We discuss a bent-ray ultrasound tomography algorithm with total-variation (TV) regularization. We have applied this algorithm to 61 in vivo breast datasets collected with our in-house clinical prototype for imaging sound-speed distributions in the breast. Our analysis showed that TV regularization could preserve sharper lesion edges than the classic Tikhonov regularization. Furthermore, the image quality of our TV bent-ray sound-speed tomograms was superior to that of the straight-ray counterparts for all types of breasts within BI-RADS density categories 1 through 4. Our analysis showed that the improvements for average sharpness (in the unit of (m · s)−1) of lesion edges in our TV bent-ray tomograms are between 2.1 to 3.4-fold compared with the straight ray tomograms. Reconstructed sound-speed tomograms illustrated that our algorithm could successfully image fatty and glandular tissues within the breast. We calculated the mean sound-speed values for fatty tissue and breast parenchyma as 1422±9 m/s (mean±SD) and 1487±21 m/s, respectively. Based on 32 lesions in a cohort of 61 patients, we also found that the mean sound-speed for malignant breast lesions 1548±17 m/s was higher, on average, than that of benign ones (1513±27 m/s) (one-sided p < 0.001). These results suggest that, clinically, sound-speed tomograms can be used to assess breast density (and therefore, breast cancer risk), as well as detect and help differentiate breast lesions. Finally, our sound-speed tomograms may also be a useful tool to monitor the clinical response of breast cancer patients to neo-adjuvant chemotherapy.

Objective and motivation-Time-of-flight (TOF) tomography used by a clinical ultrasound tomography device can efficiently and reliably produce sound-speed images of the breast for cancer diagnosis. Accurate picking of TOFs of transmitted ultrasound signals is extremely important to ensure high-resolution and high-quality ultrasound sound-speed tomograms. Since manually picking is time-consuming for large datasets, we developed an improved automatic TOF picker based on the Akaike information criterion (AIC), as described in this paper.Methods-We make use of an approach termed multi-model inference (model averaging), based on the calculated AIC values, to improve the accuracy of TOF picks. By using multi-model inference, our picking method incorporates all the information near the TOF of ultrasound signals. Median filtering and reciprocal pair comparison are also incorporated in our AIC picker to effectively remove outliers.Results-We validate our AIC picker using synthetic ultrasound waveforms, and demonstrate that our automatic TOF picker can accurately pick TOFs in the presence of random noise with absolute amplitudes up to 80% of the maximum absolute signal amplitude. We apply the new method to 1160 in vivo breast ultrasound waveforms, and compare the picked TOFs with manual picks and amplitude threshold picks. The mean value and standard deviation between our TOF picker and manual picking are 0.4 μs and 0.29 μs, while for amplitude threshold picker the values are 1.02 μs and 0.9 μs, respectively. Tomograms for in vivo breast data with high signal-to-noise ratio (SNR) (~25 dB) and low SNR (~18 dB) clearly demonstrate that our AIC picker is much less sensitive to the SNRs of the data, compared to the amplitude threshold picker.Discussion and conclusions-The picking routine developed here is aimed at determining reliable quantitative values, necessary for adding diagnostic information to our clinical ultrasound tomography device -CURE. It has been successfully adopted into CURE, and allows us to generate such values reliably. We demonstrate that in vivo sound-speed tomograms with our TOF picks significantly improve the reconstruction accuracy and reduce image artifacts.

No abstract

No abstract

A migration approach based on a local application of the Born approximation within each extrapolation interval contains a singularity that can make direct application unstable. Previous authors have suggested adding an imaginary part to the vertical wavenumber to eliminate the singularity. However, their method requires that the reference slowness must be the maximum slowness of a given layer; consequently, the slowness perturbations are larger than those when the average slowness is selected as a reference slowness. Therefore, its applicability is limited. We develop an extended local Born Fourier migration method that circumvents the singularity problem of the local Born solution and makes it possible to choose the average slowness as a reference slowness. It is computationally efficient because of the use of a fast Fourier transform algorithm. It can handle wider angles (or steeper interfaces) and scattering effects of heterogeneities more accurately than the split-step Fourier (SSF) method, which accounts for only the phase change as a result of the slowness perturbations but not amplitude change. To handle large lateral slowness variations, we introduce different reference slownesses in different regions of a medium to ensure the condition of small perturbation. The migration result obtained using the extended local Born Fourier method with multiple reference slownesses demonstrates that the method can produce high-quality images of complex structures with large lateral slowness variations.

No abstract

No abstract

We develop a method to evaluate the resolution of seismic images using a wave-equation-based operator. The one-way wave-equation-based propagator is used to extrapolate the wavefields from both the source and receivers to a target region. The plane wave analysis method is used to decompose the wavefields and calculate the wavenumber domain resolution function. The method can be applied to complex models without smoothing velocity. It can also handle irregular acquisition geometry and finite frequency band. The resulting resolution function can be used to estimate image quality or correct the image distortion caused by the uneven illumination.