(57)【要約】【課題】高精度で安定した領域抽出を可能とする放射線画像領域抽出装置、放射線画像領域抽出プログラムおよび放射線撮影装置を提供することを目的とする。【解決手段】画像領域抽出部１１は、放射線画像を被検体Ｍ領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する低周波画像生成部１５と、低周波画像の輝度変化を検出して輝度変化画像を生成する輝度変化画像生成部１７と、輝度変化画像および低周波画像の画素値に基づいて関心領域の境界に位置する境界画素を選択する境界選択部１９とを備えている。低周波画像を生成することで、ノイズ成分や被検体Ｍ内部の細かい構造成分の影響を低減し、被検体Ｍの輪郭形状の認識を容易にしている。低周波画像から生成した輝度変化画像と、低周波画像とを組み合わせて境界画素を選択するので、高精度で安定した領域抽出を行うことができる。【選択図】図１
2011255033.tif 000002
【特許請求の範囲】
【請求項１】
撮影された放射線画像の中から被検体の関心領域を抽出する放射線画像領域抽出装置において、
前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する低周波画像生成部と、
前記低周波画像の輝度変化を検出して輝度変化画像を生成する輝度変化画像生成部と、
前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する境界選択部と、
を備えていることを特徴とする放射線画像領域抽出装置。
【請求項２】
請求項１に記載の放射線画像領域抽出装置において、
前記低周波画像生成部は、前記放射線画像を縮小することにより前記低周波画像を生成することを特徴とする放射線画像領域抽出装置。
【請求項３】
請求項２に記載の放射線画像領域抽出装置において、
放射線画像の縮小は、１／ｎの縮小率の場合、ｎ×ｎ画素の平均画素値を１画素の画素値とすることを特徴とする放射線画像領域抽出装置。
【請求項４】
請求項２に記載の放射線画像領域抽出装置において、
放射線画像の縮小は、１／ｎの縮小率の場合、ｎ×ｎ画素のうちの間引いた画素の平均画素値を１画素の画素値とすることを特徴とする放射線画像領域抽出装置。
【請求項５】
請求項２から４のいずれかに記載の放射線画像領域抽出装置において、
さらに、前記境界選択部で選択された境界画素を前記放射線画像のサイズに拡大する処理を行うとともに、前記放射線画像の中から前記関心領域を抽出する処理を行う領域抽出処理部を備えていることを特徴とする放射線画像領域抽出装置。
【請求項６】
請求項１から５のいずれかに記載の放射線画像領域抽出装置において、
境界選択部は、輝度変化画像の中から輝度変化の最も大きい画素である初期画素を選択し、選択した初期画素に隣接する全方向の画素から輝度変化の最も大きい画素を次の境界画素として選択し、その後、選択した画素の進行方向に隣接する画素から輝度変化の最も大きい画素を次の境界画素として選択することを繰り返して前記関心領域の境界画素を選択することを特徴とする放射線画像領域抽出装置。
【請求項７】
請求項６に記載の放射線画像領域抽出装置において、
境界選択部は、輝度変化画像の新たに選択した画素の画素値が予め設定した所定値以下の場合、輝度変化画像の選択前の画素と同位置における低周波画像の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択することを特徴とする放射線画像領域抽出装置。
【請求項８】
請求項６または７に記載の放射線画像領域抽出装置において、
境界選択部は、輝度変化画像の新たに選択した画素および選択前の画素と同位置における低周波画像の各画素の画素値を比較して予め設定した所定値よりも大きい場合、輝度変化画像の選択前の画素と同位置における低周波画像の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択することを特徴とする放射線画像領域抽出装置。
【請求項９】
撮影された放射線画像の中から被検体の関心領域を抽出する処理をコンピュータに実行させるための放射線画像領域抽出プログラムにおいて、
前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する工程と、
前記低周波数画像の輝度変化を検出して輝度変化画像を生成する工程と、
前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する工程と、を備え、
これらの工程処理をコンピュータに実行させることを特徴とする放射線画像領域抽出プログラム。
【請求項１０】
被検体を撮影する放射線撮影装置において、
放射線を検出する放射線検出器と、
前記放射線検出器から出力された放射線画像の中から被検体の関心領域を抽出する画像領域抽出部と、を備え、
前記画像領域抽出部は、前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する低周波画像生成部と、前記低周波画像の輝度変化を検出して輝度変化画像を生成する輝度変化画像生成部と、前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する境界選択部とを備えていることを特徴とする放射線撮影装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、撮影された放射線画像の中から被検体の関心領域を抽出する放射線画像領域抽出装置、放射線画像領域抽出プログラムおよび放射線撮影装置に関する。
【背景技術】
【０００２】
放射線撮影装置では、診断に適した画像を得るために、放射線（例えば、Ｘ線またはγ線等）を検出して撮影された放射線画像の中から被検体の関心領域を抽出し、抽出した被検体の関心領域に対し、種々の画像処理が行われる。被検体の胸部の撮影において、放射線画像の中で特に診断に必要とされるのは肺野領域である。しかしながら、放射線画像の中には、診断に必要な肺野領域とともに、診断に不要な腹部や腕、あるいは直接放射線が入射する領域なども含まれてしまう。これらの領域は、被検体の肺野領域である関心領域と比べて画像濃度やコントラストが大きく異なるので、より高性能な画像処理を行うために、関心領域以外を画像処理の領域から除外することが求められる（例えば、特許文献１参照）。
【０００３】
特許文献１では、まず、被検体の肺野領域を抽出するために、しきい値に基づき画像を２値化する処理を行い、２値化された画像にラベリング処理を行う。そして、ラベリング処理された複数の領域のうち、所定の面積以下の領域は、肺野領域外で高濃度を示す領域として、また、画像の上下左右の端に接する領域は、Ｘ線の「す抜け領域」として除去される。一方、除去されなかった領域を肺野領域として抽出している。また、しきい値は、例えば、す抜け領域を除去した被検体の画像領域内から最大画素（濃度）値と、この最大画素値の座標を通るプロファイルの凹部から最小画素値を算出し、最大画素値と最小画素値から推定している。
【０００４】
また、この種の他の装置として、例えば、特許文献２が開示されている。特許文献２では、まず、マンモグラフィ画像から診断情報の無い背景空間（す抜け領域）を除去するためにしきい値に基づき２進画像を生成する。そして、皮膚ライン画素を決定し、領域を抽出している。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４１４３１４９号公報
【特許文献２】特開２００２−３６９０７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、このような従来装置では、す抜け領域を除去するためのしきい値、誤抽出された所定の面積以下の領域を除外するための面積値、または、被検体の関心領域を抽出するためのしきい値などを設定しているが、例えば、撮影対象や撮影条件が変化した場合に、その都度、最適な値に設定し直さなければ、精度の高い領域抽出を行うことができないという問題がある。
【０００７】
また、設定するしきい値やノイズなどの影響によって、抽出される領域が変化し、精度にばらつきが生じてしまう問題がある。さらに、画像全体に対して２値化処理とラベリング処理で領域検出を行っているので、処理時間がかかってしまい、かつノイズなどの影響も受けやすいという問題がある。
【０００８】
本発明は、このような事情に鑑みてなされたものであって、高精度で安定した領域抽出を可能とする放射線画像領域抽出装置、放射線画像領域抽出プログラムおよび放射線撮影装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、このような目的を達成するために、次のような構成をとる。すなわち、本発明に係る放射線画像領域抽出装置は、撮影された放射線画像の中から被検体の関心領域を抽出する放射線画像領域抽出装置において、前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する低周波画像生成部と、前記低周波画像の輝度変化を検出して輝度変化画像を生成する輝度変化画像生成部と、前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する境界選択部と、を備えていることを特徴とするものである。
【００１０】
本発明に係る放射線画像領域抽出装置によれば、低周波画像生成部は、放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成し、輝度変化画像生成部は、低周波画像の輝度変化を検出して輝度変化画像を生成する。すなわち、低周波画像を生成することで、ノイズ成分や被検体内部の細かい構造成分の影響を低減し、被検体の輪郭形状の認識を容易にしている。また、低周波画像から生成された輝度変化画像は、低周波画像と同様の効果を有している。境界選択部は、そのような輝度変化画像および低周波画像を組み合わせて被検体の関心領域の境界に位置する境界画素を選択しているので、高精度で安定した領域抽出を行うことができる。
【００１１】
また、本発明に係る放射線画像領域抽出装置において、前記低周波画像生成部は、前記放射線画像を縮小することにより前記低周波画像を生成することが好ましい。通常のフィルタに比べて、より低周波な低周波画像を生成することができるので、ノイズなどの外乱の影響が少なく、高精度で安定した領域抽出を行うことができる。また、低周波画像を用いて輝度変化画像を生成したり、境界画素を選択したりしているので、比較的演算するデータ量が少なくでき、被検体の関心領域の抽出を高速に行うことができる。
【００１２】
また、本発明に係る放射線画像領域抽出装置において、放射線画像の縮小の一例は、１／ｎの縮小率の場合、ｎ×ｎ画素の平均画素値を１画素の画素値とすることである。また、放射線画像の縮小の一例は、１／ｎの縮小率の場合、ｎ×ｎ画素のうちの間引いた画素の平均画素値を１画素の画素値とすることである。これらによって放射線画像を縮小することにより低周波画像を生成することができる。また、間引きする場合は、縮小する際の演算するデータ量を少なくすることができる。
【００１３】
また、本発明に係る放射線画像領域抽出装置において、さらに、前記境界選択部で選択された境界画素を前記放射線画像のサイズに拡大する処理を行うとともに、前記放射線画像の中から前記関心領域を抽出する処理を行う領域抽出処理部を備えていることが好ましい。放射線画像を縮小して低周波画像を生成した場合、元の放射線画像のサイズに戻す拡大する処理を行うことにより、放射線画像の中から被検体の関心領域を抽出することができる。
【００１４】
また、本発明に係る放射線画像領域抽出装置において、境界選択部の一例は、輝度変化画像の中から輝度変化の最も大きい画素である初期画素を選択し、選択した初期画素に隣接する全方向の画素から輝度変化の最も大きい画素を次の境界画素として選択し、その後、選択した画素の進行方向に隣接する画素から輝度変化の最も大きい画素を次の境界画素として選択することを繰り返して前記関心領域の境界画素を選択するものである。低周波画像から生成された輝度変化画像の中から、選択した画素に隣接する輝度変化の最も大きい画素を境界画素として選択しているので、高精度で安定した領域抽出することができる。
【００１５】
また、本発明に係る放射線画像領域抽出装置において、境界選択部は、輝度変化画像の新たに選択した画素の画素値が予め設定した所定値以下の場合、輝度変化画像の選択前の画素と同位置における低周波画像の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択することが好ましい。輝度変化画像の中において、境界選択部で新たに選択した画素の画素値が予め設定した所定値以下の場合は、境界画素として信頼性が低いと判定される。すなわち、輝度変化が小さいほど境界画素として信頼性が低くなる。この場合、低周波画像から、選択前の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択する。最も値が近い画素を選択することは、同じ領域内で境界画素である可能性が高いことを示している。したがって、新たに選択した画素を判定して信頼性が低いものを振り分けし、信頼性が低いと判定された場合でも低周波画像から境界画素を選択しているので、より高精度で安定した領域抽出を行うことができる。
【００１６】
また、本発明に係る放射線画像領域抽出装置において、境界選択部は、輝度変化画像の新たに選択した画素および選択前の画素と同位置における低周波画像の各画素の画素値を比較して予め設定した所定値よりも大きい場合、輝度変化画像の選択前の画素と同位置における低周波画像の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択することが好ましい。低周波画像において、輝度変化画像で新たに選択した画素および選択前の画素と同位置の各画素の画素値を比較して予め設定した所定値よりも大きい場合は、新たに選択した画素が境界画素として信頼性が低いと判定される。すなわち、輝度変化の大きい画素であっても（新たに選択した画素の画素値と選択前の画素の画素値とが同じような大きさであっても）、低周波画像において画素値が大きく異なることがある。この場合も低周波画像から、選択前の画素の画素値と、その画素に隣接する画素のうち最も値が近い画素を次の境界画素として選択する。最も値が近い画素を選択することは、同じ領域内で境界画素である可能性が高いことを示している。したがって、新たに選択した画素を判定して信頼性が低いものを振り分けし、信頼性が低いと判定された場合でも低周波画像から境界画素を選択しているので、より高精度で安定した領域抽出を行うことができる。
【００１７】
また、本発明に係る放射線画像領域抽出プログラムは、撮影された放射線画像の中から被検体の関心領域を抽出する処理をコンピュータに実行させるための放射線画像領域抽出プログラムにおいて、前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する工程と、前記低周波数画像の輝度変化を検出して輝度変化画像を生成する工程と、前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する工程と、を備え、これらの工程処理をコンピュータに実行させることを特徴とするものである。
【００１８】
本発明の放射線画像領域抽出プログラムによれば、撮影された放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成し、この低周波画像の輝度変化を検出して輝度変化画像を生成する。生成された輝度変化画像および低周波画像の画素値に基づいて、画像濃度やコントラストが大きく異なる被検体の関心領域の境界に位置する境界画素を選択する。放射線画像から低周波画像を生成することで、ノイズ成分や被検体内部の細かい構造成分の影響を低減し、被検体の輪郭形状の認識を容易にしている。また、低周波画像から輝度変化画像を生成することで、輝度変化画像は、低周波画像と同様の効果を有している。境界選択部は、そのような輝度変化画像および低周波画像を組み合わせて、被検体の関心領域の境界に位置する境界画素を選択しているので、高精度で安定した領域抽出を行うことができる。
【００１９】
本発明に係る放射線撮影装置は、被検体を撮影する放射線撮影装置において、放射線を検出する放射線検出器と、前記放射線検出器から出力された放射線画像の中から被検体の関心領域を抽出する画像領域抽出部と、を備え、前記画像領域抽出部は、前記放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成する低周波画像生成部と、前記低周波画像の輝度変化を検出して輝度変化画像を生成する輝度変化画像生成部と、前記輝度変化画像および前記低周波画像の画素値に基づいて前記関心領域の境界に位置する境界画素を選択する境界選択部とを備えていることを特徴とするものである。
【００２０】
本発明の放射線撮影装置によれば、放射線検出器は、放射線を検出して放射線画像を出力し、放射線画像領域抽出装置は、放射線検出器から出力された放射線画像の中から被検体の関心領域を抽出する。放射線画像領域抽出装置において、低周波画像生成部は、放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって、放射線画像の空間分解能を劣化させて低周波画像を生成し、輝度変化画像生成部は、低周波画像の輝度変化を検出して輝度変化画像を生成する。生成された輝度変化画像および低周波画像の画素値に基づいて、画像濃度やコントラストが大きく異なる被検体の関心領域の境界に位置する境界画素を選択する。放射線画像から低周波画像を生成することで、ノイズ成分や被検体内部の細かい構造成分の影響を低減し、被検体の輪郭形状の認識を容易にしている。また、低周波画像から生成された輝度変化画像は、低周波画像と同様の効果を有している。境界選択部は、そのような輝度変化画像および低周波画像を組み合わせて、被検体の関心領域の境界に位置する境界画素を選択しているので、高精度で安定した領域抽出を行うことができる。
【発明の効果】
【００２１】
本発明によれば、低周波画像生成部は、放射線画像を被検体領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成し、輝度変化画像生成部は、低周波画像の輝度変化を検出して輝度変化画像を生成する。すなわち、低周波画像を生成することで、ノイズ成分や被検体内部の細かい構造成分の影響を低減し、被検体の輪郭形状の認識を容易にしている。また、低周波画像から生成された輝度変化画像は、低周波画像と同様の効果を有している。境界選択部は、そのような輝度変化画像および低周波画像を組み合わせて、被検体の関心領域の境界に位置する境界画素を選択しているので、高精度で安定した領域抽出を行うことができる。
【図面の簡単な説明】
【００２２】
【図１】実施例に係る放射線撮影装置の概略構成を示す図である。
【図２】フラットパネル型放射線検出器（ＦＰＤ）の検出面の模式図である。
【図３】低周波画像生成部の動作の説明に供する図であり、（ａ）は縮小前の画素を示し、（ｂ）は縮小後の画素を示す。
【図４】境界選択部の動作の説明に供する図である。
【図５】境界選択部の動作の説明に供する図であり、（ａ）は進行方向が縦横である場合の探索する３画素を示し、（ｂ）は進行方向が斜めである場合の探索する３画素を示す。
【図６】領域抽出処理部の動作の説明に供する図であり、（ａ）は拡大前の境界画素を示し、（ｂ）は拡大後の境界画素を示す。
【図７】放射線撮影装置の動作を示すフローチャートである。
【図８】境界選択部の動作を示すフローチャートである。
【図９】変形例に係る境界選択部の動作の説明に供する図である。（ａ）は輝度変化画像における動作を示し、（ｂ）は低周波画像における動作を示し、（ｃ）は選択された境界画素を示す。
【図１０】変形例に係る境界選択部の動作の説明に供する図であり、（ａ）は進行方向が縦横である場合の探索する５画素を示し、（ｂ）は進行方向が斜めである場合の探索する５画素を示す。
【実施例】
【００２３】
以下、図面を参照してこの発明の実施例を説明する。図１は、実施例に係る放射線撮影装置の概略構成を示す図であり、図２は、フラットパネル型放射線検出器（ＦＰＤ）の検出面の模式図である。図３は、低周波画像生成部の動作の説明に供する図であり、（ａ）は縮小前の画素を示し、（ｂ）は縮小後の画素を示す。図４は、境界選択部の動作の説明に供する図である。図５は、境界選択部の動作の説明に供する図であり、（ａ）は進行方向が縦横である場合の探索する３画素を示し、（ｂ）は進行方向が斜めである場合の探索する３画素を示す。図６は、領域抽出処理部の動作の説明に供する図であり、（ａ）は拡大前の境界画素を示し、（ｂ）は拡大後の境界画素を示す。
【００２４】
＜放射線撮影装置＞
図１を参照する。本実施例に係る放射線撮影装置１は、被検体Ｍを載置する天板２と、被検体Ｍに向けて放射線（例えばＸ線）を照射する放射線源（例えばＸ線管）３と、放射線源３と対向して配置され、被検体Ｍを透過した放射線を検出するフラットパネル型放射線検出器（以下、「ＦＰＤ」と略記する）４とを備えている。
【００２５】
放射線源３は、図示しない照射制御部で制御される。照射制御部は、例えばＸ線を照射する場合、後述する入力部８で予め設定された管電圧・管電流などの照射条件に応じて、放射線源３に対して照射に必要な制御を実行する。
【００２６】
ＦＰＤ４は、図２に示すように、その検出面に検出素子ｄが２次元マトリックス状に配列して構成されている。検出素子ｄは、被検体Ｍを透過した放射線の強度に応じて電気信号に変換して一旦蓄積し、その蓄積された電気信号を読み出すことで放射線を検出している。各々の検出素子ｄで検出された電気信号を各画素の画素値として割り当てることで放射線画像を出力している。
【００２７】
図１に戻る。ＦＰＤ４の後段には、図示しないＡ／Ｄ変換器を介して、画像処理部５が設けられている。画像処理部５は、ＦＰＤ４で撮影した放射線画像に種々の画像処理を行う。そして、種々の画像処理が行われた放射線画像は、モニタ等で構成される表示部６に表示される。
【００２８】
また、放射線撮影装置１は、主制御部７、入力部８およびメモリ部９を備えている。主制御部７は、中央演算処理装置（ＣＰＵ）などによって構成され、各種のプログラムを実行することで、放射線源３、ＦＰＤ４および画像処理部５等の各構成を適切に動作させる制御を行う。入力部８では操作者により放射線の照射条件等の入力設定が行われる。メモリ部９は、ＲＯＭ（Read-only Memory）、ＲＡＭ（Random-Access Memory）、またはハードディスク等の記憶媒体で構成され、取得した放射線画像等を記憶する。
【００２９】
画像処理部５は、ＣＰＵなどで構成されている。なお、種々の画像処理を行うためのプログラム等をＲＯＭやＲＡＭなどの記憶媒体に書き込んで記憶し、その記憶媒体からプログラム等を読み出して画像処理部５のＣＰＵが実行することで、そのプログラムに応じた画像処理を行う。画像処理部５は、撮影された放射線画像の中から被検体Ｍの関心領域を抽出する画像領域抽出部１１を備えている。なお、画像領域抽出部１１が本発明における放射線画像領域抽出装置に相当する。
【００３０】
＜画像領域抽出部＞
画像領域抽出部１１は、メモリ部１３、低周波画像生成部１５、輝度変化画像生成部１７、境界選択部１９および領域抽出処理部２１を備えている。メモリ部１３は、ＦＰＤ４で取得した放射線画像を記憶する。
【００３１】
低周波画像生成部１５は、メモリ部１３に記憶された放射線画像を入力画像として、放射線画像を被検体Ｍ内部の空間周波数より小さい周波数フィルタ処理によって、空間分解能を劣化させた（ぼやけた）低周波画像を生成する。低周波画像生成部１５は、例えば、図３に示すように、放射線画像を縮小することにより低周波画像を生成する。具体的な例として、胸部正面画像の肺野領域を抽出する場合では、低周波画像（縮小画像）のサイズが元の放射線画像のサイズの１／８〜１／１０程度であることが好ましい。例えば、図３（ａ）に示すように、１／８に縮小するときは、放射線画像中の縦横８×８画素の平均画素値を、図３（ｂ）のように、１画素として縮小する。なお、図３のハッチングが同一の箇所は、縮小前後の画素を示している。
【００３２】
輝度変化画像生成部１７は、低周波画像生成部１５で生成された低周波画像の輝度変化を検出して輝度変化画像を生成する。輝度変化画像生成部１７は、例えば、特徴量として輝度勾配を算出する。輝度勾配の演算は、画素上の座標（ｘ、ｙ）における画素値をＩ（ｘ、ｙ）とし、算出される特徴量のうち、ｘ軸方向の特徴量をＥｘ（ｘ、ｙ）とし、ｙ軸方向の特徴量をＥｙ（ｘ、ｙ）とすると、下記（１）式および（２）式のように表される。
Ｅｘ（ｘ、ｙ）＝Ｉ（ｘ＋１、ｙ）−Ｉ（ｘ−１、ｙ）      …（１）
Ｅｙ（ｘ、ｙ）＝Ｉ（ｘ、ｙ＋１）−Ｉ（ｘ、ｙ−１）      …（２）
【００３３】
輝度勾配の強度をＰ（ｘ、ｙ）とすると、Ｐ（ｘ、ｙ）は下記（３）式より算出される。
Ｐ（ｘ、ｙ）＝√〔｛Ｅｘ（ｘ、ｙ）｝２＋｛Ｅｙ（ｘ、ｙ）｝２〕      …（３）
【００３４】
境界選択部１９は、低周波画像生成部１５および輝度変化画像生成部１７で生成された輝度変化画像および低周波画像の各画素の画素値に基づいて、関心領域の境界に位置する境界画素を選択する。
【００３５】
境界選択部１９は、輝度変化画像中から輝度変化の最も大きい画素である初期画素を選択し、選択した初期画素に隣接する全方向の画素のうち輝度変化の最も大きい画素を次の画素として選択し、その後、選択した画素の進行方向に隣接する画素から輝度変化の最も大きい画素を次の画素として選択することを繰り返して領域の境界画素を選択する。
【００３６】
境界選択部１９は、まず、輝度変化画像の中から輝度変化が最も大きい画素である初期画素を選択する。例えば、図４に示すように、輝度変化画像の中から輝度変化が最も大きい画素３１を選択する。なお、輝度変化画像が放射線源３およびＦＰＤ４の放射線が図示しないコリメータによって制限されている場合は、輝度変化画像の照射野範囲の中から初期画素３１を選択する。次に、選択した初期画素３１に隣接する全方向の８画素から輝度変化の最も大きい画素を次の境界画素として選択する。例えば、図４に示すように、初期画素３１の周辺８画素を探索し、輝度変化の最も大きい画素３２を次の境界画素として選択する。
【００３７】
境界選択部１９は、その後、選択した画素の進行方向に隣接する画素から輝度変化の最も大きい画素を次の境界画素として選択する。具体的には、図５に示すように、境界画素からなる境界線を延長する進行方向に隣接する３画素を探索する。例えば、図５（ａ）に示すように、現在選択している画素を示すｎ番目に選択した画素をＰ（ｎ）とし、この１つ前に選択された画素をＰ（ｎ−１）とする場合、境界線は、Ｐ（ｎ−１）からＰ（ｎ）へ向かうので、この方向が進行方向となる。選択した画素の進行方向に隣接する３画素は、Ｐ（ｎ）を挟んでＰ（ｎ−１）の反対側の画素Ｂを中心に、画素Ｐ（ｎ）を取り囲むように画素Ｂと隣接する画素Ａおよび画素Ｃで構成される。すなわち、画素Ａ、画素Ｂおよび画素Ｃの３画素を探索する。そして、３画素のうち輝度変化の最も大きい画素を次の画素として選択する。なお、図５（ｂ）も同様である。例えば、図４に示すように、画素Ｐ（ｎ）を画素３２とし、Ｐ（ｎ−１）を画素３１とする場合、画素３２の進行方向に隣接する３画素を探索し、そのうちの輝度変化の最も大きい例えば画素３３を次の境界画素として選択する。
【００３８】
境界選択部１９は、選択した画素（例えば画素３３）が既に選択している画素であるかどうかを判定し、選択した画素でない場合には、選択した画素の進行方向に隣接する３画素から次の境界画素を選択する処理を繰り返し行う。一方、既に選択した画素である場合は、境界選択部１９による処理を終了して、例えば、図４に示すような、選択した境界画素を示す画像データを出力する。
【００３９】
なお、境界選択部１９から出力される境界画素として、図４に示すような境界画素の画像データに限定されない。すなわち、関心領域の内部と外部とを示す２値画像として出力するようにしてもよい。また、選択した境界画素の座標データを保持して座標データ群として出力するようにしてもよい。
【００４０】
領域抽出処理部２１は、境界選択部１９で選択された境界画素を、低周波画像生成部１５に入力される前の放射線画像である元の放射線画像のサイズに拡大する処理を行うとともに、元の放射線画像の中から被検体Ｍの関心領域を抽出する処理を行う。領域抽出処理部２１には、境界選択部１９で選択された境界画素と、メモリ部１３から転送された元の放射線画像と、低周波画像生成部１５から転送された、低周波画像を生成した際の画像縮小率の信号とが転送される。
【００４１】
すなわち、領域抽出処理部２１は、境界画素、元の放射線画像および画像縮小率に基づいて、元の放射線画像の中から被検体Ｍの関心領域を抽出する処理を行う。まず、境界画素と画像縮小率に基づいて、元の放射線画像のサイズに拡大する処理を行う。例えば、１／８の画像縮小率で縮小して低周波画像を生成した場合、境界画素データを８倍に拡大する。具体的には、図６（ａ）に示す境界画素の１画素分を、図６（ｂ）に示す縦横８×８画素に拡大する。このとき、８×８画素のうち、中央に位置する任意の幅（例えば、１画素分）の画素列を境界画素とする。境界画素の幅を狭く設定することで、境界画素付近の画像情報を多く残すことができる。
【００４２】
なお、この境界画素の位置は、拡大した画像の中央に限定されない。例えば、図６の右側を領域の内側とする場合、拡大された８×８画素のうち、抽出する外側の領域を示す左側の画素列を境界画素としたときは、境界とその境界の外側を含んだ領域を抽出することができるので、領域抽出が正確に行われているかを確認できる。また、拡大された８×８画素のうち、抽出する内側を示す右側の画素列を境界画素したときは、境界やその境界の外側などの余分な部分を除いた領域を抽出することができる。右端の画素列を境界画素としたときは、８×８画素全部を境界画素としたときと同様の効果がある。また、画像縮小率の信号は、低周波数画像生成部１５から転送されることに限定されず、例えば、主制御部７から与えるように構成してもよい。
【００４３】
次に、領域抽出処理部２１は、拡大された境界画素データに基づいて、元の放射線画像の中から被検体Ｍの関心領域を抽出する処理を行う。すなわち、領域抽出処理部２１は、元の放射線画像の中から拡大した境界画像の位置における領域を分割することにより、被検体Ｍの関心領域を抽出する。また、領域抽出処理部２１は、例えば、放射線画像内に照射野領域を表す線を表示するようにしたり、抽出した関心領域外の画素を暗くするように画素値を設定したりする（例えば画素値を０にする）。画像処理部５は、画像領域抽出部１１で被検体Ｍの関心領域が抽出された放射線画像に対し、コントラストの操作などその他必要な画像処理を行う。画像処理部５でその他必要な画像処理が行われた放射線画像は、表示部６に表示されたり、メモリ部９に保存されたりする。
【００４４】
＜動作の説明＞
次に、放射線撮影装置１について、特に画像領域抽出部１１の動作を図７および図８のフローチャートに沿って説明する。
【００４５】
〔ステップＳ１〕撮影
図１において図示しない撮影ボタンを押すことにより、放射線源３から放射線（例えばＸ線）を照射し、これに連動してＦＰＤ４によって撮影が開始される。放射線源３から照射された放射線は被検体Ｍを透過し、被検体Ｍを透過した放射線をＦＰＤ４で検出する。ＦＰＤ４は、検出した放射線の強度に応じて電気信号に変換して放射線画像を出力する。ＦＰＤ４から出力された放射線画像は、図示しないＡ／Ｄ変換器でディジタルデータに変換されて画像処理部５の画像領域抽出部１１に転送される。
【００４６】
〔ステップＳ２〕低周波画像の生成
画像領域抽出部１１において、先ず、ＦＰＤ４で撮影された放射線画像は、メモリ部１３に一旦記憶される。メモリ部１３に記憶された放射線画像は、読み出されて低周波画像生成部１５に転送される。低周波画像生成部１５は、放射線画像を被検体Ｍ領域内部の空間周波数より低い周波数フィルタ処理によって、放射線画像の空間分解能を劣化させた低周波画像を生成する。低周波画像生成部１５は、入力された放射線画像を縮小して縮小画像を生成することにより低周波画像を生成する。放射線画像の画像縮小率は、１／８〜１／１０程度が好ましい。低周波画像を生成することにより、ノイズ成分や被検体Ｍ内部の細かい構造（例えば肋骨等）の影響を低減しているので安定し、被検体Ｍの輪郭形状の認識、例えば、関心領域の境界に位置する境界画素の選択を容易に行うことができる。なお、低周波画像生成部１５で生成された低周波画像は、輝度変化画像生成部１７および境界選択部１９に転送される。
【００４７】
〔ステップＳ３〕輝度変化画像の生成
輝度変化画像生成部１７は、低周波画像生成部１５で生成された低周波画像の各画素の画素値に基づいて、各画素における輝度変化量を求めることで輝度変化画像を生成する。輝度変化量は、画像の輝度勾配により算出される。輝度勾配の算出は、上述した式（１）〜式（３）によって行われる。
【００４８】
〔ステップＳ４〕境界画素の選択
境界選択部１９は、ステップＳ２の低周波画像生成部１５で生成された低周波画像、およびステップＳ３の輝度変化画像生成部１７で生成された輝度変化画像の画素値に基づいて、関心領域の境界に位置する境界画素を選択する。ここで図８のフローチャートを参照して、境界選択部１９による境界画素の選択の動作について説明する。
【００４９】
〔ステップＳ４１〕初期画素選択（輝度変化画像内）
境界選択部１９は、先ず、輝度変化画像の中から（図１において図示しないコリメータ等で照射野が制限されている場合には輝度変化画像の照射野範囲の中から）、最も変化の大きい画素、すなわち、輝度変化画像の中から画素値が最も大きい画素を選択する。そして、境界選択部１９は、その選択した画素を領域の境界に位置する画素の初期画素（初期点）として設定する。
【００５０】
〔ステップＳ４２〕画素選択（輝度変化画像内）
選択している画素に隣接する周囲の画素を探索し、その隣接する画素から輝度変化の最も大きい画素を選択する。ここで、選択している画素が初期画素である場合は、初期画素に隣接する全方向の８画素を探索し、それら８画素から輝度変化の最も大きい画素を次の境界画素として選択する。また、選択している画素が初期画素でない場合は、図５に示すような、選択している画素Ｐ（ｎ）の進行方向に隣接する３画素（画素Ａ、画素Ｂおよび画素Ｃ）を探索し、それら３画素から輝度変化の最も大きい画素を次の境界画素として選択する。
【００５１】
〔ステップＳ４３〕第１判定処理
ステップＳ４２で新たに選択した画素が境界画素として信頼性があるかどうかの判定を行う。すなわち、輝度変化画像の新たに選択した画素の画素値が予め設定した所定値以下の場合は、新たに選択した画素に対し、境界画素としての信頼性が低いとして判定する（ステップＳ４５に進む）。一方、輝度変化画像の新たに選択された画素が予め設定した所定値より大きい場合は、ステップＳ４４に進む。輝度変化が小さいほど境界画素として信頼性が低いことを示している。仮に誤って輝度変化が所定値以下の画素を選択した場合に、境界画素として新たに選択した画素の信頼性が低いとして振り分けられるので、より高精度で安定した領域抽出を行うことができる。
【００５２】
〔ステップＳ４４〕第２判定処理
ステップＳ４３で判定処理を行った新たに選択した画素に対し、さらに判定を行う。すなわち、輝度変化画像の新たに選択した画素および選択前の画素と同位置における低周波画像の各画素の画素値を比較して予め設定した所定値よりも大きい場合は、新たに選択した画素に対し、境界画素としての信頼性が低いとして判定する（ステップＳ４５に進む）。一方、予め設定した所定値以下の場合は、ステップＳ４６に進む。なお、選択前の画素とは、新たに選択した画素の選択前の画素をいう。また、比較は、低周波画像において、新たに選択した画素の画素値と選択前の画素の画素値との差を求めて行われる。低周波画像において隣接する画素の画素値が近いほど、同じ領域内にある可能性が高いことを示している。また、新たに選択した画素が輝度変化の大きい画素であっても（新たに選択した画素の画素値と選択前の画素の画素値とが同じような大きさであっても）、低周波画像において画素値が大きく異なることがある。この場合、境界画素として新たに選択した画素の信頼性が低いとして振り分けられるので、より高精度で安定した領域抽出を行うことができる。
【００５３】
なお、ステップＳ４３およびＳ４４の判定処理において、予め設定する所定値は固定値であるが、これに限定されない。すなわち、画像の画素値に応じて決定される値を設定してもよい。例えば、低周波画像の画素値範囲の１０％を所定値として設定するようにしてもよい。このとき、低周波画像の画素値範囲が１２０であれば、所定値として設定される値は１２０の１０％である１２となる。
【００５４】
〔ステップＳ４５〕画素選択（低周波画像内）
ステップＳ４３およびステップＳ４４の判定処理において、信頼性が低いと判定された場合は、低周波画像の中から境界画素を選択する。すなわち、輝度変化画像の選択前の画素と同位置における低周波画像の画素の画素値と、その画素に隣接する画素から最も値が近い画素を次の境界画素として選択する。低周波画像の画素の画素値が近い値の画素である程、同じ領域内でかつ境界画素である可能性が高いことを示している。信頼性が低いとして判定された場合でも、境界画素である可能性の高い画素を選択することができる。そのため、より高精度で安定した領域抽出を行うことができる。なお、境界画素の選択は、ステップＳ４２に示すように、選択している画素の全方向または進行方向に隣接する３画素を探索して行う。
【００５５】
〔ステップＳ４６〕既に選択した画素であるかどうかの判定
ステップＳ４２、Ｓ４５において新たに選択した画素が、既に選択している画素であるかどうかを判定する。既に選択した画素でない場合には、ステップＳ４２に戻って、同様の処理を繰り返す。一方、既に選択した画素である場合は、境界選択部１９の処理を終了する。図７に戻る。図７のステップＳ４が終了し、ステップＳ５に進む。なお、このように繰り返し境界画素の選択を行うことにより、被検体Ｍの関心領域の輪郭を示す境界画素が選択される。選択後、境界選択部１９より境界画素データとして出力される。
【００５６】
〔ステップＳ５〕領域抽出処理
領域抽出処理部２１には、ステップＳ４で選択された境界画素データ、メモリ部１３に記憶された低周波画像に生成される前の放射線画像（元の放射線画像）、および低周波画像を生成した際の画像縮小率が転送される。領域抽出処理部２１は、境界画素データと画像縮小率とに基づいて元の放射線画像のサイズに拡大する処理を行い、拡大された境界画素データと元の放射線画像とに基づいて、元の放射線画像の中から関心領域を抽出する処理を行う。
【００５７】
〔ステップＳ６〕コントラスト等の調整
画像領域抽出部１１で被検体Ｍの関心領域が抽出された放射線画像は、画像処理部５により、関心領域内のコントラストの調整など、その他必要な画像処理が行われる。
【００５８】
〔ステップＳ７〕画像表示および保存
画像処理部５で画像処理された、被検体Ｍの関心領域が抽出された放射線画像は、表示部６に表示される。また、メモリ部９に保存される。
【００５９】
このような放射線撮影装置１によれば、低周波画像生成部１５は、放射線画像を被検体Ｍ領域内部の空間周波数より低い周波数フィルタ処理によって低周波画像を生成し、輝度変化画像生成部１７は、低周波画像の輝度変化を検出して輝度変化画像を生成する。すなわち、低周波画像を生成することで、ノイズ成分や被検体Ｍ内部の細かい構造成分の影響を低減し、被検体Ｍの輪郭形状の認識を容易にしている。また、低周波画像から生成された輝度変化画像は、低周波画像と同様の効果を有している。境界選択部１９は、輝度変化画像および低周波画像を組み合わせて、被検体Ｍの関心領域の境界に位置する境界画素を選択しているので、高精度で安定した領域抽出を行うことができる。
【００６０】
また、低周波画像生成部１５は、放射線画像を縮小することにより低周波画像を生成している。それにより、通常のフィルタに比べてより低周波な低周波画像を生成することができるので、ノイズなどの外乱の影響が少なく、高精度で安定した領域抽出を行うことができる。また、低周波画像および輝度変化画像のデータ量が比較的少ないので、被検体Ｍの関心領域の抽出を高速に行うことができる。
【００６１】
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。
【００６２】
（１）上述した実施例では、境界選択部１９は、初期画素を選択し、選択した画素に隣接する画素から次の境界画素を選択していたが、この方法に限定されない。すなわち、輝度変化画像の画素値が予め設定した所定値より大きい画素で、かつ低周波画像の画素の画素値により検出された画素に基づいて、関心領域の境界に位置する境界画素を選択するようにしてもよい。
【００６３】
例えば、まず、図９（ａ）に示すように、輝度変化画像から輝度変化が予め設定した所定値よりも大きい画素を、例えばラベリング処理等して検出する。所定値よりも大きい画素を太線で囲った画素で示す。次に、図９（ｂ）に示すように、輝度変化画像で検出された輝度変化が所定値より大きい画素と同位置における低周波画像の画素の画素値に基づいて境界画素を例えばラベリング処理等して選択する。例えば、肺野領域を抽出する場合は、肺野領域が周囲と比較して明るいので、低周波画像の中から予め設定した所定値よりも大きい画素を選択して境界画素を選択する。選択した境界画素を図９（ｃ）に示す。
【００６４】
なお、抽出する被検体Ｍの関心領域が周囲と比較して暗い場合は、図９（ｂ）において、低周波画像の中から予め設定した所定値以下の画素値の画素を選択して境界画素として選択する。また、低周波画像の中から予め設定した所定の範囲内の画素値の画素を境界画素として選択してもよい。
【００６５】
また、図９（ａ）において、輝度変化画像から輝度変化が予め設定した所定値よりも大きい画素を例えばラベリング処理等して検出し、図９（ｂ）において、低周波画像の中から予め設定した所定値よりも大きい画素を検出し、両方の画像から検出された画素を選択画素として選択するようにしてもよい。
【００６６】
（２）上述した実施例では、放射線画像を縮小することにより低周波画像を生成する方法（低周波画像生成部１５による縮小処理）は、例えば、放射線画像を１／８に縮小する場合、８×８画素の画素値の平均値を求めて、その平均値を縮小画像の１画素として低周波画像を生成したが、この方法に限定されない。すなわち、放射線画像の画素を間引くことにより縮小する方法でもよい。例えば、放射線画像を１／８に縮小する場合、８×８画素のうちの間引いた残りの画素（例えば、図３の符号Ｔに示すように４×４画素）の平均画素値を１画素として低周波画像を生成する。間引くことで、縮小する際の演算するデータ量を少なくすることができる。また、その他通常において用いられる縮小方法であってもよい。
【００６７】
（３）上述した実施例では、低周波画像生成部１５は、放射線画像を縮小することにより低周波画像を生成していたが、この方法に限定されない。例えば、画素を１つずつずらしながら処理するフィルタ処理による平滑化によって低周波画像を生成する方法、放射線画像を空間周波数領域に変換して、変換された周波数領域のうちの高周波領域をカットして低周波領域を実空間に変換することにより低周波画像を生成する方法、あるいは、これらの方法の組合せなど、通常において用いられる低周波画像を生成する手法であればよい。
【００６８】
上述の、放射線画像を空間周波数領域に変換して、変換された周波数領域のうちの高周波領域をカットして低周波領域を実空間に変換することにより低周波画像を生成する方法は以下のような方法がある。すなわち、放射線画像をフーリエ変換により空間周波数領域に変換して、変換された空間周波数領域から高周波成分を除去する、あるいは低周波成分のみを通過させることで低周波数領域を求め、その低周波数領域を逆フーリエ変換により実空間に変換することにより低周波画像を生成してもよい。このフーリエ変換、逆フーリエ変換を用いた方法の場合には、画像サイズを変えずに低周波数画像を生成することが可能である。また、フーリエ変換に限定されず、ウェーブレット変換やガボールフィルタのなどの処理を行ってもよい。
【００６９】
（４）上述した実施例では、輝度変化画像生成部１７は、低周波画像のそれぞれの画素位置で、その画素に隣接する上下・左右の画素値の差に基づき輝度変化量（特徴量）として放射線画像の輝度勾配を求めていたが、これに限定されない。すなわち、輝度勾配はその他公知の方法を使用しても構わない。また、輝度変化量の算出は、例えば、Sobel、Prewitt、Laplacian、Cannyなどの１次元または２次元の微分フィルタ処理によるエッジ算出、周波数処理による高周波成分抽出など、画像の輝度変化を特徴として抽出する方法なら他の処理を行ってもよい。
【００７０】
（５）上述した実施例では、境界選択部１９は、選択した画素の進行方向に隣接する３画素を探索して境界画素を選択していたが、これに限定されない。すなわち、境界選択部１９は、図１０に示すように、選択した画素の進行方向に隣接する５画素を探索して境界画素を選択するようにしてもよい。
【００７１】
（６）上述した実施例では、被検体Ｍの関心領域として胸部の肺野領域の抽出を例に採って説明したが、これに限定されない。すなわち、抽出する被検体Ｍの関心領域は、その関心領域の内側と外側で画像濃度やコントラストが大きく異なる他の部位であってもよい。
【００７２】
（７）上述した実施例では、放射線撮影装置１は、放射線検出器としてフラットパネル型放射線検出器（ＦＰＤ）４を備えていたが、これに限定されない。すなわち、放射線撮影装置１は、放射線検出器としてイメージインテンシファイアを備えていてもよい。
【００７３】
（８）上述した実施例では、放射線としてＸ線を例に採って説明したが、これに限定されない。すなわち、放射線はγ線であってもよい。例えば、放射線検出器は、放射性同位元素（ＲＩ）を投与された被検体Ｍから放射されるγ線を検出してもよい。
【００７４】
（９）上述した実施例では、放射線撮影装置１は、単に放射線源３とＦＰＤ４を備えていたが、これに限定されない。すなわち、天板２に載置された被検体Ｍの体軸を中心に放射線源３とＦＰＤ４が回転して断層画像を撮影するＸ線ＣＴ装置であってもよい。また、工業用に用いられる非破壊検査装置のように被検体Ｍ（検査の対象物）をベルト上に運搬させて撮影を行う構成のものであってもよい。
【符号の説明】
【００７５】
１    …  放射線撮影装置
３    …  放射線源
４    …  フラットパネル型放射線検出器（ＦＰＤ）
５    …  画像処理部
１１  …  画像領域抽出部
１５  …  低周波画像生成部
１７  …  輝度変化画像生成部
１９  …  境界選択部
２１  …  領域抽出処理部
【図１】
2011255033.tif 000003
【図２】
2011255033.tif 000004
【図３】
2011255033.tif 000005
【図４】
2011255033.tif 000006
【図５】
2011255033.tif 000007
【図６】
2011255033.tif 000008
【図７】
2011255033.tif 000009
【図８】
2011255033.tif 000010
【図９】
2011255033.tif 000011
【図１０】
2011255033.tif 000012