Cancer characterization using light backscattering spectroscopy and quantitative ultrasound: an ex vivo study on sarcoma subtypes

Histological analysis is the gold standard method for cancer diagnosis. However, it is prone to subjectivity and sampling bias. In response to these limitations, we introduce a quantitative bimodal approach that aims to provide non-invasive guidance towards suspicious regions. Light backscattering spectroscopy and quantitative ultrasound techniques were combined to characterize two different bone tumor types from animal models: chondrosarcomas and osteosarcomas. Two different cell lines were used to induce osteosarcoma growth. Histological analyses were conducted to serve as references. Three ultrasound parameters and intensities of the light reflectance profiles showed significant differences between chondrosarcomas and osteosarcomas at the 5% level. Likewise, variations in the same biomarkers were reported for the two types of osteosarcoma, despite their similar morphology observed in the histological examinations. These observations show the sensitivity of our techniques in probing fine tissue properties. Secondly, the ultrasound spectral-based technique identified the mean size of chondrosarcoma cells and nuclei with relative errors of about 22% and 9% respectively. The optical equivalent technique correctly extracted scatterer size distributions that encompass nuclei and cells for chondrosarcomas and osteosarcomas ([Formula: see text] and [Formula: see text] respectively). The optical scattering contributions of nuclei were estimated at 52% for the chondrosarcomas and 69% for the osteosarcomas, probably indicating the abundant and the absent extracellular matrix respectively. Thus, the ultrasound and the optical methods brought complementary parameters. They successfully estimated morphological parameters at the cellular and the nuclear scales, making our bimodal technique promising for tumor characterization.

Sensitivity Enhancement Using Chirp Transmission for an Ultrasound Arthroscopic Probe

Meniscal tear in the knee joint is a highly common injury that can require an ablation. However, the success rate of meniscectomy is highly impacted by difficulties in estimating the thin vascularization of the meniscus, which determines the healing capacities of the patient. Indeed, vascularization is estimated using arthroscopic cameras that lack of high sensitivity to blood flow. Here, we propose an ultrasound method for estimating the density of vascularization in the meniscus during surgery. This approach uses an arthroscopic probe driven by ultrafast sequences. To enhance the sensitivity of the method, we propose to use a chirp-coded excitation combined with a mismatched compression filter robust to the attenuation. This chirp approach was compared to a standard ultrafast emission and a Hadamard-coded emission using a flow phantom. The mismatched filter was also compared to a matched filter. Results show that, for a velocity of a few millimeters per second, the mismatched filter gives a 4.4-10.4-dB increase of the signal-to-noise ratio (SNR) compared to the Hadamard emission and a 3.1-6.6-dB increase compared to the matched filter. Such increases are obtained for a loss of axial resolution of 13% when comparing the point spread functions (PSFs) of the mismatched and matched filters. Hence, the mismatched filter allows increasing significantly the probe capacity to detect slow flows at the cost of a small loss in axial resolution. This preliminary study is the first step toward an ultrasensitive ultrasound arthroscopic probe able to assist the surgeon during meniscectomy.

Prognostic significance of vascular and valvular calcifications in low- and high-gradient aortic stenosis

AIMS

In low-gradient aortic stenosis (LGAS), the high valvulo-arterial impedance observed despite low valvular gradient suggests a high vascular load. Thoracic aortic calcifications (TACs) and valvular aortic calcifications (VACs) are, respectively, surrogates of aortic load and aortic valvular gradient. The aim of this study was to compare the respective contributions of TAC and VAC on 3-year cardiovascular (CV) mortality following TAVI in LGAS vs. high-gradient aortic stenosis (HGAS) patients.

METHODS AND RESULTS

A total of 1396 consecutive patients were included. TAC and VAC were measured on the pre-TAVI CT-scan. About 435 (31.2%) patients had LGAS and 961 (68.8%) HGAS. LGAS patients were more prone to have diabetes, coronary artery disease (CAD), atrial fibrillation (AF), and lower left ventricular ejection fraction (LVEF), P<0.05 for all. During the 3 years after TAVI, 245(17.8%) patients experienced CV mortality, 92(21.6%) in LGAS and 153(16.2%) in HGAS patients, P=0.018. Multivariate analysis adjusted for age, gender, diabetes, AF, CAD, LVEF, renal function, vascular access, and aortic regurgitation showed that TAC but not VAC was associated with CV mortality in LGAS, hazard ratio (HR) 1.085 confidence interval (CI) (1.019-1.156), P=0.011, and HR 0.713 CI (0.439-1.8), P=0.235; the opposite was observed in HGAS patients with VAC but not TAC being associated with CV mortality, HR 1.342 CI (1.034-1.742), P=0.027, and HR 1.015 CI (0.955-1.079), P=0.626.

CONCLUSION 

TAC plays a major prognostic role in LGAS while VAC remains the key in HGAS patients. This confirms that LGAS is a complex vascular and valvular disease.

P907Respective pronostic value of the valvular aortic calcifications and the thoracic aorta calcifications in patients with and without low gradient aortic stenosis after TAVI

Background

Both the valvular aortic calcifications (VAC) and the thoracic aorta calcifications (TAC) have a prognostic impact in patients with aortic stenosis. Their respective prognostic values in patients with and without low gradient aortic stenosis (LGAS) remain unknown after TAVI.

Objectives

To assess the prognostic significance of VAC and TAC in patients with and without LGAS regarding cardiovascular mortality after 3 years follow-up.

Methods

The CAPRI-LGAS is an ancillary study of the C4CAPRI trial (NCT02935491) including 1282 consecutive TAVI patients. Calcifications were measured on pre-TAVI CT. The primary outcome was defined as cardiovascular mortality 3 years after TAVI.

Results

Among the 1282 patients, 397 (31%) had a LGAS. Compared to the other patients, LGAS patients were more prone to be men, younger, with atrial fibrillation, and lower left ventricular ejection fraction (LVEF), p<0.05 for all. No statistically significant difference was noticed for pulmonary systolic pressures, history of diabetes, chronic respiratory disease, renal insufficiency or peripheral vascular disease. VAC was lower in LGAS compared to non-LGAS patients (1.05 cm3±0.7 vs 0.75 cm3±0.5), p<0.001, the contrary was noticed for TAC, (3.1 cm3±3 vs 3.7 cm3±3.7), p=0.011. After 3 years follow-up, 227 (17.7%) patients died from cardiovascular causes; respectively 85 (21.4%) and 142 (16.1%) patients with and without LGAS, p=0.02. In univariate analysis, in LGAS patients each increase of 1cm3 TAC was associated with cardiovascular mortality while VAC was not, respectively Hazard Ratio (HR) 1.07 and confidence interval (CI) (1.023–1.119) p=0.003, and HR 0.822 CI (0.523–1.292), p=0.39. In patients without LGAS both TAC and VAC were associated with mortality, respectively HR 1.054 CI (1.006–1.104), p=0.028 and HR 1.363 CI (1.092–1.701), p=0.006. Multivariate analysis was adjusted for TAC, VAC, age, gender, atrial fibrillation, and LVEF. In LGAS patients TAC but not VAC was still a predictor of cardiovascular mortality, respectively HR 1.092 CI (1.031–1.158), p=0.003, and HR 0.743 CI (0.464–1.191), p=0.21. In patients without LGAS TAC was no more associated with cardiovascular mortality while VAC was, respectively HR 1.306 CI (1.024–1.666), p=0.031, and HR 1.038 CI (0.985–1.094), p=0.161. When further adjusting on pulmonary systolic pressures, history of diabetes, chronic respiratory disease, renal insufficiency and peripheral vascular disease, the results remained similar ie in LGAS patients, TAC HR 1.090 CI (1.022–1.162), p=0.009 while in patients without LGAS VAC HR 1.377 CI (1.049–1.809), p=0.021.

Conclusions

The present study shows that VAC and TAC involve different prognostic information in patients with and without LGAS after TAVI. While VAC may be a marker of early and periprocedural mortality and aortic regurgitation in non-LGAS patients, TAC may continue to be harmful and increase afterload in patients with LGAS whom LVEF is often impaired.

Evaluation and comparison of current biopsy needle localization and tracking methods using 3D ultrasound

This article compares four different biopsy needle localization algorithms in both 3D and 4D situations to evaluate their accuracy and execution time. The localization algorithms were: Principle component analysis (PCA), random Hough transform (RHT), parallel integral projection (PIP) and ROI-RK (ROI based RANSAC and Kalman filter). To enhance the contrast of the biopsy needle and background tissue, a line filtering pre-processing step was implemented. To make the PCA, RHT and PIP algorithms comparable with the ROI-RK method, a region of interest (ROI) strategy was added. Simulated and ex-vivo data were used to evaluate the performance of the different biopsy needle localization algorithms. The resolutions of the sectorial and cylindrical volumes were 0.3mm×0.4mm×0.6mmand0.1mm×0.1mm×0.2mm (axial×lateral×azimuthal) respectively. In so far as the simulation and experimental results show, the ROI-RK method successfully located and tracked the biopsy needle in both 3D and 4D situations. The tip localization error was within 1.5mm and the axis accuracy was within 1.6mm. To the best of our knowledge, considering both localization accuracy and execution time, the ROI-RK was the most stable and time-saving method. Normally, accuracy comes at the expense of time. However, the ROI-RK method was able to locate the biopsy needle with high accuracy in real time, which makes it a promising method for clinical applications.

Plane-wave transverse oscillation for high-frame-rate 2-D vector flow imaging

Transverse oscillation (TO) methods introduce oscillations in the pulse-echo field (PEF) along the direction transverse to the ultrasound propagation direction. This may be exploited to extend flow investigations toward multidimensional estimates. In this paper, the TOs are coupled with the transmission of plane waves (PWs) to reconstruct high-framerate RF images with bidirectional oscillations in the pulse-echo field. Such RF images are then processed by a 2-D phase-based displacement estimator to produce 2-D vector flow maps at thousands of frames per second. First, the capability of generating TOs after PW transmissions was thoroughly investigated by varying the lateral wavelength, the burst length, and the transmission frequency. Over the entire region of interest, the generated lateral wavelengths, compared with the designed ones, presented bias and standard deviation of -3.3 ± 5.7% and 10.6 ± 7.4% in simulations and experiments, respectively. The performance of the ultrafast vector flow mapping method was also assessed by evaluating the differences between the estimated velocities and the expected ones. Both simulations and experiments show overall biases lower than 20% when varying the beam-to-flow angle, the peak velocity, and the depth of interest. In vivo applications of the method on the common carotid and the brachial arteries are also presented.

Compressed Sensing Reconstruction of 3D Ultrasound Data Using Dictionary Learning and Line-Wise Subsampling

In this paper we present a compressed sensing (CS) method adapted to 3D ultrasound imaging (US). In contrast to previous work, we propose a new approach based on the use of learned overcomplete dictionaries that allow for much sparser representations of the signals since they are optimized for a particular class of images such as US images. In this study, the dictionary was learned using the K-SVD algorithm and CS reconstruction was performed on the non-log envelope data by removing 20% to 80% of the original data. Using numerically simulated images, we evaluate the influence of the training parameters and of the sampling strategy. The latter is done by comparing the two most common sampling patterns, i.e., point-wise and line-wise random patterns. The results show in particular that line-wise sampling yields an accuracy comparable to the conventional point-wise sampling. This indicates that CS acquisition of 3D data is feasible in a relatively simple setting, and thus offers the perspective of increasing the frame rate by skipping the acquisition of RF lines. Next, we evaluated this approach on US volumes of several ex vivo and in vivo organs. We first show that the learned dictionary approach yields better performances than conventional fixed transforms such as Fourier or discrete cosine. Finally, we investigate the generality of the learned dictionary approach and show that it is possible to build a general dictionary allowing to reliably reconstruct different volumes of different ex vivo or in vivo organs.

Real-time ultrasound-tagging to track the 2D motion of the common carotid artery wall in vivo

PURPOSE

Tracking the motion of biological tissues represents an important issue in the field of medical ultrasound imaging. However, the longitudinal component of the motion (i.e., perpendicular to the beam axis) remains more challenging to extract due to the rather coarse resolution cell of ultrasound scanners along this direction. The aim of this study is to introduce a real-time beamforming strategy dedicated to acquire tagged images featuring a distinct pattern in the objective to ease the tracking.

METHODS

Under the conditions of the Fraunhofer approximation, a specific apodization function was applied to the received raw channel data, in real-time during image acquisition, in order to introduce a periodic oscillations pattern along the longitudinal direction of the radio frequency signal. Analytic signals were then extracted from the tagged images, and subpixel motion tracking of the intima-media complex was subsequently performed offline, by means of a previously introduced bidimensional analytic phase-based estimator.

RESULTS

The authors' framework was applied in vivo on the common carotid artery from 20 young healthy volunteers and 6 elderly patients with high atherosclerosis risk. Cine-loops of tagged images were acquired during three cardiac cycles. Evaluated against reference trajectories manually generated by three experienced analysts, the mean absolute tracking error was 98 ± 84 μm and 55 ± 44 μm in the longitudinal and axial directions, respectively. These errors corresponded to 28% ± 23% and 13% ± 9% of the longitudinal and axial amplitude of the assessed motion, respectively.

CONCLUSIONS

The proposed framework enables tagged ultrasound images of in vivo tissues to be acquired in real-time. Such unconventional beamforming strategy contributes to improve tracking accuracy and could potentially benefit to the interpretation and diagnosis of biomedical images.

Full 3-D transverse oscillations: a method for tissue motion estimation

We present a new method to estimate 4-D (3-D + time) tissue motion. The method used combines 3-D phase based motion estimation with an unconventional beamforming strategy. The beamforming technique allows us to obtain full 3-D RF volumes with axial, lateral, and elevation modulations. Based on these images, we propose a method to estimate 3-D motion that uses phase images instead of amplitude images. First, volumes featuring 3-D oscillations are created using only a single apodization function, and the 3-D displacement between two consecutive volumes is estimated simultaneously by applying this 3-D estimation. The validity of the method is investigated by conducting simulations and phantom experiments. The results are compared with those obtained with two other conventional estimation methods: block matching and optical flow. The results show that the proposed method outperforms the conventional methods, especially in the transverse directions.

Back-propagation beamformer design for motion estimation in echocardiography

Transverse oscillation (TO) techniques have shown their potential for improving the accuracy of local motion estimation in the transverse direction (i.e., the direction perpendicular to the beam axis). The conventional design of TOs in linear geometry, which is based on the Fraunhofer approximation, relates point spread function (PSF) and apodization function through a Fourier transform. Motivated by the adaptation of TOs in echocardiography, we propose a specific beamforming approach based on back-propagation (BP) to build TOs in sector-shaped geometry. Numerical simulations and experimental data give a comparison between proposed and conventional beamforming for TOs. The accuracy is first quantified by comparing the generated and theoretical PSF using the root mean square error (RMSE) and shows that BP-based beamforming approximates the desired TOs more closely than the conventional approach. Motion estimation is then evaluated. The axial and lateral displacements are within the range [0-0.6] mm and [0°-6.4°], respectively, which correspond to 0.8 times the axial (0.73 mm) and lateral (8°) wavelengths. The result shows that the proposed method yields a clear improvement for lateral displacements, by reducing the error by 28.6% compared with Fourier transform-based beamforming, while maintaining the same error for axial motion estimation. Experimental measurements are discussed to complete this study and confirm that BP-based beamforming leads to better controlled TO images than conventional Fourier-based beamforming.

Comparison of carotid artery blood velocity measurements by vector and standard Doppler approaches

Although severely affected by the angle dependency, carotid artery peak systolic velocity measurements are widely used for assessment of stenosis. In this study, blood peak systolic velocities in the common and internal carotid arteries of both healthy volunteers and patients with internal carotid artery stenosis were measured by two vector Doppler (VD) methods and compared with measurements obtained with the conventional spectral Doppler approach. Although the two VD techniques were completely different (using the transmission of focused beams and plane waves, respectively), the measurement results indicate that these techniques are nearly equivalent. The peak systolic velocities measured in 22 healthy common carotid arteries by the two VD techniques were very close (according to Bland-Altman analysis, the average difference was 3.2%, with limits of agreement of ± 8.6%). Application of Bland-Altman analysis to comparison of either VD technique with the spectral Doppler method provided a 21%-25% average difference with ± 13%-15% limits of agreement. Analysis of the results obtained from 15 internal carotid arteries led to similar conclusions, indicating significant overestimation of peak systolic velocity with the spectral Doppler method. Inter- and intra-operator repeatability measurements performed in a group of 8 healthy volunteers provided equivalent results for all of the methods (coefficients of variability in the range 2.7%-6.9%), even though the sonographers were not familiar with the VD methods. The results of this study suggest that the introduction of vector Doppler methods in commercial machines may finally be considered mature and capable of overcoming the angle-dependent overestimation typical of the standard spectral Doppler approach.

Experimental evaluation of spectral-based quantitative ultrasound imaging using plane wave compounding

Quantitative ultrasound (QUS) based on backscatter coefficient (BSC) estimation has shown potential for tissue characterization. Beamforming using plane wave compounding has advantages for echographic, Doppler, and elastographic imaging; however, to date, plane wave compounding has not been experimentally evaluated for the purpose of BSC estimation. In this study, two BSC-derived parameters (i.e., the BSC midband fit and intercept) were estimated from experimental data obtained using compound plane wave beamforming. For comparison, QUS parameters were also estimated from data obtained using both fixed focus and dynamic receive beamforming. An ultrasound imaging system equipped with a 9-MHz center frequency, 64-element array was used to collect data up to a depth of 45 mm. Two gelatin phantoms with randomly distributed 20-μm inclusions with a homogeneous scatterer concentration and a two-region scatterer concentration were used for assessing the precision and lateral resolution of QUS imaging, respectively. The use of plane wave compounding resulted in accurate QUS estimation (i.e., bias in the BSC parameters of less than 2 dB) and relatively constant lateral resolution (i.e., BSC midband fit 10% to 90% rise distance ranging between 1.0 and 1.5 mm) throughout a 45 mm field of view. Although both fixed focus and dynamic receive beamforming provided the same performance around the focal depth, the reduction in SNR away from the focus resulted in a reduced field of view in the homogeneous phantom (i.e., only 28 mm). The lateral resolution also degraded away from the focus, with up to a 2-fold and 10-fold increase in the rise distance at 20 mm beyond the focal depth for dynamic receive and fixed focus beamforming, respectively. These results suggest that plane wave compounding has the potential to improve the performance of spectral-based quantitative ultrasound over other conventional beamforming strategies.

High-frame-rate 2-D vector blood flow imaging in the frequency domain

Conventional ultrasound Doppler techniques estimate the blood velocity exclusively in the axial direction to produce the sonograms and color flow maps needed for diagnosis of cardiovascular diseases. In this paper, a novel method to produce bi-dimensional maps of 2-D velocity vectors is proposed. The region of interest (ROI) is illuminated by plane waves transmitted at the pulse repetition frequency (PRF) in a fixed direction. For each transmitted plane wave, the backscattered echoes are recombined offline to produce the radio-frequency image of the ROI. The local 2-D phase shifts between consecutive speckle images are efficiently estimated in the frequency domain, to produce vector maps up to 15 kHz PRF. Simulations and in vitro steady-flow experiments with different setup conditions have been conducted to thoroughly evaluate the method's performance. Bias is proved to be lower than 10% in most simulations and lower than 20% in experiments. Further simulations and in vivo experiments have been made to test the approach's feasibility in pulsatile flow conditions. It has been estimated that the computation of the frequency domain algorithm is more than 50 times faster than the computation of the reference 2-D cross-correlation algorithm.

OntoVIP: an ontology for the annotation of object models used for medical image simulation

This paper describes the creation of a comprehensive conceptualization of object models used in medical image simulation, suitable for major imaging modalities and simulators. The goal is to create an application ontology that can be used to annotate the models in a repository integrated in the Virtual Imaging Platform (VIP), to facilitate their sharing and reuse. Annotations make the anatomical, physiological and pathophysiological content of the object models explicit. In such an interdisciplinary context we chose to rely on a common integration framework provided by a foundational ontology, that facilitates the consistent integration of the various modules extracted from several existing ontologies, i.e. FMA, PATO, MPATH, RadLex and ChEBI. Emphasis is put on methodology for achieving this extraction and integration. The most salient aspects of the ontology are presented, especially the organization in model layers, as well as its use to browse and query the model repository.

A new technique for the estimation of cardiac motion in echocardiography based on transverse oscillations: a preliminary evaluation in silico and a feasibility demonstration in vivo

Quantification of regional myocardial motion and deformation from cardiac ultrasound is fostering considerable research efforts. Despite the tremendous improvements done in the field, all existing approaches still face a common limitation which is intrinsically connected with the formation of the ultrasound images. Specifically, the reduced lateral resolution and the absence of phase information in the lateral direction highly limit the accuracy in the computation of lateral displacements. In this context, this paper introduces a novel setup for the estimation of cardiac motion with ultrasound. The framework includes an unconventional beamforming technique and a dedicated motion estimation algorithm. The beamformer aims at introducing phase information in the lateral direction by producing transverse oscillations. The estimator directly exploits the phase information in the two directions by decomposing the image into two 2-D single-orthant analytic signals. An in silico evaluation of the proposed framework is presented on five ultra-realistic simulated echocardiographic sequences, where the proposed motion estimator is contrasted against other two phase-based solutions exploiting the presence of transverse oscillations and against block-matching on standard images. An implementation of the new beamforming strategy on a research ultrasound platform is also shown along with a preliminary in vivo evaluation on one healthy subject.

Automatic needle detection and tracking in 3D ultrasound using an ROI-based RANSAC and Kalman method

This article proposes a robust technique for needle detection and tracking using three-dimensional ultrasound (3D US). It is difficult for radiologists to detect and follow the position of micro tools, such as biopsy needles, that are inserted in human tissues under 3D US guidance. To overcome this difficulty, we propose a method that automatically reduces the processed volume to a limited region of interest (ROI), increasing at the same time the calculation speed and the robustness of the proposed technique. First, a line filter method that enhances the contrast of the needle against the background is used to facilitate the initialization of ROI using the tubularness information of the complete US volume. Then, the random sample consensus (RANSAC) and Kalman filter (RK) algorithm is used in the ROI to detect and track the precise position of the needle. A series of numerical inhomogeneous phantoms with a needle simulated from real 3D US volumes are used to evaluate our method. The results show that the proposed method is much more robust than the RANSAC algorithm when detecting the needle, regardless of whether or not the insertion axis corresponds to an acquisition plane in the 3D US volume. The possibility of failure is also discussed in this article.

3D strain assessment in ultrasound (Straus): a synthetic comparison of five tracking methodologies

This paper evaluates five 3D ultrasound tracking algorithms regarding their ability to quantify abnormal deformation in timing or amplitude. A synthetic database of B-mode image sequences modeling healthy, ischemic and dyssynchrony cases was generated for that purpose. This database is made publicly available to the community. It combines recent advances in electromechanical and ultrasound modeling. For modeling heart mechanics, the Bestel-Clement-Sorine electromechanical model was applied to a realistic geometry. For ultrasound modeling, we applied a fast simulation technique to produce realistic images on a set of scatterers moving according to the electromechanical simulation result. Tracking and strain accuracies were computed and compared for all evaluated algorithms. For tracking, all methods were estimating myocardial displacements with an error below 1 mm on the ischemic sequences. The introduction of a dilated geometry was found to have a significant impact on accuracy. Regarding strain, all methods were able to recover timing differences between segments, as well as low strain values. On all cases, radial strain was found to have a low accuracy in comparison to longitudinal and circumferential components.

Multi-resolution transverse oscillation in ultrasound imaging for motion estimation

Beyond all recent methods for estimating multidimensional motion vectors in ultrasound (US) image sequences, the transverse oscillation (TO) images coupled with a phase-based motion estimator has shown to be an extremely promising approach. To improve existing TO techniques, this paper proposes a multi-resolution transverse oscillation (MRTO) that measures the motion in US image sequences at different resolutions. Several TO images are created using harmonic images obtained from nonlinear US propagation or from composite emissions. Motion estimations are made with increasing resolution levels, i.e., decreasing transverse wavelength, to improve the previous measurements. When only lateral motion is applied to the tissue, the resulting estimations, quantitatively assessed in both simulations and in vitro experiments, are significantly improved in terms of mean motion error and standard deviation. With a 1-mm displacement in the chirp transmission, the error is reduced from 40.9% to 0.1% in the simulation and from 41.6% to 1.5% in the experiment. The first quasi-static elastography results, in simulation and experimentation, also confirm this improvement with an increased contrast-to-noise ratio and signal-to-noise ratio.

Pre-beamformed RF signal reconstruction in medical ultrasound using compressive sensing

Compressive sensing (CS) theory makes it possible - under certain assumptions - to recover a signal or an image sampled below the Nyquist sampling limit. In medical ultrasound imaging, CS could allow lowering the amount of acquired data needed to reconstruct the echographic image. CS thus offers the perspective of speeding up echographic acquisitions and could have many applications, e.g. triplex acquisitions for CFM/B-mode/Doppler imaging, high-frame-rate echocardiography, 3D imaging using matrix probes, etc. The objective of this paper is to study the feasibility of CS for the reconstruction of channel RF data, i.e. the 2D set of raw RF lines gathered at the receive elements. Successful application of CS implies selecting a representation basis where the data to be reconstructed have a sparse expansion. Because they consist mainly in warped oscillatory patterns, channel RF data do not easily lend themselves to a sparse representation and thus represent a specific challenge. Within this perspective, we propose to perform and assess CS reconstruction of channel RF data using the recently introduced wave atoms [1] representation, which exhibit advantageous properties for sparsely representing such oscillatory patterns. Reconstructions obtained using wave atoms are compared with the reconstruction performed with two conventional representation bases, namely Fourier and Daubechies wavelets. The first experiment was conducted on simulated channel RF data acquired from a numerical cyst phantom. The quality of the reconstructions was quantified through the mean absolute error at varying subsampling rates by removing 50-90% of the original samples. The results obtained for channel RF data reconstruction yield error ranges of [0.6-3.0]×10(-2), [0.2-2.6]×10(-2), [0.1-1.5]×10(-2), for wavelets, Fourier and wave atoms respectively. The error ranges observed for the associated beamformed log-envelope images are [2.4-20.6]dB, [1.1-12.2]dB, and [0.5-8.8dB] using wavelets, Fourier, and wave atoms, respectively. These results thus show the superiority of the wave atom representation and the feasibility of CS for the reconstruction of US RF data. The second experiment aimed at showing the experimental feasibility of the method proposed using a data set acquired by imaging a general-purpose phantom (CIRS Model 054GS) using an Ultrasonix MDP scanner. The reconstruction was performed by removing 80% of the initial samples and using wave atoms. The reconstructed image was found to reliably preserve the speckle structure and was associated with an error of 5.5dB.

Internal strain estimation for quantification of human heel pad elastic modulus: A phantom study

Shock absorption is the most important function of the human heel pad. However, changes in heel pad elasticity, as seen in e.g. long-distance runners, diabetes patients, and victims of Falanga torture are affecting this function, often in a painful manner. Assessment of heel pad elasticity is usually based on one or a few strain measurements obtained by an external load-deformation system. The aim of this study was to develop a technique for quantitative measurements of heel pad elastic modulus based on several internal strain measures from within the heel pad by use of ultrasound images. Nine heel phantoms were manufactured featuring a combination of three heel pad stiffnesses and three heel pad thicknesses to model the normal human variation. Each phantom was tested in an indentation system comprising a 7MHz linear array ultrasound transducer, working as the indentor, and a connected load cell. Load-compression data and ultrasound B-mode images were simultaneously acquired in 19 compression steps of 0.1mm each. The internal tissue displacement was for each step calculated by a phase-based cross-correlation technique and internal strain maps were derived from these displacement maps. Elastic moduli were found from the resulting stress-strain curves. The elastic moduli made it possible to distinguish eight of nine phantoms from each other according to the manufactured stiffness and showed very little dependence of the thickness. Mean elastic moduli for the three soft, the three medium, and the three hard phantoms were 89kPa, 153kPa, and 168kPa, respectively. The combination of ultrasound images and force measurements provided an effective way of assessing the elastic properties of the heel pad due to the internal strain estimation.

A virtual imaging platform for multi-modality medical image simulation

This paper presents the Virtual Imaging Platform (VIP), a platform accessible at http://vip.creatis.insa-lyon.fr to facilitate the sharing of object models and medical image simulators, and to provide access to distributed computing and storage resources. A complete overview is presented, describing the ontologies designed to share models in a common repository, the workflow template used to integrate simulators, and the tools and strategies used to exploit computing and storage resources. Simulation results obtained in four image modalities and with different models show that VIP is versatile and robust enough to support large simulations. The platform currently has 200 registered users who consumed 33 years of CPU time in 2011.

Model fitting using RANSAC for surgical tool localization in 3-D ultrasound images

Ultrasound guidance is used for many surgical interventions such as biopsy and electrode insertion. We present a method to localize a thin surgical tool such as a biopsy needle or a microelectrode in a 3-D ultrasound image. The proposed method starts with thresholding and model fitting using random sample consensus for robust localization of the axis. Subsequent local optimization refines its position. Two different tool image models are presented: one is simple and fast and the second uses learned a priori information about the tool's voxel intensities and the background. Finally, the tip of the tool is localized by finding an intensity drop along the axis. The simulation study shows that our algorithm can localize the tool at nearly real-time speed, even using a MATLAB implementation, with accuracy better than 1 mm. In an experimental comparison with several alternative localization methods, our method appears to be the fastest and the most robust one. We also show the results on real 3-D ultrasound data from a PVA cryogel phantom, turkey breast, and breast biopsy.

Transverse oscillations for tissue motion estimation

This paper gives an overview of the methods developed for tissue motion estimation using transverse oscillation images (TO). TO images are specific radiofrequency ultrasound images featuring oscillations in both spatial directions. The initial studies on TO were published in the late 1990s. This paper reviews the main ideas and applications behind this motion estimation approach. First the origin and motivation of TO is briefly reviewed. Then the beamforming methods that lead to TO images are given, detailing the receive-only approach and the transmit-and-receive approach using synthetic aperture data. The different medical applications where TO has been used are discussed (blood flow, elastography and echocardiography), showing how it can improve motion estimation. Finally, the future perspectives of TO are outlined.

Phase-based block matching applied to motion estimation with unconventional beamforming strategies

A phase-based block matching method adapted to motion estimation with unconventional beamforming strategies is presented. The unconventional beamforming technique used allows us to obtain 2-D RF images with axial and lateral modulations. Based on these images, we propose a method that uses phase images instead of amplitude images. This way of proceeding allows us to provide an analytical solution to the local displacement estimation so that no minimization of a classical cost function is used for the local estimation. For this reason, the local estimator is directly applied to signals, without the need to process a complex cross-correlation function, as is done with most of the phase shift estimators. In this paper, the method is applied to elastography. Results with simulated data show that a downsampling of axial and lateral modulated signals leads to very little change in the accuracy and in the spatial resolution of the proposed method. For example, for decimation factors of 2 in the axial direction and of 4 in the lateral direction, the mean axial absolute error is 3 mum. The same estimation with original images provides a mean axial error of 0.7 microm. The accuracy of the lateral motion is unchanged in this case. The accuracy of our method with downsampled signals is an important issue in the purpose of a real-time implementation. With experimental data, for the same level of estimation error, classical block matching using the maximum of cross correlation as a local estimator requires images that are 36 times larger (in number of pixels) and consequently a computational time roughly 10 times longer. Our phase block matching is also shown to provide 10 percent less error than a motion estimation method based on seeking the zero of the complex correlation function phase. Finally, it is shown that given the separability of the local estimator that we propose, our method can be applied on both n-D signals and classical RF ultrasound images. The phase block matching method presented was implemented in real time on an ultrasound research scanner.

Analytic estimation of subsample spatial shift using the phases of multidimensional analytic signals

In this correspondence, a method of analytic subsample spatial shift estimation based on an a priori n-D signal model is proposed. The estimation uses the linear phases of n analytic signals defined with the multidimensional Hilbert transform. This estimation proposes: i) an analytic solution to the n-D shift estimation and ii) an estimation without processing complex cross-correlation function or cross-spectra between signals contrary to most phase shift estimators. The method provides better performance in estimating subsample shifts than two classical estimators, one using the maximum of cross-correlation function and the other seeking the zero of the complex correlation function phase. Two delay estimators using the in-phase and quadrature-phase components of signals are also compared to our estimator. Like most estimators using the complex signal phases, the estimator proposed herein presents the advantage of unaltered accuracy when low sampled signals are used. Moreover, we show that this method can be applied to motion tracking with ultrasound images. Thus, included in a block-based motion estimation method and tested with ultrasound data, this estimator provides an analytical solution to the translation estimation problem.

Lateral RF image synthesis using a synthetic aperture imaging technique

The oscillating profile naturally present in ultrasound images has been shown to be extremely valuable in different applications, particularly in motion estimation. Recent studies have shown that it is possible to produce images with transverse oscillations (TOs) based on a specific type of beamforming. However, there is still a great difference between the nature of the lateral oscillations produced with current methods and the axial profile of ultrasound images. In this study, we propose to combine synthetic aperture imaging (synthetic transmit aperture, STA) using a specific beamformer in both transmit mode and receive mode combined with a heterodyning demodulation method to produce lateral radiofrequency signals (LRFs). The aim was to produce lateral signals as close as possible to conventional axial signals, which would make it possible to estimate lateral displacements with the same accuracy as in the axial direction. The feasibility of this approach was validated in simulation and experimentally on an ultrasound research platform, the Ultrasonix RP system. We show that the combination of STA and the heterodyning demodulation can divide the wavelength of the LRF signals by 4 and divide the width of the lateral envelope of the point spread function (PSF) by 2 compared with the previous approaches using beamforming in receive mode only. Finally, we also illustrate the potential of our beamforming for motion estimation compared with previous TO methods.

Parallel integral projection transform for straight electrode localization in 3-D ultrasound images

In surgical practice, small metallic instruments are frequently used to perform various tasks inside the human body. We address the problem of their accurate localization in the tissue. Recent experiments using medical ultrasound have shown that this modality is suitable for real-time visualization of anatomical structures as well as the position of surgical instruments. We propose an image-processing algorithm that permits automatic estimation of the position of a line-segment-shaped object. This method was applied to the localization of a thin metallic electrode in biological tissue. We show that the electrode axis can be found through maximizing the parallel integral projection transform that is a form of the Radon transform. To accelerate this step, hierarchical mesh-grid algorithm is implemented. Once the axis position is known, localization of the electrode tip is performed. The method was tested on simulated images, on ultrasound images of a tissue mimicking phantom containing a metallic electrode, and on real ultrasound images from breast biopsy. The results indicate that the algorithm is robust with respect to variations in electrode position and speckle noise. Localization accuracy is of the order of hundreds of micrometers and is comparable to the ultrasound system axial resolution.

PSF dedicated to estimation of displacement vectors for tissue elasticity imaging with ultrasound

This paper investigates a new approach devoted to displacement vector estimation in ultrasound imaging. The main idea is to adapt the image formation to a given displacement estimation method to increase the precision of the estimation. The displacement is identified as the zero crossing of the phase of the complex cross-correlation between signals extracted from the lateral direction of the ultrasound RF image. For precise displacement estimation, a linearity of the phase slope is needed as well as a high phase slope. Consequently, a particular point spread function (PSF) dedicated to this estimator is designed. This PSF, showing oscillations in the lateral direction, leads to synthesis of lateral RF signals. The estimation is included in a 2-D displacement vector estimation method. The improvement of this approach is evaluated quantitatively by simulation studies. A comparison with a speckle tracking technique is also presented. The lateral oscillations improve both the speckle tracking estimation and our 2-D estimation method. Using our dedicated images, the precision of the estimation is improved by reducing the standard deviation of the lateral displacement error by a factor of 2 for speckle tracking and more than 3 with our method compared to using conventional images. Our method performs 7 times better than speckle tracking. Experimentally, the improvement in the case of a pure lateral translation reaches a factor of 7. Finally, the experimental feasibility of the 2-D displacement vector estimation is demonstrated on data acquired from a Cryogel phantom.

Two-dimensional least-squares estimation for motion tracking in ultrasound elastography

This paper proposes a method of 2-D translations estimation using an a priori signal model. Two analytical signals defined with multidimensional Hilbert transform are considered and shown to have linear phases with respect to the translations to estimate. A least squares estimator (LSE) is then developed to adjust the measured phases of the complex signals to their theoretical forms. Moreover, the LSE provides an analytical solution to the 2-D translation estimation problem. The estimator is then included in a block matching method for motion tracking with ultrasound images. We compared our results with those obtained with a classical sum of absolute differences (SAD) cost function. We show that with our method there is no need of interpolating the images. Thus, for images at the original resolution level, the results obtained with the proposed estimator are largely more accurate than with SAD. Moreover, we show that using SAD on images with resolution five times higher provide roughly the same results as with our method, but the processing time is ten times higher in this case.

Estimation methods for flow imaging with high frequency ultrasound

This article proposes to estimate slow blood flow with high frequency ultrasound imaging. The proposed technique combines 2 methods. First, a statistical method, called Speckle Flow Imaging (SFI) based on the analysis of changes in the speckle pattern along time, gives an index directly related to the total velocity vector. Secondly, a block matching approach estimates the in-plane velocity components. Results on calibrated flow sequences of blood mimicking fluid have shown good agreement with the statistical model. The quantification of flow is achieved with pulsed flow and is also angle independent when the flow is perpendicular to the ultrasound beam. Speckle Tracking has been evaluated on the same data and has shown good estimation of the in-plane velocity vector when the component of velocity perpendicular to the imaging plane is inferior to 1mm/s. The results of these two methods permit the evaluation of the total 3D velocity field and the orthogonal velocity component relative to the imaging plane. This allows the quantification of blood flow (volumetric per time unit across the sequence).