(57)【要約】【課題】Ｓ／Ｎ比が高く、かつ被検体内組織輪郭が受像部領域に因らず一様に描出されるサブトラクション像を得ることができるようにする。【解決手段】被写体に放射線をパルス状に照射して撮影した造影剤注入前の原マスク像を第１〜第ｎのフレームメモリ２３−ｉに記憶して、造影剤注入後のライブ像を撮影した際に、当該ライブ画像と同じ位置から撮影した原マスク像にフィルタ２８にてフィルタ処理を施して二次マスク像を生成し、生成した二次マスク像をライブ像から減算してサブトラクション像を取得するようにして、ライブ像から減算するマスク像におけるノイズを抑制できるようにするとともに、画像領域に因らず画像の輪郭を一様に描出することができるようにする。【選択図】  図２
2005012691.tif 000002
【特許請求の範囲】
【請求項１】
被写体に対して放射線をパルス状に照射する放射線照射手段と、
上記被写体を透過した放射線を検出して放射線強度に応じた電気信号を出力する放射線撮像手段と、
上記放射線撮像手段より出力される電気信号をディジタルデータに変換する変換手段と、
上記変換手段より出力される第１の画像のディジタルデータを記憶する記憶手段と、
上記第１の画像と同じ位置から撮影された上記第１の画像とは異なる第２の画像のディジタルデータから、上記第１の画像に所定の画像処理を施して得られたディジタルデータを減算する演算手段とを備えることを特徴とする放射線撮影装置。
【請求項２】
上記画像処理は、空間周波数領域におけるフィルタ処理であることを特徴とする請求項１に記載の放射線撮影装置。
【請求項３】
上記画像処理は、上記第１の画像の各画素における当該画素及び隣接画素の画素値の加重平均処理であることを特徴とする請求項１に記載の放射線撮影装置。
【請求項４】
上記第２の画像を解析処理する画像解析手段をさらに備え、上記演算手段は、上記画像解析手段による解析結果に応じて、上記第１の画像に施す画像処理のパラメータを変更することを特徴とする請求項１〜３の何れか１項に記載の放射線撮影装置。
【請求項５】
被写体に対して放射線をパルス状に照射する放射線照射手段と、
上記被写体を透過した放射線を検出し、検出した放射線強度に対応する電気信号を出力する放射線撮像手段と、
上記放射線撮像手段より出力される電気信号をディジタルデータに変換する変換手段と、
上記変換手段より出力される造影剤注入前の状態での原マスク像のディジタルデータを記憶する記憶手段と、
上記変換手段より造影剤注入後の状態でのライブ像のディジタルデータが得られた際に、上記ライブ像に対応する原マスク像に基づいて二次マスク像を生成して上記ライブ像から減算する演算手段とを備えることを特徴とする放射線撮影装置。
【請求項６】
上記二次マスク像は、上記原マスク像に空間周波数領域におけるフィルタ処理を施した画像であることを特徴とする請求項５に記載の放射線撮影装置。
【請求項７】
上記演算手段は、上記ライブ像における造影部を識別した結果に基づいて画像処理パラメータを変更して上記二次マスク像を生成することを特徴とする請求項６に記載の放射線撮影装置。
【請求項８】
上記画像処理パラメータは、上記空間周波数領域におけるフィルタ処理のカットオフ周波数であることを特徴とする請求項７に記載の放射線撮影装置。
【請求項９】
上記二次マスク像は、上記原マスク像にコンボリューション演算を施した画像であることを特徴とする請求項５に記載の放射線撮影装置。
【請求項１０】
上記演算手段は、上記ライブ像における造影部を識別した結果に基づいて画像処理パラメータを変更して上記二次マスク像を生成することを特徴とする請求項９に記載の放射線撮影装置。
【請求項１１】
上記画像処理パラメータは、上記コンボリューション演算におけるピクセルグループマスクのインテンシティ値であることを特徴とする請求項１０に記載の放射線撮影装置。
【請求項１２】
上記画像処理パラメータは、上記コンボリューション演算におけるピクセルグループマスクのマスクサイズであることを特徴とする請求項１０又は１１に記載の放射線撮影装置。
【請求項１３】
被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段と、
上記放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して得られたディジタルデータを、上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから減算する画像信号処理手段とを備えることを特徴とする放射線撮影装置。
【請求項１４】
被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成し、
上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算することを特徴とする放射線画像処理方法。
【請求項１５】
被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成するステップと、
上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算するステップとをコンピュータに実行させるためのプログラムを記録したことを特徴とするコンピュータ読み取り可能な記録媒体。
【請求項１６】
被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成するステップと、
上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算するステップとをコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、放射線撮影装置、放射線画像処理方法、記録媒体及びプログラムに関し、特に、撮影した放射線画像をディジタルデータで取得するディジタル放射線撮影装置に関し、例えばディジタルサブトラクションアンギオグラフィ（以下、「ＤＳＡ」と称す。）装置に用いて好適なものである。
【０００２】
【従来の技術】
回転ＤＳＡは、放射線発生部と受像部とを被検体の周囲に回転させ、造影剤注入前の透過放射線像（マスク像）と造影剤注入後の透過放射線像（ライブ像）とを所定の回転角度毎にそれぞれ撮影して、同一角度から撮影したライブ像とマスク像との間で減算を行うものである。従来のＤＳＡ装置においては、同一角度から撮影した１枚ずつのライブ像とマスク像との間で減算を行い、当該角度でのサブトラクション像を得ていた。
【０００３】
【特許文献１】
特許第２８７１０３７号公報
【０００４】
【発明が解決しようとする課題】
しかしながら、同一角度から撮影した１枚ずつのライブ像とマスク像との間で減算を行う従来のＤＳＡ装置においては、得られるサブトラクション像のＳ／Ｎ比が悪いという問題があった。回転ＤＳＡでは、放射線発生部と受像部とを被検体の周囲に回転させながら撮影を行うことや、被検体の総被曝線量の点から短時間（５ｍｓ程度）のパルス放射線により１つの画像を撮影するため、マスク像及びライブ像ともにＳ／Ｎ比が悪い画像になる。マスク像とライブ像との間での加減算におけるノイズ重畳は二乗和で効いてくるので、サブトラクション像のＳ／Ｎ比はさらに悪化する。
【０００５】
また、サブトラクション像は、マスク像とライブ像との間での減算を行うことから造影剤注入部（血管等）以外の部位については完全に消去されてしまい、造影剤注入部の被検体内における位置関係の把握が困難になってしまう。
【０００６】
これらの問題点に鑑み、上記特許文献１には、ある角度でのライブ像に対して、前後の角度でのマスク像を含む複数のマスク像を加算平均して得られたマスク像を減算することによりサブトラクション像を得る方法が開示されている。上記特許文献１には、前後角度でのマスク像を加算平均することにより、演算後のマスク像はＳ／Ｎ比が向上し、且つ前後角度での画像成分がマスク像に取り入れられるので減算処理後も造影剤注入部以外の部位についての画像が完全には消去されずに輪郭が残り、位置関係の把握が容易になると示されている。
【０００７】
しかしながら、上記特許文献１に開示されている方法では、異なる角度の画像情報により被検体画像の輪郭を描出している。ここで、被検体に対する放射線の照射角度に応じた被検体内組織輪郭部の受像部上への投影位置は、受像部の領域により異なる。すなわち、２つの異なる角度での画像の輪郭部のズレ量は、受像部中心ほど小さく、受像部中心から周辺部に至るに従って（受像部中心からの距離の増加にともなって）大きくなっていく。したがって、輪郭部のズレ量が小さい受像部中心ほど被検体内組織輪郭の描出が行われなくなる。
【０００８】
ここで、サブトラクション像にて注目する関心領域は、通常、画像の中心部、すなわち受像部中心に配置させるようにすることが多いので、上記特許文献１に開示されている方法では関心領域ほど位置関係の把握が困難になるという問題がある。さらに、観察対象になるべき造影部境界においても背景部位の輪郭描出が行われるので、造影部境界においてはフィルタ処理の影響による鮮鋭度の低下や輪郭のオーバーシュートなど不必要な効果が減算処理により加えられ、診断能を劣化させるおそれがある。
【０００９】
本発明は、このような問題を解決するために成されたものであり、Ｓ／Ｎ比が高く、かつ被検体内組織輪郭が受像部領域に因らず一様に描出されるサブトラクション像を得ることができるようにすることを目的とする。
【００１０】
【課題を解決するための手段】
本発明の放射線撮影装置は、被写体に対して放射線をパルス状に照射する放射線照射手段と、上記被写体を透過した放射線を検出して放射線強度に応じた電気信号を出力する放射線撮像手段と、上記放射線撮像手段より出力される電気信号をディジタルデータに変換する変換手段と、上記変換手段より出力される第１の画像のディジタルデータを記憶する記憶手段と、上記第１の画像と同じ位置から撮影された上記第１の画像とは異なる第２の画像のディジタルデータから、上記第１の画像に所定の画像処理を施して得られたディジタルデータを減算する演算手段とを備えることを特徴とする。
【００１１】
本発明の放射線撮影装置は、被写体に対して放射線をパルス状に照射する放射線照射手段と、上記被写体を透過した放射線を検出し、検出した放射線強度に対応する電気信号を出力する放射線撮像手段と、上記放射線撮像手段より出力される電気信号をディジタルデータに変換する変換手段と、上記変換手段より出力される造影剤注入前の状態での原マスク像のディジタルデータを記憶する記憶手段と、上記変換手段より造影剤注入後の状態でのライブ像のディジタルデータが得られた際に、上記ライブ像に対応する原マスク像に基づいて二次マスク像を生成して上記ライブ像から減算する演算手段とを備えることを特徴とする。
【００１２】
本発明の放射線撮影装置は、被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段と、上記放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して得られたディジタルデータを、上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから減算する画像信号処理手段とを備えることを特徴とする。
【００１３】
本発明の放射線画像処理方法は、被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成し、上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算することを特徴とする。
【００１４】
本発明のコンピュータ読み取り可能な記録媒体は、被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成するステップと、上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算するステップとをコンピュータに実行させるためのプログラムを記録したことを特徴とする。
【００１５】
本発明のプログラムは、被写体に対して放射線をパルス状に照射し、上記被写体を透過した放射線を検出して得られる画像のディジタルデータを出力する放射線画像取得手段より出力される第１の画像のディジタルデータに所定の画像処理を施して減算ディジタルデータを生成するステップと、上記第１の画像とは異なる上記放射線画像取得手段より出力される第２の画像のディジタルデータから上記減算ディジタルデータを減算するステップとをコンピュータに実行させることを特徴とする。
【００１６】
【発明の実施の形態】
以下、本発明の実施形態を図面に基づいて説明する。
（第１の実施形態）
図１は、本発明の第１の実施形態による放射線撮影装置の一構成例を示すブロック図である。
【００１７】
図１において、放射線管球１及び受像部３がアーム１１に保持されており、アーム１１の中央部にはベッドに横たわった被検者２が配される。
放射線管球１は、放射線制御部９からの信号に基づき高圧発生部８より供給される電源により、被検者２に対してパルス状に放射線を発する。放射線管球１より曝射される放射線は、被検者２を透過して受像部３に入射される。
【００１８】
受像部３は、入射された放射線を放射線強度に応じた電気信号に変換し、受像制御部４を介して画像信号処理部１２に出力する。受像部３は、放射線強度に応じて入射される放射線を電気信号に変換するデバイス（装置）であれば良く、例えばイメージインテンシファイアとＴＶカメラとの組み合わせ、あるいは蛍光体と二次元光センサアレイとの組み合わせなど、その種別は問わない。
【００１９】
アーム１１は、駆動部１０により駆動されて被検者２の周囲を回転する。駆動部１０によるアーム１１の回転は、制御部６により制御される。さらに、制御部６は、回転角度によって放射線の曝射タイミングを制御するように放射線制御部９を制御し、ライブ像（造影剤注入後の放射線像）の撮影（取得）時に造影剤を被検者２に注入するインジェクタ７を制御する。
【００２０】
これら各機能部の動作タイミング、例えば放射線の曝射レート、放射線管球１に与える管電圧、管電流、曝射時間や、アーム１１の回転速度、インジェクタ７の造影剤注入のタイミングや注入時間等のパラメータは制御部６に接続されたキーボードなどの入力部５より設定され、これら設定に基づいて制御部６は接続された各機能部の制御を行う。
【００２１】
図２は、上記図１に示した画像信号処理部１２の構成例を示すブロック図である。図２に示すように画像信号処理部１２は、Ａ／Ｄ（アナログ／ディジタル）変換器２１、スイッチ２２、第１〜第ｎのフレームメモリ２３−ｉ（ｉは添え字であり、ｉ＝１〜ｎ）、マルチプレクサ２７、フィルタ２８及び減算器２９を有する。
【００２２】
受像部３から出力されたアナログの電気信号は、受像制御部４を介してＡ／Ｄ変換器２１に入力され、ディジタルデータに変換される。ここでスイッチ２２は、マスク像（造影剤非注入時の放射線像）の撮影時において閉じられる（オン状態にされる）。
【００２３】
これにより、マスク像の撮影時には、アーム１１が被検者２の周囲を回転しながら放射線管球１より被検者２に対して放射線を曝射して、異なる角度からの被検者２の透過放射線像を取得し、取得した放射線像のディジタルデータを第１〜第ｎのフレームメモリ２３−ｉに順次記憶する。したがって、第１〜第ｎのフレームメモリ２３−ｉには、造影剤が注入されていない状態での複数の異なる角度からの原マスク像の画像データ（以下、「原マスク像データ」と称す。）がそれぞれ記憶される。
【００２４】
原マスク像の撮影が終了すると、ライブ像の撮影が行われる。このとき、スイッチ２２は開かれる（オフ状態にされる）とともに、制御部６はインジェクタ７を動作させて被検者２の体内に造影剤を注入する。注入された造影剤が被検者２内の観察対象部に流入するのに要する時間が経過した後、再びアーム１１を回転駆動するとともに、放射線管球１による被検者２への放射線曝射が行われる。これにより、造影剤が注入された状態での透過放射線像であるライブ像が取得される。
【００２５】
ここで、ある角度からのライブ像の撮影が行われたとき、マルチプレクサ２７は、第１〜第ｎのフレームメモリ２３−ｉに記憶されている原マスク像データの中から同一角度で撮影された原マスク像データを選択してフィルタ２８に出力する。フィルタ２８は、マルチプレクサ２７より入力された原マスク像データに所定の画像処理（フィルタ処理）を施すことにより二次マスク像データを生成する。
【００２６】
フィルタ２８による二次マスク像の生成アルゴリズムとしては、画像（原マスク像）を空間周波数領域に展開してフィルタ処理を施して、再び空間領域に戻す方法や、任意の一点の画素に対して単純に周囲画素を加重平均するコンボリューション演算によりフィルタ処理を行う方法などがある。フィルタ２８でのフィルタ処理による二次マスク像の生成の詳細については後述する。
【００２７】
撮影されたライブ像データから、フィルタ２８にて生成された二次マスク像データが減算器２９にて減算処理され、その演算結果（サブトラクション像）は記録装置３０に供給されて記録されたり、表示装置３１に供給されて表示されたりする。このライブ像を撮影してサブトラクション像を取得する処理を、アーム１１を被検者２の周囲に回転させながら任意の角度で行い、それぞれの角度でのサブトラクション像を取得する。
【００２８】
二次マスク像を生成する際のフィルタ２８でのフィルタ処理の一例について説明する。
原マスク像を空間周波数領域に展開してフィルタ処理を施して、再び空間領域に戻す方法について説明する。
【００２９】
空間周波数領域におけるフィルタ処理は一般的には以下の手順で行われる。
ステップ１：空間領域における画像データｆ（ｉ，ｊ）にフーリエ変換を施してフーリエ係数Ｆ（ｕ，ｖ）を求める。
ステップ２：フーリエ係数Ｆ（ｕ，ｖ）を、周波数領域フィルタ関数Ｈ（ｕ，ｖ）を用いて補正する。ここで、フィルタ補正係数をＧ（ｕ，ｖ）であらわすと、
【００３０】
【数１】
2005012691.tif 000003
【００３１】
になる。
ステップ３：フィルタ補正係数Ｇ（ｕ，ｖ）にフーリエ逆変換を施してｆ（ｉ，ｊ）に対する処理画像ｇ（ｉ，ｊ）を得る。
ここで、ｉ，ｊ，ｕ，ｖは離散的変数を示し、大きさＮ×Ｎの画像データに対して、
【００３２】
【数２】
2005012691.tif 000004
【００３３】
と与えられる。ｕ，ｖに関するフーリエ係数Ｆ（ｕ，ｖ）の周期性からｕ，ｖを以下のように定める。
【００３４】
【数３】
2005012691.tif 000005
【００３５】
また、変数ｕ，ｖと周波数νには以下に示す（式１）の関係がある。
【００３６】
【数４】
2005012691.tif 000006
【００３７】
ここで、ｄは画像のピクセル幅を示す。
ナイキスト周波数をνｎとすると、フィルタ関数がν≦νｎ＝１／（２ｄ）を満足する周波数に対して定義されることから、大きさＮ×Ｎの周波数領域フィルタ関数Ｈ（ｕ，ｖ）は次のように定義される。
【００３８】
【数５】
2005012691.tif 000007
【００３９】
例えば、最終的に得られるサブトラクション像に背景部分の輪郭を描出させたい場合には、減算処理を行ってもライブ像における高周波成分が残存するようにマスク像における高周波成分を除去する低域通過形フィルタを使用する。この低域通過形フィルタとして、バターワースフィルタを使用する場合には、上記周波数領域フィルタ関数Ｈ（ｕ，ｖ）は以下の（式２）のように与えられる。
【００４０】
【数６】
2005012691.tif 000008
【００４１】
ここで、ｎはフィルタの次数を表し、ωｃはカットオフ周波数νｃに対応するもので、ωｃ＝νｃ／Δνである。
上記（式２）に示したフィルタ関数Ｈ（ｕ，ｖ）を用いて、上述したステップ１〜３を実行することにより、フィルタ２８により画像の高周波成分が除去された二次マスク像を生成することができる。
【００４２】
また、最終的に得られるサブトラクション像に背景部分の輪郭ではなく、その全体像をうっすらと（薄く）描出させたい場合には、減算処理によりライブ像における低周波成分が残存するようにマスク像の低周波成分を除去する高域通過形フィルタを使用すれば良い。このように、サブトラクション像に描出させたい背景の種類、あるいは原マスク像のノイズ除去を考慮して、これら諸特性を満足する帯域通過形や、帯域制限形のフィルタ、複数のフィルタから構成されるハイブリッドフィルタを用いても良い。
【００４３】
また、フィルタ２８でのフィルタ処理に、コンボリューション演算による処理を用いた場合でも、上述した周波数領域におけるフィルタ処理と同様の効果を得ることができる。最終的に得られるサブトラクション像に背景の輪郭部を描出させたい場合には、原マスク像に対してコンボリューション演算による平滑化処理を施す。平滑化処理のもっとも一般的な方法は、各画素の周囲に３×３、５×５など所定の大きさのピクセルグループマスクを設定し、設定したピクセルグループマスク内の画素について画素値の平均を求め、これを平滑値にする方法である。
【００４４】
以上、説明したように第１の実施形態によれば、撮影により得られた造影剤注入前の原マスク像にフィルタ処理を施して二次マスク像を生成し、造影剤注入後のライブ像から減算するマスク像として生成した二次マスク像を使用してサブトラクション像を取得する。これにより、ライブ像から減算するマスク像におけるノイズを抑制することができ、得られるサブトラクション像のＳ／Ｎ比の向上を図ることができる。
【００４５】
また、二次マスク像を用いてライブ像に対する減算処理を行うことにより、サブトラクション像は、造影剤流入部以外の部位（例えば、背景部位）も減算処理により完全に消失してしまうことがなく、サブトラクション像上にその情報（輪郭あるいは全体像）を残存させることができる。したがって、本来の観察対象である造影部位とその他の背景部位とが同一画像上に適度な割合で描出され、残った背景部位の画像を参照して造影部位の位置を容易に把握することができる。
【００４６】
さらに、一枚のライブ像に対して、当該ライブ像に対応する一枚の原マスク像から二次マスク像を生成することにより、回転ＤＳＡにおける従来のような前後角度からのマスク像、つまり複数のマスク像による平均化処理という限定にとらわれず、通常のＤＳＡにおいても適用することができる。
【００４７】
上記特許文献１に開示された方法で生じる受像部中心と周辺部との輪郭部のズレ量の相違に起因する背景画像の描出ムラにより画像中心部に位置する関心領域ほど背景画像の描出が少なく位置関係の把握が困難になるという問題を解消し、受像部のどの領域においても背景の描出が一様に可能となり、位置関係を容易に把握することができる。
【００４８】
（第２の実施形態）
次に、第２の実施形態について説明する
本発明の第２の実施形態による放射線撮影装置の構成については、上記図１に示した第１の実施形態による放射線撮影装置と同じ構成であるので、説明は省略する。
【００４９】
図３は、第２の実施形態における画像信号処理部１２’の構成例を示すブロック図である。なお、この図３において、図２に示したブロック等と同一の機能を有するブロック等には同一の符号を付し、重複する説明は省略する。
【００５０】
図３に示すように画像信号処理部１２’は、上記図２に示した画像信号処理部１２に、撮影により得られたライブ像の画像データを一時的に記憶するフレームバッファ４８と、当該ライブ像における造影部の解析を行う画像解析部４９とをさらに設ける。そして、フィルタ５０は、画像解析部４９での解析結果に応じて画像処理パラメータを変更して原マスク像データにフィルタ処理を施し、二次マスク像データを生成する。
【００５１】
原マスク像を順次撮影し、第１〜第ｎのフレームメモリ２３−ｉに原マスク像データを記憶する動作は上述した第１の実施形態と同様である。
原マスク像の撮影が終了すると、上述した第１の実施形態と同様にしてライブ像の撮影が行われる。第２の実施形態においては、撮影により得られたライブ像データはフレームバッファ４８に一時的に記憶される。
【００５２】
フレームバッファ４８に記憶されたライブ像データは、画像解析部４９に供給されて造影部の解析が行われる。画像中の造影部は造影剤により放射線透過量が他の部分と比較して低くなるために、ある部位に対して撮影した画像の全画素値のヒストグラムは、図４（Ａ）に示すように低画素値側に急峻な凸部を有する。
この凸部は造影部に相当するので、この凸部の前後に閾値を定めて範囲ＶＡを設定し、ライブ像を構成する各画素の画素値が設定した範囲ＶＡ内であるか否かをそれぞれ検証する。このようにしてライブ像の各画素を造影部とそれ以外の部分とに振り分けることにより、結果としてライブ像の２値化が行われる。
【００５３】
さらに別の方法として、フレームバッファ４８から入力されるライブ像と当該ライブ像に対応する原マスク像との減算処理を画像分析部４９にて行い、得られたサブトラクション像に対して造影部であるか否かの検証を行う方法も考えられる。図４（Ｂ）に示すようにライブ像の全画素値のヒストグラムにおいて、急峻な凸部以外は背景画像による部分である。したがって、ライブ像から原マスク像を減算することにより背景画像部分をキャンセルして、造影部の抽出をより明確に行い造影部の抽出精度を向上させることができる。
【００５４】
ここで、ある角度からのライブ像の撮影が行われたとき、マルチプレクサ２７は、同一角度で取得された原マスク像データを第１〜第ｎのフレームメモリ２３−ｉより選択しフィルタ５０に出力する。フィルタ５０は、マルチプレクサ２７より入力された原マスク像に基づいて、画像処理を施した二次マスク像を生成する。フィルタ５０による二次マスク像の生成アルゴリズムは上述した第１の実施形態と同様である。
【００５５】
しかしながら、造影剤が流入するのは血管等のように障壁を有する器官であり、この障壁が観察部位とそれ以外の背景部を分けていることになる。原マスク像から二次マスク像を生成する際、原マスク像領域全体に対して同様のフィルタ処理を施すと、この障壁に対してもフィルタ処理が施されることになる。したがって、二次マスク像をライブ像より減算処理して得られるサブトラクション像の造影部境界においては、フィルタ処理の影響が及ぼす鮮鋭度の低下や輪郭のオーバーシュートなどが発生する可能性がある。
【００５６】
そこで、第２の実施形態においては、造影部に対するフィルタ処理の影響を緩和するために、画像解析部４９で得られた造影部の情報に基づいてライブ像中の造影部に対応する領域については画像処理パラメータを変更する、あるいはその領域については画像処理そのものを行わないようにして、原マスク像にフィルタ処理を施して二次マスク像を生成する。
【００５７】
フィルタ５０でのフィルタ処理に空間周波数領域におけるフィルタ処理を用いた場合には、原マスク像を構成する全画素が入力になり、得られる結果はフィルタ処理が施された二次マスク像の全画素である。そこで、フィルタ５０は、画像解析部４９から供給される造影部の領域情報に基づいて、造影部に対応する領域の画素についてはフィルタ処理により得られた結果を用いずに、原マスク像の画素値を用いて二次マスク像を生成する。
【００５８】
あるいはフィルタ５０は、非造影部に対するカットオフ周波数よりも高いカットオフ周波数でのフィルタ処理を別途行い、造影部に対応する画素については、高いカットオフ周波数でのフィルタ処理により得られた画像データを採用する方法でも良い。
【００５９】
さらに別の方法として、フィルタ５０にてフィルタ処理を行う画像サイズを全画像領域、つまり１つの画像領域ではなく、複数の領域に分割して個別にフィルタ処理を行う。そして、画像解析部４９から供給される造影部の領域情報に基づいて、分割したフィルタ処理単位の領域に造影部が含まれると判断された場合には、フィルタ５０のカットオフ周波数を高周波側にシフトさせる方法もある。
【００６０】
また、フィルタ５０でのフィルタ処理に、コンボリューション演算による処理を用いた場合には、原マスク像を構成する各画素と当該画素を中心とする周囲画素とで構成されるピクセルグループマスクにおける重み係数を造影部の領域情報に基づいて設定する。
【００６１】
図５は、ピクセルグループマスクにおける重み係数の設定例を示す図である。
図５に示す例では、造影部に対応しない領域の場合には、図５（Ａ）に示すように全画素の加重平均時の重み係数であるインテンシティ値は一律に“１”とし、単なる加算平均処理が行われているものとする。
【００６２】
例えば、ピクセルグループマスク内の周囲画素のある画素が、画像解析部４９から供給される造影部の領域情報に合致したとする。このとき、フィルタ５０は、図５（Ｂ）に示すように造影部に対応する画素５１のインテンシティ値を他画素に対して相対的に低くする。あるいは、図５（Ｃ）に示すように造影部に対応する画素５２のインテンシティ値を“０”にして、加重平均処理の対象から除外し造影部に対応する画素の影響を排除することができる。
【００６３】
図５に示したようにピクセルグループマスクにおける重み係数（インテンシティ値）を変更するようにした場合の中心画素の平滑値ｓ（ｉ，ｊ）は、ピクセルグループマスクの各画素値ｆ（ｉ，ｊ）とそれぞれの画素に対応するインテンシティ値ｇ（ｋ，ｌ）とから、
【００６４】
【数７】
2005012691.tif 000009
【００６５】
となる。
なお、ピクセルグループマスクにおける中心画素が造影部に対応する画素の場合には、周囲画素のインテンシティ値を“０”にする、すなわち原マスク像の画素値をそのまま用いれば良い。
【００６６】
また、別の方法としては、ピクセルグループマスクのマスクサイズを変更する方法がある。
図６は、ピクセルグループマスクのマスクサイズの変更例を示す図である。
造影部に対応しない領域の場合には、例えば図６（Ａ）に示すようにピクセルグループマスクサイズは７×７とする。例えば、７×７のピクセルグループマスク内に造影部に対応する領域が含まれる場合には、図６（Ｂ）に示すようにマスクサイズを５×５に縮小して造影部に対応する領域の影響を排除する。
【００６７】
さらに、５×５にマスクサイズを縮小しても造影部に対応する領域を排除できない場合には、図６（Ｃ）に示すようにマスクサイズを３×３に縮小し、最終的にはマスクサイズ１×１、つまり自身の画素のみになるまでピクセルグループマスクのマスクサイズを変更する。ここで、マスクサイズが１×１になる場合には、当該画素（中心画素）は、造影部に対応する領域に隣接する画素である。なお、中心画素が造影部に対応する画素の場合には、原マスク像の画素値を用いれば良い。
【００６８】
以上説明したような空間周波数領域におけるフィルタ処理や、コンボリューション演算によるフィルタ処理等により、フィルタ５０は、造影部に対応する領域とそれ以外の背景領域とについて、それぞれの領域の影響を抑制した二次マスク像を生成することができる。
【００６９】
上述のようにしてフィルタ５０にて生成した二次マスク像データが、フレームバッファ４８から読み出されたライブ像データから減算器２９にて減算処理され、その演算結果（サブトラクション像）は記録装置３０に供給されて記録されたり、表示装置３１に供給されて表示されたりする。このライブ像を撮影してサブトラクション像を取得する処理を、アーム１１を被検者２の周囲に回転させながら任意の角度で行い、それぞれの角度でのサブトラクション像を取得する。
【００７０】
以上、説明したように第２の実施形態によれば、第１の実施形態により得られる効果に加え、撮影したライブ像から抽出した造影部の領域情報に応じて、原マスク像にフィルタ処理を施す際のパラメータを変更して二次マスク像を生成することにより、原マスク像から二次マスク像を生成する際の造影部境界におけるフィルタ処理の影響を緩和することができる。これにより、サブトラクション像における造影部境界での鮮鋭度の低下、オーバーシュート等を抑制し、造影部周辺を明確に描出することができ、診断能の劣化を防止することができる。
【００７１】
なお、上述した第１及び第２の実施形態においては、回転ＤＳＡに適用した場合を一例として説明したが、本発明はこれに限定されるものではなく、例えば、単一画像に対する通常のＤＳＡにおいても適用可能であり、第１〜第ｎのフレームメモリ２３−ｉの容量及び個数が変わり、マルチプレクサ２７が不要であるなど、必須要件の増加はない。
【００７２】
（本発明の他の実施形態）
上述した実施形態の機能を実現するべく各種のデバイスを動作させるように、該各種デバイスと接続された装置あるいはシステム内のコンピュータに対し、上記実施形態の機能を実現するためのソフトウェアのプログラムコードを供給し、そのシステムあるいは装置のコンピュータ（ＣＰＵあるいはＭＰＵ）に格納されたプログラムに従って上記各種デバイスを動作させることによって実施したものも、本発明の範疇に含まれる。
【００７３】
また、この場合、上記ソフトウェアのプログラムコード自体が上述した実施形態の機能を実現することになり、そのプログラムコード自体は本発明を構成する。また、そのプログラムコードをコンピュータに供給するための手段、例えばかかるプログラムコードを格納した記録媒体は本発明を構成する。かかるプログラムコードを記憶する記録媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００７４】
また、コンピュータが供給されたプログラムコードを実行することにより、上述の実施形態の機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）あるいは他のアプリケーションソフト等と共同して上述の実施形態の機能が実現される場合にもかかるプログラムコードは本発明の実施形態に含まれることは言うまでもない。
【００７５】
さらに、供給されたプログラムコードがコンピュータの機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに格納された後、そのプログラムコードの指示に基づいてその機能拡張ボードや機能拡張ユニットに備わるＣＰＵ等が実際の処理の一部又は全部を行い、その処理によって上述した実施形態の機能が実現される場合にも本発明に含まれることは言うまでもない。
【００７６】
例えば、第１及び第２の実施形態に示した放射線撮影装置は、図７に示すようなコンピュータ機能７００を有し、そのＣＰＵ７０１により第１及び第２の実施形態での動作が実施される。
【００７７】
コンピュータ機能７００は、上記図７に示すように、ＣＰＵ７０１と、ＲＯＭ７０２と、ＲＡＭ７０３と、キーボード（ＫＢ）７０９のキーボードコントローラ（ＫＢＣ）７０５と、表示部としてのＣＲＴディスプレイ（ＣＲＴ）７１０のＣＲＴコントローラ（ＣＲＴＣ）７０６と、ハードディスク（ＨＤ）７１１及びフレキシブルディスク（ＦＤ）７１２のディスクコントローラ（ＤＫＣ）７０７と、ネットワークインタフェースカード（ＮＩＣ）７０８とが、システムバス７０４を介して互いに通信可能に接続された構成としている。
【００７８】
ＣＰＵ７０１は、ＲＯＭ７０２あるいはＨＤ７１１に記憶されたソフトウェア、あるいはＦＤ７１２より供給されるソフトウェアを実行することで、システムバス７０４に接続された各構成部を総括的に制御する。
すなわち、ＣＰＵ７０１は、上述したような動作を行うための処理プログラムを、ＲＯＭ７０２、あるいはＨＤ７１１、あるいはＦＤ７１２から読み出して実行することで、第１及び第２の実施形態での動作を実現するための制御を行う。
【００７９】
ＲＡＭ７０３は、ＣＰＵ７０１の主メモリあるいはワークエリア等として機能する。
ＫＢＣ７０５は、ＫＢ７０９や図示していないポインティングデバイス等からの指示入力を制御する。
ＣＲＴＣ７０６は、ＣＲＴ７１０の表示を制御する。
ＤＫＣ７０７は、ブートプログラム、種々のアプリケーション、ユーザファイル、ネットワーク管理プログラム、および第１及び第２の実施形態における上記処理プログラム等を記憶するＨＤ７１１及びＦＤ７１２とのアクセスを制御する。
ＮＩＣ７０８はネットワーク７１３上の他の装置と双方向にデータをやりとりする。
【００８０】
なお、上記実施形態は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００８１】
【発明の効果】
以上説明したように、本発明によれば、被写体にパルス状に放射線を照射して撮影した第１の画像のディジタルデータを記憶し、上記第１の画像と同じ位置から撮影した上記第１の画像とは異なる第２の画像のディジタルデータから上記第１の画像に所定の画像処理を施して得られたディジタルデータを減算して、第１の画像と第２の画像とに基づいた差分画像を得る。これにより、第１の画像におけるノイズを抑制して第２の画像に対する減算処理を行うことができ、差分画像のＳ／Ｎ比を向上させることができる。また、同じ位置から撮影した第１及び第２の画像を用いるので、画像領域に関わらず画像の輪郭が一様に描出された差分画像を得ることができる。
【図面の簡単な説明】
【図１】第１の実施形態による放射線撮影装置の一構成例を示すブロック図である。
【図２】第１の実施形態における画像信号処理部の構成例を示すブロック図である。
【図３】第２の実施形態における画像信号処理部の構成例を示すブロック図である。
【図４】第２の実施形態における画像解析部での造影部抽出方法を説明するための図である。
【図５】第２の実施形態におけるコンボリューション演算に係る重み係数の設定例を示す図である。
【図６】第２の実施形態におけるコンボリューション演算に係るマスクサイズの変更例を示す図である。
【図７】第１及び第２の実施形態による放射線撮影装置を実現可能なコンピュータ機能を示すブロック図である。
【符号の説明】
１  放射線管球
２  被検者
３  受像部
４  受像制御部
５  入力部
６  制御部
７  インジェクタ
８  高圧発生部
９  放射線制御部
１０  駆動部
１２、１２’  画像信号処理部
２１  Ａ／Ｄ変換器
２２  スイッチ
２３  フレームメモリ
２７  マルチプレクサ
２８、５０  フィルタ
２９  減算器
４８  フレームバッファ
４９  画像解析部
【図１】
2005012691.tif 000010
【図２】
2005012691.tif 000011
【図３】
2005012691.tif 000012
【図４】
2005012691.tif 000013
【図５】
2005012691.tif 000014
【図６】
2005012691.tif 000015
【図７】
2005012691.tif 000016