Brevet FR3126779A1 — publié le 20230310

L’invention concerne un procédé de contrôle non destructif d’une pièce par radiographie en transmission, comportant les étapes suivantes :acquisition de N projections ( ) de la pièce, génération (E2) de N images ( ) calculées de la pièce, estimation (E3a, E3b, E3c) par itérations successives, du vecteur à partir d’un vecteur initial et du vecteur à partir d’un vecteur initial et/ou du vecteur de paramètres à partir d’un vecteur initial , par minimisation de la norme quadratique d’une différence calculée entre les projections ( ) et les images ( ), traitement des projections ( ) et/ou des images ( ), identification de défauts de la pièce par comparaison des projections ( ) traitées et des images traitées ( ). Figure pour l’abrégé : Figure. 4 Procédé, dispositi f et programme d’ordinateur de contrôle d’une pièce par radiographie aux rayons X L'invention concerne un procédé, un dispositif et un programme d’ordinateur pour le contrôle non destructif d’une pièce par radiographie aux rayons X. Le domaine de l’invention concerne les pièces aéronautiques, notamment les aubes de turbomachines, notamment leurs aubes de turbine. Le Contrôle Non Destructif (CND) des pièces aéronautiques est un élément essentiel de la sûreté de fonctionnement des avions dont l’objectif est d’éviter tout défaut qui pourrait provoquer une défaillance en vol. Parmi les méthodes de CND, la radiographie rayons X se distingue par sa capacité à visualiser l’intérieur de la pièce de manière peu/pas intrusive et à résoudre des détails jusqu’à une taille micrométrique. Ces performances l’ont rendu populaire dans le secteur aéronautique. La radiographie rayons X permet d’imager la structure interne de la pièce en transmission. La tomographie consiste en l’acquisition d’un grand nombre de radiographies lors d’une rotation, la plupart du temps complète, d’une pièce dans le but de calculer une image tridimensionnelle complète d’une pièce. L'important temps d'acquisition de ces images tomographiques conduit les industriels à ne considérer qu'un nombre limité d'images radiographiques (appelées radiographies ou projections par la suite) pour réaliser le CND de santé matière et dimensionnel (indications géométriques tridimensionnelles). Toutefois, des artefacts d’images dus au durcissement de faisceau et à la diffusion Compton, affectent les radiographies et rendent délicate voire incertaine l’analyse métrologique effectuée lors du CND. Lorsqu’un grand nombre de projections est considéré, par exemple dans le cadre de la tomographie, ces artefacts peuvent être négligés. En revanche, ils doivent être pris soigneusement en compte lorsque le CND est réalisé à partir de peu de vues. Les cabines de CND de radiographie numérique (des marques déposées telles que Yxlon, GE, Nikon) ou les logiciels d’exploitation (Vg Studio, Aviso, RX solution) et de traitement d’images obtenues par radiographie aux rayons X présentent tous des filtres pour aider le contrôleur dans sa tâche de sanction. Une radiographie (normalisée) est interprétée comme une image de l’atténuation des rayons X lors de la traversée de la pièce, atténuation elle-même reliée à l’épaisseur par une loi qui est souvent approximée par une fonction exponentielle (loi de Beer-Lambert). Les grandeurs d'intérêt dans la caractérisation des indications géométriques sont de nature tridimensionnelle, et leur estimation à partir des radiographies (images bidimensionnelles) est partielle et peu fiable. Ceci peut induire des erreurs non négligeables dans le contrôle, et rend la procédure de sanction sensible à la variation de l’angle de vue lors de l'acquisition de chaque image. L’analyse demande donc une connaissance précise de la géométrie du système cabine-pièce, et de la correction fidèle des artefacts d’images. L’utilisation de la tomographie aux rayons X implique l’acquisition d’un ou plusieurs milliers de projections, ce qui est couteux en temps et entraine également des étapes de traitement à considérer. Ainsi, il est préférable de n’acquérir qu’un faible nombre (de l’ordre d’une dizaine) de projections, soit typiquement une centaine de fois moins. La sanction de pièce par radiographie aux rayons X à partir d'un faible nombre de vue est habituellement effectuée manuellement : des contrôleurs, techniciens spécialistes formés à cette tâche, analysent les images produites par le système à la recherche de toute anomalie. Les anomalies sont représentées par une variation anormale des niveaux de gris des images. Toutefois, des artefacts d'images altèrent ces niveaux de gris, rendant incertaine et difficile la sanction : lorsqu’un faible nombre d’images est considéré, les artefacts d’images ont un fort poids et ne peuvent pas être négligés. Ainsi, la réduction d’artefacts d’imageries pour la sanction de pièces aéronautiques à partir d’un nombre limité de projections est un élément crucial pour cette application. Pour des tensions d’accélération autour de 350 keV (typiquement utilisées pour acquérir des images d’aube de turbine), les artefacts à traiter sont essentiellement le durcissement de faisceau et la diffusion Compton. Une variabilité inter et intra contrôleur existe, réduisant la fiabilité de la sanction. Par ailleurs, l’analyse minutieuse des images par les contrôleurs constitue un travail laborieux, fatiguant, aussi bien sur le plan physique que sur le plan psychologique. Un objectif de l’invention est d’obtenir un procédé, un dispositif et un programme d’ordinateur pour le contrôle non destructif d’une pièce par radiographie aux rayons X, qui pallient les inconvénients mentionnés ci-dessus, en minimisant un ou plusieurs types d’artefacts d’images permettant alors l’exploitation d’un faible nombre de projections acquises pour réaliser le CND. A cet effet, un premier objet de l’invention est un procédé de contrôle non destructif d’une pièce par radiographie en transmission, comportant les étapes suivantes, exécutées par un calculateur : acquisition de N projections de la pièce à partir d’un dispositif de radiographie par transmission selon N angles de vue différents et prédéterminés de la pièce, où N est un entier naturel prescrit, génération de N images calculées de la pièce à partir d’un modèle de référence de la pièce correspondant aux N angles de vue et à partir d’un vecteur de paramètres caractérisant la géométrie de projection de l’acquisition pour les N angles de vue à chacune de plusieurs itérations successives, estimation, par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres rend compte du durcissement de faisceau du rayonnement dans la pièce et le vecteur des paramètres caractérise la diffusion Compton du rayonnement dans la pièce, par minimisation du carré de la norme des différences entre les N projections acquises et les N images calculées, traitement des N projections et/ou des N images calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images calculées à partir du vecteur ayant été estimé, identification de défauts de la pièce par comparaison des N projections traitées aux N images calculées traitées. L’invention permet ainsi le contrôle non destructif de pièces aéronautiques par radiographie aux rayons X, à partir d’un nombre limité de projections acquises. L'invention permet la sanction de pièces aéronautiques en utilisant un nombre limité de projections acquises (pouvant être par exemple inférieur ou égal à 1000) de la pièce inspectée. En effet, après application du traitement algorithmique décrit dans l’invention, les images acquises par le système et celles reproduites par simulation peuvent être quantitativement comparées, et une sanction fondée sur leurs différences est effectuée. Le procédé, le dispositif et le programme de contrôle suivant l’invention aident à la sanction de la pièce inspectée. La procédure étant paramétrée systématiquement, cette sanction est alors fiabilisée. L’invention présente une méthode de CND à partir d'un faible nombre de projections. L’invention permet d’augmenter la fiabilité du procédé et dispositif de contrôle, mais également sa performance en permettant la sanction de la pièce avec un nombre limité de projections. Pour cela, un ou plusieurs types d’artefact sont estimés à partir de leur modélisation. Cette estimation est utilisée soit pour corriger les artefacts dans les images acquises, soit pour les reproduire dans des images simulées à partir du modèle de la pièce idéale (par exemple son modèle de conception assistée par ordinateur ou modèle CAO). Cela est rendu possible par la reproduction des modifications en intensité des pixels, qualifiées d’artefacts, dues aux phénomènes physiques qui sont responsable d’écart par rapport à la simple loi de Beer-Lambert, tel que le durcissement de faisceau ou la diffusion Compton. Cette prise en compte permet une comparaison plus pertinente des niveaux de gris entre acquis et simulé, et donc une sanction plus précise. Une fois ces artefacts considérés, la méthode de CND proposée permet de mettre en évidence toutes les indications quelle que soit leur nature (sante matière, géométrique), ou leur position dans la pièce, indépendamment des conditions d’acquisitions utilisées (en suivant le processus de sanction par CND défini dans un contexte industriel), ou encore de la géométrie de la pièce. Le gain de fiabilité ne s’accompagne d’aucun autre coût, si ce n’est celui, modeste, de modélisation de la projection. L'invention concerne notamment la mise en place d'un système de contrôle constitué d’un dispositif d’acquisition de radiographies aux rayons X d’une pièce composée d’un seul matériau connu, couplé à un algorithme de sanction et de validation. Le traitement algorithmique permet de minimiser un ou plusieurs types d’artefacts d’images attendus pour le matériau et la puissance du faisceau de rayons X utilisé, ou bien de les reproduire, de manière à rendre quantitativement comparables les images acquises et les images simulées. Le durcissement de faisceau peut être corrigé par le calculateur en effectuant un étalonnage de l’atténuation puis une correction des intensités des images. Suivant un mode de réalisation de l’invention, le procédé comporte à chaque itération l’estimation du vecteur des paramètres de la géométrie de projection de l’acquisition pour les N angles de vue par le calculateur, cette estimation comportant : le calcul de résidus de projection à partir du vecteur des valeurs initiales le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales , où est le vecteur colonne des paramètres de la géométrie de projection, est le vecteur colonne des valeurs initiales prescrites des paramètres de la géométrie de projection, est un entier naturel désignant le numéro de l’angle de vue et allant de 1 à N, est mis à jour selon , où est le vecteur colonne des paramètres de la géométrie de projection mis à jour, est le vecteur colonne de variation des paramètres de la géométrie de projection et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Suivant un mode de réalisation de l’invention, le procédé comporte à chaque itération l’estimation du vecteur des paramètres d’étalonnage de durcissement de faisceau du rayonnement dans la pièce par le calculateur, cette estimation comportant: le calcul de résidus de projection à partir du vecteur des valeurs initiales , le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales , où est l’image obtenue en appliquant la fonction à l’intensité de chacun des pixels de , sont les N images calculées pour les N angles de vue, sont les N projections acquises de la pièce, est le vecteur colonne des paramètres d’étalonnage de durcissement de faisceau, est un vecteur colonne des valeurs initiales prescrites des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce, est une base de fonctions de forme prescrites, est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , est mis à jour selon , où est le vecteur colonne des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce, ayant été mis à jour, est le vecteur colonne de variation des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Suivant un mode de réalisation de l’invention, le procédé comporte à chaque itération l’estimation du vecteur des paramètres de diffusion Compton du rayonnement dans la pièce par le calculateur, cette estimation comportant: le calcul de résidus de projection à partir d’un vecteur des valeurs initiales , le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales , où sont les N projections acquises de la pièce, sont les N images calculées pour les N angles de vue, est l’image obtenue en convoluant des images simulées, ayant été obtenues à partir au moins des N images calculées, avec le noyau , est le vecteur colonne des paramètres de diffusion Compton du rayonnement dans la pièce, est un vecteur colonne des valeurs initiales prescrites des paramètres de diffusion Compton du rayonnement dans la pièce, est le noyau de convolution défini par est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , un noyau gaussien bidimensionnel d’écart-type prescrit, est la fonction de Dirac au pixel , est mis à jour selon , où est le vecteur colonne des paramètres de diffusion Compton du rayonnement dans la pièce, ayant été mis à jour, est le vecteur colonne de variation des paramètres de diffusion Compton du rayonnement dans la pièce et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Suivant un mode de réalisation de l’invention, est l’image simulée, obtenue en appliquant la fonction à l’intensité de chacun des pixels de l’image calculée . Suivant un mode de réalisation de l’invention, N est inférieur ou égal à 1000. Un deuxième objet de l’invention est un programme d’ordinateur, comportant des instructions de code pour la mise en œuvre des étapes suivantes d’un procédé de contrôle non destructif d’une pièce par radiographie par transmission, lorsqu’il est exécuté par un calculateur: réception de N projections de la pièce à partir d’un dispositif de radiographie par transmission selon N angles de vue différents et prédéterminés de la pièce, où N est un entier naturel prescrit, génération de N images calculées de la pièce à partir d’un modèle de référence de la pièce correspondant aux N angles de vue et à partir d’un vecteur de paramètres caractérisant la géométrie de projection de l’acquisition pour les N angles de vue à chacune de plusieurs itérations successives, estimation, par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres rend compte du durcissement de faisceau du rayonnement dans la pièce et le vecteur des paramètres caractérise la diffusion Compton du rayonnement dans la pièce, par minimisation du carré de la norme des différences entre les N projections acquises et les N images calculées, traitement des N projections et/ou des N images calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images calculées à partir du vecteur ayant été estimé, identification de défauts de la pièce par comparaison des N projections traitées aux N images calculées traitées. Un troisième objet de l’invention est un dispositif de contrôle non destructif d’une pièce par radiographie en transmission, comportant : un dispositif de radiographie en transmission, pour l’acquisition de N projections de la pièce selon N angles de vue différents et prédéterminés de la pièce, où N est un entier naturel prescrit, un calculateur configuré pour effectuer les étapes suivantes : génération de N images calculées de la pièce à partir d’un modèle de référence de la pièce correspondant aux N angles de vue et à partir d’un vecteur de paramètres caractérisant la géométrie de projection de l’acquisition pour les N angles de vue à chacune de plusieurs itérations successives, estimation, par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres rend compte du durcissement de faisceau du rayonnement dans la pièce et le vecteur des paramètres caractérise la diffusion Compton du rayonnement dans la pièce, par minimisation du carré de la norme des différences entre les N projections acquises et les N images calculées, traitement des N projections et/ou des N images calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images calculées à partir du vecteur ayant été estimé, identification de défauts de la pièce par comparaison des N projections traitées aux N images calculées traitées. L’invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d’exemple non limitatif en référence aux figures ci-dessous des dessins annexés. représente schématiquement un synoptique modulaire d’un dispositif de contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente schématiquement un synoptique modulaire d’un dispositif de contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente schématiquement un synoptique modulaire d’un dispositif de contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente schématiquement un synoptique modulaire d’un dispositif de contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente schématiquement un organigramme d’un procédé de contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente schématiquement une courbe de l’intensité des pixels d’une projection dans le cas idéal et dans un cas pratique altéré par des artéfacts de durcissement de faisceau d’une projection obtenue par un dispositif de l’état de la technique. représente un schéma explicatif de la diffusion Compton lors de l’acquisition d’une projection aux rayons X. représente la différence entre, d’une part l’inverse de l’exponentielle de l’intensité d’une projection acquise, et d’autre part l’inverse de l’exponentielle de l’intensité d’une projection simulée, avant correction sur une partie d’une pale d’aube de turbine de turbomachine, obtenue par un dispositif de l’état de la technique. représente la différence entre, d’une part l’inverse de l’exponentielle de l’intensité d’une projection acquise, et d’autre part l’inverse de l’exponentielle de l’intensité d’une projection simulée, avant correction sur une partie d’un pied d’aube de turbine de turbomachine, obtenue par un dispositif de l’état de la technique. représente la différence entre, d’une part l’inverse de l’exponentielle de l’intensité d’une projection acquise, et d’autre part l’inverse de l’exponentielle de l’intensité d’une projection simulée, après correction sur une partie de la même pale d’aube de turbine de turbomachine qu’à la , obtenue par le procédé, un dispositif et un programme d’ordinateur pour le contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. représente la différence entre, d’une part l’inverse de l’exponentielle de l’intensité d’une projection acquise, et d’autre part l’inverse de l’exponentielle de l’intensité d’une projection simulée, après correction sur une partie du même pied d’aube de turbine de turbomachine qu’à la , obtenue par le procédé, un dispositif et un programme d’ordinateur pour le contrôle non destructif d’une pièce par radiographie aux rayons X suivant un mode de réalisation de l’invention. Ainsi qu’illustré aux figures 1 à 5, le dispositif 1 de contrôle non destructif d’une pièce réelle 200 comporte un dispositif 100 de radiographie en transmission permettant d’acquérir, au cours de l’étape E1 du procédé, des projections , également appelées projections acquises brutes, de la pièce 200. Dans ce qui suit, le rayonnement est des rayons X. Le dispositif 100 de radiographie permet d’inspecter la pièce 200 dans des conditions d’acquisition maîtrisée avec une géométrie d’acquisition connue et permet d’imager en transmission la structure interne de la pièce 200. Dans ce qui suit, les projections et les images mentionnées ci-dessous sont définies par leurs niveaux de gris (ou intensité) pour leurs différents pixels . Dans les modes de réalisation décrits ci-dessous, les projections sont des images radiographiques en transmission de rayons X, dans le procédé de contrôle non destructif d’une pièce par radiographie aux rayons X, le dispositif 1 de contrôle non destructif d’une pièce par radiographie aux rayons X et le programme d’ordinateur de contrôle non destructif d’une pièce par radiographie aux rayons X. La pièce 200 est une pièce mécanique, et peut être d’une manière non limitative une pièce aéronautique, notamment une aube de turbomachine (par exemple de turboréacteur), par exemple une de leurs aubes de turbine (pale de turbine), ou peut-être autre. La pièce 200 est en un matériau semi-transparent au rayonnement, ici aux rayons X, et peut être constituée en un seul matériau (mono-matériau). La pièce 200 peut être une pièce métallique ou une pièce composite à matrice céramique, ou autre. Le procédé, le dispositif 1 et le programme d’ordinateur pour le contrôle non destructif de la pièce 200 par radiographie aux rayons X sont employés pour contrôler la pièce 200 au cours de sa fabrication ou lors d’opérations de maintenance, afin de détecter des défauts de cette pièce, pouvant par exemple provoquer une défaillance en vol dans le cas d’une pièce aéronautique. Le dispositif 1 de contrôle non destructif de la pièce 200 comporte et le procédé de dispositif de contrôle non destructif de la pièce 200 utilise un (ou plusieurs) calculateur CAL. Le calculateur CAL peut être ou comprendre un ou plusieurs ordinateurs, un ou plusieurs serveurs, une ou plusieurs machines, un ou plusieurs processeurs, un ou plusieurs microprocesseurs, une ou plusieurs mémoires permanentes MEM, une ou plusieurs mémoires vives MEM. Le calculateur CAL peut comprendre une ou plusieurs interfaces physiques INT1 d’entrée de données, une ou plusieurs interfaces physiques INT2 de sortie de données. Cette ou ces interfaces physiques INT1 d’entrée de données peuvent être ou comprendre un ou plusieurs claviers d’ordinateur, un ou plusieurs ports physiques de communication de données, un ou plusieurs écrans tactiles, ou autres. Cette ou ces interfaces physiques INT2 de sortie de données peuvent être ou comprendre un ou plusieurs ports physiques de communication de données, un ou plusieurs écrans, ou autres. Un programme d’ordinateur peut être enregistré et exécuté sur le calculateur CAL du dispositif 1 de contrôle non destructif de la pièce 200 et comporter des instructions de code, qui lorsqu’elles sont exécutées sur celui-ci, mettent en œuvre tout ou partie du procédé de contrôle non destructif de la pièce 200 suivant l’invention (dont la réception des N projections au cours de l’étape E1). Le dispositif 100 de radiographie aux rayons X comporte une source 101 de rayons X, un support 102 sur lequel se trouve la pièce 200, un mécanisme 104 de contrôle pour faire tourner le support 102 et la source 101 l’un par rapport à l’autre 101 autour d’un axe 103 de rotation, pouvant être par exemple vertical (par exemple la source 101 est fixe et le support 102 est mis en rotation autour de l’axe 103), un détecteur 105 des rayons X traversant la pièce 200, la pièce 200 se trouvant donc sur le trajet des rayons X entre la source 101 et le détecteur 105. La source 101, le support 102 et le détecteur 105 sont disposés dans une cabine à rayons X de haute puissance. Le détecteur 105 fournit les projections de la pièce 200 au cours de la première étape E1. Le mécanisme 104 de contrôle est contrôlé pour l’acquisition aux N angles ANG (n) de vue, différents les uns des autres, de la pièce 200 vis-à-vis des rayons X, par le détecteur 105, de N projections . N est un entier naturel prescrit, supérieur ou égal à 1. L’entier naturel va de 1 à N et désigne le numéro de l’angle ANG (n) de vue respectif et donc le numéro de la projection acquise. Le dispositif 100 de radiographie permet ainsi d’acquérir N projections (c’est-à-dire N projections ) du volume de la pièce 200 sous respectivement les N angles ANG ( n ) de vue. Les N angles ANG ( n ) de vue ont été prédétermines pour être utilisés pour l’étape ultérieure E5 d’analyse de la pièce. L’important temps d'acquisition des projections acquises par rayons X conduit à ne considérer qu'un nombre limité N de projections acquises, inférieur à 1000 voire à 100, ou autre. Le dispositif 100 de radiographie aux rayons X fournit ainsi l’ensemble des projections acquises au cours de l’étape E1. Au cours d’une deuxième étape E2, postérieure à la première étape E1, le calculateur CAL génère N images calculées de la pièce 200 à partir d’un modèle MODP numérique de référence de la pièce 200, ainsi qu’illustré aux figures 1 à 5. Ces N images calculées correspondent aux N angles ANG (n) de vue, c’est-à-dire adoptent le même angle ANG (n) de vue que respectivement les N projections . Ce modèle MODP est une référence géométrique de la pièce 200, est enregistré à l’avance dans une mémoire du calculateur CAL et peut être par exemple un modèle de conception assistée par ordinateur (ou CAO en abréviation) de la pièce 200, reproduisant une pièce 200 idéale sans défauts, ayant une forme extérieure de coordonnées tridimensionnelles définies. Ce modèle MODP peut tenir compte de la composition du matériau de la pièce 200. A chaque étape E2, le calculateur CAL génère N images calculées de la pièce 200 à partir du modèle MODP de référence de la pièce 200 correspondant aux N angles ANG (n) de vue et à partir d’un vecteur de paramètres caractérisant la géométrie de projection de l’acquisition pour les N angles ANG (n) de vue à chacune de plusieurs itérations successives. Le calculateur CAL effectue une ou plusieurs ou toutes les étapes E3a, E3b, E3c d’estimation d’artefacts de projections acquises multi-vues, qui seront décrites ci-dessous. Suivant des modes de réalisation, la détermination des paramètres décrits ci-dessous se fait par une procédure d’optimisation qui exploite des champs de sensibilité pour minimiser des résidus de projection. Au cours de la troisième étape E3a, postérieure à la première étape E1, le calculateur CAL estime le vecteur des paramètres de la géométrie de projection de l’acquisition de l’étape E1 pour les N angles ANG ( n ) de vue, par minimisation du carré de la norme quadratique des différences, calculées par le calculateur CAL, entre les N projections acquises et les N images calculées . Au cours de la troisième étape E3a le calculateur CAL estime le vecteur décrivant la géométrie du dispositif 100 de radiographie par rayons X par rapport au modèle MODP de la pièce 200 idéale (par exemple son modèle de CAO). La troisième étape E3a permet de déterminer précisément par le calculateur CAL le vecteur de la géométrie de projection utilisée lors de l’acquisition des projections par le dispositif 100 de radiographie par rayons X. Cela permet d’assurer une simulation plus fiable des images. La troisième étape E3a permet d’effectuer lors de l’itération suivante de l’étape E2 un recalage géométrique des N images calculées pour les N angles ANG ( n ) de vue en calculant les N images calculées par projection du modèle MODP sur le plan du détecteur 105 selon la géométrie déterminée par le vecteur des N projections acquises. L’étape E3a peut être mise en œuvre par un simulateur SIMX, ainsi que représenté à la figure 4. Au cours de la troisième étape E3a, le calculateur CAL estime le vecteur des paramètres de la géométrie de projection de l’acquisition par les itérations successives à partir d’un vecteur initial . Le vecteur des paramètres de la géométrie de projection est utilisé pour simuler les images calculées . Ainsi, lors de l’étape E2, on utilise déjà le vecteur des paramètres des valeurs initiales prescrites des paramètres de la géométrie de projection lors de la première itération ; puis on utilise lors des itérations suivantes de l’étape E2 les valeurs estimées du vecteur . Les images calculées sont par la suite utilisées pour estimer le vecteur et/ou dans l’étape E3a et l’étape E3b et/ou E3c. Les images calculées ne sont pas corrigées uniquement à l’étape E4, mais également à l’étape E2. Ainsi, l’étape E2 inclut la régénération de l’image calculée à chaque nouvelle itération avec le nouveau vecteur , qui a été estimé lors de l’étape E3a. L’estimation du vecteur lors de l’étape E3a est une procédure itérative, et donc à chaque itération on utilise la meilleure estimation possible (jusqu’à convergence) pour générer les images calculées . Le vecteur est utilisé à l’étape E2 et à chacune des itérations de l’étape E3a. La reproduction ou correction des artéfacts repose sur la simulation numérique des projections acquises par un système rayons X. Les artéfacts sont ensuite estimés en se basant sur une modélisation paramétrique (étapes E3b, E3c) avant d’être soit reproduits (étape E4), soit corrigés (étape E4). Les paramètres de projection (tels les angles ANG ( n ) de projection, distance de la source 101 au détecteur 105) permettent de lier le modèle MODP idéal de la pièce 200 aux projections acquises de la pièce 200 et sont ainsi des paramètres de mise en correspondance (en anglais : registration parameters) du modèle MODP idéal de la pièce 200 avec les projections acquises . Suivant un mode de réalisation de l’invention, au cours de l’étape E3a, le calculateur CAL calcule (par exemple par le calculateur CRES de résidus à la figure 4) à chaque itération les résidus de projection, ces résidus étant les différences entre les projections observées et les projections simulées selon l’équation . A partir de la première approximation prescrite, notée , de ces paramètres , le calculateur CAL estime (par exemple par l’estimateur ESTV à la figure 4) le vecteur décrivant au mieux le système observé, c’est-à-dire minimisant les résidus de projection. Les champs de sensibilité sont calculés (par exemple par différences finies, ou autre) par le calculateur CAL selon l’équation suivante : et traduisent l’effet de la variation de chacun des paramètres sur l’image calculée. Les paramètres de la géométrie de projection sont calculés par le calculateur CAL par optimisation à partir des résidus de projection , de la valeur , et des champs de sensibilité . Le calculateur CAL calcule la variation optimale du vecteur de paramètres et met à jour le vecteur selon l’équation suivante : où est le vecteur colonne des paramètres de la géométrie de projection mis à jour, est le vecteur colonne de variation des paramètres de la géométrie de projection et est calculé par le calculateur CAL comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Le nombre de lignes de la matrice est le nombre de pixels exploités dans les projections , dans les images calculées et dans le résidu . Le nombre de colonnes de la matrice est le nombre de paramètres . Le calculateur CAL calcule lors de l’étape E3a l’ensemble des champs de sensibilité pour allant de 1 à N et allant de 1 à . L’estimation du vecteur des paramètres de la géométrie de projection pour les N angles ANG ( n ) de vue par le calculateur CAL est effectuée par itérations successives sur qui devient , pour : calculer à chaque itération actuelle, à partir du vecteur de l’itération précédente prenant la place du vecteur mentionné ci-dessus les champs de sensibilité selon l’équation mentionnée ci-dessus, pour chaque itération actuelle : calculer un nouveau et mettre à jour selon , avec calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, selon jusqu’à ce que converge vers une meillleure estimation (déterminée par exemple par le fait que la norme de la différence calculée entre le de l’itération actuelle et le de l’itération précédente devient inférieure à un seuil prescrit). Au cours de la quatrième étape E3b, postérieure à la première étape E1, le calculateur CAL estime le vecteur des paramètres d’étalonnage de durcissement de faisceau de rayons X dans la pièce 200 par minimisation du carré de la norme quadratique des différences entre les N projections acquises et les N images calculées, calculé par le calculateur CAL. Le durcissement de faisceau se réfère au fait que l’absorption d’un photon par le matériau constituant la pièce analysée dépend de son énergie. Les photons de faible énergie sont absorbés préférentiellement comparés à ceux de haute énergie. De l’écart entre la courbe réelle CRI donnée par le niveau de gris des pixels d’une projection mesurée et la courbe théorique CTI résultant de la modélisation utilisée dans les analyses des projections (loi de Beer-Lambert), illustrées à la figure 6, émerge des artéfacts dits de durcissement de faisceau. Le durcissement de faisceau induit une relation non linéaire entre l’intensité mesurée d’une projection et l’épaisseur de matière traversée (écart par rapport à la loi de Beer-Lambert) dans le cas d’une pièce mono-matériau. L’épaisseur est l’épaisseur traversée par le rayon avant d’atteindre le pixel . Dans le cas de pièces 200 mono-matériaux, les artéfacts de durcissement de faisceau peuvent être exactement corrigés par étalonnage de l’atténuation de durcissement de faisceau. Le calculateur CAL effectue une identification de la fonction d’étalonnage reliant l’atténuation effective des rayons X traversant la pièce 200 et la longueur de matériau traversée dans la pièce 200. L’étape E3b peut être mise en œuvre par un simulateur SIMBH de durcissement de faisceau de rayons X, ainsi que représenté à la figure 4. Lors de l’étape E3b, le calculateur CAL estime, par les itérations successives, le vecteur à partir d’un vecteur initial . On décrit ci-dessous des modes réalisation de la quatrième étape E3b. Ce qui est décrit ci-dessous est effectué à chacune des itérations de l’étape E3b. Le calculateur CAL étalonne le durcissement de faisceau en identifiant la fonction faisant passer de à , à savoir . Suivant un mode de réalisation de l’invention, l’étalonnage du durcissement de faisceau est effectuée par l’application d’une fonction linéaire par morceau reliant et . Suivant un mode de réalisation de l’invention, la fonction est discrétisée sur une base de fonctions de forme prescrites (par exemple par des éléments finis monodimensionnels, ou des fonctions spline, ou des fonctions polynômes, ou autres) : où est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , est une variable muette désignant le niveau de gris (intensité) d’un pixel d’une image de référence. Le vecteur colonne contenant les paramètres quantifie la fonction d’étalonnage de durcissement de faisceau de rayons X pour allant de 1 à . A partir d’une estimation initiale prescrite de ces paramètres, le calculateur CAL calcule le vecteur décrivant au mieux l’effet observé. Suivant un mode de réalisation de l’invention, le calculateur CAL calcule (par exemple par différences finies) les champs de sensibilité, associés à chaque degré de liberté : à partir du vecteur des valeurs initiales prescrites des paramètres d’étalonnage de durcissement de faisceau de rayons X dans la pièce 200, avec l’image simulée, obtenue en appliquant la fonction (avec les paramètres contenus dans le vecteur ) à chacun des pixels de . Le vecteur des valeurs initiales est tel que . On a ainsi : Les champs de sensibilité traduisent l’effet de la variation de chacun des paramètres sur l’image calculée. Suivant un mode de réalisation de l’invention, au cours de l’étape E3b, le calculateur CAL calcule (par exemple par le calculateur CRES de résidus à la figure 4) les résidus de projection , ces résidus étant les différences entre les projections observées et les projections simulées . A partir d’une première approximation prescrite du vecteur colonne de paramètres , le calculateur CAL estime (par exemple par l’estimateur ESTV à la figure 4) le vecteur décrivant au mieux le système observé, c’est-à-dire minimisant les résidus de projection. Le vecteur est calculé par le calculateur CAL par optimisation à partir des résidus de projection, de la valeur , et des champs de sensibilité. Le calculateur CAL calcule la variation optimale du vecteur de paramètres et met à jour le vecteur selon l’équation suivante : où est le vecteur colonne des paramètres d’étalonnage de durcissement de faisceau de rayons X dans la pièce 200 mis à jour, est le vecteur colonne de variation des paramètres d’étalonnage de durcissement de faisceau de rayons X dans la pièce 200 et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Le nombre de lignes de la matrice est le nombre de pixels exploités dans les projections , dans les images calculées et dans le résidu . Le nombre de colonnes de la matrice est le nombre de paramètres . Le calculateur CAL calcule lors de l’étape E3b l’ensemble des champs de sensibilité pour allant de 1 à N et allant de 1 à . L’estimation du vecteur des paramètres d’étalonnage de durcissement de faisceau par le calculateur CAL est effectuée par itérations successives sur qui devient , pour : calculer à chaque itération actuelle, à partir du vecteur de l’itération précédente prenant la place du vecteur mentionné ci-dessus les champs de sensibilité selon l’équation mentionnée ci-dessus, pour chaque itération actuelle : calculer un nouveau et mettre à jour selon , avec calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, selon jusqu’à ce que converge vers une meillleure estimation (déterminée par exemple par le fait que la norme de la différence calculée entre le de l’itération actuelle et le de l’itération précédente devient inférieure à un seuil prescrit). Au cours de la cinquième étape E3c, postérieure à la première étape E1, le calculateur CAL estime le vecteur des paramètres de diffusion Compton des rayons X dans la pièce 200 par minimisation du carré de la norme quadratique des différences, calculé par le calculateur CAL, entre les N projections acquises et les N images simulées, ayant été calculés à partir des N images calculées. Les paramètres de diffusion Compton représente l’effet de la diffusion Compton sur les projections acquises. Ainsi qu’illustré à la figure 7, la diffusion Compton consiste en l’interaction d’un photon du rayon X avec un électron libre (ou faiblement rattaché à son atome) de la pièce 200 inspectée, provoquant la déviation du photon de sa trajectoire initiale. En atteignant le détecteur 105, le photon ainsi diffusé et dévié PDEV induit une augmentation de l'intensité mesurée (appelée signal secondaire dans ce qui suit) dans la zone DEV située autour du point d’impact direct et central CENT du rayon primaire PCENT, et une diminution de cette intensité (appelée signal primaire dans ce qui suit) en ce point d’impact CENT, créant dans les projections un flou dégradant l’analyse métrologique. Au contraire des méthodes connues utilisant une correction directe de l’effet Compton par une opération de déconvolution qui est numériquement instable, rendant peu fiables ces radiographies retraitées, l’étape E3c prévoit de reproduire cet effet de diffusion Compton sur des images simulées, pour permettre leur comparaison avec les projections acquises. L’étape E3c peut être mise en œuvre par un estimateur ESDC de diffusion Compton, ainsi que représenté à la figure 4. Lors de l’étape E3c, le calculateur CAL estime, par les itérations successives, le vecteur à partir d’un vecteur initial . Ce qui est décrit ci-dessous est effectué pour chaque itération de l’étape E3c. Suivant un mode de réalisation de l’invention, chaque projection acquise est modélisée par le calculateur CAL comme étant la somme de deux composantes : le signal primaire provenant des photons traversant l'objet sans diffusion Compton mais subissant l'absorption par les atomes constituant le matériau (effet photoélectrique); et le signal secondaire qui correspond à la contribution des photons diffusés par diffusion Compton à partir de la pièce 200, selon l’équation suivante : Le signal secondaire, ou diffusé, est calculé par le calculateur CAL comme étant le résultat de la convolution du signal primaire avec un noyau paramétrique : Le calculateur CAL calcule le noyau paramétrique comme étant une somme pondérée de gaussiennes bidimensionnelles de largeurs différentes l’une de l’autre (écart-types différents l’un de l’autre) : où est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , un noyau gaussien bidimensionnel d’écart-type prescrit (suivant une suite ), est la fonction de Dirac au pixel . En variante, les paramètres du noyau gaussien peuvent être identifiés expérimentalement ou via des simulations de Monte-Carlo. Suivant un mode de réalisation de l’invention, en identifiant le signal primaire aux images simulées numériquement , le calculateur CAL calcule le vecteur qui approxime au mieux la diffusion Compton observée. est le vecteur colonne des paramètres de diffusion Compton des rayons X dans la pièce 200, pour allant de 1 à . A partir d’un vecteur des valeurs initiales prescrites des paramètres de diffusion Compton des rayons X dans la pièce 200, le calculateur CAL calcule des images à partir des images simulées, obtenues à partir des N images calculées , selon l’équation : où le signal diffusé est simulé et ajouté aux images de sorte à générer les images . Le calculateur CAL calcule ainsi l’ensemble des N images au cours de l’étape E3c. Suivant un mode de réalisation de l’invention, est l’image obtenue en appliquant la fonction à l’intensité (niveau de gris) de chacun des pixels de l’image calculée . Bien entendu, les images simulées peuvent être obtenues d’une autre manière à partir des N images calculées , dans le cas où l’étape E3c est mise en œuvre sans l’étape E3b. Le vecteur des valeurs initiales prescrites des paramètres de diffusion Compton valeurs initiales est tel que pour . Suivant un mode de réalisation de l’invention, au cours de l’étape E3c, le calculateur CAL calcule (par exemple par le calculateur CRES de résidus à la figure 4) les résidus de projection, ces résidus de projection étant les différences entre les projections acquises et les projections simulées selon l’équation : A partir de la première approximation prescrite du vecteur de paramètres, notée , le calculateur CAL estime (par exemple par l’estimateur ESTV à la figure 4) le vecteur décrivant au mieux le système observé, c’est-à-dire minimisant les résidus de projection. Suivant un mode de réalisation de l’invention, le calculateur CAL calcule (par exemple par différences finies) les champs de sensibilité, associés à chaque degré de liberté : à partir du vecteur colonne des valeurs initiales prescrites des paramètres de diffusion Compton des rayons X dans la pièce 200. Les paramètres de diffusion Compton des rayons X dans la pièce 200 sont calculés par le calculateur CAL par optimisation à partir des résidus de projection, de la valeur , et des champs de sensibilité. Le calculateur CAL calcule la variation optimale du vecteur de paramètres et met à jour le vecteur selon l’équation suivante : où est le vecteur colonne de variation au moins des paramètres de diffusion Compton des rayons X dans la pièce 200 et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice au moins des champs de sensibilité. Le nombre de lignes de la matrice est le nombre de pixels exploités dans les projections , dans les images calculées et dans le résidu . Le nombre de colonnes de la matrice est le nombre de paramètres . Le calculateur CAL calcule lors de l’étape E3c l’ensemble des champs de sensibilité pour allant de 1 à N et j allant de 1 à . L’estimation du vecteur des paramètres de diffusion Compton des rayons X par le calculateur CAL est effectuée par itérations successives sur qui devient , pour : calculer à chaque itération actuelle, à partir du vecteur de l’itération précédente prenant la place du vecteur mentionné ci-dessus les champs de sensibilité selon l’équation mentionnée ci-dessus, pour chaque itération actuelle : calculer un nouveau et mettre à jour selon , avec calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, selon jusqu’à ce que converge vers une meillleure estimation (déterminée par exemple par le fait que la norme de la différence calculée entre le de l’itération actuelle et le de l’itération précédente devient inférieure à un seuil prescrit). Suivant un mode de réalisation de l’invention, N est inférieur ou égal à 1000 ou à 100 ou à 50 ou à 10. Cela traduit le faible nombre de projections acquises malgré lequel l’invention est capable de fonctionner. Ce nombre est inférieur aux nombres de projections acquises nécessaires dans l’état de la technique pour la reconstruction tomographique de la pièce 200. Au cours de la septième étape E4, postérieure aux itérations de l’étape E3a et de l’étape E3b et/ou E3c, le calculateur CAL traite les N projections obtenues à l’étape E1 et/ou les N images calculées, ayant été obtenues après les itérations de l’étape E2. Le calculateur CAL effectue un premier traitement et/ou un deuxième traitement. Le premier traitement comporte une correction du durcissement de faisceau sur les N projections à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images calculées à partir du vecteur ayant été estimé. Le deuxième traitement comporte une correction de la diffusion Compton sur les N projections à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images calculées à partir du vecteur ayant été estimé. Le calculateur obtient ainsi N projections acquises , qui sont égales à ou corrigées à partir des N projections acquises de la pièce 200, et N images simulées , qui sont égales à ou générées à partir des N images calculées de la pièce 200, ces projections acquises étant directement comparables aux N images simulées , pour pouvoir ensuite effectuer l’analyse de la pièce 200 lors de l’étape E5. Par exemple, le calculateur CAL génère d’abord lors de l’étape E2 les N images calculées en leur appliquant le recalage géométrique pour les N angles ANG (n) de vue, et ce en utilisant les paramètres de projection pour simuler les N images calculées. Ce sont ces images calculées ainsi modifiées (recalées) qui sont utilisées comme images calculées dans les traitements ci-dessous. Cette modification par recalage utilisant les paramètres de projection est effectuée au préalable sur les N images calculées dans les premier et deuxième modes de réalisation de l’étape E4, décrits ci-dessous. Suivant un premier mode de réalisation de l’étape E4, le calculateur CAL corrige les projections acquises en leur appliquant l’inverse de la fonction dépendant des paramètres d’étalonnage de durcissement de faisceau de rayons X à chacun des pixels de (pour effectuer ainsi une correction du durcissement de faisceau), pour obtenir les N projections acquises corrigées selon l’équation suivante : Le calculateur CAL calcule la fonction dépendant des paramètres d’étalonnage de durcissement de faisceau de rayons X selon l’équation suivante : où désigne les paramètres d’étalonnage de durcissement de faisceau de rayons X, contenus dans le vecteur et ayant été calculés. Le calculateur CAL traite les N images calculées en les convoluant avec le noyau ayant été calculé en fonction des paramètres de diffusion Compton (reproduction de la diffusion Compton sur les N images calculées), pour obtenir les N images simulées selon les équations suivantes : où désigne les paramètres de diffusion Compton, contenus dans le vecteur et ayant été calculés. Suivant un deuxième mode de réalisation de l’étape E4, le calculateur CAL calcule les projections acquises selon l’équation suivante : Le calculateur traite les N images calculées en leur appliquant la fonction (calculée selon le calcul décrit ci-dessus) dépendant des paramètres d’étalonnage de durcissement de faisceau de rayons X (pour reproduire ainsi le durcissement de faisceau), puis en les convoluant avec le noyau ayant été calculé (calculé selon le calcul décrit ci-dessus) en fonction des paramètres de diffusion Compton (reproduction de la diffusion Compton sur les N images calculées), pour obtenir les N images simulées selon les équations suivantes : Au cours de la huitième étape E5 d’analyse de la pièce 200 réelle, postérieure à l’étape E4, le calculateur CAL identifie des défauts de la pièce 200 par comparaison des N projections traitées aux N images calculées traitées. Le calculateur CAL identifie des défauts de la pièce 200 par comparaison des N projections acquises aux N images simulées suivant l’un des modes de réalisation décrits ci-dessus. Le calculateur CAL effectue une identification de défauts de la pièce par comparaison des images acquises et simulées ayant pris en compte la correction ou le traitement soit sur les unes soit sur les autres. Le calculateur CAL calcule N résidus de projection, égaux à la différence entre les projections acquises et les images simulées suivant l’équation : Le calculateur CAL enregistre dans la mémoire MEM les N projections acquises et/ou les N images simulées et/ou les N résidus de projection. Le calculateur CAL analyse ces N résidus de projection afin d’identifier des défauts inhérents à la pièce 200 inspectée. Le calculateur CAL peut fournir sur l’interface INT2 de sortie ces N résidus de projection et/ou les défauts identifiés à partir de ces N résidus de projection. Le calculateur CAL peut fournir sur l’interface INT2 de sortie un certificat de validité ou d’invalidité de la pièce, déterminé par le calculateur CAL à partir des N résidus de projection. Le calculateur CAL ou l’utilisateur peut, sur la base des N projections acquises et/ou des N images simulées et/ou des N résidus de projection, sanctionner la pièce 200 comme étant non valide car ayant trop de défauts, soit valider la pièce 200 comme étant valide car n’ayant pas trop de défauts. Suivant un mode de réalisation de l’invention, le vecteur regroupe l’un à la suite de l’autre l’ensemble des paramètres de la géométrie de projection, des paramètres de durcissement de faisceau du rayonnement dans la pièce 200 et des paramètres de diffusion Compton du rayonnement et a donc une dimension de . D’une manière correspondante, le vecteur regroupe l’un à la suite de l’autre l’ensemble des champs de sensibilité , des champs de sensibilité et des champs de sensibilité et a donc une dimension de . D’une manière correspondante, le vecteur regroupe l’un à la suite de l’autre l’ensemble des variations des paramètres de la géométrie de projection, des paramètres de durcissement de faisceau du rayonnement dans la pièce 200 et des paramètres de diffusion Compton du rayonnement et a donc une dimension de . La figure 10 montre un exemple de différence entre, d’une part l’image d’intensité acquise corrigée ayant été obtenue suivant l’invention à partir de la projection acquise , où est un réglage prescrit à partir des projections acquises, et d’autre part l’image d’intensité simulée corrigée ayant été obtenue à partir de la projection simulée pour une pale d’aube de turbine de turbomachine. La figure 11 montre un exemple de différence entre, d’une part l’image d’intensité acquise corrigée ayant été obtenue suivant l’invention à partir de la projection acquise , où est un réglage prescrit à partir des projections acquises, et d’autre part l’image d’intensité simulée ayant été obtenue à partir de la projection simulée pour un pied d’aube de turbine de turbomachine. Les figures 10 et 11 montrent que l’analyse de la pièce 200 est facilitée par l’invention prenant en compte les artéfacts. Les figure 8 et 9 montrent , sans correction sur , et où est généré avec , et . Les figures 10 et 11 montrent , sans correction sur , et où est généré avec , et . Dans le cas de la combinaison des étapes E3a, E3b et E3c, le fait de prendre en compte les connaissances sur le système d’acquisition (cabine et pièce, étape E3a) et de le coupler à une méthode algorithmique E3a, E3b, E3c de traitement d’images adaptée permet un gain significatif sur la qualité de l’information contenue dans les images produites par le système et dans celles reproduites par simulation, et donc de leur différence. Cela permet l’estimation des artefacts (1) du durcissement de faisceau et (2) de la diffusion Compton, dans le cadre de l’analyse de pièces aéronautiques à partir d’un nombre limité de projections multi-vues. L’étape E3a permet l’estimation des paramètres géométriques permettant de lier le modèle MODP de la pièce idéale aux projections acquises de la pièce à contrôler. Cette étape E3a permet de retrouver les paramètres dans le but, notamment, de les utiliser pour simuler numériquement les projections et reproduire des artéfacts dans les images simulées. L’invention permet également à l’étape E4, soit de corriger les artefacts dus au durcissement de faisceau (E3b) dans les projections acquises, soit de les reproduire dans les images simulées ; et de reproduire les artefacts dus à la diffusion Compton (E3c) dans les images simulées ou de les corriger dans les projections acquises. Le système d’acquisition d’images et la suite algorithmique fiabilisent la validation et la sanction des indications de santé matière et des indications de la géométrie tridimensionnelle de la pièce 200 à partir d’un nombre N limité de vues radiographiques par rayons X. Suivant un mode de réalisation de l’invention, le procédé, le dispositif et le programme suivant l’invention peuvent être mis en œuvre sur une ou plusieurs régions d’intérêt, sélectionnée par un module de sélection du calculateur CAL, dans la projection . Bien entendu, les modes de réalisation, caractéristiques, possibilités et exemples décrits ci-dessus peuvent être combinés l’un avec l’autre ou être sélectionnés indépendamment l’un de l’autre. Procédé de contrôle non destructif d’une pièce (200) par radiographie en transmission, comportant les étapes suivantes, exécutées par un calculateur (CAL): acquisition (E1) de N projections ( ) de la pièce (200) à partir d’un dispositif (100) de radiographie par transmission selon N angles (ANG (n) ) de vue différents et prédéterminés de la pièce (200), où N est un entier naturel prescrit, génération (E2) de N images ( ) calculées de la pièce (200) à partir d’un modèle (MODP) de référence de la pièce (200) correspondant aux N angles (ANG (n) ) de vue et à partir d’un vecteur de paramètres ( ) caractérisant la géométrie de projection de l’acquisition pour les N angles (ANG (n) ) de vue à chacune de plusieurs itérations successives, estimation (E3a, E3b, E3c), par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres ( ) rend compte du durcissement de faisceau du rayonnement dans la pièce (200) et le vecteur des paramètres ( ) caractérise la diffusion Compton du rayonnement dans la pièce (200), par minimisation du carré de la norme des différences entre les N projections ( ) acquises et les N images ( ) calculées, traitement (E4) des N projections ( ) et/ou des N images ( ) calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images ( ) calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images ( ) calculées à partir du vecteur ayant été estimé, identification (E5) de défauts de la pièce (200) par comparaison des N projections ( ) traitées aux N images ( ) calculées traitées. Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce qu’il comporte à chaque itération l’estimation (E3a) du vecteur des paramètres de la géométrie de projection de l’acquisition pour les N angles (ANG (n) ) de vue par le calculateur (CAL), cette estimation comportant : le calcul de résidus de projection à partir du vecteur des valeurs initiales , le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales où est le vecteur colonne des paramètres de la géométrie de projection, est un vecteur colonne des valeurs initiales prescrites des paramètres de la géométrie de projection, est un entier naturel désignant le numéro de l’angle (ANG (n) ) de vue et allant de 1 à N, est mis à jour selon , où est le vecteur colonne des paramètres de la géométrie de projection mis à jour, est le vecteur colonne de variation des paramètres de la géométrie de projection et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce qu’il comporte à chaque itération l’estimation (E3b) du vecteur des paramètres d’étalonnage de durcissement de faisceau du rayonnement dans la pièce (200) par le calculateur (CAL), cette estimation comportant: le calcul de résidus de projection à partir du vecteur des valeurs initiales , le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales , où est l’image obtenue en appliquant la fonction à l’intensité de chacun des pixels de , sont les N images calculées pour les N angles (ANG (n) ) de vue, sont les N projections acquises de la pièce (200), est le vecteur colonne des paramètres d’étalonnage de durcissement de faisceau, est un vecteur colonne des valeurs initiales prescrites des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce (200), est une base de fonctions de forme prescrites, est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , est mis à jour selon , où est le vecteur colonne des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce (200), ayant été mis à jour, est le vecteur colonne de variation des paramètres d’étalonnage de durcissement de faisceau de rayonnement dans la pièce (200) et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce qu’il comporte à chaque itération l’estimation (E3c) du vecteur des paramètres de diffusion Compton du rayonnement dans la pièce (200) par le calculateur (CAL), cette estimation comportant: le calcul de résidus de projection à partir du vecteur des valeurs initiales , le calcul de champs de sensibilité selon à partir du vecteur des valeurs initiales , où sont les N projections acquises de la pièce (200), sont les N images calculées pour les N angles (ANG (n) ) de vue, est l’image obtenue en convoluant des images simulées, ayant été obtenues à partir au moins des N images calculées, avec le noyau , est le vecteur colonne des paramètres de diffusion Compton du rayonnement dans la pièce (200), est un vecteur colonne des valeurs initiales prescrites des paramètres de diffusion Compton du rayonnement dans la pièce (200) est le noyau de convolution défini par est un entier naturel prescrit supérieur ou égal à 1, est un entier naturel allant de 1 à , un noyau gaussien bidimensionnel d’écart-type prescrit, est la fonction de Dirac au pixel , est mis à jour selon , où est le vecteur colonne des paramètres de diffusion Compton du rayonnement dans la pièce (200), ayant été mis à jour, est le vecteur colonne de variation des paramètres de diffusion Compton du rayonnement dans la pièce (200) et est calculé comme le minimisant le carré de la somme des différences entre les résidus de projection et le produit de par pour allant de 1 à N, où est un vecteur colonne, est la matrice des champs de sensibilité. Procédé suivant les revendications 3 et 4 prises ensemble, caractérisé en ce que est l’image simulée, obtenue en appliquant la fonction à l’intensité de chacun des pixels de l’image calculée . Procédé suivant l’une quelconque des revendications précédentes, caractérisé en ce que N est inférieur ou égal à 1000. Programme d’ordinateur, comportant des instructions de code pour la mise en œuvre des étapes suivantes d’un procédé de contrôle non destructif d’une pièce (200) par radiographie par transmission, lorsqu’il est exécuté par un calculateur (CAL) : réception (E1) de N projections ( ) de la pièce (200) à partir d’un dispositif (100) de radiographie par transmission selon N angles (ANG (n) ) de vue différents et prédéterminés de la pièce (200), où N est un entier naturel prescrit, génération (E2) de N images ( ) calculées de la pièce (200) à partir d’un modèle (MODP) de référence de la pièce (200) correspondant aux N angles (ANG (n) ) de vue et à partir d’un vecteur de paramètres ( ) caractérisant la géométrie de projection de l’acquisition pour les N angles (ANG (n) ) de vue à chacune de plusieurs itérations successives, estimation (E3a, E3b, E3c), par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres ( ) rend compte du durcissement de faisceau du rayonnement dans la pièce (200) et le vecteur des paramètres ( ) caractérise la diffusion Compton du rayonnement dans la pièce (200), par minimisation du carré de la norme des différences entre les N projections ( ) acquises et les N images ( ) calculées, traitement (E4) des N projections ( ) et/ou des N images ( ) calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images ( ) calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images ( ) calculées à partir du vecteur ayant été estimé, identification (E5) de défauts de la pièce (200) par comparaison des N projections ( ) traitées aux N images ( ) calculées traitées. Dispositif de contrôle non destructif d’une pièce (200) par radiographie en transmission, comportant : un dispositif (100) de radiographie en transmission, pour l’acquisition (E1) de N projections ( ) de la pièce (200) selon N angles (ANG (n) ) de vue différents et prédéterminés de la pièce (200), où N est un entier naturel prescrit, un calculateur (CAL) configuré pour effectuer les étapes suivantes : génération (E2) de N images ( ) calculées de la pièce (200) à partir d’un modèle (MODP) de référence de la pièce (200) correspondant aux N angles (ANG (n) ) de vue et à partir d’un vecteur de paramètres ( ) caractérisant la géométrie de projection de l’acquisition pour les N angles (ANG (n) ) de vue à chacune de plusieurs itérations successives, estimation (E3a, E3b, E3c), par les itérations successives, du vecteur à partir d’un vecteur initial et d’au moins l’un du vecteur de paramètres à partir d’un vecteur initial et du vecteur de paramètres à partir d’un vecteur initial , où le vecteur des paramètres ( ) rend compte du durcissement de faisceau du rayonnement dans la pièce (200) et le vecteur des paramètres ( ) caractérise la diffusion Compton du rayonnement dans la pièce (200), par minimisation du carré de la norme des différences entre les N projections ( ) acquises et les N images ( ) calculées, traitement (E4) des N projections ( ) et/ou des N images ( ) calculées comprenant un premier traitement et/ou un deuxième traitement, le premier traitement comportant une correction du durcissement de faisceau sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération du durcissement de faisceau sur les N images ( ) calculées à partir du vecteur ayant été estimé, le deuxième traitement comportant une correction de la diffusion Compton sur les N projections ( ) à partir du vecteur ayant été estimé ou une génération de la diffusion Compton sur les N images ( ) calculées à partir du vecteur ayant été estimé, identification (E5) de défauts de la pièce (200) par comparaison des N projections ( ) traitées aux N images ( ) calculées traitées.