(57)【要約】【課題】メタリックインクを用いた印刷結果の画質を、総合的かつ観察者の観察条件や主観をより適切に反映した態様で評価する。【解決手段】印刷媒体にメタリックインクを付着させた印刷結果に基づく画質評価方法であって、少なくとも、上記印刷結果の読み取りにより得られた明度に基づいて算出される空間的ノイズと、上記印刷結果における特定角度の反射光の明度に応じた第１補正項と、に基づいて定義される第１評価指標と、少なくとも、上記印刷結果における正反射光の明度と、上記特定角度の反射光の明度に応じた第２補正項と、に基づいて定義される第２評価指標と、を算出する。【選択図】図２
2013213712.tif 000002
【特許請求の範囲】
【請求項１】
印刷媒体にメタリックインクを付着させた印刷結果に基づく画質評価方法であって、
少なくとも、上記印刷結果の読み取りにより得られた明度に基づいて算出される空間的ノイズと、上記印刷結果における特定角度の反射光の明度に応じた第１補正項と、に基づいて定義される第１評価指標と、少なくとも、上記印刷結果における正反射光の明度と、上記特定角度の反射光の明度に応じた第２補正項と、に基づいて定義される第２評価指標と、を算出することを特徴とする画質評価方法。
【請求項２】
上記特定角度は、上記印刷結果を複数の観察角度で測定することにより得られた明度に対する評価のうち、上記印刷結果の観察者による明度の評価と相関が相対的に高い評価に対応する角度であることを特徴とする請求項１に記載の画質評価方法。
【請求項３】
上記正反射光の明度と、上記第２補正項と、上記印刷結果からの反射輝度を複数の角度で測定することにより得られた反射輝度の分布をロジスティック分布を含む関数によりシミュレートした結果におけるピーク値に応じた第３補正項と、に基づいて上記第２評価指標を算出することを特徴とする請求項１または請求項２に記載の画質評価方法。
【請求項４】
上記空間的ノイズと、上記第１補正項と、上記印刷結果からの反射輝度を複数の角度で測定することにより得られた反射輝度の分布をロジスティック分布を含む関数によりシミュレートした結果から得られる複数のパラメーターに基づく第４補正項と、に基づいて上記第１評価指標を算出することを特徴とする請求項１〜請求項３のいずれかに記載の画質評価方法。
【請求項５】
上記複数のパラメーターは、上記関数の輝度方向の倍率および、上記ロジスティック分布の標準偏差に比例するパラメーターであることを特徴とする請求項４に記載の画質評価方法。
【請求項６】
上記第１評価指標と上記第２評価指標とを各軸とする２次元平面に上記印刷結果に対する画質評価の結果を表示させた画像を所定の画面に出力させることを特徴とする請求項１〜請求項５のいずれかに記載の画質評価方法。
【請求項７】
印刷媒体にメタリックインクを付着させた印刷結果に基づく画質評価を実行する画質評価装置であって、
少なくとも、上記印刷結果の読み取りにより得られた明度に基づいて算出される空間的ノイズと、上記印刷結果における特定角度の反射光の明度に応じた第１補正項と、に基づいて定義される第１評価指標と、少なくとも、上記印刷結果における正反射光の明度と、上記特定角度の反射光の明度に応じた第２補正項と、に基づいて定義される第２評価指標と、を算出することを特徴とする画質評価装置。
【請求項８】
印刷媒体にメタリックインクを付着させた印刷結果に基づく画質評価をコンピューターに実行させる画質評価プログラムであって、
少なくとも、上記印刷結果の読み取りにより得られた明度に基づいて算出される空間的ノイズと、上記印刷結果における特定角度の反射光の明度に応じた第１補正項と、に基づいて定義される第１評価指標と、少なくとも、上記印刷結果における正反射光の明度と、上記特定角度の反射光の明度に応じた第２補正項と、に基づいて定義される第２評価指標と、を算出する機能を実行させることを特徴とする画質評価プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画質評価方法、画質評価装置および画質評価プログラムに関する。
【背景技術】
【０００２】
金属物質が放つ特殊な光沢感や質感を印刷結果において表現するためのメタリックインクを用いた印刷が行われている。また、メタリックインクを用いた印刷結果の画質を評価する手法がいくつか提案されている。
関連技術として、観察時の測定エリアに対する光源の照射角度を異ならせることによりそれぞれ得られた第１画像および第２画像に基づいて算出した、粒子感に関する定量化された第１の評価指数および第２の評価指数から、粒子感の角度依存性に関する角度特性パラメーターを求めるメタリック塗装面評価方法が知られている（特許文献１参照）。
また、複数角度の反射明度を用いて、主成分分析により得た第一主成分を輝度感、第二主成分を白味感として扱い、第一、第二主成分平面へプロットすることでメタリック塗装のメタリック感を評価する方法が知られている（特許文献２参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００６‐２０８３２７号公報
【特許文献２】特開２００４‐２８６６７２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
メタリックインクを用いた印刷結果の画質をより高精度に評価するには、文献１のような粒子感といった１つの指標だけでなく、画質を評価するための複数の指標により総合的に評価する必要がある。また、複数の指標を用いて評価する場合であっても、それぞれの指標が、観察者の観察条件や主観をより適切に反映したものでなければ、画質評価の精度を高めることは難しい。
【０００５】
本発明は少なくとも上記課題を解決するためになされたものであり、従来に増して、メタリックインクを用いた印刷結果の画質を、総合的かつ観察者の観察条件や主観をより適切に反映した態様で評価することが可能な画質評価方法、画質評価装置および画質評価プログラムを提供する。
【課題を解決するための手段】
【０００６】
本発明の態様の一つは、印刷媒体にメタリックインクを付着させた印刷結果に基づく画質評価方法であって、少なくとも、上記印刷結果の読み取りにより得られた明度に基づいて算出される空間的ノイズと、上記印刷結果における特定角度の反射光の明度に応じた第１補正項と、に基づいて定義される第１評価指標と、少なくとも、上記印刷結果における正反射光の明度と、上記特定角度の反射光の明度に応じた第２補正項と、に基づいて定義される第２評価指標と、を算出する。
当該構成によれば、空間的ノイズを主要因とした第１評価指標および正反射光の明度を主要因とした第２評価指標によって、メタリックインクの印刷結果を総合的に評価することができる。また、第１評価指標および第２評価指標はいずれも、上記特定角度の反射光の明度に応じた補正項によって補正されているため、観察者の観察条件や主観をより反映した値となり、評価の精度が向上する。
【０００７】
本発明の態様の一つは、上記特定角度は、上記印刷結果を複数の観察角度で測定することにより得られた明度に対する評価のうち、上記印刷結果の観察者による明度の評価と相関が相対的に高い評価に対応する角度であるとしてもよい。当該構成によれば、上記特定角度は、観察者が印刷結果を観察したときに最も明度を評価し易い角度であるため、上記特定角度の反射光の明度に応じた補正項を採用することで、観察者の観察条件や主観をより反映した評価結果を得ることができる。
【０００８】
本発明の態様の一つは、上記正反射光の明度と、上記第２補正項と、上記印刷結果からの反射輝度を複数の角度で測定することにより得られた反射輝度の分布をロジスティック分布を含む関数によりシミュレートした結果におけるピーク値に応じた第３補正項と、に基づいて上記第２評価指標を算出する。当該構成によれば、上記正反射光の明度を、上記シミュレート結果におけるピーク値に応じた第３補正項にて補正することで、観察者の観察条件や主観をより反映した評価結果を得ることができる。
【０００９】
本発明の態様の一つは、上記空間的ノイズと、上記第１補正項と、上記印刷結果からの反射輝度を複数の角度で測定することにより得られた反射輝度の分布をロジスティック分布を含む関数によりシミュレートした結果から得られる複数のパラメーターに基づく第４補正項と、に基づいて上記第１評価指標を算出する。当該構成によれば、上記シミュレート結果における複数のパラメーターを用いることで、“メタリック感の強さ”（以下、メタリック強度と呼ぶ。）といったものをある程度数値化することができ、このメタリック強度を示す第４補正項にて上記空間的ノイズを補正することで、観察者の主観や感じ方をより反映した評価結果を得ることができる。なお、上記複数のパラメーターは、例えば、上記関数の輝度方向の倍率および、上記ロジスティック分布の標準偏差に比例するパラメーターである。
【００１０】
本発明の態様の一つは、上記第１評価指標と上記第２評価指標とを各軸とする２次元平面に上記印刷結果に対する画質評価の結果を表示させた画像を所定の画面に出力させる。当該構成によれば、メタリックインクを用いたある印刷結果に対する画質評価であって第１評価指標および第２評価指標による評価を、当該画面を見るユーザーに対して判りやすく提示することができる。
【００１１】
本発明にかかる技術的思想は画質評価方法という態様のみで実現されるものではなく、例えば、当該方法を実行する装置（画質評価装置）によって実現されてもよいし、あるいは、当該方法をハードウェア（コンピューター）に実行させるプログラム（画質評価プログラム）によって実現されてもよい。
【図面の簡単な説明】
【００１２】
【図１】ハードウェア構成およびソフトウェア構成を示す図である。
【図２】画質評価処理を示すフローチャートである。
【図３】対象サンプルの測定結果から得られた反射輝度の分布を示す図である。
【図４】（２）式を用いたシミュレーションで得られた反射輝度の分布を示す図である。
【図５】（３）式を用いたシミュレーションで得られた反射輝度の分布を示す図である。
【図６】パラメーターａおよびｓがプロットされた２次元平面を示す図である。
【図７】画質評価確認画面の一例を示す図である。
【発明を実施するための形態】
【００１３】
以下、図面を参照しながら本発明の実施形態について説明する。
１．装置等の概略構成
図１は、本実施形態にかかる画質評価装置１等のハードウェア構成およびソフトウェア構成を示している。画質評価装置１は主としてコンピューター１０により構成される。画質評価装置１は、画質評価方法の実行主体である。また、画質評価装置１を複数装置からなるシステムとした場合、画質評価装置１には、ディスプレー３０等が含まれる。コンピューター１０およびディスプレー３０は、ラップトップ型ＰＣあるいはタブレット型コンピューター等のように、それらが一体的に構成されていてもよいし、互いに別体の装置であってもよい。またコンピューター１０には、プリンター２０を接続することが可能であり、コンピューター１０がプリンター２０を制御する印刷制御装置として機能してもよい。
【００１４】
コンピューター１０においては、ＣＰＵ１２が、ハードディスクドライブ（ＨＤＤ）１１等で実現される記憶部に記憶されたプログラム（画質評価プログラム１１ａ）を読み込み、当該プログラムをＲＡＭ１３等の内部メモリーに展開しながらプログラムに従った演算を実行することにより、輝度シミュレーション部１４、空間的ノイズ算出部１５、明度補正項算出部１６、メタリック強度補正項算出部１７、ピーク値補正項算出部１８、評価指標算出部１９等の各機能を実現する。これらの各機能については後述する。
【００１５】
また、コンピューター１０は、図示しないビデオインターフェイスを介してディスプレー３０と接続されており、図示しない入力インターフェイスを介してキーボードやマウスやタッチパッドやタッチパネル等の操作部４０と接続されている。さらにコンピューター１０は、後述するメタリックインクのサンプルを測定（測色）可能な測定装置５０やスキャナー６０と接続する。
【００１６】
２．画質評価処理
図２は、コンピューター１０において画質評価プログラム１１ａに従って実行される画質評価処理をフローチャートにより示している。画質評価処理では、評価対象としてユーザーによって任意に用意されたサンプル（対象サンプル）の画質評価を行なう。ここで言う「サンプル」とは、基本的にはメタリックインクが印刷された印刷媒体である。つまり、メタリックインクが付着した印刷媒体の表面の画質を評価する。当該評価においては、第１評価指標と第２評価指標の算出結果をもって一つの評価結果とする。本実施形態では、第１評価指標は、メタリックインクによる質感を評価するための一つの指標であって、空間的な微視的ノイズ（空間的ノイズ）を主要因とした指標である。第１評価指標は、ユーザーの感覚的には「キラキラ度合い」等と表現することができる。言い換えると、本実施形態では「キラキラ度合い」を評価するための指標として第１評価指標を算出する。「キラキラ度合い」は、サンプル表面の場所毎の光の濃淡や光の粒状感の密度等によって変化するものである。
【００１７】
また、本実施形態では、第２評価指標は、メタリックインクによる光沢感を評価するための一つの指標であって、正反射光の明度を主要因とした指標である。第２評価指標は、ユーザーの感覚的には「ピカピカ度合い」等と表現することができる。言い換えると、本実施形態では「ピカピカ度合い」を評価するための指標として第２評価指標を算出する。なお、本実施形態でいうメタリックインクとは、金属を含んで構成されたインクであり、例えば特開２０１０‐５２２２５号公報に開示されたものを用いることができる。
【００１８】
まず、ステップＳ１００では、輝度シミュレーション部１４が、測定装置５０を制御して測定装置５０に対象サンプルを測定させる。測定装置５０は、分光変角色差計（例えば、装置名：ＧＣ５０００）である。測定装置５０では、測定対象（サンプル）の測定面に対して鉛直方向を角度θ＝０としたとき、θ＝−４５°の角度から照射光を入射させ、測定面からの反射光を受光部で受光する。このとき、任意に設定した受光角度（反射角）の位置へ受光部を自動的に移動させることで、様々な反射角による分光反射率・分光透過率等を測定でき、測定結果から複数の角度における対象サンプルの各色彩値情報（ＣＩＥ‐ＸＹＺ表色系の値や、ＣＩＥ‐Ｌ*ａ*ｂ*表色系の値等）を算出することができる。従って、当該ステップＳ１００では、輝度シミュレーション部１４は、測定装置５０から出力されるこのような色彩値情報を入力する。
【００１９】
ステップＳ１１０では、輝度シミュレーション部１４は、測定装置５０より得られた対象サンプルの反射輝度Ｙの分布を、シミュレートする。
図３は、測定装置５０による対象サンプルの測定結果から得られた反射輝度Ｙの分布を示している。図３では、縦軸に反射輝度Ｙを取り、横軸に反射角θを取っている。ここで、当該分布を観察すると、何らかの確率分布モデルを適用することで、当該分布の特性をシミュレートできると考えられる。そこで本実施形態では、種々の考察の結果、下記の（１）式で表されるロジスティック分布を用いて当該分布の特性をシミュレートするとした。
【００２０】
【数１】
2013213712.tif 000003
ロジスティック分布は、期待値がｘp、標準偏差がπｓ／√３、となる確率密度関数である。上記（１）式においては、ｘは反射角θ、ｓは標準偏差に比例するパラメーターである。また、図３に示したような反射輝度Ｙの分布のシミュレートにあたっては、輝度方向の倍率を調整する係数ａ、輝度方向のオフセット量を調整する係数ｂ、角度依存の補正関数ｃｏｓθを用いて、正反射角（θ＝４５°）周辺を除いて反射輝度Ｙを下記の（２）式で表すことができると考えられる。
【数２】
2013213712.tif 000004
【００２１】
なお、補正関数ｃｏｓθは、反射面における反射強度の減衰を示す関数である。つまり、受光部で受光可能な反射強度は反射角θの変化に応じて変化すると考えられるため、ここでは反射輝度Ｙをｃｏｓθにより補正することで、予測される反射強度の減衰を加味している。上記（２）式における各パラメーター（ａ，ｂ，ｓ，ｘp）を様々に変化させるシミュレーションを行なうことで、図４に示すような結果が得られた。
【００２２】
図４では、上記（２）式における各パラメーター（ａ，ｂ，ｓ，ｘp）を変化させるシミュレーションの結果得られた反射輝度Ｙの分布を実線で示しており、図３に示した対象サンプルの実測結果から得られた反射輝度Ｙの分布を鎖線で併せて示している。図４から判るように、上記（２）式を採用することで、概ね良好なシミュレーションができる。なお、正反射角近傍（θ＝４０°，４５°，５０°）では、実測の輝度と図４に示したシミュレーション結果との差分を埋めるために、単純に輝度の差分量を加算する補正を行なう。
【００２３】
以上の検討を踏まえ、本実施形態では反射輝度Ｙを、下記の（３）式によってシミュレートするとした。
【数３】
2013213712.tif 000005
上記（３）式において、θ1は、４０°，４５°，５０°を表し、それ以外の反射角におけるｐ（θ1）は０であるとする。このような（３）式を用いるシミュレーションにおいては、いわゆる多変量非線形最適化を行なうことで各パラメーター（ａ，ｂ，ｓ，ｘp，ｐ（４０°），ｐ（４５°），ｐ（５０°））を決定する。その際、輝度シミュレーション部１４は、シミュレーションによる輝度を明度Ｌ*に変換し、測定装置５０による対象サンプルの実測明度Ｌ*との二乗誤差の積算が最小化するように各パラメーターを決定する。なお、光源方向近傍である反射角θ＝−４０°，−４５°，−５０°での明度Ｌ*は実測値の取得が困難であるため、これら反射角における上記二乗誤差の積算は除外する。
【００２４】
図５では、上記（３）式を用いるシミュレーションの結果得られた反射輝度Ｙの分布を実線で示しており、図３に示したサンプルの実測結果から得られた反射輝度Ｙの分布を鎖線で併せて示している。図５から判るように、上記（３）式を採用することで、正反射角近傍（θ＝４０°，４５°，５０°）も含めて良好なシミュレーションが行なわれた。
次に、ステップＳ１２０では、輝度シミュレーション部１４は、上記シミュレーションの結果得られた複数のパラメーターに基づいて、対象サンプルの“メタリック強度”を決定する。
【００２５】
図６は、上記複数のパラメーターのうち、ａおよびｓを２次元平面にプロットした様子を示している。図６では、異なる複数のサンプルの測定結果に基づいてそれぞれ行なわれた上述のシミュレーションにより得られたパラメーターａ，ｓで特定される点の分布を示している。ただし、図６において黒丸で示した点は、メタリックインクが付着された印刷媒体の測定結果に基づいて得られたパラメーターの位置を示しているが、白丸で示した点は、メタリックインクが付着されていない印刷媒体の測色結果に基づいて得られたパラメーターの位置を示している。このような黒丸および白丸の点を観察すると、黒丸の点については一定の特徴を有する分布を形成しており、白丸の点と大きく位置が離れていることから、黒丸の分布と白丸との間において、“メタリック感の強さ”が変化していると考察できる。そこで、図６ではこれら黒丸の点を近似した直線を実線にて示しており、当該直線に垂直であって右上を向く方向に略沿って、“メタリック感の強さ”すなわちメタリック強度が変化していると考察できる（図６における鎖線の矢印参照）。
【００２６】
上記近似直線の傾きは、計測の結果、１／−２６．８であるため、ベクトル（ａ，ｓ）＝（−２６．８，１）の方向が上記近似直線の方向となる。上記近似直線に垂直で右上を向いたベクトルは（１，２６．８）であるため、これとパラメーター（ａ，ｓ）との内積により、メタリック強度と想定される向きの程度を求めることができる。そこで本実施形態では、かかる内積の結果を、下記（４）式に示すようにメタリック強度ｘと定義する。
【数４】
2013213712.tif 000006
なお、（４）式では、上記内積の結果を定数２６．８で割った状態で示している。ただし、（４）式において定数ｕ＝−２６．８である。つまり、輝度シミュレーション部１４は、上記シミュレーションの結果得られた複数のパラメーターに基づいて、対象サンプルのメタリック強度ｘを算出する。
【００２７】
ステップＳ１３０では、空間的ノイズ算出部１５が、スキャナー６０を制御してスキャナー６０に対象サンプルの読み取りを実行させる。この結果、空間的ノイズ算出部１５は、対象サンプルの読取画像としてのＲＧＢデータをスキャナー６０から入力する。
ステップＳ１４０では、空間的ノイズ算出部１５は、スキャナー６０による読取画像に基づいて、対象パッチにおける空間的ノイズを算出する。この場合、空間的ノイズ算出部１５は、上記読取画像としてのＲＧＢデータを、ＩＣＣプロファイル等を用いてＬ*ａ*ｂ*表色系の画像データに変換する。そして、当該画像データのうちＬ*チャンネル（つまり明度画像）について、フーリエ変換等によって空間周波数へ変換し、ウィナースペクトルＷＳ（ｕ）を取得する。そして、ウィナースペクトルに、下記（５）式で表される視覚の空間周波数特性ＶＴＦ（ｕ）を適用する。
【数５】
2013213712.tif 000007
なお、ウィナースペクトルＷＳ（ｕ）および空間周波数特性ＶＴＦ（ｕ）におけるｕは周波数であり、また（５）式における観察距離ｌは、例えば３００ｍｍとする。
【００２８】
そして、空間的ノイズ算出部１５は、当該空間周波数特性によるフィルターを施したウィナースペクトルを積分することによって、下記（６）式に示すような空間的ノイズＬ*_Ｎが得られる。
【数６】
2013213712.tif 000008
【００２９】
次に、ステップＳ１５０では、メタリック強度補正項算出部１７が、上記メタリック強度ｘに基づく補正項ｆk（ｘ）を求める。補正項ｆk（ｘ）は、上記空間的ノイズＬ*_Ｎに対する補正項であり、第４補正項に該当する。上述したように、「キラキラ度合い」を評価するための第１評価指標は、空間的ノイズを主たる要素とするが、この空間的ノイズは評価対象における粒状性を示す基本的な情報であるため、非メタリックのサンプルにおいても空間的ノイズが大きい場合がある。一方、「キラキラ度合い」は、メタリック強度と想定される方向に依存して感じ易かったり感じ難かったりすると推測できる。そこで本実施形態では、上記メタリック強度ｘに基づく補正項ｆk（ｘ）を、空間的ノイズＬ*_Ｎに対する補正項として求めるようにしている。
【００３０】
メタリック強度ｘが高い場合に、当該補正項が大きくなる必要があるため、本実施形態では一例として、下記（７）式に示すような補正項ｆk（ｘ）を定義する。このような補正項ｆk（ｘ）は、メタリック強度ｘに対して単調増加となる。
【数７】
2013213712.tif 000009
上記（７）式におけるαはメタリック強度ｘを調整するための定数であり、例えば、α＝０．３８という値を用いる。
【００３１】
ステップＳ１６０では、明度補正項算出部１６が、上記測定装置５０による測定結果として得られる明度Ｌ*であって、特定角度（特定反射角）の明度Ｌ*に応じた補正項を算出する。当該ステップＳ１６０で算出する補正項は、第１補正項に該当する。ここでいう特定反射角とは、サンプルを観察するユーザーの観察角度として想定される値であって予め決められた値である。より詳しくは、特定反射角は以下のように決められる。つまり、複数の異なるサンプルをユーザーに自由に目視で観察させ、例えば、明るく感じられる順にサンプル間に順位を付けさせる。一方、測定装置５０によってもこれら複数のサンプルを測定し、反射角毎に、測定されたサンプル毎の明度Ｌ*の順位を付ける。そして、ユーザーが付けた順位と、測定されたサンプル毎の明度Ｌ*の順位との相関が最も高い反射角が、上記特定反射角として決定される。本実施形態では、このような相関は、反射角θ＝２５°〜３５°程度で高くなることが確認されたため、一例として、特定反射角を２５°としている。以下では、特定反射角の明度Ｌ*を、Ｌ*25と表記する。
【００３２】
このような特定反射角の明度Ｌ*25は、ユーザーの主観的評価に大きな影響を与える明度であると言える。また、対象サンプルの明るさに応じて、ユーザーが認識する「キラキラ度合い」も変動すると言える。そこで、本実施形態では、特定反射角の明度Ｌ*25に基づく補正項ｇk（Ｌ*25）を、空間的ノイズＬ*_Ｎに対する補正項として求めるようにしている。本実施形態では一例として、下記（８）式に示すようなｇk（Ｌ*25）を定義する。
【数８】
2013213712.tif 000010
上記（８）式におけるβkは定数であり、例えば、βk＝２．７２という値を用いる。
【００３３】
ステップ１７０では、評価指標算出部１９が、ここまでに算出された空間的ノイズＬ*_Ｎ、メタリック強度ｘに基づく補正項ｆk（ｘ）、および特定反射角の明度Ｌ*25基づく補正項ｇk（Ｌ*25）を掛け合わせて、第１評価指標Ｉｎ_１を算出する。第１評価指標Ｉｎ_１は、下記（９）式により表される。
【数９】
2013213712.tif 000011
【００３４】
ステップＳ１８０では、ピーク値補正項算出部１８が、上記ロジスティック分布を含む関数によりシミュレートした反射輝度の分布におけるピーク値ｒに応じた補正項を算出する。当該補正項は、第３補正項に該当し、第２評価指標の主要素たる正反射光の明度に対する補正項である。当該ピーク値ｒは、上記（３）式を参照すると、θ＝ｘpの場合であるため、係数ａを加味して、下記（１０）式で表すことができる。なおここでは、簡単のため、パラメーターｂ、ｐ（θ1）およびｃｏｓθは無視している。
【数１０】
2013213712.tif 000012
【００３５】
上述したように第２評価指標は、対象サンプルの「ピカピカ度合い」を評価するための指標であり、そのような光沢感は正反射の観察位置近傍で主に感じやすいものであると言える。そのため、本実施形態では、上記ロジスティック分布において正反射角で得られる上記ピーク値を、第２評価指標における補正項として採用し、ユーザーが実際に感じる「ピカピカ度合い」が第２評価指標として表されるようにしている。本実施形態では一例として、下記（１１）式に示すようなピーク値ｒに応じた補正項ｆp（ｒ）を定義する。
【数１１】
2013213712.tif 000013
上記（１１）式におけるγは定数であり、例えば、γ＝０．１２５という値を用いる。
【００３６】
ステップＳ１９０では、明度補正項算出部１６が、特定反射角の明度Ｌ*25に基づく補正項ｇp（Ｌ*25）を求める。上述したように、特定反射角の明度Ｌ*25は、ユーザーの主観的評価に大きな影響を与える明度であると言える。また、対象サンプルの明るさに応じて、ユーザーが認識する「ピカピカ度合い」も変動すると言える。そこで、本実施形態では、特定反射角の明度Ｌ*25に基づく補正項ｇp（Ｌ*25）を、第２評価指標の主要素たる正反射光の明度に対する補正項として求めるようにしている。ここで求める補正項は、第２補正項に該当する。本実施形態では一例として、下記（１２）式に示すようなｇp（Ｌ*25）を定義する。
【数１２】
2013213712.tif 000014
上記（１２）式におけるβpは定数であり、例えば、βp＝−０．３６という値を用いる。
【００３７】
ステップ２００では、評価指標算出部１９が、上記測定装置５０による測定結果として得られている正反射光の明度（反射角θ＝４５°の明度）Ｌ*45と、これまでに算出されたピーク値（ｒ）に応じた補正項ｆp（ｒ）および特定反射角の明度Ｌ*25基づく補正項ｇp（Ｌ*25）を掛け合わせて、第２評価指標Ｉｎ_２を算出する。第２評価指標Ｉｎ_２は、下記（１３）式により表される。
【数１３】
2013213712.tif 000015
なお、上記（１３）式の分母としての１００は、第２評価指標Ｉｎ_２のスケールを第１評価指標Ｉｎ_１に略合わせるために便宜的に入れた値であるため、このような分母は無くても良い（分母は１でもよい）し、１００以外の数値であってもよい。
【００３８】
ステップＳ２１０では、評価指標算出部１９は、上記算出した、第１評価指標Ｉｎ_１および第２評価指標Ｉｎ_２を、ディスプレー３０に表示させる処理を行なう。これら評価指標の表示態様は特に限られないが、例えば、ディスプレーにおいて、第１評価指標と第２評価指標とを各軸とする２次元平面に上記対象サンプルについて算出された第１評価指標Ｉｎ_１および第２評価指標Ｉｎ_２を表示させた画像を、出力させる。
【００３９】
図７は、ディスプレー３０に出力された画質評価確認画面３１の一例を示している。画質評価確認画面３１においては、縦軸を第１評価指標すなわち「キラキラ度合い」とし、横軸を第２評価指標すなわち「ピカピカ度合い」とした２次元平面が表されており、当該平面内における、上記算出された第１評価指標Ｉｎ_１および第２評価指標Ｉｎ_２に応じて特定される位置に、画質評価結果を表す特定形状の画像（例えば丸印Ｃ）がプロットされている。むろん画質評価確認画面３１は、このような表示態様に限られず、上記算出された第１評価指標Ｉｎ_１、第２評価指標Ｉｎ_２をそれぞれ別に表示してもよいし、これまで画質評価済みの複数の対象サンプルに関する第１評価指標Ｉｎ_１、第２評価指標Ｉｎ_２を同時に表示してもよい。また、画質評価確認画面３１においては、第１評価指標Ｉｎ_１、第２評価指標Ｉｎ_２を単に表示するだけでなく、第１評価指標Ｉｎ_１、第２評価指標Ｉｎ_２の各値に応じて予め用意された具体的な画質評価のコメント等を表示するとしてもよい。
【００４０】
このように本実施形態によれば、コンピューター１０は、少なくとも、対象サンプルの読取により得られた明度に基づいて算出される空間的ノイズＬ*_Ｎと、ユーザーの主観的評価との相関が高い特定反射角の明度Ｌ*25に応じた第１補正項ｇk（Ｌ*25）と、に基づいて定義される「キラキラ度合い」としての第１評価指標Ｉｎ_１を算出する。また、コンピューター１０は、少なくとも、対象サンプルにおける正反射光の明度Ｌ*45と、上記特定反射角の明度Ｌ*25に応じた第２補正項ｇp（Ｌ*25）と、に基づいて定義される第２評価指標Ｉｎ_２を算出する。従って、ユーザーが主観評価において感じる「キラキラ度合い」や「ピカピカ度合い」に極めて近い画質評価を、第１評価指標Ｉｎ_１および第２評価指標Ｉｎ_２の提示をもって実現することができる。
【００４１】
また、第１評価指標Ｉｎ_１は、対象サンプルの反射輝度の分布をシミュレートしたロジスティック分布を含む関数における複数のパラメーターの線形結合で表されたメタリック強度ｘに基づく補正項ｆk（ｘ）によって補正された上で提示される。そのため、実際にユーザーが感じると予測される「キラキラ度合い」に極めて近い画質評価を、第１評価指標Ｉｎ_１の提示をもって実現することができる。また、第２評価指標Ｉｎ_２は、上記ロジスティック分布を含む関数におけるピーク値ｒに基づく補正項ｆp（r）によって補正された上で提示される。そのため、実際にユーザーが感じると予測される「ピカピカ度合い」に極めて近い画質評価を、第２評価指標Ｉｎ_２の提示をもって実現することができる。
【００４２】
また本実施形態によれば、対象サンプルのメタリック感を、「キラキラ度合い」および「ピカピカ度合い」といった複数の評価指標に基づいて把握することができるため、対象サンプルについて総合的かつ効率的な画質評価を行なうことができる。
【００４３】
なお、上述した画質評価処理を実行する具体的製品は、コンピューター１０に限られない。例えば、プリンター２０が自機に搭載されたＣＰＵやメモリーといった資源を用いて内部のプログラムに従って、図２の処理の全部あるいは一部を実行し、また、プリンター２０が備える液晶パネル等のディスプレーに上記画質評価確認画面３１を表示したりしてもよい。
【符号の説明】
【００４４】
１…画質評価装置、１０…コンピューター、１１…ＨＤＤ、１１ａ…画質評価プログラム、１２…ＣＰＵ、１３…ＲＡＭ、１４…輝度シミュレーション部、１５…空間的ノイズ算出部、１６…明度補正項算出部、１７…メタリック強度補正項算出部、１８…ピーク値補正項算出部、１９…評価指標算出部、２０…プリンター、３０…ディスプレー、３１…画質評価確認画面、５０…測定装置、６０…スキャナー
【図１】
2013213712.tif 000016
【図２】
2013213712.tif 000017
【図３】
2013213712.tif 000018
【図４】
2013213712.tif 000019
【図５】
2013213712.tif 000020
【図６】
2013213712.tif 000021
【図７】
2013213712.tif 000022