Magnetic-resonance-based electrical properties tomography: a review

In-Vivo Electrical Properties Estimation of Biological Tissues by Means of a Multi-Step Microwave Tomography Approach

The accurate quantitative estimation of the electromagnetic properties of tissues can serve important diagnostic and therapeutic medical purposes. Quantitative microwave tomography is an imaging modality that can provide maps of the in-vivo electromagnetic properties of the imaged tissues, i.e. both the permittivity and the electric conductivity. A multi-step microwave tomography approach is proposed for the accurate retrieval of such spatial maps of biological tissues. The underlying idea behind the new imaging approach is to progressively add details to the maps in a step-wise fashion starting from single-frequency qualitative reconstructions. Multi-frequency microwave data is utilized strategically in the final stage. The approach results in improved accuracy of the reconstructions compared to inversion of the data in a single step. As a case study, the proposed workflow was tested on an experimental microwave data set collected for the imaging of the human forearm. The human forearm is a good test case as it contains several soft tissues as well as bone, exhibiting a wide range of values for the electrical properties.

MR-based electrical property tomography using a physics-informed network at 3 and 7 T

Magnetic resonance electrical propert tomography promises to retrieve electrical properties (EPs) quantitatively and non-invasively in vivo, providing valuable information for tissue characterization and pathology diagnosis. However, its clinical implementation has been hindered by, for example, B1 measurement accuracy, reconstruction artifacts resulting from inaccuracies in underlying models, and stringent hardware/software requirements. To address these challenges, we present a novel approach aimed at accurate and high-resolution EPs reconstruction based on water content maps by using a physics-informed network (PIN-wEPT). The proposed method utilizes standard clinical protocols and conventional multi-channel receive arrays that have been routinely equipped in clinical settings, thus eliminating the need for specialized RF sequence/coil configurations. Compared with the original wEPT method, the network generates accurate water content maps that effectively eliminate the influence ofB → 1 + $$ {\overrightarrow{B}}_1^{+} $$ andB → 1 - $$ {\overrightarrow{B}}_1^{-} $$ by incorporating data mismatch with electrodynamic constraints derived from the Helmholtz equation. Subsequent regression analysis develops a broad relationship between water content and EPs across various types of brain tissue. A series of numerical simulations was conducted at 7 T to assess the feasibility and performance of the method, which encompassed four normal head models and models with tumorous tissues incorporated, and the results showed normalized mean square error below 1.0% in water content, below 11.7% in conductivity, and below 1.1% in permittivity reconstructions for normal brain tissues. Moreover, in vivo validations conducted over five healthy subjects at both 3 and 7 T showed reasonably good consistency with empirical EPs values across the white matter, gray matter, and cerebrospinal fluid. The PIN-wEPT method, with its demonstrated efficacy, flexibility, and compatibility with current MRI scanners, holds promising potential for future clinical application.

Correlation between standardized uptake value of 18F-FDG PET/CT and conductivity with pathologic prognostic factors in breast cancer

We investigated the correlation between standardized uptake value (SUV) of 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET)/computed tomography (CT) and conductivity parameters in breast cancer and explored the feasibility of conductivity as an imaging biomarker. Both SUV and conductivity have the potential to reflect the tumors' heterogeneous characteristics, but their correlations have not been investigated until now. Forty four women diagnosed with breast cancer who underwent breast MRI and 18F-FDG PET/CT at the time of diagnosis were included. Among them, 17 women received neoadjuvant chemotherapy followed by surgery and 27 women underwent upfront surgery. For conductivity parameters, maximum and mean values of the tumor region-of-interests were examined. For SUV parameters, SUVmax, SUVmean, and SUVpeak of the tumor region-of-interests were examined. Correlations between conductivity and SUV were evaluated, and among them, the highest correlation was observed between mean conductivity and SUVpeak (Spearman's correlation coefficient = 0.381). In a subgroup analysis for 27 women with upfront surgery, tumors with lymphovascular invasion (LVI) showed higher mean conductivity than those without LVI (median: 0.49 S/m vs 0.06 S/m, p < 0.001). In conclusion, our study shows a low positive correlation between SUVpeak and mean conductivity in breast cancer. Furthermore, conductivity showed a potential to noninvasively predict LVI status.

Evaluating Compressed SENSE (CS) MRI Metal Artifact Reduction Using Pig L-Spine Phantom and Transplant Patients: Focused on the CS-SEMAC (SPIR), mDixon(O-MAR) and STIR Techniques

This study evaluates the clinical usefulness of the images obtained after applying mDixon (O-MAR), CS-SEMAC (SPIR), and STIR techniques to Pig L-Spine Phantom and transplant patients according to the difference in the reduction in metal artifacts and provides the optimal MAR image technique. This study was conducted with Phantom and 30 transplant patients who had an implant on the L-Spine (22 men, 8 women, mean age: 64.2 ± 12.98). All data analyzed were evaluated, using Philips Ingenia 3.0T CX. As pulse sequences, applied to the analysis, mDixon (O-MAR), CS-SEMAC (SPIR), and STIR were used. As the coil used to obtain data, the dStream Head Spine Coil was used. When tested directly applying to the transplant patients in the conditions the same as for the Phantom, as for the MAR effect of T1 and T2 images, the SNR value showed the highest effect on the increase in the signal in T1, T2 CS-SEMAC (SPIR), followed by mDixon (O-MAR) and STIR, which was the same result as the Phantom (p < 0.05). In addition, in the results of the histogram measurement in both of the subjects, Phantom and transplant patients, the count of T1, the T2 Sagittal image was the highest in T1, T2 STIR, followed by T1, T2 mDixon (O-MAR) and T1, and T2 CS-SEMAC (SPIR). As a result of the qualitative analysis, the quality was the best in T2 CS-SEMAC(SPIR) (c), followed by mDixon (O-MAR) (b) and T2 STIR (a). In conclusion, when the MAR effect on the Pig L-spine Phantom and Transplant patients was compared, it was noted that the CS-SEMAC (SPIR) technique was the most excellent in the following order: STIR < mDixon (O-MAR) < CS-SEMAC (SPIR).

Correlation analysis between the complex electrical permittivity and relaxation time of tissue mimicking phantoms in 7 T MRI

Dielectric relaxation theory describes the complex permittivity of a material in an alternating field; in particular, Debye theory relates the time it takes for an applied field to achieve the maximum polarization and the electrical properties of the material. Although, Debye’s equations were proposed for electrical polarization, in this study, we investigate the correlation between the magnetic longitudinal relaxation time T1 and the complex electrical permittivity of tissue-mimicking phantoms using a 7 T magnetic resonance scanner. We created phantoms that mimicked several human tissues with specific electrical properties. The electrical properties of the phantoms were measured using bench-test equipment. T1 values were acquired from phantoms using MRI. The measured values were fitted with functions based on dielectric estimations, using relaxation times of electrical polarization, and the mixture theory for dielectrics. The results show that, T1 and the real permittivity are correlated; therefore, the correlation can be approximated with a rational function in the case of water-based phantoms. The correlation between index loss and T1 was determined using a fitting function based on the Debye equation and mixture theory equation, in which the fraction of the materials was taken into account. This phantom study and analysis provide an insight into the application relaxation times used for estimating dielectric properties. Currently, the measurement of electrical properties based on dielectric relaxation theory is based on an antenna, sometimes invasive, that irradiates an electric field into a small sample; thus, it is not possible to create a map of electrical properties for a complex structure such as the human body. This study could be further used to compute the electrical properties maps of tissues by scanning images and measuring T1 maps.

Measurement and image-based estimation of dielectric properties of biological tissues —past, present, and future—

The dielectric properties of biological tissues are fundamental pararmeters that are essential for electromagnetic modeling of the human body. The primary database of dielectric properties compiled in 1996 on the basis of dielectric measurements at frequencies from 10 Hz to 20 GHz has attracted considerable attention in the research field of human protection from non-ionizing radiation. This review summarizes findings on the dielectric properties of biological tissues at frequencies up to 1 THz since the database was developed. Although the 1996 database covered general (normal) tissues, this review also covers malignant tissues that are of interest in the research field of medical applications. An intercomparison of dielectric properties based on reported data is presented for several tissue types. Dielectric properties derived from image-based estimation techniques developed as a result of recent advances in dielectric measurement are also included. Finally, research essential for future advances in human body modeling is discussed.

On the measurement of skeletal muscle anisotropic permittivity property with a single cross-shaped needle insertion

Application of minimally invasive methods to enable the measurement of tissue permittivity in the neuromuscular clinic remain elusive. This paper provides a theoretical and modeling study on the measurement of the permittivity of two-dimensional anisotropic tissues such as skeletal muscle with a multi-electrode cross-shaped needle. For this, we design a novel cross-shaped needle with multiple-electrodes and analyse apparent impedance corresponding to the measured impedance. In addition, we propose three methods of estimate anisotropic muscle permittivity. Compared to existing electrical impedance-based needle methods that we have developed, the new needle design and numerical methods associated enable estimating in vivo muscle permittivity values with only a single needle insertion. Being able to measure muscle permittivity directly with a single needle insertion could open up an entirely new area of research with direct clinical application, including using these values to assist in neuromuscular diagnosis and to assess subtle effects of therapeutic intervention on muscle health.

Combining deep learning and 3D contrast source inversion in MR-based electrical properties tomography

Magnetic resonance electrical properties tomography (MR-EPT) is a technique used to estimate the conductivity and permittivity of tissues from MR measurements of the transmit magnetic field. Different reconstruction methods are available; however, all these methods present several limitations, which hamper the clinical applicability. Standard Helmholtz-based MR-EPT methods are severely affected by noise. Iterative reconstruction methods such as contrast source inversion electrical properties tomography (CSI-EPT) are typically time-consuming and are dependent on their initialization. Deep learning (DL) based methods require a large amount of training data before sufficient generalization can be achieved. Here, we investigate the benefits achievable using a hybrid approach, that is, using MR-EPT or DL-EPT as initialization guesses for standard 3D CSI-EPT. Using realistic electromagnetic simulations at 3 and 7 T, the accuracy and precision of hybrid CSI reconstructions are compared with those of standard 3D CSI-EPT reconstructions. Our results indicate that a hybrid method consisting of an initial DL-EPT reconstruction followed by a 3D CSI-EPT reconstruction would be beneficial. DL-EPT combined with standard 3D CSI-EPT exploits the power of data-driven DL-based EPT reconstructions, while the subsequent CSI-EPT facilitates a better generalization by providing data consistency.

In vivo electrical conductivity measurement of muscle, cartilage, and peripheral nerve around knee joint using MR-electrical properties tomography

This study aimed to investigate whether in vivo MR-electrical properties tomography (MR-EPT) is feasible in musculoskeletal tissues by evaluating the conductivity of muscle, cartilage, and peripheral nerve around the knee joint, and to explore whether these measurements change after exercise. This prospective study was approved by the institutional review board. On February 2020, ten healthy volunteers provided written informed consent and underwent MRI of the right knee using a three-dimensional balanced steady-state free precession (bSSFP) sequence. To test the effect of loading, the subjects performed 60 squatting exercises after baseline MRI, immediately followed by post-exercise MRI with the same sequences. After reconstruction of conductivity map based on the bSSFP sequence, conductivity of muscles, cartilages, and nerves were measured. Measurements between the baseline and post-exercise MRI were compared using the paired t-test. Test–retest reliability for baseline conductivity was evaluated using the intraclass correlation coefficient. The baseline and post-exercise conductivity values (mean ± standard deviation) [S/m] of muscles, cartilages, and nerves were 1.73 ± 0.40 and 1.82 ± 0.50 (p = 0.048), 2.29 ± 0.47 and 2.51 ± 0.37 (p = 0.006), and 2.35 ± 0.57 and 2.36 ± 0.57 (p = 0.927), respectively. Intraclass correlation coefficient for the baseline conductivity of muscles, cartilages, and nerves were 0.89, 0.67, and 0.89, respectively. In conclusion, in vivo conductivity measurement of musculoskeletal tissues is feasible using MR-EPT. Conductivity of muscles and cartilages significantly changed with an overall increase after exercise.

Visible Human Project® female surface based computational phantom (Nelly) for radio-frequency safety evaluation in MRI coils

Quantitative modeling of specific absorption rate and temperature rise within the human body during 1.5 T and 3 T MRI scans is of clinical significance to ensure patient safety. This work presents justification, via validation and comparison, of the potential use of the Visible Human Project (VHP) derived Computer Aided Design (CAD) female full body computational human model for non-clinical assessment of female patients of age 50–65 years with a BMI of 30–36 during 1.5 T and 3 T based MRI procedures. The initial segmentation validation and four different application examples have been identified and used to compare to numerical simulation results obtained using VHP Female computational human model under the same or similar conditions. The first application example provides a simulation-to-simulation validation while the latter three application examples compare with measured experimental data. Given the same or similar coil settings, the computational human model generates meaningful results for SAR, B1 field, and temperature rise when used in conjunction with the 1.5 T birdcage MRI coils or at higher frequencies corresponding to 3 T MRI. Notably, the deviation in temperature rise from experiment did not exceed 2.75° C for three different heating scenarios considered in the study with relative deviations of 10%, 25%, and 20%. This study provides a reasonably systematic validation and comparison of the VHP-Female CAD v.3.0–5.0 surface-based computational human model starting with the segmentation validation and following four different application examples.

High-frequency electrical properties tomography at 9.4T as a novel contrast mechanism for brain tumors

PURPOSE

To establish high-frequency magnetic resonance electrical properties tomography (MREPT) as a novel contrast mechanism for the assessment of glioblastomas using a rat brain tumor model.

METHODS

Six F98 intracranial tumor bearing rats were imaged longitudinally 8, 11 and 14 days after tumor cell inoculation. Conductivity and mean diffusivity maps were generated using MREPT and Diffusion Tensor Imaging. These maps were co-registered with T2 -weighted images and volumes of interests (VOIs) were segmented from the normal brain, ventricles, edema, viable tumor, tumor rim, and tumor core regions. Longitudinal changes in conductivity and mean diffusivity (MD) values were compared in these regions. A correlation analysis was also performed between conductivity and mean diffusivity values.

RESULTS

The conductivity of ventricles, edematous area and tumor regions (tumor rim, viable tumor, tumor core) was significantly higher (P < .01) compared to the contralateral cortex. The conductivity of the tumor increased over time while MD from the tumor did not change. A marginal positive correlation was noted between conductivity and MD values for tumor rim and viable tumor, whereas this correlation was negative for the tumor core.

CONCLUSION

We demonstrate a novel contrast mechanism based on ionic concentration and mobility, which may aid in providing complementary information to water diffusion in probing the microenvironment of brain tumors.

Metamaterial-Inspired Radiofrequency (RF) Shield With Reduced Specific Absorption Rate (SAR) and Improved Transmit Efficiency for UHF MRI

To prevent the interferences between radiofrequency (RF) coils and other components in the magnetic resonance imaging (MRI) system such as gradient coils, it is essential to place an RF shield between the RF coils and gradient coils. However, the induced currents on conventional RF shields have negative influences on the RF coil performance. To reduce these influences, metamaterial absorbers (MA), a class of metamaterials exhibiting nearly unity absorption rate for the incident electromagnetic fields, can be employed for the design of a novel RF shield. However, the adoption of metamaterials in MRI systems is usually problematic because of the bulkiness of the metamaterial structure. In this work, capacitors and metallic interconnectors are used to miniaturize the MA so that the unit MA cell can operate at the Larmor frequencies of 7T and 9.4T MRI and stay compact. This MA-RF shield is used to improve the transmit efficiency of RF surface coils and reduce the specific absorption rate (SAR) in the region of interest (ROI). It is successfully demonstrated by simulations and experiments that, compared with conventional RF shield structure, the transmit efficiency can be enhanced by more than 32% and the peak SAR value can be reduced by 22% using the MA-RF shield. Moreover, it is observed that the transmit field penetration is improved when the surface coil is used with the MA-RF shield. This proof-of-concept study suggests a new practical way for the utilization of metamaterials in ultra-high field MRI applications.

Radio-Frequency Vector Magnetic Field Mapping in Magnetic Resonance Imaging

A method is presented to measure the radio-frequency (RF) vector magnetic field inside an object using magnetic resonance imaging (MRI). Conventional " [Formula: see text] mapping" in MRI can measure the proton co-rotating component ( [Formula: see text] of the RF field produced by a transmit coil. Here we show that by repeating [Formula: see text] mapping on the same object and coil at multiple (8) specific orientations with respect to the main magnet, the magnitudes and relative phases of all (x, y, z) Cartesian components of the RF field can be determined unambiguously. We demonstrate the method on a circularly polarized volume coil and a loop coil tuned at 123.25 MHz in a 3 Tesla MRI scanner, with liquid phantoms. The volume coil measurement showed the axial component of the RF field, which is normally unmeasurable in MRI, away from the center of the coil. The measured RF vector field maps of both coils compared favorably with numerical simulation, with volumetric normalized root-mean-square difference in the range of 7~20%. While the proposed method cannot be applied to human imaging at present, applications to phantoms and small animals could provide a useful experimental tool to validate RF simulation and verify certain assumptions in [Formula: see text] map-based electrical properties tomography (EPT).

A Survey on State-of-the-art Denoising Techniques for Brain Magnetic Resonance Images

The accuracy of the magnetic resonance (MR) image diagnosis depends on the quality of the image, which degrades mainly due to noise and artifacts. The noise is introduced because of erroneous imaging environment or distortion in the transmission system. Therefore, denoising methods play an important role in enhancing the image quality. However, a tradeoff between denoising and preserving the structural details is required. Most of the existing surveys are conducted on a specific MR image modality or on limited denoising schemes. In this context, a comprehensive review on different MR image denoising techniques is inevitable. This survey suggests a new direction in categorizing the MR image denoising techniques. The categorization of the different image models used in medical image processing serves as the basis of our classification. This study includes recent improvements on deep learning-based denoising methods alongwith important traditional MR image denoising methods. The major challenges and their scope of improvement are also discussed. Further, many more evaluation indices are considered for a fair comparison. An elaborate discussion on selecting appropriate method and evaluation metric as per the kind of data is presented. This study may encourage researchers for further work in this domain.

A comparative analysis framework of 3T and 7T TOF-MRA based on automated cerebrovascular segmentation

PURPOSE

High field strength 3T and 7T Time-Of-Flight Magnetic Resonance Angiography (TOF- MRA) achieves better visualization of intracranial vessels, so it attracts much attention. However, quantitative comparison between 3T and 7T MRA is lacking in the aspects of image quality and the practical application of cerebrovascular diseases.

METHODS

In this paper, a quantitative framework of 3T and 7T TOF-MRA comparison is proposed, which contains two steps including the automated cerebrovascular segmentation and statistical analysis. Firstly, the whole vascular structures on both 3T and 7T TOF-MRA images are segmented automatically, especially those small blood vessels in 7T MRA. The skeleton extraction-based automatic seed point detection is implemented to ensure the segmented vascular structure complete and precise. Secondly, the statistical analysis of the differences between 3T and 7T MRA is carried out in the aspects of image quality and the characteristics of some important vessels. The objects of statistical analysis are achieved and analyzed automatically without needing the time- consuming human beings' participation, therefore, it is efficient and objective.

RESULTS

The comparison experiments on seven pairs of 3T and 7T TOF MRA images validated that about image quality, the contrast-to-noise ratio of 7T MRA was about 4.53 ± 0.95 times as much as that of 3T MRA. About the cerebrovascular information, small vessels were more abundant in 7T MRA compared with 3T MRA (branches number: 462.0 ± 58.5 vs 393.1 ± 63.3).

CONCLUSIONS

The proposed framework can segment the whole cerebrovascular structure automatically and compare TOF-MRA with different field strengths objectively and quantitatively. It is helpful for clinical cerebrovascular disease, especially cerebral small vessel diseases.

Electrical Properties Tomography: A Methodological Review

Electrical properties tomography (EPT) is an imaging method that uses a magnetic resonance (MR) system to non-invasively determine the spatial distribution of the conductivity and permittivity of the imaged object. This manuscript starts by providing clear definitions about the data required for, and acquired in, EPT, followed by comprehensively formulating the physical equations underlying a large number of analytical EPT techniques. This thorough mathematical overview of EPT harmonizes several EPT techniques in a single type of formulation and gives insight into how they act on the data and what their data requirements are. Furthermore, the review describes machine learning-based algorithms. Matlab code of several differential and iterative integral methods is available upon request.

Configuration-based electrical properties tomography

PURPOSE

To introduce phase-based conductivity mapping from a configuration space analysis.

METHODS

The frequency response function of balanced SSFP (bSSFP) is used to perform a configuration space analysis. It is shown that the transceive phase for conductivity mapping can be directly obtained by a simple fast Fourier transform of a series of phase-cycled bSSFP scans. For validation, transceive phase and off-resonance mapping with fast Fourier transform is compared with phase estimation using a recently proposed method, termed PLANET. Experiments were performed in phantoms and for in vivo brain imaging at 3 T using a quadrature head coil.

RESULTS

For fast Fourier transform, aliasing can lead to systematic phase errors. This bias, however, decreases rapidly with increasing sampling points. Interestingly, Monte Carlo simulations revealed a lower uncertainty for the transceive phase and the off-resonance using fast Fourier transform as compared with PLANET. Both methods, however, essentially retrieve the same phase information from a set of phase-cycled bSSFP scans. As a result, configuration-based conductivity mapping was successfully performed using eight phase-cycled bSSFP scans in the phantoms and for brain tissues. Overall, the retrieved values were in good agreement with expectations. Conductivity estimation and mapping of the field inhomogeneities can therefore be performed in conjunction with the estimation of other quantitative parameters, such as relaxation, using configuration theory.

CONCLUSIONS

Phase-based conductivity mapping can be estimated directly from a simple Fourier analysis, such as in conjunction with relaxometry, using a series of phase-cycled bSSFP scans.

Magnetic resonance electrical property mapping at 21.1 T: a study of conductivity and permittivity in phantoms, ex vivo tissue and in vivo ischemia

Electrical properties (EP), namely conductivity and permittivity, can provide endogenous contrast for tissue characterization. Using electrical property tomography (EPT), maps of EP can be generated from conventional MRI data. This report investigates the feasibility and accuracy of EPT at 21.1 T for multiple RF coils and modes of operation using phantoms. Additionally, it demonstrates the EP of the in vivo rat brain with and without ischemia. Helmholtz-based EPT was implemented in its Full-form, which demands the complex [Formula: see text] field, and a simplified form requiring either just the [Formula: see text] field phase for conductivity or the [Formula: see text] field magnitude for permittivity. Experiments were conducted at 21.1 T using birdcage and saddle coils operated in linear or quadrature transceive mode, respectively. EPT approaches were evaluated using a phantom, ex and in vivo Sprague-Dawley rats under naïve conditions and ischemic stroke via transient middle cerebral artery occlusion. Different conductivity reconstruction approaches applied to the phantom displayed average errors of 12%-73% to the target acquired from dielectric probe measurements. Permittivity reconstructions showed higher agreement and an average 3%-8% error to the target depending on reconstruction approach. Conductivity and permittivity of ex and in vivo rodent brain were measured. Elevated EP in the ischemia region correlated with the increased sodium content and the influx of water intracellularly following ischemia in the lesion were detected. The Full-form technique generated from the linear birdcage provided the best accuracy for EP of the phantom. Phase-based conductivity and magnitude-based permittivity mapping provided reasonable estimates but also demonstrated the limitations of Helmholtz-based EPT at 21.1 T. Permittivity reconstruction was improved significantly over lower fields, suggesting a novel metric for in vivo brain studies. EPT applied to ischemic rat brain proved sensitivity to physiological changes, motivating the future application of more advanced reconstruction approaches.

Role of Simulations in the Treatment Planning of Radiofrequency Hyperthermia Therapy in Clinics

Hyperthermia therapy is a treatment modality in which tumor temperatures are elevated to higher temperatures to cause damage to cancerous tissues. Numerical simulations are integral in the development of hyperthermia treatment systems and in clinical treatment planning. In this study, simulations in radiofrequency hyperthermia therapy are reviewed in terms of their technical development and clinical aspects for effective clinical use. This review offers an overview of mathematical models and the importance of tissue properties; locoregional mild hyperthermia therapy, including phantom and realistic human anatomy models; phase array systems; tissue damage; thermal dose analysis; and thermoradiotherapy planning. This review details the improvements in numerical approaches in treatment planning and their application for effective clinical use. Furthermore, the modeling of thermoradiotherapy planning, which can be integrated with radiotherapy to provide combined hyperthermia and radiotherapy treatment planning strategies, are also discussed. This review may contribute to the effective development of thermoradiotherapy planning in clinics.

Parallel Transmission for Ultrahigh Field MRI

Magnetic resonance imaging (MRI) has been driven toward ultrahigh magnetic fields (UHF) in order to benefit from correspondingly higher signal-to-noise ratio and spectral resolution. Technological challenges associated with UHF, such as increased radiofrequency (RF) energy deposition and RF excitation inhomogeneity, limit realization of the full potential of these benefits. Parallel RF transmission (pTx) enables decreases in the inhomogeneity of RF excitations and in RF energy deposition by using multiple-transmit RF coils driven independently and operating simultaneously. pTx plays a fundamental role in UHF MRI by bringing the potential applications of UHF into reality. In this review article, we review the recent developments in pTx pulse design and RF safety in pTx. Simultaneous multislice imaging and inner volume imaging using pTx are reviewed with a focus on UHF applications. Emerging pTx design approaches using improved pTx design frameworks and calibrations are reviewed together with calibration-free approaches that remove the necessity of time-consuming calibrations necessary for successful pTx. Lastly, we focus on the safety of pTx that is improved by using intersubject variability analysis, proactively managing pTx and temperature-based pTx approaches.

IMPULSE: A scalable algorithm for design of minimum specific absorption rate parallel transmit RF pulses

PURPOSE

Managing local specific absorption rate (SAR) in parallel transmission requires ensuring that the peak SAR over a large number of voxels (> 10 5 ) is below the regulatory limit. The safety risk to the patient depends on cumulative (not instantaneous) SAR thus making a joint design of all RF pulses in a sequence desirable. We propose the Iterative Minimization Procedure with Uncompressed Local SAR Estimate (IMPULSE), an efficient optimization formulation and algorithm that can handle uncompressed SAR matrices and optimize pulses for all slices jointly within a practical time frame.

THEORY AND METHODS

IMPULSE optimizes parallel transmit pulses for small-tip-angle slice selective excitation to minimize a single cost function incorporating multiple quantities (local SAR, global SAR, and per-channel power) averaged over the entire multislice scan subject to a strict constraint on excitation accuracy. Pulses for an 8-channel 7T head coil were designed with IMPULSE and compared with pulses designed using generic optimization algorithms and VOPs to assess the computation time and SAR performance benefits.

RESULTS

IMPULSE achieves lower SAR and shorter computation time compared with a VOP approach. Compared with the generic sequential quadratic programming algorithm, computation time is reduced by a factor of 5-6 by using IMPULSE. Using as many as 6 million local SAR terms, up to 120 slices can be designed jointly with IMPULSE within 45 s.

CONCLUSIONS

IMPULSE can handle significantly larger number of SAR matrices and slices than conventional optimization algorithms, enabling the use of uncompressed or partially compressed SAR matrices to design pulses for a multislice scan in a practical time frame.

Quantitative Estimation of the Equivalent Radiation Dose Escalation using Radiofrequency Hyperthermia in Mouse Xenograft Models of Human Lung Cancer

Hyperthermia is a potent radiosensitizer, and its effect varies according to the different types of cancer cells. In the present study, the radiosensitizing effect of hyperthermia on lung cancer cell lines A549 and NCI-H1299 was determined based on the equivalent radiation dose escalation. In vitro cell experiments were conducted using lung cancer cell lines A549 and NCI-H1299 to determine thermal radiosensitivity. In vivo experiments were conducted using mouse heterotopic xenograft models to determine the treatment response and increase in the temperature of tumors using a 13.56 MHz radiofrequency (RF) hyperthermia device. Using the α and β values of the linear–quadratic equations of cell survival curves, numerical simulations were performed to calculate the equivalent radiation dose escalations. The dielectric properties of tumors were measured, and their effect on the calculated equivalent radiation dose was analyzed. Hyperthermia increased the equivalent radiation dose of lung cancer xenografts and a higher escalation was found in NCI-H1299 cells compared with that observed in A549 cells. An underestimation of the calculated equivalent radiation dose was observed when the dielectric property of the tumor was varied. This study may contribute to the effective planning of thermoradiotherapy in clinics.

Developments in Electrical-Property Tomography Based on the Contrast-Source Inversion Method

The main objective of electrical-property tomography (EPT) is to retrieve dielectric tissue parameters from B ^ 1 + data as measured by a magnetic-resonance (MR) scanner. This is a so-called hybrid inverse problem in which data are defined inside the reconstruction domain of interest. In this paper, we discuss recent and new developments in EPT based on the contrast-source inversion (CSI) method. After a short review of the basics of this method, two- and three-dimensional implementations of CSI-EPT are presented along with a very efficient variant of 2D CSI-EPT called first-order induced current EPT (foIC-EPT). Practical implementation issues that arise when applying the method to measured data are addressed as well, and the limitations of a two-dimensional approach are extensively discussed. Tissue-parameter reconstructions of an anatomically correct male head model illustrate the performance of two- and three-dimensional CSI-EPT. We show that 2D implementation only produces reliable reconstructions under very special circumstances, while accurate reconstructions can be obtained with 3D CSI-EPT.

An MR technique for simultaneous quantitative imaging of water content, conductivity and susceptibility, with application to brain tumours using a 3T hybrid MR-PET scanner

Approaches for the quantitative mapping of water content, electrical conductivity and susceptibility have been developed independently. The purpose of this study is to develop a method for simultaneously acquiring quantitative water content, electrical conductivity and susceptibility maps based on a 2D multi-echo gradient echo sequence. Another purpose is to investigate the changes in these properties caused by brain tumours. This was done using a 3T hybrid magnetic resonance imaging and positron emission tomography (MR-PET) scanner. Water content maps were derived after performing T₂* and transmit-receive field bias corrections to magnitude images essentially reflecting only the H₂O content contrast. Phase evolution during the multi-echo train was used to generate field maps and derive quantitative susceptibility, while the conductivity maps were retrieved from the phase value at zero echo time. Performance of the method is demonstrated on phantoms and two healthy volunteers. In addition, the method was applied to three patients with brain tumours and a comparison to maps obtained from PET using O-(2-[18 F]fluoroethyl)-L-tyrosine and clinical MR images is presented. The combined information of the water content, conductivity and susceptibility may provide additional information about the tissue viability. Future studies can benefit from the evaluation of these contrasts with shortened acquisition times.

Mapping electrical properties heterogeneity of tumor using boundary informed electrical properties tomography (BIEPT) at 7T

PURPOSES

To develop and evaluate a boundary informed electrical properties tomography (BIEPT) technique for high-resolution imaging of tumor electrical properties (EPs) heterogeneity on a rodent tumor xenograft model.

METHODS

Tumor EP distributions were inferred from a reference area external to the tumor, as well as internal EP spatial variations derived from a plurality of relative transmit B1 measurements at 7T. Edge sparsity constraint was enforced to enhance numerical stability. Phantom experiments were performed to determine the imaging accuracy and sensitivity for structures of various EP values, as well as geometrical sizes down to 1.5 mm. Numerical simulation of a realistic rodent model was used to quantify the algorithm performance in the presence of noise. Eleven athymic rats with human breast cancer xenograft were imaged in vivo, and representative pathological samples were acquired for comparison.

RESULTS

Reconstructed EPs of the phantoms correspond well to the ground truth acquired from dielectric probe measurements, with the smallest structure reliably detectable being 3 mm. EPs heterogeneity inside a tumor is successfully retrieved in both simulated and experimental cases. In vivo tumor imaging results demonstrate similar local features and spatial patterns to anatomical MRI and pathological slides. The imaged conductivity of necrotic tissue is higher than that of viable tissues, which agrees with our expectation.

CONCLUSION

BIEPT enables robust detection of tumor EPs heterogeneity with high accuracy and sensitivity to small structures. The retrieved quantitative EPs reflect tumor pathological features (e.g., necrosis). These results provide strong rationale to further expand BIEPT studies toward pathological conditions where EPs may yield valuable, non-invasive biomarkers.

Electrical properties tomography: Available contrast and reconstruction capabilities

MR-based electrical properties tomography converts the MRI transmit/receive RF field measurements to tissue electrical property maps through dedicated reconstruction algorithms. Recent reports showed that despite limitations, electrical properties tomography holds promise for generating additional contrast for tumor detection and patient-specific modeling of tissue-RF field interactions. This review summarizes the available tissue electrical property contrasts and compares them with the capabilities of the most commonly used electrical properties tomography reconstruction method. Future directions and prospects of clinical translation are discussed.

An intraoperative motor tract positioning method in brain tumor surgery: technical note

Intraoperative 3D recognition of the motor tract is indispensable to avoiding neural fiber injury in brain tumor surgery. However, precise localization of the tracts is sometimes difficult with conventional mapping methods. Thus, the authors developed a novel brain mapping method that enables the 3D recognition of the motor tract for intrinsic brain tumor surgeries. This technique was performed in 40 consecutive patients with gliomas adjacent to motor tracts that have a risk of intraoperative pyramidal tract damage. Motor tracts were electrically stimulated and identified by a handheld brain-mapping probe, the NY Tract Finder (NYTF). Sixteen-gauge plastic tubes were mounted onto the NYTF and inserted in the estimated direction of the motor tract with reference to navigational information. Only the NYTF was removed, leaving the plastic tubes in their places, immediately after muscle motor evoked potentials were recorded at the minimum stimulation current. Motor tracts were electrically identified in all cases. Three-dimensional information on the position of motor tracts was given by plastic tubes that were neurophysiologically placed. Tips of tubes showed the resection limit during tumor removal. Safe tumor resection with an arbitrary safety margin can be performed by adjusting the length of the plastic tubes. The motor tract positioning method enabled the 3D recognition of the motor tract by surgeons and provided for safe resection of tumors. Tumor resections were performed safely before damaging motor tracts, without any postoperative neurological deterioration.

MR-based electrical property tomography using a modified finite difference scheme

Magnetic resonance electrical property tomography (MR-EPT) reconstructs electrical properties (EPs) from measured magnetic fields in magnetic resonance imaging (MRI) systems. In this study, an MR-EPT method was proposed that utilized a new finite difference approximation of the involved differential wave equation. Compared with existing MR-EPT approaches, the construction of the system matrix involves applying the first derivative twice based on a larger number of neighbouring finite-difference grids, which is different from a standard Laplacian operator on a regular grid structure, leading to a better conditioned linear inverse problem. With improved noise robustness, more faithful EPs can be obtained by the proposed method, particularly at tissue boundaries and regions with a poorly measured magnetic field (low signal-to-noise ratio). Numerical simulations with a specially designed multi-slice phantom and an anatomically accurate head model (Duke) have demonstrated that the proposed method can provide a more faithful reconstruction of EPs compared to existing methods, which usually offer unreliable solutions associated with traditional finite difference approximation of the central wave equation and unrealistic assumptions. Experiments on a 9.4 T MRI system have been conducted to validate the simulations.

Computational method for estimating boundary of abdominal subcutaneous fat for absolute electrical impedance tomography

Abdominal fat accumulation is considered an essential indicator of human health. Electrical impedance tomography has considerable potential for abdominal fat imaging because of the low specific conductivity of human body fat. In this paper, we propose a robust reconstruction method for high-fidelity conductivity imaging by abstraction of the abdominal cross section using a relatively small number of parameters. Toward this end, we assume homogeneous conductivity in the abdominal subcutaneous fat area and characterize its geometrical shape by parameters defined as the ratio of the distance from the center to boundary of subcutaneous fat to the distance from the center to outer boundary in 64 equiangular directions. To estimate the shape parameters, the sensitivity of the noninvasively measured voltages with respect to the shape parameters is formulated for numerical optimization. Numerical simulations are conducted to demonstrate the validity of the proposed method. A 3-dimensional finite element method is used to construct a computer model of the human abdomen. The inverse problems of shape parameters and conductivities are solved concurrently by iterative forward and inverse calculations. As a result, conductivity images are reconstructed with a small systemic error of less than 1% for the estimation of the subcutaneous fat area. A novel method is devised for estimating the boundary of the abdominal subcutaneous fat. The fidelity of the overall reconstructed image to the reference image is significantly improved. The results demonstrate the possibility of realization of an abdominal fat scanner as a low-cost, radiation-free medical device.

In vivo imaging of electrical properties of an animal tumor model with an 8-channel transceiver array at 7 T using electrical properties tomography

PURPOSE

To develop and evaluate a technique for imaging electrical properties ((EPs), conductivity and permittivity) of an animal tumor model in vivo using MRI.

METHODS

Electrical properties were reconstructed from the calculated EP gradient, which was derived using two sets of measured transmit B1 magnitude and relative phase maps with the sample and radiofrequency (RF) coil oriented in the positive and negative z-directions, respectively. An eight-channel transceiver microstrip array RF coil fitting the size of the animal was developed for generating and mapping B1 fields to reconstruct EPs. The technique was evaluated at 7 tesla using a physical phantom and in vivo on two Copenhagen rats with subcutaneously implanted AT-1 rat prostate cancer on a hind limb.

RESULTS

The reconstructed EPs in the phantom experiment was in good agreement with the target EP map determined by a dielectric probe. Reconstructed conductivity map of the animals revealed the boundary between tumor and healthy tissue consistent with the boundary indicated by T1 -weighted MRI.

CONCLUSION

A technique for imaging EP of an animal tumor model using MRI has been developed with high sensitivity, accuracy, and resolution, as demonstrated in the phantom experiment. Further animal experiments are needed to demonstrate its translational value for tumor diagnosis. Magn Reson Med 78:2157-2169, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Error analysis of helmholtz-based MR-electrical properties tomography

PURPOSE

MR electrical properties tomography (MR-EPT) aims to measure tissue electrical properties by computing spatial derivatives of measured B1+ data. This computation is very sensitive to spatial fluctuations caused, for example, by noise and Gibbs ringing. In this work, the error arising from the computation of spatial derivatives using finite difference kernels (FD error) has been investigated. In relation to this FD error, it has also been investigated whether mitigation strategies such as Gibbs ringing correction and Gaussian apodization can be beneficial for conductivity reconstructions.

METHODS

Conductivity reconstructions were performed on a phantom (by means of simulations and MR measurements at 3T) and on a human brain model. The accuracy was evaluated as a function of image resolution, FD kernel size, k-space windowing, and signal-to-noise ratio. The impact of mitigation strategies was also investigated.

RESULTS

The adopted small FD kernel is highly sensitive to spatial fluctuations, whereas the large FD kernel is more noise-robust. However, large FD kernels lead to extended numerical boundary error propagation, which severely hampers the MR-EPT reconstruction accuracy for highly spatially convoluted tissue structures such as the human brain. Mitigation strategies slightly improve the accuracy of conductivity reconstructions. For the adopted derivative kernels and the investigated scenario, MR-EPT conductivity reconstructions show low accuracy: less than 37% of the voxels have a relative error lower than 30%.

CONCLUSION

The numerical error introduced by the computation of spatial derivatives using FD kernels is one of the major causes of limited accuracy in Helmholtz-based MR-EPT reconstructions. Magn Reson Med 80:90-100, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Analytic Modeling of Neural Tissue: I. A Spherical Bidomain

Presented here is a model of neural tissue in a conductive medium stimulated by externally injected currents. The tissue is described as a conductively isotropic bidomain, i.e. comprised of intra and extracellular regions that occupy the same space, as well as the membrane that divides them, and the injection currents are described as a pair of source and sink points. The problem is solved in three spatial dimensions and defined in spherical coordinates [Formula: see text]. The system of coupled partial differential equations is solved by recasting the problem to be in terms of the membrane and a monodomain, interpreted as a weighted average of the intra and extracellular domains. The membrane and monodomain are defined by the scalar Helmholtz and Laplace equations, respectively, which are both separable in spherical coordinates. Product solutions are thus assumed and given through certain transcendental functions. From these electrical potentials, analytic expressions for current density are derived and from those fields the magnetic flux density is calculated. Numerical examples are considered wherein the interstitial conductivity is varied, as well as the limiting case of the problem simplifying to two dimensions due to azimuthal independence. Finally, future modeling work is discussed.

In Vivo Electrical Conductivity Contrast Imaging in a Mouse Model of Cancer Using High-Frequency Magnetoacoustic Tomography With Magnetic Induction (hfMAT-MI)

Cancerous tissues have electrical-conductivity signatures different from normal tissues, which contain potentially useful information for early detection. Despite recent advancements in electrical-conductivity imaging and its applications, imaging electrical conductivities with high spatial resolution remains a challenge for non-invasive diagnosis of early-stage cancer. Among the various electrical-conductivity imaging methods, magnetoacoustic tomography with magnetic induction (MAT-MI) is a promising technology for non-invasive detection of breast cancer. However, previous efforts to use MAT-MI for cancer imaging have suffered due to insufficient spatial resolution. In this work, we have developed a high-frequency MAT-MI (hfMAT-MI) system with a 2-D spatial resolution of 1 mm, a significant improvement over previous methods. Furthermore, we demonstrated the performance of this method using an in vivo cancer model in nude mice with human breast xenograft hindlimb tumors. hfMAT-MI was able to resolve not only the boundaries between cancerous and healthy tissues, but also the tumors' internal structures. Importantly, we were able to track a growing tumor using our hfMAT-MI method for the first time in an in vivo mouse model, demonstrating the promise of this magneto-acoustic imaging system for effective detection and diagnosis of early-stage breast cancer.

Analysis of bipolar external excitation of spherical tissue by spatially opposed current source and sink points

The recently increasing role in medical imaging that electrophysiology plays has spurned the need for its quantitative analysis at all scales-ions, cells, tissues, organs, etc.; so, here is presented a model of nerve tissue in a spherical volume excited by a point current source at one pole and a point current sink at the opposite pole. The sphere of tissue is described as an isotropic bidomain, consisting of the intra- and extra-cellular regions and the membrane that separates them, and is immersed in an infinite isotropic conductive bath. The system of coupled differential equations is solved by redefining the domains to be in terms of a monodomain and a membrane. The solution takes the form of an infinite sum of the product of certain transcendental functions. The study concludes with a numeric example in which the boundary conditions are shown to be satisfied, validating this analysis, paving the way for more sophisticated models of excitable tissue.

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.

Parallel transmission pulse design with explicit control for the specific absorption rate in the presence of radiofrequency errors

PURPOSE

A new framework for the design of parallel transmit (pTx) pulses is presented introducing constraints for local and global specific absorption rate (SAR) in the presence of errors in the radiofrequency (RF) transmit chain.

METHODS

The first step is the design of a pTx RF pulse with explicit constraints for global and local SAR. Then, the worst possible SAR associated with that pulse due to RF transmission errors ("worst-case SAR") is calculated. Finally, this information is used to re-calculate the pulse with lower SAR constraints, iterating this procedure until its worst-case SAR is within safety limits.

RESULTS

Analysis of an actual pTx RF transmit chain revealed amplitude errors as high as 8% (20%) and phase errors above 3° (15°) for spokes (spiral) pulses. Simulations show that using the proposed framework, pulses can be designed with controlled "worst-case SAR" in the presence of errors of this magnitude at minor cost of the excitation profile quality.

CONCLUSION

Our worst-case SAR-constrained pTx design strategy yields pulses with local and global SAR within the safety limits even in the presence of RF transmission errors. This strategy is a natural way to incorporate SAR safety factors in the design of pTx pulses. Magn Reson Med 75:2493-2504, 2016. © 2015 Wiley Periodicals, Inc.

Magneto acoustic tomography with short pulsed magnetic field for in-vivo imaging of magnetic iron oxide nanoparticles

Nanoparticles are widely used as contrast and therapeutic agents. As such, imaging modalities that can accurately estimate their distribution in-vivo are actively sought. We present here our method Magneto Acoustic Tomography (MAT), which uses magnetomotive force due to a short pulsed magnetic field to induce ultrasound in the magnetic nanoparticle labeled tissue and estimates an image of the distribution of the nanoparticles in-vivo with ultrasound imaging resolution. In this study, we image the distribution of superparamagnetic iron oxide nanoparticles (IONP) using MAT method. In-vivo imaging was performed on live, nude mice with IONP injected into LNCaP tumors grown subcutaneously within the hind limb of the mice. Our experimental results indicate that the MAT method is capable of imaging the distribution of IONPs in-vivo. Therefore, MAT could become an imaging modality for high resolution reconstruction of MNP distribution in the body.

FROM THE CLINICAL EDITOR

Many magnetic nanoparticles (MNPs) have been used as contrast agents in magnetic resonance imaging. In this study, the authors investigated the use of ultrasound to detect the presence of MNPs by magneto acoustic tomography. In-vivo experiments confirmed the imaging quality of this new approach, which hopefully would provide an alternative method for accurate tumor detection.