Modelling of the acoustic field of a multi-element HIFU array scattered by human ribs

Rib detection using pitch-catch ultrasound and classification algorithms for a novel ultrasound therapy device

BACKGROUND

Noninvasive ultrasound (US) has been used therapeutically for decades, with applications in tissue ablation, lithotripsy, and physical therapy. There is increasing evidence that low intensity US stimulation of organs can alter physiological and clinical outcomes for treatment of health disorders including rheumatoid arthritis and diabetes. One major translational challenge is designing portable and reliable US devices that can be used by patients in their homes, with automated features to detect rib location and aid in efficient transmission of energy to organs of interest. This feasibility study aimed to assess efficacy in rib bone detection without conventional imaging, using a single channel US pitch-catch technique integrated into an US therapy device to detect pulsed US reflections from ribs.

METHODS

In 20 healthy volunteers, the location of the ribs and spleen were identified using a diagnostic US imaging system. Reflected ultrasound signals were recorded at five positions over the spleen and adjacent ribs using the therapy device. Signals were classified as between ribs (intercostal), partially over a rib, or fully over a rib using four models: threshold-based time domain classification, threshold-based frequency domain classification, logistic regression, and support vector machine (SVM).

RESULTS

SVM performed best overall on the All Participants cohort with accuracy up to 96.25%. All models' accuracies were improved by separating participants into two cohorts based on Body Mass Index (BMI) and re-fitting each model. After separation into Low BMI and High BMI cohorts, a simple time-thresholding approach achieved accuracies up to 100% and 93.75%, respectively.

CONCLUSION

These results demonstrate that US reflection signal classification can accurately provide low complexity, real-time automated onboard rib detection and user feedback to advance at-home therapeutic US delivery.

Accelerating frequency-domain numerical methods for weakly nonlinear focused ultrasound using nested meshes

The numerical simulation of weakly nonlinear ultrasound is important in treatment planning for focused ultrasound (FUS) therapies. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough to resolve the highest harmonic present in the problem, leading to a very large number of degrees of freedom. This paper proposes a more efficient strategy in which each harmonic is approximated on a separate mesh, the size of which is proportional to the wavelength of the harmonic. The increase in resolution required to resolve a smaller wavelength is balanced by a reduction in the domain size. This nested meshing is feasible owing to the increasingly localised nature of higher harmonics near the focus. Numerical experiments are performed for FUS transducers in homogeneous media to determine the size of the meshes required to accurately represent the harmonics. In particular, a fast volume potential approach is proposed and employed to perform convergence experiments as the computation domain size is modified. This approach allows each harmonic to be computed via the evaluation of an integral over the domain. Discretising this integral using the midpoint rule allows the computations to be performed rapidly with the FFT. It is shown that at least an order of magnitude reduction in memory consumption and computation time can be achieved with nested meshing. Finally, it is demonstrated how to generalise this approach to inhomogeneous propagation domains.

Calculating the Effect of Ribs on the Focus Quality of a Therapeutic Spherical Random Phased Array

The overlaying rib cage is a major hindrance in treating liver tumors with high intensity focused ultrasound (HIFU). The problems caused are overheating of the ribs due to its high ultrasonic absorption capability and degradation of the ultrasound intensity distribution in the target plane. In this work, a correction method based on binarized apodization and geometric ray tracing approach was employed to avoid heating the ribs. A detailed calculation of the intensity distribution in the focus plane was undertaken to quantify and avoid the effect on HIFU beam generated by a 1-MHz 256-element random phased array after the ultrasonic beam passes through the rib cage. Focusing through the ribs was simulated for 18 different idealized ribs-array configurations and 10 anatomically correct ribs-array configurations, to show the effect of width of the ribs, intercostal spacing and the relative position of ribs and array on the quality of focus, and to identify the positions that are more effective for HIFU applications in the presence of ribs. Acoustic simulations showed that for a single focus without beam steering and for the same total acoustic power, the peak intensity at the target varies from a minimum of 211 W/cm² to a maximum of 293 W/cm² for a nominal acoustic input power of 15 W, whereas the side lobe level varies from 0.07 I_peak to 0.28 I_peak and the separation between the main lobe and side lobes varies from 2.5 mm to 6.3 mm, depending on the relative positioning of the array and ribs and the beam alignment. An increase in the side lobe level was observed by increasing the distance between the array and the ribs. The parameters of focus splitting and the deterioration of focus quality caused by the ultrasonic propagation through the ribs were quantified in various possible different clinical scenarios. In addition to idealized rib topology, anatomical realistic ribs were used to determine the focus quality of the HIFU beam when the beam is steered both in axial and transverse directions and when the transducer is positioned at different depths from the rib cage.

Robot-assisted ultrasound navigation platform for 3D HIFU treatment planning: Initial evaluation for conformal interstitial ablation

Interstitial Ultrasound-guided High Intensity Focused Ultrasound (USgHIFU) therapy has the potential to deliver ablative treatments which conform to the target tumor. In this study, a robot-assisted US-navigation platform has been developed for 3D US guidance and planning of conformal HIFU ablations. The platform was used to evaluate a conformal therapeutic strategy associated with an interstitial dual-mode USgHIFU catheter prototype (64 elements linear-array, measured central frequency f = 6.5 MHz), developed for the treatment of HepatoCellular Carcinoma (HCC). The platform included a 3D navigation environment communicating in real-time with an open research dual-mode US scanner/HIFU generator and a robotic arm, on which the USgHIFU catheter was mounted. 3D US-navigation was evaluated in vitro for guiding and planning conformal HIFU ablations using a tumor-mimic model in porcine liver. Tumor-mimic volumes were then used as targets for evaluating conformal HIFU treatment planning in simulation. Height tumor-mimics (ovoid- or disc-shaped, sizes: 3-29 cm³) were created and visualized in liver using interstitial 2D US imaging. Robot-assisted spatial manipulation of these images and real-time 3D navigation allowed reconstructions of 3D B-mode US images for accurate tumor-mimic volume estimation (relative error: 4 ± 5%). Sectorial and full-revolution HIFU scanning (angular sectors: 88-360°) could both result in conformal ablations of the tumor volumes, as soon as their radii remained ≤ 24 mm. The presented US navigation-guided HIFU procedure demonstrated advantages for developing conformal interstitial therapies in standard operative rooms. Moreover, the modularity of the developed platform makes it potentially useful for developing other HIFU approaches.

Patient-Specific 3-Dimensional Model for High-Intensity Focused Ultrasound Treatment Through the Rib Cage: A Preliminary Study

OBJECTIVES

The purpose of this study was to develop a patient-specific 3-dimensional model for high-intensity focused ultrasound (HIFU) treatment through the rib cage using patient data.

METHODS

Experimental testing to derive parameters used in defining the amount of energy and alteration needed in treatment protocols for upper abdominal disorders under the rib cage was performed. Reconstructed rib cage models based on patient data, tissue-mimicking material phantoms, and magnetic resonance imaging-guided HIFU using a multielement phased array transducer were used in the experiments. Changes in the focal temperature, acoustic power, and acoustic pressure distribution were investigated with and without the presence of the rib cage model. An ExAblate system (InSightec Ltd, Tirat Carmel, Israel) was used to sonicate phantoms by varying the target phantom or rib cage model location.

RESULTS

The effect of the rib cage on the acoustic pressure distribution and acoustic power was closely related to the anatomic structures of the ribs. Thermometry revealed that heating at the focus could be controlled by changing either the power or duration of HIFU application to improve the focal temperature change. The focal temperature change was found to be related to the distance between the rib cage model and focus and the shadow area on the transducer elements covered by the rib cage model in the beam path.

CONCLUSIONS

Experimental results suggest that the rib cage model is a valuable and useful tool that can provide realistic human anatomic structures and properties for evaluating the effects of the rib cage on ultrasound propagation.

The Optimization of Transcostal Phased Array Refocusing Using the Semidefinite Relaxation Method

Tumors in organs partially obscured by the rib cage represent a challenge for high-intensity focused ultrasound (HIFU) therapy. The ribs distort the HIFU beams in a manner that reduces the focusing gain at the target, which could result in treatment-limiting collateral damage. In fact, skin burns are a common complication during the ablation of hepatic tumors. This problem can be addressed by employing optimal refocusing algorithms that are designed to achieve a specified focusing gain at the target while controlling the exposure to the ribs in the path of the HIFU beam. However, previously proposed optimal refocusing algorithms did not allow for the controlled transmission through the ribs. In this article, we introduce a new approach for refocusing that can more efficiently steer power toward the target while limiting the power deposition on the ribs. The approach utilizes the semidefinite relaxation (SDR) technique to approximate the original (nonconvex) optimization problem. An important advantage of the SDR-based method over previously proposed optimization methods is the control of the side lobes in the focal plane. The method also allows for specifying an acceptable level of exposure to the ribs. Simulation results using a 1-MHz spherical concave phased array focused on an inhomogeneous medium are presented to demonstrate the performance of the SDR refocusing approach. A finite-difference time-domain propagation model was used to model the propagation in the inhomogeneous tissues, including the ribs. Temperature simulations based on the inhomogeneous transient bioheat transfer equation (tBHTE) demonstrate the significance of the improvements in the focusing gain when using the limited power deposition (LPD) method. The results also demonstrate that the LPD method yields well-behaved array excitation vectors, realizable by currently existing drivers.

Determination of Acoustic Cavitation Probabilities and Thresholds Using a Single Focusing Transducer to Induce and Detect Acoustic Cavitation Events: I. Method and Terminology

A method to determine acoustic cavitation probabilities in tissue-mimicking materials (TMMs) is described that uses a high-intensity focused ultrasound (HIFU) transducer for both inducing and detecting the acoustic cavitation events. The method was evaluated by studying acoustic cavitation probabilities in agar-based TMMs with and without scatterers and for different sonication modes like continuous wave, single pulses (microseconds to milliseconds) and repeated burst signals. Acoustic cavitation thresholds (defined here as the peak rarefactional in situ pressure at which the acoustic cavitation probability reaches 50%) at a frequency of 1.06 MHz were observed between 1.1 MPa (for 1 s of continuous wave sonication) and 4.6 MPa (for 1 s of a repeated burst signal with 25-cycle burst length and 10-ms burst period) in a 3% (by weight) agar phantom without scatterers. The method and its evaluation are described, and general terminology useful for standardizing the description of insonation conditions and comparing results is provided. In the accompanying second part, the presented method is used to systematically study the acoustic cavitation thresholds in the same material for a range of sonication modes.

Blocked Elements in 1-D and 2-D Arrays-Part I: Detection and Basic Compensation on Simulated and In Vivo Targets

During a transcostal ultrasound scan, ribs and other highly attenuating and/or reflective tissue structures can block parts of the array. Blocked elements tend to limit the acoustic window and impede visualization of structures of interest. Here, we demonstrate a method to detect blocked elements and we measure the loss of image quality they introduce in simulation and in vivo. We utilize a fullwave simulation tool and a clinical ultrasound scanner to obtain element signals from fully sampled matrix arrays during simulated and in vivo transcostal liver scans, respectively. The elements that were blocked by a rib showed lower average signal amplitude and lower average nearest-neighbor cross correlation than the elements in the remainder of the 2-D aperture. The growing receive-aperture B-mode images created from the element data indicate that the signals on blocked elements are dominated by noise and that turning them OFF has a potential to improve visibility of liver vasculature. Adding blocked elements to the growing receive apertures for five in vivo transcostal acquisitions resulted in average decrease in vessel contrast and contrast to noise ratio of 19% and 10%, respectively.

Feasibility of an electromagnetic compatibility method for MRgFUS using a wire mesh screen

This study evaluated an electromagnetic compatibility method for high-intensity focused ultrasound (HIFU) and magnetic resonance (MR) imaging in an MR-guided focused ultrasound surgery using a conductive wire mesh screen. This screen has a good ultrasound transmission and shielding effectiveness. A hybrid acoustic simulation method was developed to analyze the effects of mesh parameters and the HIFU working frequency on the acoustic field. Experiments were performed to measure both acoustic pressure profile and radiated electromagnetic noise. With the proposed mesh screen, the electromagnetic radiation emission was reduced by 14dB at 128MHz while the acoustic focal intensity was reduced by less than 11% using one screen. This shielding method is easy to implement and requires no additional phase correction method. This method also improves the quality of MR images.

Simulating Focused Ultrasound Transducers Using Discrete Sources on Regular Cartesian Grids

Accurately representing the behavior of acoustic sources is an important part of ultrasound simulation. This is particularly challenging in ultrasound therapy where multielement arrays are often used. Typically, sources are defined as a boundary condition over a 2-D plane within the computational model. However, this approach can become difficult to apply to arrays with multiple elements distributed over a nonplanar surface. In this paper, a grid-based discrete source model for single- and multielement bowl-shaped transducers is developed to model the source geometry explicitly within a regular Cartesian grid. For each element, the source model is defined as a symmetric, simply connected surface with a single grid point thickness. Simulations using the source model with the open-source k-Wave toolbox are validated using the Rayleigh integral, O'Neil solution, and experimental measurements of a focused bowl transducer under both quasi-continuous wave and pulsed excitations. Close agreement is shown between the discrete bowl model and the axial pressure predicted by the O'Neil solution for a uniform curved radiator, even at very low grid resolutions. Excellent agreement is also shown between the discrete bowl model and experimental measurements. To accurately reproduce the near-field pressure measured experimentally, it is necessary to derive the drive signal at each grid point of the bowl model directly using holography. However, good agreement is also obtained in the focal region using uniformly radiating monopole sources distributed over the bowl surface. This allows the response of multielement transducers to be modeled, even where measurement of an input plane is not possible.

A fast boundary element method for the scattering analysis of high-intensity focused ultrasound

High-intensity focused ultrasound (HIFU) techniques are promising modalities for the non-invasive treatment of cancer. For HIFU therapies of, e.g., liver cancer, one of the main challenges is the accurate focusing of the acoustic field inside a ribcage. Computational methods can play an important role in the patient-specific planning of these transcostal HIFU treatments. This requires the accurate modeling of acoustic scattering at ribcages. The use of a boundary element method (BEM) is an effective approach for this purpose because only the boundaries of the ribs have to be discretized instead of the standard approach to model the entire volume around the ribcage. This paper combines fast algorithms that improve the efficiency of BEM specifically for the high-frequency range necessary for transcostal HIFU applications. That is, a Galerkin discretized Burton-Miller formulation is used in combination with preconditioning and matrix compression techniques. In particular, quick convergence is achieved with the operator preconditioner that has been designed with on-surface radiation conditions for the high-frequency approximation of the Neumann-to-Dirichlet map. Realistic computations of acoustic scattering at 1 MHz on a human ribcage model demonstrate the effectiveness of this dedicated BEM algorithm for HIFU scattering analysis.

Ultrasound-guided therapeutic focused ultrasound: current status and future directions

This paper reviews ultrasound imaging methods for the guidance of therapeutic focused ultrasound (USgFUS), with emphasis on real-time preclinical methods. Guidance is interpreted in the broadest sense to include pretreatment planning, siting of the FUS focus, real-time monitoring of FUS-tissue interactions, and real-time control of exposure and damage assessment. The paper begins with an overview and brief historical background of the early methods used for monitoring FUS-tissue interactions. Current imaging methods are described, and discussed in terms of sensitivity and specificity of the localisation of the FUS effects in both therapeutic and sub-therapeutic modes. Thermal and non-thermal effects are considered. These include cavitation-enhanced heating, tissue water boiling and cavitation. Where appropriate, USgFUS methods are compared with similar methods implemented using other guidance modalities, e.g. magnetic resonance imaging. Conclusions are drawn regarding the clinical potential of the various guidance methods, and the feasibility and current status of real-time implementation.

A comparison of methods for focusing the field of a HIFU array transducer through human ribs

A forward model, which predicts the scattering by human ribs of a multi-element high-intensity focused ultrasound transducer, was used to investigate the efficacy of a range of focusing approaches described in the literature. This forward model is based on the boundary element method and was described by Gélat et al (2011 Phys. Med. Biol. 56 5553-81; 2012 Phys. Med. Biol. 57 8471-97). The model has since been improved and features a complex surface impedance condition at the surface of the ribs. The inverse problem of focusing through the ribs was implemented on six transducer array-rib topologies and five methods of focusing were investigated, including spherical focusing, binarized apodization based on geometric ray tracing, phase conjugation and the decomposition of the time-reversal operator method. The excitation frequency was 1 MHz and the array was of spherical-section type. Both human and idealized rib topologies were considered. The merit of each method of focusing was examined. It was concluded that the constrained optimization approach offers greater potential than the other focusing methods in terms of maximizing the ratio of acoustic pressure magnitudes at the focus to those on the surface of the ribs whilst taking full advantage of the dynamic range of the phased array.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39-44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from 'model' to 'clinic'. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

The role of acoustic nonlinearity in tissue heating behind a rib cage using a high-intensity focused ultrasound phased array

The goal of this study was to investigate theoretically the effects of nonlinear propagation in a high-intensity focused ultrasound (HIFU) field produced by a therapeutic phased array and the resultant heating of tissue behind a rib cage. Three configurations of focusing were simulated: in water, in water with ribs in the beam path and in water with ribs backed by a layer of soft tissue. The Westervelt equation was used to model the nonlinear HIFU field, and a 1 MHz phased array consisting of 254 circular elements was used as a boundary condition to the model. The temperature rise in tissue was modelled using the bioheat equation, and thermally necrosed volumes were calculated using the thermal dose formulation. The shapes of lesions predicted by the modelling were compared with those previously obtained in in vitro experiments at low-power sonications. Intensity levels at the face of the array elements that corresponded to the formation of high-amplitude shock fronts in the focal region were determined as 10 W cm(-2) in the free field in water and 40 W cm(-2) in the presence of ribs. It was shown that exposures with shocks provided a substantial increase in tissue heating, and its better spatial localization in the main focal region only. The relative effects of overheating ribs and splitting of the focus due to the periodic structure of the ribs were therefore reduced. These results suggest that utilizing nonlinear propagation and shock formation effects can be beneficial for inducing confined HIFU lesions when irradiating through obstructions such as ribs. Design of compact therapeutic arrays to provide maximum power outputs with lower intensity levels at the elements is necessary to achieve shock wave regimes for clinically relevant sonication depths in tissue.

The optimization of acoustic fields for ablative therapies of tumours in the upper abdomen

High intensity focused ultrasound (HIFU) enables highly localized, non-invasive tissue ablation and its efficacy has been demonstrated in the treatment of a range of cancers, including those of the kidney, prostate and breast. HIFU offers the ability to treat deep-seated tumours locally, and potentially bears fewer side effects than more invasive treatment modalities such as resection, chemotherapy and ionizing radiation. There remains however a number of significant challenges which currently hinder its widespread clinical application. One of these challenges is the need to transmit sufficient energy through the ribcage to ablate tissue at the required foci whilst minimizing the formation of side lobes and sparing healthy tissue. Ribs both absorb and reflect ultrasound strongly. This sometimes results in overheating of bone and overlying tissue during treatment, leading to skin burns. Successful treatment of a patient with tumours in the upper abdomen therefore requires a thorough understanding of the way acoustic and thermal energy is deposited. Previously, a boundary element approach based on a Generalized Minimal Residual (GMRES) implementation of the Burton-Miller formulation was developed to predict the field of a multi-element HIFU array scattered by human ribs, the topology of which was obtained from CT scan data (Gélat et al 2011 Phys. Med. Biol. 56 5553-81). The present paper describes the reformulation of the boundary element equations as a least-squares minimization problem with nonlinear constraints. The methodology has subsequently been tested at an excitation frequency of 1 MHz on a spherical multi-element array in the presence of ribs. A single array-rib geometry was investigated on which a 50% reduction in the maximum acoustic pressure magnitude on the surface of the ribs was achieved with only a 4% reduction in the peak focal pressure compared to the spherical focusing case. This method was then compared with a binarized apodization approach based on ray tracing and against the decomposition of the time-reversal operator (DORT). In both cases, the constrained optimization provided a superior ratio of focal peak pressure to maximum pressure magnitude on the surface of the ribs.

Theoretical study for safe and efficient energy transfer to deeply implanted devices using ultrasound

The goal of this paper is to prove that a safe and efficient energy transfer is possible between an external transducer located on the patient's skin and a device deeply implanted in the abdomen. An ultrasound propagation model based on the Rayleigh-Sommerfeld diffraction integral is coupled with the data from the Visible Human Project to account for the geometry of the organs in the body. The model is able to predict the amount of acoustic power received by the device for different acoustic paths. The acoustic model is validated by comparison with measurements in water and in heterogeneous liquid phantoms. Care is taken to minimize adverse bioeffects-mainly temperature rise and cavitation in tissues. Simulations based on the bio-heat transfer equation are performed to check that thermal effects are indeed small.