Brevet FR3156526A1 — publié le 20250613

Méthode d’imagerie par tomographie à faisceau conique avec correction d’artefacts de reconstruction générés par la géométrie du faisceau Un aspect de l’invention concerne une méthode mise en œuvre par ordinateur, de détermination d’une représentation 3D d’une pièce par tomographie à faisceau conique émis depuis une source. Figure à publier avec l’abrégé : Figure 3 Méthode d’imagerie par tomographie à faisceau conique avec correction d’artefacts de reconstruction générés par la géométrie du faisceau DOMAINE TECHNIQUE DE L’INVENTION Le domaine technique de l’invention est celui du contrôle non-destructif. La présente invention concerne une méthode d’imagerie par tomographie à faisceau conique, notamment pour l’émission de rayons X ou de rayons non-collimatés dans le domaine du visible. ARRIERE-PLAN TECHNOLOGIQUE DE L’INVENTION La tomographie , par exemple à rayons X (RX) ou à rayons non-collimatés à des longueurs d’ondes visibles, est une méthode d’imagerie classiquement utilisée pour réaliser une inspection non destructive d’un corps ou d’un objet (typiquement une pièce ou 'une partie de la pièce). Cette technique utilise l'émission de RX pour obtenir une représentation 3D de l’objet. La conformité de l’objet peut ensuite être contrôlée en comparant la représentation 3D obtenue à des données de référence, telle une représentation 3D de référence de la pièce inspectée (dimensions externes, épaisseurs des parois, densité de matière, volume, etc.). Typiquement, dans un cadre industriel, la pièce à inspecter est positionnée sur un plateau tournant, situé entre la source et le détecteur de photons. A chaque rotation du plateau (tous les 0,1° par exemple), une image radiographique, également appelée « radiographie », est enregistrée. Une fois toutes les radiographies obtenues, l’ensemble des images est chargé dans un algorithme de reconstruction, lequel sert à reconstruire le volume numérique, i.e., la représentation 3D, de la pièce. Ce volume numérique est ensuite utilisé pour contrôler la pièce inspectée en la comparant au modèle numérique de référence. La qualité du contrôle dépend donc directement de la qualité de la reconstruction du volume numérique. Parmi les techniques de tomographie existantes, la tomographie par faisceau conique (en anglais, « Cone Beam »), aussi appelée « imagerie volumétrique par faisceau conique », repose sur l’utilisation d’un faisceau divergent (RX ou dans le domaine visible), formant un cône, pour illuminer la pièce et acquérir les radiographies. Cette approche permet d’acquérir les images radiographiques rapidement et avec une haute résolution, grâce au grandissement de l’objet sur le détecteur par la divergence du faisceau. Elle engendre, cependant, des lacunes de mesure dans des régions éloignées de l’axe optique du faisceau. Ces lacunes génèrent, dans ces régions, des artefacts lors de la reconstruction (en anglais, « cone beam artefacts ») qui dégradent la qualité du contrôle et mènent à une évaluation erronée des propriétés de la pièce. Ces artefacts peuvent tout autant altérer cette évaluation des propriétés de la pièce à l’échelle globale, par exemple une mauvaise estimation du volume total de la pièce, qu’à l’échelle locale, par exemple une mauvaise estimation d’un diamètre d’un cylindre interne. Outre ceci, si un artefact est situé dans une zone de recalage, avec le modèle numérique de référence, l’analyse peut se trouver totalement faussée. Afin de supprimer la présence des artefacts liés à la divergence du faisceau, il est connu d’appliquer, lors de l’acquisition des radiographies, une trajectoire à la pièce dans le champ de vision du faisceau, en plus de sa rotation sur le plateau. Par exemple, il est connu de compenser les lacunes de mesure par application d’une trajectoire hélicoïdale, élevant et/ou abaissant la pièce et permettant de l’illuminer sous plusieurs angles du faisceau. Toutefois, cette approche accroît significativement le temps d’acquisition et requiert un équipement de mesure adapté ainsi qu’un algorithme de reconstruction spécifique. En conséquence, la tomographie par faisceau conique est peu utilisée en production pour vérifier la conformité métrologique des pièces dans leur intégralité, d’une part, car la présence de ces artefacts ne permet pas d’obtenir une résolution suffisante du volume 3D reconstruit et, d’autre part, parce que les temps d’acquisition des approches existantes pour compenser ces artefacts ne sont pas compatibles avec les cadences de production. Il existe donc un besoin pour corriger ou diminuer les artefacts liés à la divergence du faisceau afin d’évaluer de façon précise et rapide la conformité d’une pièce. L’invention offre une solution aux problèmes évoqués précédemment, en permettant la mise en œuvre d’une méthode d’imagerie par tomographie par faisceau conique permettant la diminution ou la suppression d’artefacts liés à la géométrie conique des faisceaux de rayons et engendrés lors de la reconstruction du volume numérique à partir des radiographies acquises. Un premier aspect de l’invention concerne une méthode mise en œuvre par ordinateur, de détermination d’une représentation 3D d’une pièce par tomographie à faisceau conique émis depuis une source, la méthode comprenant : Recevoir un modèle numérique 3D de la pièce comprenant une pluralité de mailles ; Obtenir un ensemble de configurations d’irradiation de la pièce, chaque configuration d’irradiation étant associée une pluralité respective de positions prédéterminées successives de la pièce par rapport à la source ; Pour chaque configuration d’irradiation de l’ensemble de configurations d’irradiation : Déterminer une pluralité de projections du modèle numérique de la pièce, chaque projection étant respectivement associée à une position prédéterminée de la pluralité de positions prédéterminées successives associée à ladite configuration d’irradiation, chaque projection comprenant un ensemble de portions, chaque portion étant associée à une maille respective du modèle numérique 3D ; Pour chaque maille du modèle numérique 3D : Déterminer un nombre de projections parmi la pluralité de projections déterminée pour lesquelles la portion associée à la maille respecte un critère de qualité prédéfini, le critère de qualité étant lié à un risque de présence d’un artefact dans une représentation 3D de la pièce construite à partir des projections déterminées, le critère de qualité étant une condition de Tuy ; Si le nombre de projections déterminé dépasse un seuil prédéfini, associer à la maille une première valeur d’un indicateur de qualité ; sinon, associer à la maille une deuxième valeur de l’indicateur de qualité ; Déterminer une première configuration d’irradiation et une deuxième configuration d’irradiation parmi l’ensemble de configurations d’irradiation en minimisant un nombre de mailles pour lesquelles la valeur de l’indicateur de qualité calculé pour la première configuration et la valeur de l’indicateur de qualité calculé pour la deuxième configuration sont égales à la deuxième valeur ; Recevoir une première pluralité d’images radiographiques 2D de la pièce lorsque la pièce est dans la première configuration d’irradiation et une deuxième pluralité de d’images radiographiques 2D de la pièce lorsque la pièce est dans la deuxième configuration d’irradiation ; Construire la représentation 3D de la pièce à partir de la première pluralité d’images radiographiques 2D et de la deuxième pluralité d’images radiographiques 2D. On entend par « configuration d’irradiation » un ensemble de paramètres décrivant, notamment, un emplacement initial de la pièce. L’emplacement initial de la pièce est défini par une position initiale, une inclinaison initiale et/ou une orientation initiale. L’ensemble de paramètres décrit également le déplacement de la pièce inspectée, depuis l’emplacement initial ou autour de celui-ci, par rapport à l’appareil servant à l’imagerie de la pièce par tomographie, en particulier par rapport à la source de photons. Ce déplacement sert à modifier le positionnement relatif de la pièce par rapport à la source pour réaliser une mesure dans des conditions spécifiques. Les conditions spécifiques de la configuration d’irradiation sont définies par la valeur de chaque paramètre de son ensemble de paramètres. Une mesure est donc associée à une configuration d’irradiation et comprend l’acquisition d’une ou plusieurs images radiographiques de la pièce qui se déplace par rapport à la source depuis son emplacement initial selon les conditions spécifiques spécifiées par ladite configuration d’irradiation. Dans certains modes de réalisation, le déplacement est appliqué à la source plutôt qu’à la pièce. L’ensemble de paramètres définissant la configuration d’inspection peut également comprendre des paramètres de la source (tels que sa puissance, sa taille, la différence de potentiel du tube et l’énergie des photons émis), la composition de la pièce cible, des paramètres de filtrages, des paramètres de l’écran servant à la détection (tels que sa taille en pixels, le type de détection, l’épaisseur et la composition de la zone sensible), etc. Dans la suite, le terme « image radiographique » se rapporte à une image radiographique acquise au cours d’une expérimentation, par exemple acquise par un appareil d’imagerie de tomographie à faisceau conique par rayons X ou par rayons émis dans le domaine du visible (c’est-à-dire dont la longueur d’onde d’émission dans le vide de chaque photon est typiquement comprise entre 380 nm et 780 nm). On entend par « modèle numérique 3D » une représentation numérique en trois dimensions (3D) de la pièce, par exemple sous forme d’une conception assistée par ordinateur (CAO), construite au moyen d’un outil de simulation. Le modèle numérique 3D comprend un maillage formé d’une pluralité de mailles adjacentes. Chaque maille est délimitée par au moins trois nœuds, qui correspondant aux « sommets » de la maille. On entend par « projection » un projeté sur un plan parallèle à l’écran de détection se trouvant entre la source et l’écran de détection, préférentiellement entre la pièce et l’écran de détection, typiquement sur l’écran de détection ou à l’emplacement d’un barycentre de la pièce. Cette projection comprend le projeté des mailles du maillage sur ce plan, pour une position prédéterminée de la pièce au cours de son déplacement par rapport à la source depuis son emplacement initial, pour la configuration d’irradiation associée. Il s’agit donc d’une représentation géométrique d’une image radiographique qui serait acquise par l’écran de détection dans la configuration d’irradiation concernée et pour une position spécifique de la pièce par rapport à la source. Une projection est donc construite pour chaque position prédéterminée d’une pluralité de positions prédéterminée successives de la pièce au cours de son déplacement. Chaque position prédéterminée permet d’irradier la pièce suivant un angle de vue différent depuis la source. On entend par « artefact » une anomalie ou une erreur de reconstruction dans une représentation 3D de la pièce faussant la lecture de la représentation 3D. La représentation 3D de la pièce est une reconstruction numérique de la pièce à partir de d’images radiographiques acquises de la pièce pour former une reproduction numérique en trois dimensions de la pièce imagée. Le modèle numérique sert également au recalage de la reconstruction volume afin de l’analyser dans un espace connu. On entend par « critère de qualité » une propriété associée à chaque projection et pour une maille donnée. Cette propriété estime le risque de produire un artefact dans la reconstruction de la représentation 3D pour cette maille, lorsque la pièce est dans une configuration d’irradiation spécifique. On entend par « indicateur de qualité » une propriété associée à chaque maille. Cette propriété indique de façon binaire si un artefact est présent de manière quasi-certaine pour cette maille dans la reconstruction numérique. Grâce à l’invention, il est possible d’améliorer la tomographie de la pièce en minimisant le nombre d’artefacts de reconstruction présents dans la tomographie et liés à la divergence du faisceau. En effet, la méthode ci-dessus permet de déterminer les configurations d’irradiation optimales pour inspecter la pièce de sorte à minimiser le nombre d’artefacts lors de la reconstruction dans la tomographie. Des images radiographiques sont ensuite obtenues pour ces configurations optimales déterminées. La reconstruction de la tomographie est ainsi réalisée en extrayant, pour chaque maille, une donnée ou plusieurs données des images radiographiques associées à l’une des configurations optimales. Il est ainsi possible de maitriser et d’améliorer la fiabilité et la précision de la reconstruction en chaque maille, et d’augmenter ainsi la qualité globale de la représentation 3D de l’ensemble de la pièce. A cet effet, l’utilisation de la condition de Tuy permet de garantir une estimation fiable de la présence d’un artefact dans la tomographie. La méthode est, en outre, facile et rapide à mettre en œuvre puisqu’elle ne requiert, a minima, que la détermination de projections selon deux configurations d’irradiation. La méthode est donc compatible avec les cadences élevées des services de production. De plus, la méthode est robuste, fiable et produit des tomographies avec une précision élevée. Par ailleurs, la méthode n’impose pas de restriction quant aux propriétés de la pièce inspectée, notamment concernant ses dimensions géométriques. Avantageusement, la méthode ci-dessus peut être utilisée afin de déterminer plus de deux configurations optimales, par exemple trois ou plus configurations optimales, dans le but de réduire encore le nombre d’artefacts dans la représentation 3D, voire de garantir l’absence d’artefact. Outre les caractéristiques qui viennent d’être évoquées, la méthode selon le premier aspect de l’invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. Dans un mode de réalisation, la représentation 3D de la pièce comprend le modèle numérique 3D et une pluralité de valeurs d’amplitude, dites valeurs de mailles, chaque valeur de maille de la pluralité de valeurs de mailles étant respectivement associée à une maille de la pluralité de mailles, chaque valeur de maille étant déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image radiographique 2D de la première pluralité d’images radiographiques 2D et/ou de la deuxième pluralité d’images radiographiques 2D. On entend par « valeur d’amplitude » une valeur numérique affectée à un niveau de couleur, typiquement un niveau d’une échelle en niveaux de gris. Cette valeur correspond à une grandeur physique qui peut être directement ou indirectement déterminée à partir de la mesure effectuée via la tomographie mise en œuvre. La nature de cette grandeur physique dépend de l’application concernée et peut ainsi être obtenue par des méthodes connues en soi. Typiquement, la valeur d’amplitude à la dimension d’une énergie déposée, laquelle est fonction, notamment, du nombre de photons détectés en chaque pixel, de l’énergie de chaque photon détecté, de la fonction de transfert du capteur, etc. En d’autres termes, chaque maille de la représentation 3D est affectée à un niveau de l’échelle de couleur choisie, où ce niveau correspond à une valeur d’amplitude, laquelle dépend des propriétés du matériau telles que sa densité et son absorption. Dans un mode de réalisation, la méthode selon le premier aspect comprend en outre : Détection d’une non-conformité de la pièce par comparaison de la représentation 3D de la pièce avec une représentation 3D de référence ; Emission d’une alerte relative à la non-conformité détectée. On entend par « représentation 3D de référence » une représentation 3D comprenant le modèle numérique 3D et une pluralité de valeurs de référence d’amplitude de laquelle chacune des valeurs de référence est associée à une maille de la pluralité de mailles du modèle numérique 3D. Ces valeurs d’amplitude de référence sont des valeurs d’amplitude provenant, par exemple, d’une simulation et/ou d’une expérimentation sur une pièce saine, c’est-à-dire sans défaut. Dans un mode de réalisation, la détermination de la première configuration d’irradiation et de la deuxième configuration d’irradiation comprend : Déterminer, pour chaque paire de configurations d’irradiation d’une pluralité de paires de configurations d’irradiation, le nombre de mailles pour lesquelles la valeur de l’indicateur de qualité calculé de chacune des configurations d’irradiation de ladite paire de configurations d’irradiation est égale à la deuxième valeur ; Sélectionner, parmi la pluralité de paires de configurations d’irradiation, la paire de configurations d’irradiation pour laquelle le nombre de mailles est minimal, la première configuration d’irradiation étant l’une des configurations d’irradiation de la paire de configuration d’irradiation sélectionnée, et la deuxième configuration d’irradiation étant l’autre configuration d’irradiation de la paire de configuration d’irradiation sélectionnée. Il est ainsi possible de déterminer le meilleur couple de configurations d’irradiation, parmi tous les couples possibles ou parmi un sous-ensemble de ces couples, pour maximiser la qualité de la représentation 3D. Ce couple est donc celui qui minimise le nombre d’artefacts de reconstruction liés à la divergence du faisceau. Dans un mode de réalisation, pour chaque maille de la pluralité de mailles, la détermination de la valeur de maille comprend : Lorsqu’au moins un nombre parmi le nombre de projections déterminé pour la première configuration d’irradiation et le nombre de projections déterminé pour la deuxième configuration dépasse le seuil prédéfini, sélectionner, parmi la première configuration d’irradiation et la deuxième configuration d’irradiation, la configuration d’irradiation pour laquelle le nombre de projections déterminé dépassant le seuil prédéfini pour ladite maille est le plus grand, la valeur de maille étant déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image radiographique 2D de la pluralité d’images radiographiques 2D associée à la configuration d’irradiation sélectionnée ; Lorsque le nombre de projections déterminé pour la première configuration d’irradiation et le nombre de projections déterminé pour la deuxième configuration ne dépassent pas le seuil prédéfini, la valeur de maille est déterminée à partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D et d’au moins une valeur d’amplitude provenant d’au moins une image de la deuxième pluralité d’images radiographiques 2D. Il est ainsi garanti que de la représentation 3D est construite de façon optimale, en utilisant les données acquises de l’une ou l’autre des première et deuxième pluralité d’images radiographiques qui garantissent l’absence d’artefact. Dans un mode de réalisation, la détermination de la première configuration d’irradiation et de la deuxième configuration d’irradiation comprend : Sélectionner comme première configuration d’irradiation la configuration d’irradiation ayant un plus petit nombre de mailles pour lesquelles la valeur de l’indicateur de qualité calculé est égale à la deuxième valeur parmi l’ensemble de configurations d’irradiation ; Sélectionner comme deuxième configuration d’irradiation la configuration d’irradiation ayant un plus grand nombre de mailles corrigées parmi l’ensemble de configurations d’irradiation, le nombre de mailles corrigées étant le nombre de mailles pour lesquelles : La valeur de l’indicateur de qualité calculé pour la première configuration d’irradiation est égale à la deuxième valeur ; et La valeur de l’indicateur de qualité calculé pour la deuxième configuration d’irradiation est égale à la première valeur. Il est ainsi possible de déterminer la deuxième configuration de sorte à maximiser la qualité de la représentation 3D, pour chaque forme géométrique, en fonction de la qualité de la première configuration d’irradiation. Ce mode de réalisation fournit une solution sous-optimale mais permet de réduire le coût calculatoire de la méthode proposée. En effet, il n’est pas besoin de tester tous les couples de configurations de l’ensemble de configurations d’irradiation. La première configuration est donc la « meilleure » configuration, et la deuxième configuration est celle qui « complète » le mieux la première configuration. Dans un mode de réalisation, la valeur de maille est déterminée : A partir d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D si le nombre de projections déterminé pour ladite maille dépasse le seuil prédéfini lorsque la pièce est dans la première configuration d’irradiation ; A partir d’au moins une valeur d’amplitude provenant d'au moins une image de la deuxième pluralité d’images radiographiques 2D si le nombre de projections déterminé pour ladite maille ne dépasse pas le seuil prédéfini lorsque la pièce est dans la première configuration d’irradiation et si le nombre de projections déterminé pour ladite maille dépasse le seuil prédéfini lorsque la pièce est dans la deuxième configuration d’irradiation ; A partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D et d’au moins une valeur d’amplitude provenant d'au moins une image de la deuxième pluralité d’images radiographiques 2D si le nombre de projections déterminé pour ladite maille ne dépasse pas le seuil prédéfini lorsque la pièce est dans la première configuration d’irradiation et si le nombre de projections déterminé pour ladite maille ne dépasse pas le seuil prédéfini lorsque la pièce est dans la deuxième configuration d’irradiation. Dans un mode de réalisation, la méthode selon le premier aspect comprend en outre : recaler la première pluralité d’images radiographiques 2D et la deuxième pluralité d’images radiographiques 2D avec le modèle numérique 3D de la pièce à partir d’une pluralité de mailles de recalage, la construction de la représentation 3D étant effectuée à partir des images recalées. Dans un mode de réalisation, la première configuration d’irradiation et de la deuxième configuration d’irradiation sont déterminés parmi un sous-ensemble de configurations d’irradiation, le sous-ensemble de configurations d’irradiation comprenant une ou plusieurs configurations d’irradiation de l’ensemble de configurations d’irradiations pour lesquelles la valeur de l’indicateur de qualité calculé est égale à la première valeur pour une ou plusieurs mailles de recalage de la pluralité de mailles de recalage. Il est ainsi possible de garantir la fiabilité du recalage puisque les première et deuxième configurations d’irradiation sont choisies parmi des configurations d’irradiation qui ne génèrent pas d’artefact pour les mailles du modèle numérique 3D utilisées pour le recalage. Un deuxième aspect de l’invention concerne un système de détermination d’une représentation 3D d’une pièce par tomographie à faisceau conique émis depuis une source, le système étant configuré pour mettre en œuvre la méthode selon le premier aspect. Un autre aspect de l’invention concerne un produit-programme d’ordinateur comprenant des instructions qui, quand le programme est exécuté sur un ordinateur, conduisent celui-ci à mettre en œuvre les étapes de la méthode selon le premier aspect. Un autre aspect de l’invention concerne un support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre les étapes de la méthode selon le premier aspect. L’invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l’examen des figures qui l’accompagnent. BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l’invention. La FIG. 1 est un schéma synoptique illustrant l’enchaînement des étapes de la méthode selon l’invention. La FIG. 2 est une représentation schématique d’un appareil d’imagerie de tomographie à faisceau conique et d’un ordinateur pour mettre en œuvre la méthode de la FIG. 1 . La FIG. 3 est une représentation schématique de l’acquisition d’une image radiographique d’une pièce par l’appareil d’imagerie de la FIG. 2 . La FIG. 4 est une représentation schématique d’un modèle numérique 3D de la pièce, selon l’invention. La FIG. 5 est une représentation schématique d’une image radiographique selon l’invention. La FIG. 6 illustre les grandeurs physiques utilisées pour déterminer le critère de qualité de l’imagerie selon l’invention. DESCRIPTION DETAILLEE Sauf précision contraire, un même élément apparaissant sur des figures différentes présente une référence unique. L’invention concerne une méthode d’imagerie par tomographie à faisceau conique, par exemple de rayons X (RX) ou de rayons émis dans le domaine visible, l’image obtenue pouvant être par exemple utilisée pour l’inspection d’une pièce. La méthode sert à améliorer la qualité de l’imagerie en réduisant ou supprimant les artefact liés à la divergence du faisceau. A cet effet, plusieurs configurations expérimentales, par la suite appelées configurations d’irradiation, sont évaluées afin de déterminer, avant d’effectuer l’acquisition, des configurations d’irradiation, dites « optimales », maximisant la qualité de l’inspection, en minimisant le nombre d’artefacts de reconstruction. Une acquisition est ensuite effectuée pour chacune des configurations d’irradiation optimales et une représentation 3D de la pièce est construite à partir des images radiographiques acquises pour ces configurations d’irradiation optimales. La description ci-dessous présente le cas particulier où deux configurations optimales sont déterminées, mais l’invention ne se limite pas à ce cas de figure et s’étend à plus de deux configurations optimales. Un premier aspect de l’invention concerne donc une méthode de détermination d’une représentation 3D d’une pièce par tomographie à faisceau conique. Le faisceau conique comprend un axe principal et un angle solide, lequel définit la divergence du faisceau. Cette méthode 100 est représentée de façon schématique sur la FIG. 1 . La méthode 100 est mise en œuvre par un ordinateur 10, représenté sur la FIG. 2 . L’ordinateur 10 est, par exemple, intégré à un appareil d’imagerie 20 par tomographie à faisceau conique (RX ou dans le visible), ou relié (par un support de communication 30 filaire ou sans fil, par exemple) à un tel appareil 20 afin de lui transmettre des données de configurations optimales à l’appareil 20 ou de recevoir des d’images acquises au moyen de l’appareil d’imagerie 20. L’ordinateur 10 comprend dans sa mémoire 11 des instructions qui, lorsqu’exécutées par son processeur 12, mettent en œuvre les étapes de la méthode 100. L’appareil d’imagerie 20 comprend ou est relié à une source 21 configurée pour émettre un flux de photons en un faisceau conique. L’appareil 20 comprend ou est relié à un écran 22 configuré pour capter le faisceau conique de photons émis par la source. L’appareil 20, illustré sur la FIG. 3 , est donc configuré pour piloter la source 21 de sorte à émettre le faisceau conique, notamment vers la pièce 40 à imager, et pour acquérir des échantillons liés au flux de photons capté par l’écran 22, notamment après que le faisceau a traversé la pièce 40. Les échantillons du faisceau forment sur l’écran 22 le projeté 40’ de la pièce, que l’on appelle également l’image radiographique 40’. L’appareil d’imagerie 20 est en outre configuré pour modifier la position et/ou l’inclinaison de la pièce 40 par rapport à la source 21. L’appareil 20 est donc configuré pour translater et/ou incliner la pièce 40 ou la source 21 par rapport à la source 21 ou la pièce 40, respectivement. Cette modification peut également ou alternativement être définie par rapport à l’écran 22. La modification de la position et/ou de l’inclinaison est, par exemple, comprise dans la ou les configurations optimales transmises par l’ordinateur 10. La pièce 40 à imager est n’importe quelle pièce non-opaque aux photons émis, par exemple non-opaque aux RX ou aux photons émis dans le domaine visible, dont on veut obtenir une représentation 3D, en particulier pour en contrôler la conformité, par exemple après fabrication ou au cours de sa durée de vie. La première étape 110 de la méthode 100 est une étape de réception d’un modèle numérique 3D de la pièce 40. Le modèle numérique 3D 50, illustré sur la FIG. 4 , comprend une pluralité de mailles 51. La pluralité de mailles 51 définit la géométrie interne et externe de la pièce. C’est-à-dire que la géométrie de la pièce est reproduite par la pluralité de mailles 51. Le modèle numérique 3D 50 est, par exemple, une conception assistée par ordinateur (CAO) de la pièce 40 produite par un opérateur au moyen d’un outil de simulation. Chaque maille 51 est un polygone, par exemple, un triangle ou un parallélépipède. Chaque sommet d’une maille 51 est associé à un nœud 52 de la pluralité de maille et un nœud 52 peut être commun à l’un des sommets de plusieurs mailles 52. La deuxième étape 120 est ensuite une étape d’obtention d’un ensemble de configurations d’irradiation de la pièce 40. Chaque configuration d’irradiation comprend un emplacement initial de la pièce par rapport à la source 21 et un déplacement de la pièce 40 par rapport à la source 21 depuis cet emplacement initial ou autour de cet emplacement initial. Le déplacement de la pièce 40 sert à modifier le positionnement relatif de la pièce 40 par rapport à la source 21 de sorte que la pièce 40 soit irradiée par la source 21 suivant plusieurs angles de vue différents, à partir de l’emplacement initial. Un intérêt est d’être en d’acquérir une ou plusieurs images radiographiques de la pièce 40 dans chacun de ces différents angles de vue. Le terme « initial » se réfère donc à l’emplacement relatif à la source 21 dans lequel se trouve la pièce 40 avant que le déplacement ne soit appliqué. Le déplacement comprend une translation et/ou une rotation de la pièce 40 par rapport à la source 21. Alternativement, la pièce 40 reste à son emplacement initial et c’est la source 21 qui est déplacée autour de la pièce 40 pour acquérir les images suivant les différents angles de vues. Le déplacement de la configuration d’irradiation concerne alors le mouvement de la source 21 par rapport à la pièce 40 depuis l’angle de vue dans lequel la pièce 40 se trouvait à son emplacement initial. La pièce et le dispositif peuvent être disposés dans un espace associé à un repère à trois dimensions, classiquement notées x, y et z, comme représenté sur la FIG. 2 . Dans certains modes de réalisation, l’origine du repère peut être la source 21. L’emplacement initial de la pièce 40 est défini par une position initiale et une inclinaison initiale de la pièce 40. La position initiale indique la position où se trouve la pièce 40, avant d’appliquer le déplacement, suivant un ou plusieurs des axes du repère. L’inclinaison initiale indique l’angle d’inclinaison de la pièce 40, ou d’un axe de la pièce 40, par rapport à un ou plusieurs des axes x, y et z. Préférentiellement, pour chaque axe, la position initiale de chaque configuration d’irradiation est comprise dans un intervalle de valeurs. De même, pour chaque axe, l’inclinaison initiale de chaque configuration d’irradiation est comprise dans un intervalle de valeurs. Ces intervalles de valeurs définissent l’ensemble des positions et inclinaisons admissibles pour la pièce à imager. Chaque intervalle est défini par deux bornes comprenant une valeur supérieure et une valeur inférieure. Les valeurs des bornes dépendent notamment des contraintes d’accessibilités de l’appareil d’imagerie 20, en particulier de la distance entre la source 21 et l’écran 22, mais également des caractéristiques de la pièce 40 telles que ses dimensions et sa masse. Chaque configuration d’irradiation indique donc une position et une inclinaison spécifiques dans lesquelles est initialisé l’emplacement de la pièce 40 aux fins de l’imagerie. L’ensemble des configurations d’irradiation est, par exemple, construit par un opérateur, au moyen d’un outil numérique configuré pour générer lesdites configurations d’irradiation. Ces configurations peuvent être générées par incrémentations successives de la valeur de la donnée de position et/ou de la donnée d’inclinaison dans un ou plusieurs des intervalles de valeurs. Ces configurations peuvent également être automatiquement générées par l’outil numérique générant ces configurations en fonction de paramètres renseignés par l’opérateur. Le déplacement comprend une ou plusieurs rotations de la pièce par rapport à la source, formant ainsi une pluralité de positions prédéterminées successives, elle-même formant une trajectoire. Le positionnement relatif de la pièce 40 par rapport à la source 21 est ainsi modifié par cette trajectoire. La pluralité de positions prédéterminées successives comprend l’emplacement initial. Dans un mode de réalisation, ainsi qu’illustré sur la FIG. 3 , le déplacement comprend une rotation R40 de la pièce 40 autour d’un axe rotation R. Le déplacement peut alternativement ou également comprendre une rotation R21 de la source autour de l’axe de rotation R. L’axe de rotation R est compris dans un plan qui intersecte au moins une fois la pièce 40. Dans ce mode de réalisation, chaque position prédéterminée du déplacement est définie par un angle de rotation, formant ainsi une pluralité d’angles de rotation. Chaque angle de rotation correspond à l’angle formé par rapport à l’axe de rotation R entre les positions respectives d’un point de la pièce pour deux positions prédéterminées successives. La pièce, ou la source, est donc, dans chaque configuration d’irradiation, incrémentalement orientée de façon à successivement changer l’angle de vue de la pièce 40 depuis la source 21. Dans certains modes de réalisation, le déplacement est une rotation pure, telle que la rotation R40 de la pièce 40 autour d’un axe rotation R. Les valeurs des angles de la pluralité d’angles de rotation sont comprises dans un intervalle angulaire délimité par un angle maximal et un angle minimal. L’angle minimal est, préférentiellement, supérieur ou égal à 0°, et l’angle maximal, plus grand que l’angle minimal, est préférentiellement inférieur ou égal à 360°. A titre d’exemple, la pluralité d’angles d’incidence comprend des angles compris entre l’angle minimal et l’angle maximal et répartis dans l’intervalle angulaire par un pas angulaire. Le pas angulaire est typiquement compris entre 0° et 90°. Par exemple, le pas angulaire est égal à 0,1°, 1°, 2°, 5°, etc. Dans certains cas, la pièce 40 est solidaire d’un plateau mobile 24 servant à déplacer la pièce 40 en translation selon un ou plusieurs des axes du repère, ainsi qu’en rotation autour d’un axe de rotation du plateau mobile 24. La position initiale et l’inclinaison initiale s’appliquent alors au plateau 24. La pluralité d’angles de rotation sert au pilotage du plateau mobile 24 et indiquent l’écart angulaire à appliquer à la rotation R24 du plateau 24 autour de son axe de rotation afin de successivement déplacer la pièce 40 d’une position prédéterminée à une autre parmi la pluralité de positions prédéterminées successives. La troisième étape 130 est ensuite une étape de détermination d’une pluralité de projections du modèle numérique de la pièce. La troisième étape 130 est mise en œuvre pour chaque configuration d’irradiation de l’ensemble de configurations d’irradiation. Chaque pluralité de projections peut être obtenue par une approche connue en soi, typiquement via un outil numérique dédié, par exemple via une construction géométrique. Chaque projection est construite en positionnant la pièce 40 dans l’une des positions prédéterminées de la pluralité de positionnées prédéterminées successives, pour la configuration d’irradiation concernée. Chaque projection est associée à l’angle de rotation correspondant à la position prédéterminée dans laquelle la pièce se trouve. Lorsque la pièce est à l’emplacement initial, l’angle de rotation correspondant est égal à 0. Chaque projection appartient, par exemple, au plan formé par la surface réceptrice de l’écran et orthogonal à la l’axe principal C du faisceau. Chaque projection est donc une image à deux dimensions. Ainsi qu’illustré sur la FIG. 5 , chaque projection 50’ comprend un ensemble de portions 51’ dont une portion 51’ correspond à une partie de ladite projection. Les portions 51’ sont définies telles que chaque portion 51’ de la projection correspond à au moins une maille 51 du modèle numérique 3D 50. Chaque portion 51’ est donc associée à une maille 51 respective du modèle numérique 3D 50. La quatrième étape 140 est ensuite une étape de détermination d’un nombre de projections qui respectent un critère de qualité prédéfini. La quatrième étape 140 est mise en œuvre pour chaque configuration d’irradiation de l’ensemble de configurations d’irradiation et pour chaque maille 51 du modèle numérique 3D 50. Pour la maille concernée, une projection 50’ respecte le critère de qualité donnée lorsque la portion 51’ de ladite projection qui est associée à la maille concernée respecte ledit critère de qualité. Le critère de qualité indique si, pour la maille 51 concernée et pour la configuration d’inspection concernée, la représentation 3D a un risque de comprendre un artefact lié à la divergence du faisceau conique. Ce critère permet donc de prédire la fiabilité, et donc la qualité, de la représentation 3D dans le cas où elle est construite à partir d’images radiographiques 40’ acquises lorsque la pièce 40 est dans la configuration concernée. Dans certains modes de réalisation, le critère de qualité est basé sur une condition de Tuy. Par exemple, le critère de qualité peut être une condition de Tuy-Smith pour tout angle formé par le plan source, c’est-à-dire un plan comprenant le rayon incident au point d’impact de ce rayon sur la pièce, et la tangente à la face de la pièce illuminée par le faisceau, telle que , où : , et sont les coordonnées de la position de la source selon les trois dimensions du repère, respectivement ; est la distance entre la source et l’axe de rotation de la pièce évaluée dans le plan horizontal ; , et sont les coordonnées de la position de la face, i.e., la maille considérée, de la pièce illuminée par le faisceau selon les trois dimensions du repère, par exemple les coordonnées du barycentre de ladite face ; , et sont les coordonnées d’un vecteur normal à la face de la pièce illuminée par le faisceau selon les trois dimensions du repère. Ces différentes grandeurs sont définies sur la FIG. 6 . La condition de Tuy exprime le risque de produire un artefact dans la représentation 3D pour la portion 51’ de projection 50’ concernée en fonction de son éloignement par rapport à l’axe principal du faisceau. Plus la portion 51’ associée à la maille 51 concernée est éloignée de l’axe principal, et plus le risque de former un artefact dans la représentation 3D est élevé. Plus précisément, la condition de Tuy est satisfaite lorsque chaque plan sécant de la pièce 40 imagée est aussi sécant d’un rayon du faisceau au moins une fois. La portion 51’ respecte donc le critère de qualité lorsqu’elle respecte le critère de Tuy. Respectivement, lorsque la portion 51’ ne respecte pas le critère de Tuy, elle ne respecte pas le critère de qualité. La cinquième étape 150 est ensuite une étape d’association de chaque maille 51 avec une valeur d’un indicateur de qualité. La cinquième étape 150 est donc mise en œuvre pour chaque maille 51 du modèle numérique 3D et pour chaque configuration d’irradiation. L’indicateur de qualité prend une première valeur ou une deuxième valeur. Par exemple, l’indicateur de qualité est binaire. La maille 51 concernée pour la configuration concernée est associée à la première valeur ou la deuxième valeur de la manière suivante : La maille 51 est associée à la première valeur lorsque le nombre de projections 50’ qui respectent le critère de qualité prédéfini, pour la maille 51 concernée et pour la configuration concernée, est supérieur à un seuil prédéfini ; ou La maille 51 est associée à la deuxième valeur lorsque le nombre de projections 50’ qui respectent le critère de qualité prédéfini, pour la maille 51 concernée et pour la configuration concernée, est inférieur au seuil prédéfini. Autrement dit, la maille 51 est associée à la première valable lorsque, pour cette maille, la représentation 3D est « acceptable » en termes de nombre de projections respectant le critère de qualité. Une maille 51 est donc acceptable lorsque le nombre de projections respectant le critère de qualité pour cette maille est suffisamment élevé. Respectivement, la maille 51 est associée à la deuxième valeur lorsque, pour cette maille, la représentation 3D est « non-acceptable » en termes de nombre de projections respectant le critère de qualité. Une maille 51 est donc non-acceptable lorsque le nombre de projections respectant le critère de qualité pour cette maille n’est pas suffisamment élevé. Le seuil prédéfini est un nombre strictement positif, préférentiellement entier, défini en fonction des propriétés de la pièce, de contraintes expérimentales et/ou d’exigence du cahier des charges pour la fabrication et/ou l’utilisation de ladite pièce. Alternativement, le seuil prédéfini correspond à un pourcentage du nombre total de projections pour la maille 51 concernée et pour la configuration concernée. Par exemple ce seuil prédéfini est égal à 90 % de ce nombre total de projections. Le seuil prédéfini est, par exemple défini par un opérateur. Le seul prédéfini peut également être défini en fonction de l’expérience métier de l’opérateur. Lorsque le nombre de projections qui respectent le critère de qualité prédéfini dépasse ce seuil, il est quasiment garanti que la reconstruction 3D ne comporte pas d’artefact pour la maille concernée et pour la configuration concernée. La première et la deuxième valeur peuvent être un nombre, un caractère ou une suite de caractères. Par exemple, la première valeur peut être 1 ou « True », et la deuxième valeur peut être 0 ou « False », ou inversement. D’autres valeurs ou caractères peuvent être utilisés. La sixième étape 160 est ensuite une étape de détermination d’une première configuration d’irradiation et d’une deuxième configuration d’irradiation, sélectionnées parmi l’ensemble de configurations d’irradiation. La première et la deuxième configurations d’irradiation sont sélectionnées pour maximiser la qualité de la reconstruction 3D, c’est-à-dire à minimiser le nombre d’artefacts de reconstruction liés à la géométrique du cône de diffraction. Chaque maille 51 peut être : Acceptable pour la première configuration et pour la deuxième configuration ; Non-acceptable pour la première configuration et pour la deuxième configuration ; Acceptable pour l’une des deux configurations et non-acceptable pour l’autre de ces deux configurations. Ainsi, lorsque la maille 51 est acceptable pour au moins une des deux configurations, la représentation 3D ne présentera pas d’artefact pour cette maille 51 pour ladite configuration. Au contraire, lorsque la maille 51 est non-acceptable pour les deux configurations, la représentation 3D présentera un artefact pour cette maille 51 qu’importe la configuration. Les première et deuxième configurations sont donc déterminées de sorte à être complémentaires l’une de l’autre afin de minimiser le nombre de mailles 51 qui sont non-acceptables aussi bien lorsque la pièce 40 est dans la première configuration que lorsque la pièce 40 est dans la deuxième configuration. Ces deux configurations sont donc sélectionnées de sorte à minimiser le nombre de mailles du modèle numérique 3D associées à la deuxième valeur de l’indicateur de qualité lorsque : La pièce 40 est positionnée dans la première configuration d’irradiation ; et La pièce est positionnée dans la deuxième configuration d’irradiation. En d’autres termes, les deux configurations sont déterminées de sorte à maximiser le nombre de mailles 51 acceptables lorsque la pièce 40 est dans au moins une des première et deuxième configurations. Il est donc possible que la première configuration et la deuxième configuration présentent individuellement un nombre de mailles non-acceptables élevé (il s’agit donc de configurations non optimales), mais que leur complémentarité fasse que le nombre de mailles non-acceptables en commun soit petit, voire minimal. En effet, il n’est en rien garanti que le nombre de mailles non-acceptables en commun soit minimal pour deux configurations individuellement optimales car elles peuvent ne pas être complémentaires. La septième étape 170 est ensuite une étape de réception d’images radiographiques 2D acquises pour chacune des première et deuxième configurations d’irradiation. Sont ainsi reçues une première pluralité d’images radiographiques 2D acquises et une deuxième pluralité de d’images radiographiques 2D acquises de la pièce 40 associées à la première configuration d’irradiation et à la deuxième configuration d’irradiation, respectivement. Une image radiographique 2D comprend une pluralité de pixels (ou de régions), chacun associé à un niveau de gris, indiquant une valeur d’amplitude, liée au flux de photons reçu L’image radiographique 2D représente le projeté de la pièce 40 sur l’écran 22 dans un plan coplanaire à la surface de l’écran 22. La première pluralité d’images radiographiques est acquise lorsque la pièce 40 est dans la première configuration d’irradiation, c’est-à-dire lorsque la pièce 40 se déplace selon le déplacement et depuis l’emplacement initial définis dans la première configuration d’irradiation. Chaque image de la première pluralité d’images radiographiques est associée à l’une des positions prédéterminées de la pluralité de positionnées prédéterminées successives de la première configuration d’irradiation. C’est-à-dire que pour chaque position prédéterminée de la pluralité de positionnées prédéterminées successives de la première configuration d’irradiation, au moins une image radiographique est acquise. La deuxième pluralité d’images radiographiques est acquise lorsque la pièce 40 est dans la deuxième configuration d’irradiation, c’est-à-dire lorsque la pièce 40 se déplace selon le déplacement et depuis l’emplacement initial définis dans la deuxième configuration d’irradiation. Chaque image de la deuxième pluralité d’images radiographiques est associée à l’une des positions prédéterminées de la pluralité de positionnées prédéterminées successives de la deuxième configuration d’irradiation. C’est-à-dire que pour chaque position prédéterminée de la pluralité de positionnées prédéterminées successives de la deuxième configuration d’irradiation, au moins une image radiographique est acquise. La huitième étape 180 est une étape optionnelle de recalage de la première pluralité d’images radiographiques 2D 40’ et de la deuxième d’images radiographiques 2D 40’. Le recalage est mis en œuvre au moyen d’un outil numérique de recalage, tel que par recalage de volumes, et utilisant, par exemple, une méthode des moindres carrés ou une distance de Chebyshev. Chaque image radiographique 40’ de chaque pluralité d’images radiographiques 40’ comprend une pluralité de portions 41’. Une portion 41’ d’une image radiographiques 40’ correspond, par exemple, à un ou plusieurs pixels de cette image. Le recalage est effectué par association d’une ou plusieurs portions 41’ de la pluralité de portions 41’ de chaque image radiographique 40’ de chaque pluralité d’images radiographiques 2D avec une maille 51 de recalage d’une pluralité de mailles 51 de recalage. La pluralité de mailles 51 de recalage est comprise dans la pluralité de mailles 51 du modèle numérique 3D 50. Une maille 51 de recalage est une maille 51 du modèle numérique 3D 50 servant au recalage de chaque image radiographique 40’ acquise avec le modèle numérique 50 et/ou avec d’autres images radiographiques 40’ acquises ou simulées. Le recalage associe alors chaque maille 51 de recalage du modèle numérique 3D 50 avec une portion 41’ correspondante dans chaque image radiographique 40’ de chacune des pluralités d’images radiographiques 40’. La portion 41’ correspondante comprend au moins une valeur d’amplitude, mesurée par l’appareil d’imagerie 20, associée à la partie de la pièce 40 qui correspond à ladite maille 51 de recalage. La portion 41’ correspondante correspond donc au projeté de la maille 51 de recalage dans l’image radiographique 40’. Les autres mailles 51 de la pluralité de mailles 51 du modèle numérique 3D 50, qui ne sont pas comprises dans la pluralité de mailles 51 de recalage, sont ensuite associées à une ou plusieurs portions 41’ de la pluralité de portions 41’ de chaque images de chaque pluralité d’images radiographiques 40’, typiquement suivant des méthodes connues en soi. Chaque portion 41’ comprend au moins une valeur d’amplitude mesurée par l’appareil d’imagerie 20 associée à la partie de la pièce 40 qui correspond à ladite maille 51 de la pluralité de mailles 51 du modèle de recalage 50. La neuvième étape 190 est une étape de construction de la représentation 3D de la pièce 40. La représentation 3D est construite à partir de la première pluralité d’images radiographiques 2D et de la deuxième d’images radiographiques 2D. La représentation 3D de la pièce comprend le modèle numérique 3D 50 auquel est associée une pluralité de valeurs d’amplitude, également appelées valeurs de maille. Chaque valeur de maille est déterminée à partir des valeurs d’amplitude provenant des images de la première pluralité d’images radiographiques 2D ou de la deuxième pluralité d’images radiographiques 2D. Chaque valeur de maille peut être représentée, dans la représentation 3D, par une couleur sur une échelle de couleur. Les nuances de couleur de l’échelle de couleur sont, par exemple, en niveau de gris et prennent des valeurs entre 0 et 255. Préférentiellement, à chaque maille 51 du modèle numérique 3D est associée une valeur de maille de la pluralité de valeurs d’amplitude. La valeur de maille est déterminée, pour chaque maille 51, à partir d’au moins une valeur d’amplitude d’au moins une image de la première pluralité d’images radiographiques 2D ou de la deuxième pluralité d’images radiographiques 2D. La valeur de maille peut être égale à une valeur d’amplitude directement mesurée par l’appareil d’imagerie 20, autrement appelée « valeur brute » ou être calculée à partir de cette valeur brute ou de plusieurs valeurs brutes, par exemple de valeurs brutes provenant d’images différentes de la première pluralité d’images radiographiques 2D et/ou de la deuxième pluralité d’images radiographiques 2D. La valeur de maille peut être déterminée à partir de valeurs brutes correspondant à des grandeurs physiques différentes. La ou leurs valeurs d’amplitude servant à la détermination de la valeur de maille sont extraites de la portion 41’ d’une ou plusieurs images radiographiques 40’ associée à ladite maille 51. La représentation 3D peut être déterminée par combinaison d’une ou plusieurs données issues de la première pluralité d’images radiographiques 40’ avec une plusieurs données de la deuxième pluralité d’images radiographiques 40’. La dixième étape 200 est ensuite une étape optionnelle de détection d’une non-conformité de la pièce 40. La non-conformité est détectée par comparaison de la représentation 3D de la pièce 40 avec une représentation 3D de référence. La non-conformité peut être détectée via n’importe quelle méthode de l’état de l’art. Typiquement, il peut s’agir d’une méthode qui évalue un écart de position des différentes surfaces de la pièce, ou d’une analyse géométrique plus élaborée et connue en soi. La représentation 3D de référence est, par exemple, une conception assistée par ordinateur de la pièce. La onzième étape 210 est enfin une étape optionnelle d’émission d’une alerte relative à la non-conformité détectée. L’alerte émise contient donc une information portant sur la non-conformité détectée. L’alerte est émise vers un opérateur, dans un cadre de production ou de maintenance, pour indiquer une anomalie de fabrication ou une anomalie engendrée par l’utilisation de la pièce 40. L’alerte est, par exemple, un message comprenant des caractères alphanumériques et émis vers l’opérateur via un module d’affichage. L’alerte peut être, en outre, mémorisée dans une mémoire. Dans un premier et un deuxième modes de réalisation, la neuvième étape 160 de détermination des première et deuxième configurations d’irradiations comprend deux sous-étapes. Dans le premier mode de réalisation, la première sous-étape est une étape de détermination d’un nombre de mailles 51 pour lesquelles la qualité de la représentation 3D est communément dégradée pour les configurations d’irradiation de chaque paire de configurations d’irradiation d’une pluralité de paires de configurations d’irradiation. Une pluralité de paires de configurations d’irradiation est donc formée à partir de l’ensemble de configurations d’irradiation. La pluralité de paires de configurations comprend toutes les paires possibles de configurations parmi l’ensemble de configurations. Alternativement, il est possible d’utiliser un algorithme de recherche itératif dont la mise en œuvre permet de rechercher de façon itérative une paire de configurations qui satisfait un critère d’optimisation. Le critère d’optimisation porte sur une minimisation du nombre de mailles pour lesquelles la valeur de l’indicateur de qualité calculé pour la première configuration et la valeur de l’indicateur de qualité calculé pour la deuxième configuration sont égales à la deuxième valeur. Dans ce cas, la recherche itérative est effectuée sur seulement une partie des paires possibles parmi les combinaisons de l’ensemble de configurations, et non sur toute les paires possibles. Ceci permet de déterminer plus rapidement la première et la deuxième configurations d’irradiation. Ceci permet également de déterminer plus finement la première et la deuxième configurations d’irradiation que dans le cas où la recherche est effectuée parmi toutes les configurations possibles, c’est-à-dire suivant une approche de type « force brute » (par exemple suivant un incrément fixe dans un intervalle)t Cet algorithme de recherche itérative permet, ainsi, de mettre en œuvre itérativement les étapes 120 à 160. L’algorithme est, par exemple, un algorithme de recherche locale, telle une descente de gradient, ou globale (algorithme génétique, essaim particulaire, algorithme des fourmis, etc.). Pour chaque paire de la pluralité de paires est calculé le nombre de mailles 51 pour lesquelles la qualité de la reconstruction 3D est dégradée par la présence d’un artefact à la fois par l’une et par l’autre des deux configurations d’irradiation de la paire de configurations d’irradiation. C’est-à-dire que, pour chaque paire de la pluralité de paires, est déterminé le nombre de mailles 51 dont la valeur de l’indicateur de qualité est égale à la deuxième valeur à la fois lorsque la pièce 40 est dans l’une des configurations d’irradiation de la paire de configurations concernée et lorsque la pièce 40 est dans l’autre configuration d’irradiation de la paire de configuration concernée. Autrement dit, pour chaque paire est calculé le nombre de mailles 51 qui sont communément non-acceptables pour les deux configurations d’irradiation de la paire concernée. Dans le premier mode de réalisation, la deuxième sous-étape est ensuite une étape de sélection de la paire de configurations d’irradiation, parmi la pluralité de paires de configurations d’irradiation, qui minimise le nombre de mailles 51 pour lesquelles la qualité de la reconstruction 3D est dégradée par la présence d’un artefact. Il s’agit donc de la paire de configurations minimisant le nombre de mailles 51 communément non-acceptables. La paire de configurations d’irradiation sélectionnée est celle pour laquelle le nombre de mailles 51 déterminé à la première sous-étape est le plus petit. La première configuration d’irradiation est alors l’une des configurations d’irradiation de la paire de configuration d’irradiation sélectionnée, et la deuxième configuration d’irradiation est l’autre configuration d’irradiation de la paire de configuration d’irradiation sélectionnée. Dans le premier mode de réalisation, pour chaque maille, la valeur de maille associée à ladite maille 51 est déterminée, à la neuvième étape 190 d’association, à partir de la pluralité d’images 40’ pour laquelle la qualité de la représentation 3D est la meilleure pour ladite maille 51. C’est-à-dire que la valeur d’amplitude servant à la détermination de la valeur de maille considérée provient de la pluralité d’images radiographiques 2D 40’ pour laquelle la configuration d’irradiation respective est la configuration d’irradiation pour laquelle le nombre de projections 50’ déterminé, pour ladite maille 51, qui dépasse le seuil prédéfini est le plus grand, parmi la première pluralité d’images radiographiques 2D 40’ et la deuxième pluralité d’images radiographiques 2D 40’. Lorsqu’aucune des première et deuxième configurations n’a son nombre de projections 50’ déterminé qui dépasse le seuil prédéfini, la valeur de maille est déterminée à partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d’au moins une image de la pluralité d’images 40’ associée à la première configuration d’irradiation et d’au moins un valeur d’amplitude provenant d’au moins une image de la pluralité d’images 40’ associée à la deuxième configuration d’irradiation. Par exemple, la valeur de maille est égale à une moyenne des valeurs d’amplitudes provenant des deux configuration d’irradiation. Alternativement, dans une telle situation, la valeur de maille est une valeur arbitrairement choisie ou une valeur issue d’une interpolation des mailles voisines de ladite maille. L’intérêt d’assigner une valeur à ladite maille lorsqu’aucune des deux configurations d’irradiation ne permet d’imager la pièce sans produire un artefact pour ladite maille est d’éviter d’engendrer des erreurs lors de l’analyse de la reproduction 3D, par des approches connues en soi. Le fait d’affecter une valeur issue d’une combinaison linéaire des deux acquisitions, plutôt qu’une valeur arbitraire ou provenant d’une interpolation inter-mailles, permet de conserver le réalisme de la représentation 3D. Dans le deuxième mode de réalisation, la première sous-étape est une étape de sélection de la première configuration d’irradiation parmi l’ensemble de configurations d’irradiation. La première configuration d’irradiation est sélectionnée de sorte à être la configuration pour laquelle la représentation 3D aura la meilleure qualité, c’est-à-dire celle minimisant le nombre d’artefact à la reconstruction, autrement dit celle pour laquelle il y a le plus de mailles 51 acceptables. La première configuration d’irradiation est donc la configuration d’irradiation de l’ensemble de configurations d’irradiation pour laquelle le nombre de mailles 51 pour lesquelles la valeur de l’indicateur de qualité calculé est égale à la deuxième valeur est le plus petit. Dans le deuxième mode de réalisation, la deuxième sous-étape est ensuite une étape de sélection de la deuxième configuration d’irradiation parmi l’ensemble de configurations d’irradiation, excepté la première configuration sélectionnée. La deuxième configuration d’irradiation est sélectionnée de sorte maximiser la qualité de la représentation 3D associée pour les mailles 51 pour lesquelles un artefact est présent dans la représentation 3D associée à la première configuration d’inspection. C’est-à-dire que la deuxième configuration permet de corriger la représentation 3D pour les mailles 51 pour lesquelles la représentation 3D est dégradée par la présence d’un artefact lorsque la pièce est dans la première configuration. Autrement dit, la deuxième configuration est celle pour laquelle le nombre de mailles 51 qui sont non-acceptables lorsque la pièce 40 est dans la première configuration et qui sont acceptables lorsque la pièce 40 est dans la deuxième configuration est le plus grand. La deuxième configuration d’irradiation est donc la configuration d’irradiation de l’ensemble de configurations d’irradiation pour laquelle un nombre de mailles 51 corrigées est le plus grand. Une maille 51 corrigée est une maille 51 du modèle numérique 3D pour laquelle : Un artefact est présent dans la représentation 3D lorsque la pièce 40 est dans la première configuration d’irradiation ; et Aucun artefact n’est présent dans la représentation 3D lorsque la pièce 40 est dans la deuxième configuration d’irradiation. Le nombre de mailles 51 corrigées est alors le nombre de mailles 51 pour lesquelles la valeur de l’indicateur de qualité calculé est égale à : La deuxième valeur lorsque la pièce 40 est dans la première configuration d’irradiation ; et La première valeur lorsque la pièce 40 est dans la deuxième configuration d’irradiation. Dans le deuxième mode de réalisation, pour chaque maille 51, la valeur de maille, à la neuvième étape 190 d’association, est : Déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image de la première pluralité d’images radiographiques 40’ si aucun artefact n’est généré pour ladite maille 51 lorsque la pièce 40 est dans la première configuration d’irradiation ; ou Déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image de la deuxième pluralité d’images radiographiques 40’ si aucun artefact n’est généré pour ladite maille 51 lorsque la pièce 40 est dans la deuxième configuration d’irradiation mais qu’un artefact est généré pour ladite maille 51 lorsque la pièce 40 est dans la première configuration d’irradiation ; Egale à une valeur arbitraire, est déterminée à partir d’une valeur issue d’une interpolation des mailles voisines de ladite maille ou déterminée à partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D 40’ et d’au moins une valeur d’amplitude provenant d'au moins une image de la deuxième pluralité d’images radiographiques 2D 40’ si le nombre de projections 50’ déterminé pour ladite maille 51 ne dépasse pas le seuil prédéfini lorsque la pièce 40 est dans la première configuration d’irradiation et si le nombre de projections 50’ déterminé pour ladite maille 51 ne dépasse pas le seuil prédéfini lorsque la pièce 40 est dans la deuxième configuration d’irradiation. En d’autres termes, pour chaque maille 51, la valeur de maille est : Déterminée à partir de la première pluralité d’images radiographiques 2D 40’ si le nombre de projections 50’ déterminé pour ladite maille 51 dépasse le seuil prédéfini lorsque la pièce 40 est dans la première configuration d’irradiation ; Déterminée à partir de la deuxième pluralité d’images radiographiques 2D 40’ si le nombre de projections 50’ déterminé pour ladite maille 51 ne dépasse pas le seuil prédéfini lorsque la pièce 40 est dans la première configuration d’irradiation et si le nombre de projections 50’ déterminé pour ladite maille 51 dépasse le seuil prédéfini lorsque la pièce 40 est dans la deuxième configuration d’irradiation ; Egale à une valeur arbitraire, est déterminée à partir d’une valeur issue d’une interpolation des mailles voisines de ladite maille ou est déterminée à partir d’une combinaison linéaire de valeurs d’amplitude extraites de la première et de la deuxième pluralité d’images radiographiques 2D 40’, lorsque le nombre de projections 50’ déterminé pour ladite maille 51 ne dépasse pas le seuil prédéfini alors que la pièce 40 est dans la première configuration d’irradiation ni alors que la pièce 40 est dans la première configuration d’irradiation. Dans un troisième mode de réalisation, compatible avec les modes de réalisation précédents, la première configuration d’irradiation et de la deuxième configuration d’irradiation sont déterminées, à la sixième étape 160, parmi un sous-ensemble de l’ensemble de configurations d’irradiation. Le sous-ensemble de configurations d’irradiation comprend une ou plusieurs configurations d’irradiation pour lesquelles aucun artefact n’est généré lors de la construction de la représentation 3D pour une ou plusieurs des mailles 51 de recalage de la pluralité de mailles 51 de recalage. En d’autres termes, la première et la deuxième configurations d’irradiation sont choisies de sorte que les mailles de recalage ne comportent pas d’artefacts à la reconstruction selon au moins l’une des première ou deuxième configurations d’irradiation. Méthode (100), mise en œuvre par ordinateur, de détermination d’une représentation 3D d’une pièce (40) par tomographie à faisceau conique émis depuis une source (21), la méthode (100) comprenant : Recevoir (110) un modèle numérique 3D de la pièce (40) comprenant une pluralité de mailles (51) ; Obtenir (120) un ensemble de configurations d’irradiation de la pièce (40), chaque configuration d’irradiation étant associée une pluralité respective de positions prédéterminées successives de la pièce (40) par rapport à la source (21) ; Pour chaque configuration d’irradiation de l’ensemble de configurations d’irradiation : Déterminer (130) une pluralité de projections (50’) du modèle numérique de la pièce (40), chaque projection (50’) étant respectivement associée à une position prédéterminée de la pluralité de positions prédéterminées successives associée à ladite configuration d’irradiation, chaque projection (50’) comprenant un ensemble de portions (41’), chaque portion (41’) étant associée à une maille (51) respective du modèle numérique 3D ; Pour chaque maille (51) du modèle numérique 3D : Déterminer (140) un nombre de projections (50’) parmi la pluralité de projections (50’) déterminée pour lesquelles la portion (41’) associée à la maille (51) respecte un critère de qualité prédéfini, le critère de qualité étant lié à un risque de présence d’un artefact dans une représentation 3D de la pièce (40) construite à partir des projections (50’) déterminées, le critère de qualité étant une condition de Tuy ; Si le nombre de projections (50’) déterminé dépasse un seuil prédéfini, associer (150) à la maille (51) une première valeur d’un indicateur de qualité ; sinon, associer (150) à la maille (51) une deuxième valeur de l’indicateur de qualité ; Déterminer (160) une première configuration d’irradiation et une deuxième configuration d’irradiation parmi l’ensemble de configurations d’irradiation en minimisant un nombre de mailles (51) pour lesquelles la valeur de l’indicateur de qualité calculé pour la première configuration et la valeur de l’indicateur de qualité calculé pour la deuxième configuration sont égales à la deuxième valeur ; Recevoir (170) une première pluralité d’images radiographiques 2D (40’) de la pièce (40) lorsque la pièce (40) est dans la première configuration d’irradiation et une deuxième pluralité de d’images radiographiques 2D (40’) de la pièce (40) lorsque la pièce (40) est dans la deuxième configuration d’irradiation ; Construire (190) la représentation 3D de la pièce (40) à partir de la première pluralité d’images radiographiques 2D (40’) et de la deuxième pluralité d’images radiographiques 2D (40’). Méthode (100) selon la revendication 1, dans laquelle la représentation 3D de la pièce (40) comprend le modèle numérique 3D et une pluralité de valeurs d’amplitude, dites valeurs de mailles, chaque valeur de maille de la pluralité de valeurs de mailles étant respectivement associée à une maille (51) de la pluralité de mailles (51), chaque valeur de maille étant déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image radiographique 2D (40’) de la première pluralité d’images radiographiques 2D (40’) et/ou de la deuxième pluralité d’images radiographiques 2D (40’). Méthode (100) selon l’une des revendications précédentes, comprenant en outre : Détection (200) d’une non-conformité de la pièce (40) par comparaison de la représentation 3D de la pièce (40) avec une représentation 3D de référence ; Emission (210) d’une alerte relative à la non-conformité détectée. Méthode (100) selon l’une des revendications précédentes, dans laquelle la détermination de la première configuration d’irradiation et de la deuxième configuration d’irradiation comprend : Déterminer, pour chaque paire de configurations d’irradiation d’une pluralité de paires de configurations d’irradiation, le nombre de mailles (51) pour lesquelles la valeur de l’indicateur de qualité calculé de chacune des configurations d’irradiation de ladite paire de configurations d’irradiation est égale à la deuxième valeur ; Sélectionner, parmi la pluralité de paires de configurations d’irradiation, la paire de configurations d’irradiation pour laquelle le nombre de mailles (51) est minimal, la première configuration d’irradiation étant l’une des configurations d’irradiation de la paire de configuration d’irradiation sélectionnée, et la deuxième configuration d’irradiation étant l’autre configuration d’irradiation de la paire de configuration d’irradiation sélectionnée. Méthode (100) selon la revendication 4 et la revendication 2, dans laquelle, pour chaque maille (51) de la pluralité de mailles (51), la détermination de la valeur de maille comprend : Lorsqu’au moins un nombre parmi le nombre de projections (50’) déterminé pour la première configuration d’irradiation et le nombre de projections (50’) déterminé pour la deuxième configuration dépasse le seuil prédéfini, sélectionner, parmi la première configuration d’irradiation et la deuxième configuration d’irradiation, la configuration d’irradiation pour laquelle le nombre de projections (50’) déterminé dépassant le seuil prédéfini pour ladite maille (51) est le plus grand, la valeur de maille étant déterminée à partir d’au moins une valeur d’amplitude provenant d’au moins une image radiographique 2D (40’) de la pluralité d’images radiographiques 2D (40’) associée à la configuration d’irradiation sélectionnée ; Lorsque le nombre de projections (50’) déterminé pour la première configuration d’irradiation et le nombre de projections (50’) déterminé pour la deuxième configuration ne dépassent pas le seuil prédéfini, la valeur de maille est déterminée à partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D (40’) et d’au moins une valeur d’amplitude provenant d’au moins une image de la deuxième pluralité d’images radiographiques 2D (40’). Méthode (100) selon l’une des revendications 1 à 3, dans laquelle la détermination de la première configuration d’irradiation et de la deuxième configuration d’irradiation comprend : Sélectionner comme première configuration d’irradiation la configuration d’irradiation ayant un plus petit nombre de mailles (51) pour lesquelles la valeur de l’indicateur de qualité calculé est égale à la deuxième valeur parmi l’ensemble de configurations d’irradiation ; Sélectionner comme deuxième configuration d’irradiation la configuration d’irradiation ayant un plus grand nombre de mailles (51) corrigées parmi l’ensemble de configurations d’irradiation, le nombre de mailles (51) corrigées étant le nombre de mailles (51) pour lesquelles : La valeur de l’indicateur de qualité calculé pour la première configuration d’irradiation est égale à la deuxième valeur ; et La valeur de l’indicateur de qualité calculé pour la deuxième configuration d’irradiation est égale à la première valeur. Méthode (100) selon la revendication 6 et la revendication 3, dans laquelle la valeur de maille est déterminée : A partir d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D (40’) si le nombre de projections (50’) déterminé pour ladite maille (51) dépasse le seuil prédéfini lorsque la pièce (40) est dans la première configuration d’irradiation ; A partir d’au moins une valeur d’amplitude provenant d'au moins une image de la deuxième pluralité d’images radiographiques 2D (40’) si le nombre de projections (50’) déterminé pour ladite maille (51) ne dépasse pas le seuil prédéfini lorsque la pièce (40) est dans la première configuration d’irradiation et si le nombre de projections (50’) déterminé pour ladite maille (51) dépasse le seuil prédéfini lorsque la pièce (40) est dans la deuxième configuration d’irradiation ; A partir d’une combinaison linéaire d’au moins une valeur d’amplitude provenant d'au moins une image de la première pluralité d’images radiographiques 2D (40 ‘) et d’au moins une valeur d’amplitude provenant d'au moins une image de la deuxième pluralité d’images radiographiques 2D (40’) si le nombre de projections (50’) déterminé pour ladite maille (51) ne dépasse pas le seuil prédéfini lorsque la pièce (40) est dans la première configuration d’irradiation et si le nombre de projections (50’) déterminé pour ladite maille (51) ne dépasse pas le seuil prédéfini lorsque la pièce (40) est dans la deuxième configuration d’irradiation. Méthode (100) selon l’une des revendications précédentes, comprenant en outre : recaler (180) la première pluralité d’images radiographiques 2D (40’) et la deuxième pluralité d’images radiographiques 2D (40’) avec le modèle numérique 3D de la pièce (40) à partir d’une pluralité de mailles (51) de recalage, la construction de la représentation 3D étant effectuée à partir des images recalées. Méthode (100) selon la revendication 8, dans laquelle la première configuration d’irradiation et de la deuxième configuration d’irradiation sont déterminés parmi un sous-ensemble de configurations d’irradiation, le sous-ensemble de configurations d’irradiation comprenant une ou plusieurs configurations d’irradiation de l’ensemble de configurations d’irradiations pour lesquelles la valeur de l’indicateur de qualité calculé est égale à la première valeur pour une ou plusieurs mailles (51) de recalage de la pluralité de mailles (51) de recalage. Système de détermination d’une représentation 3D d’une pièce (40) par tomographie à faisceau conique émis depuis une source (21), le système étant configuré pour mettre en œuvre la méthode (100) selon l’une des revendications 1 à 9.