MRI-based transfer function determination through the transfer matrix by jointly fitting the incident and scattered B1+ field

A workflow for predicting radiofrequency-induced heating around bilateral deep brain stimulation electrodes in MRI

BACKGROUND

Heating around deep brain stimulation (DBS) in magnetic resonance imaging (MRI) occurs when the time-varying electromagnetic (EM) fields induce currents in the electrodes which can generate heat and potentially cause tissue damage. Predicting the heating around the electrode contacts is important to ensure the safety of patients with DBS implants undergoing an MRI scan. We previously proposed a workflow to predict heating around DBS contacts and introduced a parameter, equivalent transimpedance, that is independent of electrode trajectories, termination, and radiofrequency (RF) excitations. The workflow performance was validated in a unilateral DBS system.

PURPOSE

To predict RF heating around the contacts of bilateral (DBS) electrodes during an MRI scan in an anthropomorphic head phantom.

METHODS

Bilateral electrodes were fixed in a skull phantom filled with hydroxyethyl cellulose (HEC) gel. The electrode shafts were suspended extracranially, in a head and torso phantom filled with the same gel material. The current induced on the electrode shaft was experimentally measured using an MR-based technique 3 cm above the tip. A transimpedance value determined in a previous offline calibration was used to scale the shaft current and calculate the contact voltage. The voltage was assigned as a boundary condition on the electrical contacts of the electrode in a quasi-static (EM) simulation. The resulting specific absorption rate (SAR) distribution became the input for a transient thermal simulation and was used to predict the heating around the contacts. RF heating experiments were performed for eight different lead trajectories using circularly polarized (CP) excitation and two linear excitations for one trajectory. The measured temperatures for all experiments were compared with the simulated temperatures and the root-mean-squared errors (RMSE) were calculated.

RESULTS

The RF heating around the contacts of both bilateral electrodes was predicted with ≤ 0.29°C of RMSE for 20 heating scenarios.

CONCLUSION

The workflow successfully predicted the heating for different bilateral DBS trajectories and excitation patterns in an anthropomorphic head phantom.

Towards an integrated radiofrequency safety concept for implant carriers in MRI based on sensor-equipped implants and parallel transmission

To protect implant carriers in MRI from excessive radiofrequency (RF) heating it has previously been suggested to assess that hazard via sensors on the implant. Other work recommended parallel transmission (pTx) to actively mitigate implant-related heating. Here, both ideas are integrated into one comprehensive safety concept where native pTx safety (without implant) is ensured by state-of-the-art field simulations and the implant-specific hazard is quantified in situ using physical sensors. The concept is demonstrated by electromagnetic simulations performed on a human voxel model with a simplified spinal-cord implant in an eight-channel pTx body coil at 3 T . To integrate implant and native safety, the sensor signal must be calibrated in terms of an established safety metric (e.g., specific absorption rate [SAR]). Virtual experiments show that E -field and implant-current sensors are well suited for this purpose, while temperature sensors require some caution, and B 1 probes are inadequate. Based on an implant sensor matrix Q s , constructed in situ from sensor readings, and precomputed native SAR limits, a vector space of safe RF excitations is determined where both global (native) and local (implant-related) safety requirements are satisfied. Within this safe-excitation subspace, the solution with the best image quality in terms of B 1 + magnitude and homogeneity is then found by a straightforward optimization algorithm. In the investigated example, the optimized pTx shim provides a 3-fold higher mean B 1 + magnitude compared with circularly polarized excitation for a maximum implant-related temperature increase ∆ T imp ≤ 1 K . To date, sensor-equipped implants interfaced to a pTx scanner exist as demonstrator items in research labs, but commercial devices are not yet within sight. This paper aims to demonstrate the significant benefits of such an approach and how this could impact implant-related RF safety in MRI. Today, the responsibility for safe implant scanning lies with the implant manufacturer and the MRI operator; within the sensor concept, the MRI manufacturer would assume much of the operator's current responsibility.

Simplified modeling of implanted medical devices with metallic filamentary closed loops exposed to low or medium frequency magnetic fields

BACKGROUND AND OBJECTIVES

Electric currents are induced in implanted medical devices with metallic filamentary closed loops (e.g., fixation grids, stents) when exposed to time varying magnetic fields, as those generated during certain diagnostic and therapeutic biomedical treatments. A simplified methodology to efficiently compute these currents, to estimate the altered electromagnetic field distribution in the biological tissues and to assess the consequent biological effects is proposed for low or medium frequency fields.

METHODS

The proposed methodology is based on decoupling the handling of the filamentary wire and the anatomical body. To do this, a circuital solution is adopted to study the metallic filamentary implant and this solution is inserted in the electromagnetic field solution involving the biological tissues. The Joule losses computed in the implant are then used as a forcing term for the thermal problem defined by the bioheat Pennes' equation. The methodology is validated against a model problem, where a reference solution is available.

RESULTS

The proposed simplified methodology is proved to be in good agreement with solutions provided by alternative approaches. In particular, errors in the amplitude of the currents induced in the wires result to be always lower than 3%. After the validation, the methodology is applied to check the interactions between the magnetic field generated by different biomedical devices and a skull grid, which represents a complex filamentary wire implant.

CONCLUSIONS

The proposed simplified methodology, suitable to be applied to closed loop wires in the low to intermediate frequency range, is found to be sufficiently accurate and easy to apply in realistic exposure scenarios. This modeling tool allows analyzing different types of small implants, from coronary and biliary duct stents to orthopedic grids, under a variety of exposure scenarios.

MRI-based transfer function determination through the transfer matrix by jointly fitting the incident and scattered B1+ field

PURPOSE

A purely experimental method for MRI-based transfer function (TF) determination is presented. A TF characterizes the potential for radiofrequency heating of a linear implant by relating the incident tangential electric field to a scattered electric field at its tip. We utilize the previously introduced transfer matrix (TM) to determine transfer functions solely from the MR measurable quantities, that is, the B 1 + and transceive phase distributions. This technique can extend the current practice of phantom-based TF assessment with dedicated experimental setup toward MR-based methods that have the potential to assess the TF in more realistic situations.

THEORY AND METHODS

An analytical description of the B 1 + magnitude and transceive phase distribution around a wire-like implant was derived based on the TM. In this model, the background field is described using a superposition of spherical and cylindrical harmonics while the transfer matrix is parameterized using a previously introduced attenuated wave model. This analytical description can be used to estimate the transfer matrix and transfer function based on the measured B 1 + distribution.

RESULTS

The TF was successfully determined for 2 mock-up implants: a 20-cm bare copper wire and a 20-cm insulated copper wire with 10 mm of insulation stripped at both endings in respectively 4 and 3 different trajectories. The measured TFs show a strong correlation with a reference determined from simulations and between the separate experiments with correlation coefficients above 0.96 between all TFs. Compared to the simulated TF, the maximum deviation in the estimated tip field is 9.4% and 12.2% for the bare and insulated wire, respectively.

CONCLUSIONS

A method has been developed to measure the TF of medical implants using MRI experiments. Jointly fitting the incident and scattered B 1 + distributions with an analytical description based on the transfer matrix enables accurate determination of the TF of 2 test implants. The presented method no longer needs input from simulated data and can therefore, in principle, be used to measure TF's in test animals or corpses.