FR 2902527 A1 20071221 FR 0605460 A 20060620 DISPOSITIF DE LOCALISATION TRIDIMENSIONNELLE DE SOURCES DE RAYONNEMENT _DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne un dispositif de localisation tridimensionnelle de sources de rayonnement. Par rayonnement , on entend un rayonnement X et/ou gamma. L'invention s'applique notamment aux domaines suivants : - médecine nucléaire, - surveillance des centrales nucléaires, -surveillance des aéroports, et - radiographie et contrôle non destructif. ÉTAT DE LA TECHNIQUE ANTÉRIEURE pourra se reporter aux documents suivants . [1] Mise en oeuvre et étude des propriétés spectrales de la gamma-caméra ISGRI, thèse présentée et soutenue par Olivier Limousin le 27 novembre 2001 [2] The basic component of the ISGRI CdTe gamma-ray camera for space telescope IBIS on board the INTEGRAL satellite, O. Limousin et al., Nuclear 25 Instruments and Methods in Physics Research A 428 (1999) 216-222 [3] Polycell: the elementary detection unit of the ISGRI CdTe & gamma-ray camera, O. Limousin et al., Nuclear Instruments and Methods in Physics Research A 458 1-2 (2001) 551-557 [4] A basic component for ISGRI, the CdTe gamma camera on board the Integral satellite, M. Arques et al., IEEE Transactions on nuclear science, Vol. 46, N 3, 1999, 181-186 [5] The ISGRI CdTe gamma-ray camera: first steps, O. Limousin et al., Nuclear Instruments and Methods in Physics Research A 442 (2000) 244-249 [6] Qualification model of the space ISGRI CdTe gamma-ray camera, O. Limousin et al., Nuclear Instruments and Methods in Physics Research A 471 (2000) 174-178. Pour l'imagerie ou la spectro-imagerie dans le domaine des rayons gamma d'environ 20 keV à 2 MeV, les systèmes de détection connus utilisent généralement l'une des trois techniques suivantes : imagerie à masque codé, imagerie Compton et focalisation gamma. L'imagerie à masque codé consiste à enregistrer, sur un détecteur de position qui est sensible aux photons gamma, l'ombre d'un masque qui est partiellement opaque à ces photons et dont le motif est connu, ce masque étant éclairé par la source de rayonnement gamma que l'on cherche à détecter et à localiser. Le détecteur de position échantillonne l'ombre portée. La position de la source est reconstruite par calcul. Il n'est pas nécessaire d'enregistrer l'énergie déposée par les photons gamma pour constituer l'image, et cette technique ne donne pas d'information sur l'énergie. En revanche, elle peut être appliquée sur de grands domaines spectraux. Toutefois, elle est généralement limitée aux basses énergies, inférieures à 1 MeV. Ses principaux intérêts résident dans sa bonne résolution angulaire, et dans la soustraction automatique du bruit de fond, améliorant la sensibilité (rapport signal (source) sur bruit (environnement)). La résolution angulaire peut être typiquement de l'ordre de quelques minutes d'arc. Elle est fixe et dépend de la taille des éléments du masque, de la capacité du détecteur de position à échantillonner ces éléments, et de la distance du masque au plan de détection. Les systèmes d'imagerie à masque codé ont une hauteur de quelques centimètres à quelques mètres. Leur champ de vue est limité à quelques degrés, voire quelques dizaines de degrés. La technique de l'imagerie Compton utilise au moins deux détecteurs de rayonnements gamma qui sont capables d'enregistrer les positions des interactions d'un photon gamma incident avec ces détecteurs et les énergies déposées par ce photon en ces positions. Le principe de cette technique repose sur le fait qu'un photon gamma, ayant une énergie incidente suffisamment élevée, présente une forte probabilité de ne céder qu'une partie de cette énergie à un électron qui est lié aux atomes du premier détecteur remontré, à la suite d'un choc inélastique, puis de diffuser vers le second détecteur où il cède tout ou partie de son énergie résiduelle. L'angle de diffusion est lié à l'énergie du photon incident et à l'intensité du dépôt d'énergie lors de l'interaction avec le premier détecteur. Un système dit Compton est un spectroimageur. Dans un tel système, la direction du photon incident et donc la position de la source émettrice sont reconstruites par le calcul. La résolution angulaire du système et la précision de localisation dépendent de la précision de mesure de l'énergie et de la statistique de comptage (temps de pose et flux de la source). Dans un système Compton idéal, dont la mesure spectrale est parfaite , la résolution angulaire est systématiquement limitée par l'effet d'élargissement Doppler : l'angle de diffusion est affecté d'une incertitude liée au fait que l'électron sur lequel diffuse le photon n'est pas au repos au moment du choc. La résolution angulaire est limitée à environ 2 à 3 . En outre, les systèmes Compton ont un champ de vue très étendu, potentiellement 4nt stéradians. La limitation de la sensibilité de ce type d'imageur vient du traitement des coïncidences fortuites qui ne sont pas dues à la diffusion Compton. De plus, ces systèmes sont compacts et peuvent avoir des dimensions réduites à quelques centimètres. Le principe d'un imageur spectrométrique utilisant la focalisation gamma est fondé sur l'utilisation d'une lentille gamma qui est équipée de cristaux sur lesquels les photons gamma sont diffusés (diffusion de Bragg). La résolution angulaire d'un tel imageur est excellente - elle vaut environ 1 minute d'arc -- mais l'imageur est limité par une distance focale gigantesque qui dépend fortement de l'énergie des photons incidents et vaut environ 15 m à 122 keV et environ 80 m à 511 keV. En outre, l'imagerie est directe mais limitée à un domaine spectral étroit et le champ de vue de ce type de système est minuscule : il vaut quelques minutes d'arc. Les applications de la technique de focalisation gamma sont donc plutôt liées à l'astronomie. Les trois techniques que l'on vient de rappeler permettent une résolution en énergie et une localisation angulaire de sources qui sont situées à grande distance, quasiment à l'infini. Ainsi les systèmes de détection correspondants trouvent-t-ils leurs principales applications dans l'astronomie à haute énergie. De plus, tous ces systèmes sont soit encombrants, soit limités quant à leur bande spectrale, soit peu performants quant à leur résolution angulaire à haute énergie. De plus, ils ne sont bien adaptés qu'à la localisation bidimensionnelle de sources de rayonnement gamma que l'on considère comme à l'infini , c'est-à-dire à une distance égale à au moins huit fois la hauteur du télescope incorporant un système de l'un des trois types considérés plus haut. EXPOSÉ DE L'INVENTION La présente invention a pour but de remédier aux inconvénients précédents et propose un dispositif qui permet de localiser des sources de 10 rayonnement dans l'espace. De façon précise, la présente invention a pour objet un dispositif de localisation tridimensionnelle de sources de rayonnement, ce dispositif comprenant : 15 a) un ensemble de détection comprenant : - au moins deux détecteurs de rayonnement, à savoir des premier et deuxième détecteurs, qui sont plans, parallèles et sensibles à la position, chacun des premier et deuxième détecteurs étant apte à fournir 20 une impulsion électrique lorsqu'il reçoit un photon du rayonnement, cette impulsion électrique étant représentative de l'énergie de ce photon et de la position d'impact de ce dernier sur le détecteur qui le reçoit, et 25 - un masque codé comportant un matériau qui est opaque au rayonnement et qui est percé de trous selon un motif prédéfini, ce masque codé étant placé en regard du premier détecteur, ce dernier étant compris entre le masque codé et le deuxième détecteur, ce 30 positionnement relatif permettant d'obtenir respectivement, à partir d'une source que l'on veut5 localiser, des première et deuxième ombres portées sur les premier et deuxième détecteurs qui engendrent alors respectivement des première et deuxième images sous forme électrique de ces première et deuxième ombres portées, et b) des moyens électroniques de traitement des impulsions électriques fournies par les premier et deuxième détecteurs, pour localiser la source, ce dispositif étant caractérisé en ce que 10 les moyens électroniques de traitement sont prévus pour : - sélectionner parmi les impulsions électriques qui sont respectivement fournies par les premier et deuxième détecteurs celles qui sont en 15 coïncidence temporelle, - effectuer la déconvolution de la première image ou de la deuxième image ou des première et deuxième images sur les impulsions électriques sélectionnées, 20 - déterminer la position angulaire, c'est-à-dire la direction, de la source dans l'espace à partir du résultat de cette déconvolution, et - déterminer la distance de la source par rapport à l'ensemble de détection, par analyse de 25 l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs. L'invention peut également être définie comme un dispositif comportant : 30 - un masque codé qui est opaque au rayonnement d'une source à observer, - un spectro-imageur Compton comportant au moins deux détecteurs de rayonnement, à savoir des premier et deuxième détecteurs, qui sont plans, parallèles et sensibles à la position, chacun des premier et deuxième détecteurs étant apte à fournir une impulsion électrique lorsqu'il reçoit un photon du rayonnement, cette impulsion électrique étant représentative de l'énergie de ce photon et de la position d'impact de ce dernier sur le détecteur qui le reçoit, le masque codé étant placé entre la source et ce spectro-imageur, ce dernier étant orienté de manière à ce que la source à observer, que l'on recherche, projette une ombre du masque codé sur le premier détecteur et au moins une partie du second détecteur, et des moyens électroniques de calcul permettant de déterminer d'une part la direction de la source par une déconvolution déterminant une même partie du motif du masque codé, cette même partie apparaissant sur les images fournies par les premier et deuxième détecteurs, et d'autre part la distance de la source, par analyse de l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs. De préférence, les premier et deuxième détecteurs sont des détecteurs semiconducteurs pixellisés. Selon un mode de réalisation préféré du 30 dispositif objet de l'invention, les moyens électroniques de traitement sont en outre prévus pour déterminer, après la détermination de la distance de la source, le spectre énergétique et le flux de la source par analyse des impulsions électriques sélectionnées et des première et deuxième images, sur la dynamique spectrale totale du dispositif. Selon un premier mode de réalisation particulier du dispositif objet de l'invention, les moyens électroniques de traitement sont prévus pour : - rechercher la direction approximative de la source à partir des impulsions électriques sélectionnées, ce qui permet de déterminer une région d'intérêt, - mesurer la distance de la source par rapport à l'ensemble de détection, à partir des 15 impulsions électriques sélectionnées, - après avoir orienté l'ensemble de détection vers la région d'intérêt, cet ensemble étant ainsi orienté suivant une direction, reconstruire la direction de la source dans la région d'intérêt, par 20 décodage du motif prédéfini du masque codé, projeté sur l'un au moins des premier et deuxième détecteurs, - déterminer la distance de la source par rapport à l'ensemble de détection, par analyse de l'image du masque codé et du grandissement de cette 25 image sur l'un au moins des premier et deuxième détecteurs, et - déterminer le spectre énergétique et le flux de la source, par analyse des impulsions électriques sélectionnées et des première et deuxième 30 images, sur la dynamique spectrale totale du dispositif. Selon un deuxième mode de réalisation particulier du dispositif objet de l'invention, lorsque l'ensemble de détection est orienté vers une région d'intérêt, cet ensemble étant ainsi orienté suivant une direction, les moyens électroniques de traitement sont prévus pour : - reconstruire la direction de la source dans la région d'intérêt, par décodage du motif prédéfini du masque codé, projeté sur l'un au moins des premier et deuxième détecteurs, - déterminer la distance de la source par rapport à l'ensemble de détection, par analyse de l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs, - déterminer la distance de la source par rapport à l'ensemble de détection, par analyse des impulsions électriques sélectionnées, et -déterminer le spectre énergétique et le flux de la source, par analyse des impulsions électriques sélectionnées et des première et deuxième images, sur la dynamique spectrale totale du dispositif. De préférence, le masque codé comporte un ensemble de pavés, chacun des premier et deuxième détecteurs comporte un ensemble de pixels et la taille des plus petits pavés du masque est au moins égale à trois fois la taille des pixels des premiers et deuxième détecteurs. Selon un mode de réalisation préféré de l'invention, les moyens électroniques de traitement sont en outre prévus pour calculer la proportion d'événements fortuits qui sont présents dans des données constituées par les première et deuxième images et les impulsions électriques sélectionnées, pour générer des pseudo-événements fortuits afin de les soustraire aux données. De préférence, le dispositif objet de l'invention comprend en outre des moyens de limitation de bruit en dehors du champ visuel de ce dispositif. Les moyens de limitation de bruit peuvent comprendre des parois en plomb ou en tantale. Les premier et deuxième détecteurs peuvent être identiques. BRÈVE DESCRIPTION DES DESSINS La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés ci-après, à titre purement indicatif et nullement limitatif, en faisant référence aux dessins annexés sur lesquels : - la figure 1 est une vue schématique d'un mode de réalisation particulier du dispositif objet de l'invention, - la figure 2 est une vue en perspective 25 schématique du masque codé et d'un détecteur faisant partie du dispositif de la figure 1, - la figure 3 illustre schématiquement le principe de l'imagerie Compton qui est utilisé dans le dispositif de la figure 1, et - la figure 4 est une vue schématique du masque codé et de l'un des détecteurs qui font partie du dispositif de la figure 1. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Un dispositif conforme à l'invention combine la technique d'imagerie à masque codé et la technique d'imagerie Compton . Ce dispositif utilise simultanément un système de détection Compton , comportant au moins deux détecteurs de rayonnement, et un masque codé placé devant la surface sensible de l'un des détecteurs. La réponse du dispositif objet de l'invention offre les avantages de la bonne résolution angulaire du masque codé sur toute la dynamique spectrale du système de détection Compton , conserve une résolution angulaire moyenne sur un champ de vue très étendu et, de surcroît, permet de mesurer la distance entre la source à localiser et ce dispositif. La présente invention s'affranchit des limitations des dispositifs connus et permet de : - localiser avec précision en trois dimensions une (ou des) source(s) ponctuelle(s) de rayonnement sur une gamme spectrale étendue (20 keV à 10 MeV pour des rayonnements gamma), - localiser une source ou simultanément des sources du champ de vue du masque codé en trois dimensions, c'est-à-dire déterminer la position de sources de rayonnement qui sont situées à des distances finies et sont éventuellement différentes (le dispositif fournit une image bidimensionnelle des sources et la distance à laquelle il se trouve des sources), - localiser avec une résolution limitée à 2 ou 3 une (ou des) source(s) ponctuelle(s) ou étendue(s), sur une gamme spectrale étendue allant d'environ 50 keV à 10 MeV pour des rayonnements gamma (par source étendue , on entend une source dont les dimensions angulaires sont supérieures à la résolution du masque codé), - déterminer la distance entre le dispositif et une (ou plusieurs) source(s) ponctuelle(s) ou étendue(s), qui sont éventuellement situées à des distances différentes du dispositif, - mesurer l'énergie de chaque photon incident individuellement, - mesurer le spectre en énergie de la (ou des) source(s), ceci dans un vaste domaine allant de quelques keV à quelques MeV pour des rayonnements gamma, - déterminer la date d'arrivée de chaque photon, de déterminer le flux de la (ou des) source(s) ponctuelle(s) détectée(s), de déterminer précisément le flux de la (ou des) source(s) ponctuelle(s) ou étendue(s) par le traitement approprié des coïncidences fortuites pour la reconstruction des photons détectés lors des diffusions Compton, et - déterminer la variation de flux de la (ou 30 des) source(s) détectée(s). L'invention tire avantage de la compacité des systèmes Compton et des systèmes à masque codé. Un dispositif conforme à l'invention comporte au moins deux détecteurs plans qui ont de préférence des résolutions spatiale et spectrale élevées, travaillent en coïncidence et sont associés à un masque codé qui permet d'affiner l'image d'une prise de vue Compton et de déterminer la position de la source et sa distance au dispositif. On obtient ainsi une cartographie en trois dimensions des sources observées. On aboutit, en définitive, à la réalisation d'une caméra X et/ou gamma tridimensionnelle, plus simplement appelée caméra 3D par la suite. Par ailleurs, pour fabriquer les détecteurs utilisés dans l'invention, on emploie de préférence une technologie qui permet leur utilisation à température ambiante, ce qui d'une part simplifie leur mise en oeuvre et d'autre part permet d'optimiser le poids et l'encombrement de la caméra 3D pour qu'elle soit transportable. L'association de ces détecteurs à des circuits électroniques fortement intégrés, de type ASIC (pour Application Specific Integrated Circuit) et de consommation modérée, permet de rendre la caméra 3D portative et autonome. De plus, cette caméra peut être rendue capable de travailler dans un milieu bruité, par exemple une centrale nucléaire, sans perte notable de sensibilité. Ceci peut être fait de manière classique, grâce à l'adjonction de blindages (par exemple en plomb, en tungstène ou en d'autres métaux lourds), mais aussi par un traitement qui est spécifique aux mesures Compton et qui sera exposé ultérieurement. Les blindages ont pour effet de réduire l'angle d'observation de la caméra, alors que le traitement spécifique aux mesures Compton n'a pas cet inconvénient. L'invention peut utiliser l'une ou l'autre de ces protections contre le bruit ou l'association des deux. La figure 1 est une vue schématique d'un dispositif conforme à l'invention. Ce dispositif est destiné à la localisation tridimensionnelle de sources de rayonnement. Ce dispositif comprend un ensemble de détection 2 et des moyens électroniques de traitement 4 qui sont munis de moyens 5 d'affichage des résultats qu'ils fournissent. L'ensemble de détection 2 comprend un premier détecteur 6, apte à détecter un rayonnement, et un deuxième détecteur 8, également apte à détecter ce rayonnement. Les détecteurs 6 et 8 sont plans, parallèles et sensibles à la position. Chacun de ces détecteurs est apte à fournir une impulsion électrique lorsqu'il reçoit un photon du rayonnement, cette impulsion électrique étant représentative de l'énergie de ce photon et de la position d'impact de ce dernier sur le détecteur qui le reçoit. - L'ensemble de détection 2 comprend aussi un masque codé 10 qui est fait d'un matériau opaque au rayonnement. Ce masque codé 10 est percé de trous 12 selon un motif prédéfini ; il est placé en regard du premier détecteur 6 et ce dernier est compris entre le masque codé 10 et le deuxième détecteur 8. Le dispositif comprend aussi un blindage 14 sur lequel on reviendra par la suite. Au moyen du dispositif, on veut localiser une source 16 qui émet un rayonnement (X et/ou gamma) 18. Pour ce faire, on oriente le dispositif vers cette source. On reviendra sur cette orientation par la suite. L'axe Z du dispositif est alors dirigé vers 15 la zone où se trouve la source. La source 16 engendre une ombre portée sur chacun des premier et deuxième détecteurs 6 et 8 par l'intermédiaire du masque codé 10. Les premier et deuxième détecteurs 6 et 8 engendrent alors 20 respectivement des première et deuxième images sous forme électrique de cette ombre portée. Les moyens électroniques de traitement 4 sont prévus pour traiter les impulsions électriques qui sont fournies par les premier et deuxième détecteurs, 25 pour localiser la source 16. Conformément à l'invention, les moyens électroniques de traitement 4 sélectionnent les impulsions électriques qui sont respectivement fournies par les premier et deuxième détecteurs 6 et 8 et qui 30 coïncident temporellement. Puis les moyens électroniques de traitement 4 effectuent la déconvolution de la première image ou de la deuxième image ou des deux. Cette déconvolution est effectuée sur les impulsions électriques sélectionnées, de manière à rechercher tout ou partie du motif du masque 10 dans les images des impulsions sélectionnées des détecteurs 6 et 8. Compte tenu du motif choisi pour le masque, à une position de la source 16 ne correspond qu'un seul motif possible d'ombre portée 22. La meilleure corrélation entre tout ou partie du motif du masque 10 et l'ombre portée 22 enregistrée dans les détecteurs 6 et 8 donne la direction de la source 16. Puis les moyens électroniques de traitement 4 déterminent la position de la source 16 dans l'espace à partir de la première image ainsi déconvoluée ou de la deuxième image ainsi déconvoluée ou des première et deuxième images ainsi déconvoluées. Cette détermination consiste à trouver la direction de la source 16 puis la distance D de la source au premier détecteur 6 par analyse de l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs. On reviendra sur cette question par la suite. Les moyens électroniques de traitement 4 sont en outre prévus pour déterminer le spectre énergétique de la source 16 à partir des impulsions électriques sélectionnées. Pour ce faire, on procède de la façon 30 suivante : on classe les impulsions sélectionnées par intervalle d'énergie, et l'on représente le résultat sous forme d'histogramme. Puis on compare le résultat de ce classement à celui des sources radiatives connues afin d'identifier leur présence. La présence et la mesure d'un ou plusieurs pics dans cette distribution permet l'identification d'une source radioactive donnée 16. Le mode de fonctionnement du dispositif conforme à l'invention dépend de la connaissance préalable de la scène observée. On distingue deux cas. 1. Premier cas : aucune position de source n'est connue a priori dans le champ de détection. La procédure de recherche est la suivante : on recherche approximativement la position de la (ou des) source(s) avec le mode Compton , c'est-à-dire que l'on détermine statistiquement la (ou les) direction(s) d'origine des photons par analyse de leur angle de diffusion et de leur dépôt d'énergie dans les deux détecteurs ; on reviendra sur cette question par la suite ; cette détermination est faite sur un champ de vue très étendu qui va jusqu'à 4n stéradians pour un détecteur non blindé et 2n stéradians pour un détecteur blindé ; puis la (ou les) source(s) sont repérées à 2 ou 3 près ; et l'on détermine une région spatiale d'intérêt ; on mesure la distance de la (ou des) sources, ponctuelle(s) ou étendue(s) à partir des données Compton ; ces données sont obtenues par la recherche en mode Compton qui consiste à calculer la projection, sur l'axe Z du dispositif, de l'intersection des directions d'origine de chaque photon, déterminées à partir des impulsions sélectionnées ; on repointe le dispositif dans la région spatiale d'intérêt, on reconstruit la position bidimensionnelle de la (du des) source(s) dans le champ repointé, par décodage du motif du masque, qui est projeté sur le (ou les) plan(s) imageur(s), c'est-à-dire le (ou les) plan (s) du (ou des) détecteur (s) ; le décodage du motif du masque consiste à rechercher la meilleure corrélation entre l'image enregistrée sur le premier ou le deuxième détecteur à partir des impulsions sélectionnées, et le motif du masque codé, une direction de la source correspondant à une portion donnée du masque et donc à un motif unique sur les images ; par position bidimensionnelle de la source dans le champ , on entend la position de la source dans un plan perpendiculaire à l'axe Z du dispositif, sans information sur la distance de la source au dispositif ; on détermine précisément la (ou les) distance(s) de la (ou des) source(s) ponctuelle(s), par analyse de l'image du masque 10 et du grandissement de l'ombre portée de ce masque sur le ou les plans détecteurs, par exemple en déterminant la (ou les) distance(s) de la (ou des) source(s) ponctuelle(s) suivant leur énergie, c'est-à-dire en déterminant la distance de l'une des sources en sélectionnant les photons dans une bande d'énergie particulière, caractéristique de la source ; cette bande caractéristique peut être déterminée par l'analyse préalable du spectre ; en ce qui concerne la distance, on précise qu'elle peut être calculée par rapport au premier détecteur ou au deuxième détecteur ou au masque, ce qui revient au même si l'on connaît la distance entre les détecteurs 6 et 8 et entre le masque 10 et le détecteur 6 ; on détermine le spectre et le flux (c'est-à-dire la distribution de l'énergie des photons et le nombre de photons émis par unités de temps et d'énergie) de la (ou des) sources étendue(s) ou ponctuelle(s), par analyse des données Compton et des données résultant du masque codé, sur la dynamique spectrale totale du dispositif. 2. Deuxième cas : l'existence d'une (ou de) 1.5. source(s) est connue dans une région spatiale d'intérêt et la direction de visée du dispositif est fixe. La procédure de recherche est la suivante : on pointe le dispositif directement sur la région d'intérêt ; 20 - on reconstruit la position bidimensionnelle de la (ou des) source(s) dans le champ repointé, par décodage du motif de masque projeté sur le (ou les) plan (s) imageur(s) ; on détermine la (ou les) distance(s) de 25 la (ou des) source(s), par analyse de l'image du masque et de son grandissement sur le ou les plans détecteurs (suivant l'énergie et la distance) ; - on détermine la (ou les) distance(s) de la (ou des) source(s) ponctuelle(s) ou étendue(s), par 30 analyse de l'image Compton , en particulier si des sources ne sont pas dans le champ de vue du masque codé ; û on détermine le spectre et le flux de la (ou des) source(s), par analyse des données Compton et des données résultant du masque codé, sur la dynamique spectrale totale du dispositif. On indique ci-après le principe de base 10 d'utilisation du masque codé 10 en faisant référence à la figure 1 déjà décrite et à la figure 2. Cette figure 2 est une vue schématique en perspective du masque codé, pourvu des perçages 12, et du premier détecteur plan 6. 15 La position angulaire d'une source ponctuelle lointaine (non représentée) par rapport au détecteur 4 est définie dans un espace à deux dimensions. On appelle source ponctuelle , une 20 source dont la dimension est inférieure à la résolution de l'appareillage de détection constitué par le détecteur 6 et le masque 10. La résolution est un paramètre qui dépend de la taille des trous 12 du masque, de la taille des pixels du détecteur 6 et de la 25 distance d du détecteur 6 au masque 10. Le masque codé est constitué d'une plaque d'un matériau qui est opaque au rayonnement (par exemple le plomb, le tantale ou le tungstène) et qui est percé des trous 12 selon un motif prédéfini. 30 Lorsqu'il est exposé à une source de rayonnement, une partie de ce rayonnement est arrêtée par le matériau opaque du masque alors qu'une autre partie traverse les trous et crée sur le plan détecteur un motif composé de zones claires 20 et de zones obscures 22. Ce motif est directement lié à la direction de la source. A partir du codage connu du masque et du motif reçu sur le plan détecteur, on peut déterminer précisément la direction d'une source ponctuelle dans le champ de vue par un traitement informatique approprié. A ce sujet, on se reportera aux documents suivants . [7] the INTEGRAL IBIS/ISGRI. System point 15 spread function and source location accuracy, A. Gros et al., Astronomy and Astrophysics 411 (2003) L179 [8] INTEGRAL/IBIS scientific data analysis, A. Goldwim et al., Astronomy and Astrophysics 411 (2003) L223. 20 On précise que le quadrillage présent sur le masque 10 de la figure 2 correspond au motif du masque codé 10 (les éléments opaques au rayonnement 25 sont représentés en noir e_t, les éléments 12 transparents au rayonnement en blanc) et que le quadrillage présent sur le détecteur 6 de cette figure 2 correspond à l'ombre portée 22 par la source 16 de rayonnement de ce masque codé sur le détecteur 6. 30 On explique ci-après le principe de localisation d'une source de rayonnement en utilisant 10 deux plans détecteurs. Une telle localisation est une reconstruction Compton . On fait référence à la figure 3. Les plans détecteurs ont encore les références 6 et 8. L'axe du dispositif a encore la référence Z. Cet axe Z est orthogonal aux plans détecteurs 6 et 8. On définit un repère O1X1Y1Z où le point O1 appartient au plan détecteur 6 et les axes X1 et Y1 sont situés dans ce plan détecteur 6 et sont perpendiculaires l'un à l'autre et à l'axe Z. On définit aussi un repère O2X2Y2Z où le point 02 appartient aux plan détecteur 8 et les axes X2 et Y2 sont situés dans ce plan détecteur 8 et sont perpendiculaires l'un à l'autre et à l'axe Z. La localisation d'une source est déterminée par l'étude des évènements en coïncidence entre les deux détecteurs plans 6 et 8. Cette étude est basée sur l'effet Compton selon lequel un photon du rayonnement dépose de l'énergie dans le premier détecteur 6, est dévié, puis dépose le reste de son énergie dans le deuxième détecteur 8. Dans un premier temps, on déduit l'angle 02 entre la direction 24 du photon incident et la direction 26 du photon diffusé entre les deux détecteurs plans, grâce à la mesure précise du dépôt d'énergie dans chacun des deux détecteurs plans. Au sujet du calcul de l'angle 02, on se reportera au document suivant : [9] Gamma-ray burst detection and localization capabilities of the IBIS/INTEGRAL telescope Compton mode, Marcinkowski R. et al., Il Nuovo Cimento 28 (2005) 845. Puis on mesure l'angle 01 entre les plans des détecteurs du dispositif et la direction 26 du photon diffusé, grâce à la connaissance précise des points respectifs d'impact dans les deux détecteurs plans. Dans chaque détecteur, le point d'impact correspond à un pixel de ce détecteur. L'angle 01 est donné par la formule suivante : 81 = arctan ([((x2-x1) 2 + (y2-yl) 2) 1/2]/d) où x1 et yl sont les coordonnées du point d'impact dans le détecteur 6 et x2 et 172 sont les coordonnées du point d'impact dans le détecteur 8. La direction 24 de la source 16 de photons par rapport à l'axe Z du dispositif, qui permet de reconstruire l'image bidimensionnelle, et la distance entre le dispositif et la source s'obtiennent par l'étude statistique des angles 01 et 02 pour chaque photon sélectionné. La direction d'un photon incident étant sur le cône d'axe 26 et de demiangle au sommet 02, on déttermine la direction 24 de la source en calculant l'intersection de tous les cônes associés à chaque photon sélectionné. Pour obtenir la distance entre le dispositif et la source, on procède de la façon suivante : on calcule la projection, sur l'axe Z du dispositif, de l'intersection des cônes associés à chaque photon, et correspondant à des impulsions sélectionnées. On explique ci-après le principe de la détermination de la position en trois dimensions de la (ou des) source (s) dans le champ de vue du masque codé 10, en faisant référence à la figure 4. Le motif induit par le masque codé 10 sur le (ou les) détecteur(s) (seul le détecteur 6 est considéré sur la figure 4), dépend de la distance h de la source 16 au masque 10, du fait que le flux de photons provenant d'une source proche du détecteur n'est pas parallèle. Le motif est plus ou moins grossi selon que le dispositif est plus ou moins proche de la source. Le traitement informatique du motif permet d'accéder à la distance du dispositif à la source. Cette méthode fonctionne même s'il y a plusieurs sources dans le champ de vue, y compris si ces sources se trouvent à des distances différentes. Grâce à l'information sur la position angulaire, qui est déterminée comme nous l'avons vu plus haut, et grâce à l'information sur la distance entre la source et le dispositif, on aboutit à une image tridimensionnelle des sources proches qui sont situées dans le champ du dispositif. L'exploitation des grossissements différents des motifs de masque, motifs qui sont enregistrés sur les deux (ou plus) détecteurs plans et parallèles les uns aux autres, permet d'améliorer la précision de mesure de la distance de la source, étant donné qu'il s'agit de mesures indépendantes et multiples du même paramètre. Selon un mode de réalisation préféré, après avoir recherché et identifié une partie connue du motif du masque sur l'un des deux détecteurs, on évalue tout d'abord la direction de la source. Puis, connaissant les dimensions des motifs du masque, on mesure la dimension de leur ombre sur la surface détectrice considérée, qui est à une distance connue de ce masque, pour en déduire par triangulation la distance de la source. On effectue alors une mesure analogue avec l'autre surface détectrice, afin d'obtenir une autre évaluation de la distance de la source. Les deux évaluations sont alors moyennées afin d'améliorer la précision. Selon un autre mode de réalisation possible, qui permet une détection sur un angle d'ouverture plus grand, la direction de la source est déterminée comme ci-dessus, par recherche et identification d'une partie connue du motif du masque sur l'un des deux détecteurs, on évalue tout d'abord la direction de la source. Puis on utilise uniquement la différence de grossissement des images formées sur les deux surfaces détectrices pour évaluer la distance de la source par triangulation, la distance entre ces deux surfaces détectrices étant connue. On précise ci-après la technologie des deux détecteurs plans. A partir du principe de localisation que l'on a vu précédemment, on comprend que la résolution spatiale et la résolution spectrale (précision de mesure de la position d'interaction et précision de mesure de l'énergie des photons) sont des paramètres qui déterminent les performances du dispositif au premier ordre. Du fait de l'utilisation de données Compton, la résolution des mesures fournies par ce dispositif dépend beaucoup de la résolution en énergie des détecteurs. On utilise donc de préférence des détecteurs semi-conducteurs pixellisés, dont les résolutions spectrale et spatiale sont élevées. Ces détecteurs peuvent être en CdTe, en silicium ou en germanium. On peut par exemple utiliser un détecteur en Si et un détecteur en CdTe ou deux détecteurs en CdTe. A titre d'exemple, on a développé des détecteurs en CdTe, qui ont une haute résolution en énergie et peuvent être utilisés dans le dispositif de l'invention ou caméra 3D. Leurs performances spectrales sont suffisantes pour trouver des sources à 1 mètre avec une précision de 1 cm, et ceci dans une bande spectrale qui s'étend d'environ 200 keV à environ 2 MeV pour les rayonnements gamma. Au sujet de ces détecteurs, on peut se reporter aux documents suivants . Documents [1] à [6] [10] 3D modeling of Cd(Zn)Te detectors for the SIMBOL-X space mission, B.P.F. Dirks et al., High Energy Detectors in Astronomy, Proceedings of SPIE vol. 5501, Glasgow (juin 2004) [11] Leakage current measurements on pixeltade CdZnTe detectors, B.P.F. Dirks et al., Nuclear Instruments and Methods in Physics Research A (2006). A titre d'exemple, on peut utiliser des détecteurs dont les pixels forment des carrés de 0,5 mm à 10 mm de côté. On considère ci-après la minimisation du 10 bruit du dispositif. L'utilisation d'évènements en coïncidence, que l'on a explicitée plus haut, permet non seulement la localisation mais aussi une première sélection des évènements correspondant aux sources parmi de nombreux 15 évènements de bruit. Elle permet donc au dispositif de détection de travailler dans un environnement potentiellement hostile et fortement bruité. On considère ci-après la suppression des 20 événements fortuits. Certains évènements détectés par la caméra 3D peuvent être fortuitement en coïncidence entre les deux détecteurs 6 et 8. Du fait de leur caractère fortuit, ils n'ont rien à voir avec la source à étudier 25 et diminuent la sensibilité de la caméra lors de la production de l'image Compton . Cet effet peut être minimisé par un traitement informatique spécifique qui met en valeur les spécificités de la mesure Compton : le traitement utilise la sélection d'impulsions 30 arrivant sur les première et seconde surfaces détectrices de façon simultanée vis-à-vis de la résolution temporelle de ces surfaces détectrices. Le traitement spécifique consiste alors à considérer comme du bruit, et donc à éliminer, toutes les impulsions n'apparaissant pas simultanément sur les première et seconde surfaces détectrices. On a développé un logiciel de traitement permettant de corriger l'effet mentionné ci-dessus à partir des caractéristiques des deux détecteurs 6 et 8. Ce logiciel calcule la proportion d'évènements fortuits, qui sont présents dans les données de la caméra 3D, et génère des pseudo-événements fortuits afin de les soustraire aux données reconstruites à partir des impulsions électriques sélectionnées en coïncidence. Cette méthode est indispensable pour l'estimation précise du flux de sources faibles (une source faible étant une source dont on détecte la présence bien qu'elle ne contribue au taux de comptage du dispositif que de manière comparable, voire inférieure, au bruit de fond). On précise ci-après quelques éléments quantitatifs au sujet de la caméra 3D. - La distance séparant les deux détecteurs plans est un paramètre qui intervient dans la sensibilité de la détection. On peut montrer que la sensibilité est optimale lorsque la distance séparant les deux détecteurs plans est du même ordre de grandeur que la diagonal de ces détecteurs plans. - La précision de séparation s du masque codé pour deux sources situées dans un plan perpendiculaire à l'axe Z de la caméra 3D est telle que : s=arctan(m/d) qui est peu différent de m/d (en radians). Et, pour un rapport signal sur bruit N des photons de la source, la précision de localisation 1 est telle que : 1=s/N1/2 (en radians). Dans ces formules, m est la taille des plus petits éléments opaques du masque 10 et d est la distance entre le masque 10 et le premier détecteur. - Pour avoir un bon échantillonnage des pavés du masque avec le détecteur 6 ou les détecteurs 6 et 8, il est préférable que la taille des plus petits pavés du masque soit au moins égale à trois fois la taille des pixels de ces détecteurs 6 et 8. On donne ci-après un exemple d'application 20 numérique à titre purement indicatif et nullement limitatif : - la taille des pixels des détecteurs plans vaut 0,5 mm û m vaut au moins 1,5 mm ; û d vaut 30 cm et l'on obtient s = 0,005 radian ; û on peut donc distinguer deux sources qui sont séparées de 5 mm l'une de l'autre et se trouvent à 1 m du dispositif ; û la précision de localisation 1 vaut 1 mm pour un rapport signal sur bruit égal à 25. 25 30 en utilisant la reconstruction Compton , on peut déterminer la distance d'une source, qui est située à 1 m du dispositif, avec une précision de 5 cm. On donne ci-après un exemple de caméra 3D portable, conforme à l'invention. On utilise deux détecteurs, l'un en Si et l'autre en CdTe, avec des pixels de 0,5 mm de coté. Il faut alors utiliser un masque dont les trous ont une taille minimale de 1,5 mm. Pour avoir une caméra portable, on peut envisager deux détecteurs de 10 cm x 10 cm de coté, et une distance de 30 cm entre le masque 10 et le premier détecteur 6. Un tel dispositif permet une précision de localisation de 1 cm à une distance de 1 m en mode imagerie Compton (grand champ) et de 1 mm à une distance de 1 m en mode imagerie à masque codé. En ajoutant des moyens permettant de limiter le bruit en dehors du champ de vue de la caméra 3D, on obtient une caméra d'environ 10 à 15 kilogrammes qui est facilement transportable. Cette caméra est du genre de celle qui est schématiquement représentée sur la figure 1. Dans l'exemple représenté sur cette figure, les moyens permettant de limiter le bruit en dehors du champ de vue de la caméra sont des parois 14 en plomb ou en tantale, qui entourent l'ensemble formé par les détecteurs 6 et 8 et le masque 10 et recouvrent aussi celle des faces du détecteur 8 qui n'est pas en regard du détecteur 6, comme on le voit. - La réalisation des détecteurs ayant de multiples petits pixels impose l'utilisation de circuits électroniques de lecture très fortement intégrés (ASICS), ce qui tend à limiter la consommation électrique de manière radicale. Une seule voie indépendante de lecture d'un pixel peut consommer seulement de 1 mW à 3 mW. Pour le dispositif donné en exemple, cela correspond à un total de 40 W à 120 W. - Une caméra dont les dimensions sont changées de manière homothétique par rapport à la caméra décrite ci-dessus a les mêmes propriétés optiques. Une plus grande surface des détecteurs permet d'obtenir une meilleure sensibilité. Ce qui précède n'est donc qu'un exemple possible de caméra pour illustrer l'invention. D'autres dimensions sont envisageables, pour répondre aux différentes contraintes de précision de localisation, de sensibilité, de poids et/ou de consommation électrique. On précise en outre que les détecteurs pixellisés peuvent comporter des pixels individuels ou des électrodes segmentées. De plus, on précise que le dispositif permet de localiser précisément et simultanément, en trois dimensions, une (ou des) source(s) présentes dans le champ de vue du masque codé. En outre, le dispositif permet de mesurer l'énergie de chaque photon incident. Pour ce faire, on procède de la façon suivante : chaque pixel du dispositif est associé à une chaîne électronique analogique, capable de produire un signal électrique dont la hauteur de l'impulsion est proportionnelle à l'énergie déposée par le photon dans le pixel. Et ce dispositif permet de déterminer la date d'arrivée de chaque photon sur chaque détecteur. A chaque fois qu'un photon interagit dans l'un des détecteurs, son interaction est signalée à l'électronique de lecture par un signal numérique rapide. Ce signal est synchronisé sur une horloge qui permet d'en faire une datation précise (typiquement de l'ordre de 1 à 10 us. En outre, le dispositif permet de déterminer le flux de la (ou des) source(s) détectée(s). Et le dispositif permet de déterminer précisément le flux de la (ou des) source(s) ponctuelle (s) ou étendue(s), au delà de 50 keV, par un traitement approprié des coïncidences fortuites pour la reconstruction des photons détectés lors des diffusions Compton. Ce traitement est par exemple le suivant : On génère par calcul une liste de pseudoévènements fortuits en associant aléatoirement deux à deux tous les évènements enregistrés par les deux détecteurs. Ces pseudo-évènements subissent le même traitement numérique que les évènements sélectionnés. On calcule ensuite la proportion d'évèriements fortuits dans les évènements sélectionnés. Enfin, on soustrait cette proportion de pseudo-évènements aux évènements solutionnés. Le dispositif permet aussi de déterminer la variation du flux de la (ou des) source(s) détectée(s). Dans les exemples de l'invention que l'on a donnés, on a utilisé deux détecteurs de rayonnement. Cependant, plus de deux détecteurs peuvent être utilisés dans l'invention afin d'étendre la gamme d'énergies de cette dernière jusqu'à 10 MeV voire plus. De plus, dans les exemples donnés, on a surtout considéré la localisation de sources de rayonnement gamma. Mais l'invention permet aussi la localisation de sources de rayonnement X et de sources de rayonnement X et gamma, en choisissant des détecteurs adaptés à de tels rayonnements. REVENDICATIONS 1. Dispositif de localisation tridimensionnelle de sources de rayonnement, ce dispositif comprenant : a) un ensemble de détection (2) comprenant : -au moins deux détecteurs de rayonnement, à savoir des premier (6) et deuxième (8) détecteurs, qui sont plans, parallèles et sensibles à la position, chacun des premier et deuxième détecteurs étant apte à fournir une impulsion électrique lorsqu'il reçoit un photon du rayonnement (18), cette impulsion électrique étant représentative de l'énergie de ce photon et de la position d'impact de ce dernier sur le détecteur qui le reçoit, et un masque codé (10) comportant un matériau qui est opaque au rayonnement et qui est percé de trous (12) selon un motif prédéfini, ce masque codé étant placé en regard du premier détecteur, ce dernier étant compris entre le masque codé et le deuxième détecteur, ce positionnement relatif permettant d'obtenir respectivement, à partir d'une source (16) que l'on veut localiser, des première et deuxième ombres portées sur les premier et deuxième détecteurs qui engendrent alors respectivement des première et deuxième images sous forme électrique de ces première et deuxième ombres portées, et b) des moyens électroniques (4) de traitement des impulsions électriques fournies par les premier et deuxième détecteurs, pour localiser la source,ce dispositif étant caractérisé en ce que les moyens électroniques de traitement (4) sont prévus pour -sélectionner parmi les impulsions 5 électriques qui sont respectivement fournies par les premier et deuxième détecteurs celles qui sont en coïncidence temporelle, - effectuer la déconvolution de la première image ou de la deuxième image ou des première et 10 deuxième images sur les impulsions électriques sélectionnées, - déterminer la position angulaire de la source dans l'espace à partir du résultat de cette déconvolution, et - déterminer la distance de la source par rapport à l'ensemble de détection, par analyse de l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs. 2. Dispositif selon la revendication 1, dans lequel les premier et deuxième détecteurs (6, 8) sont des détecteurs semi-conducteurs pixellisés. 3. Dispositif selon l'une quelconque des revendications 1 et 2, dans lequel les moyens 25 électroniques de traitement (4) sont en outre prévus pour déterminer le spectre énergétique et le flux de la source (16) par analyse des impulsions électriques sélectionnées et des première et deuxième images, sur la dynamique spectrale totale du dispositif. 15 20 4. Dispositif selon l'une quelconque des revendications 1 à 3, dans lequel les moyens électroniques de traitement (4) sont prévus pour - rechercher la direction approximative de la source (16) à partir des impulsions électriques sélectionnées, ce qui permet de déterminer une région d'intérêt, - mesurer la distance de la source (16) par rapport à l'ensemble de détection, à partir des impulsions électriques sélectionnées, - après avoir orienté l'ensemble de détection (2) vers la région d'intérêt, cet ensemble étant ainsi orienté suivant une direction, reconstruire la direction de la source dans la région d'intérêt, par décodage du motif prédéfini du masque codé (10), projeté sur l'un au moins des premier et deuxième détecteurs (6, 8), - déterminer la distance de la source (16) par rapport à l'ensemble de détection (2), par analyse de l'image du masque codé et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs (6, 8), et - déterminer le spectre énergétique et le flux de la source (16), par analyse des impulsions électriques sélectionnées et des première et deuxième images, sur la dynamique spectrale totale du dispositif. 5. Dispositif selon l'une quelconque des revendications 1 à 3, dans lequel l'ensemble de détection (2) est orienté vers une région d'intérêt, cet ensemble étant ainsi orienté suivant une direction, et les moyens électroniques de traitement (4) sont prévus pour : reconstruire la direction de la source (16) dans la région d'intérêt, par décodage du motif prédéfini du masque codé (10), projeté sur l'un au moins des premier et deuxième détecteurs (6, 8), - déterminer la distance de la source (16) par rapport à l'ensemble de détection (2), par analyse de l'image du masque codé (10) et du grandissement de cette image sur l'un au moins des premier et deuxième détecteurs (6, 8), - déterminer la distance de la source {16) par rapport à l'ensemble de détection (2), par analyse des impulsions électriques sélectionnées, et - déterminer le spectre énergétique et le flux de la source (16), par analyse des impulsions électriques sélectionnées et des première et deuxième images, sur la dynamique spectrale totale du dispositif. 6. Dispositif selon l'une quelconque des revendications 1 à 5, dans lequel le masque codé (10) comporte un ensemble de pavés, chacun des premier et deuxième détecteurs (6, 8) comporte un ensemble de pixels et la taille des plus petits pavés du masque est au moins égale à trois fois la taille des pixels des premiers et deuxième détecteurs. 7. Dispositif selon l'une quelconque des revendications 1 à 6, dans lequel les moyens électroniques de traitement (4) sont en outre prévus pour calculer la proportion d'événements fortuits qui sont présents dans des données constituées par les première et deuxième images et les impulsions électriques sélectionnées, pour générer des pseudoévénements fortuits afin de les soustraire aux données. 8. Dispositif selon l'une quelconque des revendications 1 à 7, comprenant en outre des moyens (14) de limitation de bruit en dehors du champ visuel du dispositif. 9. Dispositif selon la revendication 8, dans lequel les moyens de limitation de bruit comprennent des parois (14) en plomb ou en tantale. 10. Dispositif selon l'une quelconque des revendications 1 à 9, dans lequel les premier et deuxième détecteurs (6, 8) sont identiques.