JP4086669B2 - 画像処理方法及び装置 - Google Patents

画像処理方法及び装置 Download PDF

Info

Publication number
JP4086669B2
JP4086669B2 JP2003016713A JP2003016713A JP4086669B2 JP 4086669 B2 JP4086669 B2 JP 4086669B2 JP 2003016713 A JP2003016713 A JP 2003016713A JP 2003016713 A JP2003016713 A JP 2003016713A JP 4086669 B2 JP4086669 B2 JP 4086669B2
Authority
JP
Japan
Prior art keywords
value
image data
noise
quantization
quantization error
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.)
Expired - Fee Related
Application number
JP2003016713A
Other languages
English (en)
Other versions
JP2004229133A5 (ja
JP2004229133A (ja
Inventor
尚 石川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Inc
Original Assignee
Canon Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Canon Inc filed Critical Canon Inc
Priority to JP2003016713A priority Critical patent/JP4086669B2/ja
Priority to US10/611,948 priority patent/US8009327B2/en
Publication of JP2004229133A publication Critical patent/JP2004229133A/ja
Publication of JP2004229133A5 publication Critical patent/JP2004229133A5/ja
Application granted granted Critical
Publication of JP4086669B2 publication Critical patent/JP4086669B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Color, Gradation (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は画像処理方法及び装置に関し、特に入力画像濃度と出力画像濃度等の差を誤差拡散法等により保存しつつ、入力データを2値又は多値データに量子化処理する画像処理方法及び装置に関する。
【0002】
【従来の技術】
従来より、入力された多階調画像データを2値画像データに変換して出力する画像処理装置では、多階調画像を2値画像に変換する方法として、例えば、誤差拡散法が用いられている。この誤差拡散法は、注目画素の階調値と、この階調値を2値化処理するための固定の2値化閾値との量子化誤差が、注目画素の近傍画素の階調値に拡散され、これら各階調値を逐次的に2値化する方法である(例えば、非特許文献1参照)。
【0003】
【非特許文献1】
Proceeding of the S.I.D., Vol 17/2, 1976, R.Floyd and L.Steinberg, "An Adaptive Algorithm for Spatial Greyscale", pp.75-76。
【0004】
ところが、誤差拡散法を利用した2値化方法では、誤差拡散処理特有の、一般的に“テクスチャ”や“ワーム”と呼ばれる縞模様が再生画像に発生し、視覚的印象に悪影響を与える場合があった。
【0005】
また、濃度変化の大きな部分においてはドット生成の遅れによる“はき寄せ”と呼ばれる白とび(高濃度から低濃度への変化部)、黒つぶれ(低濃度から高濃度への変化部)が再生画像に発生し、濃度の再現性を損なう場合があった。
【0006】
さらに、誤差拡散法を2値以上の多値に拡張した場合、量子化代表値付近においてはドットが消失してしまうため、擬似輪郭として認識されてしまう場合があった。
【0007】
このような問題を解決するために従来からも、入力画像から一定幅の均一濃度の画素データを検出し、そこに乱数により生成したノイズを重畳することによりテクスチャパターンの発生を防止する方法が提案されている(例えば、特許文献1参照)。
【0008】
さらに、画像情報の入力値に応じて、その最大値を変化させた乱数を2画素以上の周期で生成し、この乱数を2画素以上の周期内で、その負号を正/負に反転した信号と画像情報の入力値とを加算し、加算された信号値をn値化することにより、上述した“はき寄せ”や“テクスチャ”などの発生を防止する方法も提案されている(例えば、特許文献2参照)。
【0009】
【特許文献1】
特開平11−146201号公報。
【0010】
【特許文献2】
特開平9−179974号公報。
【0011】
【発明が解決しようとする課題】
しかしながら上記従来例では、入力画像にノイズを重畳しているため、重畳したノイズの影響が出力画像に残ってしまうという欠点があった。
【0012】
このノイズの影響を出力画像に及ぼさないためには、重畳するノイズのDC成分は“0”でなければならないが、このようなノイズは“はき寄せ”や“テクスチャ”の発生を抑止することが困難であった。また、出力画像の粒状感を少なくするためには、ノイズの低周波成分を抑圧し、かつ、周期性を崩す必要があるが、このような要求を満たすノイズの生成は困難であった。
【0013】
さらに、“テクスチャ”や“ワーム”を排除するために上記のようなノイズの振幅を大きくして重畳すると、そのノイズが出力画像に残ってしまい、粒状感の目立つ、質の悪い画像が形成されてしまう。従って、最適な振幅のノイズを付与することも困難であった。
【0014】
本発明は上記従来例に鑑みてなされたもので、誤差拡散によって生じる、例えば、“テクスチャ”、“ワーム”、“はき寄せ”等の画質を劣化させる要因を抑え、画像に重畳されるノイズの影響を最小化する画像処理方法及び装置を提供することを目的としている。
【0015】
【課題を解決するための手段】
上記目的を達成するため本発明の画像処理方法は以下のような工程からなる。
【0016】
即ち、n値画像データに誤差拡散処理を施して擬似中間調画像を形成するためのm値画像データ(2≦m<n)に変換する画像処理方法であって、n値画像データを入力する入力工程と、前記入力されたn値画像データに以前の誤差拡散によって累積した量子化誤差を加算する第1加算工程と、ノイズの絶対値を発生する発生工程と、前記発生工程において発生するノイズの値を予め定められた値以下に制限する制限工程と、前記入力されたn値画像データの特定ビットの値によって前記ノイズの符号を決定する決定工程と、前記決定工程において符号が決定されたノイズを前記第1加算工程において量子化誤差が加算されたn値画像データにさらに加算する第2加算工程と、前記第2加算工程においてノイズが加算されたn値画像データをm値化する量子化工程と、前記量子化工程においてm値化されたm値データを量子化代表値に逆量子化する逆量子化工程と、前記第1加算工程において量子化誤差が加算されたn値画像データから前記逆量子化工程において得られた量子化代表値を減算して量子化誤差を生成する量子化誤差生成工程と、前記量子化誤差生成工程において生成した量子化誤差をバッファに格納する格納工程とを備えることを特徴とする。
【0017】
また、本発明は、上記画像処理方法における各工程をコンピュータによって実行するために、コンピュータが実行可能なコードによって記述されたプログラムの形態であっても良い。
【0018】
さらに、そのプログラムは、コンピュータが読み取り可能であるように、コンピュータ可読の記憶媒体に格納するように備えられても良い。
【0019】
このようにして、プログラムや記憶媒体の形で本発明を実現することも可能である。
【0020】
またさらに、本発明は、上記画像処理方法を実行できる画像処理装置の形態であっても良い。
【0021】
即ち、n値画像データに誤差拡散処理を施して擬似中間調画像を形成するためのm値画像データ(2≦m<n)に変換する画像処理装置であって、n値画像データを入力する入力手段と、前記入力手段によって入力されたn値画像データに以前の誤差拡散によって累積した量子化誤差を加算する第1加算手段と、ノイズの絶対値を発生する発生手段と、前記発生手段により発生するノイズの値を予め定められた値以下に制限する制限手段と、前記入力されたn値画像データの特定ビットの値によって前記ノイズの符号を決定する決定手段と、前記決定手段によって符号が決定されたノイズを前記第1加算手段によって量子化誤差が加算されたn値画像データにさらに加算する第2加算手段と、前記第2加算手段によってノイズが加算されたn値画像データをm値化する量子化手段と、前記量子化手段によってm値化されたm値データを量子化代表値に逆量子化する逆量子化手段と、前記第1加算手段によって量子化誤差が加算されたn値画像データから前記逆量子化手段によって得られた量子化代表値を減算して量子化誤差を生成する量子化誤差生成手段と、前記量子化誤差生成手段によって生成された量子化誤差を格納するバッファとを備えることを特徴とするものである。
【0022】
【発明の実施の形態】
以上のような解決手段の構成をさらに詳しく言えば、前記第1加算工程ではバッファに格納された量子化誤差を所定の分配係数によって複数の画素に分配して累積された量子化誤差を加算することが望ましい。
【0023】
また、前記決定工程では、入力されたn値画像データの特定ビットの値によってノイズの負号を決定することが望ましい。
【0024】
さらに、前記発生工程において発生するノイズの値を所定値以下に制限する制限工程を備えることが好ましく、この制限は、リミッタやゲイン調整器などによって実行される。
【0025】
より具体的には、前記制限工程では、特定ビットが1であるなら、入力されたn値画像データの特定ビットより下位にあるビットによって表現される値、或いは、第1の所定値の内、いずれか小さい方をその所定値とする一方、特定ビットが0であるなら、入力されたn値画像データの特定ビットより下位にあるビットを反転することにより得られるビットによって表現される値、或いは、第1の所定値の内、いずれか小さい方を所定値とするようにして、その制限を実行すると良い。また、その制限された値が第2の所定値以下の場合、第2の所定値を前記所定値とすると良い。
【0026】
さらにまた、前記発生工程は、発生するノイズのビット数より長いシフトレジスタを用いたM系列擬似ランダム負号生成回路のシフトレジスタ出力より前記ノイズのビット数相当分を選択することにより、前記ノイズを生成すると良い。
【0027】
以下添付図面を参照して本発明の好適な実施形態について詳細に説明する。
【0028】
これ以降の説明では、まず多値画像データを入力して誤差拡散処理を行なってm値化データを出力するm値化処理装置のいくつかの実施形態について説明し、最後にこのm値化処理装置が応用される画像入出力システムについて説明する。この画像入出力システムでは出力装置としてインクジェットプリンタが用いられ、インクジェット記録ヘッドがカラーインクを吐出してカラー画像を記録する。
【0029】
<第1実施形態>
図1は本発明の第1実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。図1に示す構成では、後述する拡散係数を用いて1画素8ビットの多値画像データに5値の誤差拡散(即ち、m=5)を行う例を示している。ここでは、8ビットの多値画像データを5値化するための4つの閾値を“32”、“96”、“160”、“224”とする。
【0030】
図1において、加算器1では拡散フィルタ10からの補正値(e)を入力画像データ(x)に加算する。ノイズ生成回路2では、補正値が加算された画像データに重畳する1画素5ビットのノイズの絶対値を生成する。比較器3では入力画像データと量子化における閾値を比較し、その結果を1ビットの信号(s)として負号付加回路4に入力する。
【0031】
負号付加回路4では、比較器3による比較結果に従い、入力画像データが量子化器6の閾値よりも大きい場合はノイズ生成回路2の出力を負に変換して(具体的にはビット反転して+1する)加算器5に入力する。従って、負号付加回路4から出力される信号(n)は6ビットの信号である。一方、入力画像データが量子化器6の閾値よりも小さい場合はノイズ生成回路2の出力をそのまま加算器5に入力する。
【0032】
加算器5では、上記補正値が加算された画像データに上記負号が付加されたノイズを加算して、その加算結果(x+e+n)を量子化器6に入力する。量子化器6では、ノイズが付加された補正後の入力画像データを閾値と比較し、3ビットのコードに変換する。このコードは出力される一方、逆量子化器7にも入力される。逆量子化器7はそのコードより量子化代表値を生成して、その結果(r)を減算器8に入力する。減算器8では、量子化前のノイズが付加されていない値より量子化代表値を減算して7ビットの量子化誤差を生成し、これを誤差バッファ9に格納する。
【0033】
誤差バッファ9では、その量子化誤差を約1ライン分遅延し、拡散フィルタ10の演算に必要な画素の誤差値を出力する。拡散フィルタ10では、誤差バッファ9からの前ラインの拡散対象画素及び注目画素の量子化誤差に対応する係数を乗してその総和を加算器1に出力する。この出力が補正値として次の入力画像データに加算される。
【0034】
次に、以上の構成におけるm値化処理部の具体的な動作について説明する。
【0035】
この実施形態では、誤差バッファ9の容量を削減するため、誤差拡散後の値(拡散フィルタの値)を誤差バッファに格納するのではなく、量子化誤差そのものを誤差バッファに格納する構成になっている。
【0036】
図2は拡散フィルタ10の演算に用いる拡散係数の例を示す図である。
【0037】
図2に示したように、周囲の量子化誤差より次の画素の補正値を求める方法(平均誤差最小法)では誤差拡散法の拡散係数と位置が点対称の関係になる。
【0038】
ノイズ生成回路2は、量子化器6の閾値と量子化代表値との差分未満の値、即ち、0から31の値をランダムに生成する。
【0039】
また、比較器3は入力画像データと量子化器6の閾値及び逆量子化器7の量子化代表値とを比較し、ノイズの負号を決定する。ここでは入力画像データ(x)とノイズの負号を示す信号(s)との関係を以下のようにしている。
【0040】
即ち、 x<32 であれば、s=0、
32≦x<64 であれば、s=1、
64≦x<96 であれば、s=0、
96≦x<128 であれば、s=1、
128≦x<160 であれば、s=0、
160≦x<192 であれば、s=1、
192≦x<224 であれば、s=0、
224≦x であれば、s=1である。
【0041】
即ち、ノイズの負号を示す信号(s)は入力画像データ(x)のLSBから6ビット目の情報を抽出すれば良い。
【0042】
さて、量子化器6は10ビットで表現される入力データ(x+e+n)を5値のコード(c)に変換する。即ち、補正後の入力画像データ(x+e)に対して、上記のような閾値を考慮すると、3ビットの出力コード(c)は以下のようになる。
【0043】
即ち、 x+e<32 であれば、c=0、
32≦x+e<96 であれば、c=1、
96≦x+e<160であれば、c=2、
160≦x+e<224であれば、c=3、
224≦x+e であれば、c=4である。
【0044】
ここで、量子化器6の構成を簡略化するため、負号付加回路4では最大31のオフセットを加算したノイズ(n)を加算器5において補正後の入力画像データ(x+e)を加算する。即ち、ノイズ加算後の補正入力画像データ(x+e+n)に対して、出力コード(c)は以下のようになる。
【0045】
即ち、 x+e+n<64 であれば、c=0、
64≦x+e+n<128であれば、c=1、
128≦x+e+n<192であれば、c=2、
192≦x+e+n<256であれば、c=3、
256≦x+e+n であれば、c=4である。
【0046】
即ち、出力コード(c)はノイズ加算後の10ビット補正入力画像データの負号ビットを除いた上位3ビットを抽出すれば良い。但し、後述するように、補正値(e)の値によって出力コード(c)は、範囲外の値をとることがある。具体的には、(x+e+n)<0の場合と(x+e+n)>319の場合である。
【0047】
(x+e+n)<0の場合、負号ビットを除いた上位3ビットを抽出した値は“7”であり、c=7となるが、所望の出力コード(c)は“0”なので、その値が“7”の場合は“0”に変換して出力する。なお、量子化誤差の範囲は±63のため、c=6となることないので上位2ビットの値がともに“1”なら、c=0(即ち、3ビット全てを0)にしても良い。
【0048】
一方、(x+e+n)>319の場合、負号ビットを除いた上位3ビットを抽出した値は“5”であり、c=5となるが、所望の出力コード(c)は“4”なので、c=5の場合は、これをc=4に変換して出力する(即ち、上位2ビットの値が上位側から夫々、“1”、“0”ならLSBを“0”にする)。
【0049】
逆量子化器7は量子化器6の出力コード(c)を8ビットで表現される量子化代表値(r)に変換する。具体的には、次の通りである。
【0050】
即ち、c=0であれば、r=0、
c=1であれば、r=64、
c=2であれば、r=128
c=3であれば、r=192
c=4であれば、r=255である。
【0051】
従って、c=4以外(即ち、3ビットの出力コードのMSBが“0”となる場合)では、そのMSBを除いた2ビットの下位に“0”を6ビット結合して8ビットとした値が量子化代表値となる。なお、c=4の場合(即ち、3ビットの出力コードのMSBが“1”となる場合)では、8ビット全てを“1”にした値が量子化代表値となる。
【0052】
負号付加回路4はノイズ生成回路2から発生するノイズに負号及びオフセットを付加する。
【0053】
具体的には、ノイズ生成回路2の出力を(rnd)、比較器3の1ビット信号出力(これは入力画像データ(x)のビット5(b5)の値である)を(s)とすると負号付加回路4の6ビット出力(n)は以下の通りである。
【0054】
即ち、s=0であれば、n=rnd+32、
s=1であれば、n=32−rndである。
【0055】
ここで、s=0の時とs=1の絶対値が一致する必然性はない。つまり、s=0の時は32から63のノイズが、s=1の時は0から31のノイズが出力されれば良い。従って、比較器3の出力(s)を反転し、ノイズ生成回路2の出力(rnd)の上位に結合するだけで良い。このノイズは入力画像データを閾値側に引きつける(言い換えると、閾値をまたぎやすくする)ように作用するだけなので、本来選択される可能性のないコード(量子化代表値)が選択されることはほとんどない。
【0056】
例えば、入力画像データ(x)が“95”であった場合、補正値(e)の値によって、以下に述べるように、出力コード(c)として“1”または“2”が出力される。
【0057】
x=95であるなら、b5=0、即ち、s=0であるので、ノイズを重畳した場合、正のノイズが重畳される。
【0058】
負の量子化誤差の絶対値が最大になるのは量子化代表値(r)がr=128となる場合で、(x+e+n)=128(ただし、32のオフセットがついている)の時である。上述のように、s=0であれば、32≦n≦63であるので、(x+e)の最小値は128−63=65となる。従って、負の量子化誤差の最大値は65−128=−63となる。
【0059】
一方、正の量子化誤差が最大になるのは量子化代表値(r)がr=64となる場合で。(x+e+n)=127(ただし、32のオフセットがついている)の時である。上述のように、s=0であれば、32≦n≦63であるので、(x+e)の最大値は127−32=95となる。従って、正の量子化誤差の最大値は95−64=31となる。
【0060】
なお、補正値(e)は拡散フィルタ10により周囲の量子化誤差がフィルタリングされるので上記量子化誤差の範囲より狭くなるが、この範囲を補正値として入力画像に加えても、64≦x+e+n≦158<192となり、出力コード(c)は1或いは2のいずれかのみとなる。
【0061】
量子化器6における量子化誤差の範囲はノイズの重畳により±63(量子化代表値間の距離―1)に拡大されるため、減算器8の出力ビットとしては7ビット必要になる。
【0062】
なお、図2に示す拡散係数の分母は16であるので、拡散フィルタ10での補正値を表現するには整数部7ビット、小数部4ビットの合計11ビットが必要になるが、整数部7ビットが入力画像データに加算される。
【0063】
この実施形態におけるノイズ生成回路2は上記範囲の乱数を生成する回路で、例えば、LFSR(Linear Feedback Shift Register)を用いたM系列擬似ランダム符号生成回路を用いることにする。
【0064】
図3は9ビットのM系列擬似ランダム符号生成回路の構成を示す図である。図3において、11はシフトレジスタ、12はEXOR回路である。
【0065】
このようなM系列擬似ランダム符号生成回路では、シフトレジスタ11にまず“0”以外の値をロードしておき、シフトレジスタ11の原始多項式に対応するビットをEXOR(排他的論理和)を演算した結果をLSB(ビット0)に入力し、左シフトすることで次の乱数を得る。
【0066】
この実施形態では、シフトレジスタ11の段数が9段なので、9次の原始多項式X9+X4+1を用いている。よって、シフトレジスタ11のMSB(ビット8)とビット3とをEXOR回路12においてEXORを演算した結果をLSB(ビット0)の入力とし、次のクロックで左シフトする。このような構成により1〜511の乱数が、511周期で生成される。一方、ノイズ生成回路2のノイズ生成レンジは0〜31であるので、シフトレジスタ11の任意の5ビットを用いることで、511周期で0〜31の値を発生するノイズ生成回路2を構成することができる。
【0067】
上述のように、M系列擬似ランダム符号はシフトレジスタを用いているので、連続する乱数には若干相関性がある。このため、比較的低周波にノイズパワーのピークがあり、この値をノイズとして入力画像に重畳してしまうと、出力画像にノイズが観測されやすい。一方、この実施形態のように誤差拡散のフィードバックループ内にノイズを重畳した場合は、たとえDCオフセットがあるノイズであってもフィードバックループにより濃度は保存される。また、上述のように比較的低周波にノイズパワーのピークがあっても誤差拡散のフィードバックループにより、量子化誤差(ノイズ)の低周波成分は視覚的に検知されにくい高周波に変調されるので、出力画像におけるノイズの影響はかなり小さくなる。
【0068】
従って以上説明した実施形態に従えば、ノイズを入力画像にではなく、補正後の入力画像に加算する構成をとっているため、いかなるノイズを加えても、誤差拡散のフィードバックループにより濃度は保存される。また、誤差拡散のフィードバックループにより、ノイズの低周波成分は視覚的に検知されにくい高周波に変調されるので、出力画像におけるノイズの影響を抑えることができる。
【0069】
また、入力画像データに加えるノイズの絶対値を量子化器の閾値と量子化代表値との差分値未満(即ち、32未満)に制限し、さらにノイズを加算する方向を閾値方向に合わせているので、最小のノイズ振幅で閾値をまたがせる(出力コードの位置をずらす)ことが可能となる。
【0070】
さらに、ノイズを重畳しても一定入力値に対する出力コードは各々の出力コードの配置がずれるだけで、連続する2つのコード以外のコードはほとんど出力されず、コードの比率も変わらない。よって、ノイズ重畳による粒状感の増加等の弊害は殆どなく“テクスチャ”、“ワーム”、“はき寄せ”等の画像劣化を引き起こす要因を排除することができる。
【0071】
<第2実施形態>
図4は本発明の第2実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。図4において、13はリミッタである。以下、第1の実施形態との違いを中心に説明する。
【0072】
この実施形態ではノイズ生成回路2と負号付加回路4の間にリミッタ13が挿入されている。リミッタ13はノイズ生成回路2の出力を入力画像データと量子化における閾値との差分に応じた値に制限して負号付加回路4に入力する。
【0073】
次に、具体的な動作を説明する。
【0074】
この実施形態はノイズの重畳により閾値をまたぐように補正し、反対側のコードを生成しやすくすることを特徴としているので、ノイズ生成回路2の最大値は、入力画像データと閾値との差分値程度で充分である。従って、入力画像データと閾値との差分によってリミッタの制限値を制御することにより、より最適なノイズの生成が可能である。
【0075】
具体的には、入力画像データのビット5(b5)が“1”なら、即ち、s=1なら、入力画像データ(x)の下位5ビットの値(b4〜b0(LSB))によって表現される値)または所定値のいずれか小さい方を制限値(LMT)とし、ノイズ生成回路2の出力がこの制限を越える値になった時は、この制限値を出力する。一方、入力画像データのビット5(b5)が“0”なら、即ち、s=なら、入力画像データの下位5ビット(b4〜b0)を反転した値または所定値のいずれか小さい方を制限値(LMT)とし、ノイズ生成回路2の出力がこの制限を越える値になった時は、この制限値を出力する。
【0076】
なお、このようなアルゴリズムでは入力画像データ(x)の値が閾値付近であるときは、制限値が“0”付近の非常に小さな値となり、重畳するノイズが小さいものとなって、ノイズを重畳するという効果が薄れてしまうので、この制限値が設定された最小値以下の場合は、その設定された最小値を制限値とする。
【0077】
ノイズ生成回路2が発生する最大値を量子化器6の閾値と量子化代表値との差分付近に設定した場合、即ち、この実施形態のように、その最大値“31”が差分値“32”に近い場合、その最大値が重畳された画素はほんのわずかな値の補正値(e)で入力画像データが表す値よりはるかに遠い量子化代表値が選ばれてしまう。このため、エッジ部分では入力画像データより遠いコードがある間隔で並んだり、例えば、画像形成をインクジェット記録ヘッドを走査して行う場合、その記録ヘッドの走査方向にわずかなグラデーションがかかっているような画像で量子化代表値付近の値をもつ画像入力データに対し、本来生成されるはずのないコードが生成されてしまうという弊害が生じる場合がある(これは過補償と呼ばれる)。
【0078】
しかしながら、この実施形態によれば、リミッタ13で発生ノイズの最大値を適切な値に制限することにより、たとえ補正値(e)を考慮しても、量子化代表値が入力画像データの値に近接しない値になってしまう確率は非常に少なくなるので、より弊害なく“テクスチャ”、“ワーム”、“はき寄せ”の発生を抑え、形成画像の画質劣化を防止することができる。
【0079】
なお、この実施形態では制限値を入力画像データの下位ビットの値から直接生成したが、入力画像データの下位ビットを用いてテーブルを参照することにより同様の制限値を決定するようにしても良い。
【0080】
<第3実施形態>
図5は本発明の第3実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。図5において、14はゲイン調整回路である。
【0081】
以下の説明では、第1〜第2実施形態との違いを中心に説明する。
【0082】
この実施形態ではノイズ生成回路2と負号付加回路4の間にゲイン調整回路14が挿入されている。ゲイン調整回路14はノイズ生成回路2の出力を入力画像データと量子化における閾値との差分に応じた値に制限するため、ノイズ生成回路2からの出力に所定値を乗じた値を負号付加回路4に入力する。
【0083】
次に、具体的な動作を説明する。
【0084】
第2実施形態でも説明したように、ノイズ生成回路2の最大値は、入力画像データと閾値との差分値程度で充分である。従って、入力画像データと閾値との差分によってゲイン調整回路14のゲインを制御することにより、より最適なノイズを生成する。
【0085】
具体的には、入力画像データのビット5(b5)が“1”なら、即ち、s=1なら、入力画像データ(x)の下位5ビットの値(b4〜b0(LSB))によって表現される値)または所定値のいずれか小さい方をゲインとしてノイズ生成回路2の出力に乗じ、その結果得られた信号の上位ビットをノイズとして負号付加回路4に出力する。一方、一方、入力画像データのビット5(b5)が“0”なら、即ち、s=なら、入力画像データの下位5ビット(b4〜b0)を反転した値または所定値のいずれか小さい方をゲインとしてノイズ生成回路2の出力に乗じ、その結果得られた信号の上位ビットをノイズとして負号付加回路4に出力する。
【0086】
なお、第2の実施形態と同様、このようなアルゴリズムでは入力画像データ(x)の値が閾値付近であるときは、制限値が“0”付近の非常に小さな値となり、重畳するノイズが小さいものとなって、ノイズを重畳するという効果が薄れてしまうので、そのゲインが設定された最小値以下の場合は、その設定された最小値をゲインとする。
【0087】
従って以上説明した実施形態に従えば、ノイズのゲインを調整することにより、所定値以下のノイズの発生確率が一様になり、第2実施形態と同様、より弊害の少ないノイズの重畳が可能となる。
【0088】
なお、この実施形態ではゲインを入力画像データの下位ビットの値から直接生成したが、入力画像データの下位ビットを用いてテーブルを参照することにより同様のゲインを決定するようにしても良い。
【0089】
また、この実施形態の説明ではゲイン調整回路14として乗算器を用いた構成としたが、ゲイン調整できるものであれば、乗算器でなくとも良い。例えば、ビットシフトを用いる構成や、上位ビットをマスクする構成等でも同様の効果が得られる。
【0090】
特に、ノイズのビットをマスクする構成の場合、入力画像データのビット5(b5)が“1”なら、即ち、s=1なら、入力画像データの下位5ビット(b4〜b0(LSB))をマスクデータとしてノイズと論理積(AND)を演算する一方、入力画像データのビット5(b5)が“0”なら、即ち、s=0なら、入力画像データの下位5ビット(b4〜b0(LSB))を反転した値をマスクデータとしてノイズとAND演算をすると良い。
【0091】
またこの時、下位ビットの出力を固定パターン(例えば、“011”)とすることにより、設定最小値を簡単に得ることができる。
【0092】
<画像入出力システム>
ここでは上述したいくつかの実施形態に従うm値化処理装置が応用される画像入出力システムについて説明する。ここで、2≦m<nであり、nは入力多値画像データの各色成分の画像データを表すビット数である。以上説明した実施形態ではm=5として説明した。
【0093】
図6は画像入出力システムの一般的な構成を示すブロック図である。
【0094】
カラースキャナ、デジタルカメラなどの入力装置1000で読み取った或いは生成された多値カラー画像データは、通常パーソナルコンピュータ1010を経て、CRT、PDP、或いはLCDなどのモニタ1020に出力されその画面にカラー画像が表示される。
【0095】
ユーザは、そのモニタ画面を見ながら、パーソナルコンピュータ1010に接続されたキーボード1030、スイッチ(SW)1040、ポインティングデバイス1050などを操作しながらレイアウト修正やRGB各色成分について色修正を行う。
【0096】
このようにして、修正編集されたカラー画像データはパーソナルコンピュータ1010において、これから出力する出力装置(例えば、インクジェットプリンタ)1060で用いる記録媒体(例えば、記録用紙)の種類、記録モードに合わせて輝度濃度変換が実行され、RGBデータからCcMmYK成分からなる多値濃度データに変換される。
【0097】
このCcMmYK各色成分の多値濃度データは、前述したm値化処理部に入力されてm値化処理が行われ、CMYK各色成分m値化データが生成され、これらが出力装置1060(例えば、インクジェットプリンタ)に転送されてカラー画像が記録される。
【0098】
なお、図6ではm値化処理装置がパーソナルコンピュータ1010に組み込まれた例を示しているが、その代わりに出力装置1060の側に組み込まれていても良い。
【0099】
以上説明した実施形態では、誤差拡散処理が論理回路を用いて実行されるようなものとして説明したが、本発明はこれによって限定されるものではない。例えば、近年におけるマイクロプロセッサ(MPU)の高性能化に伴って、誤差拡散処理をROMに格納されたプログラムをそのマイクロプロセッサが読みだして実行するようにして実現しても良い。
【0100】
この場合、そのマイクロプロセッサがパーソナルコンピュータ1010に内蔵されている構成でも良いし、デジタルカメラを直接接続してカラー画像を出力する出力装置(例えば、インクジェットプリンタ1060)の場合にはその出力装置にマイクロプロセッサが内蔵されるような構成でも良い。
【0101】
つまり本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体(または記録媒体)を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ(またはCPUやMPU)が記憶媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本発明を構成することになる。また、コンピュータが読み出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム(OS)などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0102】
さらに、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるCPUなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0103】
【発明の効果】
以上説明したように本発明によれば、“テクスチャ”、“ワーム”、“はき寄せ”等、誤差拡散によって発生する画質を劣化させる要因を抑え、画像に重畳されるノイズの影響を最小化し、高品位な画像を形成することができるという効果がある。
【図面の簡単な説明】
【図1】本発明の第1実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。
【図2】拡散フィルタ9の演算に用いる量子化誤差の位置と拡散係数の例を示す図である。
【図3】ノイズ生成回路2の構成を示す図である。
【図4】本発明の第2実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。
【図5】本発明の第3実施形態に従うm値化処理方法を用いたm値化処理部の構成を示すブロック図である。
【図6】画像入出力システムの一般的な構成を示すブロック図である。
【符号の説明】
1 加算器
2 ノイズ生成回路
3 比較器
4 負号付加回路
5 加算器
6 量子化器
7 逆量子化器
8 減算器
9 誤差バッファ
10 拡散フィルタ
11 シフトレジスタ
12 EXOR回路
13 リミッタ
14 ゲイン調整回路

Claims (10)

  1. n値画像データに誤差拡散処理を施して擬似中間調画像を形成するためのm値画像データ(2≦m<n)に変換する画像処理方法であって、
    n値画像データを入力する入力工程と、
    前記入力されたn値画像データに以前の誤差拡散によって累積した量子化誤差を加算する第1加算工程と、
    ノイズの絶対値を発生する発生工程と、
    前記発生工程において発生するノイズの値を予め定められた値以下に制限する制限工程と、
    前記入力されたn値画像データの特定ビットの値によって前記ノイズの符号を決定する決定工程と、
    前記決定工程において符号が決定されたノイズを前記第1加算工程において量子化誤差が加算されたn値画像データにさらに加算する第2加算工程と、
    前記第2加算工程においてノイズが加算されたn値画像データをm値化する量子化工程と、
    前記量子化工程においてm値化されたm値データを量子化代表値に逆量子化する逆量子化工程と、
    前記第1加算工程において量子化誤差が加算されたn値画像データから前記逆量子化工程において得られた量子化代表値を減算して量子化誤差を生成する量子化誤差生成工程と、
    前記量子化誤差生成工程において生成した量子化誤差をバッファに格納する格納工程とを備えることを特徴とする画像処理方法。
  2. 前記第1加算工程ではバッファに格納された量子化誤差を特定の分配係数によって複数の画素に分配して累積された量子化誤差を加算することを特徴とする請求項1に記載の画像処理方法。
  3. 前記発生工程では、発生するノイズのビット数より長いシフトレジスタを用いたM系列擬似ランダム符号生成回路のシフトレジスタ出力より前記ノイズのビット数相当分を選択することにより、前記ノイズを生成することを特徴とする請求項1に記載の画像処理方法。
  4. 前記制限工程では、特定ビットが1であるなら、入力されたn値画像データの前記特定ビットより下位にあるビットによって表現される値、或いは、第1の予め定められた値の内、いずれか小さい方を前記予め定められた値とする一方、特定ビットが0であるなら、入力されたn値画像データの前記特定ビットより下位にあるビットを反転することにより得られるビットによって表現される値、或いは、前記第1の予め定められた値の内、いずれか小さい方を前記予め定められた値とすることを特徴とする請求項に記載の画像処理方法。
  5. 前記制限工程によって制限された値が第2の予め定められた値以下の場合、前記第2の予め定められた値を前記予め定められた値とすることを特徴とする請求項に記載の画像処理方法。
  6. 前記制限工程は、リミッタによって制限が実行されることを特徴とする請求項に記載の画像処理方法。
  7. 前記制限工程は、ゲイン調整器によって制限が実行されることを特徴とする請求項に記載の画像処理方法。
  8. n値画像データに誤差拡散処理を施して擬似中間調画像を形成するためのm値画像データ(2≦m<n)に変換する画像処理装置であって、
    n値画像データを入力する入力手段と、
    前記入力手段によって入力されたn値画像データに以前の誤差拡散によって累積した量子化誤差を加算する第1加算手段と、
    ノイズの絶対値を発生する発生手段と、
    前記発生手段により発生するノイズの値を予め定められた値以下に制限する制限手段と
    前記入力されたn値画像データの特定ビットの値によって前記ノイズの符号を決定する決定手段と、
    前記決定手段によって符号が決定されたノイズを前記第1加算手段によって量子化誤差が加算されたn値画像データにさらに加算する第2加算手段と、
    前記第2加算手段によってノイズが加算されたn値画像データをm値化する量子化手段と、
    前記量子化手段によってm値化されたm値データを量子化代表値に逆量子化する逆量子化手段と、
    前記第1加算手段によって量子化誤差が加算されたn値画像データから前記逆量子化手段によって得られた量子化代表値を減算して量子化誤差を生成する量子化誤差生成手段と、
    前記量子化誤差生成手段によって生成された量子化誤差を格納するバッファとを備えることを特徴とする画像処理装置。
  9. 請求項1乃至のいずれかに記載の画像処理方法における各工程をコンピュータによって実行するプログラム。
  10. 請求項に記載のプログラムをコンピュータで実行可能なコードによって記述し、該コードを格納した記憶媒体。
JP2003016713A 2002-07-11 2003-01-24 画像処理方法及び装置 Expired - Fee Related JP4086669B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2003016713A JP4086669B2 (ja) 2003-01-24 2003-01-24 画像処理方法及び装置
US10/611,948 US8009327B2 (en) 2002-07-11 2003-07-03 Method and apparatus for image processing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003016713A JP4086669B2 (ja) 2003-01-24 2003-01-24 画像処理方法及び装置

Publications (3)

Publication Number Publication Date
JP2004229133A JP2004229133A (ja) 2004-08-12
JP2004229133A5 JP2004229133A5 (ja) 2006-03-09
JP4086669B2 true JP4086669B2 (ja) 2008-05-14

Family

ID=32904082

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003016713A Expired - Fee Related JP4086669B2 (ja) 2002-07-11 2003-01-24 画像処理方法及び装置

Country Status (1)

Country Link
JP (1) JP4086669B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5309068B2 (ja) * 2010-03-30 2013-10-09 京セラドキュメントソリューションズ株式会社 画像処理装置および画像形成装置

Also Published As

Publication number Publication date
JP2004229133A (ja) 2004-08-12

Similar Documents

Publication Publication Date Title
US5353127A (en) Method for quantization gray level pixel data with extended distribution set
JP4217657B2 (ja) 画像処理方法、プログラム、記憶媒体及び装置
US8009327B2 (en) Method and apparatus for image processing
JPH08237483A (ja) イメージデータを処理するためのシステム及び方法
JP5534787B2 (ja) 画像処理装置及び画像処理方法
JP4149269B2 (ja) 中間ドットを用いた誤差拡散方法及びシステム
JP3949931B2 (ja) 画像形成装置
JPH10112801A (ja) 入力イメージ処理方法
JP2002077623A (ja) 画像処理装置
JP4248654B2 (ja) 出力装置へ出力する文書イメージを準備する処理装置
JP4086669B2 (ja) 画像処理方法及び装置
JP3821327B2 (ja) 画像処理方法および装置
JP5504740B2 (ja) 画像処理装置
JP5312158B2 (ja) 画像処理装置及び画像処理方法
JP5471135B2 (ja) 画像処理装置
JP2008306400A (ja) 画像処理方法、画像処理装置、画像形成装置、コンピュータプログラム及び記録媒体
JP3124589B2 (ja) 画像処理装置
JP2002344742A (ja) 画像形成装置
JP3560265B2 (ja) 疑似中間調処理方法
JP2005051337A (ja) 画像処理方法、画像処理装置およびプログラム
KR100400012B1 (ko) 고속 및 고 해상도의 이미지 변환 장치 및 방법
JP3741829B2 (ja) カラー画像処理装置
JP3428822B2 (ja) 画像処理装置及び画像処理方法
JP3788913B2 (ja) 画像処理方法、画像処理装置および画像形成装置
JPH10136205A (ja) 画像処理装置

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060119

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060119

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20070723

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070727

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070925

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: 20080204

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080219

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110228

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4086669

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: 20120229

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130228

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140228

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees