JP2006167161A5

JP2006167161A5 -

Info

Publication number: JP2006167161A5
Application number: JP2004363904A
Authority: JP
Filing date: 2004-12-16
Publication date: 2007-10-25

Claims

撮影対象物の周囲を多数の方向からＸ線により投影して投影データ（Ｉ_１）を取得する投影データ取得工程と、
該投影データに対してアーチファクトを低減する処理を行って該投影データの修正投影データ（Ｉ_１’）を演算する修正投影データ演算工程と、
該修正投影データを逆投影演算して画像を再構成する画像再構成工程と
を具備するＸ線ＣＴ画像再構成方法。
前記修正投影データ演算工程は、
前記投影データ（Ｉ_１）より高Ｘ線吸収係数の物体領域を抽出する高Ｘ線吸収係数物体領域抽出工程と、
前記投影データより前記高Ｘ線吸収係数の物体領域を除去した第１の投影データ（Ｉ_１Ａ）を演算する第１の投影データ演算工程と、
該第１の投影データに対して連続的に増大変化する前記高Ｘ線吸収係数の物体領域の第２の投影データ（Ｉ_１Ｂ）を演算する第２の投影データ演算工程と、
前記第１の投影データと前記第２の投影データとの合成により前記修正投影データ（Ｉ_１’）を演算する工程と
を具備する請求項１に記載のＸ線ＣＴ画像再構成方法。
前記高Ｘ線吸収係数物体領域抽出工程は、
前記第１の投影データの各Ｘ線吸収係数と第１のしきい値（ＴＨ１）との比較により前記投影データの高Ｘ線吸収係数領域を抽出する工程と、
前記投影データの差分投影データ（ΔＩ_１）の各値（Δｆ（ｉ，ｊ））と第２のしきい値（ＴＨ２）との比較により前記投影データの不連続領域を抽出する不連続領域抽出工程と、
前記高Ｘ線吸収係数領域と前記不連続領域との合成により前記高Ｘ線吸収係数の物体領域を演算する高Ｘ線吸収係数物体領域演算工程と
を具備する請求項２に記載のＸ線ＣＴ画像再構成方法。
前記第１の投影データ演算工程は、
前記投影データの各Ｘ線吸収係数のうち前記高Ｘ線吸収係数の物体領域の各Ｘ線吸収係数を０とする工程
を具備する請求項２に記載のＸ線ＣＴ画像再構成方法。
第２の投影データ演算工程は、
前記第１の投影データ（Ｉ_１Ａ）を逆投影して再構成画像（Ｐ_０）を演算する工程と、
該再構成画像より前記高Ｘ線吸収係数の物体領域画像を抽出する工程と、
該物体領域画像を演算により投影して前記第２の投影データを取得する工程と
を具備する請求項２に記載のＸ線ＣＴ画像再構成方法。
前記物体領域画像の演算による投影は、該物体領域画像のＸ線吸収係数を任意の大きさとし、前記物体領域画像以外の画像のＸ線吸収係数を０として行われる請求項５に記載のＸ線ＣＴ画像再構成方法。
前記投影データ取得工程は前記Ｘ線の異なる投影エネルギーを用いて低Ｘ線投影データ（Ｉ_１Ａ）及び高Ｘ線投影データ（Ｉ_１Ｂ）を取得し、
前記修正投影データ演算工程は前記低Ｘ線投影データと前記高Ｘ線投影データとの合成により前記修正投影データを演算する請求項１に記載のＸ線ＣＴ画像再構成方法。
前記投影データ取得工程は前記Ｘ線の低投影エネルギーを用いて低Ｘ線投影データ（Ｉ_１Ａ）を取得すると共に、既に抽出されている高Ｘ線吸収係数の物体領域を演算により投影して高Ｘ線投影データ（Ｉ_１Ｂ）を取得し、
前記修正投影データ演算工程は前記低Ｘ線投影データと前記高Ｘ線投影データとの合成により前記修正投影データを演算する請求項１に記載のＸ線ＣＴ画像再構成方法。
前記修正投影データ演算工程は、
前記高Ｘ線投影データの各Ｘ線吸収係数と第３のしきい値（ＴＨ３）との比較により高Ｘ線吸収係数の物体領域を判別する工程と、
該高Ｘ線吸収係数の物体領域外の前記低Ｘ線投影データと前記高Ｘ線吸収係数の物体領域内の高Ｘ線投影データとの合成により前記修正投影データを演算すると共に、前記低Ｘ線投影データと前記高Ｘ線投影データとの接続点（Ｐ，Ｑ）を連続となるように合成する工程と
を具備する請求項８あるいは９に記載のＸ線ＣＴ画像再構成方法。
さらに、
前記再構成された画像（Ｐ_ｋ）を前回値画像（Ｐ_ｋ−１）として演算により投影して前回値投影データ（Ｉ_ｋ−１）を取得する前回値投影データ取得工程と、
前記修正投影データ（Ｉ_１’）と前記前回値投影データとの差投影データ（ΔＩ）を演算する差投影データ演算工程と、
該差投影データを逆投影して差画像（ΔＰ）を演算する差画像演算工程と、
該差画像を前記前回値画像に加算することにより新たな再構成画像（Ｐ_ｋ）を取得する再構成画像取得工程と、
前記前回値投影データ取得工程、前記差投影データ演算工程、前記差画像演算工程及び前記再構成画像取得工程を前記差投影データが収束するまで繰返す繰返工程と
を具備する請求項１に記載のＸ線ＣＴ画像再構成方法。
さらに、前記差投影データ演算前に、前回値投影データのスケールを前記修正投影データのスケールに調整する工程を具備する請求項１０に記載のＸ線ＣＴ画像再構成方法。
前記前回値投影データのスケールの前記修正投影データのスケールへの調整は、前回値投影データの最大値及び最小値が前記修正投影データの最大値及び最小値にそれぞれ一致するように前記前回値投影データを１次変換することにより行う請求項１１に記載のＸ線ＣＴ画像再構成方法。
撮影対象物の周囲を多数の方向からＸ線により投影して投影データ（Ｉ_１）を取得し、該投影データを逆投影演算して画像を再構成するＸ線ＣＴ画像再構成方法において、
前記再構成された画像（Ｐ_ｋ）を前回値画像（Ｐ_ｋ−１）として演算により投影して前回値投影データ（Ｉ_ｋ−１）を取得する前回値投影データ取得工程と、
前回値投影データのスケールを前記修正投影データのスケールに調整するスケール調整工程と、
前記修正投影データ（Ｉ_１’）と前記前回値画像との差投影データ（ΔＩ）を演算する差投影データ演算工程と、
該差投影データを逆投影して差画像（ΔＰ）を演算する差画像演算工程と、
該差画像を前記前回値画像に加算することにより新たな再構成画像（Ｐ_ｋ）を取得する再構成画像取得工程と、
前記前回値投影データ取得工程、前記差投影データ演算工程、前記差画像演算工程及び前記再構成画像取得工程を前記差投影データが収束するまで繰返す繰返工程と
を具備することを特徴とするＸ線ＣＴ画像再構成方法。
前記前回値投影データのスケールの前記修正投影データのスケールへの調整は、前回値投影データの最大値及び最小値が前記修正投影データの最大値及び最小値にそれぞれ一致するように前回値投影データを１次変換することにより行う請求項１３に記載のＸ線ＣＴ画像再構成方法。
撮影対象物の周囲を多数の方向からＸ線により投影して投影データ（Ｉ_１）を取得する投影データ取得手順と、
該投影データに対してアーチファクトを低減する処理を行って該投影データの修正投影データ（Ｉ_１’）を演算する修正投影データ演算手順と、
該修正投影データを逆投影演算して画像を再構成する画像再構成手順と
を実行させるためのＸ線ＣＴ画像再構成プログラム。
前記修正投影データ演算手順は、
前記投影データ（Ｉ_１）より高Ｘ線吸収係数の物体領域を抽出する高Ｘ線吸収係数物体領域抽出手順と、
前記投影データより前記高Ｘ線吸収係数の物体領域を除去した第１の投影データ（Ｉ_１Ａ）を演算する第１の投影データ演算手順と、
該第１の投影データに対して連続的に増大変化する前記高Ｘ線吸収係数の物体領域の第２の投影データ（Ｉ_１Ｂ）を演算する第２の投影データ演算手順と、
前記第１の投影データと前記第２の投影データとの合成により前記修正投影データ（Ｉ_１’）を演算する手順と
を具備する請求項１５に記載のＸ線ＣＴ画像再構成プログラム。
前記高Ｘ線吸収係数物体領域抽出手順は、
前記第１の投影データの各Ｘ線吸収係数と第１のしきい値（ＴＨ１）との比較により前記投影データの高Ｘ線吸収係数領域を抽出する手順と、
前記投影データの差分投影データ（ΔＩ_１）の各値（Δｆ（ｉ，ｊ））と第２のしきい値（ＴＨ２）との比較により前記投影データの不連続領域を抽出する不連続領域抽出手順と、
前記高Ｘ線吸収係数領域と前記不連続領域との合成により前記高Ｘ線吸収係数の物体領域を演算する高Ｘ線吸収係数物体領域演算手順と
を具備する請求項１６に記載のＸ線ＣＴ画像再構成プログラム。
前記第１の投影データ演算手順は、
前記投影データの各Ｘ線吸収係数のうち前記高Ｘ線吸収係数の物体領域の各Ｘ線吸収係数を０とする手順
を具備する請求項１６に記載のＸ線ＣＴ画像再構成プログラム。
第２の投影データ演算手順は、
前記第１の投影データ（Ｉ_１Ａ）を逆投影して再構成画像（Ｐ_０）を演算する手順と、
該再構成画像より前記高Ｘ線吸収係数の物体領域画像を抽出する手順と、
該物体領域画像を演算により投影して前記第２の投影データを取得する手順と
を具備する請求項１６に記載のＸ線ＣＴ画像再構成プログラム。
前記物体領域画像の演算による投影は、該物体領域画像のＸ線吸収係数を任意の大きさとし、前記物体領域画像以外の画像のＸ線吸収係数を０として行われる請求項１９に記載のＸ線ＣＴ画像再構成プログラム。
前記投影データ取得手順は前記Ｘ線の異なる投影エネルギーを用いて低Ｘ線投影データ（Ｉ_１Ａ）及び高Ｘ線投影データ（Ｉ_１Ｂ）を取得し、
前記修正投影データ演算手順は前記低Ｘ線投影データと前記高Ｘ線投影データとの合成により前記修正投影データを演算する請求項１５に記載のＸ線ＣＴ画像再構成プログラム。
前記投影データ取得手順は前記Ｘ線の低投影エネルギーを用いて低Ｘ線投影データ（Ｉ_１Ａ）を取得すると共に、既に抽出されている高Ｘ線吸収係数の物体領域を演算により投影して高Ｘ線投影データ（Ｉ_１Ｂ）を取得し、
前記修正投影データ演算手順は前記低Ｘ線投影データと前記高Ｘ線投影データとの合成により前記修正投影データを演算する請求項１５に記載のＸ線ＣＴ画像再構成プログラム。
前記修正投影データ演算手順は、
前記高Ｘ線投影データの各Ｘ線吸収係数と第３のしきい値（ＴＨ３）との比較により高Ｘ線吸収係数の物体領域を判別する手順と、
該高Ｘ線吸収係数の物体領域外の前記低Ｘ線投影データと前記高Ｘ線吸収係数の物体領域内の高Ｘ線投影データとの合成により前記修正投影データを演算すると共に、前記低Ｘ線投影データと前記高Ｘ線投影データとの接続点（Ｐ，Ｑ）を連続となるように合成する手順と
を具備する請求項２１あるいは２２に記載のＸ線ＣＴ画像再構成プログラム。
さらに、
前記再構成された画像（Ｐ_ｋ）を前回値画像（Ｐ_ｋ−１）として演算により投影して前回値投影データ（Ｉ_ｋ−１）を取得する前回値投影データ取得手順と、
前記修正投影データ（Ｉ_１’）と前記前回値投影データとの差投影データ（ΔＩ）を演算する差投影データ演算手順と、
該差投影データを逆投影して差画像（ΔＰ）を演算する差画像演算手順と、
該差画像を前記前回値画像に加算することにより新たな再構成画像（Ｐ_ｋ）を取得する再構成画像取得手順と、
前記前回値投影データ取得手順、前記差投影データ演算手順、前記差画像演算手順及び前記再構成画像取得手順を前記差投影データが収束するまで繰返す繰返手順と
を具備する請求項１５に記載のＸ線ＣＴ画像再構成プログラム。
さらに、前記差投影データ演算前に、前回値投影データのスケールを前記修正投影データのスケールに調整する手順を具備する請求項２４に記載のＸ線ＣＴ画像再構成プログラム。
前記前回値投影データのスケールの前記修正投影データのスケールへの調整は、前回値投影データの最大値及び最小値が前記修正投影データの最大値及び最小値にそれぞれ一致するように前回値投影データを１次変換することにより行う請求項２５に記載のＸ線ＣＴ画像再構成プログラム。
撮影対象物の周囲を多数の方向からＸ線により投影して投影データ（Ｉ_１）を取得し、該投影データを逆投影演算して画像を再構成するためのＸ線ＣＴ画像再構成プログラムにおいて、
前記再構成された画像（Ｐ_ｋ）を前回値画像（Ｐ_ｋ−１）として演算により投影して前回値投影データ（Ｉ_ｋ−１）を取得する前回値投影データ取得手順と、
前回値投影データのスケールを前記修正投影データのスケールに調整するスケール調整手順と、
前記修正投影データ（Ｉ_１’）と前記前回値画像との差投影データ（ΔＩ）を演算する差投影データ演算手順と、
該差投影データを逆投影して差画像（ΔＰ）を演算する差画像演算手順と、
該差画像を前記前回値画像に加算することにより新たな再構成画像（Ｐ_ｋ）を取得する再構成画像取得手順と、
前記前回値投影データ取得手順、前記差投影データ演算手順、前記差画像演算手順及び前記再構成画像取得手順を前記差投影データが収束するまで繰返す繰返手順と
を実行させるためのＸ線ＣＴ画像再構成プログラム。
前記前回値投影データのスケールの前記修正投影データのスケールへの調整は、前回値投影データの最大値及び最小値が前記修正投影データの最大値及び最小値にそれぞれ一致するように前記前回値投影データを１次変換することにより行う請求項２７に記載のＸ線ＣＴ画像再構成プログラム。