Brevet FR3077138A1 — publié le 20190726

PROCEDE DE CONTROLE NON DESTUCTIF POUR UNE PIECE AERONAUTIQUE DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne le domaine du contrôle non destructif, plus particulièrement appliqué à une pièce aéronautique, par exemple une pièce composite. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Parmi les méthodes de contrôle non destructif, on distingue celles utilisant la tomographie à rayons X. L'imagerie par rayons X permet d'accéder à la microstructure interne d'une pièce aéronautique. La tomographie à rayons X de la pièce aéronautique a pour résultat une image 3D. L'analyse de cette image 3D permet notamment le contrôle non destructif de la pièce aéronautique. La tomographie d'une pièce nécessite tout d'abord l'acquisition d'un grand nombre de radiographies de la pièce, typiquement un ou plusieurs milliers. Elle nécessite ensuite le calcul, appelé reconstruction, de l'image 3D de la pièce à partir des radiographies. Ces deux étapes sont couteuses en temps, aussi bien pour l'étape d'acquisition qui peut demander une ou plusieurs heures que pour l'étape de calcul. EXPOSÉ DE L'INVENTION L'invention vise à résoudre les problèmes de la technique antérieure en fournissant un procédé de contrôle non destructif pour une pièce aéronautique mettant en oeuvre un dispositif de tomographie comportant un générateur de rayons X et un détecteur de rayons X, le procédé comportant l'étape de : - Détermination d'un nombre de projections de la pièce aéronautique à effectuer, chaque projection correspondant à un angle respectif d'un 2 premier ensemble d'angles échantillonnés de manière uniforme, et à un premier pas d'échantillonnage uniforme de détecteur, caractérisé en ce qu'il comporte en outre des étapes de : - Détermination de facteurs de forme de la pièce aéronautique, en 5 fonction d'une quantité d'informations utiles de la pièce aéronautique, respectivement selon les axes d'un repère tridimensionnel, - Identification de projections à effectuer, en fonction du nombre de projections à effectuer et des facteurs de forme de la pièce aéronautique, et détermination d'un second ensemble d'angles échantillonnés de manière non uniforme 10 correspondant respectivement aux projections à effectuer, - Acquisition des projections identifiées, selon le second ensemble d'angles et le premier pas d'échantillonnage de détecteur, - Ré-échantillonnage des projections acquises, par transformation des projections selon des facteurs respectifs dépendant du nombre de projections et des 15 facteurs de forme de la pièce aéronautique, puis ajustement des projections selon un second pas respectif d'échantillonnage de détecteur dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, et - Reconstruction d'une image 3D de la pièce aéronautique à partir des projections ré-échantillonnées. 20 L'invention repose sur l'exploitation du caractère fortement anisotrope de la microstructure des matériaux. A titre d'exemple, on peut citer les matériaux composites à renfort 3D. Plus particulièrement, la propriété anisotrope d'intérêt sera la microstructure mesurée le long des différents axes. 25 Grâce à l'invention, l'acquisition d'informations non pertinentes est réduite, et par conséquent le stockage et le traitement de ces informations non pertinentes également. L'invention permet ainsi de réduire significativement les temps d'acquisition et de calcul pour l'obtention d'une image 3D de qualité comparable à ce qui est pratiqué selon la technique antérieure. 3 Les facteurs de forme sont choisis de sorte que l'image 3D de la pièce aéronautique soit isotrope. Selon une caractéristique préférée, le procédé de contrôle non 5 destructif pour une pièce aéronautique comporte en outre une étape de transformation de l'image 3D de la pièce aéronautique en une autre image 3D de la pièce aéronautique dans le repère tridimensionnel. Il est ainsi possible de revenir à une représentation de la pièce aéronautique qui soit anisotrope. 10 Selon une caractéristique préférée, la détermination de facteurs de forme de la pièce aéronautique est telle qu'un des facteurs de forme selon un axe du repère tridimensionnel, correspondant à l'axe de rotation de l'objet dans le dispositif de tomographie, est égal à 1. 15 Selon une caractéristique préférée, la détermination d'un second ensemble d'angles correspondant respectivement aux projections à effectuer comporte pour chaque indice i compris entre 1 et N, où N est le nombre de projections à effectuer, Ol=atan la détermination du ième second angle cl); par : -./..), _ --.- 20 Où i et J représentent respectivement les vecteurs unitaires des axes du repère tridimensionnel, ne correspondant pas à l'axe de rotation de l'objet dans le dispositif de tomographie, Où est un vecteur unitaire directeur défini par : A= Où A est une matrice diagonale définie par : 0 aY), dans 25 laquelle ax et ay sont les facteurs de forme selon les axes du repère tridimensionnel ne 4 correspondant pas à l'axe de rotation de l'objet dans le dispositif de tomographie et p., est un vecteur unitaire défini par : 360 a. = Où a, est un angle du premier ensemble d'angles défini par : 5 Selon une caractéristique préférée, la transformation des projections selon des facteurs respectifs dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique comporte pour chaque indice i compris entre 1 et N, où N est le nombre de projections à effectuer, la transformation de la ième projection acquise c = I Ien utilisant un facteur c défini par : det(A)1 10 Selon une caractéristique préférée, l'ajustement de la ième projection acquise selon un second pas d'échantillonnage respectif de détecteur comporte la détermination du second pas t, défini par : Où p est le pas d'échantillonnage uniforme du détecteur du dispositif de 15 tomographie à rayons X. L'invention concerne aussi un dispositif de contrôle non destructif pour une pièce aéronautique, comportant un générateur de rayons X, un détecteur de rayons X et un module de calcul apte à déterminer un nombre de projections de la pièce 20 aéronautique à effectuer, chaque projection correspondant à un angle respectif d'un premier ensemble d'angles échantillonnés de manière uniforme, et à un premier pas d'échantillonnage uniforme de détecteur, caractérisé en ce que le module de calcul est en outre apte à : - Déterminer des facteurs de forme de la pièce aéronautique, en 25 fonction d'une quantité d'informations utiles de la pièce aéronautique, respectivement selon les axes d'un repère tridimensionnel, 5 - Identifier des projections à effectuer, en fonction du nombre de projections à effectuer et des facteurs de forme de la pièce aéronautique, et déterminer un second ensemble d'angles échantillonnés de manière non uniforme correspondant respectivement aux projections à effectuer, 5 - Acquérir les projections identifiées, selon le second ensemble d'angles et le premier pas d'échantillonnage uniforme de détecteur, - Ré-échantillonner les projections acquises, par transformation des projections selon des facteurs respectifs dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, puis ajuster les projections selon un second 10 pas respectif d'échantillonnage de détecteur dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, et - Reconstruire une image 3D de la pièce aéronautique à partir des projections ré-échantillonnées. 15 Selon une caractéristique préférée, le module de calcul est en outre apte à transformer de l'image 3D de la pièce aéronautique en une autre image 3D de la pièce aéronautique dans le repère tridimensionnel. Le dispositif présente des avantages analogues à ceux précédemment présentés. 20 Dans un mode particulier de réalisation, les étapes du procédé selon l'invention sont mises en oeuvre par des instructions de programme d'ordinateur. En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un 25 ordinateur, ce programme comportant des instructions adaptées à la mise en oeuvre des étapes d'un procédé tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe 30 quelle autre forme souhaitable. 6 L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions de programme d'ordinateur adaptées à la mise en oeuvre des étapes d'un procédé tel que décrit ci-dessus. Le support d'informations peut être n'importe quelle entité ou dispositif 5 capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette ou un disque dur. D'autre part, le support d'informations peut être un support 10 transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour 15 être utilisé dans l'exécution du procédé selon l'invention. BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages apparaîtront à la lecture de la description suivante d'un mode de réalisation préféré donné à titre d'exemple non 20 limitatif, décrit en référence aux figures dans lesquelles : La figure 1 représente un dispositif de contrôle non destructif pour une pièce aéronautique selon un mode de réalisation de l'invention, La figure 2 représente un procédé de contrôle non destructif pour une pièce aéronautique selon un mode de réalisation de l'invention, 25 La figure 3 représente un exemple d'objet dans sa configuration réelle, avant traitement selon la présente invention, La figure 4 représente l'objet de la figure 3, dans une représentation isotrope, après transformation selon la présente invention. 7 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Selon un mode de réalisation préféré représenté à la figure 1, un dispositif de contrôle non destructif pour une pièce aéronautique comporte un générateur à rayons X 1 apte à délivrer un faisceau de rayons X. Dans la suite, on 5 s'intéresse plus particulièrement à l'émission de rayons X en faisceaux parallèles, mais l'émission peut également être en faisceaux divergents collimatés ou encore en faisceaux coniques. Un détecteur de rayons X 2 est placé en vis-à-vis du générateur 1, de manière à recevoir les rayons X émis par le générateur 1. 10 Un plateau 3 est disposé entre le générateur 1 et le détecteur 2. Le plateau 3 est rotatif selon un axe vertical et est destiné à recevoir un objet 4, tel qu'une pièce aéronautique. On définit un repère tridimensionnel (Ox, Oy, Oz), où l'axe Oz est l'axe de rotation du plateau 3. La taille du plateau 3 dépend de la taille des objets à traiter. Lorsque l'objet 4 est disposé sur le plateau 3, les rayons X émis par le générateur 1 15 traversent l'objet 4 avant d'atteindre le détecteur 2. Un ordinateur 5 est relié au générateur 1, au détecteur 2 et au plateau 3. L'ordinateur 5 permet à un opérateur de commander le fonctionnement de ces différents éléments, de collecter les résultats de la détection de rayons X par le détecteur 2 et de traiter les données. 20 L'ordinateur 5 comporte notamment un processeur 100, une mémoire 101, une interface d'entrée 102 et une interface de sortie 103. Ces différents éléments sont classiquement reliés par un bus 105. Le processeur 100 exécute un programme d'ordinateur mettant en oeuvre le procédé selon l'invention. Ces traitements sont réalisés sous la forme 25 d'instructions de code du programme d'ordinateur qui sont mémorisées par la mémoire 101 avant d'être exécutées par le processeur 100. L'interface de sortie 103 délivre les données qui représentent des instructions de commandes du générateur 1, du détecteur 2 et du plateau 3. 8 L'interface d'entrée 102 est reliée au détecteur 2 et est destinée à recevoir les données représentant les projections effectuées comme exposé dans la suite. Par « projection », on désigne une image h qui est calculée à partir d'une radiographie p. Il ne s'agit pas du signal brut du détecteur. Pour ce calcul, il est 5 nécessaire de connaitre l'intensité de la source seule (sans échantillon) en chaque point du détecteur. Cette radiographie supplémentaire dont on fait l'acquisition en l'absence de pièce, avant ou après la tomographie de celle-ci, est nommée « image de blancs » (en anglais : flat-field) po. Dans les applications les plus courantes, la projection se calcule selon la formule : 10 h=log(po/p). D'autres variantes existent, comme la prise en compte d'une radiographie acquise avec la source. Cette radiographie, éteinte qui qualifie le bruit d'obscurité du détecteur, est nommée « image de noir » (en anglais : dark-field) pi. Dans ce cas, la projection se calcule selon la formule : 15 h=108(Po-P1)/(P-131). Il existe de nombreuses sophistications de cette formule bien connues par l'homme de métier, tout cela étant très lié au savoir-faire des fournisseurs de tomographe. 20 Le procédé de contrôle non destructif de la figure 2 représente le fonctionnement du dispositif de contrôle non destructif de la figure 1. Le procédé comporte des étapes El à E7. On suppose que l'objet 4 a été placé sur le plateau 3 dans le but d'être tomographié. 25 L'étape El est la détermination de la valeur de facteurs de forme ax, ay et a,. Les facteurs de forme représentent la quantité d'informations essentielles de l'objet selon chaque axe du repère (0, x, y, z). Les facteurs de forme dépendent de l'objet et de l'utilisation ultérieure des résultats de la tomographie. Les facteurs de forme permettent de passer d'un objet réel dont la 30 structure est anisotrope à une représentation de cet objet, redimensionnée par 9 transformation géométrique. La transformation est choisie pour produire une représentation transformée de l'objet qui est isotrope. La transformation choisie est une affinité orthogonale d'axe Oz qui est comme on l'a vu l'axe de rotation du plateau. Le facteur de forme a, selon cet axe est égal 5 à 1, ce qui signifie qu'il n'y a pas de transformation selon cet axe. La transformation est une dilatation de facteur de forme, ou d'échelle, ax, respectivement ay, le long de l'axe Ox, respectivement Oy. Les facteurs de forme ax et ay ne sont pas nécessairement déterminés avec précision. De simples approximations peuvent suffire. 10 La transformation géométrique permet donc de passer d'un objet anisotrope à une représentation isotrope de cet objet contenant autant d'informations utiles que l'objet. La figure 3 représente un exemple d'objet dans sa configuration réelle. 15 Le matériau qui constitue l'objet est un composite tissé qui comporte des fibres dont la section transversale est allongée dans les directions des axes Ox et Oy, en raison du procédé de fabrication de l'objet. Ce matériau est anisotrope. La figure 4 représente le même objet, après transformation géométrique pour en fournir une représentation isotrope. Les facteurs de forme ax et ay 20 ont été choisis pour que la section transversale des fibres soit sensiblement circulaire. En référence à nouveau à la figure 2, le résultat de l'étape El est une matrice diagonale A3D de transformation linéaire : (a' 0 0 /43D = 0 a 0 Y 0 0 az/ où a, = 1. 25 Cette matrice relie le repère tridimensionnel (Ox, 0y, Oz) de l'objet réel, donc anisotrope, avec le repère tridimensionnel (Ou, Ov, Ow) de la représentation isotrope de l'objet : (u,v,w)T= A3D(xS,z)T. 10 Ces repères sont représentés dans les figures 3 et 4. Puisque le repère (Ou, Ov, Ow) désigne le repère dans lequel la représentation de l'objet est isotrope, dans la suite il sera simplement nommé repère isotrope. Similairement le repère (Ox, Oy, Oz), sera nommé repère anisotrope. 5 Afin de simplifier les équations suivantes, la matrice A3D est simplifiée en : A = L'étape suivante E2 est la détermination d'un nombre N de projections à 10 effectuer, chaque projection correspondant à un angle respectif d'un premier ensemble d'angles échantillonnés de manière uniforme dans le repère isotrope. On considère ici un pas p d'échantillonnage de détecteur constant, par exemple égal à 1 pixel. Le pas p d'échantillonnage de détecteur est également uniforme. Le nombre de projections à effectuer est un compromis entre qualité 15 d'image et limitation du temps d'acquisition, entre autres. Par exemple, le nombre de projections à effectuer est déterminé en fixant comme règle que la projection du point de l'objet le plus éloigné de l'axe de rotation dans le repère isotrope ne se déplace pas de plus d'un pixel sur le détecteur 2 entre deux projections successives. Le nombre total de projections est réparti de manière équiangulaire sur une rotation de D degrés de l'objet, 20 ce qui produit un incrément d'angle constant entre deux projections successives. La valeur de D est habituellement 360°, mais dans le cas d'un faisceau parallèle de rayons X, la rotation peut être de 180°. Le nombre N de projections à effectuer est donc associé à un ensemble d'angles échantillonnés de manière uniforme a, et exprimés en degrés, et à leurs vecteurs 25 unitaires directeurs respectifs = , avec i =[1,N] 11 Le résultat de l'étape E2 est donc le nombre de projections à effectuer N et l'ensemble des angles a, associés à leurs vecteurs unitaires directeurs respectifs p.,. L'étape suivante E3 est l'identification de projections à effectuer dans le 5 repère anisotrope. Cette étape comporte la détermination de positions angulaires auxquelles doivent être effectuées les projections dans le repère anisotrope, à partir de celles dans le repère isotrope. Cette détermination dépend du nombre N de projections à effectuer et des facteurs de forme de l'objet. Un ensemble de vecteurs unitaires directeurs pour i compris entre 1 10 et N, contenant les directions des projections est déterminé. L'ensemble des vecteurs unitaires directeurs est défini par : - , avec i = [1,N] Un second ensemble d'angles cl), correspondant respectivement aux 15 projections à effectuer est déterminé. Les angles cl), sont échantillonnés de manière non uniforme. L'ensemble des angles cl), est défini par : ( = a tan , avec i = [1,N] où i et j représentent respectivement les vecteurs unitaires des axes 20 Ox et Oy. Le résultat de l'étape E3 est donc l'ensemble des N projections à effectuer, définies par l'ensemble des angles cl), respectifs, pour i allant de 1 à N. 25 L'étape suivante E4 est l'acquisition des N projections identifiées à l'étape E3. Pour chacune de ces projections, l'acquisition comporte l'émission de rayons X par le générateur 1. Les rayons X atteignent le détecteur 2 en traversant l'objet 4. Pour chacune des N projections, la position de l'objet est définie par l'angle cl),, i allant de 1 à N. 12 Le résultat de l'étape E4 est l'ensemble des projections acquises g-, pour chacune des positions angulaires définies par les angles c, pour i allant de 1 à N, ce qui forme un échantillonnage angulaire non uniforme, et pour le pas d'échantillonnage de détecteur uniforme p, par exemple égal à 1 pixel. 5 L'étape suivante E5 est le ré-échantillonnage des projections acquises à l'étape précédentes. Il s'agit tout d'abord d'effectuer une transformation de chaque projection acquise en utilisant un facteur respectif c, dépendant de la matrice A 10 déterminée à l'étape El et du nombre de projections déterminé à l'étape E2. Le facteur c, est défini par : 1,441,1 C1 - det(A)1, pour i allant de 1 à N. I Il faut ensuite effectuer un ajustement de chacune des N projections 15 acquises selon un pas d'échantillonnage respectif t, de détecteur, pour i allant de 1 à N, dépendant de la matrice A déterminée à l'étape El et du nombre de projections déterminé à l'étape E2. Cette opération doit prendre en compte le fait que les projections représentent une densité à ajuster selon l'orientation. Le pas d'échantillonnage t, de détecteur est défini par : 20 t1 A./.111.p , pour i allant de là N. Le pas d'échantillonnage t, de détecteur est uniforme pour chaque orientation considérée, mais différent d'une orientation à l'autre, ce qui forme un pas d'échantillonnage t, de détecteur non uniforme. Le résultat de l'étape E5 est l'ensemble des N projections ré- 25 échantillonnées f correspondant respectivement aux positions angulaires définies par les angles c, pour i allant de 1 à N, ce qui forme un échantillonnage angulaire non uniforme, et pour un pas d'échantillonnage non uniforme de détecteur t,, pour i allant de là N: 1 kr,(P,) cl où t, u et sont des vecteurs qui regroupent respectivement les valeurs de t,, p., et pour i allant de 1 à N. L'étape suivante E6 est la reconstruction d'une image 3D de l'objet 4 à 5 partir des projections ré-échantillonnées obtenues à l'étape précédente. La reconstruction est classique en soi. La reconstruction repose soit sur l'utilisation de la transformée de Fourier soit des techniques dites algébriques. Le résultat de l'étape E6 est une image 3D f(u, v, w) de l'objet dans le repère isotrope (0, u, v, w). Cette image est une représentation isotrope de l'objet. 10 L'étape suivante E7 est optionnelle. L'étape E7 est la transformation de l'image 3D de l'objet dans le repère isotrope (0, u, v, w) obtenue à l'étape précédente en une image 3D de l'objet dans le repère anisotrope (0, x, y, z). Cette étape comporte l'application à l'image f(u, y, w) d'une transformation pour obtenir une image g(x, y, z). Cette transformation est définie par : 15 g(x,y,z) 1 f( u y 142 Idet(A)1 a ' REVENDICATIONS1. Procédé de contrôle non destructif pour une pièce aéronautique mettant en oeuvre un dispositif de tomographie comportant un générateur de rayons X et un détecteur de rayons X, le procédé comportant l'étape de : - Détermination (E2) d'un nombre (N) de projections de la pièce aéronautique à effectuer, chaque projection correspondant à un angle respectif d'un premier ensemble d'angles échantillonnés de manière uniforme, et à un premier pas d'échantillonnage uniforme de détecteur, caractérisé en ce qu'il comporte en outre des étapes de : - Détermination (El) de facteurs de forme de la pièce aéronautique, en fonction d'une quantité d'informations utiles de la pièce aéronautique, respectivement selon les axes d'un repère tridimensionnel, - Identification (E3) de projections à effectuer, en fonction du nombre de projections à effectuer et des facteurs de forme de la pièce aéronautique, et détermination d'un second ensemble d'angles échantillonnés de manière non uniforme correspondant respectivement aux projections à effectuer, - Acquisition (E4) des projections identifiées, selon le second ensemble d'angles et le premier pas d'échantillonnage de détecteur, - Ré-échantillonnage (E5) des projections acquises, par transformation des projections selon des facteurs respectifs dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, puis ajustement des projections selon un second pas respectif d'échantillonnage de détecteur dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, et 15 - Reconstruction (E6) d'une image 3D de la pièce aéronautique à partir des projections ré-échantillonnées. 2. Procédé de contrôle non destructif pour une pièce aéronautique 5 selon la revendication 1, caractérisé en ce qu'il comporte en outre une étape de transformation (E7) de l'image 3D de la pièce aéronautique en une autre image 3D de la pièce aéronautique dans le repère tridimensionnel. 3. Procédé de contrôle non destructif pour une pièce aéronautique 10 selon la revendication 1 ou 2, caractérisé en ce que la détermination (El) de facteurs de forme de la pièce aéronautique est telle qu'un des facteurs de forme selon un axe du repère tridimensionnel, correspondant à l'axe de rotation de l'objet dans le dispositif de tomographie, est égal à 1. 15 4. Procédé de contrôle non destructif pour une pièce aéronautique selon la revendication 3, caractérisé en ce que la détermination (E3) d'un second ensemble d'angles correspondant respectivement aux projections à effectuer comporte pour chaque indice i compris entre 1 et N, où N est le nombre de projections à effectuer, la détermination du ième second angle cl); par : = a tan 20 Où i et j représentent respectivement les vecteurs unitaires des axes du repère tridimensionnel, ne correspondant pas à l'axe de rotation de l'objet dans le dispositif de tomographie, Où est un vecteur unitaire directeur défini par : A.p Où A est une matrice diagonale définie par : A = , dans 25 laquelle ax et ay sont les facteurs de forme selon les axes du repère tridimensionnel ne 16 correspondant pas à l'axe de rotation de l'objet dans le dispositif de tomographie et p., est un vecteur unitaire défini par : pi = Où a, est un angle du premier ensemble d'angles défini par : 5 5. Procédé de contrôle non destructif pour une pièce aéronautique selon la revendication 4, caractérisé en ce que la transformation (E5) des projections selon des facteurs respectifs dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique comporte pour chaque indice i compris entre 1 et N, où N est le nombre de projections à effectuer, la transformation de la ième projection acquise 1A.p, 10 en utilisant un facteur c, défini par : c- Idet(A)1" 6. Procédé de contrôle non destructif pour une pièce aéronautique selon la revendication 4 ou 5, caractérisé en ce que l'ajustement (E5) de la ième projection acquise selon un second pas d'échantillonnage respectif de détecteur comporte la 15 détermination du second pas t, défini par : t1=1A1.111.p Où p est le pas d'échantillonnage uniforme du détecteur du dispositif de tomographie à rayons X. 7. Dispositif de contrôle non destructif pour une pièce aéronautique, 20 comportant un générateur de rayons X (1), un détecteur de rayons X (2) et un module de calcul (5) apte à déterminer un nombre de projections de la pièce aéronautique à effectuer, chaque projection correspondant à un angle respectif d'un premier ensemble d'angles échantillonnés de manière uniforme, et à un premier pas d'échantillonnage uniforme de détecteur, 25 caractérisé en ce que le module de calcul est en outre apte à : 17 - Déterminer des facteurs de forme de la pièce aéronautique, en fonction d'une quantité d'informations utiles de la pièce aéronautique, respectivement selon les axes d'un repère tridimensionnel, - Identifier des projections à effectuer, en fonction du nombre de 5 projections à effectuer et des facteurs de forme de la pièce aéronautique, et déterminer un second ensemble d'angles échantillonnés de manière non uniforme correspondant respectivement aux projections à effectuer, - Acquérir les projections identifiées, selon le second ensemble d'angles et le premier pas d'échantillonnage uniforme de détecteur, 10 - Ré-échantillonner les projections acquises, par transformation des projections selon des facteurs respectifs dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, puis ajuster les projections selon un second pas respectif d'échantillonnage de détecteur dépendant du nombre de projections et des facteurs de forme de la pièce aéronautique, et 15 - Reconstruire une image 3D de la pièce aéronautique à partir des projections ré-échantillonnées. 8. Dispositif de contrôle non destructif pour une pièce aéronautique selon la revendication 1, caractérisé en ce que le module de calcul est en outre apte à 20 transformer l'image 3D de la pièce aéronautique en une autre image 3D de la pièce aéronautique dans le repère tridimensionnel. 9. Programme d'ordinateur comportant des instructions pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 6 25 lorsque ledit programme est exécuté par un ordinateur. 10. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution des étapes du procédé selon l'une quelconque des revendications 1 à 6.