Ultrasound computed tomography based on full waveform inversion with source directivity calibration

Ultrasound computed tomography based on full waveform inversion has the potential to provide high-resolution images of human tissues in a quantitative manner. A successful ultrasound computed tomography system requires the decent knowledge of acquisition array, including the spatial position and the directivity of each transducer, to meet the high-level demand of clinical applications. The conventional full waveform inversion algorithm assumes a point source with the omni-directional emission. Such assumption does not hold when the directivity of emitting transducer is not negligible. For a practical implementation, an efficient and accurate self-checking evaluation of directivity is crucial prior to the reconstruction of images. We propose to measure the directivity of each emitting transducer using the full-matrix captured data obtained with a water-immersed and target-free experiment. We introduce the weighted virtual point-source array to act as the proxy of emitting transducer during the numerical simulation. The weights of different points in the virtual array can be calculated from the observed data using the gradient-based local optimization method. Although the full waveform imaging method relies on the finite-difference solver of wave equation, such directivity estimation benefits from the introduction of analytical solver. The trick significantly reduces the numerical cost, enabling an automatic directivity self-check at boot. We verify the feasibility, efficiency, and accuracy of the virtual array method through simulated and experimental tests. For the experimental test, we also illustrate that full waveform inversion with directivity calibration can reduce the artifacts introduced by the conventional point source assumption, improving the quality of reconstructed images..