JPH11234513A

JPH11234513A - 擬似階調画像処理装置

Info

Publication number: JPH11234513A
Application number: JP10032881A
Authority: JP
Inventors: Nobuhito Matsushiro; 信人松代
Original assignee: Oki Data Corp
Current assignee: Oki Electric Industry Co Ltd
Priority date: 1998-02-16
Filing date: 1998-02-16
Publication date: 1999-08-27
Anticipated expiration: 2018-02-16
Also published as: US6459817B1; JP3698541B2

Abstract

(57)【要約】（修正有）【課題】プリンタ等の出力装置の特性に合わせた、画
素配置の最適化を行うことができる画像階調画像処理送
致を提供する。【解決手段】閾値設定部により、周期の長さが互いに
異なる複数の周期信号に基づいて閾値の値を変動させ、
擬似階調画像処理装置で、上記変動させた閾値を用いて
擬似階調画像処理を行う。また、上記周期信号が水平走
査線に沿う、一次元的なものである。また、上記周期信
号が、二次元的なものである。また、上記閾値設定部
は、原画像信号の座標値に対応する上記複数の周期信号
のそれぞれの信号値に、原画像階調値に依存する重み値
を掛けて加算し、該加算結果によりｎ値化（ｎは２以上
の自然数）のための閾値を修正する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、従来原画像を表
わす多値画像入力信号を受け、それよりも低い階調（例
えば２値）の信号を用い、特定の階調の画素の（例えば
黒画素の）密度を変化させることにより、元の多値画像
と同じ濃淡を感じさせるようにする擬似階調画像処理装
置に関する。

【０００２】

【従来の技術】従来、上記のごとき擬似階調画像処理方
法として、誤差拡散法やディザ法がある。

【０００３】誤差拡散法は、量子化する入力画素の階調
値と出力画素の階調値との差を誤差とし、量子化する画
素の周辺の画素の誤差に重み付けを施した後に量子化す
る画素の階調値に加算し、その加算値を固定閾値で量子
化する方法であり、階調変換の際に発生する量子化誤差
による階調差を補正するのに効果的である。出力画素が
２値で表わされ、拡散法を一次元的に行なう場合を図１
９に示す。

【０００４】入力画素の階調Ｉに対する処理の順序は、
図の左から右である。先ず、図１９（ｂ）に示す第１の
画素Ｐ１は、固定閾値Ｔｈより大きいので、白（以後シ
ンボル１と呼ぶ）と出力される。ここで、白は、図１９
（ａ）に示す線の値（階調値の最大値）であり、これと
入力値との差（入力値−出力値）が誤差Ｅ１（負の値）
となる。この誤差Ｅ１に、該誤差を拡散させる画素が１
つなので、重み値１を掛けて右隣りの第２の画素Ｐ２に
加算する。第２の画素Ｐ２は、前記誤差Ｅ１を加算して
修正した後の値（修正階段調値）を固定閾値Ｔｈで２値
化する。ここでの出力は黒（以後シンボル０と呼ぶ）で
あり、生じた誤差Ｅ２（正の値）を隣の第３の画素Ｐ３
に加算する。このようにして、順次誤差Ｅ１，Ｅ２，Ｅ
３，．．．を加算して固定閾値Ｔｈで２値化するのが一
次元の誤差拡散法である。

【０００５】図２０は、二次元の誤差拡散法の説明図で
あり、主走査方向ｘ及び副走査方向ｙとして、注目画素
Ｐ１に対して誤差を拡散させる近傍画素Ｐ1A，Ｐ1B，Ｐ
1C，Ｐ1Dと、その近傍画素Ｐ1A，Ｐ1B，Ｐ1C，Ｐ1Dに対
する重み値７／１６，３／１６，５／１６,１／１６が
示されている。

【０００６】この二次元の誤差拡散法では、一次元の誤
差拡散法と同様に（入力値−出力値）の誤差を求め、こ
の求めた誤差に、所定の重み値を掛けて走査方向ｘ，ｙ
の画素の階調値に拡散して加算し、その加算値を固定閾
値Ｔｈで２値化している。いま２値化しようとしている
着目画素Ｐ２からみると、既に処理済みの近傍４画素Ｐ
2A，Ｐ2B，Ｐ2C，Ｐ2D から図２１に示すように拡散誤
差が集められることと等価となる。

【０００７】

【発明が解決しようとする課題】しかしながら前記従来
方法では、プリンタ等の出力装置の特性に合わせた、画
素配置の最適化を行うことができなかった。

【０００８】即ち、変換された出力データの個々の値
は、量子化誤差の値によってランダムに決定されるの
で、ある階調ではある特定の並びの出力データパターン
にするのが好ましい場合であっても、そのような操作は
不可能であった。例えば、階調が５０％の近傍のグレー
の入力であった場合、チェッカーフラッグ状のパターン
で擬似階調を表すのが好ましいが、そうはならない。

【０００９】

【課題を解決するための手段】本発明の擬似階調画像処
理装置は、周期の長さが互いに異なる複数の周期信号に
基づいて、閾値の値を変動させる閾値設定部（３０）
と、該変動させた閾値を用いて擬似階調画像処理を行う
擬似階調画像処理部（１０）とを有するものである。

【００１０】本発明を用いれば、モノクロプリンタやカ
ラープリンタで出力できる２値ないし数階調程度の階調
数で、擬似的に階調を再現することができる。

【００１１】

【発明の実施の形態】最初に本発明の基本概念を説明す
る。

【００１２】まず、簡単のために一次元の場合について
説明する。図１の上部には、一次元的に１６個の画素が
配置されている。この階調値が最大値に対して１２／１
６で、１６個の画素の内４個を黒、残りの１２個を白に
したい場合を考える。この場合、図１の下部に示すよう
な周波数が異なる２種の周期信号を固定閾値Ｔｈに重畳
して２値化すれば、図１の上部に示した画素配置が得ら
れる。このように任意の画素配置制御は、白画素と黒画
素の比率（階調値）に依存する形で、周期が互いに異な
る周期信号を重ね合わせ、これにより固定閾値（Ｔｈ＝
Imax/2：Imax：階調最大値）を周期的に修正することで
行うことができる。

【００１３】二次元的に処理する場合には、上記一次元
の周期波形の代わりに図２（ａ）〜（ｄ）に示す閾値パ
ターンマトリクスを用い、これを繰返して用いることに
より、画面に周期性を持たせる。

【００１４】図２（ａ）〜（ｄ）に示すパターンマトリ
クスＰＭ１〜ＰＭ４は、２×２、４×４、８×８、１６
×１６の４種類からなる。このパターンマトリクスを画
面全体に繰返し、画面全体をこのパターンマトリクスで
埋め尽くすことにより、２次元の周期性を持たせる。な
お、図３は、パターンマトリクスＰＭ１を画面全体で繰
返す場合の配列を示し、図４は、パターンマトリクスＰ
Ｍ２を画面全体で繰返す場合の配列を示す。図３及び図
４において、実線がパターンマトリスクの境界を示し、
点線は一つのパターンマトリクス内の画素の境界を示
す。

【００１５】２×２のパターンマトリクスＰＭ１は周波
数が高い最も周期パターンを構成し、１６×１６のパタ
ーンマトリクスＰＭ４は周波数が低い周期パターンを構
成する。パターンマトリクスＰＭ１〜ＰＭ４は、各々第
１ないし第４象限において異なる値（信号値）を有し、
これを元に各画素に対する閾値が修正される。

【００１６】本発明においては、入力信号の階調値は、
最小値（黒色）が０、最大（白色）が２５５であるもの
とする。

【００１７】この場合、誤差拡散法により２値化する場
合の閾値は図５(a)から(c)に示すように、基本的には２
５５／２＝１２８となる。

【００１８】本発明においては、この閾値を下記複数の
テーブル（パターンマトリクスの信号値を記憶するも
の）に記憶された値（信号値）に基づき定められた値
（閾値調整値）を重畳して図５(a)又は(b)のように変化
させることにより入力に合わせた規則的な閾値変化を行
う。

【００１９】これら各信号値は、最大値を１２７、最小
値を−１２８とすることにより、基本的な閾値１２８を
０〜２５５の間で変化させることができる。又、最大値
及び最小値を１２７、−１２８として割り付けたことに
より、その間を３等分する値として図５(c)に示すよう
に、４２（＝１２７−２５６／３）と−４３（１２７−
２×２５６／３）を選んだ。この４つの値を第１ないし
第４象限に割当て、この４つの値を元にして、入力信号
に対する各象限の信号値を定めている。閾値修正に当っ
ては、上記の４つの値に対し、入力画素値に応じたテー
ブル毎の重み値（重み付け係数）Ｗが乗じられ、それぞ
れのテーブルについて得られた値を加算し、さらに、加
算結果に対して修正の度合いを調整するベース係数α
（例えば、０．５程度の値が好ましい）が乗じられ、こ
れにより、基本閾値１２８に近い範囲で調整される。

【００２０】次に、各パターンマトリクスについて説明
する。パターンマトリクス２を例にとると、これは４×
４のマトリクスで、左上の象限の２×２の４つの画素に
対する信号値は全て−１２８の値であり、同様に右上の
象限の２×２の４つの画素に対する信号値は全て１２７
の値を取る。他の象限の信号値も同様に全て同じ値を取
る。他のパターンマトリクスの信号値も４つに分けた象
限毎に同じ値を取るように構成されている。

【００２１】尚、印刷する画像の総ドット数が大きい場
合は、３２×３２のマトリクス、６４×６４のマトリク
スと更にマトリクスの種類を増加していくと、より精度
が向上する。

【００２２】各テーブルにおける横方向に隣接する象限
同士で信号値を正負逆とし、絶対値を略同じ（−１２８
と１２７、４２とー４３）とする理由はフーリエ級数展
開に用いる周期関数として例えば（正の期間と負の期間
が同じ）正弦波を用いるのと同様に、＋の範囲（周期で
いえば正の期間に相当）と−の範囲（負の期間）を同じ
にする為である。又、４つの象限の信号値の合計をほぼ
０とすることにより、擬似階調化することによる画像濃
度の変化を避けることができる。逆に言えば、この信号
値の総和が０以下であれば、変換後の画像は白っぽくな
り、０以上であれば黒っぽくなる。

【００２３】尚、本発明において−１２８、１２７、４
２、−４３の値を用いたが、これを−１２８，１２７，
１２７，−１２８とすれば、横方向のみならず、縦方向
のフーリエ展開も同時に可能となる。

【００２４】重み値は、上記のように、各パターンマト
リクスの信号値に掛けられるもので、各テーブルに対応
した重み値をＷ１〜Ｗ４とすると、その値は、入力画素
値に対して例えば図６に示すように変化する。各所定の
入力画素値範囲でゼロ以外の値を取り、該所定の範囲以
外ではゼロである。また、上記所定の入力画素値範囲内
においては例えば図示のように正弦波の正側の半波の形
状である。そして、より小さな周期のパターンマトリク
スに対応する重み値の上記所定入力画素値範囲（ゼロ以
外の値となる範囲）が、より高い入力画素値の領域に設
定されている。また、パターンマトリクスを大きさ順を
並べたとき、隣合うパターンマトリクスに対応した重み
値の上記所定範囲は互に重なり合っている。より具体的
には、最も小さいパターンマトリクスに対応した重み値
Ｗ１は、上記所定範囲が１７０から２５５であり、２５
５においてピークとなる正弦半波（の左半分）の形状で
ある。２番目に小さいテーブルに対応した重み値Ｗ２
は、上記所定範囲が８５から２５５であり、１７０にお
いてピークとなる正弦半波の形状である。３番目に小さ
いパターンマトリクスに対応した重み値Ｗ３は、上記所
定範囲が０から１７０であり、８５においてピークとな
る正弦半波の形状である。最も大きいパターンマトリク
スに対応した重み値Ｗ４は、上記所定範囲が０から８５
であり、０においてピークとなる正弦半波（の右半分）
の形状である。重み値曲線のピークはすべて１であり、
互いに重なりあう範囲において、一方が次第に小さくな
るに連れ、他方が次第に大きくなっており、重なり合う
重み値の和が１に近い値になるように設定されている。
（位相が９０度異なる正弦波を重ね合わせると、１〜
１．４の範囲で変動するが、ここでは１．４も１に近い
値と見る。）このように、入力画素の画面内の位置を、図２に示す基
本周期パターンマトリクスに対応付けて、その画素の位
置に応じて信号値を読み出すとともに、入力画素値によ
って、重み値Ｗ１〜Ｗ４を決め、上記閾値調整値の元と
なる値（パターンマトリクスの信号値）にそれぞれ対応
する重み値を乗算し、乗算結果を加算することでフーリ
エ積分が実現できる。

【００２５】第１の実施の形態図７は、本発明の第１の実施の形態の画像処理装置を示
すブロック図である。図示のようにこの画像処理装置
は、入力画像メモリ２と、誤差拡散部１０と、閾値設定
部３０と、制御部１０と、アドレス生成部４とを有す
る。

【００２６】入力画像メモリ２は、外部から入力される
画素信号（０から２５５までの値を取る、２５６階調値
で表わされる）を受け、各々画面内の位置に対応したア
ドレス（ｘ’，ｙ’で表わされる）に記憶する。

【００２７】アドレス生成部４は、後述の制御部１０か
ら、画面上の位置を表わすアドレスないし座標値ｘ，ｙ
を受け、入力画像メモリ２への書込み及び読み出しのた
めのアドレス（ｘ’，ｙ’）を生成する。

【００２８】画像メモリ２に蓄えられた画素値は順に読
み出され、誤差拡散部１０、閾値設定部３０による処理
を受ける。読み出しの順序はラスタスキャンと同様、画
面上の水平走査線に沿い、上から下へと読み出される。
画面上の位置に対応したメモリ内の位置から順に読み出
される。

【００２９】誤差拡散部１０は、誤差拡散法により多値
階調データを２値階調データに変換する回路ブロックで
ある。

【００３０】閾値設定部３０は、誤差拡散部１０で用い
られる閾値を設定するである。

【００３１】一画面分の処理が終わると、画像メモリ２
の内容が書換えられる。書換えの時間中処理が中断する
のを防ぐため、２画面分の容量を持つ画像メモリを用意
し、２つの画面に対応する領域を交互に使うようにして
も良い。

【００３２】誤差拡散部１０は、加算部１１と、比較部
１２と、０／１処理部１３と、画像出力部１４と、減算
部１５と、誤差メモリ１６と、フィルタ１７とを有す
る。

【００３３】加算部１１は、画像メモリ２からの出力
（入力原画像の画素の階調値）にフィルタ１７からの出
力を加算して修正階調値を発生し、これを比較部１２に
送る。

【００３４】比較部１２は、加算部１１からの出力Ｉ’
（修正された階調値）を、閾値設定部３０からの出力
（修正された閾値）Ｔｈ’と比較し、前者が後者以上で
あれば出力をＩmax＝２５５とし、そうでなければ出力
をＩmin＝０とする。

【００３５】０／１処理部１３は、比較部１２の出力が
Ｉmax＝２５５のときはその出力を１（シンボル１）と
し、比較部１２の出力がＩmin＝０のときはその出力を
０（シンボル０）にする。

【００３６】画像出力部１４は、０／１処理部１３の出
力を、画面内の座標値を表わす情報（ｘ，ｙ）に対応づ
けて出力する。

【００３７】減算部１５は、加算部１１の出力（修正階
調値）と比較部１２の出力（＊変形例では、２値−２５
６値変換部４’の出力）の差（２値化に伴う誤差）を求
める。

【００３８】即ち、誤差＝（修正階調値）−（２値化階調値）で求められる。２値化階調値はＩmax＝２５５またはＩm
in＝０の値を取る。

【００３９】誤差メモリ１６は減算部１５で計算された
誤差（２値化によって発生した誤差）をその時の画素位
置に対応するメモリ位置に記憶する。このため、誤差メ
モリ１６は、２ライン分の記憶容量を有し、現在処理さ
れている画素と同じ水平走査線上の、それより左のすべ
ての画素、及び現在処理されている画素の属する水平走
査線より一つ上の水平走査線上のすべての画素について
誤差を記憶している。その誤差の計算法については後で
詳細に説明する。

【００４０】フィルタ１７は、誤差メモリ１６に記憶さ
れている誤差を元に、現に画像メモリ２から出力されて
いる画素（着目画素）に対して適用すべき誤差を計算し
て出力する。この誤差の計算に当っては、着目画素を基
準とし、そのそれより前に処理された画素、具体的に
は、図２１に示すように着目画素の左隣の画素、着目画
素の真上の画素、その左隣及び右隣の画素についての誤
差を読み出し、これに近さの程度に応じた重み付け係数
をを掛けて加算し、加算結果を誤差として出力する。こ
のような処理を誤差フィルタリングと言う。

【００４１】図２１は、重み付け係数の一例をも示して
いる。右上の画素から重み付け係数が左上の画素からの
重み付け係数よりも大きいのは、処理の順の関係で、着
目画素の右側の誤差を考慮に入れることができないた
め、これを補償するためである。

【００４２】以上のようにして誤差拡散部１０は、誤差
拡散法により擬似階調処理を行なう。

【００４３】本発明の特徴は、誤差拡散部１０で用いる
修正閾値Ｔｈ’の作成法、言換えると、これを設定する
閾値設定部３０の構成にある。

【００４４】閾値構成部３０は、４つのテーブルＴ１、
Ｔ２、Ｔ３、Ｔ４と、重み付けテーブル３１と、乗算部
Ｍ１〜Ｍ４と、加算部３２と、乗算部３３と加算部３４
とを有する。

【００４５】テーブル部Ｔ１〜Ｔ４は、それぞれ、周期
の異なる二次元周期信号のパターンマトリクスＰＭ１〜
ＰＭ４を記憶している。信号値のアドレスは、(x,y)座
標値に対応付けられている。

【００４６】パターンマトリクスＰＭ１〜ＰＭ４は図３
及び図４に示すように、画面全体に繰返されているが、
画面全体に対応したアドレスを持つテーブルを用意する
代りに、一つのパターンマトリクスを記憶したテーブル
を用意し、アドレス変換により同じ部分を繰返しアクセ
スすることで所望の結果を得ることができる。

【００４７】パターンマトリクスの大きさ（画素数で表
わされた大きさ）が２のべきである場合、上記アドレス
変換は、例えば図８のように、アドレスの下位ビットの
みを用いることで実現できる。例えば、テーブルＴ１に
対しては最下位ビットのみでアクセスする。テーブルＴ
２に対しては、下位２ビットでアクセスする。テーブル
Ｔ３には下位３ビットを用いる。テーブルＴ４には下位
４ビットを用いる。これにより、画像メオリ２内の全画
素が順に読み出される間に、テーブルＴ１〜Ｔ４内の各
位置が繰返しアクセスされる。

【００４８】テーブルの他の例が図９に示されている。
この例では、各テーブルとも記憶位置は４つ（２×２）
ずつであり、一方ｘ，ｙアドレスも各１ビットずつしか
入力されない。即ち、テーブルＴ１には最下位ビットの
みが供給され、テーブルＴ２には下位から２番目のビッ
トのみが供給され、テーブルＴ３には下位から３番目の
ビットのみが供給され、テーブルＴ４には下位から４番
目のビットのみが供給される。このような構成であれ
ば、テーブル内の同じ位置が繰返しアクセスされる。例
えば、テーブルＴ４の場合、画像メモリ２の８つの画素
が出力される間、テーブルＴ４内の同じ位置が繰返しア
クセスされ、同じ値が出力される。

【００４９】重み付けテーブル３１は、画像メモリから
出力される入力画素値を入力とし、それに対応して重み
値Ｗ１〜Ｗ４を出力する。重み値Ｗ１〜Ｗ４はそれぞれ
テーブルＴ１〜Ｔ４に対応する。入力画素値と重み値Ｗ
１〜Ｗ４の関係は図６に示す通りである。

【００５０】この重み値テーブル３１は入力画素の階調
値により図６のグラフに対応するのいずれかの値を出力
する。尚、重み値として図６のような正弦波を用いる場
合はテーブルに代え計算により算出するようにしてもよ
い。

【００５１】乗算部Ｍ１〜Ｍ４はテーブルＴ１〜Ｔ４の
出力と、対応する重み値Ｗ１〜Ｗ４を乗算する。即ち、
乗算部Ｍ１はテーブルＴ１の出力と重み値Ｗ１を乗算す
る。乗算部Ｍ２はテーブルＴ２の出力と重み値Ｗ２を乗
算する。乗算部Ｍ３はテーブルＴ３の出力と重み値Ｗ３
を乗算する。乗算部Ｍ４はテーブルＴ４の出力と重み値
Ｗ４を乗算する。

【００５２】加算部３２は、乗算部Ｍ１〜Ｍ４の出力を
加算する。

【００５３】各テーブルの出力である信号値は各乗算部
により重み値テーブル３１の出力である係数Ｗと掛け合
わされた後加算されることにより、フーリエ展開され
る。

【００５４】乗算部３３は、加算部３２の出力にベース
係数αを掛ける。ベース係数αの値は制御部１０から与
えられる。ベース係数の値は０．５程度に定められる。
ベース係数αとの乗算により、閾値変化の範囲が調整さ
れる。

【００５５】加算部３４は、乗算部３３の出力ΔＴｈを
固定の閾値Ｔｈに加算する。閾値Ｔｈは１２８（０と２
５５の中間）に定められる。

【００５６】加算部３４の出力が修正された閾値Ｔｈ’
として誤差拡散部１０の比較部１２に入力される。

【００５７】重み値テーブル３１の値及びベース係数α
の値は、出力装置（例えばプリンタ）の特性に応じて設
定される。

【００５８】なお、重み値テーブルの関数を一例として
余弦波、正弦波としたが、出力装置の特性に合わせて半
円あるいは三角波を用いることもできる。更に、ベース
係数も出力装置の特性に応じて規則性を減じたい場合に
は０．４等、０．５より小さい値を用い、逆に規則性を
大きくしたい場合には０．６等、０．５より大きな値を
用いるとよい。

【００５９】次に、上記の装置の動作を図１０及び図１
１のフローチャートを用いて説明する。

【００６０】step 1：最初に、入力原画像の座標値、即
ち、画面上の座標値を初期化する。この初期化処理は制
御部１０が行う(step 1)。

【００６１】step 2 ：誤差メモリ１の内容を全て
「０」で初期化する。メモリの初期化処理も制御部１０
が行う。

【００６２】step 3 ：誤差メモリ１６の内容の内、着
目画素(図２１のＰ（ｘ，ｙ）の近傍４画素の位置に対
応する誤差値（図２１参照）を読み出し、重み付き和Ｅ(x,y)＝Ａ(x-1,y-1)・Ｅ(x-1,y-1)＋Ａ(x,y-1)・Ｅ
(x,y-1)＋Ａ(x+1,y-1)・Ｅ(x+1,y-1)＋Ａ(x-1,y)・Ｅ(x
-1,y) を求め、これをフィルタリング値Ｅ(x,y)とする。ここ
でＥ(x-1,y-1)，Ｅ(x,y-1)，Ｅ(x+1,y-1)，Ｅ(x-1,y)は
それぞれ着目画素の近傍の画素の誤差であり、Ａ(x-1,y
-1)，Ａ(x,y-1)，Ａ(x+1,y-1)，Ａ(x-1,y)は近傍の画素
に与えられる重み付け計数である。このフィルタリング
は、フィルタ１７が誤差メモリ１６を使用し行う処理で
ある。

【００６３】step 4 ：入力原画像の画素の階調値Ｉ(x,
y)にＥ(x,y)を加えた修正階調値Ｉ'(x,y)を求める。こ
の加算演算は加算部１１が行う。

【００６４】step 5 ：着目画素の座標値(x,y)から４個
のテーブルＴ１〜Ｔ４を参照し、(x,y)の値に依存した
テーブル値（パターンマトリクスＰＭの信号値）Ｋｊ
(x,y)とする（但し、ｊ＝１，２，３又は４）。

【００６５】step 6 ：４個のテーブル値Ｋ(x,y))にそ
れぞれ対応する重み値Ｗ１〜Ｗ４を掛け、その積の総和Ｋ１(x,y)・Ｗ１＋Ｋ２(x,y)・Ｗ２＋Ｋ３(x,y)・Ｗ
３＋Ｋ４(x,y)・Ｗ４＝Ｑを求める。総和Ｑを求める計算は加算部３２で行なう。

【００６６】次に加算部３２での加算結果Ｑにαを乗算
し、閾値修正値ΔＴｈ(x,y)を求める。この乗算は乗算
部３３で行なう。ΔＴｈ(x,y)と固定閾値Ｔｈを加算し
て、修正閾値Ｔｈ'(x,y)を求める。この加算は加算部３
４が行う。

【００６７】step 7 ：Ｔｈ'(x,y)とI'(x,y)の大小比較
を行う。大小比較は比較部１２が行う。Ｔｈ'(x,y)＜I'
(x,y)ならば、step 8に分岐する。そうでなければstep
10に分岐する。

【００６８】step 8 ：座標(x,y)の２値化出力シンボル
を１とする。この処理は、０／１処理部１３が行う。

【００６９】step 9 ：座標(x,y)での２値化誤差をＥ
(x,y)=Ｉ'(x,y)-Ｉmaxで求め、誤差メモリ１６の所定の
位置に格納する。Ｅ(x,y)の計算（Ｉ'(x,y)-Ｉmaxの減
算）は減算部１５が行う。その後、step 12 （図１１）
に進む。

【００７０】step 10 ：座標(x,y)の２値化出力シンボ
ルを０とする。この処理は０／１処理部１３が行う。

【００７１】step 11：座標(x,y)での２値化誤差をＥ
(x,y)=Ｉ'(x,y)-Ｉminで求め、誤差メモリ１６の所定の
位置に格納する。Ｅ(x,y)の減算は減算部１５が行う。
その後step 12に進む。

【００７２】step 12：xがx_maxに達しているかどうか判
定する。ここでx_maxは画面上の右端のx座標値である。
達していれば、step 14に分岐する。達していなければs
tep 13に分岐する。

【００７３】step 13：xを１だけ大きくする。その後st
ep 2に戻る。

【００７４】step 14：yがy_maxに達しているかどうか判
定する。ここでy_maxは画面上の下端のx座標値である。
達していれば、終了する。達していなければstep 15に
分岐する。

【００７５】step 15：xを0に戻し、yを１だけ大きくす
る。その後step 2に戻る。

【００７６】以上にように、step 12〜step 15におい
て、上記step 1 〜step 11の処理が画面上のすべの画素
について行なわれたかどうかを判断し、終っていなけれ
ば、座標値制御を行ない、step 3に分岐し、次の画素に
ついて処理を行なうこととする。すべての画素について
処理が終った場合、一旦処理を終了し、次の画面につい
ての処理の開始命令を待つ。

【００７７】以上説明した実施の形態によれば、階調特
性に依存した二次元空間周波数の画素配置構造の枠組み
に拘束され、拘束位置を中心にランダム性を持って変動
する画素配置構造を持つ（言換えれば、大きな枠の規則
性を有する）画像を生成することができる。このことに
より、階調値に依存したプリンタドットのノイズ特性の
キャンセルや、階調値に依存したカラードットの重なり
合わせ特性の制御を行うことができ、プリンタデバイス
の特性に合わせた、高品質の擬似階調画像を生成でき
る。

【００７８】また、上記のようなテーブルＴ１〜Ｔ４を
用いることにより、各テーブルの第２象限、第４象限の
位置の画素は信号値に応じて白く変換され易くなり、第
１象限、第３象限の画素は信号値に応じて黒く変換され
易くなる。この様な枠組みがテーブルＴ４〜テーブルＴ
１に至る階層構造で得られこれら各象限の枠組みと言う
規則の中でランダムな画素変換がなされるようになる。

【００７９】尚、擬似階調画像処理部として、誤差拡散
法による処理を行なうものを用い、誤差拡散法の代りに
ディザ法による処理を行なうものを用いてもよい。この
場合は、周期信号の信号値がディザマトリクスの数だけ
存在するが、擬似階調処理を行う対象のディザ閾値にフ
ーリエ展開された信号値を加算した結果の閾値を用いて
擬似階調処理を行う。

【００８０】第２の実施の形態図１２は、本発明の第２の実施の形態を示すブロック図
である。図７と同一の符号は同一又は対応する部材を示
す。

【００８１】第１２の画像処理装置は図７のものと類似
である。異なるのは、誤差拡散部１０内の、減算部１５
と誤差メモリ１６の間に、累加算部１７を挿入し、フィ
ルタ１７と加算部１１の間に割算部２１を挿入し、制御
部１０とアドレス生成部４の間に座標値変換部６を挿入
したことである。

【００８２】この実施の形態は、誤差拡散が画素毎にで
はなく、ブロック毎に行なわれる。ここでブロックと
は、例えば図１３に示す縦ｍ個、横ｎ個（ｍ，ｎ≧１）
の画素マトリクスからなる矩形の領域である。画面全体
には、このような複数の矩形領域は縦横に配列されてい
る。

【００８３】誤差拡散部１０における画素の処理順は、
第１の実施の形態のように水平走査線毎にではなく、ブ
ロック毎に行なう。即ち、最初に最も上の行の、最も左
のブロックが選択され、その中の画素を順次選択する。
ブロック内における画素の選択順は、水平走査線を上か
ら下へ、各水平走査線内では左から右へと言う順であ
る。

【００８４】このような順で処理を行なうため、制御部
１０は、処理すべきブロックのアドレスxb, yb及びブロ
ック内の画素の順位i (iは0からm×nの範囲で変化す
る）を順に出力する。座標変換部６は、xb,yb,iに基づ
きx,yを出力する。

【００８５】累加算部１７は、減算部１５から供給され
る誤差を累積加算し、ブロックごとの総和を求める。こ
のため、各ブロックの始めに制御部１０によりリセット
され、そのブロックのすべての画素についての誤差の累
積加算が終ったら、これが誤差メモリ１６に記憶され
る。

【００８６】誤差メモリ１６は、各画素についての誤差
ではなく、各ブロックについての誤差を記憶する。着目
画素を含むブロックと同じ行で左側の全てのブロック
と、一つ上の行の全てのブロック全部とについて、誤差
を記憶している。（図１４参照。）割算部２１は、処理ブロックの誤差を画素数Ｌで割り、
１画素あたりの誤差を求める。

【００８７】上記のように、誤差拡散部１０においてブ
ロック毎に処理を行なうのに伴い、閾値設定部３０にお
ける画素の処理順も第１の実施の形態とは異なる。しか
し、その他の点では、第２の実施の形態は第１の実施の
形態と同じである。

【００８８】以下本実施の形態の動作を、図１５及び図
１６のフローチャートを参照して説明する。

【００８９】step 101 ：入力原画像の処理ブロックの
座標値を初期化する。初期化処理は制御部１０が行う。

【００９０】step 102 ：誤差メモリ１６の内容を全て
０で初期化する。メモリの初期化処理は制御部１３が行
う。

【００９１】step 103 ：誤差メモリ１６の内容のう
ち、着目画素を含むブロックの近傍の４つのブロック
（図１４参照）の誤差値を読み出し、重み付き和Ｅb(xb,yb)＝Ａb(xb-1,yb-1)・Ｅb(xb-1,yb-1)＋Ａb(x
b,yb-1)・Ｅb(x,yb-1)＋Ａb(xb+1,yb-1)・Ｅb(xb+1,yb-
1)＋Ａb(xb-1,yb)・Ｅb(xb-1,yb) を求め、これをフィルタリング値Ｅb(x,y)とする。ここ
でＥb(xb-1,yb-1)，Ｅ(xb,yb-1)，Ｅ(xb+1,yb-1)，Ｅ(x
b-1,yb)はそれぞれ着目画素を含むブロックの近傍のブ
ロックの誤差であり、Ａb(xb-1,yb-1)，Ａb(xb,yb-1)，
Ａb(xb+1,yb-1)，Ａ(x-1,y)は近傍の画素に与えられる
重み付け計数である。このフィルタリングは、フィルタ
１７が誤差メモリ１６を使用して行う処理である。

【００９２】step 104 ：座標(xb,yb)のブロック内の画
素インデックスｉを初期化する。この画素インデクスｉ
は、ブロック内における画素の処理順を示す番号であ
り、左への画素が１、右下の最後の画素がｍ×ｎであ
る。

【００９３】step 105 ：インデックスｉの画素の階調
値Ｉ(xb,yb,i)をstep 103でもとめたＥb(xb,yb)で修正
し、Ｉ'(xb,yb,i)とする。この修正は加算部１１が行
う。ここでＬは１ブロック内に含まれる画素数であり、
Ｅb(xb,yb)の１／Ｌが各画素での修正値となる。ｍ×ｎ
が２のべきである場合には、Ｌで割る計算は、割算部８
をシフタで構成することにより容易に行なえる。

【００９４】step 106 ：xb,yb, i から４個のテーブル
を読み、テーブル値をＫ１〜Ｋ４を得る。

【００９５】step 107 ：４個のテーブル値Ｋ１〜Ｋ４
に、乗算部Ｍ１〜Ｍ４でそれぞれ重み付け計数を掛け
て、加算部３２で総和Ｑを求める。

【００９６】step 108 ：加算部３２の出力Ｑにベース
係数αをかけ、閾値の修正値ΔＴｈ(xb,yb,i)を求め
る。次に、加算部３４で閾値修正値ΔＴｈ（xb,yb,i)を
固定閾値Ｔｈに加算して、修正閾値Ｔｈ'(xb,yb,i)とす
る。

【００９７】step 109 ：比較部１２で、Ｔｈ(xb,yb,i)
とＩ'(xb,yb,i)の大小比較を行う。Ｔｈ'(xb,yb,i)＜
Ｉ'(xb,yb,i)ならば、step 110に分岐し、そうでなけれ
ばstep112に分岐する。

【００９８】step 110：座標(xb,yb,i)の２値化出力シ
ンボルを１とする。

【００９９】step 111：２値化誤差Ｅ(xb,yb,i)をＥ(x
b,yb,i)＝Ｉ'(xb,yb,i)−Ｉ_maxで求める。その後、step
114に進む。

【０１００】step 112：座標(xb,yb,i)の２値化出力シ
ンボルを０とする。

【０１０１】step 113：２値化誤差Ｅ(xb,yb,i)をＥ(x
b,yb,i)＝Ｉ'(xb,yb,i)−Ｉminで求める。

【０１０２】step 114：インデックスｉを１つインクリ
メントし、ｉ≦Ｌならばstep 115に戻る。

【０１０３】step 115：そうでなければ、step 116に進
む。

【０１０４】step 116：１ブロック内で発生した誤差の
総和Ｅ(xb,yb)を求め、誤差メモリの所定の位置に格納
する。

【０１０５】step 117：xbがxb_maxに達しているかどう
か判定する。ここでxb_maxは画面上の右端のブロックのx
b座標値である。達していれば、step 119に分岐する。
達していなければstep 118に分岐する。

【０１０６】step 118：xbを１だけ大きくする。その後
step 103に戻る。

【０１０７】step 119：ybがyb_maxに達しているかどう
か判定する。ここでyb_maxは画面上の下端のブロックのx
b座標値である。達していれば、終了する。達していな
ければstep 120に分岐する。

【０１０８】step 120：xbを0に戻し、ybを１だけ大き
くする。その後step 103に戻る。

【０１０９】このように、step 117〜step 120では、全
てのブロックについて、処理が終了したならば処理を終
了する。そうでなければ、ブロック座標位置制御を行
い、step 113に戻る。

【０１１０】第２の実施の形態によれば、階調特性に依
存した二次元空間周波数の構造の枠組みの中でランダム
性を持って変動する画素配置構造を持つ画像を生成する
ことができる。このランダム性は、処理ブロックが大き
くなる程減少し、ランダム製の制御も可能となる。また
ランダム性も制御可能である。このことにより、階調値
に依存したプリンタドットのノイズ特性のキャンセル
や、階調値に依存したカラードットの重なり合わせ特性
の制御を行うことができ、プリンタデバイスの特性に合
わせた、高品質の擬似階調画像を生成できる。

【０１１１】なお、上記の説明では、誤差を拡散された
ブロック内で各画素に対しブロック内の画素数で割っ
て、誤差を各画素に均等に割り当てているが、各画素に
対し重みを付けて分配しても良い。この重みは、例え
ば、ブロック内の画素の位置を考慮して定めることして
も良い。

【０１１２】第３の実施の形態第１の実施の形態、第２の実施の形態では、２値化を行
う実施の形態を示したが、多値（＞２）化を行うものに
拡張することもできる。

【０１１３】図１７に多値の一次元誤差拡散の説明図を
示す。この例は４値（g0，g1，g2，g3）化を行うもので
あり、g0とg1の中間、g1とg2の中間、g2とg3の中間にそ
れぞれ閾値Ｔｈ1，Ｔｈ2，Ｔｈ3が設定されている。４
値のいずれかに量子化されたときの量子化誤差が次に処
理する画素の階調値を修正し、以下同様の処理を繰り返
す。図１７の原画像階調値に対しては、誤差の拡散によ
り、階調値が図の破線のように逐次修正され、量子化シ
ンボルは、g3，g2，g3…となっている。

【０１１４】この多値誤差拡散法についても、発明が解
決しようとする課題で述べた問題があり、第１の実施の
形態、第２の実施の形態と同様の構成で解決することが
できる。

【０１１５】第１の実施の形態、第２の実施の形態にお
いて、比較部１２をＴｈ1，Ｔｈ2，Ｔｈ3の閾値につい
ての４値化の比較機能を持つようにし、０／１処理部１
３の代りに0/1/2/3処理部を用い、４値のシンボルを生
成するようにする。他は第１の実施の形態の図７のブロ
ック図と同様の機能で、第１の実施の形態、第２の実施
の形態の多値化が行える。

【０１１６】第４の実施の形態入力画像が複数色プレーンによって構成されるカラー画
像になっている場合を考える。図１８は第１の実施の形
態において、カラー画像が入力される場合の実施の形態
のブロック図を示すものである。この場合、色プレーン
は３であり、各色プレーンに対応して別個の入力画像メ
モリ２Ｃ（シアン用）、２Ｍ（マゼンタ用）、２Ｙ（イ
エロー用）を有する。

【０１１７】３つの色の信号は順に処理される。即ち、
ある色（例えばＣ）について１画面分の処理が行なわ
れ、同じ画面について別の色（Ｍ）が処理され、最後に
第３の色（Ｙ）が処理され、その後次の画面の処理に移
る。

【０１１８】図示の装置は図７に示すものと類似である
が、入力画像メモリが３つのプレーン２Ｃ、２Ｍ、２Ｙ
を有し、その出力側に選択器２２が設けられ、制御部８
は、座標値ｘ，ｙのほか色を指定する信号ｊを出力し、
閾値設定部３０は割算器３５を有する。

【０１１９】制御部８からの座標値ｘ，ｙに基づき、ア
ドレス発生器４からアドレスｘ’，ｙ’が発生され、こ
れが３つの色プレーン２Ｃ、２Ｍ、２Ｙに同時に入力さ
れる。制御部８からの信号ｊにより選択器２２が制御さ
れる。即ち処理すべき色に応じて色プレーン２Ｃ、２
Ｍ、２Ｙの出力の何れかが選択される。選択された色信
号が符号Ｚで表わされている。

【０１２０】Ｚ＝Ｃ、Ｍ又はＹである。

【０１２１】信号Ｚに対する誤差拡散部１０内の処理は
第１の実施の形態で述べたのと同様である。

【０１２２】色プレーン２Ｃ、２Ｍ、２Ｙの出力信号及
び選択器２２の信号Ｚは割算器３５に供給される。割算
器では、Ｚ／（Ｃ＋Ｍ＋Ｙ）×Ｉｍａｘの計算を行ない、その出力を重みテーブル３１に供給す
る。重みテーブル３１では、入力されるＺ／（Ｃ＋Ｍ＋
Ｙ）×Ｉｍａｘに基づいて重み値Ｗ１〜Ｗ４を出力す
る。これ以外の点では、閾値設定部３０は第１の実施の
形態と同様に動作する。

【０１２３】Ｚ／（Ｃ＋Ｍ＋Ｙ）の計算を行なうことにより、インクの総量に示す各色の
比率で重み値を制御することができ、他の色のインクに
よって占有される面積も考慮して画素配置を制御するこ
とができる。即ち、インクの総量に示す比率が低い場
合、比較的分散して配置し、インクの総量に示し比率が
高い場合、比較的集中して配置することとなる。

【０１２４】上記の実施の形態によれば、Ｃ、Ｍ、Ｙの
各色の値が小さい場合、各色のドットをできるだけ重ね
ないように制御する（色毎に２次元空間配置の位相がず
れるように画素配置構造の枠組を設定）ことで、ドット
を分散配置し、粒状感の少ないカラー階調を再現するこ
とができる。

【０１２５】また、Ｃ、Ｍ、Ｙの各色の値が大きい場
合、色ドットをできるだけ重ねるように制御すること
で、ドットを集中配置（色毎に２次元空間の位相が合う
ように画素配置構造の枠組みを設定）し、階調特性に優
れたカラー階調を再現することができる。

【０１２６】なお、第２の実施の形態についても、テー
ブル部Ｔ１〜Ｔ４の構成を同様にすることで、色プレー
ン間の色重ね合わせの制御を行うことができる。

【０１２７】

【発明の効果】以上のように本発明によれば、周期が互
いに異なる周期信号を、階調値に依存した重み値で重畳
し、２値化固定閾値に加え固定閾値の値を修正すること
で、画素値配置構造の拘束条件の中で、誤差拡散法によ
るランダムに変動する画素配置を実現することができ、
これにより、階調値に依存したプリンタドットのノイズ
特性のキャンセルや、階調値に依存したカラードットの
重なり合わせ特性の制御を行うことができ、プリンタデ
バイスの特性に合わせた、高品質の擬似階調画像を生成
できる。

【図面の簡単な説明】

【図１】本発明における一次元の誤差拡散法で得られ
る出力画素の一例、及び本発明で用いられる周波数の異
なる２つの周期信号及びそれによる閾値の修正を示す図
である。

【図２】本発明による二次元誤差拡散法で用いられる
異なる大きさのパターンマトリクスを示す図である。

【図３】２×２パターンマトリクスの画面上における
配列を示す図である。

【図４】４×４パターンマトリクスの画面上における
配列を示す図である。

【図５】パターンの各象限の、閾値修正値の元となる
信号値の設定の仕方を示す図である。

【図６】入力画素の階調値に対する各パターンマトリ
クスに対するの重み値の変化を示す図である。

【図７】第１の実施の形態の画像処理装置を示すブロ
ック図である。

【図８】一つの構成例における、パターンマトリクス
テーブルと、該テーブルに与えるアドレス信号を示す図
である。

【図９】他の構成例における、パターンマトリクステ
ーブルと、該テーブルに与えるアドレス信号を示す図で
ある。

【図１０】図７の装置の動作を示すフローチャートで
ある。

【図１１】図７の装置の動作を示すフローチャートで
ある。

【図１２】第２の実施の形態の画像処理装置を示すブ
ロック図である。

【図１３】画素のブロックの構成及びその配列を示す
図である。

【図１４】着目画素を含むブロックと、その処理の際
誤差を参照する近傍のブロックとの位置関係を示す図で
ある。

【図１５】図１２の装置の動作を示すフローチャート
である。

【図１６】図１２の装置の動作を示すフローチャート
である。

【図１７】多値の一次元誤差拡散の説明図である。

【図１８】カラー画像を処理する装置を示すブロック
図である。

【図１９】一次元誤差拡散法を説明するための図であ
る。

【図２０】ある画素と、その画素における誤差を拡散
する相手の画素との位置関係を示す図である。

【図２１】着目画素とその処理の際誤差を参照する近
傍の画素との位置関係を示す図である。

【符号の説明】

２：入力画像メモリ；８：制御部；１０：誤差
拡散部；３０：閾値設定部；Ｔ１〜Ｔ４：テー
ブル；３１：重み値テーブル；Ｍ１〜Ｍ４，３
３：乗算部；３２、３４：加算部。

Claims

【特許請求の範囲】

【請求項１】周期の長さが互いに異なる複数の周期信
号に基づいて、閾値の値を変動させる閾値設定部（３
０）と、該変動させた閾値を用いて擬似階調画像処理を行う擬似
階調画像処理部（１０）とを有する擬似階調画像処理装
置。
【請求項２】上記周期信号が水平走査線に沿う、一次
元的なものであることを特徴とする請求項１にの擬似階
調画像処理装置。
【請求項３】上記周期信号が、二次元的なものである
ことを特徴とする請求項１に記載の擬似階調画像処理装
置。
【請求項４】上記閾値設定部は、原画像信号の座標値
に対応する上記複数の周期信号のそれぞれの信号値に、
原画像階調値に依存する重み値を掛けて加算し、該加算
結果によりｎ値化（ｎは２以上の自然数）のための閾値
を修正することを特徴とする請求項１に記載の擬似階調
画像処理装置。
【請求項５】上記閾値設定部は、上記加算結果に、調
整係数（α）を乗算し、この乗算結果（ΔＴｈ）を固定
の閾値（Ｔｈ）に加算することにより、上記ｎ値化のた
めの閾値を得ることを特徴とする請求項４に記載の擬似
階調画像処理装置。
【請求項６】上記複数の周期信号のそれぞれの信号値
に掛けられる重み値が、所定の原画像階調値の所定の範
囲においてゼロ以外の値であり、上記所定の範囲以外に
おいてゼロであることを特徴とする請求項４に記載の擬
似階調画像処理装置。
【請求項７】上記所定の範囲において、上記重み値が
正弦波状に変化することを特徴とする請求項６に記載の
擬似階調画像処理装置。
【請求項８】上記複数の周期信号のうち、周期のより
短いもののための重み値がゼロ以外の値となる上記所定
の範囲が、より高い原画像階調値を含むものであること
を特徴とする請求項７に記載の擬似階調画像処理装置。
【請求項９】上記複数の周期信号のための重み値がゼ
ロ以外の値となる上記所定の範囲が、一部互いに重なり
あっていることを特徴とする請求項８に記載の擬似階調
画像処理装置。
【請求項１０】上記擬似階調画像処理部（１０）が、
誤差拡散法による画像処理を行なうものであることを特
徴とする請求項１に記載の擬似階調画像処理装置。
【請求項１１】上記ｎが２であることを特徴とする請
求項１に記載の擬似階調画像処理装置。
【請求項１２】上記ｎが３以上の値であることを特徴
とする請求項１に記載の擬似階調画像処理装置。
【請求項１３】上記擬似階調画像処理部（１０）が、
処理画像をｍ画素×ｎライン（1≦m,n）のブロックの各
々の内部で発生したｎ値量子化誤差の総和を近傍のブロ
ックに重み値を掛けて拡散し、拡散されたブロック内で
拡散誤差を所定の比率でブロック内の各画素に分配し、
各画素では分配された誤差を原画像階調値に加算するこ
とにより、原画像の画素の階調値を修正する誤差拡散を
行うことを特徴とする請求項３に記載の擬似階調画像処
理装置。