JP2010510833A - 円状及び直線状軌道を用いた医用撮像における断層撮影法の再構成のための方法及びシステム - Google Patents

円状及び直線状軌道を用いた医用撮像における断層撮影法の再構成のための方法及びシステム Download PDF

Info

Publication number
JP2010510833A
JP2010510833A JP2009538528A JP2009538528A JP2010510833A JP 2010510833 A JP2010510833 A JP 2010510833A JP 2009538528 A JP2009538528 A JP 2009538528A JP 2009538528 A JP2009538528 A JP 2009538528A JP 2010510833 A JP2010510833 A JP 2010510833A
Authority
JP
Japan
Prior art keywords
data
circular
projection data
linear
detector
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
JP2009538528A
Other languages
English (en)
Other versions
JP5161888B2 (ja
Inventor
ザミャチン，アレキサンダー・エー
チャン，ビー−シャン・エス
シルバー，マイケル・ディー
知 中西
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.)
Toshiba Corp
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of JP2010510833A publication Critical patent/JP2010510833A/ja
Application granted granted Critical
Publication of JP5161888B2 publication Critical patent/JP5161888B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

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)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Optics & Photonics (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)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

対象物のボリュームデータ画像を再構成する方法は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、受信された円状投影データから前処理された円状投影データを生成し、対象物に関する直線状経路に沿って検出器によって収集された直線状投影データを受信し、受信された直線状投影データから前処理された直線状投影データを生成し、前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて対象物に関する再構成された円状経路ボリューム画像を生成し、前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて対象物に関する再構成された直線状経路ボリューム画像を生成し、再構成された円状経路ボリューム画像と再構成された直線状経路ボリューム画像とを結像して、対象物のボリューム画像を生成する。装置とコンピュータプログラム製品もまた記述される。
【選択図】図１

Description

本発明は、概して医用画像の再構成に関する。より詳細には、本発明は、医用画像の質と円状及び直線状軌道スキャンからの投影データに基づく３次元コーンビームコンピュータ断層撮影（ＣＴ）再構成の効率との向上のための新しい方法及びシステム装置に関する。
本発明は、下記の参考文献のリスト内で参照され記載されている様々な技術の使用を含む。これら参考文献は、括弧内の対応する参照番号により明細書中で引用され、参照によりここに本質的に組込まれる。
Figure 2010510833
コンピュータ断層撮影（ＣＴ）装置により生成された再構成画像の質と効率とは、ＣＴ装置の全般的な有効性にとって重要である。画像を再構成する際に用いられるアルゴリズムは、質と効率とに影響を与える。
それに加え、円状収集（circular acquisition）における１回転で全ての人体器官を覆うＣＴスキャナに対する臨床的な要求がある。円状軌道は、心臓／冠状動脈撮像やパーフュージョン等のような動的検査においてヘリカルよりも優位である。例えば、人の心臓を覆うことは、アキシャル方向に関して１０〜１２センチメートルの検出器サイズを要求することを明記しておく。
たった１回のスキャンでほとんどの人体器官をスキャンするために、大ボリューム範囲を備えるスキャナが要求される。このような大ボリューム範囲は、医用ＣＴにおいて広く用いられているフェルドカンプ（Feldkamp）、デイビス（Davis）、クレス（Kress）のアルゴリズム（ＦＤＫ）に基づく慣習的な近似又は非厳密な再構成法に対する挑戦である。
ＦＤＫは、高い達成可能空間解像度（high achievable spatial resolution）を有する円状軌道Ｘ線ＣＴデータを再構成するための効率的な方法である。一方、他のアルゴリズムは、リビニング（ribinning）のような空間解像度を低下させる幾つかのステップを含む。しかしながら、ＦＤＫのコーンビームジオメトリは、円状軌道の理論的な不完全性、特に体軸方向に関する急な減弱変化がおきる場所、に起因するコーンビームアーチファクトを受ける。
一方、厳密なフィルタードバックプロジェクション（ＦＢＰ）再構成アルゴリズムが提案されている［１１］。Ｘ線軌跡が単一の円を構成している場合、コーンビームの厳密な再構成は不可能であることが知られている。なぜならば、軌跡に交差しない多くの平面が存在し、これら平面の間のデータが計測できないからである。換言すれば、例えば図１に示すように、円状データ１２は、空間１０を覆う不完全なラドン空間を与える。この制限を解消するため、図２に示すように、円状スキャン軌跡２２を直線状スキャン軌跡２０で補うことが提案され、様々な円及び直線状アルゴリズムが提案されている［６〜１０］。しかしながら、［６］〜［１０］で提案されたアルゴリズムは、商業的なＣＴスキャナには全く適していなかった。すなわち、２つの要因：（ａ）それらアルゴリズムの非能率な構成に起因する画像再構成のために長時間を要求する（有効なshift-invariantなフィルタードバックプロジェクション型ではなかった）、（ｂ）製造に実用的でない、とても大きい検出器を要求する、すなわち、同等には、いわゆる“長体問題（long object problem）”（これは、人体スキャンにおいてよくあることであり、検出器の外側にまで亘る対象物（体軸方向に関して打ち切られる）は、再構成可能できないことを意味する）を受ける。
カツェビッチ（Kathevich）の厳密な直線状＋円状再構成法は、円状及び直線状データの両方に対するshift-invariant（ピクセル依存なしに独立して行なわれた、再構成された対象物のコンボルーション）と、フィルタードバックプロジェクション（ＦＢＰ）とを含む。
しかしながら、カツェビッチによるアルゴリズムの厳密な方法は、現状のＣＴスキャナのハードウェアに適用するのが容易でない。なぜならば、ヒルベルトコンボルーションは周波数領域で適用され、ヒルベルトコンボルーションは検出されたデータに対する微分を要求するからである。これは、空間解像度の低下を起こしうる。それに加え、ヒルベルト変換は、リビニングを要求する。これは、計算上の効率を相当に低下させ、解像度を多少低下させる。このように、円状及び直線状データの両方に対するそのようなアルゴリズムは、多くの時間を消費しうり、付加的な高価な処理資源を要求し、そして商業的なＣＴスキャナに適さない。
本発明に係るある目的は、再構成画像の質を向上することとコーンビームアーチファクトを除去することとである。次節において我々は、コーンビームアーチファクトを除去するための、円状スキャンと直線状スキャンとを組合せる方法を記述する。
従って、上述の技術における再構成アルゴリズムの問題を解消するため、本発明は、方法、システム装置、及びコンピュータプログラム製品を提供する。
本発明のある目的は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成するデータ収集部と、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成する円状投影データ再構成部と、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成する直線状投影データ再構成部と、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する画像ボリューム処理部と、を具備する新規のボリューム画像再構成装置を提供することにある。
本発明の他の目的は、前記直線状経路が前記円状経路の平面に対して垂直でない、新規の装置を提供することにある。
本発明の他の目的は、前記円状投影データ再構成部が、バックプロジェクション部と方向性フィルタリング部とを具備し、前記バックプロジェクション部が、フィルタリングされた検出器空間データに基づいて前記再構成された円状経路ボリューム画像を生成し、前記方向性フィルタリング部が、前記前処理された円状投影データを、予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、リビニングされたデータを生成するリサンプリング部と、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記フィルタリングされた検出器空間データを生成する逆リビニング部と、を備える、新規の装置を提供することにある。
本発明の他の目的は、前記円状投影データ再構成部が、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
本発明の他の目的は、前記円状投影データ再構成部が、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するハイブリッドフィルタリング部と、前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
本発明の他の目的は、前記データ収集部が、さらに、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された円状投影データを受信し、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、新規の装置を提供することにある。
本発明の他の目的は、前記直線状投影データ再構成部が、前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
本発明の他の目的は、前記直線状投影データ再構成部が、前記前処理された直線状投影データを、２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、ダウンサンプリングされたデータを生成するダウンサンプリング部と、前記ダウンサンプリングされたデータの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、前記フィルタリングされた検出器空間データをバックプロジェクションし、ダウンサンプリングされた直線状経路ボリューム画像を生成するバックプロジェクション部と、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングし、前記再構成された直線状経路ボリューム画像を生成するアップサンプリング部と、を備える新規の装置を提供することにある。
本発明の他の目的は、前記対象物に電流に対応する放射線を照射するものであって、前記電流は２０ｍＡよりも低い又は等しいＸ線管と、前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する検出器と、をさらに備える新規の装置を提供することにある。
本発明の他の目的は、前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整する寝台不正確さ補償部をさらに備え、前記直線状投影データ再構成部が、さらに前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、新規の装置を提供することにある。
本発明の他の目的は、前記直線状投影データ再構成部が、前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える新規の装置を提供することにある。
本発明の他の目的は、前記対象物を放射線で制御可能に照射する線源と、前記対象物が前記放射線で照射された場合、前記直線状投影データを生成する検出器と、前記線源を制御して前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射する動的コリメーション部と、をさらに備える新規の装置を提供することにある。
本発明の他の目的は、前記データ収集部が、さらに、異なる捕獲時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる捕獲時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記画像ボリューム処理部が、さらに、前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成するものであり、前記対象物の前記複数のボリューム画像を前記画像ボリューム処理部から統合し、前記複数のボリューム画像の統合されたデータを生成する統合部と、前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する合成部と、をさらに備える新規の装置を提供することにある。
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、前記合成部が、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、新規の装置を提供することにある。
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、前記データ収集部が、さらに、異なる獲得時間に円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、新規の装置を提供することにある。
本発明の他の目的は、前記データ収集部が、さらに、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、前記画像ボリューム処理部は、さらに、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記対象物の前記部分の前記複数のボリューム画像を合成し、前記対象物の前記ボリューム画像を生成する、新規の装置を提供することにある。
本発明の他の目的は、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成し、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する、対象物のボリューム画像を再構成する新規の方法を提供することにある。
本発明の他の目的は、前記直線状経路が、前記円状経路の平面に対して垂直でない、新規の方法を提供することにある。
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらに、リビニングされたデータを生成するために、前記前処理された円状投影データを予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、前記フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションする、新規の方法を提供することにある。
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらにプレ重みづけされたデータを得るために、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションする、新規の方法を提供することにある。
本発明の他の目的は、前記再構成された円状経路ボリューム画像を生成することは、さらに、プレ重みづけされたデータを得るために、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってハイブリッドフィルタリングし、前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションする、新規の方法を提供することにある。
本発明の他の目的は、前記円状投影データを受信することは、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状データを生成することは、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、新規の方法を提供することにある。
本発明の他の目的は、前記再構成された直線状経路ボリューム画像を生成することは、さらに、微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションする、新規の方法を提供することにある。
本発明の他の目的は、さらに、ダウンサンプリングされたデータを生成するために、前記前処理された直線状投影データを２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、微分係数データを生成するために、前記ダウンサンプリングされたデータの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、ダウンサンプリングされた直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するために、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングする、新規の方法を提供することにある。
本発明の他の目的は、さらに、前記対象物に、Ｘ線管内の電流であって２０ｍＡよりも低い又は等しい前記電流に対応する放射線を照射し、前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する、新規の方法を提供することにある。
本発明の他の目的は、さらに、前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整することによって寝台不正確さ補償し、前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、新規の方法を提供することにある。
本発明の他の目的は、前記直線状経路ボリューム画像を生成することは、さらに、微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションする、新規の方法を提供することにある。
本発明の他の目的は、さらに、前記対象物を放射線で制御可能に照射し、前記対象物が前記放射線で照射された場合、前記直線状投影データを生成し、前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射するように前記照射することを制御する、新規の方法を提供することにある。
本発明の他の目的は、前記円状投影データを前記受信することは、異なる獲得時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記結合することは、前記前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成することであり、さらに、前記対象物の前記複数のボリューム画像を統合し、前記複数のボリューム画像の統合されたデータを生成し、前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する、新規の方法を提供することにある。
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、さらに、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、新規の方法を提供することにある。
本発明の他の目的は、前記対象物は繰り返し位相を示すものであり、さらに、前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、前記円状経路の複数部分に沿って、異なる獲得時間に前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、新規の方法を提供することにある。
本発明の他の目的は、前記円状投影データを前記受信することは、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記合成することは、前記対象物の前記ボリューム画像を生成するために、前記対象物の前記部分の前記複数のボリューム画像を合成する、新規の方法を提供することにある。
本発明の他の目的は、コンピュータにより実行された場合に、対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成し、前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、を具備するステップをコンピュータに実行させるコンピュータプログラムインストラクションを記憶するコンピュータ読み取り可能媒体を提供することにある。
ハードウェア開発の影響とコストとを最小にするため、本発明は、現行の装置に対する最小の変化で実行できうる円状データのランプ再構成を含む。それに加え、ランプフィルタに基づく再構成アルゴリズム（例えば、ＦＤＫ、またはザミャチン（Zamyatin）、タグチ（Taguchi）、シルバー（Silver）（ＺＴＳ）ランプフィルタに基づく）は、慣習的なＣＴ医用撮像装置でよく用いられうるので、それぞれが撮像対象の体部分に基づく特定のアプリケーションのために用意されているランプ関数の確立されたライブラリーが存在する。このように、本発明は、このランプ関数の現存するライブラリーからのランプ関数を利用可能であり、異なるアプローチに関するフィルタ関数に特有なアプリケーションの開発（例えば、ヒルベルト関数に対する体部分に特有のフィルタリングアプローチに基づくアプローチ）に関連づけられる不確定とコストとを避ける。
直線状データは、不完全な円状軌道に起因するコーンビームアーチファクトを相殺する追加項を提供する。
本アプローチは、擬人ファントムのシミュレーションデータと実測のスキャンデータとに基づいて実現可能であると証明されている。本アプローチは、スキャン中における造影剤注入、体動、及びミスアライメントに関して安定している証明されている。直線状スキャンに起因して患者が照射される余分な線量もまた、円状スキャンのそれに比して有利に少ない。本発明は、コーン角がとても大きい場合でさえも、コーンビームアーチファクトの影響のない再構成をも可能とし、さらに検出器列の数が２５６を超える場合であっても使用可能である。
本発明の実施形態によれば、円状データのランプフィルタに基づく再構成（例えば、ＦＤＫ又はＺＴＣのような）が行なわれる。さらに、低線量の直線状スキャン（ヒルベルト変換を含む方法を用いて再構成される）が追加される。直線状データは、再構成を速めるため、及びさらなる処理改善（例えば、メモリ削減や処理容量）を得るため、２から４の係数によってダウンサンプリングされうる。このように、本発明は、商業的なＣＴスキャナに適している。
幾つかの臨床応用は、円状軌道に対して直線状部分を傾かせるガントリチルトを要求する。換言すれば、その直線状部分は、円状軌道の平面に対して垂直でない。このように、本発明の実施形態は、ガントリが傾いている場合、例えば、脳撮像期間中に適している。この場合、患者の目と患者の耳の中心とを結ぶ眼窩耳孔（ＯＭ）線に対して再構成面が平行するように、スキャン前にガントリが傾かれる。このＯＭ線を含む面は、基準解剖面として供される。また、ガントリを傾けることは、患者の目のレンズへの照射線量を最小にするのに役立つ。
本発明の他の方法、システム、及びコンピュータプログラム製品は、以下の図面及び好ましい実施形態の詳細な記述の検討のもと、一般人あるいは当業者にとって明らかになろう。
本発明のより完全な理解及びそれの多くの付随する利点は、以下の詳細な記述を、以下の添付図面に関連づけて考慮して参照することによってより容易に理解できるようになるであろう。
図１は、ラドン空間における円状データ範囲の一例の図である。 図２は、直線状スキャン軌跡と円状スキャン軌跡との間の関係の一例の等角投影図である。 図３は、本発明の実施のブロック図である。 図４は、本発明の実施に係る物理座標系のブロック図である。 本発明の実施に係る画像再構成の実行方法のフローチャートである。 本発明の実施に係る直線状ボリューム画像の再構成の実行方法のフローチャートである。 図７Ａは、円状スキャン中において収集されたデータに対する、検出器上のフィルタリングラインの一例である。 図７Ｂは、直線状スキャン中において収集されたデータに対する、検出器上のフィルタリングラインの一例である。 図８は、直線状スキャンの範囲の一例のブロック図である。 図９は、平面及び湾曲した検出器のジオメトリを考慮するためのさらなる用語を有するブロック図。 図１０は、本実施形態の一例に係るｗ対ρのグラフ上におけるＢＰＪエリアのプロットである。 図１１Ａは、傾いたガントリを含む本発明の実施形態の一例の等角投影図である。 図１１Ｂは、円状軌道に関連付けられた、傾いた座標の等角投影図である。 図１２Ａは、本発明の実施形態に係る座標の等角投影図である。 図１２Ｂは、本発明の実施形態に係る座標の他の図である。 図１３Ａは、慣習的な円状ＦＤＫ法を用いた、シミュレーションデータの再構成画像の一例である。 図１３Ｂは、直線状データ再構成に組合わされたランプに基づく円状データ再構成を有する本発明の実施形態を用いた、シミュレーションデータの再構成画像の一例である。 図１４Ａは、ゲート制御された円状スキャンデータをのみを用いて画像が再構成された場合の心像動の影響を示す再構成画像である。 図１４Ｂは、本発明の実施形態に係るゲート制御された円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された場合の心臓動の影響を示す再構成画像である。 図１４Ｃは、図１４Ａと図１４Ｂとの間の差分の画像である。 図１４Ｄは、ゲート制御されていない円状スキャンデータのみを用いて画像が再構成された場合のシンチレータ心臓動の影響を示す再構成画像である。 図１４Ｅは、ゲート制御されていない円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された場合の心臓動の影響を示す再構成画像である。 図１４Ｆは、図１４Ｄと図１４Ｅとの間の差分の画像である。 図１５は、擬人ファントムの画像である。 図１６Ａは、ＦＤＫ円状スキャンデータのみを用いる再構成画像の一例である。 図１６Ｂは、直線状再構成データの一例である。 図１６Ｃは、本発明の実施形態を用いる再構成画像の一例である。 図１７Ａは、ＦＤＫ円状データのみを用いて再構成された画像の一例である。 図１７Ｂは、本発明の実施形態に係る、管電流５０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。 図１７Ｃは、図１７Ｂの例において用いられた直線状データのみの一例である。 図１７Ｄは、本発明の実施形態に係る、管電流３０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。 図１７Ｅは、図１７Ｄの例において用いられた直線状データのみの一例である。 図１７Ｆは、本発明の実施形態に係る、管電流５０ｍＡでの円状及び直線状データ用いて再構成された画像の一例である。 図１７Ｇは、図１７Ｆの例において用いられた直線状データのみの一例である。 図１８Ａは、ＦＤＫ円状データ再構成のみを用いて再構成された画像である。 図１８Ｂは、本発明の実施形態に係る、正確な寝台位置情報で再構成された画像である。 図１８Ｃは、不正確な寝台位置情報で再構成された画像である。 図１９Ａは、ＦＤＫ円状スキャンのみで得られたＦＯＶの図である。 図１９Ｂは、仮想的な拡張検出器を用いて得られた拡張ＦＯＶの図である。 図２０は、本発明の実施形態に係る直線状スキャンコリメーションを用いる直線状スキャンの照射線量を低減する方法におけるＦＯＶの図である。 図２１は、動的コリメーションに対するコリメーション窓の一例である。 図２２は、本発明に係るジオメトリの更なる図である。 図２３Ａは、直線状スキャンデータ中のフィルタラインのグラフである。 図２３Ｂは、直線状スキャンデータ中のフィルタラインの数のグラフである。 図２４は、本発明の他の実施形態に係る直線状及び円状アルゴリズムの方法のブロック図である。 図２５Ａは、ショートスキャンコンボルーションにおいて用いられるカーブのファミリーの一例である。 図２５Ｂは、ショートスキャンコンボルーションにおいて用いられるカーブの他のファミリーの一例である。 図２５Ｃは、図２５Ａにおけるカーブのファミリーの他の表現である。 図２５Ｄは、図２５Ｂにおけるカーブのファミリーの他の表現である。 図２６Ａは、複数のファミリーからのラインにおけるフルスキャンコンボルーションの一例である。 図２６Ｂは、図２６Ａの２つのファミリーのラインの重なりの１次元図である。 図２７Ａは、本発明の実施形態に係る逆（inverse）リビニングカーブの一例である。 図２７Ｂは、逆（reverse）リビニング処理の一例である。 図２８Ａは、本発明の実施形態に係る拡張再構成関心ボリューム（ＶＯＩ）の一例である。 図２８Ｂは、ＦＤＫ―ＶＯＩを拡張させたＦＤＫ―ＶＯＩの形状の線図である。 図２９Ａは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。 図２９Ｂは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。 図２９Ｃは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。 図２９Ｄは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。 図２９Ｅは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。 図２９Ｆは、本発明の他の実施形態に係る、１又はそれ以上の直線状スキャンを有する複数の円状スキャンからの結合データの一例を示す図である。
では図面、より詳細には図３を参照せよ。種々の図面の全体に亘って、同一の参照符号は、同一又は対応する部分を示している。図３は、本発明に係るＸ線コンピュータ断層撮影装置の実施形態のブロック図である。撮像装置は、略錐形状のＸ線束のコーンビームを発生するＸ線源３を有するガントリ１と、２次元状に配列された複数の検出素子（すなわち、一列に配列された複数の素子が複数列に積まれている）を含む２次元アレイ型のＸ線検出器５とを含む。Ｘ線源３と２次元アレイ型のＸ線検出器５とは、回転リング２に設置され、スライド制の寝台６に載置されている被検体に対向して向き合っている。２次元アレイ型のＸ線検出器５は、回転リング２上に搭載されている。各検出素子は、検出されたデータの１チャンネルに対応する。Ｘ線源３からのＸ線は、Ｘ線フィルタ４を介して被検体に向けられる。被検体を透過したＸ線は、２次元アレイ型のＸ線検出器５によって電気信号として検出される。
Ｘ線制御装置８は、高電圧発生器７にトリガ信号を供給する。高電圧発生器７は、トリガ信号を受信したタイミングで高電圧をＸ線源３に印加し、Ｘ線をＸ線源３から放射させる。ガントリ／寝台制御器９は、ガントリ１の回転リング２の回転と寝台６のスライドシートのスライドとを同期制御する。システム制御器１０は、全システムのコントロールセンターを構成し、被検体から見てＸ線源３が円状スキャンと直線状スキャンとを実行するようにＸ線制御器８とガントリ／寝台制御器９とを制御する。
本実施形態に係る円状スキャン期間中、回転リング２（線源３と検出器５とを含む）は、スライド制の寝台６が所定の位置に留まったまま、所定の角速度で連続的に回転し、Ｘ線は、所定角度おきに連続的または間欠的にＸ線源３から放射される。このように、寝台６が静止している時、円状検出器５は、スキャンされている対象物に関する真の円状経路に沿って移動し、螺旋経路上を移動しない。直線状スキャン中、寝台６は、Ｘ線が連続的または間欠的にＸ線源３（静的な位置に保持されている、例えば、直接的に寝台６の上に）から放射されたまま、直線的にガントリ内を所定速度で移動される。
上記の直線状スキャンの実施形態において、撮像される患者又は対象物は、線源と検出器とが所定位置に保持されている間、ガントリ内を移動される。しかしながら、本発明は、対象物と検出器との両方が移動する実施形態と同様に、対象物が所定位置に保持され、検出器が対象物に沿って移動する実施形態をも含む。円状スキャンと直線状スキャンとは、あらゆる順番で交互に行なわれる。
上記の円状スキャンの実施形態において、検出器５の経路は、円であると示された。しかしながら、本発明は、円状経路に限定されず、対象物の周りを囲む他の湾曲された経路をも含む。
２次元配列型Ｘ線検出器５の出力信号は、投影データを生成するために、特に直線状スキャン中においては直線状投影データを生成し、円状スキャン中においては円状投影データを生成するために、チャンネル毎にデータ収集部１１により増幅され、デジタル信号に変換される。投影データは、当業者に知られている方法を用いて前処理され、前処理されたデータは、データ収集部１１から再構成処理部１２に出力される。再構成処理部１２は、前処理されたデータを用いて、スキャンされた対象物に関するボリューム画像を生成する。
再構成されたボリューム画像は、ディスプレイ装置１４に送信され、そこで３次元画像または断層像として視覚的に表示される。
この明細書の目的のため、我々は、ボリューム画像をスキャンされている身体物の描写であると定義するものとする。その画像は、幾つかの撮像技術によって発生される。撮像技術の例は、テレビジョン、ＣＣＤカメラ、又はＸ線、音波、あるいは超音波撮像装置を含みうる。画像が記録される初期媒体は、電子的固体装置、画像フィルム、又は輝尽性蛍光体のような幾つかの他の装置が可能である。記録された画像は、その後、電子的（ＣＣＤ信号の場合のように）又は機械的／光学的方法（写真フィルムをデジタル化する場合や輝尽性蛍光体からのデータをデジタル化する場合のように）の組合せによって、デジタル形式に変換され得る。
本実施形態に用いられる検出器５は、１２０ミリメートルのボリューム範囲を有する、２５６列の医用ＣＴスキャナ検出器を含みうる［１−４］。当業者は、本発明が他の慣習的なＣＴスキャナ検出器にも適用するということを認識するであろう。本実施形態に係る補足的なスキャナパラメータは、表１に示される。
Figure 2010510833
図４は、本撮像装置の物理座標システムの定義の一例のブロック図である。円状軌跡上の線源位置は、ｙＣ（β）＝（Ｒｃｏｓβ，Ｒｓｉｎβ，０）である。ここで、Ｒは円状軌跡の半径、βはビュー角である。直線状軌跡上の線源位置は、ｙＬ（ｈ）＝（Ｒｃｏｓθ，Ｒｓｉｎθ，ｈ）、ここでθは所定の直線角度、ｈはｚ軸に沿う直線状軌跡上の線源の垂直座標である。慣習的な完全アルゴリズム［１１］に対する方程式は、方程式（１）に示すような、円及び直線状データの積分にそれぞれ対応する２つの積分項の加算を含む。
Figure 2010510833
方程式（１）において、ｘは再構成点、γはファン角、νは垂直検出器座標である。Ｈｓｉｎγ Ｄｉ，ｉ＝１，２，は、カーネル１／ｓｉｎγを有する１次元方向ヒルベルトフィルタである。線及び円状スキャンに対するフィルタリング方向は、異なる。
一方、本発明によれば、方程式（１）の厳密なアルゴリズムは簡約されている。第１に、湾曲した検出器上のフィルタリング方向は、円状データに対して略水平である（それらは平面検出器上において厳密に水平）。従って、本発明によれば、円状データに対してリビニングが適用されない。第２に、γ微分係数に適用されるヒルベルトコンボルーションは、ランプコンボルーションに置き換え可能であることが提案されている［１２］。しかしながら、本発明は、さらに進んでおり、ヒルベルト再構成に比して速度と効率とを改善するためだけでなく、解像度を改善するために方程式（１）の第１項の全てをＦＤＫ再構成に置き換える（比較のため［１２］を参照見よ）。
本発明によれば、円状スキャンは、ＦＤＫ再構成アルゴリズム又はＺＴＳ再構成アルゴリズムのような、ランプフィルタを有する再構成アルゴリズムを用いて再構成される。さらに、本発明によれば、直線状スキャンデータは、再構成された円状スキャンデータを補足するための低周波数補正項の計算のために再構成される。
図５は、本発明の実施形態に係る画像再構成を実行する方法のフローチャートである。図５において、その方法は、例えば、前処理器又は前処理器を含む検出器からの前処理された円状投影データが受信されるステップＳ５００を含む。ステップＳ５０２において、受信されたデータは、円状投影ボリューム画像を生成するために再構成される。ステップＳ５０１において、前処理された直線状投影データは受信され、そしてステップＳ５０３において、前処理された直線状スキャン投影データは、直線状投影ボリューム画像を生成するため、フィルタードバックプロジェクション（ＦＢＰ）ヒルベルト再構成を用いて再構成される。ステップＳ５０４において、円状投影ボリューム画像と直線状投影ボリューム画像とは、検出された対象物の最終及び全ボリューム画像を生成するために組合される。
図６は、本発明に係る、例えば、図５におけるステップＳ５０３を実行する、直線状ボリューム画像の再構成を実行する方法のフローチャートである。図６において、ステップＳ６０１は、受信された前処理後の直線状投影データの微分を実行する。ステップＳ６０２において、微分されたデータは、フィルタリングラインに沿ってリビニングされ、ステップＳ６０３において、リビニングされたデータは、フィルタリングラインに沿ってヒルベルトフィルタでフィルタリングされる。ステップＳ６０４において、ヒルベルトフィルタリングされたデータは、検出器空間又は格子点に逆リビニングされ、ステップ６０５において、リビニングされたデータは、直線状投影ボリューム画像を生成するため、３次元バックプロジェクションされる。
図７Ａは、円状スキャン中に収集されたデータに対する検出器上のフィルタリングラインの一例を示す。図７Ｂは、直線状スキャン中に収集されたデータに対する検出器上のフィルタリングラインの一例を示す。円状データにおいて、フィルタリング方向は略水平であり、直線状データに対しては、フィルタリング方向は傾いている。
直線状データのコンボルーション。初めに、ＰＩラインの概念を紹介する。再構成ピクセルをｘとすれば、ＰＩラインは、円状軌跡上の数点と直線状軌道上の数点とを接続し、ピクセルｘに交差する直線である。ＰＩラインは、各再構成ピクセルｘに対して存在することが［１１］に示されている。平面検出器上のＰＩ窓（軌跡の直線状部分と円状部分との間をＰＩラインによって追跡された）は、放物線：νｆｌ＝−｛（ｈｕ２）／（２Ｒ２）｝−ｈ／２によって与えられる。ここで、ｕ及びνｆｌは、それぞれ水平及び垂直平面検出器座標である。放物線の頂点は、（０，−ｈ／２）に位置する。フィルタリングラインのパラメータとしてρ＝ｔａｎγを用い、そして平面検出器上のラインνｆｌ＝ｍｕ＋ｂが湾曲検出器上のカーブν＝ｂｃｏｓγ＋ｍＲｓｉｎγに投影されることを考慮すると、我々は、図７Ｂに示すように（図７Ｂの垂直軸はνに対応し、水平軸はγに対応する）、コンボルーションラインの下記の方程式を得る。
Figure 2010510833
フィルタリングファミリーは、検出器エリアの二重範囲（double coverage）を提供する。すなわち、ＰＩ窓内のあらゆるピクセルに対して、そのピクセルを通過する２つのフィルタリングラインを見つけられる。フルスキャンモードにおいては、両方のラインからの寄与が加えられる。ショートスキャンモードにおいては、１つのサブファミリーのみが用いられる。
微分は、線源パラメータ（ｄ／ｄｈ）の方向に関して実行され、２つの連続した投影間の単純投影減算として実行され、それらの直線距離で割り算される。あるいは、３点以上の有限差分式で実行できる。
さらにフィルタリングラインへのリビニングは、直線状データのヒルベルト再構成の一部として実行されうる。リビニングは、カーテシアン検出器座標から非カーテシアンフィルタリング座標へのリサンプリングとして実行されることを明記しておく。通常２次元内挿が要求されるが、我々はν方向に関する１次元内挿のみを用いてリサンプリングを適用し、そして我々はγ方向に関してサンプリングを同じに保つ。
直線状データのバックプロジェクション。直線状データのバックプロジェクションは、直接的にフィルタリングラインから実行されうる。このように、ＰＩ窓内のデータのみが用いられる（すなわち、ν＞νπ（γ）、直線状データの冗長性が考慮されるように）。他の方法は、検出器座標に逆リビニングを適用することである。この場合、直線状検出器データは、ＰＩ窓でマスクされなければならず、方程式（３）によって与えられる。
Figure 2010510833
図８は、直線状スキャンの範囲の一例のブロック図である。アイソセンタにおける検出器の幅をＷとする。そのとき、直線状スキャンに必要な長さは２Ｗ、すなわち、直線状スキャンの長さは、アイソセンタにおける検出器サイズの２倍、又は概ね物理検出器（physical detector）のサイズでなければならない。この場合、直線状データは、ＦＤＫによって正確に再構成されうる六角形領域（図８において影付けされた領域）を完全に覆う。すなわち、直線状データは、ＦＤＫフィルタ又はＺＴＳフィルタ［１５］又は［１６］のようなランプに基づくフィルタを含むランプフィルタに基づくアルゴリズムに応じたフィールドオブビューの厳密な再構成に十分足りている。
図９は、平面及び湾曲検出器のジオメトリを考慮するためのさらなる用語を有する図である。Ｗは、検出器幅の半分（例えば、Ｗ＝６４ｍｍ）、Ｈは１方向に関する直線状スキャンの長さ、そしてｈは直線の線源位置、ｈ＜Ｈである。ここで我々は、ｈ＞０の場合のみ記述し、そしてｈ＜０の場合は同様に扱いうるとする。（ν，γ）は、円柱状検出器座標である。νは垂直座標であり、γはファン角である。（ｗ，ρ）は、平面検出器座標である。ｗは、垂直座標であり、ｗ＝ν／ｃｏｓγである。ρは水平座標であり、ρ＝ｔａｎγである。図７Ａ及び図７Ｂのそれぞれの垂直軸はρに対応し、水平軸はγに対応する。
平面検出器。初めに、平面検出器に対する方程式が示される。そして、それら方程式は、円柱状検出器に適用されるために、平面検出器上のラインｗ＝ｍρ＋ｂが円柱状検出器上のカーブν＝ｂｃｏｓγ＋ｍｓｉｎγに対応することを示す関係を用いて変換される。
ＰＩ窓の投影境界は、方程式（４）により与えられる。
Figure 2010510833
上記方程式によれば、検出器上のコンボルーションラインは、放物線ｗＰＩ（ρ，ｈ）の接線である。放物線上の点（ｗ，ρ）の接線の傾きは、
Figure 2010510833
によって与えられる。
このように、点（ｘ０，ｙ０）を通り傾きｍを有するラインの方程式は、
Figure 2010510833
によって与えられる。
放物線ｗＰＩ（ρ，ｈ）上の与えられた点（ｗ０，ρ０）に対して、この点を通る接線の傾きはｍ＝−ｈρ０であり、ｗ０＝−ｈ（１＋ρ０ ２）／２である。この時、この点を通る接線の方程式は、ｗ＝−ｈ（１＋ρ０ ２）／２−ｈρ０（ρ−ρ０）である。この方程式を解くことによって、我々は平面検出器上の接線の方程式を得る。
Figure 2010510833
次に、ρ０の極限、すなわちρｍｉｎ及びρｍａｘは、２つの場合、すなわちフルスキャン及びハーフスキャンのそれぞれにおいて、フィルタリングラインが検出器上のＢＰＪエリアの全体を覆うように選択される。
フルスキャンの場合において、フィルタリングファミリーは、対称的である、すなわち、ρｍｉｎ＝−ρｍａｘである。従って、ρｍａｘを探せば十分である。パラメータρｍａｘは、ＦＤＫボリューム投影のピーク点を覆うフィルタリングラインによって決定される（臨界状態）。このように、点Ａの座標は、（ρ＝０,ｗ＝Ｗ−ｈ）である。これら座標を方程式（７）に代入すると、我々は、ρ０がρｍａｘに対応するということが分かる。
Figure 2010510833
これから、
Figure 2010510833
であり、ここで、
Figure 2010510833
である。
式（９）〜（１０）は、直線状スキャンの長さＨに制限を課す。座標ρＸは、負であってはならない。これは、条件：ｈ≦２Ｗを強要する。従って、ＦＤＫボリューム再構成のために、直線状スキャンの長さは、検出器範囲の厳密に２倍、Ｈ＝２Ｗでなければならない。
図１０は、Ｗ＝６４ｍｍ，ｈ＝１００ｍｍ，及びρＸ＝０．２８における本実施形態の一例に係る、ｗ対ρのグラフ上のＢＰＪエリアのプロットである。
ハーフスキャンエリアにおいて、フィルタリングファミリーは、ＢＰＪは１つの側のみを覆うべきであるので、対称でない。我々が右からＢＰＪを覆う場合、ρｍａｘは、方程式（１０）で与えられる。
パラメータρｍｉｎは、ＢＰＪエリアによって決定される。ρ１＝ｔａｎγｍａｘとする。このとき、
Figure 2010510833
ＰＩ窓の投影境界は、
Figure 2010510833
によって与えられる。
コンボルーションライン。方程式（７）から我々は、
Figure 2010510833
であることが分かる。（今我々は、簡単のため“ρ０”の代わりに“ρ”と書いた）、そして我々は、平面検出器上のフィルタリングカーブの方程式を得る。
Figure 2010510833
円柱状検出器上のフィルタリングファミリーは、平面検出器上のファミリーと同じ指標付きパラメータを用いることを明記しておく。従って我々は、ρに対する極限を見つけるために方程式（１０）〜（１１）を用いることができる。
図１１Ａは、傾いたガントリを含む本発明の実施形態の等角投影図である。図１１Ｂは、円状軌道に関連する傾いた座標の等角投影図である。傾いたガントリにおいて、アルゴリズムの構造は上記の傾いていない場合と同様のままであり、画像再構成は傾いた座標系で実行され、このことは円状データの再構成がガントリの傾きがないのと同様の状態であることを意味する。従って、この節において、我々は、傾いたガントリにおける直線状データからの再構成のみを記述する。
本実施形態によれば、図１１Ａに示すように、ξをガントリのチルト角とし、図１２Ａにおいて、μバー（μバーは、方程式上ではμの上に―を付したものとする）は、直線状軌跡μの単位方向ベクトルとする。ここで、図１２Ａ上の影付けされた領域は、再構成ボリュームを表す（すなわち、データ外挿が必要ない最大領域）。円と直線との交点ｙ０バー（ｙ０バー、方程式上ではｙ０の上に―を付したものとする）が座標（Ｒ，０，０）を有するように、軸ｘ１，ｘ２を回転すると便利である。
本実施形態において、ベクトルμバーは、座標
Figure 2010510833
を有し、直線状線源の座標は、
Figure 2010510833
により与えられる。
方程式（１８）は、以下に示すような、平面検出器上のフィルタリングラインのための、方程式（１９）は、湾曲検出器上のフィルタリングラインのためのものである。
Figure 2010510833
方程式（１８）及び（１９）においてρＬは、方程式ρＬ＝−ｃｏｔ（λ／２）に基づくフィルタリングパラメータである。ここでλは、フィルタリング面の接点への角度である。
直線状スキャンの範囲は、
Figure 2010510833
によって定義される。ガントリの傾きがゼロの場合において、方程式（２０）が上記のようにｈｍａｘ＝２Ｗにまとめられることが容易にわかる。
本発明の全実施形態は、便宜的に、慣習的で汎用なコンピュータ又は、当事者にとって明らかなように、本発明の教示に応じてプログラムされたマイクロプロセッサーを用いて実行される。適切なソフトウェアは、ソフトウェア業界において通常のスキルを有する者にとって明らかなように、本開示の教示に基づく通常のスキルを備えたプログラマによって容易に用意されうる。特に、コンピュータの筐体は、ＣＰＵ、メモリ（例えば、ＤＲＡＭや、ＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、ＳＲＡＭ、ＳＤＲＡＭ、フラッシュＲＡＭ）、及び他のオプションの専用論理装置（例えばＡＳＩＣＳ）又は適合性のある論理デバイス（例えばＧＡＬ、及び再プログラム可能なＦＰＧＡ）を備えるマザーボードを収容する。コンピュータは、複数の入力装置（例えばキーボードやマウス）、及びモニターを制御するためのディスプレイカードもまた含む。さらにコンピュータは、フロッピー（登録商標）ディスクドライブ；他の着脱可能な媒体装置（例えばコンパクトディスク、テープ、及び着脱可能な光磁気媒体）；及びハードディスク又は他の固定された高密度メディアドライブを含んでいても良く、適切なデバイスバス（例えばＳＣＳＩバス、拡張ＩＤＥバス、又はウルトラＤＭＡバス）を用いて接続されていてもよい。コンピュータは、さらにコンパクトディスク読取装置、コンパクトディスク読み取り／書き込みユニット、又はコンパクトディスク・ジュークボックスを含んでもよい。それらは同じデバイスバス又は他のデバイスバスに接続されてもよい。
本発明に関連したコンピュータ読取り可能な媒体の例は、コンパクトディスク、ハードディスク、フロッピー（登録商標）ディスク、テープ、光磁気ディスク、ＰＲＯＭ（例えばＥＰＲＯＭ、ＥＥＰＲＯＭ、フラッシュＥＰＲＯＭ）、ＤＲＡＭ、ＳＲＡＭ、ＳＤＲＡＭ等を含む。任意の一つ又はこれらコンピュータ読取り可能な媒体の組合せに格納されて、本発明は、コンピュータの両方のハードウェアを制御するため、及びコンピュータが人間のユーザーと対話できるようにするためのソフトウェアを含む。そのようなソフトウェアは、開発ツールのようなデバイス・ドライバや、オペレーティングシステム、ユーザーアプリケーションを含んでいてもよいが、これらに制限されない。本発明のコンピュータプログラム製品は、コンピュータによって実行された時にコンピュータに本発明の方法を行なわせるコンピュータプログラムインストラクション（例えばコンピュータコード装置）を格納するあらゆるコンピュータ読み取り可能媒体を含んでいる。本発明のコンピュータコード装置は、スクリプトや、インタープリタ、ダイナミックリンクライブラリ（ＤＬＬ）、Ｊａｖａ（登録商標）クラス、完全実行可能プログラムを含むあらゆる解釈可能又は実行可能なコードメカニズムでもよいが、これらに制限されない。さらに、本発明の処理部分は、よりよい性能、信頼性、及び／又はコストのために分配（例えば（１）多重ＣＰＵ間、（２）少なくとも１つのＣＰＵと少なくとも１つの配置可能な論理装置との間で）されてもよい。例えば、輪郭（outline）又は画像は、第1のコンピュータ上で選択され、遠隔診断のための第２のコンピュータへ送信されてもよい。
その発明は、当業者に容易に明らかなように、特定用途向け集積回路の準備によって又は慣習的な構成回路の適切なネットワークを相互に接続することによって実行されてもよい。
本発明への画像データの源は、Ｘ線装置またはＣＴ装置のような任意の適切な画像収集装置が可能である。さらに、収集されたデータは、既にデジタル形式でなければデジタル化されてもよい。あるいは、収集され処理される画像データの源は、画像収集装置によって生成されたデータを格納するメモリでもよい。そのメモリは、局部メモリであってもよいし、又は遠隔メモリであってもよく、その場合、ＰＡＣＳ（Picture Archiving Computer System）のようなデータ通信網は、本発明に係わる処理のために画像データにアクセスするために用いられてもよい。
本発明の実施形態の実行は、コーンビームアーチファクトの完全な削除を含み、シミュレートされた、ｚ方向に関して急な変化を有する背骨を含むノイズのない胴部を用いて評価される。
図１３Ａは、慣習的な円状ＦＤＫ方法を用いてシミュレートされたデータの再構成画像の一例である。
図１３Ｂは、上記のように、直線状データ再構成と組合された、ランプに基づく円状データ再構成を有する本発明の実施形態を用いてシミュレートされたデータの再構成画像の一例である。図１３Ａ及び図１３Ｂにおけるディスプレイウィンドウは、Ｌ／Ｗ＝０／２００である。
本発明の実施形態の実行は、人体の心臓動を現実的に模した動的な動態ファントムを用いて評価された。
図１４Ａ〜１４Ｆは、シミュレートされた心臓動の影響の結果の一例である。
図１４Ａは、ゲート制御された円状スキャンデータのみを用いて画像が再構成された時の心臓動の影響を示す再構成画像である。
図１４Ｂは、本発明の実施形態における、ゲート制御された円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
図１４Ｃは、図１４Ａと図１４Ｂとの差分を示す。
図１４Ｄは、ゲート制御されていない円状スキャンデータのみを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
図１４Ｅは、ゲート制御されていない円状スキャンデータとゲート制御されていない直線状データとを用いて画像が再構成された時の心臓動の評価を示す再構成画像である。
図１４Ｆは、図１４Ｄと図１４Ｅとの差分を示す。
画像例は、直線状スキャンデータのゲート制御を要求することなく円状スキャンデータをゲート制御することによって、再構成が改善され、コーンビームアーチファクトが除去されたことを示す。さらに、コーンビームアーチファクトの低減は、円状スキャンデータにおける器官動によって影響されない。
ゲート制御された円状データは、心臓が比較的動いていないとき、その円状データが同一の特定心位相中に獲得された（captured）画像のみを含むように制限する効果を有する。ゲート制御は、所望の心位相中のスキャンデータのみを生成する（例えば、その位相に対応する時間のみに線源で患者を放射する）、又は複数の心位相中にデータを獲得して所望の心周期に対応するデータをのみを抽出することによって実行されうる。
擬人ファントムの投影データは、円状及び直線状軌跡に沿って２５６スライススキャナで収集された。図１５は、擬人ファントムの画像である。幾つかのデータセットは、直線放射線量管電流（ｍＡ）の設定と体部分すなわち頭、胸、及び腹も依存する性能を調査するために収集された。この例に対するスキャンパラメータは、表１に与えられている。この例において、最大コーン角は、１２度を越える。直線状スキャンのための管電流は、円状スキャンのそれに比してかなり低いことを明記しておく。
図１６Ａは、ＦＤＫ円状スキャンデータのみを用いた再構成画像の一例を示す。図１６Ｂは、直線状再構成データの一例である。図１６Ｃは、上記のような、本発明の実施形態を用いた再構成画像の一例である。図１６Ａ〜１６Ｃの例において、表示ウィンドウは、Ｌ／Ｗ＝４０／２００である。この例は、本発明のアプローチは、目に見えるコーンビームアーチファクトを消去することをさらに例証している。
さらに図１６Ｂは、直線状スキャンデータがコーンビームアーチファクト補正のみに寄与し、他のいかなる解剖学的構造を含んでいないことを例証している。直線状スキャンデータは、ｚ方向のみに関して振動数が高いことを表し、ｘｙ平面に関してとても振動数が低いことを示している。
さらに、とても低い管電流（ｍＡ）に応じて、とても低い放射線量を用いて直線状スキャンを実行することが可能である。上記の本発明の実施形態の性能評価は、コーンビームアーチファクトが５０ｍＡより低い、好ましくは２０ｍＡより低い、又は約１０ｍＡの管電流を用いて低減されうるとういうことを示す。
図１７Ａは、ＦＤＫ円状データのみを用いて再構成された画像の一例を示す。図１７Ｂは、５０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｃは、図１７Ｂの例において用いられた直線状データのみの例である。図１７Ｄは、３０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｅは、図１７Ｄの例において用いられた直線状データのみの例である。図１７Ｆは、５０ｍＡの管電流を用いて、本発明の実施形態による円状及び直線状データを用いて再構成された画像の一例を示す。図１７Ｇは、図１７Ｆの例において用いられた直線状データのみの例である。従って、この例は、患者への直線状スキャンの余分な線量が最小であり、それでもコーンビームアーチファクトが低減されていることを示している。
さらに上記の実施形態は、寝台位置に関するデータの不正確さに敏感でありうる。例えば、図１８Ａは、ＦＤＫ円状データ再構成のみを用いて再構成された画像である。図１８Ｂは、正確な寝台位置情報を用いた、上記の本発明の実施形態により再構成された画像である。図１８Ｃは、不正確な寝台位置情報、特に、＋／−０．３ｍｍの不正確さを有する寝台位置情報を用いた、上記の本発明の実施形態により再構成された画像である。図１８Ｃの例は、不正確な寝台位置情報に起因するストリークアーチファクトを示す。
本発明のさらなる実施形態は、直線状再構成又は円状再構成の１つ又は両方を生成するために用いられる再構成パラメータを調整することによる、寝台に関する不正確を補償することを含む。その再構成パラメータは、円及び直線状軌道との共通点の付近にある円状投影データと円状投影データとの間の相関に基づいて決定される（この点において両方の軌道からの投影は、一致していなければならない）。
直線状スキャンの再構成ボリューム範囲は、仮想的な拡張検出器（virtual expanded detector）を用いることによってさらに低減されうる。ＦＤＫにおける円状経路のみの再構成において、ＦＯＶ、図１９Ａの例に示すように、６角形形状により制限される（例えば、１９Ａにおける６角形形状に形作られたＲＯＩ領域）。ＦＤＫ―ＦＯＶを増加させるために、仮想的な拡張検出器が用いられ得る。拡張検出器列における非計測データは、境界の検出器例における計測データを重複させることによって、又は外挿によって得られる。図１９Ｂは、仮想的な拡張検出器を用いて得られた、拡張されたＦＯＶを示す。ＦＤＫでのＦＯＶの外側にある拡張領域における再構成は、ＦＤＫでのＦＯＶの内側よりも正確さに欠けることを明記しておく。
図２０の（Ａ）〜（Ｆ）は、直線状スキャンコリメーションを用いる直線状スキャンの照射線量を低減させる方法における複数の線図である。各線図において、ｈは直線状スキャン経路上のｚ座標を示し、影付けされた領域は直線状スキャンにおけるバックプロジェクションステップにおいて用いられたデータを示す。さらに、（Ｂ）〜（Ｆ）は、ｈ＞０（直線状スキャンの上部）の場合のみを示しているが、ｈ＜０の場合は対称的であり、対応する方法で扱われうる。このように、本実施形態によれば、直線状スキャンに患者によって受けられる放射線量は、影付けされた領域の外側の放射を避けるため、動的なコリメーションを用いて低減される。
図２１は、動的コリメーションに対するコリメーション窓の一例を示す。この例によれば、方程式（９）、（１０）、（１３）、（１５）を用いて我々は、垂直トップ（vertical top）とボトム境界(bottom boundaries)とのために、νＴＯＰ（ｈ）とνＢＯＴ（ｈ）とをそれぞれ得る。
Figure 2010510833
さらに、方程式（２１）は、ｈ＞５．８１ｍｍの時におけるνＴＯＰを計算するためにのみ用いられる。ｈ＜５．８１ｍｍの時、νＴＯＰ＝Ｗである。
図２２は、本発明の係るジオメトリのさらなる図である。本発明によれば、以下のパラメータは、システムオペレータによってスキャン条件に応じて決定される。
Figure 2010510833
Ｚｓｔａｒｔ 直線状スキャンの開始ｚ位置
Ｚｅｎｄ 直線状スキャンの終了ｚ位置
Ｚｃｉｒｃ 円のｚ位置
ＴｏｔａｌＶｉｅｗｓ 直線ビューの全数
以下のパラメータは、本発明の実施形態に応じて計算される。
Figure 2010510833
Ｎｖｉｅｗ 再構成のためのビュー数
ＶｉｅｗＪｕｍｐ 再構成のための最初のビュー
Δｈ 直線ビューピッチ
ｈ 直線ビュー位置、円に相対的
さらに
Figure 2010510833
又は、
Figure 2010510833
によれば、方程式（２７）が満足される。
Figure 2010510833
そして、そのような定義は、Ｚｅｎｄ＞Ｚｓｔａｒｔ、又はＺｅｎｄ＜Ｚｓｔａｒｔの場合、独立である。さらに、
Figure 2010510833
図２３Ａは、直線状スキャンデータにおけるフィルタラインのグラフである。
図２３Ｂは、直線状スキャンデータにおけるフィルタラインの数のグラフである。詳細には、この例によれば、フィルタリングラインの数は、直線状部分ＡＣ＝Ｗ−ｈ／２の長さに比例する。好ましくは、フィルタリングラインＮｆｌの数は、以下のように計算される。
Figure 2010510833
図２３Ｂは、フィルタラインＮｆｌの数のグラフである。
図２４は、本発明の他の実施形態に係る直線及び円状アルゴリズムの方法のブロック図である。この場合において、収集された直線状データにおけるサンプル（チャンネル）の数はファクタＫＤＳによって低減されうる。ここで、例としてダウンサンプリング係数ＫＤＳが２、３、４、又は５に一致しうる。
図２５Ａと図２５Ｂとは、ショートスキャンコンボルーションにおいて用いられるカーブの２つのファミリーの一例を示す。図２５Ｃと図２５Ｄとは、カーブの２つのファミリーのさらなる表現を示す図である。フィルタリング平面（又は湾曲）は、フィルタリング平面と検出器面（平坦又は円筒形。検出器のタイプによる）との交差によって得られる。特に、図２５Ａに示されるカーブは、図２５Ｃにおける弧Ａに接するフィルタリング平面に対応し、図２５Ｂにおけるカーブは、図２５Ｄにおける弧Ｂに接するフィルタリング平面に対応する。
本発明の実施形態によれば、ショートスキャン（これは円弧の部分のみに沿ってスキャンする）は、円状スキャンの一部が用いられることに依存する、ラインの１つのファミリーのみを用いる（すなわち、図２５Ａ又は図２５Ｂの１つ）。すなわち、各検出器ピクセルに対して、１つのフィルタリングラインのみが存在する。
図２６Ａは、複数のファミリー（すなわち、図２６Ａにおける全ライン）からのラインが用いられる場合における、フルスキャンコンボルーションの一例を示す図である。このファミリーは、検出器を２回覆う。すなわち、各検出器ピクセルに対して、２つのフィルタリングラインが存在する。各フィルタリングに対応する値を得るために、２つの値の差分を得る。
Figure 2010510833
図２６Ｂは、図２６Ａの２つのファミリーのラインの重なりの１次元図である。
図２７Ａと図２７Ｂとは、本発明の実施形態に係る逆リビニングの一例を示す。図２７Ａは、チャンネル１０に対する逆リビニングカーブ２７０、チャンネル４４７（例えば、センターチャンネル）に対する逆リビニングカーブ２７１、及びチャンネル８８７に対する逆リビニングカーブ２７２示す。
図２７Ｂは、逆リビニング処理の一例を示す。この例によれば、初めに逆リビニングカーブ２７４とセグメント座標Ｓｅｇとの交点のρＡとρＢとを得る。次に、これらρＡとρＢとにおけるＣＯＮＶ値を内挿し、差分をとる。
Figure 2010510833
図２８Ａは、本発明の実施形態に係る、拡張再構成関心ボリューム（ＶＯＩ）の一例を示す。物理検出器２８２は、領域２８４に対応するデータを受信する。しかしながら、ｚ軸に関する長さ２Ｗを有する仮想的に拡張された検出器２８６は、領域２８８に対応するデータを受信する。このように、表示ボリューム２９０内において、ＦＤＫ―ＶＯＩ２９２は物理検出器に対応し、ＦＤＫ―ＶＯＩ２９４は仮想的に拡張された検出器に対応する。図２８Ｂは、分かりやすさのため、ＦＤＫ―ＶＯＩ２９２とＦＤＫ―ＶＯＩ２９４との形状の別々の線図である。
複数の画像、大きいスキャン対象物の差分部分のそれぞれは、発明的なシステムによって連続して獲得されうり、全体的な対象物の１つの画像を生成するために組合わされうる。特に、複数の円状スキャンからのデータは、１つ又はそれ以上の直線状スキャンからのデータに組合わされうる。
図２９Ａ〜Ｆは、１つ又はそれ以上の直線状スキャンを有する複数の円状スキャンからのデータの結合例を示す。図２９Ａ〜Ｆのそれぞれにおいて、各分離された水平線３００は、単一の直線状スキャンに対応し、各円３０１は、単一の円状スキャン経路（これは、この経路を通過する１つ以上の円状スキャンを含みうる）。さらに、様々な直線及び円状スキャンがあらゆる適当な順番において実行されうる。
それに加え、前処理された円状データは、受信された円状投影データ中の一重に計測された及び二重に計測されたレイに対して別個の重みが与えられうる。例として、完全回転以下のデータを用いた場合、（パーカー（Parker）、シルバー（Silver））［１７］、［１８］により記述されたような冗長性重みづけが用いられうる。
ｇ（β，γ，α）を、２つの角（γ，α）によって与えられる方向での線源位置ａ（β）からのレイに沿う線積分を表すものとする。ここで、γはファン角を、αはコーン角を表す。角γ及びαは、物理検出器によって制限される。すなわち、−γｍ≦γ≦γｍ、及び−αｍ≦α≦αｍである。例えば、幾つかのスキャナは、−３０°＜γ＜３０°、−６°＜α＜６°を有する。ファンビームジオメトリの場合には、ｇ（β，γ）＝ｇ（β＋π＋２γ，−γ）であることを明記しておく。従って、我々が全回転分のファンビームデータを得る場合（いわゆるフルスキャン）、各レイは、２回計数される。そして我々は、ｗ（β，γ）＝１／２でデータを重みづけしうる。画像を再構成するために、全回転に亘ってスキャンする必要はない。ｇ（β，γ）＝ｇ（β＋π＋２γ，−γ）という関係から、我々は、π＋２γｍの再構成範囲のみが厳密なファンビーム再構成に十分であることが分かる。ここでγｍは、検出器が取り得る最大ファン角である。
パーカー（Parker）［１７］において、不連続性ができるだけ均一に分布されるように、
最小な完全データセット（π＋２γｍ）においてデータを重みづけすることが提案されている。そして以下の重みづけ関数が提案された。
Figure 2010510833
ここから、定義されていない場合、重みはゼロであるとする。例えば、β＜０、あるいはβ＞π＋２γｍの場合、ｗＰ（β，γ）＝０である。仮想的にγｍを増加させることによって、我々は、より大きい再構成範囲を得られ、従って良好なノイズ特性が得られることも明記しておく（シルバー（Silver））［１８］。ｗＰ（β，γ）に対する方程式において、物理的な最大ファン角γｍを仮想的なΓ（Γ≧γｍ）で置き換えることにより、我々は、他の重み関数を得る。
Figure 2010510833
パーカー（Parker）の重み関数は、ＭＨＳ重みづけの特殊な事例であり、従ってｗＰ（β，γ）の代わりにｗＭＨＳ（β，γ）を用いることが好まれる。
他の重み関数がノー（Noo）によって提案された［１９］。これは、任意の再構成範囲Λ＝（β０，β１）を許可する。ここでβ０及びβ１はそれぞれ、再構成角範囲の開始及び終了点である。これは、
Figure 2010510833
によって与えられる。ここで、
Figure 2010510833
である。そして関数ｃ（β）は、
Figure 2010510833
によって与えられる。ここで、Δβは、所定数（２０°）、又は全体の再構成範囲（β０−β１）の割合（但し５０％を超えない）として選ばれ得る平滑化間隔である。
本発明の多数の改良及び変形が上記の技術に基づいて可能である。従って、添付された請求項の範囲内において、ここに特別に記述された以外のやり方で本発明が実施可能であるということが認識されるべきである。

Claims (33)

  1. 対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、前記受信された円状投影データから前処理された円状投影データを生成し、前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、前記受信された直線状投影データから前処理された直線状投影データを生成するデータ収集部と、
    前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成する円状投影データ再構成部と、
    前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成する直線状投影データ再構成部と、
    前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成する画像ボリューム処理部と、
    を具備するボリューム画像再構成装置。
  2. 前記直線状経路は、前記円状経路の平面に対して垂直でない、請求項１記載の装置。
  3. 前記円状投影データ再構成部は、バックプロジェクション部と方向性フィルタリング部とを具備し、
    前記バックプロジェクション部は、フィルタリングされた検出器空間データに基づいて前記再構成された円状経路ボリューム画像を生成し、
    前記方向性フィルタリング部は、
    前記前処理された円状投影データを、予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、リビニングされたデータを生成するリサンプリング部と、
    前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、
    前記フィルタリングされたデータを検出器格子点に逆リビニングし、前記フィルタリングされた検出器空間データを生成する逆リビニング部と、を備える、
    請求項１記載の装置。
  4. 前記円状投影データ再構成部は、
    前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、
    前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するフィルタリング部と、
    前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
    請求項１記載の装置。
  5. 前記円状投影データ再構成部は、
    前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、プレ重みづけされたデータを得るプレ重みづけ部と、
    前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、フィルタリングされたデータを生成するハイブリッドフィルタリング部と、
    前記フィルタリングされたデータをバックプロジェクションし、前記再構成された円状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
    請求項１記載の装置。
  6. 前記データ収集部は、さらに、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された円状投影データを受信し、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成する、請求項１記載の装置。
  7. 前記直線状投影データ再構成部は、
    前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、
    前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
    前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
    前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、
    前記フィルタリングされた検出器空間データをバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
    請求項１記載の装置。
  8. 前記直線状投影データ再構成部は、
    前記前処理された直線状投影データを、２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、ダウンサンプリングされたデータを生成するダウンサンプリング部と、
    前記ダウンサンプリングされたデータの微分係数を計算し、微分係数データを生成する微分器と、
    前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
    前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
    前記フィルタリングされたデータを検出器格子点に逆リビニングし、フィルタリングされた検出器空間データを生成する逆リビニング部と、
    前記フィルタリングされた検出器空間データをバックプロジェクションし、ダウンサンプリングされた直線状経路ボリューム画像を生成するバックプロジェクション部と、
    前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングし、前記再構成された直線状経路ボリューム画像を生成するアップサンプリング部と、を備える、
    請求項１記載の装置。
  9. 前記対象物に電流に対応する放射線を照射するものであって、前記電流は２０ｍＡよりも低い又は等しいＸ線管と、
    前記対象物が前記放射線に照射された場合、前記直線状投影データを生成する検出器と、をさらに備える、
    請求項１記載の装置。
  10. 前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整する寝台不正確さ補償部をさらに備え、
    前記直線状投影データ再構成部は、さらに前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成する、
    請求項１記載の装置。
  11. 前記直線状投影データ再構成部は、
    前記検出された直線状投影データの微分係数を計算し、微分係数データを生成する微分器と、
    前記微分係数データをフィルタリングカーブ上にリビニングし、リビニングされたデータを生成するリビニング部と、
    前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、フィルタリングされたデータを生成するヒルベルトフィルタリング部と、
    前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションし、前記再構成された直線状経路ボリューム画像を生成するバックプロジェクション部と、を備える、
    請求項１記載の装置。
  12. 前記対象物を放射線で制御可能に照射する線源と、
    前記対象物が前記放射線で照射された場合、前記直線状投影データを生成する検出器と、
    前記線源を制御して前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射する動的コリメーション部と、をさらに備える、
    請求項１記載の装置。
  13. 前記データ収集部は、さらに、異なる捕獲時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる捕獲時間にそれぞれ対応する複数の前処理された円状投影データを生成し、
    前記画像ボリューム処理部は、さらに、前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成するものであり、
    前記対象物の前記複数のボリューム画像を前記画像ボリューム処理部から統合し、前記複数のボリューム画像の統合されたデータを生成する統合部と、
    前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成する合成部と、をさらに備える、
    請求項１記載の装置。
  14. 前記対象物は繰り返し位相を示すものであり、
    前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、
    前記合成部は、前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成する、
    請求項１３記載の装置。
  15. 前記対象物は繰り返し位相を示すものであり、
    前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定する対象物繰り返しモード決定部をさらに備え、
    前記データ収集部は、さらに、異なる獲得時間に円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成する、
    請求項１記載の装置。
  16. 前記データ収集部は、さらに、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、
    前記画像ボリューム処理部は、さらに、前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、前記対象物の前記部分の前記複数のボリューム画像を合成し、前記対象物の前記ボリューム画像を生成する、
    請求項１記載の装置。
  17. 対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、
    前記受信された円状投影データから前処理された円状投影データを生成し、
    前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、
    前記受信された直線状投影データから前処理された直線状投影データを生成し、
    前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、
    前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、
    前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、
    を特徴とする対象物のボリューム画像を再構成する方法。
  18. 前記直線状経路は、前記円状経路の平面に対して垂直でないこと、を特徴とする請求項１７記載の方法。
  19. 前記再構成された円状経路ボリューム画像を生成することは、さらに
    リビニングされたデータを生成するために、前記前処理された円状投影データを予め特定されたフィルタリング方向が表すフィルタリングカーブ上にリサンプリングし、
    フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿って前記ランプフィルタでフィルタリングし、
    前記フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
    前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションすること、
    を特徴とする請求項１７記載の方法。
  20. 前記再構成された円状経路ボリューム画像を生成することは、さらに
    プレ重みづけされたデータを得るために、前記検出器のファン角のコサインと前記検出器のコーン角のコサインとで前記前処理された円状投影データをプレ重みづけし、
    フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをフェルドカンプ、デイビス、クレス（ＦＤＫ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってフィルタリングし、
    前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションすること、
    を特徴とする請求項１７記載の方法。
  21. 前記再構成された円状経路ボリューム画像を生成することは、さらに、
    プレ重みづけされたデータを得るために、前記前処理されたデータを前記検出器のコーン角のコサインでプレ重みづけし、
    フィルタリングされたデータを生成するために、前記プレ重みづけされたデータをザミャチン、タグチ、シルバー（ＺＴＳ）ランプフィルタで前記検出器の列に沿って又は予め特定されたフィルタリングカーブに沿ってハイブリッドフィルタリングし、
    前記再構成された円状経路ボリューム画像を生成するために、前記フィルタリングされたデータをバックプロジェクションすること、
    を特徴とする請求項１７記載の方法。
  22. 前記円状投影データを受信することは、完全円の部分を備える前記円状経路に沿って前記検出器によって収集された前記円状投影データを受信し、
    前記前処理された円状データを生成することは、前記受信された円状投影データにおける一重に計測された及び二重に計測されたレイに対して異なる重みづけを有する前記前処理された円状データを生成すること、
    を特徴とする請求項１７記載の方法。
  23. 前記再構成された直線状経路ボリューム画像を生成することは、さらに、
    微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、
    リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
    フィルタリングされたデータを生成するために、前記リビンニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、
    フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
    前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションすること、
    を特徴とする請求項１７記載の方法。
  24. さらに、
    ダウンサンプリングされたデータを生成するために、前記前処理された直線状投影データを２、３、４、又は５から選択されたサンプリング係数によってダウンサンプリングし、
    微分係数データを生成するために、前記ダウンサンプリングされたデータの微分係数を計算し、
    リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
    フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブ沿ってヒルベルトフィルタリングし、
    フィルタリングされた検出器空間データを生成するために、前記フィルタリングされたデータを検出器格子点に逆リビニングし、
    ダウンサンプリングされた直線状経路ボリューム画像を生成するために、前記フィルタリングされた検出器空間データをバックプロジェクションし、
    前記再構成された直線状経路ボリューム画像を生成するために、前記ダウンサンプリングされた直線状経路ボリューム画像を前記サンプリング係数によってアップサンプリングすること、
    を特徴とする請求項１７記載の方法。
  25. さらに、
    前記対象物に、Ｘ線管内の電流であって２０ｍＡよりも低い又は等しい前記電流に対応する放射線を照射し、
    前記対象物が前記放射線に照射された場合、前記直線状投影データを生成すること、
    を特徴とする請求項１７記載の方法。
  26. さらに、
    前記前処理された円状投影データと前記前処理された直線状投影データとの間の相関に基づいて再構成パラメータを調整することによって寝台不正確さ補償し、
    前記再構成パラメータに基づいて前記対象物の前記再構成された直線状経路ボリューム画像を生成すること、
    を特徴とする請求項１７記載の方法。
  27. 前記直線状経路ボリューム画像を生成することは、さらに、
    微分係数データを生成するために、前記検出された直線状投影データの微分係数を計算し、
    リビニングされたデータを生成するために、前記微分係数データをフィルタリングカーブ上にリビニングし、
    フィルタリングされたデータを生成するために、前記リビニングされたデータを前記フィルタリングカーブに沿ってヒルベルトフィルタリングし、
    前記再構成された直線状経路ボリューム画像を生成するために、前記フィルタリングされたデータを直接的に前記フィルタリングカーブからバックプロジェクションすること、
    を特徴とする請求項１７記載の方法。
  28. さらに、
    前記対象物を放射線で制御可能に照射し、
    前記対象物が前記放射線で照射された場合、前記直線状投影データを生成し、
    前記対象物に、フィルタリング部によってフィルタリングされるデータを受信する前記検出器の一部分に対応する照射開口で前記放射線を照射するように前記照射することを制御すること、
    を特徴とする請求項１７記載の方法。
  29. 前記円状投影データを前記受信することは、異なる獲得時間に前記円状経路の複数部分に沿って前記検出器によって収集された前記円状投影データを受信し、前記前処理された円状投影データを前記生成することは、異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、前記結合することは、前記前記複数の前処理された円状投影データに基づいて前記対象物の複数のボリューム画像を生成することであり、
    さらに、
    前記対象物の前記複数のボリューム画像を統合し、
    前記複数のボリューム画像の統合されたデータを生成し、
    前記複数のボリューム画像の前記統合されたデータから前記対象物の合成ボリューム画像を生成すること、
    を特徴とする請求項１７記載の方法。
  30. 前記対象物は繰り返し位相を示すものであり、
    さらに、
    前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、
    前記対象物が前記関心位相を示した場合、前記位相時間に対応する複数のボリューム画像のうちの少なくとも１つに基づいて前記合成ボリューム画像を生成すること、
    を特徴とする請求項２９記載の方法。
  31. 前記対象物は繰り返し位相を示すものであり、
    さらに、
    前記対象物が前記繰り返し位相中の関心位相を示した場合、位相時間を同定し、
    前記円状経路の複数部分に沿って、異なる獲得時間に前記検出器によって収集された前記円状投影データを受信し、
    異なる獲得時間にそれぞれ対応する複数の前処理された円状投影データを生成し、
    前記対象物が前記関心位相を示した場合、時間間隔に対応する前記複数の前処理された円状投影データから前記前処理された円状投影データを生成すること、
    を特徴とする請求項１７記載の方法。
  32. 前記円状投影データを前記受信することは、前記検出器によって複数の円状経路に沿って収集された前記円状投影データを受信し、
    前記前処理された円状投影データを前記生成することは、前記複数の円状経路中の異なる円状経路にそれぞれ対応する複数の前処理された円状投影データを生成し、
    前記対象物の部分の複数のボリューム画像を、前記複数の前処理された円状投影データにそれぞれ対応する複数の再構成された円状経路ボリューム画像に基づいて生成し、
    前記合成することは、前記対象物の前記ボリューム画像を生成するために、前記対象物の前記部分の前記複数のボリューム画像を合成する、
    ことを特徴とする請求項１７記載の方法。
  33. コンピュータにより実行された場合に、
    対象物に関する円状経路に沿って検出器によって収集された円状投影データを受信し、
    前記受信された円状投影データから前処理された円状投影データを生成し、
    前記対象物に関する直線状経路に沿って前記検出器によって収集された直線状投影データを受信し、
    前記受信された直線状投影データから前処理された直線状投影データを生成し、
    前記前処理された円状投影データから、ランプフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された円状経路ボリューム画像を生成し、
    前記前処理された直線状投影データから、ヒルベルトフィルタを含む再構成アルゴリズムを用いて前記対象物に関する再構成された直線状経路ボリューム画像を生成し、
    前記再構成された円状経路ボリューム画像と前記再構成された直線状経路ボリューム画像とを結合して、前記対象物のボリューム画像を生成すること、
    を具備するステップをコンピュータに実行させるコンピュータプログラムインストラクションを記憶するコンピュータ読み取り可能媒体。
JP2009538528A 2006-11-24 2007-11-26 円状及び直線状軌道を用いた医用撮像における断層撮影法の再構成のための方法及びシステム Active JP5161888B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US86716506P 2006-11-24 2006-11-24
US60/867,165 2006-11-24
PCT/US2007/085510 WO2008064367A2 (en) 2006-11-24 2007-11-26 Method and system for tomographic reconstruction in medical imaging using the circle and line trajectory

Publications (2)

Publication Number Publication Date
JP2010510833A true JP2010510833A (ja) 2010-04-08
JP5161888B2 JP5161888B2 (ja) 2013-03-13

Family

ID=39430622

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009538528A Active JP5161888B2 (ja) 2006-11-24 2007-11-26 円状及び直線状軌道を用いた医用撮像における断層撮影法の再構成のための方法及びシステム

Country Status (3)

Country Link
US (1) US8270559B2 (ja)
JP (1) JP5161888B2 (ja)
WO (1) WO2008064367A2 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019031793A1 (ko) * 2017-08-07 2019-02-14 삼성전자주식회사 워크 스테이션 이를 포함하는 의료영상 장치, 및 그 제어방법
US11158093B2 (en) 2017-08-07 2021-10-26 Samsung Electronics Co., Ltd. Workstation, medical imaging device including same, and control method therefor

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5641748B2 (ja) * 2009-03-18 2014-12-17 株式会社東芝 Ｘ線コンピュータ断層撮影装置及びデータ処理方法
US8559687B2 (en) * 2009-10-30 2013-10-15 Kabushiki Kaisha Toshiba Method for circular scan reconstruction in computed tomography and computed tomographic device
EP2330559B9 (en) * 2009-10-30 2013-11-06 Kabushiki Kaisha Toshiba Computed tomographic device
TW201215126A (en) * 2010-09-27 2012-04-01 Hon Hai Prec Ind Co Ltd Image dividing system for cameras and using method of the same
DE102011076547A1 (de) * 2011-05-26 2012-11-29 Siemens Aktiengesellschaft Verfahren zum Gewinnen eines 3D-Bilddatensatzes zu einem Bildobjekt
US9237873B2 (en) 2011-10-31 2016-01-19 General Electric Company Methods and systems for CT projection domain extrapolation
US8995735B2 (en) * 2012-05-23 2015-03-31 General Electric Company System and method for wide cone helical image reconstruction using blending of two reconstructions
US9737279B2 (en) 2012-08-31 2017-08-22 Toshiba Medical Systems Corporation X-ray CT apparatus
US8977858B1 (en) 2014-05-27 2015-03-10 Support Intelligence, Inc. Using space-filling curves to fingerprint data
US9558569B2 (en) * 2014-08-12 2017-01-31 Toshiba Medical Systems Corporation Method and system for substantially reducing cone beam artifacts based upon image domain differentiation in circular computer tomography (CT)
CN108292428B (zh) 2015-12-11 2023-03-31 上海联影医疗科技股份有限公司 图像重建的系统和方法
US10255696B2 (en) 2015-12-11 2019-04-09 Shanghai United Imaging Healthcare Co., Ltd. System and method for image reconstruction
US10249066B2 (en) * 2017-04-14 2019-04-02 Shanghai United Imaging Healthcare Co., Ltd. System and method for MPR streak reduction
US10546143B1 (en) 2017-08-10 2020-01-28 Support Intelligence, Inc. System and method for clustering files and assigning a maliciousness property based on clustering

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001327514A (ja) * 2000-05-23 2001-11-27 Toshiba Corp 放射線治療計画装置及び放射線治療計画法
JP2002282246A (ja) * 2001-01-30 2002-10-02 Toshiba Corp Ｘ線ｃｔ装置、ボリューム画像化方法、及び当該ボリューム画像化方法を実現するプログラムを格納するコンピュータが読みとり可能な記録媒体
JP2003144428A (ja) * 2001-11-13 2003-05-20 Koninkl Philips Electronics Nv Ｘ線透視コンピュータ断層撮影方法
JP2003199740A (ja) * 2001-10-25 2003-07-15 Toshiba Corp Ｘ線ｃｔ装置、ｘ線ｃｔの３次元再構成方法、及びｘ線ｃｔの重み設定方法
JP2005296469A (ja) * 2004-04-14 2005-10-27 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置及びその制御方法
US20060067457A1 (en) * 2004-09-29 2006-03-30 Kabushiki Kaisha Toshiba Image reconstruction method using hilbert transform
JP2006087921A (ja) * 2004-09-21 2006-04-06 General Electric Co <Ge> 対象領域情報を用いて連続的に多重解像度３次元画像を再構成する方法およびシステム
JP2006288472A (ja) * 2005-04-06 2006-10-26 Toshiba Corp Ｘ線ｃｔ装置及び画像再構成方法
WO2007119124A1 (en) * 2006-04-19 2007-10-25 Philips Intellectual Property & Standards Gmbh Cone beam computed tomography with multiple partial scanning trajectories
JP2008512145A (ja) * 2004-09-03 2008-04-24 ジャック イー ジュニ 核医学画像装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5170439A (en) * 1991-06-11 1992-12-08 Picker International, Inc. Cone beam reconstruction using combined circle and line orbits
US5784481A (en) * 1996-06-25 1998-07-21 General Electric Company CT cone beam image reconstruction with circle and line scan path
US6014419A (en) * 1997-11-07 2000-01-11 Hu; Hui CT cone beam scanner with fast and complete data acquistion and accurate and efficient regional reconstruction
US6148056A (en) * 1999-02-22 2000-11-14 General Electric Company Efficient cone-beam reconstruction system using circle-and-line orbit data
US6330298B1 (en) 2000-10-17 2001-12-11 Siemens Corporate Research, Inc. Method and apparatus for simplifying the correction of image inaccuracies caused by processing of masked cone beam projection data
US7197105B2 (en) * 2001-08-16 2007-03-27 Research Foundation Of The University Of Central Florida, Inc. Efficient image reconstruction algorithm for the circle and line cone beam computed tomography
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
US6771732B2 (en) 2002-02-28 2004-08-03 The Board Of Trustees Of The University Of Illinois Methods and apparatus for fast divergent beam tomography
US6917663B2 (en) 2003-06-16 2005-07-12 Kabushiki Kaisha Toshiba Cone-beam reconstruction apparatus and computed tomography apparatus
US7372937B2 (en) 2004-07-16 2008-05-13 University Of Iowa Research Foundation Systems and methods of non-standard spiral cone-beam computed tomograpy (CT)
US7359478B2 (en) * 2004-11-18 2008-04-15 Toshiba Medical Systems Corporation Method for restoring truncated helical cone-beam computed tomography data

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001327514A (ja) * 2000-05-23 2001-11-27 Toshiba Corp 放射線治療計画装置及び放射線治療計画法
JP2002282246A (ja) * 2001-01-30 2002-10-02 Toshiba Corp Ｘ線ｃｔ装置、ボリューム画像化方法、及び当該ボリューム画像化方法を実現するプログラムを格納するコンピュータが読みとり可能な記録媒体
JP2003199740A (ja) * 2001-10-25 2003-07-15 Toshiba Corp Ｘ線ｃｔ装置、ｘ線ｃｔの３次元再構成方法、及びｘ線ｃｔの重み設定方法
JP2003144428A (ja) * 2001-11-13 2003-05-20 Koninkl Philips Electronics Nv Ｘ線透視コンピュータ断層撮影方法
JP2005296469A (ja) * 2004-04-14 2005-10-27 Ge Medical Systems Global Technology Co Llc Ｘ線ｃｔ装置及びその制御方法
JP2008512145A (ja) * 2004-09-03 2008-04-24 ジャック イー ジュニ 核医学画像装置
JP2006087921A (ja) * 2004-09-21 2006-04-06 General Electric Co <Ge> 対象領域情報を用いて連続的に多重解像度３次元画像を再構成する方法およびシステム
US20060067457A1 (en) * 2004-09-29 2006-03-30 Kabushiki Kaisha Toshiba Image reconstruction method using hilbert transform
JP2006095297A (ja) * 2004-09-29 2006-04-13 Toshiba Corp スキャン対象に関するｃｔ画像内の再構成点における画像データ値を決定する再構成方法及びｘ線コンピュータ断層撮影装置
JP2006288472A (ja) * 2005-04-06 2006-10-26 Toshiba Corp Ｘ線ｃｔ装置及び画像再構成方法
WO2007119124A1 (en) * 2006-04-19 2007-10-25 Philips Intellectual Property & Standards Gmbh Cone beam computed tomography with multiple partial scanning trajectories

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019031793A1 (ko) * 2017-08-07 2019-02-14 삼성전자주식회사 워크 스테이션 이를 포함하는 의료영상 장치, 및 그 제어방법
US11158093B2 (en) 2017-08-07 2021-10-26 Samsung Electronics Co., Ltd. Workstation, medical imaging device including same, and control method therefor

Also Published As

Publication number Publication date
JP5161888B2 (ja) 2013-03-13
WO2008064367A3 (en) 2008-07-31
US8270559B2 (en) 2012-09-18
US20100283779A1 (en) 2010-11-11
WO2008064367A2 (en) 2008-05-29

Similar Documents

Publication Publication Date Title
JP5161888B2 (ja) 円状及び直線状軌道を用いた医用撮像における断層撮影法の再構成のための方法及びシステム
JP6214128B2 (ja) 画像処理装置、画像処理方法、及び記憶媒体
JP5019193B2 (ja) スキャン対象に関するｃｔ画像内の再構成点における画像データ値を決定する再構成方法及びｘ線コンピュータ断層撮影装置
US8705828B2 (en) Methods and apparatus for super resolution scanning for CBCT system and cone-beam image reconstruction
US7142633B2 (en) Enhanced X-ray imaging system and method
EP1741062B1 (en) Cone beam ct apparatus using truncated projections and a previously acquired 3d ct image
US7706497B2 (en) Methods and apparatus for noise estimation for multi-resolution anisotropic diffusion filtering
JP5199081B2 (ja) 心臓ｃｔ撮影のバンドアーチファクトの抑制
JP4740516B2 (ja) コーン・傾斜平行式のサンプリング及び再構成の方法及び装置
US10789738B2 (en) Method and apparatus to reduce artifacts in a computed-tomography (CT) image by iterative reconstruction (IR) using a cost function with a de-emphasis operator
US20090297009A1 (en) Method of reconstructing an image function from radon data
EP2030170B1 (en) Dynamic computed tomography imaging
WO2002030282A9 (en) System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit
JP2007527253A (ja) 冗長データを含むヘリカルコーンビーム・コンピュータ断層撮像法に対する正確な画像再構成システム及び方法
JP3987024B2 (ja) 横方向のフィルタリング処理を用いたトモシンセシス画像を強調する方法及びシステム
JP2007000408A (ja) Ｘ線ｃｔ装置
JP3290726B2 (ja) 透過型三次元断層撮影装置
US7027552B2 (en) High resolution CT scanner
JP4509255B2 (ja) 透視画像作成方法及び装置
US6647084B1 (en) Method and apparatus for filtering projection data of a helical scan
US20080086052A1 (en) Methods and apparatus for motion compensation
CN117425918A (zh) 用于ct重建的方法
JP4282110B2 (ja) Ｘ線コンピュータ断層撮影装置
Fuchs et al. Fast volume scanning approaches by X-ray-computed tomography
JP2009106759A (ja) コンピュータ断層撮影装置、その処理方法及び記録媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120731

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120926

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121214

R150 Certificate of patent or registration of utility model

Ref document number: 5161888

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

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350