1
|
Ali R, Mitcham TM, Brevett T, Agudo OC, Martinez CD, Li C, Doyley MM, Duric N. 2-D Slicewise Waveform Inversion of Sound Speed and Acoustic Attenuation for Ring Array Ultrasound Tomography Based on a Block LU Solver. IEEE TRANSACTIONS ON MEDICAL IMAGING 2024; 43:2988-3000. [PMID: 38564345 PMCID: PMC11294001 DOI: 10.1109/tmi.2024.3383816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Ultrasound tomography is an emerging imaging modality that uses the transmission of ultrasound through tissue to reconstruct images of its mechanical properties. Initially, ray-based methods were used to reconstruct these images, but their inability to account for diffraction often resulted in poor resolution. Waveform inversion overcame this limitation, providing high-resolution images of the tissue. Most clinical implementations, often directed at breast cancer imaging, currently rely on a frequency-domain waveform inversion to reduce computation time. For ring arrays, ray tomography was long considered a necessary step prior to waveform inversion in order to avoid cycle skipping. However, in this paper, we demonstrate that frequency-domain waveform inversion can reliably reconstruct high-resolution images of sound speed and attenuation without relying on ray tomography to provide an initial model. We provide a detailed description of our frequency-domain waveform inversion algorithm with open-source code and data that we make publicly available.
Collapse
|
2
|
Zhang N, Zhao Y, Yuan Y, Xiao Y, Qin M, Shen Y. Cross-correlation adjustment full-waveform inversion with source encoding in ultrasound computed tomography. ULTRASONICS 2024; 142:107392. [PMID: 38991429 DOI: 10.1016/j.ultras.2024.107392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Revised: 05/21/2024] [Accepted: 06/26/2024] [Indexed: 07/13/2024]
Abstract
Full-waveform inversion (FWI) is one of the leading-edge techniques in ultrasound computed tomography (USCT). FWI reconstructs the images of sound speed by iteratively minimizing the difference between the predicted and measured signals. The challenges of FWI are to improve its stability and reduce its computational cost. In this paper, a new USCT algorithm based on cross-correlation adjustment FWI with source encoding (CCAFWI-SE) is proposed. In this algorithm, the gradient is adjusted using the intermediate signals as the inversion target rather than the measured signals during iteration. The intermediate signals are generated using the travel time difference calculated by cross-correlation. In the case of conventional FWI failure, using the proposed algorithm, the estimated sound speed can converge toward the ground truth. To reduce the computational cost, an intermittent update strategy is implemented. This strategy only requires one time for the calculation of the travel time difference per stage, so that the source encoding can be used. Simulation and laboratory experiments are implemented to validate this approach. The experiment results show it has successfully recovered the sound speed model, while conventional FWI failed when the initial model greatly differed from the ground truth. This verifies that our approach improves the stability of the reconstruction in USCT. In practice, additional computational costs can be reduced by combining our approach with existing methods. The proposed approach increases the robustness of the FWI and expands its application.
Collapse
Affiliation(s)
- Nuomin Zhang
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China
| | - Yue Zhao
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China.
| | - Yu Yuan
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China
| | - Yang Xiao
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China
| | - Mengting Qin
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China
| | - Yi Shen
- Harbin Institute of Technology, No 92 Dazhi Street, Harbin, 150001, Heilongjiang Province, China
| |
Collapse
|
3
|
Wen J, Jiang C, Chen H. Detection of Multi-Layered Bond Delamination Defects Based on Full Waveform Inversion. SENSORS (BASEL, SWITZERLAND) 2024; 24:4017. [PMID: 38931799 PMCID: PMC11207941 DOI: 10.3390/s24124017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2024] [Revised: 06/17/2024] [Accepted: 06/18/2024] [Indexed: 06/28/2024]
Abstract
This study aimed to address the challenges encountered in traditional bulk wave delamination detection methods characterized by low detection efficiency. Additionally, the limitations of guided wave delamination detection methods were addressed, particularly those utilizing reflected waves, which are susceptible to edge reflections, thus complicating effective defect extraction. Leveraging the full waveform inversion algorithm, an innovative approach was established for detecting delamination defects in multi-layered structures using ultrasonic guided wave arrays. First, finite element modeling was employed to simulate guided wave data acquisition by a circular array within an aluminum-epoxy bilayer structure with embedded delamination defects. Subsequently, the full waveform inversion algorithm was applied to reconstruct both regular and irregular delamination defects. Analysis results indicated the efficacy of the proposed approach in accurately identifying delamination defects of varying shapes. Furthermore, an experimental platform for guided wave delamination defect detection was established, and experiments were conducted on a steel-cement bilayer structure containing an irregular delamination defect. The experimental results validated the exceptional imaging precision of our proposed technique for identifying delamination defects in multi-layered boards. In summary, the proposed method can accurately determine both the positions and sizes of defects with higher detection efficiency than traditional pulse-echo delamination detection methods.
Collapse
Affiliation(s)
- Jiawei Wen
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China;
- School of Electronic, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Can Jiang
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China;
| | - Hao Chen
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China;
- School of Electronic, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
4
|
Ren J, Li J, Liu C, Chen S, Liang L, Liu Y. Deep Learning With Physics-Embedded Neural Network for Full Waveform Ultrasonic Brain Imaging. IEEE TRANSACTIONS ON MEDICAL IMAGING 2024; 43:2332-2346. [PMID: 38329866 DOI: 10.1109/tmi.2024.3363144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
The convenience, safety, and affordability of ultrasound imaging make it a vital non-invasive diagnostic technique for examining soft tissues. However, significant differences in acoustic impedance between the skull and soft tissues hinder the successful application of traditional ultrasound for brain imaging. In this study, we propose a physics-embedded neural network with deep learning based full waveform inversion (PEN-FWI), which can achieve reliable quantitative imaging of brain tissues. The network consists of two fundamental components: forward convolutional neural network (FCNN) and inversion sub-neural network (ISNN). The FCNN explores the nonlinear mapping relationship between the brain model and the wavefield, replacing the tedious wavefield calculation process based on the finite difference method. The ISNN implements the mapping from the wavefield to the model. PEN-FWI includes three iterative steps, each embedding the F CNN into the ISNN, ultimately achieving tomography from wavefield to brain models. Simulation and laboratory tests indicate that PEN-FWI can produce high-quality imaging of the skull and soft tissues, even starting from a homogeneous water model. PEN-FWI can achieve excellent imaging of clot models with constant uniform distribution of velocity, randomly Gaussian distribution of velocity, and irregularly shaped randomly distributed velocity. Robust differentiation can also be achieved for brain slices of various tissues and skulls, resulting in high-quality imaging. The imaging time for a horizontal cross-sectional imag e of the brain is only 1.13 seconds. This algorithm can effectively promote ultrasound-based brain tomography and provide feasible solutions in other fields.
Collapse
|
5
|
Soydan DA, Top CB, Gençer NG. On the utilization of the adjoint method in microwave tomography. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2024; 40:e3818. [PMID: 38556485 DOI: 10.1002/cnm.3818] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 03/14/2024] [Accepted: 03/17/2024] [Indexed: 04/02/2024]
Abstract
In microwave imaging, the adjoint method is widely used for the efficient calculation of the update direction, which is then used to update the unknown model parameter. However, the utilization and the formulation of the adjoint method differ significantly depending on the imaging scenario and the applied optimization algorithm. Because of the problem-specific nature of the adjoint formulations, the dissimilarities between the adjoint calculations may be overlooked. Here, we have classified the adjoint method formulations into two groups: the direct and indirect methods. The direct method involves calculating the derivative of the cost function, whereas, in the indirect method, the derivative of the predicted data is calculated. In this review, the direct and indirect adjoint methods are presented, compared, and discussed. The formulations are explicitly derived using the two-dimensional wave equation in frequency and time domains. Finite-difference time-domain simulations are conducted to show the different uses of the adjoint methods for both single source-multiple receiver, and multiple transceiver scenarios. This study demonstrated that an appropriate adjoint method selection is significant to achieve improved computational efficiency for the applied optimization algorithm.
Collapse
Affiliation(s)
- Damla Alptekin Soydan
- Electrical and Electronics Engineering Department, Middle East Technical University, Ankara, Turkey
- ASELSAN A.Ş, Ankara, Turkey
| | | | - Nevzat G Gençer
- Electrical and Electronics Engineering Department, Middle East Technical University, Ankara, Turkey
| |
Collapse
|
6
|
Jiang X, Gan K, Wang Y, Tao C, Liu X, Yuan J, Jin Z. Demonstration Study of Reflector-Based Volumetric Speed-of-Sound Imaging With Linear Ultrasound Arrays. ULTRASONIC IMAGING 2024; 46:186-196. [PMID: 38647142 DOI: 10.1177/01617346241246807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
Conventional B-mode ultrasound imaging has difficulty in delineating homogeneous soft tissues with similar acoustic impedances, as the reflectivity depends on the acoustic impedance at the interface. As a quantitative imaging biomarker sensitive to alteration of biomechanical properties, speed-of-sound (SoS) holds promising potential for tissue and disease differentiation such as delineation of different breast tissue types with similar acoustic impedance. Compared to two-dimensional (2D) SoS images, three-dimensional (3D) volumetric SoS images achieved through a full-angle ultrasound scan can reveal more intricate morphological structures of tissues; however, they generally require a ring transducer. In this study, we introduce a 3D SoS reconstruction system that utilizes hand-held linear arrays instead. This system employs a passive reflector positioned opposite the linear arrays, serving as an echogenic reference for time-of-flight (ToF) measurements, and a high-definition camera to track the location corresponding to each group of transmit-receive data. To merge these two streams of ToF measurements and location tracking, a voxel-based reconstruction algorithm is implemented. Experimental results with gelatin phantom and ex vivo tissue have demonstrated the stability of our proposed method. Moreover, the results underscore the potential of this system as a complementary diagnostic modality, particularly in the context of diseases such as breast cancer.
Collapse
Affiliation(s)
- Xiaoyi Jiang
- School of Electronic Science and Engineering, Nanjing University, Nanjing, China
| | - Kexin Gan
- School of Electronic Science and Engineering, Nanjing University, Nanjing, China
| | - Yuxin Wang
- School of Electronic Science and Engineering, Nanjing University, Nanjing, China
| | - Chao Tao
- School of Physics, Nanjing University, Nanjing, China
| | - Xiaojun Liu
- School of Physics, Nanjing University, Nanjing, China
| | - Jie Yuan
- School of Electronic Science and Engineering, Nanjing University, Nanjing, China
| | - Zhibin Jin
- Affilated Drum Tower Hospital, Medical School of Nanjing University, Nanjing, China
| |
Collapse
|
7
|
Li F, Villa U, Duric N, Anastasio MA. A Forward Model Incorporating Elevation-Focused Transducer Properties for 3-D Full-Waveform Inversion in Ultrasound Computed Tomography. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2023; 70:1339-1354. [PMID: 37682648 PMCID: PMC10775680 DOI: 10.1109/tuffc.2023.3313549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/10/2023]
Abstract
Ultrasound computed tomography (USCT) is an emerging medical imaging modality that holds great promise for improving human health. Full-waveform inversion (FWI)-based image reconstruction methods account for the relevant wave physics to produce high spatial resolution images of the acoustic properties of the breast tissues. A practical USCT design employs a circular ring-array comprised of elevation-focused ultrasonic transducers, and volumetric imaging is achieved by translating the ring-array orthogonally to the imaging plane. In commonly deployed slice-by-slice (SBS) reconstruction approaches, the 3-D volume is reconstructed by stacking together 2-D images reconstructed for each position of the ring-array. A limitation of the SBS reconstruction approach is that it does not account for 3-D wave propagation physics and the focusing properties of the transducers, which can result in significant image artifacts and inaccuracies. To perform 3-D image reconstruction when elevation-focused transducers are employed, a numerical description of the focusing properties of the transducers should be included in the forward model. To address this, a 3-D computational model of an elevation-focused transducer is developed to enable 3-D FWI-based reconstruction methods to be deployed in ring-array-based USCT. The focusing is achieved by applying a spatially varying temporal delay to the ultrasound pulse (emitter mode) and recorded signal (receiver mode). The proposed numerical transducer model is quantitatively validated and employed in computer simulation studies that demonstrate its use in image reconstruction for ring-array USCT.
Collapse
|
8
|
Drainville RA, Chatillon S, Moore D, Snell J, Padilla F, Lafon C. A simulation study on the sensitivity of transcranial ray-tracing ultrasound modeling to skull properties. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2023; 154:1211-1225. [PMID: 37610718 DOI: 10.1121/10.0020761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 08/08/2023] [Indexed: 08/24/2023]
Abstract
In transcranial focused ultrasound therapies, such as treating essential tremor via thermal ablation in the thalamus, acoustic energy is focused through the skull using a phased-array transducer. Ray tracing is a computationally efficient method that can correct skull-induced phase aberrations via per-element phase delay calculations using patient-specific computed tomography (CT) data. However, recent studies show that variations in CT-derived Hounsfield unit may account for only 50% of the speed of sound variability in human skull specimens, potentially limiting clinical transcranial ultrasound applications. Therefore, understanding the sensitivity of treatment planning methods to material parameter variations is essential. The present work uses a ray-tracing simulation model to explore how imprecision in model inputs, arising from clinically significant uncertainties in skull properties or considerations of acoustic phenomena, affects acoustic focusing quality through the skull. We propose and validate new methods to optimize ray-tracing skull simulations for clinical treatment planning, relevant for predicting intracranial target's thermal rise, using experimental data from ex-vivo human skulls.
Collapse
Affiliation(s)
| | | | - David Moore
- Focused Ultrasound Foundation, Charlottesville, Virginia 22903, USA
| | - John Snell
- Histosonics, Ann Arbor, Michigan 48103, USA
| | - Frederic Padilla
- Focused Ultrasound Foundation, Charlottesville, Virginia 22903, USA
| | - Cyril Lafon
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, Lyon, France
| |
Collapse
|
9
|
Zhou C, Xu K, Ta D. Frequency-domain full-waveform inversion-based musculoskeletal ultrasound computed tomography. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2023; 154:279-294. [PMID: 37449785 DOI: 10.1121/10.0020151] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 06/27/2023] [Indexed: 07/18/2023]
Abstract
Recently, full-waveform inversion (FWI) has become a promising tool for ultrasound computed tomography (USCT). However, as a computationally intensive technique, FWI suffers from computational burden, especially in conventional time-domain full-waveform inversion (TDFWI). On the contrary, frequency-domain full-waveform inversion (FDFWI) provides a relatively high computational efficiency as the propagation of discrete frequencies is much cheaper than full time-domain modeling. FDFWI has already been applied in soft tissue imaging, such as breast, but for the musculoskeletal model with high impedance contrast between hard and soft tissues, there is still a lack of an effective source estimation method. In this paper, a water-referenced data calibration method is proposed to address the source estimation challenge in the presence of bones, which achieves consistency between the measured and simulated data before the FDFWI procedure. To avoid the cycle-skipping local minimum effect and facilitate the algorithm convergence, a starting frequency criterion for musculoskeletal FDFWI is further proposed. The feasibility of the proposed method is demonstrated by numerical studies on retrieving the anatomies of the leg models and different musculoskeletal lesions. The study extends the advanced FDFWI method to the musculoskeletal system and provides an alternative solution for musculoskeletal USCT imaging with high computational efficiency.
Collapse
Affiliation(s)
- Chenchen Zhou
- Center for Biomedical Engineering, School of Information Science and Technology, Fudan University, Shanghai, China
| | - Kailiang Xu
- Center for Biomedical Engineering, School of Information Science and Technology, Fudan University, Shanghai, China
| | - Dean Ta
- Center for Biomedical Engineering, School of Information Science and Technology, Fudan University, Shanghai, China
| |
Collapse
|
10
|
Yuan Y, Zhao Y, Zhang N, Xiao Y, Jin J, Feng N, Shen Y. Full-Waveform Inversion for Breast Ultrasound Tomography Using Line-Shape Modeled Elements. ULTRASOUND IN MEDICINE & BIOLOGY 2023; 49:1070-1081. [PMID: 36737306 DOI: 10.1016/j.ultrasmedbio.2022.12.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 12/06/2022] [Accepted: 12/10/2022] [Indexed: 05/11/2023]
Abstract
OBJECTIVE The objective of the work described here was to incorporate the spatial shapes of the transducer elements into the framework of the full-waveform inversion. METHODS An element is treated as its cross-section in the 2-D imaging plane, that is, a line segment. The elements are not simply modeled as a set of point sources on their surface to avoid staircasing artifacts. By use of the Fourier collocation method, an element is spatially represented as the discrete convolution between its spatial distribution and a band-limited delta function. The excitation pulses on the emitters and recorded signals on the receivers are then weighted based on the discrete convolution results. Digital and physical experiments are implemented to validate the method. DISCUSSION It is meaningful to model the shapes of the elements if their spatial sizes are similar to or larger than the acoustic wavelengths. It should, however, be noted that because this article focuses on 2-D imaging, the inter-plane effects are not considered. CONCLUSION The approach helps reduce the root mean square errors and increase the structural similarity of the reconstructed images. It also helps to improve the stability of convergence and to accelerate the convergence speed.
Collapse
Affiliation(s)
- Yu Yuan
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China
| | - Yue Zhao
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China.
| | - Nuomin Zhang
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China
| | - Yang Xiao
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China
| | - Jing Jin
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China
| | - Naizhang Feng
- Control Theory and Engineering, School of Astronautics, Harbin Institute of Technology, Harbin, China
| | - Yi Shen
- Shenzhen Engineering Lab for Medical Intelligent Wireless Ultrasonic Imaging Technology, Harbin Institute of Technology, Harbin, China
| |
Collapse
|
11
|
Wu X, Li Y, Su C, Li P, Wang X, Lin W. Ultrasound computed tomography based on full waveform inversion with source directivity calibration. ULTRASONICS 2023; 132:107004. [PMID: 37071945 DOI: 10.1016/j.ultras.2023.107004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/04/2022] [Revised: 03/06/2023] [Accepted: 03/31/2023] [Indexed: 05/03/2023]
Abstract
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..
Collapse
Affiliation(s)
- Xiaoqing Wu
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yubing Li
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China.
| | - Chang Su
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Panpan Li
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xiangda Wang
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China; Ruyuan Yao Autonomous Dongyangguang Industrial Development Co. Ltd, Shaoguan 512721, China
| | - Weijun Lin
- Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China.
| |
Collapse
|
12
|
Hu P, Li L, Wang LV. Location-Dependent Spatiotemporal Antialiasing in Photoacoustic Computed Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2023; 42:1210-1224. [PMID: 36449587 PMCID: PMC10171137 DOI: 10.1109/tmi.2022.3225565] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/09/2023]
Abstract
Photoacoustic computed tomography (PACT) images optical absorption contrast by detecting ultrasonic waves induced by optical energy deposition in materials such as biological tissues. An ultrasonic transducer array or its scanning equivalent is used to detect ultrasonic waves. The spatial distribution of the transducer elements must satisfy the spatial Nyquist criterion; otherwise, spatial aliasing occurs and causes artifacts in reconstructed images. The spatial Nyquist criterion poses different requirements on the transducer elements' distributions for different locations in the image domain, which has not been studied previously. In this research, we elaborate on the location dependency through spatiotemporal analysis and propose a location-dependent spatiotemporal antialiasing method. By applying this method to PACT in full-ring array geometry, we effectively mitigate aliasing artifacts with minimal effects on image resolution in both numerical simulations and in vivo experiments.
Collapse
|
13
|
Qu X, Ren C, Yan G, Zheng D, Tang W, Wang S, Lin H, Zhang J, Jiang J. Deep-Learning-Based Ultrasound Sound-Speed Tomography Reconstruction with Tikhonov Pseudo-Inverse Priori. ULTRASOUND IN MEDICINE & BIOLOGY 2022; 48:2079-2094. [PMID: 35922265 PMCID: PMC10448397 DOI: 10.1016/j.ultrasmedbio.2022.05.033] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 05/26/2022] [Accepted: 05/31/2022] [Indexed: 06/15/2023]
Abstract
Ultrasound sound-speed tomography (USST) is a promising technology for breast imaging and breast cancer detection. Its reconstruction is a complex non-linear mapping from the projection data to the sound-speed image (SSI). The traditional reconstruction methods include mainly the ray-based methods and the waveform-based methods. The ray-based methods with linear approximation have low computational cost but low reconstruction quality; the full wave-based methods with the complex non-linear model have high quality but high cost. To achieve both high quality and low cost, we introduced traditional linear approximation as prior knowledge into a deep neural network and treated the complex non-linear mapping of USST reconstruction as a combination of linear mapping and non-linear mapping. In the proposed method, the linear mapping was seamlessly implemented with a fully connected layer and initialized using the Tikhonov pseudo-inverse matrix. The non-linear mapping was implemented using a U-shape Net (U-Net). Furthermore, we proposed the Tikhonov U-shape net (TU-Net), in which the linear mapping was done before the non-linear mapping, and the U-shape Tikhonov net (UT-Net), in which the non-linear mapping was done before the linear mapping. Moreover, we conducted simulations and experiments for evaluation. In the numerical simulation, the root-mean-squared error was 6.49 and 4.29 m/s for the UT-Net and TU-Net, the peak signal-to-noise ratio was 49.01 and 52.90 dB, the structural similarity was 0.9436 and 0.9761 and the reconstruction time was 10.8 and 11.3 ms, respectively. In this study, the SSIs obtained with the proposed methods exhibited high sound-speed accuracy. Both the UT-Net and the TU-Net achieved high quality and low computational cost.
Collapse
Affiliation(s)
- Xiaolei Qu
- School of Instrumentation and Optoelectronics Engineering, Beihang University, Beijing, China
| | - Chujian Ren
- School of Instrumentation and Optoelectronics Engineering, Beihang University, Beijing, China
| | - Guo Yan
- School of Instrumentation and Optoelectronics Engineering, Beihang University, Beijing, China
| | - Dezhi Zheng
- Research Institute for Frontier Science, Beihang University, Beijing, China
| | - Wenzhong Tang
- School of Computer Science and Engineering, Beihang University, Beijing, China
| | - Shuai Wang
- Research Institute for Frontier Science, Beihang University, Beijing, China
| | - Hongxiang Lin
- Research Center for Healthcare Data Science, Zhejiang Lab, Hangzhou, China
| | - Jingya Zhang
- School of Instrumentation and Optoelectronics Engineering, Beihang University, Beijing, China
| | - Jue Jiang
- Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, New York, USA.
| |
Collapse
|
14
|
Sarno D, Baker C, Curtis S, Hodnett M, Zeqiri B. In Vivo Measurements of the Bulk Ultrasonic Attenuation Coefficient of Breast Tissue Using a Novel Phase-Insensitive Receiver. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:2943-2954. [PMID: 35976833 DOI: 10.1109/tuffc.2022.3198815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
This study describes the first in vivo acoustic attenuation measurements of breast tissue undertaken using a novel phase-insensitive detection technique employing a differential pyroelectric sensor. The operation of the sensor is thermal in nature, with its output signal being dictated by the acoustic power integrated over its surface. The particularly novel feature of the sensor lies in its differential principle of operation, which significantly enhances its immunity to background acoustic and vibration noise. A large area variant of the sensor was used to detect ultrasonic energy generated by an array of 14 discrete 3.2-MHz plane piston transducers, transmitted through pendent breasts in water. The transduction and reception capability represent key parts of a prototype Quantitative Ultrasound Computed Tomography Test Facility developed at the National Physical Laboratory to study the efficacy of phase-insensitive ultrasound computed tomography of breast phantoms containing a range of appropriate inclusions, in particular, the measurement uncertainties associated with quantitative reconstructions of the acoustic attenuation coefficient. For this study, attenuation coefficient measurements were made using 1-D projections on 12 nominally healthy study volunteers, whose age ranged from 19 to 65 years. Averaged or bulk attenuation coefficient values were generated in the range 1.7-4.6 dBcm -1 at 3.2 MHz and have been compared with existing literature, derived from in vivo and ex vivo studies. Results are encouraging and indicate that the relatively simple technique could be applied as a robust method for assessing the properties of breast tissue, particularly the balance of fatty (adipose) and fibroglandular components.
Collapse
|
15
|
Fogel H, Hughson M, Asefi M, Jeffrey I, LoVetri J. Generation of Prior Information in a Dual-Mode Microwave-Ultrasound Breast Imaging System. SENSORS (BASEL, SWITZERLAND) 2022; 22:7087. [PMID: 36146432 PMCID: PMC9502705 DOI: 10.3390/s22187087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Revised: 09/09/2022] [Accepted: 09/15/2022] [Indexed: 06/16/2023]
Abstract
A new breast imaging system capable of obtaining ultrasound and microwave scattered-field measurements with minimal or no movement of the breast between measurements has recently been reported. In this work, we describe the methodology that has been developed to generate prior information about the internal structures of the breast based on ultrasound data measured with the dual-mode system. This prior information, estimating both the geometry and complex-valued permittivity of tissues within the breast, is incorporated into the microwave inversion algorithm as a means of enhancing image quality. Several techniques to map reconstructed ultrasound speed to complex-valued relative permittivity are investigated. Quantitative images of two simplified dual-mode breast phantoms obtained using experimental data and the various forms of prior information are presented. Though preliminary, the results presented herein provide an understanding of the impacts of different forms of prior information on dual-mode reconstructions of the breast and can be used to inform future work on the subject.
Collapse
Affiliation(s)
- Hannah Fogel
- Department of Electrical and Computer Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
| | - Max Hughson
- Antec Controls, Winnipeg, MB R2K 3Z9, Canada
| | | | - Ian Jeffrey
- Department of Electrical and Computer Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
| | - Joe LoVetri
- Department of Electrical and Computer Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
| |
Collapse
|
16
|
Solgi R, Ghadiri H. Effects of Beamforming Techniques on Quality of Ultrasound Computed Tomography Images. J Biomed Phys Eng 2022; 12:349-358. [PMID: 36059289 PMCID: PMC9395622 DOI: 10.31661/jbpe.v0i0.2107-1367] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2021] [Accepted: 10/20/2021] [Indexed: 06/15/2023]
Abstract
BACKGROUND In body tissues, tumors generally have different speeds of sound (SOS) than normal tissues. In this respect, ultrasound computed tomography (UCT) can generate a cross-sectional SOS map as an innovative ultrasound imaging method. This technique can produce images with a resolution of millimeters and a high signal-to-noise ratio. OBJECTIVE This study aimed to improve UCT image quality without increasing breast cancer screening and diagnosis time. MATERIAL AND METHODS In this analytical study, a ring-shaped UCT breast imaging system was simulated using the K-wave toolbox of MATLAB. The system has a 20 cm diameter and 256 ultrasonic piezoelectrics placed in the ring's circumference. Different beamforming techniques imaged two designed phantoms (i.e., resolution and contrast), and the resolution and contrast to noise ratio (CNR) were calculated. RESULTS The results of resolution phantom imaging without any beamforming showed that only bars with the value of 0.125 and 0.167 lp/mm were distinguishable, and the 0.1 bars were not recognizable in the imaging. In addition, increasing the number of transmitters led to no noticeable change in resolution for 0.125 and 0.167 lp/mm bars. In all beamforming techniques for imaging the contrast phantom, the CNR parameter up to an object with a diameter of 8 mm increases with increasing diameter without any change. CONCLUSION The beamforming technique using three simultaneous transmitters improved the resolution by about 1 mm compared to the normal strategy. In addition to high-contrast images, beamforming with 9 simultaneous transmitters led to a preferable technique.
Collapse
Affiliation(s)
- Razieh Solgi
- MSc, Department of Medical Physics and Biomedical Engineering, Tehran University of Medical Sciences, Tehran, Iran
- MSc, Preclinical Lab, Core Facility, Tehran University of Medical Sciences, Tehran, Iran
| | - Hossein Ghadiri
- PhD, Department of Medical Physics and Biomedical Engineering, Tehran University of Medical Sciences, Tehran, Iran
- PhD, Advanced Medical Technologies and Equipment Institute, Research Center for Molecular and Cellular Imaging, Tehran University of Medical Sciences, Tehran, Iran
| |
Collapse
|
17
|
Rosa PTCR, Fontes-Pereira AJ, Grimal Q, Pereira WCDA. Femoral neck phantom imaging using time-domain topological energy method. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2022; 152:706. [PMID: 35931554 DOI: 10.1121/10.0012695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2021] [Accepted: 06/30/2022] [Indexed: 06/15/2023]
Abstract
Ultrasonic bone imaging is a complex task, primarily because of the low energy contained in the signals reflected from the internal bone structures. In this study, the reconstruction of a bone-mimicking phantom echographic image using time-domain topological energy (TDTE) is proposed. A TDTE image results from a combination of forward and adjoint fields. The first is a solution of a numerical model that reproduces the setup of the experimental data acquisition to the best extent possible. The second has similar characteristics, but the source term is the time-reversed residue between the forward field and signals obtained from the experiment. The acquisition-reconstruction system used a linear phased-array transducer with a 5 MHz center frequency to acquire the signals and was coupled with a k-wave toolbox to implement the numerical models and perform the image reconstruction. The results showed good agreement between the geometry of the real phantom and the ultrasonic images. However, thickness evaluation errors were observed, which may be due to incorrect assumptions about the velocity models throughout the medium, a priori assumed to be known. Thus, this method has shown promising results and should be applied to the real femoral neck as a long-term objective.
Collapse
Affiliation(s)
- Paulo Tadeu C R Rosa
- Laboratório de Ultrassom, Programa de Engenharia Biomédica-COPPE, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-914, Brazil
| | | | - Quentin Grimal
- Laboratoire d'Imagerie Biomédicale, Sorbonne Université, Paris 75006, France
| | | |
Collapse
|
18
|
Ulrich IE, Boehm C, Zunino A, Bösch C, Fichtner A. Diffuse ultrasound computed tomography. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2022; 151:3654. [PMID: 35778223 DOI: 10.1121/10.0011540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Accepted: 05/13/2022] [Indexed: 06/15/2023]
Abstract
An alternative approach to acquire transmission travel time data is proposed, exploiting the geometry of devices commonly used in ultrasound computed tomography for medical imaging or non-destructive testing with ultrasonic waves. The intent is to (i) shorten acquisition time for devices with a large number of emitters, (ii) to eliminate the calibration step, and (iii) to suppress instrument noise. Inspired by seismic ambient field interferometry, the method rests on the active excitation of diffuse ultrasonic wavefields and the extraction of deterministic travel time information by inter-station correlation. To reduce stochastic errors and accelerate convergence, ensemble interferograms are obtained by phase-weighted stacking of observed and computed correlograms, generated with identical realizations of random sources. Mimicking an imaging setup, the accuracy of the travel time measurements as a function of the number of emitters and random realizations can be assessed both analytically and with spectral-element simulations for phantoms mimicking the model parameter distribution. The results warrant tomographic reconstructions with straight- or bent-ray approaches, where the effect of inherent stochastic fluctuations can be made significantly smaller than the effect of subjective choices on regularisation. This work constitutes a first conceptual study and a necessary prelude to future implementations.
Collapse
Affiliation(s)
- Ines Elisa Ulrich
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Christian Boehm
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Andrea Zunino
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Cyrill Bösch
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Andreas Fichtner
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| |
Collapse
|
19
|
Prasad S, Almekkawy M. DeepUCT: Complex cascaded deep learning network for improved ultrasound tomography. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac5296] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2021] [Accepted: 02/07/2022] [Indexed: 11/12/2022]
Abstract
Abstract
Ultrasound computed tomography is an inexpensive and radiation-free medical imaging technique used to quantify the tissue acoustic properties for advanced clinical diagnosis. Image reconstruction in ultrasound tomography is often modeled as an optimization scheme solved by iterative methods like full-waveform inversion. These iterative methods are computationally expensive, while the optimization problem is ill-posed and nonlinear. To address this problem, we propose to use deep learning to overcome the computational burden and ill-posedness, and achieve near real-time image reconstruction in ultrasound tomography. We aim to directly learn the mapping from the recorded time-series sensor data to a spatial image of acoustical properties. To accomplish this, we develop a deep learning model using two cascaded convolutional neural networks with an encoder–decoder architecture. We achieve a good representation by first extracting the intermediate mapping-knowledge and later utilizing this knowledge to reconstruct the image. This approach is evaluated on synthetic phantoms where simulated ultrasound data are acquired from a ring of transducers surrounding the region of interest. The measurement data is acquired by forward modeling the wave equation using the k-wave toolbox. Our simulation results demonstrate that our proposed deep-learning method is robust to noise and significantly outperforms the state-of-the-art traditional iterative method both quantitatively and qualitatively. Furthermore, our model takes substantially less computational time than the conventional full-wave inversion method.
Collapse
|
20
|
Fincke J, Zhang X, Shin B, Ely G, Anthony BW. Quantitative Sound Speed Imaging of Cortical Bone and Soft Tissue: Results From Observational Data Sets. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:502-514. [PMID: 34570702 DOI: 10.1109/tmi.2021.3115790] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
This work presents the first quantitative ultrasonic sound speed images of ex vivo limb cross-sections containing both soft tissue and bone using Full Waveform Inversion (FWI) with level set (LS) and travel time regularization. The estimated bulk sound speed of bone and soft tissue are within 10% and 1%, respectively, of ground truth estimates. The sound speed imagery shows muscle, connective tissue and bone features. Typically, ultrasound tomography (UST) using FWI is applied to imaging breast tissue properties (e.g. sound speed and density) that correlate with cancer. With further development, UST systems have the potential to deliver volumetric operator independent tissue property images of limbs with non-ionizing and portable hardware platforms. This work addresses the algorithmic challenges of imaging the sound speed of bone and soft tissue by combining FWI with LS regularization and travel time methods to recover soft tissue and bone sound speed with improved accuracy and reduced soft tissue artifacts when compared to conventional FWI. The value of leveraging LS and travel time methods is realized by evidence of improved bone geometry estimates as well as promising convergence properties and reduced risk of final model errors due to un-modeled shear wave propagation. Ex vivo bulk measurements of sound speed and MRI cross-sections validates the final inversion results.
Collapse
|
21
|
Ali R, Telichko AV, Wang H, Sukumar UK, Vilches-Moure JG, Paulmurugan R, Dahl JJ. Local Sound Speed Estimation for Pulse-Echo Ultrasound in Layered Media. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2022; 69:500-511. [PMID: 34723801 PMCID: PMC9127706 DOI: 10.1109/tuffc.2021.3124479] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Our previous methodology in local sound speed estimation utilized time delays measured by the cross correlation of delayed full-synthetic aperture channel data to estimate the average speed of sound. However, focal distortions in this methodology lead to biased estimates of the average speed of sound, which, in turn, leads to biased estimates of the local speed of sound. Here, we demonstrate the bias in the previous methodology and introduce a coherence-based average sound speed estimator that eliminates this bias and is computationally much cheaper in practice. Because this coherence-based approach estimates the average sound speed in the medium over an equally spaced grid in depth rather than time, we derive a refined model that relates the local and average speeds of sound as a function of depth in layered media. A fast, closed-form inversion of this model yields highly accurate local sound speed estimates. The root-mean-square (rms) error of local sound speed reconstruction in simulations of two-layer media is 4.6 and 2.5 m/s at 4 and 8 MHz, respectively. This work examines the impact of frequency, f -number, aberration, and reverberation on sound speed estimation. Phantom and in vivo experiments in rats further validate the coherence-based sound speed estimator.
Collapse
|
22
|
Lasaygues P, Espinosa L, Bernard S, Petit P, Guillermin R. Ultrasound Computed Tomography. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2022; 1364:227-250. [DOI: 10.1007/978-3-030-91979-5_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
23
|
An Ultrasound Tomography Method for Monitoring CO 2 Capture Process Involving Stirring and CaCO 3 Precipitation. SENSORS 2021; 21:s21216995. [PMID: 34770301 PMCID: PMC8587525 DOI: 10.3390/s21216995] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 10/16/2021] [Accepted: 10/19/2021] [Indexed: 11/17/2022]
Abstract
In this work, an ultrasound computed tomography (USCT) system was employed to investigate the fast-kinetic reactive crystallization process of calcium carbonate. USCT measurements and reconstruction provided key insights into the bulk particle distribution inside the stirred tank reactor and could be used to estimate the settling rate and settling time of the particles. To establish the utility of the USCT system for dynamical crystallization processes, first, the experimental imaging tasks were carried out with the stirred solid beads, as well as the feeding and stirring of the CaCO3 crystals. The feeding region, the mixing process, and the particles settling time could be detected from USCT data. Reactive crystallization experiments for CO2 capture were then conducted. Moreover, there was further potential for quantitative characterization of the suspension density in this process. USCT-based reconstructions were investigated for several experimental scenarios and operating conditions. This study demonstrates a real-time monitoring and fault detection application of USCT for reactive crystallization processes. As a robust noninvasive and nonintrusive tool, real-time signal analysis and reconstruction can be beneficial in the development of monitoring and control systems with real-world applications for crystallization processes. A diverse range of experimental studies shown here demonstrate the versatility of the USCT system in process application, hoping to unlock the commercial and industrial utility of the USCT devices.
Collapse
|
24
|
Pei Y, Zhang G, Zhang Y, Zhang W. Breast Acoustic Parameter Reconstruction Method Based on Capacitive Micromachined Ultrasonic Transducer Array. MICROMACHINES 2021; 12:963. [PMID: 34442585 PMCID: PMC8400655 DOI: 10.3390/mi12080963] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Revised: 08/07/2021] [Accepted: 08/10/2021] [Indexed: 11/17/2022]
Abstract
Ultrasound computed tomography (USCT) systems based on capacitive micromachined ultrasonic transducer (CMUT) arrays have a wide range of application prospects. For this paper, a high-precision image reconstruction method based on the propagation path of ultrasound in breast tissue are designed for the CMUT ring array; that is, time-reversal algorithms and FBP algorithms are respectively used to reconstruct sound speed distribution and acoustic attenuation distribution. The feasibility of this reconstruction method is verified by numerical simulation and breast model experiments. According to reconstruction results, sound speed distribution reconstruction deviation can be reduced by 53.15% through a time-reversal algorithm based on wave propagation theory. The attenuation coefficient distribution reconstruction deviation can be reduced by 61.53% through FBP based on ray propagation theory. The research results in this paper will provide key technological support for a new generation of ultrasound computed tomography systems.
Collapse
Affiliation(s)
| | - Guojun Zhang
- State Key Laboratory of Dynamic Testing Technology, North University of China, Taiyuan 030051, China; (Y.P.); (Y.Z.); (W.Z.)
| | | | | |
Collapse
|
25
|
Top CB. A Generalized Split-Step Angular Spectrum Method for Efficient Simulation of Wave Propagation in Heterogeneous Media. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2021; 68:2687-2696. [PMID: 33891551 DOI: 10.1109/tuffc.2021.3075367] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Angular spectrum (AS) methods enable efficient calculation of wave propagation from one plane to another inside homogeneous media. For wave propagation in heterogeneous media such as biological tissues, AS methods cannot be applied directly. Split-stepping techniques decompose the heterogeneous domain into homogeneous and perturbation parts, and provide a solution for forward wave propagation by propagating the incident wave in both frequency-space and frequency-wavenumber domains. Recently, a split-step hybrid angular spectrum (HAS) method was proposed for plane wave propagation of focused ultrasound beams. In this study, we extend these methods to enable simulation of acoustic pressure field for an arbitrary source distribution, by decomposing the source and reflection spectra into orthogonal propagation direction components, propagating each component separately, and summing all components to get the total field. We show that our method can efficiently simulate the pressure field of arbitrary sources in heterogeneous media. The accuracy of the method was analyzed comparing the resultant pressure field with pseudospectral time domain (PSTD) solution for breast tomography and hemispherical transcranial-focused ultrasound simulation models. Eighty times acceleration was achieved for a 3-D breast simulation model compared to PSTD solution with 0.005 normalized root mean-squared difference (NRMSD) between two solutions. For the hemispherical phased array, aberrations due to skull were accurately calculated in a single simulation run as evidenced by the resultant-focused ultrasound beam simulations, which had 0.001 NRMSD with 40 times acceleration factor compared to the PSTD method.
Collapse
|
26
|
Robins T, Camacho J, Agudo OC, Herraiz JL, Guasch L. Deep-Learning-Driven Full-Waveform Inversion for Ultrasound Breast Imaging. SENSORS 2021; 21:s21134570. [PMID: 34283105 PMCID: PMC8272012 DOI: 10.3390/s21134570] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Revised: 06/22/2021] [Accepted: 06/25/2021] [Indexed: 02/07/2023]
Abstract
Ultrasound breast imaging is a promising alternative to conventional mammography because it does not expose women to harmful ionising radiation and it can successfully image dense breast tissue. However, conventional ultrasound imaging only provides morphological information with limited diagnostic value. Ultrasound computed tomography (USCT) uses energy in both transmission and reflection when imaging the breast to provide more diagnostically relevant quantitative tissue properties, but it is often based on time-of-flight tomography or similar ray approximations of the wave equation, resulting in reconstructed images with low resolution. Full-waveform inversion (FWI) is based on a more accurate approximation of wave-propagation phenomena and can consequently produce very high resolution images using frequencies below 1 megahertz. These low frequencies, however, are not available in most USCT acquisition systems, as they use transducers with central frequencies well above those required in FWI. To circumvent this problem, we designed, trained, and implemented a two-dimensional convolutional neural network to artificially generate missing low frequencies in USCT data. Our results show that FWI reconstructions using experiment data after the application of the proposed method successfully converged, showing good agreement with X-ray CT and reflection ultrasound-tomography images.
Collapse
Affiliation(s)
- Thomas Robins
- Department of Earth Science and Engineering, Faculty of Engineering, Imperial College London, London SW7 2AZ, UK; (T.R.); (O.C.A.)
| | - Jorge Camacho
- Ultrasound Systems and Technology Group (GSTU), Institute for Physical and Information Technologies (ITEFI), Spanish National Research Council (CSIC), 28006 Madrid, Spain;
| | - Oscar Calderon Agudo
- Department of Earth Science and Engineering, Faculty of Engineering, Imperial College London, London SW7 2AZ, UK; (T.R.); (O.C.A.)
| | - Joaquin L. Herraiz
- Nuclear Physics Group and IPARCOS, Faculty of Physical Sciences, University Complutense of Madrid, CEI Moncloa, 28040 Madrid, Spain;
- Health Research Institute of the Hospital Clínico San Carlos (IdISSC), 28040 Madrid, Spain
| | - Lluís Guasch
- Department of Earth Science and Engineering, Faculty of Engineering, Imperial College London, London SW7 2AZ, UK; (T.R.); (O.C.A.)
- Correspondence:
| |
Collapse
|
27
|
Robins T, Camacho J, Agudo OC, Herraiz JL, Guasch L. Deep-Learning-Driven Full-Waveform Inversion for Ultrasound Breast Imaging. SENSORS 2021. [DOI: https://doi.org/10.3390/s21134570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Ultrasound breast imaging is a promising alternative to conventional mammography because it does not expose women to harmful ionising radiation and it can successfully image dense breast tissue. However, conventional ultrasound imaging only provides morphological information with limited diagnostic value. Ultrasound computed tomography (USCT) uses energy in both transmission and reflection when imaging the breast to provide more diagnostically relevant quantitative tissue properties, but it is often based on time-of-flight tomography or similar ray approximations of the wave equation, resulting in reconstructed images with low resolution. Full-waveform inversion (FWI) is based on a more accurate approximation of wave-propagation phenomena and can consequently produce very high resolution images using frequencies below 1 megahertz. These low frequencies, however, are not available in most USCT acquisition systems, as they use transducers with central frequencies well above those required in FWI. To circumvent this problem, we designed, trained, and implemented a two-dimensional convolutional neural network to artificially generate missing low frequencies in USCT data. Our results show that FWI reconstructions using experiment data after the application of the proposed method successfully converged, showing good agreement with X-ray CT and reflection ultrasound-tomography images.
Collapse
|
28
|
Arridge SR, Ehrhardt MJ, Thielemans K. (An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200205. [PMID: 33966461 DOI: 10.1098/rsta.2020.0205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
Imaging is omnipresent in modern society with imaging devices based on a zoo of physical principles, probing a specimen across different wavelengths, energies and time. Recent years have seen a change in the imaging landscape with more and more imaging devices combining that which previously was used separately. Motivated by these hardware developments, an ever increasing set of mathematical ideas is appearing regarding how data from different imaging modalities or channels can be synergistically combined in the image reconstruction process, exploiting structural and/or functional correlations between the multiple images. Here we review these developments, give pointers to important challenges and provide an outlook as to how the field may develop in the forthcoming years. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.
Collapse
Affiliation(s)
- Simon R Arridge
- Department of Computer Science, University College London, London, UK
| | - Matthias J Ehrhardt
- Department of Mathematical Sciences, University of Bath, Bath, UK
- Institute for Mathematical Innovation, University of Bath, Bath, UK
| | - Kris Thielemans
- Institute of Nuclear Medicine, University College London, London, UK
| |
Collapse
|
29
|
Quantitative Inspection of Complex-Shaped Parts Based on Ice-Coupled Ultrasonic Full Waveform Inversion Technology. APPLIED SCIENCES-BASEL 2021. [DOI: 10.3390/app11104433] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Ultrasonic methods have been extensively developed in nondestructive testing for various materials and components. However, accurately extracting quantitative information about defects still remains challenging, especially for complex structures. Although the immersion technique is commonly used for complex-shaped parts, the large mismatch of acoustic impedance between water and metal prevents effective ultrasonic transmission and leads to a low signal-to-noise ratio(SNR). In this paper, a quantitative imaging method is proposed for complex-shaped parts based on an ice-coupled full waveform inversion (FWI) method. Numerical experiments were carried out to quantitatively inspect the various defects in a turbine blade. Firstly, the k-space pseudospectral method was applied to simulate ice-coupled ultrasonic testing for the turbine blade. The recorded full wavefields were then applied for a frequency-domain FWI based on the Limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) method. With a carefully selected iterative number and frequency, a successive-frequency FWI can well detect half wavelength defects. Extended studies on an open notch with different orientations and multiple adjacent defects proved its capability to detect different types of defects. Finally, an uncertainty analysis was conducted with inaccurate initial velocity models with a relative error of ±2%, demonstrating its robustness even with a certain inaccuracy. This study demonstrates that the proposed method has a high potential to inspect complex-shaped structures with an excellent resolution.
Collapse
|
30
|
Espinosa L, Doveri E, Bernard S, Monteiller V, Guillermin R, Lasaygues P. Ultrasonic Imaging of High-contrasted Objects Based on Full-waveform Inversion: Limits under Fluid Modeling. ULTRASONIC IMAGING 2021; 43:88-99. [PMID: 33563137 DOI: 10.1177/0161734621990011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Quantitative ultrasound techniques have been previously used to evaluate biological hard tissues, characterized by a large acoustic impedance contrast. Here, we are interested in the imaging of experimental data from different test-targets with high acoustic impedance contrast, using the Full Waveform Inversion (FWI) method to solve the inverse problem. This method is based on high-resolution numerical modeling of the forward problem of interaction between waves and medium, considering the full time series. To reduce the complexity of the numerical implementation, the model considers a fluid medium. Therefore, the aim is to evaluate the precision of the reconstruction under this assumption for materials with a different level of attenuation of shear waves, to study the limits of this hypothesis. Images of the sound speed obtained using the experimental data are presented, and the precision of the reconstruction is evaluated. Future work should include viscoelastic materials.
Collapse
Affiliation(s)
- Luis Espinosa
- Aix-Marseille Université, CNRS, Centrale Marseille, LMA, Marseille, France
| | - Elise Doveri
- Aix-Marseille Université, CNRS, Centrale Marseille, LMA, Marseille, France
| | - Simon Bernard
- Laboratoire Ondes et Milieux Complexes, Université du Havre, Le Havre, France
| | - Vadim Monteiller
- Aix-Marseille Université, CNRS, Centrale Marseille, LMA, Marseille, France
| | - Régine Guillermin
- Aix-Marseille Université, CNRS, Centrale Marseille, LMA, Marseille, France
| | - Philippe Lasaygues
- Aix-Marseille Université, CNRS, Centrale Marseille, LMA, Marseille, France
| |
Collapse
|
31
|
Jiang Y, Peng C, Zhu Y, Ma X, Xu G, Yuan J, Wang X, Carson P. Biomedical Photoacoustic Imaging With Unknown Spatially Distributed Ultrasound Sensor Array. IEEE Trans Biomed Eng 2021; 68:2948-2956. [PMID: 33534699 DOI: 10.1109/tbme.2021.3056715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
OBJECTIVE With the growth of interest in different medical study on biological function, non-invasive photoacoustic imaging of biological tissue attracts the interests for researchers. To eliminate the limited angle effect of photoacoustic imaging based on ultrasound linear array, spatially distributed ultrasound sensor array is applied. The accurate sensor array position determines the quality of the imaging results. In this study, we proposed three methods based on photoacoustic and ultrasound signals to enhance the imaging quality using a 256-element full-ring array. METHODS Groups of photoacoustic and ultrasound signals are used to regress the position of each element sensor. RESULT In phantom study and mouse brain study, photoacoustic imaging results can both yield details clearly with average error rate of less than 1% (50 [Formula: see text]). CONCLUSION The performance of our three methods have proved that they can be potentially applied to other ultrasound-based medical imaging studies with unknown distributed positions of sensor array to enhance the imaging quality. SIGNIFICANCE The proposed methods can contribute to precise biomedical imaging with unknown distributed positions of sensor array in different application scenarios.
Collapse
|
32
|
Lafci B, Merčep E, Herraiz JL, Deán-Ben XL, Razansky D. Noninvasive multiparametric characterization of mammary tumors with transmission-reflection optoacoustic ultrasound. Neoplasia 2020. [DOI: https://doi.org/10.1016/j.neo.2020.10.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
33
|
Lafci B, Merčep E, Herraiz JL, Deán-Ben XL, Razansky D. Noninvasive multiparametric charac-terization of mammary tumors with transmission-reflection optoacoustic ultrasound. Neoplasia 2020; 22:770-777. [PMID: 33142241 PMCID: PMC7644559 DOI: 10.1016/j.neo.2020.10.008] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Revised: 10/09/2020] [Accepted: 10/12/2020] [Indexed: 02/07/2023]
Abstract
Development of imaging methods capable of furnishing tumor-specific morphological, functional, and molecular information is paramount for early diagnosis, staging, and treatment of breast cancer. Ultrasound (US) and optoacoustic (OA) imaging methods exhibit excellent traits for tumor imaging in terms of fast imaging speed, ease of use, excellent contrast, and lack of ionizing radiation. Here, we demonstrate simultaneous tomographic whole body imaging of optical absorption, US reflectivity, and speed of sound (SoS) in living mice. In vivo studies of 4T1 breast cancer xenografts models revealed synergistic and complementary value of the hybrid imaging approach for characterizing mammary tumors. While neovasculature surrounding the tumor areas were observed based on the vascular anatomy contrast provided by the OA data, the tumor boundaries could be discerned by segmenting hypoechoic structures in pulse-echo US images. Tumor delineation was further facilitated by enhancing the contrast and spatial resolution of the SoS maps with a full-wave inversion method. The malignant lesions could thus be distinguished from other hypoechoic regions based on the average SoS values. The reported findings corroborate the strong potential of the hybrid imaging approach for advancing cancer research in small animal models and fostering development of new clinical diagnostic approaches.
Collapse
Affiliation(s)
- Berkan Lafci
- Institute for Biomedical Engineering and Institute of Pharmacology and Toxicology, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - Elena Merčep
- Institute for Biological and Medical Imaging, Technical University of Munich and Helmholtz Center Munich, Neuherberg, Germany; iThera Medical GmbH, Munich, Germany
| | - Joaquin L Herraiz
- Nuclear Physics Group and IPARCOS, Complutense University of Madrid, Madrid, Spain; Health Research Institute of Hospital Clínico San Carlos (IdISSC), Madrid, Spain
| | - Xosé Luís Deán-Ben
- Institute for Biomedical Engineering and Institute of Pharmacology and Toxicology, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - Daniel Razansky
- Institute for Biomedical Engineering and Institute of Pharmacology and Toxicology, University of Zurich and ETH Zurich, Zurich, Switzerland; Institute for Biological and Medical Imaging, Technical University of Munich and Helmholtz Center Munich, Neuherberg, Germany.
| |
Collapse
|
34
|
Taskin U, van Dongen KWA. Multi-parameter inversion with the aid of particle velocity field reconstruction. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2020; 147:4032. [PMID: 32611169 DOI: 10.1121/10.0001396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2020] [Accepted: 05/25/2020] [Indexed: 06/11/2023]
Abstract
Multi-parameter inversion for medical ultrasound leads to an improved tissue classification. In general, simultaneous reconstruction of volume density of mass and compressibility would require knowledge of the particle velocity field along with the pressure field. However, in practice the particle velocity field is not measured. Here, the authors propose a method for multi-parameter inversion where the particle velocity field is reconstructed from the measured pressure field. To this end, the measured pressure field is described using outward propagating Hankel functions. For a synthetic setup, it is shown that the reconstructed particle velocity field matches the forward modelled particle velocity field. Next, the reconstructed particle velocity field is used together with the synthetically measured pressure field to reconstruct density and compressibility profiles with the aid of contrast source inversion. Finally, comparing the reconstructed speed of sound profiles obtained via single-parameter versus multi-parameter inversion shows that multi-parameter outperforms single-parameter inversion with respect to accuracy and stability.
Collapse
Affiliation(s)
- Ulas Taskin
- Department of Imaging Physics, Delft University of Technology, Delft, 2628CJ, The Netherlands
| | - Koen W A van Dongen
- Department of Imaging Physics, Delft University of Technology, Delft, 2628CJ, The Netherlands
| |
Collapse
|
35
|
Perez-Liva M, Udías JM, Camacho J, Merčep E, Deán-Ben XL, Razansky D, Herraiz JL. Speed of sound ultrasound transmission tomography image reconstruction based on Bézier curves. ULTRASONICS 2020; 103:106097. [PMID: 32078843 DOI: 10.1016/j.ultras.2020.106097] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2019] [Revised: 01/09/2020] [Accepted: 01/09/2020] [Indexed: 02/07/2023]
Abstract
Speed of Sound (SoS) maps from ultrasound tomography (UST) provide valuable quantitative information for soft tissue characterization and identification of lesions, making this technique interesting for breast cancer detection. However, due to the complexity of the processes that characterize the interaction of ultrasonic waves with matter, classic and fast tomographic algorithms such as back-projection are not suitable. Consequently, the image reconstruction process in UST is generally slow compared to other more conventional medical tomography modalities. With the aim of facilitating the translation of this technique into real clinical practice, several reconstruction algorithms are being proposed to make image reconstruction in UST to be a fast and accurate process. The geometrical acoustic approximation is often used to reconstruct SoS with less computational burden in comparison with full-wave inversion methods. In this work, we propose a simple formulation to perform on-the-flight reconstruction for UST using geometrical acoustics with refraction correction based on quadratic Bézier polynomials. Here we demonstrate that the trajectories created with these polynomials are an accurate approximation to reproduce the refracted acoustic paths connecting the emitter and receiver transducers. The method is faster than typical acquisition times in UST. Thus, it can be considered a step towards real-time reconstructions, which may contribute to its future clinical translation.
Collapse
Affiliation(s)
- Mailyn Perez-Liva
- Grupo de Física Nuclear and IPARCOS, Univ. Complutense de Madrid, CEI Moncloa, Spain; Université de Paris, PARCC, INSERM, F-75015 Paris, France.
| | - José Manuel Udías
- Grupo de Física Nuclear and IPARCOS, Univ. Complutense de Madrid, CEI Moncloa, Spain; Instituto de Investigación Biomédica, Hospital General Universitario Carlos III, Madrid, Spain
| | - Jorge Camacho
- Ultrasound Systems and Technology Group (GSTU), Spanish National Research Council (CSIC), Madrid, Spain
| | - Elena Merčep
- iThera Medical GmbH, Munich, Germany; Faculty of Medicine, Technical University of Munich, Germany
| | - Xosé Luís Deán-Ben
- Faculty of Medicine and Institute of Pharmacology and Toxicology, University of Zurich, Switzerland; Institute for Biomedical Engineering and Department of Information Technology and Electrical Engineering, ETH Zurich, Switzerland
| | - Daniel Razansky
- Faculty of Medicine and Institute of Pharmacology and Toxicology, University of Zurich, Switzerland; Institute for Biomedical Engineering and Department of Information Technology and Electrical Engineering, ETH Zurich, Switzerland
| | - Joaquín L Herraiz
- Grupo de Física Nuclear and IPARCOS, Univ. Complutense de Madrid, CEI Moncloa, Spain; Instituto de Investigación Biomédica, Hospital General Universitario Carlos III, Madrid, Spain
| |
Collapse
|
36
|
Baker C, Sarno D, Eckersley RJ, Zeqiri B. Ring Artifact Correction for Phase-Insensitive Ultrasound Computed Tomography. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2020; 67:513-525. [PMID: 31634829 DOI: 10.1109/tuffc.2019.2948429] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
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.
Collapse
|
37
|
Taskin U, van der Neut J, Gemmeke H, van Dongen KWA. Redatuming of 2-D Wave Fields Measured on an Arbitrary-Shaped Closed Aperture. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2020; 67:173-179. [PMID: 31545717 DOI: 10.1109/tuffc.2019.2942453] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Whole-breast ultrasound scanning systems are used to screen a women's breast for suspicious lesions. Typically, the transducers are located at fixed positions at relatively large distances from the breast to avoid any contact with the breast. Unfortunately, these large distances give rise to large spatial domains to be imaged. These large domains hamper the applicability of imaging by inversion. To reduce the size of the spatial computational domain, we present a 2-D redatuming method based on the Hankel decomposition of the measured field. With this method, the field measured over an arbitrary-shaped closed curve can be redatumed to a new curve enclosing a smaller spatial domain. Additional advantages of the proposed method are that it allows to account for the finite size and orientation of a transducer and that it is robust to noise. The proposed method is successfully validated using the synthetic and measured data, and the results show that the recorded field can be redatumed to any position in the embedding.
Collapse
|
38
|
Korta Martiartu N, Boehm C, Hapla V, Maurer H, Balic IJ, Fichtner A. Optimal experimental design for joint reflection-transmission ultrasound breast imaging: From ray- to wave-based methods. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2019; 146:1252. [PMID: 31472544 DOI: 10.1121/1.5122291] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 07/30/2019] [Indexed: 06/10/2023]
Abstract
Ultrasound computed tomography (USCT) is an emerging modality to image the acoustic properties of the breast tissue for cancer diagnosis. With the need of improving the diagnostic accuracy of USCT, while maintaining the cost low, recent research is mainly focused on improving (1) the reconstruction methods and (2) the acquisition systems. D-optimal sequential experimental design (D-SOED) offers a method to integrate these aspects into a common systematic framework. The transducer configuration is optimized to minimize the uncertainties in the estimated model parameters, and to reduce the time to solution by identifying redundancies in the data. This work presents a formulation to jointly optimize the experiment for transmission and reflection data and, in particular, to estimate the speed of sound and reflectivity of the tissue using either ray-based or wave-based imaging methods. Uncertainties in the parameters can be quantified by extracting properties of the posterior covariance operator, which is analytically computed by linearizing the forward problem with respect to the prior knowledge about parameters. D-SOED is first introduced by an illustrative toy example, and then applied to real data. This shows that the time to solution can be substantially reduced, without altering the final image, by selecting the most informative measurements.
Collapse
Affiliation(s)
- Naiara Korta Martiartu
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Christian Boehm
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Vaclav Hapla
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | - Hansruedi Maurer
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| | | | - Andreas Fichtner
- Department of Earth Sciences, Eidgenössische Technische Hochschule Zürich, Zürich CH-8092, Switzerland
| |
Collapse
|
39
|
Giza OM, Sánchez-Parcerisa D, Sánchez-Tembleque V, Herraiz JL, Camacho J, Avery S, Udías JM. Photoacoustic dose monitoring in clinical high-energy photon beams. Biomed Phys Eng Express 2019. [DOI: 10.1088/2057-1976/ab04ed] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
40
|
Almekkawy M, Chen J, Ellis MD, Haemmerich D, Holmes DR, Linte CA, Panescu D, Pearce J, Prakash P, Zderic V. Therapeutic Systems and Technologies: State-of-the-Art Applications, Opportunities, and Challenges. IEEE Rev Biomed Eng 2019; 13:325-339. [PMID: 30951478 PMCID: PMC7341980 DOI: 10.1109/rbme.2019.2908940] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
In this review, we present current state-of-the-art developments and challenges in the areas of thermal therapy, ultrasound tomography, image-guided therapies, ocular drug delivery, and robotic devices in neurorehabilitation. Additionally, intellectual property and regulatory aspects pertaining to therapeutic systems and technologies are addressed.
Collapse
|
41
|
Gu J, Jing Y. Numerical Modeling of Ultrasound Propagation in Weakly Heterogeneous Media Using a Mixed-Domain Method. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2018; 65:1258-1267. [PMID: 29993378 PMCID: PMC6055067 DOI: 10.1109/tuffc.2018.2828316] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
A mixed-domain method (MDM) is presented in this paper for modeling one-way linear/nonlinear wave propagation in biological tissue with arbitrary heterogeneities, in which sound speed, density, attenuation coefficients, and nonlinear coefficients are all spatial varying functions. The present method is based on solving an integral equation derived from a Westervelt-like equation. One-dimensional problems are first studied to verify the MDM and to reveal its limitations. It is shown that this method is accurate for cases with small variation of sound speed. A 2-D case is further studied with focused ultrasound beams to validate the application of the method in the medical field. Results from the MATLAB toolbox k-Wave are used as the benchmark. Normalized root-mean-square (rms) error estimated at the focus of the transducer is 0.0133 when the coarsest mesh (1/3 of the wavelength) is used in the MDM. Fundamental and second-harmonic fields throughout the considered computational domains are compared and good agreement is observed. Overall, this paper demonstrates that the MDM is a computationally efficient and accurate method when used to model wave propagation in biological tissue with relatively weak heterogeneities.
Collapse
|
42
|
Nguyen LT, Modrak RT. Ultrasonic wavefield inversion and migration in complex heterogeneous structures: 2D numerical imaging and nondestructive testing experiments. ULTRASONICS 2018; 82:357-370. [PMID: 28988141 DOI: 10.1016/j.ultras.2017.09.011] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2017] [Revised: 08/24/2017] [Accepted: 09/13/2017] [Indexed: 06/07/2023]
Abstract
Delaminations, cracks and other defects in engineered structures often lie close to the theoretical resolution limit for ultrasonic waves. While ultrasonic waveform tomography has succeeded in detecting such features, recovery is difficult because it requires computationally expensive high-frequency numerical wave simulations and an accurate understanding of large-scale background variations of the engineered structure. Without such knowledge, small defects may be incorrectly imaged or go undetected altogether. To reduce computational cost and improve detection of small defects, a useful approach is to divide the waveform tomography procedure into two steps: first, a low-frequency model-building step aimed at recovering background structure, and second, a high-frequency imaging step targeting defects. The first is naturally formulated as waveform inversion for wavespeed parameters and the second as time reversal migration for reflectivity. Through synthetic test cases, we show that the two-step workflow appears more promising in most cases than a single-step inversion. In particular, we find that new workflow succeeds in the challenging scenario where the defect lies along preexisting layer interface in a composite bridge deck and in related experiments involving noisy data or inaccurate source parameters.
Collapse
Affiliation(s)
- Luan T Nguyen
- International Geothermal Centre, Bochum University of Applied Sciences, Bochum, Germany.
| | - Ryan T Modrak
- Geophysical Institute, University of Alaska Fairbanks, AK, United States.
| |
Collapse
|
43
|
Bernard S, Monteiller V, Komatitsch D, Lasaygues P. Ultrasonic computed tomography based on full-waveform inversion for bone quantitative imaging. ACTA ACUST UNITED AC 2017; 62:7011-7035. [DOI: 10.1088/1361-6560/aa7e5a] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|