- Oral presentation
Accurate method for measuring arterial pulse wave velocity by cardiovascular magnetic resonance
Journal of Cardiovascular Magnetic Resonancevolume 14, Article number: O12 (2012)
A modified version of the flow-area (QA) method for measuring pulse wave velocity (PWV) from CMR images is introduced. The new ‘flow-time-area’ (QTA) method is more robust to errors in area measurements than the QA method, which is mathematically proved and demonstrated in pulmonary hypertension study. PWV was correlated to measurements from catheterization, echo, and CMR. The effect of image resolution, inter-, and intra-observer variabilities were studied. The results showed superiority of the QTA method over QA method.
The flow-area (QA) method has potential importance for measuring pulse wave velocity (PWV) from CMR images, especially in the pulmonary artery (PA). Nevertheless, the QA method is highly affected by errors in area measurements. In this work, we introduce the ‘flow-time-area’(QTA) method, a modified version of QA, which has benign behavior to area measurements inaccuracies.
In QTA, area and flow changes are calculated separately w.r.t. time, from which PWV is calculated by dividing the two slopes. The following sections show that QTA is more accurate than than QA.
In QA, flow measurements, q, are related to area measurements, a, by linear least-squares fit:q=ua+v, where PWVQA = u. By solving the normal equation⇒ PWVQA=(m∑aiqi-∑ai∑qi)/(m∑(ai)2-(∑ai)2), where m= #measurements. In QTA, two lines are formed: q=u1t+v1 and a=u2t+v2, from which⇒ PWVQTA=(m∑tiqi-∑ti∑qi)/(m∑tiai-∑ti∑ai). Under perfect linear relationship between q and a, PWVQA=PWVQTA=u. We analyze the case of positive error in area measurement (Δa) at point j. j is assumed >m/2, which reduces PWV and the accuracy measure M=PWV/ΔPWV-1, where ΔPWV is the difference between exact and measured PWV. For zero measurement error, M→∞. By mathematical manipulation, MQA=(N1+N2)/(D1+D2) and MQTA=N1/D1, where N1=u2(m∑(ti)2-(∑ti)2), N2=D1=(mtj-∑ti)Δa, and D2=(m-1)(Δa)2/u2. It can be shown that N1D2-N2D1=(ti-tk)2≥0 (summation over i>k; i,k≠j), which proves that MQA≤MQTA.
50 human subjects (25 PA hypertension (PAH) and 25 volunteers) were scanned on Siemens 3T scanner. Two sets of velocity-encoded images were acquired with resolutions 0.7mm and 1.4mm (Fig. 1). PA boundary was semi-automatically determined by two observers. PWV was calculated by QTA and QA. The measurements were compared using t-test. Normalized least-squares errors were calculated. Correlation analysis was conducted between PWV and mean PA pressure (mPAP), pulmonary vascular-resistance (PVR), and cardiac-index (CI). Bland-Altman analysis was conducted to measure inter- and intra-observer variabilities.
PWVQTA and PWVQA showed significant differences in PAH(p<.005), and less significant in volunteers (p=.01). The following results are for PAH: Normalized fitting-errors were 0.16/0.31for QTA/QA. Correlation coefficients between PWV versus mPAP, PVR, and CI were 0.93/0.61, 0.86/0.65, -0.71/-0.62 for QTA/QA. Standard deviations of differences between high- and low-resolution measurements were 0.19/0.41 m/s for QTA/QA. Correlation coefficients between low-resolution PWV and mPAP were 0.87/0.4 for QTA/QA. Bland-Altman showed minimal inter- and intra-observer variabilities.
The QTA method is more accurate than QA, because area measurement errors are not propagated to flow. The improved QTA accuracy is imperative when changes in artery cross-sectional area are minimal(e.g. PAH), when spatial resolution is low, or automatic vessel segmentation is implemented. Thus, QTA should be adopted for accurate PWV estimates.
James & Esther King Grant # 09KN-03-23138.