Theo J.C. Faes

GREIT: A unified approach to 2D linear EIT reconstruction of lung images

Electrical impedance tomography (EIT) is an attractive method for clinically monitoring patients during mechanical ventilation, because it can provide a non-invasive continuous image of pulmonary impedance which indicates the distribution of ventilation. However, most clinical and physiological research in lung EIT is done using older and proprietary algorithms; this is an obstacle to interpretation of EIT images because the reconstructed images are not well characterized. To address this issue, we develop a consensus linear reconstruction algorithm for lung EIT, called GREIT (Graz consensus Reconstruction algorithm for EIT). This paper describes the unified approach to linear image reconstruction developed for GREIT. The framework for the linear reconstruction algorithm consists of (1) detailed finite element models of a representative adult and neonatal thorax, (2) consensus on the performance figures of merit for EIT image reconstruction and (3) a systematic approach to optimize a linear reconstruction matrix to desired performance measures. Consensus figures of merit, in order of importance, are (a) uniform amplitude response, (b) small and uniform position error, (c) small ringing artefacts, (d) uniform resolution, (e) limited shape deformation and (f) high resolution. Such figures of merit must be attained while maintaining small noise amplification and small sensitivity to electrode and boundary movement. This approach represents the consensus of a large and representative group of experts in EIT algorithm design and clinical applications for pulmonary monitoring. All software and data to implement and test the algorithm have been made available under an open source license which allows free research and commercial use.

Pulmonary vascular resistance and compliance stay inversely related during treatment of pulmonary hypertension

AIMS Pulmonary arterial compliance (C) is increasingly being recognized as an important contributor to right ventricular afterload, but for monitoring of treatment of pulmonary hypertension (PH) most often still only pulmonary vascular resistance (R) is used. We aimed at testing the hypothesis that R and C are coupled during treatment of PH and that substantial changes in both R and C would result in more haemodynamic improvement than changes in R alone. METHODS AND RESULTS Data were analysed of two right-heart catheterizations of 52 patients with pulmonary arterial hypertension and 10 with chronic-thromboembolic PH. The product of R and C (= stroke volume over pulse pressure) did not change during therapy (P = 0.320), implying an inverse relationship. Changes in cardiac index correlated significantly (P < 0.001) with changes in R (R(2) = 0.37), better with changes in C (R(2) = 0.66), and best with changes in both (R(2) = 0.74). CONCLUSION During therapy for PH, R and C remain inversely related. Therefore, changes in both R and C better explain changes in cardiac index than either of them alone. Not only resistance but also compliance plays a prominent role in PH especially in an early stage of the disease.

Correction of phase offset errors in main pulmonary artery flow quantification

To investigate whether an existing method for correction of phase offset errors in phase‐contrast velocity quantification is applicable for assessment of main pulmonary artery flow with an MR scanner equipped with a high‐power gradient system.

Modeling the Instantaneous Pressure–Volume Relation of the Left Ventricle: A Comparison of Six Models

Simulations are useful to study the heart’s ability to generate flow and the interaction between contractility and loading conditions. The left ventricular pressure–volume (PV) relation has been shown to be nonlinear, but it is unknown whether a linear model is accurate enough for simulations. Six models were fitted to the PV-data measured in five sheep and the estimated parameters were used to simulate PV-loops. Simulated and measured PV-loops were compared with the Akaike information criterion (AIC) and the Hamming distance, a measure for geometric shape similarity. The compared models were: a time-varying elastance model with fixed volume intercept (LinFix); a time-varying elastance model with varying volume intercept (LinFree); a Langewouter’s pressure-dependent elasticity model (Langew); a sigmoidal model (Sigm); a time-varying elastance model with a systolic flow-dependent resistance (Shroff) and a model with a linear systolic and an exponential diastolic relation (Burkh). Overall, the best model is LinFree (lowest AIC), closely followed by Langew. The remaining models rank: Sigm, Shroff, LinFix and Burkh. If only the shape of the PV-loops is important, all models perform nearly identically (Hamming distance between 20 and 23%). For realistic simulation of the instantaneous PV-relation a linear model suffices.

Cardiac phase-dependent time normalization reduces load dependence of time-varying elastance

The time-varying elastance concept provides a comprehensive description of the intrinsic mechanical properties of the left ventricle that are assumed to be load independent. Based on pressure-volume measurements obtained with combined pressure conductance catheterization in six open-chest anesthetized sheep, we show that the time to reach end systole (defined as maximal elastance) is progressively prolonged for increasing ventricle pressures, which challenges the original (load-independent) time-varying elastance concept. Therefore, we developed a method that takes into account load dependency by normalization of time course of the four cardiac phases (isovolumic contraction, ejection, isovolumic relaxation, filling) individually. With this normalization, isophase lines are obtained that connect points in pressure-volume loops of different beats at the same relative time in each of the four cardiac phases, instead of isochrones that share points at the same time in a cardiac cycle. The results demonstrate that pressure curves can be predicted with higher accuracy, if elastance curves are estimated using isophase lines instead of using isochrones [root-mean-square error (RMSE): 3.8 +/- 1.0 vs. 14.0 +/- 7.4 mmHg (P < 0.001), and variance accounted for (VAF): 94.8 +/- 1.3 vs. 78.6 +/- 14.8% (P < 0.001)]. Similar results were found when the intercept volume was assumed to be time varying [RMSE: 1.7 +/- 0.3 vs. 13.4 +/- 7.4 mmHg (P < 0.001), and VAF: 97.4 +/- 0.5 vs. 81.8 +/- 15.5% (P < 0.001)]. In conclusion, phase-dependent time normalization reduces cardiac load dependency of timing and increases accuracy in estimating time-varying elastance.

Evaluation of model-independent deconvolution techniques to estimate blood perfusion

This report evaluates several methods to estimate blood perfusion and residue functions in dynamic contrast enhanced (DCE) MRI. Among these are model-dependent and model-independent techniques. All methods were applied to series of Monte Carlo simulations to evaluate the accuracy in order to reproduce different underlying vascular residue functions and blood perfusions. Of the model-independent approaches the use of B-splines with Tikhonov regularization was shown to have a reasonable accuracy in blood perfusion estimations and was less biased than all model-dependent approaches. This technique seems most promising for application to experimental data.

Vascular narrowing in pulmonary arterial hypertension is heterogeneous: rethinking resistance

In idiopathic pulmonary arterial hypertension (PAH), increased pulmonary vascular resistance is associated with structural narrowing of small (resistance) vessels and increased vascular tone. Current information on pulmonary vascular remodeling is mostly limited to averaged increases in wall thickness, but information on number of vessels affected and internal diameter decreases for vessels of different sizes is limited. Our aim was to quantify numbers of affected vessels and their internal diameter decrease for differently sized vessels in PAH in comparison with non‐PAH patients. Internal and external diameters of transversally cut vessels were measured in five control subjects and six PAH patients. Resistance vessels were classified in Strahler orders, internal diameters 13 μm (order 1) to 500 μm (order 8). The number fraction, that is, percentage of affected vessels, and the internal diameter fraction, that is, percentage diameter of normal diameter, were calculated. In PAH, not all resistance vessels are affected. The number fraction is about 30%, that is, 70% of vessels have diameters not different from vessels of control subjects. Within each order, the decrease in diameter of affected vessels is variable with an averaged diameter fraction of 50–70%. Narrowing of resistance vessels is heterogeneous: not all vessels are narrowed, and the decrease in internal diameters, even within a single order, vary largely. This heterogeneous narrowing alone cannot explain the large resistance increase in PAH. We suggest that rarefaction could be an important contributor to the hemodynamic changes.