This subproject is one of many research subprojects utilizing theresources provided by a Center grant funded by NIH/NCRR. The subproject andinvestigator (PI) may have received primary funding from another NIH source,and thus could be represented in other CRISP entries. The institution listed isfor the Center, which is not necessarily the institution for the investigator.Introduction. Involuntary patient motion is still a great challenge in MRI. Specifically, in the elderly and pediatric patient population or in patients whose medical conditions (tremor, seizure, stroke) preclude them to hold still, effective strategies to compensate for motion are paramount. In this study, a variant of parallel imaging is introduced that can correct k-space inconsistencies arising from rigid body motion (rotation or translation). This motion correction scheme first identifies the degree of motion, corrects the k-space data accordingly and thereafter employs an augmented conjugate gradient based iterative image reconstruction to synthesize missing data in k-space. The method is described and verified in simulated interleaved EPI and spiral images scans as well as in vivo using bi-density spiral scanning.Materials and Methods: Reconstruction  Generally, an object rotation in image space is paralleled by a similar rotation of k-space data, whereas translations are reflected by linear phase rolls. If these motion components are known, k-space data can be corrected for but usually leading to a fragmentation of k-space. This, in turn, gives rise to significant ghost artifacts in the final image. Our correction builds upon an augmented version of an iterative SENSE reconstruction1 and is performed as follows: 1) counter-rotating k-space data by applying the corresponding rotation matrix to the k-space trajectory coordinate points of each profile/interleave prior to gridding. 2) Rotating the coil sensitivity map that enters the encoding matrix E1 for each profile/interleave. This rotation is necessary because even if the object is rotated back to its desired position, different regions of the object have been exposed to different coil sensitivities during the acquisition. 3) Correcting the altered sampling density after rotation. In this study, Voronoi tessellation has been used to derive the new sampling density from the rotated k-space trajectories. 4) Phasing the data to account for translation by applying the correction term pcorr(?) = exp{-j(2??x/FOVx) (kx(?)/[kx,max-kx,min])  j(2??y/FOVy)(ky(?)/[ky,max-ky,min])} to the original k-space data prior to gridding. Motion detection  Various methods exist to derive the extent of translational and rotational motion from MR data. In this study, the motion information was extracted from navigator echoes. The navigator information can be derived from the scan trajectory itself (i.e. self-navigating trajectories) or alternatively from a separate acquisition that provides a low resolution image. Here, a multi-grid registration approach was used that finds the maximum Pearson correlation between a reference image and individual navigator images and provided a reliable estimate of the amount of rotation and translation relative to the reference image (average over all images). To increase robustness and to improve the accuracy of co-registration this step was repeated at least twice. Experiments  Synthetic data for interleaved spiral and EPI acquisitions (8 interleaves) were generated by using inverse gridding operations2 on a motion corrupted phantom. For each of the eight interleaves a random head rotation (range 30 ) and translation (range 15mm) was generated. Prior to the inverse gridding step, each of the individually rotated and shifted images were multiplied by coil sensitivities simulating receiver coil sensitivities from six coils that were attached around the circumference of the object. In vivo validation was performed in 3 healthy volunteers using T2w spin echo scans with an interleaved spiral-in/spiral-out readout and an 8-channel head coil. The spiral-in part (3-5ms duration) provided for each interleaf data a low resolution navigator image (322). The spiral-out part was a normal interleaved spiral acquisition: TR/TE=4,000ms/85ms, slice thickness/ gap=4/1mm, 17 slices, FOV=24cm, matrix=256, interleaves = 32, and NEX=1. The receiver bandwidth for the spiral acquisition was +/- 125kHz. During each experiment the volunteers were asked to rotate and/or shift their heads at three increasing levels of motion (no, mild [~ 15 ], and moderate [~ 25 ] motion ).References: 1Pruessmann K, et al. MRM 46: 638-51, 2001; 2Rasche V, et al. IEEE TMI 18: 385-92, 1999. Acknowledgements: This work was supported in part by the NIH (1R01EB002771), the Center of Advanced MR Technology at Stanford (P41RR09784), Lucas Foundation.

National Institute of Health (NIH)
National Center for Research Resources (NCRR)
Biotechnology Resource Grants (P41)
Project #
Application #
Study Section
Special Emphasis Panel (ZRG1-SBIB-F (40))
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
Stanford University
Schools of Medicine
United States
Zip Code
Maclaren, Julian; Aksoy, Murat; Ooi, Melvyn B et al. (2018) Prospective motion correction using coil-mounted cameras: Cross-calibration considerations. Magn Reson Med 79:1911-1921
Guo, Jia; Holdsworth, Samantha J; Fan, Audrey P et al. (2018) Comparing accuracy and reproducibility of sequential and Hadamard-encoded multidelay pseudocontinuous arterial spin labeling for measuring cerebral blood flow and arterial transit time in healthy subjects: A simulation and in vivo study. J Magn Reson Imaging 47:1119-1132
Tamir, Jonathan I; Uecker, Martin; Chen, Weitian et al. (2017) T2 shuffling: Sharp, multicontrast, volumetric fast spin-echo imaging. Magn Reson Med 77:180-195
Lai, Lillian M; Cheng, Joseph Y; Alley, Marcus T et al. (2017) Feasibility of ferumoxytol-enhanced neonatal and young infant cardiac MRI without general anesthesia. J Magn Reson Imaging 45:1407-1418
Taviani, Valentina; Alley, Marcus T; Banerjee, Suchandrima et al. (2017) High-resolution diffusion-weighted imaging of the breast with multiband 2D radiofrequency pulses and a generalized parallel imaging reconstruction. Magn Reson Med 77:209-220
Uecker, Martin; Lustig, Michael (2017) Estimating absolute-phase maps using ESPIRiT and virtual conjugate coils. Magn Reson Med 77:1201-1207
Kogan, Feliks; Hargreaves, Brian A; Gold, Garry E (2017) Volumetric multislice gagCEST imaging of articular cartilage: Optimization and comparison with T1rho. Magn Reson Med 77:1134-1141
Aksoy, Murat; Maclaren, Julian; Bammer, Roland (2017) Prospective motion correction for 3D pseudo-continuous arterial spin labeling using an external optical tracking system. Magn Reson Imaging 39:44-52
Bian, W; Tranvinh, E; Tourdias, T et al. (2016) In Vivo 7T MR Quantitative Susceptibility Mapping Reveals Opposite Susceptibility Contrast between Cortical and White Matter Lesions in Multiple Sclerosis. AJNR Am J Neuroradiol 37:1808-1815
Vos, Sjoerd B; Aksoy, Murat; Han, Zhaoying et al. (2016) Trade-off between angular and spatial resolutions in in vivo fiber tractography. Neuroimage 129:117-132

Showing the most recent 10 out of 446 publications