JP2004064690A

JP2004064690A - 画像評価装置及び画像評価方法

Info

Publication number: JP2004064690A
Application number: JP2002223731A
Authority: JP
Inventors: Takeki Kutsuma; 沓間　丈輝
Original assignee: Canon Inc
Current assignee: Canon Inc
Priority date: 2002-07-31
Filing date: 2002-07-31
Publication date: 2004-02-26

Abstract

【課題】プリンタで媒体上に印刷されたグラデーション画像の品質の評価を心理評価と対応した物理評価によって好適に行うことができる画像評価装置及び画像評価方法を提供する。
【解決手段】画像入力部１０から所定のグラデーション画像が入力され、色空間変換部１１でグラデーション画像の色空間を変換して第１の明度画像が取得される。次いで、フーリエ変換処理部１２で第１の明度画像がフーリエ変換され第１の空間周波数情報を取得し、視覚特性乗算部１３で第１の空間周波数情報に人間の視覚特性に関する情報を乗算する。さらに、逆フーリエ変換処理部１４で逆フーリエ変換された第２の明度画像に基づいて、グラデーション画像の階調性を評価するために階調性評価値出力部１７で評価値が算出される。
【選択図】　図１

Description

【０００１】
【発明の属する技術分野】
本発明は、プリンタによって媒体上に印刷されたグラデーション画像の品質を評価する画像評価装置及び画像評価方法に関する。
【０００２】
【従来の技術】
プリンタによって媒体上に印刷された画像の品質を評価する方法として、評価者が視覚に感じる程度を主観的に数量化する心理評価と、画像構造の性質を客観的に測定した量で評価する物理評価とが知られている。また、画像の品質に関する重要な要素の一つとして、例えば、グラデーション画像等の滑らかさのような「画像の階調性」がある。そして、画像の階調性を物理的に表現する尺度として、画素値が一定の大きさだけ離れた一定面積を有するパッチ画像の濃度域をＲＭＳ（Ｒｏｏｔ−Ｍｅａｎ−Ｓｑｕａｒｅ）粒状度の２倍で除算する「識別階調数」が従来から用いられている。
【０００３】
【発明が解決しようとする課題】
しかしながら、従来の技術によれば、プリンタによって媒体上に印刷される画像の識別階調数を計算することは可能であったが、シャドウ（暗）部からハイライト（明）部にかけて連続的に色や明度が変化するグラデーション画像の品質を客観的に評価することは困難であった。仮に、識別階調数の計算に使用されるパッチ画像の入力画素値の刻みを細かくした場合であっても、プリンタが印刷する際に印刷領域内にノイズが重畳されたり、プリンタのシェイディング特性から印刷位置によって色や明度の変動等が生じるために、従来の技術では、正確な階調性評価をすることができなかった。
【０００４】
さらに、従来のパッチ画像を用いる技術では、グラデーション画像のように色や明度が微小な領域で漸次的に変化する画像の識別階調数と、それらの画像に対する評価者の目視による識別階調数の印象とが異なるため、それぞれによる評価値の対応がよくないという問題があった。
【０００５】
本発明は、上記課題を解決するためになされたものであり、プリンタで媒体上に印刷されたグラデーション画像の品質の評価を心理評価と対応した物理評価によって好適に行うことができる画像評価装置及び画像評価方法を提供することを目的とする。
【０００６】
【課題を解決するための手段】
上記課題を解決するために、本発明に係る画像評価装置は、所定のグラデーション画像を入力する入力手段と、前記グラデーション画像の色空間を変換して第１の明度画像を取得する色空間変換手段と、前記第１の明度画像を周波数変換して第１の空間周波数情報を取得する周波数変換手段と、前記第１の空間周波数情報に所定の第２の空間周波数情報を乗算する周波数情報乗算手段と、前記第２の空間周波数情報が乗算された前記第１の空間周波数情報を第２の明度画像に逆周波数変換する逆周波数変換手段と、前記第２の明度画像に基づいて、前記グラデーション画像の階調性を評価する階調性評価手段とを備えることを特徴とする。
【０００７】
また、本発明に係る画像評価装置は、前記第２の空間周波数情報が、人間の視覚の明度に関する空間周波数特性に基づく空間周波数情報であることを特徴とする。
【０００８】
さらに、本発明に係る画像評価装置は、前記階調性評価手段が、前記第２の明度画像の所定領域ごとに平均明度情報を算出する平均明度算出手段と、隣接する前記所定領域間の明度差情報を算出する明度差算出手段と、前記所定領域ごとに明度偏差情報を算出する明度偏差算出手段と、前記明度偏差情報に基づいて前記グラデーション画像の階調性評価値を算出する階調性評価値算出手段とを備えることを特徴とする。
【０００９】
さらにまた、本発明に係る画像評価装置は、前記所定領域が、前記グラデーション画像のグラデーション方向について、該グラデーション画像を構成する全階調数で均等に分割した矩形領域であることを特徴とする。
【００１０】
さらにまた、本発明に係る画像評価装置は、前記入力手段が、所定の画像形成装置において所定の媒体上に形成されたグラデーション画像を読み取って入力することを特徴とする。
【００１１】
さらにまた、本発明に係る画像評価装置は、前記グラデーション画像が、カラーグラデーション画像であることを特徴とする。
【００１２】
さらにまた、本発明に係る画像評価装置は、前記カラーグラデーション画像が、Ｌ^＊ａ^＊ｂ^＊表色系又はＬ^＊ｕ^＊ｖ^＊表色系に基づいて表現された色空間に変換されることを特徴とする。
【００１３】
さらにまた、本発明に係る画像評価装置は、前記カラーグラデーション画像が、明度画像、色相画像及び彩度画像を含む色空間に変換されることを特徴とする。
【００１４】
【発明の実施の形態】
以下、図面を参照して、本発明の一実施形態に係る画像評価装置について説明する。
【００１５】
図１は、本発明の一実施形態に係る画像評価装置の構成を示すブロック図である。図１に示すように、本実施形態に係る画像評価装置は、プリンタ等で紙等の媒体上に印刷されたカラーグラデーション画像を読み取って入力する画像入力部１０、入力されたカラーグラデーション画像をＬ^＊ａ^＊ｂ^＊表色系に変換する色空間変換部１１、Ｌ^＊成分についてフーリエ変換を行ってスペクトル情報を算出するフーリエ変換処理部１２、人間の明度に対する視覚特性をスペクトル情報に乗算する視覚特性乗算部１３、視覚特性が乗算されたスペクトル情報を逆フーリエ変換する逆フーリエ変換処理部１４、逆フーリエ変換されたＬ^＊成分について所定ブロックごとの平均明度を算出する平均明度算出部１５、隣接するブロック間の平均明度差の標準偏差を算出する明度偏差算出部１６及び算出された各偏差に基づいて階調性評価値を算出する階調性評価値算出部１７で構成される。
【００１６】
すなわち、本発明に係る画像評価装置は、所定のグラデーション画像の色空間を変換して第１の明度画像を取得し、第１の明度画像を周波数変換して第１の空間周波数情報を取得し、第１の空間周波数情報に所定の第２の空間周波数情報を乗算し、第２の空間周波数情報が乗算された第１の空間周波数情報を第２の明度画像に逆周波数変換し、第２の明度画像に基づいて、グラデーション画像の階調性を評価することを特徴とする。
【００１７】
また、本発明に係る画像評価装置は、第２の空間周波数情報が、人間の視覚の明度に関する空間周波数特性に基づく空間周波数情報であることを特徴とする。
【００１８】
さらに、本発明に係る画像評価装置は、階調性評価が、第２の明度画像の所定領域ごとに平均明度情報を算出し、隣接する所定領域間の明度差情報を算出し、所定領域ごとに明度偏差情報を算出し、明度偏差情報に基づいてグラデーション画像の階調性評価値を算出することを特徴とする。
【００１９】
次に、上述した構成の画像評価装置を用いてカラーグラデーション画像の階調性評価値を求める動作について説明する。図２は、本実施形態の画像評価装置の動作手順を説明するためのフローチャートである。
【００２０】
まず、画像入力部１０によって階調性評価用のカラーグラデーション画像が解像度２０００ｄｐｉで読み込まれ、ＲＧＢ画像データとして不図示の記憶部等へ保存される（ステップＳ２００）。図３は、本実施形態において画像入力部１０から入力される媒体上に印刷された階調性評価用のカラーグラデーション画像のモデルを示す図である。尚、図３に示されるカラーグラデーション画像は、解像度１０００ｄｐｉのプリンタで出力された場合のモデルを示す図である。図３に示されるモデルは、横方向の画素値が９４７２画素（２４０．５８８８ｍｍ）、縦方向の画素値が１０２４画素（２１７ｍｍ）である。また、図４は、図３に示されるモデルを画像入力部１０によって解像度２０００ｄｐｉで読み込み、予め決められた大きさで切り出されたカラーグラデーション画像を示す図である。
【００２１】
すなわち、本発明は、所定領域が、グラデーション画像のグラデーション方向について、グラデーション画像を構成する全階調数で均等に分割した矩形領域であることを特徴とする。また、本発明に係る画像評価装置は、所定の画像形成装置において所定の媒体上に形成されたグラデーション画像を読み取ることを特徴とする。さらに、本発明は、グラデーション画像が、カラーグラデーション画像であることを特徴とする。
【００２２】
次に、色空間変換部１１では、式（１）から式（５）を用いて、ＲＧＢデータをＣＩＥ　ＬＡＢに変換する（ステップ２０１）。
【００２３】
【数１】

【００２４】
ここで、Ｒ’，Ｇ’，Ｂ’≦０．０４０４５の場合、式（２）を用いる。
【００２５】
【数２】

【００２６】
一方、Ｒ’，Ｇ’，Ｂ’＞０．０４０４５の場合、式（３）を用いる。
【００２７】
【数３】

【００２８】
次いで、式（４）を用いて、Ｘ、Ｙ、Ｚを計算する。
【００２９】
【数４】

【００３０】
そして、式（５）を用いて、Ｌ^＊、ａ^＊、ｂ^＊が算出される。
【００３１】
【数５】

【００３２】
尚、Ｘ、Ｙ、Ｚは、試料のＸＹＺ表色系における３刺激値である。また、Ｘｎ、Ｙｎ、Ｚｎは、完全拡散反射面の３刺激値である。また、式（５）において、（Ｘ／Ｘｎ），（Ｙ／Ｙｎ），（Ｚ／Ｚｎ）＞０．００８８５６とする。
【００３３】
尚、（Ｘ／Ｘｎ），（Ｙ／Ｙｎ），（Ｚ／Ｚｎ）＜０．００８８５６の場合は、式（６）に示すようにそれぞれを置き換えて計算する。
【００３４】
【数６】

【００３５】
上述した手順で変換されたＬ^＊ａ^＊ｂ^＊画像は不図示の記録部に保存される。また、図４に示すように、同一画素列にある１０２４画素の平均明度Ｌ^＊が求められる（ステップＳ２０２）。
【００３６】
また、図５は、Ｌ^＊成分の画像データについて、グラデーションが変化する方向（横方向）の画素値である１８９４４画素を８ビットの全階調数である２５６で除算して、１階調ごとに７４画素幅毎の階調レベルのブロック（以下、「階調ブロック」と称す。）に分割するブロック分割の概要を示す図である。また、図６は、各階調レベルと階調ブロック内の平均明度情報との関係を示す図である。次に、フーリエ変換処理部１２では、ステップＳ２０２で求められた平均明度Ｌ^＊について式（７）に従ってフーリエ変換が行われ、スペクトルＳ（ｆｒ）が求められる（ステップＳ２０３）。尚、本実施形態では周波数変換の一例としてフーリエ変換を用いているが、その他の周波数変換であっても同様である。
【００３７】
【数７】

【００３８】
但し、式（７）において、Ｎは、図５に示される図における横方向の全画素値（１８９４４画素）であり、ｊは虚数単位である。
【００３９】
フーリエ変換が行われた後、周波数領域のカウンタｆｒが０に初期化される（ステップＳ２０４）。さらに、フーリエ変換処理部１２では、式（８）を用いて空間周波数ｆ’（ｃｙｃｌｅｓ／ｄｅｇｒｅｅ）が計算される（ステップＳ２０５）。
【００４０】
【数８】

【００４１】
但し、ｄｐｉは入力解像度であり、本実施形態では２０００ｄｐｉである。また、Ｒは明視距離であり、本実施形態ではＲ＝３００ｍｍとする。
【００４２】
次に、視覚特性乗算部１３では、ステップＳ２０３で求められた空間周波数ｆ’のスペクトル情報Ｓ（ｆ’）に対して、式（９）で定義される人間の視覚の明度に対する空間周波数特性ＶＴＦを乗算したスペクトルＳ’（ｆ）が計算される（ステップＳ２０６）。ここで、図７は、空間周波数と人間の明度に対する視覚特性である相対視覚感度との関係を示す図である。
【００４３】
【数９】

【００４４】
尚、本実施形態では、ＶＴＦの形状として式（９）で示されるものを用いているが、その他のＶＴＦの視覚特性のものを用いても同様の手順で処理が可能である。
【００４５】
次に、カウンタｆｒがＮ未満であるか否かが判定される（ステップＳ２０７）。その結果、カウンタｆｒがＮ未満の場合（Ｙｅｓ）、カウンタｆｒをインクリメントしてステップＳ２０６に戻る（ステップＳ２０８）。一方、カウンタｆｒがＮである場合（Ｎｏ）、ステップＳ２０９に進む。そして、逆フーリエ変換処理部１４においてスペクトルＳ’（ｆ）が式（１０）に従って逆フーリエ変換され、実空間の平均明度情報Ｌ’^＊が求められる（ステップＳ２０９）。
【００４６】
【数１０】

【００４７】
次に、平均明度算出部１５において、画像入力部１０により解像度２０００ｄｐｉで読み込まれた図４に示すようなカラーグラデーション画像から、グラデーションが変化する方向（横方向）について、一階調レベルごとの階調ブロックの階調幅（７４画素）が算出される（ステップ２１０）。
【００４８】
次に、階調番号を表すカウンタＩを０に初期化する（ステップＳ２１１）。そして、階調レベルＩの階調ブロック内の全画素数７４画素の平均明度Ｌ^＊（Ｉ）が求められる（ステップＳ２１２）。その後、カウンタＩが全階調数である２５６未満であるか否かが判定される（ステップＳ２１３）。その結果、カウンタＩが２５６未満の場合（Ｙｅｓ）、カウンタＩをインクリメントしてステップＳ２１２に戻る（ステップＳ２１４）。一方、カウンタＩが２５６である場合（Ｎｏ）、ステップＳ２１５に進む。
【００４９】
ステップＳ２１５では、カウンタＩが０に初期化される。そして、ＶＴＦが乗算された平均明度情報Ｌ^＊から、式（１１）によって図８のモデル図に示すような隣接する階調間の明度差が求められる（ステップＳ２１６）。図８は、本実施形態におけるグラデーション画像の隣接階調レベル間の明度差のモデルを示す図である。
【００５０】
【数１１】

【００５１】
次に、カウンタＩが２５５未満であるか否かが判定される（ステップＳ２１７）。その結果、カウンタＩが２５５未満の場合（Ｙｅｓ）、カウンタＩをインクリメントしてステップＳ２１６に戻る（ステップＳ２１８）。一方、カウンタＩが２５５の場合（Ｎｏ）、明度偏差算出部１６では全階調の明度差の標準偏差が求められ（ステップＳ２１９）、式（１２）に示すように、標準偏差に対数をとり、符号反転を行って階調性評価値を階調性評価値出力部１７で出力する（ステップＳ２２０）。このように、本実施形態では、主観評価との対応をよくするために、標準偏差に対数をとり、符号反転を行った。
【００５２】
【数１２】

【００５３】
上記手順によって、本実施形態に係る画像評価装置を用いて、入力されたカラーグラデーション画像の階調性評価値を算出することができ、その結果、プリンタ印刷画像の客観的な階調性評価を行うことが可能となった。
【００５４】
また、本実施形態ではカラーグラデーションの入力画素値Ｉを２５６階調としたが、本発明は２５６階調のグラデーション画像以外であっても同様に適用することが可能である。
【００５５】
また、本実施形態では、平均明度算出部１５により、平均明度を算出した後に１次元のＶＴＦの乗算を行ったが、平均明度の算出前にグラデーション画像に対して２次元のフーリエ変換を行い、２次元の空間周波数情報に対して２次元のＶＴＦを乗算して平均明度を算出し、その後明度偏差を求めるような方法を用いることも可能である。
【００５６】
また、上述した実施形態では、カラー画像として、Ｌ^＊ａ^＊ｂ^＊表色系について説明したが、Ｌ^＊ｕ^＊ｖ^＊表色系等を用いても同様に本発明を適用することが可能である。また、本発明においては、カラー画像が、明度画像、色相画像及び彩度画像から構成されているものであっても同様に適用することが可能である。
【００５７】
すなわち、本発明は、カラーグラデーション画像が、Ｌ^＊ａ^＊ｂ^＊表色系又はＬ^＊ｕ^＊ｖ^＊表色系に基づいて表現された色空間に変換されることを特徴とする。また、本発明は、カラーグラデーション画像が、明度画像、色相画像及び彩度画像を含む色空間に変換されることを特徴とする。
【００５８】
尚、本発明は、複数の機器（例えば、ホストコンピュータ、インタフェース機器、リーダ、プリンタ等）から構成されるシステムに適用しても、一つの機器からなる装置（例えば、複写機、ファクシミリ装置等）に適用してもよい。
【００５９】
また、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記録媒体（または記憶媒体）を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。また、コンピュータが読み出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム（ＯＳ）などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００６０】
さらに、記録媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００６１】
本発明を上記記録媒体に適用する場合、その記録媒体には、先に説明したフローチャートに対応するプログラムコードが格納されることになる。
【００６２】
【発明の効果】
以上説明したように、本発明によれば、プリンタで媒体上に印刷されたグラデーション画像の品質の評価を心理評価と対応した物理評価によって好適に行うことができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る画像評価装置の構成を示すブロック図である。
【図２】本実施形態の画像評価装置の動作手順を説明するためのフローチャートである。
【図３】本実施形態において画像入力部１０から入力される媒体上に印刷された階調性評価用のカラーグラデーション画像のモデルを示す図である。
【図４】図３に示されるモデルを画像入力部１０によって解像度２０００ｄｐｉで読み込み、予め決められた大きさで切り出されたカラーグラデーション画像を示す図である。
【図５】Ｌ^＊成分の画像データについて、グラデーションが変化する方向の画素値である１８９４４画素を８ビットの全階調数である２５６で除算して、１階調ごとに７４画素幅毎の階調レベルのブロックに分割するブロック分割の概要を示す図である。
【図６】各階調レベルと階調ブロック内の平均明度情報との関係を示す図である。
【図７】空間周波数と人間の明度に対する視覚特性である相対視覚感度との関係を示す図である。
【図８】本実施形態におけるグラデーション画像の隣接階調レベル間の明度差のモデルを示す図である。
【符号の説明】
１０　画像入力部
１１　色空間変換部
１２　フーリエ変換処理部
１３　視覚特性乗算部
１４　逆フーリエ変換処理部
１５　平均明度算出部
１６　明度偏差算出部
１７　階調性評価値出力部

Claims

所定のグラデーション画像を入力する入力手段と、
前記グラデーション画像の色空間を変換して第１の明度画像を取得する色空間変換手段と、
前記第１の明度画像を周波数変換して第１の空間周波数情報を取得する周波数変換手段と、
前記第１の空間周波数情報に所定の第２の空間周波数情報を乗算する周波数情報乗算手段と、
前記第２の空間周波数情報が乗算された前記第１の空間周波数情報を第２の明度画像に逆周波数変換する逆周波数変換手段と、
前記第２の明度画像に基づいて、前記グラデーション画像の階調性を評価する階調性評価手段と
を備えることを特徴とする画像評価装置。
前記第２の空間周波数情報が、人間の視覚の明度に関する空間周波数特性に基づく空間周波数情報であることを特徴とする請求項１記載の画像評価装置。
前記階調性評価手段が、
前記第２の明度画像の所定領域ごとに平均明度情報を算出する平均明度算出手段と、
隣接する前記所定領域間の明度差情報を算出する明度差算出手段と、
前記所定領域ごとに明度偏差情報を算出する明度偏差算出手段と、
前記明度偏差情報に基づいて前記グラデーション画像の階調性評価値を算出する階調性評価値算出手段と
を備えることを特徴とする請求項１又は２に記載の画像評価装置。
前記所定領域が、前記グラデーション画像のグラデーション方向について、該グラデーション画像を構成する全階調数で均等に分割した矩形領域であることを特徴とする請求項３記載の画像評価装置。
前記入力手段が、所定の画像形成装置において所定の媒体上に形成されたグラデーション画像を読み取って入力することを特徴とする請求項１から４までのいずれか１項に記載の画像評価装置。
前記グラデーション画像が、カラーグラデーション画像であることを特徴とする請求項１から５までのいずれか１項に記載の画像評価装置。
前記カラーグラデーション画像が、Ｌ^＊ａ^＊ｂ^＊表色系又はＬ^＊ｕ^＊ｖ^＊表色系に基づいて表現された色空間に変換されることを特徴とする請求項６記載の画像評価装置。
前記カラーグラデーション画像が、明度画像、色相画像及び彩度画像を含む色空間に変換されることを特徴とする請求項６記載の画像評価装置。
所定のグラデーション画像の色空間を変換して第１の明度画像を取得する色空間変換工程と、
前記第１の明度画像を周波数変換して第１の空間周波数情報を取得する周波数変換工程と、
前記第１の空間周波数情報に所定の第２の空間周波数情報を乗算する周波数情報乗算工程と、
前記第２の空間周波数情報が乗算された前記第１の空間周波数情報を第２の明度画像に逆周波数変換する逆周波数変換工程と、
前記第２の明度画像に基づいて、前記グラデーション画像の階調性を評価する階調性評価工程と
を有することを特徴とする画像評価方法。
前記第２の空間周波数情報が、人間の視覚の明度に関する空間周波数特性に基づく空間周波数情報であることを特徴とする請求項９記載の画像評価方法。
前記階調性評価工程が、
前記第２の明度画像の所定領域ごとに平均明度情報を算出する平均明度算出工程と、
隣接する前記所定領域間の明度差情報を算出する明度差算出工程と、
前記所定領域ごとに明度偏差情報を算出する明度偏差算出工程と、
前記明度偏差情報に基づいて前記グラデーション画像の階調性評価値を算出する階調性評価値算出工程と
を有することを特徴とする請求項９又は１０に記載の画像評価方法。
前記所定領域が、前記グラデーション画像のグラデーション方向について、該グラデーション画像を構成する全階調数で均等に分割した矩形領域であることを特徴とする請求項１１記載の画像評価方法。
所定の画像形成装置において所定の媒体上に形成されたグラデーション画像を読み取る読み取り工程をさらに有することを特徴とする請求項９から１２までのいずれか１項に記載の画像評価方法。
前記グラデーション画像が、カラーグラデーション画像であることを特徴とする請求項９から１３までのいずれか１項に記載の画像評価方法。
前記カラーグラデーション画像が、Ｌ^＊ａ^＊ｂ^＊表色系又はＬ^＊ｕ^＊ｖ^＊表色系に基づいて表現された色空間に変換されることを特徴とする請求項１４記載の画像評価方法。
前記カラーグラデーション画像が、明度画像、色相画像及び彩度画像を含む色空間に変換されることを特徴とする請求項１４記載の画像評価方法。
コンピュータに、
所定のグラデーション画像の色空間を変換して第１の明度画像を取得する色空間変換手順と、
前記第１の明度画像を周波数変換して第１の空間周波数情報を取得する周波数変換手順と、
前記第１の空間周波数情報に所定の第２の空間周波数情報を乗算する周波数情報乗算手順と、
前記第２の空間周波数情報が乗算された前記第１の空間周波数情報を第２の明度画像に逆周波数変換する逆周波数変換手順と、
前記第２の明度画像に基づいて、前記グラデーション画像の階調性を評価する階調性評価手順と
を実行させるためのプログラム。
請求項１７記載のプログラムを格納することを特徴とするコンピュータ読み取り可能な記録媒体。