Magnetic resonance electrical properties tomography for small anomalies using boundary conditions

A simulation study

Joonsung Lee, Narae Choi, Jin Keun Seo, Donghyun Kim

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Purpose: Magnetic resonance electrical property tomography (MREPT) is an emerging imaging modality using measured B1 maps from magnetic resonance imaging (MRI) to measure a distribution of electric conductivity and permittivity of the subject at the Larmor frequency. Conventional MREPT approaches at single transmit channel system using the Helmholtz equation rely on an assumption that conductivity and permittivity of the subject are locally homogeneous. For small tissue structures and tissue boundaries, in which the assumption of locally homogeneous conductivity and permittivity does not hold, the reconstructed conductivity values deviated from the actual values, so called "Boundary Artifacts." The aim of this study is to propose new reconstruction processes based on time-harmonic Maxwell's equations to reconstruct conductivity for small tissue structures and tissue boundaries. Methods: Instead of removing the electric fields from the equations as done in the Helmholtz equation, three key identities of circularly polarized and longitudinal components of electric fields, circularly polarized component of magnetic fields, and electric properties from time-harmonic Maxwell's equations are derived. Based on the three key identities, the proposed reconstruction methods determine conductivity, permittivity, and circularly polarized component and longitudinal component of electric fields using the measured H1 +. In each iterative step, estimated conductivity, permittivity, electric fields, and artifact-free mask region, Ω, where the contribution of the boundary artifacts is small, were updated. Using the estimated values in the artifact-free mask region as boundary conditions, the estimates beyond the mask region were updated. EM simulations were performed on three types of numerical phantoms with very small regions of homogeneous conductivity and permittivity. The performance of the proposed methods was evaluated using the simulated electric and magnetic fields. Results: For the numerical simulation model, the proposed methods significantly reduced the boundary artifacts compared to conventional methods using Helmholtz equations. In addition, previous methods using the Helmholtz equation could measure conductivity of only large anomalies, but the proposed method can measure the conductivity of the small compartments whose size is 2-3 voxels. The proposed approaches are compatible with spatial filtering which can be used to reduce noise. If a good image segmentation is available as a prior information, better initial boundary conditions can be estimated, and thus the proposed approach can be more accurate for small tissue structures. Conclusions: The proposed reconstruction method not only determines electrical properties, but also circularly polarized component and longitudinal component of electric fields using an iterative process. The proposed method can quantitatively detect the conductivity of the small anomalies better than conventional methods.

Original languageEnglish
Pages (from-to)4773-4785
Number of pages13
JournalMedical physics
Volume44
Issue number9
DOIs
Publication statusPublished - 2017 Sep 1

Fingerprint

Magnetic Resonance Spectroscopy
Tomography
Artifacts
Masks
Electric Conductivity
Magnetic Fields
Noise
Magnetic Resonance Imaging

All Science Journal Classification (ASJC) codes

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

@article{d3b56ed9138f4745a5c30841a463e9ae,
title = "Magnetic resonance electrical properties tomography for small anomalies using boundary conditions: A simulation study",
abstract = "Purpose: Magnetic resonance electrical property tomography (MREPT) is an emerging imaging modality using measured B1 maps from magnetic resonance imaging (MRI) to measure a distribution of electric conductivity and permittivity of the subject at the Larmor frequency. Conventional MREPT approaches at single transmit channel system using the Helmholtz equation rely on an assumption that conductivity and permittivity of the subject are locally homogeneous. For small tissue structures and tissue boundaries, in which the assumption of locally homogeneous conductivity and permittivity does not hold, the reconstructed conductivity values deviated from the actual values, so called {"}Boundary Artifacts.{"} The aim of this study is to propose new reconstruction processes based on time-harmonic Maxwell's equations to reconstruct conductivity for small tissue structures and tissue boundaries. Methods: Instead of removing the electric fields from the equations as done in the Helmholtz equation, three key identities of circularly polarized and longitudinal components of electric fields, circularly polarized component of magnetic fields, and electric properties from time-harmonic Maxwell's equations are derived. Based on the three key identities, the proposed reconstruction methods determine conductivity, permittivity, and circularly polarized component and longitudinal component of electric fields using the measured H1 +. In each iterative step, estimated conductivity, permittivity, electric fields, and artifact-free mask region, Ω, where the contribution of the boundary artifacts is small, were updated. Using the estimated values in the artifact-free mask region as boundary conditions, the estimates beyond the mask region were updated. EM simulations were performed on three types of numerical phantoms with very small regions of homogeneous conductivity and permittivity. The performance of the proposed methods was evaluated using the simulated electric and magnetic fields. Results: For the numerical simulation model, the proposed methods significantly reduced the boundary artifacts compared to conventional methods using Helmholtz equations. In addition, previous methods using the Helmholtz equation could measure conductivity of only large anomalies, but the proposed method can measure the conductivity of the small compartments whose size is 2-3 voxels. The proposed approaches are compatible with spatial filtering which can be used to reduce noise. If a good image segmentation is available as a prior information, better initial boundary conditions can be estimated, and thus the proposed approach can be more accurate for small tissue structures. Conclusions: The proposed reconstruction method not only determines electrical properties, but also circularly polarized component and longitudinal component of electric fields using an iterative process. The proposed method can quantitatively detect the conductivity of the small anomalies better than conventional methods.",
author = "Joonsung Lee and Narae Choi and Seo, {Jin Keun} and Donghyun Kim",
year = "2017",
month = "9",
day = "1",
doi = "10.1002/mp.12343",
language = "English",
volume = "44",
pages = "4773--4785",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "9",

}

Magnetic resonance electrical properties tomography for small anomalies using boundary conditions : A simulation study. / Lee, Joonsung; Choi, Narae; Seo, Jin Keun; Kim, Donghyun.

In: Medical physics, Vol. 44, No. 9, 01.09.2017, p. 4773-4785.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Magnetic resonance electrical properties tomography for small anomalies using boundary conditions

T2 - A simulation study

AU - Lee, Joonsung

AU - Choi, Narae

AU - Seo, Jin Keun

AU - Kim, Donghyun

PY - 2017/9/1

Y1 - 2017/9/1

N2 - Purpose: Magnetic resonance electrical property tomography (MREPT) is an emerging imaging modality using measured B1 maps from magnetic resonance imaging (MRI) to measure a distribution of electric conductivity and permittivity of the subject at the Larmor frequency. Conventional MREPT approaches at single transmit channel system using the Helmholtz equation rely on an assumption that conductivity and permittivity of the subject are locally homogeneous. For small tissue structures and tissue boundaries, in which the assumption of locally homogeneous conductivity and permittivity does not hold, the reconstructed conductivity values deviated from the actual values, so called "Boundary Artifacts." The aim of this study is to propose new reconstruction processes based on time-harmonic Maxwell's equations to reconstruct conductivity for small tissue structures and tissue boundaries. Methods: Instead of removing the electric fields from the equations as done in the Helmholtz equation, three key identities of circularly polarized and longitudinal components of electric fields, circularly polarized component of magnetic fields, and electric properties from time-harmonic Maxwell's equations are derived. Based on the three key identities, the proposed reconstruction methods determine conductivity, permittivity, and circularly polarized component and longitudinal component of electric fields using the measured H1 +. In each iterative step, estimated conductivity, permittivity, electric fields, and artifact-free mask region, Ω, where the contribution of the boundary artifacts is small, were updated. Using the estimated values in the artifact-free mask region as boundary conditions, the estimates beyond the mask region were updated. EM simulations were performed on three types of numerical phantoms with very small regions of homogeneous conductivity and permittivity. The performance of the proposed methods was evaluated using the simulated electric and magnetic fields. Results: For the numerical simulation model, the proposed methods significantly reduced the boundary artifacts compared to conventional methods using Helmholtz equations. In addition, previous methods using the Helmholtz equation could measure conductivity of only large anomalies, but the proposed method can measure the conductivity of the small compartments whose size is 2-3 voxels. The proposed approaches are compatible with spatial filtering which can be used to reduce noise. If a good image segmentation is available as a prior information, better initial boundary conditions can be estimated, and thus the proposed approach can be more accurate for small tissue structures. Conclusions: The proposed reconstruction method not only determines electrical properties, but also circularly polarized component and longitudinal component of electric fields using an iterative process. The proposed method can quantitatively detect the conductivity of the small anomalies better than conventional methods.

AB - Purpose: Magnetic resonance electrical property tomography (MREPT) is an emerging imaging modality using measured B1 maps from magnetic resonance imaging (MRI) to measure a distribution of electric conductivity and permittivity of the subject at the Larmor frequency. Conventional MREPT approaches at single transmit channel system using the Helmholtz equation rely on an assumption that conductivity and permittivity of the subject are locally homogeneous. For small tissue structures and tissue boundaries, in which the assumption of locally homogeneous conductivity and permittivity does not hold, the reconstructed conductivity values deviated from the actual values, so called "Boundary Artifacts." The aim of this study is to propose new reconstruction processes based on time-harmonic Maxwell's equations to reconstruct conductivity for small tissue structures and tissue boundaries. Methods: Instead of removing the electric fields from the equations as done in the Helmholtz equation, three key identities of circularly polarized and longitudinal components of electric fields, circularly polarized component of magnetic fields, and electric properties from time-harmonic Maxwell's equations are derived. Based on the three key identities, the proposed reconstruction methods determine conductivity, permittivity, and circularly polarized component and longitudinal component of electric fields using the measured H1 +. In each iterative step, estimated conductivity, permittivity, electric fields, and artifact-free mask region, Ω, where the contribution of the boundary artifacts is small, were updated. Using the estimated values in the artifact-free mask region as boundary conditions, the estimates beyond the mask region were updated. EM simulations were performed on three types of numerical phantoms with very small regions of homogeneous conductivity and permittivity. The performance of the proposed methods was evaluated using the simulated electric and magnetic fields. Results: For the numerical simulation model, the proposed methods significantly reduced the boundary artifacts compared to conventional methods using Helmholtz equations. In addition, previous methods using the Helmholtz equation could measure conductivity of only large anomalies, but the proposed method can measure the conductivity of the small compartments whose size is 2-3 voxels. The proposed approaches are compatible with spatial filtering which can be used to reduce noise. If a good image segmentation is available as a prior information, better initial boundary conditions can be estimated, and thus the proposed approach can be more accurate for small tissue structures. Conclusions: The proposed reconstruction method not only determines electrical properties, but also circularly polarized component and longitudinal component of electric fields using an iterative process. The proposed method can quantitatively detect the conductivity of the small anomalies better than conventional methods.

UR - http://www.scopus.com/inward/record.url?scp=85029281365&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85029281365&partnerID=8YFLogxK

U2 - 10.1002/mp.12343

DO - 10.1002/mp.12343

M3 - Article

VL - 44

SP - 4773

EP - 4785

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 9

ER -