JP2002247391A

JP2002247391A - 画像処理方法及びその装置

Info

Publication number: JP2002247391A
Application number: JP2001040820A
Authority: JP
Inventors: Suzuko Fukao; 珠州子深尾
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2001-02-16
Filing date: 2001-02-16
Publication date: 2002-08-30

Abstract

(57)【要約】【課題】多値画像の擬似中間調表現処理に際し、全階
調において良好な処理結果を得られる画像処理方法及び
その装置を提供する。【解決手段】入力されたカラー多階調画像を、入力画
像の階調数より少ない階調数を有する画像に変換する誤
差拡散法による２値化処理１２において、２値化処理を
行った画像データ１５の評価値を求め１６，１７、評価
値に基づいて前記２値化処理の拡散係数を選択し１３，
１４，１８、入力画素の各色成分の階調に応じて、選択
された拡散係数を使って２値化処理を行う。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は画像処理方法及びそ
の装置、特に、入力されたカラー多階調画像を入力画像
の階調数より少ない階調数を有する画像に変換する画像
処理方法及びその装置に関するものである。

【０００２】

【従来の技術】従来、多値画像データを２値画像（又は
２値以上で入力階調数より少ない階調数を有する画像）
に変換する手段として、R. Floydらによる誤差拡散
法（"An adaptive algorithm for spatila gray scal
e"，SID International Symposium Digest of Technica
l Papers, vol4.3, 1975, pp.36-37）がある。この誤差
拡散法は、ある画素で生じた２値化誤差を以降の複数画
素へ拡散することにより、擬似的に階調表現を行うもの
である。一般的な誤差拡散法では、２値化誤差を拡散す
る際の重み付け（拡散係数）を入力データの値に関わら
ず固定としていた。

【０００３】しかし、拡散係数を固定とする誤差拡散法
では、ハイライト部及びシャドー部においてドットが均
一に分散されず、鎖状に連なって発生するなどの問題点
がある。この問題点を解決するため、特開平１０−３１
２４５８号には、複数の入力階調で異なる拡散係数を用
いる誤差拡散法が開示されている。この方法では、ハイ
ライト部及びシャドー部における第１の拡散係数と、中
間階調部における第２の拡散係数を定義し、ハイライト
部から中間階調部、及び中間階調部からシャドー部への
境界領域では、この２つの拡散係数を線形補間して用い
ている。

【０００４】

【発明が解決しようとする課題】しかしながら、上記の
方法では、ハイライト部及びシャドー部、又中間階調領
域に用いる拡散係数は定義されているが、これ以外の入
力階調領域においては、２つの拡散係数を線形補間した
拡散係数が用いられており、全ての階調において最適な
拡散係数が定義されているわけではない。又、既に拡散
係数が定義されている階調においても、定義された拡散
係数が最適であるという根拠はない。更に、画像の解像
度などの諸条件によっても最適な拡散係数が異なる、と
いう問題点もある。

【０００５】本発明は、上記課題を解決するためになさ
れたものであり、多値画像の擬似中間調表現処理に際
し、全階調において良好な処理結果を得られる画像処理
方法及びその装置を提供することを目的とする。

【０００６】

【課題を解決するための手段】上記目的を達成するため
に、本発明の画像処理装置は、入力されたカラー多階調
画像を、入力画像の階調教より少ない階調数を有する画
像に変換する擬似中間調表現処理を行う画像処理装置で
あって、入力画像データの特徴に応じて、擬似中間調表
現処理のパラメータを変更することが可能な画像処理手
段と、該画像処理手段で擬似中間調表現処理を行った出
力画像の評価値を求める評価値演算手段と、該評価値に
基づいて、前記画像処理手段に設定される前記擬似中間
調表現処理のパラメータを選択するパラメータ選択手段
とを有することを特徴とする。

【０００７】ここで、前記擬似中間調表現処理は色成分
毎の誤差拡散法であって、前記パラメータは、色成分毎
の誤差拡散係数である。また、前記入力画像データの特
徴は、色成分毎の入力画素値である。また、前記画像処
理手段は、色成分毎の入力画素値に対応して誤差拡散係
数を記憶するルックアップテーブルを有し、前記パラメ
ータ選択手段は、選択されたパラメータを、色成分毎の
入力画素値に対応して該ルックアップテーブルに記憶す
る。また、前記誤差拡散法は、入力画素値と出力画素値
の誤差の初期値を０とする。また、前記誤差拡散法は、
入力画素値と出力画素値の誤差の初期値を乱数とする。
また、前記パラメータ選択手段は、色成分毎の各入力画
素値に対応する前記評価値に基づいて、誤差拡散係数を
選択する。また、前記パラメータ選択手段は、色成分毎
の複数の入力画素値に対応する前記評価値に基づいて誤
差拡散係数を選択する。また、前記パラメータ選択手段
は、異なる色成分間の全ての画素値の組み合わせにおけ
る評価値に基づいて各色成分の誤差拡散係数を選択す
る。また、前記パラメータ選択手段は、異なる色成分間
の同一の画素値における評価値に基づいて各色成分の誤
差拡散係数を選択する。また、前記評価値演算手段は、
前面が所定の画素値の入力画像に擬似中間調表現処理を
行った出力画像を周波数領域へ変換する変換手段と、周
波数領域において評価値を演算する演算手段とを含む。
また、前記変換手段は、前記出力画像の一部からパワー
スペクトルを求め、前記演算手段は、前記パワースペク
トルをフィルタ処理しながら、所定の周波数範囲にわた
って加算する。また、前記変換手段は、前記出力画像の
一部から２次元フーリエ変換により２次元パワースペク
トルを求め、更に、該２次元パワースペクトルから周波
数帯域毎に平均値を取った１次元パワースペクトルを求
め、前記演算手段は、前記１次元パワースペクトルを視
覚特性フィルタ処理しながら、所定の周波数範囲にわた
って加算する。また、前記パラメータ選択手段は、前記
演算手段の加算値が最小の誤差拡散係数を選択する。

【０００８】又、本発明の画像処理方法は、入力された
カラー多階調画像を、入力画像の階調数より少ない階調
数を有する画像に変換する擬似中間調表現処理を行う画
像処理方法であって、擬似中間調表現処理を行った画像
データの評価値を求める評価値演算工程と、該評価値に
基づいて前記擬似中間調表現処理のパラメータを選択す
るパラメータ選択工程と、入力画像データの特徴に応じ
て、選択された擬似中間調表現処理のパラメータを使っ
て擬似中間調表現処理を行う画像処理工程とを有するこ
とを特徴とする。

【０００９】ここで、前記擬似中間調表現処理は色成分
毎の誤差拡散法であって、前記パラメータは、色成分毎
の誤差拡散係数である。また、前記入力画像データの特
徴は、色成分毎の入力画素値である。また、前記誤差拡
散法は、入力画素値と出力画素値の誤差の初期値を０と
する。また、前記誤差拡散法は、入力画素値と出力画素
値の誤差の初期値を乱数とする。また、前記パラメータ
選択工程では、色成分毎の各入力画素値に対応する前記
評価値に基づいて、誤差拡散係数を選択する。また、前
記パラメータ選択工程では、色成分毎の複数の入力画素
値に対応する前記評価値に基づいて誤差拡散係数を選択
する。また、前記パラメータ選択工程では、異なる色成
分間の全ての画素値の組み合わせにおける評価値に基づ
いて各色成分の誤差拡散係数を選択する。また、前記パ
ラメータ選択工程では、異なる色成分間の同一の画素値
における評価値に基づいて各色成分の誤差拡散係数を選
択する。また、前記評価値演算工程は、前面が所定の画
素値の入力画像に擬似中間調表現処理を行った出力画像
を周波数領域へ変換する変換工程と、周波数領域におい
て評価値を演算する演算工程とを含む。また、前記変換
工程では、前記出力画像の一部からパワースペクトルを
求め、前記演算工程では、前記パワースペクトルをフィ
ルタ処理しながら、所定の周波数範囲にわたって加算す
る。また、前記変換工程では、前記出力画像の一部から
２次元フーリエ変換により２次元パワースペクトルを求
め、更に、該２次元パワースペクトルから周波数帯域毎
に平均値を取った１次元パワースペクトルを求め、前記
演算工程では、前記１次元パワースペクトルを視覚特性
フィルタ処理しながら、所定の周波数範囲にわたって加
算する。また、記パラメータ選択工程は、前記演算手段
の加算値が最小の誤差拡散係数を選択する。

【００１０】又、本発明の記憶媒体は、入力されたカラ
ー多階調画像を、入力画像の階調数より少ない階調数を
有する画像に変換する擬似中間調表現処理を行うための
プログラムを記憶するコンピュータ読取り可能な記憶媒
体であって、前記プログラムは、少なくとも、擬似中間
調表現処理を行った画像データの評価値を求める評価値
演算ステップと、該評価値に基づいて前記擬似中間調表
現処理のパラメータを選択するパラメータ選択ステップ
とを有することを特徴とする。

【００１１】ここで、前記擬似中間調表現処理は色成分
毎の誤差拡散法であって、前記パラメータは、色成分毎
の誤差拡散係数である。また、前記評価値演算ステップ
は、前面が所定の画素値の入力画像に擬似中間調表現処
理を行った出力画像を周波数領域へ変換する変換ステッ
プと、周波数領域において評価値を演算する演算ステッ
プとを含む。また、前記変換ステップでは、前記出力画
像の一部から２次元フーリエ変換により２次元パワース
ペクトルを求め、更に、該２次元パワースペクトルから
周波数帯域毎に平均値を取った１次元パワースペクトル
を求め、前記演算ステップでは、前記１次元パワースペ
クトルを視覚特性フィルタ処理しながら、所定の周波数
範囲にわたって加算する。また、前記パラメータ選択ス
テップでは、前記演算手段の加算値が最小の誤差拡散係
数を選択する。また、前記パラメータ選択ステップで
は、色成分毎の各入力画素値に対応する前記評価値に基
づいて、誤差拡散係数を選択する。また、前記パラメー
タ選択ステップでは、色成分毎の複数の入力画素値に対
応する前記評価値に基づいて誤差拡散係数を選択する。
また、前記パラメータ選択ステップでは、異なる色成分
間の全ての画素値の組み合わせにおける評価値に基づい
て各色成分の誤差拡散係数を選択する。また、前記パラ
メータ選択ステップでは、異なる色成分間の同一の画素
値における評価値に基づいて各色成分の誤差拡散係数を
選択する。また、入力画像データの特徴に応じて、選択
された擬似中間調表現処理のパラメータを使って擬似中
間調表現処理を行う画像処理ステップを更に有する。

【００１２】又、入力されたカラー多階調画像を、入力
画像の階調数より少ない階調数を有する画像に変換する
擬似中間調表現処理を行うためのデータを記憶するコン
ピュータ読取り可能な記憶媒体であって、前記データ
は、擬似中間調表現処理を行った画像データの評価値を
求め、該評価値に基づいて選択された前記擬似中間調表
現処理のパラメータを含むことを特徴とする。

【００１３】

【発明の実施の形態】以下、図面を参照しながら本発明
に係る実施の形態を詳細に説明する。

【００１４】［実施の形態１］＜本実施の形態の画像処理装置の構成例＞図１は、本発
明の実施の形態の画像処理装置のブロック図である。

【００１５】図１において、１０は画像生成ユニットで
あり、画像生成ユニット１０において生成された入力画
像データは画像メモリ１１に格納される。１２は画像２
値化ユニットであり、後述する方法により画像メモリ１
１に格納されているカラー多値画像データのシアン
（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック
（Ｋ）の各色成分毎の２値化を行う。１３は拡散係数テ
ーブルであり、画像２値化ユニット１２において画素の
２値化誤差配分処理を行う際の、色成分毎の各入力画素
値に対する重み付け（拡散係数）が格納されたテーブル
である。１４は候補拡散係数生成ユニットであり、拡散
係数テーブル１３に格納する係数の候補として検討する
拡散係数を順次生成する。

【００１６】１５は２値画像メモリであり、画像２値化
ユニット１２において色毎に２値化された画像データが
格納され、該２値化された画像データが出力画像データ
として以降の処理、例えばプリントや表示などのために
使用される。１６は画像評価値演算ユニットであり、画
像の評価値を求める。画像評価値演算ユニット１６にお
いて求めた画像評価値データは、画像評価値メモリ１７
に格納される。１８は画像評価値メモリ１７に格納され
た画像評価値データをもとに、候補拡散係数の中から各
色成分の画素値ごとに１組の拡散係数（本例では、Ｋ₁
からＫ₄)を選択する拡散係数選択ユニットであり、拡散
係数選択ユニット１８において選ばれた拡散係数が最終
的に拡散係数テーブル１３に格納される。

【００１７】（誤差拡散係数テーブル）図２は、本実施
の形態による誤差拡散係数テーブル１３の一例を示す図
である。

【００１８】図２に示すように、横に隣接する４つのま
す目の値が注目画素の色成分毎の２値化誤差を周辺画素
に分配する際の拡散係数であり、０から１の範囲内の実
数が格納される（これを計算する式については、追って
図１０で示される）。欄外左の数字は、色毎の入力画素
値（０〜２５５：８ビット）を表す。又、欄上の文字
は、Ｋが拡散係数を表し、下付きのＫ１〜Ｋ４がブラッ
クの２値化誤差を分配する画素の注目画素に対する位置
を表し、下付きのＣ１〜Ｃ４がシアンの２値化誤差を分
配する画素の注目画素に対する位置を表し、下付きのＭ
１〜Ｍ４がマゼンタの２値化誤差を分配する画素の注目
画素に対する位置を表し、下付きのＹ１〜Ｙ４がイエロ
ーの２値化誤差を分配する画素の注目画素に対する位置
を表す。１から４の位置ついては、図８及び図９に従い
後述する。

【００１９】（画像２値化ユニット）図３は、本実施の
形態による画像２値化ユニット１２の構成を示したブロ
ック図である。

【００２０】図３において、２１はデータ入力端子であ
り、０から２５５の範囲内の８ビットの入力値を読み込
む。２２は入力バッファであり、入力カラー画像の１色
の１ライン分の入力データを記憶する。２３は入力補正
ユニットであり、１画素分の入力データに処理済みの画
素からの累積誤差を加算する。２４は２値化ユニットで
あり、入力補正ユニット２３において補正された入力デ
ータの２値化を行う。２５は出力バッファであり、２値
化ユニット２４において２値化されたデータを１ライン
分記憶する。２６は差分演算ユニットであり、注目画素
に対する２値化誤差を求める。２７は誤差配分ユニット
であり、差分演算ユニット２６で求めた２値化誤差を周
辺画素に拡散する処理を行う。２８は誤差バッファであ
り、注目画素の周辺画素に拡散した誤差を記憶するＲＡ
Ｍである。２９は出力レベル０又は２５５の２値信号の
出力端子である。

【００２１】尚、ここで入力カラー画像とは、各画素が
Ｃ，Ｍ，Ｙ，Ｋ各色につき０から２５５の範囲内の値
（８ビット）をもつ多値画像データを指す。又、図示さ
れない方法により入力カラー画像の横サイズ（画素数）
Ｗ、縦サイズＨが指定されているものとする。

【００２２】図４は、本実施の形態による入力バッファ
２２を示す図である。各ます目内の値がＣＭＹＫ何れか
１色の入力データであり、０から２５５のいずれかの整
数値が記憶される。従って、入力バッファ２２はＩ
（１）からＩ（Ｗ）までのＷ画素分ある。

【００２３】図５は、本実施の形態による出力バッファ
２５を示す図である。各ます目内の値が出力されるデー
タであり、０か２５５かのどちらかの値が記憶される。
従って、出力バッファ２５はＰ（１）からＰ（Ｗ）まで
のＷ画素分ある。

【００２４】図６は、本実施の形態による誤差バッファ
２８を示す図である。各ます目内の値が２値化処理済み
の画素からの累積誤差であり、−２５５から２５５の範
囲の実数が記憶される。又、本バッファには、ライン端
の処理に対応するために入力画像の横サイズＷよりも２
画素分余裕をもたせている。従って、誤差バッファ２８
はＥ（０）からＥ（Ｗ＋１）までのＷ＋２画素分ある。
又、注目画素の次の主走査方向の画素への誤差拡散（図
８及び図９のＫ₁とＫ₄）は、現在のラインと次のライン
への誤差バッファ２８への記憶が重複するので、次のラ
インへの誤差拡散はＥ_tempに一時貯えられて、次の画素
の処理時に次のラインの同じ位置の画素への誤差拡散
（Ｋ₃）に加えられることになる。

【００２５】（画像２値化ユニットでの２値化例）図７
は、図１の画像処理装置により入力カラー画像の色成分
毎の２値化手順を示したフローチャートである。このフ
ローチャートでは、副走査方向のライン番号をｙ、主走
査方向のドット番号をｘとし、２値化の処理は画像左上
端から右下へ順に行うものとする。

【００２６】図７に従って本画像処理装置の２値化手順
を説明する。

【００２７】まず、ライン番号ｙを１に初期化し、誤差
バッファ２８を０に初期化する（Ｓ１００）。つまり、
図６の誤差バッファのＥ（０）からＥ（Ｗ＋１）とＥ
_tempの値を０に設定する。次に、ｙライン目の画素デー
タを入力バッファ２２に読み込む（Ｓ１０１）。即ち、
ｙライン目のｘドット目の入力データをＩ（ｘ）に代入
する作業をｘ＝１からｘ＝Ｗについて行う。

【００２８】こうして１ライン分の入力データを読み込
んだ後、画素位置ｘを初期値１に設定し（Ｓ１０２）、
ｘ番目の入力データＩ（ｘ）に後述する処理済み画素か
らの拡散誤差Ｅ（ｘ）を加え、補正データＩ'（ｘ）と
する処理を行う（Ｓ１０３）。この補正データＩ'
（ｘ）を閾値"１２７"と比較し（Ｓ１０４）、Ｉ'
（ｘ）が"１２７"よりも大きければ２５５を結果値Ｐ
（ｘ）に代入し（Ｓ１０５）、小さければＰ（ｘ）を０
とする（Ｓ１０６）。こうして結果値を求めた後、注目
画素の補正データＩ'（ｘ）と結果値Ｐ（ｘ）との差
を、２値化誤差Ｅｒｒ（ｘ）として求める（Ｓ１０
７）。

【００２９】次に、この２値化誤差Ｅｒｒ（ｘ）を周辺
画素に分配する処理を行う。このとき、配分の重み付け
（拡散係数）は入力データＩ（ｘ）によって異なる。ま
ず入力データの値Ｉ（ｘ）の値に対応する拡散係数を拡
散係数テーブル１３より読み込み（Ｓ１０８）、その重
み付けに従い２値化誤差Ｅｒｒ（ｘ）を周辺画素に分配
する処理を行う。本実施の形態では、拡散係数テーブル
１３より読み込んだ拡散係数Ｋ₁，Ｋ₂，Ｋ₃，Ｋ₄（ここ
では、色成分の違いは省略している）を、図８に示すよ
うに注目画素周辺の未処理画素に２値化誤差Ｅｒｒ
（ｘ）を分配する。分配した誤差は、誤差バッファに次
のように格納される（Ｓ１０９）。

【００３０】Ｅ（ｘ＋１）←Ｅ（ｘ＋１）＋Ｋ₁×Ｅｒｒ（ｘ）Ｅ（ｘ−１）←Ｅ（ｘ−１）＋Ｋ₂×Ｅｒｒ（ｘ）Ｅ（ｘ＋１）←Ｅ_temp（＝Ｋ₄×Ｅｒｒ（ｘ−１））＋
Ｋ₃×Ｅｒｒ（ｘ）Ｅ_temp←Ｋ₄×Ｅｒｒ（ｘ）次に、注目画素位置ｘがラインの端まで到達しているか
を判断し（Ｓ１１０）、到達していなければ注目画素を
主走査方向に１つ移動し（Ｓ１１１）、ステップＳ１０
３からステップＳ１１０の処理を繰り返す。端まで処理
が終了した時点で、１ライン分の結果値Ｐ（１）からＰ
（Ｗ）を出力し、Ｅ_tempを０にリセットする（Ｓ１１
２）。その後、ライン番号ｙが最大値Ｈに到達したかを
判断する（Ｓ１１３）。ｙがＨより小さければ、ライン
番号ｙをインクリメントし（Ｓ１１４）、ステップＳ１
０１からステップＳ１１３の処理を繰り返す。ｙがＨに
到達していれば、本処理を終了する。

【００３１】尚、本実施の形態では、１ラインを処理す
る方向は固定としたが、ライン毎に左右交互に、あるい
は予め定められた順序で方向を切り替えてもよい。右端
から左方向に２値化処理を行う場合は、２値化の誤差を
図９のように、図８の拡散係数の位置を左右逆にした位
置に分配する。

【００３２】（拡散係数テーブルの作成例）次に、本実
施の形態による拡散係数テーブル１３の作成方法につい
て説明する。

【００３３】図１０は、候補拡散係数生成ユニット１４
において生成される全候補拡散係数を示したものであ
り、＊は注目画素を表す。ＣＭＹＫ各色の階調データ値
毎に、この全候補拡散係数の中から１組の係数が選ば
れ、拡散係数テーブル１３に格納される。

【００３４】図１１は、図１の実施の形態による拡散係
数設定の手順を示したフローチャートである。

【００３５】まずＫ成分の階調毎の拡散係数を後述の方
法により決定し、拡散係数テーブル１３に格納する（Ｓ
２００）。次にＣ成分の拡散係数を決定し（Ｓ２０
１）、以下、Ｍ（Ｓ２０２）、Ｙ（Ｓ２０３）成分の順
で色成分毎の拡散係数を決定し、拡散係数テーブル１３
に格納する。

【００３６】（Ｋ成分の拡散係数設定例）図１２は、図
１の実施の形態によるＫ成分の拡散係数設定の手順を示
したフローチャートである。

【００３７】まず、Ｋ成分の階調データ値ｇ_Kの初期値
を０に設定する（Ｓ３００）。その後、ｇ_Kをインクリ
メントする（Ｓ３０１）。

【００３８】次に、図１０に示した全候補拡散係数につ
いて順次以下の処理を行う。まず候補拡散係数を設定し
（Ｓ３０２）、拡散係数テーブルのＫ成分のｇ_Kに対応
する位置に格納する。次に、Ｋ成分は全画素がｇ_Kの値
をもち、Ｃ、Ｍ、Ｙ各成分は白である縦２０４８画素、
横５１２画素の入力カラー画像を作成する（Ｓ３０
３）。作成された画像データのＫ成分を前述の方法によ
り２値化し、この画像から図１４に示すように、下中央
部分の縦横２５６画素分を切り取る（Ｓ３０４）。ここ
で、２値化した画像の下中央部を切り取って評価するの
は、画像の上端及び左右の端では誤差拡散法特有のドッ
ト生成遅延が発生してしまうため、これを避けるのが目
的である。しかし、ドットの生成が安定した領域で評価
を行う方法はこれに限定されるものでなく、他の方法が
適用されてもよい。

【００３９】こうして作成された縦横２５６画素のカラ
ー画像の画像評価値を、後述する方法によって計算し
（Ｓ３０５）、画像評価値メモリ１７に格納する。その
後、全候補拡散係数分の評価値演算が終了したかを確認
し（Ｓ３０６）、終了していなければステップＳ３０１
からステップＳ３０５の処理を繰り返す。全候補拡散係
数についての評価値演算が終了していれば、ｇ_Kと２５
５を比較する（Ｓ３０７）。ｇ_Kが２５５より小さけれ
ば、ステップＳ３０１からステップＳ３０７の処理を繰
り返す。ｇ_Kの値が２５５に達していれば、画像評価値
メモリ１７に格納された評価値をもとに、階調データ値
ｇ_K毎に全候補拡散係数の中から、後述の方法により１
組の拡散係数を選択し、拡散係数テーブルに格納する
（Ｓ３０８）。

【００４０】尚、階調データ値ｇ_Kが０の場合の拡散係
数はｇ_K＝１の際と同じ拡散係数、又ｇ_Kの値が２５５の
場合はｇ_K＝２５４における拡散係数を拡散係数テーブ
ル１３に格納する。

【００４１】（Ｃ成分の拡散係数設定例）図１３は、図
１の実施の形態によるＣ成分の拡散係数設定の手順を示
したフローチャートである。

【００４２】まず、Ｃ成分の階調データ値ｇ_Cの初期値
を０に設定する（Ｓ４００）。その後、ｇ_Cをインクリ
メントする（Ｓ４０１）。

【００４３】次に、図１０に示した全候補拡散係数につ
いて順次以下の処理を行う。まず候補拡散係数を設定し
（Ｓ４０２）、拡散係数テーブルのＣ成分のｇ_Cに対応
する位置に格納する。次に、Ｋ成分の階調データ値ｇ_K
と画像評価値合計Ｒを０に初期化する（Ｓ４０３）。そ
の後、入力カラー画像を作成する（Ｓ４０４）。この入
力カラー画像のＫ成分は、全画素がｇ_Kの値をもつグレ
ー画像を、図１１のステップ２００において設定された
拡散係数を用いて２値化した画像とする。同Ｃ成分は全
画素がｇ_Cの値をもち、Ｍ、Ｙ成分は白とする。こうし
て作成されたカラー画像のＣ成分の２値化を行い、画像
の下中央部を切り取る（Ｓ４０５）。上記Ｋ成分と同様
に、縦２０４８画素、横５１２画素の入力カラー画像を
作成して２値化し、この画像から図１４に示すように、
下中央部分の縦横２５６画素分を切り取る。その後、上
記Ｋ成分の拡散係数設定の際と同様に、ステップＳ４０
５によってＣ成分が２値化されたカラー画像の画像評価
値を、後述する方法によって計算し、画像評価値合計Ｒ
に加算、画像評価値メモリ１７に格納する（Ｓ４０
６）。

【００４４】次に、ｇ_Kをインクリメントし（Ｓ４０
７）、ｇ_Kと２５５を比較する（Ｓ４０８）。ｇ_Kが２５
５を超えていなければステップＳ４０４からステップＳ
４０８の処理を繰り返す。ｇ_Kが２５５を超えていた場
合は、全候補拡散係数分の評価値演算が終了したかを確
認し（Ｓ４０９）、終了していなければステップＳ４０
２からステップＳ４０９の処理を繰り返す。全候補拡散
係数についての評価値演算が終了していれば、ｇ_Cと２
５５を比較する（Ｓ４１０）。ｇ_Cが２５５以下であれ
ば、ステップＳ４０１からステップＳ４１０の処理を繰
り返す。ｇ_Cの値が２５５を超えていれば、画像評価値
メモリ１７に格納された評価値合計をもとに、階調デー
タ値ｇ_C毎に全候補拡散係数の中から後述の方法により
１組の拡散係数を選択し、拡散係数テーブルに格納する
（Ｓ４１１）。

【００４５】尚、Ｋ成分の拡散係数と同様、階調データ
値ｇ_Cが０の場合の拡散係数はｇ_C＝１の際と同じ拡散係
数、又ｇ_Cの値が２５５の場合はｇ_C＝２５４における拡
散係数を拡散係数テーブル１３に格納する。

【００４６】Ｍ、Ｙ成分の拡散係数設定の手順は、Ｃ成
分の拡散係数設定の手順と同様である。但し、Ｍ成分の
拡散係数設定の際に用いる入力カラー画像は、Ｋ成分及
びＣ成分が２値化された画像であり、Ｙ成分の場合は
Ｋ、Ｃ、Ｍ成分が２値化された画像とする。

【００４７】（画像評価値演算ユニットの評価値演算
例）次に、画像評価値演算ユニット１６における評価値
演算方法について説明する。図１５は、本実施形態にお
ける画像の評価値演算処理の手順を示すフローチャート
である。

【００４８】まず、カラー画像を読み込む（Ｓ５０
０）。その後、各画素値を明度（Ｌ^*値）に変換する
（Ｓ５０１）。２色以上のドットが重なる場合は、該当
するインクを重ねた際のＬ^*値の測定値を使用する。こ
うして各画素値をＬ^*値に変換した画像データの２次元
フーリエ変換を行い、２次元パワースペクトルＰＷ
（ｉ，ｊ）を求める（Ｓ５０２）。ｉ，ｊは−１２８か
ら１２８の範囲の値であり、２次元パワースペクトル上
の位置を表す。図１６に、２次元パワースペクトルの一
例を示す。

【００４９】こうして求めた画像データの２次元パワー
スペクトルを１次元化する（Ｓ５０３）。図１７に示す
ように、２次元パワースペクトルＰＷ（ｉ，ｊ）を同心
円状の輪帯で区切り、周波数帯域ｆ毎のパワースペクト
ルの平均値を求める。

【００５０】即ち、

【００５１】

【数１】を求めるのである。ＩＮＴは小数点以下の切り上げ処理
を表す。周波数ｆは０から１８１までの整数、Ｎ（ｆ）
はｆ＝ＩＮＴ（（ｉ²＋ｊ²)^1/2）を満たす画素位置
（ｉ，ｊ）の数である。

【００５２】図１８は、横軸に周波数ｆ、縦軸に各輪帯
内のパワースペクトル平均値ＲＡＰＳ（ｆ）をとったグ
ラフの一例である。

【００５３】こうして１次元化したパワースペクトルＲ
ＡＰＳ（ｆ）に対して、視覚特性フィルタ（ＶＴＦ）処
理を行う（Ｓ４０４）。本実施の形態で用いる視覚特性
フィルタは、次の式で表される。

【００５４】

【数２】

【００５５】ここで、ｆ'は空間周波数(radial frequen
cy)であり、視野角１度あたりの周期（cycles/degree）
として与えられている。観察距離を３００ｍｍとし、画
像の解像度を１２００ｄｐｉに設定した場合、ｆ'を前
述したｆに変換すると、上式は次のようになる。

【００５６】

【数３】

【００５７】図１９は、ｆを横軸、ＶＴＦ（ｆ）を縦軸
にとったグラフである。

【００５８】こうして求めた視覚特性フィルターＶＴＦ
（ｆ）を用いて、ＲＡＰＳ（ｆ）に視覚特性フィルタ処
理を行い、Ｆｉｌｔｅｒｅｄ＿ＲＡＰＳ（ｆ）を求め
る。即ち、０から１８１までのｆの値について、次の演
算を行うのである。

【００５９】

【数４】

【００６０】こうして求めた視覚特性フィルタ処理を行
った１次元パワースペクトルＦｉｌｔｅｒｄ＿ＲＡＰＳ
（ｆ）の総和を求め（直流成分は除く）、画像評価値と
する（Ｓ５０５）。つまり、

【００６１】

【数５】で表される値を求め、画像評価値演算処理を終了する。

【００６２】ここで演算された画像評価値は、人の目に
付きやすい低周波ノイズが少ない２値画像に対応する値
は小さく、反対にノイズの多い画像に対応する値は大き
くなる傾向にある。このため、この画像評価値が小さく
なるような拡散係数を選択してハーフトーン処理に使用
することにより、ノイズが少なく粒状性の良い２値画像
を生成することが可能となる。すなわち、全画素が同一
の値を持つ入力画像に対し、１次元パワースペクトルを
求めた場合、周波数０におけるパワー（DC成分）以外の
パワーは存在しない。ところが、上記入力画像に対して
２値化を行い、該２値画像の１次元パワースペクトルを
求めると、DC成分以外のパワーが発生し、これは２値化
により生じるノイズと考えることができる。ここで更
に、人間の視覚の感度が周波数に応じて異なるという特
性を考慮するため、上記パワースペクトルに対して、視
覚特性フィルタ処理を行い、このようにフィルタ処理さ
れたパワースペクトルの総和を画像評価値とした。

【００６３】（拡散係数選択ユニットの拡散係数選択
例）次に、拡散係数選択ユニット１８における拡散係数
選択方法について説明する。

【００６４】図２０は、本実施の形態における色成分毎
の拡散係数選択の方法を示したフローチャートである。
図２０のｇは、処理する色成分に対応してｇ_K、ｇ_C、ｇ
_M、ｇ_Yとなる。

【００６５】まず、階調データ値ｇを１に初期化する
（Ｓ６００）。次に、画像評価値メモリ１７に格納され
た全候補拡散係数についての画像評価値を読み込み、同
評価値が最小となる候補拡散係数を拡散係数テーブル１
３のｇ番目に格納する（Ｓ６０１）。ｇの値をインクリ
メントし（Ｓ６０２）、ｇと２５５を比較する（Ｓ６０
３）。ｇが２５５よりも小さければステップＳ６０１か
らＳ６０３を繰り返す。ｇが２５５に達していれば、本
処理を終了する。

【００６６】このように、画像評価値が最小となるよう
な拡散係数を選択してハーフトーン処理に使用すること
により、ノイズが少なく粒状性の良い２値画像を生成す
ることが可能となる。

【００６７】このように、実施の形態１によれば、入力
階調毎の画像評価値を用いて入力階調毎に拡散係数を選
ぶことにより、各階調に適した２値化処理を行うことが
可能となる。又、既に拡散係数を設定した色成分との組
み合わせを考慮することにより、複数色が重なった際の
最適な２値化処理を行うことが可能となる。

【００６８】［実施の形態２］上記実施の形態１では、
拡散係数選択ユニット１８において色成分毎に各階調に
おける拡散係数を決定する際に、該当する１階調の画像
評価値のみを参照するものとしたが、実施の形態２では
複数の階調における画像評価値を参照する。

【００６９】（拡散係数選択ユニットの拡散係数選択の
他例）図２１は、本実施の形態２におけるある単色の拡
散係数選択方法を示したフローチャートである。図２１
のｇも、処理する色成分に対応してｇ_K、ｇ_C、ｇ_M、ｇ_Y
となる。

【００７０】まず、評価範囲を示す値ＩＮＴＶＬを設定
し（Ｓ７００）、階調データ値ｇを１に初期化する（Ｓ
７０１）。本実施の形態では、ＩＮＴＶＬを４とした。

【００７１】次に、画像評価値メモリ１７に格納された
画像評価値を読み込み、候補拡散係数毎に階調（ｇ−Ｉ
ＮＴＶＬ）から（ｇ＋ＩＮＴＶＬ）における画像評価値
の和Ｓを求める（Ｓ６０２）。このとき、ｇがＩＮＴＶ
Ｌよりも小さい場合は、階調１から（ｇ＋ＩＮＴＶＬ）
までの評価値の和をＳとする。ｇが（２５５−ＩＮＴＶ
Ｌ）よりも大きい場合は、階調（ｇ−ＩＮＴＶＬ）から
２５４までの評価値の和をＳとする。即ち、各候補係数
について

【００７２】

【数６】を求める。

【００７３】その後、Ｓが最小となる候補拡散係数を求
め、拡散係数テーブル１３のｇ番目に格納する（Ｓ６０
３）。ｇの値をインクリメントし（Ｓ６０４）、ｇと２
５５を比較する（Ｓ６０５）。ｇが２５５よりも小さけ
ればステップＳ６０２からステップＳ６０３の処理を繰
り返す。ｇが２５５に達していれば、本処理を終了す
る。

【００７４】このように、実施の形態２によれば、１階
調だけでなく近隣の階調における画像評価値も考慮する
ことにより、選ばれる係数が１階調ごとに大きく変化す
ることがなくなる。このため、階調変化が滑らかな画像
の２値化が可能となる。

【００７５】［実施の形態３］上記実施の形態１では、
画像２値化ユニット１２において入力カラー画像の色成
分毎の２値化処理を行う際、誤差バッファ２８の初期値
を０とした（図７のＳ１００）が、実施の形態３では同
初期値を乱数とする。

【００７６】（画像２値化ユニットでの２値化の他例）
図２２は、本実施の形態における入力カラー画像の一色
分の２値化処理手順を示したフローチャートである。

【００７７】まず、誤差バッファのＥ_temp及びＥ（０）
からＥ（Ｗ＋１）の初期値を−６４から６４の範囲内の
乱数（整数）に設定する（Ｓ８００）。その後のステッ
プＳ８０１からステップＳ８１１の処理は、図７のステ
ップＳ１０１からステップＳ１１１と同じであるため、
説明を省略する。

【００７８】１ライン分のデータの２値化処理が終了し
た後、−６４から６４のいずれかの整数をランダムに選
択し、Ｅ_tempに代入する（Ｓ８１２）。ステップＳ８１
３、ステップＳ８１４は、図７のステップＳ１１３、ス
テップＳ１１４と同じであり、ライン番号ｙが入力画像
の縦サイズＨに達した時点で２値化処理を終了する。

【００７９】実施の形態３によれば、誤差バッファの初
期値を０とした場合に生じる、ドットが横一列に規則的
に並ぶ現象を緩和することが可能となる。

【００８０】［実施の形態４］上記実施の形態１では、
Ｃ成分の各階調における拡散係数を設定する際に、Ｋ成
分の全階調における２値化結果との組み合わせを考慮し
た（図１３参照）が、実施の形態４では異なる色成分の
同一階調の組み合わせのみを考慮する。

【００８１】（Ｃ成分の拡散係数設定の他例）図２３
は、本実施の形態におけるＣ成分の拡散係数設定の手順
を示したフローチャートである。

【００８２】ステップＳ９００からステップＳ９０２
は、図１３のステップＳ４００からステップＳ４０２と
同じである。候補拡散係数設定後、入力カラー画像を作
成する（Ｓ９０３）。入力カラー画像のＫ成分は、全画
素がｇ_Kの値をもつ画像を、既に設定されたＫ成分の拡
散係数を用いて２値化したものとする。同Ｃ成分は実施
の形態１と同様、全画素値がｇ_Cである画像とする。そ
の後、この入力カラー画像のＣ成分の２値化を行う（Ｓ
９０４）。ステップＳ９０５からステップＳ９０８は、
図１３のステップＳ４０６，Ｓ４０９〜Ｓ４１１と同じ
である。同様に、Ｍ成分の拡散係数設定の際も、Ｃ、Ｋ
成分の同じ階調における２値化結果との組み合わせのみ
を考慮し、Ｙ成分の拡散係数設定の際は、Ｃ、Ｍ、Ｋ成
分の同じ階調における２値化結果との組み合わせのみを
考慮する。

【００８３】実施の形態４によれば、実施の形態１と比
較し、拡散係数テーブル作成に必要となる時間を短縮す
ることが可能となる。

【００８４】尚、上記の複数の実施の形態を組み合わせ
た処理を行うことも可能である。

【００８５】又、本発明は、複数の機器（例えばホスト
コンピュータ、インタフェイス機器、リーダ、プリンタ
など）から構成されるシステムに適用しても、一つの機
器からなる装置（例えば、複写機、ファクシミリ装置な
ど）に適用してもよい。

【００８６】又、本発明の目的は前述した実施形態の機
能を実現するソフトウェアのプログラムコードを記録し
た記憶媒体を、システムあるいは装置に供給し、そのシ
ステムあるいは装置のコンピュータ（又はＣＰＵやＭＰ
Ｕ）が記憶媒体に格納されたプログラム読出し実行する
ことによっても、本発明の目的が達成されることは言う
までもない。

【００８７】この場合、記憶媒体から読出されたプログ
ラムコード自体が本発明の新規な機能を実現することに
なり、そのプログラムコードを記憶した記憶媒体は本発
明を構成することになる。

【００８８】プログラムコードを供給するための記憶媒
体としては、例えば、フロッピー（登録商標）ディス
ク、ハードディスク、光ディスク、光磁気ディスク、Ｃ
Ｄ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリ
カード、ＲＯＭなどを用いることができる。又、コンピ
ュータが読出したプログラムコードを実行することによ
り、前述した実施形態の機能が実現されるだけでなく、
そのプログラムコードの指示に基づき、コンピュータ上
で稼動しているＯＳなどが実際の処理の一部又は全部を
行い、その処理によって前述した実施形態の機能が実現
される場合も含まれることは言うまでもない。

【００８９】さらに、記憶媒体から読出されたプログラ
ムコードが、コンピュータに挿入された機能拡張ボード
や機能拡張ユニットに備わるＣＰＵなどが実際の処理の
一部又は全部を行い、その処理によって前述した実施形
態の機能が実現される場合も含まれることは言うまでも
ない。

【００９０】本発明を上記記憶媒体に適用する場合、そ
の記憶媒体には、先に説明したフローチャートに対応す
るプログラムコードが格納されることになる。

【００９１】

【発明の効果】以上説明したように、本発明は各色成分
について入力階調毎に最適な拡散係数を用いることによ
り、全階調にわたって良好な２値化処理結果を得ること
ができる。又、異なる色成分間の相関を考慮することに
より、カラー画像に適した色成分毎の２値化を行うこと
ができる。

【図面の簡単な説明】

【図１】本発明の実施の形態である画像処理装置の構成
を示すブロック図である。

【図２】実施の形態１における拡散係数テーブルの一例
を示す図である。

【図３】画像２値化ユニット１２の構成を示すブロック
図である。

【図４】入力バッファ２２の構成例を示す図である。

【図５】出力バッファ２５の構成例を示す図である。

【図６】誤差バッファ２８の構成例を示す図である。

【図７】実施の形態１の２値化処理を示すフローチャー
トである。

【図８】誤差拡散の様子を示す説明図である。

【図９】誤差拡散の処理方向を逆にした様子を示す説明
図である。

【図１０】候補拡散係数を示す図である。

【図１１】拡散係数テーブル１３の作成手順を示すフロ
ーチャートである。

【図１２】ブラック成分の拡散係数テーブル作成手順の
詳細を示すフローチャートである。

【図１３】シアン成分の拡散係数テーブル作成手順の詳
細を示すフローチャートである。

【図１４】画像評価の対象とする部分を示す図である。

【図１５】画像の評価債演算処理を示すフローチャート
である。

【図１６】２次元パワースペクトルの一例を示した写真
である。

【図１７】２次元パワースペクトルの一次元化を説明す
る図である。

【図１８】１次元パワースペクトルの一例を示す図であ
る。

【図１９】視覚特性フィルタの一例を示す図である。

【図２０】実施の形態１における拡散係数選択を示すフ
ローチャートである。

【図２１】実施の形態２における拡散係数選択を示すフ
ローチャートである。

【図２２】実施の形態３の２値化処理を示すフローチャ
ートである。

【図２３】実施の形態４におけるシアン成分の拡散係数
テーブル作成手順を示すフローチャートである。

───────────────────────────────────────────────────── フロントページの続きＦターム(参考） 2C262 AA24 AB13 BB03 BB08 BC01 BC13 CA08 GA23 5B057 AA11 CA01 CA08 CA12 CA16 CB01 CB07 CB12 CB16 CC01 CE13 CE14 CH08 5C077 LL03 LL19 MP08 NN13 PP33 PP49 PQ08 PQ12 PQ23 SS02 TT02 TT06 5C079 HB03 LA12 LA14 LA31 LA34 LC09 MA04 MA11 NA05

Claims

【特許請求の範囲】

【請求項１】入力されたカラー多階調画像を、入力画
像の階調教より少ない階調数を有する画像に変換する擬
似中間調表現処理を行う画像処理装置であって、入力画像データの特徴に応じて、擬似中間調表現処理の
パラメータを変更することが可能な画像処理手段と、該画像処理手段で擬似中間調表現処理を行った出力画像
の評価値を求める評価値演算手段と、該評価値に基づいて、前記画像処理手段に設定される前
記擬似中間調表現処理のパラメータを選択するパラメー
タ選択手段とを有することを特徴とする画像処理装置。
【請求項２】前記擬似中間調表現処理は色成分毎の誤
差拡散法であって、前記パラメータは、色成分毎の誤差
拡散係数であることを特徴とする請求項１に記載の画像
処理装置。
【請求項３】前記入力画像データの特徴は、色成分毎
の入力画素値である請求項２に記載の画像処理装置。
【請求項４】前記画像処理手段は、色成分毎の入力画
素値に対応して誤差拡散係数を記憶するルックアップテ
ーブルを有し、前記パラメータ選択手段は、選択されたパラメータを、
色成分毎の入力画素値に対応して該ルックアップテーブ
ルに記憶することを特徴とする請求項３に記載の画像処
理装置。
【請求項５】前記誤差拡散法は、入力画素値と出力画
素値の誤差の初期値を０とすることを特徴とする請求項
２に記載の画像処理装置。
【請求項６】前記誤差拡散法は、入力画素値と出力画
素値の誤差の初期値を乱数とすることを特徴とする請求
項２に記載の画像処理装置。
【請求項７】前記パラメータ選択手段は、色成分毎の
各入力画素値に対応する前記評価値に基づいて、誤差拡
散係数を選択することを特徴とする請求項３に記載の画
像処理装置。
【請求項８】前記パラメータ選択手段は、色成分毎の
複数の入力画素値に対応する前記評価値に基づいて誤差
拡散係数を選択することを特徴とする請求項３に記載の
画像処理装置。
【請求項９】前記パラメータ選択手段は、異なる色成
分間の全ての画素値の組み合わせにおける評価値に基づ
いて各色成分の誤差拡散係数を選択することを特徴とす
る請求項７又は８に記載の画像処理装置。
【請求項１０】前記パラメータ選択手段は、異なる色
成分間の同一の画素値における評価値に基づいて各色成
分の誤差拡散係数を選択することを特徴とする請求項７
又は８に記載の画像処理装置。
【請求項１１】前記評価値演算手段は、全面が所定の
画素値の入力画像に擬似中間調表現処理を行った出力画
像を周波数領域へ変換する変換手段と、周波数領域にお
いて評価値を演算する演算手段とを含むことを特徴とす
る請求項３に記載の画像処理装置。
【請求項１２】前記変換手段は、前記出力画像の一部
からパワースペクトルを求め、前記演算手段は、前記パ
ワースペクトルをフィルタ処理しながら、所定の周波数
範囲にわたって加算することを特徴とする請求項１１に
記載の画像処理装置。
【請求項１３】前記変換手段は、前記出力画像の一部
から２次元フーリエ変換により２次元パワースペクトル
を求め、更に、該２次元パワースペクトルから周波数帯
域毎に平均値を取った１次元パワースペクトルを求め、
前記演算手段は、前記１次元パワースペクトルを視覚特
性フィルタ処理しながら、所定の周波数範囲にわたって
加算することを特徴とする請求項１２に記載の画像処理
装置。
【請求項１４】前記パラメータ選択手段は、前記演算
手段の加算値が最小の誤差拡散係数を選択することを特
徴とする請求項１２又は１３に記載の画像処理装置。
【請求項１５】入力されたカラー多階調画像を、入力
画像の階調数より少ない階調数を有する画像に変換する
擬似中間調表現処理を行う画像処理方法であって、擬似中間調表現処理を行った画像データの評価値を求め
る評価値演算工程と、該評価値に基づいて前記擬似中間調表現処理のパラメー
タを選択するパラメータ選択工程と、入力画像データの特徴に応じて、選択された擬似中間調
表現処理のパラメータを使って擬似中間調表現処理を行
う画像処理工程とを有することを特徴とする画像処理方
法。
【請求項１６】前記擬似中間調表現処理は色成分毎の
誤差拡散法であって、前記パラメータは、色成分毎の誤
差拡散係数であることを特徴とする請求項１５に記載の
画像処理方法。
【請求項１７】前記入力画像データの特徴は、色成分
毎の入力画素値である請求項１６に記載の画像処理方
法。
【請求項１８】前記誤差拡散法は、入力画素値と出力
画素値の誤差の初期値を０とすることを特徴とする請求
項１６に記載の画像処理方法。
【請求項１９】前記誤差拡散法は、入力画素値と出力
画素値の誤差の初期値を乱数とすることを特徴とする請
求項１６に記載の画像処理方法。
【請求項２０】前記パラメータ選択工程では、色成分
毎の各入力画素値に対応する前記評価値に基づいて、誤
差拡散係数を選択することを特徴とする請求項１７に記
載の画像処理方法。
【請求項２１】前記パラメータ選択工程では、色成分
毎の複数の入力画素値に対応する前記評価値に基づいて
誤差拡散係数を選択することを特徴とする請求項１７に
記載の画像処理方法。
【請求項２２】前記パラメータ選択工程では、異なる
色成分間の全ての画素値の組み合わせにおける評価値に
基づいて各色成分の誤差拡散係数を選択することを特徴
とする請求項２０又は２１に記載の画像処理方法。
【請求項２３】前記パラメータ選択工程では、異なる
色成分間の同一の画素値における評価値に基づいて各色
成分の誤差拡散係数を選択することを特徴とする請求項
２０又は２１に記載の画像処理方法。
【請求項２４】前記評価値演算工程は、前面が所定の
画素値の入力画像に擬似中間調表現処理を行った出力画
像を周波数領域へ変換する変換工程と、周波数領域にお
いて評価値を演算する演算工程とを含むことを特徴とす
る請求項１７に記載の画像処理方法。
【請求項２５】前記変換工程では、前記出力画像の一
部からパワースペクトルを求め、前記演算工程では、前
記パワースペクトルをフィルタ処理しながら、所定の周
波数範囲にわたって加算することを特徴とする請求項２
４に記載の画像処理方法。
【請求項２６】前記変換工程では、前記出力画像の一
部から２次元フーリエ変換により２次元パワースペクト
ルを求め、更に、該２次元パワースペクトルから周波数
帯域毎に平均値を取った１次元パワースペクトルを求
め、前記演算工程では、前記１次元パワースペクトルを
視覚特性フィルタ処理しながら、所定の周波数範囲にわ
たって加算することを特徴とする請求項２５に記載の画
像処理方法。
【請求項２７】前記パラメータ選択工程は、前記演算
手段の加算値が最小の誤差拡散係数を選択することを特
徴とする請求項２５又は２６に記載の画像処理方法。
【請求項２８】入力されたカラー多階調画像を、入力
画像の階調数より少ない階調数を有する画像に変換する
擬似中間調表現処理を行うためのプログラムを記憶する
コンピュータ読取り可能な記憶媒体であって、前記プログラムは、少なくとも、擬似中間調表現処理を行った画像データの評価値を求め
る評価値演算ステップと、該評価値に基づいて前記擬似中間調表現処理のパラメー
タを選択するパラメータ選択ステップとを有することを
特徴とする記憶媒体。
【請求項２９】前記擬似中間調表現処理は色成分毎の
誤差拡散法であって、前記パラメータは、色成分毎の誤
差拡散係数であることを特徴とする請求項２８に記載の
記憶媒体。
【請求項３０】前記評価値演算ステップは、前面が所
定の画素値の入力画像に擬似中間調表現処理を行った出
力画像を周波数領域へ変換する変換ステップと、周波数
領域において評価値を演算する演算ステップとを含むこ
とを特徴とする請求項２９に記載の記憶媒体。
【請求項３１】前記変換ステップでは、前記出力画像
の一部から２次元フーリエ変換により２次元パワースペ
クトルを求め、更に、該２次元パワースペクトルから周
波数帯域毎に平均値を取った１次元パワースペクトルを
求め、前記演算ステップでは、前記１次元パワースペク
トルを視覚特性フィルタ処理しながら、所定の周波数範
囲にわたって加算することを特徴とする請求項３０に記
載の記憶媒体。
【請求項３２】前記パラメータ選択ステップでは、前
記演算手段の加算値が最小の誤差拡散係数を選択するこ
とを特徴とする請求項３０又は３１に記載の記憶媒体。
【請求項３３】前記パラメータ選択ステップでは、色
成分毎の各入力画素値に対応する前記評価値に基づい
て、誤差拡散係数を選択することを特徴とする請求項２
９に記載の記憶媒体。
【請求項３４】前記パラメータ選択ステップでは、色
成分毎の複数の入力画素値に対応する前記評価値に基づ
いて誤差拡散係数を選択することを特徴とする請求項２
９に記載の記憶媒体。
【請求項３５】前記パラメータ選択ステップでは、異
なる色成分間の全ての画素値の組み合わせにおける評価
値に基づいて各色成分の誤差拡散係数を選択することを
特徴とする請求項３３又は３４に記載の記憶媒体。
【請求項３６】前記パラメータ選択ステップでは、異
なる色成分間の同一の画素値における評価値に基づいて
各色成分の誤差拡散係数を選択することを特徴とする請
求項３３又は３４に記載の記憶媒体。
【請求項３７】入力画像データの特徴に応じて、選択
された擬似中間調表現処理のパラメータを使って擬似中
間調表現処理を行う画像処理ステップを更に有すること
を特徴とする記憶媒体。
【請求項３８】入力されたカラー多階調画像を、入力
画像の階調数より少ない階調数を有する画像に変換する
擬似中間調表現処理を行うためのデータを記憶するコン
ピュータ読取り可能な記憶媒体であって、前記データは、擬似中間調表現処理を行った画像データ
の評価値を求め、該評価値に基づいて選択された前記擬
似中間調表現処理のパラメータを含むことを特徴とする
記憶媒体。