Multi-dimensional spline-based non-rigid image registration

Image registration, or equivalently motion estimation, plays a central role in a broad range of ultrasound applications including elastography, estimation of blood or tissue motion, radiation force imaging, and extended field of view imaging. Because of its central significance, motion estimation accuracy, precision, and computational cost are of critical importance. Furthermore, since motion estimation is typically performed on sampled signals, while estimates are usually desired over a continuous domain, performance should be considered in conjunction with associated interpolation. We have previously presented a highly accurate, spline-based time delay estimator that directly determines sub-sample time delay estimates from sampled data. The algorithm uses cubic splines to produce a continuous time representation of a reference signal and then computes an analytical matching function between this reference and a delayed signal. The location of the minima of this function yields estimates of the time delay. In this paper we describe a MUlti-dimensional Spline-based Estimator (MUSE) that allows accurate and precise estimation of multi-dimensional displacements/strain components from multi-dimensional data sets. In this paper we describe the mathematical formulation for three-dimensional (3D) motion/strain estimation and present simulation results to assess the intrinsic bias and standard deviation of this algorithm and compare it to currently available multi-dimensional estimators. In 1,000 noise-free simulations we found that 2D MUSE exhibits maximum bias errors of 4.8nm and 297nm in range and azimuth respectively. The maximum simulated standard deviation of estimates in both dimensions was comparable at 0.0026 samples (corresponding to 54nm axially and 378nm laterally). These results are two to three orders of magnitude lower than currently used 2D tracking methods. Simulation of performance in 3D yielded similar results to those observed in 2D. We also performed experiments using 2D MUSE on an Ultrasonix Sonix RP imaging system with an L14-5/38 linear array transducer operating at 6.6MHz. With this experimental data we found that bias errors were significantly smaller than geometric errors induced by machining of the transducer mount.

[1]  M. O’Donnell,et al.  Internal displacement and strain imaging using ultrasonic speckle tracking , 1994, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  Hans Torp,et al.  Comparison between cross-correlation and auto-correlation technique in color flow imaging , 1993 .

[3]  Gregg Trahey,et al.  Acoustic radiation force impulse imaging: in vivo demonstration of clinical feasibility. , 2002, Ultrasound in medicine & biology.

[4]  Gaetano Scarano,et al.  Discrete time techniques for time delay estimation , 1993, IEEE Trans. Signal Process..

[5]  Michael Unser,et al.  A pyramid approach to subpixel registration based on intensity , 1998, IEEE Trans. Image Process..

[6]  Gaetano Giunta,et al.  Fine estimators of two-dimensional parameters and application to spatial shift estimation , 1999, IEEE Trans. Signal Process..

[7]  W.F. Walker,et al.  P1A-6 A Novel Spline-Based Algorithm for Multidimensional Displacement and Strain Estimation , 2006, 2006 IEEE Ultrasonics Symposium.

[8]  M. O’Donnell,et al.  Phase-aberration correction using signals from point reflectors and diffuse scatterers: basic principles , 1988, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[9]  W F Walker,et al.  A method of imaging viscoelastic parameters with acoustic radiation force. , 2000, Physics in medicine and biology.

[10]  J. Ophir,et al.  Methods for Estimation of Subsample Time Delays of Digitized Echo Signals , 1995 .

[11]  B. Friemel,et al.  Relative performance of two-dimensional speckle-tracking techniques: normalized correlation, non-normalized correlation and sum-absolute-difference , 1995, 1995 IEEE Ultrasonics Symposium. Proceedings. An International Symposium.

[12]  J. Ophir,et al.  Elastography: A Quantitative Method for Imaging the Elasticity of Biological Tissues , 1991, Ultrasonic imaging.

[13]  W. Walker,et al.  A comparison of the performance of time-delay estimators in medical ultrasound , 2003, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[14]  W. O’Brien,et al.  Volumetric blood flow via time-domain correlation: experimental verification , 1990, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[15]  W.F. Walker,et al.  A spline-based algorithm for continuous time-delay estimation using sampled data , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[16]  W. Walker,et al.  Radiation force imaging of viscoelastic properties with reduced artifacts , 2003, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[17]  R C Waag,et al.  Measurements of ultrasonic pulse arrival time differences produced by abdominal wall specimens. , 1991, The Journal of the Acoustical Society of America.

[18]  P. Suetens,et al.  A simulation study on the performance of different estimators for two-dimensional velocity estimation , 2002, 2002 IEEE Ultrasonics Symposium, 2002. Proceedings..

[19]  W.F. Walker,et al.  Efficient and accurate spline-based time delay estimation , 2004, IEEE Ultrasonics Symposium, 2004.

[20]  Gregg E. Trahey,et al.  A comparison of algorithms for tracking sub-pixel speckle motion , 1997, 1997 IEEE Ultrasonics Symposium Proceedings. An International Symposium (Cat. No.97CH36118).

[21]  G. Trahey,et al.  Synthetic receive aperture imaging with phase correction for motion and for tissue inhomogeneities. II. Effects of and correction for motion , 1992, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[22]  O. Bonnefous,et al.  Time Domain Formulation of Pulse-Doppler Ultrasound and Blood Velocity Estimation by Cross Correlation , 1986, Ultrasonic imaging.