# Investigations concerning the application of the cross-correlation method in cardiac output measurements

- Maciej Gawlikowski
^{1, 2}Email author and - Tadeusz Pustelny
^{2}

**11**:24

https://doi.org/10.1186/1475-925X-11-24

© Gawlikowski and Pustelny; licensee BioMed Central Ltd. 2012

**Received: **6 October 2011

**Accepted: **2 April 2012

**Published: **20 May 2012

## Abstract

### Background

In spite of numerous non-invasive examinations the “gold clinical standard” of cardiac output measurements is the invasive pulmonary artery catheterization by means of the Swan-Ganz catheter and the application of the thermodilution method to estimate the blood flow. The results obtained by means of thermodilution are sensitive to many physical and biological disturbances. The unreliability of this method amounts to 20-45% and depends on the given variant of the method. Therefore some other method, more accurate and resistant to disturbances, was looked for. This paper presents a new approach to cardiac output measurements, based on cross-correlation signal analysis. The goal of investigations was to verify experimentally the application of the cross-correlation method of cardiac output measurements.

### Results

In 99.2% of the examined cases the extreme of the cross-correlation function was easy to be estimated by numerical algorithms. In 0,8% of the remaining cases (with a plateau region adjacent to the maximum point) numerical detection of the extreme was inaccurate. The typical unreliability of the investigated method amounted o 5.1% (9.8% in the worst case). Investigations performed on a physical model revealed that the unreliability of cardiac output measurements by means of the cross-correlation method is 3–5 times better than in the case of thermodilution.

### Conclusions

The performed investigations and theoretical analysis have shown, that the cross-correlation method may be applied in cardiac output measurements. This kind of measurements seems to be more accurate and disturbance-resistant than clinically applied thermodilution.

## Keywords

## Background

The most fundamental hemodynamic parameter is the cardiac output (*CO*). It is defined as the average blood flow which is pumped through the cardiac muscle in the course of one minute [1, 2]. For over 40 years pulmonary artery catheterization (*PAC*) and *CO* measurements by means of thermodilution constituted the “gold standard” of hemodynamics diagnostics [2–6]. Nowadays, invasive examinations are slowly taken off the intensive care wards and the no-invasive one (e.g. echocardiography) or minimally invasive (e.g. PiCCO [7], LidCO [8]) have been put into practice. It should be emphasized that non-invasive methods of *CO* measurements are recognized as less accurate [2, 3, 6, 8] and minimally invasive methods must be calibrated before examination by means of thermodilution [7, 8]. Actually there are a few methods based on other physical principles [9, 10], but nowadays they are being developed. The *PAC* technique by means of the Swan-Ganz catheter allows to obtain or calculate other valuable hemodynamic parameters basing on intracardiac pressures, e.g. pulmonary artery wedge pressure, pulmonary resistance etc. [2, 6]. It has been evidently proved that invasive *PAC* examinations did not increase the mortality rate [11]. Therefore, in spite of their invasiveness, *PAC* and *CO* measurements by thermodilution still constitutes the “gold standard” and they are frequently utilized in intensive care wards.

*PAC*method requires the insertion of a special catheter (Swan-Ganz catheter) into the pulmonary artery [1–3, 6]. A miniaturized thermistor is located at the tip of the catheter. The examination runs as follows: first the indicator (mostly iced or room-temperature 0.9%NaCl) is injected into the right atrium by a duct located in the interior of the catheter. Next, the temperature variations triggered by the indicator are measured by the thermistor and the indicator dilution curve (

*IDC*) is registered by an external computer (most often by a patient monitor). The

*CO*value is calculated by means of the Stewart-Hamilton formula

**(1)**. The estimation of the unreliability of the thermodilution method is rather complicated because of the lack of other reference methods or means of significant influences on the examined object (human or experimental animal) [8, 12–14]. Therefore the accuracy is determined theoretically. Comprehensive investigations of the unreliability of the thermodilution method (theoretical and performed on a physical model) revealed a high susceptibility of this method to disturbances and a low repeatability of measurements. The typical unreliability of “iced” and “room-temperature” thermodilution proved to amount to 18% and 27%, respectively (the worst cases were 34% and 54%, respectively) [13, 15].

where: *c*
_{
b
}
*, c*
_{
i
} – specific heat of the blood and indicator, *ρ*
_{
b
}
*, ρ*
_{
i
} – mass density of the blood and indicator, *V*
_{
i
} – indicator volume, *T*
_{
b
}
*, T*
_{
i
} – temperature of the blood and indicator.

### Goal

In spite of certain imperfections, thermodilution is essentially a clinical method of *CO* measurements. Low accuracy and repeatability cause that prospecting other measurement methods, based on other physical principles, is justified.

One of the flow measurement methods (occasionally applied for atypical liquid or intricate objects) is the cross-correlation [13, 16, 17]. The goal of investigations was to verify experimentally the application of the cross-correlation method for cardiac output measurements.

## Methods

### Background of cross-correlation flow measurements

*v*

_{ AV }(see Figure 1). A disturbance (naturally or artificially inducted) of the flow morphology appears at point

*A*and moves to point

*B*without any change of its form. The disturbance is detected at the points

*A*and

*B*. The distance between the sensors amounts to

*L*. Supposing, that the stochastic process is ergodic and stationary and that the disturbance does not change its form between the points

*A*and

*B*, the signal registered in point B may be defined by

**(**2

**)**:

*τ*– transit time of disturbance.

*R*

_{ AB }calculated at

*A(t)*and

*B(t)*signals proves to achieve its extreme at

*t*=

*τ*[17]:

where: *T* – finite time of integration.

**(**4

**)**:

The relationship between the transit time of disturbance and its velocity is more complex and depends on the spatial distribution of the flow. At a laminar flow through the cylindrical vessel the spatial distribution of the flow velocity is parabolic [18]. For this kind of flow the average velocity *v*
_{
AV
} equals half the maximum velocity *v*
_{
max
}.

**(**5

**)**[19]. The value of the parameter

*m*depends on the character of the flow defined by Reynold’s number and is presented in Figure 2.

*x*– distance from the centre line of the vessel,

*R*– radius of the vessel,

*v*

_{ max }– flow velocity at the centre line of the vessel.

**(**6

**)**.

where: *Q* – volumetric flow, *v*
_{
AV
} – average flow velocity, *A* – area of the cross section of the vessel.

### The procedure of the experiment

*IDC*obtained (at the same flow rates and examination conditions) of the patient and the model (Figure 4).

Experiments were performed on the pulsating flow at the following conditions: heart rate 40–70 beats per minute, duty cycle 25-50%.

*(0)*, right ventricle

*(1)*, pulmonary artery

*(2)*and at the pulmonary artery bifurcation

*(3)*(Figure 5). The distance between the detectors

*(0)*and

*(3)*was 35 cm. The disturbance was excited in the right atrium. As the operation fluid 30% water-glycerine solution was applied (mass density

*ρ*= 1067 kg/m

^{3}, absolute viscosity

*μ*= 0.004 Pa·s). During the experiments the temperature of the operation liquid was stabilized in the range of 36.5 – 37.1°C.

After the injection of the indicator the signals originated by the sensors *(0)*…*(3)* were registered (sampling frequency 1 kHz). Next, the registered signals were processed (LabVIEW, SignalProcessing Toolkit) in the following way: filtering (digital Bessel low-pass filter, cutoff frequency 0.25 Hz), normalization, cross-correlation calculation (Fourier’s method), searching for the cross-correlation function extreme, estimation of the transit time *τ* of the disturbance. The extreme of the cross-correlation function was determined by the differential method. The last step was the calculation of the volumetric flow value by means of the algorithm described below:

● step 1: calculation of *v*
_{
max
} value from **(4)** – signals from the detectors *(0)* and *(3)* were analyzed,

*v*

_{ PTmax }was calculated by means of

**(**4

**)**. The transit time

*τ*

_{ 23 }was calculated as the difference between

*τ*

_{ 03 }and

*τ*

_{ 02 }. Those signals are received by detectors located in a vessel with a constant diameter, as opposed to the sensors

*(0)*and

*(1)*, which are located in the atrium and ventricle. For calculations of the value

*Re*the following medium parameters were assumed [13, 20]: mass density

*ρ*= 1067 kg/m

^{3}, absolute viscosity

*μ*= 0.004 Pa·s. The

*v*

_{ PTmax }is a maximum value of the continuous flow velocity. Because of the pulsatory character of the flow, the actual value of

*v*

_{ PTmax }is greater than for the continuous one and it depends on the duty cycle. Therefore

*v*

_{ PTmax }was corrected by the factor

*k*

_{ p }and the respective equations are presented below:

*k*

_{ p }– correction factor,

*L*

_{ 23 }– distance between the detectors (2) and (3),

*τ*

_{ 23 }–transit time of the disturbance between the detectors (2) and (3).

where: *T*
_{
SYS
} –ejection time of the liquid during its pulsatory flow, *T* – period of flow pulsation.

● step 3: from *Re* the value of the factor *m* in **(** 5**)** is estimated, basing on Figure 2,

● step 4: the average flow velocity *v*
_{
AV
} is calculated making use of the equation **(** 5**)**,

● step 5: the volumetric flow is calculated by means of **(** 6**)**.

## Results

In the applied physical model of pulmonary circulation the pulmonary trunk diameter was *D*
_{
PT
} = 1.3 cm. As the liquid flows through elements with different diameters (right ventricle and pulmonary trunk), the equivalent diameter of the vessel was assumed to be *A*
_{
z
} = 4 cm^{2}
[13].

The total amount of the measured flow set points was *N* = 31. For all set points *n* = 8 the transit time was measured. The measurements were performed for a flow in the range of 1.3 – 4.0 L/min (the corresponding average and maximum flow velocity was from 0.10 to 0.26 m/s and 0.20 – 0.55 m/s, respectively). Those values correspond to the maximum flow velocity through the pulmonary artery valve (0.6 m/s, refer to [1]). At these conditions the Reynolds number was in the range of 700 – 2200.

*(0)*and

*(3)*in the consecutive stages of signal processing are presented in Figures 6, 7, 8.

In 99.2% of the examined cases the extreme of the cross-correlation function was easy to be estimated by numerical algorithms. In 0,8% of the remaining cases (with a plateau region adjacent to the maximum point) numerical detection of the extreme was inaccurate.

*vs.*measured disturbance transit time) of the investigated cross-correlation method is being compared.

**Unreliability of the cross-correlation method**

Unreliability | Value | Remarks |
---|---|---|

Average | 5.1% | Standard deviation = ±3.5% |

The worst case | 9.8% | For fuzzy maximum of the cross-correlation function |

## Discussion

One of the assumptions of the cross-correlation method is the immutability of the disturbance signal form received by the detectors [16, 17]. In the performed investigations the disturbance was excited artificially by injecting the indicator. In this case the shape of the disturbance signal is modified mainly by a diffusive dilution of the indicator. This phenomenon was discussed theoretically [13] and defined in the form of the stochastic *Local Diffusion Random Walk* model of dilution [13, 15]. Actually, the effect of the modification of the shape of the disturbance signal can be observed on raw signals received by the detectors (see: Figure 6). In spite of that a sharp extreme is to be observed in the cross-correlation function (Figure 8). Only in 0.8% of the analyzed cases the maximum of the cross-correlation function was fuzzy, resulting in a poor accuracy of estimation of the transit time. It was proved [17] that the recognition of the function extreme by means of the weighted mean method allows to achieve a 30% lower systematic error than by simple max. or min. searching (e.g. by differentiation of the analysed function).

Transit times lower than 420 ms were immeasurable by means of the presented method. The reason of this effect has not been explained unambiguously (low sampling frequency and near-field of thermal disturbance were taken into account). Therefore, in order to estimate of *Re*, *τ*
_{
23
} was calculated as a *τ*
_{
03
} – *τ*
_{
02
}. The transit time *τ*
_{
03
} was well-measurable and correlated with the reference flow. The *τ = f(Q)* dependence is non-linear but it may be precisely approximated by the power function (*R*^{
2
} = 0.961). The reason of nonlinearity has not been unambiguously explained and the best fitting by power function was found experimentally.

**Unreliability of various cardiac output measurement methods comparison**

Method | Unreliability | |
---|---|---|

Typical | The worst case | |

Iced thermodilution | 18.3% | 26.9% |

Room-temperature thermodilution | 34.2% | 54.7% |

Cross-correlation method | 5.1% | 9.8% |

*CO*measurements in comparison with thermodilution are: a better accuracy and higher resistance to disturbances. The method may be carried out by a simple modification of the typical construction of the Swan-Ganz catheter (one extra thermistor must be added), which allows to retain the possibility of intracardiac pressure measurements. A disadvantage of

*CO*measurements by means of cross-correlation is the need of estimating the pulmonary trunk (

*PT*) diameter. This may be performed by transthoracic echocardiography (one of examinations performed by the authors has been presented in Figure 10) or by assuming the estimated diameter of this vessel basing on other factors (e.g. age, gender,

*BMI*, diameter of other vessels – additional investigations are required).

There are many other methods of *CO* measurements and some of them have been applied clinically. Among non-invasive methods the most significant one is echocardiography. In Doppler’s echocardiography the *CO* is calculated by means of the blood velocity [13, 22]. The spatial profile of the flow of blood through the specific vessel is unknown (in calculations a laminar flow is assumed), therefore the results of measurements are inaccurate. Another approach is based on the estimation of the left ventricle volume by means of its two-dimensional cross section [13, 23]. The estimation is realized by rough equations which influents the uncertainty of measurements. Another non-invasive method is rheocardiography [3, 5]. The method itself was designed for continuously monitoring of tthe hemodynamic parameters in healthy patients (aircraft pilots). It was clinically proved [5] that with reference to patients with cardiac a illness this method is suitable for monitoring long-term tendency rather than isolated measurement. Most of the modern minimally invasive methods based on the pulse contour analysis (e.g. PiCCO [7] and LidCO [8]) assume an immutable transfer function between excitation (blood flow generated by the heart) and response (pressure in peripheral artery). Generally this transfer function is unknown; therefore, its empirical estimation by means of other methods (most often by thermodiluion) is indispensable.

In fact, the accuracy of the methods mentioned above is unknown. In many papers only the correlation of the investigated method with the reference is estimated. Mostly, the reference constitutes thermodilution (or - hardly ever - the Fick method) with a typical accuracy in the range of 18-35% [2, 12, 13].

The application of the cross-correlation method is not confined to cardiac output measurements only. One of the most fundamental problems connected with the mechanical heart supporting therapy is the estimation of the flow generated by the heart prosthesis [9, 24–27]. The cross-correlation method with a spontaneously excited disturbance (e.g. the blood pressure variation produced by artificial valves) may be applied in ventricular assisting devices.

## Conclusions

The performed investigations and theoretical analysis have shown, that the cross-correlation method may be applied in cardiac output measurements. This kind of measurements seems to be more accurate and disturbance-resistant than clinically applied thermodilution. The idea of cross-correlation cardiac output measurements has been patented [15] (patent pending).

The presented investigations were performed on one geometry of physical models of heart and vessels. Therefore the obtained results should be treated as preliminary ones. Investigations should be continued on other geometries and – in future – on animals.

## Declarations

### Acknowledgements

This work was supported by the Polish Ministry of Science and Higher Education (grant No. N518 336135).

## Authors’ Affiliations

## References

- Guyton AC:
*Textbook of medical physiology*. Philadelphia: W.B. Saunders Company; 1991.Google Scholar - Headley JM:
*Invasive hemodynamic monitoring: physiological principles and clinical applications*. Irvine: Edwards Lifescience; 2002.Google Scholar - Guzik P, Greberski K, Wysocki H:
**The comparison of invasive and non-invasive methods of the hemodynamic parameters measurements.***Nowiny Lekarskie*2002,**6:**349–354. (in Polish).Google Scholar - Gawlikowski M, Pustelny T, Przywara-Chowaniec B, Struk P:
**The anatomic structure of pulmonary arteries as a source of unreliability in thermodilution cardiac output measurement.***Acta Phys Pol A*2008,**114:**81–89.View ArticleGoogle Scholar - Przywara-Chowaniec B, Polonski L, Gawlikowski M, Pustelny T:
**Clinical studies of rheocardiography application to hemodynamic monitoring of patients with dilated cardiomyopathy.***Acta Phys Pol A*2009,**116:**380–382.View ArticleGoogle Scholar - Lichtental PR:
*Cardiopulmonary care*. Irvine: Edwards Lifescience; 2002.Google Scholar - Goedje O, Höke K, Goetz AE,
*et al*.:**Reliability of a new algorithm for continuous cardiac output determination by pulse-contour analysis during hemodynamic instability.***Crit Care Med*2002,**30**(1):52–8. 10.1097/00003246-200201000-00008View ArticleGoogle Scholar - Linton RA, Young LE, Marlin DJ,
*et al*.:**Cardiac output measured by lithium dilution, thermodilution and tran-soesophageal Doppler echocardiography in anaesthetised horses.***Am J Vet Res*2000,**61:**731–737. 10.2460/ajvr.2000.61.731View ArticleGoogle Scholar - Karamanoglu M, Bennett T, Ståhlberg M,
*et al*.:**Estimation of cardiac output in patients with congestive heart failure by analysis of right ventricular pressure waveforms.***Biomed Eng Online*2011,**10:**36. 10.1186/1475-925X-10-36View ArticleGoogle Scholar - Sundaresan A, Chase JG, Hann HE, Shaw GM:
**Cardiac output estimation using pulmonary mechanics in mechanically ventilated patients.***Biomed Eng Online*2010,**9:**80. 10.1186/1475-925X-9-80View ArticleGoogle Scholar - Payen D, Gayat E:
**Which general intensive care unit patients can benefit from placement of the pulmonary artery catheter?***Crit Care*2006,**10**(Suppl 3):S7. 10.1186/cc4925View ArticleGoogle Scholar - Nishikawa T, Doshi S:
**Errors in the measurement of cardiac output by thermodilution.***Can J Anaesth*1993,**40:**142–153. 10.1007/BF03011312View ArticleGoogle Scholar - Gawlikowski M:
**Model investigation on selected metrological parameters of circulatory system for hemodynamics diagnostics.***Ph.D. dissertation*2011. Dept. Elect. Eng., Bialystok Univ. of Tech., Bialystok PL.Google Scholar - Gawlikowski M, Pustelny T, Przywara-Chowaniec B, Nowak-Gawlikowska J:
**Model study of cardiac output measurement by thermodilution in thermal instability.***Acta Phys Pol A*2010,**118:**1124–1026.View ArticleGoogle Scholar - Gawlikowski M, Pustelny T, Urbanczyk M:
**The method of the average flow measurement in the circulatory system, especially the cardiac output measurement, Polish Patent application P-394270.**2011.Google Scholar - Beck MS, Plaskowski A:
*Cross correlation flowmeters - their design and application*. Hilger: Bristol; 1987.Google Scholar - Zator S:
*Cross-correlation measurement of liquids volumetric flow*. Kielce PL: Szumacher Publishers; 1997. (in Polish).Google Scholar - Grybos R:
*Fluid mechanics and hydraulics*. Silesian Technical University Publishers: Gliwice PL; 1977. (in Polish).Google Scholar - Walus S:
*Ultrasound flowmeters. Methodology of application*. Silesian Technical University Publishers: Gliwice PL; 1997. (in Polish).Google Scholar - Pustelny T, Struk P, Nawrat Z, Gawlikowski M:
**Design and numerical analyses of the human greater circulatory system.***Eur Phys J Spec Top*2008,**154:**171–174. 10.1140/epjst/e2008-00539-8View ArticleGoogle Scholar - Dietrich CF:
*Uncertainty, calibration and probability. The statistics of scientific and industrial measurement*. Briston UK: A. Hilger; 1991.Google Scholar - Wilkinson JR, Ruff C, Patel SK:
**A comparison between oesophageal Doppler and continuous thermodilution for the measurement of cardiac output in critically ill patients.***Crit Care*2000,**4:**P9.View ArticleGoogle Scholar - Parisi AF,
*et al*.:**Approaches to determination of left ventricular volume and ejection fFraction by real-time two-dimensional echocardiography.***Clin Cardiol*1979,**2:**257–263. 10.1002/clc.4960020404View ArticleGoogle Scholar - Gawlikowski M, Pustelny T, Kustosz R, Struk P:
**The methods of physical parameters measurement regarding the heart supporting automation.***Eur Phys J Spec Top*2008,**154:**71–76. 10.1140/epjst/e2008-00519-0View ArticleGoogle Scholar - Gawlikowski M, Pustelny T, Kustosz R:
**The physical parameters estimation of physiologically worked heard prosthesis.***J Phys IV*2006,**137:**73–78.Google Scholar - Gawlikowski M, Pustelny T, Kustosz R:
**Selected problems of mechanical heart supporting automation.***Eur Phys J Spec Top*2008,**154:**65–69. 10.1140/epjst/e2008-00518-1View ArticleGoogle Scholar - Konieczny G, Opilski Z, Pustelny T, Maciak E:
**State of the work diagram of the artificial heart.***Acta Phys Pol A*2009,**116:**344–347.View ArticleGoogle Scholar

## Copyright

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.