CA2163882C - Procede et dispositif pour determiner la fonction d'attenuation d'un objet pour ce qui est de la transmission de l'epaisseur d'un materiau de reference - Google Patents
Procede et dispositif pour determiner la fonction d'attenuation d'un objet pour ce qui est de la transmission de l'epaisseur d'un materiau de reference Download PDFInfo
- Publication number
- CA2163882C CA2163882C CA002163882A CA2163882A CA2163882C CA 2163882 C CA2163882 C CA 2163882C CA 002163882 A CA002163882 A CA 002163882A CA 2163882 A CA2163882 A CA 2163882A CA 2163882 C CA2163882 C CA 2163882C
- Authority
- CA
- Canada
- Prior art keywords
- energy
- function
- intensity
- spectra
- reference material
- 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.)
- Expired - Lifetime
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/06—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 measuring the absorption
- G01N23/083—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 measuring the absorption the radiation being X-rays
Landscapes
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
Abstract
La fonction de transmission de l'objet est exprimée comme un développement fini en puissances de la transmission d'une épaisseur de référence d'un matériau de référence, chaque puissance du développement étant égale à un rapport entre une épaisseur du matériau de référence et l'épaisseur de référence; les coefficients du développement son t déterminés grâce aux mesures de l'intensité transmise par les différentes épaisseur, choisies, incluant l'épaisseur nulle, du matériau de référence soumis à un faisceau de mesure présentant successivement plusieurs spectres en énergie et aux mesures de l'intensité transmise par l'objet soumis à ce même faisceau de mesure à spectre variable; la fonction d'atténuation est déduite de la détermination de la fonction de transmission. Application à la radioscopie.
Description
~O 94/28398 PCTlFR94/00581 PROCEDE ET DISPOSITIF POUR DETERMINER LA FONCTION D'ATTENUATION D'UN OBJET
PAR RAPPORT A LA TRANSMISSION D'UNE EPAISSEUR D'UN MATERIAU DE REFERENCE
La présente invention concerne un procédé pour déterminer la fonction d'atténuation d'un objet et un dispositif pour la mise en oeuvre du procédé. Elle s'applique notamment mais pas uniquement à la radioscopie. D'une manière générale, elle s'applique au domaine de i o la spectroscopie dans la gamme d'énergie allant de quelques keV à
quelques MeV.
De façon usuelle, on ne détermine pas la fonction d'atténuation d'un objet mais uniquement des valeurs de cette fonction pour des i s énergies données. La fonction d'atténuation peut étre ensuite reconstituée par interpolation d'une courbe passant par ces points.
On utilise le plus souvent des faisceaux monochromatiques de rayons X ou y d'énergie stable et déterminée, provenant de sources isotopiques dont l'énergie est éventuellement diminuée par effet
PAR RAPPORT A LA TRANSMISSION D'UNE EPAISSEUR D'UN MATERIAU DE REFERENCE
La présente invention concerne un procédé pour déterminer la fonction d'atténuation d'un objet et un dispositif pour la mise en oeuvre du procédé. Elle s'applique notamment mais pas uniquement à la radioscopie. D'une manière générale, elle s'applique au domaine de i o la spectroscopie dans la gamme d'énergie allant de quelques keV à
quelques MeV.
De façon usuelle, on ne détermine pas la fonction d'atténuation d'un objet mais uniquement des valeurs de cette fonction pour des i s énergies données. La fonction d'atténuation peut étre ensuite reconstituée par interpolation d'une courbe passant par ces points.
On utilise le plus souvent des faisceaux monochromatiques de rayons X ou y d'énergie stable et déterminée, provenant de sources isotopiques dont l'énergie est éventuellement diminuée par effet
2 o Compton ou issus de tubes à rayons X à spectre large filtré par un réseau de Bragg pour obtenir un spectre de raies.
L'objet à caractériser est placé sur le parcours du faisceau et un détecteur mesure, pour les différentes énergies, l'intensité du faisceau 2s transmis.
Si l'on pose que lo(E) représente l'intensité du faisceau incident d'énergie E, que It(E) représente l'intensité du faisceau transmis après traversée de l'objet, alors la transmission Tr(E) pour l'énergie E est égale au rapport It(E)/lo(E). L'atténuation considérée pour l'énergie E,
L'objet à caractériser est placé sur le parcours du faisceau et un détecteur mesure, pour les différentes énergies, l'intensité du faisceau 2s transmis.
Si l'on pose que lo(E) représente l'intensité du faisceau incident d'énergie E, que It(E) représente l'intensité du faisceau transmis après traversée de l'objet, alors la transmission Tr(E) pour l'énergie E est égale au rapport It(E)/lo(E). L'atténuation considérée pour l'énergie E,
3 o Att(E), est égale à - Log Tr(E).
Ces techniques connues présentent de nombreux inconvénients.
Les atténuations ne sont déterminées que pour un faible nombre d'énergies. En conséquence, si l'on veut déterminer la fonction 3 s atténuation, il faut ajuster une courbe passant par ces points ce qui rend la précision dépendante de la manière dont est réalisée l'interpolation.
FEUILLE DE REMPLACEMENT (REGL E 26) WO 94/28398 ' PCT/FR94/0058~
De plus, la détermination des valeurs de l'atténuation par ces techniques est lente si l'on désire obtenir une précision de l'ordre de 1 %. En effet, il est nécessaire de moyenner un grand nombre de mesures et un temps de l'ordre d'une seconde est nécessaire pour s déterminer une valeur de l'atténuation.
Ce dernier point interdit l'utilisation de ces techniques pour la réalisation d'images, et plus encore la réalisation d'images en temps quasi-réel, impérative dans les applications touchant à l'inspection i o des bagages et la sécurité.
Une image comprend environ 500 points par ligne et environ 600 lignes; si chaque mesure nécessite une seconde, la réalisation d'une image, pour une énergie unique, demande environ 83h, à multiplier i s par le nombre voulu d'énergies.
D'autre part, la diffraction de Bragg et l'effet Compton sont dépendants de l'angle d'incidence du faisceau et ne permettent d'obtenir un faisceau à l'énergie désirée que pour un angle d'incidence donné
20 obtenu avec un faisceau parallèle. Avec ce type de faisceau collimaté, il n'est possible de faire la mesure de la transmission qu'en un point de l'objet . Pour obtenir simultanément un grand nombre de points de mesure, il faut multiplier les moyens d'émission de faisceau, car étant donnée la dépendance angulaire, l'utilisation d'un faisceau en 2s éventail est interdite. On comprend que la mise en oeuvre d'un tel montage est très complexe et coûteuse.
L'utilisation de sources isotopiques est plus simple mais ces sources présentent aussi des inconvénients. En effet, pour des raisons de 3 o sécurité, la manipulation des sources est sévèrement réglementée et l'intensité des faisceaux est faible.
Donc, les techniques connues de mesure d'atténuation ne permettent des évaluations que pour des énergies particulières dépendant 3 s fortement du dispositif utilisé et non de la volonté ou des besoins de O 94/28398 ~ ~ ~ ~ ~ ~ PCT/FR94/00581 l'utilisateur. D'autre part, elles nécessitent un temps de mesure très long, sont coûteuses et délicates à mettre en oeuvre; par conséquent, ces techniques sont difficilement utilisables en radioscopie, ou en imagerie de sécurité qui sont des domaines importants d'application s de la détermination de l'atténuation d'un objet.
La présente invention pallie ces inconvénients. Elle permet de déterminer la fonction d'atténuation de manière très rapide (en un temps de l'ordre de la dizaine de millisecondes) en utilisant des tubes i o à rayons X ou des générateurs de rayons y standards. La fonction d'atténuation est déterminée en fonction d'une variable en relation bijective avec l'énergie. Comme on va le voir plus loin, la transmission u d'une épaisseur de référence d'un matériau de référence est choisie comme variable. Si l'on désire exprimer l'atténuation de l'objet en i s fonction de l'énergie, il suffit de mesurer une fois pour toute (et non à
chaque détermination) la transmission u pour les énergies souhaitées, ceci à l'aide des techniques de l'art antérieur.
De manière plus précise, l'invention concerne un procédé pour 2o déterminer la fonction d'atténuation d'un objet, comprenant les étapes suivantes:
A) Exprimer la fonction de transmission de l'objet comme un développement fini en puissances de la transmission u d'une épaisseur de référence d'un matériau de référence, chaque 2s puissance de u étant affectée d'un coefficient, les puissances du développement étant égales à un rapport entre des épaisseurs préalablement déterminées du matériau de référence et l'épaisseur de référence;
3 o B) Déterminer les coefficients du développement:
a) en mesurant l'intensité d'un faisceau de mesure pour différents spectres en énergie, b) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par chacune des épaisseurs 3s considérées de matériau de référence, WO 94/28398 ~ ~ ~ ~ PCT/FR94/0058~
Ces techniques connues présentent de nombreux inconvénients.
Les atténuations ne sont déterminées que pour un faible nombre d'énergies. En conséquence, si l'on veut déterminer la fonction 3 s atténuation, il faut ajuster une courbe passant par ces points ce qui rend la précision dépendante de la manière dont est réalisée l'interpolation.
FEUILLE DE REMPLACEMENT (REGL E 26) WO 94/28398 ' PCT/FR94/0058~
De plus, la détermination des valeurs de l'atténuation par ces techniques est lente si l'on désire obtenir une précision de l'ordre de 1 %. En effet, il est nécessaire de moyenner un grand nombre de mesures et un temps de l'ordre d'une seconde est nécessaire pour s déterminer une valeur de l'atténuation.
Ce dernier point interdit l'utilisation de ces techniques pour la réalisation d'images, et plus encore la réalisation d'images en temps quasi-réel, impérative dans les applications touchant à l'inspection i o des bagages et la sécurité.
Une image comprend environ 500 points par ligne et environ 600 lignes; si chaque mesure nécessite une seconde, la réalisation d'une image, pour une énergie unique, demande environ 83h, à multiplier i s par le nombre voulu d'énergies.
D'autre part, la diffraction de Bragg et l'effet Compton sont dépendants de l'angle d'incidence du faisceau et ne permettent d'obtenir un faisceau à l'énergie désirée que pour un angle d'incidence donné
20 obtenu avec un faisceau parallèle. Avec ce type de faisceau collimaté, il n'est possible de faire la mesure de la transmission qu'en un point de l'objet . Pour obtenir simultanément un grand nombre de points de mesure, il faut multiplier les moyens d'émission de faisceau, car étant donnée la dépendance angulaire, l'utilisation d'un faisceau en 2s éventail est interdite. On comprend que la mise en oeuvre d'un tel montage est très complexe et coûteuse.
L'utilisation de sources isotopiques est plus simple mais ces sources présentent aussi des inconvénients. En effet, pour des raisons de 3 o sécurité, la manipulation des sources est sévèrement réglementée et l'intensité des faisceaux est faible.
Donc, les techniques connues de mesure d'atténuation ne permettent des évaluations que pour des énergies particulières dépendant 3 s fortement du dispositif utilisé et non de la volonté ou des besoins de O 94/28398 ~ ~ ~ ~ ~ ~ PCT/FR94/00581 l'utilisateur. D'autre part, elles nécessitent un temps de mesure très long, sont coûteuses et délicates à mettre en oeuvre; par conséquent, ces techniques sont difficilement utilisables en radioscopie, ou en imagerie de sécurité qui sont des domaines importants d'application s de la détermination de l'atténuation d'un objet.
La présente invention pallie ces inconvénients. Elle permet de déterminer la fonction d'atténuation de manière très rapide (en un temps de l'ordre de la dizaine de millisecondes) en utilisant des tubes i o à rayons X ou des générateurs de rayons y standards. La fonction d'atténuation est déterminée en fonction d'une variable en relation bijective avec l'énergie. Comme on va le voir plus loin, la transmission u d'une épaisseur de référence d'un matériau de référence est choisie comme variable. Si l'on désire exprimer l'atténuation de l'objet en i s fonction de l'énergie, il suffit de mesurer une fois pour toute (et non à
chaque détermination) la transmission u pour les énergies souhaitées, ceci à l'aide des techniques de l'art antérieur.
De manière plus précise, l'invention concerne un procédé pour 2o déterminer la fonction d'atténuation d'un objet, comprenant les étapes suivantes:
A) Exprimer la fonction de transmission de l'objet comme un développement fini en puissances de la transmission u d'une épaisseur de référence d'un matériau de référence, chaque 2s puissance de u étant affectée d'un coefficient, les puissances du développement étant égales à un rapport entre des épaisseurs préalablement déterminées du matériau de référence et l'épaisseur de référence;
3 o B) Déterminer les coefficients du développement:
a) en mesurant l'intensité d'un faisceau de mesure pour différents spectres en énergie, b) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par chacune des épaisseurs 3s considérées de matériau de référence, WO 94/28398 ~ ~ ~ ~ PCT/FR94/0058~
4 c) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmis par l'objet, d) pour chacun desdits spectres, en décomposant l'intensité du faisceau transmis par l'objet en une somme des intensités s mesurées à l'étape B-b) relatives au spectre considéré, ces intensités étant affectées des coefficients recherchés;
e) en résolvant le système d'équations obtenu à l'étape B-d);
C) Déduire des étapes A et B, la fonction d'atténuation de l'objet i o exprimée en fonction de u.
Lesdits spectres en énergie sont contenus dans une gamme allant de à 500 keV correspondant à la gamme des rayons X ou lesdits spectres en énergie sont contenus dans une gamme allant de 0,5 i s MeV à 20 MeV correspondant à la gamme des rayons y.
Avantageusement, lesdits spectres en énergie sont tels que les spectres successifs se recouvrent partiellement.
Pour permettre la détermination de la fonction d'atténuation, le nombre 2o d'épaisseurs de matériau de référence doit être inférieur ou égal au nombre de spectres en énergie; par exemple, N+1, N entier >_ 2, spectres en énergie sont utilisés, et N épaisseurs de matériau de rëférence sont soumises à ces spectres.
Par exemple, N peut être égal à 9.
2s De manière avantageuse, les puissances du développement augmentées de 1 /2 forment une suite géométrique de raison égale à
une racine de trois, la première puissance du dévelopment étant égale à zéro.
3 o Dans la succession des épaisseurs du matériau de référence, l'ordre de l'épaisseur de référence détermine le choix de la racine de trois.
C'est ainsi que si l'épaisseur de référence est choisie comme la troisième épaisseur dans la succession des épaisseurs du matériau de référence, la raison de la suite sera la racine cubique de trois.
3 s De cette manière, zéro et un font partie de la suite. Zéro correspond à
~WO 94/28398 ~ PCT/FR94/00581 s la puissance de u pour une transmission libre, sans matériau de référence dans le trajet du faisceau incident. Un correspond à la puissance de u pour la transmission à travers l'épaisseur de référence.
Pour la gamme 10 - 500 KeV, de préférence, le matériau de référence s est choisi parmi les corps ayant un numéro atomique effectif, Z effectif, compris entre 5 et 26. De cette manière, le matériau de référence est choisi avec un Z effectif moyen situé entre les matériaux légers et les matériaux lourds.
Pour cette gamme d'énergies comprise entre 10 et 500keV, le i o matériau de référence peut être du duralumin.
L'épaisseur de référence est alors comprise entre 1 et 5mm.
Avantageusement, l'épaisseur de référence est sensiblement égale à
4 mm.
i s Si l'on considère la gamme d'énergies comprise entre 0,5 et 20 MeV, le matériau de référence est choisi parmi les corps ayant un numéro atomique effectif, Z effectif, compris entre 13 et 60.
Dans cette gamme, le matériau de référence est avantageusement de l'acier.
2o Pour cette gamme, l'épaisseur de référence peut être comprise entre 0,5 et 2cm; avantageusement, l'épaisseur de référence est prise égale à lcm~
Pour obtenir une évaluation de la fonction d'atténuation en fonction de 2s l'énergie, les étapes supplémentaires suivantes sont nécessaires:
- A l'aide de faisceaux quasi monochromatiques d'énergie et d'intensité connues, mesurer u en fonction de l'énergie, - En déduire les valeurs de la fonction d'atténuation de l'objet pour ces énergies.
La présente invention présente aussi un dispositif pour la mise en oeuvre du procédé comprenant:
- des moyens pour émettre un faisceau de mesure possédant successivement N+1, N entier >_ 2, spectres en énergie différents, WO 94/28398 PCT/FR94/0058~
- N échantillons en un matériau de référence et possédant des épaisseurs différentes, une de ces épaisseurs étant choisie comme épaisseur de référence, ces échantillons étant aptes à être placës à volonté sur le parcours du faisceau, s - au moins un détecteur apte à la mesure de l'intensité du faisceau pour chaque spectre en énergie ainsi qu'à la mesure de l'intensité
du faisceau transmis par les échantillons et l'objet, soumis au faisceau de mesure, - pour chaque détecteur, des moyens de traitement aptes à
i o calculer, à partir de la mesure de l'intensité du faisceau ainsi que de la mesure de l'intensité du faisceau transmis par les is échantillons et l'objet, la fonction de transmission de l'objet et à en déduire la fonction atténuation exprimée en fonction de la transmission de l'épaisseur de référence du matériau de référence.
Selon une variante, le dispositif pour la mise en oeuvre du procédé
comprend:
- des moyens pour émettre un faisceau de mesure à spectre large, - N, N étant un entier >_ 2, échantillons en un matériau de référence 2o et possédant des épaisseurs différentes, une de ces épaisseurs étant choisie comme épaisseur de référence, - au moins un empilement de détecteurs, chaque détecteur étant sensible à une partie du spectre en énergie, l'empilement étant ordonné du détecteur sensible à la partie du spectre d'énergie la 2s plus faible au détecteur sensible à la partis du spectre d'énergie la plus haute, - pour chaque empilement, connectés en sortie des détecteurs de l'empilement considéré, des moyens de traitement aptes à calculer, à partir de la mesure de l'intensité des différents spectres du 3o faisceau de mesure, de la mesure de l'intensité des spectres correspondant au faisceau de mesure transmis par les échantillons et l'objet, la fonction de transmission de l'objet et en déduire la fonction d'atténuation de l'objet exprimée en fonction de la transmission de l'épaisseur de référence du matériel de référence.
~O 94/28398 PCT/FR94/00581 L'invention et ses avantages seront mieux compris à la lecture de la description qui suit, donnée à titre illustratif et non limitatif et faisant référence aux dessins annexés dans lesquels:
La figure 1 représente schématiquement un dispositif pour la mise en s oeuvre du procédé conforme à l'invention;
La figure 2 représente schématiquement une suite de spectres en énergie utilisée pour mesurer des transmissions;
La figure 3 représente schématiquement une coupe longitudinale d'une pièce en forme de marches;
i o La figure 4 représente schématiquement une variante de réalisation d'un dispositif pour la mise en oeuvre du procédé conforme à
l'invention.
En référence à la figure 1, on décrit maintenant un premier dispositif pour i s la mise en oeuvre du procédé de l'invention. Ce dispositif comprend un générateur 10 de rayons X (gamme d'énergie allant de 10 à 500 keV) ou de rayons 'y (gamme d'énergie allant de 0,5 à 20 Mev) selon la gamme d'énergie souhaitée. Dans cet exemple de réalisation, sous l'effet de la commande 12 à potentiel variable, le générateur 10 délivre un faisceau 20 14 dont l'énergie maximum est variable.
Un système 16 permet de faire défiler différents filtres 16a, en synchronisme avec la variation d'énergie du faisceau. Dans l'exemple de la figure 1, le système à filtres comprend un disque dont les 2s portions correspondent chacune à un filtre différent. Les filtres sont du type passe-haut et sont connus en eux mêmes. La rotation du disque est entraînée par un système de commande 17 connectée à la commande 12 pour régler la synchronisation.
3 o D'autre systèmes à filtres peuvent être utilisés de manière équivalente;
par exemple, les filtres, disposés successivement, peuvent défiler dans le parcours du faisceau sous l'effet d'une translation.
Sur la figure 2, on a représenté la suite des spectres obtenus 3 s successivement lors de la rotation du disque en fonction de l'énergie WO 94/28398 PCT/~'It94/0058~
a du faisceau. Les filtres sont choisis de manières à obtenir un chevauchement des spectres successifs.
Dans le spectre en énergie du faisceau 14 délivré par le générateur 10, chaque filtre supprime toutes les énergies inférieures à un seuil s caractéristique du filtre. Pour chaque spectre de la figure 2, la partie des hautes énergies correspond à l'énergie maximale délivrée par le générateur 10 lorsque le filtre associé est placé dans le parcours du faisceau.
i o Dans la pratique, les spectres en énergie ne sont pas caractérisés de manière précise, et leur allure exacte est sans incidence sur la mesure. Par contre, il est important que l'allure et l'intensité de chaque spectre ne varient pas entre le moment où l'on effectue les mesures de référence et les mesures sur l'objet à caractériser.
is Par ailleurs, les faibles discontinuités dans la succession des spectres sont, elles aussi, sans incidence sur les mesures. Toutefois, pour ne pas perturber les mesures, ces discontinuités doivent être inférieures à 5 keV
pour la plage allant de 20 à 40 keV et inférieures à 10 kev au delà.
De retour à la figure 1, on voit que le dispositif comprend une pièce 18 en forme de marches d'escalier. Chaque marche correspond à un échantillon de matériau d'épaisseur différente. Un exemple de pièce 18 est représenté schématiquement en coupe longitudinale sur la 2s figure 3.
Cette pièce 18 est réalisée dans un matériau de référence dont le Z
effectif peut être compris entre 5 et 26 pour un faisceau incident de rayons X. On choisit par exemple le Duralumin (mélange composé de 95% AI, 4,5% Cu, et 0,5%Mn) dont le Z effectif est environ égal à 13,5.
Pour un faisceau incident de rayons y, le matériau de référence est chosi parmi les matériaux ayant un Z effectif compris entre 13 et 60. On choisit par exemple l'acier dont le Z effectif est égal à 26.
3s On choisit la première marche comme étant la marche de moindre ~O 94/28398 s épaisseur, l'épaisseur des marches suivantes va croissant.
Dans cette réalisation, l'épaisseur de la seconde marche 18a est choisie comme épaisseur de référence. Bien entendu, toute autre s marche pourrait être choisie comme marche de référence.
L'épaisseur de référence pour le Duralumin peut ëtre choisie dans la gamme allant de 1 à 5 mm, par exemple 4mm.
L'épaisseur de référence pour l'acier peut être choisie dans la gamme i o allant de 0,5 à 2cm, par exemple 1 cm.
Une épaisseur de référence trop faible défavorise la sensibilité aux hautes énergies; une épaisseur de référence trop importante défavorise la sensibilité aux basses énergies.
ls D'une manière générale, un nombre N d'épaisseurs distinctes de matériau de référence et donc de marches de la pièce 18, est nécessaire lorsque N+1 spectres sont utilisés. La pièce 18 représentée sur la figure 3 comprend neuf marches et est donc 2o adaptée à un dispositif utilisant dix spectres en énergie.
A nouveau en référence à la figure 1, on voit que la pièce 18 peut ëtre placöae à volonté sur le parcours du faisceau grâce à un mouvement de translation. Pour des raisons de clareté de représentation, le 2s dispasitif permetfiant la translation de la pièce 18 n'est pas représenté
sur la figure 1.
Par ailleurs, la pièce 18 peut prendre d'autres formes, équivalentes à
celle de la pièce représentée. Par exemple, la pièce 18 peut prendre 30 la forme d'un disque dont les portions possèdent les différentes épaisseurs souhaitées et dont une portion est évidée pour laisser libre le trajet du faisceau.
Sur la figure 1, on voit que le dispositif comprend un détecteur 20 3s délivrant un signal électrique proportionnel à l'intensité du faisceau WO 94/28398 PCT/FR94/00581~
atteignant le détecteur 20.
Le détecteur 20 est relié à un système de commande et de traitement 22. Le système 22 peut être un ordinateur, par exemple. Le système s 22 gère les variations de la commande de potentiel variable 12.
On décrit maintenant le procédé conforme à l'invention, pour déterminer la fonction d'atténuation d'un objet.
Lors d'une étape préliminaire, on effectue des mesures de calibration.
i o Ces mesures de calibration consistent à mesurer l'intensité transmise par les différentes épaisseurs du matériau de référence soumises au faisceau de mesure.
Dans un premier temps, on mesure cette intensité transmise pour une épaisseur nulle, c'est à dire sans que la pièce 18 soit interposée dans i s le trajet du faisceau.
Le système 22 commande la variation du potentiel variable et le filtrage adéquat du faisceau par le système 16 à filtres. Sous l'effet conjugué de la variation de potentiel et du filtrage adéquat, le faisceau de mesure 14 présente successivement N+1 spectres en énergie 2o différents. N+1 peut être par exempte égal à 10.
A chaque valeur d'énergie maximum du faisceau correspond un filtre du système 16 pour obtenir successivement les différents spectres en énergie. Chaque variation de l'énergie maximum est synchronisée 2s avec le positionnement d'un nouveau filtre.
Le détecteur 20 mesure l'intensité intégrée du faisceau de mesure pour chacun des spectres et ces mesures sont mises en mémoire par le système de traitement 22.
3 o On refait ensuite ces mêmes mesures mais en plaçant successivement sur le trajet du faisceau les différentes épaisseurs du matériau de référence grâce à une translation de la pièce 18. Une fois ces mesures mises en mémoire par le système 22, la pièce 18 est retirée du trajet du faisceau.
3 s La calibration du dispositif est donc très rapide. Lors d'une utilisation ~O 94/28398 PCT/FR94100581 i~
continue du dispositif, elle peut être renouvelée plusieurs fois par jour, car si l'allure exacte de chacun des spectres est sans importance pour les précisions de la mesure (puisque les mesures sont intégrées en énergie), il est par contre important que d'une mesure à l'autre, l'allure s et l'intensité des spectres soient identiques.
L'objet 24 à caractériser est ensuite placé sur le trajet du faisceau.
L'abjet 24 est soumis successivement à chacun des spectres en énergie réalisés par la variation de l'énergie maximum de faisceau i o synchronisée avec l'interposition du filtre adéquat. Pour chaque spectre en énergie, l'intensité transmise après traversée de l'objet 24 est mesurée par le détecteur 20 puis mise en mémoire.
Grâce à ces mesures, le système 22 déduit une formule analytique de i s l'atténuation de l'objet au point traversé par le faisceau, cette atténuation étant exprimée en fonction d'un paramètre u qui est la transmission d'une épaisseur donnée d'un matériau de référence.
Seul le cas des rayons X, avec un matériau de référence en 2o Duralumin sera envisagé dans la suite de la description. On comprend que la transposition au cas des rayons y et de l'acier est immédiat.
Pour épaisseur de référence, on choisit l'épaisseur de la seconde marche de la pièce 18.
u est donc égale à:
2 s u - e - att(E) . 2,7 . 0,4 où att(E) est l'atténuation massique;
2,7 correspond à la densité du Duralumin exprimée en g/cm3~
0,4 correspond à l'épaisseur de référence choisie égale à 4mm et exprimée en cm.
Conformément à l'invention, le système 22 met en oeuvre le procédé
suivant.
Si on définit l'indice j comme étant le numéro du spectre auquel est soumis l'objet, j=1 correspondant au spectre d'énergie la plus basse 3 s (le spectre le plus à gauche sur la figure 2), alors l'intensité Dj WO 94/28398 PCT/FR94/00581~
transmise par l'objet s'exprime comme suit:
Dj = j Ij (E) Tr (E) dE (0) ou Dj = J Ij (E) e-att {E) dE (1 ) où Ij (E) est l'intensité du jème spectre en énergie, fonction de l'énergie E, Tr (E) et att (E) sont respectivement la fonction de transmission et d'atténuation exprimées ici en fonction de E.
i o En effectuant le changement de variable, faisant apparaître le paramètre u, associé à l'épaisseur de référence ep (en cm) du matériau choisi de densité p (en g/cm3) u = e -att(E) . ep . p i s Dj s'écrit Dj = j l'j (u) e-att' (u) du (2) où l'j (u) est l'intensité du jème spectre exprimée en fonction de u et att' (u) est la fonction d'atténuation recherchée.
La fonction de transmission de l'objet est exprimée en fonction du paramètre u . Dans une première étape, le système 22 réalise ce changement de variable puis il fait l'approximation suivante: la fonction de transmission de l'objet est exprimée comme un 2s développement polynomial fini en puissances du paramètre u; soit:
N
e-att' (u) ~ ~ a i u f(i) (3) i=0 3o I set un indice qui varie de 0 à N.
Cette approximation s'est revélée en pratique très précise puisque l'erreur introduite n'excède pas 10-4 en valeur relative.
3 s Le nombre de termes du développement {N+1 ) est inférieur ou égal ~O 94/28398 PCT/FR94/00581 ~~~~i au nombre de spectres en énergies différentes utilisé pour les mesures de transmission, il est par conséquent inférieur ou égal au nombre N d'épaisseurs de matériau de référence utilisées pour la calibration plus un. De ce fait, chaque indice i correspond à un numéro s de marche de la pièce 18, l'indice i = 0 correspondant à une épaisseur nulle c'est à dire que la pièce n'est pas dans le trajet du faisceau.
On comprend donc que, pour une même gamme d'énergie, par exemple la gamme allant de 10 keV à 150 keV pour des rayons X, i o plus le nombre de spectres réalisés sera élevé, meilleure sera l'interpolation. On peut montrer que l'utilisation d'une dizaine de spectres pour une plage d'énergie de 140 keV permet d'obtenir la fonction d'atténuation sur cette plage avec une précision de 10-4.
i 5 Paur chaque valeur de i, f (i), puissance du développement, représente un rapport entre une des épaisseurs du matériau de référence et l'épaisseur de référence. De manière avantageuse, ces rapports sont tels que les puissances f (i) augmentées de 1/2 forment une suite géométrique de raison égale à une racine de trois.
Sous forme d'expression mathématique, ce choix peut s'écrire:
f(i)+1/2=ai (f (o}+1/2);
De plus, on choisit f (o) = o , puisque i = o corresponde au trajet libre du faisceau et donc à une épaisseur nulle:
f(i}=(ai-1)/2.
De cette manière, si la seconde marche est choisie comme épaisseur de référence, on a 3 o f (2) = 1 et par conséquent a = ~I3 Si la seconde marche 18a est choisie comme marche de référence, la raison de la suite est égale à racine carrée de trois.
WO 94/28398 ' PCT/FR94/0058 ~~~~~'~2 On peut montrer que ce choix pour les valeurs des puissances du développement et donc pour les épaisseurs des marches permet d'obtenir un développement stable.
s En réalité, à cause des inévitables imprécisions d'usinage, les f (i) ne prennent pas exactement les valeurs calculées. Mais ces faibles écarts sont sans conséquence sur la validité du résultat.
On montre que l'intensité transmise par l'objet pour le spectre j peut io N
se décomposer comme suit : Dj = ~ ai Cji {4) i=0 somme dans laquelle les termes Cji correspondent pour chaque spectre j à l'intensité transmise par les différentes épaisseurs de i s matériau de référence, i = 0 correspondant à l'épaisseur nulle (la pièce 18 étant hors du trajet du faisceau).
Dans l'égalité (4) , les Dj et les Cji sont connus par leur mesure, seuls les coefficients ai sont inconnus. Le nombre d'épaisseurs de matériau 2o de référence (en comptant en plus l'épaisseur nulle) peut être choisi inférieur ou égal au nombre de spectres pour que l'on puisse déterminer les coefficients ai. Dans l'exemple décrit ici, on réalise dix spectres dififérents en énergie, la pièce 18 comporte neuf marches auxquelles il faut adjoindre les mesures supplémentaires effctuées 2s lorsque la pièce 18 est hors du trajet du faisceau.
Par conséquent, lors d'une étape suivante, à l'aide des valeurs mises en mémoire lors de la calibration et des mesures de transmission de l'objet, le système 22 détermine les ai . Ceci peut être réalisé par toute 3 o méthode connue par exemple par la méthode des moindres carrés.
L'atténuation Att est définie comme:
Att = -LogTr où Tr est la transmission.
ss Une fois, les coefficients ai déterminés, le système 22 les remplace par O 94/28398 ~ PCT/FR94/00581 leurs valeurs dans l'équation (3), en prend le logarithme népérien et multiplie par (-1 ). De cette manière, le système 22 déduit la fonction d'atténuation de l'objet 24 en fonction du paramètre u.
s La figure 4 représente schématiquement une variante de réalisation d'un dispositif pour la mise en oeuvre du procédé conforme à
l'invention.
Dans cette variante, au lieu de tester successivement les différents spectres en énergie du faisceau de mesure, un faisceau de mesure à
i o spectre fixe et large est utilisé. Les mesures d'intensité transmises par la pièce 18 ou par l'objet 24 sont effectuées grâce à un empilement 30 de détecteurs 30a, 30b, ... , et jouant chacun le rôle de filtre passe haut paur les détecteurs suivants. C'est ainsi que chaque détecteur délivre un signal électrique correspondant à l'intensité pour une partie i s déterminée du spectre large en énergie.
On comprend donc que dans le dispositif de la figure 1, les objets à
caractériser sont soumis à un faisceau de mesure présentant successivement des spectres en énergie différents de manière à
balayer la gamme en énergie désirés, alors que dans le dispositif de la figure 4, les objets sont soumis à un faisceau de mesure à spectre large correspondant à la gamme en énergie désirée, ce spectre large étant ensuite filtré par les différents détecteurs.
Le dispositif de la figure 4 comprend un générateur 26 de faisceau de 2s mesure qui peut être un générateur de rayons X ou y. Ce générateur 26 est connecié à une commande 28 à potentiel fixe telle que le faisceau de mesure présente un spectre large en énergie. Par spectre large en énergie, on entend toute la gamme d'énergie que l'on désire tester, par exemple la gamme s'étendant de 10 keV à 500 keV pour les rayons X ou 0,5 MeV à 20 MeV pour les rayons y Comme précédemment, le dispositif comprend une pièce 18 en forme de marches d'épaisseurs déterminées et réalisée en un matériau de référence. La pièce 18 peut défiler à volonté dans le faisceau ou bien s s ëtre positionnée hors du trajet du faisceau grâce à un mouvement de WO 94/28398 PCT/FR94/0058~
translation non représenté.
Chaque détecteur 30a, 30b, ... de l'empilement 30 est connecté au système de commande et de traitement 22. Dans cette variante, les mesures d'intensités transmises pour chaque spectre découpé dans le s spectre large par le filtrage dû aux détecteurs sont éffectuées simultanément. Le traitement réalisé par le système 22 pour obtenir la fonction d'atténuation est similaire à celui décrit précédement.
Pour certaines applications, il peut être souhaitable d'obtenir la i o fonction d'atténuation non pas en fonction du paramètre u mais en fonction de l'énergie. Pour cela, il suffit lors d'une étape préliminaire effectuée une fois pour toute, de mesurer la fonction de transmission de l'épaisseur de référence du matériau de référence en fonction de l'énergie. Cette mesure peut être effectuée par tout moyen connu.
i s En effet, cette mesure étant effectuée une fois pour toute, les contraintes de temps ou de difficultés de mise en oeuvre sont acceptables pour obtenir le plus grand nombre de points de mesure possible avec la meilleure précision possible.
2o En mesurant la transmission pour un grand nombre d'énergies différentes (par exemple, 10 énergies différentes comprises dans une gamme allant de 10 keV à 500 keV), la fonction de transmission de l'épaisseur de référence du matériau de référence peut être interpolée avec une précision de 10 -4.
2s Lorsque la mesure de la fonction d'atténuation d'un objet est souhaitée en fonction de l'énergie, il faut donc enregistrer, dans une mémoire du système 22, la fonction de transmission u exprimée en fonction de l'énergie. Une fois la fonction d'atténuation déterminée en fonction du paramètre u, il suffit alors de remplacer u par son 3o expression en fonction de l'énergie pour obtenir la fonction d'atténuation recherchée.
La description précédente se rapporte à la mesure d'une fonction d'atténuation d'un objet en un point de celui-ci. Pour obtenir la fonction ss d'atténuation en chaque point de l'objet, on peut par exemple utiliser ~O 94/2839.8 un faisceau en éventail et une batterie de détecteurs (autant de détecteurs que de points souhaités) disposés le long d'une ligne perpendiculaire à la direction de propagation médiane du faisceau.
Ce dispositif permet d'obtenir les fonctions d'atténuation selon une s ligne, la translation de l'objet permet de faire des mesures selon des lignes successives.
e) en résolvant le système d'équations obtenu à l'étape B-d);
C) Déduire des étapes A et B, la fonction d'atténuation de l'objet i o exprimée en fonction de u.
Lesdits spectres en énergie sont contenus dans une gamme allant de à 500 keV correspondant à la gamme des rayons X ou lesdits spectres en énergie sont contenus dans une gamme allant de 0,5 i s MeV à 20 MeV correspondant à la gamme des rayons y.
Avantageusement, lesdits spectres en énergie sont tels que les spectres successifs se recouvrent partiellement.
Pour permettre la détermination de la fonction d'atténuation, le nombre 2o d'épaisseurs de matériau de référence doit être inférieur ou égal au nombre de spectres en énergie; par exemple, N+1, N entier >_ 2, spectres en énergie sont utilisés, et N épaisseurs de matériau de rëférence sont soumises à ces spectres.
Par exemple, N peut être égal à 9.
2s De manière avantageuse, les puissances du développement augmentées de 1 /2 forment une suite géométrique de raison égale à
une racine de trois, la première puissance du dévelopment étant égale à zéro.
3 o Dans la succession des épaisseurs du matériau de référence, l'ordre de l'épaisseur de référence détermine le choix de la racine de trois.
C'est ainsi que si l'épaisseur de référence est choisie comme la troisième épaisseur dans la succession des épaisseurs du matériau de référence, la raison de la suite sera la racine cubique de trois.
3 s De cette manière, zéro et un font partie de la suite. Zéro correspond à
~WO 94/28398 ~ PCT/FR94/00581 s la puissance de u pour une transmission libre, sans matériau de référence dans le trajet du faisceau incident. Un correspond à la puissance de u pour la transmission à travers l'épaisseur de référence.
Pour la gamme 10 - 500 KeV, de préférence, le matériau de référence s est choisi parmi les corps ayant un numéro atomique effectif, Z effectif, compris entre 5 et 26. De cette manière, le matériau de référence est choisi avec un Z effectif moyen situé entre les matériaux légers et les matériaux lourds.
Pour cette gamme d'énergies comprise entre 10 et 500keV, le i o matériau de référence peut être du duralumin.
L'épaisseur de référence est alors comprise entre 1 et 5mm.
Avantageusement, l'épaisseur de référence est sensiblement égale à
4 mm.
i s Si l'on considère la gamme d'énergies comprise entre 0,5 et 20 MeV, le matériau de référence est choisi parmi les corps ayant un numéro atomique effectif, Z effectif, compris entre 13 et 60.
Dans cette gamme, le matériau de référence est avantageusement de l'acier.
2o Pour cette gamme, l'épaisseur de référence peut être comprise entre 0,5 et 2cm; avantageusement, l'épaisseur de référence est prise égale à lcm~
Pour obtenir une évaluation de la fonction d'atténuation en fonction de 2s l'énergie, les étapes supplémentaires suivantes sont nécessaires:
- A l'aide de faisceaux quasi monochromatiques d'énergie et d'intensité connues, mesurer u en fonction de l'énergie, - En déduire les valeurs de la fonction d'atténuation de l'objet pour ces énergies.
La présente invention présente aussi un dispositif pour la mise en oeuvre du procédé comprenant:
- des moyens pour émettre un faisceau de mesure possédant successivement N+1, N entier >_ 2, spectres en énergie différents, WO 94/28398 PCT/FR94/0058~
- N échantillons en un matériau de référence et possédant des épaisseurs différentes, une de ces épaisseurs étant choisie comme épaisseur de référence, ces échantillons étant aptes à être placës à volonté sur le parcours du faisceau, s - au moins un détecteur apte à la mesure de l'intensité du faisceau pour chaque spectre en énergie ainsi qu'à la mesure de l'intensité
du faisceau transmis par les échantillons et l'objet, soumis au faisceau de mesure, - pour chaque détecteur, des moyens de traitement aptes à
i o calculer, à partir de la mesure de l'intensité du faisceau ainsi que de la mesure de l'intensité du faisceau transmis par les is échantillons et l'objet, la fonction de transmission de l'objet et à en déduire la fonction atténuation exprimée en fonction de la transmission de l'épaisseur de référence du matériau de référence.
Selon une variante, le dispositif pour la mise en oeuvre du procédé
comprend:
- des moyens pour émettre un faisceau de mesure à spectre large, - N, N étant un entier >_ 2, échantillons en un matériau de référence 2o et possédant des épaisseurs différentes, une de ces épaisseurs étant choisie comme épaisseur de référence, - au moins un empilement de détecteurs, chaque détecteur étant sensible à une partie du spectre en énergie, l'empilement étant ordonné du détecteur sensible à la partie du spectre d'énergie la 2s plus faible au détecteur sensible à la partis du spectre d'énergie la plus haute, - pour chaque empilement, connectés en sortie des détecteurs de l'empilement considéré, des moyens de traitement aptes à calculer, à partir de la mesure de l'intensité des différents spectres du 3o faisceau de mesure, de la mesure de l'intensité des spectres correspondant au faisceau de mesure transmis par les échantillons et l'objet, la fonction de transmission de l'objet et en déduire la fonction d'atténuation de l'objet exprimée en fonction de la transmission de l'épaisseur de référence du matériel de référence.
~O 94/28398 PCT/FR94/00581 L'invention et ses avantages seront mieux compris à la lecture de la description qui suit, donnée à titre illustratif et non limitatif et faisant référence aux dessins annexés dans lesquels:
La figure 1 représente schématiquement un dispositif pour la mise en s oeuvre du procédé conforme à l'invention;
La figure 2 représente schématiquement une suite de spectres en énergie utilisée pour mesurer des transmissions;
La figure 3 représente schématiquement une coupe longitudinale d'une pièce en forme de marches;
i o La figure 4 représente schématiquement une variante de réalisation d'un dispositif pour la mise en oeuvre du procédé conforme à
l'invention.
En référence à la figure 1, on décrit maintenant un premier dispositif pour i s la mise en oeuvre du procédé de l'invention. Ce dispositif comprend un générateur 10 de rayons X (gamme d'énergie allant de 10 à 500 keV) ou de rayons 'y (gamme d'énergie allant de 0,5 à 20 Mev) selon la gamme d'énergie souhaitée. Dans cet exemple de réalisation, sous l'effet de la commande 12 à potentiel variable, le générateur 10 délivre un faisceau 20 14 dont l'énergie maximum est variable.
Un système 16 permet de faire défiler différents filtres 16a, en synchronisme avec la variation d'énergie du faisceau. Dans l'exemple de la figure 1, le système à filtres comprend un disque dont les 2s portions correspondent chacune à un filtre différent. Les filtres sont du type passe-haut et sont connus en eux mêmes. La rotation du disque est entraînée par un système de commande 17 connectée à la commande 12 pour régler la synchronisation.
3 o D'autre systèmes à filtres peuvent être utilisés de manière équivalente;
par exemple, les filtres, disposés successivement, peuvent défiler dans le parcours du faisceau sous l'effet d'une translation.
Sur la figure 2, on a représenté la suite des spectres obtenus 3 s successivement lors de la rotation du disque en fonction de l'énergie WO 94/28398 PCT/~'It94/0058~
a du faisceau. Les filtres sont choisis de manières à obtenir un chevauchement des spectres successifs.
Dans le spectre en énergie du faisceau 14 délivré par le générateur 10, chaque filtre supprime toutes les énergies inférieures à un seuil s caractéristique du filtre. Pour chaque spectre de la figure 2, la partie des hautes énergies correspond à l'énergie maximale délivrée par le générateur 10 lorsque le filtre associé est placé dans le parcours du faisceau.
i o Dans la pratique, les spectres en énergie ne sont pas caractérisés de manière précise, et leur allure exacte est sans incidence sur la mesure. Par contre, il est important que l'allure et l'intensité de chaque spectre ne varient pas entre le moment où l'on effectue les mesures de référence et les mesures sur l'objet à caractériser.
is Par ailleurs, les faibles discontinuités dans la succession des spectres sont, elles aussi, sans incidence sur les mesures. Toutefois, pour ne pas perturber les mesures, ces discontinuités doivent être inférieures à 5 keV
pour la plage allant de 20 à 40 keV et inférieures à 10 kev au delà.
De retour à la figure 1, on voit que le dispositif comprend une pièce 18 en forme de marches d'escalier. Chaque marche correspond à un échantillon de matériau d'épaisseur différente. Un exemple de pièce 18 est représenté schématiquement en coupe longitudinale sur la 2s figure 3.
Cette pièce 18 est réalisée dans un matériau de référence dont le Z
effectif peut être compris entre 5 et 26 pour un faisceau incident de rayons X. On choisit par exemple le Duralumin (mélange composé de 95% AI, 4,5% Cu, et 0,5%Mn) dont le Z effectif est environ égal à 13,5.
Pour un faisceau incident de rayons y, le matériau de référence est chosi parmi les matériaux ayant un Z effectif compris entre 13 et 60. On choisit par exemple l'acier dont le Z effectif est égal à 26.
3s On choisit la première marche comme étant la marche de moindre ~O 94/28398 s épaisseur, l'épaisseur des marches suivantes va croissant.
Dans cette réalisation, l'épaisseur de la seconde marche 18a est choisie comme épaisseur de référence. Bien entendu, toute autre s marche pourrait être choisie comme marche de référence.
L'épaisseur de référence pour le Duralumin peut ëtre choisie dans la gamme allant de 1 à 5 mm, par exemple 4mm.
L'épaisseur de référence pour l'acier peut être choisie dans la gamme i o allant de 0,5 à 2cm, par exemple 1 cm.
Une épaisseur de référence trop faible défavorise la sensibilité aux hautes énergies; une épaisseur de référence trop importante défavorise la sensibilité aux basses énergies.
ls D'une manière générale, un nombre N d'épaisseurs distinctes de matériau de référence et donc de marches de la pièce 18, est nécessaire lorsque N+1 spectres sont utilisés. La pièce 18 représentée sur la figure 3 comprend neuf marches et est donc 2o adaptée à un dispositif utilisant dix spectres en énergie.
A nouveau en référence à la figure 1, on voit que la pièce 18 peut ëtre placöae à volonté sur le parcours du faisceau grâce à un mouvement de translation. Pour des raisons de clareté de représentation, le 2s dispasitif permetfiant la translation de la pièce 18 n'est pas représenté
sur la figure 1.
Par ailleurs, la pièce 18 peut prendre d'autres formes, équivalentes à
celle de la pièce représentée. Par exemple, la pièce 18 peut prendre 30 la forme d'un disque dont les portions possèdent les différentes épaisseurs souhaitées et dont une portion est évidée pour laisser libre le trajet du faisceau.
Sur la figure 1, on voit que le dispositif comprend un détecteur 20 3s délivrant un signal électrique proportionnel à l'intensité du faisceau WO 94/28398 PCT/FR94/00581~
atteignant le détecteur 20.
Le détecteur 20 est relié à un système de commande et de traitement 22. Le système 22 peut être un ordinateur, par exemple. Le système s 22 gère les variations de la commande de potentiel variable 12.
On décrit maintenant le procédé conforme à l'invention, pour déterminer la fonction d'atténuation d'un objet.
Lors d'une étape préliminaire, on effectue des mesures de calibration.
i o Ces mesures de calibration consistent à mesurer l'intensité transmise par les différentes épaisseurs du matériau de référence soumises au faisceau de mesure.
Dans un premier temps, on mesure cette intensité transmise pour une épaisseur nulle, c'est à dire sans que la pièce 18 soit interposée dans i s le trajet du faisceau.
Le système 22 commande la variation du potentiel variable et le filtrage adéquat du faisceau par le système 16 à filtres. Sous l'effet conjugué de la variation de potentiel et du filtrage adéquat, le faisceau de mesure 14 présente successivement N+1 spectres en énergie 2o différents. N+1 peut être par exempte égal à 10.
A chaque valeur d'énergie maximum du faisceau correspond un filtre du système 16 pour obtenir successivement les différents spectres en énergie. Chaque variation de l'énergie maximum est synchronisée 2s avec le positionnement d'un nouveau filtre.
Le détecteur 20 mesure l'intensité intégrée du faisceau de mesure pour chacun des spectres et ces mesures sont mises en mémoire par le système de traitement 22.
3 o On refait ensuite ces mêmes mesures mais en plaçant successivement sur le trajet du faisceau les différentes épaisseurs du matériau de référence grâce à une translation de la pièce 18. Une fois ces mesures mises en mémoire par le système 22, la pièce 18 est retirée du trajet du faisceau.
3 s La calibration du dispositif est donc très rapide. Lors d'une utilisation ~O 94/28398 PCT/FR94100581 i~
continue du dispositif, elle peut être renouvelée plusieurs fois par jour, car si l'allure exacte de chacun des spectres est sans importance pour les précisions de la mesure (puisque les mesures sont intégrées en énergie), il est par contre important que d'une mesure à l'autre, l'allure s et l'intensité des spectres soient identiques.
L'objet 24 à caractériser est ensuite placé sur le trajet du faisceau.
L'abjet 24 est soumis successivement à chacun des spectres en énergie réalisés par la variation de l'énergie maximum de faisceau i o synchronisée avec l'interposition du filtre adéquat. Pour chaque spectre en énergie, l'intensité transmise après traversée de l'objet 24 est mesurée par le détecteur 20 puis mise en mémoire.
Grâce à ces mesures, le système 22 déduit une formule analytique de i s l'atténuation de l'objet au point traversé par le faisceau, cette atténuation étant exprimée en fonction d'un paramètre u qui est la transmission d'une épaisseur donnée d'un matériau de référence.
Seul le cas des rayons X, avec un matériau de référence en 2o Duralumin sera envisagé dans la suite de la description. On comprend que la transposition au cas des rayons y et de l'acier est immédiat.
Pour épaisseur de référence, on choisit l'épaisseur de la seconde marche de la pièce 18.
u est donc égale à:
2 s u - e - att(E) . 2,7 . 0,4 où att(E) est l'atténuation massique;
2,7 correspond à la densité du Duralumin exprimée en g/cm3~
0,4 correspond à l'épaisseur de référence choisie égale à 4mm et exprimée en cm.
Conformément à l'invention, le système 22 met en oeuvre le procédé
suivant.
Si on définit l'indice j comme étant le numéro du spectre auquel est soumis l'objet, j=1 correspondant au spectre d'énergie la plus basse 3 s (le spectre le plus à gauche sur la figure 2), alors l'intensité Dj WO 94/28398 PCT/FR94/00581~
transmise par l'objet s'exprime comme suit:
Dj = j Ij (E) Tr (E) dE (0) ou Dj = J Ij (E) e-att {E) dE (1 ) où Ij (E) est l'intensité du jème spectre en énergie, fonction de l'énergie E, Tr (E) et att (E) sont respectivement la fonction de transmission et d'atténuation exprimées ici en fonction de E.
i o En effectuant le changement de variable, faisant apparaître le paramètre u, associé à l'épaisseur de référence ep (en cm) du matériau choisi de densité p (en g/cm3) u = e -att(E) . ep . p i s Dj s'écrit Dj = j l'j (u) e-att' (u) du (2) où l'j (u) est l'intensité du jème spectre exprimée en fonction de u et att' (u) est la fonction d'atténuation recherchée.
La fonction de transmission de l'objet est exprimée en fonction du paramètre u . Dans une première étape, le système 22 réalise ce changement de variable puis il fait l'approximation suivante: la fonction de transmission de l'objet est exprimée comme un 2s développement polynomial fini en puissances du paramètre u; soit:
N
e-att' (u) ~ ~ a i u f(i) (3) i=0 3o I set un indice qui varie de 0 à N.
Cette approximation s'est revélée en pratique très précise puisque l'erreur introduite n'excède pas 10-4 en valeur relative.
3 s Le nombre de termes du développement {N+1 ) est inférieur ou égal ~O 94/28398 PCT/FR94/00581 ~~~~i au nombre de spectres en énergies différentes utilisé pour les mesures de transmission, il est par conséquent inférieur ou égal au nombre N d'épaisseurs de matériau de référence utilisées pour la calibration plus un. De ce fait, chaque indice i correspond à un numéro s de marche de la pièce 18, l'indice i = 0 correspondant à une épaisseur nulle c'est à dire que la pièce n'est pas dans le trajet du faisceau.
On comprend donc que, pour une même gamme d'énergie, par exemple la gamme allant de 10 keV à 150 keV pour des rayons X, i o plus le nombre de spectres réalisés sera élevé, meilleure sera l'interpolation. On peut montrer que l'utilisation d'une dizaine de spectres pour une plage d'énergie de 140 keV permet d'obtenir la fonction d'atténuation sur cette plage avec une précision de 10-4.
i 5 Paur chaque valeur de i, f (i), puissance du développement, représente un rapport entre une des épaisseurs du matériau de référence et l'épaisseur de référence. De manière avantageuse, ces rapports sont tels que les puissances f (i) augmentées de 1/2 forment une suite géométrique de raison égale à une racine de trois.
Sous forme d'expression mathématique, ce choix peut s'écrire:
f(i)+1/2=ai (f (o}+1/2);
De plus, on choisit f (o) = o , puisque i = o corresponde au trajet libre du faisceau et donc à une épaisseur nulle:
f(i}=(ai-1)/2.
De cette manière, si la seconde marche est choisie comme épaisseur de référence, on a 3 o f (2) = 1 et par conséquent a = ~I3 Si la seconde marche 18a est choisie comme marche de référence, la raison de la suite est égale à racine carrée de trois.
WO 94/28398 ' PCT/FR94/0058 ~~~~~'~2 On peut montrer que ce choix pour les valeurs des puissances du développement et donc pour les épaisseurs des marches permet d'obtenir un développement stable.
s En réalité, à cause des inévitables imprécisions d'usinage, les f (i) ne prennent pas exactement les valeurs calculées. Mais ces faibles écarts sont sans conséquence sur la validité du résultat.
On montre que l'intensité transmise par l'objet pour le spectre j peut io N
se décomposer comme suit : Dj = ~ ai Cji {4) i=0 somme dans laquelle les termes Cji correspondent pour chaque spectre j à l'intensité transmise par les différentes épaisseurs de i s matériau de référence, i = 0 correspondant à l'épaisseur nulle (la pièce 18 étant hors du trajet du faisceau).
Dans l'égalité (4) , les Dj et les Cji sont connus par leur mesure, seuls les coefficients ai sont inconnus. Le nombre d'épaisseurs de matériau 2o de référence (en comptant en plus l'épaisseur nulle) peut être choisi inférieur ou égal au nombre de spectres pour que l'on puisse déterminer les coefficients ai. Dans l'exemple décrit ici, on réalise dix spectres dififérents en énergie, la pièce 18 comporte neuf marches auxquelles il faut adjoindre les mesures supplémentaires effctuées 2s lorsque la pièce 18 est hors du trajet du faisceau.
Par conséquent, lors d'une étape suivante, à l'aide des valeurs mises en mémoire lors de la calibration et des mesures de transmission de l'objet, le système 22 détermine les ai . Ceci peut être réalisé par toute 3 o méthode connue par exemple par la méthode des moindres carrés.
L'atténuation Att est définie comme:
Att = -LogTr où Tr est la transmission.
ss Une fois, les coefficients ai déterminés, le système 22 les remplace par O 94/28398 ~ PCT/FR94/00581 leurs valeurs dans l'équation (3), en prend le logarithme népérien et multiplie par (-1 ). De cette manière, le système 22 déduit la fonction d'atténuation de l'objet 24 en fonction du paramètre u.
s La figure 4 représente schématiquement une variante de réalisation d'un dispositif pour la mise en oeuvre du procédé conforme à
l'invention.
Dans cette variante, au lieu de tester successivement les différents spectres en énergie du faisceau de mesure, un faisceau de mesure à
i o spectre fixe et large est utilisé. Les mesures d'intensité transmises par la pièce 18 ou par l'objet 24 sont effectuées grâce à un empilement 30 de détecteurs 30a, 30b, ... , et jouant chacun le rôle de filtre passe haut paur les détecteurs suivants. C'est ainsi que chaque détecteur délivre un signal électrique correspondant à l'intensité pour une partie i s déterminée du spectre large en énergie.
On comprend donc que dans le dispositif de la figure 1, les objets à
caractériser sont soumis à un faisceau de mesure présentant successivement des spectres en énergie différents de manière à
balayer la gamme en énergie désirés, alors que dans le dispositif de la figure 4, les objets sont soumis à un faisceau de mesure à spectre large correspondant à la gamme en énergie désirée, ce spectre large étant ensuite filtré par les différents détecteurs.
Le dispositif de la figure 4 comprend un générateur 26 de faisceau de 2s mesure qui peut être un générateur de rayons X ou y. Ce générateur 26 est connecié à une commande 28 à potentiel fixe telle que le faisceau de mesure présente un spectre large en énergie. Par spectre large en énergie, on entend toute la gamme d'énergie que l'on désire tester, par exemple la gamme s'étendant de 10 keV à 500 keV pour les rayons X ou 0,5 MeV à 20 MeV pour les rayons y Comme précédemment, le dispositif comprend une pièce 18 en forme de marches d'épaisseurs déterminées et réalisée en un matériau de référence. La pièce 18 peut défiler à volonté dans le faisceau ou bien s s ëtre positionnée hors du trajet du faisceau grâce à un mouvement de WO 94/28398 PCT/FR94/0058~
translation non représenté.
Chaque détecteur 30a, 30b, ... de l'empilement 30 est connecté au système de commande et de traitement 22. Dans cette variante, les mesures d'intensités transmises pour chaque spectre découpé dans le s spectre large par le filtrage dû aux détecteurs sont éffectuées simultanément. Le traitement réalisé par le système 22 pour obtenir la fonction d'atténuation est similaire à celui décrit précédement.
Pour certaines applications, il peut être souhaitable d'obtenir la i o fonction d'atténuation non pas en fonction du paramètre u mais en fonction de l'énergie. Pour cela, il suffit lors d'une étape préliminaire effectuée une fois pour toute, de mesurer la fonction de transmission de l'épaisseur de référence du matériau de référence en fonction de l'énergie. Cette mesure peut être effectuée par tout moyen connu.
i s En effet, cette mesure étant effectuée une fois pour toute, les contraintes de temps ou de difficultés de mise en oeuvre sont acceptables pour obtenir le plus grand nombre de points de mesure possible avec la meilleure précision possible.
2o En mesurant la transmission pour un grand nombre d'énergies différentes (par exemple, 10 énergies différentes comprises dans une gamme allant de 10 keV à 500 keV), la fonction de transmission de l'épaisseur de référence du matériau de référence peut être interpolée avec une précision de 10 -4.
2s Lorsque la mesure de la fonction d'atténuation d'un objet est souhaitée en fonction de l'énergie, il faut donc enregistrer, dans une mémoire du système 22, la fonction de transmission u exprimée en fonction de l'énergie. Une fois la fonction d'atténuation déterminée en fonction du paramètre u, il suffit alors de remplacer u par son 3o expression en fonction de l'énergie pour obtenir la fonction d'atténuation recherchée.
La description précédente se rapporte à la mesure d'une fonction d'atténuation d'un objet en un point de celui-ci. Pour obtenir la fonction ss d'atténuation en chaque point de l'objet, on peut par exemple utiliser ~O 94/2839.8 un faisceau en éventail et une batterie de détecteurs (autant de détecteurs que de points souhaités) disposés le long d'une ligne perpendiculaire à la direction de propagation médiane du faisceau.
Ce dispositif permet d'obtenir les fonctions d'atténuation selon une s ligne, la translation de l'objet permet de faire des mesures selon des lignes successives.
Claims (18)
1. Procédé pour déterminer la fonction d'atténuation d'un objet, caractérisé en ce qu'il comprend les étapes suivantes:
A) Exprimer la fonction de transmission de l'objet comme un développement fini en puissances de la transmission u d'une épaisseur de référence d'un matériau de référence, chaque puissance de u étant affectée d'un coefficient, les puissances du développement étant égales à un rapport entre des épaisseurs préalablement déterminées du matériau de référence et l'épaisseur de référence;
B) Déterminer les coefficients du développement:
a) en mesurant l'intensité d'un faisceau de mesure pour différents spectres en énergie, b) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par chacune des épaisseurs considérées de matériau de référence, c) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par l'objet, d) pour chacun desdits spectres en énergie, en décomposant l'intensité du faisceau transmis par l'objet en une somme des intensités mesurées à l'étape B-b) relatives au spectre considéré, ces intensités étant affectées des coefficients recherchés, e) en résolvant le système d'équations obtenu à
l'étape B-d);
C) Déduire les étapes A et B, la fonction d'atténuation de l'objet exprimée en fonction de u.
18a
A) Exprimer la fonction de transmission de l'objet comme un développement fini en puissances de la transmission u d'une épaisseur de référence d'un matériau de référence, chaque puissance de u étant affectée d'un coefficient, les puissances du développement étant égales à un rapport entre des épaisseurs préalablement déterminées du matériau de référence et l'épaisseur de référence;
B) Déterminer les coefficients du développement:
a) en mesurant l'intensité d'un faisceau de mesure pour différents spectres en énergie, b) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par chacune des épaisseurs considérées de matériau de référence, c) pour chacun desdits spectres en énergie, en mesurant l'intensité du faisceau transmise par l'objet, d) pour chacun desdits spectres en énergie, en décomposant l'intensité du faisceau transmis par l'objet en une somme des intensités mesurées à l'étape B-b) relatives au spectre considéré, ces intensités étant affectées des coefficients recherchés, e) en résolvant le système d'équations obtenu à
l'étape B-d);
C) Déduire les étapes A et B, la fonction d'atténuation de l'objet exprimée en fonction de u.
18a
2. Procédé selon la revendication 1 caractérisé en ce que lesdits spectres en énergie sont contenus dans une gamme allant de 10 à 500 KeV.
3. Procédé selon la revendication 1 caractérisé en ce que lesdits spectres en énergie sont contenus dans une gamme allant de 500 KeV à 20 MeV.
4. Procédé selon la revendication 1 caractérisé en ce que lesdits spectres en énergie sont tels que les spectres successifs se recouvrent partiellement.
5. Procédé selon la revendication 1 caractérisé en ce que N+1, N entier >= 2, spectres en énergie sont utilisés, et en ce que N
épaisseurs de matériau de référence sont soumises à ces spectres.
épaisseurs de matériau de référence sont soumises à ces spectres.
6. Procédé selon la revendication 5 caractérisé en ce que N est égal à 9.
7. Procédé selon la revendication 1 caractérisé en ce que les puissances du développement augmentées de 1/2 forment une suite géométrique de raison égale à une racine de trois, et la première puissance du développement étant égale à zéro.
8. Procédé selon la revendication 2 caractérisé en ce que le matériau de référence est choisi parmi les corps ayant un Z
effectif compris entre 5 et 26.
effectif compris entre 5 et 26.
9. Procédé selon la revendication 8 caractérisé en ce que le matériau de référence est du duralumin.
10. Procédé selon la revendication 8 caractérisé en ce que l'épaisseur de référence est comprise entre 1 mm et 5 mm.
11. Procédé selon la revendication 10 caractérisé en ce que l'épaisseur de référence est sensiblement égale à 4 mm.
12. Procédé selon la revendication 3 caractérisé en ce que le matériau de référence est choisi parmi les corps ayant un Z
effectif compris entre 13 et 60.
effectif compris entre 13 et 60.
13. Procédé selon la revendication 12 caractérisé en ce que le matériau de référence est de l'acier.
14. Procédé selon la revendication 12 caractérisé en ce que l'épaisseur de référence est comprise entre 0,5 et 2cm.
15. Procédé selon la revendication 14 caractérisé en ce que l'épaisseur de référence est sensiblement égale à 1 cm.
16. Procédé selon la revendication 1 caractérisé en ce qu'il comporte les étapes supplémentaires suivantes:
à l'aide de faisceaux quasi monochromatiques d'énergie et d'intensité connues, mesurer u en fonction de l'énergie, en déduire les valeurs de la fonction d'atténuation de l'objet pour ces énergies.
à l'aide de faisceaux quasi monochromatiques d'énergie et d'intensité connues, mesurer u en fonction de l'énergie, en déduire les valeurs de la fonction d'atténuation de l'objet pour ces énergies.
17. Dispositif pour la mise en oeuvre du procédé de la revendication 1 caractérisé en ce qu'il comprend:
des moyens (10, 16) pour émettre un faisceau de mesure (14) possédant successivement N+1, N entier >= 2, spectres en énergie différents, N échantillons en un matériau de référence et possédant des épaisseurs différentes, une de ces épaisseurs (18a) étant choisie comme épaisseur de référence, ces échantillons étant aptes à
être placés à volonté sur le parcours du faisceau (14), au moins un détecteur (20) apte à la mesure de l'intensité du faisceau (14) pour chaque spectre en énergie ainsi qu'à la mesure de l'intensité des faisceaux transmis par les échantillons et l'objet (24), soumis aux faisceaux de mesure, pour chaque détecteur (20), des moyens de traitement (22) aptes à calculer, à partir de la mesure de l'intensité du faisceau pour chaque spectre en énergie ainsi que de la mesure de l'intensité
du faisceau transmis par les échantillons et l'objet, la fonction de transmission de l'objet et à en déduire la fonction atténuation exprimée en fonction de la transmission de l'épaisseur de référence du matériau de référence.
des moyens (10, 16) pour émettre un faisceau de mesure (14) possédant successivement N+1, N entier >= 2, spectres en énergie différents, N échantillons en un matériau de référence et possédant des épaisseurs différentes, une de ces épaisseurs (18a) étant choisie comme épaisseur de référence, ces échantillons étant aptes à
être placés à volonté sur le parcours du faisceau (14), au moins un détecteur (20) apte à la mesure de l'intensité du faisceau (14) pour chaque spectre en énergie ainsi qu'à la mesure de l'intensité des faisceaux transmis par les échantillons et l'objet (24), soumis aux faisceaux de mesure, pour chaque détecteur (20), des moyens de traitement (22) aptes à calculer, à partir de la mesure de l'intensité du faisceau pour chaque spectre en énergie ainsi que de la mesure de l'intensité
du faisceau transmis par les échantillons et l'objet, la fonction de transmission de l'objet et à en déduire la fonction atténuation exprimée en fonction de la transmission de l'épaisseur de référence du matériau de référence.
18. Dispositif pour la mise en oeuvre du procédé de la revendication 1 caractérisé en ce qu'il comprend:
des moyens (26) pour émettre un faisceau de mesure à spectre large, N, N étant un entier>=2, échantillons en un matériau de référence et possédant des épaisseurs différentes, une (18a) de ces épaisseurs étant choisie comme épaisseur de référence, au moins un empilement de détecteurs (30), chaque détecteur (30 a, 30 b, ... ) étant sensible à une partie du spectre en énergie, l'empilement étant ordonné du détecteur sensible à la partie du spectre d'énergie la plus faible au détecteur sensible à la partie du spectre d'énergie la plus haute, pour chaque empilement (30), connectés en sortie des détecteurs de l'empilement considéré, des moyens de traitement (22) aptes à calculer, à partir de la mesure de l'intensité des différents spectres du faisceau de mesure, de la mesure de l'intensité des spectres correspondant au faisceau de mesure transmis par les échantillons et l'objet, la fonction de transmission de l'objet et en déduire la fonction d'atténuation de l'objet exprimée en fonction de la transmission de l'épaisseur de référence du matériel de référence.
des moyens (26) pour émettre un faisceau de mesure à spectre large, N, N étant un entier>=2, échantillons en un matériau de référence et possédant des épaisseurs différentes, une (18a) de ces épaisseurs étant choisie comme épaisseur de référence, au moins un empilement de détecteurs (30), chaque détecteur (30 a, 30 b, ... ) étant sensible à une partie du spectre en énergie, l'empilement étant ordonné du détecteur sensible à la partie du spectre d'énergie la plus faible au détecteur sensible à la partie du spectre d'énergie la plus haute, pour chaque empilement (30), connectés en sortie des détecteurs de l'empilement considéré, des moyens de traitement (22) aptes à calculer, à partir de la mesure de l'intensité des différents spectres du faisceau de mesure, de la mesure de l'intensité des spectres correspondant au faisceau de mesure transmis par les échantillons et l'objet, la fonction de transmission de l'objet et en déduire la fonction d'atténuation de l'objet exprimée en fonction de la transmission de l'épaisseur de référence du matériel de référence.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR93/06449 | 1993-05-28 | ||
| FR9306449A FR2705785B1 (fr) | 1993-05-28 | 1993-05-28 | Procédé pour déterminer la fonction d'atténuation d'un objet par rapport à la transmission d'une épaisseur de référence d'un matériau de référence et dispositif pour la mise en Óoeuvre du procédé. |
| PCT/FR1994/000581 WO1994028398A1 (fr) | 1993-05-28 | 1994-05-16 | Procede et dispositif pour determiner la fonction d'attenuation d'un objet par rapport a la transmission d'une epaisseur d'un materiau de reference |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2163882A1 CA2163882A1 (fr) | 1994-12-08 |
| CA2163882C true CA2163882C (fr) | 2006-07-18 |
Family
ID=36694005
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002163882A Expired - Lifetime CA2163882C (fr) | 1993-05-28 | 1994-05-16 | Procede et dispositif pour determiner la fonction d'attenuation d'un objet pour ce qui est de la transmission de l'epaisseur d'un materiau de reference |
Country Status (1)
| Country | Link |
|---|---|
| CA (1) | CA2163882C (fr) |
-
1994
- 1994-05-16 CA CA002163882A patent/CA2163882C/fr not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| CA2163882A1 (fr) | 1994-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0700516B1 (fr) | Procede et dispositif pour determiner la fonction d'attenuation d'un objet par rapport a la transmission d'une epaisseur d'un materiau de reference | |
| EP0701703B1 (fr) | Procede et dispositif pour la reconnaissance de materiaux determines dans la composition d'un objet | |
| EP0711419B1 (fr) | Procede et dispositif pour detecter la presence d'un objet, comportant un materiau donne, non accessible a la vue | |
| Breedt et al. | Twelve years of X-ray and optical variability in the Seyfert galaxy NGC 4051 | |
| Burud et al. | An optical time delay estimate for the double gravitational lens system B1600+ 434 | |
| EP0346181B1 (fr) | Procédé et système d'étalonnage d'un scanner à rayons x en utilisant un ou plusieurs étalons circulaires excentrés | |
| EP2798339B1 (fr) | Procédé d'analyse d'un échantillon de matériau par diffractométrie et diffractomètre associé | |
| EP3106864B1 (fr) | Caracterisation d'un echantillon par representation en base de materiaux | |
| FR2626976A1 (fr) | Procede et dispositif pour l'analyse quantitative et differentielle en profondeur d'echantillons solides en mettant en oeuvre deux faisceaux d'ions | |
| EP1852717A1 (fr) | Procédé d'estimation de rayonnement diffusé dans un détecteur bidimensionnel | |
| FR2707115A1 (fr) | Procédé pour produire des formes d'ondelettes sismiques ayant une énergie de lobes latéraux minimale. | |
| EP0762145B1 (fr) | Système de traitement d'impulsions provenant de l'interaction d'une particule gamma avec un détecteur de rayonnement CdTe | |
| EP0485265A1 (fr) | Dispositif et procédé de contrôle non destructif à acquisition simultanée de données radiographiques et de données tomographiques | |
| FR2656109A1 (fr) | Procede de calibration de la chaine de mesure d'un appareil a rayons x. | |
| EP0334762B1 (fr) | Procédé et système d'étalonnage d'un scanner à rayons X en utilisant un seul étalon non circulaire | |
| CA2163882C (fr) | Procede et dispositif pour determiner la fonction d'attenuation d'un objet pour ce qui est de la transmission de l'epaisseur d'un materiau de reference | |
| EP2439520A1 (fr) | Procédé de mesure de l'orientation et de la déformation élastique de grains dans des matériaux multicristallins | |
| Qiu et al. | Proper motion measurements for stars up to 100 kpc with Subaru HSC and SDSS Stripe 82 | |
| Christensen et al. | The host galaxy of GRB 990712 | |
| EP0415815B1 (fr) | Procédé de correction des signaux d'un détecteur linéaire de radiations et dispositif de correction mettant en oeuvre ce procédé | |
| EP1081505A1 (fr) | Dispositif numérique de corrélation et procédé de corrélation pour systèmes télémétriques | |
| CA2163884C (fr) | Procede et dispositif pour la reconnaissance de materiaux determines dans la composition d'un objet | |
| FR2476474A1 (fr) | Dispositif de tomodensitometrie et procede d'exploration et de reconstruction d'images utilisant un tel dispositif | |
| EP1058128B1 (fr) | Dispositif et procédé de discrimination d'impulsions de détecteurs de rayonnement à semi-conducteur | |
| FR2700259A1 (fr) | Dispositif automatique de collimation d'un faisceau de rayons X d'un scanner et procédé de mise en Óoeuvre. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKEX | Expiry |
Effective date: 20140516 |