(57)【要約】【課題】線状ノイズ成分を除去するときに線状ノイズ成分の抽出時間の短縮を図る。【解決手段】放射線画像Ｐ０を垂直方向または水平方向または両方向に縮小した縮小画像Ｐ１が生成され、生成した縮小画像Ｐ１（Ｐ２）の中から矢印Ｙ方向の線状ノイズ成分ＬＮを示す線状ノイズ画像ＮＰ１が抽出される。その後、抽出された線状ノイズ画像ＮＰ１が放射線画像Ｐ０の画像サイズに拡大され、拡大された線状ノイズ画像ＮＰ０を用いて放射線画像Ｐ０から線状ノイズ成分が除去される。【選択図】図１
2010009348.tif 000002
【特許請求の範囲】
【請求項１】
画像を主方向または副方向または両方向ともに縮小した縮小画像を生成する画像縮小手段と、
該画像縮小手段において生成された前記縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成するノイズ画像抽出手段と、
該ノイズ画像抽出手段により生成された前記線状ノイズ画像を前記画像の画像サイズに拡大する画像拡大手段と、
該画像拡大手段により拡大された前記線状ノイズ画像を用いて前記画像から前記線状ノイズ成分を除去する画像補正手段と
を備えたことを特徴とする画像処理装置。
【請求項２】
前記ノイズ画像抽出手段が、
前記縮小画像から前記線状ノイズと直交する方向の高周波成分を抽出する周波数処理手段と、
該周波数処理手段において抽出された前記高周波成分から前記線状ノイズに沿う方向の低周波成分を抽出することにより前記線状ノイズ成分を示す前記線状ノイズ画像を生成するノイズ画像生成手段と
を有するものであることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記ノイズ画像抽出手段が、
前記縮小画像から前記線状ノイズと直交する方向の高周波成分を抽出する周波数処理手段と、
該周波数処理手段において抽出された前記高周波成分について前記線状ノイズに沿う方向の画素ラインを複数の分割画素ラインに分割するライン分割手段と、
該ライン分割手段により分割された前記各分割画素ラインのメディアン値を算出するメディアン算出手段と、
該メディアン算出手段により算出された前記各分割画素ライン毎の前記メディアン値を用いて前記線状ノイズ成分を示す線状ノイズ画像を生成するノイズ画像生成手段と
を有するものであることを特徴とする請求項１記載の画像処理装置。
【請求項４】
前記ノイズ画像抽出手段が、
前記縮小画像から前記線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、
該周波数処理手段において抽出された前記高周波成分から前記線状ノイズに沿う方向の低周波成分を抽出することにより前記線状ノイズ成分を示す第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、
前記周波数処理手段において抽出された前記高周波成分について前記画素ラインを複数の分割画素ラインに分割するライン分割手段と、
該ライン分割手段により分割された前記各分割画素ラインのメディアン値を算出するメディアン算出手段と、
該メディアン算出手段により算出された前記各分割画素ライン毎の前記メディアン値を用いて前記線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段と
を備えたものであることを特徴とする請求項１記載の画像処理装置。
【請求項５】
前記ノイズ画像抽出手段が、
前記縮小画像から前記線状ノイズに直交する方向の高周波成分と低周波成分とを抽出する周波数処理手段と、
該周波数処理手段において抽出された前記高周波成分から前記線状ノイズに沿う方向の低周波成分を抽出することにより前記線状ノイズ成分を示す第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、
該周波数処理手段において抽出された前記低周波成分について前記線状ノイズに沿う方向の画素ラインを複数の分割画素ラインに分割するライン分割手段と、
該ライン分割手段により分割された前記各分割画素ラインのメディアン値を算出するメディアン算出手段と、
該メディアン算出手段により算出された前記各分割画素ライン毎の前記メディアン値を用いて前記線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段と
を備えたものであり、
前記画像補正手段が、前記画像拡大手段により拡大された前記第１線状ノイズ画像と前記第２線状ノイズ画像とを用いて前記画像から前記線状ノイズ成分を除去するものであることを特徴とする請求項１記載の画像処理装置。
【請求項６】
前記ノイズ画像生成手段が、前記各分割画素ライン毎に算出された複数の前記メディアン値に基づいて、前記画素ラインの画素位置と画素値との関係を示す関数を導出し、導出した該関数を用いて各画素の画素値を算出し前記線状ノイズ画像を生成するものであることを特徴とする請求項３から５のいずれか１項記載の画像処理装置。
【請求項７】
前記画像補正手段が、前記第１線状ノイズ画像と前記第２線状ノイズ画像との各画素の絶対値を比較し、絶対値の小さい値を用いて前記画像から前記線状ノイズ成分を除去するものであることを特徴とする請求項４から６のいずれか１項記載の画像処理装置。
【請求項８】
前記画像補正手段が、前記画像の性質に合わせて前記第１線状ノイズ画像もしくは前記第２線状ノイズ画像のいずれかを選択して前記画像から前記線状ノイズ成分を除去するものであることを特徴とする請求項４から７のいずれか１項記載の画像処理装置。
【請求項９】
前記画像補正手段が、前記第１線状ノイズ画像を用いて前記線状ノイズ成分を除去した前記画像と、前記第２線状ノイズ画像を用いて前記線状ノイズ成分を除去した前記画像とを生成するものであることを特徴とする請求項４から８のいずれか１項記載の画像処理装置。
【請求項１０】
前記縮小画像に前記線状ノイズに直交する方向における隣接した画素間の画素値の急激な変化を緩やかにする前処理を施す前処理手段をさらに備え、
前記ノイズ画像抽出手段が前記前処理手段により前処理された前記縮小画像から前記線状ノイズ画像を抽出するものであることを特徴とする請求項１から９のいずれか１項記載の画像処理装置。
【請求項１１】
前記ライン分割手段が、前記画像の素抜け部分を除いて前記分割画素ラインに分割し、線状ノイズ成分を生成するものであることを特徴とする請求項１から１０のいずれか１項記載の画像処理装置。
【請求項１２】
画像を主方向または副方向または両方向ともに縮小した縮小画像を生成し、
生成した前記縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成し、
生成した前記線状ノイズ画像を前記画像の画像サイズに拡大し、
拡大した前記線状ノイズ画像を用いて前記画像から前記線状ノイズ成分を除去する
ことを特徴とする画像処理方法。
【請求項１３】
コンピュータに、
画像を主方向または副方向または両方向ともに縮小した縮小画像を生成し、
生成した前記縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成し、
生成した前記線状ノイズ画像を前記画像の画像サイズに拡大し、
拡大した前記線状ノイズ画像を用いて前記画像から前記線状ノイズ成分を除去する
ことを実行させるための画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像から線状ノイズを除去するための画像処理装置および方法ならびにプログラムに関するものである。
【背景技術】
【０００２】
従来、被写体に放射線を曝射し被写体を透過した放射線量を放射線画像として検出し、放射線画像を用いて被写体の診断が行われている。この放射線画像を取得する装置として撮影する被写体の部位毎に異なる装置が用いられている。たとえば乳房を撮影する際にはマンモグラフィ装置、胸部レントゲン等の２次元の放射線画像を撮影する際には一般撮影装置がそれぞれ用いられる。
【０００３】
上述した各種装置により取得した放射線画像には線状ノイズ（スジムラ）が発生することがある。たとえばラインセンサにより放射線検出器に蓄積された画像情報を読取る場合、ラインセンサを構成する各受光部毎の特性の違いによる副方向（放射線画像の垂直方向）の線状ノイズが発生する場合がある。あるいはラインセンサへ加えられた衝撃により、主方向（放射線画像の水平方向）の線状ノイズが現れることがある。
【０００４】
上述した線状ノイズを放射線画像から除去する方法として放射線画像に対しフィルタ処理する方法が知られている（たとえば特許文献１参照）。特許文献１において、放射線画像データから線状ノイズ成分を抽出した後、抽出した線状ノイズ成分を放射線画像から差し引く方法が提案されている。特に、線状ノイズ成分を抽出するとき、画素値の急激な変化を緩やかにする前処理を行った後、フィルタ処理により、線状ノイズ成分を抽出している。
【特許文献１】特開２００５−８４９０２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、上述した特許文献１に示す線状ノイズ成分をフィルタ処理により抽出するときに必要な処理時間は放射線画像の画像サイズとフィルタのサイズとに依存する。線状ノイズ成分を抽出するには多くの計算が必要であり時間がかかる為、線状ノイズ成分の抽出処理を短くすることが望まれている。
【０００６】
そこで、本発明は、線状ノイズ成分を除去するときに線状ノイズ成分の抽出時間の短縮を図ることができる画像処理装置および方法ならびにプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００７】
本発明の画像処理装置は、画像を主方向または副方向または両方向に縮小した縮小画像を生成する画像縮小手段と、画像縮小手段において生成された縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成するノイズ画像抽出手段と、ノイズ画像抽出手段により生成された線状ノイズ画像を元画像の画像サイズに拡大する画像拡大手段と、画像拡大手段により拡大された線状ノイズ画像を用いて元画像から線状ノイズ成分を除去する画像補正手段とを備えたことを特徴とするものである。
【０００８】
本発明の画像処理方法は、画像を主方向または副方向または両方向に縮小した縮小画像を生成し、生成した縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成し、生成した線状ノイズ画像を元画像の画像サイズに拡大し、拡大した線状ノイズ画像を用いて元画像から線状ノイズ成分を除去することを特徴とするものである。
【０００９】
本発明の画像処理プログラムは、コンピュータに、画像を主方向または副方向または両方向に縮小した縮小画像を生成し、生成した縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成し、生成した線状ノイズ画像を元画像の画像サイズに拡大し、拡大した線状ノイズ画像を用いて元画像から線状ノイズ成分を除去することを実行させることを特徴とするものである。
【００１０】
ここで、画像は、なんでもよく、たとえばマンモグラフィ装置により取得された放射線画像や胸部撮影により取得された放射線画像であってもよいし、医用画像に限らず非破壊検査に用いられる放射線画像であってもよい。さらにはオフセット補正やシェーディング補正等に用いられる補正用の放射線画像であってもよい。
【００１１】
また、画像縮小手段は、主方向（たとえば画像の水平方向）または副方向（たとえば画像の垂直方向）のいずれか一方もしくは両方向に縮小するものであればよく、線状ノイズに沿う方向と、さらに線状ノイズが低周波成分のみである場合には、線状ノイズに直交する方向にさらに縮小することが可能である。
【００１２】
なお、ノイズ画像抽出手段は、画像内の線状ノイズを抽出するものであればどのような構成であってもよい。たとえば、ノイズ画像抽出手段は、縮小画像から線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、周波数処理手段により高周波成分から線状ノイズに沿う方向の低周波成分を抽出することにより線状ノイズ画像を生成するノイズ画像生成手段とを有するものであってもよい。
【００１３】
あるいは、ノイズ画像抽出手段は、縮小画像から線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、周波数処理手段により抽出された高周波成分について線状ノイズに沿う方向の画素ラインを複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す線状ノイズ画像を生成するノイズ画像生成手段とを有するものであってもよい。
【００１４】
ここで、ライン分割手段は、画素ラインを複数の分割画素ラインに分割するものであれば分割画素ラインの画素数は問わないが、たとえば画素ラインを生体構造に存在し得る直線の２倍以上の長さからなる分割画素ラインに分割するものであるのが望ましい。なお、生体構造に存在し得る直線の２倍以上の長さとは、生体構造に含まれる直線成分の２倍の長さであって、経験的もしくは統計的に求められるものであり予め所定の長さに設定されるものである。
【００１５】
さらには、ノイズ画像抽出手段は、縮小画像から線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、周波数処理手段において抽出された高周波成分から線状ノイズに沿う方向の低周波成分を抽出することにより第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、周波数処理手段により抽出された高周波成分について線状ノイズに沿う方向の画素ラインを生体構造に存在し得る直線の２倍以上の長さからなる複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段とを備えたものであってもよい。このとき、画像補正手段は、第１線状ノイズ画像と第２線状ノイズ画像を用いて作成した線状ノイズ画像を画像拡大手段により拡大し元画像から線状ノイズ成分を除去するものであってもよい。
【００１６】
また、ノイズ画像抽出手段は、縮小画像から線状ノイズに直交する方向の高周波成分と低周波成分とを抽出する周波数処理手段と、周波数処理手段において抽出された高周波成分から線状ノイズに沿う方向の低周波成分を抽出することにより第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、周波数処理手段により抽出された低周波成分について線状ノイズに沿う方向の画素ラインを生体構造に存在し得る直線の２倍以上の長さからなる複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段とを備えたものであってもよい。
【００１７】
このとき、ノイズ画像生成手段は、各分割画素ライン毎に算出された複数のメディアン値に基づいて、画素ラインの画素の位置の画素値を示す関数を求め、関数により導出された画素値を画素値とする線状ノイズ画像を生成するものであってもよい。
【００１８】
また、画像補正手段は、第１線状ノイズ画像と第２線状ノイズ画像とを用いて線状ノイズ画像を生成するものであればよく、第１線状ノイズ画像と第２線状ノイズ画像との各画素の絶対値を比較し、絶対値の小さい値を用いることにより線状ノイズ画像を作成するものであってもよいし、画像の性質に合わせて第１線状ノイズ画像もしくは第２線状ノイズ画像のいずれかを選択することにより線状ノイズ画像を作成するものであってもよい。ここで、放射線画像の性質とは、たとえば人体構造が撮影された放射線画像と、ベタ画像や幾何学パターンが撮影された補正用画像やテストパターン画像のように放射線画像に含まれる成分の性質を意味する。
【００１９】
さらに、画像補正手段は、第１線状ノイズ画像を用いて線状ノイズ成分を除去した放射線画像と、第２線状ノイズ画像を用いて線状ノイズ成分を除去した放射線画像とを生成するものであってもよい。
【００２０】
また、放射線画像処理装置は、縮小画像に対し線状ノイズに直交する方向における隣接した画素間の画素値の急激な変化を緩やかにする前処理を施す前処理手段をさらに備えていてもよい。このとき、ノイズ画像抽出手段は前処理を施した縮小画像から線状ノイズ画像の生成を行う。
【発明の効果】
【００２１】
本発明の画像処理装置および方法ならびにプログラムによれば、画像を主方向または副方向または両方向に縮小した縮小画像を生成し、生成した縮小画像から線状ノイズ成分を抽出し線状ノイズ画像を生成し、生成した線状ノイズ画像を元画像の画像サイズに拡大し、拡大した線状ノイズ画像を用いて元画像から線状ノイズ成分を除去することにより、線状ノイズは線状ノイズが延びる方向への濃度変化が少なく、縮小画像から抽出した線状ノイズと元の画像から抽出した線状ノイズとはほぼ同等の抽出結果になることを利用し、縮小画像に対しフィルタリング処理を行うことにより線状ノイズを抽出するため、線状ノイズの抽出精度を維持しながら処理時間の短縮を図ることができる。
【００２２】
また、ノイズ画像抽出手段が、縮小画像から線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、周波数処理手段により抽出された高周波成分について線状ノイズに沿う方向の画素ラインを複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す線状ノイズ画像を生成するノイズ画像生成手段とを有するものであれば、縮小画像を分割した分割画素ライン毎にメディアン値を求めればよいため、メディアン値の算出時間を短縮することができるとともに、メディアン値として人体構造を示す値が算出されることを最小限に抑えることができるため、人体のパターンの中にスジ方向に伸びている部位がある場合、その部位も線状ノイズとして除去されるのを防止し線状ノイズ成分のみ除去し、画像補正による画質の劣化が生じるのを防止することができる。
【００２３】
さらに、ノイズ画像抽出手段が、縮小画像から線状ノイズに直交する方向の高周波成分を抽出する周波数処理手段と、周波数処理手段において抽出された高周波成分から線状ノイズに沿う方向の低周波成分を抽出することにより線状ノイズ成分を示す第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、周波数処理手段により抽出された高周波成分について線状ノイズに沿う方向の画素ラインを生体構造に存在し得る直線の２倍以上の長さからなる複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段とを備えたものであるとき、２種類の第１線状ノイズ画像と第２線状ノイズ画像とを用いて線状ノイズ画像が生成されるため、補正後の画像におけるアーチファクトの発生を低減することができる。
【００２４】
また、ノイズ画像抽出手段が、縮小画像から線状ノイズに直交する方向の高周波成分と低周波成分とを抽出する周波数処理手段と、周波数処理手段において抽出された高周波成分から線状ノイズに沿う方向の低周波成分を抽出することにより第１線状ノイズ画像を生成する第１ノイズ画像生成手段と、周波数処理手段により抽出された低周波成分について線状ノイズに沿う方向の画素ラインを生体構造に存在し得る直線の２倍以上の長さからなる複数の分割画素ラインに分割するライン分割手段と、ライン分割手段により分割された各分割画素ラインのメディアン値を算出するメディアン算出手段と、メディアン算出手段により算出された各分割画素ライン毎のメディアン値を用いて線状ノイズ成分を示す第２線状ノイズ画像を生成する第２ノイズ画像生成手段とを備えたものであり、画像補正手段が、第１線状ノイズ画像と第２線状ノイズ画像とを用いて作成した線状ノイズ画像を画像拡大手段により拡大し元画像から線状ノイズ成分を除去するものであれば、画像内の人体構造を示す成分は低周波成分が多いことを利用し、補正によるアーチファクトの発生を低減することができる。
【００２５】
また、画像補正手段が、第１線状ノイズ画像と第２線状ノイズ画像との各画素の絶対値を比較し、画素の絶対値の小さい値を用いることにより放射線画像から線状ノイズ成分を除去するものであれば、線状ノイズである確率が高い成分を用いて線状ノイズ画像を作成することができるため、補正によるアーチファクトを低減させることができる。
【００２６】
さらに、画像補正手段が、画像の種類に合わせて第１線状ノイズ画像もしくは第２線状ノイズ画像のいずれかを選択することにより画像から線状ノイズ成分を除去するものであるとき、たとえば画像内に人体構造を示す短い線状成分が少ない場合には第１線状ノイズ画像を補正に用いる線状ノイズ画像として選択して確実に線状ノイズ画像を除去し、画像内に人体構造を示す短い線状成分が含まれている場合には第２線状ノイズ画像を補正に用いる線状ノイズ画像として選択してアーチファクトを低減することができる。
【００２７】
さらに、画像補正手段が、第１線状ノイズ画像を用いて線状ノイズ成分を除去した画像と、第２線状ノイズ画像を用いて線状ノイズ成分を除去した画像とを生成するものであるとき、処理内容の異なる複数の補正済みの画像を提供することができる。
【００２８】
また、縮小画像に対し線状ノイズに直交する方向における画素値の急激な変化を緩やかにする前処理を施す前処理手段をさらに備え、ノイズ画像抽出手段が前処理手段により前処理された縮小画像から線状ノイズ画像を抽出するものであるとき、周波数処理手段において周波数処理をする際に急激な画素値の変化部分にリンギング成分が発生するのを防止してアーチファクトを低減することができる。
【発明を実施するための最良の形態】
【００２９】
以下、図面を参照して本発明を放射線画像装置に実施した場合の形態を詳細に説明する。図１は本発明の放射線画像処理装置１の好ましい実施の形態を示す概略図である。なお、図１のような放射線画像処理装置１の構成は、補助記憶装置に読み込まれた放射線画像処理プログラムをコンピュータ（たとえばパーソナルコンピュータ等）上で実行することにより実現される。また、この放射線画像処理プログラムは、ＣＤ−ＲＯＭ等の情報記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされることになる。
【００３０】
放射線画像処理装置１は、画像取得手段１０、画像縮小手段１５、ノイズ画像抽出手段３０、画像拡大手段４０、画像補正手段５０を備えている。画像取得手段１０は、たとえば図２に示すような胸部撮影をしたときの放射線画像Ｐ０やその他マンモグラフィ装置により取得された放射線画像Ｐ０をデジタルデータとして取得するものである。
【００３１】
画像縮小手段１５は、放射線画像Ｐ０を主方向または副方向に縮小して縮小画像Ｐ１を生成するものである。たとえば図３は放射線画像Ｐ０を副方向（矢印Ｙ方向）に縮小した場合について例示している。なお、画像縮小手段１５は予め設定された縮小率により放射線画像Ｐ０を平均縮小する。
【００３２】
図１の前処理手段２０は、縮小画像Ｐ１に対し線状ノイズに直交する方向について画素値の急激な変化を緩やかにする前処理を施すものであって、差分画像生成手段２１、差分画像調整手段２２、画像加算手段２３を備えている。差分画像生成手段２１は、抽出する線状ノイズ成分ＬＮを抽出する方向に対し直交する方向に隣接した画素間の差分を算出することにより差分画像ＳＰ１を生成するものである。たとえば副方向（矢印Ｙ方向）の線状ノイズ成分ＬＮを抽出・除去するとき、差分画像生成手段２１は図４Ａ（図４〜７は図３のＨ−Ｈ線の抜き出した画素値を示す）に示す縮小画像Ｐ１について、主方向（矢印Ｘ方向）に隣接する画素値の差分を算出することにより、図４Ｂに示すような差分値を画素値とする差分画像ＳＰ１を生成する。
【００３３】
差分画像調整手段２２は、差分画像ＳＰ１内の画素値が設定画素しきい値（ｃｌｉｐ値）以上であるとき、その画素値を所定の設定値（たとえば０）に調整し調整差分画像ＳＰ２を生成するものである（ｃｌｉｐ処理）。すなわち、画素値が急激に変化している部分は図４Ｂの部位ＭＣのように差分画像ＳＰ１における画素値（差分値）が大きくなっている。差分画像調整手段２２は当該設定画素しきい値よりも大きい画素値を所定の値に調整して図５に示すように部位ＭＣの画素値をたとえば０にする。ここで、ＳＰ１において当該設定画素しきい値よりも大きい画素を含んで単調増加または単調減少する画素も同様に所定の値に調整することにより、急激に変化する領域全体を緩やかな変化にすることができる。
【００３４】
図１の画像加算手段２３は差分画像調整手段２２において調整された調整差分画像ＳＰ２を用いて画素値を累積加算することにより前処理済みの縮小画像Ｐ２を生成するものである。具体的には、画像加算手段２３は、差分画像ＳＰ１を生成する前の放射線画像Ｐ０を復元するように、図５の調整差分画像ＳＰ２について主方向（矢印Ｘ方向）に隣接する画素値同士の加算を繰り返し行い、図６に示すような前処理済の縮小画像Ｐ２を生成する。このとき、差分画像調整手段２２においてｃｌｉｐ処理が施されているため、縮小画像Ｐ２は元の縮小画像Ｐ１の主方向（矢印Ｘ方向）における画素値が急激に変化している部分が緩やかになったものとなる。
【００３５】
図１のノイズ画像抽出手段３０は、縮小画像Ｐ１にフィルタ処理を施すことにより矢印Ｙ方向の線状ノイズ成分ＬＮを抽出するものであって、周波数処理手段３１、ノイズ画像生成手段３２を備えている。図１の周波数処理手段３１は、前処理手段２０により前処理が施された縮小画像Ｐ２に対し直交方向（矢印Ｘ方向）の高周波成分ＨＰを抽出するものである。具体的には、周波数処理手段３１は、縮小画像Ｐ２にローパスフィルタ処理を施して図７Ａに示すような低周波成分を抽出し、抽出した低周波成分を縮小画像Ｐ２から除去することにより、図７Ｂに示すような高周波成分ＨＰを抽出する。このとき、前処理手段２０による前処理が行われているため、画素値が急激に変化する部位においてリンギングが発生するのを防止することができる。
【００３６】
図１のノイズ画像生成手段３２は、周波数処理手段３１により抽出された高周波抽出画像に対し縮小方向（矢印Ｙ方向）の低周波成分を抽出することにより、図８に示す線状ノイズ画像ＮＰ１を生成するものである。つまり、線状ノイズ画像ＮＰ１は、図９に示すように、フーリエ空間において副方向の低周波数領域であって主方向の高周波数領域の成分を抽出した画像となる。
【００３７】
図１の画像拡大手段４０は、ノイズ画像生成手段３２により生成された線状ノイズ画像ＮＰ１を放射線画像Ｐ０と同一の画像サイズまで拡大するものである。つまり、画像拡大手段４０は、画像縮小手段１５における縮小率に対応した拡大率で線状ノイズ画像ＮＰ１を拡大して図１０に示すような線状ノイズ画像ＮＰ０を生成する。
【００３８】
画像補正手段５０は、画像拡大手段４０により拡大された線状ノイズ画像ＮＰ０を用いて放射線画像Ｐ０から線状ノイズを除去する補正を行うものである。具体的には、画像補正手段５０は、放射線画像Ｐ０と線状ノイズ画像ＮＰ０との差分を算出して、図１１に示す補正後の放射線画像Ｐ１０として取得する。
【００３９】
図１２は本発明の放射線画像処理方法の好ましい実施の形態を示すフローチャートであり、図１から図１２を参照して放射線画像処理方法について説明する。まず、画像取得手段１０において放射線画像Ｐ０が取得され（ステップＳＴ１、図２参照）、画像縮小手段１５において放射線画像Ｐ０を縮小した縮小画像Ｐ１が生成される（ステップＳＴ２、図３参照）。前処理手段２０において取得した放射線画像Ｐ１に対し前処理が施される（ステップＳＴ３、図４〜図６参照）。
【００４０】
次に、ノイズ画像抽出手段３０において線状ノイズ成分ＬＮを抽出した線状ノイズ画像の生成が行われる（ステップＳＴ４、図７〜図１０参照）。具体的には、前処理された縮小画像Ｐ２から矢印Ｘ方向の高周波成分ＨＰが抽出された後（図７Ａ、７Ｂ参照）、抽出された高周波成分ＨＰの矢印Ｙ方向の低周波成分が抽出されることにより線状ノイズ画像ＮＰ１が抽出される（図８、図９参照）。その後、画像拡大手段４０において線状ノイズ画像ＮＰ1が放射線画像Ｐ０の画像サイズまで拡大されることにより線状ノイズ画像ＮＰ０が作成され（ステップＳＴ５、図１０参照）、元の放射線画像Ｐ０から線状ノイズ画像ＮＰ０を除去することにより、補正後の放射線画像Ｐ１０が取得される（ステップＳＴ６、図１１参照）。
【００４１】
このように、放射線画像Ｐ０を縮小した縮小画像Ｐ１を用いて線状ノイズ画像ＮＰ０を抽出することにより、線状ノイズ成分ＬＮの抽出精度を維持しながら処理速度の向上を図ることができる。すなわち、放射線画像Ｐ０にフィルタ処理を施して線状ノイズ成分の抽出を行うとき、放射線画像Ｐ０の画像サイズに合わせたサイズのフィルタを用いてフィルタ処理を行う必要があり、フィルタ処理に時間が掛かってしまう。そこで、線状ノイズ画像ＮＰ０の抽出を縮小画像Ｐ１から行うことにより、ノイズ画像抽出手段３０内の線状ノイズ方向のフィルタ処理をサイズの小さいフィルタを用いて行うことができるため、処理時間の短縮を図ることができる。ここで、放射線画像Ｐ０からの線状ノイズ成分ＬＮの抽出結果と、縮小画像Ｐ１からの線状ノイズ成分ＬＮの抽出結果とがほぼ同一の結果になる。
【００４２】
図１３は本発明の放射線画像処理装置におけるノイズ画像抽出手段の第２の実施形態を示すブロック図であり、図１３を参照してノイズ画像抽出手段１３０について説明する。なお、図１３のノイズ画像抽出手段１３０において図１のノイズ画像抽出手段３０と同一の構成を有する部位には同一の符号を付してその説明を省略する。図１３のノイズ画像抽出手段１３０が図１のノイズ画像抽出手段３０と異なる点は、メディアン値を用いて線状ノイズ画像ＮＰ１０を生成することである。
【００４３】
ノイズ画像抽出手段１３０は、周波数処理手段３１、ライン分割手段１３１、メディアン算出手段１３２、ノイズ画像生成手段１３３を備えている。ライン分割手段１３１は、周波数処理手段３１により抽出された高周波成分ＨＰについて矢印Ｙ方向の各画素ラインを複数の分割画素ラインＢＬ１〜ＢＬ３に分割するものである。具体的には、図１４に示すように、副方向（矢印Ｙ方向）の線状ノイズを抽出する際、ライン分割手段１３１は副方向（矢印Ｙ方向）の画素ラインＬをたとえば３等分して１画素ラインについて３つの分割画素ラインＢＬ１〜ＢＬ３を生成する。
【００４４】
このとき、ライン分割手段１３１は分割画素ラインＢＬ１〜ＢＬ３が生体構造に存在し得る直線の２倍以上の長さからなるように分割する。生体構造に存在し得る直線の２倍以上の長さとは経験的もしくは統計的に求められるものであって、予め所定の長さに設定される。
【００４５】
つまり分割画素数は、スジ除去の処理を行った放射線画像Ｐ１０が診断に問題がなく、且つ、視認されるスジが残らないように最適化されるパラメータである。言い換えれば、補正後の放射線画像Ｐ１０において、診断等において問題になるような線状ノイズ成分ＬＮを除去しながら、生体構造が線状ノイズ成分ＬＮとして抽出されることによる画質の劣化を最小限に抑えることができるように最適化する。ここで、最適な分割画素数は除去するスジの太さによって変わる。ここでは分割画素ラインＢＬ１〜ＢＬ３が生成される場合を示しているが、分割数は分割画素数によって変わり、分割が１、つまり分割されない場合もある。
【００４６】
図１３のメディアン算出手段１３２は、ライン分割手段１３１により分割された各分割画素ラインＢＬ１〜ＢＬ３毎にそれぞれメディアン値Ｍ１〜Ｍ３を算出するものである。つまり、メディアン算出手段１３２は分割画素ラインＢＬ１〜ＢＬ３の画素数を計算幅としてメディアン値Ｍを算出する。このとき、分割画素ラインＢＬ１〜ＢＬ３は、生体構造に存在し得る直線の２倍以上の長さからなっているため、メディアン値（中間値）Ｍ１〜Ｍ３として生体構造を示す画素値が算出されることがない。たとえば図１４に示すように、線状ノイズ成分ＬＮが副方向（矢印Ｙ方向）の画素ラインＬの一部に生体構造を示す画素値（濃度）が存在した場合であっても、図１５のように生体構造を示す画素値が分割画素ラインＢＬ２内の複数の画素値の中間値にならず、メディアン値Ｍ２として算出されることはない。つまり、メディアン算出手段１３２は線状ノイズ成分ＬＮを示す画素値をメディアン値Ｍ２として算出することになる。
【００４７】
図１３のノイズ画像生成手段１３３は、メディアン算出手段１３２により算出された各分割画素ラインＢＬ１〜ＢＬ３毎のメディアン値Ｍ１〜Ｍ３を用いて線状ノイズ成分ＬＮを示す線状ノイズ画像ＮＰ１１を生成するものである。具体的には、図１６に示すように、ノイズ画像生成手段１３３は、各分割画素ラインＢＬ１〜ＢＬ３の中心の画素位置の画素値をメディアン値Ｍ１〜Ｍ３とし、複数のメディアン値Ｍ１、Ｍ２およびメディアン値Ｍ２、Ｍ３をそれぞれ直線で結ぶ。そして、ノイズ画像生成手段１３３は副方向（矢印Ｙ方向）の各画素に対する画素値を示す画素関数Ｆ１（ｙ）、Ｆ２（ｙ）（ｙは副方向の画素）を画素ラインＬ毎に生成する。その後、ノイズ画像生成手段１３３は生成した画素関数Ｆ１（ｙ）、Ｆ２（ｙ）を用いて各画素毎に画素値を算出し線状ノイズ画像ＮＰ１１を生成する。
【００４８】
なお、図１６において、各メディアン値Ｍ１、Ｍ２、Ｍ３をそれぞれ直線で結ぶ関数を生成する場合について例示しているが、これに限定されず、各メディアン値Ｍ１〜Ｍ３を１つの曲線で結ぶ関数Ｆ（ｙ）を生成し、線状ノイズ画像ＮＰ１１を生成するものであってもよい。
【００４９】
このように、メディアン値Ｍ１〜Ｍ３を用いて線状ノイズ成分ＬＮの抽出を行う場合であっても、縮小画像Ｐ２を用いて行うことによりメディアン計算幅を小さくすることができるため、線状ノイズ成分ＬＮの検出精度を維持しつつ処理時間の短縮化を図ることができる。さらに、分割画素ラインＢＬ１〜ＢＬ３が人体構造に含まれる直線の長さの２倍以上であるため、線状ノイズ成分ＬＮとして人体構造が含まれることが無く、画像補正手段５０において放射線画像を補正するときに人体構造が除去されるのを防止することができる。
【００５０】
図１７は本発明の放射線画像処理装置におけるノイズ画像抽出手段の第３の実施形態を示すブロック図であり、図１７を参照してノイズ画像抽出手段２３０について説明する。なお、図１７のノイズ画像抽出手段２３０において図１および図１３のノイズ画像抽出手段３０、１３０と同一の構成を有する部位には同一の符号を付してその説明を省略する。図１７のノイズ画像抽出手段２３０が図１、図１３のノイズ画像抽出手段３０、１３０と異なる点は、図１のフィルタ処理により生成された第１線状ノイズ画像ＮＰ１と、図１２のメディアン値Ｍ１〜Ｍ３を用いて生成された第２線状ノイズ画像ＮＰ１１とを用いて補正に用いる線状ノイズ画像ＮＰを生成することである。
【００５１】
具体的には、図１７のノイズ画像抽出手段２３０は、周波数処理手段３１、第１ノイズ画像生成手段３２、ライン分割手段１３１、メディアン算出手段１３２、第２ノイズ画像生成手段１３３、ノイズ画像作成手段２３１を備えている。上述したように第１ノイズ画像生成手段３２は第１線状ノイズ画像ＮＰ１を生成し（図８参照）、第２ノイズ画像生成手段１３３は第２線状ノイズ画像ＮＰ１１を生成する（図１４〜１６参照）。そして、ノイズ画像作成手段２３１は、第１線状ノイズ画像ＮＰ１および第２線状ノイズ画像ＮＰ１１を用いて補正に用いる線状ノイズ画像ＮＰを生成する。
【００５２】
ここで、ノイズ画像作成手段２３１は、第１線状ノイズ画像ＮＰ１と第２線状ノイズ画像ＮＰ１１との各画素の絶対値を比較し、絶対値の小さい値を用いることにより線状ノイズ画像ＮＰを作成する。つまり、一般的に線状ノイズ成分ＬＮの濃度（絶対値）は小さいため、絶対値の小さい方の値がより線状ノイズ成分ＬＮであると推定することができる。そこで、２つの線状ノイズ画像ＮＰ１、ＮＰ１１の各画素における絶対値を比較し、絶対値の小さい方を用いて補正に用いる線状ノイズ画像ＮＰを生成することにより、アーチファクトを低減させることができる。
【００５３】
もしくは、ノイズ画像作成手段２３１は、放射線が曝射された被写体の種類に合わせて第１線状ノイズ画像ＮＰ１もしくは第２線状ノイズ画像ＮＰ１１のいずれか一方を選択することにより線状ノイズ画像ＮＰを作成するものであってもよい。つまり、放射線画像の性質によって、フィルタ処理もしくはメディアン処理を用いた方が適している放射線画像がある。たとえば短い線状成分が多く含まれている被写体を撮影した放射線画像Ｐ０については第１線状ノイズ画像ＮＰ１を用いて補正を行い、テスト画像で用いられる幾何学パターンを撮影した放射線画像Ｐ０については第２線状ノイズ画像ＮＰ１１を用いて補正を行うようにしてもよい。
【００５４】
さらに、ノイズ画像作成手段２３１は、１つの線状ノイズ画像ＮＰを生成する場合について例示しているが、上記２つの線状ノイズ画像ＮＰ１、ＮＰ１１をそのまま補正に用いる線状ノイズ画像ＮＰとして用い、２つの線状ノイズ画像ＮＰ１、ＮＰ１１を用いて２つの補正済みの放射線画像Ｐ１０を生成するようにしてもよい。
【００５５】
図１８は本発明の放射線画像処理装置におけるノイズ画像抽出手段の第４の実施形態を示すブロック図であり、図１８を参照してノイズ画像抽出手段３３０について説明する。なお、図１８のノイズ画像抽出手段３３０において図１７のノイズ画像抽出手段２３０と同一の構成を有する部位には同一の符号を付してその説明を省略する。図１８のノイズ画像抽出手段３３０が図１７のノイズ画像抽出手段２３０と異なる点は、縮小画像Ｐ２を各周波数成分毎に分けて当該周波数成分毎に異なる処理をすることにより各種線状ノイズ画像を生成する点である。
【００５６】
図１８の周波数処理手段３３１は、矢印Ｘ方向の高周波成分ＨＰを抽出するとともに低周波成分ＬＰをも抽出する機能を備えている。そして、第１ノイズ画像生成手段３２は高周波成分ＨＰを用いて第１線状ノイズ画像ＮＰ１を生成する。一方、ライン分割手段１３１、メディアン算出手段１３２、第２ノイズ画像生成手段１３３は低周波成分ＬＰを用いて第２線状ノイズ画像ＮＰ１１を生成する。なお、この場合であっても、ノイズ画像作成手段２３１は、図１７において説明したように、２つの線状ノイズ画像ＮＰ１、ＮＰ１１を用いてノイズ画像を作成する。画像拡大手段４０でノイズ画像を拡大し、画像補正手段５０でオリジナル画像Ｐ０より引くことで補正済みの放射線画像Ｐ１０を生成する。
【００５７】
このように、縮小画像Ｐ２から各周波数成分を抽出し、各周波数成分毎に異なる線状ノイズ成分ＬＮの抽出処理を行うことにより、各周波数成分における線状ノイズ成分ＬＮの特性に応じて線状ノイズ成分ＬＮの抽出処理を行うことができるため、アーチファクトの低減を図ることができる。つまり、高周波成分に対しては濃度がある程度変動していても線状ノイズ成分ＬＮとして抽出が可能なフィルタ処理により第１線状ノイズ画像ＮＰ１を生成するとともに、低周波成分に対しては人体構造を多く含むことを考慮し、短い線状成分を線状ノイズ成分ＬＮとして抽出するのを防止することができるメディアン値Ｍ１〜Ｍ３を用いた第２線状ノイズ画像の生成を行う。したがって、アーチファクトの低減を図ることができる。
【００５８】
上記各実施の形態によれば、放射線画像Ｐ０を垂直方向または水平方向に縮小した縮小画像Ｐ１を生成し、生成した縮小画像Ｐ１（Ｐ２）の中から矢印Ｙ方向の線状ノイズ成分ＬＮを線状ノイズ画像ＮＰとして抽出し、抽出した線状ノイズ画像ＮＰを放射線画像Ｐ０の画像サイズに拡大し、拡大した線状ノイズ画像ＮＰを用いて放射線画像Ｐ０から線状ノイズ成分ＬＮを除去することにより、線状ノイズ成分ＬＮは線状ノイズが延びる方向への濃度変化が少なく、縮小画像Ｐ１から抽出した線状ノイズ成分ＬＮと元の放射線画像Ｐ０から抽出した線状ノイズ成分ＬＮとはほぼ同等の抽出結果になる知見に基づき、抽出する方向に放射線画像Ｐ０を縮小した縮小画像Ｐ１を用いて線状ノイズ成分ＬＮを抽出するものであるため、線状ノイズ成分ＬＮの抽出精度を維持しながら処理時間の短縮を図ることができる。
【００５９】
また、図１に示すように、ノイズ画像抽出手段３０が、縮小画像Ｐ１（Ｐ２）から矢印Ｘ方向の高周波成分ＨＰを抽出する周波数処理手段３１と、周波数処理手段３１により高周波成分ＨＰから線状ノイズに直交する方向の低周波成分を抽出することにより線状ノイズ画像ＮＰ１を生成するノイズ画像抽出手段３２とを有するものであるとき、縮小画像Ｐ１の画像サイズに合わせたフィルタサイズによりフィルタ処理を行うことができるため、処理時間の短縮を図ることができる。
【００６０】
さらに、図１３〜図１６に示すように、ノイズ画像抽出手段１３０が、縮小画像Ｐ１（Ｐ２）から矢印Ｘ方向の高周波成分ＨＰを抽出する周波数処理手段３１と、周波数処理手段３１により抽出された高周波成分ＨＰについて矢印Ｘ方向の画素ラインＬを複数の分割画素ラインＢＬ１〜ＢＬ３に分割するライン分割手段１３１と、ライン分割手段１３１により分割された各分割画素ラインＢＬ１〜ＢＬ３のメディアン値Ｍ１〜Ｍ３を算出するメディアン算出手段１３２と、メディアン算出手段１３２により算出された各分割画素ラインＢＬ１〜ＢＬ３毎のメディアン値Ｍ１〜Ｍ３を用いて分割画素ラインＢＬ１〜ＢＬ３が延びる方向の線状ノイズ成分ＬＮを示す線状ノイズ画像ＮＰ１１を生成するノイズ画像生成手段１３３とを有するものであれば、メディアン計算幅を小さくして処理時間の短縮を図ることができるとともに、メディアン値Ｍ１〜Ｍ３として人体構造を示す値が算出されることを抑えることができるため、人体のパターンの中にスジ方向に伸びている部位がある場合、その部位も線状ノイズとして除去されるのを防止し線状ノイズ成分ＬＮのみ除去し、画像補正による画質の劣化が生じるのを防止することができる。
【００６１】
さらに、図１７に示すように、ノイズ画像抽出手段２３０が、縮小画像Ｐ１（Ｐ２）から矢印Ｘ方向の高周波成分ＨＰを抽出する周波数処理手段３１と、周波数処理手段３１において抽出された高周波成分ＨＰから矢印Ｙ方向の低周波成分を抽出することにより第１線状ノイズ画像ＮＰ１を生成する第１ノイズ画像生成手段３２と、周波数処理手段３１により抽出された高周波成分ＨＰについて矢印Ｙ方向の画素ラインを複数の分割画素ラインＢＬ１〜ＢＬ３に分割するライン分割手段１３１と、ライン分割手段により分割された各分割画素ラインのメディアン値Ｍ１〜Ｍ３を算出するメディアン算出手段１３２と、メディアン算出手段１３２により算出された各分割画素ラインＢＬ１〜ＢＬ３毎のメディアン値Ｍ１〜Ｍ３を用いて分割画素ラインＢＬ１〜ＢＬ３が延びる方向の線状ノイズ成分ＬＮを示す第２線状ノイズ画像ＮＰ１１を生成する第２ノイズ画像生成手段１３３とを備えたものであり、ノイズ画像作成手段２３１が第１線状ノイズ画像ＮＰ１と第２線状ノイズ画像ＮＰ１１とを用いて線状ノイズ画像ＮＰ２０を作成するものであるとき、線状ノイズの特性に合わせた最適な処理で抽出した２種類の第１線状ノイズ画像と第２線状ノイズ画像とを用いて線状ノイズ画像が生成されるため、補正後の放射線画像におけるアーチファクトの発生を低減することができる。
【００６２】
また、図１８に示すように、ノイズ画像抽出手段３３０が、縮小画像Ｐ１（Ｐ２）から矢印Ｙ方向の高周波成分ＨＰと低周波成分ＬＰとを抽出する周波数処理手段３３１と、周波数処理手段３３１において抽出された高周波成分ＨＰから矢印Ｙ方向低周波成分を抽出することにより第１線状ノイズ画像ＮＰ１を生成する第１ノイズ画像生成手段３２と、周波数処理手段３３１において抽出された低周波成分ＬＰについて線状ノイズに沿う方向（矢印Ｙ方向）の画素ラインを複数の分割画素ラインＢＬ１〜ＢＬ３に分割するライン分割手段１３１と、ライン分割手段１３１により分割された各分割画素ラインＢＬ１〜ＢＬ３のメディアン値Ｍ１〜Ｍ３を算出するメディアン算出手段１３２と、メディアン算出手段１３２により算出された各分割画素ラインＢＬ１〜ＢＬ３毎のメディアン値Ｍ１〜Ｍ３を用いて線状ノイズ成分ＬＮを示す第２線状ノイズ画像ＮＰ１１を生成する第２ノイズ画像生成手段１３３とを備えたものであり、ノイズ画像作成手段２３１が第１線状ノイズ画像ＮＰ１と第２線状ノイズ画像ＮＰ１１とを用いて線状ノイズ画像ＮＰ２０を作成するものであれば、放射線画像内の人体構造を示す成分は低周波成分が多いことを利用し、線状ノイズ成分ＬＮを確実に除去するとともに補正によるアーチファクトの発生を低減することができる。
【００６３】
また、ノイズ画像作成手段２３１が、第１線状ノイズ画像ＮＰ１と第２線状ノイズ画像ＮＰ１１との各画素の絶対値を比較し、絶対値の小さい値を用いることにより線状ノイズ画像ＮＰ２０を作成するものであれば、線状ノイズである確率が高い成分を用いて放射線画像Ｐ０を補正することができるため、アーチファクトを低減させることができる。
【００６４】
さらに、ノイズ画像作成手段２３１が、放射線画像の性質に合わせて第１線状ノイズ画像ＮＰ１もしくは第２線状ノイズ画像ＮＰ１１のいずれかを選択して放射線画像Ｐ０を補正するものであるとき、たとえば幾何学パターンが描かれた放射線画像Ｐ０の場合には第２線状ノイズ画像ＮＰ１１を用いて補正することにより確実に線状ノイズ成分ＬＮを除去し、放射線画像Ｐ０内に人体構造を示す短い線状成分が複数含まれている場合には第１線状ノイズ画像ＮＰ１を用いて補正することによりアーチファクトを低減することができる。
【００６５】
本発明の実施形態は、上記実施形態に限定されない。たとえば、上記各実施形態において、実際に被写体が撮影された放射線画像Ｐ０に対する補正を行う場合について例示しているが、たとえばオフセット補正、シェーディング補正、残像補正等を行うために必要な補正用画像（たとえばベタ画像）についても上述した線状ノイズ成分ＬＮを除去する補正を適用することができる。なお、補正用画像には、人体構造に含まれる直線成分のような線状ノイズ成分ＬＮ以外の線状のデータは含まれていないため、強力なスジ除去を行うようなフィルタパラメータ等を用いるようにしてもよい。
【００６６】
さらに、上記各実施形態において、放射線画像Ｐ０全体について各種処理を行う場合について例示しているが、たとえばマンモグラフィ撮影により取得された放射線画像は、図１９に示すようにいわゆる素抜け領域が存在することになる。そこで、画像取得手段１０において素抜け領域を抽出し、素抜け領域以外の実データ領域についてのみ上記各種処理を行うようにしてもよい。
【００６７】
以上、放射線画像について述べてきたが、可視光画像なども含めた2次元画像全てに適用可能であることは勿論である。
【図面の簡単な説明】
【００６８】
【図１】本発明の放射線画像処理装置の好ましい実施の形態を示すブロック図
【図２】図１の画像取得手段において取得される放射線画像の一例を示す図
【図３】図１の画像縮小手段において放射線画像が縮小された様子を示す図
【図４Ａ】図１の前処理手段において前処理された画像の一例を示すグラフ
【図４Ｂ】図１の前処理手段において生成された差分画像の一例を示すグラフ
【図５】図１の前処理手段において調整が行われた差分画像の一例を示すグラフ
【図６】図１の前処理手段において差分画像に対し累積加算処理が行われた様子を示すグラフ
【図７Ａ】図１の周波数処理手段において抽出された低周波成分の一例を示すグラフ
【図７Ｂ】図１の周波数処理手段において抽出された高周波成分の一例を示すグラフ
【図８】図１のノイズ画像生成手段において生成された線状ノイズ画像の一例を示すグラフ
【図９】図８の線状ノイズ画像に含まれる周波数成分をフーリエ空間上に表したグラフ
【図１０】図１の画像拡大手段において拡大された線状ノイズ画像の一例を示す図
【図１１】図１の画像補正手段により補正された放射線画像の一例を示す図
【図１２】本発明の放射線画像処理方法の好ましい実施形態を示すフローチャート
【図１３】本発明の放射線画像処理装置におけるノイズ画像抽出手段の第２の実施形態を示すブロック図
【図１４】図１３のライン分割手段により画素ラインが分割される様子を示す模式図
【図１５】図１３のメディアン算出手段により各分割画素ライン毎にメディアン値が算出される様子を示すグラフ
【図１６】図１３のノイズ画像生成手段においてメディアン値から線状ノイズ画像を生成する際に導出される関数の一例を示すグラフ
【図１７】本発明の放射線画像処理装置におけるノイズ画像抽出手段の第３の実施形態を示すブロック図
【図１８】本発明の放射線画像処理装置におけるノイズ画像抽出手段の第４の実施形態を示すブロック図
【図１９】本発明の放射線画像処理装置に適用されるマンモグラフィー撮影により取得された放射線画像の一例を示す図
【符号の説明】
【００６９】
１、１００、２００、３００     放射線画像処理装置
１０  画像取得手段
１５  画像縮小手段
２０  前処理手段
２１  差分画像生成手段
２２  差分画像調整手段
２３  画像加算手段
３０、１３０、２３０、３３０   ノイズ画像抽出手段
３１、３３１       周波数処理手段
３２  ノイズ画像生成手段
４０  画像拡大手段
５０  画像補正手段
１３１      ライン分割手段
１３２      メディアン算出手段
１３３      ノイズ画像生成手段
２３１ノイズ画像作成手段
ＢＬ１、ＢＬ２、ＢＬ３   分割画素ライン
Ｍ１、Ｍ２、Ｍ３   メディアン値
ＮＰ  線状ノイズ画像
ＮＰ１      第１線状ノイズ画像
ＮＰ１１    第２線状ノイズ画像
Ｐ０  放射線画像
Ｐ１０      補正画像
Ｐ１、Ｐ２  縮小画像
【図１】
2010009348.tif 000003
【図２】
2010009348.tif 000004
【図３】
2010009348.tif 000005
【図４Ａ】
2010009348.tif 000006
【図４Ｂ】
2010009348.tif 000007
【図５】
2010009348.tif 000008
【図６】
2010009348.tif 000009
【図７Ａ】
2010009348.tif 000010
【図７Ｂ】
2010009348.tif 000011
【図８】
2010009348.tif 000012
【図９】
2010009348.tif 000013
【図１０】
2010009348.tif 000014
【図１１】
2010009348.tif 000015
【図１２】
2010009348.tif 000016
【図１３】
2010009348.tif 000017
【図１４】
2010009348.tif 000018
【図１５】
2010009348.tif 000019
【図１６】
2010009348.tif 000020
【図１７】
2010009348.tif 000021
【図１８】
2010009348.tif 000022
【図１９】
2010009348.tif 000023