JP2001273490A - 周期的パターン抑制処理方法および装置 - Google Patents
周期的パターン抑制処理方法および装置Info
- Publication number
- JP2001273490A JP2001273490A JP2000395577A JP2000395577A JP2001273490A JP 2001273490 A JP2001273490 A JP 2001273490A JP 2000395577 A JP2000395577 A JP 2000395577A JP 2000395577 A JP2000395577 A JP 2000395577A JP 2001273490 A JP2001273490 A JP 2001273490A
- Authority
- JP
- Japan
- Prior art keywords
- grid
- image signal
- frequency
- wavelet transform
- periodic pattern
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 125000004122 cyclic group Chemical group 0.000 title abstract 3
- 230000001629 suppression Effects 0.000 claims abstract description 70
- 230000004044 response Effects 0.000 claims abstract description 17
- 230000009466 transformation Effects 0.000 claims abstract 9
- 230000000737 periodic effect Effects 0.000 claims description 65
- 230000003068 static effect Effects 0.000 claims description 39
- 230000008569 process Effects 0.000 claims description 29
- 238000003672 processing method Methods 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 238000011946 reduction process Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 description 34
- 238000010586 diagram Methods 0.000 description 23
- 238000000354 decomposition reaction Methods 0.000 description 9
- 238000001914 filtration Methods 0.000 description 9
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 5
- 229910052782 aluminium Inorganic materials 0.000 description 5
- 238000005070 sampling Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 101100339395 Hydra viridissima Hln-1 gene Proteins 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5252—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT image
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Facsimile Scanning Arrangements (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
Abstract
空間周波数成分を抑制する周期的パターン抑制方法およ
び装置において、鮮鋭度を低下させないようにする。 【解決手段】 ウェーブレット変換部32において、グ
リッドピッチに対応する空間周波数以上のレスポンスが
略ゼロとなる特性を有するローパスフィルタを用いて2
次元ウェーブレット変換を施して、ウェーブレット変換
係数信号を得る。抑制手段36は、方向判定手段35に
よるグリッド方向の判定結果に基づいて、各信号のうち
のグリッド成分を含む信号(縦グリッド時にはHL1)
に対して、さらにグリッド方向について1次元ウェーブ
レット変換を繰り返し施した後、変換係数信号のうちの
低周波側の成分をゼロにし、このゼロにした信号および
残りの信号に対して1次元逆ウェーブレット変換を施
す。逆ウェーブレット変換部34において、グリッド成
分が抑制された信号を用いて元の画像を復元する。
Description
制処理方法および装置に関し、より詳細には、例えば、
静止グリッドを使用して撮影を行うことにより得られ
る、該静止グリッドに対応する縞模様(モアレ成分を含
む)が被写体像に重畳された画像から、前記縞模様を低
減・除去する方法および装置に関するものである。
とこの放射線エネルギの一部が蓄積され、その後可視光
等の励起光を照射すると蓄積されたエネルギに応じた光
量の輝尽発光光を放射する蓄積性蛍光体(輝尽性蛍光
体)を利用して、人体などの被写体の放射線画像を一旦
シート状の蓄積性蛍光体に撮影記録し、蓄積性蛍光体シ
ートをレーザ光などの励起光で走査して輝尽発光光を発
生させ、得られた輝尽発光光をフォトマルチプライヤな
どの読取手段により光電的に読み取って画像信号を得、
この画像信号に基づいて被写体の放射線画像を写真感光
材料などの記録材料やCRTなどに可視像として出力さ
せる放射線記録再生システムが提案されている(特開昭
55-12429号、同56-11395号、同55-163472 号、同56-164
645号、同55-116340 号など)。
の放射線画像を撮影記録する際に、被写体により散乱さ
れた放射線がシートに照射されないように4本/mm程度
の細かなピッチで放射線の透過しない鉛などと透過しや
すいアルミニウムや木材などとが交互に配置された静止
グリッドを被写体とシートとの間に配置して撮影を行う
ことがある。この静止グリッドを用いて撮影を行うと被
写体により散乱された放射線がシートに照射されにくく
なるため、被写体の放射線画像のコントラストを向上さ
せることができるが、被写体像とともに静止グリッドに
対応した細かな縞模様状のグリッド像が記録される。
37号などにおいて、静止グリッドの縞模様に対応する空
間周波数成分を除去するフィルタリング処理を施すこと
により、縞模様が低減した観察しやすい画像を得る方法
を提案している。この方法においては、例えば、静止グ
リッドの格子配列ピッチ(以下グリッドピッチという)
が4本/mmである場合には、4.0cycle/mm付近の空間周
波数帯域に縞模様が現れるため、この周波数帯域のレス
ポンスを除去あるいは低減するフィルタによりフィルタ
リング処理を施して、縞模様を除去するようにしてい
る。
願平10−164737号に記載の方法で用いているフィルタ
は、使用される静止グリッドのグリッドピッチに対応す
る空間周波数付近だけでなく、それ以上の高周波成分の
レスポンスをもゼロとするものであり、静止グリッドの
縞模様だけでなく、本来の画像に含まれる高周波成分を
も抑制・除去してしまうため、鮮鋭度が低下した画像と
なるという問題がある。
り、鮮鋭度を低下させずに、例えば、静止グリッド等に
起因して画像中に表われる縞模様等の周期的パターンを
目立たなくすることができる周期的パターン抑制処理方
法および装置を提供することを目的とするものである。
ーン抑制処理方法は、画像信号中の周期的パターンを形
成する空間周波数成分を低減する周期的パターン抑制処
理方法であって、実空間領域で表された前記画像信号を
周波数領域で取り扱うことのできる変換画像信号に変換
する第1のステップと、該変換画像信号のうち、周期的
パターンの配列方向近傍のみにおける少なくとも該周期
的パターンの周波数に対応する空間周波数の成分を含む
所望の周波数範囲の画像信号を低減する第2のステップ
とを有することを特徴とするものである。
に周期的に表われるパターンを広く意味するものであ
り、例えば、静止グリッドに起因して画像中に発生する
縞模様や、テレビにおいて撮像系の特性により画像に発
生するモアレ的な像など種々のパターンを示すものであ
る。
理方法は、静止グリッドを使用した撮影により得られた
画像信号中の、前記静止グリッドに起因する空間周波数
成分を低減する周期的パターン抑制処理方法であって、
実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する第1のステッ
プと、該変換画像信号のうち、静止グリッドの格子配列
方向近傍のみにおける少なくとも該静止グリッドの格子
配列周波数に対応する空間周波数の成分を含む所望の周
波数範囲の画像信号を低減する第2のステップとを有す
るものとすることもできる。
る変換画像信号に変換」するとは、画像信号に含まれる
所望の周波数成分について所望の処理ができる信号に変
換することを意味し、特に本願発明においては、実空間
領域で表された画像信号に含まれる前記静止グリッドに
対応する周波数成分を抑制する処理ができる信号に変換
することを意味する。この変換処理としては、例えば、
実空間領域で表現された画像信号をフーリエ変換して周
波数領域表現された画像信号(周波数スペクトル)にし
たり、あるいは実空間領域表現された画像信号を所定の
帯域にサブバンド分割(ウェーブレット変換やラプラシ
アンピラミッド展開などの場合)する多重解像度変換な
どがある。
分」とは、静止グリッドそのものの空間周波数成分だけ
でなく、ナイキスト周波数以下のサンプリング周期でサ
ンプリングしたり、あるいは縮小処理することにより、
静止グリッドに起因して発生するモアレ成分も含む意味
である。なお「格子配列周波数に対応する空間周波数の
成分」としたのはこのためであって、グリッド周波数と
同じ空間周波数成分だけでなく、これと関係する前記モ
アレ周波数と同じ空間周波数成分を含む意味である。
数に対応する空間周波数の成分を含む所望の周波数範囲
の画像信号を低減する」とは、グリッド周波数に対応す
る空間周波数の成分およびその近傍の空間周波数成分を
抑制する処理を施すことを意味し、従来例のように、グ
リッド周波数に対応する空間周波数成分以上の略全ての
高周波成分を抑制するものではなく、鮮鋭度に影響を与
える、前記グリッド周波数に対応する空間周波数成分以
上の高周波成分(例えばナイキスト周波数に対応する空
間周波数成分近傍)を可能な限り抑制しないことを意味
する。
おける」とは、従来例のように、グリッド方向に拘わら
ず抑制処理を施すのとは異なり、静止グリッドの格子配
列方向あるいはその近傍方向について抑制処理を施し、
それ以外の方向については抑制処理を施さないことを意
味する。なお、ここでいう「静止グリッド」は、後述す
るように、実際に使用する静止グリッドのみに限らず、
使用が予定される(グリッド方向がそれぞれ異なる)各
静止グリッドであってもよく、これらのぞれぞれについ
ての「格子配列方向近傍のみ」であればよい。換言する
と、通常の装置において、使用される静止グリッドのグ
リッド方向はある程度決まっており、例えば主走査方向
あるいは副走査方向であるので、互いに直交する主走査
方向および副走査方向において、グリッド周波数に対応
する成分を抑制する処理を施すものとしてもよい。より
判りやすく言えば、使用が予定される静止グリッドのグ
リッド方向の略全てについて、前記抑制処理を施しても
よいということである。例えば、主走査方向をv軸、副
走査方向をu軸とするフーリエ空間上において、v軸上
の近傍で(v軸と直交するu軸の正負方向に多少の幅を
持って)且つv軸方向のグリッド周波数を含む所望の空
間周波数範囲の成分と、u軸上の近傍で(u軸と直交す
るv軸の正負方向に多少の幅を持って)且つu軸方向の
グリッド周波数を含む所望の空間周波数範囲の成分のみ
を抑制するとよい。この場合、一方の静止グリッドの格
子長さ方向は、他方の静止グリッドの格子配列方向と同
じ方向となるので、結果的には、互いに直交する静止グ
リッドの格子配列方向および格子長さ方向のいずれにつ
いても抑制処理を施すことになる。
理方法においては、実際に使用する静止グリッドのみに
ついて前記抑制処理を行なうのが好ましい。ここで、実
際に使用する静止グリッドのみについて前記抑制処理を
行なうということは、実際に使用する静止グリッドの格
子配列方向あるいはその近傍方向のみについて抑制処理
を施し、それ以外の方向(その静止グリッドの格子長さ
方向も含む)については抑制処理を施さないことを意味
する。例えば、縦方向の静止グリッドを用いた場合に
は、上述の例において、v軸上の近傍で且つv軸方向の
グリッド周波数を含む所望の空間周波数範囲の成分のみ
を抑制することを意味する。
においては、前記第1のステップを、静止グリッドの空
間周波数以上のレスポンスが略ゼロとなるように帯域分
割するローパスフィルタを用いて前記画像信号に対して
2次元ウェーブレット変換を施すことにより複数の変換
画像信号を得るものとすると共に、前記第2のステップ
を、複数の変換画像信号のうちの、格子配列周波数に対
応する空間周波数の成分を含む信号について、さらに、
所定周波数以下の成分を低減する処理を施し、その後、
逆ウェーブレット変換を行なうものとするのが望まし
い。
波数に対応する空間周波数の成分を含む変換画像信号
(ウェーブレット変換係数信号)に対して、該静止グリ
ッドの格子長さ方向(グリッド方向)について、所定の
帯域分割フィルタを用いて1次元ウェーブレット変換を
所定回数だけ再帰的に繰り返し施し、該1次元ウェーブ
レット変換により得られた複数の画像信号(ウェーブレ
ット変換係数信号)のうちの、低周波成分の変換係数を
ゼロにし、その後、1次元逆ウェーブレット変換を施す
ことにより、前記所定周波数以下の成分を低減するもの
とすればより望ましい。
ずしも静止グリッドの空間周波数以上のレスポンスが略
ゼロとなるように帯域分割するローパスフィルタである
必要はない。
て1次元ウェーブレット変換を施すに際しては、例えば
縦方向グリッドのときには2次元ウェーブレット変換に
より得られたHL成分に対して副走査方向に、横方向グ
リッドのときにはLH成分に対して主走査方向に、それ
ぞれ1次元ウェーブレット変換を施すとよい。一方、ク
ロスグリッドのときにはHH成分に対して主走査方向の
1次元ウェーブレット変換を施した後に副走査方向の1
次元ウェーブレット変換を施す(主・副の順序が逆も
可)とよい。
レット変換により得られた画像信号のうちの、静止グリ
ッドの成分を含む画像信号についてハイパスフィルタを
用いて原画像を表す成分のみを抽出することで、静止グ
リッド成分を表す所定周波数以下の成分を低減すること
もできる。
においては、第2のステップを、前記複数の変換画像信
号(ウェーブレット変換係数信号)のパワーをそれぞれ
求め、求めた各パワーがそれぞれ所定の閾値よりも大き
いか否かに基づいて、静止グリッドの格子長さ方向(グ
リッド方向)を判定し、この判定結果に基づいて、前記
所定周波数以下の成分を低減する処理を施すものとする
ことが望ましい。
理方法においては、実際に使用する静止グリッドのみに
ついてに限らず、使用が予定された各静止グリッドにつ
いて、前記2次元ウェーブレット変換により得た変換画
像信号を用いて、前記第2のステップによる処理を施す
ものとしてもよい。
においては、前記第1のステップを、実空間領域で表さ
れた画像信号に対して、静止グリッドの格子長さ方向に
ついて、所定の帯域分割フィルタを用いて1次元ウェー
ブレット変換を施すことにより複数の変換画像信号を得
るものとすると共に、前記第2のステップを、複数の変
換画像信号のうちの、静止グリッドの格子配列周波数に
対応する空間周波数の成分を含む低周波側の変換画像信
号について、さらに、所定周波数以下の成分を低減する
処理を施し、その後、逆ウェーブレット変換を行なうも
のとすることもできる。
る静止グリッドのみについてに限らず、使用が予定され
た各静止グリッドについて、前記第1および第2の各ス
テップによる処理を施すものとすることができる。
は、上記方法を実施する装置、すなわち、画像信号中の
周期的パターンを形成する空間周波数成分を低減する周
期的パターン抑制処理装置であって、実空間領域で表さ
れた画像信号を周波数領域で取り扱うことのできる変換
画像信号に変換する画像信号変換手段と、該変換画像信
号のうち、周期的パターンの配列方向近傍のみにおける
少なくとも該周期的パターンの周波数に対応する空間周
波数の成分を含む所望の周波数範囲の画像信号を低減す
る周期的パターン成分抑制手段とを備えたことを特徴と
するものである。
理装置は、静止グリッドを使用した撮影により得られた
画像信号中の、静止グリッドに起因する空間周波数成分
を低減する周期的パターン抑制処理装置であって、実空
間領域で表された画像信号を周波数領域で取り扱うこと
のできる変換画像信号に変換する画像信号変換手段と、
該変換画像信号のうち、静止グリッドの格子配列方向近
傍のみにおける少なくとも該静止グリッドの格子配列周
波数に対応する空間周波数の成分を含む所望の周波数範
囲の画像信号を低減する静止グリッド成分抑制手段とを
備えたものとすることもできる。
においては、画像信号変換手段を、静止グリッドの空間
周波数以上のレスポンスが略ゼロとなるように帯域分割
するローパスフィルタを用いて前記画像信号に対して2
次元ウェーブレット変換を施すことにより複数の変換画
像信号を得るものとし、静止グリッド成分抑制手段を、
前記複数の変換画像信号のうちの、静止グリッドの格子
配列周波数に対応する空間周波数の成分を含む画像信号
について、さらに、所定周波数以下の成分を低減する処
理を施し、その後、逆ウェーブレット変換を行なうもの
とするのが望ましい。
段を、格子配列周波数に対応する空間周波数の成分を含
む変換画像信号に対して、静止グリッドの格子長さ方向
について、所定の帯域分割フィルタを用いて1次元ウェ
ーブレット変換を所定回数だけ再帰的に繰り返し施し、
該1次元ウェーブレット変換により得られた複数の画像
信号のうちの、低周波成分の変換係数をゼロにし、その
後、1次元逆ウェーブレット変換を施すことにより、前
記所定周波数以下の成分を低減するものするととより望
ましい。
理装置においては、ウェーブレット変換を施すことによ
り得られた複数の画像信号のパワーをそれぞれ求め、求
めた各パワーがそれぞれ所定の閾値よりも大きいか否か
に基づいて、静止グリッドの格子長さ方向(グリッド方
向)を判定する静止グリッド方向判定手段をさらに備え
たものとすると共に、静止グリッド成分抑制手段を、静
止グリッド方向判定手段による判定結果に基づいて、所
定周波数以下の成分を低減する処理を施すものとするの
が望ましい。
実際に使用する静止グリッドのみについてに限らず、使
用が予定された各静止グリッドについて、所定周波数以
下の成分を低減する処理を施すものとすることができ
る。
理装置においては、画像信号変換手段を、実空間領域で
表された画像信号に対して、静止グリッドの格子長さ方
向について、所定の帯域分割フィルタを用いて1次元ウ
ェーブレット変換を施すことにより複数の変換画像信号
を得るものとすると共に、静止グリッド成分抑制手段
を、前記複数の変換画像信号のうちの、静止グリッドの
格子配列周波数に対応する空間周波数の成分を含む低周
波側の画像信号について、さらに、所定周波数以下の成
分を低減する処理を施し、その後、逆ウェーブレット変
換を行なうものとすることもできる。
用が予定された各静止グリッドの格子長さ方向につい
て、1次元ウェーブレット変換を施すものとすると共
に、静止グリッド成分抑制手段についても、使用が予定
された各静止グリッドについて低減する処理および逆ウ
ェーブレット変換を行なうものとすることができる。
法および装置によれば、周期的パターンが表われた画像
の画像信号を周波数領域で取り扱うことのできる画像信
号に変換し、該周波数領域で取り扱うことのできる画像
信号のうち、周期的パターンの略配列方向における周期
的パターンの周波数と同じ空間周波数の成分を含む所望
の周波数範囲の画像信号を低減するようにしたので、画
像中に表われた周期的パターンを目立たなくすることが
できる。
因したものとした場合には、静止グリッドを使用した撮
影により得られた実空間領域で表された画像信号を周波
数領域で取り扱うことのできる画像信号に変換し、該周
波数領域で取り扱うことのできる画像信号のうち、静止
グリッドの略格子配列方向におけるグリッド周波数と同
じ空間周波数の成分を含む所望の周波数範囲の画像信号
を低減するようにしたので、静止グリッドに起因する縞
模様が目立たなくなる。一方、静止グリッドに起因する
空間周波数成分以上の全ての高周波成分を抑制すること
はないので、その分だけ、最終的に得られる画像信号に
はナイキスト周波数に対応する空間周波数成分近傍の成
分が含まれることとなり、鮮鋭度を低下させることがな
い。
ずフーリエ空間上の全方向に抑制処理を行なうのとは異
なり、静止グリッドの格子配列方向についてのみグリッ
ド成分が除去されるので、必要以上に画像を劣化させる
ことがない。
グリッドの格子配列方向のみについて前述の抑制処理を
施すようにすれば、一層必要以上に画像を劣化させるこ
とがなくなる。
いにグリッド方向の異なるもに)の各格子配列方向につ
いて、それぞれ前述の抑制処理を施すようにすれば、実
際に使用されるグリッド方向を気にすることなく(2次
元ウェーブレット変換を行なうものの場合には特にその
効果が大きい)、上記抑制処理の効果を得ることができ
る。
スポンスが略ゼロとなるように帯域分割するローパスフ
ィルタを用いて画像信号に対して2次元ウェーブレット
変換を施し、該ウェーブレット変換により得られた画像
信号のうちの、静止グリッドの成分を含む画像信号につ
いて、さらに、所定周波数以下の成分を低減する処理を
施すようにすれば、この所定周波数以下の成分が低減さ
れた信号を用いて元の画像(解像度レベル0)を復元し
たときには、静止グリッドの成分が低減された画像を復
元することができる。
り得られる低解像度信号成分には、静止グリッドの空間
周波数成分が含まれなくなるので、以後、さらにウェー
ブレット変換を多段階に亘って施して低解像度画像信号
を得ても、各ウェーブレット変換係数信号には静止グリ
ッドの空間周波数成分が含まれない。したがって、画像
の復元に際して、解像度レベル0まで復元することな
く、途中の解像度レベルまで復元しても、レベル1以降
の低解像度信号にはグリッド成分が含まれないので、縮
小画像(低解像度画像)において、静止グリッドに起因
するモアレが発生するということもない。つまり、画像
を任意に拡大あるいは縮小処理しても、モアレを発生す
る心配がない。
波数領域で取り扱うことのできる画像信号を得る場合に
おいて、静止グリッドの成分を含む画像信号(ウェーブ
レット変換係数信号)に対して、グリッド方向につい
て、所定の帯域分割フィルタを用いて1次元ウェーブレ
ット変換を所定回数だけ再帰的に繰り返し施し、該1次
元ウェーブレット変換により得られた複数の画像信号の
うちの、低周波成分の変換係数をゼロにするようにすれ
ば、繰り返しの回数が多いほど、復元された画像におい
ては、静止グリッドの成分のみが抑制された画像にする
ことができ、一層鮮鋭度の高い画像を得ることができ
る。
り得られた複数の画像信号のパワーをそれぞれ求め、求
めた各パワーがそれぞれ所定の閾値よりも大きいか否か
に基づいてグリッド方向を判定するようにすれば、静止
グリッド成分の抑制処理に際して、どのようなグリッド
方向のものを使用して撮影を行ったかを予め知る必要が
ない。
得られた実空間領域で表された画像信号に対して、静止
グリッドの格子長さ方向について、所定の帯域分割フィ
ルタを用いて1次元ウェーブレット変換を施し、該ウェ
ーブレット変換により得られた複数の画像信号のうち
の、静止グリッドの成分を含む低周波側の画像信号につ
いて、さらに、所定周波数以下の成分を低減する処理を
施すことによっても、静止グリッドに起因する空間周波
数成分を含む所望の周波数範囲の画像信号を低減するこ
とができ、この静止グリッドに起因する成分が低減され
た信号を用いて画像を復元したときには、復元の解像度
レベルに拘わらず、グリッドの成分が低減された画像を
復元することができる。
施の形態について詳細に説明する。なお、以下の説明で
は、上記特開昭55-12492号や特開昭56-11395号などに記
載されいているように、記録シートとしての蓄積性蛍光
体シートを利用した放射線画像情報記録再生システムに
おいて、前記シートに記録された人体の放射線画像をレ
ーザビーム走査によりデジタル画像信号として読み取る
態様を例に説明する。
た図である。図1に示すように、放射線源1から放射さ
れた放射線2は、被写体3を経由しさらにグリッド4に
到達する。グリッド4は、放射線2を吸収する鉛4a
と、放射線2を透過するアルミニウム4bとが、放射線
源1から発せられた放射線2がアルミニウム4bを経由
して、シート11に真っ直ぐに入射するように、位置に
応じて多少傾きをもって(図1参照)、4本/mmのピッ
チで交互に配置されているものである。このため放射線
源1から発せられ被写体3を真っ直ぐに透過した放射線
2は、鉛4aに吸収されて遮ぎられる一方、アルミニウ
ム4bを透過してシート11を照射し、シート11には
被写体像とともに4本/mmの縞模様状のグリッド像が蓄
積記録される。一方、被写体3内で散乱された散乱放射
線2aはグリッド4の前記傾きに対して斜めに入射する
ため、アルミニウム4b部に入射したものもグリッド4
内部で鉛4aに吸収され、またはグリッド4の表面で反
射されるため、シート11には照射されず、したがって
シート11には散乱放射線2aの照射の少ない鮮明な放
射線画像が蓄積記録される。
ことによりシート11に蓄積記録された、被写体像(図
の斜線部)5と縞模様状のグリッド像6を示した図であ
る。このようにシート11には被写体像5とグリッド像
6とが重畳された放射線画像が記録される。
処理方法を実施する周期的パターン抑制処理装置の一実
施の形態を適用した、放射線画像読取装置の斜視図であ
る。
線画像が記録されたシート11は、不図示の駆動手段に
より駆動されるエンドレスベルトなどのシート搬送手段
15により、走査ピッチ10本/mmで矢印Y方向に搬送
(副走査)される。一方、レーザー光源16から発せら
れた光ビーム17はモータ24により駆動され矢印方向
に高速回転する回転多面鏡18によって反射偏向され、
fθレンズなどの集束レンズ19を通過した後、ミラー
20により光路を変えてシート11に入射し副走査の方
向(矢印Y方向)と略直角な矢印X方向に主走査する。
シート11の光ビーム17が照射された箇所からは、蓄
積記録されている放射線画像情報に応じた光量の輝尽発
光光21が発散され、この輝尽発光光21は光ガイド2
2の入射端面22aに入射し、該ガイド22内の内部を
全反射を繰り返して進み、射出端面22bから射出して
フォトマル(光電子増倍管)23に受光されることによ
り、放射線画像を表す輝尽発光光21が該フォトマル2
3によって光電的に検出されて電気信号Saに変換され
る。なお、画像信号Sには、良好な放射線可視画像を再
生出力するために必要な、所望とする範囲の空間周波数
帯のうちの最高の空間周波数(後述するナイキスト周波
数)fn=5.0cycle/mmより低い、図2に示すグリッド
像6に関する情報である4.0cycle/mmの空間周波数帯域
の情報も含まれている。このグリッド像6に関する情報
は可視画像を観察する際にその可視画像を見にくくする
原因のひとつとなるものであり、取り除くべき情報であ
る。
対数的に増幅された後、A/D変換器28において空間
周波数fs=10.0 cycle/mmに対応するサンプリング間
隔でサンプリングされてディジタル化され、読取密度の
高い、高密度画像を担持するディジタルの画像信号Sが
得られる。
ート11に対して主走査方向(横方向)にレーザビーム
を走査させながらシート11を副走査方向(縦方向)に
移動させてシート11を2次元走査して得られた画像情
報を表すことになる。なお、このようにして得られた画
像信号Sには、ナイキスト周波数fn以下の情報を担持
しており、したがって図2に示すグリッド像6の情報
(4.0cycle/mm)も含まれている。なお、この実施の形
態においてはグリッド像6の情報(4.0cycle/mm)の2
倍以上の空間周波数fsに対応するサンプリング間隔で
デジタル化しているので、エリアジングによるグリッド
像6のモアレは発生しない。
れた後、画像信号処理部30に入力され、以下のように
して画像処理が施される。
理方法を実施するための画像信号処理部(周期的パター
ン抑制処理装置の一態様)30の概略を表すブロック図
である。図5に示すように本発明による周期的パターン
抑制処理方法を実施するための画像信号処理部30は、
記憶部29から読み出した画像信号Sに対して多重解像
度分解処理を施す画像信号変換手段としての多重解像度
分解処理手段32と、多重解像度分解処理手段32にお
いて複数の周波数帯域に分解された各帯域制限画像信号
(バンドパス信号,サブバンド信号)のうち、所定の周
波数帯域の画像(帯域制限画像)に対して所望の(例え
ば強調処理など)処理を施す処理手段33と、処理手段
33により所望の処理が施された帯域制限画像信号およ
び他の周波数帯域を担持する帯域制限画像信号を復元し
て処理済画像信号を得るための復元処理手段34と、復
元処理手段34により復元された処理済画像信号に基づ
いて、可視画像を再生するための画像出力手段35を有
する。
域で表される画像信号Sを周波数領域で取り扱うことの
できる画像信号に変換するに際して、ウェーブレット変
換処理を用いた多重解像度分解処理を使用するものとす
る。したがって、前記多重解像度分解処理手段32は、
放射線画像を表す画像信号Sに対してウェーブレット変
換を施すウェーブレット変換手段として機能し、前記復
元処理手段34は逆ウェーブレット変換処理を施して処
理済画像信号を得る逆ウェーブレット変換手段として機
能する。以下の説明においては、多重解像度分解処理手
段32をウェーブレット変換部32といい、復元処理手
段34を逆ウェーブレット変換部34という。
行われる多重解像度分解処理としてのウェーブレット変
換処理を説明するための機能ブロック図であり、図7は
各ウェーブレット変換手段32aの詳細を示したブロッ
ク図である。図6に示すように、ウェーブレット変換部
32には解像度(画素密度)レベルに応じた段数分のウ
ェーブレット変換手段32aが設けられ、逆ウェーブレ
ット変換部34には、ウェーブレット変換手段32aと
同じ段数分の逆ウェーブレット変換手段34aが設けら
れる。
0は、ウェーブレット変換により得た信号HL1,LH
1,HH1のパワーをそれぞれ求め、求めた各パワーが
それぞれ所定の閾値TH1よりも大きいか否かに基づい
て、静止グリッド6の格子長さ方向、すなわちグリッド
方向を判定する静止グリッド方向判定手段(以下方向判
定手段という)35と、信号HL,LH,HHのうち、
静止グリッド4の格子配列方向(グリッド方向と直交す
る方向)における、少なくとも該静止グリッド4に起因
する空間周波数成分を含む所望の周波数範囲の画像信号
を低減する静止グリッド成分抑制手段36とを備えてい
る。
理方法について説明する。なお、本実施の形態において
は、ウェーブレット変換の各係数が直交する、2次元の
直交ウェーブレット変換を行なうものとする。
から読み出された原画像を表すデジタルの画像信号Sが
ウェーブレット変換部32に入力されると、画像信号S
を原画像信号Sorg として、該原画像信号Sorg に対し
てウェーブレット変換が施される。すなわち、原画像信
号Sorg (信号LL0と等価)の主走査方向にウェーブ
レット関数H1,G1によりフィルタリング処理を行う
とともに、主走査方向の画素を1画素おきに間引き(図
中↓2で表す)、主走査方向の画素数を1/2にする。
ここで、関数H1はハイパスフィルタであり、関数G1
はローパスフィルタである。さらに、この画素が間引か
れた信号のそれぞれに対して副走査方向に前記関数H
1,G1によりフィルタリング処理を行うとともに、副
走査方向の画素を1画素おきに間引き、副走査方向の画
素数を1/2にして、ウェーブレット変換係数信号(以
下単に信号ともいう)HH1,HL1,LH1,LL1
を得る。ここで、信号LL1は原画像の縦横を各々1/
2に縮小した1/4縮小画像を表し、それぞれ原画像の
1/4縮小画像において、信号LH1は副走査方向(縦
方向)の高周波成分(横エッジ)を表す画像、信号HL
1は主走査方向(横方向)の高周波成分(縦エッジ)を
表す画像、信号HH1は対角方向の高周波成分(斜めエ
ッジ)を表す画像、信号LL1は原画像に対して1/2
解像度の低周波成分の画像を表すことになる。帯域分割
の基準となる周波数、すなわち、低周波成分と各高周波
成分との境目の周波数は、関数H1,G1のフィルタ特
性によって決まり、この初段での関数H1,G1による
フィルタリング処理の際には、ローパスフィルタとして
の関数G1は、静止グリッド4のグリッドピッチに対応
して、4.0cycle/mm以上の空間周波数の伝達特性(レス
ポンス)が略ゼロとなる特性を有するようにし、ハイパ
スフィルタとしての関数H1は、関数G1のローパス特
性を補うハイパス特性を有するようにする。このような
特性を有するローパスフィルタとしては、例えば特願平
10−164737号に記載のモアレ除去フィルタと同じ特性を
有するもの、すなわち静止グリッド4のグリッドピッチ
に対応する空間周波数成分に対して97%以上の空間周
波数成分のレスポンスを5%以下に低減するものを用い
るとよく、例えば、表1に示す(17,7)タップ、表
2に示す(13,7)タップ、表3に示す(15,5)
タップなどのフィルタ係数を用いたウェーブレット変換
フィルタで実現することができる。なお、各ウェーブレ
ット変換フィルタにおけるローパスフィルタの周波数応
答特性を図8に示す。
用されたときにはグリッド成分が信号LH1に現れ、縦
方向グリッドが使用されたときにはグリッド成分が信号
HL1に現れ、クロスグリッドが使用されたときにはグ
リッド成分が信号HH1に現れ、グリッド方向に拘わら
ず、信号LL1においては、静止グリッド4の空間周波
数成分が十分に抑制されているので、静止グリッド4に
起因する縞模様(モアレ)が殆ど現れない。
1,LH1,LL1のうち、グリッド成分を含んでいる
可能性のある各信号HH1,HL1,LH1が方向判定
手段35および静止グリッド成分抑制手段36に入力さ
れる。
おいて、基本ウェーブレット関数H0,G0を用いて、
信号LL1に対してウェーブレット変換が施されて、信
号HH2,HL2,LH2,LL2が得られる。ここ
で、信号LL2は原画像の縦横を各々1/4に縮小した
1/16縮小画像を表し、それぞれ原画像の1/16縮
小画像において、信号HL2、LH2およびHH2は縦
エッジ、横エッジおよび斜めエッジ成分の画像を表すも
のとなる。なお、上述のように、信号LL1にはグリッ
ド成分が殆ど現れないので、この2段目で用いられるウ
ェーブレット関数H0,G0は、上記初段で用いられた
ウェーブレット関数H1,G1とは異なり、静止グリッ
ド4のグリッドピッチに対応して設定されるものである
必要はなく、例えば、表4に示すドーベッチ(Daubechi
es)の(9,7)タップのフィルタ係数を用いたウェー
ブレット変換フィルタとするとよい。なお、この(9,
7)タップウェーブレット変換フィルタにおけるローパ
スフィルタの周波数応答特性を図8に示す。
得られるウェーブレット変換係数信号LLkに対するウ
ェーブレット変換をn回繰り返すことによりウェーブレ
ット変換係数信号HH1〜HHn,HL1〜HLn,L
H1〜LHn,LL1〜LLnを得る。ここで、n回目
のウェーブレット変換により得られるウェーブレット変
換係数信号HHn,HLn,LHn,LLnは、原画像
信号Sorg と比較して主副各方向の画素数が(1/2)n
となった(1/2)2n縮小画像を表し、各ウェーブレット
変換係数信号HHn,HLn,LHn,LLnはnが大
きいほど周波数帯域が低くなる。このように、ウェーブ
レット変換係数信号HHk,HLk,LHk,LLk
(kは解像度レベルであって、k=1〜nの整数)は、
それぞれ、原画像信号Sorg の周波数範囲のうちの所定
範囲の周波数成分を担持する帯域制限画像信号となる。
また、信号HHkは原画像信号Sorg の主副両方向の周
波数の変化を表しkが大きいほど低周波信号となる。ま
た信号HLkは原画像信号Sorg の主走査方向の周波数
の変化を表すものであり、kが大きいほど低周波信号と
なる。さらに信号LHkは原画像信号Sorg の副走査方
向の周波数の変化を表すものであり、kが大きいほど低
周波信号となる。
ブレット変換により得られた各信号LH1,HL1,H
H1のパワーをそれぞれ求め、求めた各パワーがそれぞ
れ所定の閾値TH0よりも大きいか否かを確認し、その
結果に基づいて、静止グリッド4のグリッド方向を判定
する。具体的には、上述のように、グリッド成分は、横
方向グリッド時には信号LH1に現れ、縦方向グリッド
時には信号HL1に現れ、クロスグリッド時には信号H
H1に現れるので、信号LHのパワーのみが閾値TH0
よりも大きいときには横方向グリッド、信号HLのパワ
ーのみが閾値TH0よりも大きいときには縦方向グリッ
ド、信号HLおよび信号LHのいずれもが閾値TH0よ
りも大きいときにはクロスグリッドが使用されているも
のと判定する。この判定結果は、静止グリッド成分抑制
手段36に入力される。
いて、方向判定手段35によるグリッド方向の判定結果
に基づいて、信号LH1,HL1,HH1のうちの静止
グリッドの成分を含む信号に対して、さらに静止グリッ
ド4の格子長さ方向(グリッド方向)について、所定の
帯域分割フィルタを用いて1次元ウェーブレット変換を
所定回数だけ再帰的に繰り返し施す。例えば、図2に示
すように縦方向グリッドが使用されている場合には、H
L1に対して縦(副走査)方向の1次元ウェーブレット
変換を施して低周波側を担持する信号HL1(L1)と
高周波側を担持する信号HL1(H1)を得、次に低周
波側を担持する信号HL1(L1)に対して縦方向の1
次元ウェーブレット変換を施して低周波側を担持する信
号HL1(L2)と高周波側を担持する信号HL1(H
2)を得、以下同様に、低周波側を担持する信号HL1
(Lj)(j=1〜mの整数)に対して縦方向の1次元
ウェーブレット変換をm回繰り返すことにより、信号H
L1に対して複数段(m段)の縦方向の1次元ウェーブ
レット変換を施して、ウェーブレット変換係数信号HL
1(H1)〜HL1(Hm),HL1(Lm)を得る。
して各成分に分解した結果を示す図である。なお、図9
においては、1回目の2次元ウェーブレット変換を行な
い(同図(a))、その後2回目の2次元ウェーブレッ
ト変換を行なった後に、信号HL1について、さらに縦
方向の1次元ウェーブレット変換を1回施した状態まで
(同図(b))を表している。
レット変換を繰り返し施すことによる効果を説明する模
式図である。
走査方向についての低周波成分を表し、図10に示すよ
うに、信号HL1は縦方向グリッドの成分と原画像信号
自身が持つ縦方向を主成分とする画像信号を含むことに
なる。そして、この信号HL1に対して格子長さ方向
(本例では副走査方向)について1次元ウェーブレット
変換を施すと、縦方向グリッドを用いた場合のグリッド
成分は1次元の縦長の低周波信号と見られるために、H
L1(Lj)には縦方向グリッドの成分を含む低周波成
分、HL1(Hj)にはその他の高周波成分が含まれる
ように、帯域分割される。この1次元ウェーブレット変
換を繰り返し施すと、十分に絞り込まれた低周波信号H
L1(Lm)は、グリッドピッチの大きさに拘わらず、
グリッドの目の方向(格子長さ方向)成分を含む極めて
低周波数の成分を表し、mを大きくするほど信号HL1
(Lm)がグリッド成分のみを表すようになる。
平・垂直あるいは水平・垂直が1対1のクロス方向には
置かれない場合が多いので、繰返しの回数mをあまり大
きくすると、信号HL1(Hj)にもグリッド成分が現
れるようになるので、1次元ウェーブレット変換は数回
程度繰り返すに留めるのが好まい。
の1次元ウェーブレット変換を施して得られた各信号H
L1(H1)〜HL1(Hm),HL1(Lm)のう
ち、最も低周波側の成分である信号HL1(Lm)をゼ
ロにする、つまり、グリッド成分を含む信号HL1(L
m)を抑制するようにする。その後、ゼロにした信号H
L1(Lm)および各信号HL1(H1)〜HL1(H
m)に対して1次元逆ウェーブレット変換を施して、信
号HL1’を得る。ここで、前述のように、少なくとも
静止グリッド4の空間周波数成分を含む信号HL1(L
m)が抑制されているので、グリッド成分を含む所定範
囲の空間周波数成分が低減された信号HL1’が得られ
る。
て、所望の(例えば強調処理など)処理が施され、その
後、逆ウェーブレット変換部34において、信号LL
n、HLk,LHkに対して、レベルnからレベル1ま
で、順次逆ウェーブレット変換が施される。
構成を示す概略ブロック図である。図11に示すよう
に、最低周波数帯域の信号HHn,HLn,LHn,L
Lnに対して逆ウェーブレット変換手段34aにおいて
逆ウェーブレット変換を施して信号LLn−1を得る。
aにおいて行われる逆ウェーブレット変換処理を説明す
るための機能ブロック図である。図12に示すように信
号LLnおよび信号LHn(LHk)の副走査方向に対
して画素間に1画素分の間隔をあける処理を行うととも
に(図中↑2で表す)、ウェーブレット変換の際に用い
た関数G0,H0に対応する逆ウェーブレット変換関数
G0’,H0’によりフィルタリング処理を副走査方向
に施してこれらを加算し、さらに加算により得られた信
号(第1の加算信号とする)の主走査方向に対して画素
間に1画素分の間隔をあける処理を行うとともに、関数
G0’によりフィルタリング処理を主走査方向に施して
第1の信号を得る。一方、信号HLn(HLk)および
信号HHn(HHk)の副走査方向に対して画素間に1
画素分の間隔をあける処理を行うとともに、関数G
0’,H0’によりフィルタリング処理を副走査方向に
施してこれらを加算し、さらに加算により得られた信号
(第2の加算信号とする)の主走査方向に対して画素間
に1画素分の間隔をあける処理を行うとともに、関数H
0’によりフィルタリング処理を主走査方向に施して第
2の信号を得る。そして第1および第2の信号を加算し
て信号LLn−1(LLk−1)を得る。
n−1,LLn−1に対して上記と同様に逆ウェーブレ
ット変換手段34aにおいて逆ウェーブレット変換を行
って、処理済み信号LLn−2を得る。そして、以下上
記と同様にして逆ウェーブレット変換を解像度レベル1
まで繰り返すことにより信号LL1が得られる。
は、信号LL1、LH1,HH1および信号HL1’に
対して、初段におけるウェーブレット変換の際に用いた
関数G1,H1に対応する逆ウェーブレット変換関数G
1’,H1’により逆ウェーブレット変換を施して画像
を再構成する。ここで、信号HL1’はグリッド成分が
低減されたものであり、再構成画像は、静止グリッド4
の成分が軽減されている。
あって、図13(a)は特願平10−164737号(以下従来
例という)に記載の方法を適用した場合における周波数
応答特性を示した図、図13(b)は本願発明を適用し
た場合における周波数応答特性を示した図、図13
(c)は従来例および本願発明による方法におけるフー
リエ空間での空間周波数特性を示した図である。なお、
図13(c)は、主走査方向をv軸、副走査方向をu軸
とするフーリエ空間上において、縦方向の静止グリッド
を用いた場合を示したものである。
れる静止グリッドのグリッドピッチに対応する空間周波
数付近だけでなく、それ以上の高周波成分のレスポンス
もゼロとするフィルタを用いているので、図13(a)
に示すように、グリッドピッチに対応する空間周波数以
上の高周波成分までもが抑制され、静止グリッド4の縞
模様だけでなく、本来の画像に含まれる高周波成分をも
除去・低減してしまうことになり、鮮鋭度が低下した画
像となってしまう。また、方向を考慮したフィルタリン
グ処理となっていないため、図13(c)に示すよう
に、フーリエ空間上における低周波領域(図中の中央部
の空白部)の範囲外全て(図中の斜線部)が抑制され
る、つまり縦方向グリッドによる縦模様(本当に抑制し
たい成分)だけでなく、本来の画像に含まれる横模様や
斜め模様などの、本来であれば抑制したくない高周波成
分も抑制されることになる。
は、初段のウェーブレット変換の際に用いられるフィル
タとして、特願平10−164737号に記載のものと同じ特性
を示すものを使用したとしても、図13(b)に示すよ
うに、静止グリッド4の空間周波数成分近傍の所定範囲
のみを抑制(カット)し、残りの高周波成分を保存する
ことができるので、静止グリッドに起因する縞模様が目
立たたず、しかも鮮鋭度の高い高画質な画像を得ること
ができる。また、信号HL1,LH1,HH1のうちの
グリッド成分を含む信号のみが抑制対象となるようにグ
リッド方向に応じて切り替えることができるので、図1
3(c)に示すように、フーリエ空間上における、グリ
ッド方向の軸上(縦グリッドのときはv軸)の近傍で
(該軸と直交する方向に多少の幅を持って)、且つ静止
グリッド4の空間周波数成分を含む高周波側の所定範囲
のみを抑制することができるので、例えば縦方向グリッ
ドの場合には、グリッドによる縦模様だけを抑制し、本
来の画像に含まれる横模様や斜め模様などの高周波成分
を抑制することがない。
により得られた変換係数信号HH1,HL1,LH1,
LL1のうち、グリッド成分を含んでいる可能性のある
各信号HH1,HL1,LH1に基づいて静止グリッド
の方向を判定することができるので、静止グリッド成分
の抑制処理に際して、どのようなグリッド方向のものを
使用して撮影を行ったかを予め知る必要もない。
像度レベル1まで復元すれば、前述のように、信号LL
1においては静止グリッドの空間周波数成分が十分に抑
制されているので、静止グリッド4に起因する縞模様が
殆ど現れない。また、グリッド成分が抑制された信号L
L1に対してウェーブレット変換を施して得たレベル1
以降のLLkにもモアレ成分が含まれないので、全ての
縮小画像において、静止グリッドに起因するモアレが発
生するということがなく、画像を任意に拡大あるいは縮
小処理しても、モアレを発生する心配がない。したがっ
て、強調処理などの所望の画像処理を施して、画像処理
が施された画像を逆多重解像度変換することにより復元
しても、いずれの解像度レベルの画像においても、復元
された処理済画像にはモアレによるアーチファクトが含
まれず、見やすい高画質の画像を提供することができ
る。
ーブレット変換の初段においてグリッド成分が抑制さ
れ、それ以降の解像度レベルの画像にはモアレが現れな
くなり、またウェーブレット変換を用いているので、拡
大/縮小画像を得るのにも好都合であり、ネットワーク
を介して画像信号を転送する上でも都合がよい。
て説明したが、本発明は上述した実施の形態に限定され
るものではなく、本発明の要旨を逸脱しない限りにおい
て、種々の変更を行うことが可能である。
6における作用説明は、静止グリッド4が縦方向グリッ
ドであるとして説明したが、グリッド方向が縦方向以外
の場合には、取り扱う信号HL1を、そのグリッド方向
に応じて信号LH1あるは信号HH1とすればよい。す
なわち、静止グリッド4が横方向グリッドである場合に
は、信号LH1に対してさらに複数段の横方向(主走査
方向)の1次元ウェーブレット変換を施し、低周波側の
信号成分をゼロにした後、1次元逆ウェーブレット変換
して画像を再構成するとよく、静止グリッド4がクロス
グリッドである場合には、信号HLおよび信号LHそれ
ぞれに対してさらに複数段の縦方向(信号HL)/横方
向(信号LH)の1次元ウェーブレット変換を施し、低
周波側の信号成分をゼロにした後、1次元逆ウェーブレ
ット変換して画像を再構成するとよい。
対する1次元ウェーブレット変換の代わりに、グリッド
成分を低減させるようなハイパスフィルタを用いてもよ
い。
何に関わらず、使用が予定された静止グリッドの各グリ
ッド方向について上述の処理を施す、例えば図9に示し
た処理に加えて、信号LH1に対して主走査方向につい
て1次元ウェーブレット変換を繰り返し施すようにすれ
ば、図13(c)に示すv軸の近傍の静止グリッド4の
空間周波数成分を含む高周波側の所定範囲だけでなく、
u軸の近傍の静止グリッド4の空間周波数成分を含む高
周波側の所定範囲についても抑制効果を得ることができ
る。したがって、通常は主走査方向若しくは副走査方向
のいずれかと同方向のグリッドが用いられることが多い
ので、この場合、グリッド成分は信号HL1(縦グリッ
ドのとき),LH1(横グリッドのとき)のいずれかに
含まれることになり、結果的に、グリッド方向を全く気
にすることなく、上記抑制処理の効果を得ることができ
る。また、この場合においても、無制限な方向に対して
抑制処理を施すものではなく、図13(C)に示すよう
に、抑制処理の対象方向が使用が予定された静止グリッ
ドの各グリッド方向についてのみに限られるので、画質
に与える影響を従来よりも軽減することができる。な
お、信号HH1に対して主/副の両走査方向について1
次元ウェーブレット変換を繰り返し施すことにより、ク
ロスグリッドの使用が予定される場合にも対応できる。
次元ウェーブレット変換を施すものとして説明したが、
グリッド方向が予め判っている場合には、図14に示す
ように、その方向に応じた(縦グリッドのときには副走
査方向に)1次元ウェーブレット変換を原画像信号に対
して施し、グリッド成分を含むウェーブレット変換係数
信号(図14(B)の信号L1)に対してグリッド成分
を低減する処理を施すようにしてもよい。なお、1次元
ウェーブレット変換の際に用いられるウェーブレット関
数は、静止グリッド4のグリッドピッチに拘わらず、任
意のものであってもよく、上記H0,G0を用いること
ができる。
目の方向(格子長さ方向)を表す極めて周波数の成分が
ウェーブレット変換係数信号のうちの低周波側の成分
(信号L1)に現れるので、この低周波側のウェーブレ
ット変換係数信号L1について、さらにグリッド成分を
低減する処理、例えば同方向の1次元ウェーブレット変
換を繰り返し施した後(図14(C)参照)に最も低周
波のウェーブレット変換係数信号をゼロにしたり、ある
いはハイパスフィルタでグリッド成分を除く原画像を表
す信号成分のみを抽出することで、グリッド成分を表す
極めて周波数の成分を抑制することができる。
た信号を用いて1次元逆ウェーブレット変換を施して画
像を復元したときには、復元の解像度レベルに拘わら
ず、グリッドの成分が低減された画像を復元することが
できる。
際に使用されるグリッド方向の如何に関わらず、使用が
予定された静止グリッドの各グリッド方向について上述
の処理をそれぞれ施すものとすることができる。例えば
図14に示した副走査方向の処理に加えて、主走査方向
の処理を行なうものとし、最終的に得られた2つの画像
のうち、グリッド成分の少ない方の画像を使用するよう
にすることもできる。
像信号Sを周波数領域で取り扱うことのできる画像信号
に変換するに際して、ウェーブレット変換処理を用いた
多重解像度分解処理を使用するものとして説明したが、
画像信号に含まれる静止グリッド成分を含む所望範囲の
周波数成分を低減する処理を施すことができるものであ
る限り、どのような変換方法を用いてもよい。例えばラ
プラシアンピラミッド展開を用いた多重解像度分解処理
を使用することもできる。また、実空間領域表現された
画像信号をフーリエ変換して周波数領域表現された信号
(周波数スペクトル)にしてもよい。
図
態様としての画像信号処理部の概略を表すブロック図
ブレット変換処理を説明するための機能ブロック図
ク図
フィルタの周波数応答特性を示す図
分解した結果を示す図
とによる効果を説明する模式図
ロック図
逆ウェーブレット変換処理を説明するための機能ブロッ
ク図
態様) 32 ウェーブレット変換部(画像信号変換手段の一
態様) 34 逆ウェーブレット変換部 35 静止グリッド方向判定手段 36 静止グリッド成分抑制手段
Claims (16)
- 【請求項1】 画像信号中の周期的パターンを形成する
空間周波数成分を低減する周期的パターン抑制処理方法
であって、 実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する第1のステッ
プと、 該変換画像信号のうち、前記周期的パターンの配列方向
近傍のみにおける少なくとも該周期的パターンの周波数
に対応する空間周波数の成分を含む所望の周波数範囲の
画像信号を低減する第2のステップとを有することを特
徴とする周期的パターン抑制処理方法。 - 【請求項2】 静止グリッドを使用した撮影により得ら
れた画像信号中の、前記静止グリッドに起因する空間周
波数成分を低減する周期的パターン抑制処理方法であっ
て、 実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する第1のステッ
プと、 該変換画像信号のうち、前記静止グリッドの格子配列方
向近傍のみにおける少なくとも該静止グリッドの格子配
列周波数に対応する空間周波数の成分を含む所望の周波
数範囲の画像信号を低減する第2のステップとを有する
ことを特徴とする周期的パターン抑制処理方法。 - 【請求項3】 前記第1のステップが、前記静止グリッ
ドの空間周波数以上のレスポンスが略ゼロとなるように
帯域分割するローパスフィルタを用いて前記画像信号に
対して2次元ウェーブレット変換を施すことにより複数
の前記変換画像信号を得るものであり、 前記第2のステップが、前記複数の変換画像信号のうち
の、前記格子配列周波数に対応する空間周波数の成分を
含む信号について、さらに、所定周波数以下の成分を低
減する処理を施し、その後、逆ウェーブレット変換を行
なうものであることを特徴とする請求項2記載の周期的
パターン抑制処理方法。 - 【請求項4】 前記第2のステップが、前記格子配列周
波数に対応する空間周波数の成分を含む変換画像信号に
対して、該静止グリッドの格子長さ方向について、所定
の帯域分割フィルタを用いて1次元ウェーブレット変換
を所定回数だけ再帰的に繰り返し施し、該1次元ウェー
ブレット変換により得られた複数の画像信号のうちの、
低周波成分の変換係数をゼロにし、その後、1次元逆ウ
ェーブレット変換を施すことにより、前記所定周波数以
下の成分を低減するものであることを特徴とする請求項
3記載の周期的パターン抑制処理方法。 - 【請求項5】 前記第2のステップが、前記複数の変換
画像信号のパワーをそれぞれ求め、求めた各パワーがそ
れぞれ所定の閾値よりも大きいか否かに基づいて、前記
静止グリッドの格子長さ方向を判定し、この判定結果に
基づいて、前記所定周波数以下の成分を低減する処理を
施すものであることを特徴とする請求項3または4記載
の周期的パターン抑制処理方法。 - 【請求項6】 使用が予定された各静止グリッドについ
て、前記第2のステップによる処理を施すことを特徴と
する請求項3または4記載の周期的パターン抑制処理方
法。 - 【請求項7】 前記第1のステップが、前記実空間領域
で表された画像信号に対して、前記静止グリッドの格子
長さ方向について、所定の帯域分割フィルタを用いて1
次元ウェーブレット変換を施すことにより複数の前記変
換画像信号を得るものであり、 前記第2のステップが、前記複数の変換画像信号のうち
の、前記静止グリッドの格子配列周波数に対応する空間
周波数の成分を含む低周波側の変換画像信号について、
さらに、所定周波数以下の成分を低減する処理を施し、
その後、逆ウェーブレット変換を行なうものであること
を特徴とする請求項2記載の周期的パターン抑制処理方
法。 - 【請求項8】 使用が予定された各静止グリッドについ
て、前記第1および第2の各ステップによる処理を施す
ことを特徴とする請求項7記載の周期的パターン抑制処
理方法。 - 【請求項9】 画像信号中の周期的パターンを形成する
空間周波数成分を低減する周期的パターン抑制処理装置
であって、 実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する画像信号変換
手段と、 該変換画像信号のうち、前記周期的パターンの配列方向
近傍のみにおける少なくとも該周期的パターンの周波数
に対応する空間周波数の成分を含む所望の周波数範囲の
画像信号を低減する周期的パターン成分抑制手段とを備
えたことを特徴とする周期的パターン抑制処理装置。 - 【請求項10】 静止グリッドを使用した撮影により得
られた画像信号中の、前記静止グリッドに起因する空間
周波数成分を低減する周期的パターン抑制処理装置であ
って、 実空間領域で表された前記画像信号を周波数領域で取り
扱うことのできる変換画像信号に変換する画像信号変換
手段と、 該変換画像信号のうち、前記静止グリッドの格子配列方
向近傍のみにおける少なくとも該静止グリッドの格子配
列周波数に対応する空間周波数の成分を含む所望の周波
数範囲の画像信号を低減する静止グリッド成分抑制手段
とを備えたことを特徴とする周期的パターン抑制処理装
置。 - 【請求項11】 前記画像信号変換手段が、前記静止グ
リッドの空間周波数以上のレスポンスが略ゼロとなるよ
うに帯域分割するローパスフィルタを用いて前記画像信
号に対して2次元ウェーブレット変換を施すことにより
複数の前記変換画像信号を得るものであり、 前記静止グリッド成分抑制手段が、前記複数の変換画像
信号のうちの、前記静止グリッドの格子配列周波数に対
応する空間周波数の成分を含む画像信号について、さら
に、所定周波数以下の成分を低減する処理を施し、その
後、逆ウェーブレット変換を行なうものであることを特
徴とする請求項10記載の周期的パターン抑制処理装
置。 - 【請求項12】 前記静止グリッド成分抑制手段が、前
記格子配列周波数に対応する空間周波数の成分を含む変
換画像信号に対して、該静止グリッドの格子長さ方向に
ついて、所定の帯域分割フィルタを用いて1次元ウェー
ブレット変換を所定回数だけ再帰的に繰り返し施し、該
1次元ウェーブレット変換により得られた複数の画像信
号のうちの、低周波成分の変換係数をゼロにし、その
後、1次元逆ウェーブレット変換を施すことにより、前
記所定周波数以下の成分を低減するものであることを特
徴とする請求項11記載の周期的パターン抑制処理装
置。 - 【請求項13】 前記複数の変換画像信号のパワーをそ
れぞれ求め、求めた各パワーがそれぞれ所定の閾値より
も大きいか否かに基づいて、前記静止グリッドの格子長
さ方向を判定する静止グリッド方向判定手段をさらに備
え、 前記静止グリッド成分抑制手段が、該静止グリッド方向
判定手段による判定結果に基づいて、前記所定周波数以
下の成分を低減する処理を施すものであることを特徴と
する請求項11または12記載の周期的パターン抑制処
理装置。 - 【請求項14】 前記静止グリッド成分抑制手段が、使
用が予定された各静止グリッドについて前記所定周波数
以下の成分を低減する処理を施すものであることを特徴
とする請求項11または12記載の周期的パターン抑制
処理装置。 - 【請求項15】 前記画像信号変換手段が、前記実空間
領域で表された前記画像信号に対して、前記静止グリッ
ドの格子長さ方向について、所定の帯域分割フィルタを
用いて1次元ウェーブレット変換を施すことにより複数
の前記変換画像信号を得るものであり、 前記静止グリッド成分抑制手段が、前記複数の変換画像
信号のうちの、前記静止グリッドの格子配列周波数に対
応する空間周波数の成分を含む低周波側の画像信号につ
いて、さらに、所定周波数以下の成分を低減する処理を
施し、その後、逆ウェーブレット変換を行なうものであ
ることを特徴とする請求項10記載の周期的パターン抑
制処理装置。 - 【請求項16】 前記画像信号変換手段が、使用が予定
された各静止グリッドの格子長さ方向についいて、前記
1次元ウェーブレット変換を施すものであり、 前記静止グリッド成分抑制手段が、前記使用が予定され
た各静止グリッドにつて前記低減する処理および前記逆
ウェーブレット変換を行なうものであることを特徴とす
る請求項15記載の周期的パターン抑制処理装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000395577A JP4265869B2 (ja) | 2000-01-20 | 2000-12-26 | 周期的パターン抑制処理方法および装置 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000-11174 | 2000-01-20 | ||
JP2000011174 | 2000-01-20 | ||
JP2000395577A JP4265869B2 (ja) | 2000-01-20 | 2000-12-26 | 周期的パターン抑制処理方法および装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2001273490A true JP2001273490A (ja) | 2001-10-05 |
JP2001273490A5 JP2001273490A5 (ja) | 2005-11-04 |
JP4265869B2 JP4265869B2 (ja) | 2009-05-20 |
Family
ID=26583819
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000395577A Expired - Fee Related JP4265869B2 (ja) | 2000-01-20 | 2000-12-26 | 周期的パターン抑制処理方法および装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4265869B2 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006289095A (ja) * | 2005-04-14 | 2006-10-26 | Agfa Gevaert Nv | 画像の周期的パターンの抑制方法 |
JP2006318103A (ja) * | 2005-05-11 | 2006-11-24 | Fuji Photo Film Co Ltd | 画像処理装置および方法並びにプログラム |
JP2011010829A (ja) * | 2009-07-01 | 2011-01-20 | Toshiba Corp | X線画像撮影装置及び画像処理装置 |
JP2014174803A (ja) * | 2013-03-11 | 2014-09-22 | Omron Corp | パターン認識のためのモアレ除去方法およびこの方法を用いたモアレ除去装置ならびにプログラム |
-
2000
- 2000-12-26 JP JP2000395577A patent/JP4265869B2/ja not_active Expired - Fee Related
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006289095A (ja) * | 2005-04-14 | 2006-10-26 | Agfa Gevaert Nv | 画像の周期的パターンの抑制方法 |
JP2006318103A (ja) * | 2005-05-11 | 2006-11-24 | Fuji Photo Film Co Ltd | 画像処理装置および方法並びにプログラム |
US8107764B2 (en) | 2005-05-11 | 2012-01-31 | Fujifilm Corporation | Image processing apparatus, image processing method, and image processing program |
JP2011010829A (ja) * | 2009-07-01 | 2011-01-20 | Toshiba Corp | X線画像撮影装置及び画像処理装置 |
JP2014174803A (ja) * | 2013-03-11 | 2014-09-22 | Omron Corp | パターン認識のためのモアレ除去方法およびこの方法を用いたモアレ除去装置ならびにプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP4265869B2 (ja) | 2009-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1120743B1 (en) | Method and unit for suppressing a periodic pattern | |
JP5698577B2 (ja) | 画像処理装置、画像処理方法並びに画像処理プログラム | |
JPH0696200A (ja) | ノイズ減少方法および装置 | |
JP2002216126A (ja) | 画像処理装置、及びその方法、並びにプログラム、記憶媒体 | |
JPH06301766A (ja) | ピラミッド的画像分解に基づいた放射線画像の多重処理法 | |
JP2003150954A (ja) | 周期的パターン抑制処理方法および装置 | |
JP2004242285A (ja) | ノイズ抑制処理方法および装置並びにプログラム | |
JP4244094B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
EP0610603B1 (en) | Fast interactive off-line processing method for radiographic images | |
JPH10268451A (ja) | 細長対象物の放射画像の記録及び読取方法 | |
JP3706201B2 (ja) | 画像処理方法 | |
JP2000224421A (ja) | 画像処理方法および装置並びに記録媒体 | |
JP2004261514A (ja) | 画像処理装置、画像処理方法およびプログラム | |
JP2003008885A (ja) | 画像再生システム | |
JP2754068B2 (ja) | 放射線画像信号生成方法および放射線画像読取装置 | |
JP4209061B2 (ja) | 画像処理符号復号化方法および画像処理符号復号化システム、画像処理符号化装置および画像処理復号化装置、並びに記録媒体 | |
JPH0312785A (ja) | 画像処理方法 | |
JP4265869B2 (ja) | 周期的パターン抑制処理方法および装置 | |
US6614044B2 (en) | Image signal generating method, apparatus and program | |
JP3370797B2 (ja) | 画像重ね合せ方法およびエネルギーサブトラクション方法 | |
JPH0944659A (ja) | エネルギーサブトラクション処理方法および装置 | |
JPH06274614A (ja) | 画像処理方法 | |
JP2001189866A (ja) | 静止グリッド抑制処理方法および装置 | |
JP4163401B2 (ja) | 画像信号生成方法および装置並びにプログラム | |
JP2001223901A (ja) | 画像保存方法および装置並びに記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050912 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050912 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20061202 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20081112 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081118 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090113 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090210 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090212 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4265869 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120227 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120227 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130227 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140227 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |