Brevet FR3061554A1 — publié le 20180706

Procédé de traitement d'un spectre d'un rayonnement transmis par un objet Description DOMAINE TECHNIQUE Le domaine technique de l'invention est l'analyse d'un objet par une irradiation de ce dernier par un rayonnement électromagnétique ionisant, et d'une acquisition d'un spectre d'un rayonnement transmis par l'objet. ART ANTERIEUR L'analyse d'objets par irradiation à l'aide de rayonnements électromagnétiques ionisants, par exemple un rayonnement X ou gamma permet d'estimer la nature des matériaux composant cet objet. Cette analyse peut être bidimensionnelle, de type radiographie. Il peut également s'agir d'une analyse tridimensionnelle, de type tomographie. Une telle analyse peut être réalisée à des fins de diagnostic médical, de contrôle non destructif dans le domaine industriel, ou encore dans la détection de matériaux sensibles, par exemple dans le contrôle de bagages. Récemment, l'émergence de détecteurs spectralement résolus a permis d'améliorer sensiblement les performances de ces analyses. De tels détecteurs permettent d'acquérir un spectre d'un rayonnement transmis par un matériau irradié. Ce spectre comporte plusieurs canaux, chaque canal étant représentatif d'une bande d'énergie. Un tel spectre peut généralement comporter des dizaines voire des centaines de canaux. Généralement, plusieurs spectres d'un même objet sont acquis, de façon à reconstruire l'objet, c'est-à-dire obtenir une distribution spatiale de caractéristiques de l'objet. Il peut s'agir d'une estimation de la distribution spatiale de la densité et ou du numéro atomique de l'objet, ou d'une estimation de la nature de matériaux présents dans l'objet. Une pratique courante est alors de regrouper les canaux, de façon à concentrer l'information sur un nombre raisonnable de canaux, généralement compris entre 2 et 20, et typiquement 5. Un tel regroupement permet d'optimiser la reconstruction, en particulier le rapport signal sur bruit. Le regroupement est fréquemment effectué par des circuits électroniques, en sortie de détecteurs. En suivant des protocoles de traitement connus, il est alors possible de décomposer l'objet en différentes mailles élémentaires et d'estimer la nature des matériaux composant chaque maille. 3061554 2 L'invention est un perfectionnement d'un tel procédé d'analyse, à partir d'un détecteur spectrométrique, permettant d'obtenir une estimation plus précise des matériaux constituant un objet. L'invention peut être mise en oeuvre préalablement ou durant une reconstruction tomographique, auquel cas la qualité de la reconstruction est améliorée. 5 EXPOSE DE L'INVENTION Un objet de l'invention est un procédé d'analyse d'un objet comportant les étapes suivantes : a) irradiation de l'objet à l'aide d'une source de rayonnement, ladite source émettant un rayonnement électromagnétique ionisant se propageant vers l'objet ; b) obtention d'un spectre d'un rayonnement transmis par l'objet à l'aide d'un détecteur 10 de rayonnement, appartenant à un dispositif de détection; le spectre ainsi mesuré comportant plusieurs de canaux élémentaires, chaque canal élémentaire étant associée à une bande d'énergie ; c) à partir du spectre obtenu, formation d'un spectre de regroupement constitué de plusieurs canaux, dits canaux de regroupement, au moins un canal de regroupement 15 étant obtenu en regroupant, selon un critère de regroupement prédéterminé, plusieurs canaux élémentaires, la valeur du spectre de regroupement en chaque canal de regroupement étant fonction de la valeur du spectre, obtenu lors de l'étape b), mesuré en chaque canal élémentaire ainsi regroupé ; d) sélection d'une base comportant des vecteurs de base, dont chaque terme représente 20 une atténuation du rayonnement électromagnétique ionisant dans une plage d'énergie associée à chaque canal de regroupement ; e) obtention de coordonnées d'une atténuation correspondant au spectre de regroupement, formé lors de l'étape c), dans la base sélectionnée lors de l'étape d) ; f) estimation d'une caractéristique de l'objet en fonction des coordonnées obtenues lors 25 de l'étape e). L'étape c) peut comporter, pour chaque canal de regroupement, les sous étapes suivantes : - en fonction du critère de regroupement, association d'une plage de regroupement, définissant les canaux élémentaires regroupés dans le canal de regroupement ; prise en compte d'une fonction de regroupement, de telle sorte que la valeur du spectre de 30 regroupement au canal de regroupement soit déterminé en appliquant la fonction de regroupement à la valeur du spectre, obtenu lors de l'étape b), dans chaque canal élémentaire ainsi regroupé. 3061554 3 Lors de l'étape c), la valeur attribuée à chaque canal de regroupement peut être une somme ou une somme pondérée de la valeur du spectre, obtenu lors de l'étape b), en chaque canal élémentaire regroupé. Lorsqu'il s'agit d'une somme pondérée, le contenu de chaque canal élémentaire est affecté d'un coefficient de pondération dépendant de la bande d'énergie 5 associée audit canal élémentaire. Lors de l'étape c), le critère de regroupement peut être un critère d'uniformité, selon lequel les valeurs respectives de chaque canal de regroupement composant le spectre de regroupement sont égales, à une tolérance près. Lors de l'étape d), chaque terme des vecteurs de base est défini pour chaque canal de 10 regroupement et correspond, dans la plage d'énergie associée audit canal de regroupement, - soit à une atténuation d'un matériau de calibration ; - soit à une proportion relative d'atténuation photoélectrique ou de diffusion Compton dans ladite plage d'énergie; - soit à une discontinuité spectrale. 15 Le dispositif de détection peut comporter une pluralité de détecteurs de rayonnements adjacents les uns des autres, les étapes b) à e) étant mises en oeuvre pour chaque détecteur de rayonnement, de telle sorte que : chaque étape c) est mise en oeuvre en considérant un même critère ; lors de chaque étape e), chaque spectre de regroupement est décomposé en un 20 même nombre de coordonnées. lors de l'étape f), l'estimation d'une caractéristique de l'objet est réalisée à partir des coordonnées obtenues dans chaque étape e). Le procédé peut comporter une répétition des étapes b) à e), de telle sorte que chaque étape b), comporte une mesure d'un spectre dans une position différente 25 du détecteur relativement à l'objet ; chaque étape c) est mise en oeuvre en considérant un même critère ; éventuellement, lors de chaque étape e), le spectre de regroupement est associé à un même nombre de coordonnées ; lors de l'étape f), l'estimation d'une caractéristique de l'objet est réalisée à partir 30 des coordonnées obtenues dans chaque étape e) L'étape f) peut comporter : - une discrétisation de l'objet en mailles bidimensionnelles ou tridimensionnelles ; 3061554 4 une estimation d'une atténuation dans chaque maille ou de la nature de matériau composant chaque maille à partir des coordonnées obtenues lors de chaque étape e). Selon un mode de réalisation, lors de l'étape c), le nombre de canaux du spectre de regroupement est égal à 2, de telle sorte que la formation du spectre de regroupement 5 comporte une détermination d'un seuil, dit seuil de regroupement, de telle sorte que le premier canal du spectre de regroupement rassemble les canaux élémentaires du spectre obtenu lors de l'étape b), correspondant à une énergie inférieure au seuil de regroupement, et que le deuxième canal de regroupement rassemble les canaux élémentaires dudit spectre obtenu correspondant à une énergie supérieure au seuil de regroupement. 10 Le nombre de canaux du spectre obtenu lors de l'étape b) peut être faible, par exemple compris entre 3 et 10. Le nombre de canaux du spectre de regroupement formé lors de l'étape c) est supérieur ou égal à 2. Un autre objet de l'invention est un dispositif de détection d'un rayonnement ionisant comportant : 15 un détecteur de rayonnement apte à interagir avec un rayonnement électromagnétique ionisant, chaque pixel comportant un circuit électronique apte à générer un signal dont l'amplitude dépend d'une énergie déposée par le rayonnement interagissant dans ledit détecteur; un circuit dit de spectrométrie, permettant d'obtenir une répartition de signaux générés 20 par le détecteur de rayonnement en fonction de leur amplitude, de façon à former un spectre représentatif de ladite répartition ; une unité de traitement apte à mettre en oeuvre les étapes c) à f) d'un procédé tel que décrit dans cette demande à l'aide des spectres formés par le circuit de spectrométrie. Le dispositif peut comporter différents détecteurs, chaque détecteur formant un pixel, de façon 25 à former un dispositif de détection pixelisé. D'autres avantages et caractéristiques ressortiront plus clairement de la description qui va suivre de modes particuliers de réalisation de l'invention, donnés à titre d'exemples non limitatifs, et représentés sur les figures listées ci-dessous. FIGURES 30 La figure 1A représente un dispositif d'analyse d'un objet selon l'invention. La figure 1B schématise les étapes d'un procédé d'analyse d'un objet selon l'invention. 3061554 5 La figure 1C représente une distribution en énergie d'un rayonnement. La figure 1D est un spectre obtenu à partir de la distribution représentée sur la figure 1C. La figure 1E est un exemple de formation d'un spectre de regroupement selon l'art antérieur. La figure 2A représente un spectre d'un rayonnement émis par une source d'irradiation et 5 détecté par un détecteur, sans objet entre la source d'irradiation et le détecteur. La figure 2B représente un spectre d'un rayonnement transmis par un premier objet, et détecté par un détecteur. La figure 2C représente le spectre, illustré sur la figure 2B, obtenu après un regroupement en deux canaux selon l'art antérieur. La figure 2D représente un spectre obtenu après un regroupement selon deux canaux effectué selon l'invention. 10 La figure 3A représente un spectre d'un rayonnement transmis par un deuxième objet, et acquis par un détecteur. La figure 3B représente le spectre, illustré sur la figure 3A, obtenu après un regroupement en deux canaux selon l'art antérieur. La figure 3C représente un spectre obtenu après un regroupement selon deux canaux effectué selon l'invention. La figure 4A représente un spectre d'un rayonnement transmis par un troisième objet, et acquis 15 par un détecteur. La figure 4B représente le spectre, illustré sur la figure 4A, obtenu après un regroupement en deux canaux selon l'art antérieur. La figure 4C représente un spectre obtenu après un regroupement selon deux canaux effectué selon l'invention. La figure 5 montre différents spectres de rayonnement correspondant à différentes configurations. 20 Les figures 6A et 6B représentent respectivement des épaisseurs de tissus et d'os d'un fantôme étalon modélisé. La figure 6C représente une distribution spatiale bidimensionnelle d'une valeur du seuil appliqué des spectres d'un rayonnement transmis par le fantôme. La figure 6D est un profil de la figure 6C. 25 La figure 7A montre différentes épaisseurs d'os du fantôme de calibration utilisé. Les figures 7B, 7C et 7D sont des profils montrant différentes les épaisseurs d'os théoriques ainsi que des épaisseurs d'os obtenues en regroupant des spectres mesurés selon l'art antérieur ou selon l'invention. EXPOSE DE MODES DE REALISATION PARTICULIERS 30 La figure 1A représente un mode de réalisation d'un dispositif 1 mettant en oeuvre un procédé selon l'invention. Une source d'irradiation 11 émet un rayonnement électromagnétique ionisant 12, dit rayonnement incident, vers un objet 10. L'objet 10 est disposé entre la source d'irradiation 11 3061554 6 et un dispositif de détection 20. Le dispositif de détection comporte des détecteurs de rayonnement 20i prenant la forme de pixels 20i agencés selon un plan, dit plan de détection P. Sur la figure 1A, on a représenté deux pixels adjacents 201 et 202. L'indice i désigne les coordonnées de chaque pixel dans le plan de détection. Les pixels peuvent s'étendre selon une 5 ligne mais en général, ils s'étendent selon une matrice régulière bidimensionnelle. L'objet 10 peut être un tissu biologique vivant, par exemple une partie du corps d'un animal ou d'un être humain. Le dispositif est alors un dispositif d'imagerie médicale. L'objet peut également être une pièce industrielle ou un bagage, le dispositif étant alors utilisé à des fins de contrôle non destructif. 10 Le terme rayonnement électromagnétique ionisant désigne un rayonnement électromagnétique constitué de photons d'énergie supérieure à 1 keV, et de préférence inférieure à 5 MeV. La plage d'énergie du rayonnement ionisant peut être comprise entre 1 keV et 2 MeV, mais elle s'étend le plus souvent entre 1 keV et 150 keV ou 300 keV. Le rayonnement ionisant peut être un rayonnement X ou y. De préférence, la source de rayonnement ionisant 15 est poly-énergétique, le rayonnement incident étant émis selon une plage d'énergie s'étendant généralement selon plusieurs dizaines voire centaines de keV. Il s'agit notamment d'un tube émetteur de rayons X. Ainsi, l'objet 10 irradié par la source 11 transmet au dispositif de détection 20 un rayonnement 14, dit rayonnement transmis, ce dernier atteignant les pixels, chacun formant un détecteur 20i. 20 Chaque pixel 20i est un détecteur spectrométrique comportant : un matériau détecteur, apte à interagir avec les photons du rayonnement 14 transmis par l'objet, ce matériau étant de type scintillateur ou, de préférence, un matériau semiconducteur compatible avec une utilisation à la température ambiante, de type CdTe, CdZnTe ; 25 un circuit électronique, apte à générer un signal dont l'amplitude dépend, et est de préférence proportionnelle, à l'énergie déposée par chaque photon interagissant dans le matériau détecteur ; un circuit de spectrométrie, apte à établir un spectre en énergie, noté Spi, des signaux détectés pendant une période temporelle, dite période d'acquisition. 30 Le dispositif de détection comporte également une unité de traitement 22 permettant la mise en oeuvre des opérations de traitement décrites ci-après. L'unité de traitement peut comprendre un microprocesseur et/ou des microcontrôleurs électroniques. 3061554 7 Selon un autre mode de réalisation, chaque pixel est couplé à un circuit électronique comportant des comparateurs d'amplitude permettant un regroupement des impulsions collectées par chaque pixel selon différents canaux. Selon un tel mode de réalisation, la résolution spectrale, c'est-à-dire la largeur des canaux, n'est pas aussi fine que dans les détecteurs spectrométriques 5 : le nombre de canaux est généralement compris entre 2 et 10. Un avantage de ce mode de réalisation est de pouvoir être implémenté directement dans un circuit électronique basé sur des discriminateurs d'amplitude. Sous l'effet de l'irradiation par le rayonnement incident 12, l'objet 10 transmet un rayonnement 14, dit rayonnement transmis, vers le dispositif de détection 20, dont chaque pixel 20i est apte 10 à détecter ledit rayonnement transmis 14 et à en former un spectre en énergie Spi. Le terme spectre en énergie correspond à un histogramme de l'amplitude des signaux détectés au cours d'une période d'acquisition du spectre. Une relation entre l'amplitude A et l'énergie E peut être obtenue par une fonction d'étalonnage en énergie g telle que E = g (A), selon des principes connus de l'homme du métier. Un spectre d'énergie Spi est donc un vecteur, dont 15 chaque terme Spi(h) représente une quantité de rayonnement détecté par le pixel 20i dans une plage d'énergie E ± -aE, avec dE étant la largeur spectrale de chaque canal. 2 La figure 1C représente une distribution en énergie d'un rayonnement et la figure 1D représente un spectre Spi discrétisé selon H canaux h. Dans les procédés d'analyse de l'art antérieur, comme précédemment évoqué, chaque spectre Spi fait l'objet d'un regroupement de canaux, 20 de manière à former un spectre regroupé Sp' i, chaque canal k du spectre regroupé regroupant l'information répartie dans des canaux adjacents du spectre initial Spi. Selon l'art antérieur, La valeur de chaque canal Sp'1(k) du spectre regroupé est constitué à partir du contenu des canaux Spi(h) du spectre mesuré Spi s'étendant entre deux bornes hmk-hink ax prédéterminées. Ainsi, la valeur de chaque canal Sp'1(k) est telle que Sp'1(k) = 25 hk zax Spi(h) (1). min Un exemple d'un tel regroupement est représenté sur la figure 1E. Quels que soient les pixels du détecteur, les bornes du regroupement n'évoluent pas. L'invention propose une approche différente, selon laquelle le regroupement d'un spectre dépend dudit spectre, c'est-à-dire de sa forme, sans que les bornes de regroupement soient 30 déterminées a priori, mais adaptées à chaque spectre. Cette approche est désignée par le terme 3061554 8 « regroupement adaptatif ». Pour cela, le regroupement est réalisé selon un critère C de regroupement prédéterminé. Selon un premier mode de réalisation, un tel critère est l'exigence que chaque canal k du spectre regroupé ait la même valeur Sp'1(k), ou sensiblement la même valeur, le terme sensiblement 5 signifiant la prise en compte d'une tolérance, par exemple de l'ordre de 10% ou 20%, tenant compte de la discrétisation en énergie du spectre. Un tel critère permet une répartition uniforme du bruit, supposé poissonien, affectant le spectre. Ainsi, chaque canal de regroupement Sp'1(k) est affecté d'un même bruit. En considérant un tel critère, Le procédé selon l'invention comporte une définition des bornes de regroupement hink ga,relative à 10 chaque canal k de telle sorte que quel que soit k, Sp'1(k) = E - Sp1(h) (2) et que pour hi,min tout couple de canaux km, Ici, du spectre de regroupement que Sp' i(km) = Sp' i(kii) (3). Chaque canal de regroupement k est constitué en appliquant une fonction de regroupement fk à des canaux h du spectre mesuré répartis dans une plage de regroupement Ak, la fonction de regroupement étant paramétrée pour de telle sorte que le spectre de regroupement satisfasse 15 au critère C. Selon ce mode de réalisation, la fonction de regroupement fk est une somme du contenu Sp1(h) des canaux du spectre mesuré, ces canaux faisant partie de la plage de regroupement Ak. Dans cet exemple, la plage de regroupement Ak s'étend entre deux bornes de regroupement, de telle sorte que Sp'1(k) = fk( Spi(h)), h E Ak (4). Dans cet exemple : - la plage de regroupement Ak de chaque canal de regroupement est constitué de canaux 20 adjacents s'étendant entre une borne minimale hink in et une borne maximale 14n, ; - la fonction de regroupement est une somme du contenu Sp1(h) de chaque canal h du spectre mesuré situé dans la plage de regroupement Ak. La plage de regroupement Ak, c'est-à-dire chaque borne de regroupement hlic,rnin, hkmax dépend du spectre mesuré Sp1(h). De ce fait, deux spectres Sp'1(k), Sp' j(k) respectivement 25 acquis par deux pixels différents i, j d'un détecteur situé face à un même objet peuvent être regroupés en considérant des plages de regroupement différentes. De préférence, le nombre K de canaux de regroupement est le même pour l'ensemble des pixels du détecteur. Ce nombre est par exemple prédéterminé. De même, le critère selon lequel s'effectue le regroupement est préférentiellement le même pour l'ensemble des pixels du détecteur. Selon une variante, le 30 nombre de canaux de regroupement varie d'un pixel à un autre, comme décrit par la suite, en lien avec le troisième mode de réalisation. 3061554 9 On a représenté, sur la figure 2A, un spectre Spi mesuré par un pixel 201. Il s'agit d'un spectre simulé à l'aide d'un code de calcul simulant le transport de photons. Il correspond au spectre mesuré par le pixel exposé à un rayonnement produit par une source d'irradiation 11, sans objet 10 disposé entre la source et le pixel. Un tel spectre est usuellement désigné par le terme "plein 5 flux". La source d'irradiation 11 modélisée est une anode de tungstène soumise à un potentiel de 110 kV. Le pixel 201 est un parallélépipède de CdZnTe de dimension 1 mm x 1 mm et d'épaisseur 500 um. Les pics observés vers 57-58 keV et entre 65 keV et 70 keV correspondent à des pics de fluorescence X du matériau détecteur et du tungstène. Le spectre de la figure 2B a été obtenu en simulant un objet 10 constitué de 5 cm de tissu entre la source et le détecteur. 10 Le tissu est constitué d'eau. La figure 2C représente un spectre de regroupement obtenu en regroupant les canaux h du spectre de la figure 2B en se basant sur un seuil arbitraire égal à 60 keV. Les impulsions d'énergie inférieure à 60 keV sont regroupées dans un premier canal de regroupement k = 1 tandis que les impulsions d'énergie supérieure à 60 keV sont regroupées dans un deuxième canal de regroupement k = 2. On observe que les contenus de chaque canal 15 ne sont pas identiques. La figure 2D montre le résultat d'un regroupement adaptatif, selon l'invention, en se basant sur un critère d'uniformité des canaux issus du regroupement. Selon un tel critère, après regroupement, chaque canal regroupé comporte un même nombre d'impulsions, à une tolérance près, la tolérance étant par exemple de ± 10% ou ± 20%. Ce critère a été appliqué au spectre représenté sur la figure 2B. L'analyse du spectre a déterminé que pour 20 satisfaire au critère de regroupement, il faut considérer un seuil, dit seuil de regroupement, de 56 keV. Le contenu des canaux h du spectre correspondant à une énergie inférieure à ce seuil sont regroupés dans un premier canal (k = 1), tandis que les canaux h supérieurs à ce seuil sont regroupés dans un deuxième canal (k = 2). On obtient ainsi un spectre de regroupement Spi , comportant deux canaux, dont chaque canal k a sensiblement la même valeur. Les seuils 25 de 60 keV et de 56 keV sont matérialisés sur la figure 2B par des lignes verticales en pointillés. La figure 3A est un autre exemple, dans lequel l'objet 10 disposé entre la source d'irradiation 11 et le pixel 201 est constitué d'une épaisseur de 20 cm de tissu. Cela correspond à une atténuation plus importante que celle de l'exemple précédent. La figure 3B représente un spectre de regroupement obtenu en regroupant les canaux h du spectre de la figure 3A en se basant sur un 30 seuil arbitraire égal à 60 keV. Les impulsions d'énergie inférieure à 60 keV sont regroupées dans un premier canal de regroupement k = 1 tandis que les impulsions d'énergie supérieure à 60 keV sont regroupées dans un deuxième canal de regroupement k = 2. On observe que les contenus de chaque canal ne sont pas identiques. La figure 3C montre le résultat d'un 3061554 10 regroupement adaptatif, selon l'invention, en se basant sur un critère d'uniformité des canaux issus du regroupement. Ce critère a été appliqué au spectre représenté sur la figure 3A. L'analyse du spectre a déterminé que pour satisfaire au critère de regroupement, il faut considérer un seuil de regroupement de 64 keV. Le contenu des canaux h du spectre correspondant à une 5 énergie inférieure à ce seuil sont regroupés dans un premier canal (k = 1), tandis que les canaux h supérieurs à ce seuil sont regroupés dans un deuxième canal (k = 2). On obtient ainsi un spectre de regroupement Spi , comportant deux canaux, dont chaque canal k a sensiblement la même valeur. Les seuils de 60 keV et de 64 keV sont matérialisés sur la figure 3A par des lignes verticales en pointillés. 10 La figure 4A est un autre exemple, dans lequel l'objet 10 disposé entre la source d'irradiation 11 et le pixel 201 est constitué d'une épaisseur de 5 cm de tissu et de 2 cm d'os. La figure 4B représente un spectre de regroupement obtenu en regroupant les canaux h du spectre de la figure 4A en se basant sur un seuil arbitraire égal à 60 keV. Les impulsions d'énergie inférieure à 60 keV sont regroupées dans un premier canal de regroupement k = 1 tandis que les impulsions 15 d'énergie supérieure à 60 keV sont regroupées dans un deuxième canal de regroupement k = 2. On observe que les contenus de chaque canal ne sont pas identiques. La figure 4C montre le résultat d'un regroupement adaptatif, selon l'invention, en se basant sur un critère d'uniformité des canaux issus du regroupement. Ce critère a été appliqué au spectre représenté sur la figure 4A. L'analyse du spectre a déterminé que pour satisfaire au critère de regroupement, il faut 20 considérer un seuil de regroupement de 65 keV. Le contenu des canaux h du spectre correspondant à une énergie inférieure à ce seuil sont regroupés dans un premier canal (k = 1), tandis que les canaux h supérieurs à ce seuil sont regroupés dans un deuxième canal (k = 2). On obtient ainsi un spectre de regroupement Spi , comportant deux canaux, dont chaque canal k a sensiblement la même valeur. Les seuils de 60 keV et de 65 keV sont matérialisés sur 25 la figure 4A par des lignes verticales en pointillés. Les trois exemples qui précèdent, illustrent la variabilité de la valeur d'un seuil de regroupement permettant le regroupement de canaux, en fonction de la constitution de l'objet observé, en considérant un même critère de regroupement. Ils justifient l'approche adaptative sur laquelle se base l'invention. Dans ces exemples, le nombre de canaux de regroupement est de 2 (K = 30 2), mais l'invention s'applique naturellement un nombre de canaux de regroupement plus élevé. 3061554 11 Selon un mode de réalisation, le critère de regroupement est identique à celui du premier mode de réalisation, c'est-à-dire une répartition uniforme du contenu des canaux de regroupement k, mais la fonction de regroupement fk, combinant les canaux regroupés h, est différente. Cette dernière est ici une somme pondérée du contenu Spi(h) des canaux h situés dans une plage de 5 regroupement Ak, le terme de pondération étant dépendant d'une valeur d'énergie E (h) à laquelle est associé chaque canal h du spectre mesuré. Ainsi, Sp' i(k) = fk( Spi(h)), h E Ak= Eh e A k E(h)-3Spi(h) Eax E(h)-3Spi(h) (5) h i,min Une telle pondération prend en compte le fait qu'en dessous d'une certaine énergie, le coefficient d'absorption d'un matériau varie en fonction de E-3. 10 Selon un troisième mode de réalisation, le critère de regroupement dépend des pixels, en se basant sur un a priori. Cet a priori peut résulter d'une caractérisation préalable de l'objet, permettant d'obtenir une première information sur la nature des matériaux dans l'objet. Le critère de regroupement tient alors compte de cet a priori. Par exemple, en se basant sur l'a priori, certains pixels sont susceptibles de comporter une raie K caractéristique d'un matériau 15 spécifique. Dans ces pixels, on peut prévoir des canaux de regroupement disposés de part et d'autre de la raie K. Dans les applications liées au diagnostic médical, on peut ainsi cibler les pixels correspondant aux parties d'un corps analysé dans lesquelles une présence d'un produit de contraste (par exemple Gadolinium ou Iode) est suspectée. Cela permet également d'appliquer des regroupements différents à des pixels correspondant à des tissus de différentes 20 natures : par exemple, on applique un premier regroupement aux pixels correspondant aux tissus gras, et un deuxième critère de regroupement aux pixels correspondant aux tissus maigres, ou aux pixels correspondant à des os, ou encore du métal. L'invention peut être mise en oeuvre de façon itérative : au cours d'une première itération, on applique un critère de regroupement identique à chaque spectre acquis, ce qui permet de disposer d'un a priori sur la 25 nature des matériaux composant l'objet. Lors d'itérations suivantes, le critère de regroupement appliqué à chaque pixel dépend de la nature des matériaux traversés par le rayonnement entre la source d'irradiation et le pixel. La figure 5 représente quatre spectres simulés en se basant sur une source d'irradiation et sur un détecteur tel que décrit en lien avec les exemples précédents. Les spectres représentés sur 30 la figure 5 correspondent aux configurations suivantes : spectre I : objet formé de 5 cm de tissu ; 3061554 12 spectre II : objet formé de 5 cm de tissu et de 1 cm de gadolinium ; spectre III : objet formé de 5 cm de tissu et de 1 cm d'os ; spectre IV : objet formé de 5 cm de tissu et de 1 cm d'os et de 1 cm de gadolinium. On observe que la présence de gadolinium entraîne une discontinuité de part et d'autre de 5 l'énergie 50 keV. Ainsi, le procédé peut comporter les étapes suivantes : au cours d'une première itération, formation de canaux de regroupement de part et d'autre de 50 keV, par exemple selon deux canaux correspondant à des énergies [45 keV - 50 keV] et [51 keV et 56 keV] ; comparaison des contenus des canaux ainsi regroupés ; 10 en cas de différence significative du contenu, l'objet observé contient du gadolinium. Le regroupement peut être réalisé selon un critère de regroupement déterminé, tout en imposant deux canaux de regroupement de part et d'autre de 50 keV. Par exemple un premier canal de regroupement [10 keV - 50 keV], un deuxième canal de regroupement [51 keV - 75 keV], et au moins un troisième canal de regroupement au-delà de 75 keV. 15 Les canaux au-delà de 75 keV peuvent faire l'objet d'un regroupement selon un critère de regroupement, par exemple un critère d'uniformité. lorsque la comparaison réalisée n'aboutit pas à une différence significative entre les contenus des canaux à [45 keV - 50 keV] et [51 keV et 56 keV], regroupement du spectre selon un critère de regroupement, par exemple uniformité. 20 Ainsi, le spectre obtenu fait l'objet d'un regroupement itératif, la présence d'un matériau particulier étant vérifiée au cours d'une première itération, la deuxième itération tenant compte des résultats de la première itération. Dans les modes de réalisation précédemment décrits, la plage de regroupement est composée de canaux adjacents, mais une plage de regroupement peut comporter des canaux non 25 adjacents les uns des autres. Chaque spectre regroupé Sp'1 représente une atténuation Att'1 du rayonnement électromagnétique ionisant 12 par l'objet. L'atténuation Att'1 peut ensuite être décomposée dans une base BK, cette dernière étant formée par N vecteurs de base 17,2K, de dimension K, de telle sorte que : Atei ai,'17,21( où ai,' désigne une coordonnée du spectre de 30 regroupement Att'1 dans la base. L'atténuation Att'1 correspondant au spectre regroupé Sp'1 est un vecteur de même dimension que le spectre regroupé Sp'1. La détermination des coordonnées de l'atténuation dans la base peut être obtenue en utilisant une méthode 3061554 13 d'estimation statistique, de type maximum de vraisemblance, ou par une méthode paramétrique de type polynomiale. Ces méthodes permettent d'estimer les ai,n en fonction des contenus des canaux du spectre regroupé Sp'i. Le symbole désigne le fait que la décomposition de l'atténuation dans la base est réalisée par estimation, du fait que la relation 5 entre le spectre regroupé et l'atténuation ne soit pas linéaire. Chaque vecteur 17,2K de la base BK est de la même dimension que le nombre de canaux K du spectre de regroupement Sp'i. La décomposition d'un spectre selon une base de coefficients d'atténuation est décrite dans la publication Rebuffel "Dual-energy X-ray imaging : benefits and limits" Insight Vol 49 No 10, october 2007 ou dans la publication Alvarez R "Near optimal energy selective x-ray imaging 10 system performance with simple detectors", Med. Phys. 37 (2), February 2010. Chaque vecteur de base Vir peut notamment être représentatif d'un spectre obtenu en disposant un matériau de calibration calibra connu à la place de l'objet à examiner, et faisant l'objet d'un regroupement similaire au regroupement ayant abouti au spectre regroupé Sp'i. Par regroupement similaire, on entend un regroupement effectué selon des mêmes plages de 15 regroupement Ak et selon une même fonction de regroupement fk que le spectre regroupé Sp'i. Dans ce cas, une coordonnée an représente une quantité équivalente du matériau de calibration calibra dans l'objet analysé. Par exemple, on peut disposer de deux matériaux de calibrations calib1 et calib2. Dans ce cas, on mesure successivement, à l'aide du détecteur 20i, deux spectres de calibration Sn rixalibl et 20 S'11 -ri,calib2, chaque spectre de calibration étant défini selon H canaux. Après l'acquisition d'un spectre Spi d'un objet à analyser, le spectre mesuré fait l'objet d'un regroupement pour former le spectre de regroupement Sp'i. On peut alors obtenir les coordonnées correspondant à ce spectre de regroupement dans une base dont chaque vecteur de base est obtenu en appliquant le même regroupement aux spectres de calibration Sp et Spi,calib2, de façon à former 25 deux spectres regroupés Spi,calib1et Spi',calib2, ces derniers formant respectivement les vecteurs de base 171 et 17 . La décomposition du spectre regroupé Sp'i dans la base est telle que : Att'i aiyiK + a2Vf , a1 et a2 représentant respectivement une proportion équivalente de chaque matériau de calibration dans le matériau à analyser. Par même regroupement, on entend un regroupement effectué selon la même plage de regroupement 30 Ak et en appliquant la même fonction de regroupement fk. Alternativement, chaque spectre de calibration peut être normalisé par un spectre Spi,o, dit spectre de la source, du rayonnement 12 émis par la source 11 sans objet disposé entre la source 3061554 14 et le détecteur 20, Le spectre Spi d'un objet à analyser peut alors être également normalisé par le spectre de la source Spi,o, avant de faire l'objet d'un regroupement, formant le spectre de regroupement Sp' i. La décomposition du spectre, ou plus exactement de l'atténuation correspondant à ce spectre, dans la base est telle que Ate + a2Vf , où Vii_c et 17f 5 sont obtenus après normalisation respective des spectres de calibration S Pi,calibl et SPi,calibl par Spi,o, puis regroupement des spectres de calibration ainsi normalisés. Les coordonnées a1 et a2 représentent respectivement une épaisseur équivalente de chaque matériau de calibration dans le matériau à analyser. Un vecteur de base peut également être représentatif d'un phénomène physique, par exemple 10 une discontinuité spectrale. Cela permet de prendre en compte des discontinuités de type absorption K d'un matériau irradié par la source d'irradiation 11, en particulier lorsque la source d'irradiation dispose d'un collimateur. Cela permet également la prise en compte de discontinuités d'absorption dans un spectre, de type raies K. Un vecteur de base peut également être représentatif d'un type d'atténuation dans l'objet 15 analysé. Par exemple, une base peut comporter un premier vecteur de base Vii_c représentatif d'une section efficace de l'absorption photoélectrique d'un matériau représentatif de l'objet à analyser, et un deuxième vecteur de base Vf représentatif d'une section efficace de la diffusion Compton dans ledit matériau. Dans ce cas, chaque coordonnée a1, a2 représente alors une proportion relative de chaque type d'atténuation (photoélectrique ou Compton) dans le spectre 20 analysé. D'une façon générale, un vecteur de base est une fonction spectrale définie selon les mêmes canaux k que le spectre regroupé Sp'i. L'atténuation Atei correspondant au spectre regroupé peut alors être décomposée selon la base. Chaque vecteur de base peut être obtenu à partir d'une fonction spectrale initialement établie selon un nombre important de canaux, par 25 exemple en autant de canaux que le spectre mesuré Spi, puis ayant fait l'objet d'un regroupement dans des canaux correspondant respectivement aux énergies de chaque plage de regroupement k. Par exemple, un vecteur de base peut être défini par une fonction spectrale définie sur autant de canaux que le spectre mesuré, après regroupement de la fonction spectrale selon la même plage de regroupement Ak et la même fonction de regroupement fk que celles 30 appliquées au spectre mesuré. Comme précédemment indiqué, deux pixels élémentaires i,j différents placés en regard d'un même objet 10, acquièrent respectivement deux spectres élémentaires Spi, Spi, ces derniers 3061554 15 faisant ensuite l'objet d'un regroupement de façon à former respectivement deux spectres de regroupement Sp'i,Sp' j. On obtient ensuite les coordonnées de chaque spectre de regroupement Sp'i,Sp' j dans une base formée à partir de vecteurs de base. Les vecteurs de base utilisés sont représentatifs des mêmes matériaux de calibration. Cependant, les 5 coordonnées du spectre Sp'1 peuvent être obtenues selon une base Br, formée de N vecteurs 171,c,i, tandis que les coordonnées du spectre Sp' j peuvent être obtenues selon une base BK, formée de N vecteurs Vlic,j. Les bases Br et edépendent respectivement du regroupement appliqué aux spectres Spi et Spi. Les bases dans lesquelles on détermine les coordonnées les spectres de regroupement, 10 correspondant à un même objet 10, peuvent avoir la même dimension, c'est-à-dire le même nombre K de vecteurs de base. Deux spectres de regroupement représentatifs d'un même objet peuvent être décomposés selon des bases différentes, ces dernières dépendant de la fonction de regroupement fk et des plages de regroupement Ak. Autrement dit, le regroupement appliqué à chaque spectre détecté Spi, c'est-à-dire les plages de regroupement Ak et la fonction 15 de regroupement fk, conditionne la base BK dans laquelle le spectre de regroupement Sp'1 est décomposé. Toutefois, comme chaque base BK peut comporter un même nombre N de vecteurs, chaque spectre de regroupement peut être caractérisé par un même nombre de coordonnées an. Chaque coordonnée aim, représente une proportion de chaque vecteur de base gr dans le 20 l'atténuation correspondant au spectre regroupé Sp'1, et le spectre regroupé associé à chaque pixel 20i est caractérisé par le même nombre N de coordonnées ai,n, cela quel que soit le regroupement appliqué à chaque pixel i. Les différentes étapes d'un exemple de procédé de traitement d'un spectre sont illustrées en lien avec la figure 1B. L'étape 100 correspond à l'obtention d'un spectre Spi par un détecteur 25 20i. Au cours d'une étape 110, on définit, à partir du spectre mesuré Spi, des plages de regroupement Ak, satisfaisant à un critère de regroupement C (par exemple une répartition uniforme dans les canaux de regroupement), la fonction de regroupement fk étant une donnée d'entrée, au même titre que le critère C et que le nombre de canaux K du spectre de regroupement. Au cours d'une étape 120, on constitue un spectre de regroupement Sp'1 en 30 prenant en compte les plages de regroupement Ak déterminées lors de l'étape 110 sachant le critère de regroupement C et la fonction de regroupement fk. 3061554 16 Les plages de regroupement Ak conditionnent également la sélection d'une base Br dans laquelle le spectre de regroupement est décomposé (étape 130), l'obtention des coordonnées du spectre de regroupement dans la base ainsi sélectionnée faisant l'objet de l'étape 140. On obtient alors un ensemble de N coordonnées aim, du spectre Sp'1 dans la base Br, ,N 5 représentant le nombre de vecteurs de la base, N étant commun à chaque base Br, de telle sorte qu'à chaque spectre Sp'1 est associé un même ensemble de N coordonnées ai,n. Les étapes 100 à 140 peuvent être réitérées pour différents pixel 20i d'un détecteur et/ou en déplaçant le détecteur par rapport à l'objet. De cette façon, on obtient autant de coordonnées ai,n, l'indice i représentant une position respective d'un détecteur par rapport à l'objet 10 examiné. L'objet analysé peut être décomposé en mailles élémentaires tridimensionnelles, ou voxels r. L'étape 150 est une étape de caractérisation de l'objet à partir des coordonnées établies dans l'étape 140. Par exemple, une reconstruction tomographique peut être réalisée, visant à caractériser une atténuation gr) de chaque voxel. La reconstruction est mise en oeuvre selon des algorithmes de reconstruction tomographique connus de l'homme du métier, visant à 15 inverser un modèle direct de type : S = M A où : S est un vecteur (I x N, 1), dont chaque terme représente une coordonnée aim, déterminée au cours de chaque étape 140, / représentant le nombre de spectres Spi 20 considérés. La dimension du vecteur S est I x N, c'est-à-dire le nombre de spectres détectés Spi multiplié par le nombre de coordonnées aim, déterminées pour chaque spectre, après regroupement de ce dernier. M est une matrice de projection modélisée, de dimension (I x N,R), R représente le nombre de voxels. 25 A est un vecteur (R, 1) dont chaque terme représente une caractéristique de l'objet, par exemple un coefficient d'atténuation gr) d'un voxel de l'objet ; A peut également représenter une autre caractéristique de l'objet, par exemple une estimation de la nature du matériau composant cet objet. Selon un mode de réalisation, durant une première phase, le procédé est mis en oeuvre en 30 regroupant chaque spectre mesuré selon un critère prédéterminé. Les spectres regroupés sont alors décomposés selon une même base. L'objet est décomposé en mailles élémentaires r. A partir de chaque coordonnée issue de la décomposition des spectres regroupés, on estime une caractéristique de chaque maille élémentaire de l'objet, par exemple un coefficient 3061554 17 d'atténuation ou la nature du matériau composant chaque maille. A partir de cette reconstruction, on détermine un critère de regroupement C. Ce critère est par exemple un critère de séparation entre des bandes d'énergies définies grâce à la reconstruction. Le procédé se poursuit par : 5 un regroupement de chaque spectre mesuré par le critère de regroupement déterminé suite à la reconstruction précédemment réalisée ; une décomposition de chaque spectre ainsi regroupé dans une base ; une reconstruction de l'objet en fonction des coordonnées obtenues lors de la décomposition de chaque spectre dans la base. 10 L'invention peut également être mise en oeuvre en géométrie bidimensionnelle, c'est-à-dire en mode radiographie, en discrétisant l'objet en éléments de surface, ou pixels, le vecteur A obtenu représentant l'atténuation de chaque pixel. L'invention a été mise en oeuvre de façon expérimentale dans une application de radiologie. On a simulé des spectres mesurés en disposant un objet composite 10, formé par du tissu et des os, 15 entre la source d'irradiation 11 préalablement décrite et un détecteur 20 pixellisé, chaque pixel 20, étant un pavé de CdZnTe de dimension 1 mm x 1 mm et d'épaisseur 500 um. Le spectre plein flux de la source d'irradiation 11 est celui représenté sur la figure 2A. Les figures 6A et 6B montrent les composants de l'objet composite simulé. Cet objet est formé de quatre bandes de tissu, de longueur 20 cm et de largeur 5cm, et d'épaisseur respectivement égales à 5 cm, 10 cm, 20 15 cm et 20 cm. Ces bandes de tissu sont représentées sur la figure 6A, le niveau de gris indiquant leur épaisseur. L'objet comporte également des cylindres d'os, de diamètre 4 cm, et d'épaisseurs égales à 1 cm, 2 cm et 4 cm. Ces cylindres sont représentés sur la figure 6B, le niveau de gris indiquant leur épaisseur. Les cylindres d'os représentés sur la figure 6B ont été déposés sur les bandes de tissu 25 représentées sur la figure 6A. On a ensuite acquis différents spectres Spi, par chaque pixel 20, correspondant aux différentes parties de l'objet composite ainsi constitué. Chaque spectre a fait l'objet d'un regroupement pour former un spectre de regroupement SK dont chaque canal k a sensiblement la même valeur Spak). Dans cet exemple, chaque spectre de regroupement comporte 2 canaux, définis de part et d'autre d'un seuil, dit seuil de regroupement. La figure 6C 30 représente une distribution spatiale d'un tel seuil de regroupement, le niveau de gris représentant la valeur du seuil de regroupement. On observe que la valeur du seuil de regroupement varie en fonction des différentes parties de l'objet. La figure 6D représente un 3061554 18 profil montrant l'évolution du seuil de regroupement selon une ligne horizontale en pointillés représentés sur la figure 6C. Cette figure permet d'apprécier la variabilité du seuil en fonction de la composition de l'objet. On a représenté, sur la figure 6C, quelques valeurs particulières du seuil. 5 La figure 7A représente les différentes épaisseurs d'os formant le tissu composite. A partir de chaque spectre Spi acquis par chaque pixel 20, on a estimé une épaisseur d'os en considérant un spectre de regroupement SK formé : soit sur la base d'un seuil de regroupement fixe, égal à 60 keV soit sur la base d'un seuil de regroupement tel que représenté sur la figure 6C. 10 Les figures 7B, 7C et 7D représentent un profil de l'épaisseur estimée, le long des lignes horizontales représentées sur la figure 7A, respectivement du haut vers le bas. Sur chaque figure, on a représenté : un profil obtenu en mettant en oeuvre un seuil de regroupement fixe (profil 1), ce qui représente l'art antérieur. 15 un profil obtenu en mettant en oeuvre un seuil de regroupement tel que représenté sur la figure 6C (profil 2); un profil correspondant à la réalité (profil 3) On observe que sur chaque profil, les estimations d'épaisseurs mettant en oeuvre l'invention (profils 2) aboutissent à une erreur d'estimation moindre que les estimations d'épaisseur 20 obtenues sur la base d'un spectre de regroupement obtenu sur la base d'un seuil de regroupement fixe. L'invention permet d'améliorer la qualité de l'estimation quant à la nature de l'objet observé. REVENDICATIONS1. Procédé d'analyse d'un objet (10) comportant les étapes suivantes : a) irradiation de l'objet à l'aide d'une source de rayonnement (11), ladite source émettant un rayonnement électromagnétique ionisant (12) se propageant vers l'objet ; b) obtention d'un spectre (Spi) d'un rayonnement transmis (14) par l'objet à l'aide d'un détecteur de rayonnement (20i) appartenant à un dispositif de détection (20), le spectre obtenu comportant plusieurs canaux élémentaires (h), chaque canal élémentaire étant associé à une bande d'énergie (8E (h)); c) à partir du spectre obtenu (Spi), formation d'un spectre de regroupement (Sp' i) constitué de plusieurs canaux (k), dits canaux de regroupement, au moins un canal de regroupement étant obtenu en regroupant, selon un critère de regroupement prédéterminé (C), plusieurs canaux élémentaires (h), la valeur du spectre de regroupement (Spak)) en chaque canal de regroupement (k) étant fonction de la valeur du spectre (Spi(h)), obtenu lors de l'étape b), en chaque canal élémentaire (h) ainsi regroupé ; d) sélection d'une base (Br) comportant des vecteurs de base (e), dont chaque terme (e) représente une atténuation du rayonnement électromagnétique ionisant dans une plage d'énergie (AE(k)) associée à chaque canal de regroupement (k); e) obtention de coordonnées (ai,n) d'une atténuation, correspondant au spectre de regroupement (Sp') formé lors de l'étape c), dans la base (Br) sélectionnée lors de l'étape d) ; f) estimation d'une caractéristique de l'objet en fonction des coordonnées (ai,n) obtenues lors de l'étape e). 2. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape c) comporte, pour chaque canal de regroupement (k), les sous-étapes suivantes : en fonction du critère de regroupement (C), une association d'une plage de regroupement (Ak), définissant les canaux élémentaires (h) regroupés dans le canal de regroupement ; la prise en compte d'une fonction de regroupement (fk), de telle sorte que la valeur du spectre de regroupement (Spak)) au canal de regroupement (k) soit déterminé en 3061554 20 appliquant la fonction de regroupement (fk) à la valeur du spectre (Spi(h)), obtenu lors de l'étape b), dans chaque canal élémentaire (h) ainsi regroupé. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape c), la valeur attribuée à chaque canal de regroupement (Sp'(k)) est une somme ou une somme 5 pondérée de la valeur du spectre (Sp(h)), obtenu lors de l'étape b), en chaque canal élémentaire regroupé. 4. Procédé selon la revendication 3, dans lequel la valeur attribuée à chaque canal de regroupement (Sp'(k)) est une somme pondérée dans laquelle le contenu de chaque canal élémentaire est affecté d'un coefficient de pondération dépendant de la bande d'énergie (E (h)) 10 associée audit canal élémentaire (Sp(h)). 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape c), le critère de regroupement (C) est un critère d'uniformité, selon lequel les valeurs respectives de chaque canal de regroupement (Sp'(k)) composant le spectre de regroupement (Sp') sont égales, à une tolérance prédéterminée près. 15 6. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape d), chaque terme des vecteurs de base (le) est défini pour chaque canal de regroupement et correspond, dans la plage d'énergie (AE(k)) associée audit canal de regroupement (k), soit à une atténuation d'un matériau de calibration ( calibra) ; soit à une proportion relative d'atténuation photoélectrique ou de diffusion Compton 20 dans ladite plage d'énergie; soit à une discontinuité spectrale. 7. Procédé selon l'une quelconque des revendications précédentes, dans lequel le dispositif de détection (20) comporte une pluralité de détecteurs de rayonnements (20i, 20g) adjacents les uns des autres, les étapes b) à e) étant mises en oeuvre pour chaque détecteur de rayonnement, 25 de telle sorte que : chaque étape c) est mise en oeuvre en considérant un même critère (C); lors de chaque étape e), chaque spectre de regroupement (Spi Sp' i) est décomposé en un même nombre de coordonnées. lors de l'étape f), l'estimation d'une caractéristique de l'objet est réalisée à partir 30 des coordonnées obtenues dans chaque étape e). 3061554 21 8. Procédé selon l'une quelconque des revendications précédentes, comportant une répétition des étapes b) à e), de telle sorte que - chaque étape b), comporte une obtention d'un spectre (Spi) dans une position différente du détecteur relativement à l'objet ; 5 - chaque étape c) est mise en oeuvre en considérant un même critère (C) ; lors de chaque étape e), le spectre de regroupement (Sp' i) est associé à un même nombre de coordonnées ; lors de l'étape f), l'estimation d'une caractéristique de l'objet est réalisée à partir des coordonnées obtenues dans chaque étape e) 10 9. Procédé selon l'une quelconque des revendications 7 ou 8, dans lequel l'étape f) comporte : - une discrétisation de l'objet (10) en mailles bidimensionnelles (r) ou tridimensionnelles ; - une estimation d'une atténuation dans chaque maille ou de la nature de matériau composant chaque maille à partir des coordonnées (ai,n) obtenues lors de chaque étape 15 e). 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape c), le nombre de canaux du spectre de regroupement est égal à 2, de telle sorte que la formation du spectre de regroupement comporte une détermination d'un seuil, dit seuil de regroupement, de telle sorte que le premier canal du spectre de regroupement rassemble les canaux 20 élémentaires correspondant à une énergie inférieure au seuil de regroupement, et que le deuxième canal du spectre de regroupement rassemble les canaux élémentaires correspondant à une énergie supérieure au seuil de regroupement. 11. Dispositif de détection d'un rayonnement ionisant comportant : - un détecteur de rayonnement (20,) apte à interagir avec un rayonnement 25 électromagnétique ionisant, chaque pixel comportant un circuit électronique apte à générer un signal dont l'amplitude dépend d'une énergie déposée par le rayonnement interagissant dans ledit détecteur ; - un circuit dit de spectrométrie, permettant d'obtenir une répartition de signaux générés par le détecteur de rayonnement en fonction de leur amplitude, de façon à former un 30 spectre représentatif de ladite répartition ; 3061554 22 une unité de traitement (22) apte à mettre en oeuvre les étapes c) à f) d'un procédé selon l'une quelconque des revendications 1 à 10 à l'aide des spectres formés par le circuit de spectrométrie.