(57)【要約】  【課題】  画像からのスジムラ除去において、リンギングの発生を抑制し、また、構造物をスジムラとして除去することも抑制でき画像補正方法および装置を提供する。  【解決手段】  元画像にスジムラと直交する方向にメディアン処理を施し、この画像を元画像から減算し、減算して得られた画像から、濃度が閾値を超える領域の濃度を低減し、この濃度低減処理を行なった画像に、スジムラ方向にローパスフィルタによる周波数処理を行い、周波数処理を施した画像を、元画像から減算することにより、前記課題を解決する。【選択図】図３
2010239480.tif 000002
【特許請求の範囲】
【請求項１】
直線状のムラを有する元画像から、この直線状のムラを除去するに際し、
前記元画像に、前記直線状のムラの延在方向と直交する方向にメディアン処理を施してメディアン処理画像を生成し、このメディアン処理画像を前記元画像から減算して減算画像を生成し、この減算画像から、濃度が所定の閾値を超える部分に対して、濃度を低減する閾値処理を行なって閾値処理画像を生成し、この閾値処理画像に、前記直線状のムラの延在方向と同方向に一次元のローパスフィルタによる周波数処理を施して処理用画像を生成し、
前記元画像から、この処理用画像を減算することにより、前記元画像から前記直線状のムラを除去することを特徴とする画像補正方法。
【請求項２】
前記メディアン処理は、前記直線状のムラの延在方向と直交する方向の画素数に対して、２倍以上の画素のマスクを有する１次元フィルタを用いて行なう請求項１に記載の画像補正方法。
【請求項３】
前記フィルタのマスクの画素数が奇数である請求項２に記載の画像補正方法。
【請求項４】
前記フィルタが、前記直線状のムラの延在方向と直交する方向の画素数ｎに対して、［２ｎ＋１］画素のマスクを有するフィルタである請求項３に記載の画像補正方法。
【請求項５】
前記元画像が、二次元的に配列された画素配列の一方向に信号を読み取る撮像素子によって得られた画像であり、前記メディアン処理を、前記読取方向と直交する方向に行なう請求項１〜４のいずれかに記載の画像補正方法。
【請求項６】
前記元画像が、一方向に配列した読取画素を有するラインセンサによって、前記画素列と直交する方向に、前記ラインセンサと読取対象とを走査して読み取られた画像であり、前記メディアン処理を、前記走査方向と直交する方向に行なう請求項１〜４のいずれかに記載の画像補正方法。
【請求項７】
前記元画像が、静止グリッドを介して撮影された放射線画像であり、前記メディアン処理を、前記静止グリッドの延在方向と直交する方向に行なう請求項１〜４のいずれかに記載の画像補正方法。
【請求項８】
直線状のムラを有する元画像に、この直線状のムラの延在方向と直交する方向にメディアン処理を施すメディアン処理手段と、
前記元画像から、前記メディアン処理手段で処理した画像を減算する第１減算手段と、
前記第１減算手段によって得られた画像中の濃度が所定の閾値を超える部分に対して、濃度を低減する処理を行なう閾値処理手段と、
前記閾値処理手段によって得られた画像に対して、前記直線状のムラの延在方向に、一次元のローパスフィルタによる周波数処理を施すフィルタ処理手段と、
元画像から、前記フィルタ処理手段で処理した画像を減算する第２減算手段とを有することを特徴とする画像補正装置。
【請求項９】
前記メディアン処理手段は、前記直線状のムラの延在方向と直交する方向の画素数に対して、２倍以上の画素のマスクを有する１次元フィルタを用いて元画像を処理する請求項８に記載の画像補正装置。
【請求項１０】
前記フィルタのマスクの画素数が奇数である請求項９に記載の画像補正装置。
【請求項１１】
前記フィルタが、前記直線状のムラの延在方向と直交する方向の画素数ｎに対して、［２ｎ＋１］画素のマスクを有するフィルタである請求項１０に記載の画像補正装置。
【請求項１２】
二次元的に配列された画素配列の一方向に信号を読み取る撮像素子によって得られた画像を前記元画像とする際に、前記メディアン処理手段は、前記読取方向と直交する方向にメディアン処理を行なう請求項８〜１１のいずれかに記載の画像補正装置。
【請求項１３】
一方向に配列した読取画素を有するラインセンサによって、前記画素列と直交する方向に、前記ラインセンサと読取対象とを走査することにより読み取られた画像を前記元画像とする際に、前記メディアン処理手段は、前記走査方向と直交する方向にメディアン処理を行なう請求項８〜１１のいずれかに記載の画像補正装置。
【請求項１４】
静止グリッドを介して撮影された放射線画像を前記元画像とする際に、前記メディアン処理手段は、前記静止グリッドの延在方向と直交する方向にメディアン処理を行なう請求項８〜１１いずれかに記載の画像補正装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像に生じてしまったスジ状（直線状）のムラの補正に関し、詳しくは、リンギングを抑制し、かつ、必要な構造物の除去も防止して、画像のスジ状のムラを補正できる画像補正方法および画像補正装置に関する。
【背景技術】
【０００２】
従来より、医療用の診断画像の撮影や工業用の非破壊検査などに、被写体を透過した放射線（Ｘ線、α線、β線、γ線、電子線、紫外線等）を電気的な信号として取り出すことにより放射線画像を撮影する、放射線画像検出器が利用されている。
この放射線画像検出器としては、放射線を電気的な画像信号として取り出す放射線固体検出器（いわゆる「Flat Panel Detector」  以下、ＦＰＤとする）が知られている。
【０００３】
このＦＰＤには、例えば放射線の入射によってアモルファスセレンなどの光導電膜が発した電子−正孔対（ｅ−ｈペア）を収集して電気信号として読み出す、いわば放射線を直接的に電気信号に変換する直接方式と、放射線の入射によって発光（蛍光）する蛍光体で形成された蛍光体層（シンチレータ層）を有し、この蛍光体層によって放射線を可視光に変換し、この可視光を光電変換素子で読み出す、いわば放射線を可視光として電気信号に変換する間接方式との、２つの方式がある。
【０００４】
さらに、ＦＰＤではなく、放射線の照射を受けると、この放射線エネルギーの一部を蓄積し、その後、可視光等の励起光の照射を受けると、蓄積されたエネルギーに応じた輝尽発光を示す蛍光体（輝尽性蛍光体）の膜を有する放射線画像変換パネル（ＩＰ(imaging plate)）を利用するシステムも知られており、例えば、富士フイルム社製のＦＣＲ（Fuji Computed Radiography)等として実用化されている。
【０００５】
ＦＰＤでは、放射線画像の撮影条件等によっては、放射線の照射に応じた電気信号の読み出し方法に応じて、画素列の方向に延在する（長尺な）、直線状のムラ（スジムラ）が生じてしまう場合が有る。
さらに、前記ＩＰに撮影した放射線画像をラインセンサで読み取る場合に、ラインセンサに不良画素が生じた場合には、読み取りの走査方向にスジムラが生じてしまう。
【０００６】
また、このようなＦＰＤやＩＰを用いる放射線画像の撮影においては、被写体によって散乱された散乱光が撮影媒体に入射しないように、放射線を透過しない鉛などの板と、放射線を透過し易いアルミニウムなどの板が交互に配置された静止グリッドを用いて撮影を行なうことが有る。このような静止グリッドを用いて放射線画像を撮影すると、散乱光の抑制によって放射線画像のコントラストを向上できる反面、この静止グリッドに対応して、静止グリッドの鉛板の長手方向に延在するスジムラが放射線画像に乗ってしまう。
【０００７】
このような放射線画像のスジムラを除去する方法として、本件出願人にかかる特許文献１に開示される方法が知られている。この方法は、元画像（画像信号（画像データ））に、スジムラ（静止グリッドに起因する周期的パターン）と同一方向の一次元フィルタによる空間周波数処理、および、この方向と直交する方向の一次元フィルタによる空間周波数処理を行うことにより、元画像のスジムラを抽出し、元画像からスジムラを差し引くことにより、元画像からスジムラを除去するものである。
具体的には、一例として、元画像に対して、主走査方向（スジムラと直交する方向）に一次元のローパスフィルタ処理を行って低周波成分を取り出し、この低周波成分を元画像から差し引くことにより、元画像から高周波数成分を抽出し、さらに、抽出した高周波数成分に対して、副走査方向に一次元のローパスフィルタ処理を行なうことにより、スジムラのみを抽出し、このスジムラを元画像から差し引くことにより、元画像からスジムラを除去する。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００３−１５０９５４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
このようなローパスフィルタ（ハイパスフィルタ）を用いた処理を行なうことにより、元画像からスジムラを抽出して、取り除くことができる。
しかしながら、スジムラと直交する方向にローパスフィルタ等によるフィルタ処理を施すと、濃度変化の大きいスジムラやエッジの境界のように、フィルタ処理の方向で画像データの値が急激に変化する領域において、リンギングが発生してしまう。
【００１０】
例えば、元画像が、図４（Ａ）に示されるようなスジムラを有する場合に、スジムラを含んで、スジムラと直交する領域Ａのプロファイルを取ると、そのプロファイルは、図４（Ｂ）の上段に実線で示すような状態となる。
このような領域Ａに対して、スジムラと直交する方向にローパスフィルタで処理すると、図４（Ｂ）の上段に破線で示すような画像が得られる。このローパスフィルタで処理した画像（破線）を、元画像（実線）から減算すると、図４（Ｂ）の下段に示すように、濃度が急激に変動するスジムラの両側が異常に下がってしまうリンギングが発生し、此処にアーチファクトが生じてしまう。
【００１１】
また、ローパスフィルタ等で処理した画像を用いて元画像を処理するスジムラの抽出方法では、元画像中に、スジムラと平行な構造物（被写体の絵柄）が存在した場合には、この構造物（その画像信号）もスジムラとして検出してしまい、スジムラの除去によって、スジムラのみならず、元画像中の必要な構造物を除去してしまう可能性も有る。
【００１２】
本発明の目的は、前記従来技術の問題点を解決することにあり、放射線画像等の各種の画像（画像データ）に存在するスジムラを、リンギングの発生を十分に抑制して除去することができ、かつ、画像中の必要な構造物を除去することも防止できる画像補正方法、および、この画像補正方法を実施する画像補正装置を提供することにある。
【課題を解決するための手段】
【００１３】
このような目的を達成するために、本発明の画像補正方法は、直線状のムラを有する元画像から、この直線状のムラを除去するに際し、前記元画像に、前記直線状のムラの延在方向と直交する方向にメディアン処理を施してメディアン処理画像を生成し、このメディアン処理画像を前記元画像から減算して減算画像を生成し、この減算画像から、濃度が所定の閾値を超える部分に対して、濃度を低減する閾値処理を行なって閾値処理画像を生成し、この閾値処理画像に、前記直線状のムラの延在方向と同方向に一次元のローパスフィルタによる周波数処理を施して処理用画像を生成し、前記元画像から、この処理用画像を減算することにより、前記元画像から前記直線状のムラを除去することを特徴とする画像補正方法を提供する。
【００１４】
このような本発明の画像補正方法において、前記メディアン処理は、前記直線状のムラの延在方向と直交する方向の画素数に対して、２倍以上の画素のマスクを有する１次元フィルタを用いて行なうのが好ましく、特に、前記フィルタのマスクの画素数が奇数であるのが好ましく、さらに、前記フィルタが、前記直線状のムラの延在方向と直交する方向の画素数ｎに対して、［２ｎ＋１］画素のマスクを有するフィルタであるのが好ましい。
また、前記元画像が、二次元的に配列された画素配列の一方向に信号を読み取る撮像素子によって得られた画像であり、前記メディアン処理を、前記読取方向と直交する方向に行なうのが好ましく、あるいは、前記元画像が、一方向に配列した読取画素を有するラインセンサによって、前記画素列と直交する方向に、前記ラインセンサと読取対象とを走査して読み取られた画像であり、前記メディアン処理を、前記走査方向と直交する方向に行なうのが好ましく、あるいは、前記元画像が、静止グリッドを介して撮影された放射線画像であり、前記メディアン処理を、前記静止グリッドの延在方向と直交する方向に行なうのが好ましい。
【００１５】
また、本発明の画像補正装置は、直線状のムラを有する元画像に、この直線状のムラの延在方向と直交する方向にメディアン処理を施すメディアン処理手段と、前記元画像から、前記メディアン処理手段で処理した画像を減算する第１減算手段と、前記減算手段によって得られた画像中の濃度が所定の閾値を超える部分に対して、濃度を低減する処理を行なう閾値処理手段と、前記閾値処理手段によって得られた画像に対して、前記直線状のムラの延在方向に、一次元のローパスフィルタによる周波数処理を施すフィルタ処理手段と、元画像から、前記フィルタ処理手段で処理した画像を減算する第２減算手段とを有することを特徴とする画像補正装置を提供する。
【００１６】
このような本発明の画像補正装置において、前記メディアン処理手段は、前記直線状のムラの延在方向と直交する方向の画素数に対して、２倍以上の画素のマスクを有する１次元フィルタを用いて元画像を処理するのが好ましく、特に、前記フィルタのマスクの画素数が奇数であるのが好ましく、さらに、前記フィルタが、前記直線状のムラの延在方向と直交する方向の画素数ｎに対して、［２ｎ＋１］画素のマスクを有するフィルタであるのが好ましい。
また、二次元的に配列された画素配列の一方向に信号を読み取る撮像素子によって得られた画像を前記元画像とする際に、前記メディアン処理手段は、前記読取方向と直交する方向にメディアン処理を行なうのが好ましく、もしくは、一方向に配列した読取画素を有するラインセンサによって、前記画素列と直交する方向に、前記ラインセンサと読取対象とを走査することにより読み取られた画像を前記元画像とする際に、前記メディアン処理手段は、前記走査方向と直交する方向にメディアン処理を行なうのが好ましく、もしくは、静止グリッドを介して撮影された放射線画像を前記元画像とする際に、前記メディアン処理手段は、前記静止グリッドの延在方向と直交する方向にメディアン処理を行なうのが好ましい。
【発明の効果】
【００１７】
上記構成を有する本発明によれば、放射線画像等の元画像からスジムラ（直線状上のムラ）を除去するに際し、ローパスフィルタやハイパスフィルタ等ではなく、メディアン処理を行い、このメディアン処理を施した画像を元画像から減算する処理を利用してスジムラを抽出する。そのため、元画像中に、濃度変化の大きいスジムラや構造物のエッジのように、濃度（画像データ）が急激に変動する領域が存在しても、リンギングの発生を十分に抑制することができ、これに起因するアーチファクト等の画質劣化を、抑制する事が出来る。
さらに、メディアン処理を施した画像に対して、所定の閾値を超える領域に対して、画像データを小さくする処理を行い、この閾値処理を行なった画像に、スジムラの延在方向にローパスフィルタをかけてスジムラを抽出し、元画像から減算する。そのため、元画像中にスジムラと平行する構造物が存在した場合でも、この構造物をスジムラとして抽出してしまうことを防止できる。
【図面の簡単な説明】
【００１８】
【図１】本発明の画像補正装置を利用する放射線画像撮影装置の一例の概念図である。
【図２】本発明の画像補正装置の一例を概念的に示すブロック図である。
【図３】（Ａ）〜（Ｅ）は、本発明によるスジムラ除去を説明するための概念図である。
【図４】（Ａ）〜（Ｃ）は、本発明および従来のスジムラ除去を説明するための概念図である。
【図５】（Ａ）および（Ｂ）は、本発明によるスジムラ除去を説明するための概念図である。
【発明を実施するための形態】
【００１９】
以下、本発明の画像補正方法および画像補正装置について、添付の図面に示される好適実施例を基に、詳細に説明する。
【００２０】
図１に、本発明の画像補正方法を実施する本発明の画像補正装置を利用する、放射線画像撮影装置の一例を概念的に示す。
図１に示す放射線画像撮影装置１０（以下、撮影装置１０とする）は、被写体Ｈ（被検者）の放射線画像（診断画像）を撮影する放射線画像診断装置であって、放射線画像を撮影する撮影部１２と、撮影部１２が撮影した被写体画像を処理する画像処理部１４と、モニタ１６と、プリンタ１８とを有して構成させる。なお、撮影装置１０は、公知の通信ネットワークによって、管理システムや画像サーバなどの外部の装置やシステム等に接続されてもよい。
【００２１】
撮影部１２は、被写体Ｈの放射線画像（以下、被写体画像とする）を撮影する部位であり、曝射制御手段２０と、放射線源２２と、撮影台２４と、撮影手段２６とを有して構成される。
曝射制御手段２０は、放射線源２２から放射線を照射するための高圧電源を有し、放射線源２２を駆動して、設定された強度の放射線を設定された時間だけ照射させる、公知の放射線源２２の制御手段である。放射線源２２も、各種の放射線画像撮影装置に設置される通常の放射線源である。
撮影台２４も、各種の放射線画像撮影装置で利用される通常の撮影台である。なお、撮影装置１０は、必要に応じて、放射線源２２の移動手段、撮影台２４の昇降手段や水平方向への移動手段、撮影台２４を傾ける傾斜手段等を有してもよい。
【００２２】
撮影手段２６は、放射線固定検出器３０（以下、ＦＰＤ３０とする）によって、被写体画像を撮影するものである。
撮影装置１０は、通常の放射線画像撮影装置１０と同様に、放射線源２２が照射し、被写体Ｈを透過した放射線をＦＰＤ３０の受光面で受光し、ＦＰＤ３０で放射線を光電変換することにより、被写体画像を撮影する。
【００２３】
本発明において、ＦＰＤ３０は、放射線画像撮影装置に利用される通常のＦＰＤ（Flat Panel Detector(フラットパネル検出器)）である。
本発明において、ＦＰＤ３０は各種のものが全て利用可能である。従って、アモルファスセレン等の光導電膜を有し、放射線の入射によって光導電膜が発した電荷（電子−正孔対（ｅ−ｈペア））を収集して電気信号として読み出す、いわゆる直接方式のＦＰＤでも、「ＣｓＩ：Ｔｌ」などの放射線の入射によって発光（蛍光）する蛍光体で形成されたシンチレータ層とフォトダイオードとを用い、放射線の入射によるシンチレータ層の発光をフォトダイオードで光電変換して、電気信号として読み出す、いわゆる間接方式のＦＰＤでもよい。
【００２４】
図示例の撮影装置１０において、ＦＰＤ３０は、ｘ−ｙ方向の二次元的に配列された個々の画素に対応して二次元的にＴＦＴ(Thin Film Transistor)等を配置して、ｘ−ｙ方向の二次元的に配列された画素配列の１方向（一例として、ｘ方向とする）に配置した電極に沿った読取方向に、光導電膜が発した各画素の電荷を流して、下流端部で、順次、収拾することにより、個々の画素に対応して光導電膜が発した電荷を収拾する構成を有するものである。
【００２５】
なお、撮影手段２６は、ＦＰＤ３０以外にも、必要に応じて、散乱放射線がＦＰＤ３０に入射する散乱放射線を遮蔽するためのグリッド、グリッドの移動手段等、公知の放射線画像撮影装置が有する各種の部材を有してもよいのは、もちろんである。
【００２６】
撮影手段２６（ＦＰＤ３０）が撮影した放射線画像の出力信号は、画像処理部１４に出力される。
画像処理部１４は、ＦＰＤ３０が出力した出力信号を処理して、モニタ１６による表示や、プリンタ１８でのプリント出力、さらには、ネットワークや記録媒体を用いた放射線画像（データ）の出力に対応する被写体画像（画像データ（画像信号））とするものである。図示例の撮影装置１０において、画像処理部１４は、データ処理手段３２と、画像処理手段３４とを有する。
このような画像処理部１４は、一例として、１台もしくは複数台のコンピュータやワークステーションで構成されるものであり、図示した部位以外にも、必要に応じて、各種の操作や指示の入力等をするためのキーボードやマウス等を有している。
【００２７】
データ処理手段３２は、ＦＰＤ３０の出力信号に、Ａ／Ｄ変換等の所定の処理を施して、被写体画像の画像データ（画像信号）に変換し、画像処理手段３４に供給するものである。
【００２８】
画像処理手段３４は、データ処理手段３２が処理した画像データ（被写体画像）に、所定の画像処理を施して、モニタ１６などによる画像表示、プリンタ１８などによるプリント（ハードコピー）の出力、ネットワークや記憶媒体への出力、等に対応する画像データとして、モニタ１６、プリンタ１８、およびネットワーク等の指示された部位に出力するものである。
【００２９】
画像処理手段３４が行なう画像処理には、特に限定は無く、ゲイン補正（シェーディング補正）、オフセット補正、欠陥画素補正、残像補正、階調補正、濃度補正、モニタ表示用やプリント出力用のデータに画像データを変換するデータ変換など、各種の放射線画像撮影装置や画像処理装置で行なわれている画像処理が、全て実施可能である。また、これらの補正は、全て、公知の方法で行なえばよい。
さらに、画像処理手段３４は、本発明の画像補正方法を実施する、本発明の画像補正装置にかかる、スジムラ補正手段４２を有する。
【００３０】
スジムラ補正手段４２は、画像中に存在する所定方向のスジムラ（直線状のムラ）を除去する、スジムラの除去を行なう部位であり、一例として、図２に概念図で示すように、メディアン処理手段５０、第１減算手段５２、閾値処理手段５４、ＬＰＦ処理手段５６、および、第２減算手段５８を有する。
以下、図２および図３を参照して、スジムラ補正手段４２によるスジムラ除去処理の作用を説明することにより、スジムラ補正手段４２、ならびに、本発明の画像補正方法および画像補正装置について、より詳細に説明する。
【００３１】
なお、本例において、元画像は、一例として、図３（Ａ）に概念的に示すように、ＦＰＤ３０の画素配列方向（ｘ−ｙ方向）に延在する直線で形成された格子状のパターンを有するファントムを撮影したものであり、黒色の太線がファントムの画像であり、ｘ方向に延在する灰色の細線がスジムラである。
【００３２】
スジムラ補正手段４２においては、まず、元画像（原画像）に対して、メディアン処理手段５０において、メディアン処理を施す。
なお、図示例の撮影装置１０（画像処理手段３４）において、スジムラ補正手段４２によるスジムラ除去を行なうタイミングには、特に限定はなく、画像処理手段３４で施す画像処理の種類や順番等に応じて、適宜、決定すればよい。
【００３３】
メディアン処理手段５０は、元画像に対して、元画像が有するスジムラと直交する方向に、メディアン処理を施す。
前述のように、図示例の撮影装置１０のＦＰＤ３０は、ｘ−ｙ方向に二次元的に配列された画素配列において、光導電膜が発した電荷をｘ方向（読取方向）に流して、ｘ方向の下流端部で、順次、収拾することにより、個々の画素に対応して光導電膜が発した電荷を収拾する構成を有するものである。
このような構成を有するＦＰＤ３０は、例えば大線量の放射線で撮影を行なう場合など、撮影条件等によっては、各画素から電荷が溢れてしまう場合が有り、この溢れた電荷が、読取方向であるｘ方向に流れて検出されてしまい、これにより、このｘ方向に延在してスジムラが生じてしまう。すなわち、ＦＰＤ３０では、元画像のスジムラは、このｘ方向に直線状に延在して発生することが、構造上、予め分かっている。
【００３４】
従って、メディアン処理手段５０は、電荷の読取方向であるｘ方向と直交するｙ方向にメディアン処理を行なう。
具体的には、図３（Ａ）に概念的に示すように、ｙ方向の一次元マスクを有するフィルタＦを用いて、元画像に対して、このフィルタＦをｙ方向に移動して、メディアン処理を行なう。
【００３５】
ここで、本発明において、メディアン処理を行なうフィルタＦの画素数には、特に限定はなく、想定されるスジムラのサイズ（太さ＝延在方向と直交する方向すなわちｙ方向の画素数）に応じて、適宜、設定すればよいが、想定されるスジムラのサイズの２倍以上の画素数のマスクを有するフィルタＦが好ましい。このようなフィルタＦを用いることにより、より確実に、スジムラを検出することが可能となる。
好ましくは、中央値をより簡易に検出できるように、奇数の画素数のマスクを有するのが好ましく、特に、想定されるスジムラの画素数をｎとして、「２ｎ＋１」画素のマスクを有するフィルタＦを用いるのが好ましい。周知のように、メディアン処理は、中央値を検出するための画像データの並び替え等が必要なため、データ処理量が大きくなるが、このように、「２ｎ＋１」画素のマスクを有するフィルタＦを用いることにより、データ処理量を最低限として、スジムラ除去処理を迅速に行なうことが可能となる。また、図示例のように読取方向に電極を配置し、この読取方向に光導電膜が発した電荷を流すＦＰＤ３０においては、想定されるスジムラのサイズは、一般的に、数画素程度の高周波のスジムラであるので、小さなフィルタＦで十分であり、計算量が掛からずに、スジムラを抽出できる。
【００３６】
メディアン処理手段５０で、元画像に対してｙ方向にメディアン処理を行なってメディアン処理画像を生成したら、次いで、第１減算手段５２において、元画像からメディアン処理画像を減算する。これにより、ｘ方向に延在するスジムラを含め、元画像中に存在するｘ方向に延在するスジ状（直線状）の画像を抽出した、図３（Ｂ）に示すような減算画像を生成できる。
【００３７】
前述の特許文献１に示されるように、従来のスジムラの除去処理では、低周波数の画素を通す一次元のローパスフィルタ（あるいはハイパスフィルタ）を用いて、スジムラと直交する方向およびスジムラの延在方向にフィルタによる周波数処理を行なって、スジムラ
の除去（スジムラの抽出）を行なっている。
しかしながら、例えば、図４（Ａ）の領域Ａのような濃度変化の大きいスジムラのように、急激に濃度が変動する領域では、図４（Ｂ）の上段に模式的に示すように、スジムラを含んでスジムラと直交するｙ方向の領域Ａのプロファイルは実線、これをローパスフィルタ（以下、ＬＰＦとする）でｙ方向に処理した画像は同破線の様になる。その結果、ＬＰＦで処理した画像を元画像から減算すると、図４（Ｂ）の下段に示すように、スジムラの両側にリンギングが発生し、これに起因して、スジムラを除去した画像にアーチファクト等が生じてしまう。
【００３８】
これに対し、ＬＰＦを用いた処理ではなく、スジムラと直交する方向にメディアン処理を行なう本発明によれば、同様のプロファイルを有する領域をメディアン処理して得られるメディアン処理画像は、図４（Ｃ）の上段に示すように、破線で示すような直線状となる。従って、その結果、元画像からメディアン処理画像を減算しても、図４（Ｃ）の下段に示すように、リンギングが発生する抑制することができ、これに起因して、スジムラを除去した画像にアーチファクト等が生じるのを好適に抑制できる。
【００３９】
第１減算手段５２で得られた減算画像は、次いで、閾値処理手段５４で処理される。
閾値処理手段５４は、減算画像から、濃度が所定の閾値を超える画素（画像データ）を抽出し、この画素の濃度を小さくする、閾値処理を行なう。好ましくは、閾値処理では、濃度が所定の閾値を超える画素の濃度を「０（最低濃度）」にする。これにより、図３（Ｃ）に示すような、閾値処理画像が得られる。
すなわち、閾値処理手段５４で行なう閾値処理は、減算画像において、濃度が閾値を超える領域を構造物として抽出し、この構造物の濃度を低減、好ましくは「０」とする処理である。
【００４０】
なお、濃度が閾値以上の画素の濃度を低減する方法には、特に限定はなく、他に、スジムラの延在方向において、閾値を超える前の画素の濃度値（画像データ）をコピーすることによって、閾値を超えた画素の濃度を小さくしてもよい。
【００４１】
ここで、図３（Ｂ）および図５（Ａ）の左側に示すような減算画像のｘ方向の領域Ｂが、図５（Ａ）の右側に示すようなプロファイルを有している場合、濃度が閾値を超えた画素の濃度を０にするような閾値処理を行なうと、閾値処理済みの領域Ｂのプロファイルは、図５（Ｂ）の右側に示すようになり、これにより、図３（Ｃ）および図５（Ｂ）の左側に示すような、減算画像から、高濃度領域を除去した、閾値処理画像が得られる。
【００４２】
メディアン処理画像を元画像から減算した減算画像には、スジムラと同方向に延在する直線状の構造物が残っている可能性が有る。
後述するが、本発明においては、この後、閾値処理画像にローパスフィルタによる処理を行なって、ｘ方向の長さが短い構造物等を除去してスジムラを抽出し、このスジムラを元画像から減算することにより、元画像からスジムラを除去する。ここで、閾値処理をしないで、このローパスフィルタによる処理を行なうと、スジムラと同方向に延在する直線状の構造物もスジムラとして抽出してしまうので、元画像から、必要な構造物も除去してしまう。
これに対し、本発明では、濃度が所定の閾値を超えている領域は、構造物として濃度を低減、好ましくは「０」にする。これにより、元画像中の構造物をスジムラとして抽出することを抑制し、より高精度なスジムラの除去を行なうことが可能になる。
【００４３】
ＦＰＤ３０で撮影した一般的な放射線画像（各種の固体撮像素子による画像）においては、前述の画素（光導電膜）の飽和等に起因するスジムラの濃度は、曝射量等の撮影条件やＦＰＤ３０の温度等によって異なるものの、元画像の構造物よりも低濃度であるのが通常である。
すなわち、構造物とスジムラとの間には、ある程度の濃度差が有り、この間に閾値を設定することにより、閾値処理画像において、スジムラと同方向に延在する直線状の構造物を、取り除くことができる。
【００４４】
本発明において、閾値の設定方法には、特に限定はなく、一例として、実験やシミュレーションを行なって、適宜、設定すればよい。例えば、ファントム等を用い、撮影条件を、種々、変更して実験やシミュレーションを行なって、構造物を取り除くことが無く、かつ、スジムラを取り除くことができる濃度を、適宜、選択／設定すればよい。
【００４５】
なお、本発明においては、閾値処理における閾値は、撮影条件（管電圧、曝射時間、撮影線量等）、撮影部位、ＦＰＤ３０の温度、前回の撮影からの経過時間などに応じて、適宜、変更してもよい。
また、撮影装置１０においては、上記撮影条件や撮影部位等に応じた閾値を、複数、準備しておき、最適な閾値を選択して、閾値処理を行なうにしてもよい。
【００４６】
次いで、閾値処理手段５４によって得られた閾値処理画像は、次いで、ＬＰＦ処理手段５６において処理される。
ＬＰＦ処理手段５６は、スジムラと直交する方向すなわちｙ方向に延在する一次元のＬＰＦ（低周波数の空間周波数成分のみを抽出するＬＰＦ）によって、スジムラの延在方向すなわちｘ方向に処理を行ない、ｘ方向に延在する領域をスジムラとして抽出する。
これにより、図３（Ｃ）に示すような閾値処理画像から、ｘ方向の長さが短い構造物等を除去して、図３（Ｄ）に示すように、元画像からスジムラのみを抽出した処理用画像を生成できる。
ＬＰＦ処理手段５６で用いるＬＰＦは、カットオフ周波数が高すぎると、スジムラとして抽出したくない短い構造物を抽出する可能性があり、逆に、カットオフ周波数が低すぎると、抽出したいスジムラをスジムラとして抽出できなくなる。ＬＰＦの特性は、例えば、ファントム等を用い、撮影条件を、種々、変更して実験やシミュレーションを行なって、スジムラを取り除くことができ、かつ、構造物を抽出することの無い特性を、適宜、選択／設定すればよい。
【００４７】
このようにして、処理用画像を生成したら、最後に、第２減算手段５８において、元画像から処理用画像を減算して、元画像から処理用画像を減算して、図３（Ｅ）に示すように、元画像からスジムラを除去した、スジムラ補正済画像を得る。
【００４８】
以上の説明より明らかなように、本発明によれば、スジムラと直交する方向に、ＬＰＦによる処理ではなく、メディアン処理を行なうので、濃度変化の大きいスジムラや構造物のエッジ等の濃度が急激に変動する部分において、リンギングが発生することを抑制することができ、リンギングに起因するアーチファクト等を生じる事を抑制して、元画像からスジムラの除去を行なうことができる。
また、処理用画像となる、元画像からメディアン処理画像を減算して得られた減算画像に、濃度が閾値を超える領域（画素）の濃度を低減、好ましくが「０」にする閾値処理を行なうことにより、元画像中でスジムラと同方向に延在する構造物をスジムラとして抽出することが抑制でき、スジムラの除去処理によって、必要な構造物を除去してしまうことも、好適に防止できる。
【００４９】
このようにして、画像処理手段３４において、スジムラ除去を含む所定の画像処理を施された放射線画像は、モニタ１６に出力されて可視像として表示され、また、プリンタ１８に出力されてハードコピーとして出力され、また、ネットワークなどを利用して指示された場所に出力される。
【００５０】
以上の例は、本発明を、ＦＰＤ３０で撮影した放射線画像からのスジムラ除去に利用した例であるが、本発明は、これに限定はされず、延在方向が既知であるスジムラが発生する、各種の画像に応じて、スジムラの除去を行なうことが可能である。
例えば、一方向に読取画素が配列されたラインＣＣＤセンサ等のラインセンサと、ラインセンサと同方向に延在する線状光源を用い、ラインセンサの読取画素の配列方向と直交する方向に、ラインセンサおよび線状光源と原稿とを走査（相対的に移動）することで原稿の画像を読み取るスキャナでは、ラインセンサに不良画素が発生すると、この不良画素に起因して、走査方向（ラインセンサと原稿との相対移動方向）に延在するスジムラが生じる。また、線状光源の一部の光量変動が生じても、走査方向に延在するスジムラが生じる場合がある。
従って、このようなスキャナが読み取った画像を元画像とする場合には、前記走査方向と直交する方向にメディアン処理を行なえばよい。
【００５１】
また、同様のラインセンサと、同方向に延在するＬＥＤアレイ等のライン状の励起光光源とを組にした読取手段を用い、ラインセンサの延在方向と直交する方向に、ライン状の励起光光源から励起光を照射しつつ、放射線画像変換パネル（ＩＰ）と読取手段とを走査することにより、ラインセンサによってＩＰに撮影された放射線画像を読み取る、放射線画像読取装置（前記ＦＣＲ等）が知られている。
この放射線画像読取装置においては、同様に、ラインセンサに不良画素が生じた場合や、ライン状の励起光光源の或る点に輝度不足等の不良が生じた場合に、走査方向に延在するスジムラが発生する。
従って、このような放射線画像読取装置による放射線画像を元画像する場合にも、同様に、前記走査方向と直交する方向にメディアン処理を行なえばよい。
【００５２】
また、特許文献１にも示されるように、静止グリッドを用いる放射線画像の撮影においては、グリッドの延在方向に延在するスジムラが、スジムラと直交する方向に周期的に発生するので、メディアン処理手段５０において、グリッドの延在方向と直交する方向にメディアン処理を行なえばよい。
さらに、デジタルカメラ等に利用される、２次元的に画素が配列されたＣＣＤセンサでは、高光量の光が照射されて画素が飽和すると（電荷蓄積が限界を超えると）、いわゆる垂直方向（電荷の転送方向）に白抜きのようなスジムラが生じる、いわゆるスミアが発生することが知られている。従って、このようなＣＣＤセンサが撮影した画像を元画像とする場合には、ＣＣＤにおける垂直方向（電荷の転送方向）と直交する方向にメディアン処理を行なえばよい。
【００５３】
以上、本発明の画像補正方法および画像補正装置について詳細に説明したが、本発明は、上記実施例に限定はされず、本発明の要旨を逸脱しない範囲において、各種の改良や変更を行なってもよいのは、もちろんである。
【産業上の利用可能性】
【００５４】
本発明は、放射線画像撮影装置、ＣＣＤセンサ（ライン／エリア）を用いるスキャナや撮影装置等によって撮影された画像の画像処理や、これらの装置に、好適に利用可能である。
【符号の説明】
【００５５】
１０  （放射線画像）撮影装置
１２  撮影部
１４  画像処理部
１６  モニタ
１８  プリンタ
２０  曝射制御手段
２２  放射線源
２４  撮影台
２６  撮影手段
３０  ＦＰＤ（放射線固体検出器）
３２  データ処理手段
３４  画像処理手段
４２  スジムラ補正手段
５０  メディアン処理手段
５２  第１減算手段
５４  閾値処理手段
５６  ＬＰＦ処理手段
５８  第２減算手段
【図１】
2010239480.tif 000003
【図２】
2010239480.tif 000004
【図３】
2010239480.tif 000005
【図４】
2010239480.tif 000006
【図５】
2010239480.tif 000007