FR 3141245 A1 20240426 FR 2210787 A 20221019 Procédé d’analyse tomographique d’une pièce, le procédé comprenant les étapes suivantes : - acquisition d’au moins une image tomographique au moyen d’un dispositif de tomographie, - division de l’image tomographique en sous-parties élémentaires et obtention d’une pluralité d’images (32) bidimensionnelles pour chaque sous-partie, - description des images (32) par application d’au moins un outil de traitement numérique sur chaque image de la sous-partie et obtention d’une valeur d’au moins un paramètre représentatif, - détermination de variations spatiales dans chaque sous-partie, selon trois directions (X, Y, Z), des valeurs de chaque paramètre et détection de valeurs hors de plages de conformité, - calcul d’un score de conformité de chaque sous-partie, et - comparaison du score de conformité de chaque sous-partie à un seuil de conformité prédéterminé pour déterminer une conformité de la sous-partie. Figure à publier avec l’abrégé : 5 Procédé d’analyse tomographique Domaine technique de l’invention L’invention concerne un procédé d’analyse tomographique d’une pièce mécanique. La pièce fait par exemple partie d’un aéronef, et plus précisément d’une turbomachine, comme une partie d’un carter de soufflante, une aube de soufflante, ou une partie d’une structure d’aubages fixes. Etat de la technique antérieure Il est connu d’inspecter les pièces mécaniques structurelles par des méthodes de tomographie, par exemple radiographique ou acoustique, afin de détecter la présence d’éventuels défauts en surface et à l’intérieur de la pièce. Ces méthodes sont fiables, non-invasives et permettent d’inspecter l’intérieur des pièces, ce qui permet de statuer rapidement sur l’état de la pièce et décider de sa potentielle remise en fonctionnement ou son remplacement. Une mesure tomographique consiste en un balayage d’un objet observé, ici une pièce mécanique, au moyen d’un faisceaux d’ondes, et en une mesure du faisceau transmis dans toutes les directions afin de reconstituer une image tridimensionnelle de l’objet. La illustre un dispositif de tomographie 1 pendant l’inspection d’une pièce mécanique 20. Le dispositif de tomographie 1 comprend au moins un dispositif d’émission 3, configuré pour émettre un faisceau incident 5 d’impulsions d’ondes en direction de la pièce 20, par exemple des ondes radiofréquence, des rayons X ou encore des ondes acoustiques, et un moins un récepteur 7 aptes à capter un faisceau transmis 9 d’ondes, agencés de part et d’autre de la pièce 20. La pièce 20 est généralement montée sur un support 11 rotatif autour d’un axe A, afin d’être observée sous toutes les directions pendant la mesure tomographique. Il peut exister un montage inverse, où la pièce est fixe et le couple dispositif d’émission / récepteur tournent autour de celle-ci. La pièce 20 observée par tomographie peut être constituée d’un matériau composite, tel que représenté en détail sur la . Un tel matériau composite peut présenter un tissage dit bidimensionnel, qui comprend des fibres de trames généralement orientées selon une direction de trame X et des fibres de chaines généralement orientées selon une direction de chaine Y, tissées les unes avec les autres sur un plan de tissage P, ou plusieurs plans de tissage P superposés selon une direction d’épaisseur Z. Les fibres peuvent être des fibres de carbone, de verre ou un mélange des deux. Elles sont noyées dans une matrice 23 comprenant par exemple un ou plusieurs polymères et/ou résines, qui se solidifient pour former la pièce 20 définitive. Alternativement, comme c’est le cas dans l’exemple de la , le tissage des fibres 21 peut être tridimensionnel, avec des fibres 21 s’étendant selon les trois directions de l’espace X, Y, Z. Les amplitudes des ondes transmises sous les différents angles d’observations sont traduites en niveaux de gris dans un dispositif de traitement 13, et une analyse numérique permet de reconstituer une image tridimensionnelle en niveaux de gris du volume de la pièce 20. Le dispositif de traitement 13 comprend un ordinateur. L’analyse par tomographie permet alors d’obtenir un visuel sur l’extérieur et l’intérieur de la pièce. Ce type de contrôle présente l’avantage majeur de permettre la visualisation dans l’épaisseur du matériau et de la pièce étudiée, tout en étant non-destructif et fiable. Un opérateur humain peut alors évaluer la qualité de la pièce, en recherchant dans l’image tridimensionnelle la présence ou l’absence d’une anomalie ou d’un dégât, et en évaluant la conformité du matériau aux normes applicables. Par exemple, dans le cas d’une pièce en composite tissé, l’anisotropie du tissage, la continuité des fibres et celle de la matrice sont des paramètres pertinents évalués visuellement par l’opérateur. Cependant, l’analyse par un inspecteur qualifié d’une pièce comme une aube de soufflante est chronophage, cette opération « manuelle » et fastidieuse, une grande partie de ce temps consistant à faire défiler les images de tomographie. L’opération génère également une fatigue visuelle importante qui peut réduire la vigilance et la qualité d’analyse de l’inspecteur après un certain temps. De plus, d’un point de vue technique, la mise en œuvre d’une telle analyse par des êtres humains peut engendrer des biais d’appréciation et in fine potentiellement une grande variabilité de résultat. Présentation de l’invention L’invention vise à remédier à ces inconvénients, en fournissant un procédé d’analyse permettant une évaluation rapide et systématique de l’état d’une pièce mécanique sans nécessiter l’intervention d’un inspecteur humain. Cela permet de rendre l’inspection plus rapide et réduit les risques d’erreur. A cet effet, l’invention a pour objet un procédé d’analyse tomographique d’une pièce, le procédé comprenant les étapes suivantes : - acquisition d’au moins une image tomographique de la pièce au moyen d’un dispositif de tomographie, - division de l’image tomographique en sous-parties élémentaires et obtention d’une pluralité d’images bidimensionnelles pour chaque sous-partie, chaque image comprenant une pluralité de pixels, chaque pixel présentant un niveau de gris respectif, - application d’au moins un outil de traitement numérique sur chaque image de la sous-partie et obtention, pour chaque outil appliqué, d’une valeur d’au moins un paramètre représentatif pour chaque image, - détermination de variations spatiales dans chaque sous-partie, selon trois directions caractéristiques, des valeurs de chaque paramètre et détection de valeurs des paramètres hors de plages de conformité prédéterminées, - calcul d’un score de conformité de chaque sous-partie, à partir des valeurs des paramètres hors des plages de conformité, et - comparaison du score de conformité de chaque sous-partie à un seuil de conformité prédéterminé pour déterminer une conformité de la sous-partie. Un tel procédé permet une analyse de la totalité d’une pièce mécanique, la détection de non-conformités, et la formulation d’un diagnostic global de l’état de la pièce sans nécessiter la présence d’un inspecteur humain. Le procédé peut comprendre une étape de traitement des images bidimensionnelles, mise en œuvre avant l’étape d’application d’au moins un outil de traitement numérique, au cours de laquelle chaque pixel dont le niveau de gris est hors d’une plage médiane prédéterminée est remplacé par un pixel corrigé présentant un niveau de gris compris dans ladite plage médiane. Une telle étape permet d’éliminer les pixels dont les niveaux de gris présentent des valeurs extrêmes, et qui pourraient générer de forts gradients et fausser les calculs suivants, en étant considérés comme des singularités par les outils numériques. Le niveau de gris corrigé peut être déterminé aléatoirement à partir d’une distribution des niveaux de gris des pixels de l’image dont le niveau de gris est dans la plage médiane. Une telle caractéristique permet de fondre les pixels corrigés dans le reste de l’image, avec une distribution de type « bruit blanc ». L’au moins un outil de traitement numérique peut comprendre un calcul d’autocorrélation, comprenant la génération d’une image d’autocorrélation à partir de l’image, la détection d’un pic central de l’image d’autocorrélation et la détermination d’un rapport d’isotropie du pic central, le rapport d’isotropie étant choisi comme valeur de paramètre représentatif. Un tel outil permet de qualifier l’anisotropie d’une image, et par exemple de détecter les directions de tissage préférentielles d’une pièce composite. L’au moins un outil de traitement numérique peut comprendre un calcul de filtrage, comprenant l’application à l’image d’au moins un filtre de Gabor et l’obtention d’une image filtrée, la détermination, pour chaque pixel de l’image d’une somme des valeurs correspondantes des images filtrées avec chacun des filtres de Gabor, le maximum sur l’image des sommes étant choisi comme valeur de paramètre représentatif. Un tel outil permet de détecter des anomalies d’entrelacement locales. Plusieurs filtres de Gabor peuvent être appliqués, avec des bases respectives obtenues à partir des directions caractéristiques de l’image par rotation. Dans un mode de réalisation quatre filtres de Gabor peuvent être appliqués, avec des bases respectives obtenues à partir des directions caractéristiques de l’image par rotation par exemple d’un angle respectivement de 0°, 45°, 90° et 135° pour une structure tissée à 0° et 90 °. Une telle caractéristique permet de traiter toutes les orientations de l’image. L’au moins un outil de traitement numérique peut comprendre un calcul d’un écart-type des niveaux de gris sur l’ensemble des pixels de l’image, le dit écart-type étant choisi comme valeur de paramètre représentatif. Un tel outil permet de détecter une surreprésentation locale d’une phase du matériau, par exemple un amas de résine dans le cas d’une pièce en composite tissé. Trois outils de traitement numériques différents peuvent être appliqués sur chaque image, l’étape de calcul du score de conformité comprenant le calcul d’un score partiel pour chaque outil de traitement numérique et pour chaque direction, à partir de la quantité de valeurs des paramètres hors des plages de conformité obtenue pour ledit outil de traitement numérique et ladite direction, le score de conformité étant la somme des scores partiels. Une telle caractéristique permet d’observer plusieurs facteurs différents en vue d’établir un diagnostic plus précis quant à l’état de la pièce. Les plages de conformité pour les paramètres sont prédéterminées à partir d’éprouvette saines, chaque plage de conformité comprenant les 99 centiles inférieurs des valeurs du paramètre associé mesurée sur lesdites éprouvettes saines. Une telle caractéristique permet de détecter de manière fiable des valeurs non-conformes afin d’améliorer la précision de l’analyse. L’invention concerne également un produit programme d’ordinateur comprenant des instructions permettant, lorsqu’il est exécuté sur un ordinateur pilotant un dispositif d’analyse tomographique, la mise en œuvre du procédé comme plus haut. Brève description des figures la est une vue schématique de côté d’un dispositif de tomographie au cours de la mise en œuvre d’un procédé selon l’invention, la est une vue schématique de détail d’une pièce en matériau composite tissé, est une vue schématique d’une étape de division d’une image tomographique en sous-parties élémentaires, et est une vue schématique d’application d’une correction des niveaux de gris extrême par remplacement par bruit blanc est une vue schématique d’utilisation de la fonction d’autocorrélation pour mise en évidence d’une indication de composite tissé Description détaillée de l’invention Un procédé d’analyse tomographique selon l’invention va maintenant être décrit. Ce procédé met en œuvre le dispositif d’analyse 1 décrit précédemment, et vise à détecter la présence d’anomalies dans une pièce 20, réalisée par exemple en matériau composite. La pièce est ici une partie d’une pale d’aubage, par exemple d’une aube de soufflante. Ladite aube de soufflante présente une direction d’élongation principale choisie comme direction Z de l’espace, la direction X étant choisie sensiblement tangente à la surface externe de l’aube et la direction Y étant la direction d’épaisseur de l’aube. Une épaisseur de la pièce 20, mesurée selon la direction Y, est par exemple compris entre 5 et 25 mm, notamment proche de 10 mm. Le procédé comprend une première étape d’acquisition d’au moins une image tridimensionnelle de la pièce 20 au moyen du dispositif de tomographie 1. Un faisceau incident 5 d’impulsions d’ondes est émis par le dispositif d’émission 3, en direction de la pièce 20. Le faisceau incident 5 est par exemple un faisceau d’ondes radiofréquences. Les ondes traversent la pièce 20, et un faisceau transmis 9 issu de la pièce 20 est capté par le récepteur 7. La distribution d’intensité du faisceau transmis 9 obtenue est convertie en une image bidimensionnelle en niveaux de gris de la pièce 20 par le dispositif de traitement 13. La pièce 20 est mise en rotation au moyen du support 11, et des images bidimensionnelles de la pièce 20 sont acquise sous toutes les directions. Une image tridimensionnelle 30 de la pièce 20 en niveaux de gris, ou image tomographique, est alors reconstituée par traitement d’image au moyen du dispositif de traitement 13. De manière non limitative, la pièce 20 est une pièce en matériau composite à tissage tridimensionnel, telle que représentée sur la . Le matériau composite comprend des fibres 21 selon les trois directions X, Y et Z de l’espace, noyées dans une matrice 23. L’écart moyen entre deux fibres 21 voisines et de même direction est de l’ordre de 2 mm. Les fibres 21 peuvent être de matériaux identiques ou différents (verre, kevlar ou carbone par exemple). La matrice 23 comprend au moins un polymère organique et/ou au moins une résine. Les fibres de différents types et la matrice sont désignées sous le terme de phases du matériau composite composant la pièce 20. L’image tridimensionnelle 30 comprend donc une pluralité de voxels présentant chacun un niveau de gris respectif représentatif notamment de la densité de la phase du matériau correspondant dans la pièce 20. Le procédé d’analyse tomographique comprend alors une étape de division de l’image tridimensionnelle 30 tomographique en sous-parties 31 élémentaires, comme représenté sur la . Ces sous-parties 31 ont des dimensions, selon trois directions de l’espace X, Y, Z, calculées à partir des caractéristiques du matériau de la pièce 20. Les sous-parties 31 sont par exemple de forme parallélépipédique. Ces dimensions des sous-parties 31 peuvent varier en fonction de la localisation de la sous-partie dans la pièce 20, comme représenté sur la qui est une coupe plane d’une cartographie des sous-parties 31 dans un plan XZ de la pièce 20. Les dimensions des sous-parties 31 sont par exemple choisies comme suit. Selon la direction Y, la dimension de chaque sous-partie 31 est choisie égale à l’épaisseur de la pièce 20 selon cette direction Y. Selon la direction Z, la dimension de chaque sous-partie 31 est choisie égale à un nombre entier n d’espacements entre colonnes de tissage (également désigné par le terme pick spacing en anglais, PS), c’est-à-dire un nombre entier de pas entre fibres 21 voisines espacées selon la direction Z. Le nombre entier n est choisi à partir de résultats expérimentaux et de modélisation obtenus avec le matériau, et correspond par exemple à une longueur totale comprise entre 15 mm et 25 mm. Selon la direction X, la dimension de chaque sous-partie 31 est choisie afin que la sous-partie 31 constitue un volume élémentaire représentatif (VER). Les dimensions de sous-parties 31 selon les directions X et Z sont par exemple voisines l’une de l’autre (par exemple comprise entre la moitié et le double l’une de l’autre). Le procédé selon l’invention ayant pour objectif de simuler une analyse par un inspecteur, une pluralité d’images 32 bidimensionnelles est ensuite obtenue pour chaque sous-partie 31, simulant des coupes du volume de la sous-partie 31. Les images 32 sont des coupes planes en niveau de gris, et chaque image comprend ainsi une pluralité de pixels, chaque pixel présentant un niveau de gris respectif. Ces images peuvent être des coupes de l’image tomographique dans des plans perpendiculaires aux directions X, Y et Z de l’espace, ou dans des plans obliques (c’est-à-dire non orthogonaux aux directions X, Y et Z). Le nombre d’images 32 obtenues pour chaque sous-partie 31 est prédéterminé de manière empirique, selon le nombre de coupes pertinentes pour l’analyse de la conformité de la sous-partie 31. L’ensemble de ces images 32 constituent les données brutes à traiter, représentatives du matériau et de la pièce 20 à analyser, par exemple pour la mise en évidence des anomalies d’entrelacement des fibres 21. Le procédé comprend avantageusement une étape de traitement des images 32, représentée sur la , pour retirer les niveaux de gris parasites susceptible des fausser le traitement numérique postérieur. Les pixels de l’image 32 dont le niveau de gris est hors d’une plage médiane prédéterminée sont remplacés par des pixels dont le niveau de gris est compris dans cette plage médiane. La plage médiane s’étend par exemple entre 10% et 90% de la pleine échelle de niveaux de gris du dispositif d’acquisition. Les pixels présentant un niveau de gris compris entre 0 et 10% et entre 90% et 100% du maximum de l’échelle (0 correspondant généralement au noir et 100% correspondant généralement au blanc) sont donc considérés comme parasite et remplacés. Ces valeurs de 10% et 90% sont données à titre indicatif et peuvent varier selon les cas, en fonction de la nature du matériau et de la sensibilité des capteurs du dispositif d’acquisition 1, par exemple. Dans l’exemple représenté , les pixels noirs de l’image 32 non traitée, à gauche, qui correspondent à un vide, sont considérés comme présentant des niveaux de gris extrêmes. La valeur de niveau de gris de remplacement est avantageusement déterminée en fonction de la distribution de niveau de gris de l’image 32, par exemple choisie égale à la médiane de niveaux de gris des pixels de l’image. Plus avantageusement, le niveau de gris de chaque pixel remplacé de l’image 32 est déterminé aléatoirement sur la base de la distribution de niveaux de gris des pixels de l’image 32 situés dans la plage médiane. Ainsi, les pixels remplacés forment un bruit blanc dans l’image et ne modifient pas les étapes de traitement suivantes. Un tel remplacement par une valeur aléatoire de niveau de gris est représenté sur la , où l’image 32 de gauche a vu ses pixels extrêmes remplacés par un « bruit blanc » déterminé à partir de la distribution des niveaux de gris non-extrêmes de l’image 32 initiale. Le procédé comprend ensuite une étape de description des images 32 par application d’au moins un outil de traitement numérique sur chaque image 32 de la sous-partie 31 et obtention, pour chaque outil appliqué, d’une valeur d’au moins un paramètre représentatif pour chaque image 31. Avantageusement, au moins trois tels outils de traitement numériques sont appliqués sur chaque image 32, de manière indépendante, et l’obtention d’au moins une valeur d’un paramètre pour chacun de ces outils. Par exemple, le procédé comprend un calcul d’autocorrélation et la détermination d’un rapport d’anisotropie pour chaque image 32, une application d’au moins un filtre de Gabor et la détermination d’une somme des maximas des images filtrées pour chaque image 32, et un calcul d’écart-type des niveaux de gris et l’obtention d’une valeur d’écart type pour chaque image 32. Le calcul d’autocorrélation est une méthode utilisée dans la littérature afin de rendre compte de l’anisotropie, de la texture ou de la périodicité d’une image. Un tel calcul est représenté sur la . Une image d’autocorrélation 33 est générée pour chacune des images 32, comprenant un pic central 34. Un seuil est calculé à partir des valeurs des maximas locaux de l’image d’autocorrélation 33 en excluant le pic central 34, et le pic central 34 est isolé grâce à ce seuil, pour obtenir une image d’isolation 35 comprenant uniquement le pic central 34. L’anisotropie de l’image 32 peut alors être caractérisée par le calcul des valeurs propres et des vecteurs propres de la matrice d’inertie de l’image d’isolation 35 du pic central 34. Un rapport d’anisotropie est calculé comme rapport des longueurs des vecteurs propres, et caractérise la présence d’une direction préférentielle de l’image 32. Des valeurs proches de 1 du rapport d’anisotropie caractérisent des images 32 où les deux directions principales ont des poids similaires, alors que des valeurs éloignées de 1 dénotent une anisotropie. C’est cette orientation préférentielle de l’image 32 qui est mesurée car impactée dans le cas de présence d’anomalie d’entrelacement de fibres 21. Un filtre de Gabor est un filtre linéaire dont la réponse impulsionnelle est une sinusoïde modulée par une fonction gaussienne. La fréquence de la sinusoïde pourrait par exemple être choisie égale à 0,05 (en pixel -1 ), et l’écart-type σ 1 de la gaussienne est calculé selon l’équation : . Avantageusement, quatre filtres de Gabor sont appliqués successivement à l’image 32, avec des bases respectives obtenus par des rotations depuis les directions principales de l’image 32 d’angles qui pourrait être par exemple égaux à 0°, 45°, 90° et 135°. Pour chaque pixel de l’image 32, la somme des valeurs obtenues avec les quatre filtres de Gabor est calculée, et la somme la plus élevés parmi les pixels de l’image est choisie comme valeur de paramètre représentatif de l’application de l’outil sur l’image 32. Cette valeur maximale peut être comparée à des valeurs issues d’une base de données, et prend des valeurs très importantes dans le cas de présence d’anomalie d’entrelacement à signature linéaire ou non linéaire. Le calcul de l’écart type de la distribution de niveaux de gris pour chaque image 32 est calculé, et la valeur d’écart-type est la valeur de paramètre représentatif. Cette valeur est caractéristique d’un amas non négligeable d’une des phases du matériau, par exemple un amas de résine, linéaire ou non. En effet, le nombre de pixels à niveau de gris faible (correspondant à la phase résine) va avoir pour conséquence une augmentation de l’écart type des niveaux de gris dans l’image 32. Le procédé comprend alors une étape de détermination de variations spatiales dans chaque sous-partie 31, selon les trois directions X, Y, Z, des valeurs de chaque paramètre. Une telle représentation des variations spatiales permet de simuler l’inspection d’une sous-partie 3A par défilement des coupes par un inspecteur. La dynamique des images 32 joue un rôle important dans la détection des anomalies d’entrelacement. Ainsi, pour chacun des paramètres caractéristiques des outils numériques précédemment appliqués, et pour chacune des directions d’analyse X, Y et Z, une courbe d’évolution des valeurs du paramètres est tracée en fonction des positions des images dans la sous partie 31 selon la direction. Le procédé comprend ensuite la détection de valeurs non-conformes de ces paramètres, qui se situent hors de plages de conformité prédéterminées. Les plages de conformité sont obtenues, pour chaque paramètre, à partir d’éprouvette saines sorties de production, chaque plage de conformité comprenant par exemple les 99 centiles inférieurs des valeurs du paramètre associé mesurée sur lesdites éprouvettes saines. Cette étape permet de détecter et quantifier les valeurs dites « non conformes » pour chaque paramètre et selon chaque direction de l’espace, et de déterminer un score partiel de conformité correspondant à la quantité de non-conformités détectées. Les scores partiels de conformité sont par exemple classés dans une matrice à trois lignes et trois colonnes, pour chacun des trois outils de traitement numérique et chacune des trois directions. Un score de conformité global de la sous-partie 31 est obtenu en sommant les scores partiels (par exemple les neufs scores partiels dans le cas d’une matrice à trois lignes et trois colonnes), et permet de caractériser l’état de conformité de la sous-partie 31, notamment du tissage des fibres 21. Ce score de conformité est comparé à une valeur de seuil de conformité prédéterminée, qui dépend du niveau d’exigences relatif à la pièce 20, afin de décider si la pièce 20 présente une ou plusieurs potentielles anomalies. Le procédé décrit, exécuté entièrement par un ordinateur connecté au système de mesure tomographique, permet d’obtenir un diagnostic en un temps rapide, diminué d’un facteur 5, et sans nécessiter une intervention d’un inspecteur qualifié. Si nécessaire, l’inspecteur peut ensuite réexaminer la pièce. Dans ce cas, les scores de conformité de chaque sous-partie 31 et les courbes d’évolution des paramètres caractéristiques permettent de localiser rapidement les régions pertinentes de la pièce 20 à inspecter, et facilitent cette observation supplémentaire. Procédé d’analyse tomographique d’une pièce (20), le procédé comprenant les étapes suivantes : - acquisition d’au moins une image tomographique (30) de la pièce (20) au moyen d’un dispositif de tomographie (1), - division de l’image tomographique (30) en sous-parties (31) élémentaires et obtention d’une pluralité d’images (32) bidimensionnelles pour chaque sous-partie (31), chaque image (32) comprenant une pluralité de pixels, chaque pixel présentant un niveau de gris respectif, - description des images (32) par application d’au moins un outil de traitement numérique sur chaque image (31) de la sous-partie (31) et obtention, pour chaque outil appliqué, d’une valeur d’au moins un paramètre représentatif pour chaque image (32), - détermination de variations spatiales dans chaque sous-partie (31), selon trois directions (X, Y, Z) caractéristiques, des valeurs de chaque paramètre et détection de valeurs des paramètres hors de plages de conformité prédéterminées, - calcul d’un score de conformité de chaque sous-partie (31), à partir des valeurs des paramètres hors des plages de conformité, et - comparaison du score de conformité de chaque sous-partie (31) à un seuil de conformité prédéterminé pour déterminer une conformité de la sous-partie (31). Procédé selon la revendication précédente, dans lequel le procédé comprend une étape de traitement des images (31) bidimensionnelles, mise en œuvre avant l’étape d’application d’au moins un outil de traitement numérique, au cours de laquelle chaque pixel dont le niveau de gris est hors d’une plage médiane prédéterminée est remplacé par un pixel corrigé présentant un niveau de gris compris dans ladite plage médiane. Procédé selon la revendication précédente, dans lequel le niveau de gris corrigé est déterminé aléatoirement à partir d’une distribution des niveaux de gris des pixels de l’image (31) dont le niveau de gris est dans la plage médiane. Procédé selon l’une des revendications précédentes, dans lequel l’au moins un outil de traitement numérique comprend un calcul d’autocorrélation, comprenant la génération d’une image d’autocorrélation (33) à partir de l’image (31), la détection d’un pic central (34) de l’image d’autocorrélation (33) et la détermination d’un rapport d’isotropie du pic central (34), le rapport d’isotropie étant choisi comme valeur de paramètre représentatif. Procédé selon l’une des revendications précédentes, dans lequel l’au moins un outil de traitement numérique comprend un calcul de filtrage, comprenant l’application à l’image (31) d’au moins un filtre de Gabor et l’obtention d’une image filtrée, la détermination, pour chaque pixel de l’image (31) d’une somme des valeurs correspondantes des images filtrées avec chacun des filtres de Gabor, le maximum sur l’image (31) des sommes étant choisi comme valeur de paramètre représentatif. Procédé selon la revendication précédente, dans lequel quatre filtres de Gabor sont appliqués, avec des bases respectives obtenues à partir des directions (X, Y, Z) caractéristiques de l’image (31). Procédé selon l’une des revendications précédentes, dans lequel l’au moins un outil de traitement numérique comprend un calcul d’un écart-type des niveaux de gris sur l’ensemble des pixels de l’image (31), ledit écart-type étant choisi comme valeur de paramètre représentatif. Procédé selon l’une des revendications précédentes, dans lequel trois outils de traitement numériques sont appliqués sur chaque image (31), l’étape de calcul du score de conformité comprenant le calcul d’un score partiel pour chaque outil de traitement numérique et pour chaque direction (X, Y, Z), à partir de la quantité de valeurs des paramètres hors des plages de conformité obtenue pour ledit outil de traitement numérique et ladite direction, le score de conformité étant la somme des scores partiels. Procédé selon l’une des revendications précédentes, dans lequel les plages de conformité des paramètres sont prédéterminées à partir d’éprouvette saines, chaque plage de conformité comprenant les 99 centiles inférieurs des valeurs du paramètre associé mesurée sur lesdites éprouvettes saines. Produit programme d’ordinateur comprenant des instructions permettant, lorsqu’il est exécuté sur un ordinateur pilotant un dispositif d’analyse tomographique (1), la mise en œuvre du procédé selon l’une des revendications précédentes.