FR2885007A1 - Procede de compensation de la diffusion et procede de mesure - Google Patents

Procede de compensation de la diffusion et procede de mesure Download PDF

Info

Publication number
FR2885007A1
FR2885007A1 FR0601817A FR0601817A FR2885007A1 FR 2885007 A1 FR2885007 A1 FR 2885007A1 FR 0601817 A FR0601817 A FR 0601817A FR 0601817 A FR0601817 A FR 0601817A FR 2885007 A1 FR2885007 A1 FR 2885007A1
Authority
FR
France
Prior art keywords
pro
dispersion
subject
data
projection length
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.)
Withdrawn
Application number
FR0601817A
Other languages
English (en)
Inventor
Masatake Nukui
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of FR2885007A1 publication Critical patent/FR2885007A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • 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
    • G01N23/046—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 using tomography, e.g. computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/58—Testing, adjusting or calibrating thereof
    • A61B6/582—Calibration
    • A61B6/583—Calibration using calibration phantoms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/419—Imaging computed tomograph
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/612—Specific applications or type of materials biological material
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S378/00—X-ray or gamma ray systems or devices
    • Y10S378/901—Computer tomography program or processor

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

La présente invention vise à compenser la disposition qui survient pendant une radiographie multicoupe. La longueur de projection pro_ex d'entités non-sujets, à savoir des entités autres qu'un sujet, provoquant la dispersion est acquise, et la longueur de projection pro_pt d'un objet à radiographier est acquise. L'objet à radiographier est radiographié avec une épaisseur de faisceau établie à la même valeur qu'une épaisseur do de détecteur afin de mesurer des données I (do, do). L'objet à radiographier est radiographié avec une épaisseur de faisceau établie à la même valeur qu'une épaisseur do de détecteur afin de mesurer des données I (do, d). Une ampleur de dispersion S (do, d) est calculée d'après la différence entre les données I (do, do) et des données I (do, d), et est enregistrée en association avec la somme des logueurs de projection pro_pt + pro_ex. Après la radiographie d'un sujet por produire des données, la longueur de projection pro_ex des entités non sujets ayant affecté les données est calculée. L'amplitude de dispersion S (do, d) associée à la somme de la longueur de projection pro_pt et de la longueur de projection pro_ex est extraite. Ainsi, les données sont corrigées afin de compenser la dispersion. Finalement, il est possible de supprimer la balise de qualité d'image tomodensiométrique attribuable à la dispersion.

Description

PROCEDE DE COMPENSATION DE DISPERSION ET PROCEDE DE
MESURE
La présente invention est relative à un procédé de compensation de dispersion, un procédé de mesure de dispersion et un système de tomographie. Plus particulièrement, la présente invention concerne un procédé de compensation de dispersion, un procédé de mesure de dispersion et un système de tomographie pour compenser l'effet négatif de la dispersion survenant pendant une radiographie multicoupe.
On connaît dans la technique antérieure divers procédés de compensation de dispersion adaptables à des systèmes de tomographie employant un détecteur de rayons X à un seul canal (cf., par exemple, le document de brevet n 1 et le document de brevet n 2).
De plus, un procédé de compensation de dispersion adaptable à des 15 systèmes de tomographie employant un détecteur de rayons X à plusieurs canaux a été proposé (cf., par exemple, document de brevet n 3).
[Document de brevet n 1] Publication non examinée de brevet japonais n 7-213 517 [Document de brevet n 2] Publication non examinée de brevet japonais n 20 8-131 431 [Document de brevet n 3] Publication non examinée de brevet japonais n 11-299 768 Lorsqu'on utilise un détecteur de rayons X à plusieurs canaux ayant une pluralité de groupes de détecteurs pour effectuer une radiographie multicoupe, l'épaisseur d'un faisceau est supérieure à l'épaisseur de chaque détecteur (ou de chaque groupe de détecteurs). Ainsi, l'effet préjudiciable de la dispersion s'aggrave.
Cependant, un procédé de compensation de dispersion selon la technique antérieure pour des systèmes tomographiques employant un détecteur de rayons X à un seul canal ne tient pas compte du cas évoqué ci-dessus et ne peut pas être appliqué à ce cas.
En outre, un procédé de compensation de dispersion selon la technique antérieure pour des systèmes tomographiques employant un détecteur de rayons X à plusieurs canaux peut être appliqué au cas évoqué ci-dessus, mais nécessite deux reconstructions d'image. Ainsi, ce procédé impose un gros travail de traitement.
La présente invention vise à réaliser un procédé de compensation de dispersion pour compenser l'effet préjudiciable de dispersion survenant pendant une radiographie multicoupe, un procédé de mesure de dispersion pour mesurer une ampleur de dispersion et un système de tomographie.
Selon le premier aspect de la présente invention, il est proposé un procédé de compensation de dispersion comprenant les étapes consistant à : acquérir une longueur de projection pro_ex d'entités non-sujets, c'està-dire des entités autres qu'un sujet, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée lorsqu'elle survient dans un passage de rayons X au cours de la radiographie d'un sujet; acquérir une longueur de projection pro_pt d'un objet à radiographier; mesurer des données I(do,do) en effectuant une radiographie avec les entités non-sujets et l'objet à radiographier introduits dans le passage de rayons X et avec une épaisseur de faisceau établie à la même valeur qu'une épaisseur do de détecteur; mesurer des données I(do,d) en radiographiant l'objet à radiographier avec l'épaisseur de faisceau établie à une valeur d supérieure à l'épaisseur do de détecteur; calculer une ampleur de dispersion S(do,d) d'après une différence entre les données I(do,do) et les données I(do,d) ; stocker les informations nécessaires au calcul de l'ampleur de la dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier; radiographier l'objet de manière à produire des données, et utiliser les données pour calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets qui ont affecté les données; et lire les informations correspondant à la longueur de projection pro_pt du sujet et à la longueur de projection pro_ex des entités non-sujets, et utiliser les informations pour corriger les données de manière à compenser la dispersion.
Dans le procédé ci-dessus, l'objet à radiographier est normalement virtuel, mais il peut s'agir d'un sujet. De plus, les paramètres écrits entre parenthèses pour distinguer les données I(do,do), les données I(do, d) ou l'ampleur de dispersion S(do,d) désignent respectivement une épaisseur de détecteur et une épaisseur de faisceau.
Dans le procédé de compensation de dispersion selon le premier aspect, l'épaisseur do du détecteur est maintenue inchangée, seule l'épaisseur d du faisceau est modifiée. On estime par conséquent qu'un incrément des données I(do,do) aux données I(do,d) est imputable à la dispersion. Par conséquent, l'ampleur de dispersion S(do,d) peut être calculée d'après la différence entre les données I(do,do) et les données I(do,d). Par ailleurs, on ne tient pas compte de la dispersion représentée par les données I(do,do). Autrement, puisque les données I(do, do) sont considérées comme des données produites par le détecteur de rayons X à un seul canal, la dispersion représentée par les données I(do, do) peut être compensée suivant le procédé de compensation de dispersion classique.
La dispersion de rayons X survenant pendant la radiographie d'un sujet est provoquée par des entités non-sujets telles qu'un filtre et des caches et le sujet présent dans le passage de rayons X. La longueur de projection pro_ex des entités non-sujets (égale à un produit d'un coefficient d'absorption de rayons X par une longueur de transmission de rayons X) et la longueur de projection pro_pt du sujet expriment respectivement la propriété des entités non-sujets et la propriété du sujet, et on estime qu'elles sont en corrélation avec une ampleur de dispersion S(do,d). Par conséquent, la corrélation entre la longueur de projection pro_ex des entités non-sujets, la longueur de projection pro_pt de l'objet à radiographier et l'ampleur de dispersion S(do,d) est définie et les informations sur la corrélation sont enregistrées en association avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier.
Comme information, il est possible d'enregistrer directement l'ampleur de dispersion S(do,d). Autrement, on dira indirectement, une formule de fonction servant à calculer l'ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier ou des paramètres (coefficients) qui définissent la formule de fonction peuvent être enregistrés. Le rapport de l'ampleur de dispersion S(do,d) à des données ou une valeur calculées en soustrayant le rapport de 1, c'est-à-dire le rapport d'une valeur calculée en soustrayant des données l'ampleur de dispersion S(do,d), peut être enregistré.
De plus, les informations sur la corrélation peuvent être enregistrées en association qualifiée d'indirecte avec un paramètre calculé d'après la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier au lieu d'être enregistrées en association directe avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier. Le paramètre peut, par exemple, être une somme de longueurs de projection pro_pt+pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier (informations sur chaque détecteur). Autrement, le paramètre peut être une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs appartenant à tous les groupes de détecteurs et étant présents dans le canal dans lequel se trouve un détecteur concerné (informations sur chaque canal), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné et qui sont disposés dans le canal dans lequel se trouve le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs disposés dans tous les canaux et appartenant au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque groupe de détecteurs), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+ pro_ex calculées par rapport à des détecteurs disposés dans une pluralité de canaux dont le canal dans lequel se trouve le détecteur concerné et qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+ pro_ex calculées par rapport à des détecteurs disposés dans une pluralité de canaux dont le canal dans lequel se trouve le détecteur concerné et qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), ou une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à tous les détecteurs (informations partagées par tous les détecteurs).
Un sujet est radiographié afin de produire des données, et la longueur de projection pro_pt du sujet est calculée à l'aide des données. De plus, la longueur de projection pro_ex des entités non-sujets ayant affectée les données est calculée (la longueur de projection pro_ex est déterminée avec la relation géométrique entre un détecteur ayant produit les données et les entités non-sujets). Ensuite, les informations enregistrées en association avec la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets sont extraites. Par conséquent, une ampleur de dispersion peut être déduite à partir des informations, et les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le deuxième aspect de la présente invention, il est proposé un procédé de compensation de dispersion selon lequel les informations nécessaires au calcul d'une ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex d'entités non-sujets, c'est-à-dire des entités autres qu'un sujet, et de la longueur de projection pro_pt d'un objet à radiographier sont enregistrées, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée lorsqu'elle se trouve dans le passage de rayons X pendant la radiographie du sujet. Le sujet est radiographié afin de produire des données. Les données servant à calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets ayant affecté les données. Les informations associées à la longueur de projection pro_pt du sujet et à la longueur de projection pro_ex des entités non-sujets sont extraites et servent à corriger les données afin de compenser la dispersion.
Dans le procédé de compensation de dispersion selon le deuxième aspect, un sujet est radiographié afin de produire des données. Les données servant à calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets ayant affecté les données. Les informations enregistrées en association avec la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets sont extraites. Par conséquent, une ampleur de dispersion peut être déduite à partir des informations. Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le troisième aspect de la présente invention, il est proposé un procédé de compensation de dispersion identique au procédé de compensation de dispersion selon le premier ou le deuxième aspect, sauf que les informations concernent une fonction de dispersion f(pro_pt+ pro_ex) destinée à servir à calculer l'ampleur de dispersion S(do,d) d'après une somme de longueurs de projection pro_pt+pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier.
Dans le procédé de compensation de dispersion selon le troisième aspect, la fonction de dispersion f(pro_pt+pro_ex) enregistrée en association avec la somme pro_pt+pro_ex de la longueur de projection pro_pt du sujet et de la longueur de projection pro_ex des entités non-sujets est extraite. Par conséquent, l'ampleur de dispersion S(do,d) peut être calculée à l'aide de la fonction de dispersion f(pro_pt+pro_ex). Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le quatrième aspect de la présente invention, il est proposé un procédé de compensation de dispersion identique au procédé de compensation de dispersion selon le troisième aspect, sauf que la fonction de dispersion f(pro_pt+pro_ex) est une fonction d'ordre n (où n>1) de la somme de longueurs de projection pro_pt+pro_ex.
Les informations portent sur tout ou partie de coefficients aO, ... , an, inclus dans la fonction d'ordre n.
Dans le procédé ci-dessus, l'expression "partie des coefficients aO, ..., an" signifie qu'un coefficient d'une valeur 0 doit être enregistré. Par exemple, mais pour un terme constant, il n'est pas nécessaire de stocker le coefficient aO.
Dans le procédé de compensation de dispersion selon le quatrième aspect, une fois que sont extraits les coefficients aO, ..., an, enregistrés en association avec la somme pro_pt+pro_ex de la longueur de projection pro_pt du sujet et de la longueur de projection pro_ex des entités nonsujets, l'ampleur de dispersion S(do,d) peut être calculée à l'aide d'une fonction d'ordre n définie avec les coefficients d'extraction. Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le cinquième aspect de la présente invention, il est proposé un procédé de compensation de dispersion identique au procédé de compensation de dispersion selon l'un quelconque des premier à quatrième aspects, sauf que les informations servant à calculer une ampleur de dispersion par rapport à chaque détecteur, distingué à l'aide d'un numéro de groupe de détecteurs et d'un numéro de canal, et l'ampleur de dispersion est soustraite des données produites par chaque détecteur. Ainsi, les données sont corrigées afin de compenser la dispersion.
Comme l'ampleur de dispersion est soustraite des données produites, l'effet défavorable de la dispersion peut être supprimé.
Selon le sixième aspect de la présente invention, il est proposé un procédé de compensation de dispersion identique au procédé de compensation de dispersion selon l'un quelconque des premier à cinquième aspects, sauf que les informations sont enregistrées en association avec deux valeurs différentes, ou plus, d'une valeur de sortie Io de tube radiogène et/ou d'une épaisseur do de détecteur, et/ou d'une épaisseur d de faisceau, et/ou d'un numéro r de groupe de détecteurs et/ou d'un numéro ch de canal.
Dans le procédé de compensation de dispersion selon le sixième aspect, par exemple, même si le signal de sortie Io du tube radiogène, l'épaisseur do du détecteur ou l'épaisseur d du faisceau est modifiée, la dispersion peut être compensée en conséquence.
Selon le septième aspect de la présente invention, il est proposé un procédé de mesure de dispersion comprenant: une étape consistant à acquérir la longueur de 35 projection pro_ex d'entités non-sujets, c'està-dire des entités autres qu'un sujet, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée si elles sont présentes dans le passage de rayons X pendant la radiographie du sujet; une étape consistant à acquérir la longueur de projection pro_pt d'un objet à radiographier; une étape consistant à introduire les entités nonsujets et l'objet à radiographier dans le passage de rayons X, mesurer des données I(do,do) en effectuant une radiographie avec une épaisseur de faisceau établie à la même valeur que l'épaisseur do de détecteur, mesurer des données I(do,d) en radiographiant l'objet à radiographier avec l'épaisseur de faisceau établie à une valeur d supérieure à l'épaisseur do de détecteur, et calculer une ampleur de dispersion S(do,d) d'après la différence entre les données 1 o I(do,do) et les données I(do, d) ; et une étape consistant à enregistrer les informations nécessaires pour calculer l'ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier.
Dans le procédé de mesure de dispersion selon le septième aspect, l'épaisseur do de détecteur est maintenue inchangée, seule l'épaisseur d de faisceau est modifiée. Par conséquent, on estime qu'un incrément des données I(do,do) aux données I(do,d) est imputable à la dispersion. Par conséquent, l'ampleur de dispersion S(do,d) peut être calculée d'après la différence entre les données I(do,do) et les données I(do,d). Par ailleurs, il n'est pas tenu compte de la dispersion représentée par les données I(do,do). Autrement, comme les données I(do,do) peuvent être considérées comme des données produites par un détecteur à un seul canal, le procédé de compensation de dispersion selon la technique antérieure peut être employé pour corriger les informations sur la dispersion contenues dans les données I(do,do).
La dispersion de rayons X survenant pendant la radiographie d'un sujet est provoquée par des entités non-sujets telles qu'un filtre et des caches présents dans le passage de rayons X, et par le sujet. La longueur de projection pro_ex des entités non-sujets (égale à un produit d'un coefficient d'absorption de rayons X par une longueur de transmission de rayons X) et la longueur de projection pro_pt du sujet expriment respectivement la propriété des entités non-sujets et la propriété du sujet, et on estime qu'elles sont en corrélation avec l'ampleur de dispersion S(do,d). Par conséquent, la corrélation entre la longueur de projection pro_ex des entités non-sujets, la longueur de projection pro_pt de l'objet à radiographier et l'ampleur de dispersion S(do,d) définie, et les informations sur la corrélation sont enregistrées en association avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier.
Comme les informations sur la corrélation, l'ampleur de dispersion S(do,d) peut être directement enregistrée. Autrement, on dira indirectement, une formule de fonction servant à calculer l'ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier, ou des paramètres (coefficients ou analogues) définissant la formule de fonction peuvent être enregistrés. Autrement, le rapport de l'ampleur de dispersion S(do,d) à des données ou une valeur calculées en soustrayant le rapport de 1, c'est-à-dire le rapport de la valeur calculée en soustrayant des données l'ampleur de dispersion S(do,d), peut être enregistré.
De plus, les informations sur la corrélation peuvent être enregistrées en association dite indirecte avec un paramètre calculé d'après la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier au lieu d'être enregistrées en association directe avec la longueur de projection pro_ex des entités nonsujets et la longueur de projection pro_pt de l'objet à radiographier. Par exemple, le paramètre peut être une somme de longueurs de projection pro_pt+pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier (informations que chaque détecteur). Autrement, le paramètre peut être une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs appartenant à tous les groupes de détecteurs et disposés dans le canal dans lequel se trouve un détecteur concerné (informations sur chaque canal), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné et qui sont disposés dans le canal dans lequel se trouve le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs qui sont disposés dans tous les canaux et qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque groupe de détecteurs), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs qui sont disposés dans une pluralité de canaux dont le canal dans lequel se trouve le détecteur concerné et qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à des détecteurs disposés dans une pluralité de canaux dont le canal dans lequel se trouve le détecteur concerné et qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), ou une valeur calculée en additionnant les sommes de longueurs de projection pro_pt+pro_ex calculées par rapport à tous les détecteurs (informations partagées par tous les détecteurs).
Selon le huitième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon le septième aspect, sauf que l'étape d'acquisition de la longueur de projection pro_ex des entités non-sujets est une étape consistant à mesurer des données I_exout en effectuant une radiographie avec les entités non-sujets non introduites dans le passage de rayons X, mesurer des données I_exin en effectuant une radiographie avec les entités non-sujets introduites dans le passage de rayons X, et calculer la longueur de projection pro_ex des entités non-sujets à l'aide des données I_exout et des données I_exin.
Dans le procédé de mesure de dispersion selon le huitième aspect, les données I_exout sont mesurées avec les entités non-sujets non introduites dans le passage de rayons X et les données I_exin sont mesurées avec les entités non-sujets introduites dans le passage de rayons X. Par conséquent, on estime que la différence entre les données I_exout et les données I_exin est imputable à la longueur de projection pro_ex des entités non-sujets. Par conséquent, la longueur de projection pro_ex des entités non-sujets peut être calculée à l'aide des données I_exout et des données I_exin.
Selon le neuvième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon le septième aspect, sauf que l'étape d'acquisition de la longueur de projection pro_ex des entités non-sujets est une étape consistant à calculer la longueur de projection pro_ex des entités nonsujets d'après les matières et les formes des entités non-sujets.
Dans le procédé de mesure de dispersion selon le neuvième aspect, un coefficient d'absorption de rayons X est appris d'après les matières des entités non-sujets, et une longueur de transmission de rayons X est déduite d'après les formes de celles-ci. La longueur de projection pro_ex des entités non-sujets peut être calculée en multipliant le coefficient d'absorption de rayons X par la longueur de transmission de rayons X. Selon le dixième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon l'un quelconque des septième à neuvième aspects, sauf que l'étape d'acquisition de la longueur de projection pro_pt de l'objet à radiographier est une étape consistant à mesurer des données I_ptout en effectuant une radiographie avec l'objet à radiographier non introduit dans le passage de rayons X, mesurer des données I_ptin en radiographiant l'objet à radiographier avec l'objet à radiographier introduit dans le passage de rayons X, et calculer la longueur de projection pro_pt de l'objet à radiographier à l'aide des données I_ptout et des données I_ptin.
Dans le procédé de mesure de dispersion selon le dixième aspect, les données I_ptout sont mesurées avec l'objet à radiographier non introduit dans le passage de rayons X, et les données I_ptin sont mesurées avec l'objet à radiographier introduit dans celui-ci. Par conséquent, on estime que la différence entre les données I_ptout et les données I_ptin est imputable à la longueur de projection pro_pt de l'objet à radiographier. Par conséquent, la longueur de projection pro_pt de l'objet à radiographier peut être calculée à l'aide des données I_ptout et des données I_ptin.
Selon le onzième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon l'un quelconque des septième à dixième aspects, sauf qu'une fonction de dispersion f(pro_pt+pro_ex) destinée à servir à calculer l'ampleur de dispersion S(do,d) d'après la somme de longueurs de projection pro_pt+pro_ex, c'est-à-dire que la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier est définie et enregistrée comme information.
Dans le procédé de mesure de dispersion selon le onzième aspect, la fonction de dispersion f(pro_pt+pro_ex) servant à calculer l'ampleur de dispersion S(do,d) d'après la somme pro_pt+pro_ex de la longueur de projection pro_pt de l'objet et de la longueur de projection pro_ex des entités non-sujets est définie et enregistrée. Une fois que la fonction de dispersion f(pro_pt+pro_ex) est extraite, l'ampleur de dispersion S(do, d) peut être calculée. Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le douzième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon le onzième aspect, sauf que la fonction de dispersion f (pro_pt + pro_ex) est une fonction d'ordre n (où n> 1) de la somme de longueurs de projection pro_pt + pro_ex, et que tout ou partie de coefficients aO, ..., an, définissant al fonction d'ordre n sont enregistrées comme informations.
Dans le procédé ci-dessus, l'expression "partie des coefficients aO, ..., an" signifie qu'il n'est pas nécessaire d'enregistrer un coefficient à valeur O. Par exemple, mais pour un terme constant, il n'est pas nécessaire d'enregistrer le coefficient aO.
Dans le procédé de mesure de dispersion selon le douzième aspect, les coefficients aO, ..., an, définissant la fonction d'ordre n f (pro_pt + pro_ex) destinée à calculer l'ampleur de dispersion S(do,d) d'après la somme pro_pt + pro_ex de la longueur de projection pro_pt du sujet et de la longueur de projection pro_ex des entités non-sujets sont calculés et enregistrés. Une fois que les coefficients aO, ..., an sont extraits,l'ampleur de dispersion S(do,d) peut être calculée à l'aide de la fonction d'ordre n f (pro_pt + pro_ex). Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le treizième aspect de la présente invention, il est proposé un procédé de mesure de dispersion identique au procédé de mesure de dispersion selon l'un quelconque des septième à douzième aspects, sauf que la mesure est effectuée avec un signal de sortie Io de tube radiogène et/ou une épaisseur do et/ou une épaisseur d de faisceau et/ou un numéro r de groupe de détecteurs et/ou un numéro ch de canal établis à deux valeurs différentes ou plus, et les informations sont enregistrées en association avec chacune des valeurs.
Dans le procédé de mesure de dispersion selon le treizième aspect, par exemple, même si le signal de sortie Io de tube radiogène, l'épaisseur do de détecteur ou l'épaisseur d de faisceau est modifié, la dispersion peut être compensée en conséquence.
Selon le quatorzième aspect de la présente invention, il est proposé un système de tomodensitométrie comprenant: un tube radiogène; un détecteur de rayons X à plusieurs canaux; un moyen de balayage pour produire des données pendant la rotation relative du tube radiogène et/ou du détecteur de rayons X à plusieurs canaux autour d'un sujet; un moyen d'acquisition de longueur de projection d'entités non-sujets pour acquérir la longueur de projection pro_ex d'entités non-sujets, c'est-à-dire des entités autres qu'un sujet, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée pendant qu'elles se trouvent dans le passage de rayons X au cours d'une radiographie du sujet; un moyen d'acquisition de longueur de projection d'objet à radiographier servant à acquérir la longueur de projection pro_pt de l'objet à radiographier; un moyen d'acquisition d'ampleur de dispersion servant à introduire les entités non-sujets et l'objet à radiographier dans le passage de rayons X, à mesurer des données I(do,do) en effectuant une radiographie avec une épaisseur de faisceau établie à la même valeur que l'épaisseur do de détecteur, à mesurer des données I(do,d) en radiographiant l'objet à radiographier avec l'épaisseur de faisceau établie à une valeur d supérieure à l'épaisseur do de détecteur, et à calculer une ampleur de dispersion S(do,d) d'après la différence entre les données I(do,do) et les données I(do,d) ; un moyen d'enregistrement d'informations pour mémoriser des informations nécessaires au calcul de l'ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier; et un moyen de compensation de dispersion servant à irradier un faisceau de rayons X ayant l'épaisseur d de faisceau de façon à produire des données à l'aide d'un détecteur ayant l'épaisseur do de détecteur, en se servant des données pour calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets ayant affecté les données, à lire des informations associées aux longueurs de projection et à utiliser des informations pour corriger des données afin de compenser la dispersion.
Dans le système de tomodensitométrie selon le quatorzième aspect, l'épaisseur do de détecteur est maintenue inchangée, seule l'épaisseur d de faisceau étant modifiée. Par conséquent, on estime que l'incrément des données I(do,do) aux données I(do,d) est imputable à la dispersion. Par conséquent, l'ampleur de dispersion S(do,d) peut être calculée d'après la différence entre les données I(do,do) et les données I(do,d). A cet instant, on ne tient pas compte de la dispersion représentée par les données I(do,do). Autrement, comme les données I(do,do) peuvent être considérées comme des données produites par un détecteur de rayons X à un seul canal, le procédé de compensation de dispersion selon la technique antérieure peut être adopté pour corriger également les informations sur la dispersion contenues dans les données I(do,do).
Une dispersion de rayons X survenant au cours d'une radiographie d'une sujet est provoquée par des entités non-sujets telles qu'un filtre et des caches présents dans le passage de rayons X, et par un sujet. La longueur de projection proj_ex des entités non-sujets et la longueur de projection pro_pt du sujet expriment respectivement la propriété des entités nonsujets et la propriété du sujet, et on estime qu'elles sont en corrélation avec l'ampleur de dispersion S(do,d). Par conséquent, la corrélation entre la longueur de projection pro_ex des entités non-sujets, la longueur de projection pro_pt de l'objet à radiographier et l'ampleur de dispersion S(do,d) est définie et les informations sur la corrélation sont enregistrées en association avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier.
Comme informations, il est possible d'enregistrer directement l'ampleur de dispersion S(do,d) elle-même. Autrement, on dira indirectement, une formule de fonction destinée à servir au calcul de l'ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier ou un paramètre (coefficient ou analogue) définissant la formule de fonction peut être enregistré. Autrement, le rapport de l'ampleur de dispersion S(do,d) a des données ou a une valeur calculées en soustrayant le rapport de 1, c'est-à-dire le rapport d'une valeur calculé en soustrayant des données l'ampleur de dispersion S(do,d) peut être enregistré.
De plus, les informations sur la corrélation peuvent être enregistrées en association dite indirecte avec un paramètre déterminé avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier au lieu d'être enregistrées en association directe avec la longueur de projection pro_ex des entités nonsujets et la longueur de projection pro_pt de l'objet à radiographier. On peut adopter comme paramètre une somme de longueurs de projection pro_pt + pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier (informations sur chaque détecteur). Autrement, le paramètre peut être une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui sont disposés dans le canal dans lequel se trouve un détecteur concerné et qui appartiennent à tous les groupes de détecteurs (informations sur chaque canal), une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné et qui sont disposés dans le canal dans lequel se trouve le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui sont disposés dans tous les canaux et qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque groupe de détecteur), une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro-ex calculées par rapport à des détecteurs qui sont disposés dans une pluralité de canaux dont le canal dans lequel est disposé le détecteur concerné et qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui sont disposés dans une pluralité de canaux dont le canal dans lequel est disposé le détecteur concerné et qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteur auquel appartient le détecteur concerné (informations sur chaque détecteur, ou une valeur calculée en additionnant des sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à tous les détecteurs (informations partagées par tous les détecteurs).
Un sujet est radiographié afin de produire des données. Les données servent à calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex d'entités non-sujets ayant affecté les données (la longueur de projection pro_ex est déterminée avec la relation géométrique entre un détecteur ayant produit les données et les entités non-sujets). Une fois que sont extraites des informations enregistrées en association avec la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets, une ampleur de dispersion peut être déduite à partir des informations. Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le quinzième aspect de la présente invention, il est proposé un système de tomodensitométrie comprenant: un tube radiogène; un détecteur de rayons X à plusieurs canaux; un moyen de balayage pour produire des données pendant la rotation relative du tube radiogène et/ou du détecteur de rayons X à plusieurs canaux autour d'un sujet; un moyen d'enregistrement d'informations pour mémoriser des informations nécessaires afin de calculer une ampleur de dispersion S(do,d) à l'aide de la longueur de projection pro_ex d'entités non-sujets, à savoir des entités autres que le sujet, et la longueur de projection pro_pt de l'objet à radiographier, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée pendant leur présence dans le passage de rayons au cours d'une radiographie du sujet; et un moyen de compensation de dispersion servant à radiographier le sujet de manière à produire des données, à utiliser les données pour calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités nonsujets ayant affecté les données, à lire les informations associées aux deux longueurs de projection et à utiliser les informations pour corriger les données de manière à compenser la dispersion.
Dans le système de tomodensitométrie selon le quinzième aspect, un sujet est radiographié afin de produire des données. Les données servent à calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets ayant affecté les données. Une fois que sont extraites les informations enregistrées en association avec la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets, une ampleur de dispersion peut être déduite d'après les informations. Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le seizième aspect de la présente invention, il est proposé un système de tomodensitométrie identique au système de tomodensitométrie selon le quatorzième ou le quinzième aspect, sauf que les informations portent sur une fonction de dispersion f (pro_pt + pro_ex) destinée à servir au calcul de l'ampleur de dispersion S(do,d) d'après une somme de longueurs de projection pro_pt + pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt du sujet.
Dans le système de tomodensitométrie selon le seizième aspect, une fois que la fonction de dispersion f (pro_pt + pro_ex) enregistrée en association avec la somme pro_pt + pro_ex de la longueur de projection pro_pt du sujet et de la longueur de projection pro_ex des entités nonsujets est extraite, l'ampleur de dispersion S(do,d) peut être calculée à l'aide de la fonction de dispersion f (pro_pt + pro_ex). Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le dix-septième aspect de la présente invention, il est proposé un système de tomodensitométrie identique au système de tomodensitométrie selon le seizième aspect, sauf que la fonction de dispersion f (pro_pt + pro_ex) est une fonction d'ordre n (où n ?1) de la somme de longueurs de projection pro_pt + pro_ex, et les informations portent sur tout ou partie de coefficients aO, ..., an définissant la fonction d'ordre n.
Dans le système ci-dessus, l'expression "partie de coefficients aO, ..., an" signifie qu'il n'est pas nécessaire d'enregistrer un coefficient à valeur O. Par exemple, mais pour un terme constant, il n'est pas nécessaire d'enregistrer le coefficient aO.
Dans le système de tomodensitométrie selon le dix-septième aspect, les coefficients aO, ..., an définissant la fonction d'ordre n f (pro_pt + pro_ex) destinée à servir au calcul de l'ampleur de dispersion S (do,d) d'après la somme pro_pt + pro_ex de la longueur de projection pro_pt du sujet et de la longueur de projection pro_ex des entités non-sujets sont calculés et enregistrés. Une fois que sont extraits les coefficients aO, . .., an, l'ampleur de dispersion S(do,d) peut être calculée à l'aide de la fonction d'ordre n f (pro_pt + pro_ex). Finalement, les données produites peuvent être corrigées afin de compenser la dispersion.
Selon le dix-huitième aspect de la présente invention, il est proposé un système de tomodensitométrie identique au système de tomodensitométrie selon l'un quelconque des quatorzième à quinzième aspects, sauf que le moyen de compensation de dispersion utilise les informations pour calculer une ampleur de dispersion par rapport à chaque détecteur distingué par un numéro de groupe de détecteurs et un numéro de canal, et soustrait l'ampleur de dispersion des données produites par chaque détecteur. Ainsi, les données sont corrigées afin de compenser la dispersion.
En soustrayant des données l'ampleur de dispersion, on peut supprimer l'effet préjudiciable de la dispersion.
Selon le dix-neuvième aspect de la présente invention, il est proposé un système de tomodensitométrie identique au système de tomodensitométrie selon l'un quelconque des quatorzième à dix-huitième aspects, sauf que le moyen d'enregistrement d'informations mémorise les informations en association avec deux valeurs différentes, ou plus, d'un signal de sortie Io de tube radiogène et/ou d'une épaisseur do de détecteur et/ou d'une épaisseur d de faisceau et/ou d'un numéro r de groupe de détecteur et/ou d'un numéro ch de canal.
Dans le système de tomodensitométrie selon le dix-neuvième aspect, par exemple, même si le signal de sortie Io de tube radiogène, l'épaisseur do de détecteur ou le l'épaisseur d de faisceau est modifié, la dispersion peut être compensée.
Selon un procédé de compensation de dispersion et un système de tomodensitométrie dans lesquels est mise en oeuvre la présente invention, il est possible de compenser la dispersion survenant pendant une radiographie multicoupe.
De plus, selon un procédé de mesure de dispersion et un système de tomodensitométrie dans lesquels est mise en oeuvre la présente invention, il est possible de mesurer la dispersion survenant pendant une radiographie multicoupe.
L'invention sera mieux comprise à l'étude de la description détaillée d'un mode de réalisation pris à titre d'exemple non limitatif et illustré par les dessins annexés, sur lesquels: la Fig. 1 est un schéma de principe représentant un système de 5 tomodensitométrie selon la première forme de réalisation; la Fig. 2 est un organigramme décrivant une mesure de dispersion employée dans la première forme de réalisation; la Fig. 3 est un schéma explicatif présentant les circonstances dans lesquelles est réalisée une mesure de dispersion; la Fig. 4 est un schéma de principe illustrant une forme d'enregistrement de la longueur de projection d'entités non-sujets; la Fig. 5 est un schéma de principe présentant une forme d'enregistrement de coefficients de dispersion; la Fig. 6 est un schéma explicatif présentant les circonstances dans 15 lesquelles la radiographie d'un sujet est effectuée; la Fig. 7 est un organigramme décrivant la compensation de dispersion employée dans la première forme de réalisation; et la Fig. 8 est un organigramme décrivant la mesure de dispersion employée dans la deuxième forme de réalisation.
[Première forme de réalisation] La Fig. 1 est un schéma de principe représentant la configuration d'un système de tomodensitométrie 100 selon la première forme de réalisation.
Le système de tomodensitométrie 100 comprend un pupitre 1 d'opérateur, 25 une table de radiographie 10 et un portique de tomodensitométrie 20.
Le pupitre 1 d'opérateur comprend un dispositif de saisie 2 qui reçoit des données entrées par un opérateur, un dispositif de traitement central 3 qui réalise une mesure de dispersion, un traitement de balayage, une compensation de dispersion, un prétraitement, une reconstruction d'image et un post-traitement, une mémoire tampon 5 de regroupement de données qui regroupe des données acquises par le portique de tomodensitométrie 20, un dispositif d'affichage 6 sur lequel sont affichées des images tomographiques reconstruites d'après des données de projection produites en prétraitant les données regroupées, et un dispositif de mémorisation 7 dans lequel sont stockés des programmes, des données et des images tomographiques.
La table de radiographie 10 comporte un berceau 12 sur lequel est allongé un sujet et qui introduit et sort le sujet dans/de l'ouverture du portique de tomodensitométrie 20. Le berceau 12 est levé ou abaissé et déplacé de manière rectiligne par un moteur intégré dans la table de radiographie 10.
Le portique de tomodensitométrie 20 comprend un tube radiogène 21, un moyen de commande 22 de rayons X, un collimateur 23, un détecteur 24 de rayons X à plusieurs canaux, un système d'acquisition de données (SAD) 25, un moyen de commande de rotation 26 d'ensemble qui commande le tube radiogène 21 et d'autres éléments qui tournent autour d'un axe de rotation, un moyen de commande 29 qui transmet des signaux de commande et autres vers et depuis le pupitre 1 d'opérateur et la table de radiographie 10, et une bague collectrice 30 qui permet la circulation d'électricité, de signaux de commande et de données produites. De plus, un moyen de commande d'inclinaison 27 de portique de tomodensitométrie permet une inclinaison du portique de densitométrie 20 sur 30 vers l'avant ou vers l'arrière.
Un filtre 40 et des caches 41 et 42 sont disposés dans le passage de rayons X. Le filtre et les caches sont des entités non-sujets, à savoir des entités autres qu'un sujet, provoquant une dispersion de rayons X, qui ne peut pas être négligée, pendant leur présence dans le passage de rayons X au cours d'une radiographie du sujet.
La Fig. 2 est un organigramme décrivant la mesure de la dispersion selon la première forme de réalisation.
Lors de l'étape Pl, comme représenté sur la Fig. 3(a), de l'air est radiographié avec une épaisseur do de faisceau établie à la même valeur que l'épaisseur do d'un détecteur det (r, ch, do) avec le filtre 40 et les caches 41 et 42, c'est-à-dire les entités non-sujets, démontés. Par conséquent, des données I (Io, 0, r, ch, do, do, exout) sont produites.
Les paramètres r, ch, et do servant à distinguer un détecteur det, (r, ch, do) indiquent respectivement un groupe de détecteurs, un numéro de canal et une épaisseur de détecteur.
De plus, les paramètres Io, 0, r, ch, do, do et exout servant à distinguer les données I (Io, 0, r, ch, do, do, exout) désignent un signal de sortie de tube radiogène, un angle de vue, un numéro de groupe de détecteurs, un numéro de canal, une épaisseur de détecteur, une épaisseur de faisceau et l'absence d'entités non-sujets.
Lors de l'étape P2, comme représenté sur la Fig. 3(b), le filtre 40 et les caches 41 et 42, c'est-à-dire les entités non-sujets, sont montés, et l'air est radiographié avec l'épaisseur de faisceau établie à la même valeur do que l'épaisseur do du détecteur det (r, ch, do). Par conséquent, les données I (Io, 0, r, ch, do, do, exin) sont produites.
Par ailleurs, un paramètre exin figurant parmi les paramètres utilisés pour distinguer les données I (Io, 0, r, ch, do, do, exin) signifie la présence d'entités non-5 sujets.
Lors de l'étape P3, comme représenté sur la Fig. 3(c), un objet virtuel a est radiographié avec l'épaisseur de faisceau établi à la même valeur do que l'épaisseur do du détecteur det (r, ch, do) afin de produire les données I (Io, 0, r, ch, do, do, pt). Cette procédure est répétée à l'aide d'objets virtuels a de différentes dimensions.
Par ailleurs, un paramètre pt figurant parmi les paramètres utilisés pour distinguer les données I (Io, 0, r, ch, do, do, pt) signifie la présence d'un objet virtuel et d'entités non-sujets.
Lors de l'étape P4, comme représenté sur la Fig. 3(d), l'objet virtuel a est radiographié avec l'épaisseur du faisceau établie à une valeur d supérieure à l'épaisseur do du détecteur det (r, ch, do) afin de produire les données I (Io, 0, r, ch, do, do, pt). Cette procédure est répétée à l'aide d'objets virtuels a de différentes dimensions.
Lors de l'étape P5, la longueur de projection pro_ex (0, r, ch) d'entités non-sujets est calculée à l'aide des données I (Io, 0, r, ch, do, do, exout) et des données I(Io, 0, r, ch, do, do, exin), puis est enregistrée comme représenté sur la Fig. 4.
Pro_ex (0, r, ch) _ - log { I (Io, 0, r, ch, do, do, exin) / I (Io, 0, r, ch, do, do, exout) } Lors de l'étape P6, la longueur de projection pro_pt (0, r, ch) de l'objet virtuel a est calculée à l'aide des données I (Io, 0, r, ch, do, do, exin) et des données I(Io, 0, r, ch, do, do, pt). Cette procédure est répétée à l'aide de données produites en radiographiant les objets virtuels a de différentes dimensions, grâce à quoi on obtient des longueurs de projection pro_pt (0, r, ch) différentes.
pro_pt (0, r, ch) = - log { I (Io, 0, r, ch, do, do, pt) / I (Io, 0, r, ch, do, do, exin) } Lors de l'étape P7, la différence entre les données I (Io, 0, r, ch, do, do, pt) et les données I (Io, 0, r, ch, do, d, pt) est calculée comme ampleur de dispersion S (Io, 0, r, ch, do, d, pt). Cette procédure est répétée à l'aide des données produites en radiographiant les objets virtuels a de dimensions différentes, grâce à quoi les ampleurs de dispersion S sont obtenues en association avec différentes longueurs de projections pro_pt (0, r ch).
S (Io, 0, r, ch, do, d, pt) = I (Io, 0, r, ch, do, d, pt) / I (Io, 0, r, ch, do, do, pt) Lors de l'étape P8, l'ampleur de dispersion S est adaptée ou approchée d'une fonction d'ordre n d'une somme de longueurs de projection W=pro_pt (0, r, ch) + pro_ex (0, r, ch), à savoir la somme de la longueur de projection pro_ex (0, r, ch) des entités non-sujets et de la longueur de projection pro_pt (0, r, ch) de l'objet virtuel a, grâce à quoi des coefficients aO, ..., an, définissant la fonction d'ordre n sont calculés et enregistrés comme représenté sur la Fig. 5.
S=aO+al É W+a2 É W2+a3 É W3+... +an É W La Fig. 6 est un schéma de principe représentant les conditions de radiographie d'un sujet.
Pour commencer, comme représenté sur la Fig. 6(a), une radiographie est effectuée sans sujet afin de produire des données d'étalonnage et I (Io, 0, r, ch, do, d, air).
Ensuite, comme représenté sur la Fig. 6(b), un sujet (3 est introduit et radiographié afin de produire des données I (Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique.
La Fig. 7 est un organigramme décrivant la compensation de dispersion employée dans la première forme de réalisation.
Lors de l'étape Q1, la longueur de projection pro_pt (0, r, ch) du sujet f3 est calculée à l'aide des données d'étalonnage I (Io, 0, r, ch, do, d, air) et des données I(Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique.
pro_pt (0, r, ch) _ - log { I (Io, 0, r, ch, do, d, pt) / I (Io, 0, r, ch, do, d, air) } Lors de l'étape Q2, la longueur de projection pro_ex (0, r, ch) des entités non-sujets est extraite.
Lors de l'étape Q3, la somme des longueurs de projection W = pro_pt (0, r, ch) + pro_pt (0, r, ch) est calculée.
Lors de l'étape Q4, des coefficients de dispersion aO, al, ..., an, qui sont 30 associés à la somme de longueurs de projection W = pro_pt + pro_pt dans les conditions Io, 0, r, ch, do et d sont extraits.
Lors de l'étape Q5, les coefficients extraits aO, al, ..., an servent à calculer l'ampleur de dispersion S (Io, 0, r, ch, do, d, pt).
S =a0+al É W +a2 É W2+a3 É W3 + ... +an É W Lors de l'étape Q6, l'ampleur calculée de la dispersion S (Io, 0, r, ch, do, d, pt) et d'autres ampleurs de dispersion sont lissées dans la direction de groupes de détecteurs et la direction de canaux, grâce à quoi on obtient une ampleur de dispersion S' (Io, 0, r, ch, do, d, pt) lissée.
S' (Io, 0, r, ch, do, d, pt) = smooth rj_chk { S (Io, 0, r, ch, do, d, pt) } Ici, smooth rj_chk { S (Io, 0, r, ch, do, d, pt) } est une fonction (par exemple, une fonction servant à fournir une valeur moyenne) destinée à lisser des ampleurs de dispersion S (Io, 0, r -j, ch, do, d, pt) à s (Io, 0, r -i, ch, do, d, pt) calculées par rapport à des détecteurs det (r j, ch, do) à del (r+j, ch, do) qui sont disposés dans le même canal et qui sont situés dans la plage de j dans la direction de groupes de détecteurs avec, comme centre, un détecteur det (r, ch, do), par rapport auquel est calculée l'ampleur de dispersion S (Io, 0, r, ch, do, d, pt). De plus, conformément à la fonction smooth rj_chk { S (Io, 0, r, ch, do, d, pt) }, les ampleurs de dispersion S (Io, 0, r, ch-k, do, d, pt) à s (Io, 0, r, ch+k, do, d, pt) calculées par rapport aux détecteurs det (r, ch-k, do) à det (r, ch+k, do) qui appartiennent au même groupe de détecteurs et sont situées dans la plage de +k dans la direction de canaux, avec comme centre le détecteur det (r, ch, do), par rapport auquel est calculée l'ampleur de dispersion S (Io, 0r, ch, do, d, pt), sont lissées.
Lors de l'étape Q7, l'ampleur de dispersion S' (Io, 0, r, ch, do, d, pt) est soustraite des données I (Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique afin de calculer des données I' (Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique qui ont subi une compensation de dispersion. Il est ensuite mis fin à la compensation de dispersion.
I' (Io, 0, r, ch, do, d, pt) = I (Io, 0, r, ch, do, d, pt) S' (Io, 0, r, ch, do, d, pt) Selon la première forme de réalisation, une image tomographique reconstruite d'après les données I' (Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique, qui ont subi une compensation de dispersion, bénéficie d'une meilleure qualité d'image, ce qui donne un avantage par rapport à une image tomographique reconstruite d'après des données I (Io, 0, r, ch, do, d, pt) de reconstruction d'image tomographique qui n'ont pas subi une compensation de dispersion. En donnant à n une valeur 4 et à aO une valeur 0, le résultat de la compensation de dispersion est satisfaisant (dans ce cas, il n'est pas nécessaire d'enregistrer aO).
Selon la première forme de réalisation, la dispersion qui survient pendant une radiographie multicoupe est mesurée et compensée d'une manière préférée. Il en résulte une image multicoupe qui ne comporte plus d'artefacts imputables à la dispersion survenant pendant la radiographie multicoupe.
[Deuxième forme de réalisation] Selon la deuxième forme de réalisation, lalongueur de projection pro_ex des entités non-sujets est calculée mais non mesurée.
La Fig. 8 est un organigramme décrivant la mesure de dispersion employée dans la deuxième forme de réalisation.
La mesure de dispersion employée dans la deuxième forme de réalisation est différente de celle (Fig. 2) employée dans la première forme de réalisation en ce que l'étape Pl est exclue et l'étape P5' remplace l'étape P5. Seule l'étape P5' sera décrite ci-dessous.
Lors de l'étape P5', la longueur de projection pro_ex (0, r, ch) des entités non-sujets est calculée d'après un coefficient d'absorption de rayons X et une longueur de transmission L (0, r, ch) de rayons X présentés par les entités non-sujets, et elle est enregistrée comme représenté sur la Fig. 4.
Pro_ex (0, r, ch) = L (0, r, ch) Selon la deuxième forme de réalisation, il est inutile d'acquérir des données avec les entités non-sujets démontées.
[Troisième forme de réalisation] Selon la troisième forme de réalisation, les coefficients aO, ..., an, définissant la fonction d'ordre n de la somme de longueurs de projection W de laquelle est approchée l'ampleur de dispersion S sont utilisés en commun parmi toutes les vues ou tous les détecteurs.
En particulier, selon la première et la deuxième formes de réalisation, les coefficients aO, ..., an définissant la fonction d'ordre n sont calculés par rapport à chaque ensemble de conditions Io, O, r, ch,, do et d lors de l'étape P8 (Fig. 2 ou Fig. 8). Selon la troisième forme de réalisation, les coefficients a, ..., an définissant la fonction d'ordre n sont calculés par rapport à chaque ensemble de conditions Io, do et d, puis sont enregistrés.
Selon la troisième forme de réalisation, il est nécessaire d'enregistrer les coefficients aO, ..., an en association avec chaque angle de vue O et chaque détecteur det, (r, ch, do). Il en résulte une procédure simple.
[Quatrième forme de réalisation] Selon les première à troisième formes de réalisation, les coefficients aO, ..., an définissant la fonction d'ordre n de la somme de longueurs de projection W, de laquelle est approchée l'ampleur de dispersion S, sont enregistrés. En revanche, une ampleur de dispersion S (Io, 0, r, ch, do, d, pt) elle-même peut être enregistrée.
Autrement, le rapport de l'ampleur de dispersion S (Io, 0, r, ch, do, d, pt) aux données I (Io, 0, r, ch, do, d, pt) (c'est-à-dire le rapport d'une ampleur de dispersion à des données) peut être enregistré. Autrement, une valeur calculée en soustrayant de 1 le rapport de l'ampleur de dispersion S (Io, 0, r, ch, do, d, pt) aux données I (Io, 0, r, ch, do, d, pt) peut être enregistrée.
[Cinquième forme de réalisation] Selon les première à quatrième forme de réalisation, les coefficients aO, an sont enregistrés en association avec une somme de longueurs de projection pro_pt + pro_ex, à savoir la somme de la longueur de projection pro_ex d'entités non-sujets et de la longueur de projection pro_pt d'un objet à radiographier. Selon une autre possibilité, les coefficients peuvent être enregistrés en association directe avec la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier. Autrement, les coefficients peuvent être enregistrés en association avec une valeur calculée en additionnant les sommes de longueurs de projection pro_pt pro_ex calculées par rapport à des détecteurs qui sont disposés dans le canal dans lequel est présent un détecteur concerné et qui appartient à tous les groupes de détecteurs (informations sur chaque canal), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui sont disposés dans le canal dans lequel est présent le détecteur concerné et qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné et qui sont disposés dans tous les canaux (informations sur chaque groupe de détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui appartiennent au groupe de détecteurs auquel appartient le détecteur concerné et qui sont présents dans une pluralité de canaux dont le canal dans lequel est disposé le détecteur concerné (informations sur chaque détecteur), une valeur calculée en additionnant les sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à des détecteurs qui appartiennent à une pluralité de groupes de détecteurs dont le groupe de détecteurs auquel appartient le détecteur concerné et qui sont disposés dans une pluralité de canaux dont le canal dans lequel est disposé le détecteur concerné (informations sur chaque détecteur), ou une valeur calculée en additionnant les sommes de longueurs de projection pro_pt + pro_ex calculées par rapport à tous les détecteurs (informations partagées par tous les détecteurs). FIG.1
100: SYSTEME DE TOMODENSITOMETRIE 1: PUPITRE D'OPERATEUR 6: DISPOSITIF D'AFFICHAGE 2: DISPOSITIF DE SAISIE 3: DISPOSITIF DE TRAITEMENT CENTRAL 5: MEMOIRE TAMPON DE REGROUPEMENT DE DONNEES 7: DISPOSITIF DE MEMORISATION 20: PORTIQUE DE TOMODENSITOMETRIE 21: TUBE RADIOGENE 40: FILTRE 23: COLLIMATEUR 41: CACHE 26: MOYEN DE COMMANDE DE ROTATION D'ENSEMBLE 22: MOYEN DE COMMANDE DE RAYONS X 42: CACHE 24: DETECTEUR DE RAYONS X A PLUSIEURS CANAUX 30: BAGUE COLLECTRICE 29:MOYEN DE COMMANDE 20 27: MOYEN DE COMMANDE D'INCLINAISON DE PORTIQUE DE TOMODENSITOMETRIE 33: TABLEAU DE COMMANDE DE TOMODENSITOMETRIE 10: TABLE DE RADIOGRAPHIE 12: BERCEAU FIG. 2
LANCER LA MESURE DE DISPERSION
P 1: RADIOGRAPHIER L'AIR AVEC LES ENTITES NON-SUJETS DEMONTEES ET AVEC UNE EPAISSEUR DE FAISCEAU ETABLIE A LA MEME VALEUR QUE L'EPAISSEUR do DE DETECTEUR, ET PRODUIRE LES DONNEES I (Io, 0, r, ch, do, do, exout) P2: RADIOGRAPHIER L'AIR AVEC LES ENTITES NON-SUJETS MONTEES ET AVEC UNE EPAISSEUR DE FAISCEAU ETABLIE A LA MEME VALEUR QUE L'EPAISSEUR do DE DETECTEUR, ET PRODUIRE DES DONNEES I (Io, 0, r, ch, do, do, exin) P3: RADIOGRAPHIER L'OBJET A RADIOGRAPHIER AVEC UNE EPAISSEUR DE FAISCEAU ETABLIE A LA MEME VALEUR QUE L'EPAISSEUR do DE DETECTEUR, ET PRODUIRE LES DONNEES I(Io, 0, r, ch, do, do, pt) P4: RADIOGRAPHIER L'OBJET A RADIOGRAPHIER AVEC UNE EPAISSEUR d DE FAISCEAU ETABLIE A UNE VALEUR PLUS GRANDE QUE L'EPAISSEUR do DE DETECTEUR, ET PRODUIRE LES DONNEES I (Io, 0, r, ch, do, d, pt) P5: CALCULER LA LONGUEUR DE PROJECTION pro_ex (0, r, ch) DES ENTITES NON-SUJETS A L'AIDE DES DONNEES I (Io, 0, r, ch, do, do, exout) ET DES DONNEES I (Io, 0, r, ch, do, do, exin), ET ENREGISTRER LA LONGUEUR DE PROJECTION CALCULEE P6: CALCULER LA LONGUEUR DE PROJECTION pro_pt (0, r, ch) DE L'OBJET A RADIOGRAPHIER A L'AIDE DES DONNEES I (Io, 0, r, ch, do, do, 15 exin) ET DES DONNEES I (Io, 0, r, ch, do, do, pt) P7: CALCULER L'AMPLEUR DE DISPERSION S (Io, 0, r, ch, do, do, pt) A L'AIDE DES DONNEES I (Io, 0, r, ch, do, do, pt) ET DES DONNEES I (Io, 0, r, ch, do, d, pt) P8: DETERMINER LA FONCTION DE DISPERSION f (pro_pt + pro_ex) QUI DONNE L'AMPLEUR DE DISPERSION S (Io, 0, r, ch, do, d, pt) DANS DES CONDITIONS DE Io, 0, r, ch, do ET d, ET STOCKER ET ENREGISTRER LES COEFFICIENTS DE DISPERSION a0, al, ..., an,
FIN
FIG. 3 0: DIRECTION DU FAISCEAU EN EVENTAIL do: EPAISSEUR DE FAISCEAU IC: AXE GEOMETRIQUE DE ROTATION H: DIRECTION DE GROUPES DE DETECTEURS 40: FILTRE 41: CACHE 42: CACHE a: OBJET VIRTUEL FIG. 6 0: DIRECTION DE FAISCEAU EN EVENTAIL 41: FILTRE d: EPAISSEUR DE FAISCEAU IC: AXE GEOMETRIQUE DE ROTATION 41: CACHE 42: CACHE f3: SUJET H: DIRECTION DE GROUPES DE DETECTEURS FIG. 7
LANCER LA COMPENSATION DE DISPERSION
Q1: CALCULER LA LONGUEUR DE PROJECTION pro_pt (0, r, ch) DU SUJET A L'AIDE DES DONNEES I (Io, 0, r, ch, do, d, air) ET DES DONNEES I (Io, 0, r, ch, do, d, pt) Q2: LIRE LA LONGUEUR DE PROJECTION pro_ex (0, r, ch) DES ENTITES 10 NON-SUJETS Q3: CALCULER LA SOMME DES LONGUEURS DE PROJECTION W = pro_pt + pro_ex Q4: LIRE LES COEFFICIENTS DE DISPERSION a0, al, ..., an, ASSOCIES A LA SOMME DE LONGUEURS DE PROJECTION W = pro_pt + pro_ex DANS 15 LES CONDITIONS Io, 0, r, ch, do, ET d Q5: CALCULER L'AMPLEUR DE DISPERSION S (Io, 0, r, ch, do, d, pt) Q6: LISSER L'AMPLITUDE DE DISPERSION S (Io, 0, r, ch, do, d, pt) ET D'AUTRES AMPLEURS DE DISPERSION DANS LA DIRECTION DE CANAUX ET LA DIRECTION DE GROUPES DE DETECTEURS Q7: SOUSTRAIRE L'AMPLEUR DE DISPERSION S' (Io, 0, r, ch, do, d, pt) DES DONNEES I (Io, 0, r, ch, do, d, pt) DE FACON A CALCULER DES DONNEES I' (Io, 0, r, ch, do, d, pt) QUI ONT SUBI UNE COMPENSATION DE DISPERSION
FIN
FIG. 8
LANCER LA MESURE DE DISPERSION
P2: RADIOGRAPHIER L'AIR AVEC LES ENTITES NON-SUJETS MONTEES ET AVEC L'EPAISSEUR DE FAISCEAU ETABLIE A LA MEME VALEUR QUE L'EPAISSEUR do DE DETECTEUR DE MANIERE A PRODUIRE DES DONNEES I (Io, 0, r, ch, do, do, exin) P3: RADIOGRAPHIER L'OBJET A RADIOGRAPHIER AVEC UNE EPAISSEUR DE FAISCEAU ETABLIE A LA MEME VALEUR QUE L'EPAISSEUR do DE DETECTEUR DE MANIERE A PRODUIRE DES DONNEES I (Io, 0, r, ch, do, do, pt) P4: RADIOGRAPHIER L'OBJET A RADIOGRAPHIER AVEC UNE EPAISSEUR DE FAISCEAU ETABLIE A UNE VALEUR d SUPERIEURE A L'EPAISSEUR do DE DETECTEUR DE MANIERE A PRODUIRE DES DONNEES I (Io, 0, r, ch, do, d, pt) P5' : CALCULER LA LONGUEUR DE PROJECTION pro_ex (0, r, ch) DES ENTITES NON-SUJETS A L'AIDE DU COEFFICIENT D'ABSORPTION DE RAYONS X ET DE LA LONGUEUR DE TRANSMISSION DE RAYONS X L (0, r, ch) DES ENTITES NON-SUJETS P6: CALCULER LA LONGUEUR DE PROJECTION pro pt (0, r, ch) DE 10 L'OBJET A RADIOGRAPHIER A L'AIDE DES DONNEES I (Io, 0, r, ch, do, do, exin) ET DES DONNEES I (Io, 0, r, ch, do, d, pt) P7: CALCULER L'AMPLEUR DE DISPERSION S (Io, 0, r, ch, do, d, pt) A L'AIDE DES DONNEES I (Io, 0, r, ch, do, do, pt) ET DES DONNEES I (Io, 0, r, ch, do, d, pt) P8: DETERMINER LA FONCTION DE DISPERSION f (pro pt + pro_ex) QUI DONNE L'AMPLEUR DE DISPERSION S (Io, 0, r, ch, do, d, pt) DANS LES CONDITIONS Io, 0, r, ch, do ET d, ET ENREGISTRER LES COEFFICIENTS DE DISPERSION a0, al, ..., an
FIN 15

Claims (9)

REVENDICATIONS
1. Procédé de compensation de dispersion, comprenant les étapes consistant à : acquérir la longueur de projection pro_ex d'entités non- sujets, à savoir des entités autres qu'un sujet, les entités provoquant une dispersion de rayons qui ne peut pas être négligée pendant qu'elle se trouve dans le passage de rayons X au cours de la radiographie du sujet; acquérir la longueur de projection pro_pt d'un objet à radiographier; mesurer des données I (do, do) en effectuant la radiographie avec les entités non-sujets et l'objet à radiographier introduits dans le passage de rayons X et avec une épaisseur de faisceau établie à la même valeur qu'une épaisseur do de détecteur; mesurer des données I (do, d) en radiographiant l'objet à radiographier avec l'épaisseur de faisceau établie à une valeur d supérieure à l'épaisseur do de détecteur; calculer une ampleur de dispersion S (do, d) d'après la différence entre les données I (do, do) et les données I (do, d) ; enregistrer les informations nécessaires au calcule de l'ampleur de dispersion S (do, d) à l'aide de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier; utiliser les données produites en radiographiant le sujet pour calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des 25 entités non-sujets ayant affecté les données; lire les informations associées à la longueur de projection pro_pt du sujet et à la longueur de projection pro_ex des entités non-sujets; et utiliser des informations pour corriger les données de manière à compenser la dispersion.
2. Procédé de compensation de dispersion comprenant les étapes consistant à : enregistrer les informations nécessaires au calcul d'une ampleur de dispersion S (do, d) en utilisant à la fois la longueur de projection pro_ex d'entités non-sujets, à savoir des entités autres qu'un sujet, et la longueur de projection pro_pt d'un objet à radiographier, les entités provoquant une dispersion de rayons X qui ne peut être négligée pendant qu'elles se trouvent dans le passage de rayons X au cours de la radiographie du sujet; utiliser les données produites en radiographiant le sujet pour calculer la longueur de projection pro_pt du sujet et la longueur de projection pro_ex des entités non-sujets ayant affecté les données; lire les informations associées à la longueur de projection pro_pt du sujet et à la longueur de projection pro_ex des entités nonsujets; et utiliser les informations pour corriger les données de manière à compenser la dispersion.
3. Procédé de compensation de dispersion selon la revendication 1, caractérisé en ce que les informations se rapportent à une fonction de dispersion f(pro_pt + pro_ex) à utiliser pour calculer l'ampleur de dispersion S (do, d) d'après une somme de longueurs de projection pro_pt + pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités non-sujets et de la longueur de projection pro_pt de l'objet à radiographier.
4. Procédé de compensation de dispersion selon la revendication 2, dans lequel les informations se rapportent à une fonction de dispersion f (pro_pt + pro_ex) destinée à servir au calcul de l'ampleur de dispersion S (do, d) d'après une somme de longueurs de projection pro_pt + pro_ex, c'est-à-dire la somme de la longueur de projection pro_ex des entités nonsujets et de la longueur de projection pro_pt de l'objet à radiographier.
5. Procédé de compensation de dispersion selon la revendication 3, dans lequel la fonction de dispersion f (pro_pt + pro_ex) est une fonction d'ordre n (où n>1) d'une somme de longueurs de projection pro_pt + pro_ex, et les informations se rapportent à tout ou partie des coefficients aO, ..., an, définissant la fonction d'ordre n.
6. Procédé de compensation de dispersion selon la revendication 4, dans lequel la fonction de dispersion f (pro_pt + pro_ex) est une fonction d'ordre n (où n21) d'une somme de longueurs de projection pro_pt + pro_ex, et les informations se rapportent à tout ou partie des coefficients aO, . .., an, définissant la fonction d'ordre n.
7. Procédé de mesure de dispersion, comprenant: une étape consistant à acquérir la longueur de projection pro_ex des entités non-sujets, à savoir des entités autres qu'un sujet, les entités provoquant une dispersion de rayons X qui ne peut pas être négligée pendant qu'ils se trouvent dans le passage de rayons X au cours de la radiographie du sujet; une étape consistant à acquérir la longueur de projection pro_pt d'un objet à radiographier; une étape consistant à introduire les entités non- sujets et l'objet à radiographier dans le passage de rayons X, mesurer les données I (do, do) en effectuant une radiographie avec une épaisseur de faisceau établie à la même valeur que l'épaisseur do de détecteur, mesurer les données I (do, d) en radiographiant l'objet à radiographier avec l'épaisseur de faisceau établie à une valeur d supérieure à l'épaisseur do de détecteur, et calculer une ampleur de dispersion S (do, d) d'après la différence entre les données I (do, do) et les données I (do, d) ; et une étape consistant à enregistrer les informations nécessaires au calcule de l'ampleur de dispersion S (do, d) en utilisant la longueur de projection pro_ex des entités non-sujets et la longueur de projection pro_pt de l'objet à radiographier.
8. Procédé de mesure de dispersion selon la revendication 7, dans lequel l'étape d'acquisition de la longueur de projection pro_ex des entités nonsujets est une étape consistant à mesurer des données I_exout en effectuant une radiographie avec les entités non-sujets non introduites dans le passage de rayons X, à mesurer des données I_exin en effectuant une radiographie avec les entités non-sujets introduites dans le passage de rayons X, et calculer la longueur de projection pro_ex des entités nonsujets à l'aide des données 1_exout et des données I_exin.
9. Procédé de mesure de dispersion selon la revendication 7, dans lequel l'étape d'acquisition de la longueur de projection pro_ex des entités nonsujets est une étape consistant à calculer la longueur de projection pro_ex des entités non-sujets d'après les matières et les formes des entités non-sujets.
FR0601817A 2005-03-01 2006-03-01 Procede de compensation de la diffusion et procede de mesure Withdrawn FR2885007A1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005056083A JP4336661B2 (ja) 2005-03-01 2005-03-01 Ｘ線ｃｔ装置および散乱測定方法

Publications (1)

Publication Number Publication Date
FR2885007A1 true FR2885007A1 (fr) 2006-10-27

Family

ID=36848338

Family Applications (1)

Application Number Title Priority Date Filing Date
FR0601817A Withdrawn FR2885007A1 (fr) 2005-03-01 2006-03-01 Procede de compensation de la diffusion et procede de mesure

Country Status (6)

Country Link
US (1) US7336759B2 (fr)
JP (1) JP4336661B2 (fr)
KR (1) KR20060099412A (fr)
CN (1) CN100563572C (fr)
DE (1) DE102006009762A1 (fr)
FR (1) FR2885007A1 (fr)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5097384B2 (ja) * 2006-11-09 2012-12-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置及び散乱補正方法
JP5220383B2 (ja) * 2007-10-29 2013-06-26 株式会社日立メディコ 放射線撮像装置
KR20100133950A (ko) * 2008-01-30 2010-12-22 더 리전츠 오브 더 유니버시티 오브 캘리포니아 동적인 제약들에 따른 오브젝트 주변의 사용을 통한 단층 촬영에 있어서의 양 감소 및 이미지 강화
EP2276408B1 (fr) 2008-05-08 2019-07-10 Arineta Ltd. Système d'imagerie par rayons x avec correction de rayonnement de diffusion et son procédé d'utilisation
JP5121889B2 (ja) 2010-07-22 2013-01-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔ装置
JP2013059511A (ja) * 2011-09-14 2013-04-04 Fujifilm Corp 画像処理装置
KR102003197B1 (ko) * 2017-10-17 2019-07-24 연세대학교 원주산학협력단 단일 방사선 영상의 산란선 교정 시스템 및 방법
CN108095745B (zh) * 2017-12-13 2020-12-04 上海联影医疗科技股份有限公司 一种医学成像设备以及一种获取用于评价医学成像设备的性能的图像的方法
CN108606805B (zh) * 2018-03-15 2022-02-11 东软医疗系统股份有限公司 一种散射校正方法、装置及设备
EP3886707A1 (fr) * 2018-11-30 2021-10-06 Accuray, Inc. Tomodensitométrie à faisceau conique hélicoïdal avec un détecteur décentré
US11166690B2 (en) 2020-03-19 2021-11-09 Accuray, Inc. Noise and artifact reduction for image scatter correction
US11647975B2 (en) 2021-06-04 2023-05-16 Accuray, Inc. Radiotherapy apparatus and methods for treatment and imaging using hybrid MeV-keV, multi-energy data acquisition for enhanced imaging
US11605186B2 (en) 2021-06-30 2023-03-14 Accuray, Inc. Anchored kernel scatter estimate
US11794039B2 (en) 2021-07-13 2023-10-24 Accuray, Inc. Multimodal radiation apparatus and methods
US11854123B2 (en) 2021-07-23 2023-12-26 Accuray, Inc. Sparse background measurement and correction for improving imaging

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5594241A (en) * 1979-01-11 1980-07-17 Hitachi Medical Corp Xxray transverse layer device
DE3608965A1 (de) * 1986-03-18 1987-10-01 Philips Patentverwaltung Verfahren zur bestimmung der raeumlichen struktur in einer schicht eines untersuchungsbereiches
US4881251A (en) * 1986-07-31 1989-11-14 Kabushiki Kaisha Toshiba Computed tomograph apparatus
DE69321244T3 (de) * 1993-12-24 2003-03-13 Agfa-Gevaert N.V., Mortsel Verfahren mit einer teildurchsichtigen Abschirmung zum Ausgleich der Röntgenbilddarstellung von Streustrahlen in Röntgenbildern
JP3426677B2 (ja) 1994-01-27 2003-07-14 株式会社日立メディコ Ｘ線ｃｔ装置
JPH08131431A (ja) 1994-11-08 1996-05-28 Hitachi Medical Corp Ｘ線ｃｔ装置
DE19721535C2 (de) * 1997-05-22 2001-09-06 Siemens Ag Röntgen-Computertomograph zur Erzeugung von Röntgenschattenbildern
JP4218908B2 (ja) 1998-04-24 2009-02-04 株式会社東芝 Ｘ線ｃｔ装置
DE10047720A1 (de) * 2000-09-27 2002-04-11 Philips Corp Intellectual Pty Vorrichtung und Verfahren zur Erzeugung eines Röntgen-Computertomogramms mit einer Streustrahlungskorrektur
US6568851B2 (en) * 2000-10-25 2003-05-27 Kabushiki Kaisha Toshiba X-ray CT scanner
EP1521947B1 (fr) * 2002-06-17 2018-05-16 Edax Inc. Procede de spectres de diffusion pour analyse par fluorescence des rayons x avec des composantes optiques
US6789943B2 (en) * 2002-11-12 2004-09-14 Ge Medical Systems Global Technology Company, Llc Method and apparatus for scatter measurement using an occluded detector ring
JP3942178B2 (ja) * 2003-07-29 2007-07-11 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｘ線ｃｔシステム
JP3961468B2 (ja) * 2003-09-19 2007-08-22 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 放射線計算断層画像装置およびそれに用いる放射線検出器
JP3919724B2 (ja) * 2003-09-19 2007-05-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 放射線計算断層画像装置および断層像データ生成方法
JP3950855B2 (ja) * 2004-01-07 2007-08-01 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 散乱測定方法、散乱補正方法およびｘ線ｃｔ装置
US7136450B2 (en) * 2004-05-26 2006-11-14 Analogic Corporation Method of and system for adaptive scatter correction in multi-energy computed tomography

Also Published As

Publication number Publication date
US7336759B2 (en) 2008-02-26
US20060198489A1 (en) 2006-09-07
KR20060099412A (ko) 2006-09-19
CN1846624A (zh) 2006-10-18
JP4336661B2 (ja) 2009-09-30
CN100563572C (zh) 2009-12-02
JP2006239003A (ja) 2006-09-14
DE102006009762A1 (de) 2006-09-07

Similar Documents

Publication Publication Date Title
FR2885007A1 (fr) Procede de compensation de la diffusion et procede de mesure
US6421552B1 (en) Methods and apparatus for estimating cardiac motion using projection data
FR2644590A1 (fr) Procede d'acquisition de donnees radiologiques et de reconstruction de structures correspondant a ce corps
FR2830741A1 (fr) Procede et appareil pour etalonner un systeme d'imagerie
US7801264B2 (en) Method for calibrating a dual -spectral computed tomography (CT) system
FR2878729A1 (fr) Procede d'evaluation de dose et appareil de tomodensitometrie.
BE1007766A3 (nl) Werkwijze en inrichting voor computer tomografie.
JPH0798041B2 (ja) オーバレンジ補正回路およびオーバレンジ補正方法
FR2930132A1 (fr) Systeme de mammographie et procede pour le faire fonctionner.
EP0360653A1 (fr) Procédé et système de correction des défauts d'images d'un scanner dus aux déplacements de ce dernier
EP0871143B1 (fr) Système de traitement d'une séquence d'images bruitées, et appareil d'examen médical incluant un tel système
FR2875994A1 (fr) Systeme et procede d'imagerie a l'aide de sources de rayons x monoenergetiques.
FR2876268A1 (fr) Procede de regulation de dose de rayons x et appareil de tomographie aux rayons x utilisant un ordinateur
FR2632749A1 (fr) Procede et systeme d'etalonnage d'un scanner a rayons x en utilisant un ou plusieurs etalons circulaires excentres
FR2896607A1 (fr) Procede et dispositif de controle de la qualite d'images de tomosynthese.
FR2666426A1 (fr) Procede de correction des mesures de densite optique effectuees sur un film radiographique.
JP3950855B2 (ja) 散乱測定方法、散乱補正方法およびｘ線ｃｔ装置
EP1274989B1 (fr) Procede d'amelioration d'un examen radiologique et dispositif pour la mise en oeuvre de ce procede
JP3270153B2 (ja) コンピュータトモグラフ
US10595803B2 (en) Signal processing method and imaging system for scatter correction in computed tomography
JP3846576B2 (ja) コンピュータ断層撮影装置
EP1537826B1 (fr) Procédé de traitement de signal dans un appareil de radiologie dentaire
EP1014078B1 (fr) Contrôle par radiographie d'un objet possédant un réseau cristallin
FR2870447A1 (fr) Dispositif de classification de pixels en mammographie a contraste accentue
JP3857002B2 (ja) Ｘ線断層撮影方法及び装置

Legal Events

Date Code Title Description
ST Notification of lapse

Effective date: 20081125