Brevet FR3131629A1 — publié le 20230707

L’invention est un procédé de traitement d’un spectre de calibration acquis par un détecteur spectrométrique de photons X ou gamma. Le procédé comporte une prise en compte d’une forme paramétrique d’une fonction de calibration, la fonction de calibration reliant le rang d’un canal d’énergie à une valeur d’énergie. Le procédé comporte une confrontation entre les canaux des pics du spectre de calibration et des énergies d’émission d’isotopes de calibration. La confrontation permet de définir les valeurs des paramètres de la fonction de calibration. Procédé de calibration en énergie d’un détecteur spectrométrique Le domaine technique de l'invention est la spectrométrie appliquée à la détection de rayonnements. ART ANTERIEUR Les dispositifs de détection de rayonnement ionisants, basés sur des matériaux détecteurs gazeux, semi-conducteurs ou scintillateurs, permettent d'obtenir des impulsions électriques formées par des interactions du rayonnement dans le matériau détecteur. L'amplitude de chaque impulsion dépend de l'énergie déposée par le rayonnement au cours de chaque interaction. Ces dispositifs sont usuellement couplés à des circuits de mesure spectrométriques. De tels dispositifs permettent d’obtenir un spectre, qui correspond à un histogramme des amplitudes détectées au cours d’une période d’exposition. Le spectre est établi selon différents canaux d’amplitude, usuellement plusieurs centaines ou plusieurs milliers de canaux. A chaque canal correspond une bande étroite d’amplitude. Chaque canal est assigné d’une valeur, correspondant à un nombre d’impulsions, détectées durant la période d’exposition, dont l’amplitude est située dans la bande d’amplitude correspondant au canal. Les systèmes de mesures spectrométriques sont à présent largement industrialisés. Les domaines d'application sont larges, et comprennent notamment les mesures de déchets, d’équipements ou d’installations nucléaires, ou la surveillance radiologique de l’environnement. Des logiciels permettent le paramétrage du traitement des impulsions, ainsi que le pilotage et le dépouillement automatisé de mesures. Cependant, certaines étapes sont délicates et laissent encore une large place à des réglages manuels. Il s’agit par exemple de la calibration en énergie, qui permet d’établir une relation entre l’amplitude d’une impulsion et l’énergie libérée par l’interaction ayant donné lieu à l’impulsion. La relation entre l’amplitude et l’énergie dépend du détecteur utilisé, ainsi que des paramétrages du circuit de mesure spectrométrique, par exemple le gain d’amplification, des paramètres de mise en forme des impulsions ou le nombre de canaux. De plus, pour un même dispositif, sans modification des réglages, la relation entre l’amplitude et l’énergie peut subir des dérives, ce qui nécessite une calibration périodique. La calibration en énergie est usuellement effectuée en exposant un détecteur à un objet de calibration dont les énergies d’émission sont connues. L’objet de calibration peut notamment comporter un ou plusieurs isotopes dont les énergies d’émission sont connues. La calibration consiste à acquérir un spectre des photons émis par les isotopes, puis d’identifier les principaux pics présents dans le spectre. Une partie des pics détectés correspond à une énergie d’émission connue. Chaque pic détecté s’étend de part et d’autre d’un canal central. La calibration en énergie permet de déterminer différents couples canal (i-e le canal central d’un pic) - énergie d’émission. A partir de ces différents couples, on détermine une fonction analytique, dite fonction de calibration, reliant le rang de chaque canal à une énergie. La calibration en énergie est une opération fréquente, qui doit être effectuée avec précision, afin de ne pas fausser l’interprétation de mesures. L’invention décrite ci-après facilite l’automatisation et la répétabilité de calibration en énergie. Elle peut également permettre de maximiser le nombre de pics pris en compte pour déterminer la fonction de calibration. Un premier objet de l’invention est un procédé de traitement d’un spectre de calibration, formé par un dispositif de mesure spectrométrique, le dispositif comportant : un détecteur, configuré pour détecter des particules, et pour former, à chaque détection, une impulsion dont l’amplitude dépend d’une énergie d’une particule ayant interagi dans le détecteur ; un circuit de mesure spectrométrique, configuré pour former un spectre, le spectre correspondant à un nombre de particules détectées dans des canaux, à chaque canal étant assigné un rang, à chaque rang correspondant une amplitude d’impulsion; le procédé comportant les étapes suivantes : a) disposition du dispositif face à un objet de calibration, l’objet de calibration émettant, vers le détecteur, des particules de différentes énergies d’émission connues ; b) détection d’une partie des particules par le détecteur et formation d’un spectre de calibration des particules détectées, le spectre comportant différents pics, chaque pic correspondant à une énergie connue ; c) assignation d’un rang de canal à chaque pic détecté; d) à partir des rangs des canaux assignés à différents pics lors de l’étape c) et des énergies d’émission, détermination d’une fonction de calibration, la fonction de calibration déterminant une énergie à partir du rang de chaque canal, la fonction de calibration étant une fonction bijective ; le procédé étant caractérisé en ce que l’étape d) comporte les sous-étapes suivantes: d-i) prise en compte d’un modèle analytique de la fonction de calibration, le modèle analytique étant paramétré par un ensemble de paramètres ; d-ii) pour différentes valeurs de chaque paramètre : application de la fonction de calibration au rang de chaque canal assigné à chaque pic détecté, de façon à obtenir, pour chaque rang, une énergie déterminée par la fonction de calibration ; ou application de la fonction réciproque de la fonction de calibration à chaque énergie d’émission, de façon à obtenir, pour chaque énergie d’émission, un canal déterminé par la fonction réciproque ; d-iii) détermination de la valeur de chaque paramètre pour laquelle : les énergies déterminées lors de d-ii) sont les plus proches des énergies d’émission de l’objet de calibration ; les canaux déterminés lors de d-ii) sont les plus proches des canaux assignés à chaque pic détecté. Selon un mode de réalisation, la sous-étape d-iii) peut comporter : calcul d’une fonction de vraisemblance de la valeur de chaque paramètre ; estimation de la valeur de chaque paramètre maximisant la fonction de vraisemblance. L’étape c) peut comporter une détermination d’une largeur de chaque pic. La fonction de vraisemblance peut alors prendre en compte la largeur de chaque pic. La sous-étape d-iii) peut comporter: prise en compte de différents couples, chaque couple comportant une énergie d’émission et un canal assigné à un pic, suite à l’étape c); détermination de la valeur de chaque paramètre pour laquelle la fonction de calibration permet de relier plusieurs couples, chaque valeur déterminée formant une valeur initiale de chaque paramètre ; pour chaque valeur initiale, définition d’un domaine de recherche s’étendant autour de ladite valeur initiale, de façon que la valeur de chaque paramètre maximisant la fonction de vraisemblance est estimée dans les domaines de recherches définis autour des valeurs initiales dudit paramètre. La sous-étape d-iii) est préférentiellement mise en œuvre par un algorithme d’optimisation. Selon une possibilité, l’étape c) comporte une sélection des pics du spectre de calibration en fonction d’un critère de sélection du pic. Le critère de sélection peut être un nombre de photons détectés dans chaque pic ou un rapport signal sur bruit déterminé pour chaque pic. La sous-étape d-iii) peut comporter une prise en compte d’un a-priori sur la valeur d’au moins un paramètre pour déterminer la valeur dudit paramètre. Les particules peuvent être choisies parmi : photons, neutrons, particules chargées. Un deuxième objet de l’invention est un dispositif destiné à acquérir un spectre de particules émises par un objet, le dispositif comportant : un détecteur, configuré pour détecter des particules, et pour former, à chaque détection, une impulsion dont l’amplitude dépend d’une énergie libérée la particule ayant interagi dans le détecteur ; un circuit de mesure spectrométrique, configuré pour former un spectre, le spectre correspondant à une distribution des amplitudes des impulsions détectées par le détecteur ; une unité de traitement, programmée pour mettre en œuvre l’étape d) d’un procédé selon le premier objet de l’invention à partir du spectre. Un troisième objet de l’invention est un support, destiné à être connecté à un ordinateur, comportant des instructions pour la mise en œuvre de l’étape d) d’un procédé selon le premier objet de l’invention à partir d’un spectre représentatif de l’énergie de particules détectées. Le support peut être intégré dans l’ordinateur ou relié à l’ordinateur par une liaison filaire ou sans fil. Un quatrième objet de l’invention est un procédé de traitement d’un spectre de calibration formé par un dispositif de mesure spectrométrique, le dispositif comportant : un détecteur, configuré pour détecter des particules, et pour former, à chaque détection, une impulsion dont l’amplitude dépend d’une masse d’une particule ayant interagi dans le détecteur ; un circuit de mesure spectrométrique, configuré pour former un spectre, le spectre correspondant à un nombre de particules détectées dans des canaux, à chaque canal étant assigné un rang, à chaque rang correspondant une amplitude d’impulsion; le procédé comportant les étapes suivantes : a) exposition du dispositif à un objet de calibration, l’objet de calibration comportant des particules de masses connues ; b) détection de tout ou partie des particules par le détecteur et formation d’un spectre de calibration des particules détectées, le spectre comportant différents pic, chaque pic correspondant à une masse connue ; c) assignation d’un rang de canal à chaque pic ; d) à partir des rangs des canaux assignés à différents pics lors de l’étape c) et des différentes masses connues, détermination d’une fonction de calibration, la fonction de calibration déterminant une masse à partir du rang de chaque canal, la fonction de calibration étant une fonction bijective ; le procédé étant caractérisé en ce que l’étape d) comporte : d-i) prise en compte d’un modèle analytique de la fonction de calibration, le modèle analytique étant paramétré par un ensemble de paramètres ; d-ii) pour différentes valeurs de chaque paramètre : application de la fonction de calibration au rang de chaque canal assigné à chaque pic, de façon à obtenir, pour chaque rang, une masse déterminée par la fonction de calibration ; ou application de la fonction réciproque de la fonction de calibration à chaque masse, de façon à obtenir, pour chaque masse, un canal déterminé par la fonction réciproque ; d-iii) détermination des valeurs des paramètres pour lesquels : les énergies déterminées lors de d-ii) sont les plus proches des masses des particules de l’objet de calibration ; les canaux déterminés lors de d-ii) sont les plus proches des canaux assignés à chaque pic. L'invention sera mieux comprise à la lecture de l'exposé des exemples de réalisation présentés, dans la suite de la description, en lien avec les figures listées ci-dessous. FIGURES La représente un exemple de dispositif permettant une mise en œuvre de l’invention. La représente un exemple de spectre. La est un détail de la . La montre un pic d’un spectre ainsi que l’application d’un opérateur différentiel linéaire d’ordre 2 au spectre. La montre les principales étapes d’un mode de réalisation d’un procédé selon l’invention. La illustre des correspondances entre des énergies d’émission et des canaux centraux de pics de détection. La illustre des correspondances entre des canaux centraux de pics de détection (axe des abscisses) et des énergies d’émission (axe des ordonnées). Les figures 7A et 7B montrent un spectre de calibration acquis avec un détecteur Germanium. On a également représenté la position de canaux assignés à différents pics détectés. La montre la valeur d’une fonction de vraisemblance d’un paramètre de la fonction de calibration pour différentes valeurs dudit paramètre. La correspond au spectre représenté sur les figures 7A et 7B. La est un détail de la . La montre la valeur d’une fonction de vraisemblance d’un autre paramètre de la fonction de calibration pour différentes valeurs dudit paramètre. La correspond au spectre représenté sur les figures 7A et 7B. La est un détail de la . La illustre des correspondances entre des canaux assignés à des pics de détection (axe des abscisses) du spectre représenté sur la et des énergies d’émission (axe des ordonnées). La montre un spectre de calibration acquis avec un détecteur LaBr3 (Bromure de Lanthane). On a également représenté la position de canaux assignés à différents pics détectés ainsi qu’une plage d’incertitude relative à la position. La montre la valeur d’une fonction de vraisemblance d’un paramètre de la fonction de calibration pour différentes valeurs dudit paramètre. La correspond au spectre représenté sur la . La est un détail de la . La montre la valeur d’une fonction de vraisemblance d’un autre paramètre de la fonction de calibration pour différentes valeurs dudit paramètre. La correspond au spectre représenté sur la . La est un détail de la . La illustre des correspondances entre des canaux assignés à des pics de détection (axe des abscisses) du spectre représenté sur la et des énergies d’émission (axe des ordonnées). EXPOSE DE MODES DE REALISATION PARTICULIERS On a représenté, sur la , un dispositif 1 permettant la mise en œuvre d’invention. Le dispositif est un système de mesure, comportant un détecteur 10, apte à interagir avec un rayonnement 5 émis par un objet 2. L’objet 2 est ici un déchet nucléaire, pouvant comporter différents isotopes irradiants 2 j . Généralement, les isotopes irradiants susceptibles d’être présents dans un objet mesuré sont préalablement connus. On peut notamment constituer une liste comportant des isotopes irradiants potentiellement présents dans l’objet analysé. L’indice j désigne chaque isotope irradiant. Dans l’exemple représenté, le détecteur 10 est configuré pour détecter un rayonnement photonique ionisant. Par rayonnement photonique ionisant, on entend un rayonnement photonique de type X ou gamma, formé de photons dont l'énergie est par exemple comprise entre 1 keV et 2 MeV. Le procédé selon l’invention s’applique à la détection d’autres particules, par exemple des neutrons, des particules chargées (par exemple rayonnement α ou rayonnement β-), ou des photons. Dans l'exemple représenté, le détecteur comporte un matériau semi-conducteur, de type Germanium (Ge), mais il pourrait également s'agir d'un matériau semi-conducteur communément mis en œuvre pour la détection de photons ionisants, par exemple de type Si, CdTe, CdZnTe. Les photons formant le rayonnement incident interagissent dans le matériau détecteur. Le matériau détecteur est soumis à une tension de polarisation V . Chaque interaction génère des porteurs de charges, lesquels sont collectés par une électrode, généralement une anode. D'autres types de détecteurs, par exemple des scintillateurs couplés à un convertisseur photons/porteurs de charge, ou un détecteur gazeux de type chambre d'ionisation, sont utilisables, dès lors qu'ils permettent la collecte d'une quantité de charges sous l'effet d'une énergie libérée par le rayonnement ionisant au cours d'une interaction dans le détecteur 10. Parmi les détecteurs usuels de type scintillateur, on peut citer NaI(Tl) ou LaBr 3 . Le détecteur 10 est relié à un circuit électronique 12, configuré pour générer une impulsion dont l’amplitude dépend, et est de préférence proportionnelle, à la quantité de charge collectée lors d'une interaction. La quantité de charge correspond à l’énergie déposée par le rayonnement au cours de l'interaction. Dans l’exemple représenté, le détecteur 10 est relié à un cryostat 16, comportant de l’azote liquide pour maintenir le détecteur Ge à une température de fonctionnement. Le circuit électronique 12 est relié à une unité de spectrométrie 13, disposée en aval du circuit électronique, qui permet de rassembler l’ensemble des impulsions formées durant une durée d’acquisition. Chaque impulsion correspond à une interaction du rayonnement incident dans le matériau de détection. Le circuit de spectrométrie classe alors les impulsions en fonction de leur amplitude A , pour fournir un histogramme comportant le nombre d’impulsions détectées en fonction de leur amplitude. Cet histogramme est un spectre en amplitude. Il est usuellement obtenu en utilisant un analyseur multicanaux. Chaque amplitude est discrétisée selon des canaux, à chaque canal étant assignée une bande d’amplitude. La valeur de chaque canal du spectre correspond à un nombre d’impulsions dont l’amplitude est située dans la bande d’amplitude assignée au canal. Chaque bande d’amplitude correspond à une bande d’énergie, la correspondance étant bijective. Ainsi, à chaque canal est assignée une bande d’énergie ou une bande d’amplitude. La relation entre amplitude et énergie peut être réalisée en irradiant le détecteur à l’aide d’un objet de calibration, émettant un rayonnement dont on connaît l’énergie. Il s'agit en particulier d'un rayonnement présentant au moins une discontinuité, ou pic d'énergie, à une valeur connue de l'énergie. Cette opération est usuellement désignée par le terme étalonnage en énergie. Dans le cadre de la spectrométrie gamma, on expose le détecteur à une source de calibration de type 152 Eu, produisant des photons à des énergies d’émission connues. On peut également mettre en œuvre d’autres sources isotopiques, comme décrit par la suite dans un des exemples de réalisation. Il est préférable que l’isotope ou les isotopes soient tels que des photons soient émis selon au moins deux énergies d’émission différentes. Les énergies d’émission s’étendent entre une valeur minimale et une valeur maximale. Il est préférable que l’écart entre la valeur minimale et la valeur maximale couvre le domaine spectral dans lequel le détecteur est destiné à être utilisé. La calibration en énergie permet d’établir une fonction de calibration , permettant d’établir une relation analytique entre l’amplitude mesurée et l’énergie. La prise en compte de la fonction de calibration permet, par un changement de variable, d’assigner la valeur des canaux à des valeurs d’énergies à la place des amplitudes. En effet, les amplitudes dépendent du détecteur et des réglages effectués, tandis que les énergies correspondent à des grandeurs physiques fixes, indépendantes du détecteur. Le spectre y est un histogramme des amplitudes de chaque impulsion détectée, discrétisée selon des canaux d’énergie ou d’amplitude. A chaque canal est assignée une bande d’amplitude qui, par l’application de la fonction de calibration, devient une bande d’énergie. Le spectre peut être exprimé sous la forme d’un vecteur ( ), où correspond au nombre total de canaux. A chaque canal est assigné un rang , avec . Chaque canal de rang est délimité par une amplitude inférieure et une amplitude supérieure , de telle sorte qu’une impulsion détectée est affectée au canal de rang lorsque son amplitude est comprise entre et . L’amplitude inférieure de chaque canal correspond à une énergie inférieure L’amplitude supérieure de chaque canal correspond à une énergie supérieure La correspondance amplitude-énergie est établie par la fonction de calibration décrite en lien avec l’art antérieur. Ainsi, et est un ensemble de paramètres de la fonction , décrit par la suite. La grandeur correspond à une position du centre d’un canal de rang relativement au nombre de canaux du spectre. Il s’agit d’un rang normalisé de chaque canal, compris entre 1 (k = 1) et (k = ). La normalisation permet d’établir une fonction de calibration indépendante du nombre de canaux , ce dernier étant paramétrable. Le dispositif comporte une unité de traitement 14, programmée pour mettre en œuvre des étapes de l’algorithme décrit en lien avec la . L’unité de traitement 14 est reliée à l’unité de spectrométrie 13, dont elle reçoit chaque spectre mesuré. L’objectif d’une mesure par spectrométrie gamma est d’identifier les isotopes 2 j présents dans l’objet 2 et, de préférence, d’estimer leurs activités respectives. La figure 2A représente un spectre résultant d’un détecteur Germanium, exposé à une source comportant l’isotope 152 Eu. L’axe des abscisses correspond aux canaux et l’axe des ordonnées correspond au nombre d’impulsions dont l’amplitude correspond à chaque canal. La est un détail de la , correspondant à la zone rectangulaire dessinée sur la . Le spectre représenté comporte plusieurs pics, chaque pic correspondant à une énergie d’émission d’un isotope présent dans l’objet 2. Ces pics forment l’information utile de chaque spectre, à partir de laquelle on peut identifier les isotopes et quantifier leurs activités respectives. Le spectre comporte également un continuum, correspondant à des diffusions de photons dans le détecteur ou avant d’atteindre le détecteur. Le continuum correspond à la partie du spectre sous et entre chaque pic. Le spectre comporte également une composante de bruit de fond, se traduisant par des fluctuations statistiques. Une des étapes courantes d’un traitement de spectre est la détection de chaque pic. Sur la , on a représenté différents pics. L’axe des abscisses correspond aux canaux d’amplitude tandis que l’axe des ordonnées correspond au nombre d’impulsion détectées dans chaque canal. L’unité de traitement peut ensuite déterminer un canal caractéristique de chaque pic, par exemple le canal central autour duquel le pic est centré. Le canal central est usuellement le canal comportant le plus grand nombre d’impulsions détectées. Chaque pic est délimité par deux canaux, constituant respectivement les bornes inférieures et supérieures de chaque pic. Par la suite, chaque pic est assigné à un ordre , par ordre d’amplitude croissante. Chaque pic d’ordre est identifié par le canal de rang de son canal central, ce dernier étant disposé entre la borne inférieure et la borne supérieure du pic. Selon une possibilité, on peut délimiter manuellement chaque pic. La détection et la délimitation de chaque pic peuvent également être effectuées de façon automatique, en mettant par exemple en œuvre un opérateur de différenciation d’ordre 2, noté . La figure 3 représente une partie d’un spectre correspondant à un pic. Sur la figure 3, on a également représenté un opérateur de différenciation d’ordre 2 appliqué au spectre , noté . La ligne verticale, qui correspond à un minimum local de , permet d’identifier le canal de rang correspondant au centre du pic. L’opérateur de différenciation peut être un filtre de Savitzky et Golay d’ordre 3 et de largeur de fenêtre 11 La schématise les principales étapes d’un procédé de détermination d’une fonction de calibration, ces étapes étant décrites ci-après. Etape 100 : acquisition du spectre de calibration . Au cours de cette étape, le détecteur est exposé à un objet de calibration. L’objet de calibration comporte généralement un ou plusieurs isotopes de calibration dont les énergies d’émission sont connues. Etape 110 : détection des pics du spectre de calibration et assignation d’un canal à chaque pic détecté. On note les rangs des canaux assignés à chaque pic détecté sur le spectre de calibration. désigne le nombre de pics détectés sur le spectre de calibration , qui correspondent respectivement à des énergies d’émission des isotopes de calibration. Chaque pic du spectre de calibration s’étend selon une largeur à mi-hauteur . Un exemple de largeur à mi-hauteur est représenté sur la . La largeur à mi-hauteur est également déterminée lors de la détection des pics par des moyens usuels. Par exemple, lorsque la détection des pics est effectuée en appliquant un filtre de Savitsky et Golay tel que décrit en lien avec la figure 3, la largeur à mi-hauteur de chaque pic peut être estimée en déterminant le nombre de canaux entre deux passages par zéro de la fonction de part et d’autre du minimum local correspondant au pic. Sur la figure 3, on a représenté une largeur à mi-hauteur . La largeur à mi-hauteur peut être déterminée par d’autres moyens usuels de logiciels de spectrométrie, suite à une délimitation manuelle des pics. Etape 120 : prise en compte d’un modèle analytique de la fonction de calibration . Au cours de cette étape, on définit un modèle analytique de la fonction de calibration, par exemple une fonction linéaire ou un polynôme de degré prédéterminé. Le modèle est paramétré par un ensemble de paramètres . On suppose que chaque énergie correspondant à un pic détecté est située au centre du canal de rang qui rassemble les impulsions d’énergies comprises entre et . Ainsi, on suppose que (3). Compte tenu de (1), (2), et (3) . où : correspond au nombre total de canaux du spectre. et sont les amplitudes respectivement minimales et maximales du canal d’amplitude de rang correspondant à l’énergie . correspond à une amplitude médiane du canal de rang , qui est supposée correspondre à l’énergie d’émission . La fonction de calibration peut être linéaire, auquel cas l’ensemble de paramètres est ; et : Dans le cas d’un polynôme de degré 2, la fonction de calibration est de type : . L’ensemble des paramètres est ; . La fonction de calibration est nécessairement bijective. Il est usuel que les canaux soient classés selon des amplitudes croissantes, auquel cas la fonction de calibration est croissante. représente l’énergie correspondant à l’amplitude médiane du canal de rang , par l’application de la fonction de calibration . Ainsi, la fonction de calibration fait correspondre, à chaque canal de rang , une énergie . L’objet de calibration émet des photons à des énergies d’émission . correspond au nombre d’énergies d’émission. Tous les pics, correspondant à des énergies d’émission, peuvent ne pas être détectés par la détection de pics, par exemple si la quantité de photons dans un pic est trop faible. Les rangs correspondant à certains pics détectés correspondent à certaines énergies d’émission . La figure 5 représente un exemple de canaux ( = 5)correspondant à chaque pic détecté pris dans l’ordre croissant ( ), ainsi que les énergies d’émission des isotopes, également pris dans l’ordre ( croissant ( = 7). La fonction de calibration doit maximiser le nombre d’occurrences selon lesquelles . Sur la figure 5, chaque occurrence est représentée par une flèche. Lorsque deux énergies d’émission sont rapprochées, elles peuvent être associées à un canal de même rang. Sur la figure 5, ce cas particulier est illustré par les énergies et . Etape 130 : Etablissement d’une densité de probabilité A chaque énergie d’émission correspond un canal de rang par l’application de la fonction , de telle sorte que La notation désigne le fait que dépend de . Comme précédemment indiqué, un pic détecté présente une certaine largeur spectrale , induisant une incertitude sur le rang du canal . On peut considérer que pour chaque énergie d’émission , le rang du canal dans lequel le pic est détecté est une variable suivant une densité de probabilité prédéterminée , par exemple gaussienne, telle que : où : est le rang « vrai » du canal associé au pic d’énergie d’émission ; est le rang assigné au ème pic détecté ; . Partant de (7), on peut écrire : où : est un vecteur comportant l’ensemble des rangs des canaux correspondant à des pics détectés : est une probabilité qu’un pic de rang corresponde à une énergie d’émission ; est le vecteur ( ). Toutes les combinaisons sont considérées équiprobables. Ainsi, Etape 140 : établissement d’une fonction de vraisemblance et maximisation. En prenant en compte (7), (8), et (9), on peut écrire : On en déduit une fonction de vraisemblance de , notée telle que : correspond aux énergies d’émission des radionucléides utilisés pour la calibration : . On peut estimer le vecteur de paramètres qui maximise la fonction de vraisemblance : Le vecteur est celui pour lequel le nombre de correspondances entre et est maximisé. Cela revient à obtenir un maximum d’occurrences selon lesquelles . Le vecteur ainsi obtenu permet de définir la fonction de calibration . La figure 6 illustre une situation selon laquelle la fonction est linéaire. La fonction est telle que : pour chaque canal de rang . De façon réciproque, une énergie correspond à un canal de rang tel que : L’expression (13) devient : L’algorithme d’optimisation, exprimé par (16), peut avantageusement être mis en œuvre autour de valeurs de et prédéterminées. Cela permet d’applique l’algorithme de maximisation dans un domaine de recherche prédéterminé autour des valeurs initiales prédéterminées de et Cela permet d’éviter des « pièges d’optimisation », correspondant à la détermination de maximas locaux de et . Cela permet également d’optimiser le temps de calcul. Le domaine de recherche est défini autour de chaque couple de valeurs et . D’une façon plus générale, quelle que soit le type de fonction de calibration, le vecteur est avantageusement déterminé en appliquant un algorithme d’optimisation, permettant de résoudre (16), sur la base d’un domaine de recherche restreint autour de valeurs initiales, discrètes, du vecteur . Le domaine de recherche correspond à un intervalle continu s’étendant autour de chaque valeur initiale. La définition du domaine de recherche restreint est décrite ci-dessous, en lien avec la . Sur la figure 6, l’axe des abscisses correspond au rang normalisé de chaque canal. L’axe des ordonnées correspond aux énergies d’émission normalisées par l’énergie d’émission la plus élevée avec . Une énergie d’émission est détectée lorsqu’un canal lui est assigné. On a représenté, par des lignes verticales noires, les canaux ayant été identifiés par la recherche de pics sur le spectre de calibration. Comme précédemment indiqué, les canaux d’énergie sont affectés d’une incertitude , qui est matérialisée par deux lignes verticales fines de part et d’autre d’une ligne verticale noire. Les énergies d’émission normalisées correspondent à des lignes horizontales. Il s’agit de valeurs exactes, sans incertitude associée. On a tracé, en gris, des droites qui passent par au moins deux intersections entre une énergie d’émission normalisée et un canal de détection , en prenant en compte l’incertitude associée. Chaque intersection correspond à une correspondance , . On a pris en compte le fait que la fonction de calibration est nécessairement croissante lorsque les rangs des canaux sont ordonnés de façon croissante avec l’amplitude. Ainsi, on a représenté uniquement des droites pouvant correspondre à une fonction croissante ( . La droite représentée en noir est la droite présentant le nombre maximal d’intersections, en l’occurrence ( , ), ( , ), ( , ). La figure 6 illustre la définition d’un domaine de recherche restreint pour mettre en œuvre l’algorithme d’optimisation précédemment décrit. Il s’agit de déterminer des valeurs initiales des paramètres de la fonction de calibration pour lesquelles l’application de la fonction de calibration, bijective, génère des occurrences entre (ou et ou, de façon réciproque, entre et . Autrement dit, on forme différents couples, chaque couple comportant une énergie d’émission et le rang d’un canal assigné à un pic détecté. L’algorithme d’optimisation décrit dans (16) est mis en œuvre autour de valeurs initiales du vecteur de paramètres de la fonction de calibration pour lesquelles la fonction de calibration est la fonction vérifiant pour au moins deux couples dont les énergies et les rangs respectifs sont différents d’un couple à l’autre : autrement dit, pour deux couples et , et A partir de la figure 6, on peut établir différentes valeurs initiales de et avec : (17) et (18) Selon une variante, les étapes 130 et 140 sont mises en œuvre en sélectionnant les énergies d’émission les plus intenses et/ou les canaux pour lesquels les surfaces des pics sont les plus élevées, ou présentent un rapport signal sur bruit considéré comme élevé. Cela permet d’optimiser le temps de calcul. Selon cette variante, le nombre d’énergies d’émission et de canaux peuvent être limitées à quelques dizaines. On peut prendre par exemple = 30 et = 20. D’une façon plus générale, le procédé comporte une prise en compte d’un critère de sélection de chaque pic. Les pics traités dans les étapes 130 et 140 sont ceux qui satisfont au critère de sélection. Selon une variante, lors de l’étape 140, on prend en compte un a priori sur les valeurs d’un ou plusieurs paramètres formant l’ensemble . Par exemple, on définit une plage de variation pour un ou plusieurs paramètres. On impose que chaque paramètre, pour lequel un a priori est pris en compte, soit situé dans la plage de variation prédéfinie. La prise en compte d’un a priori permet d’éviter que l’algorithme d’optimisation génère des valeurs non réalistes de paramètres. L’a-priori constitue une contrainte pour la mise en œuvre de l’algorithme d’optimisation, s’ajoutant à la condition selon laquelle la fonction de calibration est monotone. Essais expérimentaux Les inventeurs ont mis en œuvre le procédé de calibration tel que précédemment décrit en utilisant deux types de détecteurs : un détecteur de type Ge (Germanium), tel que schématisé sur la , et un détecteur de type LaBr 3 (Bromure de Lanthane). Les résultats des calibrations sont à présent reportés. Au cours d’une première série d’essais, on a effectué une calibration en énergie du détecteur Germanium. On a utilisé une source étalon de 152 Eu, d’activité 35.04 KBq. On a acquis un spectre de calibration durant une durée de 3 heures. La source était placée à une distance de 40 cm du détecteur. Le spectre a été acquis selon 8192 canaux. La fonction de calibration a été considérée linéaire. Les figures 7A et 7B représentent les canaux résultant de l’algorithme de détection de pic. On observe qu’il n’y a pas de pic détecté au-delà du canal de rang 6600. Sur les figures 7A et 7B, l’axe des abscisses correspond aux canaux et l’axe des ordonnées aux nombres d’impulsions (ou coups) détectées et assignées à chaque canal. Il s’agit d’une représentation spectrale classique. La détection des pics a permis d’obtenir 179700 couples de correspondances possibles canal d’un pic détecté – énergie d’émission. En prenant en compte le fait que la fonction de calibration est croissante, on obtient 82650 couples possibles. Cela correspond à = 20 (nombre d’énergies d’émission considérées) et = 30 (nombre de pics détectés), le nombre de couples possibles, correspondant à une fonction de calibration croissante, est égal à : . (19) La figure 7C montre la valeur de la fonction de vraisemblance L (axe des ordonnées) en fonction de (axe des abscisses). La figure 7D est un détail de la figure 7B. La valeur de maximisant la fonction de vraisemblance est égale à 0.3689 keV. La figure 7E montre la valeur de la fonction de vraisemblance L (axe des ordonnées) en fonction de (axe des abscisses). La figure 7F est un détail de la figure 7E. La valeur de maximisant la fonction de vraisemblance est égale à 1907.2475 keV. La figure 7G est une représentation similaire à la figure 6 précédemment décrite. L’axe des abscisses correspond aux rangs des canaux des pics détectés, ainsi que l’incertitude associée. L’axe des ordonnées correspond aux énergies d’émission de 152 Eu (keV). L’application de la fonction de calibration prenant en compte les valeurs et optimales permet de maximiser le nombre de correspondances entre les canaux des pics détectés et les énergies d’émission. Au cours d’une deuxième série d’essais, on a effectué une calibration en énergie du détecteur LaBr 3 . On a utilisé des sources étalon de 22 Na, 57 Co, 60 Co et 133 Ba. Les énergies d’émission prépondérantes étaient 276.4 keV, 302.85 keV, 356.01 keV, 383.85 keV, 511 keV, 1173.23 keV, 1274.54 keV, 1332.49 keV. Le spectre a été acquis selon 4096 canaux. La fonction de calibration a été considérée linéaire. La représente les canaux résultant de l’algorithme de détection de pic. 8 pics ont été détectés, ce qui correspond au nombre d’énergie d’émission. La détection des pics a permis d’obtenir 1568 points en prenant en compte le fait que la fonction de calibration est croissante. La figure 8B montre la valeur de la fonction de vraisemblance L (axe des ordonnées) en fonction de (axe des abscisses). La figure 8C est un détail de la figure 8B. La valeur de maximisant la fonction de vraisemblance est égale à -120.7273 keV. La figure 8D montre la valeur de la fonction de vraisemblance L (axe des ordonnées) en fonction de (axe des abscisses). La figure 8E est un détail de la figure 8D. La valeur de maximisant la fonction de vraisemblance est égale à 1607.1911 keV. La figure 8F est similaire aux figures 7F et 6 précédemment décrites. L’axe des abscisses correspond aux rangs des canaux des pics détectés, ainsi que l’incertitude associée. L’axe des ordonnées correspond aux énergies d’émission précédemment listées. L’application de la fonction de calibration prenant en compte les valeurs et permet de maximiser le nombre de correspondances entre les canaux des pics détectés et les énergies d’émission. La fonction de calibration peut être différente d’une fonction linéaire. ll peut par exemple s’agir d’un polynôme de degré supérieur ou égal à 2. L’invention a été décrite en lien avec la spectrométrie gamma. Elle peut être généralisée à la spectrométrie de rayonnement d’autres types de photons, ou d’autres types de particules, en particulier neutrons ou particules chargées. Elle peut également être mise en œuvre dans la détection de photons non ionisants, par exemple des photons des domaines spectraux infra-rouge, visible ou proche ultraviolet. On comprend que l’invention peut être appliquée à la spectrométrie de masse. Dans ce cas, le détecteur spectrométrique est configuré pour former des impulsions dont l’amplitude dépend de la masse de particules ayant interagi dans le détecteur. L’établissement de la fonction de calibration est effectué en exposant le détecteur à un objet de calibration, comportant des particules de différentes masses, ces dernières étant connues. La fonction de calibration permet d’effectuer un lien bijectif entre l’amplitude de chaque impulsion et la masse des particules détectées par le détecteur. Procédé de traitement d’un spectre de calibration ( ) formé par un dispositif de mesure spectrométrique (1), le dispositif comportant : un détecteur (10), configuré pour détecter des particules, et pour former, à chaque détection, une impulsion dont l’amplitude dépend d’une énergie d’une particule ayant interagi dans le détecteur ; un circuit de mesure spectrométrique (12, 13), configuré pour former un spectre, le spectre correspondant à un nombre de particules détectées dans des canaux, à chaque canal étant assigné un rang, à chaque rang correspondant une amplitude d’impulsion; le procédé comportant les étapes suivantes : a) disposition du dispositif face à un objet de calibration (2), l’objet de calibration émettant, vers le détecteur, des particules de différentes énergies d’émission ( )connues ; b) détection d’une partie des particules par le détecteur et formation d’un spectre de calibration ( ) des particules détectées, le spectre comportant différents pics, chaque pic correspondant à une énergie connue ; c) assignation d’un rang de canal ( ) à chaque pic détecté; d) à partir des rangs des canaux assignés à différents pics lors de l’étape c) et des énergies d’émission ( ), détermination d’une fonction de calibration, la fonction de calibration déterminant une énergie à partir du rang de chaque canal, la fonction de calibration étant une fonction bijective ; le procédé étant caractérisé en ce que l’étape d) comporte les sous-étapes suivantes : d-i) prise en compte d’un modèle analytique de la fonction de calibration ( ), le modèle analytique étant paramétré par un ensemble de paramètres ( ); d-ii) pour différentes valeurs de chaque paramètre : application de la fonction de calibration au rang de chaque canal assigné à chaque pic détecté, de façon à obtenir, pour chaque rang, une énergie déterminée par la fonction de calibration ; ou application de la fonction réciproque de la fonction de calibration à chaque énergie d’émission, de façon à obtenir, pour chaque énergie d’émission, un canal déterminé par la fonction réciproque ; d-iii) détermination de la valeur de chaque paramètre pour laquelle : les énergies déterminées lors de d-ii) sont les plus proches des énergies d’émission de l’objet de calibration ; les canaux déterminés lors de d-ii) sont les plus proches des canaux assignés à chaque pic détecté. Procédé selon la revendication 1, dans lequel la sous-étape d-iii) comporte : calcul d’une fonction de vraisemblance de la valeur de chaque paramètre ; estimation de la valeur de chaque paramètre ( ) maximisant la fonction de vraisemblance. Procédé selon la revendication 2, dans lequel l’étape c) comporte une détermination d’une largeur ( ) de chaque pic ; la fonction de vraisemblance prend en compte la largeur de chaque pic. Procédé selon l’une quelconque des revendications 2 ou 3, dans lequel la sous-étape d-iii) comporte: prise en compte de différents couples ( ), chaque couple comportant une énergie d’émission et un canal assigné à un pic, suite à l’étape c); détermination de la valeur de chaque paramètre pour laquelle la fonction de calibration permet de relier plusieurs couples, chaque valeur déterminée formant une valeur initiale de chaque paramètre ; pour chaque valeur initiale, définition d’un domaine de recherche s’étendant autour de ladite valeur initiale, de façon que la valeur de chaque paramètre maximisant la fonction de vraisemblance est estimée dans les domaines de recherches définis autour des valeurs initiales dudit paramètre. Procédé selon l’une quelconque des revendications 2 à 4, dans lequel la sous-étape d-iii) est mise en œuvre par un algorithme d’optimisation. Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape c) comporte une sélection des pics du spectre de calibration en fonction d’un critère de sélection du pic. Procédé selon la revendication 6, dans lequel le critère de sélection est un nombre de photons détectés dans chaque pic ou un rapport signal sur bruit déterminé pour chaque pic. Procédé selon l’une quelconque des revendications précédentes, dans lequel la sous-étape d-iii) comporte une prise en compte d’un a-priori sur la valeur d’au moins un paramètre pour déterminer la valeur dudit paramètre. Procédé selon l’une quelconque des revendications précédentes, dans lequel les particules sont choisies parmi : photons, neutrons, particules chargées. Dispositif (1) destiné à acquérir un spectre de particules émises par un objet, le dispositif comportant : un détecteur (10), configuré pour détecter des particules, et pour former, à chaque détection, une impulsion dont l’amplitude dépend d’une énergie libérée la particule ayant interagi dans le détecteur ; un circuit de mesure spectrométrique (12), configuré pour former un spectre ( ), le spectre correspondant à une distribution des amplitudes des impulsions détectées par le détecteur ; une unité de traitement (14), programmée pour mettre en œuvre l’étape d) d’un procédé selon l’une quelconque des revendications précédentes à partir du spectre. Support, destiné à être connecté à un ordinateur, comportant des instructions pour la mise en œuvre de l’étape d) d’un procédé selon l’une quelconque des revendications 1 à 9 à partir d’un spectre représentatif de l’énergie de particules détectées.