Improved scatterer property estimates from ultrasound backscatter for small gate lengths using a gate-edge correction factor

Backscattered rf signals used to construct conventional ultrasound B-mode images contain frequency-dependent information that can be examined through the backscattered power spectrum. The backscattered power spectrum is found by taking the magnitude squared of the Fourier transform of a gated time segment corresponding to a region in the scattering volume. When a time segment is gated, the edges of the gated regions change the frequency content of the backscattered power spectrum due to truncating of the waveform. Tapered windows, like the Hanning window, and longer gate lengths reduce the relative contribution of the gate-edge effects. A new gate-edge correction factor was developed that partially accounted for the edge effects. The gate-edge correction factor gave more accurate estimates of scatterer properties at small gate lengths compared to conventional windowing functions. The gate-edge correction factor gave estimates of scatterer properties within 5% of actual values at very small gate lengths (less than 5 spatial pulse lengths) in both simulations and from measurements on glass-bead phantoms. While the gate-edge correction factor gave higher accuracy of estimates at smaller gate lengths, the precision of estimates was not improved at small gate lengths over conventional windowing functions.     

A theoretical comparison of attenuation measurement techniques from backscattered ultrasound echoes

Accurate characterization of tissue pathologies using ultrasonic attenuation is strongly dependent on the accuracy of the algorithm that is used to obtain the attenuation coefficient estimates. In this paper, computer simulations were used to compare the accuracy and the precision of the three methods that are commonly used to estimate the local ultrasonic attenuation within a region of interest (ROI) in tissue; namely, the spectral log difference method, the spectral difference method, and the hybrid method. The effects of the inhomgeneities within the ROI on the accuracy of the three algorithms were studied, and the optimal ROI size (the number of independent echoes laterally and the number of pulse lengths axially) was quantified for each method. The three algorithms were tested for when the ROI was homogeneous, the ROI had variations in scatterer number density, and the ROI had variations in effective scatterer size. The results showed that when the ROI was homogeneous, the spectral difference method had the highest accuracy and precision followed by the spectral log difference method and the hybrid method, respectively. Also, when the scatterer number density varied, the spectral difference method completely failed, while the log difference method and hybrid method still gave good results. Lastly, when the scatterer size varied, all of the methods failed.     

Time domain attenuation estimation method from ultrasonic backscattered signals

Ultrasonic attenuation is important not only as a parameter for characterizing tissue but also for compensating other parameters that are used to classify tissues. Several techniques have been explored for estimating ultrasonic attenuation from backscattered signals. In the present study, a technique is developed to estimate the local ultrasonic attenuation coefficient by analyzing the time domain backscattered signal. The proposed method incorporates an objective function that combines the diffraction pattern of the source/receiver with the attenuation slope in an integral equation. The technique was assessed through simulations and validated through experiments with a tissue mimicking phantom and fresh rabbit liver samples. The attenuation values estimated using the proposed technique were compared with the attenuation estimated using insertion loss measurements. For a data block size of 15 pulse lengths axially and 15 beamwidths laterally, the mean attenuation estimates from the tissue mimicking phantoms were within 10% of the estimates using insertion loss measurements. With a data block size of 20 pulse lengths axially and 20 beamwidths laterally, the error in the attenuation values estimated from the liver samples were within 10% of the attenuation values estimated from the insertion loss measurements.     

Method of improved scatterer size estimation and application to parametric imaging using ultrasound

The frequency dependence of RF signals backscattered from random media (tissues) has been used to describe the microstructure of the media. The frequency dependence of the backscattered RF signal is seen in the power spectrum. Estimates of scatterer properties (average scatterer size) from an interrogated medium are made by minimizing the average squared deviation (MASD) between the measured power spectrum and a theoretical power spectrum over an analysis bandwidth. Estimates of the scatterer properties become increasingly inaccurate as the average signal to noise ratio (SNR) over the analysis bandwidth becomes smaller. Some frequency components in the analysis bandwidth of the measured power spectrum will have smaller SNR than other frequency components. The accuracy of estimates can be improved by weighting the frequency components that have the smallest SNR less than the frequencies with the largest SNR in the MASD. A weighting function is devised that minimizes the noise effects on the estimates of the average scatterer sizes. Simulations and phantom experiments are conducted that show the weighting function gives improved estimates in an attenuating medium. The weighting function is applied to parametric images using scatterer size estimates of a rat that had developed a spontaneous mammary tumor.     

Axial coding in full-field microscopy using three-dimensional structured illumination implemented with no moving parts

We report a simple optical setup to produce both axial and lateral structured illumination through a single objective lens. With a minimum of six full-field images obtained without moving either the sample or the microscope objective, 100 nm diameter fluorescent beads can be localized axially with an accuracy of 50 nm in a 1.76-microm-thick layer. We show that this axial localization improvement can easily be combined with classical lateral structured illumination, so that lateral resolution enhancement by a factor of 2 is maintained.     

Time and frequency parameters of bottlenose dolphin whistles as predictors of surface behavior in the Mississippi Sound

In this study, an algorithm previously developed for estimating the total ultrasonic attenuation along the propagation path from the surface of the transducer to a region of interest (ROI) in tissue, was modified to make it more practical for use in clinical settings. Specifically, the algorithm was re-derived for when a tissue mimicking phantom rather than a planar reflector is used to obtain the reference power spectrum. The reference power spectrum is needed to compensate for the transfer function of the transmitted pulse, the transfer function of transducer, and the diffraction effects that result from focusing/beam forming. The modified algorithm was tested on simulated radio frequency (RF) echo lines obtained from two samples that have different scatterer sizes and different attenuation coefficient slopes, one of which was used as a reference. The mean and standard deviation of the percent errors in the attenuation coefficient estimates (ACEs) were less than 5% and 10%, respectively, for ROIs that contain more than 10 pulse lengths and more than 25 independent echo lines. The proposed algorithm was also tested on two tissue mimicking phantoms that have attenuation coefficient slopes of 0.7 dB/cm-MHz and 0.5 dB/cm-MHz respectively, the latter being the reference phantom. When a single element spherically focused source was used, the mean and standard deviation of the percent errors in the ACEs were less than 5% and 10% respectively for windows that contain more than 10 pulse lengths and more than 17 independent echo lines. When a clinical array transducer was used, the mean and standard deviation of the percent errors in the ACEs were less than 5% and 25%, respectively, for windows that contain more than 12 pulse lengths and more than 45 independent echo lines.     

Impact of local attenuation approximations when estimating correlation length from backscattered ultrasound echoes

Estimating the characteristic correlation length of tissue microstructure from the backscattered power spectrum could improve the diagnostic capability of medical ultrasound. Previously, size estimates were obtained after compensating for source focusing, the frequency-dependent attenuation along the propagation path (total attenuation), and the frequency-dependent attenuation in the scattering region (local attenuation). In this study, the impact of approximations of the local attenuation on the scatterer size estimate was determined using computer simulations and theoretical analysis. The simulations used Gaussian impedance distributions with an effective radius of 25 microm randomly positioned in a homogeneous half-space sonified by a spherically focused source (f/1 to f/4). The approximations of the local attenuation that were assessed neglected local attenuation (i.e., assume 0 dB/cm-MHz) neglected frequency dependence of the local attenuation, and assumed a finite frequency dependence (i.e., 0.5 dB/cm-MHz) independent of the true attenuation of the medium. Errors in the scatterer size estimate due to the local attenuation approximations increased with increasing window length, increasing true local attenuation and increasing f number. The most robust estimates were obtained when the local attenuation was approximated by a tissue-independent attenuation value that was greater than 70% of the largest attenuation expected in the tissue region of interest.     

Fourier-domain angle-resolved low coherence interferometry through an endoscopic fiber bundle for light-scattering spectroscopy

We present a novel endoscopic fiber bundle probe incorporated in a Fourier-domain angle-resolved low coherence interferometry system for the measurement of depth-resolved angular scattering distributions to permit the determination of scatterer size via elastic scattering properties. Depth resolution is achieved with a superluminescent diode via a Mach-Zehnder interferometer. The sample is illuminated with a collimated beam, and a Fourier plane image of the backscattered light is collected by a coherent fiber bundle. The angular scattering distribution relayed by the fiber bundle is mixed with the reference field and made to coincide with the input slit of an imaging spectrograph. The data collected are processed in real time, producing a depth-resolved angular scattering distribution in 0.37 s. The data are used to determine the sizes of polystyrene microspheres with subwavelength precision and accuracy.     

Identifying ultrasonic scattering sites from three-dimensional impedance maps

Ultrasonic backscattered signals contain frequency-dependent information that is usually discarded to produce conventional B-mode images. It is hypothesized that parametrization of the quantitative ultrasound frequency-dependent information (i.e., estimating scatterer size and acoustic concentration) may be related to discrete scattering anatomic structures in tissues. Thus, an estimation technique is proposed to extract scatterer size and acoustic concentration from the power spectrum derived from a three-dimensional impedance map (3DZM) of a tissue volume. The 3DZM can be viewed as a computational phantom and is produced from a 3D histologic data set. The 3D histologic data set is constructed from tissue sections that have been appropriately stained to highlight specific tissue features. These tissue features are assigned acoustic impedance values to yield a 3DZM. From the power spectrum, scatterer size and acoustic concentration estimates were obtained by optimization. The 3DZM technique was validated by simulations that showed relative errors of less than 3% for all estimated parameters. Estimates using the 3DZM technique were obtained and compared against published ultrasonically derived estimates for two mammary tumors, a rat fibroadenoma and a 4T1 mouse mammary carcinoma. For both tumors, the relative difference between ultrasonic and 3DZM estimates was less than 10% for the average scatterer size.     

Binary-phase spatial filter for real-time swept-source optical coherence microscopy

We report a novel scheme to optimize the focusing condition for real-time, swept-source optical coherence microscopy. The axial and lateral behaviors of four-zone binary-phase spatial filters are presented numerically. A nearly constant axial intensity distribution along an extended depth of focus of 1.5 mm and a lateral resolution of 5 microm are experimentally verified. The A-line scan rate is up to 16 kHz, yielding a frame rate of 25 Hz and 640 lines per image.     

A new method of ultrasound color flow mapping

Conventional ultrasound color flow mapping systems estimate and visualize only the axial velocity component. To obtain the transverse velocity component a modification of a multiple-beam method is proposed. The new two-dimensional color flow mapping system has a small size and consists of three transducers. The central transducer is an appodized and focused phased array. The other transducers are unfocused probes. Three transducers act as receivers and the central transducer operates as a transmitter. All receivers acquire rf scan lines that are then processed to estimate three axial velocity components using an autocorrelation method. These estimates are then combined to estimate the transverse velocity component, taking into account the geometric relationships among three transducers. Two algorithms for transverse velocity estimation are proposed. The first uses the Doppler angle estimate for calculation of the transverse velocity component. The other algorithm calculates the transverse velocity component directly from the axial components. The accuracy of the flow velocity estimators is estimated by simulations. Analysis of accuracy allows choosing the more effective algorithm for two-dimensional velocity estimation, which is insensitive to variations of the Doppler angle.     

Parabola detection using matched filtering for ultrasound B-scans

Time of flight diffraction (ToFD) outputs B-scans using an ultrasound emitter and receiver at a constant separation, scanned over a sample surface parallel to the line between the transducers. The B-scan, with time and scan position axes, contains parabolic features characteristic of a point-like scatterer. Human vision effectively detects these shapes, but this is time consuming and requires training. A parabola matched filter has been developed that is simple to implement and transforms parabolic shapes to peaks whilst reducing noise in the scan. The scan can then be displayed as depth versus lateral position.     

Axial super-resolution by mirror-reflected stimulated emission depletion microscopy

In stimulated emission depletion (STED) microscopy, the lateral resolution is in the range of tens of nanometers depending on the sample and the instrument. The axial resolution, however, is in standard systems limited by diffraction to about 500 nm. We present an approach to three-dimensional diffraction-unlimited resolution by observing the sample at two optical angles. The system is realized by using an atomic force microscope (AFM) chip as a microreflector to deflect the STED beams near the region-of-interest (ROI), thus allowing observations at an angle ∠. Consequently, the superior lateral resolution can be utilized to resolve details in the axial direction of the main optical axis of the microscope. Here, fluorescent nanoparticles 90 nm apart and biological structures 80 nm apart along axial direction were distinguished by utilizing an off-the-shelf, commercial STED microscope, coupled with an AFM and an AFM chip micro-reflector.     

远场超分辨随机光重建显微镜(STORM)研究进展

了解细胞内分子尺度的动态和结构的特征是生命科学迫切需要解决的问题,要求远场光学成像要求纳米或亚纳米量级的空间分辨率.介绍了一种实现打破衍射极限的远场荧光显微成像技术--随机光重建显微术(STORM),其分辨率可以达到横向分辨率20 nm,轴向分辨率50 nm,理论上这种方法的空间分辨率可以达到单分子定位的精度.具体介绍...     

Accelerated 3D Coronary Vessel Wall MR Imaging Based on Compressed Sensing with a Block-Weighted Total Variation Regularization

Coronary vessel wall magnetic resonance (MR) imaging is important for heart disease diagnosis but often suffers long scan time. Compressed sensing (CS) has been previously used to accelerate MR imaging by reconstructing an MR image from undersampled k-space data using a regularization framework. However, the widely used regularizations in the current CS methods often lead to smoothing effects and thus are unable to reconstruct the coronary vessel walls with sufficient resolution. To address this issue, a novel block-weighted total variation regularization is presented to accelerate the coronary vessel wall MR imaging. The proposed regularization divides the image into two parts: a region-of-interest (ROI) which contains the coronary vessel wall, and the other region with less concerned features. Different penalty weights are given to the two regions. As a result, the small details within ROI do not suffer from over-smoothing while the noise outside the ROI can be significantly suppressed. Results with both numerical simulations and in vivo experiments demonstrated that the proposed method can reconstruct the coronary vessel wall from undersampled k-space data with higher qualities than the conventional CS with the total variation or the edge-preserved total variation.     

Correlation of ultrasonic scatterer size estimates for the statistical analysis and optimization of angular compounding

Ultrasonic scatterer size estimates generally have large variances due to the inherent noise of spectral estimates used to calculate size. Compounding partially correlated size estimates associated with the same tissue, but produced with data acquired from different angles of incidence, is an effective way to reduce the variance without making dramatic sacrifices in spatial resolution. This work derives theoretical approximations for the correlation between these size estimates, and the coherence between their associated spectral estimates, as functions of ultrasonic system parameters. A Gaussian spatial autocorrelation function is assumed to adequately model scatterer shape. Both approximations compare favorably with simulation results, which consider validation near the focus. Utilization of the correlation/coherence expressions for statistical analysis and optimization is discussed. Approximations, such as the invariance of phase and amplitude terms with angle, are made to obtain closed-form solutions to the derived spectral coherence near the focus and permit analytical optimization analysis. Results indicate that recommended parameter adjustments for performance improvement generally depend upon whether, for the system under consideration, the primary source of change in total coherence with rotation is phase term variation due to the change in the relative position of scattering sites, or field amplitude term variation due to beam movement.     

Estimate of the attenuation coefficient using a clinical array transducer for the detection of cervical ripening in human pregnancy

Premature delivery is the leading cause of infant mortality in the United States. Currently, premature delivery cannot be prevented and new treatments are difficult to develop due to the inability to diagnose symptoms prior to uterine contractions. Cervical ripening is a long period that precedes the active phase of uterine contractions and cervical dilation. The changes in the microstructure of the cervix during cervical ripening suggest that the ultrasonic attenuation should decrease. The objective of this study is to use the reference phantom algorithm to estimate the ultrasonic attenuation in the cervix of pregnant human patients. Prior to applying the algorithm to in vivo human data, two homogeneous phantoms with known attenuation coefficients were used to validate the algorithm and to find the length and the width of the region of interest (ROI) that achieves the smallest error in the attenuation coefficient estimates. In the phantom data, we found that the errors in the attenuation coefficients estimates are less than 12% for ROIs that contain 40 wavelengths or more axially and 30 echo lines or more laterally. The reference phantom algorithm was then used to obtain attenuation maps of the echoes from two human pregnant cervices at different gestational ages. It was observed that the mean of the attenuation coefficient estimates in the cervix of the patient at a more advanced gestational age is smaller than the mean of the attenuation coefficient estimates in the cervix of the patient at an earlier gestational age which suggests that ultrasonic attenuation decreases with increasing gestational age. We also observed a large variance between the attenuation coefficient estimates in the different regions of the cervix due to the natural variation in tissue micro-structures across the cervix. The preliminary results indicate that the algorithm could potentially provide an important diagnostic tool for diagnosing the risk of premature delivery.     

Estimation and application of spatially variable noise fields in diffusion tensor imaging

Optimal interpretation of magnetic resonance image content often requires an estimate of the underlying image noise, which is typically realized as a spatially invariant estimate of the noise distribution. This is not an ideal practice in diffusion tensor imaging because the noise distribution is usually spatially varying due to the use of fast imaging and noise suppression techniques. A new estimation approach for spatially varying noise fields (NFs) is proposed in this article. The approach is based on a noise invariance property in scenarios in which more than one image, each with potentially different signal levels, is acquired on each slice, as in diffusion-weighted MRI. This technique leads to improved NF estimates in simulations, phantom experiments and in vivo studies when compared to traditional NF estimators that use regional variability or background intensity histograms. The proposed method reduces the NF estimation error by a factor of 100 in simulations, shows a strong linear correlation (R²=0.99) between theoretical and estimated noise changes in phantoms and demonstrates consistent (<5% variability) NF estimates in vivo. The advantages of spatially varying NF estimation are demonstrated for power analysis, outlier detection and tensor estimation.     

Errors in ultrasonic scatterer size estimates due to phase and amplitude aberration

Current ultrasonic scatterer size estimation methods assume that acoustic propagation is free of distortion due to large-scale variations in medium attenuation and sound speed. However, it has been demonstrated that under certain conditions in medical applications, medium inhomogeneities can cause significant field aberrations that lead to B-mode image artifacts. These same aberrations may be responsible for errors in size estimates and parametric images of scatterer size. This work derives theoretical expressions for the error in backscatter coefficient and size estimates as a function of statistical parameters that quantify phase and amplitude aberration, assuming a Gaussian spatial autocorrelation function. Results exhibit agreement with simulations for the limited region of parameter space considered. For large values of aberration decorrelation lengths relative to aberration standard deviations, phase aberration errors appear to be minimal, while amplitude aberration errors remain significant. Implications of the results for accurate backscatter and size estimation are discussed. In particular, backscatter filters are suggested as a method for error correction. Limitations of the theory are also addressed. The approach, approximations, and assumptions used in the derivation are most appropriate when the aberrating structures are relatively large, and the region containing the inhomogeneities is offset from the insonifying transducer.     

Ultrasound line-by-line scanning method of spatial–temporal active cavitation mapping for high-intensity focused ultrasound

This paper presented an ultrasound line-by-line scanning method of spatial–temporal active cavitation mapping applicable in a liquid or liquid filled tissue cavities exposed by high-intensity focused ultrasound (HIFU). Scattered signals from cavitation bubbles were obtained in a scan line immediately after one HIFU exposure, and then there was a waiting time of 2 s long enough to make the liquid back to the original state. As this pattern extended, an image was built up by sequentially measuring a series of such lines. The acquisition of the beamformed radiofrequency (RF) signals for a scan line was synchronized with HIFU exposure. The duration of HIFU exposure, as well as the delay of the interrogating pulse relative to the moment while HIFU was turned off, could vary from microseconds to seconds. The feasibility of this method was demonstrated in tap-water and a tap-water filled cavity in the tissue-mimicking gelatin–agar phantom as capable of observing temporal evolutions of cavitation bubble cloud with temporal resolution of several microseconds, lateral and axial resolution of 0.50 mm and 0.29 mm respectively. The dissolution process of cavitation bubble cloud and spatial distribution affected by cavitation previously generated were also investigated. Although the application is limited by the requirement for a gassy fluid (e.g. tap water, etc.) that allows replenishment of nuclei between HIFU exposures, the technique may be a useful tool in spatial–temporal cavitation mapping for HIFU with high precision and resolution, providing a reference for clinical therapy.