Array Of Transducer Elements Research Articles

A frequency-domain model-based reconstruction method for transcranial photoacoustic imaging: A 2D numerical investigation

Phase aberration caused by the skull is a major barrier to achieving high quality photoacoustic images of human and non-human primates’ brains. To address this issue, time-reversal methods have been used but they are computationally demanding and slow due to relying on solving the full-wave equation. The proposed approach is based on model-based image reconstruction in the frequency-domain to achieve near real-time image reconstruction. The relationship between an imaging region and transducer array elements can be mathematically described as a model matrix and the image reconstruction can be performed by pseudo-inverse of the model matrix. The model matrix is numerically calculated due to the lack of analytical solutions for transcranial ultrasound. However, this calculation only needs to be performed once for a given experimental setup and the same acoustic medium, and is an offline process not affecting the actual image reconstruction time. This non-iterative mode-based method demonstrates a substantial improvement in image reconstruction time, being approximately 18 times faster than the time-reversal method, all while maintaining comparable image quality.

All-digital transmit beamformer for portable high-frequency ultrasound imaging systems.

To meet the requirements of high-frequency ultrasound imaging systems, a transmit-beamforming integrated circuit with higher delay resolution than conventional transmit-beamforming circuits, which are typically implemented using field-programmable gate array chips, is presented. It also requires smaller volumes, allowing for portable applications. Its proposed design includes two all-digital delay-locked loops providing a specified digital control code for a counter-based beamforming delay chain (CBDC) to generate stable and suitable delays for exciting the array transducer elements without variations in process, voltage, and temperature. Moreover, to maintain the duty cycle of long propagation signals, this novel CBDC requires only a few delay cells, significantly reducing hardware costs and power consumption. Simulations were conducted, revealing a maximum time delay of 451.9ns with a time resolution of 652ps and a maximum lateral resolution error of 0.04mm at 6.8mm.

Compensation for Electrical Fatigue in Piezoelectric Transducers

Malfunction of ultrasonic transducers is one of the most common factors affecting the quality of medical ultrasound images. In this work, the effect of electrical fatigue in piezoelectric transducers on the image quality is investigated and subsequently compensated. To this end, through a combination of simulation and controlled experimental and clinical evaluations, the effect of fatigue and its compensation is realized. Initially, the electrical fatigue is simulated by reducing the piezoelectric permittivity and piezoelectric constant parameters based on finite element method (FEM). Next, the transducer’s electromechanical impulse response is obtained, and the quality of reconstructed images from the fatigued transducer array elements is evaluated by applying the FEM simulation results to the Field II simulation program. Also, experimental measurements are conducted using a tissue-mimicking phantom and in vivo imaging. To evaluate the image quality, various metrics such as signal to noise ratio (SNR), side lobe level (SLL), full width at half maximum (FWHM) and contrast to noise ratio (CNR) are measured for the simulated and tissue-mimicking phantoms. To compensate for negative effects of electrical fatigue, a method based on the nonblind deconvolution of the ultrasound radio frequency (RF) signals is proposed. This method utilizes the corresponded electromechanical impulse responses to the fatigue level as the known undesired impulse responses. The results show that the image quality metrics which were significantly degraded due to electrical fatigue, are appropriately improved as a result of the proposed restoration method.

Fabrication, Acoustic Characterization and Phase Reference-Based Calibration Method for a Single-Sided Multi-Channel Ultrasonic Actuator.

The ultrasonic actuator can be used in medical applications because it is label-free, biocompatible, and has a demonstrated history of safe operation. Therefore, there is an increasing interest in using an ultrasonic actuator in the non-contact manipulation of micromachines in various materials and sizes for therapeutic applications. This research aims to design, fabricate, and characterize a single-sided transducer array with 56 channels operating at 500 kHz, which provide benefits in the penetration of tissue. The fabricated transducer is calibrated using a phase reference calibration method to reduce position misalignment and phase discrepancies caused by acoustic interaction. The acoustic fields generated by the transducer array are measured in a 300 mm × 300 mm × 300 mm container filled with de-ionized water. A hydrophone is used to measure the far field in each transducer array element, and the 3D holographic pattern is analyzed based on the scanned acoustic pressure fields. Next, the phase reference calibration is applied to each transducer in the ultrasonic actuator. As a result, the homogeneity of the acoustic pressure fields surrounding the foci area is improved, and the maximum pressure is also increased in the twin trap. Finally, we demonstrate the capability to trap and manipulate micromachines with acoustic power by generating a twin trap using both optical camera and ultrasound imaging systems in a water medium. This work not only provides a comprehensive study on acoustic actuators but also inspires the next generation to use acoustics in medical applications.

Algorithm to Correct Measurement Offsets Introduced by Inactive Elements of Transducer Arrays in Ultrasonic Flow Metering

Ultrasonic flow meters (UFMs) based on transducer arrays offer several advantages. With electronic beam steering, it is possible to tune the steering angle of the beam for optimal signal-tonoise ratio (SNR) upon reception. Moreover, multiple beams can be generated to propagate through different travel paths, covering a wider section of the flow profile. Furthermore, in a clamp-on configuration, UFMs based on transducer arrays can perform self-calibration. In this manner, userinput is minimized and measurement repeatability is increased. In practice, transducer array elements may break down. This could happen due to aging, exposure to rough environments, and/or rough mechanical contact. As a consequence of inactive array elements, the measured transit time difference contains two offsets. One offset originates from non-uniform spatial sampling of the generated wavefield. Another offset originates from the ill-defined beam propagating through a travel path different from the intended one. In this paper, an algorithm is proposed that corrects for both of these offsets. The algorithm also performs a filtering operation in the frequency-wavenumber domain of all spurious (i.e., flow-insensitive) wave modes. The advantage of implementing the proposed algorithm is demonstrated on simulations and measurements, showing improved accuracy and precision of the transit time differences compared to the values obtained when the algorithm is not applied. The proposed algorithm can be implemented in both in-line and clamp-on configuration of UFMs based on transducer arrays.

FPCB as an Acoustic Matching Layer for 1D Linear Ultrasound Transducer Arrays.

An acoustic matching layer is an essential component of an ultrasound transducer to achieve maximum ultrasound transmission efficiency. Here, we develop a flexible printed circuit board (FPCB) with a composite structure consisting of multiple polyimide and copper layers and demonstrate it as a novel acoustic matching layer. With a flexible substrate and robust ACF bonding, the FPCB not only serves as an acoustic matching layer between piezoelectric elements and the surrounding medium but also as a ground for the electrical connection between the transducer array elements and the folded substrate. A 1D linear ultrasound transducer array with the FPCB matching layer exhibits larger output pressure, wider -3dB bandwidth, and higher ultrasound beam intensity compared to that of an ultrasound transducer array with the alumina/epoxy matching layer, which is one of the most commonly applied composite matching layers. The enhanced transmission performance verifies that the proposed FPCB is an excellent matching layer for 1D linear ultrasound transducer arrays.

A Generic Approach for Computation of Near-Field and Far-Field Pattern of Beamforming Acoustic Transducer Array

This paper addresses a mathematical formulation to determine the radiation performance of an acoustic beamforming array, which is valid for both near-field and far-field regions. The developed beamforming analysis is generic for application to any array configuration, with the array elements having a radiating aperture of canonical contours. The beamforming analysis presented in this paper for the acoustic transducer array is obtained by the process of vector summation of radiated fields from the elements of the acoustic transducer array with an array factor and beam steering phase factor. The initial validity of the beamforming analysis is examined by considering the far-field region as a particular case of the near-field region with the relevant distance criterion. The numerous simulation studies presented in this paper show that the proposed method produces excellent results with the distinct structure of sidelobes at a distance RNF ≥ 0.5D2/λ from the acoustic aperture of (7(1) linear and (7(7) planar array of acoustic transducers. If the distance between the transducer array and the point of observation is very low (RNF < 0.5D2/λ), simulation results indicate a less directive main beam. The two-dimensional near field beam steering method proposed in this paper finds the significant utility in many practical applications such as Underwater Acoustic Imaging (UAI) systems and Underwater Acoustic Communication (UAC).

Accurate coupled vibration analysis of a piezoelectric array element by the superposition method

In this paper, the superposition method is extended to obtain analytical solutions for the coupled vibration of a piezoelectric slender bar in a configuration corresponding to a typical ultrasonic linear array transducer element. The problem can be described mathematically by three partial differential equations with electrical and mechanical boundary conditions. To solve this, the vibrations in lateral and thickness directions are referred to as two building blocks. In each building block, the expressions of displacements and electric potential are assumed first based on their symmetry properties and then the induced dynamic responses, such as in-plane stress and electric displacements, are calculated. Finally, the vibration responses of the two building blocks are superimposed to satisfy the boundary conditions using Fourier series expansions. Electrical impedance and mode shapes, represented by the spatial distribution of displacements and electric potential, are calculated analytically and compared with the results of the finite element method. An excellent agreement is observed. The method can be applied to design and optimize piezoelectric array transducers for various applications.

Superposition model of mode shapes composed of travelling torsional guided waves excited by multiple circular transducer arrays in pipes

In pipe inspection using ultrasonic guided wave technique, the current commercial transmitters are designed for the unidirectional guided wave excitation using multiple circular piezoelectric transducers arrays in the axial direction. However, the source with many individual transducer elements in arrays has difficulty in achieving an axisymmetric loading perfectly for defect detection. Therefore, a quasi-axisymmetric wave is formed due to many undesired wave modes are launched instead of a pure axisymmetric wave at a given excitation frequency. In this paper, a realistic superposition model of axial multiple transducer arrays is proposed. The model has many potential applications; one example is investigating the source influence on the generated quasi-axisymmetric wave effect. The analytical model is employed to achieve the predictions for investigating a transmitter’s influence for the unidirectional enhancement of torsional T(0,1) guided wave mode excitation in a pipe inspection system composing of three piezoelectric transducer ring arrays. The excitation function with variable power levels among transducers in arrays is also introduced. The predictive results using the analytical model for the distribution of circumferential displacement amplitudes over time are verified using the finite element method and a CLV-3D laser vibrometry measurement on a 219.1-mm-outer-diameter steel pipe without defect. A comparison between calculated and test results has been analysed quantitatively. The respective results are in good agreement. Thus, predictions for the superposed wavefield can be used to analyse the realistic characterisation of the excitation function in axial multiple transducer arrays. Additionally, a sensitivity analysis for part-circumferential crack detection using the quasi-axisymmetric torsional modes generated is also evaluated using finite element modelling.

32-Channel transmit beamformer with high timing resolution for high-frequency ultrasound imaging systems.

This paper presents a 32-channel high timing resolution transmit-beamforming circuit for use in high-frequency ultrasound imaging systems. Conventional transmit-beamforming circuits are typically implemented using field-programmable gate array (FPGA) chips. However, it is difficult for FPGAs to provide high timing resolution to meet the beamforming requirements of high-frequency ultrasound imaging systems. The proposed transmit-beamforming design can generate stable and suitable delays to excite 32-channel array transducer elements without variations in the process, voltage, and temperature. In addition, the proposed low-complexity architecture can maintain the duty cycle of long prorogation signals with low hardware cost to meet the timing requirements of a large channel number array transducer. The proposed designed transmit beamformer uses 0.18-µm CMOS technology for a 30-MHz high-frequency linear array, and the simulation results show that the proposed transmit-beamforming application-specific integrated circuit can achieve a maximum time delay of 619.5 ns with a time resolution of 617 ps.

High-Resolution All-Digital Transmit Beamformer for High-Frequency and Wearable Ultrasound Imaging Systems

In this article, we present a high delay resolution and all-digital transmit-beamforming circuit for use in high-frequency and wearable ultrasound imaging systems. The proposed transmit beamformer comprises an all-digital delay-locked loop (ADDLL), a cycle-based beamforming delay line (CBDL), and a pulse controller. ADDLL can provide specific digital control code for the CBDL to generate stable and suitable delays to excite array transducer elements without variations in process, voltage, and temperature. CBDL not only provides a finer delay resolution for high-frequency ultrasound imaging systems but also requires a few delay cells to maintain the duty cycle of long prorogation signals, thereby significantly reducing hardware costs. The pulse controller can adjust the number and width of output pulses to meet the requirements of different ultrasonic transducers and Doppler ultrasound color imaging. The proposed transmit beamformer uses a 0.18- $\mu \text{m}$ CMOS technology for a 30-MHz high-frequency linear array with 256 elements, and the focusing depth is from 3.5 to 8.1 mm. The measurement results indicate that the proposed transmit-beamforming design can achieve a maximum time delay of 562.6 ns with an average time resolution of 535 ps. The proposed design offers a smaller volume and finer delay resolution, and it is highly suitable for wearable ultrasound imaging systems.

Ring Artifact Correction for Phase-Insensitive Ultrasound Computed Tomography.

An algorithm was developed for the correction of ring artifacts in phase-insensitive ultrasound computed tomography attenuation images. Differences in the measurement sensitivity between the ultrasound transducer array elements cause discontinuities in the sinogram which manifest as rings and arcs in the reconstructed image. The magnitudes of the discontinuities are potentially time-varying and dependent on the attenuation being measured. The algorithm dynamically determines the measurement sensitivity of each transducer in the array during the scan by comparison with both the elements to its left and the elements to its right. Elements at either end of the array are corrected, assuming a zero-attenuation path. The two estimates of sensitivity are combined using a weighted mean similar to a Kalman filter. The algorithm was tested on simulated and experimentally acquired data. It was demonstrated to reduce the root-mean-square error (RMSE) of simulated images against ground-truth images by up to a factor of 50 compared with uncorrected images and to visibly reduce artifacts on images reconstructed from the experimentally acquired data.

Numerical Study of Square Border Ultrasonic Transducer Element Arrays for Breast Imaging in Ultrasound Computed Tomography with Waveform Inversion Method

Ultrasound computed tomography (USCT) is a promising approach for early breast cancer screening. Existing studies that estimate breast images using Waveform inversion (WI) methods utilize a Circular ultrasonic transducer element array (CUTA) to collect the measurement data. However, very accurate transducer element positioning and directivity are required for signal calibration in the hardware system of these studies, thereby causing difficulties in hardware manufacturing. The purpose of this study is to estimate high-resolution USCT images using a WI method with a Square border ultrasonic transducer element array (SUTA), which reduces the difficulty of hardware manufacturing and creates accurate correspondence between the continuous and discrete forms of the transducer element positions. Therefore, in this study, an SUTA is adopted to collect the measurement data of breast imaging using the WI method. A penalized least-squares optimization problem is constructed to obtain the numerical solution of a sound speed distribution. Computer simulations are conducted to compare images reconstructed from the measurement data that are collected using an SUTA and a CUTA. The performance is evaluated using a numerical breast phantom. Results suggest that the biases of the images reconstructed are less than 1% with the evenly distributed SUTA and 1% with the CUTA, under a noise condition.

Tailoring reflected and diffracted wave fields from tessellated acoustic arrays by origami folding

Recently the physical reconfigurations of origami-inspired acoustic arrays have been leveraged for adaptive guidance of acoustic wave energy. The spatial translations and rotations of transducer elements in tessellated arrays may tailor how reflected and diffracted waves contribute to the total acoustic pressure at the field point. To date such wave phenomena have not been illuminated, which inhibits understanding on how to best leverage the concept for highly reconfigurable arrays. Here, the reflection and diffraction phenomena induced from a Miura-ori based acoustic array are examined, wherein such wave fields superpose with directly radiated waves to yield the total acoustic response at receiving locations. As tessellated arrays take on highly folded configurations, reflected waves are shown to greatly contribute to the acoustic field at points normal to the plane of the unfolded array, whereas diffracted waves primarily govern the acoustic pressure for locations at nearly grazing angles of incidence to the unfolded array plane. Supported by experimental validation, the modeling methods established here can be extended to arbitrary wave-radiating architectures and may inspire investigations for other physical domains where reconfigurable array concepts are of great interest.

Addressing spatiotemporal distortion of high-speed tissue motion in B-mode ultrasound

The sequential excitation of transducer elements in linear ultrasound arrays can produce distorted images in B-mode since the A-lines in an image frame are not captured at the same time. If the imaged tissue is moving laterally during image acquisition, substantial distortion can occur and lead to errors in tissue displacement and velocity measurements, particularly for high-speed tissue motion such as in injury biomechanics tests. A method for correcting for this distortion in tracked time-series displacement data was developed for systems in which each A-line of the image is completed at a time proportional to its lateral position across the transducer linear array. The proposed method involves estimating more precisely the time at which each point in the image was captured, rather than relying on the time elapsed between image frames to estimate timing. The developed method was evaluated using results from motion tests of human dura imaged using a Vevo 2100 Imaging System at a frame rate of 865 frames per second. At tissue speeds ranging from 0.2 m s−1 to 2.3 m s−1, tracking of tissue displacement without correction resulted in mean displacement measurement error of 17%. Correction using proposed methods reduced error to a mean of 2%. This correction concept has the potential to improve the accuracy of velocity or displacement measurements in any B-mode imaging application where tissue velocity is high relative to ultrasound sweep velocity including in injury biomechanics and cardiac biomechanics research.

Three-Dimensional Ultrasonic Needle Tip Tracking with a Fiber-Optic Ultrasound Receiver.

Ultrasound is frequently used for guiding minimally invasive procedures, but visualizing medical devices is often challenging with this imaging modality. When visualization is lost, the medical device can cause trauma to critical tissue structures. Here, a method to track the needle tip during ultrasound image-guided procedures is presented. This method involves the use of a fiber-optic ultrasound receiver that is affixed within the cannula of a medical needle to communicate ultrasonically with the external ultrasound probe. This custom probe comprises a central transducer element array and side element arrays. In addition to conventional two-dimensional (2D) B-mode ultrasound imaging provided by the central array, three-dimensional (3D) needle tip tracking is provided by the side arrays. For B-mode ultrasound imaging, a standard transmit-receive sequence with electronic beamforming is performed. For ultrasonic tracking, Golay-coded ultrasound transmissions from the 4 side arrays are received by the hydrophone sensor, and subsequently the received signals are decoded to identify the needle tip's spatial location with respect to the ultrasound imaging probe. As a preliminary validation of this method, insertions of the needle/hydrophone pair were performed in clinically realistic contexts. This novel ultrasound imaging/tracking method is compatible with current clinical workflow, and it provides reliable device tracking during in-plane and out-of-plane needle insertions.

Three-Dimensional Ultrasonic Needle Tip Tracking with a Fiber-Optic Ultrasound Receiver

Ultrasound is frequently used for guiding minimally invasive procedures, but visualizing medical devices is often challenging with this imaging modality. When visualization is lost, the medical device can cause trauma to critical tissue structures. Here, a method to track the needle tip during ultrasound image-guided procedures is presented. This method involves the use of a fiber-optic ultrasound receiver that is affixed within the cannula of a medical needle to communicate ultrasonically with the external ultrasound probe. This custom probe comprises a central transducer element array and side element arrays. In addition to conventional two-dimensional (2D) B-mode ultrasound imaging provided by the central array, three-dimensional (3D) needle tip tracking is provided by the side arrays. For B-mode ultrasound imaging, a standard transmit-receive sequence with electronic beamforming is performed. For ultrasonic tracking, Golay-coded ultrasound transmissions from the 4 side arrays are received by the hydrophone sensor, and subsequently the received signals are decoded to identify the needle tip's spatial location with respect to the ultrasound imaging probe. As a preliminary validation of this method, insertions of the needle/hydrophone pair were performed in clinically realistic contexts. This novel ultrasound imaging/tracking method is compatible with current clinical workflow, and it provides reliable device tracking during in-plane and out-of-plane needle insertions.

Portable ultrasound imaging system with super-resolution capabilities

This paper discusses an ultrasound technique where the echo signals from the array of transducer elements are compressed to as few as two RF channels while still in analog domain, with a much simplified front-end electronics. The method can achieve resolutions well beyond the diffraction limit, which is set by the excitation signal wavelength and numerical aperture of the imaging system. The fundamental principle that underlies this model based imaging technique is the preservation of the spatial frequency information content of the recorded echo signals with the help of pseudo-random apodization function followed by summation. A Verasonics V1 ultrasonic scanner is used to conduct experiments using an anechoic cyst made from gel phantom, immersed in degassed water. The estimated images were compared to those obtained using traditional B-mode delay-and-sum imaging available with the Verasonics V1 ultrasound machine. The estimated images using the proposed imaging technique showed a contrast ratio of 0.96 and Full-Width-Half-Maximum (FWHM) of about half the wavelength at a depth of 9.1 cm and at 1.875 MHz center frequency while the traditional delay and sum images had a contrast ratio of 0.62 and FWHM of about 5.5 wavelengths.

Measurements of the Relationship Between CT Hounsfield Units and Acoustic Velocity and How It Changes With Photon Energy and Reconstruction Method.

Transcranial magnetic resonance-guided focused ultrasound continues to gain traction as a noninvasive treatment option for a variety of pathologies. Focusing ultrasound through the skull can be accomplished by adding a phase correction to each element of a hemispherical transducer array. The phase corrections are determined with acoustic simulations that rely on speed of sound estimates derived from CT scans. While several studies have investigated the relationship between acoustic velocity and CT Hounsfield units (HUs), these studies have largely ignored the impact of X-ray energy, reconstruction method, and reconstruction kernel on the measured HU, and therefore the estimated velocity, and none have measured the relationship directly. In this paper, 91 ex vivo human skull fragments from two skulls are imaged by 80 CT scans with a variety of energies and reconstruction methods. The average HU from each fragment is found for each scan and correlated with the speed of sound measured using a through transmission technique in that fragment. As measured by the -squared value, the results show that CT is able to account for 23%-53% of the variation in velocity in the human skull. Both the X-ray energy and the reconstruction technique significantly alter the -squared value and the linear relationship between HU and speed of sound in bone. Accounting for these variations will lead to more accurate phase corrections and more efficient transmission of acoustic energy through the skull.

Multi-ray medical ultrasound simulation without explicit speckle modelling.

To develop a medical ultrasound (US) simulation method using T1-weighted magnetic resonance images (MRI) as the input that offers a compromise between low-cost ray-based and high-cost realistic wave-based simulations. The proposed method uses a novel multi-ray image formation approach with a virtual phased array transducer probe. A domain model is built from input MR images. Multiple virtual acoustic rays are emerged from each element of the linear transducer array. Reflected and transmitted acoustic energy at discrete points along each ray is computed independently. Simulated US images are computed by fusion of the reflected energy along multiple rays from multiple transducers, while phase delays due to differences in distances to transducers are taken into account. A preliminary implementation using GPUs is presented. Preliminary results show that the multi-ray approach is capable of generating view point-dependent realistic US images with an inherent Rician distributed speckle pattern automatically. The proposed simulator can reproduce the shadowing artefacts and demonstrates frequency dependence apt for practical training purposes. We also have presented preliminary results towards the utilization of the method for real-time simulations. The proposed method offers a low-cost near-real-time wave-like simulation of realistic US images from input MR data. It can further be improved to cover the pathological findings using an improved domain model, without any algorithmic updates. Such a domain model would require lesion segmentation or manual embedding of virtual pathologies for training purposes.