No abstract

Image registration and motion estimation play central roles in many fields, including RADAR, SONAR, light microscopy, and medical imaging. Because of its central significance, estimator accuracy, precision, and computational cost are of critical importance. 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 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 the MUlti-dimensional Spline-based Estimator (MUSE) that allows accurate and precise estimation of multidimensional displacements/strain components from multidimensional data sets. We describe the mathematical formulation for two-and three-dimensional 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 1000 noise-free simulations of ultrasound data we found that 2D MUSE exhibits maximum bias of 2.6 × 10 −4 samples in range and 2.2 × 10 −3 samples in azimuth (corresponding to 4.8 and 297 nm, respectively). The maximum simulated standard deviation of estimates in both dimensions was comparable at roughly 2.8 × 10 −3 samples (corresponding to 54 nm axially and 378 nm laterally). These results are between two and three orders of magnitude better than currently used 2D tracking methods. Simulation of performance in 3D yielded similar results to those observed in 2D. We also present experimental results obtained using 2D MUSE on data acquired by an Ultrasonix Sonix RP imaging system with an L14-5/38 linear array transducer operating at 6.6 MHz. While our validation of the algorithm was performed using ultrasound data, MUSE is broadly applicable across imaging applications.

The portability, low cost, and non-ionizing radiation associated with medical ultrasound suggest that it has potential as a superior alternative to X-ray for bone imaging. However, when conventional ultrasound imaging systems are used for bone imaging, clinical acceptance is frequently limited by artifacts derived from reflections occurring away from the main axis of the acoustic beam. In this paper, the physical source of off-axis artifacts and the effect of transducer geometry on these artifacts are investigated in simulation and experimental studies. In agreement with diffraction theory, the sampled linear-array geometry possessed increased off-axis energy compared with single-element piston geometry, and therefore, exhibited greater levels of artifact signal. Simulation and experimental results demonstrated that the linear-array geometry exhibited increased artifact signal when the center frequency increased, when energy off-axis to the main acoustic beam (i.e., grating lobes) was perpendicularly incident upon off-axis surfaces, and when off-axis surfaces were specular rather than diffusive. The simulation model used to simulate specular reflections was validated experimentally and a correlation coefficient of 0.97 between experimental and simulated peak reflection contrast was observed. In ex vivo experiments, the piston geometry yielded 4 and 6.2 dB average contrast improvement compared with the linear array when imaging the spinous process and interlaminar space of an animal spine, respectively. This work indicates that off-axis reflections are a major source of ultrasound image artifacts, particularly in environments comprising specular reflecting (i.e., bone or bone-like) objects. Transducer geometries with reduced sensitivity to off-axis surface reflections, such as a piston transducer geometry, yield significant reductions in image artifact.

Technical difficulties, poor image quality and reliance on pattern identifications represent some of the drawbacks of two-dimensional ultrasound imaging of spinal bone anatomy. To overcome these limitations, we sought to develop real-time volumetric imaging of the spine using a portable handheld device. The device measured 19.2 cm x 9.2 cm x 9.0 cm and imaged at 5 MHz center frequency. 2D imaging under conventional ultrasound and volumetric (3D) imaging in real time was achieved and verified by inspection using a custom spine phantom. Further device performance was assessed and revealed a 75-minute battery life and average frame rate of 17.7 Hz in volumetric imaging mode. Our results suggest that real-time volumetric imaging of the spine is a feasible technique for more intuitive visualization of the spine. These results may have important ramifications for a large array of neuraxial procedures.

Two-dimensional arrays present significant beamforming computational challenges because of their high channel count and data rate. These challenges are even more stringent when incorporating a 2-D transducer array into a battery-powered hand-held device, placing significant demands on power efficiency. Previous work in sonar and ultrasound indicates that 2-D array beamforming can be decomposed into two separable line-array beamforming operations. This has been used in conjunction with frequency-domain phase-based focusing to achieve fast volume imaging. In this paper, we analyze the imaging and computational performance of approximate near-field separable beamforming for high-quality delay-and-sum (DAS) beamforming and for a low-cost, phaserotation-only beamforming method known as direct-sampled in-phase quadrature (DSIQ) beamforming. We show that when high-quality time-delay interpolation is used, separable DAS focusing introduces no noticeable imaging degradation under practical conditions. Similar results for DSIQ focusing are observed. In addition, a slight modification to the DSIQ focusing method greatly increases imaging contrast, making it comparable to that of DAS, despite having a wider main lobe and higher side lobes resulting from the limitations of phase-only time-delay interpolation. Compared with non-separable 2-D imaging, up to a 20-fold increase in frame rate is possible with the separable method. When implemented on a smart-phone-oriented processor to focus data from a 60 × 60 channel array using a 40 × 40 aperture, the frame rate per C-mode volume slice increases from 16 to 255 Hz for DAS, and from 11 to 193 Hz for DSIQ. Energy usage per frame is similarly reduced from 75 to 4.8 mJ/ frame for DAS, and from 107 to 6.3 mJ/frame for DSIQ. We also show that the separable method outperforms 2-D FFT-based focusing by a factor of 1.64 at these data sizes. This data indicates that with the optimal design choices, separable 2-D beamforming can significantly improve frame rate and battery life for hand-held devices with 2-D arrays.