(57)【特許請求の範囲】
【請求項１】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理方法において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１のステップと、
該第１のステップによって得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２のステップと、
前記第１のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第３のステップを含み、該第３のステップによって得られた前記加算された空間周波数成分を前記原画像信号から減算する減算ステップと、
を含むことを特徴とする周期的パターン抑制処理方法。
【請求項２】
前記第２のステップのフィルタ処理における遮断周波数より高い遮断周波数を用いて、前記第１のステップによって得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第５のステップを更に含み、
前記第３のステップは、前記第５のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第５のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算するステップであることを特徴とする請求項１に記載の周期的パターン抑制処理方法。
【請求項３】
前記第１のステップによって得られた画像信号と前記第５のステップによって得られた画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくする第６のステップを更に含み、
前記第３のステップは、前記第５のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第５のステップによって得られた画像信号の空間周波数成分の重みを大きくした空間周波数成分に、更に前記第６のステップによって得られた空間周波数成分を加算して、該加算した空間周波数成分を前記第２のステップによって得られた画像信号の空間周波数成分に加算するものであることを特徴とする請求項２に記載の周期的パターン抑制処理方法。
【請求項４】
前記原画像は、静止グリッドを介して撮影された放射線画像と、前記静止グリッドに対応する縞模様状のグリッド像とからなるものであり、前記周期的パターンは前記グリッド像であることを特徴とする請求項１〜３の何れか一項に記載の周期的パターン抑制処理方法。
【請求項５】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理装置において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１の手段と、
該第１の手段において得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２の手段と、
前記第１の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第３の手段を備え、該第３の手段において得られた前記加算された空間周波数成分を前記原画像信号から減算する減算手段と、
を備えることを特徴とする周期的パターン抑制処理装置。
【請求項６】
前記第２の手段のフィルタ処理における遮断周波数より高い遮断周波数を用いて、前記第１の手段において得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第５の手段を更に備え、
前記第３の手段は、前記第５の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第５の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算するものであることを特徴とする請求項５に記載の周期的パターン制処理装置。
【請求項７】
前記第１の手段において得られた画像信号と前記第５の手段において得られた画像信号のコントラスト差が大きいほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくする第６の手段を更に含み、
前記第３の手段は、前記第５の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第５の手段において得られた画像信号の空間周波数成分の重みを大きくした空間周波数成分に、更に前記第６の手段において得られた空間周波数成分を加算して、該加算した空間周波数成分を前記第２の手段において得られた画像信号の空間周波数成分に加算するものであることを特徴とする請求項６に記載の周期的パターン抑制処理装置。
【請求項８】
前記原画像は、静止グリッドを介して撮影された放射線画像と、前記静止グリッドに対応する縞模様状のグリッド像とからなるものであり、前記周期的パターンは前記グリッド像であることを特徴とする請求項５〜７の何れか一項に記載の周期的パターン抑制処理装置。
【請求項９】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理方法において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１のステップと、
該第１のステップによって得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２のステップと、
前記原画像信号の濃度が高いほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第４のステップを含み、該第４のステップによって得られた前記加算された空間周波数成分を前記原画像信号から減算する減算ステップと、
を含むことを特徴とする周期的パターン抑制処理方法。
【請求項１０】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理方法において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１のステップと、
該第１のステップによって得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２のステップと、
前記第１のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第３のステップ、及び、前記原画像信号の濃度が高いほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第４のステップを含み、前記第３のステップ及び前記第４のステップによって得られた前記加算された空間周波数成分を前記原画像信号から減算する減算ステップと、
を含むことを特徴とする周期的パターン抑制処理方法。
【請求項１１】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理装置において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１の手段と、
該第１の手段において得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２の手段と、
前記原画像信号の濃度が高いほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第４の手段を備え、該第４の手段において得られた前記加算された空間周波数成分を前記原画像信号から減算する減算手段と、
を備えることを特徴とする周期的パターン抑制処理装置。
【請求項１２】
縞模様状の周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理装置において、
前記原画像信号から前記周期的パターンの並び方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第１の手段と、
該第１の手段において得られた画像信号から前記周期的パターンと平行な方向における該周期的パターンが持つ周波数を含む周波数成分を抽出する一次元フィルタ処理を行う第２の手段と、
前記第１の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第３の手段、及び、前記原画像信号の濃度が高いほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第４の手段を備え、前記第３の手段及び前記第４の手段において得られた前記加算された空間周波数成分を前記原画像信号から減算する減算手段と、
を備えることを特徴とする周期的パターン抑制処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像信号における周期的パターンに応じた空間周波数成分を抑制する周期的パターン抑制処理方法及び装置に関するものである。
【背景技術】
【０００２】
従来より、放射線（Ｘ線、α線、β線、γ線、電子線、紫外線等）が照射されることにより放射線エネルギーの一部を蓄積し、蓄積領域に可視光やレーザ光等の励起光が照射されると、蓄積された放射線エネルギーに応じて輝尽発光光を発光する蓄積性蛍光体（輝尽性蛍光体）が放射線画像読取装置等に広く利用されている。具体的には、例えば、支持体上に蓄積性蛍光体を積層した蓄積性蛍光体シートに人体などの被写体を透過した放射線が照射されるように蓄積性蛍光体シートを配置する。蓄積性蛍光体シートは被写体を透過した放射線によって形成された放射線画像情報を一旦蓄積記録する。そして、蓄積性蛍光体シートに励起光を照射して、輝尽発光光を生じさせる。この輝尽発光光を光電変換すると放射線画像情報に応じた画像信号を得ることができる。
【０００３】
上述した蓄積性蛍光体シートに被写体の放射線画像を撮影記録する際、被写体により散乱された放射線が蓄積性蛍光体シートに照射されないように、放射線の透過しない鉛等と、透過しやすいアルミニウムや木材等とが約４本／ｍｍのピッチで交互に配置された静止グリッドを被写体と蓄積性蛍光体シートの間に配置して撮影を行うことがある。この静止グリッドを用いて撮影を行うと、被写体によって散乱された放射線が蓄積性蛍光体シートに照射されにくくなるため、被写体の放射線画像のコントラストを向上させることができる。しかし、被写体による放射線画像と共に、静止グリッドに対応した縞模様状の周期的パターン（以下「グリッド像」と表記する）も記録されることになる。
【０００４】
そこで、グリッド像が混入した原画像からグリッド像を除去する方法が特許文献１に開示されている。具体的には、グリッド像の縞模様の並び方向に一次元ハイパスフィルタ処理を施し、更にグリッド像の縞模様と平行な方向に一次元ローパスフィルタ処理を施してグリッド像に応じた空間周波数成分を原画像から抽出し、抽出された空間周波数成分を原画像から差し引く。これにより、原画像に混入されたグリッド像を除去することができる。
【特許文献１】特開２００３−１５０９５４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に開示された方法では、以下のような問題点があった。原画像に対して上記した一次元ハイパスフィルタ処理と一次元ローパスフィルタ処理を施すと、元々グリッド像が存在しない部分にもグリッド像が現れてしまう。元々グリッド像が存在しない部分とは、主に原画像においてグリッド像が潰れている部分（つまり、原画像の高濃度部分）に相当する。従って、原画像に対してグリッド像を除去する処理を行った結果、高濃度部分に元々存在しない縞模様が混入してしまい問題となっていた。
【０００６】
一次元ローパスフィルタの遮断周波数を低く設定すれば、縞模様の発生を抑えることができるが、高周波数成分も遮断することになるため、処理後の画像のコントラストが低下してしまう。
【０００７】
本発明は、上記事情を鑑みてなされたものであり、画質を劣化させることなく原画像の周期的パターンを抑制する周期的パターン抑制処理方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
以上の課題を解決するために、本発明の周期的パターン抑制処理方法は、周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理方法において、前記原画像信号に対して前記周期的パターンの並び方向にフィルタ処理を行う第１のステップと、該第１のステップによって得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第２のステップと、前記第１のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第３のステップ、及び、前記原画像信号の濃度が高いほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算する第４のステップ、の少なくとも何れか一方を含み、前記第３のステップ及び／又は前記第４のステップによって得られた前記加算された空間周波数成分を前記原画像信号から減算する減算ステップと、を含むことを特徴とする。
【０００９】
ここで「周期的パターン」とは、例えば、静止グリッドを介して放射線撮影された際に原画像に含まれる静止グリッドに対応した縞模様状のグリッド像、又はデジタルカメラ等を用いた通常撮影において網戸やブラインド等を介して撮影された際に撮影画像に含まれる周期的パターン等を言う。更に、例えば、静止グリッドそのものの空間周波数成分によるグリッド像だけでなく、放射線画像を読み取る際に用いられるサンプリング周波数の数値や縮小処理等に起因して発生した縞状モアレも含む。
【００１０】
また、「周期的パターンの並び方向」とは、例えば上記グリッド像の縞模様の並び方向を示し、「周期的パターンと平行な方向」とは、例えば上記グリッド像の縞模様と平行な方向、即ち、上記グリッド像の縞模様の並び方向と直交する方向のことを示す。
【００１１】
また、「画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして」とは、単にコントラスト差に比例して第１のステップによって得られた画像信号の空間周波数成分の重みの度合いを大きくしていくこととしてもよいし、コントラスト差が所定の範囲内である場合は空間周波数成分の重みの度合いを一定とし、コントラスト差が所定の範囲を超えた場合は空間周波数成分の重みの度合いをコントラスト差に比例して大きくしていくこととしてもよい。
【００１２】
また、前記第２のステップのフィルタ処理における遮断周波数より高い遮断周波数を用いて、前記第１のステップによって得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第５のステップを更に含み、前記第３のステップは、前記第５のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第５のステップによって得られた画像信号の空間周波数成分の重みを大きくして前記第２のステップによって得られた画像信号の空間周波数成分に加算するステップであることとしてもよい。
【００１３】
また、前記第１のステップによって得られた画像信号と前記第５のステップによって得られた画像信号のコントラスト差が大きいほど前記第１のステップによって得られた画像信号の空間周波数成分の重みを大きくする第６のステップを更に含み、前記第３のステップは、前記第５のステップによって得られた画像信号と前記第２のステップによって得られた画像信号のコントラスト差が大きいほど前記第５のステップによって得られた画像信号の空間周波数成分の重みを大きくした空間周波数成分に、更に前記第６のステップによって得られた空間周波数成分を加算して、該加算した空間周波数成分を前記第２のステップによって得られた画像信号の空間周波数成分に加算するものであることとしてもよい。
【００１４】
更に、本発明における周期的パターン抑制処理装置は、周期的パターンを含む原画像の原画像信号から、フィルタリング処理により前記周期的パターンの空間周波数成分を抽出し、該空間周波数成分を前記原画像信号から減算することにより前記原画像から前記周期的パターンを抑制する周期的パターン抑制処理装置において、前記原画像信号に対して前記周期的パターンの並び方向にフィルタ処理を行う第１の手段と、該第１の手段において得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第２の手段と、前記第１の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第３の手段、及び、前記原画像信号の濃度が高いほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算する第４の手段、の少なくとも何れか一方を備え、前記第３の手段及び／又は前記第４の手段において得られた前記加算された空間周波数成分を前記原画像信号から減算する減算手段と、を備えることを特徴とする。
【００１５】
また、前記第２の手段のフィルタ処理における遮断周波数より高い遮断周波数を用いて、前記第１の手段において得られた画像信号に対して前記周期的パターンと平行な方向にフィルタ処理を行う第５の手段を更に備え、前記第３の手段は、前記第５の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第５の手段において得られた画像信号の空間周波数成分の重みを大きくして前記第２の手段において得られた画像信号の空間周波数成分に加算するものであることとしてもよい。
【００１６】
また、前記第１の手段において得られた画像信号と前記第５の手段において得られた画像信号のコントラスト差が大きいほど前記第１の手段において得られた画像信号の空間周波数成分の重みを大きくする第６の手段を更に含み、前記第３の手段は、前記第５の手段において得られた画像信号と前記第２の手段において得られた画像信号のコントラスト差が大きいほど前記第５の手段において得られた画像信号の空間周波数成分の重みを大きくした空間周波数成分に、更に前記第６の手段において得られた空間周波数成分を加算して、該加算した空間周波数成分を前記第２の手段において得られた画像信号の空間周波数成分に加算するものであることとしてもよい。
【発明の効果】
【００１７】
第１のステップによって得られた画像信号と第２のステップによって得られた画像信号のコントラスト差が大きいほど第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして第２のステップによって得られた画像信号の空間周波数成分に加算する、及び／又は、原画像信号の濃度が高いほど第１のステップによって得られた画像信号の空間周波数成分の重みを大きくして第２のステップによって得られた画像信号の空間周波数成分に加算し、第３のステップ及び／又は第４のステップによって得られた加算された空間周波数成分を原画像信号から減算して画像信号を得ることにより、原画像の画質を劣化させることなく原画像に混入されている周期的パターンを抑制し、更に本来グリッド像が存在していない部分（例えば、原画像における高濃度部分）における周期的パターンの出現をも抑えることができる。
【発明を実施するための最良の形態】
【００１８】
以下、図面を参照して本発明の実施の形態を説明する。尚、以下の実施の形態では、本発明における周期的パターン抑制処理装置を放射線画像読取装置に用いた場合について説明するが、デジタルカメラ等を用いた通常撮影において網戸やブラインド等を介して撮影した際に撮影画像に含まれる周期的パターンを抑制するための画像処理装置等に用いられてもよい。尚、放射線画像読取装置は、蓄積性蛍光体シートに記録された人体の放射線画像をレーザビーム走査によりデジタル画像信号として読み取るものである。
【００１９】
図１は放射線画像撮影装置の概略図である。放射線源１から放射された放射線２は、被写体３を透過して静止グリッド（以下、単に「グリッド」と表記する）４に到達する。グリッド４は放射線２を吸収する鉛４ａと、放射線２を透過するアルミニウム４ｂとが例えば４本／ｍｍ程度のピッチで交互に配置されているものである。また、放射線２がアルミニウム４ｂを透過して蓄積性蛍光体シート１１に入射するように鉛４ａは位置に応じて多少傾きを変化させて設置されている。
【００２０】
従って、被写体３を透過した放射線２は、鉛４ａに吸収されて蓄積性蛍光体シート１１への到達を遮られる一方で、アルミニウム４ｂを透過して蓄積性蛍光体シート１１に到達する。蓄積性蛍光体シート１１には被写体像と共に４本／ｍｍの縞模様状のグリッド像が記録される。一方、被写体３内で散乱された散乱放射線２ａは、位置に応じて傾きを持って設置された鉛４ａによって吸収され、又はグリッド４の表面で反射するため、蓄積性蛍光体シート１１には到達しない。従って、蓄積性蛍光体シート１１は散乱放射線２ａの照射の少ない鮮明な放射線画像を記録することができる。図２は、図１に示した放射線画像撮影装置を用いて撮影した際に蓄積性蛍光体シート１１に蓄積記録された被写体像５とグリッド像６の放射線画像の一例である。
【００２１】
図３は、放射線画像読取装置の斜視図と機能ブロック図を組み合わせた図である。読取部１０の所定位置にセットされた蓄積性蛍光体シート１１は、駆動手段（不図示）によって駆動されるエンドレスベルト等のシート搬送手段１５によって、例えば走査ピッチ１０本／ｍｍで矢印Ｙ方向に搬送（副走査）される。一方、レーザ光源１６から発せられた光ビーム１７は、モータ２４により駆動され矢印方向に高速回転する回転多面鏡１８によってｆθレンズ等の集束レンズ１９に向けて反射される。光ビーム１７は集束レンズ１９を通過した後、ミラー２０によって蓄積性蛍光体シート１１に向けて反射される。蓄積性蛍光体シート１１は光ビーム１７によって副走査方向（矢印Ｙ方向）と略直角な矢印Ｘ方向に主走査される。
【００２２】
蓄積性蛍光体シート１１において光ビーム１７が照射された箇所からは蓄積記録されている放射線画像情報に応じた光量の輝尽発光光２１が発散される。輝尽発光光２１は光ガイド２２の入射端面２２ａから入射し、光ガイド２２内の内部を全反射を繰り返して出射端面２２ｂから出射する。出射された輝尽発光光２２はフォトマルチプライヤ２３によって受光されて光電変換され、アナログ画像信号Ｓａに変換される。
【００２３】
アナログ画像信号Ｓａはログアンプ２６によって対数的に増幅された後、Ａ／Ｄ変換器２８において空間周波数ｆｓ＝１０ｃｙｃｌｅ／ｍｍに対応するサンプリング間隔でサンプリングされてデジタル化され、デジタル画像信号Ｓｄ（以下、単に「画像信号Ｓｄ」と表記する）を出力する。この画像信号Ｓｄは、図４に示すように、蓄積性蛍光体シート１１に対して主走査方向（横方向）に光ビーム１７を走査させながら蓄積性蛍光体シート１１を副走査方向（縦方向）に移動させることにより蓄積性蛍光体シート１１を２次元走査して得られた放射線画像情報を示す。尚、このようにして得られた画像信号Ｓｄには、被写体３に対応する放射線画像５の情報の他に、グリッド４に対応するグリッド像６の情報も含まれている。
【００２４】
画像信号Ｓｄは一旦記憶部２９に記憶された後、画像信号処理部３０に入力される。画像信号処理部３０は本発明における周期的パターン抑制処理方法を実行するための周期的パターン抑制処理部（周期的パターン抑制処理装置に相当）４０を備えている。以下、実施例を用いて、この周期的パターン抑制処理部４０について詳しく説明する。
【００２５】
（実施例１）
図５（ａ）は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ａの構成を説明する為のフローチャートである。記憶部２９から読み出された画像信号Ｓｄは、横方向ハイパスフィルタ（以下、略して「横方向ＨＰＦ」と表記する）４１によってＨＰＦ処理が施される（第１のステップ）。ＨＰＦ処理後の信号を信号Ｓｈとする。ここで「横方向ＨＰＦ」は、図４の“主走査方向（横方向）”に示す矢印の方向、つまりグリッド像の縞模様の方向にＨＰＦ処理を施す手段のことを言う。そして信号Ｓｈは更に縦方向ローパスフィルタ（以下、略して「縦方向ＬＰＦ」と表記する）４２によってＬＰＦ処理が施される（第２のステップ）。ＬＰＦ処理後の信号を信号Ｓｈｌとする。ここで「縦方向ＬＰＦ」は、図４の“副走査方向（縦方向）”に示す矢印の方向、つまりグリッド像の縞模様と平行な方向にＬＰＦ処理を施す手段のことを言う。そして減算器４３は信号Ｓｈと信号Ｓｈｌを入力して減算処理を行って信号Ｓ１を出力する。ここで、信号Ｓ１＝（信号Ｓｈ−信号Ｓｈｌ）である。この信号Ｓ１は信号Ｓｈと信号Ｓｈｌのコントラスト差（濃淡差）を示す信号であり、実施例１ではこのコントラスト差に応じて画像信号Ｓｄから減算する空間周波数成分を変化させることによって周期的パターン（グリッド像）を抑制する方法を適用した周期的パターン処理部について説明している。
【００２６】
次に、信号Ｓ１は第１関数処理部４４によって関数処理が施される。ここで、第１関数処理部４４において実行される関数のグラフを図５（ｂ）に示す。同図に示すグラフは、横軸が入力（即ち、信号Ｓ１＝Ｓｈ−Ｓｌ）であり、縦軸が出力（即ち、ｆｃ（Ｓｈ−Ｓｌ））となっている。ここで、実線で示した関数ｆｃを用いて詳しく説明する。関数ｆｃは、−ＴＨ≦信号Ｓ１≦ＴＨの場合、出力＝０とし、信号Ｓ１＞ＴＨ、及び信号Ｓ１＜−ＴＨの場合、出力は信号Ｓ１に比例して増減するグラフとなっている。即ち、ＴＨ≦信号Ｓ１≦ＴＨの範囲は出力＝０であるが、信号Ｓ１がＴＨを超えている場合、信号Ｓ１の値が大きいほど、つまり信号Ｓｈと信号Ｓｈｌのコントラスト差が大きいほど、出力が大きくなる。尚、撮影部位や撮影条件等に応じて、ＴＨの値を変化させてもよい。即ち、撮影部位や撮影条件等に応じて、複数の関数（例えば、関数ｆｃ、関数ｆｃ１、関数ｆｃ２等）から選択可能なようにしてもよい。
【００２７】
第１関数処理部４４による関数処理後の信号を信号Ｓ２＝ｆｃ（Ｓｈ−Ｓｈｌ）とする。加算器４５は信号Ｓ２と信号Ｓｈｌを入力して加算処理を行い、信号Ｓ３＝Ｓ２＋Ｓｈｌ＝ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌを出力する（第３のステップ）。更に減算器４６は信号Ｓ３と画像信号Ｓｄを入力して減算処理を行い（減算ステップ）、画像信号Ｓｐ＝Ｓｄ−Ｓ３＝Ｓｄ−｛ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。
【００２８】
ここで、従来の周期的パターン抑制処理方法について説明する。図７（ａ）は画像信号Ｓｄに対応する画像の一例である。被写体に対応した放射線画像Ｘにグリッド像Ｇが含まれている。また放射線画像Ｘにおける領域Ｘａは高濃度部分であることを示す。図７（ａ）に示す放射線画像に対して横方向ＨＰＦ処理を施した画像が図７（ｂ）に示す画像であり、図７（ｂ）の画像に対して更に縦方向ＬＰＦを施した画像が図７（ｃ）に示す画像である。即ち、図７（ｃ）の画像は、図５のフローチャートにおける信号Ｓｈｌに対応する画像である。ここで、図７（ｃ）に示す画像には、図７（ａ）の領域Ｘａに対応する領域においても、グリッド像が存在している。そして図７（ａ）の画像から図７（ｃ）の画像を減算した画像が図７（ｄ）に示す画像である。ここで、図７（ｄ）の画像は放射線画像Ｘのみが抽出されていることが望ましいが、縦方向ＬＰＦ後の画像（即ち、図７（ｃ）に示す画像）には本来グリッド像が存在しない領域Ｘａに相当する領域にもグリッド像が存在している為、図７（ｄ）に示すように領域Ｘａにグリッド像が混入されてしまっていた。
【００２９】
そこで、本発明における実施例１では、横方向ＨＰＦ処理後の画像（例えば、図７（ｂ）に示す画像）と縦方向ＬＰＦ処理後の画像（例えば、図７（ｃ）に示す画像）のコントラスト差を算出し（図５（ａ）における減算器４３に相当）、このコントラスト差が大きい領域においては横方向ＨＰＦ処理後の画像に対応する空間周波数成分を大きくして（図５（ａ）における第１関数処理部４４に相当）、縦方向ＬＰＦ処理後の信号Ｓｈｌに加算する（図５（ａ）における加算器４５に相当）することとした。これにより、画像全体に混入されているグリッド像を抑制しつつ、本来グリッド像が存在していない高濃度部分（例えば、図７（ｄ）における領域Ｘａ）におけるグリッド像の出現をも抑えることができる。
【００３０】
図５（ａ）に示したフローチャートに対応する式は以下の通りとなる。
【００３１】
Ｓｐ＝Ｓｄ−｛ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝・・・（１）
尚、本実施例における図５（ｂ）の示す関数ｆｃは、信号Ｓｈの濃度が信号Ｓｈｌの濃度より高いほど、信号Ｓｈの空間周波数成分の重みを大きくすることとして示したが、周期的パターン抑制処理部４０ａに入力される前に画像信号Ｓｄに対して階調反転処理等が施されている場合、信号Ｓｈの空間周波数成分の重みの付け方も反転する。即ち、画像信号Ｓｄに対して階調反転処理が施されている場合は、信号Ｓｈの濃度が信号Ｓｈｌの濃度より低いほど、信号Ｓｈの空間周波数成分の重みを大きくする。
【００３２】
また、第１関数処理部４４において用いられる関数は、図６（ｂ）に示すようなグラフに対応する関数でもよい。図６（ｂ）に示す関数をｆｃ１とする。関数ｆｃ１は、−ＴＨ≦信号Ｓ１≦＋ＴＨの範囲では、信号Ｓ１の大きさに比例して出力が増減し、−ＴＨ＞信号Ｓ１、及び＋ＴＨ＜信号Ｓ１の範囲では出力は一定となっている。図６（ａ）は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ｂの構成を説明する為のフローチャートである。記憶部２９から読み出された画像信号Ｓｄは、横方向ＨＰＦ４１によってＨＰＦ処理が施される。ＨＰＦ処理後の信号を信号Ｓｈとする。そして信号Ｓｈは更に縦方向ＬＰＦ４２によってＬＰＦ処理が施される。ＬＰＦ処理後の信号を信号Ｓｈｌとする。そして減算器４３は信号Ｓｈと信号Ｓｈｌを入力して減算処理を行って信号Ｓ１を出力する。次に、信号Ｓ１は第３関数処理部５４によって関数処理が施される。第１関数処理部５４において実行される関数のグラフが図７（ｂ）に示されたグラフである。
【００３３】
第３関数処理部５４による関数処理後の信号を信号Ｓ１８＝ｆｃ１（Ｓｈ−Ｓｈｌ）とする。更に減算器５３は信号Ｓ１と信号Ｓ１８を入力して減算処理を行い、信号Ｓ１９を出力する。加算器４５は信号Ｓ１９と信号Ｓｈｌを入力して加算処理を行い、信号Ｓ３＝Ｓ１９＋Ｓｈｌ＝｛（Ｓｈ−Ｓｈｌ）−ｆｃ１（Ｓｈ−Ｓｈｌ）｝＋Ｓｈｌ＝Ｓｈ−ｆｃ１（Ｓｈ−Ｓｈｌ）を出力する。更に減算器４６は信号Ｓ３と画像信号Ｓｄを入力して減算処理を行い、画像信号Ｓｐ＝Ｓｄ−Ｓ３＝Ｓｄ−｛Ｓｈ−ｆｃ１（Ｓｈ−Ｓｈｌ）｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。この関数ｆｃ１を用いる場合、画像信号Ｓｐは
Ｓｐ＝Ｓｄ−｛Ｓｈ−ｆｃ１（Ｓｈ−Ｓｈｌ）｝・・・（２）
の式で表されることとなる。
【００３４】
（実施例２）
実施例１では、信号Ｓｈと信号Ｓｈｌに対応する画像のコントラスト差の大きい領域では信号Ｓｈの空間周波数の重みを大きくして信号Ｓｈｌと加算し、加算後の信号を画像信号Ｓｄから減算することによってグリッド像を抑制する方法を説明した。実施例２では、画像信号Ｓｄの高濃度部分において信号Ｓｈの空間周波数の重みを大きくして信号Ｓｈｌと加算し、加算後の信号を画像信号Ｓｄから減算することによってグリッド像を抑制する方法について説明する。
【００３５】
図８（ａ）は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ｃの構成を説明する為のフローチャートである。記憶部２９から読み出された画像信号Ｓｄは、横方向ＨＰＦ４１によってＨＰＦ処理が施される（信号Ｓｈ）。信号Ｓｈは更に縦方向ＬＰＦ４２によってＬＰＦ処理が施される（信号Ｓｈｌ）。そして減算器４３は信号Ｓｈと信号Ｓｈｌを入力して減算処理を行い、信号Ｓ１を出力する。
【００３６】
一方、画像信号Ｓｄは第２関数処理部４７よって関数処理が施される。ここで、第２関数処理４７において実行される関数のグラフを図８（ｂ）に示す。同図に示すグラフは、横軸が入力（即ち、画像信号Ｓｄ）であり、縦軸が出力となっている。また、横軸における縦軸に近い方が高濃度であることを示し、縦軸から遠ざかるに従って低濃度であることを示す。ここで、実線で示した関数ｆｄを用いて詳しく説明する。関数ｆｄは、０≦画像信号Ｓｄ≦ＴＨ１の場合、出力＝１となり、ＴＨ２≦画像信号Ｓｄの場合、出力＝０となる。ＴＨ１＜画像信号Ｓｄ＜ＴＨ２の範囲は、画像信号Ｓｄの濃度が高いほど出力値が大きくなり、濃度が低いほど出力値が小さくなるグラフとなっている。尚、撮影部位や撮影条件等に応じて、ＴＨ１及びＴＨ２の値を変化させてもよい。即ち、撮影部位や撮影条件等に応じて、複数の関数（例えば、関数ｆｄ、関数ｆｄ１、関数ｆｄ２等）から選択可能なようにしてもよい。
【００３７】
第２関数処理部４６による関数処理後の信号を信号Ｓ６＝ｆｄ（Ｓｄ）とする。乗算器４８は信号Ｓ６と信号Ｓ１を入力して乗算処理を行い、信号Ｓ４＝ｆｄ（Ｓｄ）×Ｓ１＝ｆｄ（Ｓｄ）×（Ｓｈ−Ｓｈｌ）を出力する。そして信号Ｓ４と信号Ｓｈｌは加算器４５によって加算処理が施され（第４のステップ）、信号Ｓ５＝Ｓ４＋Ｓｈｌ＝ｆｄ（Ｓｄ）×（Ｓｈ−Ｓｈｌ）＋Ｓｈｌが出力される。ここで、第２関数処理部４６では、画像信号Ｓｄが高濃度を示す値であるほど１に近い値を出力し、画像信号Ｓｄが低濃度を示す値であるほど０に近い値を出力する。即ち画像信号Ｓｄが高濃度を示す値であれば、信号Ｓ５に含まれる信号Ｓｈの割合が大きくなり（信号Ｓｈの重みが大きくなり）、画像信号Ｓｄが低濃度を示す値であれば、信号Ｓ５に含まれる信号Ｓｈの割合が小さくなる。
【００３８】
更に減算器４６は信号Ｓ５と画像信号Ｓｄを入力して減算処理を行い、画像信号Ｓｐ＝Ｓｄ−Ｓ５＝Ｓｄ−｛ｆｄ（Ｓｄ）×（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。つまり、図８（ａ）に示したフローチャートに対応する式は以下の通りとなる。
【００３９】
Ｓｐ＝Ｓｄ−｛ｆｄ（Ｓｄ）×（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝・・・（３）
以上説明したように、画像信号Ｓｄにおける高濃度部分における信号Ｓｈの空間周波数成分の重みを大きくして信号Ｓｈｌと加算し、加算後の信号を画像信号Ｓｄから減算して画像信号Ｓｐを得ることによって、画像信号Ｓｄに対応する画像の画質を劣化させることなく画像全体に混入されているグリッド像を抑制し、更に本来グリッド像が存在していない高濃度部分におけるグリッド像の出現をも抑えることができる。
【００４０】
尚、本実施例における図８（ｂ）の示す関数ｆｄは、画像信号Ｓｄの示す画像の濃度が高いほど、信号Ｓｈの空間周波数成分の重みを大きくすることとして示したが、周期的パターン抑制処理部４０ａに入力される前に画像信号Ｓｄに対して階調反転処理等が施されている場合、信号Ｓｈの空間周波数成分の重みの付け方も反転する。即ち、画像信号Ｓｄに対して階調反転処理が施されている場合は、画像信号Ｓｄの示す画像の濃度が低いほど、信号Ｓｈの空間周波数成分の重みを大きくすることとなる。
【００４１】
（実施例３）
実施例３では実施例１及び実施例２において説明した周期的パターン抑制処理方法を組み合わせた方法について説明する。図９は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ｄの構成を説明する為のフローチャートである。記憶部２９から読み出された画像信号Ｓｄは、横方向ＨＰＦ４１によってＨＰＦ処理が施される（信号Ｓｈ）。信号Ｓｈは更に縦方向ＬＰＦ４２によってＬＰＦ処理が施される（信号Ｓｈｌ）。そして減算器４３は信号Ｓｈと信号Ｓｈｌを入力して減算処理を行って信号Ｓ１を出力し、信号Ｓ１は第１関数処理部４４に入力されて関数処理が施される。第１関数処理部４４において実行される関数処理は実施例１において説明した関数処理と同様のものである。関数処理後の信号を信号Ｓ２とする。
【００４２】
一方、画像信号Ｓｄは第２関数処理部４７よって関数処理が施される。ここで、第２関数処理部４７において実行される関数処理は実施例２において説明した関数処理と同様のものである。関数処理後の信号を信号Ｓ６とする。乗算器４８は信号Ｓ６と信号Ｓ２を入力して乗算処理を行い、信号Ｓ７＝Ｓ６×Ｓ２＝ｆｄ（Ｓｄ）×ｆｃ（Ｓｈ−Ｓｈｌ）を出力する。そして加算器４５は信号Ｓ７と信号Ｓｈｌを入力して加算処理を行い、信号Ｓ８＝Ｓ７＋Ｓｈｌ＝ｆｄ（Ｓｄ）×ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌを出力する。
【００４３】
更に減算器４６は信号Ｓ８と画像信号Ｓｄを入力して減算処理を行い、画像信号Ｓｐ＝Ｓｄ−Ｓ８＝Ｓｄ−｛ｆｄ（Ｓｄ）×ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。つまり、図９に示したフローチャートに対応する式は以下の通りとなる。
【００４４】
Ｓｐ＝Ｓｄ−｛ｆｄ（Ｓｄ）×ｆｃ（Ｓｈ−Ｓｈｌ）＋Ｓｈｌ｝・・・（４）
即ち、式（４）は、実施例１において説明した式（１）と実施例２において説明した式（３）とを組み合わせた形となっている。尚、フローチャートを示していないが、式（２）と式（３）とを組み合わせてもよいことは勿論である。
【００４５】
以上説明したように、信号Ｓｈと信号Ｓｈｌのコントラスト差と、画像信号Ｓｄの濃度とに応じて画像信号Ｓｄから減算する空間周波数成分における信号Ｓｈの重みを変化させることによって、画像全体に混入されているグリッド像を抑制しつつ、本来グリッド像が存在していない高濃度部分におけるグリッド像の出現をも抑えることができる。更に、式（４）における画像信号Ｓｐに対応する画像は、実施例１及び実施例２における画像信号Ｓｐに対応する画像比べて、信号Ｓｈと信号Ｓｈｌの空間周波数の重みが切り替わる境界部分が滑らかに表現された画像となる。
【００４６】
（実施例４）
実施例３では、実施例１と実施例２とを組み合わせた周期的パターン抑制処理方法について説明した。実施例４では、コントラスト低下等の画像劣化を低減させるために、縦方向ＬＰＦ４２に用いられている遮断周波数より高い遮断周波数を用いたＬＰＦを更に用いた場合の周期的パターン抑制処理方法について説明する。
【００４７】
図１０は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ｅの構成を説明する為のフローチャートである。記憶部２９から読み出された画像信号Ｓｄは、横方向ＨＰＦ４１によってＨＰＦ処理が施される（信号Ｓｈ）。信号Ｓｈは更に縦方向ＬＰＦ４２によってＬＰＦ処理が施される（信号Ｓｈｌ）。また、信号Ｓｈは縦方向ＬＰＦ４９によってＬＰＦ処理が施される（信号Ｓｈｌｌ）。ここで、縦方向ＬＰＦ４９において用いられる遮断周波数は、縦方向ＬＰＦ４２において用いられる遮断周波数よりも高い周波数であることとする。そして減算器５０は信号Ｓｈｌｌと信号Ｓｈｌを入力して減算処理を行って、信号Ｓ９を出力する。信号Ｓ９は第１関数処理部４４に入力されて関数処理が施される。第１関数処理部４４において実行される関数処理は実施例１において説明した関数処理と同様のものである。関数処理後の信号を信号Ｓ１０とする。
【００４８】
一方、画像信号Ｓｄは第２関数処理部４７よって関数処理が施される。ここで、第２関数処理部４７において実行される関数処理は実施例２において説明した関数処理と同様のものである。関数処理後の信号を信号Ｓ６とする。乗算器４８は信号Ｓ６と信号Ｓ１０を入力して乗算処理を行い、信号Ｓ１１＝Ｓ６×Ｓ１０＝ｆｄ（Ｓｄ）×ｆｃ（Ｓｈｌｌ−Ｓｈｌ）を出力する。そして加算器４５は信号Ｓ１１と信号Ｓｈｌを入力して加算処理を行い、信号Ｓ１２＝Ｓ１１＋Ｓｈｌ＝ｆｄ（Ｓｄ）×ｆｃ（Ｓｈｌｌ−Ｓｈｌ）＋Ｓｈｌを出力する。更に減算器４６は信号Ｓ１２と画像信号Ｓｄを入力して減算処理を行い、画像信号Ｓｐ＝Ｓｄ−Ｓ１２＝Ｓｄ−｛ｆｄ（Ｓｄ）×ｆｃ（Ｓｈｌｌ−Ｓｈｌ）＋Ｓｈｌ｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。つまり、図１０に示したフローチャートに対応する式は以下の通りとなる。
【００４９】
Ｓｐ＝Ｓｄ−｛ｆｄ（Ｓｄ）×ｆｃ（Ｓｈｌｌ−Ｓｈｌ）＋Ｓｈｌ｝・・・（５）
また、図１１は周期的パターン抑制処理部の一例である周期的パターン抑制処理部４０ｆの構成を説明する為のフローチャートであり、周期的パターン抑制処理部４０ｅの変形例である。記憶部２９から読み出された画像信号Ｓｄは、横方向ＨＰＦ４１によってＨＰＦ処理が施される（信号Ｓｈ）。信号Ｓｈは更に縦方向ＬＰＦ４２によってＬＰＦ処理が施される（信号Ｓｈｌ）。また、信号Ｓｈは縦方向ＬＰＦ４９によってＬＰＦ処理が施される（信号Ｓｈｌｌ）。そして減算器５０は信号Ｓｈｌｌと信号Ｓｈｌを入力して減算処理を行い、信号Ｓ９を出力する。信号Ｓ９は第１関数処理部４４に入力されて関数処理が施される。関数処理後の信号を信号Ｓ１０とする。
【００５０】
更に減算器５１は信号Ｓｈｌｌと信号Ｓｈを入力して減算処理を行って信号Ｓ１３を出力する。信号Ｓ１３は第１関数処理部４４に入力されて関数処理が施される。関数処理後の信号を信号Ｓ１４とする。更に加算器５２は信号Ｓ１０と信号Ｓ１４を入力して加算処理を行う。加算処理後の信号を信号Ｓ１５とする。ここで信号Ｓ１５＝ｆｃ（Ｓｈ−Ｓｈｌｌ）＋ｆｃ（Ｓｈｌｌ−Ｓｈｌ）である。
【００５１】
一方、画像信号Ｓｄは第２関数処理部４７よって関数処理が施される。関数処理後の信号を信号Ｓ６とする。乗算器４８は信号Ｓ６と信号Ｓ１５を入力して乗算処理を行い、信号Ｓ１６＝Ｓ６×Ｓ１５＝ｆｄ（Ｓｄ）×｛ｆｃ（Ｓｈ−Ｓｈｌｌ）＋ｆｃ（Ｓｈｌｌ−Ｓｈｌ）｝を出力する。そして加算器４５信号Ｓ１６と信号Ｓｈｌを入力して加算処理を行い、信号Ｓ１７＝Ｓ１６＋Ｓｈｌ＝ｆｄ（Ｓｄ）×｛ｆｃ（Ｓｈ−Ｓｈｌｌ）＋ｆｃ（Ｓｈｌｌ−Ｓｈｌ）｝＋Ｓｈｌを出力する。更に減算器４６が信号Ｓ１７と画像信号Ｓｄを入力して減算処理を行い、画像信号Ｓｐ＝Ｓｄ−Ｓ１７＝Ｓｄ−｛ｆｄ（Ｓｄ）×｛ｆｃ（Ｓｈ−Ｓｈｌｌ）＋ｆｃ（Ｓｈｌｌ−Ｓｈｌ）｝＋Ｓｈｌ｝を出力する。この画像信号Ｓｐに対応する画像がモニタ等の表示手段（不図示）に表示され、又はプリンタ等の出力手段（不図示）に出力される。つまり、図１１に示したフローチャートに対応する式は以下の通りとなる。
【００５２】
Ｓｐ＝Ｓｄ−｛ｆｄ（Ｓｄ）×｛ｆｃ（Ｓｈ−Ｓｈｌｌ）＋ｆｃ（Ｓｈｌｌ−Ｓｈｌ）｝＋Ｓｈｌ｝・・・（６）
以上説明したように、信号Ｓｈｌと信号Ｓｈｌｌのコントラスト差と、信号Ｓｈｌｌと信号Ｓｈのコントラスト差と、画像信号Ｓｄの濃度とに応じて画像信号Ｓｄから減算する空間周波数成分における信号Ｓｈの重みを変化させることによって、画像全体に混入されているグリッド像を抑制しつつ、本来グリッド像が存在していない高濃度部分におけるグリッド像の出現をも抑えることができる。更に、縦方向ＬＰＦ４２に用いられている遮断周波数よりも高い遮断周波数を用いた縦方向ＬＰＦ４９を用いて画像信号Ｓｐを求めることによって、画像ボケ等の画像劣化を低減させることができる。
【図面の簡単な説明】
【００５３】
【図１】放射線画像撮影装置の概略図
【図２】蓄積性蛍光体シートに蓄積記録された被写体像とグリッド像の放射線画像の一例
【図３】放射線画像読取装置の斜視図と機能ブロック図を組み合わせた図
【図４】蓄積性蛍光体シートの主走査方向及び副走査方向を示す図
【図５】実施例１における周期的パターン抑制処理部の構成を説明する為のフローチャートと関数ｆｃのグラフ
【図６】実施例１における周期的パターン抑制処理部の構成を説明する為のフローチャートと関数ｆｃのグラフ
【図７】従来の周期的パターン抑制処理方法を説明する為の図
【図８】実施例２における周期的パターン抑制処理部の構成を説明する為のフローチャートと関数ｆｄのグラフ
【図９】実施例３における周期的パターン抑制処理部の構成を説明する為のフローチャート
【図１０】実施例４における周期的パターン抑制処理部の構成を説明する為のフローチャート
【図１１】実施例４における周期的パターン抑制処理部の構成を説明する為のフローチャート
【符号の説明】
【００５４】
４１          横方向ＨＰＦ
４２、４９    縦方向ＬＰＦ
４３、４６、５０、５１、５３    減算器
４４          第１関数処理部
４５、５２    加算器
４７          第２関数処理部
４８          乗算器
５４          第３関数処理部
【図１】
0004823737.tif 000002
【図２】
0004823737.tif 000003
【図３】
0004823737.tif 000004
【図４】
0004823737.tif 000005
【図５】
0004823737.tif 000006
【図６】
0004823737.tif 000007
【図７】
0004823737.tif 000008
【図８】
0004823737.tif 000009
【図９】
0004823737.tif 000010
【図１０】
0004823737.tif 000011
【図１１】
0004823737.tif 000012