FR 3055728 A1 20180309 FR 1658197 A 20160902 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 de pièces aéronautiques, telles que des pièces équipant des moteurs d'aéronef comme par exemple des aubes de turbine haute-pression ou basse-pression, des distributeurs, etc. 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 par exemple 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. Parmi les méthodes de contrôle non destructif existantes, certaines s'appuient sur des images numériques fournies par des systèmes de radiographie ou de tomographie numérique. L'avantage de ces images est qu'elles fournissent des informations directement exploitables sur l'intérieur des structures ou des matériaux, et permettent ainsi de détecter des défauts internes pouvant affecter ces structures ou ces matériaux, tels que par exemple la présence d'inclusions ou de retassures. Dans la suite de la description, on désigne plus généralement par « pièce » la structure ou le matériau sur lequel on envisage un contrôle non destructif. De nombreuses variables interviennent dans le réglage des systèmes d'imagerie numérique précités pour acquérir des images numériques des pièces sur lesquelles on envisage d'effectuer un contrôle non destructif. De telles variables sont notamment le grandissement, le temps d'exposition de la pièce aux rayons électromagnétiques (ex. rayons X), la taille des pixels, l'efficacité quantique de détection du capteur équipant le système d'imagerie numérique, etc. Aujourd'hui, ces réglages sont définis et validés, pour différentes gammes de mesure, par des 3055728 2 comités de certification appelés COSAC (COmité Sectoriel Aérospatial de Certification) 2 et 3. Ils s'appuient très largement sur l'expérience des professionnels qui manipulent les systèmes d'imagerie numérique et sur la mesure de facteurs de qualité des images (ou IQI pour Indicateurs de 5 Qualité d'Images) sur une grande gamme de conditions. L'optimisation des paramètres de réglage des systèmes d'imagerie numérique tels que les systèmes de radiographie ou de tomographie implique donc actuellement un travail et un temps considérables, et est en grande partie liée à la dextérité de la personne qui l'effectue. 10 Objet et résumé de l'invention La présente invention permet notamment de remédier à ces inconvénients en proposant un procédé de contrôle non destructif réalisé sur une pièce aéronautique comprenant : 15 - une étape de réglage d'un système d'imagerie numérique, ledit système d'imagerie numérique étant un système de radiographie ou de tomographie numérique comprenant une source de rayons électromagnétiques et un détecteur apte à détecter les rayons électromagnétiques émis par la source ; 20 - une étape d'acquisition d'au moins une image numérique de la pièce aéronautique au moyen du système de radiographie ou de tomographie numérique réglé ; et - une étape de contrôle non destructif réalisé sur la pièce aéronautique en utilisant ladite au moins une image numérique acquise ; 25 Ce procédé de contrôle non destructif est remarquable en ce qu'il comprend en outre une étape préalable d'estimation de paramètres de réglage du système d'imagerie numérique, cette étape d'estimation comprenant : - une étape d'obtention d'un modèle numérique caractérisant le 30 système d'imagerie numérique, ledit modèle numérique comprenant un modèle caractérisant la source et un modèle caractérisant le détecteur ; - une étape d'obtention d'un modèle numérique caractérisant la pièce aéronautique ; - une étape d'évaluation, pour une pluralité de valeurs distinctes des 35 paramètres de réglage, d'un rapport contraste sur bruit obtenu pour le système d'imagerie numérique et ladite pièce aéronautique, cette étape 3055728 3 d'évaluation utilisant les modèles numériques obtenus caractérisant le système d'imagerie numérique et la pièce aéronautique ; et - une étape de détermination de valeurs dites optimales des paramètres de réglage optimisant le rapport contraste sur bruit, ces 5 valeurs optimales étant appliquées au système d'imagerie numérique lors de l'étape de réglage. Corrélativement, l'invention vise également un système de contrôle non destructif d'une pièce aéronautique comprenant : - un système d'imagerie numérique, ce système d'imagerie numérique 10 étant un système de radiographie ou de tomographie numérique comprenant une source de rayons électromagnétiques et un détecteur apte à détecter les rayons électromagnétiques émis par la source ; et - un dispositif de contrôle non destructif de la pièce aéronautique configuré pour utiliser au moins une image numérique de la pièce 15 aéronautique acquise par le système d'imagerie numérique ; Ce système de contrôle non destructif est remarquable en ce qu'il comprend en outre un dispositif d'estimation de valeurs dites optimales de paramètres de réglage du système d'imagerie numérique, lesdites valeurs optimales de paramètres de réglage étant destinés à être appliquées audit 20 système d'imagerie numérique avant l'acquisition de ladite au moins une image numérique de la pièce aéronautique lors d'un réglage du système d'imagerie numérique, ledit dispositif d'estimation comprenant : - un premier module d'obtention configuré pour obtenir un modèle numérique caractérisant le système d'imagerie numérique, ce modèle 25 numérique comprenant un modèle caractérisant la source et un modèle caractérisant le détecteur ; - un deuxième module d'obtention, configuré pour obtenir un modèle numérique caractérisant la pièce aéronautique ; - un module d'évaluation configuré pour évaluer, pour une pluralité de 30 valeurs distinctes des paramètres de réglage, un rapport contraste sur bruit obtenu pour le système d'imagerie numérique et ladite pièce aéronautique, ce module d'évaluation utilisant les modèles numériques obtenus caractérisant le système d'imagerie numérique et la pièce aéronautique ; et 3055728 4 - un module de détermination, configuré pour déterminer des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit. La source de rayons électromagnétiques est par exemple une 5 source de rayons X. Ainsi, l'invention propose un procédé et un système permettant de déterminer de façon automatique les paramètres de réglage d'un système d'imagerie numérique de type radiographique ou tomographique destiné à être utilisé dans un processus de contrôle non destructif d'une 10 pièce aéronautique. Ces paramètres comprennent par exemple de façon classique, au moins un paramètre parmi : - une tension d'opération de la source du système d'imagerie numérique ; - un temps d'exposition de la pièce aéronautique à un faisceau de 15 rayons émis par la source du système d'imagerie numérique ; - une valeur de courant appliquée à la source du système d'imagerie numérique ; et - une épaisseur de filtre appliqué dans le système d'imagerie numérique à la source (aussi désigné par filtre externe à la source). 20 Le procédé selon l'invention s'appuie à cet effet sur divers modèles numériques du système d'imagerie numérique et de la pièce aéronautique qui permettent la détermination de paramètres de réglage optimaux : grâce à l'invention, on obtient un réglage efficace de ces systèmes tout en supprimant toute subjectivité pouvant affecter ce 25 réglage. Ces modèles numériques peuvent être par exemple prédéterminés et stockés dans des bases de données, sous forme d'un ensemble de valeurs et/ou de courbes. L'obtention de ces modèles comprend alors l'accès à ces bases de données et l'extraction des paramètres physiques modélisés pertinents de ces bases pour mettre en 30 oeuvre l'invention. Le critère retenu par l'invention pour l'optimisation des paramètres de réglage précités est le rapport contraste sur bruit qui permet de quantifier la qualité des images fournies par le système d'imagerie numérique considéré. Ce rapport est déterminé à partir des 35 différentes modélisations réalisées conformément à l'invention : modélisation du système d'imagerie numérique, et plus particulièrement 3055728 5 de ses principaux composants, à savoir de sa source de rayons électromagnétiques et du détecteur qu'il utilise pour détecter les rayons électromagnétiques émis par la source, et modélisation de la pièce aéronautique sur laquelle est réalisée le contrôle non destructif et dont on 5 souhaite acquérir des images au moyen du système d'imagerie numérique. La démarche ainsi adoptée par l'invention permet d'offrir une solution de réglage générique pouvant s'appliquer à de nombreux systèmes d'imagerie numérique conventionnels de radiographie et de 10 tomographie. En recourant à un modèle numérique du système d'imagerie numérique utilisé lors du contrôle non destructif, l'invention permet avantageusement de prendre en compte les spécificités et l'unicité de ce système. En outre, le réglage proposé par l'invention tient compte de la 15 pièce aéronautique sur laquelle on effectue un contrôle non destructif. Il est donc adapté à cette pièce et permet d'obtenir une qualité d'image optimale pour cette pièce. L'invention offre donc une solution rapide et efficace permettant d'automatiser le paramétrage de systèmes d'imagerie numérique utilisés 20 dans le cadre de contrôle non destructif. La précision de cette solution la rend applicable aisément dans un milieu industriel à n'importe quel système d'imagerie numérique de type radiographique ou tomographique, pour n'importe quelle pièce issue de l'industrie aéronautique. En outre, elle s'applique à de nombreux défauts susceptibles 25 d'être présents dans la pièce, comme notamment une sous-épaisseur ou une surépaisseur de matériau, une porosité, une inclusion, etc. On note par ailleurs que grâce à l'invention, en optimisant les paramètres de réglage du système d'imagerie numérique, on s'assure d'une utilisation optimale de ce dernier et de ses composants : l'invention permet 30 notamment, en corollaire, de réduire l'usage de la source de rayons électromagnétiques et de limiter taux de vieillissement du détecteur utilisé par le système d'imagerie numérique. Comme mentionné précédemment, l'invention s'appuie sur la modélisation des principaux composants du système d'imagerie numérique 35 et de la pièce aéronautique sur laquelle on souhaite réaliser un contrôle non destructif. Cette modélisation permet de quantifier l'impact de ces 3055728 6 éléments sur la qualité de l'image numérique fournie par le système d'imagerie numérique. Préférentiellement, on considère pour les principaux composants du système d'imagerie numérique, une modélisation physique 5 caractérisant la réponse de ces composants en énergie. Ainsi, dans un mode particulier de réalisation, le modèle caractérisant la source fournit des spectres d'émission de la source pour différentes valeurs de tension d'opération de la source. Chaque spectre d'émission fournit l'intensité des rayons émis 10 par la source (i.e. la quantité de photons générée par unité de temps ou par émission) en fonction de la longueur d'onde. Une telle modélisation permet avantageusement de tenir compte du caractère polychromatique de la source. Elle prend en considération la géométrie de la source, et peut être déterminée au moyen de simulations (ex. de type Monte-Carlo), 15 de données empiriques et/ou de d'expressions analytiques. Dans un mode particulier de réalisation, le modèle caractérisant la source fournit en outre des estimations d'au moins deux couches de demi-atténuation de la source pour différentes valeurs de tension d'opération de la source, le procédé comprenant en outre une étape de 20 validation du modèle caractérisant la source en comparent les estimations desdites au moins deux couches de demi-atténuation de la source à des résultats expérimentaux. Cette étape de validation permet de s'assurer que le modèle numérique utilisé pour caractériser la source est bien représentatif de la 25 réalité physique de la source et donc qu'il est bien adapté. On s'assure ainsi que les paramètres de réglage extraits grâce à l'invention sont bien optimaux pour le système d'imagerie considéré. Dans un mode particulier de réalisation, le modèle caractérisant le détecteur fournit une réponse spectrale du détecteur à un faisceau de 30 rayons émis par la source. Une telle réponse caractérise la façon dont le détecteur pondère les énergies des différents photons arrivant à la source, en fonction de la longueur d'onde. Comme pour la source, le modèle considéré pour le détecteur permet donc également de prendre en considération le 35 caractère polychromatique de la source. 3055728 7 Dans un mode particulier de réalisation, le procédé comprend en outre une étape de validation du modèle caractérisant le détecteur comprenant : - une estimation d'au moins un niveau de gris théorique, de 5 préférence au moins 3 niveaux, vu par le détecteur pour au moins une tension d'opération donnée de la source, en utilisant le modèle caractérisant la source et le modèle caractérisant le détecteur ; - une obtention d'une mesure d'au moins un niveau de gris, de préférence au moins 3 niveaux, vu par le détecteur pour ladite au moins 10 une tension d'opération donnée de la source ; - une comparaison dudit au moins un niveau de gris théorique estimé avec ledit au moins un niveau de gris mesuré. De cette sorte, on s'assure que le modèle numérique utilisé pour caractériser le détecteur et extraire les paramètres de réglage optimaux 15 du système d'imagerie numérique est bien représentatif de la réalité. Ceci permet de garantir la pertinence des paramètres optimaux déterminés conformément à l'invention. Comme mentionné précédemment, l'invention s'appuie sur une modélisation numérique et physique des principaux composants du 20 système numérique mais également de la pièce aéronautique sur laquelle on souhaite effectuer un contrôle non destructif. Dans un mode particulier de réalisation, le modèle numérique caractérisant la pièce aéronautique fournit des coefficients d'atténuation massique d'un matériau composant la pièce pour différentes valeurs 25 d'énergie d'un faisceau de rayons émis par la source. En effet, chaque pièce dont on souhaite acquérir une image numérique au moyen du système d'imagerie numérique pour effectuer un contrôle non destructif, a ses propres caractéristiques physiques (liées en partie à la composition du matériau qui la compose) qui influent sur le 30 système d'imagerie numérique et sur son fonctionnement. Prendre en compte la composition de la pièce aéronautique et plus particulièrement son coefficient d'atténuation massique pour différentes valeurs d'énergie des rayons émis par la source permet de régler le système d'imagerie numérique de façon éclairée et optimale pour la pièce en question. 3055728 8 Pour des matériaux non élémentaires, l'obtention du modèle numérique caractérisant la pièce aéronautique comprend par exemple, dans un mode particulier de réalisation : - une étape d'obtention d'une composition chimique du matériau 5 composant la pièce, ledit matériau étant composé d'au moins un élément chimique ; - pour chaque élément chimique composant le matériau, une étape d'obtention d'un coefficient d'atténuation massique de cet élément pour lesdites différentes valeurs d'énergie ; et 10 - une étape de calcul des coefficients d'atténuation massiques du matériau pour lesdites différentes valeurs d'énergie à partir des coefficients d'atténuation massiques des éléments composant le matériau pour ces différentes énergies et des proportions de ces éléments dans le matériau. 15 Le coefficient d'atténuation massique de chaque élément du matériau peut être obtenu par exemple en consultant une base de données NIST (National Institute of Standards and Technology). Cette procédure permet de déterminer de façon automatique le coefficient massique d'atténuation massique d'une pièce aéronautique 20 pour chaque énergie de rayons considérée. Le critère utilisé par l'invention pour extraire les paramètres de réglage optimaux du système numérique est le rapport contraste sur bruit (ou CNR pour Contrast-to-Noise Ratio). Conformément à l'invention, ce rapport est évalué pour différentes valeurs des paramètres que l'on 25 cherche à optimiser, puis les valeurs de ces paramètres conduisant au CNR optimal sont déterminées à partir des rapports ainsi évalués. Dans un mode particulier de réalisation, lors de l'étape d'évaluation, le rapport contraste sur bruit noté CNR obtenu pour le système d'imagerie numérique et ladite pièce aéronautique est évalué 30 pour différentes valeurs des paramètres de réglage selon l'expression : Sdet(L, E) - Sdet(L) OU : Sdet(L, E) + Sdet(L) CNR = .J Sdet(L) edet(L) + Sdet(L, E) CNR = edet(L) 3055728 9 où Sdet(L), respectivement Sdet(L, E), désigne une valeur de signal vue par le détecteur du système d'imagerie numérique lorsqu'une épaisseur L, respectivement une épaisseur L - E, de la pièce aéronautique est traversée par le faisceau de rayons électromagnétiques (ex. rayons X) 5 émis par la source, E désignant une épaisseur d'un indicateur recherché durant le contrôle non destructif dans la pièce aéronautique, lesdites valeurs de signal étant définies par : Sdet(L) =110(Ei,V)A(Ej, L)R(Ei) et Sdet(L,E) = I o(Ei , V)As(Ei , L, E)R(E1) où : 10 - io(Ei, V) désigne l'intensité des photons présents dans un faisceau de rayons électromagnétiques émis par la source ayant pour énergie E; correspondant à une longueur d'onde j, ladite intensité dépendant de la tension V d'opération de la source, du temps d'exposition de la pièce aéronautique au faisceau de rayons émis par la source, et de la valeur de 15 courant appliquée à la source ; - A(Ei,L) désigne un facteur d'atténuation subie par le faisceau de rayons électromagnétiques émis par la source ayant pour énergie Ei , ce facteur d'atténuation dépendant d'un coefficient d'atténuation massique du matériau composant la pièce aéronautique et de l'épaisseur L ; 20 - Ae(Ei, L, E) désigne un facteur d'atténuation subie par le faisceau de rayons électromagnétiques émis par la source ayant pour énergie E; , ce facteur d'atténuation dépendant d'un coefficient d'atténuation massique du matériau composant la pièce aéronautique, de l'épaisseur L, et de l'épaisseur E de l'indicateur ; et 25 R(Ei) désigne la réponse spectrale du détecteur en fonction de l'énergie E3. Le calcul proposé par l'invention donne une évaluation précise et fiable de la qualité des images fournies par le système numérique pour la pièce aéronautique considérée. Il en résulte l'extraction de paramètres 30 de réglage optimaux pour le système d'imagerie numérique adaptés aux 3055728 10 principaux composants du système d'imagerie numérique et à la pièce aéronautique faisant l'objet du contrôle non destructif. Il convient de noter que le modèle proposé par l'invention pour calculer le rapport contraste-sur-bruit s'applique, moyennant une 5 définition appropriée de l'épaisseur E à différents types d'indicateurs recherchés dans la pièce aéronautique lors du contrôle destructif. Par indicateurs, on entend des structures que l'on recherche dans la pièce et qui sont représentatifs notamment de la présence d'un défaut dans la pièce comme par exemple une surépaisseur (ex. surépaisseur en cavité 10 interne dans des noyaux céramiques pour l'élaboration d'aubes de turbine hate pression), une sous-épaisseur, une inclusion, une porosité, etc. En appliquant le modèle précité, les paramètres de réglage déterminés pour le système d'imagerie numérique sont également adaptés à l'indicateur recherché lors du contrôle non destructif dans la pièce aéronautique, et 15 facilite sa détection. A titre illustratif, lorsque le défaut recherché est une sous-épaisseur ou une porosité, l'épaisseur E est choisie « négative » de sorte à traduire une variation d'épaisseur négative par rapport à l'épaisseur L de la pièce aéronautique, une sous-épaisseur ou une porosité se traduisant 20 par un manque de matière au niveau de la pièce. Lorsque le défaut recherché est une surépaisseur, l'épaisseur E est au contraire choisie « positive » pour traduire la présence d'un surplus de matière. Pour un défaut de type inclusion, l'épaisseur E traduit l'épaisseur de matériau de l'inclusion. 25 Dans un mode particulier de réalisation : le facteur d'atténuation A(EJ, L) est défini par : A(Ei, L) = Afint(Ep L fint)Afext(Ep Lfext )Amat(Ej, L) ; et - le facteur d'atténuation A(Ei,L,E) est défini par : A(E1, L, = Afint(Ei, L fint)Afext(Ei, L fext)Amat (Ei , L - E)Aind(Ei, E) OÙ Afint(EpLfint), Afext(Ej, Lfext), Amat(Ej, L) et Amat(Ej, L - désignent 30 respectivement l'atténuation subie par le faisceau de rayons électromagnétiques (ex. rayons X) ayant pour énergie Ei due à un filtrage inhérent à la source, l'atténuation subie par le faisceau de rayons électromagnétiques ayant pour énergie Ei due à un filtrage externe à la source, l'atténuation subie par le faisceau de rayons électromagnétiques 35 ayant pour énergie Ei due au matériau composant la pièce aéronautique et 3055728 11 l'atténuation subie par le faisceau de rayons électromagnétiques ayant pour énergie Ei en présence de l'indicateur recherché d'épaisseur E. Le facteur d'atténuation défini dans ce mode de réalisation permet de prendre en considération de façon très complète les différentes 5 atténuations auxquelles est soumis le faisceau de rayons émis par la source lors de l'utilisation du système d'imagerie numérique. Ces atténuations sont modélisées ici par des opérations de filtrages réalisées sur le faisceau de rayons émis par la source : filtrage inhérent à la source (et au matériau composant la source), filtrage externe à la source, mais 10 aussi filtrage lié à proprement parler au matériau de la pièce aéronautique que l'on souhaite contrôler incluant le filtrage opéré par les éléments de la pièce mais également par l'indicateur que l'on recherche dans la pièce (cet indicateur pouvant être constitué de vide par exemple dans le cas d'une porosité). 15 Dans un mode particulier de réalisation, l'étape d'estimation comprend la mise en oeuvre d'une pluralité d'itérations balayant tout ou partie de la pluralité de valeurs distinctes des paramètres de réglage, et au cours desquelles, l'évaluation du rapport contraste-sur-bruit est réalisée en accédant à des valeurs des modèles numériques caractérisant la 20 source, le détecteur et la pièce aéronautique stockées préalablement dans des bases de données. Ce mode de réalisation s'appuyant sur un stockage préalable des modèles numériques dans des bases de données et sur des boucles imbriquées balayant différentes valeurs des paramètres de réglage permet 25 d'accélérer la détermination des paramètres de réglage optimaux du système d'imagerie numérique. Dans un mode particulier de réalisation, certaines étapes du procédé de contrôle non destructif, et en particulier l'étape d'estimation des paramètres de réglage, sont déterminées par des instructions de 30 programmes d'ordinateurs. En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en oeuvre dans un dispositif d'estimation de paramètres de réglage d'un système d'imagerie numérique ou plus 35 généralement dans un ordinateur, ce programme comportant des 3055728 12 instructions adaptées à la mise en oeuvre de l'étape d'estimation d'un procédé de contrôle non destructif tel que décrit ci-dessus. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de 5 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. L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur 10 tel que mentionné ci-dessus. Le support d'informations peut être n'importe quelle entité ou 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 15 d'enregistrement magnétique, par exemple une disquette (floppy disc) ou un disque dur. 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 20 moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. Alternativement, le support d'informations peut être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé dans l'exécution du procédé en 25 question. On peut également envisager, dans d'autres modes de réalisation, que le procédé et le système de contrôle non destructif selon l'invention présentent en combinaison tout ou partie des caractéristiques précitées. 30 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 35 caractère limitatif. Sur les figures : 3055728 13 - la figure 1 représente, dans son environnement, un système de contrôle non destructif conforme à l'invention comprenant un système d'imagerie numérique, un dispositif de contrôle non destructif et un dispositif d'estimation de paramètres de réglage du système d'imagerie 5 numérique ; - la figure 2 représente, sous forme schématique, l'architecture matérielle du dispositif d'estimation de paramètres de réglage de la figure 1; - la figure 3 représente, sous forme d'ordinogramme, les principales 10 étapes d'un procédé de contrôle non destructif selon l'invention, dans un mode particulier de réalisation dans lequel il est mis en oeuvre par le système de contrôle non destructif de la figure 1 ; - les figures 4A, 4B et 4C illustrent respectivement les principales étapes mises en oeuvre par le dispositif d'estimation de paramètres de 15 réglage de la figure 1 pour obtenir un modèle numérique de la source et un modèle numérique du détecteur du système d'imagerie numérique de la figure 1, ainsi qu'un modèle numérique d'une pièce aéronautique sur lequel opère le système de contrôle non destructif de la figure 1 ; et - la figure 5 représente, sous forme d'ordinogramme, les principales 20 étapes mises en oeuvre lors de l'étape d'estimation des paramètres de réglage optimaux mise en oeuvre dans le procédé de contrôle non destructif illustré à la figure 3. Description détaillée de l'invention 25 La figure 1 représente, dans son environnement, un système 1 de contrôle non destructif conforme à l'invention, dans un mode particulier de réalisation. Dans l'exemple envisagé ici, le système 1 permet de procéder au contrôle non destructif de pièces aéronautiques telles que par exemple une aube fan 2 de turbine. Toutefois, aucune limitation n'est 30 attachée à la nature de la pièce sur laquelle le contrôle non destructif est appliquée. 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, de fusée, etc. Conformément à l'invention, le système 1 de contrôle non 35 destructif comprend : 3055728 14 - un système d'imagerie numérique 3, apte à acquérir et à fournir une ou plusieurs images numériques IM de la pièce aéronautique 2 ; et - un dispositif 4 de contrôle non destructif de la pièce aéronautique 2 configuré pour utiliser la ou les images numériques IM de la pièce 5 aéronautique acquise(s) et fournie(s) par le système d'imagerie numérique 3 et former un diagnostic d'intégrité DIAG de la pièce 2. Ce diagnostic porte notamment sur un ou plusieurs indicateurs IND ou paramètres recherchés dans la pièce. Aucune limitation n'est attachée à la nature de ces indicateurs ; il peut s'agir par exemple d'indicateurs 10 représentatifs de surépaisseur dans la pièce, de sous-épaisseur, ou de la présence de porosités ou d'inclusions, etc. Un tel dispositif de contrôle non destructif est connu en soi et n'est pas décrit plus en détail ici. Le système d'imagerie numérique 3 est ici un système de radiographie ou de tomographie numérique. Il comprend, de façon 15 connue, une source 3A de rayons électromagnétiques (par exemple une source émettant des faisceaux de rayons X) et un détecteur 3B apte à détecter les rayons électromagnétiques émis par la source 3A. La source 3A applique ici le principe du bombardement d'électrons sur une cible métallique dans un tube à rayons X : les 20 électrons sont extraits d'un filament métallique (la cathode) puis accélérés par une tension électrique importante appliquée à la source 3A (tension d'opération au sens de l'invention) dans un tube sous vide. Ce faisceau d'électrons émis par la source 3A est focalisé de manière à bombarder une cible métallique (l'anode), par exemple en tungstène. Les interactions 25 entre les électrons et la matière lorsque ceux-ci heurtent la cible produisent des photons de rayons « X ». Un tel fonctionnement est connu en soi et n'est pas décrit plus en détail ici. Le détecteur 3B du système d'imagerie numérique 3 fonctionne de manière indirecte : un scintillateur absorbe les photons des rayons X 30 produits par la source 3A, et par des mécanismes de fluorescence, émet une lumière visible. Les photons visibles émis sont dirigés sur un réseau de photodiodes, chaque élément du réseau de photodiodes étant associé à un pixel. Les photons visibles sont convertis en charge électrique au niveau du réseau de photodiodes puis en signal électrique par un 35 photomultiplicateur. Le signal électrique véhicule des niveaux de gris associés à chaque pixel du détecteur et reflétant l'absorption des photons 305 5 72 8 15 au niveau du pixel considéré. Un tel fonctionnement est connu en soi et n'est pas décrit plus en détail ici. Pour fournir des images numériques IM de la pièce aéronautique 2 permettant la mise en oeuvre d'un contrôle non destructif 5 efficace par le dispositif 4 de contrôle non destructif, il est classique d'apporter certains réglages au système d'imagerie numérique 3. Ces réglages concernent en particulier différents paramètres du système d'imagerie numérique 3, comprenant dans le mode de réalisation décrit ici : 10 - la tension V d'opération de la source 3A ; - le temps d'exposition de la pièce aéronautique 2 au faisceau de rayons émis par la source 3A ; - la valeur de courant i appliquée à la source 3A ; et - une épaisseur de filtre appliqué dans le système d'imagerie 15 numérique à la source 3A (aussi désigné par filtre externe à la source). Cette liste n'est toutefois pas exhaustive et d'autres paramètres de réglage du système d'imagerie numérique 3 peuvent être envisagés ; de même, seulement une partie des paramètres précités peut être envisagée. 20 Conformément à l'invention, le système 1 de contrôle non destructif comporte de façon avantageuse un dispositif 5 d'estimation automatique des paramètres de réglage précités. Ce dispositif 5 d'estimation est apte à estimer des valeurs optimales de ces paramètres de réglage destinées à être appliquées au système 3 d'imagerie 25 numérique lors d'un réglage de celui-ci réalisé avant l'acquisition des images numériques IM de la pièce aéronautique 2. Dans le mode de réalisation décrit ici, le dispositif 5 d'estimation a l'architecture matérielle d'un ordinateur telle que représentée schématiquement à la figure 2. 30 Il comprend, notamment un processeur 6, une mémoire vive 7, une mémoire morte 8, une mémoire flash non volatile 9 ainsi que des moyens d'entrée/sortie 10 permettant à un opérateur du système 3 d'imagerie numérique d'interagir avec le dispositif 5 d'estimation. Ces moyens d'entrée/sortie 10 comprennent par exemple un écran, un clavier, 35 une interface homme-machine permettant à l'opérateur du système 3 notamment de saisir des données dans le dispositif 5 d'estimation et 3055728 16 d'obtenir les paramètres de réglage optimaux du système 3 d'imagerie numérique déterminés par celui-ci conformément à l'invention. La mémoire morte 8 du dispositif d'estimation 5 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 6 5 et sur lequel est enregistré un programme d'ordinateur conforme à l'invention. Le programme d'ordinateur PROG définit des modules fonctionnels et logiciels ici, configurés pour mettre en oeuvre une estimation des paramètres de réglage optimaux du système 3 d'imagerie 10 numérique conformément à l'invention. Ces modules fonctionnels s'appuient sur et/ou commandent les éléments matériels 6-10 du dispositif d'estimation 5 cités précédemment. Ils comprennent notamment ici, comme illustré sur la figure 1 : - un premier module d'obtention 5A configuré pour obtenir un modèle 15 numérique MOD3 caractérisant le système d'imagerie numérique 3, ce modèle numérique comprenant un modèle MOD3A caractérisant la source 3A et un modèle MOD3B caractérisant le détecteur 3B. Le premier module d'obtention 5A est configuré ici pour extraire ces modèles numériques MOD3A et MOD3B des bases de données DB3A et DB3B du dispositif 5 20 d'estimation, dans lesquels ils ont été préalablement stockés après leur détermination. Par base de données, on entend ici au sens large tout ensemble de données permettant de stocker des informations ; - un deuxième module d'obtention 5B, configuré pour obtenir un modèle numérique MOD2 caractérisant la pièce aéronautique 2 ; le 25 deuxième module d'obtention 5B est configuré ici pour extraire le modèle numérique MOD2 d'une base de donnée DB2 du dispositif 5 d'estimation dans lequel il a été préalablement stocké après sa détermination ; - un module d'évaluation 5C configuré pour évaluer, pour une pluralité de valeurs distinctes des paramètres de réglage, un rapport contraste sur 30 bruit obtenu pour le système d'imagerie numérique 3 et la pièce aéronautique 2, ce module d'évaluation 5C utilisant les modèles numériques MOD3 (incluant les modèles MOD3A et MOD3B) et MOD2 obtenus par les modules 5A et 5B et caractérisant respectivement le système d'imagerie numérique 3 et la pièce aéronautique 2 ; et 3055728 17 - un module de détermination 5D, configuré pour déterminer des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit. Les fonctions des différents modules sont décrites plus en détail 5 en référence à l'étape d'estimation du procédé de contrôle non destructif selon l'invention mise en oeuvre par le dispositif 5 d'estimation du système 1 de contrôle non destructif de la figure 1. On note que deux modules d'obtention distincts 5A et 5B sont décrits ici, toutefois il est possible que les fonctions mises en oeuvre par 10 ces deux modules soient mises en oeuvre dans un autre mode de réalisation par un module unique. Pour mieux comprendre les avantages apportés par l'invention, nous allons maintenant décrire brièvement les principales étapes du procédé de contrôle non destructif selon l'invention. La figure 3 illustre 15 ces étapes telles qu'elles sont mises en oeuvre par le système 1 de contrôle non destructif de la figure 1. Conformément à l'invention, on suppose ici que l'on souhaite réaliser un contrôle non destructif de la pièce aéronautique 2 à partir d'une ou plusieurs images numériques de cette pièce acquises au moyen 20 du système 3 d'imagerie numérique (système radiographique ou tomographique). A cet effet, un réglage préalable de divers paramètres du système 3 d'imagerie numérique est effectué (étape D10). Dans le mode de réalisation décrit ici, ces paramètres comprennent, comme indiqué 25 précédemment : - la tension V d'opération de la source 3A ; - le temps d'exposition de la pièce aéronautique 2 au faisceau de rayons émis par la source 3A ; - la valeur de courant i appliquée à la source 3A ; et 30 - l'épaisseur de filtre appliqué dans le système d'imagerie numérique à la source. Conformément à l'invention, le réglage de ces paramètres est effectué en deux temps : - dans un premier temps, des valeurs optimales de ces paramètres sont 35 estimées automatiquement par le dispositif 5 d'estimation (étape D10- 1) , 3055728 18 - puis dans un second temps, les valeurs optimales ainsi estimées sont appliquées au système 3 d'imagerie numérique pour le paramétrer (étape D10-2). Ce réglage du système 3 d'imagerie numérique peut être effectué par exemple manuellement par un opérateur du système 5 3 au moyen des valeurs optimales fournies par le dispositif 5 d'estimation. Puis une ou plusieurs images numériques IM de la pièce aéronautique 2 sont acquises par le système 3 d'imagerie numérique ainsi réglé (étape D20), et fournies au dispositif 4 de contrôle non destructif. 10 Un contrôle non destructif de la pièce aéronautique 2 est alors réalisé via le dispositif 4 de contrôle non destructif à partir des images numériques IM (étape D30). Une telle étape peut être réalisée de diverses manières, de façon connue en soi, et n'est pas décrite en détail ici. Les images ainsi obtenues permettent un accès non destructif à la pièce à 15 contrôler. Elles peuvent par exemple être comparées par des méthodes de traitement d'images connues à des images acquises par le même procédé pour une pièce saine dans le cadre de recherche d'anomalie(s) dans la pièce à contrôler. Conformément à l'invention, l'estimation réalisée à l'étape D10-1 20 des paramètres de réglage optimaux du système 3 d'imagerie numérique pour son utilisation sur la pièce aéronautique 2 s'appuie sur divers modèles numériques caractérisant physiquement le système 3 d'imagerie numérique, et notamment ses principaux composants à savoir la source 3A et le détecteur 3B, et la pièce aéronautique 2. Ces modèles 25 numériques sont obtenus par les modules 5A et 5B respectivement du dispositif 5 d'estimation. Ils peuvent être pré-calculés ou obtenus à la volée lors de l'estimation des paramètres de réglage optimaux. Toutefois, dans le mode de réalisation décrit ici, ces modèles numériques sont pré-calculés et stockés dans des bases de données (DB3A, DB3B et DB2) du 30 dispositif d'estimation 5 pour une utilisation future. Cela permet avantageusement d'accélérer la recherche des paramètres de réglage optimaux lorsque l'on souhaite utiliser le système d'imagerie numérique 3. Nous allons dans un premier temps décrire comment ces modèles numériques sont obtenus (i.e. construits) et mis à disposition des 35 modules 5A et 5B du dispositif 5 d'estimation dans les bases de données DB3A, DB3B et DB2. 3055728 19 La figure 4A représente les principales étapes mises en oeuvre pour obtenir le modèle MOD3A caractérisant la source 3A du système d'imagerie numérique 3. On note que ces étapes peuvent être mises en oeuvre par le module 5A du dispositif d'estimation 5 ou via tout autre 5 module externe au dispositif d'estimation 5 et apte à stocker le modèle MOD3A obtenu dans la base de données DB3A. Le modèle MOD3A fournit ici des spectres d'émission de la source 3A pour différentes valeurs de tension d'opération de la source notées V1, 1/2, ...,17N où N désigne un entier supérieur à 1. Dans l'exemple 10 envisagé ici, les tensions d'opération V1, V2, ..., VN varient de 0 à 450 kV, par pas de 1kV. Dans le mode de réalisation décrit ici, le modèle MOD3A fournit en outre des spectres d'émission de la source 3A pour N valeurs de tension d'opération de la source et pour M valeurs de l'angle d'inclinaison 15 par rapport à l'horizontal de l'anode de la source 3A, M désignant un entier supérieur à 1. On note en effet que l'angle d'inclinaison de l'anode de la source influe sur l'angle d'émission de la source 3A, ce qui a une importance sur la qualité du spectre (effet Heel). Par spectre d'émission, on entend l'intensité (i.e. le nombre de 20 photons émis par unité de temps ou par émission) des faisceaux de rayons émis par la source 3A en fonction de la longueur d'onde de ces rayons. Une telle modélisation permet avantageusement de prendre en considération le caractère polychromatique de la source 3A. On note que la longueur d'onde est reliée de manière univoque à l'énergie des 25 faisceaux de rayons considérés. Autrement dit à une longueur d'onde indexé par un entier j et notée 2 correspond une unique énergie de faisceaux Ei. Pour déterminer les spectres d'émission du modèle MOD3A correspondant aux N valeurs possibles de la tension d'opération de la 30 source et M valeurs possibles de l'angle d'inclinaison de l'anode (étapes d'initialisation E10, de test E30 et E35 et d'incrémentation E40 et E45), on utilise, dans le mode de réalisation décrit ici, un logiciel de simulation permettant de modéliser le comportement de la source 3A. Ce logiciel est par exemple ici le logiciel SpekCalc, développé par V.N. Hansen et al., et 35 décrit notamment dans le document intitulé « An efficient Monte Carlobased algorithm for scatter correction in kev, cone-beam et », Physics in 3055728 20 Medicine and Biology, 54 : 3847-3864, 2009. Il utilise des processus d'interaction des électrons au niveau de l'anode, ainsi que des données issues de simulations de type Monte-Carlo et/ou d'expressions empiriques et/ou analytiques permettant d'estimer les spectres d'émission. On note 5 que le logiciel SpekCalc permet de modéliser des spectres d'émission pour des angles d'inclinaison de l'anode de la source compris entre 1 et 90 degrés, toutefois les sources de rayons électromagnétiques et en particulier de rayons X possèdent rarement des angles d'inclinaison allant au-delà de 45 degrés. Par exemple ici, on considère un angle d'inclinaison 10 de l'anode de la source variant de 5 à 45 degrés par incrément de 1 degré. Les spectres d'émission correspondant aux N valeurs de tensions d'opération et aux M valeurs d'angles d'inclinaison obtenus grâce au logiciel SpekCalc (étape E20) sont stockés dans la base de données 15 DB3A (étape E50). Dans le mode de réalisation décrit ici, pour accéder à l'un des spectres d'émission stockés dans la base de données DB3A, le module 5A fournit à la base de données DB3A (ou de manière équivalente interroge au moyen de) : - l'angle d'inclinaison noté al de l'anode de la source 3A par rapport à 20 l'horizontal ; - l'énergie maximum Emax(V1) (en EV ou kEV) correspondant à la tension Vi d'opération de la source 3A considérée pour laquelle on souhaite estimer le spectre d'émission de la source. Le module 5A obtient ainsi le spectre d'émission de la source 3A 25 pour la tension d'opération Vi et l'angle al, noté dans la suite de la description ios(Vi, cri). On note que chaque spectre d'émission /0s(Vi, ai) fournit le nombre de photons d'énergie Ej par unité de courant et par unité de temps pour différentes énergies j=1,...,K (ou de façon équivalente, pour différentes longueurs d'onde j=1,...,K), K désignant 30 un entier supérieur à 1. Dans le mode de réalisation décrit ici, la base de données D3A comprend également des estimations d'au moins deux couches de demi-atténuation de la source 3A pour deux matériaux distincts, à savoir le cuivre (Cu) et l'aluminium (AI), pour chaque tension d'opération Vi et pour 35 chaque angle d'inclinaison ai. Ces estimations ont été fournies par le 305 5 72 8 21 logiciel SpekCalc au cours de l'étape E20. Plus précisément, la base de données D3A comprend : - les couches de demi-atténuation HVL1Cu(Vi, ai) et HVL1A1(Vi, ai) désignant l'épaisseur de matière nécessaire (cuivre ou aluminium ici) 5 pour diminuer de moitié la radiation du faisceau à 1m de la fenêtre de la source ; et - les couches de demi-atténuation HVL2Cu(vi, al) et HVL2A1(Vi, al) désignant l'épaisseur de matière nécessaire (cuivre ou aluminium ici) pour diminuer de moitié la radiation du faisceau résultant 10 respectivement de l'utilisation des couches de demi-atténuation HVL1Cu(Vi, al) et HVL1AI(Vi, al) à 1m de la fenêtre de la source. Le module 5A accède ici à ces estimations en même temps qu'au spectre d'émission de la source 3A pour la tension d'opération Vi et l'angle al. 15 Dans le mode de réalisation décrit ici, on envisage une validation du modèle numérique MOD3A, autrement dit des spectres d'émission stockés dans la base de données D3A (étape E60). Cette validation peut être mise en oeuvre par exemple par le module 5A ou par un autre module externe au dispositif d'estimation 5, en accédant aux 20 spectres d'émission stockés dans la base de données D3A et aux estimations des couches de demi-atténuation associées. Dans la suite de la description, par souci de simplification, on désignera le module chargé de cette validation de façon générale par module de validation. Cette validation vise ici à s'assurer que les spectres générés 25 pour la source 3A sont bien représentatifs de la réalité physique de la source 3A. Cela contribue à permettre au dispositif d'estimation 5 d'obtenir grâce à l'invention des valeurs optimales des paramètres de réglage. A cet effet, le module de validation compare ici, pour au moins une sélection de spectres comprenant au moins deux spectres générés au 30 cours de l'étape E20 pour deux tensions d'opération différentes de la source 3A, les estimations des couches de demi-atténuation HVL1Cu, HVL1AI, HVL2Cu et HVL2AI fournies lors de l'étape E20 pour ces tensions avec des résultats expérimentaux. Les valeurs expérimentales des couches de demi-atténuation 35 HVL1Cu, HVL1AI, HVL2Cu et HVL2AI peuvent être générées préalablement, au moyen de couches de cuivre et d'aluminium de 3055728 22 différentes épaisseurs et d'un dosimètre, puis stockées dans la base de données DB3A pour être comparées avec les valeurs des couches de demi atténuation accompagnant les spectres dans la base de données DB3A. Si les valeurs HVL1Cu, HVL1AI, HVL2Cu et HVL2AI 5 correspondant à chaque spectre de la sélection de spectres considérée de la base de données correspondent aux valeurs expérimentales, alors cela signifie que les spectres générés à l'étape E20 et stockés dans la base de données DB3A sont valides. Ils constituent un modèle numérique MOD3A validé caractérisant la source 3A au sens de l'invention. On note qu'en 10 pratique il est suffisamment de procéder à cette validation sur une sélection de spectres correspondant à une plage prédéterminée de valeurs de tensions et de courants. On note également que dans le mode de réalisation décrit ici, si les valeurs HVL1Cu, HVL1AI, HVL2Cu et HVL2AI d'un spectre de la base de 15 données DB3A ne correspondent pas aux valeurs expérimentales, celui-ci est conservé dans la base de données DB3A. Toutefois le module de validation cherchera par ailleurs dans la base de données DB3A, le spectre ayant des valeurs HVL1 et HVL2 qui correspondent le mieux aux valeurs de HVLA et HVL2 trouvées expérimentalement. 20 Nous allons maintenant décrire, en référence à la figure 4B, les principales étapes mises en oeuvre pour obtenir le modèle MOD3B caractérisant le détecteur 3B du système d'imagerie numérique 3. On note que ces étapes peuvent être mises en oeuvre par le module 5A du dispositif d'estimation 5 ou via tout autre module externe au dispositif 25 d'estimation 5 et apte à stocker le modèle MOD3B obtenu dans la base de données DB3B. Le modèle MOD3B fournit ici une réponse spectrale notée R du détecteur à un faisceau de rayons émis par la source d'énergie Ej, j=1,...,K. En d'autres mots, le modèle MOD3B caractérise comment le 30 détecteur 3B pondère les différentes énergies de photons arrivant de la source. On note que cette réponse spectrale dépend en grande partie de la capacité du scintillateur du détecteur 3B à absorber les différentes énergies de photons. Par conséquent, dans le mode de réalisation décrit ici, pour obtenir le modèle MOD3B, on simule préalablement, à l'aide d'un 35 logiciel de simulation approprié, les courbes dites de réponse montrant la quantité d'énergie absorbée dans le scintillateur du détecteur 3B en 305 5 72 8 23 fonction de l'énergie des photons incidents (étape F10), en fonction de l'épaisseur et type du scintillateur du détecteur 3B. L'élaboration de telles courbes de réponse s'appuie ici sur la simulation de la dose déposée par le faisceau de photons sur le 5 scintillateur, c'est-à-dire de l'énergie absorbée par le scintillateur. A cette fin, un nombre K de sources « élémentaires » de photons indépendantes, de même intensité chacune et monochromatiques sont modélisées. Chaque source indexée par j, j=1,...,K, émet un faisceau de photons d'énergie Ei avec Ej REVENDICATIONS1. Procédé de contrôle non destructif réalisé sur une pièce aéronautique (2) comprenant : - une étape de réglage (D10) d'un système d'imagerie numérique, ledit système (3) d'imagerie numérique étant un système de radiographie ou de tomographie numérique comprenant une source (3A) de rayons électromagnétiques et un détecteur (3B) apte à détecter les rayons électromagnétiques émis par la source ; - une étape d'acquisition (D20) d'au moins une image numérique (IM) de la pièce aéronautique au moyen du système de radiographie ou de tomographie numérique réglé ; et - une étape de contrôle non destructif (D30) réalisé sur la pièce aéronautique en utilisant ladite au moins une image numérique acquise ; ledit procédé étant caractérisé en ce qu'il comprend en outre une étape préalable d'estimation (D10-1) de paramètres de réglage du système d'imagerie numérique, cette étape d'estimation comprenant : - une étape d'obtention (H20) d'un modèle numérique caractérisant le système d'imagerie numérique, ledit modèle numérique comprenant un 20 modèle (MOD3A) caractérisant la source et un modèle (MOD3B) caractérisant le détecteur ; - une étape d'obtention (H20) d'un modèle numérique (MOD2) caractérisant la pièce aéronautique ; - une étape d'évaluation (H30), pour une pluralité de valeurs distinctes 25 des paramètres de réglage, d'un rapport contraste sur bruit obtenu pour le système d'imagerie numérique et ladite pièce aéronautique, cette étape d'évaluation utilisant les modèles numériques obtenus caractérisant le système d'imagerie numérique et la pièce aéronautique ; et - une étape de détermination (H70) de valeurs dites optimales des 30 paramètres de réglage optimisant le rapport contraste sur bruit, ces valeurs optimales étant appliquées au système d'imagerie numérique lors de l'étape de réglage. 2. Procédé selon la revendication 1 dans lequel les paramètres 35 de réglage du système comprennent au moins un paramètre parmi : 305 5 72 8 34 - une tension d'opération de la source du système d'imagerie numérique ; - un temps d'exposition de la pièce aéronautique à un faisceau de rayons émis par la source du système d'imagerie numérique ; 5 - une valeur de courant appliquée à la source du système d'imagerie numérique ; et - une épaisseur de filtre appliqué dans le système d'imagerie numérique à la source. 10 3. Procédé selon la revendication 1 ou 2 dans lequel le modèle caractérisant la source (MOD3A) fournit des spectres d'émission de la source pour différentes valeurs de tension d'opération de la source. 4. Procédé selon la revendication 3 dans lequel le modèle 15 (MOD3A) caractérisant la source fournit en outre des estimations d'au moins deux couches de demi-atténuation de la source pour différentes valeurs de tension d'opération de la source, ledit procédé comprenant une étape de validation (E50) du modèle caractérisant la source en comparent les estimations desdites au moins deux couches de demi-atténuation de la source à des résultats expérimentaux. 5. Procédé selon l'une quelconque des revendications 1 à 4 dans lequel le modèle caractérisant le détecteur (MOD3B) fournit une réponse spectrale du détecteur à un faisceau de rayons émis par la source. 6. Procédé selon la revendication 5 comprenant en outre une étape de validation (F30) du modèle caractérisant le détecteur comprenant : - une estimation d'au moins un niveau de gris théorique vu par le 30 détecteur pour au moins une tension d'opération donnée de la source, en utilisant le modèle caractérisant la source et le modèle caractérisant le détecteur ; - une obtention d'une mesure d'au moins un niveau de gris vu par le détecteur pour ladite au moins une tension d'opération donnée de la 35 source ; 305 5 72 8 35 - une comparaison dudit au moins un niveau de gris théorique estimé avec ledit au moins un niveau de gris mesuré. 7. Procédé selon l'une quelconque des revendications 1 à 6 dans 5 lequel le modèle numérique (MOD2) caractérisant la pièce aéronautique (2) fournit des coefficients d'atténuation massique d'un matériau composant la pièce aéronautique pour différentes valeurs d'énergie d'un faisceau de rayons émis par la source. 10 8. Procédé selon la revendication 7 dans lequel l'obtention du modèle numérique caractérisant la pièce aéronautique comprend : - une étape d'obtention (G10) d'une composition chimique du matériau composant la pièce, ledit matériau étant composé d'au moins un élément chimique ; 15 - pour chaque élément chimique composant le matériau, une étape d'obtention (G20) d'un coefficient d'atténuation massique de cet élément pour lesdites différentes valeurs d'énergie ; et - une étape de calcul (G30) des coefficients d'atténuation massiques du matériau pour lesdites différentes valeurs d'énergie à partir des 20 coefficients d'atténuation massiques des éléments composant le matériau pour ces différentes énergies et des proportions de ces éléments dans le matériau. 9. Procédé selon la revendication 8 dans lequel le coefficient 25 d'atténuation massique de chaque élément du matériau est obtenu en consultant une base de données NIST. 10. Procédé selon l'une quelconque des revendications 1 à 9 dans lequel lors de l'étape d'évaluation (H30), le rapport contraste sur bruit noté CNR obtenu pour le système d'imagerie numérique et ladite pièce aéronautique est évalué pour différentes valeurs des paramètres de réglage selon l'expression : CNR = sdee, E) - Sdet(L) ,/sdet(L) OU : 3055728 36 Sdet(L, E) + Sdet(L) CNR = \ Sdet(L) + Sdet(L, E) où Sdet(L), respectivement Sdet(L,E), désigne une valeur de signal vue par le détecteur du système d'imagerie numérique lorsqu'une épaisseur L, respectivement une épaisseur L - E, de la pièce aéronautique est traversée par le faisceau de rayons émis par la source, E désignant une 5 épaisseur d'un indicateur recherché durant le contrôle non destructif dans la pièce aéronautique, lesdites valeurs de signal étant définies par : S det(L) = 10(E1,V)A(Ei, L)R (Ei) et Sdet-E) = I 0(Ei , V)A,(Ei , L, E)R(Ei) où : - 10 (E;, V) désigne l'intensité des photons présents dans un faisceau de 10 rayons d'énergie Ei émis par la source correspondant à une longueur d'onde ladite intensité dépendant de la tension V d'opération de la source, du temps d'exposition de la pièce aéronautique au faisceau de rayons émis par la source, et de la valeur de courant appliquée à la source ; 15 - A(Ei,L) désigne un facteur d'atténuation subie par le faisceau de rayons émis par la source d'énergie Ei , ce facteur d'atténuation dépendant d'un coefficient d'atténuation massique du matériau composant la pièce aéronautique et de l'épaisseur L ; - AE(Ei, L, E) désigne un facteur d'atténuation subie par le faisceau de 20 rayons émis par la source d'énergie Ei , ce facteur d'atténuation dépendant d'un coefficient d'atténuation massique du matériau composant la pièce aéronautique, de l'épaisseur L, et de l'épaisseur E de l'indicateur ; et - R(Ei) désigne la réponse spectrale du détecteur en fonction de 25 l'énergie 11. Procédé selon la revendication 10 dans lequel : 3055728 37 - le facteur d'atténuation A(Ei, L) est défini par : A(Ei, L) = Afint(Ei,L fint)Af't(Ei,Lf't )Amat(Ej, ; et - le facteur d'atténuation A(Ei, L, E) est défini par : A(Ei, L, = Af int(Ej, Lf int)Afext(Ej, Lfext)Amat(Ei, L - E)Aind(Ei, E) OÙ Afint(Ei, L f int), Afext(Ei, L fext), Amat(E , L) et Amat(Ei, L - E) désignent 5 respectivement l'atténuation subie par le faisceau de rayons émis par la source d'énergie E j due à un filtrage inhérent à la source, l'atténuation subie par le faisceau de rayons émis par la source d'énergie Ei due à un filtrage externe à la source, l'atténuation subie par le faisceau de rayons émis par la source d'énergie Ei due au matériau composant la pièce 10 aéronautique et l'atténuation subie par le faisceau de rayons d'énergie Ei émis par la source en présence de l'indicateur recherché d'épaisseur E. 12. Procédé selon l'une quelconque des revendications 1 à 11 dans lequel l'étape d'estimation comprend la mise en oeuvre d'une pluralité d'itérations (iter) balayant tout ou partie de la pluralité de valeurs 15 distinctes des paramètres de réglage, et au cours desquelles, l'évaluation du rapport contraste-sur-bruit est réalisée en accédant à des valeurs des modèles numériques caractérisant la source, le détecteur et la pièce aéronautique stockées préalablement dans des bases de données. 20 13. Système (1) de contrôle non destructif d'une pièce aéronautique comprenant : - un système (3) d'imagerie numérique, ce système d'imagerie numérique étant un système de radiographie ou de tomographie numérique comprenant une source de rayons électromagnétiques et un 25 détecteur apte à détecter les rayons électromagnétiques émis par la source ; et - un dispositif (4) de contrôle non destructif de la pièce aéronautique configuré pour utiliser au moins une image numérique de la pièce aéronautique acquise par le système d'imagerie numérique ; 30 le système étant caractérisé en ce qu'il comprend en outre un dispositif d'estimation (5) de valeurs dites optimales de paramètres de réglage du système d'imagerie numérique, lesdites valeurs optimales de paramètres de réglage étant destinées à être appliquées audit système d'imagerie numérique avant l'acquisition de ladite au moins une image numérique de 305 5 72 8 38 la pièce aéronautique lors d'un réglage du système d'imagerie numérique, ledit dispositif d'estimation comprenant : - un premier module d'obtention (5A) configuré pour obtenir un modèle numérique caractérisant le système d'imagerie numérique, ce 5 modèle numérique comprenant un modèle caractérisant la source et un modèle caractérisant le détecteur ; - un deuxième module d'obtention (5B), configuré pour obtenir un modèle numérique caractérisant la pièce aéronautique ; - un module d'évaluation (5C) configuré pour évaluer, pour une 10 pluralité de valeurs distinctes des paramètres de réglage, un rapport contraste sur bruit obtenu pour le système d'imagerie numérique et ladite pièce aéronautique, ce module d'évaluation utilisant les modèles numériques obtenus caractérisant le système d'imagerie numérique et la pièce aéronautique ; et 15 - un module de détermination (5D), configuré pour déterminer des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit.