Full-Matrix Phase Shift Migration Method for Transcranial Ultrasonic Imaging

Transcranial ultrafast ultrasound Doppler imaging: A phantom study

Ultrafast ultrasound Doppler imaging facilitates the assessment of cerebral hemodynamics with high spatio-temporal resolution. However, the significant acoustic impedance mismatch between the skull and soft tissue results in phase aberrations, which can compromise the quality of transcranial imaging and introduce biases in velocity and direction quantification of blood flow. This paper proposed an aberration correction method that combines deep learning-based skull sound speed modelling with ray theory to realize transcranial plane-wave imaging and ultrafast Doppler imaging. The method was validated through phantom experiments using a linear array with a center frequency of 6.25 MHz, 128 elements, and a pitch of 0.3 mm. The results demonstrated an improvement in the imaging quality of intracranial targets when using the proposed method. After aberration correction, the average locating deviation decreased from 1.40 mm to 0.27 mm in the axial direction, from 0.50 mm to 0.20 mm in the lateral direction, and the average full-width-at-half-maximum (FWHM) decreased from 1.37 mm to 0.97 mm for point scatterers. For circular inclusions, the average contrast-to-noise ratio (CNR) improved from 8.1 dB to 11.0 dB, and the average eccentricity decreased from 0.36 to 0.26. Furthermore, the proposed method was applied to transcranial ultrafast Doppler flow imaging. The results showed a significant improvement in accuracy and quality for blood Doppler flow imaging. The results in the absence of the skull were considered as the reference, and the average normalized root-mean-square errors of the axial velocity component on the five selected axial profiles were reduced from 17.67% to 8.02% after the correction.

Ultrasonic Imaging of Deeper Bone Defect Using Virtual Source Synthetic Aperture with Phased Shift Migration: A Phantom Study

Ultrasound imaging for bone is a difficult task in the field of medical ultrasound. Compared with other phase array techniques, the synthetic aperture (SA) has a better lateral resolution but a limited imaging depth due to the limited ultrasonic energy emitted by the single emitter in each transmission. In contrast, the virtual source (VS) synthetic aperture allows a simultaneous multi-element emission and could provide a higher ultrasonic incident energy in each transmission. Therefore, the VS might achieve a high imaging quality at a deeper depth for bone imaging than the traditional SA. In this study, we proposed the virtual source phase shift migration (VS-PSM) method to achieve ultrasonic imaging of the deeper bone defect featured in the multilayer structure. The proposed VS-PSM method was validated using standard soft tissue phantom and printed bone phantom with artificial defects. The image quality was evaluated in terms of contrast-to-noise ratios (CNR) and amplitudes of scatters and defects at different imaging depths. The results showed that the VS-PSM method could achieve a high imaging quality of the soft tissues with a significant improvement in the scattering amplitude and without a significant sacrifice of the lateral and axial resolution. The PSM was superior to the DAS in suppressing the background noise in the images. Compared with the traditional SA-PSM, the VS-PSM method could image deeper bone defects at different ultrasonic frequencies, with an average improvement of 50% in CNR. In conclusion, this study demonstrated that the proposed VS-PSM method could image deeper bone defects and might help the diagnosis of bone disease using ultrasonic imaging.

An efficient method for transcranial ultrasound focus correction based on the coupling of boundary integrals and finite elements

Transcranial focused ultrasound is a novel technique for the noninvasive treatment of brain diseases. The success of the treatment greatly depends on achieving precise and efficient intraoperative focus. However, compensating for aberrated ultrasound waves caused by the skull through numerical simulation-based phase corrections is a challenging task due to the significant computational burden involved in solving the acoustic wave equation. In this article, we propose a promising strategy using the coupling of the boundary integral equation method (BIEM) and the finite element method (FEM) to overcome the above limitation. Specifically, we adopt the BIEM to obtain the Robin-to-Dirichlet maps on the boundaries of the skull and then couple the maps to the FEM matrices via a dual interpolation technique, resulting in a computational domain including only the skull. Three simulation experiments were conducted to evaluate the effectiveness of the proposed method, including a convergence test and two skull-induced aberration corrections in 2D and 3D ultrasound. The results show that the method's convergence is guaranteed as the element size decreases, leading to a decrease in pressure error. The computation times for simulating a 500 kHz ultrasound field on a regular desktop computer were found to be 0.47 ± 0.01 s in the 2D case and 43.72 ± 1.49 s in the 3D case, provided that lower-upper decomposition (approximately 13 s in 2D and 2.5 h in 3D) was implemented in advance. We also demonstrated that more accurate transcranial focusing can be achieved by phase correction compared to the noncorrected results (with errors of 1.02 mm vs. 6.45 mm in 2D and 0.28 mm vs. 3.07 mm in 3D). The proposed strategy is valuable for enabling online ultrasound simulations during treatment, facilitating real-time adjustments and interventions.

Ray theory-based compounded plane wave ultrasound imaging for aberration corrected transcranial imaging: Phantom experiments and simulations

Compounded plane wave imaging (CPWI) allows high-frame-rate measurement and has been one of the most promising modalities for real-time brain imaging. However, ultrasonic brain imaging using the CPWI modality is usually performed with a worn thin or removal of the skull layer. Otherwise, the skull layer is expected to distort the ultrasonic wavefronts and significantly decrease intracranial imaging quality. The motivation of this study is to investigate a CPWI method for transcranial brain imaging with the skull layer. A coordinate transformation ray-tracing (CTRT) approach was proposed to track the distorted ultrasonic wavefronts and calculate the time delays for the ultrasound plane wave passing through the skull layer. With an accurate correction for the time delays in beamforming, the CTRT-based CPWI could achieve high-quality intracranial images with the presence of skulls. The proposed CTRT-based CPWI method was verified using a simplified three-layer transcranial model. The full-wave simulation demonstrated that CTRT could accurately (i.e., relative percentage error less than0.18%) track the distorted transmitting wavefront through skull. Compared with the CPWI without aberration correction, the CTRT-based CPWI provided high-quality intracranial imaging and could accurately localize intracranial point scatterers; specifically, positioning error decreases from 0.5 mm to 0.1 mm on average in the axial direction and from 0.7 mm to 0.1 mm on average in the lateral direction. As the compounded angles increased in the CTRT-based CPWI, the contrast improved by 16.2 dB on average for the region of interest, and the array performance indicator (representing resolution) decreased by 4.0 on average for the intracranial point scatterers. The CTRT is of low computational cost compared with full wave simulation. This study suggested that the proposed CTRT-based CPWI might have the potential for real-time and non-invasive transcranial aberration-corrected imaging.

A multiscale residual U-net architecture for super-resolution ultrasonic phased array imaging from full matrix capture data

Ultrasonic phased array imaging using full-matrix capture (FMC) has raised great interest among various communities, including the nondestructive testing community, as it makes full use of the echo space to provide preferable visualization performance of inhomogeneities. The conventional way of FMC data postprocessing for imaging is through beamforming approaches, such as delay-and-sum, which suffers from limited imaging resolution and contrast-to-noise ratio. To tackle these difficulties, we propose a deep learning (DL)-based image forming approach, termed FMC-Net, to reconstruct high-quality ultrasonic images directly from FMC data. Benefitting from the remarkable capability of DL to approximate nonlinear mapping, the developed FMC-Net automatically models the underlying nonlinear wave-matter interactions; thus, it is trained end-to-end to link the FMC data to the spatial distribution of the acoustic scattering coefficient of the inspected object. Specifically, the FMC-Net is an encoder-decoder architecture composed of multiscale residual modules that make local perception at different scales for the transmitter-receiver pair combinations in the FMC data. We numerically and experimentally compared the DL imaging results to the total focusing method and wavenumber algorithm and demonstrated that the proposed FMC-Net remarkably outperforms conventional methods in terms of exceeding resolution limit and visualizing subwavelength defects. It is expected that the proposed DL approach can benefit a variety of ultrasonic array imaging applications.

Aberration correction in diagnostic ultrasound: A review of the prior field and current directions

Medical ultrasound images are reconstructed with simplifying assumptions on wave propagation, with one of the most prominent assumptions being that the imaging medium is composed of a constant sound speed. When the assumption of a constant sound speed are violated, which is true in most in vivoor clinical imaging scenarios, distortion of the transmitted and received ultrasound wavefronts appear and degrade the image quality. This distortion is known as aberration, and the techniques used to correct for the distortion are known as aberration correction techniques. Several models have been proposed to understand and correct for aberration. In this review paper, aberration and aberration correction are explored from the early models and correction techniques, including the near-field phase screen model and its associated correction techniques such as nearest-neighbor cross-correlation, to more recent models and correction techniques that incorporate spatially varying aberration and diffractive effects, such as models and techniques that rely on the estimation of the sound speed distribution in the imaging medium. In addition to historical models, future directions of ultrasound aberration correction are proposed.

Super-resolution Ultrasound Microvascular Angiography for Spinal Cord Penumbra Imaging

OBJECTIVE

After spinal cord injury (SCI) or ischemia, timely intervention in the penumbra, such as recanalization and tissue reperfusion, is essential for preservation of its function. However, limited by imaging resolution and micro-blood flow sensitivity, golden standard angiography modalities, including magnetic resonance angiography (MRA) and digital subtraction angiography (DSA), are still not applicable for spinal cord microvascular imaging. Regarding spinal cord penumbra, to the best of authors' knowledge, currently, there is no efficient in vivo imaging modality for its evaluation. With tens-of-micrometer resolution and deep penetration, advanced ultrasound localization microscopy (ULM) could potentially meet the needs of emergent diagnosis and long-term monitoring of spinal cord penumbra.

METHODS

ULM microvasculature imaging was performed on rats with all laminae removed to obtain the blood supply in major spinal cord segments (C5-L5). For adult rats with spinal cord penumbra induced by compression injury (1 s, 10 s and 15 s), ULM imaging was conducted. The corresponding angiography results are investigated in terms of microvessel saturation, morphology, and flow velocity. The Basso/Beattie/Bresnahan (BBB) locomotor rating scale and hematoxylin and eosin staining were utilized for model validation and comparison.

RESULTS

The feasibility of ULM enabling spinal cord penumbra imaging and development monitoring was demonstrated. The focal injury core and penumbra can be clearly identified using the proposed method. Significant difference of perfusion can be observed after 1 s, 10 s and 15 s compression. Quantitative results show a high correlation between in vivo ultrasonic angiography, BBB functional evaluation and ex vivo histology assessment under different compression duration.

CONCLUSION

It is demonstrated that the super-resolution ULM micro-vasculature imaging can be used to evaluate the penumbra in spinal cord at acute and early stage of chronic phase, providing an efficient modality for micro-hemodynamics monitoring of the spinal cord.

Sound Speed Estimation for Distributed Aberration Correction in Laterally Varying Media

Spatial variation in sound speed causes aberration in medical ultrasound imaging. Although our previous work has examined aberration correction in the presence of a spatially varying sound speed, practical implementations were limited to layered media due to the sound speed estimation process involved. Unfortunately, most models of layered media do not capture the lateral variations in sound speed that have the greatest aberrative effect on the image. Building upon a Fourier split-step migration technique from geophysics, this work introduces an iterative sound speed estimation and distributed aberration correction technique that can model and correct for aberrations resulting from laterally varying media. We first characterize our approach in simulations where the scattering in the media is known a-priori. Phantom and in-vivo experiments further demonstrate the capabilities of the iterative correction technique. As a result of the iterative correction scheme, point target resolution improves by up to a factor of 4 and lesion contrast improves by up to 10.0 dB in the phantom experiments presented.

Transcranial Ultrasound Imaging With Decomposition Descent Learning-Based Full Waveform Inversion

Noninvasive brain diagnosis is extremely important because of its efficiency, low cost, and painless nature in the prediction of stroke, cerebral hemorrhage, and other brain research. At present, achieving full 3-D quantitative ultrasonic imaging of the human brain is a cutting-edge challenge due to the complex structures of the human brain and the strong scattering caused by the skulls. In this article, we achieved quantitative ultrasonic imaging of inside-brain anomalies with our proposed method, the decomposition descent learning-based full waveform inversion (DDL-FWI). The proposed method adopts a linear residual decomposing technique to greatly alleviate the computation burden in fast inversion tomography (FIT) with enhanced convergence guaranteed by residual functions. Testing results in both simulation and laboratory experiments demonstrated that our method can achieve high-quality quantitative imaging of brain soft tissues and skulls even starting from homogeneous water background in 2-D, and this method is capable of reconstructing both complex brain tissues and clots in 2-D and 3-D cases using either clean or noisy signals, with a robust 3-D clot resolution as small as 18 mm and 2-D reconstruction speed in 11.20 s. Combined with advanced ultrasonic hardware, DDL-FWI can be easily trained and used for brain imaging efficiently that frees patients from harmful influences from traditional imaging techniques, e.g., ionization radiations from X-ray computed tomography (CT).

Accelerated 2-D Real-Time Refraction-Corrected Transcranial Ultrasound Imaging

In a recent study, we proposed a technique to correct aberration caused by the skull and reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. Given a sound speed map, the arrival times were calculated using a fast marching technique (FMT), which solves the Eikonal equation and, therefore, is computationally expensive for real-time imaging. In this article, we introduce a two-point ray tracing method, based on Fermat's principle, for fast calculation of the travel times in the presence of a layered aberrator in front of the ultrasound probe. The ray tracing method along with the reconstruction technique is implemented on a graphical processing unite (GPU). The point spread function (PSF) in a wire phantom image reconstructed with the FMT and the GPU implementation was studied with numerical synthetic data and experiments with a bone-mimicking plate and a sagittally cut human skull. The numerical analysis showed that the error on travel times is less than 10% of the ultrasound temporal period at 2.5 MHz. As a result, the lateral resolution was not significantly degraded compared with images reconstructed with FMT-calculated travel times. The results using the synthetic, bone-mimicking plate, and skull dataset showed that the GPU implementation causes a lateral/axial localization error of 0.10/0.20, 0.15/0.13, and 0.26/0.32 mm compared with a reference measurement (no aberrator in front of the ultrasound probe), respectively. For an imaging depth of 70 mm, the proposed GPU implementation allows reconstructing 19 frames/s with full synthetic aperture (96 transmission events) and 32 frames/s with multiangle plane wave imaging schemes (with 11 steering angles) for a pixel size of [Formula: see text]. Finally, refraction-corrected power Doppler imaging is demonstrated with a string phantom and a bone-mimicking plate placed between the probe and the moving string. The proposed approach achieves a suitable frame rate for clinical scanning while maintaining the image quality.

Index-Rotated Fast Ultrasound Imaging of Cortical Bone Based on Predicted Velocity Model

Due to the significant acoustic impedance contrast at cortical boundaries, highly inside attenuation, and the unknown sound velocity distribution, accurate ultrasound cortical bone imaging remains a challenge, especially for the traditional pulse-echo modalities using unique sound velocity. Moreover, the large amounts of data recorded by multielement probe results in a relatively time-consuming reconstruction process. To overcome these limitations, this article proposed an index-rotated fast ultrasound imaging method based on predicted velocity model (IR-FUI-VP) for cortical cross section ultrasound tomography (UST) imaging, utilizing ray-tracing synthetic aperture (RTSA). In virtue of ring probe, the sound velocity model was predicted in advance using bent-ray inversion (BRI). With the predicted velocity model, index-rotated fast ultrasound imaging (IR-FUI) was further applied to image the cortical cross sections in the sectors corresponding to the dynamic apertures (DAs) and ring center. The final result was merged by all sector images. One cortical bone phantom and two ex vivo bovine femurs were utilized to demonstrate the performance of the proposed method. Compared to the conventional synthetic aperture (SA) imaging, the method can not only accurately image the outer cortical boundary but also precisely reconstruct the inner cortical surface. The mean relative errors of the predicted sound velocity in the region of interest (ROI) were all smaller than 7%, and the mean errors of cortical thickness are all less than 0.31 mm. The reconstructed images of bovine femurs were in good agreement with the reference images scanned by micro-computed tomography ( μ CT) with respect to the morphology and thickness. The speed of IR-FUI is about 3.73 times faster than the traditional SA. It is proved that the proposed IR-FUI-VP-based UST is an effective way for fast and accurate cortical bone imaging.

The Spectrum-Beamformer for Conventional B-Mode Ultrasound Imaging System: Principle, Validation, and Robustness

Fast and efficient imaging techniques are important for real-time ultrasound imaging. The delay and sum (DAS) beamformer is the most widely-used strategy in focused ultrasound imaging (FUI) modality. However, calculating the time delays and coherently summing the amplitude response in DAS is computationally expensive and generally require a high-performance processor to realize real-time processing. In this study, an efficient spectrum beamformer, namely full-matrix capture (FMC)-stolt, is proposed in FUI system with a linear phased array. The imaging performance of FMC-stolt was validated with the point-scatter simulation and in vitro point and cyst phantoms, and then compared with that of five beamformers, that is, Multiline acquisition (MLA), retrospective transmit beamforming (RTB) in the FUI modality, as well as DAS, Garcia's frequency-wavenumber (f-k), Lu's f-k in the coherent plane wave compounding imaging (CPWCI) modality, under specific conditions. We show that the imaging performance of FMC-stolt is better than MLA-DAS in non-transmit-focal regions, and comparable with RTB-DAS at all imaging depths. FMC-stolt also shows better discontinuity alleviation than MLA and RTB. In addition, FMC-stolt has similar imaging characteristics (e.g., off-axis resolution, computational cost) as the f-k beamformers. The computational complexity and actual computational time indicate that FMC-stolt is comparable to Garcia's f-k, Lu's f-k, and faster than RTB and CPWCI-DAS if the transmitting numbers are close for FUI and CPWCI. The study demonstrates that the proposed FMC-stolt could achieve good reconstruction speed while preserving high-quality images and thus provide a choice for software beamforming for conventional B-mode ultrasound imaging, especially for hand-held devices with limited performance processors.

Refraction-Corrected Transcranial Ultrasound Imaging Through the Human Temporal Window Using a Single Probe

Transcranial ultrasound imaging (TUI) is a diagnostic modality with numerous applications, but unfortunately, it is hindered by phase aberration caused by the skull. In this article, we propose to reconstruct a transcranial B-mode image with a refraction-corrected synthetic aperture imaging (SAI) scheme. First, the compressional sound velocity of the aberrator (i.e., the skull) is estimated using the bidirectional headwave technique. The medium is described with four layers (i.e., lens, water, skull, and water), and a fast marching method calculates the travel times between individual array elements and image pixels. Finally, a delay-and-sum algorithm is used for image reconstruction with coherent compounding. The point spread function (PSF) in a wire phantom image and reconstructed with the conventional technique (using a constant sound speed throughout the medium), and the proposed method was quantified with numerical synthetic data and experiments with a bone-mimicking plate and a human skull, compared with the PSF achieved in a ground truth image of the medium without the aberrator (i.e., the bone plate or skull). A phased-array transducer (P4-1, ATL/Philips, 2.5 MHz, 96 elements, pitch = 0.295 mm) was used for the experiments. The results with the synthetic signals, the bone-mimicking plate, and the skull indicated that the proposed method reconstructs the scatterers with an average lateral/axial localization error of 0.06/0.14 mm, 0.11/0.13 mm, and 1.0/0.32 mm, respectively. With the human skull, an average contrast ratio (CR) and full-width-half-maximum (FWHM) of 37.1 dB and 1.75 mm were obtained with the proposed approach, respectively. This corresponds to an improvement of CR and FWHM by 7.1 dB and 36% compared with the conventional method, respectively. These numbers were 12.7 dB and 41% with the bone-mimicking plate.

3D Internal Visualization of Concrete Structure Using Multifaceted Data for Ultrasonic Array Pulse-Echo Tomography

This research proposes a 3D internal visualization using ultrasonic pulse-echo tomography technique to evaluate accurately the state of concrete structures for their efficient maintenance within a limited budget. Synthetic aperture focusing technique (SAFT) is used as a post-processing algorithm to manipulate the data measured by the ultrasonic pulse-echo technique. Multifaceted measurements improve the weakness of the existing ultrasonic pulse-echo tomography technique that cannot identify the area beyond a reflector as well as the area located far away from measuring surfaces. The application of apodization factor, pulse peak delay calibration and elimination of trivial response not only complements the weaknesses of the SAFT algorithm but also improves the accuracy of the SAFT algorithm. The results show that the proposed method reduces the unnecessary surface noise and improves the expressiveness of the reflector's boundaries on the resulting images. It is expected that the proposed 3D internal visualization technique will provide a useful non-destructive evaluation tool in combination with another structure evaluation method.

Fourier-Domain Ultrasonic Imaging of Cortical Bone Based on Velocity Distribution Inversion

There is a significant acoustic impedance contrast between the cortical bone and the surrounding soft tissue, resulting in difficulty for ultrasound penetration into bone tissue with high frequency. It is challenging for the conventional pulse-echo modalities to give accurate cortical bone images using uniform sound velocity model. To overcome these limitations, an ultrasound imaging method called full-matrix Fourier-domain synthetic aperture based on velocity inversion (FM-FDSA-VI) was developed to provide accurate cortical bone images. The dual linear arrays were located on the upper and lower sides of the imaging region. After full-matrix acquisition with two identical linear array probes facing each other, travel-time inversion was used to estimate the velocity distribution in advance. Then, full-matrix Fourier-domain synthetic aperture (FM-FDSA) imaging based on the estimated velocity model was applied twice to image the cortical bone, utilizing the data acquired from top and bottom linear array, respectively. Finally, to further improve the image quality, the two images were merged to give the ultimate result. The performance of the method was verified by two simulated models and two bone phantoms (i.e., regular and irregular hollow bone phantom). The mean relative errors of estimated sound velocity in the region-of-interest (ROI) are all below 12%, and the mean errors of cortical section thickness are all less than 0.3 mm. Compared to the conventional synthetic aperture (SA) imaging, the FM-FDSA-VI method is able to accurately image cortical bone with respect to the structure. Moreover, the result of irregular bone phantom was close to the image scanned by microcomputed tomography ( μ CT) in terms of macro geometry and thickness. It is demonstrated that the proposed FM-FDSA-VI method is an efficient way for cortical bone ultrasonic imaging.

Automatic Detection and Measurement of Spinous Process Curve on Clinical Ultrasound Spine Images

The ultrasound (US) imaging technique has been applied to scoliosis assessment, and the proxy Cobb angle can be acquired on the US coronal images. The spinous process angle (SPA) is a valuable parameter to indicate 3-D deformity of spine. However, the SPA cannot be measured on US images since the spinous process (SP) is merged in the soft tissue layer and impossible to be identified on the coronal view directly. A new method based on the gradient vector flow (GVF) snake model was proposed to automatically locate SP position on the US transverse images, and the density-based spatial clustering of application with noise (DBSCAN) was used to remove the outliers out of the detected location results. With marking the SP points on the US coronal image, the SP curve was interpolated and the SPA was measured. The algorithm was evaluated on 50 subjects with various severity of scoliosis, and two raters measured the SPA on both US images and radiographs manually. The mean absolute differences (MADs) of SPAs obtained from the two modalities were 3.4° ± 2.4° and 3.6° ± 2.8° for the two raters, respectively, which were less than the clinical acceptance error (5°), and the results reported a good linear correlation ( ) between the US method and radiography. It indicates that the proposed method can be a promising approach for SPA measurement using the US imaging technique.