JP2001086368A - 画像処理装置 - Google Patents
画像処理装置Info
- Publication number
- JP2001086368A JP2001086368A JP25618199A JP25618199A JP2001086368A JP 2001086368 A JP2001086368 A JP 2001086368A JP 25618199 A JP25618199 A JP 25618199A JP 25618199 A JP25618199 A JP 25618199A JP 2001086368 A JP2001086368 A JP 2001086368A
- Authority
- JP
- Japan
- Prior art keywords
- image
- processing
- pixel
- periodicity
- periodicity index
- 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
Landscapes
- Image Processing (AREA)
- Picture Signal Circuits (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
理技術を提供する。 【解決手段】 周期画像領域抽出部10によって、画像
から周期的な画像パターンの繰り返しが存在する周期画
像領域が抽出される。また、このときに周期性指標算出
部20が画像パターンの繰り返しの有無の程度を示す周
期性指標を算出する。算出された周期性指標からは、加
工後周期性指標算出部25によって加工後周期性指標が
求められる。一方、差分画像作成部30は元の画像から
差分画像を作成する。その差分画像に基づいてコントラ
スト量算出部50がコントラスト量を算出する。その
後、ぼかし量決定部60およびシフト量決定部70が周
期性指標、加工後周期性指標およびコントラスト量に基
づいてそれぞれぼかし量およびシフト量を決定する。そ
れらぼかし量およびシフト量に従ってぼかし処理部95
およびシフト処理部90が画像のぼかし処理およびシフ
ト処理を行う。
Description
を施す画像処理装置に係り、特には、モアレの発生の抑
制などに効果がある画像処理装置に関する。
どの画像入力器で原稿または対象物体を光学的に読み取
り、デジタル化することによってデジタル画像を取得
し、そのデジタル画像に対して種々の処理を施すことに
より行われる。
は取得したデジタル画像にシャープネス処理が施され
る。また、シャープネス処理以外にも、画像中の画素の
階調値を変換する階調変換処理等が行われる。
例えば、デジタル画像に対してシャープネス処理などの
画像処理を行った場合、処理後の画像にモアレが発生す
ることがあった。
や発生原因を調査した。その結果、まず、このようなモ
アレは、画像内に布地や網模様、格子パターンなどの周
期性を有する画像パターンが存在する画像領域で発生す
ることが判明した。
発生すると考えられる。例えば、原稿内に図30に示す
ような周期性を有する画像パターン(白黒の縞模様)GP
Tが存在し、この原稿を画像入力器で読み取ってデジタ
ル画像を得る場合を考える。このとき、原稿の画像パタ
ーンGPTと画像入力器の光学読み取り素子(光電素子)
IDとの位置関係によって、得られたデジタル画像の各
画素の階調値は変化する。
素子IDとの位置関係が図30(a)に示す状態であれ
ば、得られたデジタル画像の各画素の階調値は、図31
(a)の階調プロファイルに示すように、白黒が明確に
分かれる「INフェイズ状態」となる。また、原稿の画
像パターンGPTと光電素子IDとの位置関係が図30
(b)に示す状態であれば、得られたデジタル画像の各
画素の階調値は、図31(b)の階調プロファイルに示
すように、白黒の中間値が発生する「OUTフェイズ状
態」となる。
フェイズ状態」と「OUTフェイズ状態」とで相違する
が、各々の平均階調値AKIN、AKOUTは略等しい。
シャープネス処理を施すと、「INフェイズ状態」で読
み取られた画像と「OUTフェイズ状態」で読み取られ
た画像とで平均階調値AKIN、AKOUTに違いが生じ
る。例えば、図31(a)、(b)に示す各画像にそれ
ぞれシャープネス処理を施すと、計算上、各画素の階調
値は各々図32(a)、(b)に示すようになる。ここ
で、計算上、階調値が階調レンジKR(図では、0〜2
55)から外れるサチレーションが起きると、その階調
値は、強制的に階調レンジKRの上下限値に修正され
る。すなわち、図31(a)、(b)に示す各画像にそ
れぞれシャープネス処理を施すと、処理後の各画素の階
調値は各々図33(a)、(b)の実線に示すようにな
る。
デジタル画像上で4画素周期の周期性を有しているが、
1周期ごとに、図33(a)では、白2画素だけにサチ
レーションSRが起き、図33(b)では、白1画素と
黒1画素とにサチレーションSRが起きている。その結
果、図33(a)(「INファイズ状態」)の平均階調値
AKINは、図33(b)(「OUTファイズ状態」)の平
均階調値AKOUTに比べて低くなる。
Tの白黒の幅が、光電素子IDの幅の整数倍である場合
を示したが、画像パターンGPTの白黒の幅が、光電素
子IDの幅の整数倍でない場合、原稿の画像パターンG
PTと光電素子IDとの位置関係が周期的にずれるの
で、得られた画像上では、画像パターンGPTは、「I
Nフェイズ状態」と「OUTフェイズ状態」とで読み込
まれる状態が周期的に繰り返されることになる。この場
合、このデジタル画像にシャープネス処理を施すと、平
均階調値が周期的に変動する状態が生じる。この平均階
調値が周期的に変動する状態が、モアレとなって現れる
ものと考えられる。
画像であるため、画像の晶質が著しく低下することにな
る。
ス処理を施した場合について説明したが、例えば、図3
4に示すようなルックアップテーブルによって各画素の
階調値を変換するような階調変換処理など、サチレーシ
ョンに偏りが生じるような画像処理を施した場合にも上
記のような平均階調値が周期的に変動する状態が発生
し、その結果、処理後の画像にモアレが発生すると考え
られる。
クアップテーブルによって各画素の階調値を変換する階
調変換処理などでも、「INフェイズ状態」と「OUT
フェイズ状態」とで平均階調値に違いが生じ、上記のよ
うな平均階調値が周期的に変動する状態が発生する結
果、処理後の画像にモアレが発生すると考えられる。
であり、モアレの発生の抑制などに効果がある画像処理
装置を提供することを目的とする。
め、請求項1の発明は、画像に対して処理を施す画像処
理装置であって、画像から周期的な画像パターンの繰り
返しが存在する周期画像領域を抽出する周期画像領域抽
出手段と、前記周期画像領域における注目画素につい
て、当該注目画素を含む注目領域と当該注目領域以外の
周辺領域との相関特性に基づいて画像パターンの繰り返
しの有無の程度を示す周期性指標を算出する周期性指標
算出手段と、前記周期性指標によって表現された周期性
指標画像を加工することによって加工後周期性指標画像
を取得し、当該加工後周期性指標画像における前記注目
画素の加工後周期性指標を算出する加工後周期性指標算
出手段と、前記画像における注目画素とその周辺画素と
の明暗比を示す指標であるコントラスト量を算出するコ
ントラスト量算出手段と、前記加工後周期性指標および
前記コントラスト量に基づいて、前記画像における注目
画素についてのぼかし処理の程度であるぼかし量を決定
するぼかし量決定手段と、前記周期性指標、前記加工後
周期性指標および前記コントラスト量に基づいて、前記
画像における注目画素についてのシフト処理の程度であ
るシフト量を決定するシフト量決定手段と、前記ぼかし
量決定手段により決定されたぼかし量に基づいて、前記
画像における注目画素についてのぼかし処理を行うぼか
し処理手段と、前記シフト量決定手段により決定された
シフト量に基づいて、前記画像における注目画素につい
てのシフト処理を行うシフト処理手段と、を備え、前記
加工後周期性指標算出手段に、前記周期性指標画像内を
複数の領域に分割し、前記複数の領域のそれぞれに含ま
れる画素の階調値を平均化して1つの画素の階調値とし
て出力することにより、前記周期性指標画像をバイリニ
ア縮小する縮小手段と、前記バイリニア縮小によって得
られた縮小画像に対してメディアンフィルタによるフィ
ルタ操作を施すフィルタ手段と、前記フィルタ操作後の
縮小画像に含まれる画素の階調値を補間しつつ、当該縮
小画像を前記周期性指標画像と同等の大きさに拡大する
ことにより、前記加工後周期性指標画像を得る拡大手段
と、を含ませている。
理を施す画像処理装置であって、画像から周期的な画像
パターンの繰り返しが存在する周期画像領域を抽出する
周期画像領域抽出手段と、前記周期画像領域を前記画像
パターンの繰り返しの1周期分拡張して拡張周期画像領
域を取得し、前記拡張周期画像領域における注目画素に
ついて、画像パターンの繰り返しの有無の程度を示す周
期性指標を算出する周期性指標算出手段と、前記画像に
おける注目画素とその周辺画素との明暗比を示す指標で
あるコントラスト量を算出するコントラスト量算出手段
と、前記周期性指標および前記コントラスト量に基づい
て、前記画像における注目画素についてのぼかし処理の
程度であるぼかし量を決定するぼかし量決定手段と、前
記周期性指標および前記コントラスト量に基づいて、前
記画像における注目画素についてのシフト処理の程度で
あるシフト量を決定するシフト量決定手段と、前記ぼか
し量決定手段により決定されたぼかし量に基づいて、前
記画像における注目画素についてのぼかし処理を行うぼ
かし処理手段と、前記シフト量決定手段により決定され
たシフト量に基づいて、前記画像における注目画素につ
いてのシフト処理を行うシフト処理手段と、を備えてい
る。
実施の形態について詳細に説明する。
画像処理装置のハードウェア構成の一例を示すブロック
図である。この画像処理装置はいわゆるコンピューター
システムを用いて構成されている。CPU1は、内部メ
モリ2に記憶された処理プログラムに従って、モアレの
発生を抑制する画像処理(後述する周期画像領域の抽出
処理やコントラスト量算出を含む一連の処理であり、以
下「モアレ抑制処理」とも略す)の他、シャープネス処
理などの通常の画像処理も実行する。CPU1と内部メ
モリ2とはバスラインBLを介して接続されている。内
部メモリ2には、上述した処理プログラムを記憶するプ
ログラム記憶部2a以外にも、少なくともモアレ抑制処
理を施す前のデジタル画像(以下、「処理前画像」とも
言う)を記憶する処理前画像記憶部2bやモアレ抑制処
理を施した後のデジタル画像(以下、「処理後画像」と
も言う)を記憶する処理後画像記憶部2cを有する。
スIFを介して、記録媒体用ドライバ3、外部記憶装置
4、表示装置5、入力装置6、画像入力器7などにも接
続されている。
理をCPU1に実行させるための処理プログラムが記録
されている記録媒体8が装填される。記録媒体8から読
み出された処理プログラムは内部メモリ2内のプログラ
ム記憶部2aに記憶され、CPU1に実行される。
成される外部記憶装置4は、処理前画像や処理後画像な
どの画像の保存などに用いられる。また、CRTなどで
構成される表示装置5は、処理前画像や処理後画像、シ
ャープネス処理後の画像などの表示などに用いられる。
さらに、キーボードやマウスなどで構成される入力装置
6は、オペレータからの各種の指示や処理に必要なデー
タの設定などに用いられる。
7によって読み込まれたデジタル画像(処理前画像)は、
内部メモリ2内の処理前画像記憶部2bに記憶され、こ
の処理前画像に対して、後述するような周期画像領域の
抽出処理やコントラスト量算出等の各処理が施される。
を示す機能ブロック図である。同図に示すように、この
画像処理装置は機能的に、周期画像領域抽出部10と、
周期性指標算出部20と、加工後周期性指標算出部25
と、差分画像作成部30と、コントラスト量算出部50
と、ぼかし量決定部60と、シフト量決定部70と、シ
フト量生成部80と、シフト処理部90と、ぼかし処理
部95とを備えている。さらに、差分画像作成部30
は、加工部40を含んでいる。これらの処理部はいずれ
も、CPU1内において、上記処理プログラムによって
実現される手段であり、それぞれの機能および処理内容
については後に詳述する。
装置における画像処理の処理手順の概要を示すフローチ
ャートである。画像処理に先立って、原稿または対象物
体が画像入力器7によってデジタル画像として読み込ま
れ、処理前画像記憶部2bに記憶される。そして、読み
込まれたデジタル画像は、周期画像領域抽出部10に与
えられ、画像から周期的な画像パターンの繰り返しが存
在する周期画像領域が抽出される(ステップS1)。ま
た、このときに後述する相関特性に基づいて、周期性指
標算出部20が画像パターンの繰り返しの有無の程度を
示す周期性指標を算出する(ステップS2)。算出され
た周期性指標はシフト量決定部70に伝達されるととも
に、加工後周期性指標算出部25にも付与され、加工後
周期性指標算出部25によって加工後周期性指標が算出
される(ステップS3)。
画像作成部30にも与えられ、処理前画像から差分画像
が作成される。差分画像は加工部40によって加工さ
れ、その加工後の差分画像に基づいてコントラスト量算
出部50が明暗比を示す指標であるコントラスト量を算
出する(ステップS4)。なお、ステップS1〜S3の
処理とステップS4の処理とはその処理順序が逆であっ
ても良い。
後周期性指標およびコントラスト量に基づいてぼかし処
理の程度であるぼかし量を決定する(ステップS5)。
また、シフト量決定部70が上記の周期性指標、加工後
周期性指標およびコントラスト量に基づいてシフト処理
の程度であるシフト量を決定する(ステップS6)。な
お、ステップS5の処理とステップS6の処理とはその
処理順序が逆であっても良い。
し処理部95が読み込まれた画像のぼかし処理を行い
(ステップS7)、決定されたシフト量に従ってシフト
処理部90がシフト処理を行う(ステップS8)。
行われる。なお、ステップS1からステップS7の処理
は画像中の画素単位で行われるものであり、原則として
読み込まれた画像の全画素について行われる。また、処
理が終了した画像は処理後画像記憶部2cに記憶され
る。
細に説明する。
性指標の算出>周期性を有する画像パターンが存在する
画像領域、すなわちモアレの発生が予測される周期的な
画像パターンの繰り返しが存在する周期画像領域の抽出
は、周期画像領域抽出部10が処理前画像の自己相関特
性を調べることにより行われる。
を注目画素(xc,yc)として、以下の数1によって注目画
素(xc,yc)の周辺領域の自己相関データS(a,b)を求め
る。
P(x,y)は処理前画像の画素(x,y)の階調値、m,nは差
分積算領域Eを決める定数、a,bは自己相関を比較す
るズラシ量、wx,wyは1つの中心画素(xc,yc)に対し
て自己相関特性を調べる範囲Wを決める定数である。
×3)、wx=xy=2(a=b=-2〜+2)とした場合において、a=b=
+2のときの自己相関データS(a,b)の算出形態を図4に
示す。
固定値として処理するようにしてもよいし、入力装置6
からオペレータによって適宜に変更可能に構成してもよ
い。
タを2次元的に求めているので処理量が多くなる。そこ
で、例えば、2次元の処理前画像内で互いに直交する2
つの画素列方向であるx方向、y方向それぞれに沿った
周期的な画像パターンの繰り返しの存在を調べて周期画
像領域を抽出するために、x方向、y方向それぞれに沿
った自己相関データH(a)、V(b)を以下の数2、数3に
よって求めて、処理の高速化を図るようにしてもよい。
1)、wx=2(a=-2〜+2)とした場合において、a=+2のとき
のx方向に沿った自己相関データH(a)の算出形態を図
5に、また、(xc,yc)=(4,4)、n=1(差分積算領域:1×
3)、wy=2(b=-2〜+2)とした場合において、b=+2のとき
のy方向に沿った自己相関データV(b)の算出形態を図
6にそれぞれ示す。
b)、または、H(a)、V(b))に基づき、画像中における
周期的な画像パターンの繰り返しの有無を調べる。
しが存在すると、その画像パターンの周期ごとに自己相
関が高くなり、上記数1、数2、数3で求まる自己相関
データは規則的に小さくなる。従って、まず、(A)自
己相関データの極小値を検索し、(B)それら極小値が
所定レベル以下で、かつ、(C)それら極小値が規則的
に存在していることを調べる。
5)とした場合のP(x,y)、H(a)の一例を示すデータとそ
のH(a)をグラフ化した図である。すなわち、同図は注
目画素(xc,yc)=(7,3)について、当該注目画素を含む注
目領域と当該注目領域以外の周辺領域との相関特性を示
すものである。なお、a=0は、同じ画素同志の自己相関
であるので、H(0)=0となり極小値になる。
(A)の処理は、〔(H(k-1)>H(k))and(H(k)<H(k+
1))〕の条件を満たすkを+側と−側とで求める。この条
件を満たすkについてのH(k)の値が極小値となる。
満たすH(k)、すなわち極小値となるH(k)が所定のしき
い値以下となるか否かで判定する。このしきい値は、予
め入力装置6等によって複数が設定されており、例えば
図7においては、”SL1=7.5”、”SL2=5”の2つ
が設定されている。そして、極小値となるH(k)が少な
くとも最も高い値のしきい値SL1以下となっていれ
ば、上記(B)の処理において、所定レベル以下である
と判定される。
ど、H(a)の全体レベルが上がるため、しきい値は固定
値ではなく、H(a)の最大値に対する比率(例えば、S
L1=58%,SL2=38%)で規定してもよい。
(A)の条件を満たす+側のkをkp、−側のkをkmとした
とき(ABS{kp+km}≦1)を満たすか否かで規則性の有無
を判定する。また、例えば、上記(A)の条件を満たす
各極小値に対して、隣接する各極小値間の幅の相互の差
分の絶対値が各々1以下であるか否か、すなわち、ABS
{D(d+1)−D(d)}≦1(d=0,1,2,・・・:図7では、d=O)
によって規則性の有無を判定することもできる。後者の
判定によれば、例えば、極小値が+側、−側にそれぞれ
2個以上存在する場合にも適用できる。
きかったり(最も高い値のしきい値SL1を越えていた
り)、それら極小値が不規則に存在しているような場合
は周期性が有るとは言い難いが、上記(B)によりレベ
ル判定が行え、上記(C)により極小値の規則性が判別
でき、周期性の有無を確実に判定できる。
件を全て満たす場合、注目画素(xc,yc)の周囲の範囲W
内の画像に周期的な画像パターンの繰り返しが存在する
ことになる。このことは換言すれば、注目画素(xc,yc)
の周囲の範囲W内の画像が周期画像領域として抽出され
たことを意味する。例えば、図7に示す自己相関データ
H(a)には周期性が有るため、注目画素(xc,yc)の周囲の
範囲W内の画像が周期画像領域として抽出される。これ
に対して、図8に示すような自己相関データH(a)には
周期性が無く、周期画像領域として抽出されることもな
い。
Thは((kp一km)/2))または(D(d))である。
無や周期値(Th)について説明したが、y方向に沿った周
期性の有無や周期値(Tv)も同様の処理により求めること
ができる。
y方向の周期値(Tv)が判れば、周期方向θは(arctan(Tv
/Th))であり、処理前画像上の実際の周期値(周期方向
に沿った周期値)Tは(Th×sinθ)(または、Tv×cosθ)
により求めることができる。
は、上記と同様の処理を2次元的に行うことで、周期性
の有無や処理前画像上の実際の周期方向及び周期値を得
ることができるが、上述したようにx方向、y方向それ
ぞれに沿った自己相関データH(a)、V(b)を用いても同
様の結果を高速に得ることができる。なお、自己相関デ
ータS(a,b)を用いた場合、実際の画像パターンの周期
を種々の方向から調べることになるが、このとき、周期
値が最小となる方向がその画像パターンの周期方向とな
る。
(A)、(B)、(C)の処理を行う過程において、図
7に示したような相関特性に基づいて、周期性指標算出
部20が画像パターンの繰り返しの有無の程度を示す周
期性指標Mを算出する。自己相関データH(a)は、周期
性の度合い(画像パターンの繰り返しの有無の程度)を
示すものであるため、極小値となるH(k)がより低い値
のしきい値以下となっている場合には、周期性の度合い
が大きいと言える。例えば、図7において、極小値H(-
3)はしきい値SL1としきい値SL2との間であるのに
対して、極小値H(+3)はしきい値SL2よりも小さい。
つまり、極小値H(+3)の方が周期性の度合いが大きいこ
とを示している。従って、極小値H(k)がいずれのしき
い値以下となっているかを判定することによって、画像
パターンの繰り返しの有無の程度を判別することができ
るのである。
算出部20が周期性指標Mを算出する。周期性指標Mの
算出に関与するのは、上記(A)の処理において極小値
と判定され、上記(B)の処理において最高しきい値S
L1以下と判定され、かつ上記(C)の処理において規
則性有りと判定された極小値H(k)のみである。但し、
H(0)は、注目画素自身との相関であって、必ず0にな
る値であるため、周期性指標Mの算出には関与しない。
従って、例えば、図7においては、しきい値SL1以下
の極小値H(-3)および極小値H(+3)がいずれのしきい値
以下となっているかによって周期性指標Mは算出され
る。なお、上記(A)、(B)、(C)の条件を満たす
H(k)が存在しない場合には、その注目画素(xc,yc)につ
いての周期性指標M(xc,yc)=0とすることにより、周
期性は存在しないものされる。
Mの値が0から1の範囲にて規定されており、例えば図
7の例では、しきい値SL1に”M=0”、しきい値S
L2に”M=1”が規定されている。そして、しきい値
SL2以下の極小値H(k)については”M(k)=1”とす
る。一方、しきい値SL1としきい値SL2との間の極
小値H(k)については、しきい値SL1としきい値SL
2との間において周期性指標Mが線形的に変化するもの
とみなして、算出する。例えば、図7の場合、以下の数
4に従って算出する。
されたM(k)の平均値が注目画素(xc,yc)についての周期
性指標M(xc,yc)となる。例えば、図7においては、極
小値H(+3)はしきい値SL2よりも小さいため、M(+3)
=1となる。また、極小値H(-3)はしきい値SL1とし
きい値SL2との間であるため、数4に従い、M(-3)=
(7.5-6)/(7.5-5)=0.6となる。よって、注目画素(xc,yc)
=(7,3)についての周期性指標M(7,3)=(1+0.6)/2=0.8
となる。
い値は2つに限定されるものではなく、3つ以上であっ
ても良い。最高しきい値と最低しきい値との間において
周期性指標Mが線形的に変化する場合は、2つで十分で
あるが、非線形的に変化する場合は3つ以上のしきい値
を設定した方が好ましい。
データH(a)の最大値について判定するためのしきい値
を設けるようにしても良い。自己相関データH(a)の最
大値が所定のしきい値よりも小さい場合は、均一な画像
中にノイズが存在している場合が多いと考えられ、この
ような場合は周期的な画像パターンの繰り返しであると
は言えない。従って、自己相関データH(a)の最大値が
所定のしきい値よりも小さい場合は、注目画素(xc,yc)
についての周期性指標M(xc,yc)=0とすることによ
り、周期性は存在しないものと判定するのである。
よび周期性指標算出部20による一連の処理を、処理前
画像内の全ての画素(または、一定間隔ごとの画素でも
よい)を順次注目画素(xc,yc)として行うことにより、処
理前画像内において周期性を有する画像パターンが存在
する画像領域、すなわち周期的な画像パターンの繰り返
しが存在する周期画像領域を全て抽出することができる
とともに、注目画素(xc,yc)についての周期性の度合い
(画像パターンの繰り返しの有無の程度)である周期性
指標M(xc,yc)を算出することができる。
と周期方向と周期値とを周期性データとする。
性指標算出部20によって算出された周期性指標Mは、
加工後周期性指標算出部25に伝達され、加工後周期性
指標算出部25によって加工後周期性指標が算出され
る。図10は、加工後周期性指標算出の手順を示すフロ
ーチャートである。
素について周期性指標M(xc,yc)が求められるため、周
期性指標M(xc,yc)によって表現された周期性指標画像
を作成することができる。周期性指標画像は、その画素
の値を周期性指標M(xc,yc)とする一種のグレー画像で
あり、加工後周期性指標算出部25によって作成され、
処理前画像記憶部2bに記憶される。そして、加工後周
期性指標算出部25は、その周期性指標画像に対して以
下の加工を行った後、その画素の値たる周期性指標を検
出して加工後周期性指標とするのである。
て、加工後周期性指標算出部25がバイリニア縮小を行
う(ステップS31)。バイリニア縮小とは、例えば、
周期性指標画像を1つの領域が4×4の画素群で構成さ
れる複数の領域に分割し、複数の領域のそれぞれに含ま
れる16の画素の階調値を平均して1つの画素の階調値
として出力する画像縮小処理である。このようなバイリ
ニア縮小によって、周期性指標画像は縦横ともに1/4
に縮小されるとともに、周期性指標画像の階調値が平滑
化されることとなる。
小画像に対して加工後周期性指標算出部25がメディア
ンフィルタによるフィルタ操作を行う(ステップS3
2)。メディアンフィルタは、その構成画素(例えば、
3×3=9画素)の階調値のうちの中央値(メディア
ン)を中央画素の出力値とするフィルタであり、画像中
のノイズ部分を除去するためのフィルタである。メディ
アンフィルタを上記縮小画像に対して順次に走査させる
と、当該縮小画像中のノイズ部分が除去される。周期性
指標画像を縮小した縮小画像中においてノイズ部分とな
るのは、元の周期性指標画像中において周期性指標が急
激に変化している局所領域である。
ンフィルタによるフィルタ操作を行うことは、元の周期
性指標画像に対して12×12のメディアンフィルタに
よるフィルタ操作を行うことと実質的に等価である。フ
ィルタ操作は、フィルタのサイズが小さいほどCPU1
の負担が軽く、高速に実行することができる。すなわ
ち、バイリニア縮小を行った後に小さなフィルタによる
フィルタ操作を行うことで、大きなフィルタによるフィ
ルタ操作を行ったのと同様の効果を高い処理効率にて得
ることができるのである。
加工部40がバイリニア拡大を行い、加工後の周期性指
標画像を得る(ステップS33)。バイリニア拡大と
は、フィルタ操作後の縮小画像に含まれる画素の階調値
を補間しつつ、当該縮小画像を元の周期性指標画像と同
等の大きさに拡大する処理である。
バイリニア拡大を概念的に説明するための図である。フ
ィルタ操作後の縮小画像において、図11(a)に示す
ように、周期性指標M=1.0、M=0.0、M=1.
0の順に並ぶ画素列が存在していたとする。これを縦横
ともに4倍に拡大すると、図11(b)に示すような元
の周期性指標画像と同等の大きさの画像が得られるので
あり、このときに周期性指標M=1.0の画素と周期性
指標M=0.0の画素との間に生成された画素列L1お
よび画素列L2については周期性指標M=1.0と周期
性指標M=0.0との間にて直線近似にて補間(バイリ
ニア補間)した周期性指標Mを付与する。このようにす
れば、フィルタ操作後の縮小画像において隣接する画素
の周期性指標Mの値が大きく異なっていたとしても、拡
大後の加工後周期性指標画像においては周期性指標Mの
値が滑らかに変化する状態(いわゆる「なじむ」状態)
となる。
部25による加工が実行されるのであり、バイリニア縮
小とバイリニア拡大との組み合わせによって周期性指標
画像に対する一種のぼかし処理を行うとともに、メディ
アンフィルタによるフィルタ操作を行うことによって、
周期性指標が急激に変化している局所領域が十分に除去
された加工後の周期性指標画像を取得することができ
る。
性指標算出部25が上記の加工後周期性指標画像から注
目画素の値(すなわち周期性指標)を検出することによ
って、当該注目画素(xc,yc)についての加工後周期性指
標PM(xc,yc)を算出する。ここでの注目画素とは、既
述した周期性指標の算出時の注目画素(xc,yc)と同じで
ある。
コントラスト量の算出について説明する。既述したよう
に、処理前画像は差分画像作成部30にも与えられ、差
分画像作成部30、加工部40およびコントラスト量算
出部50による処理を経て注目画素(xc,yc)についての
コントラスト量C(xc,yc)が算出される。注目画素(xc,y
c)についてのコントラスト量とは、注目画素(xc,yc)と
その周辺画素との明暗比を示す指標である。
すフローチャートである。図12に示した各処理のう
ち、ステップS41〜ステップS43の処理は差分画像
作成部30によって、ステップS44〜ステップS46
の処理は加工部40によって、ステップS47の処理は
コントラスト量算出部50によって実行されるものであ
り、以下これらについて順に説明する。
張画像および高濃度領域縮小画像を作成する(ステップ
S41,S42)。図13は、高濃度領域の拡張または
縮小について説明する図である。同図に示すように、白
地中に黒地の四角形の図形(斜線を付した部分)が描か
れた処理前画像の高濃度領域を拡張または縮小する場合
を例として説明する。
255”の高濃度領域であり、黒の図形部分が階調値P
=”0”の低濃度領域である。高濃度領域を拡張すると
きには、この画像に対して3×3の画素行列で構成され
るフィルタF1によるフィルタ操作を行う。フィルタF
1は、その構成画素(3×3=9画素)の階調値のうち
の最大値を中央画素の出力値とするフィルタである。そ
して、このフィルタF1を、画像に対して順次に走査さ
せることによって高濃度領域の拡張を行う。すなわち、
フィルタF1が黒の図形の周辺部分を走査するときに
は、高濃度(高階調値)である白地部分の階調値(P=
255)が出力され、当該周縁部分が黒から白地に置換
される。このことは、高濃度領域(白地)が拡張される
ことを意味しており、いわゆる太らせ処理が行われてい
ることとなる。
張された高濃度領域拡張画像を示す図である。高濃度領
域(白地)が拡張され、黒の図形部分が処理前画像より
も小さくなっている。
と同様に、処理前画像に対して3×3の画素行列で構成
されるフィルタF2によるフィルタ操作を行う。フィル
タF2は、その構成画素(3×3=9画素)の階調値の
うちの最小値を中央画素の出力値とするフィルタであ
る。従って、上記拡張と同様に、フィルタF2を処理前
画像に対して順次に走査させると、黒の図形の周辺部分
において、低濃度(低階調値)である黒の階調値(P=
0)が出力され、当該周辺部分が白地から黒に置換され
る。このことは、高濃度領域(白地)が縮小されること
を意味しており、いわゆる細らせ処理が行われているこ
ととなる。
小された高濃度領域縮小画像を示す図である。高濃度領
域(白地)が縮小され、黒の図形部分が処理前画像より
も大きくなっている。
張画像と高濃度領域縮小画像との差分である差分画像を
作成する(ステップS43)。本実施形態における差分
画像とは、高濃度領域拡張画像の画素の階調値から当該
画素に対応する高濃度領域縮小画像の画素の階調値を減
算して得られる画像である。すなわち、以下の数5によ
って得られる画像である。
素(x,y)の階調値、Pmax(x,y)は高濃度領域拡張画像の
画素(x,y)の階調値、Pmin(x,y)は高濃度領域縮小画像
の画素(x,y)の階調値である。
ある。図13に示した処理前画像と比較すると明らかな
ように、処理前画像においてコントラストが高い領域、
すなわち、ある画素とその周辺との明暗比が大きい領域
が差分画像においては大きな階調値を有している。この
ことは、差分画像とは、処理前画像からコントラストが
高い領域を抽出した画像であることを意味しており、差
分画像の作成とは高コントラスト領域の抽出処理である
と言える。
コントラスト領域にぼかし処理を施すことによってモア
レの発生を抑制している。ある程度面積の大きな高コン
トラスト領域(例えば、エアコンのフロントパネルにつ
いての画像領域)であれば、ぼかし処理が有効なのであ
るが、高コントラストであっても面積の小さな領域(例
えば、エアコンのスイッチについての画像領域)にぼか
し処理を施すと鮮鋭度を喪失するため適当ではない。そ
こで、本実施形態においては、加工部40によって上記
差分画像を加工し、高コントラストであっても面積の小
さな領域やノイズ部分を除去している(コントラストを
低下させている)。
〜ステップS46の3段階の工程によって行われる。こ
の加工は、上述した加工後周期性指標の算出時における
加工処理と類似の処理である。すなわち、まず、ステッ
プS43にて作成された差分画像について、加工部40
がバイリニア縮小を行う(ステップS44)。ここでの
バイリニア縮小とは、例えば、差分画像を1つの領域が
8×8の画素群で構成される複数の領域に分割し、複数
の領域のそれぞれに含まれる64の画素の階調値を平均
して1つの画素の階調値として出力する画像縮小処理で
ある。このようなバイリニア縮小によって、差分画像は
縦横ともに1/8に縮小されるとともに、差分画像の階
調値が平滑化されることとなる。
小画像に対して加工部40がメディアンフィルタによる
フィルタ操作を行う(ステップS45)。メディアンフ
ィルタは、既述した如く、その構成画素(例えば、3×
3=9画素)の階調値のうちの中央値(メディアン)を
中央画素の出力値とするフィルタであり、画像中のノイ
ズ部分を除去するためのフィルタである。メディアンフ
ィルタを上記縮小画像に対して順次に走査させると、当
該縮小画像中のノイズ部分が除去される。ここで縮小画
像中においてノイズ部分となるのは、元の処理前画像中
における面積の小さな高コントラスト領域やノイズ領域
である。
ンフィルタによるフィルタ操作を行うことは、元の差分
画像に対して24×24のメディアンフィルタによるフ
ィルタ操作を行うことと実質的に等価である。フィルタ
操作は、フィルタのサイズが小さいほどCPU1の負担
が軽く、高速に実行することができる。すなわち、バイ
リニア縮小を行った後に小さなフィルタによるフィルタ
操作を行うことで、大きなフィルタによるフィルタ操作
を行ったのと同様の効果を高い処理効率にて得ることが
できるのである。
加工部40がバイリニア拡大を行い、加工後の差分画像
を得る(ステップS46)。バイリニア拡大とは、フィ
ルタ操作後の縮小画像に含まれる画素の階調値を補間し
つつ、当該縮小画像を元の差分画像と同等の大きさに拡
大する処理である。
イリニア拡大を概念的に説明するための図である。フィ
ルタ操作後の縮小画像において、図17(a)に示すよ
うに、階調値P=255、P=0、P=255の順に並
ぶ画素列が存在していたとする。これを縦に8倍に拡大
すると、図17(b)に示すような元の差分画像と同等
の大きさの画像が得られるのであり、このときに階調値
P=255の画素と階調値P=0の画素との間に生成さ
れた画素列L3および画素列L4については階調値P=
255と階調値P=0との間にて直線近似にて補間(バ
イリニア補間)した階調値を付与する。このようにすれ
ば、フィルタ操作後の縮小画像において隣接する画素の
階調値が大きく異なっていたとしても、拡大後の差分画
像においては階調値が滑らかに変化する状態(いわゆる
「なじむ」状態)となる。このことは、後述するコント
ラスト量の急激な変化を抑制して、ぼかし処理の程度の
急激な変化による画質低下の防止につながるのである。
が実行されるのであり、バイリニア縮小とバイリニア拡
大との組み合わせによって差分画像に対する一種のぼか
し処理を行うとともに、メディアンフィルタによるフィ
ルタ操作を行うことによって、高コントラストであって
も面積の小さな領域やノイズ部分が十分に除去された加
工後の差分画像を取得することができる。
加工処理は必須の処理ではなく、処理前画像中に高コン
トラストであっても面積の小さな領域等がほとんど存在
しないような場合は、省略可能な処理である。
ト量算出部50が差分画像(加工後の差分画像または加
工を行っていない差分画像)から注目画素の階調値を検
出することによって、当該注目画素についてのコントラ
スト量を算出する。ここでの注目画素とは、既述した周
期性指標の算出時の注目画素(xc,yc)と同じである。
が高い領域を抽出した画像であり、その階調値が大きい
ということは、処理前画像においてコントラストが高い
ことを意味している。従って、差分画像における注目画
素(xc,yc)の階調値は、そのまま注目画素(xc,yc)とその
周辺画素との明暗比を示す指標となり、すなわち注目画
素(xc,yc)についてのコントラスト量C(xc,yc)となるの
である。
素についてのコントラスト量を、差分画像から算出する
ように説明したが、例えば次のようにコントラスト量を
算出してもよい。
右、上下方向の隣接画素)の階調値の差分を求め、求め
た差分値の平均値又は合計値に基づいてコントラスト量
を算出することができる。
えば、左右、上下方向の隣接画素)間の階調値の差分を
求め、求めた差分値の平均値又は合計値に基づいてコン
トラスト量を算出することができる。
して、注目画素(xc,yc)についての周期性指標M(xc,y
c)、加工後周期性指標PM(xc,yc)およびコントラスト
量C(xc,yc)が算出された後、ぼかし量決定部60が加
工後周期性指標PM(xc,yc)およびコントラスト量C(x
c,yc)に基づいてぼかし処理の程度であるぼかし量を決
定する(図3のステップS5)。本実施形態において、
ぼかし処理とは、平滑化処理のことであり、注目画素(x
c,yc)を含む所定領域に含まれる画素の階調値の加重平
均を注目画素(xc,yc)の階調値として出力する加重平均
処理である。
するが、加重平均フィルタ(いわゆる平滑化フィルタ)
を用いたフィルタ操作により行っている。そして、本実
施形態では、平滑化処理に使用する加重平均フィルタを
予め複数用意しておき、それらから加重平均フィルタを
選択することによってぼかし処理の程度であるぼかし量
を決定している。
ための図である。本実施形態においては、注目画素(xc,
yc)についての加工後周期性指標PM(xc,yc)とコントラ
スト量C(xc,yc)との積に基づいて加重平均フィルタを
選択している。すなわち、加工後周期性指標PMとコン
トラスト量Cとの積をぼかし量を規定するためのぼかし
量規定値(C×PM)とする。そして、ぼかし量規定値
(C×PM)に基づき、図18に示すようなルックアッ
プテーブル(LUT)を用いて加重平均フィルタのID
番号を決定する。同図に示すように、ID番号として6
4段階(ID=0〜63)が設定されており、ぼかし量
規定値(C×PM)の値が大きくなるほど、決定される
加重平均フィルタのID番号も大きくなる。
である。図19(a)に示す最も大きなID番号の加重
平均フィルタによれば、フィルタ内の全画素の階調値が
均等に平均化され、平滑化の程度が大きい(ぼかし量が
大きい)といえる。一方、図19(b)に示す最も小さ
なID番号の加重平均フィルタによれば、フィルタ内の
中心画素の階調値がそのまま出力されることとなり、こ
れは全く平滑化が行われない(ぼかし量が0)ことを意
味する。これら以外のID番号の加重平均フィルタは、
そのID番号の大きさに応じて上記2つの加重平均フィ
ルタによる平滑化の程度の間の平滑化を行うフィルタで
あり、換言すれば、ID番号の大きな加重平均フィルタ
ほど平滑化の程度が大きくなる。
が大きくなるほど、ルックアップテーブルによって決定
される加重平均フィルタのID番号が大きくなり、選択
される加重平均フィルタの平滑化の程度が大きく、すな
わちぼかし量が大きくなるのである。
法に限定されるものではなく、ぼかし量規定値(C×P
M)の値が大きくなるほど加重平均フィルタの平滑化の
程度が大きくなるような手法であれば良い。
処理の程度であるシフト量の決定について説明する(図
3のステップS6)。シフト量の決定は、シフト量生成
部80が生成したシフト量を基準にして、シフト量決定
部70が注目画素(xc,yc)についての周期性指標M(xc,y
c)、加工後周期性指標PM(xc,yc)およびコントラスト
量C(xc,yc)に基づいて決定する。
ト量について説明する。処理前画像の各画素(x,y)のシ
フト量は、図20に示すように、x方向のシフト量dx
(x,y)と、y方向のシフト量dy(x,y)とにより規定され
る。なお、図20では格子状の交点部分が各画素位置で
あり、処理前画像の注目画素(xc,yc)に対するx方向の
シフト量dx(xc,yc)と、y方向のシフト量dy(xc,yc)とに
より規定されるシフト方向とシフト量とを示している。
このように、シフト量生成部80は、シフト量とともに
シフト方向(以下、これらを総称して「シフトデータ」
とする)も決定する。このようなシフトデータの決定に
は以下のような方法を採用できる。
る。
の乱数を発生する関数で、kは最大シフト量(通常は0
<k≦1)である。なお、dx(x,y)を決定する際のR(-1.
0〜+1.O)とdy(x,y)を決定する際のR(-1.0〜+1.0)とは
別個であり、決定されたdx(x,y)の値とdy(x,y)の値とは
必ずしも同じではない。
に決定することにより、シフト方向とシフト量とをラン
ダムに決定することができる。生成されたシフト量は、
-k〜+kの間の値になる。
定する、
で、Fは周期関数の周波数(通常は、0<F≦1)、R(O.
0〜+1.O)はO.0〜+1.Oの間の乱数を発生する関数、RR
(x)はx座標における周期関数の初期値を規定するデー
タ、RR(y)はy座標における周期関数の初期値を規定
するデータである。
方向と直交する方向に周期性を有する関数であり、dx
(x,y)の決定においてはx座標が同じdx(x,y)は全て同じ
周期関数によって決定し、dy(x,y)の決定においてはy
座標が同じdy(x,y)は全て同じ周期関数によって決定す
る。周期関数を用いた場合であっても、生成されたシフ
ト量は、-k〜+kの間の値になる。
os関数など他の周期性を有するdx(x,y)、dy(x,y)が決定
できる関数を用いることができる。
5,+0.5]や[+0.5,-0.5]、[-0.5,+0.5]、[-0.5,-0.5]など
に固定する。この場合、生成されたシフト量は、固定さ
れた[dx(x,y),dy(x,y)]によって規定される値となる。
方法を述べたが、シフト量生成部80がシフトデータを
決定するに際して上記いずれの方法を採用するかは、既
述した周期性データ、より具体的には周期値T(図9)
を参照して決定しても良い。例えば、画像パターンが2
画素周期の周期性(周期値T=2)を有するときにはの
方法を選択するのが有効である。また、本発明者の実験
により、画像パターンが所定の周期値の場合にの方法
でシフトデータを決定すると良好な結果が得られないこ
とがあったので、例えば、周期値Tが2以外で、かつ、
上記所定の周期値以外の場合はの方法を選択し、上記
所定の周期値の場合はの方法を選択することなども考
えられる。
として処理してもよいし、入力装置6からオペレータに
より設定されたデータを用いて処理してもよいし、さら
に、周期性データに基づき決定してもよい。例えば、周
期関数の周期と、画像パターンの周期とが干渉するのを
避けるような周波数Fを決定することが好ましい。
囲)による乱数の発生頻度は乱数発生範囲内で均一にし
てもよいし、例えば、0.0付近の乱数を発生し易くした
り、+1.O付近の乱数を発生し易くしたりするなど、乱数
の発生頻度を偏重させるようにしてもよい。この乱数の
発生頻度の設定は、入力装置6からオペレータにより設
定可能に構成してもよい。
x方向と一致していれば、dy(x,y)=0として、dx(x,y)だ
けでシフトデータを決定したり、画像パターンの周期方
向がy方向と一致していれば、dx(x,y)=0として、dy(x,
y)だけでシフトデータを決定したりすることにより、画
像パターンの周期方向だけに注目画素の位置がシフトさ
れるようにシフトデータのシフト方向を決定してもよ
い。
シフトデータを決定し、シフト量が生成される。そし
て、生成されるシフト量は、上記いずれの方法によった
としても、-k〜+kの間の値になる。
フト処理時における実際のシフト量を規定するためのシ
フト量規定値を算出する。シフト量規定値は、上述した
ぼかし量規定値(C×PM)に基づいてルックアップテ
ーブルにより求められた値に周期性指標Mを乗じて得ら
れる値である。
れるルックアップテーブルを示す図である。ぼかし量規
定値(C×PM)は上述のぼかし量決定に用いたものと
同じであり、注目画素(xc,yc)についての加工後周期性
指標PM(xc,yc)とコントラスト量C(xc,yc)との積であ
る。図21に示すように、ぼかし量規定値(C×PM)
が所定の値よりも小さいときは、ルックアップテーブル
により求められる値は1.0である。そして、ぼかし量規
定値(C×PM)が所定の値よりも大きいときは、ぼか
し量規定値(C×PM)が大きくなるに従ってルックア
ップテーブルにより求められる値は小さくなり、ある値
以上になるとルックアップテーブルにより求められる値
が0.0になる。すなわち、このルックアップテーブルの
特性を定性的に表現すれば、ぼかし量規定値(C×P
M)が大きくなるに従って、小さな値を出力するもので
あると言える。
テーブルにより求められる値に周期性指標Mを乗じてシ
フト量規定値Sを得る。シフト量規定値Sも注目画素ご
とに算出されるものであり、注目画素(xc,yc)について
シフト量規定値S(xc,yc)が決定されるのである。な
お、周期性指標Mは0以上1以下の値であり、図19の
ルックアップテーブルにより求められる値も0以上1以
下の値であるため、それらの積であるシフト量規定値S
も0以上1以下の値となる。
生成部80が生成したシフト量にシフト量規定値Sを乗
ずることによって、最終的なシフト量、すなわちシフト
処理時における実際のシフト量を算出する。なお、最終
的なシフト量の算出も注目画素(xc,yc)ごとに行われる
ものである。また、シフト量生成部80が生成するシフ
ト量は−k以上k以下の値であり、シフト量規定値Sは0
以上1以下の値であるため、最終的なシフト量は−k以
上k以下の値となる。
いてのシフト量が算出され、シフト方向も決定されるの
である。
注目画素(xc,yc)についてのぼかし量およびシフト量が
上述の如く算出され、シフト方向が決定された後、その
ぼかし量に従ってぼかし処理部95が処理前画像のぼか
し処理を行い、そのシフト量およびシフト方向に従って
シフト処理部90がシフト処理をこの順序で行う(図3
のステップS7,S8)。
とは、ぼかし量決定部60が選択した加重平均フィルタ
(図19参照)を用いたフィルタ操作による処理であ
る。このフィルタ操作自体は、一般的なフィルタリング
処理であり、ぼかし量決定部60が選択した加重平均フ
ィルタを処理前画像の注目画素(xc,yc)に対して作用さ
せるものである。具体的には、ぼかし量決定部60が選
択した図19の如き加重平均フィルタを構成する各画素
の係数を図22のように表すとすると、この加重平均フ
ィルタからの出力Qb(xc,yc)は、次の数8に従って算
出され、この値が注目画素(xc,yc)の新たな階調値とな
る。
(xc,yc)であることは上述の通りである。例えば、図1
9(a)に示す加重平均フィルタが選択されたとする
と、注目画素(xc,yc)を中心とする5×5画素の領域内
の各画素の階調値のそれぞれに1/25を乗じ、これら
を積算した値がフィルタの出力Qb(xc,yc)であり、注
目画素(xc,yc)の新たな階調値となる。
0が行うシフト処理は、上記の決定されたシフト方向お
よびシフト量に従って、注目画素(xc,yc)の位置を移動
(シフト)させる処理である。このときに、シフト後の
位置に応じた注目画素(xc,yc)の新たな階調値Q(xc,yc)
は、シフト後の位置の近傍画素の階調値を補間して求め
る。この補間法は、一般的なバイリニア法などを採用す
ることができる。
図20に示すようなシフト方向およびシフト量(シフト
量規定値S=1とする)に従って移動した場合、ぼかし
処理後の階調値Qb(xc,yc)、Qb(xc+1,yc)、Qb(xc,
yc+1)、Qb(xc+1,yc+1)に基づき、これら4つ画素位置
((xc,yc)、(xc+1,yc)、(xc,yc+1)、(xc+1,yc+1))の間
で階調値が線形的に変化しているものと仮定して、移動
後の注目画素(xc+dx(xc,yc),yc+dy(xc,yc))の階調値を
算出する。より具体的には、ぼかし処理後の階調値Qb
(xc,yc)、Qb(xc+1,yc)から階調値Qb(xc+dx(xc,yc),
yc)を直線近似にて算出し、ぼかし処理後の階調値Qb
(xc,yc+1)、Qb(xc+1,yc+1)から階調値Qb(xc+dx(xc,
yc),yc+1)を直線近似にて算出する。その後、これらの
階調値Qb(xc+dx(xc,yc),yc)および階調値Qb(xc+dx
(xc,yc),yc+1)から移動後の注目画素の階調値Qb(xc+d
x(xc,yc),yc+dy(xc,yc))を直線近似にて算出するのであ
る。
において、ぼかし量規定値(C×PM)が大きい場合に
は、平滑化の程度が大きい加重平均フィルタ(図19
(a)参照)が選択されて大きなぼかし量のぼかし処理
が行われる一方、図21のルックアップテーブルにより
シフト量が0になり、シフト処理が行われないこととな
る。ぼかし量規定値(C×PM)が大きい場合とは、注
目画素の周辺において画像パターンの繰り返しの有無の
程度が大きく、かつ高コントラストであることを意味し
ており、例えばエアコンのフロントパネルについての画
像等が該当する。
中の輪郭部分(高コントラストを発生する部分)が乱れ
て画質が劣化するのであるが、本実施形態のようにすれ
ば、大きなぼかし量のぼかし処理のみが行われてシフト
処理は実質的に行われないため、画質の劣化はほとんど
生じない。なぜなら、高コントラスト部分はUSM強調
する必要がなく、ぼかし処理をおこなっても画質が劣化
しないためである。そして、本処理後の画像に対してシ
ャープネス処理等を施したとしても、サチレーションが
抑制もしくは解消されるため、先に説明した平均レベル
の変動すなわち、平均階調値が周期的に変動する状態は
生じず、モアレの発生を抑制することができる。すなわ
ち、画像パターンの繰り返しの有無の程度が大きく、か
つ高コントラストの画像であったとしても、その画質を
劣化させることなくモアレの発生を抑制することができ
るのである。
ンの繰り返しの有無の程度は大きいが、低コントラスト
である場合、例えば衣服の布地についての画像等の場合
は、コントラスト量Cの値が小さくなるため、ぼかし量
規定値(C×PM)の値も小さくなる。従って、平滑化
の程度が小さい加重平均フィルタ(図19(b)参照)
が選択されてほとんどぼかし処理が行われない一方、シ
フト量規定値Sの値は大きくなり、大きなシフト処理が
行われる。
きいが、低コントラストの画像である場合には、大きな
シフト処理を行っても画質の劣化はほとんどない。なぜ
なら、このようなテクスチャ領域はぼかし処理ではテク
スチャが消滅するなどの問題が発生するが、シフト処理
ではテクスチャは消滅しないからである。しかも「IN
フェイズ状態」と「OUTフェイズ状態」との周期性が
崩れ、「INフェイズ状態」と「OUTフェイズ状態」
とがランダムに分散される。従って、本処理後の画像に
対してシャープネス処理等を施したとしても、平均階調
値が周期的に変動する状態は生じず、モアレの発生を抑
制することができる。すなわち、画像パターンの繰り返
しの有無の程度が大きく、かつ低コントラストの画像で
あったとしても、その画質を劣化させることなくモアレ
の発生を抑制することができるのである。
度が小さな画像については、周期性指標Mおよび加工後
周期性指標PMの値が小さくなるため、そのコントラス
トの高低の如何に関わらず、ぼかし量規定値(C×P
M)の値およびシフト量規定値Sの値ともに小さくな
り、ぼかし処理もシフト処理も行われない。画像パター
ンの繰り返しの有無の程度が小さな画像についは、そも
そもモアレの発生がほとんどないものと考えられるた
め、画質の低下を回避する観点からぼかし処理もシフト
処理も行わないのが好ましいのである。
もある。文字には、例えば漢字の「日」、「口」、
「田」のように周期性を有するものもある。このような
周期性を有する文字については、画像パターンの繰り返
しであると判断され、周期性指標Mの値が大きくなる。
一方、文字については、高コントラストであるため、コ
ントラスト量Cの値も大きくなる。すなわち、文字につ
いては、周期性があり、かつ高コントラストであるた
め、例えばエアコンのフロントパネルと同様にぼかし処
理の対象にすると、いわゆる文字ぼけが生じることな
り、適当でない。なお、加工部40によって差分画像を
加工したとしても、面積の小さな領域(例えば文字が単
独で存在するような場合)であれば除去できるが、面積
の大きな領域(例えば、複数の文字からなる文章として
存在する場合)については除去することができない。
期性指標算出部25が周期性指標画像に対して加工(図
10)を行った後、周期性指標を検出して加工後周期性
指標としている。加工後周期性指標算出部25による加
工において、周期性指標が急激に変化している局所領域
が除去されるのは既述したとおりである。文章中に周期
性を有する文字が点在している(周期性を有さない文字
に周期性を有する文字が混じっている)場合には、周期
性を有する文字についての周期性指標Mの値は大きくな
り、その周辺の周期性を有さない文字についての周期性
指標Mの値は小さくなる。つまり、周期性を有する文字
部分のみが周期性指標Mが急激に上昇する局所領域とな
り、加工後周期性指標算出部25に除去されることとな
る。具体的には、周期性を有する文字部分の大きな周期
性指標Mの値が、加工によって減少され、加工後周期性
指標PMの値は小さくなるのである。
は、ぼかし量規定値(C×PM)の値が小さくなり、平
滑化の程度が小さい加重平均フィルタ(図19(b)参
照)が選択されてほとんどぼかし処理が行われない。よ
って、いわゆる「文字ぼけ」による画質の劣化を防止す
ることができるのである。
画像の内容に応じてぼかし処理およびシフト処理の程度
を変化させることにより、種々の種類の画像について画
質を劣化させることなくモアレの発生を抑制することが
できるのである。
装置が上記第1実施形態の画像処理装置と相違するの
は、抽出された周期画像領域を周期性指標算出部20が
拡張し、その拡張後の拡張周期画像領域における注目画
素について周期性指標を算出する点である。画像処理装
置のハードウェア構成等の残余の点については上記第1
実施形態と同じである。
1実施形態と同様の処理を行う場合について考える。な
お、図23において、斜線部は黒色領域であり、残余の
部分は白色領域であるものとする。
域P2に存在する場合、自己相関データH(a)は図24
にて示すようになる。すなわち、黒色領域P2の両隣に
は、黒色領域P1および黒色領域P3が存在し、画像パ
ターンの繰り返しとなっているため、所定のしきい値以
下の極小値H(a)が規則的に存在する。その結果、黒色
領域P2は周期画像領域として抽出され、注目画素(xc,
yc)には高い値(例えば”1.0”)の周期性指標M(x
c,yc)が付与される。注目画素(xc,yc)が黒色領域P3に
存在する場合も同様である。
黒色領域P1に存在する場合、自己相関データH(a)は
図25にて示すようになる。すなわち、黒色領域P1の
右隣には黒色領域P2が存在するものの、左隣には黒色
領域が存在しないため、所定のしきい値以下の極小値H
(a)が1つしか存在しない、つまり規則的に存在しない
こととなる。その結果、黒色領域P1は周期画像領域と
して抽出されず、注目画素(xc,yc)の周期性指標M(xc,y
c)は”0”とされる。注目画素(xc,yc)が黒色領域P4
に存在する場合も同様である。
23に示した画像から図26の斜線部にて示すような周
期画像領域が抽出されることとなる。すなわち、図23
の処理前画像における画像パターンの繰り返しよりも狭
い領域が周期画像領域として抽出される。
性指標Mに基づいてシフト処理等が行われたとすると、
周期性指標Mの値が大きな黒色領域P2、P3について
は大きなシフト処理が行われるのに対して、周期性指標
Mの値が小さな黒色領域P1、P4についてはシフト処
理がほとんど行われない。従って、本処理後の画像に対
してシャープネス処理等を施したとすると、鮮鋭化の程
度が黒色領域P2、P3では比較的弱いのに対して、黒
色領域P1、P4では非常に強く現れることとなる。こ
のことは、換言すれば、均質であった処理前画像が本処
理によって不均質な画像になることを意味する。
れた周期画像領域を周期性指標算出部20が画像パター
ンの繰り返しの1周期分、すなわち周期値Tだけ拡張し
て拡張周期画像領域を取得しているのである。
注目画素(xc,yc)を中心として±Tの範囲内に存在する
画素の周期性指標Mと当該注目画素(xc,yc)の周期性指
標M(xc,yc)とを比較し、±Tの範囲内に存在する画素
の周期性指標Mの方が小さければ、その画素の周期性指
標Mの値を注目画素(xc,yc)の周期性指標M(xc,yc)の値
に置き換えているのである。例えば、図26において抽
出された周期画像領域内に含まれる黒色領域P2の注目
画素(xc,yc)を中心とした±Tの範囲内には、黒色領域
P1が含まれる。黒色領域P1内の画素の周期性指標M
は、黒色領域P2の注目画素(xc,yc)の周期性指標M(x
c,yc)よりも低いため、黒色領域P1内の画素の周期性
指標Mの値が黒色領域P2の注目画素(xc,yc)の周期性
指標M(xc,yc)の値に置換されるのである。同様にし
て、黒色領域P4内の画素の周期性指標Mの値が黒色領
域P3の注目画素(xc,yc)の周期性指標M(xc,yc)の値に
置換される。
な拡張周期画像領域が取得されることとなる。この拡張
周期画像領域は、本来周期画像領域として抽出されるべ
き図23の処理前画像における画像パターンの繰り返し
の全てを含むものである。そして、拡張周期画像領域内
の画素には、元の周期画像領域の画素が有していた周期
性指標Mの値と同じ値の周期性指標Mが付与されてい
る。
標算出部20は、拡張周期画像領域の境界部分のなじま
せ処理を行う。図28は、拡張周期画像領域の境界部分
のなじませ処理を説明する図である。このなじませ処理
は、上記第1実施形態におけるぼかし処理と同様の処理
であり、図29に示すようなフィルタを用いたフィルタ
操作による処理である。具体的には、拡張後の周期性指
標Mによって表現された画像に対して図29に示すよう
なフィルタを順次に走査させる処理である。
拡張周期画像領域の境界部分において、周期性指標Mの
値が段階的に推移する状態、すなわちにじむ状態とな
る。
は、取得した拡張周期画像領域における注目画素につい
て、画像パターンの繰り返しの有無の程度を示す周期性
指標を算出しているのである。なお、上述においては、
説明を簡単にするため、一方向のみに拡張する場合につ
いて説明したが、2方向に拡張した後それを合成するこ
とによって拡張周期画像領域を取得するようにしても良
い。拡張周期画像領域における注目画素について周期性
指標が算出された後の処理は、上記第1実施形態と同じ
であるためその説明は省略する。
の効果が得られるのに加えて、図23に示したような画
像が処理対象であった場合でも、黒色領域P1、P2、
P3、P4の全てについてほぼ均一な値の周期性指標M
が付与される結果、シフト処理等の程度も同程度とな
る。従って、本処理後の画像に対してシャープネス処理
等を施したとしても、黒色領域P1、P2、P3、P4
のそれぞれについての鮮鋭化の程度も同程度となり、画
質の劣化が生じない均質な画像処理が行われることとな
る。
ませ処理を行うことによって、周期性指標Mの値が急激
に変化する状態が生じないため、シフト処理等の程度が
急減に変化することもなく、画質の著しい劣化を防止す
ることができる。
じませ処理を行うことによって、拡張周期画像領域にお
ける周期性を有するノイズ部分(例えば、上述した周期
性を有する文字が単独で存在する部分)を除去すること
もできるのである。このノイズ除去の意味においては、
加工後周期性指標算出部25による加工と同様の効果
(いわゆる「文字ぼけ」の防止)を得ることができる。
なお、拡張周期画像領域におけるノイズ除去という観点
からは、フィルタ操作に使用するフィルタとしては図2
9に示すようなものに限定されず、例えばメディアンフ
ィルタを使用するようにしても良い。
期性指標算出部25による加工および第2実施形態にて
説明した周期画像領域の拡張は、それらの双方を実施す
るようにしても良いし、それぞれを単独で実施するよう
にしてもよい。双方を実施した場合には、加工後周期性
指標算出部25による加工がぼかし処理に強く影響する
ため、周期画像領域の拡張は主としてシフト処理に影響
することとなる。なお、周期画像領域の拡張のみを行う
場合には、加工後周期性指標が算出されないため、第1
実施形態にて説明した加工後周期性指標の代わりに周期
性指標をそのまま用いてぼかし量およびシフト量を算出
することとなる。
によれば、周期性指標によって表現された周期性指標画
像を複数の領域に分割し、複数の領域のそれぞれに含ま
れる画素の階調値を平均化して1つの画素の階調値とし
て出力することにより、周期性指標画像をバイリニア縮
小し、そのバイリニア縮小によって得られた縮小画像に
対してメディアンフィルタによるフィルタ操作を施し、
フィルタ操作後の縮小画像に含まれる画素の階調値を補
間しつつ、当該縮小画像を周期性指標画像と同等の大き
さに拡大することにより、加工後周期性指標画像を取得
し、当該加工後周期性指標画像における注目画素の加工
後周期性指標を算出しているため、周期性を有する文字
部分等についてのぼかし処理の程度が小さくなり、画質
の劣化を防止しつつモアレの発生を抑制することができ
る。
領域を画像パターンの繰り返しの1周期分拡張して拡張
周期画像領域を取得し、その拡張周期画像領域における
注目画素について、画像パターンの繰り返しの有無の程
度を示す周期性指標を算出しているため、画像パターン
の繰り返しが存在する領域の全てについて均質なシフト
処理等が行われることとなり、画質の劣化を防止しつつ
モアレの発生を抑制することができる。
の一例を示すブロック図である。
ロック図である。
順の概要を示すフローチャートである。
る。
例を示す図である。
例を示す図である。
から得られたx方向に沿った自己相関データを示す図で
ある。
一例を示す図である。
期値と実際の周期方向及び周期値との関係を示す図であ
る。
ャートである。
大を概念的に説明するための図である。
ートである。
る図である。
度領域拡張画像を示す図である。
度領域縮小画像を示す図である。
を概念的に説明するための図である。
る。
す図である。
ップテーブルを示す図である。
ある。
示す図である。
例を示す図である。
の例を示す図である。
示す図である。
理を説明する図である。
示す図である。
画像パターンと光電素子との位置関係を示す図である。
ある場合の「INフェイズ状態」と「OUTフェイズ状
態」とを示す図である。
施した場合の計算上の各画素の階調値を示す図である。
施した場合の各画素の最終的な階調値を示す図である。
ルックアップテーブルの一例を示す図である。
ルックアップテーブルの他の例を示す図である。
Claims (2)
- 【請求項1】 画像に対して処理を施す画像処理装置で
あって、 画像から周期的な画像パターンの繰り返しが存在する周
期画像領域を抽出する周期画像領域抽出手段と、 前記周期画像領域における注目画素について、当該注目
画素を含む注目領域と当該注目領域以外の周辺領域との
相関特性に基づいて画像パターンの繰り返しの有無の程
度を示す周期性指標を算出する周期性指標算出手段と、 前記周期性指標によって表現された周期性指標画像を加
工することによって加工後周期性指標画像を取得し、当
該加工後周期性指標画像における前記注目画素の加工後
周期性指標を算出する加工後周期性指標算出手段と、 前記画像における注目画素とその周辺画素との明暗比を
示す指標であるコントラスト量を算出するコントラスト
量算出手段と、 前記加工後周期性指標および前記コントラスト量に基づ
いて、前記画像における注目画素についてのぼかし処理
の程度であるぼかし量を決定するぼかし量決定手段と、 前記周期性指標、前記加工後周期性指標および前記コン
トラスト量に基づいて、前記画像における注目画素につ
いてのシフト処理の程度であるシフト量を決定するシフ
ト量決定手段と、 前記ぼかし量決定手段により決定されたぼかし量に基づ
いて、前記画像における注目画素についてのぼかし処理
を行うぼかし処理手段と、 前記シフト量決定手段により決定されたシフト量に基づ
いて、前記画像における注目画素についてのシフト処理
を行うシフト処理手段と、を備え、 前記加工後周期性指標算出手段は、 前記周期性指標画像内を複数の領域に分割し、前記複数
の領域のそれぞれに含まれる画素の階調値を平均化して
1つの画素の階調値として出力することにより、前記周
期性指標画像をバイリニア縮小する縮小手段と、 前記バイリニア縮小によって得られた縮小画像に対して
メディアンフィルタによるフィルタ操作を施すフィルタ
手段と、 前記フィルタ操作後の縮小画像に含まれる画素の階調値
を補間しつつ、当該縮小画像を前記周期性指標画像と同
等の大きさに拡大することにより、前記加工後周期性指
標画像を得る拡大手段と、を含むことを特徴とする画像
処理装置。 - 【請求項2】 画像に対して処理を施す画像処理装置で
あって、 画像から周期的な画像パターンの繰り返しが存在する周
期画像領域を抽出する周期画像領域抽出手段と、 前記周期画像領域を前記画像パターンの繰り返しの1周
期分拡張して拡張周期画像領域を取得し、前記拡張周期
画像領域における注目画素について、画像パターンの繰
り返しの有無の程度を示す周期性指標を算出する周期性
指標算出手段と、 前記画像における注目画素とその周辺画素との明暗比を
示す指標であるコントラスト量を算出するコントラスト
量算出手段と、 前記周期性指標および前記コントラスト量に基づいて、
前記画像における注目画素についてのぼかし処理の程度
であるぼかし量を決定するぼかし量決定手段と、 前記周期性指標および前記コントラスト量に基づいて、
前記画像における注目画素についてのシフト処理の程度
であるシフト量を決定するシフト量決定手段と、 前記ぼかし量決定手段により決定されたぼかし量に基づ
いて、前記画像における注目画素についてのぼかし処理
を行うぼかし処理手段と、 前記シフト量決定手段により決定されたシフト量に基づ
いて、前記画像における注目画素についてのシフト処理
を行うシフト処理手段と、を備えることを特徴とする画
像処理装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP25618199A JP3734652B2 (ja) | 1999-09-09 | 1999-09-09 | 画像処理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP25618199A JP3734652B2 (ja) | 1999-09-09 | 1999-09-09 | 画像処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2001086368A true JP2001086368A (ja) | 2001-03-30 |
JP3734652B2 JP3734652B2 (ja) | 2006-01-11 |
Family
ID=17289043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP25618199A Expired - Fee Related JP3734652B2 (ja) | 1999-09-09 | 1999-09-09 | 画像処理装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3734652B2 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002300383A (ja) * | 2001-03-30 | 2002-10-11 | Ricoh Co Ltd | 画像処理装置、画像処理方法、画像処理方法をコンピュータに実行させるプログラム、およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP2003158644A (ja) * | 2001-11-20 | 2003-05-30 | Toppan Printing Co Ltd | 一定周期パターンを持つ平面状被写体の撮像装置 |
CN100386797C (zh) * | 2005-11-03 | 2008-05-07 | 西安交通大学 | 一种利于实现数字视频图像缩放插值的行场信号产生方法 |
JP2022160498A (ja) * | 2018-03-22 | 2022-10-19 | カシオ計算機株式会社 | 投影制御装置、投影装置、投影制御方法及びプログラム |
-
1999
- 1999-09-09 JP JP25618199A patent/JP3734652B2/ja not_active Expired - Fee Related
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002300383A (ja) * | 2001-03-30 | 2002-10-11 | Ricoh Co Ltd | 画像処理装置、画像処理方法、画像処理方法をコンピュータに実行させるプログラム、およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP4605931B2 (ja) * | 2001-03-30 | 2011-01-05 | 株式会社リコー | 画像処理装置、画像処理方法、画像処理方法をコンピュータに実行させるプログラム、およびそのプログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP2003158644A (ja) * | 2001-11-20 | 2003-05-30 | Toppan Printing Co Ltd | 一定周期パターンを持つ平面状被写体の撮像装置 |
CN100386797C (zh) * | 2005-11-03 | 2008-05-07 | 西安交通大学 | 一种利于实现数字视频图像缩放插值的行场信号产生方法 |
JP2022160498A (ja) * | 2018-03-22 | 2022-10-19 | カシオ計算機株式会社 | 投影制御装置、投影装置、投影制御方法及びプログラム |
JP7380771B2 (ja) | 2018-03-22 | 2023-11-15 | カシオ計算機株式会社 | 投影制御装置、投影装置、投影制御方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP3734652B2 (ja) | 2006-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6915024B1 (en) | Image sharpening by variable contrast mapping | |
JP4727720B2 (ja) | 画像処理方法および画像処理装置 | |
US7792384B2 (en) | Image processing apparatus, image processing method, program, and recording medium therefor | |
JP2010009381A (ja) | 画素補間装置、画素補間方法および画素補間プログラム | |
JP2007527567A (ja) | 領域エッジのシャープネス補正を伴う画像シャープニング | |
JP2004040247A (ja) | 画像処理装置及びその制御方法 | |
JPH07131634A (ja) | 画像処理装置 | |
JPH11195113A (ja) | 画像処理装置、方法及びコンピュータ読み取り可能な記憶媒体 | |
JP3734652B2 (ja) | 画像処理装置 | |
JP3092769B2 (ja) | 画像処理装置 | |
JPH0944128A (ja) | 2値画像多値化縮小処理装置 | |
JP3751806B2 (ja) | 画像処理装置および記録媒体 | |
JP6444049B2 (ja) | 画像処理装置、方法およびプログラム | |
KR102470242B1 (ko) | 영상 처리 장치, 영상 처리 방법, 및 프로그램 | |
JP3753900B2 (ja) | 画像処理装置 | |
JP3262425B2 (ja) | 画像処理装置 | |
JP3200351B2 (ja) | 画像処理装置及びその方法 | |
JP2004133592A (ja) | 画像を拡大するための画像処理装置および画像処理方法および画像処理プログラム | |
JP3624153B2 (ja) | 画像処理装置及び画像処理方法 | |
JP2001078039A (ja) | 画像処理装置および記録媒体 | |
JP3581460B2 (ja) | 画像処理方法とその装置 | |
JP2005293265A (ja) | 画像処理装置及び方法 | |
JP2015106318A (ja) | 画像処理装置および画像処理方法 | |
JP3714820B2 (ja) | 画像処理方法 | |
JP2006050481A (ja) | 画像処理装置、画像処理方法及びそのプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051007 |
|
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: 20051018 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20051019 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091028 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091028 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091028 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101028 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101028 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111028 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111028 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121028 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121028 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131028 Year of fee payment: 8 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |