Algorithmic processing of pressure waveforms to facilitate estimation of cardiac elastance
© Stevenson et al.; licensee BioMed Central Ltd. 2012
Received: 23 December 2011
Accepted: 28 May 2012
Published: 15 June 2012
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.
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.
The method located all points on 87 out of 88 waveforms for P pa , to within the sampling frequency. For P ao , out of 616 total points, 605 were found within 1%, 5 within 5%, 4 within 10% and 2 within 20%.
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.
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.
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 , P v (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 , 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 (E es ) [14, 23, 24] and ejection fraction . 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 E es  (although this cannot be found from the normalised TVE waveform). However, TVE features are highly correlated to clinical parameters  and contain similar information to pressure-volume (PV) loops, which are known to contain information on cardiac function  including cardiac work [28, 29], contractility [13, 30], O2 consumption [29, 31], and all the states of filling, contraction ejection and relaxation . 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 (P ao ) and pulmonary artery (P pa ) 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 P ao and P pa 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.
locate points on P ao and P pa (this paper)
correlate points of the pressure waveforms to points on the measured cardiac elastance
use these correlations to estimate the points on the cardiac elastance
create a continuous function e(t) through the estimated points with (2)-(4)
- 5.compare the estimated elastance waveform to the measured elastance waveform(2)
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.
Both local maximum, , and minimum, , are also required. exists and is the maximum point of , if and only if, there exists a maximum stationary point that does not fall at the temporal boundary of , and similarly for .
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, P 2, over a given range, t, near the first guess, P 1, and also defining a threshold time, D. If the point P 2 lies temporally within , then the correct point is the initial guess P 1, otherwise the correct point is P 2.
The choice of or , 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.
Point location method
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.
The error data for each location, grouped by number of points found per error band for each location
L S ∗
The automated, algorithmic method presented enables the mapping between aortic pressure (P ao ), pulmonary artery pressure (P pa ), and the ventricle TVE (e rv (t) and e lv (t)), by accurately processing the P ao and P pa 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 P ao and P pa 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.
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.
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.
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)
- Guyton A, Hall J: Textbook of Medical Physiology. Philadelphia: W.B. Saunders Company; 2000.Google Scholar
- Grenvik A, Ayres SM, Holbrook PR: Textbook of Critical Care. Philadelphia: W.B. Saunders Company; 1989.Google Scholar
- 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 care. Crit Care Med 2001, 29(7):1303–10. 10.1097/00003246-200107000-00002View ArticleGoogle Scholar
- Kearon C: Diagnosis of pulmonary embolism. CMAJ 2003, 168(2):183–94.Google Scholar
- Pineda LA, Hathwar VS, Grant BJB: Clinical Suspicion of Fatal Pulmonary Embolism. Chest 2001, 120(3):791–5. [ http://chestjournal.chestpubs.org/content/120/3/791.full] 10.1378/chest.120.3.791View ArticleGoogle Scholar
- Frazier S, Skinner GJ: Pulmonary artery catheters: state of the controversy. Journal of Cardiovascular Nursing 2008, 23: 113–21.View ArticleGoogle Scholar
- Chatterjee K: The Swan-Ganz Catheters: Past, Present, and Future. Circulation 2009, 119: 147–152. 10.1161/CIRCULATIONAHA.108.811141View ArticleGoogle Scholar
- Cooper A, Doig WJ, Sibbald GS: Pulmonary artery catheters in the critically ill. Crit Care Med 1996, 12: 777–94.Google Scholar
- Weber KT, Janicki JS: The heart as a muscle-pump system and the concept of heart-failure. Am Heart J 1979, 98(3):371–84. 10.1016/0002-8703(79)90051-6View ArticleGoogle Scholar
- Ross JJ: Afterload mismatch and preload reserve: a conceptual framework for the analysis of ventricular function. Prog Cardiovasc Dis 1976, 18: 255–64. 10.1016/0033-0620(76)90021-9View ArticleGoogle Scholar
- Suga H, Sagawa K: Instantaneous pressure-volume relationships and their ratio in the excised, supported canine left ventricle. Circ Res 1974, 35: 117–26. 10.1161/01.RES.35.1.117View ArticleGoogle Scholar
- Sagawa K: The end-systolic pressure-volume relation of the ventricle: definition, modifications and clinical use. Circulation 1981, 63(6):1223–7. 10.1161/01.CIR.63.6.1223View ArticleGoogle Scholar
- 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 ratio. Circ Res 1973, 32(3):314–22. 10.1161/01.RES.32.3.314View ArticleGoogle Scholar
- 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 application. Circulation 1996, 94(10):2497–506. 10.1161/01.CIR.94.10.2497View ArticleGoogle Scholar
- Sunagawa K, Sagawa K, Maughan WL: Ventricular interaction with the loading system. Ann Biomed Eng 1984, 12(2):163–89. 10.1007/BF02584229View ArticleGoogle Scholar
- Smith BW, Chase JG, Nokes RI, Shaw GM, Wake G: Minimal haemodynamic system model including ventricular interaction and valve dynamics. Med Eng Phys 2004, 26(2):131–9. 10.1016/j.medengphy.2003.10.001View ArticleGoogle Scholar
- 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 hemofiltration. Math Biosci 2008, 216(2):132–9. 10.1016/j.mbs.2008.08.014MathSciNetView ArticleGoogle Scholar
- Chung DC, Niranjan SC, Clark J, Bidani A, Johnston WE, Zwischenberger JB, Traber DL: A dynamic model of ventricular interaction and pericardial influence. Am J Physiol 1997, 272(6 Pt 2):H2942—62.Google Scholar
- 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 method. Annals of Intensive Care 2011., 1(33): [ http://www.annalsofintensivecare.com/content/1/1/33]Google Scholar
- 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 measurements. Comput Methods Programs Biomed 2011. in-press. [ http://www.sciencedirect.com/science/article/pii/S0169260711003014]Google Scholar
- Guarini M, Urzúa J, Cipriano A, González W: Estimation of caridac function from computer analysis of the arterial pressure waveform. IEEE Transactions on Biomedical Engineering 1998, 45(12):1420–8. 10.1109/10.730436View ArticleGoogle Scholar
- Swamy G, Kuiper J, Gudur MSR, Bari Oliver N, Mukkamala R: Continuous Left Ventricular Ejection Fraction Monitoring by Aortic Pressure Waveform Analysis. Ann Biomed Eng 2009, 37(6):1055–68. 10.1007/s10439-009-9675-4View ArticleGoogle Scholar
- 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 curve. Circulation 2000, 102(16):1983–9. 10.1161/01.CIR.102.16.1983View ArticleGoogle Scholar
- 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 failure. Acta Physiol (Oxf) 2010, 198: 37–46. 10.1111/j.1748-1716.2009.02040.xView ArticleGoogle Scholar
- Burkhoff D: Chasing the Elusive Pressure-Volume Relationship. J Am Coll Cardiol Img 2009, 2: 1282–4.View ArticleGoogle Scholar
- Starfinger C, Hann CE, Chase JG, Desaive T, Ghuysen A, Shaw GM: Model-based cardiac diagnosis of pulmonary embolism. Comput Methods Programs Biomed 2007, 87: 46–60. 10.1016/j.cmpb.2007.03.010View ArticleGoogle Scholar
- Sagawa K: The ventricular pressure-volume diagram revisited. Circ Res 1978, 43(5):677–87. 10.1161/01.RES.43.5.677View ArticleGoogle Scholar
- Suga H: Ventricular energetics. Physiol Rev 1990, 70(2):247–77.Google Scholar
- Burkhoff D, Sagawa K: Ventricular efficiency predicted by an analytical model. Regu Physiol 1986, 250(6):R1021—7.Google Scholar
- 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 area. J Cardiothorac Vasc Anesth 2006, 20(3):340–6. 10.1053/j.jvca.2006.01.014View ArticleGoogle Scholar
- Suga H: Total mechanical energy of a ventricle model and cardiac oxygen consumption. Heart 1979, 236(3):H498–5.Google Scholar
- 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 embolism. Singapore; 2005.Google Scholar
- 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 overloading. Shock 2008, 29(2):197–204.Google Scholar
- 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 efficiency. Cardiovasc Res 2003, 59(2):412–8. 10.1016/S0008-6363(03)00368-7View ArticleGoogle Scholar
- 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 pigs. Artif Organs 2006, 30(7):560–4. 10.1111/j.1525-1594.2006.00260.xView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.