WO1994028398A1 - 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 - Google Patents

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 Download PDF

Info

Publication number
WO1994028398A1
WO1994028398A1 PCT/FR1994/000581 FR9400581W WO9428398A1 WO 1994028398 A1 WO1994028398 A1 WO 1994028398A1 FR 9400581 W FR9400581 W FR 9400581W WO 9428398 A1 WO9428398 A1 WO 9428398A1
Authority
WO
WIPO (PCT)
Prior art keywords
energy
reference material
intensity
function
spectra
Prior art date
Application number
PCT/FR1994/000581
Other languages
English (en)
Inventor
Serge Maitrejean
Didier Perion
Original Assignee
Europ Scan, S.A.
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 Europ Scan, S.A. filed Critical Europ Scan, S.A.
Priority to AU68003/94A priority Critical patent/AU6800394A/en
Priority to US08/553,427 priority patent/US5687210A/en
Priority to CA002163882A priority patent/CA2163882C/fr
Priority to EP94916291A priority patent/EP0700516B1/fr
Priority to DE69401433T priority patent/DE69401433T2/de
Priority to JP7500297A priority patent/JPH09500719A/ja
Publication of WO1994028398A1 publication Critical patent/WO1994028398A1/fr
Priority to NO954819A priority patent/NO954819L/no
Priority to FI955716A priority patent/FI955716A/fi

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

Definitions

  • the present invention relates to a method for determining the attenuation function of an object and to a device for implementing the method. It applies in particular but not only to radioscopy. In general, it applies to the domain of 1 o spectroscopy in the energy range from a few keV to a few MeV.
  • the object to be characterized is placed on the beam path and a detector measures, for the different energies, the intensity of the beam
  • Att (E) is equal to - Log T r (E).
  • An image has about 500 dots per line and about 600 lines; if each measurement requires a second, the realization of an image, for a single energy, requires approximately 83h, to multiply 15 by the desired number of energies.
  • the Bragg diffraction and the Compton effect are dependent on the angle of incidence of the beam and only allow a beam to be obtained at the desired energy for a given angle of incidence 0 obtained with a parallel beam.
  • this type of collimated beam it is only possible to measure the transmission at one point of the object. To obtain a large number of measurement points simultaneously, it is necessary to multiply the beam emission means, because given the angular dependence, the use of a fan beam is prohibited. It is understood that the implementation of such an arrangement is very complex and costly.
  • isotopic sources are simpler but these sources also have drawbacks. In fact, for safety reasons, the handling of sources is strictly regulated and the intensity of the beams is low.
  • known attenuation measurement techniques allow assessments only for particular energies strongly dependent on the device used and not on the will or needs of the user.
  • they require a very long measurement time, are expensive and difficult to implement; consequently, these techniques are difficult to use in radioscopy, or in security imaging which are important fields of application of the determination of the attenuation of an object.
  • the present invention overcomes these drawbacks. It allows the attenuation function to be determined very quickly (in a time of the order of ten milliseconds) using standard X-ray tubes or ⁇ -ray generators.
  • the attenuation function is determined according to a variable in bijective relation with energy. As will be seen below, the transmission u of a reference thickness of a reference material is chosen as variable. If one wishes to express the attenuation of the object as a function of energy, it suffices to measure once and for all (and not at each determination) the transmission u for the desired energies, this using techniques of the prior art.
  • the invention relates to a method for determining the attenuation function of an object, comprising the following steps:
  • A) Express the transmission function of the object as a finite development in powers of the transmission u of a reference thickness of a reference material, each power of u being assigned a coefficient, the powers of the development being equal to a ratio between previously determined thicknesses of the reference material and the reference thickness;
  • B) Determine the development coefficients: a) by measuring the intensity of a measurement beam for different energy spectra, b) for each of said energy spectra, by measuring the intensity of the beam transmitted by each of the thicknesses considered of reference material, c) for each of said energy spectra, by measuring the intensity of the beam transmitted by the object, d) for each of said spectra, by decomposing the intensity of the beam transmitted by the object into a sum of the intensities measured at step Bb) relating to the spectrum considered, these intensities being affected by the desired coefficients; e) by solving the system of equations obtained in step Bd);
  • Said energy spectra are contained in a range from 10 to 500 keV corresponding to the range of X-rays or said energy spectra are contained in a range ranging from 0.5 MeV to 20 MeV corresponding to the range of ⁇ rays.
  • said energy spectra are such that the successive spectra partially overlap.
  • the number of reference material thicknesses must be less than or equal to the number of energy spectra; for example, N + 1, N integer> 2, energy spectra are used, and N thicknesses of reference material are subjected to these spectra.
  • N can be equal to 9.
  • the powers of development increased by 1/2 form a geometric sequence of reason equal to a root of three, the first power of development being equal to zero.
  • the order of the reference thickness determines the choice of the root of three.
  • the reason for the sequence will be the cube root of three. In this way, zero and one are part of the sequence.
  • Zero corresponds to the power of u for free transmission, without reference material in the path of the incident beam.
  • Un corresponds to the power of u for the transmission through the reference thickness.
  • the reference material 5 is chosen from bodies having an effective atomic number, effective Z, between 5 and 26. In this way, the reference material is chosen with an effective Z medium located between light materials and heavy materials.
  • the reference material can be duralumin.
  • the reference thickness is then between 1 and 5mm.
  • the reference thickness is substantially equal to 4 mm.
  • the reference material is chosen from bodies with an effective atomic number, effective Z, between 13 and 60.
  • the reference material is advantageously steel. 20
  • the reference thickness can be between
  • the reference thickness is taken equal to 1 cm.
  • the present invention also presents a device for implementing the method comprising:
  • a measurement beam having successively N + 1, N integer> 2, different energy spectra, - N samples of a reference material and having different thicknesses, one of these thicknesses being chosen as the reference thickness, these samples being able to be placed at will on the path of the beam, - at least one detector suitable for measuring the intensity of the beam for each energy spectrum as well as the measurement of the intensity of the beam transmitted by the samples and the object, subjected to the measurement beam,
  • processing means able to calculate, from the measurement of the intensity of the beam as well as from the measurement of the intensity of the beam transmitted by the samples and the object, the transmission function of the object and to deduce the attenuation function expressed as a function of the transmission of the reference thickness of the reference material.
  • the device for implementing the method comprises:
  • N N being an integer> 2 samples of a reference material and having different thicknesses, one of these thicknesses being chosen as the reference thickness,
  • each detector being sensitive to part of the energy spectrum
  • the stack being ordered from the detector sensitive to the weakest part of the energy spectrum to the detector sensitive to the part of the energy spectrum the highest
  • FIG. 1 schematically represents a device for implementing the process in accordance to the invention
  • FIG. 2 schematically represents a series of energy spectra used to measure transmissions
  • FIG. 3 schematically represents a longitudinal section of a part in the form of steps
  • î o Figure 4 schematically shows an alternative embodiment of a device for implementing the method according to the invention.
  • This device comprises an X-ray generator 10 (energy range ranging from 10 to 500 keV) or ⁇ rays (energy range ranging from 0.5 to 20 Mev) depending on the desired energy range.
  • the generator 10 delivers a beam 20 14 whose maximum energy is variable.
  • a system 16 makes it possible to cycle through different filters 16a, in synchronism with the energy variation of the beam.
  • the filter system comprises a disc, the 5 portions of which each correspond to a different filter.
  • the filters are of the high pass type and are known in themselves.
  • the rotation of the disc is driven by a control system 17 connected to the control 12 to adjust the synchronization.
  • filters can be used equivalently; for example, the filters, arranged successively, can scroll in the path of the beam under the effect of a translation.
  • FIG. 2 the sequence of spectra obtained successively during the rotation of the disc as a function of the energy has been represented. of the beam.
  • the filters are chosen so as to obtain an overlap of the successive spectra.
  • each filter suppresses all the energies below a characteristic threshold of the filter.
  • the part of the high energies corresponds to the maximum energy delivered by the generator 10 when the associated filter is placed in the path of the beam.
  • the energy spectra are not precisely characterized, and their exact shape has no effect on the measurement.
  • the low discontinuities in the succession of the spectra are also without effect on the measurements.
  • these discontinuities must be less than 5 keV for the range from 20 to 40 keV and less than 10 kev beyond.
  • the device comprises a part 18 in the form of stair treads. Each step corresponds to a sample of material of different thickness.
  • An example of part 18 is shown diagrammatically in longitudinal section in FIG. 3.
  • This part 18 is made of a reference material whose effective Z can be between 5 and 26 for an incident beam of X-rays.
  • Duralumin is chosen (mixture composed of 95% Al, 4.5% Cu, and 0.5% Mn) whose effective Z is approximately equal to 13.5.
  • the reference material is chosen from among the materials having an effective Z of between 13 and 60.
  • steel is chosen whose effective Z is equal to 26.
  • the first step is chosen as being the least step thickness, the thickness of the following steps increases.
  • the thickness of the second step 18a is chosen as the reference thickness.
  • any other step could be chosen as the reference step.
  • the reference thickness for Duralumin can be chosen from the range from 1 to 5 mm, for example 4 mm.
  • the reference thickness for the steel can be chosen from the range from 0.5 to 2 cm, for example 1 cm.
  • Too thin a reference thickness disadvantages sensitivity to high energies; too large a reference thickness disadvantages the sensitivity to low energies.
  • N a number N of distinct thicknesses of reference material and therefore of steps of the piece 18, is necessary when N + 1 spectra are used.
  • the part 18 represented in FIG. 3 comprises nine steps and is therefore 0 suitable for a device using ten energy spectra.
  • the part 18 can be placed at will on the path of the beam thanks to a translational movement.
  • the device allowing the translation of the part 18 is not shown in FIG. 1.
  • the part 18 can take other forms, equivalent to that of the part shown.
  • the part 18 can take the form of a disk whose portions have the various desired thicknesses and of which a portion is hollowed out to leave the path of the beam free.
  • the device comprises a detector 20 5 delivering an electrical signal proportional to the intensity of the beam reaching detector 20.
  • the detector 20 is connected to a control and processing system 22.
  • the system 22 can be a computer, for example.
  • the system 22 manages the variations of the variable potential command 12.
  • calibration measurements are made. î o These calibration measurements consist in measuring the intensity transmitted by the different thicknesses of the reference material subjected to the measurement beam.
  • this transmitted intensity is measured for zero thickness, that is to say without the part 18 being interposed in
  • the system 22 controls the variation of the variable potential and the adequate filtering of the beam by the system 16 with filters. Under the combined effect of the potential variation and the adequate filtering, the measurement beam 14 successively presents N + 1 energy spectra
  • N + 1 can for example be equal to 10.
  • Each maximum energy value of the beam corresponds to a system filter 16 to successively obtain the different energy spectra.
  • Each variation of the maximum energy is synchronized with the positioning of a new filter.
  • the detector 20 measures the integrated intensity of the measurement beam for each of the spectra and these measurements are stored in memory by the processing system 22.
  • the device calibration is therefore very quick. When in use continuous of the device, it can be renewed several times a day, because if the exact shape of each of the spectra is not important for the precision of the measurement (since the measurements are integrated in energy), it is on the other hand important that one measurement to another, the pace and intensity of the spectra are identical.
  • the object 24 to be characterized is then placed on the path of the beam.
  • the object 24 is subjected successively to each of the energy spectra produced by the variation of the maximum beam energy synchronized with the interposition of the appropriate filter.
  • the intensity transmitted after crossing the object 24 is measured by the detector 20 and then stored.
  • the system 22 deduces an analytical formula of the attenuation of the object at the point crossed by the beam, this attenuation being expressed as a function of a parameter u which is the transmission of a given thickness of a reference material.
  • the thickness of the second step of the piece 18 is chosen.
  • 2.7 corresponds to the density of Duralumin expressed in g / cm 3 - 0.4 corresponds to the chosen reference thickness equal to 4mm and expressed in cm.
  • the system 22 implements the following method.
  • j (u) is the intensity of the j th spectrum expressed as a function of u and att '(u) is the sought attenuation function.
  • the transmission function of the object is expressed as a function of the parameter u.
  • the system 22 performs this change of variable then it makes the following approximation: the transmission function of the object is expressed as a finite polynomial expansion in powers of the parameter u; is:
  • f (i) For each value of i, f (i), power of development, represents a relationship between one of the thicknesses of the reference material and the reference thickness.
  • these ratios are such that the powers f (i) increased by 1/2 form a geometric sequence of reason equal to a root of three.
  • the reason for the sequence is equal to the square root of three.
  • the Dj and the Cji are known by their measurement, only the coefficients ai are unknown.
  • the number of thicknesses of reference material (counting in addition the zero thickness) can be chosen less than or equal to the number of spectra so that we can determine the coefficients aj.
  • the part 18 comprises nine steps to which must be added the additional measurements carried out when the part 18 is out of the beam path.
  • the system 22 determines the aj. This can be achieved by any known method, for example by the method of least squares.
  • Att -LogTr where Tr is the transmission.
  • the system 22 replaces them with their values in equation (3), takes the natural logarithm and multiplies by (-1). In this way, the system 22 deduces the attenuation function of the object 24 according to the parameter u.
  • FIG. 4 schematically represents an alternative embodiment of a device for implementing the method according to the invention.
  • a measurement beam with i o fixed and wide spectrum is used instead of successively testing the different energy spectra of the measurement beam.
  • the intensity measurements transmitted by the part 18 or by the object 24 are carried out thanks to a stack 30 of detectors 30a, 30b, ..., and each playing the role of high pass filter for the following detectors. This is how each detector delivers an electrical signal corresponding to the intensity for a part
  • the objects to be characterized are subjected to a measurement beam having successively different energy spectra so as to scan the range of energy desired, while in the device of 0 in FIG. 4 , the objects are subjected to a broad spectrum measurement beam corresponding to the desired energy range, this broad spectrum then being filtered by the various detectors.
  • the device of FIG. 4 comprises a generator 26 of a measurement beam which can be an X or ⁇ ray generator.
  • This generator 26 is connected to a control 28 with fixed potential such that the measurement beam has a broad energy spectrum.
  • broad energy spectrum is meant the whole range of energy which it is desired to test, for example the range extending from 10 keV to 500 keV for 0 X-rays or 0.5 MeV to 20 MeV for rays ⁇ .
  • the device comprises a part 18 in the form of steps of determined thicknesses and made of a reference material.
  • the piece 18 can scroll at will in the beam or else be positioned outside the path of the beam by means of a movement of translation not shown.
  • Each detector 30a, 30b, ... of the stack 30 is connected to the control and processing system 22.
  • the intensity measurements transmitted for each spectrum cut out from the broad spectrum by the filtering due to the detectors are performed simultaneously.
  • the processing carried out by the system 22 to obtain the attenuation function is similar to that described above.
  • the i o attenuation function may be desirable to obtain not as a function of the parameter u but as a function of energy.
  • This measurement can be carried out by any known means. In fact, this measurement being carried out once and for all, the constraints of time or of difficulties of implementation are acceptable in order to obtain the greatest number of measurement points possible with the best possible accuracy.
  • the transmission function of the reference thickness of the reference material can be interpolated with an accuracy of 10 "4. 5
  • the transmission function u expressed in function of the energy
  • the preceding description relates to the measurement of an attenuation function of an object at a point thereof.
  • an attenuation function of an object at a point thereof.
  • This device makes it possible to obtain the attenuation functions along a line, the translation of the object makes it possible to make measurements along successive lines.

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 sont déterminés grâce aux mesures de l'intensité transmise par les différentes épaisseurs 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

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 1 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
15 é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 γ d'énergie stable et déterminée, provenant de sources isotopiques dont l'énergie est éventuellement diminuée par effet
20 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
25 transmis.
Si l'on pose que lo(E) représente l'intensité du faisceau incident d'énergie E, que lt(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 lt(E)/l0(E). L'atténuation considérée pour l'énergie E,
30 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 35 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. 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 5 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 î 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 15 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é 0 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 5 é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 0 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 35 fortement du dispositif utilisé et non de la volonté ou des besoins de 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 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 à rayons X ou des générateurs de rayons γ 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 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 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 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 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 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 exprimée en fonction de u.
Lesdits spectres en énergie sont contenus dans une gamme allant de 10 à 500 keV correspondant à la gamme des rayons X ou lesdits spectres en énergie sont contenus dans une gamme allant de 0,5 MeV à 20 MeV correspondant à la gamme des rayons γ.
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 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.
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 development étant égale à zéro. 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. De cette manière, zéro et un font partie de la suite. Zéro correspond à 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 5 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 î 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.
15 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. 20 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 à 1cm.
Pour obtenir une évaluation de la fonction d'atténuation en fonction de 5 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. 0
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, - 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, - 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 à calculer, à partir de la mesure de l'intensité du faisceau 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.
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 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 plus faible au détecteur sensible à la partie 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 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. 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 5 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; î 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 15 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 γ (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 5 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.
30 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 35 successivement lors de la rotation du disque en fonction de l'énergie 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 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.
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.
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 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% Al, 4,5% Cu, et 0,5%Mn) dont le Z effectif est environ égal à 13,5.
Pour un faisceau incident de rayons γ, 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. On choisit la première marche comme étant la marche de moindre é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 5 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 î 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.
1 5
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 0 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ée à volonté sur le parcours du faisceau grâce à un mouvement de translation. Pour des raisons de clareté de représentation, le 5 dispositif permettant 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 0 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 5 délivrant un signal électrique proportionnel à l'intensité du faisceau 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 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. î 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
15 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
20 différents. N+1 peut être par exemple é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 5 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.
30 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.
35 La calibration du dispositif est donc très rapide. Lors d'une utilisation 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 et l'intensité des spectres soient identiques.
L'objet 24 à caractériser est ensuite placé sur le trajet du faisceau. L'objet 24 est soumis successivement à chacun des spectres en énergie réalisés par la variation de l'énergie maximum de faisceau 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 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 Duralumin sera envisagé dans la suite de la description. On comprend que la transposition au cas des rayons γ 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 à: 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 (le spectre le plus à gauche sur la figure 2), alors l'intensité Dj transmise par l'objet s'exprime comme suit: Dj = i lj (E) Tr (E) dE (0) ou
Dj = J lj (E) e-att (E) dE (1 )
5 où lj (E) est l'intensité du jènπe 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
15 Dj s'écrit :
Dj = j l'j (u) e-att' (u) du (2)
où l'j (u) est l'intensité du jeme spectre exprimée en fonction de u et att' (u) est la fonction d'atténuation recherchée.
20
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 25 développement polynomial fini en puissances du paramètre u; soit:
N e-att' (u) ~ ∑ a j U f(i) (3) i=0
30 I set un indice qui varie de 0 à N.
Cette approximation s'est révélée en pratique très précise puisque l'erreur introduite n'excède pas 10"4 en valeur relative.
35 Le nombre de termes du développement (N+1 ) est inférieur ou égal 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 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, 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.
Pour 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 = a' (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) = (a' - 1) / 2 .
De cette manière, si la seconde marche est choisie comme épaisseur de référence, on a : f (2) = 1 et par conséquent a = 3
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. 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.
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 N se décomposer comme suit : Dj = ∑ aj Cji (4) i=0 somme dans laquelle les termes Cji correspondent pour chaque spectre j à l'intensité transmise par les différentes épaisseurs de 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 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 aj. Dans l'exemple décrit ici, on réalise dix spectres différents en énergie, la pièce 18 comporte neuf marches auxquelles il faut adjoindre les mesures supplémentaires effctuées 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 aj . Ceci peut être réalisé par toute 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.
Une fois, les coefficients aj déterminés, le système 22 les remplace par 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.
5 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 pour les détecteurs suivants. C'est ainsi que chaque détecteur délivre un signal électrique correspondant à l'intensité pour une partie
15 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ée, alors que dans le dispositif de 0 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 5 mesure qui peut être un générateur de rayons X ou γ. Ce générateur 26 est connecté à 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 0 les rayons X ou 0,5 MeV à 20 MeV pour les rayons γ.
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 5 être positionnée hors du trajet du faisceau grâce à un mouvement de 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 5 spectre large par le filtrage dû aux détecteurs sont effectué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. 15 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.
0 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. 5 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 0 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 35 d'atténuation en chaque point de l'objet, on peut par exemple utiliser 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 ligne, la translation de l'objet permet de faire des mesures selon des lignes successives.

Claims

REVENDICATIONS
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; c) en résolvant le système d'équations obtenu à l'étape B-d); C) Déduire des étapes 1 et B, la fonction d'atténuation de l'objet exprimée en fonction de u.
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 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.
i o
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
15 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 0 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.
5 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. 0
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.
5
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.
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.
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.
PCT/FR1994/000581 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 WO1994028398A1 (fr)

Priority Applications (8)

Application Number Priority Date Filing Date Title
AU68003/94A AU6800394A (en) 1993-05-28 1994-05-16 Method and device for determining the attenuation function of an object with respect to the transmission of a reference material thickness
US08/553,427 US5687210A (en) 1993-05-28 1994-05-16 Method and device for determining the attenuation function of an object with respect to the transmission of a reference material thickness
CA002163882A 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
EP94916291A EP0700516B1 (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
DE69401433T DE69401433T2 (de) 1993-05-28 1994-05-16 Verfahren und gerät zur bestimmung der abschwaechungsfunktion eines objektes in verhältnis zur transmission einer referenzmaterialdicke
JP7500297A JPH09500719A (ja) 1993-05-28 1994-05-16 基準材料の厚みの透過率との関係における物体の減衰関数を決定するための方法及び装置
NO954819A NO954819L (no) 1993-05-28 1995-11-27 Fremgangsmåte og anordning for å bestemme svekningsfunksjonen for en gjenstand med hensyn på transmisjon gjennom en tykkelse av et referansemateriale
FI955716A FI955716A (fi) 1993-05-28 1995-11-27 Menetelmä ja laite kohteen vaimennusfunktion määrittämiseksi vertailumateriaalipaksuuden läpäisyyn verrattuna

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
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é.
FR93/06449 1993-05-28

Publications (1)

Publication Number Publication Date
WO1994028398A1 true WO1994028398A1 (fr) 1994-12-08

Family

ID=9447564

Family Applications (1)

Application Number Title Priority Date Filing Date
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

Country Status (10)

Country Link
US (1) US5687210A (fr)
EP (1) EP0700516B1 (fr)
JP (1) JPH09500719A (fr)
AT (1) ATE147510T1 (fr)
AU (1) AU6800394A (fr)
DE (1) DE69401433T2 (fr)
FI (1) FI955716A (fr)
FR (1) FR2705785B1 (fr)
NO (1) NO954819L (fr)
WO (1) WO1994028398A1 (fr)

Families Citing this family (96)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2708751B1 (fr) * 1993-07-30 1995-10-06 Schlumberger Ind Sa Procédé et dispositif pour détecter la présence d'un objet, comportant un matériau donné, non accessible à la vue.
US8480754B2 (en) 2001-05-25 2013-07-09 Conformis, Inc. Patient-adapted and improved articular implants, designs and related guide tools
US8545569B2 (en) 2001-05-25 2013-10-01 Conformis, Inc. Patient selectable knee arthroplasty devices
US8735773B2 (en) 2007-02-14 2014-05-27 Conformis, Inc. Implant device and method for manufacture
US8617242B2 (en) 2001-05-25 2013-12-31 Conformis, Inc. Implant device and method for manufacture
US9603711B2 (en) 2001-05-25 2017-03-28 Conformis, Inc. Patient-adapted and improved articular implants, designs and related guide tools
US8771365B2 (en) 2009-02-25 2014-07-08 Conformis, Inc. Patient-adapted and improved orthopedic implants, designs, and related tools
US8556983B2 (en) 2001-05-25 2013-10-15 Conformis, Inc. Patient-adapted and improved orthopedic implants, designs and related tools
US8882847B2 (en) 2001-05-25 2014-11-11 Conformis, Inc. Patient selectable knee joint arthroplasty devices
US6064716A (en) * 1997-09-05 2000-05-16 Cyberlogic, Inc. Plain x-ray bone densitometry apparatus and method
US5917877A (en) * 1997-09-05 1999-06-29 Cyberlogic, Inc. Plain x-ray bone densitometry apparatus and method
US7184814B2 (en) 1998-09-14 2007-02-27 The Board Of Trustees Of The Leland Stanford Junior University Assessing the condition of a joint and assessing cartilage loss
US7239908B1 (en) 1998-09-14 2007-07-03 The Board Of Trustees Of The Leland Stanford Junior University Assessing the condition of a joint and devising treatment
ATE439806T1 (de) 1998-09-14 2009-09-15 Univ Leland Stanford Junior Zustandsbestimmung eines gelenks und schadenvorsorge
US6173038B1 (en) 1999-12-01 2001-01-09 Cyberlogic, Inc. Plain x-ray bone densitometry apparatus and method
JP4049312B2 (ja) 2000-08-29 2008-02-20 イメージング セラピューティクス，インコーポレーテッド Ｘ線画像の定量解析法および装置
US6904123B2 (en) 2000-08-29 2005-06-07 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
US7050534B2 (en) * 2000-08-29 2006-05-23 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
US7467892B2 (en) * 2000-08-29 2008-12-23 Imaging Therapeutics, Inc. Calibration devices and methods of use thereof
ATE426357T1 (de) 2000-09-14 2009-04-15 Univ Leland Stanford Junior Beurteilung des zustandes eines gelenkes und planung einer behandlung
US8639009B2 (en) 2000-10-11 2014-01-28 Imatx, Inc. Methods and devices for evaluating and treating a bone condition based on x-ray image analysis
US20070047794A1 (en) * 2000-10-11 2007-03-01 Philipp Lang Methods and devices for analysis of x-ray images
US7660453B2 (en) 2000-10-11 2010-02-09 Imaging Therapeutics, Inc. Methods and devices for analysis of x-ray images
US8000766B2 (en) 2001-05-25 2011-08-16 Imatx, Inc. Methods to diagnose treat and prevent bone loss
JP2005504563A (ja) 2001-05-25 2005-02-17 イメージング セラピューティクス，インコーポレーテッド 関節を再表面化する方法および組成物
US6570955B1 (en) 2002-01-08 2003-05-27 Cyberlogic, Inc. Digital x-ray material testing and bone densitometry apparatus and method
US8275091B2 (en) 2002-07-23 2012-09-25 Rapiscan Systems, Inc. Compact mobile cargo scanning system
US8503605B2 (en) * 2002-07-23 2013-08-06 Rapiscan Systems, Inc. Four sided imaging system and method for detection of contraband
US7783004B2 (en) 2002-07-23 2010-08-24 Rapiscan Systems, Inc. Cargo scanning system
US9958569B2 (en) 2002-07-23 2018-05-01 Rapiscan Systems, Inc. Mobile imaging system and method for detection of contraband
US7369643B2 (en) * 2002-07-23 2008-05-06 Rapiscan Security Products, Inc. Single boom cargo scanning system
US7322745B2 (en) * 2002-07-23 2008-01-29 Rapiscan Security Products, Inc. Single boom cargo scanning system
US7486768B2 (en) 2002-07-23 2009-02-03 Rapiscan Security Products, Inc. Self-contained mobile inspection system and method
US7963695B2 (en) 2002-07-23 2011-06-21 Rapiscan Systems, Inc. Rotatable boom cargo scanning system
US7840247B2 (en) 2002-09-16 2010-11-23 Imatx, Inc. Methods of predicting musculoskeletal disease
US8965075B2 (en) 2002-09-16 2015-02-24 Imatx, Inc. System and method for predicting future fractures
EP1555962B1 (fr) 2002-10-07 2011-02-09 Conformis, Inc. Implant articulaire par chirurgie non effractive a geometrie tridimensionnelle correspondant aux surfaces articulaires
JP2006505366A (ja) 2002-11-07 2006-02-16 コンフォーミス・インコーポレイテッド 半月板サイズおよび形状の決定および工夫した処置の方法
FR2849241B1 (fr) * 2002-12-20 2005-06-24 Biospace Instr Procede et dispositif d'imagerie radiographique
CA2519187A1 (fr) 2003-03-25 2004-10-14 Imaging Therapeutics, Inc. Procedes de compensation de technique d'imagerie dans le traitement d'images radiographiques
GB0525593D0 (en) 2005-12-16 2006-01-25 Cxr Ltd X-ray tomography inspection systems
US9113839B2 (en) 2003-04-25 2015-08-25 Rapiscon Systems, Inc. X-ray inspection system and method
US20050058242A1 (en) 2003-09-15 2005-03-17 Peschmann Kristian R. Methods and systems for the rapid detection of concealed objects
US8243876B2 (en) 2003-04-25 2012-08-14 Rapiscan Systems, Inc. X-ray scanners
US8451974B2 (en) 2003-04-25 2013-05-28 Rapiscan Systems, Inc. X-ray tomographic inspection system for the identification of specific target items
US7949101B2 (en) 2005-12-16 2011-05-24 Rapiscan Systems, Inc. X-ray scanners and X-ray sources therefor
US8837669B2 (en) 2003-04-25 2014-09-16 Rapiscan Systems, Inc. X-ray scanning system
US8223919B2 (en) 2003-04-25 2012-07-17 Rapiscan Systems, Inc. X-ray tomographic inspection systems for the identification of specific target items
US6928141B2 (en) * 2003-06-20 2005-08-09 Rapiscan, Inc. Relocatable X-ray imaging system and method for inspecting commercial vehicles and cargo containers
US7856081B2 (en) 2003-09-15 2010-12-21 Rapiscan Systems, Inc. Methods and systems for rapid detection of concealed objects using fluorescence
US8290564B2 (en) 2003-09-19 2012-10-16 Imatx, Inc. Method for bone structure prognosis and simulated bone remodeling
AU2004274003A1 (en) 2003-09-19 2005-03-31 Imaging Therapeutics, Inc. Method for bone structure prognosis and simulated bone remodeling
CA2580726A1 (fr) 2004-09-16 2006-03-30 Imaging Therapeutics, Inc. Systeme et procede de prediction de futures fractures
GB2423687B (en) 2005-02-25 2010-04-28 Rapiscan Security Products Ltd X-ray security inspection machine
US7471764B2 (en) 2005-04-15 2008-12-30 Rapiscan Security Products, Inc. X-ray imaging system having improved weather resistance
US8213570B2 (en) 2006-02-27 2012-07-03 Rapiscan Systems, Inc. X-ray security inspection machine
US7526064B2 (en) * 2006-05-05 2009-04-28 Rapiscan Security Products, Inc. Multiple pass cargo inspection system
GB0609744D0 (en) * 2006-05-16 2006-06-28 Oxford Instr Analytical Ltd Method of determining the feasibility of a proposed x-ray structure analysis process
US7822253B2 (en) * 2006-07-12 2010-10-26 General Electric Company Methods and apparatus for BMD measuring
US7742568B2 (en) * 2007-06-09 2010-06-22 Spectrum San Diego, Inc. Automobile scanning system
GB0803641D0 (en) 2008-02-28 2008-04-02 Rapiscan Security Products Inc Scanning systems
US9036779B2 (en) 2008-02-28 2015-05-19 Rapiscan Systems, Inc. Dual mode X-ray vehicle scanning system
GB0803642D0 (en) 2008-02-28 2008-04-02 Rapiscan Security Products Inc Drive-through scanning systems
GB0803640D0 (en) 2008-02-28 2008-04-02 Rapiscan Security Products Inc Scanning systems
GB0803643D0 (en) 2008-02-28 2008-04-02 Rapiscan Security Products Inc Mobile scanning systems
US8682052B2 (en) 2008-03-05 2014-03-25 Conformis, Inc. Implants for altering wear patterns of articular surfaces
EP2303193A4 (fr) 2008-05-12 2012-03-21 Conformis Inc Dispositifs et procédés pour le traitement de facette et d autres articulations
GB0809109D0 (en) 2008-05-20 2008-06-25 Rapiscan Security Products Inc Scanner systems
GB0809107D0 (en) * 2008-05-20 2008-06-25 Rapiscan Security Products Inc Scannign systems
GB0809110D0 (en) 2008-05-20 2008-06-25 Rapiscan Security Products Inc Gantry scanner systems
GB0810638D0 (en) 2008-06-11 2008-07-16 Rapiscan Security Products Inc Photomultiplier and detection systems
US8963094B2 (en) 2008-06-11 2015-02-24 Rapiscan Systems, Inc. Composite gamma-neutron detection system
US8939917B2 (en) 2009-02-13 2015-01-27 Imatx, Inc. Methods and devices for quantitative analysis of bone and cartilage
WO2010099231A2 (fr) 2009-02-24 2010-09-02 Conformis, Inc. Systèmes automatisés de fabrication d'implants orthopédiques spécifiques au patient et instrumentation
US9310323B2 (en) 2009-05-16 2016-04-12 Rapiscan Systems, Inc. Systems and methods for high-Z threat alarm resolution
DE102010038544A1 (de) * 2009-10-19 2011-04-21 Robert Bosch Gmbh Sensorvorrichtung für eine Verpackungsmaschine
US8314394B1 (en) 2009-11-04 2012-11-20 Science Applications International Corporation System and method for three-dimensional imaging using scattering from annihilation coincidence photons
EP2509539B1 (fr) 2009-12-11 2020-07-01 ConforMIS, Inc. Implants orthopédiques mis au point pour un patient et spécifiques à un patient
WO2012109273A2 (fr) 2011-02-08 2012-08-16 Rapiscan Systems, Inc. Surveillance secrète utilisant une détection à modalité multiple
EP2754419B1 (fr) 2011-02-15 2024-02-07 ConforMIS, Inc. Implants orthopédiques améliorés et adaptés au patient
US9218933B2 (en) 2011-06-09 2015-12-22 Rapidscan Systems, Inc. Low-dose radiographic imaging system
IN2014DN06514A (fr) 2012-02-03 2015-06-12 Rapiscan Systems Inc
US10670740B2 (en) 2012-02-14 2020-06-02 American Science And Engineering, Inc. Spectral discrimination using wavelength-shifting fiber-coupled scintillation detectors
EP2941775A4 (fr) 2013-01-07 2016-08-24 Rapiscan Systems Inc Dispositif de balayage à rayons x à réseau de détecteurs de discrimination d'énergie partielle
PL2952068T3 (pl) 2013-01-31 2021-07-26 Rapiscan Systems, Inc. Przenośny system kontroli bezpieczeństwa
FR3015681B1 (fr) * 2013-12-20 2016-04-22 Commissariat Energie Atomique Methode de mesure du numero atomique effectif d'un materiau
US9557427B2 (en) 2014-01-08 2017-01-31 Rapiscan Systems, Inc. Thin gap chamber neutron detectors
US9903964B2 (en) 2014-05-02 2018-02-27 Koninklike Philips N.V. Method to calibrate a photon detector, absorption filter assembly and imaging apparatus
EP3271709B1 (fr) 2015-03-20 2022-09-21 Rapiscan Systems, Inc. Système d'inspection par rétrodiffusion portable portatif
US10345479B2 (en) 2015-09-16 2019-07-09 Rapiscan Systems, Inc. Portable X-ray scanner
US10634797B2 (en) 2017-07-07 2020-04-28 International Business Machines Corporation Real time X-ray dosimeter using diodes with variable thickness degrader
EP3811117A4 (fr) 2018-06-20 2022-03-16 American Science & Engineering, Inc. Détecteurs de scintillation couplés à une feuille à décalage de longueur d'onde
US11175245B1 (en) 2020-06-15 2021-11-16 American Science And Engineering, Inc. Scatter X-ray imaging with adaptive scanning beam intensity
US11340361B1 (en) 2020-11-23 2022-05-24 American Science And Engineering, Inc. Wireless transmission detector panel for an X-ray scanner
JP7335284B2 (ja) * 2021-01-15 2023-08-29 アンリツ株式会社 Ｘ線検査装置
GB2618277A (en) 2021-02-23 2023-11-01 Rapiscan Systems Inc Systems and methods for eliminating cross-talk in scanning systems having multiple x-ray sources

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2004437A (en) * 1977-09-14 1979-03-28 Philips Nv Device for computer aided tomography
US4400827A (en) * 1981-11-13 1983-08-23 Spears James R Method and apparatus for calibrating rapid sequence radiography
EP0334762A1 (fr) * 1988-03-25 1989-09-27 General Electric Cgr S.A. Procédé et système d'étalonnage d'un scanner à rayons X en utilisant un seul étalon non circulaire
EP0402244A2 (fr) * 1989-06-09 1990-12-12 General Electric Cgr S.A. Procédé d'étalonnage d'un système radiologique et de mesure de l'épaisseur équivalente d'un objet
WO1992002892A1 (fr) * 1990-08-10 1992-02-20 Vivid Technologies, Inc. Dispositif et procede d'inspection de bagages et d'autres objets
AT394654B (de) * 1988-08-08 1992-05-25 Inst Med Physik Vet Med Uni Wi Testphantom zur laufenden kontrolle der bildguete in der roentgendiagnostik

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5044002A (en) * 1986-07-14 1991-08-27 Hologic, Inc. Baggage inspection and the like
NZ237767A (en) * 1992-04-09 1994-09-27 Inst Geolog Nuclear Sciences Luggage scanning by fast neutrons and gamma radiation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2004437A (en) * 1977-09-14 1979-03-28 Philips Nv Device for computer aided tomography
US4400827A (en) * 1981-11-13 1983-08-23 Spears James R Method and apparatus for calibrating rapid sequence radiography
EP0334762A1 (fr) * 1988-03-25 1989-09-27 General Electric Cgr S.A. Procédé et système d'étalonnage d'un scanner à rayons X en utilisant un seul étalon non circulaire
AT394654B (de) * 1988-08-08 1992-05-25 Inst Med Physik Vet Med Uni Wi Testphantom zur laufenden kontrolle der bildguete in der roentgendiagnostik
EP0402244A2 (fr) * 1989-06-09 1990-12-12 General Electric Cgr S.A. Procédé d'étalonnage d'un système radiologique et de mesure de l'épaisseur équivalente d'un objet
WO1992002892A1 (fr) * 1990-08-10 1992-02-20 Vivid Technologies, Inc. Dispositif et procede d'inspection de bagages et d'autres objets

Also Published As

Publication number Publication date
DE69401433D1 (de) 1997-02-20
FR2705785B1 (fr) 1995-08-25
NO954819D0 (no) 1995-11-27
NO954819L (no) 1996-01-26
FI955716A (fi) 1996-01-26
FR2705785A1 (fr) 1994-12-02
EP0700516A1 (fr) 1996-03-13
AU6800394A (en) 1994-12-20
EP0700516B1 (fr) 1997-01-08
DE69401433T2 (de) 1997-06-12
FI955716A0 (fi) 1995-11-27
US5687210A (en) 1997-11-11
JPH09500719A (ja) 1997-01-21
ATE147510T1 (de) 1997-01-15

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
EP2791620B1 (fr) Mesure par interferometrie atomique
EP3106864B1 (fr) Caracterisation d'un echantillon par representation en base de materiaux
FR2961904A1 (fr) Procede d'identification de materiaux a partir de radiographies x multi energies
EP0470909B1 (fr) Appareil de détection nucléaire, notamment du genre gamma-caméra, à filtres de déconvolution
FR2953298A1 (fr) Procede de correction du phenomene d'empilement applique a des spectres de rayonnement x acquis a l'aide d'un capteur spectrometrique
FR2738693A1 (fr) Systeme de traitement d'impulsions provenant de l'interaction d'une particule gamma avec un detecteur de rayonnement cdte
EP0334762B1 (fr) Procédé et système d'étalonnage d'un scanner à rayons X en utilisant un seul étalon non circulaire
EP0414587A1 (fr) Procédé et chaîne de spectrométrie gamma
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
Christensen et al. The host galaxy of GRB 990712
FR2695730A1 (fr) Méthode et appareil pour traiter un signal numérique amélioré utilisant une transformée de Fourier.
EP3617751A1 (fr) Procédé d'auto-calibration d'un dispositif de détection de rayonnements ionisants
FR3038722B1 (fr) Procede d’identification d’un materiau
EP2225588B1 (fr) Procede de traitement de signaux de mesure de chambres a fission
CA2163884C (fr) Procede et dispositif pour la reconnaissance de materiaux determines dans la composition d'un objet
WO2021122772A1 (fr) Procédé d'évaluation de la teneur en uranium par spectrométrie gamma dans un puits de forage et dispositif associé
EP1105751A1 (fr) Procede de simulation de la reponse d'un detecteur de rayonnements emis par des objets radioactifs et procede de controle d'elements de combustible nucleaire utilisant cette simulation
FR2787573A1 (fr) Controle par radiographie d'un objet possedant un reseau cristallin
FR2476474A1 (fr) Dispositif de tomodensitometrie et procede d'exploration et de reconstruction d'images utilisant un tel dispositif
WO2011036282A1 (fr) Procede de determination du rapport isotopique de matiere fissile contenue dans une chambre a fission
CA2166985C (fr) Procede et dispositif pour detecter la presence d'un objet, comportant un materiau donne, non accessible a la vue
EP2737339A1 (fr) Dispositif et procede de traitement pour une mesure spectrometrique d'un flux de photons

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AU BB BG BR BY CA CN CZ FI HU JP KP KR KZ LK MG MN MW NO NZ PL RO RU SD SK UA US VN

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN ML MR NE SN TD TG

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2163882

Country of ref document: CA

Ref document number: 955716

Country of ref document: FI

WWE Wipo information: entry into national phase

Ref document number: 1994916291

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1994916291

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 08553427

Country of ref document: US

WWG Wipo information: grant in national office

Ref document number: 1994916291

Country of ref document: EP