JP2007236942A - 画像再構成方法およびｘ線ｃｔ装置 - Google Patents

画像再構成方法およびｘ線ｃｔ装置 Download PDF

Info

Publication number
JP2007236942A
JP2007236942A JP2007054179A JP2007054179A JP2007236942A JP 2007236942 A JP2007236942 A JP 2007236942A JP 2007054179 A JP2007054179 A JP 2007054179A JP 2007054179 A JP2007054179 A JP 2007054179A JP 2007236942 A JP2007236942 A JP 2007236942A
Authority
JP
Japan
Prior art keywords
ray
image reconstruction
image
system matrix
reconstruction method
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
JP2007054179A
Other languages
English (en)
Other versions
JP5162143B2 (ja
Inventor
Yu Zhou
チョウ・ユー
Wu Li
リ・ウー
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 JP2007236942A publication Critical patent/JP2007236942A/ja
Application granted granted Critical
Publication of JP5162143B2 publication Critical patent/JP5162143B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/04—Positioning of patients; Tiltable beds or the like
    • 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/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/424—Iterative
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Optics & Photonics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Algebra (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Biophysics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

【課題】画像再構成時間が短く、メモリの消費量が小さく、画像の品質が良い、イテラティブ・リコンストラクションによる画像再構成方法、および、そのような画像再構成を行うＸ線ＣＴ装置を実現する。
【解決手段】それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成し、前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成し、前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成し、前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行い、得られた画像を矩形ピクセルの画像に変換する。
【選択図】図６

Description

本発明は、画像再構成方法およびＸ線ＣＴ(Computed Tomography)装置に関し、特に、イテラティブ・リコンストラクション(iterative reconstruction)による画像再構成方法、および、そのような画像再構成を行うＸ線ＣＴ装置に関する。
Ｘ線ＣＴ装置やＰＥＴ(Positron Emission Tomography)装置等によって得られた投影データ(data)から画像を再構成する技法の一種として、イテラティブ・リコンストラクションが知られている。
これは、システムモデル(system model)上の係数の推定分布から作られる投影データを実測の投影データと比較して、両者が一致するまで係数の推定分布を逐次修正するもので、スタティスティカル・リコンストラクション(statistical reconstruction)とも呼ばれる。そのような技法の代表例として、オーダード・サブセッツ・エクスペクテーション・マキシマイゼーション(ordered subsets Expectation Maximization: OSEM)がある。（例えば、非特許文献１参照）。
エイ ラーミン、ティー ジェー ラス および ヴイ ソッシ（A, Rahmin, T. J. Ruth, and V. Sossi)、スタデイ オブ ア コンバージェント サブセッタイズド リストモード イー エム リコンストラクション アルゴリズム（Study of a Convergent Subsetized List-mode EM Reconstruction Algorithm)、「ニュークリア サイエンス シンポジューム コンファレンス レコード、 ２００４ アイ イー イー イー(Nuclear Science Symposium Conference Record, 2004 IEEE)」、２００４年１０月、第６巻、p. 3978-3982
上記の技法においては、システムモデルとして、係数の空間的配列が直角格子状になっているものが採用される。そのようなシステムモデルに対する画像再構成用のシステムマトリクス(system matrix)は、構造が複雑でデータ量が膨大なものとなるので、その計算はコンピュータ(computer)にとって大きな負荷となる。
このため、図１５の（ａ）に示すように、予め計算したストアド・マトリクス(stored matrix)を使用するときはメモリ(memory)の消費量が膨大になり、（ｂ）に示すように 、システムマトリクスをそのつど計算しながら画像再構成を行うときは画像再構成時間 が長くなる。また、画像再構成時間を短縮するためにシステムマトリクスを単純化する と画像の品質が低下する。
そこで、本発明の課題は、画像再構成時間が短く、メモリの消費量が小さく、画像の品質が良い、イテラティブ・リコンストラクションによる画像再構成方法、および、そのような画像再構成を行うＸ線ＣＴ装置を実現することである。
上記の課題を解決するための第１の観点での発明は、それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成し、前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成し、前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成し、前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行い、得られた画像を矩形ピクセルの画像に変換する、ことを特徴とする画像再構成方法である。
上記の課題を解決するための第２の観点での発明は、被検体をＸ線でスキャンして投影データを収集するデータ収集装置と、収集された投影データに基づいて画像を再構成する画像再構成装置とを有するＸ線ＣＴ装置であって、前記画像再構成装置は、それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成するシステムモデル構成手段と、前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成するシステムマトリクス構成手段と、前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成するシステムマトリクス作成手段と、前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行う画像再構成手段と、
得られた画像を矩形ピクセルの画像に変換する変換手段と、を具備することを特徴とするＸ線ＣＴ装置である。
前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影線の長さであることが、システムマトリクスを簡易化する点で好ましい。
前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影ストライプの部分的エリアであることが、システムマトリクスを精密化する点で好ましい。
前記部分的エリアは正規化されたエリアであることが、エレメントを規格化する点で好ましい。
前記他のビューのシステムマトリクスの作成をイテレーションごとに行うことが、システムマトリクスを適時に作成する点で好ましい。
前記イテラティブ・リコンストラクションはオーダード・サブセッツ・エクスペクテーション・マキシマイゼーション(OSEM)であることが、画像再構成を適切に行う点で好ましい。
前記矩形ピクセルの画像への変換をバイリニア・インターポレーションによって行うことが、変換を適切に行う点で好ましい。
前記投影線はファンビームＸ線に対応する投影線であることが、ファンビームＸ線による投影データに適応する点で好ましい。
前記投影線はパラレルビームＸ線に対応する投影線であることが、パラレルビームＸ線による投影データに適応する点で好ましい。
本発明によれば、それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成し、前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成し、前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成し、前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行い、得られた画像を矩形ピクセルの画像に変換するので、画像再構成時間が短く、メモリの消費量が小さく、画像の品質が良い、イテラティブ・リコンストラクションによる画像再構成方法、および、そのような画像再構成を行うｘＸ線ＣＴ装置を実現することができる。
以下、図面を参照して発明を実施するための最良の形態を説明する。なお、本発明は、発明を実施するための最良の形態に限定されるものではない。図１にＸ線ＣＴ装置の模式的構成を示す。本装置は本発明を実施するための最良の形態の一例である。本装置の構成によって、Ｘ線ＣＴ装置に関する発明を実施するための最良の形態の一例が示される。
本装置は、ガントリ(gantry)１００、テーブル(table)２００およびオペレータコンソール(operator console)３００を有する。ガントリ１００は、テーブル２００によって搬入される被検体１０を、Ｘ線照射・検出装置１１０でスキャンして複数ビューの透過Ｘ線信号（投影データ）を収集し、オペレータコンソール３００に入力する。ガントリ１００は、本発明におけるデータ収集装置の一例である。
オペレータコンソール３００は、内蔵のコンピュータにより、ガントリ１００から入力された投影データに基づいて画像再構成を行い、再構成した画像をディスプレイ(display)３０２に表示する。オペレータコンソール３００は、本発明における画像再構成装置の一例である。
オペレータコンソール３００は、ガントリ１００とテーブル２００の動作を制御する。オペレータコンソール３００による制御の下で、ガントリ１００は所定のスキャン条件でスキャンを行い、テーブル２００は所定の部位がスキャンされるように、被検体１０の位置決めを行う。位置決めは、内蔵する位置調節機構により、天板２０２の高さおよび天板上のクレードル(cradle)２０４の水平移動距離を調節することによって行われる。
クレードル２０４を停止させた状態でスキャンすることにより、アキシャルスキャン(axial scan)を行うことができる。クレードル２０４を連続的に移動させながら複数回のスキャンを連続的に行うことにより、ヘリカルスキャン(helical scan)を行うことができる。クレードル２０４を間欠的に移動させながら停止位置ごとにスキャンすることによりクラスタスキャン(cluster scan)を行うことができる。
天板２０２の高さ調節は、支柱２０６をベース(base)２０８への取り付け部を中心としてスイング(swing)させることによって行われる。支柱２０６のスイングによって、天板２０２は垂直方向および水平方向に変位する。クレードル２０４は天板２０２上で水平方向に移動して天板２０２の水平方向の変位を相殺する。スキャン条件によっては、ガントリ１００をチルト(tilt)させた状態でスキャンが行われる。ガントリ１００のチルトは、内蔵のチルト機構によって行われる。
なお、テーブル２００は、図２に示すように、天板２０２がベース２０８に対して垂直に昇降する方式のものであってよい。天板２０２の昇降は内蔵の昇降機構によって行われる。このテーブル２００においては、昇降に伴う天板２０２の水平移動は生じない。
図３に、Ｘ線照射・検出装置１１０の構成を模式的に示す。Ｘ線照射・検出装置１１０は、Ｘ線管１３０の焦点１３２から放射されたＸ線１３４をＸ線検出器１５０で検出するようになっている。
Ｘ線１３４は、図示しないコリメータ(collimator)で成形されて左右対称なコーンビーム(cone beam)またはファンビーム(fan beam)のＸ線となっている。Ｘ線検出器１５０は、Ｘ線の広がりに対応して２次元的に広がるＸ線入射面１５２を有する。Ｘ線入射面１５２は円筒の一部を構成するように湾曲している。円筒の中心軸は焦点１３２を通る。
Ｘ線照射・検出装置１１０は、撮影中心すなわちアイソセンタ(isocenter)Ｏを通る中心軸の周りを回転する。中心軸は、Ｘ線検出器１５０が形成する部分円筒の中心軸に平行である。
回転の中心軸の方向をｚ方向とし、アイソセンタＯと焦点１３２を結ぶ方向をｙ方向とし、ｚ方向およびｙ方向に垂直な方向をｘ方向とする。これらｘ，ｙ，ｚ軸はｚ軸を中心軸とする回転座標系の３軸となる。
図４に、Ｘ線検出器１５０のＸ線入射面１５２の平面図を模式的に示す。Ｘ線入射面１５２は検出セル(cell)１５４がｘ方向とｚ方向に２次元的に配置されたものとなっている。すなわち、Ｘ線入射面１５２は検出セル１５４の２次元アレイ(array)となっている。なお、ファンビームＸ線を用いる場合は、Ｘ線入射面１５２は検出セル１５４の１次元アレイとしてよい。
個々の検出セル１５４はＸ線検出器１５０の検出チャンネル(channel)を構成する。これによって、Ｘ線検出器１５０は多チャンネルＸ線検出器となる。検出セル１５４は、例えばシンチレータ(scintillator)とフォトダイオード(photo diode)の組合せによって構成される。
図５に、オペレータコンソール３００による画像再構成のプロセス(process)の概要を示す。本プロセスによって、画像再構成方法に関する本発明を実施するための最良の形態の一例が示される。各プロセスは、内蔵のコンピュータによって実行される。
図５に示すように、プロセス１０１でイテラティブ・リコンストラクションを行い、プロセス１０３で補間を行って画像を出力する。プロセス１０１でのイテラティブ・リコンストラクションは、例えばオーダード・サブセッツ・エクスペクテーション・マキシマイゼーション(OSEM)によって行われる。これにより、画像再構成を適切に行うことができる。
なお、イテラティブ・リコンストラクションは、OSEMに限らず、マルティプリカティブ・アルジェブレイック・リコンストラクション(multiplicative algebraic reconstruction: MART)やパラボロイダル・サロゲート・コーディネイト・ディセント(paraboloidal surrogate coordinate descent: PSCD)等によって行ってもよい。以下、OSEMによる画像再構成について説明するが、MARTやPSCDによる画像再構成も同様になる。
イテラティブ・リコンストラクションには、投影データとプロセス１０５で計算されたシステムマトリクスが用いられる。システムマトリクスの計算には、予め記憶されている部分マトリクスが用いられる。
部分マトリクスについて説明する。その前に、部分マトリクスの前提となるシステムモデルについて説明する。図６にシステムモデルの模式的構成を示す。図６に示すように、システムモデルは複数の同心円環で構成される。
各円環は等角度のセクタ(sector)に分割されている。分割の角度はビューの角度ステップ(step)に等しいかまたはその整数倍である。なお、分割数は全円環を通じて同一である必要はなく、外側のものほど多くしてもよい。このようにすることにより、セクタの長さをシステムモデル全体にわたってほぼ均一にすることができる。また、円環の幅も同一である必要はなく個々に異ならせてもよい。
このようなシステムモデルに関する１つのビュー（ビュー番号をｉとする）の投影は、ソース(source)から放射状に伸びる複数の投影線によって得られる。これら投影線は、例えばファンビーム(fan beam)Ｘ線を形成する個別Ｘ線に相当する。ファンビームＸ線の中心ビームに相当する投影線は、システムモデルの中心を通り、かつ、円環を分割する複数の半径のうちの１つに垂直である。このような投影線を用いることにより、ファンビームＸ線による投影データに適応することができる。
ファンビームＸ線の代わりにパラレルビーム(parallel beam)Ｘ線が用いられるときは、複数の投影線はすべて平行となり、かつ、円環を分割する半径のうちの１つに垂直となる。このような投影線を用いることにより、パラレルビームＸ線による投影データに適応することができる。以下、放射状の投影線の例で説明するが、平行投影線の場合も同様である。
各投影線は、システムモデルにおける特定のセクタをそれぞれ通過する。ここで、ｊ番目の投影線が通過する任意の１つのセクタ（強調して表示）に着目すると、部分マトリクスは、１つのビューの全投影線を通じてのこのようなセクタの集まりとして定義される。これはまた、ビューｉについてのシステムマトリクスともなる。このような部分マトリクスの、全ビューを通じての集まりがシステムマトリクスである。以下、部分マトリクスを単にマトリクスともいう。
システムマトリクスのエレメント(element)は、１つのセクタによって切り取られる投影線ｊの長さａijで与えられる。このようにすることにより、システムマトリクスを簡易化することができる。
エレメントは、投影ストライプ(strip)ｋが１つのセクタ（強調して表示）によって切り取られる部分的エリアａikで与えてもよい。このようにすることにより、システムマトリクスを精密化することができる。
図７に部分的エリアの詳細を示す。網掛けした部分が、セクタによって切り取られた投影ストライプの部分的エリアである。このエリアの大きさは、この部分の面積Ｓikを
のように正規化して表すのがよい。ここで、ｌtopおよびｌbottomは、それぞれ、セクタに関するストライプｋのソース側およびその反対側の弧の長さである。このようにすることにより、エレメントを規格化することができる。
図８によって部分マトリクスの特質を説明する。図８の（ａ），（ｂ）は、引き続く２つのビュー１，２において、投影線ｊがシステムモデルを通過する状態をそれぞれ示す。ただし、ここでは、簡単のためにビューの全数を８とする。
図８の（ａ）、（ｂ）に示すように、投影線ｊは、ビュー１ではセクタ１を通過し、ビュー２ではセクタ２を通過する。システムモデルにおけるセクタの分割角度がビューの角度ステップに等しいので、ビュー１，２における投影線ｊ上のセクタ１，２の位置は全く同じになる。
したがって、ビュー２において投影線ｊがセクタ２を通過する状態は、ビュー１において投影線ｊがセクタ１を通過する状態と全く同じになる。このため、ビュー２についての部分マトリクスのエレメントａ'ijは、ビュー１についての部分マトリクスのエレメントａijと同一の値を持つ。
同様な関係がビュー１，２の間で全ての投影線について成立し、かつ、他の全てのビュー３，４，５，６，７，８に対しても成立する。これによって、部分マトリクスの構造は、ビュー番号以外は全てのビューについて同一となるので、１つのビューの部分マトリクスを予め計算して記憶しておけは、それを変換することによって残りの全ての部分マトリクスを作成することができる。セクタによって切り取られた投影ストライプの部分的エリアをエレメントとする場合も同様である。
図９に、画像再構成のさらに詳細なプロセスを示す。図９に示すように、プロセス５０１で、初期化並びに投影データとマトリクスの読み込みを行う。マトリクスは予め計算されてメモリに記憶されている。このマトリクスは部分マトリクスすなわち１ビュー分のシステムマトリクスである。このため、その計算はコンピュータにとって大きな負担にはならず、かつ、それを記憶するにも大してメモリを消費しない。
プロセス５０３で投影データのサブセットを選択し、プロセス５０５でサブセット用のマトリクスを求める。サブセット用のマトリクスは、メモリから読み出した部分マトリクスからの計算によって求められる。この計算は単純な変換であるからコンピュータの負担は軽少である。
プロセス５０７で、OSEMの１回のイテレーションを行い、プロセス５０９でイテレーション継続の要否を判定し、必要と判定したときはプロセス５０３に戻る。そして、プロセス５０３における投影データのサブセットの選択、プロセス５０５におけるマトリクス計算、プロセス５０７におけるイテレーションを順次行う。このようにして、イテラティブ・リコンストラクションによる画像再構成が行われる。サブセットのマトリクスの計算をイテレーションごとに行うので、システムマトリクスを適時に作成することができる。
プロセス５０９でイテレーション不要と判定したときは、プロセス５１１で補間演算を行う。補間演算は、再構成画像が図１０の左側に示すように、セクタ形のピクセル(pixel)を持つ画像となるので、それを右側のように矩形のピクセルを持つ画像に変換するために行われる。補間演算には、例えばバイリニア・インターポレーション(bilinear interpolation)が利用される。これによって、画像変換を適切に行うことができる。なお、補間演算は、バイリニア・インターポレーションに限らず、適宜の補間アルゴリズム(algorithm)を用いてよい。
図１１に、上記のような画像再構成を行うときの所要時間を、従来方法による画像再構成の所要時間と対比して示す。図１１に示すように、本方法を使用することにより、画像再構成の速度が従来の約４倍ないし２０倍に向上する。
図１２に、ファントム(phantom)の再構成画像を、従来方法と本方法について対比して示す。（ａ）はファントムの実像、（ｂ）は従来方法による再構成画像、（ｃ）は本方法による再構成画像である。本方法により従来方法と変わらない再構成画像が得られていることが明らかである。
図１３に、ファントムの実像、従来方法による再構成画像および本方法による再構成画像の１次元プロファイル(profile)を示す。破線がファントム実像のプロファイルであり、二点鎖線が従来方法による再構成画像のプロファイルであり、実線が本方法による再構成画像のプロファイルである。従来方法による再構成画像および本方法による再構成画像のプロファイルに差がないので二点鎖線と実線は完全に重なっている。
図１４に、３種類のファントムの本方法による再構成画像を、フィルタード・バックプロジェクション(filtered back projection: FBP)による再構成画像と対比して示す。（ａ１）（ａ２）（ａ３）がFBPによる再構成画像、（ｂ１）（ｂ２）（ｂ３）が本方法による再構成画像、（ｃ１）、（ｃ２）、（ｃ３）が差分画像である。本方法によりFBPとほとんど差がない画像が得られていることが明らかである。
本発明を実施するための最良の形態の一例のＸ線ＣＴ装置の構成を示す図である。 本発明を実施するための最良の形態の一例のＸ線ＣＴ装置の構成を示す図である。 Ｘ線照射・検出装置の構成を示す図である。 Ｘ線検出器のＸ線入射面の構成を示す図である。 本発明を実施するための最良の形態の一例のプロセスの概要を示す図である。 システムモデルとその投影線を示す図である。 投影ストライプの一部を示す図である。 引き続く２つのビューにおけるシステムモデルと投影線を示す図である。 本発明を実施するための最良の形態の一例のプロセスを示す図である。 補間による画像の変換を示す図である。 本方法と従来方法の再構成時間の対比を示す図である。 本方法と従来方法の再構成画像の対比を中間調の写真によって示す図であ る。 本方法と従来方法の再構成画像プロファイルの対比を示す図である。 本方法とFBPの再構成画像の対比を中間調の写真によって示す図である。

従来方法のプロセスを示す図である。
符号の説明
１０ ： 被検体
１００ ： ガントリ
１１０ ： Ｘ線照射・検出装置
１３０ ： Ｘ線管
１３２ ： 焦点
１３４ ： Ｘ線
１５０ ： Ｘ線検出器
１５２ ： Ｘ線入射面
１５４ ： 検出セル
２００ ： テーブル
２０２ ： 天板
２０４ ： クレードル
２０６ ： 支柱
２０８ ： ベース
３００ ： オペレータコンソール
３０２ ： ディスプレイ
１０１−１０３ プロセス
５０１−５１１ プロセス

Claims (18)

  1. それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成し、
    前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成し、
    前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成し、
    前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行い、
    得られた画像を矩形ピクセルの画像に変換する、
    ことを特徴とする画像再構成方法。
  2. 前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影線の長さである、
    ことを特徴とする請求項１に記載の画像再構成方法。
  3. 前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影ストライプの部分的エリアである、
    ことを特徴とする請求項１に記載の画像再構成方法。
  4. 前記部分的エリアは正規化されたエリアである、
    ことを特徴とする請求項３に記載の画像再構成方法。
  5. 前記他のビューのシステムマトリクスの作成をイテレーションごとに行う、
    ことを特徴とする請求項１ないし請求項４のうちのいずれか１つに記載の画像再構成方法。
  6. 前記イテラティブ・リコンストラクションはオーダード・サブセッツ・エクスペクテーション・マキシマイゼーション(OSEM)である、
    ことを特徴とする請求項１ないし請求項５のうちのいずれか１つに記載の画像再構成方法。
  7. 前記矩形ピクセルの画像への変換をバイリニア・インターポレーションによって行う、
    ことを特徴とする請求項１ないし請求項６のうちのいずれか１つに記載の画像再構成方法。
  8. 前記投影線はファンビームＸ線に対応する投影線である、
    ことを特徴とする請求項１ないし請求項７のうちのいずれか１つに記載の画像再構成方法。
  9. 前記投影線はパラレルビームＸ線に対応する投影線である、
    ことを特徴とする請求項１ないし請求項７のうちのいずれか１つに記載の画像再構成方法。
  10. 被検体をＸ線でスキャンして投影データを収集するデータ収集装置と、収集された投影データに基づいて画像を再構成する画像再構成装置とを有するＸ線ＣＴ装置であって、
    前記画像再構成装置は、
    それぞれ複数のセクタに分割された複数の同心円環によってシステムモデルを構成するシステムモデル構成手段と、
    前記システムモデルの個々のセクタのうち１つのビューの各投影線上に位置するセクタによってそのビューのシステムマトリクスを構成するシステムマトリクス構成手段と、
    前記１つのビューのシステムマトリクスを変換して他のビューのシステムマトリクスを作成するシステムマトリクス作成手段と、
    前記システムマトリクスと投影データを用いてイテラティブ・リコンストラクションにより画像再構成を行う画像再構成手段と、
    得られた画像を矩形ピクセルの画像に変換する変換手段と、
    を具備することを特徴とするＸ線ＣＴ装置。
  11. 前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影線の長さである、
    ことを特徴とする請求項１０に記載のＸ線ＣＴ装置。
  12. 前記システムマトリクスのエレメントは、前記セクタによって切り取られる投影ストライプの部分的エリアである、
    ことを特徴とする請求項１０に記載のＸ線ＣＴ装置。
  13. 前記部分的エリアは正規化されたエリアである、
    ことを特徴とする請求項１２に記載のＸ線ＣＴ装置。
  14. 前記システムマトリクス作成手段は、前記他のビューのシステムマトリクスの作成をイテレーションごとに行う、
    ことを特徴とする請求項１０ないし請求項１３のうちのいずれか１つに記載のＸ線ＣＴ装置。
  15. 前記イテラティブ・リコンストラクションはオーダード・サブセッツ・エクスペクテーション・マキシマイゼーション(OSEM)である、
    ことを特徴とする請求項１０ないし請求項１４のうちのいずれか１つに記載のＸ線ＣＴ装置。
  16. 前記変換手段は、前記矩形ピクセルの画像への変換をバイリニア・インターポレーションによって行う、
    ことを特徴とする請求項１０ないし請求項１５のうちのいずれか１つに記載のＸ線ＣＴ装置。
  17. 前記投影線はファンビームＸ線に対応する投影線である、
    ことを特徴とする請求項１０ないし請求項１６のうちのいずれか１つに記載のＸ線ＣＴ装置。
  18. 前記投影線はパラレルビームＸ線に対応する投影線である、
    ことを特徴とする請求項１０ないし請求項１６のうちのいずれか１つに記載のＸ線ＣＴ装置。
JP2007054179A 2006-03-10 2007-03-05 画像再構成方法およびｘ線ｃｔ装置 Active JP5162143B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CNB2006100739625A CN100565586C (zh) 2006-03-10 2006-03-10 图像重建方法和x射线ct设备
CN200610073962.5 2006-03-10

Publications (2)

Publication Number Publication Date
JP2007236942A true JP2007236942A (ja) 2007-09-20
JP5162143B2 JP5162143B2 (ja) 2013-03-13

Family

ID=38478939

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007054179A Active JP5162143B2 (ja) 2006-03-10 2007-03-05 画像再構成方法およびｘ線ｃｔ装置

Country Status (3)

Country Link
US (1) US7656991B2 (ja)
JP (1) JP5162143B2 (ja)
CN (1) CN100565586C (ja)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8503968B2 (en) * 2007-01-19 2013-08-06 Samsung Electronics Co., Ltd. Method and system for power saving in wireless communications
US8699421B2 (en) * 2007-01-19 2014-04-15 Samsung Electronics Co., Ltd. Method and system for wireless communication using channel selection and bandwidth reservation
WO2011148277A2 (en) * 2010-05-27 2011-12-01 Koninklijke Philips Electronics N.V. Improved reconstruction for cone-beam computed tomography imaging with off-center flat panel detector
CN102376097A (zh) * 2010-08-25 2012-03-14 东软飞利浦医疗设备系统有限责任公司 非对称检测器ct迭代重建方法
US8620054B2 (en) * 2010-09-08 2013-12-31 SOCPRA—Sciences Sante et Humaines, s.e.c. Image reconstruction based on accelerated method using polar symmetries
CN101980302A (zh) * 2010-10-22 2011-02-23 南方医科大学 投影数据恢复导引的非局部平均低剂量ct重建方法
US8503753B2 (en) * 2010-12-02 2013-08-06 Kabushiki Kaisha Toshiba System and method for triangular interpolation in image reconstruction for PET
CN102682199B (zh) * 2012-04-28 2015-08-26 昆明理工大学 基于二维索引计算投影系数构造系统矩阵的简易方法
CN104183286B (zh) * 2013-05-23 2017-11-28 环境保护部核与辐射安全中心 用于堆芯熔融物状态监测的图像重建方法、装置和系统
US9406154B2 (en) 2013-06-14 2016-08-02 General Electric Company Iterative reconstruction in image formation
BE1024052B1 (nl) * 2013-12-03 2017-11-08 Layerwise N.V. Werkwijze en inrichting voor het kalibreren van meerdere energiestralen voor het additief vervaardigen van een object
CN104361615B (zh) * 2014-10-22 2017-07-11 南方医科大学 一种采用圆轨道扇形束x射线ct扫描机快速重建断层图像的方法
CN108428253A (zh) * 2018-03-12 2018-08-21 武汉大学 一种构造虚拟doi及相应系统矩阵提高pet图像重建质量的方法
CN110264537B (zh) * 2019-06-13 2023-07-18 上海联影医疗科技股份有限公司 Pet图像重建方法、系统、可读存储介质和设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060013354A1 (en) * 2004-07-16 2006-01-19 Bjoern Heismann Method for reconstructing tomograms from detector measured data of a tomography unit

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5483567A (en) * 1994-09-22 1996-01-09 Wisconsin Alumni Research Foundation Computerized method for converting polar tomographic data to a cartesian image format
US5909476A (en) * 1997-09-22 1999-06-01 University Of Iowa Research Foundation Iterative process for reconstructing cone-beam tomographic images
DE10211485A1 (de) * 2002-03-15 2003-09-25 Philips Intellectual Property Verfahren zur Bestimmung einer Objektfunktion
DE102004034500A1 (de) * 2004-07-16 2006-02-09 Siemens Ag Verfahren zur Rekonstruktion von Schnittbildern aus Detektormessdaten eines Tomographiegerätes
US7127095B2 (en) 2004-10-15 2006-10-24 The Brigham And Women's Hospital, Inc. Factor analysis in medical imaging
US7173248B2 (en) 2004-10-20 2007-02-06 General Electric Company Methods and systems for positron emission tomography data correction
US7272205B2 (en) * 2004-11-17 2007-09-18 Purdue Research Foundation Methods, apparatus, and software to facilitate computing the elements of a forward projection matrix
US7251306B2 (en) * 2004-11-17 2007-07-31 General Electric Company Methods, apparatus, and software to facilitate iterative reconstruction of images
US7381959B2 (en) 2005-08-17 2008-06-03 General Electric Company Technique for reconstructing PET scan images
US7606348B2 (en) * 2006-02-09 2009-10-20 L-3 Communications Security and Detection Systems Inc. Tomographic imaging systems and methods
US8314796B2 (en) * 2006-02-24 2012-11-20 The Board Of Trustees Of The Leland Stanford Junior University Method of reconstructing a tomographic image using a graphics processing unit
CN101071111B (zh) * 2006-05-08 2011-05-11 清华大学 一种多视角航空集装箱安全检查系统及方法
US7949172B2 (en) * 2007-04-27 2011-05-24 Siemens Medical Solutions Usa, Inc. Iterative image processing

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060013354A1 (en) * 2004-07-16 2006-01-19 Bjoern Heismann Method for reconstructing tomograms from detector measured data of a tomography unit

Also Published As

Publication number Publication date
US20070211846A1 (en) 2007-09-13
CN101034479A (zh) 2007-09-12
US7656991B2 (en) 2010-02-02
CN100565586C (zh) 2009-12-02
JP5162143B2 (ja) 2013-03-13

Similar Documents

Publication Publication Date Title
JP5162143B2 (ja) 画像再構成方法およびｘ線ｃｔ装置
US20200196965A1 (en) Systems and methods for controlling motion of detectors having moving detector heads
USRE43272E1 (en) X-ray computerized tomographic apparatus
EP2053972B1 (en) Computed tomography image acquisition
JP6246936B2 (ja) Ｘ線撮像装置および画像再構成方法
US9613440B2 (en) Digital breast Tomosynthesis reconstruction using adaptive voxel grid
JP4585167B2 (ja) Ｘ線コンピュータ断層撮影システム
RU2534612C2 (ru) Система визуализации с массивом из множества детекторов
JP6571313B2 (ja) 医用画像診断装置及び制御方法
WO2014192775A1 (ja) Ｘ線ｃｔ装置及び画像診断装置
JP2013022455A (ja) 放射線画像診断装置及び撮像方法
CN114067013A (zh) 用于经由单应重采样变换进行重投影和反投影的系统和方法
US9895113B2 (en) Systems and methods for planar imaging with detectors having moving detector heads
JP2004160218A (ja) Ｘ線コンピュータ断層装置、ｘ線コンピュータ断層装置制御方法、及びｘ線コンピュータ断層撮影プログラム
JP2007151806A (ja) Ｘ線ｃｔ撮像方法およびｘ線ｃｔ装置
US5457321A (en) Spect apparatus utilizing fan-parallel data conversion
CN115005847A (zh) 基于可寻址纳米线冷阴极x射线平面源的静态ct系统
JP4433487B2 (ja) Ｘ線ｃｔ装置
JP2006110332A (ja) Ｘ線コンピュータ断層撮影装置
JP6129914B2 (ja) Ｘ線コンピュータ断層撮影装置
US20180259659A1 (en) Systems and methods for image quality enhancement for out of focus regions for multi-head camera
JP2009118950A (ja) Ｘ線ｃｔ装置
US20240169531A1 (en) Medical image processing apparatus, medical image processing method, and model generation method
JP2003225230A (ja) コンピュータ断層撮影装置
JP7466401B2 (ja) 医用画像診断装置

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20100118

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20111114

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120904

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121026

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121217

R150 Certificate of patent or registration of utility model

Ref document number: 5162143

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20151221

Year of fee payment: 3

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