Towards clinical prostate ultrasound elastography using full inversion approach

On the soft tissue ultrasound elastography using FEM based inversion approach

Elastography is a medical imaging modality that enables visualization of tissue stiffness. It involves quasi-static or harmonic mechanical stimulation of the tissue to generate a displacement field which is used as input in an inversion algorithm to reconstruct tissue elastic modulus. This paper considers quasi-static stimulation and presents a novel inversion technique for elastic modulus reconstruction. The technique follows an inverse finite element framework. Reconstructed elastic modulus maps produced in this technique do not depend on the initial guess, while it is computationally less involved than iterative reconstruction approaches. The method was first evaluated using simulated data (in-silico) where modulus reconstruction's sensitivity to displacement noise and elastic modulus was assessed. To demonstrate the method's performance, displacement fields of two tissue mimicking phantoms determined using three different motion tracking techniques were used as input to the developed elastography method to reconstruct the distribution of relative elastic modulus of the inclusion to background tissue. In the next stage, the relative elastic modulus of three clinical cases pertaining to liver cancer patient were determined. The obtained results demonstrate reasonably high elastic modulus reconstruction accuracy in comparison with similar direct methods. Also it is associated with reduced computational cost in comparison with iterative techniques, which suffer from convergence and uniqueness issues, following the same formulation concept. Moreover, in comparison with other methods which need initial guess, the presented method does not require initial guess while it is easy to understand and implement.

Enhanced full-inversion-based ultrasound elastography for evaluating tumor response to neoadjuvant chemotherapy in patients with locally advanced breast cancer

PURPOSE

An enhanced ultrasound elastography technique is proposed for early assessment of locally advanced breast cancer (LABC) response to neoadjuvant chemotherapy (NAC).

METHODS

The proposed elastography technique inputs ultrasound radiofrequency data obtained through tissue quasi-static stimulation and adapts a strain refinement algorithm formulated based on fundamental principles of continuum mechanics, coupled with an iterative inverse finite element method to reconstruct the breast Young's modulus (E) images. The technique was explored for therapy response assessment using data acquired from 25 LABC patients before and at weeks 1, 2, and 4 after the NAC initiation (100 scans). The E ratio of tumor to the surrounding tissue was calculated at different scans and compared to the baseline for each patient. Patients' response to NAC was determined many months later using standard clinical and histopathological criteria.

RESULTS

Reconstructed E ratio changes obtained as early as one week after the NAC onset demonstrate very good separation between the two cohorts of responders and non-responders to NAC. Statistically significant differences were observed in the E ratio changes between the two patient cohorts at weeks 1 to 4 after treatment (p-value < 0.001; statistical power greater than 97%). A significant difference in axial strain ratio changes was observed only at week 4 (p-value = 0.01; statistical power = 76%). No significant difference was observed in tumor size changes at weeks 1, 2 or 4.

CONCLUSION

The proposed elastography technique demonstrates a high potential for chemotherapy response monitoring in LABC patients and superior performance compared to strain imaging.

Second-Order Ultrasound Elastography With L1-Norm Spatial Regularization

Time delay estimation (TDE) between two radio-frequency (RF) frames is one of the major steps of quasi-static ultrasound elastography, which detects tissue pathology by estimating its mechanical properties. Regularized optimization-based techniques, a prominent class of TDE algorithms, optimize a nonlinear energy functional consisting of data constancy and spatial continuity constraints to obtain the displacement and strain maps between the time-series frames under consideration. The existing optimization-based TDE methods often consider the L2 -norm of displacement derivatives to construct the regularizer. However, such a formulation over-penalizes the displacement irregularity and poses two major issues to the estimated strain field. First, the boundaries between different tissues are blurred. Second, the visual contrast between the target and the background is suboptimal. To resolve these issues, herein, we propose a novel TDE algorithm where instead of L2 -, L1 -norms of both first- and second-order displacement derivatives are taken into account to devise the continuity functional. We handle the non-differentiability of L1 -norm by smoothing the absolute value function's sharp corner and optimize the resulting cost function in an iterative manner. We call our technique Second-Order Ultrasound eLastography (SOUL) with the L1 -norm spatial regularization ( L1 -SOUL). In terms of both sharpness and visual contrast, L1 -SOUL substantially outperforms GLobal Ultrasound Elastography (GLUE), tOtal Variation rEgulaRization and WINDow-based time delay estimation (OVERWIND), and SOUL, three recently published TDE algorithms in all validation experiments performed in this study. In cases of simulated, phantom, and in vivo datasets, respectively, L1 -SOUL achieves 67.8%, 46.81%, and 117.35% improvements of contrast-to-noise ratio (CNR) over SOUL. The L1 -SOUL code can be downloaded from http://code.sonography.ai.

Combining First- and Second-Order Continuity Constraints in Ultrasound Elastography

Ultrasound elastography is a prominent noninvasive medical imaging technique that estimates tissue elastic properties to detect abnormalities in an organ. A common approximation to tissue elastic modulus is tissue strain induced after mechanical stimulation. To compute tissue strain, ultrasound radio frequency (RF) data can be processed using energy-based algorithms. These algorithms suffer from ill-posedness to tackle. A continuity constraint along with the data amplitude similarity is imposed to obtain a unique solution to the time-delay estimation (TDE) problem. Existing energy-based methods exploit the first-order spatial derivative of the displacement field to construct a regularizer. This first-order regularization scheme alone is not fully consistent with the mechanics of tissue deformation while perturbed with an external force. As a consequence, state-of-the-art techniques suffer from two crucial drawbacks. First, the strain map is not sufficiently smooth in uniform tissue regions. Second, the edges of the hard or soft inclusions are not well-defined in the image. Herein, we address these issues by formulating a novel regularizer taking both first- and second-order derivatives of the displacement field into account. The second-order constraint, which is the principal novelty of this work, contributes both to background continuity and edge sharpness by suppressing spurious noisy edges and enhancing strong boundaries. We name the proposed technique: Second-Order Ultrasound eLastography (SOUL). Comparative assessment of qualitative and quantitative results shows that SOUL substantially outperforms three recently developed TDE algorithms called Hybrid, GLUE, and MPWC-Net++. SOUL yields 27.72%, 62.56%, and 81.37% improvements of the signal-to-noise ratio (SNR) and 72.35%, 54.03%, and 65.17% improvements of the contrast-to-noise ratio (CNR) over GLUE with data pertaining to simulation, phantom, and in vivo tissue, respectively. The SOUL code can be downloaded from code.sonography.ai.

Supervised Classification of the Accuracy of the Time Delay Estimation in Ultrasound Elastography

The accuracy of time-delay estimation (TDE) in ultrasound elastography is usually measured by calculating the value of normalized cross correlation (NCC) at the estimated displacement. NCC value, however, could be very high at a displacement estimate with large error, a well-known problem in TDE referred to as peak-hopping. Furthermore, NCC value could suffer from jitter error, which is due to electric noise and signal decorrelation. Herein, we propose a novel method to assess the accuracy of TDE by investigating the NCC profile around the estimated time-delay. We extract several features from the NCC profile, and utilize support vector machine to classify peak-hopping and jitter error. The results on simulation, phantom, and in vivo data show the significant improvement of the proposed algorithm compared to the state of the art techniques.

The Role of Affordable, Point-of-Care Technologies for Cancer Care in Low- and Middle-Income Countries: A Review and Commentary

As the burden of non-communicable diseases such as cancer continues to rise in low- and middle-income countries (LMICs), it is essential to identify and invest in promising solutions for cancer control and treatment. Point-of-care technologies (POCTs) have played critical roles in curbing infectious disease epidemics in both high- and low-income settings, and their successes can serve as a model for transforming cancer care in LMICs, where access to traditional clinical resources is often limited. The versatility, cost-effectiveness, and simplicity of POCTs warrant attention for their potential to revolutionize cancer detection, diagnosis, and treatment. This paper reviews the landscape of affordable POCTs for cancer care in LMICs with a focus on imaging tools, in vitro diagnostics, and treatment technologies and aspires to encourage innovation and further investment in this space.

Ultrasound Elastography of the Prostate Using an Unconstrained Modulus Reconstruction Technique: A Pilot Clinical Study

A novel full-inversion-based technique for quantitative ultrasound elastography was investigated in a pilot clinical study on five patients for non-invasive detection and localization of prostate cancer and quantification of its extent. Conventional-frequency ultrasound images and radiofrequency (RF) data (~5 MHz) were collected during mechanical stimulation of the prostate using a transrectal ultrasound probe. Pre and post-compression RF data were used to construct the strain images. The Young's modulus (YM) images were subsequently reconstructed using the derived strain images and the stress distribution estimated iteratively using finite element (FE) analysis. Tumor regions determined based on the reconstructed YM images were compared to whole-mount histopathology images of radical prostatectomy specimens. Results indicated that tumors were significantly stiffer than the surrounding tissue, demonstrating a relative YM of 2.5 ± 0.8 compared to normal prostate tissue. The YM images had a good agreement with the histopathology images in terms of tumor location within the prostate. On average, 76% ± 28% of tumor regions detected based on the proposed method were inside respective tumor areas identified in the histopathology images. Results of a linear regression analysis demonstrated a good correlation between the disease extents estimated using the reconstructed YM images and those determined from whole-mount histopathology images (r² = 0.71). This pilot study demonstrates that the proposed method has a good potential for detection, localization and quantification of prostate cancer. The method can potentially be used for prostate needle biopsy guidance with the aim of decreasing the number of needle biopsies. The proposed technique utilizes conventional ultrasound imaging system only while no additional hardware attachment is required for mechanical stimulation or data acquisition. Therefore, the technique may be regarded as a non-invasive, low cost and potentially widely-available clinical tool for prostate cancer diagnosis.

Estimation of Strain Elastography from Ultrasound Radio-Frequency Data by Utilizing Analytic Gradient of the Similarity Metric

Most strain imaging techniques follow a pipeline strategy: in the first step, tissue displacement is estimated from radio-frequency (RF) frames, and in the second step, a spatial derivative operation is applied. There are two main issues that arise from this framework. First, the gradient operation amplifies noise, and therefore, smoothing techniques have to be adopted. Second, strain estimation does not exploit the original RF data. It rather relies solely on the noisy displacement field. In this paper, a novel technique is proposed that utilizes both the displacement field and the RF frames to accurately obtain the strain estimates. The normalized cross correlation (NCC) metric between two corresponding windows around the samples of the pre- and post-compressed images is employed to generate a dissimilarity measurement. The derivative of NCC with respect to the strain is analytically derived using the chain rule. This allows an efficient minimization of the dissimilarity metric with respect to the strain using the gradient descent optimization technique. The effectiveness of the proposed method is investigated through simulation data, phantom experiments, and in vivo patient data. The experimental results show that exploiting the information in RF data significantly improves the strain estimates.

Relative Elastic Modulus Imaging Using Sector Ultrasound Data for Abdominal Applications: An Evaluation of Strategies and Feasibility

We reconstruct the elastic modulus distribution for one tissue mimicking (TM) phantom and two in vivo biopsy-confirmed liver tumors using curvilinear ultrasound echo data. Spatial distribution of the relative elastic modulus values is determined by solving an inverse problem within a region of interest (ROI). This inverse problem solution requires knowledge of the ultrasonically measured displacement field in a uniform rectilinear grid to ensure that the resolution on the resultant relative elastic modulus elastogram will be uniform over the entire ROI. Taking advantage of a new speckle tracking algorithm, two different displacement tracking strategies are investigated: 1) sector-shaped ultrasound data were converted to ultrasound data on a rectilinear grid prior to speckle tracking and 2) axial and lateral displacements directly obtained from sector-shaped data were converted to vertical and horizontal displacements on a rectilinear grid after speckle tracking. Compared with strain elastography (SE), TM phantom results show that relative elastic modulus imaging (REMI) using Strategy 2 provided higher contrast-to-noise ratios (>300% and 25% increases compared with SE and REMI using Strategy 1, respectively). Furthermore, in phantoms, REMI using Strategy 2 more accurately (a 1.3% difference to shear wave elastography measurements) estimated the elastic contrast ratio between the target and the background, compared with both SE (20%-25%) and REMI using Strategy 1 (4.1%). It was also observed that relative modulus elastograms were more consistent with anatomical structures visualized on corresponding B-mode images for the two in vivo liver cases. Overall, we conclude that applying REMI is feasible for abdominal organs such as the liver. Strategy 2 offered improved and consistent results for the data investigated.

A regularization-free elasticity reconstruction method for ultrasound elastography with freehand scan

BACKGROUND

In ultrasound elastography, reconstruction of tissue elasticity (e.g., Young's modulus) requires regularization and known information of forces and/or displacements on tissue boundaries. In practice, it is challenging to choose an appropriate regularization parameter; and the boundary conditions are difficult to obtain in vivo. The purpose of this study is to develop a more applicable algorithm that does not need any regularization or boundary force/displacement information.

METHODS

The proposed method adopts the bicubic B-spline as the tissue motion model to estimate the displacement fields. Then the estimated displacements are input to the finite element inversion scheme to reconstruct the Young's modulus of each element. In the inversion, a modulus boundary condition is used instead of force/displacement boundary conditions. Simulation and experiments on tissue-mimicking phantoms are carried out to test the proposed method.

RESULTS

The simulation results demonstrate that Young's modulus reconstruction of the proposed method has a relative error of -3.43 ± 0.43% and root-squared-mean error of 16.94 ± 0.25%. The phantom experimental results show that the target hardening artifacts in the strain images are significantly reduced in the Young's modulus images. In both simulation and phantom studies, the size and position of inclusions can be accurately depicted in the modulus images.

CONCLUSIONS

The proposed method can reconstruct tissue Young's modulus distribution with a high accuracy. It can reduce the artifacts shown in the strain image and correctly delineate the locations and sizes of inclusions. Unlike most modulus reconstruction methods, it does not need any regularization during the inversion procedure. Furthermore, it does not need to measure the boundary conditions of displacement or force. Thus this method can be used with a freehand scan, which facilitates its usage in the clinic.