Accurate T1 mapping in patients with Pulmonary Hypertension and age matched volunteers using synthetic image based registration
Laura Claire Saunders1, Neil J Stewart1, Charlotte Hammerton1, David Capener1, Valentina O Puntmann2, David G Kiely3, Martin J Graves4, Andy Swift1, and Jim M Wild1

1Academic Unit of Radiology, The University of Sheffield, Sheffield, United Kingdom, 2Department of Cardiovascular Imaging, Kings College London, London, United Kingdom, 3The University of Sheffield, Sheffield, United Kingdom, 4University of Cambridge School of Clinical Medicine, University of Cambridge, Cambridge, United Kingdom

Synopsis

Patients with suspected pulmonary hypertension (n=94) and healthy volunteers (n=26) underwent T1 mapping of the right ventricle with a Modified Look Locker inversion recovery (MOLLI) sequence at 1.5T. MOLLI images were registered using pairwise registration to synthetic images produced using a simplified inversion recovery model to correct cardiac or respiratory motion. 89% of patients and 100% of healthy volunteers were successfully registered, with mean T1s of 1.00±0.10s and 0.97±0.06s (septal), 1.05±0.11s and 0.97±0.06s (right ventricular insertion point) and 1.02±0.11s and 1.04±0.13s (right ventricular free wall) respectively.

Background

Pulmonary hypertension (PH) is a condition characterised by increased blood pressure in the pulmonary artery, with a diagnostic criteria of mean pulmonary artery pressure ≥25mmHg as measured by right heart catheter. Patients with pulmonary hypertension may experience shortness of breath, leading to difficulty maintain breath hold during a CMR (cardiac magnetic resonance) scan.

Right ventricular remodelling and failure occur as a result of prolonged afterload. Accurate measurement of the changes in RV free wall tissue would be of clinical value. T1 elevation has been measured in animal models with pulmonary hypertension in the right ventricular insertion point and linked to severity of disease [1] In comparison to the in the septum and insertion points the right ventricular free wall is less spatially resolved on an MRI and it is unknown whether accurate measurement of the T1 can be made. It is therefore of interest to perform accurate T1 mapping in the right ventricular free wall and insertion points, in patients with pulmonary hypertension in comparison with age and sex matched volunteers.

Purpose

The aim of this study was to measure accurate myocardial T1, in the right ventricular free wall, septum and insertion points in patients with pulmonary hypertension and age and sex matched volunteers using a custom image registration method to correct for respiratory and cardiac motion.

Methods

All patients and volunteers underwent MOLLI imaging on a 1.5 T GE scanner, using a 3-3-5 image acquisition pattern and a flip angle of 35o. 94 consecutive patients referred from a PH referral centre with suspected PH were analysed. 26 age and sex matched volunteers were recruited for this study.

All acquired MOLLI images were registered using pair-wise image registration to synthetic images, based upon the work of Xue et al [2]. However, in this work synthetic images were created using a non-iterative simplified inversion recovery model (see Eq. 1) with an input of three pre-registered acquired images: the first image, the fourth image and the last image. The fourth image is included based on empirical optimisation to provide a good model fitting at the myocardial null point. The model assumes: the last image (image l) is fully recovered; that all signal is fully inverted, and that T1 is equivalent to T1*. This method is less time and computation intensive than an iterative energy minimisation solution, whilst still producing spatially aligned synthetic images with similar contrast to acquired images, enabling registration. See figure 1.

Image registration was performed using a non-rigid, intensity based, registration toolkit ShIRT [3] which was integrated into the custom T1 mapping algorithm in Matlab. $$ s(t)=s(l)-2s(l)-e^{\frac{-t}{T1}} Eq 1 $$
Acquired MOLLI images were overlaid pre and post registration to demonstrate areas of motion correction. T1 maps were produced pre and post registration to highlight areas of T1 correction due to motion correction. See Figure 2. T1 analysis regions of interest (see Figure 2) covered the entire septum, and as much RVFW as was visible. See Figure 3.

Results

84/94 of all acquired MOLLI images from patients with suspected PH were free from motion post registration, of which 83/94 were also free from artefact and therefore able to produce accurate T1 maps, and were used in the results section. See Figure 4.

Right ventricular free wall was identifiable in 73/83 patients, and 22/26 healthy volunteers.

T1s in the myocardium of patients were found to be: septal 1.00±0.10s; right ventricular insertion point 1.05±0.11s; right ventricular free wall 1.02±0.11s).

Healthy volunteer myocardial T1’s were consistent with known values: septal T1 = 0.97±0.06s; right ventricular insertion point T1 =0.97±0.06s [4], right ventricular free wall was found to be 1.04±0.13s.

Discussion and conclusion

Image registration to synthetic images, using a simplified inversion recovery model, created accurate (motion free) T1 maps were in 90% of all patients with pulmonary hypertension and 100% of healthy volunteers.

Right ventricular free wall was partially to fully spatially resolved in 88% of patients and 85% of healthy volunteers. Right ventricular free wall T1 was measured to be 1.02±0.11s and 1.04±0.13s in patients and volunteers respectively. The elevated T1 in volunteer right ventricular free wall, compared to septal T1, indicates partial volume effects within the right ventricular free wall measurements, and therefore inaccurate RVFW T1 values.

The use of image registration can produce accurate T1 maps in both healthy volunteers, and patients with pulmonary hypertension who have difficulty maintaining breath hold. T1 values measured healthy volunteers within the septum and right ventricular insertion point are consistent with known values.

Acknowledgements

This work was funded by the National Institute of Health Research (NIHR).

References

1. Garcia-Alvarez, A., et al., Association of Myocardial T1-Mapping CMR With Hemodynamics and RV Performance in Pulmonary Hypertension. Jacc-Cardiovascular Imaging, 2015. 8(1): p. 76-82.
2. Xue, H., et al., Motion correction for myocardial T1 mapping using image registration with synthetic image estimation. Magn Reson Med, 2012. 67(6): p. 1644-55.
3. Barber, D.C. and D.R. Hose, Automatic segmentation of medical images using image registration: diagnostic and simulation applications. J Med Eng Technol, 2005. 29(2): p. 53-63.
4. Rogers, T., et al., Standardization of T1 measurements with MOLLI in differentiation between health and disease--the ConSept study. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance, 2013. 15: p. 78.

Figures

Work flow of image analysis process.

Pre and post registration acquired images were overlaid (green: pre-registration, pink: post-registration). T1 maps were produced pre-registration (A) and View Figure 1 post-registration (B). Resultant T1 maps were overlaid to show areas of T1 correction due to registration (C, D).


T1 maps, post registration, of a healthy volunteer (A) and patient with PH (B).

Demonstration of region of interest selection.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0460