(57)【特許請求の範囲】
【請求項１】
造影剤の流入前後の被写体のＸ線画像から複数のフレームからなるサブトラクション動画像を生成する減算手段と、
前記複数のフレームからなるサブトラクション動画像の各フレームについて、血管領域における評価値を算出する評価手段と、
前記各フレームに対する評価値の比較に基づいて、前記複数のフレームからなるサブトラクション動画像から少なくとも１つのフレームを選択し、前記選択されたフレームを参照画像として出力する選択手段と、
を有することを特徴とする画像処理装置。
【請求項２】
前記減算手段は、造影剤の流入後の複数フレームからなるＸ線画像であるライブ像から造影剤の流入前の前記被写体のＸ線画像から生成されるマスク像を減算処理して前記サブトラクション動画像を生成することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記評価手段は、前記サブトラクション動画像の各フレームから、閾値に対して所定範囲の画素値を持つ画素の数を用いた評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
前記評価手段は、前記サブトラクション動画像の各フレームから、閾値に対して所定範囲の画素値を持ち、かつ連続した配置の画素の数を用いた評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項５】
前記評価手段は、前記サブトラクション動画像に含まれるモーションアーチファクトから、閾値に対して所定範囲の画素値を持つ画素の数を用いて評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項６】
前記評価手段は、前記サブトラクション動画像の各フレームから、前記造影剤と背景領域との境界画素数を用いて評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項７】
前記評価手段は、前記サブトラクション動画像を生成した際の撮影部位、撮影手技、撮影条件、及び造影剤の注入からの経過時間に基づく評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項８】
前記評価手段は、前記撮影条件から前記サブトラクション動画像のフレーム範囲を求め、フレーム範囲内で前記評価値を算出することを特徴とする請求項７に記載の画像処理装置。
【請求項９】
前記評価手段は、前記サブトラクション動画像の各フレーム間で対応する各画素のうち、最小の画素値を有する画素によって構成される最小値投影像と、前記サブトラクション動画像との比較から評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項１０】
前記評価手段は、前記サブトラクション動画像の各フレームから複数の特徴量を算出し、算出した複数の特徴量の加算に基づき評価値を算出することを特徴とする請求項１または２に記載の画像処理装置。
【請求項１１】
前記評価手段は、前記サブトラクション動画像の各フレームから求めた複数の特徴量に重み係数を乗じた値を線形結合により加算した値を評価値として算出することを特徴とする請求項１０記載の画像処理装置。
【請求項１２】
前記参照画像を表示する画像表示手段を更に備えることを特徴とする請求項１乃至１１のいずれか１項に記載の画像処理装置。
【請求項１３】
前記画像表示手段は、前記参照画像が複数フレームについて選択された場合、所定のフレームレートで切り替えて前記参照画像を表示することを特徴とする請求項１２に記載の画像処理装置。
【請求項１４】
前記参照画像を記憶する記憶手段を更に備えることを特徴とする請求項１乃至１３のいずれか１項に記載の画像処理装置。
【請求項１５】
前記評価手段は、前記複数のフレームからなるサブトラクション動画像の各フレームについて、前記造影剤の流入状況を評価することにより前記評価値を算出することを特徴とする請求項１乃至１４のいずれか１項に記載の画像処理装置。
【請求項１６】
前記評価手段は、前記複数のフレームからなるサブトラクション動画像の各フレームについて、血管領域を評価することにより前記評価値を算出することを特徴とする請求項１乃至１４のいずれか１項に記載の画像処理装置。
【請求項１７】
前記選択手段は、前記複数のフレームのうち血管領域の面積が最大となるフレームを、前記参照画像として選択することを特徴とする請求項１６に記載の画像処理装置。
【請求項１８】
画像処理装置の減算手段が、造影剤の流入前後の被写体のＸ線画像から複数のフレームからなるサブトラクション動画像を生成する減算工程と、
前記画像処理装置の評価手段が、前記複数のフレームからなるサブトラクション動画像の各フレームについて、血管領域における評価値を算出する評価工程と、
前記画像処理装置の選択手段が、前記各フレームに対する評価値の比較に基づいて、前記複数のフレームからなるサブトラクション動画像から少なくとも１つのフレームを選択し、前記選択されたフレームを参照画像として出力する選択工程と、
を有することを特徴とする画像処理方法。
【請求項１９】
コンピュータを、請求項１乃至１７のいずれか１項に記載の画像処理装置の各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法およびプログラムに関する。特に複数のＸ線ディジタルサブトラクションアンギオグラフィ像から診断や手術に好適な画像を選択する画像処理技術に関するものである。
【背景技術】
【０００２】
近年のディジタル技術の進歩により、医療の分野でも画像にディジタル処理を施すことが一般的になっている。従来のＸ線診断用のフィルムを用いたＸ線撮影に代わり、Ｘ線画像をディジタル画像として出力する２次元Ｘ線センサの普及も進み、２次元Ｘ線センサが出力するディジタル画像へのディジタル画像処理の応用範囲は広い。
【０００３】
このディジタル画像処理の好適な適用例として、ディジタルサブトラクションアンギオグラフィ像（以下、ＤＳＡ像）を取得するＤＳＡ処理が挙げられる。ＤＳＡ像は、被写体への造影剤流入の前後で画像を取得し、造影剤流入後の画像（以下、ライブ像）から造影剤流入前の画像（以下、マスク像）を減算処理して得られるサブトラクション像である。ＤＳＡ像は、診断上の関心領域である血管領域を造影剤流入の前後で取得した画像間の変化領域として保持する一方、変化領域以外の不要な領域を背景領域として取り除いて均一な領域として表現する。即ち、ＤＳＡ像は、ライブ像では背景領域によって埋もれることで視認性を落とす造影血管領域を、サブトラクション像として明瞭に描出する。このＤＳＡ像は、一般に血管インターベンショナル・ラジオロジー（以下、ＩＶＲ）において用いられる。血管ＩＶＲは、Ｘ線透視像を見ながら体内の目的臓器にカテーテルを挿入する手法である。
【０００４】
ＤＳＡ処理に関しては、例えば、特許文献１に、外科手術等において血管走行状態を容易に把握できるオリエンテーション画像を得るのに好適なＸ線画像処理装置としての発明の開示がある。
【０００５】
通常、Ｘ線透視像では、Ｘ線透過率の高い血管を描出することはできない。そこで、血管を描出するために造影剤が投与される。しかし、血管ＩＶＲにおいてカテーテルを目的の臓器へ到達させるまで造影剤を流し続けると大量の造影剤を投与することになり、患者への負担が大きい。そこで、一般に血管ＩＶＲでは、透視の前に所定量の造影剤を投与して連続撮影を行い、造影後の連続画像から造影前の画像をそれぞれ減算することでサブトラクション像を動画として取得し、ＤＳＡ動画像とする。そして、ＤＳＡ動画像の中からカテーテルを進めるための血管が明瞭に描出されたフレーム（以下、参照画像）を選択する。この参照画像は、透視像を表示するメインモニタとは別のセカンドモニタ上に表示され、透視時にカテーテルを進める際の手掛かりとして用いられる。
【０００６】
一般的な血管ＩＶＲでは、目的の臓器にカテーテルが到達するまで、透視によるカテーテル操作、造影剤を投与しながらの連続撮影、得られた画像からのＤＳＡ動画像の生成、参照画像として用いるフレームの選択を繰り返し行う。透視化で造影剤を投与しながらカテーテルを操作するのではなく、透視前に取得した参照画像を用いることで、造影剤投与を抑えることができる。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開平１−２２１１４２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、従来、ＤＳＡ動画像からの参照画像の選択作業は、術者がＤＳＡ動画を再生、コマ送り、コマ戻しなどして目視で確認し、最適なフレームを探していた。一般に血管ＩＶＲでは、一度のＤＳＡ動画像作成、一度の透視によるカテーテル操作で目的臓器に到達するものではない。即ち、ＤＳＡ動画像作成及び参照画像の選択は、カテーテルの挿入状況に応じて複数回繰り返される。従って、この参照画像の選択作業は、術者にとって非常に煩雑な作業であり、処置の長時間化にもつながる。
【課題を解決するための手段】
【０００９】
本発明は、上記の問題を考慮して、血管ＩＶＲにおけるセカンドモニタに表示すべき参照画像の選択作業の労力を削減することを可能にする画像処理技術の提供を目的とする。
【００１０】
本発明の一つの側面に係る画像処理装置は、造影剤の流入前後の被写体のＸ線画像から複数のフレームからなるサブトラクション動画像を生成する減算手段と、
前記複数のフレームからなるサブトラクション動画像の各フレームについて、血管領域における評価値を算出する評価手段と、
前記各フレームに対する評価値の比較に基づいて、前記複数のフレームからなるサブトラクション動画像から少なくとも１つのフレームを選択し、前記選択されたフレームを参照画像として出力する選択手段と、
を有することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、血管ＩＶＲにおけるセカンドモニタに表示すべき参照画像の選択作業の労力を削減することが可能になる。
【図面の簡単な説明】
【００１２】
【図１】実施形態にかかる画像処理装置の構成を説明する図。
【図２】実施形態にかかる画像処理方法の流れを説明する図。
【図３】サブトラクション像を模式的に説明する図。
【図４】最小値投影像を模式的に説明する図。
【図５】実施形態にかかる情報処理装置の構成を説明する図。
【発明を実施するための形態】
【００１３】
以下、図面を参照して、本発明の好適な実施形態を例示的に説明する。
【００１４】
（第１実施形態）
本発明の第１実施形態にかかる画像処理装置（以下、「Ｘ線画像処理装置１００」）を図１の参照により説明する。Ｘ線画像処理装置１００は、Ｘ線発生部１０１、２次元Ｘ線センサ１０４、前処理部１０５、画像記憶部１０６、画像間減算部１０７、評価部１０８、参照画像選択部１０９、および画像表示部１１０を有する。
【００１５】
Ｘ線発生部１０１は、毎秒３〜３０パルスのＸ線パルス（Ｘ線）を発生することが可能である。２次元Ｘ線センサ１０４は、動画を撮影することが可能であり、被写体１０２を透過したＸ線１０３を受けてＸ線パルスに同期した動画像を撮像する。２次元Ｘ線センサ１０４は、Ｘ線照射された被写体１０２の動画像を撮像する撮像部として機能する。
【００１６】
前処理部１０５は、２次元Ｘ線センサ１０４が出力する動画像の各フレームに対して前処理を行う。また、画像記憶部１０６は、前処理部１０５によって前処理された動画像のうち少なくとも１フレームを造影剤流入前のマスク像として記憶する。マスク像として記憶するフレームは、例えば、動画像撮像開始直後のフレームや、造影剤の流入を動画像から検出して自動取得した造影剤流入の直前フレームである。あるいは操作者が造影剤流入開始時に記憶タイミングを指示して選択したフレームである。また、画像記憶部１０６に複数のフレームを記憶しておき、記憶してある複数のフレームの中からマスク像として用いるフレームを適宜選択するか、あるいはマスク像として複数のフレームを合成したものでも良い。
【００１７】
画像間減算部１０７は、造影剤の流入前後の被写体のＸ線画像から複数のフレームからなる画像を生成する。すなわち、画像間減算部１０７は、前処理部１０５が出力する造影剤流入後の画像（フレーム）であるライブ像から、画像記憶部１０６が記憶するマスク像を減算したサブトラクション像を生成し、ＤＳＡ動画像の各フレームとして出力する。
【００１８】
評価部１０８は、画像間減算部１０７が出力するＤＳＡ動画像の各フレームから、各フレームの造影剤流入状況を示す評価値を算出し、出力する。
【００１９】
参照画像選択部１０９は、評価部１０８が出力する各フレームの評価値に基づいて、ＤＳＡ動画像から少なくとも１フレームを選択し、参照画像として出力する。
【００２０】
画像表示部１１０は、参照画像選択部１０９が出力する参照画像を表示する。ここで画像表示部１１０は、参照画像が複数フレームについて選択された場合、所定のフレームレートで複数の参照画像の表示を切り替えながらループ再生することが可能である。また、参照画像の出力先として、参照画像選択部１０９は参照画像を、不図示の記憶部に保存することも可能である。画像表示部１１０は、記憶部に記憶された参照画像を表示することも可能である。
【００２１】
以上の構成を持つＸ線画像処理装置１００に関して、本実施形態の特徴的な構成である評価部１０８及び画像表示部１１０の動作を、図２に示すフローチャートを用いて説明する。
【００２２】
ステップＳ２０１において、画像間減算部１０７はＤＳＡ動画像の各フレームを評価部１０８に対して出力し、評価部１０８はＤＳＡ動画像の各フレームの造影剤流入状況を評価した評価値を出力する。造影剤流入状況の評価は、様々な手法を適用することができる。ここでは、例として、ＤＳＡ動画像の各フレーム（ＤＳＡ像）の血管領域の面積に基づいて評価する手法を説明する。
【００２３】
ライブ像からマスク像を減算処理したサブトラクション像の模式図を図３に示す。図３に示すサブトラクション像は、−１．０〜１．０の範囲に画素値を持つ画像であるとし、１．０に近い画素値を白色、−１．０に近い画素値を黒色で表わしている。ＤＳＡ像は造影剤の流入後の複数フレームからなるＸ線画像であるライブ像から造影剤の流入前の被写体のＸ線画像から生成されるマスク像を減算処理したサブトラクション像である。従って、理想的なサブトラクション像では、造影剤が投入（流入）された血管領域３０１を除いた背景領域３０２は、減算処理によりキャンセルされ、画素値は０となる。さらに血管領域３０１は、背景領域３０２の画素値に対して片側（この例では、負値の側）に集中する。
【００２４】
従って、ＤＳＡ像の血管領域３０１の面積により造影剤流入状況を評価する手法では、画素値として負値を持つ画素の数を数え上げる手法をとる。ここで、総フレーム数Ｔ、各フレームの縦方向画素数Ｈ、横方向画素数ＷのＤＳＡ動画像のｔ(０≦ｔ＜Ｔ)フレーム目を、サブトラクション像Ｉtとする。更に、サブトラクション像Ｉtの画素（ｘ,ｙ） (０≦ｘ＜Ｗ,０≦ｙ＜Ｈ)の画素値をＩt(ｘ,ｙ)とする。更に、画素値を０とした背景領域３０２と血管領域３０１とを分離する閾値を一般化してthreshとする。以上の関係を用いると、ＤＳＡ動画像のｔフレーム目の評価値Ｓｖtは、次式で表わすことができる。
【００２５】
【数１】
0006021420.tif 000002
【００２６】
ステップＳ２０２において、評価部１０８は評価値を参照画像選択部１０９に対して出力し、参照画像選択部１０９は、評価部１０８から出力された評価値に基づいて、参照画像として好適なＤＳＡ像を少なくとも１フレーム選択する。
【００２７】
上述した例では、ステップＳ２０１において、評価部１０８は造影剤流入状況をＤＳＡ像の血管領域の面積で評価し、ＤＳＡ動画像のｔ(０≦ｔ＜Ｔ)フレーム目の評価値として評価値Ｓｖtを算出した。
【００２８】
造影剤が血管領域において広範囲にいきわたった状態を示す画像が参照画像として好適である。従って、参照画像選択部１０９が選択するフレームｔrefは、評価部１０８から出力された各フレームの評価値Ｓｖのうち、最大値となるフレームであり、評価部１０８から出力された評価値を用いると次式で表わすことができる。
【００２９】
【数２】
0006021420.tif 000003
【００３０】
上式は、評価値Ｓｖtが最大値となるフレームｔを、参照画像選択部１０９が選択する参照画像（フレームｔref）とすることを意味している。上式において、参照画像選択部１０９は、評価値Ｓｖtが最大値となるｔ(０≦ｔ＜Ｔ)フレーム目の１フレームを参照画像として選択することを想定している。しかし、参照画像選択部１０９が選択する参照画像は評価値Ｓｖの最大値に従って１フレームの画像に限定されるものではなく、評価値Ｓｖを用いて複数のフレームの画像を選択することも可能である。例えば、参照画像を２枚以上の複数フレームとし、所定のフレームレートで画像を切り替えながらループ再生する場合は、評価値Ｓｖが最大となるフレームの前後の複数フレームを選択するよう構成しても良い。または、評価値Ｓｖが高い方から順番に評価値Ｓｖが低くなる複数フレームを選択するよう構成しても良い。
【００３１】
本実施形態によれば、評価部１０８は、ＤＳＡ動画像の各フレーム（ＤＳＡ像）を、ＤＳＡ像の血管領域の面積で評価し、参照画像選択部１０９は、血管領域の面積が最大となるフレームを参照画像として選択（出力）する。この参照画像を画像表示部１１０に出力することで、Ｘ線画像処理装置の操作者によるフレーム選択の手間を削減することができる。
【００３２】
上述の例では、各フレーム（ＤＳＡ像）の画素の有する画素値と閾値との比較により、所定範囲の画素値を持つ画素の数を用いた評価値を説明したが、ＤＳＡ像には造影血管らしい画素値を持つ細かなノイズが存在する。従って、評価部１０８は、所定範囲の画素値を持つ画素数に加えて、更に、画素の連続性を考慮した値を評価値としてもよい。例えば、閾値との比較により所定範囲の画素値を持ち、かつ、その条件を満たした画素が連続して配置される数を評価値としてもよい。造影血管領域はカテーテルの先端から連続した領域であると考えられるため、この評価値を用いることで、参照画像の選択精度の更なる向上を図ることができる。
【００３３】
また、図３に示す血管領域３０１と背景領域３０２の例では、理想的なサブトラクション像に対する評価値として、ＤＳＡ像の血管領域３０１の面積を評価値とする例を説明した。しかし、実際のサブトラクション像では、被写体動や挿入されたカテーテルの動きによりモーションアーチファクト３０３が生じ得る。モーションアーチファクト３０３の画素値は、例えば図３に示すとおり、正負に渡って生じるため、ＤＳＡ像の画質に大きく影響する。参照画像としてはモーションアーチファクト３０３が小さいことが望ましい。モーションアーチファクト３０３の存在を加味して参照画像を選択するために、ステップＳ２０１において評価部１０８が算出する評価値に、モーションアーチファクトを考慮した値を加えても良い。
【００３４】
モーションアーチファクトを考慮した評価値の算出は様々な手法が考えられる。画素値が負値側に集中する血管領域３０１とは異なり、モーションアーチファクト３０３の画素値は、図３に示すように、正負両方の値をとり得る。モーションアーチファクト３０３には、例えば、１．０に近い画素値を示す領域（白色で示す領域）と、−１．０に近い画素値を示す領域（黒色で示す領域）とが含まれ得る。ここでは簡易な例として、画素値として正の値を持つ画素の数を数え上げる手法を説明する。これは、数式１と同様の記号を用い、背景領域３０２とモーションアーチファクト３０３とを分離する閾値をthresh２とすれば、モーションアーチファクト３０３に基づく評価値Ｓｍtは、次式で表わすことができる。
【００３５】
【数３】
0006021420.tif 000004
【００３６】
このモーションアーチファクトに基づく評価値Ｓｍに血管領域３０１の面積に基づく評価値Ｓｖを加算した評価値を、最終的に算出すべき評価値Ｓ（＝Ｓｍ＋Ｓｖ）とする。さらに、モーションアーチファクト３０３に基づく評価値Ｓｍを用いる場合、モーションアーチファクト３０３は小さいほうが参照画像として好適である。これを踏まえ、評価値Ｓｖに対する重みとして０より大きい係数ａと、評価値Ｓｍに対する重みとして０より小さい係数ｂを用いて、評価値Ｓｖと評価値Ｓｍとの和として算出される評価値Ｓを次式で表わす。
【００３７】
【数４】
0006021420.tif 000005
【００３８】
ＤＳＡ動画像のｔ(０≦ｔ＜Ｔ)フレーム目の評価値として評価値Ｓtを用いると、ステップＳ２０２における参照画像選択部１０９が選択するフレームｔrefは、次式で表わすことができる。
【００３９】
【数５】
0006021420.tif 000006
【００４０】
上式は、式２と同様、評価値Ｓtが最大値となるフレームtを、参照画像選択部１０９が選択する参照画像（ｔref）とすることを意味している。即ち、参照画像選択部１０９は、血管領域の面積が大きく、かつモーションアーチファクトが小さいフレームを参照画像として選択する。
【００４１】
また、評価値には、造影剤と背景領域の境界画素数を用いても良い。造影剤と背景領域の境界は、例えば、サブトラクション像にＣａｎｎｙ法などによるエッジ検出処理を施し、検出したエッジ位置での画素値を調べれば、造影剤と背景領域との境界を判別することができる。この境界画素の数を評価値に加味すれば、上述した血管領域の面積だけでは評価できない細い血管構造の複雑度を評価値に加えることができ、より好適な参照画像を選択することが可能となる。
【００４２】
更に、評価部１０８は、ＤＳＡ動画像撮影時の撮影部位、撮影手技、撮影条件、造影剤の注入からの経過時間に基づく評価値を算出することも可能である。参照画像として選択される画像は造影剤が被写体にいきわたったタイミングにおけるＤＳＡ像である。従って造影剤の注入直後や造影剤注入から長時間が経過したＤＳＡ像は参照画像としては使えない可能性が高い。例えば、撮影部位や撮影手技で定まる造影剤の投与量をＱ、撮影部位による血流の速さをＤ、撮影条件により定まるフレームレートをＲとする。これらの条件を変えた場合の参照画像として利用可能なフレーム範囲[tstart,tend]を、予め実験的に求めてテーブルとして保持する。
【００４３】
【数６】
0006021420.tif 000007
【００４４】
実際の撮影では、評価部１０８は、このテーブルから撮影条件によって定められるフレーム範囲[tstart,tend]を求め、そのフレーム範囲内で評価値Ｓtを算出する。そして、参照画像選択部１０９は、評価部１０８によって算出された評価値に基づき参照画像としてフレームｔrefを選択する。この構成によれば、フレームの探索範囲を制限できるため処理の高速化が可能となる。また、テーブルに保持するデータとしては、利用可能なフレーム範囲に限定されるものではない。例えば、造影剤の投与量Ｑ、撮影部位による血流の速さＤ、撮影条件により定まるフレームレートＲの条件に応じた重みｗtとして、０．０〜１．０の値を保持するようにテーブルを構成してもよい。例えば、以下の式に示すように、造影剤の投与量Ｑ、血流の速さＤ、フレームレートＲが定められると、これらに対応した重みｗtを求めることができる。
【００４５】
【数７】
0006021420.tif 000008
【００４６】
この重みｗtを用いると、ステップＳ２０２における参照画像選択部１０９が選択するフレームｔrefは次式で表わすことができる。
【００４７】
【数８】
0006021420.tif 000009
【００４８】
以上の構成によると、造影剤の投与量Ｑ、血流の速さＤ、フレームレートＲに応じてフレーム毎に評価値Ｓtに重みを与えることで、より精度の高い参照画像の選択が可能になる。
【００４９】
更に、ＤＳＡ動画像から最小値投影像Ｍを作成し、この最小値投影像とＤＳＡ動画像の各フレームを比較した結果を評価値として用いても良い。最小値投影像Ｍは、ＤＳＡ動画像の各フレームをＩt、画素(ｘ,ｙ)の値をＩt (ｘ,ｙ)で表わした時、次式で表わされる画像である。
【００５０】
【数９】
0006021420.tif 000010
【００５１】
最小値投影像Ｍの各画素(ｘ,ｙ)の画素値Ｍ(ｘ,ｙ)は、ＤＳＡ動画像の全フレームにおける各フレーム間で対応する画素の最小値である。図４は、ＤＳＡ像４０１、４０２、４０３から作成した最小値投影像４００を模式的に示した図である。図４に示す通り、ＤＳＡ動画像の各フレームであるＤＳＡ像４０１、４０２、４０３は、時間経過に伴い、造影血管領域の描画位置が移動するため、それぞれのＤＳＡ像（フレーム）では全血管領域の一部ずつしか描出されていない。これらのフレームから構成される最小値投影像４００は、各フレーム間で対応する各画素のうち、最小の画素値を有する画素を取得するによって構成される。ＤＳＡ動画像を構成するフレームのどれかに造影血管領域がある場合、その造影血管領域を持つ画素値が最小値として取得されるため、最小値投影像４００は全血管領域を描出することができる。この最小値投影像ＭとＤＳＡ動画像の各フレームとの比較による評価値Ｓｃtを次式で算出する。
【００５２】
【数１０】
0006021420.tif 000011
【００５３】
評価値Ｓｃtは、最小値投影像Ｍに最も近い、即ち、造影剤が広い範囲にいきわたっているＤＳＡ動画像のフレームにおいて最小になる。これを踏まえ、モーションアーチファクトに基づく評価値Ｓｍと同様に０より小さい所定の係数を乗算し、評価値に加えることで、より好適な参照画像を選択することが可能となる。
【００５４】
以上、いくつかの評価値の算出を例として説明したが、本発明で用いることができる評価値は、上述の例に限定されるものではない。例えば、評価部１０８は、サブトラクション画像の各フレームから複数の特徴量を算出し、算出した複数の特徴量の加算に基づき評価値を算出することが可能である。例えば、画像Iからｎ種類の特徴量（Ｓ1,Ｓ2,…,Ｓn）を求め、重み係数（ｋ１,ｋ２,…,ｋｎ）を用いて線形結合により加算した値を評価値Ｓとする場合、次式で表わすことができる。
【００５５】
【数１１】
0006021420.tif 000012
【００５６】
重み係数ｋの値は、例えば、対応する特徴量が大きいときにその特徴量が参照画像を示すものには０より大（Ｓ＞０）とし、対応する特徴量が小さいときにその特徴量が参照画像を示すものには０より小（Ｓ＜０）とする。
【００５７】
このとき、参照画像選択部１０９は、評価値Ｓに応じて（例えば、評価値Ｓが最大となる）フレームを選択することで複数の特徴量に対応した参照画像の選択が可能となる。
【００５８】
ここで、図１に示した各部は専用のハードウェアにより構成しても良いが、ハードウェアの機能構成をソフトウェアにより実現することも可能である。この場合、図１に示す各部の機能は、情報処理装置にソフトウェアをインストールし、ソフトウェアの実行による画像処理方法を情報処理装置の演算機能を利用して実現することができる。ソフトウェアの実行により、例えば、２次元Ｘ線センサ１０４が出力する動画像の各フレームに対して減算工程、評価工程および選択工程が実行される。
【００５９】
図５は、情報処理装置のハードウェア構成、及びその周辺機器の構成を示すブロック図である。情報処理装置１０００は、撮像装置２０００と接続しており、互いにデータ通信が可能な構成となっている。
【００６０】
（情報処理装置）
ＣＰＵ１０１０は、ＲＡＭ１０２０やＲＯＭ１０３０に格納されているプログラムやデータを用いて情報処理装置１０００の全体的な制御を行うことが可能である。また、ＣＰＵ１０１０は、プログラムの実行により予め定められた画像処理に関する演算処理を実行することが可能である。
【００６１】
ＲＡＭ１０２０は、光磁気ディスク１０６０やハードディスク１０５０からロードされたプログラムやデータを一時的に記憶するためのエリアを備える。更に、ＲＡＭ１０２０は、撮像装置２０００から取得したＸ線透視の動画像データを一時的に記憶するためのエリアを備える。また、ＲＡＭ１０２０は、ＣＰＵ１０１０が各種の処理を実行する際に用いるワークエリアも備える。    
ＲＯＭ１０３０は、情報処理装置１０００の設定データや、ブートプログラムなどを格納する。
【００６２】
ハードディスク１０５０は、情報記憶媒体としての一例であり、ＯＳ（オペレーティングシステム）や、図１に示した各部が行う各処理をＣＰＵ１０１０に実行させるためのプログラムやデータを保持している。そして、これらはＣＰＵ１０１０による制御に従って適宜、ＲＡＭ１０２０にロードされ、ＣＰＵ１０１０による処理対象となる。また、Ｘ線透視の動画像のデータをこのハードディスク１０５０に保存することも可能である。
【００６３】
光磁気ディスク１０６０は、情報記憶媒体としての一例であり、ハードディスク１０５０に保存されているプログラムやデータの一部若しくは全部をこの光磁気ディスク１０６０に格納することが可能である。
【００６４】
マウス１０７０、キーボード１０８０は、情報処理装置１０００の操作者が操作することで、各種の指示をＣＰＵ１０１０に対して入力することができる。
【００６５】
プリンタ１０９０は、画像表示部１１０に表示された画像を記録媒体上に印刷出力することが可能である。
【００６６】
表示装置1１００は、ＣＲＴや液晶画面などにより構成されており、ＣＰＵ１０１０による処理結果を画像や文字などでもって表示することが可能である。例えば、図１に示した各部により処理され、最終的に画像表示部１１０から出力された画像を表示することが可能である。この場合、画像表示部１１０は、表示装置1１００に画像を表示するための表示制御手段として機能する。バス１０４０は、情報処理装置１０００内の各部を繋ぎ、各部の間でデータの送受信を行うことを可能にする。
【００６７】
（撮像装置）
次に、撮像装置２０００について説明する。撮像装置２０００は、例えば、Ｘ線透視装置のような、動画像を撮像することが可能であり、撮像された画像データは情報処理装置１０００に送信される。なお、画像データは複数をまとめて情報処理装置１０００に送信するようにしても良いし、撮像の都度、順次、画像データの送信を行うようにしてもよい。
【００６８】
（その他の実施例）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【図１】
0006021420.tif 000013
【図２】
0006021420.tif 000014
【図５】
0006021420.tif 000015
【図３】
0006021420.tif 000016
【図４】
0006021420.tif 000017