Les systèmes environnementaux présentent une grande variabilité spatio-temporelle. En raison des influences extérieures et de la stochasticité introduite par le processus de reproduction, la dynamique des populations est également caractérisée par une grande variabilité. La recherche de lois universelles en écologie implique souvent une distribution de la forme loi de puissance, ces distributions peuvent survenir dans la dynamique des populations ou dans les patrons spatiaux de distribution d'abondance et de richesse spécifique. En prenant les polychètes, groupe colonisant un grand nombre d'habitats des sédiments meubles et durs et considérés comme indicateurs des principales conditions qui contrôlent la structure et le fonctionnement des communautés macrobenthiques, nous essayons d’identifier les changements spatio-temporels de la biodiversité de ce groupe. A partir de séries biologiques à long terme de trois communautés de fonds meubles, nous étudions la dynamique des populations des polychètes, nous caractérisons des événements extrêmes pour des données d’abondance et nous introduisons des méthodes de quantification pour des séries biologiques très intermittentes. Ensuite, nous abordons la distribution spatiale des espèces de polychètes ayant comme objectifs d'identifier les patrons latitudinaux, longitudinaux et bathymétriques du plateau continental de l’Atlantique nord-est européen; et de tester l’existence de patrons généraux, voir même universels, caractérisant la biodiversité (i.e. l'augmentation de la diversité avec l’aire échantillonnée, sa décroissance de l’équateur vers les pôles et l'augmentation de la richesse en espèces avec l’abondance totale d’individus). / One of the key features of environmental field studies is their high variability at many different time and space scales. Because of these external influences and of the stochasticity introduced by the reproduction, population dynamics are also characterised by high variability over time and space. The search for universal scaling laws in ecology often involves considering a form of power-law distribution, power laws can emerge in population dynamics or in patterns of abundance, distribution, and richness. Using the polychaetes, group that colonises a large range of soft and hard marine sediment habitats, from intertidal to hadal zones, and are considered to be good surrogates to identify the main environmental conditions that control the structure and functioning of benthic communities, we try to identify the spatiotemporal changes in biodiversity for this characteristic benthic group. First, we discuss the dynamics of polychaete populations. Based on long-term series of three soft-bottom communities, we study the dynamics of polychaete populations using different statistical techniques; we characterise extreme events in abundance data and we show how to apply some quantification methods to highly erratic and intermittent biological series. Then, we discuss the spatial distribution of polychaete species aiming to: identify latitudinal, longitudinal and bathymetric patterns on the European northeast Atlantic continental shelf; and test the existence of general, perhaps universal, patterns for characterising biodiversity i.e. increasing diversity with sampled area, its decay from the equator to the poles and the increase in richness with the total abundance of individuals.

Cette thèse porte sur la modélisation 0D/1D de la combustion Diesel dans les moteurs récents. L'objectif est d'augmenter la précision des modèles tout en limitant les temps de calcul associés afin d'utiliser la simulation comme un outil dédié à la mise au point. Dans une première partie, le développement d'un modèle 0D orienté simulation système est présenté. La prise en compte de l'ensemble des phénomènes physico-chimiques se déroulant dans la chambre de combustion confère au modèle un niveau de prédictivité conséquent. Un nouveau modèle de combustion de prémélange est proposé, permettant une modélisation détaillée des combustions fortement diluées et des combustions relatives aux injections précoces. Une approche innovante permettant de quantifier les interactions entre les jets pour la multi injection est également proposée. Après calibration sur un nombre restreint d'essais moteur, les résultats du modèle global sont comparés à des mesures expérimentales pour toute la plage de fonctionnement du moteur. La seconde partie de ce travail porte sur la modélisation 1D de la combustion Diesel. Un modèle de jet Diesel est d'abord développé et validé sur des mesures expérimentales. Ce modèle est ensuite étendu à des conditions réactionnelles à l'aide d'un couplage avec un modèle de combustion. Ce dernier s'appuie sur une tabulation des mécanismes de cinétique chimique, ainsi que sur une approche Eddy Break-Up permettant de modéliser le taux de réaction lié au micro mélange. Ce modèle est ensuite intégré à un modèle de chambre de combustion et une première validation du modèle sur des essais moteur réels est entreprise.

Cette thèse porte sur la modélisation 0D/1D de la combustion Diesel dans les moteurs récents. L’objectif est d’augmenter la précision des modèles tout en limitant les temps de calcul associés afin d’utiliser la simulation comme un outil dédié à la mise au point. Dans une première partie, le développement d’un modèle 0D orienté simulation système est présenté. La prise en compte de l’ensemble des phénomènes physico-chimiques se déroulant dans la chambre de combustion confère au modèle un niveau de prédictivité conséquent. Un nouveau modèle de combustion de prémélange est proposé, permettant une modélisation détaillée des combustions fortement diluées et des combustions relatives aux injections précoces. Une approche innovante permettant de quantifier les interactions entre les jets pour la multi injection est également proposée. Après calibration sur un nombre restreint d’essais moteur, les résultats du modèle global sont comparés à des mesures expérimentales pour toute la plage de fonctionnement du moteur. La seconde partie de ce travail porte sur la modélisation 1D de la combustion Diesel. Un modèle de jet Diesel est d’abord développé et validé sur des mesures expérimentales. Ce modèle est ensuite étendu à des conditions réactionnelles à l’aide d’un couplage avec un modèle de combustion. Ce dernier s’appuie sur une tabulation des mécanismes de cinétique chimique, ainsi que sur une approche Eddy Break-Up permettant de modéliser le taux de réaction lié au micro mélange. Ce modèle est ensuite intégré à un modèle de chambre de combustion et une première validation du modèle sur des essais moteur réels est entreprise. / The present thesis focuses on the 0D/1D Diesel combustion modeling of recent engines. The goal is to improve models accuracy while minimizing computation times in order to use simulation as a tool for engine pre-mapping. In the first part, a 0D model designed as a system simulation-oriented tool is proposed. The main contribution of this study is the modeling of the premixed part of the Diesel combustion. This model allows a detailed modeling of highly diluted combustion and combustion related to early injections. A new approach to quantify interactions between each spray in the case of multi injection strategies is also proposed. After calibration using a very small number of engine tests, results for the global combustion chamber model are compared with experimental measurements for the overall engine operating conditions. The second part of this work deals with the 1D Diesel combustion modeling. A Diesel spray model is at first developed and validated on experimental measurements. This model is then extended to reaction conditions using the coupling with a combustion model. The combustion model makes use of tabulated local reaction rates of fuel and is based on the Eddy Break-Up approach to describe the reaction rate related to the turbulent mixing process. The next step is the integration of the burning spray model into a Diesel engine combustion chamber model. A first validation using experimental results for a recent Diesel engine is done.

Ce travail porte sur les aspects numériques de la décomposition en ondelettes et quelques applications aux tourbillons océaniques ainsi qu'a l'estimation de densité de probabilité. On présente une synthèse des différents types de transformée en ondelettes. Nous proposons une étude de régularité et de décroissance a l'infini des ondelettes complexes progressives issues d'ondelettes réelles. On s'intéresse au calcul numérique de la transformée en ondelettes dans les cas unidimensionnels et bidimensionnels. Les aspects de l'implémentation, à l'aide de la fft (fast Fourier transform), de l'algorithme général de la transformée en ondelettes sont présentés. Nous montrons une application de l'analyse par ondelettes a des signaux électriques industriels ainsi qu'a des tourbillons océaniques engendres par un modèle numérique du gulf stream. On introduit une approche originale d'analyse spectrale locale par ondelettes fondée sur l'utilisation de la transformée en ondelettes continue. Nous présentons également une methode de filtrage numérique par ondelettes orthogonales qui permet d'extraire une structure spécifique dans un signal. La dernière partie est relative a l'application des ondelettes a l'estimation de densité de probabilité. Nous développons un algorithme de calcul efficace, fonde sur l'utilisation des ondelettes a support compact de Daubechies, pour construire un estimateur d'une densité inconnue de probabilité

La cristallisation en solution est un procédé largement utilisé dans l'industrie comme opération de séparation et de purification qui a pour but de produire des solides avec des propriétés spécifiques. Les propriétés concernant la taille et la forme ont un impact considérable sur la qualité finale des produits. Il est donc primordial de pouvoir déterminer la distribution granulométrique (DG) des cristaux en formation. En utilisant une caméra in situ, il est possible de visualiser en temps réel les projections 2D des particules 3D présentes dans la suspension. La projection d'un objet 3D sur un plan 2D entraîne nécessairement une perte d'informations : déterminer sa taille et sa forme à partir de ses projections 2D n'est donc pas aisé. C'est tout l'enjeu de ce travail: caractériser géométriquement et morphométriquement des objets 3D à partir de leurs projections 2D. Tout d'abord, une méthode basée sur le maximum de vraisemblance des fonctions de densité de probabilité de mesures géométriques projetées a été développée pour déterminer la taille d'objets 3D convexes. Ensuite, un descripteur de forme stéréologique basé sur les diagrammes de forme a été proposé. Il permet de caractériser la forme d'un objet 3D convexe indépendamment de sa taille et a notamment été utilisé pour déterminer les facteurs d'anisotropie des objets 3D convexes considérés. Enfin, une combinaison des deux études précédentes a permis d'estimer à la fois la taille et la forme des objets 3D convexes. Cette méthode a été validée grâce à des simulations, comparée à une méthode de la littérature et utilisée pour estimer des DGs d'oxalate d'ammonium qui ont été comparées à d'autres méthodes granulométriques.

Un des dangers chimiques les plus importants relevés par l'OMS concerne la contamination des céréales par les mycotoxines et notamment les aflatoxines et les fumonisines. La réglementation recommande des contaminations maximales d'aflatoxines dans les céréales inférieures à 20 mg/kg ; cependant on relève couramment des taux supérieurs à 200 mg/kg dans le maïs au Mexique. Bien qu'il ait été évalué que le processus de nixtamalisation détruit plus de 90 % des fumonisines et de 80 à 95 % des aflatoxines, le taux résiduel peut encore être élevé : certaines publications rapportent des concentrations jusqu'à 100 mg/kg dans les tortillas, ce qui représente un risque avéré vue la grande consommation de tortillas au Mexique (325g/j). Le JECFA (2001) a établi une dose maximale acceptable de 2μg/kg pc/j pour la fumonisine et recommande de réduire l'exposition aux aflatoxines au plus faible niveau possible en considérant que le seuil de 1 ng/kg pc/j ne devrait pas être atteint. Au cours de cette thèse 3 enquêtes aléatoires et représentatives ont été menées dans 40 tortillerias de la ville de Veracruz. La consommation de maïs de la population a été évaluée à partir d'un questionnaire de consommation. L'analyse des mycotoxines a été réalisée par HPLC-FD par utilisation de colonnes à immunoaffinité selon la réglementation européenne (CIRAD-Montpellier). L'analyse des données obtenues a été effectuée selon une méthode probabiliste permettant de construire une fonction de distribution de probabilités à partir de la méthode de Monte Carlo (UBO). La représentativité de la population a été validée à partir d'évaluation de quotas de population après échantillonnage aléatoire initial. La contamination des tortillas a été mesurée à 0.54-1.96 mg/kg pour les aflatoxines et à 65-136 mg/kg pour les fumonisines. La consommation moyenne de tortillas a été mesurée à 148 g de maïs par jour. L'exposition de la population aux aflatoxines apparaît alors comprise entre 0,94 et 3,14 ng/kg pc/j et celle aux fumonisines entre 146 et 315 ng/kg pc/j. Les échantillons les plus contaminés proviennent des tortillerias réalisant elles-mêmes leur procédure de nixtamalisation. L'analyse des résultats montre que 60 % de la population de Veracruz serait à risque selon les préconisations du JECFA. L'exposition aux fumonisines atteint 5 % de la dose maximale acceptable, du fait d'une relativement faible contamination du maïs à cette mycotoxine. Les résultats montrent donc un risque sanitaire pour la population de la ville de Veracruz. Une extension de ce travail à la totalité de l’Etat de Veracruz, incluant la population rurale, devrait être menée du fait du risque probablement accru de cette dernière catégorie de population en lien avec sa plus forte consommation de maïs. / One of the chemical hazards that WHO has reported more frequently is cereals contamination with mycotoxins, mainly aflatoxins and fumonisins. NOM-188-SSA1-2002 establishes that aflatoxin concentration in grain should not exceed 20 mg kg-1 ; however, there are reported concentrations > 200 mg kg-1 in maize. Although it has been documented that nixtamalizacion removes more than 90% of fumonisins and between 80 and 95% of aflatoxins, the residual amount could be important, finding reports concentrations higher than 100 mg kg-1 of aflatoxin in tortilla, representing a risk due to the high consumption of tortillas in Mexico (325 g d-1). The JECFA (2001) establishes a maximum intake of 2 mg kg-1 pc d-1 for fumonisin and aflatoxin recommends reducing “as low as reasonably achievable” levels. 3 random and representative sampling in Veracruz city, each in 40 tortillerias, were made. Corn intake and weight of the population were estimated using a consumption questionnaire. Mycotoxins analysis were performed by HPLC-FD using immunoaffinity columns according to European standard UNE-EN ISO 14123 : 2008 for aflatoxins and UNE-EN 13585 : 2002 for fumonisin in the CIRAD (Montpellier, France). Statistical analysis were performed under a probabilistic approach in collaboration with the University of Bretagne Occidentale (Brest, France), building probability density function (PDF) and using the Monte Carlo method. PDF parameters of the weight of the population was 74.15kg for men (which coincides with reported by CANAIVE) and 65.83kg for women ; the pollution aflatoxin tortilla was 0.54 – 1.96mg kg-1 and fumonisin from 65.46 – 136.00mg kg-1 ; the tortilla consumption was 148.3g of corn per person per day ; the daily intake of aflatoxins was 0.94 – 3.14ng kg-1 bw d-1 and fumonisin of 146.24 – 314.99ng kg-1 bw d-1. Samples with higher aflatoxin contamination came from tortillerias that make the nixtamalization in situ. In assessing exposure it was found that up to 60% of the population could be consuming more than the recommended by JECFA (2001) for aflatoxin dose (1ng kg-1 bw d-1). Exposure to fumonisins intake was < 5% due to low contamination by these mycotoxins. The results suggest that the population of the city of Veracruz could be in food risk by eating contaminated corn tortillas AFT. It is advisable to extend this study in rural communities, where risk factors could increase.

Dans cette thèse, une approche combinant des données issues de capteurs hétérogènes pour la reconnaissance d'activités des personnes âgées à domicile est proposée. Cette approche consiste à combiner les données fournies par des capteurs vidéo avec des données fournies par des capteurs environnementaux pour suivre l'interaction des personnes avec l'environnement. La première contribution est un nouveau modèle de capteur capable de donner une représentation cohérente et efficace des informations fournies par différents types de capteurs physiques. Ce modèle inclue l'incertitude sur la mesure. La deuxième contribution est une approche, basée sur une fusion multicapteurs, pour la reconnaissance d'activités. Cette approche consiste à détecter la personne, suivre ses mouvements, reconnaître ses postures et ses activités d'intérêt, par une analyse multicapteurs et une reconnaissance d'activités humaines. Pour résoudre le problème de la présence de capteurs hétérogènes, nous avons choisi de réaliser la fusion à haut niveau (niveau événement) des différentes données issues des différents capteurs, en combinant les événements vidéo avec les événements environnementaux. La troisième contribution est l'extension d'un langage de description qui permet aux utilisateurs (ex. le corps médical) de décrire les activités d'intérêt dans des modèles formels. Les résultats de cette approche sont montrés pour la reconnaissance des AVQ pour de vraies personnes agées évoluant dans un appartement expérimental appelé GERHOME équipé de capteurs vidéo et de capteurs environnementaux. Les résultats obtenus de la reconnaissance des différentes AVQ sont encourageants.

L’Analyse Linéaire de Stabilité a été largement utilisée afin de décrire l’évolution du mode dominant des instabilités pour différentes configurations. Toutefois, plusieurs observations expérimentales et numériques ont démontré l’existence de distributions de longueurs inter-strictions, prémices aux distributions de tailles de fragments. Une extension de cette utilisation classique de l’Analyse Linéaire de Stabilité dédiée à l’extension d’un barreau en traction a été proposée afin de rendre compte de la contribution de l’ensemble des modes de perturbation sur le profil final des strictions. Cette approche, correspondant physiquement au cas des anneaux fins en expansion, est à même de déterminer une distribution de longueurs inter-strictions. En se basant sur les mêmes outils mathématiques que ceux ayant permis d’étendre l’analyse linéaire de stabilité, une nouvelle approche numérique a été suivie afin de déterminer le temps d’apparition et le nombre de strictions en fonction de la vitesse de chargement, de l’amplitude et de la taille de cellules de perturbations. Une comparaison, en termes de distribution de longueurs inter-strictions, a été par la suite obtenue entre résultats analytiques et numériques. De bonnes corrélations ont été observées. Le développement d’un montage expérimental d’expansion d’anneaux par forces électromagnétiques a été entamé au cours de ce travail de thèse. Lors de sa finalisation, il devrait permettre une validation expérimentale des approches développées dans les domaines numériques et analytiques / Linear Stability Analysis has been widely used in order to describe the evolution of the dominant necking pattern in different configurations. However, experimental and numerical observations have demonstrated that a distribution of internecking lengths is obtained instead of a unique dominant pattern. In this work, an extension of the classical Linear Stability Analysis applied to the dynamic extension of a round bar case has been developed to take into account the contribution of all perturbation modes on the final necking pattern. This approach, corresponding physically to the case of thin ring expansion, is able to determine a distribution on internecking distances. Based on the same mathematical tools that was developed for the extension of the linear stability analysis, a new numerical approach has been followed to determine the number and the onset time of necking with respect to loading velocity, amplitude and cell size of the perturbations. A comparison, in terms of distributions of inter-necking distance, has been carried out between analytical and numerical results. Good correlations have been observed. An electromagnetic ring expansion experimental device has been proposed during this work. Whith some additional development, it should enable to validate experimentally the working assumptions developed in the numerical and analytical studies

La législation sur les émissions de polluants des Moteurs à Combustion Interne (ICEs) est de plus en plus contraignante et représente un gros défi pour les constructeurs automobiles. De nouvelles stratégies de combustion telles que la Combustion à Allumage par Compression Homogène (HCCI) et l’exploitation de stratégies d’injections multiples sont des voies prometteuses qui permettent de respecter les normes sur les émissions de NOx et de suies, du fait que la combustion a lieu dans un mélange très dilué et par conséquent à basse température. Ces aspects demandent la création d’outils numériques adaptés à ces nouveaux défis. Cette thèse présente le développement d’un nouveau modèle 0D de combustion Diesel HCCI : le dual Combustion Model (dual - CM). Le modèle dual-CM a été basé sur l’approche PCM-FPI utilisée en Mécanique des Fluides Numérique (CFD) 3D, qui permet de prédire les caractéristiques de l’auto-allumage et du dégagement de chaleur de tous les modes de combustion Diesel. Afin d’adapter l’approche PCM-FPI à un formalisme 0D, il est fondamental de décrire précisément le mélange à l’intérieur du cylindre. Par consequent, des modèles d’évaporation du carburant liquide, de formation de la zone de mélange et de variance de la fraction de mélange, qui permettent d’avoir une description détaillée des proprietés thermochimiques locales du mélange y compris pour des configurations adoptant des stratégies d’injections multiples, sont proposés. Dans une première phase, les résultats du modèle ont été comparés aux résultats du modèle 3D. Ensuite, le modèle dual-CM a été validé sur une grande base de données expérimentales; compte tenu du bon accord avec l’expérience et du temps de calcul réduit, l’approche présentée s’est montrée prometteuse pour des applications de type simulation système. Pour conclure, les limites des hypothèses utilisées dans dual-CM ont été investiguées et des perspectives pour les dévélopements futurs ont été proposées. / More and more stringent restrictions concerning the pollutant emissions of Internal Combustion Engines (ICEs) constitute a major challenge for the automotive industry. New combustion strategies such as Homogeneous Charge Compression Ignition (HCCI) and the implementation of complex injection strategies are promising solutions for achieving the imposed emission standards as they permit low NOx and soot emissions, via lean and highly diluted combustions, thus assuring low combustion temperatures. This requires the creation of numerical tools adapted to these new challenges. This Ph.D presents the development of a new 0D Diesel HCCI combustion model : the dual Combustion Model (dual−CM ). The dual-CM is based on the PCM-FPI approach used in 3D CFD, which allows to predict the characteristics of Auto-Ignition and Heat Release for all Diesel combustion modes. In order to adapt the PCM-FPI approach to a 0D formalism, a good description of the in-cylinder mixture is fundamental. Consequently, adapted models for liquid fuel evaporation, mixing zone formation and mixture fraction variance, which allow to have a detailed description of the local thermochemical properties of the mixture even in configurations adopting multiple injection strategies, are proposed. The results of the 0D model are compared in an initial step to the 3D CFD results. Then, the dual-CM is validated against a large experimental database; considering the good agreement with the experiments and low CPU costs, the presented approach is shown to be promising for global engine system simulations. Finally, the limits of the hypotheses made in the dual-CM are investigated and perspectives for future developments are proposed.

L'imagerie intravasculaire ultrasonore (IVUS) est une technologie médicale par cathéter qui produit des images de coupe des vaisseaux sanguins. Elle permet de quantifier et d'étudier la morphologie de plaques d'athérosclérose en plus de visualiser la structure des vaisseaux sanguins (lumière, intima, plaque, média et adventice) en trois dimensions. Depuis quelques années, cette méthode d'imagerie est devenue un outil de choix en recherche aussi bien qu'en clinique pour l'étude de la maladie athérosclérotique. L'imagerie IVUS est par contre affectée par des artéfacts associés aux caractéristiques des capteurs ultrasonores, par la présence de cônes d'ombre causés par les calcifications ou des artères collatérales, par des plaques dont le rendu est hétérogène ou par le chatoiement ultrasonore (speckle) sanguin. L'analyse automatisée de séquences IVUS de grande taille représente donc un défi important. Une méthode de segmentation en trois dimensions (3D) basée sur l'algorithme du fast-marching à interfaces multiples est présentée. La segmentation utilise des attributs des régions et contours des images IVUS. En effet, une nouvelle fonction de vitesse de propagation des interfaces combinant les fonctions de densité de probabilité des tons de gris des composants de la paroi vasculaire et le gradient des intensités est proposée. La segmentation est grandement automatisée puisque la lumière du vaisseau est détectée de façon entièrement automatique. Dans une procédure d'initialisation originale, un minimum d'interactions est nécessaire lorsque les contours initiaux de la paroi externe du vaisseau calculés automatiquement sont proposés à l'utilisateur pour acceptation ou correction sur un nombre limité d'images de coupe longitudinale. La segmentation a été validée à l'aide de séquences IVUS in vivo provenant d'artères fémorales provenant de différents sous-groupes d'acquisitions, c'est-à-dire pré-angioplastie par ballon, post-intervention et à un examen de contrôle 1 an suivant l'intervention. Les résultats ont été comparés avec des contours étalons tracés manuellement par différents experts en analyse d'images IVUS. Les contours de la lumière et de la paroi externe du vaisseau détectés selon la méthode du fast-marching sont en accord avec les tracés manuels des experts puisque les mesures d'aire sont similaires et les différences point-à-point entre les contours sont faibles. De plus, la segmentation par fast-marching 3D s'est effectuée en un temps grandement réduit comparativement à l'analyse manuelle. Il s'agit de la première étude rapportée dans la littérature qui évalue la performance de la segmentation sur différents types d'acquisition IVUS. En conclusion, la segmentation par fast-marching combinant les informations des distributions de tons de gris et du gradient des intensités des images est précise et efficace pour l'analyse de séquences IVUS de grandes tailles. Un outil de segmentation robuste pourrait devenir largement répandu pour la tâche ardue et fastidieuse qu'est l'analyse de ce type d'images. / Intravascular ultrasound (IVUS) is a catheter based medical imaging technique that produces cross-sectional images of blood vessels. These images provide quantitative assessment of the vascular wall, information about the nature of atherosclerotic lesions as well as the plaque shape and size. Over the past few years, this medical imaging modality has become a useful tool in research and clinical applications, particularly in atherosclerotic disease studies. However, IVUS imaging is subject to catheter ring-down artifacts, missing vessel parts due to calcification shadowing or side-branches, heterogeneously looking plaques and ultrasonic speckle from blood. The automated analysis of large IVUS data sets thus represents an important challenge. A three-dimensional segmentation algorithm based on the multiple interface fast-marching method is presented. The segmentation is based on region and contour features of the IVUS images: a new speed fonction for the interface propagation that combines the probability density functions (PDFs) of the vessel wall components and the intensity gradients is proposed. The segmentation is highly automated with the detection of the lumen boundary that is fully automatic. Minimal interactions are necessary with a novel initialization procedure since initial contours of the external vessel wall border are also computed automatically on a limited number of longitudinal images and then proposed to the user for acceptance or correction. The segmentation method was validated with in-vivo IVUS data sets acquired from femoral arteries. This database contained 3 subgroups: pullbacks acquired before balloon angioplasty, after the intervention and at a 1 year follow-up examination. Results were compared with validation contours that were manually traced by different experts in IVUS image analysis. The lumen and external wall boundaries detected with the fast-marching method are in agreement with the experts' manually traced contours with similarly found area measurements and small point-to-point contour differences. In addition, the 3D fast-marching segmentation method dramatically reduced the analysis time compared to manual tracing. Such a valdiation study, with comparison between pre- and post-intervention data, has never been reported in the IVUS segmentation literature. In conclusion, the fast-marching method combining the information on the gray level distributions and intensity gradients of the images is precise and efficient to analyze large IVUS sequences. It is hoped that the fast-marching method will become a widely used tool for the fastidious and difficult task of IVUS image processing.