Variable-Density Sampling and Radial View-Ordering (VDRad)
For 3D Cartesian imaging, we focus the design of the data acquisition scheme on achieving motion-robustness and ideal undersampling patterns for accelerated imaging techniques. The technique that we have developed is the Variable Density sampling and Radial view-ordering, VDRad for short.
Because of the advantageous motion properties of radial imaging, we mimic the effects for 3D Cartesian imaging by grouping (ky,kz)-views into radial spokes. These spokes are ordered according to the Golden Ratio . Lastly, we incorporate variable-density sampling for an ideal undersampling pattern for accelerated imaging. An important feature in our method is that we provide flexibility to support arbitrary (both isotropic and anisotropic) sampling patterns to cater to different hardware configurations.
The generated view-ordering/sampling pattern is also ideal for time-resolved acquisition techniques where it is important to have unique sampling patterns at each time point.
The MATLAB scripts/functions are included here to provide a way to illustrate the VDRad algorithm. Different sampling patterns and (both isotropic and anisotropic) accelerations can be prescribed to cater to arbitrary coil array configurations.
Downloads – MATLAB scripts/functions demonstrating VDRad
(Updated 16 October 2014)
Nonrigid Motion Correction in 3D: Autofocusing
Because of the lengthy duration of typical MR scans, data acquired are sensitive to motion effects. There has been a lot of research on methods to minimize these motion artifacts – from prospective to retrospective methods. We wanted to focus on a method that can be generalized to any sequence. Thus, we worked on correcting for these artifacts retrospectively during the reconstruction process.
In order to do so, we looked into autofocusing (autocorrection). Some work has been published on this idea for the application of motion correction. By localizing the autofocusing process, we extend the concept to nonrigid motion artifacts. In order to make the algorithm computationally feasible, we limit the search space to motion paths to be linear paths collected with multi-channel receivers. The motion paths were estimated from navigator data collected using the Butterfly modification.
More details of the pulse sequence design and algorithm can be found in the references.
A complete MATLAB-based implementation can be downloaded. It is much slower than the OMP,CUDA,C++ version, but it demonstrates the key concepts of the algorithm. For more information, see the informative README file.
You can also download each file separately:
Alternatively, the entire package can be downloaded in one ZIP file. The additional MATLAB tools needed are also included.
(Added 24 Sept 2014 with old version)
In order to run some of the example code, you will need some additional MATLAB functions. You can download the additional tools here:
In the archived file, you will find some debugging tools and common functions needed for typical MR reconstructions