(57)【特許請求の範囲】
【請求項１】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データから求める前記出力手段における前記画像処理信号のＤＣ成分ごとの画像再現特性に基づいて前記処理手段の前記空間周波数処理係数を決定する係数決定手段を備え、該係数決定手段は、前記テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定する測定手段と、該測定手段で測定された前記領域の空間周波数成分の測定結果と前記テストパターン画像データの空間周波数成分とを比較する比較手段と、前記比較手段の比較結果から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記処理手段の前記空間周波数処理係数を決定する空間周波数処理係数決定手段を有し、前記処理手段は、該空間周波数処理係数に基づいて前記画像処理信号のＤＣ成分ごとに前記注目画素の値に加える変更の程度を変化させることを特徴とする画像処理装置。
【請求項２】
前記画像処理信号のＤＣ成分ごとに付与すべき前記所定の特性は、前記出力装置の持つ画像再現特性を前記画像処理信号のＤＣ成分ごとに変化しないとする理想特性にするためのものであることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データから求める前記出力手段の画像再現特性に基づいて前記処理手段の前記空間周波数処理係数を決定する係数決定手段を備え、前記空間周波数処理係数は、基本処理係数と補助処理係数を含み、前記係数決定手段は、前記テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定する測定手段と、該測定手段で測定された前記領域の空間周波数成分の測定結果と前記テストパターン画像データの空間周波数成分とを比較する比較手段と、前記比較手段の比較結果から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記処理手段の前記空間周波数処理係数を決定する空間周波数処理係数決定手段を有し、前記空間周波数処理係数決定手段は、前記出力装置に入力する前記画像処理信号の持つ空間周波数特性に所定の特性を付与するために前記基本処理係数を決定する基本処理係数決定手段と、該基本係数決定手段で決定される基本処理係数で変化する画像処理信号の持つ空間周波数特性にさらに加えるべき処理の特性を付与するための前記補助処理係数を決定する補助処理係数決定手段を有していることを特徴とする画像処理装置。
【請求項４】
前記補助処理係数は、補助空間周波数処理係数と補助ゲイン係数を含み、前記処理係数決定手段は、前記画像処理信号の持つ空間周波数特性に所定の特性を付与するための前記補助空間周波数処理係数を決定する補助空間周波数処理係数決定手段と、該補助空間周波数処理係数決定手段で決定される補助空間周波数処理係数に基づき変更する画像処理信号に掛け合わせるべきゲインを得るための補助ゲイン係数を決定するゲイン係数決定手段を含むことを特徴とする請求項３に記載の画像処理装置。
【請求項５】
前記基本処理係数は、基本空間周波数処理係数と基本ゲイン係数を含み、前記処理係数決定手段は、前記画像処理信号の持つ空間周波数特性に所定の特性を付与するための前記基本空間周波数処理係数を決定する基本空間周波数処理係数決定手段と、該基本空間周波数処理係数決定手段で決定される基本空間周波数処理係数に基づき変更する画像処理信号に掛け合わせるべきゲインを得るための基本ゲイン係数を決定するゲイン係数決定手段を含むことを特徴とする請求項３に記載の画像処理装置。
【請求項６】
前記ゲイン係数決定手段は、前記画像処理信号の局所的なコントラストに対応する基本ゲイン係数または補助ゲイン係数を決定することを特徴とする請求項４または請求項５に記載の画像処理装置。
【請求項７】
前記ゲイン係数決定手段は、前記画像処理信号の局所的な平均値に対応する基本ゲイン係数または補助ゲイン係数を決定することを特徴とする請求項４または請求項５に記載の画像処理装置。
【請求項８】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項９】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１０】
前記空間周波数処理係数は、基本処理係数と補助処理係数を含み、前記基本処理係数は、前記出力装置に入力する前記画像処理信号の持つ空間周波数特性に所定の特性を付与するためのものであり、前記補助処理係数は、前記基本処理係数で変化する画像処理信号の持つ空間周波数特性にさらに加えるべき処理の特性を付与するためのものであり、前記処理手段は、前記画像処理信号の前記注目画素の値を前記基本処理係数に基づき変更する基本処理手段と、前記注目画素の値を前記補助処理係数に基づき変更する補助処理手段を含み、それぞれの変更された注目画素の値に基づいて前記出力手段に与える画像処理信号を作成することを特徴とする請求項８または請求項９に記載の画像処理装置。
【請求項１１】
前記補助処理係数は、補助空間周波数処理係数と補助ゲイン係数を含み、前記補助処理手段は、前記補助空間周波数処理係数に基づき画像処理信号に変更を加え、さらに前記特徴量抽出手段により求めた特徴と前記補助ゲイン係数から得られるゲインに基づいた変更を加えることを特徴とする請求項１０に記載の画像処理装置。
【請求項１２】
前記基本処理係数は、基本空間周波数処理係数と基本ゲイン係数を含み、前記基本処理手段は、前記基本空間周波数処理係数に基づき画像処理信号に変更を加え、さらに前記特徴量抽出手段により求めた特徴と前記基本ゲイン係数から得られるゲインに基づいた変更を加えることを特徴とする請求項１０に記載の画像処理装置。
【請求項１３】
前記空間周波数処理係数は、前記テストパターン画像を実際に測定した空間周波数成分と前記テストパターン画像データの空間周波数成分との比較結果とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて算出されたものであることを特徴とする請求項９に記載の画像処理装置。
【請求項１４】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記画像処理信号の空間周波数成分ごとのゲインを決定するゲイン決定手段と、前記入力手段に入力された画像処理信号の空間周波数成分ごとに前記ゲイン決定手段で決定されたそれぞれの前記ゲインに基づいて前記注目画素の値に変更を加え得られたそれぞれの値から変更された画像処理信号を得る処理手段と、該処理手段によって得られた前記変更された画像処理信号の画像を出力する出力手段を備えており、前記画像処理信号の空間周波数成分ごとの前記ゲインは、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１５】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記画像処理信号の空間周波数成分ごとのゲインを決定するゲイン決定手段と、前記入力手段に入力された画像処理信号の空間周波数成分ごとに前記ゲイン決定手段で決定されたそれぞれの前記ゲインに基づいて前記注目画素の値に変更を加え得られたそれぞれの値から変更された画像処理信号を得る処理手段と、該処理手段によって得られた前記変更された画像処理信号の画像を出力する出力手段を備えており、前記画像処理信号の空間周波数成分ごとの前記ゲインは、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１６】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、前記入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記処理手段の前記空間周波数処理係数を設定する処理係数決定手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１７】
注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、前記入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記処理手段の前記空間周波数処理係数を設定する処理係数決定手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１８】
均等色空間の色座標系で表現された値を含む画像処理信号を入力する第１の入力手段と、該第１の入力手段に入力された画像処理信号を出力手段の制御に適した前記色空間と異なる第２の色空間の色座標系に変換する色変換手段と、該色変換手段により色変換された注目画素および該注目画像に対する周辺画素の画像信号の値を含む画像処理信号を入力する第２の入力手段と、該第２の入力手段に入力された前記画像信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、前記第２の入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像のデータの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項１９】
均等色空間の色座標系で表現された値を含む画像処理信号を入力する第１の入力手段と、該第１の入力手段に入力された画像処理信号を出力手段の制御に適した前記色空間と異なる第２の色空間の色座標系に変換する色変換手段と、該色変換手段により色変換された注目画素および該注目画像に対する周辺画素の画像信号の値を含む画像処理信号を入力する第２の入力手段と、該第２の入力手段に入力された前記画像信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、前記第２の入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像のデータの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とする画像処理装置。
【請求項２０】
前記空間周波数処理係数は、前記テストパターン画像を実際に測定した空間周波数成分と前記テストパターン画像データの空間周波数成分との比較結果とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて算出されたものであることを特徴とする請求項１６ないし請求項１９のいずれか１項に記載の画像処理装置。
【請求項２１】
特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段により出力し出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的なコントラストを特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程を有することを特徴とする画像処理方法。
【請求項２２】
特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段により出力し出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的な平均値を特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程を有することを特徴とする画像処理方法。
【請求項２３】
前記空間周波数処理係数を決定する工程は、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記空間周波数処理係数を算出することを特徴とする請求項２１または請求項２２に記載の画像処理方法。
【請求項２４】
特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段に出力させ出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的なコントラストを特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程をコンピュータに実行させるプログラムを記憶した前記コンピュータが読み取り可能な記憶媒体。
【請求項２５】
特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段に出力させ出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的な平均値を特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程をコンピュータに実行させるプログラムを記憶した前記コンピュータが読み取り可能な記憶媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、画像データに対して空間周波数処理を施して出力する画像処理に関するものである。
【０００２】
【従来の技術】
複数の出力機器が接続されたデスクトップパブリシング（ＤＴＰ）システム等の画像処理システムを用いてカラー画像の編集および出力を行なっている現状において、出力機器によらず同一の出力画像を得たいという要望がある。しかし、画像出力装置の再現特性は、各装置ごとにばらつきが大きく、そのまま出力したのでは同一の特性を有する出力画像が得られない。
【０００３】
このような個々の画像出力装置の再現特性のばらつきに対して、色の再現に関しては、画像出力装置によらず同一の色再現出力画像が得られるようになりつつある。例えば、カラーマネージメントシステム（ＣＭＳ）を用いて、装置に依存しないカラー信号から各装置固有の色信号への色変換を行なうことにより、出力機器によらず同一の色再現出力画像が得られる。この際の色変換に用いられる色変換係数は、例えば、特開平２−２８９３６７号公報に述べられているように、複数のカラーパッチの出力画像を測色することにより算出することが可能である。算出された色変換係数は、装置特性情報として、例えばＩＣＣプロファイルの形成で記述され、システムに保持される。
【０００４】
しかし、それぞれの出力機器に応じた色変換を施した画像データを作成し、出力しても、出力機器ごとに鮮鋭度の異なる出力画像となってしまう。鮮鋭度を比較する尺度として、画像処理の分野では各周波数ごとにコントラストをどのくらい伝達できるかを示すＣＴＦ（コントラストトランスファーファンクション）を用いる。以下の説明では、このＣＴＦを各出力装置の空間周波数特性とする。
【０００５】
現状においては、各画像出力装置の空間周波数特性の影響による出力画像の変化をカンと経験から定性的に予測し、画像処理アプリケーションソフトウェアエアを用いて、画像データに対してシャープネス処理や平滑化処理などの空間周波数処理を繰り返し行ない、希望する出力画像を得ている。
【０００６】
また、あらかじめプリンタ内部に空間周波数処理部を内蔵しているカラープリンタにおいては、その開発段階において、空間周波数処理部の処理係数を決定している。その際には、処理係数をカンと経験に基づいて変化させながら、複数種類の基準画像データの出力画像をカラープリンタで繰り返し作成し、出力画像の鮮鋭度を観察することにより、処理係数を決定している。
【０００７】
【発明が解決しようとする課題】
本発明は、上述した事情に鑑みてなされたもので、画像出力装置の空間周波数伝達特性を補正するために行なう画像データに対する空間周波数処理において、その処理係数をカンや経験によらずに決定できる処理係数算出手段によって各画像出力装置ごとに算出された空間処理係数を用い、画像出力装置によらず同一の鮮鋭度の出力画像を得ることが可能な画像処理装置及び画像処理方法と、そのような画像処理プログラムを記憶した記憶媒体を提供することを目的とするものである。
【０００８】
【課題を解決するための手段】
請求項１に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データから求める前記出力手段における前記画像処理信号のＤＣ成分ごとの画像再現特性に基づいて前記処理手段の前記空間周波数処理係数を決定する係数決定手段を備え、該係数決定手段は、前記テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定する測定手段と、該測定手段で測定された前記領域の空間周波数成分の測定結果と前記テストパターン画像データの空間周波数成分とを比較する比較手段と、前記比較手段の比較結果から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記処理手段の前記空間周波数処理係数を決定する空間周波数処理係数決定手段を有し、前記処理手段は、該空間周波数処理係数に基づいて前記画像処理信号のＤＣ成分ごとに前記注目画素の値に加える変更の程度を変化させることを特徴とするものである。
【０００９】
請求項２に記載の発明は、請求項１に記載の画像処理装置において、前記画像処理信号のＤＣ成分ごとに付与すべき前記所定の特性は、前記出力装置の持つ画像再現特性を前記画像処理信号のＤＣ成分ごとに変化しないとする理想特性にするためのものであることを特徴とするものである。
【００１０】
請求項３に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データから求める前記出力手段の画像再現特性に基づいて前記処理手段の前記空間周波数処理係数を決定する係数決定手段を備え、前記空間周波数処理係数は、基本処理係数と補助処理係数を含み、前記係数決定手段は、前記テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定する測定手段と、該測定手段で測定された前記領域の空間周波数成分の測定結果と前記テストパターン画像データの空間周波数成分とを比較する比較手段と、前記比較手段の比較結果から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記処理手段の前記空間周波数処理係数を決定する空間周波数処理係数決定手段を有し、前記空間周波数処理係数決定手段は、前記出力装置に入力する前記画像処理信号の持つ空間周波数特性に所定の特性を付与するために前記基本処理係数を決定する基本処理係数決定手段と、該基本係数決定手段で決定される基本処理係数で変化する画像処理信号の持つ空間周波数特性にさらに加えるべき処理の特性を付与するための前記補助処理係数を決定する補助処理係数決定手段を有していることを特徴とするものである。
【００１１】
請求項４に記載の発明は、請求項３に記載の画像処理装置において、前記補助処理係数は、補助空間周波数処理係数と補助ゲイン係数を含み、前記処理係数決定手段は、前記画像処理信号の持つ空間周波数特性に所定の特性を付与するための前記補助空間周波数処理係数を決定する補助空間周波数処理係数決定手段と、該補助空間周波数処理係数決定手段で決定される補助空間周波数処理係数に基づき変更する画像処理信号に掛け合わせるべきゲインを得るための補助ゲイン係数を決定するゲイン係数決定手段を含むことを特徴とするものである。
【００１２】
請求項５に記載の発明は、請求項３に記載の画像処理装置において、前記基本処理係数は、基本空間周波数処理係数と基本ゲイン係数を含み、前記処理係数決定手段は、前記画像処理信号の持つ空間周波数特性に所定の特性を付与するための前記基本空間周波数処理係数を決定する基本空間周波数処理係数決定手段と、該基本空間周波数処理係数決定手段で決定される基本空間周波数処理係数に基づき変更する画像処理信号に掛け合わせるべきゲインを得るための基本ゲイン係数を決定するゲイン係数決定手段を含むことを特徴とするものである。
【００１３】
請求項６に記載の発明は、請求項４または請求項５に記載の画像処理装置において、前記ゲイン係数決定手段は、前記画像処理信号の局所的なコントラストに対応する基本ゲイン係数または補助ゲイン係数を決定することを特徴とするものである。
【００１４】
請求項７に記載の発明は、請求項４または請求項５に記載の画像処理装置において、前記ゲイン係数決定手段は、前記画像処理信号の局所的な平均値に対応する基本ゲイン係数または補助ゲイン係数を決定することを特徴とするものである。
【００１５】
請求項８に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００１６】
請求項９に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００１７】
請求項１０に記載の発明は、請求項８または請求項９に記載の画像処理装置において、前記空間周波数処理係数は、基本処理係数と補助処理係数を含み、前記基本処理係数は、前記出力装置に入力する前記画像処理信号の持つ空間周波数特性に所定の特性を付与するためのものであり、前記補助処理係数は、前記基本処理係数で変化する画像処理信号の持つ空間周波数特性にさらに加えるべき処理の特性を付与するためのものであり、前記処理手段は、前記画像処理信号の前記注目画素の値を前記基本処理係数に基づき変更する基本処理手段と、前記注目画素の値を前記補助処理係数に基づき変更する補助処理手段を含み、それぞれの変更された注目画素の値に基づいて前記出力手段に与える画像処理信号を作成することを特徴とするものである。
【００１８】
請求項１１に記載の発明は、請求項１０に記載の画像処理装置において、前記補助処理係数は、補助空間周波数処理係数と補助ゲイン係数を含み、前記補助処理手段は、前記補助空間周波数処理係数に基づき画像処理信号に変更を加え、さらに前記特徴量抽出手段により求めた特徴と前記補助ゲイン係数から得られるゲインに基づいた変更を加えることを特徴とするものである。
【００１９】
請求項１２に記載の発明は、請求項１０に記載の画像処理装置において、前記基本処理係数は、基本空間周波数処理係数と基本ゲイン係数を含み、前記基本処理手段は、前記基本空間周波数処理係数に基づき画像処理信号に変更を加え、さらに前記特徴量抽出手段により求めた特徴と前記基本ゲイン係数から得られるゲインに基づいた変更を加えることを特徴とするものである。
【００２０】
請求項１３に記載の発明は、請求項９に記載の画像処理装置において、前記空間周波数処理係数は、前記テストパターン画像を実際に測定した空間周波数成分と前記テストパターン画像データの空間周波数成分との比較結果とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて算出されたものであることを特徴とするものである。
【００２１】
請求項１４に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記画像処理信号の空間周波数成分ごとのゲインを決定するゲイン決定手段と、前記入力手段に入力された画像処理信号の空間周波数成分ごとに前記ゲイン決定手段で決定されたそれぞれの前記ゲインに基づいて前記注目画素の値に変更を加え得られたそれぞれの値から変更された画像処理信号を得る処理手段と、該処理手段によって得られた前記変更された画像処理信号の画像を出力する出力手段を備えており、前記画像処理信号の空間周波数成分ごとの前記ゲインは、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２２】
請求項１５に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、該入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記画像処理信号の空間周波数成分ごとのゲインを決定するゲイン決定手段と、前記入力手段に入力された画像処理信号の空間周波数成分ごとに前記ゲイン決定手段で決定されたそれぞれの前記ゲインに基づいて前記注目画素の値に変更を加え得られたそれぞれの値から変更された画像処理信号を得る処理手段と、該処理手段によって得られた前記変更された画像処理信号の画像を出力する出力手段を備えており、前記画像処理信号の空間周波数成分ごとの前記ゲインは、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２３】
請求項１６に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、前記入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記処理手段の前記空間周波数処理係数を設定する処理係数決定手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２４】
請求項１７に記載の発明は、画像処理装置において、注目画素および該注目画素に対する周辺画素の画像信号の値を含む画像処理信号を入力する入力手段と、前記入力手段に入力された画像処理信号と空間周波数処理係数に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段と、前記入力手段に入力された前記画像処理信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、該特徴量抽出手段で求めた特徴に基づいて前記処理手段の前記空間周波数処理係数を設定する処理係数決定手段を備えており、前記空間周波数処理係数は、テストパターン画像データに基づいて前記出力手段で作成されたテストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像データの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２５】
請求項１８に記載の発明は、画像処理装置において、均等色空間の色座標系で表現された値を含む画像処理信号を入力する第１の入力手段と、該第１の入力手段に入力された画像処理信号を出力手段の制御に適した前記色空間と異なる第２の色空間の色座標系に変換する色変換手段と、該色変換手段により色変換された注目画素および該注目画像に対する周辺画素の画像信号の値を含む画像処理信号を入力する第２の入力手段と、該第２の入力手段に入力された前記画像信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的なコントラストを特徴として求める特徴量抽出手段と、前記第２の入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像のデータの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２６】
請求項１９に記載の発明は、画像処理装置において、均等色空間の色座標系で表現された値を含む画像処理信号を入力する第１の入力手段と、該第１の入力手段に入力された画像処理信号を出力手段の制御に適した前記色空間と異なる第２の色空間の色座標系に変換する色変換手段と、該色変換手段により色変換された注目画素および該注目画像に対する周辺画素の画像信号の値を含む画像処理信号を入力する第２の入力手段と、該第２の入力手段に入力された前記画像信号の注目画素および周辺画素の画像信号の値から該画像処理信号の局所的な平均値を特徴として求める特徴量抽出手段と、前記第２の入力手段に入力された画像処理信号と空間周波数処理係数と前記特徴量抽出手段で求めた特徴に基づいて前記注目画素の値に変更を加える処理手段と、該処理手段によって値の変更を受けた画像処理信号の画像を出力する出力手段を備えており、前記空間周波数処理係数は、テストパターン画像の実際の空間周波数成分を測定し前記テストパターン画像のデータの空間周波数成分との比較によって前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように決定されたものであることを特徴とするものである。
【００２７】
請求項２０に記載の発明は、請求項１６ないし請求項１９のいずれか１項に記載の画像処理装置において、前記空間周波数処理係数は、前記テストパターン画像を実際に測定した空間周波数成分と前記テストパターン画像データの空間周波数成分との比較結果とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて算出されたものであることを特徴とするものである。
【００２８】
請求項２１に記載の発明は、画像処理方法において、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段により出力し出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的なコントラストを特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程を有することを特徴とするものである。
【００２９】
請求項２２に記載の発明は、画像処理方法において、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段により出力し出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的な平均値を特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程を有することを特徴とするものである。
【００３０】
請求項２３に記載の発明は、請求項２１または請求項２２に記載の画像処理方法において、前記空間周波数処理係数を決定する工程は、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係とともに前記出力手段が実効的に許容する入力値の範囲とから目標空間周波数伝達特性を決定し、該目標空間周波数伝達特性に基づいて前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように前記空間周波数処理係数を算出することを特徴とするものである。
【００３１】
請求項２４に記載の発明は、コンピュータが読み取り可能な記憶媒体において、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段に出力させ出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的なコントラストを特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程をコンピュータに実行させるプログラムを記憶したことを特徴とするものである。
【００３２】
請求項２５に記載の発明は、コンピュータが読み取り可能な記憶媒体において、特定の複数の空間周波数成分を有する領域を含むテストパターン画像データを出力手段に出力させ出力画像を作成する工程と、前記出力画像の空間周波数特性を測定する工程と、前記テストパターン画像データの空間周波数特性と前記出力画像の空間周波数特性の測定結果との対応関係から前記テストパターン画像の空間周波数特性に前記出力手段の空間周波数特性を近づけるように空間周波数処理係数を決定する工程と、入力された注目画素および周辺画素の画像信号の値を含む画像処理信号の局所的な平均値を特徴として求める工程と、前記画像処理信号と前記特徴と前記空間周波数処理係数に基づいて前記注目画素の値に変更を加える工程と、変更を加えた前記画像処理信号の画像を出力する工程をコンピュータに実行させることを特徴とするものである。
【００５７】
【発明の実施の形態】
図１は、本発明の実施の一形態を示すブロック図である。図中、１は処理演算部、２は出力部、３は測定部、４は係数算出部、５はテストパターン保持部である。処理演算部１は、入力された画像データに対して各種の処理を施し、出力画像データを作成する。ここでは、少なくとも空間周波数処理を行ない、画像出力装置に応じた空間周波数となるように、入力画像データに対して変換処理を行なう。この時行なわれる空間周波数処理は、係数算出部４で算出された処理係数を用いて行なう。また、処理係数の算出時には、テストパターン保持部５に保持されているテストパターン画像データに対して空間周波数処理以外の必要な処理を施し、出力画像データを作成する。出力部２は、処理演算部１で作成された出力画像データを出力する。
【００５８】
測定部３は、空間周波数処理で用いる処理係数を算出する際に、テストパターン保持部５に保持されているテストパターン画像データに基づいて出力部２から出力されたテストパターン画像の空間周波数特性を測定する。
【００５９】
係数算出部４は、テストパターン保持部５に保持されているテストパターン画像データの空間周波数特性を求め、また、測定部３で測定したテストパターン画像の空間周波数特性を得て、両者の対応関係から空間周波数処理で用いる処理係数を算出する。
【００６０】
テストパターン保持部５は、テストパターン画像データを保持する。テストパターン画像データは、例えば空間周波数の異なる１以上のパッチ画像により構成することができる。パッチ画像の１つとして、ＤＣ成分、すなわち空間周波数が０のパッチ画像を含めることができる。
【００６１】
空間周波数処理で用いる処理係数を算出する際には、テストパターン保持部５に保持されているテストパターン画像データに対して処理演算部１で空間周波数処理以外の必要な処理、例えばカラー処理等を行なった後、出力部２から出力する。出力されたテストパターン画像の空間周波数特性を測定部３で測定する。係数算出部４において、測定部３で測定した空間周波数特性と、テストパターン保持部５に保持されているテストパターン画像データの空間周波数特性とを対応づけ、２つの空間周波数特性を近づけるような空間周波数処理を行なうための処理係数を算出する。算出された処理係数は処理演算部１に渡される。
【００６２】
画像データの出力処理は、このようにして空間周波数処理を行なうための処理係数が求められた後に行なう。処理演算部１は、入力された画像に対して係数算出部４で算出された処理係数を用いた空間周波数処理を含め、各種の画像処理を行なって出力画像データを作成し、出力部２から出力する。このとき、空間周波数処理のための処理係数が既に算出されて処理演算部１に設定されているので、処理演算部１および出力部２のみで動作でき、測定部３、係数算出部４、テストパターン保持部５はなくてもよい。
【００６３】
図２は、処理演算部の構成の一例を示すブロック図である。図中、１１は色空間変換部、１２は空間周波数処理部、１３は階調変換部である。ここでは、具体例として、カラー画像データに対する処理を示している。入力されるカラー画像データは、一例として、入出力装置に依存しないＬ* ａ* ｂ* 各８Ｂｉｔのデータであるものとする。また、出力部２は、ＹＭＣ各８Ｂｉｔのカラー画像を受け取り、画像を出力するものとする。もちろん、入力、出力とも、他の表色系を用いたり、ビット数が異なっていてもよい。
【００６４】
色空間変換部１１は、Ｌ* ａ* ｂ* 各８Ｂｉｔの画像データから、ＹＭＣ各８Ｂｉｔの画像データへの変換を行なう。このＹＭＣカラー信号は、ａ* ＝ｂ* ＝０のときＹ＝Ｍ＝Ｃ、ａ* ＝ｂ* ＝０かつＬ* ＝１００のときＹ，Ｍ，Ｃ＝０、ａ* ＝ｂ* ＝０かつＬ* ＝０のときＹ，Ｍ，Ｃ＝２５５であり、ａ* ＝ｂ* ＝０のときＹ，Ｍ，Ｃの値はＬ* 値にリニアな関係となるように階調設定がされている。このような色変換処理は、例えばダイレクト・ルックアップ・テーブルを用いて変換を行なうことができる。
【００６５】
次に、空間周波数処理部１２においてＹ，Ｍ，Ｃの各カラー画像データに対して空間周波数処理を行なう。空間周波数処理としては、デジタルフィルタ処理が一般的である。図３は、デジタルフィルタ処理の一例の説明図である。いま、図３（Ｂ）に示すような３×３のフィルタ係数マトリクスを考える。空間周波数処理の場合、このフィルタ係数マトリクスの係数の和が１となるように設定されている。図３（Ａ）に示す画像データが入力されたとき、各画素を順次注目画素として、周囲の画素とともに３×３画素について積和演算を行なう。図３に示した例では、中央の値３の画素を注目画素とし、その注目画素を中心とする３×３画素を用いて演算を行なう。例えば、１×（−０．２）＋２×（−０．１）＋１×（−０．２）＋３×２．２＋２×（−０．２）の計算を行ない、値５．８を得る。得られた値を注目画素の新たな値とする。このような演算処理を、各画素を注目画素として行なうことにより、図３（Ｃ）に示すような空間周波数処理後の画像データが得られる。
【００６６】
ここで、図３（Ｂ）に示すフィルタ係数マトリクスは、係数算出部４において、出力部２の再現特性に応じて算出され、設定される。フィルタ係数マトリクスは、ＹＭＣ各軸ごとに、それぞれ１つ以上設定される。後述するように、入力画像データの平均値やコントラスト等に応じた複数のフィルタ係数マトリクスを設定し、これらを入力された画像データの空間周波数特性に応じて切り換えて用いるように構成することができる。
【００６７】
最後に、階調変換部１３において各々ルックアップ・テーブルにより階調特性を調整し、例えば、面積変調により画像を形成するカラー画像形成装置においては、出力カバレッジにリニアな画像信号である、Ｙ’Ｍ’Ｃ’各８Ｂｉｔに変換して出力部２に出力する。階調変換部１３は、ＹＭＣ各軸に対応したそれぞれのルックアップ・テーブルで構成することができる。
【００６８】
上述の色空間変換部１１におけるダイレクト・ルックアップ・テーブルの処理係数、および、階調変換部１３におけるルックアップ・テーブルの処理係数は、出力装置色変換情報としてあらかじめ出力部２の出力装置ごとに保持されているものとする。使用する出力装置に対応した出力装置色変換情報を用いて処理を行なうことにより、入出力装置に依存しない形式で表現された画像データから出力装置によらず同一色再現の出力画像を得ることが可能となっている。また、空間周波数処理部１２におけるデジタルフィルタの処理係数も、係数算出部４で出力装置に対応づけて算出しておくことにより、出力装置によらず同一の空間周波数特性を有する出力画像を得ることが可能となっている。なお、空間周波数処理部１２におけるデジタルフィルタの処理係数も、係数算出部４で算出後、出力装置ごとに保持しておけばよい。
【００６９】
図２に示した構成では、テストパターン保持部５に保持されているテストパターン画像データはＹＭＣ各８Ｂｉｔの画像データとして保持されていることを仮定している。そのため、空間周波数処理部１２で用いる処理係数を算出する際には、階調変換部１３の処理を行なって、出力部２からテストパターン画像を出力している。テストパターン画像データがＹＭＣ色空間以外の色空間のデータであれば、色空間の変換を行なってから階調変換部１３に入力することになる。例えば、Ｌ* ａ* ｂ* 各８Ｂｉｔのデータであれば、色空間変換部１１の処理後、空間周波数処理部１２の処理を行なわずに階調変換部１３の処理を行なうように構成すればよい。
【００７０】
空間周波数処理部１２についてさらに説明する。図４は、種々の出力装置における空間周波数応答特性の一例の説明図である。図中、実線は入力された画像データの信号を示し、破線は出力された画像の信号を示す。出力装置は、それぞれの出力方式の特徴や機器の構成などによって、入力された画像データの空間周波数特性と出力された画像の空間周波数特性は、線形に対応しない場合がある。例えば、図４（Ａ）に示すように、同じ空間周波数、同じ振幅であっても、平均の明るさ（以下ＤＣ成分と呼ぶ）の異なる入力画像データに対して出力画像の振幅が異なる場合がある。また、図４（Ｂ）に示すように、同じ空間周波数、同じＤＣ成分で異なる振幅の場合、例えば、小さな振幅の入力画像データに対して振幅の大きな出力画像となったり、大きな振幅の入力画像データに対して振幅の小さな出力画像が得られたりする。この時の入力画像データの振幅と出力画像の振幅は線形に変化するとは限らない。以下の説明では、画像信号からＤＣ成分を取り除いた信号をＡＣ成分と呼ぶことにする。
【００７１】
このような異なる空間周波数応答をする出力装置が存在するため、これらに対応した空間周波数処理を行なうと有効である。図４（Ａ）に示す場合に対応するため、同じ空間周波数で異なるＤＣ成分の入力画像データに対応して補正量を調整し、また、図４（Ｂ）に示す場合に対応するため、同じ空間周波数で異なるＡＣ成分の入力信号に対応して補正量を調整すればよい。
【００７２】
図５は、空間周波数処理部の構成の第１の例を示すブロック図である。図中、２１はハイパスフィルタ、２２はローパスフィルタ、２３，２４はルックアップテーブル、２５は乗算部、２６は加算部、３１は入力画像信号、３２は補正量信号、３３はＡＣ成分ゲイン信号、３４は調整補正量信号、３５はＤＣ成分信号、３６はＤＣ成分ゲイン信号、３７は出力画像信号である。図２に示した空間周波数処理部１２では、ＹＭＣ各画像信号に対してデジタルフィルタ処理を行なう構成として、例えば、図５に示す構成を適用することができる。
【００７３】
ハイパスフィルタ２１、ローパスフィルタ２２は、例えば２次元ＦＩＲフィルタ（ＦＩＲＦ：Ｆｉｎｉｔｅ  Ｉｍｐｕｌｓｅ  Ｒｅｓｐｏｎｓｅ  Ｆｉｌｔｅｒ：有限インパルス応答フィルタ）で構成することができる。２次元ＦＩＲフィルタは、画像データ中の各注目画素を注目画素値およびその周辺画素値の重み付加により算出される値に変換するものであり、例えば３×３フィルタにおいては、入力データに対して図３に示したごとくの動作を行なう。もちろん、フィルタサイズはもっと大きくてもよい。ハイパスフィルタ２１は、フィルタ係数の和が０である２次元ＦＩＲフィルタであり、入力画像信号３１に対して空間周波数処理を行ない、注目画素に対して補正すべき補正量信号３２を出力する。また、ローパスフィルタ２２は、フィルタ係数の和が１かつ各係数が等しい２次元ＦＩＲフィルタであり、入力画像信号３１に対して空間周波数処理を行ない、ＤＣ成分信号３５を出力する。
【００７４】
ルックアップテーブル２３は、ハイパスフィルタ２１から補正量信号３２を受け取り、ＡＣ成分ゲイン信号３３を生成し出力するルックアップテーブルである。また、ルックアップテーブル２４は、ローパスフィルタ２２からＤＣ成分信号３５を受け取り、ＤＣ成分ゲイン信号３６を生成し出力するルックアップテーブルである。
【００７５】
乗算部２５は、ルックアップテーブル２３から入力されるＡＣ成分ゲイン信号３３およびルックアップテーブル２４から入力されるＤＣ成分ゲイン信号３６により、ハイパスフィルタ２１から入力される補正量信号３２を増幅調整し、調整補正量信号３４として出力する。
【００７６】
加算部２６は、入力画像信号３１と調整補正量信号３４とを加算演算し、空間周波数処理を施した出力画像信号３７を出力する。
【００７７】
上述の構成において、空間周波数処理部１２に入力された入力画像信号３１は、加算部２６とハイパスフィルタ２１とローパスフィルタ２２に入力される。ハイパスフィルタ２１はフィルタ係数の和が０の２次元ＦＩＲフィルタであるので、このハイパスフィルタ２１によって、入力画像信号３１からＡＣ成分のうちの補正量のみが取り出され、補正量信号３２として出力される。そして、ルックアップテーブル２３によって補正量信号３２からＡＣ成分のゲイン調整を行なうためのＡＣ成分ゲイン信号３３が生成され、乗算部２５に出力される。このルックアップテーブル２３で行なわれる演算が、図４（Ｂ）に示す同じ空間周波数で異なるＡＣ成分の入力信号に対応するための補正である。
【００７８】
一方、ローパスフィルタ２２は、入力画像信号３１からＤＣ成分を抽出し、ＤＣ成分信号３５として出力される。そして、ルックアップテーブル２４によってＤＣ成分信号３５からＤＣ成分のゲイン調整を行なうためのＤＣ成分ゲイン信号が生成され、乗算部２５に出力される。このルックアップテーブル２４で行なわれる演算が、図４（Ａ）に示す同じ空間周波数で異なるＤＣ成分の入力信号に対応するための補正である。
【００７９】
ハイパスフィルタ２１で得られた補正量信号３２に対し、ＡＣ成分の相違による調整量であるＡＣ成分ゲイン信号３３と、ＤＣ成分の相違による調整量であるＤＣ成分ゲイン信号３５とを乗算部２５で乗算し、最終的な調整量を表わす調整補正量信号３４を生成する。そして、入力画像信号３１と調整補正量信号３４とを加算部２６で加算し、出力画像信号３７を出力する。
【００８０】
このような構成において、決定すべき処理係数は、ハイパスフィルタ２１の２次元ＦＩＲフィルタのフィルタ係数、および、ルックアップテーブル２３，ルックアップテーブル２４の変換テーブルである。これらの値は出力装置の特性に従って決定されるものであり、以下、これらの値の決定方式について順次説明して行く。なお、図４に示したようなＤＣ成分、ＡＣ成分の相違による空間周波数特性の変化を無視できる出力装置のみが接続されたシステムであれば、ローパスフィルタ２２、ルックアップテーブル２３，２４、乗算部２５は不要であり、ハイパスフィルタ２１の出力を入力画像信号３１に加算部２６で直接加算すれば補正を行なうことができる。この場合、以下の説明中、ハイパスフィルタ２１のフィルタ係数のみを求めればよいことになる。もちろん、ＤＣ成分の相違を考慮しない構成や、ＡＣ成分の相違を考慮しない構成も可能である。
【００８１】
図６は、空間周波数処理部の第１の例において用いる処理係数を決定する処理の一例を示すフローチャートである。図６に示した空間周波数処理で用いる処理係数を決定する処理は、画像データ出力工程４０、出力画像の測定工程５０、空間周波数処理係数の決定工程６０の３工程からなる。まず、画像データ出力工程４０において、特定の空間周波数成分からなるテストパターン画像データを、処理係数を決定したい出力装置で出力し、その出力画像を作成する。次いで、出力画像測定工程５０において、その出力画像の空間周波数特性を測定する。最後に、空間周波数処理係数の決定工程６０において、テストパターン画像データの空間周波数成分と、出力画像測定工程５０で測定したテストパターンの出力画像の空間周波数特性との関係から、出力装置の空間周波数伝達特性を求め、それを補正するための各空間周波数処理係数を算出する。
【００８２】
図６に示す処理係数を決定する処理について具体例を用いながらさらに説明する。以下の説明では、具体的な出力装置の例として、入力データの解像度が１０ｄｏｔ／ｍｍのカラープリンタを想定する。
【００８３】
まず、Ｓ４１において、テストパターン画像データのレイアウトを決定する。図７は、空間周波数処理部の第１の例において用いるテストパターン画像データのレイアウトの一例を示す説明図である。図中、７１は基本ＡＣパッチ画像群、７２は補助ＡＣパッチ画像群、７３はＤＣパッチ画像群である。図７に示したテストパターン画像データのレイアウトは、基本ＡＣパッチ画像群７１、補助ＡＣパッチ画像群７２、ＤＣパッチ画像群７３から構成されている。各パッチ画像群には、複数のパッチ画像が含まれている。図７では、合計で２９個のパッチ画像が示されている。各パッチ画像のサイズは１０ｄｏｔ／ｍｍのプリンタで出力した際に２０ｍｍ四方となるように、２００画素×２００画素としている。
【００８４】
ここで、基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２における各パッチ画像は、０ｌｐ／ｍｍ以外の空間周波数成分を有するパッチ画像、すなわち濃淡の変化のあるパッチ画像である。基本ＡＣパッチ画像群７１は、出力装置の基本的な空間周波数伝達特性を求め、図５に示したハイパスフィルタ２１のフィルタ係数を決定するためのものである。基本ＡＣパッチ画像群７１のパッチ画像には、Ｖ１〜Ｖ８で示される垂直方向の縞のパッチ画像と、Ｈ１〜Ｈ８で示される水平方向の縞のパッチ画像が設けられる。Ｖ１とＨ１，Ｖ２とＨ２，・・・，Ｖ８とＨ８はそれぞれ空間周波数が同じで縞の方向が異なるパッチ画像である。
【００８５】
補助ＡＣパッチ画像群７２は、図５に示したルックアップテーブル２３，２４の変換テーブルの係数を決定するためのものである。基本ＡＣパッチ画像群７１のパッチ画像と同様に、Ｖ２．１〜Ｖ２．４で示されるパッチ画像は垂直方向の縞であり、Ｈ２．１〜Ｈ２．４で示されるパッチ画像は水平方向の縞である。ここでは基本ＡＣパッチ画像群７１中のＶ２，Ｈ２のパッチを基本として、図４（Ａ）に示すようにＤＣ成分の異なるパッチ画像をＶ２．１，Ｖ２．２，Ｈ２．１，Ｈ２．２にレイアウトする。また、図４（Ｂ）に示すようにＡＣ成分の異なるパッチ画像としてＶ２．３，Ｖ２．４．Ｈ２．３，Ｈ２．４をレイアウトしている。なお、Ｖ２．１とＨ２．１，・・・，Ｖ２．４とＨ２．４はそれぞれ方向が異なるのみでＡＣ成分、ＤＣ成分とも同じパッチ画像である。
【００８６】
ＤＣパッチ画像群７３における各パッチ画像は、０ｌｐ／ｍｍ以外の空間周波数成分を含まないパッチ画像、すなわち濃度変化のない平坦な濃度のパッチ画像であり、後述するように、基本ＡＣパッチ画像群７１、補助ＡＣパッチ画像群７２の測定値を規格化するためのものである。Ｄ０〜Ｄ２５５はそれぞれ濃度の異なるパッチ画像である。
【００８７】
図７に示したテストパターン画像データのレイアウトは一例であって、どの位置にどのようなパッチ画像が配置されているかがわかれば他のレイアウトであってもかまわない。また、それぞれのパッチ画像群が別々の複数枚のテストパターン画像として形成されるものであってもよい。各パッチ画像群に含まれるパッチ画像の数や特性、大きさ等は任意である。特にこの例では補助ＡＣパッチ群７２は基本ＡＣパッチ群７１中のパッチ画像Ｖ２，Ｈ２をもとに作成しているが、他の基本ＡＣパッチ群７１中のパッチ画像をもとにしたパッチ画像を追加することによって精度を向上させることが可能である。
【００８８】
なお、例えば、図７に示したようなテストパターン画像データのレイアウトは、予め作成してテストパターン保持部５などに記憶させておいてもよい。複数のテストパターン画像データのレイアウトを記憶させておき、選択して用いるように構成することもできる。
【００８９】
Ｓ４２において、例えば、図７に示したようなテストパターン画像データのうちの基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２のパッチ画像の画像データを作成する。上述のように、図７において、基本ＡＣパッチ画像群７１の各パッチに付された記号Ｖ１〜Ｖ８は垂直方向のライン、Ｈ１〜Ｈ８は水平方向のラインを表わしており、記号ＶまたはＨの後ろの数字は各ラインの幅を画素数で表わしている。図８は、空間周波数処理部の第１の例において用いる基本ＡＣパッチ画像群のパッチ画像の一例の説明図である。図中、Ｂは黒画素を表わし、Ｗは白画素を表わしている。図８（Ａ）はＶ２のパッチ画像の一部を示しており、幅が２画素の白と黒の垂直方向のラインが交互に配列されている。また、図８（Ｂ）はＨ３のパッチ画像の一部を示しており、幅が３画素の白と黒の水平方向のラインが交互に配列されている。ここでは、黒画素のデータ値を１６０、白画素のデータ値を９６としている。
【００９０】
次に、図７に示したテストパターン画像データのうちの補助ＡＣパッチ画像群７２について、各パッチ画像Ｖ２．１〜Ｖ２．４，Ｈ２．１〜Ｈ２．４の画像データを作成する。上述のように、基本ＡＣパッチ画像群７１中のＶ２，Ｈ２のパッチを基本として、ＤＣ成分の異なるパッチ画像Ｖ２．１，Ｖ２．２，Ｈ２．１，Ｈ２．２と、ＡＣ成分の異なるパッチ画像Ｖ２．３，Ｖ２．４．Ｈ２．３，Ｈ２．４を作成する。
【００９１】
図９は、空間周波数処理部の第１の例において用いる基本ＡＣパッチ画像群と補助ＡＣパッチ画像群の各パッチ画像の特性値の一例の説明図である。図９には、基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２に含まれるパッチ画像について、その空間周波数、黒画素のデータ値、白画素のデータ値、ＤＣ成分値、ＡＣ成分値、信号コントラスト値を示している。空間周波数は、１ｍｍあたり何本のラインが含まれるかを示したものであり、Ｖ１，Ｈ１であれば１画素おきに白と黒のラインが交互に現われるため、想定しているプリンタの解像度である１ｍｍあたり１０ｄｏｔ中に５本のラインが含まれることを示している。ＤＣ成分値は平均濃度を示しており、ＡＣ成分値は黒画素のデータ値および白がそのデータ値とＤＣ成分値の差分を示している。なお、濃度は、最も黒の濃度が濃い場合を２５５、最も黒の濃度が薄い白部分の場合を０としている。信号コントラスト値は後述する。
【００９２】
図９にも示したように、基本ＡＣパッチ画像群７１の各パッチ画像では、黒画素のデータ値、白画素のデータ値、ＤＣ成分値、ＡＣ成分値はすべて同じであり、空間周波数のみが異なる。また、補助ＡＣパッチ画像群７２の各パッチ画像では、図９に示すように空間周波数が一定で、黒画素のデータ値、白画素のデータ値としてそれぞれの値を有する。
【００９３】
例えば、図９に示した各パラメータに従って、各パッチ画像の画像データを作成し、作成した各パッチ画像データを例えば図７に示すようなテストパターン画像データのレイアウトに従って埋め込む。
【００９４】
次に、Ｓ４３において、基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２の各パッチ画像の画像データから特性値を算出する。基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２における各パッチ画像の各画素のデータ値は、例えば基本ＡＣパッチ画像群７１においては、黒画素値＝１６０，白画素値＝９６というように設定した。これは、別の表現をすると、ＤＣ成分値＝１２８，ＡＣ成分値＝±３２と表わすことができる。ＤＣ成分値，ＡＣ成分値で表現した各パッチの設定値は図９に示した。
【００９５】
ここで、基本ＡＣパッチ画像群７１および補助ＡＣパッチ画像群７２の各パッチ画像の画像データの特性値として、ＤＣ成分値と信号のコントラスト値を考える。ＤＣ成分値は、図９に示した各ＤＣ成分値である。信号のコントラスト値は、上述した各ＡＣ成分値をもとに補正を行なう。画像の場合には信号は正弦波というよりも矩形波であるため、実効コントラスト値は正弦波のときよりも向上する。このことを考慮してコントラスト値の補正を行なう。実際の補正は、
信号コントラスト値＝ＡＣ成分値×（４／３．１４）
によって計算することができる。このような補正を行なって求めた各パッチ画像のコントラスト値を図９の信号コントラスト値の欄に示している。
【００９６】
Ｓ４４において、ＤＣパッチ画像群７３のパッチ画像の画像データを作成する。図７において、ＤＣパッチ画像群７２の各パッチ画像に付された記号Ｄ０〜Ｄ２５５において、記号ＤはＤＣ成分すなわち０ｌｐ／ｍｍのみの空間周波数成分からなるパッチ画像であることを示し、後ろの数字はＤＣ成分のデータ値を表わしている。たとえば、Ｄ６４のパッチ画像は、２００画素×２００画素の全画素のデータ値が６４である均一なパッチ画像を意味する。
【００９７】
ここでは、データ値＝０からデータ値＝２５５までを、ほぼ均等に分割し、９個のパッチ画像を用意している。各パッチ画像の画像データを作成し、作成した各パッチ画像の画像データを、例えば図７に示したようなテストパターン画像データのレイアウトに従って埋め込む。これにより、テストパターン画像データが完成する。完成したテストパターン画像データは、テストパターン保持部５に保持される。同一のテストパターン画像データが、Ｙ，Ｍ，Ｃ用に３色分用意される。あるいは１つのテストパターン画像データが３色共通に使用される。
【００９８】
画像データ出力工程４０の最後として、Ｓ４５においてテストパターン画像データを出力装置から出力させ、テストパターン画像を作成する。テストパターン画像の出力に際しては、まず出力装置を単色モードに設定し、テストパターン画像の出力を行なう。例えば、Ｙ単色モードに設定されているとき、テストパターン保持部５からはＹ用画像信号としてテストパターン画像データが出力される。テストパターン保持部５から出力されたテストパターン画像データは、図２に示す階調変換部１３においてルックアップテーブルを用いて階調変換処理が施された後、出力装置から出力される。出力装置では、Ｙ用の階調変換処理が施されたテストパターン画像データを受け取り、Ｙ単色モードにてその出力の画像を出力する。Ｍ，Ｃに対しても同様に、出力装置をＭ単色モード、Ｃ単色モードに設定しながら、それぞれのテストパターン画像を得、計３枚のテストパターン画像を得る。
【００９９】
得られたテストパターン画像は例えば図７と同様のレイアウトの画像である。ただし、各パッチ画像は、出力装置の特性を反映しており、テストパターン画像データを忠実に反映した画像とは限らない。以下、得られた３枚のテストパターン画像のそれぞれに対して、出力画像の測定工程５０および空間周波数処理係数の決定工程６０を行なう。説明の都合上、そのうちの１色についてのみ説明するが、３色それぞれについて同様の処理を行なう。
【０１００】
出力画像の測定工程５０では、テストパターン画像データを出力装置から出力したテストパターン画像の測定を行ない、測定結果を解析する。まずＳ５１において、テストパターン画像中の基本ＡＣパッチ画像群および補助ＡＣパッチ画像群について、微少領域ごとに反射率を測定する。反射率の測定には、例えばミクロ反射計などを用いることができる。図１０は、反射率の測定の具体例の説明図である。図中、８１はアパーチャ、８２はパッチ画像である。反射率の測定の具体例として、図１０に示すように、例えば、０．１ｍｍ×０．０２ｍｍのアパーチャ８１を使用し、アパーチャの長辺がパッチ画像８２中のライン構造と平行になるようにセットする。そして、アパーチャの短辺方向（図１０中の矢印方向）に測定位置を０．０２ｍｍづつ移動しながら各位置における反射光量を測定し、測定結果を得る。この際、補色フィルタを使用して測定するとよい。
【０１０１】
続いてＳ５２において、テストパターン画像中のＤＣパッチ画像群の各パッチ画像についても微少領域ごとに反射率を測定し、平均反射率を求める。この場合も、例えばミクロ反射計を用いることができる。ＤＣパッチ画像群は上述のようにＤＣ成分のみの均一なパッチ画像なので、アパーチャを任意な方向にセットし、アパーチャの短辺方向に測定位置を移動しながら各位置における反射光量を各パッチに対して複数点で測定し、測定結果の平均値を得る。この際、基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の測定に際してフィルタを使用した場合には、ここでも同一のフィルタを使用する。
【０１０２】
このようにして得られた測定値に対して解析を行なう。まずＳ５３において、基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の各パッチ画像における測定値を規格化する。図１１は、ＤＣパッチ画像群の各パッチ画像の平均反射率のデータとテストパターン画像データにおけるＤＣパッチ画像群の各パッチ画像のデータ値との対応関係の一例を示すグラフである。Ｓ５２で測定して得られたＤＣパッチ画像群の各パッチ画像の平均反射率のデータと、テストパターン画像データにおけるＤＣパッチ画像群の各パッチ画像のデータ値とを対応づけると、例えば図１１に示すような関係が得られる。まずこの対応関係ｆdc（ｘ）を近似値により求める。基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の各パッチ画像に対するミクロ反射計での測定値を、この対応関係ｆdc（ｘ）を用いて変換する。これにより、測定結果を処理演算部１で取り扱うＹＭＣ画像信号と同等に扱うことが可能となる。
【０１０３】
図１２は、ＡＣパッチ画像の測定値の規格化および成形処理による波形の変化の一例の説明図である。基本ＡＣパッチ画像群または補助ＡＣパッチ画像群のあるパッチ画像に対する測定データ列を図１２（Ａ）に示している。この測定データ列を例えば図１１に示すように得られた対応関係ｆdc（ｘ）を用いて出力装置に依存した画像信号値に変換し、規格化データ列を作成する。作成された規格化データ列を図１２（Ｂ）に示す。この規格化データ列は、図２において色空間変換部１１から出力されるＹＭＣ信号と同等のものである。ここではＹＭＣ信号は８Ｂｉｔ２５６階調としているが、これに限られるものではない。なお、上述のように、規格化データ列において値が小さい方が白に近く、値が大きい方が色濃度が濃くなることを示している。そのため、反射率の値の大小とは逆になっている。
【０１０４】
次にＳ５４において、Ｓ５３で規格化した測定値に対して整形処理を行ない、成形データ列を作成する。図１３は、規格化データ列に対する成形処理の一例を示すフローチャートである。例えば、図１２（Ｂ）に示すような波形の規格化データ列は、上述のように２０ｍｍ四方のパッチ画像を０．０２ｍｍ間隔で測定した場合には１０００個の測定データからなっている。Ｓ９１において、これら１０００個のデータからなる規格化データ列の中央部から５１２個のデータを切り出す。次にＳ９２において、Ｓ９１で切り出した５１２個のデータを離散フーリエ変換（ＤＦＴ）処理する。これにより、空間周波数成分が算出される。Ｓ９３において、Ｓ９２で算出された空間周波数成分のうち、不要な高周波成分を除去し、その後、Ｓ９４において離散フーリエ逆変換処理する。このような成形処理によって、ノイズなどによる高調波成分が除去されて、例えば、図１２（Ｃ）にしめすような、Ｓｉｎ波形を基本とした整形データ列が得られる。
【０１０５】
Ｓ５５において、基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の各パッチ画像に対応する整形データ列をもとに、それらの特性値としてコントラスト値を算出する。図１４は、空間周波数処理部の第１の例における整形データ列に対する特性値の測定手法の一例の説明図である。図１４に示すように、整形データ列に対して、破線で示したようなウインドウ領域を設定し、このウインドウ領域内での最大値と最小値を検出し、（最大値−最小値）を算出する。ウインドウ領域を移動しながら、（最大値−最小値）を繰り返し算出し、その平均値をコントラストとする。また、最大値と最大値の間、あるいは最小値と最小値の間など、特徴点の間のデータ数を係数することによって周期を求めることができる。測定ピッチは０．２ｍｍであったから、０．２×データ数で周期が求められる。また、周期の逆数を計算することにより、空間周波数が求められる。ここで、ウインドウ領域のサイズは、基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の各パッチ画像に含まれる縞状の線の２〜３本分程度に設定すればよい。
【０１０６】
図１５は、空間周波数処理部の第１の例においてテストパターン画像から得られた基本ＡＣパッチ画像群の各パッチ画像の周期および空間周波数の一例の説明図、図１６は、同じく得られたコントラスト値の一例を示すグラフである。図１５には、上述のようにして算出された基本ＡＣパッチ画像群の各パッチ画像に含まれる縞状の線の周期を長さとデータ数で示すとともに、算出される空間周波数と、使用したウインドウ領域のサイズをデータ数で示している。また、各パッチ画像ごとに算出された空間周波数とコントラストの関係の一例を図１６に示している。図１６において、黒丸はＨ１〜Ｈ８の各パッチ画像における空間周波数とコントラストの関係を示し、白丸はＶ１〜Ｖ８の各パッチ画像における空間周波数とコントラストの関係を示している。
【０１０７】
なお、図１５、図１６には基本ＡＣパッチ画像群の各パッチ画像について、周期、空間周波数、コントラスト等を示したが、補助ＡＣパッチ画像群の各パッチ画像についても同様に周期を求めて空間周波数を算出し、また、コントラストを算出する。
【０１０８】
空間周波数処理係数の決定工程６０では、テストパターン画像データの基本ＡＣパッチ画像群の各パッチ画像のデータのコントラスト値と、それに対応するテストパターン画像のコントラスト値との関係から出力装置の空間周波数伝達特性を求め、それを補正するための空間周波数処理部１２における処理係数を算出する。
【０１０９】
Ｓ６１において、出力装置の空間周波数伝達特性を算出する。空間周波数伝達特性として各空間周波数のコントラスト伝達率を考える。出力装置の空間周波数伝達率は、（出力画像のコントラスト値／入力データのコントラスト値）により求まる。ここで、基本ＡＣパッチ画像群のうち、記号Ｖで示した垂直方向の縞からなるパッチ画像からは、それに直交する水平方向の空間周波数伝達率が求まり、同様に、記号Ｈで示した水平方向の縞からなるパッチ画像からからは、それに直交する垂直方向の空間周波数伝達率が求まる。なお、空間周波数伝達率の算出は、基本ＡＣパッチ画像群および補助ＡＣパッチ画像群の両者に関して行なう。
【０１１０】
図１７は、基本ＡＣパッチ画像群における算出された空間周波数伝達特性の一例を示すグラフである。この例において、出力画像のコントラスト値は、上述の図１６に示すグラフの値であり、入力データのコントラスト値は図９に示す信号コントラスト値である。ただし、図９に示した信号コントラスト値はＤＣ成分からの差分で示しているので、示されている値を２倍して（最大値−最小値）の値に直すか、あるいは出力画像のコントラスト値を１／２にしてから計算する必要がある。
【０１１１】
図１７において、黒丸は図１６で黒丸で示した水平方向の縞のパッチ画像におけるコントラストから得られた垂直方向の空間周波数伝達率を表わしており、白丸は図１６で白丸で示した垂直方向の縞のパッチ画像におけるコントラストから得られた水平方向の空間周波数伝達率を表わしている。
【０１１２】
次にＳ６２において、基本的な空間周波数伝達特性を決定する。ここで、基本的な空間周波数伝達特性とは、図４で説明したようなＤＣ成分の違いやＡＣ成分の違いを考慮しない空間周波数伝達特性のことを指す。
【０１１３】
図２に示したように、処理演算部１に入力された画像データは、空間周波数処理部１２にて空間周波数伝達率Ｆ１の空間周波数処理を施された後に、出力部２に送られ、空間周波数伝達率Ｆ２である画像形成工程を経て、出力画像が作成される。したがって、入力された画像データの有する空間周波数成分をＳ１，出力画像の有する空間周波数成分をＳ２とすると、Ｓ２＝Ｓ１×Ｆ１×Ｆ２となる。
【０１１４】
上述のごとく、出力装置によらず同一の鮮鋭度を得るためには、空間周波数処理部１２における空間周波数伝達率Ｆ１は、画像形成工程の空間周波数伝達率Ｆ２を補正するものであればよい。すなわち、空間周波数処理部１２における空間周波数伝達率Ｆ１＝１／Ｆ２となるように、空間周波数処理部１２に設定する処理係数を設計すればよい。
【０１１５】
図１８は、基本ＡＣパッチ画像群の各パッチ画像から求められた空間周波数伝達特性に基づいて設計された空間周波数処理部の第１の例における空間周波数伝達率の一例を示すグラフである。図１８において、破線は図１７に示した出力部２の空間周波数伝達率の逆数のグラフであり、実線は実際に設計した空間周波数処理部の空間周波数伝達率のグラフである。例えば、出力部２の基本ＡＣパッチ画像群の各パッチ画像から求まる空間周波数伝達特性が図１７に示すものであるとき、空間周波数処理部１２の空間周波数伝達率の目標値は、図１８に破線で示すものとなる。このままでは高い空間周波数帯域においてノイズが極端に強調されてしまう。ノイズ成分の強調を抑えるため、解像度１０ｄｏｔ／ｍｍのカラープリンタでは、そのナイキスト周波数である５ｌｐ／ｍｍにおいて伝達率が１になるように設定している。
【０１１６】
このようにして設計された空間周波数処理部１２の空間周波数伝達率を実現するように、Ｓ６３において、基本的な空間周波数処理における処理係数を算出する。空間周波数処理部１２として例えば図５に示す構成を用いる場合、そのうち基本的な空間周波数処理はハイパスフィルタ２１が行なう。このＳ６３における処理では、このハイパスフィルタ２１のフィルタ係数を決定する。
【０１１７】
まず、先に基本ＡＣパッチ画像群より求めた、垂直、水平の２方向の空間周波数伝達率の設計値をもとに、設計値を２次元面補間し、その逆フーリエ変換を施す。次に逆フーリエ変換をした結果を総和が１となるように規格化することで、目標とする空間周波数伝達特性のフィルタ係数が算出される。このようなフィルタ係数の算出方法は種々開発されており、適宜用いればよい。ただし、ハイパスフィルタ２１は、上述したとおり補正信号のみを出力するフィルタである。そのため、算出されたフィルタ係数のうち、中心画素に対する係数から１を減算し、ハイパスフィルタ２１のフィルタ係数として決定する。
【０１１８】
次にＳ６４において、例えば、図４（Ｂ）に示したような、ＡＣ成分による非線形性の補正係数を算出する。上述したように、ハイパスフィルタ２１は画像信号のコントラストが±４０．５である基本ＡＣパッチ画像群を基にフィルタ係数を算出した。いま、水平方向の縞のパッチ画像における空間周波数が２．５ｌｐ／ｍｍの入力画像データに対するハイパスフィルタ２１の応答を考える。
【０１１９】
図９に示すように、テストパターン画像データにおける記号Ｖ２のパッチ画像は、ＤＣ成分＝１２８、信号コントラスト＝±４０．５、空間周波数＝２．５ｌｐ／ｍｍである。このパッチ画像で算出した出力装置の空間周波数伝達率をＴ０とすると、信号コントラスト＝Ｃ１である入力信号に対してハイパスフィルタ２１から出力される補正量信号３２のコントラストＣｈ０は、Ｃｈ０＝Ｃ１×（１／Ｔ０−１）となる。
【０１２０】
ここで、信号コントラスト＝Ｃ１である入力信号に対する出力装置の空間周波数伝達率をＴ１とすると、本来望まれる補正後の信号のコントラストＣｗ１は、Ｃｗ１＝Ｃ１×（１／Ｔ１−１）となる。したがって、入力信号の信号コントラストがＣ１である場合には、ハイパスフィルタ２１から出力される補正量信号３２に対して増幅処理をする必要があり、その増幅率Ｇａｉｎ１は、
Ｇａｉｎ１＝Ｃｗ１／Ｃｈ０＝（１／Ｔ１−１）／（１／Ｔ０−１）
となる。
【０１２１】
ルックアップテーブル２３は、この増幅率Ｇａｉｎ１を得るためのものであり、補正量信号３２によりルックアップテーブルを参照して、対応する増幅率Ｇａｉｎ１を出力するものである。ルックアップテーブル２３の処理係数は、基本ＡＣパッチ画像群の中のＶ２、Ｈ２パッチ画像、補助ＡＣパッチ画像群中のＶ２．３，Ｈ２．３，Ｖ２．４，Ｈ２．４パッチ画像に対応する出力装置の空間周波数伝達率を算出し、その算出結果を用いて決定する。
【０１２２】
具体例として、信号コントラストＣ０が±４０．５であるＶ２パッチ画像に対応する空間周波数伝達率Ｔ０が０．８であり、信号コントラストＣ１が±２０．５であるＶ２．３パッチ画像に対応する空間周波数伝達率Ｔ１が０．９、信号コントラストＣ２が±８１であるＶ２．４パッチ画像に対応する空間周波数伝達率Ｔ２が０．７であるとする。このとき、信号コントラスト＝Ｃ０，Ｃ１，Ｃ２である入力信号に対してハイパスフィルタ２１から出力される補正量信号３２のコントラストをＣｈ０，Ｃｈ１，Ｃｈ２とすると、
Ｃｈ０＝Ｃ０×（１／Ｔ０−１）＝±１０．１
Ｃｈ１＝Ｃ１×（１／Ｔ０−１）＝±５．１
Ｃｈ２＝Ｃ２×（１／Ｔ０−１）＝±２０．３
となる。またルックアップテーブル２３における各々に対する増幅率をＧ０，Ｇ１，Ｇ２とすると、
Ｇ０＝（１／Ｔ０−１）／（１／Ｔ０−１）＝１
Ｇ１＝（１／Ｔ１−１）／（１／Ｔ０−１）＝０．４４
Ｇ２＝（１／Ｔ２−１）／（１／Ｔ０−１）＝１．７１
となる。
【０１２３】
図１９は、補正量信号のコントラストと増幅率の関係を示すグラフである。上述のようにしてＣｈ０とＧ０、Ｃｈ１とＧ１、Ｃｈ２とＧ２の組の値が求められた。これらをグラフ上にプロットし、これらの点を補間することにより、図１９に示すようなグラフが得られる。このグラフと同等の変換テーブルを作成する。同様に、Ｈ２、Ｈ２．３、Ｈ２．４パッチ画像を用いて変換テーブルを作成し、両者の平均をルックアップテーブル２３の係数とすればよい。
【０１２４】
次にＳ６５において、例えば、図４（Ａ）に示したような、ＤＣ成分による非線形性の補正係数を算出する。ルックアップテーブル２３の処理係数算出時と同様に、ハイパスフィルタ２１の応答を考える。具体例として、Ｖ２パッチ画像（ＤＣ成分＝１２３，データコントラスト＝±４０．５、空間周波数＝２．５ｌｐ／ｍｍ）で算出した出力装置の空間周波数伝達率をＴ０とすると、データコントラスト＝Ｃ０である入力信号に対してハイパスフィルタ２１から出力される補正量信号３２のコントラストＣｈ０は、
Ｃｈ０＝Ｃ１×（１／Ｔ０−１）
となる。
【０１２５】
ここで、ＤＣ成分＝ｄ１である入力信号に対する出力装置の空間周波数伝達率をＴ’１とすると、本来望まれる補正後の信号のコントラストＣｗ’１は、
Ｃｗ’１＝Ｃ０×（１／Ｔ’１−１）
となる。したがって、入力信号のＤＣ成分＝ｄ１である場合には、ハイパスフィルタ２１から出力される補正量信号３２に対して増幅処理をする必要があり、その増幅率Ｇａｉｎ２は、
Ｇａｉｎ２＝Ｃｗ’１／Ｃｈ０＝（１／Ｔ’１−１）／（１／Ｔ０−１）
となる。
【０１２６】
ルックアップテーブル２４は、この増幅率Ｇａｉｎ２を得るためのものであり、入力であるＤＣ成分信号３５によりルックアップテーブルを参照して対応する増幅率Ｇａｉｎ２を出力するものである。ルックアップテーブル２４の処理係数は、基本ＡＣパッチ画像群７１中のＶ２，Ｈ２パッチ画像、補助ＡＣパッチ画像群７３中のＶ２．１，Ｈ２．１，Ｖ２．２，Ｈ２．２パッチ画像に対応する出力装置の空間周波数伝達率の算出結果を用いて決定する。
【０１２７】
具体例として、ＤＣ成分ｄ０が１２８であるＶ２パッチ画像に対応する空間周波数伝達率Ｔ０が０．８であり、ＤＣ成分ｄ１が６４であるＶ２．１パッチ画像に対応する空間周波数伝達率Ｔ１が０．７であり、ＤＣ成分ｄ２が１９２であるＶ２．２パッチ画像に対応する空間周波数伝達率Ｔ２が０．９であるとする。ＤＣ成分＝ｄ０，ｄ１，ｄ２である入力信号に対するルックアップテーブル２４における各々に対応する増幅率をＧ’０，Ｇ’１，Ｇ’２，とすると、
Ｇ’０＝（１／Ｔ０−１）／（１／Ｔ０−１）＝１
Ｇ’１＝（１／Ｔ’１−１）／（１／Ｔ０−１）＝１．７１
Ｇ’２＝（１／Ｔ’２−１）／（１／Ｔ０−１）＝０．４４
となる。
【０１２８】
図２０は、ＤＣ成分信号値と増幅率の関係を示すグラフである。上述のようにしてｄ０とＧ’０、ｄ１とＧ’１、ｄ２とＧ’２の組の値が求められた。これらをグラフ上にプロットし、これらの点を補間することにより、図２０に示すようなグラフが得られる。このグラフと同等の変換テーブルを作成する。同様に、Ｈ２、Ｈ２．１、Ｈ２．２パッチ画像を用いて変換テーブルを作成し、両者の平均をルックアップテーブル２４の係数とすればよい。
【０１２９】
以上により、ハイパスフィルタ２１のフィルタ係数、ルックアップテーブル２３および２４の係数が得られた。これらの係数をそれぞれ設定し、空間周波数処理部１２で空間周波数処理を行なう。空間周波数処理部１２では、出力部２に応じて設定された処理係数に従って空間周波数処理を行なうことにより、出力部２の特性によらない空間周波数特性を有する出力画像を得ることができる。
【０１３０】
上述のテストパターン画像データの基本ＡＣパッチ画像群および補助ＡＣパッチ画像群中の各パッチ画像に含まれる空間周波数成分は、各々単一空間周波数成分とその高調波成分であった。しかしこれに限らず、各パッチ画像に含まれる空間周波数成分は各々複数の空間周波数成分であってもよい。図２１は、ＡＣパッチ画像群のパッチ画像の別の例の説明図である。図２１（Ａ）に示すパッチ画像の例は、黒画素＝１６０，白画素＝９６が横方向に各々２画素づつ繰り返される縦縞のパターンである。また、図２１（Ｂ）に示すパッチ画像の例は、黒画素＝１６０，白画素＝９６が横方向に各々３画素づつ繰り返される縦縞のパターンである。例えば、この図２１（Ａ）に示すパッチ画像のパターンと図２１（Ｂ）に示すパッチ画像のパターンとを合成し、図２１（Ｃ）に示すようなパターンをテストパターン画像データのパッチ画像として用いることもできる。この場合、図６のＳ５４において説明したＡＣパッチ規格化測定値の整形処理において、図１３のＳ９３で不要な高周波成分を除去する際に、空間周波数データをＤＣ成分と合成前の各々の画像データパターンに対応した空間周波数周波数成分に分解し、各々に対してＳ９４で離散フーリエ逆変換を行なうことにより、同様に空間周波数処理のための処理係数を算出することができる。
【０１３１】
また、上述の例では補助ＡＣパッチ画像としてＶ２，Ｈ２パッチ画像と同じ空間周波数のパッチ画像のみを用意したが、これに限らず、種々の空間周波数に対応する補助ＡＣパッチ画像を用意してもよい。このとき、図１９、図２０に示したような変換テーブルは空間周波数ごとに得られることになる。ルックアップテーブル２３を参照する際には、それぞれ、補正量信号３２のコントラストとともに空間周波数を用いて、また、ルックアップテーブル２４を参照する際にはＤＣ成分信号３５とともにＡＣ成分の空間周波数を用いて、それぞれ参照するように構成すればよい。この場合、空間周波数の違いによるＡＣ成分、ＤＣ成分による補正量の違いを補正することが可能となる。
【０１３２】
図２２は、空間周波数処理部の構成の第２の例を示すブロック図である。図中、１１１は基本補正部、１１２は補助補正部、１１３は増幅率決定部、１１４は乗算部、１１５は加算部、１２１は１２１は入力画像信号、１２２は基本補正信号、１２３は補助補正基本信号、１２４は補助補正量増幅率信号、１２５は補助補正信号、１２６は出力画像信号である。
【０１３３】
基本補正部１１１および補助補正部１１２は、入力画像信号１２１の空間周波数特性を変更するためのものである。基本補正部１１１および補助補正部１１２としては、例えば２次元ＦＩＲフィルタなどのデジタルフィルタを用いることができる。デジタルフィルタを用いることにより、再現性の高い処理系の設計が容易になる。また、既存の専用チップを用いることにより高速な処理を安価で実現できる。また、ＦＩＲフィルタにすることにより、入力画像信号１２１に対して処理系を容易に安定化することができる。基本補正部１１１は、例えばフィルタ係数の和が１である２次元フィルタとすることができ、入力画像信号１２１に対して空間周波数処理を行ない、基本補正信号１２２を出力する。補助補正部１１２は、例えばフィルタ係数の和が０である２次元フィルタとすることができ、入力画像信号１２１に対して空間周波数処理を行ない、補助補正基本信号１２３を出力する。
【０１３４】
増幅率決定部１１３は、入力画像信号１２１から補助補正基本信号１２３の増幅率を決定し、補助補正量増幅率信号１２４として出力する。詳しい内部動作については後述する。
【０１３５】
乗算部１１４は、増幅率決定部１１３より入力される補助補正量増幅率信号１２４により、補助補正部１１２から入力される補助補正基本信号１２３を増幅調整し、補助補正信号１２５として出力する。
【０１３６】
加算部１１５は、基本補正信号１２２と補助補正信号１２５とを加算演算し、空間周波数処理を施した出力画像信号１２６を出力する。
【０１３７】
上述の構成において、空間周波数処理部１２に入力された入力画像信号１２１は、基本補正部１１１と補助補正部１１２と増幅率決定部１１３に入力される。基本補正部１１１は、入力画像信号１２１に対して基本的な補正を施し、基本補正信号１２２として加算部１１５に出力する。
【０１３８】
一方、補助補正部１１２は、入力画像信号１２１より基本補正部１１１を用いた補正の過不足分を調整するための予備信号である補助補正信号基本信号１２３を算出し、乗算部１１４に出力する。乗算部１１４で補助補正基本信号１２３は、入力画像信号１２１より増幅率決定部１１３において算出された補助補正量増幅率信号１２４に基づき増幅調整され、補助補正信号１２５として加算部１１５に出力される。そして、基本補正信号１２２と補助補正信号１２５を加算部１１５で加算し、出力画像信号１２６を出力する。
【０１３９】
この補助補正部１１２と増幅率決定部１１３と乗算部１１４で行なわれる演算が、上述の図４で説明した同じ空間周波数で異なるＤＣ成分および／またはＡＣ成分の入力に対応するための補正である。
【０１４０】
図２３は、増幅率決定部の構成の一例を示すブロック図である。図中、１３１はＡＣ成分抽出器、１３２はＤＣ成分抽出器、１３３は２次元ルックアップテーブル、１３４はＡＣ成分信号、１３５はＤＣ成分信号である。ＡＣ成分抽出器１３１は、注目画素周辺の局所的なＡＣ成分を抽出するためのもので、入力画像信号１２１に対して、例えば注目画素を中心に１１×１１の領域の最大値と最小値の差を算出し、ＡＣ成分信号１３４として出力する。ＡＣ成分抽出器１３１を設けることにより、異なるＡＣ成分を持つ入力画像信号１２１に対して効率よく対応することができる。
【０１４１】
ＤＣ成分抽出器１３２は、注目画素周辺の局所的なＤＣ成分を抽出するためのもので、入力画像信号１２１に対して、例えば注目画素を中心に１１×１１の領域を平均値を算出し、ＤＣ成分信号１３５として出力する。ＤＣ成分抽出器１３２を設けることにより、異なるＤＣ成分を持つ入力画像１２１に対して効率よく対応することができる。
【０１４２】
２次元ルックアップテーブル１３３は、ＡＣ成分抽出器１３１から入力されるＡＣ成分信号１３４と、ＤＣ成分抽出器１３２から入力されるＤＣ成分信号１３５を受け取り、補助補正量増幅率信号１２４を生成し出力する。２次元ルックアップテーブル１３３を用いることにより、入力信号に応じた任意の出力を得ることができる。また、２次元のルックアップテーブルとすることにより、入力のＡＣ成分信号とＤＣ成分信号の任意の組合せに柔軟に応じた出力信号を得ることができる。
【０１４３】
このような構成において、決定すべき処理係数は、基本補正部１１１と補助補正部１１２の例えば２次元ＦＩＲフィルタのフィルタ係数、および、２次元ルックアップテーブル１３３の変換テーブルである。これらの値は出力装置の特性に従って決定されるものであり、以下、これらの値の決定方式について順次説明していく。なお、図４に示したようなＤＣ成分、ＡＣ成分の相違による空間周波数特性の変化を無視できる出力装置のみが接続されたシステムであれば、補助補正部１１２、増幅率決定部１１３、乗算部１１４、加算部１１５は不要であり、基本補正部１１１の出力を直接出力画像信号とすることができる。またこの場合は、補助補正部１１２のフィルタ係数と増幅率決定部１１３内の２次元ルックアップテーブル１３３の変換テーブルのいずれか、または両方の係数を全て０にしてもよい。この場合、以下の説明中、基本補正部１１１のフィルタ係数のみを求めればよいことになる。もちろん、ＡＣ成分の相違を考慮しない構成や、ＡＣ成分の相違を考慮しない構成も可能である。
【０１４４】
図２４は、空間周波数処理部の第２の例において用いる処理係数を決定する処理の一例を示すフローチャートである。図２４に示すフローチャートの概略は上述の図６に示す流れとほぼ同様である。図２４における画像データ出力工程１４０は図６の画像データ出力工程４０に対応し、同様に出力画像の測定工程１５０は出力画像の測定工程５０に対応し、空間周波数処理係数の決定工程１６０は空間周波数処理係数の決定工程６０に対応する。図６に示す空間周波数処理部１２の処理と重複する部分もあるが、以下に空間周波数処理部１２の第２の例における処理係数の決定処理の流れを順に説明して行く。なお、以下の説明では、具体的な出力装置の例として、入力データの解像度が１０ｄｏｔ／ｍｍの場合を想定する。
【０１４５】
まず、Ｓ１４１において、テストパターン画像データのレイアウトを決定する。図２５は、空間周波数処理部の第２の例において用いるテストパターン画像データのレイアウトの一例を示す説明図である。図中、１７１はＡＣパッチ画像群、１７２はＤＣパッチ画像群である。図２５に示したテストパターン画像のレイアウトは、ＡＣパッチ画像群１７１、ＤＣパッチ画像群１７２から構成されている。各パッチ画像群には、複数のパッチ画像群が含まれている。図２５では、合計で８１個のパッチ画像が示されている。各パッチ画像のサイズは１０ｄｏｔ／ｍｍの解像度のプリンタで出力した際に２０ｍｍ四方となるように、２００画像×２００画素としている。
【０１４６】
ここで、ＡＣパッチ画像群１７１における各パッチ画像は、０ｌｐ／ｍｍ以外の空間周波数成分を有するパッチ画像、すなわち濃淡の変化のあるパッチ画像である。ＡＣパッチ画像群１７１は、出力装置の空間周波数伝達特性を求めるためのものである。ＡＣパッチ画像群１７１を用いて求めた出力装置の空間周波数伝達特性に基づき、図２２に示した基本補正部１１１および補助補正部１１２のフィルタ係数と図２３に示した２次元ルックアップテーブル１３３の変換テーブルを決定する。ＡＣパッチ画像群１７１のパッチ画像には、Ｖ１．１〜Ｖ９．８で示される垂直方向の縞のパッチ画像と、Ｈ１．１〜Ｈ９．８で示される水平方向の縞のパッチ画像が設けられる。各パッチ画像を示す記号のうち、Ｖは垂直方向の縞、Ｈは水平方向の縞を表わしている。ＶまたはＨの後ろの２つの数字は、１番目の数字は各パッチ画像に含まれているＡＣ成分とＤＣ成分の組合せ条件を示す通し番号で、２番目の数字は各パッチ画像の濃淡の空間周波数を表わす通し番号である。すなわち、Ｈ１．１〜Ｈ１．８はＡＣ成分とＤＣ成分が同じで濃淡の空間周波数のみが異なるパッチ画像である。Ｈ２．１〜Ｈ２．８・・・、Ｖ９．１〜Ｖ９．８も各々同様である。また、Ｖ１．１とＨ１．１・・・、Ｖ９．８とＨ９．８はそれぞれＡＣ成分とＤＣ成分と濃淡の空間周波数が同じで、縞の方向のみが異なるパッチである。
【０１４７】
ＤＣパッチ画像群１７２における各パッチ画像は、０ｌｐ／ｍｍ以外の空間周波数成分を含まないパッチ画像、すなわち濃度変化のない平坦な濃度のパッチ画像であり、後述するように、ＡＣパッチ画像群１７１の測定値を規格化するするためのものである。Ｄ０〜Ｄ２５５はそれぞれ濃度の異なるパッチである。
【０１４８】
図２５に示したテストパターン画像データのレイアウトは一例であって、どの位置にどのようなパッチ画像が配置されているかがわかれば他のレイアウトであってもかまわない。また、それぞれのパッチ画像群が別々の複数枚のテストパターン画像として形成されるものであってもよい。各パッチ画像群に含まれるパッチ画像群の数や特性、大きさ等は任意である。特にこの例では、９種類のＡＣ成分／ＤＣ成分の組合せを持つＡＣパッチ画像を用いているが、ＡＣパッチ画像を追加することにより精度を向上させることが可能である。
【０１４９】
なお、例えば、図２５に示したようなテストパターン画像データのレイアウトは、あらかじめ作成してテストパターン保持部５などに記憶させておいてもよい。複数のテストパターン画像データのレイアウトを記憶させておき、選択して用いるように構成することもできる。テストパターン画像データをテストパターン保持部５に記憶させておくことにより、テストパターン画像データのレイアウトは１度行なえばよいことになり、２回目以降にテストパターン画像を出力する際の高速化を図ることができる。これにより、画像データ出力工程１４０、出力画像の測定工程１５０、空間周波数処理係数の決定工程１６０の一連の作業時間を短縮することができ、画像出力装置の補正を簡単に行なうことができるようになる。
【０１５０】
Ｓ１４２において、例えば、図２５に示したようなテストパターン画像データのうちＡＣパッチ画像群１７１の画像データを作成する。上述のように図２５において、ＡＣパッチ画像群１７１の各パッチに付与された記号Ｖ１．１〜Ｖ１．８、Ｖ２．１〜Ｖ１．８・・・Ｖ９．１〜Ｖ９．８は垂直方向のライン、Ｈ１．１〜Ｈ１．８、Ｈ２．１〜Ｈ１．８、・・・、Ｈ９．１〜Ｈ９．８は水平方向のラインを表わしており、ＶまたはＨの後ろの数字はＡＣ成分／ＤＣ成分の組合せ条件の種類と各ラインの幅を画素数で表わしている。これらの各パッチ画像は、上述の空間周波数処理部の第１の例において説明した図８に示す例と同様に構成することができる。なお、図８（Ａ）はＶ１．２のパッチ画像の１部、図８（Ｂ）はＨ１．３のパッチ画像の一部に相当する。ここでは、黒画素のデータ値を１６０、白画素のデータ値を９６としている。もちろん、図２１に示したような複数の空間周波数成分が含まれるパッチを使用してもよい。
【０１５１】
図２６、図２７は、空間周波数処理部の第２の例において用いるＡＣ画像群の各パッチ画像の特性値の一例である。図２６、図２７には、ＡＣパッチ画像群１７１について、その空間周波数、黒画素のデータ値、白画素のデータ値、信号コントラスト値を示している。空間周波数は、１ｍｍあたり何本のラインが含まれるかを示したものであり、単位はｌｐ／ｍｍで示している。Ｖ１．１、Ｈ１．１であれば１画素おきに白と黒のラインが交互に現れるため、想定しているプリンタの解像度である１ｍｍあたり１０ｄｏｔ中に５本のラインが含まれていることを示している。ＤＣ成分値は平均濃度を示しており、ＡＣ成分値は黒画素のデータ値および白画素のデータ値とＤＣ成分値の差分を示している。なお、濃度は最も黒の濃度が濃い場合を２５５、最も黒の濃度が薄い白部分の場合を０としている。信号コントラスト値は後述する。
【０１５２】
例えば、図２６、図２７に示した各パラメータに従って、ＡＣパッチ画像群１７１の各パッチ画像の画像データを作成し、作成した各パッチ画像データを例えば図２５に示すようなテストパターン画像データのレイアウトに従って、埋め込む。
【０１５３】
次に、Ｓ１４３において、ＡＣパッチ画像群１７１の各パッチ画像の画像データから特性値を算出する。ＡＣパッチ画像群１７１における各パッチ画像の各画素のデータ値は、例えばＶ１．１においては、黒画素値＝１６０、白画素値＝９６というように設定した。これは、別の表現をすると、ＤＣ成分値＝１２３、ＡＣ成分値＝±３２と表わすことができる。ＤＣ成分値、ＡＣ成分値で表現したパッチの設定値は図２６、図２７に示した。
【０１５４】
ここで、ＡＣパッチ画像群１７１の各パッチ画像の画像データの特性値として、ＤＣ成分値と信号のコントラスト値を考える。ＤＣ成分値は、図２６、図２７に示した各ＤＣ成分値である。信号のコントラスト値は、上述した各ＡＣ成分値をもとに補正を行なう。画像の場合には信号は正弦波というよりも矩形波であるため、実効コントラスト値は正弦波のときよりも向上する。このことを考慮してコントラスト値の補正を行なう。実際の補正は、
信号コントラスト値＝ＡＣ成分値×（４／３．１４）
によって計算することができる。このような補正を行なって決めた各パッチ画像のコントラスト値を図２６、図２７の信号コントラスト値の欄に示している。
【０１５５】
Ｓ１４４において、ＤＣパッチ画像群１７２のパッチ画像の画像データを作成する。図２５において、ＤＣパッチ画像群１７２の各パッチ画像に付された記号Ｄ０〜Ｄ２５５において、記号ＤはＤＣ成分すなわち０ｌｐ／ｍｍのみの空間周波数成分からなるパッチ画像であることを示し、後ろの数字はＤＣ成分のデータ値を表わしている。例えば、ＤＣ６４のパッチ画像は、２００画素×２００画素の全画素のデータ値が６４である均一なパッチ画像を意味する。
【０１５６】
ここでは、データ値＝０からデータ値＝２５５までを、ほぼ均等に分割し、９個のパッチ画像を用意している。各パッチ画像の画像データを作成し、作成した各パッチ画像のデータを、例えば図２５に示したようなテストパターン画像データのレイアウトにしたがって埋め込む。これにより、テストパターン画像データが完成する。完成したテストパターン画像データは、テストパターン保持部５により保持される。カラープリンタでは、このような同一のテストパターン画像データが、Ｙ、Ｍ、Ｃ用に３色分、さらにはＫ用にもう１色分用意される。あるいは１つのテストパターン画像データが３色あるいは４色共通に使用される。
【０１５７】
画像データ出力工程１４０の最後として、Ｓ１４５においてテストパターン画像データを出力装置から出力させ、テストパターン画像を作成する。テストパターン画像の出力に際しては、まず出力装置を単色モードに設定し、テストパターン画像の出力を行なう。例えば、Ｙ単色モードに設定されているとき、テストパターン保持部５からはＹ用画像信号としてテストパターン画像データが出力される。テストパターン保持部５から出力されたテストパターン画像データは、図２に示す階調変換部１３においてルックアップテーブルを用いて階調変換処理が施された後、出力装置から出力される。出力装置では、Ｙ用の階調変換処理が施されたテストパターン画像データを受け取り、Ｙ単色モードにてその出力の画像を出力する。Ｍ、Ｃに対しても同様に、出力装置をＭ単色モード、Ｃ単色モードに設定しながら、それぞれのテストパターン画像を得、計３枚のテストパターン画像を得る。
【０１５８】
得られたテストパターン画像は、例えば図２５と同様のレイアウトの画像である。ただし、各パッチ画像は、出力装置の特性を反映しており、テストパターン画像データを忠実に反映した画像とは限らない。以下、得られた３枚のテストパターン画像のそれぞれに対して、出力画像の測定工程１５０および空間周波数処理係数の決定工程１６０を行なう。説明の都合上、そのうちの１色についてのみ説明するが、３色それぞれについて同様の処理を行なう。
【０１５９】
出力画像の測定工程１５０では、テストパターン画像データを出力装置から出力したテストパターン画像の測定を行ない、測定結果を解析する。まずＳ１５１において、テストパターン画像中のＡＣパッチ画像群について、微小領域ごとに反射率を測定する。反射率の測定には、例えばデジタルカメラなどを用いることができる。デジタルカメラを用いた場合は同時に広い領域の反射率の測定が可能である。図２８は、反射率の測定の別の具体例の説明図である。図中、１８１はデジタルカメラ、１８２はパッチ画像である。例えば、テストパターン画像の２０ｍｍ四方を１０００画素×１０００画素の画像データとして取り込めるようにデジタルカメラ１８１の位置を調整して、テストパターン画像の撮像を行なう。これは、テストパターン画像の反射率を０．０２ｍｍ間隔で測定していることと等価である。この際、撮像を行なうときに、デジタルカメラの軸に対して、テストパターン画像が傾かないようにする。また、補色フィルタを使用して撮像するとよい。もちろん、図１０に示したようなアパーチャ８１による走査を行なう方式を用いたり、ラインタイプの撮像素子で走査してもよい。また、図２８に示す構成を上述の空間周波数処理部の第１の例において用いることもできる。
【０１６０】
続いてＳ１５２において、テストパターン画像中のＤＣパッチ画像群の各パッチ画像についても微小領域ごとに反射率を測定し、平均反射率を求める。この場合も、例えば図２８に示すようにデジタルカメラ１８１を用いたり、図１０に示すような走査方式を用いることができる。ＤＣパッチ画像群は上述のようにＤＣ成分のみの均一なパッチ画像なので、ＤＣパッチ画像を撮像し、撮像した測定データ画像の平均値を得る。この際、ＡＣパッチ画像群の測定に際してフィルタを使用した場合には、ここでも同一のフィルタを使用する。
【０１６１】
このようにして得られた測定値に対して解析を行なう。まずＳ１５３において、ＡＣパッチ画像群の各パッチ画像における測定値を規格化する。上述の空間周波数処理部の第１の例と同様に、Ｓ１５２で測定して得られたＤＣパッチ画像群の各パッチ画像の平均反射率のデータと、テストパターン画像データにおけるＤＣパッチ画像群の各パッチ画像のデータ値とを対応づけると、例えば図１１に示すような関係が得られる。まずこの対応関係ｆdc（ｘ）を近似値により求める。ＡＣパッチ画像群の各パッチ画像に対するデジタルカメラ１８１による測定値を、この対応関係ｆdc（ｘ）を用いて変換する。これにより、測定結果を処理演算部１で取り扱うＹＭＣ画像信号と同等に扱うことが可能となる。
【０１６２】
図２９は、空間周波数処理部の第２の例において測定された測定データ列の並び方向の説明図である。測定データ列は、水平方向の縞のＡＣパッチ画像に対する測定データ画像の場合は、図２９（Ａ）に示すように垂直方向の断面の測定データ列であり、同様に垂直方向の縞のＡＣパッチ画像に対する測定データ画像の場合は、図２９（Ｂ）に示すように水平方向の断面の測定データ列である。ＡＣパッチ画像に対する測定データ画像の断面の測定データ列は、例えば上述した図１２（Ａ）に示すようになる。測定データ画像を例えば図１１に示すように得られた対応関係ｆdc（ｘ）を用いて出力装置に依存した画像信号値に変換して、規格化データ画像を作成する。規格化データ画像の断面の規格化データ列は、例えば図１２（Ｂ）に示すようになる。断面の方向は、図２９に示すもので、測定データ画像に対するものと同様である。この規格化データ画像は、図２において色空間変換部１１から出力されるＹＭＣ信号と同等のものである。ここではＹＭＣ信号は８ｂｉｔ２５６階調としているが、これに限られるものではない。なお、上述のように、規格化データ画像において値が小さい方が白に近く、値が大きい方が色濃度が濃くなることを示している。そのため、反射率の値の大小とは逆になっている。
【０１６３】
次にＳ１５４において、Ｓ１５３で規格化した測定データ画像に対して整形処理を行ない、整形データ画像を作成する。この処理は、空間周波数処理部の第１の例における図１３に示した処理と基本的に同じである。ここでは、図２８に示したようにデジタルカメラ１８１によって上述のように２０ｍｍ四方のパッチ画像を０．０２ｍｍ間隔で測定するものとし、２次元の処理を行なう場合について図１３を用いて説明する。この場合には、規格化画像データは、例えば１０００画素×１０００画素の測定データからなっている。図１３のＳ９１において、これら１０００画素×１０００画素のデータからなる規格化データ画像の中央部から５１２画素×５１２画素のデータを切り出す。次にＳ９２において、Ｓ９１で切り出した５１２画素×５１２画素のデータを２次元離散フーリエ変換（２Ｄ−ＤＦＴ）処理する。これにより、空間周波数成分が算出される。Ｓ９３において、Ｓ９２で算出された空間周波数成分のうち、不要な高周波成分を除去し、その後、Ｓ９４において２次元離散フーリエ逆変換処理（２Ｄ−ＩＤＦＴ）する。このような整形処理によって、ノイズおよび高周波成分が除去されて、例えば、図１２（Ｃ）に示すような、断面においてＳｉｎ波形を基本とした整形データ列を持つ整形データ画像が得られる。ここで、断面の方向は図２９に示すもので、測定データ画像に対するものと同様である。
【０１６４】
Ｓ１５５において、ＡＣパッチ画像群に対応する整形データ画像をもとに、それらの特性値としてコントラスト値を算出する。図３０は、空間周波数処理部の第２の例における整形データ画像に対する特性値の算出手法の一例の説明図である。まず、整形データ画像より、図３０（Ａ）に示すように縞と直交する方向に沿って整形データ列を取り出す。整形データ列は複数取り出すことで特性値の算出精度を上げることができる。図３０（Ａ）では、整形データ列として、１個の整形データ画像から６本の整形データ列を取り出している。整形データ列は、上述した図１２（Ｃ）に示すようなＳｉｎ波形を基本としたデータ列である。図３０（Ｂ）に示すように、１本の整形データ列に対して、破線で示したようなウインドウ領域を設定し、このウインドウ領域内での最大値と最小値を検出し、（最大値−最小値）を算出する。ウインドウ領域を移動しながら、（最大値−最小値）を繰り返し算出する。１個の整形データ画像から取り出した、全ての整形データ列につて繰り返し算出した（最大値−最小値）の平均値をコントラストとする。整形データ列からは、極大値と極大値の間、あるいは極小値と極小値の間など、特徴の間のデータ数を計数することによって周期を求めることができる。測定ピッチは、０．０２ｍｍであったから、０．０２ｍｍ×データ数で周期が求められる。また、周期の逆数を計算することにより、空間周波数が求められる。ここで、ウインドウ領域のサイズは、ＡＣパッチ画像群の各パッチ画像に含まれる縞状の線の２〜３本分程度に設定すればよい。
【０１６５】
図３１は、空間周波数処理部の第２の例においてテストパターン画像から得られたＡＣパッチ画像のうちＶ１．１〜Ｖ１．８およびＨ１．１〜Ｈ１．８の各パッチ画像の周期および空間周波数の一例の説明図、図３２は、同じく得られたコントラスト値の一例を示すグラフである。図３１には、上述のようにして算出されたＡＣパッチ画像群の各パッチ画像に含まれる縞状の線の周期を長さとデータ数で示すとともに、算出される空間周波数と、使用したウインドウ領域のサイズをデータ数で示している。また、各パッチ画像ごとに算出された空間周波数とコントラストの関係の一例を図３２に示している。図３２において、黒丸はＨ１．１〜Ｈ１．８の各パッチ画像における空間周波数とコントラストの関係を示し、白丸はＶ１．１〜Ｖ１．８の各パッチ画像における空間周波数とコントラストの関係を示している。
【０１６６】
図３１、図３２にはＡＣパッチ画像群のうちＶ１．１〜Ｖ１．８およびＨ１．１〜Ｈ１．８の各パッチについて、周期、空間周波数、コントラスト等を示したが、ＡＣパッチ画像群中の他の各パッチ画像についても同様にコントラストを算出する。
【０１６７】
空間周波数処理係数の決定工程１６０では、テストパターン画像データのＡＣパッチ画像群の各パッチ画像のデータのコントラスト値と、それに対応するテストパターン画像のコントラスト値との関係から出力装置の空間周波数伝達特性を求め、それを補正するための空間周波数処理部１２における処理係数を算出する。
【０１６８】
Ｓ１６１において、出力装置の入出力対応特性を算出する。上述のように、出力装置の空間周波数伝達特性を求めて補正処理系の目標空間周波数伝達特性を決定してもよいが、ここでは、出力装置の入出力コントラスト特性から補正処理系の目標空間周波数伝達特性を決定する例を示す。
【０１６９】
図３３は、出力装置における入出力コントラストの関係の一例を示すグラフである。例えばＤＣおよび空間周波数を固定し、コントラストを示すＡＣ成分について、その入出力コントラストの対応づけを考える。例えば、ＡＣパッチ画像群中から、水平方向の縞を持ち、ＤＣ＝１２８で濃淡の空間周波数＝５ｌｐ／ｍｍについて、ＡＣパッチ画像Ｈ１．１（入力信号コントラスト値＝±４０．５）、Ｈ２．１（入力信号コントラスト値＝±２０．４）、Ｈ３．１（入力信号コントラスト値＝±８１．５）を用いて行なうと、例として図３３に示す対応関係ｆout-in（Ｘ）が得られる。このような対応関係ｆout-in（Ｘ）を、各ＤＣ（＝６４，１２８，１９２）および各空間周波数（＝５，２．５，１．２５，０．６２５）の組み合わせごとに求める。
【０１７０】
次にＳ１６２において、各ＡＣ成分とＤＣ成分の組合せ条件における、出力装置の空間周波数伝達特性を補正する目標空間周波数伝達特性を決定する。図３３に示すようなＤＣおよび空間周波数を特定したときの入出力コントラストの対応関係ｆout-in（Ｘ）から、ＡＣパッチ画像群中に含まれる９通りの各々のＡＣ成分とＤＣ成分の組合せ条件下での目標空間周波数伝達特性を算出する。
【０１７１】
図２に示したように、処理演算部１に入力された画像データは、空間周波数処理部１２にて空間周波数伝達率Ｆ１の空間周波数処理を施された後に、出力部２に送られ、空間周波数伝達率Ｆ２である画像形成工程を経て、出力画像が作成される。したがって、入力された画像データの有する空間周波数成分をＳ１、出力画像の有する空間周波数成分をＳ２とすると、Ｓ２＝Ｓ１×Ｆ１×Ｆ２となる。
【０１７２】
上述のごとく、出力装置によらず同一の鮮鋭度を得るためには、空間周波数処理部１２における空間周波数伝達率Ｆ１は、画像形成工程の空間周波数伝達率Ｆ２を補正するものであればよい。すなわち、空間周波数処理部１２における空間周波数伝達率Ｆ１＝１／Ｆ２となるように、空間周波数処理部１２に設定する処理係数を設計すればよい。
【０１７３】
入出力コントラストの対応関係ｆout-in（Ｘ）を用いると、例えば出力コントラスト±４０．５を得るために必要な入力コントラスト値を求めることができる。ここで求まる入力コントラスト値を±Ｄ1 とする。ＡＣパッチ画像Ｈ１．１はＤＣ＝１２８、信号コントラスト値＝±４０．５、空間周波数＝５ｌｐ／ｍｍのパッチ画像なので、ＡＣパッチ画像を出力装置の特性の影響を受けずに出力するのには、ＡＣ成分をあらかじめＤ1 ／４０．５倍しておけばよい。この時のＤ1 ／４０．５がＤＣ＝１２８、入力コントラスト＝±４０．５、空間周波数＝５ｌｐ／ｍｍにおける補正系の目標空間周波数伝達特性である。同様にして、ＤＣ＝１２８、空間周波数＝２．５ｌｐ／ｍｍ、１．２５ｌｐ／ｍｍ、０．６２５ｌｐ／ｍｍについても、それぞれの対応関係ｆout-in（Ｘ）から信号コントラスト値が±４０．５の時の入力信号コントラスト値Ｄ2 〜Ｄ4 を求める。
【０１７４】
図３４は、基本ＡＣパッチ画像群の各パッチ画像から求められた空間周波数伝達特性に基づいて設計された空間周波数処理部の第２の例における空間周波数伝達率の一例を示すグラフである。上述の手順に従って求めた、ＤＣ＝１２８、入力コントラスト値＝±４０．５におけるＤn ／４０．５の値を伝達率とし、空間周波数と対応づけてグラフ化すると、例えば図３４の破線で示すようになる。これが垂直方向の目標周波数伝達特性である。このままでは高い空間周波数帯域においてノイズが極端に強調されてしまう。ノイズ成分の強調を抑えるため、解像度１０ｄｏｔ／ｍｍのカラープリンタでは、そのナイキスト周波数である５ｌｐ／ｍｍにおいて伝達率が０になるように設定している。このようにして実際に設計した補正系の空間周波数伝達特性のグラフを図３４において実線で示している。
【０１７５】
ＤＣ＝１２８、入力コントラスト値＝±４０．５における垂直方向の目標空間周波数伝達特性の設計値の求め方について説明した。同様にして、入力コントラスト値＝±２０．４、±８１．５の場合についても目標空間周波数伝達特性の設計値を算出する。また、同様に水平方向の目標空間周波数伝達特性の設計値を算出する。さらに、ＤＣ＝６４、１９２についても、各入力コントラスト値について図３４に示すような補正系の空間周波数伝達特性を算出する。図３５は、空間周波数処理部の第２の例において算出された垂直方向の目標空間周波数伝達特性の一例を示すグラフである。図中の各線は、特定のＡＣ成分とＤＣ成分の組合せ条件下での目標空間周波数伝達特性の設計値を表わしている。このような処理を、垂直方向および水平方向についてそれぞれ行なう。
【０１７６】
上述の例において、図３３に示すようなｆout-in（Ｘ）から求めた入力コントラスト値を用いるとき、ＤＣ成分の値が大きい場合や小さい場合に入力コントラスト値が大きいと、画像信号の最大値が非常に大きな値となったり、あるいはマイナスの値となる場合が発生する。一般に、出力部２で受け付ける入力信号の値は制限されている。例えば、８ビット入力のプリンタの場合には、入力信号値は０以上２５５以下の整数値に限られており、それ以外の入力信号値を出力部２に入力することはできない。入力信号がこのような範囲を越えたあるいは下回った値を含む場合、一般のプリンタなどでは入力する前に不適当な値を取り除く処理が行なわれている。広く利用されている処理の例として、例えば打ち切り処理がある。打ち切り処理を行なうと、入力信号の値が２５５を越える場合には２５５とされ、マイナスの場合には０とされる。
【０１７７】
図３６は、入力信号値が制限範囲を越える場合の画像の出力特性の説明図、図３７は、入力信号値が制限範囲を下回る場合の画像の出力特性の説明図である。ここでは、出力部２における入力信号値の制限範囲が０〜２５５であるものとして示している。例えば図３６（Ａ）に示すようにＤＣ成分および入力コントラスト値が大きく、入力信号値の一部が制限範囲を越えている。制限範囲を越えた部分について、例えば打ち切り処理を行なうと図３６（Ｂ）に示すようになり、このような入力信号値が出力部２に与えられることになる。出力部２は、図３６（Ｂ）に示すような入力信号値に応じた画像を形成する。形成された画像の反射率は図３６（Ｃ）に示すようになる。ここで、入力信号値が大きいと濃度が高くなるため反射率は低下する。この例では、入力信号値が２５５のとき反射率は０になっている。
【０１７８】
また、図３７（Ａ）に示す例では、ＤＣ成分が小さく、入力コントラスト値が大きいため、入力信号値の一部が制限範囲を下回り、マイナスとなっている。打ち切り処理を行なえば、入力信号値は、図３７（Ｂ）に示すようにマイナスの部分が０に変換される。図３７（Ｂ）に示すような入力信号値が出力部２に与えられ、画像が形成される。この時形成された画像の反射率は、図３７（Ｃ）に示すようになる。
【０１７９】
図３６（Ｃ）および図３７（Ｃ）からわかるように、入力信号値が出力部２の制限範囲を越えあるいは下回る部分が存在する場合、実際に出力された画像においては、反射率の振幅が制限されて小さくなっており、所望のコントラストを実現できていないことがわかる。さらに、出力された画像の局所的な反射率の平均値も変化していることになる。これは、出力画像の色が入力信号から期待される色と異なることを意味し、このような入力信号は適切ではない。
【０１８０】
図３８は、出力部において入力信号値が出力された画像の反射率に反映されない場合の説明図である。実際の画像出力装置では、高コントラスト、高周波数の入力信号に対して、入力信号の上限や下限近傍に対する応答が出力に反映されない現象が起こることがある。例えば図３８（Ａ）に示すように入力信号としては許可されている範囲内の入力信号が出力部２に与えられたとき、入力信号が小さく、淡い色の部分で、ある一定の反射率よりも明るい画像を出力できず、図３８（Ｂ）に示すような反射率の画像が出力されてしまう場合がある。このような現象は、例えばＣＲＴディスプレイ等で見られる。ＣＲＴディスプレイにおいては、特に高コントラストを再現する際に、図３８（Ｂ）とは逆に暗部の輝度が上昇することが実験から確認されている。このような現象が発生する原因は、画像出力装置の応答が高コントラスト、高周波数の入力信号に追従できなくなることによる。このような現象は、図３８に示すように、上述の図３６、図３７に示す現象と同様に、出力された画像の反射率の平均値を変化させてしまう。そのため、出力画像の色が期待される色と異なってしまうことがある。このような現象が発生する場合、入力信号が出力に反映される範囲が実際の制限範囲となる。
【０１８１】
このように、入力信号値が制限範囲を越えたり下回った場合には、出力画像における反射率の平均値が維持され、出力画像の色が期待される色と異ならないように、入力信号値を調整する必要がある。図３９は、入力信号値の調整処理の概念図である。例えば図３６（Ａ）に示したように、入力信号値が制限範囲を越える場合には、図３９（Ａ）に示すように、入力信号値が制限範囲を越えず、平均値が変化しないように、振幅を変化させる。同様に図３７（Ａ）に示すように入力信号値が制限範囲を下回る場合には、図３９（Ｂ）に示すように、入力信号値が制限範囲を下回らず、平均値が変化しないように、振幅を変化させればよい。図３８（Ａ）に示すように出力画像に反映できない入力信号値が存在する場合にも、実際の制限範囲を考慮して図３９（Ａ）または（Ｂ）に示すように調整処理を行なえばよい。この実施の形態では、出力部２としてカラープリンタを想定しているが、このような調整処理を行なうことによって、ＣＲＴディスプレイや液晶ディスプレイのような画像表示機器に関しても、画像の色の一致がとれ、鮮鋭度を保存することができる。
【０１８２】
図４０は、入力信号値の制限範囲を考慮した入出力コントラストの関係の一例を示すグラフである。一例として図３３に示した入力コントラストと出力コントラストとの関係を示すグラフｆout-in（Ｘ）について、出力部２における実際の制限範囲も考慮に入れ、このグラフを修正する。図３３はＤＣ成分の値と空間周波数を一定（ＤＣ＝１２８、空間周波数＝５ｌｐ／ｍｍ）としている。ＤＣ成分の値を一定とした場合、振幅すなわちコントラスト値が小さい場合には入力信号値は制限範囲内に収まるので調整を施す必要はないが、コントラスト値がある閾値以上となると、一部が制限範囲を越えたりあるいは下回ったりすることになる。そのため、コントラスト値が閾値以上では、図３９に示すように調整処理を行なってコントラスト値を一定値に制限する。すなわち調整処理は、閾値を境として出力コントラスト値を一定にすれば図３９に示したような調整を行なうことができる。このような調整処理により、例えば図４０に示すｆ’out-in（Ｘ）を得る。このｆ’in-out（Ｘ）を用いて出力コントラスト値から入力コントラスト値を求めると、求められた入力コントラスト値は０から閾値までの値に収まることになる。
【０１８３】
次にこの閾値を求める方法について一例を説明する。図４１は、入力信号コントラストに対する出力画像の実際の平均反射率の一例を示すグラフである。図中、破線は目標とする平均反射率であり、ここでは一例としてＳ１５２において測定したＤＣ＝１２８のパッチ画像の平均反射率を表わしている。また一点鎖線は平均反射率から５％のズレを表わしている。Ｓ１５４において求めたＤＣ＝１２８、空間周波数＝５ｌｐ／ｍｍのＡＣパッチ画像Ｈ１．１、Ｈ２．１、Ｈ３．１の測定データ画像に整形処理を施した画像の平均反射率をプロットすることによって図４１に示すようなグラフが得られる。
【０１８４】
上述のように、入力コントラストを大きくしてゆくと、入力信号の一部が制限範囲を越え、あるいは下回り、平均反射率が変化する。図４１に示した例では、入力コントラストが大きくなると平均反射率が低下し、平均的な色合いが変化することがわかる。入力コントラストが１０２の時に平均反射率から５％ずれる。これを越える場合、平均的な色合いをこれ以上変化させないようにするため、調整処理を行なうこととする。すなわち、この時の入力コントラスト値を閾値として用い、上述のｆin-out（Ｘ）を修正すればよい。図４０では、図４１から求められた入力コントラスト値＝１０２を閾値とし、図３３に示すｆin-out（Ｘ）を修正してｆ’in-out（Ｘ）を得ている。なお、上述の例では平均反射率が５％以上ずれないように調整処理を行なったが、このずれ量の設定は任意に設定可能である。
【０１８５】
このようにして図３３に示すＤＣ＝１２８、空間周波数＝５ｌｐ／ｍｍの時の入力コントラストと出力コントラストとの対応関係ｆin-out（Ｘ）を、入力コントラストが閾値（＝１０２）を越えないように修正した対応関係ｆ’in-out（Ｘ）が得られた。同様にして、各ＤＣ（＝６４，１２８，１９２）および各空間周波数（＝５，２．５，１．２５，０．６２５）の組み合わせごとに修正した対応関係ｆ’in-out（Ｘ）を求める。
【０１８６】
以下同様にして、Ｓ１６２において各ＡＣ成分とＤＣ成分における目標空間周波数伝達特性を算出すればよい。例えば上述のようにＤＣ＝６４、信号コントラスト値±４０．５の場合の目標空間周波数伝達特性は、ＤＣ＝６４で空間周波数＝５，２．５，１．２５，０．６２５の時の４つの対応関係ｆ’in-out（Ｘ）を用い、それぞれのグラフから出力コントラストが４０．５の時の入力コントラスト値Ｄ1 〜Ｄ4 を求め、Ｄ1 ／４０．５，Ｄ2 ／４０．５，Ｄ3 ／４０．５，Ｄ4 ／４０．５を目標空間周波数伝達特性とすればよい。このようにして求められた各空間周波数ごとの目標空間周波数伝達特性をグラフ化すると、すべての空間周波数で制限範囲を越えない場合には図３４に破線で示すグラフとなる。
【０１８７】
しかし、コントラスト値が大きい場合には、調整処理の影響を受ける。この例では各ＡＣパッチの画素値は図２６、図２７に示すように０〜２５５に収まるように設計しているが、例えばＤＣ＝６４で信号コントラスト値＝±８１．５の場合や、ＤＣ＝１９２で信号コントラスト値＝±８１．５の場合などでは、出力部２の特性や仕様によっては制限範囲を越えあるいは下回る場合がある。このような入力信号値が制限範囲を越えあるいは下回るＤＣ成分、信号コントラスト値、空間周波数の組み合わせの目標空間周波数伝達特性は、対応関係ｆ’in-out（Ｘ）から求められる入力コントラスト値Ｄが、図４１に示すようにして求められた閾値となるので、（閾値／信号コントラスト値）となる。図３３に示す対応関係ｆin-out（Ｘ）と図４０に示すｆ’in-out（Ｘ）とを比較してわかるように、調整処理を行なった方が出力コントラスト値から得られる入力コントラスト値Ｄは小さくなるので、目標空間周波数伝達特性の値も小さくなる。
【０１８８】
図４２は、空間周波数処理部の第２の例における空間周波数伝達率の別の例を示すグラフである。図中、破線は調整処理を行なわない入力コントラスト値と出力コントラスト値との対応関係ｆin-out（Ｘ）を用いて求めた目標空間周波数伝達特性である。図４２に示した例では、空間周波数が高くなるに従って入力信号値に対する応答が出力に反映されにくくなる場合を想定し、空間周波数が５ｍｍ／ｌｐのときに対応関係ｆ’in-out（Ｘ）のうち調整処理を行なった部分を使用した場合を示している。図４２に示すように、調整処理を行なった対応関係ｆ’in-out（Ｘ）を用いることによって、空間周波数特性は調整処理を施した部分で値が抑えられる。そのため、上述のように空間周波数が高い部分で発散することはなく、ある程度の値に落ちつくことになる。
【０１８９】
上述のように、図３４に破線で示した目標空間周波数伝達特性では、高い空間周波数帯域においてノイズが極端に強調されてしまうため、図３４に実線で示したような目標空間周波数伝達特性に修正する。調整処理を行なった対応関係ｆ’in-out（Ｘ）を用いて得た、例えば図４２に示すような目標空間周波数伝達特性についても、同様に修正を施せばよい。しかし図４２に示すように、調整処理を行なった対応関係ｆ’in-out（Ｘ）を用いた場合、高い空間周波数帯域における空間周波数伝達率が小さくなっているため、極端にノイズが強調されることはない。そのため、図４２に示すような目標空間周波数伝達特性をそのまま用いるようにしてもよく、目標空間周波数伝達特性の修正処理を省略して高速化することが可能である。
【０１９０】
このような目標空間周波数伝達特性を、ＤＣ成分、ＡＣ成分の組み合わせごとに求め、さらに垂直方向および水平方向について求める。以下の説明では、一例として図３４に実線で示したように目標空間周波数伝達特性に修正を加え、図３５に示すようにそれぞれのＤＣ成分、ＡＣ成分の組み合わせごとに求めたものとし、さらに垂直方向および水平方向について求めたものとして説明を続ける。
【０１９１】
このようにして設計された空間周波数処理部１２の空間周波数伝達特性を実現するように、Ｓ１６３、Ｓ１６４、Ｓ１６５において、空間周波数処理の処理係数を算出する。空間周波数処理部１２として例えば図２２に示す構成を用いる場合、Ｓ１６３において基本補正部１１１のフィルタ係数を、Ｓ１６４において補助補正部１１２のフィルタ係数をそれぞれ決定する。また、Ｓ１６５において図２３中の２次元ルックアップテーブル１３３の変換テーブルを決定する。
【０１９２】
Ｓ１６３において、まず、基本補正部１１１の空間周波数伝達特性を決める。基本補正部１１１の空間周波数伝達特性は、例えばＳ１６２において各ＡＣパッチ画像から求めた、垂直方向と水平方向の目標空間周波数伝達特性の設計値の各々の平均特性をもとに、２次元面補間したものとする。次に基本補正部１１１の空間周波数伝達特性に逆フーリエ変換を施し、逆フーリエ変換をした結果を総和が１となるように規格化することで、目標とする空間周波数伝達特性のフィルタ係数を算出することができる。このようなフィルタ係数の算出方法は種々開発されており、適宜用いればよい。
【０１９３】
Ｓ１６４において、補助補正部１１２のフィルタ係数を決定する。まず、各ＡＣ成分とＤＣ成分の組合せ条件下における空間周波数処理部１２の空間周波数伝達特性の設計値から、Ｓ１６３において設定した基本補正部１１１の空間周波数伝達特性を引き、各ＡＣ成分とＤＣ成分の組合せ条件下における差分特性を算出する。図４３は、空間周波数処理部の第２の例における差分特性の一例を示すグラフである。例えば図３５に示した垂直方向の目標空間周波数伝達特性の設計値をもとに、Ｓ１６３において設定した基本補正部１１１の空間周波数伝達特性を引き、差分特性を求めると、図４３に示すようなものとなる。図中の黒丸は、ＡＣ成分とＤＣ成分と空間周波数の特定の組合せを示しており、ＡＣパッチ画像群中の各パッチ画像と対応している。
【０１９４】
この差分特性を代表する１つの空間周波数伝達特性を求め、それに基づき補助補正部１１２のフィルタ係数を決定する。図４４は、水平方向および垂直方向の空間周波数を伝達率の関係の説明図である。各ＡＣ成分とＤＣ成分の組合せ条件下における、水平方向の空間周波数が０．６２５、１．２５、２．５、５ｌｐ／ｍｍおよび垂直方向の空間周波数が０．６２５、１．２５、２．５、５ｌｐ／ｍｍの時の目標空間周波数伝達特性は、図４４における垂直方向の空間周波数軸と水平方向の空間周波数軸上の曲線として現われる。これらの曲線上の値をそれぞれの次元として、図４４に示すように８次元空間中の１点として表わすことができる。この例では、各ＡＣ成分とＤＣ成分の組合せの条件について行なうと、８次元空間中の点が９点求まる。この際、空間周波数に応じて人間の視覚特性（ＶＴＦ：Ｖｉｓｕａｌ  Ｔｒａｎｓｆｅｒ  Ｆｕｎｃｔｉｏｎ）を掛けて重み付けを行なってもよい。また、図４４の対応関係を用いると、８次元空間中の１点から水平方向の空間周波数０．６２５、１．２５、２．５、５ｌｐ／ｍｍおよび垂直方向の空間周波数０．６２５、１．２５、２．５、５ｌｐ／ｍｍにおける空間周波数伝達特性に展開することも可能である。
【０１９５】
次に、各ＡＣ成分とＤＣ成分の組合せ条件下で求めた８次元空間中の９点の最近傍を通り原点を通過する直線を算出する。ここで、８点の最近傍を通る直線とは、各点と直線の間の距離の自乗和が最も小さいものである。上述した問題は、方程式を解くことで求めることができる。また、例えば最急降下法を用いて反復計算を行ない、解を求めてもよい。求めた直線の、長さを１に規格化した方向ベクトルを８次元中の１点と見なし、図４４に示した対応関係を用いて水平方向および垂直方向における空間周波数が０．６２５、１．２５、２．５、５ｌｐ／ｍｍの場合の空間周波数伝達特性に展開する。このようにして直線の方向ベクトルより求めた空間周波数伝達特性を２次元面補間し、逆フーリエ変換を施したものを補助補正部１１２のフィルタ係数とする。ただし、２次元面補間を行なう際に、水平方向および垂直方向の空間周波数＝０の空間周波数伝達特性を０として面補間を行なう。
【０１９６】
Ｓ１６５において、各ＡＣ成分とＤＣ成分の組合せ条件下における空間周波数処理部１２の空間周波数伝達特性の設計値とＳ１６３において設定した基本補正部１１１の空間周波数伝達特性の差に、Ｓ１６４において算出した補助補正部１１２の空間周波数伝達特性を定数倍することによってさらに近づける。ここで最も近いとは両者の差の自乗の和を最小にすることである。各条件についてこの補助補正部１１２の空間周波数伝達特性の倍率を求め、さらに面補間を行なう。図４５は、空間周波数処理部の第２の例における補助補正部の空間周波数伝達特性に対する各条件における倍率の一例を示すグラフである。上述のようにして補助補正部１１２の空間周波数伝達特性の倍率を求めて面補間し、ＡＣ成分とＤＣ成分を軸にとった３次元グラフを作成すると、図４５に示すようになる。この３次元グラフを２次元ルックアップテーブル１３３の変換テーブルとすればよい。
【０１９７】
以上により、基本補正部１１１および補助補正部１１２のフィルタ係数、２次元ルックアップテーブル１３３の変換テーブルが得られた。これらの係数をそれぞれ設定し空間周波数処理部１２で空間周波数処理を行なう。空間周波数処理部１２では、入力画像信号１２１に対してＡＣ成分抽出器１３１でＡＣ成分を求め、ＤＣ成分抽出器１３２でＤＣ成分を求めて、求めたＡＣ成分とＤＣ成分から図４５に示すグラフをテーブル化した２次元ルックアップテーブル１３３を引いて倍率すなわち補助補正量増幅率信号１２４を得る。これを補助補正部１１２によるフィルタ処理によって得られる補助補正基本信号１２３と掛け合わせることによって、補助補正信号１２５が得られる。この補助補正信号１２５は、入力画像信号１２１のＡＣ成分、ＤＣ成分によって決定される例えば垂直方向のみを取れば図４３に示すような差分特性に近い補助特性による補正量を表わしている。加算部１１５においてこの補正量を基本補正信号１２２に加え、補正後の出力画像信号１２６を得る。このようにして、出力部２に応じて設定された処理係数にしたがって空間周波数処理を行なうことにより、出力部２の特性によらない空間周波数特性を有する出力画像を得ることができる。
【０１９８】
図４６は、空間周波数処理部の構成の第３の例を示すブロック図である。図中、図２２および図２３と同様の部分には同じ符号を付して説明を省略する。１１６は乗算部、１２７は基本補正量増幅率信号、１３６は２次元ルックアップテーブルである。この例では、基本補正信号１２２にも増幅率を掛け合わせる例を示している。
【０１９９】
２次元ルックアップテーブル１３６は、２次元ルックアップテーブル１３３と同様のテーブルであり、ＡＣ成分抽出器１３２から入力されるＡＣ成分信号１３５と、ＤＣ成分抽出器１３１から入力されるＤＣ成分信号１３４を受け取り、基本補正量増幅率信号１２７を生成して乗算部１１６へ出力する。乗算部１１６は、基本補正部１１から出力される基本補正信号１２２に２次元ルックアップテーブル１３６から出力される基本補正量増幅率信号１２７で与えられる増幅率を掛け合わせて出力する。加算部１１５は、乗算部１１６から出力される増幅率が掛けられた基本補正信号１２２と乗算部１１４から出力される補助補正信号１２５とを加算して出力画像信号１２６を出力する。
【０２００】
ここで新たに設定が必要となるのは、２次元ルックアップテーブル１３６である。上述の空間周波数処理部の第２の例において例えば図３５に示すような目標空間周波数伝達特性の設計値から基本補正部１１１のフィルタ係数を求める際に、すべての条件の目標空間周波数伝達特性の設計値を最もよく表わす特性とその倍率を例えば最小２乗法等を用いて算出する。ここで得られた特性をもとに、上述のようにして基本補正部１１１のフィルタ係数を決定する。また、各条件について得られた倍率に適当な補間処理等を施し、２次元ルックアップテーブル１３６の内容とすればよい。
【０２０１】
補助補正部１１２のフィルタ係数を決定する際には、各ＡＣ成分とＤＣ成分の組み合わせ条件における目標空間周波数伝達特性の設計値から、基本補正部１１１の空間周波数特性に上述のようにして２次元ルックアップテーブル１３６に設定した各条件ごとの倍率を掛け合わせたものを減算し、差分特性とする。得られた差分特性から補助補正部１１２のフィルタ係数および２次元ルックアップテーブル１３３の内容を決定する処理は、上述の空間周波数処理部の第２の例と同様である。
【０２０２】
このように、基本補正信号１２２にも増幅率を掛け合わせる構成とすることによって、さらに補正精度の向上が期待できる。
【０２０３】
図４７は、空間周波数処理部の構成の第４の例を示すブロック図である。図中、図２２、図２３と同様の部分には同じ符号を付して説明を省略する。１９１はデジタルフィルタ、１９２は乗算部、１９３は２次元ルックアップテーブル、１９４は加算部である。この例では、目標空間周波数伝達特性の設計値を空間周波数帯域ごとに分割し、各空間周波数帯域ごとにデジタルフィルタおよび増幅率によって補正し、全体として目標空間周波数伝達特性の設計値を実現しようとするものである。
【０２０４】
デジタルフィルタ１９１は、各空間周波数帯域ごとに設けられる。これらのデジタルフィルタ１９１は、既存の帯域分離用のデジタルフィルタ群を用いればよい。２次元ルックアップテーブル１９３は、図２３に示す２次元ルックアップテーブルと同様のテーブルであり、各空間周波数帯域ごとに設けられる。内容として、各ＡＣ成分およびＤＣ成分の組み合わせごとに、対応する空間周波数帯域における増幅率を有している。乗算部１９２は、各空間周波数帯域ごとに設けられ、対応する空間周波数帯域のデジタルフィルタ１９１の出力と２次元ルックアップテーブル１９３の出力を掛け合わせる。加算部１９４は、各乗算部１９２の出力をすべて加算して、出力画像信号１２６として出力する。
【０２０５】
図４８は、空間周波数処理部の構成の第５の例における空間周波数分割による目標空間周波数伝達特性の実現の説明図である。各デジタルフィルタ１９１の空間周波数特性は、図４８（Ａ）に示すようにそれぞれの特性を有している。２次元ルックアップテーブル１９３から得られる増幅率は、これらの特性を図中矢印の方向へそれぞれ変形させる。例えば、図３４に示したような目標空間周波数伝達特性が設計された場合、図４８（Ｂ）に示すように各空間周波数帯域における伝達率の和として、設計された目標空間周波数伝達特性を近似することができる。すなわち、所定の特性を有するデジタルフィルタ１９１の出力のそれぞれに対して、所定の増幅率を乗算し、すべてを加算することによって、図４８（Ｂ）に示すように目標空間周波数伝達特性を得ることができる。このような手法は、例えばウェーブレット等によって実現されている。
【０２０６】
この例の場合の２次元ルックアップテーブル１９３の内容は、図３４に示すような目標空間周波数伝達特性を各ＡＣ成分およびＤＣ成分ごとに水平方向および垂直方向とも求め、各空間周波数領域におけるデジタルフィルタ１９１の特性に応じて倍率を算出し、テーブルに設定してゆけばよい。
【０２０７】
図４９は、空間周波数処理部の構成の第６の例を示すブロック図である。図中、図２２、図２３と同様の部分には同じ符号を付して説明を省略する。２０１はデジタルフィルタ、２０２はパラメータデータベースである。この例では、図２２における基本補正部１１１、補助補正部１１２、乗算部１１４、加算部１１５、増幅率決定部１１３中の図２３に示す２次元ルックアップテーブル１３３による演算を、１つのデジタルフィルタ２０１によって行なうものである。上述のように空間周波数伝達特性は各ＡＣ成分およびＤＣ成分によって異なるため、固定のフィルタ係数では対応できない。そのため、各ＡＣ成分およびＤＣ成分の組み合わせに応じたフィルタ係数をパラメータデータベース２０２に蓄積しておく。蓄積しておくフィルタ係数は、各ＡＣ成分およびＤＣ成分の組み合わせに対応する目標空間周波数伝達特性から設計すればよい。
【０２０８】
ＡＣ成分抽出器１３１で抽出されたＡＣ成分、および、ＤＣ成分抽出器１３２で抽出されたＤＣ成分によってパラメータデータベース２０２からフィルタ係数を取り出す。そして、取り出したフィルタ係数をデジタルフィルタ２０１にセットして、入力画像信号１２１に対して補正処理を行ない、出力画像信号１２６を生成して出力する。
【０２０９】
この第６の例によれば、目標空間周波数伝達特性を近似する処理は含まれていないので、精度よく補正を行なうことができる。
【０２１０】
上述の第１の例におけるルックアップテーブル２３，２４や、第２ないし第５の例における各２次元ルックアップテーブルは、例えばニューラルネットワークで置き換えることも可能である。その場合には、ルックアップテーブルの作成工程はニューラルネットワークの学習工程として実現すればよい。
【０２１１】
また、上述の各例では、図２に示すように、ＹＭＣ空間において空間周波数処理を行なっているが、これに限らず、例えばＬ* ａ* ｂ* 色空間など、他の色空間のカラー画像に対して空間周波数処理を行なってもよい。この場合、出力はフルカラーモードで行ない、測定値をＬ* ａ* ｂ* 色空間に変換し、Ｌ* ，ａ* ，ｂ* のそれぞれの伝達特性を算出すればよい。また、カラー画像に限らず、モノトーンの画像についても本発明を適用することができる。
【０２１２】
図５０は、本発明を適用した画像編集・出力システムの一例を示す構成図である。図中、１０１はホストコンピュータ、１０２はカラープリンタ、１０３は画像データベース、１０４はフォトＣＤである。現在、デスクトップカラーパブリッシングの分野においては、例えば図５０に示すように、ホストコンピュータ１０１、カラープリンタ１０２、画像データベース１０３、フォトＣＤ１０４、などの装置による画像編集／出力システム１０を用いて、画像の編集出力処理が行なわれている。このようなシステムにおいて、それぞれの装置は専用信号線もしくはネットワークを介して接続されている。ホストコンピュータ１０１では、画像データベース１０３もしくはフォトＣＤ１０４、図示しないカラースキャナなどから入力された画像データや、ホストコンピュータ１０１上で作成されたコンピュータグラフィックス画像やテキスト情報を、画像編集ソフトウェアやレイアウトソフトウェアなどを用いて編集することによりディジタル画像データを作成し、作成したディジタル画像データをカラープリンタ１０２に転送し、カラープリンタ１０２により出力画像を得ている。
【０２１３】
このようなシステムにおいて、画像データベース１０３やフォトＣＤ１０４等から得られたカラー画像データは、カラープリンタ１２の特性によらずに再現したい。このような場合に、デジタル画像データをカラープリンタ１０２に転送する際に、カラープリンタ１０２に応じた出力処理を施す。この出力処理を行なう構成として、図２に示したような処理演算部１を適用することが可能である。
【０２１４】
例えば、図２に示すような処理演算部１をホストコンピュータ１０１内に設ける場合には、色空間変換部１１におけるダイレクト・ルックアップ・テーブルの処理係数、本発明で求める空間周波数処理部１２における処理係数、階調変換部１３におけるルックアップ・テーブルの処理係数は、出力装置情報としてあらかじめホストコンピュータ１０１内の図示せぬ出力装置情報保持部に出力装置ごとに保持させておく。そしてデジタル画像データを出力するカラープリンタ１０２を決定すると、そのカラープリンタ１０２に対応した出力装置情報を取り出して色空間変換部１１、空間周波数処理部１２、階調変換部１３に設定し、処理を行なうことにより、出力装置によらない空間周波数特性および色特性を有する出力画像を得ることが可能となっている。カラープリンタ１０２が１台だけ、ホストコンピュータ１０１に固定的に接続されている場合には、出力装置情報の選択的な使用は必要なく、各係数を最初から設定しておいて出力処理を行なえばよい。
【０２１５】
このシステムでは、画像データベース１３やフォトＣＤ１４から供給されるカラー画像データは入出力装置に依存しない形式で表現された画像データであってよく、望ましくはシステム内で、あるいはシステム間で共通した形式で表現されるデータであるとよい。
【０２１６】
ホストコンピュータ１０１には、係数算出部４を設けておくことができ、測定部３に対応する測定器を接続して、上述のような処理により空間周波数処理のための処理係数を求めることができる。測定部３をも予め設けておいてもよい。
【０２１７】
あるいは、係数算出部４および測定部３はホストコンピュータ１０１とは別の構成としておき、係数算出部４の出力をホストコンピュータ１０１に入力するように構成してもよい。この場合、係数算出部４の出力である空間周波数処理のための処理係数を、直接、ホストコンピュータ１０１に入力する構成のほか、ファイルの形式でホストコンピュータ１０１内の記憶装置にコピーしたり、ネットワークを介して配信することもできる。
【０２１８】
処理演算部１は、ホストコンピュータ１０１に設ける構成のほか、カラープリンタ１０２側に設けることもできる。例えば、ネットワークなどに接続されたカラープリンタでは、複数台のホストコンピュータからの出力要求を受け付けて画像を出力する。また、出力要求を行なうホストコンピュータ１０１側もどのプリンタに出力するかを決定しないで出力要求を出す場合もある。このような場合、カラープリンタ１０２側で、そのカラープリンタ１０２で固有の出力処理を行なえばよい。この場合、係数算出部４、測定部３はカラープリンタ１０２とは別の構成とし、係数算出部４の出力である空間周波数処理に使用する処理係数をカラープリンタ１０２に予め組み込んでおくことができ、係数算出部４、測定部３をカラープリンタ１０２に組み込む場合に比べてカラープリンタ１０２のコストを下げることができる。
【０２１９】
【発明の効果】
以上の説明から明らかなように、本発明によれば、カンや経験によらずに空間周波数処理の処理係数を決定することが可能となり、決定した処理係数を用いて画像データに空間周波数を施すことにより、出力画像装置の空間周波数伝達特性の影響を受けることなく、出力装置の種類によらず、同一の鮮鋭度の出力画像を得ることができるという効果がある。
【図面の簡単な説明】
【図１】  本発明の実施の一形態を示すブロック図である。
【図２】  処理演算部の構成の一例を示すブロック図である。
【図３】  デジタルフィルタ処理の一例の説明図である。
【図４】  種々の出力装置における空間周波数応答特性の一例の説明図である。
【図５】  空間周波数処理部の構成の第１の例を示すブロック図である。
【図６】  空間周波数処理部の第１の例において用いる処理係数を決定する処理の一例を示すフローチャートである。
【図７】  空間周波数処理部の第１の例において用いるテストパターン画像データのレイアウトの一例を示す説明図である。
【図８】  空間周波数処理部の第１の例において用いる基本ＡＣパッチ画像群のパッチ画像の一例の説明図である。
【図９】  空間周波数処理部の第１の例において用いる基本ＡＣパッチ画像群と補助ＡＣパッチ画像群の各パッチ画像の特性値の一例の説明図である。
【図１０】  反射率の測定の具体例の説明図である。
【図１１】  ＤＣパッチ画像群の各パッチ画像の平均反射率のデータとテストパターン画像データにおけるＤＣパッチ画像群の各パッチ画像のデータ値との対応関係の一例を示すグラフである。
【図１２】  ＡＣパッチ画像の測定値の規格化および成形処理による波形の変化の一例の説明図である。
【図１３】  規格化データ列に対する成形処理の一例を示すフローチャートである。
【図１４】  空間周波数処理部の第１の例における整形データ列に対する特性値の測定手法の一例の説明図である。
【図１５】  空間周波数処理部の第１の例においてテストパターン画像から得られた基本ＡＣパッチ画像群の各パッチ画像の周期および空間周波数の一例の説明図である。
【図１６】  空間周波数処理部の第１の例においてテストパターン画像の基本ＡＣパッチ画像群から得られたコントラスト値の一例を示すグラフである。
【図１７】  基本ＡＣパッチ画像群における算出された空間周波数伝達特性の一例を示すグラフである。
【図１８】  基本ＡＣパッチ画像群の各パッチ画像から求められた空間周波数伝達特性に基づいて設計された空間周波数処理部における空間周波数伝達率の一例を示すグラフである。
【図１９】  補正量信号のコントラストと増幅率の関係を示すグラフである。
【図２０】  ＤＣ成分信号値と増幅率の関係を示すグラフである。
【図２１】  ＡＣパッチ画像群のパッチ画像の別の例の説明図である。
【図２２】  空間周波数処理部の構成の第２の例を示すブロック図である。
【図２３】  増幅率決定部の構成の一例を示すブロック図である。
【図２４】  空間周波数処理部の第２の例において用いる処理係数を決定する処理の一例を示すフローチャートである。
【図２５】  空間周波数処理部の第２の例において用いるテストパターン画像データのレイアウトの一例を示す説明図である。
【図２６】  空間周波数処理部の第２の例において用いるＡＣ画像群の各パッチ画像の特性値の一例である。
【図２７】  空間周波数処理部の第２の例において用いるＡＣ画像群の各パッチ画像の特性値の一例（続き）である。
【図２８】  反射率の測定の別の具体例の説明図である。
【図２９】  空間周波数処理部の第２の例において測定された測定データ列の並び方向の説明図である。
【図３０】  空間周波数処理部の第２の例における整形データ画像に対する特性値の算出手法の一例の説明図である。
【図３１】  空間周波数処理部の第２の例においてテストパターン画像から得られたＡＣパッチ画像のうちＶ１．１〜Ｖ１．８およびＨ１．１〜Ｈ１．８の各パッチ画像の周期および空間周波数の一例の説明図である。
【図３２】  空間周波数処理部の第２の例においてテストパターン画像から得られたコントラスト値の一例を示すグラフである。
【図３３】  出力装置における入出力コントラストの関係の一例を示すグラフである。
【図３４】  基本ＡＣパッチ画像群の各パッチ画像から求められた空間周波数伝達特性に基づいて設計された空間周波数処理部の第２の例における空間周波数伝達率の一例を示すグラフである。
【図３５】  空間周波数処理部の第２の例において算出された垂直方向の目標空間周波数伝達特性の一例を示すグラフである。
【図３６】  入力信号値が制限範囲を越える場合の画像の出力特性の説明図である。
【図３７】  入力信号値が制限範囲を下回る場合の画像の出力特性の説明図である。
【図３８】  出力部において入力信号値が出力された画像の反射率に反映されない場合の説明図である。
【図３９】  入力信号値の調整処理の概念図である。
【図４０】  入力信号値の制限範囲を考慮した入出力コントラストの関係の一例を示すグラフである。
【図４１】  入力信号コントラストに対する出力画像の実際の平均反射率の一例を示すグラフである。
【図４２】  空間周波数処理部の第２の例における空間周波数伝達率の別の例を示すグラフである。
【図４３】  空間周波数処理部の第２の例における差分特性の一例を示すグラフである。
【図４４】  水平方向および垂直方向の空間周波数を伝達率の関係の説明図である。
【図４５】  空間周波数処理部の第２の例における補助補正部の空間周波数伝達特性に対する各条件における倍率の一例を示すグラフである。
【図４６】  空間周波数処理部の構成の第３の例を示すブロック図である。
【図４７】  空間周波数処理部の構成の第４の例を示すブロック図である。
【図４８】  空間周波数処理部の構成の第５の例における空間周波数分割による目標空間周波数伝達特性の実現の説明図である。
【図４９】  空間周波数処理部の構成の第６の例を示すブロック図である。
【図５０】  本発明を適用した画像編集・出力システムの一例を示す構成図である。
【符号の説明】
１…処理演算部、２…出力部、３…測定部、４…係数算出部、５…テストパターン保持部、１１…色空間変換部、１２…空間周波数処理部、１３…階調変換部、２１…ハイパスフィルタ、２２…ローパスフィルタ、２３，２４…ルックアップテーブル、２５…乗算部、２６…加算部、３１…入力画像信号、３２…補正量信号、３３…ＡＣ成分ゲイン信号、３４…調整補正量信号、３５…ＤＣ成分信号、３６…ＤＣ成分ゲイン信号、３７…出力画像信号、１１１…基本補正部、１１２…補助補正部、１１３…増幅率決定部、１１４…乗算部、１１５…加算部、１１６…乗算部、１２１…１２１…入力画像信号、１２２…基本補正信号、１２３…補助補正基本信号、１２４…補助補正量増幅率信号、１２５…補助補正信号、１２６…出力画像信号、１２７…基本補正量増幅率信号、１３１…ＡＣ成分抽出器、１３２…ＤＣ成分抽出器、１３３…２次元ルックアップテーブル、１３４…ＡＣ成分信号、１３５…ＤＣ成分信号、１３６…２次元ルックアップテーブル、１９１…デジタルフィルタ、１９２…乗算部、１９３…２次元ルックアップテーブル、１９４…加算部、２０１…デジタルフィルタ、２０２…パラメータデータベース。
【図１】
0004035668.tif 000002
【図２】
0004035668.tif 000003
【図３】
0004035668.tif 000004
【図４】
0004035668.tif 000005
【図５】
0004035668.tif 000006
【図６】
0004035668.tif 000007
【図７】
0004035668.tif 000008
【図８】
0004035668.tif 000009
【図９】
0004035668.tif 000010
【図１０】
0004035668.tif 000011
【図１１】
0004035668.tif 000012
【図１２】
0004035668.tif 000013
【図１３】
0004035668.tif 000014
【図１４】
0004035668.tif 000015
【図１５】
0004035668.tif 000016
【図１６】
0004035668.tif 000017
【図１７】
0004035668.tif 000018
【図１８】
0004035668.tif 000019
【図１９】
0004035668.tif 000020
【図２０】
0004035668.tif 000021
【図２１】
0004035668.tif 000022
【図２２】
0004035668.tif 000023
【図２３】
0004035668.tif 000024
【図２４】
0004035668.tif 000025
【図２５】
0004035668.tif 000026
【図２６】
0004035668.tif 000027
【図２７】
0004035668.tif 000028
【図２８】
0004035668.tif 000029
【図２９】
0004035668.tif 000030
【図３０】
0004035668.tif 000031
【図３１】
0004035668.tif 000032
【図３２】
0004035668.tif 000033
【図３３】
0004035668.tif 000034
【図３４】
0004035668.tif 000035
【図３５】
0004035668.tif 000036
【図３６】
0004035668.tif 000037
【図３７】
0004035668.tif 000038
【図３８】
0004035668.tif 000039
【図３９】
0004035668.tif 000040
【図４０】
0004035668.tif 000041
【図４１】
0004035668.tif 000042
【図４２】
0004035668.tif 000043
【図４３】
0004035668.tif 000044
【図４４】
0004035668.tif 000045
【図４５】
0004035668.tif 000046
【図４６】
0004035668.tif 000047
【図４７】
0004035668.tif 000048
【図４８】
0004035668.tif 000049
【図４９】
0004035668.tif 000050
【図５０】
0004035668.tif 000051