Evaluation of the Accuracy of Liver Lesion DCEUS Quantification With Respiratory Gating

A Comprehensive and Repeatable Contrast-Enhanced Ultrasound Quantification Approach for Clinical Evaluations of Tumor Blood Flow

OBJECTIVE

The aim of this study is to define a comprehensive and repeatable contrast-enhanced ultrasound (CEUS) imaging protocol and analysis method to quantitatively assess lesional blood flow. Easily repeatable CEUS evaluations are essential for longitudinal treatment monitoring. The quantification method described here aims to provide a structure for future clinical studies.

MATERIALS AND METHODS

This retrospective analysis study included liver CEUS studies in 80 patients, 40 of which contained lesions (primarily hepatocellular carcinoma, n = 28). Each patient was given at least 2 injections of a microbubble contrast agent, and 60-second continuous loops were acquired for each injection to enable evaluation of repeatability. For each bolus injection, 1.2 mL of contrast was delivered, whereas continuous, stationary scanning was performed. Automated respiratory gating and motion compensation algorithms dealt with breathing motion. Similar in size regions of interest were drawn around the lesion and liver parenchyma, and time-intensity curves (TICs) with linearized image data were generated. Four bolus transit parameters, rise time ( RT ), mean transit time ( MTT ), peak intensity ( PI ), and area under the curve ( AUC ), were extracted either directly from the actual TIC data or from a lognormal distribution curve fitted to the TIC. Interinjection repeatability for each parameter was evaluated with coefficient of variation. A 95% confidence interval was calculated for all fitted lognormal distribution curve coefficient of determination ( R2 ) values, which serves as a data quality metric. One-sample t tests were performed between values obtained from injection pairs and between the fitted lognormal distribution curve and direct extraction from the TIC calculation methods to establish there were no significant differences between injections and measurement precision, respectively.

RESULTS

Average interinjection coefficient of variation with both the fitted curve and direct calculation of RT and MTT was less than 21%, whereas PI and AUC were less than 40% for lesion and parenchyma regions of interest. The 95% confidence interval for the R2 value of all fitted lognormal curves was [0.95, 0.96]. The 1-sample t test for interinjection value difference showed no significant differences, indicating there was no relationship between the order of the repeated bolus injections and the resulting parameters. The 1-sample t test between the values from the fitted lognormal distribution curve and the direct extraction from the TIC calculation found no statistically significant differences (α = 0.05) for all perfusion-related parameters except lesion and parenchyma PI and lesion MTT .

CONCLUSIONS

The scanning protocol and analysis method outlined and validated in this study provide easily repeatable quantitative evaluations of lesional blood flow with bolus transit parameters in CEUS data that were not available before. With vital features such as probe stabilization ideally performed with an articulated arm and an automated respiratory gating algorithm, we were able to achieve interinjection repeatability of blood flow parameters that are comparable or surpass levels currently established for clinical 2D CEUS scans. Similar values and interinjection repeatability were achieved between calculations from a fitted curve or directly from the data. This demonstrated not only the strength of the protocol to generate TICs with minimal noise, but also suggests that curve fitting might be avoided for a more standardized approach. Utilizing the imaging protocol and analysis method defined in this study, we aim for this methodology to potentially assist clinicians to assess true perfusion changes for treatment monitoring with CEUS in longitudinal studies.

Ultrasound fusion imaging for improving diagnostic and therapeutic strategies of focal liver lesions: A preliminary study

PURPOSE

To assess the effect of ultrasound (US) fusion imaging on the clinical diagnostic and therapeutic strategies of focal liver lesions, which are difficult to detect or diagnose by conventional US.

METHODS

From November 2019 to June 2022, 71 patients with invisible or undiagnosed focal liver lesions who underwent fusion imaging combining US with CT or MR were included in this retrospective study. The reasons for US fusion imaging were as follows: (1) lesions that were undetectable or inconspicuous on B-mode US; (2) post-ablation lesions that could not be assessed accurately by B-mode US; (3) to evaluate whether the lesions detected by B-mode US that were consistent with those presented on MRI/CT images.

RESULTS

Of the 71 cases, 43 cases were single lesions, and 28 cases were multiple lesions. Among the 46 cases which were invisible on conventional US, the display rate of lesions using US-CT/MRI fusion imaging was 30.8%, and that combined with CEUS was 76.9%. US-guided biopsy was performed in 30 patients after the detection and localization determined by fusion imaging, with a positive rate of 73.3%. Six patients with recurrence after ablation therapy were all detected and located accurately after fusion imaging, and 4 of them successfully underwent ablation therapy again.

CONCLUSION

Fusion imaging contributes to the understanding of the anatomical relationship between lesion location and blood vessels. Additionally, fusion imaging can improve the diagnostic confidence, be helpful to guide interventional operations, and hence be conducive to clinical therapeutic strategies.

Irregular Respiratory Motion Compensation for Liver Contrast-Enhanced Ultrasound via Transport-Based Motion Estimation

Contrast-enhanced ultrasound (CEUS) imaging has been widely applied for the detection and characterization of focal liver lesions (FLLs), for its ability to visualize the blood flow in real time. However, cyclic liver motion poses a great challenge to the recovery of perfusion curves as well as quantitative kinetic parameters estimation. Recently, a few gating methods have been proposed to eliminate unexpected intensity fluctuations by the breathing phase estimation, with the assumption that each breathing phase corresponds to a specific lesion position strictly. While practical liver motion tends to be irregular due to changes in the patient's underlying physiologic status, that is, the same phase might not correspond to the same position. To tackle the challenge of motion irregularity, we present a novel motion estimation-based respiratory compensation method, named RCME, which first estimates salient motion information through the framework of optimal transport (OT) by jointly modeling pixel intensity as well as their locations and then employs sparse subspace clustering (SSC) to identify the subset of frames acquired at the same position. Our proposed method is evaluated on 15 clinical CEUS sequences in both qualitative and quantitative manners. Experimental results demonstrate good performance on irregular liver motion compensation.

Respiratory motion correction for liver contrast-enhanced ultrasound by automatic selection of a reference image

PURPOSE

Respiratory motion correction is necessary for the quantitative analysis of liver contrast-enhanced ultrasound (CEUS) image sequences. Most respiratory motion correction methods are based on the dual mode of CEUS image sequences, including contrast and grayscale image sequences. Due to free-breathing motion, the acquired two-dimensional (2D) ultrasound cine might show the in-plane and out-of-plane motion of tumors. The registration of an entire 2D ultrasound contrast image sequence based on out-of-plane images is ineffective. For the respiratory motion correction of CEUS sequences, the reference image is usually considered the standard for the deletion of any out-of-plane images. Most methods used for the selection of the reference image are subjective in nature. Here, a quantitative selection method for an optimal reference image from CEUS image sequences in the B mode and contrast mode was explored.

METHODS

The original high-dimensional ultrasound grayscale image data were mapped into a two-dimensional space using Laplacian Eigenmaps (LE), and K-means clustering was adopted. The center image of the larger cluster with a near-peak contrast intensity was considered the optimal ultrasound reference image. In the ultrasound grayscale image sequence, the images with the maximum correlations to the reference image in the same time interval were selected as the corrected image sequence. The effectiveness of this proposed method was then validated on 18 CEUS cases of VX2 tumors in rabbit livers.

RESULTS

Correction smoothed the time-intensity curves (TICs) extracted from the region of interest of the CEUS image sequences. Before correction, the average of the total mean structural similarity (TMSSIM) and the average of the mean correlation coefficient (MCC) from the image sequences were 0.45 ± 0.11 and 0.67 ± 0.16, respectively, and after correction, the average TMSSIM and MCC increased (P < 0.001) by 31% to 0.59 ± 0.11 and by 21% to 0.81 ± 0.11, respectively. The average deviation value (DV) index of the TICs from the image sequences prior to correction was 92.16 ± 18.12, and correction reduced the average to 31.71 ± 7.31. The average TMSSIM and MCC values after correction using the mean frame of the reference image (MBMFRI) were clearly lower than those after correction using the proposed method (P < 0.001). Moreover, the average DV after correction using the MBMFRI was obviously higher than that after correction using the proposed method (P < 0.001).

CONCLUSIONS

The breathing frequency of rabbits is notably faster than that of human beings, but the proposed correction method could reduce the effect of the respiratory motion in the CEUS image sequences. The reference image was selected quantitatively, which could improve the accuracy of the quantitative analysis of rabbit liver CEUS sequences using the reference image method based on the current standard of manual selection and the MBMFRI. This easy-to-operate method can potentially be used in both animal studies and clinical applications.

Automatic Respiratory Gating Hepatic DCEUS-based Dual-phase Multi-parametric Functional Perfusion Imaging using a Derivative Principal Component Analysis

Purpose: Angiogenesis in liver cancers can be characterized by hepatic functional perfusion imaging (FPI) on the basis of dynamic contrast-enhanced ultrasound (DCEUS). However, accuracy is limited by breathing motion which results in out-of-plane image artifacts. Current hepatic FPI studies do not correct for these artifacts and lack the evaluation of correction accuracy. Thus, a hepatic DCEUS-based dual-phase multi-parametric FPI (DM-FPI) scheme using a derivative principal component analysis (PCA) respiratory gating is proposed to overcome these limitations.

Materials and Methods: By considering severe 3D out-of-plane respiratory motions, the proposed scheme's accuracy was verified with in vitro DCEUS experiments in a flow model mimicking a hepatic vein. The feasibility was further demonstrated by considering in vivo DCEUS measurements in normal rabbit livers, and hepatic cavernous hemangioma and hepatocellular carcinoma in patients. After respiratory kinetics was extracted through PCA of DCEUS sequences under free-breathing condition, dual-phase respiratory gating microbubble kinetics was identified by using a derivative PCA zero-crossing dual-phase detection, respectively. Six dual-phase hemodynamic parameters were estimated from the dual-phase microbubble kinetics and DM-FPI was then reconstructed via color-coding to quantify 2.5D angiogenic hemodynamic distribution for live tumors.

Results: Compared with no respiratory gating, the mean square error of respiratory gating DM-FPI decreased by 1893.9 ± 965.4 (p < 0.05), and mean noise coefficients decreased by 17.5 ± 7.1 (p < 0.05), whereas correlation coefficients improved by 0.4 ± 0.2 (p < 0.01). DM-FPI observably removed severe respiratory motion artifacts on PFI and markedly enhanced the accuracy and robustness both in vitro and in vivo.

Conclusions: DM-FPI precisely characterized and distinguished the heterogeneous angiogenic hemodynamics about perfusion volume, blood flow and flow rate within two anatomical sections in the normal liver, and in benign and malignant hepatic tumors. DCEUS-based DM-FPI scheme might be a useful tool to help clinicians diagnose and provide suitable therapies for liver tumors.

In-vitro evaluation of accuracy of dynamic contrast-enhanced ultrasound (DCEUS)-based parametric perfusion imaging with respiratory motion-compensation

PURPOSE

The accuracy of multi-parametric perfusion imaging (PPI) based on dynamic contrast-enhanced ultrasound is disturbed by the respiratory motion in some cases, especially during characterizing hemodynamic features of abdominal tumor angiogenesis. This study aimed to effectively remove those disturbances on PPI and evaluate its accuracy.

METHOD

The respiratory motion-compensation (rMoCo) strategy in PPI was modified by employing non-negative matrix factorization combined with phase-by-phase compensation. According to the known and controllable ground truths in in-vitro perfusion experiments, the accuracy of the modified rMoCo strategy was further evaluated from multiple perspectives in a simulated dual-vessel flow phantom.

RESULTS

Compared with that of PPIs without rMoCo, the mean correlation coefficient between six PPIs with rMoCo and the corresponding static PPIs was up to 0.98 ± 0.01 and improved by 0.17 ± 0.04 (P < 0.05). The estimated error of vascular diameter decreased from 87.85% (P < 0.05) to 7.25% (P < 0.05) after rMoCo. PPIs with rMoCo were significantly consistent with static PPIs without respiratory motion disturbances.

CONCLUSIONS

These quantitative results illustrated the disturbances induced by respiratory motion were effectively removed and the accuracy of PPIs was significantly improved. The partial parabolic and bimodal hemodynamic characteristics and the anatomical structures and sizes were accurately quantified and depicted by PPIs with rMoCo. The modified method can benefit physicians in providing accurate diagnoses and in developing appropriate therapeutic strategies for abdominal diseases.

Automatic Identification of the Optimal Reference Frame for Segmentation and Quantification of Focal Liver Lesions in Contrast-Enhanced Ultrasound

Post-examination interpretation of contrast-enhanced ultrasound (CEUS) cineloops of focal liver lesions (FLLs) requires offline manual assessment by experienced radiologists, which is time-consuming and generates subjective results. Such assessment usually starts by manually identifying a reference frame, where FLL and healthy parenchyma are well-distinguished. This study proposes an automatic computational method to objectively identify the optimal reference frame for distinguishing and hence delineating an FLL, by statistically analyzing the temporal intensity variation across the spatially discretized ultrasonographic image. Level of confidence and clinical value of the proposed method were quantitatively evaluated on retrospective multi-institutional data (n = 64) and compared with expert interpretations. Results support the proposed method for facilitating easier, quicker and reproducible assessment of FLLs, further increasing the radiologists' confidence in diagnostic decisions. Finally, our method yields a useful training tool for radiologists, widening CEUS use in non-specialist centers, potentially leading to reduced turnaround times and lower patient anxiety and healthcare costs.