Comparing a Portable Motion Analysis System against the Gold Standard for Potential Anterior Cruciate Ligament Injury Prevention and Screening.

Knee kinematics during a drop vertical jump, measured by the Kinect V2 (Microsoft, Redmond, WA, USA), have been shown to be associated with an increased risk of non-contact anterior cruciate ligament injury. The accuracy and reliability of the Microsoft Kinect V2 has yet to be assessed specifically for tracking the coronal and sagittal knee angles of the drop vertical jump. Eleven participants performed three drop vertical jumps that were recorded using both the Kinect V2 and a gold standard motion analysis system (Vicon, Los Angeles, CA, USA). The initial coronal, peak coronal, and peak sagittal angles of the left and right knees were measured by both systems simultaneously. Analysis of the data obtained by the Kinect V2 was performed by our software. The differences in the mean knee angles measured by the Kinect V2 and the Vicon system were non-significant for all parameters except for the peak sagittal angle of the right leg with a difference of 7.74 degrees and a p-value of 0.008. There was excellent agreement between the Kinect V2 and the Vicon system, with intraclass correlation coefficients consistently over 0.75 for all knee angles measured. Visual analysis revealed a moderate frame-to-frame variability for coronal angles measured by the Kinect V2. The Kinect V2 can be used to capture knee coronal and sagittal angles with sufficient accuracy during a drop vertical jump, suggesting that a Kinect-based portable motion analysis system is suitable to screen individuals for the risk of non-contact anterior cruciate ligament injury.

Clutter suppression in ultrasound: performance evaluation and review of low-rank and sparse matrix decomposition methods.

Vessel diseases are often accompanied by abnormalities related to vascular shape and size. Therefore, a clear visualization of vasculature is of high clinical significance. Ultrasound color flow imaging (CFI) is one of the prominent techniques for flow visualization. However, clutter signals originating from slow-moving tissue are one of the main obstacles to obtain a clear view of the vascular network. Enhancement of the vasculature by suppressing the clutters is a significant and irreplaceable step for many applications of ultrasound CFI. Currently, this task is often performed by singular value decomposition (SVD) of the data matrix. This approach exhibits two well-known limitations. First, the performance of SVD is sensitive to the proper manual selection of the ranks corresponding to clutter and blood subspaces. Second, SVD is prone to failure in the presence of large random noise in the dataset. A potential solution to these issues is using decomposition into low-rank and sparse matrices (DLSM) framework. SVD is one of the algorithms for solving the minimization problem under the DLSM framework. Many other algorithms under DLSM avoid full SVD and use approximated SVD or SVD-free ideas which may have better performance with higher robustness and less computing time. In practice, these models separate blood from clutter based on the assumption that steady clutter represents a low-rank structure and that the moving blood component is sparse. In this paper, we present a comprehensive review of ultrasound clutter suppression techniques and exploit the feasibility of low-rank and sparse decomposition schemes in ultrasound clutter suppression. We conduct this review study by adapting 106 DLSM algorithms and validating them against simulation, phantom, and in vivo rat datasets. Two conventional quality metrics, signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR), are used for performance evaluation. In addition, computation times required by different algorithms for generating clutter suppressed images are reported. Our extensive analysis shows that the DLSM framework can be successfully applied to ultrasound clutter suppression.

The effect of low back pain and lower limb injury on lumbar multifidus muscle morphology and function in university soccer players.

BACKGROUND: The lumbar multifidus muscle (LMM) plays a critical role to stabilize the spine. While low back pain (LBP) is a common complaint in soccer players, few studies have examined LMM characteristics in this athletic population and their possible associations with LBP and lower limb injury. Therefore, the purpose of this study was to 1) investigate LMM characteristics in university soccer players and their potential association with LBP and lower limb injury; 2) examine the relationship between LMM characteristics and body composition measurements; and 3) examine seasonal changes in LMM characteristics. METHODS: LMM ultrasound assessments were acquired in 27 soccer players (12 females, 15 males) from Concordia University during the preseason and assessments were repeated in 18 players at the end of the season. LMM cross-sectional area (CSA), echo-intensity and thickness at rest and during contraction (e.g. function) were assessed bilaterally in prone and standing positions, at the L5-S1 spinal level. A self-reported questionnaire was used to assess the history of LBP and lower limb injury. Dual-energy x-ray absorptiometry (DEXA) was used to acquire body composition measurements. RESULTS: Side-to-side asymmetry of the LMM was significantly greater in males (p = 0.02). LMM thickness when contracted in the prone position (p = 0.04) and LMM CSA in standing (p = 0.02) were also significantly greater on the left side in male players. The LMM % thickness change during contraction in the prone position was significantly greater in players who reported having LBP in the previous 3-months (p < 0.001). LMM CSA (r = - 0.41, p = 0.01) and echo-intensity (r = 0.69, p < 0.001) were positively correlated to total % body fat. There was a small decrease in LMM thickness at rest in the prone position over the course of the season (p = 0.03). CONCLUSIONS: The greater LMM contraction in players with LBP may be a maladaptive strategy to splint and project the spine. LMM morphology measurements were correlated to body composition. The results provide new insights with regards to LMM morphology and activation in soccer players and their associations with injury and body composition measurements.

LUMINOUS database: lumbar multifidus muscle segmentation from ultrasound images.

BACKGROUND: Among the paraspinal muscles, the structure and function of the lumbar multifidus (LM) has become of great interest to researchers and clinicians involved in lower back pain and muscle rehabilitation. Ultrasound (US) imaging of the LM muscle is a useful clinical tool which can be used in the assessment of muscle morphology and function. US is widely used due to its portability, cost-effectiveness, and ease-of-use. In order to assess muscle function, quantitative information of the LM must be extracted from the US image by means of manual segmentation. However, manual segmentation requires a higher level of training and experience and is characterized by a level of difficulty and subjectivity associated with image interpretation. Thus, the development of automated segmentation methods is warranted and would strongly benefit clinicians and researchers. The aim of this study is to provide a database which will contribute to the development of automated segmentation algorithms of the LM. CONSTRUCTION AND CONTENT: This database provides the US ground truth of the left and right LM muscles at the L5 level (in prone and standing positions) of 109 young athletic adults involved in Concordia University's varsity teams. The LUMINOUS database contains the US images with their corresponding manually segmented binary masks, serving as the ground truth. The purpose of the database is to enable development and validation of deep learning algorithms used for automatic segmentation tasks related to the assessment of the LM cross-sectional area (CSA) and echo intensity (EI). The LUMINOUS database is publicly available at http://data.sonography.ai . CONCLUSION: The development of automated segmentation algorithms based on this database will promote the standardization of LM measurements and facilitate comparison among studies. Moreover, it can accelerate the clinical implementation of quantitative muscle assessment in clinical and research settings.

Population-averaged MRI atlases for automated image processing and assessments of lumbar paraspinal muscles.

PURPOSE: Growing evidence suggests an association between lumbar paraspinal muscle degeneration and low back pain (LBP). Currently, time-consuming and laborious manual segmentations of paraspinal muscles are commonly performed on magnetic resonance imaging (MRI) axial scans. Automated image analysis algorithms can mitigate these drawbacks, but they often require individual MRIs to be aligned to a standard "reference" atlas. Such atlases are well established in automated neuroimaging analysis. Our aim was to create atlases of similar nature for automated paraspinal muscle measurements. METHODS: Lumbosacral T2-weighted MRIs were acquired from 117 patients who experienced LBP, stratified by gender and age group (30-39, 40-49, and 50-59 years old). Axial MRI slices of the L4-L5 and L5-S1 levels at mid-disc were obtained and aligned using group-wise linear and nonlinear image registration to produce a set of unbiased population-averaged atlases for lumbar paraspinal muscles. RESULTS: The resulting atlases represent the averaged morphology and MRI intensity features of the corresponding cohorts. Differences in paraspinal muscle shapes and fat infiltration levels with respect to gender and age can be visually identified from the population-averaged data from both linear and nonlinear registrations. CONCLUSION: We constructed a set of population-averaged atlases for developing automated algorithms to help analyze paraspinal muscle morphometry from axial MRI scans. Such an advancement could greatly benefit the fields of paraspinal muscle and LBP research. These slides can be retrieved under Electronic Supplementary Material.

Evaluation of an automated thresholding algorithm for the quantification of paraspinal muscle composition from MRI images.

BACKGROUND: The imaging assessment of paraspinal muscle morphology and fatty infiltration has gained considerable attention in the past decades, with reports suggesting an association between muscle degenerative changes and low back pain (LBP). To date, qualitative and quantitative approaches have been used to assess paraspinal muscle composition. Though highly reliable, manual thresholding techniques are time consuming and not always feasible in a clinical setting. The tedious and rater-dependent nature of such manual thresholding techniques provides the impetus for the development of automated or semi-automated segmentation methods. The purpose of the present study was to develop and evaluate an automated thresholding algorithm for the assessment of paraspinal muscle composition. The reliability and validity of the muscle measurements using the new automated thresholding algorithm were investigated through repeated measurements and comparison with measurements from an established, highly reliable manual thresholding technique. METHODS: Magnetic resonance images of 30 patients with LBP were randomly selected cohort of patients participating in a project on commonly diagnosed lumbar pathologies in patients attending spine surgeon clinics. A series of T2-weighted MR images were used to train the algorithm; preprocessing techniques including adaptive histogram equalization method image adjustment scheme were used to enhance the quality and contrast of the images. All muscle measurements were repeated twice using a manual thresholding technique and the novel automated thresholding algorithm, from axial T2-weigthed images, at least 5 days apart. The rater was blinded to all earlier measurements. Inter-method agreement and intra-rater reliability for each measurement method were assessed. The study did not received external funding and the authors have no disclosures. RESULTS: There was excellent agreement between the two methods with inter-method reliability coefficients (intraclass correlation coefficients) varying from 0.79 to 0.99. Bland and Altman plots further confirmed the agreement between the two methods. Intra-rater reliability and standard error of measurements were comparable between methods, with reliability coefficient varying between 0.95 and 0.99 for the manual thresholding and 0.97-0.99 for the automated algorithm. CONCLUSION: The proposed automated thresholding algorithm to assess paraspinal muscle size and composition measurements was highly reliable, with excellent agreement with the reference manual thresholding method.

Alternating direction method of multipliers for displacement estimation in ultrasound strain elastography.

BACKGROUND: Ultrasound strain imaging, which delineates mechanical properties to detect tissue abnormalities, involves estimating the time delay between two radio-frequency (RF) frames collected before and after tissue deformation. The existing regularized optimization-based time-delay estimation (TDE) techniques suffer from at least one of the following drawbacks: (1) The regularizer is not aligned with the tissue deformation physics due to taking only the first-order displacement derivative into account; (2) The L 2 $L2$ -norm of the displacement derivatives, which oversmooths the estimated time-delay, is utilized as the regularizer; (3) The modulus function defined mathematically should be approximated by a smooth function to facilitate the optimization of L 1 $L1$ -norm. PURPOSE: Our purpose is to develop a novel TDE technique that resolves the aforementioned shortcomings of the existing algorithms. METHODS: Herein, we propose employing the alternating direction method of multipliers (ADMM) for optimizing a novel cost function consisting of L 2 $L2$ -norm data fidelity term and L 1 $L1$ -norm first- and second-order spatial continuity terms. ADMM empowers the proposed algorithm to use different techniques for optimizing different parts of the cost function and obtain high-contrast strain images with smooth backgrounds and sharp boundaries. We name our technique ADMM for totaL variaTion RegUlarIzation in ultrasound STrain imaging (ALTRUIST). ALTRUIST's efficacy is quantified using absolute error (AE), Structural SIMilarity (SSIM), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and strain ratio (SR) with respect to GLUE, OVERWIND, and L 1 $L1$ -SOUL, three recently published energy-based techniques, and UMEN-Net, a state-of-the-art deep learning-based algorithm. Analysis of variance (ANOVA)-led multiple comparison tests and paired t $t$ -tests at 5 % $5\%$ overall significance level were conducted to assess the statistical significance of our findings. The Bonferroni correction was taken into account in all statistical tests. Two simulated layer phantoms, three simulated resolution phantoms, one hard-inclusion simulated phantom, one multi-inclusion simulated phantom, one experimental breast phantom, and three in vivo liver cancer datasets have been used for validation experiments. We have published the ALTRUIST code at http://code.sonography.ai. RESULTS: ALTRUIST substantially outperforms the four state-of-the-art benchmarks in all validation experiments, both qualitatively and quantitatively. ALTRUIST yields up to 573 % ∗ ${573\%}^{*}$ , 41 % ∗ ${41\%}^{*}$ , and 51 % ∗ ${51\%}^{*}$ SNR improvements and 443 % ∗ ${443\%}^{*}$ , 53 % ∗ ${53\%}^{*}$ , and 15 % ∗ ${15\%}^{*}$ CNR improvements over L 1 $L1$ -SOUL, its closest competitor, for simulated, phantom, and in vivo liver cancer datasets, respectively, where the asterisk (*) indicates statistical significance. In addition, ANOVA-led multiple comparison tests and paired t $t$ -tests indicate that ALTRUIST generally achieves statistically significant improvements over GLUE, UMEN-Net, OVERWIND, and L 1 $L1$ -SOUL in terms of AE, SSIM map, SNR, and CNR. CONCLUSIONS: A novel ultrasonic displacement tracking algorithm named ALTRUIST has been developed. The principal novelty of ALTRUIST is incorporating ADMM for optimizing an L 1 $L1$ -norm regularization-based cost function. ALTRUIST exhibits promising performance in simulation, phantom, and in vivo experiments.

Displacement Tracking Techniques in Ultrasound Elastography: From Cross Correlation to Deep Learning.

Ultrasound elastography is a noninvasive medical imaging technique that maps viscoelastic properties to characterize tissues and diseases. Elastography can be divided into two classes in a broad sense: strain elastography (SE), which relies on Hooke's law to delineate strain as a surrogate for elasticity, and shear-wave elastography (SWE), which tracks the propagation of shear waves (SWs) in tissues to estimate the elasticity. As tracking the displacement field in the temporal or spatial domain is an inevitable step of both SE and SWE, the success is contingent on the displacement estimation accuracy. Recent reviews mostly focused on clinical applications of elastography, disregarding advances in displacement tracking algorithms. Here, we comprehensively review the recently proposed displacement estimation algorithms applied to both SE and SWE. In addition to cross correlation, block-matching-based (i.e., window-based), model-based, energy-based, and deep learning-based tracking techniques, we review large and lateral displacement tracking, adaptive beamforming, data enhancement, and noise-suppression algorithms facilitating better displacement estimation. We also discuss the simulation models for displacement tracking validation, clinical translation and validation of displacement tracking methods, performance evaluation metrics, and publicly available codes and data for displacement tracking in elastography. Finally, we provide experiential opinions on different tracking algorithms, list the limitations of the current state of elastographic tracking, and comment on possible future research.

Denoising Plane Wave Ultrasound Images Using Diffusion Probabilistic Models.

Ultrasound plane wave imaging is a cutting-edge technique that enables high frame-rate imaging. However, one challenge associated with high frame-rate ultrasound imaging is the high noise associated with them, hindering their wider adoption. Therefore, the development of a denoising method becomes imperative to augment the quality of plane wave images. Drawing inspiration from Denoising Diffusion Probabilistic Models (DDPMs), our proposed solution aims to enhance plane wave image quality. Specifically, the method considers the distinction between low-angle and high-angle compounding plane waves as noise and effectively eliminates it by adapting a DDPM to beamformed radiofrequency (RF) data. The method underwent training using only 400 simulated images. In addition, our approach employs natural image segmentation masks as intensity maps for the generated images, resulting in accurate denoising for various anatomy shapes. The proposed method was assessed across simulation, phantom, and in vivo images. The results of the evaluations indicate that our approach not only enhances image quality on simulated data but also demonstrates effectiveness on phantom and in vivo data in terms of image quality. Comparative analysis with other methods underscores the superiority of our proposed method across various evaluation metrics. The source code and trained model will be released along with the dataset at: http://code.sonography.ai.

Mitigating Aberration-Induced Noise: A Deep Learning-Based Aberration-to-Aberration Approach.

One of the primary sources of suboptimal image quality in ultrasound imaging is phase aberration. It is caused by spatial changes in sound speed over a heterogeneous medium, which disturbs the transmitted waves and prevents coherent summation of echo signals. Obtaining non-aberrated ground truths in real-world scenarios can be extremely challenging, if not impossible. This challenge hinders the performance of deep learning-based techniques due to the domain shift between simulated and experimental data. Here, for the first time, we propose a deep learning-based method that does not require ground truth to correct the phase aberration problem and, as such, can be directly trained on real data. We train a network wherein both the input and target output are randomly aberrated radio frequency (RF) data. Moreover, we demonstrate that a conventional loss function such as mean square error is inadequate for training such a network to achieve optimal performance. Instead, we propose an adaptive mixed loss function that employs both B-mode and RF data, resulting in more efficient convergence and enhanced performance. Finally, we publicly release our dataset, comprising over 180,000 aberrated single plane-wave images (RF data), wherein phase aberrations are modeled as near-field phase screens. Although not utilized in the proposed method, each aberrated image is paired with its corresponding aberration profile and the non-aberrated version, aiming to mitigate the data scarcity problem in developing deep learning-based techniques for phase aberration correction. Source code and trained model are also available along with the dataset at http://code.sonography.ai/main-aaa.

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.

SCANED: Siamese collateral assessment network for evaluation of collaterals from ischemic damage.

This study conducts collateral evaluation from ischemic damage using a deep learning-based Siamese network, addressing the challenges associated with a small and imbalanced dataset. The collateral network provides an alternative oxygen and nutrient supply pathway in ischemic stroke cases, influencing treatment decisions. Research in this area focuses on automated collateral assessment using deep learning (DL) methods to expedite decision-making processes and enhance accuracy. Our study employed a 3D ResNet-based Siamese network, referred to as SCANED, to classify collaterals as good/intermediate or poor. Utilizing non-contrast computed tomography (NCCT) images, the network automates collateral identification and assessment by analyzing tissue degeneration around the ischemic site. Relevant features from the left/right hemispheres were extracted, and Euclidean Distance (ED) was employed for similarity measurement. Finally, dichotomized classification of good/intermediate or poor collateral is performed by SCANED using an optimal threshold derived from ROC analysis. SCANED provides a sensitivity of 0.88, a specificity of 0.63, and a weighted F1 score of 0.86 in the dichotomized classification.

Homodyned K-Distribution Parameter Estimation in Quantitative Ultrasound: Autoencoder and Bayesian Neural Network Approaches.

Quantitative ultrasound (QUS) analyzes the ultrasound (US) backscattered data to find the properties of scatterers that correlate with the tissue microstructure. Statistics of the envelope of the backscattered radio frequency (RF) data can be utilized to estimate several QUS parameters. Different distributions have been proposed to model envelope data. The homodyned K-distribution (HK-distribution) is one of the most comprehensive distributions that can model US backscattered envelope data under diverse scattering conditions (varying scatterer number density and coherent scattering). The scatterer clustering parameter ( α ) and the ratio of the coherent to diffuse scattering power ( k ) are the parameters of this distribution that have been used extensively for tissue characterization in diagnostic US. The estimation of these two parameters (which we refer to as HK parameters) is done using optimization algorithms in which statistical features such as the envelope point-wise signal-to-noise ratio (SNR), skewness, kurtosis, and the log-based moments have been utilized as input to such algorithms. The optimization methods minimize the difference between features and their theoretical value from the HK model. We propose that the true value of these statistical features is a hyperplane that covers a small portion of the feature space. In this article, we follow two approaches to reduce the effect of sample features' error. We propose a model projection neural network based on denoising autoencoders to project the noisy features into this space based on this assumption. We also investigate if the noise distribution can be learned by the deep estimators. We compare the proposed methods with conventional methods using simulations, an experimental phantom, and data from an in vivo animal model of hepatic steatosis. The network weight and a demo code are available online at ht.tp://code.sonography.ai.

MixTURE: L1-Norm-Based Mixed Second-Order Continuity in Strain Tensor Ultrasound Elastography.

Energy-based displacement tracking of ultrasound images can be implemented by optimizing a cost function consisting of a data term, a mechanical congruency term, and first- and second-order continuity terms. This approach recently provided a promising solution to two-dimensional axial and lateral displacement tracking in ultrasound strain elastography. However, the associated second-order regularizer only considers the unmixed second derivatives and disregards the mixed derivatives, thereby providing suboptimal noise suppression and limiting possibilities for total strain tensor imaging. We propose to improve axial, lateral, axial shear, and lateral shear strain estimation quality by formulating and optimizing a novel L1-norm-based second-order regularizer that penalizes both mixed and unmixed displacement derivatives. We name the proposed technique L1-MixTURE, which stands for L1-norm Mixed derivative for Total UltRasound Elastography. When compared to simulated ground truth results, the mean structural similarity (MSSIM) obtained with L1-MixTURE ranged 0.53 to 0.86 and the mean absolute error (MAE) ranged 0.00053 to 0.005. In addition, the mean elastographic signal-to-noise ratio (SNR) achieved with simulated, experimental phantom, and in vivo breast datasets, ranged 1.87 to 52.98, and the mean elastographic contrast-to-noise ratio (CNR) ranged 7.40 to 24.53. When compared to a closely related existing technique that does not consider the mixed derivatives, L1-MixTURE generally outperformed the MSSIM, MAE, SNR, and CNR by up to 37.96%, 67.82%, and 25.53% in the simulated, experimental phantom, and in vivo datasets, respectively. These results collectively highlight the ability of L1-MixTURE to deliver highly accurate axial, lateral, axial shear, and lateral shear strain estimates and advance the state-of-the-art in elastography-guided diagnostic and interventional decisions.

Open access segmentations of intraoperative brain tumor ultrasound images.

PURPOSE: Registration and segmentation of magnetic resonance (MR) and ultrasound (US) images could play an essential role in surgical planning and resectioning brain tumors. However, validating these techniques is challenging due to the scarcity of publicly accessible sources with high-quality ground truth information. To this end, we propose a unique set of segmentations (RESECT-SEG) of cerebral structures from the previously published RESECT dataset to encourage a more rigorous development and assessment of image-processing techniques for neurosurgery. ACQUISITION AND VALIDATION METHODS: The RESECT database consists of MR and intraoperative US (iUS) images of 23 patients who underwent brain tumor resection surgeries. The proposed RESECT-SEG dataset contains segmentations of tumor tissues, sulci, falx cerebri, and resection cavity of the RESECT iUS images. Two highly experienced neurosurgeons validated the quality of the segmentations. DATA FORMAT AND USAGE NOTES: Segmentations are provided in 3D NIFTI format in the OSF open-science platform: https://osf.io/jv8bk. POTENTIAL APPLICATIONS: The proposed RESECT-SEG dataset includes segmentations of real-world clinical US brain images that could be used to develop and evaluate segmentation and registration methods. Eventually, this dataset could further improve the quality of image guidance in neurosurgery.

ULTRA-SR Challenge: Assessment of Ultrasound Localization and TRacking Algorithms for Super-Resolution Imaging.

With the widespread interest and uptake of super-resolution ultrasound (SRUS) through localization and tracking of microbubbles, also known as ultrasound localization microscopy (ULM), many localization and tracking algorithms have been developed. ULM can image many centimeters into tissue in-vivo and track microvascular flow non-invasively with sub-diffraction resolution. In a significant community effort, we organized a challenge, Ultrasound Localization and TRacking Algorithms for Super-Resolution (ULTRA-SR). The aims of this paper are threefold: to describe the challenge organization, data generation, and winning algorithms; to present the metrics and methods for evaluating challenge entrants; and to report results and findings of the evaluation. Realistic ultrasound datasets containing microvascular flow for different clinical ultrasound frequencies were simulated, using vascular flow physics, acoustic field simulation and nonlinear bubble dynamics simulation. Based on these datasets, 38 submissions from 24 research groups were evaluated against ground truth using an evaluation framework with six metrics, three for localization and three for tracking. In-vivo mouse brain and human lymph node data were also provided, and performance assessed by an expert panel. Winning algorithms are described and discussed. The publicly available data with ground truth and the defined metrics for both localization and tracking present a valuable resource for researchers to benchmark algorithms and software, identify optimized methods/software for their data, and provide insight into the current limits of the field. In conclusion, Ultra-SR challenge has provided benchmarking data and tools as well as direct comparison and insights for a number of the state-of-the art localization and tracking algorithms.

DiffeoRaptor: diffeomorphic inter-modal image registration using RaPTOR.

PURPOSE: Diffeomorphic image registration is essential in many medical imaging applications. Several registration algorithms of such type have been proposed, but primarily for intra-contrast alignment. Currently, efficient inter-modal/contrast diffeomorphic registration, which is vital in numerous applications, remains a challenging task. METHODS: We proposed a novel inter-modal/contrast registration algorithm that leverages Robust PaTch-based cOrrelation Ratio metric to allow inter-modal/contrast image alignment and bandlimited geodesic shooting demonstrated in Fourier-Approximated Lie Algebras (FLASH) algorithm for fast diffeomorphic registration. RESULTS: The proposed algorithm, named DiffeoRaptor, was validated with three public databases for the tasks of brain and abdominal image registration while comparing the results against three state-of-the-art techniques, including FLASH, NiftyReg, and Symmetric image Normalization (SyN). CONCLUSIONS: Our results demonstrated that DiffeoRaptor offered comparable or better registration performance in terms of registration accuracy. Moreover, DiffeoRaptor produces smoother deformations than SyN in inter-modal and contrast registration. The code for DiffeoRaptor is publicly available at https://github.com/nimamasoumi/DiffeoRaptor .

Lateral Strain Imaging Using Self-Supervised and Physically Inspired Constraints in Unsupervised Regularized Elastography.

Convolutional Neural Networks (CNN) have shown promising results for displacement estimation in UltraSound Elastography (USE). Many modifications have been proposed to improve the displacement estimation of CNNs for USE in the axial direction. However, the lateral strain, which is essential in several downstream tasks such as the inverse problem of elasticity imaging, remains a challenge. The lateral strain estimation is complicated since the motion and the sampling frequency in this direction are substantially lower than the axial one, and a lack of carrier signal in this direction. In computer vision applications, the axial and the lateral motions are independent. In contrast, the tissue motion pattern in USE is governed by laws of physics which link the axial and lateral displacements. In this paper, inspired by Hooke's law, we, first propose Physically Inspired ConsTraint for Unsupervised Regularized Elastography (PICTURE), where we impose a constraint on the Effective Poisson's ratio (EPR) to improve the lateral strain estimation. In the next step, we propose self-supervised PICTURE (sPICTURE) to further enhance the strain image estimation. Extensive experiments on simulation, experimental phantom and in vivo data demonstrate that the proposed methods estimate accurate axial and lateral strain maps.

Exploiting Mechanics-Based Priors for Lateral Displacement Estimation in Ultrasound Elastography.

Tracking the displacement between the pre- and post-deformed radio-frequency (RF) frames is a pivotal step of ultrasound elastography, which depicts tissue mechanical properties to identify pathologies. Due to ultrasound's poor ability to capture information pertaining to the lateral direction, the existing displacement estimation techniques fail to generate an accurate lateral displacement or strain map. The attempts made in the literature to mitigate this well-known issue suffer from one of the following limitations: 1) Sampling size is substantially increased, rendering the method computationally and memory expensive. 2) The lateral displacement estimation entirely depends on the axial one, ignoring data fidelity and creating large errors. This paper proposes exploiting the effective Poisson's ratio (EPR)-based mechanical correspondence between the axial and lateral strains along with the RF data fidelity and displacement continuity to improve the lateral displacement and strain estimation accuracies. We call our techniques MechSOUL (Mechanically-constrained Second-Order Ultrasound eLastography) and L1 -MechSOUL ( L1 -norm-based MechSOUL), which optimize L2 - and L1 -norm-based penalty functions, respectively. Extensive validation experiments with simulated, phantom, and in vivo datasets demonstrate that MechSOUL and L1 -MechSOUL's lateral strain and EPR estimation abilities are substantially superior to those of the recently-published elastography techniques. We have published the MATLAB codes of MechSOUL and L1 -MechSOUL at https://code.sonography.ai.

Robust landmark-based brain shift correction with a Siamese neural network in ultrasound-guided brain tumor resection.

PURPOSE: In brain tumor surgery, tissue shift (called brain shift) can move the surgical target and invalidate the surgical plan. A cost-effective and flexible tool, intra-operative ultrasound (iUS) with robust image registration algorithms can effectively track brain shift to ensure surgical outcomes and safety. METHODS: We proposed to employ a Siamese neural network, which was first trained using natural images and fine-tuned with domain-specific data to automatically detect matching anatomical landmarks in iUS scans at different surgical stages. An efficient 2.5D approach and an iterative re-weighted least squares algorithm are utilized to perform landmark-based registration for brain shift correction. The proposed method is validated and compared against the state-of-the-art methods using the public BITE and RESECT datasets. RESULTS: Registration of pre-resection iUS scans to during- and post-resection iUS images were executed. The results with the proposed method shows a significant improvement from the initial misalignment ([Formula: see text]) and the method is comparable to the state-of-the-art methods validated on the same datasets. CONCLUSIONS: We have proposed a robust technique to efficiently detect matching landmarks in iUS and perform brain shift correction with excellent performance. It has the potential to improve the accuracy and safety of neurosurgery.