Brevet FR2968403A1 — publié le 20120608

L'invention porte sur un procédé de traitement d'un spectre de rayonnement diffusé par un matériau exposé à un rayonnement incident, en vue de l'extraction du spectre de rayonnement diffusé premier. Elle porte également sur un dispositif de traitement associé et sur un programme d'ordinateur pour un tel procédé de traitement. Le domaine d'application de l'invention s'étend en tout premier lieu à la spectrométrie de diffusion de rayons X ou gamma, avec notamment l'analyse des matériaux, mais comprend également d'autres spectrométries de diffusion. Une telle spectrométrie peut être utilisée dans l'étude des matériaux, par exemple la détection des explosifs. La spectrométrie de diffusion de rayons X est basée sur l'exposition d'un matériau à un rayonnement X incident d'énergie de quelques dizaines à quelques centaines de keV. Les photons X, quand ils rencontrent le matériau sur lequel ils sont projetés, induisent différents types d'interactions avec celui-ci : fluorescence ou conversion interne (effet photoélectrique au cours duquel le photon transfert toute son énergie au matériau qui la restitue ensuite), diffusion inélastique (ou effet Compton qui inclut un changement de la direction du photon ainsi qu'une diminution de son énergie), création de paires positon - électron (uniquement pour des rayons X de très grande énergie non considérés dans la présente invention), ou diffusion Rayleigh (ou diffusion élastique, minoritaire aux énergies considérées). L'invention s'applique préférentiellement au traitement de signaux produits par des tubes X délivrant des photons dont l'énergie est comprise entre zéro et 300 keV. De tels photons diffusent dans la matière. Certaines caractéristiques des matériaux étudiés (coefficient d'absorption µ (E), densité, rapport Â entre le numéro atomique Z et le nombre de masse atomique A, ainsi que la composition chimique notamment) peuvent être déterminées sur la base de connaissances théoriques et de l'obtention du spectre diffusé premier du matériau exposé à un rayonnement X, c'est-à-dire le spectre de diffusion qui serait obtenu dans une situation où chaque photon n'interagirait qu'une seule fois avec le matériau. La demande de brevet FR 09 55011, déposée le 20 juillet 2009, montre par exemple comment on peut déterminer certaines caractéristiques de matériau en analysant le spectre en énergie d'un rayonnement diffusé par ce matériau, ce dernier étant soumis à une irradiation par un rayonnement X. Or, les spectres de diffusion X comportent une composante 10 importante de photons diffusés ayant interagi plusieurs fois avec le matériau. Cette composante est appelée le spectre de diffusé multiple. Certaines informations, notamment la densité du matériau, peuvent être obtenues sur la base du spectre de diffusé total car le caractère atténuant du matériau affecte d'une manière similaire les deux composantes, première et 15 multiple, du spectre de diffusion. Pour l'obtention d'une meilleure estimation de cette densité ou d'autres informations physico-chimiques, on constate que l'utilisation du spectre de diffusion total mène à des résultats imparfaits. Le document WO2007/007247 enseigne d'utiliser un système où l'objet étudié est placé entre une source de rayonnement et une matrice de 20 détecteurs pixellisés. Une calibration du détecteur est effectuée sur la base des lignes caractéristiques de la source, et une évaluation du diffusé multiple est faite pour chacune de ces lignes caractéristiques avant d'être extrapolée à tout le spectre. L'invention se place dans un autre contexte, puisqu'elle s'applique en 25 premier lieu à un système d'analyse avec une source de rayonnement fortement collimatée, et un détecteur placé dans le même demi-espace que la source vis-à-vis de la surface de matériau étudiée, et également fortement collimaté. L'invention a pour objet de remédier, dans le contexte de cet 30 agencement analytique, au problème évoqué ci-dessus, en proposant un procédé de correction du spectre de rayonnement diffusé total pour en extraire le rayonnement diffusé premier, c'est-à-dire un rayonnement n'ayant diffusé qu'une seule fois dans le matériau analysé avant d'être détecté. Pour ce faire, il est proposé un procédé d'extraction d'un spectre de rayonnement diffusé premier hors d'un spectre de diffusion de rayonnement diffusé, selon un angle de diffusion, issu d'un matériau exposé à un rayonnement incident à travers une surface, le procédé comprenant au moins l'application d'une fonction spectrale de réponse correspondant au spectre du rayonnement diffusé multiple lorsqu'on détecte un photon d'énergie donnée appartenant au spectre de diffusion premier. Grâce à ce procédé, on extrait aisément le spectre de diffusé premier et on accède ainsi à de nombreuses informations physico-chimiques. Selon un mode de réalisation particulièrement avantageux, ladite étape d'application comprend l'application d'au moins deux fonctions spectrales de réponse, les au moins deux fonctions spectrales de réponse étant organisées sous forme d'une matrice (M), dont chaque valeur ai correspond au nombre de photons détectés, d'énergie Ei, constituant le rayonnement diffusé multiple, lorsqu'on détecte un photon, d'énergie Ei, du rayonnement diffusé premier. Ce mode de calcul offre l'avantage de l'efficacité, notamment en termes de programmation. Selon un aspect de l'invention, la fonction spectrale de réponse fournit un spectre discrétisé sur un nombre fini de plages d'énergie. Selon un aspect avantageux de l'invention, le procédé comprend un processus itératif dont chaque pas comprend une estimation du spectre de diffusion multiple d'après une estimation précédente du spectre de diffusion premier, et une nouvelle estimation du spectre diffusion premier, par soustraction dudit spectre de diffusion multiple estimé du spectre de diffusion détecté, ledit processus itératif étant poursuivi jusqu'à la satisfaction d'un critère de convergence portant sur les spectres diffusés multiples estimés successifs. Selon un mode de réalisation, la fonction spectrale de réponse est le résultat d'une interpolation d'une pluralité de fonctions spectrales de réponse à une exposition monochromatique de différents matériaux pour la profondeur donnée en considérant uniquement la densité du matériau. Selon un autre mode de réalisation, la fonction spectrale de réponse est le résultat d'une interpolation d'une pluralité de fonctions spectrales de réponse du matériau pour différentes profondeurs. Il est aussi proposé que la fonction spectrale de réponse est préalablement obtenue par une acquisition expérimentale d'une réponse du matériau pour la profondeur donnée. Selon un aspect de l'invention, on procède aussi à une étape préliminaire d'estimation d'un premier spectre de diffusé multiple estimé, par multiplication d'un spectre de diffusion total mesuré par un coefficient multiplicatif dépendant de la profondeur d'inspection. L'invention consiste également en un dispositif d'extraction d'un spectre de rayonnement diffusé premier hors d'un spectre de diffusion de rayonnement diffusé, selon un angle de diffusion, issu d'un matériau exposé à un rayonnement incident à travers une surface, le dispositif comprenant au moins des moyens d'application d'une fonction spectrale de réponse correspondant au spectre du rayonnement diffusé multiple lorsqu'on détecte un photon d'énergie donnée appartenant au spectre de diffusion premier, et un dispositif d'analyse spectrométrique comprenant un tel dispositif d'extraction de spectre de rayonnement diffusé premier. Il est également proposé un programme d'ordinateur comprenant une suite d'instructions aptes, lorsqu'elles sont exécutées par un microprocesseur, à mettre en oeuvre un procédé selon l'invention. Celle-ci va maintenant être décrite en relation avec les figures annexées, à titre d'illustrations non limitatives. La figure 1 représente un mode de réalisation d'un dispositif selon l'invention. La figure 2 représente la décomposition générale des spectres de diffusion. La figure 3 représente une réponse spectrale à une exposition monochromatique d'un matériau telle qu'envisagée dans l'invention. La figure 4 représente une étape d'un mode de réalisation d'un procédé selon l'invention. La figure 5 représente une étape ultérieure du même procédé. La figure 6 représente une étape d'un premier mode de réalisation particulier de l'invention. Les figures 7 et 8 représentent des résultats obtenus avec le premier mode de réalisation particulier de l'invention. La figure 9 représente deux fonctions mises en oeuvre dans un deuxième mode de réalisation particulier de l'invention. La figure 10 représente des résultats obtenus pour le deuxième mode de réalisation particulier de l'invention. La figure 11 représente la confrontation entre des résultats obtenus avec le procédé selon l'invention et des mesures obtenues par ailleurs. La figure 1 représente l'agencement matériel mis en oeuvre au cours d'un procédé selon l'invention. Une source 100 émet un rayonnement X 150 qui est collimaté par un collimateur 110. Le rayonnement 150 forme un faisceau dont l'angle d'ouverture est très faible, typiquement quelques degrés, autour de l'axe central du collimateur de la source. La source 100 est monoénergétique, ou polyénergétique, et émet des photons de manière continue sur toute une gamme d'énergie. Le faisceau est projeté sur un volume de matériau 200 par une surface du matériau 210. Les photons X interagissent avec le matériau suivant les différents phénomènes physiques évoqués en introduction. Certains des photons interagissent très près de la surface alors que d'autres s'enfoncent profondément dans le matériau. Un détecteur 300 équipé d'un collimateur 310 observe la surface 210 du matériau 200 et mesure l'énergie des photons que son angle de collimation lui permet de recevoir. Le positionnement relatif de la source 100 et du détecteur 300 ainsi que leur positionnement vis-à-vis du matériau 200 et leurs angles de collimation respectifs définissent une profondeur d'analyse dans le matériau 200, qui est telle que les photons mesurés par le détecteur 300 proviennent d'un volume du matériau 200 dit volume d'inspection. Ce volume est déterminé par l'intersection entre les angles solides sous lesquels la source et le détecteur voient le volume de matériau 200. Les axes de collimation respectifs de la source et du détecteur se croisent dans le matériau analysé. La profondeur d'analyse, ou profondeur d'inspection, correspond à la distance entre la surface du matériau et cette intersection. L'angle entre les deux axes de collimation est noté O. Le détecteur détecte les photons émis par la source et ayant diffusé selon un angle 0' _ - 0 dans le volume d'inspection. Lorsque 0' = 0, il n'y a pas de diffusion inélastique : le photon suit une trajectoire rectiligne. La source 100 et le détecteur 300 étant disposés du même côté de la surface 210, on parle de spectre de diffusion en arrière ou spectre de rétrodiffusion. En figure 2, on a représenté le spectre de diffusion mesuré pour un matériau 200 donné. En abscisses est figurée l'énergie des photons qui s'étale typiquement de 20 à 100 keV, et en ordonnées, on a représenté l'intensité du rayonnement détecté pour une énergie donnée en unité arbitraire. Le spectre de diffusion totale 1000 se décompose en un spectre de diffusion premier 1100 et un spectre de diffusion multiple 1200 comme évoqué en introduction. En figure 1, on a représenté la trajectoire d'un photon diffusé premier 410. Par définition, sa trajectoire ne change qu'une fois : il ne subit qu'une seule diffusion avant d'atteindre le détecteur. Celle-ci se fait selon l'angle de diffusion 0', et se produit dans le volume d'inspection. Les photons tels que celui précédemment évoqué forment le rayonnement de diffusion premier selon l'angle 0', constitué de photons n'ayant diffusé qu'une fois, selon cet angle de diffusion 0', dans le volume d'inspection. On a également représenté la trajectoire d'un photon diffusé multiple qui fait l'objet de plusieurs interactions avec le matériau 200 et qui change de trajectoire à chacune de ces interactions. Il subit alors de multiples diffusions avant d'atteindre le détecteur. On parle alors d'un rayonnement diffusé multiple, car constitué de photons ayant diffusé plus d'une fois avant d'atteindre le détecteur. On précise que d'autres photons ressortent du volume de matériau 200, à l'instar des photons 410 et 420, sous d'autres angles et ne sont pas mesurés par le détecteur 300. Ils n'apparaissent donc pas sur les spectres de la figure 2. Une partie du rayonnement diffusé premier et du rayonnement diffusé multiple est détectée par le détecteur. Ce dernier n'est pas capable de distinguer ces deux composantes. Ainsi, le détecteur détermine un spectre de rayonnement en énergie, du rayonnement diffusé total, noté DT. Ce dernier est composé d'une composante correspondant au spectre de diffusion premier, noté DP, et précédemment défini et d'une composante correspondant au spectre de diffusion multiple, noté DM, précédemment défini. Ainsi, le spectre de diffusion total est la somme du spectre de diffusion premier et du spectre de diffusion multiple. DT = DP + DM En figure 3, on a représenté la réponse spectrale de rayonnement diffusé en réaction à une exposition monochromatique du matériau 200 dans l'agencement de la figure 1. Le matériau utilisé pour l'exemple est du Delrin (marque déposée) et l'excitation monochromatique 2000 est effectuée à 150 keV. Le volume d'inspection est à une profondeur de 4 cm, l'angle entre le faisceau incident et l'axe de collimation du détecteur étant de 120°. On constate que, comme attendu, la réponse spectrale 2100 est composée uniquement d'énergies inférieures à celle de l'excitation monochromatique 2000. Dans un système d'analyse tel qu'en figure 1, c'est-à-dire utilisant une mesure de rétrodiffusion et une collimation à la fois de la source et du détecteur, le spectre de diffusé premier en réponse à une excitation monochromatique est lié de manière bijective à l'énergie de l'activation, et est essentiellement monochromatique. Il apparaît sous la référence 2200 en figure 3. Ainsi, en dehors de l'énergie mesurée à très grande proximité d'une valeur d'énergie correspondant à l'ensemble du spectre de diffusion premier 2200, qui en figure 3 est de 104,14 keV, l'ensemble du spectre de diffusion 2100 est constitué de la composante de diffusé multiple DM. On représente ce spectre DM sous une forme discrétisée sur un (a 1 nombre limité n de canaux d'énergie de détection, par un vecteur colonne a. \an / a; étant l'énergie du rayonnement diffusé multiple à l'énergie j. Le vecteur dépend naturellement du matériau et de la profondeur d'analyse. Le matériau peut être irradié selon un grand nombre d'énergies monochromatiques Eirradiation-i, chaque énergie Eirradiation-i étant telle que le rayonnement diffusé premier, selon l'angle 0', détecté par le détecteur a une énergie Ei. L'ensemble des réponses du matériau à des excitations monochromatiques, discrétisées pour un nombre déterminé de canaux d'énergie de détection est présenté sous forme de la matrice de corrélation. (all al. alp M= ail ail anl a \ np / Dans cette matrice, chaque valeur ai, correspond à la proportion de photons détectés constituant le rayonnement diffusé multiple, d'énergie Ei lorsqu'on détecte un photon, d'énergie Ej, du rayonnement diffusé premier. De cette façon, chaque somme de coefficients ai, selon une colonne j de cette matrice correspond au nombre de photons diffusés multiples détectés Nmultiples_j lorsqu'un photon diffusé premier à l'énergie Ej est détecté. Nmultiplea-; _ L agi, i Chaque colonne Ci de la matrice M représente le spectre en énergie du rayonnement diffusé multiple détecté lorsqu'on détecte un photon diffusé premier d'énergie Ej. Cette fonction est appelée fonction de réponse spectrale correspondant à la détection d'un photon diffusé premier d'énergie Ej. Le procédé selon l'invention utilise la matrice de corrélation M pour extraire de manière itérative la composante de diffusés multiples dans le spectre de diffusion mesuré expérimentalement. Il s'agit d'un procédé itératif, permettant une estimation du spectre de diffusion premier DP. Ainsi, à chaque itération k, on réalise une estimation du spectre de diffusion premier notée DPi. Cette estimation est réalisée en utilisant la matrice de corrélation précédemment définie, selon l'algorithme suivant : - DMk : estimation du spectre de diffusion multiple lors de l'itération i ; - DPk_l : estimation du spectre de diffusion premier établie à l'itération k-1 lorsque k 1. Exceptionnellement, lorsque k = 1, DPk_l est l'estimation initiale du spectre de diffusion premier ; - DT : spectre de diffusion total, c'est-à-dire spectre de diffusion mesuré par le détecteur. Les itérations se poursuivent soit selon un nombre d'itérations prédéterminé, soit jusqu'à ce qu'un critère de convergence soit atteint, comme 15 cela sera décrit plus bas. Le procédé itératif comporte ainsi - une estimation du spectre diffusé premier initial (étape d'initialisation) - une estimation du spectre de diffusion multiple DMk correspondant 20 au spectre de diffusion premier de l'itération précédente DPk_,, ou, lors de la première itération, de l'estimation initiale du spectre diffusé premier (DP0). Cette estimation est réalisée à l'aide de la matrice de corrélation M, - une nouvelle estimation du spectre de diffusion premier DPk, effectuée en soustrayant le spectre diffusé multiple estimé DMk au spectre de 25 diffusion mesuré DT, - une nouvelle étape d'itération au cours de laquelle le spectre de diffusé multiple est évalué sur la base de l'estimation du spectre de diffusé premier. Ainsi, suivant la méthode d'inversion du point fixe, on sépare les 30 deux composantes dans le spectre mesuré physiquement. Le spectre de diffusé premier initial, noté DPk-o est tout d'abord estimé sur la simple base du spectre de diffusion total DT mesuré, multiplié par DMk = M * DPk_l (1) DPk = DT - DMk (2) Avec : un facteur d'atténuation globale, noté a, compris entre 0.3 et 1 en fonction de la profondeur d'inspection, pour des matériaux de densité comprise entre 0.7 et 2.2. Pour une profondeur d'inspection de 2 cm, une valeur de 0.6 est utilisée, et pour une profondeur de 4 cm, la valeur est de 0.45. Ces valeurs peuvent être ajustées à chaque mise en oeuvre. Ainsi, DPk=o = a DT (3) Connaissant le spectre diffusé total DT (il s'agit du spectre mesuré par le détecteur), et la valeur initiale d'estimation du spectre diffusé premier DPk-o, on réalise les étapes (1) et (2) précédemment décrites, de façon itérative, jusqu'à ce que le critère de convergence soit atteint ou jusqu'à ce que le nombre d'itérations atteigne le nombre prédéterminé. Un critère de convergence possible est une distance entre les estimations successives de deux spectres de diffusion premier, DPk_ et DPk, 15 cette distance pouvant être exprimée selon l'expression suivante. DPk = .\iE (DPk (E) - DPk_I (E))2 E (4) E Dans cette expression, le terme E désigne l'énergie et le terme s désigne le critère de convergence. A chaque itération, le spectre de diffusé premier est représenté sous dpp 20 forme d'un vecteur . dp2 , et le diffusé multiple est également représenté sous \dpp /k dml dm2 forme d'un vecteur . Ce dernier se déduit sous forme de calcul matriciel \dmn /k dmn a~ clip dpp par l'opération dm2 a21 a22 dp2 \dmn /k .. anp dpp Les pas de discrétisation en énergie ainsi que la dimension des vecteurs précédemment décrits sont de préférence identiques. Egalement, le nombre de canaux n de discrétisation du spectre de diffusion multiple est, de préférence, choisi égal au nombre de canaux p du spectre de diffusion premier. Dans ce cas, la matrice M est carrée. Le procédé itératif est illustré aux figures 4 et 5. La figure 4 représente la multiplication d'une intensité dp; du spectre de diffusé premier estimé 3000 à l'énergie i = 40 keV avec une colonne ai, j variant de 1 à n, représentée sous la référence 3100 en partie basse de la figure, le résultat étant une fraction 3200 du spectre de diffusé multiple estimé 3300, visible sur la (dmi dm2 partie droite de la figure 4 et exprimée sous la forme du vecteur , que l'on ~dmn i additionne aux fractions de spectre calculées pour les valeurs d'énergie inférieures à i pour obtenir un total intermédiaire 3300. La figure 5 représente la multiplication d'une autre intensité dp; du spectre de diffusé premier à l'énergie i avec une colonne a;,, référencée 3110, le résultat étant une deuxième fraction du spectre de diffusé multiple estimé 3210, les différentes fractions étant ensuite additionnées pour obtenir le spectre de diffusé multiple estimé 3300. Un élément important de l'invention est l'utilisation de la matrice de corrélation M, permettant d'estimer le spectre de diffusion multiple DMk correspondant à un spectre de diffusion premier DPk donné. Cette matrice doit être établie pour une position donnée de la source et du détecteur, par rapport à l'objet analysé. En outre, elle dépend des paramètres suivants : - l'angle 0 entre les axes de collimation du détecteur et de la source, qui conditionne l'angle de diffusion 0«0 = - 0» - la nature du matériau constituant l'objet analysé - la profondeur du volume d'inspection, laquelle est déterminée par l'intersection des axes de collimation du détecteur et de la source. La matrice de corrélation M peut être construite expérimentalement, en utilisant, successivement, différentes sources de rayons X ou gamma monoénergétiques d'énergie Eirradiation-j, comme des sources radioactives dont les raies sont connues, ou un rayonnement synchrotron, qui a de plus l'avantage d'être de grande énergie. On mesure alors un spectre de diffusion total, tel que celui apparaissant en figure 6, et référencé 4100. On procède alors selon trois étapes : Tout d'abord, on détermine le pic de diffusion premier. A partir du spectre de diffusion total 4100, on isole le pic 4200 correspondant au rayonnement de diffusion premier détecté. Ce pic s'étend de part et d'autre d'une énergie connue, Ej correspondant à l'énergie de diffusion Eirradiation_j du rayonnement monochromatique émis par la source selon l'angle de diffusion 0'. Ce pic est d'autant plus fin que la source et le détecteur sont collimatés. De préférence, la largeur à mi-hauteur de ce pic correspond à une bande énergétique inférieure à quelques keV. Ce pic 4200 correspond à la composante du rayonnement de diffusion premier dans le spectre de diffusion total détecté. Au cours d'une deuxième étape, on détermine le spectre de diffusion multiple. En retirant la composante de diffusion premier 4200 du spectre de diffusion total 4100, on obtient le spectre de diffusion multiple 4300. Pour cela, on détermine le fond "sous le pic", c'est-à-dire la courbe continue, par exemple linéaire, joignant les deux extrémités du spectre 4300 de part et d'autre de la base du pic 4200. On peut utiliser des algorithmes de détermination du signal sous un pic, couramment mis en oeuvre dans le domaine de la spectrométrie X ou gamma. Enfin, au cours d'une troisième étape, on procède à une normalisation. On divise le spectre de diffusion multiple 4300 par la quantité de photons correspondant au pic de diffusion premier 4200. Le spectre de diffusion multiple ainsi normalisé correspond alors au spectre de diffusion multiple obtenu lorsqu'un photon du rayonnement diffusé premier est détecté. Ce spectre correspond aussi à une colonne Ci de la matrice de corrélation M, dont l'indice j correspond à l'énergie Ej du pic de diffusion premier observé, cette dernière étant déterminée par l'énergie Eirradiation-j de la source de rayonnement et l'angle de diffusion 0'. Les trois étapes ainsi décrites peuvent être réalisées non pas de façon expérimentale, mais par simulation, en utilisant par exemple un code de transport de particules de type MonteCarlo connu de l'homme du métier, et notamment MCNP, GEANT ou PENELOPE. En outre, la matrice de corrélation peut être constituée en combinant simulation et expérimentation.. En renouvelant cette opération en utilisant, ou en simulant, successivement différentes sources monoénergétiques, il est possible d'obtenir plusieurs colonnes de la matrice, correspondant à autant d'énergies de diffusion premier Ei, ces dernières étant calculées d'après la connaissance de l'angle de diffusion 0' et de l'énergie d'émission Eirradiation-i de la source mise en oeuvre. Ces différentes colonnes peuvent être interpolées entre elles, par exemple de façon linéaire, afin de constituer une matrice complète, les pas de discrétisation en énergie entre chaque ligne et chaque colonne correspondant respectivement au pas de discrétisation en énergie du vecteur représentant un spectre de diffusion multiple et au pas de discrétisation en énergie du vecteur représentant le spectre de diffusion premier. Une telle matrice est réalisée pour différentes configurations de 20 mesure, en faisant varier les paramètres suivants : - nature du matériau étalon - position, dans le matériau étalon, de l'intersection des axes de collimation du détecteur et de la source, cette position correspondant à la profondeur d'inspection. 25 Ainsi, pour un angle de diffusion donné, on obtient une pluralité de matrices de corrélation M(nature, profondeur). Par exemple, lorsque l'angle de diffusion 0' est de 120°, correspondant à un angle séparant les axes de collimation 0 de 60°, on peut simuler le spectre de diffusion multiple correspondant à une source 30 monoénergétique dont l'énergie varie entre 10 et 150 keV, par pas de 10 keV. Quinze simulations sont alors nécessaires. Cela correspond à des photons de diffusion premier d'énergies comprise entre 9.7 keV (lorsque l'énergie de la source est égale à 10 keV) et 104.14 keV (lorsque l'énergie de la source est égale à 150 keV). La matrice est ensuite complétée par interpolation. On constitue de telles matrices pour différents matériaux étalons, ces derniers pouvant être le polyéthylène, l'eau, le Plexiglas (marque déposée), le nylon (marque déposée), le Delrin, la cellulose, le Kynar (marque déposée), le nitrate d'ammonium ou le Téflon (marque déposée). Pour chaque matériau, on détermine une matrice de corrélation par profondeur d'inspection, cette dernière variant par exemple de 0 à 7 cm, par pas de 0.5 cm. Pour un même matériau, on peut également interpoler des matrices selon la profondeur. La figure 7 présente un exemple de mise en oeuvre du procédé itératif d'estimation du spectre de diffusion premier, précédemment décrit, en utilisant un matériau en plexiglas et une profondeur de 2 cm. On constate une excellente convergence de cette estimation. La convergence est ici estimée par calcul d'une distance normalisée entre deux spectres de diffusé premier estimés successifs et comparaison entre cette distance et une valeur de seuil s suivant une formule du type DPk-DPk_i =.\jE(DPk(E)-DPk_I(E))2 matériaux à étudier à leur densité, qui détermine l'atténuation sur l'ensemble de la gamme spectrale, pour la préparation de la matrice de corrélation. La figure 9 présente la répartition spectrale de la diffusion multiple pour un diffusé premier à 104.14 keV, qui ne dépend pas du matériau, et la figure 10 présente la fonction d'atténuation pour la même valeur d'énergie de diffusé premier. Des interpolations peuvent être également réalisées pour les valeurs d'énergie incidente et la profondeur, de manière à adapter la matrice M utilisée effectivement aux valeurs précises de la situation étudiée. Le procédé de traitement de spectre décrit dans cette demande peut être appliqué sur un matériau inconnu, mais de densité connue, ou estimée. On sélectionne une matrice déjà constituée, correspondant à un matériau dont la densité est proche de celle du matériau inconnu. L'extraction du spectre de diffusion premier peut permettre de déterminer plus précisément la nature de ce matériau, par exemple en suivant le procédé décrit dans la demande de brevet FR 09 55011, déposée le 20 juillet 2009, précédemment citée. Il est alors possible de refaire une extraction du spectre de diffusion premier à l'aide d'une matrice tenant compte de la nature précise du matériau, telle que celle-ci a été précédemment déterminée. La figure 11 présente l'excellente convergence constatée pour la détermination du spectre de diffusé premier dans ce mode de réalisation, dans le cas de l'étude d'un volume de plexiglas à 3 cm de profondeur. La convergence est atteinte en moins de 10 itérations, avec un phénomène d'oscillations autour du spectre final. La figure 12 présente les résultats de la décomposition en diffusé premier et diffusé multiple d'un spectre de diffusion total obtenu avec du Delrin étudié à 3 cm de profondeur. Les résultats sont confrontés avec des données simulées du spectre de diffusion total, du spectre de diffusion premier, du spectre de diffusion multiple. L'invention ne se limite pas au mode de réalisation décrit mais s'étend à toutes les variantes à la portée de l'homme du métier. REVENDICATIONS1. Procédé d'extraction d'un spectre de rayonnement diffusé premier (1100) hors d'un spectre de diffusion (1000) de rayonnement diffusé, selon un angle de diffusion, issu d'un matériau exposé (200) à un rayonnement incident (150) à travers une surface (210), le procédé comprenant au moins l'application d'une fonction spectrale de réponse correspondant au spectre du rayonnement diffusé multiple lorsqu'on détecte un photon d'énergie donnée appartenant au spectre de diffusion premier. 2. Procédé d'extraction selon la revendication 1, caractérisé en ce ladite étape d'application comprend l'application d'au moins deux fonctions spectrales de réponse (3100, 3110), les au moins deux fonctions spectrales de réponse étant organisées sous forme d'une matrice (M), dont chaque valeur au correspond au nombre de photons détectés, d'énergie Ei, constituant le rayonnement diffusé multiple, lorsqu'on détecte un photon, d'énergie Ei, du rayonnement diffusé premier. 3. Procédé d'extraction selon l'une des revendications précédentes, 20 caractérisé en ce que ladite fonction spectrale de réponse fournit un spectre discrétisé sur un nombre fini de plages d'énergie. 4. Procédé d'extraction selon l'une des revendications précédentes, caractérisé en ce qu'il comprend un processus itératif dont chaque pas 25 comprend une estimation du spectre de diffusion multiple d'après une estimation précédente du spectre de diffusion premier et une nouvelle estimation du spectre diffusion premier, par soustraction dudit spectre de diffusion multiple estimé du spectre de diffusion détecté, ledit processus itératif étant poursuivi jusqu'à la satisfaction d'un critère de convergence portant sur 30 les spectres diffusés multiples estimés (3300) successifs. 5. Procédé d'extraction selon l'une des revendications précédentes, caractérisé en ce que la fonction spectrale de réponse (3100, 3110) est le résultat d'une interpolation d'une pluralité de fonctions spectrales de réponse à une exposition monochromatique de différents matériaux pour la profondeur donnée en considérant uniquement la densité du matériau. 6. Procédé d'extraction selon l'une des revendications précédentes, caractérisé en ce que la fonction spectrale de réponse (3100, 3110) est le résultat d'une interpolation d'une pluralité de fonctions spectrales de réponse à une exposition monochromatique du matériau pour différentes profondeurs. 7. Procédé d'extraction selon l'une des revendications précédentes, caractérisé en ce que la fonction spectrale de réponse (3100, 3110) est préalablement obtenue par une acquisition expérimentale d'une réponse à une exposition monochromatique du matériau pour la profondeur donnée. 8. Procédé d'extraction selon l'une des revendications précédentes, comprenant de plus une étape préliminaire d'estimation d'un premier spectre de diffusé multiple estimé, par multiplication d'un spectre de diffusion total mesuré par un coefficient multiplicatif dépendant de la profondeur d'inspection. 9. Procédé d'extraction selon l'une des revendications précédentes, caractérisé en ce que le rayonnement est un rayonnement X. 10. Dispositif d'extraction d'un spectre de rayonnement diffusé premier (1100) hors d'un spectre de diffusion (1000) de rayonnement diffusé, selon un angle de diffusion, issu d'un matériau exposé (200) à un rayonnement incident (150) à travers une surface (210), le dispositif comprenant au moins des moyens d'application d'une fonction spectrale de réponse correspondant au spectre du rayonnement diffusé multiple lorsqu'on détecte un photon d'énergie donnée appartenant au spectre de diffusion premier.30 11. Dispositif d'analyse spectrométrique comprenant un dispositif d'extraction de spectre de rayonnement diffusé premier selon la revendication 10. 12. Programme d'ordinateur comprenant une suite d'instructions aptes, lorsqu'elles sont exécutées par un microprocesseur, à mettre en oeuvre un procédé selon l'une des revendications de procédé 1 à 9.10