Adaptive projection method applied to three-dimensional ultrasonic focusing and steering through the 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.

An endoluminal cylindrical sectored-ring ultrasound phased-array applicator for minimally-invasive therapeutic ultrasound

BACKGROUND

The size of catheter-based ultrasound devices for delivering ultrasound energy to deep-seated tumors is constrained by the access pathway which limits their therapeutic capabilities.

PURPOSE

To devise and investigate a deployable applicator suitable for minimally-invasive delivery of therapeutic ultrasound, consisting of a 2D cylindrical sectored-ring ultrasound phased array, integrated within an expandable paraboloid-shaped balloon-based reflector. The balloon can be collapsed for compact delivery and expanded close to the target position to mimic a larger-diameter concentric-ring sector-vortex array for enhanced dynamic control of focal depth and volume.

METHODS

Acoustic and biothermal simulations were employed in 3D generalized homogeneous and patient-specific heterogeneous models, for three-phased array transducers with 32, 64, and 128 elements, composed of sectored 4, 8, and 16 tubular ring transducers, respectively. The applicator performance was characterized as a function of array configuration, focal depth, phasing modes, and balloon reflector geometry. A 16-element proof-of-concept phased array applicator assembly, consisting of four tubular transducers each divided into four sectors, was fabricated, and characterized with hydrophone measurements along and across the axis, and ablations in ex vivo tissue.

RESULTS

Simulation results indicated that transducer arrays (1.5 MHz, 9 mm OD × 20 mm long), balloon sizes (41-50 mm expanded diameter, 20-60 mm focal depth), phasing mode (0-4) and sonication duration (30 s) can produce spatially localized acoustic intensity focal patterns (focal length: 3-22 mm, focal width: 0.7-8.7 mm) and ablative thermal lesions (width: 2.7-16 mm, length: 6-46 mm) in pancreatic tissue across a 10-90 mm focal depth range. Patient-specific studies indicated that 0.1, 0.46, and 1.2 cm³ volume of tumor can be ablated in the body of the pancreas for 120 s sonications using a single axial focus (Mode 0), or four, and eight simultaneous foci in a toroidal pattern (Mode 2 and 4, respectively). Hydrophone measurements demonstrated good agreement with simulation. Experiments in which chicken meat was thermally ablated indicated that volumetric ablation can be produced using single or multiple foci.

CONCLUSIONS

The results of this study demonstrated the feasibility of a novel compact ultrasound applicator design capable of focusing, deep penetration, electronic steering, and volumetric thermal ablation. The proposed applicator can be used for compact endoluminal or laparoscopic delivery of localized ultrasound energy to deep-seated targets.

Deployable ultrasound applicators for endoluminal delivery of volumetric hyperthermia

PURPOSE

To investigate the design of an endoluminal deployable ultrasound applicator for delivering volumetric hyperthermia to deep tissue sites as a possible adjunct to radiation and chemotherapy.

METHOD

This study considers an ultrasound applicator consisting of two tubular transducers situated at the end of a catheter assembly, encased within a distensible conical shaped balloon-based reflector that redirects acoustic energy distally into the tissue. The applicator assembly can be inserted endoluminally or laparoscopically in a compact form and expanded after delivery to the target site. Comprehensive acoustic and biothermal simulations and parametric studies were employed in generalized 3D and patient-specific pancreatic head and body tumor models to characterize the acoustic performance and evaluate heating capabilities of the applicator by investigating the device at a range of operating frequencies, tissue acoustic and thermal properties, transducer configurations, power modulation, applicator positioning, and by analyzing the resultant 40, 41, and 43 °C isothermal volumes and penetration depth of the heating volume. Intensity distributions and volumetric temperature contours were calculated to define moderate hyperthermia boundaries.

RESULTS

Parametric studies demonstrated the frequency selection to control volume and depth of therapeutic heating from 62 to 22 cm³ and 4 to 2.6 cm as frequency ranges from 1 MHz to 4.7 MHz, respectively. Width of the heating profile tracks closely with the aperture. Water cooling within the reflector balloon was effective in controlling temperature to 37 °C maximum within the luminal wall. Patient-specific studies indicated that applicators with extended OD in the range of 3.6-6.2 cm with 0.5-1 cm long and 1 cm OD transducers can heat volumes of 1.1-7 cm³, 3-26 cm³, and 3.3-37.4 cm³ of pancreatic body and head tumors above 43, 41, and 40 °C, respectively.

CONCLUSION

In silico studies demonstrated the feasibility of combining endoluminal ultrasound with an integrated expandable balloon reflector for delivering volumetric hyperthermia in regions adjacent to body lumens and cavities.

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.

3D synthetic aperture imaging with a therapeutic spherical random phased array for transcostal applications

Experimental validation of a synthetic aperture imaging technique using a therapeutic random phased array is described, demonstrating the dual nature of imaging and therapy of such an array. The transducer is capable of generating both continuous wave high intensity beams for ablating the tumor and low intensity ultrasound pulses to image the target area. Pulse-echo data is collected from the elements of the phased array to obtain B-mode images of the targets. Since therapeutic arrays are optimized for therapy only with concave apertures having low f-number and large directive elements often coarsely sampled, imaging can not be performed using conventional beamforming. We show that synthetic aperture imaging is capable of processing the acquired RF data to obtain images of the field of interest. Simulations were performed to compare different synthetic aperture imaging techniques to identify the best algorithm in terms of spatial resolution. Experimental validation was performed using a 1 MHz, 256-elements, spherical random phased array with 130 mm radius of curvature. The array was integrated with a research ultrasound scanner via custom connectors to acquire raw RF data for variety of targets. Imaging was implemented using synthetic aperture beamforming to produce images of a rib phantom and ex vivo ribs. The array was shown to resolve spherical targets within ±15 mm of either side of the axis in the focal plane and obtain 3D images of the rib phantom up to ±40 mm of either side of the central axis and at a depth of 3-9 cm from the array surface. The lateral and axial full width half maximum was 1.15 mm and 2.75 mm, respectively. This study was undertaken to emphasize that both therapy and image guidance with a therapeutic random phased array is possible and such a system has the potential to address some major limitations in the existing high intensity focused ultrasound (HIFU) systems. The 3D images obtained with a therapeutic array can be used to identify and locate strong scattering objects aiding to image guidance and treatment planning of the HIFU procedure.

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.

Self-adaptive ultrasonic beam amplifiers: application to transcostal shock wave therapy

Ultrasound shock wave therapy is increasingly used for non-invasive surgery. It requires the focusing of very high pressure amplitude in precisely controlled focal spots. In transcostal therapy of the heart or the liver, the high impedance mismatch between the bones and surrounding tissues gives rise to strong aberrations and attenuation of the therapeutic wavefront, with potential risks of injury at the tissue-bone interface. An adaptive propagation of the ultrasonic beam through the intercostal spaces would be required. Several solutions have been developed so far, but they require a prior knowledge of the patient's anatomy or an invasive calibration process, not applicable in clinic. Here, we develop a non-invasive adaptive focusing method for ultrasound therapy through the ribcage using a time reversal cavity (TRC) acting as an ultrasonic beam amplifier. This method is based on ribcage imaging through the TRC and a projection orthogonally to the strongest identified reflectors. The focal pressure of our device was improved by up to 30% using such self-adaptive processing, without degrading the focal spots size and shape. This improvement allowed lesion formation in an Ultracal^® phantom through a ribcage without invasive calibration of the device. This adaptive method could be particularly interesting to improve the efficiency and the safety of pulsed cavitational therapy of the heart or the liver.

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.

Computational study on the propagation of strongly focused nonlinear ultrasound in tissue with rib-like structures

The presence of a rib cage is a significant hindrance to the potential applications of focused ultrasound as a noninvasive extracorporeal surgery modality for various internal organs. Here the influence of ribs on the propagation of strongly focused high-intensity nonlinear ultrasound beam inside the body is studied. Based on the spheroidal beam equation, a three-dimensional numerical algorithm is developed to solve the nonlinear acoustic field generated by a focused ultrasonic transducer with a large aperture angle. Idealized ribs, of rectangular cross sectional, with high absorption and impedance, and various dimensions, are used to simulate human anatomical configurations. The changes in the spatial distribution of acoustic intensity and the reduction of the acoustic pressure amplitude and heat deposition rate due to the presence of "ribs" are investigated. It is somewhat surprising that in some cases, the axial peak positions shift less than 2 mm and more than 80% of the sound energy can propagate through the space of the rib cage in the strongly focused sound field. This study also includes quantitative analyses of the effects of different rib configurations and transducers of various f-numbers. The results can be used as reference information for further study and clinical applications.

Pre-clinical study of in vivo magnetic resonance-guided bubble-enhanced heating in pig liver

Bubble-enhanced heating (BEH) can be exploited to increase heating efficiency in treatment of liver tumors with non-invasive high-intensity focused ultrasound (HIFU). The objectives of this study were: (i) to demonstrate the feasibility of increasing the heating efficiency of sonication exploiting BEH in pig liver in vivo using a clinical platform; (ii) to determine the acoustic threshold for such effects with real-time, motion-compensated magnetic resonance-guided thermometry; and (iii) to compare the heating patterns and thermal lesion characteristics resulting from continuous sonication and sonication including a burst pulse. The threshold acoustic power for generation of BEH in pig liver in vivo was determined using sonication of 0.5-s duration ("burst pulse") under real-time magnetic resonance thermometry. In a second step, experimental sonication composed of a burst pulse followed by continuous sonication (14.5 s) was compared with conventional sonication (15 s) of identical energy (1.8 kJ). Modification of the heating pattern at the targeted region located at a liver depth between 20 and 25 mm required 600-800 acoustic watts. The experimental group exhibited near-spherical heating with 40% mean enhancement of the maximal temperature rise as compared with the conventional sonication group, a mean shift of 7 ± 3.3 mm toward the transducer and reduction of the post-focal temperature increase. Magnetic resonance thermometry can be exploited to control acoustic BEH in vivo in the liver. By use of experimental sonication, more efficient heating can be achieved while protecting tissues located beyond the focal point.

The road to clinical use of high-intensity focused ultrasound for liver cancer: technical and clinical consensus

Clinical use of high-intensity focused ultrasound (HIFU) under ultrasound or MR guidance as a non-invasive method for treating tumors is rapidly increasing. Tens of thousands of patients have been treated for uterine fibroid, benign prostate hyperplasia, bone metastases, or prostate cancer. Despite the methods' clinical potential, the liver is a particularly challenging organ for HIFU treatment due to the combined effect of respiratory-induced liver motion, partial blocking by the rib cage, and high perfusion/flow. Several technical and clinical solutions have been developed by various groups during the past 15 years to compensate for these problems. A review of current unmet clinical needs is given here, as well as a consensus from a panel of experts about technical and clinical requirements for upcoming pilot and pivotal studies in order to accelerate the development and adoption of focused ultrasound for the treatment of primary and secondary liver cancer.

Transcostal high-intensity focused ultrasound treatment using phased array with geometric correction

In the high-intensity focused ultrasound treatment of liver tumors, ultrasound propagation is affected by the rib cage. Because of the diffraction and absorption of the bone, the sound distribution at the focal plane is altered, and more importantly, overheating on the rib surface might occur. To overcome these problems, a geometric correction method is applied to turn off the elements blocked by the ribs. The potential of steering the focus of the phased-array along the propagation direction to improve the transcostal treatment was investigated by simulations and experiments using different rib models and transducers. The ultrasound propagation through the ribs was computed by a hybrid method including the Rayleigh-Sommerfeld integral, k-space method, and angular spectrum method. A modified correction method was proposed to adjust the output of elements based on their relative area in the projected "shadow" of the ribs. The simulation results showed that an increase in the specific absorption rate gain up to 300% was obtained by varying the focal length although the optimal value varied in each situation. Therefore, acoustic simulation is required for each clinical case to determine a satisfactory treatment plan.

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.