(57)【要約】【課題】体動等を検出するための解析点を放射線画像上に効率よく設定する。【解決手段】被写体領域抽出部５１が放射線画像から被写体領域を抽出する。解析点抽出部５２が、放射線画像を多重解像度変換することにより生成した低周波成分において、被写体領域の画素値を特定のしきい値と比較する。画素値がしきい値を超える画素を、例えば体動検出のための解析点として抽出する。【選択図】図４
2014064857.tif 000002
【特許請求の範囲】
【請求項１】
放射線画像から対象となる被写体の領域を抽出する被写体領域抽出手段と、
前記被写体の領域における空間周波数の低周波成分が特定のしきい値を超える点を、体動を検出するための処理を行うための解析点として抽出する解析点抽出手段とを備えたことを特徴とする解析点設定装置。
【請求項２】
前記しきい値は、前記被写体の領域における画素値の標準偏差に基づいて決定されることを特徴とする請求項１記載の解析点設定装置。
【請求項３】
請求項１または２記載の解析点設定装置と、
前記解析点において、前記放射線画像の高周波成分のコントラストおよび低周波成分のコントラストを算出するコントラスト算出手段と、
前記高周波成分のコントラストと前記低周波成分のコントラストとの比率を算出する比率算出手段と、
該比率に基づいて前記放射線画像に含まれる被写体の体動の有無を判定する判定手段とを備えたことを特徴とする体動検出装置。
【請求項４】
放射線画像から対象となる被写体の領域を抽出し、
前記被写体の領域における空間周波数の低周波成分が特定のしきい値を超える点を、体動を検出するための処理を行うための解析点として抽出することを特徴とする解析点設定方法。
【請求項５】
請求項３記載の解析点設定方法により前記放射線画像に解析点を設定し、
前記解析点において、前記放射線画像の高周波成分のコントラストおよび低周波成分のコントラストを算出し、
前記高周波成分のコントラストと前記低周波成分のコントラストとの比率を算出し、
該比率に基づいて前記放射線画像に含まれる被写体の体動の有無を判定することを特徴とする体動検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線画像に体動を検出するための解析点を設定する解析点設定装置および方法、並びに設定した解析点を用いた体動検出装置および方法に関するものである。
【背景技術】
【０００２】
被写体にＸ線等の放射線を照射して被写体を撮影し、被写体の放射線画像を取得する際に、撮影中に被写体が動いてしまう体動により、取得される放射線画像にボケが生じ、画質が劣化してしまう場合がある。とくに、撮影対象が、頚部、胸部、腹部および腰部等の放射線の照射時間が長い部位であったり、小児を撮影したりする場合には、体動が生じやすいことから、放射線画像の劣化が問題となる。
【０００３】
このような体動の有無は、撮影を行う技師が、撮影中の被写体の様子を観察し、かつ取得された放射線画像を直ちに確認することにより、見極めることができる。しかしながら、撮影室にあるモニタは解像度が低く、放射線画像が縮小して表示され、かつ観察環境が比較的明るい場合が多いため、表示された放射線画像を観察して体動の有無を判断することは困難である。このような場合、放射線画像を拡大することが考えられるが、放射線画像を拡大する作業は、技師のワークフローを低下させることとなる。
【０００４】
このため、放射線画像に基づいて体動を自動で検出する手法が提案されている。例えば、放射線画像に含まれるエッジの劣化を、縦方向横方向の２方向の特徴量として学習しておき、放射線画像からエッジを抽出し、学習した特徴量と抽出したエッジとに基づいて、体動によるブレが放射線画像に生じているか否かを判定する手法が提案されている（特許文献１参照）。また、頭部の側面という特定の部位についての特徴点の学習を行い、学習結果を用いて取得した放射線画像の体動を検出する手法も提案されている（非特許文献１参照）。
【０００５】
このような体動を検出するためには、体動を検出するための解析点を放射線画像上に設定する必要がある。このため、技師は表示された放射線画像上にマニュアル操作により解析点を設定して、体動を検出することができる。また、体動検出のための高空間周波数を持つ部位の部分画素を解析点として、自動で抽出する手法も提案されている（特許文献２参照）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許第７８９９２２９号明細書
【特許文献２】特開２００５−２０５０８６号公報
【非特許文献１】Proc. Of SPIE Vol.6914, 69140U, (2008)
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、撮影中に体動が生じた場合、取得される放射線画像においては、体動ブレにより高周波成分のコントラストが低下するため、特許文献２に記載された手法では、体動が生じている場合には、精度良く解析点を抽出することができない。
【０００８】
本発明は上記事情に鑑みなされたものであり、体動を検出するための解析点を放射線画像上に精度良く設定することを目的とする。
【課題を解決するための手段】
【０００９】
本発明による解析点設定装置は、放射線画像から対象となる被写体の領域を抽出する被写体領域抽出手段と、
前記被写体の領域における空間周波数の低周波成分が特定のしきい値を超える点を、体動を検出する処理を行うための解析点として抽出する解析点抽出手段とを備えたことを特徴とするものである。
【００１０】
なお、本発明による解析点設定装置においては、しきい値を、前記被写体の領域における画素値の標準偏差に基づいて決定してもよい。
【００１１】
本発明による体動検出装置は、本発明の解析点設定装置と、
前記解析点において、前記放射線画像の高周波成分のコントラストおよび低周波成分のコントラストを算出するコントラスト算出手段と、
前記高周波成分のコントラストと前記低周波成分のコントラストとの比率を算出する比率算出手段と、
該比率に基づいて前記放射線画像に含まれる被写体の体動の有無を判定する判定手段とを備えたことを特徴とするものである。
【００１２】
本発明による解析点設定方法は、放射線画像から対象となる被写体の領域を抽出し、
前記被写体の領域における空間周波数の低周波成分が特定のしきい値を超える点を、体動を検出するための処理を行うための解析点として抽出することを特徴とするものである。
【００１３】
本発明による体動検出方法は、本発明の解析点設定方法により前記放射線画像に解析点を設定し、
前記解析点において、前記放射線画像の高周波成分のコントラストおよび低周波成分のコントラストを算出し、
前記高周波成分のコントラストと前記低周波成分のコントラストとの比率を算出し、
該比率に基づいて前記放射線画像に含まれる被写体の体動の有無を判定することを特徴とするものである。
【００１４】
なお、本発明による解析点設定方法および体動検出方法をコンピュータに実行させるためのプログラムとして提供してもよい。
【発明の効果】
【００１５】
本発明の解析点設定装置および方法によれば、放射線画像から対象となる被写体の領域が抽出され、被写体の領域における空間周波数の低周波成分が特定のしきい値を超える点が、体動を検出する処理を行うための解析点として抽出される。ここで、撮影中に体動が生じた場合、取得される放射線画像においては、高周波成分のコントラストは低下するが、低周波成分のコントラストは殆ど低下しない。このため、本発明によれば、体動に影響されることなく、体動を検出する処理を行うための解析点を、放射線画像に精度良く設定することができ、その結果、精度良く体動を検出することができる。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施形態による解析点設定装置および体動検出装置を適用した放射線画像撮影システムの構成を示す概略図
【図２】本実施形態において行われる処理を示すフローチャート
【図３】多重解像度変換を説明するための図
【図４】解析点設定部の構成を示す概略ブロック図
【図５】放射線画像の例を示す図
【図６】設定した解析点の例を示す図
【図７】解析点を中心とする８方向を示す図
【図８】解析領域および水平方向のラインを説明するための図
【図９】体動ありと判定された場合の判定結果の表示画面の例を示す図
【図１０】体動ありと判定された場合の判定結果の表示画面の例を示す図
【図１１】体動ありと判定された場合の判定結果の表示画面の例を示す図
【図１２】横軸にコントラスト比の値を、縦軸に各コントラスト比の値において体動ありと判定した画像の数、および体動なしと判定した画像の数をプロットしたヒストグラムを示す図
【発明を実施するための形態】
【００１７】
以下、本発明の実施形態について図面を参照して詳細に説明する。図１は、本発明の実施形態による解析点設定装置および体動検出装置を適用した放射線画像撮影システムの構成を示す概略図である。本実施形態による放射線画像撮影システムは、被写体を撮影することにより、被写体の放射線画像Ｇ０を取得する撮影装置１０と、放射線画像Ｇ０に基づいて体動を検出する、本実施形態の解析点設定装置および体動検出装置を備えた信号処理部２０と、信号処理部２０に各種指示を与えるための入力部３０と、撮影により取得された放射線画像等を表示する表示部４０とから構成されている。
【００１８】
撮影装置１０は、被写体Ｈに放射線Ｒを照射する放射線管球１２、放射線管球１２の駆動を制御する撮影制御部１３および被写体Ｈが載置される撮影台１４を備える。撮影台１４には、被写体Ｈの放射線検出信号を出力する放射線検出器１１が設けられている。放射線検出器１１は、マトリクス状に配置された画素毎に照射放射線のエネルギーレベルに対応した放射線検出信号を出力するものであり、この検出信号はＡ／Ｄ変換処理され、被写体Ｈの放射線画像Ｇ０を表すデジタル画像信号として出力される。
【００１９】
なお、放射線検出器１１としては、例えば特開平７−７２２５３号公報に記載されているように、放射線の照射を受けて可視光を発するシンチレータと、その可視光を検出する固体光検出素子とが積層されてなるもの、あるいは例えば特開２０１０−２０６０６７号公報に記載されているように、放射線の照射を受けてそのエネルギーに対応した電気信号を出力する放射線光導電層を有してなるものを適用することができる。
【００２０】
信号処理部２０は、放射線画像Ｇ０を表すデジタル画像信号が入力される前処理部２１、多重解像度変換部２２、解析点設定部２３、コントラスト算出部２４、比率算出部２５、判定部２６および表示制御部２７を備える。
【００２１】
入力部３０は、例えばキーボード３１およびマウス３２等からなるものであり、放射線技師等の使用者による各種指示を信号処理部２０に入力する。
【００２２】
表示部４０は、例えば液晶ディスプレイやＣＲＴディスプレイ等からなるものであり、体動の有無の判定結果や、撮影された被写体の放射線画像等を必要に応じて表示する。
【００２３】
以上述べた信号処理部２０、入力部３０および表示部４０は、例えば一般的なパーソナルコンピュータ等のコンピュータシステムから構成することができる。
【００２４】
次いで、放射線画像の撮影について説明する。放射線画像の撮影に際しては、放射線検出器１１が撮影装置１０の撮影台１４の上に載置され、その上に被写体Ｈが置かれる。この状態で撮影制御部１３が操作されることにより放射線管球１２が駆動されて、放射線Ｒが被写体Ｈを透過して放射線検出器１１に照射される。撮影が終了すると、放射線検出器１１から放射線画像Ｇ０を表すデジタル画像信号が取得される。この放射線画像Ｇ０は表示部４０において再生表示可能である。
【００２５】
放射線画像Ｇ０は、信号処理部２０に入力される。以下、信号処理部２０が行う処理について説明する。図２は信号処理部２０において行われる処理を示すフローチャートである。信号処理部２０に入力された放射線画像Ｇ０は、まず前処理部２１において、放射線の照射ムラおよび放射線検出器１１の検出ムラ等に起因する信号値の変動を補正する処理、放射線画像Ｇ０の画質を向上させるための濃度、コントラスト、周波数成分等を補正する画像処理、並びに必要に応じて適宜なされるその他の処理を受ける（前処理、ステップＳＴ１）。
【００２６】
次いで、多重解像度変換部２２が、前処理済みの放射線画像Ｇ０を多重解像度変換し、放射線画像Ｇ０の高周波成分および低周波成分を抽出する（ステップＳＴ２）。図３は多重解像度変換を説明するための図である。まず、多重解像度変換部２２は、放射線画像Ｇ０に対してσ＝１のガウシアンフィルタによりフィルタリング処理を行って、放射線画像Ｇ０を１／２に縮小してガウシアン成分Ｇ１を生成する。ガウシアン成分Ｇ１は放射線画像Ｇ０を１／２に縮小したものとなる。次いで、多重解像度変換部２２は、例えば３次Ｂスプライン補間等の補間演算を行って、ガウシアン成分Ｇ１を放射線画像Ｇ０と同一サイズに拡大し、拡大したガウシアン成分Ｇ１を放射線画像Ｇ０から減算して、最高周波数帯域のラプラシアン成分Ｌ０を生成する。本実施形態においては、このラプラシアン成分Ｌ０を高周波成分として用いる。なお、本実施形態では、最高周波数帯域を便宜上第０の周波数帯域と称する。
【００２７】
次いで、多重解像度変換部２２は、ガウシアン成分Ｇ１に対してσ＝１のガウシアンフィルタによりフィルタリング処理を行って、ガウシアン成分Ｇ１をさらに１／２に縮小してガウシアン成分Ｇ２を生成し、ガウシアン成分Ｇ２をガウシアン成分Ｇ１と同一サイズに拡大し、拡大したガウシアン成分Ｇ２をガウシアン成分Ｇ１から減算して、第１の周波数帯域のラプラシアン成分Ｌ１を生成する。さらに、所望とする周波数帯域のラプラシアン成分が生成されるまで上記の処理を繰り返すことにより、複数の周波数帯域のラプラシアン成分Ｌｊ（ｊ＝０〜ｎ）を生成する。
【００２８】
本実施形態においては、第３の周波数帯域のラプラシアン成分Ｌ３が得られるまで、上記の処理を繰り返し、第１周波数帯域のラプラシアン成分Ｌ０と同一サイズとなるように、ラプラシアン成分Ｌ３を８倍に拡大して、ラプラシアン成分Ｌ３′を生成する。そして、このラプラシアン成分Ｌ３′を低周波成分として用いる。
【００２９】
ここで、ガウシアン成分の各画素の信号値はその画素の濃度を表し、ラプラシアン成分の各画素の信号値は、その画素におけるその周波数帯域の周波数成分の大きさを表すものとなる。
【００３０】
なお、ウェーブレット変換等の他の多重解像度変換の手法を用いることにより、周波数帯域が異なる複数の帯域画像を生成して最高周波数帯域の帯域画像を高周波成分、上記と同様の第３の周波数帯域の帯域画像を８倍に拡大したものを低周波成分として用いてもよい。なお、低周波成分の周波数帯域は第３の周波数帯域に限定されるものではなく、最高周波数帯域よりも低い周波数帯域であれば、任意の周波数帯域の周波数成分を低周波成分として用いることができる。また、高周波成分としては、放射線画像Ｇ０そのものを用いるようにしてもよい。
【００３１】
次いで、解析点設定部２３が、後述するコントラストを算出するための解析点を、前処理が施された放射線画像Ｇ０に設定する（ステップＳＴ３）。具体的には、低周波成分Ｌ３′を用いて解析点を設定する。以下、解析点の設定について説明する。図４は、解析点設定部２３の構成を示す概略ブロック図である。図４に示すように、解析点設定部２３は、被写体領域抽出部５１および解析点抽出部５２を備える。被写体領域抽出部５１は、放射線画像Ｇ０に含まれるスキンラインや照射野枠の影響を除外して、体動を検出する対象を放射線画像Ｇ０に含まれる骨部を含む被写体領域に限定すべく、放射線画像Ｇ０から被写体領域を抽出する。ここで、撮影中に被写体の体動が生じているか否かを判断するためには、放射線画像Ｇ０における高コントラストな領域を抽出する必要がある。このため、本実施形態においては、被写体領域抽出部５１は、骨部を含む領域を被写体領域として抽出する。
【００３２】
具体的には、対象となる被写体についての学習結果を用いる手法、被写体についてのルールベースの手法、グラフカット法等の各種セグメンテーション技術を用いて、放射線画像Ｇ０から被写体領域を抽出する。例えば、図５に示すように、腰椎を含む放射線画像Ｇ０から腰椎の部分を抽出する場合、腰椎の形状について学習を行った学習結果を用いて、放射線画像Ｇ０から腰椎を含むあらかじめ定められた大きさの領域を被写体領域Ａ０として抽出する。
【００３３】
解析点抽出部５２は、被写体領域Ａ０から体動の有無を判定する処理の基準となる解析点を抽出する。解析点は、体動の有無を判定しやすい部分から抽出することが好ましい。具体的には、骨部のエッジ上に解析点を抽出することが好ましい。このため、解析点抽出部５２は、多重解像度変換部２２が生成した低周波成分Ｌ３′における被写体領域Ａ０から比較的高い画素値となる画素を解析点として抽出する。具体的には、解析点抽出部５２は、低周波成分Ｌ３′の被写体領域Ａ０内の各画素の画素値をしきい値Ｔｈ１と比較し、画素値がしきい値Ｔｈ１を超える画素位置を解析点として抽出する。なお、しきい値Ｔｈ１は、例えば、Ｔｈ１＝σ×αにより算出する。ここで、σは低周波成分Ｌ３′の被写体領域Ａ０内の画素値の標準偏差、αは係数である。係数αは、低周波成分Ｌ３′の各画素の画素値に応じて適宜設定すればよい。このようにして設定した解析点の例を図６に示す。
【００３４】
次いで、コントラスト算出部２４が、解析点における高周波成分のコントラストおよび低周波成分のコントラストを算出する（ステップＳＴ４）。ここで、被写体に放射線を照射するわずかな時間（数十ミリ秒〜数百ミリ秒）においては、体動は一方向のブレとして現れる。体動が発生した場合、放射線画像Ｇ０は体動の方向と同一方向に劣化するが、体動の方向とは垂直な方向には劣化しない。このため、コントラスト算出部２４は、まず各解析点におけるエッジの勾配方向を算出する。具体的には、コントラスト算出部２４は、低周波成分Ｌ３′上に、図７に示すように、解析点Ｐ０を中心とする８つの方向を設定する。例えば、紙面の下から上に向かう方向を基準（０度）とした場合、−６７．５度、−４５度、−２２．５度、０度、２２．５度、４５度、６７．５度および９０度の８つの方向を設定する。なお、８つの方向をそれぞれθ１〜θ８とする。また、設定する方向は８に限定されるものではなく、２方向、４方向あるいは１６方向等、任意の数とすることができる。
【００３５】
そして、各方向毎に、解析点Ｐ０を中心とする特定サイズの解析領域を設定する。解析領域は、例えば９×９画素の矩形の領域である。そして、解析領域内に複数のラインを仮想的に設定し、各ラインにおける画素値の最大値と最小値との差分値を各ラインにおけるコントラストとして算出する。ここで、ラインの方向は、その解析領域を設定した方向と一致させる。例えば、図７に示す方向θ１に解析領域を設定した場合、その解析領域に設定するラインの方向はθ１となる。
【００３６】
図８は解析領域および水平方向のラインを説明するための図である。図８に示すように解析点Ｐ０のある方向θ１に設定した解析領域Ａ２に複数のラインが設定される。なお、ライン数は、解析領域Ａ２の垂直方向の画素数と同一とすればよい。すなわち、解析領域Ａ２のサイズが９×９画素であれば、ライン数は９とすればよい。しかしながら、ライン数を適宜間引いて少なくしてもよい。
【００３７】
そして、コントラスト算出部２４は、すべてのラインについてのコントラストの平均値を、低周波成分Ｌ３′においてその解析領域を設定した方向についてのコントラストとして算出する。そして、コントラスト算出部２４は、８つの方向θ１〜θ８のすべてに設定した解析領域のそれぞれについて、コントラストの平均値を各解析領域についてのコントラストＣＬθｉ（ｉ＝１〜８）として算出する。さらに、コントラスト算出部２４は、算出した８つのコントラストＣＬθｉのうち、最も値が高いコントラストを算出した解析領域を設定した方向を、その解析点についてのエッジの勾配方向に決定する。
【００３８】
同様に、コントラスト算出部２４は、高周波成分Ｌ０についても、低周波成分Ｌ３′と同様に、すべての解析点において、８つの方向θ１〜θ８のすべてについての８つのコントラストＣＨθｉ（ｉ＝１〜８）を算出する。そして、算出した８つのコントラストＣＨθｉのうち、最も値が高いコントラストを算出した解析領域を設定した方向を、その解析点についてのエッジの勾配方向に決定する。
【００３９】
なお、各解析点において、コントラストが変化する方向は低周波成分Ｌ３′であっても高周波成分Ｌ０であっても同一となることから、各解析点におけるエッジの勾配方向は一致することとなる。このため、各解析点において同一のエッジ勾配方向についての、それぞれ１つの高周波成分のコントラストＣＨθｉおよび低周波成分のコントラストＣＬθｉが算出される。
【００４０】
ここで、撮影中に体動が生じた場合、取得される放射線画像Ｇ０においては、高周波成分のコントラストは低下するが、低周波成分のコントラストは殆ど低下しない。また、撮影条件あるいは被写体の種類によって、放射線画像全体のコントラストは変動するが、その変動は低周波成分のコントラストの差として現れる。したがって、高周波成分のコントラストと低周波成分のコントラストとの比率は、体動の有無を表す指標として用いることができる。
【００４１】
このため、比率算出部２５は、各解析点における高周波成分のコントラストＣＨθｉおよび低周波成分のコントラストＣＬθｉの比率、より具体的には、低周波成分のコントラストＣＬθｉに対する高周波成分のコントラストＣＨθｉの比率ＣＲθｉ（＝ＣＨθｉ／ＣＬθｉ）を算出する（ステップＳＴ５）。さらに、方向θｉ毎に、比率ＣＲθｉの平均値ＣＲｍθｉを算出し、最も低い平均値ＣＲｍθｉを体動の有無を判定するための指標値として算出する（ステップＳＴ６）。
【００４２】
なお、体動の有無を判定するための指標値の算出は、上記手法に限定されるものではなく、各解析点における８つの方向のすべてについての高周波成分のコントラストＣＨθｉおよび低周波成分のコントラストＣＬθｉの比率ＣＲθｉを算出し、方向θｉ毎にすべての解析点についての高周波成分のコントラストＣＨθｉおよび低周波成分のコントラストＣＬθｉの比率ＣＲθｉの平均値ＣＲｍθｉを算出し、この平均値のうちで最も低い平均値ＣＲｍＬを体動の有無を判定するための指標値として算出してもよい。
【００４３】
上述したように、撮影中に体動が生じた場合、取得される放射線画像Ｇ０においては、高周波成分のコントラストは低下するが、低周波成分のコントラストは殆ど低下しない。このため、体動が発生しなかった場合、指標値は１に近い値となるが、体動が発生した場合、指標値は１より小さく、かつ体動が大きいほど小さい値となる。
【００４４】
このため、判定部２６は、指標値をしきい値Ｔｈ２と比較し、指標値がしきい値Ｔｈ２以上か否かを判定する（ステップＳＴ７）。判定部２６は、指標値がしきい値Ｔｈ２以上の場合に体動なしと判定し（ステップＳＴ８）、しきい値Ｔｈ２未満の場合に体動ありと判定する（ステップＳＴ９）。なお、しきい値Ｔｈ２は１よりも小さいあらかじめ定められた値とする。
【００４５】
表示制御部２７は、判定部２６による判定結果を表示部４０に表示し（ステップＳＴ１０）、処理を終了する。図９は体動ありと判定された場合の判定結果の表示画面を示す図である。図９に示すように、判定結果の表示画面には、放射線画像Ｇ０および体動の方向および大きさを表す矢印が解析点に付与されている。なお、矢印を付与する解析点の数は、各解析点について算出した比率ＣＲθｉを小さい順に並べた場合における上位所定数（図９においては４つ）であってもよく、指標値がしきい値Ｔｈ２未満となった解析点のすべてであってもよい。また、矢印は比率ＣＲθｉの大きさに応じて色を変更してもよく、図９に示すように、比率ＣＲθｉの大きさに応じて線種を変更してもよい。なお、本実施形態においては、比率ＣＲθｉが小さいほど体動が大きいことから、矢印は比率ＣＲθｉが小さいほど体動が大きいことを表すものとなる。また、図１０に示すように比率ＣＲθｉの大きさに応じて矢印の長さを変更してもよい。ここで、図１０においては、比率が小さいほど矢印は長くなる。このように判定結果を表示することにより、放射線画像Ｇ０における体動が生じている部分および体動の大きさを確認することが容易となる。
【００４６】
なお、上述したように矢印を付与した場合であっても、放射線画像Ｇ０において体動を目視により確認したい場合がある。しかしながら、放射線画像の撮影室にある表示部４０の解像度は比較的小さく、また観察環境も比較的明るいことから、目視による体動の確認は容易ではない。このため、図１１に示すように、放射線画像Ｇ０上に所定サイズの関心領域Ｗを設定できるようにし、関心領域Ｗの画像を拡大して放射線画像Ｇ０に並べて表示するようにしてもよい。なお、関心領域Ｗの位置は、矢印を付与した解析点を含むようにあらかじめ設定しておいてもよいが、使用者が入力部３０を用いて放射線画像Ｇ０上の所望とする位置に移動できるようにしてもよい。
【００４７】
さらに、指標値を表示部４０に表示してもよい。この場合、単に数値のみを表示したのでは体動の大きさが判断しにくいため、今までに体動の有無を判定した結果を用いて、現在解析中の放射線画像Ｇ０における体動がどの程度あるかを分かりやすくしてもよい。例えば、図１２に示すように、横軸に指標値であるコントラストの比率（コントラスト比）の値を、縦軸に各コントラスト比の値において体動ありと判定した画像の数、および体動なしと判定した画像の数をプロットしたヒストグラムを表示し、そのヒストグラムにおける現在解析中の放射線画像Ｇ０の指標値を矢印にて指し示すようにしてもよい。なお、図１２において、ヒストグラムにおけるしきい値Ｔｈ２より左側の領域が体動あり、右側の領域が体動なしをそれぞれ表すものとなる。
【００４８】
一方、体動なしと判定された場合には、「体動はありません」というメッセージを表示したり、音声を出力したりすればよい。また、体動ありと判定された場合にも、「体動が見られます。撮影をし直してください」というメッセージを表示したり、音声を出力したりしてもよい。
【００４９】
このように、本実施形態においては、放射線画像Ｇ０から被写体領域Ａ０を抽出し、被写体領域Ａ０における低周波成分Ｌ３′がしきい値Ｔｈ１を超える画素を、体動検出のための解析点として抽出するようにしたものである。ここで、撮影中に体動が生じた場合、取得される放射線画像Ｇ０においては、高周波成分のコントラストは低下するが、低周波成分のコントラストは殆ど低下しない。このため、本実施形態によれば、体動に影響されることなく、体動を検出する処理を行うための解析点を放射線画像に精度良く設定することができ、その結果、精度良く体動を検出することができる。
【００５０】
なお、上記実施形態においては、多重解像度変換部２２において放射線画像Ｇ０を多重解像度変換して、放射線画像Ｇ０の高周波成分および低周波成分を取得しているが、解析点を含む所定範囲の領域の放射線画像Ｇ０をフーリエ変換して、解析点についてのパワースペクトルを求め、パワースペクトルにおけるあらかじめ定められた高周波数帯域の値および低周波数帯域の値を、解析点における高周波成分のコントラストおよび低周波成分のコントラストとして算出してもよい。この場合、高周波成分のコントラストと低周波成分のコントラストとの比率は、実質的には、パワースペクトルにおける高周波成分と低周波成分との比率となる。
【００５１】
また、上記実施形態においては、放射線検出器１１を用いて被写体の放射線画像を撮影する撮影装置１０において取得した放射線画像を用いて体動を判定しているが、特開平８−２６６５２９号公報、特開平９−２４０３９号公報等に示される放射線検出体としての蓄積性蛍光体シートに被写体の放射線画像情報を蓄積記録し、蓄積性蛍光体シートから光電的に読み取ることにより取得した放射線画像を用いて体動を判定するに際しても、本発明を適用できることはもちろんである。
【符号の説明】
【００５２】
１０    撮影装置
１１    放射線検出器
１２    放射線管球
１３    撮影制御部
１４    撮影台
２０    信号処理部
２１    前処理部
２２    多重解像度変換部
２３    解析点設定部
２４    コントラスト算出部
２５    比率算出部
２６    判定部
２７    表示制御部
３０    入力部
４０    表示部
５１    被写体領域抽出部
５２    解析点抽出部
【図１】
2014064857.tif 000003
【図２】
2014064857.tif 000004
【図３】
2014064857.tif 000005
【図４】
2014064857.tif 000006
【図５】
2014064857.tif 000007
【図６】
2014064857.tif 000008
【図７】
2014064857.tif 000009
【図８】
2014064857.tif 000010
【図９】
2014064857.tif 000011
【図１０】
2014064857.tif 000012
【図１１】
2014064857.tif 000013
【図１２】
2014064857.tif 000014