Electrode displacement strain imaging of thermally-ablated liver tissue in an in vivo animal model

Ultrasound Entropy Imaging for Detection and Monitoring of Thermal Lesion During Microwave Ablation of Liver

Ultrasonic B-mode imaging offers non-invasive and real-time monitoring of thermal ablation treatment in clinical use, however it faces challenges of moderate lesion-normal contrast and detection accuracy. Quantitative ultrasound imaging techniques have been proposed as promising tools to evaluate the microstructure of ablated tissue. In this study, we introduced Shannon entropy, a non-model based statistical measurement of disorder, to quantitatively detect and monitor microwave-induced ablation in porcine livers. Performance of typical Shannon entropy (TSE), weighted Shannon entropy (WSE), and horizontally normalized Shannon entropy (hNSE) were explored and compared with conventional B-mode imaging. TSE estimated from non-normalized probability distribution histograms was found to have insufficient discernibility of different disorder of data. WSE that improves from TSE by adding signal amplitudes as weights obtained area under receiver operating characteristic (AUROC) curve of 0.895, whereas it underestimated the periphery of lesion region. hNSE provided superior ablated area prediction with the correlation coefficient of 0.90 against ground truth, AUROC of 0.868, and remarkable lesion-normal contrast with contrast-to-noise ratio of 5.86 which was significantly higher than other imaging methods. Data distributions shown in horizontally normalized probability distribution histograms indicated that the disorder of backscattered envelope signal from ablated region increased as treatment went on. These findings suggest that hNSE imaging could be a promising technique to assist ultrasound guided percutaneous thermal ablation.

Differential Imaging of Liver Tumors before and after Microwave Ablation with Electrode Displacement Elastography

Liver cancer is a leading cause of cancer-related deaths; however, primary treatment options such as surgical resection and liver transplant may not be viable for many patients. Minimally invasive image-guided microwave ablation (MWA) provides a locally effective treatment option for these patients with an impact comparable to that of surgery for both cancer-specific and overall survival. MWA efficacy is correlated with accurate image guidance; however, conventional modalities such as B-mode ultrasound and computed tomography have limitations. Alternatively, ultrasound elastography has been used to demarcate post-ablation zones, yet has limitations for pre-ablation visualization because of variability in strain contrast between cancer types. This study attempted to characterize both pre-ablation tumors and post-ablation zones using electrode displacement elastography (EDE) for 13 patients with hepatocellular carcinoma or liver metastasis. Typically, MWA ablation margins of 0.5-1.0 cm are desired, which are strongly correlated with treatment efficacy. Our results revealed an average estimated ablation margin inner quartile range of 0.54-1.21 cm with a median value of 0.84 cm. These treatment margins lie within or above the targeted ablative margin, indicating the potential to use EDE for differentiating index tumors and ablated zones during clinical ablations. We also obtained a high correlation between corresponding segmented cross-sectional areas from contrast-enhanced computed tomography, the current clinical gold standard, when compared with EDE strain images, with r2 values of 0.97 and 0.98 for pre- and post-ablation regions.

Elastography for characterization of focal liver lesions: current evidence and future perspectives

Focal liver lesions (FLLs) are a common finding during routine abdominal ultrasound (US). The differential diagnosis between diverse types of FLLs, especially between benign and malignant ones, is extremely important and can often be particularly challenging. Radiological techniques with contrast administration and/or liver biopsy are mostly necessary for establishing diagnosis, but they have several contraindications or complications. Due to limitations of these tools, there is urgent and still unmet need to develop a first line, non-invasive and simple method to diagnose FLLs. Elastography is an US-based imaging modality that provides information about the physical parameter corresponding to the tissue stiffness and can be considered a virtual biopsy. Several elastographic approaches have been developed, such as transient elastography, strain imaging and share wave imaging, which include point shear wave elastography and 2D shear wave elastography. These tools are already in use for evaluating liver fibrosis and in the assessment of focal lesions in other organs, like breast and thyroid gland. This review aims to assess the current evidence of different techniques based on elastography in the setting of FLLs, in order to evaluate accuracy, limitations and future perspectives. In particular, we focused on two contexts: the ability of discriminating between benign and malignant lesions, especially hepatocellular carcinoma and liver metastasis, and the surveillance after percutaneous therapy. This could have a high clinical impact making elastography crucial to identify the appropriate management of FLLs.

Adaptation of Dictionary Learning for Electrode Displacement Elastography. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2020;2020:2023-2026. [PMID: 33018401 PMCID: PMC7538652 DOI: 10.1109/embc44109.2020.9175319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Microwave ablation has become a common treatment method for liver cancers. Unfortunately, microwave ablation success is correlated with clinician's ability for proper electrode placement and assess ablative margins, requiring accurate imaging of liver tumors and ablated zones. Conventionally, ultrasound and computed tomography are utilized for this purpose, yet both have their respective drawbacks. As an alternate approach, electrode displacement elastography offers promise but is still plagued by decorrelation artifacts reducing lesion depiction and visualization. A recent filtering method, namely dictionary representation, has improved contrast-to-noise ratios without reducing delineation contrast. As a supplement to this recent work, this paper evaluates adaptations on this initial dictionary-learning algorithm and applies them to an EDE phantom and 15 in-vivo patient datasets. Two new adaptations of dictionary representations were evaluated, namely a combined dictionary and magnitude-based dictionary representation. When comparing numerical results, the combined dictionary representation algorithm outperforms the previous developed dictionary representation in signal-to-noise (1.54 dB) and contrast-to-noise (0.67 dB) ratios, while a magnitude dictionary representation produces higher noise levels, but improves visualized strain tensor resolution.

Physiological Motion Reduction Using Lagrangian Tracking for Electrode Displacement Elastography

Minimally invasive treatments such as microwave ablation (MWA) have been growing in popularity for extending liver cancer survival rates in patients, when surgery is not an option. As a non-ionizing, real-time alternative to contrast-enhanced computed tomography, electrode displacement elastography (EDE) has shown promise as an imaging modality for MWA. Despite imaging efficacy, motion artifacts caused by physiological motion result in unintended speckle pattern variance, thereby inhibiting consistent and accurate ablated region visualization. To combat these unavoidable motion artifacts, a Lagrangian deformation tracking (LDT) approach based on freehand EDE was developed to track tissue movement and better define tissue properties. For validating LDT efficacy, a spherical inclusion phantom as well as seven in vivo data sets were processed, and strain tensor images were compared with identical time sampled images estimated using a traditional Eulerian approach. In vivo results revealed greater consistency among visualized LDT strain tensor images, with segmented ablated regions exhibiting standard deviation reductions of up to 98% when compared with Eulerian strain tensor images. Additionally, Lagrangian strain tensor images provided Dice coefficient improvements up to 25%, and success rates improved from approximately 50% to nearly 100% for ablated region visualization.

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.

Comparison of Displacement Tracking Algorithms for in Vivo Electrode Displacement Elastography

Hepatocellular carcinoma and liver metastases are common hepatic malignancies presenting with high mortality rates. Minimally invasive microwave ablation (MWA) yields high success rates similar to surgical resection. However, MWA procedures require accurate image guidance during the procedure and for post-procedure assessments. Ultrasound electrode displacement elastography (EDE) has demonstrated utility for non-ionizing imaging of regions of thermal necrosis created with MWA in the ablation suite. Three strategies for displacement vector tracking and strain tensor estimation, namely coupled subsample displacement estimation (CSDE), a multilevel 2-D normalized cross-correlation method, and quality-guided displacement tracking (QGDT) have previously shown accurate estimations for EDE. This paper reports on a qualitative and quantitative comparison of these three algorithms over 79 patients after an MWA procedure. Qualitatively, CSDE presents sharply delineated, clean ablated regions with low noise except for the distal boundary of the ablated region. Multilevel and QGDT contain more visible noise artifacts, but delineation is seen over the entire ablated region. Quantitative comparison indicates CSDE with more consistent mean and standard deviations of region of interest within the mass of strain tensor magnitudes and higher contrast, while Multilevel and QGDT provide higher CNR. This fact along with highest success rates of 89% and 79% on axial and lateral strain tensor images for visualization of thermal necrosis using the Multilevel approach leads to it being the best choice in a clinical setting. All methods, however, provide consistent and reproducible delineation for EDE in the ablation suite.

Dictionary Representations for Electrode Displacement Elastography

Ultrasound electrode displacement elastography (EDE) has demonstrated the potential to monitor ablated regions in human patients after minimally invasive microwave ablation procedures. Displacement estimation for EDE is commonly plagued by decorrelation noise artifacts degrading displacement estimates. In this paper, we propose a global dictionary learning approach applied to denoising displacement estimates with an adaptively learned dictionary from EDE phantom displacement maps. The resulting algorithm is one that represents displacement patches sparsely if they contain low noise and averages remaining patches thereby denoising displacement maps while retaining important edge information. The results of dictionary-represented displacements presented with a higher signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) with improved contrast, as well as improved phantom inclusion delineation when compared to initial displacements, median-filtered displacements, and spline smoothened displacements, respectively. In addition to visualized noise reduction, dictionary-represented displacements presented with the highest SNR, CNR, and improved contrast with values of 1.77, 4.56, and 4.35 dB, respectively, when compared to axial strain tensor images estimated using the initial displacements. Following EDE phantom imaging, we utilized dictionary representations from in vivo patient data, further validating efficacy. Denoising displacement estimates are a newer application for dictionary learning producing strong ablated region delineation with little degradation from denoising.

Delineation of Post-Procedure Ablation Regions with Electrode Displacement Elastography with a Comparison to Acoustic Radiation Force Impulse Imaging

We compared a quasi-static ultrasound elastography technique, referred to as electrode displacement elastography (EDE), with acoustic radiation force impulse imaging (ARFI) for monitoring microwave ablation (MWA) procedures on patients diagnosed with liver neoplasms. Forty-nine patients recruited to this study underwent EDE and ARFI with a Siemens Acuson S2000 system after an MWA procedure. On the basis of visualization results from two observers, the ablated region in ARFI images was recognizable on 20 patients on average in conjunction with B-mode imaging, whereas delineable ablation boundaries could be generated on 4 patients on average. With EDE, the ablated region was delineable on 40 patients on average, with less imaging depth dependence. Study of tissue-mimicking phantoms revealed that the ablation region dimensions measured on EDE and ARFI images were within 8%, whereas the image contrast and contrast-to-noise ratio with EDE was two to three times higher than that obtained with ARFI. This study indicated that EDE provided improved monitoring results for minimally invasive MWA in clinical procedures for liver cancer and metastases.

Monitoring Microwave Ablation of Ex Vivo Bovine Liver Using Ultrasonic Attenuation Imaging

Thermal ablation of soft tissue changes the tissue microstructure and, consequently, induces changes in its acoustic properties. Although B-mode ultrasound provides high-resolution and high-frame-rate images of ablative therapeutic procedures, it is not particularly effective at delineating boundaries of ablated regions because of poor contrast in echogenicity between ablated and surrounding normal tissue. Quantitative ultrasound techniques can provide quantitative estimates of acoustic properties, such as backscatter and attenuation coefficients, and differentiate ablated and unablated regions more effectively, with the potential for monitoring minimally invasive thermal therapies. In this study, a previously introduced attenuation estimation method was used to create quantitative attenuation coefficient maps for 11 microwave ablation procedures performed on refrigerated ex vivo bovine liver. The attenuation images correlate well with the pathologic images of the ablated region. The mean attenuation coefficient for regions of interest drawn inside and outside the ablated zones were 0.9 (±0.2) and 0.45 (±0.15) dB/cm/MHz, respectively. These estimates agree with reported values in the literature and establish the usefulness of non-invasive attenuation imaging for monitoring therapeutic procedures in the liver.

Comparison of three dimensional strain volume reconstructions using SOUPR and wobbler based acquisitions: A phantom study

PURPOSE

Ultrasound strain imaging is a relatively low cost and portable modality for monitoring percutaneous thermal ablation of liver neoplasms. However, a 3D strain volume reconstruction from existing 2D strain images is necessary to fully delineate the thermal dose distribution. Tissue mimicking (TM) phantom experiments were performed to validate a novel volume reconstruction algorithm referred to as sheaf of ultrasound planes reconstruction (SOUPR), based on a series of 2D rotational imaging planes.

METHODS

Reconstruction using SOUPR was formulated as an optimization problem with constraints on data consistency with 2D strain images and data smoothness of the volume data. Reconstructed ablation inclusion dimensions, volume, and elastographic signal to noise ratio (SNRe) and contrast to noise ratio (CNRe) were compared with conventional 3D ultrasound strain imaging based on interpolating a series of quasiparallel 2D strain images with a wobbler transducer.

RESULTS

Volume estimates of the phantom inclusion were in a similar range for both acquisition approaches. SNRe and CNRe obtained with SOUPR were significantly higher on the order of 250% and 166%, respectively. The mean error of the inclusion dimension reconstructed with a wobbler transducer was on the order of 10.4%, 3.5%, and 19.0% along the X, Y, and Z axes, respectively, while the error with SOUPR was on the order of 2.6%, 2.8%, and 9.6%. A qualitative comparison of SOUPR and wobbler reconstruction was also performed using a thermally ablated region created in ex vivo bovine liver tissue.

CONCLUSIONS

The authors have demonstrated using experimental evaluations with a TM phantom that the reconstruction results obtained with SOUPR were superior when compared with a conventional wobbler transducer in terms of inclusion shape preservation and detectability.

Post-Procedure Evaluation of Microwave Ablations of Hepatocellular Carcinomas Using Electrode Displacement Elastography

Microwave ablation has been used clinically as an alternative to surgical resection. However, lack of real-time imaging to assess treated regions may compromise treatment outcomes. We previously introduced electrode displacement elastography (EDE) for strain imaging and verified its feasibility in vivo on porcine animal models. In this study, we evaluated EDE on 44 patients diagnosed with hepatocellular carcinoma, treated using microwave ablation. The ablated region was identified on EDE images for 40 of the 44 patients. Ablation areas averaged 13.38 ± 4.99 cm² on EDE, compared with 7.61 ± 3.21 cm² on B-mode imaging. Contrast and contrast-to-noise ratios obtained with EDE were 232% and 98%, respectively, significantly higher than values measured on B-mode images (p < 0.001). This study indicates that EDE is feasible in patients and provides improved visualization of the ablation zone compared with B-mode ultrasound.

A coupled subsample displacement estimation method for ultrasound-based strain elastography

Obtaining accurate displacement estimates along both axial (parallel to the acoustic beam) and lateral (perpendicular to the beam) directions is an important task for several clinical applications such as shear strain imaging, modulus reconstruction and temperature imaging, where a full description of the two or three-dimensional (2D/3D) deformation field is required. In this study we propose an improved speckle tracking algorithm where axial and lateral motion estimations are simultaneously performed to enhance motion tracking accuracy. More specifically, using conventional ultrasound echo data, this algorithm first finds an iso-contour in the vicinity of the peak correlation between two segments of the pre- and post-deformation ultrasound radiofrequency echo data. The algorithm then attempts to find the center of the iso-contour of the correlation function that corresponds to the unknown (sub-sample) motion vector between these two segments of echo data. This algorithm has been tested using computer-simulated data, studies with a tissue-mimicking phantom, and in vivo breast lesion data. Computer simulation results show that the method improves the accuracy of both lateral and axial tracking. Such improvements are more significant when the deformation is small or along the lateral direction. Results from the tissue-mimicking phantom study are consistent with findings observed in computer simulations. Using in vivo breast lesion data we found that, compared to the 2D quadratic subsample displacement estimation methods, higher quality axial strain and shear strain images (e.g. 18.6% improvement in contrast-to-noise ratio for shear strain images) can be obtained for large deformations (up to 5% frame-to-frame and 15% local strains) in a multi-compression technique. Our initial results demonstrated that this conceptually and computationally simple method could improve the image quality of ultrasound-based strain elastography with current clinical equipment.

A PDE-Based Regularization Algorithm Toward Reducing Speckle Tracking Noise: A Feasibility Study for Ultrasound Breast Elastography

Obtaining accurate ultrasonically estimated displacements along both axial (parallel to the acoustic beam) and lateral (perpendicular to the beam) directions is an important task for various clinical elastography applications (e.g., modulus reconstruction and temperature imaging). In this study, a partial differential equation (PDE)-based regularization algorithm was proposed to enhance motion tracking accuracy. More specifically, the proposed PDE-based algorithm, utilizing two-dimensional (2D) displacement estimates from a conventional elastography system, attempted to iteratively reduce noise contained in the original displacement estimates by mathematical regularization. In this study, tissue incompressibility was the physical constraint used by the above-mentioned mathematical regularization. This proposed algorithm was tested using computer-simulated data, a tissue-mimicking phantom, and in vivo breast lesion data. Computer simulation results demonstrated that the method significantly improved the accuracy of lateral tracking (e.g., a factor of 17 at 0.5% compression). From in vivo breast lesion data investigated, we have found that, as compared with the conventional method, higher quality axial and lateral strain images (e.g., at least 78% improvements among the estimated contrast-to-noise ratios of lateral strain images) were obtained. Our initial results demonstrated that this conceptually and computationally simple method could be useful for improving the image quality of ultrasound elastography with current clinical equipment as a post-processing tool.

A survey of ultrasound elastography approaches to percutaneous ablation monitoring

Percutaneous thermal ablation has been widely used as a minimally invasive treatment for tumors. Treatment monitoring is essential for preventing complications while ensuring treatment efficacy. Mechanical testing measurements on tissue reveal that tissue stiffness increases with temperature and ablation duration. Different types of imaging methods can be used to monitor ablation procedures, including temperature or thermal strain imaging, strain imaging, modulus imaging, and shear modulus imaging. Ultrasound elastography demonstrates the potential to become the primary imaging modality for monitoring percutaneous ablation. This review briefly presented the state-of-the-art ultrasound elastography approaches for monitoring radiofrequency ablation and microwave ablation. These techniques were divided into four groups: quasi-static elastography, acoustic radiation force elastography, sonoelastography, and applicator motion elastography. Their advantages and limitations were compared and discussed. Future developments were proposed with respect to heat-induced bubbles, tissue inhomogeneities, respiratory motion, three-dimensional monitoring, multi-parametric monitoring, real-time monitoring, experimental data center for percutaneous ablation, and microwave ablation monitoring.

Three-dimensional sheaf of ultrasound planes reconstruction (SOUPR) of ablated volumes

This paper presents an algorithm for 3-D reconstruction of tumor ablations using ultrasound shear wave imaging with electrode vibration elastography. Radio-frequency ultrasound data frames are acquired over imaging planes that form a subset of a sheaf of planes sharing a common axis of intersection. Shear wave velocity is estimated separately on each imaging plane using a piecewise linear function fitting technique with a fast optimization routine. An interpolation algorithm then computes velocity maps on a fine grid over a set of C-planes that are perpendicular to the axis of the sheaf. A full 3-D rendering of the ablation can then be created from this stack of C-planes; hence the name "Sheaf Of Ultrasound Planes Reconstruction" or SOUPR. The algorithm is evaluated through numerical simulations and also using data acquired from a tissue mimicking phantom. Reconstruction quality is gauged using contrast and contrast-to-noise ratio measurements and changes in quality from using increasing number of planes in the sheaf are quantified. The highest contrast of 5 dB is seen between the stiffest and softest regions of the phantom. Under certain idealizing assumptions on the true shape of the ablation, good reconstruction quality while maintaining fast processing rate can be obtained with as few as six imaging planes suggesting that the method is suited for parsimonious data acquisitions with very few sparsely chosen imaging planes.

Considering angle selection when using ultrasound electrode displacement elastography to evaluate radiofrequency ablation of tissues

Percutaneous radiofrequency ablation (RFA) is a minimally invasive treatment to thermally destroy tumors. Ultrasound-based electrode-displacement elastography is an emerging technique for evaluating the region of RFA-induced lesions. The angle between the imaging probe and the RFA electrode can influence electrode-displacement elastography when visualizing the ablation zone. We explored the angle effect on electrode-displacement elastography to measure the ablation zone. Phantoms embedded with meatballs were fabricated and then ablated using an RFA system to simulate RFA-induced lesions. For each phantom, a commercial ultrasound scanner with a 7.5 MHz linear probe was used to acquire raw image data at different angles, ranging from 30° to 90° at increments of 10°, to construct electrode-displacement images and facilitate comparisons with tissue section images. The results revealed that the ablation regions detected using electrode-displacement elastography were highly correlated with those from tissue section images when the angle was between 30° and 60°. However, the boundaries of lesions were difficult to distinguish, when the angle was larger than 60°. The experimental findings suggest that angle selection should be considered to achieve reliable electrode-displacement elastography to describe ablation zones.

Mean scatterer spacing estimation in normal and thermally coagulated ex vivo bovine liver

The liver has been hypothesized to have a unique arrangement of microvasculature that presents as an arrangement of quasiperiodic scatterers to an interrogating ultrasound pulse. The mean scatterer spacing (MSS) of these quasiperiodic scatterers has been proposed as a useful quantitative ultrasound biomarker for characterizing liver tissue. Thermal ablation is an increasingly popular method for treating hepatic tumors, and ultrasonic imaging approaches for delineating the extent of thermal ablation are in high demand. In this work, we examine the distribution of estimated MSS in thermally coagulated bovine liver and normal untreated bovine liver ex vivo. We estimate MSS by detecting local maxima in the spectral coherence function of radio frequency echoes from a clinical transducer, the Siemens VFX 9L4 transducer operating on an S2000 scanner. We find that normal untreated bovine liver was characterized by an MSS of approximately 1.3 mm. We examined regions of interest 12 mm wide laterally, and ranging from 12 mm to 18 mm axially, in 2 mm increments. Over these parameters, the mode of the MSS estimates was between 1.25 and 1.37 mm. On the other hand, estimation of MSS in thermally coagulated liver tissue yields a distribution of MSS estimates whose mode varied between 0.45 and 1.0 mm when examining regions of interest over the same sizes. We demonstrate that the estimated MSS in thermally coagulated liver favors small spacings because the randomly positioned scatterers in this tissue are better modeled as aperiodic scatterers. The submillimeter spacings result from the fact that this was the most probable spacing to be estimated if the discretely sampled spectral coherence function was a uniformly random two-dimensional function.

Visualizing ex vivo radiofrequency and microwave ablation zones using electrode vibration elastography

PURPOSE

Electrode vibration elastography is a new shear wave imaging technique that can be used to visualize thermal ablation zones. Prior work has shown the ability of electrode vibration elastography to delineate radiofrequency ablations; however, there has been no previous study of delineation of microwave ablations or radiological-pathological correlations using multiple observers.

METHODS

Radiofrequency and microwave ablations were formed in ex vivo bovine liver tissue. Their visualization was compared on shear wave velocity and maximum displacement images. Ablation dimensions were compared to gross pathology. Elastographic imaging and gross pathology overlap and interobserver variability were quantified using similarity measures.

RESULTS

Elastographic imaging correlated with gross pathology. Correlation of area estimates was better in radiofrequency than in microwave ablations, with Pearson coefficients of 0.79 and 0.54 on shear wave velocity images and 0.90 and 0.70 on maximum displacement images for radiofrequency and microwave ablations, respectively. The absolute relative difference in area between elastographic imaging and gross pathology was 18.9% and 22.9% on shear wave velocity images and 16.0% and 23.1% on maximum displacement images for radiofrequency and microwave ablations, respectively.

CONCLUSIONS

Statistically significant radiological-pathological correlation was observed in this study, but correlation coefficients were lower than other modulus imaging techniques, most notably in microwave ablations. Observers provided similar delineations for most thermal ablations. These results suggest that electrode vibration elastography is capable of imaging thermal ablations, but refinement of the technique may be necessary before it can be used to monitor thermal ablation procedures clinically.

Characterizing the compression-dependent viscoelastic properties of human hepatic pathologies using dynamic compression testing

Recent advances in elastography have provided several imaging modalities capable of quantifying the elasticity of tissue, an intrinsic tissue property. This information is useful for determining tumour margins and may also be useful for diagnosing specific tumour types. In this study, we used dynamic compression testing to quantify the viscoelastic properties of 16 human hepatic primary and secondary malignancies and their corresponding background tissue obtained following surgical resection. Two additional backgrounds were also tested. An analysis of the background tissue showed that F4-graded fibrotic liver tissue was significantly stiffer than F0-graded tissue, with a modulus contrast of 4:1. Steatotic liver tissue was slightly stiffer than normal liver tissue, but not significantly so. The tumour-to-background storage modulus contrast of hepatocellular carcinomas, a primary tumour, was approximately 1:1, and the contrast decreased with increasing fibrosis grade of the background tissue. Ramp testing showed that the background stiffness increased faster than the malignant tissue. Conversely, secondary tumours were typically much stiffer than the surrounding background, with a tumour-to-background contrast of 10:1 for colon metastases and 10:1 for cholangiocarcinomas. Ramp testing showed that colon metastases stiffened faster than their corresponding backgrounds. These data have provided insights into the mechanical properties of specific tumour types, which may prove beneficial as the use of quantitative stiffness imaging increases.

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.

Bayesian regularization applied to ultrasound strain imaging

Noise artifacts due to signal decorrelation and reverberation are a considerable problem in ultrasound strain imaging. For block-matching methods, information from neighboring matching blocks has been utilized to regularize the estimated displacements. We apply a recursive Bayesian regularization algorithm developed by Hayton et al. [Artif. Intell., vol. 114, pp. 125-156, 1999] to phase-sensitive ultrasound RF signals to improve displacement estimation. The parameter of regularization is reformulated, and its meaning examined in the context of strain imaging. Tissue-mimicking experimental phantoms and RF data incorporating finite-element models for the tissue deformation and frequency-domain ultrasound simulations are used to compute the optimal parameter with respect to nominal strain and algorithmic iterations. The optimal strain regularization parameter was found to be twice the nominal strain and did not vary significantly with algorithmic iterations. The technique demonstrates superior performance over median filtering in noise reduction at strains 5% and higher for all quantitative experiments performed. For example, the strain SNR was 11 dB higher than that obtained using a median filter at 7% strain. It has to be noted that for applied deformations lower than 1%, since signal decorrelation errors are minimal, using this approach may degrade the displacement image.

Shear wave velocity imaging using transient electrode perturbation: phantom and ex vivo validation

This paper presents a new shear wave velocity imaging technique to monitor radio-frequency and microwave ablation procedures, coined electrode vibration elastography. A piezoelectric actuator attached to an ablation needle is transiently vibrated to generate shear waves that are tracked at high frame rates. The time-to-peak algorithm is used to reconstruct the shear wave velocity and thereby the shear modulus variations. The feasibility of electrode vibration elastography is demonstrated using finite element models and ultrasound simulations, tissue-mimicking phantoms simulating fully (phantom 1) and partially ablated (phantom 2) regions, and an ex vivo bovine liver ablation experiment. In phantom experiments, good boundary delineation was observed. Shear wave velocity estimates were within 7% of mechanical measurements in phantom 1 and within 17% in phantom 2. Good boundary delineation was also demonstrated in the ex vivo experiment. The shear wave velocity estimates inside the ablated region were higher than mechanical testing estimates, but estimates in the untreated tissue were within 20% of mechanical measurements. A comparison of electrode vibration elastography and electrode displacement elastography showed the complementary information that they can provide. Electrode vibration elastography shows promise as an imaging modality that provides ablation boundary delineation and quantitative information during ablation procedures.