Brevet FR2632749A1 — publié le 19891215

PROCEDE ET SYSTEME D'ETALONNAGE D'UN SCANNER A RAYONS X EN UTILISANT UN OU PLUSIEURS ETALONS CIRCULAIRES EXCENTRES L'invention concerne les scanners à rayons X et, plus particulièrement un procédé d'étalonnage d'appareils de ce type qui utilise un ou plusieurs étalons circulaires excentrés par rapport à l'axe de-rotation du scanner; elle concerne également un système qui permet de mettre en oeuvre ce procédé. Pour examiner un patient, on utilise de plus en plus des appareils à rayons X appelés "scanners" qui réalisent des images de coupes transversales du patient. Ces appareils sont basés sur le phénomène physique d'absorption des rayons X par le corps humain. Pour un faisceau monochromatique, cette absorption est directement liée à la distance parcourue x. des rayons X dans un corps homogène selon la formule: I - ioe - bx formule dans laquelle: Io est l'intensité du rayonnement entrant dans le corps humain, I est l'intensité du rayonnement sortant du corps humain, b est un coefficient d'atténuation qui dépend du corps traversé. Dans une échelle de mesure logarithmique, l'atténuation I/Io est égale à bx, c'est-à-dire qu'elle est proportionnelle à la distance x. Ces appareils sont constitués essentiellement, comme le montre la figure 1, d'une source de rayons X référencée qui est associée à un dispositif de détection, 11, ces deux éléments étant disposés l'un par rapport à l'autre dans une relation géométrique fixe de manière à pouvoir intercaler entre eux le corps à examiner. En outre, ils sont supportés par une structure (non représentée) qui peut tourner autour du corps à examiner de manière à irradier le corps suivant des angles différents. La source de rayons X, qui est commandée par un dispositif 13, émet ses rayons suivant un secteur angulaire ayant une largeur suffisante pour illuminer toute la section transversale du corps. Le dispositif de détection 11 a la forme d'un secteur annulaire dont la longueur est adaptée à la largeur du faisceau de rayons X et est constitué d'un grand nombre de détecteurs élémentaires 12 juxtaposés les uns à côté des autres. Pour obtenir une image de la section transversale du corps humain traversé par le faisceau de rayons X, on fait tourner la structure de support de la source 10 et du dispositif de détection 11 autour du corps et on mesure les signaux de sortie des détecteurs élémentaires 12 pour les traiter de manière appropriée selon des procédés connus afin d'en tirer une image représentative de la section transversale. Pour ce traitement, les détecteurs élémentaires 12, également appelés canaux, sont connectés à un dispositif électronique 14 qui effectue en premier lieu un calcul du logarithme des signaux reçus de manière à obtenir un signal dont l'amplitude est proportionnelle à l'atténuation des rayons X. Par suite de l'effet de différents phénomènes qui ne seront pas expliqués ici, l'amplitude de ce signal pour chaque détecteur élémentaire ou canal n'est pas proportionnelle à l'atténuation effectivement subie. En conséquence, pour remédier à cet inconvénient, il a été imaginé divers procédés qui consistent par exemple à relever les signaux de sortie des canaux en présence de corps dont on connaît les dimensions et le coefficient d'absorption, de manière à calculer les attenuations (calculs de logarithmes) et à comparer ces atténuations mesurées à des valeurs calculées en fonction des dimensions et du coefficient d'absorption du corps ou étalon. Ces comparaisons permettent de déduire une loi de correspondance ou une loi modificatrice entre les valeurs mesurées et les valeurs que l'on devrait obtenir. Cette loi peut être sous la forme de fichiers de correspondance ou de formules mathématiques traduisant cette correspondance pour chaque canal de détection. Les étalons qui sont utilisés pour effectuer ces mesures dites d'étalonnage sont par exemple des cales d'épaisseurs différentes qui sont introduites à proximité de la source de rayons X, ce qui implique des manipulations au niveau de la source pour introduire et retirer ces cales. Dans le brevet US N 4352020, il est proposé d'utiliser des cales ou étalons circulaires 15, 16 et 17, de différents diamètres, qui sont disposées au centre de rotation de la structure du support. Ceci permet de se rapprocher des conditions des mesures qui seront effectuées sur le corps à examiner. Ce brevet propose également d'utiliser un étalon en forme de cône à section circulaire qui est déplacé transversalement par rapport au faisceau de manière à obtenir différentes longueurs d'atténuation. Avec les étalons décrits, les mesures sont effectuées pour une position déterminée de la structure de support et par étalon. La figure 2 montre l'allure de trois courbes de réponse , 21 et 22 de l'atténuation en fonction de la position des canaux dans le cas -de mesures sur trois étalons de forme circulaire. Les valeurs mesurées sont représentées par les pointillés et varient autour d'une valeur moyenne qui représente la valeur théorique dans un système linéaire. Ces courbes peuvent être utilisées de la manière suivante; lorsque le signal mesuré correspond à un point A, on en déduira que le signal linéaire est le point A' de la courbe moyenne 20. Lorsque le signal mesuré correspond e un point B situé entre les courbes 20 et 21, on en déduira le signal linéaire par interpolation entre les courbes 20 et 21. Cette interpolation peut être calculée suivant une loi linéaire ou plus généralement une loi polynomiale. Les courbes 23 et 24 de la figure 3 montrent sous une autre forme le principe de l'étalonnage au niveau d'un canal. Il s'agit de courbes décrivant dans un canal donné l'atténuation en fonction de l'épaisseur x, pour des valeurs mesurées (courbe 23) et pour des valeurs calculées (droite 24): en fait, les valeurs mesurées donnent des points, points que l'on relie entre eux selon une loi choisie: linéaire ou polynomiale, de manière à obtenir une courbe continue. Lorsqu'on mesure une atténuation, cela correspond par exemple au point C de la courbe 23 et on en déduit la valeur linéaire correspondant au point C' de la courbe 24. Le brevet américain précité décrit un appareil dans lequel la correspondance entre les valeurs mesurées et les valeurs réelles d'atténuation est effectuée par un système de fichiers créés pendant l'opération d'étalonnage. En ce qui concerne l'interpolation, le brevet propose des interpolations linéaires, cubiques et biquadratiques mais seule l'interpolation linéaire est décrite de manière détaillée. Les procédés d'étalonnage qui ont été succinctement décrits ci-dessus présentent l'inconvénient majeur de mettre en oeuvre plusieurs étalons, ce qui conduit à de nombreuses manipulations; en outre, ces manipulations doivent être précises, notamment dans le.cas d'étalons circulaires dont les différents centres doivent coïncider avec le centre de rotation de la structure. Il est à noter que le brevet américain propose d'utiliser un étalon unique qui aurait la forme donnée par la figure 12 dudit brevet, et de faire tourner la structure autour de cet étalon,.ce qui conduit à obtenir des trajets d'absorption de longueur différente selon la position angulaire de la structure. Cependant, une telle manière de faire n'est que citée et n'indique pas le procédé ni les moyens de le mettre en oeuvre dans ce cas. En résumé, le but des procédés d'étalonnage est d'homogénéiser la réponse des canaux entre eux et à tendre vers une réponse linéaire. Pour obtenir cette réponse linéaire, une correction polynomiale est appliquée à chaque canal et le calcul de cette correction, c'est-à-dire la détermination des coefficients du polynôme, s'effectue à partir de mesures effectuées avec des valeurs différentes du trajet d'absorption. C'est ainsi que pour déterminer les coefficients d'un polynôme d'ordre n d'un canal déterminé, il est nécessaire et suffisant d'effectuer (n+l) mesures de trajets différents pour ce canal.' Dans les exemples de réalisation de l'art antérieur décrits ci-dessus, il a été montré que l'on utilise par exemple (n+1) cales d'épaisseurs différentes ou (n+l) cales circulaires de diamètres différents, les valeurs des épaisseurs et des diamètres étant choisis pour couvrir la dynamique requise. Il en résulte donc de nombreuses manipulations pour passer d'une cale à une autre, manipulations qui doivent être précises dana le cas de cales circulaires dont les centres doivent coïncider avec le centre de la structure. Un but de la présente invention est donc de mettre en oeuvre un procédé d'étalonnage qui utilise un seul étalon circulaire dont la position est excentrée par rapport à l'axe de rotation du scanner. Un autre but de la présente invention est de réaliser un système pour mettre en oeuvre ledit procédé. L'invention se rapporte à un procédé d'étalonnage d'un scanner à rayons X, qui comporte une source de rayons X et un dispositif de détection à N détecteurs ou canaux portés par une structure qui tourne autour d'un axe, en utilisant un étalon circulaire, caractérisé en ce qu'il comprend les opérations suivantes: - la mise en place d'un étalon circulaire dont le centre est placé à une position excentrée par rapport à l'axe de rotation du scanner, - le déplacement du scanner pour réaliser m mesures d'atténuation Yij par canal de rang i, chaque mesure correspondant à une position angulaire dj du scanner, - le calcul pour chaque mesure Yij de la distance Xij parcourue par le rayonnement dans l'étalon, et - le calcul pour chaque canal de rang i, des coefficients Aik du polynôme Pn de degré n (n correspondantes b*Xij. - des moyens pour comparer les valeurs de Yij avec les valeurs calculées des atténuations b.Xij pour le même canal et calculer les coefficients Aik de l'approximation polynomiale, et, - des moyens pour enregistrer lesdits coefficients Aik ainsi calculées de manière à les utiliser ensuite comme coefficients multiplicateurs des atténuations Yij subies par le rayonnement X dans le corps du patient. D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description suivante d'un exemple particulier de réalisation, ladite description étant faite en relation avec les dessins joints dans lesquels: - la figure 1 est un schéma de principe d'un scanner à rayons X qui montre un certain nombre d'étalons circulaires utilisés pour l'étalonnage, selon l'art antérieur, - la figure 2 est un diagramme montrant différentes courbes d'atténuation en fonction de la position des détecteurs ou canaux et pour plusieurs valeurs du diamètre; - la figure 3 est un diagramme montrant les courbes d'atténuation théorique et mesurée pour un canal déterminé en fonction du trajet d'absorption, - la figure 4 est un schéma essentiellement géométrique qui permet de comprendre le calcul des trajets du rayonnement X pour un détecteur déterminé lorsque, selon l'invention, on utilise un étalon circulaire excentré, et - la figure 5 est un schéma fonctionnel d'un système de traitement des signaux de sortie des détecteurs mettant en oeuvre le procédé selon l'invention. Les figures 1, 2 et 3 qui ont permis de définir l'art antérieur dans le préambule ne seront pas décrites à nouveau. Sur la figure 4, on a représenté de manière schématique mais précise les positions d'une source de rayonnement X matérialisé par le point F, d'un dispositif de détection 31 et d'un étalon circulaire 32 de rayon R disposé entre les deux premiers éléments 30 et 31. La source 30 et le dispositif de détection 31 sont portés par une structure, non représentée, qui tourne autour d'un axe perpendiculaire au plan du dessin au point O. Le centre E de l'étalon circulaire 32 est situé à une distance r du point O et la demi-droite OE fait un angle { avec l'axe Ox d'un système d'axes perpendiculaires x'x-y'y. La position d'ordre j de la structure de support est définie par l'angle Pj entre l'axe Oy et la demi-droite OF (position de la mesure ou de la vue) tandis que la position d'un détecteur I d'ordre i du dispositif de détection 31 est définie par l'angle 6i entre FO et FI. Ce schéma géométrique permet de calculer le trajet AB appelé Xij d'un rayonnement X dans l'étalon circulaire excentré 32 selon la formule Xij = AB = 2 /R2-. r sinSi - r cos ( formule (2). On remarquera que le nombre m est largement supérieur à n+l afin de tenir compte du bruit des mesures. Pour mettre en oeuvre ce procédé d'étalonnage, l'invention propose un système qui sera maintenant décrit en relation avec la figure 5. Il comprend un convertisseur analogique/numérique 50 auquel sont appliqués les signaux de sortie CI à CN des N détecteurs ou canaux du dispositif de détection. Les N codes numériques correspondant à une position ou vue de rang j sont appliqués à un circuit de calcul de logarithme 51 qui fournit, pour chaque canal i et chaque vue j, un code représentatif de l'atténuation subie dans l'étalon. Les N codes résultant de cette opération de calcul logarithmique sont appliqués à un circuit de soustraction 52 dans lequel on leur soustrait une valeur REF représentative de l'atténuation subie par le rayonnement X à l'extérieur de l'étalon, c'est-à-dire dans l'air. Cette valeur est obtenue à l'aide d'un détecteur appelé moniteur dont la position sur le dispositif de détection 11 est telle qu'il reçoit le rayonnement X sans atténuation par l'étalon. Ce sont les codes qui résultent de cette soustraction qui constituent les mesures Yij définies ci-dessus. Les N valeurs correspondant à une vue de rang j sont enregistrées dans deux mémoires, l'une référencée 53 pour le calcul des coefficients Aik conformément au procédé d'étalonnage et l'autre référencée 62, pour être corrigée en dehors des opérations d'étalonnage, à l'aide des coefficients Aik- Chaque vue correspondant à un angle Pj différent, N valeurs Yij sont enregistrées dans les mémoires 53 et 62. La mémoire 53 est par exemple prévue pour enregistrer les N.m codes correspondant à m-vues si 1l'on souhaite effectuer une approximation polynomiale d'ordre n (n analogue à la mémoire 53. Les codes contenus dans les mémoires 53 et 61 sont utilisés dans un circuit 58 pour calculer les coefficients polynomiaux Aik pour chaque canal de rang i en appliquant la formule (2) afin d'obtenir une valepr minimale pour Ci. Ces coefficients sont enregistrés dans une mémoire 59 de manière à être utilisés dans un circuit 63 pour corriger les valeurs mesurées en présence du corps du patient. Le calcul des coordonnées r et { du centre E de l'étalon peut être réalisé de différentes manières. L'une d'entre elles est basée sur le fait que l'atténuation la plus forte due à l'étalon correspond à un trajet diamétral, c'est-à-dire passant par le centre E. Ainsi, pour la- vue d'ordre j correspondant à Pj=O (OF confondu avec l'axe Oy), on détermine le canal i dans lequel l'atténuation est maximale, ce qui conduit à connaître l'angle a' entre l'axe Oy et la demi-droite FE. De la même manière, on peut déterminer l'angle a" entré l'axe Ox' et la demi- droite FE pour la vue correspondant à pj=90' (ou o). La connaissance de ces deux angles a' et a" associée à celle de la distance FO permet de calculer par les relations trigonométriques les coordonnées cartésiennes du centre E de l'étalon par rapport aux axes x'x et y'y. Une transformation de coordonnées cartésiennes en coordonnées polaires permet d'obtenir r et À. REVENDICATIONS 1. Procédé d'étalonnage d'un scanner à rayons X qui comporte une source de rayons X et un dispositif de détection à N détecteurs ou canaux portés par une structure qui tourne autour d'un axe, en utilisant un étalon circulaire, caractérisé en ce qu'il comprend les opérations suivantes: la mise en place d'un étalon circulaire dont le centre est placé à une position excentrée par rapport à l'axe de rotation du scanner, - le déplacement du scanner pour réaliser n mesures d'atténuation Yij par canal de rang i, chaque mesure correspondant à une position angulaire hj du scanner, - le calcul pour chaque mesure Yij de la distance Xij parcourue par le rayonnement X dans l'étalon,-et - le calcul pour chaque canal de rang i, des i coefficients Aik du polynôme P de degré n tel que: n ij k=n k P (Yij) = E Aik Yij n k=o de manière à obtenir une approximation polynomiale des valeurs calculées Xij suivant une loi d'écart déterminée. 2. Procédé selon la revendication 1, caractérisé en ce que la distance Xij est calculéepar la formule: JR2-( r sin6i - r cos (I-j-6i))2 (1) 3. Procédé selon la revendication 1 ou 2 caractérisé en ce que le calcul des coefficients Aik du polynme P (Y) ce que le calcul des coefficients Aik du polynôme P (Y) n est effectué pour que l'écart quadratique: m n k 2 Ci= ú { (Aik Yij) - b- Xij) (2) j=l k=o soit aussi faible que possible. 4. Système pour mettre en oeuvre le procédé selon l'une quelconque des revendications 1, 2 ou 3 caractérisé en ce qu'il comprend: - des moyens pour disposer au moins un étalon circulaire entre la source de rayons X et le dispositif de détection dans une position excentrée par rapport à l'axe de rotation de la structure tournante, - des moyens pour déplacer le scanner à m positions angulaires différentes Pj proches les unes des autres, - des moyens (51,52) pour calculer, à partir des signaux reçus sur les N détecteurs, les atténuations Yij subies par le rayonnement X dans l'étalon - des moyens (53) pour enregistrer les N valeurs Yij pour chacune des m positions angulaires pj, - des moyens (54) pour calculer à partir des N.m valeurs Yij les éléments du décentrage r et I de l'étalon par rapport à l'axe de rotation du scanner, - des moyens (55) pour calculer, à partir des éléments r et e du décentrage, les longueurs Xij des trajets du rayonnement X dans l'étalon ainsi que les atténuations correspondantes, b'Xij, - des moyens pour comparer les valeurs de Yij avec les valeurs calculées des atténuations b-Xij pour le même canal, et calculer les coefficients Aik de l'approximation polynomiale, et - des moyens pour enregistrer lesdits coefficients Aik -ainsi calculés de manière à les utiliser ensuite comme coefficients multiplicateurs des atténuations Yij subies par le rayonnement X dans le corps du patient.