Fast conductivity imaging in magnetic resonance electrical impedance tomography (MREIT) for RF ablation monitoring

Bias correction for phase-based cr-MREPT using low resolution B1+ magnitude

ObjectiveFull-form Magnetic Resonance Electrical Properties Tomography (MREPT) requires bothB1+magnitude and phase information. SinceB1+phase can be obtained faster and with higher SNR compared toB1+magnitude, several phase-based methods have been developed for conductivity imaging. However, phase-based methods suffer from a concave bias due to the assumption that∇|B1+|is negligible in the ROI.ApproachIn this paper, we re-derive the central equation of phase-based cr-MREPT without assuming that∇|B1+|is negligible and thus propose a correction method directly integrated into the equation system.Main resultsProposed method successfully corrects the concave bias on both simulated and experimental data and significantly increases image quality.SignificanceThe proposed correction method depends on a very low-resolution|B1+|map, and therefore the imaging time does not increase significantly for obtainingB1+magnitude. Moreover, correction can be achieved using simulatedB1+magnitude, hence completely removing the additional imaging requirement.

Use of dielectric padding to eliminate low convective field artifact in cr-MREPT conductivity images

PURPOSE

Convection-reaction equation-based magnetic resonance electrical properties tomography (cr-MREPT) provides conductivity images that are boundary artifact-free and robust against noise. However, these images suffer from the low convective field (LCF) artifact. We propose to use dielectric pads to alter the transmit magnetic field (B₁ ⁺ ), shift the LCF region, and eliminate the LCF artifact.

METHODS

Computer simulations were conducted to analyze the effects of pad electrical properties, pad thickness, pad height, arc angle, and thickness of the pad-object gap. In 3T MR experiments, water pads and BaTiO₃ pads were used with agar-saline phantoms. Two data sets (e.g., with the pad located on the left or on the right of the object [phantom]) were acquired, and the corresponding linear systems were simultaneously solved to get LCF artifact-free conductivity images.

RESULTS

A pad needed to have 180° arc angle and the same height with the phantom for maximum benefit. Increasing the pad thickness and/or the relative permittivity of the pad increased the LCF shift, whereas excessive amounts of these parameters caused errors in conductivity reconstructions because the effect of neglected B_z terms became noticeable. Conductivity of the pad, on the other hand, had minimal effect on elimination of the LCF artifact. Combining 2 data sets (i.e., with 2 different dielectric pad positions) resulted in more accurate conductivity maps (low L² -errors) as opposed to no pad or single pad cases in experiments and simulations.

CONCLUSIONS

Using the proposed technique, LCF artifact is significantly removed, and the reconstructed conductivity values are improved.

Accelerating acquisition strategies for low-frequency conductivity imaging using MREIT

We sought to improve efficiency of magnetic resonance electrical impedance tomography data acquisition so that fast conductivity changes or electric field variations could be monitored. Undersampling of k-space was used to decrease acquisition times in spin-echo-based sequences by a factor of two. Full MREIT data were reconstructed using continuity assumptions and preliminary scans gathered without current. We found that phase data were reconstructed faithfully from undersampled data. Conductivity reconstructions of phantom data were also possible. Therefore, undersampled k-space methods can potentially be used to accelerate MREIT acquisition. This method could be an advantage in imaging real-time conductivity changes with MREIT.

Feasibility of conductivity imaging using subject eddy currents induced by switching of MRI gradients

PURPOSE

To investigate the feasibility of low-frequency conductivity imaging based on measuring the magnetic field due to subject eddy currents induced by switching of MRI z-gradients.

METHODS

We developed a simulation model for calculating subject eddy currents and the magnetic fields they generate (subject eddy fields). The inverse problem of obtaining conductivity distribution from subject eddy fields was formulated as a convection-reaction partial differential equation. For measuring subject eddy fields, a modified spin-echo pulse sequence was used to determine the contribution of subject eddy fields to MR phase images.

RESULTS

In the simulations, successful conductivity reconstructions were obtained by solving the derived convection-reaction equation, suggesting that the proposed reconstruction algorithm performs well under ideal conditions. However, the level of the calculated phase due to the subject eddy field in a representative object indicates that this phase is below the noise level and cannot be measured with an uncertainty sufficiently low for accurate conductivity reconstruction. Furthermore, some artifacts other than random noise were observed in the measured phases, which are discussed in relation to the effects of system imperfections during readout.

CONCLUSION

Low-frequency conductivity imaging does not seem feasible using basic pulse sequences such as spin-echo on a clinical MRI scanner. Magn Reson Med 77:1926-1937, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Hyperthermia treatment planning for cervical cancer patients based on electrical conductivity tissue properties acquired in vivo with EPT at 3 T MRI

Introduction The reliability of hyperthermia treatment planning (HTP) is strongly dependent on the accuracy of the electric properties of each tissue. The values currently used are mostly based on ex vivo measurements. In this study, in vivo conductivity of human muscle, bladder content and cervical tumours, acquired with magnetic resonance-based electric properties tomography (MR-EPT), are exploited to investigate the effect on HTP for cervical cancer patients. Methods Temperature-based optimisation of five different patients was performed using literature-based conductivity values yielding certain antenna settings, which are then used to compute the temperature distribution of the patient models with EPT-based conductivity values. Furthermore, the effects of altered bladder and muscle conductivity were studied separately. Finally, the temperature-based optimisation was performed with patient models based on EPT conductivity values. Results The tumour temperatures for all EPT-based dielectric patient models were lower compared to the optimal tumour temperatures based on literature values. The largest deviation was observed for patient 1 with ΔT90 = -1.37 °C. A negative impact was also observed when the treatment was optimised based on the EPT values. For four patients ΔT90 was less than 0.6 °C; for one patient it was 1.5 °C. Conclusions Electric conductivity values acquired by EPT are higher than commonly used from literature. This difference has a substantial impact on cervical tumour temperatures achieved during hyperthermia. A higher conductivity in the bladder and in the muscle tissue surrounding the tumour leads to higher power dissipation in the bladder and muscle, and therefore to lower tumour temperatures.

Gradient-based electrical conductivity imaging using MR phase

PURPOSE

To develop a fast, practically applicable, and boundary artifact free electrical conductivity imaging method that does not use transceive phase assumption, and that is more robust against the noise.

THEORY

Starting from the Maxwell's equations, a new electrical conductivity imaging method that is based solely on the MR transceive phase has been proposed. Different from the previous phase based electrical properties tomography (EPT) method, a new formulation was derived by including the gradients of the conductivity into the equations.

METHODS

The governing partial differential equation, which is in the form of a convection-reaction-diffusion equation, was solved using a three-dimensional finite-difference scheme. To evaluate the performance of the proposed method numerical simulations, phantom and in vivo human experiments have been conducted at 3T.

RESULTS

Simulation and experimental results of the proposed method and the conventional phase-based EPT method were illustrated to show the superiority of the proposed method over the conventional method, especially in the transition regions and under noisy data.

CONCLUSION

With the contributions of the proposed method to the phase-based EPT approach, a fast and reliable electrical conductivity imaging appears to be feasible, which is promising for clinical diagnoses and local SAR estimation. Magn Reson Med 77:137-150, 2017. © 2016 Wiley Periodicals, Inc.

Current density imaging sequence for monitoring current distribution during delivery of electric pulses in irreversible electroporation

Background

Electroporation is gaining its importance in everyday clinical practice of cancer treatment. For its success it is extremely important that coverage of the target tissue, i.e. treated tumor, with electric field is within the specified range. Therefore, an efficient tool for the electric field monitoring in the tumor during delivery of electroporation pulses is needed. The electric field can be reconstructed by the magnetic resonance electric impedance tomography method from current density distribution data. In this study, the use of current density imaging with MRI for monitoring current density distribution during delivery of irreversible electroporation pulses was demonstrated.

Methods

Using a modified single-shot RARE sequence, where four 3000 V and 100 μs long pulses were included at the start, current distribution between a pair of electrodes inserted in a liver tissue sample was imaged. Two repetitions of the sequence with phases of refocusing radiofrequency pulses 90° apart were needed to acquire one current density image. For each sample in total 45 current density images were acquired to follow a standard protocol for irreversible electroporation where 90 electric pulses are delivered at 1 Hz.

Results

Acquired current density images showed that the current density in the middle of the sample increased from first to last electric pulses by 60%, i.e. from 8 kA/m² to 13 kA/m² and that direction of the current path did not change with repeated electric pulses significantly.

Conclusions

The presented single-shot RARE-based current density imaging sequence was used successfully to image current distribution during delivery of short high-voltage electric pulses. The method has a potential to enable monitoring of tumor coverage by electric field during irreversible electroporation tissue ablation.

Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures

The application of supraphysiological temperatures (>40°C) to biological tissues causes changes at the molecular, cellular, and structural level, with corresponding changes in tissue function and in thermal, mechanical and dielectric tissue properties. This is particularly relevant for image-guided thermal treatments (e.g. hyperthermia and thermal ablation) delivering heat via focused ultrasound (FUS), radiofrequency (RF), microwave (MW), or laser energy; temperature induced changes in tissue properties are of relevance in relation to predicting tissue temperature profile, monitoring during treatment, and evaluation of treatment results. This paper presents a literature survey of temperature dependence of electrical (electrical conductivity, resistivity, permittivity) and thermal tissue properties (thermal conductivity, specific heat, diffusivity). Data of soft tissues (liver, prostate, muscle, kidney, uterus, collagen, myocardium and spleen) for temperatures between 5 to 90°C, and dielectric properties in the frequency range between 460 kHz and 3 GHz are reported. Furthermore, perfusion changes in tumors including carcinomas, sarcomas, rhabdomyosarcoma, adenocarcinoma and ependymoblastoma in response to hyperthmic temperatures up to 46°C are presented. Where appropriate, mathematical models to describe temperature dependence of properties are presented. The presented data is valuable for mathematical models that predict tissue temperature during thermal therapies (e.g. hyperthermia or thermal ablation), as well as for applications related to prediction and monitoring of temperature induced tissue changes.