Elastographic image quality vs. tissue motion in vivo

Noncontact longitudinal shear wave imaging for the evaluation of heterogeneous porcine brain biomechanical properties using optical coherence elastography

High-resolution quantification of heterogeneous brain biomechanical properties has long been an important topic. Longitudinal shear waves (LSWs) can be used to assess the longitudinal Young's modulus, but contact excitation methods have been used in most previous studies. We propose an air-coupled ultrasound transducer-based optical coherence elastography (AcUT-OCE) technique for noncontact excitation and detection of LSWs in samples and assessment of the nonuniformity of the brain's biomechanical properties. The air-coupled ultrasonic transducer (AcUT) for noncontact excitation of LSWs in the sample has a center frequency of 250 kHz. Phase-resolved Doppler optical coherence tomography (OCT) was used to image and reconstruct the propagation behavior of LSWs and surface ultrasound waves at high resolution. An agar phantom model was used to verify the feasibility of the experimental protocol, and experiments with ex vivo porcine brain samples were used to assess the nonuniformity of the brain biomechanical properties. LSWs with velocities of 0.83 ± 0.11 m/s were successfully excited in the agar phantom model. The perivascular elastography results in the prefrontal cortex (PFC) of the ex vivo porcine brains showed that the Young's modulus was significantly higher in the longitudinal and transverse directions on the left side of the cerebral vessels than on the right side and that the Young's modulus of the PFC decreased with increasing depth. The AcUT-OCE technique, as a new scheme for LSW applications in in vivo elastography, can be used for noncontact excitation of LSWs in brain tissue and high-resolution detection of heterogeneous brain biomechanical properties.

An unsupervised learning approach to ultrasound strain elastography with spatio-temporal consistency

Quasi-static ultrasound elastography (USE) is an imaging modality that measures deformation (i.e. strain) of tissue in response to an applied mechanical force. In USE, the strain modulus is traditionally obtained by deriving the displacement field estimated between a pair of radio-frequency data. In this work we propose a recurrent network architecture with convolutional long-short-term memory decoder blocks to improve displacement estimation and spatio-temporal continuity between time series ultrasound frames. The network is trained in an unsupervised way, by optimising a similarity metric between the reference and compressed image. Our training loss is also composed of a regularisation term that preserves displacement continuity by directly optimising the strain smoothness, and a temporal continuity term that enforces consistency between successive strain predictions. In addition, we propose an open-access in vivo database for quasi-static USE, which consists of radio-frequency data sequences captured on the arm of a human volunteer. Our results from numerical simulation and in vivo data suggest that our recurrent neural network can account for larger deformations, as compared with two other feed-forward neural networks. In all experiments, our recurrent network outperformed the state-of-the-art for both learning-based and optimisation-based methods, in terms of elastographic signal-to-noise ratio, strain consistency, and image similarity. Finally, our open-source code provides a 3D-slicer visualisation module that can be used to process ultrasound RF frames in real-time, at a rate of up to 20 frames per second, using a standard GPU.

Fast Strain Estimation and Frame Selection in Ultrasound Elastography Using Machine Learning

Ultrasound elastography aims to determine the mechanical properties of the tissue by monitoring tissue deformation due to internal or external forces. Tissue deformations are estimated from ultrasound radio frequency (RF) signals and are often referred to as time delay estimation (TDE). Given two RF frames I₁ and I₂ , we can compute a displacement image, which shows the change in the position of each sample in I₁ to a new position in I₂ . Two important challenges in TDE include high computational complexity and the difficulty in choosing suitable RF frames. Selecting suitable frames is of high importance because many pairs of RF frames either do not have acceptable deformation for extracting informative strain images or are decorrelated and deformation cannot be reliably estimated. Herein, we introduce a method that learns 12 displacement modes in quasi-static elastography by performing principal component analysis (PCA) on displacement fields of a large training database. In the inference stage, we use dynamic programming (DP) to compute an initial displacement estimate of around 1% of the samples and then decompose this sparse displacement into a linear combination of the 12 displacement modes. Our method assumes that the displacement of the whole image could also be described by this linear combination of principal components. We then use the GLobal Ultrasound Elastography (GLUE) method to fine-tune the result yielding the exact displacement image. Our method, which we call PCA-GLUE, is more than 10× faster than DP in calculating the initial displacement map while giving the same result. This is due to converting the problem of estimating millions of variables in DP into a much simpler problem of only 12 unknown weights of the principal components. Our second contribution in this article is determining the suitability of the frame pair I₁ and I₂ for strain estimation, which we achieve by using the weight vector that we calculated for PCA-GLUE as an input to a multilayer perceptron (MLP) classifier. We validate PCA-GLUE using simulation, phantom, and in vivo data. Our classifier takes only 1.5 ms during the testing phase and has an F1-measure of more than 92% when tested on 1430 instances collected from both phantom and in vivo data sets.

Thermal strain imaging in vivo using interpolated IQ-images

Thermal strain imaging (TSI) is a promising technique for ultrasonic thermometry, especially in the applications of thermal therapies. The accuracy of TSI is dependent on the sampling rate and line density of B-Scan images, and the prevalent IQ-demodulated ultrasound data outputted from low- and middle-end machines are therefore insufficient. Here, the feasibility of using interpolated IQ images for TSI (based on the "infinitesimal echo strain filter" model) is studied through in vivo experiments targeting the perirenal fat of pigs. It is demonstrated that, axial interpolations, especially those using the zero-padding algorithm, can recover the capabilities of the low-sampling-rate complex IQ images in TSI, and make their performances comparable to those of RF/IQ complex images with higher sample rate. Meanwhile, interpolations along the lateral direction can increase the line density of IQ images, reduce TSI errors, and reveal more details in the temperature maps. In the experiments, the variation in the thermometry coefficient (the k-value) is well below 3%. The findings here bring down the requirement of high sampling rate as well as high line density of US images in TSI, making it possible to be applied on common US machines.

Effect of Calcifications on Shear-Wave Elastography in Evaluating Breast Lesions

This study aimed to investigate the effect of calcifications on shear-wave elastography in evaluating breast lesions. We retrospectively reviewed ultrasound images of 673 breast lesions and compared the elasticity between lesions with and without calcifications in three subgroups: benign lesions, in situ carcinomas and invasive carcinomas. Breast lesions were confirmed histologically (n = 401) or by follow-up images for more than 2 y (n = 272). Calcifications were present in 25.3% (170/673) of the lesions. The E_mean values with and without calcifications, respectively, were as follows: 62.8 and 29.8 kPa in benign lesions (p = 0.000), 114.6 and 52.8 kPa in in situ carcinomas (p = 0.037) and 171.9 and 146.4 kPa in invasive carcinomas (p = 0.018). The presence of calcifications significantly increased the E_mean of breast lesions. Shear-wave elastography should be carefully interpreted in benign lesions with calcifications and in situ carcinomas without calcifications.

Progress in the Application of Ultrasound Elastography for Brain Diseases

Ultrasound (US) can be used to evaluate the brain structure and nervous system damage. Patients with neurologic symptoms need rapid, noninvasive imaging with high spatial resolution and tissue contrast. Magnetic resonance imaging is currently the most sensitive and specific imaging method for evaluating neuropathologic conditions. This approach does present some challenges, such as the need to transport patients who may be seriously ill to the magnetic resonance imaging suite and the need for patients to remain for a considerable time. Cranial US provides a very valuable imaging method for clinicians, which can make a rapid diagnosis and evaluation without ionizing radiation. The main disadvantage of cranial US is its low sensitivity and specificity for subtle/early lesions. In recent years, with the rapid development of anatomic and functional US technology, the practicability of US diagnosis and intervention has been greatly improved. Ultrasound elastography may have the potential to improve the sensitivity and specificity of various cranial nerve conditions. Ultrasound elastography has received considerable critical attention, and an increasing number of studies have recognized its critical role in evaluating brain diseases. At present, US elastography has been applied to the evaluation of traumatic brain injury, ischemic stroke, intraoperative brain tumors, and hypoxic ischemic encephalopathy. The latest animal experiments and human clinical trial developments in the applications of US elastography for brain diseases are summarized in this review.

A Spline Interpolation-based Data Reconstruction Technique for Estimation of Strain Time Constant in Ultrasound Poroelastography

Ultrasound poroelastography is a cost-effective and noninvasive imaging technique, which can be used to reconstruct mechanical parameters of tissues such as Young's modulus, Poisson's ratio, interstitial permeability, and vascular permeability. To estimate interstitial permeability and vascular permeability using poroelastography, accurate estimation of the strain time constant (TC) is required. This can be a challenging task due to the nonlinearity of the exponential strain curve and noise affecting the experimental data. Due to motion artifacts caused by the sonographer, animal/patient, and/or the environment, noise affecting some strain frames can be significantly higher than the strain signal. If these frames are used for the computation of the strain TC, the resulting TC estimate can be highly inaccurate, which, in turn, can cause high errors in the reconstructed mechanical parameters. In this paper, we introduce a cubic spline-based interpolation method, which allows to use only good quality strain frames (i.e., frames with sufficiently high signal-to-noise ratio [SNR]) to estimate the strain TC. Using finite element simulations, we demonstrate that the proposed interpolation method can improve the estimation accuracy of the strain TC by 46% with respect to the case where no interpolation and filtering are used and by 37% with respect to the case where the strain frames are Kalman filtered before TC estimation (at an SNR of 30 dB). We also prove the technical feasibility of the proposed technique using in vivo experimental data. While detecting the bad frames in both simulations and experiments, we assumed the lower limit SNR to be below 10 dB. Based on our results, the proposed technique may be of great help in applications relying on the accurate assessment of the temporal behavior of strain data.

Benign lesion evaluation: Factors causing the "stiff rim" sign in breast tissues using shear-wave elastography

OBJECTIVE:

To investigate the factors causing the "stiff rim" sign in breast lesions using shear-wave elastography.

METHODS:

A total of 907 patients with 907 lesions were included retrospectively in this study. Traditional ultrasound and shear-wave elastography imaging were both performed. Patients age, maximum diameter, depth, distance, echogenicity, shape, boundary, margin, internal components, CDFI, calcification, echogenicity attenuation and longitudinal growth of lesions were observed and calculated by both univariate and multivariate analyses.

RESULTS:

Univariate analyses indicated that the age, depth, shape, margin, internal components, CDFI, calcification and pathology showed significant difference between the benign lesions with and without a "stiff rim", whereas there was no correlation of "stiff rim" with maximum diameter, distance, boundary, echogenicity, echo attenuation and longitudinal growth of the lesions. Multivariate analysis expressed that CDFI, margin, internal components, depth and age were significantly associated with the "stiff rim" sign in breast benign lesions, whereas there was no correlation with the pathology, shape or calcification of the lesions.

CONCLUSIONS:

The "stiff rim" sign can be helpful for differentiation between benign and malignant lesions. Older patients with a "stiff rim" sign whose benign masses are deep, poorly defined, heterogeneous and have a positive CDFI should be examined more closely to avoid unnecessary false-positives.

ADVANCES IN KNOWLEDGE:

The "stiff rim" sign can be helpful for differentiation between benign and malignant lesions. Positive CDFI, poorly defined margin, heterogeneous internal components, deep depth and older age were significantly associated with the "stiff rim" sign in benign breast lesions.

Clinical application of sonoelastography in thyroid, prostate, kidney, pancreas, and deep venous thrombosis

This article reviews the clinical applications of current ultrasound elastography methods in non-hepatic conditions including thyroid nodules, prostate cancer, chronic kidney disease, solid renal lesions, pancreatic lesions, and deep vein thrombosis. Pathophysiology alters tissue mechanical properties via ultrastructural changes including fibrosis, increased cellularity, bleeding, and necrosis, creating a target biomarker, which can be imaged qualitatively or quantitatively with US elastography. US elastography methods can add information to conventional US methods and improve the diagnostic performance of conventional US in a range of disease processes.

Spatial and Temporal Control of Hyperthermia Using Real Time Ultrasonic Thermal Strain Imaging with Motion Compensation, Phantom Study

Mild hyperthermia has been successfully employed to induce reversible physiological changes that can directly treat cancer and enhance local drug delivery. In this approach, temperature monitoring is essential to avoid undesirable biological effects that result from thermal damage. For thermal therapies, Magnetic Resonance Imaging (MRI) has been employed to control real-time Focused Ultrasound (FUS) therapies. However, combined ultrasound imaging and therapy systems offer the benefits of simple, low-cost devices that can be broadly applied. To facilitate such technology, ultrasound thermometry has potential to reliably monitor temperature. Control of mild hyperthermia was previously achieved using a proportional-integral-derivative (PID) controller based on thermocouple measurements. Despite accurate temporal control of heating, this method is limited by the single position at which the temperature is measured. Ultrasound thermometry techniques based on exploiting the thermal dependence of acoustic parameters (such as longitudinal velocity) can be extended to create thermal maps and allow an accurate monitoring of temperature with good spatial resolution. However, in vivo applications of this technique have not been fully developed due to the high sensitivity to tissue motion. Here, we propose a motion compensation method based on the acquisition of multiple reference frames prior to treatment. The technique was tested in the presence of 2-D and 3-D physiological-scale motion and was found to provide effective real-time temperature monitoring. PID control of mild hyperthermia in presence of motion was then tested with ultrasound thermometry as feedback and temperature was maintained within 0.3°C of the requested value.

Spatial and Temporal Control of Hyperthermia Using Real Time Ultrasonic Thermal Strain Imaging with Motion Compensation, Phantom Study

Mild hyperthermia has been successfully employed to induce reversible physiological changes that can directly treat cancer and enhance local drug delivery. In this approach, temperature monitoring is essential to avoid undesirable biological effects that result from thermal damage. For thermal therapies, Magnetic Resonance Imaging (MRI) has been employed to control real-time Focused Ultrasound (FUS) therapies. However, combined ultrasound imaging and therapy systems offer the benefits of simple, low-cost devices that can be broadly applied. To facilitate such technology, ultrasound thermometry has potential to reliably monitor temperature. Control of mild hyperthermia was previously achieved using a proportional-integral-derivative (PID) controller based on thermocouple measurements. Despite accurate temporal control of heating, this method is limited by the single position at which the temperature is measured. Ultrasound thermometry techniques based on exploiting the thermal dependence of acoustic parameters (such as longitudinal velocity) can be extended to create thermal maps and allow an accurate monitoring of temperature with good spatial resolution. However, in vivo applications of this technique have not been fully developed due to the high sensitivity to tissue motion. Here, we propose a motion compensation method based on the acquisition of multiple reference frames prior to treatment. The technique was tested in the presence of 2-D and 3-D physiological-scale motion and was found to provide effective real-time temperature monitoring. PID control of mild hyperthermia in presence of motion was then tested with ultrasound thermometry as feedback and temperature was maintained within 0.3°C of the requested value.

Observation of a cavitation cloud in tissue using correlation between ultrafast ultrasound images

The local application of ultrasound is known to improve drug intake by tumors. Cavitating bubbles are one of the contributing effects. A setup in which two ultrasound transducers are placed confocally is used to generate cavitation in ex vivo tissue. As the transducers emit a series of short excitation bursts, the evolution of the cavitation activity is monitored using an ultrafast ultrasound imaging system. The frame rate of the system is several thousands of images per second, which provides several tens of images between consecutive excitation bursts. Using the correlation between consecutive images for speckle tracking, a decorrelation of the imaging signal appears due to the creation, fast movement, and dissolution of the bubbles in the cavitation cloud. By analyzing this area of decorrelation, the cavitation cloud can be localized and the spatial extent of the cavitation activity characterized.

A survey of ultrasound elastography approaches to percutaneous ablation monitoring

Percutaneous thermal ablation has been widely used as a minimally invasive treatment for tumors. Treatment monitoring is essential for preventing complications while ensuring treatment efficacy. Mechanical testing measurements on tissue reveal that tissue stiffness increases with temperature and ablation duration. Different types of imaging methods can be used to monitor ablation procedures, including temperature or thermal strain imaging, strain imaging, modulus imaging, and shear modulus imaging. Ultrasound elastography demonstrates the potential to become the primary imaging modality for monitoring percutaneous ablation. This review briefly presented the state-of-the-art ultrasound elastography approaches for monitoring radiofrequency ablation and microwave ablation. These techniques were divided into four groups: quasi-static elastography, acoustic radiation force elastography, sonoelastography, and applicator motion elastography. Their advantages and limitations were compared and discussed. Future developments were proposed with respect to heat-induced bubbles, tissue inhomogeneities, respiratory motion, three-dimensional monitoring, multi-parametric monitoring, real-time monitoring, experimental data center for percutaneous ablation, and microwave ablation monitoring.

Freehand elastography for determination of breast cancer size: comparison with B-mode sonography and histopathologic measurement

OBJECTIVES

Elastography assesses the strain of soft tissues and is used to enhance diagnostic accuracy in evaluating breast tumors, but minimal data exist on its ability to accurately assess tumor size. This study was performed to assess the preoperative accuracy of measuring the size of biopsyproven breast cancer lesions with elastography and conventional B-mode sonography compared with the reference standard size measured by histopathologic examination.

METHODS

Elastography and conventional B-mode sonography were performed on 69 women with histologically proven breast cancer, and tumor sizes on both modalities were recorded. These measurements were compared with the final pathologic size, which was used as the reference standard. The sizes and differences between sonographic, elastographic, and pathologic measurements were statistically tested, and an analysis of equivalence to the reference standard was performed using Bland-Altman plots.

RESULTS

There was a significant difference between sizes on elastography and pathologic examination, with elastography overestimating the tumor size (P = .0187). Sonography slightly underestimated the tumor size, but this finding was not significant (P = .36). Bland-Altman plots confirmed that sonography but not elastography was an acceptable standard compared with the pathologic size.

CONCLUSIONS

Breast elastography but not B-mode sonography overestimates the size of breast tumors compared with the final pathologic size.

Diagnostic performance of ultrasound and ultrasound elastography with respect to physician experience

The aim of this study was to compare the diagnostic performance of gray-scale ultrasound (US), elastography and a combination of gray-scale ultrasound and elastography (US-E) in differentiating benign and malignant thyroid nodules with respect to the level of physician experience. Three hundred fifty-eight patients with 367 thyroid nodules who underwent both gray-scale US and elastography, from November 2011 to January 2012, were included in this study. The diagnostic performance of US performed by experienced and less experienced physicians was compared. Comparisons of the diagnostic performance of US, elastography and US-E were evaluated for each group separately. Of 367 nodules, 121 were malignant and 246 were benign. When we compared the diagnostic performance of the experienced and less experienced physician groups, specificity was statistically higher in the experienced physician group for both US alone (p = 0.001) and US-E (p = 0.048). However, the experienced and less experienced physician groups did not differ significantly on other measures of diagnostic performance, regardless of modality. For the experienced physicians, the specificity and positive predictive value US were 88.0% and 76.8%, respectively; both of them were significantly higher than the corresponding values for US-E. For the less experienced physicians, specificity was significantly higher on elastography (93.8%) than on US (71.4%) (p < 0.001). However, diagnostic performance did not differ significantly between US and US-E for the less experienced physicians. Experienced physicians had superior specificity compared with less experienced physicians. The diagnostic performance of elastography and US-E was inferior compared with that of US alone, irrespective of the level of experience of the physician.

Averaging improves strain images of the biceps brachii using quasi-static ultrasound elastography

OBJECTIVE

Quasi-static ultrasound elastography is a technique for measuring tissue deformation (strain) under externally applied loading and can be used to identify the presence of abnormalities. The objective of this study was to demonstrate the efficacy of averaging strain images from repeated compression cycles in mitigating user-induced error using quasi-static ultrasound elastography.

METHODS

Freehand compressions were performed with an ultrasound transducer on the biceps brachii of nine participants (five males and four females), as well as with a custom automated compression system. Sets of strain images from the freehand techniques were averaged to create single representative images and compared against strain images from the automated compressions using both qualitative and quantitative metrics.

RESULTS

Significant improvements in intra-operator repeatability and interoperator reproducibility can be achieved by averaging strain images from four to eight repeated compressions. The resulting strain images did not lose significant image data compared with strain images from single automated compressions.

CONCLUSION

Averaging is introduced as a feasible and appropriate technique to improve strain image quality without sacrificing important image data.

ADVANCES IN KNOWLEDGE

Simple averaging of multiple freehand elastography measures can achieve a similar degree of accuracy, repeatability and reproducibility as that of more awkward and expensive automated methods. The resulting elastograms can be used to obtain a more accurate and complete diagnosis without additional cost to the doctor or the patient.

Real-time monitoring of high-intensity focused ultrasound treatment using axial strain and axial-shear strain elastograms

Axial strain elastograms (ASEs) have been found to help visualize sonographically invisible thermal lesions. However, in most studies involving high-intensity focused ultrasound (HIFU)-induced thermal lesions, elastography imaging was performed separately later, after the lesion was formed. In this article, the feasibility of monitoring, in real time, tissue elasticity variation during HIFU treatment and immediately thereafter is explored using quasi-static elastography. Further, in addition to ASEs, we also explore the use of simultaneously acquired axial-shear strain elastograms (ASSEs) for HIFU lesion visualization. Experiments were performed on commercial porcine liver samples in vitro. The HIFU experiments were conducted at two applied acoustic power settings, 35 and 20 W. The experimental setup allowed us to interrupt the HIFU pulse momentarily several different times during treatment to perform elastographic compression and data acquisition. At the end of the experiments, the samples were cut along the imaging plane and photographed to compare size and location of the formed lesion with those visualized on ASEs and ASSEs. Single-lesion and multiple-lesion experiments were performed to assess the contribution of ASEs and ASSEs to lesion visualization and treatment monitoring tasks. At both power settings, ASEs and ASSEs provided accurate location information during HIFU treatment. At the low-power setting case, ASEs and ASSEs provide accurate lesion size in real-time monitoring. Lesion appearance in ASEs and ASSEs was affected by the cavitation bubbles produced at the high-power setting. The results further indicate that the cavitation bubbles influence lesion appearance more in ASEs than in ASSEs. Both ASEs and ASSEs provided accurate size information after a waiting period that allowed the cavitation bubbles to disappear. The results indicate that ASSEs not only improve lesion visualization and size measurement of a single lesion, but, under certain conditions, also help to identify untreated gaps between adjacent lesions with high contrast.

In vivo thermal ablation monitoring using ultrasound echo decorrelation imaging

Previous work indicated that ultrasound echo decorrelation imaging can track and quantify changes in echo signals to predict thermal damage during in vitro radiofrequency ablation (RFA). In the in vivo studies reported here, the feasibility of using echo decorrelation imaging as a treatment monitoring tool was assessed. RFA was performed on normal swine liver (N = 5), and ultrasound ablation using image-ablate arrays was performed on rabbit liver implanted with VX2 tumors (N = 2). Echo decorrelation and integrated backscatter were computed from Hilbert transformed pulse-echo data acquired during RFA and ultrasound ablation treatments. Receiver operating characteristic (ROC) curves were employed to assess the ability of echo decorrelation imaging and integrated backscatter to predict ablation. Area under the ROC curves (AUROC) was determined for RFA and ultrasound ablation using echo decorrelation imaging. Ablation was predicted more accurately using echo decorrelation imaging (AUROC = 0.832 and 0.776 for RFA and ultrasound ablation, respectively) than using integrated backscatter (AUROC = 0.734 and 0.494).

A freehand ultrasound elastography system with tracking for in vivo applications

Ultrasound transducers are commonly tracked in modern ultrasound navigation/guidance systems. In this article, we demonstrate the advantages of incorporating tracking information into ultrasound elastography for clinical applications. First, we address a common limitation of freehand palpation: speckle decorrelation due to out-of-plane probe motion. We show that by automatically selecting pairs of radio frequency frames with minimal lateral and out-of-plane motions, combined with a fast and robust displacement estimation technique, greatly improves in vivo elastography results. We also use tracking information and image-quality measures to fuse multiple images with similar strains that are taken from roughly the same location so as to obtain a high-quality elastography image. Finally, we show that tracking information can be used to give the user partial control over the rate of compression. Our methods are tested on a tissue-mimicking phantom, and experiments have been conducted on intraoperative data acquired during animal and human experiments involving liver ablation. Our results suggest that in challenging clinical conditions, our proposed method produces reliable strain images and eliminates the need for a manual search through the ultrasound data in order to find radio frequency pairs suitable for elastography.

Ultrasonic imaging of 3D displacement vectors using a simulated 2D array and beamsteering

Most quasi-static ultrasound elastography methods image only the axial strain, derived from displacements measured in the direction of ultrasound propagation. In other directions, the beam lacks high resolution phase information and displacement estimation is therefore less precise. However, these estimates can be improved by steering the ultrasound beam through multiple angles and combining displacements measured along the different beam directions. Previously, beamsteering has only considered the 2D case to improve the lateral displacement estimates. In this paper, we extend this to 3D using a simulated 2D array to steer both laterally and elevationally in order to estimate the full 3D displacement vector over a volume. The method is tested on simulated and phantom data using a simulated 6-10MHz array, and the precision of displacement estimation is measured with and without beamsteering. In simulations, we found a statistically significant improvement in the precision of lateral and elevational displacement estimates: lateral precision 35.69μm unsteered, 3.70μm steered; elevational precision 38.67μm unsteered, 3.64μm steered. Similar results were found in the phantom data: lateral precision 26.51μm unsteered, 5.78μm steered; elevational precision 28.92μm unsteered, 11.87μm steered. We conclude that volumetric 3D beamsteering improves the precision of lateral and elevational displacement estimates.

On the advantages of imaging the axial-shear strain component of the total shear strain in breast tumors

Axial-shear strain elastography was described recently as a method to visualize the state of bonding at an inclusion boundary. Although total shear strain elastography was initially proposed for this purpose, it did not evolve beyond the initial reported finite element model (FEM) and simulation studies. One of the major reasons for this was the practical limitation in estimating the tissue motion perpendicular (lateral) to the ultrasound (US) beam as accurately as the motion along the US beam (axial). Nevertheless, there has been a sustained effort in developing methods to improve the lateral motion tracking accuracy and thereby obtain better quality total shear strain elastogram (TSSE). We hypothesize that in some cases, even if good quality TSSE becomes possible, it may still be advantageous to utilize only the axial-shear strain (one of the components of the total shear strain) elastogram (ASSE). Specifically, we show through FEM and corroborating tissue-mimicking gelatin phantom experiments that the unique "fill-in" discriminant feature that was introduced recently for asymmetric breast lesion classification is depicted only in the ASSE and not in the TSSE. Note that the presence or conspicuous absence of this feature in ASSE was shown to characterize asymmetric inclusions' boundaries as either loosely-bonded or firmly-bonded to the surrounding, respectively. This might be an important observation because the literature suggests that benign breast lesions tend to be loosely-bonded, while malignant tumors are usually firmly-bonded. The results from the current study demonstrate that the use of shear strain lesion "fill-in" as a discriminant feature in the differentiation between asymmetric malignant and benign breast lesions is only possible when using the ASSEs and not the TSSEs.

A pilot study to evaluate assisted freehand ultrasound elasticity imaging in the sizing of early breast cancer: a comparison of B-mode and AFUSON elasticity ultrasound with histopathology measurements

OBJECTIVE

This pilot study investigates the role of assisted-freehand ultrasound (AFUSON) elasticity imaging of the breast in assessing the contour, size and area of 23 early breast cancers by making comparison of AFUSON with the equivalent B-mode ultrasound images and gold standard histopathology slides.

METHODS

The B-mode, AFUSON and digitised histopathology slides of three early breast cancers were compared for contour, size and area with histopathology scans. AFUSON features that corresponded to areas of known malignant change on the histopathology slides were regarded as diagnostic. These diagnostic criteria were then applied to the B-mode and AFUSON elasticity images of all 23 breast cancers in the pilot study without having the availability of the histopathology scans for reference. Corresponding diameters were measured and the results were compared with the equivalent measurements on the scans of the histology slides. The results were tabulated in histogram form. Diagnostic confidence levels were evaluated.

RESULTS

Size dimension accuracy increased from 66% using B-mode alone to 82% using combined B-mode and AFUSON elasticity images. Tumour area accuracy was also increased. A small number of cases had a striking visual similarity of shape on AFUSON elasticity scans and histopathology slides.

CONCLUSION

In spite of the shortfalls in this study, AFUSON elasticity imaging was capable of acquiring some high-quality images that showed strong correlation between AFUSON elasticity and scans of histology slides. Further studies will be carried out to refine the technique and determine if it has a role in the diagnosis and management of breast cancer.

Dynamic frame selection for in vivo ultrasound temperature estimation during radiofrequency ablation

Minimally invasive therapies such as radiofrequency ablation have been developed to treat cancers of the liver, prostate and kidney without invasive surgery. Prior work has demonstrated that ultrasound echo shifts due to temperature changes can be utilized to track the temperature distribution in real time. In this paper, a motion compensation algorithm is evaluated to reduce the impact of cardiac and respiratory motion on ultrasound-based temperature tracking methods. The algorithm dynamically selects the next suitable frame given a start frame (selected during the exhale or expiration phase where extraneous motion is reduced), enabling optimization of the computational time in addition to reducing displacement noise artifacts incurred with the estimation of smaller frame-to-frame displacements at the full frame rate. A region of interest that does not undergo ablation is selected in the first frame and the algorithm searches through subsequent frames to find a similarly located region of interest in subsequent frames, with a high value of the mean normalized cross-correlation coefficient value. In conjunction with dynamic frame selection, two different two-dimensional displacement estimation algorithms namely a block matching and multilevel cross-correlation are compared. The multi-level cross-correlation method incorporates tracking of the lateral tissue expansion in addition to the axial deformation to improve the estimation performance. Our results demonstrate the ability of the proposed motion compensation using dynamic frame selection in conjunction with the two-dimensional multilevel cross-correlation to track the temperature distribution.

Intra-operative ultrasound hand-held strain imaging for the visualization of ablations produced in the liver with a toroidal HIFU transducer: first in vivo results

The use of hand-held ultrasound strain imaging for the intra-operative real-time visualization of HIFU (high-intensity focused ultrasound) ablations produced in the liver by a toroidal transducer was investigated. A linear 12 MHz ultrasound imaging probe was used to obtain radiofrequency signals. Using a fast cross-correlation algorithm, strain images were calculated and displayed at 60 frames s(-1), allowing the use of hand-held strain imaging intra-operatively. Fourteen HIFU lesions were produced in four pigs. Intra-operative strain imaging of HIFU ablations in the liver was feasible owing to the high frame rate. The correlation between dimensions measured on gross pathology and dimensions measured on B-mode images and on strain images were R = 0.72 and R = 0.94 respectively. The contrast between ablated and non-ablated tissue was significantly higher (p < 0.05) in the strain images (22 dB) than in the B-mode images (9 dB). Strain images allowed equivalent or improved definition of ablated regions when compared with B-mode images. Real-time intra-operative hand-held strain imaging seems to be a promising complement to conventional B-mode imaging for the guidance of HIFU ablations produced in the liver during an open procedure. These results support that hand-held strain imaging outperforms conventional B-mode ultrasound and could potentially be used for the assessment of thermal therapies.

Axial-shear strain distributions in an elliptical inclusion model: experimental validation and in vivo examples with implications to breast tumor classification

Recently, we reported on the axial-shear strain fill-in of the interior of loosely bonded stiff elliptical inclusions in a soft background at non-normal orientations, and the lack of fill-in in firmly bonded inclusions at any orientation. In this paper, we report on the experimental validation of the simulation studies using tissue-mimicking gelatin-based phantoms. We also show a few confirmatory examples of the existence of these phenomena in benign vs. malignant breast lesions in vivo. Phantom experiments showed that axial-shear strain zones caused by firmly bonded elliptical inclusions occurred only outside of the inclusion, as predicted by the simulation. By contrast, the axial-shear strain zones filled in the interior of loosely bonded elliptical inclusions at non-normal orientations. The axial-shear strain elastograms obtained from the in vivo cases appeared to be in general agreement with our experimental results. The results reported in this paper may have important clinical implications. Specifically, axial-shear strain fill-in inside an inclusion may be a unique signature of stiff, loosely bonded, ellipsoidal or elongated inclusions at non-normal orientations. Thus, it may be useful as a marker of benignity of benign breast lesions (e.g., fibroadenomas) that are generally stiff, elongated and loosely bonded to the host tissues.

Assisted-freehand ultrasound elasticity imaging

Good-quality elasticity imaging requires highly controlled compressions of the breast, which are often challenging to obtain with freehand, even by an experienced radiologist. This paper presents assisted-freehand ultrasound (AFUSON): a fusion of freehand and automated ultrasound systems designed to assisted elasticity imaging acquisition while remaining as flexible as freehand. In the form of a hand-held device, this semi-automatic solution delivers both increased acquisition precision and control. Compared with freehand acquisitions, it reduces out-of-plane motion decorrelation by one-half and lateral motion by one-third, increases within-scan repeatability by 50%, and does so across operators.

[Elastosonography of thyroid lesions]

While ultrasound is the imaging modality of choice for diagnosis of thyroid lesions, characterization remains limited and tissue diagnosis frequently is required for management. The availability of additional tools such as elastography may improve lesion characterization and direct management.

MATERIALS AND METHODS

A total of 96 patients (11 males and 85 females; 58+/-24 years) referred for fine needle aspiration (FNA) of mainly solid thyroid nodules 9-32 mm in diameter underwent conventional US and elastosonography. Results on elastography were correlated with histological results from FNA and classified as follows: suspected malignant lesion, suspected benign lesion, suspicious, indeterminate.

RESULTS

The nodules were classified as follows: 95 nodules were soft (classes I and II) and 13 nodules were hard (classes III and IV). No cancers were detected in class and II lesions and 6 cancers were detected in class III and IV lesions. FNA provided insufficient cellular material for diagnosis in 5 class I-II nodules and 2 class III-IV nodules.

CONCLUSION

Real-time elastosonography may be a useful adjunct to conventional US in the evaluation and characterization of thyroid nodules allowing identification of patients at high risk of malignancy for whom tissue diagnosis and/or close follow-up is required.

A novel motion compensation algorithm for acoustic radiation force elastography

A novel method of physiological motion compensation for use with radiation force elasticity imaging has been developed. The method utilizes a priori information from finite element method models of the response of soft tissue to impulsive radiation force to isolate physiological motion artifacts from radiation force-induced displacement fields. The new algorithmis evaluated in a series of clinically realistic imaging scenarios, and its performance is compared to that achieved with previously described motion compensation algorithms. Though not without limitations, the new model-based motion compensation algorithm performs favorably in many circumstances and may be a logical choice for use with in vivo abdominal imaging.

The feasibility of using poroelastographic techniques for distinguishing between normal and lymphedematous tissues in vivo

Lymphedema is a common condition involving an abnormal accumulation of lymphatic fluid in the interstitial space that causes swelling, most often in the arm(s) and leg(s). Lymphedema is a significant lifelong concern that can be congenital or develop following cancer treatment or cancer metastasis. Common methods of evaluation of lymphedema are mostly qualitative making it difficult to reliably assess the severity of the disease, a key factor in choosing the appropriate treatment. In this paper, we investigate the feasibility of using novel elastographic techniques to differentiate between lymphedematous and normal tissues. This study represents the first step of a larger study aimed at investigating the combined use of elastographic and sonographic techniques for the detection and staging of lymphedema. In this preliminary study, poroelastographic images were generated from the leg (8) and arm (4) subcutis of five normal volunteers and seven volunteers having lymphedema, and the results were compared using statistical analyses. The preliminary results reported in this paper suggest that it may be feasible to perform poroelastography in different lymphedematous tissues in vivo and that poroelastography techniques may be of help in differentiating between normal and lymphedematous tissues.

Myocardial elastography at both high temporal and spatial resolution for the detection of infarcts

Myocardial elastography is a novel method for noninvasively assessing regional myocardial function, with the advantages of high spatial and temporal resolution and high signal-to-noise ratio (SNR). In this paper, in-vivo experiments were performed in anesthetized normal and infarcted mice (one day after left anterior descending coronary artery [LAD] ligation) using a high-resolution (30 MHz) ultrasound system (Vevo 770, VisualSonics Inc., Toronto, ON, Canada). Radiofrequency (RF) signals of the left ventricle (LV) in longitudinal (long-axis) view and the associated electrocardiogram (ECG) were simultaneously acquired. Using a retrospective ECG gating technique, 2-D full field-of-view RF frames were acquired at an extremely high frame rate (8 kHz) that resulted in high-quality incremental displacement and strain estimation of the myocardium. The incremental results were further accumulated to obtain the cumulative displacements and strains. Two-dimensional and M-mode displacement images and strain images (elastograms), as well as displacement and strain profiles as a function of time, were compared between normal and infarcted mice. Incremental results clearly depicted cardiac events including LV contraction, LV relaxation and isovolumetric phases in both normal and infarcted mice, and also evidently indicated reduced motion and deformation in the infarcted myocardium. The elastograms indicated that the infarcted regions underwent thinning during systole rather than thickening, as in the normal case. The cumulative elastograms were found to have higher elastographic SNR (SNR(e)) than the incremental elastograms (e.g., 10.6 vs. 4.7 in a normal myocardium, and 6.0 vs. 2.4 in an infarcted myocardium). Finally, preliminary statistical results from nine normal (m = 9) and seven infarcted (n = 7) mice indicated the capability of the cumulative strain in differentiating infracted from normal myocardia. In conclusion, myocardial elastography could provide regional strain information at simultaneously high temporal (>/=0.125 ms) and spatial ( approximately 55 microm) resolution as well as high precision ( approximately 0.05 microm displacement). This technique was thus capable of accurately characterizing normal myocardial function throughout an entire cardiac cycle, at the same high resolution, and detecting and localizing myocardial infarction in vivo.

The impact of physiological motion on tissue tracking during radiation force imaging

The effect of physiological motion on the quality of radiation force elasticity images has been investigated. Experimental studies and simulated images were used to investigate the impact of motion effects on image quality metrics over a range of clinically realistic velocity and acceleration magnitudes. Evaluation criteria included motion filter effectiveness, image signal-to-noise ratio (SNR) and the contrast-to-noise ratio (CNR) of a stiff inclusion embedded in a homogeneous background material. Two transmit frequencies (2.5 and 4.4 MHz) were analyzed and contrasted in terms of image quality over a range of target motions. Results indicate that situations may exist where liver and cardiac motion magnitudes lead to poor image quality, but optimized transducer orientations may help suppress motion artifacts if some a priori information concerning target motion characteristics is known. In the presence of significant target motion, utilizing a lower transmit frequency can improve SNR and CNR in elasticity images.

Assessing image quality in effective Poisson's ratio elastography and poroelastography: II

Poroelastography is a novel elastographic technique for imaging the time variation of the mechanical behaviour of poroelastic materials. Poroelastograms are generated as a series of time-sequenced effective Poisson's ratio (EPR) elastograms, obtained from the imaged material under sustained compression. In the companion report (Righetti et al 2007 Phys. Med. Biol. 52 1303), we investigated image quality of EPR elastography starting from a theoretical analysis of the performance limitations of axial strain elastography and lateral strain elastography. In this report, we extend this analysis to poroelastography. The theoretical analysis reported in these two companion papers allows understanding the performance limitations of these novel techniques and identifying the fundamental parameters that control their signal-to-noise ratio, contrast-to-noise ratio and resolution. The results of these studies also indicate that EPR elastograms and poroelastograms of reasonable image quality can be generated in practical applications that may be of clinical interest provided that advanced elastographic techniques in combination with other commonly employed imaging methods to increase signal-to-noise and contrast-to-noise ratios are used.

Signal-to-noise ratio, contrast-to-noise ratio and their trade-offs with resolution in axial-shear strain elastography

In axial-shear strain elastography, the local axial-shear strain resulting from the application of quasi-static axial compression to an inhomogeneous material is imaged. In this paper, we investigated the image quality of the axial-shear strain estimates in terms of the signal-to-noise ratio (SNR(asse)) and contrast-to-noise ratio (CNR(asse)) using simulations and experiments. Specifically, we investigated the influence of the system parameters (beamwidth, transducer element pitch and bandwidth), signal processing parameters (correlation window length and axial window shift) and mechanical parameters (Young's modulus contrast, applied axial strain) on the SNR(asse) and CNR(asse). The results of the study show that the CNR(asse) (SNR(asse)) is maximum for axial-shear strain values in the range of 0.005-0.03. For the inclusion/background modulus contrast range considered in this study (<10), the CNR(asse) (SNR(asse)) is maximum for applied axial compressive strain values in the range of 0.005%-0.03%. This suggests that the RF data acquired during axial elastography can be used to obtain axial-shear strain elastograms, since this range is typically used in axial elastography as well. The CNR(asse) (SNR(asse)) remains almost constant with an increase in the beamwidth while it increases as the pitch increases. As expected, the axial shift had only a weak influence on the CNR(asse) (SNR(asse)) of the axial-shear strain estimates. We observed that the differential estimates of the axial-shear strain involve a trade-off between the CNR(asse) (SNR(asse)) and the spatial resolution only with respect to pitch and not with respect to signal processing parameters. Simulation studies were performed to confirm such an observation. The results demonstrate a trade-off between CNR(asse) and the resolution with respect to pitch.