(57)【特許請求の範囲】
【請求項１】
画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得手段と、
前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択手段と、
前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復手段とを有し、
前記選択手段は、前記画像形成装置によって形成される画像の空間周波数に対する明るさの特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けるパラメータを選択することを特徴とする画像処理装置。
【請求項２】
前記第１の周波数は、画像の前記鮮鋭性の認識が困難である空間周波数であることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記パラメータは、フィルタの係数であって、前記回復処理は、前記フィルタの畳み込み処理であることを特徴とする請求項１又は請求項２に記載の画像処理装置。
【請求項４】
前記パラメータは、周波数特性であって、前記回復処理は、周波数空間上で前記画像データの周波数特性と前記パラメータが示す周波数特性との積を算出する処理であることを特徴とする請求項１又は請求項２に記載の画像処理装置。
【請求項５】
前記第２の周波数は、前記画像形成装置の出力解像度に対応する周波数であることを特徴とする請求項１乃至４の何れか一項に記載の画像処理装置。
【請求項６】
画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得手段と、
前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択手段と、
前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復手段とを有し、
前記パラメータの周波数特性は、前記画像形成装置によって形成される画像の周波数特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けることを特徴とする画像処理装置。
【請求項７】
前記パラメータの周波数特性が、前記第１の周波数から、前記第２の周波数における周波数帯において、略線形に変化することを特徴とする請求項６に記載の画像処理装置。
【請求項８】
前記パラメータの周波数特性が、前記第１の周波数から、前記第２の周波数における周波数帯において、前記第１の周波数における振幅ゲインを保持することを特徴とする請求項７に記載の画像処理装置。
【請求項９】
前記パラメータの周波数特性は、前記第２の周波数における明るさが、前記第１の周波数における明るさと、略同一となるように決定されることを特徴とする請求項７に記載の画像処理装置。
【請求項１０】
前記パラメータの周波数特性は、回復処理の前後における前記第２の周波数における明るさの差が、前記第１の周波数における明るさの差と、略同一となるように決定されることを特徴とする請求項７に記載の画像処理装置。
【請求項１１】
前記第２の周波数は、前記画像形成装置の出力解像度に対応する周波数であることを特徴とする請求項６乃至１０の何れか一項に記載の画像処理装置。
【請求項１２】
コンピュータを、請求項１乃至請求項１１の何れか一項に記載された画像処理装置の各手段として機能させることを特徴とするプログラム。
【請求項１３】
画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得工程と、
前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択工程と、
前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復工程と
を有し、
前記選択工程は、前記画像形成装置によって形成される画像の空間周波数に対する明るさの特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けるパラメータを選択することを特徴とする画像処理方法。
【請求項１４】
画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得工程と、
前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択工程と、
前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復工程と
を有し、
前記パラメータの周波数特性は、前記画像形成装置によって形成される画像の周波数特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けること
を特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は画像の鮮鋭性を回復する処理に関する。
【背景技術】
【０００２】
インクジェットプリンタなどの画像形成装置より出力される画像は、インクの着弾位置ずれやインクの滲み（メカニカルドットゲイン）、光学的暈け（オプティカルドットゲイン）等により、入力画像に比べ鮮鋭性が低下することが知られている。このとき、出力画像の周波数特性をあらかじめ取得できれば、入力画像をその逆特性を持つフィルタで畳み込み処理（以下、鮮鋭回復処理と呼ぶ）をすることで、鮮鋭低下を補償することができる。
【０００３】
しかしながら、インクジェットプリンタなどの画像形成装置では、入力画像に対して逆特性を持つフィルタで鮮鋭回復処理を行うと、高周波部で明るさの低下が発生することが知られている。
【０００４】
特許文献１では、鮮鋭回復処理の前後における明るさ変化量をあらかじめ計測し、入力画像の各画素に対して計測された変化量で補正することで、明るさの低下を抑制する技術が記載されている。
【０００５】
また、特許文献２では、ディジタル画像中にデルタ関数を埋め込み、鮮明化、印刷、及び、走査のような幾つかの中間処理段階を通じて処理されたディジタル画像からデルタ関数を引き出す。取り出されたデルタ関数からＭＴＦ（周波数特性）を推定して補正する方法が記載されている。中間処理装置のノイズ特性が既知の場合、鮮明化（回復）処理を行うウィーナーフィルタの特性は、ＭＴＦの二乗及びノイズの二乗の和でＭＴＦを割算することで鮮明化処理の信号対ノイズ比を最大化する技術が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１１−２４０４９
【特許文献２】特開２００１−１８８９０３
【発明の概要】
【発明が解決しようとする課題】
【０００７】
特許文献１に記載の技術では、画像の特徴量（周波数成分ごとの振幅値、平均輝度（明度）値）に応じた輝度変化量をあらかじめ測定し、保持しておく必要があった。さらに、鮮鋭回復処理時には入力画像の画素ごとに、前記特徴量（周波数成分ごとの振幅値、平均輝度値）を算出しなければならず、保持コスト、回復処理の測定コスト、計算コストが膨大であった。
特許文献２に記載の技術は、ノイズ量に応じて周波数特性を抑制したものであり、明るさ低下には対応できなかった。
そこで本発明では、明度低下が知覚されにくい鮮鋭回復処理を従来に比べて簡易に行うことを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明に係る画像処理装置は、画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得手段と、前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択手段と、前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復手段とを有し、前記選択手段は、前記画像形成装置によって形成される画像の空間周波数に対する明るさの特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けるパラメータを選択することを特徴とする。
【０００９】
また、本発明に係る画像処理装置は、画像形成装置が形成する画像の鮮鋭性が低下する特性に対応する出力条件を取得する取得手段と、前記出力条件に基づき、画像の鮮鋭性を回復するためのパラメータを選択する選択手段と、前記選択されたパラメータを用いて鮮鋭性の回復処理を画像データに施す回復手段とを有し、前記パラメータの周波数特性は、前記画像形成装置によって形成される画像の周波数特性が、第１の周波数から前記画像形成装置の出力解像度以下の解像度に対応する周波数である第２の周波数まで、変曲点または不連続点なく、一定もしくは減少し続けることを特徴とする。
【発明の効果】
【００１０】
本発明によれば、明度低下が知覚されにくい鮮鋭回復処理を従来に比べて簡易に行うことが可能となる。
【図面の簡単な説明】
【００１１】
【図１】視覚特性の一例を示す図である。
【図２】画像形成装置の周波数特性と、その逆特性の一例を示す図である。
【図３】逆特性で回復処理をおこなった場合の明るさの特性の一例を示す図である。
【図４】空間位置に対して周波数が滑らかに変化する入力画像に対して、逆特性で回復処理を行った場合の明るさ低下の一例を示す図である。
【図５】周波数特性Ｆａ（ｕ）の一例を示す図である。
【図６】周波数特性Ｆａ（ｕ）で回復処理を行った場合の明るさの特性Ｌａ（ｕ）の一例を示す図である。
【図７】図４（ａ）に示す入力画像に対して、周波数特性Ｆａ（ｕ）で回復処理した場合の明るさ低下の一例を示す図である。
【図８】周波数特性Ｆｂ（ｕ）と、周波数特性Ｆｂ（ｕ）で回復処理した場合の明るさの特性Ｌｂ（ｕ）の一例を示す図である。
【図９】実施例１における画像処理装置の構成の一例を示す図である。
【図１０】実施例１における処理フローの一例を示す図である。
【図１１】連携ＬＵＴの一例を示す図である。
【図１２】回復フィルタと連携ＬＵＴの作成フローの一例を示す図である。
【図１３】鮮鋭性計測チャートの一例を示す図である。
【図１４】「フィルタ番号」列が空欄である連携ＬＵＴの一例を示す図である。
【図１５】周波数ｕ＿ａ，ｕ＿ｂ，ｕ＿ｃにおけるフィルタの周波数特性に対する明るさの変化の一例を示す図である。
【図１６】周波数特性Ｆｃ（ｕ）の一例を示す図である。
【図１７】周波数特性Ｆｃ（ｕ）で回復処理を行った場合の明るさの特性Ｌｃ（ｕ）を示す図である。
【図１８】周波数ｕ１，ｕ２における回復前後の明るさの差分が略一致するようにした場合の明るさの特性の例を示す図である。
【図１９】実施例２における回復フィルタの作成フローの一例を示す図である。
【図２０】実施例２における変形例を示す図である。
【発明を実施するための形態】
【００１２】
以下、本発明の実施形態について、図面を参照して説明する。なお、以下の実施形態は本発明を限定するものではなく、また、本実施形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。なお、同一の構成については、同じ符号を付して説明する。
【実施例１】
【００１３】
図１に公知の方法（Ｄｏｏｌｅｙの近似式）で計算された視覚特性を示す。ただし、観察距離は３００ｍｍとする。図１に示すように、空間周波数の上昇に伴い人間の視覚の感度特性（以下、視覚特性）は低下する。例えば４サイクル／ｍｍを超える空間周波数域（以下、周波数域）において視覚特性は０．２５以下になり、鮮鋭性の低下が発生したとしても、その認識は難しい。以下では、鮮鋭性低下の認識が難しくなる空間周波数（例えば４サイクル／ｍｍ）を「視覚限界周波数」と呼ぶことにする。また、以下では、空間周波数を単に「周波数」と記載する場合がある。
【００１４】
また、以下の実施例においては、ＣＩＥ  Ｌａｂに含まれる明度Ｌ＊の周波数に対する特性を、明るさの特性として説明する。しかしながら、たとえば、ＸＹＺ表色系におけるＹの特性を明るさの特性として用いてもよい。あるいは、紙面上の光学濃度の特性を明るさの特性として用いてもよい。
（鮮鋭回復処理による明るさ低下）
出力画像の周波数特性Ｐ（ｕ）を取得できれば、Ｐ（ｕ）の逆特性Ｒ（ｕ）＝１／Ｐ（ｕ）を持つフィルタで、入力画像に対して鮮鋭回復処理を行うことで、出力画像の鮮鋭低下を補償ことができる。たとえば、出力画像の周波数特性が図２の曲線Ｐ（ｕ）に示す特性とすると、逆特性１／Ｐ（ｕ）は、図２の曲線Ｒ（ｕ）に示す特性となる。このとき、逆特性Ｒ（ｕ）を逆フーリエ変換して得られるフィルタを用いて鮮鋭回復処理した画像を出力すれば、鮮鋭低下を抑制した画像を得ることができる。
【００１５】
しかしながら、画像形成装置の入力画像に対して、逆特性を持つフィルタで鮮鋭回復処理を行うと、高周波領域において平均の明るさが低下する。逆特性Ｒ（ｕ）で回復処理を行った場合の明るさ低下の例を、図３の曲線Ｌ（ｕ）で示す。Ｌ（ｕ）は、周波数ｕ、平均明度Ｌａｖｅの振幅を持つ正弦波を入力画像した場合の、周波数ｕに対する出力画像の平均の明るさを示している。
【００１６】
図３に示すように、逆特性Ｒ（ｕ）で鮮鋭回復処理を行うと、周波数ｕに応じて明るさが低下する。例えば、空間位置に応じて周波数が滑らかに変化する入力画像に逆特性Ｒ（ｕ）で鮮鋭回復処理を行った場合、明るさ低下が、空間位置に応じて発生する。図４を用いてより具体的に説明する。図４（ａ）は、正弦波の空間周波数ｕが、空間位置ｘ［ｍｍ］に対して線形に変化する画像の一例を示している。つまり、画素の位置が空間周波数に対応し、図の左側は周波数が低く、図の右側に行くほど周波数が高い。このような画像は下記式（１）を用いて生成することができる。
【００１７】
Ｉ（ｘ，ｙ）＝Ｌｈ×ｓｉｎ（π×ｘ２／ｗ）＋Ｌａｖｅ      （１）
ただし、０≦ｘ＜ｗとし、Ｉ（ｘ，ｙ）は座標（ｘ，ｙ）における画素値、ｗ［ｍｍ］は画像の横幅、Ｌｈは正弦波の振幅、Ｌａｖｅは平均明度、πは円周率とする。このとき、位置ｘと空間周波数ｕとの関係は、比例定数Ａを用いると、ｕ＝Ａ×ｘであらわされる。ここで、ｕ＿ｍａｘを出力する画像形成装置が表現できる最高の解像度（ｄｐｉ）とし、位置ｘと空間周波数ｕの比例定数Ａ＝（ｕ＿ｍａｘ／（２×２５．４×ｗ））とすると、図４（ａ）に示すような、低周波から画像形成装置が表現できる最高の解像度まで遷移する画像となる。この図４（ａ）を入力画像とし、図２（ｂ）に示す逆特性Ｒ（ｕ）で回復処理を行うと、出力画像の平均の明るさの特性は、Ｌ（ｕ）＝Ｌ（Ａ×ｘ）となる。したがって、Ｌ（ｕ）が図３に示す特性の場合、ｘが大きくなるほど明るさ低下が生じる。よって、出力画像は図４（ｂ）のようになる。
【００１８】
図４（ｂ）に示すような明るさの低下は、鮮鋭回復処理によって高周波成分の振幅が拡大されることに起因する。そのため、回復処理に用いるフィルタの高周波領域での周波数特性を抑制すれば、明るさ低下が知覚されにくくなる。たとえば、鮮鋭回復処理に用いるフィルタの周波数特性Ｆａ（ｕ）の高周波部を１．０に抑制すれば、高周波数部の振幅が拡大されず、鮮鋭回復処理に起因する明るさ低下は生じない。しかしながら、Ｆａ（ｕ）の高周波部を１．０に抑制すると、高周波部の鮮鋭回復処理の効果も失われる。そのため、鮮鋭性低下が知覚困難となる視覚限界周波数よりも低周波な領域では、鮮鋭回復の効果を優先して鮮鋭回復特性を抑制せず、視覚限界周波数よりも高周波な領域では、明るさ低下を低減するために鮮鋭回復特性を抑制することが好ましい。すなわち、明るさの低下は、回復処理によって振幅が拡大されることに起因するため、回復処理に用いるフィルタの高周波領域での周波数特性を抑制することで、より低コストに明るさ低下を抑制することができる。従って、視覚限界周波数をｕ＿ｖｉｓとしたとき、たとえば、図５に示すような周波数特性Ｆａ（ｕ）で鮮鋭回復処理を行えばよい。Ｆａ（ｕ）は、以下の式（２）で定義される。
【００１９】
Ｆａ（ｕ）  ＝  Ｒ（ｕ）      （０≦ｕ＜ｕ＿ｖｉｓ）
＝１．０        （ｕ＿ｖｉｓ≦ｕ＜ｕ＿ｍａｘ）      （２）
周波数特性Ｆａ（ｕ）の鮮鋭回復フィルタを用いることで、視覚的に重要な領域において、フィルタの効果を残しつつ、明るさ低下を抑制することができる。即ち、周波数特性Ｆａ（ｕ）を逆フーリエ変換してフィルタを設計すればよく、画像の特徴量（周波数成分ごとの振幅値、平均輝度（明度）値）に応じた輝度変化量を測定する必要がない。また、鮮鋭回復処理としては、通常の鮮鋭回復フィルタと同様の処理のみでよく、入力画素ごとの前記特徴量（周波数成分ごとの振幅値、平均輝度（明度）値）を算出する必要もない。したがって、測定コスト、計算コストの増加なく、鮮鋭回復処理における明るさ低下を抑制することができる。
【００２０】
ところが、周波数特性Ｆａ（ｕ）で鮮鋭回復処理を行うと、明るさの特性は図６に示すＬａ（ｕ）のような下に凸の不連続な形状となる。そのため、図４（ａ）のような空間位置に対して周波数が滑らかに変化する入力画像に周波数特性Ｆａ（ｕ）で鮮鋭回復処理を行うと、図７に示すような帯状の不自然な明るさ変動が知覚される。このような不自然な明るさ変動は、鮮鋭回復処理後の出力画像の画質を著しく損なう。よって、出力画像において、明るさ変動が知覚されにくくなるような制御が必要となる。
【００２１】
本実施例では、鮮鋭回復処理に用いるフィルタの周波数特性を所定の周波数ｕ１から振幅ゲインを保持することで、鮮鋭回復処理による明るさの変動を知覚されにくくする。
【００２２】
例えば式（３）に示すように、周波数ｕ１までは、鮮鋭回復の効果を優先して鮮鋭回復特性（振幅ゲイン）を抑制せず、周波数ｕ１よりも高周波な領域の周波数帯では、明るさ低下を低減するために鮮鋭回復特性（振幅ゲイン）をＲ（ｕ１）に保持する。こうすれば、明るさ変動が知覚されにくい鮮鋭回復処理が実現できる。
【００２３】
Ｆｂ（ｕ）  ＝  Ｒ（ｕ）            （０≦ｕ＜ｕ１）
＝  Ｒ（ｕ１）          （ｕ１≦ｕ）      （３）
上記を図を用いてより詳細に説明する。図８（ａ）の曲線Ｆｂ（ｕ）は、式（３）で定義されるフィルタ特性Ｆｂ（ｕ）である。Ｆｂ（ｕ）は、周波数ｕ１〜ｕ２において一定である。実験の結果、フィルタ特性Ｆｂ（ｕ）による鮮鋭回復処理後の明るさＬｂ（ｕ）は、周波数ｕ１〜ｕ２において線形に低下することがわかった。このときの明るさＬｂ（ｕ）を、図８（ｂ）に示す。
【００２４】
一般的な画像形成装置における出力画像の周波数特性は、図２の周波数特性Ｐ（ｕ）に示すように、周波数が大きくなると振幅ゲインは小さくなるため、逆特性Ｒ（ｕ）は曲線２０２に示すように大きくなり、明るさの低下量が増加する。そこで、図８（ａ）に示すように周波数ｕ１から鮮鋭回復フィルタの周波数特性（振幅ゲイン）を保持するようにすると、周波数ｕ１よりも高周波な領域において、明るさの低下が抑制される。すなわち、鮮鋭回復処理の周波数特性として前記Ｆｂ（ｕ）を用いれば、周波数ｕ１以降の明るさの低下を抑制することができる。このときの線形回復処理後の明るさＬｂ（ｕ）は、図８（ｂ）に示すように周波数ｕ１〜ｕ２において、周波数に対してほぼ線形に低下する。そのため、鮮鋭回復処理を行っても、図７に示すような不自然な明るさ変動が生じない。ここで、周波数ｕ１を視覚限界周波数ｕ＿ｖｉｓとし、周波数ｕ２を画像形成装置の出力解像度ｕ＿ｍａｘとすれば、フィルタ特性Ｆｂ（ｕ）で鮮鋭回復処理を行うことで、明るさ変動が知覚されにくい鮮鋭回復処理を実現できる。
【００２５】
（装置の構成と処理フロー）
図９は、本実施例における画像処理装置１００の構成例を示すブロック図である。また、図１０は本実施例における画像処理装置１００の処理を示すフローチャートである。以下、図９のブロック図と図１０のフローチャートを用いて、画像処理装置１００における処理の概要を説明する。
【００２６】
まず、回復処理部１０７は、データ入出力（ｉ／ｏ）部１０４を介して、情報処理装置１５０などから画像形成対象の画像データｉを取得する。さらに、取得した画像データｉを例えばＲＡＭなどのメモリ部１０１に格納する（Ｓ１００１）。
【００２７】
次に、条件取得部１０３は、ユーザインタフェース（ＵＩ）部１０２を介して、あるいは、情報処理装置１５０から、画像形成部１０８における出力条件Ｏを取得する（Ｓ１００２）。出力条件Ｏは、パス数、キャリッジ速度、印刷方向（双方向または単方向）、ハーフトーン設定、印刷ヘッドと記録媒体の間の距離、クリアインク使用の有無、色設定、サイレント設定の有無、記録媒体の種類などを含む。なお、クリアインク（以下、ＣＬインク）とは、無色透明または実質的に無色透明の顔料を含むインクである。
【００２８】
次に、フィルタ選択部１０６は、パラメータ保持部１０５にアクセスして、出力条件Ｏに対応する鮮鋭回復フィルタｒを取得する（Ｓ１００３）。鮮鋭回復フィルタの詳細、および出力条件Ｏと鮮鋭回復フィルタｒの対応については後述する。
【００２９】
次に、鮮鋭回復処理部１０７は、入力画像データｉに、鮮鋭回復フィルタｒによる鮮鋭回復処理を実行し、鮮鋭回復処理後の画像データｉ´をメモリ部１０１に格納する（Ｓ１００４）。そして、画像形成部１０８は、出力条件Ｏに基づき、鮮鋭回復処理後の画像データｉ´が表す画像を記録媒体２０８上に形成する（Ｓ１００５）。
【００３０】
（画像形成部）
画像形成部１０８の色分解部２０１は、色分解テーブルを参照して、例えばＲＧＢ形式の画像データｉ´を画像形成部１０８が備えるインク色（ＣＬインクを含む）に対応する色材量データに色分解する。ハーフトーン（ＨＴ）処理部２０２は、色分解部２０１が出力する色材量データを量子化してドットパターンに変換し、記録データ（以下、ＨＴデータ）として出力する。
【００３１】
画像形成部１０８は、熱転写方式、インクジェット方式などの記録装置であり、記録ヘッド２０５を記録媒体２０８に対して相対的に縦横に移動して、バンド単位に入力されるＨＴデータが示す画像を記録媒体２０８上に形成する。その際、インク色選択部２０３は、入力されるＨＴデータに対応するインク色を、記録ヘッド２０５に搭載されたインク色の中から選択する。
【００３２】
記録ヘッド２０５は、一つ以上の記録素子（インクジェット方式であればノズル）を有する。記録ヘッド２０５の相対的な縦横の移動は、ヘッド制御部２０４が移動部２０６を制御して記録ヘッド２０５をＸ方向（主走査方向）に移動し、ヘッド制御部２０４が搬送部２０７を制御して記録媒体２０８をＹ方向（副走査方向）に搬送することで実現される。
【００３３】
（パラメータ保持部）
パラメータ保持部１０５は、例えばＥＥＰＲＯＭやフラッシュメモリなどの不揮発性メモリであり、前述した出力条件の組み合わせと鮮鋭回復フィルタとの対応を示すルックアップテーブル（以下、対応ＬＵＴ）を保持する。図１１により対応ＬＵＴの一例を示す。なお、出力条件に対応する鮮鋭回復フィルタの作成方法は後述する。
【００３４】
（情報処理装置）
画像処理装置１００は、パーソナルコンピュータなどの情報処理装置１０５に接続される。画像処理装置１００と情報処理装置１５０の接続には、ＵＳＢなどのシリアルバス、有線または無線ネットワークなどを利用することができる。
【００３５】
情報処理装置１５０のＣＰＵ１５１は、ＲＡＭ１５４をワークメモリとして、ＲＯＭ１５３や記憶部１５５に格納されたＯＳや各種プログラムを実行し、情報処理装置１５０の動作を制御する。ＨＤＤやＳＳＤなどの記憶部１５５には、プリンタドライバや、対応ＬＵＴおよび鮮鋭回復フィルタを作成する作成プログラムなどが格納されている。プリンタドライバを実行するＣＰＵ１５１により、ＵＳＢなどの汎用インタフェイス１５２や有線または無線ネットワークを介して、画像データが画像処理装置１００へ供給される。また、作成プログラムをＣＰＵ１５１で実行して対応ＬＵＴや対応回復フィルタを作成し、汎用インタフェイス１５２や有線または無線ネットワークを介して、パラメータ保持部１０５に保持される。
【００３６】
また、画像処理装置１００がスキャナを備える複合機の場合など、画像処理装置１００の図示しないＣＰＵに対応ＬＵＴおよび鮮鋭回復フィルタを作成するための作成プログラムを実行させて、対応ＬＵＴおよび鮮鋭回復フィルタの作成を行ってもよい。
【００３７】
（鮮鋭回復フィルタと対応ＬＵＴの作成）
図１２に示すフローチャートを用いて、鮮鋭回復フィルタと連携対応ＬＵＴの作成手順を説明する。
【００３８】
ＣＰＵ１５１は、カウンタｉを「１」に初期化し（Ｓ１２０１）、カウンタｉに対応する鮮鋭回復フィルタを生成する出力条件Ｏｉを取得する（Ｓ１２０２）。たとえば、図１１に出力条件Ｏ１として示す、記録媒体「Ａ」、パス数「３２」、キャリッジ速度「遅」、印刷方向「単方向」、ＨＴ処理「誤差拡散」、ヘッド距離「短」、クリアインク「無」、色設定「カラー」が取得される。
【００３９】
次に、ＣＰＵ１５１は、画像形成部１０８を制御して、出力条件Ｏｉにより鮮鋭性の計測チャートを出力する（Ｓ１２０３）。なお、計測チャートの出力時は回復処理を行わない。図１３に計測チャートの一例を示す。計測チャートは、周波数が異なる複数の正弦波パターンと、左下の均一パターン（例えば白べたと黒べた）を含む画像チャートである。
【００４０】
次に、ＣＰＵ１５１は、汎用インタフェイス１５２に接続された不図示の測定装置を用いて、周波数特性の取得に必要な情報を計測チャートから取得する。測定装置には、例えば、スキャナ、ディジタルカメラ、顕微鏡、マイクロデンシトメータなどを用いることができる。ＣＰＵ１５１は、取得した情報に基づき、出力条件Ｏｉに対する画像形成装置の周波数応答値Ｐ（ｕ）を算出する（Ｓ１２０４）。
【００４１】
周波数応答値Ｐ（ｕ）は、例えば下式を用いて算出される光学伝達関数（ＭＴＦ）を用いることができる。
【００４２】
Ｐ（ｕ）  ＝  Ｃ（ｕ）／Ｃ’      （４）
ここで、ｕは正弦波の周波数、
Ｃ（ｕ）＝｛Ｍａｘ（ｕ）−Ｍｉｎ（ｕ）｝／
｛Ｍａｘ（ｕ）＋Ｍｉｎ（ｕ）｝、
Ｃ’＝（Ｗｈｉｔｅ−Ｂｌａｃｋ）／（Ｗｈｉｔｅ＋Ｂｌａｃｋ）、
Ｍａｘ（ｕ）は周波数ｕで変化する正弦波パターンの最大明度、
Ｍｉｎ（ｕ）は周波数ｕで変化する正弦波パターンの最小明度、
ＷｈｉｔｅとＢｌａｃｋはそれぞれ均一パターンの明度。
【００４３】
勿論、光学伝達関数の算出は式（４）に限定されず、例えば下式を用いてもよい。
【００４４】
Ｐ（ｕ）＝｛Ｍａｘ（ｕ）−Ｍｉｎ（ｕ）｝／（Ｗｈｉｔｅ−Ｂｌａｃｋ）
（４’）
なお、上記では、Ｍａｘ（ｕ）とＭｉｎ（ｕ）、ＷｈｉｔｅとＢｌａｃｋを明度として周波数応答値Ｐ（ｕ）を算出しているが、例えば輝度や濃度、測定装置のデバイスＲＧＢ値等を用いて算出してもよい。
【００４５】
また、Ｓ１２０３で出力する計測チャートとして、図１３に示すような正弦波パターンではなく、矩形波パターンを用いて周波数特性Ｐ（ｕ）を取得してもよい。その場合、矩形波パターンに対して式（４）を適用することにより算出されるコントラスト伝達関数（ＣＴＦ）の値を周波数特性Ｐ（ｕ）として用いる。もしくは、ＣＴＦ値を公知のコルトマン補正式を用いて変換したＭＴＦ値を周波数特性Ｐ（ｕ）に用いてもよい。
【００４６】
次に、ＣＰＵ１５１は、Ｓ１２０４で取得した空間周波数特性Ｐ（ｕ）に基づき、鮮鋭回復フィルタの周波数特性Ｆｂ（ｕ）を算出する。このとき、前述の式（３）に従って鮮鋭回復フィルタを生成するため、まずＣＰＵ１５１は、ｕとｕ１との大小を比較する（Ｓ１２０５）。そして、Ｓ１２０５における比較の結果、ｕ＜ｕ１であれば、ＣＰＵ１５１はＦｂ（ｕ）＝Ｒ（ｕ）＝１／Ｐ（ｕ）より周波数特性Ｆｂ（ｕ）を算出する（Ｓ１２０６）。また、Ｓ１２０５における比較の結果、ｕ≧ｕ１であれば、ＣＰＵ１５１はＦｂ（ｕ）＝Ｒ（ｕ１）＝１／Ｐ（ｕ１）より周波数特性Ｆｂ（ｕ）を算出する（Ｓ１２０７）。
【００４７】
次に、ＣＰＵ１５１は、画像形成装置の出力解像度まで空間周波数特性Ｐ（ｕ）を取得したかを判断するため、解像度ｕとｕ２との大小を比較する（Ｓ１２０８）。比較の結果、ｕ＜ｕ２であれば、出力解像度まで取得していないと判断し、ｕを更新し（Ｓ１２１０）、鮮鋭回復フィルタの周波数特性Ｆｂ（ｕ）の算出を続ける。比較の結果、ｕ≧ｕ２であれば、出力条件Ｏｉにおける鮮鋭回復フィルタの周波数特性Ｆｂ（ｕ）の算出を完了したと判断する。そして、周波数特性Ｆｂ（ｕ）を逆フーリエ変換することで、鮮鋭回復フィルタの係数ｆｂ（ｘ）を算出する（Ｓ１２０９）。
【００４８】
次に、ＣＰＵ１５１は、鮮鋭回復フィルタ係数ｆｂ（ｘ）を、出力条件Ｏｉと対応付けてパラメータ保持部１０５に格納する（Ｓ１２１１）。このとき、出力条件Ｏｉと鮮鋭回復フィルタが対応するように、対応ＬＵＴを作成する。
【００４９】
たとえば、図１４に示すように、フィルタ番号列が空欄である対応ＬＵＴを、あらかじめパラメータ保持部１０５に保持しておく。そして、対応ＬＵＴに格納されたフィルタ番号のうち、最も大きなフィルタ番号に「１」を加えたフィルタ番号を新規のフィルタ番号ｒとして、対応する出力条件Ｏｉの行に書き込めばよい。このとき、算出した鮮鋭回復フィルタの係数は、フィルタ番号ｒに関連付けてパラメータ保持部１０５に格納される。
【００５０】
次に、ＣＰＵ１５１は、全ての出力条件で鮮鋭回復フィルタの周波数特性を取得したかを判定する（Ｓ１２１２）。未所得の出力条件がある場合、ＣＰＵ１５１は、カウンタｉをインクリメントし（Ｓ１２１３）、処理をステップＳ１２０３に戻す。これにより、次の出力条件へのフィルタ番号の割り当てが行われる。また、全ての出力条件で鮮鋭回復フィルタの周波数特性を取得した場合、ＣＰＵ１５１は処理を終了する。
【００５１】
なお、上記では、ＣＰＵ１５１がパラメータ保持部１０５の対応ＬＵＴにアクセスしてフィルタ番号の登録処理を行う例を説明した。しかし、ＣＰＵ１５１がパラメータ保持部１０５、記憶部１５５または図示しないサーバ装置などからフィルタ番号が未登録の対応ＬＵＴを取得し、フィルタ番号の登録処理を行った後、登録処理後の対応ＬＵＴをパラメータ保持部１０５に格納するようにしても良い。
【００５２】
また、必ずしも全出力条件に対して固有の鮮鋭回復フィルタを用意する必要はなく、図１１に示すように、複数の出力条件に対して、同一のフィルタを割り当ててもよい。
【００５３】
以上、説明したように本実施例によれば、周波数特性Ｆｂ（ｕ）は、周波数ｕ１までは、鮮鋭回復の効果を優先して鮮鋭回復の周波数特性（振幅ゲイン）を抑制しない。一方、周波数ｕ１よりも高周波な領域では、明るさ低下を低減するために鮮鋭回復の周波数特性（振幅ゲイン）をＲ（ｕ１）に抑制する。このような周波数特性Ｆｂ（ｕ）を逆フーリエ変換することによって得られるフィルタ係数ｆｂを用いることで、逆特性Ｒ（ｕ）を用いた鮮鋭回復処理に比べ、周波数ｕ１以降の明るさの低下を抑制することができる。また、このときの鮮鋭回復処理後の明るさＬｂ（ｕ）は、周波数に対してほぼ線形に低下する。そのため、鮮鋭回復処理を行っても不自然な明るさ変動が生じにくく、明るさ変動が知覚されにくい鮮鋭回復処理を実現できる。
【００５４】
なお、本実施例では鮮鋭回復フィルタにて鮮鋭回復処理を行う例を示したが、これに限らず、例えば、鮮鋭回復フィルタの代りに周波数特性Ｆｂ（ｕ）をパラメータ保持部に格納し、鮮鋭回復処理に用いるようにしてもよい。その場合の鮮鋭回復処理は、まず、入力画像ｉに対してフーリエ変換を行い、入力画像ｉの周波数特性Ｉ（ｕ）を得る。次に、周波数空間上で、各周波数ｕにおいて入力画像ｉの周波数特性Ｉ（ｕ）と鮮鋭回復の周波数特性Ｆｂ（ｕ）の積を求め、鮮鋭回復処理後の周波数特性Ｉ´を得る。そして、鮮鋭回復処理後の周波数特性Ｉ´を逆フーリエ変換することで鮮鋭回復処理後の画像ｉ´を得ればよい。
【００５５】
また、上記実施例では、所定の周波数ｕ１に着目し、所定の周波数ｕ１より高周波な領域で振幅ゲインを保持する例を示した。しかしながら、所定の周波数ではなく、所定の応答値（振幅ゲイン）に着目し、所定の応答値（振幅ゲイン）Ｆ＿ｌｉｍとなる周波数ｕ＿ｌｉｍを算出し、周波数ｕ＿ｌｉｍより高周波な領域で、周波数特性を保持するようにしてもよい。あるいは、所定の応答値（振幅ゲイン）Ｆ＿ｌｉｍ以下となるように制限するようにしても良い。通常、出力画像の周波数特性Ｐ（ｕ）は周波数が高くなるに従い減衰するので、上述したように所定の応答値（振幅ゲイン）Ｆ＿ｌｉｍ以下となるように制限すると、回復処理の周波数特性は、所定の周波数より高周波な領域で振幅ゲインを保持する場合と同様の特性となる。
【実施例２】
【００５６】
以下、実施例２の画像処理および情報処理を説明する。本実施例においては、鮮鋭回復処理に用いるフィルタの周波数特性Ｆｃ（ｕ）が、所定の周波数領域ｕ１〜ｕ２において、目標の周波数特性に対して略線形に減少するように抑制することで、鮮鋭回復処理による明るさ変動を知覚されにくくする。
【００５７】
なお、本実施例において、上記実施例１と略同様の構成については、同一の符号を付して、その詳細な説明を省略する。
【００５８】
本願の発明者は、ある周波数の振幅画像の鮮鋭回復処理後の明るさ（１周期分の平均明度）の低下量は、鮮鋭回復フィルタの振幅ゲインに対して略線形に変化することを導き出した。この周波数と明るさの関係に基づき、各周波数における入力平均値（明るさ）を一定にしたときの所定周波数ｕ１における鮮鋭回復処理後の明るさと、出力解像度に対応する周波数ｕ２おける鮮鋭回復処理後の明るさが略一致する振幅ゲインを求める。さらに、周波数ｕ１〜ｕ２において、略線形に減少するように（例えばｕ１−ｕ２間の振幅ゲインを直線で結ぶ等で）鮮鋭回復フィルタの周波数特性を設計することで、明るさ変動が知覚されにくい鮮鋭回復処理を実現する。
【００５９】
以下、図１５、図１６を用いてより詳細に説明する。図１５（ａ）に示す直線ｕ＿ａ，ｕ＿ｂ，ｕ＿ｃは、それぞれ周波数ｕ＝ｕ＿ａ，ｕ＿ｂ，ｕ＿ｃの矩形波を入力した場合の、鮮鋭回復処理の振幅ゲインにおける出力画像の平均の明るさを示している。ただし、ｕ＿ａ＜ｕ＿ｂ＜ｕ＿ｃとする。図１５（ａ）に示す通り、フィルタの振幅ゲインに対して、出力画像の明るさＬ＊は略線形に変化する。一方、出力画像の明るさがＬ１となる振幅ゲインを周波数でプロットしたものを同図（ｂ）に示す。このようにある周波数以上の領域においては、同一明るさとなるゲインは略線形に変化する。よって、所定周波数以上の領域においては、振幅ゲインをコントロールすることで鮮鋭回復処理後の明るさをコントロールすることができる。
【００６０】
図１６の曲線は、本実施例における鮮鋭回復フィルタの周波数特性Ｆｃ（ｕ）の一例を示している。ここで、周波数ｕ１において、振幅ゲインβ＝Ｒ（ｕ１）で鮮鋭回復処理した場合の出力画像の平均明るさをＬ１´とする。さらに、周波数ｕ１より高い周波数ｕ２において、振幅ゲインβ＝Ｒ（ｕ１）で鮮鋭回復処理した場合の出力画像の平均明るさをＬ２´とする。また、回復処理を行わずに出力した場合（振幅ゲインα＝１．０）の周波数ｕ１における出力画像の平均明るさをＬ１、周波数ｕ２における出力画像の平均明るさをＬ２とする。このとき、Ｌ２´´がＬ１´と略一致する振幅ゲインχは、以下の式（５）によって求めることができる。ただし、α＜χ＜βとする。
【００６１】
χ＝α＋（β−α）×（Ｌ１´−Ｌ２）／（Ｌ２´−Ｌ２）      （５）
但し、Ｌ２≧Ｌ１´＞Ｌ２´
このときの鮮鋭回復処理後の出力画像の平均明るさの周波数特性Ｌｃ（ｕ）の一例を図１７に示す。このように、（ｕ１，β）、（ｕ２，χ）を直線で結んだ周波数特性Ｆｃ（ｕ）で鮮鋭回復処理することで、図１７の実線で示すように、周波数領域ｕ１〜ｕ２において出力画像の平均明るさＬｃ（ｕ）は保持される。その結果、明るさ変動が知覚されにくくなる。このようなフィルタの周波数特性Ｆｃ（ｕ）は、たとえば、以下の式（６）で定義される。
【００６２】
Ｆｃ（ｕ）＝Ｒ（ｕ）          （０≦ｕ＜ｕ１）
＝β−（β−χ）×（ｕ−ｕ１）／（ｕ２−ｕ１）  （ｕ１≦ｕ）
（６）
なお、式（５）では、Ｌ１´とＬ２´´が略一致するように振幅ゲインχを定めた。しかし、図１８に実線で示すように、Ｌ１とＬ１´の差分ΔＬ１を取得し、Ｌ２とＬ２´´の差分ΔＬ２が、ΔＬ１と略一致するようにＬ２´´を定め、線形補間により振幅ゲインχを算出してもよい。
【００６３】
また、上記の説明では、図１６に示すように周波数領域ｕ１〜ｕ２における振幅ゲインをβ＝Ｒ（ｕ１）としてＬ１´、Ｌ２´を求め、次に、周波数ｕ２の振幅ゲインα＝１．０としてＬ２を求めた。しかし、前述の通り、α≠βであれば、式（５）（６）により、Ｆｃ（ｕ）を算出できる。たとえば、α＝０．８としてＦｃ（ｕ）を求めてもよい。αを１より小さくすると、出力画像の明るさの低下を抑えることができるので、黒つぶれを起こしやすいメディア（ドットのにじみが大きい場合）に有効である。なお、α＝１．０の場合は、鮮鋭性回復前の周波数特性を測定するために使用するチャートと同一のチャートを用いることができるため、測定のコストの面から、α＝１．０とすることが望ましい。また、鮮鋭性回復後の画素値が入力画素値の範囲を超えてしまう（以下、この状態をクリッピングが発生と称する）と、図１５に示すような線形性が崩れるため、測定チャートの平均値は入力画素値の範囲の中間レベル付近で振幅ゲインβはクリッピングが発生しない範囲とすることが望ましい。
【００６４】
装置の構成と回復処理の処理フローは、それぞれ図９ならびに図１０で示した実施例１の構成と処理フローと同様であるため、説明を省略する。
【００６５】
以下、図１９を用いて、フィルタの作成手順を説明する。なお、図１２に示した実施例１における処理と同等の処理については、説明を省略する。
【００６６】
ＣＰＵ１５１は画像形成部１０８を制御し、たとえば図１３に示すような周波数特性計測チャートを、出力条件Ｏｉにより出力する。さらに。不図示の測定装置を用いて、周波数特性の取得に必要な情報を計測チャートから取得し、たとえば式（４）´を用いて、出力条件Ｏｉに対する画像形成装置の周波数応答値Ｐ（ｕ）を算出する（Ｓ１９０１）。なお、計測チャートの出力時は鮮鋭回復処理を行わない。
【００６７】
次に、ＣＰＵ１５１は、汎用インタフェイス１５２に接続された不図示の測定装置を用いて、周波数ｕ１およびｕ２における、鮮鋭回復処理を行わない場合の明るさ低下量（図１６、図１７のＬ１，Ｌ２）を取得する（Ｓ１９０２）。たとえば、マイクロデンシトメータを用い、Ｓ１９０１で出力した周波数特性計測チャートのうち、周波数ｕ１およびｕ２光学濃度を周波数毎に算出し明度Ｌに変換すればよい。
【００６８】
次に、ＣＰＵ１５１は、Ｓ１９０１で取得された画像形成装置の周波数応答値Ｐ（ｕ）に基づき、振幅ゲインβ＝Ｒ（ｕ１）＝１／Ｐ（ｕ１）を算出する。さらに、式（３）に従い、周波数ｕ１より高周波な領域で周波数特性がβで一定となる周波数特性Ｆｂ（ｕ）を求め、逆フーリエ変換により鮮鋭回復フィルタｆｂを作成する（Ｓ１９０３）。
【００６９】
次に、ＣＰＵ１５１は、周波数特性計測チャートに対し、鮮鋭回復フィルタｆｂを用いた回復処理を、回復処理部１０７を制御して実行する。さらに、画像形成部１０８を制御し、回復処理後の画像を出力条件Ｏｉにより出力する（Ｓ１９０４）。
【００７０】
次に、ＣＰＵ１５１は、周波数ｕ１およびｕ２における回復処理後の明るさ（図１６、図１７のＬ１´，Ｌ２´）を取得する（Ｓ１９０５）。明るさの取得は、Ｓ１９０２と同様に行えばよい。
【００７１】
次に、ＣＰＵ１５１は、Ｌ１´，Ｌ２，Ｌ２´を用い、式（５）から、Ｌ２´´≒Ｌ１´となる振幅ゲインχを算出する（Ｓ１９０６）。
【００７２】
次に、ＣＰＵ１５１は、Ｓ１９０６で算出したχを用いて、式（６）に従い、鮮鋭回復フィルタの周波数特性Ｆｃ（ｕ）を算出する。そして、Ｆｃ（ｕ）に逆フーリエ変換を行うことで、明るさ低下が検知されにくい鮮鋭回復フィルタｆｃを作成する（Ｓ１９０７）。
【００７３】
なお、実施例１と同様に、本実施例でも各出力条件Ｏｉに対してＳ１９０１〜Ｓ１９０７までの処理を行い、作成したフィルタは対応ＬＵＴに対応付けてパラメータ保持部１０５に格納する。
（変形例）
上記実施例１では振幅ゲインに着目し、所定周波数の振幅ゲインの値を出力解像度に対応する周波数まで保持することにより、鮮鋭回復処理後の明るさの低下を抑制する方法を示した。また、上記実施例２では、各周波数における振幅ゲインと鮮鋭回復処理後の明るさの関係に着目し、所定周波数における明るさから出力解像度に対応する周波数における明るさを決定し、明るさが低下しないように所定周波数以上の領域の振幅ゲインを制御することで明るさ変動が知覚されにくくする方法を示した。本実施例では振幅ゲインと明るさの両方を制御することで、鮮鋭性と明るさ変動のバランスをとる方法を示す。
【００７４】
鮮鋭回復処理後の明るさの周波数特性Ｌ（ｕ）が、逆フィルタよりも明るさの低下が少なく、かつ、明るさ変動が知覚されにくい特性となっていれば、所定周波数以上の領域でフィルタの特性を必ずしも略線形とする必要はない。ただし、回復処理に用いるフィルタの特性に変曲点（２次微分の符号が変化する点）、不連続点（微分不可能な点）が存在すると、それらの点の近傍で不自然な明るさの変動が発生する。そのため、変曲点や不連続点を含まない周波数特性を用いて回復処理を行うことが望ましい。
【００７５】
図２０を用いて、回復処理に適用可能な周波数特性について、より詳細に説明する。
【００７６】
図２０（Ａ）、（Ｂ）、（Ｃ）中の曲線（Ａ−１）〜（Ａ−４）、（Ｂ−１）〜（Ｂ−３）、（Ｃ−１）〜（Ｃ−２）は、ｕ１＜ｕ＜ｕ２におけるフィルタの周波数特性を示したものである。それぞれの詳細については後述する。一方、図中の曲線（ａ−１）〜（ａ−４）、（ｂ−１）〜（ｂ−３）、（ｃ−１）〜（ｃ−２）は、それぞれ前記フィルタ特性（Ａ−１）〜（Ｃ−２）を用いて鮮鋭回復処理した後の明るさの特性Ｌ（ｕ）を示している。
【００７７】
フィルタの周波数特性（Ａ−１）は、ｕ１＜ｕ＜ｕ２における周波数特性Ｆ（ｕ）が１．０と一定の場合である。この場合、鮮鋭回復処理による明るさの低下は生じないが、鮮鋭回復処理による効果も得られない。
【００７８】
フィルタの周波数特性（Ａ−２）は、逆特性Ｒ（ｕ）を示している。明るさの特性（ａ−２）が示すように、逆特性Ｒ（ｕ）を用いて鮮鋭回復処理を行うと、高周波領域において、明るさの大きな低下が発生する。
【００７９】
フィルタの周波数特性（Ａ−３）は、実施例１における鮮鋭回復フィルタｆｂの周波数特性Ｆｂ（ｕ）を示している。周波数特性（Ａ−３）は、ｕ１＜ｕ＜ｕ２において、周波数特性がＦ（ｕ１）に保持されている。明るさの特性（ａ−３）が示すように、保持された周波数領域では、明るさは略線形に緩やかに変化する。そのため、周波数特性（Ａ−３）を用いることで、明るさ変動が知覚されにくくなる。
【００８０】
フィルタの周波数特性（Ａ−４）は、実施例２における回復フィルタの周波数特性Ｆｃ（ｕ）を示している。周波数特性（Ａ−４）の周波数ｕ２における振幅ゲインは、式（６）により、回復処理後の周波数ｕ２における明るさと、回復処理後の周波数ｕ１における明るさが一致するように決定されている。この時、周波数特性（Ａ−４）は、ｕ１＜ｕ＜ｕ２において、略線形に変化する。周波数特性を略線形に変化させると、回復処理後の明るさも、略線形に変化することが本願発明者の実験によって明らかとなっている。回復処理後の周波数ｕ２における明るさと周波数ｕ１における明るさが一致しすることから、周波数特性（Ａ−４）を用いると、（ａ−４）に示すように、ｕ１＜ｕ＜ｕ２において、明るさが略一致する。そのため、周波数特性（Ａ−４）を用いることで、明るさの変動が極めて知覚されにくい処理が実現できる。
【００８１】
上述の通り、明るさ変動が知覚されにくい特性となっていれば、周波数ｕ１とｕ２における回復処理後の明るさを、必ずしも略一致させる必要はない。加えて、必ずしもフィルタの特性を略線形に変化させる必要もない。たとえば、周波数特性（Ｂ−１）〜（Ｂ−３）を回復フィルタの周波数特性として用いてもよい。
【００８２】
周波数特性（Ｂ−１）は、ｕ１＜ｕ＜ｕ２において、略線形に変化するが、周波数ｕ２における周波数特性（Ｂ−１）の振幅ゲインは、周波数特性（Ａ−４）の周波数ｕ２における振幅ゲインより大きい。したがって、回復処理後の明るさ（ｂ−１）は、周波数に対して略線形に低下するが、周波数特性（Ａ−４）を用いた回復処理よりも、回復効果は大きい。そのため、明るさの低下の抑制より、処理の効果を優先する場合には、周波数特性（Ｂ−１）が好ましい。なお、不図示であるが、周波数ｕ２における振幅ゲインを、周波数特性（Ａ−４）の周波数ｕ２における振幅ゲインより小さくした周波数特性（Ｂ−１）´を回復処理に用いた場合、明るさの低下をより抑制できるが、周波数ｕ１付近で明るさの特性が下に凸となり不自然である。よって、周波数特性（Ｂ−１）´より、周波数特性（Ａ−４）を用いて回復処理したほうが、ｕ１における明るさ変動が知覚されにくく、より好ましい。
【００８３】
周波数特性（Ｂ−２）、（Ｂ−３）は、ｕ１＜ｕ＜ｕ２において、非線形に変化する。このとき、明るさの特性は、（ｂ−２）、（ｂ−３）に示す通り、周波数に対してなだらかに低下する。そのため、明るさの変動として知覚されにくい。
【００８４】
一方、周波数特性（Ｃ−１）に示すように、回復処理に用いるフィルタの特性に変曲点（２次微分の符号が変化する点）があると、明るさの特性（ｃ−１）にも変曲点が生まれる。明るさの特性（ｃ−１）の変曲点近傍では、マッハバンドとして知られる明るさの変動として知覚されやすくなる。そのため、回復処理に用いるフィルタの特性は変曲点がないことが望ましい。加えて、マッハバンドは、明るさの不連続点でも知覚されやすい。そのため、たとえば、周波数特性（Ｃ−２）に示すような不連続な周波数特性で回復処理を行うと、明るさの特性（ｃ−２）にも不連続な点が生じ、明るさの変動が知覚されやすくなる。したがって、回復処理に用いるフィルタの特性には、所定周波数以上の領域に不連続点がないことが望ましい。
【００８５】
上記（Ａ）、（Ｂ）に示すような、変曲点や不連続点を含まない周波数特性を用いて回復処理を行うことで、明るさの変動が知覚されにくい鮮鋭回復処理が実現できる。
【００８６】
なお、上記実施例の説明では、簡潔な説明のため、出力画像の周波数特性を１次元の関数Ｐ（ｕ）として扱っているが、プリンタの周波数特性は異方性を持つことが知られており、実際には、周波数ｕと、ｕに直交する方向の周波数ｖの２次元の関数Ｐ（ｕ，ｖ）のとして、回復フィルタの周波数特性を設計することが望ましい。
【００８７】
また、上記実施例においては、所定周波数ｕ１として視覚限界周波数ｕ＿ｖｉｓを用いたが、ｕ１はｕ＿ｖｉｓに限定されるものではない。たとえば、視覚感度ではなく、処理前後の明度の低下明るさの変化量からｕ１を決定してもよい。すなわち、処理前後の明るさ変化量が、離間比較で差の弁別が困難とされるＬ＊＝２となる周波数をｕ１として用いてもよい。
【００８８】
さらに、上記実施例においては、周波数ｕ２として画像形成装置の出力解像度に対応する周波数ｕ＿ｍａｘを用いた。しかし、画像形成装置と出力条件Ｏの組み合わせによっては、画像形成装置が出力解像度に相当する画像を形成できない場合がある。たとえば、画像形成装置としてインクジェットプリンタを用い、記録媒体として専用用紙ではなく普通紙を用いるような場合、インクのにじみが大きく、高周波な入力信号は画像として形成されない。なお、周波数ｕ＿ｍａｘ´までの信号しか画像形成装置が形成できない場合、周波数ｕ＿ｍａｘ´より高い周波数では、信号の平均明度に基づく略均一な画像として形成される。したがって周波数ｕ＿ｍａｘ´より高周波な領域では、明るさ低下の度合いが周波数ｕ＿ｍａｘ´と略一致する。従って、周波数ｕ＿ｍａｘ´より高周波な領域では図８（ｂ）に示すような線形性が成り立たない。よって、このような場合には、画像形成装置の出力解像度に対応する周波数ｕ＿ｍａｘではなく、前記周波数ｕ＿ｍａｘ´、もしくは周波数ｕ＿ｍａｘ´より低い適切な周波数をｕ２として用い、緩やかに変化するようにするとよい。
【００８９】
また、本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【図１】
0006632222.tif 000002
【図２】
0006632222.tif 000003
【図３】
0006632222.tif 000004
【図４】
0006632222.tif 000005
【図５】
0006632222.tif 000006
【図６】
0006632222.tif 000007
【図７】
0006632222.tif 000008
【図８】
0006632222.tif 000009
【図９】
0006632222.tif 000010
【図１０】
0006632222.tif 000011
【図１１】
0006632222.tif 000012
【図１２】
0006632222.tif 000013
【図１３】
0006632222.tif 000014
【図１４】
0006632222.tif 000015
【図１５】
0006632222.tif 000016
【図１６】
0006632222.tif 000017
【図１７】
0006632222.tif 000018
【図１８】
0006632222.tif 000019
【図１９】
0006632222.tif 000020
【図２０】
0006632222.tif 000021