Document Detail

Algorithmic processing of pressure waveforms to facilitate estimation of cardiac elastance.
Jump to Full Text
MedLine Citation:
PMID:  22703604     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
BACKGROUND: Cardiac elastances are highly invasive to measure directly, but are clinically useful due to the amount of information embedded in them. Information about the cardiac elastance, which can be used to estimate it, can be found in the downstream pressure waveforms of the aortic pressure (P(ao)) and the pulmonary artery (P(pa)). However these pressure waveforms are typically noisy and biased, and require processing in order to locate the specific information required for cardiac elastance estimations. This paper presents the method to algorithmically process the pressure waveforms.
METHODS: A shear transform is developed in order to help locate information in the pressure waveforms. This transform turns difficult to locate corners into easy to locate maximum or minimum points as well as providing error correction.
RESULTS: The method located all points on 87 out of 88 waveforms for Ppa, to within the sampling frequency. For Pao, out of 616 total points, 605 were found within 1%, 5 within 5%, 4 within 10% and 2 within 20%.
CONCLUSIONS: The presented method provides a robust, accurate and dysfunction-independent way to locate points on the aortic and pulmonary artery pressure waveforms, allowing the non-invasive estimation of the left and right cardiac elastance.
Authors:
David Stevenson; James Revie; J Geoffrey Chase; Christopher E Hann; Geoffrey M Shaw; Bernard Lambermont; Alexandre Ghuysen; Philippe Kolh; Thomas Desaive
Related Documents :
12243374 - Do ancillary properties of antihypertensive drugs explain the outcome results of recent...
7490144 - Hypertension prevalence, awareness, treatment, and control in egypt. results from the e...
6753894 - Long-term experience of captopril in the treatment of primary (essential) hypertension.
24814814 - Comparison of isoflurane and propofol for maintenance of anesthesia in dogs with intrac...
3181284 - Reproducibility of heart rate changes following adenosine infusion in man.
12144004 - Effect of a primary care physician-focused, population-based approach to blood pressure...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2012-06-15
Journal Detail:
Title:  Biomedical engineering online     Volume:  11     ISSN:  1475-925X     ISO Abbreviation:  Biomed Eng Online     Publication Date:  2012  
Date Detail:
Created Date:  2013-01-01     Completed Date:  2013-04-23     Revised Date:  2013-07-12    
Medline Journal Info:
Nlm Unique ID:  101147518     Medline TA:  Biomed Eng Online     Country:  England    
Other Details:
Languages:  eng     Pagination:  28     Citation Subset:  IM    
Affiliation:
Department of Mechanical Engineering, Centre for Bio Engineering at the University of Canterbury, Christchurch, New Zealand.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Algorithms*
Aorta / physiology
Blood Pressure*
Heart / physiology*
Humans
Models, Statistical*
Pulmonary Artery / physiology
Wavelet Analysis*
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Biomed Eng Online
Journal ID (iso-abbrev): Biomed Eng Online
ISSN: 1475-925X
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Stevenson et al.; licensee BioMed Central Ltd.
open-access:
Received Day: 23 Month: 12 Year: 2011
Accepted Day: 28 Month: 5 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 15 Month: 6 Year: 2012
Volume: 11First Page: 28 Last Page: 28
PubMed Id: 22703604
ID: 3533753
Publisher Id: 1475-925X-11-28
DOI: 10.1186/1475-925X-11-28

Algorithmic processing of pressure waveforms to facilitate estimation of cardiac elastance
David Stevenson1 Email: david.stevenson@pg.canterbury.ac.nz
James Revie1 Email: Jamejames.revie@pg.canterbury.ac.nz
J Geoffrey Chase1 Email: geoff.chase@canterbury.ac.nz
Christopher E Hann1 Email: chris.hann@canterbury.ac.nz
Geoffrey M Shaw2 Email: Geoff.Shaw@cdhb.govt.nz
Bernard Lambermont3 Email: b.lambermont@chu.ulg.ac.be
Alexandre Ghuysen3 Email: A.Ghuysen@chu.ulg.ac.be
Philippe Kolh3 Email: Philippe.Kolh@chu.ulg.ac.be
Thomas Desaive3 Email: tdesaive@ulg.ac.be
1Department of Mechanical Engineering, Centre for Bio Engineering at the University of Canterbury, Christchurch, New Zealand
2Department of Intensive Care, Christchurch Hospital, Christchurch, New Zealand
3Cardiovascular Research Center, University of Liege, Belgium

Background

In an Intensive Care Unit (ICU), cardiac disturbances are difficult to diagnose and treat, which can lead to poor management [1,2]. Inadequate diagnosis can be common, and plays a significant role in increased length of stay and death [3-5], despite access to many different cardiac measurements and metrics. Currently, internal measurements are only possible at the locations where catheters are placed. This limited set of data can severely restrict clinical diagnostic capability, and thus these catheters are not necessarily associated with improved outcomes [6-8]. Overall, a lot of data currently available to ICU clinicians, that could have significant clinical value, is under utilised.

Using modelling techniques, this limited set of data can be expanded to estimate a much greater set of clinically relevant data to enable more accurate diagnosis. For example, acute cardiovascular dysfunction, like pulmonary embolism (PE) and septic shock, severely alter \ system (CVS) hemodynamics around the heart. These changes can be seen by catheter measurements as a change in the balance of preload and afterload, resulting in an altered cardiac energetic state [9,10]. Detailed cardiac energetics are too invasive to measure in an ICU setting. However, if the relevant energetics could be captured from a nearby catheter, the clinical potential of such measurements could be realised. To date, no such method achieves this aim.

Time-varying cardiac elastance (TVE) is defined [11]:

[Formula ID: bmcM1]
(1) 
e(t)=Pv(t)Vv(t)−Vd

where, V d is assumed to be equal to V 0 for simplicity, V 0 is the intercept of the end-systolic pressure-volume relation (ESPRV) with the volume axis [12], Pv(t) is the ventricle pressure and V v(t) is the ventricle volume. It thus provides a measure of heart function and energetics [13-15]. The waveform e(t) is typically normalised to a value of 1.0 [13], and can also be used as the input or driver function in lumped parameter CVS models [16-20].

There have been several attempts to estimate TVE [14,21-24]. However, none have estimated it for its own sake. Most studies present a method using the TVE to estimate a specific parameter, most commonly end-systolic elastance (Ees) [14,23,24] and ejection fraction [22]. However, their validation is based on these metrics, not on the resulting TVE waveform.

This research is unique in that the end goal is to produce the TVE function in its own right, validating the TVE waveform on its own accuracy for eventual use as a diagnostic tool. It is unclear to date how much specific information can be obtained from the TVE waveform, other than the highly sought-after Ees[25] (although this cannot be found from the normalised TVE waveform). However, TVE features are highly correlated to clinical parameters [26] and contain similar information to pressure-volume (PV) loops, which are known to contain information on cardiac function [27] including cardiac work [28,29], contractility [13,30], O2 consumption [29,31], and all the states of filling, contraction ejection and relaxation [1]. Thus, a TVE waveform reflects cardiac state, cardiac output or blood volume, and net preload and afterload, all of which change with cardiac dysfunction. Hence, the ability to easily and non-invasively obtain TVE waveforms could enable clinically useful diagnostics and metrics.

This paper presents the first step in estimating TVE from already available measurements, namely aortic (Pao) and pulmonary artery (Ppa) pressure. This research develops algorithms to process these pressure waveforms to extract specific features and points, which, in turn, allow the estimation of the end goal, the TVE waveform. The Pao and Ppa pressure waveforms, as typically measured, are noisy and/or biased, which can significantly effect this process. Hence, a method is presented for automatically processing the pressure waveforms to robustly and accurately locate the points required by the correlations, that enable accurate, not additionally invasive cardiac elastance estimation and construction.


Methods
Concept

The specific approach presented defines a representative set of points on the TVE waveforms. These points allow a continuous waveform to be constructed that adequately follows the shape of a true (invasively measured) TVE, capturing the necessary dynamics. Correlations between these representative points and the points or properties on the typically measured waveforms, Pao and Ppa, enable construction of continuous beat-to-beat estimations of TVE with knowledge only of the Pao and Ppa waveforms and the pre-defined correlations. It thus uses typically available data to construct what would otherwise require a highly invasive added test. A high level view of this approach is presented in Figure 1, which shows the formation and use of the correlations to generate and estimated cardiac elastance waveform. This paper is focused on the left half of Figure 1, that of producing the correlations for later use. A brief overview of the whole method is given for clarity:

1. locate points on Pao and Ppa (this paper)

2. correlate points of the pressure waveforms to points on the measured cardiac elastance

3. use these correlations to estimate the points on the cardiac elastance

4. create a continuous function e(t) through the estimated points with (2)-(4)

5. compare the estimated elastance waveform to the measured elastance waveform

[Formula ID: bmcM2]
(2) 
e(t)=Fα(t)0<t<cα(1−x2)(t−cα)cβ−cα+x2cα<t<cβFβ(t)cβ<t<period

where:

[Formula ID: bmcM3]
(3) 
Fi=ai·e−bi(t−ci)2

and the coefficients of (3), also seen in (2), are fitted for a specific waveform, and are defined:

[Formula ID: bmcM4]
(4) 
aα=x2bα=−log(x1/x2)exp(log(−log(x1/x2)·2·(x1/x˙1))·2)cα=−log(x1/x2)·2·x1−x˙1·t1x˙1

where aβ, bβ and cβ are similarly defined by replacing subscript 1 with 3 and setting x2 = 1.

Figure 2 shows an illustrative mapping between points on Pao and TVE. However, this approach is useful if and only if it is possible to automate the detection of the required points, defined in (5), on the Pao and Ppa waveforms, shown in Figures 3 and 4.

This paper focuses on the robust capturing of the points on the pressure waveforms and leaves the specific correlations and methodology of creating the TVE waveforms to a paper in review. The details contained in this paper about the formulation of the correlations and there use are illustrative only, and are assumed to be correct for the purposes of demonstrating a potential use for this research.

For describing the methods in this paper, a naming convention is defined:

[Formula ID: bmcM5]
(5) 
e(t)≡time varying cardiac elastancePao≡aortic pressurePpa≡pulmonary artery pressureDMPG≡driver maximum positive gradientMN≡minimum pointMPG≡maximum positive gradientLS≡left shoulderMX≡maximumRS≡right shoulderMNG≡maximum negative gradientDN≡dicrotic notch

Shear Transform

This paper uses a shear transform to extract features from the Pao and Ppa waveforms, defined:

[Formula ID: bmcM6]
(6) 
S≡(t,X(t))→t,ϕshear(X(t))

where:

[Formula ID: bmcM7]
(7) 
ϕshearX(t)=X(t)+mt+c,t0<t<tend

[Formula ID: bmcM8]
(8) 
X(t)≡descrete, time valued data

and the parameters t0 and tend are set depending on the region or period of interest. The parameters m and c are chosen such that:

[Formula ID: bmcM9]
(9) 
ϕshear(X(t0))=X(t0)=ϕshear(X(tend))

Equation (9) leads to:

[Formula ID: bmcM10]
(10) 
X(t0)+m·t0+c=X(t0)X(tend)+m·tend+c=X(t0)

Solving (10), for m and c yields:

[Formula ID: bmcM11]
(11) 
m=X(t0)−X(tend)tend−t0c=−m·t0

To better visualize how this transformation operates, imagine a line from A to C in Figure 5, representing a portion of the waveform in Figure 3, rotated about A so that end points align horizontally, while time remains unchanged. Hence, it is a rotation and contraction that projects the line onto a horizontal axis (time). The effect of this transformation is to transform the difficult to find “shoulder” point B into an easily found peak of a curve, or for the reflection of A B in x, a minimum or the curve. A “shoulder” is defined as a point at which two near linear lines with different slopes meet, such as the point LS in Figure 3. Thus, the use of this transform makes it far easier to, algorithmically, locate aspects of the waveforms which can be otherwise difficult to find.

The transform, S, is used in two ways. First it is used to locate a point of interest based on the maximum or minimum point of S. This use is demonstrated in Figure 5. Hence, a desired point P is defined:

[Formula ID: bmcM12]
(12) 
P=SmaxorSmin

Both local maximum, Smax(local), and minimum, Smin(local), are also required. Smax(local) exists and is the maximum point of S, if and only if, there exists a maximum stationary point that does not fall at the temporal boundary of S, and similarly for Smin(local).

The second way the shear transform of (6) is used relates to the verification of a particular point given an initial guess. For example, the first guess of the point MN is the global minimum of the waveform, after which this point is verified using the shear transformation, resulting in confirmation of the point or a new point to use instead.

This works by locating the maximum of minimum point of the shear transform, P2, over a given range, t, near the first guess, P1, and also defining a threshold time, D. If the point P2 lies temporally within tP1±D, then the correct point is the initial guess P1, otherwise the correct point is P2.

The choice of Smax or Smin, the range of time, t, and the threshold time, D, are defined for the type of point under consideration. The specific values, listed in Section Point location method, are chosen empirically, based on what features that appear close to the point of interest and the temporal variation that has been observed in these features.

These two situations are graphically shown in Figure 6 and Figure 7 for positive values of D. However the same applies for negative values, for which the real point lies before the point P1, instead of after it. The complete process is defined:

[Formula ID: bmcM13]
(13) 
P=P1iftP2lies temporally withintP1±DP2otherwise

where:

[Formula ID: bmcM14]
(14) 
P1≡initial point to be checkedP2≡SmaxorSmin(chosen separately)D≡threshold timet∈{t:t0<t<tend}

Combined these two methods of use, shown in Figure 5-7, create a robust and computationally fast method for locating certain hard to find points on a waveform.

Point location method

The method for finding the points is described in Figure 8 along with the following two sections (Finding DMPG and Finding DN). Figure 8 gives the full text, reproducible and ordered method except for the points DMPG and DN (which are described in the next two sections) along with a graphical illustration. The graphical illustrations are the out-working of the method for a representative waveform, and are only intended to aid the reader in their understanding of the method, and not to formally describe the method itself. Due to the complexity of the method for the points DMPG and DN, these two have been described in Figure 8 only for the simplest (as well as the and most common) case, with the full method described in separate sections with relevant figures.

Finding DMPG

There are some cases, where DMPG, see Figure 3, is equivalent to MN. These cases occur when:

[Formula ID: bmcM15]
(15) 
(tMX−tMN)<period·0.25

When (15) is not true, DMPG is found as the Smax(local) from 34·tMN to MN, see Figure 9. However, there are a few cases, for both sepsis and pulmonary embolism, in which a local maximum of S does not exist except at the boundaries of the region, which is not acceptable if an automated detection method is desired. In this case, a point P2, is defined as the Smin from 34·tMN to MN. If Smax(local) from P2 to MN exists, this is taken as DMPG, see Figure 10, otherwise DMPG is defined as Smax(local) from 34·tMN to P2, see Figure 11. If this final local maximum does not exist, DMPG is defined as 34·tMN. These cases occur due to noise, variability and dysfunction and are part of what makes robust algorithmic or automated processing difficult.

This process is defined:

[Formula ID: bmcM16]
(16) 
P=P1if∃P1P3if∃P3P4if∃P434·tMNotherwise

where:

[Formula ID: bmcM17]
(17) 
P1≡Smax(local),t∈{t:34·tMN<t<tMN}P2≡Smin(local),t∈{t:34·tMN<t<tMN}P3≡Smax(local),t∈{t:tP2<t<tMN}P4≡Smax(local),t∈{t:34·tMN<t<tP2}

Finding DN

The general approach to find the point DN, see Figure 3, is to find Smin between tMX and tend (or period). However, in a number of cases this approach fails due to oscillations towards the end of the waveform, see Figure 12. Also, using only the first local minimum (as is the case in Figure 12) works only in a few cases and therefore is not a robust solution either. Hence a more specific algorithm is required.

A second maximum point is defined:

[Formula ID: bmcM18]
(18) 
MX2≡max point ofPao(t),t∈t:tMN+period5<t<period

and two more intermediate points are defined:

[Formula ID: bmcM19]
(19) 
DN1≡lowestSmin(local),t∈{t:tMX2<t<period}

[Formula ID: bmcM20]
(20) 
DN2≡firstSmin(local),t∈{t:tMX2<t<period}

From DN1 and DN2, the real DN is chosen, defined:

[Formula ID: bmcM21]
(21) 
DN=DN1if(CA∧CB)∨(CC∧CD)DN2otherwise

where:

[Formula ID: bmcM22]
(22) 
CA≡DN2>MX2CB≡(tDN2−tMX2)<0.15·periodCC≡DN2−DN1Pao(tMX)>0.02CD≡(tDN2−tMX2)<0.11·period

While using a measured Ppa waveform, DN is defined:

[Formula ID: bmcM23]
(23) 
DN=DN1

Validation Test

The method presented was developed on a set of five pigs (51 waveforms) that were induced with pulmonary embolism [32,33], and then independently tested on a further five pigs (37 waveforms) induced with septic shock, and treated with haemofiltration [34,35].

The points for all waveforms (see Figures 3-4) were identified or checked individually by eye. The two gradients (MPG and MNG) were first located through simple computation, and then individually checked by eye and corrected were necessary. Due to the nature and location of these points (maximum gradient of a sigmoidal function), they are the two easiest and most reliable to find algorithmically, and in fact the algorithmic approach is more accurate than hand selection. The two shoulders (LS and RS) were first located through the algorithm developed prior to that which is described in this paper, after which each point was individually checked and corrected. Because there is no formal definition for the location of these “shoulder” points, it was left to an algorithmic definition. For a validation test this definition is self fulfilling. However, as LS (RS is only used to aid in finding MNG , and is hence not included in the validation results) is found as an intermediate step to the estimation of the cardiac elastance, its full and more formal validation would be the results of the cardiac estimation which is not in the scope of this paper. All the remaining points were hand selected.

The automated method was applied to the waveforms and the identified points assessed against the known points for accuracy in time. The use of separate data with different cardiac dysfunction to design and test the method ensures the robustness of the validated method.


Results

For the points (MN and DN) required when using Ppa (Figure 4), the method located both points in 87 of the 88 waveforms to within the sample frequency of 200Hz (0.005 sec), missing DN, from one waveform. This missed point is in a waveform at the start of the third pig of the sepsis cohort and is unique to the data set, both in the measured TVE and Ppa, as shown in Figure 13, compared to the more typical Ppa waveform in Figure 4. The failure is due to the unusual second peak of Ppa, and the early decay of the TVE.

For each Pao waveform, the method locates eight points, MX, MN, DMPG, LS, MPG, DN, RS, MNG. However, RS and MN are only used to aid the location of other points. These two points were both located sufficiently to enable the method to progress in all 88 waveforms. Results for the other six points are shown in Table 1. Of 616 total points, 605 were found within 1%, 5 within 5%, 4 within 10% and 2 within 20%.


Discussion

The automated, algorithmic method presented enables the mapping between aortic pressure (Pao), pulmonary artery pressure (Ppa), and the ventricle TVE (erv(t) and elv(t)), by accurately processing the Pao and Ppa waveforms to identify specific points. Once combined, they enable a very useful tool, for clinicians to obtain very accurate TVE without further invasive or risky sensors or procedures.

There are other ways to locate the points on the pressure waveforms, most notably a derivative and second derivative method. However, this becomes problematic in practice due the noise inherent in the waveforms. The method that has been developed in this paper, was designed to work with the level of noise that is typically seen on these measurements, and is therefore more involved than a simple derivative method.

The method presented was robust to the typical and significant variation and noise in Pao and Ppa waveforms. The method was developed on five pigs induced with pulmonary embolism, and then tested independently with data from another set of five pigs induced with septic shock. The results give confidence that this method will generalize to a wider set of disease states and to human data.

However, while the results were very good, this research needs further validation on a wider cohort of pigs and types of dysfunction to further quantify the limits and accuracy of this approach. Direct validation on humans is the ultimate goal. However, the results appear robust, and justify and enable a wide range of further more in-depth validation studies of both the method and its potential uses when reconstructing TVE for monitoring and diagnosis.

Clinically, it must be noted that for this method to work a Swan-Ganz catheter is assumed. If radial artery pressure was measured instead, there would be more oscillations in the waveform, potentially requiring modifications. However, Swan-Ganz catheters are still commonly used, and this application would add value to their use, which is otherwise sometimes contested [6-8].

The method developed in this paper shows promise for gaining clinical insight and improving diagnosis. It can enable clinicians to get more information about the current patient state, without the use of more invasive measurements, as well as beat-to-beat tracking of this information. This level of detail is far more than currently available and could potentially lead to better and earlier diagnosis of dysfunction, as well as better knowledge of response to treatment, non-invasively, as it needs no further procedures or sensors required.


Conclusions

This paper has presented a robust, potentially dysfunction-independent method to find the waveform points necessary to use proven methods to non-invasively and automatically estimate the otherwise unavailable left and right ventricle TVEs with accuracy well within measurement error. This capability is enabled using standard measurements that are already commonly used in an intensive care setting, thus involving no additional risk to the patient. The results thus justify prospective validation of these conclusions.


Competing interests

The authors declare that they have no competing interests.


Author’s contributions

DS drafted the manuscript and developed the algorithm. JR and JC participated in the algorithm development with added input from TD. CH participated in the initial mathematical formulation. GS provided physiological understanding and clinical input at all stages. BL, AG, PK and TD provided the porcine data and further clinical input and relevance. JC and TD edited and aided the writing of the manuscript and revisions with DS. All authors read and approved the final manuscript.


Acknowledgements

This work was supported by the Foundation for Research Science and Technology (FRST) of New Zealand, the Fonds de la Recherche Scientifique (FNRS, Belgium) and the French Community of Belgium (Actions de Recherches Concertées Académie Wallonie-Europe)


References
Guyton A,Hall J,Textbook of Medical PhysiologyYear: 2000Philadelphia: W.B. Saunders Company
Grenvik A,Ayres SM,Holbrook PR,Textbook of Critical CareYear: 1989Philadelphia: W.B. Saunders Company
Angus DC,Linde-Zwirble WT,Lidicker J,Clermont G,Carcillo J,Pinsky MR,Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of careCrit Care MedYear: 200129713031010.1097/00003246-200107000-0000211445675
Kearon C,Diagnosis of pulmonary embolismCMAJYear: 200316821839412538548
Pineda LA,Hathwar VS,Grant BJB,Clinical Suspicion of Fatal Pulmonary EmbolismChestYear: 200112037915 [ http://chestjournal.chestpubs.org/content/120/3/791.full]. 10.1378/chest.120.3.79111555511
Frazier S,Skinner GJ,Pulmonary artery catheters: state of the controversyJournal of Cardiovascular NursingYear: 2008231132118382253
Chatterjee K,The Swan-Ganz Catheters: Past, Present, and FutureCirculationYear: 200911914715210.1161/CIRCULATIONAHA.108.81114119124674
Cooper A,Doig WJ,Sibbald GS,Pulmonary artery catheters in the critically illCrit Care MedYear: 19961277794
Weber KT,Janicki JS,The heart as a muscle-pump system and the concept of heart-failureAm Heart JYear: 19799833718410.1016/0002-8703(79)90051-6382811
Ross JJ,Afterload mismatch and preload reserve: a conceptual framework for the analysis of ventricular functionProg Cardiovasc DisYear: 1976182556410.1016/0033-0620(76)90021-9128034
Suga H,Sagawa K,Instantaneous pressure-volume relationships and their ratio in the excised, supported canine left ventricleCirc ResYear: 1974351172610.1161/01.RES.35.1.1174841253
Sagawa K,The end-systolic pressure-volume relation of the ventricle: definition, modifications and clinical useCirculationYear: 19816361223710.1161/01.CIR.63.6.12237014027
Suga H,Sagawa K,Shoukas AA,Load independence of the instantaneous pressure-volume ratio of the canine left ventricle and effects of epinephrine and heart rate on the ratioCirc ResYear: 19733233142210.1161/01.RES.32.3.3144691336
Senzaki H,Chen CH,Kass DA,Single-beat estimation of end-systolic pressure-volume relation in humans. A new method with the potential for noninvasive applicationCirculationYear: 19969410249750610.1161/01.CIR.94.10.24978921794
Sunagawa K,Sagawa K,Maughan WL,Ventricular interaction with the loading systemAnn Biomed EngYear: 19841221638910.1007/BF025842296507965
Smith BW,Chase JG,Nokes RI,Shaw GM,Wake G,Minimal haemodynamic system model including ventricular interaction and valve dynamicsMed Eng PhysYear: 2004262131910.1016/j.medengphy.2003.10.00115036180
Starfinger C,Chase JG,Hann CE,Shaw GM,Lambermont B,Ghuysen A,Kolh P,Dauby PC,Desaive T,Model-based identification and diagnosis of a porcine model of induced endotoxic shock with hemofiltrationMath BiosciYear: 20082162132910.1016/j.mbs.2008.08.01418817788
Chung DC,Niranjan SC,Clark J,Bidani A,Johnston WE,Zwischenberger JB,Traber DL,A dynamic model of ventricular interaction and pericardial influenceAm J PhysiolYear: 19972726 Pt 2H2942—629227574
Revie J,Stevenson D,Chase G,Hann CE,Lambermont B,Ghuysen A,Kolh P,Morimont P,Shaw GM,Desaive T,Clinical detection and monitoring of acute pulmonary embolism: proof of concept of a computer-based methodAnnals of Intensive CareYear: 2011133 [ http://www.annalsofintensivecare.com/content/1/1/33].
Revie J,Stevenson D,Chase JG,Hann CE,Lambermont B,Ghuysen A,Kolh P,Shaw GM,Heldmann S,Desaive T,Validation of subject-specific cardiovascular system models from porcine measurementsComput Methods Programs BiomedYear: 2011 in-press. [ http://www.sciencedirect.com/science/article/pii/S0169260711003014].
Guarini M,Urzúa J,Cipriano A,González W,Estimation of caridac function from computer analysis of the arterial pressure waveformIEEE Transactions on Biomedical EngineeringYear: 199845121420810.1109/10.7304369835191
Swamy G,Kuiper J,Gudur MSR,Bari Oliver N,Mukkamala R,Continuous Left Ventricular Ejection Fraction Monitoring by Aortic Pressure Waveform AnalysisAnn Biomed EngYear: 200937610556810.1007/s10439-009-9675-419308732
Shishido T,Hayashi K,Shigemi K,Sato T,Sugimachi M,Sunagawa K,Single-beat estimation of end-systolic elastance using bilinearly approximated time-varying elastance curveCirculationYear: 2000102161983910.1161/01.CIR.102.16.198311034949
ten Brinke EA,Klautz RJ,Verwey HF,van der Wall EE,Dion RA,Steendijk P,Single-beat estimation of the left ventricular end-systolic pressure-volume relationship in patients with heart failureActa Physiol (Oxf)Year: 2010198374610.1111/j.1748-1716.2009.02040.x19735484
Burkhoff D,Chasing the Elusive Pressure-Volume RelationshipJ Am Coll Cardiol ImgYear: 2009212824
Starfinger C,Hann CE,Chase JG,Desaive T,Ghuysen A,Shaw GM,Model-based cardiac diagnosis of pulmonary embolismComput Methods Programs BiomedYear: 200787466010.1016/j.cmpb.2007.03.01017478006
Sagawa K,The ventricular pressure-volume diagram revisitedCirc ResYear: 19784356778710.1161/01.RES.43.5.677361275
Suga H,Ventricular energeticsPhysiol RevYear: 1990702247772181496
Burkhoff D,Sagawa K,Ventricular efficiency predicted by an analytical modelRegu PhysiolYear: 19862506R1021—7
Broscheit JA,Weidemann F,Strotmann J,Steendijk P,Karle H,Roewer N,Greim CA,Time-varying elastance concept applied to the relation of carotid arterial flow velocity and ventricular areaJ Cardiothorac Vasc AnesthYear: 2006203340610.1053/j.jvca.2006.01.01416750733
Suga H,Total mechanical energy of a ventricle model and cardiac oxygen consumptionHeartYear: 19792363H4985
Desaive T,Dutron S,Lambermont B,Kolh P,Hann CE,Chase JG,Dauby PC,Ghuysen A,Close-loop model of the cardiovascular system including ventricular interaction and valve dynamics: application to pulmonary embolismYear: 2005Singapore
Ghuysen A,Lambermont B,Kolh P,Tchana-Sato V,Magis D,Gerard P,Mommens V,Janssen N,Desaive T,D’Orio V,Alteration of right ventricular-pulmonary vascular coupling in a porcine model of progressive pressure overloadingShockYear: 200829219720417693928
Lambermont B,Ghuysen A,Kolh P,Tchana-Sato V,Segers P,Gerard P,Morimont P,Magis D,Dogne JM,Masereel B,D’Orio V,Effects of endotoxic shock on right ventricular systolic function and mechanical efficiencyCardiovasc ResYear: 2003592412810.1016/S0008-6363(03)00368-712909324
Lambermont B,Delanaye P,Dogne JM,Ghuysen A,Janssen N,Dubois B,Desaive T,Kolh P,D’Orio V,Krzesinski JM,Large-pore membrane hemofiltration increases cytokine clearance and improves right ventricular-vascular coupling during endotoxic shock in pigsArtif OrgansYear: 2006307560410.1111/j.1525-1594.2006.00260.x16836739

Figures

[Figure ID: F1]
Figure 1 

Overview. The figure shows a conceptualised overview of the process described in this paper and further implications. From the many measured left cardiac elastance (e(t)) waveforms, alone with many aortic pressure (Pao) waveforms, correlations are derived (the information flow is shown through the large grey arrow). Once these correlations are known, they can be used along with the aortic pressure waveform (from a patient), to arrive at an estimation of their cardiac elastance waveform. The equivalent for the right cardiac elastance is also shown, with the pulmonary artery pressure (Ppa).



[Figure ID: F2]
Figure 2 

Illustrative elastance estimation. An example of what can be done with the identified points on the aortic pressure, and an example of the formation of the estimated cardiac elastance is shown here, while the terms are defined in (5). This figure is not part of the method of this paper, rather as a illustration of what the method as a whole leads to.



[Figure ID: F3]
Figure 3 

Aortic pressure waveform and relevant points. A representative aortic pressure waveform over one heart beat with relevant points (defined in (5)) marked on it. The two dashed circles, MN and RS are used only in locating other points.



[Figure ID: F4]
Figure 4 

Pulmonary artery pressure waveform and relevant points. A representative pulmonary artery pressure waveform over one heart beat with relevant points (defined in (5)) marked on it. The dashed circle, MX is only used to help find other points.



[Figure ID: F5]
Figure 5 

Shear transform. An illustration of the shear transformation of (6), turning a hard to locate “shoulder” ( ) points into an easily found maximum point ( ).



[Figure ID: F6]
Figure 6 

Use of the shear transform, A. The desired point for MN is P. However, in this example the global minimum of the waveform is P1, which is the initial guess for MN. A shear transform of the pressure waveform between P1 and MX reveals a minimum (P2) outside the range of D, and hence the time of P2 is taken as the time of MN.



[Figure ID: F7]
Figure 7 

Use of the shear transform, B. This example is the other situation in the process of finding MN to Figure 6, i.e. the initial guess of the global minimum for MN is correct. Here, the minimum of the shear transform from P1 to MX falls within the range of D and hence the P1 is taken as MN.



[Figure ID: F8]
Figure 8 

The method. The step by step method for finding the points on Pao and Ppa, as labelled on the right. The graphics beside each step are for illustration only and are not meant to be part of the definition of the method, rather to see the method in operation on a representative Pao waveform. Note that the methods described here for DMPG and DN are note complete as these require a more complex method, refer to Sections Finding DMPG and Finding DN for the complete method for these two points.



[Figure ID: F9]
Figure 9 

Finding DMPG, A. A straight forward case for finding DMPG, where P1 of (17) exists, hence DMPG P1.



[Figure ID: F10]
Figure 10 

Finding DMPG, B. A less common case for finding DMPG, where P1 of (17) does not exist, but P3 does, hence DMPG P3.



[Figure ID: F11]
Figure 11 

Finding DMPG, C. A less common case for finding DMPG, where P1, and P3 of (17) do not exist, but P4 does, hence DMPG P4.



[Figure ID: F12]
Figure 12 

Finding DN. An example of where the first local minimum of the shear transform is the correct time for the point DN.



[Figure ID: F13]
Figure 13 

Where the method fails.Ppa alongside the matching TVE. The automatic or algorithmic method failed to capture the correct DN point (circle), the real DN and associated MX are marked by squares.



Tables
[TableWrap ID: T1] Table 1 

The error data for each location, grouped by number of points found per error band for each location


  <1% 15% 510% 1020%
DMPG
84
1
3
0
MPG
88
0
0
0
LS
86
1
1
0
MX
88
0
0
0
MNG
87
0
0
1
DN
86
2
0
0
TOTAL 605 5 4 2

Accuracy of the method: number of points grouped by absolute error (of 88 total points). * note that the validation LS here is partially self fulfilling and only included here for completeness.



Article Categories:
  • Research


Previous Document:  Incidence and prevalence of idiopathic inflammatory myopathies among commercially insured, Medicare ...
Next Document:  Subclinical impairment of ovarian reserve in juvenile systemic lupus erythematosus after cyclophosph...