(57)【要約】【課題】バンディングを抑制しながらテクスチャの無い一様な画像を出力することが可能な量子化処理方法を提供する。【解決手段】処理対象画素で発生した量子化誤差ΔＥを、処理対象画素の周辺に位置する量子化が行われていない周辺画素ｉに拡散するための拡散係数として、第１拡散係数ｈ１＿ｉと、第２拡散係数ｈ２＿ｉとを用意する。量子化誤差ΔＥは、処理対象画素（ｘ、ｙ）に対応づけて決定されたブレンド比率αに基づいて、第１拡散係数ｈ１＿ｉと第２拡散係数ｈ２＿ｉとに従って前記周辺画素に拡散される。【選択図】図４
2020124835.tif 000002
【特許請求の範囲】
【請求項１】
画素が有する画素値を、誤差拡散処理によって量子化するための画像処理装置であって、
処理対象画素の画素値を量子化する際に発生した量子化誤差を取得する手段と、
前記処理対象画素の周辺に位置する量子化が行われていない周辺画素に前記量子化誤差を拡散するための第１拡散係数と、前記第１拡散係数とは異なる第２拡散係数を保持する保持手段と、
前記処理対象画素に対応づけて、前記第１拡散係数と前記第２拡散係数とのブレンド比率を決定するブレンド比率決定手段と、
前記ブレンド比率に基づき、前記量子化誤差を、前記第１拡散係数と前記第２拡散係数とに従って前記周辺画素に拡散する拡散手段と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記ブレンド比率決定手段は、前記ブレンド比率を、画素が配列する方向に対し連続性を有するように、前記処理対象画素の画素位置に対応づけて決定することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記画素が配列する方向は、前記誤差拡散処理によって生成された記録データに従って画像を記録する記録ヘッドにおいて複数の記録素子が配列する方向であり、
前記ブレンド比率決定手段は、前記ブレンド比率を、前記複数の記録素子のぞれぞれの記録特性に基づいて決定することを特徴とする請求項２に記載の画像処理装置。
【請求項４】
前記複数の記録素子のぞれぞれの記録特性を記憶する記憶手段を更に備え、
前記記録特性は、前記記録ヘッドが記録したテストチャートを読み取って得られる画像データに基づいて導出されることを特徴とする請求項３に記載の画像処理装置。
【請求項５】
前記ブレンド比率は、前記画素が配列する方向に対し、低周波成分と部分的な高周波成分とを有するように連続的に変化することを特徴とする請求項２から４のいずれか１項に記載の画像処理装置。
【請求項６】
前記ブレンド比率決定手段は、前記ブレンド比率を、画素値に対し連続性を有するように、前記処理対象画素の画素値に対応づけて決定することを特徴とする請求項１に記載の画像処理装置。
【請求項７】
前記拡散手段は、前記ブレンド比率に基づいて前記第１拡散係数と前記第２拡散係数とをブレンドすることによってブレンド拡散係数を導出し、当該ブレンド拡散係数に従って前記量子化誤差を前記周辺画素に拡散することを特徴とする請求項１から６のいずれか１項に記載の画像処理装置。
【請求項８】
前記拡散手段は、前記周辺画素において、前記ブレンド拡散係数の総和が、前記第１拡散係数の総和および前記第２拡散係数の総和のそれぞれと等しくなるように、前記ブレンド拡散係数を導出することを特徴とする請求項７に記載の画像処理装置。
【請求項９】
前記拡散手段は、前記第１拡散係数に従う前記量子化誤差の第１拡散値と、前記第２拡散係数に従う前記量子化誤差の第２拡散値とを、前記ブレンド比率に基づいてブレンドすることにより、前記周辺画素に拡散する拡散値を導出することを特徴とする請求項１から６のいずれか１項に記載の画像処理装置。
【請求項１０】
前記拡散手段は、前記周辺画素における前記拡散値の総和が前記量子化誤差と等しくなるように、前記第１拡散値と前記第２拡散値とをブレンドすることを特徴とする請求項９に記載の画像処理装置。
【請求項１１】
前記第１拡散係数のみに従って前記量子化誤差を前記周辺画素に拡散した場合には、ブルーノイズ特性を有するドットパターンが得られ、
前記第２拡散係数のみに従って前記量子化誤差を前記周辺画素に拡散した場合には、グリーンノイズ特性を有するドットパターンが得られることを特徴とする請求項１から１０のいずれか１項に記載の画像処理装置。
【請求項１２】
前記量子化誤差を、前記周辺画素に拡散するための、前記第１拡散係数および前記第２拡散係数とは異なる第３拡散係数を保持する手段を更に備え、
前記ブレンド比率決定手段は、前記処理対象画素に対応づけて、前記第１拡散係数、前記第２拡散係数、および前記第３拡散係数のブレンド比率を決定し、
前記拡散手段は、前記ブレンド比率に基づき、前記量子化誤差を、前記第１拡散係数と前記第２拡散係数と前記第３拡散係数とに従って前記周辺画素に拡散する
ことを特徴とする請求項１から１１のいずれか１項に記載の画像処理装置。
【請求項１３】
画素が有する画素値を、誤差拡散処理によって量子化するための画像処理方法であって、
処理対象画素の画素値を量子化する際に発生した量子化誤差を取得する工程と、
前記処理対象画素の周辺に位置する量子化が行われていない周辺画素に前記量子化誤差を拡散するための第１拡散係数と、前記第１拡散係数とは異なる第２拡散係数を保持する保持工程と、
前記処理対象画素に対応づけて、前記第１拡散係数と前記第２拡散係数とのブレンド比率を決定するブレンド比率決定工程と、
前記ブレンド比率に基づき、前記量子化誤差を、前記第１拡散係数と前記第２拡散係数とに従って前記周辺画素に拡散する拡散工程と、
を有することを特徴とする画像処理方法。
【請求項１４】
コンピュータを請求項１から１２のいずれか１項に記載の画像処理装置の各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法、及びプログラムに関する。特に、画素が有する画素値を量子化するための量子化方法に関する。
【背景技術】
【０００２】
インクジェット記録装置においては、記録ヘッドに配列する複数のノズルのそれぞれが記録データに従ってインクを吐出し、吐出されたインクが記録媒体に着弾することによって、画像が記録される。個々のノズルにおいては、吐出量や吐出方向などの吐出特性に、どうしてもある程度のばらつきが含まれるが、このばらつきが大きいと、画像上にスジや濃度ムラが確認されることがある。
【０００３】
特に、記録ヘッドと記録媒体との１回の相対走査によって画像を完成させるシングルパス記録では、各ノズルが記録するドットが記録媒体上で一列に配列するため、上記濃度ムラは目立ちやすい傾向にある。加えて、シングルパス記録では、隣接するドットの着弾時間差が短く、記録媒体に吸収される前のインク滴同士が接触および統合しやすいため、上記濃度ムラはさらに強調される。以下、このような記録ヘッドや記録装置固有の濃度ムラを、本明細書ではバンディングと称する。
【０００４】
特許文献１には、バンディングを目立たなくするための量子化方法が開示されている。具体的には、バンディング現象に関与するノズルと、関与しないノズルとで個別の誤差拡散マトリクスを用意し、画素ごとにこれら２つの誤差拡散マトリクスを切り替えながら誤差拡散処理を行っている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００６−２３１９０３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１の方法では、通常の誤差拡散マトリクスと、バンディングを抑えるための誤差拡散マトリクスとが、画素単位で切り替えられながら用いられる。そのため、記録媒体上では、誤差拡散マトリクスの切り替わりに伴うテクスチャが目立ってしまう場合がある。
【０００７】
本発明は上記問題点を解消するためになされたものである。よってその目的とするところは、インクジェット方式で画像を記録する構成において、バンディングを抑制した高品位な画像を出力することが可能な量子化処理方法を提供することである。
【課題を解決するための手段】
【０００８】
そのために本発明は、画素が有する画素値を、誤差拡散処理によって量子化するための画像処理装置であって、処理対象画素の画素値を量子化する際に発生した量子化誤差を取得する手段と、前記処理対象画素の周辺に位置する量子化が行われていない周辺画素に前記量子化誤差を拡散するための第１拡散係数と、前記第１拡散係数とは異なる第２拡散係数を保持する保持手段と、前記処理対象画素に対応づけて、前記第１拡散係数と前記第２拡散係数とのブレンド比率を決定するブレンド比率決定手段と、前記ブレンド比率に基づき、前記量子化誤差を、前記第１拡散係数と前記第２拡散係数とに従って前記周辺画素に拡散する拡散手段と、を備えることを特徴とする。
【発明の効果】
【０００９】
本発明によれば、バンディングを抑制した高品位な画像を出力することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】画像処理システムの一例を示すブロック図である。
【図２】記録装置と記録ヘッドの概略構成図である。
【図３】画像処理の工程を説明するためのフローチャートである。
【図４】第１の実施形態の量子化処理部の機能構成を示すブロック図である。
【図５】第１の実施形態の量子化処理の工程を説明するためのフローチャートである。
【図６】吐出特性情報Ｎ（ｙ）とブレンド比率αの関係を示す図である。
【図７】第１拡散係数と第２拡散係数のマトリクスを示す図である。
【図８】吐出特性情報更新モードを説明するためのフローチャートである。
【図９】位置ズレ量Ｘ（ｙ）とフィルタ処理を施した結果を示す。
【図１０】第１拡散係数および第２拡散係数と、ドットパターンを示す図である。
【図１１】ブレンド比率αを様々に振った場合のドットパターンを比較する図である。
【図１２】ブレンド比率とドットパターンの関係を示す図である。
【図１３】ブレンド比率とドットパターンの関係を示す図である。
【図１４】第２の実施形態の量子化処理部の機能構成を示すブロック図である。
【図１５】第２の実施形態の量子化処理の工程を説明するためのフローチャートである。
【図１６】第３の実施形態の量子化処理部の機能構成を示すブロック図である。
【図１７】第３の実施形態の量子化処理の工程を説明するためのフローチャートである。
【図１８】第４の実施形態の量子化処理部の機能構成を示すブロック図である。
【図１９】第４の実施形態の量子化処理の工程を説明するためのフローチャートである。
【図２０】階調値Ｉに対するブレンド比率αの例を示す図である。
【図２１】階調値Ｉとブレンド比率αの対応関係を求めるためのフローチャートである。
【図２２】テストチャートの一例を示す図である。
【図２３】演算誤差を説明するための模式図である。
【図２４】補正処理を説明するためのフローチャートである。
【発明を実施するための形態】
【００１１】
以下、本発明の実施形態について、図面を参照しながら詳細に説明する。  
【００１２】
（第１の実施形態）
図１は、第１の実施形態で使用可能な画像処理システムの一例を示すブロック図である。カラースキャナ、デジタルカメラなどの入力装置１０００で取得された画像データは、パーソナルコンピュータのような画像処理装置１１００で所定の画像処理が施された後、記録装置１２００に出力される。
【００１３】
本実施形態において、入力装置１０００から画像処理装置１１００に入力される画像データは、ＲＧＢの８ｂｉｔ（２５６階調）データとする。ユーザは、モニタ１１３０で表示されるカラー画像を確認しながら、キーボード１１４０、マウス１１６０を用いて、画像のレイアウトや色補正などを行い、修正後の画像に対し、記録コマンドを発生する。記録コマンドが発生すると、画像処理装置１１００のＭＰＵ１１１０は、修正後の画像データに対し、ＲＯＭ１１２０に記憶されているプログラムに従って所定の画像処理を施し、記録装置１２００が記録可能な２値の記録データを生成する。この際、量子化処理部１は、８ビットの画像データを、記録装置１２００が記録可能な１ビット２値の記録データに変換する量子化処理を行う。
【００１４】
記録装置１２００は、ＭＰＵ１２０１、ＲＡＭ１２０２などを有している。ＭＰＵ１２０１は、画像処理装置１１００より受信した記録データをＲＡＭ１２０２に展開し、各種機構を制御して記録媒体に画像を記録する。
【００１５】
図２（ａ）および（ｂ）は、本実施形態で使用する記録装置１２００と記録ヘッド２０の概略構成図である。本実施形態の記録装置１２００は、フルラインタイプのインクジェット記録装置である。
【００１６】
図２（ａ）に示すように、装置内に固定されている記録ヘッド２０には、ブラック、シアン、マゼンタおよびイエローのインクをそれぞれ吐出する４つのノズル列（２０Ｋ、２０Ｃ、２０Ｍ、２０Ｙ）がｘ方向に並列配置されている。それぞれのノズル列（２０Ｋ、２０Ｃ、２０Ｍ、２０Ｙ）には、同種類のインクを吐出する複数のノズル１２０３が、記録媒体Ｓの幅に相当する距離だけ１２００ｄｐｉのピッチでｙ方向に配列されている。記録ヘッド２０と対向する位置には平板から成るプラテン１２０４が配され、記録中の記録媒体Ｓを背面から支持する。
【００１７】
以上の構成のもと、搬送ローラ１２０３によって所定の速度でｘ方向に搬送される記録媒体Ｓに対し、記録ヘッド２０が記録データに従ってインクを吐出することにより、記録媒体Ｓに所望の画像が記録される。
【００１８】
図３は記録コマンドが発生した際に画像処理装置１１００のＭＰＵ１１１０が実行する画像処理の工程を説明するためのフローチャートである。本処理は、ＭＰＵ１１１０が、ＲＯＭ１１２０に記憶されたプログラムに従って量子化処理部１などを利用しながら実行するものである（図１参照）。以下の説明において、各工程（ステップ）を符号の前に「Ｓ」と表記する。
【００１９】
本処理が開始されると、ＭＰＵ１１１０は、まず、記録すべき画像データを取得する。
本実施形態では、個々の画素についてレッド（Ｒ）、グリーン（Ｇ）、ブルー（Ｂ）それぞれの画素値を８ｂｉｔ（２５６階調）で表す画像データを取得する。
【００２０】
Ｓ３０１において、ＭＰＵ１１１０は取得した画像データに対し、色補正処理を実行する。色補正処理とは、ｓＲＧＢ等の規格化された色空間を、記録装置１２００が記録可能な色空間に対応づけるための処理である。具体的には、ＭＰＵ１１１０は、予め記憶されたＲＧＢの３次元ルックアップテーブルを参照し、８ｂｉｔ（２５６階調）のＲＧＢデータを、同じく８ｂｉｔ（２５６階調）のＲ´Ｇ´Ｂ´データに変換する。
【００２１】
Ｓ３０２において、ＭＰＵ１１１０はＳ３０１で生成したＲ´Ｇ´Ｂ´データに対し、色変換処理を実行する。色変換処理とは、輝度情報を示すＲ´Ｇ´Ｂ´データを、記録装置１２００が使用するインク色に対応した、Ｋ、Ｃ、Ｍ、Ｙのそれぞれの階調値を示す階調データに変換するための処理である。具体的には、ＭＰＵ１１１０は、予め記録された３次元ルックアップテーブルを参照し、各画素の８ｂｉｔ（２５６階調）のＲ´Ｇ´Ｂ´データを、８ｂｉｔ（２５６階調）のＫＣＭＹデータに変換する。
【００２２】
Ｓ３０３において、ＭＰＵ１１１０はＳ３０２で生成したＫＣＭＹデータに対し、量子化処理を行う。量子化処理については後に詳しく説明するが、この量子化処理によって、８ｂｉｔ（２５６階調）のＫＣＭＹデータは、１ｂｉｔのＫ´Ｃ´Ｍ´Ｙ´データに変換される。本実施形態において、Ｋ´Ｃ´Ｍ´Ｙ´データは、ＫＣＭＹの各インク色について、ドットの記録（１）又は非記録（０）を画素ごとに指定する２値データである。
【００２３】
Ｓ３０４において、ＭＰＵ１１１０はＳ３０３で生成したＫ´Ｃ´Ｍ´Ｙ´データを記録データとして記録装置１２００に出力する。以上で本処理が終了する。
【００２４】
記録データを受信した記録装置１２００のＭＰＵ１２０１は、記録データをＲＡＭ１２０２に展開する。そして、搬送ローラ１２０３を所定の速度で回転させながら、ＲＡＭ１２０２に展開した記録データに従って記録ヘッド２０を駆動し、記録媒体に画像を記録する。
【００２５】
図４は、本実施形態の量子化処理部１の機能構成を説明するためのブロック図である。本実施形態において、図３のＳ３０３で実行される量子化処理は、ＭＰＵ１１１０の指示のもと量子化処理部１によって実行されるものである。
【００２６】
量子化処理部１は、ＫＣＭＹデータの各色（Ｋデータ、Ｃデータ、Ｍデータ、Ｙデータ）に対し、同じ処理を並行して行う。よって、以下ではブラックの階調データ（Ｋデータ）の量子化処理部１について説明する。
【００２７】
色変換処理によって生成された階調データＩ（ｘ、ｙ）は、入力部２より量子化処理部１に入力される。ここで、（ｘ、ｙ）は処理対象となる画素の記録媒体上での画素位置（ｘｙ座標）を示し、Ｉ（ｘ、ｙ）は画素位置（ｘ、ｙ）の画素が有する２５６階調の画素値（以下、階調値と称す）を示す。
【００２８】
加算器３は、入力部２より入力された階調値Ｉ（ｘ、ｙ）に、誤差拡散値Ｅ（ｘ、ｙ）を加算し比較値Ｉ´（ｘ、ｙ）を生成する。比較器４は、比較値Ｉ´（ｘ、ｙ）を予め定められた閾値Ｔと比較し、量子化代表値Ｒ（ｘ、ｙ）を出力する。出力部５は、量子化代表値Ｒ（ｘ、ｙ）に基づいて、「１」または「０」の記録データＯ（ｘ、ｙ）を出力する。
【００２９】
減算器６は、比較値Ｉ´（ｘ、ｙ）と、量子化代表値Ｒ（ｘ、ｙ）との差分を算出し、画素位置（ｘ、ｙ）における量子化誤差ΔＥとして、拡散処理部１５に提供する。拡散処理部１５は、画素位置（ｘ、ｙ）に対応づけられたブレンド拡散係数Ｈｉに従って、量子化誤差ΔＥを周辺の画素に拡散する。
【００３０】
ブレンド比率決定部７は、吐出特性情報格納部８を参照し、処理対象画素の画素位置（ｘ、ｙ）に対応するブレンド比率αを出力する。
【００３１】
減算器１０は、１．０よりブレンド比率αを減算した値（１．０−α）を第１乗算器９に出力する。第１乗算器９は、第１拡散係数格納部１１に保持されている拡散係数ｈ１＿ｉと減算器１０より受信した値（１．０−α）を乗算する。ここで、ｉは、処理対象画素の周囲に位置する４つの画素の、処理対象画素に対する相対位置を示す。一方、第２乗算器１２は、第２拡散係数格納部１３に保持されている拡散係数ｈ２＿ｉとブレンド比率αを乗算する。
【００３２】
加算器１４は、第１乗算器９の出力値Ｈ１＿ｉと、第２乗算器１２の出力値Ｈ２＿ｉとを加算し、得られた結果をブレンド拡散係数Ｈｉとして、拡散処理部１５に提供する。
【００３３】
なお、図４では、量子化処理部１をハードウェアとして示しているが、量子化処理部１の各機能は、ＭＰＵ１１１０がＲＯＭ１１２０に格納されたプログラムに従って実現する形態であってもよい。
【００３４】
図５は、図４で説明した量子化処理部１を用いて、画像処理装置１１００のＭＰＵ１１１０が実行する量子化処理の工程を説明するためのフローチャートである。以下、図４のブロック図を参照しながら、図５の各工程について説明する。
【００３５】
本処理が開始されると、ＭＰＵ１１１０は、まずＳ１００において、１ページに含まれる全画素についての階調データを取得する。
【００３６】
Ｓ１０１において、ＭＰＵ１１１０は、全画素位置の中から処理対象画素（ｘ、ｙ）を設定し、続くＳ１０２において、処理対象画素（ｘ、ｙ）の階調値Ｉ（ｘ、ｙ）を取得する。
【００３７】
Ｓ１０３において、ＭＰＵ１１１０は、処理対象画素の階調値Ｉ（ｘ、ｙ）に処理対象画素の誤差拡散値Ｅ（ｘ、ｙ）を加算し、比較値Ｉ´（ｘ、ｙ）を得る。
Ｉ´（ｘ、ｙ）＝Ｉ（ｘ、ｙ）＋Ｅ（ｘ、ｙ）
ここで、処理対象画素が最初の画素である場合、誤差拡散値Ｅ（ｘ、ｙ）は０であるのでそのままＳ１０４に進む。
【００３８】
Ｓ１０４において、ＭＰＵ１１１０は、予め定められている閾値Ｔと比較値Ｉ´（ｘ、ｙ）とを比較し、量子化代表値Ｒ（ｘ、ｙ）を取得する。閾値Ｔについては、固定値であってもよいが、画素ごとに異なる値が設定されてもよい。ここで、例えば閾値をＴ＝１２８とすると、
Ｉ´（ｘ、ｙ）≧１２８  のとき  Ｒ（ｘ、ｙ）＝２５５
Ｉ´（ｘ、ｙ）＜１２８  のとき  Ｒ（ｘ、ｙ）＝０
となる。更に、ＭＰＵ１１１０は、得られた量子化代表値Ｒ（ｘ、ｙ）に従って２値の記録データＯ（ｘ、ｙ）を求め、記録装置１２００に出力する。
Ｒ（ｘ、ｙ）＝２５５  のとき  Ｏ（ｘ、ｙ）＝１
Ｒ（ｘ、ｙ）＝０      のとき  Ｏ（ｘ、ｙ）＝０
【００３９】
Ｓ１０５において、ＭＰＵ１１１０は、比較値Ｉ´（ｘ、ｙ）と量子化代表値Ｒ（ｘ、ｙ）との差分を求め、これを処理対象画素の量子化誤差ΔＥとする。
ΔＥ＝Ｉ´（ｘ、ｙ）−Ｒ（ｘ、ｙ）
【００４０】
Ｓ１０６において、ＭＰＵ１１１０は、吐出特性情報格納部８をアクセスし、処理対象画素（ｘ、ｙ）に対応づけられた吐出特性情報Ｎ（ｙ）を取得する。本実施形態において、吐出特性情報Ｎ（ｙ）とは、個々のノズル（ｙ方向の画素位置）に対応づけて予め記憶されている、スジムラの程度を示す値である。本実施形態の吐出特性情報Ｎ（ｙ）は０〜５に規格化されており、値が大きいほど±ｙ方向へのずれが大きく、画像上でスジムラが目立ち易いことを意味している。
【００４１】
Ｓ１０７にて、ＭＰＵ１１１０は、取得した吐出特性情報Ｎ（ｙ）に基づいて、処理対象画素（ｘ、ｙ）のブレンド比率αを設定する。
【００４２】
図６は、吐出特性情報Ｎ（ｙ）とブレンド比率αの関係を示す図である。ブレンド比率αは０．０〜１．０の値を有する実数であり、ここでは、０〜５の吐出特性情報Ｎ（ｙ）に対し線形関係を有するものとする。すなわち、ブレンド比率αは、定数ｋを用い、以下の式で求めることができる。
α＝ｋ×Ｎ（ｙ）
【００４３】
但し、吐出特性情報Ｎ（ｙ）とブレンド比率αは、必ずしも線形関係を有していなくてもよい。吐出特性情報Ｎ（ｙ）とブレンド比率αは、互いに１対１で対応づけられていればよい。
【００４４】
図５のフローチャートに戻る。Ｓ１０７でブレンド比率αが設定されると、ＭＰＵ１１１０は、Ｓ１０８において、第１拡散係数格納部１１より第１拡散係数ｈ１＿ｉを読み出し、ブレンド比率αに従って第１係数Ｈ１＿ｉを求める。また、Ｓ１０９において、第２拡散係数格納部１３より第２拡散係数ｈ２＿ｉを読み出し、ブレンド比率αに従って第２係数Ｈ２＿ｉを求める。
【００４５】
図７（ａ）および（ｂ）は、第１拡散係数格納部１１および第２拡散係数格納部１３のそれぞれに格納されている第１拡散係数ｈ１＿ｉと第２拡散係数ｈ２＿ｉのマトリクスを示す図である。拡散係数とは、誤差拡散処理を行う際に、処理対象画素（＊）で発生した量子化誤差ΔＥを、量子化処理（２値化処理）が未処理の周辺画素ｉに分配するための係数である。処理対象画素（＊）に対する未処理の周辺画素の位置は、量子化処理を進めていく方向によって異なるため、拡散係数は量子化処理を進めていく方向に対応づけて用意される。
【００４６】
例えば、図７（ａ）は、量子化処理を＋ｙ方向に進めて行く場合の、処理対象画素（＊）と周辺画素ｉの相対位置、及び各周辺画素ｉにおける第１の分散係数ｈ１＿ｉと第２の分散係数ｈ２＿ｉを示している。また、図７（ｂ）は、量子化処理を−ｙ方向に進めて行く場合の、処理対象画素（＊）と周辺画素ｉの相対位置、及び各周辺画素ｉにおける第１の分散係数ｈ１＿ｉと第２の分散係数ｈ２＿ｉを示している。量子化処理は、図７（ａ）で示す方向で行ってもよいし、図７（ｂ）で示す方向で行ってもよいし、図７（ａ）で示す方向と図７（ｂ）で示す方向とを、ラインごとに交互に切替えながら行ってもよい。
【００４７】
更に加えると、図では示していないが、＋ｘ方向や−ｘ方向に量子化処理を進めて行ってもよい。いずれにせよ、図７（ａ）および（ｂ）で示すような拡散係数のマトリクスが、量子化処理が行われる方向に対応づけて用意されていればよい。
【００４８】
以下では、簡単のため、量子化処理は＋ｙ方向に進められ、Ｓ１０８及びＳ１０９では、図７（ａ）に示す分散係数が読み出されるものとする。具体的に説明すると、図５のＳ１０８において、ＭＰＵ１１１０は、分散係数ｈ１＿１、ｈ１＿２、ｈ１＿３、ｈ１＿４を読み出し、それぞれに（１−α）を乗算して４つの第１係数、Ｈ１＿１、Ｈ１＿２、Ｈ１＿３、Ｈ１＿４を得る。
Ｈ１＿ｉ＝（１−α）×ｈ１＿ｉ
【００４９】
また、Ｓ１０９において、ＭＰＵ１１１０は、分散係数ｈ２＿１、ｈ２＿２、ｈ２＿３、ｈ２＿４を読み出し、それぞれにαを乗算して４つの第２係数、Ｈ２＿１、Ｈ２＿２、Ｈ２＿３、Ｈ２＿４を得る。
Ｈ２＿ｉ＝α×ｈ２＿ｉ
【００５０】
Ｓ１１０において、ＭＰＵ１１１０は、Ｓ１０８及びＳ１０９で求めた第１係数Ｈ１＿ｉと第２係数Ｈ２＿ｉを加算して、ブレンド拡散係数Ｈｉを得る。
Ｈｉ＝Ｈ１＿ｉ＋Ｈ２＿ｉ
【００５１】
Ｓ１１１において、ＭＰＵ１１１０は、ブレンド拡散係数Ｈｉに基づいて、Ｓ１０５で求めた量子化誤差ΔＥを周辺画素ｉ＝１〜４に拡散する。量子化処理が未処理である各周辺画素においては、処理対象画素が変わる度に新たなブレンド拡散係数Ｈｉが求められることになるが、ＭＰＵ１１１０は処理対象画素が変わる度に新たに拡散された誤差ΔＥ×Ｈｉを画素位置（ｘ、ｙ）に対応づけて加算しておく。そして、以降のＳ１０３では、新たに処理対象となった画素（ｘ、ｙ）に対し、加算管理していた誤差ΔＥ×Ｈｉの和を、当該画素の誤差拡散値Ｅ（ｘ、ｙ）として使用する。
【００５２】
Ｓ１１２においてＭＰＵ１１１０は、Ｓ１００で取得した全ての画素について量子化処理が終了したか否かを判定する。まだ処理すべき画素が残っている場合はＳ１０１に戻り、本例の場合は＋ｙ方向に隣接する新たな処理対象画素を設定する。一方、全ての画素について量子化処理が終了したと判定した場合は、本処理を終了する。
【００５３】
なお、以上では、Ｓ１００において１ページに含まれる全画素についての階調データを取得し、その後、取得した全画素の中から処理対象画素を１つずつ設定して行く内容で説明したが、本実施形態はこれに限定されない。Ｓ１００で取得する階調データの単位としては、ｙ方向に並ぶ１画素ラインの単位であってもよいし、複数の画素ラインで構成されるバンド単位であってもよい。
【００５４】
次に、吐出特性情報格納部８（図４参照）に格納する吐出特性情報の取得方法について説明する。本実施形態の記録システムにおいては、吐出特性情報を取得するための吐出特性情報更新モードを適切なタイミングで実行し、吐出特性情報を更新可能とする。
【００５５】
図８は、吐出特性情報更新モードを説明するためのフローチャートである。本処理も画像処理装置のＭＰＵ１１１０が、ＲＯＭ１１２０に記憶されたプログラムに従って実行するものとする。
【００５６】
本処理が開始されると、ＭＰＵ１１１０は、まずＳ４０１において、記録装置１２００にコマンドを送り、所定のテストチャートを記録媒体に記録させる。本実施形態で使用する記録装置は、図２（ａ）および（ｂ）で示したようなフルラインタイプのインクジェット記録装置である。テストチャートには、記録媒体に記録されたドットの特徴（大きさ、位置）が、各色の各ノズルに対応づけて確認できるようなパターンが含まれている。
【００５７】
Ｓ４０２において、ＭＰＵ１１１０は、テストチャートの読み取り処理を実行する。ここでは、１２００ｄｐｉで配列するノズルのそれぞれの吐出量や記録位置ズレ量を把握できるよう、２４００ｄｐｉの解像度で１６ｂｉｔの輝度値を取得するような読取処理を行う。
【００５８】
本実施形態において、テストチャートの読み取り処理を実際に行うデバイスは特に限定されない。例えば、記録装置１２００から出力されたテストチャートをユーザが入力装置１０００にセットし、画像処理装置１１００が、入力装置１０００に対して読み取り処理を行わせるためのコマンドを発信してもよい。また、記録装置１２００にスキャナが備わっている場合は、出力したテストチャートを記録装置１２００のスキャナにそのまま読み取らせてもよい。更に、図１に示すシステムとは別の、オフラインで用意したスキャナなどでテストチャートを読み取る形態としてもよい。
【００５９】
Ｓ４０３において、ＭＰＵ１１１０は、Ｓ４０３の読み取り処理で取得した画像データに基づいて、位置ズレ量Ｘ（ｙ）を導出する。具体的には、個々の画素位置ｙに対応する複数のドットのｙ方向の重心位置を求め、理想位置との差の絶対値を位置ズレ量Ｘ（ｙ）とする。理想位置については、Ｓ４０１でテストチャートを出力する際に用意した画像データから求めることができる。但し、位置ズレ量Ｘ（ｙ）については、必ずしも理想位置に基づいて導出しなくてもよく、例えばｙ方向に隣接する２つのドットやラインの距離に基づいて導出してもよい。いずれにせよ、位置ズレ量Ｘ（ｙ）は、その値が大きいほど白スジ又は黒スジが目立ち易いことを示す値であればよい。
【００６０】
図９（ａ）は、Ｓ４０３で取得される位置ズレ量Ｘ（ｙ）の一例を示す図である。ここでは、１２００ｄｐｉのピッチ（約２１μｍ）で配列する７６８個のノズルのそれぞれについて、理想位置からの位置ズレ量（絶対値）が示されている。
【００６１】
図８のフローチャートに戻る。Ｓ４０４において、ＭＰＵ１１１０は、位置ズレ量Ｘ（ｙ）に対し所定のフィルタ処理および正規化処理を施し、吐出特性情報Ｎ（ｙ）を導出する。フィルタ処理は、ノズル位置毎の位置ズレ量から低周波成分を抽出するための処理である。すなわちフィルタ処理は、図９（ａ）に示す位置ズレ量のノズル位置毎の変化から、ゆるやかな特性を抽出する。フィルタ処理において適切なフィルタサイズは、ノズル数および、生じた位置ズレ量のレンジや変動の急峻さに応じて異なる。例えば複数の異なるフィルタサイズを用いてフィルタ処理を実行し、記録ヘッド固有の位置ズレ量の変化の低周波成分を効果的に抽出することができるフィルタ処理結果を採用すればよい。ここでは、フィルタサイズは２００を用いるものとする。図９（ｂ）は、図９（ａ）の位置ズレ量Ｘ（ｙ）に対し、所定のフィルタ処理を施した結果を示す。ＭＰＵ１１１０は、これらの値を更に０〜５の範囲に規格化し、吐出特性情報Ｎ（ｙ）とする。
【００６２】
Ｓ４０５において、ＭＰＵ１１１０は、Ｓ４０４で導出した吐出特性情報Ｎ（ｙ）を画素位置ｙに対応づけて吐出特性情報格納部８に格納する。以上で本処理は終了する。  
以上説明したＳ４０１〜Ｓ４０５の各工程は、記録ヘッド２０に配列する複数のノズル列２０Ｋ、２０Ｃ、２０Ｍ及び２０Ｙのそれぞれについて、独立且つ並行に行う。すなわち、Ｓ４０１では、ブラック、シアン、マゼンタ、及びイエローのそれぞれについて個々のノズルの吐出状態が確認できるテストチャートを出力する。そして、Ｓ４０２〜Ｓ４０５では、各色のテストチャートを読み取って得られた画像データを基に、位置ズレ量Ｘ（ｙ）及び吐出特性情報Ｎ（ｙ）を色ごとに導出し、吐出特性情報格納部８に格納する。つまり、吐出特性情報格納部８には、ノズル列２０Ｋ、２０Ｃ、２０Ｍ、２０Ｙのそれぞれについて、ノズル数分の吐出特性情報Ｎ（ｙ）が個別に記憶される。
【００６３】
一般に、記録ヘッドには、図９（ａ）に示すような高周波の位置ズレも存在するが、図９（ｂ）に示すような低周波の位置ズレも存在する。そして、記録媒体に記録された画像において視覚的に目立ち易いのは、図９（ｂ）に示すような低周波の位置ズレに起因する濃度ムラである。よって本実施形態では、図９（ｂ）に示すような低周波の位置ズレに起因する濃度ムラを積極的に抑制するための誤差拡散処理を行う。具体的には、分散性に優れた第１拡散係数に加え、位置ズレに伴う濃度ムラへの耐性に優れた第２拡散係数を用意し、第１拡散格納部１１と第２拡散格納部１３にそれぞれ格納しておく（図４参照）。その上で、実際に誤差拡散処理を行う際には、第１拡散係数と第２拡散係数のブレンドによって個々の画素の拡散係数を画素ごとに設定することにより、２つの拡散係数の切り替わりに伴うテクスチャの発生を抑える。
【００６４】
図１０（ａ）および（ｂ）は、第１拡散係数および第２拡散係数の内容と、それぞれの拡散係数を用いて誤差拡散処理を行った場合のドットパターンを比較する図である。ここでは、図６（ａ）で示す処理方向（＋ｙ方向）で量子化処理を進める場合について示している。
【００６５】
図１０（ａ）の左側には、第１拡散係数格納部１１に記憶されている第１拡散係数のマトリクスを示す。図６（ａ）に示す４つの第１拡散係数に対し、ｈ１＿１＝７、ｈ１＿２＝３、ｈ１＿３＝５、ｈ１＿４＝１が設定され、全係数の総和は１６となっている。
【００６６】
このような第１拡散係数を用いて誤差拡散処理を行った場合、処理対象画素（＊）で発生した量子化誤差ΔＥは、周辺に位置する未処理の画素に等方的に拡散される。具体的には、＋ｙ方向に隣接する未処理の画素には、量子化誤差ΔＥの７／１６が拡散される。また、＋ｘ方向に隣接する未処理の画素には、量子化誤差ΔＥの６／１６が拡散され、更にその画素の±ｙ方向に隣接する画素には、量子化誤差ΔＥの３／１６と１／１６がそれぞれ拡散される。
【００６７】
図１０（ａ）の右側には、全ての画素についてＩ（ｘ、ｙ）＝６４（２５％デューティ）が入力され、第１拡散係数を用いて誤差拡散処理を行った場合のドットパターンを示している。量子化誤差ΔＥが等方的に拡散されるため、ドットも一様且つ等方的に分布した状態となり、ブルーノイズ特性を有するドットパターンが得られている。
【００６８】
但し、この場合、ｙ方向に配列するノズルのいずれかのドットが±ｙ方向にずれると、記録媒体においてはドットの粗密が現れ、濃度ムラとして認識されるおそれが生じる。
【００６９】
一方、図１０（ｂ）の左側には、第２拡散格納部１３に記憶されている第２拡散係数のマトリクスを示す。図６（ａ）に示す４つの第２拡散係数に対し、ｈ２＿１＝０、ｈ２＿２＝０、ｈ２＿３＝１６、ｈ２＿４＝０が設定され、全係数の総和は１６である。このような第２拡散係数を用いて誤差拡散処理を行った場合、処理対象画素（＊）で発生した量子化誤差ΔＥは、＋ｘ方向に隣接する画素のみに拡散され、他の画素には拡散されないことになる。
【００７０】
図１０（ｂ）の右側には、全ての画素についてＩ（ｘ、ｙ）＝６４（２５％デューティ）が入力され、第２拡散係数を用いて誤差拡散処理を行った場合のドットパターンを示している。量子化誤差ΔＥが±ｙ方向に拡散させないため、ｙ方向に並ぶ個々の画素は互いに影響を与えず、ｘ座標のほぼ同じ位置で閾値を超える画素が現れる。その結果、複数のドットがｙ方向に連続して形成されるラインが、ｘ方向に周期的に現れるような、グリーンノイズ特性を有するドットパターンが得られている。
【００７１】
この場合、ｙ方向に配列する複数のドットのうちいくつかのドットが、±ｙ方向にずれたとしても、ドットの粗密は現れ難く、濃度ムラが認識されるおそれも少ない。すなわち、第２拡散係数を用いて誤差拡散処理を行うと、ｙ方向の記録位置ズレに対する耐性を向上させることができる。但し、グリーンノイズのパターンはブルーノイズに比べ粒状感が強いため、分散性については第１拡散係数よりも劣る状況となる。
【００７２】
本実施形態では、このような特徴の異なる第１拡散係数と第２拡散係数とを予め用意しておく。そして、実際に誤差拡散処理を行う際には、個々の画素位置の吐出特性情報（すなわち位置ズレの程度）によって、第１拡散係数と第２拡散係数をブレンドするブレンド比率を画素ごとに設定する。具体的には、位置ズレ量が比較的小さな画素位置に対しては濃度ムラの抑制よりも分散性を優先し、第１拡散係数の特徴が第２拡散係数の特徴よりも強く現れるように、ブレンド比率を設定する。一方、位置ズレ量が比較的大きな画素位置に対しては、分散性よりも濃度ムラの抑制を優先し、第２拡散係数の特徴が第１拡散係数の特徴よりも強く現れるように、ブレンド比率を設定する。そして、設定されたブレンド比率に従って導出した拡散係数Ｈｉを用いて量子化誤差ΔＥを周囲の画素に拡散する。
【００７３】
図１１は、図１０（ａ）および（ｂ）に示す第１拡散係数と第２拡散係数において、ブレンド比率αを様々に振った場合の量子化後のドットパターンを比較する図である。ここでも、全ての画素についてＩ（ｘ、ｙ）＝６４（２５％デューティ）が入力された場合を示している。ブレンド比率αの値が０に近いほど第１拡散係数の特徴が大きく現れ、１．０に近いほど第２拡散係数の特徴が大きく現れている。
【００７４】
このような状況において、例えば、α＝０としたドットパターンと、α＝１．０としたドットパターンを隣接して配置すると、これらの違いは違和感を招きテクスチャとして感知されやすい。しかし、ブレンド比率αをｙ方向に連続的に且つ緩やかに変化させれば、ドットパターンは、ブルーノイズ特性を有する状態とグリーンノイズ特性を有する状態との間で連続的に変化し、テクスチャは認識され難くなる。
【００７５】
図１２（ａ）および（ｂ）は、画素位置ｙに応じてブレンド比率を連続的に変化させた場合の、ブレンド比率とドットパターンを示す図である。図１２（ａ）は画素位置ｙとブレンド比率αの関係を示し、図１２（ｂ）は同図（ａ）のブレンド比率のもと、全画素に対しＩ（ｘ、ｙ）＝６４が入力された場合のドットパターンを示している。
【００７６】
図１２（ｂ）によれば、ブレンド比率αが０に近い両端部では第１拡散係数の特徴（ブルーノイズ）が現れ、ブレンド比率が１．０に近い中央部では第２拡散係数の特徴（グリーンノイズ）が現れている。ドットパターンは違和感の無い状態でｙ方向に徐々に変化しており、目立ったテクスチャも出現していない。このように、ブレンド比率αがｙ方向に対して連続的に変化する状態であれば、テクスチャのない滑らかな画像を出力することができる。
【００７７】
本実施形態においては、図９（ａ）に示すような位置ズレ量Ｘ（ｙ）に対し所定のフィルタ処理を施すことにより、図９（ｂ）に示すような滑らかに変化する吐出特性情報Ｎ（ｙ）を求めている。このため、ブレンド比率αもｙ方向に対して連続的に変化し、ドットパターンは違和感の無い状態でｙ方向に徐々に変化する状態となる。
【００７８】
但し、図１２（ｂ）においては、ブレンド比率が１．０に近い中央部で、部分的にではあるが、僅かなテクスチャが確認される。本発明者らの検討によれば、このようなテクスチャの出現は、ブレンド比率αの傾きが急激に変化することに起因すると考えられる。
【００７９】
図１３（ａ）および（ｂ）は、図１２（ｂ）に現れる部分的なテクスチャを抑制する方法を示す図である。図１３（ａ）のブレンド比率αでは、図１２（ａ）のブレンド比率αに対し、ブレンド比率α＝１．０の近傍に比較的振幅の小さい高周波成分の振動を部分的に与え、中央部におけるブレンド比率αの傾きの変化を目立たせないようにしている。その結果、図１３（ｂ）に示すドットパターンにおいても、図１２（ｂ）で確認されるテクスチャは緩和されている。
【００８０】
なお、図１２（ｂ）に示すような部分的なテクスチャを緩和するための手法としては、上記以外にも、様々な手法を採用することができる。例えば、テクスチャが目立つ部分に対してローパスフィルタ処理を行ってブレンド比率αの傾きが急激に変化する箇所を緩やかな傾きに補正してもよいし、テクスチャが目立つ部分に対して高周波のノイズを付加してもよい。
【００８１】
以上説明したように、本実施形態によれば、分散性に優れたドットパターンが得られる第１拡散係数と、記録位置ズレへの耐性が高いドットパターンが得られる第２拡散係数とを用意する。そして、個々の画素位置においては、当該画素位置の位置ズレ量に応じたブレンド比率で、第１拡散係数と第２拡散係数とをブレンドし、得られた拡散係数を用いて誤差拡散処理を行う。これにより、記録ヘッド固有のバンディングを抑制しつつ、拡散係数の切替えに伴うテクスチャのない一様な画像を出力することが可能となる。
【００８２】
（第２の実施形態）
本実施形態も図１、２で示した画像処理システムおよび記録装置を用い、図３で示したフローチャートに従って画像処理を行うものとする。
【００８３】
図１４は、本実施形態の量子化処理部１の機能構成を説明するためのブロック図である。以下、第１の実施形態で説明した図４と異なる点について説明する。
【００８４】
減算器６で生成された量子化誤差ΔＥは、第１拡散処理部２１と第２拡散処理部２２に提供される。第１拡散処理部２１は、量子化誤差ΔＥに対し第１拡散係数格納部１１に記憶されている第１拡散係数ｈ１＿ｉを乗算し、第１拡散値ΔＥ１＿ｉを求める。そして、第１乗算器２３は、第１拡散値ΔＥ１＿ｉに減算器２４より受信した値（１．０−α）を乗算し、第１ブレンド値ΔＥ１＿ｉ´を求める。
【００８５】
一方、第２拡散処理部２２は、量子化誤差ΔＥに第２拡散係数格納部１３に記憶されている第２拡散係数ｈ２＿ｉを乗算し、第２拡散値ΔＥ２＿ｉを求める。そして、第２乗算器２５は、第２拡散値ΔＥ２＿ｉにブレンド比率α決定部より受信したブレンド比率αを乗算し、第２ブレンド値ΔＥ２＿ｉ´を求める。
【００８６】
加算器２６は、第１ブレンド値ΔＥ１＿ｉ´と第２ブレンド値ΔＥ２＿ｉ´とを加算し、拡散値ΔＥｉを求める。
【００８７】
図１５は、図１４で説明した量子化処理部１を用いて、画像処理装置１１００のＭＰＵ１１１０が実行する量子化処理の工程を説明するためのフローチャートである。以下、図１４のブロック図を参照しながら、図１５の各工程について説明する。なお、Ｓ１００〜Ｓ１０７およびＳ１１２の工程は、第１の実施形態で説明した図５のフローチャートと同じ内容であるので、ここでの説明は省略する。
【００８８】
Ｓ２０１において、ＭＰＵ１１１０は、第１拡散係数格納部１１より第１拡散係数ｈ１＿ｉを読み出す。そして、Ｓ１０５で取得した量子化誤差ΔＥに各画素ｉの拡散係数ｈ１＿ｉを乗算し、周辺画素ｉにおける第１拡散値ΔＥ１＿ｉを求める。
ΔＥ１＿ｉ＝ｈ１＿ｉ×ΔＥ
【００８９】
Ｓ２０２において、ＭＰＵ１１１０は、第２拡散係数格納部１３より第２拡散係数ｈ２＿ｉを読み出す。そして、Ｓ１０５で取得した量子化誤差ΔＥに各画素ｉの拡散係数ｈ２＿ｉを乗算し、周辺画素ｉにおける第２拡散値ΔＥ２＿ｉを求める。
ΔＥ２＿ｉ＝ｈ２＿ｉ×ΔＥ
【００９０】
Ｓ２０３において、ＭＰＵ１１１０は、Ｓ２０１で求めた第１拡散値ΔＥ１＿ｉに１．０からブレンド比率αを減算した値（１．０−α）を乗算し、第１ブレンド値ΔＥ１＿ｉ´を求める。
ΔＥ１＿ｉ´＝（１．０−α）×ΔＥ１＿ｉ
【００９１】
Ｓ２０４において、ＭＰＵ１１１０は、Ｓ２０２で求めた第２拡散値ΔＥ２＿ｉにブレンド比率αを乗算し、第２ブレンド値ΔＥ２＿ｉ´を求める。
ΔＥ２＿ｉ´＝α  ×ΔＥ２＿ｉ
【００９２】
Ｓ２０５において、ＭＰＵ１１１０は、Ｓ２０３で求めた第１ブレンド値ΔＥ１＿ｉ´とＳ２０４で求めた第２ブレンド値ΔＥ２＿ｉ´を加算し、周辺画素ｉ＝１〜４のそれぞれについての拡散値ΔＥｉを求める。
【００９３】
量子化処理が未処理である各周辺画素においては、処理対象画素が変わる度に新たな拡散値ΔＥｉが求められることになるが、ＭＰＵ１１１０は、処理対象画素が変わる度に、新たな拡散値ΔＥｉを画素位置（ｘ、ｙ）に対応づけて加算しておく。そして、以降のＳ１０３では、新たに処理対象となった画素（ｘ、ｙ）に対し、加算管理していた拡散値ΔＥｉの和を、当該画素の誤差拡散値Ｅ（ｘ、ｙ）として使用する。
【００９４】
以上説明したように、本実施形態によれば、分散性に優れたドットパターンが得られる第１拡散係数と、記録位置ズレへの耐性が高いドットパターンが得られる第２拡散係数とを用意する。そして、個々の画素位置においては、第１拡散係数に基づいて拡散された第１拡散値と第２拡散係数に基づいて拡散された第２拡散値とを、当該画素位置の位置ズレ量に応じたブレンド比率でブレンドし、得られた拡散値を用いて誤差拡散処理を行う。
【００９５】
このような本実施形態においても、誤差拡散値Ｅ（ｘ、ｙ）を求めるための算出の手順は異なるものの、実質的に第１の実施形態と同様の効果を得ることができる。すなわち、目立ったテクスチャを出現させること無く、第１拡散係数の特徴と第２拡散係数の特徴とが適切にブレンドされたドットパターンを、違和感の無い状態で連続的に配置させることができる。その結果、記録ヘッド固有のバンディングを抑制しつつ一様な画像を出力することが可能となる。
【００９６】
（第３の実施形態）
本実施形態も図１、２で示した画像処理システムおよび記録装置を用い、図３で示したフローチャートに従って画像処理を行うものとする。
【００９７】
図１６は、本実施形態の量子化処理部１の機能構成を説明するためのブロック図である。以下、第１の実施形態で説明した図４と異なる点について説明する。減算器６で生成された量子化誤差ΔＥは、第１乗算器３３と第２乗算器３５に提供される。
【００９８】
第１乗算器３３は、量子化誤差ΔＥに減算器３２より受信した値（１．０−α）を乗算し、得られた第１誤差ΔＥ１を第１拡散処理部３１に提供する。第２乗算器３５は、量子化誤差ΔＥにブレンド比率αを乗算し、得られた第２誤差ΔＥ２を第２拡散処理部３４に提供する。
【００９９】
第１拡散処理部３１は、第１乗算部３３から取得した第１誤差ΔＥ１に第１拡散係数格納部１１に記憶されている第１拡散係数ｈ１＿ｉを乗算し、第１拡散値ΔＥ１＿ｉを求める。一方、第２拡散処理部３４は、第２乗算部３５から取得した第２誤差ΔＥ２に第２拡散係数格納部１３に記憶されている第２拡散係数ｈ２＿ｉを乗算し、第２拡散値ΔＥ２＿ｉを求める。
【０１００】
加算器３６は、第１拡散値ΔＥ１＿ｉと第２拡散値ΔＥ２＿ｉとを加算して、拡散値ΔＥｉを求める。
【０１０１】
図１７は、図１６で説明した量子化処理部１を用いて、画像処理装置１１００のＭＰＵ１１１０が実行する量子化処理の工程を説明するためのフローチャートである。以下、図１６のブロック図を参照しながら、図１７の各工程について説明する。なお、Ｓ１００〜Ｓ１０７およびＳ１１２の工程は、第１の実施形態で説明した図５のフローチャートと同じ内容であるので、ここでの説明は省略する。
【０１０２】
Ｓ３０１において、ＭＰＵ１１１０は、第１拡散係数格納部１１より第１拡散係数ｈ１＿ｉを読み出す。そして、Ｓ１０５で取得した量子化誤差ΔＥに拡散係数ｈ１＿ｉと、（１．０−α）を乗算し、第１拡散値ΔＥ１＿ｉを求める。
【０１０３】
Ｓ３０２において、ＭＰＵ１１１０は、第２拡散係数格納部１３より第２拡散係数ｈ２＿ｉを読み出す。そして、Ｓ１０５で取得した量子化誤差ΔＥに拡散係数ｈ２＿ｉとブレンド比率αを乗算し、第２拡散値ΔＥ２＿ｉを求める。
【０１０４】
Ｓ３０３において、ＭＰＵ１１１０は、Ｓ３０１で求めた第１拡散値ΔＥ１＿ｉとＳ３０２で求めた第２拡散値ΔＥ２＿ｉを加算し、拡散値ΔＥｉを求める。
【０１０５】
量子化処理が未処理である各周辺画素においては、処理対象画素が変わる度に新たな拡散値ΔＥｉが得られることになるが、ＭＰＵ１１１０は、処理対象画素が変わる度に、新たな拡散値ΔＥｉを画素位置（ｘ、ｙ）に対応づけて加算しておく。そして、以降のＳ１０３では、新たに処理対象となった画素（ｘ、ｙ）に対し、管理していた拡散値ΔＥｉの和を、当該画素の誤差拡散値Ｅ（ｘ、ｙ）として使用する。
【０１０６】
以上説明したように、本実施形態によれば、分散性に優れたドットパターンが得られる第１拡散係数と、記録位置ズレへの耐性が高いドットパターンが得られる第２拡散係数とを用意する。そして、個々の画素位置においては、当該画素位置の位置ズレ量に応じたブレンド比率で量子化誤差ΔＥを分配し、分配されたそれぞれの量子化誤差を、第１拡散係数と第２拡散係数によって、周辺画素に拡散する。
【０１０７】
このような本実施形態においても、誤差拡散値Ｅ（ｘ、ｙ）を求めるための算出の手順は異なるものの、実質的に第１、第２の実施形態と同様の効果を得ることができる。すなわち、目立ったテクスチャを出現させること無く、第１拡散係数の特徴と第２拡散係数の特徴とが適切にブレンドされたドットパターンを、違和感の無い状態で連続的に配置させることができる。その結果、記録ヘッド固有のバンディングを抑制しつつ一様な画像を出力することが可能となる。
【０１０８】
（第４の実施形態）
上記実施形態では、図９（ｂ）に示すように、ノズル並び方向に対してブレンド比率αが連続性を有するように、ｙ方向の個々の画素位置に対応づけてブレンド比率αを設定した。これに対し本実施形態では、入力階調値Ｉに対してブレンド比率αが連続性を有するように、０〜２５５の個々の階調値に対応づけてブレンド比率αを設定する。
【０１０９】
なお、本実施形態も図１、２で示した画像処理システムおよび記録装置を用い、図３で示したフローチャートに従って画像処理を行うものとする。
【０１１０】
図１８は、本実施形態の量子化処理部１の機能構成を説明するためのブロック図である。以下、第１の実施形態で説明した図４と異なる点について説明する。
【０１１１】
ブレンド比率格納部４２には、０〜２５５の階調値Ｉのそれぞれに対応づけて、０．０〜１．０のブレンド比率αが予め記憶されている。プレンド比率決定部４１は、ブレンド比率格納部４２にアクセスし、入力部２に入力された階調値Ｉに対応するブレンド比率αを出力する。
【０１１２】
図１９は、図１８で説明した量子化処理部１を用いて、画像処理装置１１００のＭＰＵ１１１０が実行する量子化処理の工程を説明するためのフローチャートである。図１９のフローチャートにおいて、第１の実施形態で説明した図５のフローチャートと異なる点は、Ｓ１０６とＳ１０７の工程が、Ｓ４０１に置き換わっていることである。よって、Ｓ４０１以外の工程についての説明は省略する。
【０１１３】
Ｓ４０１において、ＭＰＵ１１１０は、ブレンド比率格納部４２をアクセスし、処理対象画素（ｘ、ｙ）の階調値Ｉ（ｘ、ｙ）に対応するブレンド比率αを設定する。
【０１１４】
図２０（ａ）および（ｂ）は、ブレンド比率格納部４２に格納された、階調値Ｉに対するブレンド比率αの例を示す図である。図２０（ａ）は、０〜２５５の階調値Ｉのうち、Ｉ＝０近傍のハイライト部とＩ＝２５５近傍の高濃度部ではブレンド比率αが低く、Ｉ＝１２８近傍の中間濃度部でブレンド比率αが高く設定される例を示している。
【０１１５】
一般に、ドットの記録位置ずれが発生した場合、中間濃度部では、記録媒体に対するドット被覆率が比較的大きく変動し、濃度ムラが目立ち易い傾向にある。これに対し、ハイライト部や高濃度部では、記録媒体に対するドットの被覆率は比較的小さく変動し、濃度ムラは目立ち難い傾向にある。
【０１１６】
図２０（ａ）のようにブレンド拡散係数αが階調値Ｉに対応づけて設定されていれば、中間濃度部では第２拡散係数の影響を強めて記録位置ズレへの耐性を高め、ハイライト部や高濃度部では第１拡散係数の影響を強めて分散性を向上させることができる。更に、図２０（ａ）においては、ハイライト部〜中間濃度部〜高濃度部にわたる全階調領域において、ブレンド比率αを０〜１．０の間で滑らかに変化させている。よって、ハイライト部から高濃度部まで階調が様々に変化する実画像を記録する場合においても、これらドットパターンを違和感の無い状態で連続的に変化させることができる。
【０１１７】
一方、図２０（ｂ）は、ブレンド比率αが階調値Ｉによらず一定の値α＝０．３に設定された場合を示している。記録位置ズレが発生し難い場合や、記録位置ズレに伴う濃度ムラが階調値Ｉによらず安定しているような場合には、図２０（ａ）のような対応づけが好ましいと言える。
【０１１８】
図２１は、上記のような階調値Ｉとブレンド比率αの対応づけを求める手順を説明するためのフローチャートである。本処理についても、画像処理装置のＭＰＵ１１１０が、ＲＯＭ１１２０に記憶されたプログラムに従って実行するものとする。
【０１１９】
本処理が開始されると、ＭＰＵ１１１０は、まずＳ１６０１において、記録装置１２００にコマンドを送り、所定のテストチャートを記録媒体に記録させる。
【０１２０】
図２２は、Ｓ１６０１で出力されるテストチャートの一例を示す図である。ここでは、０〜２５５の階調値Ｉのうち、Ｉ＝１７、３４・・・２３８の１４段階の階調値Ｉのそれぞれについて、ブレンド比率αを０．１刻みで１０段階に振った場合のドットパターンを、個々のパッチとして記録している。ドットパターンの生成方法は、第１の実施形態で説明した方法を採用する。
【０１２１】
図２２のようなテストチャートが出力されると、Ｓ１６０２において、ＭＰＵ１１１０は、テストチャートの読み取り処理を行う。ここでは、１２００ｄｐｉの解像度で１６ｂｉｔの輝度値を取得するような読取処理を行う。なお、本実施形態においても、テストチャートの読み取りを行わせるデバイスは特に限定されない。読み取り処理は入力装置１０００に行わせてもよいし、記録装置１２００に備えられたスキャナに行わせてもよいし、他のデバイスに読み取らせてもよい。
【０１２２】
Ｓ１６０３において、ＭＰＵ１１１０は、Ｓ１６０２の読み取り処理で取得した画像データに基づいて、個々のパッチを評価する。ここでは、粒状性、バンディングおよびテクスチャの３項目についての評価を行う。粒状性については、公知のＲＭＳ粒状性評価法等を用いて評価することができる。バンディングについては、個々のパッチにおけるｙ方向の濃度変動値を累積して評価してもよいし、スキャン画像を周波数変換しｙ方向の低周波成分のパワーを累積して評価してもよい。テクスチャについては、スキャン画像に含まれる低周波なモアレ成分を抽出することによって評価することができる。
【０１２３】
どのような方法を採用するにせよ、３項目のそれぞれの評価値について、好ましい状態であるほど高い評価値が得られることが求められる。その上で、３項目それぞれの評価値を規格化し、総合的に好ましい状態の場合に高い総合評価値が得られるように、３項目の評価値の加重平均処理が行われればよい。
【０１２４】
Ｓ１６０４において、ＭＰＵ１１１０は、Ｓ１６０３で求めた総合評価値に基づいて、０〜２５５の階調値Ｉのそれぞれに対応づけてブレンド比率αを決定する。決定方法は特に限定されるものではないが、例えば、１４段階の階調値Ｉのそれぞれについて、総合評価値の最も高いパッチのブレンド比率αを階調値Ｉに対応付けて仮設定する。更に、仮設定されたブレンド比率αに対し階調値Ｉが増減する方向のローパスフィルタ処理を施すことにより、０〜２５５の階調値Ｉと０〜１．０のブレンド比率αとが、より連続的で滑らかに変化するような曲線を導出する。そして、この曲線に基づいて、０〜２５５の階調値のそれぞれに、０〜１．０のブレンド比率αを対応づけた情報を、図１８に示すブレンド比率格納部４２に格納すればよい。
【０１２５】
また、例えば図２０（ｂ）のようにブレンド比率αを一律にすることが決められている場合には、１４段階の階調値Ｉのそれぞれで総合評価値の最も高いブレンド比率αを選び、最も多く選出されたブレンド比率αを全階調値に対応づけて一律に設定してもよい。以上で本処理を終了する。
【０１２６】
以上説明した本実施形態によれば、分散性に優れたドットパターンが得られる第１拡散係数と、記録位置ズレへの耐性が高いドットパターンが得られる第２拡散係数とを用意する。そして、個々の画素位置においては、当該画素位置の階調値に応じたブレンド比率で、第１拡散係数と第２拡散係数とをブレンドし、得られた拡散係数を用いて誤差拡散処理を行う。これにより、各階調においてバンディングを抑制しつつテクスチャのない一様な画像を出力することが可能となる。
【０１２７】
なお、本実施形態では、ブレンド比率αを入力値に対し連続的に変化させる方法を第１の実施形態の変形例として説明したが、本実施形態の方法は、第２、第３の実施形態の変形例にも応用することもできる。すなわち、第２、第３の実施形態においても、吐出特性格納部８をブレンド比率情報格納部４２に変更し、ブレンド比率決定部４１が階調値Ｉに対応するブレンド比率αを決定するようにすれば、本実施形態の効果を得ることができる。
【０１２８】
（第５の実施形態）
本実施形態では、第１拡散係数と第２拡散係数とをブレンドする際に発生する演算誤差を補正するための処理を、第１、第４の実施形態に対して更に追加した誤差拡散処理について説明する。
【０１２９】
図２３は、上記演算誤差の具体例を説明するための模式図である。ここでは、第１または第４の実施形態において、第１拡散係数１９１および第２拡散係数１９２が０〜１６の５ビットで管理され、ブレンド比率α及び減算値（１−α）が、０〜１．０の８ビットで管理される場合を示している。この場合、第１および第２の拡散係数は１／１６＝０．０６２５の分解能を有し、ブレンド比率α及び減算値（１−α）は１／２５５＝０．０３９０６２５の分解能を有することになる。このような制限の下、例えばブレンド比率をα＝０．９とし、第１の実施形態で説明した方法で第１拡散係数１９１と第２拡散係数１９２をブレンドすると、図の右側に示すようなブレンド拡散係数１９３が得られる。
【０１３０】
図において、第１拡散係数１９１と第２拡散係数１９２における、係数ｈ１＿ｉ及びｈ２＿ｉそれぞれの総和は「１６」である。これに対し、ブレンド後のブレンド拡散係数１９３における各係数Ｈｉの総和は「１４」であり、拡散係数の総和が減少してしまっている。この場合、ブレンド拡散係数１９３をそのまま用いて誤差拡散処理を行うと、処理対象画素で発生した量子化誤差ΔＥのうち、１４／１６に相当する誤差しか周囲画素に分散されないことになり、実画像の濃度が入力画像の濃度よりも低くなってしまう。
【０１３１】
図２３では、ブレンド後係数Ｈｉの総和が分散係数の総和よりも小さい場合を示したが、ブレンド後係数Ｈｉの総和は分散係数の総和よりも大きくなってしまう場合もある。この場合は、処理対象画素で発生した量子化誤差ΔＥよりも大きな誤差が周辺画素に分散されことになり、実画像の濃度は入力画像の濃度よりも高くなってしまう。すなわち、拡散係数を保存するためのビット数（拡散係数の分解能）がブレンド比率αを保存するためのビット数（ブレンド比率の分解能）よりも小さい場合、ブレンドする際に演算誤差が発生し、拡散係数のＤＣ（寸法基準）が損なわれてしまう場合がある。
【０１３２】
このため、本実施形態では、拡散係数の総和がブレンド後にも維持されるように、ブレンド拡散係数Ｈｉを補正するための処理を行う。
【０１３３】
図２４は、上記補正処理を説明するためのフローチャートである。本補正処理は、図５または図１９で説明したフローチャートのＳ１１０とＳ１１１の間で実行すればよい。以下、各工程について具体的に説明する。
【０１３４】
Ｓ５０１において、ＭＰＵ１１１０は、Ｓ１１０で求められた４つのブレンド拡散係数Ｈｉの総和ΣＨｉを求める。また、Ｓ５０２において、ＭＰＵ１１１０は、既知数である第１、第２拡散係数の総和Ｓと、Ｓ５０１で求めたブレンド拡散係数Ｈｉの総和ΣＨｉの差分ΔＨを求める。
【０１３５】
図２３の場合、ΔＨ＝Ｓ−ΣＨｉ＝１６−１４＝２となる。
【０１３６】
Ｓ５０３において、ＭＰＵ１１１０は、Ｓ５０２で求めた差分ΔＨを周辺画素ｉに再分配する。ここで、ΔＨ＞０であれば、周辺画素ｉのブレンド拡散係数Ｈｉは、より大きな値に補正されることになる。また、ΔＨ＜０であれば、周辺画素ｉのブレンド拡散係数Ｈｉは、より小さな値に補正されることになる。
【０１３７】
以上のような補正処理を行えば、第１、第４実施形態のブレンド処理において演算誤差が発生しても、拡散係数のＤＣを維持し、入力画像の濃度を実画像でも保存することが可能となる。
【０１３８】
なお、以上では第１、第４実施形態に採用可能な補正処理について説明したが、演算誤差に伴うＤＣの損失は、他の実施形態でも起こりうる現象である。例えば、階調値Ｉのビット数がブレンド比率αのビット数よりも小さい場合には、第２の実施形態や第３の実施形態でも、処理対象画素の量子化誤差ΔＥの全てが、周辺画素ｉに拡散誤差ΔＥｉとして拡散されない場合がある。いずれにせよ、演算誤差に伴うＤＣの損失が懸念される算出工程が存在する場合は、その工程の後に、ＤＣの損失が懸念されるパラメータについてＤＣの損失量を求め、再分配すれば入力画像の濃度は保存される。
【０１３９】
（その他の実施形態）
以上では、図１に示す画像処理システムを用い、記録装置１２００で記録するための画像処理を画像処理装置１１００が実行する内容で説明したが、本発明はこのような形態に限定されるものではない。図３に示す画像処理の一部又は全ての工程を記録装置１２００のＭＰＵ１２０１が行う形態としてもよい。記録装置１２００が上記実施形態で説明した誤差拡散処理を行う場合、図５、図１５、図１７、図１９および図２４の全ての工程は、記録装置１２００のＭＰＵ１２０１が実行することとなる。
【０１４０】
また、以上の実施形態では図２を用い、記録装置１２００がフルライン型の記録装置である場合を例に説明したが、シリアル型の記録装置にも適用可能である。シリアル型の記録装置であっても、記録ヘッド特有のバンディングは存在し、シングルパス記録は勿論、マルチパス記録であってもマルチパス数が少ない場合に、同様の効果を発揮させることができる。
【０１４１】
また、上述の実施形態は、例えば電子写真方式や感熱方式の記録装置にも適用できる。電子写真方式や感熱方式の記録装置においても記録素子の記録特性のばらつきなど、装置の様々な要因で装置固有の濃度ムラが出現することはあり、上記実施形態で説明した同様の効果を発揮させることができる。
【０１４２】
また、以上の実施形態では８ビット２５６値の階調データを１ビット２値の記録データに量子化する場合について説明してきたが、他の例を採用することもできる。例えば階調データの階調数は更に多くても少なくても良いし、量子化値（Ｎ値）は３値以上であってもよい。階調値Ｉを、Ｎ値（Ｎは３以上の整数）に量子化する場合、比較器４（図４参照）は比較値Ｉ´（ｘ、ｙ）と予め用意された（Ｎ−１）個の閾値Ｔとの大小関係から、Ｎ値の量子化値を求めればよい。この場合であっても、Ｎ段階の量子化代表値Ｒ（ｘ、ｙ）と比較値Ｉ´（ｘ、ｙ）との量子化誤差ΔＥは発生するため、これを上記実施形態と同様の方法で周辺画素に拡散すれば、上記実施形態と同様の効果をえることができる。
【０１４３】
また、第１〜第３の実施形態では、吐出特性情報格納部８において、記録位置ズレの程度を示す吐出特性情報Ｎ（ｙ）をノズル位置（画素位置ｙ）に対応づけて保存したが、記録位置ズレの情報は、ｘ座標の位置やｘｙ座標の位置に対応づけて保存してもよい。例えば、記録装置においては、搬送ローラの駆動に起因する濃度ムラが、記録媒体の搬送方向に周期的に現れる場合がある。また、シリアル型の記録装置の場合は、キャリッジの速度変動に伴う濃度ムラがキャリッジの進行方向に周期的に現れる場合がある。すなわち、±ｙ方向の記録位置ズレだけでなく±ｘ方向の記録位置ズレが濃度ムラを引き起こす場合も有る。
【０１４４】
±ｘ方向の記録位置ズレに伴う濃度ムラが発生する場合は、第２拡散係数として、処理対象画素（＊）で発生した量子化誤差ΔＥを、ｙ方向に隣接する画素に積極的に拡散するような拡散係数を用意することが好ましい。この場合、複数のドットがｘ方向に連続して形成されるラインがｙ方向に周期的に現れるような、すなわち図１０（ｂ）のドットパターンと直交関係にあるような、グリーンノイズ特性を有するドットパターンが得られることになる。そして、濃度ムラ情報Ｎ（ｘ、ｙ）を記録媒体のｘｙ座標に対応づけて保存しておき、それぞれの位置で最適なブレンド比率αで第１拡散誤差と第２拡散誤差のブレンド処理を行えば、ｙ方向の濃度ムラもｙ方向の濃度ムラも適切に抑制することができる。無論この場合であっても、ブレンド比率αをｘｙ方向に連続的に且つ緩やかに変化させることにより、テクスチャの発生を抑えることができる。
【０１４５】
また、以上の実施形態では、第１拡散係数と第２拡散係数の２つを用意し、これらを一定のブレンド比率αでブレンドする形態を例に説明してきたが、ブレンドする拡散係数は３種類以上用意してもよい。この場合、例えば、第１拡散係数のブレンド比率をα、第２拡散係数のブレンド比率をβ、第３拡散係数のブレンド比率をγとしたとき、α、β及びγの和が１．０になるように、ブレンド比率α、β、γが設定されればよい。
【０１４６】
更に、第１〜３の実施形態では、ブレンド比率決定部７が吐出特性情報格納部８より読み出した吐出特性情報Ｎ（ｙ）に基づいてブレンド比率αを決定したが、吐出特性情報格納部８には画素位置ｙに対応するブレンド比率α（ｙ）が直接格納されていてもよい。
【０１４７】
また、図８に示すフローチャートでは、テストチャートを読み取って取得した画像データに基づいて、位置ズレ量Ｘ（ｙ）を求め、更にフィルタ処理を施すことによって吐出特性情報Ｎ（ｙ）を導出したが、吐出特性情報Ｎ（ｙ）の導出方法はこれに限らない。例えば、テストチャートに一低階調のパッチを含め、このパッチの読取画像のｙ方向における濃度変位に基づいて吐出特性情報Ｎ（ｙ）を導出してもよい。
【０１４８】
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１４９】
１    量子化処理部
７    ブレンド比率決定部
１１    第１拡散係数格納部
１３    第２拡散係数格納部
１５    拡散処理部
１１００    画像処理装置
【図１】
2020124835.tif 000003
【図２】
2020124835.tif 000004
【図３】
2020124835.tif 000005
【図４】
2020124835.tif 000006
【図５】
2020124835.tif 000007
【図６】
2020124835.tif 000008
【図７】
2020124835.tif 000009
【図８】
2020124835.tif 000010
【図９】
2020124835.tif 000011
【図１０】
2020124835.tif 000012
【図１１】
2020124835.tif 000013
【図１２】
2020124835.tif 000014
【図１３】
2020124835.tif 000015
【図１４】
2020124835.tif 000016
【図１５】
2020124835.tif 000017
【図１６】
2020124835.tif 000018
【図１７】
2020124835.tif 000019
【図１８】
2020124835.tif 000020
【図１９】
2020124835.tif 000021
【図２０】
2020124835.tif 000022
【図２１】
2020124835.tif 000023
【図２２】
2020124835.tif 000024
【図２３】
2020124835.tif 000025
【図２４】
2020124835.tif 000026