Brevet FR3032281A1 — publié le 20160805

PROCÉDÉ ET SYSTÈME DE STABILISATION D'UNE CHAÎNE DE DÉTECTION SPECTROMÈTRIQUE D'IMPULSIONS PAR CORRÉLATION ENTRE SPECTRES DOMAINE TECHNIQUE Le domaine de l'invention est celui de la mesure de l'activité d'un radioélément. L'invention correspond plus particulièrement à un procédé et à un système visant à s'affranchir de la dérive du gain d'une chaîne de détection spectrométrique d'impulsions issues d'une source de rayonnements ionisants, de lumière ou de vibrations mécaniques lorsque le spectre est invariant et que la mesure utilise l'intensité du signal. Elle trouve notamment application dans le domaine de la médecine nucléaire pour permettre de connaître avec précision l'activité d'un radio élément contenu dans une seringue. ÉTAT DE LA TECHNIQUE ANTÉRIEURE Un système de production de radiotraceurs «à la demande» pour l'imagerie moléculaire TEP (Tomographie par Emission de Positons) intégrant une chaîne complète et automatisée de production et de préparation des radiotraceurs est actuellement en cours de développement. Un tel système est destiné à une installation in situ sans contraintes majeures de radioprotection au sein de centres de recherche et d'hôpitaux. Un tel système nécessite une chaîne de détection de spectrométrie impulsionnelle à même de réaliser un comptage d'impulsions issues d'une source de radioéléments, par exemple pour fournir une information sur l'activité d'un radioélément contenu dans une seringue. Une telle chaîne de détection comprend un scintillateur, un photomultiplicateur et un amplificateur. L'utilisation d'un scintillateur en bromure de lanthane associé à une photodiode à avalanche en silicium permet de fournir un détecteur compact, rapide, disposant d'une bonne résolution énergétique et un bon rendement (densité et numéro atomiques élevés). Mais un tel détecteur est très sensible à la température. A titre d'exemple, un spectre réalisé à 20,6°C présente un pic à 662 keV. Ce pic se retrouve à 780 keV lorsque le spectre est réalisé à 18,4°C, soit une augmentation de 18% quand la température baisse de 2,2°C. On relèvera que les variations de gain et de taux de comptage sont à peu près linéaires en fonction de la température Dès lors, afin de rendre possible une mesure d'activité avec une précision de l'ordre du pourcent, une stabilisation efficace de ce détecteur s'avère nécessaire pour compenser les variations de gain et le rendre notamment indépendant de la température. On connaît du brevet FR 2 626 121 un dispositif électronique de compensation de dérive du gain par égalisation d'un comptage des impulsions au sein de deux fenêtres adjacentes dont une première est délimitée entre un premier et un deuxième seuil homothétiques et dont la deuxième a pour borne inférieure le deuxième seuil. Un tel dispositif ne s'avère pas suffisamment performant, notamment en termes de précision et de robustesse aux bruits, pour répondre aux besoins actuels. EXPOSÉ DE L'INVENTION L'invention vise à répondre au besoin suscité d'une stabilisation efficace d'une chaîne de détection de spectrométrie impulsionnelle, et propose pour ce faire un procédé de compensation de dérive du gain d'une chaîne de détection d'impulsions issues d'une source, comprenant les étapes suivantes : acquisition d'un premier spectre d'impulsions issues de la source ; détermination d'un degré de corrélation du premier spectre avec un spectre de référence fixation d'une variation de gain de la chaîne de détection ; détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence, dans lequel le second spectre correspond directement ou indirectement à une acquisition d'un spectre d'impulsions issues de la source après modification du gain de la chaîne de détection de ladite variation de gain. Certains aspects préférés mais non limitatifs de ce procédé sont les suivants : les étapes de fixation d'une variation de gain et de détermination d'un degré de corrélation du second spectre avec le spectre de référence sont réitérées jusqu'à identifier une variation de gain permettant d'atteindre un degré de corrélation donné, typiquement le degré de corrélation maximal, du second spectre avec le spectre de référence ; l'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend la modification du gain de la chaîne de détection de ladite variation de gain et l'acquisition du second spectre à l'aide de la chaîne de détection à gain modifié ; l'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend la détermination d'un coefficient d'expansion à partir de ladite variation de gain et le calcul du second spectre par compression ou dilation du premier spectre conformément au coefficient d'expansion. L'invention porte également sur un module d'étalonnage d'un système de détection spectroscopique configuré pour mettre en oeuvre le procédé, ainsi que sur un produit programme d'ordinateur apte à mettre en oeuvre le procédé. BRÈVE DESCRIPTION DES DESSINS D'autres aspects, buts, avantages et caractéristiques de l'invention apparaîtront mieux à la lecture de la description détaillée suivante de formes de réalisation préférées de celle-ci, donnée à titre d'exemple non limitatif, et faite en référence aux dessins annexés sur lesquels : les figures la et lb sont des ordinogrammes représentant deux modes de réalisation possible de l'invention ; la figure 2 est un schéma d'un module d'étalonnage configuré pour mettre en oeuvre le procédé selon l'invention ; les figures 3a et 3b illustrent une dilatation, respectivement une compression, de spectre ; les figures 4a, 5a, 6a et 7a illustrent, pour différentes valeurs de gain de la chaîne de détection, le second spectre et le spectre de référence ; les figures 4b, 5b, 6b et 7b illustrent, pour les différentes valeurs de gain de la chaîne de détection, le degré de corrélation entre le second spectre et le spectre de référence ; la figure 8 est un schéma représentant le degré de corrélation en fonction du gain, autour du gain de référence ; la figure 9 est un schéma illustrant la variation de comptage en réponse à une variation brutale de gain ; la figure 10 illustre un filtrage de spectre pouvant être mis en oeuvre dans un mode de réalisation possible de l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention porte sur un procédé de compensation de dérive du gain d'une chaîne de détection d'impulsions issues d'une source. Elle s'étend à un système de détection spectroscopique des impulsions issues de la source, le système comprenant un module d'étalonnage configuré pour mettre en oeuvre le procédé. Elle porte également sur un produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé, lorsque ledit programme est exécuté sur un ordinateur L'invention n'est pas limitée à un système mono-détecteur, mais trouve également application dans des systèmes multi-détecteurs, en permettant la mise en coïncidence des différentes chaînes de détection, ce qui permet d'en réaliser la sommation afin d'être plus robuste au bruit. Le système de détection est plus précisément configuré pour réaliser un comptage des impulsions selon leur amplitude en correspondance avec l'énergie détectée. Le signal incident est collecté sur plusieurs canaux à la fois, chaque canal correspondant à une division de l'échelle d'énergie. Le nombre d'impulsions pour chaque canal est appelé « nombre de coups », le taux de comptage correspondant au nombre de coups par seconde. D'une manière générale, l'invention propose d'avoir recours à un spectre de référence et lors des mesures suivantes de corriger le gain de la chaîne de détection de telle sorte que le spectre acquis vienne se superposer au spectre de référence. Les dérives liées notamment à la température sont ainsi compensées par modification du gain de manière à ce que les deux spectres, réalisés à des températures différentes, se superposent. La corrélation est plus particulièrement utilisée dans le cadre de l'invention pour mettre en coïncidence les deux spectres décalés l'un de l'autre du fait de la dérive du gain de la chaîne de détection. La méthode proposée utilise plus particulièrement la corrélation entre deux spectres homothétiques selon l'axe des amplitudes (correspondant aux énergies). Il est ainsi possible de déterminer le coefficient d'homothétie (gain de la chaîne de détection) le plus probable par le maximum de corrélation avec le coefficient d'homothétie comme paramètre. Généralement, la corrélation est utilisée pour mettre en coïncidence deux signaux semblables mais décalés dans le temps. Il est ainsi possible de déterminer le décalage en temps le plus probable entre eux par la recherche du maximum de corrélation avec le temps comme paramètre. Une mesure d'un degré de corrélation pouvant être utilisée est le coefficient de corrélation Bravais-Pearson. Il a l'avantage d'être adimensionnel (ne dépend pas de l'amplitude des signaux) et de varier entre 0 et 1 (ou -1). Les deux séries de mesures sont d'autant mieux corrélées que le coefficient est loin de zéro. L'expression de ce coefficient rp vaut : r = où o-xy désigne la covariance entre les variables x et y, et o-x ,ay leur écart type. Pour deux séries de mesures x, et y, avec 1 N : r = Avec la moyenne de X: et la moyenne de y: Dans le cadre de l'invention, la même méthode est utilisée avec, pour séries de mesures, les deux spectres et comme paramètre le gain de la chaîne de détection à la place du temps. Selon l'invention, on vient ainsi acquérir un spectre d'impulsions issues d'une source, et modifier le gain de la chaîne de détection pour atteindre un degré de corrélation donné, typiquement le degré de corrélation maximal, du spectre acquis avec un spectre de référence d'impulsions issues de la source. A cet effet, et en référence aux figures 1 et 2 qui en illustrent deux modes de réalisation possibles, le procédé selon l'invention comprend les étapes suivantes : acquisition « ACQ » d'un premier spectre d'impulsions issues de la source ; 15 détermination d'un degré de corrélation « CORR » du premier spectre avec un spectre de référence ; fixation d'une variation de gain « MOD-GAIN » de la chaîne de détection ; détermination d'un degré de corrélation « CORR » d'un second spectre d'impulsions avec le spectre de référence, dans lequel le second spectre 20 correspond directement ou indirectement à une acquisition « ACQ » d'un spectre d'impulsions issues de la source après modification du gain de la chaîne de détection de ladite variation de gain. La comparaison des degrés de corrélation ainsi déterminés permet de déterminer lequel du premier spectre ou du second spectre est statistiquement le plus proche du spectre de référence, et donc de vérifier si l'effet de la variation de gain est positif en ce sens qu'il permet d'augmenter le degré de corrélation avec le spectre de référence. On a représenté sur la figure 2 un module d'étalonnage 1 configuré pour mettre en oeuvre le procédé selon l'invention. Un tel module comprend une unité de corrélation 2 permettant de déterminer le degré de corrélation du premier spectre Si, respectivement du second spectre S2, avec le spectre de référence, une unité de fixation d'une variation de gain 3 et une unité de détermination du second spectre S2 soit par acquisition, soit par calcul comme cela sera décrit par la suite. La fixation de la variation de gain peut consister à réaliser une variation incrémentale, par exemple selon un pas de 2.10-4 pour des spectres de 105 impulsions. La variation de gain peut également être mise en oeuvre selon une méthode dichotomique qui présente l'avantage d'être plus rapide, mais l'inconvénient d'être moins robuste (courbe en cloche). Les figures 4a, 5a, 6a et 7a illustrent, pour différentes valeurs de gain G de la chaîne de détection, le second spectre S2 et le spectre de référence Sref (représentés sous la forme d'un nombre de coups Co par canal Ca) tandis que les figures 4b, 5b, 6b et 7b illustrent, pour ces différentes valeurs de gain G de la chaîne de détection, le degré de corrélation Dc entre le second spectre S2 et le spectre de référence Sref. Le gain G est plus précisément de : 0,5 sur les figures 4a et 4b, 0,9 sur les figures 5a et 5b, 1 sur les figures 6a et 6b, 2 sur les figures 7a et 7b. Dans cet exemple, c'est le gain G égal à 1 qui correspond au plus haut degré de corrélation Dc (figure 6a) et qui permet de mettre en coïncidence les spectres. C'est ce gain qui sera retenu pour les mesures suivantes. Dans un mode de réalisation, les étapes de fixation d'une variation de gain et de détermination d'un degré de corrélation du second spectre avec le spectre de référence sont réitérées jusqu'à identifier une variation de gain permettant d'atteindre un degré de corrélation donné du second spectre avec le spectre de référence. Il s'agit par exemple d'atteindre un degré de corrélation supérieur à un seuil, ou encore d'atteindre le degré de corrélation maximal. La réitération de ces étapes forme une boucle qui, dans un premier mode de réalisation possible illustré par la figure la, est mise en oeuvre en faisant varier le gain et en venant acquérir un nouveau spectre. Cette boucle consiste à répéter les étapes suivantes jusqu'à atteindre le degré de corrélation donné « MAX? »: détermination du degré de corrélation « CORR » du spectre acquis avec le spectre de référence ; modification « MOD-GAIN » du gain de la chaîne de détection; acquisition « ACQ » d'un nouveau spectre (second spectre) d'impulsions issues de la source à l'aide de la chaîne de détection à gain modifié. Une fois le degré de corrélation donné atteint, tel que vérifié à l'opération « MAX? », le procédé s'achève lors d'une étape « END ». Dans ce premier mode de réalisation, le second spectre correspond donc directement à une acquisition « ACQ » d'un spectre d'impulsions issues de la source après modification du gain de la chaîne de détection de ladite variation de gain. L'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend ainsi la modification du gain de la chaîne de détection de ladite variation de gain et l'acquisition d'un nouveau spectre d'impulsions issues de la source à l'aide de la chaîne de détection à gain modifié. Dans un second mode de réalisation possible illustré par la figure lb, la réitération de ces étapes forme une boucle qui est mise en oeuvre en venant compresser/dilater le premier spectre, l'effet de cette compression/dilatation étant semblable à une variation du gain de la chaîne de détection. La modification du gain « MOD- GAIN » est ainsi réalisée après mise en oeuvre d'une boucle consistant à répéter les étapes suivantes : - calcul « CAL)> d'un spectre à partir du spectre acquis lors de l'étape « ACQ » (premier spectre), le spectre calculé correspondant à une compression ou une dilatation du spectre acquis conformément à un coefficient d'expansion dont l'effet est semblable à une modification de gain de la chaîne de détection ; détermination du degré de corrélation « CORR » du spectre calculé avec le spectre de référence ; modification du coefficient d'expansion « MOD-EXP »; jusqu'à atteindre le degré de corrélation donné « MAX ? » entre le spectre calculé et le spectre de référence. Dans ce second mode de réalisation, le second spectre est un spectre calculé pour simuler une modification du gain et qui correspond donc indirectement à une acquisition « ACQ » d'un spectre d'impulsions issues de la source après modification du gain de la chaîne de détection de ladite variation de gain. L'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend ainsi la détermination d'un coefficient d'expansion à partir de ladite variation de gain et le calcul du second spectre par compression ou dilation du premier spectre conformément au coefficient d'expansion. Le coefficient d'expansion correspond typiquement au ratio du nombre de canaux du second spectre et du nombre de canaux du premier spectre. Ainsi à partir d'un spectre d'origine (canaux Ni, N2, etc.), on vient créer un spectre calculé (canaux n1, n2, etc.) qui comporte plus (dilatation) ou moins (compression) de canaux que le spectre d'origine et dont la somme des canaux est identique à celle du spectre d'origine. Les figures 3a et 3b illustrent à ce propos une dilatation de spectre de coefficient d'expansion 0,75, respectivement une compression de coefficient d'expansion 1,25. L'effet d'une telle dilatation/compression est semblable à une variation de gain de la chaîne de détection spectrométrique, comme cela est d'ailleurs est visible sur les figures 4a, 5a, 6a et 7a où une réduction de gain est synonyme de compression et où une augmentation de gain est synonyme de dilatation. Une méthode de calcul du contenu des canaux du second spectre résultant de la compression ou dilation du premier spectre selon un coefficient d'expansion a est la suivante. Le contenu des canaux du spectre dilaté ou compressé est la somme (ou une partie) du contenu des canaux du spectre d'origine. Pour un canal i du spectre calculé, on calcule le début et la fin des canaux concernés dans le spectre d'origine. Le début du canal i se situe à (i-1).a et la fin à i.a (dans l'échelle des canaux du spectre d'origine). Dans le spectre d'origine, le canal de début correspond à la partie entière de (i-1),a et le canal de fin à la partie entière de isa. La somme des canaux se décompose en : - le reste du premier canal, soit le contenu de ce canal multiplié par [1-partie décimale de (i-1).a] ; 10 - la somme des contenus d'un certain nombre de canaux entiers. Ce nombre est déterminé par la différence entre le canal de fin et celui du début ; - Le début du dernier canal, soit le contenu de ce canal multiplié par la partie décimale de /*a. Ainsi, le contenu ni d'un canal d'indice i du second spectre correspond à: 15 - ni = a.ND si = 0 - n = ND. (1 - partie décimale [i - 1] ) + NF. (partie décimale i. a) si à = 1 - ni = ND. (1 - partie décimale [i - 1]. a) + ND+k 20 N. (partie décimale i. a) si A k 2 avec : i un entier compris entre 1 et la partie entière de N*a, où N désigne le nombre de canaux du premier spectre ; Ni le contenu du canal d'indice j du premier spectre ; 25 - D le canal du premier spectre au sein duquel débute le canal d'indice i du second spectre : il s'agit du canal du premier spectre dont l'indice correspond à la partie entière de (i-1).a ; - F le canal du premier spectre au sein duquel s'achève le canal d'indice i du second spectre: il s'agit du canal du premier spectre dont l'indice correspond à 30 la partie entière de i.a ; - =F-D. On a vu précédemment qu'une variation incrémentale sur le gain selon un pas de 2.104 pouvait être mise en oeuvre. Un tel pas est significatif car il entraîne sur le coefficient de corrélation une variation importante, très supérieure au bruit. La figure 8 illustre à ce propos le degré de corrélation en fonction du gain, autour du gain de référence. La forme de cette courbe est très pointue. De sorte qu'une variation de gain de 1/1000 par exemple entraînerait une variation de 2% sur le degré de corrélation, soit 20 fois plus que la variation du gain. C'est ce qui donne à la méthode selon l'invention une si grande sensibilité. Dans un mode de réalisation possible, la durée d'acquisition du spectre de référence est supérieure à celle du premier spectre. On peut par exemple prévoir une durée d'acquisition du spectre de référence de 100s et une durée d'acquisition du premier spectre de 10s. On bénéficie ainsi d'un spectre de référence moins bruité, ce qui permet d'améliorer les résultats de la corrélation. On relèvera que les degrés de corrélation couramment obtenus par les mesures réalisées par les inventeurs sont supérieurs à 0,999. La conséquence d'une telle sensibilité est une stabilité extrêmement élevée sur les taux de comptages obtenus. On a représenté sur la figure 9 la variation de comptage Vco en réponse à une variation brutale de température se traduisant par une variation sur le gain G de 27%. Après l'acquisition d'un spectre de durée 10 secondes et un calcul de quelques secondes, le gain est recalé et le système est prêt pour réaliser un comptage. La figure 9 illustre que la variation brutale de gain n'entraîne pas de variation de comptage Vco perceptible sur les différentes mesures Nm réalisées. On a par ailleurs représenté sur la figure 10 un filtrage de spectre pouvant être mis en oeuvre dans un mode de réalisation possible de l'invention qui comprend avant détermination du degré de corrélation avec le spectre de référence, une étape de filtrage du ou des spectres (spectre de référence, premier spectre, et seconds spectres dans le premier mode de réalisation) acquis par un filtre en rampe. Ce filtrage vise à accroître la robustesse du procédé au bruit de la chaîne de spectrométrie à basse énergie et diminuer l'influence du seuil de coupure dans le calcul du coefficient de corrélation. L'application d'un filtre en rampe (passe haut) se traduit par la multiplication du contenu de chaque canal par son numéro de canal. Ensuite, pour normaliser, l'ensemble est divisé par la moitié du nombre de canaux du spectre. En effet, la partie haute du spectre apporte à l'unité de corrélation plus d'informations que la partie basse. Le bruit en partie basse est sérieusement affaibli. On relèvera que sur la figure 10, la normalisation n'a pas été réalisée telle que précédemment décrite, mais pour mettre en coïncidence le pic présent autour du canal 890 dans le spectre acquis Sa et le spectre acquis filtré Saf. Un seuillage peut par ailleurs être mis en oeuvre pour ignorer le contenu des canaux de plus basses énergies. Une étape préalable de détermination du gain de la chaîne de détection peut par ailleurs être mise en oeuvre pour gagner du temps et diminuer la dynamique du procédé. Ce gain peut tout simplement correspondre au gain recalé lors de la précédente utilisation du système de détection spectroscopique. Il peut également s'agir d'un gain déterminé selon un autre procédé de mesure, tel que par exemple suivant la méthode de suiveur de pics présentée dans le brevet suscité FR 2 626 121 et qui consiste à réaliser une égalisation du comptage des impulsions au sein de deux fenêtres adjacentes dont une première est délimitée entre un premier et un deuxième seuil homothétiques, et la deuxième a pour borne inférieure le deuxième seuil. Une comparaison des performances du procédé selon l'invention avec deux autres méthodes de stabilisation a été réalisée. Ces deux méthodes sont la méthode de suiveur de pic, et la méthode du barycentre du pic d'absorption. Selon cette dernière méthode, on cherche dans le pic d'absorption totale le canal qui a le comptage le plus élevé et nous calculons le barycentre dans une fenêtre qui se situe à ±5 % de la position de ce canal (voir spectre ci-dessous). Le numéro du canal correspondant à la position du EY(c-xx.) barycentre est la suivante : X = ' où N est le nombre de canaux dans la fenêtre Ei' ci du pic, x, le numéro d'un canal et c, le comptage qui lui est associé. La position du barycentre est proportionnelle au gain. Variation du comptage pour une variation de deux fois l'écart-type sur le gain Temps de mesure (s) 0,1 1 10 Nombre de photons comptés 11300 13000 130000 Méthode par barycentre du pic instable 0.04% 0.02% Méthode par seuils suiveurs de pic 2% 0.65% 0.20% Méthode par corrélation entre spectres 0.11% 0.06% 0.03% On relèvera qu'une variation de gain entraîne une variation trois fois plus petite sur le comptage (avec le seuil bas de la fenêtre de comptage placé à 20% de la position du pic). Par ailleurs, la variation du temps d'acquisition permet de faire varier le bruit des comptages du spectre. Ce bruit suit la loi de Poisson et est inversement proportionnel à la racine carrée du temps. La méthode par corrélation est peu sensible au nombre de canaux du spectre, ce qui est un avantage car il est possible d'utiliser un spectromètre peu résolutif. La méthode par corrélation apparaît être la meilleure car elle est très précise et très peu sensible aux bruits (parasites, bruits électroniques et bruits statistiques de comptage).15 REVENDICATIONS1. Procédé de compensation de dérive du gain d'une chaîne de détection d'impulsions issues d'une source, comprenant les étapes suivantes : - acquisition (ACQ) d'un premier spectre d'impulsions issues de la source ; détermination d'un degré de corrélation (CORR) du premier spectre avec un spectre de référence - fixation d'une variation de gain (MOD-GAIN, MOD-EXP) de la chaîne de détection ; détermination d'un degré de corrélation (CORR) d'un second spectre d'impulsions avec le spectre de référence, dans lequel le second spectre correspond directement ou indirectement à une acquisition d'un spectre d'impulsions issues de la source après modification du gain de la chaîne de détection de ladite variation de gain. 2. Procédé selon la revendication 1, dans lequel les étapes de fixation d'une variation de gain (MOD-GAIN, MOD-EXT) et de détermination d'un degré de corrélation (CORR) du second spectre avec le spectre de référence sont réitérées jusqu'à identifier une variation de gain (MAX?) permettant d'atteindre un degré de corrélation donné du second spectre avec le spectre de référence. 3. Procédé selon la revendication 2, dans lequel l'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend la modification du gain de la chaîne de détection (MOD-GAIN) de ladite variation de gain et l'acquisition (ACQ) du second spectre à l'aide de la chaîne de détection à gain modifié. 4. Procédé selon la revendication 2, dans lequel l'étape de détermination d'un degré de corrélation d'un second spectre d'impulsions avec le spectre de référence comprend la détermination d'un coefficient d'expansion (MOD-EXP) à partir de laditevariation de gain et le calcul (CAL) du second spectre par compression ou dilation du premier spectre conformément au coefficient d'expansion. 5. Procédé selon la revendication 4, dans lequel le coefficient d'expansion correspond au ratio d'un nombre de canaux du second spectre et d'un nombre de canaux du premier spectre. 6. Procédé selon la revendication 5, dans lequel le contenu ri, d'un canal d'indice du second spectre correspond à : o = a.ND si A = 0 o ni = ND. (1 - partie décimale [i -1].a) + NF. (partie décimale i. a) si A = 1 o ni = ND. (1 - partie décimale [i a) + EL1 N D + k NF. (partie décimale i. a) si A > 2 avec : i un entier compris entre 1 et la partie entière de N*a, où N désigne le nombre de canaux du premier spectre ; Nj le contenu du canal d'indice] du premier spectre ; D le canal du premier spectre au sein duquel débute le canal d'indice i du second spectre : il s'agit du canal du premier spectre dont l'indice correspond à la partie entière de (i-1)*a; - F le canal du premier spectre au sein duquel s'achève le canal d'indice i du second spectre: il s'agit du canal du premier spectre dont l'indice correspond à la partie entière de i*a ; - A = F-D. 7. Procédé selon l'une des revendications 1 à 6, dans lequel la durée d'acquisition du spectre de référence est supérieure à celle du premier spectre.30 8. Procédé selon l'une des revendications 1 à 7, comprenant, avant détermination du degré de corrélation avec le spectre de référence, une étape de filtrage du premier spectre et du spectre de référence par un filtre en rampe. 9. Procédé selon l'une des revendications 1 à 8, comprenant une étape préalable de détermination du gain de la chaîne de détection par égalisation d'un comptage des impulsions au sein de deux fenêtres adjacentes dont une première est délimitée entre un premier et un deuxième seuils homothétiques, et la deuxième a pour borne inférieure le deuxième seuil. 10. Système de détection spectroscopique d'impulsions issues d'une source, comprenant un module d'étalonnage (1) configuré pour mettre en oeuvre le procédé selon l'une des revendications 1 à 9. 11. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution des étapes du procédé selon l'une des revendications 1 à 9, lorsque ledit programme est exécuté sur un ordinateur.20