JP4350206B2 - Ｘ線ｃｔ装置 - Google Patents
Ｘ線ｃｔ装置 Download PDFInfo
- Publication number
- JP4350206B2 JP4350206B2 JP13920099A JP13920099A JP4350206B2 JP 4350206 B2 JP4350206 B2 JP 4350206B2 JP 13920099 A JP13920099 A JP 13920099A JP 13920099 A JP13920099 A JP 13920099A JP 4350206 B2 JP4350206 B2 JP 4350206B2
- Authority
- JP
- Japan
- Prior art keywords
- contrast
- image
- axial
- scout image
- scout
- 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.)
- Expired - Fee Related
Links
- 239000002872 contrast media Substances 0.000 claims description 30
- 238000000034 method Methods 0.000 description 36
- 210000004204 blood vessel Anatomy 0.000 description 15
- 238000010586 diagram Methods 0.000 description 6
- 239000008186 active pharmaceutical agent Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、閾値決定方法および造影領域３次元データ作成方法に関し、さらに詳しくは、造影領域を抽出する基となる複数の造影アキシャル画像の体軸方向の長さが大きいときでも適切な閾値を容易に決定できる閾値決定方法および当該閾値を用いて造影領域３次元データを作成する造影領域３次元データ作成方法に関する。
【０００２】
【従来の技術】
従来のＸ線ＣＴ装置で、被検体内の血管を３次元表示する手順について説明する。
（１）図１０に示すように、被検体Ｈの造影対象の血管αに造影剤Ｍを注入した後、息止めした被検体Ｈの体軸方向のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応したアキシャル画像（造影アキシャル画像）Ａ１，Ａ２，…，Ａ６を連続して撮影する。
（２）造影アキシャル画像Ａ１，Ａ２，…，Ａ６における血管αのＣＴ値（CT１'，CT2'，CT3u'，CT3d'，…，CT6u'，CT6d'）に基づいて、造影領域を抽出するための１つの閾値Ｔｈを決定する。
（３）造影アキシャル画像Ａ１，Ａ２，…，Ａ６から、前記閾値Ｔｈに基づいて造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。
（４）図１１に示すように、作成した３次元データに基づいて、画面上に、血管αに相当する造影領域３次元画像Ｉを表示する。
【０００３】
【発明が解決しようとする課題】
旧式のＸ線ＣＴ装置では、スキャン速度が遅いので、１回の息止めで撮影できる体軸方向の長さＤ’（図１０参照）はせいぜい２０ｃｍ程度であった。この程度の長さＤ’では、血管αの全域で造影剤Ｍの濃度をほぼ均一と見なすことができ、１つの閾値Ｔｈでも造影領域の抽出を支障なく行えた。
【０００４】
これに対し、最近のＸ線ＣＴ装置では、多層型検出器の採用等によりスキャン速度が高速化されているので、１回の息止めで撮影できる体軸方向の長さＤ’が、例えば４０ｃｍ以上に増大してきている。ところが、このように長さＤ’が増大すると、血管αの区間によって造影剤Ｍの濃度が著しく異なる（例えば最上流側と最下流端とで異なる）ため、１つの閾値Ｔｈでは、本来の造影領域を正しく抽出できない問題点がある。すなわち、血管αの最小ＣＴ値に合わせて閾値Ｔｈを設定すると、閾値Ｔｈが過小となって、臓器などの非造影対象が誤って抽出されてしまう。また、血管αの最大ＣＴ値に合わせて閾値Ｔｈを設定すると、ＣＴ値の小さい血管αの区間を正しく抽出できなくなってしまう。
【０００５】
なお、上記問題点を解消するため、操作者が、造影アキシャル画像Ａ１，Ａ２，…，Ａ６ごとに血管αのＣＴ値を把握して、造影アキシャル画像ごとに閾値を個別に手動設定することも考えられるが、手間と時間がかかる問題点がある。
【０００６】
そこで、本発明の目的は、造影領域を抽出する基となる複数の造影アキシャル画像の体軸方向の長さが大きいときでも適切な閾値を容易に決定できる閾値決定方法および当該閾値を用いて造影領域３次元データを作成する造影領域３次元データ作成方法を提供することにある。
【０００７】
【課題を解決するための手段】
第１の観点では、本発明は、造影剤を注入しないで被検体の体軸方向のスカウト画像（プレーンスカウト画像）を撮影するプレーンスカウト画像撮影過程と、造影剤を注入して被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影過程と、前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成過程と、前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成過程と、前記差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定過程とを有することを特徴とする閾値決定方法を提供する。
上記第１の観点による閾値決定方法では、造影剤を注入しないで被検体の体軸方向のスカウト（scout）画像を撮影することで、プレーンスカウト（plain scout）画像を得る。また、造影剤を注入して各スライス位置のアキシャル（axial）画像を撮影して造影アキシャル画像を得て、それらの造影アキシャル画像を基にしてプレーンスカウト画像に相当する造影スカウト画像を作成する。
プレーンスカウト画像の画素値と，造影スカウト画像の画素値は、造影剤の存在する領域すなわち造影領域のみで実質的に異なり、造影領域以外の領域で同等と見なせる。したがって、差分画像として、スカウト画像内の造影領域のみが描出される。
この結果、差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定することで、造影領域の区間によって造影剤の濃度が著しく異なる場合でも、適切な閾値をきめ細かく決定することが出来る。
【０００８】
第２の観点では、本発明は、造影剤を注入しないで被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（プレーンアキシャル画像）を撮影するプレーンアキシャル画像撮影過程と、造影剤を注入して前記被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影過程と、前記プレーンアキシャル画像を基にして演算により被検体の体軸方向のスカウト画像（プレーンスカウト画像）を作成するプレーンスカウト画像作成過程と、前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成過程と、前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成過程と、前記差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定過程とを有することを特徴とする閾値決定方法を提供する。
上記第２の観点による閾値決定方法では、造影剤を注入しないで被検体の体軸方向の各スライス位置でアキシャル画像を撮影することでプレーンスアキシャル画像を得て、それらのプレーンスアキシャル画像を基にしてプレーンスカウト画像を得る。また、造影剤を注入して各スライス位置のアキシャル画像を撮影して造影アキシャル画像を得て、それらの造影アキシャル画像を基にしてプレーンスカウト画像に相当する造影スカウト画像を作成する。
プレーンスカウト画像の画素値と，造影スカウト画像の画素値は、造影剤の存在する領域すなわち造影領域のみで実質的に異なり、造影領域以外の領域で同等と見なせる。したがって、差分画像として、スカウト画像内の造影領域のみが描出される。なお、プレーンスカウト画像と，造影スカウト画像のいずれもが、画質の高いアキシャル画像を基にして作成されているので、造影領域を正確に抽出することが出来る。
この結果、差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定することで、造影領域の区間によって造影剤の濃度が著しく異なる場合でも、適切な閾値をきめ細かく決定することが出来る。
【０００９】
第３の観点では、本発明は、上記構成の閾値決定方法により決定した閾値を用いて前記造影アキシャル画像のそれぞれから造影領域を抽出する造影領域抽出過程と、前記抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する３次元データ作成過程とを有することを特徴とする造影領域３次元データ作成方法を提供する。
上記第３の観点による造影領域３次元データ作成方法では、上記閾値決定方法により決定した閾値を用いて、造影アキシャル画像のそれぞれから造影領域を過不足なく抽出することが可能となり、精度の高い造影領域３次元データを作成することが出来る。
【００１０】
【発明の実施の形態】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
−第１の実施形態−
図１は、本発明の第１の実施形態にかかる造影領域３次元データ作成方法を実施する高速スキャン型のＸ線ＣＴ装置１００を示すブロック図である。
このＸ線ＣＴ装置１００は、操作コンソール１と、撮影テーブル１０と、走査ガントリ２０とを具備している。
前記操作コンソール１は、操作者の指示や情報などを受け付ける入力装置２と、スキャン処理や補間処理や画像再構成処理やレイサム（ray-summation）演算処理などを実行する中央処理装置３と、制御信号などを撮影テーブル１０や走査ガントリ２０へ出力する制御インタフェース４と、走査ガントリ２０で取得したデータを収集するデータ収集バッファ５と、画像などを表示するＣＲＴ６と、各種のデータやプログラムを記憶する記憶装置７とを具備している。
前記撮影テーブル１０は、被検体を乗せて体軸方向に移動させる。
前記走査ガントリ２０は、ファンビーム方式のＸ線管２１，コリメータ２２および検出器２３と、被検体の体軸の回りにＸ線管２１や検出器２３などを回転させる回転コントローラ２４と、Ｘ線照射のタイミングや強度を調整するＸ線コントローラ２５と、データ収集部２６とを具備している。
【００１１】
図２は、このＸ線ＣＴ装置１００による造影領域３次元データ作成処理を示すフロー図である。
ステップＳＴ１では、造影剤を注入しないで被検体の体軸方向のスカウト画像（プレーンスカウト画像）を撮影する。一般に、前記スカウト画像は、コロナル（coronal）画像またはサジタル（sagital）画像である。すなわち、図３に示すように、Ｘ線管２１および検出器２３を回転させずに被検体Ｈの体軸方向に長さＤだけ直線的に相対移動させ、被検体Ｈを透過したファンビームＢのＸ線強度を検出器２３の各チャネルで検出し、当該検出強度に対応した画素値をスカウト画像平面Ｆ上に投影した如きスカウト画像を作成する。前記長さＤは、造影領域の区間に相当する長さであり、例えば４０ｃｍ程度である。図４に、プレーンスカウト画像ＰＳを例示する。造影剤を注入しない状態では、プレーンスカウト画像ＰＳ内に血管αは現れない。
【００１２】
ステップＳＴ２では、図５に示すように、被検体Ｈの血管αに造影剤Ｍを注入した後、息止めした被検体Ｈの体軸方向のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応したアキシャル画像（造影アキシャル画像）Ａ１，Ａ２，…，Ａ６を連続して撮影する。図示の都合上、スライス数を“６”としたが、スライス数は任意である。
ステップＳＴ３では、造影アキシャル画像Ａ１，Ａ２，…，Ａ６を基にして、演算により前記プレーンスカウト画像ＰＳに相当する造影スカウト画像を作成する。すなわち、図６に示すように、造影アキシャル画像Ａ１，Ａ２，…，Ａ６の各スライス面内で、ファンビームＢに相当する投影線Ｂ’を前記スカウト画像平面Ｆ上に投影したときの各投影線上にある画素に対応するＣＴ値を積算するレイサム演算処理を行い、その結果を新たな画素値とした如きレイサム画像を造影スカウト画像ＲＳとする。図７に、前記造影スカウト画像ＲＳを例示する。血管αに相当する造影領域α’の画素値は、造影剤Ｍの濃度を反映して変化する。
【００１３】
ステップＳＴ４では、プレーンスカウト画像ＰＳと造影スカウト画像ＲＳの差分をとって、差分画像を作成する。図８に、差分画像ＤＳを例示する。プレーンスカウト画像ＰＳの画素値と，造影スカウト画像ＲＳの画素値は、造影領域α’のみで実質的に異なり、造影領域α’以外の領域で同等と見なせる。したがって、臓器や骨格などの領域の画素値が相殺され、造影領域α’のみが描出される。ステップＳＴ５では、差分画像ＤＳ中のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応する位置での造影領域α’のＣＴ値に基づいて造影領域を抽出するための閾値を決定する。すなわち、図８に示すように、スライス位置Ｌ１に対応する位置でのＣＴ値CT1に基づいて閾値Ｔｈ１を決定し、スライス位置Ｌ２に対応するＣＴ値CT2に基づいて閾値Ｔｈ２を決定し、スライス位置Ｌ３に対応するＣＴ値CT3u，CT3dに基づいて閾値Ｔｈ３を決定し、スライス位置Ｌ４に対応するＣＴ値CT4u，CT4dに基づいて閾値Ｔｈ４を決定し、スライス位置Ｌ５に対応するＣＴ値CT5u，CT5dに基づいて閾値Ｔｈ５を決定し、スライス位置Ｌ６に対応するＣＴ値CT6u，CT6dに基づいて閾値Ｔｈ６を決定する。
【００１４】
ステップＳＴ６では、上記ステップＳＴ５で決定した閾値Ｔｈ１，Ｔｈ２，…，Ｔｈ６を用いて造影アキシャル画像Ａ１，Ａ２，…，Ａ６のそれぞれから造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。その後、作成した造影領域３次元データに基づいて、画面上に、血管αに相当する造影領域３次元画像（図１１のＩに相当）を表示する。
【００１５】
以上のＸ線ＣＴ装置１００によれば、造影領域α’内で造影剤Ｍが不均一に分布した場合でも、スライス位置ごとに適切な閾値Ｔｈ１，Ｔｈ２，…，Ｔｈ６を自動決定することが出来る。
【００１６】
−第２の実施形態−
図９は、本発明の第２の実施形態にかかる造影領域３次元データ作成処理を示すフロー図である。
ステップＵ１では、造影剤を注入しないで被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（プレーンアキシャル画像）を撮影する。
ステップＳＴ２では、造影剤を注入して被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する。
ステップＵ２では、前記プレーンアキシャル画像を基にして演算により被検体の体軸方向のスカウト画像（プレーンスカウト画像）を作成する。この作成は、先に図６を参照して説明したように、レイサム演算処理により行う（プレーンアキシャル画像を基にする点は、上記第１の実施形態と異なる）。
ステップＳＴ３では、造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する。
ステップＳＴ４では、プレーンスカウト画像と造影スカウト画像の差分をとって、差分画像を作成する。
ステップＳＴ５では、差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する。
ステップＳＴ６では、造影アキシャル画像のそれぞれから造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。
以上の造影領域３次元データ作成処理によれば、画質の高いプレーンアキシャル画像を基にしてプレーンスカウト画像を作成するので、差分画像内に造影領域（図８のα’に相当）をいっそう正確に描出でき、造影領域３次元データの精度をいっそう高めることが出来る。
【００１７】
【発明の効果】
本発明の閾値決定方法によれば、造影剤を被検体に注入しない状態で得られたプレーンスカウト画像（またはプレーンアキシャル画像を基にして作成されたプレーンスカウト画像）と、造影剤を被検体に注入した状態で得られた造影アキシャル画像を基にして作成された造影スカウト画像との差分画像を作成し、当該差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定するので、造影領域内で造影剤が不均一に分布した場合でも、スライス位置ごとに適切な閾値を自動決定することが出来る。
また、本発明の造影領域３次元データ作成方法によれば、上記閾値決定方法により決定した閾値を用いて、各造影アキシャル画像から本来の造影領域を好適に抽出でき、高画質の３次元画像を表示し得る造影領域３次元データを作成することが出来る。
【図面の簡単な説明】
【図１】本発明の第１の実施形態にかかるＸ線ＣＴ装置を示すブロック図である。
【図２】図１のＸ線ＣＴ装置による造影領域３次元データ作成処理を示すフロー図である。
【図３】プレーンスカウト画像を撮影する状態を示す説明図である。
【図４】プレーンスカウト画像の例示図である。
【図５】造影アキシャル画像を撮影する状態を示す説明図である。
【図６】造影アキシャル画像を基にしてレイサム演算処理により造影スカウト画像を作成する状態を示す説明図である。
【図７】造影スカウト画像の例示図である。
【図８】差分画像の例示図である。
【図９】本発明の第２の実施形態にかかる造影領域３次元データ作成処理を示すフロー図である。
【図１０】従来のＸ線ＣＴ装置で造影アキシャル画像を撮影する状態を示す説明図である。
【図１１】造影領域３次元画像の例示図である。
【符号の説明】
１００ Ｘ線ＣＴ装置
１ 操作コンソール
２ 入力装置
３ 中央処理装置
４ 制御インタフェース
１０ テーブル装置
２０ 走査ガントリ
２１ Ｘ線管
２３ 検出器
Ａ１，Ａ２，…，Ａ６ 造影アキシャル画像
ＤＳ 差分画像
Ｈ 被検体
Ｌ１，Ｌ２，…、Ｌ６ スライス位置
Ｍ 造影剤
ＲＳ 造影スカウト画像
Ｔｈ１，Ｔｈ２，…，Ｔｈ６ 閾値
α 血管
α’ 造影領域
【発明の属する技術分野】
本発明は、閾値決定方法および造影領域３次元データ作成方法に関し、さらに詳しくは、造影領域を抽出する基となる複数の造影アキシャル画像の体軸方向の長さが大きいときでも適切な閾値を容易に決定できる閾値決定方法および当該閾値を用いて造影領域３次元データを作成する造影領域３次元データ作成方法に関する。
【０００２】
【従来の技術】
従来のＸ線ＣＴ装置で、被検体内の血管を３次元表示する手順について説明する。
（１）図１０に示すように、被検体Ｈの造影対象の血管αに造影剤Ｍを注入した後、息止めした被検体Ｈの体軸方向のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応したアキシャル画像（造影アキシャル画像）Ａ１，Ａ２，…，Ａ６を連続して撮影する。
（２）造影アキシャル画像Ａ１，Ａ２，…，Ａ６における血管αのＣＴ値（CT１'，CT2'，CT3u'，CT3d'，…，CT6u'，CT6d'）に基づいて、造影領域を抽出するための１つの閾値Ｔｈを決定する。
（３）造影アキシャル画像Ａ１，Ａ２，…，Ａ６から、前記閾値Ｔｈに基づいて造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。
（４）図１１に示すように、作成した３次元データに基づいて、画面上に、血管αに相当する造影領域３次元画像Ｉを表示する。
【０００３】
【発明が解決しようとする課題】
旧式のＸ線ＣＴ装置では、スキャン速度が遅いので、１回の息止めで撮影できる体軸方向の長さＤ’（図１０参照）はせいぜい２０ｃｍ程度であった。この程度の長さＤ’では、血管αの全域で造影剤Ｍの濃度をほぼ均一と見なすことができ、１つの閾値Ｔｈでも造影領域の抽出を支障なく行えた。
【０００４】
これに対し、最近のＸ線ＣＴ装置では、多層型検出器の採用等によりスキャン速度が高速化されているので、１回の息止めで撮影できる体軸方向の長さＤ’が、例えば４０ｃｍ以上に増大してきている。ところが、このように長さＤ’が増大すると、血管αの区間によって造影剤Ｍの濃度が著しく異なる（例えば最上流側と最下流端とで異なる）ため、１つの閾値Ｔｈでは、本来の造影領域を正しく抽出できない問題点がある。すなわち、血管αの最小ＣＴ値に合わせて閾値Ｔｈを設定すると、閾値Ｔｈが過小となって、臓器などの非造影対象が誤って抽出されてしまう。また、血管αの最大ＣＴ値に合わせて閾値Ｔｈを設定すると、ＣＴ値の小さい血管αの区間を正しく抽出できなくなってしまう。
【０００５】
なお、上記問題点を解消するため、操作者が、造影アキシャル画像Ａ１，Ａ２，…，Ａ６ごとに血管αのＣＴ値を把握して、造影アキシャル画像ごとに閾値を個別に手動設定することも考えられるが、手間と時間がかかる問題点がある。
【０００６】
そこで、本発明の目的は、造影領域を抽出する基となる複数の造影アキシャル画像の体軸方向の長さが大きいときでも適切な閾値を容易に決定できる閾値決定方法および当該閾値を用いて造影領域３次元データを作成する造影領域３次元データ作成方法を提供することにある。
【０００７】
【課題を解決するための手段】
第１の観点では、本発明は、造影剤を注入しないで被検体の体軸方向のスカウト画像（プレーンスカウト画像）を撮影するプレーンスカウト画像撮影過程と、造影剤を注入して被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影過程と、前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成過程と、前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成過程と、前記差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定過程とを有することを特徴とする閾値決定方法を提供する。
上記第１の観点による閾値決定方法では、造影剤を注入しないで被検体の体軸方向のスカウト（scout）画像を撮影することで、プレーンスカウト（plain scout）画像を得る。また、造影剤を注入して各スライス位置のアキシャル（axial）画像を撮影して造影アキシャル画像を得て、それらの造影アキシャル画像を基にしてプレーンスカウト画像に相当する造影スカウト画像を作成する。
プレーンスカウト画像の画素値と，造影スカウト画像の画素値は、造影剤の存在する領域すなわち造影領域のみで実質的に異なり、造影領域以外の領域で同等と見なせる。したがって、差分画像として、スカウト画像内の造影領域のみが描出される。
この結果、差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定することで、造影領域の区間によって造影剤の濃度が著しく異なる場合でも、適切な閾値をきめ細かく決定することが出来る。
【０００８】
第２の観点では、本発明は、造影剤を注入しないで被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（プレーンアキシャル画像）を撮影するプレーンアキシャル画像撮影過程と、造影剤を注入して前記被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影過程と、前記プレーンアキシャル画像を基にして演算により被検体の体軸方向のスカウト画像（プレーンスカウト画像）を作成するプレーンスカウト画像作成過程と、前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成過程と、前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成過程と、前記差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定過程とを有することを特徴とする閾値決定方法を提供する。
上記第２の観点による閾値決定方法では、造影剤を注入しないで被検体の体軸方向の各スライス位置でアキシャル画像を撮影することでプレーンスアキシャル画像を得て、それらのプレーンスアキシャル画像を基にしてプレーンスカウト画像を得る。また、造影剤を注入して各スライス位置のアキシャル画像を撮影して造影アキシャル画像を得て、それらの造影アキシャル画像を基にしてプレーンスカウト画像に相当する造影スカウト画像を作成する。
プレーンスカウト画像の画素値と，造影スカウト画像の画素値は、造影剤の存在する領域すなわち造影領域のみで実質的に異なり、造影領域以外の領域で同等と見なせる。したがって、差分画像として、スカウト画像内の造影領域のみが描出される。なお、プレーンスカウト画像と，造影スカウト画像のいずれもが、画質の高いアキシャル画像を基にして作成されているので、造影領域を正確に抽出することが出来る。
この結果、差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定することで、造影領域の区間によって造影剤の濃度が著しく異なる場合でも、適切な閾値をきめ細かく決定することが出来る。
【０００９】
第３の観点では、本発明は、上記構成の閾値決定方法により決定した閾値を用いて前記造影アキシャル画像のそれぞれから造影領域を抽出する造影領域抽出過程と、前記抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する３次元データ作成過程とを有することを特徴とする造影領域３次元データ作成方法を提供する。
上記第３の観点による造影領域３次元データ作成方法では、上記閾値決定方法により決定した閾値を用いて、造影アキシャル画像のそれぞれから造影領域を過不足なく抽出することが可能となり、精度の高い造影領域３次元データを作成することが出来る。
【００１０】
【発明の実施の形態】
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
−第１の実施形態−
図１は、本発明の第１の実施形態にかかる造影領域３次元データ作成方法を実施する高速スキャン型のＸ線ＣＴ装置１００を示すブロック図である。
このＸ線ＣＴ装置１００は、操作コンソール１と、撮影テーブル１０と、走査ガントリ２０とを具備している。
前記操作コンソール１は、操作者の指示や情報などを受け付ける入力装置２と、スキャン処理や補間処理や画像再構成処理やレイサム（ray-summation）演算処理などを実行する中央処理装置３と、制御信号などを撮影テーブル１０や走査ガントリ２０へ出力する制御インタフェース４と、走査ガントリ２０で取得したデータを収集するデータ収集バッファ５と、画像などを表示するＣＲＴ６と、各種のデータやプログラムを記憶する記憶装置７とを具備している。
前記撮影テーブル１０は、被検体を乗せて体軸方向に移動させる。
前記走査ガントリ２０は、ファンビーム方式のＸ線管２１，コリメータ２２および検出器２３と、被検体の体軸の回りにＸ線管２１や検出器２３などを回転させる回転コントローラ２４と、Ｘ線照射のタイミングや強度を調整するＸ線コントローラ２５と、データ収集部２６とを具備している。
【００１１】
図２は、このＸ線ＣＴ装置１００による造影領域３次元データ作成処理を示すフロー図である。
ステップＳＴ１では、造影剤を注入しないで被検体の体軸方向のスカウト画像（プレーンスカウト画像）を撮影する。一般に、前記スカウト画像は、コロナル（coronal）画像またはサジタル（sagital）画像である。すなわち、図３に示すように、Ｘ線管２１および検出器２３を回転させずに被検体Ｈの体軸方向に長さＤだけ直線的に相対移動させ、被検体Ｈを透過したファンビームＢのＸ線強度を検出器２３の各チャネルで検出し、当該検出強度に対応した画素値をスカウト画像平面Ｆ上に投影した如きスカウト画像を作成する。前記長さＤは、造影領域の区間に相当する長さであり、例えば４０ｃｍ程度である。図４に、プレーンスカウト画像ＰＳを例示する。造影剤を注入しない状態では、プレーンスカウト画像ＰＳ内に血管αは現れない。
【００１２】
ステップＳＴ２では、図５に示すように、被検体Ｈの血管αに造影剤Ｍを注入した後、息止めした被検体Ｈの体軸方向のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応したアキシャル画像（造影アキシャル画像）Ａ１，Ａ２，…，Ａ６を連続して撮影する。図示の都合上、スライス数を“６”としたが、スライス数は任意である。
ステップＳＴ３では、造影アキシャル画像Ａ１，Ａ２，…，Ａ６を基にして、演算により前記プレーンスカウト画像ＰＳに相当する造影スカウト画像を作成する。すなわち、図６に示すように、造影アキシャル画像Ａ１，Ａ２，…，Ａ６の各スライス面内で、ファンビームＢに相当する投影線Ｂ’を前記スカウト画像平面Ｆ上に投影したときの各投影線上にある画素に対応するＣＴ値を積算するレイサム演算処理を行い、その結果を新たな画素値とした如きレイサム画像を造影スカウト画像ＲＳとする。図７に、前記造影スカウト画像ＲＳを例示する。血管αに相当する造影領域α’の画素値は、造影剤Ｍの濃度を反映して変化する。
【００１３】
ステップＳＴ４では、プレーンスカウト画像ＰＳと造影スカウト画像ＲＳの差分をとって、差分画像を作成する。図８に、差分画像ＤＳを例示する。プレーンスカウト画像ＰＳの画素値と，造影スカウト画像ＲＳの画素値は、造影領域α’のみで実質的に異なり、造影領域α’以外の領域で同等と見なせる。したがって、臓器や骨格などの領域の画素値が相殺され、造影領域α’のみが描出される。ステップＳＴ５では、差分画像ＤＳ中のスライス位置Ｌ１，Ｌ２，…，Ｌ６に対応する位置での造影領域α’のＣＴ値に基づいて造影領域を抽出するための閾値を決定する。すなわち、図８に示すように、スライス位置Ｌ１に対応する位置でのＣＴ値CT1に基づいて閾値Ｔｈ１を決定し、スライス位置Ｌ２に対応するＣＴ値CT2に基づいて閾値Ｔｈ２を決定し、スライス位置Ｌ３に対応するＣＴ値CT3u，CT3dに基づいて閾値Ｔｈ３を決定し、スライス位置Ｌ４に対応するＣＴ値CT4u，CT4dに基づいて閾値Ｔｈ４を決定し、スライス位置Ｌ５に対応するＣＴ値CT5u，CT5dに基づいて閾値Ｔｈ５を決定し、スライス位置Ｌ６に対応するＣＴ値CT6u，CT6dに基づいて閾値Ｔｈ６を決定する。
【００１４】
ステップＳＴ６では、上記ステップＳＴ５で決定した閾値Ｔｈ１，Ｔｈ２，…，Ｔｈ６を用いて造影アキシャル画像Ａ１，Ａ２，…，Ａ６のそれぞれから造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。その後、作成した造影領域３次元データに基づいて、画面上に、血管αに相当する造影領域３次元画像（図１１のＩに相当）を表示する。
【００１５】
以上のＸ線ＣＴ装置１００によれば、造影領域α’内で造影剤Ｍが不均一に分布した場合でも、スライス位置ごとに適切な閾値Ｔｈ１，Ｔｈ２，…，Ｔｈ６を自動決定することが出来る。
【００１６】
−第２の実施形態−
図９は、本発明の第２の実施形態にかかる造影領域３次元データ作成処理を示すフロー図である。
ステップＵ１では、造影剤を注入しないで被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（プレーンアキシャル画像）を撮影する。
ステップＳＴ２では、造影剤を注入して被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する。
ステップＵ２では、前記プレーンアキシャル画像を基にして演算により被検体の体軸方向のスカウト画像（プレーンスカウト画像）を作成する。この作成は、先に図６を参照して説明したように、レイサム演算処理により行う（プレーンアキシャル画像を基にする点は、上記第１の実施形態と異なる）。
ステップＳＴ３では、造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する。
ステップＳＴ４では、プレーンスカウト画像と造影スカウト画像の差分をとって、差分画像を作成する。
ステップＳＴ５では、差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する。
ステップＳＴ６では、造影アキシャル画像のそれぞれから造影領域を抽出し、抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する。
以上の造影領域３次元データ作成処理によれば、画質の高いプレーンアキシャル画像を基にしてプレーンスカウト画像を作成するので、差分画像内に造影領域（図８のα’に相当）をいっそう正確に描出でき、造影領域３次元データの精度をいっそう高めることが出来る。
【００１７】
【発明の効果】
本発明の閾値決定方法によれば、造影剤を被検体に注入しない状態で得られたプレーンスカウト画像（またはプレーンアキシャル画像を基にして作成されたプレーンスカウト画像）と、造影剤を被検体に注入した状態で得られた造影アキシャル画像を基にして作成された造影スカウト画像との差分画像を作成し、当該差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定するので、造影領域内で造影剤が不均一に分布した場合でも、スライス位置ごとに適切な閾値を自動決定することが出来る。
また、本発明の造影領域３次元データ作成方法によれば、上記閾値決定方法により決定した閾値を用いて、各造影アキシャル画像から本来の造影領域を好適に抽出でき、高画質の３次元画像を表示し得る造影領域３次元データを作成することが出来る。
【図面の簡単な説明】
【図１】本発明の第１の実施形態にかかるＸ線ＣＴ装置を示すブロック図である。
【図２】図１のＸ線ＣＴ装置による造影領域３次元データ作成処理を示すフロー図である。
【図３】プレーンスカウト画像を撮影する状態を示す説明図である。
【図４】プレーンスカウト画像の例示図である。
【図５】造影アキシャル画像を撮影する状態を示す説明図である。
【図６】造影アキシャル画像を基にしてレイサム演算処理により造影スカウト画像を作成する状態を示す説明図である。
【図７】造影スカウト画像の例示図である。
【図８】差分画像の例示図である。
【図９】本発明の第２の実施形態にかかる造影領域３次元データ作成処理を示すフロー図である。
【図１０】従来のＸ線ＣＴ装置で造影アキシャル画像を撮影する状態を示す説明図である。
【図１１】造影領域３次元画像の例示図である。
【符号の説明】
１００ Ｘ線ＣＴ装置
１ 操作コンソール
２ 入力装置
３ 中央処理装置
４ 制御インタフェース
１０ テーブル装置
２０ 走査ガントリ
２１ Ｘ線管
２３ 検出器
Ａ１，Ａ２，…，Ａ６ 造影アキシャル画像
ＤＳ 差分画像
Ｈ 被検体
Ｌ１，Ｌ２，…、Ｌ６ スライス位置
Ｍ 造影剤
ＲＳ 造影スカウト画像
Ｔｈ１，Ｔｈ２，…，Ｔｈ６ 閾値
α 血管
α’ 造影領域
Claims (3)
- 造影剤を注入しないで被検体の体軸方向のスカウト画像（プレーンスカウト画像）を撮影するプレーンスカウト画像撮影手段と、
造影剤を注入して被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影手段と、
前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成手段と、
前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成手段と、
前記差分画像中の前記スライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定手段と
を有することを特徴とするＸ線ＣＴ装置。 - 造影剤を注入しないで被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（プレーンアキシャル画像）を撮影するプレーンアキシャル画像撮影手段と、
造影剤を注入して前記被検体の体軸方向の複数の異なるスライス位置でアキシャル画像（造影アキシャル画像）を撮影する造影アキシャル画像撮影手段と、
前記プレーンアキシャル画像を基にして演算により被検体の体軸方向のスカウト画像（プレーンスカウト画像）を作成するプレーンスカウト画像作成手段と、
前記造影アキシャル画像を基にして演算により前記プレーンスカウト画像に相当する造影スカウト画像を作成する造影スカウト画像作成手段と、
前記プレーンスカウト画像と前記造影スカウト画像の差分をとって差分画像を作成する差分画像作成手段と、
前記差分画像中の前記造影アキシャル画像のスライス位置に対応する位置での造影領域の画素値に基づいて造影領域を抽出するための閾値を決定する閾値決定手段と
を有することを特徴とするＸ線ＣＴ装置。 - 前記閾値決定手段により決定した閾値を用いて前記造影アキシャル画像のそれぞれから造影領域を抽出する造影領域抽出手段と、
前記抽出した造影領域のデータを基に３次元データ（造影領域３次元データ）を作成する３次元データ作成手段と
を有することを特徴とする請求項１または請求項２に記載のＸ線ＣＴ装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13920099A JP4350206B2 (ja) | 1999-05-19 | 1999-05-19 | Ｘ線ｃｔ装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13920099A JP4350206B2 (ja) | 1999-05-19 | 1999-05-19 | Ｘ線ｃｔ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000325340A JP2000325340A (ja) | 2000-11-28 |
| JP4350206B2 true JP4350206B2 (ja) | 2009-10-21 |
Family
ID=15239897
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP13920099A Expired - Fee Related JP4350206B2 (ja) | 1999-05-19 | 1999-05-19 | Ｘ線ｃｔ装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4350206B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6470070B2 (en) * | 2000-12-20 | 2002-10-22 | Cedara Software Corp. | Image reconstruction using multiple X-ray projections |
| JP4647360B2 (ja) * | 2004-04-05 | 2011-03-09 | 富士フイルム株式会社 | 差分画像作成装置、差分画像作成方法、及び、そのプログラム |
| JP5670025B2 (ja) * | 2009-02-04 | 2015-02-18 | 株式会社東芝 | Ｘ線コンピュータ断層撮影装置 |
| JP2013223809A (ja) * | 2013-08-07 | 2013-10-31 | Toshiba Corp | Ｘ線コンピュータ断層撮影装置 |
-
1999
- 1999-05-19 JP JP13920099A patent/JP4350206B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2000325340A (ja) | 2000-11-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4965433B2 (ja) | トランケートされた投影と事前に取得した３次元ｃｔ画像を用いる円錐ビームｃｔ装置 | |
| US6487432B2 (en) | Method and system for selecting and displaying medical image data | |
| US5852646A (en) | X-ray imaging method | |
| EP2501290B1 (en) | Scan plan field of view adjustor, determiner, and/or quality assessor | |
| Zijp et al. | Extraction of the respiratory signal from sequential thorax cone-beam X-ray images | |
| JP2009034503A (ja) | トモシンセシス画像を表示するための方法及びシステム | |
| US20210104040A1 (en) | System and method for automated angiography | |
| US10748263B2 (en) | Medical image processing apparatus, medical image processing method and medical image processing system | |
| US9532750B2 (en) | Object identification in dual energy contrast-enhanced CT images | |
| JP2005323628A (ja) | 画像再構成方法およびｘ線ｃｔ装置 | |
| CN111184523A (zh) | 基于dr设备的三维图像重建方法及*** | |
| US5561695A (en) | Methods and apparatus for reducing image artifacts | |
| US7260175B2 (en) | CT image producing method and X-ray CT apparatus | |
| US11771390B2 (en) | Method and device for determining the contour of anatomical structures in a digital X-ray-based fluoroscopic image | |
| JP4350206B2 (ja) | Ｘ線ｃｔ装置 | |
| JPH07148155A (ja) | コンピュータ断層撮影装置 | |
| US6834096B2 (en) | Apparatus and method for providing a three-dimensional radiographic image of an object | |
| WO2009004571A1 (en) | Method and apparatus for image reconstruction | |
| Mol et al. | Application of digital image analysis in dental radiography for the description of periapical bone lesions: a preliminary study | |
| JP2008154680A (ja) | Ｘ線ｃｔ装置 | |
| JP3647318B2 (ja) | ツイン非対称スキャン用スライス厚設定方法およびｘ線ｃｔ装置 | |
| JP2002200073A (ja) | Ｘ線ｃｔ装置 | |
| JP2005323627A (ja) | 撮影計画作成方法およびｘ線ｃｔ装置 | |
| EP1542166A1 (en) | Three-dimension back-projection method and apparatus, and x-ray CT apparatus | |
| JPH04249746A (ja) | 断層像再構成法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060417 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20090413 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090421 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090618 |
|
| 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: 20090714 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090722 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120731 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130731 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130731 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |