Motion compensated iterative reconstruction for cardiac X-ray tomography

Within this Ph.D. project, three-dimensional reconstruction methods for moving objects (with a focus on the human heart) from cone-beam X-ray projections using iterative reconstruction algorithms were developed and evaluated. This project was carried in collaboration with the Digital Imaging Group of Philips Research Europe – Hamburg. In cardiac cone-beam computed tomography (CT) a large effort is continuously dedicated to increase scanning speed in order to minimize patient or organ motion during acquisition. In particular, motion causes severe artifacts such as blurring and streaks in tomographic images. While for a large class of applications the current scanning speed is sufficient, in cardiac CT image reconstruction improvements are still required. Whereas it is currently feasible to achieve stable image quality in the resting phases of the cardiac cycle, in the phase of fast motion data acquisition is too slow. A variety of algorithms to reduce or compensate for motion artifacts have been proposed in literature. Most of the correction methods address the calculation of consistent projection data belonging to the same motion state (gated CT reconstruction). Even if gated CT leads to better results, not only with respect to the processing time but also regarding the image quality, it is also limited in its temporal and spatial resolution due to the mechanical movement of the gantry. This can lead to motion blurring, especially in the phases of fast cardiac motion during the RR interval. A motion-compensated reconstruction method for CT can be used to improve the resolution of the reconstructed image and to suppress motion blurring. Iterative techniques are a promising approach to solve this problem, since no direct inversion methods are known for arbitrarily moving objects. In this work, we therefore introduced motion compensation into image reconstruction. In order to determine the unknown cardiac motion, 3 different cardiac-motion estimation methodologies were implemented. Visual and quantitative assessment of the method in a number of applications, including: phantoms; cardiac CT reconstructions; Region of Interest (ROI) CT reconstructions of left and right coronaries of several clinical patients, confirmed its potential

Recalage préservant la topologie des vaisseaux: application à la cardiologie interventionnelle

In percutaneous coronary interventions, integrating into the live fluoroscopic image vessel calcifications and occlusion information that are revealed in the pre-operative Computed Tomography Angiography can greatly improve guidance of the clinician. Fusing pre- and intra-operative information into a single space aims at taking advantage of two complementary modalities and requires a step of registration that must provide good alignment and relevant correspondences between them. Most of the existing 3D/2D vessel registration algorithms do not take into account the particular topology of the vasculature to be matched, resulting into pairings that may be topologically inconsistent along the vasculature.A first contribution consisted in a registration framework dedicated to curve matching, denoted the Iterative Closest Curve (ICC). Its main feature is to preserve the topological consistency along curves by taking advantage of the Frechet distance that not only computes the distance between two curves but also builds ordered pairings along them. A second contribution is a pairing procedure designed for the matching of a vascular tree structure that endorses its particular topology and that can easily take advantage of the ICC-framework. Centerlines of the 3D tree are matched to curves extracted from the 2D vascular graph while preserving the connectivity at 3D bifurcations. The matching criterion used to build the pairings takes into account the geometric distance and the resemblance between curves both based on a global formulation using the Frechet distance.To evaluate our approach we run experiments on a database composed of 63 clinical cases, measuring accuracy on real conditions and robustness with respect to a simulated displacement. Quantitative results have been obtained using two complementary measures that aim at assessing the results both geometrically and topologically, and quantify the resulting alignment error as well as the pairing error. The proposed method exhibits good results both in terms of pairing and alignment and demonstrates to be low sensitive to the rotations to be compensated (up to 30 degrees).Cette thèse s’inscrit dans le cadre de la cardiologie interventionnelle. Intégrer des informations telles que la position des calcifications ainsi que la taille et forme d’une occlusion dans les images fluoroscopiques constituerait un bénéfice pour le praticien. Ces informations, invisibles dans les images rayons-X pendant la procédure, sont présentes au sein du scanner CT préopératoire. La fusion de cette modalité avec la fluoroscopie apporterait une aide précieuse au guidage temps réel des outils interventionnels en bénéficiant des informations fournies par le CT. Cette fusion requiert une étape de recalage qui vise à aligner au mieux les deux modalités et fournir des correspondances pertinentes entre elles. La plupart des algorithmes de recalage 3D/2D de vaisseaux rencontrent des difficultés à construire des appariements anatomiquement pertinents, essentiellement à cause du manque de cohérence topologique le long du réseau vasculaire.Afin de résoudre ce problème, nous proposons dans cette thèse un cadre générique pour le recalage de structures curvilinéaires. L’algorithme qui en découle préserve la structure des courbes appariées. Les artères coronaires pouvant être représentées par un ensemble de courbes arrangées en arbre, nous proposons aussi une procédure d’appariement qui respecte cette structure. Le recalage d’un arbre 3D sur un graphe 2D est ainsi réalisé en assurant la préservation des connectivités aux bifurcations. Le choix de l’appariement est basé sur un critère prenant en compte la distance géométrique ainsi que la ressemblance entre courbes. Ce critère est évalué grâce à une forme modifiée de la distance de Fréchet.Une base de données de 63 cas cliniques a été utilisée à travers différentes expériences afin de prouver la robustesse et la précision de notre approche. Nous avons proposé deux mesures complémentaires visant à quantifier la qualité de l’alignement d’une part et des appariements engendrés d’autre part. La méthode proposée se montre précise pour les alignements de la projection du modèle CT et des artères coronaires observées dans les images angiographiques. De plus, les appariements obtenus sont anatomiquement pertinents et lálgorithme a prouvé sa robustesse face aux perturbations de la position initiale. Nous attribuons cette robustesse à la qualité des appariements construits au fur et à mesure des itérations

Traitement et exploration d'images TDM pour l'évaluation des bioprothèses valvulaires aortiques

Le but de cette étude est d évaluer la faisabilité de l analyse tomodensitométrique 3D des bioprothèses aortiques pour faciliter leur évaluation morphologique durant le suivi et d aider la sélection de cas et améliorer la planification d une procédure valvein-valve. Le challenge était représenté par le rehaussement des feuillets valvulaires, en raison d images très bruitées. Un angio-scanner synchronisé était réalisé chez des patients porteurs d une bioprotèses aortique dégénérée avant réintervention (images in-vivo). Différentes méthodes pour la réduction du bruit étaient proposées. La reconstruction tridimensionnelle des bioprothèses était réalisée en utilisant des méthodes de segmentation de régions par "sticks". Après réopération ces méthodes étaient appliquées aux images scanner des bioprothèses explantées (images ex-vivo) et utilisées comme référence. La réduction du bruit obtenue par le filtre stick modifié montrait meilleurs résultats en rapport signal/bruit en comparaison aux filtres de diffusion anisotropique. Toutes les méthodes de segmentation ont permis une reconstruction 3D des feuillets. L analyse qualitative a montré une bonne concordance entre les images obtenues in-vivo et les altérations des bioprothèses. Les résultats des différentes méthodes étaient comparés par critères volumétriques et discutés. Les bases d'une première approche de visualisation spatio-temporelle d'images TDM 3D+T de la prothèse valvulaire ont été proposés. Elle implique des techniques de rendu volumique et de compensation de mouvement. Son application à la valve native a aussi été envisagée. Les images scanner des bioprothèses aortiques nécessitent un traitement de débruitage et de réduction des artéfacts de façon à permettre le rehaussement des feuillets prothétiques. Les méthodes basées sticks semblent constituer une approche pertinente pour caractériser morphologiquement la dégénérescence des bioprothèses.The aim of the study was to assess the feasibility of CT based 3D analysis of degenerated aortic bioprostheses to make easier their morphological assessment. This could be helpful during regular follow-up and for case selection, improved planning and mapping of valve-in-valve procedure. The challenge was represented by leaflets enhancement because of highly noised CT images. Contrast-enhanced ECG-gated CT scan was performed in patients with degenerated aortic bioprostheses before reoperation (in-vivo images). Different methods for noise reduction were tested and proposed. 3D reconstruction of bioprostheses components was achieved using stick based region segmentation methods. After reoperation, segmentation methods were applied to CT images of the explanted prostheses (exvivo images). Noise reduction obtained by improved stick filter showed best results in terms of signal to noise ratio comparing to anisotropic diffusion filters. All segmentation methods applied to the best phase of in-vivo images allowed 3D bioprosthetic leaflets reconstruction. Explanted bioprostheses CT images were also processed and used as reference. Qualitative analysis revealed a good concordance between the in-vivo images and the bioprostheses alterations. Results from different methods were compared by means of volumetric criteria and discussed. A first approach for spatiotemporal visualization of 3D+T images of valve bioprosthesis has been proposed. Volume rendering and motion compensation techniques were applied to visualize different phases of CT data. Native valve was also considered. ECG-gated CT images of aortic bioprostheses need a preprocessing to reduce noise and artifacts in order to enhance prosthetic leaflets. Stick based methods seems to provide an interesting approach for the morphological characterization of degenerated bioprostheses.RENNES1-Bibl. électronique (352382106) / SudocSudocFranceF

Traitement et exploration d'images TDM pour l'évaluation des bioprothèses valvulaires aortiques

Le but de cette étude est d évaluer la faisabilité de l analyse tomodensitométrique 3D des bioprothèses aortiques pour faciliter leur évaluation morphologique durant le suivi et d aider la sélection de cas et améliorer la planification d une procédure valvein-valve. Le challenge était représenté par le rehaussement des feuillets valvulaires, en raison d images très bruitées. Un angio-scanner synchronisé était réalisé chez des patients porteurs d une bioprotèses aortique dégénérée avant réintervention (images in-vivo). Différentes méthodes pour la réduction du bruit étaient proposées. La reconstruction tridimensionnelle des bioprothèses était réalisée en utilisant des méthodes de segmentation de régions par "sticks". Après réopération ces méthodes étaient appliquées aux images scanner des bioprothèses explantées (images ex-vivo) et utilisées comme référence. La réduction du bruit obtenue par le filtre stick modifié montrait meilleurs résultats en rapport signal/bruit en comparaison aux filtres de diffusion anisotropique. Toutes les méthodes de segmentation ont permis une reconstruction 3D des feuillets. L analyse qualitative a montré une bonne concordance entre les images obtenues in-vivo et les altérations des bioprothèses. Les résultats des différentes méthodes étaient comparés par critères volumétriques et discutés. Les bases d'une première approche de visualisation spatio-temporelle d'images TDM 3D+T de la prothèse valvulaire ont été proposés. Elle implique des techniques de rendu volumique et de compensation de mouvement. Son application à la valve native a aussi été envisagée. Les images scanner des bioprothèses aortiques nécessitent un traitement de débruitage et de réduction des artéfacts de façon à permettre le rehaussement des feuillets prothétiques. Les méthodes basées sticks semblent constituer une approche pertinente pour caractériser morphologiquement la dégénérescence des bioprothèses.The aim of the study was to assess the feasibility of CT based 3D analysis of degenerated aortic bioprostheses to make easier their morphological assessment. This could be helpful during regular follow-up and for case selection, improved planning and mapping of valve-in-valve procedure. The challenge was represented by leaflets enhancement because of highly noised CT images. Contrast-enhanced ECG-gated CT scan was performed in patients with degenerated aortic bioprostheses before reoperation (in-vivo images). Different methods for noise reduction were tested and proposed. 3D reconstruction of bioprostheses components was achieved using stick based region segmentation methods. After reoperation, segmentation methods were applied to CT images of the explanted prostheses (exvivo images). Noise reduction obtained by improved stick filter showed best results in terms of signal to noise ratio comparing to anisotropic diffusion filters. All segmentation methods applied to the best phase of in-vivo images allowed 3D bioprosthetic leaflets reconstruction. Explanted bioprostheses CT images were also processed and used as reference. Qualitative analysis revealed a good concordance between the in-vivo images and the bioprostheses alterations. Results from different methods were compared by means of volumetric criteria and discussed. A first approach for spatiotemporal visualization of 3D+T images of valve bioprosthesis has been proposed. Volume rendering and motion compensation techniques were applied to visualize different phases of CT data. Native valve was also considered. ECG-gated CT images of aortic bioprostheses need a preprocessing to reduce noise and artifacts in order to enhance prosthetic leaflets. Stick based methods seems to provide an interesting approach for the morphological characterization of degenerated bioprostheses.RENNES1-Bibl. électronique (352382106) / SudocSudocFranceF

Analysis of contrast-enhanced medical images.

Early detection of human organ diseases is of great importance for the accurate diagnosis and institution of appropriate therapies. This can potentially prevent progression to end-stage disease by detecting precursors that evaluate organ functionality. In addition, it also assists the clinicians for therapy evaluation, tracking diseases progression, and surgery operations. Advances in functional and contrast-enhanced (CE) medical images enabled accurate noninvasive evaluation of organ functionality due to their ability to provide superior anatomical and functional information about the tissue-of-interest. The main objective of this dissertation is to develop a computer-aided diagnostic (CAD) system for analyzing complex data from CE magnetic resonance imaging (MRI). The developed CAD system has been tested in three case studies: (i) early detection of acute renal transplant rejection, (ii) evaluation of myocardial perfusion in patients with ischemic heart disease after heart attack; and (iii), early detection of prostate cancer. However, developing a noninvasive CAD system for the analysis of CE medical images is subject to multiple challenges, including, but are not limited to, image noise and inhomogeneity, nonlinear signal intensity changes of the images over the time course of data acquisition, appearances and shape changes (deformations) of the organ-of-interest during data acquisition, determination of the best features (indexes) that describe the perfusion of a contrast agent (CA) into the tissue. To address these challenges, this dissertation focuses on building new mathematical models and learning techniques that facilitate accurate analysis of CAs perfusion in living organs and include: (i) accurate mathematical models for the segmentation of the object-of-interest, which integrate object shape and appearance features in terms of pixel/voxel-wise image intensities and their spatial interactions; (ii) motion correction techniques that combine both global and local models, which exploit geometric features, rather than image intensities to avoid problems associated with nonlinear intensity variations of the CE images; (iii) fusion of multiple features using the genetic algorithm. The proposed techniques have been integrated into CAD systems that have been tested in, but not limited to, three clinical studies. First, a noninvasive CAD system is proposed for the early and accurate diagnosis of acute renal transplant rejection using dynamic contrast-enhanced MRI (DCE-MRI). Acute rejection–the immunological response of the human immune system to a foreign kidney–is the most sever cause of renal dysfunction among other diagnostic possibilities, including acute tubular necrosis and immune drug toxicity. In the U.S., approximately 17,736 renal transplants are performed annually, and given the limited number of donors, transplanted kidney salvage is an important medical concern. Thus far, biopsy remains the gold standard for the assessment of renal transplant dysfunction, but only as the last resort because of its invasive nature, high cost, and potential morbidity rates. The diagnostic results of the proposed CAD system, based on the analysis of 50 independent in-vivo cases were 96% with a 95% confidence interval. These results clearly demonstrate the promise of the proposed image-based diagnostic CAD system as a supplement to the current technologies, such as nuclear imaging and ultrasonography, to determine the type of kidney dysfunction. Second, a comprehensive CAD system is developed for the characterization of myocardial perfusion and clinical status in heart failure and novel myoregeneration therapy using cardiac first-pass MRI (FP-MRI). Heart failure is considered the most important cause of morbidity and mortality in cardiovascular disease, which affects approximately 6 million U.S. patients annually. Ischemic heart disease is considered the most common underlying cause of heart failure. Therefore, the detection of the heart failure in its earliest forms is essential to prevent its relentless progression to premature death. While current medical studies focus on detecting pathological tissue and assessing contractile function of the diseased heart, this dissertation address the key issue of the effects of the myoregeneration therapy on the associated blood nutrient supply. Quantitative and qualitative assessment in a cohort of 24 perfusion data sets demonstrated the ability of the proposed framework to reveal regional perfusion improvements with therapy, and transmural perfusion differences across the myocardial wall; thus, it can aid in follow-up on treatment for patients undergoing the myoregeneration therapy. Finally, an image-based CAD system for early detection of prostate cancer using DCE-MRI is introduced. Prostate cancer is the most frequently diagnosed malignancy among men and remains the second leading cause of cancer-related death in the USA with more than 238,000 new cases and a mortality rate of about 30,000 in 2013. Therefore, early diagnosis of prostate cancer can improve the effectiveness of treatment and increase the patient’s chance of survival. Currently, needle biopsy is the gold standard for the diagnosis of prostate cancer. However, it is an invasive procedure with high costs and potential morbidity rates. Additionally, it has a higher possibility of producing false positive diagnosis due to relatively small needle biopsy samples. Application of the proposed CAD yield promising results in a cohort of 30 patients that would, in the near future, represent a supplement of the current technologies to determine prostate cancer type. The developed techniques have been compared to the state-of-the-art methods and demonstrated higher accuracy as shown in this dissertation. The proposed models (higher-order spatial interaction models, shape models, motion correction models, and perfusion analysis models) can be used in many of today’s CAD applications for early detection of a variety of diseases and medical conditions, and are expected to notably amplify the accuracy of CAD decisions based on the automated analysis of CE images