Preliminary Results: Cardiac Cine “Watermark” MRI provides both Anatomical Function via Magnitude Cine and 2D Myocardial Strain via Spatially Modulated Phase
Ronald J Beyers1, Davis M Vigneault2, Dean Schwartz3, Nouha Salibi1,4, David A Bluemke2, and Thomas Denney1

1MRI Research Center, Auburn University, Auburn University, AL, United States, 2Radiology and Imaging Sciences, National Institutes of Health, Bethesda, MD, United States, 3Anatomy, Physiology and Pharmacology, Auburn University, Auburn University, AL, United States, 4MR R&D, Siemens Healthcare, Malvern, PA, United States

Synopsis

We developed a Cine Watermark (CWM) cine sequence that produces normal cine magnitude images, plus a grid pattern of tags added only in the phase for quantitative cine strain, while requiring no extra operator effort. Using spatial cosine modulation combined with k-space sum/differencing produced separate normal magnitude cine and unique phase-only grid-tags for strain calculation. In vivo rat and human scans demonstrated good magnitude cine and phase-only quantified displacement. Calculated by Farneback optical flow algorithm, the peak principle strain, averaged around the LV for rat = -16.5±2.4 % and human = -17.8±6.2 % (mean±StdDev).

Purpose

Myocardial tagging, introduced over 25 years ago, is universally implemented on all manufacturers’ MRI scanners1-3. Despite extensive validation as the best method for clinical quantitative measure of myocardial wall deformation, tagging is limited to a research context due to the time intensive analysis and need to acquire additional sequences. Cardiac MRI (CMR) sequences must provide maximum diagnostic information, and maximum patient throughput, with minimum operator effort and time (cost)4. We developed a CMR cine sequence with an added “hidden” capability to acquire cardiac strain data while requiring no extra operator effort. Called Cine Watermark (CWM), this method produces normal (qualitative) cine magnitude images, plus a grid pattern of tags added only in the phase (watermark) information for quantitative cine strain capability.

Methods

Fig 1 presents the ECG-triggered CWM sequence that acquires both normal and complimentary spatial modulation of magnetization (SPAMM/CSPAMM) prepared cine image sets. A classic “1-1” saturation SPAMM was used were each non-selective Rect RF pulse had a flip angle of 45°. The sequence alternated the acquisition of SPAMM or CSPAMM on every other ECG trigger. The cine readout used a standard, low flip angle, spoiled GRE method. CWM acquisition runs were acquired with the cosine modulation in the readout (RO) axis and then phase encode (PE) axis directions to create a grid-tag pattern. Two rats were scanned in a 7T human-size scanner with a birdcage 1-chan small animal RF coil and two healthy human volunteers, with informed consent, were scanned in a 3T Verio scanner with a 32-chan anterior/posterior RF coil array (both scanners: Siemens, Erlangen, Germany). All CWM scan slice locations were mid-left ventricular (LV), short-axis. Scan parameters for rats (humans): field-of-view = 64x40(264x198) mm, slice thickness = 2(8) mm, matrix = 256x160(176x132), in-plane pixel size = 0.25x0.25(1.5x1.5) mm, flip angle = 8(10)°, BW = 488(406) Hz/pix, averages = 4(1), cine phases = 15(20), tag spacing = 1.5(5) mm, ECG-triggered, free-breathing(11 breath holds). All CWM image reconstruction was performed offline using customized Matlab programs (Mathworks, Natick, MA). During image reconstruction, the sum of the SPAMM + CSPAMM raw k-space data eliminated the cosine modulation components leaving only the conventional (central k-space) T1-weighted cine image component. The magnitude images of this sum’s inverse Fourier transform provided the pure cardiac cine image set without tag pattern lines. Conversely, the difference of SPAMM – CSPAMM raw k-space eliminated the central k-space T1-weighted component leaving only the cosine modulation components. The phase images of this difference’s inverse Fourier transform provided strong and persistent (0°/180°) phase-only tag lines that tracked cardiac displacement and allowed subsequent strain calculation. The myocardial strain analysis was implemented in C++ using the ITK5 and OpenCV6 libraries. Magnitude and phase images were weighted to create synthetic tag images. Dense displacement fields were obtained between successive frames using the Farneback optical flow algorithm (5 pyramids, 7px window). Principal strain values were then calculated pixelwise for each image, and the median value within each anatomical segment was calculated. A weighted sum of forward and reverse registrations was calculated to avoid propagation of error.

Results

Both rat and human CWM scans gave cine magnitude and displacement cine phase image sets suitable for analysis. Fig 2 shows the in vivo rat SPAMM/CSPAMM k-space sum effectively isolated the central k-space component (magnitude cine) and the k-space difference isolated the simultaneous RO and PE axes cosine modulation components (for strain). Fig 3 presents example rat CWM images at end-diastole (ED) and end-systole (ES) time points. The rat CWM exhibited some signal loss in the myocardium that is common at 7T due to susceptibility B0 inhomogeneities. Fig 4 presents similar human CWM ED and ES image results. Both Figs 3&4 show observable displacement of the phase grid-tag patterns evident at ES compared to the ED starting point, with regionalized displacement (quiver maps) displayed at far-right. The calculated peak principle strain, averaged around the LV, for rat = -16.5±2.4 % and human = -17.8±6.2 % (mean±StdDev).

Discussion

These preliminary results provided an in vivo demonstration that CWM can be used to obtain anatomical cine and tagged phase images in a single scan, although with some short-comings that need correction and optimization. The phase-tag displacements are much more obvious in the phase cine loop playback (not shown). Proper strain calculation was problematic and attributed to be from phase tags discontinuities. This CWM Encode method is likely compatible with SSFP cine and such an integration should be investigated.

Conclusions

We developed and demonstrated a CMR cine sequence with added “hidden” capability to acquire strain data. Future development effort will improve its performance.

Acknowledgements

No acknowledgement found.

References

1. Bolster BD, et al. “Myocardial tagging in polar coordinates with striped tags”, Radiology 1990;177:769–72.

2. Axel L, Dougherty L. “MR imaging of motion with spatial modulation of magnetization”, Radiology.1989;171:841–9.

3. Axel L, Dougherty L. “Heart wall motion: Improved method of spatial modulation…”, Radiology 1989;172:349–50.

4. Edelstein WA, et al. “MRI: Time Is Dose—and Money and Versatility, J Am Coll Radiol 2010;7(8): 650–2.

5. Insight Segmentation and Registration Toolkit (ITK), www.itk.org

6. Open Source Computer Vision (OpenCV), www.opencv.org

Figures

Fig 1 - Cine Watermark sequence is ECG-triggered and captures an entire RR cardiac period. In the CWM Encode module, the second Rect RF flip angle is set +45° for normal SPAMM and -45° for CSPAMM.

Fig 2 - Example in vivo raw k-space magnitude images showing that SPAMM + CSPAMM leaves only the central k-space T1-weighted component and that SPAMM – CSPAMM leaves only the orthogonal cosine modulation tagging components.

Fig 3 - CWM rat result images with added segmentation lines over the left ventricle myocardium. The end systole (ES) phase image clearly displays tags pattern displacement that is further visually quantified by the green quiver vectors in the ES Displacement image.

Fig 4-CWM human result images with added segmentation lines over the left ventricle myocardium. Similar to Fig 3 rat results, the end systole (ES) phase image clearly displays tags pattern displacement that is further visually quantified by the green quiver vectors in the ES Displacement image.



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