EP3507593A1 - Non-destructive inspection method and system carried out on an aeronautical part - Google Patents

Non-destructive inspection method and system carried out on an aeronautical part

Info

Publication number
EP3507593A1
EP3507593A1 EP17771806.1A EP17771806A EP3507593A1 EP 3507593 A1 EP3507593 A1 EP 3507593A1 EP 17771806 A EP17771806 A EP 17771806A EP 3507593 A1 EP3507593 A1 EP 3507593A1
Authority
EP
European Patent Office
Prior art keywords
source
imaging system
digital imaging
digital
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
EP17771806.1A
Other languages
German (de)
French (fr)
Inventor
Edward Romero
Clément REMACHA
Maximilian MELISSAS
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Safran SA
Original Assignee
Safran SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Safran SA filed Critical Safran SA
Publication of EP3507593A1 publication Critical patent/EP3507593A1/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60—Testing or inspecting aircraft components or systems

Landscapes

  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Health & Medical Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The non-destructive inspection method is carried out on an aeronautical part by means of a digital imaging system comprising a source of electromagnetic rays and a detector, and comprises a prior step of estimating adjustment parameters for adjusting the digital imaging system, comprising: — a step (H20) of obtaining a digital model characterising the digital imaging system, and comprising models (MOD3A, MOD3B) characterising the source and the detector; — a step (H20) of obtaining a digital model (MOD2) characterising the aeronautical part; — a step (H30) of evaluating, for a plurality of distinct values of the adjustment parameters, a contrast-to-noise ratio for the digital imaging system and the aeronautical part, by using the numerical model characterising the digital imaging system and the aeronautical part; and — a step (H70) of automatically determining so-called optimum values of adjustment parameters optimising the contrast-to-noise ratio and applied to the digital imaging system during an adjustment step.

Description

Procédé et système de contrôle non destructif réalisé sur une pièce aéronautique  Non destructive testing method and system performed on an aeronautical part
Arrière-plan de l'invention Background of the invention
L'invention concerne le domaine général de l'aéronautique.  The invention relates to the general field of aeronautics.
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.  It relates more particularly to the non-destructive testing of aeronautical parts, such as parts fitted to aircraft engines such as high pressure or low pressure turbine blades, distributors, etc. However, no limitation is attached to the type of aeronautical part considered or the material in which the piece is composed. It may be for example a composite material or not.
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.  In a known manner, non-destructive testing means a set of methods that make it possible to characterize the state of integrity and / or the quality of structures or materials without degrading them. Non-destructive testing has a preferred but non-limiting application in the field of aeronautics, and more generally in any field in which the structures whose state or quality are to be characterized are expensive and / or their reliability. operation is critical. Non-destructive testing can be advantageously performed on the structure or material considered both during production and during use or 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 dinclusions 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.  Among the existing non-destructive testing methods, some rely on digital images provided by X-ray or digital tomography systems. The advantage of these images is that they provide information directly exploitable on the inside of the structures or materials, and thus make it possible to detect internal defects that may affect these structures or materials, such as for example the presence of inclusions or blowholes. In the remainder of the description, the term "part" is generally used to designate the structure or material on which non-destructive testing is envisaged.
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 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 ia mesure de facteurs de qualité des images (ou IQI pour Indicateurs de Qualité d'Images) réalisée sur une grande gamme de conditions à partir d'images numériques en deux ou en trois dimensions issues des systèmes d'imagerie numérique considérés. L'optimisation « manuelle » 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. Many variables are involved in tuning the aforementioned digital imaging systems to acquire digital images of parts on which non-destructive testing is contemplated. Such variables include the magnification, the exposure time of the part to electromagnetic radiation (eg radii X), the size of the pixels, the quantum detection efficiency of the sensor equipping the digital imaging system, etc. Today, these settings are defined and validated, for different measurement ranges, by certification committees called COSAC (Sectorial Aerospace Certification Committee) 2 and 3. They are based largely on the experience of professionals who handle digital imaging systems and the measurement of image quality factors (or IQIs for Image Quality Indicators) performed over a wide range of conditions from two-dimensional or three-dimensional digital images from digital imaging systems. digital imaging considered. The "manual" optimization of the adjustment parameters of digital imaging systems such as X-ray or tomography systems therefore currently involves considerable work and time, and is largely related to the dexterity of the person who uses it. performs.
Objet et résumé de l'invention Object and summary of the 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 :  The present invention notably makes it possible to remedy these drawbacks by proposing a non-destructive testing method carried out on an aeronautical part comprising:
— 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 ; A step of setting up a digital imaging system, said digital imaging system being an X-ray or digital tomography system comprising an electromagnetic radiation source and a detector able to detect the electromagnetic rays emitted by the source;
— 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 A step of acquiring at least one digital image of the aeronautical part by means of the set radiography or digital tomography system; and
— 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 ;  A non-destructive testing step performed on the aeronautical part using said at least one acquired digital image;
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 : This non-destructive testing method is remarkable in that it further comprises a prior step of estimating adjustment parameters of the digital imaging system, this estimation step comprising:
— une étape d'obtention d'un modèle numérique caractérisant le 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 ; A step of obtaining a digital model characterizing the digital imaging system, said digital model comprising a model characterizing the source and a model characterizing the detector; A step of obtaining a digital model characterizing the aeronautical part;
— une étape d'évaluation, pour une pluralité de valeurs distinctes 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  A step of evaluating, for a plurality of distinct values of the adjustment parameters, a contrast-to-noise ratio obtained for the digital imaging system and said aeronautical part, this evaluation step using the numerical models obtained characterizing the digital imaging system and the aeronautical part; and
— une étape de détermination automatique de valeurs dites optimales des 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.  A step of automatically determining so-called optimal values of the adjustment parameters optimizing the contrast-to-noise ratio, these optimum values being applied to the digital imaging system during the adjustment step.
Corrélativement, l'invention vise également un système de contrôle non destructif d'une pièce aéronautique comprenant :  Correlatively, the invention also relates to a non-destructive control system of an aeronautical part comprising:
— un système 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 détecteur apte à détecter les rayons électromagnétiques émis par la source ; et A digital imaging system, this digital imaging system being a radiography or digital tomography system comprising a source of electromagnetic rays and a detector able to detect the electromagnetic rays emitted by the source; and
— 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 aéronautique acquise par le système d'imagerie numérique ; A non-destructive control device for the aeronautical part configured to use at least one digital image of the aeronautical part acquired by the digital imaging system;
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, lesdltes valeurs optimales de paramètres de réglage étant destinés à être appliquées audit 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 :  This non-destructive control system is remarkable in that it furthermore comprises a device for estimating so-called optimal values of adjustment parameters of the digital imaging system, the high optimum values of adjustment parameters being intended to be applied to said system. digital imaging device before acquiring said at least one digital image of the aeronautical part during an adjustment of the digital imaging system, said estimating device comprising:
— 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 numérique comprenant un modèle caractérisant la source et un modèle caractérisant le détecteur ;  A first obtaining module configured to obtain a digital model characterizing the digital imaging system, this digital model comprising a model characterizing the source and a model characterizing the detector;
— un deuxième module d'obtention, configuré pour obtenir un modèle numérique caractérisant la pièce aéronautique ;  A second obtaining module configured to obtain a digital model characterizing the aeronautical part;
— un module d'évaluation configuré pour évaluer, pour une 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 An evaluation module configured to evaluate, for a plurality of distinct values of the adjustment parameters, a contrast-to-noise ratio obtained for the digital imaging system and said piece aeronautical, this evaluation module using the numerical models obtained characterizing the digital imaging system and the aeronautical part; and
— un module de détermination automatique, configuré pour déterminer automatiquement des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit.  An automatic determination module, configured to automatically determine optimum values of the adjustment parameters optimizing the contrast-to-noise ratio.
La source de rayons électromagnétiques est par exemple une source de rayons X.  The source of electromagnetic radiation is, for example, an X-ray source.
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 pièce aéronautique. Ces paramètres comprennent par exemple de façon classique, au moins un paramètre parmi :  Thus, the invention proposes a method and a system for automatically determining the adjustment parameters of a radiographic or tomographic type digital imaging system intended to be used in a non-destructive inspection process of an aeronautical part. . These parameters comprise, for example conventionally, at least one parameter among:
— une tension d'opération de la source du système d'imagerie numérique ; An operating voltage of the source of the digital imaging system;
— 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 ;  A time of exposure of the aeronautical part to a ray beam emitted by the source of the digital imaging system;
— une valeur de courant appliquée à la source du système d'imagerie numérique ; et  A current value applied to the source of the digital imaging system; and
— une épaisseur de filtre appliqué dans le système d'imagerie numérique à la source (aussi désigné par filtre externe à la source).  A filter thickness applied in the digital imaging system at the source (also referred to as an external filter at the source).
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 automatique 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 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 œuvre l'invention.  The method according to the invention is based for this purpose on various numerical models of the digital imaging system and the aeronautical part which allow the automatic determination of optimal adjustment parameters: thanks to the invention, an effective adjustment of these systems while removing any subjectivity that may affect this setting. These numerical models may for example be predetermined and stored in databases, in the form of a set of values and / or curves. The obtaining of these models then includes access to these databases and the extraction of the relevant modeled physical parameters from these bases to implement the 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é, préférentiellement de manière analytique, à partir des différentes modélisations réalisées conformément à l'invention : modélisation du système d'imagerie numérique, et plus particulièrement 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 souhaite acquérir des images au moyen du système d'imagerie numérique. The criterion chosen by the invention for optimizing the aforementioned adjustment parameters is the contrast-to-noise ratio which makes it possible to quantify the quality of the images provided by the system. digital imaging considered. This ratio is determined, preferably analytically, from the different modelizations made according to the invention: modeling of the digital imaging system, and more particularly of its main components, namely its source of electromagnetic radiation and the detector it uses to detect the electromagnetic rays emitted by the source, and modeling of the aeronautical part on which is carried out the non-destructive control and of which one wishes to acquire images by means of the digital imaging system.
La démarche ainsi adoptée par l'invention permet d'offrir une solution de réglage générique et automatique pouvant s'appliquer à de nombreux systèmes d'imagerie numérique conventionnels de radiographie et de tomographie. En recourant à un modèle numérique du système d'imagerie numérique utilisé lors du contrôle non destructif tenant compte de ses paramètres physiques, l'invention permet avantageusement de prendre en compte les spécificités et l'unicité de ce système.  The approach thus adopted by the invention provides a generic and automatic adjustment solution that can be applied to many conventional digital imaging systems for radiography and tomography. By using a digital model of the digital imaging system used during the nondestructive testing taking into account its physical parameters, the invention advantageously allows to take into account the specificities and uniqueness of this system.
En outre, le réglage proposé par l'invention tient compte de la 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.  In addition, the setting proposed by the invention takes into account the aeronautical part on which a non-destructive test is carried out. It is therefore adapted to this piece and allows to obtain an optimal image quality for this piece.
L'invention offre donc une solution rapide et efficace permettant d'automatiser le paramétrage de systèmes d'imagerie numérique utilisés dans le cadre de contrôle non destructif. L'invention ne requiert aucun traitement d'images numériques en deux ou trois dimensions issues du système d'imagerie que l'on chercher à régler n'est requis, qui peut s'avérer lourd et compliqué et s'oppose en tout état de cause à la détermination automatique de paramètres de réglage optimaux du système d'imagerie numérique. De plus, la précision de la solution proposée par l'invention 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.  The invention thus provides a quick and efficient solution for automating the parameterization of digital imaging systems used in the context of non-destructive testing. The invention does not require any digital image processing in two or three dimensions from the imaging system that is to be tuned is required, which can be cumbersome and complicated and opposes in any state of because of the automatic determination of optimal adjustment parameters of the digital imaging system. In addition, the accuracy of the solution proposed by the invention makes it easily applicable in an industrial environment to any digital imaging system of radiographic or tomographic type, for any part from the aeronautical industry.
En outre, elle s'applique à de nombreux défauts susceptibles 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 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. In addition, it applies to many defects that may be present in the room, such as in particular a sub-thickness or an extra thickness of material, porosity, inclusion, etc. It should also be noted that, thanks to the invention, by optimizing the parameters of adjustment of the digital imaging system, it ensures an optimal use of the latter and its components: the invention makes it possible in particular to reduce the use of the source of electromagnetic radiation and to limit the rate of aging of the detector used by the digital imaging system.
Comme mentionné précédemment, l'invention s'appuie sur la modélisation des principaux composants du système d'imagerie numérique 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 éléments sur la qualité de l'image numérique fournie par le système dlmagerie numérique.  As mentioned above, the invention is based on the modeling of the main components of the digital imaging system and the aeronautical part on which it is desired to carry out non-destructive testing. This modeling makes it possible to quantify the impact of these elements on the quality of the digital image provided by the digital imaging system.
Préférentiellement, on considère pour les principaux composants du système d'imagerie numérique, une modélisation physique caractérisant la réponse de ces composants en énergie.  For the main components of the digital imaging system, it is preferable to consider a physical model characterizing the response of these components to energy.
Par exemple, le modèle caractérisant la source et le modèle caractérisant le détecteur fournissent des réponses en énergie respectivement de la source et du détecteur en fonction d'une longueur d'onde. Ceci permet avantageusement de prendre en considération le caractère polychromatique de la source.  For example, the model characterizing the source and the model characterizing the detector provide energy responses respectively of the source and the detector as a function of a wavelength. This advantageously makes it possible to take into consideration the polychromatic nature of the source.
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.  Thus, in a particular embodiment, the model characterizing the source provides emission spectra of the source for different operation voltage values of the source.
Chaque spectre d'émission fournit l'intensité des rayons émis 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 comme mentionné précédemment. 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), de données empiriques et/ou de d'expressions analytiques.  Each emission spectrum provides the intensity of the rays emitted by the source (i.e. the amount of photons generated per unit time or emission) as a function of the wavelength. Such modeling advantageously allows to take into account the polychromatic nature of the source as mentioned above. It takes into account the geometry of the source, and can be determined by means of simulations (eg Monte-Carlo type), empirical data and / or analytical expressions.
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 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. In a particular embodiment, the model characterizing the source further provides estimates of at least two half-attenuation layers of the source for different operating voltage values of the source, the method further comprising a step of validation of the model characterizing the source by comparing the estimates said at least two half-attenuation layers of the source to experimental results.
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 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é.  This validation step makes it possible to ensure that the numerical model used to characterize the source is well representative of the physical reality of the source and therefore that it is well adapted. This ensures that the adjustment parameters extracted by the invention are optimal for the imaging system considered.
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 rayons émis par la source.  In a particular embodiment, the model characterizing the detector provides a spectral response of the detector to a ray beam emitted by the 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 caractère polychromatique de la source.  Such a response characterizes the way in which the detector weights the energies of the different photons arriving at the source, as a function of the wavelength. As for the source, the model considered for the detector also allows to take into consideration the polychromatic nature of the source.
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 :  In a particular embodiment, the method further comprises a validation step of the model characterizing the detector comprising:
— une estimation d'au moins un niveau de gris théorique, de 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 ;  An estimate of at least one theoretical gray level, preferably at least 3 levels, seen by the detector for at least one given operating voltage of the source, by using the model characterizing the source and the model characterizing the detector ;
— 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 une tension d'opération donnée de la source ;  Obtaining a measurement of at least one gray level, preferably at least 3 levels, seen by the detector for said at least one given operating voltage of the source;
— une comparaison dudit au moins un niveau de gris théorique estimé avec ledit au moins un niveau de gris mesuré.  A comparison of said at least one estimated theoretical gray level with said at least one measured gray level.
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 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.  In this way, it is ensured that the numerical model used to characterize the detector and extract the optimal adjustment parameters of the digital imaging system is well representative of reality. This makes it possible to guarantee the relevance of the optimal parameters determined in accordance with the invention.
Comme mentionné précédemment, l'invention s'appuie sur une modélisation numérique et physique des principaux composants du 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 d'énergie d'un faisceau de rayons émis par la source. As mentioned above, the invention relies on a numerical and physical modeling of the main components of the digital system but also of the aeronautical part on which it is desired to perform a non-destructive inspection. In a particular embodiment, the numerical model characterizing the aeronautical part provides mass attenuation coefficients of a component material for different energy values of a ray beam emitted by the 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 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.  Indeed, each part of which it is desired to acquire a digital image by means of the digital imaging system to carry out a non-destructive inspection, has its own physical characteristics (partly related to the composition of the material which composes it) which influence the system. digital imaging and its operation. Taking into account the composition of the aeronautical part and more particularly its mass attenuation coefficient for different energy values of the rays emitted by the source makes it possible to adjust the digital imaging system in a lighted and optimal way for the part in question.
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 :  For non-elementary materials, the obtaining of the numerical model characterizing the aeronautical part comprises for example, in a particular embodiment:
— une étape d'obtention d'une composition chimique du matériau composant la pièce, ledit matériau étant composé d'au moins un élément chimique ;  A step of obtaining a chemical composition of the material composing the part, said material being composed of at least one chemical element;
— 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 For each chemical element composing the material, a step of obtaining a mass attenuation coefficient of this element for said different energy values; and
— 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.  A step of calculating the mass attenuation coefficients of the material for said different energy values from the mass attenuation coefficients of the elements composing the material for these different energies and of the proportions of these elements in the material.
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).  The mass attenuation coefficient of each element of the material can be obtained for example by consulting a National Institute of Standards and Technology (NIST) database.
Cette procédure permet de déterminer de façon automatique le coefficient massique d'atténuation massique d'une pièce aéronautique pour chaque énergie de rayons considérée.  This procedure makes it possible to automatically determine the mass attenuation mass coefficient of an aeronautical part for each ray energy considered.
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 cherche à optimiser, puis les valeurs de ces paramètres conduisant au CNR optimal sont déterminées à partir des rapports ainsi évalués. The criterion used by the invention to extract the optimal adjustment parameters of the digital system is the Contrast-to-Noise Ratio (CNR). According to the invention, this The ratio is evaluated for different values of the parameters that are to be optimized, and then the values of these parameters leading to the optimal CNR are determined from the reports thus evaluated.
Dans un mode particulier de réalisation, le rapport contraste- sur-bruit est évalué, pour ladite pluralité de valeurs distinctes des paramètres de réglage, à partir de signaux caractérisés par des paramètres physiques du système d'imagerie numérique, autrement dit, du système réel considéré.  In a particular embodiment, the contrast-on-noise ratio is evaluated, for said plurality of distinct values of the adjustment parameters, from signals characterized by physical parameters of the digital imaging system, in other words, of the real system. considered.
Le rapport contraste-sur-bruit peut notamment être évalué pour différentes valeurs des paramètres que l'on cherche à optimiser, à partir de valeurs de signaux vus par le détecteur du système d'imagerie numérique lorsque des épaisseurs déterminées de la pièce aéronautique sont traversées par le faisceau de rayons électromagnétiques émis par la source, ces valeurs de signaux étant obtenues de manière analytique (i.e. au moyen d'une formule analytique) à partir des modèles numériques caractérisant le système d'imagerie numérique et la pièce aéronautique.  The contrast-to-noise ratio can in particular be evaluated for different values of the parameters that are to be optimized, based on signal values seen by the detector of the digital imaging system when determined thicknesses of the aeronautical part are crossed. by the beam of electromagnetic rays emitted by the source, these signal values being obtained analytically (ie by means of an analytical formula) from the numerical models characterizing the digital imaging system and the aeronautical part.
Cette évaluation analytique du rapport contraste-sur-bruit tenant compte des paramètres physiques du système numérique facilite la détermination automatique des paramètres de réglage optimaux du système d'imagerie numérique. L'invention est en effet facile à mettre en œuvre car elle s'appuie sur de simples calculs numériques qui plus est scalaires : le calcul du rapport contraste-sur-bruit ne nécessite en effet que la connaissance des valeurs des signaux vus par le détecteur du système d'imagerie numérique lorsque des épaisseurs déterminées de la pièce aéronautique sont traversées par le faisceau de rayons électromagnétiques émis par la source (ces valeurs reflétant par exemple un niveau de gris) ; en outre, pour des paramètres de réglage donnés du système d'imagerie numérique, ces valeurs ne varient elles-mêmes qu'en fonction d'une seule grandeur, à savoir l'épaisseur considérée de la pièce traversée par le faisceau de rayons émis par la source.  This analytical evaluation of the contrast-to-noise ratio taking into account the physical parameters of the digital system facilitates the automatic determination of the optimal adjustment parameters of the digital imaging system. The invention is indeed easy to implement because it relies on simple numerical calculations which are more scalar: the calculation of the contrast-on-noise ratio requires only the knowledge of the values of the signals seen by the detector. of the digital imaging system when determined thicknesses of the aeronautical part are traversed by the beam of electromagnetic rays emitted by the source (these values reflecting for example a gray level); in addition, for given adjustment parameters of the digital imaging system, these values themselves vary only as a function of a single quantity, namely the considered thickness of the part traversed by the ray beam emitted by source.
L'invention ne requiert donc pas de traitement lourd et compliqué d'images numériques en deux ou trois dimensions issues du système d'imagerie que l'on chercher à régler.  The invention therefore does not require heavy and complicated processing of digital images in two or three dimensions from the imaging system that we seek to adjust.
En outre, grâce aux étapes de validation des modèles numériques mises en œuvre le cas échéant, on garantit que la formule analytique utilisée pour déterminer automatiquement les paramètres de réglage du système d'imagerie numérique reflète de façon précise les caractéristiques techniques du système d'imagerie numérique. In addition, thanks to the validation steps of the numerical models implemented where appropriate, it is ensured that the analytical formula used to automatically determine the parameters of Setting the digital imaging system accurately reflects the technical characteristics of the digital imaging system.
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é pour différentes valeurs des paramètres de réglage selon l'expression :  In a particular embodiment, during the evaluation step, the noise-contrast ratio noted CNR obtained for the digital imaging system and said aeronautical part is evaluated for different values of the adjustment parameters according to the expression:
ou : or :
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 - ε, de la pièce aéronautique est traversée par le faisceau de rayons électromagnétiques (ex. rayons X) émis par la source, ε 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 ar : where S det (L), S det L, e), denotes a signal value seen by the detector of the digital imaging system when a thickness L, respectively a thickness L - ε, of the aeronautical part is traversed by the beam of electromagnetic rays (eg X-rays) emitted by the source, ε denoting a thickness of a desired indicator during non-destructive testing in the aeronautical part, said signal values being defined by:
et and
Yj l0{Ej, V)AE{Ej, L, e)R{Ej ') Y j l 0 {E j , V) A E {E j , L, e) R {E j ' )
j  j
où : or :
— /<>(£/, V) désigne l'intensité des photons présents dans un faisceau de rayons électromagnétiques émis par la source ayant pour énergie Ej 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 ; - / <> (£ /, V) denotes the intensity of the photons present in a beam of electromagnetic rays emitted by the source having for energy Ej corresponding to a wavelength λ ; · Said intensity depending on the operating voltage V of the source, the exposure time of the aeronautical part to the ray beam emitted by the source, and the current value applied to the source;
— A{Ejt L) désigne un facteur d'atténuation subie par le faisceau de rayons électromagnétiques émis par la source ayant pour énergie Ej , 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 ; - A {E jt L) denotes an attenuation factor experienced by the beam of electromagnetic rays emitted by the source whose energy Ej, this attenuation factor dependent on a mass attenuation coefficient of the material composing the aeronautical part and the thickness L;
— AE(Ej, L, E) désigne un facteur d'atténuation subie par le faisceau de rayons électromagnétiques émis par la source ayant pour énergie £) , 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 ε de l'indicateur ; et - A E (Ej, L, E) designates an attenuation factor experienced by the beam of electromagnetic rays emitted by the source having energy)), this attenuation factor depending on a mass attenuation coefficient of the component material the aeronautical part, of the thickness L, and the thickness ε of the indicator; and
— R(Ej) désigne la réponse spectrale du détecteur en fonction de l'énergie E}. - R (Ej) denotes the spectral response of the detector as a function of energy E } .
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 de réglage optimaux pour le système d'imagerie numérique adaptés aux principaux composants du système d'imagerie numérique et à la pièce aéronautique faisant l'objet du contrôle non destructif.  The calculation proposed by the invention gives an accurate and reliable evaluation of the quality of the images provided by the digital system for the aeronautical part under consideration. This results in the extraction of optimal adjustment parameters for the digital imaging system adapted to the main components of the digital imaging system and to the aeronautical part under non-destructive testing.
Il convient de noter que le modèle proposé par l'invention pour calculer le rapport contraste-sur-bruit s'applique, moyennant une définition appropriée de l'épaisseur ε à 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é 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 facilite sa détection.  It should be noted that the model proposed by the invention for calculating the contrast-to-noise ratio applies, with an appropriate definition of the thickness ε to different types of indicators sought in the aeronautical part during the destructive inspection. By indicators, we mean structures that are searched in the room and which are representative of the presence of a defect in the part such as for example an extra thickness (eg internal cavity extra thickness in ceramic cores for the elaboration high pressure turbine blades), under-thickness, inclusion, porosity, etc. By applying the aforementioned model, the adjustment parameters determined for the digital imaging system are also adapted to the indicator sought during non-destructive testing in the aeronautical part, and facilitates its detection.
A titre illustratif, lorsque le défaut recherché est une sous- épaisseur ou une porosité, l'épaisseur ε 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 par un manque de matière au niveau de la pièce. Lorsque le défaut recherché est une surépaisseur, l'épaisseur ε 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 ε traduit l'épaisseur de matériau de l'inclusion. As an illustration, when the desired defect is a sub-thickness or a porosity, the thickness ε is chosen to be "negative" so as to translate a variation of negative thickness with respect to the thickness L of the aeronautical part, a sub-thickness -thickness or porosity resulting in a lack of material at the room level. When the defect sought is an excess thickness, the thickness ε is on the contrary chosen "positive" to reflect the presence of a surplus of material. For a inclusion type defect, the thickness ε reflects the material thickness of the inclusion.
Dans un mode particulier de réalisation :  In a particular embodiment:
— le facteur d'atténuation A(Ej( L) est défini par : The attenuation factor A (E j ( L) is defined by:
A(Ej, L) = Afint(Ej, Lfint)Afext(Ej, Lfext )Amat(Ej, L) ; etA (E j, L) = A fint (Ej, L fint) A fext (Ej, Lfext) A mat (Ej, L); and
— le facteur d'atténuation A(Ej, L, ε) est défini par : The attenuation factor A (E j , L, ε) is defined by:
A(Ej, L, ε) = Afint(Ej, Lfint)Afext(Ej, Lfext)Amat(Ej, L - e)Aind(Ef, £) où Annt(Ej, Lfmt), AfextfEj. Liéx Amat(Ej, L) et Amat(Ej, L - e) désignent respectivement l'atténuation subie par le faisceau de rayons électromagnétiques (ex. rayons X) ayant pour énergie Ej due à un filtrage inhérent à la source, l'atténuation subie par le faisceau de rayons électromagnétiques ayant pour énergie Ej due à un filtrage externe à la source, l'atténuation subie par le faisceau de rayons électromagnétiques ayant pour énergie Ej due au matériau composant la pièce aéronautique et l'atténuation subie par le faisceau de rayons électromagnétiques ayant pour énergie Ej en présence de l'indicateur recherché d'épaisseur ε. A (E j , L, ε) = A fint (Ej, L fint ) A fext (Ej, L fext ) A mat (Ej, L - e) A ind (E f , £) where Annt (Ej, Lfmt) , AfextfEj. Liéx A mat (Ej, L) and A mat (Ej, L - e) respectively denote the attenuation undergone by the beam of electromagnetic rays (eg X-rays) whose energy E j due to filtering inherent to the source, the attenuation undergone by the electromagnetic ray beam whose energy E j is due to external filtering at the source, the attenuation undergone by the electromagnetic ray beam having energy E j due to the material composing the aeronautical part and the attenuation undergone by the electromagnetic ray beam having energy E j in the presence of the desired thickness indicator ε.
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 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 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é).  The attenuation factor defined in this embodiment makes it possible to take into consideration in a very complete manner the different attenuations to which the ray beam emitted by the source is subjected when using the digital imaging system. These attenuations are modeled here by filtering operations carried out on the beam of rays emitted by the source: filtering inherent to the source (and the material composing the source), filtering external to the source, but also filtering strictly related to the material of the aeronautical part which one wishes to control including the filtering carried out by the elements of the part but also by the indicator which one seeks in the room (this indicator being able to be made up of emptiness for example in the case of a porosity).
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 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 d'accélérer la détermination des paramètres de réglage optimaux du système d'imagerie numérique. In a particular embodiment, the estimation step comprises the implementation of a plurality of iterations scanning all or part of the plurality of distinct values of the adjustment parameters, and during which the evaluation of the ratio contrast-on-noise is achieved by accessing numerical model values characterizing the source, the detector and the aeronautical part previously stored in databases. This embodiment, based on prior storage of the digital models in databases and on nested loops scanning different values of the adjustment parameters, makes it possible to accelerate the determination of the optimal adjustment parameters of the digital imaging system.
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 programmes d'ordinateurs.  In a particular embodiment, certain steps of the non-destructive inspection method, and in particular the step of estimating the adjustment parameters, are determined by computer program instructions.
En conséquence, l'invention vise aussi un programme d'ordinateur sur un support d'informations, ce programme étant susceptible d'être mis en œuvre dans un dispositif d'estimation de paramètres de réglage d'un système d'imagerie numérique ou plus généralement dans un ordinateur, ce programme comportant des 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.  Accordingly, the invention also relates to a computer program on an information carrier, this program being capable of being implemented in a device for estimating adjustment parameters of a digital imaging system or more generally in a computer, this program comprising instructions adapted to the implementation of the step of estimating a non-destructive testing method as described above.
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.  This program can use any programming language, and be in the form of source code, object code, or intermediate code between source code and object code, such as in a partially compiled form, or in any other form desirable shape.
L'invention vise aussi un support d'informations lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus.  The invention also relates to a computer-readable information medium, comprising instructions of a computer program as mentioned above.
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 d'enregistrement magnétique, par exemple une disquette (floppy dise) ou un disque dur.  The information carrier may be any entity or device capable of storing the program. For example, the medium may comprise storage means, such as a ROM, for example a CD ROM or a microelectronic circuit ROM, or a magnetic recording medium, for example a floppy disk or a disk. hard.
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. 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. On the other hand, the information medium may be a transmissible medium such as an electrical or optical signal, which may be conveyed via an electrical or optical cable, by radio or by other means. The program according to the invention can be downloaded in particular on an Internet type network. Alternatively, the information carrier may be an integrated circuit in which the program is incorporated, the circuit being adapted to execute or to be used in the execution of the method in 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. Brève description des dessins  It can also be envisaged, in other embodiments, that the method and the non-destructive testing system according to the invention present in combination all or part of the aforementioned characteristics. Brief description of the drawings
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 :  Other features and advantages of the present invention will emerge from the description given below, with reference to the accompanying drawings which illustrate an embodiment having no limiting character. In the figures:
— 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 numérique ; FIG. 1 represents, in its environment, a non-destructive control system according to the invention comprising a digital imaging system, a non-destructive inspection device and a device for estimating adjustment parameters of the imaging system. digital;
— 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 l ; FIG. 2 represents, in schematic form, the hardware architecture of the adjustment parameter estimation device of FIG. 1;
— la figure 3 représente, sous forme d'ordinogramme, les principales é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 œuvre par le système de contrôle non destructif de la figure 1 ;  FIG. 3 represents, in the form of a flow chart, the main steps of a non-destructive testing method according to the invention, in a particular embodiment in which it is implemented by the non-destructive control system of the figure 1 ;
— les figures 4A, 4B et 4C illustrent respectivement les principales étapes mises en œuvre par le dispositif d'estimation de paramètres de 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 FIGS. 4A, 4B and 4C respectively illustrate the main steps implemented by the adjustment parameter estimation device of FIG. 1 to obtain a digital model of the source and a digital model of the detector of the digital imaging system. of Figure 1, and a digital model of an aeronautical part on which operates the non-destructive control system of Figure 1; and
— la figure 5 représente, sous forme d'ordinogramme, les principales étapes mises en œuvre lors de l'étape d'estimation des paramètres de réglage optimaux mise en œuvre dans le procédé de contrôle non destructif illustré à la figure 3. Description détaillée de l'invention FIG. 5 represents, in the form of a flow chart, the main steps implemented during the step of estimating the optimal adjustment parameters implemented in the non-destructive inspection method illustrated in FIG. 3. Detailed description of the invention
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 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.  FIG. 1 represents, in its environment, a non-destructive testing system 1 according to the invention, in a particular embodiment. In the example envisaged here, the system 1 makes it possible to carry out non-destructive testing of aeronautical parts such as, for example, a fan blade 2 of a turbine. However, no limitation is attached to the nature of the part on which non-destructive testing is applied. It may be more generally any type of parts, preferably aeronautical, such as for example a part equipping an aircraft engine, rocket, etc.
Conformément à l'invention, le système 1 de contrôle non destructif comprend :  According to the invention, the non-destructive control system 1 comprises:
— 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 A digital imaging system 3, able to acquire and supply one or more digital images IM of the aeronautical part 2; and
— 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 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 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. A non-destructive testing device 4 for the aeronautical part 2 configured to use the digital IM image (s) of the aeronautical part acquired and supplied by the digital imaging system 3 and form an integrity diagnosis DIAG of part 2. This diagnosis concerns in particular one or more IND indicators or parameters searched in the part. No limitation is attached to the nature of these indicators; these may be, for example indicators representative of thickness in the room, of under-thickness, or the presence of porosities or inclusions, etc. Such a non-destructive testing device is known per se and is not described in more detail here.
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 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.  The digital imaging system 3 is here a radiography or digital tomography system. It comprises, in a known manner, a source 3A of electromagnetic rays (for example a source emitting beams of X-rays) and a detector 3B able to detect the electromagnetic rays emitted by the source 3A.
La source 3A applique ici le principe du bombardement d'électrons sur une cible métallique dans un tube à rayons X : les é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 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. The source 3A applies here the principle of the bombardment of electrons on a metal target in an X-ray tube: the electrons are extracted from a metallic filament (the cathode) and then accelerated by a large electrical voltage applied to the source 3A (voltage operating in the sense of the invention) in a vacuum tube. This electron beam emitted by the source 3A is focused so as to bombard a metal target (the anode), for example tungsten. Interactions between electrons and matter when they hit the target produce 'X' ray photons. Such operation is known per se and is not described in more detail here.
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 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 photomultiplicateur. Le signal électrique véhicule des niveaux de gris associés à chaque pixel du détecteur et reflétant l'absorption des photons au niveau du pixel considéré. Un tel fonctionnement est connu en soi et n'est pas décrit plus en détail ici.  The detector 3B of the digital imaging system 3 operates indirectly: a scintillator absorbs the X-ray photons produced by the source 3A, and by fluorescence mechanisms, emits visible light. The visible photons emitted are directed on a photodiode array, each element of the photodiode array being associated with a pixel. The visible photons are converted into electrical charge at the level of the photodiode array and then into an electrical signal by a photomultiplier. The electrical signal conveys gray levels associated with each pixel of the detector and reflecting the absorption of photons at the pixel considered. Such operation is known per se and is not described in more detail here.
Pour fournir des images numériques IM de la pièce aéronautique 2 permettant la mise en œuvre d'un contrôle non destructif 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 :  To provide IM digital images of the aeronautical part 2 enabling the implementation of effective non-destructive testing by the non-destructive testing device 4, it is conventional to make certain adjustments to the digital imaging system 3. These adjustments relate in particular to different parameters of the digital imaging system 3, comprising in the embodiment described here:
— la tension V d'opération de la source 3A ;  The operating voltage V of the source 3A;
— le temps d'exposition de la pièce aéronautique 2 au faisceau de rayons émis par la source 3A ;  The exposure time of the aeronautical part 2 to the beam of rays emitted by the source 3A;
— la valeur de courant i appliquée à la source 3A ; et  The current value i applied to the source 3A; and
— une épaisseur de filtre appliqué dans le système d'imagerie numérique à la source 3A (aussi désigné par filtre externe à la source).  A filter thickness applied in the digital imaging system to the source 3A (also referred to as an external filter at the 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.  This list is however not exhaustive and other adjustment parameters of the digital imaging system 3 can be envisaged; likewise, only a part of the aforementioned parameters can be envisaged.
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 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. According to the invention, the non-destructive control system 1 advantageously comprises a device 5 for automatically estimating the aforementioned adjustment parameters. This estimation device 5 is able to estimate optimal values of these parameters. to be applied to the digital imaging system 3 during an adjustment thereof made before the acquisition of the digital images IM of the aeronautical part 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.  In the embodiment described here, the estimating device 5 has the hardware architecture of a computer as diagrammatically shown in FIG. 2.
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 dlmagerie numérique d'interagir avec le dispositif 5 d'estimation. Ces moyens d'entrée/sortie 10 comprennent par exemple un écran, un clavier, 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 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.  It comprises, in particular a processor 6, a random access memory 7, a read only memory 8, a nonvolatile flash memory 9 as well as input / output means 10 enabling an operator of the digital imaging system 3 to interact with the device 5 estimation. These input / output means 10 comprise, for example, a screen, a keyboard, a man-machine interface enabling the operator of the system 3 in particular to enter data in the estimation device 5 and to obtain the adjustment parameters. optimal digital imaging system 3 determined by it in accordance with the invention.
La mémoire morte 8 du dispositif d'estimation 5 constitue un support d'enregistrement conforme à l'invention, lisible par le processeur 6 et sur lequel est enregistré un programme d'ordinateur conforme à l'invention.  The read-only memory 8 of the estimation device 5 constitutes a recording medium in accordance with the invention, readable by the processor 6 and on which a computer program according to the invention is recorded.
Le programme d'ordinateur PROG définit des modules fonctionnels et logiciels ici, configurés pour mettre en œuvre une estimation des paramètres de réglage optimaux du système 3 d'imagerie 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 :  The PROG computer program defines functional and software modules herein, configured to implement an estimate of the optimal setting parameters of the digital imaging system 3 according to the invention. These functional modules rely on and / or control the hardware elements 6-10 of the estimation device 5 mentioned above. They include in particular here, as illustrated in FIG.
— un premier module d'obtention 5A configuré pour obtenir un modèle 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 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 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 ; A first obtaining module 5A configured to obtain a numerical model MOD3 characterizing the digital imaging system 3, this numerical model comprising a MOD3A model characterizing the source 3A and a MOD3B model characterizing the detector 3B. The first obtaining module 5A is configured here to extract these MOD3A and MOD3B digital models from the DB3A and DB3B databases of the estimation device, in which they were previously stored after their determination. By database is meant in the broad sense any set of data for storing information; A second obtaining module 5B, configured to obtain a numerical model MOD2 characterizing the aeronautical part 2; the second obtaining module 5B is configured here to extract the digital model MOD2 from a database DB2 of the estimating device in which it was previously stored after its determination;
— un module d'évaluation 5C configuré pour évaluer, pour une pluralité de valeurs distinctes des paramètres de réglage, un rapport contraste sur 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 5 A et 5B et caractérisant respectivement le système d'imagerie numérique 3 et la pièce aéronautique 2 ; et  A evaluation module 5C configured to evaluate, for a plurality of distinct values of the adjustment parameters, a contrast-to-noise ratio obtained for the digital imaging system 3 and the aeronautical part 2, this evaluation module 5C using the numerical models MOD3 (including the models MOD3A and MOD3B) and MOD2 obtained by the modules 5A and 5B and respectively characterizing the digital imaging system 3 and the aeronautical part 2; and
— un module de détermination automatique 5D, configuré pour déterminer automatiquement des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit.  An automatic determination module 5D, configured to automatically determine optimum values of the adjustment parameters optimizing the contrast-to-noise ratio.
Les fonctions des différents modules sont décrites plus en détail en référence à l'étape d'estimation du procédé de contrôle non destructif selon l'invention mise en œuvre par le dispositif 5 d'estimation du système 1 de contrôle non destructif de la figure 1.  The functions of the various modules are described in more detail with reference to the estimation step of the non-destructive testing method according to the invention implemented by the estimation device of the non-destructive testing system 1 of FIG. .
On note que deux modules d'obtention distincts 5A et 5B sont décrits ici, toutefois il est possible que les fonctions mises en œuvre par ces deux modules soient mises en œuvre dans un autre mode de réalisation par un module unique.  Note that two separate obtaining modules 5A and 5B are described here, however it is possible that the functions implemented by these two modules are implemented in another embodiment by a single module.
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 ces étapes telles qu'elles sont mises en œuvre par le système 1 de contrôle non destructif de la figure 1.  To better understand the advantages provided by the invention, we will now briefly describe the main steps of the non-destructive testing method according to the invention. FIG. 3 illustrates these steps as they are implemented by the non-destructive testing system 1 of FIG. 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 du système 3 d'imagerie numérique (système radiographique ou tomographique).  According to the invention, it is assumed here that it is desired to carry out non-destructive testing of the aeronautical part 2 from one or more digital images of this part acquired by means of the digital imaging system 3 (radiographic system or tomographic).
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é précédemment : For this purpose, a preliminary adjustment of various parameters of the digital imaging system 3 is performed (step D10). In the mode embodiment described here, these parameters include, as indicated above:
— la tension V d'opération de la source 3A ;  The operating voltage V of the source 3A;
— le temps d'exposition de la pièce aéronautique 2 au faisceau de rayons émis par la source 3A ;  The exposure time of the aeronautical part 2 to the beam of rays emitted by the source 3A;
— la valeur de courant i appliquée à la source 3A ; et  The current value i applied to the source 3A; and
— l'épaisseur de filtre appliqué dans le système d'imagerie numérique à la source.  - the filter thickness applied in the digital imaging system at the source.
Conformément à l'invention, le réglage de ces paramètres est effectué en deux temps :  According to the invention, the adjustment of these parameters is carried out in two stages:
— dans un premier temps, des valeurs optimales de ces paramètres sont estimées automatiquement par le dispositif 5 d'estimation (étape D10- i) ;  Firstly, optimal values of these parameters are automatically estimated by the estimation device (step D10-i);
— 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 And then in a second step, the optimal values thus estimated are applied to the digital imaging system 3 to parameterize it.
(é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 3 au moyen des valeurs optimales fournies par le dispositif 5 d'estimation. (step D10-2). This setting of the digital imaging system 3 can be performed for example manually by an operator of the system 3 by means of the optimum values provided by the estimating device.
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.  Then one or more IM digital images of the aeronautical part 2 are acquired by the digital imaging system 3 and adjusted (step D20), and provided to the non-destructive testing device 4.
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 à 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.  Non-destructive testing of the aeronautical part 2 is then performed via the non-destructive testing device 4 from the digital images IM (step D30). Such a step can be carried out in various ways, in a manner known per se, and is not described in detail here. The images thus obtained allow non-destructive access to the part to be controlled. For example, they can be compared by known image processing methods to images acquired by the same method for a sound part in the search for anomaly (s) in the part to be controlled.
Conformément à l'invention, l'estimation réalisée à l'étape D10-1 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 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 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. In accordance with the invention, the estimate made in step D10-1 of the optimal adjustment parameters of the digital imaging system 3 for its use on the aeronautical part 2 is based on various numerical models physically characterizing the system 3. digital imaging, and in particular its main components namely the source 3A and the detector 3B, and the aeronautical part 2. These numerical models are obtained by the modules 5A and 5B respectively of the estimation device 5. They can be pre-calculated or obtained on the fly when estimating optimal tuning parameters. However, in the embodiment described here, these numerical models are precalculated and stored in databases (DB3A, DB3B and DB2) of the estimation device 5 for future use. This advantageously makes it possible to accelerate the search for optimum adjustment parameters when it is desired to use the digital imaging system 3.
Nous allons dans un premier temps décrire comment ces modèles numériques sont obtenus (i.e. construits) et mis à disposition des modules 5A et 5B du dispositif 5 d'estimation dans les bases de données DB3A, DB3B et DB2.  We will first describe how these numerical models are obtained (i.e. constructed) and made available the modules 5A and 5B of the estimation device 5 in the DB3A, DB3B and DB2 databases.
La figure 4A représente les principales étapes mises en uvre pour obtenir le modèle OD3A caractérisant la source 3A du système d'imagerie numérique 3. On note que ces étapes peuvent être mises en œuvre par le module 5A du dispositif d'estimation 5 ou via tout autre module externe au dispositif d'estimation 5 et apte à stocker le modèle MOD3A obtenu dans la base de données DB3A.  FIG. 4A represents the main steps implemented to obtain the OD3A model characterizing the source 3A of the digital imaging system 3. It is noted that these steps can be implemented by the module 5A of the estimation device 5 or via any another module external to the estimation device 5 and able to store the model MOD3A obtained in the database 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 i, v2, w où N désigne un entier supérieur à 1. Dans l'exemple envisagé ici, les tensions d'opération Vlt V2, ... , VN varient de 0 à 450 kV, par pas de lkV. The model MOD3A here provides emission spectra of the source 3A for different operating voltage values of the source denoted i, v 2 , w where N denotes an integer greater than 1. In the example considered here, the voltages V Lt V 2 , ..., V N range from 0 to 450 kV, in steps of lkV.
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 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).  In the embodiment described here, the MOD3A model further provides emission spectra of the source 3A for N operating voltage values of the source and for M values of the tilt angle with respect to horizontal of the anode of the source 3A, M denoting an integer greater than 1. It is noted that the angle of inclination of the anode of the source influences the emission angle of the source 3A; which has an importance on the quality of the spectrum (Heel effect).
Par spectre d'émission, on entend l'intensité (i.e. le nombre de 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 faisceaux de rayons considérés. Autrement dit à une longueur d'onde indexé par un entier j et notée Λ, correspond une unique énergie de faisceaux Es. By emission spectrum is meant the intensity (ie the number of photons emitted per unit of time or emission) beams of rays emitted by the source 3A as a function of the wavelength of these rays. Such modeling advantageously makes it possible to take into consideration the polychromatic nature of the source 3A. We take note that the wavelength is uniquely related to the energy of the beams of rays considered. In other words, at a wavelength indexed by an integer j and denoted by Λ, there corresponds a single energy of beams E s .
Pour déterminer les spectres d'émission du modèle MOD3A correspondant aux N valeurs possibles de la tension d'opération de la 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 3 A. Ce logiciel est par exemple ici le logiciel SpekCalc, développé par V.N. Hansen et al., et décrit notamment dans le document intitulé « An efficient Monte Carlo- based algorithm for scatter correction in kev, cone-beam et », Physics in 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 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 de l'anode de la source variant de 5 à 45 degrés par incrément de 1 degré.  To determine the emission spectra of the model MOD3A corresponding to the N possible values of the operating voltage of the source and M possible values of the angle of inclination of the anode (initialization steps E10, test E30 and E35 and of incrementation E40 and E45), one uses, in the embodiment described here, a software of simulation making it possible to model the behavior of the source 3 A. This software is for example here software SpekCalc, developed by VN Hansen et al., and described in particular in the document entitled "An efficient Monte Carlo-based algorithm for scatter correction in kev, cone-beam and", Physics in Medicine and Biology, 54: 3847-3864, 2009. It uses electron interaction at the anode, as well as data from Monte-Carlo simulations and / or empirical and / or analytic expressions for estimating emission spectra. It is noted that the SpekCalc software makes it possible to model emission spectra for angles of inclination of the source anode of between 1 and 90 degrees, however the sources of electromagnetic rays and in particular of X-rays rarely have angles. tilt beyond 45 degrees. For example here, we consider an angle of inclination of the anode of the source varying from 5 to 45 degrees in increments of 1 degree.
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 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) :  The emission spectra corresponding to the N operating voltage values and the M tilt angle values obtained using the SpekCalc software (step E20) are stored in the DB3A database (step E50). In the embodiment described here, to access one of the transmission spectra stored in the DB3A database, the module 5A provides the database DB3A (or equivalent query with):
— l'angle d'inclinaison noté at de l'anode de la source 3A par rapport à l'horizontal ; The angle of inclination noted a t of the anode of the source 3A relative to the horizontal;
— l'énergie maximum Emax(½) (en EV ou kEV) correspondant à la tension vt 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 pour la tension d'opération ^ et l'angle α,, noté dans la suite de la description ^ Υί, ^. On note que chaque spectre d'émission ΐ05(Υί. αι) fournit le nombre de photons d'énergie Ei par unité de courant et par unité de temps pour différentes énergies j=l,...,K (ou de façon équivalente, pour différentes longueurs d'onde A/J=1/...I ), K désignant un entier supérieur à 1. The maximum energy Emax (½) (in EV or kEV) corresponding to the operating voltage v t of the source 3A considered for which it is desired to estimate the emission spectrum of the source. The module 5A thus obtains the emission spectrum of the source 3A for the operating voltage et and the angle α ,, noted in the following description Υ ί ^. We note that each emission spectrum ΐ 05 (Υ ί . Αι) provides the number of energy photons E i per unit of current and per time unit for different energies j = 1, ..., K (or of equivalent way, for different wavelengths A / J = 1 / ... I ), K denoting an integer greater than 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 (Al), pour chaque tension d'opération vt et pour chaque angle d'inclinaison ,. Ces estimations ont été fournies par le logiciel SpekCalc au cours de l'étape E20. Plus précisément, la base de données D3A comprend : In the embodiment described here, database D3A also includes estimates of at least two half-attenuation layers of source 3A for two distinct materials, namely copper (Cu) and aluminum (Al). for each operating voltage v t and for each inclination angle, These estimates were provided by the SpekCalc software during step E20. Specifically, the D3A database includes:
— les couches de demi-atténuation HVLlCu^. a,) et HVLlAI(Vif af) désignant l'épaisseur de matière nécessaire (cuivre ou aluminium ici) pour diminuer de moitié la radiation du faisceau à lm de la fenêtre de la source ; et The half-attenuation layers HVL1Cu1. a,) and HVLlAI (V if a f ) designating the thickness of material required (copper or aluminum here) to halve the radiation of the beam at 1 m from the source window; and
— les couches de demi-atténuation HVUCu^a,) et HVL2AI( , a£ désignant l'épaisseur de matière nécessaire (cuivre ou aluminium ici) pour diminuer de moitié la radiation du faisceau résultant respectivement de l'utilisation des couches de demi-atténuation HVLlCu^, α,) et HVLlAI(VÉ, a,) à lm de la fenêtre de la source. - The half-attenuation layers HVUCu ^ a,) and HVL2AI (, a £ designating the thickness of material required (copper or aluminum here) to halve the radiation of the beam resulting respectively from the use of half layers attenuation HVLlCu ^, α,) and HVLlAI (V É , a) to 1m from the source window.
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 Vt et l'angle at. The module 5A here accesses these estimates at the same time as the emission spectrum of the source 3A for the operating voltage V t and the angle at t .
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 œuvre par exemple par le module 5A ou par un autre module externe au dispositif d'estimation 5, en accédant aux 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 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. In the embodiment described here, a validation of the numerical model MOD3A, ie emission spectra stored in the database D3A (step E60), is envisaged. This validation can be implemented for example by the module 5A or by another module external to the estimation device 5, by accessing the emission spectra stored in the database D3A and the estimates of the associated half-attenuation layers. . In the remainder of the description, for the sake of simplification, the module responsible for this validation will generally be designated by validation module. This validation aims here to ensure that the spectra generated for the source 3A are well representative of the physical reality of the source 3A. This contributes to allow the estimation device 5 to obtain, thanks to the invention, optimal values of the adjustment parameters.
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 cours de l'étape E20 pour deux tensions d'opération différentes de la source 3 A, les estimations des couches de demi-atténuation HVLICu, HVL1AI, HVL2Cu et HVL2AI fournies lors de l'étape E20 pour ces tensions avec des résultats expérimentaux.  For this purpose, the validation module compares here, for at least one selection of spectra comprising at least two spectra generated during the step E20 for two different operating voltages of the source 3 A, the estimates of the half layers. HVLICu, HVL1AI, HVL2Cu and HVL2AI attenuators provided in step E20 for these voltages with experimental results.
Les valeurs expérimentales des couches de demi-atténuation HVLICu, HVL1AI, HVL2Cu et HVL2AI peuvent être générées préalablement, au moyen de couches de cuivre et d'aluminium de 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.  The experimental values of the half-attenuation layers HVLICu, HVL1AI, HVL2Cu and HVL2AI can be generated beforehand, by means of layers of copper and aluminum of different thicknesses and a dosimeter, then stored in the DB3A database to be compared with the values of the half-attenuation layers accompanying the spectra in the DB3A database.
Si les valeurs HVLICu, HVL1AI, HVL2Cu et HVL2AI 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 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.  If the values HVLICu, HVL1AI, HVL2Cu and HVL2AI corresponding to each spectrum of the relevant spectral selection of the database correspond to the experimental values, then this means that the spectra generated in step E20 and stored in the database DB3A are valid. They constitute a validated MOD3A numerical model characterizing the source 3A within the meaning of the invention. Note that in practice it is sufficient to carry out this validation on a selection of spectra corresponding to a predetermined range of values of voltages and currents.
On note également que dans le mode de réalisation décrit ici, si les valeurs HVLICu, HVL1AI, HVL2Cu et HVL2AI d'un spectre de la base de 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 HVLl et HVL2 qui correspondent le mieux aux valeurs de HVLA et HVL2 trouvées expérimentalement.  Note also that in the embodiment described here, if the HVLICu, HVL1AI, HVL2Cu and HVL2AI values of a spectrum of the DB3A database do not match the experimental values, it is stored in the DB3A database. . However, the validation module will also search the DB3A database for the spectrum having HVL1 and HVL2 values that best correspond to the values of HVLA and HVL2 found experimentally.
Nous allons maintenant décrire, en référence à la figure 4B, les principales étapes mises en œuvre 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 œuvre par le module 5A du dispositif d'estimation 5 ou via tout autre module externe au dispositif d'estimation 5 et apte à stocker le modèle MOD3B obtenu dans la base de données DB3B. We will now describe, with reference to FIG. 4B, the main steps implemented to obtain the MOD3B model characterizing the detector 3B of the digital imaging system 3. It is noted that these steps can be implemented by the module 5A of FIG. estimation device 5 or via any other module external to the estimation device 5 and able to store the model MOD3B obtained in the database 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 £}, j=l,...,K. En d'autres mots, le modèle MOD3B caractérise comment le 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 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 fonction de l'énergie des photons incidents (étape F10), en fonction de l'épaisseur et type du scintillateur du détecteur 3B.  The model MOD3B provides here a spectral response denoted R of the detector to a beam of rays emitted by the energy source £}, j = l, ..., K. In other words, the MOD3B model characterizes how the detector 3B weights the different energies of photons arriving from the source. It should be noted that this spectral response depends to a large extent on the ability of the scintillator of the detector 3B to absorb the different photon energies. Therefore, in the embodiment described here, to obtain the MOD3B model, the so-called response curves showing the amount of energy absorbed in the detector scintillator are simulated beforehand with the aid of appropriate simulation software. 3B as a function of the energy of the incident photons (step F10), as a function of the thickness and type of the scintillator of the detector 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 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=l,...,K, émet un faisceau de photons d'énergie E} avec Ej < eVmax où Vmax désigne la tension d'opération maximal de la source 3A. Le faisceau émis par chaque source est non divergent et comprend un nombre suffisant de photons (i.e. une intensité), typiquement supérieur à 106. Les K sources sont superposées au même endroit. The development of such response curves is based here on the simulation of the dose deposited by the photon beam on the scintillator, that is to say the energy absorbed by the scintillator. To this end, a number K of "elementary" sources of independent photons of the same intensity each and monochromatic are modeled. Each source indexed by j, j = l, ..., K, emits a beam of energy photons E } with Ej <eV max where V max denotes the maximum operating voltage of the source 3A. The beam emitted by each source is non-divergent and comprises a sufficient number of photons (ie an intensity), typically greater than 10 6 . The K sources are superimposed on the same place.
On note que la dose du faisceau absorbée par le scintillateur dépend directement de son épaisseur. Plusieurs épaisseurs sont donc évaluées. On obtient ainsi, à l'issue des simulations, plusieurs courbes de réponse correspondant à plusieurs épaisseurs des scintillateurs (étape F20). Ceci permet par la suite au module 5A de choisir celle qui représente le mieux le détecteur 3B. Par exemple, pour un scintillateur en iodure de césium (Csl), on peut envisager des épaisseurs allant de 350 à 750 microns par pas de 10 microns, tandis que pour un scintillateur en GdOS (Gadolinium Oxygène Sulfure), on peut envisager des épaisseurs allant de 100 à 400 microns par pas de 10 microns. Dans le mode de réalisation décrit ici, les simulations réalisées à l'étape F10 sont effectuées avec le module GateDoseSpecùvmActor\mp]émertié dans la plateforme GATE/Geant 4 (Géant 4 Application for Tomography Emission), qui permet de calculer la réponse du détecteur 3B à chaque énergie. De façon connue, la plateforme GATE/Geant 4 met en œuvre des simulations de type Monte- Carlo pour modéliser la tomographie par rayons X et la radiographie. It is noted that the dose of the beam absorbed by the scintillator depends directly on its thickness. Several thicknesses are evaluated. Thus, at the end of the simulations, several response curves corresponding to several thicknesses of the scintillators (step F20) are obtained. This then allows the module 5A to choose the one that best represents the detector 3B. For example, for a cesium iodide (CsI) scintillator, thicknesses ranging from 350 to 750 microns can be envisaged in steps of 10 microns, whereas for a GdOS (Gadolinium Oxygen Sulfide) scintillator, it is possible to envisage thicknesses ranging from from 100 to 400 microns in steps of 10 microns. In the embodiment described here, the simulations carried out in step F10 are performed with the GateDoseSpecUvmActor \ mp] module emertized in the GATE / Geant 4 platform (Giant 4 Application for Tomography Emission), which makes it possible to calculate the detector response. 3B with each energy. In a known manner, the GATE / Geant platform 4 implements Monte Carlo simulations for modeling X-ray tomography and radiography.
Les réponses spectrales du détecteur 3B R(£)), j=l,...,K obtenues lors des simulations pour les différentes épaisseurs de scintillateur sont stockées dans la base de données D3B (étape F30). Elles constituent un modèle numérique MOD3B caractérisant le détecteur 3B au sens de l'invention. Elles deviennent ainsi accessibles par le module 5A du dispositif d'estimation 5 pour déterminer les paramètres de réglage optimaux du système d'imagerie numérique, comme décrit ultérieurement.  The spectral responses of the detector 3B R ()), j = 1,..., K obtained during the simulations for the different scintillator thicknesses are stored in the database D3B (step F30). They constitute a numerical model MOD3B characterizing the detector 3B within the meaning of the invention. They thus become accessible by the module 5A of the estimation device 5 to determine the optimal adjustment parameters of the digital imaging system, as described later.
Dans le mode de réalisation décrit ici, les réponses spectrales obtenues pour le détecteur 3B sont par ailleurs validées pour s'assurer de sa conformité à la réalité physique du détecteur 3B (étape F40). Cette validation peut être mise en œuvre par le module 5A ou par tout autre module de validation, par exemple externe au dispositif d'estimation 5. Dans la suite de la description, par souci de simplification, on fera référence de façon générale à un module de validation.  In the embodiment described here, the spectral responses obtained for the detector 3B are furthermore validated to ensure that it conforms to the physical reality of the detector 3B (step F40). This validation can be implemented by the module 5A or by any other validation module, for example external to the estimation device 5. In the remainder of the description, for the sake of simplification, reference will be made generally to a module of confirmation.
Pour mettre en œuvre cette validation, les inventeurs se sont appuyés sur les constats suivants :  To implement this validation, the inventors relied on the following observations:
— l'énergie moyenne et l'énergie maximale des photons dans un spectre sont des caractéristiques qui dépendent de la tension d'opération de la source, et une variation de cette dernière se manifeste dans le niveau de gris moyen noté GL vu par le détecteur ; et  The average energy and the maximum energy of the photons in a spectrum are characteristics which depend on the operating voltage of the source, and a variation of this latter is manifested in the average gray level noted GL seen by the detector ; and
— la variation de niveau de gris moyen GL en fonction de l'énergie moyenne des faisceaux de photons émis par la source dépend de la réponse spectrale du détecteur.  The average gray level variation GL as a function of the average energy of the photon beams emitted by the source depends on the spectral response of the detector.
Pour valider la réponse spectrale obtenue à l'étape F20, les inventeurs proposent donc de s'appuyer sur la méthode proposée par P.V. To validate the spectral response obtained in step F20, the inventors therefore propose to rely on the method proposed by P.V.
GRANTON et al. dans le document intitulé « A combined dose calculation and vérification method for a small animal précision irradiator based on onboard 290 imaging », Med. Phys., vol. 39 n°7, page 4155, 2012. Cette méthode, appliquée dans le contexte de l'invention, consiste à observer de manière empirique les variations de niveaux de gris du détecteur 3B en fonction de la tension d'opération de la source 3A, et à comparer ces niveaux de gris avec une estimation théorique de ces niveaux de gris évaluée en utilisant le modèle numérique MOD3A et le modèle numérique OD3B. GRANTON et al. in the document entitled "A combined dose calculation and verification method for a small animal irradiator-based precision on onboard imaging", Med. Phys., Vol. 39, No. 7, page 4155, 2012. This method, applied in the context of the invention, consists in observing empirically the gray level variations of the detector 3B as a function of the operating voltage of the source 3A, and comparing these gray levels with a theoretical estimate of these gray levels evaluated using the numerical model MOD3A and the model Digital OD3B.
Plus particulièrement, lors de l'étape de validation F40, le module de validation obtient tout d'abord une pluralité de mesures du niveau de gris moyen noté GLfe, i=l,...I, I désignant un entier supérieur à 1, vu par le détecteur 3B pour une pluralité de tensions d'opérations V V2, ... , Vi de la source 3A variant entre une tension minimum Vmin et la tension maximale Vma . Par souci de simplification ici, on suppose que I=N, et que les tensions d'opération VX, V2, ... , V} coïncident avec les tensions d'opération V V2, ... , VN précédemment considérées pour obtenir le module MOD3A. More particularly, during the validation step F40, the validation module first obtains a plurality of measurements of the average gray level denoted GLf e , where i = 1, ... I, I denoting an integer greater than 1 seen by the detector 3B for a plurality of operating voltages VV 2 ,..., Vi of the source 3A varying between a minimum voltage V min and the maximum voltage V ma . For simplicity here, it is assumed that I = N, and the operating voltages V X, V 2, ..., V} coincide with the operating voltages VV 2, ..., V N previously considered to obtain the MOD3A module.
On note que lors de cette étape, les mesures peuvent être effectuées directement par le module de vérification, ou en variante être obtenues par ceux-ci d'une base de données dans laquelle on a stocké de telles mesures effectuées préalablement.  Note that during this step, the measurements can be performed directly by the verification module, or alternatively be obtained by them from a database in which such measurements previously made have been stored.
Le module de validation estime par ailleurs des niveaux de gris théoriques GL†eo, ί=1,...Ν, pour chacune des tensions d'opération Vx. Vz vN The validation module also estimates theoretical gray levels GL † eo , ί = 1, ... Ν, for each of the operating voltages Vx. Vz v N
Dans le mode de réalisation décrit ici, un niveau de gris théorique GL e° (E}), \=l,...N correspondant à une tension d'opération Vi est estimé par le module de vérification pour un faisceau d'énergie E} en multipliant la réponse R(£}) obtenue pour le détecteur 3B (i.e. extraite du modèle numérique MOD3B) par la valeur /0s( , Ej) du spectre d'émission losiVô de 'a source 3A associé à la tension de voltage ½ (extraite au modèle numérique MOD3A et correspondant à l'angle d'inclinaison de l'anode de la source 3A). Cette opération est répétée pour toutes les énergies de faisceaux Ejt j=l,...,K et pour toutes les épaisseurs de scintillateur pour lesquelles on dispose d'une réponse spectrale dans la base de données D3B. On note que par souci de simplification, on considère ici un unique nombre K pour toutes les tensions d'opération. Toutefois en pratique les énergies considérées varient entre une valeur minimale (par exemple IkeV) et une valeur maximale qui dépend de la tension d'opération considérée. Puis le niveau de gris théorique moyen est obtenu par le module de vérification en calculant l'expression : In the embodiment described here, a theoretical gray level GL e ° (E } ), \ = 1, ... N corresponding to an operation voltage Vi is estimated by the verification module for a beam of energy E } by multiplying the response R (£}) obtained for the detector 3B (ie extracted from the numerical model MOD3B) by the value / 0s (, Ej) of the emission spectrum losiV0 of a source 3A associated with the voltage voltage ½ (extracted from the numerical model MOD3A and corresponding to the angle of inclination of the anode of the source 3A). This operation is repeated for all the beam energies E jt j = l, ..., K and for all the scintillator thicknesses for which a spectral response is available in the database D3B. Note that for the sake of simplicity, here we consider a single number K for all operating voltages. However in practice the energies considered vary between a minimum value (for example IkeV) and a maximum value which depends on the operating voltage considered. Then the average theoretical gray level is obtained by the verification module by calculating the expression:
N NOT
CL eo = ^ ff(£-y) /0s(„ 7) pour i=l,...N. CL eo = ff (£ - y ) / 0s (" 7 ) for i = 1, ... N.
Le module de vérification compare ensuite cette valeur aux valeurs empiriques obtenues précédemment. Dans le mode de réalisation décrit ici, pour effectuer cette comparaison, le module de vérification calcule tout d'abord un coefficient de normalisation pour chaque i=l,...N : The verification module then compares this value with the empirical values previously obtained. In the embodiment described here, to perform this comparison, the verification module first calculates a normalization coefficient for each i = 1, ... N:
_ Gif puis il évalue des niveaux de gris théoriques normalisés selon :  _ Gif then he evaluates theoretical theoretical gray levels according to:
s> r theo.norm Γ i theo „ s> r theo.norm Γ i theo "
ÙLi,j = ÙLi ■ Cj Ù L i, j = ÙL i ■ C j
pour i,j=l,...N. for i, j = l, ... N.
Ensuite le module de vérification calcule l'erreur suivante :  Then the verification module calculates the following error:
N NOT
SJ = S J =
1=1  1 = 1
et sélectionne le coefficient de normalisation cj correspondant à l'erreur minimale sj. On note cLfeo,norm le niveau de gris théorique correspondant à ce coefficient de normalisation. and selects the normalization coefficient cj corresponding to the minimum error sj. Note cLf eo, norm the theoretical gray level corresponding to this coefficient of normalization.
Le module de vérification compare ensuite les valeurs de CL m et de Gi eo,norm par exemp|e en traçant la courbe Si les courbes iog(CL eo*norm) et log(CL†m) coïncident, cela veut dire que le modèle MOD3B de réponse du détecteur est validé. Le module de vérification réalise ces opérations pour toutes les épaisseurs de scintillateur pour lesquelles une réponse spectrale est stockée dans la base D3B, et sélectionne la réponse spectrale qui correspond le mieux aux valeurs de niveaux de gris mesurées. The verification module then compares the values of CL m and Gi eo, norm for example p | e by plotting the curve If the iog curves (CL eo * norm ) and log (CL † m ) coincide, this means that the MOD3B response model of the detector is validated. The verification module performs these operations for all scintillator thicknesses for which a spectral response is stored in the D3B base, and selects the spectral response that best matches the measured gray level values.
On note que les étapes de validation des modèles MOD3A et MOD3B, bien que recommandées, sont toutefois optionnelles.  It should be noted that the validation steps of MOD3A and MOD3B models, although recommended, are optional.
Enfin, comme mentionné précédemment, l'invention s'appuie également sur un troisième modèle numérique MOD2 caractérisant la pièce aéronautique. Ce modèle est stocké dans la base de données DB2 de sorte à être utilisé par le deuxième module d'obtention 5B. Nous allons maintenant décrire, en référence à la figure 4C, les principales étapes mises en œuvre pour obtenir ce modèle MOD2. Ces étapes sont mises en œuvre ici par le module 5B mais elles peuvent être mises en œuvre par tout autre module, par exemple externe au dispositif d'estimation 5. Finally, as mentioned above, the invention is also based on a third numerical model MOD2 characterizing the aeronautical part. This model is stored in the database DB2 so as to be used by the second obtaining module 5B. We will now describe, with reference to FIG. 4C, the main steps implemented to obtain this MOD2 model. These steps are implemented here by the module 5B but they can be implemented by any other module, for example external to the estimation device 5.
Dans le mode de réalisation décrit ici, le modèle numérique MOD2 caractérisant la pièce aéronautique 2 fournit des coefficients d'atténuation massique du matériau composant la pièce aéronautique pour différentes valeurs d'énergie £} , j=l,...,K du faisceau de rayons émis par la source 3A.  In the embodiment described here, the numerical model MOD2 characterizing the aeronautical part 2 provides mass attenuation coefficients of the material composing the aeronautical part for different energy values £}, j = 1, ..., K of the beam of rays emitted by the source 3A.
Pour obtenir ce modèle numérique MOD2, le deuxième module d'obtention 5B obtient dans un premier temps la composition chimique du matériau composant la pièce aéronautique 2 (étape G10). Par composition chimique on entend les différents éléments chimiques composant le matériau ainsi que leurs proportions massiques dans le matériau. Cette composition chimique peut être fournie par exemple au deuxième module 5B par l'opérateur du système d'imagerie numérique 3 ou par toute autre personne via les moyens d'entrée/sortie 10 du dispositif d'estimation 5, ou lue par le module 5B dans un fichier préalablement renseigné.  In order to obtain this numerical model MOD2, the second obtaining module 5B first obtains the chemical composition of the material composing the aeronautical part 2 (step G10). By chemical composition is meant the different chemical elements composing the material as well as their mass proportions in the material. This chemical composition can be supplied for example to the second module 5B by the operator of the digital imaging system 3 or by any other person via the input / output means 10 of the estimation device 5, or read by the module 5B in a previously informed file.
Puis le module 5B obtient pour chaque élément chimique composant le matériau de la pièce aéronautique 2, un coefficient d'atténuation massique de cet élément pour lesdites différentes valeurs d'énergie _¾ , j=l,...,K (étape G20).  Then the module 5B obtains for each chemical element composing the material of the aeronautical part 2, a mass attenuation coefficient of this element for said different energy values _¾, j = 1, ..., K (step G20).
Dans le mode de réalisation décrit ici, l'étape G20 comprend la consultation d'une base de données NIST, comprenant, pour chaque élément chimique de la table périodique des éléments (aussi connu sous le nom de tableau de Mendeleïev), le coefficient d'atténuation massique de cet élément.  In the embodiment described here, step G20 comprises consulting a NIST database, comprising, for each chemical element of the periodic table of elements (also known as the Mendeleev table), the coefficient of mass attenuation of this element.
Le module 5B calcule ensuite, pour lesdites différentes valeurs d'énergie £) , j=l,...,K, le coefficient d'atténuation massique du matériau de la pièce aéronautique 2 selon l'expression (étape G30) : The module 5B then calculates, for said different energy values)), j = 1,..., K, the mass attenuation coefficient of the material of the aeronautical part 2 according to the expression (step G30):
où M est un entier supérieur ou égal à 1 désignant le nombre d'éléments chimiques composant le matériau de la pièce 2, wm est la proportion massique de l'élément chimique indexé par m dans le matériau de la pièce 2 et μ{πι, ΕΪ) est le coefficient d'atténuation massique de l'élément chimique indexé par m pour l'énergie £}. where M is an integer greater than or equal to 1 denoting the number of chemical elements composing the material of the part 2, w m is the proportion mass of the chemical element indexed by m in the material of the part 2 and μ {πι, ΕΪ) is the mass attenuation coefficient of the chemical element indexed by m for the energy £}.
Les coefficients d'atténuation massique μ(£)), j=l,...,N ainsi obtenus constituent le modèle numérique MOD2 et sont stockés dans la base de données DB2 (étape G40).  The mass attenuation coefficients μ ()), j = 1,..., N thus obtained constitute the numerical model MOD2 and are stored in the database DB2 (step G40).
Nous allons maintenant décrire comment les trois modèles MOD3A, MOD3B et MOD2 sont utilisés conformément à l'invention par le dispositif d'estimation 5 pour déterminer les paramètres de réglage optimaux du système 3 d'imagerie numérique en vue d'obtenir une ou plusieurs images de la pièce aéronautique 2 pour réaliser sur celle-ci un contrôle non destructif. On note que ces modèles précités peuvent être calculés préalablement à leur utilisation et stockés, comme décrit précédemment dans les bases de données DB3A, DB3B et DB2 accessibles par le dispositif d'estimation 5. L'accès aux bases de données en vue d'obtenir ces modèles (ou tout ou partie des valeurs stockées dans ces modèles) constitue alors une étape d'obtention de ces modèles au sens de l'invention.  We will now describe how the three models MOD3A, MOD3B and MOD2 are used according to the invention by the estimation device 5 to determine the optimal adjustment parameters of the digital imaging system 3 in order to obtain one or more images. of the aeronautical part 2 to perform on it a non-destructive test. It is noted that these models can be calculated prior to their use and stored, as previously described in the databases DB3A, DB3B and DB2 accessible by the estimation device 5. Access to databases to obtain these models (or all or part of the values stored in these models) then constitute a step of obtaining these models within the meaning of the invention.
La figure 5 représente les principales étapes mises en oeuvre par le dispositif d'estimation 5 lors de l'étape D10-1 d'estimation automatique des paramètres de réglage optimaux du système 3 d'imagerie numérique.  FIG. 5 represents the main steps implemented by the estimation device 5 during the step D10-1 of automatic estimation of the optimal adjustment parameters of the digital imaging system 3.
Conformément à l'invention, les paramètres de réglage optimaux du système 3 d'imagerie numérique sont estimés automatiquement par le dispositif d'estimation 5 en optimisant le rapport contraste sur bruit des images numériques IM de la pièce aéronautique 2 acquises par le système 3 d'imagerie numérique. Cette optimisation est conduite sur une pluralité de valeurs distinctes des différents paramètres de réglage, les plages sur lesquelles varient les valeurs des paramètres de réglage étant configurées préalablement au niveau du dispositif d'estimation 5 (par exemple via les moyens d'entrée/sortie 10).  According to the invention, the optimal adjustment parameters of the digital imaging system 3 are automatically estimated by the estimation device 5 by optimizing the contrast-to-noise ratio of the digital images IM of the aeronautical part 2 acquired by the system 3. digital imaging. This optimization is carried out on a plurality of distinct values of the various adjustment parameters, the ranges on which the values of the adjustment parameters are previously configured at the level of the estimation device 5 (for example via the input / output means 10 ).
Comme chaque paramètre de réglage peut prendre une pluralité de valeurs distinctes, dans le mode de réalisation décrit ici, pour optimiser le rapport contraste sur bruit, le dispositif d'estimation 5 met en uvre une pluralité d'itérations balayant tout ou partie de cette pluralité de valeurs (cf. étape H 10 de mise en oeuvre d'une nouvelle itération et étape test H60 de fin du processus itératif). Cette pluralité d'itérations est mise en œuvre via la prévision de plusieurs boucles (quatre boucles ici), imbriquées les unes dans les autres et balayant respectivement chacune différentes valeurs d'un paramètre de réglage distinct, i.e. une boucle balayant les valeurs de la tension d'opération, une boucle balayant les valeurs du courant, une boucle balayant les valeurs du temps d'exposition et une boucle balayant les valeurs de l'épaisseur du filtre externe. Since each setting parameter can take a plurality of distinct values, in the embodiment described here, to optimize the contrast-to-noise ratio, the estimation device 5 implements a plurality of iterations scanning all or part of this plurality. of values (see step H 10 of implementation of a new iteration and step H60 end of the iterative process). This plurality of iterations is implemented via the prediction of several loops (four loops here), nested in each other and respectively scanning each different values of a separate adjustment parameter, ie a loop scanning the values of the voltage of operation, a loop scanning the current values, a loop scanning the exposure time values and a loop scanning the values of the external filter thickness.
A chaque itération /ter correspondant à une valeur vt de la tension d'opération, à une valeur i du courant appliquée à la source 3A, à un temps d'exposition t de la pièce aéronautique et à une épaisseur Mfext de filtre du système 3 d'imagerie numérique (étape H 10), le dispositif d'estimation 5 accède via ses modules d'obtention 5A et 5B aux modèles numériques MOD3A, OD3B et MOD2 stockés respectivement dans les bases de données DB3A, DB3B et DB2. Plus précisément, les modules 5A et 5B extraient ici des bases de données DB3A, DB3B et DB2 les valeurs de ces modèles correspondant à la valeur de tension d'opération de la source 3A (étape H20), à la valeur de l'angle d'inclinaison de l'anode de la source 3A et à l'épaisseur et au type du scintillateur du détecteur 3B. At each iteration / ter corresponding to a value v t of the operating voltage, to a value i of the current applied to the source 3A, to an exposure time t of the aeronautical part and to a filter thickness M fext of digital imaging system 3 (step H 10), the estimation device 5 accesses via its obtaining modules 5A and 5B to the numerical models MOD3A, OD3B and MOD2 stored respectively in the databases DB3A, DB3B and DB2. More precisely, the modules 5A and 5B extract here DB3A, DB3B and DB2 databases the values of these models corresponding to the operating voltage value of the source 3A (step H20), to the value of the angle d tilt of the anode of the source 3A and the thickness and type of the scintillator of the detector 3B.
Puis il évalue de façon analytique, via son module d'évaluation 5C, à partir des valeurs extraites, un rapport contraste sur bruit noté CNR pour l'itération courante iter selon l'expression (étape H30) :  Then it evaluates analytically, via its evaluation module 5C, from the extracted values, a contrast-to-noise ratio noted CNR for the iteration current iter according to the expression (step H30):
CNR(iter) = . : ou alternativement selon l'expression : CNR (iter) =. : or alternatively according to the expression:
Sdet(.Vi, L,£) + Sdet(yi, L) S det (.V i, L, £) + S det (y i, L)
CNR(iter) = où SaetiVi, L), respectivement Sdet(yit L, E), désigne la valeur de signal vue par le détecteur 3B du système d'imagerie numérique 3 lorsqu'une épaisseur L, respectivement une épaisseur L - ε, de la pièce aéronautique 2 est traversée par le faisceau de rayons émis par la source 3A. Cette valeur de signal reflète, comme mentionné précédemment, un niveau de gris vu par le détecteur 3B. Autrement dit, l'invention s'appuie sur un simple calcul scalaire du rapport contraste sur bruit CNR. Ce calcul ne dépend que d'une seule grandeur, à savoir ici une valeur de signal représentative d'un niveau de gris vu par le détecteur 3B et correspondant à une épaisseur donnée de la pièce aéronautique traversée par le faisceau de rayons émis par la source. Cette valeur de signal dépend uniquement de l'épaisseur considérée (pour des paramètres de réglage donnés du système d'imagerie, ex. à Vi constante). L'épaisseur ε est prédéterminée et correspond à l'épaisseur d'un indicateur recherché durant le contrôle non destructif dans la pièce aéronautique 2. On note que dans le mode de réalisation décrit ici, cette épaisseur est un nombre réel qui en fonction de l'indicateur recherché peut être tantôt positif tantôt négatif, comme décrit précédemment. Un tel indicateur est par exemple une surépaisseur, une sous-épaisseur, une porosité ou encore une inclusion. L'épaisseur ε est fournie au préalable au dispositif d'estimation 5 par exemple par l'opérateur du système d'imagerie numérique 3 ou celui du dispositif de contrôle non destructif 4. CNR (iter) = where SaetiVi, L), respectively S det (y it L, E), designates the signal value seen by the detector 3B of the digital imaging system 3 when a thickness L, respectively a thickness L - ε, the aeronautical part 2 is crossed by the beam of rays emitted by the source 3A. This signal value reflects, as mentioned above, a gray level seen by the detector 3B. In other words, the invention is based on a simple scalar calculation of the contrast-to-noise ratio CNR. This calculation depends only on a single quantity, namely here a signal value representative of a gray level seen by the detector 3B and corresponding at a given thickness of the aeronautical part traversed by the beam of rays emitted by the source. This signal value depends solely on the thickness considered (for given adjustment parameters of the imaging system, eg at constant Vi). The thickness ε is predetermined and corresponds to the thickness of an indicator sought during the non-destructive testing in the aeronautical part 2. It is noted that in the embodiment described here, this thickness is a real number which, depending on the The desired indicator may be sometimes positive or negative, as described above. Such an indicator is, for example, an oversize, a sub-thickness, a porosity or an inclusion. The thickness ε is supplied in advance to the estimation device 5, for example by the operator of the digital imaging system 3 or that of the non-destructive inspection device 4.
Les signaux considérés dans le calcul du rapport contraste-sur- bruit sont caractérisés par les paramètres physiques du système d'imagerie numérique 3, autrement dit, du système réel.  The signals considered in the calculation of the contrast-on-noise ratio are characterized by the physical parameters of the digital imaging system 3, in other words, of the real system.
Plus particulièrement, les valeurs de signal vues par le détecteur sont obtenues ici analytiquement par le module d'évaluation 5C via les expressions analyti ues suivantes :  More particularly, the signal values seen by the detector are obtained here analytically by the evaluation module 5C via the following analytic terms:
— lQ(Ej, Vi) désigne l'intensité des photons présents dans un faisceau de rayons d'énergie E émis par la source 3A (correspondant à une longueur d'onde λ}) pour la tension d'opération vt ; - the Q (Ej, Vi) denotes the intensity of the photons present in an energy beam E emitted by source 3A (corresponding to a wavelength λ}) to the operating voltage v t;
— A(Ej. L) et Ae(Ej, L, s) désignent des facteurs d'atténuation subie par le faisceau de rayons d'énergie £) émis par la source 3A ; et A (Ej.L) and A e (Ej, L, s) denote attenuation factors experienced by the beam of energy rays ( e ) emitted by the source 3A; and
— R Ej) désigne la réponse spectrale du détecteur pour l'énergie E} (extraite du modèle MOD3B lors de l'étape H20). - R Ej) denotes the spectral response of the detector for energy E } (extracted from model MOD3B during step H20).
L'intensité i0(,Ej, Vt) est calculée par le module d'évaluation 5C pour j=l,...,N à partir des valeurs extraites du modèle MOD3A à l'étape H20. Plus précisément le module d'évaluation 5C calcule la valeur /0(£), vt) selon : The intensity i 0 ( , Ej, V t ) is calculated by the evaluation module 5C for j = l, ..., N from the values extracted from the model MOD3A at step H20. More precisely, the evaluation module 5C calculates the value / 0 (£), v t ) according to:
Io(Ej, Vi) = L t. l0s(yi, Ej) Io (Ej, V i ) = L t. l 0s (y i , Ej)
où /0s (Vi. Ej) est extraite du modèle MOD3A lors de l'étape H20, i désigne la valeur du courant appliquée à la source 3A correspondant à l'itération courante iter, et t désigne le temps d'exposition de la pièce aéronautique 2 correspondant à l'itération courante iter. where / 0s (Vi.Ej) is extracted from the MOD3A model in step H20, i denotes the value of the current applied to the source 3A corresponding to the current iteration iter, and t denotes the exposure time of the part aeronautical 2 corresponding to the current iteration iter.
Par ailleurs, dans le mode de réalisation décrit ici, le module d'évaluation 5C évalue le facteur d'atténuation A(£J( z,) de la façon suivante : On the other hand, in the embodiment described here, the evaluation module 5C evaluates the attenuation factor A (£ J ( z)) as follows:
A(£), L) = Afint(Ej, Lfint)AfexC(Ej, Lfext Amat(Ej, L) OÙ : A (£), L) = Afi nt (Ej, Lf int ) Af exC (Ej, Lf ext A mat (Ej, L) where:
— Afint{ j, lfint) désigne l'atténuation subie par le faisceau de rayons d'énergie E} émis par la source 3A due à un filtrage inhérent de la source 3A ; - Afint {j, l f int) designates the attenuation undergone by the beam of energy beams E } emitted by the source 3A due to an inherent filtering of the source 3A;
— Afext{Ej, Lfext) désigne l'atténuation subie par le faisceau de rayons d'énergie E} émis par la source 3A due à un filtrage externe à la source et correspondant à l'épaisseur de filtre ALfext correspondant à l'itération courante iter ; - Afext {Ej, Lfext) denotes the attenuation undergone by the beam of energy beams E } emitted by the source 3A due to an external filtering at the source and corresponding to the filter thickness AL fext corresponding to the iteration current iter;
— Amat{Ej, L) désigne l'atténuation subie par le faisceau de rayons d'énergie Ej émis par la source 3A due au matériau composant la pièce aéronautique 2.  - Amat {Ej, L) denotes the attenuation undergone by the beam of energy rays Ej emitted by the source 3A due to the material composing the aeronautical part 2.
Le module d'évaluation 5C évalue en outre le facteur d'atténuation A(Ej, L, ε) de la façon suivante : The evaluation module 5C further evaluates the attenuation factor A (E j , L, ε) as follows:
A(Ej, L, ε) = Afint(Ej, Lfint}Afext(Ejt Lfext)Amat(Ej, L - e)Aind(Ej, ε) dans lequel Aini E}, e) désigne l'atténuation subie par le faisceau de rayons d'énergie £) émis par la source 3A due à l'indicateur recherché d'épaisseur ε . A (E j , L, ε) = A fint (Ej, L fint } A fext (E jt L fext ) A mat (Ej, L - e) A ind (Ej, ε) in which A ini E } , e ) designates the attenuation undergone by the beam of energy rays (e) emitted by the source 3A due to the desired thickness indicator ε.
Chacune des atténuations précitées est modélisées ici par un filtrage opéré par un matériau distinct (élémentaire ou multi-éléments), à savoir le matériau constituant le filtre inhérent à la source 3A, le matériau constituant le filtre externe à la source 3A, le matériau de la pièce aéronautique 2 et le matériau correspondant à l'indicateur recherché (ce matériau pouvant être du vide, par exemple dans le cas d'une porosité).  Each of the abovementioned attenuations is modeled here by a filtering operated by a distinct material (elementary or multi-element), namely the material constituting the filter inherent in the source 3A, the material constituting the external filter at source 3A, the material of the aeronautical part 2 and the material corresponding to the desired indicator (this material may be vacuum, for example in the case of porosity).
Pour calculer ces différentes atténuations, le module d'évaluation 5C utilise la loi de Beer Lambert, définie par : To calculate these different attenuations, the 5C evaluation module uses Beer Lambert's law, defined by:
où : or :
— I désigne un indice prenant les valeurs fint.fext. mat ou ind suivant l'atténuation considérée ;  I denotes an index taking the values fint.fext. mat or ind depending on the attenuation considered;
— i désigne la densité du matériau associé à l'atténuation considérée ;I denotes the density of the material associated with the attenuation considered;
— L, désigne l'épaisseur du matériau associé à l'atténuation considérée (pour le matériau de la pièce aéronautique 2, cette épaisseur est prise égale à L ou à L- ε, et pour le filtre externe à la source 3A, cette épaisseur est prise égale à la valeur ALfext correspondant à l'itération iter courante) ; et - L, denotes the thickness of the material associated with the attenuation considered (for the material of the aeronautical part 2, this thickness is taken equal to L or L-ε, and for the external filter at the source 3A, this thickness is taken equal to the value AL fext corresponding to the iteration iter current); and
— iii(Ej) désigne le coefficient d'atténuation massique du matériau associé à l'atténuation considérée en fonction de l'énergie £). Dans le cas du matériau de la pièce aéronautique 2, ce coefficient d'atténuation massique est obtenu du modèle numérique MOD2 lors de l'étape H20. Pour les autres matériaux, on suppose que leurs coefficients d'atténuation massique sont fournis au dispositif 5 d'estimation, de même que leurs densités.  - iii (Ej) denotes the mass attenuation coefficient of the material associated with the attenuation considered as a function of the energy)). In the case of the material of the aeronautical part 2, this mass attenuation coefficient is obtained from the numerical model MOD2 during the step H20. For the other materials, it is assumed that their mass attenuation coefficients are provided to the estimator as well as their densities.
Une fois le rapport contraste-sur-bruit CNR(/te/) calculé par le module d'évaluation 5C pour l'itération courante iter, le module 5D de détermination 5D du dispositif 5 d'estimation compare le rapport CNR(/fe/) avec un seuil CNROPT (étape test H40), initialisé à 0 lors d'une étape préliminaire d'initialisation H00 mise en œuvre avant de démarrer les itérations.  Once the contrast-to-noise ratio CNR (/ te /) calculated by the evaluation module 5C for the current iteration iter, the 5D determination module 5D of the estimation device 5 compares the ratio CNR (/ fe / ) with a CNROPT threshold (test step H40), initialized to 0 during a preliminary H00 initialization step implemented before starting the iterations.
Si le rapport contraste-sur-bruit est supérieur au seuil CNROPT courant (réponse oui à l'étape test H40), le seuil CNROPT est mis à jour avec la valeur du rapport CNR(/te/) (étape H50). L'itération courante iter est par ailleurs stockée dans un indicateur iteropt  If the contrast-on-noise ratio is greater than the current CNROPT threshold (yes answer to the test step H40), the CNROPT threshold is updated with the value of the CNR ratio (/ te /) (step H50). The current iteration is also stored in an iteropt indicator
Si le rapport contraste-sur-bruit est supérieur au seuil CNROPT courant (réponse oui à l'étape test H40), le seuil CNROPT est gardé inchangé.  If the contrast-on-noise ratio is greater than the current CNROPT threshold (yes answer to the test step H40), the CNROPT threshold is kept unchanged.
Puis le module de détermination 5D vérifie si toutes les itérations ont été mises en œuvre (étape test H60).  Then the determination module 5D checks whether all the iterations have been implemented (test step H60).
Si il reste des itérations (réponse oui à l'étape test H60), une nouvelle itération correspondant à de nouvelles valeurs des paramètres de réglage vl r '\, t et Mfext est mise en œuvre et les étapes H20 à H60 sont réitérées pour ces nouvelles valeurs. Sinon (réponse non à l'étape test H60), le module de détermination 5D détermine que les valeurs optimales des paramètres de réglage du système 3 d'imagerie numérique correspondent aux valeurs des paramètres de réglage testées lors de l'itération iteropt (étape H70). If there are still iterations (yes answer to the test step H60), a new iteration corresponding to new values of the adjustment parameters v lr ' \, t and M fext is implemented and the steps H20 to H60 are repeated for these new values. If not (answer no to the test step H60), the determination module 5D determines that the optimum values of the adjustment parameters of the digital imaging system 3 correspond to the values of the adjustment parameters tested during the iteropt iteration (step H70 ).
Ces valeurs optimales sont alors fournies à l'opérateur du système 3 d'imagerie numérique pour être appliquées sur le système 3 d'imagerie numérique lors de l'étape D10-2.  These optimal values are then provided to the operator of the digital imaging system 3 to be applied to the digital imaging system 3 in step D10-2.
On note que les exemples numériques donnés dans la description ne sont donnés qu'à titre illustratif. Bien entendu, d'autres plages de valeurs peuvent être envisagées, en fonction notamment du système 3 d'imagerie numérique et de la pièce aéronautique 2 considérée.  Note that the numerical examples given in the description are given for illustrative purposes only. Of course, other ranges of values can be envisaged, in particular according to the digital imaging system 3 and the aeronautical part 2 considered.
L'invention fournit ainsi un processus automatique permettant de déterminer automatiquement les paramètres optimaux d'un système d'imagerie numérique tel un système radiographique ou tomographique pour l'acquisition d'images d'une pièce aéronautique. Il convient de noter que bien que décrite en référence à une pièce aéronautique, l'invention s'applique à d'autres pièces pouvant faire l'objet de contrôle non destructif.  The invention thus provides an automatic process for automatically determining the optimal parameters of a digital imaging system such as a radiographic or tomographic system for acquiring images of an aeronautical part. It should be noted that although described with reference to an aeronautical part, the invention applies to other parts that can be subject to non-destructive testing.

Claims

REVENDICATIONS
1. Procédé de contrôle non destructif réalisé sur une pièce aéronautique (2) comprenant : A non-destructive testing method performed on an aeronautical part (2) comprising:
— 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 ; An adjustment step (D10) of a digital imaging system, said digital imaging system (3) being a radiography or digital tomography system comprising a source (3A) of electromagnetic rays and a detector (3B) able to detect the electromagnetic rays emitted by the source;
— une étape d'acquisition (D20) d'au moins une image numérique (I ) de la pièce aéronautique au moyen du système de radiographie ou de tomographie numérique réglé ; et A step of acquiring (D20) at least one digital image (I) of the aeronautical part by means of the set radiography or digital tomography system; and
— 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 :  A non-destructive testing step (D30) performed on the aeronautical part using said at least one acquired digital image; said method being characterized in that it further comprises a preliminary step of estimating (D10-1) adjustment parameters of the digital imaging system, this estimation step comprising:
— 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 modèle (MOD3A) caractérisant la source et un modèle (MOD3B) caractérisant le détecteur ;  A step of obtaining (H20) a digital model characterizing the digital imaging system, said digital model comprising a model (MOD3A) characterizing the source and a model (MOD3B) characterizing the detector;
— une étape d'obtention (H20) d'un modèle numérique (MOD2) caractérisant la pièce aéronautique ;  A step of obtaining (H20) a numerical model (MOD2) characterizing the aeronautical part;
— une étape d'évaluation (H30), pour une pluralité de valeurs distinctes 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  An evaluation step (H30), for a plurality of distinct values of the adjustment parameters, of a contrast-to-noise ratio obtained for the digital imaging system and said aeronautical part, this evaluation step using the numerical models obtained characterizing the digital imaging system and the aeronautical part; and
— une étape de détermination automatique (H70) de valeurs dites optimales des 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.  A step of automatic determination (H70) of so-called optimum values of the adjustment parameters optimizing the contrast-to-noise ratio, these optimum values being applied to the digital imaging system during the adjustment step.
2. Procédé selon la revendication 1 dans lequel le rapport contraste-sur-bruit est évalué, pour ladite pluralité de valeurs distinctes des paramètres de réglage, à partir de signaux caractérisés par des paramètres physiques du système d'imagerie numérique. The method of claim 1 wherein the contrast-to-noise ratio is evaluated for said plurality of distinct values. setting parameters, based on signals characterized by physical parameters of the digital imaging system.
3. Procédé selon la revendication 1 ou 2 dans lequel le modèle (MOD3A) caractérisant la source et le modèle (MOD3B) caractérisant le détecteur fournissent des réponses en énergie respectivement de la source et du détecteur en fonction d'une longueur d'onde. 3. Method according to claim 1 or 2 wherein the model (MOD3A) characterizing the source and the model (MOD3B) characterizing the detector provide energy responses respectively of the source and the detector as a function of a wavelength.
4. Procédé selon l'une quelconque des revendications 1 à 3 dans lequel les paramètres de réglage du système comprennent au moins un paramètre parmi : 4. Method according to any one of claims 1 to 3 wherein the system setting parameters comprise at least one of:
— une tension d'opération de la source du système d'imagerie numérique ;  An operating voltage of the source of the digital imaging system;
— 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 ;  A time of exposure of the aeronautical part to a ray beam emitted by the source of the digital imaging system;
— une valeur de courant appliquée à la source du système d'Imagerie numérique ; et  A current value applied to the source of the digital imaging system; and
— une épaisseur de filtre appliqué dans le système d'imagerie numérique à la source.  - a filter thickness applied in the digital imaging system at the source.
5. Procédé selon l'une quelconque des revendications 1 à 4 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. The method of any one of claims 1 to 4 wherein the model characterizing the source (MOD3A) provides emission spectra of the source for different operation voltage values of the source.
6. Procédé selon la revendication 5 dans lequel le modèle ( OD3A) 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. The method of claim 5 wherein the source model (OD3A) further provides estimates of at least two half-attenuation layers of the source for different operating voltage values of the source, said method comprising a validation step (E50) of the source characterizing model compares the estimates of said at least two half-attenuation layers of the source with experimental results.
7. Procédé selon l'une quelconque des revendications 1 à 6 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. 7. Method according to any one of claims 1 to 6 wherein the model characterizing the detector (MOD3B) provides a spectral response of the detector to a ray beam emitted by the source.
8. Procédé selon la revendication 7 comprenant en outre une étape de validation (F30) du modèle caractérisant le détecteur comprenant : 8. The method of claim 7 further comprising a validation step (F30) of the model characterizing the detector comprising:
— une estimation d'au moins un niveau de gris théorique 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 ;  An estimate of at least one theoretical gray level seen by the detector for at least one given operating voltage of the source, by using the model characterizing the source and the model characterizing the detector;
— 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 source ;  Obtaining a measurement of at least one gray level seen by the detector for said at least one given operating voltage of the source;
— une comparaison dudit au moins un niveau de gris théorique estimé avec ledit au moins un niveau de gris mesuré.  A comparison of said at least one estimated theoretical gray level with said at least one measured gray level.
9. Procédé selon l'une quelconque des revendications 1 à 8 dans lequel le modèle numérique (MOD2) caractérisant la pièce aéronautique9. Method according to any one of claims 1 to 8 wherein the numerical model (MOD2) characterizing the aeronautical part
(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. (2) provides mass attenuation coefficients of a material composing the aeronautical part for different energy values of a ray beam emitted by the source.
10. Procédé selon la revendication 8 dans lequel l'obtention du modèle numérique caractérisant la pièce aéronautique comprend : 10. The method of claim 8 wherein obtaining the numerical model characterizing the aeronautical part comprises:
— 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 ;  A step of obtaining (G10) a chemical composition of the material composing the part, said material being composed of at least one chemical element;
— 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 Iesdites différentes valeurs d'énergie ; et For each chemical element composing the material, a step of obtaining (G20) a mass attenuation coefficient of this element for said different energy values; and
— une étape de calcul (G30) des coefficients d'atténuation massiques du matériau pour Iesdites 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.  A step of calculating (G30) the mass attenuation coefficients of the material for said different energy values from the mass attenuation coefficients of the elements composing the material for these different energies and the proportions of these elements in the material.
11. Procédé selon la revendication 10 dans lequel le coefficient d'atténuation massique de chaque élément du matériau est obtenu en consultant une base de données NIST. 11. The method of claim 10 wherein the mass attenuation coefficient of each element of the material is obtained by consulting a NIST database.
12. Procédé selon l'une quelconque des revendications 1 à 11 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 : 12. Method according to any one of claims 1 to 11 wherein during the evaluation step (H30), the contrast ratio noise noted CNR obtained for the digital imaging system and said aeronautical part is evaluated for different values of the adjustment parameters according to the expression:
_ Sdet(L, )— Sdet(L) ou : S det (L,) - S det (L) or:
V¾«a) + S**(L, e) V¾ "a) + S ** (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 Ι - ε, de la pièce aéronautique est traversée par le faisceau de rayons émis par la source, ε 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 analytiquement par : where s det (L), respectively s det (L, e), denotes a signal value seen by the detector of the digital imaging system when a thickness L, respectively a thickness Ι - ε, of the aeronautical part is traversed by the beam of rays emitted by the source, ε denoting a thickness of a desired indicator during non-destructive testing in the aeronautical part, said signal values being analytically defined by:
= V)A(Ej, L)R(Ej) = V) A (Ej, L) R (Ej)
j  j
et and
I0(Ej, V)A£(Ej, L, e)R ÎEj) I 0 (Ej, V) A £ (Ej, L, e) R EEj)
j  j
où : or :
— l0(Ej, V) désigne l'intensité des photons présents dans un faisceau de rayons d'énergie Ej émis par la source correspondant à une longueur d'onde Ay, 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 ; - l 0 (Ej, V) designates the intensity of photons present in a beam of energy rays Ej emitted by the source corresponding to a wavelength A y , said intensity depending on the voltage V of operation of the source, the exposure time of the aeronautical part to the ray beam emitted by the source, and the current value applied to the source;
— A{Ej, L) désigne un facteur d'atténuation subie par le faisceau de rayons émis par la source d'é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 et de l'épaisseur L ; — AE(Ej, L, e) 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, de l'épaisseur L, et de l'épaisseur ε de l'indicateur ; et - A {Ej, L) denotes an attenuation factor experienced by the beam of rays emitted by the energy source E } -, this attenuation factor depending on a mass attenuation coefficient of the material constituting the part aeronautical and thickness L; - A E (Ej, L, e) denotes an attenuation factor experienced by the beam of rays emitted by the energy source E i , this attenuation factor depending on a mass attenuation coefficient of the material composing the aeronautical part, the thickness L, and the thickness ε of the indicator; and
— R(Ej) désigne la réponse spectrale du détecteur en fonction de l'énergie Ej.  - R (Ej) denotes the spectral response of the detector as a function of the energy Ej.
13. Procédé selon la revendication 12 dans lequel : 13. The method of claim 12 wherein:
— le facteur d'atténuation A(Ej, L) est défini par : The attenuation factor A (Ej, L) is defined by:
A(Ej, L)  A (Ej, L)
— le facteur d'atténuation A(Ej( L, ε) est défini par : The attenuation factor A (E j ( L, ε) is defined by:
A(Ej, L. ε) = Afint(Ej, Lfint)AfexC(Ej, Lfext)Amat(Ej, L— ε)ν ίη£ΐ(£), ε) où Arint(Ej, Lfint), Afext(Ejt Lfext), Amat(Ej, L) et Amat{Ej, L— ε) désignent respectivement l'atténuation subie par le faisceau de rayons émis par la source d'énergie Ej due à un filtrage inhérent à la source, l'atténuation subie par le faisceau de rayons émis par la source d'énergie £} due à un filtrage externe à la source, l'atténuation subie par le faisceau de rayons émis par la source d'énergie F,- due au matériau composant la pièce aéronautique et l'atténuation subie par le faisceau de rayons d'énergie E} émis par la source en présence de l'indicateur recherché d'épaisseur ε. A (Ej, L. ε) = Af int (Ej, Lf int ) Af exC (Ej, Lf ext ) A mat (Ej, L- ε) ν ίη £ ΐ (£), ε) where A rint (Ej, L fint ), A fext (E jt L fext ), A mat (Ej, L) and A mat {E j , L- ε) respectively denote the attenuation undergone by the beam of rays emitted by the energy source Ej due to a filtering inherent to the source, the attenuation undergone by the beam of rays emitted by the source of energy £} due to a filtering external to the source, the attenuation undergone by the beam of rays emitted by the source of energy F, - due to the material constituting the aeronautical piece and the attenuation undergone by the energy beam E} emitted by the source in the presence of the desired indicator ε thickness.
14. Procédé selon l'une quelconque des revendications 1 à 13 dans lequel l'étape d'estimation comprend la mise en uvre d'une pluralité d'itérations (iter) 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 source, le détecteur et la pièce aéronautique stockées préalablement dans des bases de données. The method according to any one of claims 1 to 13 wherein the estimating step comprises implementing a plurality of iterations scanning all or part of the plurality of distinct values of the setting parameters. , and during which, the evaluation of the contrast-on-noise ratio is carried out by accessing numerical model values characterizing the source, the detector and the aeronautical part previously stored in databases.
15. Système (1) de contrôle non destructif d'une pièce aéronautique comprenant : 15. System (1) for non-destructive testing of an aeronautical part comprising:
— 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 détecteur apte à détecter les rayons électromagnétiques émis par la source ; et A digital imaging system (3), this digital imaging system being an X-ray or digital tomography system comprising a source of electromagnetic radiation and a detector capable of detecting the electromagnetic rays emitted by the source; and
— 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 ;  A non-destructive testing device (4) for the aeronautical part configured to use at least one digital image of the aeronautical part acquired by the digital imaging system;
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 la pièce aéronautique lors d'un réglage du système d'imagerie numérique, ledit dispositif d'estimation comprenant : the system being characterized in that it further comprises a device (5) for estimating said optimum values of adjustment parameters of the digital imaging system, said optimum values of adjustment parameters being intended to be applied to said digital system. digital imaging prior to acquiring said at least one digital image of the aeronautical part during an adjustment of the digital imaging system, said estimating device comprising:
— un premier module d'obtention (5A) configuré pour obtenir un modèle numérique caractérisant le système d'imagerie numérique, ce modèle numérique comprenant un modèle caractérisant la source et un modèle caractérisant le détecteur ;  A first obtaining module (5A) configured to obtain a digital model characterizing the digital imaging system, this digital model comprising a model characterizing the source and a model characterizing the detector;
— un deuxième module d'obtention (5B), configuré pour obtenir un modèle numérique caractérisant la pièce aéronautique ;  A second obtaining module (5B) configured to obtain a digital model characterizing the aeronautical part;
— un module d'évaluation (5C) configuré pour évaluer, pour une 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  An evaluation module (5C) configured to evaluate, for a plurality of distinct values of the adjustment parameters, a contrast-to-noise ratio obtained for the digital imaging system and said aeronautical part, this evaluation module using the models obtained numerals characterizing the digital imaging system and the aeronautical part; and
— un module de détermination automatique (5D), configuré pour déterminer automatiquement des valeurs optimales des paramètres de réglage optimisant le rapport contraste sur bruit An automatic determination module (5D), configured to automatically determine optimum values of the adjustment parameters optimizing the contrast-to-noise ratio
EP17771806.1A 2016-09-02 2017-08-31 Non-destructive inspection method and system carried out on an aeronautical part Pending EP3507593A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1658197A FR3055728B1 (en) 2016-09-02 2016-09-02 METHOD AND SYSTEM FOR NON-DESTRUCTIVE CONTROL ON AERONAUTICAL COMPONENT
PCT/FR2017/052313 WO2018042132A1 (en) 2016-09-02 2017-08-31 Non-destructive inspection method and system carried out on an aeronautical part

Publications (1)

Publication Number Publication Date
EP3507593A1 true EP3507593A1 (en) 2019-07-10

Family

ID=57539392

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17771806.1A Pending EP3507593A1 (en) 2016-09-02 2017-08-31 Non-destructive inspection method and system carried out on an aeronautical part

Country Status (4)

Country Link
US (1) US10724970B2 (en)
EP (1) EP3507593A1 (en)
FR (1) FR3055728B1 (en)
WO (1) WO2018042132A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113205497B (en) * 2021-04-30 2022-06-07 江苏迪业检测科技有限公司 Image processing method, device, equipment and medium for double-wire type image quality meter
CN115061519B (en) * 2022-06-28 2022-11-25 中国飞机强度研究所 Multi-thermal-field simulation high-precision control system and method for testing strength of aerospace plane
CN116953356B (en) * 2023-09-20 2023-12-26 南京纳特通信电子有限公司 Ground-air integrated three-dimensional space radio frequency spectrum monitoring method and system

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2425415A1 (en) * 2000-10-11 2003-04-10 Randolf Hanke Method and apparatus for the representation of an object via a transmission as well as for reconstruction
US20070076842A1 (en) * 2005-09-30 2007-04-05 Tkaczyk John E Adaptable energy discriminating computed tomography system
DE102014008671A1 (en) * 2014-06-14 2015-12-17 Microvista Gmbh Method for optimally arranging an object in a device and device for displaying an internal spatial structure of the object

Also Published As

Publication number Publication date
US10724970B2 (en) 2020-07-28
FR3055728B1 (en) 2019-09-06
WO2018042132A1 (en) 2018-03-08
FR3055728A1 (en) 2018-03-09
US20190204241A1 (en) 2019-07-04

Similar Documents

Publication Publication Date Title
EP3149705B1 (en) Method and device for estimation of a quality index of a 3-d image of a piece of composite material
WO2018042132A1 (en) Non-destructive inspection method and system carried out on an aeronautical part
FR2948481A1 (en) IMAGING METHOD FOR THE PRODUCTION OF A TRIPLE ENERGY MODELING, AND DEVICE FOR IMPLEMENTING SUCH A METHOD
FR2668829A1 (en) DEVICE AND METHOD FOR NON - DESTRUCTIVE CONTROL WITH SIMULTANEOUS ACQUISITION OF RADIOGRAPHIC DATA AND TOMOGRAPHIC DATA.
EP2930540B1 (en) Method for linearisation of attenuation measurements taken by a spectrometric sensor
WO2013068687A1 (en) Method and device for estimating a degree of porosity of a sample of material on the basis of at least one image coded by grey levels
FR2966716A1 (en) METHOD OF ESTIMATING AND CORRECTION OF DIFFUSION IN MAMMOGRAPHY.
EP2726815B1 (en) Method and device for identifying a material by the spectral analysis of electromagnetic radiation passing through said material
FR2968403A1 (en) METHOD FOR EXTRACTING A PREDICTED DIFFUSION SPECTRUM
FR3071611A1 (en) METHOD AND DEVICE FOR NON-DESTRUCTIVE CONTROL OF AN AERONAUTICAL WORKPIECE
FR2927719A1 (en) METHOD FOR PROCESSING IMAGES OBTAINED BY TOMOSYNTHESIS AND APPARATUS THEREFOR
FR3106409A1 (en) Method for establishing a spectral response function of an X-ray or gamma-ray measurement system
EP3153888A1 (en) Method for correction of a spectrum
EP0965946B1 (en) Method of improving signal-to-noise ratio in the radiographic image of a moving object
EP2769209B1 (en) Method and device for characterising a material by electromagnetic radiation scattering
FR3091359A1 (en) Method for calibrating a spectrometric detector
EP3864624B1 (en) Method for correcting a spectral image
WO2023037078A1 (en) Method, device and computer program for monitoring a part by x-ray
FR3080700A1 (en) METHOD AND DEVICE FOR NON-DESTRUCTIVE CONTROL OF A PART
FR3064776A1 (en) METHOD FOR LEARNING THE DYNAMIC BEHAVIOR OF AN AIRCRAFT ENGINE
FR3118520A1 (en) Process for processing volumic images by principal component analysis
FR3136855A1 (en) Method for generating synthetic spectral data
FR3117597A1 (en) Non-destructive testing method for an aeronautical part
FR3108196A1 (en) Method and system for non-destructive testing of an aeronautical part
FR3090128A1 (en) METHOD FOR CHARACTERIZING A SPECTROMETER, COMPUTER PROGRAM PRODUCT AND RELATED COMPUTER

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20190326

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20211206