1506

Improved MR thermometry for laser-induced thermal therapy – tradeoffs between imaging approaches
Henrik Odéen1 and Dennis L Parker1

1Radiology and Imaging Sciences, University of Utah, Salt Lake City, UT, United States

Synopsis

MR thermometry is often used to monitor thermal therapies such as focused ultrasound and laser induced thermal therapies (LITT). As in MRI in general, there is an inherent tradeoff between measurement accuracy, precision, and spatial and temporal resolution in MR thermometry. In this work we present improved acquisition protocols for 2D and 3D MR thermometry for LITT applications. We investigate and compare image quality and temperature precision for 8 different 2D and 3D GRE, and 3D segmented EPI protocols. Experiments are performed in a healthy volunteer (non-heating) and tissue-mimicking gel (with heating).

Introduction

MRI has become the image modality of choice to monitor many non- and minimally-invasive interventional procedures such as focused ultrasound and laser-inducer thermal therapy (LITT) because it combines high-resolution anatomical imaging for procedure guidance with real-time MR temperature imaging (MRTI) for procedure monitoring. The proton resonance frequency shift (PRFS) method (1,2) is the most commonly used MRTI approach as it has been shown to work for a wide range of (aqueous) tissues and is linear over a relatively large temperature range (3). As in MR in general, there are some inherent trade-offs between field of view, spatial resolution, temporal resolution, and measurement precision in MRTI. In this work we compare different 2D and 3D approaches for performing MRTI for LITT applications in the brain using the PRFS approach. We investigate trade-offs in MRTI precision, image quality (distortion) and off-resonance effects for 2D and 3D gradient recalled echo (GRE), and 3D echo-planar imaging (EPI) approaches.

Methods

Table 1 lists the 8 protocol approaches investigated in this work. #1 is the 2DGRE protocol currently used clinically with our LITT system (Visualase™ MRI-Guided Laser Ablation System, Louisville, CO), where 2 orthogonal 2D slices (usually oriented along the laser fiber in Sagittal+Coronal planes) are being sampled. Using the same TR and BW (23 ms, 260 Hz/px), 3 mono- or 4 bi-polar contrasts can be acquired (#2-#3) on our scanner (80mT/m gradient, 200T/m/s slew rates). This increases the total sampling time from 3.85ms to 11.54 and 15.38ms, respectively, improving MRTI precision by sqrt(sampling time) (1.7-2.0x) with no scan time penalty. #4 is 2DGRE with increased BW (810Hz/px) to reduce off-resonance shift from temperature change/fat. #5 is 3DGRE acquiring 12 slices in the same scan time as the 2DGRE protocols. #6-#8 are 3D segmented-EPI with mono- (i.e., flyback) and bi-polar readouts and varying ETL to investigate differences in off-resonance effects and distortion.

Scans were performed on a 3T scanner (PrismaFIT, Siemens, Erlangen, Germany) with a 4-channel wrap-around/flex-coil. In-vivo scans were performed on one healthy volunteer under IRB-approved informed consent. LITT sonications were performed in a human-sized plastic skull (3B Scientific, Tucker, GA) filled with tissue-mimicking gellan gum (Fisher Scientific, Pittsburgh, PA) at 3W for 8 dynamic measurements (~47-50 s), and were repeated 3 times/protocol. All data was zero-filled interpolated (for 2D only in-plane, and for 3D in all 3 directions) to 0.5-mm voxel spacing to minimize partial volume effects. MRTI utilized referenceless reconstruction (4). Measured T1 (in vivo/phantom ~1000/~1800 ms) was used to calculate and apply the Ernst angle for all protocols.

Results

Figure 1 shows magnitude and temperature standard deviation (SD) through time maps from a central sagittal slice in the volunteer for all 8 protocols. Figure 2 shows two orthogonal views of temperature maps at time of maximum heating, and Figure 3 shows temperature vs. time plots for the hottest voxel (mean±SD of three heatings), both for the phantom experiment. Table 2 lists the precision of the measurements, defined as the SD through time of an (un-heated) background region.

Discussion and Conclusions

Table 2 shows that the best MRTI precision was achieved with multi-contrast 2D protocols, followed by 3DEPI, single-contrast 2DGRE, and last 3DGRE. This agrees well with theory and sampling time duration for the different protocols (Table 1). The precision of 3DGRE is negatively affected by the parallel imaging (resulting in at least sqrt(R) SNR decrease), needed to achieve similar scan time. The higher maximum temperature rise measured in 3D protocols compared to 2D protocols (Figure 3) is likely due to reduced partial-volume effects due to the ZFI in three-dimensions vs only two-dimensions for 2D protocols. All 2D protocols show good SNR magnitude images with minimal distortions. The high-BW protocol (#4) shows the least fat-water shift (~1 voxel) and the 4 contrast bi-polar protocol (#3) show the largest shift (~3 voxels) – seen comparing the skin-fat layer around the skull (RO direction is H->F for sagittal images). 3D protocols all have high RO-BW and experience the smallest fat-water shift (all <1 voxel), but EPI sequences have rather low PE-BW and experience larger shifts in this direction (A->P). This is most clearly seen in the flyback protocol which has the lowest PE-BW.

To maximize MRTI precision, the available TR should be filled with multiple contrasts. The already high RO-BW allows for bi-polar readout without large off-resonance shifts. 3D EPI imaging with good precision, albeit somewhat more distortions, can be achieved in the same scan time. 3D GRE achieves small distortions, but for similar precision scan time increases are needed.

Acknowledgements

This work was supported by Medtronic and the Mark H. Huntsman endowed chair.

References

1. De Poorter J. Noninvasive MRI thermometry with the proton resonance frequency method Study of susceptibility effects. Magn. Reson. Med. 1995;43:359–67.

2. Ishihara Y, Calderon A, Watanabe H, Okamoto K, Suzuki Y, Kuroda K. A precise and fast temperature mapping using water proton chemical shift. Magn. Reson. Med. 1995;34:814–23.

3. Rieke V, Butts Pauly K. MR thermometry. J. Magn. Reson. Imaging 2008;27:376–90. doi: 10.1002/jmri.21265.

4. Rieke V, Vigen KK, Sommer G, Daniel BL, Pauly JM, Butts K. Referenceless PRF shift thermometry. Magn. Reson. Med. 2004;51:1223–1231. doi: 10.1002/mrm.20090.

Figures

Table 1. Scan parameters. All protocols used a 240 mm in-plane field of view with 50% phase encoding, i.e., 256x128 (readout x phase encoding) in-plane matrix size, for voxel spacing of 0.94x1.88 mm. TR/TE=Repetition/Echo time, BW=bandwidth, RO/PE=Readout-/Phase-encoding direction, ETL=Echo train length, ES=Echo spacing PF=Partial Fourier, PE/SE=Phase-/Slice-encoding direction, PPI=Partial-parallel imaging, Tacq=Acquisition time.

Figure 1. In-vivo scans in healthy volunteer. a) – h) show sagittal magnitude (top row) and temperature standard deviation through time maps (bottom row) for the 8 protocols listed in Table 1, respectively. Multi-contrast 2D GRE approaches show the best precision, followed by 3D EPI and single-contrast 2D GRE.

Figure 2. Phantom MRTI maps. a) – h) show sagittal (top row) and coronal (bottom row) MRTI maps for the 8 protocols listed in Table 1, respectively. It should be noted that the 2D protocols acquire 2 orthogonal 2D slices, while the 3D protocols acquire a contiguous volume of 12 parallel slices. The in-plane FOV is cropped to better show details.

Figure 3. Temperature rise vs. time curves for phantom data. Shown is mean ± SD of three repeated and identical heatings; 3W for 8 dynamic measurements (corresponding to 47 - 50 s depending on protocol). a) shows the hottest voxel vs. time for the 4 2D protocols, and b) shows corresponding data for the 4 3D protocols.

Table 2. MRTI precision. The measured SD through time for in-vivo and phantom data is listed for all 8 protocols. For the in-vivo data an ROI in the high-SNR region in the posterior-head direction was used, and for the phantom experiment a background ROI which did not experience any heating was used to calculate the SD.

Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
1506