JP4668161B2 - 直線軌跡走査を採用した画像復元システム及び方法 - Google Patents

直線軌跡走査を採用した画像復元システム及び方法 Download PDF

Info

Publication number
JP4668161B2
JP4668161B2 JP2006313359A JP2006313359A JP4668161B2 JP 4668161 B2 JP4668161 B2 JP 4668161B2 JP 2006313359 A JP2006313359 A JP 2006313359A JP 2006313359 A JP2006313359 A JP 2006313359A JP 4668161 B2 JP4668161 B2 JP 4668161B2
Authority
JP
Japan
Prior art keywords
projection data
projection
scan
detector
image
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.)
Active
Application number
JP2006313359A
Other languages
English (en)
Other versions
JP2007139783A (ja
Inventor
麗 張
河偉 高
志強 陳
克軍 康
建平 程
元景 李
以農 劉
宇翔 ▲シン▼
自然 趙
永順 肖
Original Assignee
ツィンファ ユニバーシティ
同方威視技術股▲フン▼有限公司
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 ツィンファ ユニバーシティ, 同方威視技術股▲フン▼有限公司 filed Critical ツィンファ ユニバーシティ
Publication of JP2007139783A publication Critical patent/JP2007139783A/ja
Application granted granted Critical
Publication of JP4668161B2 publication Critical patent/JP4668161B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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]
    • A61B6/032—Transmission computed tomography [CT]
    • 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/30—Accessories, mechanical or electrical features
    • G01N2223/33—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts
    • G01N2223/3303—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts object fixed; source and detector move
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/421—Filtered back projection [FBP]

Landscapes

  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Pulmonology (AREA)
  • Algebra (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Mathematical Optimization (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pure & Applied Mathematics (AREA)
  • Image Processing (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)

Description

本発明は、輻射結像技術に関し、具体的には、三次元結像中の走査速度を向上させるための、直線軌跡走査を採用した画像復元システム及び方法に関する。
ＣＴ技術の進展に伴い、走査角度が制限されている場合、及びデータカットの場合にも、ある程度の品質の断層画像を復元することが可能であり、これは不完全走査復元画像の実際の応用を可能にする。理論的には、走査経路が直線である結像システムに対して、直線の長さが無限大であると、断層画像を精細に復元できることが分かっている。走査経路の長さが有限であれば、制限された角度（Ｌｉｍｉｔｅｄ−Ａｎｇｌｅ）のＣＴ走査モードと等価である。そこで、不完全復元アルゴリズムを応用して、直線走査の結像システムで採集したデータを復元することで、断層画像を獲得し、立体結像を実現できる。
実際の安全検査において、高速に移動する物体の立体結像と大型物体の回転が難題となる。これは、従来のＣＴ結像システムにおいて、被検査物体に対して円軌道走査を実行し、又は被検査物体を回転し、又は被検査物体を固定し円軌跡走査で検査しなければならないためである。つまり、走査システムと被検査物体との何れかが回転される。回転しなければならないので、大型の被検査物体、例えば列車やトラックなどのような車両に対しては困難である。なお、回転式走査には大錐角（large cone-angle）についての問題もある。
前記問題に鑑みて、直線軌跡の結像システムであるコンピュータ断層結像（Ｃｏｍｐｕｔｅｄ Ｌａｍｉｎｏｇｒａｐｈｙ）システムが提案されている。当該システムは、放射線放射角が小さく、復元アルゴリズムとしてトモグラフ（tomograph）方式を採用したものであるが、三次元立体結像の能力が不足であるので、現在まで安全検査に適用されていない。他の技術としては、再配列アルゴリズムを採用したものがある。それは、直線走査による投影データを平行束の場合のデータを再配列して、画像の復元を行う。しかし、このアルゴリズムは、空間解像度が低いという問題がある。なぜなら、これは直線軌跡走査データを円軌道平行束走査データに再配列するときに、角度方向と検知器方向での補間がシステム解像度の低下を招いているためである。
本発明は、前記課題に鑑みてなされたものである。本発明は、投影データを平行束に再配列する必要をなくし、データ再配列において角度方向での補間と検知器方向での補間による画像解像度の低下を避けることにより、復元断層の品質を向上させた、直線軌跡走査を採用した画像復元システムと方法を提供することを目的とする。
本発明の一つの形態は、直線軌跡走査を採用した画像復元システムであって、直線軌跡走査での投影データを擬平行束走査での投影データに変換するための投影データ変換部と、予め定めた重畳積分関数核で擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタ部と、フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影部とを含む画像復元システムを提供するものである。
また、本発明の一実施例では、画像復元システムは、検知器ユニットを複数備え、放射線源から放射された放射線が被検査物体を透過した透過信号を受信して透過信号を投影データに変換するための検知器アレーを更に含む。
また、本発明の一実施例では、前記複数の検知器ユニットが等間隔で配列される。
また、本発明の一実施例では、前記投影データ変換部が、投影データｐ（ｌ，ｔ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）を獲得する。ここで、投影データｐ（ｌ，ｔ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の座標位置がｔである位置への投影値を示す。前記フィルタ部は、予め定めた重畳積分関数核でｌ方向に沿って、擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）に対して一次元重畳積分を実行し、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を獲得する。前記逆投影部は、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
また、本発明の一実施例では、前記複数の検知器ユニットは、放射線源に対して等角度に配列される。
また、本発明の一実施例では、前記投影データ変換部が、投影データｐ（ｌ，γ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，γ，ｚ）を獲得する。ここで、投影データｐ（ｌ，γ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の角度位置がγである位置への投影値を示す。前記フィルタ部は、予め定めた重畳積分関数核でｌ方向に沿って、擬平行束走査での投影データｑ（ｌ，γ，ｚ）に対して一次元重畳積分を実行し、フィルタリングされた投影データＱ（ｌ´，γ，ｚ）を獲得する。前記逆投影部は、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，γ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
本発明の他の形態は、直線軌跡走査を採用した画像復元方法であって、直線軌跡走査での投影データを擬平行束走査での投影データに変換するための投影データ変換ステップと、予め定めた重畳積分関数核で擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタステップと、フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影ステップとを含む画像復元方法を提供するものである。
また、本発明の一実施例では、本発明の方法は、検知器ユニットを複数備える検知器アレーを利用して、放射線源から放射された放射線が被検査物体を透過した透過信号を受信し、透過信号を投影データに変換するステップを更に含む。
また、本発明の一実施例では、前記複数の検知器ユニットが等間隔で配列される。
また、本発明の一実施例では、前記投影データ変換ステップは、投影データｐ（ｌ，ｔ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）を獲得する。ここで、投影データｐ（ｌ，ｔ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときの、検知器アレーのｚ層目の座標位置がｔである位置への投影値を示す。前記フィルタステップは、予め定めた重畳積分関数核でｌ方向に沿って、擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）に対して一次元重畳積分を実行し、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を獲得する。前記逆投影ステップは、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
また、本発明の一実施例では、前記複数の検知器ユニットは、放射線源に対して等角度に配列される。
また、本発明の一実施例では、前記投影データ変換ステップは、投影データｐ（ｌ，γ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，γ，ｚ）を獲得する。但し、投影データｐ（ｌ，γ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときの、検知器アレーのｚ層目の角度位置がγである位置への投影値を示す。前記フィルタステップは、予め定めた重畳積分関数核でｌ方向に沿って、擬平行束走査での投影データｑ（ｌ，γ，ｚ）に対して一次元重畳積分を実行し、フィルタリングされた投影データＱ（ｌ´，γ，ｚ）を獲得する。前記逆投影ステップは、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，γ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
本発明は、再配列アルゴリズムに比べて、画像復元解像度が向上し、復元画像に対するデータカットの影響が低減される。本発明は、フィルタリング逆投影の形式を採用しているので、一般的なフィルタリング逆投影の利点、例えば復元速度が速い、容易な並行動作、容易な加速動作などの利点を有する。
これによって、従来の結像システムと比べて、本発明のシステムは、透視画像に加えて、断層画像を獲得することも可能であり、透視結像における物体オーバーラップの問題を解決し、安全検査システムにおける高速立体結像を実現することが出来る。従来のＣＴ結像と比べて、本発明のシステムは、設計が簡単で、検査速度が速く、回転の必要がなく、円軌道錐束ＣＴにおける「大錐角」についての問題がないという利点がある。
以下で、本発明の実施例を図面に参照しつつ詳しく説明する。
図１は、立体結像安全検査における本発明の結像システムの直線走査の実行を示す模式図である。図２は本発明の結像システムの構成を示す模式図である。
図１に示すように、被検査物体は、直線に沿って放射線源Ａと検知器との間を移動し、移動の過程において、放射線源Ａは制御システムの命令に従って放射線を放射し、被検査物体を透過する。検知器は透過信号を受信して、制御システムの制御によって投影データを採取し、投影データをメモリに格納する。
図２に示すようなシステムは、被検査物体１２０のサイズと応用領域に応じて選択される。システムは、放射線源１１０と、移送デバイス１３０と、検知器アレー１４０と、制御・データ信号線１５０と、制御・画像処理装置１６０とを含む。放射線源１１０は、例えばＸ線管、加速器放射線源又は同位体源などの放射線源であり、かつ水平方向の放射角が９０度より大きい角度、例えば９０〜１８０度の範囲内の角度である。移送デバイス１３０は、被検査物体１２０を搭載し、安定して移送して、検査過程において搭載された被検査物体１２０を直線に沿って移動させる。検知器アレー１４０は、エリアアレー検知器又はラインアレー検知器であり、移送デバイス１３０を挟んで放射線源１１０の対面に位置し、物体を透過した放射線を検知するのに充分な大きさを持つ。制御・データ信号線１５０は制御・データ信号を伝送する。制御・画像処理装置１６０は、放射線源１１０と移送デバイス１３０と検知器アレー１４０とに制御・データ信号線１５０を介して接続され、走査過程において移送デバイス１３０が直線に沿って移動するよう制御し、放射線源１１０に対して放射線を放射するように命令し、検知器アレー１４０が透過信号の受信を開始するように制御し、投影データを生成し、獲得した投影データに対して後続処理を行う。
これによって、被検査物体１２０は、図１に示す直線軌跡に従って等速度で移動させられ、検知器アレー１４０は、同期して、等時間間隔でサンプリングして、投影データを獲得する。
図３は、図２に示すような結像システムにおける制御・画像処理装置１６０の機能ブロック図である。図３に示すように、制御・画像処理装置１６０は、メモリ１６１と、入力ユニット１６２と、コントローラ１６３と、内部バス１６４と、画像復元ユニット１６５とを含む。メモリ１６１は、データを格納するための、例えばハードディスクなどのような記憶媒体からなる。入力ユニット１６２は、例えばキーボードなどのようなユーザがパラメータや命令を入力するための入力装置である。コントローラ１６３は、ユーザが入力ユニット１６２を介して命令してから、移送デバイス１３０に指示して被検査物体１２０を直線に沿って等速度で移動させることを開始させ、放射線源１１０と検知器アレー１４０の動作を開始させて、投影データを獲得する。内部バス１６４は、各部を接続し制御信号とデータとを伝送する。画像復元ユニット１６５は、検知器アレー１４０により獲得された投影データを復元する。
次に、図４を用いて、画像復元ユニット１６５において画像復元を行う過程を詳しく述べる。図４は等価検知器と復元物体点との間のＺ方向での幾何学関係を示す模式図である。
被検査物体ｆ（ｒ，φ，ｚ）のある近似推定を
とすると、下記の式
となる。
ここで、
である。
ここで、検知器は等間隔に配列されたものであり、データｐ（ｌ，ｔ，ｚ）は、物体１２０が直線上での座標ｌである位置に移動したときに、検知器アレー１４０のｚ層目の座標位置がｔである位置への投影値を示す。注意すべきことは、ｔ及びｚの何れもが検知器アレー１４０の各検知器ユニットを物体直線運動の中心線上に置換した数値である、ということである。
なお、式（１）〜（４）において、Ｄは放射線源１１０から直線運動中心線までの距離を示し、±ｔｍは被検査物体の移動方向での検知器アレー１４０の最小位置と最大位置を示し、ｈは重畳積分関数核を示し、その理論値は
であり、一般的にはＳ−Ｌフィルタ関数を採用し、当該関数の離散形式は、
となる。
そこで、画像復元ユニット１６５において、投影データ変換部１６５１が投影データｐ（ｌ，ｔ，ｚ）を逆転シフトしてｑ（ｌ，ｔ，ｚ）を獲得する。ここでのｑ（ｌ，ｔ，ｚ）は擬平行束走査での投影データを示す。ここでの「擬平行束走査」とは、各角度での検知器の等価サンプル間隔が異なり走査角度サンプルも不均一な可能性があることを指す。
そして、フィルタ部１６５２は、重畳積分関数核ｈでｌ方向に沿って、擬平行束走査での投影データｑ（ｌ,ｔ,ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を獲得する。
続いて、逆投影部１６５３は、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
ここで、説明すべきことは、逆転・シフトを実行する目的が、直線走査での投影データを、擬平行束走査での投影データに変換することであり、ここでの擬平行束走査は標準のＣＴ復元中の平行束ではない、ということである。これは、各々の走査角度での検知器による等価サンプル間隔が異なり、角度サンプルが不均一であることも可能であるからである。
なお、重畳積分関数核ｈでフィルタリングする目的は、標準ＦＢＰ復元アルゴリズム中のフィルタリングと同様であり、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を重み付け逆投影することで復元画像を獲得できる。
データ採集方向ｌに沿ってフィルタリングし、放射線方向に沿って逆投影され、平行束に再配列するアルゴリズムと比べて、本発明は、各有効データを充分利用し、画像の解像度を向上し、且つデータカットに対する影響を再配列アルゴリズムの場合より小さくすることが出来る。
以下に、図１，４及び５を参照して、前記の式（１）を導出する。導出する前に、まず、直線走査データを円軌道平行束走査に再配列する過程を述べる。
図１に示す走査方式のように、各検出器が一つの角度走査の角度に対応し、物体ｆ（ｘ，ｙ）の移動過程における、当該角度での平行束走査に相当する。図５の投影模式図を参照すれば、等間隔に配列された検知器アレーに対しては、直線走査データを円軌道平行束走査に再配列する再配列数式は、
となる。ここで、
となり、円軌道平行束走査における、走査角度がθで、回転中心からの距離がｓである位置への投影データを示している。ｐ（ｌ，ｔ）は、アレー検知器は物体が直線上の座標ｌである位置に移動したときの、検知器アレー座標位置がｔである位置への投影値を示している。
数式（６）を利用して、直線軌跡走査投影データを円軌道平行束走査での投影データに再配列することが実現できる。しかし、実際のシステムにおいては直線の長さが無限大であることはあり得ないので、再配列されたデータも、円軌道上で１８０度分の平行束走査データとはなり得ない。つまり、ここでのＣＴ復元について言えば、データは不完全であると言える。
直線走査においては、ｌとｔのサンプリングは均一なものであるが、円軌道平行束走査における角度θと検知器位置ｓに対するサンプリングの何れもが不均一なものである。したがって、再配列のためには、角度方向及び検知器方向において補間を行うことが必要であり、復元画像の解像度は低下してしまう。
次に、本発明による直線走査データを直接フィルタリング逆投影して画像を復元する過程を詳しく説明する。
円軌道平行束走査でのフィルタリング逆投影復元数式は、
である。無限長の直線軌跡、等間隔検知器に対して、数式（７）を利用し、パラメータ（ｌ，ｔ）で（θ，ｓ）を置換すると、
となる。ここで、
である。以下のように証明する。
ここで、
である。
直線軌跡走査では、ｐ（ｌ，ｔ）で
を置換する。同時に、図５の幾何学構造に基づいて、
となる。ここで、
であり、これは直線走査での空間サンプル位置における点（ｒ，φ）及び第ｔ検知器を経過した投影データを示している。
（１０）を（９）に代入し、
を利用し、
ｌ＝ｌ−ｔ, ｌ´＝ｌ´−ｔを代入して、直ちに復元数式（８）を得る。
数式（８）に対しては、ｑ(ｌ,ｔ)＝ｐ(−ｌ＋ｔ,ｔ)を数式（８）に代入すれば、
となる。ここで、Ｑ（ｌ´，ｔ）＝ｑ（ｌ，ｔ）＊ｈ（ｌ）である。
実際には、直線軌跡が［−Ｌ，Ｌ］であり、検知器の総長が［−ｔｍ，ｔｍ］であれば、数式（８）に従って復元される画像は精確なｆ（ｘ，ｙ）ではなく、ある近似したものである。さらに、三次元の場合を考えると、被検査物体ｆ（ｘ，ｙ，ｚ）のある近似したものである
は式（１）に示すもので表わすことができる。
図６に本発明による直線フィルタリング逆投影方法のフローチャートを示す。図６に示すように、検知器アレー１４０が投影データを獲得しメモリ１６１に格納してから、画像を復元しようとするときに、ステップＳ１１０にて、メモリから直線走査で獲得された投影データｐ（ｌ，ｔ，ｚ）を読み取る。
そして、ステップＳ１２０において、投影データ変換部１６５１で投影データｐ（ｌ，ｔ，ｚ）を逆転・シフトしてｑ（ｌ，ｔ，ｚ）を取得する。ここでのｑ（ｌ，ｔ，ｚ）は、擬平行束走査での投影データを示している。
そして、ステップＳ１３０にて、フィルタ部１６５２は、重畳積分関数核ｈでｌ方向に沿って擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）に対して一次元重畳積分を実行し、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を取得する。
次に、ステップ１４０にて、逆投影部１６５３は、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
ステップＳ１５０にて、復元した画像をメモリ１６１に格納し、又はディスプレー１７０のスクリーンに表示する。
等間隔検知器アレーの場合、直線軌跡走査のフィルタリング逆投影復元数式（１）及び本発明による方法の詳細の実行過程を前記のように導出する。実際には、検知器アレー１４０における検知器ユニット同士は、放射線源１１０に対して等角度の方式に従って配列されることも可能である。検知器が等角度に配列されれば、前記の導出過程と類似し、そのフィルタリング逆投影復元数式は、
となる。ここで、
である。ここで、検知器ユニットは等角度に配列されたものであり、データｐ（ｌ，γ，ｚ）は、物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の角度位置がγである位置への投影値を表している。注意すべきことは、γ及びｚの何れもが、検知器アレーを物体の直線運動の中心線に変換した数値である、ということである。±γｍは、被検査物体の移動方向での検知器アレーの最小の角度と最大の角度を表している。
したがって、等角度に配列された検知器の場合には、直線フィルタリング逆投影の復元過程は上述のように行われ、そこでは、逆転・シフト操作が数式（１４）に従って行われ、重畳積分の意味は等間隔の場合のそれと同じである。
言い換えれば、投影データ変換部１６５１において、投影データｐ（ｌ，γ，ｚ）を逆転・シフトすることで擬平行束走査での投影データｑ（ｌ，γ，ｚ）を獲得する。ここで、投影データｐ（ｌ，γ，ｚ）は、物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の角度位置がγである位置への投影値を表している。
フィルタ部１６５２において、予め定めた重畳積分関数核ｈでｌ方向に沿って、擬平行束走査での投影データｑ（ｌ，γ，ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，γ，ｚ）を獲得する。
逆投影部１６５３において、放射線方向に沿ってフィルタリングされた投影データＱ（ｌ´，γ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得する。
精確な画像復元を実現するために、Ｘ線結像システムは以下のようなシステムパラメータを精確に検知／設定できるべきである。当該システムパラメータは、Ｘ線源点から検知器までの距離Ｔと、Ｘ線源点から直線運動中心までの距離Ｄと、移送装置の直線移動速度νと、アレー検知器のサンプリング間隔Δｔと、検知器の物理サイズとを含む。但し、検知器の物理サイズは、単一検知器ユニットの物理サイズと検知器アレーの物理サイズなどを含む。
図７に、本発明の方法を採用した場合と再配列平行束アルゴリズムを採用した場合とで、結像システムで採集したデータを復元して獲得した断層画像（Ｘ−Ｙ平面）を比較した結果を示す。ここで、模擬放射線放射角は１５０度であり、検知器の幅は６ｍｍで等間隔に配列され、検知器アレー等価空間サンプル間隔は３ｍｍである。復元画像は３００＊３００画素であり、画素毎のサイズは３ｍｍ＊３ｍｍである。同図は、荷物モデルを利用して行われた模擬実験の結果であり、復元されたものは中心層の画像である。但し、（Ａ）は原始模擬画像を、（Ｂ）は円軌道平行束走査で精確に復元された画像を、（Ｃ）は直線軌跡走査で獲得された投影データを使用し、平行束に再配列する方法によって復元された画像を、（Ｄ）は本発明の方法によって復元された画像を、（Ｅ）は、データを検知器方向で切断した場合の、平行束に再配列する方法によって復元されたもう一つの画像を、（Ｆ）は、データを検知器方向で切断した場合の、本発明の方法によって復元されたもう一つの画像を、それぞれ表している。図７の比較から分かるように、本発明の方法では、平行束再配列方法に比べて、解像度が大きく向上し、データカットによる影響が大きく低減している。
これまで述べてきたものは、本発明における具体的な実施の形態に過ぎず、本発明の保護範囲はこれに限定されるものではなく、本発明中に示唆された技術範囲内に入る当業者が容易に想到可能な変換や置換もまた、本発明の保護範囲に含まれる。したがって、本発明の保護範囲は、クレームにおける保護範囲を基準として定められるべきである。
本発明による結像システムにおける直線軌跡走査の実行を示す平面模式図である。 本発明の実施例による結像システムの構成を示す模式図である。 図２に示す結像システムにおける制御・画像処理装置の機能ブロック図である。 等価検知器と復元物体点との間のＺ方向での幾何学関係を示す模式図である。 本発明の一実施例による直線フィルタリング逆投影方法の幾何学関係を説明するための模式図である。 本発明による直線フィルタリング逆投影方法を示すフローチャートである。 本発明の方法復元と再配列平行束アルゴリズムとを使用し、立体結像システムで採取したデータを復元することで獲得した模擬画像（Ｘ−Ｙ平面）の効果の比較を示したものである。
符号の説明
１６１ メモリ
１６２ 入力ユニット
１６３ コントローラ
１６４ 内部バス
１６５ 画像復元ユニット
１６５１ 投影データ変換部
１６５２ フィルタ部
１６５３ 逆投影部

Claims (8)

  1. 回転動作を含まない直線軌跡走査を採用した画像復元システムであって、
    直線軌跡走査での投影データを、各走査角度での検知器の等価サンプル間隔が異なり走査角度サンプルも不均一でありうる擬平行束走査での投影データに変換するための投影データ変換部と、
    予め定めた重畳積分関数核で前記擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタ部と、
    フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影部と、
    検知器ユニットを複数備え、放射線源から放射された放射線が被検査物体を透過した透過信号を受信して透過信号を投影データに変換するための検知器アレーと、
    を含み、
    前記投影データ変換部が、投影データｐ（ｌ，ｔ，ｚ）を逆転・シフトして前記擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）を獲得し、但し、前記投影データｐ（ｌ，ｔ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、前記検知器アレーのｚ層目の座標位置がｔである位置への投影値を示し、
    前記フィルタ部が、予め定めた重畳積分関数核でｌ方向に沿って、前記擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を獲得し、
    前記逆投影部が、放射線方向に沿って前記フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得し、
    前記逆転・シフトは、下記式に従って行なわれ、
    前記等価サンプル間隔とは、走査角度ごとに被検査物体が二回の走査の間に移動する距離を、放射線源と検知器ユニットとの接続線に垂直する方向への距離に置換した値である、ことを特徴とする画像復元システム。
  2. 前記複数の検知器ユニットが等間隔で配列されることを特徴とする請求項１に記載する画像復元システム。
  3. 回転動作を含まない直線軌跡走査を採用した画像復元システムであって、
    直線軌跡走査での投影データを、各走査角度での検知器の等価サンプル間隔が異なり走査角度サンプルも不均一でありうる擬平行束走査での投影データに変換するための投影データ変換部と、
    予め定めた重畳積分関数核で前記擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタ部と、
    フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影部と、
    検知器ユニットを複数備え、放射線源から放射された放射線が被検査物体を透過した透過信号を受信して透過信号を投影データに変換するための検知器アレーと、
    を含み、
    前記投影データ変換部が、投影データｐ（ｌ，γ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，γ，ｚ）を獲得し、但し、前記投影データｐ（ｌ，γ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の角度位置がγである位置への投影値を示し、
    前記フィルタ部が、予め定めた重畳積分関数核でｌ方向に沿って、前記擬平行束走査での投影データｑ（ｌ，γ，ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，γ，ｚ）を獲得し、
    前記逆投影部が、放射線方向に沿って前記フィルタリングされた投影データＱ（ｌ´，γ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得し、
    前記逆転・シフトは、下記式に従って行なわれ、
    前記等価サンプル間隔とは、走査角度ごとに被検査物体が二回の走査の間に移動する距離を、放射線源と検知器ユニットとの接続線に垂直する方向への距離に置換した値である、ことを特徴とする画像復元システム。
  4. 前記複数の検知器ユニットは、放射線源に対して等間隔で配列されることを特徴とする請求項３に記載する画像復元システム。
  5. 回転動作を含まない直線軌跡走査を採用した画像復元方法であって、
    直線軌跡走査での投影データを、各走査角度での検知器の等価サンプル間隔が異なり走査角度サンプルも不均一でありうる擬平行束走査での投影データに変換するための投影データ変換ステップと、
    予め定めた重畳積分関数核で前記擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタステップと、
    フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影ステップと、
    検知器ユニットを複数備える検知器アレーを利用して、放射線源から放射された放射線が被検査物体を透過した透過信号を受信し、透過信号を投影データに変換するステップとを含み、
    前記投影データ変換ステップが、投影データｐ（ｌ，ｔ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）を獲得し、但し、前記投影データｐ（ｌ，ｔ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の座標位置がｔである位置への投影値を示し、
    前記フィルタステップが、予め定めた重畳積分関数核でｌ方向に沿って、前記擬平行束走査での投影データｑ（ｌ，ｔ，ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）を獲得し、
    前記逆投影ステップが、放射線方向に沿って前記フィルタリングされた投影データＱ（ｌ´，ｔ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得し、
    前記逆転・シフトは、下記式に従って行なわれ、
    前記等価サンプル間隔とは、走査角度ごとに被検査物体が二回の走査の間に移動する距離を、放射線源と検知器ユニットとの接続線に垂直する方向への距離に置換した値である、ことを特徴とする画像復元方法。
  6. 前記複数の検知器ユニットが等間隔で配列されることを特徴とする請求項５に記載する画像復元方法。
  7. 回転動作を含まない直線軌跡走査を採用した画像復元方法であって、
    直線軌跡走査での投影データを、各走査角度での検知器の等価サンプル間隔が異なり走査角度サンプルも不均一でありうる擬平行束走査での投影データに変換するための投影データ変換ステップと、
    予め定めた重畳積分関数核で前記擬平行束走査での投影データを重畳積分することで、フィルタリングされた投影データを獲得するためのフィルタステップと、
    フィルタリングされた投影データに対して重み付け逆投影を実行することで画像を復元する逆投影ステップと、
    検知器ユニットを複数備える検知器アレーを利用して、放射線源から放射された放射線が被検査物体を透過した透過信号を受信し、透過信号を投影データに変換するステップとを含み、
    前記投影データ変換ステップが、投影データｐ（ｌ，γ，ｚ）を逆転・シフトして擬平行束走査での投影データｑ（ｌ，γ，ｚ）を獲得し、但し、前記投影データｐ（ｌ，γ，ｚ）は、被検査物体が直線上の座標ｌである位置に移動したときに、検知器アレーのｚ層目の角度位置がγである位置への投影値を示し、
    前記フィルタステップが、予め定めた重畳積分関数核でｌ方向に沿って、前記擬平行束走査での投影データｑ（ｌ，γ，ｚ）に対して一次元重畳積分を実行して、フィルタリングされた投影データＱ（ｌ´，γ，ｚ）を獲得し、
    前記逆投影ステップが、放射線方向に沿って前記フィルタリングされた投影データＱ（ｌ´，γ，ｚ）に対して重み付け逆投影操作を実行することで、復元画像を獲得し、
    前記逆転・シフトは、下記式に従って行なわれ、
    前記等価サンプル間隔とは、走査角度ごとに被検査物体が二回の走査の間に移動する距離を、放射線源と検知器ユニットとの接続線に垂直する方向への距離に置換した値である、ことを特徴とする画像復元方法。
  8. 前記複数の検知器ユニットは、放射線源に対して等角度に配列されることを特徴とする請求項７に記載する画像復元方法。
JP2006313359A 2005-11-21 2006-11-20 直線軌跡走査を採用した画像復元システム及び方法 Active JP4668161B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005101235880A CN100495439C (zh) 2005-11-21 2005-11-21 采用直线轨迹扫描的图像重建系统和方法

Publications (2)

Publication Number Publication Date
JP2007139783A JP2007139783A (ja) 2007-06-07
JP4668161B2 true JP4668161B2 (ja) 2011-04-13

Family

ID=38048282

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006313359A Active JP4668161B2 (ja) 2005-11-21 2006-11-20 直線軌跡走査を採用した画像復元システム及び方法

Country Status (10)

Country Link
US (1) US7424089B2 (ja)
EP (1) EP1953700B1 (ja)
JP (1) JP4668161B2 (ja)
CN (1) CN100495439C (ja)
AT (1) ATE556393T1 (ja)
AU (1) AU2006315013B2 (ja)
DE (1) DE102006054591A1 (ja)
PL (1) PL1953700T3 (ja)
RU (1) RU2334973C1 (ja)
WO (1) WO2007056922A1 (ja)

Families Citing this family (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100565336C (zh) * 2005-11-21 2009-12-02 清华大学 成像系统
CN101681521B (zh) * 2007-06-11 2012-06-27 皇家飞利浦电子股份有限公司 用于对感兴趣区域成像的成像系统和成像方法
CN101561405B (zh) 2008-04-17 2011-07-06 清华大学 一种直线轨迹扫描成像系统和方法
EP2130494A1 (en) * 2008-06-06 2009-12-09 Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt (GmbH) Scanner device and method for computed tomography imaging
DE102009020400B4 (de) * 2009-05-08 2016-04-21 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Bildbestimmung aus beim Durchlaufen einer Trajektorie aufgenommenen Röntgenprojektionen
CN101900694B (zh) * 2009-05-27 2012-05-30 清华大学 基于直线轨迹扫描的双能欠采样物质识别系统和方法
CN102590234B (zh) * 2009-05-27 2014-04-09 清华大学 基于直线轨迹扫描的双能欠采样物质识别系统和方法
EP3569289B1 (en) 2010-02-24 2020-12-09 Accuray, Inc. Gantry image guided radiotherapy system and related target tracking methods
US9687200B2 (en) 2010-06-08 2017-06-27 Accuray Incorporated Radiation treatment delivery system with translatable ring gantry
WO2011126283A2 (en) 2010-04-05 2011-10-13 Samsung Electronics Co., Ltd. Method and apparatus for encoding video based on internal bit depth increment, and method and apparatus for decoding video based on internal bit depth increment
WO2011126282A2 (en) 2010-04-05 2011-10-13 Samsung Electronics Co., Ltd. Method and apparatus for encoding video by using transformation index, and method and apparatus for decoding video by using transformation index
US9369736B2 (en) 2010-04-05 2016-06-14 Samsung Electronics Co., Ltd. Low complexity entropy-encoding/decoding method and apparatus
KR20110112168A (ko) * 2010-04-05 2011-10-12 삼성전자주식회사 내부 비트뎁스 확장에 기반한 비디오 부호화 방법 및 그 장치, 내부 비트뎁스 확장에 기반한 비디오 복호화 방법 및 그 장치
US8559596B2 (en) 2010-06-08 2013-10-15 Accuray Incorporated Target Tracking for image-guided radiation treatment
CN102279970B (zh) * 2010-06-13 2013-02-27 清华大学 基于gpu的螺旋锥束ct重建方法
CN102004111B (zh) * 2010-09-28 2012-09-19 北京航空航天大学 一种倾斜多锥束直线轨迹ct成像方法
JP6214128B2 (ja) * 2010-11-22 2017-10-18 キヤノン株式会社 画像処理装置、画像処理方法、及び記憶媒体
JP5913351B2 (ja) * 2010-11-25 2016-04-27 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 順投影装置
US8536547B2 (en) 2011-01-20 2013-09-17 Accuray Incorporated Ring gantry radiation treatment delivery system with dynamically controllable inward extension of treatment head
EP2745272A1 (en) * 2011-09-12 2014-06-25 Koninklijke Philips N.V. Method of image reconstruction for a filtered back projection in limited angle tomography
CN104144650B (zh) * 2012-03-27 2016-08-31 株式会社日立制作所 放射线拍摄装置和图像处理方法
US9959631B2 (en) 2014-02-21 2018-05-01 Samsung Electronics Co., Ltd. Tomography apparatus and method for reconstructing tomography image thereof
US9865066B2 (en) 2014-05-06 2018-01-09 Astrophysics Inc. Computed tomography system for cargo and transported containers
CN105094725B (zh) * 2014-05-14 2019-02-19 同方威视技术股份有限公司 图像显示方法
CN104142353B (zh) * 2014-06-06 2017-01-11 清华大学 可见光成像辅助定位局部ct扫描方法及系统
CN105785462B (zh) * 2014-06-25 2019-02-22 同方威视技术股份有限公司 一种定位三维ct图像中的目标的方法和安检ct系统
US10353191B2 (en) 2015-01-13 2019-07-16 Massachusetts Institute Of Technology Circular scanning technique for large area inspection
WO2018057517A1 (en) 2016-09-20 2018-03-29 Massachusetts Institute Of Technology Circular scanning technique for large area imaging
CN106934110B (zh) * 2016-12-14 2021-02-26 北京信息科技大学 一种由聚焦堆栈重建光场的反投影方法和装置
CN108242066B (zh) * 2016-12-26 2023-04-14 通用电气公司 Ct图像的空间分辨率增强装置和方法以及ct成像系统
CN106709915B (zh) * 2017-01-06 2020-02-07 天津大学 一种图像重采样操作检测方法
CN106845540B (zh) * 2017-01-06 2020-02-07 天津大学 一种图像重采样操作插值类型识别方法
US11175242B2 (en) * 2018-03-20 2021-11-16 Fei Company Geometric alignment, sample motion correction, and intensity normalization of computed tomography projections using pi-line optimization
CN109685867B (zh) * 2018-12-29 2023-05-30 上海联影医疗科技股份有限公司 一种ct系统、图像重建方法、装置、设备及介质
CN112274169A (zh) * 2020-09-18 2021-01-29 南昌大学 基于直线轨迹投影数据的pet成像系统及方法
CN112504441B (zh) * 2020-12-15 2022-12-13 西安热工研究院有限公司 一种基于重要信息重构的振动加速度信号分段积分方法
CN112649452A (zh) * 2020-12-17 2021-04-13 苏州斯玛维科技有限公司 工业x射线系统成像方法和装置
CN113303820B (zh) * 2021-05-31 2022-06-24 南京安科医疗科技有限公司 紧凑型ct探测器图像重建方法、重建装置和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002267622A (ja) * 2001-03-12 2002-09-18 Shimadzu Corp 断層撮影装置
JP2005265618A (ja) * 2004-03-18 2005-09-29 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2006071472A (ja) * 2004-09-02 2006-03-16 Yukihiro Nishikawa Ｃｔ法およびｃｔ装置

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5053042A (en) * 1990-01-16 1991-10-01 Bidwell Clifford D Biopsy needle guide for use with CT scanner
US5170439A (en) * 1991-06-11 1992-12-08 Picker International, Inc. Cone beam reconstruction using combined circle and line orbits
US6560308B1 (en) * 2001-10-26 2003-05-06 Kabushiki Kaisha Toshiba Method and system for approximating missing data in cone beam x-ray CT reconstruction
JP3950782B2 (ja) * 2002-05-22 2007-08-01 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー ３次元逆投影方法およびｘ線ｃｔ装置
US6775347B2 (en) * 2002-05-29 2004-08-10 Ge Medical Systems Global Technology Company Methods and apparatus for reconstructing an image of an object
US7245755B1 (en) * 2002-07-10 2007-07-17 Xiaochuan Pan Algorithm for image reconstruction and image noise analysis in computed tomography
DE10252661A1 (de) * 2002-11-11 2004-05-27 Philips Intellectual Property & Standards Gmbh Exaktes Computertomographie-Verfahren und Computertomograph
US6862337B2 (en) * 2003-06-25 2005-03-01 General Electric Company Linear track based digital tomosynthesis system and method
US7444011B2 (en) * 2004-02-10 2008-10-28 University Of Chicago Imaging system performing substantially exact reconstruction and using non-traditional trajectories
CN100522063C (zh) * 2004-03-19 2009-08-05 深圳安科高技术股份有限公司 一种双排或多排螺旋ct中的图象重建方法
WO2006073584A2 (en) * 2004-11-24 2006-07-13 Wisconsin Alumni Research Foundation Cone-beam filtered backprojection image reconstruction method for short trajectories
CN100565336C (zh) * 2005-11-21 2009-12-02 清华大学 成像系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002267622A (ja) * 2001-03-12 2002-09-18 Shimadzu Corp 断層撮影装置
JP2005265618A (ja) * 2004-03-18 2005-09-29 Toshiba It & Control Systems Corp コンピュータ断層撮影装置
JP2006071472A (ja) * 2004-09-02 2006-03-16 Yukihiro Nishikawa Ｃｔ法およびｃｔ装置

Also Published As

Publication number Publication date
US20070116175A1 (en) 2007-05-24
RU2334973C1 (ru) 2008-09-27
EP1953700A1 (en) 2008-08-06
JP2007139783A (ja) 2007-06-07
RU2006140975A (ru) 2008-05-27
AU2006315013B2 (en) 2011-03-10
PL1953700T3 (pl) 2013-01-31
EP1953700A4 (en) 2009-08-12
CN100495439C (zh) 2009-06-03
AU2006315013A1 (en) 2007-05-24
WO2007056922A1 (fr) 2007-05-24
CN1971620A (zh) 2007-05-30
ATE556393T1 (de) 2012-05-15
EP1953700B1 (en) 2012-05-02
US7424089B2 (en) 2008-09-09
DE102006054591A1 (de) 2007-06-14

Similar Documents

Publication Publication Date Title
JP4668161B2 (ja) 直線軌跡走査を採用した画像復元システム及び方法
JP4675869B2 (ja) 結像システム
US5598453A (en) Method for X-ray fluoroscopy or radiography, and X-ray apparatus
US6292530B1 (en) Method and apparatus for reconstructing image data acquired by a tomosynthesis x-ray imaging system
US6014419A (en) CT cone beam scanner with fast and complete data acquistion and accurate and efficient regional reconstruction
US6256370B1 (en) Method and apparatus for performing tomosynthesis
US8284892B2 (en) System and method for image reconstruction
JP2005504571A (ja) 多機能コーンビーム結像装置とその方法
US7602879B2 (en) Method for increasing the resolution of a CT image during image reconstruction
JP2004160218A (ja) Ｘ線コンピュータ断層装置、ｘ線コンピュータ断層装置制御方法、及びｘ線コンピュータ断層撮影プログラム
JP4106251B2 (ja) ３次元逆投影方法およびｘ線ｃｔ装置
US5611026A (en) Combining a priori data with partial scan data to project three dimensional imaging of arbitrary objects with computerized tomography
JP2005137389A (ja) Ｘ線ｃｔ撮像方法およびｘ線ｃｔ装置
JP2007159878A (ja) Ｘ線ｃｔ装置およびそのｘ線ｃｔ画像再構成方法
JP4201686B2 (ja) Ｘ線ｃｔ装置
JP3913798B2 (ja) 画像再構成処理装置
JP3921971B2 (ja) 断層再構成方法と断層撮影装置
JP2007198866A (ja) 広義サドルコーンビームｃｔ装置および３次元再構成法
CN2919379Y (zh) 一种采用直线轨迹扫描的图像重建装置
JP3825492B2 (ja) 画像再構成処理装置及びｘ線ｃｔ装置
JP3802650B2 (ja) Ｘ線ｃｔ装置
JPH1043178A (ja) 再構成面設定方法、再構成面設定装置、再構成画像作成方法およびｘ線ｃｔ装置
JPH11347027A (ja) Ｘ線コンピュータ断層撮影装置
JPH08112276A (ja) Ｘ線ｃｔ装置
JP2002204796A (ja) ３次元ｘ線ｃｔ装置

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090812

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090825

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20100223

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100428

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20100729

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101102

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101116

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: 20101214

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: 20110112

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

Free format text: PAYMENT UNTIL: 20140121

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4668161

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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