Speckle reduction for medical ultrasound images with an expectation-maximization framework

In view of inherent speckle noise in medical images,a speckle reduction method was proposed based on an expectation-maximization(EM) framework.First,the real component of the in-phase/quadrature(I/Q) ultrasound image is extracted.Then,it is used to blindly estimate the point spread function(PSF) of the imaging system.Finally,based on the EM framework,an iterative algorithm alternating between the Wiener Filter and the anisotropic diffusion(AD) is exploited to produce despeckled images.The comparison experiment is carried out on both simulated and in vivo ultrasound images.It is shown that,with respect to the I/Q image,the proposed method averagely improves the speckle-signal-to-noise ratio(S-SNR) and the edge preservation index(β) of images by the factor of 1.94 and 7.52.Meanwhile,it averagely reduces the normalized mean-squared error(NMSE) by the factor of 3.95.The simulation and in vivo results indicates that the proposed method has a better overall performance than exited ones.     

Analysis of speckle in ultrasound images using fractional order statistics and the homodyned k-distribution

It is necessary to identify speckled regions in ultrasound images to control adaptive speckle suppression algorithms, for tissue characterisation, and to estimate the elevational separation of B-scans by speckle decorrelation. Previous authors have proposed classification techniques based on second order powers of the homodyned k-distribution, or lower order powers of the more limited k-distribution. In this paper we explore the speckle discrimination properties of statistics based on arbitrary powers of the ultrasound echo envelope signal using a combination of simulations and theoretical results from the homodyned k-distribution. We conclude that statistics based on powers less than one are surprisingly less effective than some higher powers. A simple discriminant function for speckle is evaluated quantitatively in simulation and qualitatively on sample B-scan images.     

Effects of red blood cell aggregates dissociation on the estimation of ultrasound speckle image velocimetry

Ultrasound speckle image of blood is mainly attributed by red blood cells (RBCs) which tend to form RBC aggregates. RBC aggregates are separated into individual cells when the shear force is over a certain value. The dissociation of RBC aggregates has an influence on the performance of ultrasound speckle image velocimetry (SIV) technique in which a cross-correlation algorithm is applied to the speckle images to get the velocity field information. The present study aims to investigate the effect of the dissociation of RBC aggregates on the estimation quality of SIV technique. Ultrasound B-mode images were captured from the porcine blood circulating in a mock-up flow loop with varying flow rate. To verify the measurement performance of SIV technique, the centerline velocity measured by the SIV technique was compared with that measured by Doppler spectrograms. The dissociation of RBC aggregates was estimated by using decorrelation of speckle patterns in which the subsequent window was shifted as much as the speckle displacement to compensate decorrelation caused by in-plane loss of speckle patterns. The decorrelation of speckles is considerably increased according to shear rate. Its variations are different along the radial direction. Because the dissociation of RBC aggregates changes ultrasound speckles, the estimation quality of SIV technique is significantly correlated with the decorrelation of speckles. This degradation of measurement quality may be improved by increasing the data acquisition rate. This study would be useful for simultaneous measurement of hemodynamic and hemorheological information of blood flows using only speckle images.     

An adaptive squared-distance-weighted interpolation for volume reconstruction in 3D freehand ultrasound

Volume reconstruction is a key procedure in 3D ultrasound imaging. An algorithm named as squared-distance-weighted (SDW) interpolation has been earlier proposed to reduce the blurring effect in the 3D ultrasonic images caused by the conventional distance weighted (DW) interpolation. However, the SDW parameter alpha, which controls the weight distribution, is a constant assigned by operators so that the interpolation effect is invariant for both sharp edges and speckle noises. In this paper, we introduced a new adaptive algorithm based on SDW interpolation for volume reconstruction of 3D freehand ultrasound. In the algorithm, the local statistics of pixels surrounding each voxel grid were used to adaptively adjust the parameter alpha in SDW. The voxel grids with a higher ratio of local variance and mean in their neighbourhoods would have a smaller alpha to make the image details sharper, while the voxel grids locating in regions with a lower ratio of local variance and mean would have a larger alpha to smooth image content in homogeneous regions, where speckle noise is usually observed and damages the image quality. By comparing the simulation results using the SDW and new adaptive algorithm, it was demonstrated that this new algorithm worked well in both edge preservation and speckle reduction.     

Speckle noise reduction by division and digital processing of a hologram

One of the common problems in classical and digital holography is presence of speckle in reconstructed images. Some methods have been predicted for speckle reduction. Speckle has a statistical distribution in size and intensity. Thus complete elimination of them is impossible. But, one can use this statistical characters to reduce its effects.In this paper a method is proposed for reduction of speckle effect by recording a hologram and then dividing it to several sub-holograms. Each sub-hologram is used to produce an image. Final image is constructed by superposition of these images. The final image has a smaller amount of speckle noise. Experimental results verify effectiveness of this technique.     

Speckle detection for 3D ultrasound

Positioning accuracy of 3D ultrasound slice planes can be enhanced by speckle decorrelation of the B-scans. A suitable method requires accurate identification of zones of non-specular reflection in the image, a feature of parenchymal tissue. This has traditionally been achieved with statistics such as signal-to-noise ratios and measures of autocorrelation curves. Co-occurrence matrices and structural approaches have also been applied to ultrasound images. In this paper we assess the performance of these methods when only B-scan display data is available. We also report improved results with in vivo images using a proposed structural algorithm that detects speckle kernels. All detection algorithms yielded detection error rates between 28% and 38%, suggesting that user intervention may be required for clinically useful performance.     

Reconstruction of 3D ultrasound images based on Cyclic Regularized Savitzky-Golay filters

This paper presents a new three-dimensional (3D) ultrasound reconstruction algorithm for generation of 3D images from a series of two-dimensional (2D) B-scans acquired in the mechanical linear scanning framework. Unlike most existing 3D ultrasound reconstruction algorithms, which have been developed and evaluated in the freehand scanning framework, the new algorithm has been designed to capitalize the regularity pattern of the mechanical linear scanning, where all the B-scan slices are precisely parallel and evenly spaced. The new reconstruction algorithm, referred to as the Cyclic Regularized Savitzky-Golay (CRSG) filter, is a new variant of the Savitzky-Golay (SG) smoothing filter. The CRSG filter has been improved upon the original SG filter in two respects: First, the cyclic indicator function has been incorporated into the least square cost function to enable the CRSG filter to approximate nonuniformly spaced data of the unobserved image intensities contained in unfilled voxels and reduce speckle noise of the observed image intensities contained in filled voxels. Second, the regularization function has been augmented to the least squares cost function as a mechanism to balance between the degree of speckle reduction and the degree of detail preservation. The CRSG filter has been evaluated and compared with the Voxel Nearest-Neighbor (VNN) interpolation post-processed by the Adaptive Speckle Reduction (ASR) filter, the VNN interpolation post-processed by the Adaptive Weighted Median (AWM) filter, the Distance-Weighted (DW) interpolation, and the Adaptive Distance-Weighted (ADW) interpolation, on reconstructing a synthetic 3D spherical image and a clinical 3D carotid artery bifurcation in the mechanical linear scanning framework. This preliminary evaluation indicates that the CRSG filter is more effective in both speckle reduction and geometric reconstruction of 3D ultrasound images than the other methods.     

基于Rician分布的散斑噪声的超声图像模拟算法*

为了衡量各种去噪算法的性能,在干净无噪声的图像上添加接近真实且可控的散斑噪声是非常重要的,提出一种基于Rician分布的不完全发育的斑点噪声的超声图像模拟算法。该算法考虑到了声波区域中包含孤立的强散射体的情况,同时结合了超声成像的扫描过程。以合成图像和肾脏图像为体模进行了模拟实验,并对最终生成的伪超声图像进行了噪声分布统计及拟合检验。实验结果表明该算法生成的伪超声图像在视觉上和理论上都接近真实的超声图像。     

Speckle noise attenuation in optical coherence tomography by compounding images acquired at different positions of the sample

This study demonstrates a simple method for attenuating the speckle noise generated by coherent multiple-scattered photons in optical-coherence tomography images. The method could be included among the space-diversity techniques used for speckle reduction. It relies on displacing the sample along a weakly focused beam in the sample arm of the interferometer, acquiring a coherent image for each sample position and adding the individual images to form a compounded image. It is proven that the compounded image displays a reduction in the speckle noise generated by multiple scattered photons and an enhancement in the intensity signal caused by single-backscattered photons. To evaluate its potential biomedical applications, the method is used to investigate in vitro a caries lesion affecting the enamel layer of a wisdom tooth. Because of the uncorrelated nature of the speckle noise the compounded image provides a better mapping of the lesion compared to a single (coherent) image.     

Mechanism of speckle reduction in laser-microscope images using a rotating optical fiber

The reduction mechanism of speckles appearing in laser-microscope images is studied theoretically and experimentally when an object is illuminated by laser light through a rotating multimode optical fiber. The principle of the speckle reduction is based on independent addition of microscope images with boiling-like speckles as a result of the rotational motion of the optical fiber used for illumination. Especially, the speckle reduction is evaluated from first-order statistics of the speckle intensity in which its saturation effect is taken into consideration.     

Superresolution imaging of scatterers in ultrasound B-scan imaging

A number of imaging systems exhibit speckle, which is caused by the interaction of a coherent pulse reflecting off of random reflectors. The limitations of these systems are quite serious because the speckle phenomenon creates a pattern of nulls and peaks from subresolvable scatterers or targets that are difficult to interpret. Another limitation of these pulse-echo imaging systems is that their resolution is dependent on the full spatial extent of the propagating pulse, usually several wavelengths in the axial or propagating dimension and typically longer in the transverse direction. This limits the spatial resolution to many multiples of the wavelength. This paper focuses on the particular case of ultrasound B-scan imaging and develops an inverse filter solution that eliminates both the speckle phenomenon and the poor resolution dependency on the pulse length and width to produce super-resolution ultrasound (SURUS) images. The key to the inverse filter is the creation of pulse shapes that have stable inverses. This is derived by use of the standard Z-transform and related properties. Although the focus of this paper is on examples from ultrasound imaging systems, the results are applicable to other pulse-echo imaging systems that also can exhibit speckle statistics.     

Real-time implementation of coherent speckle suppression in B-scan images

Early experiences of new forms of adaptive filtering for ultrasound speckle reduction and parametric imaging, using off-line conventional digital processing, have been sufficiently encouraging to warrant examining the feasibility of implementing specific algorithms in real-time. A hardware two-dimensional real-time filter is described which consists of a hybrid digital/analogues system in which the video signal from any scanner is sampled to 256 points per line and passed sequentially through a series of shift registers, in order to derive a 5 x 5 window of values which surrounds the image point currently being processed. These 25 video signals are then used as inputs to an analogue processor, which provides the filtered output. The real-time processed images show clear evidence of speckle smoothing without blurring of tissue structural information but possess limited pixel resolution.     

Adaptive filtering for reduction of speckle in ultrasonic pulse-echo images

Current medical ultrasonic scanning instrumentation permits the display of fine image detail (speckle) which does not transfer useful information but degrades the apparent low contrast resolution in the image. An adaptive two-dimensional filter has been developed which uses local features of image texture to recognize and maximally low-pass filter those parts of the image which correspond to fully developed speckle, while substantially preserving information associated with resolved-object structure. A first implementation of the filter is described which uses the ratio of the local variance and the local mean as the speckle recognition feature. Preliminary results of applying this form of display processing to medical ultrasound images are very encouraging; it appears that the visual perception of features such as small discrete structures, subtle fluctuations in mean echo level and changes in image texture may be enhanced relative to that for unprocessed images.     

Improvement of ultrasound speckle image velocimetry using image enhancement techniques

Ultrasound-based techniques have been developed and widely used in noninvasive measurement of blood velocity. Speckle image velocimetry (SIV), which applies a cross-correlation algorithm to consecutive B-mode images of blood flow has often been employed owing to its better spatial resolution compared with conventional Doppler-based measurement techniques. The SIV technique utilizes speckles backscattered from red blood cell (RBC) aggregates as flow tracers. Hence, the intensity and size of such speckles are highly dependent on hemodynamic conditions. The grayscale intensity of speckle images varies along the radial direction of blood vessels because of the shear rate dependence of RBC aggregation. This inhomogeneous distribution of echo speckles decreases the signal-to-noise ratio (SNR) of a cross-correlation analysis and produces spurious results. In the present study, image-enhancement techniques such as contrast-limited adaptive histogram equalization (CLAHE), min/max technique, and subtraction of background image (SB) method were applied to speckle images to achieve a more accurate SIV measurement. A mechanical sector ultrasound scanner was used to obtain ultrasound speckle images from rat blood under steady and pulsatile flows. The effects of the image-enhancement techniques on SIV analysis were evaluated by comparing image intensities, velocities, and cross-correlation maps. The velocity profiles and wall shear rate (WSR) obtained from RBC suspension images were compared with the analytical solution for validation. In addition, the image-enhancement techniques were applied to in vivo measurement of blood flow in human vein. The experimental results of both in vitro and in vivo SIV measurements show that the intensity gradient in heterogeneous speckles has substantial influence on the cross-correlation analysis. The image-enhancement techniques used in this study can minimize errors encountered in ultrasound SIV measurement in which RBCs are used as flow tracers instead of exogenous contrast agents.     

Laser speckle reduction using motionless image conduits

We have demonstrated a speckle reduction method using motionless image conduits (MICs). Different experimental conditions by introducing the high-coherence HeNe laser and the low-coherence laser diode (LD) as the illumination light sources, by employing the straight MIC and the curved MIC as the speckle reduction components, and by recording speckle images without (objective speckle) and with (subjective speckle) the imaging lens mounted on the CCD camera are conducted, respectively. The most efficient speckle reduction condition is found by the combination of using the LD and the curved MIC, where the objective speckle contrast ratio is reduced from 0.7378 to 0.1725. Experimental results are discussed, and the causes for these speckle reduction efficiency changes are given.     

Image features in medical vibro-acoustography: in vitro and in vivo results

Vibro-acoustography is an imaging method based on audio-frequency harmonic vibrations induced in the object by the radiation force of focused ultrasound. The purpose of this study is to investigate features of vibro-acoustography images and manifestation of various tissue structures and calcifications in such images. Our motivation for this study is to pave the way for further in vitro and in vivo applications of vibro-acoustography. Here, vibro-acoustography images of excised prostate and in vivo breast are presented and compared with images obtained with other modalities. Resulting vibro-acoustography images obtained with a 3 MHz ultrasound transducer and at a vibration frequency of 50-60 kHz show soft tissue structures, tissue borders, and microcalcifications with high contrast, high resolution, and no speckle. It is concluded that vibro-acoustography offers features that may be valuable for diagnostic purposes.     

Ultrasound speckle and equivalent scatterers

B-mode ultrasound images are characterised by the speckle artefact, which introduces fine-false structures whose apparent resolution is beyond the imaging system capabilities. Speckle presence is due to interference effects between overlapping echoes and its occurrence is related to a great number of randomly distributed structure scatterers within a resolution cell. Basing our analysis on linear system theory, we show that a dense random set of scatterers can be substituted by an equivalent one with a much smaller number of periodic scatterers. This new structure with regularly distributed scatterers is able to give rise to the same B-mode image and the same speckle pattern, for a given ultrasound pulse. This new approach helps the understanding of the deterministic nature of speckle and may reduce drastically the computing time in numerical simulations. Additionally, it can contribute to periodicity analysis used in tissue characterisation.     

Optimally discriminant moments for speckle detection in real B-scan images

The paper presents and evaluates a speckle detection method for B-scan images. This is a fully automatic method and does not require information about the sensor parameters, which is often missing in retrospective studies.The characterization and posterior detection of speckle noise in ultrasound (US) has been regarded as an important research topic in US imaging, for improving signal-to-noise ratio by removing speckle noise and for exploiting speckle correlation information. Most of the existing methods require either manual intervention, the need to know sensor parameters or are based on statistical models which often do not generalize well to B-scans of different imaging areas. The proposed method aims to overcome those limitations.The main novelty of this work is to show that speckle detection can be improved based on finding optimally discriminant low order speckle statistics. In addition, and in contrast with other approaches the presented method is fully automatic and can be efficiently implemented to B-scan images.The method detects speckle patches using an ellipsoid discriminant function which classifies patches based on features extracted from optimally discriminant low order moments of the uncompressed intensity B-scan information. In addition, if the uncompressed signal is not available, we propose and evaluate a method for the estimation of this factor.The computation of low order moments using an optimality criteria, the decompression factor estimation and other key aspects of the method are quantitatively evaluated using both simulated and real (phantom and in vivo) data. Speckle detection results are obtained using again phantom and in vivo studies which show the validity of our approach. In addition, speckle probability images (SPI) are presented which provide valuable information about the distribution of speckle and non-speckle areas in an image.The presented evaluation and results show the effectiveness of our approach. In particular, the need for using discriminant analysis to determine the optimal discriminant power of the statistical moments and that this optimal value strongly depends on the characteristics and imaged tissues in the B-scan data.     

Speckle noise reduction in digital holography by use of multiple polarization holograms

A novel speckle reduction technique for digital holography is proposed. Multiple off-axis holograms are recorded using a circularly polarized illumination beam and a rotating linearly polarized reference beam. The speckle noise in the reconstructed images is suppressed by averaging these fields. We demonstrate the effectiveness of this technique experimentally and conduct additional statistical evaluation.     

Understanding the exposure-time effect on speckle contrast measurements for laser displays

To evaluate the influence of exposure time on speckle noise for laser displays, speckle contrast measurement method was developed observable at a human eye response time using a high-sensitivity camera which has a signal multiplying function. The nonlinearity of camera light sensitivity was calibrated to measure accurate speckle contrasts, and the measuring lower limit noise of speckle contrast was improved by applying spatial-frequency low pass filter to the captured images. Three commercially available laser displays were measured over a wide range of exposure times from tens of milliseconds to several seconds without adjusting the brightness of laser displays. The speckle contrast of raster-scanned mobile projector without any speckle-reduction device was nearly constant over various exposure times. On the contrary to this, in full-frame projection type laser displays equipped with a temporally-averaging speckle-reduction device, some of their speckle contrasts close to the lower limits noise were slightly increased at the shorter exposure time due to the noise. As a result, the exposure-time effect of speckle contrast could not be observed in our measurements, although it is more reasonable to think that the speckle contrasts of laser displays, which are equipped with the temporally-averaging speckle-reduction device, are dependent on the exposure time. This discrepancy may be attributed to the underestimation of temporal averaging factor. We expected that this method is useful for evaluating various laser displays and clarify the relationship between the speckle noise and the exposure time for a further verification of speckle reduction.