Brevet FR2629214A1 — publié le 19890929

PROCEDE ET SYSTEME D ETALONNAGE D UN SCANNER A RAYONS X EN UTILISANT UN SEUL ETALON NON CIRCULAIRE L'invention concerne les scanners à rayons X et plus particulièrement un procédé d'étalonnage d'appareils de ce type qui utilise un seul étalon non circulaire et un système de mise en oeuvre dudit 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. Cette absorption est directement liée à la distance parcourue x des rayons X dans le corps 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 10 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 à rayons X, qui est commandée par un dispositif 13, émet ses rayons suivant un secteur angulaire qui a 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 applelé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 connait les dimensions et le coefficient d'absorption, de manière à calculer les atténuations (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. En outre, la forme de ces cales et leur position sont éloignées de la forme et de la position du corps du malade à examiner, ce qui augmente la non-linéarité du système. Dans le brevet USP N 4352020, il est proposé d'utiliser des cales circulaires 15 à 18, 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 & 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'apaisseur 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, celà 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 coincider 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 mettre en oeuvre dans ce cas. Un but de la présente invention est de mettre en oeuvre un procédé d'étalonnage qui utilise un étalon unique ayant une forme autre que circulaire, notamment une forme elliptique ou une forme qui serait symétrique par rapport à deux axes coplanaires. 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 rayonnement X et un dispositif de détection à N canaux, en utilisant un seul étalon de forme autre que circulaire, caractérisé en ce qu'il comprend les opérations suivantes: a) La mise en place de l'étalon entre la source de rayonnement X et le dispositif de détection; b) la mesure des atténuations dans les N canaux du dispositif de détection pour P positions angulaires principales de la structure tournante et pour n positions élémentaires proches les unes des autres autour de chaque position angulaire principale; c) le calcul de la moyenne des n atténuations pour chaque canal et pour chacune des P positions angulaires principales; d) le lissage des N valeurs moyennes des atténuations d'un canal au suivant pour chacune des P positions angulaires principales de manière à obtenir une courbe de réponse de l'atténuation en fonction de la position du canal dans laquelle les composantes haute fréquence ont été éliminées, ladite courbe de réponse ainsi obtenue étant la courbe réelle de variation de l'atténuation introduite par l'étalon en fonction de la position respective des N canaux. L'invention se rapporte également à un système pour mettre en oeuvre le procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend: - des premiers moyens pour calculer le logarithme des - N signaux fournis par les N canaux, - des deuxièmes moyens pour calculer, pour chaque canal, l'atténuation introduite par l'étalon, - des troisièmes moyens pour calculer, pour chaque canal, la valeur moyenne des n atténuations mesurées pour une position angulaire principale donnée, - des quatrièmes moyens pour calculer la transformée de Fourier des N atténuations moyennes correspondant à une position angulaire principale, des cinquièmes moyens pour éliminer les composantes haute fréquence dans le spectre résultant de la transformée de Fourier, - des sixièmes moyens pour calculer la transformée de Fourier inverse des composantes basse fréquence du spectre, - des sixièmes moyens pour calculer une approximation polynomiale de la courbe résultant de la transformée de Fourier inverse, - et des moyens pour garder en mémoire les valeurs représentant l'approximation polynomiale de chaque courbe de réponse relative à une position angulaire principale. 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 dans lequel l'étalonnage est obtenu à l'aide d'étalons circulaires, - la figure 2 est un diagramme montrant différentes courbes d'atténuation en fonction de la position des détecteurs ou canaux et du diamètre de l'étalon; - 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 de principe d'un- scanner à rayons X dans lequel l'étalonnage est obtenu selon des caractéristiques de l'invention; - la figure 5 est un diagramme montrant différentes courbes d'atténuation en fonction de la position du détecteur ou canal et de la position de la source à rayons X et des détecteurs par rapport à l'étalon unique et, - la figure 6 est un schéma fonctionnel d'un système de traitement des signaux de sortie des détecteurs. 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. Le schéma de la figure 4 est analogue à celui de la figure 1 en ce sens qu'il montre un scanner à rayons X qui comprend une source à rayons X 30 et un dispositif de détection 31 à plusieurs détecteurs ou canaux 35, entre lesquels est interposé, en fonctionnement normal, le corps d'un patient ou un étalon 32 de forme élliptique lors des opérations d'étalonnage. Bien que non représentée, une structure est prévue pour supporter la source 30 et le dispositif 31 et faire tourner l'ensemble autour du corps du patient ou de l'étalon, l'axe de rotation étant matérialisé par le point 33 situé dans l'étalon mais pas nécessairement en son centre géométrique. La source à rayons X 30 est commandée par un dispositif de commande 34 tandis que les différents détecteurs ou canaux 35 du dispositif 31 sont connectés à un système électronique 36 qui sera décrit en relation avec la figure 6. Le procédé d'étalonnage selon l'invention consiste à disposer la structure de support dans différentes positions angulaires dites principales, autour de l'étalon et à effectuer, à chaque position angulaire et de part et d'autre de cette position, une série de mesures de l'atténuation. Par exemple, on peut choisir seize positions angulaires principales régulièrement espacées de 22'30' et effectuer par exemple dix mesures élémentaires sur tous les canaux de détection autour de chaque position angulaire principale. A titre d'exemple, si le nombre de détecteurs est de 1024, le nombre total de vues sera de 16 x 10: 160 et le nombre total de mesures sera de 1024 x 160: 163840. A titre d'exemple, chaque vue autour d'une position angulaire principale est séparée d'un angle d'environ un tiers de degré. Ainsi donc, pour chaque position angulaire principale, il y a dix vues sur chacun des 1024 canaux, c'est-à-dire dix mesures par canal sur un secteur angulaire de trois degrés environ, et chaque mesure correspond à un trajet légèrement différent dans l'étalon et donc à une atténuation légèrement différente. Selon l'invention, dans chaque canal, ces dix mesures sont utilisées pour calculer une moyenne ce qui permet d'améliorer le rapport signal/bruit par un coefficient fîo. Ces valeurs moyennes au nombre de seize permettent de tracer, pour chaque position angulaire principale, une courbe de variation de l'atténuation en fonction de la position des canaux. En fait, ce nombre de courbes peut être ramené à quatre car les seize positions angulaires principales se réduisent à quatre si l'on tient compte des symétries entre les positions angulaires principales et celles de l'étalon. Ces courbes de variation sont montrées sur la figure 5 pour trois des quatre positions angulaires principales après prise en compte des symétries. Sur cette figure, les points 40, 41 ou 42 représentent, les valeurs moyennes d'atténuation pour certains des 1024 canaux, et la courbe 43 résulte d'une opération de lissage conformément à la présente invention, opération qui sera décrite ci-après en relation avec la figure 6. La courbe 43 correspond à la position angulaire de la figure 4, c'est-àdire que l'étalon elliptique 32 est irradié suivant le grand axe de l'ellipse. La courbe 44 correspond à la position angulaire représentée en pointillés sur la figure 4, c'est-à-dire que l'étalon 32 est irradié suivant le petit axe de l'ellipse. Enfin, la courbe 45 correspond à une position intermédiaire entre les deux précédentes. Par ailleurs, comme indiqué ci-dessus les seize positions angulaires principales sont symétriques deux à deux, les deux positions symétriques étant séparées deux à deux d'un angle de 180 . Ces deux positions symétriques conduisent à des courbes d'atténuation identiques et il est suffisant de n'en retenir qu'une seule. En fait, selon l'invention, il est proposé de combiner les mesures provenant de ces deux positions symétriques en en effectuant la moyenne de la manière qui sera expliquée plus en détail en relation avec la figure 6. Le rapport signal/bruit est alors amélioré par un facteur a/ soit 25 au total. En outre, si un étalon de forme symétrique par rapport à deux axes, par exemple orthogonaux, est utilisé, on peut grouper les mesures correspondant à des positions symétriques par rapport à ces axes. Il en résulte que les seize positions principales, qui ont été ramenées à huit par la première symétrie peuvent à nouveau être groupées deux par deux en en effectuant la moyenne, ce qui permet d'aboutir à quatre positions angulaires principales. Le rapport signal/bruit est à nouveau amélioré de /-, soit /40 au total. Les aspects principaux du procédé d'étalonnage selon l'invention sont: l'utilisation d'un étalon non circulaire; - la mesure des signaux de sortie des N canaux pour positions angulaires principales P de la structure source-détecteur et pour n positions élémentaires ou vues proches les unes des autres autour de chaque position angulaire principale; - le calcul de la moyenne des n mesures pour chaque canal et pour chacune des P positions angulaires principales; - le lissage des valeurs moyennes d'un canal au suivant pour chaque position angulaire principale de manière à obtenir une courbe de réponse dans laquelle les composantes haute fréquence ont été éliminées. Dans le cas de positions angulaires principales espacées régulièrement et en nombre pair il est prévu un calcul de la valeur moyenne des mesures correspondant & des positions principales espacées de 180 . En outre, si l'étalon a une forme symétrique par rapport à deux axes, le calcul de la valeur moyenne est effectué également sur les mesures obtenues pour des positions symétriques par rapport auxdits axes. Pour mettre en oeuvre le procédé, il est proposé de réaliser un système dont le schéma fonctionnel est donné par la figure 6. Sur cette figure, les signaux fournis par les 1024 détecteurs sont appliqués à un circuit de codage 50 qui fournit des valeurs ou codes numériques. Ce sont sur ces codes que seront effectuées toutes les opérations qui vont être décrites ci-après, ces opérations pouvant être réalisées par un calculateur convenablement programme. Les N codes résultant d'une mesure à une position élémentaire ou vue sont appliqués à un circuit de calcul de logarithme 51. Les valeurs logarithmiques sont soustraites dans un circuit 52 à une valeur logarithmique de référence REF fournie par un détecteur appelé moniteur qui reçoit directement le rayonnement X sans atténuation. Cette différence donne donc la valeur de l'atténuation introduite par l'étalon- par rapport à un trajet dans l'air. Ces N valeurs différentielles, correspondant chacune à un canal, sont ensuite soustraites dans un circuit 53 à N valeurs logarithmiques d'atténuation mesurées en l'absence d'étalon, c'est-à-dire selon des trajets dans l'air. Ces valeurs ont été mesurées et calculées préalablement à ce procédé d'étalonnage et sont utilisées ensuite, comme il est connu, lors des mesures sur le patient. Elles sont pour cette raison enregistrées dans une mémoire 54. Cette deuxième opération de soustraction permet de tenir compte d'une partie des disparités entre les canaux et donc de supprimer leur influence d'un canal au suivant. Les N valeurs résultant de cette deuxième soustraction sont enregistrées dans une mémoire 55 qui est prévue pour recevoir les N.P.n. codes qui résulteront des n mesures effectuées (ou n vues prises) à n positions élémentaires autour d'une position angulaire principale. Ce sont les n.N. codes qu'est effectué le calcul de la moyenne sur chaque canal par l'intermédiaire d'un circuit 56. Les codes des N valeurs moyennes sont enregistrés dans une mémoire 57. On comprend que la mémoire 55 et le circuit de calcul 56 peuvent être réalisés sous la forme d'un circuit additionneur qui réaliserait neuf opérations successives d'addition. Pour bénéficier des effets de la symétrie des P positions angulaires principales et de l'étalon, plusieurs manières de faire peuvent être mises en oeuvre. L'une d'entre elles consiste à effectuer successivement les mesures pour quatre positions angulaires principales symétriques et à effectuer la moyenne pour un canal, non pas sur dix valeurs mais sur quarante valeurs. Ceci implique une mémoire 55 prévue pour recevoir 40.N codes. Une autre manière consiste à utiliser P mémoires 57, une par position angulaire principale, et à effectuer un deuxième calcul de la moyenne.sur les valeurs contenues dans les quatre mémoires correspondant aux quatre positions angulaires principales symétriques. Après ces opérations de calcul de la moyenne qui ont pour but principal de réduire le bruit et donc d'améliorer le rapport signal/bruit, interviennent les opérations de lissage qui sont de deux ordres. D'abord une opération d'élimination des hautes fréquences dans des circuits 58, 59 et 60 et une opération d'approximation polynomiale dans un circuit 61. L'opération d'élimination des hautes fréquences comporte d'abord une transformée de Fourier dans le circuit 58 pour obtenir le spectre du signal sur les N canaux, puis un filtrage passe-bas dans le circuit 59 qui élimine les composantes haute fréquence, enfin une transformée de Fourier inverse de manière à revenir au signal initial mais sans composante haute fréquence. L'opération d'approximation polynomiale dans le circuit 61 peut être limitée à la partie centrale des courbes 43, 44 et 45 de la figure 5. Elle est effectuée par tous procédés et moyens connus. Les résultats de cette approximation polynomiale sont enregistrés dans une mémoire 62 qui peut contenir les valeurs numériques des P, P/2 ou P/4 courbes lissées selon que l'on tient compte ou non des symétries mentionnées ci-dessus, les courbes étant calculées l'une d'après l'autre. Les valeurs enregistrées dans la mémoire 62 à la fin des opérations pour P positions angulaires principales constituent des valeurs d'étalonnage qui sont utilisées de manière connue, lors des opérations de mesure sur le patient, pour calculer la valeur réelle de l'atténuation à partir de la valeur mesurée. Chaque courbe ainsi lissée de la mémoire 62, correspond à une courbe non lissée de la mémoire 57 et la différence entre les valeurs de ces deux courbes pour un même canal permet de calculer une valeur d'écart. Ainsi, lorsque pour ce même canal, la différence entre la valeur mesurée avec le patient et la valeur lissée est comprise dans cette valeur d'écart, on en déduit que la valeur réelle est celle de la courbe lissée. Lorsque cette différence dépasse cette valeur d'écart, il est nécessaire d'effectuer un calcul d'interpolation de type linéaire ou polynomial pour déterminer la valeur réelle comprise entre deux courbes lissées. Le système de mise en oeuvre du procédé a été décrit sous la forme de circuits fonctionnels mais il est clair qu'il peut être réalisé pour un calculateur qui réaliserait par programme les différentes opérations décrites. REVENDICATIONS 1 - Procédé d'étalonnage d'un scanner à rayons X qui comporte une source de rayonnement X et un dispositif de détection à N canaux, en utilisant un seul étalon de forme autre que circulaire, caractérisé en ce qu'il comprend les opérations suivantes: a) la mise en place de l'étalon entre la source de rayonnement X et le dispositif de détection; b) la mesure des atténuations dans les N canaux du dispositif de détection pour P positions angulaires principales de la structure tournante et pour n positions élémentaires ou vues proches les unes des autres autour de chaque position angulaire principale; c) le calcul de la moyenne des n atténuations pour chaque canal et pour chacune des P positions angulaires principales; d) le lissage des N valeurs moyennes des atténuations d'un canal au suivant pour chacune des P positions angulaires principales de manière à obtenir une courbe de réponse de l'atténuation en fonction de la position du canal dans laquelle les composantes haute fréquence ont été éliminées, ladite courbe de réponse ainsi obtenue étant la courbe réelle de variation de l'atténuation introduite par l'étalon en fonction de la position respective des N canaux. 2 - Procédé d'étalonnage selon la revendication 1, caractérisé en ce que l'opération de lissage comprend une opération de filtrage. 3 - Procédé d'étalonnage selon la revendication 2 caractérisé en ce que l'opération de lissage comprend en outre une opération d'approximation polynomiale qui est effectuée sur la courbe résultant du filtrage. 4 - Procédé d'étalonnage selon la revendication 2 ou 3, caractérisé en ce que l'opération de filtrage consiste à réaliser la transformée de Fourier des N valeurs moyennes correspondant à une des P positions angulaires principales, à filtrer le signal obtenu pour éliminer les composantes haute fréquence et à réaliser ensuite la transformée de Fourier inverse. - Procéde selon l'une quelconque des revendications précédentes, caractérisé en ce que les P positions angulaires principales sont en nombre pair, et sont choisies de manière à être symétriques deux à deux par rapport au centre de rotation et en ce que les valeurs des atténuations correspondant à des positions élémentaires symétriques sont moyennées. 6 - Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étalon a une forme qui est symétrique par rapport à deux axes et en ce que les valeurs des atténuations correspondant à des positions symétriques de l'étalon par rapport auxdits axes sont moyennées. 7 - Système pour mettre en oeuvre le procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend: - des premiers moyens pour calculer le logarithme des N signaux fournis par les N détecteurs, - des deuxièmes moyens pour calculer, pour chaque canal, l'atténuation introduite par l'étalon, - des troisièmes moyens pour calculer, pour chaque canal, la valeur moyenne des n atténuations mesurées pour une position angulaire principale * donnée, - des quatrièmes moyens pour calculer la transformée de Fourier des N atténuations moyennes correspondant à une position angulaire principale, - des cinquièmes moyens pour-éliminer les composantes haute fréquence dans le spectre résultant de la transformée de Fourier, - des sixièmes moyens pour calculer la transformée de Fourier inverse des composantes basse fréquence du spectre, - des sixièmes moyens pour calculer une approximation polynomiale de la courbe résultant de la transformée de Fourier inverse, - et des moyens pour garder en mémoire les valeurs représentant l'approximation polynomiale de chaque courbe de réponse relative à une position angulaire principale.