Brevet FR3090901A1 — publié le 20200626

Description Titre de l'invention : PROCEDE DE COMPTAGE DE PHOTONS ASSOCIES A UNE SOURCE DE RAYONNEMENT DOMAINE DE L'INVENTION ET ETAT DE LA TECHNIQUE [0001] L'invention s'inscrit dans le domaine (le l'acquisition et du traitement (le données spectrométriques pour la caractérisation et la localisation de sources d'émission de photons, et notamment de sources radioactives ou de rayons X. [0002] On connaît des techniques de localisation en temps réel d'une source de photons, par exemple enfouie dans un terrain, par acquisition de données (le spectre d'émission de photons, avec éventuellement application de méthodes d'amplification (les signaux faibles. On utilise à cette fin un détecteur de photons qui peut être fixe ou mobile. Un tel détecteur est typiquement configuré pour détecter des rayons gamma ou des rayons X, et est capable de fournir un taux de comptage de photons reçus par unité de temps, pour plusieurs niveaux d'énergie donnés. Le taux de comptage pour un niveau d'énergie correspond à une quantité d'événements détectés par unité de temps, les événements détectés étant des impacts de photons ayant ledit niveau d'énergie. Les taux de comptage et les niveaux d'énergie associés permettent de caractériser les sources d'énergie présentes dans le terrain. La mesure du taux de comptage de photons permet en effet de calculer l'activité d'une source de photons présente dans le sol et positionnée de sorte à émettre des photons en direction de la tête du détecteur. 100031 Il est connu de quadriller une surface au sol en une pluralité de points de mesure, et d'effectuer des mesures de taux de comptage de photons en chacun desdits points, avec éventuellement une corrélation entre les mesures réalisées en des points voisins, pour repérer des sources de photons dans le sol et réaliser une cartographie. 100041 Cependant, les détecteurs de rayonnement connus n'opèrent pas de discrimination parmi les événements unitaires de détection de photons, aussi dits « coups », entre les coups qui sont liés à la présence d'une source radioactive, et les coups liés à l'activité naturelle du terrain, formant un « bruit de fond ». Typiquement, sur un spectre d'énergie acquis en un point de mesure au sol, un pic de rayonnement peut être identifié pour une raie d'énergie à laquelle émet une source, mais les coups détectés pour ladite raie d'énergie et des raies d'énergie voisines sur le spectre ne sont pas uniquement dus à la source. 100051 La détermination ultérieure de l'activité de la source s'en trouve faussée. PRESENTATION GENERALE DE L'INVENTION 100061 11 existe ainsi un besoin pour un procédé de détermination de taux de comptage de photons qui permette de mieux isoler l'activité d'une source de photons émettant à un 2 niveau d'énergie précis, par rapport au bruit de fond qui n'est pas caractéristique de la présence d'une source, et ce même lorsque la contribution de la source est faible par rapport au bruit de fond. [0007] Il existe également un besoin pour un procédé de détermination du taux de comptage qui soit rapide et précis. [0008] A cc titre, l'invention fournit un procédé de détermination d'une part de taux de comptage de photons attribuable à une source de photons, à partir d'un spectre de rayonnement comprenant un taux de comptage de photons associé à chaque niveau d'une pluralité de niveaux d'énergie, [0009] le procédé comportant les étapes suivantes : [0010] détection d'un maximum local de taux de comptage, ledit maximum constituant un pic de rayonnement, [0011] détermination d'un premier niveau d'énergie et d'un deuxième niveau d'énergie, ces deux niveaux délimitant une première fenêtre d'énergie qui comprend le pic de rayonnement, [0012] détermination d'un niveau haut supérieur au deuxième niveau, délimitant avec le deuxième niveau une fenêtre de haute énergie, le niveau haut étant déterminé en fonction d'une région d'intérêt de la première fenêtre d'énergie, [0013] et, de préférence, détermination d'un niveau bas inférieur au premier niveau, dé- limitant avec le premier niveau une fenêtre de basse énergie, le niveau bas étant alors déterminé en fonction de la région d'intérêt de la première fenêtre d'énergie, [0014] calcul d'une part du taux de comptage total de photons dans la première fenêtre d'énergie qui est attribuable à une source de photons émettant selon le pic de rayonnement, en retranchant dudit taux de comptage total le taux de comptage total dans la fenêtre de haute énergie, et éventuellement le taux de comptage total dans la fenêtre de basse énergie. [0015] Un avantage du procédé de l'invention est de permettre la prise en compte du bruit de fond. Le nombre de coups détectés dus au bruit de fond pour le niveau d'énergie correspondant au pic de rayonnement identifié dans le spectre est obtenu à partir du nombre de coups totaux détectés pour la fenêtre de haute énergie, et le cas échéant pour la fenêtre de basse énergie, lesdites fenêtres étant adjacentes à la première fenêtre d'énergie contenant le pic. [0016] On peut ainsi obtenir un taux de comptage de photons qui est réellement attribuable la source, indépendamment du bruit de fond. [0017] Les développements techniques récents ayant fourni des détecteurs de rayonnement ayant une sensibilité accrue aux « coups », il est possible de balayer plus rapidement une pluralité de points d'un terrain pour rechercher des sources, sans perdre en précision pour la détection du taux de comptage. 3 [0018] Le procédé de l'invention présente en outre l'avantage d'être de mise en oeuvre simple par une unité de traitement. Le procédé peut s'intégrer dans un contexte de localisation en temps réel de sources radioactives dans une zone et de cartographie de la radioactivité environnementale, pour une situation de crise ou non. [0019] Un procédé selon l'invention peut présenter les caractéristiques suivantes, non li- mitatives, prises seules ou dans l'une quelconque des combinaisons techniquement possibles : [0020] - le procédé comprend la détermination du niveau bas et de la fenêtre de basse énergie, l'étape de calcul comprenant la soustraction du taux de comptage total dans la fenêtre de basse énergie au taux de comptage total dans la première fenêtre d'énergie, [0021] - la région d'intérêt de la première fenêtre d'énergie est définie à l'aide de la largeur à mi-hauteur du pic de rayonnement détecté ; [0022] - la région d'intérêt de la fenêtre de haute énergie est égale à la région d'intérêt de la fenêtre de basse énergie ; [0023] - les régions d'intérêt des fenêtres de haute et de basse énergie sont toutes deux égales à la moitié de la région d'intérêt de la première fenêtre d'énergie ; [0024] - pour la détermination de la part du taux de comptage total attribuable à la source, on retranche du taux de comptage total les taux de comptage totaux dans la fenêtre de haute énergie et dans la fenêtre de basse énergie ; [0025] - lors de l'étape de détection, deux pics de rayonnement associés à deux niveaux d'énergie distincts sont détectés dans le spectre, la première fenêtre d'énergie comprenant les deux dits pics ; [0026] - la région d'intérêt de la fenêtre de haute énergie est égale à la région d'intérêt de la première fenêtre d'énergie ; [0027] - pour la détermination de la part du taux de comptage attribuable à la source, on retranche uniquement le taux de comptage total dans la fenêtre de haute énergie ; [0028] - le procédé comprend une étape ultérieure de détermination d'une activité d'une source de photons qui serait située aux coordonnées géographiques du point de mesure ; [0029] - le procédé comprend en outre une étape préliminaire d'acquisition du spectre de rayonnement au point de mesure. [0030] L'invention vise en second lieu un système d'identification de source de photons, comprenant : [0031] un détecteur de rayonnement configuré pour mesurer un spectre de rayonnement, [0032] une unité de traitement et une mémoire associée à l'unité de traitement, le détecteur étant configuré pour transmettre des mesures de spectre à l'unité de traitement, [0033] l'unité de traitement étant configurée pour déterminer, à partir d'un spectre acquis par le détecteur, une part de taux de comptage de photons attribuable à une source de photons par la mise en oeuvre d'un procédé de détermination de taux de comptage comme défini ci-avant. [0034] Le système de l'invention peut, de manière optionnelle et non limitative, présenter les caractéristiques suivantes prises seules ou en combinaison : [0035] - le système est placé sur un vecteur de transport, le système étant adapté pour être déplacé sur la surface d'un terrain, le détecteur de rayonnement étant adapté pour acquérir des spectres de rayonnement à une pluralité de points de la surface ; [0036] - le système comprend en outre une structure de support du détecteur, et une unité de propulsion adaptée pour suspendre la structure de support au-dessus de la surface d'un terrain, le détecteur de rayonnement étant adapté pour acquérir des spectres de rayonnement à une pluralité de points de la surface ; [0037] - le détecteur de rayonnement est non colli maté ; [0038] - le détecteur de rayonnement est configure pour acquérir des spectres de rayonnement gamma. [0039] Selon un autre aspect, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour la mise en oeuvre du procédé ci-avant lorsque lesdites instructions sont exécutées par une unité de traitement. PRESENTATION GENERALE DES FIGURES [0040] D'autres caractéristiques, buts et avantages de l'invention ressortiront de la des- cription qui suit, qui est donnée à titre d'exemple non limitatif, accompagnée des dessins annexés suivants parmi lesquels : [0041] [fig.l] [0042] La Figure 1 représente schématiquement un système d'identification de source de photons de l'invention. [0043] [fig.2] [0044] La Figure 2 illustre les étapes d'un procédé de détermination de taux de comptage de photons selon un mode de réalisation. [0045] [fig.3a] [0046] [fig.3b] [0047] [fig.3c] [0048] Les Figures 3a, 3b et 3c représentent des exemples théoriques de spectres de rayonnement (coups par seconde en fonction du niveau d'énergie), correspondant respectivement à un bruit de fond, à l'énergie émise par un radionucléide enfoui au sol, et au cumul de l'énergie dudit bruit et dudit radionucléide. [0049] [fig.4a] [0050] La Figure 4a représente un premier spectre de rayonnement gamma réel (coups par seconde en fonction du niveau d'énergie) acquis sur une durée de dix minutes. [0051] [fig.4b] 100521 1 fie.4c1 [0053] [fig.4d] [0054] Les Figures 4b, 4c et 4d illustrent la mise en oeuvre de certaines étapes successives d'un procédé de détermination de taux de comptage selon la Figure 2 à partir d'un deuxième spectre de rayonnement gamma acquis sur une durée d'une seconde. [0055] [fig.5.1] [0056] [fig.511] [0057] [fig.5c] [0058] Les Figures 5a, 5b et 5c illustrent la mise en oeuvre de certaines étapes successives pour la détermination d'un taux de comptage à partir d'un troisième spectre de rayonnement gamma réel acquis sur une durée de dix minutes. DESCRIPTION DETAILLEE DE MODES DE REALISATION [0059] Dans tout cc qui suit, le terme « source de photons » peut désigner une source ra- dioactive, émettant par désexcitation atomique des rayonnements gamma et de l'énergie, ou une source émettant un rayonnement infrarouge, ou encore des rayons X. On se place dans la suite dans le cas d'une source radioactive émettant des photons gamma. [0060] Lc « taux de comptage » ou « taux de comptage de photons » désigne un nombre d'événements de détection, par un détecteur de rayonnement, de photons gamma dans une certaine raie d'énergie, par unité de temps. On parle usuellement de « coups » pour désigner les événements de détection. Le taux de comptage peut par exemple être donné en coups par seconde. [0061] On parle de « région d'intérêt (ou ROI) » pour désigner la largeur d'une fenêtre d'énergie, une fenêtre d'énergie pouvant être définie comme un ensemble de raies d'énergie contiguës. On parlera enfin de « pic de rayonnement » pour un maximum local de taux de comptage dans un spectre de rayonnement. [0062] Système pour identifier des sources de photons [0063] En référence à la Figure 1, un système 1 de détection de source de photons comprend un vecteur de transport 10. Le vecteur de transport est ici un engin aéroporté tel qu'un hélicoptère. En alternative, le vecteur de transport est une voiture, un avion, un diable, etc., sur lequel sont montés un ou plusieurs détecteurs de rayonnement. Le vecteur de transport peut parcourir la surface d'un sol, permettant de parcourir plusieurs zones de mesure à la suite à la surface du sol. Le détecteur est ici configuré pour acquérir et enregistrer des spectres de rayonnement gamma. [0064] Dans le cas de la localisation d'une source radioactive, le détecteur est avanta- geusement un détecteur NaI, car il présente une sensibilité importante, ce qui permet d'acquérir des mesures d'une seconde qui sont représentatives de la zone mesurée. De plus aucun refroidissement n'est à prévoir avec ce type de détecteur. 6 [0065] Pour la caractérisation de zones et la mesure d'activités et/ou de localisation d'une source (qui ne seront pas décrites de manière détaillée dans la présente description), la distance entre le détecteur et un point de mesure situé au sol est un paramètre à prendre en compte afin de corriger la mesure et de la ramener à la surface du sol. [0066] Le système 1 comprend en outre une mémoire 13 pour stocker les résultats de mesure. De manière optionnelle, le système 1 comprend également une unité de calcul et de traitement 12. Le détecteur est alors pourvu d'une interface 15 lui permettant de transmettre des données à l'unité de calcul et de traitement. [0067] Le système peut en outre comprendre un dispositif de positionnement, par exemple un système de géo-positionnement par satellite (GPS), connecté à l'unité de calcul et de traitement 12, pour associer aux données spectrométriques acquises en un point de mesure les coordonnées géographiques dudit point. [0068] Enfin, le système peut, de façon optionnelle, être doté d'un capteur de température apte à fournir à l'unité de traitement 12 une mesure de température du détecteur. [0069] En alternative, le système 1 ne comprend pas d'unité de calcul configure° pour réaliser l'analyse des données de mesure. La mémoire 13 du système 1 peut alors être reliée via un réseau de communication à un serveur configure pour réaliser l'analyse des données. [0070] Procédé de détermination du taux de comptage pour une source [0071] On a représenté en Figure 2 les étapes d'un procédé 20 de détermination, à partir d'un spectre de rayonnement acquis à un point de mesure, d'une part de taux de comptage de photons attribuable à une source. [0072] Le procédé 20 est par exemple mis en oeuvre par une unité de traitement du système 1 décrit ci-avant, ladite unité de traitement ayant été préconfigurée avec des instructions de code informatique permettant la mise en oeuvre du procédé, ou bien mis en oeuvre par un serveur distant. [0073] Optionnellement, un utilisateur du système 1 peut, via une interface, paramétrer les analyses réalisées sur chaque spectre unitaire pour cibler une fenêtre d'énergie donnée. Par exemple, si un opérateur recherche une source de césium radioactif, il peut sélectionner une fenêtre d'énergie comprenant la raie d'énergie 661,6 keV, l'énergie de désintégration du césium 137 étant de 661,6 keV. L'opérateur peut sélectionner en outre d'autres paramètres d'analyse tels que d'autres régions d'intérêt, détections de pics, indicateurs de radionucléides artificiels. Il est également possible d'extraire différentes informations liées au détecteur et au GPS telles que l'élévation par rapport au sol, le type de détecteur, la hauteur. [0074] Dans le cadre du procédé 20 de détermination de taux de comptage, on cherche à éliminer la contribution de bruits de deux types : un premier bruit qualifié de « bruit de fond » et un deuxième bruit qualifié de « bruit du radionucléide ». 7 [0075] En effet, le taux de comptage au niveau du pic de rayonnement, dans la fenêtre d'énergie ciblée par l'utilisateur du détecteur, ne correspond en général pas uniquement à des photons émis directement, c'est-à-dire sans diffusion, par une même source à identifier. Un bruit de fond et un bruit de radionucléide se superposent en général à ces photons émis directement. Pour la raie d'énergie correspondant au pic, certains événements correspondent au bruit de fond lié à l'activité naturelle du terrain, d'autres encore peuvent correspondre à d'autres sources émettant à des niveaux d'énergie plus élevés, mais ayant subi une diffusion Compton et étant détectés à des niveaux d'énergie plus faibles que leur énergie d'émission. [0076] Le principe de l'élimination des bruits de mesure est explicité en relation aux Figures 3a à 3c. [0077] Le spectre de la Figure 3a correspond au bruit de fond, le spectre de la Figure 3h correspond au cumul de l'activité liée au pic de rayonnement à environ 520 keV d'un radionucléide et du bruit dudit radionucléide, et le spectre de la Figure 3c correspond à l'addition des deux spectres des Figures 3a et 3h. [0078] On comprendra que le détecteur de rayonnement n'est pas nécessairement configuré pour discriminer le bruit de fond, le bruit du radionucléide et l'activité liée au pic de rayonnement. Le détecteur mesure donc le spectre de la Figure 3c. [0079] Dans ces trois figures, la partie gauche donne un spectre de rayonnement matérialisé par des points successifs. Le nombre de coups par seconde est donné en ordonnée et l'énergie est donnée en keV en abscisse. On a également représenté sur la partie gauche des trois dites figures une courbe de régression superposée au spectre de rayonnement, en pointillés. [0080] La partie droite des trois dites figures représente uniquement lesdites courbes de ré- gression sur lesquelles on a repéré une première fenêtre correspondant à une région d'intérêt, ainsi que des fenêtres de basse et de haute énergie. [0081] Sur la Figure 3a, le bruit de fond DI du spectre Cl correspond à une courbe de ré- gression linéaire. [0082] Ainsi, si la région d'intérêt est L, la somme de l'aire 30 sous la courbe Dl sur la fenêtre de basse énergie d'amplitude L/2 et de l'aire 32 sous la courbe Dl sur la fenêtre de haute énergie d'amplitude L/2 est égale à l'aire 31 sous la courbe Dl dans la région d'intérêt. [0083] Sur la Figure 3b, on peut extraire du spectre C2 une première courbe de régression D2a correspondant au pic et une deuxième courbe de régression D2b con-espondant au bruit du radionucléide. [0084] On note que la courbe D2b présente une symétrie par rapport à l'axe vertical associé à l'énergie du pic (environ 520 keV). Ainsi, l'aire 40 sous la courbe D2b dans la fenêtre de basse énergie est égale à la somme des ah-es 41 et 42 sous la courbe D2b 8 dans la région d'intérêt. [0085] La Figure 3c représente un spectre C3 qui est la somme des spectres Cl et C2 et qui correspond au spectre observé par le détecteur. La courbe de régression D3h est la somme des courbes de régression Dl et D2h. [0086] Au vu de ce qui précède, il est équivalent de retrancher l'aire 51 sous la courbe D3a dans la région d'intérêt, et de retrancher les aires 50 et 52 qui correspondent respectivement aux fenêtres de basse et de haute énergie. [0087] Or, les aires 50 et 52 sont directement observables dans le spectre de rayonnement mesuré, contrairement à l'aire 51. [0088] Au cours du procédé 20, il est ainsi proposé d'obtenir, à partir d'un taux de comptage total sur une région d'intérêt, le taux de comptage attribuable à une source donnant lieu à un pic de rayonnement, en retranchant dudit taux de comptage total des taux de comptage observés sur des fenêtres de haute énergie, et éventuellement de basse énergie. [0089] De retour à la Figure 2, à une étape 100, une mesure de spectre de rayonnement S est réalisée en un point de mesure par le détecteur 11. Les mesures acquises sont transmises à l'unité de traitement 12. [0090] L'acquisition d'un spectre de rayonnement a lieu par exemple lors d'une calibration, afin de valider la fonction de réponse du détecteur aux rayonnements incidents, en fonction de mesures fournies par des irradiations de sources à activités connues. En effet, chaque détecteur possède son efficacité absolue (également appelée rendement du détecteur) et le calcul réalisé pour obtenir l'activité (les radionucléides dépend de ladite efficacité et varie en fonction de l'énergie desdits radionucléides. [0091] En alternative, la mesure est acquise à l'étape 100 lors d'une campagne de mesure, afin de caractériser une source de rayonnement et/ou d'évaluer une distribution de source(s) dans un territoire contaminé. [0092] Les coordonnées géographiques du point de mesure sont également transmises à l'unité de traitement, ces coordonnées pouvant être acquises par un système de géopositionnement. [0093] Ensuite, l'unité de traitement exécute les étapes 200 à 600 ci-après sur la base du spectre S en un point. [0094] A une étape 200, un pic de rayonnement P est déterminé sur le spectre S. Ce pic correspond à une raie d'énergie. Le pic de rayonnement peut être détecté automatiquement, par des méthodes connues, sur le spectre S (par détection d'un maximum local de taux de comptage). En alternative, la raie d'énergie associée au pic de rayonnement peut être prédéterminée, notamment si l'on connaît déjà une nature atomique des sources radioactives recherchées. [0095] On a représenté en Figure 4a un spectre de rayonnement mesuré. Ce spectre re- 9 présente un taux de comptage par photons en coups par seconde en ordonnée, et des niveaux d'énergie en kiloélectronvolt (ou keV) en abscisse. Cc spectre a été acquis sur une durée de dix minutes. [0096] On a représenté en Figure 4b un spectre S de rayonnement acquis sur une durée d'une seconde. Ce spectre présente une statistique de comptage moins importante que le spectre de la Figure 4a. Ce spectre est par exemple acquis par le détecteur 1 représenté en Figure 1. La plage des abscisses est comprise entre 0 et 1000 keV et est plus restreinte que la plage du graphique de la Figure 4a. [0097] Sur le spectre S, on a réalisé une étape 200 de détection de pic. Un pic P est détecté à la raie d'énergie à 513 keV. Cc pic constitue un maximum local de taux de comptage. [0098] A une étape 300, un premier et un deuxième niveau d'énergie sont déterminés, dé- limitant une première fenêtre d'énergie contenant la raie d'énergie correspondant au pic de rayonnement. Ces niveaux d'énergie sont déterminés pour couvrir une proportion importante des événements de détection liés à la présence de la source à identifier. Par exemple, la région d'intérêt de la première fenêtre d'énergie est définie à l'aide de la largeur à mi-hauteur du pic P de rayonnement détecté. [0099] Sur la Figure 4c, on a représenté une première fenêtre d'énergie Bp autour du pic P déterminée à l'issue de l'étape 300, la première fenêtre étant délimitée par les niveaux d'énergie NI et Nn. La largeur pour le pic P est de 109 keV. Les niveaux d'énergie NI et N2 sont pris respectivement à 458.5 keV et 567.5 keV. [0100] La fenêtre Bp correspond, de manière avantageuse, à une région d'intérêt égale à la largeur à mi-hauteur du pic P. [0101] La quantité d'événements détectés dans la fenêtre d'énergie ainsi déterminée est utile pour déterminer l'activité de la source à identifier. [0102] Cependant, tous les événements ne sont pas attribuables directement à ladite source, selon les principes décrits ci-avant en relation aux Figures 3a à 3c. [0103] Pour évaluer ultérieurement l'activité de la source à identifier, il est souhaitable de donner une estimation de la part du taux de comptage, sur la première fenêtre d'énergie Bp, qui est attribuable à la source prise isolément, en éliminant les contributions d'autres sources et/ou d'un bruit de fond. [0104] A une étape 400, l'unité de traitement détermine un troisième niveau d'énergie N1 formant avec le niveau NI une fenêtre de basse énergie BE correspondant à des énergies plus basses que la raie d'énergie du pic P, ainsi qu'un quatrième niveau N4 formant avec le niveau N, une fenêtre de haute énergie BH correspondant à des énergies plus hautes que la raie d'énergie du pic P. La fenêtre Bp est donc encadrée par les fenêtres B B et BH. [0105] Les niveaux N3 et N4 sont déterminés en fonction de la région d'intérêt de la fenêtre principale Bp. [0106] De manière préférentielle, ces niveaux sont déterminés de sorte que la région d'intérêt des fenêtres BE et B' corresponde à une proportion prédéterminée de la région d'intérêt de la fenêtre B'. [0107] Par exemple, et notamment dans le cas d'un pic P qui n'est pas situé à proximité d'autres pics de rayonnement dans le spectre S, les fenêtres BE et BE sont prises avec une région d'intérêt égale à la moitié de la région d'intérêt de la fenêtre B'. [0108] On a représenté des fenêtres BE et BE déterminées à l'issue de l'étape 400 sur la Figure 4d. Les coups détectés sur la fenêtre BE correspondent à un bruit de fond, pouvant comprendre des photons ayant subi une diffusion Compton et ayant perdu de l'énergie après leur émission. [0109] A une étape 500 du procédé 20, on utilise le bruit de fond observé sur les fenêtres BE et BE pour interpoler le bruit sur la fenêtre B' contenant le pic P. On utilise le principe explicité ci-avant en relation aux Figures 3a à 3c, pour supprimer les bruits de mesure. [0110] On calcule des sommes de taux de comptage respectivement sur les fenêtres BE, B, et BE. Ces sommes correspondent respectivement à l'intégrale de la courbe de taux de comptage de la Figure 4d sur les ensembles [N3 ; Ni], [N1 ; N2] et [N2 ; N4]. [0111] Le taux de comptage de photons attribuable à la source liée au pic de rayonnement P sur la fenêtre B, est alors déterminé comme étant égal au taux de comptage total sur la fenêtre B, auquel on soustrait le taux de comptage total sur la fenêtre BE et le taux de comptage total sur la fenêtre BH. [0112] Dans cet exemple, en ajoutant les régions d'intérêt des fenêtres BE et BH, on obtient une région d'intérêt de la bande B. [0113] Ainsi, le taux de comptage attribuable à la source obtenu à l'issue de l'étape 500 correspond à l'intégrale sur la fenêtre B, à laquelle on soustrait la moyenne des bruits de fond sur les fenêtres de haute énergie B' et de basse énergie B. [0114] Le résultat ainsi obtenu peut être avantageusement utilisé à une étape ultérieure 600 pour calculer l'activité de la source à identifier et qui est liée au pic de rayonnement P. En effet, l'activité de la source est fonction du taux de comptage sur une fenêtre d'émission de la source, et d'autres paramètres tels que la distribution de radionucléides dans le terrain, l'efficacité intrinsèque du détecteur, des informations sur la topographie du terrain... [0115] Le taux de comptage de photons attribuable à la source à identifier et l'activité de ladite source sont par exemple reliés selon l'expression mathématique suivante : [0116] A = N / (Ps*T), [0117] avec A l'activité du radionucléide en becquerel (ou en becquerel par unité de surface, ou en becquerel par unité de masse), N le nombre de coups directement attribuable au radionucléide déterminé à l'issue du procédé décrit ci-avant, I la probabilité d'émission de la raie énergétique associée au radionucléide en %, T le temps d'acquisition de la 11 mesure en secondes (600 secondes dans l'exemple qui précède) et E l'efficacité intrinsèque du détecteur à un niveau d'énergie correspondant à la raie du radionucléide. [0118] Il est possible d'effectuer des traitements ultérieurs au calcul de cette activité pour déterminer une localisation de source radioactive dans le terrain. Des modèles prédéterminés pour le calcul de probabilité de présence d'une source en un point, comme ceux de la demande de brevet FR 2 991 781 au nom de la Demanderesse, permettent d'obtenir des probabilités d'angle d'incidence d'un photon détecté à un point donné par un détecteur non collimaté. [0119] On a représenté en Figures Sa et 5h un spectre S' de rayonnement obtenu en un point, auquel on applique une variante de réalisation de procédé de détermination de part de taux de comptage attribuable à une source. Les mêmes références sont utilisées sur ces figures et sur les Figures 3h pour des éléments similaires. [0120] Sur le spectre S', on note la présence de deux pics de rayonnement Pl et P2 rapprochés. Les deux pics correspondent à deux raies d'énergie séparées entre elles de seulement 36 keV. L'application des étapes 200 et 300 du procédé 20 donne, dans le cas d'une fenêtre BE, de région d'intérêt égale à la largeur à mi-hauteur pour le pic P2, une fenêtre d'énergie dans laquelle est également présent le pic Pi. [0121] Des fenêtres de basse et haute énergie qui seraient définies comme précédemment, avec une région d'intérêt égale à la moitié de la région d'intérêt de la fenêtre B., sont susceptibles de donner un taux de comptage de la fenêtre de basse énergie B. très proche du taux de comptage moyen sur la fenêtre B. En retranchant le taux de comptage de la fenêtre BR, on risquerait d'ôter une proportion trop importante et de sous-estimer le taux de comptage attribuable aux sources à identifier. [0122] Ainsi, dans ce cas avec deux pics, à l'étape 400 on détermine, comme représenté sur la Figure 4a, le niveau N4 de sorte que la fenêtre de haute énergie B, soit de la même région d'intérêt que la fenêtre Bp. Il n'est pas nécessaire de déterminer une fenêtre de basse énergie. Ensuite, à l'étape 500, pour obtenir le taux de comptage attribuable à la source pour les pics PI et P2, on retranche à la somme de taux de comptage sur la fenêtre Bp la somme de taux de comptage sur la fenêtre BH, comme représenté en Figure 4b. [0123] Ce mode de réalisation présente, pour le cas de deux pics de rayonnement rapprochés, l'avantage de permettre une estimation de taux de comptage qui ne surestime pas le bruit de fond. Cette variante de réalisation peut également être utilisée pour le cas de trois pics de rayonnement rapprochés, ou pour plus de trois pics. 12