Modeling and phantom studies of ultrasonic wall shear rate measurements using coded pulse excitation

Wall shear rate (WSR) is the derivative of blood velocity with respect to vessel radius at the endothelial cell (EC) surface. The product of WSR and blood viscosity is the wall shear stress (WSS) that has been identified as an important factor for atherosclerosis development. High echo signal-to-noise ratio (eSNR) and high spatial resolution are crucial for minimizing the errors in WSR estimates. By transmitting coded pulses with time-bandwidth product greater than one, high eSNR from weak blood scatter can be achieved without increasing instantaneous power or sacrificing spatial resolution. This paper summarizes a series of measurements in a straight tube (5-mm diameter), constant velocity flow phantom using a 10 MHz transducer (60% bandwidth, f/1.5) imaged with a 72 degrees Doppler angle, 125 MHz sampling frequency and 1 kHz pulse repetition frequency. Measurements were made using a frequency-modulated (FM) code, phase-modulated (PM) codes, and uncoded broadband and narrow band pulse transmissions. Both simulation and experimental results show that coded-pulse excitation increases accuracy and precision in WSR estimation for laminar flow over a broad range of peak velocity values when compared to standard pulsing techniques in noise-limited conditions (eSNR < 30 dB). The code sequence and its length are selected to balance range lobe suppression with eSNR and echo coherence enhancements to minimize WSR errors. In our study, the combination of an eight bit Optimal coded pulse with a Wiener compression filter yielded the highest WSR estimation performance.     

Coded excitation for diagnostic ultrasound: a system developer's perspective

Resolution and penetration are primary criteria for clinical image quality. Conventionally, high bandwidth for resolution was achieved with a short pulse, which results in a tradeoff between resolution and penetration. Coded excitation extends the bounds of this tradeoff by increasing signal-to-noise ratio (SNR) through appropriate coding on transmit and decoding on receive. Although used for about 50 years in radar, coded excitation was successfully introduced into commercial ultrasound scanners only within the last 5 years. This delay is at least partly due to practical implementation issues particular to diagnostic ultrasound, which are the focus of this paper. After reviewing the basics of biphase and chirp coding, we present simulation results to quantify tradeoffs between penetration and resolution under frequency-dependent attenuation, dynamic focusing, and nonlinear propagation. Next we compare chirp and Golay code performance with respect to image quality and system requirements, then we show clinical images that illustrate the current applications of coded excitation in B-mode, harmonic, and flow imaging.     

Coded pulse excitation for ultrasonic strain imaging

Decorrelation strain noise can be significantly reduced in low echo-signal-to-noise (eSNR) conditions using coded excitation. Large time-bandwidth-product (>30) pulses are transmitted into tissue mimicking phantoms with 2.5-mm diameter inclusions that mimic the elastic properties of breast lesions. We observed a 5-10 dB improvement in eSNR that led to a doubling of the depth of focus for strain images with no reduction of spatial resolution. In high eSNR conditions, coded excitation permits the use of higher carrier frequencies and shorter correlation windows to improve the attainable spatial resolution for strain relative to that obtained with conventional short pulses. This paper summarizes comparative studies of strain imaging in noise-limited conditions obtained by short pulses and four common aperiodic codes (chirp, Barker, suboptimal, and Golay) as a function of attenuation, eSNR and applied strain. Imaging performance is quantified using SNR for displacement (SNRd), local modulation transfer function (LMTF), and contrast-to-noise ratio for strain (CNRepsilon). We found that chirp and Golay codes are the most robust for imaging soft tissue deformation using matched filter decoding. Their superior performance is obtained by balancing the need for low-range lobes, large eSNR improvement, and short-code duration.     

巴克码激励的超声气体流量测量研究

提出了一种用“编码激励+相关测时”的超声气体流量测量方法。理论分析和数值仿真均表明,相同信噪比条件下,与传统的门限电平法及猝发脉冲激励相关法相比,巴克编码激励相关法的时延估计量方差具有更低的克拉美罗下界。采用标准表法标定实验,在200~1 200 m3/h范围内对比研究了巴克序列二进制相位调制激励相关和猝发脉冲激励相关的超声流量测量性能,前者测得流速的标准偏差和测量误差明显优于后者。     

Binary code design for high-frequency ultrasound

This paper proposes an approach to designing binary codes suitable for high-frequency applications of coded excitation in medical ultrasound. For a high-frequency ultrasound system, transmitting well-designed binary codes with a low sampling ratio (i.e., the bit rate divided by the transducer center frequency) is a practical way to improve the signal-to-noise ratio (SNR) because the challenge of implementing arbitrary-waveform generators for transmitting nonbinary codes increases with the frequency and the switching speed of square-wave pulsers are limited. One conventional approach designs codes using a base sequence that modulates wideband sequences up to the transducer passband. Because a major portion of codes is excluded as a candidate, codes designed using this approach typically need long compression filters for restoring the axial resolution, and they do not improve the SNR efficiently. In contrast, the approach proposed here searches all the codes that match the transducer passband; hence, the resultant codes exhibit better performance. The technique was tested using a bit rate of 50 MHz and a sampling ratio of 2. For a transducer with an ideal Gaussian frequency response with a center frequency of 25 MHz and a -6 dB bandwidth of 15 MHz, the SNR for the same side-lobe extent was 1 to 6 dB higher for the codes designed using the proposed approach compared with those designed using the conventional approach. When a real transducer response with a center frequency of 26.4 MHz and a one-way -6 dB bandwidth of 20.7 MHz was considered, the codes designed using the proposed approach were superior by 0.5 to 5 dB. Therefore, our approach is better than the conventional approach for designing binary codes for high-frequency ultrasound, with the results indicating that the moderate bit rate of 50 MHz will suffice when the ultrasonic center frequency is 25 MHz.     

超声编码激励在瞬时弹性成像检测中的应用

瞬时弹性(Transient Elastography,TE)成像广泛应用于肝硬化研究。然而,在临床应用中,对于肥胖病人,该方法很难实现对深度位置的瞬时剪切波进行检测。研究了将超声编码激励应用于瞬时弹性成像剪切波检测的可行性,选用7位巴克码进行编码检测研究。通过剪切波信噪比和检测穿透力两个指标对编码检测与传统短脉冲检测结果进行量化和对比。弹性仿体实验表明,编码检测可以提供比传统短脉冲检测更高的剪切波信噪比和检测深度。离体猪肝实验表明将编码激励应用于组织检测时同样可以实现高信噪比剪切波检测。这些结果表明编码检测应用于瞬时弹性成像检测是可行的,可以增加其检测深度。     

用于医学超声内窥成像的编码激励技术

为了提高医学超声内窥系统中图像的信噪比,将编码激励技术引入超声成像系统.仿真研究了编码长度和换能器相对带宽对成像信噪比（SNR）提升的影响.发现在采用长度较短的2-5位Barker编码激励时,成像信噪比的提升在编码长度为4位时达到峰值,且信噪比的提升随换能器相对带宽的增加而增加.基于超声内窥成像系统的编码激励实验表明,在峰值激励电压为25 V和换能器相对带宽为20%的情况下,采用编码激励技术能够获得1.85 dB的成像信噪比提升.     

High frequency coded imaging system with RF

Coded transmission is an approach to solve the inherent compromise between penetration and resolution required in ultrasound imaging. Our goal was to examine the applicability of the coded excitation to HF (20-35 MHz) ultrasound imaging. A novel real-time imaging system for research and evaluation of the coded transmission was developed. The digital programmable coder- digitizer module based on the field programmable gate array (FPGA) chip supports arbitrary waveform coded transmission and RF echo sampling up to 200 megasamples per second, as well as real-time streaming of digitized RF data via a high-speed USB interface to the PC. All RF and image data processing were implemented in the software. A novel balanced software architecture supports real-time processing and display at rates up to 30 frames/sec. The system was used to acquire quantitative data for sine burst and 16-bit Golay code excitation at 20 MHz fundamental frequency. SNR gain close to 14 dB was obtained. The example of the skin scan clearly shows the extended penetration and improved contrast when a 35-MHz Golay code is used. The system presented is a practical and low-cost implementation of a coded excitation technique in HF ultrasound imaging that can be used as a research tool as well as to be introduced into production.     

Use of modulated excitation signals in medical ultrasound. Part I: Basic concepts and expected benefits

This paper, the first from a series of three papers on the application of coded excitation signals in medical ultrasound, discusses the basic principles and ultrasound-related problems of pulse compression. The concepts of signal modulation and matched filtering are given, and a simple model of attenuation relates the matched filter response with the ambiguity function, known from radar. Based on this analysis and the properties of the ambiguity function, the selection of coded waveforms suitable for ultrasound imaging is discussed. It is shown that linear frequency modulation (FM) signals have the best and most robust features for ultrasound imaging. Other coded signals such as nonlinear FM and binary complementary Golay codes also have been considered and characterized in terms of signal-to-noise ratio (SNR) and sensitivity to frequency shifts. Using the simulation program Field II, it is found that in the case of linear FM signals, a SNR improvement of 12 to 18 dB can be expected for large imaging depths in attenuating media, without any depth-dependent filter compensation. In contrast, nonlinear FM modulation and binary codes are shown to give a SNR improvement of only 4 to 9 dB when processed with a matched filter. Other issues, such as depth-dependent matched filtering and use of filters other than the matched filter (inverse and Wiener filters) also are addressed.     

Directional velocity estimation using a spatio-temporal encoding technique based on frequency division for synthetic transmit aperture ultrasound

This paper investigates the possibility of flow estimation using spatio-temporal encoding of the transmissions in synthetic transmit aperture imaging (STA). The spatial encoding is based on a frequency division approach. In STA, a major disadvantage is that only a single transmitter (denoting single transducer element or a virtual source) is used in every transmission. The transmitted acoustic energy will be low compared to a conventional focused transmission in which a large part of the aperture is used. By using several transmitters simultaneously, the total transmitted energy can be increased. However, to focus the data properly, the signals originating from the different transmitters must be separated. To do so, the pass band of the transducer is divided into a number of subbands with disjoint spectral support. At every transmission, each transmitter is assigned one of the subbands. In receive, the signals are separated using a simple filtering operation. To attain high axial resolution, broadband spectra must be synthesized for each of the transmitters. By multiplexing the different waveforms on different transmitters over a number of transmissions, this can be accomplished. To further increase the transmitted energy, the waveforms are designed as linear frequency modulated signals. Therefore, the full excitation amplitude can be used during most of the transmission. The method has been evaluated for blood velocity estimation for several different velocities and incident angles. The program Field II was used. A 128-element transducer with a center frequency of 7 MHz was simulated. The 64 transmitting elements were used as the transmitting aperture and 128 elements were used as the receiving aperture. Four virtual sources were created in every transmission. By beamforming lines in the flow direction, directional data were extracted and correlated. Hereby, the velocity of the blood was estimated. The pulse repetition frequency was 16 kHz. Three different setups were investigated with flow angles of 45, 60, and 75 degrees with respect to the acoustic axis. Four different velocities were simulated for each angle at 0.10, 0.25, 0.50, and 1.00 m/s. The mean relative bias with respect to the peak flow for the three angles was less than 2%, 2%, and 4%, respectively.     

High-resolution ultrasound displacement measurement using coded excitations

Resolution of displacement measurements based on ultrasound pulse-echo techniques is limited by the center frequency of the transmitted wave, echo sampling rate, quantization errors, and electronic noises in the measurement system. We developed a new method utilizing the clutter signal in coded excitations to determine the displacement of an object or a desired region of an object with much improved resolution. The method includes transmitting a pair of Golay complementary sequences, receiving echoes from the object or a region of the object, compressing the pulse, eliminating the main lobe, and determining the object displacement between the two transmissions from the residual clutter signal around the main lobe of the compressed pulse. Results of computer simulations showed that the new method improved the resolution by several orders of magnitude and was more robust to noise than traditional pulse-echo methods. The new method was also evaluated using an experimental ultrasound system (10 MHz center frequency, 100 MHz sampling rate, and 8-bit sampling precision). A high precision in the displacement measurement was achieved with a measurement error of ?5.76 nm ±36.27 nm (mean ± standard deviation). The method has the potential to be applied in biomedical and industrial measurements of distance, displacement, and thickness.     

Coded excitation for synthetic aperture ultrasound imaging

Peak acoustic power limits the signal-to-noise ratio (SNR) of real-time ultrasound images. For most conventional scan formats, however, the average power is well below heating limits. This means the SNR can be significantly increased using coded excitation. A coded system transmits a broadband, temporally elongated excitation pulse with a finite time-bandwidth product. The received signal must be decoded to produce an imaging pulse with improved SNR resulting from the higher average power in the elongated excitation. Decoding can produce significant range side lobes, however, greatly reducing image quality. All practical coding designs, therefore, represent a trade-off between SNR gain and range side lobes. A specific coding scheme appropriate for synthetic aperture imaging is presented. A 14.5 dB SNR improvement with acceptable range side lobes is demonstrated on a forward-looking imaging system appropriate for intravascular applications.     

Echography using correlation techniques: choice of coding signal

Theoretical studies made in the early 1980's suggest that ultrasonic imaging using correlation technique can overcome some of the drawbacks of classical pulse echography. Indeed by transmitting a continuous coded signal and then compressing it into a short, high resolution pulse at the receiver the total signal to noise ratio (SNR) is improved. The target location is determined by cross correlation of the emitted and the received signal. The band compression allows, by increasing SNR, the retrieval of echo signals buried in the receiver noise. Thus in medical-type echography, where the signal attenuation at fixed depth is proportional to the frequency, the SNR improvement allows the use of higher frequency signals and leads to improved resolution. We report here the results of comparative experimental studies of simple echo B type images as obtained by the classical pulse echo and correlation techniques. Because the optimisation of the coded signal plays a crucial role in the performance of the correlation technique we will also present a comparative study of the performances of the most common codes (m-sequences and complementary series). In particular we shall emphasise the following points: the relative importance of the central lobe as compared to the side lobes of the correlation function, which is directly related to the dynamic of the imaging system, the width of the correlation peak which is directly related to the axial resolution of the system, the facility of the realisation. The merit of B-mode images obtained with the coded signals will be discussed showing that as far as signal modulation is used the best results are obtained with periodic m-sequences     

Golay-encoded excitation for dual-frequency harmonic detection of ultrasonic contrast agents

Golay-encoded excitation in combination with the third harmonic (3f?) transmit phasing is examined for both signal-to-noise ratio (SNR) and contrast-to-tissue ratio (CTR) improvements in harmonic imaging of contrast microbubbles. To produce the cancellation pair of tissue harmonic signal in 3f? transmit phasing, the phase of the bit waveform is properly designed for both the fundamental and the 3f? transmit signals to provide the Golay encoding of the received harmonic responses. Results indicate that the proposed Golay excitation can effectively suppress the tissue harmonic amplitude to increase CTR. Meanwhile, the SNR of the contrast harmonic signal also improves because of the elongated waveform of Golay excitation. Nevertheless, the generation of marked range side-lobes of the bubble region would degrade the achievable SNR improvement and the image contrast, especially when the bit of Golay excitation increases. The range side-lobes could result from the nonlinear resonance of the microbubbles that interferes with the phase modulation of the Golay encoding.     

Spatial encoding using a code division technique for fast ultrasound imaging

This paper describes a method for spatial encoding in synthetic transmit aperture ultrasound imaging. This allows several ultrasonic sources to be active simultaneously. The method is based on transmitting pseudorandom sequences to spatially encode the transmitters. The data can be decoded after only one transmission using the knowledge of the transmitted code sequences as opposed to other spatial encoding techniques, such as Hadamard or Golay encoding. This makes the method less sensitive to motion, and data can be acquired using fewer transmissions. The aim of this paper is to analyze the underlying theory and to test the feasibility in a physical system. The method has been evaluated in simulations using Field II in which the point-spread functions were simulated for different depths for a 7 MHz linear array transducer. A signal-to-noise ratio (SNR) simulation also was included in the study in which an improvement in SNR of approximately 1.5 dB was attained compared to the standard synthetic transmit aperture (STA) firing scheme. Considering the amount of energy transmitted, this value is low. A plausible explanation is given that is verified in simulation. The method also was tested in an experimental ultrasound scanner and compared to a synthetic transmit aperture ultrasound imaging scheme using a sinusoidal excitation. The performance of the proposed method was comparable to the reference with respect to axial and lateral resolution, but it displayed poorer contrast with sidelobe levels at approximately - 40 dB compared to the mainlobe.     

Harmonic chirp imaging method for ultrasound contrast agent

Coded excitation is currently used in medical ultrasound to increase signal-to-noise ratio (SNR) and penetration depth. We propose a chirp excitation method for contrast agents using the second harmonic component of the response. This method is based on a compression filter that selectively compresses and extracts the second harmonic component from the received echo signal. Simulations have shown a clear increase in response for chirp excitation over pulse excitation with the same peak amplitude. This was confirmed by two-dimensional (2-D) optical observations of bubble response with a fast framing camera. To evaluate the harmonic compression method, we applied it to simulated bubble echoes, to measured propagation harmonics, and to B-mode scans of a flow phantom and compared it to regular pulse excitation imaging. An increase of approximately 10 dB in SNR was found for chirp excitation. The compression method was found to perform well in terms of resolution. Axial resolution was in all cases within 10% of the axial resolution from pulse excitation. Range side-lobe levels were 30 dB below the main lobe for the simulated bubble echoes and measured propagation harmonics. However, side-lobes were visible in the B-mode contrast images.     

Basic considerations in the use of coded excitation for color flow imaging applications

Coded excitation is now a well-established technique in medical ultrasound for B-mode imaging applications. It enables a gain in depth of penetration, without sacrificing the spatial resolution and maintaining an acceptable peak intensity for patient safety. The rationale of this technique for velocity estimation applications still has to be formulated in more precise terms. In particular, differences in the situation that arise in color flow imaging (CFI) applications from typical B-mode imaging conditions, such as signal-to-noise ratio conditions, pulsing strategy, and safety requirements, need to be specifically addressed to assess more quantitatively the potential of this technique. This paper discusses the potential improvement in sensitivity, resolution, and statistical performance provided by coded excitation for CFI applications from theoretical considerations and simulations.     

Design and performance of Huffman sequences in medical ultrasound coded excitation

This paper deals with coded-excitation techniques for ultrasound medical echography. Specifically, linear Huffman coding is proposed as an alternative approach to other widely established techniques, such as complementary Golay coding and linear frequency modulation. The code design is guided by an optimization procedure that boosts the signal-to-noise ratio gain (GSNR) and, interestingly, also makes the code robust in pulsed-Doppler applications. The paper capitalizes on a thorough analytical model that can be used to design any linear coded-excitation system. This model highlights that the performance in frequency-dependent attenuating media mostly depends on the pulse-shaping waveform when the codes are characterized by almost ideal (i.e., Kronecker delta) autocorrelation. In this framework, different pulse shapers and different code lengths are considered to identify coded signals that optimize the contrast resolution at the output of the receiver pulse compression. Computer simulations confirm that the proposed Huffman codes are particularly effective, and that there are scenarios in which they may be preferable to the other established approaches, both in attenuating and non-attenuating media. Specifically, for a single scatterer at 150 mm in a 0.7-dB/(MHz·cm) attenuating medium, the proposed Huffman design achieves a main-to-side lobe ratio (MSR) equal to 65 dB, whereas tapered linear frequency modulation and classical complementary Golay codes achieve 35 and 45 dB, respectively.     

Estimation of ultrasound attenuation from broadband echo-signals using bandpass filtering

A problem with video signal analysis for estimating frequency-dependent ultrasonic attenuation is that relative echogenicity versus depth curves are distorted when broadband pulses are used. In this correspondence, we present results that demonstrate improved accuracy of attenuation estimates computed from B-mode or envelope data derived after narrowband (1 MHz BW) filtering at different frequencies around the center frequency of the broadband echo signal. Based on the premise that transducer center frequencies are selected in part on penetration or imaging depth requirements, simulation and experimental results for a typical ultrasound imaging system show that narrowband video signal analysis for frequencies lower than or at the center frequency of the broadband pulse provide unbiased attenuation estimation over this depth. Filtered signals above the center frequency of the pulse yield accurate results only at shallow depths.     

Ultrasound Transducer and System for Real-Time Simultaneous Therapy and Diagnosis for Noninvasive Surgery of Prostate Tissue

For noninvasive treatment of prostate tissue using high-intensity focused ultrasound this paper proposes a design of an integrated multifunctional confocal phased array (IMCPA) and a strategy to perform both imaging and therapy simultaneously with this array. IMCPA is composed of triplerow phased arrays: a 6-MHz array in the center row for imaging and two 4-MHz arrays in the outer rows for therapy. Different types of piezoelectric materials and stack configurations may be employed to maximize their respective functionalities, i.e., therapy and imaging. Fabrication complexity of IMCPA may be reduced by assembling already constructed arrays. In IMCPA, reflected therapeutic signals may corrupt the quality of imaging signals received by the center-row array. This problem can be overcome by implementing a coded excitation approach and/or a notch filter when B-mode images are formed during therapy. The 13-bit Barker code, which is a binary code with unique autocorrelation properties, is preferred for implementing coded excitation, although other codes may also be used. From both Field II simulation and experimental results, we verified whether these remedial approaches would make it feasible to simultaneously carry out imaging and therapy by IMCPA. The results showed that the 13-bit Barker code with 3 cycles per bit provided acceptable performances. The measured -6 dB and -20 dB range mainlobe widths were 0.52 mm and 0.91 mm, respectively, and a range sidelobe level was measured to be -48 dB regardless of whether a notch filter was used. The 13-bit Barker code with 2 cycles per bit yielded -6 dB and -20 dB range mainlobe widths of 0.39 mm and 0.67 mm. Its range sidelobe level was found to be -40 dB after notch filtering. These results indicate the feasibility of the proposed transducer design and system for real-time imaging during therapy.