Resolution of synthetic-aperture imaging through turbulence

A theory is developed for the resolution of an optical synthetic-aperture imaging system viewing an object through an inhomogeneous refractive medium. The inhomogeneities of the propagation medium create errors in the phase history data with resultant space-variant image effects, including geometric distortions and broadening of the impulse response or point-spread function. I relate the intensity-impulse response to the usual wave structure function. I determine the modulation transfer function for synthetic apertures of any size and exposure time, valid whenever the optical bandwidth is small compared with the carrier frequency, and derive the resolution for monostatic and bistatic synthetic apertures, valid whenever the real sampling aperture is small compared with the medium's coherence length. The results take the same form as the well-known turbulence-limited resolution of incoherent, real-aperture imaging with short exposure. Turbulence-limited synthetic-aperture resolution is somewhat better than incoherent real-aperture resolution under the same conditions. Autofocus processing improves synthetic-aperture resolution beyond this limit, and adaptive correction of higher-order phase history errors would improve it further.     

Imaging with concave large-aperture therapeutic ultrasound arrays using conventional synthetic-aperture beamforming

Several dual-mode ultrasound array (DMUA) systems are being investigated for potential use in image- guided surgery. In therapeutic mode, DMUAs generate pulsed or continuous-wave (CW) high-intensity focused ultrasound (HIFU) beams capable of generating localized therapeutic effects within the focal volume. In imaging mode, pulse-echo data can be collected from the DMUA elements to obtain B-mode images or other forms of feedback on the state of the target tissue before, during, and after the application of the therapeutic HIFU beam. Therapeutic and technological constraints give rise to special characteristics of therapeutic arrays. Specifically, DMUAs have concave apertures with low f-number values and are typically coarsely sampled using directive elements. These characteristics necessitate pre- and post-beamforming signal processing of echo data to improve the spatial and contrast resolution and maximize the image uniformity within the imaging field of view (IxFOV). We have recently developed and experimentally validated beamforming algorithms for concave large-aperture DMUAs with directive elements. Experimental validation was performed using a 1 MHz, 64-element, concave spherical aperture with 100 mm radius of curvature. The aperture was sampled in the lateral direction using elongated elements 1-lambda x 33.3-lambda with 1.333-lambda center-to-center spacing (lambda is the wavelength). This resulted in f-number values of 0.8 and 2 in the azimuth and elevation directions, respectively. In this paper, we present a new DMUA design approach based on different sampling of the shared concave aperture to improve image quality while maintaining therapeutic performance. A pulse-wave (PW) simulation model using a modified version of the Field II program is used in this study. The model is used in generating pulse-echo data for synthetic-aperture (SA) beamforming for forming images of a variety of targets, e.g., wire arrays and speckle-generating cyst phantoms. To provide validation for the simulation model and illustrate the improvements in image quality, we show SA images of similar targets using pulse-echo data acquired experimentally using our existing 64-element prototype. The PW simulation model is used to investigate the effect of transducer bandwidth as well as finer sampling of the concave DMUA aperture on the image quality. The results show that modest increases in the sampling density and transducer bandwidth result in significant improvement in spatial and contrast resolutions in addition to extending the DMUA IxFOV.     

Optimum synthetic-aperture imaging of extended astronomical objects

In optical aperture-synthesis imaging of stellar objects, different beam combination strategies are used and proposed. Coaxial Michelson interferometers are very common and a homothetic multiaxial interferometer is recently realized in the Large Binocular Telescope. Laboratory experiments have demonstrated the working principles of two new approaches: densified pupil imaging and wide field-of-view (FOV) coaxial imaging using a staircase-shaped mirror. We develop a common mathematical formulation for direct comparison of the resolution and noise sensitivity of these four telescope configurations for combining beams from multiple apertures for interferometric synthetic aperture, wide-FOV imaging. Singular value decomposition techniques are used to compare the techniques and observe their distinct signal-to-noise ratio behaviors. We conclude that for a certain chosen stellar object, clear differences in performance of the imagers are identifiable.     

存在反射面时的虚拟环绕声系统鲁棒性研究

1引言虚拟环绕声系统在当今的多媒体计算机、家庭影院、音乐的环绕声等声重放系统中被广泛利用,实现环绕声系统的关键是串音消除系统(CCS)的设计。CCS的传输特性取决于扬声器到双耳的传输函数,因而严格来说对一定的CCS,只能在一个特定     

Computer-aided B-mode ultrasound diagnosis of hepatic steatosis: a feasibility study

Fatty liver (steatosis) occurs in obese patients, among others, and is related to the development of diabetes type-2. Timely diagnosis of steatosis is therefore of great importance. Steatosis is also the most common liver disease of high-yielding dairy cattle during early lactation. This makes it a suitable animal model for studying liver steatosis. Furthermore, reference of derived ultrasound parameters against a "gold standard" is possible in cattle by taking a liver biopsy for the assessment of fat concentration. The authors undertook this pilot study to investigate the hypothesis that quantitative, computer-aided B-mode ultrasound enables the noninvasive detection of hepatic steatosis. Echographic images were obtained postpartum from dairy cows (n = 12) in transcutaneous and direct (intraoperative) applications using a convex array transducer at 4.2 MHz. During surgery, a biopsy was taken from the caudate lobe to assess the liver fat content (fat score). A custom-designed software package for computer-aided ultrasound diagnosis (CAUS) was developed. After linearizing the post-processing look-up-table (LUT), the image gray levels were transferred into echo levels in decibels relative to the mean echo level in a tissue-mimicking phantom. The quantitative comparison of transcutaneous and intraoperative images enabled the correction for the attenuation effect of skin and subcutaneous fat layer on the mean echo level in the liver, as well as for the effects of the beam formation and attenuation of liver tissue on the echo level vs. depth. The residual attenuation coefficient (dB/cm) in fatty liver vs. normal liver was estimated and compensated for. Finally, echo level was estimated relative to the phantom used for calibration, and echo texture was characterized by the mean axial and lateral speckle size within the regions of interest. In the no fat/low fat group (n = 5) skin plus fat layer attenuation was 3.4 dB/cm. A correlation of skin layer thickness vs. fat score of r = 0.48 was found. The mean transcutaneous liver tissue echo level correlated well with fat score: r = 0.80. A residual liver attenuation coefficient of 0.76 dB/cm and 1.19 dB/cm was found in medium and high fat liver, respectively. In transcutaneous images, correlation of residual attenuation coefficient with fat score was r = 0.69. Axial and lateral speckle sizes were on the order of 0.2 and 1.0 cm, respectively, and no correlation was found with liver fat content. Results for transcutaneous and intraoperative images were similar. The authors conclude that this pilot study shows the feasibility of calibrated, computer-aided ultrasound for noninvasively diagnosing, possibly even screening, steatosis of the liver.     

High-frequency rapid B-mode ultrasound imaging for real-time monitoring of lesion formation and gas body activity during high-intensity focused ultrasound ablation

The goal of this study was to examine the ability of high-frame-rate, high-resolution imaging to monitor tissue necrosis and gas-body activities formed during high-intensity focused ultrasound (HIFU) application. Ex vivo porcine cardiac tissue specimens (n = 24) were treated with HIFU exposure (4.33 MHz, 77 to 130 Hz pulse repetition frequency (PRF), 25 to 50% duty cycle, 0.2 to 1 s, 2600 W/cm(2)). RF data from Bmode ultrasound imaging were obtained before, during, and after HIFU exposure at a frame rate ranging from 77 to 130 Hz using an ultrasound imaging system with a center frequency of 55 MHz. The time history of changes in the integrated backscatter (IBS), calibrated spectral parameters, and echo-decorrelation parameters of the RF data were assessed for lesion identification by comparison against gross sections. Temporal maximum IBS with +12 dB threshold achieved the best identification with a receiver-operating characteristic (ROC) curve area of 0.96. Frame-to-frame echo decorrelation identified and tracked transient gas-body activities. Macroscopic (millimetersized) cavities formed when the estimated initial expansion rate of gas bodies (rate of expansion in lateral-to-beam direction) crossed 0.8 mm/s. Together, these assessments provide a method for monitoring spatiotemporal evolution of lesion and gas-body activity and for predicting macroscopic cavity formation.     

Four-dimensional ultrasound current source density imaging of a dipole field

Ultrasound current source density imaging (UCSDI) potentially transforms conventional electrical mapping of excitable organs, such as the brain and heart. For this study, we demonstrate volume imaging of a time-varying current field by scanning a focused ultrasound beam and detecting the acoustoelectric (AE) interaction signal. A pair of electrodes produced an alternating current distribution in a special imaging chamber filled with a 0.9% NaCl solution. A pulsed 1 MHz ultrasound beam was scanned near the source and sink, while the AE signal was detected on remote recording electrodes, resulting in time-lapsed volume movies of the alternating current distribution.     

医学超声成像的模拟研究

简述医学超声成像模拟研究的主要方法,探讨其研究意义和应用前景。将现有的超声成像模拟算法分为基于线性假设和基于非线性假设的两大类,介绍每类中代表方法的基本原理,并讨论两类方法的差异,其中基于数值方法求解超声波动方程的模拟算法更符合实际超声成像过程,为更好地理解超声图像中的各种现象、寻求组织特性和图像特征间的匹配关系提供了有效的研究工具。     

Fabrication and characterization of transducer elements in two-dimensional arrays for medical ultrasound imaging

Some of the problems of developing a two-dimensional (2-D) transducer array for medical imaging are examined. The fabrication of a 2-D array material consisting of lead zirconate titanate (PZT) elements separated by epoxy is discussed. Ultrasound pulses and transmitted radiation patterns from individual elements in the arrays are measured. A diffraction theory for the continuous wave pressure field of a 2-D array element is generalized to include both electrical and acoustical cross-coupling between elements. This theory can be fit to model the measured radiation patterns of 2-D array elements, giving an indication of the level of cross-coupling in the array, and the velocity of the acoustic cross-coupling wave. Improvements in bandwidth and cross-coupling resulting from the inclusion of a front acoustic matching layer are demonstrated, and the effects of including a lossy backing material on the array are discussed. A broadband electrical matching network is described, and pulse-echo waveforms and insertion loss from a 2-D array element are measured.     

The design of protection circuitry for high-frequency ultrasound imaging systems

Transmission line lengths in the protection circuitry of a high-frequency (>20-MHz) ultrasound imaging system have an important effect on the frequency, amplitude, and bandwidth of the pulse-echo response of the system. A model that includes the transmission line lengths between the pulser, transducer, and receiver and the electromechanical properties of high-frequency transducers is used to illustrate the importance of correctly choosing these line lengths. An iterative optimization procedure for designing the protection circuitry for a broadband system is proposed. A theoretical and experimental analysis of the validity of this approach is reported for a 45-MHz PVDF transducer.     

Cystic resolution: a performance metric for ultrasound imaging systems

This paper describes a metric that can be used to characterize the resolution of arbitrary broadband coherent imaging systems. The metric is particularly suited to medical ultrasound because it characterizes scanner performance using the contrast obtained by imaging anechoic cysts of various sizes that are embedded in a speckle-generating background, accounting for the effect of electronic noise. We present the theoretical derivation of the metric and provide simulation examples that demonstrate its utility. We use the metric to compare a low-cost, handheld, C-scan system under development in our laboratory to conventional ultrasound scanners. We also present the results of simulations that were designed to evaluate and optimize various parameters in our system, including the f/# and apodization windows. We investigate the impact of electronic noise on our system and quantify the tradeoffs associated with quantization in the analog to digital converter. Results indicate that an f/1 receive aperture combined with 10-bit precision and a signal-to-noise ratio (SNR) of 0 dB per channel would result in adequate image quality.     

Cardiac activation mapping using ultrasound current source density imaging (UCSDI)

We describe the first mapping of biological current in a live heart using ultrasound current source density imaging (UCSDI). Ablation procedures that treat severe heart arrhythmias require detailed maps of the cardiac activation wave. The conventional procedure is time-consuming and limited by its poor spatial resolution (5-10 mm). UCSDI can potentially improve on existing mapping procedures. It is based on a pressure-induced change in resistivity known as the acousto-electric (AE) effect, which is spatially confined to the ultrasound focus. Data from 2 experiments are presented. A 540 kHz ultrasonic transducer (f/# = 1, focal length = 90 mm, pulse repetition frequency = 1600 Hz) was scanned over an isolated rabbit heart perfused with an excitation-contraction decoupler to reduce motion significantly while retaining electric function. Tungsten electrodes inserted in the left ventricle recorded simultaneously the AE signal and the low-frequency electrocardiogram (ECG). UCSDI displayed spatial and temporal patterns consistent with the spreading activation wave. The propagation velocity estimated from UCSDI was 0.25 plusmn 0.05 mm/ms, comparable to the values obtained with the ECG signals. The maximum AE signal-to-noise ratio after filtering was 18 dB, with an equivalent detection threshold of 0.1 mA/ cm². This study demonstrates that UCSDI is a potentially powerful technique for mapping current flow and biopotentialsin the heart.     

Synthetic aperture techniques with a virtual source element

A new imaging technique has been proposed that combines conventional B-mode and synthetic aperture imaging techniques to overcome the limited depth of field for a highly focused transducer. The new technique improves lateral resolution beyond the focus of the transducer by considering the focus a virtual element and applying synthetic aperture focusing techniques. In this paper, the use of the focus as a virtual element is examined, considering the issues that are of concern when imaging with an array of actual elements: the tradeoff between lateral resolution and sidelobe level, the tradeoff between system complexity (channel count/amount of computation) and the appearance of grating lobes, and the issue of signal to noise ratio (SNR) of the processed image. To examine these issues, pulse-echo RF signals were collected for a tungsten wire in degassed water, monofilament nylon wires in a tissue-mimicking phantom, and cyst targets in the phantom. Results show apodization lowers the sidelobes, but only at the expense of lateral resolution, as is the case for classical synthetic aperture imaging. Grating lobes are not significant until spatial sampling is more than one wavelength, when the beam is not steered. Resolution comparable to the resolution at the transducer focus can be achieved beyond the focal region while obtaining an acceptable SNR. Specifically, for a 15-MHz focused transducer, the 6-dB beamwidth at the focus is 157 mum, and with synthetic aperture processing the 6-dB beamwidths at 3, 5, and 7 mm beyond the focus are 189 mum, 184 mum, and 215 mum, respectively. The image SNR is 38.6 dB when the wire is at the focus, and it is 32.8 dB, 35.3 dB, and 38.1 dB after synthetic aperture processing when the wire is 3, 5, and 7 mm beyond the focus, respectively. With these experiments, the virtual source has been shown to exhibit the same behavior as an actual transducer element in response to synthetic aperture processing techniques.     

Position and time-delay calibration of transducer elements in a sparse array for underwater ultrasound imaging

This paper describes a novel method for the calibration of the position and time delay of transducer elements in a large, sparse array used for underwater, high-resolution ultrasound imaging. This method is based on the principles used in the global positioning system (GPS). However, unlike GPS, in which the wave propagation speed is generally assumed known, the sound propagation speed in the water usually is unknown and it is calibrated simultaneously in this method to achieve high calibration accuracy. In this method, a high-precision positioning system is used to scan a single hydrophone (used for transmission) in the imaging field of the array. The hydrophone transmits pulses at selected positions, and the transducer elements in the sparse array receive the transmitted signals. Time of flight (TOF) values between transducer elements and hydrophone positions then are measured. From a series of measured TOF values, the position and time delay values for each transducer element as well as the propagation speed can be calibrated. The performances of the calibration algorithm are theoretically analyzed and evaluated with numerical calculations and simulation studies. It is found that this method is capable of calibrating the positions and time delays of transducer elements with high accuracy.     

A study of the influence of damage in equitype elements on vibration of regular systems

The results of an experimental study of the effect of operation-induced damage in equitype elements on vibration of regular system are illustrated by the example of testing a fork type specimen. The in-service frequency detuning of equitype elements is found to have an influence on the formation of vibrations of regular systems similar to their technological detuning. __________ Translated from Problemy Prochnosti, No. 2, pp. 39–46, March–April, 2006.     

Vibrations in systems with dry-friction damping elements

Forced vibrations of a mechanical system having one degree of freedom with dry and viscous friction are examined. It is established that dry friction acts as a filter that allows the passage of forced vibrations of a precisely determined frequency range in relation to the ratio of the disturbing force amplitude to amount of dry friction. Addition of viscous friction to dry friction narrows the vibration-passing sections of this filter.Translated from Problemy Prochnosti, No. 4, pp. 33–36, April, 1992.     

Elasticity imaging using conventional and high-frame rate ultrasound imaging: experimental study

High-frame rate ultrasound imaging is necessary to track fast deformation in ultrasound elasticity imaging, but the image quality may be degraded. Previously, we investigated the performance of strain imaging using numerical models of conventional and ultrafast ultrasound imaging techniques. In this paper, we performed experimental studies to quantitatively evaluate the strain images and elasticity maps obtained using conventional and high frame rate ultrasound imaging methods. The experiments were carried out using point target and tissue mimicking phantoms. The experimental results were compared with the results of numerical simulation. Our experimental studies confirm that the signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and axial/lateral resolution of the displacement and strain images acquired using high-frame rate ultrasound imaging are slightly lower but comparable with those obtained using conventional imaging. Furthermore, the quality of elasticity images also exhibits similar trends. Thus, high-frame rate ultrasound imaging can be used reliably for static elasticity imaging to capture the internal tissue motion if the frame rate is critical.     

一种新的超声成像方法振动声成像

文章介绍一种新的超声成像方法———振动声成像。首先从理论上分析了振动声成像的基本原理;然后分别介绍了该方法在钙化动脉血管成像、乳房中微钙化灶的检测以及受热疗和HIFU治疗的肿瘤组织焦斑检测等方面的应用;最后比较了振动声成像与传统超声成像不同点并指出了振动声成像独特的优势。     

A Fourier transform-based sidelobe reduction method in ultrasound imaging

Focusing is widely used to increase the resolution in medical ultrasound imaging systems. Focusing increases signal levels returning from the mainlobe direction and decreases those from sidelobe directions. The sidelobes, when not completely cancelled, deteriorate the resulting image resolution. This paper proposes a method of improving the resolution by scaling the received signal according to the ratio between the mainlobe and the sidelobe levels computed in the frequency domain by the use of Fourier transform. The proposed method is verified by computer simulation and experiment and is shown to be highly effective in narrowing the mainlobe width and decreasing the sidelobe levels at the same time.     

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)