"記録媒体の搬送方向に直交する第1方向に所定の配列周期で配置されるノズルを有する印刷手段を備え、その印刷手段で印刷する画像データを誤差拡散処理によって生成する画像処理装置であって、前記ノズルの配列周期に相当する周波数を±2n倍nは0以上の整数した周波数である対象周波数のうち、少なくとも1つ以上の対象周波数についての周波数スペクトルを、又は、前記対象周波数の近傍の周波数である近傍周波数のうち少なくとも1つ以上の近傍周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数を用いて、前記誤差拡散処理で使用するしきい値を変動する変動手段を備えていることを特徴とする画像処理装置。前記対象周波数、又は、前記近傍周波数は、±10cycle/mmの範囲に含まれ、前記乱数は、ホワイトノイズから生成された乱数であることを特徴とする請求項1に記載の画像処理装置。前記近傍周波数は、前記対象周波数の±0.2cycle/mmの範囲内であることを特徴とする請求項1又は2にに記載の画像処理装置。前記乱数は、前記第1方向と、その第1方向に直交する第2方向との2次元状のホワイトノイズから、前記第2方向の周波数が±10cycle/mmの範囲に含まれる周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数であることを特徴とする請求項1から3のいずれかに記載の画像処理装置。前記乱数を、テーブル化して記憶する記憶手段を備え、前記変動手段は、前記記憶手段にテーブル化して記憶された乱数を用いて変動することを特徴とする請求項1から4のいずれかに記載の画像処理装置。画像処理プログラムであって、記録媒体の搬送方向に直交する第1方向に所定の配列周期で配置されるノズルを有する印刷手段で印刷する画像データを誤差拡散処理によって生成するコンピュータを、前記ノズルの配列周期に相当する周波数を±2n倍nは0以上の整数した周波数である対象周波数のうち、少なくとも1以上の対象周波数についての周波数スペクトルを、又は、前記対象周波数の近傍の周波数である近傍周波数のうち少なくとも1以上の近傍周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数を用いて、前記誤差拡散処理で使用するしきい値を変動する変動手段として機能させることを特徴とする画像処理プログラム。記録媒体の搬送方向に直交する第1方向に所定の配列周期で配置されるノズルを有する印刷手段で印刷する画像データを誤差拡散処理によって生成する場合に用いられる乱数の生成方法であって、ホワイトノイズを生成する生成ステップと、その生成ステップで生成されたホワイトノイズをフーリエ変換するフーリエ変換ステップと、そのフーリエ変換ステップでフーリエ変換した出力から、前記ノズルの配列周期に相当する周波数を±2n倍nは0以上の整数した周波数である対象周波数のうち、少なくとも1以上の対象周波数についての周波数スペクトルを、又は、前記対象周波数の近傍の周波数である近傍周波数のうち少なくとも1以上の近傍周波数についての周波数スペクトルを、0にする、又は、削減するフィルタステップと、そのフィルタステップでフィルタした後の出力を逆フーリエ変換する逆フーリエ変換ステップとを備えていることを特徴とする乱数の生成方法。しかしながら、誤差拡散処理における省メモリ化、高速化の要請から、乱数をテーブル化し、サイクリックにブルーノイズを擬似発生させると、ブルーノイズでは、テーブルの繰り返しにより、高周波成分よりも低周波成分が高くなり、結局、テーブルごとのむらがバンディングとして発生してしまうという問題点があった。本発明は、上記の問題点を解決するためになされたものであり、バンディングを抑制することを目的とする。この目的を達成するために、本発明の装置は、記録媒体の搬送方向に直交する第1方向に所定の配列周期で配置されるノズルを有する印刷手段を備え、その印刷手段で印刷する画像データを誤差拡散処理によって生成するものであって、前記ノズルの配列周期に相当する周波数を±2n倍nは0以上の整数した周波数である対象周波数のうち、少なくとも1つ以上の対象周波数についての周波数スペクトルを、又は、前記対象周波数の近傍の周波数である近傍周波数のうち少なくとも1つ以上の近傍周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数を用いて、前記誤差拡散処理で使用するしきい値を変動する変動手段を備えている。よって、バンディングを抑制することができる。上記の装置において、前記対象周波数、又は、前記近傍周波数は、±10cycle/mmの範囲に含まれ、前記乱数は、ホワイトノイズから生成された乱数である。よって、ホワイトノイズ特性を極力維持した乱数を用いて、誤差拡散処理で使用するしきい値を変動することができる。即ち、周波数スペクトルを0、又は、削減する周波数対象周波数、近傍周波数を、視覚上無視できる±10cycle/mmの範囲に限定することで、ホワイトノイズの高周波成分が低下するのを抑制できる。よって、バンディングを抑制することができる。上記の装置において、前記近傍周波数は、前記対象周波数の±0.2cycle/mmの範囲内である。よって、対象周波数と一致する周波数スペクトルを、演算上、0又は、削減できなくても、対象周波数に対し、かかる範囲の近傍周波数の周波数スペクトルを操作することで、バンディングを抑制できる。上記の装置において、前記乱数は、前記第1方向と、その第1方向に直交する第2方向との2次元状のホワイトノイズから、前記第2方向の周波数が±10cycle/mmの範囲に含まれる周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数である。よって、2次元状のホワイトノイズ特性を極力維持した乱数を用いて、誤差拡散処理で使用するしきい値を変動することができる。即ち、第2方向についても、第1方向と同様に、0、又は、削減する周波数スペクトルを、視覚上無視できる±10cycle/mmの範囲に限定することで、2次元状のホワイトノイズの高周波成分が低下するのを抑制できる。よって、2次元の乱数テーブルを用いても、バンディングを抑制できる。上記の装置において、前記乱数を、テーブル化して記憶する記憶手段を備え、前記変動手段は、前記記憶手段にテーブル化して記憶された乱数を用いて変動する。よって、高速に、且つ、安価に、バンディングが抑制された画像データを生成できる。なお、本発明は、画像処理装置を制御する制御装置、画像処理方法、画像処理装置を制御する画像処理プログラム、該プログラムを記録する記録媒体等の種々の態様で実現可能である。以下、本発明の実施形態について、添付図面を参照して説明する。図1は、本発明の一実施形態であるプリンタ1の断面図である。プリンタ1は、視覚的にバンディングが抑制された高品質な画質で画像を印刷できるよう構成されている。プリンタ1はラインプリンタであって、制御部10、印刷ヘッド2、用紙Pを搬送する搬送機構21、搬送機構21に用紙Pを給紙する給紙機構30、印刷ヘッド2により画像が形成された用紙Pを排紙する排紙部90を主に備えている。搬送機構21は、ベルトローラ6,7と、ベルトローラ6,7に巻き掛けられた搬送ベルト8とを有する。ベルトローラ7には、搬送モータ22図3が接続され、搬送モータ22からの回転力によって図1中矢印R方向に回転する。このとき、搬送ベルト8が用紙Pを搬送方向Tに搬送するように走行し、従動ローラであるベルトローラ6も回転する。搬送ベルト8を挟んでベルトローラ6と対向する位置には、ニップローラ4が配置される。ニップローラ4は、給紙機構30により送られてきた用紙Pを、搬送ベルト8の外周面8aに押さえ付けて、搬送ベルト8の外周面8aに保持させる。一方、搬送ベルト8の内周側には、プラテン19が配置され、内周側から搬送ベルト8を支持する。印刷ヘッド2は、4つのインク色C,M,Y,Kシアン、マゼンタ、イエロー、ブラックごとに設けられており、枠体のフレーム3に支持されて、用紙Pの搬送方向Tに併設されている。各印刷ヘッド2は、それぞれチューブを介してインク貯留部60と接続される。インク貯留部60には、インクを色別に貯留するインクカートリッジ61が設けられ、各印刷ヘッド2には、それぞれ対応する色のインクが供給される。搬送ベルト8に保持された用紙Pが、4つの各印刷ヘッド2を順番に通過することで用紙Pに画像が形成される。各印刷ヘッド2は、インク吐出面2aのノズル2d図2からインク滴を吐出し、ドットとして用紙Pに着弾させることにより、ドット列を用紙Pに形成する。プリンタ1は、ドット列の形成を繰返し、用紙Pに多数のドット列を形成することにより、画像を印刷する。画像が印刷された用紙Pは、搬送機構21によって更に下流側の排紙部90へ搬送される。排出機構90は、搬送機構21から搬送されてきた用紙Pを上方に搬送し、排紙する。図2aは、印刷ヘッド2の配置を示した図である。1色の印刷ヘッド2は、2つのヘッドユニット2bを含む。一方のヘッドユニット2bは、搬送方向Tに平行な方向における第1の位置L1に設けられ、他方のヘッドユニット2bは、搬送方向Tに平行な方向における第2の位置L2に設けられる。各ヘッドユニット2bのインク吐出面2aには、台形状の部分ヘッド2cが複数設けられる。複数の部分ヘッド2cは、搬送方向Tに直交する方向Sに沿って並ぶ。各部分ヘッド2cには、図示しない多数のノズルが形成される。尚、以下の説明では、印刷ヘッド2の長手方向をX方向、X方向に直交する方向をY方向とする。理想的には、図2aに示すように、印刷ヘッド2のX方向と方向Sとが平行となるように、印刷ヘッド2がフレーム3図1に取り付けられる。ただし、実際には、製造段階における取り付けの誤差や経時変動を考慮し、プリンタ1の実用上許容される範囲であれば、印刷ヘッド2のX方向が方向Sに対して傾いていても良い。図2bは、ノズル配列周期を説明するための図である。X方向において互いに隣り合うノズルは、Y方向において互いに位置が異なるような配置されている。X方向のノズル間隔は、600dpi、即ち、1/600×25.4×1000=42.333μmの間隔で配置されている。また、X方向のノズルは、16ノズル周期で配置されている。更に、各ノズルは、サイズが互いに異なる3種類のドット大ドット、中ドット、小ドットを、用紙P上に形成可能に構成されている。図3は、プリンタ1の主な電気的構成を模式的に示すブロック図である。図3に示すように、プリンタ1は、制御部10、インターフェイス16I/F16、印刷ヘッド2、搬送モータ22を有する。制御部10は、CPU11、ROM12、RAM13、フラッシュメモリ14、及びASIC15を備え、これらは、バスラインを介して互いに接続される。また、制御部10、I/F16、印刷ヘッド2、搬送モータ22は、ASIC15を介して互いに接続される。CPU11は、ROM12やフラッシュメモリ14に記憶される固定値やプログラムに従って、プリンタ1が有している各機能の制御や、ASIC15と接続された各部を制御する。ROM12は、制御プログラム12a、誤差拡散パラメータメモリ12b、カラープロファイルメモリ12cなどを格納した書換不能なメモリである。制御プログラム12aは、図9に示す誤差拡散処理を含む印刷処理をCPU11に実行させるためのプログラムである。誤差拡散パラメータメモリ12bは、3種類のドット毎に誤差拡散パラメータを格納する。カラープロファイルメモリ12cは、R赤,G緑,B青の各色の画素値を、C,M,Y,Kの各色の画素値に色変換するための予め定められた対応関係カラープロファイルを格納する。カラープロファイルメモリ12cにも、3種類のカラープロファイルが記憶されている。RAM13は、書換可能な揮発性のメモリであり、画像データメモリ13a、誤差拡散パラメータ設定メモリ13b、カラープロファイル設定メモリ13c、乱数カウンタ13dが設けられる。画像データメモリ13aには、処理対象の画像データが格納される。誤差拡散パラメータ設定メモリ13bには、誤差拡散パラメータメモリ12bに記憶された3種類の誤差拡散パラメータのうち、選択された1種類の誤差拡散パラメータが設定される。カラープロファイル設定メモリ13cには、カラープロファイルメモリ12cに記憶された3種類のカラープロファイルのうち、誤差拡散パラメータ設定メモリ13bに設定された誤差拡散パラメータと組み合わせて使用されるべき1種類のカラープロファイルが設定される。乱数カウンタ13dは、乱数の読み出し位置を示すカウンタである。乱数カウンタ13dの値は、1画素分の処理を終了する度に1ずつカウントアップされ、その値が、後述する乱数テーブル14aに格納された乱数の個数に達すると、1に設定される。制御部10は、乱数テーブル14aから読み出した乱数を用いて誤差拡散処理を実行する。フラッシュメモリ14は書換可能な不揮発性のメモリであり、乱数テーブル14aが格納されている。乱数テーブル14aは、図4bに示す乱数テーブル作成処理によって生成され、例えば、512個の乱数で構成されている。制御部10は、画像データに対して、誤差拡散処理を含む画像処理を施して印刷データを生成し、ASIC15に供給する。ASIC15は、印刷データに基づいて、印刷ヘッド2および搬送モータ22などを制御することにより、用紙Pに画像を印刷する。図4aは、印刷ヘッド2の周波数スペクトルと、VTF曲線とを示す図である。図4aに示す通り、印刷ヘッド2の周波数スペクトルは、ア1.476cycle/mmを中心とした±0.2cycle/mmの範囲、イ2.952cycle/mm中心とした±0.2cycle/mmの範囲、ウ5.904cycle/mm中心とした±0.2cycle/mmの範囲、エ11.808cycle/mm中心とした±0.2cycle/mmの範囲、オ23.616cycle/mm中心とした±0.2cycle/mmの範囲において、他の周波数スペクトルと比べ高い周波数スペクトルを示していることが分かる。即ち、印刷ヘッド2では、かかる周波数で、縦筋状のむら所謂、バンディングが発生していることが分かる。印刷ヘッド2に形成されているノズルは、上述した通り、X方向に600dpi、即ち、1/600×25.4×1000=42.333μmの間隔で配置されている。また、X方向のノズルは、16ノズル周期で配置されている。よって、X方向のノズルの配列周期は、600/16=37.5dpiであり、このノズルの配列周期を周波数に換算すると、1/1/600/16×25.4=1.476cycle/mmとなる。また、X方向のノズルの配列周期に相当する周波数である1.476cycle/mmの2倍は2.952cycle/mm、4倍は5.904cycle/mm、8倍は11.808cycle/mm、16倍は23.616cycle/mmとなる。即ち、バンディングは、X方向のノズルの配列周期に相当する周波数である1.476cycle/mmを、±2n倍nは0以上の整数した周波数以下、対象周波数と、その対象周波数を中心とした±0.2cycle/mmの範囲とで発生していることが分かる。また、図4bに示す観察距離を30cmとするVTF曲線観察距離を30cmに示す通り、10cycle/mm以上の周波数は、周波数スペクトルがほぼ0となり、視覚的に無視できる範囲である。本実施形態は、対象周波数X方向のノズルの配列周期に相当する周波数である1.476cycle/mmを、±2n倍nは0以上の整数した周波数であって、±10cycle/mmの範囲内の周波数について、その周波数スペクトルを0にした周波数特性を持つ乱数を、誤差拡散処理の閾値に加えるものである。図4bは、乱数テーブル作成処理を示すフローチャートである。乱数テーブル14aに記憶される乱数テーブルは、この処理によって作成される乱数テーブルが記憶される。即ち、乱数テーブル14aに記憶される乱数テーブルは、プリンタ1の工場出荷前に、図示しない乱数テーブル生成装置以下、単に、生成装置と称すによって作成される。尚、乱数テーブル生成装置としては、例えば、パーソナルコンピュータが該当する。生成装置は、公知のアルゴリズムであるメルセンヌツイスターを使って、一様分布乱数であるホワイトノイズを生成しS401、その生成したホワイトノイズを離散フーリエ変換するS402。尚、一様分布乱数であるホワイトノイズを生成可能なアルゴリズムであれば、例えば、線形合同法など、メルセンヌツイスターに限定されない。生成装置は、離散フーリエ変換された信号値から、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数対象周波数であって、±10cycle/mmの範囲の周波数について、その周波数スペクトルを0にするS403。生成装置は、S403の処理から生成される信号値を、逆離散フーリエ変換しS404、虚数値の絶対値をとり、実数値に変換するS405。また、S405の処理を経て出力される信号値の最小値が0、最大値が1となる範囲で線形変換しS406、更に、S406の処理を経て出力される信号値について、整数値変換処理を実行しS407、本処理を終了する。ここで、S407の整数値変換処理について、図5を参照して説明する。図5は、整数値変換処理を示すフローチャートである。この処理は、実数を-nからn-1の整数値に変換する処理である。生成装置は、信号値を取得しS502、その信号値が0以上1/2n未満かを判断しS502、0以上1/2n未満であればS502:Yes、取得した信号値に-nを設定しS503、S513の処理に移行する。生成装置は、取得した信号値が0以上1/2n未満でなければS502:No、取得した信号値が1/2n以上2/2n未満かを判断しS504、1/2n以上2/2n未満であればS504:Yes、取得した信号値に-n1を設定しS505、S513の処理に移行する。生成装置は、取得した信号値が、1/2n以上2/2n未満でなければS504:No、取得した信号値が2/2n以上3/2n未満かを判断しS506、2/2n以上3/2n未満であればS506:Yes、取得した信号値に-n2を設定しS507、S513の処理に移行する。その後は、上述した処理に相当する処理を繰り返し、生成装置は、取得した信号値が、1-3/2n以上1-2/2n未満かを判断しS508、1-3/2n以上1-2/2n未満であればS508:Yes、取得した信号値にn-3を設定しS509、S513の処理に移行する。生成装置は、取得した信号値が、1-3/2n以上1-2/2n未満でなければS508:No、取得した信号値が1-2/2n以上1-1/2n未満かを判断しS510、1-2/2n以上1-1/2n未満であればS510:Yes、取得した信号値にn-2を設定しS511、S513の処理に移行する。生成装置は、1-2/2n以上1-1/2n未満でなけばS510:No、取得した信号値にn-1を設定しS512、S513の処理に移行する。生成装置は、S513の処理において、全データについて、上述した処理を実行したかを判断しS513、全データについて、上述した処理を実行していなければS513:No、S501からの処理を繰り返し、全データについて、上述した処理を実行していればS513:Yes、本処理を終了する。こうして、整数値変換処理によって、図4bに示すS406の処理を経て出力される信号値の全データてについて、実数値を-nからn-1の整数値に変換することができる。図6乃至図8は、図4で説明した乱数テーブル作成方法において生成される各信号を説明するための図である。図6aは、一様分布乱数信号を示す図であり、図4のS401で生成される信号を示している。図6bは、図6aに示す一様分布乱数信号を離散フーリエ変換した信号を示す図であり、図4のS402で生成される信号を示している。図6bから、図6aに示す一様分布乱数信号を離散フーリエ変換した信号は、ホワイトノイズ特性を有していることが分かる。図6cは、図6aの一様分布乱数信号のヒストグラムを示し、乱数値が一様に分布していることが分かる。図6d1は、図6bに示す離散フーリエ変換した信号に対して、ハイパスフィルタをかけた信号を示す図であり、この信号は低周波成分が減少したブルーノイズ特性を有していることが分かる。図6d2は、図6bに示す離散フーリエ変換した信号に対して、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数対象周波数であって、且つ、±10cycle/mmの範囲の周波数について、その周波数スペクトルを0にした信号を示す図である。即ち、図4のS403で生成される信号を示している。図6d3は、図6bに示す離散フーリエ変換した信号に対して、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数対象周波数であって、且つ、±10cycle/mmの範囲の周波数について、その周波数スペクトルを50削減した信号を示す図である。即ち、図6d2は、周波数スペクトルを0にしたが、図6d3では、50削減した信号を示している。図7e1は、図6d1に示す信号に対して、逆離散フーリエ変換した信号を示す図である。図7e2は、図6d2に示す信号に対して、逆離散フーリエ変換した信号を示す図であり、図4のS404で生成される信号を示している。図7e3は、図6d3に示す信号に対して、逆離散フーリエ変換した信号を示す図である。図7f1は、図7e1に示す信号のヒストグラムを示す図である。図7f2は、図7e2に示す信号のヒストグラムを示す図である。図7f3は、図7e3に示す信号のヒストグラムを示す図である。図7f1~f3から、図7f1に示す度数は、図7f2、図7f3に示す度数に比べ、小さい乱数値に偏りが生じていることが分かる。一方、図7f2、図7f3に示す度数は、図6cに示す度数よりも乱数の一様分布性は劣るものの、図7f1に示す度数よりも、乱数の一様分布性を保っており、乱数として好ましいことが分かる。図8g1は、図7e1に示す信号をテーブル化して100回繰り返した場合のフーリエ変換した信号を示す図である。図8g2は、図7e2に示す信号をテーブル化して100回繰り返した場合のフーリエ変換した信号を示す図である。図8g3は、図7e3に示す信号をテーブル化して100回繰り返した場合のフーリエ変換した信号を示す図である。図8hは、図6aに示す信号をテーブル化して100回繰り返した場合のフーリエ変換した信号を示す図である。図8g1~g3から、図8g1に示す信号は、図8g2、図8g3に示す信号に比べ、低周波成分マイナスの場合は絶対値のスペクトルが高くなっていることが分かる。即ち、乱数の低周波成分は、画像のムラとして知覚されるので、この場合には、画像不良を招くことが分かる。一方、図8g2、図8g3に示す信号は、図8hに示す信号よりも、ホワイトノイズ特性としては劣るものの、図8g1に示す信号よりも、ホワイトノイズ特性を保っており、乱数として好ましいことが分かる。このように、図6bに示す離散フーリエ変換した信号に対して、対象周波数であって、且つ、±10cycle/mmの範囲の周波数について、その周波数スペクトルを0にした信号を、逆離散フーリエ変換して生成される信号は、図7f2に示す通り、乱数の一様分布性を保っている。また、この逆離散フーリエ変換して生成される信号は、図8g2に示す通り、ホワイトノイズ特性を保っている。よって、乱数をテーブル化して使用しても、また、記録媒体の搬送誤差や、印刷ヘッドの取り付け誤差があっても、バンディングが抑制された高品質な画質で印刷できる。また、図6bに示す離散フーリエ変換した信号に対して、対象周波数であって、且つ、±10cycle/mmの範囲の周波数について、その周波数スペクトルを50削減した信号を、逆離散フーリエ変換して生成される信号は、図7f3に示す通り、乱数の一様分布性を保っている。また、この逆離散フーリエ変換して生成される信号は、図8g3に示す通り、ホワイトノイズ特性を保っている。よって、乱数をテーブル化して使用しても、また、記録媒体の搬送誤差や、印刷ヘッドの取り付け誤差があっても、バンディングが抑制された高品質な画質で印刷できる。図9は、CPU11が実行する誤差拡散処理を示すフローチャートである。この処理は外部装置から印刷指示が入力された場合に、CPU11が実行する印刷処理に含まれる処理である。具体的には、この誤差拡散処理を実行する前に、CPU11は、入力された画像データを画像データメモリ13aに格納する。また、CPU11は、誤差拡散パラメータメモリ12bから誤差拡散パラメータを読み出し、誤差拡散パラメータ設定メモリ13bに設定する。更に、CPU11は、カラープロファイルメモリ12cからカラープロファイルを読み出し、カラープロファイル設定メモリ13cに設定する。CPU11は、カラープロファイル設定メモリ13cに設定されたカラープロファイルを用いて、画像データメモリ13aに格納されている画像データに色変換処理を施し、画像データに含まれるR,G,Bの画素値を、C,M,Y,Kの画素値に変換する。そして、CPU11は、色変換済みの画像データに対して、誤差拡散パラメータ設定メモリ13bに設定した誤差拡散パラメータを用いて誤差拡散処理を施す。尚、CPU11は、誤差拡散処理を施して生成した印刷データを、ASIC15に供給し、画像を印刷させる。この誤差拡散処理において、CPU11は、注目画素を選択しS902、補正入力値I’を決定するS904。補正入力値I’は、注目画素の画素値に、周辺画素の誤差eを加算した値である。尚、誤差拡散処理における補正入力値I’の算出方法は公知であるため詳細な説明は省略する。次に、CPU11は、基本閾値T設定処理を実行するS906。ここで、図10aを参照して、基本閾値T設定処理について説明する。図10aは、基本閾値T設定処理を示すフローチャートである。CPU11は、乱数テーブル14aからカウンタ値nが示す乱数を取得しS1001、閾値に乱数テーブル14aから取得した乱数を加算するS1002。そして、CPU10は、カウンタ値nに1を換算しS1003、カウンタ値nが乱数テーブルのサイズを超えたかを判断しS1004、超えていなければS1004:No、本処理を終了し、超えていればS1004:Yes、カウンタ値nを1に設定し、本処理を終了する。こうして、基本閾値T設定処理によってS906、閾値に乱数テーブル14aで設定されている乱数が加算され、基本閾値Tが設定される。再び、図9のフローチャートに戻り、説明を続ける。CPU11は、誤差拡散パラメータ設定メモリ13bに設定された誤差拡散パラメータのうち、threLに乱数Tを加算した値を、大ドット閾値として決定するS908。そして、CPU11は、補正入力値I’が大ドット閾値より大きいか否かを判断するS912。その結果、CPU11は、補正入力値I’が大ドット閾値より大きいと判断した場合S912:Yes、注目画素の値を、大ドットの形成を示す値とし、補正入力値I’から大ドット相対濃度DotLを減算した値を、注目画素の誤差eとして算出するS914。CPU11は、算出した誤差eを誤差バッファ図示せずに記憶しS930、最終画素か否かを判断するS932。S932の判断が否定される場合S932:No、CPU11はS902から処理を繰り返す。一方、CPU11は、補正入力値I’が大ドット閾値以下であると判断される場合S912:No、誤差拡散パラメータ設定メモリ13bに設定された誤差拡散パラメータのうち、threMに乱数Tを加算した値を、中ドット閾値として決定するS916。そして、CPU11は、補正入力値I’が中ドット閾値より大きいか否かを判断するS918。その結果、CPU11は、補正入力値I’が中ドット閾値より大きいと判断される場合S918:Yes、注目画素の値を、中ドットの形成を示す値とし、補正入力値I’から中ドット相対濃度DotMを減算した値を、注目画素の誤差eとして算出するS920。CPU11は、算出した誤差eを誤差バッファ図示せずに記憶しS930、最終画素か否かを判断するS932。S932の判断が否定される場合S932:No、CPU11はS902から処理を繰り返す。一方、CPU11は、補正入力値I’が中ドット閾値以下であると判断される場合S918:No、誤差拡散パラメータ設定メモリ13bに設定された誤差拡散パラメータのうち、threSに乱数Tを加算した値を、小ドット閾値として決定するS922。そして、CPU11は、補正入力値I’が小ドット閾値より大きいか否かを判断するS924。その結果、CPU11は、補正入力値I’が小ドット閾値より大きいと判断される場合S924:Yes、注目画素の値を小ドットの形成を示す値とし、補正入力値I’から小ドット相対濃度DotSを減算した値を、注目画素の誤差eとして算出するS926。CPU11は、算出した誤差eを誤差バッファ図示せずに記憶しS930、最終画素か否かを判断するS932。S932の判断が否定される場合S932:No、CPU11はS902から処理を繰り返す。一方、CPU11は、補正入力値I’が小ドット閾値以下であると判断される場合S924:No、注目画素の値をドットの非形成を示す値とし、補正入力値I’を誤差eとするS928。CPU11は、算出した誤差eを誤差バッファ図示せずに記憶しS930、最終画素か否かを判断するS932。S932の判断が否定される場合S932:No、CPU11はS902から処理を繰り返す。このようにして処理を繰り返すうちに、最終画素まで処理を終了するとS930:Yes、CPU11は誤差拡散処理を終了する。このように、上述した誤差拡散処理を含む印刷処理では、図4bに示す乱数テーブル作成処理によって作成され、乱数テーブル14aに記憶されている乱数テーブルが使用されて画像データを生成するので、乱数をテーブル化しても、また、記録媒体の搬送誤差や、印刷ヘッドの取り付け誤差があっても、バンディングが抑制された高品質な画質で印刷できる。次に、図10bを参照して本発明の第2実施形態について説明する。第2実施形態では、図10bに示す2次元テーブル作成処理で作成したテーブルを乱数テーブル14aに記憶する。上述した第1実施形態では、X方向用紙Pの搬送方向と直交する方向に関する1次元乱数テーブルを作成する場合について説明したが、第2実施携帯では、X方向と、X方向に直交するY方向との2次元乱数テーブルを作成する。図10bは、2次元乱数テーブル作成処理を示すフローチャートである。生成装置は、2次元の一様分布乱数を生成するS1010。尚、この2次元の一様分布乱数は、第1実施形態の場合と同様に、公知のアルゴリズムであるメルセンヌツイスターを使って、2次元の一様分布乱数であるホワイトノイズを生成する。そして、CPU11は、その生成したホワイトノイズを2次元離散フーリエ変換するS1011。生成装置は、2次元離散フーリエ変換された信号値から、X方向とY方向との空間周波数が±10cycle/mmの範囲で、且つ、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数について、その周波数スペクトルを0にするS1012。生成装置は、S1012の処理から生成される信号値を、2次元逆離散フーリエ変換しS1013、虚数値の絶対値をとり、実数値に変換するS1014。また、S1014の処理を経て出力される信号値の最小値が0、最大値が1となる範囲で線形変換しS1015、更に、S1015の処理を経て出力される信号値について、図5で説明したのと同様に、整数値変換処理を実行しS1016、本処理を終了する。この2次元乱数テーブル作成処理で作成した2次元乱数テーブルは、第1実施形態と同様に、乱数テーブル14aに記憶され、2次元乱数テーブルに記憶されている乱数が誤差拡散処理の閾値加えられる。このように、2次元乱数テーブル作成処理で作成した2次元乱数テーブルに記憶されている乱数を誤差拡散処理の閾値に加えても、第1実施形態の場合と同様に、誤差拡散処理の閾値に加えられる乱数は、X方向と、Y方向との2次元状のホワイトノイズから、Y方向の空間周波数が±10cycle/mmの範囲に含まれ周波数についての周波数スペクトルを、0にして、又は、削減して生成された乱数である。よって、2次元状のホワイトノイズ特性を極力維持した乱数を用いて、誤差拡散処理で使用する閾値を変動することができる。即ち、Y方向についても、X方向と同様に、0、又は、削減する周波数スペクトルを、視覚上無視できる±10cycle/mmの範囲に限定することで、2次元状のホワイトノイズの高周波成分が低下するのを抑制できる。よって、2次元乱数テーブルを用いる場合でも、バンディングが抑制された高品質な画質で印刷できる。上記の実施形態において、プリンタ1が画像処理装置の一例である。印刷ヘッド2が印刷手段の一例である。CPU11が実行する図9の誤差拡散処理が誤差拡散処理の一例である。CPU11が実行する図9のS906の処理が変動手段の一例である。乱数テーブル14aが記憶手段の一例である。制御プログラム12aが画像処理プラグラムの一例である。CPU11がコンピュータの一例である。生成装置が実行する図4bのS401、図10bのS1010が生成ステップの一例である。生成装置が実行する図4bのS402、図10bのS1011がフーリエ変換ステップの一例である。生成装置が実行する図4bのS403、図10bのS1013がフィルタステップの一例である。生成装置が実行する図4bのS404、図10bのS1014が逆フーリエ変換ステップの一例である。以上、実施形態に基づき本発明を説明したが、本発明は上述した実施形態に何ら限定されるものではなく、本発明の趣旨を逸脱しない範囲内で種々の改良変更が可能であることは容易に推察できるものである。上述した実施形態では、周波数スペクトルを0、又は、50に削減する周波数を、±10cycle/mmの範囲に含まれる周波数に制限する場合について説明したが、必ずしも、かかる制限をする必要はない。例えば、図4に示すS403の処理において、離散フーリエ変換された信号値から、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数対象周波数について、その周波数スペクトルを0にしても良い。かかる場合には、上述した実施形態よりは、ホワイトノイズ特性が低減するものの、上述した実施形態と同様な効果を奏することができる。また、この場合は、周波数スペクトルを0にする周波数は、1つでも、複数であっても良い。また、上述した実施形態では、対象周波数の周波数スペクトルを0にする、または、50に削減する場合について説明したが、対象周波数の周波数スペクトルを、削減さえすれば、上述した実施形態と同様な効果を奏することができる。また、上述した実施形態では、X方向のノズルの配列周期に相当する周波数を、±2n倍nは0以上の整数した周波数対象周波数について、その周波数スペクトルを0にする場合について説明したが、周波数スペクトルを0、または、削減する周波数としては、対象周波数を中心とした±0.2cycle/mmの範囲に含まれる全部、または、1以上の周波数であっても良い。"