L'imagerie cone-beam computed tomography (CBCT) est une méthodologie de contrôle non destructif permettant l'obtention d'images volumiques d'un objet. Le système d'acquisition se compose d'un tube à rayons X et d'un détecteur plan numérique. La recherche développée dans ce manuscrit se déroule dans le contexte industriel. L'objet est placé sur une platine de rotation et une séquence d'images 2D est acquise. Un algorithme de reconstruction procure des données volumiques de l'atténuation de l'objet. Ces informations permettent de réaliser une étude métrologique et de valider ou non la conformité de la pièce imagée. La qualité de l'image 3D est dégradée par différents artefacts inhérents à la plateforme d'acquisition. L'objectif de cette thèse est de mettre au point une méthode de correction adaptée à une plateforme de micro-tomographie par rayons X d'objets manufacturés poly-matériaux. Le premier chapitre décrit les bases de la physique et de l'algorithmie propres à la technique d'imagerie CBCT par rayons X ainsi que les différents artefacts nuisant à la qualité de l'image finale. Le travail présenté ici se concentre sur deux types d'artefacts en particulier: les rayonnements secondaires issus de l'objet et du détecteur et le durcissement de faisceau. Le second chapitre présente un état de l'art des méthodes visant à corriger le rayonnement secondaire. Afin de quantifier le rayonnement secondaire, un outil de simulation basé sur des techniques de Monte Carlo hybride est développé. Il permet de caractériser le système d'acquisition installé au laboratoire de façon réaliste. Le troisième chapitre détaille la mise en place et la validation de cet outil. Les calculs Monte Carlo étant particulièrement prohibitifs en terme de temps de calcul, des techniques d'optimisation et d'accélération sont décrites. Le comportement du détecteur est étudié avec attention et il s'avère qu'une représentation 2D suffit pour modéliser le rayonnement secondaire. Le modèle de simulation permet une reproduction fidèle des projections acquises avec le système réel. Enfin, le dernier chapitre présente la méthodologie de correction que nous proposons. Une première reconstruction bruitée de l'objet imagé est segmentée afin d'obtenir un modèle voxélisé en densités et en matériaux. L'environnement de simulation fournit alors les projections associées à ce volume. Le volume est corrigé de façon itérative. Des résultats de correction d'images tomographiques expérimentales sont présentés dans le cas d'un objet mono-matériaux et d'un objet poly-matériaux. Notre routine de correction réduit les artefacts de cupping et améliore la description du volume reconstruit. / Cone-beam computed tomography (CBCT) is a standard nondestructive imaging technique related to the acquisition of three-dimensional data. This methodology interests a wide range of applications. An industrial CBCT system comprises an X-ray source and a flat-panel detector. Radiographic images are acquired during a rotation of the object of interest. A reconstruction algorithm leads to a volumic representation of the object and a post-processing routine assesses its validity. Accurate quantitative reconstruction is needed to perform an efficient diagsnotic. However, it is challenged by the presence of different artefacts coming from the acquisition itself. This thesis aims at analyzing and correcting those artefacts in a context of industrial micro-tomography. After an introduction to the physical and algorithmic background of CBCT, the artefacts are presented. Our study adresses two major artefacts: beam hardening and scatter radiations coming from the object and the detector. The second chapter reports on the state of the art in secondary radiation correction. A simulation model of the CBCT imaging chain is developed in a Monte Carlo environment. This model is designed to be realistic in order to get an accurate insight on the processes contributing to the final image formation. The third chapter focuses on the built and validation of the simulation tool. Monte Carlo methods are exact but prohibitively slow. Consequently, acceleration and optimization techniques are used to speed-up the calculations without loss of accuracy. A layer model of the flat-panel detector gives some insight on its secondary radiation behavior. More specifically, we demonstrate that a 2D description of the detector would be sufficient to compute its contribution. Our projection tool fits well with the real system. Finally, the last chapter describes our iterative correction method. The noisy initial reconstruction is segmented into different materials and densities and fed to the simulation framework. Beam hardening and secondary radiations are corrected via the volume reconstructed from the difference between acquired and simulated projections. This correction method is shown to be effective on both mono-material and poly-material objects.

Le travail réalisé durant cette thèse a été motivé par l'étude des mécanismes microscopiques à l'origine du fluage dans les matériaux granulaires.%En particulier, on cherche à explorer des techniques expérimentales et numériques pour l'étude d'un tel phénomène.Dans une première partie, on cherche à mesurer les déplacements des grains dans un matériau granulaire par observations en micro-tomographie X. Une telle identification ne peut être efficacement réalisée pour des phénomènes rapides avec les méthodes classiques de corrélation d'images numériques. Une nouvelle méthode nommée emph{corrélation discrète des projections numériques} qui contourne cette difficulté est développée dans cette thèse. Cette méthode, basée sur la corrélation des projections de tomographie, permet de mesurer les déplacements avec un nombre réduit de projections (100 fois moins que les méthodes classiques), ce qui diminue énormément le temps d'acquisition nécessaire pour la mesure. La méthode, appliquée à des données expérimentales, donne une précision comparable à celles des méthodes classiques tandis que le temps d'acquisition nécessaire est réduit à quelques minutes. Une étude portant sur l'analyse des sources d'erreurs affectant la précision des résultats est également présentée.Le but de la deuxième partie est de réaliser des simulations numériques pour fournir une caractérisation de l'essai oedométrique. Différents assemblages de billes de verre légèrement poly-disperses interagissant à travers des contacts élastiques de Hertz-Mindlin et frottement de Coulomb ont été utilisés. Ces simulations ont permis d'étudier l'évolution de certains paramètres structuraux du matériau modèle, préparant ainsi le terrain pour de futures études sur le fluage. Il a été particulièrement souligné que les contacts élastiques utilisés dans ces simulations ne reproduisent pas l'irréversibilité des déformations observée dans les expériences sur des sables. Cependant, l'irréversibilité est bien visible sur le nombre de coordination et l'anisotropie. Alors que les paramètres élastiques peuvent exprimer la réponse pour des petits incréments de déformations, la compression oedometrique est belle et bien anélastique, principalement à cause de la mobilisation du frottement. Le rapport entre les contraintes horizontales et verticales (coefficient du sol au repos) n'est particulièrement constant que lorsque l'anisotropie de structure est instaurée dans l'état initial de l'assemblage. Il est par ailleurs relié à l'anisotropie interne de la structure par une formule simple. Finalement, les coefficients du tenseur élastique dépendent principalement du nombre de coordination et son anisotropie est plus liée à l'anisotropie des contacts qu'à celle des forces / The present work is motivated by the study of creep in granular materials at the microscopic scale.The first part of this thesis deals with displacement measurements by microtomography. Classical digital image correlation fails to catch time-dependent (possibly fast) phenomena such as short-term creep. A new method named emph{Discrete Digital Projection Correlation} is developed to overcome this limitation. This method requires very few projections (about 100 times less than classical methods) of the deformed state to perform the correlation and retrieve grain displacements. Therefore, the acquisition time is remarkably reduced, which allows to study time-dependent phenomena.The method is tested on experimental data. While its accuracy compares favorably to that of conventional methods, it only requires acquisition times of a few minutes. The origins of measurement errors are tracked by numerical means, on simulated grain displacements and rotations.The second part is a numerical simulation study, by the Discrete Element Method (DEM), of oedometric compression in model granular materials, carried out with a simple model material: assemblies of slightly polydisperse spherical beads interacting by Hertz-Mindlin contact elasticity and Coulomb friction. A wide variety of initial states are subject to compression, differing in density, coordination number and fabric anisotropy. Despite apparently almost reversible strains, oedometric compression proves an essentially anelastic and irreversible process,due to friction, with important internal state changes affecting coordination number and anisotropy. Elastic moduli only describe the response to very small stress increments about well equilibrated configurations. The ratio of horizontal stress to vertical stress (or coefficient of earth pressure at rest, commonly investigated in soil mechanics) only remains constant for initially anisotropic assemblies. A simple formula relates it to force and fabric anisotropy parameters, while elastic moduli are mainly sensitive to the latter. Further studies of contact network instabilities and rearrangements should pave the way to numerical investigations of creep behavior