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