FR 3047339 A1 20170804 FR 1650787 A 20160201 1 Procédé de contrôle non-destructif par redressement DOMAINE TECHNIQUE GENERAL L'invention concerne le domaine du contrôle non destructif (CND) sur des pièces industrielles, notamment dans le domaine de l'aéronautique, à l'aide de volumes tomographiques numériques. Plus précisément, l'invention concerne l'amélioration du contrôle visuel effectué par un opérateur. Le CND est primordial pour contrôler la santé matière des matériaux. Par exemple, les aubes de soufflante, qui sont en composite de fibres de carbone tissées en trois dimensions, sont des pièces critiques qui doivent être contrôlées en intégralité. Mais le carter fan, le redresseur, les pales, les aubes, etc. peuvent aussi être concernés. Les volumes tomographiques numériques sont obtenues à l'aide d'un tomographe dont le générateur de rayons X émet un faisceau traversant la pièce à explorer, avant d'être analysé, après atténuation, par un système de détection. L'image intermédiaire ainsi obtenue est appelée une « projection ». En acquérant une pluralité de projections dans différents plans de l'espace (avec un éventuel prétraitement) et en les recombinant, on obtient un volume tridimensionnel de la pièce avec une valeur de densité d'absorption des rayons X à chaque voxel. Ces images tridimensionnelles permettent de visualiser l'intérieur de la pièce d'une manière non destructive. ETAT DE L'ART Les pièces à inspecter 10 sont fabriquées en matériaux composites à l'aide de fibres de carbone tissées en trois dimensions et qui sont ensuite placés dans un moule puis injectés. La référence 10 s'applique 3047339 2 aussi bien à la pièce physique qu'à sa représentation tridimensionnelle dans le volume tomographique V. Après injection et usinage, certaines pièces présentent une double cambrure, appelée généralement « twist ». Les aubes de turbines 5 présentent généralement cette forme (voir figure 1) : une courbure liée aux profils intrados-extrados, en forme de lunule en croissant, et une courbure en hélice, qui s'enroule sur elle-même. A titre d'exemple, la torsion hélicoïdale s'étend sur 70°. Ces deux déformations s'étendent selon une direction, dite direction principale Z. 10 Pour l'analyse visuelle, l'opérateur utilise un logiciel d'affichage du volume tomographique V en trois dimensions qui lui permet de balayer l'intégralité de la pièce. Les figures 2a à 2c représentent les coupes bidimensionnelles de la pièce qui sont affichées à l'écran, chacune étant 15 orientée selon un axe de l'image brute. La figure 2d illustre ces trois plans sur une vue en trois dimension de la pièce 10. Les plans présentés en figures 2a et 2b sont facilement exploitables pour réaliser la recherche des anomalies de tissage puisque l'opérateur peut parcourir la pièce 10 et avoir une vision globale. En revanche, le 20 troisième plan, présenté en figure 2c, est plus difficile à analyser puisque le twist de la pièce ne permet pas d'avoir de vision globale et l'opérateur doit balayer cette vue à de multiples reprises. Lorsqu'une indication suspecte est détectée dans l'image, une analyse plus fine doit être menée pour confirmer ou infirmer la présence de cette indication. 25 Cela consiste à ajouter un plan qui est cette fois ci aligné en fonction de la pièce et non suivant les axes de l'image. Il s'agit d'une opération fastidieuse et très longue. Le document US 6 041 132 présente une méthode utilisant un modèle 30 simulé informatiquement, appelé modèle CAO, pour permettre de redresser le twist de la pièce. Néanmoins, l'utilisation du modèle CAO pour redresser le twist est un procédé lourd et contraignant : lorsque le modèle est disponible, il n'est 3047339 3 pas défini dans le même référentiel que le volume tomographique et il est nécessaire d'effectuer un « recalage ». Le calcul du recalage est une opération non triviale qui fait intervenir des méthodes complexes et coûteuses en ressources. En outre, il se peut que la pièce fabriquée ne 5 corresponde pas au modèle CAO : le redressement effectué ne sera alors pas conforme. La mise en oeuvre d'une telle méthode est donc complexe. PRESENTATION DE L'INVENTION 10 Pour palier les inconvénients de l'art antérieur, notamment en s'affranchissant de l'utilisation d'un modèle CAO, l'invention concerne un procédé de contrôle non-destructif d'une pièce aéronautique, comprenant une étape préliminaire d'obtention par tomographie 15 numérique d'un volume dans lequel se trouve une représentation de ladite pièce aéronautique, le volume étant défini par des voxels, et caractérisé en ce qu'il comprend un traitement de redressement comprenant les étapes suivantes : - (E21) Traitement de la représentation tomographique, ledit traitement 20 consistant à extraire dudit volume les voxels se rapportant uniquement à la pièce, - (E22) Génération d'une surface médiane de la pièce à l'intérieur du volume tomographique à partir dudit volume extrait se rapportant à la pièce, ladite surface médiane divisant la pièce en deux parties, 25 - (E23) Détermination d'un champ de vecteurs normaux à ladite surface médiane, - (E24) Aplanissement de ladite surface médiane pour former une surface médiane aplanie dans un plan, - (E25) Reconstruction d'un volume redressé contenant la pièce 30 redressée autour la surface médiane aplanie, ladite reconstruction étant effectuée par correspondance entre des voxels du volume redressé le long de directions orthogonales à la surface médiane aplanie et entre 3047339 4 des voxels du volume le long des vecteurs normaux de la surface médiane associés à ces directions orthogonales respectives, - (E3) Analyse du volume redressé ainsi obtenu pour identifier les anomalies de la pièce. 5 La représentation tridimensionnelle redressée est à présent orientée selon les axes principaux du repère, ce qui facilite la recherche d'anomalie sur un moniteur d'affichage, à l'aide des plans de coupe bidimensionnels présentés précédemment. Les avantages pour les 10 opérateurs sont multiples : analyse plus rapide, moins fatigante, moins d'opérations à effectuer, etc. La variabilité intra- et inter-opérateur est diminuée. Enfin, puisque la représentation tridimensionnelle est redressée, le volume redressé est moins gourmand en ressource informatique 15 d'affichage (mémoire vive, etc.), puisque le nombre de voxels liés à l'air est diminué (la taille du volume global est moins importante). L'invention peut comprendre les caractéristiques suivantes, prises seules ou en combinaison : 20 - l'étape de traitement de la représentation tomographique (E21) est réalisée par segmentation en deux classes des voxels en fonction d'un niveau de gris, 25 - la pièce s'étend selon une direction principale, et la surface médiane s'étend aussi selon cette direction principale, - l'étape de génération de la surface médiane est suivie d'une étape de filtrage passe-bas afin de supprimer certaines irrégularités, 30 - l'étape d'aplanissement (E24) se fait par projection de ladite surface médiane sur un plan, 3047339 5 - le plan de l'étape d'aplanissement (E24) est choisi de telle sorte que les distorsions liées à la projection sont minimisées, - le voxel (et même chaque voxel) du volume redressé est 5 caractérisé par sa distance, dite distance orthogonale, à un point de ladite surface médiane aplanie par projection orthogonale, dans lequel au point de la surface médiane aplanie correspond un point de la surface médiane, et dans lequel la valeur du voxel du volume situé à une distance le long de la direction définie par le 10 vecteur normal passant par le point est récupérée, ladite distance étant fonction de la distance dite orthogonale, - l'étape de reconstruction (E25) comprend les sous-étapes suivantes : 15 (E251) choix d'un voxel dans l'espace de la surface médiane aplanie, (E252) calcul de la distance par projection orthogonale, dite distance orthogonale, de ce voxel à la surface médiane aplanie, et détermination du point projeté sur la surface médiane aplanie, (E253) détermination du point de la surface médiane non-redressée 20 correspondant au point projeté de la surface non-redressée, (E254) récupération du vecteur normal passant par le point projeté correspondant sur la surface non-redressée, (E255) détermination de la valeur du voxel situé dans la direction du vecteur normal et à une distance fonction de la valeur de la distance 25 orthogonale, (E256) report de la valeur du voxel au voxel choisi dans la représentation redressée, - la fonction entre les deux distances est une fonction linéaire, 30 préférablement l'identité, - aucune modélisation informatique de ladite pièce aéronautique n'est utilisée, 3047339 6 L'invention propose aussi une unité de calcul comprenant des moyens de traitement de données et une mémoire, configurée pour mettre en oeuvre le traitement de redressement tel que défini précédemment. En 5 outre, un moniteur d'affichage peut être prévu, ledit moniteur affichant les données reçues par l'unité de calcul, et notamment les données du volume redressé. L'invention propose aussi un produit programme d'ordinateur, configuré 10 pour être exécuté par une unité de calcul et mettre en oeuvre le traitement de redressement tel que défini précédemment. PRESENTATION DES FIGURES 15 D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés, sur lesquels : - La figure 1 représente une aube twistée, - Les figures 2a à 2d, présentées précédemment, illustrent 20 différentes vues en coupe d'une pièce twistée, - Les figures 3 et 4 représentent différentes étapes d'un procédé selon des modes de réalisation de l'invention, - Les figures 5 et 6 représentent des illustrations d'axes médians et de surface médiane, respectivement en deux et trois dimensions, 25 - La figure 7 représente un champ de vecteurs normaux sur une surface médiane selon l'invention, - La figure 8 représente un plan de projection pour aplanir la surface médiane, - Les figures 9a à 9c représentent différentes étapes pour 30 reconstruire le volume redressé à partir de l'aplanissement de la surface médiane, - La figure 10 représente une reconstruction du volume redressé en cours, 3047339 7 - Les figures lia à lld représentent un volume redressé, avec une pièce redressée, - La figure 12 représente une unité de calcul et un moniteur d'affichage. 5 DESCRIPTION DETAILLEE A présent, en référence aux figures 3 et 4, une description détaillée du procédé de contrôle non-destructif par redressement va être donnée. 10 La pièce 10 sujette à l'inspection est destinée à une utilisation dans le domaine de l'aéronautique, et plus précisément dans le cadre de turbomachines. La description et les figures sont illustrées pour une aube telle que présentée en introduction, qui comprend une direction 15 principale Z. Néanmoins, le procédé est applicable à toutes formes de pièce, de forme non-torique. En effet, une telle forme pourrait générer des difficultés de génération de la surface médiane ou des difficultés de projection, comme cela sera explicité par la suite. 20 Afin de pouvoir effectuer un contrôle non-destructif à l'intérieur de la pièce 10, une étape préliminaire El d'acquisition à l'aide d'un système d'imagerie de tomographie numérique et de moyens de calcul (non représenté sur les figures) est effectuée. Cette étape permet d'obtenir une représentation tridimensionnelle de la pièce 10 pour l'aéronautique. 25 Le volume tomographie V obtenu est constitué de voxels Vx possédant chacun un niveau de gris. La pièce 10 est donc définie à l'intérieur du volume V par des voxels Vx ayant des niveaux de gris particulier. A présent, le procédé plus spécifique de redressement E2 va être décrit. 30 Ce procédé ne fait pas intervenir l'utilisation d'un modèle CAO. Par abus de langage, on parlera de la « pièce 10 » pour désigner la représentation tomographique de la pièce 10 qui possède une existence matérielle. La référence 10' s'applique à la « pièce 10 » redressée. 3047339 8 L'étape E2 est typiquement effectuée à l'aide d'une unité de calcul 100 comprenant des moyens de traitement de données 110 et une mémoire 120. Un moniteur d'affichage 130 peut être prévu, configuré pour 5 afficher des données notamment issues de l'unité de calcul 100. Génération de la surface médiane L'étape de génération de la surface médiane peut comprendre deux 10 étapes. Tout d'abord, dans une première étape E21, un traitement de la représentation est effectué pour extraire du volume les voxels correspondant à la pièce. On parlera de masque M de la pièce 10 dans le 15 volume tomographique V. Le masque M consiste à rassembler tous les voxels représentant uniquement la pièce. . On parle aussi d'extraction de la forme de la pièce 10. Le masque M est une image binaire pour laquelle chaque voxel qui la compose à un label 0 ou 1 : si le voxel appartient à la pièce 10, il a une valeur de 1 et si le voxel n'appartient 20 pas à la pièce 10, il a une valeur de 0. Un mode de réalisation consiste, à partir de la représentation obtenue par l'imagerie, à segmenter les voxels en deux catégories : air ou matière. Les voxels classés en « matière » sont ceux appartenant effectivement à la pièce 10 pour aéronautique, et les voxels classés en 25 « air » (ou « vide ») sont ceux n'appartenant pas à la pièce pour aéronautique, mais forme un volume d'air, complémentaire au volume de la pièce. Cette séparation peut être effectuée à l'aide des niveaux de gris, en choisissant un seuil de segmentation. 30 Une autre méthode consiste à calculer l'histogramme des niveaux de gris du volume tomographique et à maximiser la variance interclasse de cet histogramme afin de définir le seuil optimal qui sépare les voxels de 3047339 9 l'air et ceux de la matière (voir référence [1]). D'autres méthodes de seuillage automatique sont envisageables (voir référence [2]). Une fois le masque obtenu de l'étape E21 précédente, le procédé 5 comprend une étape de génération E22 d'une surface médiane SM à partir dudit masque M. La surface médiane SM est une surface qui divise la pièce 10 contenu dans le volume tomographique V en deux parties, c'est-à-dire plus exactement qui divise le masque M en deux volumes tridimensionnelles 10 (puisque la surface médiane SM est calculée à partir du masque M). Par surface médiane SM, on entend une surface qui, dans un plan transversal orthogonal, est à équidistance des deux extrémités de la pièce. La figure 5 illustre en deux dimensions une courbe médiane. La surface médiane SM est composée par l'ensemble des courbes 15 médianes, comme représentée en figure 6. En d'autres termes, la surface médiane SM passe par le milieu de la pièce 10. Afin d'avoir une surface médiane SM bien définie, il est préférable que la 20 pièce 10 présente un allongement dans une direction, ici la direction Z. La surface médiane SM s'étend alors selon cette direction Z. Une pièce torique pourrait générer une surface médiane torique ce qui complexifie les étapes suivantes du procédé. 25 La surface médiane SM peut être calculée de différentes manières. Des méthodes d'analyse automatique de la géométrie de la pièce (la surface médiane SM doit couper la pièce 10 en deux) conviennent, ou bien des algorithmes plus complexes permettant d'extraire l'axe médian en deux dimensions ou la surface médiane en trois dimensions (voir référence 30 [3] et [4]). Avantageusement, une étape postérieure Ef de filtrage est appliquée pour lisser et régulariser la surface médiane SM et limiter les 3047339 10 irrégularités. Un filtrage passe-bas, par exemple de type gaussien, peut être employé. L'importance d'une surface lisse réside dans l'extraction des vecteurs normaux, qui forme l'étape suivante. 5 Détermination d'un champ de vecteurs normaux Après l'étape E22 de génération de la surface médiane, une étape E23 de détermination d'un champ de vecteurs normaux N est effectuée (voir 10 figure 7). Par vecteur normal N, on entend un vecteur orienté et orthogonal à la surface médiane SM à l'emplacement de l'origine du vecteur. On privilégie des vecteurs unitaires, même si n'importe quelle norme peut convenir. Le champ de vecteurs N a pour but de connaitre la direction initiale de 15 recherche des différents voxels du volume V à insérer dans un volume redressé Vr, une fois la surface médiane SM aplanie. Plus d'explication seront données par la suite. Un exemple de calcul d'un vecteur normal N en un point H donné sur la 20 surface consiste à extraire un ensemble de points de la surface très proches du point H et de calculer l'inertie de ce nuage de points, c'est-à-dire les variations dans l'espace dans un voisinage. Après avoir calculé les valeurs propres et les vecteurs propres de la matrice d'inertie, on peut en déduire le vecteur normal. 25 En effet, comme le nuage de point extrait correspond à une surface qui localement s'apparente à un plan, les valeurs (À) et vecteurs propres (v,) avec i E {1,2,3} : (v1,À1), (v2,À2) et (v3,À3) sont organisés de la manière suivante : À1'k,'À2>>À3. Le vecteur v3 correspond donc à la normale à la surface recherchée au point H. Le même processus est 30 réalisé pour tous les points H appartenant à la surface médiane et le champ de vecteurs N est ainsi déterminé. 3047339 11 Les données relatives au champ de vecteurs sont stockées dans une mémoire et seront utilisées par la suite. Aplanissement de la surface médiane 5 A présent, la pièce 10 subit des transformations pour effectuer le redressement. Dans une étape E24 d'aplanissement, la surface médiane SM qui est sous forme tridimensionnelle est aplanie dans un plan P, c'est-à-dire 10 sous forme bidimensionnelle (voir figure 8, dans laquelle la pièce 10 est représentée, avec le plan P, et non pas la surface médiane SM). La surface médiane SM issue de cette étape E24 est appelée surface médiane aplanie SMa. Pour cela, une méthode consiste à choisir un plan, qui devient le plan P, 15 et à projeter la surface médiane SM tridimensionnelle dans ce nouveau référentiel. La référence à la forme non-torique de la pièce 10 mentionnée précédemment se manifeste ici par le fait qu'il est préférable que chaque point projeté dans le plan P n'ait qu'un seul antécédent par côté du plan P ainsi que chaque point du volume n'ait 20 qu'une seule projection possible. Dans le cas contraire, une incertitude serait créée lors de la reconstruction du volume redressé, comme cela sera explicité par la suite. Avantageusement, le plan est choisi de sorte qu'il minimise les distorsions liées à la projection ; complémentaire, ledit plan peut aussi 25 être choisi pour qu'il minimise l'ensemble des distances projetées. On parle désormais de surface médiane aplanie SMa pour désigner la surface obtenue suite à l'étape E24, et de surface médiane SM non aplanie obtenue suite à l'étape E22. 30 Redressement du volume 3047339 12 Enfin, à partir de la surface médiane aplanie SMa, le volume tomographique doit être reconstruit tout autour. On parle alors de volume redressé Vr. 5 Cette étape E25 de reconstruction fait intervenir les étapes E22 de génération de la surface médiane, E23 de génération du champ de vecteurs N et E24 d'aplanissement. Dans cette étape E25, une correspondance est établie entre un voxel Vx' 10 du volume redressé Vr, le long d'une direction orthogonale à la surface médiane aplanie SMa, et un voxel Vx du volume tomographique V associé, le long du vecteur normal N correspondant à cette direction orthogonale. Par correspondance, on entend que la valeur du niveau de gris du voxel Vx est implémentée dans à la valeur du voxel Vx'. 15 On parle ainsi de volume tomographique redressé Vr pour désigner le volume tomographique V obtenu à l'issue de l'étape E25, et de volume tomographique non redressé pour désigner le volume tomographique obtenu à l'issue de l'étape El. 20 En d'autres termes, une correspondance est établie entre un voxel Vx' du volume redressé Vr et un voxel Vx du volume V non redressé. Pour cela, le voxel Vx' est caractérisé par sa distance h' à un point A' de ladite surface médiane aplanie SMa par projection orthogonale (c'est à dire le long d'une direction orthogonale à la surface médiane aplanie 25 SMa et passant par le point A). A ce point A' correspond un point A de la surface médiane SM et la valeur du voxel Vx disposé à une distance h fonction de la distance h' le long de la direction définie par le vecteur normal N passant par le point A est récupérée. La définition du point A' et de sa correspondance A est typiquement un 30 intermédiaire pour pouvoir associer la direction orthogonale selon laquelle la distance h' est mesurée à un vecteur normal N selon lequel la distance h est mesurée. 3047339 13 Comme mentionné précédemment, s'il existe plusieurs points dont la projection correspond aux points A ou A', des incertitudes sont créées qui complexifient la reconstruction du volume. 5 On parle désormais de volume tomographique redressé pour désigner le volume tomographique obtenu à l'issue de l'étape E25, et de volume tomographique non redressé pour désigner le volume tomographique obtenu à l'issue de l'étape El. 10 En balayant le volume redressé Vr à l'aide des voxels Vx', l'ensemble des valeurs de niveau de gris du volume tomographique V est réintégrée voxel après voxel dans le volume redressé, ce qui génère in fine une image tridimensionnelle redressée. 15 Il se peut que la correspondance d'un voxel Vx' du volume redressé ne tombe pas exactement sur un voxel Vx du volume non redressé mais entre huit voxels Vx. En pratique, une interpolation trilinéaire est effectuée à chaque correspondance. 20 Les figures 9a à 9c et 10 illustrent en deux dimensions un mode de réalisation de l'étape de reconstruction. Dans une sous-étape E251, un voxel Vx' dans le volume redressé Vr de la surface médiane aplanie SMa est sélectionné. 25 Dans une sous-étape E252, la distance h' par projection orthogonale sur ladite surface médiane aplanie est calculée. Cette distance h' est appelée distance orthogonale. Parallèlement, le point A' de la surface médiane aplanie correspondant à la projection de Vx' est calculé. Dans une sous-étape E253, le point A de la surface médiane SM non 30 redressée correspondant au point A' de la surface médiane aplanie SMa est déterminé. Pour cela, il est possible d'effectuer une étape intermédiaire de correspondance via le point correspondant de la surface médiane SM non aplanie. 3047339 14 La figure 9a illustre les étapes précitées, qui se déroulent dans le volume redressé Vr. Dans une sous-étape E254, le vecteur normal N passant par le point A 5 est récupéré. Dans une sous-étape E255, la valeur du voxel Vx (dans le volume tomographique V non redressé) situé le long de la direction dudit vecteur normal N à une distance h fonction de la distance orthogonale h' est récupérée. 10 La figure 9b illustre les étapes précitées qui se déroulent dans le volume V. Enfin, dans une sous-étape E256, la valeur du voxel Vx du volume V est assigné au voxel Vx' du volume redressé Vr. 15 La figure 9c illustre cette étape qui se déroule à nouveau dans le volume redressé. Préférablement, pour limiter les déformations, la fonction liant h et h', c'est-à-dire liant les distances AVx et A'Vx', est une fonction linéaire de 20 type h=bh'. Enfin, encore préférablement, on a b=1, ce qui signifie que la fonction est l'identité et que h'=h : de la sorte, il suffit de reporter la même distance. Plus généralement, la fonction doit permettre de conserver la topographie des défauts. 25 Comme indiqué précédemment, les étapes E251 à E256 sont itératives pour que le volume V soit redressé en intégralité. La figure 10 illustre cette reconstruction par itération : la pièce redressée 10' est reconstruite voxel Vx' par voxel Vx', par récupération des valeurs des 30 voxels Vx respectivement associés. Pour effectuer cette itération, les étapes E251 et E252 peuvent être répétées pour tous les points situés sur la même droite, de sorte que les étapes E253 et E254 ne sont effectuées qu'une seule fois. En revanche, 3047339 15 les étapes E255 et E256 sont répétées autant de fois que les étapes E251 et E252 l'ont été. Lors de la reconstruction du volume, il se peut que le point Vx corresponde en fait dans le volume non redressé à un voxel Vx' d'air. 5 Une valeur nulle peut par exemple être accordée à Vx. Pour éviter des calculs inutiles, le procédé peut établir que les correspondances suivantes seront effectuées pour des distances h inférieures. Alternativement, le choix d'un voxel Vx peut être fait en fonction d'une direction orthogonale prédéterminée auparavant. Cela signifie que dans 10 une étape E251', il s'agit du point A qui est fixé et que la droite orthogonale à la surface médiane redressée passant par A est calculée. Ensuite, dans une étape E252', une distance h est choisie, déterminant ainsi un voxel Vx. Dans ce mode de réalisation, il est possible d'itérer le procédé en 15 répétant l'étape E252', préférablement en faisant augmenter la valeur de h, de façon à balayer le volume en direction de l'extérieur de la pièce 10. Dès lors qu'un voxel Vx est associé à un voxel Vx' d'air, le procédé peut considérer que la pièce a été intégralement reconstruite selon cette direction et les étapes E251', etc. sont ensuite répétées pour un autre 20 point A. A l'issue de l'étape de reconstruction E25, un volume tomographique dit redressé Vr a été généré. Les figures lia à 11c représentent trois différents plans en coupe orthogonaux deux à deux et la figure iid 25 représentent la pièce redressée traversées par ces trois plans. En comparaison avec la figure 3c, ce plan issu du redressement est dorénavant exploitable sans complication particulière. Du fait des transformations impliquées dans le procédé (projection et report de distance), la topographie des défauts est respectée, ce qui 30 permet d'assurer qu'aucun défaut ne disparait lors de l'application du précédé : l'analyse par un opérateur sur le volume tomographique redressé se trouve plus fiable, plus stabilisé et plus rapide. 3047339 16 En outre, du fait du redressement de la pièce 10 dans le volume tomographique, la taille du fichier est diminuée puisque le nombre de voxels d'air est diminué. 5 Enfin, une étape E3 d'analyse telle que définie en introduction est avantageusement effectuée par un opérateur, à l'aide en particulier du moniteur 130 qui reçoit les données de l'unité de calcul 100 (figure 12). Cette étape est désormais plus rapide, plus efficace, plus fiable, moins gourmande en ressource de calcul et possède une répétabilité 10 inter-opérateur accrue. On définit aussi un produit programme d'ordinateur qui peut être exécuté par les moyens de traitement de données 110 de l'unité de calcul 100 et qui permet de mettre en oeuvre le traitement de 15 redressement décrit précédemment. En particulier, l'image tridimensionnelle redressée est stockée dans la mémoire 120. Références : 20 [1] Nobuyuki Otsu, « A threshold selection method from gray-level histograms », IEEE Trans. Sys., Man., Cyber., vol. 9, 1979, p. 62-66 (DOI 10.1109/TSMC.1979.4310076). [2] M. Sezgin and B. Sankur, « Survey over image thresholding 25 techniques and quantitative performance evaluation », Journal of Electronic Imaging, vol. 13, no 1, 2003, p. 146-165 (DOI 10.1117/1.1631315). [3] Hesselink, W. H., & Roerdink, J. B. (2008). Euclidean skeletons of 30 digital image and volume data in linear time by the integer medial axis transform. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 30(12), 2204-2217. 3047339 17 [4] Chaussard, J., Couprie, M., & Talbot, H. (2009, January). A discrete À-medial axis. In Discrete Geometry for Computer Imagery (pp. 421433). Springer Berlin Heidelberg. REVENDICATIONS1. Procédé de contrôle non-destructif d'une pièce aéronautique (10), comprenant une étape préliminaire (El) d'obtention par tomographique numérique d'un volume (V) dans lequel se trouve une représentation de ladite pièce aéronautique (10), le volume (V) étant défini par des voxels (Vx), et caractérisé en ce qu'il comprend un traitement de redressement (E2) comprenant les étapes suivantes : - (E21) Traitement de la représentation tomographique, ledit traitement consistant à extraire dudit volume les voxels se rapportant uniquement à la pièce - (E22) Génération d'une surface médiane (SM) de la pièce (10) à l'intérieur du volume tomographique à partir dudit volume extrait se rapportant à la pièce, ladite surface médiane (SM) divisant la pièce (10) en deux parties, - (E23) Détermination d'un champ de vecteurs normaux (N) à ladite surface médiane (SM), - (E24) Aplanissement de ladite surface médiane (SM) pour former une surface médiane aplanie (SMa) dans un plan (P), - (E25) Reconstruction d'un volume redressé (Vr) contenant la pièce (10, 10') redressée autour la surface médiane aplanie (SMa), ladite reconstruction étant effectuée par correspondance entre des voxels (Vx') du volume redressé (Vr) le long de directions orthogonales à la surface médiane aplanie (SMa) et entre des voxels (Vx) du volume (V) le long des vecteur normaux (N) de la surface médiane (SM) associés à ces directions orthogonales respectives, le procédé comprenant en outre une étape (E3) d'analyse du volume redressé ainsi obtenu pour identifier les anomalies de la pièce. 2. Procédé de contrôle non-destructif selon la revendication précédente dans lequel l'étape (E22) de génération de la surface médiane (SM) est suivie d'une étape (Ef) de filtrage passe-bas afin de supprimer certaines irrégularités. 3047339 19 3. Procédé de contrôle non-destructif selon l'une quelconque des revendications précédentes, dans lequel l'étape d'aplanissement (E24) se fait par projection de ladite surface médiane sur le plan (P). 5 4. Procédé de contrôle non-destructif selon la revendication précédente, dans lequel le plan de l'étape d'aplanissement (E24) est choisi de telle sorte que les distances des projections sont minimisées. 10 5. Procédé de contrôle non-destructif selon l'une quelconque des revendications précédentes, dans lequel le voxel (Vx') du volume redressé (Vr) est caractérisé par sa distance (h'), dite distance orthogonale, à un point (A') de ladite surface médiane aplanie (SMa) par projection orthogonale, dans lequel au point (A') de la surface médiane 15 aplanie (SMa) correspond un point (A) de la surface médiane (SM), et dans lequel la valeur du voxel (Vx) du volume (V) situé à une distance (h) le long de la direction définie par le vecteur normal (N) passant par le point (A) est récupérée, ladite distance (h) étant fonction de la distance dite orthogonale (h'). 20 6. Procédé de contrôle non-destructif selon la revendication précédente, dans lequel la fonction liant les deux distances (h, h') est une fonction linéaire, préférablement l'identité. 25 7. Procédé de contrôle non-destructif selon l'une quelconque des revendications précédentes, dans lequel aucune modélisation informatique de ladite pièce aéronautique (10) n'est utilisée. 8. Unité de calcul (100) comprenant des moyens de traitement de 30 données (110) et une mémoire (120), configurée pour mettre en oeuvre le traitement de redressement tel que défini dans l'une quelconque des revendications 1 à 7. 3047339 20 9. Produit programme d'ordinateur, configuré pour être exécuté par une unité de calcul (100) et mettre en oeuvre le traitement de redressement tel que défini dans l'une quelconque des revendications 1 à 7.