JP2006000226A - X線ct装置 - Google Patents
X線ct装置 Download PDFInfo
- Publication number
- JP2006000226A JP2006000226A JP2004177406A JP2004177406A JP2006000226A JP 2006000226 A JP2006000226 A JP 2006000226A JP 2004177406 A JP2004177406 A JP 2004177406A JP 2004177406 A JP2004177406 A JP 2004177406A JP 2006000226 A JP2006000226 A JP 2006000226A
- Authority
- JP
- Japan
- Prior art keywords
- reconstruction
- range
- ray
- superabsorber
- distance
- 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.)
- Pending
Links
- 239000006096 absorbing agent Substances 0.000 claims description 20
- 230000000694 effects Effects 0.000 claims 1
- 230000000116 mitigating effect Effects 0.000 claims 1
- 239000002184 metal Substances 0.000 abstract description 8
- 229910052751 metal Inorganic materials 0.000 abstract description 8
- 239000002250 absorbent Substances 0.000 abstract description 6
- 230000002745 absorbent Effects 0.000 abstract description 5
- 230000007423 decrease Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 23
- 238000003384 imaging method Methods 0.000 description 16
- 238000002591 computed tomography Methods 0.000 description 13
- 238000003745 diagnosis Methods 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 230000001603 reducing effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
【課題】金属などのX線高吸収体があってもアーチファクトの少ない断層画像を得る。
【解決手段】 再構成関数は高吸収体Mが含まれる範囲R1内では、図5に示すような高周波成分を低下させる再構成関数を用い、通常の範囲R3ではShepp&Logan関数のような標準的な再構成関数を用いる。移行範囲R2においては、高吸収体Mが含まれる範囲からの距離に応じて、高吸収体範囲の再構成関数と通常範囲の再構成関数を補間して漸次移行する再構成関数を用いる。
【選択図】図4
【解決手段】 再構成関数は高吸収体Mが含まれる範囲R1内では、図5に示すような高周波成分を低下させる再構成関数を用い、通常の範囲R3ではShepp&Logan関数のような標準的な再構成関数を用いる。移行範囲R2においては、高吸収体Mが含まれる範囲からの距離に応じて、高吸収体範囲の再構成関数と通常範囲の再構成関数を補間して漸次移行する再構成関数を用いる。
【選択図】図4
Description
本発明は、X線高吸収体によるアーチファクトを軽減した断層画像を得るコーンビーム型のX線CT装置に関するものである。
X線CT(Computer Tomography)装置においては、被検体の内部に金属等のX線高吸収体が存在すると、X線強度の減弱に伴う高周波成分が投影データに含まれるようになり、再構成画像にストリークアーチファクトと呼ばれる直線状の偽像が生じ、診断の妨げとなる。
このようなアーチファクトを除去或いは低減するため、例えば特許文献1では、先ず投影データを再構成して断層画像を作成し、この断層画像から金属などに関する関心領域を設定する。そして、投影データ上における関心領域に対応する領域を求め、この領域について移動平均処理を行う。この移動平均処理された投影データを再構成することによって、アーチファクトの低減された断層画像を得る。
コーンビーム型X線CT装置などで、1人の患者に対する複数の断層画像を再構成する場合に同一の再構成関数が用いられる。このため、体内外に金属を有する患者を撮影して断層画像を作成した場合に、金属が存在する断面の断層画像は、アーチファクトが発生して診断に不向きな画像となってしまう。
また、このアーチファクトを除去するための前述の特許文献1の方法では、金属等の高吸収体の存在或いはその存在領域を断層画像より発見するために、最初に再構成処理を行わなければならない。そして、投影データを補正した後に、再度再構成処理を行うことになる。
つまり、2度の再構成処理をする必要があるが、再構成処理は計算量が極めて膨大で、時間のかかる処理である。特に、X線検出器に二次元平面検出器を用いたコーンビーム型CT装置の場合には、体軸方向のスライスピッチも画素ピッチと同等の細かさで再構成できるため、再構成処理にかかる時間は多大なものとなる。
また、アーチファクトを除去或いは低減するために、投影データに移動平均処理を行えば処理時間が増大し、更に移動平均処理では精度の良い補正を行うことができない。
本発明の目的は、上述の問題点を解消し、計算時間の損失がなく撮影サイクルを維持することができ、効率の良い撮影、診断に有効なX線CT装置を提供することにある。
上記目的を達成するための本発明に係るX線CT装置は、投影画像から高吸収体領域を抽出し再構成する断面が高吸収体を含むか否かを決定する決定手段と、前記高吸収体を含まない断面について通常の第1の再構成関数を使用して再構成する第1の再構成手段と、前記高吸収体を含む断面についてより高周波をカットするアーチファクト軽減効果の高い第2の再構成関数を用いて再構成する第2の再構成手段と、前記高吸収体を含む断面の近傍の断面についてその距離に応じて前記第2の再構成関数から前記第1の再構成関数に漸次移行する第3の再構成関数を用いて再構成する第3の再構成手段とを備えたことを特徴とする。
また、本発明に係るX線CT装置は、投影画像から高吸収体領域を抽出した該高吸収体領域と再構成する点との距離を求め、該距離によって逆投影するデータを求めるための範囲を決定する決定手段と、前記範囲内の投影データを補間することによって求めたデータを逆投影して再構成を行う再構成手段とを有することを特徴とする。
本発明に係るX線CT装置によれば、再構成関数を変更することによって金属アーチファクトの少ない断層画像を得ることができ、診断上有効であり、再構成関数を変更するだけであるので、計算時間の損失がなく撮影サイクルを維持することが可能となり、効率の良い撮影、診断が有効に実施できる。
また、逆投影するデータ範囲を高吸収体に合わせて可変にすることで、アーチファクトの少ない自然な断層画像を得ることができる。
本発明を図示の実施例に基づいて詳細に説明する。
図1は本発明によるX線CTシステムの構成図である。患者である被写体Sを挟んでX線源1とX線検出器2が配置され、X線源1には撮影システム制御部3の出力が、X線発生制御部4を介して接続されている。また、X線検出器2の出力は画像入力部5を介して撮影システム制御部3に接続されている。更に、撮影システム制御部3には、画像処理部6、画像保存部7、診断モニタ8、操作部9、ネットワーク10が接続されており、ネットワーク10にはプリンタ11、診断ワークステーション12、画像データベース13が接続されている。
X線発生装置制御部4により制御されたX線源1から発生したX線は、患者である被写体Sを透過してX線検出器2により検知され、検知されたX線は投影画像として画像入力部5に入力される。このX線源1とX線検出器2は被写体Sを回転中心として回転を行いながら所定の回転角度毎に投影画像の収集を行う。或いは、X線源1とX線検出器2の位置関係を保持しながら、回転テーブルに固定した被写体Sを回転させてもよい。
入力された各回転角度の投影画像は、撮影システム制御部3を介して画像処理部6によってX線検出器2の補正、対数変換を含めた前処理や再構成処理等の画像処理がなされ、断層画像群が作成される。作成された断層画像群は診断モニタ8に表示、画像保存部7に保存、ネットワーク10を介してプリンタ11、診断ワークステーション12、画像データベース13に出力される。表示のウインドウ操作や体軸方向の断層画像の切換表示操作や断面変換操作、三次元表面表示操作などの種々の操作は操作部9によって行われる。
図2はこのシステムにおけるCT撮影のフローチャート図である。最初に、ステップS201で管電圧、管電流、曝射時間などの撮影条件や患者名、年齢、性別などの患者情報や、検査IDなど検査情報が設定される。設定された撮影条件に合わせてステップS202で撮影が行われると、撮影により得られた投影画像に対してステップS203で前処理が実行される。
そして、前処理された投影画像をステップS204で再構成処理することによって断層画像が作成される。作成した断層画像をステップS205で表示して、ステップS206でこの表示画像を確認することによって撮影の成否の判断を行う。
撮影に失敗していればステップS202に戻って再撮影を行い、撮影に成功していれば、ステップS207で作成した断層画像を画像サーバーなどに転送し、次の撮影に移ってゆく。
図3は処理のフローチャート図を示し、先ずステップS301で撮影によって得られた投影画像から高吸収体の抽出を行う。高吸収体ではX線の減衰が大きいため、投影画像上で高吸収体領域の画素値は相対的に小さな値となっている。
従って、閾値処理によって2値化して概略の形状を抽出し、2値画像処理である膨張(expansion)処理、侵食(erosion)処理を繰り返すことにより、余分な小成分や穴などを取り除くことができ、無駄なく鮮明に高吸収体領域を抽出することができる。高吸収体領域を抽出すると、ステップS302で抽出した高吸収体領域に対応するスライス位置を同定する。
これは、図4に示すように投影画像Pの縦方向がスライス方向に相当するので、被写体S中の金属などの高吸収体Mの存在する縦方向範囲を同定し、これを高吸収体Mの存在するスライス位置範囲とする。そして、ステップS303でこのスライス位置範囲の上下に移行範囲を設定する。この範囲の幅はスライス位置範囲の幅を、所定の割合で拡大又は縮小した幅を用いればよい。
このようにして、ステップS304で高吸収体Mが存在するスライス位置範囲、移行範囲、及びその他の通常範囲が決まるので、それぞれの範囲内で再構成時の再構成関数を設定する。
再構成関数は図4に示すように高吸収体Mが含まれる範囲R1内では、図5に示すような高周波成分を低下させる再構成関数を用い、通常の範囲R3ではShepp&Logan関数のような標準的な再構成関数を用いる。移行範囲R2においては、高吸収体Mが含まれる範囲からの距離に応じて、高吸収体範囲の再構成関数と通常範囲の再構成関数を補間して漸次移行する再構成関数を用いる。
このように設定された再構成関数を、図2のステップS204で再構成処理する際に、そのスライス位置に応じて再構成関数を選択して再構成する。
このようにして、高吸収体Mを含まない断面については通常通りの再構成画像が得られ、高吸収体Mを含む断面についてはアーチファクトの軽減された再構成画像が得られる。そして、高吸収体Mを含む断面から高吸収体Mを含まない断面への移行領域の断面では、断面間の変化を滑らかなものとすることができる。
実施例2では、実施例1と同様に図2のフローチャート図のステップS203の前処理とステップS204での再構成処理が実行される。図6はこの部分における実施例2のフローチャート図を示している。
ここでは、実施例1と同様に撮影によって得られた投影画像から高吸収体の抽出を行う。そして、再構成処理ではコンボリューション処理の後に、各再構成点ごとに全投影画像において対応するデータを逆投影する。この対応するデータを高吸収体との距離に応じて適応的に算出するために、ステップS602〜ステップS604を実行する。先ず、ステップS602では投影画像上における再構成点と高吸収体との距離を求める。
このために、再構成点の投影画像上の座標を求める。座標系は図7に示すように、三次元の絶対座標系を(X,Y,Z)とし、再構成する断層面を三次元絶対座標系の(X,Y)平面としてスライス方向をZ軸にとる。
X線源1の回転中心は三次元絶対座標系の原点O(0,0,0)であるとして、X線源1の回転開始の位置の座標を(−D0,0,0)であるとする。X線源1が角度θだけ回転したときのX線源1の位置の座標は(−D0 cosθ,−D0 sinθ,0)となる。そして、再構成点を(rx,ry,rZ)として、X線源1から原点Oを通って、反対側にD0+D1の個所に垂直に投影面があるとして、投影面上にO’を原点とする局所的な二次元座標系(x,y)をとる。
この投影面と局所二次元座標系における再構成点の座標を(Px,Py)とすれば、次のように現すことができる。
Px={rx−t×(D0 cosθ+rx)−D1 cos}× cosθ−{ry−t×(D0 sinθ+ry)−D1 sinθ}×sinθ
Py=rz−t
ただし、t={(rx−D1 cosθ)×D1 cosθ+(ry−D1 sinθ)×D1 sinθ}/{(rx−D0 cosθ)×D1 cosθ+(ry−D0 sinθ)×D1 sinθ}
Px={rx−t×(D0 cosθ+rx)−D1 cos}× cosθ−{ry−t×(D0 sinθ+ry)−D1 sinθ}×sinθ
Py=rz−t
ただし、t={(rx−D1 cosθ)×D1 cosθ+(ry−D1 sinθ)×D1 sinθ}/{(rx−D0 cosθ)×D1 cosθ+(ry−D0 sinθ)×D1 sinθ}
この再構成点(Px,Py)と高吸収体との距離を求める。高吸収体の領域は前述のように2値画像として求めることができ、この高吸収体領域の境界線を求める。境界線を抽出する方法は、例えば非特許文献1に記載されている画像処理アルゴリズムを用いればよい。これにより、高吸収体の境界線の点列の座標Qn(Xn,Yn)が求まる。
画像処理アルゴリズムの最新動向:高木幹雄他編(株)新技術コミュニケーションズの5.幾何学的特徴の処理(2)
再構成点と高吸収体との距離PQは、再構成点と境界線の点列との最小距離によって求める。つまり、PQ=min√{(Px−Xn)2+(Py−Yn)2)である。
このようにして、図6のステップS602で再構成点と高吸収体との距離が求まると、ステップS603でこの距離に応じた逆投影データ範囲を決定する。この決め方は、例えば画素サイズで定義された距離と逆投影データ範囲を横画素数×縦画素数で表した次の表1のテーブルを用いればよい。
表1
距離 範囲
0〜1 7×7
1〜2 6×6
2〜3 5×5
3〜4 4×4
4〜5 3×3
距離 範囲
0〜1 7×7
1〜2 6×6
2〜3 5×5
3〜4 4×4
4〜5 3×3
また、逆投影データ範囲の横画素数と縦画素数を違う大きさにすることもできる。これは例えば図8に示すように、再構成点の最近傍の高吸収体の境界線の部分線分、例えば最短距離の点を中心にした所定の数の点列が縦長であれば、逆投影データ範囲を横長にとるようにする。また、この逆に点列が横長であれば、図9に示すように逆投影データ範囲を縦長にとるようにする。
逆投影範囲が求まると、ステップS604でこの範囲内の投影画像データを補間して逆投影データを求める。補間には、スプライン関数等の補間関数を用いてもよいし、或いは再構成点からの距離の逆数を正規化した重み付けで足し込んだデータを用いることもできる。
このようにして求めた逆投影データを逆投影することによって、アーチファクトで低減された再構成画像を作成することができる。
投影画像データから高吸収体領域を抽出した後に、再構成処理の中の逆投影処理においてアーチファクトを低減する方法を説明すると、X線CT装置の再構成処理として最も一般的に用いられている手法は、コンボリューション(演算)・バックプロジェクション(逆投影演算)法である。
コーンビームCTにおいても、このコンボリューション・バックプロジェクション法を、コーンビームに拡張したFeldkampの方法が用いられている。このFeldkampの方法は、image=BP[w1×{w2=proj*func}]Λと表すことができる。
ここで、imageは再構成画像、projは投影データ、funcは再構成関数と呼ばれるコンボリューション関数データであり、w1、w2は重みを表し、*はコンボリューション、BPはバックプロジェクションを表している。
図10は従来のバックプロジェクションの説明図を示し、バックプロジェクションは再構成点Aを通るX線源1からのX線経路の投影画像P上における点のデータを1周(ハーフスキャンでは半周+ファン角度)に渡って、全て足し込むことによって行われる。ただし、X線経路の投影画像P上における点は、必ずしも投影画像Pの画素の中心点に一致するとは限らないため、その点の近傍4点を補間したデータを逆投影データとして足し込んでいる。
本発明では、この逆投影に用いるデータを全て近傍4点とするのではなく、図11に示すように投影画像P上で再構成点Aと抽出された高吸収体Mとの距離に応じて範囲を可変にするものである。つまり、高吸収体Mとの距離が近ければ、より広い範囲にわたるデータを補間して用い、高吸収体Mとの距離が遠ければ、その範囲を狭めるものである。これによって、高吸収体Mの近辺における高周波成分の影響を緩和することができ、アーチファクトが低減された再構成画像を得ることができる。
1 X線源
2 X線検出器
3 撮影システム制御部
4 X線発生制御部
5 画像入力部
6 画像処理部
7 画像保存部
8 診断モニタ
9 操作部
2 X線検出器
3 撮影システム制御部
4 X線発生制御部
5 画像入力部
6 画像処理部
7 画像保存部
8 診断モニタ
9 操作部
Claims (4)
- 投影画像から高吸収体領域を抽出し再構成する断面が高吸収体を含むか否かを決定する決定手段と、前記高吸収体を含まない断面について通常の第1の再構成関数を使用して再構成する第1の再構成手段と、前記高吸収体を含む断面についてより高周波をカットするアーチファクト軽減効果の高い第2の再構成関数を用いて再構成する第2の再構成手段と、前記高吸収体を含む断面の近傍の断面についてその距離に応じて前記第2の再構成関数から前記第1の再構成関数に漸次移行する第3の再構成関数を用いて再構成する第3の再構成手段とを備えたことを特徴とするX線CT装置。
- 投影画像から高吸収体領域を抽出した該高吸収体領域と再構成する点との距離を求め、該距離によって逆投影するデータを求めるための範囲を決定する決定手段と、前記範囲内の投影データを補間することによって求めたデータを逆投影して再構成を行う再構成手段とを有することを特徴とするX線CT装置。
- 前記決定手段は予め前記高吸収体領域再構成点間距離と逆投影データ範囲のテーブルを用いて、前記高吸収体領域再構成点間距離から前記逆投影データ範囲を決定することを特徴とする請求項2に記載のX線CT装置。
- 前記再構成点の最近傍の前記高吸収体の境界線の部分線分が縦長であれば、前記逆投影データ範囲を横長に、前記部分線分が横長であれば前記逆投影データ範囲を縦長にとることを特徴とする請求項3に記載のX線CT装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004177406A JP2006000226A (ja) | 2004-06-15 | 2004-06-15 | X線ct装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004177406A JP2006000226A (ja) | 2004-06-15 | 2004-06-15 | X線ct装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2006000226A true JP2006000226A (ja) | 2006-01-05 |
Family
ID=35769132
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2004177406A Pending JP2006000226A (ja) | 2004-06-15 | 2004-06-15 | X線ct装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2006000226A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100734884B1 (ko) * | 2006-07-11 | 2007-07-03 | 한국전자통신연구원 | Ieee 802.16/와이브로 망상의 네이버 탐색프로토콜 메시지 전송 방법 |
EP2605211A1 (en) | 2011-12-16 | 2013-06-19 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and program |
-
2004
- 2004-06-15 JP JP2004177406A patent/JP2006000226A/ja active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100734884B1 (ko) * | 2006-07-11 | 2007-07-03 | 한국전자통신연구원 | Ieee 802.16/와이브로 망상의 네이버 탐색프로토콜 메시지 전송 방법 |
EP2605211A1 (en) | 2011-12-16 | 2013-06-19 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and program |
US9014450B2 (en) | 2011-12-16 | 2015-04-21 | Canon Kabushiki Kaisha | Method and apparatus for filtering projection images |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3950782B2 (ja) | 3次元逆投影方法およびx線ct装置 | |
US8284892B2 (en) | System and method for image reconstruction | |
JP4576032B2 (ja) | 2回パス式コーン・ビーム画像再構成の方法及び装置 | |
JP5918374B2 (ja) | X線ct装置およびx線ct画像の処理方法 | |
CN103156629B (zh) | 图像处理设备和图像处理方法 | |
JP6497912B2 (ja) | 画像処理装置、放射線撮影システム、制御方法、及びプログラム | |
US8805037B2 (en) | Method and system for reconstruction of tomographic images | |
JP2006043431A (ja) | ヘリカルマルチスライスctのための回復ノイズを伴うヘリカルウィンドミルアーチファクトを低減する方法 | |
JP6026214B2 (ja) | 連続マルチスケール再構成において詳細画像を補うx線コンピュータ断層撮像装置(x線ct装置)、医用画像処理装置及び医用画像処理方法 | |
JP2007307417A (ja) | 像データの処理方法及び像データの処理装置 | |
US20070116343A1 (en) | Methods and systems for improving quality of an image | |
EP1649421A2 (en) | Metal artifact correction in computed tomography | |
JP2007319674A (ja) | 対象の画像再構成方法およびその方法を実施するための装置 | |
JP6482934B2 (ja) | 画像処理装置、放射線検出装置および画像処理方法 | |
JP2007050259A (ja) | ボリュームデータ再構成後の断層撮影3d画像のフィルタリング方法 | |
US6751284B1 (en) | Method and system for tomosynthesis image enhancement using transverse filtering | |
JP4266683B2 (ja) | マルチスライス型画像再構成のための方法及び装置 | |
JP2007029735A (ja) | コンピュータ断層撮影システムによるコンピュータ断層画像の作成方法およびコンピュータ断層撮影システム | |
US6775347B2 (en) | Methods and apparatus for reconstructing an image of an object | |
JP6321405B2 (ja) | 画像生成装置、放射線断層撮影装置及び画像生成方法並びにプログラム | |
Maeda et al. | Maximum a posteriori X-ray computed tomography using graph cuts | |
JP2016198504A (ja) | 画像生成装置、x線コンピュータ断層撮影装置及び画像生成方法 | |
JP2006000226A (ja) | X線ct装置 | |
US6999550B2 (en) | Method and apparatus for obtaining data for reconstructing images of an object | |
JP2006167161A (ja) | X線ct画像再構成方法、装置及びx線ct画像再構成プログラムを記憶した記憶媒体 |