Simulating federated learning for steatosis detection using ultrasound images.

We aimed to implement four data partitioning strategies evaluated with four federated learning (FL) algorithms and investigate the impact of data distribution on FL model performance in detecting steatosis using B-mode US images. A private dataset (153 patients; 1530 images) and a public dataset (55 patient; 550 images) were included in this retrospective study. The datasets contained patients with metabolic dysfunction-associated fatty liver disease (MAFLD) with biopsy-proven steatosis grades and control individuals without steatosis. We employed four data partitioning strategies to simulate FL scenarios and we assessed four FL algorithms. We investigated the impact of class imbalance and the mismatch between the global and local data distributions on the learning outcome. Classification performance was assessed with area under the receiver operating characteristic curve (AUC) on a separate test set. AUCs were 0.93 (95% CI 0.92, 0.94) for source-based partitioning scenario with FedAvg, 0.90 (95% CI 0.89, 0.91) for a centralized model, and 0.83 (95% CI 0.81, 0.85) for a model trained in a single-center scenario. When data was perfectly balanced on the global level and each site had an identical data distribution, the model yielded an AUC of 0.90 (95% CI 0.88, 0.92). When each site contained data exclusively from one single class, irrespective of the global data distribution, the AUC fell in the range of 0.34-0.70. FL applied to B-mode US images provide performance comparable to a centralized model and higher than single-center scenario. Global data imbalance and local data heterogeneity influenced the learning outcome.

Benchmarking the CoW with the TopCoW Challenge: Topology-Aware Anatomical Segmentation of the Circle of Willis for CTA and MRA.

The Circle of Willis (CoW) is an important network of arteries connecting major circulations of the brain. Its vascular architecture is believed to affect the risk, severity, and clinical outcome of serious neuro-vascular diseases. However, characterizing the highly variable CoW anatomy is still a manual and time-consuming expert task. The CoW is usually imaged by two angiographic imaging modalities, magnetic resonance angiography (MRA) and computed tomography angiography (CTA), but there exist limited public datasets with annotations on CoW anatomy, especially for CTA. Therefore we organized the TopCoW Challenge in 2023 with the release of an annotated CoW dataset. The TopCoW dataset was the first public dataset with voxel-level annotations for thirteen possible CoW vessel components, enabled by virtual-reality (VR) technology. It was also the first large dataset with paired MRA and CTA from the same patients. TopCoW challenge formalized the CoW characterization problem as a multiclass anatomical segmentation task with an emphasis on topological metrics. We invited submissions worldwide for the CoW segmentation task, which attracted over 140 registered participants from four continents. The top performing teams managed to segment many CoW components to Dice scores around 90%, but with lower scores for communicating arteries and rare variants. There were also topological mistakes for predictions with high Dice scores. Additional topological analysis revealed further areas for improvement in detecting certain CoW components and matching CoW variant topology accurately. TopCoW represented a first attempt at benchmarking the CoW anatomical segmentation task for MRA and CTA, both morphologically and topologically.

CT analysis of aortic calcifications to predict abdominal aortic aneurysm rupture.

BACKGROUND: Abdominal aortic aneurysm (AAA) rupture prediction based on sex and diameter could be improved. The goal was to assess whether aortic calcification distribution could better predict AAA rupture through machine learning and LASSO regression. METHODOLOGY: In this retrospective study, 80 patients treated for a ruptured AAA between January 2001 and August 2018 were matched with 80 non-ruptured patients based on maximal AAA diameter, age, and sex. Calcification volume and dispersion, morphologic, and clinical variables were compared between both groups using a univariable analysis with p = 0.05 and multivariable analysis through machine learning and LASSO regression. We used AUC for machine learning and odds ratios for regression to measure performance. RESULTS: Mean age of patients was 74.0 ± 8.4 years and 89% were men. AAA diameters were equivalent in both groups (80.9 ± 17.5 vs 79.0 ± 17.3 mm, p = 0.505). Ruptured aneurysms contained a smaller number of calcification aggregates (18.0 ± 17.9 vs 25.6 ± 18.9, p = 0.010) and were less likely to have a proximal neck (45.0% vs 76.3%, p < 0.001). In the machine learning analysis, 5 variables were associated to AAA rupture: proximal neck, antiplatelet use, calcification number, Euclidian distance between calcifications, and standard deviation of the Euclidian distance. A follow-up LASSO regression was concomitant with the findings of the machine learning analysis regarding calcification dispersion but discordant on calcification number. CONCLUSION: There might be more to AAA calcifications that what is known in the present literature. We need larger prospective studies to investigate if indeed, calcification dispersion affects rupture risk. CLINICAL RELEVANCE STATEMENT: Ruptured aneurysms are possibly more likely to have their calcification volume concentrated in a smaller geographical area. KEY POINTS: • Abdominal aortic aneurysm (AAA) rupture prediction based on sex and diameter could be improved. • For a given calcification volume, AAAs with well-distributed calcification clusters could be less likely to rupture. • A machine learning model including AAA calcifications better predicts rupture compared to a model based solely on maximal diameter and sex alone, although it might be prone to overfitting.

Comparison of Radiologists and Deep Learning for US Grading of Hepatic Steatosis.

Background Screening for nonalcoholic fatty liver disease (NAFLD) is suboptimal due to the subjective interpretation of US images. Purpose To evaluate the agreement and diagnostic performance of radiologists and a deep learning model in grading hepatic steatosis in NAFLD at US, with biopsy as the reference standard. Materials and Methods This retrospective study included patients with NAFLD and control patients without hepatic steatosis who underwent abdominal US and contemporaneous liver biopsy from September 2010 to October 2019. Six readers visually graded steatosis on US images twice, 2 weeks apart. Reader agreement was assessed with use of κ statistics. Three deep learning techniques applied to B-mode US images were used to classify dichotomized steatosis grades. Classification performance of human radiologists and the deep learning model for dichotomized steatosis grades (S0, S1, S2, and S3) was assessed with area under the receiver operating characteristic curve (AUC) on a separate test set. Results The study included 199 patients (mean age, 53 years ± 13 [SD]; 101 men). On the test set (n = 52), radiologists had fair interreader agreement (0.34 [95% CI: 0.31, 0.37]) for classifying steatosis grades S0 versus S1 or higher, while AUCs were between 0.49 and 0.84 for radiologists and 0.85 (95% CI: 0.83, 0.87) for the deep learning model. For S0 or S1 versus S2 or S3, radiologists had fair interreader agreement (0.30 [95% CI: 0.27, 0.33]), while AUCs were between 0.57 and 0.76 for radiologists and 0.73 (95% CI: 0.71, 0.75) for the deep learning model. For S2 or lower versus S3, radiologists had fair interreader agreement (0.37 [95% CI: 0.33, 0.40]), while AUCs were between 0.52 and 0.81 for radiologists and 0.67 (95% CI: 0.64, 0.69) for the deep learning model. Conclusion Deep learning approaches applied to B-mode US images provided comparable performance with human readers for detection and grading of hepatic steatosis. Published under a CC BY 4.0 license. Supplemental material is available for this article. See also the editorial by Tuthill in this issue.

Radiomics using computed tomography to predict CD73 expression and prognosis of colorectal cancer liver metastases.

BACKGROUND: Finding a noninvasive radiomic surrogate of tumor immune features could help identify patients more likely to respond to novel immune checkpoint inhibitors. Particularly, CD73 is an ectonucleotidase that catalyzes the breakdown of extracellular AMP into immunosuppressive adenosine, which can be blocked by therapeutic antibodies. High CD73 expression in colorectal cancer liver metastasis (CRLM) resected with curative intent is associated with early recurrence and shorter patient survival. The aim of this study was hence to evaluate whether machine learning analysis of preoperative liver CT-scan could estimate high vs low CD73 expression in CRLM and whether such radiomic score would have a prognostic significance. METHODS: We trained an Attentive Interpretable Tabular Learning (TabNet) model to predict, from preoperative CT images, stratified expression levels of CD73 (CD73High vs. CD73Low) assessed by immunofluorescence (IF) on tissue microarrays. Radiomic features were extracted from 160 segmented CRLM of 122 patients with matched IF data, preprocessed and used to train the predictive model. We applied a five-fold cross-validation and validated the performance on a hold-out test set. RESULTS: TabNet provided areas under the receiver operating characteristic curve of 0.95 (95% CI 0.87 to 1.0) and 0.79 (0.65 to 0.92) on the training and hold-out test sets respectively, and outperformed other machine learning models. The TabNet-derived score, termed rad-CD73, was positively correlated with CD73 histological expression in matched CRLM (Spearman's ρ = 0.6004; P < 0.0001). The median time to recurrence (TTR) and disease-specific survival (DSS) after CRLM resection in rad-CD73High vs rad-CD73Low patients was 13.0 vs 23.6 months (P = 0.0098) and 53.4 vs 126.0 months (P = 0.0222), respectively. The prognostic value of rad-CD73 was independent of the standard clinical risk score, for both TTR (HR = 2.11, 95% CI 1.30 to 3.45, P < 0.005) and DSS (HR = 1.88, 95% CI 1.11 to 3.18, P = 0.020). CONCLUSIONS: Our findings reveal promising results for non-invasive CT-scan-based prediction of CD73 expression in CRLM and warrant further validation as to whether rad-CD73 could assist oncologists as a biomarker of prognosis and response to immunotherapies targeting the adenosine pathway.

Interpretable clinical phenotypes among patients hospitalized with COVID-19 using cluster analysis.

Background: Multiple clinical phenotypes have been proposed for coronavirus disease (COVID-19), but few have used multimodal data. Using clinical and imaging data, we aimed to identify distinct clinical phenotypes in patients admitted with COVID-19 and to assess their clinical outcomes. Our secondary objective was to demonstrate the clinical applicability of this method by developing an interpretable model for phenotype assignment. Methods: We analyzed data from 547 patients hospitalized with COVID-19 at a Canadian academic hospital. We processed the data by applying a factor analysis of mixed data (FAMD) and compared four clustering algorithms: k-means, partitioning around medoids (PAM), and divisive and agglomerative hierarchical clustering. We used imaging data and 34 clinical variables collected within the first 24 h of admission to train our algorithm. We conducted a survival analysis to compare the clinical outcomes across phenotypes. With the data split into training and validation sets (75/25 ratio), we developed a decision-tree-based model to facilitate the interpretation and assignment of the observed phenotypes. Results: Agglomerative hierarchical clustering was the most robust algorithm. We identified three clinical phenotypes: 79 patients (14%) in Cluster 1, 275 patients (50%) in Cluster 2, and 203 (37%) in Cluster 3. Cluster 2 and Cluster 3 were both characterized by a low-risk respiratory and inflammatory profile but differed in terms of demographics. Compared with Cluster 3, Cluster 2 comprised older patients with more comorbidities. Cluster 1 represented the group with the most severe clinical presentation, as inferred by the highest rate of hypoxemia and the highest radiological burden. Intensive care unit (ICU) admission and mechanical ventilation risks were the highest in Cluster 1. Using only two to four decision rules, the classification and regression tree (CART) phenotype assignment model achieved an AUC of 84% (81.5-86.5%, 95 CI) on the validation set. Conclusions: We conducted a multidimensional phenotypic analysis of adult inpatients with COVID-19 and identified three distinct phenotypes associated with different clinical outcomes. We also demonstrated the clinical usability of this approach, as phenotypes can be accurately assigned using a simple decision tree. Further research is still needed to properly incorporate these phenotypes in the management of patients with COVID-19.

Deep Learning: An Update for Radiologists.

Deep learning is a class of machine learning methods that has been successful in computer vision. Unlike traditional machine learning methods that require hand-engineered feature extraction from input images, deep learning methods learn the image features by which to classify data. Convolutional neural networks (CNNs), the core of deep learning methods for imaging, are multilayered artificial neural networks with weighted connections between neurons that are iteratively adjusted through repeated exposure to training data. These networks have numerous applications in radiology, particularly in image classification, object detection, semantic segmentation, and instance segmentation. The authors provide an update on a recent primer on deep learning for radiologists, and they review terminology, data requirements, and recent trends in the design of CNNs; illustrate building blocks and architectures adapted to computer vision tasks, including generative architectures; and discuss training and validation, performance metrics, visualization, and future directions. Familiarity with the key concepts described will help radiologists understand advances of deep learning in medical imaging and facilitate clinical adoption of these techniques. Online supplemental material is available for this article. ©RSNA, 2021.

Predicting the Response to FOLFOX-Based Chemotherapy Regimen from Untreated Liver Metastases on Baseline CT: a Deep Neural Network Approach.

In developed countries, colorectal cancer is the second cause of cancer-related mortality. Chemotherapy is considered a standard treatment for colorectal liver metastases (CLM). Among patients who develop CLM, the assessment of patient response to chemotherapy is often required to determine the need for second-line chemotherapy and eligibility for surgery. However, while FOLFOX-based regimens are typically used for CLM treatment, the identification of responsive patients remains elusive. Computer-aided diagnosis systems may provide insight in the classification of liver metastases identified on diagnostic images. In this paper, we propose a fully automated framework based on deep convolutional neural networks (DCNN) which first differentiates treated and untreated lesions to identify new lesions appearing on CT scans, followed by a fully connected neural networks to predict from untreated lesions in pre-treatment computed tomography (CT) for patients with CLM undergoing chemotherapy, their response to a FOLFOX with Bevacizumab regimen as first-line of treatment. The ground truth for assessment of treatment response was histopathology-determined tumor regression grade. Our DCNN approach trained on 444 lesions from 202 patients achieved accuracies of 91% for differentiating treated and untreated lesions, and 78% for predicting the response to FOLFOX-based chemotherapy regimen. Experimental results showed that our method outperformed traditional machine learning algorithms and may allow for the early detection of non-responsive patients.

Deep learning workflow in radiology: a primer.

Interest for deep learning in radiology has increased tremendously in the past decade due to the high achievable performance for various computer vision tasks such as detection, segmentation, classification, monitoring, and prediction. This article provides step-by-step practical guidance for conducting a project that involves deep learning in radiology, from defining specifications, to deployment and scaling. Specifically, the objectives of this article are to provide an overview of clinical use cases of deep learning, describe the composition of multi-disciplinary team, and summarize current approaches to patient, data, model, and hardware selection. Key ideas will be illustrated by examples from a prototypical project on imaging of colorectal liver metastasis. This article illustrates the workflow for liver lesion detection, segmentation, classification, monitoring, and prediction of tumor recurrence and patient survival. Challenges are discussed, including ethical considerations, cohorting, data collection, anonymization, and availability of expert annotations. The practical guidance may be adapted to any project that requires automated medical image analysis.

Acoustic radiation force induced resonance elastography of coagulating blood: theoretical viscoelasticity modeling and ex-vivo experimentation.

Deep vein thrombosis is a common vascular disease that can lead to pulmonary embolism and death. The early diagnosis and clot age staging are important parameters for reliable therapy planning. This article presents an acoustic radiation force induced resonance elastography method for the viscoelastic characterization of clotting blood. The physical concept of this method relies on the mechanical resonance of the blood clot occurring at specific frequencies. Resonances are induced by focusing ultrasound beams inside the sample under investigation. Coupled to an analytical model of wave scattering, the ability of the proposed method to characterize the viscoelasticity of a mimicked venous thrombosis in the acute phase is demonstrated. Experiments with a gelatin-agar inclusion sample of known viscoelasticity are performed for validation and establishment of the proof of concept. In addition, an inversion method is applied in-vitro for the kinetic monitoring of the blood coagulation process of six human blood samples obtained from two volunteers. The computed elasticity and viscosity values of blood samples at the end of the 90 min kinetics were estimated at 411 ± 71 Pa and 0.25 ± 0.03 Pa.s for volunteer #1, and 387 ± 35 Pa and 0.23 ± 0.02 Pa.s for volunteer #2, respectively. The proposed method allowed reproducible time-varying thrombus viscoelastic measurements from samples having physiological dimensions.

Validation and application of a nondestructive and contactless method for rheological evaluation of biomaterials.

Hydrogels are extensively used for tissue engineering, cell therapy or controlled release of bioactive factors. Nondestructive techniques that can follow their viscoelastic properties during polymerization, remodeling, and degradation are needed, since these properties are determinant for their in vivo efficiency. In this work, we proposed the viscoelastic testing of bilayered materials (VeTBiM) as a new method for nondestructive and contact-less mechanical characterization of soft materials. The VeTBiM method measures the dynamic displacement response of a material, to a low amplitude vibration in order to characterize its viscoelastic properties. We validated VeTBiM by comparing data obtained on various agar and chitosan hydrogels with data from rotational rheometry, and compression tests. We then investigated its potential to follow the mechanical properties of chitosan hydrogels during gelation and in the presence of papain and lysozyme that induce fast or slow enzymatic degradation. Due to this nondestructive and contactless approach, samples can be removed from the instrument and stored in different conditions between measurements. VeTBiM is well adapted to follow biomaterials alone or with cells, over long periods of time. This new method will help in the fine tuning of the mechanical properties of biomaterials used for cell therapy and tissue engineering. © 2016 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 105B: 2565-2573, 2017.

Ultrasound viscoelasticity assessment using an adaptive torsional shear wave propagation method.

PURPOSE: Different approaches have been used in dynamic elastography to assess mechanical properties of biological tissues. Most techniques are based on a simple inversion based on the measurement of the shear wave speed to assess elasticity, whereas some recent strategies use more elaborated analytical or finite element method (FEM) models. In this study, a new method is proposed for the quantification of both shear storage and loss moduli of confined lesions, in the context of breast imaging, using adaptive torsional shear waves (ATSWs) generated remotely with radiation pressure. METHODS: A FEM model was developed to solve the inverse wave propagation problem and obtain viscoelastic properties of interrogated media. The inverse problem was formulated and solved in the frequency domain and its robustness to noise and geometric constraints was evaluated. The proposed model was validated in vitro with two independent rheology methods on several homogeneous and heterogeneous breast tissue-mimicking phantoms over a broad range of frequencies (up to 400 Hz). RESULTS: Viscoelastic properties matched benchmark rheology methods with discrepancies of 8%-38% for the shear modulus G' and 9%-67% for the loss modulus G″. The robustness study indicated good estimations of storage and loss moduli (maximum mean errors of 19% on G' and 32% on G″) for signal-to-noise ratios between 19.5 and 8.5 dB. Larger errors were noticed in the case of biases in lesion dimension and position. CONCLUSIONS: The ATSW method revealed that it is possible to estimate the viscoelasticity of biological tissues with torsional shear waves when small biases in lesion geometry exist.

A Robotic Ultrasound Scanner for Automatic Vessel Tracking and Three-Dimensional Reconstruction of B-Mode Images.

Locating and evaluating the length and severity of a stenosis is very important for planning adequate treatment of peripheral arterial disease (PAD). Conventional ultrasound (US) examination cannot provide maps of entire lower limb arteries in 3-D. We propose a prototype 3D-US robotic system with B-mode images, which is nonionizing, noninvasive, and is able to track and reconstruct a continuous segment of the lower limb arterial tree between the groin and the knee. From an initialized cross-sectional view of the vessel, automatic tracking was conducted followed by 3D-US reconstructions evaluated using Hausdorff distance, cross-sectional area, and stenosis severity in comparison with 3-D reconstructions with computed tomography angiography (CTA). A mean Hausdorff distance of 0.97 ± 0.46 mm was found in vitro for 3D-US compared with 3D-CTA vessel representations. To evaluate the stenosis severity in vitro, 3D-US reconstructions gave errors of 3%-6% when compared with designed dimensions of the phantom, which are comparable to 3D-CTA reconstructions, with 4%-13% errors. The in vivo system's feasibility to reconstruct a normal femoral artery segment of a volunteer was also investigated. These results encourage further ergonomic developments to increase the robot's capacity to represent lower limb vessels in the clinical context.

Generation of remote adaptive torsional shear waves with an octagonal phased array to enhance displacements and reduce variability of shear wave speeds: comparison with quasi-plane shear wavefronts.

A method based on adaptive torsional shear waves (ATSW) is proposed to overcome the strong attenuation of shear waves generated by a radiation force in dynamic elastography. During the inward propagation of ATSW, the magnitude of displacements is enhanced due to the convergence of shear waves and constructive interferences. The proposed method consists in generating ATSW fields from the combination of quasi-plane shear wavefronts by considering a linear superposition of displacement maps. Adaptive torsional shear waves were experimentally generated in homogeneous and heterogeneous tissue mimicking phantoms, and compared to quasi-plane shear wave propagations. Results demonstrated that displacement magnitudes by ATSW could be up to 3 times higher than those obtained with quasi-plane shear waves, that the variability of shear wave speeds was reduced, and that the signal-to-noise ratio of displacements was improved. It was also observed that ATSW could cause mechanical inclusions to resonate in heterogeneous phantoms, which further increased the displacement contrast between the inclusion and the surrounding medium. This method opens a way for the development of new noninvasive tissue characterization strategies based on ATSW in the framework of our previously reported shear wave induced resonance elastography (SWIRE) method proposed for breast cancer diagnosis.

Frequency adaptation for enhanced radiation force amplitude in dynamic elastography.

In remote dynamic elastography, the amplitude of the generated displacement field is directly related to the amplitude of the radiation force. Therefore, displacement improvement for better tissue characterization requires the optimization of the radiation force amplitude by increasing the push duration and/or the excitation amplitude applied on the transducer. The main problem of these approaches is that the Food and Drug Administration (FDA) thresholds for medical applications and transducer limitations may be easily exceeded. In the present study, the effect of the frequency used for the generation of the radiation force on the amplitude of the displacement field was investigated. We found that amplitudes of displacements generated by adapted radiation force sequences were greater than those generated by standard nonadapted ones (i.e., single push acoustic radiation force impulse and supersonic shear imaging). Gains in magnitude were between 20 to 158% for in vitro measurements on agar-gelatin phantoms, and 170 to 336% for ex vivo measurements on a human breast sample, depending on focus depths and attenuations of tested samples. The signal-to-noise ratio was also improved more than 4-fold with adapted sequences. We conclude that frequency adaptation is a complementary technique that is efficient for the optimization of displacement amplitudes. This technique can be used safely to optimize the deposited local acoustic energy without increasing the risk of damaging tissues and transducer elements.

Rheological assessment of a polymeric spherical structure using a three-dimensional shear wave scattering model in dynamic spectroscopy elastography.

With the purpose of assessing localized rheological behavior of pathological tissues using ultrasound dynamic elastography, an analytical shear wave scattering model was used in an inverse problem framework. The proposed method was adopted to estimate the complex shear modulus of viscoelastic spheres from 200 to 450 Hz. The inverse problem was formulated and solved in the frequency domain, allowing assessment of the complex viscoelastic shear modulus at discrete frequencies. A representative rheological model of the spherical obstacle was determined by comparing storage and loss modulus behaviors with Kelvin-Voigt, Maxwell, Zener, and Jeffrey models. The proposed inversion method was validated by using an external vibrating source and acoustic radiation force. The estimation of viscoelastic properties of three-dimensional spheres made softer or harder than surrounding tissues did not require a priori rheological assumptions. The proposed method is intended to be applied in the context of breast cancer imaging.

Viscoelastic characterization of elliptical mechanical heterogeneities using a semi-analytical shear-wave scattering model for elastometry measures.

This paper presents a semi-analytical model of shear wave scattering by a viscoelastic elliptical structure embedded in a viscoelastic medium, and its application in the context of dynamic elastography imaging. The commonly used assumption of mechanical homogeneity in the inversion process is removed introducing a priori geometrical information to model physical interactions of plane shear waves with the confined mechanical heterogeneity. Theoretical results are first validated using the finite element method for various mechanical configurations and incidence angles. Secondly, an inverse problem is formulated to assess viscoelastic parameters of both the elliptic inclusion and its surrounding medium, and applied in vitro to characterize mechanical properties of agar-gelatin phantoms. The robustness of the proposed inversion method is then assessed under various noise conditions, biased geometrical parameters and compared to direct inversion, phase gradient and time-of-flight methods. The proposed elastometry method appears reliable in the context of estimating confined lesion viscoelastic parameters.

Shear wave induced resonance elastography of venous thrombi: a proof-of-concept.

Shear wave induced resonance elastography (SWIRE) is proposed for deep venous thrombosis (DVT) elasticity assessment. This new imaging technique takes advantage of properly polarized shear waves to induce resonance of a confined mechanical heterogeneity. Realistic phantoms (n = 9) of DVT total and partial clot occlusions with elasticities from 406 to 3561 Pa were built for in vitro experiments. An ex vivo study was also performed to evaluate the elasticity of two fresh porcine venous thrombi in a pig model. Transient shear waves at 45-205 Hz were generated by the vibration of a rigid plate (plane wavefront) or by a needle to simulate a radiation pressure on a line segment (cylindrical wavefront). Induced propagation of shear waves was imaged with an ultrafast ultrasound scanner and a finite element method was developed to simulate tested experimental conditions. An inverse problem was then formulated considering the first resonance frequency of the DVT inclusion. Elasticity agreements between SWIRE and a reference spectroscopy instrument (RheoSpectris) were found in vitro for total clots either in plane (r(2) = 0.989) or cylindrical (r(2) = 0.986) wavefront configurations. For total and partial clots, elasticity estimation errors were 9.0 ±4.6% and 9.3 ±11.3%, respectively. Ex vivo, the blood clot elasticity was 498 ±58 Pa within the inferior vena cava and 436 ±45 Pa in the right common iliac vein (p = 0.22). To conclude, the SWIRE technique seems feasible to quantitatively assess blood clot elasticity in the context of DVT ultrasound imaging.

Stolt's f-k migration for plane wave ultrasound imaging.

Ultrafast ultrasound is an emerging modality that offers new perspectives and opportunities in medical imaging. Plane wave imaging (PWI) allows one to attain very high frame rates by transmission of planar ultrasound wave-fronts. As a plane wave reaches a given scatterer, the latter becomes a secondary source emitting upward spherical waves and creating a diffraction hyperbola in the received RF signals. To produce an image of the scatterers, all the hyperbolas must be migrated back to their apexes. To perform beamforming of plane wave echo RFs and return high-quality images at high frame rates, we propose a new migration method carried out in the frequency-wavenumber (f-k) domain. The f-k migration for PWI has been adapted from the Stolt migration for seismic imaging. This migration technique is based on the exploding reflector model (ERM), which consists in assuming that all the scatterers explode in concert and become acoustic sources. The classical ERM model, however, is not appropriate for PWI. We showed that the ERM can be made suitable for PWI by a spatial transformation of the hyperbolic traces present in the RF data. In vitro experiments were performed to outline the advantages of PWI with Stolt's f-k migration over the conventional delay-and-sum (DAS) approach. The Stolt's f-k migration was also compared with the Fourier-based method developed by J.-Y. Lu. Our findings show that multi-angle compounded f-k migrated images are of quality similar to those obtained with a stateof- the-art dynamic focusing mode. This remained true even with a very small number of steering angles, thus ensuring a highly competitive frame rate. In addition, the new FFT-based f-k migration provides comparable or better contrast-to-noise ratio and lateral resolution than the Lu's and DAS migration schemes. Matlab codes for the Stolt's f-k migration for PWI are provided.