Gradient nonlinearity calibration and correction for a compact, asymmetric magnetic resonance imaging gradient system

S. Tao, Joshua D Trazasko, J. L. Gunter, P. T. Weavers, Yunhong Shu, John III Huston, S. K. Lee, E. T. Tan, Matthew A Bernstein

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

Due to engineering limitations, the spatial encoding gradient fields in conventional magnetic resonance imaging cannot be perfectly linear and always contain higher-order, nonlinear components. If ignored during image reconstruction, gradient nonlinearity (GNL) manifests as image geometric distortion. Given an estimate of the GNL field, this distortion can be corrected to a degree proportional to the accuracy of the field estimate. The GNL of a gradient system is typically characterized using a spherical harmonic polynomial model with model coefficients obtained from electromagnetic simulation. Conventional whole-body gradient systems are symmetric in design; typically, only odd-order terms up to the 5th-order are required for GNL modeling. Recently, a high-performance, asymmetric gradient system was developed, which exhibits more complex GNL that requires higher-order terms including both odd- and even-orders for accurate modeling. This work characterizes the GNL of this system using an iterative calibration method and a fiducial phantom used in ADNI (Alzheimer's Disease Neuroimaging Initiative). The phantom was scanned at different locations inside the 26 cm diameter-spherical-volume of this gradient, and the positions of fiducials in the phantom were estimated. An iterative calibration procedure was utilized to identify the model coefficients that minimize the mean-squared-error between the true fiducial positions and the positions estimated from images corrected using these coefficients. To examine the effect of higher-order and even-order terms, this calibration was performed using spherical harmonic polynomial of different orders up to the 10th-order including even- and odd-order terms, or odd-order only. The results showed that the model coefficients of this gradient can be successfully estimated. The residual root-mean-squared-error after correction using up to the 10th-order coefficients was reduced to 0.36 mm, yielding spatial accuracy comparable to conventional whole-body gradients. The even-order terms were necessary for accurate GNL modeling. In addition, the calibrated coefficients improved image geometric accuracy compared with the simulation-based coefficients.

Original languageEnglish (US)
Pages (from-to)N18-N31
JournalPhysics in Medicine and Biology
Volume62
Issue number2
DOIs
StatePublished - Jan 21 2017

Fingerprint

Calibration
Magnetic Resonance Imaging
Computer-Assisted Image Processing
Electromagnetic Phenomena
Statistical Models
Neuroimaging
Alzheimer Disease

Keywords

  • asymmetric gradient
  • compact 3T
  • gradient nonlinearity
  • head-only MRI system
  • image geometric distortion

ASJC Scopus subject areas

  • Radiological and Ultrasound Technology
  • Radiology Nuclear Medicine and imaging

Cite this

Gradient nonlinearity calibration and correction for a compact, asymmetric magnetic resonance imaging gradient system. / Tao, S.; Trazasko, Joshua D; Gunter, J. L.; Weavers, P. T.; Shu, Yunhong; Huston, John III; Lee, S. K.; Tan, E. T.; Bernstein, Matthew A.

In: Physics in Medicine and Biology, Vol. 62, No. 2, 21.01.2017, p. N18-N31.

Research output: Contribution to journalArticle

@article{46525ad7e7c54ff2929af9076237867d,
title = "Gradient nonlinearity calibration and correction for a compact, asymmetric magnetic resonance imaging gradient system",
abstract = "Due to engineering limitations, the spatial encoding gradient fields in conventional magnetic resonance imaging cannot be perfectly linear and always contain higher-order, nonlinear components. If ignored during image reconstruction, gradient nonlinearity (GNL) manifests as image geometric distortion. Given an estimate of the GNL field, this distortion can be corrected to a degree proportional to the accuracy of the field estimate. The GNL of a gradient system is typically characterized using a spherical harmonic polynomial model with model coefficients obtained from electromagnetic simulation. Conventional whole-body gradient systems are symmetric in design; typically, only odd-order terms up to the 5th-order are required for GNL modeling. Recently, a high-performance, asymmetric gradient system was developed, which exhibits more complex GNL that requires higher-order terms including both odd- and even-orders for accurate modeling. This work characterizes the GNL of this system using an iterative calibration method and a fiducial phantom used in ADNI (Alzheimer's Disease Neuroimaging Initiative). The phantom was scanned at different locations inside the 26 cm diameter-spherical-volume of this gradient, and the positions of fiducials in the phantom were estimated. An iterative calibration procedure was utilized to identify the model coefficients that minimize the mean-squared-error between the true fiducial positions and the positions estimated from images corrected using these coefficients. To examine the effect of higher-order and even-order terms, this calibration was performed using spherical harmonic polynomial of different orders up to the 10th-order including even- and odd-order terms, or odd-order only. The results showed that the model coefficients of this gradient can be successfully estimated. The residual root-mean-squared-error after correction using up to the 10th-order coefficients was reduced to 0.36 mm, yielding spatial accuracy comparable to conventional whole-body gradients. The even-order terms were necessary for accurate GNL modeling. In addition, the calibrated coefficients improved image geometric accuracy compared with the simulation-based coefficients.",
keywords = "asymmetric gradient, compact 3T, gradient nonlinearity, head-only MRI system, image geometric distortion",
author = "S. Tao and Trazasko, {Joshua D} and Gunter, {J. L.} and Weavers, {P. T.} and Yunhong Shu and Huston, {John III} and Lee, {S. K.} and Tan, {E. T.} and Bernstein, {Matthew A}",
year = "2017",
month = "1",
day = "21",
doi = "10.1088/1361-6560/aa524f",
language = "English (US)",
volume = "62",
pages = "N18--N31",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "2",

}

TY - JOUR

T1 - Gradient nonlinearity calibration and correction for a compact, asymmetric magnetic resonance imaging gradient system

AU - Tao, S.

AU - Trazasko, Joshua D

AU - Gunter, J. L.

AU - Weavers, P. T.

AU - Shu, Yunhong

AU - Huston, John III

AU - Lee, S. K.

AU - Tan, E. T.

AU - Bernstein, Matthew A

PY - 2017/1/21

Y1 - 2017/1/21

N2 - Due to engineering limitations, the spatial encoding gradient fields in conventional magnetic resonance imaging cannot be perfectly linear and always contain higher-order, nonlinear components. If ignored during image reconstruction, gradient nonlinearity (GNL) manifests as image geometric distortion. Given an estimate of the GNL field, this distortion can be corrected to a degree proportional to the accuracy of the field estimate. The GNL of a gradient system is typically characterized using a spherical harmonic polynomial model with model coefficients obtained from electromagnetic simulation. Conventional whole-body gradient systems are symmetric in design; typically, only odd-order terms up to the 5th-order are required for GNL modeling. Recently, a high-performance, asymmetric gradient system was developed, which exhibits more complex GNL that requires higher-order terms including both odd- and even-orders for accurate modeling. This work characterizes the GNL of this system using an iterative calibration method and a fiducial phantom used in ADNI (Alzheimer's Disease Neuroimaging Initiative). The phantom was scanned at different locations inside the 26 cm diameter-spherical-volume of this gradient, and the positions of fiducials in the phantom were estimated. An iterative calibration procedure was utilized to identify the model coefficients that minimize the mean-squared-error between the true fiducial positions and the positions estimated from images corrected using these coefficients. To examine the effect of higher-order and even-order terms, this calibration was performed using spherical harmonic polynomial of different orders up to the 10th-order including even- and odd-order terms, or odd-order only. The results showed that the model coefficients of this gradient can be successfully estimated. The residual root-mean-squared-error after correction using up to the 10th-order coefficients was reduced to 0.36 mm, yielding spatial accuracy comparable to conventional whole-body gradients. The even-order terms were necessary for accurate GNL modeling. In addition, the calibrated coefficients improved image geometric accuracy compared with the simulation-based coefficients.

AB - Due to engineering limitations, the spatial encoding gradient fields in conventional magnetic resonance imaging cannot be perfectly linear and always contain higher-order, nonlinear components. If ignored during image reconstruction, gradient nonlinearity (GNL) manifests as image geometric distortion. Given an estimate of the GNL field, this distortion can be corrected to a degree proportional to the accuracy of the field estimate. The GNL of a gradient system is typically characterized using a spherical harmonic polynomial model with model coefficients obtained from electromagnetic simulation. Conventional whole-body gradient systems are symmetric in design; typically, only odd-order terms up to the 5th-order are required for GNL modeling. Recently, a high-performance, asymmetric gradient system was developed, which exhibits more complex GNL that requires higher-order terms including both odd- and even-orders for accurate modeling. This work characterizes the GNL of this system using an iterative calibration method and a fiducial phantom used in ADNI (Alzheimer's Disease Neuroimaging Initiative). The phantom was scanned at different locations inside the 26 cm diameter-spherical-volume of this gradient, and the positions of fiducials in the phantom were estimated. An iterative calibration procedure was utilized to identify the model coefficients that minimize the mean-squared-error between the true fiducial positions and the positions estimated from images corrected using these coefficients. To examine the effect of higher-order and even-order terms, this calibration was performed using spherical harmonic polynomial of different orders up to the 10th-order including even- and odd-order terms, or odd-order only. The results showed that the model coefficients of this gradient can be successfully estimated. The residual root-mean-squared-error after correction using up to the 10th-order coefficients was reduced to 0.36 mm, yielding spatial accuracy comparable to conventional whole-body gradients. The even-order terms were necessary for accurate GNL modeling. In addition, the calibrated coefficients improved image geometric accuracy compared with the simulation-based coefficients.

KW - asymmetric gradient

KW - compact 3T

KW - gradient nonlinearity

KW - head-only MRI system

KW - image geometric distortion

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

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

U2 - 10.1088/1361-6560/aa524f

DO - 10.1088/1361-6560/aa524f

M3 - Article

C2 - 28033119

AN - SCOPUS:85010064965

VL - 62

SP - N18-N31

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 2

ER -