JP2006239003A - 散乱補正方法、散乱測定方法およびｘ線ｃｔ装置 - Google Patents

散乱補正方法、散乱測定方法およびｘ線ｃｔ装置 Download PDF

Info

Publication number
JP2006239003A
JP2006239003A JP2005056083A JP2005056083A JP2006239003A JP 2006239003 A JP2006239003 A JP 2006239003A JP 2005056083 A JP2005056083 A JP 2005056083A JP 2005056083 A JP2005056083 A JP 2005056083A JP 2006239003 A JP2006239003 A JP 2006239003A
Authority
JP
Japan
Prior art keywords
pro
data
scattering
projection length
ray
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
Application number
JP2005056083A
Other languages
English (en)
Other versions
JP4336661B2 (ja
Inventor
Masayasu 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
Priority to JP2005056083A priority Critical patent/JP4336661B2/ja
Priority to FR0601817A priority patent/FR2885007A1/fr
Priority to US11/365,249 priority patent/US7336759B2/en
Priority to CNB2006100820951A priority patent/CN100563572C/zh
Priority to DE200610009762 priority patent/DE102006009762A1/de
Priority to KR1020060020080A priority patent/KR20060099412A/ko
Publication of JP2006239003A publication Critical patent/JP2006239003A/ja
Application granted granted Critical
Publication of JP4336661B2 publication Critical patent/JP4336661B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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

【課題】マルチスライス撮影時における散乱を補正する。
【解決手段】散乱を生じさせる物体であって被検体以外の物体である被検体外物（４０〜４２）の投影長pro_exを求め、撮影対象（α）の投影長pro_ptを求め、検出器厚ｄｏと同等のビーム厚で撮影対象（α）を撮影してデータＩ(do,do)を測定し、検出器厚ｄｏより厚いビーム厚ｄで撮影対象（α）を撮影してデータＩ(do,d)を測定し、両データＩ(do,do)，Ｉ(do,d)の差に基づいて散乱量Ｓ(do,d)を求め、投影長和「pro_pt＋pro_ex」に対応付けて記憶しておく。被検体を撮影してデータを収集した後、被検体の投影長pro_ptを求めると共に対応する被検体外物の投影長pro_exを求め、投影長pro_ptと投影長pro_exの和に対応する散乱量Ｓ(do,d)を読み出し、データを散乱補正する。
【効果】散乱に起因するＣＴ画像の画質の低下を抑制できる。
【選択図】図３

Description

本発明は、散乱補正方法、散乱測定方法およびＸ線ＣＴ（Computer Tomography）装置に関し、さらに詳しくは、マルチスライス撮影時における散乱の影響を補正するための散乱補正方法、散乱量を測定するための散乱測定方法およびＸ線ＣＴ装置に関する。
従来、単列Ｘ線検出器を用いたＸ線ＣＴ装置における各種の散乱補正方法が知られている（例えば、特許文献１、特許文献２参照。）。
また、多列Ｘ線検出器を用いたＸ線ＣＴ装置に適用可能な散乱補正方法が提案されている（例えば、特許文献３参照。）。
特開平７−２１３５１７号公報 特開平８−１３１４３１号公報 特開平１１−２９９７６８号公報
複数の検出器列を有する多列Ｘ線検出器を用いてマルチスライス撮影を行う場合、個々の検出器厚（または検出器列厚）よりもビーム厚が大きくなるため、散乱の影響が大きくなる。
しかし、単列Ｘ線検出器を用いたＸ線ＣＴ装置における従来の散乱補正方法は、このような場合について考慮しておらず、適用できない問題点がある。
また、多列Ｘ線検出器を用いたＸ線ＣＴ装置における従来の散乱補正方法は、このような場合にも適用可能であるが、画像再構成を２回行う必要があり、処理の負担が大きい問題点がある。
そこで、本発明の目的は、マルチスライス撮影時における散乱の影響を補正するための散乱補正方法、散乱量を測定するための散乱測定方法およびＸ線ＣＴ装置を提供することにある。
第１の観点では、本発明は、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求め、撮影対象の投影長pro_ptを求め、前記被検体外物および前記撮影対象をＸ線通過経路に入れた状態で、検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を測定し、検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し、前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求め、前記被検体外物の投影長pro_exおよび前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶しておき、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め、前記被検体の投影長pro_ptと前記被検体外物の投影長pro_exに対応する前記情報を読み出し該情報を用いて前記データを散乱補正することを特徴とする散乱補正方法を提供する。
上記構成において、撮影対象は、通常はファントムであるが、被検体であってもよい。また、データＩ(do,do)，データＩ(do,d)，散乱量Ｓ(do,d)中の（,）は、(検出器厚，ビーム厚)を表している。
上記第１の観点による散乱補正方法では、検出器厚ｄｏは変えずにビーム厚ｄだけを変えるから、データＩ(do,do)からデータＩ(do,d)への増分はもっぱら散乱に起因すると考えられる。よって、データＩ(do,do)とデータＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求めることが出来る。なお、データＩ(do,do)に含まれる散乱は無視する。あるいは、データＩ(do,do)は単列Ｘ線検出器のデータと見なせるから、従来公知の散乱補正方法を適用してデータＩ(do,do)に含まれる散乱をさらに補正してもよい。
被検体撮影時のＸ線の散乱は、Ｘ線通過経路に存在しているフィルタやカバーなどの被検体外物および被検体によって生じる。そして、被検体外物の投影長（＝Ｘ線吸収係数×Ｘ線透過長）proj_exおよび被検体の投影長pro_ptは被検体外物の特性および被検体の特性を表しており、散乱量Ｓ(do,d)と相関があると考えられる。そこで、被検体外物の投影長pro_exと撮影対象の投影長pro_ptと散乱量Ｓ(do,d)の相関関係を求め、その情報を被検体外物の投影長pro_exと撮影対象の投影長pro_ptに対応付けて記憶しておく。
この情報としては、直接的に散乱量Ｓ(do,d)そのものを記憶してもよいし、いわば間接的に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を算出する関数式またはその関数式を規定するパラメータ（係数など）を記憶してもよいし、データ中の散乱量Ｓ(do,d)の割合またはその割合を「１」から引いた値つまりデータから散乱量Ｓ(do,d)を引いた量の割合を記憶してもよい。
また、被検体外物の投影長pro_exと撮影対象の投影長pro_ptとに直接的に対応付けて記憶しておく他に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから求められるパラメータにいわば間接的に対応付けて記憶しておいてもよい。このパラメータとしては、例えば、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」（検出器毎の情報となる）や、その投影長和「pro_pt＋pro_ex」を、対応する検出器と同一チャネルで全検出器列について加算した値（チャネル毎の情報となる）や、同一チャネルで同一検出器列を含む複数検出器列について加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列で全チャネルについて加算した値（検出器列毎の情報となる）や、同一検出器列で同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列を含む複数検出器列および同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、全検出器について加算した値（全検出器で一つの情報を共有することになる）が挙げられる。
被検体を撮影してデータを収集し、そのデータから被検体の投影長pro_ptを求めると共にそのデータに対応する被検体外物の投影長pro_exを求める（そのデータを得た検出器と被検体外物の幾何学的関係から対応する投影長pro_exが求まる）。そして、被検体の投影長pro_ptと被検体外物の投影長pro_exとに対応して記憶しておいた情報を読み出せば、その情報から散乱量を推定できるから、収集したデータを散乱補正することが出来る。
第２の観点では、本発明は、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を求めるための情報を記憶しておき、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め、前記被検体の投影長pro_ptと前記被検体外物の投影長pro_exに対応する前記情報を読み出し該情報を用いて前記データを散乱補正することを特徴とする散乱補正方法を提供する。
上記第２の観点による散乱補正方法では、被検体を撮影してデータを収集し、そのデータから被検体の投影長pro_ptを求めると共にそのデータに対応する被検体外物の投影長pro_exを求める。そして、被検体の投影長pro_ptと被検体外物の投影長pro_exとに対応して記憶しておいた情報を読み出せば、その情報から散乱量を推定できるから、収集したデータを散乱補正することが出来る。
第３の観点では、本発明は、上記第１または上記第２の観点による散乱補正方法において、前記情報が、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)であることを特徴とする散乱補正方法を提供する。
上記第３の観点による散乱補正方法では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」に対応して記憶しておいた散乱関数ｆ(pro_pt＋pro_ex)を読み出せば、その散乱関数ｆ(pro_pt＋pro_ex)から散乱量Ｓ(do,d)を算出できるから、収集したデータを散乱補正することが出来る。
第４の観点では、本発明は、上記第３の観点による散乱補正方法において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報が前記ｎ次関数の係数ａ０，…，ａｎの全部または一部であることを特徴とする散乱補正方法を提供する。
上記構成において、「係数ａ０，…，ａｎの一部」とは、値が「０」の係数は記憶しておく必要がない、という意味である。例えば、定数項がなければ、係数ａ０を記憶しておく必要がない。
上記第４の観点による散乱補正方法では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」に対応して記憶しておいた係数ａ０，…，ａｎを読み出せば、その係数によるｎ次関数から散乱量Ｓ(do,d)を算出できるから、収集したデータを散乱補正することが出来る。
第５の観点では、本発明は、上記第１から上記第４のいずれかの観点による散乱補正方法において、検出器列とチャネルによって特定される検出器毎に前記情報を用いて散乱量を求め、検出器毎の前記データから前記散乱量を減算することにより、前記データを散乱補正することを特徴とする散乱補正方法を提供する。
データから散乱量を減算することにより、散乱の影響を抑制することが出来る。
第６の観点では、本発明は、上記第１から上記第５のいずれかの観点による散乱補正方法において、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについての２種類以上の異なる値について前記情報を記憶しておくことを特徴とする散乱補正方法を提供する。
上記第６の観点による散乱補正方法では、例えばＸ線管出力Ｉｏを変えた場合や、検出器厚ｄｏを変えた場合や、ビーム厚ｄを変えた場合などに対応できる。
第７の観点では、本発明は、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求めるステップと、撮影対象の投影長pro_ptを求めるステップと、前記被検体外物および前記撮影対象をＸ線通過経路に入れ、検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を測定し、検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し、前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求めるステップと、前記被検体外物の投影長pro_exと前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶するステップとを有することを特徴とする散乱測定方法を提供する。
上記第７の観点による散乱測定方法では、検出器厚ｄｏは変えずにビーム厚ｄだけを変えるから、データＩ(do,do)からデータＩ(do,d)への増分はもっぱら散乱に起因すると考えられる。よって、データＩ(do,do)とデータＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求めることが出来る。なお、データＩ(do,do)に含まれる散乱は無視する。あるいは、データＩ(do,do)は単列Ｘ線検出器のデータと見なせるから、従来公知の散乱補正方法を適用してデータＩ(do,do)に含まれる散乱をさらに補正してもよい。
被検体撮影時のＸ線の散乱は、Ｘ線通過経路に存在しているフィルタやカバーなどの被検体外物および被検体によって生じる。そして、被検体外物の投影長（＝Ｘ線吸収係数×Ｘ線透過長）proj_exおよび被検体の投影長pro_ptは被検体外物の特性および被検体の特性を表しており、散乱量Ｓ(do,d)と相関があると考えられる。そこで、被検体外物の投影長pro_exと撮影対象の投影長pro_ptと散乱量Ｓ(do,d)の相関関係を求め、その情報を被検体外物の投影長pro_exと撮影対象の投影長pro_ptに対応付けて記憶しておく。
この情報としては、直接的に散乱量Ｓ(do,d)そのものを記憶してもよいし、いわば間接的に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を算出する関数式またはその関数式を規定するパラメータ（係数など）を記憶してもよいし、データ中の散乱量Ｓ(do,d)の割合またはその割合を「１」から引いた値つまりデータから散乱量Ｓ(do,d)を引いた量の割合を記憶してもよい。
また、被検体外物の投影長pro_exと撮影対象の投影長pro_ptとに直接的に対応付けて記憶しておく他に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから求められるパラメータにいわば間接的に対応付けて記憶しておいてもよい。このパラメータとしては、例えば、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」（検出器毎の情報となる）や、その投影長和「pro_pt＋pro_ex」を、対応する検出器と同一チャネルで全検出器列について加算した値（チャネル毎の情報となる）や、同一チャネルで同一検出器列を含む複数検出器列について加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列で全チャネルについて加算した値（検出器列毎の情報となる）や、同一検出器列で同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列を含む複数検出器列および同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、全検出器について加算した値（全検出器で一つの情報を共有することになる）が挙げられる。
第８の観点では、本発明は、上記第７の観点による散乱測定方法において、前記被検体外物の投影長pro_exを求めるステップが、前記被検体外物をＸ線通過経路に入れないで撮影してデータＩ_exoutを測定し、前記被検体外物をＸ線通過経路に入れて撮影してデータＩ_exinを測定し、前記データＩ_exoutおよび前記データＩ_exinより前記被検体外物の投影長pro_exを求めるステップであることを特徴とする散乱測定方法を提供する。
上記第８の観点による散乱測定方法では、被検体外物をＸ線通過経路に入れないときのデータＩ_exoutと入れたときのデータＩ_exinを測定するから、データＩ_exoutとデータＩ_exinの差はもっぱら被検体外物の投影長pro_exに起因すると考えられる。よって、データＩ_exoutとデータＩ_exinから被検体外物の投影長pro_exを求めることが出来る。
第９の観点では、本発明は、上記第７の観点による散乱測定方法において、前記被検体外物の投影長pro_exを求めるステップが、前記被検体外物の材質と形状より前記被検体外物の投影長pro_exを算出するステップであることを特徴とする散乱測定方法を提供する。
上記第９の観点による散乱測定方法では、被検体外物の材質よりＸ線吸収係数が判り、被検体外物の形状よりＸ線透過長が判るから、Ｘ線吸収係数×Ｘ線透過長により被検体外物の投影長pro_exを算出することが出来る。
第１０の観点では、本発明は、上記第７から上記第９のいずれかの観点による散乱測定方法において、前記撮影対象の投影長pro_ptを求めるステップが、撮影対象をＸ線通過経路に入れないで撮影してデータＩ_ptoutを測定し、前記撮影対象をＸ線通過経路に入れて撮影してデータＩ_ptinを測定し、前記データＩ_ptoutおよび前記データＩ_ptinより前記撮影対象の投影長pro_ptを求めるステップであることを特徴とする散乱測定方法を提供する。
上記第１０の観点による散乱測定方法では、撮影対象をＸ線通過経路に入れないときのデータＩ_ptoutと入れたときのデータＩ_ptinを測定するから、データＩ_ptoutとデータＩ_ptinの差はもっぱら撮影対象の投影長pro_ptに起因すると考えられる。よって、データＩ_ptoutとデータＩ_ptinから撮影対象の投影長pro_ptを求めることが出来る。
第１１の観点では、本発明は、上記第７から上記第１０のいずれかの観点による散乱測定方法において、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)を求め、前記情報として記憶することを特徴とする散乱測定方法を提供する。
上記第１１の観点による散乱測定方法では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」から散乱量Ｓ(do,d)を算出する散乱関数ｆ(pro_pt＋pro_ex)を求めて記憶しておく。この散乱関数ｆ(pro_pt＋pro_ex)を読み出すことにより、散乱量Ｓ(do,d)を算出でき、収集したデータを散乱補正することが出来る。
第１２の観点では、本発明は、上記第１１の観点による散乱測定方法において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報として前記ｎ次関数の係数ａ０，…，ａｎの全部または一部を記憶することを特徴とする散乱測定方法を提供する。
上記構成において、「係数ａ０，…，ａｎの一部」とは、値が「０」の係数は記憶しておく必要がない、という意味である。例えば、定数項がなければ、係数ａ０を記憶しておく必要がない。
上記第１２の観点による散乱測定方法では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」から散乱量Ｓ(do,d)を算出するｎ次関数ｆ(pro_pt＋pro_ex)の係数ａ０，…，ａｎを求めて記憶しておく。この係数ａ０，…，ａｎを読み出すことにより、ｎ次関数ｆ(pro_pt＋pro_ex)から散乱量Ｓ(do,d)を算出でき、収集したデータを散乱補正することが出来る。
第１３の観点では、本発明は、上記第７から上記第１２のいずれかの観点による散乱測定方法において、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについて２種類以上の異なる値で測定を行い、それぞれについて前記情報を記憶することを特徴とする散乱測定方法を提供する。
上記第１３の観点による散乱測定方法では、例えばＸ線管出力Ｉｏを変えた場合や、検出器厚ｄｏを変えた場合や、ビーム厚ｄを変えた場合などに対応できる。
第１４の観点では、本発明は、Ｘ線管と、多列Ｘ線検出器と、前記Ｘ線管および前記多列Ｘ線検出器の少なくとも一方を被検体の周りに相対回転させながらデータを収集するスキャン手段と、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求める被検体外物投影長取得手段と、撮影対象の投影長pro_ptを求める撮影対象投影長取得手段と、前記被検体外物および前記撮影対象をＸ線通過経路に入れて検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を収集すると共に検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求める散乱量取得手段と、前記被検体外物の投影長pro_exと前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶する情報記憶手段と、ビーム厚ｄのＸ線ビームで且つ検出器厚ｄｏの検出器でデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め両者に対応する情報を読み出して該情報を用いて前記データを散乱補正する散乱補正手段とを具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１４の観点によるＸ線ＣＴ装置では、検出器厚ｄｏは変えずにビーム厚ｄだけを変えるから、データＩ(do,do)からデータＩ(do,d)への増分はもっぱら散乱に起因すると考えられる。よって、データＩ(do,do)とデータＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求めることが出来る。なお、データＩ(do,do)に含まれる散乱は無視する。あるいは、データＩ(do,do)は単列Ｘ線検出器のデータと見なせるから、従来公知の散乱補正方法を適用してデータＩ(do,do)に含まれる散乱をさらに補正してもよい。
被検体撮影時のＸ線の散乱は、Ｘ線通過経路に存在しているフィルタやカバーなどの被検体外物および被検体によって生じる。そして、被検体外物の投影長proj_exおよび被検体の投影長pro_ptは被検体外物の特性および被検体の特性を表しており、散乱量Ｓ(do,d)と相関があると考えられる。そこで、被検体外物の投影長pro_exと撮影対象の投影長pro_ptと散乱量Ｓ(do,d)の相関関係を求め、その情報を被検体外物の投影長pro_exと撮影対象の投影長pro_ptに対応付けて記憶しておく。
この情報としては、直接的に散乱量Ｓ(do,d)そのものを記憶してもよいし、いわば間接的に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を算出する関数式またはその関数式を規定するパラメータ（係数など）を記憶してもよいし、データ中の散乱量Ｓ(do,d)の割合またはその割合を「１」から引いた値つまりデータから散乱量Ｓ(do,d)を引いた量の割合を記憶してもよい。
また、被検体外物の投影長pro_exと撮影対象の投影長pro_ptとに直接的に対応付けて記憶しておく他に、被検体外物の投影長pro_exと撮影対象の投影長pro_ptから求められるパラメータにいわば間接的に対応付けて記憶しておいてもよい。このパラメータとしては、例えば、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」（検出器毎の情報となる）や、その投影長和「pro_pt＋pro_ex」を、対応する検出器と同一チャネルで全検出器列について加算した値（チャネル毎の情報となる）や、同一チャネルで同一検出器列を含む複数検出器列について加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列で全チャネルについて加算した値（検出器列毎の情報となる）や、同一検出器列で同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列を含む複数検出器列および同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、全検出器について加算した値（全検出器で一つの情報を共有することになる）が挙げられる。
被検体を撮影してデータを収集し、そのデータから被検体の投影長pro_ptを求めると共にそのデータに対応する被検体外物の投影長pro_exを求める（そのデータを得た検出器と被検体外物の幾何学的関係から対応する投影長pro_exが求まる）。そして、被検体の投影長pro_ptと被検体外物の投影長pro_exとに対応して記憶しておいた情報を読み出せば、その情報から散乱量を推定できるから、収集したデータを散乱補正することが出来る。
第１５の観点では、本発明は、Ｘ線管と、多列Ｘ線検出器と、前記Ｘ線管および前記多列Ｘ線検出器の少なくとも一方を被検体の周りに相対回転させながらデータを収集するスキャン手段と、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を求めるための情報を記憶する情報記憶手段と、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め両者に対応する前記情報を読み出し該情報を用いて前記データを散乱補正する散乱補正手段とを具備したことを特徴とするＸ線ＣＴ装置を提供する。
上記第１５の観点によるＸ線ＣＴ装置では、被検体を撮影してデータを収集し、そのデータから被検体の投影長pro_ptを求めると共にそのデータに対応する被検体外物の投影長pro_exを求める。そして、被検体の投影長pro_ptと被検体外物の投影長pro_exとに対応して記憶しておいた情報を読み出せば、その情報から散乱量を推定できるから、収集したデータを散乱補正することが出来る。
第１６の観点では、本発明は、上記第１４または上記第１５の観点によるＸ線ＣＴ装置において、前記情報が、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)であることを特徴とするＸ線ＣＴ装置を提供する。
上記第１６の観点によるＸ線ＣＴ装置では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」に対応して記憶しておいた散乱関数ｆ(pro_pt＋pro_ex)を読み出せば、その散乱関数ｆ(pro_pt＋pro_ex)から散乱量Ｓ(do,d)を算出できるから、収集したデータを散乱補正することが出来る。
第１７の観点では、本発明は、上記第１６の観点によるＸ線ＣＴ装置において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報が前記ｎ次関数の係数ａ０，…，ａｎの全部または一部であることを特徴とするＸ線ＣＴ装置を提供する。
上記構成において、「係数ａ０，…，ａｎの一部」とは、値が「０」の係数は記憶しておく必要がない、という意味である。例えば、定数項がなければ、係数ａ０を記憶しておく必要がない。
上記第１７の観点によるＸ線ＣＴ装置では、被検体の投影長pro_ptと被検体外物の投影長pro_exの和「pro_pt＋pro_ex」から散乱量Ｓ(do,d)を算出するｎ次関数ｆ(pro_pt＋pro_ex)の係数ａ０，…，ａｎを求めて記憶しておく。この係数ａ０，…，ａｎを読み出すことにより、ｎ次関数ｆ(pro_pt＋pro_ex)から散乱量Ｓ(do,d)を算出でき、収集したデータを散乱補正することが出来る。
第１８の観点では、本発明は、上記第１４から上記第１７のいずれかの観点によるＸ線ＣＴ装置において、前記散乱補正手段は、検出器列とチャネルによって特定される検出器毎に前記情報を用いて散乱量を求め、検出器毎の前記データから前記散乱量を減算することにより、前記データを散乱補正することを特徴とするＸ線ＣＴ装置を提供する。
データから散乱量を減算することにより、散乱の影響を抑制することが出来る。
第１９の観点では、本発明は、上記第１４から上記第１８のいずれかの観点によるＸ線ＣＴ装置において、前記情報記憶手段は、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについての２種類以上の異なる値について前記情報を記憶しておくことを特徴とするＸ線ＣＴ装置を提供する。
上記第１９の観点によるＸ線ＣＴ装置では、例えばＸ線管出力Ｉｏを変えた場合や、検出器厚ｄｏを変えた場合や、ビーム厚ｄを変えた場合などに対応できる。
本発明の散乱補正方法およびＸ線ＣＴ装置によれば、マルチスライス撮影時における散乱を補正することが出来る。
また、本発明の散乱測定方法およびＸ線ＣＴ装置によれば、マルチスライス撮影時における散乱を測定することが出来る。
以下、図に示す実施例により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
図１は、実施例１に係るＸ線ＣＴ装置１００を示す構成ブロック図である。
このＸ線ＣＴ装置１００は、操作コンソール１と、撮影テーブル１０と、走査ガントリ２０とを具備している。
操作コンソール１は、操作者の入力を受け付ける入力装置２と、散乱測定処理，スキャン処理，散乱補正処理，前処理，画像再構成処理，後処理などを実行する中央処理装置３と、走査ガントリ２０で取得したデータを収集するデータ収集バッファ５と、収集したデータを前処理して求めた投影データから再構成した断層像を表示する表示装置６と、プログラムやデータや断層像を記憶する記憶装置７とを具備している。
撮影テーブル１０は、被検体を乗せて走査ガントリ２０の開口部に入れ出しするクレードル１２を具備している。クレードル１２は、撮影テーブル１０に内蔵するモータで昇降および直線移動される。
走査ガントリ２０は、Ｘ線管２１と、Ｘ線コントローラ２２と、コリメータ２３と、多列Ｘ線検出器２４と、ＤＡＳ（Data Acquisition System）２５と、回転中心軸の周りを回転するＸ線管２１などを制御する回転部コントローラ２６と、制御信号などを操作コンソール１や撮影テーブル１０とやり取りする制御コントローラ２９と、電源，制御信号，収集したデータを転送するスリップリング３０とを具備している。また、走査ガントリ傾斜コントローラ２７により、走査ガントリ２０は前方または後方に±３０゜ほど傾斜させることが出来る。
Ｘ線通過経路には、フィルタ４０およびカバー４１，４２が存在している。これらは、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物である。
図２は、実施例１に係る散乱測定処理を示すフロー図である。
ステップＰ１では、図３の（ａ）に示すように、フィルタ４０およびカバー４１，４２すなわち被検体外物を外した状態で、検出器det(r,ch,do)の検出器厚ｄｏと同等のビーム厚ｄｏで空気を撮影し、データＩ(Io,θ,r,ch,do,do,exout)を収集する。
検出器det(r,ch,do)中の(r,ch,do)は、（検出器列番号，チャネル番号，検出器厚)を表している。
また、データＩ(Io,θ,r,ch,do,do,exout)中の(Io,θ,r,ch,do,do,exout)は、（Ｘ線管出力，ビュー角度，検出器列番号，チャネル番号，検出器厚，ビーム厚，被検体外物無し)を表している。
ステップＰ２では、図３の（ｂ）に示すように、フィルタ４０およびカバー４１，４２すなわち被検体外物を装着し、検出器det(r,ch,do)の検出器厚ｄｏと同等のビーム厚ｄｏで空気を撮影して撮影し、データＩ(Io,θ,r,ch,do,do,exin)を収集する。
なお、データＩ(Io,θ,r,ch,do,do,exin)中の(…,exin)は、（…，被検体外物有り)を表している。
ステップＰ３では、図３の（ｃ）に示すように、検出器det(r,ch,do)の検出器厚ｄｏと同等のビーム厚ｄｏでファントムαを撮影してデータＩ(Io,θ,r,ch,do,do,pt)を収集する。これを、異なる大きさのファントムαについて繰り返す。
なお、データＩ(Io,θ,r,ch,do,do,pt)中の(…,pt)は、（…，ファントム及び被検体外物有り)を表している。
ステップＰ４では、図３の（ｄ）に示すように、検出器det(r,ch,do)の検出器厚ｄｏより厚いビーム厚ｄでファントムαを撮影してデータＩ(Io,θ,r,ch,do,d,pt)を収集する。これを、異なる大きさのファントムαについて繰り返す。
ステップＰ５では、データＩ(Io,θ,r,ch,do,do,exout)とデータＩ(Io,θ,r,ch,do,do,exin)から被検体外物の投影長pro_ex(θ,r,ch)を算出する。そして、図４に示すように記憶する。
pro_ex(θ,r,ch)＝−log{Ｉ(Io,θ,r,ch,do,do,exin)／Ｉ(Io,θ,r,ch,do,do,exout)}
ステップＰ６では、データＩ(Io,θ,r,ch,do,do,exin)とデータＩ(Io,θ,r,ch,do,do,pt)からファントムαの投影長pro_pt(θ,r,ch)を算出する。これを、異なる大きさのファントムαについて収集したデータについて繰り返し、異なる投影長pro_pt(θ,r,ch)を得る。
pro_pt(θ,r,ch)＝−log{Ｉ(Io,θ,r,ch,do,do,pt)／Ｉ(Io,θ,r,ch,do,do,exin)}
ステップＰ７では、データＩ(Io,θ,r,ch,do,do,pt)とデータＩ(Io,θ,r,ch,do,d,pt)の差を散乱量Ｓ(Io,θ,r,ch,do,d,pt)として算出する。これを、異なる大きさのファントムαについて収集したデータについて繰り返し、異なる投影長pro_pt(θ,r,ch)についての散乱量Ｓを得る。
Ｓ(Io,θ,r,ch,do,d,pt)＝Ｉ(Io,θ,r,ch,do,d,pt)−Ｉ(Io,θ,r,ch,do,do,pt)
ステップＰ８では、被検体外物の投影長pro_ex(θ,r,ch)とファントムαの投影長pro_pt(θ,r,ch)の和である投影長和Ｗ＝pro_pt(θ,r,ch)＋pro_ex(θ,r,ch)のｎ次関数として散乱量Ｓを近似し、ｎ次関数の係数ａ０，…，ａｎを求める。そして、図５に示すように記憶する。
Ｓ＝ａ０＋ａ１・Ｗ＋ａ２・Ｗ2＋ａ３・Ｗ3＋…＋ａｎ・Ｗn
図６は、被検体の撮影状態を示す概念図である。
まず、図６の（ａ）に示すように、被検体を入れないで撮影し、キャリブレーション用データＩ(Io,θ,r,ch,do,d,air)を収集する。
次に、図６の（ｂ）に示すように、被検体βを入れて撮影し、断層像再構成用データＩ(Io,θ,r,ch,do,d,pt)を収集する。
図７は、実施例１に係る散乱補正処理を示すフロー図である。
ステップＱ１では、キャリブレーション用データＩ(Io,θ,r,ch,do,d,air)と断層像再構成用データＩ(Io,θ,r,ch,do,d,pt)から被検体βの投影長pro_pt(θ,r,ch)を算出する。
pro_pt(θ,r,ch)＝−log{Ｉ(Io,θ,r,ch,do,d,pt)／Ｉ(Io,θ,r,ch,do,d,air)}
ステップＱ２では、被検体外物の投影長pro_ex(θ,r,ch)を読み出す。
ステップＱ３では、投影長和Ｗ＝pro_pt(θ,r,ch)＋pro_pt(θ,r,ch)を算出する。
ステップＱ４では、条件［Io,θ,r,ch,do,d］での投影長和Ｗ＝pro_pt＋pro_ptに対応する散乱係数ａ０，ａ１，…，ａｎを読み出す。
ステップＱ５では、読み出した係数ａ０，ａ１，…，ａｎを用いて散乱量Ｓ(Io,θ,r,ch,do,d,pt)を算出する。
Ｓ＝ａ０＋ａ１・Ｗ＋ａ２・Ｗ2＋ａ３・Ｗ3＋…＋ａｎ・Ｗn
ステップＱ６では、算出した散乱量Ｓ(Io,θ,r,ch,do,d,pt)を検出器列方向およびチャネル方向に平滑化処理し、平滑化した散乱量Ｓ'(Io,θ,r,ch,do,d,pt)を算出する。
Ｓ'(Io,θ,r,ch,do,d,pt)＝smooth_rj_chk{Ｓ(Io,θ,r,ch,do,d,pt)}
ここで、smooth_rj_chk{Ｓ(Io,θ,r,ch,do,d,pt)}は、散乱量Ｓ(Io,θ,r,ch,do,d,pt)に対応する検出器det(r,ch,do)を中心として同一チャネルで検出器列方向に±ｊ列の検出器det(r-j,ch,do)〜det(r+j,ch,do)に対応する散乱量Ｓ(Io,θ,r-j,ch,do,d,pt)〜Ｓ(Io,θ,r+j,ch,do,d,pt)および散乱量Ｓ(Io,θ,r,ch,do,d,pt)に対応する検出器det(r,ch,do)を中心として同一検出器列でチャネル方向に±ｋ個の検出器det(r,ch-k,do)〜det(r,ch+k,do)に対応する散乱量Ｓ(Io,θ,r,ch-k,do,d,pt)〜Ｓ(Io,θ,r,ch+k,do,d,pt)を平滑化処理する関数（例えば平均値を求める関数）である。
ステップＱ７では、断層像再構成用データＩ(Io,θ,r,ch,do,d,pt)から散乱量Ｓ'(Io,θ,r,ch,do,d,pt)を減算し、散乱補正した断層像再構成用データＩ'(Io,θ,r,ch,do,d,pt)を算出する。そして、散乱補正処理を終了する。
Ｉ'(Io,θ,r,ch,do,d,pt)＝Ｉ(Io,θ,r,ch,do,d,pt)−Ｓ'(Io,θ,r,ch,do,d,pt)
実施例１で散乱補正した断層像再構成用データＩ'(Io,θ,r,ch,do,d,pt)から再構成した断層像は、散乱補正しない断層像再構成用データＩ(Io,θ,r,ch,do,d,pt)から再構成した断層像よりも画質が向上した。なお、ｎ＝４，ａ０＝０として良好な散乱補正結果が得られることを確認した（この場合、ａ０は記憶しておく必要がない）。
実施例１によれば、マルチスライス撮影時における散乱を好適に測定し、補正することが出来る。これにより、マルチスライス撮影時における散乱に起因するアーチファクトを抑制したマルチスライス画像を得ることが出来る。
実施例２では、被検体外物の投影長pro_exを、実測でなく、計算により求める。
図８は、実施例２に係る散乱測定処理を示すフロー図である。
この実施例２に係る散乱測定処理は、実施例１の散乱測定処理（図２）からステップＰ１を省略し、ステップＰ５の代わりにステップＰ５’を行う点のみが相違している。そこで、ステップＰ５’のみ説明する。
ステップＰ５’では、被検体外物のＸ線吸収係数μとＸ線透過長Ｌ(θ,r,ch)より被検体外物の投影長pro_ex(θ,r,ch)を算出し、図４に示すように記憶する。
pro_ex(θ,r,ch)＝μ・Ｌ(θ,r,ch)
実施例２によれば、被検体外物を外してデータを収集する必要がなくなる。
実施例３では、投影長和Ｗのｎ次関数として散乱量Ｓを近似するｎ次関数の係数ａ０，…，ａｎを全ビュー，全検出器で共通とする。
すなわち、実施例１，２のステップＰ８（図２，図８）では条件［Io,θ,r,ch,do,d］毎にｎ次関数の係数ａ０，…，ａｎを求めたが、実施例３では条件［Io,do,d］毎にｎ次関数の係数ａ０，…，ａｎを求め、記憶する。
実施例３によれば、ビュー角度θ毎に且つ検出器det(r,ch,do)毎に係数ａ０，…，ａｎを記憶しないで済み、構成を簡単化できる。
実施例１〜３では、投影長和Ｗのｎ次関数として散乱量Ｓを近似するｎ次関数の係数ａ０，…，ａｎを記憶したが、その代わりに、散乱量Ｓ(Io,θ,r,ch,do,d,pt)そのものを記憶しておいてもよい。あるいは、データＩ(Io,θ,r,ch,do,d,pt)に対する散乱量Ｓ(Io,θ,r,ch,do,d,pt)の割合（つまり、データ中の散乱量の割合）を記憶しておいてもよい。あるいは、「１」からデータＩ(Io,θ,r,ch,do,d,pt)に対する散乱量Ｓ(Io,θ,r,ch,do,d,pt)の割合を引いた値を記憶しておいてもよい。
実施例１〜４では、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」に対応付けて係数ａ０，…，ａｎなどを記憶したが、被検体外物の投影長pro_exと撮影対象の投影長pro_ptとに直接的に対応付けて記憶してもよい。あるいは、投影長和「pro_pt＋pro_ex」を、対応する検出器と同一チャネルで全検出器列について加算した値（チャネル毎の情報となる）や、同一チャネルで同一検出器列を含む複数検出器列について加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列で全チャネルについて加算した値（検出器列毎の情報となる）や、同一検出器列で同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、対応する検出器と同一検出器列を含む複数検出器列および同一チャネルを含む複数チャネルについて加算した値（検出器毎の情報となる）や、全検出器について加算した値（全検出器で一つの情報を共有することになる）に対応付けて記憶してもよい。
本発明の散乱補正方法、散乱測定方法およびＸ線ＣＴ装置は、マルチスライス撮影時における散乱に起因するＣＴ画像の画質の低下を抑制するのに利用できる。
実施例１に係るＸ線ＣＴ装置を示すブロック図である。 実施例１に係る散乱測定処理を示すフロー図である。 散乱測定状況を示す説明図である。 被検体外物の投影長の記憶形式を示す概念図である。 散乱係数の記憶形式を示す概念図である。 被検体撮影状況を示す説明図である。 実施例１に係る散乱補正処理を示すフロー図である。 実施例２に係る散乱測定処理を示すフロー図である。
符号の説明
３ 中央処理装置
５ データ収集バッファ
７ 記憶装置
２４ 多列Ｘ線検出器
２５ ＤＡＳ
１００ Ｘ線ＣＴ装置

Claims (19)

  1. 被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求め、撮影対象の投影長pro_ptを求め、前記被検体外物および前記撮影対象をＸ線通過経路に入れた状態で、検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を測定し、検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し、前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求め、前記被検体外物の投影長pro_exと前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶しておき、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め、前記被検体の投影長pro_ptと前記被検体外物の投影長pro_exに対応する前記情報を読み出し該情報を用いて前記データを散乱補正することを特徴とする散乱補正方法。
  2. 被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を求めるための情報を記憶しておき、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め、前記被検体の投影長pro_ptと前記被検体外物の投影長pro_exに対応する前記情報を読み出し該情報を用いて前記データを散乱補正することを特徴とする散乱補正方法。
  3. 請求項１または請求項２に記載の散乱補正方法において、前記情報が、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)であることを特徴とする散乱補正方法。
  4. 請求項３に記載の散乱補正方法において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報が前記ｎ次関数の係数ａ０，…，ａｎの全部または一部であることを特徴とする散乱補正方法。
  5. 請求項１から請求項４のいずれかに記載の散乱補正方法において、検出器列とチャネルによって特定される検出器毎に前記情報を用いて散乱量を求め、検出器毎の前記データから前記散乱量を減算することにより、前記データを散乱補正することを特徴とする散乱補正方法。
  6. 請求項１から請求項５のいずれかに記載の散乱補正方法において、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについての２種類以上の異なる値について前記情報を記憶しておくことを特徴とする散乱補正方法。
  7. 被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求めるステップと、
    撮影対象の投影長pro_ptを求めるステップと、
    前記被検体外物および前記撮影対象をＸ線通過経路に入れ、検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を測定し、検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し、前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求めるステップと、
    前記被検体外物の投影長pro_exと前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶するステップと
    を有することを特徴とする散乱測定方法。
  8. 請求項７に記載の散乱測定方法において、前記被検体外物の投影長pro_exを求めるステップが、前記被検体外物をＸ線通過経路に入れないで撮影してデータＩ_exoutを測定し、前記被検体外物をＸ線通過経路に入れて撮影してデータＩ_exinを測定し、前記データＩ_exoutおよび前記データＩ_exinより前記被検体外物の投影長pro_exを求めるステップであることを特徴とする散乱測定方法。
  9. 請求項７に記載の散乱測定方法において、前記被検体外物の投影長pro_exを求めるステップが、前記被検体外物の材質と形状より前記被検体外物の投影長pro_exを算出するステップであることを特徴とする散乱測定方法。
  10. 請求項７から請求項９のいずれかに記載の散乱測定方法において、前記撮影対象の投影長pro_ptを求めるステップが、撮影対象をＸ線通過経路に入れないで撮影してデータＩ_ptoutを測定し、前記撮影対象をＸ線通過経路に入れて撮影してデータＩ_ptinを測定し、前記データＩ_ptoutおよび前記データＩ_ptinより前記撮影対象の投影長pro_ptを求めるステップであることを特徴とする散乱測定方法。
  11. 請求項７から請求項１０のいずれかに記載の散乱測定方法において、被検体外物の投影長pro_exと撮影対象の投影長pro_ptの和である投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)を求め、前記情報として記憶することを特徴とする散乱測定方法。
  12. 請求項１１に記載の散乱測定方法において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報として前記ｎ次関数の係数ａ０，…，ａｎの全部または一部を記憶することを特徴とする散乱測定方法。
  13. 請求項７から請求項１２のいずれかに記載の散乱測定方法において、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについて２種類以上の異なる値で測定を行い、それぞれについて前記情報を記憶することを特徴とする散乱測定方法。
  14. Ｘ線管と、多列Ｘ線検出器と、前記Ｘ線管および前記多列Ｘ線検出器の少なくとも一方を被検体の周りに相対回転させながらデータを収集するスキャン手段と、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exを求める被検体外物投影長取得手段と、撮影対象の投影長pro_ptを求める撮影対象投影長取得手段と、前記被検体外物および前記撮影対象をＸ線通過経路に入れて検出器厚ｄｏと同等のビーム厚で撮影してデータＩ(do,do)を収集すると共に検出器厚ｄｏより厚いビーム厚ｄで撮影対象を撮影してデータＩ(do,d)を測定し前記データＩ(do,do)と前記データＩ(do,d)の差に基づいて散乱量Ｓ(do,d)を求める散乱量取得手段と、前記被検体外物の投影長pro_exと前記撮影対象の投影長pro_ptから前記散乱量Ｓ(do,d)を求めるための情報を記憶する情報記憶手段と、ビーム厚ｄのＸ線ビームで且つ検出器厚ｄｏの検出器でデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め両者に対応する情報を読み出して該情報を用いて前記データを散乱補正する散乱補正手段とを具備したことを特徴とするＸ線ＣＴ装置。
  15. Ｘ線管と、多列Ｘ線検出器と、前記Ｘ線管および前記多列Ｘ線検出器の少なくとも一方を被検体の周りに相対回転させながらデータを収集するスキャン手段と、被検体撮影時にＸ線通過経路に存在して無視できないＸ線散乱を生じさせる物体であって且つ被検体以外の物体である被検体外物の投影長pro_exと撮影対象の投影長pro_ptから散乱量Ｓ(do,d)を求めるための情報を記憶する情報記憶手段と、被検体を撮影してデータを収集し前記データから被検体の投影長pro_ptを求めると共に前記データに対応する被検体外物の投影長pro_exを求め両者に対応する前記情報を読み出し該情報を用いて前記データを散乱補正する散乱補正手段とを具備したことを特徴とするＸ線ＣＴ装置。
  16. 請求項１４または請求項１５に記載のＸ線ＣＴ装置において、前記情報が、前記投影長和「pro_pt＋pro_ex」から前記散乱量Ｓ(do,d)を算出するための散乱関数ｆ(pro_pt＋pro_ex)であることを特徴とするＸ線ＣＴ装置。
  17. 請求項１６に記載のＸ線ＣＴ装置において、前記散乱関数ｆ(pro_pt＋pro_ex)が投影長和「pro_pt＋pro_ex」のｎ（≧１）次関数であり、前記情報が前記ｎ次関数の全部または一部であることを特徴とするＸ線ＣＴ装置。
  18. 請求項１４から請求項１７のいずれかに記載のＸ線ＣＴ装置において、前記散乱補正手段は、検出器列とチャネルによって特定される検出器毎に前記情報を用いて散乱量を求め、検出器毎の前記データから前記散乱量を減算することにより、前記データを散乱補正することを特徴とするＸ線ＣＴ装置。
  19. 請求項１４から請求項１８のいずれかに記載のＸ線ＣＴ装置において、前記情報記憶手段は、Ｘ線管出力Ｉｏと検出器厚ｄｏとビーム厚ｄと検出器列ｒとチャネルｃｈのうちの少なくとも１つについての２種類以上の異なる値について前記情報を記憶しておくことを特徴とするＸ線ＣＴ装置。
JP2005056083A 2005-03-01 2005-03-01 Ｘ線ｃｔ装置および散乱測定方法 Active JP4336661B2 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
JP2005056083A JP4336661B2 (ja) 2005-03-01 2005-03-01 Ｘ線ｃｔ装置および散乱測定方法
FR0601817A FR2885007A1 (fr) 2005-03-01 2006-03-01 Procede de compensation de la diffusion et procede de mesure
US11/365,249 US7336759B2 (en) 2005-03-01 2006-03-01 Scattering compensating method and scattering measuring method
CNB2006100820951A CN100563572C (zh) 2005-03-01 2006-03-01 散射补偿方法和散射测量方法
DE200610009762 DE102006009762A1 (de) 2005-03-01 2006-03-01 Streukompensationsverfahren und Streumessverfahren
KR1020060020080A KR20060099412A (ko) 2005-03-01 2006-03-02 산란 보정 방법 및 산란 측정 방법

Applications Claiming Priority (1)

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

Publications (2)

Publication Number Publication Date
JP2006239003A true JP2006239003A (ja) 2006-09-14
JP4336661B2 JP4336661B2 (ja) 2009-09-30

Family

ID=36848338

Family Applications (1)

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

Country Status (6)

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

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008119095A (ja) * 2006-11-09 2008-05-29 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置及び散乱補正方法
JP2009106433A (ja) * 2007-10-29 2009-05-21 Hitachi Medical Corp 放射線撮像装置
JP2011510757A (ja) * 2008-01-30 2011-04-07 ザ リージェンツ オブ ザ ユニヴァーシティー オブ カリフォルニア 対象の周囲を動的な束縛として使用することによる、トモグラフィにおける放射線量低減および画質向上
JP2022509306A (ja) * 2018-11-30 2022-01-20 アキュレイ インコーポレイテッド イメージングにおける散乱評価および散乱補正を向上させるための方法および機器
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

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2276408B1 (en) 2008-05-08 2019-07-10 Arineta Ltd. X ray imaging system with scatter radiation correction and method of using same
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 东软医疗系统股份有限公司 一种散射校正方法、装置及设备
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

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
EP0660599B2 (en) * 1993-12-24 2002-08-14 Agfa-Gevaert Partially-transparent-shield-method for scattered radiation compensation in x-ray imaging
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
US6845147B2 (en) * 2002-06-17 2005-01-18 Edax Inc. Scatter spectra method for x-ray fluorescent analysis with optical components
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

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008119095A (ja) * 2006-11-09 2008-05-29 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置及び散乱補正方法
JP2009106433A (ja) * 2007-10-29 2009-05-21 Hitachi Medical Corp 放射線撮像装置
JP2011510757A (ja) * 2008-01-30 2011-04-07 ザ リージェンツ オブ ザ ユニヴァーシティー オブ カリフォルニア 対象の周囲を動的な束縛として使用することによる、トモグラフィにおける放射線量低減および画質向上
US8532350B2 (en) 2008-01-30 2013-09-10 The Regents Of The University Of California Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints
JP2022509306A (ja) * 2018-11-30 2022-01-20 アキュレイ インコーポレイテッド イメージングにおける散乱評価および散乱補正を向上させるための方法および機器
US11890125B2 (en) 2018-11-30 2024-02-06 Accuray, Inc. Multimodal radiation apparatus and methods
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

Also Published As

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

Similar Documents

Publication Publication Date Title
JP4336661B2 (ja) Ｘ線ｃｔ装置および散乱測定方法
US7391844B2 (en) Method and apparatus for correcting for beam hardening in CT images
US7801264B2 (en) Method for calibrating a dual -spectral computed tomography (CT) system
JP4578675B2 (ja) 投影データを用いて心運動を評価する方法及び装置
JP5336364B2 (ja) 散乱放射線の割合に応じたｘ線検出器の利得校正
US7747057B2 (en) Methods and apparatus for BIS correction
JP4817736B2 (ja) Ｘ線コンピュータ断層撮影システム及びｘ線コンピュータ断層撮影システムのデータ補正方法
US20040228451A1 (en) Method and apparatus for calibrating detector spectral response
JP2007181623A (ja) Ｘ線ｃｔ装置
JP5216680B2 (ja) Ｘ線ｃｔ装置およびこれを用いたデータ取得方法
JP3950855B2 (ja) 散乱測定方法、散乱補正方法およびｘ線ｃｔ装置
WO2005011502A1 (ja) 放射線断層撮影装置
CN107533766B (zh) 用于来自牙齿图像产生系统的图像数据的图像改进方法
US7835485B2 (en) Method for scattered radiation correction in x-ray imaging devices
JP3908993B2 (ja) Ｘ線ｃｔ装置
JP2005143759A (ja) Ｘ線ｃｔ装置
JP2022113115A (ja) ビームハードニング補正方法、ｘ線ｃｔ装置及びビームハードニング補正プログラム
US6411671B2 (en) Method for reducing line artifacts in a CT image and device for implementing the method
US6904120B2 (en) Method and apparatus for correcting bone induced spectral artifacts
JP4397513B2 (ja) Ｘ線ｃｔ装置
JP4703221B2 (ja) Ｘ線ｃｔ装置
JP2007125129A (ja) Ｘ線ｃｔ装置
JP4064541B2 (ja) レファレンス信号生成方法および装置並びに放射線断層撮影装置
JP3857002B2 (ja) Ｘ線断層撮影方法及び装置
JP2006068227A (ja) Ｃｔ画像の再構成方法及びｃｔ装置並びにプログラム

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081118

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090216

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090310

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090513

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20090609

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090629

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120703

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4336661

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120703

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120703

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120703

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120703

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130703

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130703

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130703

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250