FR2664708A1 - Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. - Google Patents
Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. Download PDFInfo
- Publication number
- FR2664708A1 FR2664708A1 FR9008733A FR9008733A FR2664708A1 FR 2664708 A1 FR2664708 A1 FR 2664708A1 FR 9008733 A FR9008733 A FR 9008733A FR 9008733 A FR9008733 A FR 9008733A FR 2664708 A1 FR2664708 A1 FR 2664708A1
- Authority
- FR
- France
- Prior art keywords
- sep
- butterfly filter
- channel
- signals
- coefficients
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05G—X-RAY TECHNIQUE
- H05G1/00—X-ray apparatus involving X-ray tubes; Circuits therefor
- H05G1/08—Electrical details
- H05G1/26—Measuring, controlling or protecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
L'invention concerne, dans les scanners à rayons X, un procédé pour corriger les effets dus au durcissement du faisceau de rayons X polychromatique. La correction qui est effectuée est du type polynomial dont les coefficients alpha2 à alphan sont déterminés par le procédé de l'invention. Ce procédé consiste à mesurer les signaux des N canaux 12 pour n fantômes circulaires de rayons R1 à Rn et à calculer les rapports des intégrales des signaux pour les fantômes pris deux à deux et les comparer aux rapports correspondants des carrés des rayons. Ce système d'équations permet de déterminer alpha2 à alphan .
Description
PROCEDE DE CORRECTION DU DURCISSEMENT
D'UN FAISCEAU DE RAYONS X DANS UN SCANNER
L'invention concerne les scanners à rayons X et, plus particulièrement dans de tels scanners, un procédé pour corriger les effets dus au durcissement du faisceau de rayons X polychromatique lorsque ledit faisceau traverse un corps d'épaisseur variable.
D'UN FAISCEAU DE RAYONS X DANS UN SCANNER
L'invention concerne les scanners à rayons X et, plus particulièrement dans de tels scanners, un procédé pour corriger les effets dus au durcissement du faisceau de rayons X polychromatique lorsque ledit faisceau traverse un corps d'épaisseur variable.
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 d des rayons X dans le corps selon la formule
I = IOe~d (1) formule dans laquelle
Io est l'intensité du rayonnement entrant dans le
corps humain,
I est l'intensité du rayonnement sortant du corps
humain,
est un coefficient d'atténuation qui dépend du
corps traversé.
I = IOe~d (1) formule dans laquelle
Io est l'intensité du rayonnement entrant dans le
corps humain,
I est l'intensité du rayonnement sortant du corps
humain,
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 à d, c'est-à-dire qu'elle est proportionnelle au trajet d dans le corps.
I/Io est égale à d, c'est-à-dire qu'elle est proportionnelle au trajet d dans le corps.
Ces appareils sont constitués essentiellement, comme le montre la figure 1, d'une source 10 de rayons X 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 un corps 15 à 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 15.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 15 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 15 et on mesure les signaux de sortie des détecteurs élémentaires 12 pour les traiter de manière appropriée dans un dispositif électronique 14 selon des procédés connus afin d'en tirer une image représentative de la section transversale.
Lorsque le faisceau de rayons X est monochromatique, l'atténuation est directement proportionnelle au trajet d du rayon dans le corps 15, ce qui correspond à un coefficient A fixe. Cependant, le faisceau de rayons X est polychromatique, ce qui a pour conséquence que l'atténuation varie aussi avec la longueur d'onde du faisceau. Ainsi les rayons de faible énergie sont très atténués, jusqu'à être complètement absorbés par le corps 15, tandis que les rayons d'énergie plus forte sont moins atténués, traversent le corps 15 pour être détectés par le dispositif de détection 11. Ceci a pour effet que l'énergie moyenne du faisceau de rayons X à la sortie du corps 15 est plus grande qu'à l'entrée, phénomène que l'on appelle durcissement du faisceau polychromatique.
Ce phénomène de durcissement du faisceau polychromatique se traduit dans la formule (1) par une baisse apparente du coefficient , coefficient qui, dans le cas d'un faisceau polychromatique, correspond à une valeur moyenne. Il en résulte une atténuation moins importante des rayons de haute énergie et donc une augmentation de la valeur des signaux détectés par les détecteurs élémentaires 12 par rapport à la valeur qu'ils auraient en l'absence de ce phénomène.
Par ailleurs, ce phénomène de durcissement est d'autant plus accentué que le trajet d dans le corps 15 est plus grand. I1 y a donc aussi une baisse apparente plus importante du coefficient y pour un point se trouvant au centre 0 du corps 15 que pour un point du bord du même corps. Sur l'image de la section d'un corps homogène qui est obtenue après reconstruction, il y a baisse progressive de la densité vers le centre.
Pour tenir compte de ce phénomène, il a été proposé d'effectuer une correction des signaux détectés par chaque détecteur élémentaire. La correction la plus répandue est du type polynomial tel que
formule dans laquelle Xj est le signal d'atténuation du canal de rang j parmi les n canaux ou détecteurs élémentaires 11; Xj est donné par
Xj = Kj log + ssj (3) Io avec Kj et Pj des valeurs individuelles du canal j de rang j après étalonnage.
formule dans laquelle Xj est le signal d'atténuation du canal de rang j parmi les n canaux ou détecteurs élémentaires 11; Xj est donné par
Xj = Kj log + ssj (3) Io avec Kj et Pj des valeurs individuelles du canal j de rang j après étalonnage.
Les coefficients ai de la formule (2) sont choisis pour que l'image d'un corps homogène ne présente pas de variation de densité.
En l'absence de durcissement du faisceau, on devrait avoir ydj = Xj; du fait de la correction polynomiale indiquée ci-dessus, on a :
dj = K'Ycor + #n (4) en étant l'erreur provenant d'un développement d'ordre n de la formule (2) et K' un coefficient d'adaptation.
dj = K'Ycor + #n (4) en étant l'erreur provenant d'un développement d'ordre n de la formule (2) et K' un coefficient d'adaptation.
On peut trouver un jeu de coefficients polynomiaux ai de sorte que En soit négligeable pour un trajet (ou épaisseur) dj compris entre zéro et une valeur correspondant à un patient de corpulence moyenne.
Si #n est négligeable, on peut écrire :
dj = K' Ycor (5) et 'dj = Ycor (6) avec ' = /K' (7)
En utilisant les équations (2) et (6), on obtient :
dj = K' Ycor (5) et 'dj = Ycor (6) avec ' = /K' (7)
En utilisant les équations (2) et (6), on obtient :
Sur la figure 2, la droite 16 représente la variation du signal d'atténuation X en fonction de 'd dans le cas d'un faisceau monochromatique, c'est-à-dire ai= , et la courbe 17 représente X en fonction de 'd dans le cas d'un faisceau polychromatique.Ces courbes montrent que pour une épaisseur d et un coefficient ' donnés, le signal mesuré XA est inférieur au signal théorique XB que l'on aurait obtenu avec un faisceau monochromatique, ce qui correspond bien à une diminution apparente du coefficient '.
Dans les scanners, il est connu d'utiliser un filtre-papillon d'épaisseur variable d' qui est placé à la sortie de la source. Ce filtre est de faible épaisseur au centre du faisceau qui correspond au trajet le plus long dans le corps et de forte épaisseur sur les bords du faisceau qui correspondent à des trajets plus courts dans le corps. Ce filtre a pour effets bénéfiques de réduire la dose de rayonnement, notamment sur les côtés, à ce qui est nécessaire et à uniformiser les signaux au niveau des détecteurs, ce qui permet une amplification plus grande des signaux. Cependant, il a aussi pour effet de durcir le faisceau de rayons X qui irradie le patient et de changer son spectre d'énergie d'un canal à l'autre.
I1 en résulte qu'il faut en tenir compte pour déterminer les coefficients ai du polynôme de correction et trouver pour chaque canal un jeu de coefficients ai de sorte que 1 'épaisseur du filtre-papillon "vue" par le canal soit prise en compte et que la densité d'un fantôme homogène soit constante.
Le but de la présente invention est de mettre en oeuvre un procédé de correction du durcissement du faisceau de rayons X d'un scanner et, plus particulièrement, un procédé pour déterminer les coefficients du polynôme de correction qui tient compte de la présence d'un filtre-papillon.
La présente invention concerne un procédé de correction du durcissement du faisceau de rayons X dans un scanner comportant un filtre-papillon d'épaisseur variable d' qui consiste à effectuer une correction de type polynomial tel que le signal corrigé Ycor sera de la forme
avec Xj le signal du canal de rang j et Qi le coefficient de correction polynomial d'ordre i, caractérisé par un changement de coordonnées pour passer d'un système d'axes X1 01 Y1 correspondant à une courbe d'atténuation du signal d'un canal avec filtre-papillon d'épaisseur variable d' à un système d'axes X2 02 Y2 correspondant à une courbe d'atténuation du signal d'un canal sans filtre-papillon tel que le point 02 a pour coordonnées XO et YO telles que
X2 = X1 - Xo (9)
Y2 = Y1 - Yo (10) avec Yo = IL'd' (11) et XO la racine réelle et positive de l'équation
et caractérisé par le fait que la correction qui sera effectuée pour un canal donné qui voit le chemin optique IL'd' a travers le filtre-papillon sera donnée par
avec Xj le signal du canal de rang j et Qi le coefficient de correction polynomial d'ordre i, caractérisé par un changement de coordonnées pour passer d'un système d'axes X1 01 Y1 correspondant à une courbe d'atténuation du signal d'un canal avec filtre-papillon d'épaisseur variable d' à un système d'axes X2 02 Y2 correspondant à une courbe d'atténuation du signal d'un canal sans filtre-papillon tel que le point 02 a pour coordonnées XO et YO telles que
X2 = X1 - Xo (9)
Y2 = Y1 - Yo (10) avec Yo = IL'd' (11) et XO la racine réelle et positive de l'équation
et caractérisé par le fait que la correction qui sera effectuée pour un canal donné qui voit le chemin optique IL'd' a travers le filtre-papillon sera donnée par
Les (n-l) coefficients ai et le coefficient ' sont obtenus par les opérations suivantes (a) enlèvement du filtre-papillon; (b) mise en place successive de n fantômes cylindriques
de rayons R1 à Rn et mesure pour chaque fantôme des
signaux Xj des canaux; (c) calcul pour chaque fantôme de l'intégrale S1 à Sn
des signaux des N canaux; (d) calcul des rapports
de manière à obtenir un système de (n-l) équations dont les (n-1) inconnues sont les coefficients a2 à an; (e) calcul des coefficients a2 à an; (f) mise en place du filtre-papillon et d'un fantôme de
faible épaisseur d1 et mesure des signaux
des N canaux;; (g) calcul de ' par la formule
Y d1 (6')
D'autres buts, caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description suivante d'un exemple particulier de mise en oeuvre du procédé, ladite description étant faite en relation avec les dessins joints dans lesquels - la figure 1 est une vue schématique d'un scanner à
rayons X auquel s'applique le procédé selon
l'invention, - la figure 2 est un diagramme de courbes montrant
l'effet du durcissement du faisceau sur la mesure du
signal d'atténuation, et, - la figure 3 est un diagramme montrant le changement de
coordonnées qui est effectué dans le procédé selon
l'invention.
de rayons R1 à Rn et mesure pour chaque fantôme des
signaux Xj des canaux; (c) calcul pour chaque fantôme de l'intégrale S1 à Sn
des signaux des N canaux; (d) calcul des rapports
de manière à obtenir un système de (n-l) équations dont les (n-1) inconnues sont les coefficients a2 à an; (e) calcul des coefficients a2 à an; (f) mise en place du filtre-papillon et d'un fantôme de
faible épaisseur d1 et mesure des signaux
des N canaux;; (g) calcul de ' par la formule
Y d1 (6')
D'autres buts, caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description suivante d'un exemple particulier de mise en oeuvre du procédé, ladite description étant faite en relation avec les dessins joints dans lesquels - la figure 1 est une vue schématique d'un scanner à
rayons X auquel s'applique le procédé selon
l'invention, - la figure 2 est un diagramme de courbes montrant
l'effet du durcissement du faisceau sur la mesure du
signal d'atténuation, et, - la figure 3 est un diagramme montrant le changement de
coordonnées qui est effectué dans le procédé selon
l'invention.
Comme on l'a exposé dans le préambule, la plupart des scanners mobiles sont pourvus d'un filtre-papillon qui change le spectre d'énergie du faisceau d'un canal à l'autre. Cela signifie qu'il faut trouver pour chaque canal un jeu de coefficients de correction ai de sorte que l'épaisseur du filtre-papillon "vue" par lui soit prise en compte et que la densité d'un fantôme homogène soit constante.
La correction doit être telle qu'un canal, qui a été étalonné en présence d'un filtre-papillon, doit donner, par définition, un signal d'atténuation nul quand il ne voit que le filtre-papillon. Aussi, la correction proposée correspond à un changement de coordonnées de manière à passer d'un premier système réel où le corps à examiner est balayé à travers un filtre-papillon à un deuxième système où le filtre-papillon est apparemment absent, mais chaque canal a une autre courbe de correction en fonction du point de coupure du filtre-papillon.
Ce changement de coordonnées est montré par la figure 3 pour la courbe 18 d'un canal qui voit son chemin optique IL'd' à travers le filtre-papillon dans les deux systèmes d'axes X1 O1 Y1 et X2 O2 Y2, Le point O2 ayant les coordonnées XO et YO tel que
X2 = X1 - Xo (9)
Y2 = Y1 - Yo (10) avec Yo = 'd' (11) et Xo qui est donné par la racine réelle et positive de l'équation (8) tel que
IL'd' étant le chemin optique vu par un canal à travers le filtre-papillon.
X2 = X1 - Xo (9)
Y2 = Y1 - Yo (10) avec Yo = 'd' (11) et Xo qui est donné par la racine réelle et positive de l'équation (8) tel que
IL'd' étant le chemin optique vu par un canal à travers le filtre-papillon.
Des équations (9), (10) et (11), on tire :
X1 = X2 + Xo (13)
Y1 = Y2 + 'd' (14) puis en reportant ces valeurs dans l'équation (8), on obtient
X1 = X2 + Xo (13)
Y1 = Y2 + 'd' (14) puis en reportant ces valeurs dans l'équation (8), on obtient
La correction qui sera effectuée pour un canal donné qui voit un chemin optique IL'd' à travers le filtre-papillon sera
formule dans laquelle une partie des termes s'annulera compte-tenu de l'équation (12). La correction définie par l'équation (16) dépend de n paramètres - (n-l) paramètres pour ai - 1 paramètre pour '
Par ailleurs, d' dépend uniquement de la formule du filtre-papillon et XO dépend (formule 12) de IL'd' et ai.
formule dans laquelle une partie des termes s'annulera compte-tenu de l'équation (12). La correction définie par l'équation (16) dépend de n paramètres - (n-l) paramètres pour ai - 1 paramètre pour '
Par ailleurs, d' dépend uniquement de la formule du filtre-papillon et XO dépend (formule 12) de IL'd' et ai.
Afin de déterminer ces n paramètres, l'invention propose un procédé basé sur les remarques suivantes.
Dans le cas de deux corps homogènes de forme cylindrique dont les sections ont des rayons R1 et R2 , l'intégrale des signaux d'atténuation de tous les N canaux sont tels que
S étant égal à
à condition que le rayonnement X soit monochromatique et que les corps cylindriques soient centrés autour de l'isocentre O du scanner.
à condition que le rayonnement X soit monochromatique et que les corps cylindriques soient centrés autour de l'isocentre O du scanner.
Comme ces coefficients ai sont identiques en présence ou non du filtre-papillon, on détermine ces paramètres sans filtre-papillon.
Pour déterminer ces paramètres, le procédé selon l'invention consiste à utiliser n fantômes cylindriques différents de rayons respectifs R1 à Rn et à mesurer les signaux Xj des n canaux pour chacun des n fantômes.
Afin d'améliorer le rapport signal/bruit, on peut également prendre la moyenne des signaux Xj pour un tour du scanner.
L'étape suivante consiste à mettre en oeuvre l'équation (17) en calculant les (n-l) rapports suivants
<SEP> - <SEP> N <SEP> Xj <SEP> i <SEP>
<tb> <SEP> Xi <SEP> + <SEP> i <SEP> Xj <SEP> 2
<tb> <SEP> sl <SEP> g <SEP> i=2 <SEP> fantôme <SEP> 1 <SEP> R1
<tb> <SEP> = <SEP> =, <SEP> (l82)
<tb> <SEP> N <SEP> Xi <SEP> t <SEP> - <SEP> 2 <SEP> 2
<tb> <SEP> z <SEP> X;<SEP> aI
<tb> <SEP> j <SEP> Xj <SEP> 1=2 <SEP> fantôme <SEP> 2 <SEP> R2
<tb> <SEP> N <SEP> i
<tb> <SEP> z <SEP> Xi <SEP> z <SEP> ai <SEP> Xj <SEP> 2
<tb> <SEP> 3 <SEP> i <SEP> 3 <SEP> i=2 <SEP> fantôme <SEP> 2 <SEP> R2
<tb> = <SEP> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ <SEP> = <SEP> ~~~ <SEP> (183)
<tb> <SEP> N <SEP> 2
<tb> <SEP> z <SEP> Xj <SEP> xj <SEP> fantôme <SEP> 2
<tb> <SEP> S3 <SEP> i <SEP> Xj <SEP> i=2 <SEP> fantôme <SEP> 3
<tb> et ainsi de suite jusqu'au nième fantôme
<tb> <SEP> Xi <SEP> + <SEP> i <SEP> Xj <SEP> 2
<tb> <SEP> sl <SEP> g <SEP> i=2 <SEP> fantôme <SEP> 1 <SEP> R1
<tb> <SEP> = <SEP> =, <SEP> (l82)
<tb> <SEP> N <SEP> Xi <SEP> t <SEP> - <SEP> 2 <SEP> 2
<tb> <SEP> z <SEP> X;<SEP> aI
<tb> <SEP> j <SEP> Xj <SEP> 1=2 <SEP> fantôme <SEP> 2 <SEP> R2
<tb> <SEP> N <SEP> i
<tb> <SEP> z <SEP> Xi <SEP> z <SEP> ai <SEP> Xj <SEP> 2
<tb> <SEP> 3 <SEP> i <SEP> 3 <SEP> i=2 <SEP> fantôme <SEP> 2 <SEP> R2
<tb> = <SEP> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ <SEP> = <SEP> ~~~ <SEP> (183)
<tb> <SEP> N <SEP> 2
<tb> <SEP> z <SEP> Xj <SEP> xj <SEP> fantôme <SEP> 2
<tb> <SEP> S3 <SEP> i <SEP> Xj <SEP> i=2 <SEP> fantôme <SEP> 3
<tb> et ainsi de suite jusqu'au nième fantôme
<tb> 5n-1 <SEP> S.<SEP> Xj <SEP> N <SEP> Xj <SEP> xj <SEP> fantôme <SEP> n-1 <SEP> Rn-l
<tb> 2
<tb> <SEP> 3
<tb> <SEP> = <SEP> =
<tb> 5n <SEP> N <SEP> +.N <SEP> j <SEP> + <SEP> 2 <SEP> ti <SEP> j <SEP> ] <SEP> fantôme <SEP> n <SEP> - <SEP> 2 <SEP> Rn
<tb> <SEP> 3
<tb>
Ce système de (n-l) équations (182) à (18n) permet de déterminer les (n-l) coefficients α2 à αn.
<tb> 2
<tb> <SEP> 3
<tb> <SEP> = <SEP> =
<tb> 5n <SEP> N <SEP> +.N <SEP> j <SEP> + <SEP> 2 <SEP> ti <SEP> j <SEP> ] <SEP> fantôme <SEP> n <SEP> - <SEP> 2 <SEP> Rn
<tb> <SEP> 3
<tb>
Ce système de (n-l) équations (182) à (18n) permet de déterminer les (n-l) coefficients α2 à αn.
Le coefficient ' peut être déterminé en mettant en place le filtre-papillon et en mesurant les signaux Xj des canaux avec un fantôme de faible épaisseur dl tel que Ycor est sensiblement égal à Y.
Alors ' est donné par application de la formule (6)
' = Y (6')
a1
Le procédé de correction du durcissement d'un faisceau de rayons X dans un scanner, c'est-à-dire le procédé qui permet de déterminer les coefficients a2 à an de la correction polynomiale de degré n et le coefficient ', comprend les étapes suivantes (a) enlèvement du filtre-papillon; (b) mise en place successive de n fantômes cylindriques
de rayons R1 à Rn et mesure pour chaque fantôme des
signaux Xj des canaux; (c) calcul pour chaque fantôme de l'intégrale S1 à Sn
des signaux des N canaux.
' = Y (6')
a1
Le procédé de correction du durcissement d'un faisceau de rayons X dans un scanner, c'est-à-dire le procédé qui permet de déterminer les coefficients a2 à an de la correction polynomiale de degré n et le coefficient ', comprend les étapes suivantes (a) enlèvement du filtre-papillon; (b) mise en place successive de n fantômes cylindriques
de rayons R1 à Rn et mesure pour chaque fantôme des
signaux Xj des canaux; (c) calcul pour chaque fantôme de l'intégrale S1 à Sn
des signaux des N canaux.
(d) calcul des rapports
de manière à obtenir un système de (n-l) équations dont les (n-l) inconnues sont les coefficients a2 à an; (e) calcul des coefficients a2 à an (f) mise en place du filtre-papillon et d'un fantôme de
faible épaisseur dl et mesure des signaux
des N canaux; (g) calcul de ' par la formule
' = Y (6')
d1 I1 est possible de mettre en oeuvre le procédé en maintenant en place le filtre-papillon mais les calculs des coefficients ai sont plus compliqués car il faut tenir compte de IL'd'.
de manière à obtenir un système de (n-l) équations dont les (n-l) inconnues sont les coefficients a2 à an; (e) calcul des coefficients a2 à an (f) mise en place du filtre-papillon et d'un fantôme de
faible épaisseur dl et mesure des signaux
des N canaux; (g) calcul de ' par la formule
' = Y (6')
d1 I1 est possible de mettre en oeuvre le procédé en maintenant en place le filtre-papillon mais les calculs des coefficients ai sont plus compliqués car il faut tenir compte de IL'd'.
Claims (2)
1. Procédé de correction du durcissement d'un faisceau de rayons X dans un scanner comportant un filtre-papillon d'épaisseur variable d' qui consiste à effectuer une correction de type polynomial tel que le signal corrigé Ycor sera de la forme
avec Xj le signal du canal de rang j et ai le coefficient de correction polynomial d'ordre i, caractérisé par un changement de coordonnées pour passer d'un système d'axes X1 01 Y1 correspondant à une courbe d'atténuation du signal d'un canal avec filtre-papillon d'épaisseur variable d' à un système d'axes X2 02 Y2 correspondant à une courbe d'atténuation du signal d'un canal sans filtre-papillon tel que le point 02 a pour coordonnées XO et YO telles que
X2 = X1 - Xo (9)
Y2 = Y1 - Yo (10) avec YO = IL'd' (11) et XO la racine réelle et positive de l'équation
et caractérisé par le fait que la correction qui sera effectuée pour un canal donné qui voit le chemin optique IL'd' à travers le filtre-papillon sera donnée par
2. Procédé de correction selon la revendication 1, caractérisé en ce que les (n-l) coefficients ai et le coefficient ' sont obtenus en effectuant les opérations suivantes (a) enlèvement du filtre-papillon; (b) mise en place successive de n fantômes cylindriques
de rayons R1 à Rn et mesure pour chaque fantôme des
signaux Xj des canaux; (c) calcul pour chaque fantôme de l'intégrale S1 à Sn
des signaux des N canaux.
(d) calcul des rapports
d1
Y ' = Y (6')
des N canaux; (g) calcul de ' par la formule
faible épaisseur dl et mesure des signaux
a2 à an; (e) calcul des coefficients a2 à an (f) mise en place du filtre-papillon et d'un fantôme de
dont les (n-l) inconnues sont les coefficients
de manière à obtenir un système de (n-l) équations
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9008733A FR2664708B1 (fr) | 1990-07-10 | 1990-07-10 | Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9008733A FR2664708B1 (fr) | 1990-07-10 | 1990-07-10 | Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2664708A1 true FR2664708A1 (fr) | 1992-01-17 |
| FR2664708B1 FR2664708B1 (fr) | 1992-09-25 |
Family
ID=9398533
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR9008733A Expired - Fee Related FR2664708B1 (fr) | 1990-07-10 | 1990-07-10 | Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. |
Country Status (1)
| Country | Link |
|---|---|
| FR (1) | FR2664708B1 (fr) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1296011C (zh) * | 2003-07-16 | 2007-01-24 | 东软飞利浦医疗设备系统有限责任公司 | 一种ct机射束硬化的校正方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3755672A (en) * | 1970-11-30 | 1973-08-28 | Medinova Ab So | Exposure compensating device for radiographic apparatus |
| JPS6157840A (ja) * | 1984-08-29 | 1986-03-24 | Toshiba Corp | 放射線断層検査装置 |
| JPS62127657A (ja) * | 1985-11-29 | 1987-06-09 | Yokogawa Medical Syst Ltd | Ｘ線ｃｔ装置の１次フイルタ |
-
1990
- 1990-07-10 FR FR9008733A patent/FR2664708B1/fr not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3755672A (en) * | 1970-11-30 | 1973-08-28 | Medinova Ab So | Exposure compensating device for radiographic apparatus |
| JPS6157840A (ja) * | 1984-08-29 | 1986-03-24 | Toshiba Corp | 放射線断層検査装置 |
| JPS62127657A (ja) * | 1985-11-29 | 1987-06-09 | Yokogawa Medical Syst Ltd | Ｘ線ｃｔ装置の１次フイルタ |
Non-Patent Citations (3)
| Title |
|---|
| PATENT ABSTRACTS OF JAPAN, vol. 10, no. 220 (P-482)[2276], 31 juillet 1986; & JP-A-61 57 840 (TOSHIBA CORP.) 24-03-1986 * |
| PATENT ABSTRACTS OF JAPAN, vol. 11, no. 345 (P-636)[2792], 12 novembre 1987; & JP-A-62 127 657 (YOKOGAWA MEDICAL SYST. LTD) 09-06-1987 * |
| RADIOLOGY, vol. 99, no. 3, juin 1971, pages 694-696; P.R. EDHOLM et al.: "Primary X-ray dodging" * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1296011C (zh) * | 2003-07-16 | 2007-01-24 | 东软飞利浦医疗设备系统有限责任公司 | 一种ct机射束硬化的校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| FR2664708B1 (fr) | 1992-09-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0346181B1 (fr) | Procédé et système d'étalonnage d'un scanner à rayons x en utilisant un ou plusieurs étalons circulaires excentrés | |
| EP0409698A1 (fr) | Procédé de correction de la mesure de la densité osseuse dans un scannner | |
| EP1852717A1 (fr) | Procédé d'estimation de rayonnement diffusé dans un détecteur bidimensionnel | |
| FR2705785A1 (fr) | Procédé pour déterminer la fonction d'atténuation d'un objet par rapport à la transmission d'une épaisseur de référence d'un matériau de référence et dispositif pour la mise en Óoeuvre du procédé. | |
| EP0360653A1 (fr) | Procédé et système de correction des défauts d'images d'un scanner dus aux déplacements de ce dernier | |
| EP0871143A1 (fr) | Système de traitement d'une séquence d'images bruitées, et appareil d'examen médical incluant un tel système | |
| EP0334762B1 (fr) | Procédé et système d'étalonnage d'un scanner à rayons X en utilisant un seul étalon non circulaire | |
| Fusco et al. | Myopic deconvolution method for adaptive optics images of stellar fields | |
| FR2484824A1 (fr) | Dispositif pour la visualisation de couches du corps au moyen d'un rayonnement monoenergetique | |
| EP0488889A1 (fr) | Procédé et dispositif de reconstruction d'images tridimensionnelles d'un objet en utilisant deux trajectoires circulaires d'acquisition | |
| FR2519772A1 (fr) | Dispositif de detection de reference pour tomodensitometre multidetecteur et tomodensitometre comportant un tel dispositif | |
| FR2664708A1 (fr) | Procede de correction du durcissement d'un faisceau de rayons x dans un scanner. | |
| EP0197582B1 (fr) | Procédé et appareil d'exploration de milieux par échographie ultrasonore | |
| EP0800655B1 (fr) | Methode et dispositif pour determiner au moins un parametre caracteristique d'un corps | |
| EP3153888A1 (fr) | Procédé de correction d'un spectre | |
| FR2656697A1 (fr) | Procede de mesure de la densite osseuse a l'aide d'un scanner. | |
| FR2507466A1 (fr) | Procede et appareil de tomographie a faisceau en eventail continument emis utilisant un filtre d'optimisation | |
| EP3635377B1 (fr) | Procédé de détermination de propriétés physiques d'un échantillon | |
| EP2981213A1 (fr) | Procede et systeme de caracterisation d'un tissu osseux | |
| FR2700259A1 (fr) | Dispositif automatique de collimation d'un faisceau de rayons X d'un scanner et procédé de mise en Óoeuvre. | |
| EP1058128B1 (fr) | Dispositif et procédé de discrimination d'impulsions de détecteurs de rayonnement à semi-conducteur | |
| FR2780183A1 (fr) | Procede d'amelioration du rapport signal/bruit de l'image d'un objet en mouvement | |
| FR2838852A1 (fr) | Procede de correction du rayonnement diffuse de photons x dans un dispositif de radiographie comportant des lames de collimation | |
| WO2013014132A1 (fr) | Dispositif et procede de traitement pour une mesure spectrometrique d'un flux de photons | |
| JP2826297B2 (ja) | 光ｃｔの再構成方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |