JP2012090702A - X線ct装置、x線ct装置の画像再構成方法 - Google Patents
X線ct装置、x線ct装置の画像再構成方法 Download PDFInfo
- Publication number
- JP2012090702A JP2012090702A JP2010239234A JP2010239234A JP2012090702A JP 2012090702 A JP2012090702 A JP 2012090702A JP 2010239234 A JP2010239234 A JP 2010239234A JP 2010239234 A JP2010239234 A JP 2010239234A JP 2012090702 A JP2012090702 A JP 2012090702A
- Authority
- JP
- Japan
- Prior art keywords
- ray
- image
- rotation center
- subject
- projection
- 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
Links
- 238000000034 method Methods 0.000 title claims description 29
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000003786 synthesis reaction Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 230000002194 synthesizing effect Effects 0.000 claims description 6
- 230000001678 irradiating effect Effects 0.000 claims 2
- 238000003384 imaging method Methods 0.000 abstract description 50
- 230000005540 biological transmission Effects 0.000 description 25
- 238000001914 filtration Methods 0.000 description 7
- 239000002131 composite material Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000002872 contrast media Substances 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 241001669679 Eleotris Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Abstract
【解決手段】 複数の撮影の回転中心を設定し、複数の回転中心で撮影された投影データを平行投影像に変換し、複数回転分の前記平行投影像同士の重複領域を計算し、重複領域について設定上の回転中心位置と実際の撮影における回転中心位置に基づき複数回転分の平行投影像を合成する。
【選択図】 図5
Description
まず、図1及び図2に基づいて、本発明が適用されるコーンビームX線CT装置の概略的な構成について説明する。図1は、本発明が適用されるコーンビームX線CT装置(Cアーム方式)1を示す概略構成図である。図2は、本発明が適用される移動型X線装置に搭載されたCアーム方式コーンビームX線CT装置1aを示す概略構成図である。
撮影部10は、被検体2が横臥する寝台17と、被検体2にX線を照射するX線源11と、X線源11に対向して設置され被検体2を透過したX線を検出することによりX線透過像111を出力する2次元X線検出器12と、X線源11及び2次元X線検出器12を機械的に支持するC型アーム13と、C型アーム13を保持するC型アーム保持体14と、C型アーム保持体14をX線撮影室の天井に取り付ける天井支持体15と、天井支持体15を前後左右の2次元方向に移動可能に支持する天井レール16と、被検体2に造影剤を注入するインジェクタ18と、を備える。
制御演算部20は、撮影部10を制御する撮影部制御部100と、撮影部10が出力したX線透過像111を収集して格納する画像収集部110と、収集されたX線透過像111に基づいて3次元CT像を再構成する画像再構成部200と、回転中心位置入力部310、回転幾何学計算部320と、再構成部200が生成した3次元CT像を表示する画像表示部280と、を備える。
撮影部制御部100は、C型アーム13の、回転中心軸4の回りの回転移動を制御する撮影系回転制御部101と、天井支持体15の天井レール16上での位置を制御してC型アーム13の被検体2に対する位置を2次元的に制御する撮影系位置制御部102と、X線管11tに流す管電流のON、OFF等を制御するX線照射制御部103と、インジェクタ18が被検体2に注入する造影剤の注入量及び注入タイミングを制御するインジェクタ制御部104と、寝台17の位置を制御して被検体2の位置を調整するための寝台制御部105と、2次元X線検出器12によるX線透過像111の撮影を制御する検出系制御部107と、を備える。
画像再構成部200は、前処理部210と、フィルタリング部250と、逆投影部260と、平行投影像生成部220と、投影像合成部230と、を備える。
回転中心位置入力部310は、情報入力装置70から、2つの回転撮影の回転中心位置を入力する。ここで、回転中心位置は入力した回転中心位置と実際の回転中心位置との間に数mmの位置誤差を生じる。この位置誤差を補償する平行投影像の位置併せ処理は、後述する図5のステップS231において詳細に説明する。
回転幾何学計算部320は、投影データ211の幾何学系(回転中心(C1x,C1y,C1z))、投影データ212の幾何学系(回転中心(C2x,C2y,C2z))から、次のように定義される合成平行投影像の回転中心(Cx,Cy,Cz)、投影角度β、回転中心からの距離を設定する。
投影角度:β0−90°≦β≦β0+90°
回転中心(Cx,Cy,Cz)からの距離:X-Y面方向 R、Z方向 v
次に、コーンビームX線CT装置1の各種設定条件について次のとおり説明する。
X線管11tから照射されたX線は、被検体2を透過した後、2次元X線検出器12で検出され電気信号に変換される。2次元X線検出器12からの電気信号は、図示しないA/D変換を経た後、デジタル信号からなるX線透過像111として画像収集部110に収集される。2次元X線検出器12は毎秒30フレームで、回転撮影における投影角度間隔は1.33度で、5秒間に150枚のX線透過像111を収集する。200度の回転撮影が完了すると、X線照射制御部103はX線管11tのX線照射を終了し、撮影系回転制御部101は回転減速期間を経たのち回転を停止する。
操作者は、情報入力装置70を用い回転中心位置入力部310に、第1の撮影の回転中心(C1x,C1y,C1z)と、第2の撮影の回転中心(C2x,C2y,C2z)を入力する。回転幾何学計算部320は、回転中心位置入力部310への第1の撮影の回転中心(C1x,C1y,C1z)と、第2の撮影の回転中心(C2x,C2y,C2z)の入力を受けて、第1の撮影の回転中心(C1x,C1y,C1z)と、第2の撮影の回転中心(C2x,C2y,C2z)を記憶する。
回転幾何学計算部320は、第1の回転撮影データと第2の回転撮影データから別々に再構成演算した場合に、2つのX線CT像に重複領域があるかどうかを判定する。
重複領域が存在する(Yes)の場合、回転幾何学計算部320は、合成平行投影像の幾何学(回転中心(Cx=(C1x+C2x)/2,Cy=(C1y+C2y)/2,Cz=(C1z+C2z)/2)と、投影角度β0が{数4}のとおり計算される。ステップS331に進む。
(β0−90°≦β≦β0+90°)
(ステップS331)
制御演算部20は、第1の撮影の回転中心において第1の回転撮影を行う。画像収集部110が収集したX線透過像111は、前処理手段210により前処理が施され、第1の投影データ211を生成する。
制御演算部20は、第2の撮影の回転中心において第2の回転撮影を行う。画像収集部110が収集したX線透過像111は、前処理手段210により前処理が行われ、第2の投影データ212を生成する。
平行投影像変換部220は、合成平行投影像を生成する平行投影角度βを設定する。平行投影角度βは、ステップS320で計算した:β0−90°≦β≦β0+90°の角度範囲をとり、等角度間隔Δβで変化させる。Δβとして、例えば、180°を撮影枚数で割算した値を使用する。
平行投影像変換部220は、第1の撮影処理(ステップS331)が出力した第1の投影データ211から、合成平行投影像の幾何学に変換した、第1の平行投影角度βの平行投影像(221p)を生成する。平行投影変換は、図7に示したサブルーチンより実施される。
合成平行投影像上の画素は(R,v)と定義する。Rは、検出器面左右方向(X-Y面)における合成平行投影像の回転中心(Cx,Cy,Cz)からの距離であり、vは検出器上下方向(Z軸)の距離である。平行投影像変換部220は、合成平行投影像上の画素(R,v)から、ステップS220で設定した平行投影角度βの方向に直線を引く演算を行う。
平行投影像変換部220は、画素は(R,v)と第1の投影データ211の再構成領域1との交点があるかを計算し、交点有り(YES)の場合は、次のステップS227に進む。交点無し(NO)の場合は、ステップS228へ進む。
平行投影像変換部220は、第1の投影データ211(投影角度β1と、検出器座標(u1,v1)で定義されている対応する投影値を2次元X線検出器12から読み出す。
平行投影像変換部220は、現在の平行投影角度βについて、すべての平行投影像の画素(R,v)についてステップS226の交点の有無を計算したか確認し、すべての平行投影像の画素(R,v)について計算していない(NO)の場合は、ステップS225〜ステップS227の処理を繰り返す。すべての平行投影像の画素(R,v)について計算した(YES)の場合は、第1の平行投影角度βの平行投影像(221p)を出力し、ステップS222に進む。
平行投影像変換部220は、第2の撮影処理(ステップS332)が出力した第2の投影データ212から、第2の平行投影角度βの平行投影像(222p)を生成する。
平行投影像変換部220は、Δu方向(横方向)に検出器ピッチの数分の1(0.1mm程度)、Δv方向(縦方向)には検出器ピッチ(0.5mm)程度の精度で、位置合わせをする。
(ステップS233)
投影像合成部230は、ステップS231が計算した補正量(ΔR,Δv)を用いて、第1の平行投影像(221p)と第2の平行投影像(222p)を合成し、平行投影角度βの合成投影像233を生成する。
フィルタリング手段250は、平行投影角度βの合成投影像233に対するフィルタリング処理を行う。
逆投影手段260は、フィルタリング処理された平行投影角度βの合成投影像の逆投影処理を行う。
画像再構成部200は、すべての平行投影角度βについて処理したかを確認し、すべての平行投影角度βについて処理していない(NO)の場合は、平行投影角度βをΔβ増分し、ステップS220〜ステップS260の処理を繰り返す。すべての平行投影角度βについて処理した(YES)の場合は、ステップS280に進む。
画像再構成部200は、すべての平行投影角度βについて逆投影した再構成像を、いわゆるHounsfield値として定義されたCT値に変換し、X線CT像を出力する。
平行投影像変換部220は、ステップS231の{数4}において相関値を計算する際に計算したROI平均値を使用する。
ROI平均値は次のようにそれぞれ示される。
各ROI平均値中間値を求め、第1の平行投影像(221p)を各ROI平均値の加算値を前者のROI平均値の2倍で除した数を乗じ、第2の平行投影像(222p)を各ROI平均値の加算値を後者のROI平均値の2倍で除した数を乗じて、2つの投影データの輝度を等しくしてから、ステップS233に進む。
平行投影像変換部220は、回転中心軸変位・検出器傾き角変位テーブルを参照し、第1の撮影の回転中心軸方向ベクトル(e1x,e1y,e1z)、検出器傾き角θ1と、第2の撮影の回転中心軸方向ベクトル(e2x,e2y,e2z)、検出器傾き角θ2を計算する。
平行投影像変換部220は、合成平行投影像上の画素(R,v)から、ステップS220で設定した平行投影角度βの方向に直線を引く演算を行う。
平行投影像変換部220は、画素(R,v)が第1の投影データ211の再構成領域と交点があるかを計算し、画素(R,v)が第1の投影データ211の再構成領域と交点がある(YES)の場合は、次のステップS227に進む。画素(R,v)が第1の投影データ211の再構成領域と交点がない(NO)の場合は、ステップS228へ進む。
平行投影像変換部220は、回転中心軸方向ベクトル(e1x,e1y,e1z)と、平行投影角度βの検出器面、回転軌道面方向の単位ベクトル(cosβ,−sinβ,0)との内積tを計算する。第1の投影データ211を、−θ1−arctan(t)だけ回転し、合成平行投影像の画素(R,v)に対応する投影値を読み出す。
平行投影像変換部220は、現在の平行投影角度βについて、すべての(R,v)についてスキャンしたかを確認し、NOの場合は、ステップS225〜ステップS227の処理を繰り返す。YESの場合は、第1の平行投影角度βの平行投影像(221p)を出力し、ステップS222に進む。
実施例1〜3で説明した構成は一例であり、本発明は、技術思想を逸脱しない範囲内で適宜変更が可能である。
Claims (6)
- 被検体にX線を照射するX線源と、
前記X線源に対向配置され、前記被検体を透過したX線を検出して前記被検体の投影データを出力するX線検出器と、
前記X線源と前記X線検出器を前記被検体の回転中心で回転させる回転部と、
前記被検体の投影データに基づいて前記被検体のX線CT像を再構成する画像再構成演算部と、
前記X線CT像を表示する表示部と、を備えたX線CT装置であって、
複数の撮影の回転中心を設定する回転中心設定部を備え、
前記画像再構成演算部は、
複数の回転中心で撮影された投影データを平行投影像に変換する平行投影像変換部と、
複数回転分の前記平行投影像同士の重複領域を計算し、前記重複領域について設定上の回転中心位置と実際の撮影における回転中心位置に基づき複数回転分の平行投影像を合成する投影像合成部と、
を具備したことを特徴とするX線CT装置。 - 前記投影像合成部は、前記設定上の回転中心位置と実際の撮影における回転中心位置の誤差値を計算し、該誤差値に基づいて設定上の回転中心位置と実際の回転中心位置の誤差を補正し、前記複数回転分の平行投影像を合成することを特徴とする請求項1に記載のX線CT装置。
- 前記投影像合成部は、さらに、前記重複領域の投影データの出力レベルを合わせて前記複数回転分の平行投影像を合成することを特徴とする請求項1に記載のX線CT装置。
- 前記投影像合成部は、前記X線検出器の傾き変位に基づく誤差を補正することを特徴とする請求項1に記載のX線CT装置。
- 前記投影像合成部は、複数回転の回転中心軸の傾き変位に基づく誤差を補正することを特徴とする請求項1に記載のX線CT装置。
- 被検体にX線を照射するX線源と、
前記X線源に対向配置され、前記被検体を透過したX線を検出して前記被検体の投影データを出力するX線検出器と、
前記X線源と前記X線検出器を前記被検体の回転中心で回転させる回転部と、
前記被検体の投影データに基づいて前記被検体のX線CT像を再構成する画像再構成演算部と、
前記X線CT像を表示する表示部と、を備えたX線CT装置の画像再構成方法であって、
複数の撮影の回転中心を設定するステップと、
複数の回転中心で撮影された投影データを平行投影像に変換するステップと、
複数回転分の前記平行投影像同士の重複領域を計算し、前記重複領域について設定上の回転中心位置と実際の撮影における回転中心位置に基づき複数回転分の平行投影像を合成するステップと、
を含むことを特徴とするX線CT装置の画像再構成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010239234A JP5572521B2 (ja) | 2010-10-26 | 2010-10-26 | X線ct装置、x線ct装置の画像再構成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010239234A JP5572521B2 (ja) | 2010-10-26 | 2010-10-26 | X線ct装置、x線ct装置の画像再構成方法 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2012090702A true JP2012090702A (ja) | 2012-05-17 |
JP2012090702A5 JP2012090702A5 (ja) | 2013-12-05 |
JP5572521B2 JP5572521B2 (ja) | 2014-08-13 |
Family
ID=46384794
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010239234A Expired - Fee Related JP5572521B2 (ja) | 2010-10-26 | 2010-10-26 | X線ct装置、x線ct装置の画像再構成方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5572521B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014533168A (ja) * | 2011-11-11 | 2014-12-11 | コーニンクレッカ フィリップス エヌ ヴェ | 拡張された視野を備えるcアーム型システム |
CN112461165A (zh) * | 2019-09-09 | 2021-03-09 | 株式会社三丰 | X射线测量装置的校正方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61181443A (ja) * | 1985-01-03 | 1986-08-14 | ゼネラル・エレクトリツク・カンパニイ | 回転式扇形ビ−ムct装置の回転中心のずれを補正する方法と装置 |
JPH10243941A (ja) * | 1997-03-05 | 1998-09-14 | Toshiba Corp | 画像再構成処理装置 |
JPH11188030A (ja) * | 1997-12-25 | 1999-07-13 | Hitachi Medical Corp | X線ct装置 |
JP2000237182A (ja) * | 1999-02-17 | 2000-09-05 | Toshiba Corp | X線ct装置 |
JP2004065706A (ja) * | 2002-08-08 | 2004-03-04 | Ge Medical Systems Global Technology Co Llc | 投影データ補正方法、画像生成方法およびx線ct装置 |
JP2005144167A (ja) * | 2003-11-11 | 2005-06-09 | Ge Medical Systems Global Technology Co Llc | コンピュータ断層撮影イメージングシステムにおけるアーチファクトを低減するための方法及び装置 |
JP2008055171A (ja) * | 2006-09-01 | 2008-03-13 | Siemens Ag | 3次元画像ボリュームの再構成方法およびx線装置 |
JP2009082250A (ja) * | 2007-09-28 | 2009-04-23 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
-
2010
- 2010-10-26 JP JP2010239234A patent/JP5572521B2/ja not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61181443A (ja) * | 1985-01-03 | 1986-08-14 | ゼネラル・エレクトリツク・カンパニイ | 回転式扇形ビ−ムct装置の回転中心のずれを補正する方法と装置 |
JPH10243941A (ja) * | 1997-03-05 | 1998-09-14 | Toshiba Corp | 画像再構成処理装置 |
JPH11188030A (ja) * | 1997-12-25 | 1999-07-13 | Hitachi Medical Corp | X線ct装置 |
JP2000237182A (ja) * | 1999-02-17 | 2000-09-05 | Toshiba Corp | X線ct装置 |
JP2004065706A (ja) * | 2002-08-08 | 2004-03-04 | Ge Medical Systems Global Technology Co Llc | 投影データ補正方法、画像生成方法およびx線ct装置 |
JP2005144167A (ja) * | 2003-11-11 | 2005-06-09 | Ge Medical Systems Global Technology Co Llc | コンピュータ断層撮影イメージングシステムにおけるアーチファクトを低減するための方法及び装置 |
JP2008055171A (ja) * | 2006-09-01 | 2008-03-13 | Siemens Ag | 3次元画像ボリュームの再構成方法およびx線装置 |
JP2009082250A (ja) * | 2007-09-28 | 2009-04-23 | Ge Medical Systems Global Technology Co Llc | X線ct装置 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014533168A (ja) * | 2011-11-11 | 2014-12-11 | コーニンクレッカ フィリップス エヌ ヴェ | 拡張された視野を備えるcアーム型システム |
CN112461165A (zh) * | 2019-09-09 | 2021-03-09 | 株式会社三丰 | X射线测量装置的校正方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5572521B2 (ja) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6580777B1 (en) | X-ray CT apparatus | |
RU2394494C2 (ru) | Система сканирования х-кт | |
JP3378401B2 (ja) | X線装置 | |
JP5331940B2 (ja) | 放射線撮影システム及び放射線画像生成方法 | |
JP5384521B2 (ja) | 放射線撮像装置 | |
JP3548339B2 (ja) | X線撮影装置 | |
US20080008372A1 (en) | A method and system for reducing artifacts in a tomosynthesis imaging system | |
JP5375655B2 (ja) | 放射線撮影装置 | |
JP5878119B2 (ja) | X線ct装置、及びその制御方法 | |
KR101158619B1 (ko) | 디지털 엑스-레이 장치(변형들)를 교정하기 위한 방법 | |
US9532755B2 (en) | X-ray computed tomographic imaging apparatus and method for same | |
JP3897925B2 (ja) | コーンビームct装置 | |
WO2006028085A1 (ja) | X線ct装置、画像処理プログラム、及び画像処理方法 | |
JP3540916B2 (ja) | 3次元x線ct装置 | |
JP4408664B2 (ja) | コーンビームx線ct装置及びそれに用いるファントム | |
JP3583554B2 (ja) | コーンビームx線断層撮影装置 | |
JP5060862B2 (ja) | 断層撮影装置 | |
JP4573593B2 (ja) | X線画像補正方法及び装置 | |
JP5572521B2 (ja) | X線ct装置、x線ct装置の画像再構成方法 | |
JP2008012129A (ja) | X線ct装置 | |
JP5808672B2 (ja) | X線ct装置及びx線ct撮影方法 | |
JP2009279301A (ja) | コーンビームx線ct装置 | |
JP3583567B2 (ja) | X線断層撮像装置 | |
JP2008036272A (ja) | コーンビームx線ct装置 | |
JP3484283B2 (ja) | X線撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131018 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20131018 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140328 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140507 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140526 |
|
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: 20140616 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140630 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5572521 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140808 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |