Magnetic resonance imaging (MRI) is one of the most important imaging modalities in clinical routine today. Especially, the superior tissue contrast and the absence of ionizing radiation make MRI a vital tool for diagnosis. These upsides come with a longer time for image acquisition compared to computed tomography, which makes MRI particularly susceptible to motion artifacts . In particular in thorax imaging, there is body motion due to respiration and cardiac motion. Breath holding is one way to avoid respiratory motion artifacts  but it is inconvenient and not suitable for all patients and imaging schemes. This and cardiac motion, which is inevitable, make it necessary to deal with the problem of motion artifacts on the level of data acquisition or data processing. A way to do so are prospective and retrospective gating. In both cases, additional information to the complex k-space data is needed.
It is desirable that motion corruption could be estimated just with the complex MRI raw data without additional measurements or trigger signals. One approach for such a blind motion compensation algorithm is presented in . To test this blind motion compensation algorithm on real reproducible measurement data, a motion phantom needs to be utilized. Even though, there are commercially available MRI motion phantoms, none of them can be used in a small animal MRI system.
In this work, an MRI phantom for rigid motion is presented, which can be used in a small animal MRI system. Furthermore, results of motion corrupted imaging and retrospective gating are displayed.
2 Material and methods
For all experiments, an 1 T ICON TM MRI system (Bruker, Ettlingen, Germany) was used. It is designed for imaging small rodents and is based on a permanent magnet. The gradient system provides gradients up to 360–390 mT/m. All data were acquired with a mouse body transmit-receive coil.
2.1 Imaging sequence
For all images, data were acquired using a fast low angle shot sequence (FLASH) . For a good signal-to-noise ratio (SNR) a strong signal with a high echo amplitude is needed. The echo amplitude for basic gradient echo sequences like FLASH can be described as
with Mz0 the longitudinal magnetization in thermal equilibrium, TR the repetition time, TE the echo time, α the flip angle, T1 the longitudinal relaxation time, and the transverse relaxation time as shown in .
In this work, repetition times of 50 ms, echo times of 6 ms and flip angles of 15° were used.
The direction of object motion was chosen as phase-encoding direction. For periodic motion, this may cause motion artifacts if
with Tacq the acquisition time for all phase-encoding steps, the mth harmonic of the periodic movement and f0 the frequency of the object motion . For the first harmonic and the chosen 0.43 Hz motion frequency, (2) limits the acquisition time to 1.16 s if motion artifacts are to be avoided. With 3.2 s the acquisition time was chosen in a way that motion artifacts will appear.
2.2 Phase drift correction
In case of retrospective gating, one has to deal with longer times for imaging, since multiple repetitions are needed. For this work, a data acquisition time of 10 min for 200 repetitions was needed. In these dimensions, phase-drifts become relevant. Over time the total transverse magnetization drifts in phase but stays constant in magnitude for the same time after excitation. The oscillating signal here in complex notation can be described as
with S0 being the magnitude of the signal, ωrf the precession frequency and ϕ(t) the initial phase. For short imaging times, ϕ(t) can be assumed as constant. However, over longer times, ϕ(t) is not constant and results in a translation in image space. For combining data which is spread over several minutes in time, as needed for retrospective gating, a phase correction has to be performed as a phase subtraction so that
with ϕin the phase for the initial repetition, ϕac(t) the actual phase for a chosen repetition, and ϕco(t) the phase-correction value. Therefore, the sampled data was corrected by multiplication with a time dependent term
The repetition-dependent ϕco(t) is determined by fitting the phase for all phase-encoding steps with the same number and the same motor position with a 4th degree polynomial (see Figure 1).
If the phase-drift remains uncorrected, wave-like artifacts appear due to not corresponding phase values of the combined data of different repetitions. This effect is shown subsequently in the results section.
2.3 Motion phantom
The motion phantom is designed in a way that rigid object motion can be imaged with an MRI system.
It consists of an inner part, placed inside the MRI system, and an outer part, placed right next to the MRI system. The outer part consist of an electro motor RE50 (maxon motor GmbH, München, Germany) and a translation stage. It is used to generate a force to deform the inner part of the phantom. Both parts are connected via a polyethylene line. The electro motor rotates periodically and the rotation is transformed to a linear movement by a translation stage. This movement is transferred to the inner phantom part via the line and causes a translational relative movement between a small diameter test tube (length/inner diameter/outer diameter: 100/20/22 mm) and a fixed large diameter test tube (90/13/18 mm). To reverse the translation of the inner phantom part, an elastic band is used. The imaged object, an agarose cylinder, is placed in the small diameter test tube. A schematic overview of the phantom is shown in Figure 2. The agarose has a concentration of 0.75% and is prepared under laboratory conditions. Therefor, agarose powder (Natura, Hannover, Germany) was boiled in pure water for five minutes.
2.4 Triggering for retrospective gating
For retrospective gating, trigger signals of the phantom and the MRI system are needed to assign the measured data to the corresponding phantom motion status. The trigger signals are rectangular voltage changes and denoted as pulses in the following.
The motor control sends a trigger signal with a pulse width of 100 μs every 13.5° of motor rotation which corresponds to a linear movement of the motion phantom of 75 ms. With these pulses the motion of the phantom can be approximated. Therefore, the movement between two trigger signals is assumed to be linear. The result is shown in Figure 3.
The MRI system sends a trigger signal with a pulse width of 1 ms before every phase-coding step. The time for frequency coding is approximated as two times the echo time. With the used motor rotation speed of 100 rpm and an echo time of 6 ms this results in an object motion of 10 μm during frequency coding. Compared to a voxel size of 781 μm in the direction of movement, the motion during frequency encoding is considered negligible.
2.5 Image reconstruction and retrospective gating
During data acquisition, the k-space is fully Cartesian-sampled but the contained data is phase-drift- and motion-corrupted. Therefore, the phase-drift correction has to be carried out prior to reconstruction.
To show the influence of object motion on the reconstructed images, a 2D inverse fast Fourier transform is applied to the k-space data. In case of the retrospectively gated images, before applying a 2D inverse fast Fourier transform, data binning was executed as follows.
The motor position was discretized in intervals with an interval width of 0.2 mm. With a total moving distance of 2.5 mm this makes 13 bins for the gating process. Now the data for individual k-space lines are sorted into bins according to the motor position at the time of their echo middle. All bins contain data for all rows of k-space. The difference between data of different repetitions of a specified phase-encoding step in one bin can be considered negligible, so the repetition to be used can be chosen randomly within one bin.
3 Results and discussion
As it can be seen in Figure 4A, image reconstruction for an static agarose cylinder was performed successfully. With the chosen imaging parameters, an SNR of 32 could be achieved. It was calculated manually by the ratio of the means of object and background values.
In Figure 4B, an image of the first repetition which was corrupted due to rigid object motion can be seen. The mean square error (MSE) between the motion-corrupted image and the image without motion amounts to 140. Compared to a mean signal value in the region of interest of 150, the MSE indicates artifacts in the motion-corrupted image. Visual examination reveals noticeable artifacts in the direction of motion. The object appears to be shortened in the motion-corrupted image.
The results of the gating process can be seen in Figure 4D–F. The movement was split in 13 bins for gating and the positions 0 mm, 1 mm and 2 mm were chosen for visualization. The position 0 mm corresponds to a pause in the periodic movement which can be seen in Figure 3. Therefore, the MSE between the gated image (pos. 0 mm) and the image with no object motion should be zero for a perfect gating process. In fact the MSE amounts to 39. Visually, there can not be identified much difference. The reason for that is a difference in the mean values in the region of interest. That is caused by the not entirely perfect equilibrium state of the longitudinal magnetization. The motionless image was acquired within one repetition, whereas the gated image consists of the first seven repetitions. An observation of the raw data revealed that the signal at the beginning of a data acquisition is slightly higher than in the following repetitions.
For perfect gating, the image at position 1 mm should be a shifted version of Figure 4A and Figure 4D. Because of that, the image (pos. 1 mm) was manipulated applying a shift in the opposite direction and than the MSE to the image with no motion was calculated. The result is an MSE of 62. The fact of the small error reveals that a correction has been performed, but it seems as it has been less successful than in the case for the position 0 mm. This is probably caused by the phase-drift. A phase-drift correction was performed, which had a big impact that can be recognized by comparing Figure 4D–F with Figure 4C. Nevertheless, experiments with the data that is chosen for the gating process revealed that the gating result becomes better if the chosen data belongs to successive repetitions (data not shown).
For the gated image at position 2 mm, the results are similar to the position 1 mm.
3.2 Motion phantom
Analysis of the gated images revealed that there was an object motion of about 3 pixels. With a voxel width in the direction of movement of 0.78 mm this corresponds to 2.3 mm of object motion. Therefore, it can be assumed that almost the full movement of the translation stage could be transferred to the inner phantom. This validates the chosen phantom design.
In this work, a compact and easy to use phantom for rigid object motion was designed. It was shown that retrospective gating based on motor trigger signals can be performed. The object movement was reconstructed and can be used as the ground truth for motion estimation. The preprocessing of the raw data, e.g. phase-drift correction, was carried out succesfully. In further examinations, the gating process as implemented needs to be examined with elastic motion and therefore the motion phantom has to be modified. If successful, the data can be used as reference for blind motion compensation algorithms as presented in .
Research funding: The author state no funding involved. Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent has been obtained from all individuals included in this study. Ethical approval: The conducted research is not related to either human or animal use.
Tavallaei MA, Johnson PM, Liu J, Drangova M. Design and evaluation of an MRI-compatible linear motion stage. Med Phys. 2016;43:62–71. Google Scholar
Zaitsev M, Mclaren P, Herbst M. Motion artifacts in MRI: a complex problem with many partial solutions. J Magn Reson Imaging 2015;42:887–901. Google Scholar
Möller A, Maaß M, Mertins A. Blind sparse motion MRI with linear subpixel interpolation. Bildverarbeitung für die Medizin. 2015:510–5. Google Scholar
Nitz WR, Runge VR, Schmeets SH, Faulkner W, Desai NK. Praxiskurs MRT. 1st ed. Stuttgart: Georg Thieme Verlag KG 2007. Google Scholar
Liang ZP, Lauterbur PC. Principles of magnetic resonance imaging. 1st ed. New York: IEEE Press; 2000. Google Scholar
About the article
Published Online: 2016-09-30
Published in Print: 2016-09-01
Citation Information: Current Directions in Biomedical Engineering, Volume 2, Issue 1, Pages 471–474, ISSN (Online) 2364-5504, DOI: https://doi.org/10.1515/cdbme-2016-0104.
©2016 Max Schmiedel et al., licensee De Gruyter.. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0