Spatio-temporal tracking of myocardial deformations with a 4-D B-spline model from tagged MRI

Myocardial kinematics from tagged MRI based on a 4-D B-spline model

Current research investigating the modeling of left ventricular dynamics for accurate clinical assessment of cardiac function is extensive. Magnetic resonance (MR) tagging is a functional imaging method which allows for encoding of a grid of signal voids on cardiac MR images, providing a mechanism for noninvasive measurement of intramural tissue deformations, in vivo. We present a novel technique of employing a four-dimensional (4-D) B-spline model which permits concurrent determination of myocardial beads and myocardial strains. The method entails fitting the knot planes of the 4-D B-spline model for fixed times to a sequence of triplets of orthogonal sets of tag surfaces for all imaged volumetric frames within the constraints of the model's spatio-temporal internal energy. From a three-dimensional (3-D) displacement field, the corresponding long and short-axis Lagrangian normal, shear, and principal strain maps are produced. As an important byproduct, the points defined by the 3-D intersections of the triplets of orthogonal tag planes, which we refer to as myocardial beads, can easily be determined by our model. Displaying the beads as a movie loop allows for the visualization of the nonrigid movement of the left ventricle in 3-D.

Construction of a statistical model for cardiac motion analysis using nonrigid image registration

In this paper we present a new technique for tracking the movement of the myocardium using a statistical model derived from the motion fields in the hearts of several healthy volunteers. To build the statistical model we tracked the motion of the myocardium in 17 volunteers using a nonrigid registration technique based on free-form deformations and mapped the motion fields obtained into a common reference coordinate system. A principal component analysis (PCA) was then performed on the motion fields to extract the major modes of variation in the fields between the successive time frames. The modes of variation obtained were then used to parametrize the free-form deformations and build our statistical model. The results of using our model to track the motion of the heart in normal volunteers are also presented.

Unsupervised reconstruction of a three-dimensional left ventricular strain from parallel tagged cardiac images

A new algorithm, called the Unsupervised Tag ExTraction and Heart strain(E) Reconstruction (UNTETHER) algorithm, is presented for quantifying three-dimensional (3D) myocardial strain in tagged cardiac MR images. Five human volunteers and five postinfarct patients were imaged. 3D strains measured by UNTETHER and a user-supervised technique were compared. Each study was analyzed in 49 +/- 8 min with UNTETHER, compared to approximately 4 hr with the user-supervised technique. For pooled human data, the correlation coefficient between the two methods for circumferential shortening (E(cc)) was r = 0.91 at the mid-wall (P < 0.0005). UNTETHER is capable of measuring wall motion abnormalities resulting from coronary artery disease, and has the potential to overcome the main limitations (time and user-supervision requirements) to routine clinical use of tagged cardiac MRI.

A discrete curvature-based deformable surface model with application to segmentation of volumetric images

In this paper, we present a new curvature-based three-dimensional (3-D) deformable surface model. The model deforms under defined force terms. Internal forces are calculated from local model curvature, using a robust method by a least-squares error (LSE) approximation to the Dupin indicatrix. External forces are calculated by applying a step expansion and restoration filter (SEF) to the image data. A solution for one of the most common problems associated with deformable models, self-cutting, has been proposed in this work. We use a principal axis analysis and reslicing of the deformable model, followed by triangulation of the slices, to remedy self-cutting. We use vertex resampling, multiresolution deformation, and refinement of the mesh grid to improve the quality of the model deformation, which leads to better results. Examples of the model application to different cases (simulation, magnetic resonance imaging (MRI), computerized tomography (CT), and ultrasound images) are presented, showing diversity and flexibility of the model.

Biomechanical dynamics of the heart with MRI

Magnetic resonance imaging (MRI) provides a noninvasive way to evaluate the biomechanical dynamics of the heart. MRI can provide spatially registered tomographic images of the heart in different phases of the cardiac cycle, which can be used to assess global cardiac function and regional endocardial surface motion. In addition, MRI can provide detailed information on the patterns of motion within the heart wall, permitting calculation of the evolution of regional strain and related motion variables within the wall. These show consistent patterns of spatial and temporal variation in normal subjects, which are affected by alterations of function due to disease. Although still an evolving technique, MRI shows promise as a new method for research and clinical evaluation of cardiac dynamics.

A MAP framework for tag line detection in SPAMM data using Markov random fields on the B-spline solid

Magnetic resonance (MR) tagging is a technique for measuring heart deformations through creation of a stripe grid pattern on cardiac images. In this paper, we present a maximum a posteriori (MAP) framework for detecting tag lines using a Markov random field (MRF) defined on the lattice generated by three-dimensional (3-D) and four-dimensional (4-D) (3-D + t) uniform sampling of B-spline models. In the 3-D case, MAP estimation is cast for detecting present tag features in the current image given an initial solid from the previous frame (the initial undeformed solid is manually positioned by clicking on corner points of a cube). The method also allows the parameters of the solid model, including the number of knots and the spline order, to be adjusted within the same framework. Fitting can start with a solid with less knots and lower spline order and proceed to one with more knots and/or higher order so as to achieve more accuracy and/or higher order of smoothness. In the 4-D case, the initial model is considered to be the linear interpolation of a sequence of optimal solids obtained from 3-D tracking. The same framework proposed for the 3-D case can once again be applied to arrive at a 4-D B-spline model with a higher temporal order.

Virtual tagging: numerical considerations and phantom validation

This paper presents a virtual tagging framework for measuring, as well as visualising, myocardial deformation using magnetic resonance (MR) velocity imaging. Tagging grids are allocated artificially according to the deformation gradient with varying shapes and densities. The control points are then deformed such that the difference between the induced deformation velocity and that of actually measured MR data is minimum. A full three-dimensional implementation of the technique combined with the mass conservation constraint is provided. Numerical considerations of applying the proposed framework and different optimization strategies have been investigated with both simulated and phantom experiments. The accuracy of the technique in terms of following material deformation is compared with that of conventional tagging technique.

Four-dimensional processing of deformable cardiac PET data

A four-dimensional deformable motion algorithm is described for use in the motion compensation of gated cardiac positron emission tomography. The algorithm makes use of temporal continuity and a non-uniform elastic material model to provide improved estimates of heart motion between time frames. Temporal continuity is utilized in two ways. First, incremental motion fields between adjacent time frames are calculated to improve estimation of long-range motion between distant time frames. Second, a consistency criterion is used to insure that the image match between distant time frames is consistent with the deformations used to match adjacent time frames. The consistency requirement augments the algorithm's ability to estimate motion between noisy time frames, and the concatenated incremental motion fields improve estimation for large deformations. The estimated motion fields are used to establish a voxel correspondence between volumes and to produce a motion-compensated composite volume.

Fast LV motion estimation using subspace approximation techniques

Cardiac motion estimation is very important in understanding cardiac dynamics and in noninvasive diagnosis of heart disease. Magnetic resonance (MR) imaging tagging is a technique for measuring heart deformations. In cardiac tagged MR images, a set of dark lines are noninvasively encoded within myocardial tissue providing the means for measurement of deformations of the heart. The points along tag lines measured in different frames and in different directions carry important information for determining the three-dimensional nonrigid movement of left ventricle. However, these measurements are sparse and, therefore, multidimensional interpolation techniques are needed to reconstruct a dense displacement field. In this paper, a novel subspace approximation technique is used to accomplish this task. We formulate the displacement estimation as a variational problem and then project the solution into spline subspaces. Efficient numerical methods are derived by taking advantages of B-spline properties. The proposed technique significantly improves our previous results reported in [3] with respect to computational time. The method is applied to a temporal sequence of two-dimensional images and is validated with simulated and in vivo heart data.

Tag surface reconstruction and tracking of myocardial beads from SPAMM-MRI with parametric B-spline surfaces

Magnetic resonance imaging (MRI) is unique in its ability to noninvasively and selectively alter tissue magnetization, and create tag planes intersecting image slices. The resulting grid of signal voids allows for tracking deformations of tissues in otherwise homogeneous-signal myocardial regions. In this paper, we propose a specific spatial modulation of magnetization (SPAMM) imaging protocol together with efficient techniques for measurement of three-dimensional (3-D) motion of material points of the human heart (referred to as myocardial beads) from images collected with the SPAMM method. The techniques make use of tagged images in orthogonal views by explicitly reconstructing 3-D B-spline surface representation of tag planes (tag planes in two orthogonal orientations intersecting the short-axis (SA) image slices and tag planes in an orientation orthogonal to the short-axis tag planes intersecting long-axis (LA) image slices). The developed methods allow for viewing deformations of 3-D tag surfaces, spatial correspondence of long-axis and short-axis image slice and tag positions, as well as nonrigid movement of myocardial beads as a function of time.

Three-dimensional modeling for functional analysis of cardiac images: a review

Three-dimensional (3-D) imaging of the heart is a rapidly developing area of research in medical imaging. Advances in hardware and methods for fast spatio-temporal cardiac imaging are extending the frontiers of clinical diagnosis and research on cardiovascular diseases. In the last few years, many approaches have been proposed to analyze images and extract parameters of cardiac shape and function from a variety of cardiac imaging modalities. In particular, techniques based on spatio-temporal geometric models have received considerable attention. This paper surveys the literature of two decades of research on cardiac modeling. The contribution of the paper is three-fold: 1) to serve as a tutorial of the field for both clinicians and technologists, 2) to provide an extensive account of modeling techniques in a comprehensive and systematic manner, and 3) to critically review these approaches in terms of their performance and degree of clinical evaluation with respect to the final goal of cardiac functional analysis. From this review it is concluded that whereas 3-D model-based approaches have the capability to improve the diagnostic value of cardiac images, issues as robustness, 3-D interaction, computational complexity and clinical validation still require significant attention.

Investigating intrinsic myocardial mechanics: the role of MR tagging, velocity phase mapping, and diffusion imaging

Assessment of myocardial mechanics is an integral part of understanding and predicting heart disease. This review covers the two most common magnetic resonance (MR) methods used to measure myocardial motion: myocardial tagging and myocardial velocity mapping. Myocardial tagging has been well established in clinical research, despite its time-consuming postprocessing procedure. Myocardial velocity mapping uses the phase shifts of the spins to encode the velocity into the MR signal. This means that once the myocardial contours have been segmented, the data can be automatically processed to obtain quantitative measurements. Diffusion MR also has found applications in cardiac imaging, with preliminary results of myocardial fiber architecture being obtained recently. These three different MR techniques have provided valuable insights into the assessment of intrinsic cardiac mechanics. J. Magn. Reson. Imaging 2000;12:873-883.