Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter August 13, 2013

Twenty-fold acceleration of 3D projection reconstruction MPI

Justin J. Konkle, Patrick W. Goodwill, Emine Ulku Saritas, Bo Zheng, Kuan Lu and Steven M. Conolly

Abstract

We experimentally demonstrate a 20-fold improvement in acquisition time in projection reconstruction (PR) magnetic particle imaging (MPI) relative to the state-of-the-art PR MPI imaging results. We achieve this acceleration in our imaging system by introducing an additional Helmholtz electromagnet pair, which creates a slow shift (focus) field. Because of magnetostimulation limits in humans, we show that scan time with three-dimensional (3D) PR MPI is theoretically within the same order of magnitude as 3D MPI with a field free point; however, PR MPI has an order of magnitude signal-to-noise ratio gain.


Corresponding author: Justin J. Konkle, Department of Bioengineering, University of California, Berkeley, 340 Hearst Memorial Mining Building, Berkeley, CA 94720-1762, USA, E-mail:

J.J. Konkle and P.W. Goodwill contributed equally to this work. This work was supported in part by the CIRM Tools and Technology Grant RT2-01893, the University of California Discovery Grant, the NIH National Institute of Biomedical Imaging and Bioengineering grant 1R01EB013689, and the National Science Foundation Graduate Research Fellowship grant DGE 1106400. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of CIRM or any other agency of the State of California, the National Institute of Biomedical Imaging and Bioengineering, or the National Institutes of Health.

Appendices

Appendix 1 Projection reconstruction SNR gain calculation

We calculate the SNR gain of 3D PR with an FFL as compared to 3D imaging with an FFP with equal imaging time. Here, we assume that the two systems have identical noise characteristics, pulse sequence, resolution, and gradient strength. In a simulation study, Weizenecker et al. [40] noted that the SNR of a 3D PR image is approximately proportional to

where N is the number of projections acquired. Here, we show how the SNR is affected by the FBP operation. We begin by calculating the noise variance of FBP:

where σ0 is the standard deviation of noise per pixel assuming Gaussian random noise, k is spatial frequency, G(k) is the filter frequency response, φ is the angle in radians, and Kmax is the maximum radial spatial frequency of the acquired image (i.e., the total extent is 2Kmax). Similar calculations have been done in CT [33]. Note that filtering takes place after discretization of the acquired projection images. Hence, the spatial-frequency domain is “normalized” such that the maximum frequency is one-half cycle per pixel, i.e., 2Kmax=1. For a ramp filter G(k)=k, the noise variance becomes

We then calculate the SNR for a PR image after FBP with a ramp filter:

where μ is the signal mean or expected value. Similarly, the noise variance in an image acquired with an FFP can be calculated:

Again, we used 2Kmax=1. Thus,

Finally, we find

In general, the exact value of the constant multiplier before

depends on the shape of the filter G(k).

Appendix 2 MPI acquisition time calculation

We can calculate the optimum imaging time based on specific absorption rate (SAR) and magnetostimulation (dB/dt) limits, the two primary safety concerns when imaging human subjects using time-varying magnetic fields. Owing to the frequency range in which MPI operates, magnetostimulation (and not SAR) is the dominant limitation for scanning speed in MPI [11, 35]. For human-sized MPI scanners, magnetostimulation will restrict the amplitude of both the excitation (drive) field and the slow shift (focus) fields.

The drive field in MPI is typically a ~25 kHz frequency sinusoidal field. Between 5 and 50 kHz, the magnetostimulation threshold in the human torso is extrapolated as approximately Bth=7 mT [35]. With a G=10 T/m gradient strength, which will produce 1 mm native (i.e., no deconvolution) resolution with Resovist [11–13, 16], the FOV can be calculated as FOV=2Bth/G=1.4 mm.

To address the limited FOV coverage of the drive subsystem, slow shift magnets [12] or focus field magnets [37] are used. These slow shift magnets slowly raster the mean position of the FFP or FFL (see Figure 8), expanding the FOV beyond what is covered by the drive field alone. In a system with slow shift magnets, the space covered solely by the drive field (with the slow shift field disabled) is termed a “partial FOV” (pFOV) [12] or “imaging station” with multistation reconstruction [37]. Slow shift fields also limit the acquisition time for an MPI system owing to magnetostimulation limits. The ICNIRP define a maximal magnetic field slew rate of 20 T/s for pulse durations longer than a couple of milliseconds [18]. As we determine below, the slew rate causes the slow shift fields to become binding constraints on imaging time in addition to the drive field.

Figure 8 Trajectory distance calculation for the Lissajous pattern.The Lissajous pattern is created with two sinusoidal drive fields and is slowly rastered through the field of view with two slow shift (focus) fields.

Figure 8

Trajectory distance calculation for the Lissajous pattern.

The Lissajous pattern is created with two sinusoidal drive fields and is slowly rastered through the field of view with two slow shift (focus) fields.

Generally, the total imaging time, T [s], in MPI can be calculated as

where D [m] is the total distance traveled by the field free region and vs [m/s] is the slow shift field scanning rate. The total distance traveled depends on the type of drive field and slow shift field pulse sequences. Commonly, a linear or Lissajous trajectory is used for the drive field. These trajectories provide near ideal spatial coverage and are easy to calculate [32, 39].

The total distance traveled by the mean position of the field free region dominates the total imaging time assuming that (1) the time to cover the pFOV is nearly instantaneous, (2) magnetostimulation thresholds from the drive field and the slow shift fields do not affect each other, and (3) the pFOV is small relative to the total FOV (see Figure 8). The total distance is then calculated as

where

Here, 2D refers to a single slice in FFP imaging or a single projection using an FFL, and 3D refers to imaging using an FFP. FOVi is the size of the field of view along axis i, pFOVi is the size of the partial FOV along i, pi is the number of pFOVs along i, and β is a factor (>1) determining the overlap extent of the pFOVs required for baseline recovery [29]. We can calculate the maximal size of the pFOV using the simple relation pFOV=2Bth/G [m]. Accordingly, the imaging time can be estimated as

where Smax=20 T/s is the maximum slew rate as described above and μ0Gab [T/m] is the partial derivative of the magnetic field in the a direction with respect to b.

References

[1] Achenbach S, Marwan M, Ropers D, et al. Coronary computed tomography angiography with a consistent dose below 1 mSv using prospectively electrocardiogram-triggered high-pitch spiral acquisition. Eur Heart J 2010; 31: 340–346.10.1093/eurheartj/ehp470Search in Google Scholar PubMed

[2] Barrett T, Brechbiel M, Bernardo M, Choyke PL. MRI of tumor angiogenesis. J Magn Reson Imaging 2007; 26: 235–249.10.1002/jmri.20991Search in Google Scholar PubMed

[3] Bulte JWM, Walczak P, Gleich B, et al. MPI cell tracking: what can we learn from MRI? Imaging 2011; 7965: 79650Z–79650Z-4.10.1117/12.879844Search in Google Scholar PubMed PubMed Central

[4] Chandarana H, Godoy MCB, Vlahos I, et al. Abdominal aorta: evaluation with dual-source dual-energy multidetector CT after endovascular repair of aneurysms – initial observations. Radiology 2008; 249: 692–700.10.1148/radiol.2492080359Search in Google Scholar PubMed

[5] Erbe M, Knopp T, Sattel TF, Biederer S, Buzug TM. Experimental generation of an arbitrarily rotated field-free line for the use in magnetic particle imaging. Med Phys 2011; 38: 5200–5207.10.1118/1.3626481Search in Google Scholar PubMed

[6] Ferguson RM, Khandhar AP, Krishnan KM. Tracer design for magnetic particle imaging (invited). J Appl Phys 2012; 111: 7B318–7B3185.10.1063/1.3676053Search in Google Scholar PubMed PubMed Central

[7] Ferguson RM, Minard KR, Khandhar AP, Krishnan KM. Optimizing magnetite nanoparticles for mass sensitivity in magnetic particle imaging. Med Phys 2011; 38: 1619.10.1118/1.3554646Search in Google Scholar PubMed PubMed Central

[8] Fleischmann D, Boas FE. Computed tomography – old ideas and new technology. Eur Radiol 2011; 21: 510–517.10.1007/s00330-011-2056-zSearch in Google Scholar PubMed

[9] Gleich B, Weizenecker J. Tomographic imaging using the nonlinear response of magnetic particles. Nature 2005; 435: 1214–1217.10.1038/nature03808Search in Google Scholar PubMed

[10] Goldman LW. Principles of CT and CT technology. J Nucl Med Technol 2007; 35: 115–28.10.2967/jnmt.107.042978Search in Google Scholar PubMed

[11] Goodwill PW, Conolly SM. The x-space formulation of the magnetic particle imaging process: 1-D signal, resolution, bandwidth, SNR, SAR, and magnetostimulation. IEEE Trans Med Imaging 2010; 29: 1851–1859.10.1109/TMI.2010.2052284Search in Google Scholar

[12] Goodwill PW, Conolly SM. Multidimensional x-space magnetic particle imaging. IEEE Trans Med Imaging 2011; 30: 1581–1590.10.1109/TMI.2011.2125982Search in Google Scholar

[13] Goodwill PW, Conolly SM. Experimental demonstration of x-space magnetic particle imaging. In: Proceedings of SPIE, vol. 7965, 2011: 79650U–79650U-6.10.1117/12.878659Search in Google Scholar

[14] Goodwill PW, Konkle JJ, Zheng B, Saritas EU, Conolly SM. Projection x-space magnetic particle imaging. IEEE Trans Med Imag 2012; 31: 1076–1085.10.1109/TMI.2012.2185247Search in Google Scholar

[15] Goodwill PW, Lu K, Zheng B, Conolly SM. An x-space magnetic particle imaging scanner. Rev Sci Instrum 2012; 83: 033708.10.1063/1.3694534Search in Google Scholar

[16] Goodwill PW, Tamrazian A, Croft LR, et al. Ferrohydrodynamic relaxometry for magnetic particle imaging. Appl Phys Lett 2011; 98: 262502.10.1063/1.3604009Search in Google Scholar

[17] Holdsworth DW, Thornton MM. Micro-CT in small animal and specimen imaging. Trends Biotechnol 2002; 20: S34–S39.10.1016/S0167-7799(02)02004-8Search in Google Scholar

[18] International Commission on Non-Ionizing Radiation Protection (ICNIRP). Medical magnetic resonance (MR) procedures: protection of patients. Health Phys 2004; 87: 197–216.10.1097/00004032-200408000-00008Search in Google Scholar PubMed

[19] Jorgensen SM, Demirkaya O, Ritman EL. Three-dimensional imaging of vasculature and parenchyma in intact rodent organs with X-ray micro-CT. Am J Physiol 1998; 275: H1103–H1114.10.1152/ajpheart.1998.275.3.H1103Search in Google Scholar PubMed

[20] Knopp T, Biederer S, Sattel TF, Erbe M, Buzug TM. Prediction of the spatial resolution of magnetic particle imaging using the modulation transfer function of the imaging process. IEEE Trans Med Imaging 2011; 30: 1284–1292.10.1109/TMI.2011.2113188Search in Google Scholar PubMed

[21] Knopp T, Erbe M, Biederer S, Sattel TF, Buzug TM. Efficient generation of a magnetic field-free line. Med Phys 2010; 37: 3538.10.1118/1.3447726Search in Google Scholar PubMed

[22] Knopp T, Erbe M, Sattel TF, Biederer S, Buzug TM. Generation of a static magnetic field-free line using two Maxwell coil pairs. Appl Phys Lett 2010; 97: 092505.10.1063/1.3486118Search in Google Scholar

[23] Knopp T, Erbe M, Sattel TF, Biederer S, Buzug TM. A Fourier slice theorem for magnetic particle imaging using a field-free line. Inverse Probl 2011; 27: 095004.10.1088/0266-5611/27/9/095004Search in Google Scholar

[24] Konkle J, Goodwill P, Carrasco-Zevallos O, Conolly S. Experimental 3D X-space magnetic particle imaging using projection reconstruction. In: Buzug TM, Borgert J, editors. Magnetic particle imaging. Lubeck, Germany: Springer Proceedings in Physics 2012: 243–247.10.1007/978-3-642-24133-8_39Search in Google Scholar

[25] Konkle JJ, Goodwill PW, Carrasco-Zevallos OM, Conolly SM. Projection reconstruction magnetic particle imaging. IEEE Trans Med Imaging 2013; 32: 338–347.10.1109/TMI.2012.2227121Search in Google Scholar PubMed PubMed Central

[26] Konkle J, Goodwill P, Conolly S. Development of a field free line magnet for projection MPI. In: Pelc NJ, Dawant BM, Summers RM, et al., editors. Proc. SPIE 7965, medical imaging 2011: biomedical applications in molecular, structural, and functional imaging, vol. 7965, 2011: 79650X–79650X-7.Search in Google Scholar

[27] Krauss B, Schmidt B, Flohr T. Dual source CT. In: Johnson T, Fink C, SchÖnberg SO, Reiser MF, editors. Dual energy CT in Clinical Practice. Berlin Heidelberg: Springer 2011: 11–20.10.1007/174_2010_44Search in Google Scholar

[28] Lu K, Goodwill P, Saritas E, Zheng B, Conolly S. Linearity and shift-invariance for quantitative magnetic particle imaging. IEEE Trans Med Imaging 2013; DOI: 10.1109/TMI.2013.2257177.10.1109/TMI.2013.2257177Search in Google Scholar PubMed PubMed Central

[29] Lu K, Goodwill P, Zheng B, Conolly S. The impact of filtering direct-feedthrough on the x-space theory of magnetic particle imaging. In: SPIE proceedings, vol. 7965, 2011: 79652I–79652I-7.10.1117/12.878446Search in Google Scholar

[30] Macovski A. Medical imaging systems, 1st ed. Kailath T, editor. Englewood Cliffs, NJ: Prentice Hall 1983.Search in Google Scholar

[31] Nishimura DG. Principles of magnetic resonance imaging. Stanford, CA: Stanford University 2010.Search in Google Scholar

[32] Rahmer J, Weizenecker J, Gleich B, Borgert J. Signal encoding in magnetic particle imaging: properties of the system function. BMC Med Imaging 2009; 9: 4.10.1186/1471-2342-9-4Search in Google Scholar PubMed PubMed Central

[33] Riederer SJ, Pelc NJ, Chesler DA. The noise power spectrum in computed X-ray tomography. Phys Med Biol 1978; 23: 446–454.10.1088/0031-9155/23/3/008Search in Google Scholar PubMed

[34] Saritas EU, Goodwill PW, Croft LR, et al. Magnetic particle imaging (MPI) for NMR and MRI researchers. J Magn Reson 2012; 229: 116–126.10.1016/j.jmr.2012.11.029Search in Google Scholar PubMed PubMed Central

[35] Saritas EU, Goodwill PW, Zhang GZ, Conolly SM. Magnetostimulation limits in magnetic particle imaging. IEEE Trans Med Imaging 2013; DOI: 10.1109/TMI.2013.2260764.10.1109/TMI.2013.2260764Search in Google Scholar PubMed

[36] Sattel TF, Knopp T, Biederer S. Single-sided device for magnetic particle imaging. J Phys D Appl Phys 2009; 42: 022001.10.1088/0022-3727/42/2/022001Search in Google Scholar

[37] Schmale I, Rahmer J, Gleich B, et al. First phantom and in vivo MPI images with an extended field of view. Imaging 2011; 7965: 796510–796510-6.10.1117/12.877339Search in Google Scholar

[38] Schmitz SA, Albrecht T, Wolf KJ. MR angiography with superparamagnetic iron oxide: feasibility study. Radiology 1999; 213: 603–607.10.1148/radiology.213.2.r99oc24603Search in Google Scholar PubMed

[39] Weizenecker J, Borgert J, Gleich B. A simulation study on the resolution and sensitivity of magnetic particle imaging. Phys Med Biol 2007; 52: 6363–6374.10.1088/0031-9155/52/21/001Search in Google Scholar PubMed

[40] Weizenecker J, Gleich B, Borgert J. Magnetic particle imaging using a field free line. J Phys D Appl Phys 2008; 41: 105009.10.1088/0022-3727/41/10/105009Search in Google Scholar

[41] Weizenecker J, Gleich B, Rahmer J, Dahnke H, Borgert J. Three-dimensional real-time in vivo magnetic particle imaging. Phys Med Biol 2009; 54: L1–L10.10.1088/0031-9155/54/5/L01Search in Google Scholar PubMed

Received: 2013-1-10
Accepted: 2013-7-8
Published Online: 2013-08-13
Published in Print: 2013-12-01

©2013 by Walter de Gruyter Berlin Boston

Scroll Up Arrow