FR 3117597 A1 20220617 FR 2013391 A 20201216 Procédé de contrôle non destructif d’une pièce, par exemple une pièce aéronautique, système, produit programme d’ordinateur, pièce aéronautique et aéronef, le procédé comprenant une étape d’acquisition de plusieurs images de la pièce à différents instants d’acquisition respectifs, chaque image comprenant des pixels associés respectivement à des zones de la pièce, une valeur du pixel étant représentative d’une caractéristique de la zone, une étape d'association de chaque zone à un vecteur ayant des composantes comprenant les valeurs du pixel aux différents instants d’acquisition, une étape de détermination de paramètres de référence, associés à au moins une zone, par moyennage du vecteur de zones voisines de la zone d’une pièce de référence et une étape de comparaison du vecteur de la zone avec les paramètres de référence pour détecter une anomalie sur la zone. Figure pour l’abrégé : Fig. 1 P r oc é d é de contrôle non destructif pour une pièce aéronautique Domaine technique général et art antérieur La présente invention concerne les procédés de détection d’anomalies dans des pièces industrielles. Ces procédés sont aussi connus sous l’expression de contrôle non-destructif (CND). L’invention peut être utilisée par exemple pour la détection d’anomalies dans des pièces industrielles de grande taille. Ces pièces industrielles peuvent être des pièces aéronautiques, comme par exemple des pièces de moteur d’aéronefs de type structures internes fixes (IFS pour Inner Fixed Structures en anglais). Le contrôle non destructif permet de caractériser l’état d’intégrité et la qualité de structure d’une pièce sans avoir à dégrader cette pièce. Pour réaliser ce contrôle non destructif il est connu d’utiliser un procédé par imagerie active. Un procédé par imagerie active s’appuie sur une excitation contrôlée de la pièce à analyser, afin de faire évoluer une grandeur physique de la pièce, et sur l’observation, au moyen d’un dispositif d’acquisition, de l’évolution temporelle de la grandeur physique. Cette observation de la pièce à analyser est réalisée pendant puis après l’excitation. Elle peut aussi commencer avant l’excitation. L’analyse de l’évolution temporelle de la grandeur physique est alors utilisée pour détecter une éventuelle anomalie dans la pièce. On connait par exemple la thermographie infrarouge active par excitation flash ou par induction dans laquelle l’excitation provoque une augmentation instantanée et éphémère de plusieurs degrés de la température localisée, par exemple à la surface ou sous la surface, de la pièce et dans laquelle une caméra infrarouge permet de déterminer l’évolution de la température surfacique de la pièce à la suite de cette excitation initiale. Les procédés connus, en particulier dans une phase d’apprentissage statistique, ne sont généralement pas parallélisables ce qui rend pratiquement le processus d’apprentissage irrémédiablement lent. En effet, dans les procédés connus, les éléments de donnée, typiquement des pixels ou des patchs de pixels, alimentant l’apprentissage, statistique sont parcourus séquentiellement et la complexité calculatoire des procédés connus augmente au moins proportionnellement à la surface des pièces à contrôler. Il y a donc un besoin pour un procédé de contrôle non destructif qui permette la détection d’une indication sur la présence ou l’absence d’une anomalie dans une pièce, qui soit parallélisable et qui permette un passage à l’échelle. Présentation générale de l’invention Dans ce cadre, au moyen de la présente invention, on propose de résoudre le problème de la détection d’anomalies dans des pièces industrielles. Il est ainsi proposé selon un aspect de l’invention un procédé de contrôle non destructif d’une pièce, par exemple une pièce aéronautique, le procédé comprenant une étape d’acquisition de plusieurs images de la pièce à différents instants d’acquisition respectifs, chaque image comprenant des pixels associés respectivement à des zones de la pièce, une valeur du pixel étant représentative d’une caractéristique de la zone, une étape d'association de chaque zone à un vecteur ayant des composantes comprenant les valeurs du pixel aux différents instants d’acquisition, une étape de détermination de paramètres de référence, associés à au moins une zone, par moyennage du vecteur de zones voisines de la zone d’une pièce de référence et une étape de comparaison du vecteur de la zone avec les paramètres de référence pour détecter une anomalie sur la zone. Le procédé de contrôle non destructif présenté ci-dessus permet donc la détection automatique des anomalies dans une pièce et ne nécessite pas nécessairement l’intervention d’un opérateur humain pour analyser toutes les données et déterminer si une anomalie est présente. De plus le traitement réalisé pour une zone de la pièce est indépendant des traitements réalisés pour les autres zones, et le nombre de traitements à réaliser pour analyser une pièce augmente linéairement avec le nombre de zones et donc la surface de la pièce. Donc le procédé est parallélisable et permet un passage à l’échelle avec des calculateurs à architecturés pour le calcul parallèle, à la fois pour le procédé de détermination des statistiques de référence et pour le procédé de comparaison. Le procédé peut aussi prendre en compte des informations en provenance d’un opérateur humain sur la pertinence des anomalies détectées. Sur la base de la pertinence des anomalies le procédé permet la mise à jour des paramètres statistiques de référence. Le procédé présente également les avantages suivants : Il permet de détecter des anomalies dans de grandes pièces industrielles. Il permet de détecter des anomalies résolues arbitrairement petites ou grandes. Il permet de détecter des anomalies dans des pièces dont la profondeur est hétérogène. Il permet de détecter des anomalies dans des pièces hétérogènes radialement. Il permet d’avoir une approche non-supervisée. En effet dans le procédé décrit il n’y a pas besoin de connaitre de manière exhaustive les configurations anormales. Le procédé de contrôle non destructif peut être mis en œuvre de la manière suivante. Dans un mode de réalisation le procédé de contrôle non destructif comprend en outre une étape d’excitation de la pièce, produisant une perturbation de la caractéristique. Dans un mode de réalisation le procédé de contrôle non destructif comprend en outre, lorsque l’anomalie n’a pas été détectée sur la zone une étape d’estimation des paramètres statistiques unitaires de la zone et une étape de mise à jour des paramètres statistiques de référence de la zone à partir des paramètres statistiques unitaires de la zone, lorsque l’anomalie a été détectée sur la zone une étape de validation d’une présence de l’anomalie, lorsque la présence de l’anomalie n’est pas validée une étape d’estimation des paramètres statistiques unitaires de la zone et une étape de mise à jour des paramètres statistiques de référence de la zone à partir des paramètres statistiques unitaires de la zone. Ce mode de réalisation permet d’avoir un fonctionnement dynamique du procédé. En particulier il permet de mettre à jour les paramètres de référence qui seront utilisés pour le contrôle des pièces. Cela est rendu possible par l’aspect incrémental de la mise à jour des paramètres de référence. Dans un mode de réalisation l’étape de comparaison comprend la détermination d’un test statistique d’une conformité du vecteur de la zone avec une distribution paramétrée par les paramètres de référence, le procédé de contrôle non destructif comprenant en outre une étape d’émission d’un message indiquant la détection d’une anomalie lorsqu’une valeur du test statistique est supérieure à un seuil. Dans un mode de réalisation la caractéristique est une température surfacique de la pièce à contrôler, l’étape d’excitation est réalisée par une lampe flash ou un générateur d’un champ magnétique pour provoquer une augmentation d’une température de la pièce et l’étape d’acquisition est réalisée par une caméra infrarouge et les composantes du vecteur de la zone représentent la température surfacique de la pièce. Dans un mode de réalisation le procédé de contrôle non destructif comprend en outre une étape de pré-traitement, l’étape de pré-traitement étant réalisée entre l’étape d’acquisition et l’étape de comparaison. Dans un mode de réalisation, les paramètres de référence comprennent une matrice de covariance de référence associée à la zone, et un vecteur moyen de référence, associé à la zone. On détermine la matrice de covariance de référence à partir d’une matrice de moments d’ordre 2 unitaire, la matrice de moments d’ordre 2 par récurrence à partir de paramètres statistiques unitaires, le vecteur moyen de référence par récurrence à partir des paramètres statistiques unitaires, les paramètres statistiques unitaires à partir du vecteur des zones voisines de la zone d’une pièce de référence, par convolution du vecteur des zones voisines de la zone d’une pièce de référence avec un noyau de convolution, la zone de la pièce de référence ne présentant pas d’anomalie. Un autre aspect de l’invention concerne un procédé de détermination de paramètres de référence, le procédé comprenant les étapes suivantes réalisées pour une pluralité de pièces de référence et dont une zone ne comporte pas d’anomalie : - excitation de la pièce de référence, produisant une perturbation d’une caractéristique de la pièce de référence, - acquisition de plusieurs images de la pièce de référence à différents instants d’acquisition respectifs, chaque image comprenant des pixels associés respectivement à des zones de la pièce, une valeur du pixel étant représentative d’une caractéristique de la zone, - association de chaque zone à un vecteur ayant des composantes comprenant les valeurs du pixel aux différents instants d’acquisition, - détermination, à partir des images, de paramètres statistiques unitaires de la zone et - détermination de paramètres de référence de la zone, obtenus par récurrence à partir des paramètres statistiques unitaires associés aux pièces de références. Ainsi la détermination des paramètres de référence est réalisée de manière incrémentale ce qui permet de limiter le temps de calcul nécessaire aux calculs des paramètres de référence. Dans un mode de réalisation l’étape de détermination comprend le moyennage des paramètres statistiques unitaires associés aux pièces de références. Dans un mode de réalisation le procédé de détermination comprend en outre une étape de pré-traitement, l’étape de pré-traitement étant réalisée entre l’étape d’acquisition et l’étape de détermination. Un autre aspect de cette invention concerne un système de contrôle non destructif d’une pièce, par exemple une pièce aéronautique. Le système de contrôle comprend un dispositif d’excitation, apte à perturber une caractéristique physique de la pièce, un dispositif d’acquisition d’une pluralité d’images de la pièce et une unité de traitement apte à la mise en œuvre du procédé de contrôle non destructif et/ou la mise en œuvre du procédé de détermination de paramètres de référence. Un autre aspect de cette invention concerne un produit programme d'ordinateur comportant des instructions logicielles qui, lorsqu’elles sont mises en œuvre par un équipement informatique, mettent en œuvre le procédé de contrôle non destructif et/ou le procédé de détermination de paramètres de référence. Un autre aspect de cette invention concerne une pièce aéronautique qui a été contrôlée par le procédé de contrôle non destructif. Un autre aspect de cette invention concerne un aéronef comprenant la pièce aéronautique. D’autres caractéristiques et avantages de l’invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des figures annexées : représente un système de contrôle non destructif selon un aspect de l’invention, et représentent de deux autres manières le système de contrôle non destructif selon un aspect de l’invention, représente un procédé de contrôle non destructif selon un aspect de l’invention, représente un procédé de détermination de paramètres statistiques de référence selon un aspect invention, représente un noyau gaussien utilisable pour le calcul des paramètres statistiques de cette invention. Description d’un ou plusieurs modes de réalisation La représente un système S de contrôle non destructif d’une pièce P à contrôler, par exemple une pièce industrielle de type pièce aéronautique. Le système S comprend un dispositif d’imagerie active 101. Ce dispositif d’imagerie active 101 comprend un dispositif d’excitation 101-a et un dispositif d’acquisition 101-b. Le dispositif d’excitation 101-a est apte à modifier une caractéristique physique de la pièce P. Le dispositif d’acquisition 101-b est apte à mesurer une caractéristique physique de la pièce P en différentes zones, à la suite de cette excitation. Le dispositif d’excitation 101-a est, par exemple, une lampe flash provoquant une perturbation éphémère de l’état thermique de la pièce P. Ainsi l’état thermique de la pièce est perturbé de manière transitoire, ce qui a pour effet de rendre le contrôle de la pièce non destructif. Le dispositif d’acquisition 101-b est, par exemple, une caméra infrarouge permettant de déterminer l’évolution de la température à la surface de la pièce. Lorsque la source d’excitation 101-a est une lampe flash et le dispositif d’acquisition 101-b est une caméra infrarouge, on parle alors de thermographie infrarouge active par excitation flash. La lampe flash est utilisée avantageusement avec des pièces en matériaux composites ou du moins non-métallique et permet pratiquement de contrôler la pièce à une profondeur supérieure à celle possible avec d’autres types de dispositifs d’excitation 101-a. On peut aussi utiliser la thermographie par excitation induction qui est avantageusement choisie pour le contrôle de pièces en matériaux métalliques pour détecter des défauts dans la matière à la surface ou proche de la surface de la pièce. Dans ce cas le dispositif d’excitation 101-a permet la génération d’un champ magnétique qui provoque par induction une modification d’état de la pièce, à savoir un réchauffement, en une zone située sous la surface de la pièce. Le dispositif d’acquisition 101-b est une caméra qui mesure la variation thermique en surface de la pièce. On peut aussi utiliser la radiographie. Dans ce cas le dispositif d’excitation 101-a comprend une source de rayons X et le dispositif d’acquisition 101-b comprend un écran-capteur. Les rayons X sont atténués par la quantité de matière traversée, comme dans le cas d’une radiographie médicale. Le dispositif d’acquisition 101-b fournit une séquence temporelle d’images numériques. Les images numériques comprennent un ensemble de pixels. Chaque pixel est représenté par une valeur. Chaque zone est associée à un pixel. La valeur de chaque pixel de chaque image est représentative de la caractéristique physique en une zone de la pièce et en l’instant d’acquisition de l’image. On associe à chaque zone un vecteur. Les composantes de ce vecteur comprennent les valeurs, dans les différentes images, du pixel associé à la zone. Dans le cas d’une image unique à valeurs scalaires, on peut appliquer des filtres texturels qui permettent d’obtenir des images d’indicateurs texturels. On peut associer un vecteur à chaque zone, les composantes du vecteur étant les valeurs du pixel associé à la zone pour les différentes images-textures obtenues par les différents filtrages. Dans le cas d’une vidéo, la dimension du vecteur correspond typiquement au temps échantillonné. Ainsi les composantes du vecteur d’une zone sont les valeurs des pixels des différentes images de la vidéo. Dans la suite par vecteur d’un pixel ou d’une zone on comprend le vecteur dont les composantes sont les valeurs prises par le pixel sur les différentes images acquises ou pre-traitées. Ainsi les images numériques représentent l’évolution temporelle de la caractéristique physique en différentes zones de la pièce. Dans le cas où le dispositif d’acquisition est une caméra infrarouge, chaque image numérique peut être une image codée en niveau de gris. Le niveau de gris de chaque pixel est représentatif de la température de la surface de la pièce en la zone associée au pixel à un instant déterminé de la trame. Il est préférable que la pièce et le dispositif d’acquisition 101-b soient fixes l’un par rapport à l’autre au moment de l’acquisition de la pluralité des images. Ainsi un pixel représentera la caractéristique physique d’une même zone de la pièce évoluant dans le temps. Il est aussi possible que la position relative de la pièce P et du dispositif d’acquisition 101-b varie de manière connue durant l’acquisition des images. Dans ce cas une étape de recalage des images est nécessaire, afin qu’un même pixel de différentes images représente la même zone de la pièce. Généralement l’enregistrement réalisé par le dispositif d’acquisition 101-b débute un peu avant l’instant de début de l’excitation et continue durant un temps approprié, permettant, si l’acquisition est suffisamment longue, le retour de la zone de la pièce excitée à l’état préalable à l’excitation. Le système comprend une unité de traitement 102 comprenant une entrée 102-a pour recevoir des images du dispositif d’acquisition 101-b et une sortie 102-b pour indiquer la détection anomalies. L’unité de traitement 102 est configurée pour la mise en œuvre d’un procédé de contrôle non destructif détaillé ci-après. L’unité de traitement 102 est également configurée pour la mise en œuvre d’un procédé de détermination de paramètres statistiques de référence détaillé ci-après. Ces paramètres statistiques de référence sont aussi appelés plus simplement paramètres de référence. Les paramètres statistiques de référence sont des paramètres d’un modèle de l’évolution temporelle de la caractéristique physique en une zone ne présentant pas d’anomalie. Les figures 2-a et 2-b, représentent deux modes de réalisation du système S de contrôle non destructif d’une pièce P. Sur ces figures le système utilise la thermographie infrarouge active par excitation flash. La figure 2-a comprend également un dispositif de synchronisation SYNC permettant la synchronisation de l’excitation avec la synchronisation de l’acquisition. La représente le procédé de contrôle non destructif de la pièce P. Ce procédé est mis en œuvre par l’unité de traitement 102 et permet la détection automatique d’anomalie dans la pièce P. Cette détection est réalisée en utilisant les images fournies par le dispositif d’acquisition 101-b. Le procédé de détection d’une anomalie de la comporte : - une étape d’excitation 301 de la pièce P, permettant une perturbation éphémère de la caractéristique physique de la pièce, - une étape d’acquisition 302 d’une pluralité d’images de la pièce P, les images étant acquises à une pluralité d’instants temporels, les images étant constituées d’un ensemble de pixels, les pixels étant associés respectivement à des zones de la pièce, on associe à chaque zone un vecteur dont les composantes représentent l’évolution temporelle de la caractéristique physique de la pièce en la zone associée, - une étape de comparaison 304 des vecteurs associés aux zones avec des paramètres statistiques de référence pour détecter une anomalie sur la zone. Dans un mode de réalisation l’étape d’excitation 301 est réalisée par une lampe flash ou un générateur d’un champ magnétique pour provoquer une augmentation d’une température, par exemple une température surfacique ou sub-surfacique de la pièce P et l’étape d’acquisition 302 est réalisée par une caméra infrarouge et les composantes du vecteur de la zone représentent la température surfacique de la pièce P. Ce procédé peut fonctionner dans deux modes : un mode statique et un mode dynamique. Le mode statique est celui dans lequel les paramètres statistiques de référence n’évoluent pas lors de l’analyse des pièces. Le mode dynamique est un mode, décrit en détail par la suite, dans lequel les paramètres statistiques de référence sont mis à jour au fur et à mesure que les pièces sont analysées. Les paramètres statistiques de référence comprennent, pour chaque zone, un vecteur moyen de référence, notée dans la suite m i pour la zone i. Ce vecteur moyen de référence est la moyenne locale des vecteurs des zones voisines au de la zone i, sans anomalie, estimée sur un ensemble de pièces équivalentes. Les paramètres statistiques de référence comprennent également, pour chaque zone, une matrice de covariance de référence, notée dans la suite Γ i pour la zone i. Ces paramètres statistiques de référence sont déterminés en utilisant un ensemble de pièces de référence équivalentes à la pièce à contrôler et sur lesquelles les zones ne comportant pas d’anomalie sont connues. Un tel procédé de détermination des paramètres statistiques de référence est présenté ultérieurement. L’étape de comparaison 304 se base sur la conformité entre le vecteur x i de la zone i et la distribution gaussienne attendue par ce vecteur, représentée par les paramètres statistiques de référence, lorsque la zone ne comprend pas d’anomalie. Pour déterminer cette conformité on peut utiliser la distance entre le vecteur x i d’une zone i et une distribution associée à ce vecteur, par exemple gaussienne, paramétrée par les paramètres statistiques de référence. Cette distance est déterminée au moyen de la distance de Mahalanobis calculée à partir du vecteur x i , du vecteur moyen de référence de la zone i et de la matrice de covariance de référence de la zone i. L’étape de comparaison 304 comprend, pour chaque zone i, la détermination d’une distance de Mahalanobis (notée dm(x i )) entre le vecteur de cette zone et la distribution gaussienne attendue et caractérisée par les paramètres statistiques de référence. Cette distance de Mahalanobis permet d’estimer si les données reçues suivent un modèle gaussien paramétré par le vecteur moyen de référence et la matrice de covariance de référence. Cette distance prend en compte la variabilité locale et inter-pièces via la matrice de covariance de référence. Elle est calulée en utilisant la formule suivante : Où, m i est le vecteur moyen de référence et Γ i est la matrice de covariance de référence. La distance de Mahalanobis est comprise en 0 et ∞ et est un indicateur d’anormalité. Dans un mode de réalisation, pour chaque zone la distance de Mahalanobis est comparée à un seuil prédéterminé. Si cette distance dépasse le seuil prédéterminé, une anomalie est détectée sur la zone. Dans un mode de réalisation le seuil est déterminé de manière empirique. Ce seuil empirique est un compromis entre le nombre de fausses alarmes et la probabilité de détection effective. Dans un mode de réalisation on utilise une approche quantitative avec le calcul de la fonction de survie et la possibilité de fixer un seuil qui correspond explicitement à la probabilité théorique acceptable de fausses alarmes. Enfin, dans un mode de réalisation on peut présenter à l’opérateur humain une carte de distances de Mahalanobis puis l’opérateur décide d’investiguer sur telle ou telle zone qui semble anormale car trop claire, dans ce cas il n’y a pas de comparaison avec le seuil prédéterminé. De plus le procédé peut comprendre une étape d’émission d’un message indiquant la détection d’une anomalie si la distance de Mahalanobis est supérieure au seuil prédéterminé. Optionnellement le procédé peut comporter une étape de prétraitement 303. Cette étape de prétraitement 303 est réalisée après l’étape d’acquisition 302 et avant l’étape de comparaison 304. Cette étape de prétraitement 303 comprend deux sous-étapes. Une première sous-étape de filtrage des images acquises en utilisant une pluralité de filtres permettant d’obtenir pour chaque image acquise une pluralité d’images filtrées. Chaque image filtrée ayant typiquement la même dimension spatiale que l’image acquise initiale. Une deuxième sous-étape de concaténation des valeurs des pixels des différentes images filtrées. Cela augmente la dimension du vecteur de la zone et, en particulier, si l’image d’entrée est à valeurs scalaires (en particulier si ce n’est pas une vidéo), on peut ainsi obtenir pour chaque zone un vecteur associé. Cette étape de prétraitement 303 peut aussi comprendre la transformation des images acquises selon un échantillonnage temporel d’acquisition et une durée d’acquisition fixée en des images dont les vecteurs des pixels sont les dérivées du logarithmique du vecteur des pixels par rapport au logarithme du temps passé depuis l’excitation, grandeur dite « dérivée logarithmique », évaluée selon un échantillonnage temporel dit « échantillonnage d’estimation ». Cette transformation est décrite dans la demande de brevet FR3071611A1. Cette étape de prétraitement 303 peut aussi comprendre le filtrage des images acquises en utilisant des filtres texturaux à vecteurs mono ou multidimensionnels. À la suite de ceci les images filtrées sont concaténées pour former une image à vecteurs qui constitue l’image prétraitée. Enfin, cette étape de prétraitement 303 peut aussi comprendre, dans le cas où la pièce comprend une absence de matière (trou), un filtre de détection des trous, afin d’obtenir une carte des trous, et ainsi une détermination des zones pour lesquelles la détection d’une anomalie n’est pas pertinente. La distance de Mahalanobis peut donc être considérée comme un score de détection et on peut optionnellement afficher à l’utilisateur pour chaque zone la valeur de cette distance. Cet affichage est connu sous l’expression de carte de scores ou carte de détection. Le procédé peut aussi comprendre une étape de post-traitement 306 qui prend en entrée la valeur des distances de Mahalanobis des zones calculées par l’étape de comparaison 304. Cette étape de post-traitement 306 réalise une régularisation de ces valeurs à partir d’informations de type métier, afin d’améliorer la détection des anomalies. Quelques exemples de post-traitements : - morphologie mathématique (érosion ou ouverture) afin de supprimer les anomalies inférieures à une certaine taille, - changement de dynamique, - application d’un seuil de décision sur chaque zone pour obtenir un masque binaire de détection, - visualisation, affichage, mosaïquage, juxtaposition, projection des cartes et masques de détection par projection sur une pièce simulée numériquement (par exemple une modélisation 3D de la pièce). La représente un procédé de détermination de paramètres statistiques de référence pour une zone ne comportant pas d’anomalie. Ce procédé utilise une pluralité de pièces de référence équivalentes à la pièce à contrôler pour lesquelles les zones ne comportant pas d’anomalie sont connues. Ces pièces de référence sont par exemple obtenues en analysant les zones des pièces et en déterminant si ces zones présentent ou non une anomalie. Ensuite seules les zones qui ne présentent pas d’anomalie sont utilisées pour déterminer les paramètres statistiques de référence. Les zones qui présentent une anomalie ne sont pas utilisées pour déterminer les paramètres statistiques de référence. Ce procédé comprend pour chaque pièce de référence : - l’excitation 301’ de la pièce de référence, permettant une perturbation éphémère de la caractéristique physique de la pièce de référence, - l’acquisition 302’ d’une pluralité d’images de la pièce de référence et - l’estimation 401, à partir des images, de paramètres statistiques unitaires du vecteur associé à la zone ne comportant pas d’anomalie. À la suite de l’obtention des paramètres statistiques unitaires pour l’ensemble des pièces de référence, le procédé de détermination de paramètres statistiques de référence comprend une étape de détermination ou mise à jour 402 des paramètres statistiques de référence. Cette étape de détermination 402 est réalisée par moyennage des paramètres statistiques unitaires. Ce moyennage permet d’obtenir un vecteur moyen de référence et une matrice de moments statistiques d’ordre 2 de référence, associées au vecteur de chaque zone. L’étape excitation 301’ et l’étape d’acquisition 302’ du procédé de détermination des paramètres statistiques de référence sont réalisées de la même manière que l’étape excitation 301 et l’étape d’acquisition 302 du procédé de contrôle non destructif. Les paramètres statistiques unitaires comprennent un vecteur moyen unitaire, notée dans la suite m i pour la zone i. Ce vecteur moyen unitaire est la moyenne des vecteurs des zones voisines d’une zone i sans anomalie. Les paramètres statistiques unitaires comprennent également une matrice des moments d’ordre 2 unitaire (notée dans la suite ) des vecteurs des zones voisines de la zone i sans anomalie. Ainsi l’étape d’estimation 401 comprend l’estimation du vecteur moyen unitaire, pour chaque zone, des vecteurs des zones voisines de cette zone. Ce vecteur moyen unitaire est estimé pour chaque zone sans anomalie. On peut utiliser la formule suivante pour réaliser l’estimation de cette moyenne unitaire : où est le vecteur moyen unitaire de la zone i pour la pièce d, J représente l’ensemble de décalage d’indice pour obtenir les zones voisines de la zone i, pour la pièce d, x d [j-i] représente le vecteur de la zone d’indice j-i, pour la pièce d, k[j] représente un facteur de pondération de la zone voisine j. Ce facteur de pondération peut être 1/N, où N est le nombre de zones voisines de la zone i. L’étape d’estimation 401 des paramètres statistiques unitaires comprend aussi l’estimation de la matrice des moments d’ordre 2 du vecteur des zones voisines de la zone i. Cette matrice des moments d’ordre 2 est estimée pour chaque zone sans anomalie. On peut utiliser la formule suivante pour réaliser l’estimation de cette matrice de covariance de la zone i pour la pièce d: où est une variable intermédiaire unitaire, dite matrice des moments statistiques d’ordre 2. Avantageusement, on calcule pour chaque nouvelle pièce de référence seulement et , puis la matrice de covariance est calculée à la fin des récurrences. L’étape d’estimation 402 des paramètres statistiques de référence comprend le moyennage des vecteurs moyens unitaires associés à la zone i obtenus pour l’ensemble des pièces de référence. Cette moyenne est obtenue en utilisant la formule suivante : Où est le vecteur moyen unitaire de la zone i pour la pièce d, D est le nombre de pièces de référence, m i est le vecteur moyen de référence pour la zone i. Le calcul des vecteurs moyens unitaires se fait par un produit de convolution spatial entre les vecteurs de la zone de la pièce de référence et un noyau de convolution. Avantageusement, ce noyau de convolution est une imagette à valeurs scalaires dont les valeurs sont les facteurs de pondérations k[ ] associés aux zones voisines. Avantageusement le noyau de convolution k représente une fonction positive, normée, centrée, compacte et symétrique. Ce qui veut dire que: Positive: ∀i,k[ ]≥0 Normée : ∑_i1,i2 k[ ]=1 Centrée, compacte : il existe deux entiers positifs J1 et J2 tels que ∀j vérifiant |j1|>J1 et |j2|>J2, k[j1,j2]=0 Symétrique : k[j1,j2]=k[-j1,-j2] Pour alléger l’écriture, on utilise une notation 1D qui se généralise naturellement à la 2D : le pixel d’indice i a pour coordonnées ; il peut être déterminé par son indice ou par ses coordonnées et en ce sens, k[ ] s’écrit de manière équivalente k[i]. De même, l’étape d’estimation 402 des paramètres statistiques de référence comprend le moyennage des variables intermédiaires unitaires associées à la zone i et obtenues pour l’ensemble des pièces de référence. Cette moyenne est obtenue en utilisant la formule suivante D est le nombre de pièces de référence, q i est une variable intermédiaire unitaire de référence pour la zone i. Ensuite la matrice de covariance de référence est obtenue en utilisant la formule : L’étape d’estimation 402 a l’effet technique de pouvoir prendre en compte la variabilité dans la production des pièces, par le moyennage des paramètres statistiques sur plusieurs pièces connues comme ne présentant pas d’anomalie. On peut aussi réaliser le moyennage des paramètres statistiques estimés, pour obtenir les paramètres statistiques de référence, en utilisant les formules de récurrence suivantes : Et Où est le vecteur moyen de référence pour la zone i calculé en utilisant les D premières pièces, est le vecteur moyen de référence pour la zone i calculé en utilisant les D-1 premières pièces, est le vecteur moyen de la zone i pour la pièce D, est la variable intermédiaire de référence pour la zone i calculée en utilisant les D premières pièces, est la variable intermédiaire de référence pour la zone i calculée en utilisant les D-1 premières pièces, est la variable intermédiaire unitaire de la zone i pour la pièce D. A la suite de la mise à jour de la variable intermédiaire de référence, la matrice de covariance de référence pour la zone i calculée en utilisant D+1 pièces est déterminée en utilisant la formule : Avantageusement, cette dernière formule, qui n’est pas une formule de récurrence, n’est appliquée que lorsque toutes les pièces de références sont passées dans le processus de récurrence. On calcule alors la matrice de covariance et son inverse. C’est son inverse qui sera utilisée pour la détection d’anomalies dans de nouvelles images, dans la formule de la distance de Mahalanobis. Ces formules de récurrence ont pour effet technique de pouvoir mettre à jour les vecteurs moyens de référence et les matrices de covariance sur la base des vecteurs moyens de référence et des matrices des moments d’ordre 2 de référence estimées précédemment et sur la base des vecteurs moyens unitaires et des matrices de moments d’ordre 2 unitaires estimées sur la dernière pièce analysée. Elles offrent donc l’effet technique de ne pas avoir à sauvegarder l’ensemble des images associées à l’ensemble des pièces de référence et de ne pas avoir à reparcourir toutes les données si l’on souhaite mettre à jour les paramètres de référence avec de nouvelles données ce qui peut représenter un gain calculatoire conséquent. En effet il est nécessaire de ne sauvegarder que les derniers paramètres statistiques de référence, le nombre D d’appels de la formule de récurrence ayant servi à les déterminer et les images associées à la dernière pièce. Ainsi les procédés de détermination des paramètres statistiques de référence et de mise à jour de ces paramètres sont incrémentaux. Les formules de récurrence précédentes affectent à chaque acquisition d’images le même poids, ce qui est avantageux dans le cas d’utilisation où la seule variabilité provient de la production et du processus d’acquisition. On peut prendre en compte un vieillissement progressif du système d’acquisition, en modifiant les formules de récurrence précédentes pour oublier progressivement les contributions passées afin de choisir une distribution de poids non-uniforme privilégiant les observations récentes. Cette propriété est atteinte en considérant les moyennes mobiles exponentielles des statistiques instantanées : Où τ représente une constante de temps comprise entre 0 et 1 permettant de régler l’importance relative d’une nouvelle donnée dans la mise à jour et donc la capacité d’évaluer progressivement la contribution des données passées. Si l’on considère que les N dernières pièces apportent une information significative, on prend typiquement τ=2/(N+1). Ainsi le procédé de détermination de paramètres statistiques de référence permet d’obtenir pour chaque zone i un vecteur moyen de référence mi et une matrice des moments statistiques d’ordre 2 de référence qi. Cette matrice qi représente une matrice de paramètres intermédiaire utilisée pour le calcul de la matrice de covariance Γi. et est un paramètre de la distribution supposée gaussienne de la zone. L’estimation des paramètres de référence mi et qi associés à la zone i est réalisée au moyen d’un moyennage local calculé au moyen d’un produit de convolution de l’image. Le produit de convolution est caractérisé par un noyau de convolution à plusieurs dimensions (autant que le nombre de dimensions spatiales de l’image, à savoir 2 dans notre cas). Le noyau de convolution peut avantageusement être un noyau gaussien centré de paramètre σ, exprimé en fonction de la taille du voisinage, qui est typiquement le maximum entre l’incertitude de position δ et la taille du voisinage d’une zone dans lequel la caractéristique physique est homogène. Un tel noyau gaussien est illustré en . Un tel produit de convolution est connu de l’homme du métier. Dans le cas plus général d’un noyau issu d’une famille non-gaussienne, le paramètre σ représente le vecteur de paramètres définissant le noyau à partir de sa famille d’appartenance. Il est avantageux de choisir un noyau positif, normé à support compact (typiquement dans le cas d’une image à deux dimensions : définit la taille du voisinage et vaut entre quelques zones et quelques dizaines de zones). Le procédé de contrôle non destructif représenté peut aussi comporter optionnellement des étapes visant à la mise à jour des paramètres statistiques de référence à partir de nouvelles données contrôlées par le procédé. Dans ce cas le procédé de contrôle non destructif comprend pour chaque zone : - si le vecteur de la zone est statistiquement anormal (distance de Mahalanobis supérieure à un seuil), une étape de validation 305 pour confirmer ou infirmer le label d’anomalie statistique; typiquement, si l’anomalie statistique correspond bien à un défaut elle est confirmée, sinon, elle est infirmée. - et si la détection de l’anomalie est infirmée une étape de mise à jour 306 des paramètres statistiques de référence, en utilisant le vecteur de la zone. L’étape de validation 305 est réalisée par un opérateur humain. Celui-ci utilise son expertise pour déterminer si la zone qui a été détectée comme comportant une anomalie présente réellement une anomalie. Typiquement, l’opérateur humain analyse le signal brut ou prétraité par l’étape 303 ou prétraitée de manière spécifique à la visualisation : il examine qualitativement ou quantitativement l’allure (dimension, contraste local, signature spécifique…) du signal dans la zone suspecte en appliquant une procédure prévue à cet effet. Si la détection de l’anomalie n’est pas correcte, il est avantageux d’utiliser le vecteur de cette zone pour mettre à jour le vecteur moyen de référence et la matrice de covariance de référence, car c’est une zone qui n’est pas caractérisée correctement par les paramètres statistiques de référence. Il est aussi possible de ne pas réaliser l’étape de validation 305 si l’on a la certitude qu’aucune anomalie ne se trouve dans les données. Dans ce cas l’étape de mise à jour 306 est réalisée pour l’ensemble des zones. Cette étape de mise à jour 306 comprend la détermination de paramètres statistiques unitaires du vecteur de la zone. Cette détermination est réalisée de manière identique à l’étape 401 de détermination de paramètres de statistiques unitaires du procédé de détermination de paramètres statistiques de référence de la . Ceci permet l’estimation d’un vecteur moyen unitaire, associée à la zone. Ce vecteur moyen est estimé pour chaque zone. On peut utiliser la formule suivante pour réaliser l’estimation de ce vecteur moyen unitaire : où est le vecteur moyen unitaire associé à la zone i, J représente les décalages d’indice pour obtenir l’ensemble des zones voisines de la zone i x[j-i] représente le vecteur de la zone d’indice j-i k[j] représente un facteur de pondération, qui peut de manière avantageuse être 1/N où N est le nombre de zones voisines de la zone i. L’étape de mise à jour 306 des paramètres statistiques de référence comprend l’estimation d’une matrice de moments d’ordre 2 unitaire associée à la distribution des vecteurs des zones voisines: On peut utiliser la formule suivante pour réaliser l’estimation de la matrice de covariance : Cependant et avantageusement, n’est pas calculée si l’on souhaite ne mettre à jour la matrice de covariance qu’après le passage de plusieurs pièces, les seules statistiques unitaires calculées sont et . Ensuite l’étape de mise à jour 306 des paramètres statistiques de référence utilise les formules de récurrence décrites précédemment afin de mettre à jour les paramètres statistiques de référence. Ces formules de récurrences sont : Et Où est le vecteur moyen de référence pour la zone i calculée en utilisant D premières pièces, est le vecteur moyen de référence pour la zone i calculée en utilisant D-1 premières pièces, est le vecteur moyen de la zone i pour la pièce D, est la matrice des moments statistiques d’ordre 2 de référence pour la zone i calculée en utilisant D premières pièces, est la matrice des moments statistiques d’ordre 2 de référence pour la zone i calculée en utilisant D-1 premières pièces et est la matrice des moments statistiques d’ordre 2 unitaire de la zone i pour la pièce D. A la suite de la mise à jour de la variable intermédiaire de référence, la matrice de covariance de référence pour la zone i est calculée en utilisant la formule suivante : On peut aussi utiliser les formules de récurrence suivantes : Où τ représente une constante de temps comprise entre 0 et 1 permettant de régler l’importance relative d’une nouvelle donnée dans la mise à jour et donc la capacité d’évaluer progressivement la contribution des données passées. Cette mise à jour permet d’ajuster le procédé de contrôle non destructif à des changements locaux dans les pièces sur lesquelles le contrôle est réalisé. Procédé de contrôle non destructif d’une pièce (P), par exemple une pièce aéronautique, le procédé comprenant : - une étape d’acquisition (302) de plusieurs images de la pièce (P) à différents instants d’acquisition respectifs, chaque image comprenant des pixels associés respectivement à des zones de la pièce, une valeur du pixel étant représentative d’une caractéristique de la zone, - une étape d'association de chaque zone à un vecteur ayant des composantes comprenant les valeurs du pixel aux différents instants d’acquisition, - une étape de détermination (402) de paramètres de référence, associés à au moins une zone, par moyennage du vecteur de zones voisines de la zone d’une pièce de référence et - une étape de comparaison (304) du vecteur de la zone avec les paramètres de référence pour détecter une anomalie sur la zone. Procédé de contrôle non destructif selon la revendication 1 comprenant en outre : - une étape d’excitation (301) de la pièce (P), produisant une perturbation de la caractéristique. Procédé de contrôle non destructif selon la revendication 1 ou 2, comprenant en outre, lorsque l’anomalie n’a pas été détectée sur la zone : - une étape d’estimation (401) de paramètres statistiques unitaires de la zone et - une étape de mise à jour (306) des paramètres de référence de la zone à partir des paramètres statistiques unitaires de la zone, lorsque l’anomalie a été détectée sur la zone : - une étape de validation (305) d’une présence de l’anomalie, et lorsque la présence de l’anomalie n’est pas validée : - une étape d’estimation (401) des paramètres statistiques unitaires de la zone et - une étape de mise à jour (306) des paramètres de référence de la zone à partir des paramètres statistiques unitaires de la zone. Procédé de contrôle non destructif, selon l’une des revendications précédentes, dans lequel : - l’étape de comparaison (304) comprend la détermination d’une valeur d’un test statistique représentatif d’une conformité du vecteur de la zone avec une distribution paramétrée par les paramètres de référence, le procédé de contrôle non destructif comprenant en outre : - une étape d’émission d’un message indiquant la détection d’une anomalie lorsque la valeur du test statistique est supérieure à un seuil. Procédé de contrôle non destructif, selon l’une des revendications précédentes, dans lequel : - la caractéristique est une température surfacique de la pièce (P), - l’étape d’excitation (301) est réalisée par une lampe flash ou un générateur d’un champ magnétique pour provoquer une augmentation d’une température de la pièce (P) et - l’étape d’acquisition (302) est réalisée par une caméra infrarouge et les composantes du vecteur de la zone représentent la température surfacique de la pièce (P). Procédé de contrôle non destructif, selon l’une des revendications précédentes, comprenant en outre une étape de pré-traitement (303), l’étape de pré-traitement (303) étant réalisée entre l’étape d’acquisition (302) et l’étape de comparaison (304). Procédé de contrôle non destructif, selon l’une des revendications précédentes, les paramètres de référence comprenant : - une matrice de covariance de référence associée à la zone, et - un vecteur moyen de référence, associé à la zone, on détermine : - la matrice de covariance de référence à partir d’une matrice de moments d’ordre 2, - la matrice de moments d’ordre 2 par récurrence à partir de paramètres statistiques unitaires, - le vecteur moyen de référence par récurrence à partir des paramètres statistiques unitaires et - les paramètres statistiques unitaires à partir du vecteur des zones voisines de la zone d’une pièce de référence, par convolution du vecteur des zones voisines de la zone de la pièce de référence avec un noyau de convolution, la zone de la pièce de référence ne présentant pas d’anomalie. Procédé de détermination de paramètres de référence utilisés dans le procédé de contrôle non destructif d’une pièce (P) de la revendication 1, le procédé comprenant les étapes suivantes, réalisées pour une pluralité de pièces de référence et dont une zone ne comporte pas d’anomalie : - excitation (301’) de la pièce de référence, produisant une perturbation d’une caractéristique de la pièce de référence, - acquisition (302’) de plusieurs images de la pièce de référence à différents instants d’acquisition respectifs, chaque image comprenant des pixels associés respectivement à des zones de la pièce, une valeur du pixel étant représentative d’une caractéristique de la zone, - association de chaque zone à un vecteur ayant des composantes comprenant les valeurs du pixel aux différents instants d’acquisition, - détermination (401), à partir des images, de paramètres statistiques unitaires de la zone et - détermination (402) de paramètres de référence de la zone, obtenus par récurrence à partir des paramètres statistiques unitaires associés aux pièces de référence. Procédé de détermination de paramètres de référence, selon la revendication 8, l’étape de détermination (402) comprend le moyennage des paramètres statistiques unitaires associés aux pièces de référence. Procédé de détermination de paramètres de référence, selon la revendication 8 ou 9, comprenant en outre une étape de pré-traitement (303), l’étape de pré-traitement (303) étant réalisée entre l’étape d’acquisition (302’) et l’étape de détermination (401). Système (S) de contrôle non destructif d’une pièce (P), par exemple une pièce aéronautique, le système (S) de contrôle non destructif comprenant : - un dispositif d’excitation (101-a), apte à modifier une caractéristique physique de la pièce (P), - un dispositif d’acquisition (101-b) d’une pluralité d’images de la pièce (P) et - une unité de traitement (102) apte à la mise en œuvre du procédé de contrôle non destructif selon l’une des revendications 1 à 7 et/ou la mise en œuvre du procédé de détermination de paramètres de référence selon l’une des revendications 8 à 10. Produit programme d'ordinateur comportant des instructions logicielles qui, lorsqu’elles sont mises en œuvre par un équipement informatique, mettent en œuvre le procédé de contrôle non destructif selon l’une des revendications 1 à 7 et/ou le procédé de détermination de paramètres de référence selon l’une des revendications 8 à 10. Pièce aéronautique qui a été contrôlée par le procédé de contrôle non destructif selon l’une des revendications 1 à 7. Aéronef comprenant la pièce aéronautique selon la revendication 13.