A noncoherent correlation technique and focused beamforming forultrasonic underwater imaging: a comparative analysis

The purpose of this paper is to describe a noncoherent correlation technique for the formation of 3-D acoustic images, mainly for underwater applications. This technique is applied to the envelope of signals received by a 2-D array. The performances of the technique are compared with those of the coherent technique of focused beamforming in terms of angular resolution, depth of field, beam pattern, and accuracy in range determination. Moreover, special attention is given to the modalities for an efficient implementation on a digital computer. The complete absence of grating lobes, for whatever arrangement of the transducers, and the much higher speed of the imaging process, as compared with that of beamforming in the time domain, are the two major advantages of this technique. A reduction in the speckle effect, the possibility of creating a virtual depth of field, the feasibility of using square-law transducers, and the implementation simplicity are additional interesting characteristics. The results of some simulations are reported and compared with those obtained by the beamforming process     

Imaging of trace species distributions by degenerate four-wave mixing: diffraction effects, spatial resolution, and image referencing

We describe the theory of imaging by degenerate four-wave mixing (DFWM) using a standard diffraction theory of imaging by coherent light. We demonstrate that, even with the phase-conjugating geometry, no aberration correction can be achieved by DFWM imaging. We demonstrate the coherent nature of DFWM image formation using spatially modulated signals generated in flame OH in the phase-conjugating geometry. The intensity distribution in the Fourier plane of a telecentric lens system is shown to be the spatial Fourier transform of the object distribution characteristic of coherent imaging. The brightness of the DFWM signals exceeds that of similar laser-induced fluorescence signals that can be discriminated by restricting the aperture of the imaging system while still allowing a spatial resolution of approximately 70 ?m. DFWM imaging with the forward-folded boxcars geometry is demonstrated and used in a simple referencing scheme to compensate for structure on the images imposed by nonuniformity of the laser beams employed. Images formed in NO are used to illustrate that structure on a scale of less than 100 ?m arising from beam inhomogeneity can be removed by this referencing technique.     

Beamspace adaptive beamforming for ultrasound imaging

Applying the Capon adaptive beamformer in medical ultrasound imaging results in enhanced resolution by improving the interference-suppressing capabilities of the array. This improvement comes at the expense of an increased computational complexity. We have investigated the application of a beamspace adaptive beamformer for medical ultrasound imaging, which can be used to achieve reduced computational complexity with performance comparable to that of the Capon beamformer. The idea behind beamspace beamforming is that, instead of using the spatial statistics of the elements in the array to differentiate between signals and interference, we use the spatial statistics of a set of orthogonal beams, which are formed in different directions. This represents a shift from element space to beamspace. Because the majority of interference in medical ultrasound imaging is constrained to a limited spatial interval due to the focused transmit beam, this latter space can be reduced to a dimension that is lower than that of element space. We show, using simulations and experimental data, that this dimension can be selected as low as 3 while still achieving performance comparable to its element space counterpart.     

Phase-aberration correction using signals from point reflectors and diffuse scatterers: basic principles

Methods for correction of phase aberrations induced by near-field variations in the index of refraction are explored. Using signals obtained from a sampled aperture (i.e. transducer array), phase aberrations can be accurately measured with a correlation approach similar to methods used in adaptive optics and radar. However, the method presented here has no need for a beacon or an ideal point reflector to act as a source for estimating phase errors. It uses signals from random collections of scatterers to determine phase aberrations accurately. Because there is no longer a need for a beacon signal, the method is directly applicable not only to medical ultrasound imaging but also to any coherent imaging system with a sampled aperture, such as radar and sonar.     

Ultrasound speckle reduction using harmonic oscillator models

A speckle reduction algorithm called the harmonic imaging (HI) algorithm is presented. It is based on a multicomponent scattering model for medical ultrasonics. The backscattered ultrasound quadrature signal is modeled as the sum of three components after demodulation. The first component represents nonresolvable diffuse scatterers, while the second component represents subresolvable quasi-periodic scatterers. The third component represents resolvable quasi-periodic scatterers and mirroring surfaces. Since the second component gives rise to the most long range destructive interference effects it is eliminated in the HI algorithm to reduce speckle. Due to its slow spatial variation, it can be almost completely eliminated simply by differentiating the backscattered demodulated quadrature signal. Lissajous-like figures are observed in complex plots of the signals from ultrasound beams going through tissues with quasi-periodic components and sometimes in areas with only diffuse scatterers. Therefore the sum of the complex signals from the resolvable and nonresolvable scatterers within a resolution cell is modeled by two orthogonal and independent harmonic oscillators. The estimated, total energy of these two oscillators determines the gray level value of the HI image within the resolution cell. The HI images produced using radio frequency data from a phantom and from tissues in vivo are more blurred than ordinary envelope images, but the signal to noise ratio and tissue contrast were higher for the HI images     

Spaser Nanoparticles for Ultranarrow Bandwidth STED Super-Resolution Imaging

Super-resolution microscopy, as a powerful tool of seeing abundant spatial details, typically can only distinguish a few distinct targets at a time due to the spectral crosstalk between fluorophores. Spaser (i.e., surface plasmon laser) nanoprobes, which confine lasing emission into nanoscale, offer an opportunity to eliminate such obstacle. Here, realized is narrow band stimulated emission depletion (STED) nanoscopy on spaser nanoparticles by collecting the coherent spasing signals. Demonstrated are the physics concept and feasibility of erasing spaser emission by using a depletion beam to suppress the population inversion, which lays the foundation of spaser-based STED super-resolution. Thanks to the small size (47 nm) and narrow spectral linewidth (3.8 nm) of the spaser nanoparticles, a 74 nm spatial resolution in STED imaging within an acquisition bandwidth of 10 nm is finally obtained. These spaser nanoparticles, if multiplexing with different wavelengths, in principle, allow for spectral-multiplexed imaging, sensing, cytometry, and light operation of a large number of targets all at once.     

Interpolation methods for time-delay estimation using cross-correlation method for blood velocity measurement

The cross-correlation method (CCM) for blood flow velocity measurement using Doppler ultrasound is based on time delay estimation of echoes from pulse-to-pulse. The sampling frequency of the received signal is usually kept as low as possible in order to reduce computational complexity, and the peak in the correlation function is found by interpolating the correlation function. The parabolic-fit interpolation method introduces a bias at low sampling rate to the ultrasound center frequency ratio. In this study, four different methods are suggested to improve the estimation accuracy: (1) Parabolic interpolation with bias-compensation, derived from a theoretical signal model. (2) Parabolic interpolation combined with linear filter interpolation of the correlation function. (3) Parabolic interpolation to the complex correlation function envelope. (4) Matched filter interpolation applied to the correlation function. The new interpolation methods are analyzed both by computer simulated signals and RF-signals recorded from a patient with time delay larger than 1/f(0), where f(0) is the center frequency. The simulation results show that these methods are more accurate than the parabolic-fit method. From the simulation, the worst estimation accuracy is about 1.25% of 1/f(0) for the parabolic-fit interpolation, and it is improved by the above methods to less than 0.5% of 1/f(0) when the sampling rate is 10 MHz, the center frequency is 2.5 MHz and the bandwidth is 1 MHz. This improvement also can be observed in the experimental data. Furthermore, the matched filter interpolation gives the best performance when signal-to-noise ratio (SNR) is low. This is verified both by simulation and experimentation.     

A novel envelope detector for high-frame rate, high-frequency ultrasound imaging

This paper proposes a novel design of envelope detectors capable of supporting a small animal cardiac imaging system requiring a temporal resolution of more than 150 frames per second. The proposed envelope detector adopts the quadrature demodulation and the lookup table (LUT) method to compute the magnitude of the complex baseband components of received echo signals. Because the direct use of the LUT method for a square root function is not feasible due to a large memory size, this paper presents a new LUT strategy dramatically reducing its size by using binary logarithmic number system (BLNS). Due to the nature of BLNS, the proposed design does not require an individual LOG-compression functional block. In the implementation using a field programmable gate array (FPGA), a total of 166.56 Kbytes memories were used for computing the magnitude of 16-bit in-phase and quadrature components instead of 4 Gbytes in the case of the direct use of the LUT method. The experimental results show that the proposed envelope detector is capable of generating LOG-compressed envelope data at every clock cycle after 32 clock cycle latency, and its maximum error is less than 0.5 (i.e., within the rounding error), compared with the arithmetic results of square root function and LOG compression.     

Use of modulated excitation signals in medical ultrasound. Part III: High frame rate imaging

This paper, the last from a series of three papers on the application of coded excitation signals in medical ultrasound, investigates the possibility of increasing the frame rate in ultrasound imaging by using modulated excitation signals. Linear array-coded imaging and sparse synthetic transmit aperture imaging are considered, and the trade-offs between frame rate, image quality, and SNR are discussed. It is shown that FM codes can be used to increase the frame rate by a factor of two without a degradation in image quality and by a factor of 5, if a slight decrease in image quality can be accepted. The use of synthetic transmit aperture imaging is also considered, and it is here shown that Hadamard spatial encoding in transmit with FM emission signals can be used to increase the frame rate by 12 to 25 times with either a slight or no reduction in signal-to-noise ratio and image quality. By using these techniques a complete ultrasound-phased array image can be created using only two emissions.     

Improved synthetic aperture focusing technique with applications in high-frequency ultrasound imaging

Synthetic aperture focusing using a virtual source was used previously to increase the penetration and to extend the depth of focus in high-frequency ultrasonic imaging. However, the performance of synthetic aperture focusing is limited by its high sidelobes. In this paper, an adaptive weighting technique based on a focusing-quality index is introduced to suppress the sidelobes. The focusing-quality index is derived from the spatial spectrum of the scan-line data along the mechanical scan direction (i.e., the synthetic aperture direction) after focusing delays relative to the virtual source have been applied. The proposed technique is of particular value in high-frequency ultrasound in which dynamic focusing using array transducers is not yet possible. Experimental ultrasound data from a 50-MHz imaging system with a single-crystal transducer (f-number=2) are used to demonstrate the efficacy of the proposed technique on both wire targets and speckle-generating objects. An in vivo experiment also is performed on a mouse to further demonstrate the effectiveness. Both 50-MHz fundamental imaging and 50-MHz tissue harmonic imaging are tested. The results clearly demonstrate the effectiveness in sidelobe reduction and background-noise suppression for both imaging modes. The principles, experimental results, and implementation issues of the new technique are described in this paper.     

Soot volume fraction imaging

A technique is described for the determination of soot volume fractions by laser light extinction measurements. This technique differs from previously reported pointwise methods in that a two-dimensional array (i.e., image) of data is acquired simultaneously. In this fashion the net data rate is increased and allows the study of time-dependent phenomena and the investigation of spatial and temporal correlations. A telecentric imaging configuration is employed to provide depth invariant magnification and to permit the specification of the collection angle for scattered light. A method is also employed to suppress undesirable coherent imaging effects. A discussion of the tomographic inversion process is also provided, including the results obtained from numerical simulations.     

Lateral RF image synthesis using a synthetic aperture imaging technique

The oscillating profile naturally present in ultrasound images has been shown to be extremely valuable in different applications, particularly in motion estimation. Recent studies have shown that it is possible to produce images with transverse oscillations (TOs) based on a specific type of beamforming. However, there is still a great difference between the nature of the lateral oscillations produced with current methods and the axial profile of ultrasound images. In this study, we propose to combine synthetic aperture imaging (synthetic transmit aperture, STA) using a specific beamformer in both transmit mode and receive mode combined with a heterodyning demodulation method to produce lateral radiofrequency signals (LRFs). The aim was to produce lateral signals as close as possible to conventional axial signals, which would make it possible to estimate lateral displacements with the same accuracy as in the axial direction. The feasibility of this approach was validated in simulation and experimentally on an ultrasound research platform, the Ultrasonix RP system. We show that the combination of STA and the heterodyning demodulation can divide the wavelength of the LRF signals by 4 and divide the width of the lateral envelope of the point spread function (PSF) by 2 compared with the previous approaches using beamforming in receive mode only. Finally, we also illustrate the potential of our beamforming for motion estimation compared with previous TO methods.     

An Experimental Data Acquisition System for Ultrasound Imaging

An experimental digital data acquisition system for computerized ultrasound imaging studies is developed. The system requirements for the digital reconstruction of tomographic images from the envelope of ultrasound signals of several megahertz are described. To effectively digitize ultrasound signals, a new digital averager is implemented. Several consecutive ultrasound signals are averaged on a real-time base to enhance the signal-to-noise ratio, to reduce the volume of the data, and to achieve a match of data flow rates between the ultrasound signals and the usual computer input devices. The positions of the ultrasound transducers and the target are controlled by a microprocessor controller. These positions are sensed and also digitized by the averager. The averager transfers the digitized data, both the ultrasound signals and position signals, to a general-purpose computer for further data processing. Experimental data on reconstructing the cross section of a simple target from the ultrasound signals acquired by this system will be presented as an illustration.     

Quantitative blood speed imaging with intravascular ultrasound

Previously, we presented a method of real-time arterial color flow imaging using an intravascular ultrasound (IVUS) imaging system, where real-time RF A-scans were processed with an FIR (finite-impulse response) filter bank to estimate relative blood speed. Although qualitative flow measurements are clinically valuable, realizing the full potential of blood flow imaging requires quantitative flow speed and volume measurements in real time. Unfortunately, the rate of RF echo-to-echo decorrelation is not directly related to scatterer speed in a side-looking IVUS system because the elevational extent of the imaging slice varies with range. Consequently, flow imaging methods using any type of decorrelation processing to estimate blood speed without accounting for spatial variation of the radiation pattern will have estimation errors that prohibit accurate comparison of speed estimates from different depths. The FIR filter bank approach measures the rate of change of the ultrasound signal by estimating the slow-time spectrum of RF echoes. A filter bank of M bandpass filters is applied in parallel to estimate M components of the slow-time DFT (discrete Fourier transform). The relationship between the slow-time spectrum, aperture diffraction pattern, and scatterer speed is derived for a simplified target. Because the ultimate goal of this work is to make quantitative speed measurements, we present a method to map slow time spectral characteristics to a quantitative estimate. Results of the speed estimator are shown for a simulated circumferential catheter array insonifying blood moving uniformly past the array (i.e., plug flow) and blood moving with a parabolic profile (i.e., laminar flow)     

Parallel beamforming using synthetic transmit beams

Parallel beamforming is frequently used to increase the acquisition rate of medical ultrasound imaging. However, such imaging systems will not be spatially shift invariant due to significant variation across adjacent beams. This paper investigates a few methods of parallel beam-forming that aims at eliminating this flaw and restoring the shift invariance property. The beam-to-beam variations occur because the transmit and receive beams are not aligned. The underlying idea of the main method presented here is to generate additional synthetic transmit beams (STB) through interpolation of the received, unfocused signal at each array element prior to beamforming. Now each of the parallel receive beams can be aligned perfectly with a transmit beam--synthetic or real--thus eliminating the distortion caused by misalignment. The proposed method was compared to the other compensation methods through a simulation study based on the ultrasound simulation software Field II. The results have been verified with in vitro experiments. The simulations were done with parameters similar to a standard cardiac examination with two parallel receive beams and a transmit-line spacing corresponding to the Rayleigh criterion, wavelength times f-number (lambda x f#). From the results presented, it is clear that straightforward parallel beamforming reduces the spatial shift invariance property of an ultrasound imaging system. The proposed method of using synthetic transmit beams seems to restore this important property, enabling higher acquisition rates without loss of image quality.     

Comparison of conventional parallel beamforming with plane wave and diverging wave imaging for cardiac applications: a simulation study

When imaging the heart, good temporal resolution is beneficial for capturing the information of short-lived cardiac phases (in particular, the isovolumetric phases). To increase the frame rate, parallel beamforming is a commonly used technique for fast cardiac imaging. Conventionally, a 4 multiple-line-acquisition (4MLA) system increases the frame rate by a factor of 4, making use of a broadened transmit beam to reduce block-like artifacts. As an alternative, it has been proposed to transmit an unfocused beam (i.e., plane wave or diverging wave) for which a large number of parallel receive beams (i.e., 16) can be formed for each transmit. However, to keep the spatial resolution acceptable in these approaches, spatial compounding of overlapping successive transmits is required. As a result, the effective gain in frame rate is similar to that of a 4MLA system. To date, it remains unclear how conventional 4MLA compares to plane-wave or diverging-wave imaging when operating at similar frame rate. The goal of this study was therefore to directly contrast the performance of these beamforming methods by computer simulation. In this study, the performance of 4 different imaging systems was investigated by quantitatively evaluating the characteristics of their beam profiles. The results showed that the conventional 4MLA and plane wave imaging were very competitive imaging strategies when operating at a similar frame rate. 4MLA performed better in the near field (i.e., 10 to 50 mm), whereas plane-wave imaging had better beam profiles in the far field (i.e., 50 to 90 mm). Although diverging-wave imaging had the poorest performance in the present study, it could potentially be improved by optimizing the settings.     

A new pulsed wave Doppler ultrasound system to measure bloodvelocities beyond the Nyquist limit

Conventional pulsed wave (PW) Doppler ultrasound systems measure blood velocities unambiguously as long as the velocities are less than the Nyquist velocity. Here, a new PW Doppler ultrasound system is described, which can measure velocities up to 4.5 times higher than the Nyquist velocity. The new PW system presented here utilizes two ultrasound carrier frequencies, such that an appropriate processing of the Doppler frequencies results in an extended mean-velocity range. With this extended mean-velocity information, the complex Doppler signal is interpolated to reconstruct aliased Doppler spectra. The interpolation algorithm has been tested by computer simulations demonstrating the capabilities and limits. The integration into a commercial Doppler scanner allowed the evaluation of the new PW system in real time. In vitro and in vivo tests demonstrate the effective reconstruction of highly aliased Doppler spectra     

A new formalism for the quantification of tissue perfusion by the destruction-replenishment method in contrast ultrasound imaging

A new formalism is presented for the destruction-replenishment perfusion quantification approach at low mechanical index. On the basis of physical considerations, best-fit methods should be applied using perfusion functions with S-shape characteristics. These functions are first described for the case of a geometry with a single flow velocity, then extended to the case of vascular beds with blood vessels having multiple flow velocity values and directions. The principles guiding the analysis are, on one hand, a linearization of video echo signals to overcome the log-compression of the imaging instrument, and, on the other hand, the spatial distribution of the transmit-receive ultrasound beam in the elevation direction. An in vitro model also is described; it was used to confirm experimentally the validity of the approach using a commercial contrast agent. The approach was implemented in the form of a computer program, taking as input a sequence of contrast-specific images, as well as parameters related to the ultrasound imaging equipment used. The generated output is either flow-parameter values computed in regions-of-interest, or parametric flow-images (e.g., mean velocity, mean transit time, mean flow, flow variance, or skewness). This approach thus establishes a base for extracting information about the morphology of vascular beds in vivo, and could allow absolute quantification provided that appropriate instrument calibration is implemented.     

Phase-aberration correction using signals from point reflectors and diffuse scatterers: measurements

A method for phase-aberration correction of phased-array images is tested using a model of near-field velocity inhomogeneities. A set of grooved room-temperature vulcanizing plates was constructed to simulate near-field aberrations encountered in clinical ultrasound imaging. As expected, large image distortion was experienced when grooved plates producing significant aberrations were placed near the surface of the array. An iterative aberration correction procedure based on cross-correlation measures between neighboring elements in a phased array, using signals reflected from diffuse scatterers, significantly reduced the effects of these aberrations, producing images nearly identical to those generated in the absence of aberrations. The results suggest that a practical phase-aberration correction system can be constructed for medical ultrasound imaging and possibly all coherent imaging systems by using a sampled aperture.     

血栓多普勒信号的多参数提取及分类

血栓的准确检测可以用于早期脑血管疾病的诊断，超声多普勒是一种无损的血栓检测技术。文章使用三种信号处理方法：传统的声谱分析法、小波分析法、renyi信息量分析法对血栓多普勒信号进行分析，提取出相应的特征参数，然后对敏感的特征参数采用反向传输(Back-Propagation，简称BP)神经网络进行分类，建立起血栓、干扰噪声和正常血流信号的自动判别系统。通过对300例仿真多普勒信号和163例临床采集的大脑中动脉多普勒信号进行分析，结果表明：本文建立的系统对血栓的检测率高于传统的方法，有望可用于血栓多普勒信号的自动检测。