Motion estimation of tagged cardiac magnetic resonance images using variational techniques

A Novel Hierarchical Template Matching Model for Cardiac Motion Estimation

Cardiovascular disease diagnosis and prognosis can be improved by measuring patient-specific in-vivo local myocardial strain using Magnetic Resonance Imaging. Local myocardial strain can be determined by tracking the movement of sample muscles points during cardiac cycle using cardiac motion estimation model. The tracking accuracy of the benchmark Free Form Deformation (FFD) model is greatly affected due to its dependency on tunable parameters and regularisation function. Therefore, Hierarchical Template Matching (HTM) model, which is independent of tunable parameters, regularisation function, and image-specific features, is proposed in this article. HTM has dense and uniform points correspondence that provides HTM with the ability to estimate local muscular deformation with a promising accuracy of less than half a millimetre of cardiac wall muscle. As a result, the muscles tracking accuracy has been significantly (p < 0.001) improved (30%) compared to the benchmark model. Such merits of HTM provide reliably calculated clinical measures which can be incorporated into the decision-making process of cardiac disease diagnosis and prognosis.

Stochastic optimization with randomized smoothing for image registration

Image registration is typically formulated as an optimization process, which aims to find the optimal transformation parameters of a given transformation model by minimizing a cost function. Local minima may exist in the optimization landscape, which could hamper the optimization process. To eliminate local minima, smoothing the cost function would be desirable. In this paper, we investigate the use of a randomized smoothing (RS) technique for stochastic gradient descent (SGD) optimization, to effectively smooth the cost function. In this approach, Gaussian noise is added to the transformation parameters prior to computing the cost function gradient in each iteration of the SGD optimizer. The approach is suitable for both rigid and nonrigid registrations. Experiments on synthetic images, cell images, public CT lung data, and public MR brain data demonstrate the effectiveness of the novel RS technique in terms of registration accuracy and robustness.

Myocardial motion estimation of tagged cardiac magnetic resonance images using tag motion constraints and multi-level b-splines interpolation

Myocardial motion estimation of tagged cardiac magnetic resonance (TCMR) images is of great significance in clinical diagnosis and the treatment of heart disease. Currently, the harmonic phase analysis method (HARP) and the local sine-wave modeling method (SinMod) have been proven as two state-of-the-art motion estimation methods for TCMR images, since they can directly obtain the inter-frame motion displacement vector field (MDVF) with high accuracy and fast speed. By comparison, SinMod has better performance over HARP in terms of displacement detection, noise and artifacts reduction. However, the SinMod method has some drawbacks: 1) it is unable to estimate local displacements larger than half of the tag spacing; 2) it has observable errors in tracking of tag motion; and 3) the estimated MDVF usually has large local errors. To overcome these problems, we present a novel motion estimation method in this study. The proposed method tracks the motion of tags and then estimates the dense MDVF by using the interpolation. In this new method, a parameter estimation procedure for global motion is applied to match tag intersections between different frames, ensuring specific kinds of large displacements being correctly estimated. In addition, a strategy of tag motion constraints is applied to eliminate most of errors produced by inter-frame tracking of tags and the multi-level b-splines approximation algorithm is utilized, so as to enhance the local continuity and accuracy of the final MDVF. In the estimation of the motion displacement, our proposed method can obtain a more accurate MDVF compared with the SinMod method and our method can overcome the drawbacks of the SinMod method. However, the motion estimation accuracy of our method depends on the accuracy of tag lines detection and our method has a higher time complexity.

Hybrid intensity- and phase-based optical flow tracking of tagged MRI

Accurate tracking of the myocardium tissues in tagged Magnetic Resonance Images (MRI) is essential for evaluating the cardiac function. Current tracking methods utilize either the image intensity or the image phase as landmarks that can be tracked. In either case, the performance is vulnerable to the image quality and the fading of the tag lines. In this work, we propose a hybrid optical flow tracking method that combines both the intensity and the phase features of the image. The method is validated using numerical cardiac phantom as well as real MRI data experiments. Both experiments showed that the proposed method outperforms current intensity-based optical flow tracking and the phase-based HARP method with maximum error of 1 pixel at extreme conditions of tag fading.

Estimating dense cardiac 3D motion using sparse 2D tagged MRI cross-sections

In this work, we describe a new method, an extension of the Large Deformation Diffeomorphic Metric Mapping to estimate three-dimensional deformation of tagged Magnetic Resonance Imaging Data. Our approach relies on performing non-rigid registration of tag planes that were constructed from set of initial reference short axis tag grids to a set of deformed tag curves. We validated our algorithm using in-vivo tagged images of normal mice. The mapping allows us to compute root mean square distance error between simulated tag curves in a set of long axis image planes and the acquired tag curves in the same plane. Average RMS error was 0.31 ± 0.36(SD) mm, which is approximately 2.5 voxels, indicating good matching accuracy.

Trans-dimensional MCMC methods for fully automatic motion analysis in tagged MRI

Tagged magnetic resonance imaging (tMRI) is a well-known noninvasive method allowing quantitative analysis of regional heart dynamics. Its clinical use has so far been limited, in part due to the lack of robustness and accuracy of existing tag tracking algorithms in dealing with low (and intrinsically time-varying) image quality. In this paper, we propose a novel probabilistic method for tag tracking, implemented by means of Bayesian particle filtering and a trans-dimensional Markov chain Monte Carlo (MCMC) approach, which efficiently combines information about the imaging process and tag appearance with prior knowledge about the heart dynamics obtained by means of non-rigid image registration. Experiments using synthetic image data (with ground truth) and real data (with expert manual annotation) from preclinical (small animal) and clinical (human) studies confirm that the proposed method yields higher consistency, accuracy, and intrinsic tag reliability assessment in comparison with other frequently used tag tracking methods.

Reversible jump MCMC methods for fully automatic motion analysis in tagged MRI

Tagged magnetic resonance imaging (tMRI) is a well-known noninvasive method for studying regional heart dynamics. It offers great potential for quantitative analysis of a variety of kine(ma)tic parameters, but its clinical use has so far been limited, in part due to the lack of robustness and accuracy of existing tag tracking algorithms in dealing with low (and intrinsically time-varying) image quality. In this paper, we evaluate the performance of four frequently used concepts found in the literature (optical flow, harmonic phase (HARP) magnetic resonance imaging, active contour fitting, and non-rigid image registration) for cardiac motion analysis in 2D tMRI image sequences, using both synthetic image data (with ground truth) and real data from preclinical (small animal) and clinical (human) studies. In addition we propose a new probabilistic method for tag tracking that serves as a complementary step to existing methods. The new method is based on a Bayesian estimation framework, implemented by means of reversible jump Markov chain Monte Carlo (MCMC) methods, and combines information about the heart dynamics, the imaging process, and tag appearance. The experimental results demonstrate that the new method improves the performance of even the best of the four previous methods. Yielding higher consistency, accuracy, and intrinsic tag reliability assessment, the proposed method allows for improved analysis of cardiac motion.

Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging.