FR 3071611 A1 20190329 FR 1759022 A 20170928 Arrière-plan de l'invention L'invention concerne le domaine général de l'aéronautique. Elle se rapporte plus particulièrement au contrôle non destructif par thermographie infrarouge active de pièces aéronautiques, telles que des pièces à géométrie complexe équipant des moteurs d'aéronef comme par exemple des structures internes fixes (ou IFS pour Inner Fixed Structures) ou autres structures sandwiches composites d'inverseurs de poussée. Aucune limitation n'est attachée toutefois au type de pièce aéronautique considérée ni au matériau dans laquelle cette pièce est composée, il peut s'agir indifféremment d'un matériau composite ou non. De façon connue, on désigne par contrôle non destructif un ensemble de méthodes qui permettent de caractériser l'état d'intégrité et/ou la qualité de structures ou de matériaux sans les dégrader. Le contrôle non destructif a une application privilégiée mais non limitative dans le domaine de l'aéronautique, et plus généralement dans n'importe quel domaine dans lequel les structures dont on souhaite caractériser l'état ou la qualité sont coûteuses et/ou leur fiabilité de fonctionnement est critique. Le contrôle non destructif peut être avantageusement réalisé sur la structure ou le matériau considéré aussi bien en cours de production qu'en cours d'utilisation ou de maintenance. La thermographie active, ou thermographie infrarouge active, est une des techniques connues utilisées pour réaliser des contrôles non destructifs dans le domaine de l'aéronautique. Elle s'appuie en premier lieu sur une excitation contrôlée de la pièce considérée. Cette excitation conduit à un changement d'état thermique de la pièce (par exemple à une augmentation instantanée de la température surfacique de la pièce de plusieurs degrés dans le cas où la source d'excitation est une lampe flash). Puis dans un second temps, on observe au moyen d'une caméra infrarouge, l'évolution de la température à la surface de la pièce (par exemple la décroissance progressive de la température à la surface de la pièce dans le cas d'une excitation par lampe flash). La caméra infrarouge fournit une séquence temporelle d'images thermiques numériques représentant à différents instants la température en différents points de la surface de la pièce. Par exemple, chaque image numérique est une image codée en niveaux de gris sur laquelle chaque pixel associé à un point de la pièce a un niveau de gris représentatif d'une fonction croissante de la 2 température en ce point à l'instant de l'acquisition de l'image numérique. L'observation de la séquence temporelle d'images thermiques numériques fournies par la caméra infrarouge permet de détecter la présence d'anomalies dans l'évolution de la température à la surface de la pièce qui 5 seraient causées par des défauts inhérents à la pièce (ex. vides, inclusions, etc.) et qui perturberaient la diffusion de la chaleur de la surface vers l'intérieur de la pièce. Le document WO 2001/041421 décrit une méthode appelée méthode TSR (Thermographic Signal Reconstruction), permettant 10 d'exploiter les informations contenues dans les images thermiques fournies par la caméra infrarouge pour détecter des défauts affectant une pièce de matériau. Cette méthode s'appuie sur le constat que la réponse thermique en tout point d'une pièce de matériau décroit en fonction du temps suite à son excitation jusqu'à un état d'équilibre correspondant 15 sensiblement à l'état thermique de la pièce avant son excitation, et que, pour une excitation impulsionnelle (flash), cette décroissance est telle que le logarithme naturel de la perturbation en température (différence entre la température surfacique courante et la température à l'équilibre) varie de façon linéaire en fonction du logarithme naturel du temps écoulé après 20 l'excitation impulsionnelle dans le cas d'un matériau homogène d'épaisseur semi-infinie. Partant de ce constat, l'état de la technique propose un procédé permettant d'estimer les dérivées d'ordre 1 et/ou d'ordre 2 du logarithme naturel de la température, et d'utiliser les dérivées ainsi estimées pour déterminer la présence ou non d'un défaut dans la pièce de 25 matériau considérée. Pour estimer les dérivées susmentionnées, l'état de la technique convertit tout d'abord, au moyen d'une fonction logarithme naturel, les écarts d'amplitudes des pixels (par rapport à l'amplitude à l'équilibre évaluée à partir des images acquises avant l'excitation) des images 30 numériques fournies par la caméra infrarouge ainsi que les instants auxquels ces images numériques ont été acquises, l'origine du temps étant attribuée à l'instant de l'excitation. Puis, pour chaque pixel, il détermine par une régression au sens des moindres carrées effectuée sur l'ensemble des images numériques fournies par la caméra infrarouge, un 35 polynôme d'ordre n approximant l'évolution du logarithme de l'écart d'amplitude de ce pixel en fonction du logarithme du temps. Le polynôme 3 ainsi déterminé est alors dérivé pour estimer les dérivées d'ordre 1 et/ou d'ordre 2 recherchées. Malheureusement, la méthode proposée dans l'état de la technique est peu adaptée en pratique au contrôle non destructif de 5 grandes pièces souvent hétérogènes et/ou à géométrie complexe telles qu'elles existent couramment dans l'aéronautique, comme par exemple les structures IFS mentionnées précédemment. Le contrôle non destructif complet d'une pièce de grande dimension (plusieurs mètres carrés typiquement pour une structure IFS 10 par exemple) nécessite en effet de mener plusieurs acquisitions « unitaires » de séquences temporelles d'images thermiques pour couvrir l'ensemble de la surface de celle-ci, chaque acquisition unitaire portant sur une zone d'acquisition dite unitaire différente de la pièce. La juxtaposition des acquisitions unitaires permet ainsi d'obtenir une représentation de la 15 pièce dans son ensemble, moyennant la prise en compte d'une zone de recouvrement éventuelle entre chaque acquisition unitaire. La figure 1 illustre à titre d'exemple une pièce aéronautique telle qu'une structure IFS décomposée en une pluralité de zones d'acquisition unitaires rectangulaires, chacune étant destinée à faire l'objet d'une acquisition 20 unitaire. Les pièces aéronautiques peuvent en outre être hétérogènes, dans le sens où les propriétés thermiques du matériau qui les composent peuvent varier spatialement notamment avec la présence de défauts, et peuvent avoir des géométries complexes, c'est-à-dire que : 25 - l'épaisseur de la pièce peut varier en fonction de la position considérée à la surface de la pièce. En particulier, l'épaisseur maximale d'une zone d'acquisition unitaire peut varier d'une zone d'acquisition unitaire à l'autre ; et/ou - la surface de la pièce n'est pas nécessairement plane. 30 Plusieurs paramètres d'acquisition des images thermiques de telles pièces peuvent ainsi être amenés à varier, notamment d'une zone d'acquisition unitaire à l'autre lorsque plusieurs acquisitions unitaires sont nécessaires pour recouvrir une grande pièce dans son ensemble. Ces paramètres d'acquisition susceptibles de varier sont notamment : 35 - l'amplitude de l'excitation (i.e. énergie déposée) au niveau de la surface de la pièce, en raison par exemple de la non-uniformité de la 4 source d'excitation, de la non-uniformité de l'aspect de surface de la pièce, de sa géométrie complexe, ou de la distance variable entre la source d'excitation et la zone de la pièce considérée ; - le temps d'acquisition de chaque séquence temporelle d'images 5 thermiques de la pièce et l'échantillonnage (temporel, et voire spatial) appliqué (notamment le pas d'échantillonnage), qui peuvent varier en fonction de l'épaisseur de la zone considérée lors de l'acquisition. Le contrôle non destructif d'une pièce de grande dimension, hétérogène et/ou de géométrie complexe, doit donc s'appuyer, pour être 10 pertinent, sur des caractéristiques qui sont « raisonnablement » invariantes par rapport aux conditions d'acquisition des images thermiques sur la pièce dans son ensemble. Or, dans l'état de la technique, les dérivées d'ordre 1 et d'ordre 2 utilisées pour le contrôle non destructif sont estimées à partir d'une 15 approximation polynomiale déterminée sur toute la séquence temporelle d'images thermiques acquise. Le modèle polynomial utilisé est donc partagé par tous les échantillons temporels acquis pour un pixel. Il en résulte un biais d'estimation qui dépend, par construction, des conditions d'acquisition de la séquence d'images utilisée et notamment de la durée 20 d'acquisition de la séquence d'images. Ainsi, si les conditions d'acquisition sont modifiées entre deux zones d'acquisition unitaire adjacentes, ce qui peut s'avérer nécessaire lorsque ces zones ont des épaisseurs différentes, l'état de la technique peut conduire à des polynômes distincts et incidemment à des dérivées estimées distinctes pour une même structure 25 de matériau vue via deux acquisitions unitaires distinctes. Ceci peut mener à des résultats difficiles à interpréter en matière de contrôle non destructif (i.e. interprétation des dérivées obtenues). Cette difficulté d'interprétation est d'autant plus dommageable lorsqu'on envisage d'automatiser l'interprétation des dérivées obtenues 30 lors du contrôle non destructif de la pièce. Objet et résumé de l'invention La présente invention a pour but principal de pallier les inconvénients précités et de proposer une méthode pouvant être 35 appliquée à tout type de pièce notamment aéronautique, y compris aux 5 pièces aéronautiques de grande dimension, hétérogène et/ou de géométrie complexe. Elle propose à cet effet un procédé de contrôle non destructif d'une pièce aéronautique comprenant : 5 - une étape d'obtention d'une pluralité d'images numériques de la pièce aéronautique acquises à une pluralité d'instants d'acquisition définis sur une période de temps déterminée dite période d'acquisition au moyen d'un système de thermographie infrarouge active, chaque pixel d'une image numérique acquise à un instant d'acquisition ayant une 10 amplitude représentant une caractéristique de la pièce aéronautique à cet instant d'acquisition en un point de la pièce aéronautique ; - une étape d'estimation, à partir des images numériques acquises, d'au moins une dérivée d'une première fonction de ladite caractéristique par rapport à une deuxième fonction du temps ; et 15 - une étape de vérification de la pièce aéronautique utilisant ladite au moins une dérivée estimée. Ce procédé est remarquable en ce que l'étape d'estimation comprend, pour chaque dérivée estimée : - une étape de sélection d'une pluralité d'instants d'estimation de cette 20 dérivée sur ladite période d'acquisition ; - une étape d'application de la première fonction aux amplitudes des pixels des images numériques, ladite étape d'application résultant en des amplitudes de pixels converties ; - une étape d'application de la deuxième fonction aux instants 25 d'acquisition des images numériques résultant en des instants d'acquisition convertis et aux instants d'estimation de la dérivée résultant en des instants d'estimation convertis, lesdits instants d'estimation étant sélectionnés lors de l'étape de sélection de sorte que les instants d'estimation convertis soient uniformément espacés ; et 30 - une étape d'estimation, pour chaque instant d'estimation converti, de ladite dérivée à partir des amplitudes converties correspondant aux instants d'acquisition convertis appartenant à un intervalle d'estimation déterminé autour de l'instant d'estimation converti. Corrélativement, l'invention propose également un système de 35 contrôle non destructif d'une pièce aéronautique comprenant : 6 - un système de thermographie active comprenant un dispositif d'excitation et un dispositif d'acquisition configuré pour acquérir à une pluralité d'instants d'acquisition définis sur une période de temps déterminée dite période d'acquisition, des images numériques de la 5 pièce aéronautique, chaque pixel d'une image numérique acquise à un instant d'acquisition ayant une amplitude représentant une caractéristique de la pièce aéronautique à cet instant d'acquisition en un point de la pièce aéronautique après excitation de la pièce aéronautique par le dispositif d'excitation ; 10 - un dispositif d'estimation, configuré pour estimer à partir des images numériques acquises, au moins une dérivée d'une première fonction de ladite caractéristique par rapport à une deuxième fonction du temps ; et - un dispositif de vérification configuré pour vérifier la pièce 15 aéronautique en utilisant ladite au moins une dérivée estimée. Ce système est remarquable en ce que le dispositif d'estimation comprend une pluralité de modules activés pour chaque dérivée estimée et comprenant : - un module de sélection, configuré pour sélectionner une pluralité 20 d'instants d'estimation de cette dérivée sur ladite période d'acquisition ; - un premier module d'application, configuré pour appliquer la première fonction aux amplitudes des pixels des images numériques, ledit premier module d'application fournissant des amplitudes de pixels 25 converties ; - un deuxième module d'application, configuré pour appliquer la deuxième fonction aux instants d'acquisition des images numériques fournissant des instants d'acquisition convertis et pour appliquer la deuxième fonction aux instants d'estimation de la dérivée fournissant 30 des instants d'estimation convertis, lesdits instants d'estimation ayant étant sélectionnés par le module de sélection de sorte que les instants d'estimation convertis soient uniformément espacés ; et un module d'estimation, configuré pour estimer pour chaque instant d'estimation converti, ladite dérivée à partir des amplitudes converties 35 correspondant aux instants d'acquisition convertis appartenant à un 7 intervalle d'estimation déterminé autour de l'instant d'estimation converti. Pour une image numérique codée en niveaux de gris, chaque pixel est associé à un niveau de gris dont l'amplitude représente la 5 luminance spectrale vue par le capteur du dispositif d'acquisition du système de thermographie active et qui est, de façon connue en soi, une fonction croissante de la température de la pièce aéronautique au point correspondant au pixel considéré de la pièce aéronautique. Ainsi, l'invention propose un contrôle non destructif basé sur 10 une estimation de dérivée(s) temporelle(s) des amplitudes converties du signal acquis par thermographie infrarouge active robuste aux conditions d'acquisition des images thermiques de la pièce aéronautique. Cette robustesse est assurée par une estimation locale des dérivées (par opposition à une estimation globale de celles-ci sur l'ensemble du signal 15 acquis comme dans l'état de la technique), pour différents instants d'estimation choisis judicieusement préalablement et sur des intervalles d'estimation déterminés autour de chaque instant d'estimation pour estimer localement les dérivées. Chaque intervalle d'estimation associé à un instant d'estimation auquel on souhaite estimer la dérivée contient des 20 échantillons mesurés temporellement voisins de cet instant d'estimation. Cet intervalle d'estimation peut varier d'un instant d'estimation à l'autre. Si l'on considère que la période d'acquisition est délimitée par un instant dit initial et un instant dit final, on choisit préférentiellement pour chaque intervalle d'estimation déterminé autour d'un instant 25 d'estimation converti une taille inférieure à la différence entre l'instant final converti résultant de l'application de la deuxième fonction à l'instant final et l'instant initial converti résultant de l'application de la deuxième fonction à l'instant initial. En d'autres mots, chaque intervalle d'estimation déterminé conformément à l'invention définit avantageusement un 30 voisinage limité d'échantillons autour de l'instant d'estimation considéré pour calculer localement les dérivées. On note que conformément à l'invention, les instants d'estimation sont choisis de sorte que les instants d'estimation convertis soient uniformément espacés. Ceci permet de s'adapter à la base 35 temporelle dans laquelle on représente les dérivées (en particulier pour l'exploitation visuelle de ces dérivées qui peut être mise en oeuvre par 8 exemple lors de l'étape de vérification). En outre, les instants d'estimation peuvent être choisis de façon indépendante des instants d'acquisition (et donc différents des instants d'acquisition), et il est possible de façon avantageuse de choisir un nombre d'instants d'estimation de la dérivée 5 inférieur (voire bien inférieur) au nombre d'instants d'acquisition définis sur la période d'acquisition. En effet, pour caractériser les structures de matériau considérées, il est nécessaire de choisir un échantillonnage d'acquisition suffisamment dense (i.e. des instants d'acquisition suffisamment 10 rapprochés). Cet échantillonnage est généralement constant dans le domaine temporel et le passage dans un domaine temporel converti (i.e. via l'application de la deuxième fonction aux instants d'acquisition) peut engendrer, suivant la nature de la deuxième fonction appliquée, un sur-échantillonnage dans cet espace de représentation. 15 En considérant un échantillonnage temporel d'estimation des dérivées moins dense que l'échantillonnage temporel d'acquisition et en répartissant les instants d'estimation comme proposé dans l'invention, il est possible de contrôler la redondance d'information produite, et ainsi le temps de calcul nécessaire et la quantité d'information à stocker. 20 Dans un cas d'application de l'invention dans lequel l'étape de vérification comprend la détection d'indications de défauts ou d'hétérogénéités dans la structure de matériau considérée, une pluralité d'images représentant chacune une cartographie sur la pièce aéronautique d'une dérivée estimée à un instant déterminé peut être présentée à un 25 opérateur pour détecter si la pièce comprend ou non un ou des défauts. Si comme mentionné précédemment l'application de la deuxième fonction induit un sur-échantillonnage, il suffit de sélectionner un nombre d'instants d'estimation de la dérivée relativement parcimonieux sur le support temporel d'intérêt contenu dans la période d'acquisition, ce qui revient à 30 ne montrer à l'opérateur qu'un nombre relativement faible d'images pour lui permettre d'effectuer une vérification de la pièce. Typiquement, dans l'exemple d'une deuxième fonction appliquée de type logarithme d'une différence (ou d'un écart), une séquence de plus de mille images peut être réduite à une centaine d'images tout en permettant une qualité d'analyse 35 équivalente. Il suffit donc d'estimer la dérivée sur un faible nombre d'instants d'estimation ce qui réduit considérablement les temps de 9 calculs, l'espace de stockage des dérivées calculées et le temps d'analyse de ces dérivées. La complexité de mise en oeuvre de l'invention liée à l'estimation locale des dérivées est ainsi limitée tout en s'assurant de disposer d'une 5 bonne représentation des dérivées estimées (on peut s'affranchir par ce biais de réaliser un nombre trop important d'estimations de dérivées tout en conservant suffisamment d'informations pertinentes pour procéder efficacement à la vérification de la pièce aéronautique). L'estimation locale des dérivées proposée par l'invention permet 10 de s'affranchir de l'existence d'un biais d'estimation dépendant des conditions d'acquisition de la séquence d'images thermiques. Un biais d'estimation peut subsister avec une telle estimation locale, mais celui-ci est invariant aux conditions d'acquisition. En particulier, si une même structure de matériau est observée en différentes zones d'une grande 15 pièce vues par différentes acquisitions unitaires et que les conditions d'acquisition unitaires diffèrent (par exemple en termes de durée d'acquisition), alors les dérivées estimées par le procédé décrit sont identiques à une variance d'estimation près. Le procédé selon l'invention fournit donc des informations (valeurs de dérivées) aisément 20 interprétables sur la pièce dans son intégralité, facilitant l'analyse de la pièce sur la base de ces informations par un opérateur humain ou de façon automatique. On obtient en effet grâce à l'invention des dérivées qui caractérisent localement la structure de matériau de la pièce aéronautique considérée et non plus l'ensemble constitué de la structure 25 de matériau et des conditions d'acquisition. La méthode proposée par l'invention assure une meilleure stabilité des dérivées estimées. En outre, à performance de détection identique, elle permet une réduction du temps d'acquisition des échantillons temporels nécessaires pour l'estimation des dérivées par 30 rapport à l'état de la technique (gain d'environ 30% constaté par les inventeurs). On note que bien que l'invention ait une application privilégiée pour les pièces aéronautiques de grande taille, hétérogènes et/ou de géométrie complexe, elle s'applique néanmoins à toute pièce 35 aéronautique, y compris aux pièces qui ne nécessite qu'une unique acquisition d'images thermiques. L'invention peut également s'appliquer 10 avantageusement à d'autres domaines industriels que le domaine aéronautique. Aucune limitation n'est attachée à l'ordre des dérivées estimées. Ainsi, par exemple ladite au moins une dérivée peut comprendre une 5 dérivée d'ordre 1 et/ou une dérivée d'ordre 2 de la première fonction de la caractéristique. Le procédé selon l'invention permet toutefois d'estimer des dérivées d'ordres supérieurs, comme par exemple d'ordre 3. Par ailleurs, conformément à l'invention, on estime à partir des images numériques fournies par le système de thermographie active une 10 dérivée d'une première fonction de la caractéristique représentée sur les images numériques par rapport à une deuxième fonction du temps. La première et la deuxième fonctions peuvent être diverses. Elles peuvent être choisies en fonction de différents paramètres, comme notamment de la nature de l'excitation utilisée par le système de thermographie active, 15 et/ou le type de matériau composant la pièce aéronautique. Par exemple, dans un mode de réalisation particulier dans lequel : - la caractéristique reflétée par les images numériques est une fonction croissante de la température surfacique de la pièce aéronautique ; et 20 - le système de thermographie infrarouge active utilise un mode d'excitation impulsionnelle par réflexion de la pièce aéronautique ; et - la première fonction et la deuxième fonction sont des fonctions logarithmiques (ex. logarithme naturel) d'un écart entre une valeur et une valeur de référence déterminée. 25 En effet, dans un tel contexte, il est connu des lois de la physique, que la dérivée d'ordre 1 du logarithme naturel de l'écart entre la caractéristique acquise par le système de thermographie active après excitation (caractéristique qui est assimilée à une fonction linéaire croissante de la température) et la valeur de cette caractéristique au repos 30 (i.e. avant excitation) par rapport au logarithme du temps (si on prend comme temps de référence l'instant auquel a lieu l'excitation), est égale à la valeur -0.5 (pour un matériau composant la pièce aéronautique homogène d'épaisseur semi-infinie, et en négligeant le phénomène de convexion). Partant de ce constat, l'utilisation dans le procédé de 35 l'invention de première et deuxième fonctions telles que susmentionnées permet de faciliter la détection d'un défaut dans la pièce aéronautique. Il 11 convient de noter que pour d'autres types d'excitation appliquée à la pièce aéronautique, d'autres fonctions qu'un logarithme d'un écart entre deux valeurs peuvent être envisagées et plus pertinentes. Le procédé de contrôle non destructif selon l'invention bénéficie 5 donc d'une grande flexibilité en matière de configuration et peut être paramétré à différents niveaux, ce qui lui permet d'être appliqué dans de nombreux contextes : choix des première et deuxième fonctions selon les conditions d'acquisition, choix des instants d'estimation, choix de l'intervalle d'estimation de la dérivée, etc. 10 Comme mentionné ci-avant, l'intervalle d'estimation des dérivées est paramétrable conformément à l'invention. Son choix dépend de plusieurs facteurs, comme typiquement des caractéristiques de la pièce aéronautique, des conditions d'acquisition (incluant le mode d'excitation de la pièce par le système de thermographie active), etc. 15 Ainsi, notamment, les inventeurs ont constaté que le bruit affectant les amplitudes converties des pixels des images numériques utilisées pour estimer les dérivées peut être impacté par la nature de la première fonction appliquée pour convertir les amplitudes du signal acquis par le système de thermographie active. Pour réduire l'erreur d'estimation 20 résultant de ce bruit, il peut être judicieux d'envisager des intervalles d'estimation de taille variable en fonction du temps, c'est-à-dire contenant un nombre d'échantillons (i.e. d'amplitudes acquises par le système de thermographie active) qui varie dans le temps. Dans l'exemple précité d'une première fonction et d'une 25 deuxième fonction égales à des fonctions logarithmiques d'écarts, le fait de considérer une fenêtre temporelle (i.e. un voisinage temporel) autour de chaque instant d'estimation constant en fonction du temps (et donc d'utiliser un nombre d'échantillons acquis constant pour le calcul de la dérivée) entraîne typiquement une augmentation dans le temps de la 30 variance d'estimation de la dérivée. Pour pallier cet inconvénient, les inventeurs proposent, dans un mode particulier de réalisation de l'invention, de sélectionner des fenêtres temporelles comprenant un nombre d'échantillons variable afin de contrôler la variance d'estimation en fonction du temps 35 Ainsi, dans un mode particulier de réalisation, l'intervalle d'estimation déterminé autour de chacun des instants d'estimation 12 convertis peut avoir la même taille pour tous les instants d'estimation convertis. Si la deuxième fonction est une fonction logarithmique, et donc croissante dans le temps, cela revient bien à avoir une fenêtre temporelle d'estimation (correspondant dans le domaine temporel à l'intervalle 5 d'estimation) qui augmente dans le temps (mais qui est constante par rapport au logarithme du temps, autrement dit dans le référentiel du temps converti). On obtient de cette sorte une erreur d'estimation indépendante du logarithme du temps. Dans un mode particulier de réalisation, le procédé de contrôle 10 non destructif selon l'invention comprend, préalablement à l'étape d'estimation de la dérivée, une étape de réduction du bruit affectant les amplitudes converties. Une telle étape permet de réduire le bruit affectant les amplitudes converties et ainsi d'améliorer la qualité des dérivées estimées 15 localement. Les inventeurs ont constaté qu'une telle étape de débruitage peut typiquement conduire à une division par 5 voire 10 de l'erreur d'estimation affectant le calcul de la dérivée. Dans un mode particulier de réalisation, certaines étapes du procédé de contrôle non destructif selon l'invention sont déterminées par 20 des instructions de programmes d'ordinateurs. C'est le cas par exemple des étapes mises en oeuvre lors de l'étape de détermination du procédé de contrôle non destructif. En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations ou d'enregistrement, ce 25 programme étant susceptible d'être mis en oeuvre dans un dispositif d'estimation ou plus généralement dans un ordinateur, ce programme comportant des instructions adaptées à la mise en oeuvre de l'étape de détermination du procédé de contrôle non destructif tel que décrit ci-dessus. 30 Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable. 13 L'invention vise aussi un support d'informations ou d'enregistrement lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou 5 dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore un moyen d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. 10 D'autre part, le support d'informations peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. 15 Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en question. 20 Brève description des dessins D'autres caractéristiques et avantages de la présente invention ressortiront de la description faite ci-dessous, en référence aux dessins annexés qui en illustrent un exemple de réalisation dépourvu de tout caractère limitatif. Sur les figures : 25 - la figure 1, déjà décrite, représente une pièce aéronautique telle qu'une structure IFS sur laquelle sont définies une pluralité de zones d'acquisition unitaire pour mener une pluralité d'acquisitions unitaires par un système de thermographie active ; la figure 2 représente un système de contrôle non destructif conforme 30 à l'invention, dans un mode particulier de réalisation ; - la figure 3 représente l'architecture matérielle d'un dispositif d'estimation de dérivées du système de contrôle non destructif de la figure 2 dans un mode particulier de réalisation ; - la figure 4 illustre les principales étapes d'un procédé de contrôle non 35 destructif selon l'invention ; 14 - la figure 5 représente plusieurs dérivées estimées conformément à l'invention et pouvant être utilisées pour réaliser un contrôle non destructif selon l'invention ; et - la figure 6 représente les principales étapes mises en oeuvre par le 5 dispositif d'estimation de la figure 2 conformément à l'invention. Description détaillée de l'invention La figure 2 représente, dans son environnement, un système 1 de contrôle non destructif conforme à l'invention, dans un mode particulier 10 de réalisation. Dans l'exemple envisagé ici, le système 1 permet de procéder au contrôle non destructif de grandes pièces aéronautiques, hétérogènes et de géométrie complexe, telles que par exemple une structure interne fixe (IFS) 2 d'un inverseur de poussée équipant un aéronef. 15 Toutefois, aucune limitation n'est attachée à la nature de la pièce sur laquelle le contrôle non destructif est appliqué. Il peut s'agir plus généralement de tout type de pièces, préférentiellement aéronautiques, comme par exemple une pièce équipant un moteur d'aéronef, etc., ces pièces pouvant être de taille et de géométrie quelconques, hétérogènes 20 ou non, etc. L'invention s'applique en outre également à d'autres domaines industriels qu'au domaine aéronautique. Conformément à l'invention, le système 1 de contrôle non destructif comprend : un système de thermographie active 3 ; 25 un dispositif d'estimation 4 conforme à l'invention ; et un dispositif de vérification 5 configuré pour sanctionner (i.e. valider ou rejeter) la pièce aéronautique 2 à partir des données estimées par le dispositif d'estimation 4. Dans l'exemple envisagé à la figure 2, le système de 30 thermographie active 3 comprend une ou plusieurs sources d'excitation 3A destinées à exciter en transmission et/ou en réflexion la pièce aéronautique 2 alors positionnée sur un support physique prévu à cet effet (non représenté sur la figure), et un dispositif d'acquisition 3B apte à mesurer la réponse (thermique ici) de tout ou partie de la pièce 35 aéronautique 2 aux excitations des sources 4 sur une période déterminée Tacq. 15 Dans le mode de réalisation décrit ici, les sources d'excitation 3A sont ici des sources de chaleur aptes à provoquer une excitation impulsionnelle, de type flash par exemple, de la pièce aéronautique 2. On s'intéresse ici à une excitation par réflexion de la pièce aéronautique 2. 5 En réponse à cette excitation impulsionnelle, la température à la surface de la pièce aéronautique 2 augmente instantanément (par exemple de 10-15 degrés). Puis, la diffusion du flux de chaleur dans la pièce aéronautique 2 se traduit par une décroissance en température à la surface de la pièce aéronautique 2 jusqu'à atteindre ou quasiment la 10 température initiale de la pièce aéronautique 2 avant excitation (autrement dit, la température ambiante). Le dispositif d'acquisition 3B est une caméra thermique infrarouge apte à capter (au moyen d'un capteur approprié), sur une période (durée) de temps dite d'acquisition déterminée, la réponse 15 thermique en surface de la pièce aéronautique 2 à l'excitation impulsionnelle qui lui est appliquée par la ou les sources d'excitation 3A. La caméra thermique infrarouge 3B est en outre apte à fournir une séquence temporelle d'images numériques thermiques reflétant cette réponse sur la période d'acquisition. Chaque image numérique IMj, 20 j=1,...,3 fournie par la caméra infrarouge correspond à un instant d'acquisition tj défini sur la période d'acquisition, J désignant un entier supérieur à 1. Les instants tj sont par exemple espacés de manière uniforme sur la période d'acquisition. On considère ici par souci de simplification que l'origine du temps est l'instant d'excitation de la pièce 25 aéronautique 2 par le système de thermographie infrarouge 3. Chaque image numérique IMj reflète la réponse thermique en surface de la pièce aéronautique 2 à l'excitation par la ou les sources 3A, ou lorsqu'il s'agit d'une pièce aéronautique de grande taille, d'une partie seulement de cette pièce dite zone d'acquisition unitaire de la pièce. La 30 réponse thermique sur l'ensemble de la pièce aéronautique 2 est obtenue en réalisant une pluralité d'acquisitions unitaires, chaque acquisition unitaire ayant potentiellement une période d'acquisition différente, et chaque acquisition unitaire étant réalisée sur une zone d'acquisition unitaire dédiée de la pièce aéronautique 2. La dimension des zones 35 d'acquisition unitaire est préférentiellement choisie en tenant compte des caractéristiques de la caméra thermique infrarouge 3B et des exigences 16 sur les mesures réalisées par la caméra thermique (ex. résolution spatiale, résolution temporelle, rapport signal à bruit) induites par les exigences sur la précision d'estimation des dérivées. En outre, les zones d'acquisition unitaire sont choisies de sorte que leur juxtaposition couvre l'intégralité de 5 la pièce aéronautique 2 (par exemple comme sur la figure 1 précédemment décrite) ou tout du moins une partie choisie préalablement de cette pièce aéronautique 2 sur laquelle on souhaite baser plus particulièrement le contrôle non destructif. Pour réaliser cette pluralité d'acquisitions unitaires, la caméra 10 infrarouge 3B peut être placée par exemple sur un bras apte à se déplacer autour de la pièce aéronautique 2 et à positionner la caméra infrarouge 3B en vis-à-vis de chaque zone d'acquisition unitaire. Dans la suite de la description, on note Z1, Z2,...,ZN, N désignant un entier supérieur ou égal à 1, les zones d'acquisition unitaire 15 considérées lors des acquisitions unitaires (quand N=1, une seule zone d'acquisition unitaire couvrant l'intégralité ou une partie choisie de la pièce aéronautique 2 est considérée), et IMj(Zn), j=1,...,3(n) les images thermiques fournies par la caméra infrarouge 3B de la zone Zn, 3(n) désignant un entier supérieur à 1. Chaque acquisition unitaire correspond 20 à une période ou durée d'acquisition Tacqn, n=1,...,N. Comme mentionné précédemment, chaque image numérique IMj(Zn) de la pièce aéronautique 2 reflète la réponse thermique de la zone d'acquisition unitaire Zn de la pièce aéronautique 2 à l'instant d'acquisition tj, suite à l'excitation (impulsionnelle ici) appliquée par le biais de la ou 25 des sources d'excitation 3A. Elle comprend une pluralité de pixels correspondant à un échantillonnage spatial de la zone d'acquisition unitaire Zn, autrement dit, chaque pixel est associé à un point de la zone d'acquisition unitaire Zn. Chaque pixel est associé sur l'image IMj(Zn) à une amplitude représentant une caractéristique de la pièce aéronautique 2 30 à l'instant d'acquisition tj, cette caractéristique étant ici une fonction croissante déterminée de la température surfacique de la pièce. Dans le mode de réalisation décrit ici, les images numériques thermiques fournies par la caméra infrarouge 3B sont des images en deux dimensions codées en niveaux de gris : chaque pixel d'une image a donc 35 une amplitude qui est ici un niveau de gris reflétant la température surfacique ou plus précisément une fonction croissante de la température 17 surfacique en le point de la pièce aéronautique 2 représenté par le pixel. Cette fonction croissante résulte, de façon connue, de la combinaison de la loi de Planck en monochromatique et de la réponse supposée linéaire du capteur de la caméra infrarouge 3B : elle traduit la conversion en un 5 niveau de gris ici de la luminance spectrale captée par le capteur de la caméra infrarouge 3B en le point représenté par le pixel, cette luminance spectrale reflétant la température surfacique de la pièce aéronautique en ce point (i.e. plus la température est élevée plus la luminance est élevée). Conformément à l'invention, les images numériques acquises 10 par le dispositif d'acquisition 3 et notamment par la caméra infrarouge 3B sont fournies au dispositif d'estimation 4. Dans le mode de réalisation décrit ici, le dispositif d'estimation 4 est configuré pour estimer à partir des images numériques fournies par la caméra infrarouge 3B, au moins une dérivée d'une première fonction mathématique fi de la température par 15 rapport à une deuxième fonction mathématique f2 du temps. Dans l'exemple envisagé ici d'une excitation impulsionnelle de la pièce aéronautique 2 par le dispositif d'acquisition 3, on considère pour fl et f2 la même fonction mathématique, à savoir une fonction logarithme naturel aussi appelée logarithme népérien d'un écart entre une valeur et 20 une valeur de référence déterminée, soit fl(y)=f2(y)=In(y-yref), où y est un nombre réel et yref la valeur de référence considérée. Pour fi, la valeur de référence correspond ici à l'amplitude de niveau de gris du pixel considéré associée à une température au repos du point de la pièce aéronautique 2 représenté par ce pixel. 25 Pour f2, la valeur de référence correspond à l'instant d'excitation de la pièce aéronautique 2 par le système de thermographie infrarouge. Comme mentionné précédemment, on considère ici que les instants d'acquisition ont pour origine l'instant d'excitation de la pièce aéronautique 2 de sorte qu'on applique ici une valeur de référence nulle (i.e. elle est 30 déjà implicitement contenue dans les instants d'acquisition). Toutefois, ces hypothèses ne sont pas limitatives en soi. D'autres fonctions mathématiques peuvent être envisagées pour f1 et f2. Le choix des fonctions fi et f2 peut être guidé par différents paramètres, comme notamment de la nature de l'excitation utilisée par le système de 35 thermographie active, etc. 18 Dans le mode de réalisation décrit ici, le dispositif d'estimation 4 a l'architecture matérielle d'un ordinateur telle que représentée schématiquement à la figure 3. Il comprend, notamment un processeur 6, une mémoire vive 7, 5 une mémoire morte 8, une mémoire flash non volatile 9 ainsi que des moyens de communication 10 permettant notamment au dispositif d'estimation 4 de communiquer avec le dispositif de vérification 5 pour lui fournir par exemple la ou les dérivées estimées à partir des images numériques de la pièce aéronautique 2 fournies par la caméra infrarouge 10 3B. Ces moyens de communication comprennent par exemple un bus de données numériques ou si le dispositif d'estimation 4 et le dispositif de vérification 5 sont reliés via un réseau de télécommunications (local ou autre, filaire ou sans fil, etc.), une carte réseau ou une interface permettant de communiquer sur ce réseau. 15 Dans une variante de réalisation, le dispositif d'estimation 4 et le dispositif de vérification 5 peuvent être localisés au sein d'un même équipement, par exemple un ordinateur tel que celui représenté à la figure 3. La mémoire morte 8 du dispositif d'estimation 4 constitue un 20 support d'enregistrement conforme à l'invention, lisible par le processeur 6 et sur lequel est enregistré un programme d'ordinateur PROG conforme à l'invention. Le programme d'ordinateur PROG définit des modules fonctionnels et logiciels ici, configurés pour mettre en oeuvre une 25 estimation d'une ou de plusieurs dérivées telles que proposées par l'invention à partir des images numériques de la pièce aéronautique 2 fournies par la caméra infrarouge 3B. Ces modules fonctionnels s'appuient sur et/ou commandent les éléments matériels 6-10 du dispositif d'estimation 4 cités précédemment. Ils comportent notamment ici, comme 30 illustré sur la figure 2, une pluralité de modules activés pour chaque dérivée estimée et pour chaque zone d'acquisition unitaire Zn, n=1,...,N considérée, lesdits modules comprenant ici : - un module de sélection 4A, configuré pour sélectionner une pluralité K(n) d'instants d'estimation tck, k=1,...,K(n) de cette dérivée sur la 35 période d'acquisition Tacqn des images numériques (période de temps 19 déterminée au sens de l'invention) pour l'acquisition unitaire indexée par n, n=1,...,N, K(n) désignant un entier supérieur à 1 ; un premier module d'application 4B, configuré pour appliquer la première fonction f1 aux amplitudes des pixels des images 5 numériques, ce premier module d'application fournissant des amplitudes de pixels dites converties ; un deuxième module d'application 4C, configuré pour appliquer la deuxième fonction f2 aux instants d'acquisition tj, j=1,...,3(n) des images numériques pour la zone d'acquisition unitaire Z(n) et aux 10 instants d'estimation tck, k=1,...,K(n) de la dérivée, ce deuxième module d'application 4C fournissant des instants d'acquisition dits convertis notés Tl j=1,...,3(n) et en des instants d'estimation dits convertis notés Tck, k=1,...,K(n) ; un module de débruitage 4D, configuré pour réduire le bruit présent 15 dans les amplitudes des pixels converties ; un module de détermination 4E, configuré pour déterminer pour chaque instant d'estimation converti tck, un intervalle d'estimation (i.e. voisinage) Wk autour de cet instant d'estimation converti. Si on considère que la période d'acquisition Tacqn pour l'acquisition unitaire 20 réalisée sur la zone d'acquisition unitaire Zn est délimitée par un instant initial t_init et un instant final t_fin, chaque intervalle d'estimation Wk déterminé par le module de détermination 4E autour d'un instant d'estimation converti Tck a une taille inférieure à (est inclus dans) la période d'acquisition convertie définie par la différence entre 25 l'instant final converti résultant de l'application de la deuxième fonction f2 à l'instant final t_fin et l'instant initial converti résultant de l'application de la deuxième fonction f2 à l'instant initial t_init (on considère ici la même convention que précédemment pour l'origine des instants t_init et t_fin) ; et 30 un module d'estimation 4F, configuré pour estimer pour chaque instant d'estimation converti tck, la dérivée à partir des amplitudes converties correspondant aux instants d'acquisition convertis tj appartenant à l'intervalle d'estimation Wk déterminé par le module de détermination 4E pour l'instant d'estimation converti. 35 Les fonctions de ces différents modules sont décrites plus en détail ultérieurement en référence à l'étape d'estimation du procédé de 20 contrôle non destructif selon l'invention mise en oeuvre par le dispositif 4 d'estimation du système 1 de contrôle non destructif de la figure 2. Comme mentionné précédemment, le système 1 de contrôle non destructif comprend également un dispositif de vérification 5 de la 5 pièce aéronautique 2. Ce dispositif de vérification 5 est configuré pour permettre de vérifier l'intégrité de la pièce aéronautique 2 (c'est-à-dire de détecter la présence d'anomalies ou d'indications de telles anomalies dans la pièce aéronautique 2) en analysant la ou des dérivées évaluées par le dispositif d'estimation 4. 10 Dans le mode de réalisation décrit ici, le dispositif de vérification 5 comprend un écran et une interface graphique permettant de visualiser des images représentant la ou des dérivées produites par le dispositif d'estimation 4. Ceci permet à un opérateur visualisant ces données de détecter la présence éventuelle d'anomalie sur la pièce aéronautique 2, de 15 positionner le cas échéant cette anomalie sur la surface de la pièce et éventuellement d'évaluer la profondeur du défaut. Ce dispositif de vérification 5 est configuré pour donner des informations intelligibles à l'opérateur pour l'aider à une prise de décision quant à l'intégrité de la pièce aéronautique 2, comme décrit plus en détail ultérieurement. 20 Dans un autre mode de réalisation, le dispositif de vérification 5 peut comprendre un module de détection automatique de défauts sur la pièce aéronautique 2 configuré pour utiliser la ou les dérivées fournies par le dispositif d'estimation 4 et poser un diagnostic quant à l'intégrité de la pièce aéronautique 2 (ex. détection ou non d'indications dans la structure 25 de matériau laissant suspecter un défaut ou une anomalie), sans nécessiter l'intervention d'un opérateur extérieur à proprement parler. La figure 4 représente les différentes étapes du procédé mises en oeuvre par le système 1 de contrôle non destructif illustré à la figure 2 pour opérer un contrôle non destructif de la pièce aéronautique 2 30 conformément à l'invention. Conformément à l'invention, le contrôle non destructif de la pièce aéronautique 2 s'appuie sur un procédé de thermographie infrarouge active. Ce procédé, connu en soi et brièvement rappelé précédemment, est mis en oeuvre par le système de thermographie active 35 3 (étape E10). Plus particulièrement, le système de thermographie active 3 applique, par le biais de sources d'excitation 3A, une excitation 21 impulsionnelle par réflexion ici sur la pièce aéronautique 2 positionnée sur son support prévu à cet effet. La réponse thermique de la pièce aéronautique 2 à cette excitation impulsionnelle est captée par la caméra infrarouge 3B durant 5 une période d'acquisition déterminée, à une pluralité d'instants d'acquisition. Dans l'exemple envisagé ici, la pièce aéronautique 2 étant une structure IFS de grande taille, plusieurs (i.e. N, N désignant un entier supérieur à 1) acquisitions unitaires sont successivement réalisées sur la pièce aéronautique 2 comme mentionné précédemment. Chaque 10 acquisition unitaire est réalisée sur une période d'acquisition unitaire égale à Tacqn, chaque acquisition unitaire visant une zone d'acquisition unitaire Zn distincte, n=1,...,N, les N zones d'acquisition unitaire Zn juxtaposées permettant, comme mentionné précédemment, de couvrir la pièce aéronautique 2 dans son intégralité. 15 La caméra infrarouge 3B génère, pour chaque acquisition unitaire indexée par n, une pluralité d'images numériques IMj, j=1,...,J(n) de la pièce aéronautique 2 représentant la réponse ainsi captée de la pièce aéronautique aux différents instants d'acquisition tj, j=1,...,3(n) (étape E20). 20 Dans le mode de réalisation décrit ici, on suppose que les instants d'acquisition tj, j=1,...,3(n) définis sur la période de temps Tacqn sont espacés de façon uniforme avec un pas suffisamment faible pour avoir une bonne représentation de la réponse thermique de la zone d'acquisition unitaire Zn de la pièce aéronautique 2 à l'excitation 25 impulsionnelle appliquée à la pièce, et une redondance d'observation permettant une estimation précise en présence de bruit statistique d'observation. Les images numériques IMj(Zn), 1=1,...,3, acquises sur chacune des zones d'acquisition unitaire Zn, n=1,...,N considérées représentent 30 l'évolution de la température à la surface de la pièce aéronautique 2 suite à l'excitation impulsionnelle qui lui a été appliquée au moyen de la ou des sources d'excitation 3A. Plus particulièrement, chaque pixel noté PIXi, i=1,...,I, d'une image numérique IMj(Zn) correspond à un point de la zone d'acquisition unitaire Zn résultant d'un échantillonnage spatial de cette 35 dernière, I désignant un entier supérieur à 1 correspondant au nombre de pixels dans chaque image IMj(Zn). Chaque pixel PIXi de l'image IMj(Zn) a 22 une amplitude notée amp(PIXi) qui représente une fonction croissante de la température surfacique de la pièce aéronautique 2 (caractéristique de la pièce aéronautique au sens de l'invention), à l'instant d'acquisition tj, en le point de la zone d'acquisition unitaire Zn correspondant à ce pixel. 5 Les images numériques IMj(Zn), j=1,...,J, n=1,...,N acquises par la caméra numérique 3B sont ensuite fournies par le système de thermographie active 3 au dispositif d'estimation 4 (étape E30). Conformément à l'invention, le dispositif 4 d'estimation estime, à partir des images numériques thermiques qui lui sont fournies, au moins 10 une dérivée d'ordre x (entier) notée d(x) d'une première fonction f1 de la température par rapport à une deuxième fonction f2 du temps (étape E40). Dans l'exemple envisagé ici, comme mentionné précédemment, les fonctions fi et f2 sont toutes les deux prises égales à une fonction logarithme naturel (i.e. népérien) d'un écart (i.e. d'une différence) entre 15 une valeur et une valeur de référence déterminée (i.e. d'un écart d'amplitudes pour fi et d'un écart temporel pour f2). Ces fonctions sont particulièrement bien adaptées à une excitation de type impulsionnelle, notamment par réflexion, comme expliqué ci-avant. Aucune limitation n'est attachée à l'ordre x de la dérivée 20 estimée par le dispositif 4 d'estimation. Il s'agit par exemple d'une dérivée d'ordre 1 ou de tout ordre supérieur à 1. Par ailleurs, plusieurs dérivées d'ordres différents peuvent être estimées par le dispositif 4 d'estimation. Par souci de simplification, dans l'exemple décrit ici, on se limite à l'estimation d'une unique dérivée d'ordre 1, autrement dit d(x)=d(1). 25 Conformément à l'invention, la dérivée cl(x) estimée par le dispositif 4 d'estimation est estimée localement pour chaque pixel de chaque zone d'acquisition unitaire Zn, en plusieurs instants d'estimation tck, k=1,...,K, sur un voisinage temporel autour de ces instants d'estimation. Ce voisinage temporel est réduit par rapport à la période 30 d'acquisition Tacqn, comme décrit plus en détail ultérieurement. On obtient de cette sorte une estimation (locale) de la dérivée d(x) recherchée robuste aux conditions d'acquisition des images thermiques. Le dispositif d'estimation 4 fournit la dérivée estimée d(x) formée des dérivées estimées localement en chacun des instants d'estimation tck, 35 k=1,...,K(n) et en chaque pixel de chacune des zones d'acquisition unitaire Zn au dispositif de vérification 5 (étape E50). 23 Dans le mode de réalisation décrit ici, le dispositif de vérification 5 est configuré pour afficher sur son écran, via une interface graphique prévue à cet effet, les dérivées ainsi estimées localement. Cet affichage peut être déclenché sous l'action par exemple d'un opérateur qui, à partir 5 des dérivées affichées, est en mesure via sa propre expertise, de réaliser un diagnostic de la pièce aéronautique 2, c'est-à-dire de vérifier l'intégrité de la structure de matériau par rapport à un état attendu et de déterminer si la pièce aéronautique 2 comporte ou non un défaut (étape E60). Un exemple d'un tel affichage est donné en figure 5. Sur cette 10 figure, les différentes courbes représentées en traits discontinus correspondent aux dérivées d(1) d'ordre 1 estimées conformément à l'invention pour différents pixels associés à des points distincts d'une même zone d'acquisition unitaire. En abscisse sont représentés les indices k des différents instants d'estimation, tandis qu'en ordonnée sont 15 représentées les amplitudes des dérivées estimées. Pour rendre compte de la précision des dérivées estimées, les dérivées premières vraies (autrement dit « théoriques » ou « attendues ») sont également représentées sur la figure à titre de référence en traits continus. L'utilisation de telles informations par un opérateur pour réaliser 20 un diagnostic de la pièce aéronautique 2 (aussi appelé vérification ou « sanction » de la pièce) est connue en soi. En examinant l'allure de la dérivée d(1) d'ordre 1 du logarithme naturel de l'écart d'une fonction croissante déterminée de la température par rapport au logarithme naturel du temps écoulé après l'excitation impulsionnelle fournie par le dispositif 25 d'estimation 4, l'opérateur peut déterminer si celle-ci a une allure conforme à la structure de matériau attendue ou si elle présente des différences notables par rapport à celle-ci, laissant suspecter la présence d'hétérogénéités ou de défauts dans la pièce aéronautique 2. A titre illustratif, la figure 5 représente les dérivées d(1) estimées obtenues pour 30 différents pixels présentant tous des défauts, lors d'un contrôle en réflexion avec une excitation par lampe flash. On note que plus le défaut apparaît de manière précoce dans le temps, plus il est proche de la surface de la pièce aéronautique 2 située du côté de la caméra. Plus particulièrement, dans le cas d'application de la détection 35 d'indications par un opérateur, la sanction de la pièce aéronautique 2 peut comprendre par exemple l'observation successive pour chaque 24 instant d'estimation, d'images dites spatiales représentant une cartographie spatiale des dérivées estimées à cet instant d'acquisition. Un défaut tel qu'une lame d'air en profondeur dans la structure du matériau constitue une hétérogénéité et se traduit par un contraste entre la zone 5 présentant le défaut et une zone voisine saine dans au moins une image spatiale de dérivées. Grâce à l'invention, on s'assure de l'absence de sauts de niveaux de gris entre les images juxtaposées de deux zones d'acquisition unitaire juxtaposées (autrement dit, l'invention permet de s'affranchir de l'existence d'un biais d'estimation dépendant des conditions 10 d'acquisition, comme par exemple de la période d'acquisition unitaire qui peut varier d'une acquisition unitaire à l'autre). On s'assure également que le contraste visuel entre deux zones voisines, l'une saine et l'autre à défaut, toutes deux de structure de matériau fixée, est indépendant des conditions d'observation de la zone d'acquisition unitaire contenant cette 15 zone saine et cette zone à défaut. En variante, le vecteur de dérivées de chaque pixel de la pièce aéronautique 2 peut être comparé à des vecteurs de référence, préalablement appris et couvrant la variabilité des « possibles » attendus pour la structure de matériau. Cela permet d'avoir une caractéristique de 20 la structure de matériau en chaque point de la pièce et d'évaluer la conformité de la structure de matériau inspectée. L'invention permet de réduire la dimension de la base des vecteurs de référence considérée par rapport à l'état de la technique, en garantissant une estimation des dérivées indépendante des conditions d'acquisition. En effet, pour établir 25 la base des vecteurs de référence, il est nécessaire de tenir compte uniquement de la variabilité des structures de matériau (et non de la variabilité des structures de matériau couplée à la variabilité des conditions d'acquisition). De plus, les conditions d'acquisition peuvent être avantageusement modifiées après constitution de la base, ce qui offre une 30 souplesse d'utilisation. Le document WO 2001/041421 donne d'autres exemples d'exploitation de la ou des dérivées d(x) estimées par le dispositif d'estimation 4 pour détecter la présence d'un défaut dans une pièce aéronautique. 25 Ces exemples ne sont bien entendu donnés qu'à titre illustratif, et dépendent notamment des dérivées cl(x) estimées par le dispositif d'estimation 4 et de la structure de matériau inspectée. Dans un autre mode de réalisation, le dispositif de vérification 5 5 comprend un module de détection automatique de défauts sur la pièce aéronautique 2 configuré pour utiliser la ou les dérivées d(x) fournies par le dispositif d'estimation 4 et poser un diagnostic quant à l'intégrité de la pièce aéronautique 2, sans nécessiter l'intervention d'un opérateur externe à proprement parler. Le document EP 1 910 786 décrit par exemple une 10 méthode automatique permettant d'exploiter des informations contenues dans des dérivées cl(x) pouvant être calculées par le dispositif d'estimation 4. Cet exemple n'est toutefois donné qu'à titre illustratif. Comme mentionné précédemment, l'invention se distingue principalement de l'état de la technique existant par la façon dont elle 15 propose d'estimer, au cours de l'étape E40 du procédé de contrôle non destructif précédemment décrit, la ou les dérivées cl(x) utilisées pour réaliser le contrôle non destructif de la pièce aéronautique 2. L'invention propose en effet une estimation de cette ou de ces dérivées à la fois moins biaisée et de biais insensible aux variations usuelles des conditions 20 d'acquisitions, telles que la période d'acquisition. Cette estimation est réalisée sur un échantillonnage temporel dédié, possiblement et avantageusement différent de l'échantillonnage d'acquisition considéré par le dispositif d'acquisition 3 pour obtenir les images thermiques (i.e. les instants d'estimation tck, k=1,...,K peuvent être différents des instants 25 d'acquisition tj, j=1,...,J). En outre, l'invention s'appuie sur une estimation locale de la ou des dérivées d(x) à partir d'échantillons mesurés dans un voisinage temporel de taille limitée par rapport à la période d'acquisition Tacq des images thermiques, ce qui garantit, par construction, l'invariance de l'estimation obtenue par rapport à la période de l'acquisition et à la 30 zone d'acquisition unitaire considérée. Nous allons maintenant décrire plus en détail la façon dont la ou les dérivées d(x) sont estimées par le dispositif d'estimation 4 du système de contrôle non destructif 1 lors de l'étape E40, conformément à l'invention. On se limite ici par souci de simplification au calcul d'une dérivée unique d'ordre 1 (i.e. d(x)= dm). 35 La figure 6 illustre les principales étapes mises en oeuvre lors de l'étape E40 par le dispositif d'estimation 4 pour estimer la dérivée d(x) à 26 partir des images IMj(Zn), j=1,...,), n=1,...,N fournies par la caméra infrarouge 3B pour chaque zone d'acquisition unitaire Zn. Dans le mode de réalisation décrit ici, le dispositif d'estimation 4, via son module de sélection 4A, sélectionne tout d'abord la pluralité 5 d'instants d'estimation tck, k=1,...,K(n) auxquels il va estimer la dérivée cl(x) en chaque pixel de chaque zone d'acquisition unitaire sur la période d'acquisition Tacqn (étape F10). A cet effet, le module de sélection 4A utilise ici une règle de sélection prédéterminée, fixée préalablement par un opérateur. 10 Conformément à l'invention, cette règle fixe des instants d'estimation convertis de la dérivée uniformément espacés dans le graphe logarithmique résultant de la conversion du domaine temporel correspondant à la période d'acquisition Tacqn par la fonction f2. Autrement dit, en reprenant les notations introduites précédemment, si 15 on considère que la période d'acquisition unitaire Tacqn est délimitée par un instant initial t_init et un instant final tfin, les instants d'estimation convertis sont uniformément espacés (même écart en les instants d'estimation convertis) sur la période d'acquisition convertie définie par la différence entre l'instant final converti résultant de l'application de la 20 deuxième fonction f2 à l'instant final t_fin et l'instant initial converti résultant de l'application de la deuxième fonction f2 à l'instant initial t_init. En outre, selon la règle de construction appliquée, le nombre d'instants d'estimation sélectionnés K(n) pour la zone d'acquisition unitaire Zn est inférieur au nombre d'instants d'acquisition tj, j=1,...,3(n). 25 Une telle règle de sélection appliquée par le module de sélection 4A est par exemple la suivante : « à partir d'un instant TO défini sur la période d'acquisition Tacq et pour un paramètre entier a positif, tant que tck Tacq : tc(k+1)= or.tck avec tc1=TO ». L'instant TO et le paramètre a sont fournis par l'opérateur. TO peut par exemple être choisi de sorte à 30 coïncider avec le premier instant d'acquisition tl de la période d'acquisition Tacqn. Bien entendu, cet exemple n'est donné qu'à titre illustratif et d'autres règles de sélection ou de façons de déterminer les instants d'estimation tck, k=1,...,K sur la période d'acquisition Tacq peuvent être 35 envisagées. Par exemple, on peut envisager de définir des instants d'estimation coïncidant avec certains instants d'acquisition tj (ex. tc1=t1, 27 tc2=t3, tc3=t5, etc.). Préférentiellement, on choisit un nombre K(n) d'instants d'estimation inférieur au nombre J(n) d'instants d'acquisition sur la période de temps Tacqn pour limiter la complexité de mise en oeuvre de l'invention et notamment le temps de calcul nécessaire. 5 Puis, le dispositif d'estimation 4 procède à un formatage des données dont il dispose et qui sont contenues dans les images numériques qui lui ont été fournies par la caméra infrarouge 3B (étape F20). Plus précisément, lors de ce formatage, il détermine tout d'abord les instants d'acquisition tj, j=1,...,3(n) des images numériques qui 10 lui sont fournies pour chaque zone d'acquisition unitaire Zn, n=1,...,N de la pièce aéronautique 2. On suppose ici qu'à cet effet, les paramètres d'acquisition de ces images lui ont été fournis, à savoir notamment l'instant t_init(Zn) auquel l'acquisition de la réponse thermique de la pièce aéronautique 2 à l'excitation impulsionnelle à laquelle elle a été soumise a 15 débuté (fixée ici comme mentionné précédemment à l'instant d'excitation de la pièce aéronautique 2), pour chaque zone d'acquisition unitaire Zn considérée, ainsi que le pas d'échantillonnage d'acquisition retenu pour chaque acquisition unitaire (i.e. durée entre chaque trame acquise). Puis le dispositif d'estimation 4, via son premier module 20 d'application 4B, convertit, pour chaque zone d'acquisition unitaire Zn, n=1,...,N, les amplitudes des pixels PIXi des images numériques correspondant à cette zone d'acquisition unitaire, en leur appliquant la première fonction f1. Il génère ainsi ici des amplitudes de pixels, notées y(PIXi) dites converties, pour chaque pixel PIXi, i=1,...,I de chaque image 25 IMj(Zn), j=1,...,3(n) de la zone d'acquisition unitaire Zn, telles que : y(PIXi)=f1(amp(PIX0)=In(amp(PIX1)-amp0) où amp0 désigne le niveau de gris associé au pixel PIXi à une température de la pièce aéronautique 2 au repos (c'est-à-dire avant excitation) (valeur de référence au sens de l'invention pour la fonction fi). La valeur de 30 référence amp0 peut être déterminée aisément et fournie au dispositif d'estimation 4. Ceci permet de tenir compte uniquement de la perturbation en amplitude induite par l'excitation. Suite à ce formatage, le dispositif d'estimation 4 procède à un changement de variable temporelle (étape F30), afin de s'adapter 35 notamment au type d'excitation appliquée à la pièce aéronautique par le système de thermographie active 3. 28 Plus précisément, au cours de cette étape, il applique, via son deuxième module d'application 4C, la deuxième fonction f2 aux instants d'acquisition tj, j=1,...,J(n) précédemment déterminés pour chaque zone d'acquisition unitaire Zn et chaque image numérique IMj(Zn), ainsi qu'aux 5 instants d'estimation tck, k=1,...,K(n) de la dérivée cl(x). Il obtient ainsi des instants d'acquisition dits convertis notés zj, j=1,...,J(n) et des instants d'estimation dits convertis notés tck, k=1,...,K(n) tels que : tj=f2(tj)=In(tj), j=1,...,J(n) wk=f2(tck)=In(tck), k=1,...,K(n) 10 On note que le logarithme naturel In est ici appliqué directement aux instants d'acquisition et d'estimation (i.e. valeur de référence nulle pour la fonction f2), ceux-ci ayant pour même origine par hypothèse l'excitation de la pièce aéronautique 2 (autrement dit, la valeur de référence est en quelque sorte intégrée déjà dans la valeur des instants 15 d'acquisition et d'estimation considérée). Dans le mode de réalisation décrit ici, le dispositif d'estimation 4, via son module de débruitage 4D, applique ensuite un traitement de débruitage (i.e. de réduction du bruit) à l'ensemble des amplitudes de pixels converties obtenues à l'issue de l'étape F20 (étape F40). On note 20 que l'étape F40 peut être réalisée indifféremment avant, après ou même en parallèle de l'étape de changement de variable temporelle F30. Cette étape est optionnelle. L'étape de débruitage F40 (ou aussi désignée par étape de pré-restauration) vise à réduire, préalablement à l'estimation de la dérivée d(x), 25 le bruit affectant les amplitudes de pixels converties en tenant compte des propriétés de redondance du signal utile et des statistiques du bruit. Un tel débruitage peut être réalisé à partir de nombreuses techniques connues en soi appliquées successivement sur les amplitudes converties y(PIXi) correspondant à chaque image IMj, j=1,...,J(n) d'une même zone 30 d'acquisition unitaire Zn. De telles techniques consistent par exemple en un filtrage spatial appliqué sur chaque image, en un débruitage par patchs (aussi connu sous le nom anglais de « non-local means »), en une technique dite de morphologie mathématique, et ou une approche variationnelle régularisée (Tikhonov, variation totale, etc.), etc. 35 Aucune limitation n'est attachée à la technique de débruitage utilisée. Dans le mode de réalisation décrit ici, on suppose que le module 29 de débruitage 4D met en oeuvre un filtrage spatial gaussien 2D isotrope sur les amplitudes converties des pixels PIXi, image par image. Les amplitudes converties débruitées sont désignées dans la suite de la description par y'(PIXi). On note que les paramètres du filtrage réalisé (et 5 notamment le nombre de pixels sur lequel il est réalisé) dépendent notamment de la dimension des anomalies que l'on recherche sur la pièce. Le dispositif d'estimation 4 estime alors la dérivée d(x) à partir des amplitudes converties débruitées y'(PIXi) obtenues pour chaque image IMj et pour chaque zone d'acquisition unitaire Zn. Comme mentionné 10 précédemment, conformément à l'invention, la dérivée d(x) est estimée localement, sur des voisinages sélectionnés autour de chaque instant d'estimation tck, ou de façon équivalente autour de chaque instant d'estimation converti Tck, k=1,...,K(n). Le dispositif d'estimation 4, via son module de détermination 4E, 15 détermine donc en premier lieu, pour chaque instant d'estimation converti tck (et pour chaque zone d'acquisition unitaire Zn), un intervalle d'estimation Wk autour de cet instant d'estimation converti, cet intervalle d'estimation ayant une durée inférieure à la période d'acquisition convertie précédemment définie (étape F50). Les mêmes intervalles d'estimation Wk 20 peuvent être sélectionnés pour les différentes zones d'acquisition unitaire Zn, n=1,...,N. Pour déterminer les intervalles d'estimation Wk, le module de détermination 4E utilise ici une règle de détermination prédéfinie, fixée préalablement en fonction des caractéristiques de la pièce aéronautique 2 25 et/ou des conditions d'acquisition des images numériques thermiques (ex. de la méthode d'excitation appliquée par le système de thermographie infrarouge active). Dans l'exemple envisagé ici, cette règle de détermination consiste à sélectionner un intervalle d'estimation Wk centré autour de l'instant d'estimation converti tck et de taille A supérieure à 0, 30 identique pour tous les instants d'estimation convertis Tck, k=1,...,K. Ceci revient à considérer, dans l'exemple de la fonction f2, un voisinage temporel autour des instants d'estimation tck (avant changement de variable temporelle, i.e. avant conversion) qui croit dans le temps, de manière proportionnelle. Un tel choix conduit, pour une acquisition à un 35 échantillonnage régulier dans le temps (i.e. tj régulièrement espacés sur la période Tacqn) et un bruit d'observation blanc sur les amplitudes 30 mesurées, à une variance d'estimation de la dérivée cl(x) pratiquement indépendante des instants d'acquisition convertis Tj, j=1,...,3(n). En variante, d'autres règles de détermination des intervalles d'estimation Wk peuvent être envisagées. Typiquement, les intervalles 5 d'estimation Wk ne sont pas nécessairement centrés autour des instants d'estimation convertis. Ils ne sont en outre pas nécessairement identiques pour chacun des instants d'estimation convertis. Suite à l'étape de détermination F50, le dispositif d'estimation 4 via son module d'estimation 4F, estime pour chaque pixel PIXi, i=1,...,I 10 correspondant à l'échantillonnage spatial de chaque zone d'acquisition unitaire Zn, et pour chaque instant d'estimation converti Tck, k=1,...,K(n), la dérivée notée d(x)(n,i,k) en cet instant à partir des amplitudes converties débruitées correspondant aux instants d'acquisition convertis Tj appartenant à l'intervalle d'estimation Wk (autrement dit à partir des 15 amplitudes converties débruitées extraites des images IMj de la zone d'acquisition unitaire Zn associées aux instants d'acquisition tj correspondant aux instants d'acquisition convertis tj) (étape F60). Il procède de la même façon pour chaque zone d'acquisition unitaire Zn considérée. 20 Différentes techniques peuvent être utilisées par le module d'estimation 4F pour estimer localement la dérivée d(x)(n,i,k) sur l'intervalle d'estimation Wk. Ainsi, par exemple, le module d'estimation 4F peut estimer la dérivée d(x)(n,i,k) par une technique de dérivation d'un filtre gaussien temporel ayant un écart-type constant. 25 De façon connue, une telle technique de dérivation s'appuie ur le principe suivant. Si on suppose une fonction yi(r) (yi est ici homogène à l'amplitude convertie du pixel indexé par i) continue par rapp9rt à une , 1 f variable réelle T et une fonction gaussienne hc(r) = *,-2-1 e-;) d'écart-type a, alors la fonction y(r) = y * h0.(r) (où « * » désigne un produit de 30 convolution) est une version lissée de y(T) à une échelle définie par Q. La convolution par la fonction gaussienne revient en effet à réaliser une moyenne locale pondérée par la fonction gaussienne donnant plus de poids aux valeurs de variable proches de T. On peut par exemple fixer = 0.125 x tailleW où tailleW est la taille d'un voisinage à prendre en 35 compte. 31 La dérivée nième de la fonction continue lissée évaluée en r est alors donnée par : d5j 7()n) (r) = e [y * mer) = Y * eh,(,). Dans le mode de réalisation décrit ici, le module d'estimation 4F 5 applique ce principe (convolution par un filtre gaussien et dérivation) en utilisant des filtres gaussiens pour chaque dérivée estimée dont l'implémentation doit tenir compte d'une part, de la nature discrète des amplitudes converties des pixels, et d'autre part de la nature non-uniforme de l'échantillonnage des instants d'estimation. 10 En variante, le module d'estimation 4F peut réaliser une régression paramétrique locale sur une famille de fonctions dérivables (ex. des polynômes), puis dériver analytiquement la fonction paramétrée résultant de la régression. Bien entendu ces exemples ne sont donnés qu'à titre illustratif 15 et d'autres techniques d'estimation locale des dérivées peuvent être envisagées en variante. Comme mentionné précédemment les dérivées ainsi estimées par une approche locale sont fournies ensuite au dispositif de vérification 5 pour permettre leur interprétation, soit par un opérateur externe, soit 20 par le dispositif de vérification 5 lui-même lorsque cette interprétation est automatisée. Dans le mode de réalisation décrit ici, on a considéré une excitation impulsionnelle par réflexion, des dérivées d'ordre 1 ou 2, et des fonctions f1 et f2 identiques, égales au logarithme naturel d'un écart entre 25 une valeur et une valeur de référence. Comme mentionné, ces hypothèses ne sont pas limitatives en soi, et l'invention s'applique à tout type d'excitation, tout ordre de dérivée, et tout type de fonctions f1 et f2. La prise en compte d'autres ordres de dérivées et/ou d'autres fonctions f1 et f2 permet de s'adapter à d'autres profils d'excitation, tels que des 30 excitations créneaux par flux d'air REVENDICATIONS1. Procédé de contrôle non destructif d'une pièce aéronautique (2) comprenant : - une étape d'obtention (E20) d'une pluralité d'images numériques de la pièce aéronautique acquises à une pluralité d'instants d'acquisition définis sur une période de temps déterminée dite période d'acquisition au moyen d'un système (3) de thermographie infrarouge active, chaque pixel d'une image numérique acquise à un instant d'acquisition ayant une amplitude représentant une caractéristique de la pièce aéronautique à cet instant d'acquisition en un point de la pièce aéronautique ; une étape d'estimation (E40), à partir des images numériques acquises, d'au moins une dérivée d'une première fonction de ladite caractéristique par rapport à une deuxième fonction du temps ; et une étape de vérification (E60) de la pièce aéronautique utilisant ladite au moins une dérivée estimée ; ledit procédé étant caractérisé en ce que l'étape d'estimation (E40) comprend, pour chaque dérivée estimée : - une étape de sélection (F10) d'une pluralité d'instants d'estimation de cette dérivée sur ladite période d'acquisition ; - une étape d'application (F20) de la première fonction aux amplitudes des pixels des images numériques, ladite étape d'application résultant en des amplitudes de pixels converties ; - une étape d'application (F30) de la deuxième fonction aux instants d'acquisition des images numériques résultant en des instants d'acquisition convertis, et aux instants d'estimation de la dérivée résultant en des instants d'estimation convertis, lesdits instants d'estimation étant sélectionnés lors de l'étape de sélection de sorte que les instants d'estimation convertis soient uniformément espacés ; et - une étape d'estimation (F60), pour chaque instant d'estimation converti, de ladite dérivée à partir des amplitudes converties correspondant aux instants d'acquisition convertis appartenant à un intervalle d'estimation déterminé (F50) autour de l'instant d'estimation converti. 33 2. Procédé de contrôle non destructif selon la revendication 1 dans lequel la période d'acquisition est délimitée par un instant dit initial et un instant dit final, et chaque intervalle d'estimation déterminé autour d'un instant d'estimation converti a une taille inférieure à la différence 5 entre l'instant final converti résultant de l'application de la deuxième fonction à l'instant final et l'instant initial converti résultant de l'application de la deuxième fonction à l'instant initial. 3. Procédé de contrôle non destructif selon la revendication 1 ou 10 2 dans lequel : - ladite caractéristique est une fonction croissante de la température surfacique de la pièce aéronautique ; - le système de thermographie infrarouge (3) active utilise un mode d'excitation impulsionnelle par réflexion de la pièce aéronautique ; et 15 - la première fonction et la deuxième fonction sont des fonctions logarithmiques d'un écart entre une valeur et une valeur de référence déterminée. 4. Procédé de contrôle non destructif selon l'une quelconque des 20 revendications 1 à 3 dans lequel ladite au moins une dérivée comprend une dérivée d'ordre 1 et/ou une dérivée d'ordre 2 de la première fonction de ladite caractéristique par rapport à la deuxième fonction du temps. 5. Procédé de contrôle non destructif selon l'une quelconque des 25 revendications 1 à 4 dans lequel l'intervalle d'estimation déterminé autour de chacun des instants d'estimation convertis a la même taille pour tous les instants d'estimation convertis. 6. Procédé de contrôle non destructif selon l'une quelconque des 30 revendications 1 à 5 dans lequel le nombre d'instants d'estimation de la dérivée est choisi lors de l'étape de sélection inférieur au nombre d'instants d'acquisition sur la période d'acquisition. 7. Procédé de contrôle non destructif selon l'une quelconque des 35 revendications 1 à 6 comprenant en outre, préalablement à l'étape 34 d'estimation de la dérivée, une étape (F40) de réduction du bruit affectant les amplitudes converties. 8. Système (1) de contrôle non destructif d'une pièce 5 aéronautique comprenant : - un système de thermographie active (3) comprenant un dispositif d'excitation (3A) et un dispositif d'acquisition (3B) configurés pour acquérir à une pluralité d'instants d'acquisition définis sur une période de temps déterminée dite période d'acquisition, des images 10 numériques de la pièce aéronautique, chaque pixel d'une image numérique acquise à un instant d'acquisition ayant une amplitude représentant une caractéristique de la pièce aéronautique à cet instant d'acquisition en un point de la pièce aéronautique après excitation de la pièce aéronautique par le dispositif d'excitation ; 15 - un dispositif d'estimation (4), configuré pour estimer à partir des images numériques acquises, au moins une dérivée d'une première fonction de ladite caractéristique par rapport à une deuxième fonction du temps ; et - un dispositif de vérification (5) configuré pour vérifier la pièce 20 aéronautique en utilisant ladite au moins une dérivée estimée ; ledit système étant caractérisé en ce que le dispositif d'estimation comprend une pluralité de modules activés pour chaque dérivée estimée et comprenant : - un module de sélection (4A), configuré pour sélectionner une pluralité 25 d'instants d'estimation de cette dérivée sur ladite période d'acquisition ; - un premier module d'application (4B), configuré pour appliquer la première fonction aux amplitudes des pixels des images numériques, ledit premier module d'application fournissant des amplitudes de pixels 30 converties ; - un deuxième module d'application (4C), configuré pour appliquer la deuxième fonction aux instants d'acquisition des images numériques fournissant des instants d'acquisition convertis et pour appliquer la deuxième fonction aux instants d'estimation de la dérivée fournissant 35 des instants d'estimation convertis, lesdits instants d'estimation ayant 35 été sélectionnés par le module de sélection (4A) de sorte que les instants d'estimation convertis soient uniformément espacés ; - un module d'estimation (4F), configuré pour estimer pour chaque instant d'estimation converti, ladite dérivée à partir des amplitudes 5 converties correspondant aux instants d'acquisition convertis appartenant à un intervalle d'estimation déterminé autour de l'instant d'estimation converti.