(57)【要約】  ＦＰＤ型の放射線画像撮影装置で撮影された画像データに重畳されている筋状のムラを的確に除去可能な放射線画像処理装置を提供する。  放射線画像処理装置７０は、放射線画像撮影装置１で撮影された画像データＤが二次元状に配列された領域Ｑを各区分領域Ｓｋに分割する分割手段７１と、各区分領域Ｓｋ内の同一の走査線５の延在方向に並ぶ画像データＤの平均値Ｄａ(n)を算出する平均値算出手段７２と、撮影された対象物とその周囲との境界部分に生じる平均値Ｄａ(n)の差分を圧縮するエッジ圧縮手段７３と、差分が圧縮されたプロファイルに対して適応フィルタを施すフィルタ処理手段７４と、適応フィルタが施されたプロファイル上の各データｄｓ(n)に基づいて補正データＤｓを作成する補正データ作成手段７５とを備え、画像データＤから対応する補正データＤｓを減算処理してノイズが除去された画像データＤ＊を生成する。
2011145377.tif 000002
【特許請求の範囲】
【請求項１】
互いに交差するように配設された複数の走査線および複数の信号線により区画された各区画に各放射線検出素子が設けられ、前記各放射線検出素子が二次元状に配列された検出部を備える放射線画像撮影装置により撮影された前記各放射線検出素子ごとの各画像データに対して画像処理を行う放射線画像処理装置であって、
前記各画像データを前記二次元状に配列された複数の放射線検出素子にそれぞれ対応するように二次元状に配列した場合に、前記各画像データが二次元状に配列されて形成された領域を、少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分して、前記領域を各区分領域に分割する分割手段と、
前記各区分領域内の前記各画像データのうち、同一の前記走査線の延在方向に並ぶ前記各画像データの平均値を算出する処理を前記各走査線ごとに行い、前記各区分領域ごとに行う平均値算出手段と、
前記各区分領域ごとに算出した前記各走査線ごとの前記各画像データの平均値の前記信号線の延在方向のプロファイルにおいて、撮影された対象物とその周囲との境界部分に生じる前記各画像データの平均値の差分を圧縮した前記プロファイルを前記各区分領域ごとに作成するエッジ圧縮手段と、
前記差分が圧縮された前記プロファイルに対して適応フィルタを施す処理を前記各区分領域ごとに行うフィルタ処理手段と、
前記適応フィルタが施された前記プロファイル上の各データに基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成する補正データ作成手段と、
を備え、
前記二次元状に配列された前記各画像データから、対応する位置の前記補正データを減算処理して、前記各放射線検出素子ごとにノイズが除去された前記各画像データを生成することを特徴とする放射線画像処理装置。
【請求項２】
前記フィルタ処理手段は、前記差分が圧縮された前記プロファイルに対する前記適応フィルタとして、ウィナーフィルタを施すことを特徴とする請求の範囲第１項に記載の放射線画像処理装置。
【請求項３】
前記フィルタ処理手段は、前記ウィナーフィルタを、下記（１）式で表される形で前記差分が圧縮された前記プロファイルに適用し、算出された下記の値Ｄｆ(n)を、前記差分が圧縮された前記プロファイル上の前記平均値から減算した値を、前記適応フィルタが施された前記プロファイル上のデータとすることを特徴とする請求の範囲第２項に記載の放射線画像処理装置。
【数１】
2011145377.tif 000004
ただし、Ｄｆ(n)はフィルタ処理後のデータの値を表し、ＤＡaveは前記差分が圧縮された前記プロファイル上の前記各平均値の各区分領域ごとの平均値を表し、σｖ２は予め設定された分散の値を表し、σ２(n)は前記差分が圧縮された前記プロファイル上の前記平均値の分散を表す。
【請求項４】
前記フィルタ処理手段は、前記差分が圧縮された前記プロファイル上の前記各平均値にローパスフィルタ処理を施して、前記ウィナーフィルタを施すことを特徴とする請求の範囲第２項または第３項に記載の放射線画像処理装置。
【請求項５】
前記分割手段で前記領域を少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分する際の前記所定個数は、前記領域における前記走査線の延在方向の前記画像データの総数を自然数で除算した個数に設定されることを特徴とする請求の範囲第１項から第４項のいずれか一項に記載の放射線画像処理装置。
【請求項６】
互いに隣接する所定個数の前記各区分領域について、前記フィルタ処理手段により前記適応フィルタが施された前記プロファイル上の各データのうち、同一の前記走査線の延在方向に存在する前記各データの平均値を算出する処理を前記各走査線ごとに行う再平均化手段を備え、
前記補正データ作成手段は、前記再平均化手段が算出した前記各データの平均値に基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成することを特徴とする請求の範囲第１項から第５項のいずれか一項に記載の放射線画像処理装置。
【請求項７】
前記検出部上で、少なくとも前記複数の走査線が走査線の延在方向において分割されて構成されている前記放射線画像撮影装置により撮影された前記各画像データに対しては、前記複数の走査線が分割された前記検出部の各領域でそれぞれ得られた前記各画像データについてそれぞれ別々に前記分割手段から前記補正データ作成手段までの各処理を行い、それぞれ別々に二次元状に配列された前記各画像データから、対応する位置の前記補正データをそれぞれ別々に減算処理して、ノイズが除去された前記各画像データを生成することを特徴とする請求の範囲第１項から第６項のいずれか一項に記載の放射線画像処理装置。
【請求項８】
前記放射線画像撮影装置により取得された前記各放射線検出素子ごとの各オフセット補正値に対しても、前記分割手段から前記補正データ作成手段までの各処理を行い、二次元状に配列された前記各オフセット補正値から、対応する位置の前記補正データを減算処理して、ノイズが除去された前記各オフセット補正値を生成し、前記各画像データからノイズが除去された前記各オフセット補正値を減算することを特徴とする請求の範囲第１項から第７項のいずれか一項に記載の放射線画像処理装置。
【請求項９】
前記放射線画像撮影装置の前記各放射線検出素子ごとに基準となる前記各オフセット補正値を予め有しておき、前記放射線画像撮影装置により取得された前記各オフセット補正値から前記基準となる各オフセット補正値を減算した各差分に対して、前記分割手段から前記補正データ作成手段までの各処理を行い、二次元状に配列された前記各差分から、対応する位置の前記補正データを減算処理して、ノイズが除去された前記各差分を生成し、ノイズが除去された前記各差分を前記基準となる各オフセット補正値に加算して、前記放射線画像撮影装置の前記各放射線検出素子ごとのノイズが除去された前記各オフセット補正値を作成し、前記各画像データからノイズが除去された前記各オフセット補正値を減算することを特徴とする請求の範囲第１項から第７項のいずれか一項に記載の放射線画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線画像処理装置に係り、特に、放射線画像撮影装置で撮影された放射線画像から筋状のムラを除去する画像処理装置に関する。
【背景技術】
【０００２】
病気診断等を目的として、Ｘ線画像に代表される放射線を用いて撮影された放射線画像が広く用いられている。こうした医療用の放射線画像は、従来からスクリーンフィルムを用いて撮影されていたが、放射線画像のデジタル化を図るために輝尽性蛍光体シートを用いたＣＲ（Computed Radiography）装置が開発され、最近では、照射された放射線をフォトダイオード（光電変換素子）等の放射線検出素子で検出してデジタル画像データとして取得する放射線画像撮影装置が開発されている。【０００３】
このタイプの放射線画像撮影装置はＦＰＤ（Flat Panel Detector）として知られており、従来は、支持台等と一体的に形成された、いわゆる専用機として開発された（例えば特許文献１、２参照）。また、近年、放射線検出素子等をハウジングに収納して持ち運び可能とした可搬型の放射線画像撮影装置が開発され、実用化されている（例えば特許文献３、４参照）。【０００４】
放射線画像撮影装置としては、照射されたＸ線等の放射線の線量に応じて検出素子で電荷を発生させて電気信号に変換するいわゆる直接型の放射線画像撮影装置や、照射された放射線をシンチレータ等で可視光等の他の波長の電磁波に変換した後、変換され照射された電磁波のエネルギに応じてフォトダイオード等の光電変換素子で電荷を発生させて電気信号に変換するいわゆる間接型の放射線画像撮影装置が知られている。なお、本発明では、直接型の放射線画像撮影装置における検出素子や、間接型の放射線画像撮影装置における光電変換素子を、あわせて放射線検出素子という。【０００５】
ところで、特許文献５、６では、従来のＣＲ装置で得られた放射線画像に対して、放射線画像中に生じているムラを除去して、ムラがない放射線画像を生成するための放射線画像処理装置等が提案されている。この放射線画像処理装置等では、撮影された放射線画像の中からムラの成分を抽出し、抽出したムラの成分を元の放射線画像から除去することにより、ムラのない放射線画像を生成する。【０００６】
ＣＲ装置の場合、輝尽性蛍光体層を有するプレート等に被写体を介して放射線を照射して放射線画像撮影を行い、プレート等が放射線画像読取装置に搬送されて、放射線画像読取装置でプレートから放射線画像が読み出される。【０００７】
その際、例えば図２６に示すように、放射線画像読取装置１００では、レーザダイオード１０１から出射したレーザ光をコリメートレンズ１０２で平行光とし、その一部をビームスプリッタ１０３で反射してディテクタ１０４に送り、ディテクタ１０４によるモニタの結果がレーザドライバ回路１０５にフィードバックされてレーザダイオード１０１からのレーザ出力が調整される。【０００８】
また、シリンドリカルレンズからなる結像レンズ１０６で平行光の大半は屈折し、ポリゴンミラー１０７の鏡面上で反射し、反射光をｆθレンズ１０８、シリンドリカルミラー１０９を経てプレートＰの像面に照射する。その際、ポリゴンミラー１０９が回転することにより、プレートＰの像面に入射するレーザ光は励起光として像面上を主走査方向Ｘに移動し、読取線Ｚに沿って像面を走査する。【０００９】
プレートＰの像面に励起光であるレーザ光が入射すると、プレートＰの輝尽性蛍光体層のレーザ光が入射された位置に蓄積されている放射線エネルギが蛍光として放射される。そして、放射された蛍光は集光器１１０の導光管１１０ａを通って集光器１１０で集光され、集光された蛍光がフォトマルやフォトダイオード１１１に導かれる。フォトダイオード１１１等は、集光された蛍光に感応して出力電流を出力する。【００１０】
このようにして、プレートＰ上の各画素ごとの画像データが取得される。なお、主走査方向Ｘに直交する方向が副走査方向Ｙであり、プレートＰは副走査方向Ｙに少しずつ移動されながら、上記のように励起光であるレーザ光の走査を受ける。このようにして、プレートＰ上の各画素の画像データが二次元状に読み取られて、画像データが読み取られる。【００１１】
そして、例えば特許文献６に記載された放射線画像処理装置等では、図２７に示すように、読み取られた画像データを二次元状に並べ、例えば副走査方向Ｙに延在する短冊状の各領域Ｓａに分割し、短冊状の各領域Ｓａごとに、副走査方向Ｙに並ぶ各画素の画像データの分布を得る。そして、この画像データの分布に対してハイパスフィルタ処理やローパスフィルタ処理を施す等して、もともとの放射線画像ｐの中からムラの成分を抽出し、抽出したムラの成分を元の放射線画像ｐから除去することで、ムラのない放射線画像ｐを生成することが提案されている。【００１２】
ＣＲ装置で撮影された放射線画像ｐに生じるムラは、主に、放射線画像読取装置１００のポリゴンミラー１０７が回転する際に振動したりシリンドリカルミラー１０９が振動するなど、放射線画像読取装置１００の各部材の物理的な振動によって発生するものである。【００１３】
そこで、上記の特許文献６に記載の発明を適用する場合、放射線画像読取装置１００の各部材を調整して、それらの振動の周期が例えば高周波側で発生するように構成する等の措置がとられることが多い。放射線画像読取装置１００をこのように構成すると、放射線画像ｐに生じているムラの成分の周期を、例えば高周波側に寄せることができる。【００１４】
そして、人体内の臓器や骨等が撮影された画像データでは、画像データの分布はより低周波側（すなわち長周期側）に現れる。そのため、上記の画像データの分布に対して例えばハイパスフィルタ処理を施すことで、ムラの成分は残り、人体内の臓器や骨等が撮影された画像データの成分はカットオフされる。そのため、この場合は、ハイパスフィルタ処理により放射線画像ｐからムラの成分のみを抽出することができる。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】特許第３８９０１６３号公報
【特許文献２】特開平９−７３１４４号公報
【特許文献３】特開２００６−０５８１２４号公報
【特許文献４】特開平６−３４２０９９号公報
【特許文献５】特開平１０−２７６３３０号公報
【特許文献６】特開２００５−３０３４４８号公報
【発明の概要】
【発明が解決しようとする課題】
【００１６】
しかしながら、放射線画像撮影を上記のＦＰＤ型の放射線画像撮影装置で行う場合、放射線画像撮影装置を構成する各装置や部材を調整して、それらの振動の周期を例えば高周波側で発生するように構成する等の措置をとることは、事実上、困難である。【００１７】
また、放射線画像撮影装置では、上記の放射線画像読取装置の場合とは異なる原因で、撮影された放射線画像ｐ上にムラが生じることが分かっている。【００１８】
すなわち、放射線画像撮影装置では、通常、例えば後述する図５や図７に示すように、基板４上に複数の走査線５と複数の信号線６とが互いに交差するように設けられ、走査線５と信号線６により区画された各区画ｒに複数の放射線検出素子７が二次元状に配列されて構成される。【００１９】
そして、被写体を介して放射線画像撮影装置に放射線が照射されて放射線画像撮影が行われた後、走査駆動手段１５のゲートドライバ１５ｂから、例えば図２８に示すように、走査線５の各ラインＬ１〜Ｌｘにオン電圧が順次印加され、放射線の照射により各放射線検出素子７に蓄積された電荷がそれぞれ各信号線６を介して各読み出し回路１７で読み出されることで、各放射線検出素子７から画像データが順次読み出される。【００２０】
その際、走査駆動手段１５の電源回路１５ａ（図７参照）やバイアス電源１４等ではそれぞれ各走査線５やバイアス線９に印加する電圧にノイズが生じている。そして、走査線５のあるラインＬｎにオン電圧が印加されて、走査線５の当該ラインＬｎに接続されている各放射線検出素子７から読み出される各電荷には、読み出し処理の瞬間に電源回路１５ａやバイアス電源１４等で生じているノイズが電荷ノイズとして重畳される。【００２１】
そのため、走査線５の当該ラインＬｎに接続されている各放射線検出素子７から読み出される各電荷には、同じ電荷ノイズが重畳される。また、走査線５の次のラインＬn+1にオン電圧が印加されて、走査線５の当該ラインＬn+1に接続されている各放射線検出素子７から読み出される各電荷にも、同様にして、同じ電荷ノイズが重畳される。【００２２】
このように、各放射線検出素子７から読み出される各画像データｄには、走査線５の同一のラインＬに接続されている各放射線検出素子７については、同じ電荷ノイズが重畳される。しかし、図２８に示したように、走査線５の各ラインＬ１〜Ｌｘにはそれぞれ異なるタイミングでオン電圧が印加され、各タイミングで、電源回路１５ａやバイアス電源１４等で生じているノイズの大きさが異なる。【００２３】
そのため、各画像データそのものではなく、各画像データに重畳されるノイズ成分だけを見ると、同じ走査線５に接続されている各放射線検出素子７では同じノイズが重畳されているが、走査線５が異なると重畳されるノイズの大きさが異なる値になり、図２９に拡大して示すように、放射線画像ｐ上に、走査線５の延在方向に延びる筋状のムラが現れる。【００２４】
このように、ＦＰＤ型の放射線画像撮影装置で画像データを取得すると、それらの画像データに基づいて生成された放射線画像ｐ上に、走査線５の延在方向に延びる筋状のムラが現れるという特徴がある。なお、図２９では、走査線５の延在方向が矢印Ａで、それに直交する信号線６の延在方向が矢印Ｂで、それぞれ表されている。【００２５】
そこで、ＦＰＤ型の放射線画像撮影装置で取得された画像データに対する画像処理においては、このように特徴的に現れる筋状のムラを的確に除去することが望まれる。【００２６】
本発明は、上記の点を鑑みてなされたものであり、ＦＰＤ型の放射線画像撮影装置で撮影された画像データに重畳されている筋状のムラを的確に除去可能な放射線画像処理装置を提供することを目的とする。
【課題を解決するための手段】
【００２７】
前記の問題を解決するために、本発明の放射線画像処理装置は、
互いに交差するように配設された複数の走査線および複数の信号線により区画された各区画に各放射線検出素子が設けられ、前記各放射線検出素子が二次元状に配列された検出部を備える放射線画像撮影装置により撮影された前記各放射線検出素子ごとの各画像データに対して画像処理を行う放射線画像処理装置であって、
前記各画像データを前記二次元状に配列された複数の放射線検出素子にそれぞれ対応するように二次元状に配列した場合に、前記各画像データが二次元状に配列されて形成された領域を、少なくとも前記検出部上での前記走査線の延在方向において所定個数の前記画像データごとに区分して、前記領域を各区分領域に分割する分割手段と、
前記各区分領域内の前記各画像データのうち、同一の前記走査線の延在方向に並ぶ前記各画像データの平均値を算出する処理を前記各走査線ごとに行い、前記各区分領域ごとに行う平均値算出手段と、
前記各区分領域ごとに算出した前記各走査線ごとの前記各画像データの平均値の前記信号線の延在方向のプロファイルにおいて、撮影された対象物とその周囲との境界部分に生じる前記各画像データの平均値の差分を圧縮した前記プロファイルを前記各区分領域ごとに作成するエッジ圧縮手段と、
前記差分が圧縮された前記プロファイルに対して適応フィルタを施す処理を前記各区分領域ごとに行うフィルタ処理手段と、
前記適応フィルタが施された前記プロファイル上の各データに基づいて、前記領域における他の位置の各データを復元して、二次元状に配列された各補正データを作成する補正データ作成手段と、
を備え、
前記二次元状に配列された前記各画像データから、対応する位置の前記補正データを減算処理して、前記各放射線検出素子ごとにノイズが除去された前記各画像データを生成することを特徴とする。
【発明の効果】
【００２８】
本発明のような方式の放射線画像処理装置によれば、ＦＰＤ型の放射線画像撮影装置で撮影された画像データに対して、分割手段や平均値算出手段、エッジ圧縮手段、フィルタ処理手段、補正データ作成手段の各手段での各処理を行い、画像データから、ＦＰＤ型の放射線画像撮影装置特有の走査線の延在方向に延びる筋状のムラ（図２９参照）の成分のみを的確に抽出して、筋状のムラの成分のみからなる補正データを作成することが可能となる。【００２９】
そのため、画像データからそれぞれ対応する補正データを減算処理することによって、画像データから、それに重畳されている筋状のムラを的確に除去することが可能となり、筋状のムラが除去された画像データを的確に生成することが可能となる。そして、筋状のムラが除去された画像データに基づいて放射線画像を生成することで、放射線画像上から筋状のムラが的確に除去することが可能となり、筋状のムラのない放射線画像を生成することが可能となる。
【図面の簡単な説明】
【００３０】
【図１】放射線画像処理装置が処理する放射線画像を撮影するための放射線画像撮影システムの全体構成の例を示す図である。
【図２】放射線画像撮影装置の外観斜視図である。
【図３】図２の放射線画像撮影装置を反対側から見た外観斜視図である。
【図４】図２におけるＸ−Ｘ線に沿う断面図である。
【図５】放射線画像撮影装置の基板の構成を示す平面図である。
【図６】ＣＯＦやＰＣＢ基板等が取り付けられた基板を説明する側面図である。
【図７】放射線画像撮影装置の等価回路を表すブロック図である。
【図８】放射線画像撮影装置のコネクタとブッキー装置のコネクタとが接続された状態を表す外観斜視図である。
【図９】本実施形態における放射線画像処理装置での画像データに対する画像処理の手順を示すフローチャートである。
【図１０】放射線画像処理装置の機能的な構成を示すブロック図である。
【図１１】二次元状に配列された元画像データおよびその領域、走査線の延在方向、信号線の延在方向を説明する図である。
【図１２】元画像データの領域を分割して作成された各区分領域を説明する図である。
【図１３】区分領域内の同一の走査線の延在方向に並ぶ各元画像データの平均値を算出する処理を説明する図である。
【図１４】算出した平均値を信号線の延在方向に並べた場合の平均値の信号線の延在方向のプロファイルの例を示す図である。
【図１５】区分領域内の同一の信号線の延在方向に並ぶ各元画像データの平均値を算出する場合を説明する図である。
【図１６】エッジ圧縮処理が行われた後の平均値の信号線の延在方向のプロファイルの例を示すグラフである。
【図１７】フィルタ処理後のデータの信号線の延在方向のプロファイルの例を示すグラフである。
【図１８】区分領域ごとの補正データの信号線の延在方向のプロファイルの例を示すグラフである。
【図１９】フィルタ処理後のプロファイル上のデータを各区分領域の走査線方向の中央の位置に信号線の延在方向に並べた状態等を説明する図である。
【図２０】再平均化手段を設けた場合の放射線画像処理装置の機能的な構成を示すブロック図である。
【図２１】各走査線が分割されて構成されている放射線画像撮影装置の基板の構成を示す平面図である。
【図２２】図２１に示した放射線画像撮影装置で撮影され、前処理が施されただけの元画像データに基づいて生成された放射線画像を拡大して示す写真である。
【図２３】図２１に示した放射線画像撮影装置で撮影され、検出部の各領域を区別せずに処理された元画像データに基づいて生成された放射線画像を拡大して示す写真である。
【図２４】図２１に示した放射線画像撮影装置で撮影され、検出部の各領域で別々に処理された元画像データに基づいて生成された放射線画像を拡大して示す写真である。
【図２５】オフセット補正値から筋状のムラをより厳密に除去するための処理の手順を示すフローチャートである。
【図２６】放射線画像読取装置の構成を示す図である。
【図２７】二次元状に並べた画像データを短冊状の各領域に分割することを説明する図である。
【図２８】ＦＰＤ型の放射線画像撮影装置での画像データの読み出し処理で各走査線にオン電圧を印加するタイミングを表すタイミングチャートである。
【図２９】ＦＰＤ型の放射線画像撮影装置で撮影された画像データに基づいて生成された放射線画像上に現れる走査線の延在方向に延びる筋状のムラを表す図である。
【発明を実施するための形態】
【００３１】
以下、本発明に係る放射線画像処理装置の実施の形態について、図面を参照して説明する。ただし、本発明は以下の図示例のものに限定されるものではない。
【００３２】
以下では、放射線画像処理装置について説明する前に、放射線画像処理装置が処理する放射線画像を撮影するための放射線画像撮影システムについて説明する。本実施形態では、放射線画像撮影システム５０は、病院や医院等に設けられることが想定されており、放射線画像撮影装置１を用いて、患者の身体等を撮影する場合が想定されている。
【００３３】
なお、以下では、放射線画像撮影装置１が、シンチレータ等を備え、照射された放射線を可視光等の他の波長の電磁波に変換し光電変換素子で電気信号を得るいわゆる間接型の放射線画像撮影装置である場合について説明するが、直接型の放射線画像撮影装置であってもよい。また、以下では、放射線画像撮影装置が可搬型である場合について説明するが、支持台等と一体的に形成された放射線画像撮影装置であってもよい。
【００３４】
また、以下では、後述するように、放射線画像処理装置７０が、撮影室Ｒａと１対１に対応付けられて前室Ｒｂ内に設けられており、放射線画像処理装置７０が撮影室Ｒａ内等の各装置を制御するコンピュータ（いわゆるコンソール）の機能を兼用する場合について説明する。
【００３５】
しかし、例えば、撮影室Ｒａに対応付けられ、撮影室Ｒａ内等の各装置を制御するコンソールは別に設けておき、放射線画像処理装置７０を撮影室Ｒａや前室Ｒｂとは別の場所に設置し、放射線画像処理装置７０とコンソールとをＬＡＮ（Local Area Network）等のネットワークを介して接続するように構成することも可能である。
【００３６】
さらに、例えば、放射線画像処理装置７０と複数の撮影室Ｒａや、複数の放射線画像処理装置７０と複数の撮影室Ｒａとをネットワークを介して接続するように構成することも可能であり、放射線画像撮影システム５０は、放射線画像処理装置７０で放射線画像ｐを適切に画像処理できるように適宜構成される。
【００３７】
撮影室Ｒａは、患者の身体の一部である被写体（すなわち患者の撮影部位）に放射線を照射して放射線画像撮影を行う部屋であり、被写体に放射線を照射するための放射線照射装置の放射線発生装置５７の放射線源５２等が配置されている。また、撮影室Ｒａには、放射線画像撮影装置１を装填可能なブッキー装置５１が設けられている。
【００３８】
ここで、放射線画像撮影装置１について説明する。図２は、本実施形態に係る放射線画像撮影装置の外観斜視図であり、図３は、放射線画像撮影装置を反対側から見た外観斜視図である。また、図４は、図２のＸ−Ｘ線に沿う断面図である。放射線画像撮影装置１は、図２〜図４に示すように、筐体状のハウジング２内にシンチレータ３や基板４等で構成されるセンサパネルＳＰが収納されている。
【００３９】
図２や図３に示すように、本実施形態では、筐体２のうち、放射線入射面Ｒを有する中空の角筒状のハウジング本体部２Ａは、放射線を透過するカーボン板やプラスチック等の材料で形成されており、ハウジング本体部２Ａの両側の開口部を蓋部材２Ｂ、２Ｃで閉塞することで筐体２が形成されている。なお、筐体２をこのようないわゆるモノコック型として形成する代わりに、例えば、フレーム板とバック板とで形成された、いわゆる弁当箱型とすることも可能である。
【００４０】
図２に示すように、筐体２の一方側の蓋部材２Ｂには、電源スイッチ３７や選択スイッチ３８、コネクタ３９、バッテリ状態や放射線画像撮影装置１の稼働状態等を表示するＬＥＤ等で構成されたインジケータ４０等が配置されている。
【００４１】
また、図３に示すように、筐体２の反対側の蓋部材２Ｃには、画像データ等を放射線画像処理装置７０に無線で転送するための通信手段であるアンテナ装置４１が埋め込まれている。なお、画像データ等を放射線画像処理装置７０に有線方式で転送するように構成することも可能であり、その場合、例えば、前述したコネクタ３９にケーブル等を接続して送受信するように構成される。また、アンテナ装置４１を設ける場合には、アンテナ装置４１の筐体２上の配置場所や配置する個数は適宜決められる。
【００４２】
筐体２の内部には、図４に示すように、センサパネルＳＰの基板４の下方側に基台３１が配置され、基台３１には、電子部品３２等が配設されたＰＣＢ基板３３や緩衝部材３４等が取り付けられている。本実施形態では、基板４やシンチレータ３の放射線入射面Ｒ側には、それらを保護するためのガラス基板３５が配設されており、センサパネルＳＰと筐体２の側面との間には、それらがぶつかり合うことを防止するための緩衝材３６が設けられている。
【００４３】
シンチレータ３は、基板４の後述する検出部Ｐに貼り合わされるようになっている。本実施形態では、シンチレータ３は、例えば、蛍光体を主成分とし、放射線の入射を受けると３００〜８００ｎｍの波長の電磁波、すなわち可視光を中心とした電磁波に変換して出力するものが用いられる。
【００４４】
基板４は、本実施形態では、ガラス基板で構成されており、図５に示すように、基板４のシンチレータ３に対向する側の面４ａ上には、複数の走査線５と複数の信号線６とが互いに交差するように配設されている。基板４の面４ａ上の複数の走査線５と複数の信号線６により区画された各区画ｒには、放射線検出素子７がそれぞれ設けられている。
【００４５】
このように、走査線５と信号線６で区画された各区画ｒに二次元状に配列された複数の放射線検出素子７が設けられた区画ｒの全体、すなわち図５に一点鎖線で示される領域が検出部Ｐとされている。
【００４６】
本実施形態では、放射線検出素子７としてフォトダイオードが用いられているが、この他にも例えばフォトトランジスタ等を用いることも可能である。各放射線検出素子７は、後述する図７に示すように、スイッチ手段であるＴＦＴ８のソース電極８ｓに接続されている。また、ＴＦＴ８のドレイン電極８ｄは信号線６に接続されている。
【００４７】
そして、ＴＦＴ８は、後述する走査駆動手段１５から走査線５を介してゲート電極８ｇにオン電圧が印加されるとオン状態となり、ソース電極８ｓやドレイン電極８ｄを介して放射線検出素子７内に蓄積されている電荷を信号線６に放出させるようになっている。また、ＴＦＴ８は、接続された走査線５を介してゲート電極８ｇにオフ電圧が印加されるとオフ状態となり、放射線検出素子７から信号線６への電荷の放出を停止して、放射線検出素子７内に電荷を保持するようになっている。
【００４８】
本実施形態では、列状に配置された複数の放射線検出素子７にそれぞれバイアス線９が接続されており、図５に示すように、各バイアス線９は、基板４の検出部Ｐの外側の位置で１本の結線１０に結束されている。
【００４９】
また、各走査線５や各信号線６、バイアス線９の結線１０は、それぞれ基板４の端縁部付近に設けられた入出力端子（パッドともいう）１１に接続されている。各入出力端子１１には、図６に示すように、ＩＣ１２ａ等のチップがフィルム上に組み込まれたＣＯＦ（Chip On Film）１２が異方性導電接着フィルム（Anisotropic Conductive Film）や異方性導電ペースト（Anisotropic Conductive Paste）等の異方性導電性接着材料１３を介して接続されている。
【００５０】
また、ＣＯＦ１２は、基板４の裏面４ｂ側に引き回され、裏面４ｂ側で前述したＰＣＢ基板３３に接続されるようになっている。このようにして、放射線画像撮影装置１のセンサパネルＳＰの基板４部分が形成されている。なお、図６では、電子部品３２等の図示が省略されている。
【００５１】
ここで、図７を用いて放射線画像撮影装置１の回路構成について説明する。各放射線検出素子７の一方の電極にはそれぞれバイアス線９が接続されており、各バイアス線９は結線１０に結束されてバイアス電源１４に接続されている。バイアス電源１４は、結線１０および各バイアス線９を介して各放射線検出素子７の電極にそれぞれバイアス電圧（本実施形態では逆バイアス電圧）を印加するようになっている。
【００５２】
また、各放射線検出素子７の他方の電極はＴＦＴ８のソース電極８ｓ（図７中ではＳと表記されている。）に接続されており、各ＴＦＴ８のゲート電極８ｇ（図７中ではＧと表記されている。）は、走査駆動手段１５のゲートドライバ１５ｂから延びる走査線５の各ラインＬ１〜Ｌｘにそれぞれ接続されている。また、各ＴＦＴ８のドレイン電極８ｄ（図７中ではＤと表記されている。）は各信号線６にそれぞれ接続されている。
【００５３】
走査駆動手段１５は、ゲートドライバ１５ｂにオン電圧やオフ電圧を供給する電源回路１５ａと、走査線５の各ラインＬ１〜Ｌｘに印加する電圧をオン電圧とオフ電圧の間で切り替えるゲートドライバ１５ｂとを備えている。ゲートドライバ１５ｂは、前述したように、走査線５の各ラインＬ１〜Ｌｘを介してＴＦＴ８のゲート電極８ｇに印加する電圧をオン電圧とオフ電圧との間で切り替えて、各ＴＦＴ８のオン状態とオフ状態とを制御するようになっている。
【００５４】
そして、被写体を介して放射線画像撮影装置１に放射線が照射されて放射線画像撮影が行われた後の読み出し処理では、走査駆動手段１５のゲートドライバ１５ｂは、例えば図２８に示したように走査線５の各ラインＬ１〜Ｌｘにオン電圧を順次印加して、放射線の照射により各放射線検出素子７内に蓄積された電荷をそれぞれ各信号線６に放出させるようになっている。
【００５５】
各信号線６は、読み出しＩＣ１６内に形成された各読み出し回路１７にそれぞれ接続されている。読み出し回路１７は、増幅回路１８と、相関二重サンプリング（Correlated Double Sampling）回路１９と、アナログマルチプレクサ２１と、Ａ／Ｄ変換器２０とで構成されている。
【００５６】
放射線画像撮影の際に、被写体を介して放射線画像撮影装置１に放射線が照射されると、シンチレータ３で放射線が他の波長の電磁波に変換され、変換された電磁波がその直下の放射線検出素子７に照射される。そして、照射された放射線の線量（すなわち電磁波の光量）に応じて放射線検出素子７内で電荷が発生する。
【００５７】
そして、各放射線検出素子７からの画像データＤの読み出し処理において、上記のように走査駆動手段１５のゲートドライバ１５ｂから走査線５の所定のラインＬｎにオン電圧が印加されると、走査線５の当該ラインＬｎを介してそれに接続されている各ＴＦＴ８のゲート電極８ｇにオン電圧が印加されて各ＴＦＴ８がオン状態となり、オン状態となった各ＴＦＴ８と接続されている放射線検出素子７から各ＴＦＴ８を介して信号線６に電荷が放出される。
【００５８】
そして、放射線検出素子７から放出された電荷量に応じて増幅回路１８から電圧値が出力され、それを相関二重サンプリング回路１９で相関二重サンプリングしてアナログ値の画像データＤがマルチプレクサ２１に出力される。マルチプレクサ２１から順次出力された画像データＤは、Ａ／Ｄ変換器２０で順次デジタル値の画像データＤに変換され、記憶手段２３に出力されて順次保存されるようになっている。
【００５９】
制御手段２２は、ＣＰＵ（Central Processing Unit）やＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、入出力インターフェース等がバスに接続されたコンピュータや、ＦＰＧＡ（Field Programmable Gate Array）等により構成されている。専用の制御回路で構成されていてもよい。
【００６０】
制御手段２２は、放射線画像撮影装置１の走査駆動手段１５や読み出し回路１７等の各機能部の動作等を制御するようになっている。制御手段２２には、ＤＲＡＭ（Dynamic ＲＡＭ）等で構成される記憶手段２３や、放射線画像撮影装置１の各機能部に電力を供給するバッテリ２４が接続されている。また、制御手段２２には、前述したアンテナ装置４１が接続されている。
【００６１】
制御手段２２は、放射線画像撮影が行われると、上記のようにして、各放射線検出素子７から電荷すなわち画像データＤを読み出す読み出し処理を行うようになっている。また、本実施形態では、制御手段２２は、各放射線検出素子７から読み出して記憶手段２３に保存した各画像データＤに基づいて、各画像データＤを所定の割合で間引いてプレビュー用の間引きデータＤｔを作成するようになっている。
【００６２】
本実施形態では、制御手段２２は、放射線画像撮影が終了して各放射線検出素子７から画像データＤが読み出して記憶手段２３に保存すると、速やかに間引きデータＤｔを作成し、間引きデータＤｔを後述する放射線画像処理装置７０に送信するようになっている。
【００６３】
なお、本実施形態では、後述するように、放射線画像撮影装置１がブッキー装置５１に装填されている場合は、制御手段２２はブッキー装置５１を介して間引きデータＤｔや画像データＤ等を放射線画像処理装置７０に送信し、放射線画像撮影装置１がブッキー装置５１に装填されておらず単独で用いられている場合は、制御手段２２はアンテナ装置４１を介して間引きデータＤｔ等を送信するようになっている。
【００６４】
間引きデータＤｔは、例えば、二次元状に配列された各放射線検出素子７に対応して各画像データＤを配列した場合に３×３画素や４×４画素ごとに１画素分の画像データＤを抽出するようにして作成してもよく、或いは、走査線５の各ラインＬ１、Ｌ４、Ｌ７、…にそれぞれ接続された各放射線検出素子７からの画像データＤのように、走査線５の所定の間隔ごとの各ラインＬｎに接続された各放射線検出素子７からの画像データＤを抽出して作成するように構成することも可能である。
【００６５】
また、本実施形態では、制御手段２２は、間引きデータＤｔを放射線画像処理装置７０に送信した後、当該間引きデータＤｔの作成の基となった対応する画像データＤを放射線画像処理装置７０に自動的に送信するようになっている。
【００６６】
一方、本実施形態では、制御手段２２は、１回の放射線画像撮影が終了するごとに、或いは、一連の放射線画像撮影が終了した時点で、これらの放射線画像撮影で得られた画像データＤに重畳されている暗電荷等に起因するオフセット分を補正するためのオフセット補正値Ｏを得るために、いわゆるダーク読取処理を自動的に行うようになっている。なお、単独或いは一連の放射線画像撮影が開始される前に、ダーク読取処理を行うように構成することも可能である。
【００６７】
ダーク読取処理では、放射線画像撮影装置１の各ＴＦＴ８をオフ状態とし、放射線画像撮影装置１に放射線が照射されない状態で所定時間放置した後、上記の読み出し処理と同様にして、各放射線検出素子７から、それらに蓄積された暗電荷等がいわゆるダーク読取値ｄとして読み出され、記憶手段２３に保存される。
【００６８】
そして、制御手段２２は、読み出した各放射線検出素子７ごとのダーク読取値ｄをオフセット補正値Ｏとしたり、或いは、ダーク読取処理を複数回行って、各放射線検出素子７ごとに得られた複数のダーク読取値ｄを平均化する等して、オフセット補正値Ｏを算出して放射線画像処理装置７０に送信するようになっている。
【００６９】
次に、放射線画像撮影システム５０における各装置等について説明する。
【００７０】
図１に示すように、ブッキー装置５１は、本実施形態では、カセッテ保持部（カセッテホルダともいう。）５１ａに可搬型の放射線画像撮影装置１を装填して用いることができるようになっている。本実施形態では、図１に示すように、撮影室Ｒａには、ブッキー装置５１として立位撮影用と臥位撮影用のブッキー装置５１Ａ、５１Ｂが設置されているが、いずれか一方のみが設けられていてもよい。
【００７１】
本実施形態では、放射線画像撮影装置１は、ブッキー装置５１に装填して用いられる場合には、図８に示すように、放射線画像撮影装置１をブッキー装置５１に装填する前に、ブッキー装置５１から延びるケーブルの先端に設けられたコネクタ５１ｂを放射線画像撮影装置１のコネクタ３９に接続するようになっている。
【００７２】
なお、例えば、ブッキー装置５１のカセッテ保持部５１ａの内部にコネクタ５１ｂを設けておき、放射線画像撮影装置１をカセッテ保持部５１ａに装填すると放射線画像撮影装置１のコネクタ３９とコネクタ５１ｂが自動的に接続されるように構成することも可能である。
【００７３】
図１に示すように、撮影室Ｒａには、被写体に放射線を照射する放射線源５２が少なくとも１つ設けられており、本実施形態では、放射線源５２のうち、１つの放射線源５２Ａは、例えば撮影室Ｒａの天井からつり下げられて配置されており、その照射方向等を変えることにより、立位撮影用のブッキー装置５１Ａや臥位撮影用のブッキー装置５１Ｂに装填された放射線画像撮影装置１に対して放射線を照射することができるようになっている。
【００７４】
また、本実施形態では、ポータブルの放射線発生装置５２Ｂも設けられており、ポータブルの放射線発生装置５２Ｂは、撮影室Ｒａ内のいかなる場所にも持ち運びでき、任意の方向に放射線を照射できるようになっている。そして、ブッキー装置５１に装填されていない単独の状態の放射線画像撮影装置１を被写体である患者の身体の部分にあてがったり、臥位撮影用のブッキー装置５１Ｂや図示しないベッドと患者の身体との間に差し込んだりした状態で、ポータブルの放射線発生装置５２Ｂから放射線を照射することができるようになっている。
【００７５】
撮影室Ｒａ内には、放射線画像撮影装置１と放射線画像処理装置７０との通信を中継する無線アンテナ５３を備えた基地局５４が設けられている。なお、本実施形態では、ブッキー装置５１と基地局５４とはケーブル等で接続されており、放射線画像撮影装置１をブッキー装置５１に装填して用いる場合には、放射線画像撮影装置１から出力された間引きデータＤｔ等は、ブッキー装置５１や基地局５４等を介して放射線画像処理装置７０に送信される。
【００７６】
また、基地局５４には、クレードル５５が接続されている。クレードル５５は撮影室Ｒａと前室Ｒｂのいずれに設置されてもよく、撮影室Ｒａに設置される場合には、放射線発生装置５２から照射される放射線が到達しない位置、すなわち、例えば撮影室Ｒａのコーナーの位置等に設置される。
【００７７】
図１に示すように、前室（操作室等ともいう。）Ｒｂには、放射線源５２に対して放射線の照射開始等を指示するための曝射スイッチ５６等を備えた放射線発生装置の操作卓５７が設けられている。そして、放射線技師等の操作者は、曝射スイッチ５６を操作することにより放射線源５２から放射線画像撮影装置１に放射線を照射するようになっている。
【００７８】
また、前述したように、本実施形態では、放射線画像処理装置７０は、前室Ｒｂに設けられており、図示しないＣＰＵやＲＯＭ、ＲＡＭ、入出力インターフェース等がバスに接続されたコンピュータ等で構成されている。ＲＯＭには所定のプログラムが格納されており、放射線画像処理装置７０は、必要なプログラムを読み出してＲＡＭの作業領域に展開してプログラムに従って各種処理を実行するようになっている。
【００７９】
放射線画像処理装置７０には、ＬＣＤ（Liquid Crystal Display）等からなる表示部７０ａが設けられており、その他、キーボードやマウス等の図示しない入力手段等が接続されている。また、放射線画像処理装置７０には、ハードディスク等で構成された記憶手段７０ｂが接続されており、その他、図示を省略するが、他のコンピュータや、放射線画像処理装置７０から出力された画像データに基づいて放射線画像をフィルム等の画像記録媒体に記録して出力するイメージャ等の外部機器が、ＬＡＮ等を介して接続されている。
【００８０】
本実施形態では、放射線画像処理装置７０は、前述したように、放射線画像撮影装置１から間引きデータＤｔが送信されてくると、間引きデータＤｔに基づいて表示部７０ｂ上にプレビュー画像を表示するようになっている。このプレビュー画像の表示は、放射線技師等の操作者がプレビュー画像を見て放射線画像撮影が正常に行われたか否かを確認して再撮影の要否等を判断するために行われるものであり、操作者は、再撮影が必要であると判断すると、当該放射線画像撮影を再度行う。
【００８１】
また、操作者が再撮影は不要であると判断して、放射線画像処理装置７０に対して当該放射線画像撮影について画像処理を行う旨の入力が行われると、放射線画像処理装置７０は、当該放射線画像撮影で得られた画像データＤに対して画像処理を行うようになっている。
【００８２】
本実施形態では、放射線画像処理装置７０は、図９に示すフローチャートに従って画像データＤに対する画像処理を行うようになっている。
【００８３】
そして、放射線画像処理装置７０は、図１０に示すように、機能的に、少なくとも分割手段７１、平均値算出手段７２、エッジ圧縮手段７３、フィルタ処理手段７４および補正データ作成手段７５を備えて構成されている。
【００８４】
そして、図９のムラ除去処理（ステップＳ３）では、これらの各手段７１〜７５で、後述するように前処理が行われた画像データＤ（すなわち後述する元画像データＤori）から、前述した筋状のムラ（図２９参照）の成分のみからなる補正データＤｓを作成して、各元画像データＤoriからそれぞれ対応する補正データＤｓを減算処理することによって、筋状のムラが除去された各画像データＤ＊を生成するようになっている。
【００８５】
以下、放射線画像処理装置７０における画像データＤに対する画像処理について詳しく説明する。また、本実施形態に係る放射線画像処理装置７０の作用について説明する。
【００８６】
なお、本発明に特徴的な処理は、上記のようにムラ除去処理（ステップＳ３）であり、その他の処理は、放射線画像撮影で得られた画像データＤに対する通常の画像処理で行われる公知の処理であるため、簡単に説明したり、或いは説明を省略する。
【００８７】
放射線画像処理装置７０は、上記のように、放射線画像撮影装置１から各放射線検出素子７ごとの画像データＤ（すなわちいわゆるrawデータ）とオフセット補正値Ｏが送信されてくると、それらの各データを記憶手段７０ｂに保存して、まず、ゲイン・オフセット補正処理（ステップＳ１）を行う。
【００８８】
ゲイン・オフセット補正処理（ステップＳ１）では、放射線画像処理装置７０は、記憶手段７０ｂに記憶されている当該放射線画像撮影装置１の各放射線検出素子７ごとのゲイン補正値Ｇを読み出して、下記（１）式に従って、画像データＤに基づいて元画像データＤoriを各放射線検出素子７ごとに算出する。
Ｄori＝Ｇ×（Ｄ−Ｏ）    …（１）
【００８９】
なお、元画像データＤoriとは、下記のムラ除去処理（ステップＳ３）で、補正データＤｓを作成する基となる画像データという意味であり、作成された補正データが減算処理されて筋状のムラが除去された各画像データＤ＊を生成する基となる画像データという意味である。
【００９０】
次に、放射線画像処理装置７０は、欠陥補正処理（ステップＳ２）では、記憶手段７０ｂに記憶されている欠陥画素マップを参照して、欠陥画素の補正を行う。ここで、欠陥画素とは、異常な画像データＤが読み出される放射線検出素子７、或いはその放射線検出素子７から読み出された異常な画像データＤのことをいう。
【００９１】
欠陥補正処理（ステップＳ２）では、放射線画像処理装置７０は、欠陥画素マップを参照して、欠陥画素としてマップ上に示されている異常な放射線検出素子７の元画像データＤoriを無効として廃棄し、当該異常な放射線検出素子７に検出部Ｐ上で隣接する正常な放射線検出素子７の元画像データＤoriで線形補間する等して、異常な放射線検出素子７の元画像データＤoriとする。
【００９２】
放射線画像処理装置７０は、続いて、本発明に特徴的なムラ除去処理（ステップＳ３）を行う。なお、このムラ除去処理（ステップＳ３）は、上記のようにゲイン・オフセット補正処理（ステップＳ１）や欠陥補正処理（ステップＳ２）を行った後に行うことが好ましい。
【００９３】
また、後述する対数変換処理等（ステップＳ４）が行われた後に行うと、上記（１）式に従って算出された元画像データＤoriのうち、値が相対的に小さい元画像データＤoriは値の差異が増幅され、値が相対的に大きい元画像データＤoriは値の差異が圧縮されるため、筋状のムラを補正しづらくなる。そのため、このムラ除去処理（ステップＳ３）は、対数変換処理等（ステップＳ４）を行う前に行うことが好ましい。
【００９４】
ムラ除去処理（ステップＳ３）は、分割手段７１から補正データ作成手段７５までの各手段で、前述したように元画像データＤoriから、前述した筋状のムラ（図２９参照）の成分のみからなる補正データＤｓを作成する。そして、各元画像データＤoriからそれぞれ対応する補正データＤｓを減算処理することによって、筋状のムラが除去された各画像データＤ＊を生成する処理である。
【００９５】
まず、補正データＤｓを作成するための前提として、図１１に示すように、各元画像データＤoriを、検出部Ｐ（図５や図７参照）上で二次元状に配列された複数の放射線検出素子７にそれぞれ対応するように、二次元状に配列した場合を考える。
【００９６】
なお、検出部Ｐ上で、ｍ本目の信号線６とｎ本目の走査線５（すなわち走査線５のラインＬｎ）との交点の位置にある放射線検出素子７（ｍ，ｎ）に対応する各元画像データＤoriをＤori（ｍ，ｎ）と表す。また、図１１中では、元画像データＤori（ｍ，ｎ）が単に（ｍ，ｎ）で表されている。
【００９７】
また、元画像データＤori（ｍ，ｎ）を図１１に示したように配列させた場合、検出部Ｐ上での走査線５の延在方向（以下、走査線方向という。）は図中の矢印Ａで表され、信号線６の延在方向（以下、信号線方向という。）は図中の矢印Ｂで表される。
【００９８】
放射線画像処理装置７０の分割手段７１は、図１１に示したように各元画像データＤori（ｍ，ｎ）が二次元状に配列されて形成された領域Ｑを、図１２に示すように、少なくとも走査線方向Ａにおいて所定個数の元画像データＤoriごとに区分して、領域Ｑを各区分領域Ｓｋに分割するようになっている。
【００９９】
区分領域Ｓｋの走査線方向Ａにおける元画像データＤoriの所定個数は、適宜の個数に設定され、例えば２０個等に設定される。また、図１２に示したように、区分領域Ｓｋの信号線方向Ｂは、領域Ｑの信号線方向Ｂの全ての元画像データＤoriとすることも可能であり、区分領域Ｓｋを信号線方向Ｂにおいても所定個数の元画像データＤoriごとに区分して、区分領域Ｓｋを信号線方向Ｂにも分割するように構成してもよい。
【０１００】
放射線画像処理装置７０の平均値算出手段７２は、図１３に示すように、分割手段７１が領域Ｑを分割して作成した各区分領域Ｓｋ内の各元画像データＤori（ｍ，ｎ）のうち、同一の走査線５のラインＬｎの延在方向に並ぶ各元画像データＤori（ｍ，ｎ）の平均値Ｄａを算出する。そして、この処理を走査線５の各ラインＬ１〜Ｌｘ（図７参照）ごとに行う。
【０１０１】
このように、各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａを算出すると、各平均値Ｄａは、走査線５の各ラインＬ１〜Ｌｘごとの値となり、ライン番号ｎの変数になる。そのため、以下、走査線５のラインＬｎについての平均値ＤａをＤａ(n)と表す。また、それらの平均値Ｄａ(n)を信号線方向Ｂに並べると、図１４に示すように、各平均値Ｄａ(n)の信号線方向Ｂのプロファイルが得られる。なお、各平均値Ｄａ(n)の算出は、各区分領域Ｓｋごとに行われる。
【０１０２】
図２９に示したように、本実施形態のようなＦＰＤ型の放射線画像撮影装置１では、走査線方向Ａに延びる筋状のムラが特徴的に現れる。そのような状況において、前述したように、各元画像データＤori（ｍ，ｎ）の平均値を算出する際に、仮に、図１５に示すように、各元画像データＤori（ｍ，ｎ）の信号線方向Ｂの平均値Ｄｂ(m)を算出するように構成してしまうと、平均値Ｄｂ(m)の算出の際に走査線方向Ａに延びる筋状のムラの成分が互いに打ち消しあってしまうため、元画像データＤoriから筋状のムラの成分を抽出することができない。
【０１０３】
それに対し、図１３に示したように、各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａを算出すると、平均値Ｄａ(n)の算出の際に、走査線方向Ａに延びる筋状のムラの成分が互いに打ち消しあうことはなく、算出された平均値Ｄａ(n)に、筋状のムラの成分がそのまま残存するため、元画像データＤoriから筋状のムラの成分を抽出することが可能となる。
【０１０４】
なお、ゲイン・オフセット処理（図９のステップＳ１）や欠陥補正処理（ステップＳ２）を施されただけの元画像データＤoriには、筋状のムラの成分だけでなく、種々のノイズ成分が重畳されているが、上記のように、各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａを算出することにより、筋状のムラの成分は残存するが、その他の種々のノイズ成分は互いに打ち消しあうため低減される。
【０１０５】
そのため、各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａを算出することにより、平均値Ｄａは、ノイズ成分として、主に筋状のムラの成分を含むデータとなる。
【０１０６】
また、このようにして、各平均値Ｄａ(n)の信号線方向Ｂのプロファイルを得ると、例えば患者の人体内の臓器や骨等が撮影された部分とそれ以外の部分とで、平均値Ｄａ(n)が異なる値になる。なお、図１４では、臓器や骨等が撮影されている部分の平均値Ｄａ(n)がその周囲の部分の平均値Ｄａ(n)よりも大きな値とされる場合が示されている。
【０１０７】
次に、放射線画像処理装置７０のエッジ圧縮手段７３は、上記のように、各区分領域Ｓｋごとに算出した平均値Ｄａ(n)の信号線方向Ｂのプロファイルにおいて、撮影された患者の人体内の臓器や骨等の対象物とその周囲との境界部分に生じる上記の平均値Ｄａ(n)の差分を圧縮したプロファイルを、各区分領域Ｓｋごとに作成するようになっている。
【０１０８】
具体的には、エッジ圧縮手段７３は、例えば、上記のように、各区分領域Ｓｋごとに算出した平均値Ｄａ(n)の信号線方向Ｂのプロファイルを、例えばｎが小さい方から順に見ていった場合に、ある平均値Ｄａ(n)からその直前の平均値Ｄａ(n-1)を減算して算出した差分Ｄａ(n)−Ｄａ(n-1)の絶対値｜Ｄａ(n)−Ｄａ(n-1)｜が、予め設定された閾値ΔＤａth以上に大きい場合には、当該平均値Ｄａ(n)を含むそれ以降の各平均値Ｄａ(n)、Ｄａ(n+1)、Ｄａ(n+2)、…から前記差分を減算して、差分を圧縮したプロファイルを作成するように構成される。なお、以下、この差分を圧縮した後の平均値を、平均値Ｄａ(n)と区別するために、平均値ＤＡ(n)と表す。
【０１０９】
すなわち、区分領域Ｓｋの各平均値Ｄａ(n)の信号線方向Ｂのプロファイルが、例えば図１４に示したように算出された場合、そのプロファイルを例えばｎが小さい方から順に見ていった場合、ｎ＝ａやｎ＝ｂの所で平均値Ｄａ(n)が大きく変化している。そして、例えばｎ＝ａの所では、平均値Ｄａ(a)からその直前の平均値Ｄａ(a-1)を減算して算出した差分Ｄａ(a)−Ｄａ(a-1)の絶対値｜Ｄａ(a)−Ｄａ(a-1)｜が、予め設定された閾値ΔＤａth以上に大きくなる。
【０１１０】
そのため、エッジ圧縮手段７３は、平均値Ｄａ(a)とそれ以降の各平均値Ｄａ(a+1)、Ｄａ(a+2)、…から、それぞれ差分Ｄａ(a)−Ｄａ(a-1)を減算して、圧縮後の平均値ＤＡ(a)を算出する。
【０１１１】
この場合、例えばｎ＝ａでの平均値Ｄａ(a)から差分Ｄａ(a)−Ｄａ(a-1)を減算すると、圧縮後の平均値ＤＡ(a)は、
ＤＡ(a)＝Ｄａ(a)−（Ｄａ(a)−Ｄａ(a-1)）
∴ＤＡ(n)＝Ｄａ(a-1)    …（２）
となり、直前のｎ＝ａ−１における平均値Ｄａ(a-1)に等しくなる。
【０１１２】
また、平均値Ｄａ(a+1)、Ｄａ(a+2)、…についても、
ＤＡ(a+1)＝Ｄａ(a+1)−（Ｄａ(a)−Ｄａ(a-1)）  …（３）
ＤＡ(a+2)＝Ｄａ(a+2)−（Ｄａ(a)−Ｄａ(a-1)）  …（４）
等の各減算処理が行われて、圧縮後の平均値ＤＡ(a+1)、ＤＡ(a+2)、…が算出される。また、平均値Ｄａ(b)を含むそれ以降の各平均値Ｄａ(b)、Ｄａ(b+1)、Ｄａ(b+2)、…についても同様の処理が行われて、圧縮後の平均値ＤＡ(b)、ＤＡ(b+1)、ＤＡ(b+2)、…が算出される。
【０１１３】
エッジ圧縮手段７３でこのように処理が行われることで、例えば図１４に示した区分領域Ｓｋの平均値Ｄａ(n)のプロファイルから、図１６に示すような差分を圧縮後の平均値ＤＡ(n)のプロファイルが作成される。エッジ圧縮手段７３は、この処理を、各区分領域Ｓｋごとに行うようになっている。
【０１１４】
上記のように患者の人体内の臓器や骨等を撮影した画像では、臓器や骨等が撮影された部分と、それ以外の部分が撮影された部分との境界、すなわちいわゆるエッジの部分で、隣接する平均値Ｄａ(n)同士の値が大きく変化する。一方、筋状のムラの成分でも隣接する平均値Ｄａ(n)同士で異なる値になるが、これらの臓器や骨等のエッジ部分における隣接する平均値Ｄａ(n)同士の値の差分は、筋状のムラの成分の場合の差分より格段に大きくなる。
【０１１５】
そのため、上記の閾値ΔＤａthとして、臓器や骨等のエッジ部分における隣接する平均値Ｄａ(n)同士の値の差分と、筋状のムラの成分における隣接する平均値Ｄａ(n)同士の値の差分とを的確に区別できるような適切な値が設定される。
【０１１６】
そして、このように設定することで、筋状のムラの成分については、上記のエッジ圧縮処理が行われずに圧縮後の平均値ＤＡ(n)のプロファイル中に残り、臓器や骨等のエッジ部分については、上記のエッジ圧縮処理が行われて、圧縮後の平均値ＤＡ(n)のプロファイル中から除去することが可能となる。
【０１１７】
このようにして、上記のようなエッジ圧縮処理を行うことにより、圧縮後の平均値ＤＡ(n)のプロファイルに基づいて作成される後述する補正データＤｓ中に、臓器や骨等のエッジ部分が残らないようにすることが可能となり、補正データＤｓを、筋状のムラの成分のみのデータとすることが可能となる。
【０１１８】
なお、放射線画像処理装置７０のエッジ圧縮手段７３におけるエッジ圧縮処理は、上記のように、ある平均値Ｄａ(n)からその直前の平均値Ｄａ(n-1)を減算して算出した差分Ｄａ(n)−Ｄａ(n-1)に基づいて行う代わりに、例えば、平均値Ｄａ(n)のプロファイルにおける低周波成分を圧縮する、いわゆるダイナミック圧縮処理によって行うように構成することも可能である。
【０１１９】
また、例えば、注目するｎを含む例えばｎ−５〜ｎ＋５等の所定の範囲における平均値Ｄａ(n)の平均値を算出し、当該ｎにおける平均値Ｄａ(n)と当該所定範囲における平均値Ｄａ(n)の平均値との差分を算出し、その差分の絶対値が予め設定された閾値を越えた場合に、上記のようにエッジを圧縮するように構成することも可能である。
【０１２０】
このように、エッジ圧縮手段７３におけるエッジ圧縮処理は、上記のように、各区分領域Ｓｋごとに算出した平均値Ｄａ(n)の信号線方向Ｂのプロファイルにおいて、撮影された患者の人体内の臓器や骨等の対象物とその周囲との境界部分に生じる上記の平均値Ｄａ(n)の差分を圧縮したプロファイルを適切に作成することができるものであればよく、特定の手法に限定されない。
【０１２１】
次に、放射線画像処理装置７０のフィルタ処理手段７４は、上記のように作成された圧縮後の平均値ＤＡ(n)のプロファイルに対して、適応フィルタを施す処理を、各区分領域Ｓｋごとに行うようになっている。
【０１２２】
本実施形態では、上記の適応フィルタとしてウィナーフィルタが用いられるようになっている。ウィナーフィルタは、よく知られているように、ノイズが重畳されていない真のデータとノイズが重畳されているデータとの二乗誤差を最小にする適応フィルタであり、本実施形態の場合に適用する場合には、ウィナーフィルタは、例えば下記（５）式で表されるような形になる。
【０１２３】
【数１】
2011145377.tif 000003
【０１２４】
ここで、Ｄｆ(n)はフィルタ処理後のデータの値を表し、ＤＡaveは各区分領域Ｓｋごとの圧縮後の平均値ＤＡ(n)の平均値を表し、σｖ２は予め設定された分散の値を表し、σ２(n)はＤＡ(n)の分散を表す。すなわち、ＤＡ(n)の分散σ２(n)は、ｎを含む微小範囲（例えばｎ−５〜ｎ＋５の範囲）の（ＤＡ(n)−ＤＡave）２の総和Σ（ＤＡ(n)−ＤＡave）２である。
【０１２５】
そして、ＤＡ(n)の分散σ２(n)が設定された分散σｖ２に比べて非常に大きい場合（すなわちσ２(n)＞＞σｖ２の場合）には、上記（５）式は、
Ｄｆ(n)≒ＤＡ(n)    …（６）
となり、フィルタ処理後のデータＤｆ(n)は、元の圧縮後の平均値ＤＡ(n)に近い値になる。
【０１２６】
また、ＤＡ(n)の分散σ２(n)が設定された分散σｖ２とほぼ等しい場合（すなわちσ２(n)≒σｖ２の場合）には、上記（５）式は、
Ｄｆ(n)≒ＤＡave    …（７）
となり、フィルタ処理後のデータＤｆ(n)は、圧縮後の平均値ＤＡ(n)の平均値ＤＡaveに近い値になる。
【０１２７】
上記（５）式に従って算出されるフィルタ処理後のデータＤｆ(n)は、上記のように各区分領域Ｓｋにおける筋状のムラ等のノイズが重畳されていない真のデータであり、図１７に示すようなプロファイルになる。
【０１２８】
一方、このムラ除去処理（ステップＳ３）で得たいのは、逆に、前述した筋状のムラ（図２９参照）の成分のみからなる補正データＤｓである。
【０１２９】
そこで、放射線画像処理装置７０のフィルタ処理手段７４は、下記（８）式に従って、図１６に示した圧縮後の平均値ＤＡ(n)から、上記のようにして算出した筋状のムラ等のノイズが重畳されていない真のデータであるフィルタ処理後のデータＤｆ(n)（図１７参照）を減算することで、当該区分領域Ｓｋにおける補正データＤｓ作成の基となるデータｄｓ(n)を算出するようになっている。
ｄｓ(n)＝ＤＡ(n)−Ｄｆ(n)    …（８）
以下、このｄｓ(n)を区分領域Ｓｋごとの補正データという。
【０１３０】
そして、この区分領域Ｓｋごとの補正データｄｓ(n)を筋状のムラの成分のみのデータとするために、上記の（５）式において予め設定する分散σｖ２の値を的確に設定することが必要となる。そして、分散σｖ２の値を、筋状のムラに対応する分散の値と同程度の値に予め設定しておくと、この目的を達成することが可能となる。
【０１３１】
すなわち、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定した場合、上記（５）式によれば、ある圧縮後の平均値ＤＡ(n)の分散σ２(n)が分散σｖ２（すなわち筋状のムラに対応する分散）に比べて非常に大きい場合（すなわちσ２(n)＞＞σｖ２の場合）には、上記（６）式に示したようにＤｆ(n)≒ＤＡ(n)となり、これを上記（８）式に代入すると、区分領域Ｓｋごとの補正データｄｓ(n)はほぼ０になる。
【０１３２】
つまり、前述したようにｎを含む例えばｎ−５〜ｎ＋５の微小範囲について算出されるＤＡ(n)の分散σ２(n)が、分散σｖ２すなわち筋状のムラに対応する分散に比べて非常に大きい場合とは、例えばその微小範囲において圧縮後の平均値ＤＡ(n)が圧縮後の平均値ＤＡ(n)の平均値ＤＡaveよりも全体的に大きい、或いは全体的に小さい場合を表している。これは、その微小範囲では、筋状のムラの成分よりも、例えば患者の人体内の臓器や骨等に由来する長周期の成分が支配的であることを意味する。
【０１３３】
そして、上記のように区分領域Ｓｋごとの補正データｄｓ(n)を筋状のムラの成分のみのデータとするためには、そのような筋状のムラの成分以外の長周期の成分が区分領域Ｓｋごとの補正データｄｓ(n)中に残らないことが望まれる。そして、上記のように、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定すればこの微小範囲では区分領域Ｓｋごとの補正データｄｓ(n)はほぼ０になり、筋状のムラの成分以外の長周期の成分が区分領域Ｓｋごとの補正データｄｓ(n)中に残らない。このように、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定することが望ましい。
【０１３４】
また、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定すると、上記（５）式によれば、ある圧縮後の平均値ＤＡ(n)の分散σ２(n)が分散σｖ２（すなわち筋状のムラに対応する分散）とほぼ等しい場合（すなわちσ２(n)≒σｖ２の場合）には、上記（７）式に示したようにＤｆ(n)≒ＤＡaveとなり、これを上記（８）式に代入すると、区分領域Ｓｋごとの補正データｄｓ(n)はほぼＤＡ(n)−ＤＡaveになる。
【０１３５】
つまり、この場合には、区分領域Ｓｋごとの補正データｄｓ(n)が０ではなく、圧縮後の平均値ＤＡ(n)に含まれる筋状のムラの成分とほぼ等しい値となる。このように、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定すれば、上記の微小範囲が筋状のムラを反映した部分であれば、区分領域Ｓｋごとの補正データｄｓ(n)が筋状のムラの成分とほぼ等しい値となる。
【０１３６】
上記のように区分領域Ｓｋごとの補正データｄｓ(n)を筋状のムラの成分のみのデータとするためには、筋状のムラの成分が反映されたプロファイルの部分では筋状のムラの成分が的確に区分領域Ｓｋごとの補正データｄｓ(n)中に残ることが必要になるが、上記のように、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定すれば、筋状のムラの成分が反映されたプロファイルの部分で、筋状のムラの成分を的確に区分領域Ｓｋごとの補正データｄｓ(n)中に残すことが可能となる。従って、この点においても、分散σｖ２の値を筋状のムラに対応する分散の値と同程度の値に設定することが望ましいのである。
【０１３７】
なお、筋状のムラの成分がどの程度の大きさで出現するかは、個々の放射線画像撮影装置１で異なるため、分散σｖ２の値は、放射線画像撮影装置１ごとに、上記の目的を達成することが可能な適宜の値に設定される。
【０１３８】
以上のように、ウィナーフィルタを用いることで、前述した従来の場合のように、圧縮後の平均値ＤＡ(n)のプロファイルに対してハイパスフィルタ処理やローパスフィルタ処理を施さなくても、図１８に示すように、区分領域Ｓｋごとの補正データｄｓ(n)を的確に抽出することが可能となる。
【０１３９】
なお、上記の（５）式では、圧縮後の平均値ＤＡ(n)をそのまま用いる場合を示したが、例えば、圧縮後の平均値ＤＡ(n)にＳＩＮＣフィルタを施す等してローパスフィルタ処理を行い、そのように処理された圧縮後の平均値ＤＡ(n)に対して、上記の（５）式を適用するように構成することも可能である。
【０１４０】
その場合、上記（５）式中のＤＡ(n)にはローパスフィルタ処理が施された圧縮後の平均値ＤＡ(n)が代入され、ＤＡaveはローパスフィルタ処理が施された圧縮後の平均値ＤＡ(n)の平均値となる。
【０１４１】
次に、放射線画像処理装置７０の補正データ作成手段７５は、各区分領域Ｓｋごとに得られた区分領域Ｓｋごとの補正データｄｓ(n)のプロファイルに基づいて、補正データＤｓ（ｍ，ｎ）を作成するようになっている。
【０１４２】
具体的には、補正データ作成手段７５は、図１９に示すように、例えば、区分領域Ｓｋごとの補正データｄｓ(n)を、各区分領域Ｓｋにおける走査線方向Ａの例えば中央の位置にそれぞれ信号線方向Ｂに並べ、隣接する区分領域Ｓｋの対応する補正データｄｓ(n)同士の間の位置の補正データＤｓ（ｍ，ｎ）を例えば線形補間等により作成して、二次元状に配列された各補正データＤｓ（ｍ，ｎ）を作成する。
【０１４３】
このようにして、放射線画像処理装置７０の分割手段７１から補正データ作成手段７５における各処理で、元画像データＤori（ｍ，ｎ）から、筋状のムラ（図２９参照）の成分のみからなる補正データＤｓ（ｍ，ｎ）が抽出される。
【０１４４】
そして、放射線画像処理装置７０は、図１０に示したように、元画像データＤori（ｍ，ｎ）から、それぞれ対応する補正データＤｓ（ｍ，ｎ）を減算処理することによって、元画像データＤori（ｍ，ｎ）から筋状のムラを除去し、筋状のムラが除去された各画像データＤ＊を生成する。
【０１４５】
放射線画像処理装置７０は、筋状のムラが除去された各画像データＤ＊を生成すると、続いて、生成された各画像データＤ＊に対して正規化処理や対数変換処理、階調処理等の公知の画像処理を行って、生成された各画像データＤ＊を最終的な画像データ（すなわち診断用画像データ）に変換する（図９のステップＳ４）。
【０１４６】
なお、放射線画像撮影装置１を用いて放射線画像撮影を行う際、放射線画像撮影装置１の放射線入射面Ｒ（図２等参照）にグリッドを取り付ける等して撮影が行われる場合がある。このような場合には、グリッドにより放射線画像上にムラ（以下、上記の筋状のムラと区別するためグリッドモアレという。）が現れるが、このグリッドモアレは、上記の筋状のムラに比べて長周期で画像上に現れるため、上記のムラ除去処理（ステップＳ３）ではグリッドモアレを除去し切れず、或いは除去できない。
【０１４７】
そのため、このような場合には、放射線画像処理装置７０は、生成された各画像データＤ＊から生成した最終的な画像データに対して、公知のモアレ除去のための画像処理を施して、画像データからグリッドモアレを除去する処理を行う（ステップＳ５）。
【０１４８】
そして、放射線画像処理装置７０は、以上のようにして生成した最終的な画像データを記憶手段７０ｂに保存する（ステップＳ６）。また、必要に応じて、最終的な画像データに基づく放射線画像をイメージャでフィルム等の画像記録媒体に記録して出力したり、病院や医院のメインコンピュータに送信して格納させたりする。
【０１４９】
以上のように、本実施形態に係る放射線画像処理装置７０によれば、分割手段７１や平均値算出手段７２、エッジ圧縮手段７３、フィルタ処理手段７４、補正データ作成手段７５の各手段で、画像データＤに対して前処理が行われて得られた元画像データＤoriから筋状のムラの成分のみを的確に抽出し、筋状のムラの成分のみからなる補正データＤｓを作成する。
【０１５０】
前述したように、ＦＰＤ型の放射線画像撮影装置で放射線画像撮影を行う場合、ゲイン・オフセット補正処理や欠陥補正処理等の通常の前処理を行っただけの各元画像データＤoriに基づいて放射線画像を生成すると、放射線画像上に、ＦＰＤ型の放射線画像撮影装置特有の、走査線５の延在方向に延びる筋状のムラ（図２９参照）が生じる。
【０１５１】
しかし、平均値算出手段７２で、各区分領域Ｓｋ内の、同一の走査線５の延在方向に並ぶ各元画像データＤoriの平均値を算出することで、筋状のムラの成分が互いに打ち消しあうことなく、筋状のムラの成分を抽出することができる。また、フィルタ処理手段７４で、例えば上記（５）式に示したようなウィナーフィルタを適用することで、筋状のムラの成分が強く現れているデータを的確に抽出することが可能となる。
【０１５２】
そのため、分割手段７１から補正データ作成手段７５の各手段で前述した各処理を施すことにより、元画像データＤoriから走査線５の延在方向に延びる筋状のムラの成分のみを的確に抽出して、筋状のムラの成分のみからなる補正データＤｓを作成することが可能となる。
【０１５３】
そして、各元画像データＤoriからそれぞれ対応する補正データＤｓを減算処理することによって、各元画像データＤoriから、それに重畳されている筋状のムラを的確に除去することが可能となり、筋状のムラが除去された各画像データＤ＊を的確に生成することが可能となる。
【０１５４】
そのため、筋状のムラが除去された各画像データＤ＊に基づいて放射線画像を生成すれば、放射線画像上から筋状のムラが的確に除去され、筋状のムラのない放射線画像を生成することが可能となる。
【０１５５】
なお、前述したように、放射線画像処理装置７０の平均値算出手段７２で、各区分領域Ｓｋごとに、各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａ(n)を算出することにより、筋状のムラの成分以外のノイズ成分が打ち消しあって低減され、平均値Ｄａ(n)を、筋状のムラの成分をノイズ成分の主成分とするデータとすることができる。
【０１５６】
そのため、これを拡張すれば、図１２に示したように、分割手段７１で、各元画像データＤori（ｍ，ｎ）が二次元状に配列されて形成された領域Ｑを、少なくとも走査線方向Ａにおいて所定個数の元画像データＤoriごとに区分して、領域Ｑを各区分領域Ｓｋに分割せずに、領域Ｑの全域に対して各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａ(n)を算出するように構成すればよいとも考えられる。
【０１５７】
しかし、領域Ｑの全域に対して平均値Ｄａ(n)を算出するように構成すると、図１４にａ、ｂで示したようなエッジ部分が平均値Ｄａ(n)のプロファイル中に多数現れるようになり、エッジ圧縮手段７３で、上記のようにエッジ部分を適切に圧縮することが困難になる。
【０１５８】
また、上記の所定個数を小さい個数とすると（すなわち簡単に言えば、各区分領域Ｓｋの走査線方向Ａの幅を狭くすると）、上記のように平均値算出手段７２で各元画像データＤori（ｍ，ｎ）の走査線方向Ａの平均値Ｄａ(n)を算出する際に、筋状のムラの成分以外のノイズ成分が打ち消しあって低減される程度が低くなり、平均値Ｄａ(n)中に、筋状のムラの成分以外のノイズ成分が多く含まれる状態になる。
【０１５９】
そのため、そのような平均値Ｄａ(n)に対してエッジ圧縮手段７３以降の各手段で各処理を行っても、筋状のムラの成分を的確に抽出することができず、筋状のムラ（図２９参照）の成分のみからなる補正データＤｓを的確に作成することができない。そのため、分割手段７１で、領域Ｑを、走査線方向Ａにおいて所定個数の元画像データＤoriごとに区分する際の所定個数は、適宜の個数に設定される。
【０１６０】
その際、上記の所定個数を、領域Ｑにおける走査線方向Ａの元画像データＤoriの総数を自然数で除算した個数に設定すれば、平均値算出手段７２で、各区分領域Ｓｋごとに上記の平均値Ｄａ(n)を算出する際の平均化の対象となる元画像データＤoriの数が同じになり、演算処理の処理構成が単純になる。また、算出された各平均値Ｄａ(n)を均等に扱うことが可能となる。
【０１６１】
一方、上記のように、領域Ｑを分割する際の上記の所定個数を大きく設定すると、平均値Ｄａ(n)のプロファイル中にエッジ部分が多数現れるようになり、エッジ圧縮手段７３で、上記のようにエッジ部分を適切に圧縮することが困難になるが、エッジ部分を的確に圧縮して除去できないと、補正データＤｓ中にエッジ部分が残る。そして、元画像データＤoriから補正データＤｓを減算処理すると、元画像データＤori中のエッジ部分から補正データＤｓ中のエッジ部分が減算されてしまうため、筋状のムラが除去された各画像データＤ＊は、エッジ部分がぼやけたデータになってしまう。
【０１６２】
そして、このような各画像データＤ＊に基づいて放射線画像を生成すると、放射線画像は、例えば人体内の臓器や骨等が撮影された部分とそれ以外の部分との境界がぼやけた画像となり、画質が劣化する。このように、エッジ圧縮手段７３で、エッジ部分を適切に検出して的確に圧縮することが必要となる。
【０１６３】
そこで、上記の分割手段７１では、領域Ｑを細かく分割し、すなわち上記の所定個数を小さく設定して、エッジ部分が含まれない区分領域Ｓｋの数を多くし、また、区分領域Ｓｋにエッジ部分が含まれてもエッジ圧縮手段７３でエッジ部分を適切に検出できるようにして、エッジ圧縮手段７３やフィルタ処理手段７４での各処理を行う。
【０１６４】
そして、図２０に示すように、フィルタ処理手段７４と補正データ作成手段７５との間に再平均化手段７６を設け、上記のように細かく分割された各区分領域Ｓｋの各元画像データＤoriの走査線方向Ａの平均値Ｄａ(n)のプロファイルに対して、エッジ圧縮処理やフィルタ処理を行った後、再平均化手段７６で、隣接する所定個数の各区分領域Ｓｋについて、区分領域Ｓｋごとの補正データｄｓ(n)の走査線方向Ａの平均値ｄｓ(n)aveを算出するように構成することが可能である。
【０１６５】
補正データ作成手段７５は、この再平均化手段７６が算出した各平均値ｄｓ(n)aveに基づいて、例えば、図１９に示したように各平均値ｄｓ(n)aveを、上記のように隣接する所定個数の各区分領域Ｓｋの走査線方向Ａの中央の位置にそれぞれ信号線方向Ｂに並べ、線形補間する等して、領域Ｑにおける他の位置の各データを復元して、二次元状に配列された補正データＤｓを作成する。
【０１６６】
このように構成すれば、領域Ｑを細かく分割して生成された各区分領域Ｓｋのうち、エッジ部分が含まれない区分領域Ｓｋの数が多くなり、また、エッジ部分が含まれる区分領域Ｓｋでは、区分領域Ｓｋ中に含まれるエッジ部分の数が少なくなり、かつ明確に現れるようになるため、エッジ圧縮手段７３でエッジ部分を的確に検出して圧縮することが可能となる。
【０１６７】
また、エッジ部分を含まない区分領域Ｓｋや、エッジ部分が的確に圧縮された区分領域Ｓｋの各元画像データＤoriの走査線方向Ａの平均値ＤＡ(n)のプロファイルに対してフィルタ処理を行った後、再平均化手段７６で、隣接する所定個数の各区分領域Ｓｋについて、区分領域Ｓｋごとの補正データｄｓ(n)の走査線方向Ａの平均値ｄｓ(n)aveを算出することで、上記のように、筋状のムラ以外のノイズ成分が互いに打ち消しあって低減される。そのため、算出された平均値ｄｓ(n)aveが、筋状のムラの成分をノイズ成分の主成分とするデータになる。
【０１６８】
ところで、上記のように、エッジ圧縮手段７３で、各区分領域Ｓｋの各元画像データＤoriの走査線方向Ａの平均値Ｄａ(n)のプロファイルにおけるエッジ部分を圧縮する場合、仮にエッジ部分に筋状のムラの成分があったとしても、平均値Ｄａ(n)を圧縮することでその部分の筋状のムラの成分が消されてしまう。
【０１６９】
すなわち、圧縮前の平均値Ｄａ(n)に筋状のムラの情報が残っていたとしても、圧縮後の平均値ＤＡ(n)中からはその情報が失われてしまう。そのため、筋状をムラの成分を抽出するという観点から言えば、エッジ部分の平均値ＤＡ(n)の値は必ずしも信頼がおけるデータとは言えない。
【０１７０】
そこで、例えば、エッジ圧縮手段７３が、各区分領域Ｓｋのプロファイル中からエッジ部分を検出した場合、その位置、すなわち、例えば図１４に示した区分領域Ｓｋではｎ＝ａやｎ＝ｂを、各区分領域Ｓｋごとに記憶手段７０ｂに保存する。また、フィルタ処理手段７４によりフィルタ処理が行われた後で、各区分領域Ｓｋごとに得られた区分領域Ｓｋごとの補正データｄｓ(n)のプロファイルに基づいて、全ての区分領域Ｓｋについて、各区分領域Ｓｋごとの補正データｄｓ(n)の走査線方向Ａの平均値ｄｓ(n)allのプロファイルを作成するように構成する。
【０１７１】
そして、各区分領域Ｓｋのエッジ部分の区分領域Ｓｋごとの補正データｄｓ(n)については、各区分領域Ｓｋごとに算出したデータｄｓ(n)を破棄して全区分領域Ｓｋについての平均値ｄｓ(n)allで置換することにより、エッジ圧縮手段７３におけるエッジ圧縮処理で失われた各区分領域Ｓｋのエッジ部分の筋状のムラの成分を復元することが可能となる。
【０１７２】
［第２の実施の形態］
前述したように、画像データＤや元画像データＤori上に重畳される筋状のムラは、前述したように、走査駆動手段１５の電源回路１５ａ（図７参照）やバイアス電源１４等で発生する電圧のノイズが走査線５やバイアス線９を介して各放射線検出素子７やＴＦＴ８に伝わることに起因して生じるものである。
【０１７３】
そして、画像データＤの読み出し処理の際に、走査線５の各ラインＬ１〜Ｌｘにオン電圧が印加されるタイミングが異なるため、同じ走査線５に接続されている各放射線検出素子７には同じノイズが重畳されるが、走査線５が異なるとノイズの大きさが異なる値になる。そのため、図２９に拡大して示したように、放射線画像ｐすなわち各画像データＤ上に、走査線方向Ａに延びる筋状のムラとなって現れるものである。
【０１７４】
そして、上記の第１の実施形態では、図５等に示したように、少なくとも各走査線５が検出部Ｐの全域（すなわち図５の場合は左側から右側までの全ての領域）にわたって配列されているため、ゲイン・オフセット補正処理等が施された各元画像データＤoriを、検出部Ｐ上で二次元状に配列された複数の放射線検出素子７にそれぞれ対応するように、二次元状に配列して、分割手段７１以降の各手段で各処理を行った。
【０１７５】
しかし、ＦＰＤ型の放射線画像撮影装置の中には、例えば図２１に示すように、検出部Ｐ上で、少なくとも各走査線５が、走査線方向Ａにおいてそれぞれ走査線５ａと走査線５ｂに分割されて構成されているような放射線画像撮影装置１＊もある。
【０１７６】
このような放射線画像撮影装置１＊では、分割された各走査線５ａ、５ｂがそれぞれ存在する検出部Ｐの各領域Ｐａ、Ｐｂで、走査線５ａ、５ｂ等を介して各放射線検出素子７等に伝わるノイズが各領域Ｐａ、Ｐｂで異なる場合があり、また、各領域Ｐａ、Ｐｂで必ずしも同じタイミングで読み出し処理等が行われない場合もある。
【０１７７】
このような場合、仮に、放射線画像撮影装置１＊で撮影された各画像データＤに前処理が施されて作成された各元画像データＤoriに基づいて放射線画像ｐ１を生成すると、図２２に拡大して示すように、検出部Ｐの各領域Ｐａ、Ｐｂに対応する各領域で、走査線方向Ａに延びる筋状のムラがそれぞれ異なる状態で現れる。
【０１７８】
なお、図２２や後述する図２３、図２４において、放射線画像中に信号線方向Ｂに延びる筋状のムラは、各信号線６ごとに接続されている各読み出し回路１７の読み出し特性に由来するものであり、本発明における走査線方向Ａに延びる筋状のムラとは異なるものである。この信号線方向Ｂに延びる筋状のムラについては、別途、除去する処理が適切に行われる。
【０１７９】
そして、上記の第１の実施形態と同様に、このような元画像データＤoriを、検出部Ｐの全域で二次元状に配列された複数の放射線検出素子７にそれぞれ対応するように、二次元状に配列して、分割手段７１以降の各手段で各処理を行うと、図２３に拡大して示すように、放射線画像ｐ２上の、特に各領域Ｐａ、Ｐｂの境界Ｂｏ付近で、各領域Ｐａ、Ｐｂで筋状のムラが異なる状態で現れることの影響が排除し切れず、筋状のムラが除去し切れずに残存するという現象が生じる場合がある。
【０１８０】
そこで、このように、検出部Ｐ上で少なくとも各走査線５が走査線方向Ａにおいて分割されて構成されている放射線画像撮影装置１＊により撮影された各画像データＤに対しては、分割された走査線５ａ、５ｂがそれぞれ存在する検出部Ｐの各領域Ｐａ、Ｐｂでそれぞれ得られた各画像データＤについてそれぞれ別々に第１の実施形態で説明した画像処理を行うように構成することが望ましい。
【０１８１】
このように構成すれば、検出部Ｐの各領域Ｐａ、Ｐｂでそれぞれ得られた各画像データＤに基づく各元画像データＤoriから、それぞれ別々に補正データＤｓが抽出されて作成され、それらの補正データＤｓがそれぞれ別々に各元画像データＤoriから減算処理される。
【０１８２】
そのため、検出部Ｐの領域Ｐａで得られた各画像データＤに基づく各元画像データＤoriからはそれから抽出されて作成された補正データＤｓが減算処理され、検出部Ｐの領域Ｐｂで得られた各画像データＤに基づく各元画像データＤoriからはそれから抽出されて作成された補正データＤｓが減算処理されるため、検出部Ｐの各領域Ｐａ、Ｐｂで、筋状のムラが除去された各画像データＤ＊をそれぞれ的確に生成することが可能となる。
【０１８３】
そのため、図２４に拡大して示すように、上記のように的確に生成された各画像データＤ＊に基づいて生成された放射線画像ｐ３では、筋状のムラが的確に除去され、各領域Ｐａ、Ｐｂの境界Ｂｏ付近においても、筋状のムラが残存していない状態とすることが可能となる。
【０１８４】
［第３の実施の形態］
第１および第２の実施形態では、放射線画像撮影装置１で読み出された各画像データＤについて、ゲイン・オフセット補正処理等の前処理を行って各元画像データＤoriを作成し、各元画像データＤoriに対して分割手段７１から補正データ作成処理７５までの各処理を行って、補正データＤｓを作成し、各元画像データＤoriから対応する補正データＤｓを減算処理して、筋状のムラが除去された各画像データＤ＊を生成する場合について説明した。
【０１８５】
一方、ゲイン・オフセット補正処理（図９のステップＳ１）で各画像データＤから減算されるオフセット補正値Ｏを取得するためのダーク読取処理は、前述したように、放射線画像撮影装置１に放射線を照射しないことを除いて放射線画像撮影と同様の処理が行われて読み出されるものであるため、オフセット補正値Ｏにも、画像データＤと同様の走査線方向Ａに延びる筋状のムラが重畳されている。
【０１８６】
また、ダーク読取処理は、放射線画像撮影と同時に行うことができず、放射線画像撮影の前や後に行われるため、ダーク読取処理で読み出される各オフセット補正値Ｏに重畳される筋状のムラと、放射線画像撮影における読み出し処理で読み出される各画像データＤに重畳される筋状のムラは、全く異なる状態となる。
【０１８７】
しかし、ゲイン・オフセット補正処理で、上記（１）式に従って各画像データＤから対応する各オフセット補正値Ｏを減算すると、減算して得られる値Ｄ−Ｏには、各画像データＤに重畳されている筋状のムラから各オフセット補正値Ｏに重畳されている筋状のムラが減算されて生成される筋状のムラが重畳される状態になる。
【０１８８】
従って、この筋状のムラが重畳された値Ｄ−Ｏにゲイン補正値Ｇを乗算して各元画像データＤoriを作成し、これらの各元画像データＤoriに対して画像処理を行うように構成することも可能である。これは、第１の実施形態や第２の実施形態で説明した場合に相当する。
【０１８９】
しかし、筋状のムラが重畳された各画像データＤから、別の筋状のムラが重畳された各オフセット補正値Ｏを減算すると、減算して得られる値Ｄ−Ｏに重畳されている筋状のムラは、その大きさが、２の平方根倍すなわち約１．４倍に大きくなる。そして、画像データＤに重畳される筋状のムラがもともと比較的大きいような場合には、減算して得られる値Ｄ−Ｏに重畳される筋状のムラの大きさが異常に大きくなり、元画像データＤoriから筋状のムラを的確に除去できなくなる場合がある。
【０１９０】
そこで、ゲイン・オフセット補正処理（図９のステップＳ１）で、上記（１）式に従って各画像データＤからオフセット補正値Ｏを減算する前に、このオフセット補正値Ｏについて、第１および第２の実施形態で示した各元画像データＤoriに対する各処理と同様に処理を行い、筋状のムラが除去された各オフセット補正値Ｏを生成しておく。
【０１９１】
そして、ゲイン・オフセット補正処理（図９のステップＳ１）で、この筋状のムラが除去された各オフセット補正値Ｏを、筋状のムラが重畳されている画像データＤから減算するように構成することが可能である。
【０１９２】
このように構成すれば、元画像データＤoriに重畳されている筋状のムラが異常に大きくならず、元の画像データＤに重畳されていた筋状のムラの状態のままで各処理を行うことが可能となり、各元画像データＤoriから筋状のムラを的確に除去することが可能となる。
【０１９３】
また、オフセット補正値Ｏは、前述したように、放射線画像撮影装置１に放射線を照射しない状態で行われるダーク読取処理で取得される。そのため、各元画像データＤoriの場合と異なり、オフセット補正値Ｏには、人体内の臓器や骨等の情報は含まれておらず、上記のようなエッジ部分が存在しない。
【０１９４】
そのため、オフセット補正値Ｏに対する処理では、各オフセット補正値Ｏを、検出部Ｐ（図５や図７参照）上で二次元状に配列された複数の放射線検出素子７にそれぞれ対応するように、二次元状に配列した場合の領域Ｑを、分割手段７１で分割せず、領域Ｑの全体に対して平均値算出手段７２から補正データ作成手段７５までの各手段での各処理を行うように構成することも可能である。
【０１９５】
一方、オフセット補正値Ｏから筋状のムラをより厳密に除去するために、以下のように構成することも可能である。
【０１９６】
すなわち、事前に、放射線画像撮影装置１でダーク読取処理を複数回行って、各放射線検出素子７ごとに複数のオフセット補正値Ｏを取得しておき、それらの平均値として各放射線検出素子７ごとのオフセット補正値Ｏを算出して、放射線画像処理装置７０の記憶手段７０ｂに保存させておく。
【０１９７】
以下、この各放射線検出素子７ごとに予め取得されたオフセット補正値Ｏを、各回の放射線画像撮影ごとに得られるオフセット補正値Ｏと区別するため、基準となるオフセット補正値Ｏｃという。
【０１９８】
また、上記のように、基準となるオフセット補正値Ｏｃを取得するためのダーク読取処理は複数回行われ、各回のダーク読取処理ごとにオフセット補正値Ｏに重畳される筋状のムラの状態が異なるため、すなわち、図２９に示したように、筋状のムラが強く出る走査線５のラインＬと弱く出る走査線５のラインＬが各回のダーク読取処理ごとに異なるため、各回のダーク読取処理で取得された各オフセット補正値Ｏを平均化することで、各オフセット補正値Ｏに重畳されている筋状のムラが互いに打ち消しあう。そのため、基準となるオフセット補正値Ｏｃは、筋状のムラが重畳されていない状態となる。
【０１９９】
放射線画像処理装置７０では、このような基準となるオフセット補正値Ｏｃを有する状態で、図２５に示すフローチャートに従って、各回の放射線画像撮影の前や後に行われるダーク読取処理で取得された実際のオフセット補正値Ｏから筋状のムラを除去して、筋状のムラが除去されたオフセット補正値Ｏを作成する。
【０２００】
放射線画像処理装置７０は、まず、各回のダーク読取処理で取得された実際のオフセット補正値Ｏから基準となるオフセット補正値Ｏｃを各放射線検出素子７ごとに減算して、各放射線検出素子７ごとの各差分ΔＯを算出する（ステップＳ１１）。上記のように、基準となるオフセット補正値Ｏｃには筋状のムラが重畳されていないため、各差分ΔＯには、実際のオフセット補正値Ｏに重畳されていた筋状のムラが重畳された状態となる。
【０２０１】
そして、算出した各差分ΔＯに対して、前述した分割手段７１から補正データ作成手段７５における各処理を実施する（ステップＳ１２）。なお、この場合も、各差分ΔＯには、人体内の臓器や骨等の情報は含まれておらず、上記のようなエッジ部分が存在しないため、必ずしも分割手段７１における分割処理を行う必要はない。
【０２０２】
そして、補正データ作成手段７５で、各差分ΔＯに関する補正データΔＯｓを作成すると、元の差分ΔＯからそれぞれ対応する補正データΔＯｓを減算処理して、筋状のムラが除去された各差分ΔＯ＊を生成する。
【０２０３】
次に、放射線画像処理装置７０は、生成した筋状のムラが除去された各差分ΔＯ＊を、基準となる各オフセット補正値Ｏｃに加算する（ステップＳ１３）。前述したように、基準となるオフセット補正値Ｏｃには筋状のムラが重畳されておらず、各差分ΔＯ＊からも筋状のムラが除去されているため、これらの加算値も筋状のムラが除去されたものとなる。
【０２０４】
そのため、放射線画像処理装置７０は、筋状のムラが除去された各差分ΔＯ＊と基準となる各オフセット補正値Ｏｃの加算値を、筋状のムラが除去されたオフセット補正値Ｏ＊として記憶手段７０ｂに保存する。放射線画像処理装置７０は、このようにして、実際のオフセット補正値Ｏから、筋状のムラが除去されたオフセット補正値Ｏ＊を作成する（ステップＳ１４）。
【０２０５】
なお、このようにして作成された筋状のムラが除去されたオフセット補正値Ｏ＊を、上記（１）式のＯに代入して元画像データＤoriが作成され、この元画像データＤoriに対して、再度、分割手段７１から補正データ作成手段７５における各処理が施されて、作成された補正データＤｓを元画像データＤoriから減算処理して、筋状のムラが除去された各画像データＤ＊が生成されることは、改めて説明するまでもない。
【０２０６】
このように構成すれば、オフセット補正値Ｏから筋状のムラがより的確に除去されるため、作成された筋状のムラが除去されたオフセット補正値Ｏ＊を画像データＤから減算して得られる元画像データＤoriに対して適切に筋状のムラを除去する各処理を施すことが可能となり、筋状のムラが除去された各画像データＤ＊を的確に生成することが可能となる。
【産業上の利用可能性】
【０２０７】
放射線画像撮影を行う分野（特に医療分野）において利用可能性がある。
【符号の説明】
【０２０８】
１  放射線画像撮影装置
５、Ｌ１〜Ｌｘ  走査線
６  信号線
７  放射線検出素子
７０  放射線画像処理装置
７１  分割手段
７２  平均値算出手段
７３  エッジ圧縮手段
７４  フィルタ処理手段
７５  補正データ作成手段
７６  再平均化手段
Ａ  走査線方向（走査線の延在方向）
Ｂ  信号線方向（信号線の延在方向）
Ｄ  画像データ
Ｄ＊  ノイズが除去された画像データ
ＤＡave  差分が圧縮されたプロファイル上の各平均値の各区分領域ごとの平均値
Ｄａ(n)  平均値
Ｄａ(n)−Ｄａ(n-1)  差分
Ｄori  元画像データ（画像データ）
Ｄｓ、Ｄｓ（ｍ，ｎ）  補正データ
ｄｓ(n)  区分領域ごとの補正データ（適応フィルタが施されたプロファイル上の各データ）
ｄｓ(n)ave  区分領域ごとの補正データの平均値
Ｏ  オフセット補正値
Ｏ＊  ノイズが除去されたオフセット補正値
Ｏｃ  基準となるオフセット補正値
Ｐ  検出部
Ｐａ、Ｐｂ  検出部の各領域
Ｑ  領域
ｒ  区画
Ｓｋ  各区分領域
ΔＤａth  閾値
ΔＯ  差分
ΔＯ＊  ノイズが除去された差分
ΔＯｓ  補正データ
σ２(n)  差分が圧縮されたプロファイル上の平均値の分散
σｖ２  予め設定された分散の値
【図１】
2011145377.tif 000005
【図２】
2011145377.tif 000006
【図３】
2011145377.tif 000007
【図４】
2011145377.tif 000008
【図５】
2011145377.tif 000009
【図６】
2011145377.tif 000010
【図７】
2011145377.tif 000011
【図８】
2011145377.tif 000012
【図９】
2011145377.tif 000013
【図１０】
2011145377.tif 000014
【図１１】
2011145377.tif 000015
【図１２】
2011145377.tif 000016
【図１３】
2011145377.tif 000017
【図１４】
2011145377.tif 000018
【図１５】
2011145377.tif 000019
【図１６】
2011145377.tif 000020
【図１７】
2011145377.tif 000021
【図１８】
2011145377.tif 000022
【図１９】
2011145377.tif 000023
【図２０】
2011145377.tif 000024
【図２１】
2011145377.tif 000025
【図２２】
2011145377.tif 000026
【図２３】
2011145377.tif 000027
【図２４】
2011145377.tif 000028
【図２５】
2011145377.tif 000029
【図２６】
2011145377.tif 000030
【図２７】
2011145377.tif 000031
【図２８】
2011145377.tif 000032
【図２９】
2011145377.tif 000033
【国際調査報告】
2011145377.tif 000034
2011145377.tif 000035
2011145377.tif 000036
2011145377.tif 000037