1
|
Rosnitskiy PB, Khokhlova TD, Schade GR, Sapozhnikov OA, Khokhlova VA. Treatment Planning and Aberration Correction Algorithm for HIFU Ablation of Renal Tumors. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2024; 71:341-353. [PMID: 38231825 PMCID: PMC11003458 DOI: 10.1109/tuffc.2024.3355390] [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: 01/19/2024]
Abstract
High-intensity focused ultrasound (HIFU) applications for thermal or mechanical ablation of renal tumors often encounter challenges due to significant beam aberration and refraction caused by oblique beam incidence, inhomogeneous tissue layers, and presence of gas and bones within the beam. These losses can be significantly mitigated through sonication geometry planning, patient positioning, and aberration correction using multielement phased arrays. Here, a sonication planning algorithm is introduced, which uses the simulations to select the optimal transducer position and evaluate the effect of aberrations and acoustic field quality at the target region after aberration correction. Optimization of transducer positioning is implemented using a graphical user interface (GUI) to visualize a segmented 3-D computed tomography (CT)-based acoustic model of the body and to select sonication geometry through a combination of manual and automated approaches. An HIFU array (1.5 MHz, 256 elements) and three renal cell carcinoma (RCC) cases with different tumor locations and patient body habitus were considered. After array positioning, the correction of aberrations was performed using a combination of backpropagation from the focus with an ordinary least squares (OLS) optimization of phases at the array elements. The forward propagation was simulated using a combination of the Rayleigh integral and k-space pseudospectral method (k-Wave toolbox). After correction, simulated HIFU fields showed tight focusing and up to threefold higher maximum pressure within the target region. The addition of OLS optimization to the aberration correction method yielded up to 30% higher maximum pressure compared to the conventional backpropagation and up to 250% higher maximum pressure compared to the ray-tracing method, particularly in strongly distorted cases.
Collapse
|
2
|
Gao P, Sun Y, Zhang G, Li C, Wang L. A transducer positioning method for transcranial focused ultrasound treatment of brain tumors. Front Neurosci 2023; 17:1277906. [PMID: 37904813 PMCID: PMC10613465 DOI: 10.3389/fnins.2023.1277906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2023] [Accepted: 09/28/2023] [Indexed: 11/01/2023] Open
Abstract
Purpose As a non-invasive method for brain diseases, transcranial focused ultrasound (tFUS) offers higher spatial precision and regulation depth. Due to the altered path and intensity of sonication penetrating the skull, the focus and intensity in the skull are difficult to determine, making the use of ultrasound therapy for cancer treatment experimental and not widely available. The deficiency can be effectively addressed by numerical simulation methods, which enable the optimization of sonication modulation parameters and the determination of precise transducer positioning. Methods A 3D skull model was established using binarized brain CT images. The selection of the transducer matrix was performed using the radius positioning (RP) method after identifying the intracranial target region. Simulations were performed, encompassing acoustic pressure (AP), acoustic field, and temperature field, in order to provide compelling evidence of the safety of tFUS in sonication-induced thermal effects. Results It was found that the angle of sonication path to the coronal plane obtained at all precision and frequency models did not exceed 10° and 15° to the transverse plane. The results of thermal effects illustrated that the peak temperatures of tFUS were 43.73°C, which did not reach the point of tissue degeneration. Once positioned, tFUS effectively delivers a Full Width at Half Maximum (FWHM) stimulation that targets tumors with diameters of up to 3.72 mm in a one-off. The original precision model showed an attenuation of 24.47 ± 6.13 mm in length and 2.40 ± 1.42 mm in width for the FWHM of sonication after penetrating the skull. Conclusion The vector angles of the sonication path in each direction were determined based on the transducer positioning results. It has been suggested that when time is limited for precise transducer positioning, fixing the transducer on the horizontal surface of the target region can also yield positive results for stimulation. This framework used a new transducer localization method to offer a reliable basis for further research and offered new methods for the use of tFUS in brain tumor-related research.
Collapse
Affiliation(s)
- Penghao Gao
- Artificial Intelligence Laboratory, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
| | - Yue Sun
- Department of Biomedical Engineering, Shenyang University of Technology, Shenyang, Liaoning, China
| | - Gongsen Zhang
- Artificial Intelligence Laboratory, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
| | - Chunsheng Li
- Department of Biomedical Engineering, Shenyang University of Technology, Shenyang, Liaoning, China
| | - Linlin Wang
- Artificial Intelligence Laboratory, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
- Department of Radiation Oncology, Shandong Cancer Hospital and Institute, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan, Shandong, China
| |
Collapse
|
3
|
Cao E, Greillier P, Loyet R, Chavrier F, Robert J, Bessière F, Dillenseger JL, Lafon C. Development of a Numerical Model of High-Intensity Focused Ultrasound Treatment in Mobile and Elastic Organs: Application to a Beating Heart. ULTRASOUND IN MEDICINE & BIOLOGY 2022; 48:1215-1228. [PMID: 35430101 DOI: 10.1016/j.ultrasmedbio.2022.02.017] [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/29/2021] [Revised: 02/14/2022] [Accepted: 02/16/2022] [Indexed: 06/14/2023]
Abstract
High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.
Collapse
Affiliation(s)
- Elodie Cao
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France..
| | - Paul Greillier
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France
| | - Raphaël Loyet
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France
| | - Françoise Chavrier
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France
| | - Jade Robert
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France
| | - Francis Bessière
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France.; Hospices Civils de Lyon, Hôpital Cardiovasculaire Louis Pradel, Lyon, France
| | | | - Cyril Lafon
- LabTAU, INSERM, Centre Léon Bérard, Université Lyon 1, Univ Lyon, F-69003, LYON, France
| |
Collapse
|
4
|
Papež P, Praprotnik M. Dissipative Particle Dynamics Simulation of Ultrasound Propagation through Liquid Water. J Chem Theory Comput 2022; 18:1227-1240. [PMID: 35001631 PMCID: PMC8830050 DOI: 10.1021/acs.jctc.1c01020] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Ultrasound is widely
used as a noninvasive method in therapeutic
and diagnostic applications. These can be further optimized by computational
approaches, as they allow for controlled testing and rational optimization
of the ultrasound parameters, such as frequency and amplitude. Usually,
continuum numerical methods are used to simulate ultrasound propagating
through different tissue types. In contrast, ultrasound simulations
using particle description are less common, as the implementation
is challenging. In this work, a dissipative particle dynamics model
is used to perform ultrasound simulations in liquid water. The effects
of frequency and thermostat parameters are studied and discussed.
We show that frequency and thermostat parameters affect not only the
attenuation but also the computed speed of sound. The present study
paves the way for development and optimization of a virtual ultrasound
machine for large-scale biomolecular simulations.
Collapse
Affiliation(s)
- Petra Papež
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, Ljubljana, SI-1001, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana, SI-1000, Slovenia
| | - Matej Praprotnik
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, Ljubljana, SI-1001, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana, SI-1000, Slovenia
| |
Collapse
|
5
|
Zhang J, Lay RJ, Roberts SK, Chauhan S. Towards real-time finite-strain anisotropic thermo-visco-elastodynamic analysis of soft tissues for thermal ablative therapy. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2021; 198:105789. [PMID: 33069033 DOI: 10.1016/j.cmpb.2020.105789] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2020] [Accepted: 10/05/2020] [Indexed: 06/11/2023]
Abstract
BACKGROUND AND OBJECTIVES Accurate and efficient prediction of soft tissue temperatures is essential to computer-assisted treatment systems for thermal ablation. It can be used to predict tissue temperatures and ablation volumes for personalised treatment planning and image-guided intervention. Numerically, it requires full nonlinear modelling of the coupled computational bioheat transfer and biomechanics, and efficient solution procedures; however, existing studies considered the bioheat analysis alone or the coupled linear analysis, without the fully coupled nonlinear analysis. METHODS We present a coupled thermo-visco-hyperelastic finite element algorithm, based on finite-strain thermoelasticity and total Lagrangian explicit dynamics. It considers the coupled nonlinear analysis of (i) bioheat transfer under soft tissue deformations and (ii) soft tissue deformations due to thermal expansion/shrinkage. The presented method accounts for anisotropic, finite-strain, temperature-dependent, thermal, and viscoelastic behaviours of soft tissues, and it is implemented using GPU acceleration for real-time computation. RESULTS The presented method can achieve thermo-visco-elastodynamic analysis of anisotropic soft tissues undergoing large deformations with high computational speeds in tetrahedral and hexahedral finite element meshes for surgical simulation of thermal ablation. We also demonstrate the translational benefits of the presented method for clinical applications using a simulation of thermal ablation in the liver. CONCLUSION The key advantage of the presented method is that it enables full nonlinear modelling of the anisotropic, finite-strain, temperature-dependent, thermal, and viscoelastic behaviours of soft tissues, instead of linear elastic, linear viscoelastic, and thermal-only modelling in the existing methods. It also provides high computational speeds for computer-assisted treatment systems towards enabling the operator to simulate thermal ablation accurately and visualise tissue temperatures and ablation zones immediately.
Collapse
Affiliation(s)
- Jinao Zhang
- Department of Mechanical and Aerospace Engineering, Monash University, Clayton, Victoria, Australia.
| | - Remi Jacob Lay
- Department of Mechanical and Aerospace Engineering, Monash University, Clayton, Victoria, Australia
| | - Stuart K Roberts
- Department of Gastroenterology, The Alfred Hospital, Melbourne, Victoria, Australia
| | - Sunita Chauhan
- Department of Mechanical and Aerospace Engineering, Monash University, Clayton, Victoria, Australia.
| |
Collapse
|
6
|
Zhang J, Chauhan S. Fast computation of soft tissue thermal response under deformation based on fast explicit dynamics finite element algorithm for surgical simulation. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 187:105244. [PMID: 31805458 DOI: 10.1016/j.cmpb.2019.105244] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Revised: 11/03/2019] [Accepted: 11/26/2019] [Indexed: 06/10/2023]
Abstract
BACKGROUND AND OBJECTIVES During thermal heating surgical procedures such as electrosurgery, thermal ablative treatment and hyperthermia, soft tissue deformation due to surgical tool-tissue interaction and patient movement can affect the distribution of thermal energy induced. Soft tissue temperature must be obtained from the deformed tissue for precise delivery of thermal energy. However, the classical Pennes bio-heat transfer model can handle only the static non-moving state of tissue. In addition, in order to enable a surgeon to visualise the simulated results immediately, the solution procedure must be suitable for real-time thermal applications. METHODS This paper presents a formulation of bio-heat transfer under the effect of soft tissue deformation for fast or near real-time tissue temperature prediction, based on fast explicit dynamics finite element algorithm (FED-FEM) for transient heat transfer. The proposed thermal analysis under deformation is achieved by transformation of the unknown deformed tissue state to the known initial static state via a mapping function. The appropriateness and effectiveness of the proposed formulation are evaluated on a realistic virtual human liver model with blood vessels to demonstrate a clinically relevant scenario of thermal ablation of hepatic cancer. RESULTS For numerical accuracy, the proposed formulation can achieve a typical 10-3 level of normalised relative error at nodes and between 10-4 and 10-5 level of total errors for the simulation, by comparing solutions against the commercial finite element analysis package. For computation time, the proposed formulation under tissue deformation with anisotropic temperature-dependent properties consumes 2.518 × 10-4 ms for one element thermal loads computation, compared to 2.237 × 10-4 ms for the formulation without deformation which is 0.89 times of the former. Comparisons with three other formulations for isotropic and temperature-independent properties are also presented. CONCLUSIONS Compared to conventional methods focusing on numerical accuracy, convergence and stability, the proposed formulation focuses on computational performance for fast tissue thermal analysis. Compared to the classical Pennes model that handles only the static state of tissue, the proposed formulation can achieve fast thermal analysis on deformed states of tissue and can be applied in addition to tissue deformable models for non-linear heating analysis at even large deformation of soft tissue, leading to great translational potential in dynamic tissue temperature analysis and thermal dosimetry computation for computer-integrated medical education and personalised treatment.
Collapse
Affiliation(s)
- Jinao Zhang
- Department of Mechanical and Aerospace Engineering, Monash University, Wellington Road, Clayton, VIC 3800, Australia.
| | - Sunita Chauhan
- Department of Mechanical and Aerospace Engineering, Monash University, Wellington Road, Clayton, VIC 3800, Australia
| |
Collapse
|
7
|
Neural network methodology for real-time modelling of bio-heat transfer during thermo-therapeutic applications. Artif Intell Med 2019; 101:101728. [DOI: 10.1016/j.artmed.2019.101728] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 07/30/2019] [Accepted: 09/26/2019] [Indexed: 12/26/2022]
|
8
|
Cheema MN, Nazir A, Sheng B, Li P, Qin J, Kim J, Feng DD. Image-Aligned Dynamic Liver Reconstruction Using Intra-Operative Field of Views for Minimal Invasive Surgery. IEEE Trans Biomed Eng 2018; 66:2163-2173. [PMID: 30507524 DOI: 10.1109/tbme.2018.2884319] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
During hepatic minimal invasive surgery (MIS), 3-D reconstruction of a liver surface by interpreting the geometry of its soft tissues is achieving attractions. One of the major issues to be addressed in MIS is liver deformation. Moreover, it severely inhibits free sight and dexterity of tissue manipulation, which causes its intra-operative morphology and soft tissue motion altered as compared to its pre-operative shape. While many applications focus on 3-D reconstruction of rigid or semi-rigid scenes, the techniques applied in hepatic MIS must be able to cope with a dynamic and deformable environment. We propose an efficient technique for liver surface reconstruction based on the structure from motion to handle liver deformation. The reconstructed liver will assist surgeons to visualize liver surface more efficiently with better depth perception. We use the intra-operative field of views to generate 3-D template mesh from a dense keypoint cloud. We estimate liver deformation by finding best correspondence between 3-D templates and reconstruct a liver image to calculate translation and rotational motions. Our technique then finely tunes deformed surface by adding smoothness using shading cues. Up till now, this technique is not used for solving the human liver deformation problem. Our approach is tested and validated with synthetic as well as real in vivo data, which reveal that the reconstruction accuracy can be enhanced using our approach even in challenging laparoscopic environments.
Collapse
|
9
|
Bousselham A, Bouattane O, Youssfi M, Raihani A. 3D brain tumor localization and parameter estimation using thermographic approach on GPU. J Therm Biol 2018; 71:52-61. [DOI: 10.1016/j.jtherbio.2017.10.014] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2017] [Revised: 10/14/2017] [Accepted: 10/19/2017] [Indexed: 10/18/2022]
|
10
|
Schwenke M, Strehlow J, Demedts D, Haase S, Barrios Romero D, Rothlübbers S, von Dresky C, Zidowitz S, Georgii J, Mihcin S, Bezzi M, Tanner C, Sat G, Levy Y, Jenne J, Günther M, Melzer A, Preusser T. A focused ultrasound treatment system for moving targets (part I): generic system design and in-silico first-stage evaluation. J Ther Ultrasound 2017; 5:20. [PMID: 28748092 PMCID: PMC5523151 DOI: 10.1186/s40349-017-0098-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2017] [Accepted: 05/11/2017] [Indexed: 12/25/2022] Open
Abstract
Background Focused ultrasound (FUS) is entering clinical routine as a treatment option. Currently, no clinically available FUS treatment system features automated respiratory motion compensation. The required quality standards make developing such a system challenging. Methods A novel FUS treatment system with motion compensation is described, developed with the goal of clinical use. The system comprises a clinically available MR device and FUS transducer system. The controller is very generic and could use any suitable MR or FUS device. MR image sequences (echo planar imaging) are acquired for both motion observation and thermometry. Based on anatomical feature tracking, motion predictions are estimated to compensate for processing delays. FUS control parameters are computed repeatedly and sent to the hardware to steer the focus to the (estimated) target position. All involved calculations produce individually known errors, yet their impact on therapy outcome is unclear. This is solved by defining an intuitive quality measure that compares the achieved temperature to the static scenario, resulting in an overall efficiency with respect to temperature rise. To allow for extensive testing of the system over wide ranges of parameters and algorithmic choices, we replace the actual MR and FUS devices by a virtual system. It emulates the hardware and, using numerical simulations of FUS during motion, predicts the local temperature rise in the tissue resulting from the controls it receives. Results With a clinically available monitoring image rate of 6.67 Hz and 20 FUS control updates per second, normal respiratory motion is estimated to be compensable with an estimated efficiency of 80%. This reduces to about 70% for motion scaled by 1.5. Extensive testing (6347 simulated sonications) over wide ranges of parameters shows that the main source of error is the temporal motion prediction. A history-based motion prediction method performs better than a simple linear extrapolator. Conclusions The estimated efficiency of the new treatment system is already suited for clinical applications. The simulation-based in-silico testing as a first-stage validation reduces the efforts of real-world testing. Due to the extensible modular design, the described approach might lead to faster translations from research to clinical practice.
Collapse
Affiliation(s)
- Michael Schwenke
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Jan Strehlow
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Daniel Demedts
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Sabrina Haase
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Diego Barrios Romero
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Sven Rothlübbers
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany.,Mediri, Heidelberg, Germany
| | - Caroline von Dresky
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Stephan Zidowitz
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Joachim Georgii
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany
| | - Senay Mihcin
- Institute for Medical Science and Technology, Dundee, Scotland
| | - Mario Bezzi
- Universita Degli Studi Di Roma La Sapienza, Rome, Italy
| | - Christine Tanner
- Computer Vision Laboratory, Eidgenössische Technische Hochschule, Zurich, Switzerland
| | - Giora Sat
- GE Medical Systems Israel, Haifa, Israel
| | | | - Jürgen Jenne
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany.,Mediri, Heidelberg, Germany
| | - Matthias Günther
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany.,Mediri, Heidelberg, Germany
| | - Andreas Melzer
- Institute for Medical Science and Technology, Dundee, Scotland.,Innovation Center Computer Assisted Surgery, Leipzig, Germany
| | - Tobias Preusser
- Fraunhofer Institute for Medical Image Computing MEVIS, Am Fallturm 1, Bremen, 28359 Germany.,Jacobs University, Bremen, Germany
| |
Collapse
|