1
|
He Z, Lefebvre PM, Soullié P, Doguet M, Ambarki K, Chen B, Odille F. Phantom evaluation of electrical conductivity mapping by MRI: Comparison to vector network analyzer measurements and spatial resolution assessment. Magn Reson Med 2024; 91:2374-2390. [PMID: 38225861 DOI: 10.1002/mrm.30009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Revised: 12/21/2023] [Accepted: 12/22/2023] [Indexed: 01/17/2024]
Abstract
PURPOSE To evaluate the performance of various MR electrical properties tomography (MR-EPT) methods at 3 T in terms of absolute quantification and spatial resolution limit for electrical conductivity. METHODS Absolute quantification as well as spatial resolution performance were evaluated on homogeneous phantoms and a phantom with holes of different sizes, respectively. Ground-truth conductivities were measured with an open-ended coaxial probe connected to a vector network analyzer (VNA). Four widely used MR-EPT reconstruction methods were investigated: phase-based Helmholtz (PB), phase-based convection-reaction (PB-cr), image-based (IB), and generalized-image-based (GIB). These methods were compared using the same complex images from a 1 mm-isotropic UTE sequence. Alternative transceive phase acquisition sequences were also compared in PB and PB-cr. RESULTS In large homogeneous phantoms, all methods showed a strong correlation with ground truth conductivities (r > 0.99); however, GIB was the best in terms of accuracy, spatial uniformity, and robustness to boundary artifacts. In the resolution phantom, the normalized root-mean-squared error of all methods grew rapidly (>0.40) when the hole size was below 10 mm, with simplified methods (PB and IB), or below 5 mm, with generalized methods (PB-cr and GIB). CONCLUSION VNA measurements are essential to assess the accuracy of MR-EPT. In this study, all tested MR-EPT methods correlated strongly with the VNA measurements. The UTE sequence is recommended for MR-EPT, with the GIB method providing good accuracy for structures down to 5 mm. Structures below 5 mm may still be detected in the conductivity maps, but with significantly lower accuracy.
Collapse
Affiliation(s)
- Zhongzheng He
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
| | | | - Paul Soullié
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
| | - Martin Doguet
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- BioSerenity, Paris, France
| | | | - Bailiang Chen
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- CIC-IT 1433, INSERM, Université de Lorraine and CHRU Nancy, Nancy, France
| | - Freddy Odille
- IADI U1254, INSERM and Université de Lorraine, Nancy, France
- CIC-IT 1433, INSERM, Université de Lorraine and CHRU Nancy, Nancy, France
| |
Collapse
|
2
|
Bilgic B, Costagli M, Chan KS, Duyn J, Langkammer C, Lee J, Li X, Liu C, Marques JP, Milovic C, Robinson SD, Schweser F, Shmueli K, Spincemaille P, Straub S, van Zijl P, Wang Y. Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group. Magn Reson Med 2024; 91:1834-1862. [PMID: 38247051 PMCID: PMC10950544 DOI: 10.1002/mrm.30006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 10/31/2023] [Accepted: 12/14/2023] [Indexed: 01/23/2024]
Abstract
This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.
Collapse
Affiliation(s)
- Berkin Bilgic
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts, USA
| | - Mauro Costagli
- Department of Neuroscience, Rehabilitation, Ophthalmology, Genetics, Maternal and Child Sciences (DINOGMI), University of Genoa, Genoa, Italy
- Laboratory of Medical Physics and Magnetic Resonance, IRCCS Stella Maris, Pisa, Italy
| | - Kwok-Shing Chan
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, Massachusetts, USA
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Jeff Duyn
- Advanced MRI Section, NINDS, National Institutes of Health, Bethesda, Maryland, USA
| | | | - Jongho Lee
- Department of Electrical and Computer Engineering, Seoul National University, Seoul, Republic of Korea
| | - Xu Li
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Chunlei Liu
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California, USA
- Helen Wills Neuroscience Institute, University of California, Berkeley, California, USA
| | - José P Marques
- Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| | - Carlos Milovic
- School of Electrical Engineering (EIE), Pontificia Universidad Catolica de Valparaiso, Valparaiso, Chile
| | - Simon Daniel Robinson
- High Field MR Centre, Department of Biomedical Imaging and Image-Guided Therapy, Medical University of Vienna, Vienna, Austria
- Centre of Advanced Imaging, University of Queensland, Brisbane, Australia
| | - Ferdinand Schweser
- Buffalo Neuroimaging Analysis Center, Department of Neurology, Jacobs School of Medicine and Biomedical Sciences at the University at Buffalo, Buffalo, New York, USA
- Center for Biomedical Imaging, Clinical and Translational Science Institute at the University at Buffalo, Buffalo, New York, USA
| | - Karin Shmueli
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
| | - Pascal Spincemaille
- MRI Research Institute, Department of Radiology, Weill Cornell Medicine, New York, New York, USA
| | - Sina Straub
- Department of Radiology, Mayo Clinic, Jacksonville, Florida, USA
| | - Peter van Zijl
- Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, Maryland, USA
| | - Yi Wang
- MRI Research Institute, Departments of Radiology and Biomedical Engineering, Cornell University, New York, New York, USA
| |
Collapse
|
3
|
Özdemir S, Ilicak E, Zapp J, Schad LR, Zöllner FG. Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging. Magn Reson Med 2024; 91:1567-1575. [PMID: 38044757 DOI: 10.1002/mrm.29952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Revised: 10/09/2023] [Accepted: 11/13/2023] [Indexed: 12/05/2023]
Abstract
PURPOSE To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques. METHODS Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method. RESULTS The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m. CONCLUSION A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.
Collapse
Affiliation(s)
- Safa Özdemir
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Efe Ilicak
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Jascha Zapp
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Lothar R Schad
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| | - Frank G Zöllner
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
- Mannheim Institute for Intelligent Systems in Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, Germany
| |
Collapse
|
4
|
Zheng M, Lou F, Huang Y, Pan S, Zhang X. MR-based electrical property tomography using a physics-informed network at 3 and 7 T. NMR IN BIOMEDICINE 2024:e5137. [PMID: 38439522 DOI: 10.1002/nbm.5137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 01/29/2024] [Accepted: 02/11/2024] [Indexed: 03/06/2024]
Abstract
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.
Collapse
Affiliation(s)
- Mengxuan Zheng
- Interdisciplinary Institute of Neuroscience and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, China
- MOE Frontier Science Center for Brain Science and Brain-Machine Integration, Zhejiang University, Hangzhou, China
| | - Feiyang Lou
- MOE Frontier Science Center for Brain Science and Brain-Machine Integration, Zhejiang University, Hangzhou, China
- School of Medicine, Zhejiang University, Hangzhou, China
| | - Yiman Huang
- MOE Frontier Science Center for Brain Science and Brain-Machine Integration, Zhejiang University, Hangzhou, China
- College of Electrical Engineering, Zhejiang University, Hangzhou, China
| | - Sihong Pan
- MOE Frontier Science Center for Brain Science and Brain-Machine Integration, Zhejiang University, Hangzhou, China
- College of Electrical Engineering, Zhejiang University, Hangzhou, China
| | - Xiaotong Zhang
- Interdisciplinary Institute of Neuroscience and Technology, Key Laboratory for Biomedical Engineering of Ministry of Education, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, China
- MOE Frontier Science Center for Brain Science and Brain-Machine Integration, Zhejiang University, Hangzhou, China
- School of Medicine, Zhejiang University, Hangzhou, China
- College of Electrical Engineering, Zhejiang University, Hangzhou, China
- Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou, China
| |
Collapse
|
5
|
Jung K, Mandija S, Cui C, Kim J, Al‐masni MA, Meerbothe TG, Park M, van den Berg CAT, Kim D. Data-driven electrical conductivity brain imaging using 3 T MRI. Hum Brain Mapp 2023; 44:4986-5001. [PMID: 37466309 PMCID: PMC10502651 DOI: 10.1002/hbm.26421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 06/14/2023] [Accepted: 07/03/2023] [Indexed: 07/20/2023] Open
Abstract
Magnetic resonance electrical properties tomography (MR-EPT) is a non-invasive measurement technique that derives the electrical properties (EPs, e.g., conductivity or permittivity) of tissues in the radiofrequency range (64 MHz for 1.5 T and 128 MHz for 3 T MR systems). Clinical studies have shown the potential of tissue conductivity as a biomarker. To date, model-based conductivity reconstructions rely on numerical assumptions and approximations, leading to inaccuracies in the reconstructed maps. To address such limitations, we propose an artificial neural network (ANN)-based non-linear conductivity estimator trained on simulated data for conductivity brain imaging. Network training was performed on 201 synthesized T2-weighted spin-echo (SE) data obtained from the finite-difference time-domain (FDTD) electromagnetic (EM) simulation. The dataset was composed of an approximated T2-w SE magnitude and transceive phase information. The proposed method was tested three in-silico and in-vivo on two volunteers and three patients' data. For comparison purposes, various conventional phase-based EPT reconstruction methods were used that ignoreB 1 + magnitude information, such as Savitzky-Golay kernel combined with Gaussian filter (S-G Kernel), phase-based convection-reaction EPT (cr-EPT), magnitude-weighted polynomial-fitting phase-based EPT (Poly-Fit), and integral-based phase-based EPT (Integral-based). From the in-silico experiments, quantitative analysis showed that the proposed method provides more accurate and improved quality (e.g., high structural preservation) conductivity maps compared to conventional reconstruction methods. Representatively, in the healthy brain in-silico phantom experiment, the proposed method yielded mean conductivity values of 1.97 ± 0.20 S/m for CSF, 0.33 ± 0.04 S/m for WM, and 0.52 ± 0.08 S/m for GM, which were closer to the ground-truth conductivity (2.00, 0.30, 0.50 S/m) than the integral-based method (2.56 ± 2.31, 0.39 ± 0.12, 0.68 ± 0.33 S/m). In-vivo ANN-based conductivity reconstructions were also of improved quality compared to conventional reconstructions and demonstrated network generalizability and robustness to in-vivo data and pathologies. The reported in-vivo brain conductivity values were in agreement with literatures. In addition, the proposed method was observed for various SNR levels (SNR levels = 10, 20, 40, and 58) and repeatability conditions (the eight acquisitions with the number of signal averages = 1). The preliminary investigations on brain tumor patient datasets suggest that the network trained on simulated dataset can generalize to unforeseen in-vivo pathologies, thus demonstrating its potential for clinical applications.
Collapse
Affiliation(s)
- Kyu‐Jin Jung
- Department of Electrical and Electronic EngineeringYonsei UniversitySeoulRepublic of Korea
| | - Stefano Mandija
- Computational Imaging Group for MR Therapy and DiagnosticsUniversity Medical Center UtrechtUtrechtThe Netherlands
- Department of RadiotherapyUniversity Medical Center UtrechtUtrechtThe Netherlands
| | - Chuanjiang Cui
- Department of Electrical and Electronic EngineeringYonsei UniversitySeoulRepublic of Korea
| | - Jun‐Hyeong Kim
- Department of Electrical and Electronic EngineeringYonsei UniversitySeoulRepublic of Korea
| | - Mohammed A. Al‐masni
- Department of Artificial IntelligenceCollege of Software & Convergence Technology, Daeyang AI Center, Sejong UniversitySeoulRepublic of Korea
| | - Thierry G. Meerbothe
- Computational Imaging Group for MR Therapy and DiagnosticsUniversity Medical Center UtrechtUtrechtThe Netherlands
- Department of RadiotherapyUniversity Medical Center UtrechtUtrechtThe Netherlands
| | - Mina Park
- Department of Radiology, Gangnam Severance HospitalYonsei University College of MedicineSeoulRepublic of Korea
| | - Cornelis A. T. van den Berg
- Computational Imaging Group for MR Therapy and DiagnosticsUniversity Medical Center UtrechtUtrechtThe Netherlands
- Department of RadiotherapyUniversity Medical Center UtrechtUtrechtThe Netherlands
| | - Dong‐Hyun Kim
- Department of Electrical and Electronic EngineeringYonsei UniversitySeoulRepublic of Korea
| |
Collapse
|
6
|
Jeong H, Andersson J, Hess A, Jezzard P. Effect of subject-specific head morphometry on specific absorption rate estimates in parallel-transmit MRI at 7 T. Magn Reson Med 2023; 89:2376-2390. [PMID: 36656151 PMCID: PMC10952207 DOI: 10.1002/mrm.29589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Revised: 12/02/2022] [Accepted: 12/31/2022] [Indexed: 01/20/2023]
Abstract
PURPOSE To assess the accuracy of morphing an established reference electromagnetic head model to a subject-specific morphometry for the estimation of specific absorption rate (SAR) in 7T parallel-transmit (pTx) MRI. METHODS Synthetic T1 -weighted MR images were created from three high-resolution open-source electromagnetic head voxel models. The accuracy of morphing a "reference" (multimodal image-based detailed anatomical [MIDA]) electromagnetic model into a different subject's native space (Duke and Ella) was compared. Both linear and nonlinear registration methods were evaluated. Maximum 10-g averaged SAR was estimated for circularly polarized mode and for 5000 random RF shim sets in an eight-channel transmit head coil, and comparison made between the morphed MIDA electromagnetic models and the native Duke and Ella electromagnetic models, respectively. RESULTS The averaged error in maximum 10-g averaged SAR estimation across pTx MRI shim sets between the MIDA and the Duke target model was reduced from 17.5% with only rigid-body registration, to 11.8% when affine linear registration was used, and further reduced to 10.7% when nonlinear registration was used. The corresponding figures for the Ella model were 16.7%, 11.2%, and 10.1%. CONCLUSION We found that morphometry accounts for up to half of the subject-specific differences in pTx SAR. Both linear and nonlinear morphing of an electromagnetic model into a target subject improved SAR agreement by better matching head size, morphometry, and position. However, differences remained, likely arising from details in tissue composition estimation. Thus, the uncertainty of the head morphometry and tissue composition may need to be considered separately to achieve personalized SAR estimation.
Collapse
Affiliation(s)
- Hongbae Jeong
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical NeurosciencesUniversity of OxfordOxfordUK
- Athinoula A. Martinos Center for Biomedical Imaging, Department of RadiologyMassachusetts General HospitalBostonMassachusettsUSA
| | - Jesper Andersson
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical NeurosciencesUniversity of OxfordOxfordUK
| | - Aaron Hess
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical NeurosciencesUniversity of OxfordOxfordUK
- Centre for Clinical Magnetic Resonance Research, Department of Cardiovascular MedicineUniversity of OxfordOxfordUK
- British Heart Foundation Centre for Research ExcellenceOxfordUK
| | - Peter Jezzard
- Wellcome Centre for Integrative Neuroimaging, FMRIB Division, Nuffield Department of Clinical NeurosciencesUniversity of OxfordOxfordUK
| |
Collapse
|
7
|
Arduino A, Pennecchi F, Katscher U, Cox M, Zilberti L. Repeatability and Reproducibility Uncertainty in Magnetic Resonance-Based Electric Properties Tomography of a Homogeneous Phantom. Tomography 2023; 9:420-435. [PMID: 36828386 PMCID: PMC9961522 DOI: 10.3390/tomography9010034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Revised: 02/10/2023] [Accepted: 02/14/2023] [Indexed: 02/19/2023] Open
Abstract
Uncertainty assessment is a fundamental step in quantitative magnetic resonance imaging because it makes comparable, in a strict metrological sense, the results of different scans, for example during a longitudinal study. Magnetic resonance-based electric properties tomography (EPT) is a quantitative imaging technique that retrieves, non-invasively, a map of the electric properties inside a human body. Although EPT has been used in some early clinical studies, a rigorous experimental assessment of the associated uncertainty has not yet been performed. This paper aims at evaluating the repeatability and reproducibility uncertainties in phase-based Helmholtz-EPT applied on homogeneous phantom data acquired with a clinical 3 T scanner. The law of propagation of uncertainty is used to evaluate the uncertainty in the estimated conductivity values starting from the uncertainty in the acquired scans, which is quantified through a robust James-Stein shrinkage estimator to deal with the dimensionality of the problem. Repeatable errors are detected in the estimated conductivity maps and are quantified for various values of the tunable parameters of the EPT implementation. The spatial dispersion of the estimated electric conductivity maps is found to be a good approximation of the reproducibility uncertainty, evaluated by changing the position of the phantom after each scan. The results underpin the use of the average conductivity (calculated by weighting the local conductivity values by their uncertainty and taking into account the spatial correlation) as an estimate of the conductivity of the homogeneous phantom.
Collapse
Affiliation(s)
- Alessandro Arduino
- Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135 Torino, Italy
- Correspondence:
| | - Francesca Pennecchi
- Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135 Torino, Italy
| | | | - Maurice Cox
- National Physical Laboratory, Teddington TW11 0LW, UK
| | - Luca Zilberti
- Istituto Nazionale di Ricerca Metrologica, Strada delle Cacce 91, 10135 Torino, Italy
| |
Collapse
|
8
|
Sasaki K, Porter E, Rashed EA, Farrugia L, Schmid G. Measurement and image-based estimation of dielectric properties of biological tissues —past, present, and future—. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac7b64] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 06/22/2022] [Indexed: 12/23/2022]
Abstract
Abstract
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.
Collapse
|
9
|
Zhu D, Schär M, Qin Q. Ultrafast B1 mapping with RF-prepared 3D FLASH acquisition: Correcting the bias due to T 1 -induced k-space filtering effect. Magn Reson Med 2022; 88:757-769. [PMID: 35381114 PMCID: PMC9232926 DOI: 10.1002/mrm.29247] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 03/06/2022] [Accepted: 03/07/2022] [Indexed: 01/25/2023]
Abstract
Purpose The traditional radiofrequency (RF)‐prepared B1 mapping technique consists of one scan with an RF preparation module for flip angle‐encoding and a second scan without this module for normalizing. To reduce the T1‐induced k‐space filtering effect, this method is limited to 2D FLASH acquisition with a two‐parameter method. A novel 3D RF‐prepared three‐parameter method for ultrafast B1‐mapping is proposed to correct the T1‐induced quantification bias. Theory The point spread function analysis of FLASH shows that the prepared longitudinal magnetization before the FLASH acquisition and the image signal obeys a linear (not proportional) relationship. The intercept of the linear function causes the quantification bias and can be captured by a third saturated scan. Methods Using the 2D double‐angle method (DAM) as the reference, a 3D RF‐prepared three‐parameter protocol with 9 s duration was compared with the two‐parameter method, as well as the saturated DAM (SDAM) method, the dual refocusing echo acquisition mode (DREAM) method, and the actual flip‐angle imaging (AFI) method, for B1 mapping of brain, breast, and abdomen with different orientations and shim settings at 3T. Results The 3D RF‐prepared three‐parameter method with complex‐subtraction delivered consistently lower RMS error, error mean, error standard deviation, and higher concordance correlation coefficients values than the two‐parameter method, the three‐parameter method with magnitude‐subtraction, the multi‐slice DREAM and the 3D AFI, and were close to the results of 2D or multi‐slice SDAM. Conclusion The proposed ultrafast 3D RF‐prepared three‐parameter method with complex‐subtraction was demonstrated with high accuracy for B1 mapping of brain, breast, and abdomen.
Collapse
Affiliation(s)
- Dan Zhu
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, USA.,The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Michael Schär
- The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Qin Qin
- F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, USA.,The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| |
Collapse
|
10
|
Leijsen R, van den Berg C, Webb A, Remis R, Mandija S. Combining deep learning and 3D contrast source inversion in MR-based electrical properties tomography. NMR IN BIOMEDICINE 2022; 35:e4211. [PMID: 31840897 PMCID: PMC9285035 DOI: 10.1002/nbm.4211] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2019] [Revised: 10/09/2019] [Accepted: 10/09/2019] [Indexed: 05/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Reijer Leijsen
- Department of Radiology, C.J. Gorter Center for High Field MRILeiden University Medical CenterLeidenThe Netherlands
| | - Cornelis van den Berg
- Department of Radiotherapy, Division of Imaging & OncologyUniversity Medical Center UtrechtUtrechtThe Netherlands
- Computational Imaging Group for MR Diagnostics & Therapy, Center for Image SciencesUtrecht UniversityUtrechtThe Netherlands
| | - Andrew Webb
- Department of Radiology, C.J. Gorter Center for High Field MRILeiden University Medical CenterLeidenThe Netherlands
| | - Rob Remis
- Circuits and Systems Group, Faculty of Electrical Engineering, Mathematics and Computer ScienceDelft University of TechnologyDelftThe Netherlands
| | - Stefano Mandija
- Department of Radiotherapy, Division of Imaging & OncologyUniversity Medical Center UtrechtUtrechtThe Netherlands
- Computational Imaging Group for MR Diagnostics & Therapy, Center for Image SciencesUtrecht UniversityUtrechtThe Netherlands
| |
Collapse
|
11
|
Katscher U, Minhas AS, Katoch N. Magnetic Resonance Electrical Properties Tomography (MREPT). ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2022; 1380:185-202. [DOI: 10.1007/978-3-031-03873-0_8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
12
|
Weingärtner S, Desmond KL, Obuchowski NA, Baessler B, Zhang Y, Biondetti E, Ma D, Golay X, Boss MA, Gunter JL, Keenan KE, Hernando D. Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group. Magn Reson Med 2021; 87:1184-1206. [PMID: 34825741 DOI: 10.1002/mrm.29084] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 10/20/2021] [Accepted: 10/27/2021] [Indexed: 12/26/2022]
Abstract
On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T1 . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.
Collapse
Affiliation(s)
- Sebastian Weingärtner
- Department of Imaging Physics, Delft University of Technology, Delft, The Netherlands
| | - Kimberly L Desmond
- Brain Health Imaging Centre, Centre for Addiction and Mental Health, Toronto, Ontario, Canada.,Department of Psychiatry, University of Toronto, Toronto, Ontario, Canada
| | - Nancy A Obuchowski
- Department of Quantitative Health Sciences, Cleveland Clinic, Cleveland, Ohio, USA
| | - Bettina Baessler
- Institute of Diagnostic and Interventional Radiology, University Hospital Zurich, Zurich, Switzerland
| | - Yuxin Zhang
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin, USA.,Department of Radiology, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | - Emma Biondetti
- Department of Neuroscience, Imaging and Clinical Sciences, D'Annunzio University of Chieti and Pescara, Chieti, Italy
| | - Dan Ma
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio, USA
| | - Xavier Golay
- Brain Repair & Rehabilitation, Institute of Neurology, University College London, United Kingdom.,Gold Standard Phantoms Limited, Rochester, United Kingdom
| | - Michael A Boss
- Center for Research and Innovation, American College of Radiology, Philadelphia, Pennsylvania, USA
| | | | - Kathryn E Keenan
- National Institute of Standards and Technology, Boulder, Colorado, USA
| | - Diego Hernando
- Department of Medical Physics, University of Wisconsin-Madison, Madison, Wisconsin, USA.,Department of Radiology, University of Wisconsin-Madison, Madison, Wisconsin, USA
| | | |
Collapse
|
13
|
Katoch N, Choi BK, Park JA, Ko IO, Kim HJ. Comparison of Five Conductivity Tensor Models and Image Reconstruction Methods Using MRI. Molecules 2021; 26:5499. [PMID: 34576970 PMCID: PMC8467711 DOI: 10.3390/molecules26185499] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Revised: 08/30/2021] [Accepted: 09/07/2021] [Indexed: 11/16/2022] Open
Abstract
Imaging of the electrical conductivity distribution inside the human body has been investigated for numerous clinical applications. The conductivity tensors of biological tissue have been obtained from water diffusion tensors by applying several models, which may not cover the entire phenomenon. Recently, a new conductivity tensor imaging (CTI) method was developed through a combination of B1 mapping, and multi-b diffusion weighted imaging. In this study, we compared the most recent CTI method with the four existing models of conductivity tensors reconstruction. Two conductivity phantoms were designed to evaluate the accuracy of the models. Applied to five human brains, the conductivity tensors using the four existing models and CTI were imaged and compared with the values from the literature. The conductivity image of the phantoms by the CTI method showed relative errors between 1.10% and 5.26%. The images by the four models using DTI could not measure the effects of different ion concentrations subsequently due to prior information of the mean conductivity values. The conductivity tensor images obtained from five human brains through the CTI method were comparable to previously reported literature values. The images by the four methods using DTI were highly correlated with the diffusion tensor images, showing a coefficient of determination (R2) value of 0.65 to 1.00. However, the images by the CTI method were less correlated with the diffusion tensor images and exhibited an averaged R2 value of 0.51. The CTI method could handle the effects of different ion concentrations as well as mobilities and extracellular volume fractions by collecting and processing additional B1 map data. It is necessary to select an application-specific model taking into account the pros and cons of each model. Future studies are essential to confirm the usefulness of these conductivity tensor imaging methods in clinical applications, such as tumor characterization, EEG source imaging, and treatment planning for electrical stimulation.
Collapse
Affiliation(s)
- Nitish Katoch
- Department of Biomedical Engineering, Kyung Hee University, Seoul 02447, Korea; (N.K.); (B.-K.C.)
| | - Bup-Kyung Choi
- Department of Biomedical Engineering, Kyung Hee University, Seoul 02447, Korea; (N.K.); (B.-K.C.)
| | - Ji-Ae Park
- Division of Applied RI, Korea Institute of Radiological and Medical Science, Seoul 01812, Korea;
| | - In-Ok Ko
- Division of Applied RI, Korea Institute of Radiological and Medical Science, Seoul 01812, Korea;
| | - Hyung-Joong Kim
- Department of Biomedical Engineering, Kyung Hee University, Seoul 02447, Korea; (N.K.); (B.-K.C.)
| |
Collapse
|
14
|
Gavazzi S, van Lier ALHMW, Zachiu C, Jansen E, Lagendijk JJW, Stalpers LJA, Crezee H, Kok HP. Advanced patient-specific hyperthermia treatment planning. Int J Hyperthermia 2021; 37:992-1007. [PMID: 32806979 DOI: 10.1080/02656736.2020.1806361] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Hyperthermia treatment planning (HTP) is valuable to optimize tumor heating during thermal therapy delivery. Yet, clinical hyperthermia treatment plans lack quantitative accuracy due to uncertainties in tissue properties and modeling, and report tumor absorbed power and temperature distributions which cannot be linked directly to treatment outcome. Over the last decade, considerable progress has been made to address these inaccuracies and therefore improve the reliability of hyperthermia treatment planning. Patient-specific electrical tissue conductivity derived from MR measurements has been introduced to accurately model the power deposition in the patient. Thermodynamic fluid modeling has been developed to account for the convective heat transport in fluids such as urine in the bladder. Moreover, discrete vasculature trees have been included in thermal models to account for the impact of thermally significant large blood vessels. Computationally efficient optimization strategies based on SAR and temperature distributions have been established to calculate the phase-amplitude settings that provide the best tumor thermal dose while avoiding hot spots in normal tissue. Finally, biological modeling has been developed to quantify the hyperthermic radiosensitization effect in terms of equivalent radiation dose of the combined radiotherapy and hyperthermia treatment. In this paper, we review the present status of these developments and illustrate the most relevant advanced elements within a single treatment planning example of a cervical cancer patient. The resulting advanced HTP workflow paves the way for a clinically feasible and more reliable patient-specific hyperthermia treatment planning.
Collapse
Affiliation(s)
- Soraya Gavazzi
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | | | - Cornel Zachiu
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Eric Jansen
- Amsterdam UMC, Department of Radiation Oncology, Cancer Center Amsterdam, University of Amsterdam, Amsterdam, The Netherlands
| | - Jan J W Lagendijk
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lukas J A Stalpers
- Amsterdam UMC, Department of Radiation Oncology, Cancer Center Amsterdam, University of Amsterdam, Amsterdam, The Netherlands
| | - Hans Crezee
- Amsterdam UMC, Department of Radiation Oncology, Cancer Center Amsterdam, University of Amsterdam, Amsterdam, The Netherlands
| | - H Petra Kok
- Amsterdam UMC, Department of Radiation Oncology, Cancer Center Amsterdam, University of Amsterdam, Amsterdam, The Netherlands
| |
Collapse
|
15
|
Stijnman PRS, Stefano Mandija, Fuchs PS, van den Berg CAT, Remis RF. Transceive phase corrected 2D contrast source inversion-electrical properties tomography. Magn Reson Med 2021; 85:2856-2868. [PMID: 33280166 PMCID: PMC7898605 DOI: 10.1002/mrm.28619] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Revised: 10/08/2020] [Accepted: 11/05/2020] [Indexed: 11/11/2022]
Abstract
PURPOSE To remove the necessity of the tranceive phase assumption for CSI-EPT and show electrical properties maps reconstructed from measured data obtained using a standard 3T birdcage body coil setup. METHODS The existing CSI-EPT algorithm is reformulated to use the transceive phase rather than relying on the transceive phase assumption. Furthermore, the radio frequency (RF)-shield is numerically implemented to accurately model the RF fields inside the MRI scanner. We verify that the reformulated two-dimensional (2D) CSI-EPT algorithm can reconstruct electrical properties maps given 2D electromagnetic simulations. Afterward, the algorithm is tested with three-dimensional (3D) FDTD simulations to investigate if the 2D CSI-EPT can retrieve the electrical properties for 3D RF fields. Finally, an MR experiment at 3T with a phantom is performed. RESULTS From the results of the 2D simulations, it is seen that CSI-EPT can reconstruct the electrical properties using MRI accessible quantities. For 3D simulations, it is observed that the electrical properties are underestimated, nonetheless, CSI-EPT has a lower standard deviation than the standard Helmholtz based methods. Finally, the first CSI-EPT reconstructions based on measured data are presented showing comparable accuracy and precision to reconstructions based on simulated data, and demonstrating the feasibility of CSI-EPT. CONCLUSIONS The CSI-EPT algorithm was rewritten to use MRI accessible quantities. This allows for CSI-EPT to fully exploit the benefits of the higher static magnetic field strengths with a standard quadrature birdcage coil setup.
Collapse
Affiliation(s)
- Peter R. S. Stijnman
- Computational Imaging Group for MRI Diagnostics and TherapyCentre for Image Sciences UMC UtrechtUtrechtThe Netherlands
- Department of Biomedical EngineeringEindhoven University of TechnologyEindhovenThe Netherlands
| | - Stefano Mandija
- Computational Imaging Group for MRI Diagnostics and TherapyCentre for Image Sciences UMC UtrechtUtrechtThe Netherlands
| | - Patrick S. Fuchs
- Circuit & Systems Group of the Electrical EngineeringDelft University of TechnologyDelftThe Netherlands
| | - Cornelis A. T. van den Berg
- Computational Imaging Group for MRI Diagnostics and TherapyCentre for Image Sciences UMC UtrechtUtrechtThe Netherlands
| | - Rob F. Remis
- Circuit & Systems Group of the Electrical EngineeringDelft University of TechnologyDelftThe Netherlands
| |
Collapse
|
16
|
Jang H, Ma Y, Carl M, Jerban S, Chang EY, Du J. Ultrashort echo time Cones double echo steady state (UTE-Cones-DESS) for rapid morphological imaging of short T 2 tissues. Magn Reson Med 2021; 86:881-892. [PMID: 33755258 DOI: 10.1002/mrm.28769] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 02/17/2021] [Accepted: 02/18/2021] [Indexed: 12/12/2022]
Abstract
PURPOSE In this study, we aimed to develop a new technique, ultrashort echo time Cones double echo steady state (UTE-Cones-DESS), for highly efficient morphological imaging of musculoskeletal tissues with short T2 s. We also proposed a novel, single-point Dixon (spDixon)-based approach for fat suppression. METHODS The UTE-Cones-DESS sequence was implemented on a 3T MR system. It uses a short radiofrequency (RF) pulse followed by a pair of balanced spiral-out and spiral-in readout gradients separated by an unbalanced spoiling gradient in-between. The readout gradients are applied immediately before or after the RF pulses to achieve a UTE image (S+ ) and a spin/stimulated echo image (S- ). Weighted echo subtraction between S+ and S- was performed to achieve high contrast specific to short T2 tissues, and spDixon was applied to suppress fat by using the intrinsic complex signal of S+ and S- . Six healthy volunteers and five patients with osteoarthritis were recruited for whole-knee imaging. Additionally, two healthy volunteers were recruited for lower leg imaging. RESULTS The UTE-Cones-DESS sequence allows fast volumetric imaging of musculoskeletal tissues with excellent image contrast for the osteochondral junction, tendons, menisci, and ligaments in the knee joint as well as cortical bone and aponeurosis in the lower leg within 5 min. spDixon yields efficient fat suppression in both S+ and S- images without requiring any additional acquisitions or preparation pulses. CONCLUSION The rapid UTE-Cones-DESS sequence can be used for high contrast morphological imaging of short T2 tissues, providing a new tool to assess their association with musculoskeletal disorders.
Collapse
Affiliation(s)
- Hyungseok Jang
- Department of Radiology, University of California, San Diego, San Diego, CA, USA
| | - Yajun Ma
- Department of Radiology, University of California, San Diego, San Diego, CA, USA
| | | | - Saeed Jerban
- Department of Radiology, University of California, San Diego, San Diego, CA, USA
| | - Eric Y Chang
- Department of Radiology, University of California, San Diego, San Diego, CA, USA.,Radiology Service, Veterans Affairs San Diego Healthcare System, San Diego, CA, USA
| | - Jiang Du
- Department of Radiology, University of California, San Diego, San Diego, CA, USA
| |
Collapse
|
17
|
Leijsen R, Brink W, van den Berg C, Webb A, Remis R. Electrical Properties Tomography: A Methodological Review. Diagnostics (Basel) 2021; 11:176. [PMID: 33530587 PMCID: PMC7910937 DOI: 10.3390/diagnostics11020176] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 01/18/2021] [Accepted: 01/22/2021] [Indexed: 11/25/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Reijer Leijsen
- Department of Radiology, C.J. Gorter Center for High Field MRI, Leiden University Medical Center, 2333ZA Leiden, The Netherlands; (R.L.); (W.B.); (A.W.)
| | - Wyger Brink
- Department of Radiology, C.J. Gorter Center for High Field MRI, Leiden University Medical Center, 2333ZA Leiden, The Netherlands; (R.L.); (W.B.); (A.W.)
| | - Cornelis van den Berg
- Computational Imaging Group for MRI Diagnostics and Therapy, Centre for Image Sciences, University Medical Centre Utrecht, 3508GA Utrecht, The Netherlands;
| | - Andrew Webb
- Department of Radiology, C.J. Gorter Center for High Field MRI, Leiden University Medical Center, 2333ZA Leiden, The Netherlands; (R.L.); (W.B.); (A.W.)
| | - Rob Remis
- Circuits and Systems Group, Faculty of Electrical Engineering, Mathematics and Computes Science, Delft University of Technology, 2628CD Delft, The Netherlands
| |
Collapse
|
18
|
Liu C, Guo L, Li M, Chen H, Jin J, Chen W, Liu F, Crozier S. Divergence-Based Magnetic Resonance Electrical Properties Tomography. IEEE Trans Biomed Eng 2021; 68:192-203. [DOI: 10.1109/tbme.2020.3003460] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
19
|
Iyyakkunnel S, Schäper J, Bieri O. Configuration-based electrical properties tomography. Magn Reson Med 2020; 85:1855-1864. [PMID: 33107082 DOI: 10.1002/mrm.28542] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 09/09/2020] [Accepted: 09/14/2020] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Santhosh Iyyakkunnel
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Jessica Schäper
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| | - Oliver Bieri
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Basel, Switzerland.,Department of Biomedical Engineering, University of Basel, Basel, Switzerland
| |
Collapse
|
20
|
Soullié P, Missoffe A, Ambarki K, Felblinger J, Odille F. MR electrical properties imaging using a generalized image-based method. Magn Reson Med 2020; 85:762-776. [PMID: 32783236 DOI: 10.1002/mrm.28458] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2020] [Revised: 07/08/2020] [Accepted: 07/10/2020] [Indexed: 11/12/2022]
Abstract
PURPOSE To develop a fast and easy-to-use electrical properties tomography (EPT) method based on a single MR scan, avoiding both the need of a B1 -map and transceive phase assumption, and that is robust against noise. THEORY Derived from Maxwell's equations, conductivity, and permittivity are reconstructed from a new partial differential equation involving the product of the RF fields and its derivatives. This also allows us to clarify and revisit the relevance of common assumptions of MREPT. METHODS Our new governing equation is solved using a 3D finite-difference scheme and compared to previous frameworks. The benefits of our method over selected existing MREPT methods are demonstrated for different simulation models, as well as for both an inhomogeneous agar phantom gel and in vivo brain data at 3T. RESULTS Simulation and experimental results are illustrated to highlight the merits of the proposed method over existing methods. We show the validity of our algorithm in versatile configurations, with many transition regions notably. Complex admittivity maps are also provided as a complementary MR contrast. CONCLUSION Because it avoids time-consuming RF field mapping and generalizes the use of standard MR image for electrical properties reconstruction, this contribution is promising as a new step forward for clinical applications.
Collapse
Affiliation(s)
- Paul Soullié
- IADI, INSERM U1254, Université de Lorraine, Nancy, France
| | | | | | - Jacques Felblinger
- IADI, INSERM U1254, Université de Lorraine, Nancy, France.,CIC-IT 1433, INSERM, Université de Lorraine and CHRU de Nancy, Nancy, France
| | - Freddy Odille
- IADI, INSERM U1254, Université de Lorraine, Nancy, France.,CIC-IT 1433, INSERM, Université de Lorraine and CHRU de Nancy, Nancy, France
| |
Collapse
|
21
|
Rashed EA, Gomez-Tames J, Hirata A. Deep Learning-Based Development of Personalized Human Head Model With Non-Uniform Conductivity for Brain Stimulation. IEEE TRANSACTIONS ON MEDICAL IMAGING 2020; 39:2351-2362. [PMID: 31995479 DOI: 10.1109/tmi.2020.2969682] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Electromagnetic stimulation of the human brain is a key tool for neurophysiological characterization and the diagnosis of several neurological disorders. Transcranial magnetic stimulation (TMS) is a commonly used clinical procedure. However, personalized TMS requires a pipeline for individual head model generation to provide target-specific stimulation. This process includes intensive segmentation of several head tissues based on magnetic resonance imaging (MRI), which has significant potential for segmentation error, especially for low-contrast tissues. Additionally, a uniform electrical conductivity is assigned to each tissue in the model, which is an unrealistic assumption based on conventional volume conductor modeling. This study proposes a novel approach for fast and automatic estimation of the electric conductivity in the human head for volume conductor models without anatomical segmentation. A convolutional neural network is designed to estimate personalized electrical conductivity values based on anatomical information obtained from T1- and T2-weighted MRI scans. This approach can avoid the time-consuming process of tissue segmentation and maximize the advantages of position-dependent conductivity assignment based on the water content values estimated from MRI intensity values. The computational results of the proposed approach provide similar but smoother electric field distributions of the brain than that provided by conventional approaches.
Collapse
|
22
|
Hampe N, Katscher U, van den Berg CAT, Tha KK, Mandija S. Investigating the challenges and generalizability of deep learning brain conductivity mapping. ACTA ACUST UNITED AC 2020; 65:135001. [DOI: 10.1088/1361-6560/ab9356] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
23
|
Gavazzi S, van den Berg CAT, Savenije MHF, Kok HP, de Boer P, Stalpers LJA, Lagendijk JJW, Crezee H, van Lier ALHMW. Deep learning-based reconstruction of in vivo pelvis conductivity with a 3D patch-based convolutional neural network trained on simulated MR data. Magn Reson Med 2020; 84:2772-2787. [PMID: 32314825 PMCID: PMC7402024 DOI: 10.1002/mrm.28285] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2019] [Revised: 03/25/2020] [Accepted: 03/26/2020] [Indexed: 12/21/2022]
Abstract
PURPOSE To demonstrate that mapping pelvis conductivity at 3T with deep learning (DL) is feasible. METHODS 210 dielectric pelvic models were generated based on CT scans of 42 cervical cancer patients. For all dielectric models, electromagnetic and MR simulations with realistic accuracy and precision were performed to obtain B 1 + and transceive phase (ϕ± ). Simulated B 1 + and ϕ± served as input to a 3D patch-based convolutional neural network, which was trained in a supervised fashion to retrieve the conductivity. The same network architecture was retrained using only ϕ± in input. Both network configurations were tested on simulated MR data and their conductivity reconstruction accuracy and precision were assessed. Furthermore, both network configurations were used to reconstruct conductivity maps from a healthy volunteer and two cervical cancer patients. DL-based conductivity was compared in vivo and in silico to Helmholtz-based (H-EPT) conductivity. RESULTS Conductivity maps obtained from both network configurations were comparable. Accuracy was assessed by mean error (ME) with respect to ground truth conductivity. On average, ME < 0.1 Sm-1 for all tissues. Maximum MEs were 0.2 Sm-1 for muscle and tumour, and 0.4 Sm-1 for bladder. Precision was indicated with the difference between 90th and 10th conductivity percentiles, and was below 0.1 Sm-1 for fat, bone and muscle, 0.2 Sm-1 for tumour and 0.3 Sm-1 for bladder. In vivo, DL-based conductivity had median values in agreement with H-EPT values, but a higher precision. CONCLUSION Anatomically detailed, noise-robust 3D conductivity maps with good sensitivity to tissue conductivity variations were reconstructed in the pelvis with DL.
Collapse
Affiliation(s)
- Soraya Gavazzi
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands.,Computational Imaging Group for MR diagnostics and therapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Mark H F Savenije
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands.,Computational Imaging Group for MR diagnostics and therapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - H Petra Kok
- Department of Radiation Oncology, Amsterdam University Medical Center, Amsterdam, The Netherlands
| | - Peter de Boer
- Radiotherapy Institute Friesland, Leeuwarden, The Netherlands
| | - Lukas J A Stalpers
- Department of Radiation Oncology, Amsterdam University Medical Center, Amsterdam, The Netherlands
| | - Jan J W Lagendijk
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Hans Crezee
- Department of Radiation Oncology, Amsterdam University Medical Center, Amsterdam, The Netherlands
| | | |
Collapse
|
24
|
Rashed EA, Diao Y, Hirata A. Learning-based estimation of dielectric properties and tissue density in head models for personalized radio-frequency dosimetry. Phys Med Biol 2020; 65:065001. [PMID: 32023556 DOI: 10.1088/1361-6560/ab7308] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Radio-frequency dosimetry is an important process in assessments for human exposure safety and for compliance of related products. Recently, computational human models generated from medical images have often been used for such assessment, especially to consider the inter-subject variability. However, a common procedure to develop personalized models is time consuming because it involves excessive segmentation of several components that represent different biological tissues, which is a major obstacle in the inter-subject variability assessment of radiation safety. Deep learning methods have been shown to be a powerful approach for pattern recognition and signal analysis. Convolutional neural networks with deep architecture are proven robust for feature extraction and image mapping in several biomedical applications. In this study, we develop a learning-based approach for fast and accurate estimation of the dielectric properties and density of tissues directly from magnetic resonance images in a single shot. The smooth distribution of the dielectric properties in head models, which is realized using a process without tissue segmentation, improves the smoothness of the specific absorption rate (SAR) distribution compared with that in the commonly used procedure. The estimated SAR distributions, as well as that averaged over 10 g of tissue in a cubic shape, are found to be highly consistent with those computed using the conventional methods that employ segmentation.
Collapse
Affiliation(s)
- Essam A Rashed
- Department of Electrical and Mechanical Engineering, Nagoya Institute of Technology, Nagoya 466-8555, Japan. Department of Mathematics, Faculty of Science, Suez Canal University, Ismailia 41522, Egypt. Author to whom any correspondence should be addressed
| | | | | |
Collapse
|
25
|
Amouzandeh G, Mentink-Vigier F, Helsper S, Bagdasarian FA, Rosenberg JT, Grant SC. Magnetic resonance electrical property mapping at 21.1 T: a study of conductivity and permittivity in phantoms, ex vivo tissue and in vivo ischemia. Phys Med Biol 2020; 65:055007. [PMID: 31307020 PMCID: PMC7223161 DOI: 10.1088/1361-6560/ab3259] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ghoncheh Amouzandeh
- Department of Physics, Florida State University, Tallahassee, FL, USA
- The National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL, USA
| | | | - Shannon Helsper
- The National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL, USA
- Department of Chemical & Biomedical Engineering, FAMU-FSU College of Engineering, Tallahassee, FL, USA
| | - F. Andrew Bagdasarian
- The National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL, USA
- Department of Chemical & Biomedical Engineering, FAMU-FSU College of Engineering, Tallahassee, FL, USA
| | - Jens T. Rosenberg
- The National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL, USA
| | - Samuel C. Grant
- The National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL, USA
- Department of Chemical & Biomedical Engineering, FAMU-FSU College of Engineering, Tallahassee, FL, USA
| |
Collapse
|
26
|
Gavazzi S, Shcherbakova Y, Bartels LW, Stalpers LJA, Lagendijk JJW, Crezee H, van den Berg CAT, van Lier ALHMW. Transceive phase mapping using the PLANET method and its application for conductivity mapping in the brain. Magn Reson Med 2019; 83:590-607. [PMID: 31483520 PMCID: PMC6900152 DOI: 10.1002/mrm.27958] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2019] [Revised: 07/25/2019] [Accepted: 07/30/2019] [Indexed: 12/23/2022]
Abstract
Purpose To demonstrate feasibility of transceive phase mapping with the PLANET method and its application for conductivity reconstruction in the brain. Methods Accuracy and precision of transceive phase (ϕ±) estimation with PLANET, an ellipse fitting approach to phase‐cycled balanced steady state free precession (bSSFP) data, were assessed with simulations and measurements and compared to standard bSSFP. Measurements were conducted on a homogeneous phantom and in the brain of healthy volunteers at 3 tesla. Conductivity maps were reconstructed with Helmholtz‐based electrical properties tomography. In measurements, PLANET was also compared to a reference technique for transceive phase mapping, i.e., spin echo. Results Accuracy and precision of ϕ± estimated with PLANET depended on the chosen flip angle and TR. PLANET‐based ϕ± was less sensitive to perturbations induced by off‐resonance effects and partial volume (e.g., white matter + myelin) than bSSFP‐based ϕ±. For flip angle = 25° and TR = 4.6 ms, PLANET showed an accuracy comparable to that of reference spin echo but a higher precision than bSSFP and spin echo (factor of 2 and 3, respectively). The acquisition time for PLANET was ~5 min; 2 min faster than spin echo and 8 times slower than bSSFP. However, PLANET simultaneously reconstructed T1, T2, B0 maps besides mapping ϕ±. In the phantom, PLANET‐based conductivity matched the true value and had the smallest spread of the three methods. In vivo, PLANET‐based conductivity was similar to spin echo‐based conductivity. Conclusion Provided that appropriate sequence parameters are used, PLANET delivers accurate and precise ϕ± maps, which can be used to reconstruct brain tissue conductivity while simultaneously recovering T1, T2, and B0 maps.
Collapse
Affiliation(s)
- Soraya Gavazzi
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Yulia Shcherbakova
- Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lambertus W Bartels
- Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands.,Image Sciences Institute, Imaging Division, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Lukas J A Stalpers
- Department of Radiotherapy, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Jan J W Lagendijk
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Hans Crezee
- Department of Radiotherapy, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands.,Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | | |
Collapse
|
27
|
Automated gradient-based electrical properties tomography in the human brain using 7 Tesla MRI. Magn Reson Imaging 2019; 63:258-266. [PMID: 31425805 DOI: 10.1016/j.mri.2019.08.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Revised: 08/01/2019] [Accepted: 08/15/2019] [Indexed: 12/18/2022]
Abstract
Electrical properties of the brain tissues may yield useful biomarkers for neurological disorders and diseases, as well as contribute to safety assurance of ultra-high-field MRI. It has been reported that using B1 maps from a multi-channel RF coil, the spatial variation of the electrical properties can be robustly retrieved. The absolute electrical property values can then be obtained by spatial integration, given that an integration seed point is assigned. In this study, we propose to exploit automatically detected seed points based on tissue piece-wise homogeneity (Helmholtz equation) for spatial integration. Numerical simulations of a numerical brain model and experiments involving 12 healthy volunteers were performed to demonstrate its feasibility and robustness in various noisy conditions and head positions. For in vivo imaging, we consistently observed higher conductivity and permittivity values in the white and gray matter compared to tabulated ex vivo probe measurement results found in the literature, a discrepancy that may be attributed to ex vivo experimental constraints. Our results suggest that the proposed technique produces consistent brain electrical properties in vivo that may contribute to improving diagnostic and therapeutic decisions.
Collapse
|
28
|
Gavazzi S, van den Berg CAT, Sbrizzi A, Kok HP, Stalpers LJA, Lagendijk JJW, Crezee H, van Lier ALHMW. Accuracy and precision of electrical permittivity mapping at 3T: the impact of three B 1 + mapping techniques. Magn Reson Med 2019; 81:3628-3642. [PMID: 30737816 PMCID: PMC6593818 DOI: 10.1002/mrm.27675] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Revised: 01/07/2019] [Accepted: 01/08/2019] [Indexed: 12/29/2022]
Abstract
Purpose To investigate the sequence‐specific impact of B1+ amplitude mapping on the accuracy and precision of permittivity reconstruction at 3T in the pelvic region. Methods B1+ maps obtained with actual flip angle imaging (AFI), Bloch–Siegert (BS), and dual refocusing echo acquisition mode (DREAM) sequences, set to a clinically feasible scan time of 5 minutes, were compared in terms of accuracy and precision with electromagnetic and Bloch simulations and MR measurements. Permittivity maps were reconstructed based on these B1+ maps with Helmholtz‐based electrical properties tomography. Accuracy and precision in permittivity were assessed. A 2‐compartment phantom with properties and size similar to the human pelvis was used for both simulations and measurements. Measurements were also performed on a female volunteer’s pelvis. Results Accuracy was evaluated with noiseless simulations on the phantom. The maximum B1+ bias relative to the true B1+ distribution was 1% for AFI and BS and 6% to 15% for DREAM. This caused an average permittivity bias relative to the true permittivity of 7% to 20% for AFI and BS and 12% to 35% for DREAM. Precision was assessed in MR experiments. The lowest standard deviation in permittivity, found in the phantom for BS, measured 22.4 relative units and corresponded to a standard deviation in B1+ of 0.2% of the B1+ average value. As regards B1+ precision, in vivo and phantom measurements were comparable. Conclusions Our simulation framework quantitatively predicts the different impact of B1+ mapping techniques on permittivity reconstruction and shows high sensitivity of permittivity reconstructions to sequence‐specific bias and noise perturbation in the B1+ map. These findings are supported by the experimental results.
Collapse
Affiliation(s)
- Soraya Gavazzi
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Cornelis A T van den Berg
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands.,Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Alessandro Sbrizzi
- Center for Image Sciences, University Medical Center Utrecht, Utrecht, The Netherlands
| | - H Petra Kok
- Department of Radiation Oncology, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Lukas J A Stalpers
- Department of Radiation Oncology, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | - Jan J W Lagendijk
- Department of Radiotherapy, University Medical Center Utrecht, Utrecht, The Netherlands
| | - Hans Crezee
- Department of Radiation Oncology, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands
| | | |
Collapse
|
29
|
Wang Y, Van de Moortele PF, He B. CONtrast Conformed Electrical Properties Tomography (CONCEPT) Based on Multi- Channel Transmission and Alternating Direction Method of Multipliers. IEEE TRANSACTIONS ON MEDICAL IMAGING 2019; 38:349-359. [PMID: 30106715 PMCID: PMC6372102 DOI: 10.1109/tmi.2018.2865121] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
In magnetic resonance-based electrical properties tomography (EPT), circularly polarized magnetic field B1 from a transmit radiofrequency (RF) coil is measured and utilized to infer the electrical conductivity and permittivity of biological tissues. Compared with a quadrature RF coil, a multi-channel transmit coil provides a plurality of unique transmit B1 patterns that help to alleviate the under-determinedness of EPT reconstruction problem, and it also allows to circumvent the "transceive phase assumption" that fails at ultra-high-field MRI. Here, a new approach, contrast conformed electrical properties tomography or CONCEPT, is proposed based on the multi-channel transmission that retrieves electrical properties (EPs) by solving a linear partial differential equation with discriminated L1 and L2 norm regularization informed by intermediate EP gradient. The theory of CONCEPT and a fast reconstruction algorithm based on the alternating direction method of multipliers are described and evaluated using numerical simulations, phantom experiment, and analysis of in vivo human brain data at 7 T MRI. Compared with the multi-channel gradient-based EPT (gEPT) method, this new technology does not require receive- B1 sensitivity profiles and does not rely on symmetry assumption regarding RF coil design and imaged target. Moreover, it is not dependent on external prior information, such as integration seed point or anatomical MRI, which can be sources of bias in reconstructed EP values. By deriving EPs from transmit B1 profiles only, CONCEPT can be used with RF coils that include receive-only arrays with large channel count which can, in turn, offer substantial gains in signal-to-noise ratio. It also holds potentials to image unsymmetrical body organs and diseased brain. CONCEPT provides solutions for the practical problems during the implementation of gEPT, thus representing a more generalized framework in the context of multi-channel RF transmission.
Collapse
Affiliation(s)
- Yicun Wang
- Department of Biomedical Engineering, University of Minnesota, MN 55455, USA
| | | | - Bin He
- Department of Biomedical Engineering, University of Minnesota, MN 55455, USA; and is with the Department of Biomedical Engineering, Carnegie Mellon University, PA 15213, USA. ()
| |
Collapse
|
30
|
Yildiz G, Ider YZ. Use of dielectric padding to eliminate low convective field artifact in cr-MREPT conductivity images. Magn Reson Med 2019; 81:3168-3184. [PMID: 30693565 DOI: 10.1002/mrm.27648] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Revised: 12/05/2018] [Accepted: 12/05/2018] [Indexed: 11/07/2022]
Abstract
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 (B1 + ), 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 BaTiO3 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 Bz 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 L2 -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.
Collapse
Affiliation(s)
- Gulsah Yildiz
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey
| | - Yusuf Ziya Ider
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey
| |
Collapse
|
31
|
Cardiac Segmentation from LGE MRI Using Deep Neural Network Incorporating Shape and Spatial Priors. LECTURE NOTES IN COMPUTER SCIENCE 2019. [DOI: 10.1007/978-3-030-32245-8_62] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
32
|
Wang Y, Shao Q, Van de Moortele PF, Racila E, Liu J, Bischof J, He B. Mapping electrical properties heterogeneity of tumor using boundary informed electrical properties tomography (BIEPT) at 7T. Magn Reson Med 2019; 81:393-409. [PMID: 30230603 PMCID: PMC6258314 DOI: 10.1002/mrm.27414] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2018] [Revised: 05/30/2018] [Accepted: 05/31/2018] [Indexed: 11/06/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yicun Wang
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455
| | - Qi Shao
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455
| | | | - Emilian Racila
- Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, MN 55455
| | - Jiaen Liu
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455
- National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD 20892
| | - John Bischof
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455
- Department of Mechanical Engineering, University of Minnesota, Minneapolis, MN 55455
| | - Bin He
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN 55455
- Department of Biomedical Engineering, Carnegie Mellon University, PA 15213, USA
| |
Collapse
|
33
|
Guo L, Jin J, Li M, Wang Y, Liu CY, Liu F, Crozier S. Reference-Based Integral MR-EPT:Simulation and Experiment Studies on the 9.4T MRI. IEEE Trans Biomed Eng 2018; 66:1832-1843. [PMID: 30403619 DOI: 10.1109/tbme.2018.2879667] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Current integral-equation (IE) based MR electrical properties tomography (EPT) methods utilize simulated incident radio-frequency (RF) fields, which are inaccurate and lead to reconstruction errors. To improve the accuracy and practicability of IE-based MR-EPT methods, a new approach is presented that obtains the incident fields using reference subjects and RF field mapping techniques. The Incident field approximation (IFA) is first demonstrated in this paper. This approximation assumes that two imaged subjects with similar coil/subject interactions will have similar incident RF fields, thus one can feed the estimation of the incident fields within the imaged subject into the calculation of those within a homogeneous subject (reference subject). This is done by measuring the total RF fields ( field) of the reference using field mapping techniques, using the known EPs of the reference subject and by rearranging Ampere's Law and the integral equations. The calculated incident RF fields are then used to reconstruct the EPs' distribution with a three-dimensional (3D) integral-based MR-EPT method. Numerical simulation results indicated that the incident RF fields obtained from the reference subject provide accurate 3D reconstruction of EPs with less than 16% root mean square error (RMSE) in noise-free scenario while the conventional IE method had more than 28% RMSE. The phantom-based experiments at 9.4T MRI system have also been conducted to evaluate the performance of the proposed method and the results indicated that the proposed method achieved desirable robustness against the noise in practical scenario with less than 21% RMSE while the conventional differential equation-based method showed worse than 37% RMSE.
Collapse
|
34
|
Ozdemir S, Ider YZ. bSSFP phase correction and its use in magnetic resonance electrical properties tomography. Magn Reson Med 2018; 81:934-946. [PMID: 30357891 DOI: 10.1002/mrm.27446] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2018] [Revised: 06/06/2018] [Accepted: 06/12/2018] [Indexed: 11/05/2022]
Abstract
PURPOSE Balanced steady-state free precession (bSSFP) sequence is widely used because of its high SNR and high speed. However, bSSFP images suffer from "banding artifact" caused by B0 inhomogeneity. In this article, we propose a method to remove this artifact in bSSFP phase images and investigate the usage of the corrected phase images in phase-based magnetic resonance electrical properties tomography (MREPT). THEORY AND METHODS Two bSSFP phase images, obtained with different excitation frequencies, are collaged to get rid of the regions containing banding artifacts. Phase of the collaged bSSFP image is the sum of the transceive phase of the RF system and an error term that depends on B0 and T2 . By using B0 and T2 maps, this error is eliminated from bSSFP phase images by using pixel-wise corrections. Conductivity maps are obtained from the uncorrected and the corrected phase images using the phase-based cr-MREPT method. RESULTS Phantom and human experiment results of the proposed method are illustrated for both phase images and conductivity maps. It is shown that uncorrected phase images yield unacceptable conductivity images. When only B0 information is used for phase correction conductivity, reconstructions are substantially improved, and yet T2 information is still needed to fully recover accurate and undistorted conductivity images. CONCLUSIONS With the proposed technique, B0 sensitivity of the bSSFP phase images can be removed by using B0 and T2 maps. It is also shown that corrected bSSFP phase images are of sufficient quality to be used in conductivity imaging.
Collapse
Affiliation(s)
- Safa Ozdemir
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey
| | - Yusuf Ziya Ider
- Department of Electrical and Electronics Engineering, Bilkent University, Ankara, Turkey
| |
Collapse
|
35
|
Hancu I, Liu J, Hua Y, Lee SK. Electrical properties tomography: Available contrast and reconstruction capabilities. Magn Reson Med 2018; 81:803-810. [PMID: 30325052 DOI: 10.1002/mrm.27453] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Revised: 06/10/2018] [Accepted: 06/24/2018] [Indexed: 12/29/2022]
Abstract
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.
Collapse
Affiliation(s)
| | - Jiaen Liu
- National Institute of Neurological, Disorders and Stroke, National Institutes of Health, Bethesda, Maryland
| | - Yihe Hua
- GE Global Research, Niskayuna, New York
| | - Seung-Kyun Lee
- IBS Center for Neuroscience Imaging Research, Suwon, Republic of Korea
- Department of Biomedical Engineering, Sungkyunkwan University, Suwon, Republic of Korea
| |
Collapse
|
36
|
Hampe N, Herrmann M, Amthor T, Findeklee C, Doneva M, Katscher U. Dictionary-based electric properties tomography. Magn Reson Med 2018; 81:342-349. [PMID: 30246342 DOI: 10.1002/mrm.27401] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Revised: 05/18/2018] [Accepted: 05/21/2018] [Indexed: 12/24/2022]
Abstract
PURPOSE To develop and validate a new algorithm called "dictionary-based electric properties tomography" (dbEPT) for deriving tissue electric properties from measured B1 maps. METHODS Inspired by Magnetic Resonance fingerprinting, dbEPT uses a dictionary of local patterns ("atoms") of B1 maps and corresponding electric properties distributions, derived from electromagnetic field simulations. For reconstruction, a pattern from a measured B1 map is compared with the B1 atoms of the dictionary. The B1 atom showing the best match with the measured B1 pattern yields the optimum electric properties pattern that is chosen for reconstruction. Matching was performed through machine learning algorithms. Two dictionaries, using transmit and transceive phases, were evaluated. The spatial distribution of local matching distance between optimal atom and measured pattern yielded a reconstruction reliability map. The method was applied to reconstruct conductivity of 4 volunteers' brains. A conventional, Helmholtz-based Electric properties tomography (EPT) reconstruction was performed for reference. Noise performance was studied through phantom simulations. RESULTS Quantitative values of conductivity agree with literature values. Results of the 2 dictionaries exhibit only minor differences. Somewhat larger differences are visible between dbEPT and Helmholtz-based EPT. Quantified by the correlation between conductivity and anatomic images, dbEPT depicts brain details more clearly than Helmholtz-based EPT. Matching distance is minimal in homogeneous brain ventricles and increases with tissue heterogeneity. Central processing unit time was approximately 2 minutes per dictionary training and 3 minutes per brain conductivity reconstruction using standard hardware equipment. CONCLUSION A new, dictionary-based approach for reconstructing electric properties is presented. Its conductivity reconstruction is able to overcome the EPT transceive-phase problem.
Collapse
Affiliation(s)
| | - Max Herrmann
- University of Applied Sciences, Hamburg, Germany
| | | | | | | | | |
Collapse
|
37
|
Leijsen RL, Brink WM, van den Berg CAT, Webb AG, Remis RF. 3-D Contrast Source Inversion-Electrical Properties Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2018; 37:2080-2089. [PMID: 29994520 DOI: 10.1109/tmi.2018.2816125] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Contrast source inversion-electrical properties tomography (CSI-EPT) is an iterative reconstruction method to retrieve the electrical properties (EPs) of tissues from magnetic resonance data. The method is based on integral representations of the electromagnetic field and has been shown to allow EP reconstructions of small structures as well as tissue boundaries with compelling accuracy. However, to date, the CSI-EPT has been implemented for 2-D configurations only, which limits its applicability. In this paper, a full 3-D extension of the CSI-EPT method is presented, to enable CSI-EPT to be applied to realistic 3-D scenarios. Here, we demonstrate a proof-of-principle of 3-D CSI-EPT and present the reconstructions of a 3-D abdominal body section and a 3-D head model using different settings of the transmit coil. Numerical results show that the full 3-D approach yields accurate reconstructions of the EPs, even at tissue boundaries and is most accurate in regions where the absolute value of the electric field is highest.
Collapse
|
38
|
Ariturk G, Ider YZ. Optimal multichannel transmission for improved cr-MREPT. ACTA ACUST UNITED AC 2018; 63:045001. [DOI: 10.1088/1361-6560/aaa732] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|