Full waveform inversion using frequency shift envelope-based global correlation norm for ultrasound computed tomography.

Many studies have been carried out on ultrasound computed tomography (USCT) for its ability to offer quantitative measurements of tissue sound speed. Full waveform inversion (FWI) is a technique for reconstructing high-resolution sound speed images by iteratively minimizing the difference between the observed ultrasound data and the synthetic data based on the waveform equation. However, FWI suffers from cycle-skipping, which usually causes FWI convergence at a local minimum. Cycle-skipping occurs when the phase difference between the observed data and the synthetic data exceeds half a cycle. The simplest way to avoid cycle-skipping is to use low-frequency information for reconstruction. Nevertheless, in imaging systems, the response bandwidth of the probe is limited, and reliable low-frequency information often exceeds the response band. Therefore, it is a challenge to perform FWI imaging and avoid cycle-skipping problems without low-frequency information. In this paper, we propose a frequency shift envelope-based global correlation norm (FSEGCN), where an artificial source wavelet with a lower frequency is adopted to calculate synthetic data. FSEGCN compared with FWI, envelope inversion (EI), global correlation norm (GCN), envelope-based global correlation norm (EGCN) through concentric circle phantom without low-frequency information. The experimental results demonstrated the capability of the proposed method to recover the sound speed close to the exact model in the absence of low-frequency information, whereas FWI, EI, GCN, and EGCN cannot. Experiments on phantoms of the human head and calf show that artificial source wavelets can reduce image artifacts and enhance reconstruction robustness, when original low-frequency information is absent.

A coupling, stabilizing, and shaping strategy for breast ultrasound computed tomography (USCT) with a ring array transducer.

Breast ultrasound computed tomography (USCT) has been gradually promoted to clinical application after years of rapid development. Compared with the traditional handheld ultrasound scanning method, the scanning plane of USCT is fixed at the coronal plane, and the scanning path is designed in advance; the acoustic window is not in direct contact with the breast, a lot of coupling medium (usually degassed water is used to fill the gaps between the probe and breast. The clinical application of breast USTC faces challenges: (1) the processes of water degassing, heating, filling, draining, and cleaning prolong the entire scan cycle and reduce patient throughput. (2) The breast is not stabilized and slight movements of the breast may cause motion artifacts in the USCT images. (3) The non-normal incidence of ultrasound into the breast causes reflected and transmitted signals received with a low signal-to-noise ratio (SNR) or even unable to be detected. This article proposes a coupling, stabilizing, and shaping strategy for the clinical application of USCT with a ring array transducer. The solid gel coupling agent (SGCA) is applied for coupling, and a set of SGCA moldings is designed to stabilize and shape the breast during scanning, the breast shape and size which vary from person to person are simplified into several models. The preparation time is reduced to less than 1 min by replacing disposable moldings. The results show that the breast after shaping is close to round in the coronal plane, and slopes of the breast skin are limited in the sagittal and transverse planes, the breast subcutaneous tissue (fat and glands) has a better contrast-to-noise ratio (CNR) and can be better distinguished in the reflection images than that of the breast without shaping. The mean value of the raw beamformed data which represents the reflection signal amplitude of breast subcutaneous tissue after shaping shows 1.5 times that of the breast without shaping, the signal-to-noise ratio (SNR) of the raw transmission signal data after breast shaping is overall higher than that of the breast without shaping. The application of SGCA moldings for breast coupling, stabilizing, and shaping also benefits establishing a standardized scanning process, the standardized diagnosis of the breast lesion, and the localization of breast lesions.

Effect of spatial-domain pulse width on the resolution of scattering images in ultrasound computed tomography.

Full-aperture tomography (FAT) is the major image reconstruction method for a circular ring array (CRA)-based ultrasound computed tomography (USCT) system. The FAT technique requires transferring the reconstruction process from the temporal domain to the spatial domain, during which the imaging resolution of the USCT is degraded by the spatial-domain pulse width (SDPW) of backprojection areas. To tackle this challenge, this study investigates the characteristics of the SDPW and how it degrades the image resolution. We show that the SDPW depends on the frequency of the ultrasound and the position of the transmitting elements, receiving elements and the imaging point. To quantify the deterioration of image resolution associated with the position of the transmitting and receiving elements, a SDPW broadening factor (SDPWBF) is introduced. The results of numerical simulation show a smaller SDPWBFprovides a better reflection image resolution, and the distribution of SDPWBFshows that a shorter distance between the receiving element and the transmitting element yields a smaller SDPWBF. The SDPWBFis therefore able to be an indicator of selecting the signals acquired from the transmitting and receiving elements to perform optimal image resolution. Single-scatterer phantom andinvivoexperiments demonstrate how the SDPWBFaffects the USCT image spatial resolution and signal-to-noise ratio (SNR), and the results agree well with the theoretical predictions.

Design and Implementation of a Modular and Scalable Research Platform for Ultrasound Computed Tomography.

Increasing attention has been attracted to the research of ultrasound computed tomography (USCT). This article reports the design considerations and implementation details of a novel USCT research system named UltraLucid, which aims to provide a user-friendly platform for researchers to develop new algorithms and conduct clinical trials. The modular design strategy is adopted to make the system highly scalable. A prototype has been assembled in our laboratory, which is equipped with a 2048-element ring transducer, 1024 transmit (TX) channels, 1024 receive (RX) channels, two servers, and a control unit. The prototype can acquire raw data from 1024 channels simultaneously using a modular data acquisition and a transfer system, consisting of 16 excitation and data acquisition (EDAQ) boards. Each EDAQ board has 64 independent TX and RX channels and 4-Gb Ethernet interfaces for raw data transmission. The raw data can be transferred to two servers at a theoretical rate of 64 Gb/s. Both servers are equipped with a 10.9-TB solid-state drive (SSD) array that can store raw data for offline processing. Alternatively, after processing by onboard field-programmable gate arrays (FPGAs), the raw data can be processed online using multicore central processing units (CPUs) and graphics processing units (GPUs) in each server. Through control software running on the host computer, the researchers can configure parameters for transmission, reception, and data acquisition. Novel TX-RX scheme and coded imaging can be implemented. The modular hardware structure and the software-based processing strategy make the system highly scalable and flexible. The system performance is evaluated with phantoms and in vivo experiments.

Zone-Shrinking Fresnel Zone Travel-Time Tomography for Sound Speed Reconstruction in Breast USCT.

Many studies have been carried out on ultrasound computed tomography (USCT) for its potential application in breast imaging. The sound speed (SS) image modality in USCT can help doctors diagnose the breast cancer, as the tumor usually has a higher sound speed than normal tissues. Travel time is commonly used to reconstruct SS image. Raypath travel-time tomography (RTT) assumes that the sound wave travels through a raypath. RTT is computationally efficient but with low contrast to noise ratio (CNR). Fresnel zone travel-time tomography (FZTT) is based on the assumption that the sound wave travels through an area called the Fresnel zone. FZTT can provide SS image with high CNR but low accuracy due to the wide Fresnel zone. Here, we propose a zone-shrinking Fresnel zone travel-time tomography (ZSFZTT), where a weighting factor is adopted to shrink the Fresnel zone during the inversion process. Numerical phantom and in vivo breast experiments were performed with ZSFZTT, FZTT, and RTT. In the numerical experiment, the reconstruction biases of size by ZSFZTT, FZTT, and RTT were 0.2%~8.3%, 2.3%~31.7%, and 1.8%~25%; the reconstruction biases of relative SS value by ZSFZTT, FZTT, and RTT were 24.7%~42%, 53%~60.8%, and 30.3%~47.8%; and the CNR by ZSFZTT, FZTT, and RTT were 67.7~96.6, 68.5~98, and 1.7~2.7. In the in vivo breast experiment, ZSFZTT provided the highest CNR of 8.6 compared to 8.1 by FZTT and 1.9 by RTT. ZSFZTT improved the reconstruction accuracy of size and the relative reconstruction accuracy of SS value compared to FZTT and RTT while maintaining a high CNR similar to that of FZTT.

Segmentation-free x-ray energy spectrum estimation for computed tomography using dual-energy material decomposition.

An x-ray energy spectrum plays an essential role in computed tomography (CT) imaging and related tasks. Because of the high photon flux of clinical CT scanners, most of the spectrum estimation methods are indirect and usually suffer from various limitations. In this study, we aim to provide a segmentation-free, indirect transmission measurement-based energy spectrum estimation method using dual-energy material decomposition. The general principle of this method is to minimize the quadratic error between the polychromatic forward projection and the raw projection to calibrate a set of unknown weights, which are used to express the unknown spectrum together with a set of model spectra. The polychromatic forward projection is performed using material-specific images, which are obtained using dual-energy material decomposition. The algorithm was evaluated using numerical simulations, experimental phantom data, and realistic patient data. The results show that the estimated spectrum matches the reference spectrum quite well and the method is robust. Extensive studies suggest that the method provides an accurate estimate of the CT spectrum without dedicated physical phantom and prolonged workflow. This paper may be attractive for CT dose calculation, artifacts reduction, polychromatic image reconstruction, and other spectrum-involved CT applications.

A variation-based ring artifact correction method with sparse constraint for flat-detector CT.

The reconstructed slice quality of flat-detector computed tomography (CT) is often disturbed by concentric-ring artifacts. Since concentric rings in CT slices appear as straight stripes when transformed into polar coordinates, a variation-based model is proposed to suppress the stripes. The method is motivated by two observations about stripes in polar coordinates: (1) ring artifacts attenuate gradually along the radial direction, leading to a sparse distribution of stripes and (2) stripes greatly distort the image gradient across the stripes, while slightly affecting the image gradient along the stripes. Thus, a [Formula: see text]-norm-based data fidelity term and a [Formula: see text]-norm/[Formula: see text]-norm unidirectional variation-based regularization term are presented to characterize the stripes. The alternating direction method of multipliers is introduced to solve the resulting minimization problem. Moreover, we discuss the interpolation methods used in coordinate transformation and find that the nearest neighbor interpolation is optimal. Experimental results on simulated and real data demonstrate that our method can correct ring artifacts effectively compared with state-of-the-art coordinate transformation-based methods, as well as preserve the structures and details of slices.