Brevet FR2620829A1 — publié le 19890324

SYSTEME DE SUPPRESSION DU BRUIT ET DE SES VARIATIONS POUR LA DETECTION D'UN SIGNAL PUR DANS UN SIGNAL DISCRET BRUITE MESURE DESCRIPTION La présente invention a pour objet un système de suppression du bruit et de ses variations pour la détection d'un signal pur dans un signal discret bruité mesuré. Ce système évalue l'information apportée par le signal mesuré. Cette information se traduit par une probabilité de fausse détection, c'est-à-dire une probabilité de détecter un signal pur là o il n'y en a pas. L'information peut être comparée à un seuil si l'on désire uniquement détecter la présence ou l'absence d'un signal pur, ce seuil correspondant à un taux de fausse détection prédéterminé; l'information peut aussi être évaluée si on désire connaître le degré de certitude de l'existence d'un signal pur. L'invention concerne la détection des signaux discrets, c'est-à-dire les signaux produits par des événements discontinus. L'invention s'applique notamment à la détection des rayonnements corpusculaires, tels que les rayons neutroniques ou les rayons gamma. La détection d'un signal pur dans un signal bruité mesuré suppose que l'on peut supprimer le bruit dans ledit signal mesuré. Deux cas se présentent: le taux moyen de bruit est constant et le taux moyen de bruit est fluctuant. L'invention concerne un système de détection lorsque le taux moyen de bruit est fluctuant. Un taux moyen de bruit constant n'étant qu'un cas particulier d'un taux moyen de bruit fluctuant, le système de l'invention peut a fortiori être utilisé pour la détection lorsque le taux moyen de bruit est constant. La détection d'un rayonnement corpusculaire particulier, par exemple nucléaire, se fait en présence d'un bruit corpusculaire ayant plusieurs causes. Dans le cas d'un rayonnement nucléaire, une partie du bruit est due au rayonnement cosmique et l'autre partie est due à la radioactivité naturelle 2 2620829 des terrains et des roches. Il existe aussi un bruit dû au détecteur utilisé. Le rayonnement cosmique peut varier au cours du temps et le bruit corpusculaire dû aux roches et au terrain est essentiellement variable, selon le lieu o se fait la mesure. Un rayonnement particulier ne peut donc être détecté que si l'on connait le bruit, ou bien si on en a une estimation suffisante. Le stockage des informations est effectué pendant une durée de temps limitée AT. Ce stockage des informations se fait de façon discrète: lorsqu'une particule est détectée, le contenu d'une mémoire est augmenté d'une unité. En général, un ou plusieurs paramètres définissent l'événement nucléaire détecté. Dans le cas o cet événement ne dépend que d'un paramètre unique, celuici peut être par exemple l'énergie; c'est le cas typique de la détection d'un rayonnement gamma à l'aide d'un germanium intrinsèque ou d'un scintillateur NaI associé à un photomultiplicateur. Ce paramètre unique peut aussi être le temps, par exemple lorsqu'on mesure le temps de vol c. e particule. La condition de détection peut être également une fonction logique complexe dépendant de l'état de plusieurs détecteurs durant un intervalle de temps très court, en général de l'ordre de quelques dizaines de nanosecondes, notamment dans le cas de la détection de particules dans une chambre à fils. La difficulté de la détection d'un signal discret bruité réside dans le fait que le bruit de fond a un niveau essentiellement variable, en fonction du temps et en fonction du lien o se fait la détection. Le demandeur a décrit dans le document FR-A-2579767 un système pour la détection d'un signal pur dans un signal discret bruité mesuré. Dans ce système, on utilise une première fenêtre énergétique, appelée voie de signal, correspondant à l'énergie des événements à détecter (rayonnement corpusculaire particulier ou autre) et une seconde fenêtre énergétique, appelée voie de référence pour évaluer le bruit de fond. Pendant des durées AT égales, on compte le nombre N 3 2620829 d'événements détectés dans la voie de signal et le nombre N d'événements détectés dans la voie de référence. Les N événements détectés dans la voie de signal comprennent N événements dus au rayonnement corpusculaire à détecter et des événéments dus au bruit de fond. Le critère de détection d'un signal pur dans la voie de signal, c'est-à-dire N >0, est défini par l'inégalité N >N 3/, o X est le rapport des taux moyens de bruit M et M 2 3 3 2 respectivement dans la voie de référence et dans la voie de signal. Lorsque le critère de détection est satisfait, il existe une information I sur l'existence d'un signal pur dans le signal mesuré, et cette information est d'autant plus importante que la différence N2-N3/ est élevée. Cette information I peut être quantifiée, notamment selon l'expression I=log(PF/(1-PF)) dans laquelle PF est la probabilité de fausse détection, c'est-à-dire la probabilité de détecter un événement du au bruit et que l'on attribue au rayonnement corpusculaire particulier. Le demandeur a constaté que le système décrit dans le document FR-A2579767 ne pouvait pas être utilisé efficacement darns certaines situations. En particulier la détection est faussée si une source de rayonnement se superpose au bruit dans la voie de référence, car le bruit de fond est alors évalué de manière incorrecte. Une autre limitation apparait dans le cas de la détection d'un rayonnement gamma pour mesurer la radioactivité d'un matériau particulier. En effet, le détecteur est souvent constitué d'un cristal de NaI dont la faible résolution en énergie peut dans certains cas provoquer un recouvrement partiel de la voie de signal et de la voie de référence, alors qu'elles sont théoriquement disjointes. L'invention a pour but de remédier aux limitations du système décrit dans le document FR-A-2579767. Ce but est atteint par l'utilisation d'un système 4 2620829 comprenant une voie unique dans laquelle sont mesurés successivement le bruit de fond et le signal complet (signal pur et bruit de fond). De manière précise, l'invention a pour objet un système de suppression du bruit et de ses variations pour la détection d'un signal pur dans un signal bruité discret mesuré comprenant en série un détecteur pour convertir un rayonnement reçu en impulsions électriques, un moyen de mise en forme, un convertisseur analogique-numérique, un moyen de traitement comprenant un processeur, au moins une mémoire et un circuit d'entréesortie. Dans ce système, le processeur mesure un taux de comptage N dans une fenêtre énergétique prédéterminée pendant une durée AT; mesure un taux de comptage N dans la même s5 3T fenêtre énergétique pendant une autre durée AT, égale à q. 4Ts, o q est un nombre réel positif, pendant laquelle aucun événement du au signal pur attendu ne se produit; commèe l'information I, correspondant à la probabilité de fausse détection associée aux comptages N2, N et au coefficient q, à au moins un seuil I, chaque seuil I correspondant à un taux de m m fausse détection prédéterminé; et commande le moyen de signalisation en fonction de la valeur de l'information I par rapport aux seuils I. m Ce dispositif est utilisable efficacement lorsque les durées de comptage sont faibles devant la constante de temps de variation du bruit de fond en un lieu donné, ou si les comptages se font en des emplacements o le bruit de fond moyen est identique. Le dispositif de l'invention est ainsi particulièrement bien adapté à la détection de rayonnements gamma ou neutroniques au moyen de détecteurs portatifs déplacés à pied, et à la détection de sources localisées. La valeur du coefficient q exprimant le rapport entre les durées àT et LT peut être choisie de manière quelconque. C db s Cependant, une valeur non-entière entraÂne une complication du 2620829 traitement qui est en général inutile. Dans le cas de la détection d'un rayonnement par comptages successifs de durées identiques 6T1,..., T,..., on peut prendre le comptage fait pendant la durée &T i1 comme comptage de bruit N associé au comptage N fait pendant la durée 3 2 AT.. Dans ce cas, on a q=l. On peut également choisir comme i comptage N la somme de q comptages précédant le comptage N dont 3 2 on veut savoir s'il contient un signal pur. Dans ce cas, q est un entier quelconque; en pratique, sa valeur est de quelques unités. Il convient de noter que le procédé de l'invention consistant à prendre pour N la somme des q comptages précédant le comptage N est totalement différent du procédé classique consistant à essayer d'estimer le bruit en prenant la valeur moyenne des q comptages précédant le comptage N que l'on teste. Dans le procédé classique, on se contente de comparer le comptage N au bruit estimé, obtenu par une moyenne glissante sur q valeurs, o q est de l'ordre de quelques dizaines d'unités. L'information I apportée par la mesure peut être exploitée de deux façons:elle peut être simplement comparée au seuil I pour discriminer la présence ou l'absence d'un signal m pur; on peut aussi mesurer la différence I-I pour évaluer le m degré de certitude de la présence d'un signal pur. La comparaison entre les informations I et I peut être m faite indirectement sur la probabilité de fausse détection et le taux de fausse détection associés à ces informations. L'information I est cependant utilisée de préférence pour l'exploitation, par exemple par visualisation d'un signal représentant l'information en fonction des comptages. Selon l'invention, la détection d'un signal pur consiste à comparer l'information I apportée par les comptages N et N pour la valeur q à au moins un seuil d'information I 3 m prédéterminé. Chaque seuil est défini par un taux de fausse détection TF. Selon une première variante préférée du système de 6 2620829 l'invention, les comparaisons entre l'information I et chaque seuil I sont remplacées par des comparaisons entre la m probabilité de fausse détection PF, associée aux comptages N et N pour la valeur q, et les taux de fausse détection TF associés aux seuils I. m De préférence, l'information I est définie par la relation I=-log(PF/1-PF)); on définit de la même manière I par m I =-log(TF/1-TF)). m Selon une seconde variante préférée du système de l'invention, dans le cas o on connait une relation entre le taux moyen de signal pur M, le taux moyen de bruit M pendant les I 2 durées de comptage T, le taux moyen de bruit M pendant les s 3 durées de comptage AT, et le taux de fausse détection TF, on compare ledit taux moyen de signal pur M au taux estimé par N - i 2 N3/q. Cette comparaison produit le même résultat que la comparaison directe entre la probabilité de fausse détection PF et le taux de fausse détection TF, mais présente souvent l'avantage d'être plus facile à effectuer en pratique. Dans la variante qui vient d'être décrite, on suppose connue une relation du type f(TF, M, M2, M)=0. Compte tenu de 1 2 3 q=M3/M2, ceci revient à connaître une relation du type g(TF, M, M, q)=0. Pour les comptages mesurés N et N3, cette relation ? 2 3' devient g(PF, N -N /q, N /q, q)=0. L'utilisation pratique de la relation g=O pour décider si un signal pur est présent peut conduire, dans certains cas, à des calculs longs et compliqués. La durée du calcul direct de la probabilité de fausse détection peut notamment être fonction de la valeur des comptages mesurés. La durée de ce calcul peut être en particulier incompatible avec une détection en temps réel, si la durée de chaque comptage est de une seconde par exemple. Selon une caractéristique préférée du système de détection de l'invention, la fonction g est tabulée. Ceci permet de calculer rapidement le taux moyen de signal pur M pour le comparer à N2-N3/q, ou la probabilité de fausse détection PF pour comparer à N -N /q, ou la probabilité de fausse détection PF pour 2 3 7 2620829 la comparer au taux de fausse détection TF. La tabulation de la fonction g présente l'avantage important que la durée du calcul du taux moyen de signal pur M1, o de la probabilité de fausse détection PF, est pratiquement indépendante des valeurs de comptage N et N mesurées. 2 3 La fonction g peut être tabulée pour plusieurs valeurs différentes du taux de fausse détection TF. La comparaison de la valeur estimée du signal N2- N /q aux valeurs de M correspondant à chaque taux de fausse déection TF permet de mesurer l'information apportée par les comptages. Dans le cas o la fonction g est tabulée pour une valeur unique du taux de fausse détection TF, la comparaison entre N2-N3/q et M1 situe l'information I par rapport au seuil 23 / I I, mais ne détermine pas la valeur de l'information I. m La mémorisation sous forme de tables de la fonction g est particulièrement facile dans le cas o des relations linéaires peuvent être mises en évidence entre les variables de ladite fonction g ou entre des fonctions simples desdites variables. Les rayonnements corpusculaires, tels que les rayons neutroniques ou les rayons gamma, sont des rayonnements pour lesquels le nombre de particules détectées pendant un intervalle de temps & T est une variable aléatoire dont la densité de probabilité est celle de la loi de Poisson. On constate que, dans ce cas, une relation sensiblement linéaire existe entre Log(N) et Log(N2), o N =N 2-N /q, pour une valeur fixée de q et pour un taux de fausse détection TF fixé. Une autre possibilité de tabulation de la fonction g, qui peut être mise en oeuvre même lorsqu'il n'existe pas de relation linéaire entre les variables de la fonction g consiste à créer, pour chaque valeur du taux de fausse détection TF, une table comprenant un ensemble de couples (N, C(N)), o C(N) est (N22 l'entier le plus proche du taux moyen de bruit M dans le comptage N qui engendre une probabilité de fausse détection PF égale au taux de fausse détection TF. 8 2620829 Ainsi, pour chaque comptage N mesuré, on compare la valeur C(N) lue dans la table à l'estimation N /q du taux moyen 2 3 de bruit dans le comptage N et on décide Qu'un signal pur est présent, avec une probabilité de fausse détection inférieure à TF, si C(N)>N 3/Q. 2 3 Selon une variante, on peut choisir pour C(N) le taux moyen de bruit M égal à q.M. Pour chaque comptage, on compare 3 2 alors C(N) à N 2 3 Les caractéristiques et avantages de l'invention ressortiront mieux de la description qui va suivre, donnée à titre illustratif mais non limitatif, en référence aux dessins annexés, sur lesquels: - la figure 1 représente schématiquement le système de l'invention, - la figure 2 représente graphiquement pour q=1 la relation entre Log(M) et Log(M) pour différentes valeurs du i 2 taux de fausse détection TF, dans le cas o la densité de probabilité du nombre de particules détectées est celle de la Loic de Poisson, - la figure 3 illustre graphiquement la détermination de Log(M) en fonction du comptage N, au moyen des relations 1 2 linéaires mémorisées, et - la figure 4 est un graphe illustrant, pour différentes valeurs du taux de fausse détection prédéterminé TF, la présence ou l'absence de signal pur dans le signal bruité mesuré en fonction du comptage N et du comptage N. 2 3 L'invention s'applique à tous les rayonnements discrets, tels que les rayonnements corpusculaires, et concerne plus précisément l'analyse du signal délivré par un détecteur pour décider, avec une probabilité de fausse détection inférieure à un taux de fausse détection prédéterminé, si le signal analysé contient un signal pur dû à la présence d'une source particulière ou si ce signal se compose uniquement de bruit. A titre d'exemple, le procédé de l'invention s'applique notamment à la détection des rayons gamma et des rayons neutroniques. Le système représenté sur la figure 1 comprend un 9 2620829 détecteur 2, un moyen de mise en forme 4, un convertisseur analogique numérique 6, un moyen de traitement 8 et un moyen de signalisation 10. Le détecteur 2 dépend du rayonnement à détecter. Pour la détection de rayons gamma, ce détecteur peut être un germanium intrinsèque ou un scintillateur NaI associé à un photomultiplicateur; pour la détection de neutrons, le détecteur est par exemple un compteur à He3. Le détecteur 2 délivre une impulsion de tension pour chaque événement détecté. Cette impulsion est reçue par le moyen de mise en forme 4. Il peut comprendre notamment en série un amplificateur, un circuit de mise en forme et un discriminateur. Le convertisseur analogique-numérique 6 est disposé en aval du moyen de mise en forme 4. Il convertit l'amplitude de l'impulsion reçue en un signal binaire, par exemple sur 8 bits dans le cas d'un scintillateur NaI ou 12 bits dans le cas d'un germanium intrinsèque. Ce signal binaire est transmis au moyen de traitement 8 par une voie 12. Le moyen de traitement 8 comprend un processeur 14, au moins une mémoire 16 et un circuit d'entrée-sortie 18. Ces éléments sont reliés par une voie 20 comprenant des lignes d'adresse, des lignes de données et des lignes de commande. Le moyen de signalisation 10 est commandé par le processeur 14 à travers le circuit d'entrée-sortie 18 et une voie 22. Il peut comprendre notamment des moyens d'alerte visuels et sonores. Ceux-ci sont déclenchés lorsqu'un signal pur est détecté, dans le cas o l'information I est comparée à un seuil unique I. Ils délivrent la valeur de l'information I si celle-ci m est estimée par comparaison à plusieurs seuils différents. Le fonctionnement du système est le suivant. Des comptages N et N sont effectués pendant des durées a T 2 3 S et L T indépendantes pour compter le nombre d'impulsions reçues b par le détecteur dans une fenêtre énergétique prédéterminée correspondant à l'énergie du signal à détecter. Le comptage N est fait pendant une durée 4 T pendant 3 b 2620829 laquelle on peut supposer qu'il n'y a pas de signal, le comptage N représente donc uniquement le bruit de fond dans la fenêtre énergétique prédéterminée. Le comptage N, qui sert de référence, peut *tre effectué indifféremment avant ou après le comptage N. Les durées 4T et ET ne sont pas nécessairement b s égales. Il peut être avantageux, dans certains cas, de choisir un rapport q=4T Tb/ T différent de 1. Par exemple, lorsque le b s système est utilisé lors d'une expérience auprès d'un accélérateur de particules, il est important -compte tenu du coût de ce genre d'expériences- d'optimiser le rapport q pour obtenir la quantité d'informations désirée en un temps minimal. En revanche, dans le cas d'une détection séquentielle o les comptages ont lieu de manière répétitive pour des durées consécutives AT1. T2,..., TT,... aucun problème d'optimisation ne se pose plus. La solution la plus simple consiste alors à choisir pour le comptage N_, associé à -- comptage N effectué pendant la durée T, la somme des q 2i comptages effectués pendant les durées AT, AT, ..., i-q i-q+1 AT. Le rapport q est alors un nombre entier. Il est avantageux i -1 de choisir en pratique une valeur comprise entre 1 et 5. Une variante possible, lorsque q=l consiste à prendre comme référence N le comptage fait pendant la durée T pour 3 i-1 tous les comptages N faits pendant les durées T, T i+l, 2 i i+ LT i+2',... tant que l'information apportée par les comptages N2 et le comptage N est supérieure à un seuil prédéterminé, puis à reprendre l'avant-dernier comptage N comme comptage de référence N lorsque l'information devient inférieure au seuil prédéterminé. Pour chaque couple de comptage (N, N), le système 2 3 compare la probabilité de fausse détection PF associée à ces comptages à au moins un taux de fausse détection TF, ou une fonction de ladite probabilité de fausse détection à une fonction des taux de fausse détection prédéterminés, cette fonction pouvant *tre l'information I. l 2620829 On va maintenant décrire, à titre d'exemple, le fonctionnement du système de l'invention pour la détection d'un rayonnement corpusculaire.La densité de probabilité du nombre de particules détectées est celle de la loi de Poisson. La probabilité pour qu'un rayonnement existe, que le stockage N contienne ou non un événement qui lui soit dû est, à M constant J = N -M J N2 M 2 2J Pl e P =e.'> ---- _ P1 / J! J = 0 Par ailleurs, la probabilité pour que M soit compris entre M3 et M3+ dM3 est N3 -M3 (M3) e. N!. d.3 P 3 2N cc N3 o-M3 (M3) t e * N! -.dM3 ! La probabilité d'absence d'un rayonnement pur, ou probabilité de fausse détection, est alors: PF = f p P À {q.N21 o signifie "partie entière de", c'est-à-dire, en utilisant qM 2=M, 12 2620829 CPF N +1 (J+N3 - PF (q33 v PFJ+ J J!.-N z33 J:N2+1 K=J+N3 -(+q)'N 2 3 avec Yt 2 ((1+.N2)K K=D K! Le taux de fausse détection TF s'obtient en remplaçant dans l'expression de PF: N2 par {M1 + {M2 N3 par {M31 o sifnifie "partie entière de". Compte tenu de q=M /M, cette relation est de la forme g(TF, M1, M2, q)=0. Ainsi qu'on l'a indiqué plus haut, lorsque le densité de probabilité du nombre d'événements détectés suit une loi de Poisson, on peut mettre en évidence des relations sensiblement linéaires entre des fonctions simples des variables TF, M1, M et q. Ceci permet de mémoriser sous forme de tables la fonction g et de déterminer pratiquement sans calcul le taux moyen de signal M correspondant à un taux de fausse détection prédéterminé TF, pour un taux de bruit M estimé et un rapport q fixé, et de comparer ledit taux moyen de signal M au taux estimé égal à N -N /q. 1 2 3 La figure 2 est un graphe montrant la relation sensiblement linéaire entre Log(M) et Log(M), pour la valeur I 2 particulière q=1. Ce graphe comprend un ensemble de droites correspondant chacune à une valeur de TF. 13 2620829 La mémorisation de toutes les droites n'est pas nécessaire. On pourra ne mémoriser que les droites correspondant, -3 -4 -5 -6 par exemple, à TF=10, 10, 10 et 10. Pour chacune des droites à mémoriser, deux couples de points au moins doivent être stockés. Ce nombre pourra être augmenté si l'on estime que la linéarité est insuffisante. Pour la détection, conformément à l'invention, on construit une table comprenant k.k2 couples de points (Log(M1), Log(M2)), o k1 est le nombre de valeurs de TF mémorisées et k2 le nombre de points par droite mémorisée. La détection d'un signal pur est alors faite de la manière suivante. Des comptages N2 et N3 sont réalisés pendant des durées àT et &T, o LTb =q. AT, le comptage N étant s b b s 3 fait à un moment o on suppose qu'il n'y a pas de signal pur. On estime ensuite le taux moyen de bruit M pendant la durée NT par le rapport N /q, et on en déduit le taux moyen de s 3 signal pur M pendant la durée AT au moyen de la table. 1 s Pour la valeur TF choisie, le taux moyen de signal pur M s'obtient comme on l'a représenté sur la figure 3, par interpolation des couples (Log(M'1), Log(M')) et (Log(M", 2i Log(M")), et en fonction du taux moyen de bruit M2. 2 2 L'interpolation peut être remplacée par une extrapolation si Log(M) n'est pas compris entre Log(M') et Log(M"). 2 2 2 Il ne reste alors plus qu'à comparer la valeur M à son i estimation égale à N2-N /q. Si M est supérieur à N -N /q, on 23 I 2 3 considère que le comptage N ne contient que des événements dus au bruit; dans le cas contraire, on considère qu'un signal pur a été détecté. La comparaison entre la valeur M du taux moyen de signal donnée par la table et la valeur estimée N -N /q est 2 3 identique à la comparaison entre le taux de fausse détection TF et la probabilité de fausse détection PF, ou entre le seuil d'information I et l'information I apportée par les comptages. m La comparaison entre la valeur estimée N -N /q et une 2 3 unique valeur de M 1, définie par un taux de fausse détection 14 2620829O prédéterminé TF, réalise une détection par seuil: un signal pur est présent si le seuil est franchi, absence de signal pur si le seuil n'est Pas franchi. Comme on l'a déjà indiqué, l'information peut tre évaluée de manière plus quantitative. Cette évaluation de l'information utilise la suite des opérations qui viennent d'être décrites, celles-ci étant répétées pour plusieurs valeurs de TF, -3 -4 -5 -6 par exemple 10, 10, 10, 10. Le système détermine ainsi plusieurs valeurs du taux moyen de signal M, chacune associée à un taux de fausse détection, donc à un niveau d'information. L'information I apportée par les comptages est comprise entre les informations I' et I" correspondant aux taux moyens de bruit M' et M" définis par M' >N et (N /q) I 1 1-2 3 - 1 L'information est définie conformément à l'invention par la relation I=Log(PF/(1-PF)). La probabilité de fausse -3 détection étant en général inférieure à 10, cette relation se traduit en pratique par I=-LogPF. L'information I a donc une valeur usuelle que quelc;s. unités, pour une base de logarithme égale à 10. La représentation graphique de l'information I en fonction du numéro du comptage, c'est-àdire finalement en fonction du temps, constitue alors pour un opérateur une visualisation particulièrement claire et représentative. Le procédé de détection décrit en référence aux figures 2 et 3 ne peut être mis en oeuvre que si des relations linéaires existent entre des fonctions simples des variables TF, M, M et q. Dans le cas contraire, la détection peut être faire de la manière suivante. Dans une première étape, on crée pour chaque taux de fausse détection TF choisi une table constituée d'un ensemble de couples (N2, C(N2)), o C(N2) est l'entier le plus proche du taux de bruit moyen M du comptage N pour lequel la probabilité de 2 2 fausse détection PF est égale au taux de fausse détection TF. Ensuite, pour chaque couple de comptage (N, N) 2620829 mesuré, on estime le taux moyen de bruit M par la valeur N /q et 2 3 on compare cette estimation à la valeur C(N) Lue dans la table correspondant à TF. Si C(N 2) est supérieur à N3/q, un signal pur est présent dans le comptage N2; si C(N2) est inférieur à N3/q, le comptage N ne contient que du bruit. Suivant une variante, on prend pour C(N2) la valeur q.- M2, c'est-à-dire M3, au lieu de M2. On compare alors C(N) à N3, au lieu de N /q. 2 3 ' 3 On a représenté sur la figure 4, un graphique illustrant les tables construites pour différentes valeurs de TF. La détection selon l'invention consiste à comparer la position du point de coordonnée (N, N) à la courbe du taux de fausse 2 3 détection prédéterminé TF choisi. Si le point est au-dessus de la courbe, le comptage N ne contient que du bruit; dans le cas contraire, le comptage N contient un signal pur. Comme dans le procédé décrit en référence aux figures 2 et 3, le procédédécrit en référence à la figure 4 permet d'évaluer quantitativement l'information I lorsque plusieurs tables, associées chacune à une valeur particulière de TF, sont utilisées. 16 2620829 REVENDICATIONS 1. Système de suppression du bruit et de ses variations pour la détection d'un signal pur dans un signal bruité discret mesuré comprenant en série un détecteur (2) pour convertir un rayonnement recu en impulsions électriques, un moyen de mise en forme (4) desdites impulsions électriques, un convertisseur analogique-numérique, un moyen de traitement (8) et un moyen de signalisation (10), ledit moyen de traitement comprenant un processeur (14), au moins une mémoire (16) et un circuit d'entrée-sortie (18), ledit système étant caractérisé en ce que ledit processeur: -mesure un taux de comptage N dans une fenêtre énergétique prédéterminée pendant une durée T, s - mesure un taux de comptage N dans ladite fenêtre énergétique prédéterminée pendant une autre durée T, égale à q. AT o q est un nombre réel positif, pendant laquelle aL. s événement correspondant au signal pur attendu n'est présent, - compare l'information I, correspondant à la probabilité de fausse détection PF associée aux comptages N2, N - 3 et au coefficient q, à au moins un seuil I, chaque seuil I m m correspondant à un taux de fausse détection prédéterminée, - commande le moyen de signalisation en fonction de la valeur de l'information I par rapport aux seuils I. m 2. Système selon la revendication 1, caractérisé en ce que des comptages sont effectués séquentiellement et de manière continue pendant des durées successives égales tT, AT,... ET,..., et dans lequel le comptage N, associé à un comptage N 1 2 effectué pendant la durée aT, est pris égal à la somme des q comptages effectués pendant les durées tT, T, ..., T i-q i -q+1' AT i-1 3. Système selon la revendication 1, caractérisé en ce que les comptages sont effectués séquentiellement et de manière continue pendant des durées successives égales et en ce qu'un i-1 même comptage de référence N effectué pendant la durée e T 17 2620829 est utilisé pour les comptages N2 effectués pendant les durées T, T, ..., tant que les couples (N2, N) o N est fixé, i i13 3 apportent une information I supérieure à un seuil prédéterminé I. m 4. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que le processeur compare l'information I à un seuil unique I, un signal pur étant détecté si le seuil est m franchi. 5. Système selon l'une quelconque des revendications 1 à 3, caractérisé en ce que le processeur estime l'information I par comparaison à une pluralité de seuils I. m 6. Système selon l'une quelconque des revendications 4 et 5, caractérisé en ce que le processeur compare la probabilité de fausse détection PF, associée à l'information I, au taux de fausse détection TF, associé à chaque seuil I, l'information I m étant définie par I=-Log(PF/(1-PF)) et chaque seuil I par I - m m Log(TF/(1-TF)). 7. Système selon l'une quelconque des revendication 4 à 6 dans lequel on connaÂt une relation g(TF, M, M2, q)=0, o M est le taux moyen de signal pur et M le taux moyen de bruit, caractérisé en ce que la détection est opérée de manière indirecte sur le taux moyen de signal pur, le processeur comparant, pour au moins une valeur du taux de fausse détection TF, la valeur M au taux moyen de signal pur estimé N -N /q, o I 2 3 N3/q est une estimation de M. 8. Système selon la revendication 7, caractérisé en ce que la relation g(TF, M, M, q)=0 est stockée sous forme de table dans la mémoire (16). 9. Système selon la revendication 8, pour la détection d'un signal pour lequel Log(M) est une fonction sensiblement linéaire de Log(M), caractérisé en ce que le processeur: - crée, préalablement à toute mesure, une table comprenant, pour chaque valeur de TF prise parmi un ensemble de valeurs prédéterminées, au moins deux couples (LogM, LogM), - calcule, après chaque mesure et pour au moins une 18 2620829 valeur de TF, Log(M1) associé à l'estimation LogN 3/q) de Log(M2), par interpolation ou extrapolation des valeurs contenues dans la table, - compare le taux moyen de signal brut estimé N2-N /q à au moins un taux moyen de signal pur M calculé. 10. Système selon la revendication 8, caractérisé en ce que le processeur: - crée, préalablement à toute mesure, pour chaque valeur de TF pris parmi un ensemble de valeurs prédéterminées, une table comprenant une suite de couples (N2, C(N)), o C(N2) est l'entier le plus proche du taux moyen de bruit M pendant la durée à T du comptage N qui engendre une probabilité de fausse s 2 détection PF égale au taux de fausse détection TF, Ladite table étant stockée dans la mémoire, et - compare après chaque mesure la valeur C(NI de la table associée au comptage N2, au comptage N 3/q, un signal pur étant détecté si C(N 2) est supérieur à N3/q. 11. Système selon la revendication 8, caractérisé en que le processeur: crée, préalablement à toute mesure, pour chaque valeur de TF prise parmi un ensemble de valeurs prédéteminées, une table comprenant une suite de couples (N, C(N)), o C(N) 2 2' est l'entier le plus proche du taux moyen de bruit M pendant la durée T du comptage N qui engendre une probabilité de fausse b3 détection PF égale au taux de fausse détection TF, ladite table étant stockée dans la mémoire, et - compare après chaque mesure la valeur C(N) de la table, associée au comptage N, au comptage N, un signal pur étant détecté si C(N) est supérieur à N. 2 3 12. Système selon l'une quelconque des revendications 9 à 11, caractérisé en ce que la mémoire contient des tables pour plusieurs valeurs de TF, le processeur calculant le taux moyen de signal pur M pour plusieurs valeurs de TF et évaluant la probabilité de fausse détection associée aux comptages N et N en l'encadrant par les deux valeurs de TF correspondant aux deux valeurs de mi encadrant l'estimation N N /G. 1 2 3