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
|
Ö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
|
3
|
Cao J, Ball I, Humburg P, Dokos S, Rae C. Repeatability of brain phase-based magnetic resonance electric properties tomography methods and effect of compressed SENSE and RF shimming. Phys Eng Sci Med 2023; 46:753-766. [PMID: 36995580 DOI: 10.1007/s13246-023-01248-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2022] [Accepted: 03/19/2023] [Indexed: 03/31/2023]
Abstract
Magnetic resonance electrical properties tomography (MREPT) is an emerging imaging modality to noninvasively measure tissue conductivity and permittivity. Implementation of MREPT in the clinic requires repeatable measurements at a short scan time and an appropriate protocol. The aim of this study was to investigate the repeatability of conductivity measurements using phase-based MREPT and the effects of compressed SENSE (CS), and RF shimming on the precision of conductivity measurements. Conductivity measurements using turbo spin echo (TSE) and three-dimensional balanced fast field echo (bFFE) with CS factors were repeatable. Conductivity measurement using bFFE phase showed smaller mean and variance that those measured by TSE. The conductivity measurements using bFFE showed minimal deviation with CS factors up to 8, with deviation increasing at CS factors > 8. Subcortical structures produced less consistent measurements than cortical parcellations at higher CS factors. RF shimming using full slice coverage 2D dual refocusing echo acquisition mode (DREAM) and full coverage 3D dual TR approaches further improved measurement precision. BFFE is a more optimal sequence than TSE for phase-based MREPT in brain. Depending on the area of the brain being measured, the scan can be safely accelerated with compressed SENSE without sacrifice of precision, offering the potential to employ MREPT in clinical research and applications. RF shimming with better field mapping further improves precision of the conductivity measures.
Collapse
Affiliation(s)
- Jun Cao
- Neuroscience Research Australia, 139 Barker St, Randwick, NSW, 2031, Australia
- School of Biomedical Sciences, The University of New South Wales, Kensington, NSW, 2052, Australia
| | - Iain Ball
- Philips Australia & New Zealand, North Ryde, NSW, 2113, Australia
| | - Peter Humburg
- Neuroscience Research Australia, 139 Barker St, Randwick, NSW, 2031, Australia
- Mark Wainwright Analytical Centre, Stats Central, The University of New South Wales, Kensington, NSW, 2052, Australia
| | - Socrates Dokos
- Graduate School of Biomedical Engineering, The University of New South Wales, Kensington, NSW, 2052, Australia
| | - Caroline Rae
- Neuroscience Research Australia, 139 Barker St, Randwick, NSW, 2031, Australia.
- School of Psychology, The University of New South Wales, Kensington, NSW, 2052, Australia.
| |
Collapse
|
4
|
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
|
5
|
Eda N, Fushimi M, Hasegawa K, Nara T. A Method for Electrical Property Tomography Based on a Three-Dimensional Integral Representation of the Electric Field. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:1400-1409. [PMID: 34968176 DOI: 10.1109/tmi.2021.3139455] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Magnetic resonance electrical properties tomography (MREPT) noninvasively reconstructs high-resolution electrical property (EP) maps using MRI scanners and is useful for diagnosing cancerous tissues. However, conventional MREPT methods have limitations: sensitivity to noise in the numerical Laplacian operation, difficulty in reconstructing three-dimensional (3D) EPs and convergence not guaranteed in the iterative process. We propose a novel, iterative 3D reconstruction MREPT method without a numerical Laplacian operation. We derive an integral representation of the electric field using its Helmholtz decomposition with Maxwell's equations, under the assumption that the EPs are known on the boundary of the region of interest with the approximation that the unmeasurable magnetic field components are zero. Then, we solve the simultaneous equations composed of the integral representation and Ampere's law using a convex projection algorithm whose convergence is theoretically guaranteed. The efficacy of the proposed method was validated through numerical simulations and a phantom experiment. The results showed that this method is effective in reconstructing 3D EPs and is robust to noise. It was also shown that our proposed method with the unmeasurable component H- enhances the accuracy of the EPs in a background and that with all the components of the magnetic field reduces the artifacts at the center of the slices except when all the components of the electric field are close to zero.
Collapse
|
6
|
Keskin K, Yilmaz U, Cukur T. Constrained Ellipse Fitting for Efficient Parameter Mapping With Phase-Cycled bSSFP MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:14-26. [PMID: 34351856 DOI: 10.1109/tmi.2021.3102852] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Balanced steady-state free precession (bSSFP) imaging enables high scan efficiency in MRI, but differs from conventional sequences in terms of elevated sensitivity to main field inhomogeneity and nonstandard [Formula: see text]-weighted tissue contrast. To address these limitations, multiple bSSFP images of the same anatomy are commonly acquired with a set of different RF phase-cycling increments. Joint processing of phase-cycled acquisitions serves to mitigate sensitivity to field inhomogeneity. Recently phase-cycled bSSFP acquisitions were also leveraged to estimate relaxation parameters based on explicit signal models. While effective, these model-based methods often involve a large number of acquisitions (N ≈ 10-16), degrading scan efficiency. Here, we propose a new constrained ellipse fitting method (CELF) for parameter estimation with improved efficiency and accuracy in phase-cycled bSSFP MRI. CELF is based on the elliptical signal model framework for complex bSSFP signals; and it introduces geometrical constraints on ellipse properties to improve estimation efficiency, and dictionary-based identification to improve estimation accuracy. CELF generates maps of [Formula: see text], [Formula: see text], off-resonance and on-resonant bSSFP signal by employing a separate [Formula: see text] map to mitigate sensitivity to flip angle variations. Our results indicate that CELF can produce accurate off-resonance and banding-free bSSFP maps with as few as N = 4 acquisitions, while estimation accuracy for relaxation parameters is notably limited by biases from microstructural sensitivity of bSSFP imaging.
Collapse
|
7
|
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
|
8
|
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
|
9
|
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
|