Elastographic versus x-ray CT imaging of radio frequency ablation coagulations: Anin vitrostudy

The Role of Artificial Intelligence and Texture Analysis in Interventional Radiological Treatments of Liver Masses: A Narrative Review

Liver lesions, including both benign and malignant tumors, pose significant challenges in interventional radiological treatment planning and prognostication. The emerging field of artificial intelligence (AI) and its integration with texture analysis techniques have shown promising potential in predicting treatment outcomes, enhancing precision, and aiding clinical decision-making. This comprehensive review aims to summarize the current state-of-the-art research on the application of AI and texture analysis in determining treatment response, recurrence rates, and overall survival outcomes for patients undergoing interventional radiological treatment for liver lesions. Furthermore, the review addresses the challenges associated with the implementation of AI and texture analysis in clinical practice, including data acquisition, standardization of imaging protocols, and model validation. Future directions and potential advancements in this field are discussed. Integration of multi-modal imaging data, incorporation of genomics and clinical data, and the development of predictive models with enhanced interpretability are proposed as potential avenues for further research. In conclusion, the application of AI and texture analysis in predicting outcomes of interventional radiological treatment for liver lesions shows great promise in augmenting clinical decision-making and improving patient care. By leveraging these technologies, clinicians can potentially enhance treatment planning, optimize intervention strategies, and ultimately improve patient outcomes in the management of liver lesions.

Subregion Radiomics Analysis to Display Necrosis After Hepatic Microwave Ablation-A Proof of Concept Study

OBJECTIVES

The aim of this study was to improve the visualization of coagulation necrosis after computed tomography (CT)-guided microwave ablation (MWA) in routine postablational imaging.

MATERIALS AND METHODS

Ten MWAs were performed in 8 pigs under CT guidance. After each ablation, we obtained contrast-enhanced CT scans in venous phase. Ablations were then resected as a whole, and histologic slices were obtained orthogonally through the ablation center. Subsequently, a vital stain was applied to the sections for visualization of coagulation necrosis. Computed tomography images were reformatted to match the histologic slices. Afterwards, quantitative imaging features were extracted from the subregions of all images, and binary classifiers were used to predict the presence of coagulation necrosis for each subregion. From this, heatmaps could be created, which visually represented the extent of necrosis in each CT image. Two independent observers evaluated the extent of coagulative necrosis between the heat maps and histological sections.

RESULTS

We applied 4 different classifiers, including a generalized linear mixed model (GLMM), a stochastic gradient boosting classifier, a random forest classifier, and a k-nearest neighbor classifier, out of which the GLMM showed the best performance to display coagulation necrosis. The GLMM resulted in an area under the curve of 0.84 and a Jaccard index of 0.6 between the generated heat map and the histologic reference standard as well as a good interobserver agreement with a Jaccard index of 0.9.

CONCLUSIONS

Subregion radiomics analysis may improve visualization of coagulation necrosis after hepatic MWA in an in vivo porcine model.

Improved Visualization of the Necrotic Zone after Microwave Ablation Using Computed Tomography Volume Perfusion in an In Vivo Porcine Model

After hepatic microwave ablation, the differentiation between fully necrotic and persistent vital tissue through contrast enhanced CT remains a clinical challenge. Therefore, there is a need to evaluate new imaging modalities, such as CT perfusion (CTP) to improve the visualization of coagulation necrosis. MWA and CTP were prospectively performed in five healthy pigs. After the procedure, the pigs were euthanized, and the livers explanted. Orthogonal histological slices of the ablations were stained with a vital stain, digitalized and the necrotic core was segmented. CTP maps were calculated using a dual-input deconvolution algorithm. The segmented necrotic zones were overlaid on the DICOM images to calculate the accuracy of depiction by CECT/CTP compared to the histological reference standard. A receiver operating characteristic analysis was performed to determine the agreement/true positive rate and disagreement/false discovery rate between CECT/CTP and histology. Standard CECT showed a true positive rate of 81% and a false discovery rate of 52% for display of the coagulation necrosis. Using CTP, delineation of the coagulation necrosis could be improved significantly through the display of hepatic blood volume and hepatic arterial blood flow (p < 0.001). The ratios of true positive rate/false discovery rate were 89%/25% and 90%/50% respectively. Other parameter maps showed an inferior performance compared to CECT.

Two-dimensional ultrasound-computed tomography image registration for monitoring percutaneous hepatic intervention

PURPOSE

Deformable registration of ultrasound (US) and contrast enhanced computed tomography (CECT) images are essential for quantitative comparison of ablation boundaries and dimensions determined using these modalities. This comparison is essential as stiffness-based imaging using US has become popular and offers a nonionizing and cost-effective imaging modality for monitoring minimally invasive microwave ablation procedures. A sensible manual registration method is presented that performs the required CT-US image registration.

METHODS

The two-dimensional (2D) virtual CT image plane that corresponds to the clinical US B-mode was obtained by "virtually slicing" the 3D CT volume along the plane containing non-anatomical landmarks, namely points along the microwave ablation antenna. The initial slice plane was generated using the vector acquired by rotating the normal vector of the transverse (i.e., xz) plane along the angle subtended by the antenna. This plane was then further rotated along the ablation antenna and shifted along with the direction of normal vector to obtain similar anatomical structures, such as the liver surface and vasculature that is visualized on both the CT virtual slice and US B-mode images on 20 patients. Finally, an affine transformation was estimated using anatomic and non-anatomic landmarks to account for distortion between the colocated CT virtual slice and US B-mode image resulting in a final registered CT virtual slice. Registration accuracy was measured by estimating the Euclidean distance between corresponding registered points on CT and US B-mode images.

RESULTS

Mean and SD of the affine transformed registration error was 1.85 ± 2.14 (mm), computed from 20 coregistered data sets.

CONCLUSIONS

Our results demonstrate the ability to obtain 2D virtual CT slices that are registered to clinical US B-mode images. The use of both anatomical and non-anatomical landmarks result in accurate registration useful for validating ablative margins and comparison to electrode displacement elastography based images.

Monitoring Radiofrequency Ablation Using Ultrasound Envelope Statistics and Shear Wave Elastography in the Periablation Period: An In Vitro Feasibility Study

Radiofrequency ablation (RFA) is a minimally invasive method for treating tumors. Shear wave elastography (SWE) has been widely applied in evaluating tissue stiffness and final ablation size after RFA. However, the usefulness of periablation SWE imaging in assessing RFA remains unclear. Therefore, this study investigated the correlation between periablation SWE imaging and final ablation size. An in vitro porcine liver model was used for experimental validation (n = 36). During RFA with a power of 50 W, SWE images were collected using a clinical ultrasound system. To evaluate the effects of tissue temperature and gas bubbles during RFA, changes in the ablation temperature were recorded, and image echo patterns were measured using B-mode and ultrasound statistical parametric images. After RFA, the gross pathology of each tissue sample was compared with the region of change in the corresponding periablation SWE image. The experimental results showed that the tissue temperature at the ablation site varied between 70°C and 100°C. Hyperechoic regions and changes were observed in the echo amplitude distribution induced by gas bubbles. Under this condition, the confounding effects (including the temperature increase, tissue stiffness increase, and presence of gas bubbles) resulted in artifacts in the periablation SWE images, and the corresponding region correlated with the estimated final ablation size obtained from the gross pathology (r = 0.8). The findings confirm the feasibility of using periablation SWE imaging in assessing RFA.

Quantifying local stiffness variations in radiofrequency ablations with dynamic indentation

Elastographic imaging can be used to monitor ablation procedures; however, confident and clear determination of the ablation boundary is essential to ensure complete treatment of the pathological target. To investigate the potential for ablation boundary representation on elastographic images, local variations in the viscoelastic properties in radiofrequency-ablated regions that were formed in vivo in porcine liver tissue were quantified using dynamic indentation. Spatial stiffness maps were then correlated to stained histology, the gold standard for the determination of the ablation periphery or boundary. Regions of interest in 11 radiofrequency ablation samples were indented at 18-24 locations each, including the central zone of complete necrosis and more peripheral transition zones including normal tissue. Storage modulus and the rate of stiffening were both greatest in the central ablation zone and decreased with radial distance away from the center. The storage modulus and modulus contrast at the ablation outer transition zone boundary were 3.1 ± 1.0 kPa and 1.6 ± 0.4, respectively, and 36.2 ± 9.1 kPa and 18.3 ± 5.5 at the condensation boundary within the ablation zone. Elastographic imaging modalities were then compared to gross pathology in ex vivo bovine liver tissue. Area estimated from strain, shear-wave velocity, and gross pathology images were 470, 560, and 574 mm(2), respectively, and ablation widths were 19.4, 20.7, and 23.0 mm. This study has provided insights into spatial stiffness distributions within radiofrequency ablations and suggests that low stiffness contrast on the ablation periphery leads to the observed underestimation of ablation extent on elastographic images.

Ultrasound-based relative elastic modulus imaging for visualizing thermal ablation zones in a porcine model

The feasibility of using ultrasound-based elastic modulus imaging to visualize thermal ablation zones in an in vivo porcine model is reported. Elastic modulus images of soft tissues are estimated as an inverse optimization problem. Ultrasonically measured displacement data are utilized as inputs to determine an elastic modulus distribution that provides the best match to this displacement field. A total of 14 in vivo thermal ablation zones were investigated in this study. To determine the accuracy of delineation of each thermal ablation zone using elastic modulus imaging, the dimensions (lengths of long and short axes) and the area of each thermal ablation zone obtained from an elastic modulus image were compared to the corresponding gross pathology photograph of the same ablation zone. Comparison of elastic modulus imaging measurements and gross pathology measurements showed high correlation with respect to the area of thermal ablation zones (Pearson coefficient = 0.950 and p < 0.0001). The radiological-pathological correlation was slightly lower (correlation = 0.853, p < 0.0001) for strain imaging among these 14 in vivo ablation zones. We also found that, on average, elastic modulus imaging can more accurately depict thermal ablation zones, when compared to strain imaging (14.7% versus 22.3% absolute percent error in area measurements, respectively). Furthermore, elastic modulus imaging also provides higher (more than a factor of 2) contrast-to-noise ratios for evaluating these thermal ablation zones than those on corresponding strain images, thereby reducing inter-observer variability. Our preliminary results suggest that elastic modulus imaging might potentially enhance the ability to visualize thermal ablation zones, thereby improving assessment of ablative therapies.

Multilevel hybrid 2D strain imaging algorithm for ultrasound sector/phased arrays

Two-dimensional (2D) cross-correlation algorithms are necessary to estimate local displacement vector information for strain imaging. However, most of the current two-dimensional cross-correlation algorithms were developed for linear array transducers. Although sector and phased array transducers are routinely used for clinical imaging of abdominal and cardiac applications, strain imaging for these applications has been performed using one-dimensional (1D) cross-correlation analysis. However, one-dimensional cross-correlation algorithms are unable to provide accurate and precise strain estimation along all the angular insonification directions which can range from -45 degrees to 45 degrees with sector and phased array transducers. In addition, since sector and phased array based images have larger separations between beam lines as the pulse propagates deeper into tissue, signal decorrelation artifacts with deformation or tissue motion are more pronounced. In this article, the authors propose a multilevel two-dimensional hybrid algorithm for ultrasound sector and phased array data that demonstrate improved tracking and estimation performance when compared to the traditional 1D cross-correlation or 2D cross-correlation based methods. Experimental results demonstrate that the signal-to-noise and contrast-to-noise ratio estimates improve significantly for smaller window lengths for the hybrid method when compared to the currently used one-dimensional or two-dimensional cross-correlation algorithms. Strain imaging results on ex vivo thermal lesions created in liver tissue and in vivo on cardiac short-axis views demonstrate the improved image quality obtained with this method.

Young's modulus reconstruction for radio-frequency ablation electrode-induced displacement fields: a feasibility study

Radio-frequency (RF) ablation is a minimally invasive treatment for tumors in various abdominal organs. It is effective if good tumor localization and intraprocedural monitoring can be done. In this paper, we investigate the feasibility of using an ultrasound-based Young's modulus reconstruction algorithm to image an ablated region whose stiffness is elevated due to tissue coagulation. To obtain controllable tissue deformations for abdominal organs during and/or intermediately after the RF ablation, the proposed modulus imaging method is specifically designed for using tissue deformation fields induced by the RF electrode. We have developed a new scheme under which the reconstruction problem is simplified to a 2-D problem. Based on this scheme, an iterative Young's modulus reconstruction technique with edge-preserving regularization was developed to estimate the Young's modulus distribution. The method was tested in experiments using a tissue-mimicking phantom and on ex vivo bovine liver tissues. Our preliminary results suggest that high contrast modulus images can be successfully reconstructed. In both experiments, the geometries of the reconstructed modulus images of thermal ablation zones match well with the phantom design and the gross pathology image, respectively.

Three-dimensional electrode displacement elastography using the Siemens C7F2 fourSight four-dimensional ultrasound transducer

Because ablation therapy alters the elastic modulus of tissues, emerging strain imaging methods may enable clinicians for the first time to have readily available, cost-effective, real-time guidance to identify the location and boundaries of thermal lesions. Electrode displacement elastography is a method of strain imaging tailored specifically to ultrasound-guided electrode-based ablative therapies (e.g., radio-frequency ablation). Here tissue deformation is achieved by applying minute perturbations to the unconstrained end of the treatment electrode, resulting in localized motion around the end of the electrode embedded in tissue. In this article, we present a method for three-dimensional (3D) elastographic reconstruction from volumetric data acquired using the C7F2 fourSight four-dimensional ultrasound transducer, provided by Siemens Medical Solutions USA, Inc. (Issaquah, WA, USA). Lesion reconstruction is demonstrated for a spherical inclusion centered in a tissue-mimicking phantom, which simulates a thermal lesion embedded in a normal tissue background. Elastographic reconstruction is also performed for a thermal lesion created in vitro in canine liver using radio-frequency ablation. Postprocessing is done on the acquired raw radio-frequency data to form surface-rendered 3D elastograms of the inclusion. Elastographic volume estimates of the inclusion compare reasonably well with the actual known inclusion volume, with 3D electrode displacement elastography slightly underestimating the true inclusion volume.

Focusing cross-fire applicator for ultrasonic hyperthermia of tumors

An improved concept for ultrasonic hyperthermia of tumors is presented. This concept is based on past experience of a German government supported project , which ended in 1984. It offers a low cost alternative to common RF- and microwave methods for hyperthermia of tumors with volumes between 1 and 40 ml at treatment times between 30 and 60 min. Our new version of the system considerably improves the temperature suppression in the healthy tissue around the target area and enables the adjustment of the beam width to the actual tumor size and the field geometry to the depth and shape of the tumor. The applicator can be used for moderate hyperthermia with tissue overheating up to 10K or for ablation therapy with short high temperature pulses. Its central area is free for the integration of a commercial ultrasonic diagnostic sector scanner or a Doppler flow sensor in order to support the adjustment of the transducer and to monitor the whole area during the therapy.

Liver ablation guidance with acoustic radiation force impulse imaging: challenges and opportunities

Previous studies have established the feasibility of monitoring radiofrequency (RF) ablation procedures with acoustic radiation force impulse (ARFI) imaging. However, questions remained regarding the utility of the technique in clinically realistic scenarios and at scanning depths associated with abdominal imaging in adults. We address several of these issues and detail recent progress towards the clinical relevance of the ARFI technique. Results from in vitro bovine tissues and an in vivo ovine model are presented. Additional experiments were conducted with a tissue-mimicking phantom and parallel receive tracking techniques in order to further support the clinical feasibility of the method. Thermal lesions created during RF ablation are visualized with high contrast in both in vitro and in vivo hepatic tissues, and radial lesion growth can be monitored throughout the duration of the procedure. ARFI imaging is implemented on a diagnostic ultrasonic scanner, and thus may be a convenient option to guide RF ablation procedures, particularly when electrode insertion is also performed with sonographic guidance.

An investigation of the use of transmission ultrasound to measure acoustic attenuation changes in thermal therapy

The potential of using a commercial ultrasound transmission imaging system to quantitatively monitor tissue attenuation changes after thermal therapy was investigated. The ultrasound transmission imaging system used, the AcoustoCam (Imperium Inc., MD) allows ultrasonic images to be captured using principles similar to that of a CCD-type camera that collects light. Ultrasound energy is focused onto a piezoelectric array by an acoustic lens system, creating a gray scale acoustic image. In this work, the pixel values from the acoustic images were assigned acoustic attenuation values by imaging polyacrylamide phantoms of varying known attenuation. After the calibration procedure, data from heated polyacrylamide/bovine serum albumin (BSA) based tissue-mimicking (TM) phantoms and porcine livers were acquired. Samples were heated in water at temperatures of 35, 45, 55, 65, and 75 degrees C for 1 h. Regions of interest were chosen in the images and acoustic attenuation values before and after heating were compared. An increase in ultrasound attenuation was found in phantoms containing BSA and in porcine liver. In the presence of BSA, attenuation in the TM phantom increased by a factor of 1.5, while without BSA no significant changes were observed. The attenuation of the porcine liver increased by up to a factor of 2.4, consistent with previously reported studies. The study demonstrates the feasibility of using a quantitative ultrasound transmission imaging system for monitoring thermal therapy.

Segmentation of elastographic images using a coarse-to-fine active contour model

Delineation of radiofrequency-ablation-induced coagulation (thermal lesion) boundaries is an important clinical problem that is not well addressed by conventional imaging modalities. Elastography, which produces images of the local strain after small, externally applied compressions, can be used for visualization of thermal coagulations. This paper presents an automated segmentation approach for thermal coagulations on 3-D elastographic data to obtain both area and volume information rapidly. The approach consists of a coarse-to-fine method for active contour initialization and a gradient vector flow, active contour model for deformable contour optimization with the help of prior knowledge of the geometry of general thermal coagulations. The performance of the algorithm has been shown to be comparable to manual delineation of coagulations on elastograms by medical physicists (r = 0.99 for volumes of 36 radiofrequency-induced coagulations). Furthermore, the automatic algorithm applied to elastograms yielded results that agreed with manual delineation of coagulations on pathology images (r = 0.96 for the same 36 lesions). This algorithm has also been successfully applied on in vivo elastograms.

Temperature dependence of ultrasonic backscattered energy in motion-compensated images

Noninvasive temperature imaging would enhance the ability to uniformly heat tumors at therapeutic levels. Ultrasound is an attractive modality for this purpose. Previously, we predicted monotonic changes in backscattered energy (CBE) of ultrasound with temperature for certain subwavelength scatterers. We also measured CBE values similar to our predictions in bovine liver, turkey breast, and pork muscle in one dimension (1-D). Those measurements were corrected manually for changes in the axial position of echo signals with temperature. To investigate the effect of temperature on CBE in 2-D, we imaged 1-cm thick samples of bovine liver, turkey breast, and pork muscle during heating in a water bath. Images were formed by a phased-array imager with a 7 MHz linear probe. Using radio frequency (RF) signals permitted the use of cross correlation as a similarity measure for automatic tracking of feature displacement as a function of temperature. Feature displacement across the specimen was nonuniform with typical total displacements of 0.5 mm in both axial and lateral directions. Apparent movement in eight image regions in each specimen was tracked from 37 to 50 degrees C in 0.5 degrees C steps. Envelopes of motion-compensated image regions were found then smoothed with a 3 x 3 running average filter before forming the backscattered energy at each pixel. Our measure of CBE compared means of both the positive and negative changes in the backscattered energy (BE) images. CBE was monotonic and differed by about 4 dB at 50 degrees C from its value at 37 degrees C. Relatively noise-free CBE curves from tissue volumes of less than 1 cm3 supports the use of CBE for temperature estimation.

Monitoring stiffness changes in lesions after radiofrequency ablation at different temperatures and durations of ablation

The variations in the stiffness or stiffness contrast of lesions resulting from radiofrequency (RF) ablation of canine liver tissue at different temperatures and for different ablation durations at a specified temperature are analyzed. Tissue stiffness, in general, increases with temperature; however, an anomaly exists around 80 degrees C, where the stiffness of the lesion is lower than that of the lesion ablated at 70 degrees C. On the other hand, the stiffness increases monotonically with the duration of ablation. Plots illustrating the ratio of mean strains in normal canine liver tissue to mean strains in ablated thermal lesions demonstrate the variation in the stiffness contrast of the thermal lesions. The contrast-to-noise ratio (CNRe) of the lesions, which serves as an indicator of the detectability of the lesions under the different experimental imaging conditions described above, is also presented. The results presented in this paper show that the elastographic depiction of stiffer thermal lesions is better, in terms of the CNRe parameter. An important criterion in the elastographic depiction of RF-ablated regions of tissue is the trade-off between ablation temperature and duration of ablation. Tissue necrosis can occur either by ablating tissue to high temperatures for short durations or to lower temperatures for longer durations. In this paper, we attempt to characterize the elastographic depiction of thermal lesions under these different experimental conditions. This paper provides results that may be utilized by practitioners of RF ablation to decide the ablation temperature and duration, on the basis of the strain images of normal liver tissue and ablated thermal lesions discussed in this paper.