Automated 3D Motion Tracking Using Gabor Filter Bank, Robust Point Matching, and Deformable Models

Tagged magnetic resonance imaging (tagged MRI or tMRI) provides a means of directly and noninvasively displaying the internal motion of the myocardium. Reconstruction of the motion field is needed to quantify important clinical information, e.g., the myocardial strain, and detect regional heart functional loss. In this paper, we present a three-step method for this task. First, we use a Gabor filter bank to detect and locate tag intersections in the image frames, based on local phase analysis. Next, we use an improved version of the robust point matching (RPM) method to sparsely track the motion of the myocardium, by establishing a transformation function and a one-to-one correspondence between grid tag intersections in different image frames. In particular, the RPM helps to minimize the impact on the motion tracking result of (1) through-plane motion and (2) relatively large deformation and/or relatively small tag spacing. In the final step, a meshless deformable model is initialized using the transformation function computed by RPM. The model refines the motion tracking and generates a dense displacement map, by deforming under the influence of image information, and is constrained by the displacement magnitude to retain its geometric structure. The 2D displacement maps in short and long axis image planes can be combined to drive a 3D deformable model, using the moving least square method, constrained by the minimization of the residual error at tag intersections. The method has been tested on a numerical phantom, as well as on in vivo heart data from normal volunteers and heart disease patients. The experimental results show that the new method has a good performance on both synthetic and real data. Furthermore, the method has been used in an initial clinical study to assess the differences in myocardial strain distributions between heart disease (left ventricular hypertrophy) patients and the normal control group. The final results show that the proposed method is capable of separating patients from healthy individuals. In addition, the method detects and makes possible quantification of local abnormalities in the myocardium strain distribution, which is critical for quantitative analysis of patients' clinical conditions. This motion tracking approach can improve the throughput and reliability of quantitative strain analysis of heart disease patients, and has the potential for further clinical applications.

[1]  Xiaoxu Wang,et al.  Meshless deformable models for LV motion analysis , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[2]  Scott Schaefer,et al.  Image deformation using moving least squares , 2006, ACM Trans. Graph..

[3]  Jerry L. Prince,et al.  Imaging heart motion using harmonic phase MRI , 2000, IEEE Transactions on Medical Imaging.

[4]  Dimitris N. Metaxas Physics-Based Deformable Models , 1996 .

[5]  Mark A Fleming,et al.  Meshless methods: An overview and recent developments , 1996 .

[6]  Alistair A. Young,et al.  Model Tags: Direct 3D Tracking of Heart Wall Motion from Tagged Magnetic Resonance Images , 1998, MICCAI.

[7]  Anand Rangarajan,et al.  A new point matching algorithm for non-rigid registration , 2003, Comput. Vis. Image Underst..

[8]  Leon Axel,et al.  Dense Myocardium Deformation Estimation for 2D Tagged MRI , 2005, FIMH.

[9]  L. Axel,et al.  MR imaging of motion with spatial modulation of magnetization. , 1989, Radiology.

[10]  Grace Wahba,et al.  Spline Models for Observational Data , 1990 .

[11]  Leon Axel,et al.  2D Motion Analysis of Long Axis Cardiac Tagged MRI , 2007, MICCAI.

[12]  James S. Duncan,et al.  Boundary element method-based regularization for recovering of LV deformation , 2007, Medical Image Anal..

[13]  Marc Alexa,et al.  Point based animation of elastic, plastic and melting objects , 2004, SCA '04.

[14]  Petia Radeva,et al.  Tag surface reconstruction and tracking of myocardial beads from SPAMM-MRI with parametric B-spline surfaces , 2001, IEEE Transactions on Medical Imaging.

[15]  Steven Gold,et al.  A Graduated Assignment Algorithm for Graph Matching , 1996, IEEE Trans. Pattern Anal. Mach. Intell..

[16]  Nicholas Ayache,et al.  Definition of a four-dimensional continuous planispheric transformation for the tracking and the analysis of left-ventricle motion , 1998, Medical Image Anal..

[17]  I.E. Magnin,et al.  A realistic anthropomorphic dynamic heart phantom , 2005, Computers in Cardiology, 2005.

[18]  Demetri Terzopoulos,et al.  A dynamic finite element surface model for segmentation and tracking in multidimensional medical images with application to cardiac 4D image analysis. , 1995, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[19]  Jinah Park,et al.  Volumetric deformable models with parameter functions: A new approach to the 3D motion analysis of the LV from MRI-SPAMM , 1995, Proceedings of IEEE International Conference on Computer Vision.

[20]  L. Axel,et al.  Using Gabor Filter Banks and Temporal-Spatial Constraints to Compute 3D Myocardium Strain , 2006, 2006 International Conference of the IEEE Engineering in Medicine and Biology Society.

[21]  P. Lancaster,et al.  Surfaces generated by moving least squares methods , 1981 .

[22]  N. Tustison,et al.  Biventricular myocardial strains via nonrigid registration of anatomical NURBS model [corrected]. , 2006, IEEE transactions on medical imaging.

[23]  Daniel Rueckert,et al.  Analysis of 3-D myocardial motion in tagged MR images using nonrigid image registration , 2004, IEEE Transactions on Medical Imaging.

[24]  Alejandro F. Frangi,et al.  Automatic construction of multiple-object three-dimensional statistical shape models: application to cardiac modeling , 2002, IEEE Transactions on Medical Imaging.

[25]  Dennis Gabor,et al.  Theory of communication , 1946 .