FR 3023001 A1 20160101 FR 1456212 A 20140630 La présente invention concerne un procédé d'analyse d'un objet par diffractométrie, et un système de détection associé. La diffractométrie correspond à l'irradiation d'un matériau ou objet à l'aide 5 d'une source de rayonnement électromagnétique, de type X ou gamma, suivie de l'analyse du rayonnement de diffusion élastique à faible angle. L'expression « rayonnement de diffusion élastique à faible angle » désigne le rayonnement diffusé par le matériau ou objet, de manière cohérente, selon un angle inférieur à 15°, voire même inférieur à 10°, par rapport à la direction du rayonnement incident au matériau 10 ou objet. Au-delà en 10°, en effet, la diffusion élastique, ou diffusion Rayleigh, devient progressivement négligeable. Il est connu d'utiliser la diffractométrie pour détecter certaines substances cristallines telles que la plupart des explosifs ou de nombreuses autres structures dangereuses ou illégales. L'invention trouve ainsi une application dans le domaine de la sécurité, 15 plus particulièrement dans le domaine de la détection de matériaux explosifs dans un bagage. Elle est également utile dans le domaine médical, par exemple, pour la localisation d'une tumeur dans un sein. En effet, une publication de l'UCL (Pani, S. et al. « Characterization of breast tissue using energy-dispersive X-ray diffraction 20 computed tomography ». Applied Radiation and Isotopes 68, n° 10 (2010) : 1980-1987) a pu démontrer la possibilité de différencier (tissus adipeux, fibreux, tumeurs bénignes, adénofibromes, carcinomes...) des objets de tissus biologiques issus de biopsies du sein grâce à la mesure en diffusion de ces tissus. L'analyse du rayonnement diffusé selon un angle faible (à noter que le 25 terme « diffracté » est généralement utilisé pour un matériau cristallin tandis que le terme « diffusé » est généralement utilisé pour un matériau amorphe, mais ces deux termes sont ici indifféremment utilisés, de même que les termes diffusion et diffraction) par un matériau est une méthode d'analyse physico-chimique qui fournit des informations sur la structure de la matière permettant ainsi une meilleure 30 caractérisation des matériaux. Il est connu que l'analyse du spectre du rayonnement diffusé à faible angle, ou spectre de diffusion, permet d'établir une signature du matériau examiné. Pour les matériaux cristallins par exemple, lorsque la longueur d'onde des rayons X irradiants est du même ordre de grandeur que les distances inter- réticulaires (quelques angstdims), les rayons diffusés génèrent des interférences constructives ou destructives selon leur énergie et leur angle de diffusion. Les conditions pour lesquelles les interférences sont constructives sont déterminées par la loi de Bragg. Pour un matériau cristallin, cette loi lie la distance inter-réticulaire, l'énergie du rayonnement diffusé et l'angle de diffusion, selon l'équation suivante : E hkl n 2dhi, sin(8/ 2) Avec : dhk, : distance inter-réticulaire entre les plans cristallographiques du cristal irradié ; {hk1} : indices de Miller ; 0 : angle de diffusion, c'est-à-dire l'angle formé entre le rayonnement diffusé analysé et le faisceau incident au cristal irradié ; h : constante de Planck ; c: vitesse de la lumière ; n : l'ordre de l'interférence. On peut repérer les pics de Bragg par le transfert de quantité de mouvement défini par l'équation suivante : sin(9/ 2) n x = 2dhki L'intérêt d'exprimer les profils de diffusion (intensité mesurée) en fonction de x est dû au fait qu'un pic d'intensité peut être mesuré pour différents couples (À,13) mais pour une seule valeur de x (n fixe). A titre d'exemples, la figure 1 annexée montre les pics de Bragg de deux cristaux, à savoir le TNT (trinitrotoluène) et le sel (NaCI). Dans le cas de matériaux non cristallins, le spectre de diffusion à angle 25 faible est également représentatif du matériau examiné. A l'instar des interférences déterminées par la loi de Bragg pour un matériau cristallin, des phénomènes d'interférence peuvent également se produire entre les atomes et/ou molécules d'un matériau amorphe tel qu'un liquide, faisant cette fois intervenir une distribution continue de distances (fonction d'interférence 30 moléculaire, notée FIM). En effet, beaucoup de matériaux amorphes ont des arrangements réguliers sur des distances nanométriques (on parle d'ordre à courte distance). Ce type d'ordre est déterminé par des liaisons chimiques fortes pour les liaisons covalentes et ioniques. Cet ordre à courte distance provoque des hc interférences intramoléculaires et intermoléculaires. La figure 2 annexée illustre des exemples de fonctions d'interférence moléculaire, à savoir la FIM normalisée de l'eau (H20), la FIM normalisée de l'eau oxygénée (H202), la FIM normalisée de l'acétone(C2H6C0), et la FIM normalisée d'un matériau connu sous le nom commercial plexiglas® ((C502H8)n). Les diffractomètres les plus courants sont dits à dispersion angulaire ADXRD (acronyme anglais de « Anguier Dispersive X-ray Diffraction »). L'énergie est fixée par l'utilisation d'un rayonnement monochromatique et le nombre de photons diffractés est mesuré en fonction de l'angle. Bien que ces dispositifs soient très précis, ils nécessitent l'utilisation d'une source monochromatique puissante et ne peuvent être utilisés pour l'imagerie du fait de leur encombrement. Développée plus récemment, la technique EDXRD (acronyme anglais de « Energy Dispersive X-Ray Diffraction ») permet de pallier ces difficultés. Elle consiste cette fois à travailler à angle fixe en utilisant un ensemble de collimateurs et à illuminer l'objet avec un faisceau polychromatique pour mesurer les photons diffractés avec un détecteur spectrométrique résolu en énergie. Les pics de diffraction apparaissent alors à certaines énergies dans le spectre mesuré. La technique EDXRD, et de façon plus générale, toute technique d'analyse par spectrométrie, nécessite la mise en oeuvre d'un détecteur spectrométrique suffisamment résolu en énergie pour permettre la séparation et l'identification des différents pics caractéristiques du matériau constitutif de l'objet à analyser. Les détecteurs connus présentant la meilleure résolution en énergie sont du type Germanium. Mais ce type de détecteur doit être refroidi à des températures très faibles, par des procédés complexes et/ou coûteux (refroidissement thermoélectrique ou par un réservoir d'azote liquide). Aussi, les dispositifs d'analyse mettant en oeuvre un tel détecteur sont très encombrants. L'émergence récente de détecteurs spectrométriques pouvant être utilisés à température ambiante, tels des détecteurs de type CdTe, CdZnTe, ou des matériaux scintillateurs, offre une alternative intéressante aux détecteurs Germanium. En effet, ces détecteurs sont compacts, non refroidis et moins onéreux. Cependant, leurs performances en termes de résolution en énergie sont encore inférieures à celles obtenues avec les détecteurs Germanium, bien que correctes. Pour savoir si une substance cristalline ou amorphe donnée est contenue dans un objet, il est donc connu de : / irradier l'objet à l'aide d'un faisceau incident, émis par une source de rayonnement ionisant, de préférence collimaté par un collimateur primaire, / capter le rayonnement diffracté à l'aide d'un dispositif de détection comprenant . un détecteur, dit détecteur spectrométrique, apte à établir un spectre en énergie (ou spectre énergétique) du rayonnement diffusé selon un angle de diffusion donné, c'est-à-dire un détecteur comprenant o un matériau détecteur apte à interagir avec un rayonnement diffusé par l'objet et qui présente, du côté de l'objet, un plan dit plan de détection, o des moyens de mesure spectrométrique, aptes à mesurer une énergie libérée par chaque interaction d'un photon avec le matériau détecteur et à établir au moins un spectre en énergie. . un collimateur, dit collimateur de détection, associé au détecteur, le détecteur et le collimateur de détection étant agencés de façon à 15 présenter un axe de détection D formant, avec l'axe central Z du faisceau incident, un angle de diffusion 0, / analyser le ou les spectres mesurés par comparaison avec le spectre en énergie de la substance recherchée. D'une façon générale, un spectre en énergie illustre la distribution de 20 l'énergie d'un rayonnement sous la forme d'un histogramme représentant le nombre d'interactions de photon dans l'objet (en ordonnée) en fonction de l'énergie libérée (en abscisse). Souvent, l'axe des énergies est discrétisé en canaux de largeur 2 bE, un canal Ci, centré sur l'énergie Ei correspondant aux énergies comprises entre EiÔE et Ei-FbE. 25 On sait que l'identification d'un matériau par élaboration puis analyse d'un spectre de diffusion peut-être fastidieuse lorsqu'on souhaite caractériser des objets de volume important. En effet, cette méthode est d'autant plus précise que le volume d'objet observé est faible. L'invention vise à pallier ces inconvénients, en proposant notamment un 30 procédé d'analyse d'un objet plus rapide et plus efficace. Pour ce faire, l'invention propose un procédé d'analyse d'un objet caractérisé en ce qu'il se déroule en deux parties, à savoir une première partie comprenant les étapes suivantes : / irradiation de l'objet par un rayonnement photonique incident, / mesure d'un rayonnement transmis par l'objet, dit rayonnement en transmission mesuré, à l'aide d'un détecteur placé en transmission, / détermination d'au moins une première caractéristique de l'objet à partir du rayonnement en transmission mesuré, V vérification de la réalisation d'au moins un critère de suspicion portant sur la première caractéristique de l'objet et traduisant le fait que l'objet contient un matériau potentiellement suspect pour une application considérée ; ainsi par exemple, dans le cas d'une application relevant du domaine de la sécurité telle que la recherche d'explosifs dans des bagages, le critère de suspicion traduit le fait que l'objet contient un matériau potentiellement dangereux, c'est-à-dire un matériau dont la première caractéristique s'approche de celle d'un matériau explosif ; dans le cas d'une application médicale telle que la recherche d'une tumeur dans un sein, le critère de suspicion traduit le fait que l'objet contient un matériau potentiellement malin, c'est-à-dire un matériau dont la première caractéristique est proche de celle d'un tissu cancéreux. Le procédé selon l'invention comprend une seconde partie qui est exécutée uniquement lorsque le critère de suspicion est réalisé. Ainsi, lorsque le critère de suspicion n'est pas réalisé, il est mis fin au procédé d'analyse de l'objet à la fin de la première partie : l'objet en cours d'analyse 20 est considéré comme inoffensif, est déplacé et éventuellement remplacé (dans la zone de réception d'objet du système de détection) par un nouvel objet à analyser. A l'inverse, lorsque le critère de suspicion est réalisé, la seconde partie est exécutée, laquelle seconde partie comprend : / acquisition d'un spectre énergétique, dit spectre en diffusion mesuré, 25 diffusé par l'objet selon un angle de diffusion 0 non nul inférieur à 15° à l'aide d'un détecteur spectrométrique placé en diffusion, / détermination d'une seconde caractéristique de l'objet à partir -soit du spectre en diffusion mesuré, - soit du spectre en diffusion mesuré et du rayonnement en transmission mesuré, 30 V comparaison avec des caractéristiques de matériaux étalons mémorisées dans une base de données, soit de la seconde caractéristique de l'objet, soit non seulement de la seconde caractéristique de l'objet mais aussi de la première caractéristique de l'objet, aux fins d'identification d'un matériau constitutif de l'objet. En résumé, la définition du procédé fournie ci-dessus couvre quatre cas : - un premier cas dans lequel . la deuxième caractéristique (par exemple une signature en diffusion) est obtenue à partir du seul spectre en diffusion mesuré, le rayonnement en transmission ne servant alors que pour éliminer les objets non suspects (dans la première partie du procédé), . la deuxième caractéristique est utilisée seule pour identifier le matériau, la deuxième caractéristique pouvant alors être qualifiée de « signature déterminante » du matériau en ce sens qu'elle permet à elle seule d'identifier de façon fiable et rapide le matériau, comme par exemple une signature en diffusion (Pics de Bragg, Fonction d'interférence Moléculaire), - un deuxième cas dans lequel . la deuxième caractéristique (par exemple signature en diffusion) est obtenue à partir du spectre en diffusion mesuré et du rayonnement en transmission mesuré (qui dans ce cas est de préférence un spectre en transmission), le spectre en transmission mesuré servant alors à deux reprises, une première fois dans la première partie du procédé pour éliminer les objets non suspects et une deuxième fois dans la deuxième partie du procédé pour la détermination de la deuxième caractéristique par exemple pour améliorer la précision avec laquelle cette deuxième caractéristique est déterminée, . la deuxième caractéristique est utilisée seule pour identifier le matériau, la deuxième caractéristique étant alors une signature déterminante telle que définie supra, - un troisième cas dans lequel . la deuxième caractéristique (par exemple transferts de quantité de mouvement) est obtenue à partir du seul spectre en diffusion mesuré, le rayonnement en transmission ne servant alors que pour éliminer les objets non suspects (dans la première partie du procédé), . à la fois la deuxième caractéristique (par exemple transferts de quantité de mouvement) et la première caractéristique (par exemple numéro atomique effectif) sont utilisées pour identifier le matériau. - un quatrième cas dans lequel . la deuxième caractéristique (par exemple le spectre en diffusion lui-même, ou la position de certains pics de ce spectre en diffusion et éventuellement leur intensité relative) est obtenue à partir du spectre en diffusion mesuré et du rayonnement en transmission mesuré (qui dans ce cas est de préférence un spectre en transmission), le spectre en transmission mesuré servant alors à deux reprises, une première fois dans la première partie du procédé pour éliminer les objets non suspects et une deuxième fois dans la deuxième partie du procédé pour la détermination de la deuxième caractéristique par exemple pour améliorer la précision avec laquelle cette deuxième caractéristique est déterminée, . à la fois la deuxième caractéristique (par exemple le spectre en diffusion) et la première caractéristique (par exemple un coefficient d'atténuation, l'épaisseur d'un matériau constitutif de l'objet) sont utilisées pour identifier le matériau. A noter qu'il est possible d'utiliser un seul et même détecteur spectrométrique à la fois pour la mesure du rayonnement en transmission dans la première partie du procédé, lequel rayonnement mesuré est alors avantageusement un spectre énergétique en transmission, et pour l'acquisition du spectre en diffusion dans la deuxième partie du procédé. Il convient, dans ce cas, de déplacer le détecteur au cours du procédé, entre les deux positions d'acquisition (position en transmission, sur l'axe du rayonnement incident / position en diffusion, décalée d'un angle 0 de cet axe). En variante, on utilise un système de détection pourvu de deux détecteurs qui ne sont pas destinés à être déplacés, à savoir un détecteur placé en transmission, qui peut être un détecteur spectrométrique (version préférée) ou un détecteur autre comme par exemple un détecteur à intégration, et un détecteur spectrométrique placé en diffusion. Avantageusement et selon l'invention, lorsque le détecteur placé en transmission est un imageur, la première caractéristique peut correspondre à une forme géométrique de l'objet, ou à une de ces dimensions, comme par exemple son épaisseur suivant l'axe du rayonnement incident. Avantageusement et selon l'invention, le détecteur placé en transmission est un détecteur spectrométrique, ce qui rend possible l'estimation, en tant que première caractéristique de l'objet, outre celles citées ci-dessus, d'un coefficient d'atténuation pour une ou plusieurs énergies prédéterminées, ou pour une plage d'énergies donnée, ou encore du numéro atomique effectif Zeff d'un matériau constitutif de l'objet. Il peut également s'agir d'une indication quant à la nature d'un matériau constitutif de l'objet, comme par exemple : - une famille de matériaux à laquelle appartiendrait le matériau constitutif de l'objet, la seconde partie du procédé permettant alors de préciser au sein de cette famille 5 quelle est la nature exacte du matériau et donc d'affiner l'hypothèse fournie par la première partie, - une hypothèse sur la nature-même du matériau, la seconde partie du procédé permettant alors de confirmer l'hypothèse fournie par la première étape quant à la nature du matériau, 10 ces indications étant obtenues à partir de l'analyse du spectre en transmission. Le numéro atomique effectif Zeff fait partie des paramètres pouvant caractériser un matériau. On rappelle que le numéro atomique d'un corps simple est défini comme le nombre de protons présents dans le noyau d'un atome de ce corps. En revanche, lorsque l'on considère un composé chimique, on doit faire appel à la 15 notion de numéro atomique effectif. Celui-ci est défini comme le numéro atomique d'un corps simple qui conduirait au même spectre de transmission dans une bande d'énergie donnée. De manière générale, le numéro atomique effectif d'un composé chimique est obtenu au moyen d'une combinaison des numéros atomiques des corps simples constitutifs du composé, chaque numéro atomique étant affecté d'un 20 coefficient de pondération dépendant de la fraction massique ou atomique du corps simple dans le composé. Ainsi, en pratique, le numéro atomique effectif eff d'un Zeff = (Zieff y composé de N corps simples, vérifie où zieff est le numéro = atomique du corps simple i1,...,Net P est une constante relié à l'effet , photoélectrique (p- 462) 25 Le critère de suspicion définit de préférence une plage de valeurs de la première caractéristique à l'intérieur de laquelle il est considéré que l'objet est potentiellement suspect eu égard à l'application considérée, et à l'extérieur de laquelle l'objet est considéré comme n'étant pas suspect. L'étape de vérification de la réalisation du critère de suspicion comprend alors une comparaison entre la 30 première caractéristique déterminée pour l'objet et la plage de valeurs prédéfinie eu égard à l'application considérée. JI/p Un avantage à déterminer la première caractéristique à partir d'un spectre en transmission mesuré est qu'un tel spectre nécessite un temps d'acquisition relativement court, de sorte que la première partie du procédé qui ne permet pas d'identifier précisément le matériau mais permet d'éliminer un grand nombre d'objets inoffensifs est rapide. La deuxième partie du procédé, qui permet une meilleure discrimination du matériau mais qui est plus longue en raison de l'acquisition d'un spectre en diffusion, n'est effectuée que sur les objets suspects. Un gain de temps considérable est ainsi réalisé. La deuxième caractéristique de l'objet, déterminée à partir au moins d'un spectre en diffusion mesuré, est par exemple une fonction appelée signature en diffusion, qui représente des pics de diffusion (pics de Bragg pour un matériau cristallin ou fonction d'interférence moléculaire pour un matériau amorphe). En variante, il peut s'agir d'une grandeur physique extraite d'une telle signature en diffusion (distance inter-réticulaire pour un matériau cristallin par exemple). La deuxième caractéristique de l'objet peut également être un transfert de quantité de mouvement. Comme défini précédemment, cette deuxième caractéristique est déterminée à partir au moins du spectre en diffusion mesuré. Dans un mode de réalisation préférentiel, elle est déterminée à partir non seulement du spectre en diffusion mesuré lors de la deuxième partie du procédé mais aussi du rayonnement en transmission mesuré lors de la première partie procédé. L'utilisation combinée d'un spectre en diffusion et d'un spectre en transmission pour reconstruire une signature en diffusion d'un matériau est décrite plus en détail plus loin dans le cadre d'un exemple de réalisation de l'invention. Cette 25 combinaison permet d'améliorer l'identification du matériau constitutif de l'objet. Dans une version préférée, le détecteur spectrométrique placé en diffusion est configuré de façon à présenter un axe de détection formant, avec un axe central du rayonnement incident, un angle de diffusion 0 Inférieur à 10° et de préférence inférieur à 5°. 30 Dans un mode de réalisation préférentiel, le procédé d'analyse selon l'invention comprend une première étape préalable de calibrage d'une matrice de réponse du détecteur spectrométrique placé en diffusion (laquelle étape de calibrage peut être omise notamment si la réponse du détecteur est jugée satisfaisante, c'est-à-dire si la résolution en énergie du détecteur utilisé est particulièrement bonne), une deuxième étape préalable de calibrage d'une matrice de réponse du détecteur placé en transmission lorsque celui-ci est un détecteur spectrométrique (laquelle étape de calibrage peut être omise notamment si la réponse du détecteur est jugée satisfaisante, c'est-à-dire si la résolution en énergie du détecteur utilisé est particulièrement bonne), et une troisième étape de calibrage d'une matrice de réponse angulaire du système de détection, ces première, deuxième et troisième étapes préalables de calibrage pouvant être effectuées dans un ordre ou dans un autre, avant toute irradiation d'un objet à analyser, c'est-à-dire « off-line ». L'invention s'étend à un système de détection apte à mettre en oeuvre le 10 procédé selon l'invention, et en particulier à un système de détection pour l'analyse d'un objet comprenant : V une source de rayonnement photonique ; V une zone de réception d'un objet à analyser ; Y des moyens de détection, en aval de la zone de réception d'objet, pour la 15 mesure d'un rayonnement transmis par l'objet, dit spectre en transmission mesuré et pour l'acquisition d'un spectre énergétique, dit spectre en diffusion mesuré, diffusé par l'objet selon un angle de diffusion 0 non nul inférieur à 15° ; ces moyens de détection comprenant un détecteur spectrométrique ; Y des moyens informatiques de traitement comprenant 20 . des moyens de détermination d'au moins une première caractéristique de l'objet à partir du rayonnement en transmission mesuré, la première caractéristique pouvant être une forme, une dimension de l'objet, un coefficient d'atténuation pour une ou plusieurs énergies prédéterminées, ou pour une plage d'énergies prédéterminée, ou un numéro atomique effectif 25 Zeff d'un matériau constitutif de l'objet, ou une indication quant à la nature d'un matériau constitutif de l'objet (première estimation de la nature du matériau -à confirmer à l'aide du spectre de diffusion dans la seconde partie du procédé- ou identification d'une famille ou d'un type de matériaux), et éventuellement son épaisseur. 30 . des moyens de vérification d'au moins un critère de suspicion portant sur la première caractéristique précédemment déterminée pour l'objet et traduisant le fait que l'objet contient un matériau potentiellement suspect pour une application considérée, par exemple par comparaison de la première caractéristique avec une plage de valeurs prédéfinie eu égard à l'application concernée ; . des moyens de détermination d'une deuxième caractéristique de l'objet, à partir soit du spectre en diffusion mesuré, soit à la fois du spectre en diffusion mesuré et du rayonnement en transmission mesuré, et des moyens de comparaison à des caractéristiques de matériaux étalons mémorisées dans une base de données, soit de la deuxième caractéristiques, soit des première et deuxième caractéristiques, aux fins d'identification d'un matériau constitutif de l'objet. Selon un premier mode de réalisation possible, les moyens de détection comprennent, outre le détecteur spectrométrique précédemment défini, des moyens de positionnement pour le positionnement successif dudit détecteur spectrométrique en une position de transmission pour la mesure du rayonnement en transmission mesuré, et en une position de diffusion pour l'acquisition du spectre en diffusion mesuré. Selon un deuxième mode de réalisation possible, les moyens de détection comprennent, outre le détecteur spectrométrique précédemment défini, qui est alors placé en diffusion, un autre détecteur, placé en transmission, lequel détecteur peut être un détecteur spectrométrique ou un détecteur autre tel un détecteur à intégration. Le système de détection selon l'invention présente de plus, de préférence, 20 une ou plusieurs des caractéristiques suivantes : / la première caractéristique est choisie parmi une forme de l'objet, une dimension de l'objet, un coefficient d'atténuation à une ou plusieurs énergies prédéterminées d'un matériau constitutif de l'objet, un numéro atomique effectif Zeff d'un matériau constitutif de l'objet., une indication quant à la nature ou 25 éventuellement une épaisseur d'un matériau constitutif de l'objet, / la deuxième caractéristique est choisie parmi une signature de diffusion représentant des pics de Bragg ou une fonction d'interférence moléculaire, des transferts de quantité de mouvement. / Les moyens informatiques de traitement sont aptes à mettre en oeuvre 30 une méthode basée sur une approche de type problème inverse et à estimer une matrice de réponse globale du système de détection. / Les moyens informatiques de traitement sont aptes à estimer un spectre incident atténué par l'objet en utilisant le spectre en transmission mesuré. / Les moyens informatiques de traitement sont aptes à construire une matrice de réponse globale du système de détection à partir d'un spectre incident atténué estimé, d'une matrice de résolution angulaire calibrée du système de détection ainsi que, éventuellement, d'une matrice de réponse calibrée du détecteur spectrométrique placé en diffusion, d'une matrice de réponse calibrée du détecteur placé en transmission. / La source de rayonnement est polychromatique. / La source de rayonnement est une source de rayons X. / Le détecteur spectrométrique placé en diffusion est configuré de façon à présenter un axe de détection formant, avec un axe central du rayonnement incident, 10 un angle de diffusion 0 Inférieur à 10° et de préférence inférieur à 5°. / Le système de détection comprend un premier collimateur, dit collimateur de source, placé entre la source de rayonnement photonique et la zone de réception d'objet, de façon à collimater le rayonnement émis par la source en un faisceau incident d'axe central Z, 15 V Le système de détection comprend un deuxième collimateur, dit collimateur de diffusion, entre la zone de réception d'objet et le détecteur spectrométrique placé en diffusion, le collimateur de diffusion ayant un axe central de collimation formant un angle de diffusion 0 avec l'axe central incident Z, lequel angle de diffusion 0 est de préférence inférieur à 5°. 20 v Les détecteurs spectrométriques utilisés contiennent des capteurs à conversion directe, à savoir que les photons X incidents sur le capteur sont absorbés par un semi-conducteur (CdTe par exemple) et créent un nuage de charges électroniques (typiquement 10000 électrons pour un photon X de 60 keV). Ces charges sont ensuite collectées par des électrodes et forment un signal électrique 25 transitoire appelé impulsion. Si la collecte est complète, l'intégrale de l'impulsion mesurée pour chaque particule détectée est proportionnelle à son énergie déposée dans le semi-conducteur. Un circuit électronique permet de mesurer cette intégrale. Après numérisation, les différentes mesures de photons de même énergie sont sommées, ce qui permet de reconstruire le spectre des rayonnements X ayant 30 interagi avec l'objet irradié. L'invention s'étend à un système de détection et à un procédé d'analyse d'un objet caractérisés en combinaison par tout ou partie des caractéristiques mentionnées ci-avant et ci-après. D'autres détails et avantages de la présente invention apparaîtront à la lecture de la description suivante, qui se réfère aux dessins schématiques annexés et porte sur des modes de réalisation préférentiels, fournis à titre d'exemples non limitatifs. Sur ces dessins : / La figure 1 est un graphique représentant les pics de Bragg du TNT (trinitrotoluène) et du sel (NaCI), lesquels pics illustrent l'intensité normalisée (axe des ordonnées), c'est-à-dire le nombre relatif de photons détectés au cours d'une opération d'acquisition, en fonction du transfert de quantité de mouvement x (axe des abscisses) en nm-ldes photons détectés. / La figure 2 est un graphique représentant la fonction d'interférence moléculaire normalisée de l'eau (H20), la fonction d'interférence moléculaire normalisée de l'eau oxygénée (H202), la fonction d'interférence moléculaire normalisée de l'acétone (C2H6C0), et la fonction d'interférence moléculaire normalisée du plexiglas® ((C502H8)n), avec, en ordonnées, le nombre de photons détectés et, en abscisses, le transfert de quantité de mouvement x en nm-1. v La figure 3 est une vue schématique d'un système de détection selon l'invention. / La figure 4 est un graphique représentant des spectres en transmission pour des échantillons cylindriques de 40mm de diamètre d'eau (H20) et d'acétone (C2H6C0), tels que mesurés par le détecteur spectrométrique placé en transmission 20 du système de détection de la figure 3. / La figure 5 est un graphique représentant des spectres en diffusion pour des échantillons cylindriques de 40mm de diamètre d'eau (H20) et d'acétone (C2H6C0), tels que mesurés par le détecteur spectrométrique placé en diffusion du système de détection de la figure 3. 25 v La figure 6 est un graphique représentant la signature de diffusion (FIM estimée) d'un échantillon cylindrique d'eau de 40 mm de diamètre et celle d'un échantillon cylindrique d'acétone de 40 mm de diamètre, telles que reconstruites à partir d'un procédé et d'un système de détection selon l'invention, / La figure 7 est un graphique représentant d'une part la signature de 30 diffusion (« FIM reconstruite ») d'un échantillon cylindrique d'eau de 40 mm de diamètre et celle d'un échantillon cylindrique d'acétone de 40 mm de diamètre, telles que reconstruites à partir d'un procédé et d'un système de détection selon l'invention, et d'autres part les signatures de diffusion théoriques (« FIM théorique ») de ces échantillons, / La figure 8 est un graphique représentant une matrice de réponse du détecteur spectrométrique placé en transmission du système de détection de la figure 3. / La figure 9 est un graphique représentant une matrice de réponse du 5 détecteur spectrométrique placé en diffusion du système de détection de la figure 3. / La figure 10 est un graphique représentant la distribution angulaire du détecteur spectrométrique placé en diffusion du système de détection de la figure 3. / La figure 11 est un graphique représentant une matrice de réponse angulaire du système de détection de la figure 3. 10 V la figure 12 est un graphique représentant des spectres incidents atténués estimés pour des échantillons cylindriques de 40mm de diamètre d'acétone et d'eau, tels qu'estimés selon l'invention à partir des spectres en transmission mesurés de la figure 4. / La figure 13 est un digramme illustrant une opération de construction 15 d'une matrice de réponse globale du système de détection de la figure 3. / la figure 14 est un graphique représentant la matrice de réponse globale A du système de détection de la figure 3 pour un échantillon cylindrique d'eau de 40 mm de diamètre. / la figure 15 est un graphique représentant la matrice de réponse globale 20 A du système de détection de la figure 3 pour un échantillon cylindrique d'acétone de 40 mm de diamètre. Le système de détection selon l'invention illustré à la figure 3 comprend : . une source 1 polychromatique de rayonnements ionisants, telle qu'un tube à rayons X, 25 . un collimateur de source 2, qui permet de canaliser les rayonnements de la source 1 en un faisceau incident d'axe central Z, - une zone 4 de réception d'un objet à analyser, . un collimateur de diffusion 5 ayant un axe de collimation D, . un détecteur spectrométrique 6 placé en diffusion, associé au collimateur de 30 diffusion 5 de sorte que le détecteur 6 détecte un rayonnement diffusé selon un angle de diffusion 0 (angle entre l'axe incident Z et l'axe de collimation et de détection D) par exemple égal à 2,5° (l'échelle n'étant pas respectée sur la figure 3) ; le détecteur spectrométrique 6 placé en diffusion est apte à établir un spectre (énergétique) en diffusion mesuré, c'est-à-dire un spectre en énergie du rayonnement diffusé par l'objet selon la direction D ; le détecteur spectrométrique 6 est, en l'exemple, un détecteur à matériau semi-conducteur, tel un détecteur à CdTe ou CdZnTe ; - un détecteur spectrométrique 7 placé en transmission, apte à établir un spectre (énergétique) en transmission mesuré, c'est-à-dire un spectre en énergie du rayonnement transmis par l'objet selon la direction Z; le détecteur spectrométrique 7 placé en transmission est, en l'exemple, un détecteur à matériau semi-conducteur, tel un détecteur à CdTe ou CdZnTe ; - des moyens informatiques 8 de traitement des spectres mesurés fournis par les détecteurs spectrométriques 6 et 7. Les termes « rayonnement transmis » désignent le rayonnement constitué de photons n'ayant subi aucune interaction dans l'objet examiné. Par « spectre en transmission », on entend le spectre du rayonnement transmis selon l'axe du faisceau incident à l'objet, constitué par les photons qui n'ont subi aucune interaction dans l'objet. L'expression « placé en transmission » désigne un détecteur apte à détecter le rayonnement transmis par le matériau. Ainsi, un détecteur placé en transmission est situé sur l'axe du rayonnement incident à l'objet, l'objet étant placé entre le détecteur et la source de rayonnement. Par « détecteur spectrométrique », on entend un détecteur apte à générer 20 un spectre énergétique du rayonnement détecté. Le procédé selon l'invention vise à déterminer si un objet à analyser est dangereux (explosif, tumeur...) ou non. Pour ce faire, le procédé se déroule en deux parties. L'objectif de la première partie du procédé n'est pas l'identification précise du matériau constitutif de l'objet ; cette première partie vise seulement à déterminer si 25 l'objet est suspect ou non, c'est-à-dire à déterminer si l'objet mérite qu'une analyse plus poussée soit menée. Cette première partie permet ainsi d'écourter l'analyse pour un grand nombre d'objets qui s'avèrent manifestement inoffensifs. Le procédé décrit ici a été mis en oeuvre par les inventeurs pour étudier deux échantillons cylindriques de 40mm de diamètre constitués respectivement 30 d'acétone (C2H6C0) et d'eau (H20). Pour chaque échantillon inspecté, la première partie du procédé exploite un spectre en transmission mesuré tel que fourni par le détecteur spectrométrique 7, et dont le modèle direct est le suivant : h = RE, - (S,' x An) Avec - h: le vecteur du spectre en transmission mesuré de taille (NbEit x 1) - REt : la matrice de réponse du détecteur spectrométrique placé en transmission, de taille (NbEit x NbE,). Dans le cas d'un détecteur parfait, cette matrice est une matrice diagonale. Chaque terme REt (j,i) de la matrice représente la probabilité de détection d'une valeur d'énergie égale à j lorsque le photon incident au détecteur a une énergie i. - : le vecteur du spectre incident du tube 1 à rayons X de taille (1 x NbEi) - Att : le vecteur d'atténuation de taille (1 x NbE,) qui prend en compte les effets d'atténuation dans l'objet. - Le vecteur (Sinc x Att) représente le spectre de la source de rayonnement atténuée par l'objet. Dans l'exemple décrit ici, et qui fournit des signatures très précises, ce vecteur est pris en compte dans la construction d'une matrice de réponse globale A du système, comme on le comprendra plus loin. - NbEit et NbE, correspondent respectivement au nombre de canaux du spectre en transmission mesuré (c'est-à-dire au nombre de canaux du spectre de l'énergie détectée par le détecteur placé en transmission) et au nombre de canaux du spectre de l'énergie incidente. - Le symbole x correspond à un produit terme à terme (S,n, et Att sont multipliés terme à terme, on obtient alors un vecteur qui fait la même taille). - Le symbole . correspond au produit matriciel classique. Les spectres en transmission mesurés avec le détecteur spectrométrique 7 pour les échantillons d'acétone (C2H6C0) et d'eau (H20) susmentionnés peuvent 25 être observés à la figure 4. Plusieurs caractéristiques peuvent être extraites d'un tel spectre en transmission : coefficient d'atténuation, numéro atomique effectif Zeff d'un matériau constitutif de l'objet, une estimation de la nature, voire de la nature et de l'épaisseur, d'un matériau constitutif de l'objet. 30 Le numéro atomique effectif Zeff du matériau peut par exemple être obtenu par une méthode statistique qui compare le spectre en transmission mesuré pour l'objet à des spectres d'échantillons de calibrage (matériaux étalons) mesurés dans les mêmes conditions expérimentales, et qui utilise un modèle physique d'évolution du spectre en fonction du numéro atomique effectif et de l'épaisseur du matériau. Cette technique est décrite dans la demande de brevet français N°1363174 d'Alexia Gorecki au nom du CEA, incorporée ici par référence. Le spectre en transmission peut également être mis en oeuvre pour estimer la nature du matériau examiné. L'estimation de la nature d'un matériau, ainsi qu'éventuellement son épaisseur, à partir d'un spectre en transmission est décrite dans les demandes W02012000993 ou W02011069748, au nom du CEA, incorporées ici par référence. Une fois cette ou ces caractéristiques de l'objet estimées, une décision est prise. Si l'objet inspecté semble suspect, c'est-à-dire si un ou plusieurs critères de suspicion sont réalisés (forme douteuse, sous réserve d'utiliser un imageur en tant que détecteur placé en transmission, ou numéro atomique appartenant à une plage de valeurs critique de matériaux potentiellement dangereux définie au préalable en fonction de l'application concernée), alors un spectre en diffusion du matériau sera acquis et étudié pour permettre une meilleure discrimination du matériau. Cette seconde partie n'est pas faite de façon systématique sur tous les matériaux car, bien que très discriminante, elle nécessite un temps d'exposition aux rayons X relativement long. S'il n'y a aucun doute sur l'innocuité de l'objet, alors la première partie du procédé est menée sur un nouvel objet. Le numéro atomique effectif des deux échantillons inspectés (eau et acétone) a été évalué à partir des spectres de la figure 4. Les résultats suivants ont été obtenus : V pour l'échantillon d'eau : Zeff estimé = 7.21. En théorie, un tel échantillon 25 possède un numéro atomique effectif : Zeff théorique = 7.22. On constate donc que le résultat obtenu est très proche de la valeur théorique de l'échantillon. V Pour l'échantillon d'acétone : Zeff estimé = 6.05, pour un numéro atomique théorique Zeff théorique = 6.03. Une zone critique a été définie pour cette application (détection 30 d'explosifs) ; elle correspond à des matériaux dont le Zeff est entre 6 et 8.5. En d'autres termes, l'étape de vérification de la réalisation du critère de suspicion consiste à comparer le numéro atomique effectif déterminé, estimé à partir du spectre en transmission mesuré, à la plage de valeurs [6 ;8,5]. A noter que plusieurs zones critiques discontinues (plusieurs plages de valeurs) peuvent être définies pour une même première caractéristique (numéro atomique par ex.). Il est également possible de définir plusieurs critères de suspicion portant sur des caractéristiques distinctes ; dans ce cas, on peut considérer soit que l'objet est suspect (et donc éligible pour la seconde partie du procédé) dès lors qu'au moins l'un de ces critères de suspicion est réalisé, soit que l'objet n'est suspect que si tous ces critères de suspicion sont réalisés. Les deux échantillons inspectés ici (eau et acétone) ont des Zeff qui sont dans la zone critique ; ils sont donc considérés comme suspects et la deuxième partie du procédé selon l'invention, à savoir une étude de ces deux échantillons en 10 diffusion petit angle, doit donc être exécutée. La deuxième partie du procédé vise à identifier précisément le matériau constitutif de l'objet lorsque celui-ci a été considéré comme suspect à l'issue de la première partie. Cette deuxième partie consiste par exemple à fournir une signature en diffusion dudit matériau à partir d'un spectre en diffusion mesuré (fourni par le 15 détecteur spectrométrique 6), les termes « signature en diffusion » désignant une fonction représentative du matériau et montrant des pics de diffusion (pics de Bragg ou Fonction d'Interférence Moléculaire selon la nature -cristalline ou amorphe- du matériau). Les spectres en diffusion mesurés par le détecteur 6 pour les deux 20 échantillons peuvent être observés à la figure 5. Dans l'exemple décrit ici, il a été choisi d'exploiter, de façon originale, les spectres en diffusion mesurés en combinaison avec les spectres en transmission mesurés au cours de la première partie du procédé. En variante, il est possible d'exploiter les spectres en diffusion mesurés de la figure 5 selon une autre technique telle une technique connue utilisée 25 en EDXRD. Dans l'exemple choisi ici, la deuxième partie du procédé utilise le modèle suivant, décrivant la relation entre la signature en diffusion fdu matériau et le spectre en diffusion mesuré g: g = (RF(] x S11. X Au). Re -f=i4.f 30 Avec - g: le vecteur du spectre en diffusion (cohérente) mesuré, de taille (NbEid x 1) - REd : matrice de réponse du détecteur spectrométrique placé en diffusion, de taille (Nbqd x NbE,). Dans le cas d'un détecteur parfait, cette matrice est une - - - 25 - 30 matrice diagonale. Chaque terme REd (j,i) de la matrice représente la probabilité de détection d'une valeur d'énergie égale à j sachant que le rayonnement incident au détecteur a une énergie égale à i. D'une façon générale, la matrice de réponse d'un détecteur spectrométrique établit une relation probabiliste entre une énergie détectée par le détecteur et l'énergie du rayonnement incident à ce détecteur. Chaque colonne i de REd (j,i), avec j variant de 1 à NbEjd, correspond à la densité de probabilité d'énergie détectée par le détecteur lorsque ce dernier est soumis à un rayonnement incident d'énergie i. : vecteur du spectre incident du tube à rayons X de taille (1 x NbE,) Att: vecteur d'atténuation de taille (1 x NbE,) qui prend en compte les effets d'atténuation dans l'objet. Du fait de la faible valeur de l'angle 0, inférieur à 15° et de préférence inférieur à 10° voire 5°,on fera l'approximation que l'atténuation par l'objet le long du trajet vers les deux détecteurs (en diffusion et en transmission) est la même. R0: matrice de réponse angulaire du système de détection, de taille (N bE, x Nbx). Chaque terme Ro (j,k) de la matrice Ro correspond à une probabilité que l'énergie d'un photon détecté à l'énergie j corresponde à un transfert de quantité de mouvement égal à k . Autrement dit, Ro (j, k) correspond à la probabilité qu'un transfert de quantité de mouvement k donne lieu à la détection d'un photon d'énergie j. Chaque colonne k Ro(j,k), avec j variant de 1 à NbEjd, correspond à la densité de probabilité d'énergie détectée par le détecteur lors d'un transfert de quantité de mouvement égal à k. D'une façon plus générale, la matrice de réponse angulaire Ro permet d'établir une relation probabiliste entre l'énergie détectée par le détecteur placé en diffusion et un paramètre de diffusion élastique d'un matériau constituant l'objet, en particulier le transfert de quantité de mouvement. f: signature en diffusion, de taille (1 x Nbx), propre au matériau constituant l'objet, et qui permet de décrire soit les pics de Bragg théoriques du matériau dans le cas d'un matériau cristallin, soit la fonction d'interférence moléculaire dans le cas d'un matériau amorphe, - A: matrice de réponse globale du système en diffusion, de taille (NbEid x Nbx). Chaque terme A(j,k) de A correspond à une probabilité que l'énergie d'un photon détecté, par le détecteur en diffusion, à l'énergie j corresponde à un transfert de quantité de mouvement égal à k . Autrement dit, A(j,k) correspond à la probabilité qu'un transfert de quantité de mouvement k donne lieu à la détection d'un photon à l'énergie j. - Le symbole x correspond à un produit terme à terme (S1 et Att sont multipliés terme à terme, on obtient alors un vecteur qui fait la même taille). - Le symbole . correspond au produit matriciel classique. - NbEid, NbE, et Nb x correspondent respectivement au nombre de canaux du spectre en diffusion mesuré (c'est-à-dire au nombre de canaux du spectre de l'énergie détectée par le détecteur placé en diffusion), au nombre de canaux du spectre de l'énergie incidente et au nombre de canaux du vecteur décrivant le transfert de quantité de mouvement. A noter que le nombre de photons détectés dans chaque canal du vecteur g suit une loi de Poisson de paramètre le nombre moyen de photons dans ce canal. Le procédé selon l'invention comprend alors une opération de construction d'une matrice de réponse globale A du système de détection, utilisant le modèle ci-dessus. Pour ce faire, il convient au préalable de déterminer les termes Ro et (Sinc x 20 Att), et éventuellement REd.. Chacune de ces étapes est individuellement décrite plus loin. L'étape préalable de calibrage de la matrice de réponse REd du détecteur spectrométrique placé en diffusion n'est pas nécessaire mais avantageuse, car elle prend en compte la dégradation des spectres du fait de la réponse du détecteur. 25 Cependant, cette étape est optionnelle voire inutile, en particulier pour des détecteurs suffisamment résolus en énergie et lorsque la réponse d'un détecteur est jugée satisfaisante. Il en va de même pour la matrice de réponse REt du détecteur spectrométrique placé en transmission, utilisée plus loin. Une fois la matrice de réponse globale A construite à l'aide du modèle 30 susmentionnée, le procédé selon l'invention reconstruit la signature f (fonction d'interférence moléculaire pour les amorphes, distribution des dhk, pour les polycristallins) à partir du modèle g = A f (où A et g sont alors connus) en mettant en oeuvre une méthode basée sur une approche de type problème inverse. L'algorithme Maximum Likelihood - Expectation Maximization (noté MLEM) se propose d'estimer le spectre à calculer par maximisation itérative de la fonction de log-vraisemblance. Ce type de calcul est très fréquent dès lors qu'il s'agit 5 d'estimer un maximum de vraisemblance, et fait appel à un algorithme plus général, appelé Expectation - Maximization (EM). Cette méthode a l'avantage de prendre en compte la nature Poissonienne des données mesurées. Les coefficients de la matrice de réponse globale A du système sont notés aii. On veut maximiser la probabilité que l'estimée f de dimension Nbx génère les 10 mesures g. On sait de plus que les données mesurées obéissent à une statistique de Poisson, de par leur nature physique. On peut donc écrire la fonction de vraisemblance de l'estiméef : Sa log-vraisemblance s'exprime alors: cp(f) Io 15 On cherche ensuite à maximiser cette fonction, en annulant sa dérivée açb(f) f k a f k La résolution itérative de ce problème s'écrit alors : 20 En initialisant le vecteur I'°~k avec des valeurs positives, on a la garantie d'avoir des résultats non-négatifs. Ainsi, à partir d'une estimation de A et de la mesure de g, on peut reconstruire f en itérant l'algorithme MLEM. Les résultats de cette reconstruction (notés « FIM reconstruite » sur la fig.6) pour les échantillons d'eau et d'acétone 25 inspectés, en mettant en oeuvre un algorithme MLEM avec 100 itérations, sont illustrés à la figure 6. En d'autres termes, à partir de mesures réalisées en transmission et en diffusion sur un objet inconnu, on peut reconstruire une fonction relative à la structure Nbx ici=1"j,la j j,k d'un matériau constitutif de l'objet. Les valeurs de cette fonction sont représentées dans la matrice A. Ce matériau étant inconnu, l'objectif est alors de l'identifier. Pour ce faire, on utilise un ensemble de matériaux étalons (de type explosifs et non explosifs dans 5 le cas d'une application d'analyse de bagages par exemple ; de type tissus biologiques sains et tissus tumoraux dans le cas d'une application d'analyse médicale) dont les numéros atomiques effectifs ou autre première caractéristique et dont les signatures en diffusion (fonction d'interférence moléculaire ou pics de Bragg) ou autre seconde caractéristique sont tabulées et mémorisées dans une base de 10 données, et le procédé d'analyse selon l'invention consiste ensuite à comparer les valeurs obtenues pour l'objet à analyser avec celles de la base de données, pour identifier l'objet inconnu. En variante, certains paramètres permettant de remonter à des paramètres structurels du matériau sont extraits de la signature reconstruite pour 15 l'objet ; par exemple, dans le cas d'un matériau cristallin, l'extraction de la position des pics présents dans la signature obtenue permet de remonter aux distances inter-réticulaires du cristal. En comparant d'une part les valeurs de Zeff estimés à une base de valeurs de Zeff théoriques pour un ensemble de matériau et d'autre part la forme des 20 fonctions d'interférence moléculaire reconstruites à une base de forme de fonctions d'interférence moléculaire pour un ensemble de matériaux, il est possible d'identifier la nature des deux échantillons analysés (eau et acétone). La figure 7 permet de comparer les fonctions d'interférence moléculaire reconstruites des deux échantillons analysés (eau et acétone) et leurs fonctions d'interférence moléculaire théoriques. 25 Sont maintenant décrites les diverses étapes de l'opération de construction de la matrice de réponse globale A du système. Préalablement à toute analyse d'un objet, c'est-à-dire « off-line », des opérations de calibrage sont réalisées pour déterminer certaines spécifications du système de détection, qui dépendent notamment des détecteurs utilisés et de la 30 géométrie du système, et qui, contrairement aux vecteurs d'atténuation, ne dépendent pas de l'objet à analyser. Ces spécifications sont REt, REd, Ro. Elles sont ensuite mémorisées dans les moyens informatiques 8 de traitement. La matrice de réponse REt du détecteur spectrométrique 7 placé en transmission peut être obtenue à partir du logiciel de simulation Monte-Carlo Tasmania, qui permet de simuler toute la chaîne de détection d'un détecteur semiconducteur (interactions photon, transit des porteurs de charges...). De préférence, cette simulation est de plus confrontée avec des données expérimentales acquises par exemple avec des sources gamma. Ceci permet d'ajuster la résolution en énergie obtenue en simulation. La figure 8 montre la matrice de réponse REt calibrée pour le détecteur spectrométrique 7 placé en transmission. Cette matrice définit la probabilité de détecter un photon à l'énergie Ej alors que l'énergie incidente du photon est Ei. Cette probabilité est indiquée sur la figure 8 par une nuance de gris (dont l'échelle est reportée à droite du graphique), l'axe des abscisses du graphique représentant l'énergie incidente Ei exprimée en keV, l'axe des ordonnées correspondant à l'énergie détectée Ej exprimée en keV. Dans le cas d'un détecteur parfait, la matrice est diagonale (si elle est carrée). De façon similaire, une opération préalable de calibrage d'une matrice de 15 réponse REd du détecteur spectrométrique 6 placé en diffusion est exécutée « offline », par simulation à l'aide du logiciel de simulation Monte-Carlo et/ou par expérimentation. La matrice de réponse calibrée REd obtenue est illustrée à la figure 9. Là encore, l'axe des abscisses représente l'énergie incidente Ei en keV, l'axe des 20 ordonnées correspond à l'énergie détectée Ej en keV, la probabilité du couple (Ei, Ej) étant traduite par une nuance de gris. Est aussi exécutée « off-line », une opération préalable de calibrage d'une matrice de réponse angulaire Ro du système de détection. Cette réponse angulaire dépend de deux facteurs. Le premier correspond à la géométrie du système 25 d'acquisition et plus précisément à l'ouverture du collimateur de source 2 et à l'ouverture du collimateur de diffusion 5, sachant qu'on suppose que l'objet remplit l'intersection de deux cônes, à savoir un cône d'irradiation et un cône d'observation. Le cône d'irradiation délimite l'angle solide sous lequel l'objet voit la source, tandis que le cône d'observation correspond à l'angle solide sous lequel le détecteur voit 30 l'objet. On évalue tout d'abord une distribution angulaire 1D du système, soit à partir de simulations, soit à partir de calibrages. Utilisant la relation liant x (transfert de quantité de mouvement), E (énergie incidente au détecteur placé en diffusion) et 0 (angle de diffusion), on déduit la matrice de réponse angulaire fonction de Ei (Energie incidente) et de x à partir de la distribution angulaire 1D fonction de 0. La figure 10 montre un exemple de distribution angulaire 1D du système de diffusion lorsque l'axe de collimation D du collimateur de diffusion définit un angle de diffusion 0 égal à 2,5°, la distribution angulaire exprimant la quantité relative de photons incidents au détecteur placé en diffusion (axe des ordonnées) en fonction de 5 l'angle de détection en degré (axe des abscisses). Il s'agit d'un exemple car l'utilisation de collimateurs présentant des configurations différentes (en particulier largeur et longueur de l'ouverture des collimateurs) conduirait à l'obtention d'un autre graphique. Une matrice de réponse angulaire Ro du système de détection peut être observée à la figure 11, dont l'axe des abscisses représente le transfert de quantité 10 de mouvement x en nm-1, tandis que l'axe des ordonnées représente l'énergie E incidente au détecteur placé en diffusion en keV, des nuances de gris traduisant la quantité relative de photons incidents. Comme précédemment évoqué, cette matrice définit une relation probabiliste entre le nombre de photons incidents au détecteur placé en diffusion, à une énergie donnée, et le transfert de quantité de mouvement. 15 La construction de la matrice de réponse globale A du système de détection utilisant le modèle A =(REd x S x An). Ro nécessite encore une étape d'estimation d'un spectre incident atténué par l'objet (Sin, x Att). Avantageusement, cette étape d'estimation du spectre incident atténué utilise à nouveau le spectre en transmission mesuré lors de la première partie du 20 procédé par le détecteur spectrométrique 7 placé en transmission. Ce spectre en transmission h peut s'écrire : h = RE, - (s111 x Att) En d'autres termes, il est considéré dans cette version de l'invention que le terme (S,' x Att) dans l'expression de la matrice A est égal au terme (S,' x Att) 25 dans l'expression du spectre en transmission h. Les inventeurs ont démontré que cette approximation est tout à fait acceptable pour une diffusion à de petits angles (inférieurs à 15°) et qu'elle permet d'obtenir des signatures f d'une résolution et d'une précision excellentes pour des angles de diffusion inférieurs à 5°. Pour estimer (Sin, x Att) à partir du spectre en transmission mesuré h et de 30 la matrice réponse du détecteur REt calibrée, le système selon l'invention utilise avantageusement à nouveau une technique de type MLEM. La figure 12 montre les différents spectres (Sin, x Att) obtenus après inversion MLEM des données expérimentales enregistrées sur le détecteur spectrométrique 7 placé en transmission pour les deux échantillons inspectés (eau et acétone), lesquelles données (spectres en transmission mesurés) sont illustrées à la figure 4. Tous les termes de la matrice de réponse globale A du système ont été 5 calibrés, le procédé selon l'invention pris en exemple consiste ensuite à les combiner selon la formule A = (REd x S x Att)- Ro . Cette combinaison est résumée à la figure 13 où la référence numérique 10 désigne une multiplication ligne à ligne et où la référence numérique 11 désigne une multiplication matricielle. Les matrices de réponse globale A obtenues peuvent être observées à la figure 14 pour l'échantillon 10 d'eau et à la figure 15 pour l'échantillon d'acétone. L'invention peut faire l'objet de nombreuses variantes vis-à-vis du mode de réalisation préférentiel précédemment décrit. Ainsi par exemple, la seconde partie du procédé peut consister en une analyse par diffractométrie classique utilisant la technique EDXRD par exemple. Par ailleurs, le détecteur placé en transmission peut- 15 être un simple détecteur à intégration, auquel cas la première caractéristique de l'objet évaluée aux fins de vérification d'un critère de suspicion n'est pas le numéro atomique effectif du matériau. REVENDICATIONS1. Procédé d'analyse d'un objet caractérisé en ce qu'il se déroule en deux parties, à savoir une première partie comprenant les étapes suivantes : v irradiation de l'objet par un rayonnement photonique incident, / mesure d'un rayonnement transmis par l'objet, dit rayonnement en transmission mesuré, à l'aide d'un détecteur placé en transmission, / détermination d'au moins une première caractéristique de l'objet à partir du rayonnement en transmission mesuré, V vérification de la réalisation d'au moins un critère de suspicion portant sur la première caractéristique de l'objet et traduisant le fait que l'objet contient un matériau potentiellement suspect pour une application considérée ; et une seconde partie exécutée uniquement lorsque le critère de suspicion est réalisé, laquelle seconde partie comprend : v acquisition d'un spectre énergétique, dit spectre en diffusion mesurée, diffusé par l'objet selon un angle de diffusion (0) non nul inférieur à 15° à l'aide d'un détecteur spectrométrique placé en diffusion, / détermination d'une seconde caractéristique de l'objet à partir - soit du spectre en diffusion mesuré, - soit à la fois du spectre en diffusion mesuré et du rayonnement en transmission mesuré, / comparaison - soit de la seconde caractéristique de l'objet, - soit à la fois des première et seconde caractéristiques de l'objet, 25 avec des caractéristiques de matériaux étalons mémorisées dans une base de données, aux fins d'identification d'un matériau constitutif de l'objet. 2. Procédé d'analyse selon la revendication 1, caractérisé en ce que le détecteur placé en transmission est un détecteur spectrométrique. 30 3. Procédé d'analyse selon la revendication 2, caractérisé en ce que la première caractéristique est choisie parmi une forme, une dimension, un coefficient d'atténuation pour une ou plusieurs énergies prédéterminées ou pour une plage d'énergies prédéterminée, une indication quant à la nature d'un matériau constitutifde l'objet, un numéro atomique effectif. 4. Procédé d'analyse selon l'une des revendications précédentes, caractérisé en ce que l'étape de vérification de la réalisation du critère de suspicion comprend la 5 comparaison de la première caractéristique déterminée pour l'objet à une plage de valeurs prédéfinie eu égard à l'application considérée. 5. Procédé d'analyse selon l'une des revendication 1 à 4, caractérisé en ce que la deuxième caractéristique est choisie parmi un transfert de quantité de mouvement 10 et une fonction (f), dite signature en diffusion, représentative de l'objet et montrant des pics de diffusion, reconstruite à partir à la fois du spectre en diffusion mesuré et du spectre en transmission mesuré dans le cas où le détecteur placé en transmission est un détecteur spectrométrique. 15 6. Procédé d'analyse selon la revendication 5, caractérisé en ce que la reconstruction de la signature en diffusion de l'objet met en oeuvre une méthode basée sur une approche de type problème inverse et comprend une opération de construction d'une matrice de réponse globale (A) du système de détection, laquelle matrice de réponse globale (A) établit une relation entre une énergie détectée par le 20 détecteur placé en diffusion et un transfert de quantité de mouvement dans l'objet . 7. Procédé d'analyse selon la revendication 6, caractérisé en ce que l'opération de construction de la matrice de réponse globale (A) du système de détection est faite à partir d'un spectre incident atténué estimé (Sinc x Att) et d'une matrice de 25 réponse angulaire calibrée (R9) du système de détection. 8. Procédé d'analyse selon la revendication 7, caractérisé en ce que l'opération de construction de la matrice de réponse globale (A) du système de détection est faite à partir, de plus, d'une matrice de réponse (REd) du détecteur spectrométrique 30 placé en diffusion (6) et d'une matrice de réponse (REt) du détecteur spectrométrique placé en transmission (7). 9. Système de détection pour l'analyse d'un objet comprenant : ,( une source (1) de rayonnement photonique,/ une zone (4) pour la réception d'un objet à analyser, / des moyens de détection (6, 7), en aval de la zone de réception d'objet (4), pour la mesure d'un rayonnement transmis par l'objet, dit spectre en transmission mesuré et pour l'acquisition d'un spectre énergétique, dit spectre en diffusion mesuré, diffusé par l'objet selon un angle de diffusion 0 non nul inférieur à 15° ; ces moyens de détection comprenant un détecteur spectrométrique (6). / des moyens informatiques (8) de traitement comprenant . des moyens de détermination d'au moins une première caractéristique de l'objet à partir du spectre en transmission mesuré, . des moyens de vérification de la réalisation d'au moins un critère de suspicion portant sur la première caractéristique de l'objet et traduisant le fait que l'objet contient un matériau potentiellement suspect pour une application considérée, . des moyens de détermination d'une deuxième caractéristique de l'objet à partir soit du spectre en diffusion mesuré, soit à la fois du spectre en diffusion mesuré et du rayonnement en transmission mesuré, . et des moyens de comparaison à des caractéristiques de matériaux étalons mémorisées dans une base de données, soit de la seconde caractéristique de l'objet, soit des première et seconde caractéristiques de l'objet, aux fins d'identification d'un matériau constitutif de l'objet. 10. Système de détection selon la revendication 9, caractérisé en ce que les moyens de détection comprennent de plus un autre détecteur (7), placé en transmission, et en ce que le détecteur spectrométrique (6) est placé en diffusion. 25 11. Système de détection selon la revendication 9, caractérisé en ce que les moyens de détection comprennent de plus des moyens de positionnement pour le positionnement successif dudit détecteur spectrométrique (6) en une position de transmission pour la mesure du rayonnement en transmission mesuré, et en une 30 position de diffusion pour l'acquisition du spectre en diffusion mesuré. 12. Système de détection selon l'un des revendications 9 à 11, caractérisé en ce que les moyens de vérification de la réalisation du critère de suspicion sont aptes à comparer la première caractéristique déterminée pour l'objet à une plage de valeursprédéfinie eu égard à l'application concernée. 13. Système de détection selon l'une des revendications 9 à 12, caractérisé en ce que : .7 la première caractéristique est choisie parmi une forme, une dimension, un coefficient d'atténuation pour une ou plusieurs énergies ou pour une plage d'énergies prédéterminée, une indication quant à la nature d'un matériau constitutif de l'objet, un numéro atomique effectif, et V la deuxième caractéristique est choisie parmi un transfert de quantité de mouvement et une fonction (f), dite signature en diffusion, représentative de l'objet et montrant des pics de diffusion, reconstruite à partir à la fois du spectre en diffusion mesuré et du spectre en transmission mesuré dans le cas où le détecteur placé en transmission est un détecteur spectrométrique. 14. Système de détection selon l'une des revendications 9 à 13, caractérisé en ce que le détecteur placé en transmission (7) est un détecteur spectrométrique et en ce que les moyens informatiques (8) de traitement sont aptes à mettre en oeuvre une méthode basée sur une approche de type problème inverse et à construire une matrice de réponse globale (A) du système de détection, laquelle matrice de réponse globale (A) établit une relation entre une énergie détectée par le détecteur placé en diffusion et un transfert de quantité de mouvement dans l'objet, la matrice de réponse globale (A) étant construite à partir d'un spectre incident atténué (Sinc x Att) estimé en utilisant le spectre en transmission mesuré, et d'une matrice de réponse angulaire calibrée (R0) du système de détection. 15. Système de détection selon l'une des revendications 9 à 14, caractérisé en ce que le détecteur spectrométrique (6) placé en diffusion est configuré de façon à présenter un axe de détection (D) formant, avec un axe central (Z) du rayonnement incident, un angle de diffusion 0 Inférieur à 5°.30