JP2001273490A - 周期的パターン抑制処理方法および装置 - Google Patents

Abstract

(57)【要約】 【課題】 静止グリッド等に起因する周期的パターンの
空間周波数成分を抑制する周期的パターン抑制方法およ
び装置において、鮮鋭度を低下させないようにする。 【解決手段】 ウェーブレット変換部３２において、グ
リッドピッチに対応する空間周波数以上のレスポンスが
略ゼロとなる特性を有するローパスフィルタを用いて２
次元ウェーブレット変換を施して、ウェーブレット変換
係数信号を得る。抑制手段３６は、方向判定手段３５に
よるグリッド方向の判定結果に基づいて、各信号のうち
のグリッド成分を含む信号（縦グリッド時にはＨＬ１）
に対して、さらにグリッド方向について１次元ウェーブ
レット変換を繰り返し施した後、変換係数信号のうちの
低周波側の成分をゼロにし、このゼロにした信号および
残りの信号に対して１次元逆ウェーブレット変換を施
す。逆ウェーブレット変換部３４において、グリッド成
分が抑制された信号を用いて元の画像を復元する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、周期的パターン抑
制処理方法および装置に関し、より詳細には、例えば、
静止グリッドを使用して撮影を行うことにより得られ
る、該静止グリッドに対応する縞模様（モアレ成分を含
む）が被写体像に重畳された画像から、前記縞模様を低
減・除去する方法および装置に関するものである。

【０００２】

【従来の技術】従来より、Ｘ線などの放射線を照射する
とこの放射線エネルギの一部が蓄積され、その後可視光
等の励起光を照射すると蓄積されたエネルギに応じた光
量の輝尽発光光を放射する蓄積性蛍光体（輝尽性蛍光
体）を利用して、人体などの被写体の放射線画像を一旦
シート状の蓄積性蛍光体に撮影記録し、蓄積性蛍光体シ
ートをレーザ光などの励起光で走査して輝尽発光光を発
生させ、得られた輝尽発光光をフォトマルチプライヤな
どの読取手段により光電的に読み取って画像信号を得、
この画像信号に基づいて被写体の放射線画像を写真感光
材料などの記録材料やＣＲＴなどに可視像として出力さ
せる放射線記録再生システムが提案されている（特開昭
55-12429号、同56-11395号、同55-163472 号、同56-164
645号、同55-116340 号など）。

【０００３】上述した蓄積性蛍光体シートなどに被写体
の放射線画像を撮影記録する際に、被写体により散乱さ
れた放射線がシートに照射されないように４本／mm程度
の細かなピッチで放射線の透過しない鉛などと透過しや
すいアルミニウムや木材などとが交互に配置された静止
グリッドを被写体とシートとの間に配置して撮影を行う
ことがある。この静止グリッドを用いて撮影を行うと被
写体により散乱された放射線がシートに照射されにくく
なるため、被写体の放射線画像のコントラストを向上さ
せることができるが、被写体像とともに静止グリッドに
対応した細かな縞模様状のグリッド像が記録される。

【０００４】このため、本願出願人は、特願平10−1647
37号などにおいて、静止グリッドの縞模様に対応する空
間周波数成分を除去するフィルタリング処理を施すこと
により、縞模様が低減した観察しやすい画像を得る方法
を提案している。この方法においては、例えば、静止グ
リッドの格子配列ピッチ（以下グリッドピッチという）
が４本／mmである場合には、4.0cycle／mm付近の空間周
波数帯域に縞模様が現れるため、この周波数帯域のレス
ポンスを除去あるいは低減するフィルタによりフィルタ
リング処理を施して、縞模様を除去するようにしてい
る。

【０００５】

【発明が解決しようとする課題】しかしながら、上記特
願平10−164737号に記載の方法で用いているフィルタ
は、使用される静止グリッドのグリッドピッチに対応す
る空間周波数付近だけでなく、それ以上の高周波成分の
レスポンスをもゼロとするものであり、静止グリッドの
縞模様だけでなく、本来の画像に含まれる高周波成分を
も抑制・除去してしまうため、鮮鋭度が低下した画像と
なるという問題がある。

【０００６】本発明は上記事情に鑑みなされたものであ
り、鮮鋭度を低下させずに、例えば、静止グリッド等に
起因して画像中に表われる縞模様等の周期的パターンを
目立たなくすることができる周期的パターン抑制処理方
法および装置を提供することを目的とするものである。

【０００７】

【課題を解決するための手段】本発明による周期的パタ
ーン抑制処理方法は、画像信号中の周期的パターンを形
成する空間周波数成分を低減する周期的パターン抑制処
理方法であって、実空間領域で表された前記画像信号を
周波数領域で取り扱うことのできる変換画像信号に変換
する第１のステップと、該変換画像信号のうち、周期的
パターンの配列方向近傍のみにおける少なくとも該周期
的パターンの周波数に対応する空間周波数の成分を含む
所望の周波数範囲の画像信号を低減する第２のステップ
とを有することを特徴とするものである。

【０００８】ここで、「周期的パターン」とは、画像中
に周期的に表われるパターンを広く意味するものであ
り、例えば、静止グリッドに起因して画像中に発生する
縞模様や、テレビにおいて撮像系の特性により画像に発
生するモアレ的な像など種々のパターンを示すものであ
る。

【０００９】また、本発明による周期的パターン抑制処
理方法は、静止グリッドを使用した撮影により得られた
画像信号中の、前記静止グリッドに起因する空間周波数
成分を低減する周期的パターン抑制処理方法であって、
実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する第１のステッ
プと、該変換画像信号のうち、静止グリッドの格子配列
方向近傍のみにおける少なくとも該静止グリッドの格子
配列周波数に対応する空間周波数の成分を含む所望の周
波数範囲の画像信号を低減する第２のステップとを有す
るものとすることもできる。

【００１０】ここで「周波数領域で取り扱うことのでき
る変換画像信号に変換」するとは、画像信号に含まれる
所望の周波数成分について所望の処理ができる信号に変
換することを意味し、特に本願発明においては、実空間
領域で表された画像信号に含まれる前記静止グリッドに
対応する周波数成分を抑制する処理ができる信号に変換
することを意味する。この変換処理としては、例えば、
実空間領域で表現された画像信号をフーリエ変換して周
波数領域表現された画像信号（周波数スペクトル）にし
たり、あるいは実空間領域表現された画像信号を所定の
帯域にサブバンド分割（ウェーブレット変換やラプラシ
アンピラミッド展開などの場合）する多重解像度変換な
どがある。

【００１１】「静止グリッドに起因する空間周波数成
分」とは、静止グリッドそのものの空間周波数成分だけ
でなく、ナイキスト周波数以下のサンプリング周期でサ
ンプリングしたり、あるいは縮小処理することにより、
静止グリッドに起因して発生するモアレ成分も含む意味
である。なお「格子配列周波数に対応する空間周波数の
成分」としたのはこのためであって、グリッド周波数と
同じ空間周波数成分だけでなく、これと関係する前記モ
アレ周波数と同じ空間周波数成分を含む意味である。

【００１２】「少なくとも静止グリッドの格子配列周波
数に対応する空間周波数の成分を含む所望の周波数範囲
の画像信号を低減する」とは、グリッド周波数に対応す
る空間周波数の成分およびその近傍の空間周波数成分を
抑制する処理を施すことを意味し、従来例のように、グ
リッド周波数に対応する空間周波数成分以上の略全ての
高周波成分を抑制するものではなく、鮮鋭度に影響を与
える、前記グリッド周波数に対応する空間周波数成分以
上の高周波成分（例えばナイキスト周波数に対応する空
間周波数成分近傍）を可能な限り抑制しないことを意味
する。

【００１３】「静止グリッドの格子配列方向近傍のみに
おける」とは、従来例のように、グリッド方向に拘わら
ず抑制処理を施すのとは異なり、静止グリッドの格子配
列方向あるいはその近傍方向について抑制処理を施し、
それ以外の方向については抑制処理を施さないことを意
味する。なお、ここでいう「静止グリッド」は、後述す
るように、実際に使用する静止グリッドのみに限らず、
使用が予定される（グリッド方向がそれぞれ異なる）各
静止グリッドであってもよく、これらのぞれぞれについ
ての「格子配列方向近傍のみ」であればよい。換言する
と、通常の装置において、使用される静止グリッドのグ
リッド方向はある程度決まっており、例えば主走査方向
あるいは副走査方向であるので、互いに直交する主走査
方向および副走査方向において、グリッド周波数に対応
する成分を抑制する処理を施すものとしてもよい。より
判りやすく言えば、使用が予定される静止グリッドのグ
リッド方向の略全てについて、前記抑制処理を施しても
よいということである。例えば、主走査方向をｖ軸、副
走査方向をｕ軸とするフーリエ空間上において、ｖ軸上
の近傍で（ｖ軸と直交するｕ軸の正負方向に多少の幅を
持って）且つｖ軸方向のグリッド周波数を含む所望の空
間周波数範囲の成分と、ｕ軸上の近傍で（ｕ軸と直交す
るｖ軸の正負方向に多少の幅を持って）且つｕ軸方向の
グリッド周波数を含む所望の空間周波数範囲の成分のみ
を抑制するとよい。この場合、一方の静止グリッドの格
子長さ方向は、他方の静止グリッドの格子配列方向と同
じ方向となるので、結果的には、互いに直交する静止グ
リッドの格子配列方向および格子長さ方向のいずれにつ
いても抑制処理を施すことになる。

【００１４】なお、本発明による周期的パターン抑制処
理方法においては、実際に使用する静止グリッドのみに
ついて前記抑制処理を行なうのが好ましい。ここで、実
際に使用する静止グリッドのみについて前記抑制処理を
行なうということは、実際に使用する静止グリッドの格
子配列方向あるいはその近傍方向のみについて抑制処理
を施し、それ以外の方向（その静止グリッドの格子長さ
方向も含む）については抑制処理を施さないことを意味
する。例えば、縦方向の静止グリッドを用いた場合に
は、上述の例において、ｖ軸上の近傍で且つｖ軸方向の
グリッド周波数を含む所望の空間周波数範囲の成分のみ
を抑制することを意味する。

【００１５】本発明による周期的パターン抑制処理方法
においては、前記第１のステップを、静止グリッドの空
間周波数以上のレスポンスが略ゼロとなるように帯域分
割するローパスフィルタを用いて前記画像信号に対して
２次元ウェーブレット変換を施すことにより複数の変換
画像信号を得るものとすると共に、前記第２のステップ
を、複数の変換画像信号のうちの、格子配列周波数に対
応する空間周波数の成分を含む信号について、さらに、
所定周波数以下の成分を低減する処理を施し、その後、
逆ウェーブレット変換を行なうものとするのが望まし
い。

【００１６】この場合、第２のステップを、格子配列周
波数に対応する空間周波数の成分を含む変換画像信号
（ウェーブレット変換係数信号）に対して、該静止グリ
ッドの格子長さ方向（グリッド方向）について、所定の
帯域分割フィルタを用いて１次元ウェーブレット変換を
所定回数だけ再帰的に繰り返し施し、該１次元ウェーブ
レット変換により得られた複数の画像信号（ウェーブレ
ット変換係数信号）のうちの、低周波成分の変換係数を
ゼロにし、その後、１次元逆ウェーブレット変換を施す
ことにより、前記所定周波数以下の成分を低減するもの
とすればより望ましい。

【００１７】ここで「所定の帯域分割フィルタ」は、必
ずしも静止グリッドの空間周波数以上のレスポンスが略
ゼロとなるように帯域分割するローパスフィルタである
必要はない。

【００１８】また、静止グリッドの格子長さ方向につい
て１次元ウェーブレット変換を施すに際しては、例えば
縦方向グリッドのときには２次元ウェーブレット変換に
より得られたＨＬ成分に対して副走査方向に、横方向グ
リッドのときにはＬＨ成分に対して主走査方向に、それ
ぞれ１次元ウェーブレット変換を施すとよい。一方、ク
ロスグリッドのときにはＨＨ成分に対して主走査方向の
１次元ウェーブレット変換を施した後に副走査方向の１
次元ウェーブレット変換を施す（主・副の順序が逆も
可）とよい。

【００１９】なお、このような方法に限らず、ウェーブ
レット変換により得られた画像信号のうちの、静止グリ
ッドの成分を含む画像信号についてハイパスフィルタを
用いて原画像を表す成分のみを抽出することで、静止グ
リッド成分を表す所定周波数以下の成分を低減すること
もできる。

【００２０】本発明による周期的パターン抑制処理方法
においては、第２のステップを、前記複数の変換画像信
号（ウェーブレット変換係数信号）のパワーをそれぞれ
求め、求めた各パワーがそれぞれ所定の閾値よりも大き
いか否かに基づいて、静止グリッドの格子長さ方向（グ
リッド方向）を判定し、この判定結果に基づいて、前記
所定周波数以下の成分を低減する処理を施すものとする
ことが望ましい。

【００２１】なお、本発明による周期的パターン抑制処
理方法においては、実際に使用する静止グリッドのみに
ついてに限らず、使用が予定された各静止グリッドにつ
いて、前記２次元ウェーブレット変換により得た変換画
像信号を用いて、前記第２のステップによる処理を施す
ものとしてもよい。

【００２２】本発明による周期的パターン抑制処理方法
においては、前記第１のステップを、実空間領域で表さ
れた画像信号に対して、静止グリッドの格子長さ方向に
ついて、所定の帯域分割フィルタを用いて１次元ウェー
ブレット変換を施すことにより複数の変換画像信号を得
るものとすると共に、前記第２のステップを、複数の変
換画像信号のうちの、静止グリッドの格子配列周波数に
対応する空間周波数の成分を含む低周波側の変換画像信
号について、さらに、所定周波数以下の成分を低減する
処理を施し、その後、逆ウェーブレット変換を行なうも
のとすることもできる。

【００２３】なお、この場合においても、実際に使用す
る静止グリッドのみについてに限らず、使用が予定され
た各静止グリッドについて、前記第１および第２の各ス
テップによる処理を施すものとすることができる。

【００２４】本発明による周期的パターン抑制処理装置
は、上記方法を実施する装置、すなわち、画像信号中の
周期的パターンを形成する空間周波数成分を低減する周
期的パターン抑制処理装置であって、実空間領域で表さ
れた画像信号を周波数領域で取り扱うことのできる変換
画像信号に変換する画像信号変換手段と、該変換画像信
号のうち、周期的パターンの配列方向近傍のみにおける
少なくとも該周期的パターンの周波数に対応する空間周
波数の成分を含む所望の周波数範囲の画像信号を低減す
る周期的パターン成分抑制手段とを備えたことを特徴と
するものである。

【００２５】また、本発明による周期的パターン抑制処
理装置は、静止グリッドを使用した撮影により得られた
画像信号中の、静止グリッドに起因する空間周波数成分
を低減する周期的パターン抑制処理装置であって、実空
間領域で表された画像信号を周波数領域で取り扱うこと
のできる変換画像信号に変換する画像信号変換手段と、
該変換画像信号のうち、静止グリッドの格子配列方向近
傍のみにおける少なくとも該静止グリッドの格子配列周
波数に対応する空間周波数の成分を含む所望の周波数範
囲の画像信号を低減する静止グリッド成分抑制手段とを
備えたものとすることもできる。

【００２６】本発明による周期的パターン抑制処理装置
においては、画像信号変換手段を、静止グリッドの空間
周波数以上のレスポンスが略ゼロとなるように帯域分割
するローパスフィルタを用いて前記画像信号に対して２
次元ウェーブレット変換を施すことにより複数の変換画
像信号を得るものとし、静止グリッド成分抑制手段を、
前記複数の変換画像信号のうちの、静止グリッドの格子
配列周波数に対応する空間周波数の成分を含む画像信号
について、さらに、所定周波数以下の成分を低減する処
理を施し、その後、逆ウェーブレット変換を行なうもの
とするのが望ましい。

【００２７】また、この場合、静止グリッド成分抑制手
段を、格子配列周波数に対応する空間周波数の成分を含
む変換画像信号に対して、静止グリッドの格子長さ方向
について、所定の帯域分割フィルタを用いて１次元ウェ
ーブレット変換を所定回数だけ再帰的に繰り返し施し、
該１次元ウェーブレット変換により得られた複数の画像
信号のうちの、低周波成分の変換係数をゼロにし、その
後、１次元逆ウェーブレット変換を施すことにより、前
記所定周波数以下の成分を低減するものするととより望
ましい。

【００２８】また、本発明による周期的パターン抑制処
理装置においては、ウェーブレット変換を施すことによ
り得られた複数の画像信号のパワーをそれぞれ求め、求
めた各パワーがそれぞれ所定の閾値よりも大きいか否か
に基づいて、静止グリッドの格子長さ方向（グリッド方
向）を判定する静止グリッド方向判定手段をさらに備え
たものとすると共に、静止グリッド成分抑制手段を、静
止グリッド方向判定手段による判定結果に基づいて、所
定周波数以下の成分を低減する処理を施すものとするの
が望ましい。

【００２９】なお、上記静止グリッド成分抑制手段を、
実際に使用する静止グリッドのみについてに限らず、使
用が予定された各静止グリッドについて、所定周波数以
下の成分を低減する処理を施すものとすることができ
る。

【００３０】また、本発明による周期的パターン抑制処
理装置においては、画像信号変換手段を、実空間領域で
表された画像信号に対して、静止グリッドの格子長さ方
向について、所定の帯域分割フィルタを用いて１次元ウ
ェーブレット変換を施すことにより複数の変換画像信号
を得るものとすると共に、静止グリッド成分抑制手段
を、前記複数の変換画像信号のうちの、静止グリッドの
格子配列周波数に対応する空間周波数の成分を含む低周
波側の画像信号について、さらに、所定周波数以下の成
分を低減する処理を施し、その後、逆ウェーブレット変
換を行なうものとすることもできる。

【００３１】なお、この場合、画像信号変換手段を、使
用が予定された各静止グリッドの格子長さ方向につい
て、１次元ウェーブレット変換を施すものとすると共
に、静止グリッド成分抑制手段についても、使用が予定
された各静止グリッドについて低減する処理および逆ウ
ェーブレット変換を行なうものとすることができる。

【００３２】

【発明の効果】本発明による周期的パターン抑制処理方
法および装置によれば、周期的パターンが表われた画像
の画像信号を周波数領域で取り扱うことのできる画像信
号に変換し、該周波数領域で取り扱うことのできる画像
信号のうち、周期的パターンの略配列方向における周期
的パターンの周波数と同じ空間周波数の成分を含む所望
の周波数範囲の画像信号を低減するようにしたので、画
像中に表われた周期的パターンを目立たなくすることが
できる。

【００３３】また、周期的パターンを静止グリッドに起
因したものとした場合には、静止グリッドを使用した撮
影により得られた実空間領域で表された画像信号を周波
数領域で取り扱うことのできる画像信号に変換し、該周
波数領域で取り扱うことのできる画像信号のうち、静止
グリッドの略格子配列方向におけるグリッド周波数と同
じ空間周波数の成分を含む所望の周波数範囲の画像信号
を低減するようにしたので、静止グリッドに起因する縞
模様が目立たなくなる。一方、静止グリッドに起因する
空間周波数成分以上の全ての高周波成分を抑制すること
はないので、その分だけ、最終的に得られる画像信号に
はナイキスト周波数に対応する空間周波数成分近傍の成
分が含まれることとなり、鮮鋭度を低下させることがな
い。

【００３４】また、従来のようにグリッド方向に拘わら
ずフーリエ空間上の全方向に抑制処理を行なうのとは異
なり、静止グリッドの格子配列方向についてのみグリッ
ド成分が除去されるので、必要以上に画像を劣化させる
ことがない。

【００３５】また、グリッド方向が判る場合には、静止
グリッドの格子配列方向のみについて前述の抑制処理を
施すようにすれば、一層必要以上に画像を劣化させるこ
とがなくなる。

【００３６】また、使用が予定される静止グリッド（互
いにグリッド方向の異なるもに）の各格子配列方向につ
いて、それぞれ前述の抑制処理を施すようにすれば、実
際に使用されるグリッド方向を気にすることなく（２次
元ウェーブレット変換を行なうものの場合には特にその
効果が大きい）、上記抑制処理の効果を得ることができ
る。

【００３７】また、静止グリッドの空間周波数以上のレ
スポンスが略ゼロとなるように帯域分割するローパスフ
ィルタを用いて画像信号に対して２次元ウェーブレット
変換を施し、該ウェーブレット変換により得られた画像
信号のうちの、静止グリッドの成分を含む画像信号につ
いて、さらに、所定周波数以下の成分を低減する処理を
施すようにすれば、この所定周波数以下の成分が低減さ
れた信号を用いて元の画像（解像度レベル０）を復元し
たときには、静止グリッドの成分が低減された画像を復
元することができる。

【００３８】また、この初段のウェーブレット変換によ
り得られる低解像度信号成分には、静止グリッドの空間
周波数成分が含まれなくなるので、以後、さらにウェー
ブレット変換を多段階に亘って施して低解像度画像信号
を得ても、各ウェーブレット変換係数信号には静止グリ
ッドの空間周波数成分が含まれない。したがって、画像
の復元に際して、解像度レベル０まで復元することな
く、途中の解像度レベルまで復元しても、レベル１以降
の低解像度信号にはグリッド成分が含まれないので、縮
小画像（低解像度画像）において、静止グリッドに起因
するモアレが発生するということもない。つまり、画像
を任意に拡大あるいは縮小処理しても、モアレを発生す
る心配がない。

【００３９】また、ウェーブレット変換を行なって、周
波数領域で取り扱うことのできる画像信号を得る場合に
おいて、静止グリッドの成分を含む画像信号（ウェーブ
レット変換係数信号）に対して、グリッド方向につい
て、所定の帯域分割フィルタを用いて１次元ウェーブレ
ット変換を所定回数だけ再帰的に繰り返し施し、該１次
元ウェーブレット変換により得られた複数の画像信号の
うちの、低周波成分の変換係数をゼロにするようにすれ
ば、繰り返しの回数が多いほど、復元された画像におい
ては、静止グリッドの成分のみが抑制された画像にする
ことができ、一層鮮鋭度の高い画像を得ることができ
る。

【００４０】また、ウェーブレット変換を施すことによ
り得られた複数の画像信号のパワーをそれぞれ求め、求
めた各パワーがそれぞれ所定の閾値よりも大きいか否か
に基づいてグリッド方向を判定するようにすれば、静止
グリッド成分の抑制処理に際して、どのようなグリッド
方向のものを使用して撮影を行ったかを予め知る必要が
ない。

【００４１】また、静止グリッドを使用した撮影により
得られた実空間領域で表された画像信号に対して、静止
グリッドの格子長さ方向について、所定の帯域分割フィ
ルタを用いて１次元ウェーブレット変換を施し、該ウェ
ーブレット変換により得られた複数の画像信号のうち
の、静止グリッドの成分を含む低周波側の画像信号につ
いて、さらに、所定周波数以下の成分を低減する処理を
施すことによっても、静止グリッドに起因する空間周波
数成分を含む所望の周波数範囲の画像信号を低減するこ
とができ、この静止グリッドに起因する成分が低減され
た信号を用いて画像を復元したときには、復元の解像度
レベルに拘わらず、グリッドの成分が低減された画像を
復元することができる。

【００４２】

【発明の実施の形態】以下、図面を参照して本発明の実
施の形態について詳細に説明する。なお、以下の説明で
は、上記特開昭55-12492号や特開昭56-11395号などに記
載されいているように、記録シートとしての蓄積性蛍光
体シートを利用した放射線画像情報記録再生システムに
おいて、前記シートに記録された人体の放射線画像をレ
ーザビーム走査によりデジタル画像信号として読み取る
態様を例に説明する。

【００４３】図１は、放射線画像撮影装置の概略を示し
た図である。図１に示すように、放射線源１から放射さ
れた放射線２は、被写体３を経由しさらにグリッド４に
到達する。グリッド４は、放射線２を吸収する鉛４ａ
と、放射線２を透過するアルミニウム４ｂとが、放射線
源１から発せられた放射線２がアルミニウム４ｂを経由
して、シート１１に真っ直ぐに入射するように、位置に
応じて多少傾きをもって（図１参照）、４本／mmのピッ
チで交互に配置されているものである。このため放射線
源１から発せられ被写体３を真っ直ぐに透過した放射線
２は、鉛４ａに吸収されて遮ぎられる一方、アルミニウ
ム４ｂを透過してシート１１を照射し、シート１１には
被写体像とともに４本／mmの縞模様状のグリッド像が蓄
積記録される。一方、被写体３内で散乱された散乱放射
線２ａはグリッド４の前記傾きに対して斜めに入射する
ため、アルミニウム４ｂ部に入射したものもグリッド４
内部で鉛４ａに吸収され、またはグリッド４の表面で反
射されるため、シート１１には照射されず、したがって
シート１１には散乱放射線２ａの照射の少ない鮮明な放
射線画像が蓄積記録される。

【００４４】図２は、グリッド４を使用して撮影を行う
ことによりシート１１に蓄積記録された、被写体像（図
の斜線部）５と縞模様状のグリッド像６を示した図であ
る。このようにシート１１には被写体像５とグリッド像
６とが重畳された放射線画像が記録される。

【００４５】図３は、本発明による周期的パターン抑制
処理方法を実施する周期的パターン抑制処理装置の一実
施の形態を適用した、放射線画像読取装置の斜視図であ
る。

【００４６】読取部１０の所定位置にセットされた放射
線画像が記録されたシート１１は、不図示の駆動手段に
より駆動されるエンドレスベルトなどのシート搬送手段
１５により、走査ピッチ１０本／mmで矢印Ｙ方向に搬送
（副走査）される。一方、レーザー光源１６から発せら
れた光ビーム１７はモータ２４により駆動され矢印方向
に高速回転する回転多面鏡１８によって反射偏向され、
ｆθレンズなどの集束レンズ１９を通過した後、ミラー
２０により光路を変えてシート１１に入射し副走査の方
向（矢印Ｙ方向）と略直角な矢印Ｘ方向に主走査する。
シート１１の光ビーム１７が照射された箇所からは、蓄
積記録されている放射線画像情報に応じた光量の輝尽発
光光２１が発散され、この輝尽発光光２１は光ガイド２
２の入射端面２２ａに入射し、該ガイド２２内の内部を
全反射を繰り返して進み、射出端面２２ｂから射出して
フォトマル（光電子増倍管）２３に受光されることによ
り、放射線画像を表す輝尽発光光２１が該フォトマル２
３によって光電的に検出されて電気信号Ｓａに変換され
る。なお、画像信号Ｓには、良好な放射線可視画像を再
生出力するために必要な、所望とする範囲の空間周波数
帯のうちの最高の空間周波数（後述するナイキスト周波
数）ｆｎ＝5.0cycle／mmより低い、図２に示すグリッド
像６に関する情報である4.0cycle／mmの空間周波数帯域
の情報も含まれている。このグリッド像６に関する情報
は可視画像を観察する際にその可視画像を見にくくする
原因のひとつとなるものであり、取り除くべき情報であ
る。

【００４７】アナログ出力信号Ｓａはログアンプ２６で
対数的に増幅された後、Ａ／Ｄ変換器２８において空間
周波数ｆｓ＝10.0 cycle／mmに対応するサンプリング間
隔でサンプリングされてディジタル化され、読取密度の
高い、高密度画像を担持するディジタルの画像信号Ｓが
得られる。

【００４８】この画像信号Ｓは、図４に示すように、シ
ート１１に対して主走査方向（横方向）にレーザビーム
を走査させながらシート１１を副走査方向（縦方向）に
移動させてシート１１を２次元走査して得られた画像情
報を表すことになる。なお、このようにして得られた画
像信号Ｓには、ナイキスト周波数ｆｎ以下の情報を担持
しており、したがって図２に示すグリッド像６の情報
（4.0cycle／mm）も含まれている。なお、この実施の形
態においてはグリッド像６の情報（4.0cycle／mm）の２
倍以上の空間周波数ｆｓに対応するサンプリング間隔で
デジタル化しているので、エリアジングによるグリッド
像６のモアレは発生しない。

【００４９】画像信号Ｓは一旦記憶部２９に一旦記憶さ
れた後、画像信号処理部３０に入力され、以下のように
して画像処理が施される。

【００５０】図５は本発明による周期的パターン抑制処
理方法を実施するための画像信号処理部（周期的パター
ン抑制処理装置の一態様）３０の概略を表すブロック図
である。図５に示すように本発明による周期的パターン
抑制処理方法を実施するための画像信号処理部３０は、
記憶部２９から読み出した画像信号Ｓに対して多重解像
度分解処理を施す画像信号変換手段としての多重解像度
分解処理手段３２と、多重解像度分解処理手段３２にお
いて複数の周波数帯域に分解された各帯域制限画像信号
（バンドパス信号，サブバンド信号）のうち、所定の周
波数帯域の画像（帯域制限画像）に対して所望の（例え
ば強調処理など）処理を施す処理手段３３と、処理手段
３３により所望の処理が施された帯域制限画像信号およ
び他の周波数帯域を担持する帯域制限画像信号を復元し
て処理済画像信号を得るための復元処理手段３４と、復
元処理手段３４により復元された処理済画像信号に基づ
いて、可視画像を再生するための画像出力手段３５を有
する。

【００５１】なお、本実施の形態においては、実空間領
域で表される画像信号Ｓを周波数領域で取り扱うことの
できる画像信号に変換するに際して、ウェーブレット変
換処理を用いた多重解像度分解処理を使用するものとす
る。したがって、前記多重解像度分解処理手段３２は、
放射線画像を表す画像信号Ｓに対してウェーブレット変
換を施すウェーブレット変換手段として機能し、前記復
元処理手段３４は逆ウェーブレット変換処理を施して処
理済画像信号を得る逆ウェーブレット変換手段として機
能する。以下の説明においては、多重解像度分解処理手
段３２をウェーブレット変換部３２といい、復元処理手
段３４を逆ウェーブレット変換部３４という。

【００５２】図６はウェーブレット変換部３２において
行われる多重解像度分解処理としてのウェーブレット変
換処理を説明するための機能ブロック図であり、図７は
各ウェーブレット変換手段３２ａの詳細を示したブロッ
ク図である。図６に示すように、ウェーブレット変換部
３２には解像度（画素密度）レベルに応じた段数分のウ
ェーブレット変換手段３２ａが設けられ、逆ウェーブレ
ット変換部３４には、ウェーブレット変換手段３２ａと
同じ段数分の逆ウェーブレット変換手段３４ａが設けら
れる。

【００５３】また図６に示すように、画像信号処理部３
０は、ウェーブレット変換により得た信号ＨＬ１，ＬＨ
１，ＨＨ１のパワーをそれぞれ求め、求めた各パワーが
それぞれ所定の閾値ＴＨ１よりも大きいか否かに基づい
て、静止グリッド６の格子長さ方向、すなわちグリッド
方向を判定する静止グリッド方向判定手段（以下方向判
定手段という）３５と、信号ＨＬ，ＬＨ，ＨＨのうち、
静止グリッド４の格子配列方向（グリッド方向と直交す
る方向）における、少なくとも該静止グリッド４に起因
する空間周波数成分を含む所望の周波数範囲の画像信号
を低減する静止グリッド成分抑制手段３６とを備えてい
る。

【００５４】次いで本発明による周期的パターン抑制処
理方法について説明する。なお、本実施の形態において
は、ウェーブレット変換の各係数が直交する、２次元の
直交ウェーブレット変換を行なうものとする。

【００５５】図６および図７に示すように、記憶部２９
から読み出された原画像を表すデジタルの画像信号Ｓが
ウェーブレット変換部３２に入力されると、画像信号Ｓ
を原画像信号Ｓorg として、該原画像信号Ｓorg に対し
てウェーブレット変換が施される。すなわち、原画像信
号Ｓorg （信号ＬＬ０と等価）の主走査方向にウェーブ
レット関数Ｈ１，Ｇ１によりフィルタリング処理を行う
とともに、主走査方向の画素を１画素おきに間引き（図
中↓２で表す）、主走査方向の画素数を１／２にする。
ここで、関数Ｈ１はハイパスフィルタであり、関数Ｇ１
はローパスフィルタである。さらに、この画素が間引か
れた信号のそれぞれに対して副走査方向に前記関数Ｈ
１，Ｇ１によりフィルタリング処理を行うとともに、副
走査方向の画素を１画素おきに間引き、副走査方向の画
素数を１／２にして、ウェーブレット変換係数信号（以
下単に信号ともいう）ＨＨ１，ＨＬ１，ＬＨ１，ＬＬ１
を得る。ここで、信号ＬＬ１は原画像の縦横を各々１／
２に縮小した１／４縮小画像を表し、それぞれ原画像の
１／４縮小画像において、信号ＬＨ１は副走査方向（縦
方向）の高周波成分（横エッジ）を表す画像、信号ＨＬ
１は主走査方向（横方向）の高周波成分（縦エッジ）を
表す画像、信号ＨＨ１は対角方向の高周波成分（斜めエ
ッジ）を表す画像、信号ＬＬ１は原画像に対して１／２
解像度の低周波成分の画像を表すことになる。帯域分割
の基準となる周波数、すなわち、低周波成分と各高周波
成分との境目の周波数は、関数Ｈ１，Ｇ１のフィルタ特
性によって決まり、この初段での関数Ｈ１，Ｇ１による
フィルタリング処理の際には、ローパスフィルタとして
の関数Ｇ１は、静止グリッド４のグリッドピッチに対応
して、4.0cycle／mm以上の空間周波数の伝達特性（レス
ポンス）が略ゼロとなる特性を有するようにし、ハイパ
スフィルタとしての関数Ｈ１は、関数Ｇ１のローパス特
性を補うハイパス特性を有するようにする。このような
特性を有するローパスフィルタとしては、例えば特願平
10−164737号に記載のモアレ除去フィルタと同じ特性を
有するもの、すなわち静止グリッド４のグリッドピッチ
に対応する空間周波数成分に対して９７％以上の空間周
波数成分のレスポンスを５％以下に低減するものを用い
るとよく、例えば、表１に示す（１７，７）タップ、表
２に示す（１３，７）タップ、表３に示す（１５，５）
タップなどのフィルタ係数を用いたウェーブレット変換
フィルタで実現することができる。なお、各ウェーブレ
ット変換フィルタにおけるローパスフィルタの周波数応
答特性を図８に示す。

【００５６】

【表１】

【表２】

【表３】 これにより、静止グリッド４として横方向グリッドが使
用されたときにはグリッド成分が信号ＬＨ１に現れ、縦
方向グリッドが使用されたときにはグリッド成分が信号
ＨＬ１に現れ、クロスグリッドが使用されたときにはグ
リッド成分が信号ＨＨ１に現れ、グリッド方向に拘わら
ず、信号ＬＬ１においては、静止グリッド４の空間周波
数成分が十分に抑制されているので、静止グリッド４に
起因する縞模様（モアレ）が殆ど現れない。

【００５７】ウェーブレット変換係数信号ＨＨ１，ＨＬ
１，ＬＨ１，ＬＬ１のうち、グリッド成分を含んでいる
可能性のある各信号ＨＨ１，ＨＬ１，ＬＨ１が方向判定
手段３５および静止グリッド成分抑制手段３６に入力さ
れる。

【００５８】さらに、ウェーブレット変換手段３２ａに
おいて、基本ウェーブレット関数Ｈ０，Ｇ０を用いて、
信号ＬＬ１に対してウェーブレット変換が施されて、信
号ＨＨ２，ＨＬ２，ＬＨ２，ＬＬ２が得られる。ここ
で、信号ＬＬ２は原画像の縦横を各々１／４に縮小した
１／１６縮小画像を表し、それぞれ原画像の１／１６縮
小画像において、信号ＨＬ２、ＬＨ２およびＨＨ２は縦
エッジ、横エッジおよび斜めエッジ成分の画像を表すも
のとなる。なお、上述のように、信号ＬＬ１にはグリッ
ド成分が殆ど現れないので、この２段目で用いられるウ
ェーブレット関数Ｈ０，Ｇ０は、上記初段で用いられた
ウェーブレット関数Ｈ１，Ｇ１とは異なり、静止グリッ
ド４のグリッドピッチに対応して設定されるものである
必要はなく、例えば、表４に示すドーベッチ（Daubechi
es）の（９，７）タップのフィルタ係数を用いたウェー
ブレット変換フィルタとするとよい。なお、この（９，
７）タップウェーブレット変換フィルタにおけるローパ
スフィルタの周波数応答特性を図８に示す。

【００５９】

【表４】 以下、上記２段目と同様にして、各周波数帯域において
得られるウェーブレット変換係数信号ＬＬｋに対するウ
ェーブレット変換をｎ回繰り返すことによりウェーブレ
ット変換係数信号ＨＨ１〜ＨＨｎ，ＨＬ１〜ＨＬｎ，Ｌ
Ｈ１〜ＬＨｎ，ＬＬ１〜ＬＬｎを得る。ここで、ｎ回目
のウェーブレット変換により得られるウェーブレット変
換係数信号ＨＨｎ，ＨＬｎ，ＬＨｎ，ＬＬｎは、原画像
信号Ｓorg と比較して主副各方向の画素数が(1/2）ｎ
となった(1/2）２ｎ縮小画像を表し、各ウェーブレット
変換係数信号ＨＨｎ，ＨＬｎ，ＬＨｎ，ＬＬｎはｎが大
きいほど周波数帯域が低くなる。このように、ウェーブ
レット変換係数信号ＨＨｋ，ＨＬｋ，ＬＨｋ，ＬＬｋ
（ｋは解像度レベルであって、ｋ＝１〜ｎの整数）は、
それぞれ、原画像信号Ｓorg の周波数範囲のうちの所定
範囲の周波数成分を担持する帯域制限画像信号となる。
また、信号ＨＨｋは原画像信号Ｓorg の主副両方向の周
波数の変化を表しｋが大きいほど低周波信号となる。ま
た信号ＨＬｋは原画像信号Ｓorg の主走査方向の周波数
の変化を表すものであり、ｋが大きいほど低周波信号と
なる。さらに信号ＬＨｋは原画像信号Ｓorg の副走査方
向の周波数の変化を表すものであり、ｋが大きいほど低
周波信号となる。

【００６０】次に、方向判定手段３５において、ウェー
ブレット変換により得られた各信号ＬＨ１，ＨＬ１，Ｈ
Ｈ１のパワーをそれぞれ求め、求めた各パワーがそれぞ
れ所定の閾値ＴＨ０よりも大きいか否かを確認し、その
結果に基づいて、静止グリッド４のグリッド方向を判定
する。具体的には、上述のように、グリッド成分は、横
方向グリッド時には信号ＬＨ１に現れ、縦方向グリッド
時には信号ＨＬ１に現れ、クロスグリッド時には信号Ｈ
Ｈ１に現れるので、信号ＬＨのパワーのみが閾値ＴＨ０
よりも大きいときには横方向グリッド、信号ＨＬのパワ
ーのみが閾値ＴＨ０よりも大きいときには縦方向グリッ
ド、信号ＨＬおよび信号ＬＨのいずれもが閾値ＴＨ０よ
りも大きいときにはクロスグリッドが使用されているも
のと判定する。この判定結果は、静止グリッド成分抑制
手段３６に入力される。

【００６１】次に、静止グリッド成分抑制手段３６にお
いて、方向判定手段３５によるグリッド方向の判定結果
に基づいて、信号ＬＨ１，ＨＬ１，ＨＨ１のうちの静止
グリッドの成分を含む信号に対して、さらに静止グリッ
ド４の格子長さ方向（グリッド方向）について、所定の
帯域分割フィルタを用いて１次元ウェーブレット変換を
所定回数だけ再帰的に繰り返し施す。例えば、図２に示
すように縦方向グリッドが使用されている場合には、Ｈ
Ｌ１に対して縦（副走査）方向の１次元ウェーブレット
変換を施して低周波側を担持する信号ＨＬ１（Ｌ１）と
高周波側を担持する信号ＨＬ１（Ｈ１）を得、次に低周
波側を担持する信号ＨＬ１（Ｌ１）に対して縦方向の１
次元ウェーブレット変換を施して低周波側を担持する信
号ＨＬ１（Ｌ２）と高周波側を担持する信号ＨＬ１（Ｈ
２）を得、以下同様に、低周波側を担持する信号ＨＬ１
（Ｌｊ）（ｊ＝１〜ｍの整数）に対して縦方向の１次元
ウェーブレット変換をｍ回繰り返すことにより、信号Ｈ
Ｌ１に対して複数段（ｍ段）の縦方向の１次元ウェーブ
レット変換を施して、ウェーブレット変換係数信号ＨＬ
１（Ｈ１）〜ＨＬ１（Ｈｍ），ＨＬ１（Ｌｍ）を得る。

【００６２】図９は、原画像信号をウェーブレット変換
して各成分に分解した結果を示す図である。なお、図９
においては、１回目の２次元ウェーブレット変換を行な
い（同図（ａ））、その後２回目の２次元ウェーブレッ
ト変換を行なった後に、信号ＨＬ１について、さらに縦
方向の１次元ウェーブレット変換を１回施した状態まで
（同図（ｂ））を表している。

【００６３】図１０は、格子長さ方向に１次元ウェーブ
レット変換を繰り返し施すことによる効果を説明する模
式図である。

【００６４】ウェーブレット変換係数信号ＨＬ１は、副
走査方向についての低周波成分を表し、図１０に示すよ
うに、信号ＨＬ１は縦方向グリッドの成分と原画像信号
自身が持つ縦方向を主成分とする画像信号を含むことに
なる。そして、この信号ＨＬ１に対して格子長さ方向
（本例では副走査方向）について１次元ウェーブレット
変換を施すと、縦方向グリッドを用いた場合のグリッド
成分は１次元の縦長の低周波信号と見られるために、Ｈ
Ｌ１（Ｌｊ）には縦方向グリッドの成分を含む低周波成
分、ＨＬ１（Ｈｊ）にはその他の高周波成分が含まれる
ように、帯域分割される。この１次元ウェーブレット変
換を繰り返し施すと、十分に絞り込まれた低周波信号Ｈ
Ｌ１（Ｌｍ）は、グリッドピッチの大きさに拘わらず、
グリッドの目の方向（格子長さ方向）成分を含む極めて
低周波数の成分を表し、ｍを大きくするほど信号ＨＬ１
（Ｌｍ）がグリッド成分のみを表すようになる。

【００６５】なお、実際には静止グリッドは完全には水
平・垂直あるいは水平・垂直が１対１のクロス方向には
置かれない場合が多いので、繰返しの回数ｍをあまり大
きくすると、信号ＨＬ１（Ｈｊ）にもグリッド成分が現
れるようになるので、１次元ウェーブレット変換は数回
程度繰り返すに留めるのが好まい。

【００６６】次に、信号ＨＬ１について複数段の縦方向
の１次元ウェーブレット変換を施して得られた各信号Ｈ
Ｌ１（Ｈ１）〜ＨＬ１（Ｈｍ），ＨＬ１（Ｌｍ）のう
ち、最も低周波側の成分である信号ＨＬ１（Ｌｍ）をゼ
ロにする、つまり、グリッド成分を含む信号ＨＬ１（Ｌ
ｍ）を抑制するようにする。その後、ゼロにした信号Ｈ
Ｌ１（Ｌｍ）および各信号ＨＬ１（Ｈ１）〜ＨＬ１（Ｈ
ｍ）に対して１次元逆ウェーブレット変換を施して、信
号ＨＬ１’を得る。ここで、前述のように、少なくとも
静止グリッド４の空間周波数成分を含む信号ＨＬ１（Ｌ
ｍ）が抑制されているので、グリッド成分を含む所定範
囲の空間周波数成分が低減された信号ＨＬ１’が得られ
る。

【００６７】次に、処理手段３３により、必要に応じ
て、所望の（例えば強調処理など）処理が施され、その
後、逆ウェーブレット変換部３４において、信号ＬＬ
ｎ、ＨＬｋ，ＬＨｋに対して、レベルｎからレベル１ま
で、順次逆ウェーブレット変換が施される。

【００６８】図１１は、逆ウェーブレット変換部３４の
構成を示す概略ブロック図である。図１１に示すよう
に、最低周波数帯域の信号ＨＨｎ，ＨＬｎ，ＬＨｎ，Ｌ
Ｌｎに対して逆ウェーブレット変換手段３４ａにおいて
逆ウェーブレット変換を施して信号ＬＬｎ−１を得る。

【００６９】図１２は各逆ウェーブレット変換手段３４
ａにおいて行われる逆ウェーブレット変換処理を説明す
るための機能ブロック図である。図１２に示すように信
号ＬＬｎおよび信号ＬＨｎ（ＬＨｋ）の副走査方向に対
して画素間に１画素分の間隔をあける処理を行うととも
に（図中↑２で表す）、ウェーブレット変換の際に用い
た関数Ｇ０，Ｈ０に対応する逆ウェーブレット変換関数
Ｇ０’，Ｈ０’によりフィルタリング処理を副走査方向
に施してこれらを加算し、さらに加算により得られた信
号（第１の加算信号とする）の主走査方向に対して画素
間に１画素分の間隔をあける処理を行うとともに、関数
Ｇ０’によりフィルタリング処理を主走査方向に施して
第１の信号を得る。一方、信号ＨＬｎ（ＨＬｋ）および
信号ＨＨｎ（ＨＨｋ）の副走査方向に対して画素間に１
画素分の間隔をあける処理を行うとともに、関数Ｇ
０’，Ｈ０’によりフィルタリング処理を副走査方向に
施してこれらを加算し、さらに加算により得られた信号
（第２の加算信号とする）の主走査方向に対して画素間
に１画素分の間隔をあける処理を行うとともに、関数Ｈ
０’によりフィルタリング処理を主走査方向に施して第
２の信号を得る。そして第１および第２の信号を加算し
て信号ＬＬｎ−１（ＬＬｋ−１）を得る。

【００７０】次に、信号ＨＨｎ−１，ＨＬｎ−１，ＬＨ
ｎ−１，ＬＬｎ−１に対して上記と同様に逆ウェーブレ
ット変換手段３４ａにおいて逆ウェーブレット変換を行
って、処理済み信号ＬＬｎ−２を得る。そして、以下上
記と同様にして逆ウェーブレット変換を解像度レベル１
まで繰り返すことにより信号ＬＬ１が得られる。

【００７１】そして、原画像を表す解像度レベル０で
は、信号ＬＬ１、ＬＨ１，ＨＨ１および信号ＨＬ１’に
対して、初段におけるウェーブレット変換の際に用いた
関数Ｇ１，Ｈ１に対応する逆ウェーブレット変換関数Ｇ
１’，Ｈ１’により逆ウェーブレット変換を施して画像
を再構成する。ここで、信号ＨＬ１’はグリッド成分が
低減されたものであり、再構成画像は、静止グリッド４
の成分が軽減されている。

【００７２】図１３は、本願発明の効果を説明する図で
あって、図１３（ａ）は特願平10−164737号（以下従来
例という）に記載の方法を適用した場合における周波数
応答特性を示した図、図１３（ｂ）は本願発明を適用し
た場合における周波数応答特性を示した図、図１３
（ｃ）は従来例および本願発明による方法におけるフー
リエ空間での空間周波数特性を示した図である。なお、
図１３（ｃ）は、主走査方向をｖ軸、副走査方向をｕ軸
とするフーリエ空間上において、縦方向の静止グリッド
を用いた場合を示したものである。

【００７３】従来例の方法を適用した場合には、使用さ
れる静止グリッドのグリッドピッチに対応する空間周波
数付近だけでなく、それ以上の高周波成分のレスポンス
もゼロとするフィルタを用いているので、図１３（ａ）
に示すように、グリッドピッチに対応する空間周波数以
上の高周波成分までもが抑制され、静止グリッド４の縞
模様だけでなく、本来の画像に含まれる高周波成分をも
除去・低減してしまうことになり、鮮鋭度が低下した画
像となってしまう。また、方向を考慮したフィルタリン
グ処理となっていないため、図１３（ｃ）に示すよう
に、フーリエ空間上における低周波領域（図中の中央部
の空白部）の範囲外全て（図中の斜線部）が抑制され
る、つまり縦方向グリッドによる縦模様（本当に抑制し
たい成分）だけでなく、本来の画像に含まれる横模様や
斜め模様などの、本来であれば抑制したくない高周波成
分も抑制されることになる。

【００７４】これに対して、本願発明を適用した場合に
は、初段のウェーブレット変換の際に用いられるフィル
タとして、特願平10−164737号に記載のものと同じ特性
を示すものを使用したとしても、図１３（ｂ）に示すよ
うに、静止グリッド４の空間周波数成分近傍の所定範囲
のみを抑制（カット）し、残りの高周波成分を保存する
ことができるので、静止グリッドに起因する縞模様が目
立たたず、しかも鮮鋭度の高い高画質な画像を得ること
ができる。また、信号ＨＬ１，ＬＨ１，ＨＨ１のうちの
グリッド成分を含む信号のみが抑制対象となるようにグ
リッド方向に応じて切り替えることができるので、図１
３（ｃ）に示すように、フーリエ空間上における、グリ
ッド方向の軸上（縦グリッドのときはｖ軸）の近傍で
（該軸と直交する方向に多少の幅を持って）、且つ静止
グリッド４の空間周波数成分を含む高周波側の所定範囲
のみを抑制することができるので、例えば縦方向グリッ
ドの場合には、グリッドによる縦模様だけを抑制し、本
来の画像に含まれる横模様や斜め模様などの高周波成分
を抑制することがない。

【００７５】また、初段での２次元ウェーブレット変換
により得られた変換係数信号ＨＨ１，ＨＬ１，ＬＨ１，
ＬＬ１のうち、グリッド成分を含んでいる可能性のある
各信号ＨＨ１，ＨＬ１，ＬＨ１に基づいて静止グリッド
の方向を判定することができるので、静止グリッド成分
の抑制処理に際して、どのようなグリッド方向のものを
使用して撮影を行ったかを予め知る必要もない。

【００７６】また、レベル０まで復元することなく、解
像度レベル１まで復元すれば、前述のように、信号ＬＬ
１においては静止グリッドの空間周波数成分が十分に抑
制されているので、静止グリッド４に起因する縞模様が
殆ど現れない。また、グリッド成分が抑制された信号Ｌ
Ｌ１に対してウェーブレット変換を施して得たレベル１
以降のＬＬｋにもモアレ成分が含まれないので、全ての
縮小画像において、静止グリッドに起因するモアレが発
生するということがなく、画像を任意に拡大あるいは縮
小処理しても、モアレを発生する心配がない。したがっ
て、強調処理などの所望の画像処理を施して、画像処理
が施された画像を逆多重解像度変換することにより復元
しても、いずれの解像度レベルの画像においても、復元
された処理済画像にはモアレによるアーチファクトが含
まれず、見やすい高画質の画像を提供することができ
る。

【００７７】このように、多重解像度変換としてのウェ
ーブレット変換の初段においてグリッド成分が抑制さ
れ、それ以降の解像度レベルの画像にはモアレが現れな
くなり、またウェーブレット変換を用いているので、拡
大／縮小画像を得るのにも好都合であり、ネットワーク
を介して画像信号を転送する上でも都合がよい。

【００７８】以上、本発明の好ましい実施の形態につい
て説明したが、本発明は上述した実施の形態に限定され
るものではなく、本発明の要旨を逸脱しない限りにおい
て、種々の変更を行うことが可能である。

【００７９】例えば、上記静止グリッド成分抑制手段３
６における作用説明は、静止グリッド４が縦方向グリッ
ドであるとして説明したが、グリッド方向が縦方向以外
の場合には、取り扱う信号ＨＬ１を、そのグリッド方向
に応じて信号ＬＨ１あるは信号ＨＨ１とすればよい。す
なわち、静止グリッド４が横方向グリッドである場合に
は、信号ＬＨ１に対してさらに複数段の横方向（主走査
方向）の１次元ウェーブレット変換を施し、低周波側の
信号成分をゼロにした後、１次元逆ウェーブレット変換
して画像を再構成するとよく、静止グリッド４がクロス
グリッドである場合には、信号ＨＬおよび信号ＬＨそれ
ぞれに対してさらに複数段の縦方向（信号ＨＬ）／横方
向（信号ＬＨ）の１次元ウェーブレット変換を施し、低
周波側の信号成分をゼロにした後、１次元逆ウェーブレ
ット変換して画像を再構成するとよい。

【００８０】また、信号ＨＬおよび／または信号ＬＨに
対する１次元ウェーブレット変換の代わりに、グリッド
成分を低減させるようなハイパスフィルタを用いてもよ
い。

【００８１】また、実際に使用されるグリッド方向の如
何に関わらず、使用が予定された静止グリッドの各グリ
ッド方向について上述の処理を施す、例えば図９に示し
た処理に加えて、信号ＬＨ１に対して主走査方向につい
て１次元ウェーブレット変換を繰り返し施すようにすれ
ば、図１３（ｃ）に示すｖ軸の近傍の静止グリッド４の
空間周波数成分を含む高周波側の所定範囲だけでなく、
ｕ軸の近傍の静止グリッド４の空間周波数成分を含む高
周波側の所定範囲についても抑制効果を得ることができ
る。したがって、通常は主走査方向若しくは副走査方向
のいずれかと同方向のグリッドが用いられることが多い
ので、この場合、グリッド成分は信号ＨＬ１（縦グリッ
ドのとき），ＬＨ１（横グリッドのとき）のいずれかに
含まれることになり、結果的に、グリッド方向を全く気
にすることなく、上記抑制処理の効果を得ることができ
る。また、この場合においても、無制限な方向に対して
抑制処理を施すものではなく、図１３（Ｃ）に示すよう
に、抑制処理の対象方向が使用が予定された静止グリッ
ドの各グリッド方向についてのみに限られるので、画質
に与える影響を従来よりも軽減することができる。な
お、信号ＨＨ１に対して主／副の両走査方向について１
次元ウェーブレット変換を繰り返し施すことにより、ク
ロスグリッドの使用が予定される場合にも対応できる。

【００８２】また、上記説明は、画像信号Ｓに対して２
次元ウェーブレット変換を施すものとして説明したが、
グリッド方向が予め判っている場合には、図１４に示す
ように、その方向に応じた（縦グリッドのときには副走
査方向に）１次元ウェーブレット変換を原画像信号に対
して施し、グリッド成分を含むウェーブレット変換係数
信号（図１４（Ｂ）の信号Ｌ１）に対してグリッド成分
を低減する処理を施すようにしてもよい。なお、１次元
ウェーブレット変換の際に用いられるウェーブレット関
数は、静止グリッド４のグリッドピッチに拘わらず、任
意のものであってもよく、上記Ｈ０，Ｇ０を用いること
ができる。

【００８３】図１４に示した方法とすると、グリッドの
目の方向（格子長さ方向）を表す極めて周波数の成分が
ウェーブレット変換係数信号のうちの低周波側の成分
（信号Ｌ１）に現れるので、この低周波側のウェーブレ
ット変換係数信号Ｌ１について、さらにグリッド成分を
低減する処理、例えば同方向の１次元ウェーブレット変
換を繰り返し施した後（図１４（Ｃ）参照）に最も低周
波のウェーブレット変換係数信号をゼロにしたり、ある
いはハイパスフィルタでグリッド成分を除く原画像を表
す信号成分のみを抽出することで、グリッド成分を表す
極めて周波数の成分を抑制することができる。

【００８４】したがって、このグリッド成分が低減され
た信号を用いて１次元逆ウェーブレット変換を施して画
像を復元したときには、復元の解像度レベルに拘わら
ず、グリッドの成分が低減された画像を復元することが
できる。

【００８５】なお、図１４に示した態様においても、実
際に使用されるグリッド方向の如何に関わらず、使用が
予定された静止グリッドの各グリッド方向について上述
の処理をそれぞれ施すものとすることができる。例えば
図１４に示した副走査方向の処理に加えて、主走査方向
の処理を行なうものとし、最終的に得られた２つの画像
のうち、グリッド成分の少ない方の画像を使用するよう
にすることもできる。

【００８６】また上記説明は、実空間領域で表される画
像信号Ｓを周波数領域で取り扱うことのできる画像信号
に変換するに際して、ウェーブレット変換処理を用いた
多重解像度分解処理を使用するものとして説明したが、
画像信号に含まれる静止グリッド成分を含む所望範囲の
周波数成分を低減する処理を施すことができるものであ
る限り、どのような変換方法を用いてもよい。例えばラ
プラシアンピラミッド展開を用いた多重解像度分解処理
を使用することもできる。また、実空間領域表現された
画像信号をフーリエ変換して周波数領域表現された信号
（周波数スペクトル）にしてもよい。

【図面の簡単な説明】

【図１】放射線画像撮影装置の概略を示した図

【図２】グリッド撮影により得られた放射線画像を示す
図

【図３】放射線画像読取装置の一例を示す斜視図

【図４】走査方向と読取画像との関係を示した図

【図５】本発明による周期的パターン抑制処理装置の一
態様としての画像信号処理部の概略を表すブロック図

【図６】ウェーブレット変換部において行われるウェー
ブレット変換処理を説明するための機能ブロック図

【図７】ウェーブレット変換手段の詳細を示したブロッ
ク図

【図８】ウェーブレット変換フィルタにおけるローパス
フィルタの周波数応答特性を示す図

【図９】原画像信号をウェーブレット変換して各成分に
分解した結果を示す図

【図１０】１次元ウェーブレット変換を繰り返し施すこ
とによる効果を説明する模式図

【図１１】逆ウェーブレット変換部の構成を示す概略ブ
ロック図

【図１２】逆ウェーブレット変換手段において行われる
逆ウェーブレット変換処理を説明するための機能ブロッ
ク図

【図１３】本願発明の効果を説明する図

【図１４】本願発明の他の態様を説明する図

【符号の説明】

４ 静止グリッド ５ 被写体像 ６ グリッド像 １１ 蓄積性蛍光体シート ３０ 画像信号処理部（周期的パターン抑制装置の一
態様） ３２ ウェーブレット変換部（画像信号変換手段の一
態様） ３４ 逆ウェーブレット変換部 ３５ 静止グリッド方向判定手段 ３６ 静止グリッド成分抑制手段

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁷ 識別記号 ＦＩ テーマコート゛(参考） Ｈ０４Ｎ 1/409 Ｈ０４Ｎ 1/40 １０１Ｄ

Claims (16)

【特許請求の範囲】
1. 【請求項１】 画像信号中の周期的パターンを形成する
   空間周波数成分を低減する周期的パターン抑制処理方法
   であって、 実空間領域で表された前記画像信号を周波数領域で取り
   扱うことのできる変換画像信号に変換する第１のステッ
   プと、 該変換画像信号のうち、前記周期的パターンの配列方向
   近傍のみにおける少なくとも該周期的パターンの周波数
   に対応する空間周波数の成分を含む所望の周波数範囲の
   画像信号を低減する第２のステップとを有することを特
   徴とする周期的パターン抑制処理方法。
2. 【請求項２】 静止グリッドを使用した撮影により得ら
   れた画像信号中の、前記静止グリッドに起因する空間周
   波数成分を低減する周期的パターン抑制処理方法であっ
   て、 実空間領域で表された前記画像信号を周波数領域で取り
   扱うことのできる変換画像信号に変換する第１のステッ
   プと、 該変換画像信号のうち、前記静止グリッドの格子配列方
   向近傍のみにおける少なくとも該静止グリッドの格子配
   列周波数に対応する空間周波数の成分を含む所望の周波
   数範囲の画像信号を低減する第２のステップとを有する
   ことを特徴とする周期的パターン抑制処理方法。
3. 【請求項３】 前記第１のステップが、前記静止グリッ
   ドの空間周波数以上のレスポンスが略ゼロとなるように
   帯域分割するローパスフィルタを用いて前記画像信号に
   対して２次元ウェーブレット変換を施すことにより複数
   の前記変換画像信号を得るものであり、 前記第２のステップが、前記複数の変換画像信号のうち
   の、前記格子配列周波数に対応する空間周波数の成分を
   含む信号について、さらに、所定周波数以下の成分を低
   減する処理を施し、その後、逆ウェーブレット変換を行
   なうものであることを特徴とする請求項２記載の周期的
   パターン抑制処理方法。
4. 【請求項４】 前記第２のステップが、前記格子配列周
   波数に対応する空間周波数の成分を含む変換画像信号に
   対して、該静止グリッドの格子長さ方向について、所定
   の帯域分割フィルタを用いて１次元ウェーブレット変換
   を所定回数だけ再帰的に繰り返し施し、該１次元ウェー
   ブレット変換により得られた複数の画像信号のうちの、
   低周波成分の変換係数をゼロにし、その後、１次元逆ウ
   ェーブレット変換を施すことにより、前記所定周波数以
   下の成分を低減するものであることを特徴とする請求項
   ３記載の周期的パターン抑制処理方法。
5. 【請求項５】 前記第２のステップが、前記複数の変換
   画像信号のパワーをそれぞれ求め、求めた各パワーがそ
   れぞれ所定の閾値よりも大きいか否かに基づいて、前記
   静止グリッドの格子長さ方向を判定し、この判定結果に
   基づいて、前記所定周波数以下の成分を低減する処理を
   施すものであることを特徴とする請求項３または４記載
   の周期的パターン抑制処理方法。
6. 【請求項６】 使用が予定された各静止グリッドについ
   て、前記第２のステップによる処理を施すことを特徴と
   する請求項３または４記載の周期的パターン抑制処理方
   法。
7. 【請求項７】 前記第１のステップが、前記実空間領域
   で表された画像信号に対して、前記静止グリッドの格子
   長さ方向について、所定の帯域分割フィルタを用いて１
   次元ウェーブレット変換を施すことにより複数の前記変
   換画像信号を得るものであり、 前記第２のステップが、前記複数の変換画像信号のうち
   の、前記静止グリッドの格子配列周波数に対応する空間
   周波数の成分を含む低周波側の変換画像信号について、
   さらに、所定周波数以下の成分を低減する処理を施し、
   その後、逆ウェーブレット変換を行なうものであること
   を特徴とする請求項２記載の周期的パターン抑制処理方
   法。
8. 【請求項８】 使用が予定された各静止グリッドについ
   て、前記第１および第２の各ステップによる処理を施す
   ことを特徴とする請求項７記載の周期的パターン抑制処
   理方法。
9. 【請求項９】 画像信号中の周期的パターンを形成する
   空間周波数成分を低減する周期的パターン抑制処理装置
   であって、 実空間領域で表された前記画像信号を周波数領域で取り
   扱うことのできる変換画像信号に変換する画像信号変換
   手段と、 該変換画像信号のうち、前記周期的パターンの配列方向
   近傍のみにおける少なくとも該周期的パターンの周波数
   に対応する空間周波数の成分を含む所望の周波数範囲の
   画像信号を低減する周期的パターン成分抑制手段とを備
   えたことを特徴とする周期的パターン抑制処理装置。
10. 【請求項１０】 静止グリッドを使用した撮影により得
    られた画像信号中の、前記静止グリッドに起因する空間
    周波数成分を低減する周期的パターン抑制処理装置であ
    って、 実空間領域で表された前記画像信号を周波数領域で取り
    扱うことのできる変換画像信号に変換する画像信号変換
    手段と、 該変換画像信号のうち、前記静止グリッドの格子配列方
    向近傍のみにおける少なくとも該静止グリッドの格子配
    列周波数に対応する空間周波数の成分を含む所望の周波
    数範囲の画像信号を低減する静止グリッド成分抑制手段
    とを備えたことを特徴とする周期的パターン抑制処理装
    置。
11. 【請求項１１】 前記画像信号変換手段が、前記静止グ
    リッドの空間周波数以上のレスポンスが略ゼロとなるよ
    うに帯域分割するローパスフィルタを用いて前記画像信
    号に対して２次元ウェーブレット変換を施すことにより
    複数の前記変換画像信号を得るものであり、 前記静止グリッド成分抑制手段が、前記複数の変換画像
    信号のうちの、前記静止グリッドの格子配列周波数に対
    応する空間周波数の成分を含む画像信号について、さら
    に、所定周波数以下の成分を低減する処理を施し、その
    後、逆ウェーブレット変換を行なうものであることを特
    徴とする請求項１０記載の周期的パターン抑制処理装
    置。
12. 【請求項１２】 前記静止グリッド成分抑制手段が、前
    記格子配列周波数に対応する空間周波数の成分を含む変
    換画像信号に対して、該静止グリッドの格子長さ方向に
    ついて、所定の帯域分割フィルタを用いて１次元ウェー
    ブレット変換を所定回数だけ再帰的に繰り返し施し、該
    １次元ウェーブレット変換により得られた複数の画像信
    号のうちの、低周波成分の変換係数をゼロにし、その
    後、１次元逆ウェーブレット変換を施すことにより、前
    記所定周波数以下の成分を低減するものであることを特
    徴とする請求項１１記載の周期的パターン抑制処理装
    置。
13. 【請求項１３】 前記複数の変換画像信号のパワーをそ
    れぞれ求め、求めた各パワーがそれぞれ所定の閾値より
    も大きいか否かに基づいて、前記静止グリッドの格子長
    さ方向を判定する静止グリッド方向判定手段をさらに備
    え、 前記静止グリッド成分抑制手段が、該静止グリッド方向
    判定手段による判定結果に基づいて、前記所定周波数以
    下の成分を低減する処理を施すものであることを特徴と
    する請求項１１または１２記載の周期的パターン抑制処
    理装置。
14. 【請求項１４】 前記静止グリッド成分抑制手段が、使
    用が予定された各静止グリッドについて前記所定周波数
    以下の成分を低減する処理を施すものであることを特徴
    とする請求項１１または１２記載の周期的パターン抑制
    処理装置。
15. 【請求項１５】 前記画像信号変換手段が、前記実空間
    領域で表された前記画像信号に対して、前記静止グリッ
    ドの格子長さ方向について、所定の帯域分割フィルタを
    用いて１次元ウェーブレット変換を施すことにより複数
    の前記変換画像信号を得るものであり、 前記静止グリッド成分抑制手段が、前記複数の変換画像
    信号のうちの、前記静止グリッドの格子配列周波数に対
    応する空間周波数の成分を含む低周波側の画像信号につ
    いて、さらに、所定周波数以下の成分を低減する処理を
    施し、その後、逆ウェーブレット変換を行なうものであ
    ることを特徴とする請求項１０記載の周期的パターン抑
    制処理装置。
16. 【請求項１６】 前記画像信号変換手段が、使用が予定
    された各静止グリッドの格子長さ方向についいて、前記
    １次元ウェーブレット変換を施すものであり、 前記静止グリッド成分抑制手段が、前記使用が予定され
    た各静止グリッドにつて前記低減する処理および前記逆
    ウェーブレット変換を行なうものであることを特徴とす
    る請求項１５記載の周期的パターン抑制処理装置。