Model-based Reconstruction of Myocardial Perfusion SPECT and PET Images

Myocardial perfusion imaging is an important noninvasive tool in the diagnosis and prognosis of coronary artery disease, the leading cause of death in the United States. Electrocardiographically (ECG) gated acquisition allows combined evaluation of perfusion and left ventricular function within a single study. However, the accuracy of perfusion quantification and functional analysis is reduced by a number of image degrading factors. Particularly, partial volume effects (PVE) resulting from finite spatial resolution cause activity spillover between tissue classes and blur region boundaries. High-resolution anatomical images, such as contrast CT or MRI, can be used for partial volume compensation (PVC), but they are generally not available in clinical practice. The objective of this research is to develop and evaluate a model-based reconstruction method for emission computed tomography applied to myocardial perfusion imaging to improve perfusion quantification and functional assessment without the presence of anatomical images. The idea is to model the left ventricle (LV) using a geometry model and an activity distribution model instead of modeling them using voxels. The geometry model parameterizes the endocardial and epicardial surfaces using a set of rays originating from the long axis of LV. The rays sample the surfaces cylindrically in the basal and mid-ventricular regions and spherically in the apex. The surfaces are obtained by interpolating the corresponding intersection points with the rays using a cubic-spline function. The activity distribution model divides the myocardium into segments similar to those used in standardized myocardial quantitative analysis. The model assumes uniform activity concentrations in the segments as well as the blood pool and body background. The method estimates the parameters of the geometry and activity models instead of the intensities of all voxels, which greatly reduces the number of unknowns to be estimated. The goal of model-based reconstruction method was to estimate the parameters that give the best match between the image generated by the model and the measured data. The input image is contaminated by noise, so the metric for the goodness of fit was a statistical criterion based on the likelihood, i.e., the probability that the image resulted from the given set of model parameters. The image generated by the model includes the effects of resolution and other image degrading factors. A shape constraint was also incorporated into the objective function to regularize the ill-posed reconstruction problem and increase the robustness to perfusion defects and noise. The model parameters were optimized by seeking the maximum of the objective function using a group-wise alternating scheme along with dedicated initial parameter estimation. The hypothesis underlying the work is that resulting geometry parameters produce an accurate segmentation of the LV while the activity parameters are PVE compensated representations of the true activities in the segments. The proposed method integrates prior knowledge about the targeted object and the imaging system into one framework and allows simultaneous LV segmentation and PVC. In the evaluation with simulated myocardial perfusion SPECT images, it improved accuracy and precision in delineating the myocardium in comparison with typical segmentation methods. In addition, it recovered the myocardial activity more effectively compared to deconvolution-based PVC, which also does not require coregistered anatomical images to define regions of interest

Cardiologie nucléaire du 21ième siècle : nouveautés et réalités

Les maladies cardio-vasculaires demeurent une cause majeure de mortalité et morbidité dans les sociétés développées. La recherche de déterminants prédictifs d’évènements vasculaires représente toujours un enjeu d’actualité face aux coûts croissants des dépenses reliées aux soins médicaux et à l’élargissement des populations concernées, notamment face à l’occidentalisation des pays émergeants comme l’Inde, le Brésil et la Chine. La cardiologie nucléaire occupe depuis trente ans, une place essentielle dans l’arsenal des méthodes diagnostiques et pronostiques des cardiopathies. De plus, de nouvelles percées permettront de dépister d’une façon plus précoce et précise, la maladie athérosclérotique cardiaque et périphérique chez les populations atteintes ainsi qu’en prévention primaire. Nous présenterons dans cette thèse, deux approches nouvelles de la cardiologie nucléaire. La dysfonction endothéliale est considérée comme le signal pathologique le plus précoce de l’athérosclérose. Les facteurs de risques cardiovasculaires traditionnels atteignent la fonction endothéliale et peuvent initier le processus d’athérosclérose même en l’absence de lésion endothéliale physique. La quantification de la fonction endothéliale coronarienne comporte donc un intérêt certain comme biomarqueur précoce de la maladie coronarienne. La pléthysmographie isotopique, méthodologie développée lors de ce cycle d’étude, permet de quantifier la fonction endothéliale périphérique, cette dernière étant corrélée à la fonction endothéliale coronarienne. Cette méthodologie est démontrée dans le premier manuscrit (Harel et. al., Physiol Meas., 2007). L’utilisation d’un radiomarquage des érythrocytes permet la mesure du flot artériel au niveau du membre supérieur pendant la réalisation d’une hyperémie réactive locale. Cette nouvelle procédure a été validée en comparaison à la pléthysmographie par jauge de contrainte sur une cohorte de 26 patients. Elle a démontré une excellente reproductibilité (coefficient de corrélation intra-classe = 0.89). De plus, la mesure du flot artérielle pendant la réaction hyperémique corrélait avec les mesure réalisées par la méthode de référence (r=0.87). Le deuxième manuscrit expose les bases de la spectroscopie infrarouge comme méthodologie de mesure du flot artériel et quantification de la réaction hyperémique (Harel et. al., Physiol Meas., 2008). Cette étude utilisa un protocole de triples mesures simultanées à l’aide de la pléthysmographie par jauge de contrainte, radio-isotopique et par spectroscopie infrarouge. La technique par spectroscopie fut démontrée précise et reproductible quant à la mesure des flots artériels au niveau de l’avant-bras. Cette nouvelle procédure a présenté des avantages indéniables quant à la diminution d’artéfact et à sa facilité d’utilisation. Le second volet de ma thèse porte sur l’analyse du synchronisme de contraction cardiaque. En effet, plus de 30% des patients recevant une thérapie de resynchronisation ne démontre pas d’amélioration clinique. De plus, ce taux de non-réponse est encore plus élevé lors de l’utilisation de critères morphologiques de réponse à la resynchronisation (réduction du volume télésystolique). Il existe donc un besoin urgent de développer une méthodologie de mesure fiable et précise de la dynamique cardiaque. Le troisième manuscrit expose les bases d’une nouvelle technique radio-isotopique permettant la quantification de la fraction d’éjection du ventricule gauche (Harel et. al. J Nucl Cardiol., 2007). L’étude portant sur 202 patients a démontré une excellente corrélation (r=0.84) avec la méthode de référence (ventriculographie planaire). La comparaison avec le logiciel QBS (Cedar-Sinai) démontrait un écart type du biais inférieur (7.44% vs 9.36%). De plus, le biais dans la mesure ne démontrait pas de corrélation avec la magnitude du paramètre pour notre méthodologie, contrairement au logiciel alterne. Le quatrième manuscrit portait sur la quantification de l’asynchronisme intra-ventriculaire gauche (Harel et. al. J Nucl Cardiol, 2008). Un nouveau paramètre tridimensionnel (CHI: contraction homogeneity index) (médiane 73.8% ; IQ 58.7% - 84.9%) permis d’intégrer les composantes d’amplitude et du synchronisme de la contraction ventriculaire. La validation de ce paramètre fut effectuée par comparaison avec la déviation standard de l’histogramme de phase (SDΦ) (médiane 28.2º ; IQ 17.5º - 46.8º) obtenu par la ventriculographie planaire lors d’une étude portant sur 235 patients. Ces quatre manuscrits, déjà publiés dans la littérature scientifique spécialisée, résument une fraction des travaux de recherche que nous avons effectués durant les trois dernières années. Ces travaux s’inscrivent dans deux axes majeurs de développement de la cardiologie du 21ième siècle.Cardiovascular diseases remain a major cause of mortality and morbidity in developed countries. The search for predictive determinants of vascular events represents a relevant and timely goal, considering the increasing costs of medical care and the progress in developing countries such as India, Brazil and China. Nuclear cardiology has, for 30 years, played an essential role in the diagnosis and prognosis of various cardiac and vascular diseases. Moreover, new developments will allow earlier and more specific detection of cardiac and peripheral atherosclerosis disease in affected individuals and in primary prevention. In this thesis, we will focus on advances in two major themes of nuclear cardiology. Endothelial dysfunction is regarded as the earliest pathological markers of atherosclerosis. Traditional cardiovascular risks factors impair endothelial function and can initiate the atherosclerosis process, even in the absence of overt endothelial disruption. Quantification of coronary endothelial function is, therefore, of considerable interest as an early biomarker for coronary disease. The radionuclide plethysmography methodology developed during the course of my doctoral studies allows the quantification of peripheral endothelial function, which has been correlated with coronary endothelial function. This methodology is detailed in the first manuscript (Harel et. al., Physiol Meas., 2007). The use of red blood cell radio-labeling permits arterial flow to be measured in the upper limb during local reactive hyperemia. This new procedure was validated against strain gauge plethysmography in a cohort of 26 patients with excellent reproducibility (intraclass coefficient of correlation = 0.89). Moreover, the arterial measurements of flow during the hyperemic reaction correlated well with the reference method (r=0.87). The second manuscript exposes the basis of infrared spectroscopy as a method for measuring arterial flow and quantifying the hyperemic reaction (Harel et. al., Physiol Meas., 2008). The study protocol consisted of simultaneous measurements by strain gauge, radionuclide and infrared spectroscopy plethysmography. The spectroscopy technique was shown to be precise and reproducible for forearm measurement of arterial blood flow. This novel procedure came major advantages in reducing artifacts and in its ease of use. The second axis of my thesis relates to the analysis of cardiac contraction synchrony. Indeed, more than 30% of patients receiving resynchronization therapy do not show clinical improvement. Moreover, this non-response rate is even higher if we consider morphological criteria of resynchronization (end-systolic volume reduction). There is therefore, an urgent need to improve a methodology to reliably and precisely measure cardiac dynamics so as to identify and monitor potential responders. The third manuscript exposes the basis of a new radionuclide technique to quantify left ventricle ejection fraction (Harel et. al. J Nucl Cardiol., 2007). The study of 202 patients showed an excellent correlation (r=0.84) with the reference method (planar ventriculography). The comparison with QBS software (Cedar-Sinai), showed a lower standard deviation of bias (7.44% vs 9.36%). Moreover, unlike the alternative software, the bias did not correlate with the magnitude of the ejection fraction. The fourth manuscript relates to the quantification of the left intra-ventricular synchronism (Harel et. al. J Nucl Cardiol, 2008). A new three-dimensional parameter (CHI: contraction homogeneity index) (median 73.8%; IQ 58.7% - 84.9%) was defined to allow the integration of amplitude and synchrony components of ventricular contraction. Validation of this parameter was undertaken out by comparing the standard deviation of the histogram of phase (SDΦ) (median 28.2º; IQ 17.5º- 46.8º) obtained by planar ventriculography in a study of 235 patients. These four manuscripts, already published in the specialized scientific literature, summarize a fraction of the research tasks that we have carried out during the three last years, representing two major axes of nuclear cardiology advancement in the 21st century

In Vivo Quantitative Assessment of Myocardial Structure, Function, Perfusion and Viability Using Cardiac Micro-computed Tomography

The use of Micro-Computed Tomography (MicroCT) for in vivo studies of small animals as models of human disease has risen tremendously due to the fact that MicroCT provides quantitative high-resolution three-dimensional (3D) anatomical data non-destructively and longitudinally. Most importantly, with the development of a novel preclinical iodinated contrast agent called eXIA160, functional and metabolic assessment of the heart became possible. However, prior to the advent of commercial MicroCT scanners equipped with X-ray flat-panel detector technology and easy-to-use cardio-respiratory gating, preclinical studies of cardiovascular disease (CVD) in small animals required a MicroCT technologist with advanced skills, and thus were impractical for widespread implementation. The goal of this work is to provide a practical guide to the use of the high-speed Quantum FX MicroCT system for comprehensive determination of myocardial global and regional function along with assessment of myocardial perfusion, metabolism and viability in healthy mice and in a cardiac ischemia mouse model induced by permanent occlusion of the left anterior descending coronary artery (LAD)

Análise funcional do ventrículo esquerdo em angio-TC coronária

Doutoramento em Engenharia InformáticaCoronary CT angiography is widely used in clinical practice for the assessment of coronary artery disease. Several studies have shown that the same exam can also be used to assess left ventricle (LV) function. LV function is usually evaluated using just the data from end-systolic and end-diastolic phases even though coronary CT angiography (CTA) provides data concerning multiple cardiac phases, along the cardiac cycle. This unused wealth of data, mostly due to its complexity and the lack of proper tools, has still to be explored in order to assess if further insight is possible regarding regional LV functional analysis. Furthermore, different parameters can be computed to characterize LV function and while some are well known by clinicians others still need to be evaluated concerning their value in clinical scenarios. The work presented in this thesis covers two steps towards extended use of CTA data: LV segmentation and functional analysis. A new semi-automatic segmentation method is presented to obtain LV data for all cardiac phases available in a CTA exam and a 3D editing tool was designed to allow users to fine tune the segmentations. Regarding segmentation evaluation, a methodology is proposed in order to help choose the similarity metrics to be used to compare segmentations. This methodology allows the detection of redundant measures that can be discarded. The evaluation was performed with the help of three experienced radiographers yielding low intraand inter-observer variability. In order to allow exploring the segmented data, several parameters characterizing global and regional LV function are computed for the available cardiac phases. The data thus obtained is shown using a set of visualizations allowing synchronized visual exploration. The main purpose is to provide means for clinicians to explore the data and gather insight over their meaning, as well as their correlation with each other and with diagnosis outcomes. Finally, an interactive method is proposed to help clinicians assess myocardial perfusion by providing automatic assignment of lesions, detected by clinicians, to a myocardial segment. This new approach has obtained positive feedback from clinicians and is not only an improvement over their current assessment method but also an important first step towards systematic validation of automatic myocardial perfusion assessment measures.A angiografia coronária por TC (angio-TC) é prática clínica corrente para a avaliação de doença coronária. Alguns estudos mostram que é também possível utilizar o exame de angio-TC para avaliar a função do ventrículo esquerdo (VE). A função ventricular esquerda (FVE) é normalmente avaliada considerando as fases de fim de sístole e de fim de diástole, apesar de a angio-TC proporcionar dados relativos a diferentes fases distribuídas ao longo do ciclo cardíaco. Estes dados não considerados, devido à sua complexidade e à falta de ferramentas apropriadas para o efeito, têm ainda de ser explorados para que se perceba se possibilitam uma melhor compreensão da FVE. Para além disso, podem ser calculados diferentes parâmetros para caracterizar a FVE e, enquanto alguns são bem conhecidos dos médicos, outros requerem ainda uma avaliação do seu valor clínico. No âmbito de uma utilização alargada dos dados proporcionados pelos angio- TC, este trabalho apresenta contributos ao nível da segmentação do VE e da sua análise funcional. É proposto um método semi-automático para a segmentação do VE de forma a obter dados para as diferentes fases cardíacas presentes no exame de angio- TC. Foi também desenvolvida uma ferramenta de edição 3D que permite aos utilizadores a correcção das segmentações assim obtidas. Para a avaliação do método de segmentação apresentado foi proposta uma metodologia que permite a detecção de medidas de similaridade redundantes, a usar no âmbito da avaliação para comparação entre segmentações, para que tais medidas redundantes possam ser descartadas. A avaliação foi executada com a colaboração de três técnicos de radiologia experientes, tendo-se verificado uma baixa variabilidade intra- e inter-observador. De forma a permitir explorar os dados segmentados, foram calculados vários parâmetros para caracterização global e regional da FVE, para as diversas fases cardíacas disponíveis. Os resultados assim obtidos são apresentados usando um conjunto de visualizações que permitem uma exploração visual sincronizada dos mesmos. O principal objectivo é proporcionar ao médico a exploração dos resultados obtidos para os diferentes parâmetros, de modo a que este tenha uma compreensão acrescida sobre o seu significado clínico, assim como sobre a correlação existente entre diferentes parâmetros e entre estes e o diagnóstico. Finalmente, foi proposto um método interactivo para ajudar os médicos durante a avaliação da perfusão do miocárdio, que atribui automaticamente as lesões detectadas pelo médico ao respectivo segmento do miocárdio. Este novo método obteve uma boa receptividade e constitui não só uma melhoria em relação ao método tradicional mas é também um primeiro passo para a validação sistemática de medidas automáticas da perfusão do miocárdio

The Assessment of left ventricular Function in MRI using the detection of myocardial borders and optical flow approaches: A Review

The evaluation of left ventricular wall motion in Magnetic Resonance Imaging (MRI) clinical practice is based on a visual assessment of cine-MRI sequences. In fact, clinical interpreters (radiologists) proceed with a global visual evaluation of multiple cine-MRI sequences acquired in the three standard views. In addition, some functional parameters are quantified following a manual or a semi-automatic contouring of the myocardial borders. Although these parameters give information about the functional state of the left ventricle, they are not able to provide the location and the extent of wall motion abnormalities, which are associated with many cardiovascular diseases. In the past years, several approaches were developed to overcome the limitations of the classical evaluation techniques of left ventricular function. The aim of this article is to present an overview of the different methods and to summarize the relevant techniques based on myocardial contour detection and optical flow for regional assessment of left ventricular abnormalities

Quantification of ventricular mechanical dyssynchrony under stress

L'évaluation de l'asynchronisme mécanique ventriculaire sous stress a soulevé une attention importante en tant que facteur prédictif de la réponse au traitement de resynchronisation cardiaque (CRT). De plus, il semble exister une relation significative entre le devenir du patient et la présence d’asynchronisme au repos. Plusieurs méthodes échocardiographiques peuvent être utilisées pour évaluer l’asynchronisme. Cependant, parmi toutes les différentes méthodologies ou index existant dans ce domaine, aucun critère ne fait l’unanimité. Cette thèse étudie l'importance des techniques d'imagerie nucléaire dans le cadre de l’évaluation de l’asynchronisme cardiaque induit par le stress en utilisant trois différents modèles canins expérimentaux. Le premier chapitre vise à examiner les effets du stress sur le synchronisme de la contraction du ventricule gauche (VG) en utilisant l'imagerie synchronisée de perfusion myocardique dans une cohorte canine normale. Le stress a été induit par différents niveaux d’infusion de dobutamine sur six sujets sains. Les paramètres hémodynamiques et l’asynchronisme ont été évalués par des mesures de pressions ventriculaires. L'analyse de phase sur l’imagerie s’est effectuée en utilisant un logiciel commercialement disponible (QGS) et un logiciel interne (MHI4MPI), basée sur le déplacement et l’épaississement des parois ventriculaires. L’augmentation de la concentration de dobutamine a démontré une amélioration de la capacité fonctionnelle et une réduction de l’asynchronisme ventriculaire. L’analyse de l’asynchronisme calculée à partir de l’épaississement de la paroi semble plus robuste et plus sensible que l’utilisation du déplacement des parois. (Salimian et. al., J Nucl Cardiol., 2014) Le second chapitre étudie les différents paramètres d’asynchronisme au repos et à différents niveaux de stress dans un modèle de cardiomyopathie dilatée et à QRS étroit. Ce modèle a été créé sur dix chiens par tachycardie via stimulation de l'apex du ventricule droit pendant 3-4 semaines, permettant d’atteindre une fraction d'éjection cible de 35% ou moins. Le stress a ensuite été induit par une perfusion de dobutamine jusqu'à un maximum de 20 μg/kg/min. Les données hémodynamiques et l’asynchronisme ont été analysés par des mesures de pression ventriculaire et l’analyse de l’imagerie dynamique du compartiment sanguin. L’importante variabilité individuelle des sujets inclus dans notre cohorte empêche toute conclusion définitive sur la mesure de l’asynchronisme interventriculaire. Cependant, les différents niveaux de stress, même dans des intervalles rapprochés, ont démontré un effet significatif sur les paramètres hémodynamiques et l’asynchronisme. (Salimian et. al., J Nucl Cardiol., 2015) La troisième section vise à déterminer si l’estimation du mode de stimulation optimal effectuée au repos demeure le choix optimal lorsque le niveau d’activité cardiaque s’intensifie pour des sujets avec bloc auriculo-ventriculaire (AV) et fonction ventriculaire normale. Cinq chiens ont été soumis à une ablation du nœud AV et des sondes de stimulation ont été insérées dans l'oreillette droite pour la détection, l’apex du ventricule droit (VD) et une veine postérolatérale du VG pour la stimulation. Cinq modes de stimulation ont été utilisés : LV pur, biventriculaire (BiV) avec pré-activation de 20 ms du LV (LVRV20), BiV pur, BiV avec pré-activation de 20 ms du VD (RVLV20), VD pur. Des niveaux jusqu’à 20 μg/kg/min de dobutamine ont été atteints. Le stress a modifié l’étendue de l’asynchronisme de base et ce, pour tous les modes de stimulation. De plus, les effets physiologiques intrinsèques du stress permettent une évaluation plus précise de l’asynchronisme ventriculaire, diminuant la variabilité inter-sujet. Le mode de stimulation LVRV20 semble le mode optimal dans ce modèle, supportant l’utilisation de la stimulation bi-ventriculaire.Assessment of ventricular mechanical dyssynchrony (MD) under stress has attracted a large amount of attention as a stronger predictor of response to cardiac resynchronization therapy (CRT) and as a parameter whose variation bears a greater relationship to clinical outcomes than resting-MD either in CRT candidates or another subset of patients. Several echocardiographic methods can be used to assess stress-MD. However, no standardized approach is currently used to explore stress-induced variations in inter- and intraventricular MD. This dissertation studies the importance of nuclear imaging techniques in assessing stress-induced MD variations by providing three different experimental canine models. The first chapter sought to examine the impacts of stress on the left ventricular (LV) synchrony with phase analysis of gated SPECT myocardial perfusion imaging (GMPS) within a normal canine cohort. Stress was induced by different levels of dobutamine infusion in six healthy subjects. Hemodynamic and LV MD parameters were assessed by LV pressure measurements and phase analysis of GMPS using commercially available QGS software and in-house MHI4MPI software with thickening- and displacement-based methodology. The increase of dobutamine level was shown to be in accordance with the improvement of LV functional capacity and reduction of MD parameters. MD analysis based on wall thickening was more robust and sensitive than the global wall displacement. (Salimian et. al., J Nucl Cardiol., 2014) The second chapter investigated the range of difference in inter- and intraventricular MD parameters from rest to various levels of stress in a dilated cardiomyopathy (DCM) and narrow QRS complex model. Ten large dogs were submitted to tachycardia-induced DCM by pacing the right ventricular apex for 3-4 weeks to reach a target ejection fraction of 35% or less. Stress was then induced by infusion of dobutamine up to a maximum of 20 μg/kg/min. Hemodynamic and MD data were analyzed by LV pressure measurements and gated-blood pool SPECT (GBPS) imaging. Individual differences in the magnitude and pattern of change in the various levels of stress precluded any definitive conclusion about interventricular MD. However, different levels of stress, even in close intervals, showed a significant positive impact on hemodynamic and intraventricular MD parameters. (Salimian et. al., J Nucl Cardiol., 2015) The third chapter sought to examine if the optimal pacing mode at rest could be the best one during the maximum stress level in terms of MD parameters in subjects with an atrioventricular (AV) block and normal function. Five dogs were submitted to AV node ablation and pacing leads were placed in the right atrium for sensing, in right ventricular (RV) apex, and in posterolateral LV vein for pacing in five modes of LV, biventricular (BiV) with 20 ms of LV pre-activation (LVRV20), BiV, BiV with 20 ms of RV pre-activation (RVLV20) and RV pacing. Stress was induced by dobutamine infusion up to a maximum of 20 μg/kg/min. Data analyses were the same as chapter one. Dobutamine stress changed the extent of resting-LV MD at all pacing modes. Intrinsic physiologic effects of stress resulted in more accurate MD assessment with lesser variability in subjects who underwent pacing. LVRV20 was the preferred site of stimulation in this model rather than single-site pacing

State of the Art of Level Set Methods in Segmentation and Registration of Medical Imaging Modalities

Segmentation of medical images is an important step in various applications such as visualization, quantitative analysis and image-guided surgery. Numerous segmentation methods have been developed in the past two decades for extraction of organ contours on medical images. Low-level segmentation methods, such as pixel-based clustering, region growing, and filter-based edge detection, require additional pre-processing and post-processing as well as considerable amounts of expert intervention or information of the objects of interest. Furthermore the subsequent analysis of segmented objects is hampered by the primitive, pixel or voxel level representations from those region-based segmentation. Deformable models, on the other hand, provide an explicit representation of the boundary and the shape of the object. They combine several desirable features such as inherent connectivity and smoothness, which counteract noise and boundary irregularities, as well as the ability to incorporate knowledge about the object of interest. However, parametric deformable models have two main limitations. First, in situations where the initial model and desired object boundary differ greatly in size and shape, the model must be re-parameterized dynamically to faithfully recover the object boundary. The second limitation is that it has difficulty dealing with topological adaptation such as splitting or merging model parts, a useful property for recovering either multiple objects or objects with unknown topology. This difficulty is caused by the fact that a new parameterization must be constructed whenever topology change occurs, which requires sophisticated schemes. Level set deformable models, also referred to as geometric deformable models, provide an elegant solution to address the primary limitations of parametric deformable models. These methods have drawn a great deal of attention since their introduction in 1988. Advantages of the contour implicit formulation of the deformable model over parametric formulation include: (1) no parameterization of the contour, (2) topological flexibility, (3) good numerical stability, (4) straightforward extension of the 2D formulation to n-D. Recent reviews on the subject include papers from Suri. In this chapter we give a general overview of the level set segmentation methods with emphasize on new frameworks recently introduced in the context of medical imaging problems. We then introduce novel approaches that aim at combining segmentation and registration in a level set formulation. Finally we review a selective set of clinical works with detailed validation of the level set methods for several clinical applications