PubMed Health⌕ Search

Biomedical subjects

Jonathan Ophir

Publications and source records attributed to Jonathan Ophir.

At least 19 recordsLinked to original sources

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.

Acoustics↗

Resolution of axial shear strain elastography.

The technique of mapping the local axial component of the shear strain due to quasi-static axial compression is defined as axial shear strain elastography. In this paper, the spatial resolution of axial shear strain elastography is investigated through simulations, using an elastically stiff cylindrical lesion embedded in a homogeneously softer background. Resolution was defined as the smallest size of the inclusion for which the strain value at the inclusion/background interface was greater than the average of the axial shear strain values at the interface and inside the inclusion. The resolution was measured from the axial shear strain profile oriented at 45 degrees to the axis of beam propagation, due to the absence of axial shear strain along the normal directions. The effects of the ultrasound system parameters such as bandwidth, beamwidth and transducer element pitch along with signal processing parameters such as correlation window length (W) and axial shift (DeltaW) on the estimated resolution were investigated. The results show that the resolution (at 45 degrees orientation) is determined by the bandwidth and the beamwidth. However, the upper bound on the resolution is limited by the larger of the beamwidth and the window length, which is scaled inversely to the bandwidth. The results also show that the resolution is proportional to the pitch and not significantly affected by the axial window shift.

Algorithms↗

Analysis of a hybrid spectral strain estimation technique in elastography.

Conventional spectral elastographic techniques estimate strain using cross-correlation methods. Despite promising results, decorrelation effects compromise the accuracy of these techniques and, subsequently, the tissue strain estimates. Since tissue compression in the time-domain corresponds to upscaling in the frequency-domain, decorrelation effects become more pronounced as tissue strains increase and are a fundamental concern in spectral cross-correlation elastography. In this paper, a two-stage hybrid spectral elastographic technique is introduced. For the first stage, an approximated spectral scaling factor (i.e. initial strain estimate) is employed to compensate for bandwidth broadening (due to tissue compression) between pre- and post-compression power spectra pairs. The second stage then estimates any residual strain information using spectral cross-correlation methods due to improper scaling factor selection in the first stage. This novel hybrid spectral elastographic technique was compared to both conventional spectral and adaptive temporal elastographic methods in simulation and experimentation. In addition to demonstrating enhancement in performance over the conventional spectral elastographic technique, the hybrid spectral-based method introduced in this paper is shown to outperform the adaptive temporal-based elastographic approach.

Algorithms↗

Enhancing the performance of model-based elastography by incorporating additional a priori information in the modulus image reconstruction process.

Model-based elastography is fraught with problems owing to the ill-posed nature of the inverse elasticity problem. To overcome this limitation, we have recently developed a novel inversion scheme that incorporates a priori information concerning the mechanical properties of the underlying tissue structures, and the variance incurred during displacement estimation in the modulus image reconstruction process. The information was procured by employing standard strain imaging methodology, and introduced in the reconstruction process through the generalized Tikhonov approach. In this paper, we report the results of experiments conducted on gelatin phantoms to evaluate the performance of modulus elastograms computed with the generalized Tikhonov (GTK) estimation criterion relative to those computed by employing the un-weighted least-squares estimation criterion, the weighted least-squares estimation criterion and the standard Tikhonov method (i.e., the generalized Tikhonov method with no modulus prior). The results indicate that modulus elastograms computed with the generalized Tikhonov approach had superior elastographic contrast discrimination and contrast recovery. In addition, image reconstruction was more resilient to structural decorrelation noise when additional constraints were imposed on the reconstruction process through the GTK method.

Algorithms↗

Comparison of shift estimation strategies in spectral elastography.

This paper compares the performance of various spectral shift estimators for use in spectral elastography, namely, the normalized cross-correlation (NCC), sum squared difference (SSD) and sum absolute difference (SAD). Simulation and experimental results demonstrate that the spectral SSD-based elastographic method exhibits no marked difference in performance compared to the more computationally costly NCC-based approach, which has conventionally been the preferred estimator in spectral elastography. The spectral SAD-based strain estimator, despite being computationally less burdening, failed to exhibit performance comparable to that of the NCC- and SSD-based techniques. Furthermore, though spectral subsample estimation techniques using a cosine-fit interpolation method outperformed that of the parabolic-fit method in terms of both reduced bias errors and standard deviations, the latter was analyzed in this study due to computational simplicity. The role of spectral density was evaluated without and with parabolic-based subsample interpolation. Based on minimizing computational complexity, it is concluded that a (low density) spectral SSD strain estimator coupled with parabolic-based subsample estimation is the preferred choice for spectral elastography.

Bias↗

In vivo elastographic investigation of ethanol-induced hepatic lesions.

Ethanol-induced hepatic lesions were investigated in swine for in vivo use as a strain imaging animal model. Lesions (n = 25) were induced by injecting ethanol (doses 0.33 to 2.0 mL) directly into the surgically exposed liver at depths of 12, 15 or 25 mm. Lesions were imaged with a modified HDI 1000 scanner (Philips Medical Systems, Bothell, WA, USA). The elastograms (n = 91) characterized lesions as being areas harder than the surrounding soft hepatic tissue. Elastographic lesion sizes and the corresponding injected ethanol dose used to induce the lesions were shown to be statistically significant (r(2) = 0.22; p = 0.029) using a linear regression analysis. Additionally, lesion depth was shown to be statistically insignificant (r(2) < 0.12; p > 0.10) when regressed against elastographic lesion size. An analysis of elastographic and gross pathology lesion sizes indicated no correlation (r(2) < 0.01; p = 0.973). Subsequently, lesion types were sorted by size and regression lines were computed from quasilinear regions of the corresponding run charts. Trend lines indicate a four-to-three size relationship between the selected elastographic and pathology lesion sizes. Comparison of elastogram lesion sizes from two independent observers using a paired t-test resulted in no statistically significant difference (p = 0.14). In conclusion, ethanol-induced hepatic lesions in swine is a suitable animal model for evaluation of strain-based imaging systems, due to the ease of generation and repeatability.

Animals↗

A method for generating permeability elastograms and Poisson's ratio time-constant elastograms.

The feasibility of imaging the permeability and Poisson's ratio time-constant of porous media was investigated. The study involved the following steps. First, poroelastograms were generated from porous tofu phantoms under sustained compression. The sample materials used for the experiments were previously characterized through independent mechanical measurements. Second, corresponding Poisson's ratio time-constant elastograms were generated by calculating and displaying the decay time-constants of the local Poisson's ratios over the time interval in which the poroelastograms were acquired. Finally, for homogeneous samples, permeability elastograms were generated using the poroelastograms in combination with a previously proposed biphasic theoretical model, under very specific conditions. A comparison between the results obtained using poroelastography and those obtained through independent mechanical measurements suggests that poroelastography may be used for imaging the local time-dependent behavior of poroelastic media.

Computer Simulation↗

Comparative evaluation of strain-based and model-based modulus elastography.

Elastography based on strain imaging currently endures mechanical artefacts and limited contrast transfer efficiency. Solving the inverse elasticity problem (IEP) should obviate these difficulties; however, this approach to elastography is often fraught with problems because of the ill-posed nature of the IEP. The aim of the present study was to determine how the quality of modulus elastograms computed by solving the IEP compared with those produced using standard strain imaging methodology. Strain-based modulus elastograms (i.e., modulus elastograms computed by simply inverting strain elastograms based on the assumption of stress uniformity) and model-based modulus elastograms (i.e., modulus elastograms computed by solving the IEP) were computed from a common cohort of simulated and gelatin-based phantoms that contained inclusions of varying size and modulus contrast. The ensuing elastograms were evaluated by employing the contrast-to-noise ratio (CNR(e)) and the contrast transfer efficiency (CTE(e)) performance metrics. The results demonstrated that, at a fixed spatial resolution, the CNR(e) of strain-based modulus elastograms was statistically equivalent to those computed by solving the IEP. At low modulus contrast, the CTE(e) of both elastographic imaging approaches was comparable; however, at high modulus, the CTE(e) of model-based modulus elastograms was superior.

Algorithms↗

Investigation of parametric spectral estimation techniques for elasticity imaging.

Several autoregressive (AR) and autoregressive moving average (ARMA) parametric spectral estimators were evaluated for use in tissue strain estimation. Using both 1-D simulations and in vitro phantom experiments, the performance of these parametric spectral strain estimators were compared against both a nonparametric discrete Fourier transform (DFT) spectral strain estimator and a coherent elastographic technique. Parametric spectral estimator model orders were selected based on a modified strain filter approach. This technique illustrated the trade-offs between different signal-processing parameters and a strain estimator performance measure, namely the area under the strain filter (using applied strain dynamic range of 0.1 to 50%). The Yule-Walker AR spectral strain estimator outperformed all other parametric methods evaluated, but failed to outperform the DFT-based approach. Furthermore, both these spectral strain-estimation techniques exhibit an elastographic signal-to-noise ratio (SNR(e)) and strain estimation dynamic range not achievable using conventional elastography without global stretching.

Algorithms↗

Noise performance and signal-to-noise ratio of shear strain elastograms.

In this paper, we develop a theoretical expression for the signal-to-noise ratio (SNR) of shear strain elastograms. The previously-developed ideas for the axial strain filter (ASF) and lateral strain filter (LSF) are extended to define the concept of the shear strain filter (SSF). Some of our theoretical results are verified using simulations and phantom experiments. The results indicate that the signal-to-noise ratio of shear-strain elastograms (SNRsse) improves with increasing shear strain and with improvements in system parameters such as the sonographic signal-to-noise ratio (SNRs) beamwidth, center frequency and fractional bandwidth. The results also indicate that the amount of axial strain present along with the shear strain is an important parameter that determines the upper bound on SNRsse. The SNRsse will be higher in the absence of additional deformation due to axial strain.

Computer Simulation↗

A new method for generating poroelastograms in noisy environments.

Poroelastography has been recently introduced as a new elastographic technique that may be used to describe the spatial and temporal behavior of poroelastic materials. The experimental methodology proposed thus far for phantoms and tissues in vitro requires the acquisition of a precompression rf frame, the application of a unit step strain compression to the sample and the acquisition of subsequent post-compression frames from the material. Elastograms and poroelastograms are generated by cross-correlating the sequentially-acquired postcompression frames with the reference precompression frame. The application of poroelastography to tissues in vivo must address the echo decorrelation problems that are encountered due to uncontrolled tissue motion, which may become significant shortly after the acquisition of the precompression frame. In this paper, we investigate the feasibility of performing poroelastography experiments using an alternative experimental scheme. In the proposed experimental methodology, the reference precompression frame is continuously moved while the time interval between the frames that are correlated is kept short. This allows long data acquisition times with simultaneous minimization of the decorrelation due to undesired tissue motion in vivo. We validated this new method using both a step and a ramp compression functions. We performed poroelastographic simulations and experiments in phantoms and in tissues in vivo. The results were compared to those obtained using the traditional acquisition methodology. This study shows that the two methods yield similar results in vitro and suggests that the new method may be more robust to decorrelation noise in applications in vivo.

Arm↗

Comparing elastographic strain images with modulus images obtained using nanoindentation: preliminary results using phantoms and tissue samples.

Conventional elastography involves quasistatic mechanical compression (external or internal) of the tissue under ultrasonic insonification to obtain radiofrequency (RF) A-lines before and after compression. Cross-correlation of the pre- and postcompression A-lines results in displacement images with axial gradients that produce the strain images (strain elastograms). Though the strain elastograms show structural similarities to the modulus images, they are not related in a simple way to the modulus images because the strains depend on both modulus and geometry of the materials being deformed. Therefore, a quantification of the similarities between the strain and modulus images may enhance the interpretation confidence of strain elastograms in depicting tissue structure. To demonstrate similarities between modulus images and strain elastograms, a feasibility study of using nanoindentation to obtain modulus images of thin slices of tissue and tissue-mimicking phantoms (agar-gelatin mixtures) was performed first, with encouraging results. This was followed by a comparison of modulus images and strain elastograms obtained from the same sample slices. The experimental results indicated that, under certain experimental conditions, it is feasible to perform quantitative comparisons between strain images (using elastography) and modulus images. A good visual, as well as quantitative, correspondence between structures in the modulus and strain images could be obtained at a 3-mm scale.

Agar↗

The feasibility of using elastography for imaging the Poisson's ratio in porous media.

The feasibility of using elastography for experimentally estimating and imaging the Poisson's ratio of porous media under drained and undrained conditions was investigated. Using standard elastographic procedures, static and time-sequenced poroelastograms (strain ratio images) of homogeneous cylindrical gelatin and commercially available tofu samples were generated under sustained applied axial strain. The experimental data show similar trends to those that were observed in finite-elements simulations, and to those that were calculated from classical theoretical models proposed for biphasic materials with similar mechanical properties. To demonstrate the applicability of elastography to monitor time-dependent changes in nonhomogeneous porous structures as well, preliminary time-sequenced poroelastograms were obtained from two-layer porous phantoms and porcine muscle samples in vitro. The results suggest that elastography may have significant potential for quantitatively mapping the time-dependent mechanical behavior of poroelastic media, which is related to the dynamics of fluid flow and to the elasticity and permeability parameters of the media.

Animals↗

Tissue displacements during acupuncture using ultrasound elastography techniques.

Acupuncture needle manipulation has been previously shown to result in measurable changes in connective tissue architecture in animal experiments. In this study, we used a novel in vivo ultrasound (US)-based technique to quantify tissue displacement during acupuncture manipulation in humans. B-scan ultrasonic imaging was performed on the thighs of 12 human subjects at different stages of needle motion, including varying amounts of rotation, downward and upward movement performed with a computer-controlled acupuncture needling instrument. Tissue displacements, estimated using cross-correlation techniques, provided successful mapping and quantitative analysis of spatial and temporal tissue behavior during acupuncture needle manipulation. Increasing amounts of rotation had a significant linear effect on tissue displacement during downward and upward needle motion, as well as on rebound tissue displacement after downward needle movement. In addition to being a valuable tool for studies of acupuncture's mechanism of action, this technique may have applications to other types of needling including biopsies.

Acupuncture Therapy↗

Dynamic measurement of internal solid displacement in articular cartilage using ultrasound backscatter.

Mechanics of articular cartilage can be represented using poroelastic theories where fluid and solid displacements are viscously coupled to create a time-dependent spatially heterogeneous behavior. In recent models of this tissue, finite element methods have been used to predict tissue deformation as a function of time for adult articular cartilage bearing a characteristic depth-dependent structure and composition. However, current experimental methods are limited in providing verification of these predictions. The current study presents an apparatus for imaging the radial displacement profile of cartilage in unconfined compression using an ultrasound technique called elastography. We acquired ultrasound A-scans across the lateral diameter of full-thickness cartilage disks containing a thin layer of underlying bone, during axial compression. Elastography was then applied to correlate temporally sequential A-scans to estimate the solid radial displacement profile in articular cartilage while it undergoes compression and stress-relaxation. Both time-dependent and depth-dependent solid radial displacement profiles were obtained with a precision better than 0.2 micro The results generally agree with predictions of poroelastic models, demonstrating lateral expansion with an effective Poisson's ratio just after completion of the compression phase of the mechanical tests reaching values from 0.18 to 0.4 (depending on compression speed), followed by contraction to lower values. A more restricted movement was observed at both the articular surface and near to the subchondral bone than at regions midway between these two locations.

Animals↗

Lateral resolution in elastography.

The factors that control the lateral resolution in elastography were investigated using a simulation study. The lateral resolution was estimated from the simulated axial strain elastograms as the smallest measurable distance between two equally stiff lesions embedded in a homogeneously softer background. The lesions were symmetrically positioned lateral to the center of the target, at the focus of the transducer. Ultrasound (US) systems with different transducer frequencies, bandwidths and f-numbers were simulated. The effects of the ultrasonic parameters, the lateral spacing between adjacent echo signals, the cross-correlation window length, the lesion/background elastic contrast and the lateral motion of scatterers on the estimated lateral resolution were investigated. The results show that the lateral resolution in elastography is proportional to the beam width of the US system used to acquire the data, and is on the same order as the sonographic lateral resolution.

Computer Simulation↗

Visualisation of HIFU lesions using elastography of the human prostate in vivo: preliminary results.

An imaging system was developed for prostate elastography in vivo using a transrectal ultrasound (US) probe to guide high-intensity focused US (HIFU) therapy of prostate cancer. Uniform compression was applied using a balloon, while a sector image was acquired. Strain was calculated from the gradient of the displacements obtained from the ultrasonic signal using the cross-correlation technique. Elastograms were acquired on a total of 31 patients undergoing HIFU therapy for localised prostate cancer. For two patients, only part of the prostate was treated and posttherapy magnetic resonance imaging (MRI) confirmed the size and position of the HIFU lesions seen in the elastograms as low strain areas, with a strain contrast ratio between 1.6 and 3.2. The whole prostate was treated for the next 29 patients. After treatment, the whole prostate appeared to be stiff in the elastograms and a 40% to 60% (mean 50%) decrease in average strain was observed when compared to strains measured before HIFU application. Tumours identified by biopsies and sonograms could occasionally be seen in the preoperative elastograms. Decorrelation effects occurred mainly because of low sonographic signal-to-noise ratio (SNR) and of out-of-plane motion induced by respiration.

Adenocarcinoma↗

Elastography imaging of small animal oncology models: a feasibility study.

To test the feasibility of applying ultrasonic elastography on small animal oncology models, experiments were performed in vitro and in situ on murine mammary lesions induced exogenously by tumor cell line 66.3. In vitro studies involved three 1-week-old excised tumors embedded in a phantom block with ultrasonic properties similar to those of soft biologic tissues. In situ studies involved five mice whose bodies were embedded in pure gelatin blocks. The data were acquired from the blocks with a clinical scanner modified to have an automated compressor assembly and processed to construct the elastograms at various imaging planes within each block. The results were analyzed both qualitatively and quantitatively to assess the merits of the elastographic imaging and its limitations for in vivo serial studies of tumors in small animal oncology models.

Animals↗