Percentage of change in the calculated SAR values in human head during 3T MRI of patients with deep brain stimulation implants: A computational study of realistic vs. simplified lead trajectories
Laleh Golestanirad1, Maria Ida Iacono2, Leonardo M Angelone2, and Giorgio Bonmassar1

1Radiology, Massachusetts General Hospital, Charlestown, MA, United States, 2Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, US Food and Drug Administration, Silver Spring, MD, United States

Synopsis

Each year approximately 300,000 patients with medical implants including deep brain stimulation (DBS) devices are denied magnetic resonance imaging (MRI) examination due to safety concerns. One of the major contraindications of MRI for DBS patient population is due to the potential for permanent injuries from excessive tissue heating. One open question when evaluating RF-induced heating with DBS is the effect of the lead path and the need for patient-specific information. Using finite element method, we report results of calculated SAR maps for patient-specific lead paths based on CT images, and compare them to simplified path trajectories.

Introduction

Each year approximately 300,000 patients with medical implants including deep brain stimulation (DBS) devices are denied magnetic resonance imaging (MRI) examination due to safety concerns 1 . One of the major contraindications of MRI for DBS imaging is due to the potential for permanent injuries from excessive tissue heating. Consequently, post-operative MRI of DBS patients is currently approved by the FDA with specific and limited conditions of use (e.g., 1.5T only) 2. One open question when evaluating RF-induced heating with DBS is the effect of the lead path and the need for patient-specific information. We report here, preliminary results of SAR maps using patient-specific lead paths extracted from post-operative CT images, and compare them to simplified lead trajectories.

Methods

Finite Element Time Domain technique (FDTD) and Finite Element Method (FEM) were used to calculate local SAR values around the DBS implants at 3T MRI with a transmit head coil in two models based on patient-specific data (see Fig.1). Realistic DBS lead trajectories were extracted from post-operative CT images, and were used to construct detailed numerical models of DBS lead trajectories, and to calculate 1g- and 10g-averaged SAR. Electrically homogeneous and electrically heterogeneous head models were used for the study. Case A: FDTD simulations were performed on a 1x1x1 mm MRI-based head model. The acquisition, processing and segmentation of the head model are described in a previous work 3. Two bilateral implants were modeled. Implants were manually segmented from the CT image and adjusted onto the multi-resolution head model by applying a transformation obtained from non-linear registration (CT image to head model). Each implant was modeled by a polyline containing 500 wire segments. The first segment of wire started from subthalamic nucleus (STN) and continued up to the skull. The modeled wire was insulated except for the 2 mm in the first segment contacting the STN, to model the electrical contact. Case B: FEM simulations were performed on an electrically homogeneous head model and a uni-lateral DBS implant, built from post-operative CT images of another patient. The lead trajectory was manually segmented from CT images and a detailed model of DBS lead composed of four electrode contacts, connected through a solid core and surrounded by hollow insulation, was constructed similar to those described in previous studies 4. For both cases, we also constructed simplified DBS lead models in which loops where eliminated from lead trajectories. Simulation results were normalized to produce a total head averaged SAR of 3.2 W/kg in the head model without the implant. For both studies, B1+ fields were calculated as B1+=0.5(B1x+jB1y) and SAR ratios between the simplified models and realistic models were defined as SAR_Ratio=10log10(Simplified SAR/Realistic SAR).

Results

B1+ fields, averaged over whole head without the implant, were comparable in both studies (3.5 µT in Case A and 3.0 µT in Case B). In both cases, we found that the simplified path, without loops, led to a much higher local SAR near the lead. Specifically, the simplified lead model in Case B showed an increase of 1g-avg.SAR and 10g-avg. SAR of 84% compared to the realistic path. The simplified lead model in Case A overestimated the peak of 1g-avg. SAR and 10g-avg. SAR by approximately 60%.

Conclusion and Future Work

To the authors’ knowledge, this work is the first attempt to evaluate quantitative values of local SAR in the brain tissue based on patient-derived realistic models of DBS lead geometries. Our results have two significant implications: Firstly, if verified in larger patient cohorts, similar studies might justify the use of simplified lead trajectories in numerical simulations that aim to assess safety of MR Conditional devices. Second, we found a significant difference between SAR values predicted from simplified models and those predicted from realistic DBS paths. Our preliminary results suggest that actual SAR values during MRI of DBS patients might be less than those estimated by simple models. However, additional validation of these results is needed. We are now in the process of evaluating a cohort of 20 DBS patients to complement this study.

Acknowledgements

Preparation of this work was supported by grants 1R21EY020961-01, 1R43NS071988-01A, and 15P41RR014075-13 from the National Institutes of Health and, in part, by grants from NIDA 1R01DA027804-01 and NIMH 1R21MH084041-01A1.

References

1. Shellock FG, Spinazzi A. MRI safety update 2008: part 2, screening patients for MRI. American Journal of Roentgenology 2008;191(4):1140-1149.

2. Larson PS, Richardson RM, Starr PA, Martin AJ. Magnetic resonance imaging of implanted deep brain stimulators: experience in a large series. Stereotactic and functional neurosurgery 2008;86(2):92-100.

3. Makris N, Angelone L, Tulloch S, Sorg S, Kaiser J, Kennedy D, Bonmassar G. MRI-based anatomical model of the human head for specific absorption rate mapping. Medical & biological engineering & computing 2008;46(12):1239-1251.

4. Elwassif MM, Kong Q, Vazquez M, Bikson M. Bio-heat transfer model of deep brain stimulation-induced temperature changes. Journal of neural engineering 2006;3(4):306.

Figures

Fig.1: Top row: post-operative CT (Case A), Numerical model with realistic (left) and simplified (right) implant path.

Bottom row: post-operative CT (Case B), Numerical model with realistic implant (left) and simplified path.


Fig.2: Distribution of SAR_Ratio for 1g-averaged SAR. Top row: SAR_Ratio in the homogenous head model (Case B).

Bottom row: SAR_Ratio in the heterogeneous head model (Case A).


Table1: SAR values for realistic, simplified and no-lead cases.



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