(57)【要約】  この発明の放射線画像処理装置は、分離手段と調整手段と除去手段とを備えている。被検体を介して撮影された放射線画像を用いて補正画像を生成するので、グリッド補正のためのデータを予め取得せずに１回の撮影で各画像を取得することができる。また、被検体を介して撮影された放射線画像をグリッド画像と非グリッド画像とに分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、補正画像を生成するので、１回の撮影で精度の高い補正処理を実現することができる。
2009044452.tif 000002
【特許請求の範囲】
【請求項１】
散乱放射線を除去するグリッドを介して撮影された放射線画像を処理する放射線画像処理装置であって、前記グリッドによるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに前記放射線画像を分離する分離手段と、前記非グリッド画像に基づいて前記グリッド画像の強度を実空間上で調整して、調整画像を生成する調整手段と、前記非グリッド画像から前記調整画像を減算して、グリッドの影響を除去した補正画像を生成する除去手段とを備えていることを特徴とする放射線画像処理装置。
【請求項２】
請求項１に記載の放射線画像処理装置において、前記分離手段は、前記グリッド画像と非グリッド画像とを空間周波数領域上で分離することを特徴とする放射線画像処理装置。
【請求項３】
請求項１または請求項２に記載の放射線画像処理装置において、前記分離手段は、前記放射線画像における周期的なパターンを抽出することにより画像成分を分離することを特徴とする放射線画像処理装置。
【請求項４】
請求項３に記載の放射線画像処理装置において、前記分離手段は、特定の周波数帯域を通過するバンドパスフィルタ処理により画像成分を分離することを特徴とする放射線画像処理装置。
【請求項５】
請求項３に記載の放射線画像処理装置において、前記分離手段は、前記放射線画像の実空間上の強度に関するプロファイルに関して三角関数と類似する成分を抽出することにより画像成分を分離することを特徴とする放射線画像処理装置。
【請求項６】
請求項１から請求項５のいずれかに記載の放射線画像処理装置において、前記非グリッド画像から前記調整画像を減算した画像に前記グリッド画像と同様の成分が含まれないように、前記調整手段はグリッド画像の強度を調整して、調整画像を生成することを特徴とする放射線画像処理装置。
【請求項７】
請求項６に記載の放射線画像処理装置において、前記非グリッド画像のグリッド位置の局所的な強度変化が単調でないときに、前記グリッド画像のグリッド位置の強度が、非グリッド画像のグリッド位置の局所的な強度変化に等しくなるように、前記調整手段はグリッド画像の強度を調整して、調整画像を生成することを特徴とする放射線画像処理装置。
【請求項８】
請求項６に記載の放射線画像処理装置において、前記グリッド画像のグリッド位置の強度が、前記非グリッド画像のグリッド位置の強度と非グリッド画像を平滑化した画像のグリッド位置の強度との差分に等しくなるように、前記調整手段はグリッド画像の強度を調整して、調整画像を生成することを特徴とする放射線画像処理装置。
【請求項９】
散乱放射線を除去するグリッドを介して撮影された放射線画像の処理をコンピュータに実行させるための放射線画像処理プログラムであって、前記グリッドによるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに前記放射線画像を分離する分離工程と、前記非グリッド画像に基づいて前記グリッド画像の強度を実空間上で調整して、調整画像を生成する調整工程と、前記非グリッド画像から前記調整画像を減算して、グリッドの影響を除去した補正画像を生成する除去工程とを備え、これらの工程での処理をコンピュータに実行させることを特徴とする放射線画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、散乱放射線を除去するグリッドを介して撮影された放射線画像を処理する放射線画像処理装置および放射線画像処理プログラムに関する。
【背景技術】
【０００２】
放射線画像処理装置は、被検体を透過する放射線エネルギーを「フラットパネル型放射線検出器（FPD: Flat Panel Detector）」と呼ばれるイメージセンサに蓄積し、電気信号に変換して可視化する。このとき、ＦＰＤは、被検体を透過した直線的な放射線エネルギーを蓄積するとともに、被検体を透過する際に散乱するエネルギーの低い放射線も同時に蓄積してしまう。したがって、撮影された放射線画像のコントラストが散乱放射線によって低下するなど、画質が低下してしまう。通常、この散乱放射線を除去するために、放射線吸収率の高い鉛などを使用した放射線吸収層と、放射線吸収率の低いアルミニウムや紙などを使用した放射線透過層とを交互に配置して構成されたグリッドを介して撮影している。【０００３】
しかし、グリッドを介して被検体を撮影する場合、放射線吸収層の影響による縞模様の像（グリッド像）が撮影画像である放射線画像に重畳してしまう。従来、このグリッド像を低減させるために、グリッドを揺動させながら撮影する揺動システムがある。しかし、この揺動システムを実現することで、装置が大型化する、コストがかかるなどの問題がある。【０００４】
そこで、グリッドに起因するバラツキを持ったグリッド像を除去するために、イメージセンサを構成する複数の画素の感度バラツキに基づいて、被検体の放射線画像を補正するための補正用画像を、被検体を介さずグリッドを介して予め撮影された画像データから作成し、その補正用画像を用いてグリッドに起因する画像成分を除去する手法がある（例えば、特許文献１参照）。その他にも、複数種類のグリッドや様々なイメージセンサに対応させるために、グリッドに関するグリッド周波数とイメージセンサのサンプリング周波数との関係に基づいてフィルタ処理を行う手法（例えば、特許文献２−４参照）などがある。
【特許文献１】特開２００５−３２４０６９号公報
【特許文献２】特許第２５０７６５９号公報（第１，３−５頁、第１−３，７図）
【特許文献３】特開２００４−３４４６７０号公報（第７頁、図７、表１）
【特許文献４】特開２００３−２６００５３号公報（第６頁、表１）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、上述した特許文献１の場合には、グリッド像除去のための基準となる補正用画像を、上述したように被検体を介さずに、かつグリッドを介して予め撮影する必要がある。実際には、使用するたびに被検体を介さない補正用画像を予め取得するのは煩雑であり、実用上現実的でない。また、精度の高い補正結果を得るためには、補正用画像を取得するとき、および実際に被検体を撮影するときの撮影環境や条件などを完全に一致させることが必要となるが、実際の使用状況で撮影環境や条件を完全に一致させることは実用上困難である場合が多い。【０００６】
また、上述した特許文献１の場合には、グリッド縞が画像全体にわたって定常であると仮定して、補正の対象となる対象画像から空間フィルタリングによってグリッド縞成分を抽出している。実際には、グリッド内の放射線吸収層の傾きは画像位置によって異なるので、画像中央部と端部とではグリッド縞の周波数帯域が異なる場合もあり、定常とは限らない。したがって、実際には周波数帯域が定常でない補正用画像を用いて、周波数帯域が画像全体にわたって定常であるとの仮定の下で補正を行うので、補正結果は精度の低いものとなる。また、補正用画像を予め生成することで、感度のバラツキを低減させることはできるが、補正用画像を作成するときの補正そのものの精度のバラツキに関しては考慮されていないので、必ずしも精度の高い補正を実現することができるとは限らない。【０００７】
この発明は、このような事情に鑑みてなされたものであって、グリッド補正のためのデータを予め取得せずに１回の撮影で精度の高い補正処理を実現することができる放射線画像処理装置および放射線画像処理プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、この発明の放射線画像処理装置は、散乱放射線を除去するグリッドを介して撮影された放射線画像を処理する放射線画像処理装置であって、前記グリッドによるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに前記放射線画像を分離する分離手段と、前記非グリッド画像に基づいて前記グリッド画像の強度を実空間上で調整して、調整画像を生成する調整手段と、前記非グリッド画像から前記調整画像を減算して、グリッドの影響を除去した補正画像を生成する除去手段とを備えていることを特徴とするものである。【０００９】
この発明の放射線画像処理装置によれば、グリッドによるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに、分離手段は放射線画像を分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整手段は調整して、調整画像を生成する。そして、非グリッド画像から調整画像を除去手段は減算して、グリッドの影響を除去した補正画像を生成する。このように被検体を介して撮影された放射線画像を用いて補正画像を生成するので、グリッド像除去のための基準となる補正用画像を、被検体を介さずに、かつグリッドを介して予め撮影する必要がなく、１回の撮影で各画像を取得することができ、撮影環境や条件などを完全に一致させる必要がなく制約もない。また、被検体を介して撮影された放射線画像をグリッド画像と非グリッド画像とに分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、補正画像を生成するので、調整によるグリッド画像の各画素での強度により各画素のバラツキを考慮しており、精度の高い補正処理を実現することができる。その結果、グリッド補正のためのデータを予め取得せずに１回の撮影で精度の高い補正処理を実現することができる。【００１０】
上述した分離手段の一例は、グリッド画像と非グリッド画像とを空間周波数領域上で分離することである。空間周波数領域上で分離することでグリッド像の成分とその他の成分とをより確実に分離することができる。もちろん、分離手段の他の一例は、グリッド画像と非グリッド画像とを実空間上で分離してもよい。【００１１】
また、分離手段は、放射線画像における周期的なパターンを抽出することにより画像成分を分離する。特に、グリッド像の成分は、ある特定の範囲の周波数帯域を持つ周期的なパターンであることが多いので、グリッド像に相当する特定の周波数帯域を通過するバンドパスフィルタ処理により画像成分を分離するのが好ましい。また、グリッド像が重畳した放射線画像が実空間上で周期的なパターンとなっているのを利用して、分離手段は、放射線画像の実空間上の強度に関するプロファイルに関して三角関数と類似する成分を抽出することにより画像成分を分離してもよい。【００１２】
また、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整する際には、非グリッド画像から調整画像を減算した画像にグリッド画像と同様の成分が含まれないように、調整手段はグリッド画像の強度を調整して、調整画像を生成する。具体的には、非グリッド画像のグリッド位置の局所的な強度変化が単調でないときに、グリッド画像のグリッド位置の強度が、非グリッド画像のグリッド位置の局所的な強度変化に等しくなるように、調整手段はグリッド画像の強度を調整して、調整画像を生成する。または、グリッド画像のグリッド位置の強度が、非グリッド画像のグリッド位置の強度と非グリッド画像を平滑化した画像のグリッド位置の強度との差分に等しくなるように、調整手段はグリッド画像の強度を調整して、調整画像を生成する。このように調整することで、非グリッド画像から調整画像を減算した画像（すなわち補正画像）にグリッド画像と同様の成分が含まれなくなり、より一層精度の高い補正を実現することができる。【００１３】
また、この発明の放射線画像処理プログラムは、散乱放射線を除去するグリッドを介して撮影された放射線画像の処理をコンピュータに実行させるための放射線画像処理プログラムであって、前記グリッドによるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに前記放射線画像を分離する分離工程と、前記非グリッド画像に基づいて前記グリッド画像の強度を実空間上で調整して、調整画像を生成する調整工程と、前記非グリッド画像から前記調整画像を減算して、グリッドの影響を除去した補正画像を生成する除去工程とを備え、これらの工程での処理をコンピュータに実行させることを特徴とするものである。【００１４】
この発明の放射線画像処理プログラムによれば、被検体を介して撮影された放射線画像を用いて補正画像を生成するので、グリッド補正のためのデータを予め取得せずに１回の撮影で各画像を取得することができる。また、被検体を介して撮影された放射線画像をグリッド画像と非グリッド画像とに分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、補正画像を生成するので、１回の撮影で精度の高い補正処理を実現することができる。
【発明の効果】
【００１５】
この発明に係る放射線画像処理装置および放射線画像処理プログラムによれば、被検体を介して撮影された放射線画像を用いて補正画像を生成するので、グリッド補正のためのデータを予め取得せずに１回の撮影で各画像を取得することができる。また、被検体を介して撮影された放射線画像をグリッド画像と非グリッド画像とに分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、補正画像を生成するので、１回の撮影で精度の高い補正処理を実現することができる。
【図面の簡単な説明】
【００１６】
【図１】実施例に係る放射線画像処理装置のブロック図である。
【図２】フラットパネル型放射線検出器（ＦＰＤ）の検出面の模式図である。
【図３】（ａ）〜（ｃ）は、分離の説明に供する各画像の概略図である。
【図４】一連の調整の流れを示したフローチャートである。
【図５】（ａ）〜（ｃ）は、調整の説明に供する各プロファイルの概略図である。
【符号の説明】
【００１７】
４  …  画像処理部
３１  …  グリッド
４１  …  分離部
４２  …  調整部
４３  …  除去部
【実施例】
【００１８】
以下、図面を参照してこの発明の実施例を説明する。図１は、実施例に係る放射線画像処理装置のブロック図であり、図２は、フラットパネル型放射線検出器（ＦＰＤ）の検出面の模式図であり、図３は、分離の説明に供する各画像の概略図であり、図４は、一連の調整の流れを示したフローチャートであり、図５は、調整の説明に供する各プロファイルの概略図である。【００１９】
本実施例に係る放射線画像処理装置は、図１に示すように、被検体Ｍを載置した天板１と、被検体Ｍに向けて放射線（例えばＸ線）を照射する放射線源２（例えばＸ線管）と、放射線源２から照射されて被検体Ｍを透過した放射線を検出するフラットパネル型放射線検出器（以下、「ＦＰＤ」と略記する）３と、ＦＰＤ３によって検出された放射線に基づいて画像処理を行う画像処理部４と、画像処理部４によって各種の画像処理された放射線画像を表示する表示部５とを備えている。表示部５はモニタやテレビジョンなどの表示手段で構成されている。また、ＦＰＤ３の検出面側にはグリッド３１を配設している。【００２０】
画像処理部４は、中央演算処理装置（ＣＰＵ）などで構成されている。なお、各種の画像処理を行うためのプログラム等をＲＯＭ（Read-only Memory）などに代表される記憶媒体に書き込んで記憶し、その記憶媒体からプログラム等を読み出して画像処理部４のＣＰＵが実行することでそのプログラムに応じた画像処理を行う。特に、画像処理部４の後述する分離部４１や調整部４２や除去部４３は、放射線画像の分離やグリッド画像の強度調整やグリッドの影響の除去に関するプログラムを実行することで、そのプログラムに応じた放射線画像の分離やグリッド画像の強度調整やグリッドの影響の除去をそれぞれ行う。放射線画像の分離やグリッド画像の強度調整やグリッドの影響の除去に関するプログラムは、この発明における放射線画像処理プログラムに相当する。【００２１】
画像処理部４は、グリッド３１によるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに放射線画像を分離する分離部４１と、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、調整画像を生成する調整部４２と、非グリッド画像から調整画像を減算して、グリッドの影響を除去した補正画像を生成する除去部４３とを備えている。分離部４１は、この発明における分離手段に相当し、調整部４２は、この発明における調整手段に相当し、除去部４３は、この発明における除去手段に相当する。【００２２】
ＦＰＤ３は、図２に示すように、その検出面には放射線に有感な複数の検出素子ｄを２次元マトリックス状に配列して構成されている。検出素子ｄは、被検体Ｍを透過した放射線を電子信号に変換して一旦蓄積して、その蓄積された電気信号を読み出すことで、放射線を検出する。各々の検出素子ｄでそれぞれ検出された電気信号を、電気信号に応じた画素値に変換して、検出素子ｄの位置にそれぞれ対応した画素にその画素値を割り当てることで放射線画像を出力して、画像処理部４の分離部４１（図１を参照）に放射線画像を送り込む。【００２３】
実際の撮影および各画像に関するデータの流れについて説明する。図示を省略する撮影ボタンを押下することで、放射線源２から放射線を発生させて、その放射線を被検体Ｍに向けて照射し、それに連動して撮影を開始する。発生した放射線は被検体Ｍを透過して、グリッド３１を介してＦＰＤ３に入射することで、ＦＰＤ３は放射線を検出して放射線画像を出力して撮影を行う。上述したように、ＦＰＤ３の検出面側にグリッド３１を配設しているので、放射線が被検体Ｍを透過する際に発生した散乱放射線をグリッド３１が除去する。したがって、ＦＰＤ３から出力された放射線画像にはグリッド３１の影響によるグリッド像が重畳している。【００２４】
このグリッド像が重畳した放射線画像を画像処理部４の分離部４１に送り込んで、グリッド像の成分を含むグリッド画像と、その他の成分を含む非グリッド画像とに、分離部４１は放射線画像を分離する。本実施例では、成分の分離には周波数解析を用いる。すなわち、グリッド像の成分は、ある特定の範囲の周波数帯域を持つ周期的なパターンであることが多いので、グリッド像に相当する特定の周波数帯域を通過するバンドパスフィルタ処理により成分を分離する。【００２５】
グリッド３１が定格（規格通りに設計されたグリッド）の場合には、既知である定格のグリッド３１のグリッド密度Ｄｇと同じく既知であるＦＰＤ３のイメージセンサのサンプリングピッチＰＳおよびサンプリング周波数ｆＳとを用いると、グリッドに関するグリッド周波数ｆｇをｆｇ＝ＰＳ・Ｄｇ・ｆＳの式で求めることができる。一方、ナイキスト周波数ｆｎはサンプリング周波数ｆＳの１／２、すなわちｆｎ＝ｆＳ／２で表される。サンプリング定理より、もしｆｇ＞ｆｎならば、グリッド周波数ｆｇがグリッド像の周波数であり、もしｆｇ≦ｆｎならば、２ｆｎ−ｆｇ（すなわちｆＳ−ｆｇ）がグリッド像の周波数である。そして、バンドパスフィルタ処理によってグリッド像の周波数のみを通過して、通過したグリッド像の周波数成分をグリッド画像として、残りの周波数成分を非グリッド画像として、放射線画像を分離する。【００２６】
グリッド３１が定格でなく（すなわち画像全体にわたって定常でなく）、上述したように画像中央部と端部とではグリッド縞の周波数帯域が異なる場合には、上述したサンプリング定理を用いずに、グリッド画像と非グリッド画像とを空間周波数領域上で分離する。具体的には、放射線画像に対してフーリエ変換してパワースペクトルを求め、パワースペクトルが極大値を示す範囲の周波数帯域を、グリッド像の周波数とみなすことで、グリッド画像と非グリッド画像とを分離する。【００２７】
より具体的に説明する。放射線画像の縦横をｘ方向、ｙ方向として、各画素値をｆ（ｘ、ｙ）とすると、空間周波数領域上でのフーリエ変換後のＦ（μ、ν）は、ｆ（ｘ、ｙ）を用いて下記（１）式のように表される。【００２８】
Ｆ（μ、ν）＝ΣΣｆ（ｘ、ｙ）・ｅ−ｊ２π（μｘ＋νｙ）／Ｎ    …（１）
ただし、上記（１）式中のｊは複素数であり、Σはｙ＝０〜Ｎ−１までの総和、ｘ＝０〜Ｎ−１までの総和である。【００２９】
上記（１）式によって放射線画像ｆ（ｘ、ｙ）に対してフーリエ変換された空間周波数領域上の画像Ｆ（μ、ν）は、図３（ａ）に示す通りである。つまり、図３（ａ）は、放射線画像ｆ（ｘ、ｙ）のパワースペクトルを可視化した画像Ｆ（μ、ν）である。図３（ａ）の画像を、予め設定した所定の大きさのブロックに分割し、各々のブロックにおいてブロック内のスペクトル強度の平均で置き換えたスペクトルデータを図３（ｂ）に示す。図３（ｂ）のスペクトルデータを走査し、極大値を示す位置のブロックを抽出した結果を図３（ｃ）に示す。図３（ｃ）では、図中の黒塗りで示した位置のブロックを、極大値として抽出している。【００３０】
（フーリエ変換後の）空間周波数領域上の画像Ｆ（μ、ν）は、グリッドの影響がなければ図３の画像中心から端部に向かって強度が単調に減少する。もし、単調に減少せずに極大値をとる箇所があれば、その箇所（ブロック）がグリッド像の周波数帯域であるとみなせる。したがって、図３の画像中心から走査を行う。後述するグリッド画像の強度調整における「単調」と同様に、本明細書中では、強度が極値をとらずに変化することを「単調」であるとし、強度が極値をとって変化することを「単調」でないとする。【００３１】
フーリエ変換された成分を、抽出されたブロック位置に相当する周波数帯域の成分（図３（ｃ）中の黒塗りのブロック）と、それ以外の位置に相当する周波数帯域の成分（図３（ｃ）中の白部分）とにバンドパスフィルタ処理により分離し、それぞれに対して逆フーリエ変換を行って画像化することにより、グリッド画像と非グリッド画像とを生成する。逆フーリエ変換後のｆ（ｘ、ｙ）は、Ｆ（μ、ν）を用いて下記（２）式のように表される。【００３２】
ｆ（ｘ、ｙ）＝ΣΣＦ（μ、ν）・ｅｊ２π（μｘ＋νｙ）／Ｎ      …（２）
ただし、上記（２）式中のΣはν＝０〜Ｎ−１までの総和、μ＝０〜Ｎ−１までの総和である。【００３３】
このように、グリッド３１が定格の場合には、サンプリング定理によってグリッド像の周波数を求めて、そのグリッド像の周波数を通過するバンドパスフィルタ処理により画像成分を分離し、グリッド３１が定格でない場合には、フーリエ変換による空間周波数領域上でパワースペクトルが極大値となるブロック位置に相当する周波数をグリッド像の周波数として、そのグリッド像の周波数を通過するバンドパスフィルタ処理により画像成分を分離することで、グリッド画像と非グリッド画像とを空間周波数領域上で分離する。【００３４】
なお、分離の手法については、上述した内容に限定されない。例えば、予め設定した所定の値よりも大きい周波数帯域と、それ以外の周波数帯域とに分離することで、グリッド画像と非グリッド画像とを分離してもよい。また、上述したようなグリッド画像と非グリッド画像とを空間周波数領域上で分離するのに限定されずに、グリッド像が重畳した放射線画像が実空間上で周期的なパターンとなっているのを利用して、放射線画像の実空間上の強度に関するプロファイルに関して三角関数と類似する成分（例えば正弦波、余弦波、それらの積・和、あるいは別の関数との積・和など）を抽出することにより画像成分を分離してもよい。すなわち、放射線画像の実空間上のプロファイルであって、三角関数と類似する成分を用いて、グリッド像に直交する（グリッド像の断面）プロファイルから類似する成分（極値）を抽出して、その抽出された成分をグリッド画像としてもよい。また、周波数解析についてもフーリエ変換に限定されず、ウェーブレット(Wavelet)変換やガボールフィルタなどのその他の周波数解析を用いてもよい。分離部４１による放射線画像の分離は、この発明における分離工程に相当する。【００３５】
分離部４１で生成されたグリッド画像と非グリッド画像とを調整部４２に送り込むとともに、非グリッド画像のみを除去部４３に送り込む。調整部４２は、非グリッド画像から後述する調整画像を減算した画像にグリッド画像と同様の成分が含まれないように、調整部４２はグリッド画像の強度を調整して、調整画像を生成する。本実施例では、図４のフローで調整部４２は調整を行う。【００３６】
（ステップＳ１）極値？
グリッド画像を、グリッド像と直交する方向（グリッド像の断面方向）に走査し、注目画素の強度が極値を示す画素であるか否かを調べる。具体的には、グリッド像に直交する（グリッド像の断面）強度に関するプロファイルについて、グリッド画像および非グリッド画像間で比較する。図５（ａ）は、グリッド画像のプロファイルであり、図５（ｂ）は、図５（ａ）と同じ位置における単調でないときの非グリッド画像のプロファイルであり、図５（ｃ）は、図５（ａ）と同じ位置における単調のときの非グリッド画像のプロファイルである。【００３７】
（ステップＳ２）次の画素
図５（ａ）に示すプロファイルで、注目画素の強度が極値を示す画素でなければ、グリッド像と直交する方向に沿って次の画素に進む。そして、ステップＳ１に戻って、極値を示す画素が見つかるまで、ステップＳ１、Ｓ２を繰り返して、グリッド像と直交する方向に沿って走査することになる。【００３８】
（ステップＳ３）単調？
図５（ａ）に示すプロファイルで、注目画素の強度が極値を示す画素であれば、図５（ａ）に示すように、強度が極大または極小を示している位置（すなわちグリッド画像のグリッド位置）において、非グリッド画像が局所的に単調に変化しているか否かを調べる。すなわち、非グリッド画像のグリッド位置の局所的な強度変化が単調であるか否かを調べる。非グリッド画像が局所的に単調に変化しているか否かについては、グリッド画像の強度の極大と極小とが反転する位置（図５（ａ）中のＡ，Ｃ）と、その内部の極値の位置（図５（ａ）中のＢ）とにそれぞれ対応する点の強度を比較することにより決定する。グリッド画像のＡ，Ｂ，Ｃ点に対応する非グリッド画像を、図５（ｂ）、図５（ｃ）に示すようにＡ´，Ｂ´，Ｃ´とすると、Ａ´の強度およびＣ´の強度の平均値とＢ´の強度とを比較して、所定の範囲を超える差異がある場合には、単調でないと判断する。【００３９】
図５（ｂ）に示すように、非グリッド画像のＢ´が極値となっている場合には、Ａ´の強度およびＣ´の強度の平均値とＢ´の強度との差分が、所定の範囲を超えているので、単調でないと判断する。このＢ´の強度は、グリッド画像のＢの強度に影響されて極値になっているとみなすことができる。すなわち、非グリッド画像は、分離部４１によって、グリッド像の成分を含むグリッド画像から分離されているが、グリッド画像のＢの強度に影響された成分が僅かながら非グリッド画像に含まれており、それによって非グリッド画像のグリッド位置の局所的な強度変化が単調でないとみなすことができる。単調でなければ次のステップＳ４に進む。【００４０】
逆に、図５（ｃ）に示すように、非グリッド画像のＢ´が極値となっていない場合には、強度が小さいＡ´の強度および強度が大きいＣ´の強度の平均値が、Ａ´およびＣ´の強度の相殺により、強度が中ほどのＢ´の強度と変わらなくなる。したがって、Ａ´の強度およびＣ´の強度の平均値とＢ´の強度との差分が、所定の範囲内に収まっており、単調であると判断する。単調であれば次のステップＳ４をスキップしてさらなる次のステップＳ５に進む。【００４１】
（ステップＳ４）調整
ステップＳ３で非グリッド画像のグリッド位置の局所的な強度変化が単調でないと判断された場合には、グリッド画像のグリッド位置（すなわちＢ）の強度が、非グリッド画像のグリッド位置（すなわちＢ´）の局所的な強度変化（Ａ´の強度およびＣ´の強度の平均値とＢ´の強度との差分）に等しくなるように、グリッド画像の強度を定数倍に調整する。【００４２】
（ステップＳ５）全ての画素？
ステップＳ１〜Ｓ４（あるいは単調であればステップＳ１〜Ｓ３）を１つの注目画素について行い、全ての画素について走査したか否かを調べる。全ての画素について走査した場合には、調整部４２によるグリッド画像の強度調整を全ての画素について行ったとして一連の調整を終了する。全ての画素について走査していない場合には、ステップＳ１に戻って、全ての画素について行うまで、ステップＳ１〜Ｓ４（あるいは単調であればステップＳ１〜Ｓ３）を繰り返して、グリッド像と直交する方向に沿って走査することになる。【００４３】
このように、ステップＳ１〜Ｓ５を行うことで一連の調整を終了する。なお、調整の手法については、上述した内容に限定されない。非グリッド画像のグリッド位置の局所的な強度変化が単調であるか否かについて調べたが、それ以外に、調整部４２において、非グリッド画像を平滑化した画像を作成し、グリッド画像のグリッド位置（Ｂ）の強度が、非グリッド画像のグリッド位置（Ｂ´）の強度と上述した平滑化した画像のグリッド位置の強度との差分に等しくなるように、調整部４２がグリッド画像の強度を調整して、調整画像を生成してもよい。調整部４２によるグリッド画像の強度調整は、この発明における調整工程に相当する。【００４４】
調整部４２で調整されて生成された調整画像を除去部４３に送り込むとともに、上述したように分離部４１で生成された非グリッド画像を除去部４３に送り込む。除去部４３は、非グリッド画像から調整画像を減算して、グリッドの影響を除去した補正画像を生成する。除去部４３によるグリッドの影響の除去は、この発明における除去工程に相当する。【００４５】
除去部４３で除去されて生成された補正画像を表示部５に送り込んで、表示部５は補正画像を表示出力する。また、除去部４３で除去されて生成された補正画像を、ＲＡＭ（Random-Access Memory）などに代表される記憶媒体に書き込んで記憶して、適宜必要に応じて読み出してもよいし、プリンタなどに代表される印刷手段に印刷出力してもよい。【００４６】
本実施例に係る放射線画像処理装置によれば、グリッド３１によるグリッド像の成分を含むグリッド画像とその他の成分を含む非グリッド画像とに、分離部４１は放射線画像を分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整部４２は調整して、調整画像を生成する。そして、非グリッド画像から調整画像を除去部４３は減算して、グリッドの影響を除去した補正画像を生成する。このように被検体Ｍを介して撮影された放射線画像を用いて補正画像を生成するので、グリッド像除去のための基準となる補正用画像を、被検体Ｍを介さずに、かつグリッド３１を介して予め撮影する必要がなく、１回の撮影で各画像を取得することができ、撮影環境や条件などを完全に一致させる必要がなく制約もない。また、被検体Ｍを介して撮影された放射線画像をグリッド画像と非グリッド画像とに分離し、非グリッド画像に基づいてグリッド画像の強度を実空間上で調整して、補正画像を生成するので、調整によるグリッド画像の各画素での強度により各画素のバラツキ（例えば、ＦＰＤ３のイメージセンサの感度のバラツキや補正のバラツキなど）を考慮しており、精度の高い補正処理を実現することができる。その結果、グリッド補正のためのデータを予め取得せずに１回の撮影で精度の高い補正処理を実現することができる。【００４７】
本実施例では、グリッド３１が定格でない場合には、分離部４１は、グリッド画像と非グリッド画像とを空間周波数領域上で分離している。空間周波数領域上で分離することでグリッド像の成分とその他の成分とをより確実に分離することができる。もちろん、放射線画像の実空間上の強度に関するプロファイルに関して三角関数と類似する成分を抽出する変形実施でも述べたように、グリッド画像と非グリッド画像とを実空間上で分離してもよい。【００４８】
分離部４１は、放射線画像における周期的なパターンを抽出することにより画像成分を分離する。上述したようにグリッド像の成分は、ある特定の範囲の周波数帯域を持つ周期的なパターンであることが多いので、グリッド像に相当する特定の周波数帯域を通過するバンドパスフィルタ処理により画像成分を分離するのが好ましい。【００４９】
非グリッド画像に基づいてグリッド画像の強度を実空間上で調整する際には、非グリッド画像から調整画像を減算した画像にグリッド画像と同様の成分が含まれないように、調整部４２はグリッド画像の強度を調整して、調整画像を生成する。具体的には、非グリッド画像のグリッド位置の局所的な強度変化が単調でないときに、グリッド画像のグリッド位置の強度が、非グリッド画像のグリッド位置の局所的な強度変化に等しくなるように、調整部４２はグリッド画像の強度を調整して、調整画像を生成する。このように調整することで、非グリッド画像から調整画像を減算した画像（すなわち補正画像）にグリッド画像と同様の成分が含まれなくなり、より一層精度の高い補正を実現することができる。【００５０】
この発明は、上記実施形態に限られることはなく、下記のように変形実施することができる。【００５１】
（１）上述した実施例では、放射線としてＸ線を例に採って説明したが、Ｘ線以外の放射線（例えばγ線など）に適用してもよい。【００５２】
（２）上述した実施例では、放射線画像処理装置は、医用等に用いられる、図１に示すような天板１に被検体を載置して撮影を行う構造であったが、これに限定されない。例えば、工業用等に用いられる非破壊検査装置のように被検体（この場合には検査の対象物が被検体）をベルト上に運搬させて撮影を行う構造であってもよいし、医用等に用いられるＸ線ＣＴ装置などのような構造であってもよい。また、撮影を行う装置を外部装置として別装置で構成し、分離手段と調整手段と除去手段とを単に備えた装置であってもよい。【００５３】
（３）グリッド像が現れやすいのは、画素にそれぞれ対応して配置された検出素子から構成されたフラットパネル型放射線検出器（ＦＰＤ）のようなイメージセンサを用いた場合である。したがって、この発明がＦＰＤに適用することで、精度の高い補正を実現することができて、特に有用であるが、イメージセンサはＦＰＤに限定されない。通常において用いられるイメージセンサであれば、この発明は適用することができる。
【図１】
2009044452.tif 000003
【図２】
2009044452.tif 000004
【図３】
2009044452.tif 000005
【図４】
2009044452.tif 000006
【図５】
2009044452.tif 000007
【国際調査報告】
2009044452.tif 000008
2009044452.tif 000009