Improving IVUS palpography by incorporation of motion compensation based on block matching and optical flow

Intravascular ultrasound (IVUS) strain imaging of the luminal layer in coronary arteries, coined as IVUS palpography, utilizes conventional radio frequency (RF) signals acquired at 2 different levels of a compressional load. The signals are cross-correlated to obtain the microscopic tissue displacements, which can be directly translated into local strain of the vessel wall. However, (apparent) tissue motion and nonuniform deformation of the vessel wall, due to catheter wiggling, reduce signal correlation and result in invalid strain estimates. Implications of probe motion were studied on the tissue-mimicking phantom. The measured circumferential tissue displacement and level of the speckle decorrelation amounted to 12° and 0.58, respectively, for the catheter displacement of 456 μm. To compensate for the motion artifacts in IVUS palpography, a novel method based on the feature-based scale-space optical flow (OF), and classical block matching (BM) algorithm, were employed. The computed OF vector and BM displacement fields quantify the amount of local tissue misalignment in consecutive frames. Subsequently, the extracted circumferential displacements are used to realign the signals before strain computation. Motion compensation reduces the RF signal decorrelation and increases the number of valid strain estimates. The advantage of applying the motion correction in IVUS palpography was demonstrated in a midscale validation study on 14 in vivo pullbacks. Both methods substantially increase the number of valid strain estimates in the partial and compounded palpograms. Mean relative improvement in the number of valid strain estimates with motion compensation in comparison to one without motion compensation amounts to 28% and 14%, respectively. Implementation of motion compensation methods boosts the diagnostic value of IVUS palpography.     

An integrated compliant balloon ultrasound catheter for intravascular strain imaging

An integrated compliant balloon ultrasound catheter was developed to allow greater deformations in strain imaging with intravascular ultrasound. A 64-element circumferential array was placed inside a compliant silicone balloon catheter to capture real-time, phase-sensitive radio frequency (RF) data during deformation experiments. Strains over 40% could be applied to normal arterial wall tissue with intracatheter pressures as low as 200 kPa (2 atm). Strain images of a hard-soft rubber phantom, thrombus, and fibrotic plaque were produced using the integrated balloon ultrasound catheter. Results show that this catheter can apply large deformations at low pressures and image various vascular pathologies ex vivo. Potentially, it can serve as a multifunctional, intravascular therapeutic device to guide angioplasty and stent deployment.     

Nonlinear elasticity imaging: theory and phantom study

In tissue the Young's modulus cannot be assumed constant over a wide deformation range. For example, direct mechanical measurements on human prostate show up to a threefold increase in Young's modulus over a 10% deformation. In conventional elasticity imaging, these effects produce strain-dependent elastic contrast. Ignoring these effects generally leads to suboptimal contrast (stiffer tissues at lower strain are contrasted against softer tissues at higher strain), but measuring the nonlinear behavior results in enhanced tissue differentiation. To demonstrate the methods extracting nonlinear elastic properties, both simulations and measurements were performed on an agar-gelatin phantom. Multiple frames of phase-sensitive ultrasound data are acquired as the phantom is deformed by 12%. All interframe displacement data are brought back to the geometry of the first frame to form a three-dimensional (3-D) data set (depth, lateral, and preload dimensions). Data are fit to a 3-D second order polynomial model for each pixel that adjusts for deformation irregularities. For the phantom geometry and elastic properties considered in this paper, reconstructed frame-to-frame strain images using this model result in improved contrast to noise ratios (CNR) at all preload levels, without any sacrifice in spatial resolution. From the same model, strain hardening at all preload levels can be extracted. This is an independent contrast mechanism. Its maximum CNR occurs at 5.13% preload, and it is a 54% improvement over the best case (preload 10.6%) CNR for frame-to-frame strain reconstruction. Actual phantom measurements confirm the essential features of the simulation. Results show that modeling of the nonlinear elastic behavior has the potential to both increase detectability in elasticity imaging and provide a new independent mechanism for tissue differentiation.     

2-D high-frame-rate dynamic elastography using delay compensated and angularly compounded motion vectors: preliminary results

This paper describes a new ultrasound-based system for high-frame-rate measurement of periodic motion in 2-D for tissue elasticity imaging. Similarly to conventional 2-D flow vector imaging, the system acquires the RF signals from the region of interest at multiple steering angles. A custom sector subdivision technique is used to increase the temporal resolution while keeping the total acquisition time within the range suitable for real-time applications. Within each sector, 1-D motion is estimated along the beam direction. The intra- and inter-sector delays are compensated using our recently introduced delay compensation algorithm. In-plane 2-D motion vectors are then reconstructed from these delay-compensated 1-D motions. We show that Young's modulus images can be reconstructed from these 2-D motion vectors using local inversion algorithms. The performance of the system is validated quantitatively using a commercial flow phantom and a commercial elasticity phantom. At the frame rate of 1667 Hz, the estimated flow velocities with the system are in agreement with the velocity measured with a pulsed-wave Doppler imaging mode of a commercial ultrasound machine with manual angle correction. At the frame rate of 1250 Hz, phantom Young's moduli of 29, 6, and 54 kPa for the background, the soft inclusion, and the hard inclusion, are estimated to be 30, 11, and 53 kPa, respectively.     

Elasticity imaging using conventional and high-frame rate ultrasound imaging: experimental study

High-frame rate ultrasound imaging is necessary to track fast deformation in ultrasound elasticity imaging, but the image quality may be degraded. Previously, we investigated the performance of strain imaging using numerical models of conventional and ultrafast ultrasound imaging techniques. In this paper, we performed experimental studies to quantitatively evaluate the strain images and elasticity maps obtained using conventional and high frame rate ultrasound imaging methods. The experiments were carried out using point target and tissue mimicking phantoms. The experimental results were compared with the results of numerical simulation. Our experimental studies confirm that the signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and axial/lateral resolution of the displacement and strain images acquired using high-frame rate ultrasound imaging are slightly lower but comparable with those obtained using conventional imaging. Furthermore, the quality of elasticity images also exhibits similar trends. Thus, high-frame rate ultrasound imaging can be used reliably for static elasticity imaging to capture the internal tissue motion if the frame rate is critical.     

Elasticity reconstruction from displacement and confidence measures of a multi-compressed ultrasound RF sequence

Ultrasound elasticity imaging shows promise as a new way for early detection of cancers by assessing the elastic characteristics of soft tissue. So far the commonly used approach involves solving the so-called inverse elasticity problem of recovering elastic parameters from displacement measurements. We propose a finite-elementbased nonlinear scheme to estimate the elasticity distribution of soft tissue from multi-compressed ultrasound radio frequency (RF) data. An experimental ultrasound workstation has been developed to acquire multi-compressed data. A composite probe was employed as the compression plate. The contact forces and torques were acquired at the same time as imaging. Axial displacements under different static loads are estimated from the RF data before and after deformation using a cross-correlation technique. The confidence of displacement estimates is employed as a weighting factor in solving the objective function describing the inverse elasticity reconstruction problem. A novel splitand- merge strategy is employed over the image sequence in which strain images are used to provide a priori knowledge of the relative stiffness distribution of the tissue to constrain the inverse problem solution. The experimental study has allowed us to investigate the performance of our approach in the controlled environment of simulated and phantom data. For a simulated single inclusion model with 5% axial displacement estimation error, the L2-error between the target and the reconstructed Young's modulus was found to be about 1%. In vivo validation of the proposed method has been carried out and some preliminary results are presented.     

Strain imaging using conventional and ultrafast ultrasound imaging: numerical analysis

In elasticity imaging, the ultrasound frames acquired during tissue deformation are analyzed to estimate the internal displacements and strains. If the deformation rate is high, high-frame-rate imaging techniques are required to avoid the severe decorrelation between the neighboring ultrasound images. In these high-frame-rate techniques, however, the broader and less focused ultrasound beam is transmitted and, hence, the image quality is degraded. We quantitatively compared strain images obtained using conventional and ultrafast ultrasound imaging methods. The performance of the elasticity imaging was evaluated using custom-designed, numerical simulations. Our results demonstrate that signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and spatial resolutions in displacement and strain images acquired using conventional and ultrafast ultrasound imaging are comparable. This study suggests that the high-frame-rate ultrasound imaging can be reliably used in elasticity imaging if frame rate is critical     

Intravascular Raman spectroscopic catheter for molecular diagnosis of atherosclerotic coronary disease

An intravascular catheter for Raman spectroscopic detection and analysis of coronary atherosclerotic disease has been developed. The catheter, having an outer diameter of 2 mm, consisted of a side-view-type micro-Raman probe, an imaging fiber bundle, a working channel (injection drain), and a balloon. By inflating the balloon, the probe was brought close to the inner wall of a modeled blood flow system and detected a phantom target buried in the wall. Results obtained demonstrate the possibility of using the spectroscopic catheter for molecular diagnosis of coronary lesions.     

Frequency-domain-based strain estimation and high-frame-rate imaging for quasi-static elastography

In freehand elastography, quasi-static tissue compression is applied through the ultrasound probe, and the corresponding axial strain is estimated by calculating the time shift between consecutive echo signals. This calculation typically suffers from a poor signal-to-noise ratio or from the decorrelation between consecutive echoes resulting from an erroneous axial motion impressed by the operator. This paper shows that the quality of elastograms can be improved through the integration of two distinct techniques in the strain estimation procedure. The first technique evaluates the displacement of the tissue by analyzing the phases of the echo signal spectra acquired during compression. The second technique increases the displacement estimation robustness by averaging multiple displacement estimations in a high-frame-rate imaging system, while maintaining the typical elastogram frame-rate. The experimental results, obtained with the Ultrasound Advanced Open Platform (ULA-OP) and a cyst phantom, demonstrate that each of the proposed methods can independently improve the quality of elastograms, and that further improvements are possible through their combination.     

An elasticity microscope. Part II: Experimental results

For pt.I see ibid., vol.44, no.6, pp.1304-19 (1997). Initial experimental results from a 50 MHz elasticity microscope are shown. Using methods discussed previously, we present measured displacement and normal axial strain fields from a tissue mimicking phantom. Results from this homogenous gel are compared to a finite element simulation of the deformation experiment. The spatial resolution is estimated to be approximately 52 μm for axial displacements, and 71 μm for normal axial strains. These estimates were further tested by imaging a phantom containing a hard cylindrical inclusion with cross-sectional diameter of 265 μm. By examining the strain transition between regions in this image, the spatial resolution of the normal axial strain was verified to be at most 88 μm. A typical experiment produces peak normal axial strain around 3%. These experiments demonstrate the potential of high frequency ultrasound as a means for elasticity microscopy. Preliminary deformation experiments are presented on porcine epidermis     

Reduced peak-hopping artifacts in ultrasonic strain estimation using the Viterbi algorithm

Internal strain resulting from tissue deformation can be estimated by correlation processing of speckle patterns within complex (i.e., radio frequency) ultrasound images acquired during deformation. At large deformations, the magnitude of the correlation coefficient peak can be significantly lower than unity, so that random speckle correlations will exceed the true peak. This effect is called ?peak hopping? and causes significant errors in displacement and deformation estimates. Here we investigate the Viterbi algorithm, a dynamic programming procedure, to overcome peak-hopping artifacts by finding the most likely sequence of hidden states in a sequence of observed events. It is well suited to motion estimation in elasticity- imaging studies because adjacent tissue elements remain adjacent following deformation. Particularly, tissue elements along an ultrasonic beam in one image lie along a 3-D continuous curve in the next image instant. The observed event in this case is the correlation coefficient of a pixel at a certain displacement. Radio-frequency data were generated before and after deformation with an average strain of 6%. Simulations were performed for a homogenous medium and for a medium with a stiffer inclusion. Results show that Viterbi processing of speckle-tracking outputs can significantly reduce peak-hopping artifacts.     

An elasticity microscope. Part I: Methods

An elasticity microscope images tissue stiffness at fine resolution. Possible applications include dermatology, ophthalmology, pathology, and tissue engineering. In addition, if the resolution approaches cellular dimensions, then this system may be very useful in understanding tissue micromorphology. Elasticity images can be reconstructed from displacement and strain fields measured throughout the specimen during controlled external loading. High frequency ultrasound is used to obtain these images by tracking coherent speckle motion during deformation. In this paper, methods are presented to track speckle in two dimensions with near unity correlation coefficients using a high frequency, single element focused transducer. These techniques include improved means for speckle tracking. Procedures to control boundary conditions for consistent specimen deformation and scanning techniques required to obtain a plane-strain state in the imaging plane are also discussed. To test these methods, a 50 MHz elasticity microscope was constructed     

Clutter filters adapted to tissue motion in ultrasound color flow imaging

The quality of ultrasound color flow images is highly dependent on sufficient attenuation of the clutter signals originating from stationary and slowly moving tissue. Without sufficient clutter rejection, the detection of low velocity blood flow will be poor, and the velocity estimates will have a large bias. In some situations, e.g., when imaging the coronary arteries or when the operator moves the probe in search for small vessels, there is considerable movement of tissue. It has been suggested that clutter rejection can be improved by mixing down the signal with an estimate of the mean frequency prior to high pass filtering. In this paper, we compare this algorithm with several other adaptive clutter filtering algorithms using both experimental data and simulations. We found that realistic accelerations of the tissue have a large effect on the clutter rejection. The best results were obtained by mixing down the signal with non-constant phase increments estimated from the signal. This adapted the filter to a possibly accelerated tissue motion and produced a significant improvement in clutter rejection     

冠状动脉造影图像序列中估计二维运动及变形的弹性配准算法

提出一种由单面血管造影图像序列分析心脏冠状动脉运动和变形的方法，即利用弹性配准算法，将动脉运动的估计简化为对连续两帧图像的骨架像素的匹配．采用动态规划的方法进行配准，同时引入自回归建模，二者结合起来互相作用：动态规划为自回归模型参数的估计提供样本值；自回归模型为动态规划提供适当的成本函数。采用临床得到的单面冠状动脉造影图像序列对该算法进行了验证。实验结果表明应用该方法计算血管骨架图像中各点的光滑移动向量场，能够得到精确的结果。     

Filter-based compounded delay estimation with application to strain imaging

Ultrasonic wave interference produces local fluctuations in both the envelope, known as speckle, and phase of echoes. Furthermore, such fluctuations are correlated in space, and subsequent motion estimation from the envelope and/or phase signal produces patterned, correlated errors. Compounding, or combining information from multiple decorrelated looks, reduces such effects. We propose using a filter bank to create multiple looks to produce a compounded motion estimate. In particular, filtering in the lateral direction is shown to preserve delay estimation accuracy in the filtered sub-bands while creating decorrelation between sub-bands at the expense of some lateral resolution. For Gaussian apodization, we explicitly compute the induced signal decorrelation produced by Gabor filters. Furthermore, it is shown that lateral filtering is approximately equivalent to steering, in which filtered sub-bands correspond to signals extracted from shifted sub-apertures. Field II simulation of a point spread function verifies this claim. We use phase zero and its variants as displacement estimators for our compounded result. A simplified deformation model is used to provide computer simulations of deforming an elastic phantom. Simulations demonstrate root mean square error (RMSE) reduction in both displacement and strain of the compounded result over conventional and its laterally blurred versions. Then we apply the methods to experimental data using a commercial elastic phantom, demonstrating an improvement in strain SNR.     

一种改进的图像序列运动小目标检测方法

针对目标背景运动的情况,在随机哈夫变换运动检测方法(MDRHT)的基础上,提出了基于中心偏移的哈夫变换运动检测方法(Center-Biased MDHT)。它利用图像中显著的边界信息,有效地估计图像序列帧间图像背景的整体运动;利用图像背景运动的一致性,预先设定帧间位移的阈值,减少了预处理阶段的计算量;对前后两帧相邻图像进行差分处理,充分利用小目标在差分图像中成对出现的特点,达到了在目标检测阶段减少计算量的目的。     

High-frame rate, full-view myocardial elastography with automated contour tracking in murine left ventricles in vivo

Myocardial elastography is a novel method for noninvasively assessing regional myocardial function, with the advantages of high resolution and high precision. The purpose in this paper was to isolate the left ventricle from other structures for better displacement and strain visualization. Using a high-resolution (30 MHz) ultrasound system and a retrospective electrocardiogram (ECG)-gating technique, an extremely high frame rate (up to 8 kHz) was previously shown achievable for full-view (12-mm times 12-mm) myocardial elastography in the murine left ventricle. In vivo experiments were performed in anesthetized normal and infarcted mice [one day after left anterior descending (LAD) coronary artery ligation]. Radio frequency (RF) signals of the left ventricle (LV) in the long-axis view and the associated ECG were simultaneously acquired, with the ECG allowing gating of the RF signals. Incremental axial displacement of the myocardium was estimated using a one-dimensional (1-D) cross-correlation function. The cumulative displacement and strain then were calculated from the incremental displacement. In this paper, after manual selection of 40-50 points along the endo-and epicardial borders in the first frame of the cine-loop, myocardial contour was automatically tracked across the entire LV throughout a full cardiac cycle, which correctly determined the region of interest (ROI) for better interpretation. The cine-loop of the cumulative displacement and strain in one cardiac cycle, in both the normal and infarcted cases, showed that motion and deformation in the infarcted myocardium were significantly reduced, and that the infarcted region underwent thinning, rather than thickening, during systole. High precision of the displacement estimation, due to high frequency (30 MHz) and high frame rate (up to 8 kHz) available with this system, allowed for automated tracking of a manually-initialized myocardial contour over an entire cardiac cycle. High frame rate, full-view myocardial elastography with automated contour tracking could provide regional strain information of the LV throughout an entire cardiac cycle, and characterize normal as well as detect abnormal myocardial function, such as an infarction. The method of automated contour tracking can further enhance the capability of the elastographic technique with minimal user intervention while providing accurate functional information for the detection of disease throughout the entire cardiac cycle.     

Coronary arterial dynamics computation with medical-image-based time-dependent anatomical models and element-based zero-stress state estimates

We propose a method for coronary arterial dynamics computation with medical-image-based time-dependent anatomical models. The objective is to improve the computational analysis of coronary arteries for better understanding of the links between the atherosclerosis development and mechanical stimuli such as endothelial wall shear stress and structural stress in the arterial wall. The method has two components. The first one is element-based zero-stress (ZS) state estimation, which is an alternative to prestress calculation. The second one is a “mixed ZS state” approach, where the ZS states for different elements in the structural mechanics mesh are estimated with reference configurations based on medical images coming from different instants within the cardiac cycle. We demonstrate the robustness of the method in a patient-specific coronary arterial dynamics computation where the motion of a thin strip along the arterial surface and two cut surfaces at the arterial ends is specified to match the motion extracted from the medical images.     

Motion artifact reduction for IVUS-based thermal strain imaging

Thermal strain imaging (TSI) using intravascular ultrasound (IVUS) has the potential to identify lipid pools within rupture-prone arterial plaques and serve as a valuable supplement to current IVUS systems in diagnosing acute coronary syndromes. The major challenge for in vivo application of TSI will be cardiac motion, including bulk motion and tissue deformation. Simulations based on an artery model, including a lipid-filled plaque, demonstrate that effective bulk motion compensation can be achieved within a certain motion range using spatial interpolation. We also propose a practical imaging scheme to minimize mechanical strains caused by tissue deformation based on a linear least squares fitting strategy. This scheme was tested on clinical data by artificially superimposing thermal displacements corresponding to different temperature rises. Results suggest a 1-2 degrees C temperature rise is required to detect lipids in an atherosclerotic plaque in vivo.