In vitro estimation of mean sound speed based on minimum average phase variance in medical ultrasound imaging

Effective receive beamforming in medical ultrasound imaging is important for enhancing spatial and contrast resolution. In current ultrasound receive beamforming, a constant sound speed (e.g., 1540 m/s) is assumed. However, the variations of sound speed in soft tissues could introduce phase distortions, leading to degradation in spatial and contrast resolution. This degradation becomes even more severe in imaging fatty tissues (e.g., breast) and with obese patients. In this paper, a mean sound speed estimation method where phase variance of radio-frequency channel data in the region of interest is evaluated is presented for improving spatial and contrast resolution. The proposed estimation method was validated by the Field II simulation and the tissue mimicking phantom experiments. In the simulation, the sound speed of the medium was set to 1450 m/s and the proposed method was capable of capturing this value correctly. From the phantom experiments, the −18-dB lateral resolution of the point target at 50 mm obtained with the estimated mean sound speed was improved by a factor of 1.3, i.e., from 3.9 mm to 2.9 mm. The proposed estimation method also provides an improvement of 0.4 in the contrast-to-noise ratio, i.e., from 2.4 to 2.8. These results indicate that the proposed mean sound speed estimation method could enhance the spatial and contrast resolution in the medical ultrasound imaging systems.     

Compensation of fat layer effects in ultrasound imaging using an inclined-fat-layer model derived from magnetic resonance images

When cutaneous fat layers are in the ultrasound imaging region, the phase aberration caused by the fat layers induce image distortion as well as spatial resolution degradation. The phase aberration may complicate clinical procedures particularly when ultrasound imaging is employed for spatial positioning of medical devices like a biopsy needle or HIFU. To compensate the fat layer effects more precisely in beamforming, an inclined-fat-layer model has been established from the magnetic resonance images of the same imaging region as in the ultrasound scanning. We have verified utility of the fat layer model by taking images of a metal needle put into an inclined-fat-layer mimicking phantom. The ultrasound images taken with a 128-element linear phase array operating at 6 MHz have shown better resolution and less distortion when receive beamforming was performed with the phase delay data derived from the inclined-fat-layer model.     

Gaussian wavelet based dynamic filtering (GWDF) method for medical ultrasound systems

In this paper, a novel dynamic filtering method using Gaussian wavelet filters is proposed to remove noise from ultrasound echo signal. In the proposed method, a mother wavelet is first selected with its central frequency (CF) and frequency bandwidth (FB) equal to those of the transmitted signal. The actual frequency of the received signal at a given depth is estimated through the autocorrelation technique. Then the mother wavelet is dilated using the ratio between the transmitted central frequency and the actual frequency as the scale factor. The generated daughter wavelet is finally used as the dynamic filter at this depth. Frequency-demodulated Gaussian wavelet is chosen in this paper because its power spectrum is well-matched with that of the transmitted ultrasound signal. The proposed method is evaluated by simulations using Field II program. Experiments are also conducted out on a standard ultrasound phantom using a 192-element transducer with the center frequency of 5 MHz. The phantom contains five point targets, five circular high scattering regions with diameters of 2, 3, 4, 5, 6 mm respectively, and five cysts with diameters of 6, 5, 4, 3, 2 mm respectively. Both simulation and experimental results show that optimal signal-to-noise ratio (SNR) can be obtained and useful information can be extracted along the depth direction irrespective of the diagnostic objects.     

Design and experiment of 256-element ultrasound phased array for noninvasive focused ultrasound surgery

A 256-element phased array high intensity focused ultrasound system has been designed and constructed in our laboratory. The 256-element spherical-section ultrasound phased array made from piezocomposite material operates at 1.1 MHz with 11-cm radius of curvature, 14-cm outer diameter, and 3.4-cm diameter central hole for mounting diagnostic ultrasound imaging probe. First, the explicit sound field calculation approach for the spherical-section phased array and the genetic optimal algorithm are briefly introduced as the optimal focus pattern control methods. Then, the design guidelines of 256-element spherical-section focused ultrasound phased array and 256-channel driver system are given. The results of single on and off axial focus, multiple on and off axial foci, half sub-array focus pattern provide further evidence for the 3D focus steering and sub-array mode for avoiding obstacles in focused ultrasound surgery. The multi-foci pattern can enlarge the treatment volume to 22 times larger than that of a single focus in one sonication. Finally, in vitro and transparent tissue-mimicking phantom experiment results confirm the ability of 256-element spherical-section phased array system to induce thermal lesions for noninvasive ultrasound surgery.     

Tissue velocity imaging of coronary artery by rotating-type intravascular ultrasound

Intravascular ultrasound (IVUS) provides not only the dimensions of coronary artery but the information of tissue components. In catheterization laboratory, soft and hard plaques are classified by visual inspection of echo intensity. So-called soft plaque contains lipid core or thrombus and it is believed to be more vulnerable than a hard plaque. However, it is not simple to analyze the echo signals quantitatively. When we look at a reflection signal, the intensity is affected by the distance of the object, the medium between transducer and objects and the fluctuation caused by rotation of IVUS probe. The time of flight is also affected by the sound speed of the medium and Doppler shift caused by tissue motion but usually those can be neglected. Thus, the analysis of RF signal in time domain can be more quantitative than intensity of RF signal. In the present study, a novel imaging technique called "intravascular tissue velocity imaging" was developed for searching a vulnerable plaque. Radio-frequency (RF) signal from a clinically used IVUS apparatus was digitized at 500 MSa/s and stored in a workstation. First, non-uniform rotation was corrected by maximizing the correlation coefficient of circumferential RF signal distribution in two consecutive frames. Then, the correlation and displacement were calculated by analyzing the radial difference of RF signal. Tissue velocity was determined by the displacement and the frame rate. The correlation image of normal and atherosclerotic coronary arteries clearly showed the internal and external borders of arterial wall. Soft plaque with low echo area in the intima showed high velocity while the calcified lesion showed the very low tissue velocity. This technique provides important information on tissue character of coronary artery.     

Sound speed correction in ultrasound imaging

A constant sound speed of 1.54 mm/micros is generally used by ultrasound imaging systems for delay and timing. However, the body's sound speed in-homogeneity can lead to defocusing and increased clutter. To provide an improvement using standard transducers, the sound speed used in delay and timing was computed using different sound speeds. We observed improvement in lateral resolution and clutter in phantom, OB, abdominal, and breast imaging. In OB and abdominal imaging using a 4 MHz curved array, 1.48 mm/micros provided higher image quality in many situations. In breast with an 8 MHz linear array, 1.44 mm/micros provided better images in some cases. To provide an automated way to determine and adjust the sound speed used by the imaging system, an algorithm was developed that determines the sound speed that produces the best overall lateral image quality by analyzing the spatial frequency content in a single B-mode frame of channel data using images reconstructed using various trial sound speeds. The metric produced correlates well with the observed best lateral image quality.     

A novel power spectrum calculation method using phase-compensation and weighted averaging for the estimation of ultrasound attenuation

An estimation of ultrasound attenuation in soft tissues is critical in the quantitative ultrasound analysis since it is not only related to the estimations of other ultrasound parameters, such as speed of sound, integrated scatterers, or scatterer size, but also provides pathological information of the scanned tissue. However, estimation performances of ultrasound attenuation are intimately tied to the accurate extraction of spectral information from the backscattered radiofrequency (RF) signals. In this paper, we propose two novel techniques for calculating a block power spectrum from the backscattered ultrasound signals. These are based on the phase-compensation of each RF segment using the normalized cross-correlation to minimize estimation errors due to phase variations, and the weighted averaging technique to maximize the signal-to-noise ratio (SNR). The simulation results with uniform numerical phantoms demonstrate that the proposed method estimates local attenuation coefficients within 1.57% of the actual values while the conventional methods estimate those within 2.96%. The proposed method is especially effective when we deal with the signal reflected from the deeper depth where the SNR level is lower or when the gated window contains a small number of signal samples. Experimental results, performed at 5 MHz, were obtained with a one-dimensional 128 elements array, using the tissue-mimicking phantoms also show that the proposed method provides better estimation results (within 3.04% of the actual value) with smaller estimation variances compared to the conventional methods (within 5.93%) for all cases considered.     

The-effective directivity characteristic of a pulsed ultrasound transducer and its measurement by semi-automatic means

A simple electronic apparatus for plotting the effective directivity characteristics of pulsed, or continuously excited, ultrasound (1–10 MHz) transducers has been built. Beam profiles are plotted by recording a dc signal, derived from the peak amplitude of either the echo from a steel ball or the signal from a small hydrophone, as a function of position in the acoustic field. An additional capability of direct plotting of isoamplitude contours provides a particularly easy and informative means of recording and displaying the acoustic field.The apparatus has been used to investigate the relationship, in grey scale imaging systems, between dynamic range compression characteristics and effective display resolution. It is shown that, so long as the gain characteristics of the system remain constant, the display resolution capability is inversely related to the strength of the scatterers to be visualized.     

Non-contact sound speed measurement by optical probing of beam deflection due to sound wave

We report a non-contact and non-invasive method of sound speed measurement by optical probing of deflected laser beam due to normally incident degenerated shock wave. In this study the shock wave from an exploding wire was degenerated to an ordinary sound wave at the distance exceeding 0.23 m. Temporal resolution of the deflected beam signal was improved by passing through an adequate electronic high-pass filter, as a result we obtained a better temporal resolution than that of the acoustic pressure detection by PZT transducer in terms of rising time. The spatial resolution was improved by passing the refracted beam signal into the edge of focusing lens to make a larger deflection angle. Sound speed was calculated by monitoring the time of flight of transient deflected signal at the predetermined position. Sound speed has been measured in air, distilled water and acryl, agreed well with the published values. The sound speed measured in the solution of glycerin, magnesium sulfate (MgSO4), and dimethylformamide with various mole fractions also agrees within 3% of relative error with those measured in the present work by ultrasonic pulse echo method. The results suggest that the method proposed is to be reliable and reproducible.     

Lens multielement acoustic microscope in the mode for measuring the parameters of layered objects

An acoustic microscope with a cylindrical lens and ultrasound transducer have been considered, as well as the method based on it for the measuring of longitudinal and transverse wave velocities, the thickness and density of the investigated layer. A theoretical model of the microscope has been constructed, and the relation between the spatiotemporal output signal of the transducer and the angular dependence of the sample reflection coefficient has been found. It has been shown that the velocities of body waves and the thickness can be determined by the delays of ultrasound responses reflected from the layer boundaries measured by the transducer elements, and the density, by the amplitudes of these responses. The method was tested experimentally using a 20-element transducer with a central frequency of 15 MHz and a period of 0.8 mm. The example of a duralumin plate has shown that the error in measuring the thickness and velocity of longitudinal waves error does not exceed 1%; the velocity of transverse waves, 2%; and the density can be estimated with an accuracy of about 5%.     

压缩感知在合成发射孔径医学超声成像中的应用

为了解决合成发射孔径技术在医学超声成像实现中面临的数据量大及接收通道多的问题,提出一种超声成像系统频率域稀疏性模型的压缩感知成像算法。首先对超声系统频率域稀疏性模型进行了验证;然后根据稀疏性模型利用压缩感知理论对回波信号进行压缩采样,并使用最优化方法完成回波信号重建;最终通过合成发射孔径技术完成超声成像。针对医学成像中常用的点目标及模拟胎儿目标进行成像仿真实验,对重建图像在均方误差、分辨率及成像质量等方面与常规成像结果对比分析。实验结果表明在保证成像质量的同时,仅使用30%原始数据量及50%总接收通道数目可完成成像;频率域稀疏性模型的压缩感知成像算法可以大幅度减少合成发射孔径成像所需数据量及接收通道数,极大地降低了系统复杂度。      

A digital image analysis method for diagnostic ultrasound calibration

Acoustic test objects are commonly used for quality assurance testing of diagnostic ultrasound machines. However, the accompanying calibration protocols rely heavily on the judgment of the sonographer, are dependent on machine settings and are semi-quantitative. In the current study, two unique test objects and protocols were designed to quantitatively determine diagnostic ultrasound parameters, namely axial resolution and geometric uniformity, and lateral resolution and geometric uniformity of the ultrasound field. The effect of focal zone, signal gain, and distance from the ultrasound probe on these parameters was assessed. The investigation was performed using a typical low-frequency diagnostic unit equipped with a 7.5 MHz linear pulse–echo probe. Results underline the need to ensure that sensitivity of routine testing regimes is adequate for the measurements to be made. This study is a preliminary part of a larger project developing an ultrasound technique to be used as an engineering design tool in a non-clinical industrial setting for quality assurance testing of total knee replacements immersed in water.     

Sensitivity to point-spread function parameters in medical ultrasound image deconvolution

Clinical ultrasound images are often perceived as difficult to interpret due to image blurring and speckle inherent in the ultrasound imaging. But the image quality can be improved by deconvolution using an estimate of the point-spread function. However, it is difficult to obtain a sufficiently accurate estimate of the point-spread function in vivo because of the unknown properties of the soft tissue in clinical applications. Local variations in the speed of sound and attenuation change the pulse and beam shape. These in turn affect the point-spread function. The purpose and novelty of this paper is therefore to explore the sensitivity of a state-of-the-art deconvolution algorithm to uncertainty in the point-spread function. The point-spread function in our restoration algorithm is made shift invariant in the lateral dimension but shift dependent in the axial direction, and is modelled to match a 128-element 1D linear array often found in clinical use. We present simulated and in vitro sensitivity analyses of two-dimensional deconvolution while varying six parameters on which the point-spread function depends. Uncertainty in the ultrasound machine is analysed by varying the axial depths of lateral and elevational foci alongside height and width of transducer elements. Sensitivity to tissue influence is investigated by varying the speed of sound and frequency-dependent attenuation of the electro-mechanical impulse response. The results are analysed both quantitatively and in terms of the perceived image quality. First, the assessment of deconvolution using the logarithmic image amplitude is found to be a better indicator of the perceived improvement in the restoration. Secondly, the two most critical parameters for two-dimensional deconvolution are discovered to be the lateral focus and the speed of sound, because the success of deconvolution is perceived primarily in terms of deblurring. We also observed similar patterns for the simulation and in vitro experiment. Finally, we show that it is possible to restore in vivo ultrasound images using an assumed point-spread function and hence conclude that an exact point-spread function is not necessary for enhancing ultrasound image quality by deconvolution.     

Reflectivity function reconstruction using adaptive lattice deconvolution technique

Deconvolution of ultrasonic echo signals improves resolution and quality of ultrasonic images. The problem of reconstructing the reflectivity of a biological tissue is examined by adaptive lattice deconvolution of the echo ultrasound signals. The simulation of the signal formation process in an ultrasonic-echo scan line in noisy conditions is estimated. The reflectivity of a biological tissue is estimated as cross-correlation coefficients between forward and backward predication errors in each stage of the adaptive lattice filter.     

经颅超声刺激人脑海马的数值仿真研究

颞叶癫痫、阿尔兹海默症等神经系统疾病与大脑海马的异常放电相关。近年兴起的经颅超声刺激治疗脑疾病具有无创、能深入脑组织和可调控海马神经元放电的特点。该文基于128阵元相控阵换能器、人头颅CT数据、水体建立超声经颅刺激海马数值仿真模型,数值仿真优选适应于海马刺激的换能器结构参数,并结合电容模型,探究超声声学参数对海马神经元放电的影响。结果表明曲率半径90 mm、开口半径56 mm、阵元半径2.0 mm、频率0.9 MHz的128阵元相控超声换能器可形成适应于刺激人脑海马大小的焦域;且在超声频率为0.9 MHz时对海马神经元放电有较好的抑制作用;在较小的占空比和较低的空间峰值时间平均声强条件下对神经元放电有较好的抑制作用。     

Influence of the ultrasound transducer bandwidth on selection of the complementary Golay bit code length

In contrast to previously published papers [A. Nowicki, Z. Klimonda, M. Lewandowski, J. Litniewski, P.A. Lewin, I. Trots, Comparison of sound fields generated by different coded excitations – Experimental results, Ultrasonics 44 (1) (2006) 121–129; J. Litniewski, A. Nowicki, Z. Klimonda, M. Lewandowski, Sound fields for coded excitations in water and tissue: experimental approach, Ultrasound Med. Biol. 33 (4) (2007) 601–607], which examined the factors influencing the spatial resolution of coded complementary Golay sequences (CGS), this paper investigates the effect of ultrasound imaging transducer’s fractional bandwidth on the gain of the compressed echo signal for different spectral widths of the CGS. Two different bit lengths were considered, specifically one and two cycles. Three transducers having fractional bandwidth of 25%, 58% and 80% and operating at frequencies 6, 4.4 and 6 MHz, respectively were examined (one of the 6 MHz sources was focused and made of composite material). The experimental results have shown that by increasing the code length, i.e. decreasing the bandwidth, the compressed echo amplitude could be enhanced. The smaller the bandwidth was the larger was the gain; the pulse-echo sensitivity of the echo amplitude increased by 1.88, 1.62 and 1.47, for 25%, 58% and 80% bandwidths, respectively. These results indicate that two cycles bit length excitation is more suitable for use with bandwidth limited commercially available imaging transducers. Further, the time resolution is retained for transducers with two cycles excitation providing the fractional bandwidth is lower than approximately 90%. The results of this work also show that adjusting the code length allows signal-to-noise-ratio (SNR) to be enhanced while using limited (less that 80%) bandwidth imaging transducers. Also, for such bandwidth limited transducers two cycles excitation would not decrease the time resolution, obtained with “conventional” spike excitation. Hence, CGS excitation could be successfully implemented with the existing, relatively narrow band imaging transducers without the need to use usually more expensive wideband, composite ones.     

Estimation of speed of sound in dual-layered media using medical ultrasound image deconvolution

The speed of sound in soft tissues is usually assumed to be 1540 m/s in medical pulse-echo ultrasound imaging systems. When the true speed is different, the mismatch can lead to distortions in the acquired images, and so reduce their clinical value. Previously we reported a new method of sound-speed estimation in the context of image deconvolution. Unlike most other sound-speed estimation methods, this enables the use of unmodified ultrasound machines and a normal scanning pattern. Our approach was validated for largely homogeneous media with single sound speeds. In this article, we demonstrate that sound speeds of dual-layered media can also be estimated through image deconvolution. An ultrasound simulator has been developed for layered media assuming that, for moderate speed differences, the reflection at the interface may be neglected. We have applied our dual-layer algorithm to simulations and in vitro phantoms. The speed of the top layer is estimated by our aforesaid method for homogeneous media. Then, when the layer boundary position is known, a series of deconvolutions are carried out with dual-layered point-spread functions having different lower-layer speeds. The best restoration is selected using a correlation metric. The error level (e.g., a mean error of −9 m/s with a standard deviation of 16 m/s) for in vitro phantoms is found to be not as good as that of our single-speed algorithm, but is comparable to other local speed estimation methods where the data acquisition may not be as simple as in our proposed method.     

圆合成孔径声呐多点定位运动补偿

圆合成孔径声呐(CSAS)的成像性能受平台运动误差影响而下降,利用单侧回波可估计CSAS基阵的斜距误差,但单侧回波在小测绘带时无法估计升沉误差,针对此问题,提出了一种利用单侧回波信号的声呐平台三维运动估计和补偿方法。首先,对CSAS在不同观测角度的目标回波取极大值获得目标回波的到达时间;其次,基于多个点目标的到达时间建立CSAS三维定位模型;然后利用列文伯格-马夸尔特方法对声呐三维坐标进行估计;最后将位置估计结果与时域反投影成像方法结合实现对目标的成像.仿真结果表明:该方法能精确估计声呐平台运动误差,其空间坐标的估计误差小于仿真信号波长的1/8,从而精确补偿了CSAS在不同空间采样点上的阵元回波时间差,显著提高了目标成像质量。湖上试验结果表明,该算法能够实现对CSAS的运动误差补偿。仿真和试验结果均验证了方法的可行性和有效性。      

一种基于相对声程差的高精度逐点聚焦实现方法

针对超声成像逐点聚焦延时参数容量大、工程实现困难的问题,提出一种基于相对声程差的高精度逐点聚焦实现方法。该方法对逐点聚焦延时参数进行分解,压缩存储,并在逐点聚焦过程中实时解压生成聚焦延时参数,用于高精度的逐点聚焦中。为了验证该方法,本文以128阵元16通道的凸阵探头为例,进行了相关的数学推导和证明。在此基础上,给出了基于FPGA(Field Programmable Gate Array)的高精度逐点聚焦模块的硬件实现原理框图,并对该逐点聚焦方法的性能进行了分析讨论,验证了该方法的优越性。      

基于声学扫描振镜的超声/光声双模态成像技术

超声/光声双模态成像技术因其同时兼具超声的高分辨率结构成像和光声的高对比度功能成像优势,极大地推动了光声成像技术的临床应用推广.传统超声/光声双模态成像技术多基于超声成像所用阵列探头同时收集光声信号,系统结构紧凑且无需图像配准,操作便捷.但该类设备使用阵列探头和多通道数据采集,使得其成本较高;且成像结果易受通道一致性差异影响.本文提出了一种基于声学扫描振镜的超声/光声双模态成像技术,该技术采用单个超声换能器结合一维声学扫描振镜进行快速声束扫描,实现超声/光声双模态成像,是一种小型化、低成本的双模态快速成像技术.本文开展了系列仿体和活体成像研究,实验结果表明:系统有效成像范围为15.6 mm,超声和光声成像B扫描速度分别为1.0 s^–1和0.1 s^–1 (光声成像速度主要受制于脉冲激光器重复频率).基于本文所提技术研究,有助于进一步推动超声/光声双模态成像技术的临床转化和普及;也为基于超声信号检测的多模态成像技术提供了一种低成本、小型化和快速声信号检测的参考方案.