Martino Alessandrini

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 workίow 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.

Creaseg: A free software for the evaluation of image segmentation algorithms based on level-set

This paper describes a free open source software in Matlab (named Creaseg, http://www.creatis.insa-lyon. fr/∼bernard/creaseg) for the evaluation of the performance of different level-set based algorithms in the context of 2D image segmentation. The platform gives access to the implementation of six level-set methods that have been chosen in order to cover a wide range of data attachment terms (contour, region and localized approaches). The software also gives the possibility to compare the performance of the proposed algorithms on any kind of images. The performance can be evaluated either visually, or from similarity measurements between a reference and the results of the segmentation.

Myocardial Motion Estimation From Medical Images Using the Monogenic Signal

We present a method for the analysis of heart motion from medical images. The algorithm exploits monogenic signal theory, recently introduced as an N-dimensional generalization of the analytic signal. The displacement is computed locally by assuming the conservation of the monogenic phase over time. A local affine displacement model is considered to account for typical heart motions as contraction/expansion and shear. A coarse-to-fine B-spline scheme allows a robust and effective computation of the models parameters, and a pyramidal refinement scheme helps to handle large motions. Robustness against noise is increased by replacing the standard point-wise computation of the monogenic orientation with a robust least-squares orientation estimate. Given its general formulation, the algorithm is well suited for images from different modalities, in particular for those cases where time variant changes of local intensity invalidate the standard brightness constancy assumption. This paper evaluates the methods feasibility on two emblematic cases: cardiac tagged magnetic resonance and cardiac ultrasound. In order to quantify the performance of the proposed method, we made use of realistic synthetic sequences from both modalities for which the benchmark motion is known. A comparison is presented with state-of-the-art methods for cardiac motion analysis. On the data considered, these conventional approaches are outperformed by the proposed algorithm. A recent global optical-flow estimation algorithm based on the monogenic curvature tensor is also considered in the comparison. With respect to the latter, the proposed framework provides, along with higher accuracy, superior robustness to noise and a considerably shorter computation time.

Detection of the whole myocardium in 2D-echocardiography for multiple orientations using a geometrically constrained level-set

The segmentation of the myocardium in echocardiographic images is an important task for the diagnosis of heart disease. This task is difficult due to the inherent problems of echographic images (i.e. low contrast, speckle noise, signal dropout, presence of shadows). In this article, we propose a method to segment the whole myocardium (endocardial and epicardial contours) in 2D echographic images. This is achieved using a level-set model constrained by a new shape formulation that allows to model both contours. The novelty of this work also lays in the fact that our framework allows to segment the whole myocardium for the four main views used in clinical routine. The method is validated on a dataset of clinical images and compared with expert segmentation.

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.

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

In this paper we propose 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. Such dictionaries allow for much sparser representations of the signals since they are optimized for a particular class of images such as US images. We will investigate two undersampling patterns of the 3D US imaging: a spatially uniform random acquisition and a line-wise random acquisition. The latter being extremely interesting for 3D imaging: it would indeed allow skipping the acquisition of many lines among the several thousands required in 3D acquisitions, thus, speeding up the whole acquisition process and incrementing the imaging rate. In this study, the dictionary was learned using the K-SVD algorithm on patches extracted from a training dataset constituted of simulated 3D non-log envelope US volumes. Experiments were performed on a testing dataset made of a simulated 3D US log-envelope volume not included in the testing dataset. CS reconstruction was performed by removing 20% to 80% of the original samples according to the two undersampling patterns. Reconstructions using a K-SVD dictionary previously trained dictionary indicate minimal information loss, thus showing the potential of the overcomplete dictionaries.

Simulation of realistic echocardiographic sequences for ground-truth validation of motion estimation

We present a framework for the simulation of realistic cardiac ultrasound sequences. Both the visual aspect and the synthesized motion mimic a real echocardiography sequence used as template. The resulting simulation appears virtually indistinguishable from a real scan. As the true tissue motion is known, these synthetic sequences can provide a trustful benchmark for the analysis of heart motion. This possibility is illustrated by comparing the performance of two well known motion estimation algorithms.

Predictive Deconvolution and Hybrid Feature Selection for Computer-Aided Detection of Prostate Cancer

Computer-aided detection (CAD) schemes are decision making support tools, useful to overcome limitations of problematic clinical procedures. Trans-rectal ultrasound image based CAD would be extremely important to support prostate cancer diagnosis. An effective approach to realize a CAD scheme for this purpose is described in this work, employing a multi-feature kernel classification model based on generalized discriminant analysis. The mutual information of feature value and tissue pathological state is used to select features essential for tissue characterization. System-dependent effects are reduced through predictive deconvolution of the acquired radio-frequency signals. A clinical study, performed on ground truth images from biopsy findings, provides a comparison of the classification model applied before and after deconvolution, showing in the latter case a significant gain in accuracy and area under the receiver operating characteristic curve.

A restoration framework for ultrasonic tissue characterization

Ultrasonic tissue characterization has become an area of intensive research. This procedure generally relies on the analysis of the unprocessed echo signal. Because the ultrasound echo is degraded by the non-ideal system point spread function, a deconvolution step could be employed to provide an estimate of the tissue response that could then be exploited for a more accurate characterization. In medical ultrasound, deconvolution is commonly used to increase diagnostic reliability of ultrasound images by improving their contrast and resolution. Most successful algorithms address deconvolution in a maximum a posteriori estimation framework; this typically leads to the solution of ℓ2-norm or ℓ1-norm constrained optimization problems, depending on the choice of the prior distribution. Although these techniques are sufficient to obtain relevant image visual quality improvements, the obtained reflectivity estimates are, however, not appropriate for classification purposes. In this context, we introduce in this paper a maximum a posteriori deconvolution framework expressly derived to improve tissue characterization. The algorithm overcomes limitations associated with standard techniques by using a nonstandard prior model for the tissue response. We present an evaluation of the algorithm performance using both computer simulations and tissue-mimicking phantoms. These studies reveal increased accuracy in the characterization of media with different properties. A comparison with state-of-the-art Wiener and ℓ1-norm deconvolution techniques attests to the superiority of the proposed algorithm.

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.