Search

Objective: Statistical Parametric Mapping (SPM) is widely used for the quantitative analysis of brain images from ¹⁸F fluorodeoxyglucose positron emission tomography (FDG PET). SPM requires an initial step of spatial normalization to align all images to a standard anatomic model (the template), but this may lead to image distortion and artifacts, especially in cases of marked brain abnormalities. This study aimed at assessing a block-matching (BM) normalization algorithm, where most transformations are not directly computed on the overall brain volume but through small blocks, a principle that is likely to minimize artifacts., Methods: Large and/or small hypometabolic areas were artificially simulated in initially normal FDG PET images to compare the results provided by statistical tests computed after either SPM or BM normalization., Results: Results were enhanced by BM, compared with SPM, with regard to (i) errors in the estimation of large defects volumes (about 2-fold lower) because of a lower image distortion, and (ii) rates of false-positive foci when numerous or extended abnormalities were simulated. These observations were strengthened by analyses of FDG PET examinations from epileptic patients., Conclusions: Results obtained with the BM normalization of brain FDG PET appear more precise and robust than with SPM normalization, especially in cases of numerous or extended abnormalities.

EMPIRE10 (Evaluation of Methods for Pulmonary Image REgistration 2010) is a public platform for fair and meaningful comparison of registration algorithms which are applied to a database of intrapatient thoracic CT image pairs. Evaluation of nonrigid registration techniques is a nontrivial task. This is compounded by the fact that researchers typically test only on their own data, which varies widely. For this reason, reliable assessment and comparison of different registration algorithms has been virtually impossible in the past. In this work we present the results of the launch phase of EMPIRE10, which comprised the comprehensive evaluation and comparison of 20 individual algorithms from leading academic and industrial research groups. All algorithms are applied to the same set of 30 thoracic CT pairs. Algorithm settings and parameters are chosen by researchers expert in the configuration of their own method and the evaluation is independent, using the same criteria for all participants. All results are published on the EMPIRE10 website (http://empire10.isi.uu.nl). The challenge remains ongoing and open to new participants. Full results from 24 algorithms have been published at the time of writing. This paper details the organization of the challenge, the data and evaluation methods and the outcome of the initial launch with 20 algorithms. The gain in knowledge and future work are discussed.

Radiotherapy planning requires accurate delineations of the critical structures. To avoid manual contouring, atlas-based segmentation can be used to get automatic delineations. However, the results strongly depend on the chosen atlas, especially for the head and neck region where the anatomical variability is high. To address this problem, atlases adapted to the patient's anatomy may allow for a better registration, and already showed an improvement in segmentation accuracy. However, building such atlases requires the definition of a criterion to select among a database the images that are the most similar to the patient. Moreover, the inter-expert variability of manual contouring may be high, and therefore bias the segmentation if selecting only one image for each region. To tackle these issues, we present an original method to design a piecewise most similar atlas. Given a query image, we propose an efficient criterion to select for each anatomical region the K most similar images among a database by considering local volume variations possibly induced by the tumor. Then, we present a new approach to combine the K images selected for each region into a piecewise most similar template. Our results obtained with 105 CT images of the head and neck show that our method reduces the over-segmentation seen with an average atlas while being robust to inter-expert manual segmentation variability.

In emission tomography imaging, respiratory motion causes artifacts in lungs and cardiac reconstructed images, which lead to misinterpretations, imprecise diagnosis, impairing of fusion with other modalities, etc. Solutions like respiratory gating, correlated dynamic PET techniques, list-mode data based techniques and others have been tested, which lead to improvements over the spatial activity distribution in lungs lesions, but which have the disadvantages of requiring additional instrumentation or the need of discarding part of the projection data used for reconstruction. The objective of this study is to incorporate respiratory motion compensation directly into the image reconstruction process, without any additional acquisition protocol consideration. To this end, we propose an extension to the maximum likelihood expectation maximization (MLEM) algorithm that includes a respiratory motion model, which takes into account the displacements and volume deformations produced by the respiratory motion during the data acquisition process. We present results from synthetic simulations incorporating real respiratory motion as well as from phantom and patient data.

This manuscript tackles the reconstruction of 3-D volumes via mono-modal registration of series of 2-D biological images (histological sections, autoradiographs, cryosections, etc.). The process of acquiring these images typically induces composite transformations that we model as a number of rigid or affine local transformations embedded in an elastic one. We propose a registration approach closely derived from this model. Given a pair of input images, we first compute a dense similarity field between them with a block matching algorithm. We use as a similarity measure an extension of the classical correlation coefficient that improves the consistency of the field. A hierarchical clustering algorithm then automatically partitions the field into a number of classes from which we extract independent pairs of sub-images. Our clustering algorithm relies on the Earth mover's distribution metric and is additionally guided by robust least-square estimation of the transformations associated with each cluster. Finally, the pairs of sub-images are, independently, affinely registered and a hybrid affine/non-linear interpolation scheme is used to compose the output registered image. We investigate the behavior of our approach on several batches of histological data and discuss its sensitivity to parameters and noise.

Cardiovascular diseases remain the primary cause of death in developed countries. In most cases, exploration of possibly underlying coronary artery pathologies is performed using X-ray coronary angiography. Current clinical routine in coronary angiography is directly conducted in two-dimensional projection images from several static viewing angles. However, for diagnosis and treatment purposes, coronary artery reconstruction is highly suitable. The purpose of this study is to provide physicians with a three-dimensional (3-D) model of coronary arteries, e.g., for absolute 3-D measures for lesion assessment, instead of direct projective measures deduced from the images, which are highly dependent on the viewing angle. In this paper, we propose a novel method to reconstruct coronary arteries from one single rotational X-ray projection sequence. As a side result, we also obtain an estimation of the coronary artery motion. Our method consists of three main consecutive steps: 1) 3-D reconstruction of coronary artery centerlines, including respiratory motion compensation; 2) coronary artery four-dimensional motion computation; 3) 3-D tomographic reconstruction of coronary arteries, involving compensation for respiratory and cardiac motions. We present some experiments on clinical datasets, and the feasibility of a true 3-D Quantitative Coronary Analysis is demonstrated.

In this paper we propose a new two step method to register the liver from two acquisitions. This registration helps experts to make an intra-patient follow-up for hepatic tumors. Firstly, an original and efficient tree matching is applied on different segmentations of the vascular system of a single patient. These vascular systems are segmented from CT-scan images acquired (every six months) during disease treatement, and then modeled as trees. Our method matches common bifurcations and vessels. Secondly, an estimation of liver deformation is computed from the results of the first step. This approach is validated on a large synthetic database containing cases with various deformation and segmentation problems. In each case, after the registration process, the liver recovery is very accurate (around 95%) and the mean localization error for 3D landmarks in liver is small (around 4 mm).

Purpose: Brain tumor radiotherapy requires the volume measurements and the localization of several individual brain structures. Any tool that can assist the physician to perform the delineation would then be of great help. Among segmentation methods, those that are atlas-based are appealing because they are able to segment several structures simultaneously, while preserving the anatomy topology. This study aims to evaluate such a method in a clinical context., Methods and Materials: The brain atlas is made of two three-dimensional (3D) volumes: the first is an artificial 3D magnetic resonance imaging (MRI); the second consists of the segmented structures in this artificial MRI. The elastic registration of the artificial 3D MRI against a patient 3D MRI dataset yields an elastic transformation that can be applied to the labeled image. The elastic transformation is obtained by minimizing the sum of the square differences of the image intensities and derived from the optical flow principle. This automatic delineation (AD) enables the mapping of the segmented structures onto the patient MRI. Parameters of the AD have been optimized on a set of 20 patients. Results are obtained on a series of 6 patients' MRI. A comprehensive validation of the AD has been conducted on performance of atlas-based segmentation in a clinical context with volume, position, sensitivity, and specificity that are compared by a panel of seven experimented physicians for the brain tumor treatments., Results: Expert interobserver volume variability ranged from 16.70 cm(3) to 41.26 cm(3). For patients, the ratio of minimal to maximal volume ranged from 48% to 70%. Median volume varied from 19.47 cm(3) to 27.66 cm(3) and volume of the brainstem calculated by AD varied from 17.75 cm(3) to 24.54 cm(3). Medians of experts ranged, respectively, for sensitivity and specificity, from 0.75 to 0.98 and from 0.85 to 0.99. Median of AD were, respectively, 0.77 and 0.97. Mean of experts ranged, respectively, from 0.78 to 0.97 and from 0.86 to 0.99. Mean of AD were, respectively, 0.76 and 0.97., Conclusions: Results demonstrate that the method is repeatable, provides a good trade-off between accuracy and robustness, and leads to reproducible segmentation and labeling. These results can be improved by enriching the atlas with the rough information of tumor or by using different laws of deformation for the different structures. Qualitative results also suggest that this method can be used for automatic segmentation of other organs such as neck, thorax, abdomen, pelvis, and limbs.

The recent development of 3-D medical imaging devices has given access to the 3-D imaging of in vivo tissues, from an anatomical (MR, CT) or even functional point of view (fMRI, PET, SPECT). However, the resolution of these images is still not sufficient to image anatomical or functional details, that can only be revealed by in vitro imaging (e.g. histology, autoradiography). The deep motivation of this work is the comparison of activations detected by fMRI series analysis to the ones that can be observed in autoradiographic images. The aim of the presented work is to fuse the autoradiographic data with the pre-mortem anatomical MR image, to facilitate the above-mentioned comparison. We show that this fusion can be achieved by using only simple global transformations (rigid and affine), yielding a very satisfactory result.

Abstract: Purpose: The aim of the present study was to quantitatively assess the performance of a block matching-based automatic registration algorithm integrated within the commercial treatment planning system designated ISOgray™ from Dosisoft. The accuracy of the process was evaluated by a phantom study on computed tomography (CT), magnetic resonance (MR) and positron emission tomography (PET) images. Materials and methods: Two phantoms were used to carry out this study: the cylindrical Jaszczak® phantom and the anthropomorphic Liqui-Phil™ Head Phantom (the Phantom Laboratory), containing fillable spheres. External fiducial markers were used to quantify the accuracy of 41 CT/CT, MR/CT and PET/CT automatic registrations with images of the rotated and tilted phantoms. Results: The study first showed that a cylindrical phantom was not adapted for the evaluation of the performance of a block matching-based registration software. Secondly, the Liqui-Phil™ Head Phantom study showed that the algorithm was able to perform automatic registrations of CT/CT and MR/CT images with differences of up to 40° in phantom rotation and of up to 20–30° for PET/CT with accuracy below the image voxel size. Conclusion: The study showed that the block matching-based automatic registration software under investigation was robust, reliable and yielded very satisfactory results. This phantom-based test can be integrated into a periodical quality assurance process and used for any commissioning of image registration software for radiation therapy. [Copyright &y& Elsevier]

Abstract: Intra-subject and inter-subject nonlinear registration based on dense transformations requires the setting of many parameters, mainly for regularization. This task is a major issue, as the global quality of the registration will depend on it. Setting these parameters is, however, very hard, and they may have to be tuned for each patient when processing data acquired by different centers or using different protocols. Thus, we present in this article a method to introduce more coherence in the registration by using fewer degrees of freedom than with a dense registration. This is done by registering the images only on user-defined areas, using a set of affine transformations, which are optimized together in a very efficient manner. Our framework also ensures a smooth and coherent transformation thanks to a new regularization of the affine components. Finally, we ensure an invertible transformation thanks to the Log-Euclidean polyaffine framework. This allows us to get a more robust and very efficient registration method, while obtaining good results as explained below. We performed a qualitative and quantitative evaluation of the obtained results on two applications: first on atlas-based brain segmentation, comparing our results with a dense registration algorithm. Then the second application for which our framework is particularly well suited concerns bone registration in the lower-abdomen area. We obtain in this case a better positioning of the femoral heads than with a dense registration. For both applications, we show a significant improvement in computation time, which is crucial for clinical applications. [Copyright &y& Elsevier]