Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2933924

Formats

Article sections

Authors

Related links

Magn Reson Med. Author manuscript; available in PMC 2011 January 1.

Published in final edited form as:

PMCID: PMC2933924

NIHMSID: NIHMS228021

Julian Maclaren,^{1} Oliver Speck,^{2} Daniel Stucht,^{2} Peter Schulze,^{2} Jürgen Hennig,^{1} and Maxim Zaitsev^{1}

Corresponding author: J. Maclaren, Medical Physics, Department of Diagnostic Radiology, University Hospital Freiburg, Hugstetterstr. 55, 79106, Freiburg, Germany. Phone: +49 761 2709376, Fax: +49 761 2709379, Email: julian.maclaren/at/uniklinik-freiburg.de

The publisher's final edited version of this article is available at Magn Reson Med

See other articles in PMC that cite the published article.

Prospective motion correction in MR imaging is becoming increasingly popular to prevent the image artefacts that result from subject motion. Navigator information is used to update the position of the imaging volume before every spin excitation so that lines of acquired *k*-space data are consistent. Errors in the navigator information, however, result in residual errors in each *k*-space line. This paper presents an analysis linking noise in the tracking system to the power of the resulting image artefacts. An expression is formulated for the required navigator accuracy based on the properties of the imaged object and the desired resolution. Analytical results are compared with computer simulations and experimental data.

Subject motion remains a major problem in MRI. For high resolution brain imaging, measurement times of about 10 to 30 minutes are required to achieve an acceptable signal-to-noise ratio (SNR). During this time, involuntarily movements of the order of millimetres are unavoidable (1). The problem of motion is not limited to high-resolution imaging: in the case of fMRI, motion can produce false-positive activations (2); in routine clinical imaging, motion results in many scans having to be repeated.

Numerous motion-correction methods exist. For abdominal imaging, traditional navigator techniques can be used to selectively acquire data at only a particular point in the breathing cycle (3, 4). Rigid body head motion, however, requires a different approach. PROPELLER (5) has proved popular, although it has disadvantages, such as redundantly oversampling the *k*-space centre, and being sequence-specific and limited to in-plane motion correction. Other similar techniques (6, 7) suffer from similar drawbacks to PROPELLER, as do approaches that use post-processing correction, such as the autofocus method proposed by Atkinson et al. (8). Image-based navigation and correction methods, such as that proposed by Thesen et al. (9), can only correct for inter-scan motion.

Prospective motion correction is becoming a popular avenue of research, as it has advantages over the above approaches. Navigator data are used to update the scanner gradients, and therefore the position of the imaging volume, before every spin excitation. Although the term ‘prospective’ implies prediction, normally the latest estimate of the pose of the object is used to set the position and orientation of the imaging volume for the next excitation (1). Motion tracking of a rigid body in a full six degrees of freedom (6DOF) is possible and the technique can be applied to any existing pulse sequence to correct for head motion.

There are numerous examples of prospective correction being used during brain imaging in MRI. A variety of navigator systems have been used: reflecting laser beams off a marker attached to the subject (10, 11); optical tracking with multiple cameras (1, 12); optical tracking with a single camera (13, 14); or by tracking using small coils, known as ‘active markers’ (15, 16). The evaluation of these different systems is the subject of current research.

Regardless of the navigation system used, tracking accuracy is limited by noise in the navigator data. When imaging a stationary object, these pose errors generate artefacts in the image.

Despite the popularity of navigators, to our knowledge no analysis of the required navigator accuracy is available in literature (although accuracy has been quantified in the case of navigator echoes (17)). Hence, we develop a formalism to statistically analyse the image artefacts introduced by the prospective motion correction of each *k*-space line in the presence of tracking noise. It is then possible to predict artefact levels given knowledge of the tracking system parameters and the imaged object. The reverse is also true: one can calculate the required navigator accuracy of a tracking system given a typical target image, scanning parameters and image-quality prerequisites. This will assist the specification and development of such tracking systems.

The *k*-space signal intensity in traditional spin warp imaging, ignoring saturation, relaxation and dephasing effects can be given as

where *ρ* (**r**) is the signal density of the imaged object. Prospective motion correction using a noisy tracking system in the absence of motion adds a small uncertainty, **Δ** = **i**Δ* _{x}* +

(1)

where the approximation is valid only for **Δ** 1 voxel.

Applying the inverse Fourier transform to Eq. (1) to generate an image yields

(2)

where *ρ* (**r**) is the original signal density and *η* (**r**) describes the image-space artefacts that result from the navigator noise.

Eq. (2) can be further analysed statistically to derive the ‘artefact power’, defined here as |*η* (**r**)|^{2}, the expected value of the mean-squared artefact intensity in the image. Assuming zero-mean white uncorrelated noise with standard deviation *σ* and
and Δ* _{x}* (

(3)

where is the gradient operator.

Expanding Eq. (3) and noting that the artefact power is a function of *x*, but not *y* (the direction of phase-encoding), gives

(4)

Finally, incorporating the relative size of the FOV gives an expression for the expected value of the mean-squared artefact intensity in terms of the navigator noise in millimetres:

(5)

where *N _{y}* is the image matrix size in the phase-encoding direction and FOV

Eq. (5) shows that the mean artefact power is proportional to the integral edge power in the image and the navigator noise variance. Additionally, it is clear that the mean artefact power can be separated into the noise generated by the *x*- and *y*-components of navigator noise.

To verify the theoretical results presented here, we performed simulations and MR measurements. Both methods produce images corrupted with artefacts caused by navigator noise, which are then compared directly to the predictions made by Eq. (5).

We use the analytical version of the Shepp-Logan phantom for simulation purposes (Fig. 1(a)). The grey levels in the original phantom were chosen to represent the radiographic density variations in the human brain (18). In MRI, radiographic density is not the measured parameter; therefore, the original grey levels are inappropriate and have been adjusted here, as was done by Devaney (19).

The Shepp-Logan phantom in (a) the reference position and (b) after the introduction of a pose perturbation. The computed *k*-space data set (c) uses different pose parameters of the phantom for each *k*-space line, *i*, denoted by *x*_{i}, *y*_{i} and *θ*_{i}. As **...**

Motion simulations are performed by computing each line in *k*-space individually, where the computed values correspond to a Shepp-Logan phantom in a specified position. This position varies between *k*-space lines, according to parameters specified by a random number generator. The equation for the Fourier transform of a rotated, shifted ellipse (20) is applied to our situation. The *k*-space value at coordinates *k _{x}* and

(6)

where *J*_{1} is a first-order Bessel function of the first kind and *x _{i}*,

Phantom experiments were carried out on a 3T Magnetom TRIO system (Siemens Healthcare, Erlangen, Germany) using a single-channel birdcage coil. A standard 2D gradient echo sequence was modified to enable real-time motion position input. A stationary phantom was imaged with an in-plane resolution of 0.5 × 0.5 mm^{2} and slice thickness of 4 mm. Position noise was simulated using a random number generator with a white Gaussian distribution (21). The advantage of using this method, rather than using noise data from an actual tracking system, is that the navigator noise parameters can be tightly specified. Images were reconstructed with no filtering or other post-processing.

An *in vivo* experiment was performed on a 7T high-field system (Siemens Medical Solutions, Erlangen, Germany). A stereoscopic tracking system (ARTtrack3, Advanced Realtime Tracking GmbH, Weilheim, Germany) was used for prospective motion correction in conjunction with a gradient echo sequence modified for the purpose. Imaging resolution was 0.5 × 0.5 mm with a slice thickness of 4 mm and an axial slice orientation. The field-of-view was deliberately chosen to be much larger than the object, so that any artefacts could be clearly seen outside the support region. The navigator error standard deviations at the centre of the tracking target were *σ _{x}* = 0.08 mm and

Given an *N _{x}* ×

(7)

where the vertical bars indicate the complex magnitude. It is also useful to define a signal-to-artefact ratio, given by

(8)

Eqs (7) and (8) describe the computation of artefact power for measured data. This can now be compared to the predicted values from Eq. (5), where the theoretical value for *AP* (*x*) is equivalent to |*η* (*x*)|^{2}. The evaluation of Eq. (5) requires the computation of the partial derivatives in the image,
and
. This can conveniently be performed by frequency multiplication in the Fourier domain, which avoids errors introduced by a finite difference derivative calculation.

Fig. 2(a) and (b) show images reconstructed from simulations where navigator noise is applied in the *x*- and *y*-directions, respectively. The navigator noise standard deviation, *σ _{x}* and

Simulations with navigator noise of *σ* = 0.1 pixels in (a) the *x*-direction and (b) the *y*-direction. After subtraction of a reference image, the resulting artefacts from (a) and (b) are seen more clearly in (c) and (d), respectively. Over 100 iterations, **...**

To explore the relationship of the navigator noise variance to total artefact power, denoted here Σ*AP* (*x*) (the denominator in Eq. (8)), the above simulation was carried out for a range of values for *σ _{x}* (Fig. 3). A line of constant slope
is plotted for comparison. As predicted by Eq. (5), artefact power increases linearly with the navigator noise variance for small values of
. When

Simulations were also performed for rotation noise with *σ _{θ}* = 0.2°, and the resulting artefacts and artefact power are shown in Fig. 4. Although the extension of the theory to consider rotations is left to future work, it is reasonable to expect that edges lying in a radial direction would be a source of artefacts in this scenario.

(a) Artefacts resulting from a rotation simulation using the Shepp-Logan phantom with *σ*_{θ} = 0.2° and (b) the artefact power generated by summing down the columns in (a).

The phantom imaging experiments produce similar results to the simulations. Fig. 5(a) shows an artefact-free image. Fig. 5(b)–(d) show sample results with 0.5 × 0.5 mm in-plane spatial resolution and navigator noise with standard deviation of *σ* = 0.5 mm in (b) both directions, (c) the *x*-direction only and (d) the *y*-direction only. Again, vertical edges in the image generate artefacts when navigator noise in the *x*-direction is present, while horizontal edges in the image generate artefacts when navigator noise in the *y*-direction is present.

Phantom images acquired (a) without navigator noise, (b) with isotropic noise of *σ*_{x}, *σ*_{y} = 0.5 mm, (c) with the same noise in only the *x*- and (d) *y*-directions. It is apparent that vertical edges generate artefacts in figure (c) and horizontal **...**

Fig. 6 shows a more quantitative comparison between theory and measurements. The calculated mean artefact power based on the experimental results is plotted against that predicted by the theory based on Eq. (5). Fig. 6(a) and (b) show simulations results; Fig. 6(c) and (d) show results from the imaged phantom. For the simulations in (a) and (b), the total measured signal-to-artefact ratio is almost identical to the predicted value. For the MR data, there is a slight discrepancy: for navigator noise applied in the *x*-direction, signal-to-artefact ratio is 9.52 × 10^{5} compared to the predicted value of 8.97 × 10^{5}; for noise applied in the *y*-direction, signal-to-artefact ratio is 9.77 × 10^{5} compared to the predicted value of 9.29 × 10^{5}. This discrepancy can be explained by system imperfections increasing the difference between the images.

Mean artefact power generated by the *x*-component of the navigator noise for (a) simulated data (*σ*_{x} = 0.1 pixels) and (c) MR data (*σ*_{x} = 0.25 mm). Mean artefact power generated by the *y*-component of the navigator noise for (b) simulated **...**

It is now possible to go beyond artefact prediction to specifying the tracking accuracy requirements for prospective motion correction. Assuming that the navigator noise in the *x*- and *y*-directions has the same variance, then we define
. From Eq. (5),

(9)

The procedure to compute the required image accuracy is then

- Select an acceptable artefact level (i.e., a value for signal-to-artefact ratio). This can be done by taking a fraction of the native noise level, such that signal-to-artefact ratio > SNR.
- Compute Σ
*AP*, based on the desired signal-to-artefact ratio and Eq. (8). - Compute the edge power of the object in both the
*x*- and*y*-directions. - Evaluate Eq. (9) using the above information and the FOV size and imaging resolution.

As an example, the above is applied to the brain image in Fig. 7(a). For practical applications it is desirable for the signal-to-artefact ratio to be significantly better than the native image SNR. Here we decide that the artefact level caused by navigator noise is to be 1/4 that of the native image noise, so that signal-to-artefact ratio = 4 SNR. A typical high-resolution image, with 0.5 × 0.5 mm resolution, has an SNR of about 25. Thus, we set signal-to-artefact ratio = 100. Given this, and the edge energy in our target image, Σ*AP* = 9.0 × 10^{−6}.

(a) A high resolution brain image obtained on a 7T system using prospective motion correction with tracking noise *σ*_{x} = 0.08 mm, *σ*_{y} = 0.05 mm. (b) The same image raised to the power of 0.2, to enhance artefacts. (c) The predicted *in vivo* **...**

Finally, using the FOV size and imaging resolution, we evaluate Eq. (9) and obtain *σ _{xy}* = 0.18 mm. Thus, to meet our target value for signal-to-artefact ratio, tracking noise standard deviation must be 0.18 mm or better in both

The navigator noise constraints calculated here are well above the measured parameters from the ARTtrack3 tracking system of *σ _{x}* = 0.08 mm and

We have shown that noise in the tracking data used for prospective motion correction always results in artefacts, which can be quantified in terms of a ‘signal-to-artefact’ ratio. Mean artefact power is proportional to both the integral of the square of the derivatives of the image intensity and to the variance of navigator noise. For small values of navigator noise variance, artefacts are not always visible in the reconstructed image (as is the case in Fig. 2(a) and (b)). Artefact power increases linearly, however, with navigator noise variance. Several tracking systems mentioned in the literature (such as (22), for example) report a level of navigator noise that means artefacts will certainly be visible in the image.

By inverting the derived expression for artefact power, the required navigator accuracy can be computed. This is a function of the object properties and the acceptable level of residual artefacts. In general, however, navigator accuracy should be several times better than the target pixel size for most imaging situations.

Although these requirements seem stringent, some image quality may be recovered in postprocessing, which would relax requirements on the tracking data accuracy to achieve the same final image quality. Techniques for post-processing motion correction, such as that proposed by Bammer et al. (23) or Atkinson and Hill (24), could be applied for reduction of artefacts caused by tracking noise. However, this is only applicable if tracking noise can be estimated after data acquisition. Similarly, approaches such as Kalman filtering, recently used for cardiac navigators (25), are being investigated for head motion tracking (26). This would also relax the accuracy requirements of the tracking system itself.

This analysis only accounts for the effects of in-plane translation errors. Omitting errors in the remaining degrees of freedom is justified if the effect of these errors on artefact power is much smaller than in-plane translation errors. This applies in the case of the ARTtrack3 tracking system used in this study. Provided the target is large relative to the location errors of individual markers, the rotation of the target is accurately estimated. Our simulations show that acquisitions with realistic values for rotation noise do not display visible artefacts. The omission of through-plane translation is justified for 2D imaging, as the slice thickness used is normally around a factor of four greater than the in-plane resolution and is much smaller than navigator noise (Fig. 8). Additionally, the non-ideal slice profile results in a ‘smoothing effect’ if through-plane navigator noise moves the excited plane over a tissue boundary.

Effects of through-plane navigator noise can be neglected due to the relatively poor through-plane resolution and the non-ideal slice profile, typical in 2D Fourier imaging. (a) A sharp tissue boundary representing the ‘worst case’ in **...**

The theory presented here makes several assumptions. We assume that such navigator noise is normally distributed: visually, this appears to be the case for noise from the ARTtrack3 tracking system and a Jarque-Bera test (27) applied to the navigator noise data for each DOF was unable to reject the null-hypothesis that the data came from a normal distribution. However, the assumption of Gaussianity is a limitation to this work, as it may not apply to some navigation methods. Additionally, we assume that the imaging properties are independent of the position and orientation of the object. This will not be the case in prospective motion correction with large head movements, due to field distortions and non-uniform coil sensitivity profiles.

Another limitation of the approach described here is that mean noise power does not correlate exactly with human perception of artefacts. If artefacts are evenly distributed in a ‘noise-like’ fashion, then they will be less visible than if they appear concentrated in one area, as occurs in the case of ghosting. Therefore, this work could be usefully extended through the incorporation of a ‘human perception’ factor. This would likely incorporate the maximum value of generated artefacts, rather than their sum over the image.

Here, we compute the image gradient in the *x*- and *y*-directions using frequency multiplication in *k*-space. The result is not identical to the gradient of the imaged object, due to the truncation of *k*-space. However, we show that it is the edge power in the image *after* sampling that generates image artefacts, rather than sharp edges present in the object that are lost during imaging. Due to Fourier encoding, the acquired data in MRI is inherently band-limited; thus, errors in frequency components of the object that are not present in the acquired signal, have no effect on the final image. This concept has been recognised in the application of super-resolution to Fourier-encoded MRI data (28) where unsampled spatial components cannot be retrieved using conventional super-resolution techniques involving shifts in the FOV.

It is worth noting that the theory and results presented here can be applied in a more general sense to motion artefacts in MRI. The situation presented is identical to a line-by-line acquisition of *k*-space, without prospective correction but with Gaussian-distributed position noise caused by patient motion. Although patient motion parameters will not necessarily be Gaussian distributed, a similar expression for artefact generation may apply. It is reasonable to expect that image edges are crucial to artefact generation.

This work derives the effect of navigator noise on image artefact power. Theoretical results closely match simulations and MR data. For small translation errors, artefact power is directly proportional to the variance of the navigator noise in the tracking system. The complete expression for mean artefact power enables the specification of the required accuracy for tracking systems in prospective motion correction. Knowledge of this relationship will aid in the design of tracking systems for prospective motion correction. Tracking system accuracy must be a fraction of a pixel to be effective and to avoid the introduction of extra artefacts when imaging a stationary object.

This work is a part of the INUMAC project supported by the German Federal Ministry of Education and Research (grant 01EQ0605) and the NIH/NIDA project ‘Real-time motion tracking and correction for MR Imaging and Spectroscopy’ (grant 1R01 DA021146). We thank Dr. Valerij Kiselev for helpful discussions.

1. Zaitsev M, Dold C, Sakas G, Hennig J, Speck O. Magnetic resonance imaging of freely moving objects: prospective real-time motion correction using an external optical motion tracking system. Neuroimage. 2006;31(3):1038–1050. [PubMed]

2. Hajnal JV, Myers R, Oatridge A, Schwieso JE, Young IR, Bydder GM. Artifacts due to stimulus correlated motion in functional imaging of the brain. Magn Reson Med. 1994;31(3):283–91. [PubMed]

3. Ehman RL, Felmlee JP. Adaptive technique for high-definition MR imaging of moving structures. Radiology. 1989;173(1):255–263. [PubMed]

4. Wang Y, Rossman PJ, Grimm RC, Riederer SJ, Ehman RL. Navigator-echo-based real-time respiratory gating and triggering for reduction of respiration effects in three-dimensional coronary MR angiography. Radiology. 1996;198(1):55–60. [PubMed]

5. Pipe JG. Motion correction with PROPELLER MRI: application to head motion and free-breathing cardiac imaging. Magn Reson Med. 1999;42(5):963–969. [PubMed]

6. Welch EB, Rossman PJ, Felmlee JP, Manduca A. Self-navigated motion correction using moments of spatial projections in radial MRI. Magn Reson Med. 2004;52(2):337–345. [PubMed]

7. Maclaren JR, Bones PJ, Millane RP, Watts R. MRI with TRELLIS: a novel approach to motion correction. Magn Reson Imaging. 2008;26(4):474–83. [PubMed]

8. Atkinson D, Hill DLG, Stoyle PNR, Summers PE, Clare S, Bowtell R, Keevil SF. Automatic compensation of motion artifacts in MRI. Magn Reson Med. 1999;41(1):163–170. [PubMed]

9. Thesen S, Heid O, Mueller E, Schad LR. Prospective acquisition correction for head motion with image-based tracking for real-time fMRI. Magn Reson Med. 2000;44(3):457–463. [PubMed]

10. Eviatar H, Saunders J, Hoult D. Motion compensation by gradient adjustment. Proceedings of the International Society for Magnetic Resonance in Medicine; Vancouver. 1997. p. 1898.

11. Eviatar H, Schattka B, Sharp J, Rendell J, Alexander M. Real time head motion correction for functional MRI. Proceedings of the International Society for Magnetic Resonance in Medicine; Philadelphia. 1999. p. 269.

12. Speck O, Hennig J, Zaitsev M. Prospective real-time slice-by-slice motion correction for fMRI in freely moving subjects. Magma. 2006;19(2):55–61. [PubMed]

13. Aksoy M, Newbould R, Straka M, Holdsworth S, Skare S, Santos J, Bammer R. A Real Time Optical Motion Correction System Using a Single Camera and 2D Marker. Proceedings 16th Scientific Meeting, International Society for Magnetic Resonance in Medicine; 2008. p. 3120.

14. Qin L, Gelderen P, Zwart J, Jin F, Tao Y, Duyn J. Head Movement Correction for MRI With a Single Camera. Proceedings 16th Scientific Meeting, International Society for Magnetic Resonance in Medicine; 2008. p. 1467.

15. Derbyshire JA, Wright GA, Henkelman RM, Hinks RS. Dynamic scan-plane tracking using MR position monitoring. Journal of Magnetic Resonance Imaging. 1998;8(4):924–932. [PubMed]

16. Krueger S, Schaeffter T, Weiss S, Nehrke K, Rozijn T, Boernert P. Prospective Intra-Image Compensation for Non-Periodic Rigid Body Motion Using Active Markers. Proceedings 14th Scientific Meeting, International Society for Magnetic Resonance in Medicine; 2006. p. 3196.

17. Felmlee JP, Ehman RL, Riederer SJ, Korin HW. Adaptive motion compensation in MRI: accuracy of motion measurement. Magn Reson Med. 1991;18(1):207–13. [PubMed]

18. Shepp LA, Logan BF. Fourier reconstruction of a head section. IEEE Transactions on Nuclear Science. 1974;NS21(3):21–43.

19. Devaney AJ. A computer simulation study of diffraction tomography. IEEE Trans Biomed Eng. 1983;30(7):377–86. [PubMed]

20. Pan S, Kak A. A computational study of reconstruction algorithms for diffraction tomography: interpolation versus filtered-backpropagation. IEEE Transactions on Acoustics, Speech, and Signal Processing. 1983;31(5):1262–1275.

21. Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in C. Cambridge University Press; Cambridge, UK: 1992.

22. Nevo E, Hushek S. An Electromagnetic 3D Locator System For Use In MR Scanners. Proceedings of the International Society for Magnetic Resonance in Medicine; Hawaii. 2002. p. 334.

23. Bammer R, Aksoy M, Liu C. Augmented generalized SENSE reconstruction to correct for rigid body motion. Magnetic Resonance in Medicine. 2007;57(1):90–102. [PubMed]

24. Atkinson D, Hill DLG. Reconstruction after rotational motion. Magn Reson Med. 2003;49:183–187. [PubMed]

25. Spincemaille P, Nguyen TD, Prince MR, Wang Y. Kalman filtering for real-time navigator processing. Magn Reson Med. 2008;60(1):158–168. [PubMed]

26. Maclaren J, Speck O, Hennig J, Zaitsev M. A Kalman filtering framework for prospective motion correction. Proceedings 17th Scientific Meeting, International Society for Magnetic Resonance in Medicine; Honolulu, Hawaii. 2009. p. 4602.

27. Jarque CM, Bera AK. A test for normality of observations and regression residuals. International Statistical Review. 1987;55(2):163–172.

28. Peled S, Yeshurun Y. Superresolution in MRI–perhaps sometimes. Magnetic Resonance in Medicine. 2002;48(2):409.

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |