(57)【要約】【課題】  画像が取得された条件を操作者が認識しやすいようにする画像表示技術を提供すること。【解決手段】  画像表示装置は、放射線を検出するセンサを用いて画像を取得し、その画像の行方向と列方向との少なくともいずれかにおけるデータについて、所定の状態で画像が取得された場合にピークが出現する特定の周波数の成分が所定の周波数の成分となるように、そのデータの周波数をシフトし、周波数がシフトされたデータに基づく画像を表示する。【選択図】  図１
2014150843.tif 000002
【特許請求の範囲】
【請求項１】
放射線を検出するセンサを用いて画像を取得する取得手段と、
前記画像の行方向と列方向との少なくともいずれかにおけるデータについて、所定の状態で前記画像が取得された場合にピークが出現する特定の周波数の成分が所定の周波数の成分となるように、前記データの周波数をシフトするシフト手段と、
周波数がシフトされた前記データに基づく画像を表示する表示手段と、
を有することを特徴とする画像表示装置。
【請求項２】
前記所定の状態は、前記センサにグリッドが装着された状態であり、
前記特定の周波数は、前記グリッドに起因する縞目の周波数である、
ことを特徴とする請求項１に記載の画像表示装置。
【請求項３】
前記所定の状態は、前記センサにグリッドが装着された状態であり、
前記シフト手段は、行方向と列方向のうち、前記グリッドに起因する縞目と直交する方向について、前記データの周波数をシフトする、
ことを特徴とする請求項１に記載の画像表示装置。
【請求項４】
前記取得手段で取得した画像から、前記グリッドに起因する縞目の有無と、前記特定の周波数と、前記グリッドに起因する縞目の方向とのうち、少なくともいずれかを解析する解析手段をさらに有する、
ことを特徴とする請求項２または３に記載の画像表示装置。
【請求項５】
前記シフト手段は、前記解析手段がグリッドに起因する縞目がないと判定した場合は、センサに装着されうるグリッドの候補に起因する縞目の周波数を、前記特定の周波数として、前記データの周波数をシフトする、
ことを特徴とする請求項４に記載の画像表示装置。
【請求項６】
前記解析手段で解析した結果を通知する通知手段をさらに有する、
ことを特徴とする請求項５に記載の画像表示装置。
【請求項７】
前記所定の周波数は、視覚的にコントラストに対する感度が高い周波数である、
ことを特徴とする請求項１から６のいずれか１項に記載の画像表示装置。
【請求項８】
前記取得手段が取得する画像は、被写体を含むＸ線画像または被写体を含まないキャリブレーションに用いる画像である、
ことを特徴とする請求項１から７のいずれか１項に記載の画像表示装置。
【請求項９】
取得手段が、放射線を検出するセンサを用いて画像を取得する取得工程と、
シフト手段が、前記画像の行方向と列方向との少なくともいずれかにおけるデータについて、所定の状態で前記画像が取得された場合にピークが出現する特定の周波数の成分が所定の周波数の成分となるように、前記データの周波数をシフトするシフト工程と、
表示手段が、周波数がシフトされた前記データに基づく画像を表示する表示工程と、
を有することを特徴とする画像表示方法。
【請求項１０】
請求項９に記載の画像表示方法の各工程をコンピュータに実行させるためのプログラム。
【請求項１１】
請求項９に記載の画像表示方法の各工程をコンピュータに実行させるためのプログラムを記憶したコンピュータが読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線画像の撮影時の条件を確認するための画像表示方法に関する。
【背景技術】
【０００２】
近年、医療用のＸ線撮影装置では、さまざまな方式を用いたデジタルＸ線撮影装置が普及してきている。例えば、蛍光体と大面積アモルファスシリコン（a-Si）センサとを密着させたＸ線検出器を使用し、光学系等を介さずにＸ線像を直接デジタル化する方式が実用化されている。また、アモルファスセレン（a-Se）等を使用してＸ線を直接光電変換して電子に変換し、該電子を大面積アモルファスシリコンセンサで検出する方式も同様に実用化されている。
【０００３】
上述したＸ線検出器を用いたデジタルＸ線撮影装置の場合、光電変換素子毎の感度のバラツキや、読出し回路内でのゲインのバラツキに起因するセンサの特性を補正（以後、ゲイン補正と呼ぶ）するために、一般的には、キャリブレーションが行われる。ここで、キャリブレーションとは、被検体を除外した状態でセンサ全面に略均一なＸ線を照射し撮影を行うことで、各画素のゲイン情報を表す画像（以後、この画像をキャリブレーションデータと呼ぶ）を取得するものである。また、ゲイン補正は、実際に被検体を撮影した画像（以後、この画像を撮影画像と呼ぶ）に対してキャリブレーションデータを除算（あるいは、対数変換後に減算）することで行われる。
【０００４】
しかしながら、上述したキャリブレーションにおいて、適切なキャリブレーションデータが取得されていない場合、ゲイン補正が正しく行えない可能性がある。例えば、散乱線除去用のグリッドが装着されない状態でのキャリブレーションデータを取得すべき場合に、グリッドが装着された状態で誤ってキャリブレーションデータを取得した場合が考えられる。この場合、キャリブレーションデータにグリッドに起因する縞目が写り込み、適切なゲイン補正が行えない可能性がある。散乱線除去用のグリッドが装着された状態でのキャリブレーションデータを取得すべき場合に、グリッドが装着されない状態で誤ってキャリブレーションデータを取得した場合も同様に、適切なゲイン補正が行えない可能性がある。すなわち、グリッドの装着の有無のような、キャリブレーションデータを取得する際の条件を、認識できることが重要となる。
【０００５】
そこで、この課題を解決する方法として例えば特許文献１の方法がある。この方法では、撮影したキャリブレーションデータから予め既知であるグリッドに起因する縞目に相当する周波数帯域のスペクトル値を閾値と比較することで、グリッドの写り込みを検査し、キャリブレーションデータが適切かを判断するものである。また、その他の方法として特許文献２の方法がある。この方法では、グリッドに起因する縞目に直交するラインをフーリエ変換して、得られたスペクトルにピークが有るか否かを判定し、ピークがあった場合は、グリッドが装着された状態であることを通知する。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第４７３１６９８号公報
【特許文献２】特開２０１０−２９２３８号公報
【特許文献３】特許第３４４５２５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
特許文献１及び特許文献２に記載の技術は、何れも撮影したキャリブレーションデータからグリッドの装着の有無を自動的に解析するものである。ただし、このような自動解析による方法は、誤った判定を行う場合があるため、判定結果と合わせて、キャリブレーションデータを表示し、操作者に目視確認させることが一般的に行われている。しかしながら、グリッドの密度、センサの画素ピッチ、表示サイズ等により、表示された画像上のグリッドに起因する縞目の周波数が変化する。したがって、グリッドの装着の有無のような、操作者が表示された画像で目視確認する場合、場合によってはグリッドが表示画面上で分り難く、画像が取得された条件の確認が困難となる場合がある。
【０００８】
本発明は上記課題に鑑みなされたものであり、画像が取得された条件を操作者が認識しやすいようにする画像表示技術を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、本発明による画像表示装置は、放射線を検出するセンサを用いて画像を取得する取得手段と、前記画像の行方向と列方向との少なくともいずれかにおけるデータについて、所定の状態で前記画像が取得された場合にピークが出現する特定の周波数の成分が所定の周波数の成分となるように、前記データの周波数をシフトするシフト手段と、周波数がシフトされた前記データに基づく画像を表示する表示手段と、を有する。
【発明の効果】
【００１０】
本発明によれば、画像が取得された条件を操作者が認識しやすくなる。
【図面の簡単な説明】
【００１１】
【図１】実施例１の放射線撮影装置の構成例を示す図。
【図２】実施例１の画像表示部における処理手順を示すフローチャート。
【図３】実施例２の放射線撮影装置の構成例を示す図。
【図４】実施例２の画像表示部における処理手順を示すフローチャート。
【図５】コントラスト感度を説明する図。
【図６】グリッドに起因する縞目に応じたラインデータの取得方法を説明する図。
【図７】グリッド周波数のシフト変換を説明する図。
【発明を実施するための形態】
【００１２】
以下、添付図面を参照して本発明の実施の形態を詳細に説明する。
【００１３】
＜＜実施形態１＞＞
（放射線撮影装置の構成）
本実施例に係る放射線撮影装置１００の構成例を図１に示す。放射線撮影装置１００は、画像表示機能を有し、例えば、放射線発生部１０１、放射線検出器１０４、データ収集部１０５、前処理部１０６、ＣＰＵ１０８、メインメモリ１０９、操作部１１０、画像処理部１１１、及び画像表示部１１２を有する。これらの機能部は、例えばＣＰＵバス１０７を介して、互いにデータ授受が可能なように接続される。また、画像表示部１１２は、放射線検出器１０４で撮影したキャリブレーションデータを表示する機能を有し、画像取得部１１３、周波数シフト部１１４、画像サイズ変換部１１５、及び表示部１１６を有する。これらの各構成部もまた、ＣＰＵバス１０７に接続される。
【００１４】
メインメモリ１０９は、ＣＰＵ１０８での処理に必要な各種のデータを記憶すると共に、ＣＰＵ１０８のワーキング・メモリとして機能する。ＣＰＵ１０８は、メインメモリ１０９を用いて、操作部１１０からの操作に従った装置全体の動作制御等を行う。これにより放射線撮影装置１００は、以下のように動作する。
【００１５】
まず、操作部１１０を介してユーザから撮影指示を取得すると、ＣＰＵ１０８は、この撮影指示をデータ収集部１０５伝達する。ＣＰＵ１０８は、撮影指示を受けると、放射線発生部１０１及び放射線検出器１０４を制御して放射線撮影を実行させる。
【００１６】
放射線撮影では、まず放射線発生部１０１が、被検体１０３に対して放射線ビーム１０２を照射する。放射線発生部１０１から照射された放射線ビーム１０２は、被検体１０３を減衰しながら透過して、放射線検出器１０４に到達する。そして、放射線検出器１０４は到達した放射線強度に応じた信号を出力する。なお、本実施例では、被検体１０３を人体とする。よって、放射線検出器１０４から出力される信号は人体を撮影したデータとなる。
【００１７】
データ収集部１０５は、放射線検出器１０４から出力された信号を所定のデジタル信号に変換して、画像データとして前処理部１０６に供給する。前処理部１０６は、データ収集部１０５から供給された画像データに対して、オフセット補正を行う。また、前処理部１０６は、後述するキャリブレーションによって取得したキャリブレーションデータ（各画素についてのゲイン情報）をメインメモリ１０９から読み出して、ゲイン補正を行う。前処理部１０６で前処理が行われた画像データは、ＣＰＵ１０８を介して、メインメモリ１０９、そして画像処理部１１１に順次転送される。画像処理部１１１は、画像データに対して各種画像処理、例えば、一般的な階調変換処理、鮮鋭化処理等を行い、処理後の画像を不図示のプリンタ等に出力し、一連の撮影動作が終了する。
【００１８】
（キャリブレーションの動作）
続いて、上述の構成を有する放射線撮影装置１００において実行される、本実施例に係るキャリブレーションの動作について、図２に示すフローチャートを用いて具体的に説明する。
【００１９】
まず、操作者は被検体１０３のない状態で、操作部１１０を介してキャリブレーションの撮影指示を行う。この撮影指示はＣＰＵ１０８によりデータ収集部１０５に伝えられる。ＣＰＵ１０８は、キャリブレーションの撮影指示を受けると、放射線発生部１０１及び放射線検出器１０４を制御して放射線撮影を行う。
【００２０】
データ収集部１０５は、放射線検出器１０４から出力された信号を所定のデジタル信号に変換して画像データとして前処理部１０６に供給する。前処理部１０６は、データ収集部１０５から供給された画像データに対して、オフセット補正を行う。この前処理部１０６で前処理が行われた画像データをキャリブレーションデータ（各画素についてのゲイン情報）として、メインメモリ１０９、画像取得部１１３に順次転送する（Ｓ２０１）。
【００２１】
なお、上の説明では、１回の撮影でキャリブレーションデータを取得しているが、複数回の放射線撮影を実行して複数枚のキャリブレーションデータを取得し、加算平均することで１枚のキャリブレーションデータをメインメモリ１０９に記憶してもよい。
【００２２】
次に、画像取得部１１３で取得したキャリブレーションデータに対し、周波数シフト部１１４において、Ｓ２０２〜Ｓ２０４の各ステップを実行する。これにより、キャリブレーションデータにおけるグリッドに起因する縞目に関する特定の周波数（以後、この特定の周波数をグリッド周波数と呼ぶ）の成分が所定の周波数の成分となるような、周波数変換が行われる。
【００２３】
なお、本実施例では、所定の周波数は、視覚的なコントラスト感度が高い周波数とする。ここで、図５（ａ）のように、モニタに表示された画像を観察した場合の視角１度あたりのサイクル数を周波数ｒ（cycle／degree）とした場合、視覚伝達関数ＶＴＦ（ｒ）はＤｏｏｌｅｙの近似式を用いると下記式にて表わされる。
2014150843.tif 000003
上記式の視覚伝達関数ＶＴＦ（ｒ）は、数値が大きい周波数ｒほど視覚的なコントラスト感度が高い周波数であることを示す関数である。ＶＴＦ（ｒ）は、図５（ｂ)のように周波数が約６(cycle/deree)で視覚的なコントラスト感度が最大となる。すなわち、図５（ｂ）によれば、視覚的に最も見やすい周波数は約６(cycle/degree)であることが分かる。したがって、所定の周波数ｒtを６（cycle／degree）に設定すれば、観察者にとって、最も視認性が良い周波数で表示を行うことができる。なお、視覚伝達関数ＶＴＦ（ｒ）の値が所定値以上であるような周波数の範囲から、任意に所定の周波数が選択されてもよい。すなわち、観察者の視認性が十分となる周波数であれば、どのような周波数が所定の周波数として選択されてもよい。
【００２４】
続いて、キャリブレーションデータにおけるグリッド周波数の成分が、上述した所定の周波数ｒt（cycle／degree）の成分となるように変換するためのシフト量を決定する（Ｓ２０２）。
【００２５】
まず、下記式にて所定の周波数ｒt（cycle／degree）を、モニタ表示上の１（mm）あたりのサイクル数で表わされる周波数ｈt（cycle／mm）の形に変換する。
2014150843.tif 000004
ただし、Ｌは図５（ａ）で示した観察距離（mm）である。一般的な値として例えばＬ＝１０００（mm）としてもよい。なお、各施設に適した観察距離が、Ｌの値として予め設定されてもよい。
【００２６】
次に、所定の周波数ｈt（cycle／mm）を、キャリブレーションデータの１（pixel）あたりのサイクル数で表わされる周波数ｆt（cycle／pixel）の形に変換する。ここで、表示するモニタの画素ピッチをＰm（mm）、モニタ表示する際の拡大率をＡ（倍）とすれば、下記式にて周波数ｆt（cycle／pixel）を求めることができる。
2014150843.tif 000005
ここで、拡大率Ａは、キャリブレーションデータの１（pixel）をモニタのＡ（pixel）に割り当てることを意味する。すなわち、拡大率Ａが１より大きい場合は拡大表示、拡大率Ａが１未満の場合は縮小表示、拡大率が１のときは、キャリブレーションデータの１（pixel）をモニタの１（pixel）に割り当てる、所謂ピクセル等倍表示を、それぞれ意味する。なお、拡大率Ａは、所定の周波数ｆt（cycle／pixel）がナイキスト周波数０．５（cycle／pixel）を超えない範囲で、使用するモニタのサイズおよび解像度に応じて任意に設定してもよい。
【００２７】
次に、キャリブレーションデータ上のグリッド周波数ｆg（cycle／pixel）を算出する。ここで、グリッドの密度をＤ（本／cm）とし、キャリブレーションデータの画素ピッチをＰc（mm）とすれば、グリッド周波数ｆg（cycle／pixel）は下記式にて求めることができる。
2014150843.tif 000006
ただし、ｎは、条件式
2014150843.tif 000007
を満たす整数である。ここで、使用されるグリッドの密度Ｄ（本／cm）は、各施設で既知であるため、その値を事前に設定すればよい。
【００２８】
なお、複数の異なる密度のグリッドを使用する場合は、上記式で求めた各々のグリッド周波数のうち、最大ｆgmax（cycle／pixel）および最小ｆgmin（cycle／pixel）となる値を用いて、グリッド周波数を算出してもよい。すなわち、
2014150843.tif 000008
のようにして、グリッド周波数ｆg（cycle／pixel）を算出してもよい。このようにすることで、使用されるすべてのグリッド周波数を、所定の周波数近傍の周波数へ変換することができる。
【００２９】
その後、求めたグリッド周波数ｆg（cycle／pixel）の成分を、上述の所定の周波数ｆt（cycle／pixel）の成分にシフトするためのシフト量Δｆを、
2014150843.tif 000009
のように算出する。
【００３０】
次に、キャリブレーションデータに対してグリッドに起因する縞目と直交する方向に１次元の離散フーリエ変換を行う（Ｓ２０３）。
【００３１】
具体的には、図６（ａ）のように、グリッドに起因する縞目が縦縞の場合は、キャリブレーションデータの各行を１ラインとしたラインデータｌ[ｎ]のすべてのラインに対して、１次元の離散フーリエ変換を行う。同様に、図６（ｂ）のように、グリッドに起因する縞目が横縞の場合は、キャリブレーションデータの各列を１ラインとしたラインデータｌ[ｎ]のすべてのラインに対して、１次元の離散フーリエ変換を行う。すなわち、行方向と列方向の少なくとも一方向についてのデータについて、１次元の離散フーリエ変換を行えばよい。なお、１次元の離散フーリエ変換は、
2014150843.tif 000010
のように表される。ただし、ｉは虚数を表し、Ｎはラインデータのサンプル数を表す。なお、グリッドに起因する縞目の方向は、通常、各施設で既知であり、その値を事前に設定すればよい。
【００３２】
次に、各ラインデータｌ[ｎ]から求めた離散フーリエ変換の結果Ｌ[ｋ]に対して、周波数のシフト変換を行う（Ｓ２０４）。
【００３３】
具体的には、Ｌ[ｋ]に対して、図７のように周波数スペクトル上のグリッド周波数ｆg（cycle／pixel）が所定の周波数ｆt（cycle／pixel）となるようにΔｆだけＬ[ｋ]の値をシフトさせたＬc[ｋ]を求める。ここで、Ｎが偶数でかつΔｆが負の場合は、
2014150843.tif 000011
のようにシフト変換を行う。
【００３４】
また、Ｎが偶数でかつΔｆが正の場合は、
2014150843.tif 000012
のようにシフト変換を行う。
【００３５】
また、Ｎが奇数でかつΔｆが負の場合は、
2014150843.tif 000013
のようにシフト変換を行う。
【００３６】
また、Ｎが奇数でかつΔｆが正の場合は、
2014150843.tif 000014
のようにシフト変換を行う。なお、
2014150843.tif 000015
は、ｘを超える最小の整数を示す。
【００３７】
以上のシフト変換をすべてのラインに対して実行し、処理を終了する。なお、周波数シフトについては、上述のように、フーリエ変換と逆フーリエ変換とを用いて実行してもよいし、例えば、データに対してΔｆに対応する正弦波を乗じることにより、周波数シフトを実行してもよい。すなわち、いずれかの手法によりデータの周波数シフトが行われればよく、フーリエ変換と逆フーリエ変換は行われなくてもよい。
【００３８】
次に、シフト変換を行ったすべてのラインのＬc[ｋ]に対して、
2014150843.tif 000016
のように、１次元の逆離散フーリエ変換を行う（Ｓ２０５）。ただし、ｉは虚数を表し、Ｎはラインデータのサンプル数を表す。
【００３９】
そして、求めたｌc[ｎ]のラインを、対応するキャリブレーションデータの各行または各列として、キャリブレーションデータを構成する。これにより、グリッド周波数ｆg（cycle／pixel）が、所定の周波数ｆc（cycle／pixel）に変換されたキャリブレーションデータを生成することができる。
【００４０】
グリッドが装着された状態で撮影されたキャリブレーションデータでは、上述の変換により、変換前のグリッド周波数ｆg（cycle／pixel）にあらわれる線スペクトルが所定の周波数ｆc（cycle／pixel）にシフトされる。その結果、変換後のキャリブレーションデータでは、画像上で周波数ｆc（cycle／pixel）の周期信号が出現する。なお、グリッドが装着されていない状態で撮影されたキャリブレーションデータでは、グリッド周波数ｆg（cycle／pixel）に線スペクトルが存在しない。このため、上述の変換を行っても、画像上に上述の周期信号は現れない。
【００４１】
なお、本実施例では、グリッドに起因する縞目の方向は各施設で既知であるとして説明したが、縞目の方向が未知の場合は、各行および各列の両方向に対して上述の処理を実行することで対応可能である。具体的には、キャリブレーションデータの各行を１ラインとしたラインデータに対して上述の処理を実行した後、処理後データに対してさらに各列を１ラインとしたラインデータに同様の処理を実行すれば両方向のグリッド周波数ｆg（cycle／pixel）を所定の周波数ｆc（cycle／pixel）に変換することが可能である。
【００４２】
次に、画像サイズ変換部１１５において、求めたキャリブレーションデータを拡大率Ａで拡大し表示サイズに変換する（Ｓ２０６）。ここで、拡大方法はどのようなものであってもよく、例えば、線形補間やＣｕｂｉｃ  Ｃｏｎｖｏｌｕｔｉｏｎ等の技術を用いてもよい。
【００４３】
次に、拡大したキャリブレーションデータを表示部１１６へ表示する（Ｓ２０７）。なお、表示部１１６で表示されたキャリブレーションデータは、グリッド周波数ｆg（cycle／pixel）が、モニタ表示上の所定の周波数ｈt（cycle／mm）に変換された形で表示されるものである。このとき、表示されるキャリブレーションデータに出現する周期信号の周波数ｆc（cycle／pixel）は操作者が視認しやすい周波数が選択されている。このため、操作者は、キャリブレーションデータの画像において周期信号を視認できるか否かに応じて、グリッドが装着されているか否かを確認することができる。
【００４４】
なお、キャリブレーションデータを、グリッドが装着されない状態で取得すべき場合は、グリッドが装着されてキャリブレーションデータが取得された場合については、グリッドを取り外して再度キャリブレーションを行う（Ｓ２０８）。一方、操作者が表示されたキャリブレーションデータを確認し、グリッドに起因する周期信号が確認できない場合はキャリブレーションを終了する。
【００４５】
また、キャリブレーションデータを、グリッドが装着された状態で取得すべき場合には、グリッドが装着されずにキャリブレーションデータが取得された場合については、グリッドを取り外して再度キャリブレーションを行うようにしてもよい。この場合、操作者が表示されたキャリブレーションデータを確認し、グリッドに起因する周期信号が確認できた場合はキャリブレーションを終了する。
【００４６】
以上、実施例１では、キャリブレーションデータを表示する際に、グリッドの密度、センサの画素ピッチ、表示サイズ等によらずグリッド周波数を常に一定の範囲の周波数の周期信号に変換する。これにより、グリッド周波数を視覚的にコントラスト感度の高い周波数に変換することで、人間の目で知覚し易く、操作者が目視確認しやすくすることができる。この結果、グリッドが装着された状態で撮影されたか否かを、操作者が目視確認しやすくなる。
【００４７】
＜＜実施形態２＞＞
本実施例に係る放射線撮影装置３００の構成例を図３に示す。また、本実施形態における画像表示部１１２の処理手順を図４に示す。この放射線撮影装置３００は、図１の放射線撮影装置１００に、グリッド解析部３０１と通知部３０２とを追加した構成を有する。なお、図３の放射線撮影装置３００において、図１の放射線撮影装置１００の機能ブロックと同様に動作する機能部については同じ符号を付し、その詳細な説明を省略する。また、図４に示すフローチャートにおいても同様に、図２に示したフローチャートと同様の処理を実行するステップは同じ符号を付して説明を省略する。
【００４８】
本実施例に係る放射線撮影装置３００では、まず、画像取得部１１３において、実施例１と同様にキャリブレーションデータを取得する（Ｓ２０１）。次に、グリッド解析部３０１において、グリッドに起因する縞目の有無、グリッド周波数、及び縞目の方向の少なくともいずれかをキャリブレーションデータから解析する（Ｓ４０１）。ここで、グリッド解析の具体的な方法は特に限定するものではないが、例えば、特許文献３に記載の方法を用いることができる。
【００４９】
特許文献３に記載の方法は、画像の複数行（または、複数列）のラインから平均的なパワースペクトルを求め、パワースペクトル上にグリッドに起因する縞目に起因する有意なピークが存在するか否かによってグリッドに起因する縞目の有無を判定する。そして、当該方法では、そのピークに対応する周波数をグリッド周波数とする。詳細については特許文献３を参照されたい。
【００５０】
本実施例では、同様の方法を用いてキャリブレーションデータの複数行および複数列のそれぞれのラインに対してグリッドに起因する縞目の有無を判定する。判定の結果、有意なピークが存在しない場合はグリッドに起因する縞目がないと判定し、有意なピークが存在する場合、縞目の方向と、そのピークに対応するグリッド周波数とを合わせて記憶する。なお、ここで記憶される縞目の方向とは、例えば、複数行のラインから求めたパワースペクトル上に有意なピークがある場合は縦縞、複数列のラインから求めたパワースペクトル上に有意なピークがある場合は横縞である。
【００５１】
次に、求めたグリッドに起因する縞目の方向とグリッド周波数を用いて、Ｓ２０２〜Ｓ２０７の各処理を実施例１と同様に実行する。なお、グリッドに起因する縞目がないと判定した場合は、放射線撮影装置３００は、各施設で使用されうるグリッドの候補についてのグリッド周波数を用いて、実施例１と同様の処理を行ってもよい。なお、装着されうるグリッドの候補については、複数存在する場合はその複数のグリッドの全てについて実施例１と同様の処理を実行し、全ての結果を表示するようにしてもよい。また、グリッドの候補が複数ある場合は、例えば、使用される確率が最も高い、確率が所定値以上、または確率が所定値以上の中からランダムに、１つ以上のグリッドを選択し、そのグリッドに関するグリッド周波数を用いて上述の処理を実行してもよい。また、グリッドに起因する縞目がないと判定した場合は、放射線撮影装置３００は、Ｓ２０２以降の処理を省略して処理を終了してもよい。
【００５２】
そして、通知部３０２において、グリッド解析部３０１で解析した結果をキャリブレーションデータと共にあわせて表示する（Ｓ４０２）。なお、図４においては、上述の実施例１のＳ２０２〜Ｓ２０７の処理後に解析結果の通知が実行されるように記載されているが、これに限られない。解析結果の通知は、例えば、Ｓ４０１の解析の直後に、またはＳ２０２〜Ｓ２０７の処理と並行して、実行されてもよい。
【００５３】
以上、実施例２では、キャリブレーションデータを目視確認する際に、グリッドを自動で解析した結果をセカンドオピニオンとして利用でき、操作者の確認の手間を低減することができる。
【００５４】
＜＜実施形態３＞＞
上述の実施例１及び実施例２では、キャリブレーションデータの周波数を変換することで、グリッドの装着を判定できる画像を生成する手法を説明した。本実施例３では、キャリブレーションデータの周波数変換は行わず、バンドパスフィルタを用いてグリッド装着の有無を判定可能とする手法について説明する。
【００５５】
本実施例３においては、キャリブレーションデータに対して、グリッドに起因する縞目と直交する方向においてバンドパスフィルタを適用し、グリッド周波数（及びその近傍）の要素を抜き出した画像を取得して表示する。これにより、グリッドが装着された状態で取得されたキャリブレーションデータであれば、縞目が強調された画像が出力される。一方、グリッドが装着されていない状態で取得されたキャリブレーションデータについても、グリッド周波数とその近傍についての要素は０ではないことが想定されるが、その要素はランダムに発生するため、縞目は現れない。したがって、操作者は、バンドパスフィルタを適用後のキャリブレーションデータにより、グリッドが装着されて撮影が行われたかを判定することができる。
【００５６】
なお、上述の実施例１〜３では、放射線撮影装置の機構（ハードウェア）によるグリッド装着の検出がないものとして説明したが、このようなグリッド装着を検出する機構が存在する場合にも上述の手法を適用することができる。この場合、放射線撮影装置は、例えば、機構と上述の手法を実行した結果とが共にグリッド装着を検出した場合に、グリッドが装着されていると判定してもよい。同様に、放射線撮影装置は、機構と上述の手法を実行した結果とが共にグリッドの不装着を検出した場合に、グリッドが装着されていないと判定してもよい。このようにすることにより、二重にグリッドの装着をチェックすることができ、グリッド装着／非装着の判断の誤りを防ぐことができる。
【００５７】
また、上述の実施例１〜３では、グリッドの装着があるかを特定するものとして説明したが、他の用途にも使用されうる。すなわち、所定の状態で放射線撮影を行うことにより取得されたキャリブレーションデータに対して、特定の周波数成分にピークが生じるような場合であれば、どのような場合であっても上述の手法を適用することができる。この場合、周波数シフトされたキャリブレーションデータを用いることにより、操作者は、その所定の状態で撮影されたかを判定することができる。
【００５８】
また、上述の実施例１〜３では、キャリブレーションデータについての画像表示処理について説明したが、これに限られない。すなわち、被写体が映り込んでいるＸ線画像（放射線画像）に対しても適用可能である。この場合、操作者は、例えば、被写体が映り込んでいない領域について、グリッドに起因する縞目が表示されているかを判断することで、Ｘ線画像の撮影時にグリッドが装着されていたかを判定することができる。この他、例えば、グリッド強調画像と、元のＸ線画像と、それらの差分画像とをＧＵＩを通じてボタンの押下等によって切り替えて表示可能としてもよい。特に、バンドパスフィルタでグリッド周波数（とその近傍）を抽出して得られる画像であれば、元のＸ線画像との差分画像を表示することで、撮影したＸ線画像に含まれるグリッドに起因する縞目を除去した画像を表示することができる。さらに、バンドパスフィルタ適用後の画像を用いて、Ｘ線画像にグリッドに起因する縞目を強調した画像を重ねて表示することにより、Ｘ線画像のどの部分にグリッドに起因する縞目の影響が出ているかを操作者に示すことが可能となる。
【００５９】
以上、本発明の好ましい実施例について説明したが、本発明はこれらの実施例に限定されず、その要旨の範囲内で種々の変形及び変更が可能である。
【００６０】
＜＜その他の実施形態＞＞
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【図１】
2014150843.tif 000017
【図２】
2014150843.tif 000018
【図３】
2014150843.tif 000019
【図４】
2014150843.tif 000020
【図５】
2014150843.tif 000021
【図６】
2014150843.tif 000022
【図７】
2014150843.tif 000023