Quantification of left ventricular deformation fields from undersampled radial, real-time cardiac MRI

Background Regional ventricular elasticity is altered in ischemic heart disease and LV hypertrophy, potentially contributing to and sustaining a ventricular remodeling process. Conventional methods to measure regional ventricular elasticity in preclinical animal models require high fidelity ventricular pressure measurement, sonomicrometry transducers and dynamic modulation of preload. These procedures are highly invasive and provide only limited regional information; patient studies are prohibitive. We dynamically measure global LV elastance using undersampled radial real-time MRI and seek to find regional LV elastance in these studies. This requires tracking LV motion derived from wall thickness measurements in a large number of real-time image frames at high accuracy. We assessed several motion tracking algorithms and characterize their fidelity with respect to true motion labeled by experts and sensitivity to residual radial streak artifact. Methods A series of short axis MR slices were acquired in 9 patients using a golden-angle radial bSSFP trajectory (total views = 4000, radial views per image = 34, shared profiles = 31, reconstructed frame rate = 98 fps) and reconstructed via conjugate gradient-SENSE based iterative reconstruction. Two independent researchers manually segmented the datasets and the movies were restricted to three full heartbeats to limit the number of images. Image registrations were performed between the initial end-diastolic frame and all subsequent frames. The image registrations were performed using b-spline and symmetric diffeomorphic (SyN) image registration across a range of spatial scales (shown in Table 1). Four separate image registration metrics were investigated (mutual information, meansquared overlap, cross-correlation, and demons). Although we expected mean squares and cross-correlation to perform better in this intramodality registration, radial streaking could potentially corrupt data integrity. The resulting deformation fields were evaluated using an array of metrics (dice coefficient, area overlap, false positive, and false negative rate) for all registrations, image frames and patients.