印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出工程と、観察光源での基準彩度と前記観察光源下における正反射での前記印刷物の測色値の彩度との差分を算出する彩度差算出工程と、前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出工程と、を有することを特徴とするブロンジング指標値算出方法。前記指標値算出工程では、低彩度において前記ブロンジング指標値が低くなるよう重みを設定することを特徴とする請求項1に記載のブロンジング指標値算出方法。前記指標値算出工程では、前記色相角の差分と前記彩度の差分とを正規化することを特徴とする請求項1に記載のブロンジング指標値算出方法。前記ブロンジング指標値は、前記印刷物毎に取得されることを特徴とする請求項1~請求項3のいずれか一項に記載のブロンジング指標値算出方法。印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出手段と、観察光源下での基準彩度と前記観察光源下における正反射での前記印刷物の測色値の彩度との差分を算出する彩度差算出工程と、前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出手段と、を有することを特徴とするブロンジング指標値算出装置。印刷物のブロンジングを評価する値として、印刷物の色味を実際に測定して評価を行うことはあっても、印刷物のブロンジングの程度を定量的に示す指標値を算出することは開示されていない。本発明は、上記課題にかんがみてなされたもので、印刷物に対してそのブロンジングの程度を定量的に示す指標値であるブロンジング指標値の、算出方法の提供を目的とする。上記課題を解決するために、本発明では、印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出工程と、観察光源での基準彩度と前記観察光源下における正反射での前記記録剤の測色値の彩度との差分を算出する彩度差算出工程と、前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出工程と、を有する構成としてある。本発明にかかる具体的な構成及び作用は明細書中にて開示する。以下、下記の順序に従って本発明の実施形態を説明する。1.ブロンジング指標値算出方法1:2.ブロンジング指標値算出方法2:3.LUT作成方法:4.その他の実施形態:1.ブロンジング指標値算出方法1:以下、図を参照して、発明に係るブロンジング指標値算出方法を具体化した実施の形態について説明する。図1は、指標値取得システム100の構成を説明するためのブロック構成図である。指標値取得システム100は、コンピューター10と、プリンター20と、測色機2122を備えて構成されている。指標値取得システム100は、コンピューター10の統合的な制御のもと、測色機2122及びプリンター20を駆動させて、プリンタープリンター20に限定されないが使用するインク量Ijに対してブロンジング指標ブロンジング指標値BI、ブロンジング指標値BIIjとも記載する。を算出する。コンピューター10は、制御部11、記憶部12、USBインターフェイスI/F23、ビデオI/F31、入力I/F41などを備えている。そして、コンピューター10が備える制御部11が、記憶部12に記憶されたプログラムを読み込み、このプログラムに沿った演算を実行することにより、所定の処理を実行する。さらに、コンピューター10は、インク量Ijにおける指標値を算出する処理を行なうブロンジング指標値算出装置としても機能する。記憶部12は、例えばハードディスクドライブHDDであり、ブロンジングに係る指標値が記録された指標値テーブルTa等が記憶される。また、記録部12には、制御部11に、対してプリンター20の駆動制御、又は測色機21に測色値を取得させるための各種ドライバープログラムが記録されている。コンピューター10は、プログラムに従った演算を実行し、プリンター20をUSBI/F23等を介して制御することにより、印刷制御装置としても機能する。具体的にはコンピューター10は、印刷対象の画像データを取得し、当該画像データを色変換プロファイルによって画素単位で色変換し、色変換後の画像データにハーフトーン処理やマイクロウェーブ処理を行って印刷データを生成する。そして、当該印刷データをUSBIF23等を介してプリンター20に出力する。これによりプリンター20は、当該印刷データに基づく印刷を実行する。また、コンピューター10は、ビデオI/F31を介してディスプレイ30と接続されており、入力I/F41を介してキーボードやマウス等の操作部40と接続されている。測色機21は、測色対象物を測色してコンピューター10に対して測色値を供給する。測色機21は、キセノンランプ等により構成される光源例えば、標準のD50光源や、測色対象物を反射した光源の反射強度分光反射率とも記載する。を測定する測色部、更には測色部により取得された分光反射率を用いて測色値本実施形態では、Lab値を算出する測色値算出部を備えて構成される。ここで、測色機21の測色部は、測色対象物を拡散反射した光拡散反射光とも記載する。を受光して、その反射強度を取得するよう構成されている。測色機22は、測色機21同様、光源、測色部、更には測色値算出部を備えて構成され、コンピューター10に対して測色値を供給する。また、本実施形態では、測色機22の測色部は、測色対象物を正反射した光正反射光とも記載する。を受光して、その反射強度を取得するよう構成されている。プリンター20は、コンピューター10による制御のもと、メディア用紙にインクを噴射して画像を形成する。本実施形態ではシアンC、マゼンタM、イエローY、ブラックK、の4種類のインクを利用可能なカラープリンターを想定している。無論、プリンター20が使用するインクの種類は4種類に限定されず、これに淡シアンLc、淡マゼンダLm、淡ブラックLkを加えるものであってもよい。図2は、指標値取得システム100により実行される処理を示すフローチャートである。図2に示す処理により、コンピューター10は、各インク量Ijに対してブロンジング指標値BIを作成する。本処理により算出されるブロンジング指標値BIにより、あるインク量Ijに対してブロンジングの程度が指標化される。即ち、ブロンジング指標値BIが大きいインク量Ijほど、ブロンジングの程度が悪くなることを示し、値が小さいインク量Ijほど、ブロンジングの程度が良いことを示す。ここで、ブロンジング指標値は、その値が大きいほど程度が悪いものと定義している。利用者がコンピューター10の操作部40を操作して、図示しないUI画面において入力を行うと、コンピューター10は、ステップS1において、この入力を受付ける。また上記入力にて、ブロンジング指標値BIが作成される用紙の種別が選択される。例えば、このUI画面上において選択可能な用紙の一例として、普通紙、厚紙マット紙、コート紙光沢紙、特殊紙がある。ここで、ブロンジング指標値BIを作成する用紙を選択する理由は、発生するブロンジングの程度が用紙の種別によっても変化するからである。ステップS2では、コンピューター10は、プリンター20にカラーチャートを形成させる。図3は、プリンター20が形成するカラーチャートを示す図である。本実施形態に係るカラーチャートは、複数のパッチiが縦横に配列して構成される。また、各パッチiは、それぞれプリンター20が噴射する1次色インクC、M、Y、Kを、そのインク量Ij化させて所定の階調値が形成されている。なお、各パッチiを形成する1次色インクの階調値は、このインクを搭載するプリンターが使用する階調値に応じて設定される。なお、本実施形態では、インク量Ijを1次色インクにより構成するが、1次色インクの組合せ等とするものであってもよい。また、本明細書においてカラーチャートという用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。ステップS3では、コンピューター10は、測色機21にカラーチャートを測色させて、各パッチiの測色値を取得する。図4は、測色機におけるカラーチャートの測色方法を示す図である。ここで、図4Aは、測色機21における測色方法を示す。また、図4Bは、測色機22における測色方法を示す。図4Aに示すように、測色機21は、光源をカラーチャートに対して垂直0度に配置し、測色部をカラーチャートに対して45度に配置し、各パッチiを測色する。そのため、測色機21により、各パッチiに対して拡散反射での測色値が取得される。そして、測色機21は、カラーチャートを形成するすべてのパッチiに対して測色値を取得する。ステップS4では、コンピューター10は、測色機22にカラーチャートを測色させて各パッチiの測色値を取得する。図4Bに示すように、測色機21は、光源をカラーチャートに対して45度に配置し、測色部をカラーチャートに対して-45度に配置し、各パッチiを測色する。そのため、測色機21により、各パッチiに対して正反射での測色値が取得される。ステップS5では、コンピューター10は、ステップS3及びS4で取得した各測色値からブロンジング指標値BIを取得する。一般にブロンジングとして認識される現象として以下の2点が挙げられる。現象1:観察光源に対して正反射である角度で観察した色と、拡散反射で観察した色との、色味色相が異なる。現象2:観察光源に対して正反射である角度で観察した色と、観察光源の色との、色味色相が異なる。そのため、本実施形態では、上記現象1及び現象2を定量化し、ブロンジング指標値BIを算出する。図5は、ステップS5の処理を詳細に示すフローチャートである。また、図6は、ブロンジング指標値を説明するための図である。図5のステップS51では、コンピューター10は、拡散反射における色相角Huestdと、正反射における色相角Huespecularとの差分となるΔHueを算出する。即ち、本ステップでは、現象1として示した、正反射で観察した色味と、通常の観察拡散反射での色味の違いを、色相角Hueの差として求める。図6Aにおいて、色相角は、色相環上の角度として算出される。ここで、正反射における色味と拡散反射における色味が異なる場合、色相環上の色相角Hueも異なるため、両色味の差が大きくなるほど、差となる角度も大きくなる。そのため、本ステップでは、現象1を色相角Hueの差として定量化する。そのため、コンピューター10は、ステップS3により取得した拡散反射におけるパッチiの測色値Lab値から色相角Huestdを算出する。また、コンピューター10は、ステップS4により取得した正反射におけるパッチiの測色値Lab値から色相角Huespecularを算出する。一般に、色相角Hueは、下記式1により求めることができる。ここで、aはbは、測色値Lab値におけるa、bの値である。そうして求めた角色相角Huestd、Huespecularもとに下記式2を用いて、パッチiインク階調値毎の色相角の差ΔHueを算出する。ただし、上記ΔHueの値が取りうる範囲は、0≦ΔHue≦2πRadであるものとする。上記ステップS51の処理により、本発明の差分算出工程、差分算出手段を実現する。ステップS52では、コンピューター10は、正反射における基準彩度からの差ΔChromaを算出する。即ち、本ステップでは、現象2として示した、観察光源の色味に対して、正反射で異なる色味が見える現象を、正反射における彩度の差ΔChromaとして定量化する。図6Bに示すように、同一観察光源で正反射での色味が異なる場合、色相環上で観察光源下での彩度基準彩度とも記載する。に対して、正反射で対象物を観察した彩度Chromaspecularが異なる。ここで、彩度Chromaは、一般には、下記式3により算出することができる。本実施形態では、測色機1及び測色機2において観察光源を白色光とするため、基準彩度は限りなく小さい値≒0となる。そのため、正反射における彩度を彩度の差ΔChromaとして使用することができる。上記ステップS52の処理により、本発明の彩度差算出工程、彩度差算出手段を実現する。ステップS53では、コンピューター10は、現象1現象2を定量化した各値ΔHue、ΔChromaを乗算してブロンジング指標値BIIjを算出する。即ち、ブロンジング指標値BIIjは、下記式4により算出することができる。ここで、weightは、低彩度においてブロンジング指標値BIIjが低くなるよう調整する重みである。本実施形態では、重みweightは、下記式5に示すようにシグモイド関数として算出される。図7は、本実施形態に係る重みweightを示す図である。重みweightは、シグモイド関数であるため、拡散反射での彩度Chromaspecularの値に応じて単調増加し、1つの変曲点をもつ。即ち、低彩度域においては、重みweightが限りなく小さく、ある彩度Chromamaxを境に重みweightが増加し、やがて重みweightが1付近に飽和する。そのため、この重みweightをブロンジング指標値BIIjの一部として付与することにより、低彩度域ではブロンジング指標値BIIjの値が小さくなり即ち、指標値として良い値となる。、高彩度域ではブロンジング指標値BIIjの値が大きくなる即ち、指標値として悪い値となる。。ブロンジング指標値BIIjを彩度に応じたシグモイド関数により重み付けする理由として、色相角の差ΔHueと彩度の差ΔChromaとはその取りうる値の範囲が異なるため、影響度合いを均一にする必要があるからである。即ち、色相角の差ΔHueは、0≦ΔHue≦2πなお、単位はradの範囲での値であるのに対して、彩度の差ΔChromaは、0≦Chroma≦50の範囲での値となる。そのため、単純に両者ΔHue、ΔChromaを乗算すると、値が小さい色相角の差ΔHueが支配的となってしまう。また、ブロンジング指標値BIIjは、彩度が低彩度であるほど、その値が良くなる指標値が低くなるよう設定される必要がある。即ち、ΔChromaの算出において正反射での彩度Chromaspecularが小さくなるほどブロンジングが発生し難くなるよう値を調整する必要がある。そのため、低彩度ほど、ブロンジング指標値BIIjが低い値を取るよう、重みweightを付与している。上記ステップS53の処理により、本発明の指標値算出工程、指標値算出手段を実現する。図2に戻り、ステップS6では、コンピューター10は、ステップS5で取得した指標値を指標値テーブルTaに記録する。図8は、指標値テーブルTaを説明する図である。図8に示すように、指標値テーブルTaには、プリンターが使用する用紙毎本実施形態では、普通紙、厚紙、コート紙、特殊紙に、ブロンジング指標値BIIjが設定される。即ち、用紙毎に、インク量Ij階調値と、このインク量Ijに対応するブロンジング指標値BIIjが記録される。そのため、後述するプリンター20におけるLUTの作成において、この指標値テーブルTaが参照され、最適なインク量Ijの探索が行われる。無論、指標値テーブルTaに記録される指標値は用紙毎の形式に限定されるものではなく、コンピューター10における印刷制御情報例えば、印刷時に噴射されるインク滴重量の情報などを加えてもよい。2.ブロンジング指標値算出方法2:ブロンジング指標値BIを算出する過程において、色相角の差ΔHueと彩度の差ΔChromaをそれぞれ正規化してもよい。上記のように、色相角の差ΔHueは、彩度の差ΔChromaに比べて値が小さいため、単純に乗算するとブロンジング指標値BIにおいて、色相角の差ΔHueが支配的になってしまう。そのため、色相角の差ΔHueと彩度の差ΔChromaをそれぞれ0以上、1以下の値に正規化し、正規化後の値によりブロンジング指標値BIIjを算出してもよい。即ち、図5のステップS53において、コンピューター10は、正規化後の色相角の差ΔHueと彩度の差ΔChromaを用いてブロンジング指標値BIIjを算出する。そして、図2のステップS6において、コンピューター10は、算出されたブロンジング指標値BIIjを指標値テーブルTaに記録する。以上説明したように、本ブロンジング指標値算出方法により、所定印刷条件下での印刷物におけるブロンジングの程度を指標値として算出するため、印刷物に発生しているブロンジングの程度の評価等を定量的に扱うことができる。3.LUT作成方法:以下、本指標値作成方法により作成されたブロンジング指標値BIを用いたLUT作成方法を説明する。図9は、LUT作成システム110を示すブロック構成図である。本LUT作成システム110は、コンピューター10と、プリンター20と、測色機不図示により実現される。LUT作成方法に係るコンピューター10の記録部12には、プログラムや色変換プロファイル例えば、後述するインバースモデル初期LUT410、色補正なしLUT510、色補正LUT610、ブロンジングに係る指標値が記録された指標値テーブルTa等が記憶される。なお、LUTは、ルックアップテーブルの略語である。色補正なしLUT510は、所定の入力色空間例えばRGB色空間の階調値をプリンター20で使用される複数種類のインクの組み合わせインク量Ijに変換するための色変換テーブルである。色補正なしLUT510の入力色空間であるRGB色空間は、いわゆる機器依存色空間では無く、特定のデバイスとは無関係に設定された仮想の色空間あるいは抽象的な色空間である。この色補正なしLUT510は、例えば色補正LUT610を作成する際に使用される。色補正LUT610は、標準的な機器依存の色空間例えばsRGB色空間やJAPANCOLOR2001色空間を、特定のプリンターのインク量Ijに変換するためのルックアップテーブルである。インバースモデル初期LUT410については後述する。また、コンピューター10のLUT作成モジュールは、初期値設定部14と、スムージング処理部15と、テーブル作成部16と、ガマット予測部17と、フォワードモデルコンバーター18と、を有している。これらの各部の機能については後述する。<<プロファイルの作成について>>図10は、色補正なしLUT作成時の全体処理手順を示すフロー図である。ステップS100では、インクデータの設定が行われる。ここで、インクデータとは、プリンター20で使用されるインクの色や種類、デューティ制限値などである。また、デューティ制限値とは、メディア用紙の単位面積の中に、打ち込むことのできるインク量Ijを示すものである。このデューティ制限値は、メディアの種類や形成する画像の解像度に依存する。このため、デューティ制限値は、メディアの種類毎および印刷解像度毎に定められる。本実施形態では、写真用紙に1440×1440dpiの解像度で印刷を行うこととする。また、デューティ制限値として、インク単色によるデューティ制限値や、複数色の合計によるデューティ制限値合計デューティ制限値がある。例えば、ある紙のある解像度におけるデューティ制限値として、単色では80、各色の合計では120のように定められる。ステップS200では、カラーチャートの印刷が行われる。印刷されるカラーチャートは、図3に示すものと同様であり、複数のパッチiが配列して構成されたものである。また、各パッチiは、所定のインク量Ijを用いて形成されている。ここではデューティ制限値や色変換テーブルを用いずに、簡単なインク量Ijの組み合わせによりカラーチャートの印刷を行なう。ステップS300では、例えば、不図示の測色機やスキャナなどを用いて、カラーチャートを構成する各パッチの測色が行われる。この測色結果は、コンピューター10に取り込まれ、コンピューター10においてインク量Ijと測色値Lab値とを対応させた仮想サンプルが形成される。そして、ステップS400では、カラーチャートの測色結果仮想サンプルに基づいてフォワードモデル以下FMともいうの作成が行われる。なお、フォワードモデルとは、インク量Ijを機器独立色空間の色彩値測色値に変換する変換モデルを意味する。逆に、インバースモデルとは、機器独立色空間の色彩値をインク量Ijに変換する変換モデルを意味する。本実施形態では、機器独立色空間としてCIELab色空間を使用する。なお、以下では、CIELab色空間の色彩値を、単にLab値又はLab値とも呼ぶ。図11は、FMコンバーターとインバースモデル初期LUTの説明図である。図11に示すように、FMコンバーター18の前段を構成する分光プリンティングコンバーターRCは、複数種類のインクのインク量Ijを、そのインク量Ijに応じて印刷されるパッチの分光反射率Rλに変換する。本実施形態では、前述した4色のインクを使用するプリンター20を想定しており、分光プリンティングモデルコンバーターRCもこの4種類のインクのインク量Ijを入力としている。但し、プリンター20で使用する複数種類のインクとしては、任意のインク量Ijを利用することが可能である。色算出部CCは、分光反射率RλからLab色空間の色彩値を算出する。この色彩値の算出には、予め選択された光源例えば標準の光D50がカラーチャートの観察条件として使用される。なお、分光プリンティングモデルコンバーターRCを作成する方法としては、例えば特表2007-511175号公報に記載された方法を採用することが可能である。インバースモデル初期LUT410は、Lab値を入力とし、インク量Ijを出力とするルックアップテーブルである。この初期LUT410は、例えば、Lab色空間を複数の小セルに区分し、小セル毎に最適なインク量Ijを選択して登録したものである。この選択は、例えば、そのインク量Ijで印刷されるカラーチャートの画質を考慮して行われる。一般に、或る1つのLab値を再現するインク量Ijの組み合わせは多数存在する。そこで、初期LUT410では、ほぼ同じLab値を再現する多数のインク量Ijの組み合わせの中から、画質等の所望の観点から最適なインク量Ijを選択したものが登録されている。この初期LUT410の入力値であるLab値は各小セルの代表値である。一方、出力値であるインク量Ijはそのセル内のいずれかのLab値を再現するものである。従って、この初期LUT410では、入力値であるLabと出力値であるインク量Ijとが厳密に対応したものとなっておらず、出力値のインク量IjをFMコンバーター18でLab値に変換すると、初期LUT410の入力値とは多少異なる値が得られる。但し、初期LUT410として、入力値と出力値とが完全に対応するものを利用してもよい。また、初期LUT410を用いずに色補正なしLUT510を作成することも可能である。なお、小セル毎に最適なインク量Ijを選択して初期LUT410を作成する方法としては、例えば上記特表2007-511175号公報に記載された方法を採用することが可能である。図10のステップS500では、フォワードモデルFMとインクのデューティ制限値を用いてガマットの予測が行われる。ステップS600では、ステップS500で予測されたガマット以下、予測ガマットともいうの評価が行われる。なお、この判断は、例えばユーザーが行う。予測ガマットの評価がNGであればS600でNO、ステップS100に戻る。一方、予測ガマットの評価がOKであればS600:YES、ステップS700で、色補正なしLUT510作成のための初期入力値がユーザーによって設定される。色補正なしLUT510の入力値としては、RGBの各値として予め定められたほぼ等間隔の値が設定される。1組のRGB値はRGB色空間内の点を表していると考えられるので、1組のRGB値を入力格子点とも呼ぶ。ステップS700においては、複数の入力格子点のうちから予め選択された幾つかの少数の入力格子点に対するインク量Ijの初期値がユーザーによって入力される。この初期入力値が設定される入力格子点としては、RGB色空間における3次元色立体の頂点に相当する入力格子点を少なくとも選択することが好ましい。この3次元色立体の頂点では、RGBの各値がその定義範囲の最小値又は最大値を取る。具体的には、RGBの各値を8ビットで表現した場合には、、、、、、、、である8つの入力格子点に関してインク量Ijの初期入力値が設定される。図10のステップS800では、スムージング処理部15図9が、ステップS700で設定された初期入力値に基づいてスムージング処理平滑化処理を実行する。図12は、スムージング処理における入力格子点とLab色空間の座標値との対応関係を示す図である。なお、図12Bにはスムージング処理前の状態における複数の色点の分布が2重丸と白丸とで示されている。これらの座標値は、Lab色空間における3次元色立体CSを構成している。なお、各座標値のLab値は、色補正なしLUT510の複数の入力格子点におけるインク量Ijを、FMコンバーター18を用いてLab値に変換した値である。なお、他の入力格子点に対するインク量Ijの初期値は、初期入力値から初期値設定部14図9参照によって設定される。また、図12Cには、スムージング処理後の色点の分布が示されている。スムージング処理は、Lab色空間における複数の座標値を移動させて、それらの座標値の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各座標値のLab値を再現するために最適なインク量Ijも決定される。このとき、入力格子点に割り当てられる最適なインク量Ijを目的関数Eをもとに算出する。Lab色空間の3次元色立体CS図12B、図12Cの上記各頂点は、色補正なしLUT510の入力色空間RGB色空間の3次元色立体図12Aの頂点と一対一に対応している。また、各頂点を結ぶ辺稜線も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前図12BのLab色空間の各座標値は、色補正なしLUT510の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後図12CのLab色空間の各色点も色補正なしLUT510の入力格子点にそれぞれ対応付けられる。なお、色補正なしLUT510の入力格子点はスムージング処理によって変化しない。色補正なしLUT510を作成する際に、Lab空間においてスムージング処理を行う理由は以下の通りである。色補正なしLUT510では、なるべく大きな色域を再現できるように出力色空間のインク量Ijを設定したいという要望がある。特定のインクセットで再現可能な色域は、インクデューティ制限値一定面積に吐出可能なインク量Ijの制限値などの所定の制限条件を考慮して決定される。一方、上述したフォワードモデルFMは、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にインクデューティ制限等の制限条件を考慮してLab空間内の座標値の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、座標値の移動を行うアルゴリズムとしては、例えば、力学モデルを使用したものが利用される力学モデルについては、例えば特開2006-197080号公報参照。図13は、スムージング処理図10のステップS800の処理手順の一例を示すフロー図である。ステップS810では、初期値設定部14図9が、スムージング処理の対象とする複数の座標値を初期設定する。具体的には、まず、インク量Ijの初期入力値から、スムージング処理の対象となる各座標値の仮インク量Ijが決定される。次に、FMコンバーター18を用いて、仮インク量Ijに対応するLab値を求める。そして、得られたLab値を、インバースモデル初期LUT410図9を用いてインク量Ijに再度変換する。上述のステップS810の処理の結果、スムージング処理の対象となる座標値について、以下の初期値が決定される。i色補正なしLUTの入力格子点の値:ii各入力格子点に対応するLab色空間の座標値の初期座標値:Liii各入力格子点に対応する初期インク量Ij:I以上の説明から理解できるように、初期値設定部14は、代表的な入力格子点に関する入力初期値から他の入力格子点に関する初期値を設定する機能を有している。なお、初期値設定部14は、スムージング処理部15に含まれるものとしてもよい。ステップS820では、スムージング処理部15は、Lab空間内の座標値を移動させる。図14は、図13のステップS820~S850の処理内容を示す説明図である。図14Aに示すように、スムージング処理前には、座標値の分布にはかなりの偏りがある。図14Bは、微少時間経過後の各色点の位置を示している。この移動後のLab値をターゲット値LABtと呼ぶ。ターゲットという修飾語は、この値LABtが、以下で説明するインク量Ijの最適値の探索処理の際の目標値として使用されるからである。ステップS830では、スムージング処理部15は、以下の式6で表される目的関数Eを用いて、ターゲット値LABtに対するインク量Ijの最適値を探索する図14C参照。ここで、Ijは目的関数Eを最小とするインク量Ij、LABFMはインク量IjをFMコンバーター18で変換して得られるLab値、GIIjはインク量Ijで印刷されるカラーチャートの粒状性指数、BIIjはブロンジング指標、α、βは定数である。6式の右辺第1項は、2つのベクトルLABtLABFMで表される2つの座標値の距離の2乗を求める演算を意味している。上記式6は、ステップS820の処理において微少量だけ移動した後の色点の座標値LABtに近いLab値を再現するインク量Ijの中で、粒状性指数GI及びブロンジング指標値BIのより小さいインク量Ijが最適なインク量Ijとして決定されることを意味している。この最適なインク量Ijの探索は、例えば準ニュートン法などの最適化手法を用いて実行される。また、インク量Ijの探索は、ステップS810で設定された各入力格子点の初期インク値から開始される。従って、探索で得られるインク量Ijは、この初期インク量Ijを修正した値となる。なお、一般に、最適化処理によるインク量Ijの探索は以下の条件の下で実行される。条件i目的関数Eを最小とする条件ii使用可能として予め指定されたインクのみを使用する条件iの目的関数Eを与える式としては、上記6式以外の種々の式を利用可能である。また、条件iiにおける使用可能なインクとしては、プリンター20で利用可能な複数種類のインクのうちの一部のインクのみが指定される場合と、全種類のインクが指定される場合とがある。使用可能なインクの種類は、スムージング処理の対象となる色点毎にユーザーによって予め設定されることが好ましい。これらの最適化条件の具体例については更に後述する。ここで、ブロンジング指標値BIは、図8に示す指標値テーブルTaに記録された各値が取得される。このとき、カラーチャートが形成されたメディアの種別に応じて、インク量Ijに対応するブロンジング指標値BIIjが選択される。また、粒状性指数GIIjは、例えば以下の式7で与えられる。ここで、aLは明度補正係数、WSuはカラーチャートの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、VTFuは視覚の空間周波数特性、uは空間周波数である。粒状性指数GIの計算方法としては、例えば、特表2007-511161に記載されている。図13のステップS840では、ステップS830で探索されたインク量Ijに対応するLab値が、FMコンバーター18で再計算される図14D。ここで、Lab値を再計算する理由は、探索されたインク量Ijが目的関数Eを最小とするインク量Ijなので、そのインク量Ijで再現されるLab値は、最適化処理のターゲット値LABtから多少ずれているからである。こうして再計算されたLab値が、各色点の移動後の座標値として使用される。ステップS850では、各色点の座標値の移動量の平均値aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値aveが閾値εよりも大きい場合には、ステップS820に戻りステップS820~S850のスムージング処理が継続される。一方、移動量の平均値aveが閾値ε以下の場合には、色点の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。このように、本実施形態のスムージング処理平滑化処理では、各色点を微少時間間隔毎に移動させつつ、移動後の色点に対応する最適なインク量Ijを最適化手法で探索する。そして、色点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図14Cに示したように、平滑な色点分布を得ることが可能である。図10に戻り、ステップS900では、スムージング処理の結果を用いて、テーブル作成部16が色補正なしLUT510を作成する。すなわち、テーブル作成部16は、各入力格子点に対応付けられたLab色空間の色点を再現するための最適なインク量Ijを色補正なしLUT510の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、色補正なしLUT510の入力格子点の一部のみに対応する色点のみを処理対象として選択することも可能である。例えば、色補正なしLUT510の入力格子点におけるRGB値の間隔が16である場合に、スムージング処理の対象となる入力格子点におけるRGB値の間隔を32に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成部16は、スムージング処理結果を補間することによって色補正なしLUT510のすべての入力格子点に対するインク量Ijを決定して登録する。そして、図10のステップS1000では、色補正なしLUT510の評価が行われる。評価がNGであればS1000でNO、再度ステップS100に戻る。一方、評価がOKであれば、色補正なしLUT510を記憶部12に登録する。なお、LUT作成モジュールは、色補正なしLUT510を用いて、色補正LUT610の作成も行う。図15は、色補正なしLUTを用いた色補正LUTの作成方法を示す説明図である。上記のように色補正LUT610は、標準的な機器依存の色空間を、特定のプリンターのインク量Ijに変換するためのルックアップテーブルであり、本実施形態では、機器依存の色空間としてsRGB色空間を採用する。以下、図15に示すように、色補正なしLUTが規定するRGB値とインク量Ijの関係をもとに、色補正LUT610の作成が行われる。まず、色補正なしLUT510により、RGB値がインク量Ijに変換される。変換後のインク量Ijは、FMコンバーター18によってLab値に変換される。一方、色補正LUT610に対応付けられるsRGB値は、既知の変換式に従ってLab値に変換される。そして、変換後のLab値は、その色域が、FMコンバーター18で変換されたLab値の色域と一致するようにガマットマッピングされる。一方、色補正なしLUT510とFMコンバーター18を通じて、RGB値から変換したLab値を、逆方向ルックアップテーブルとして、逆変換LUT511が作成される。即ち、逆変換LUT511は、Lab値とRGB値とを対応付けるLUTである。さらに、ガマットマッピングされたLab値は、この逆変換LUT511によってRGB値に変換される。このRGB値は、さらに、色補正なしLUT510によってインク量Ijに再度変換される。そして、最後のインク量Ijと最初のsRGB値の対応関係をルックアップテーブルに登録することによって、色補正LUT610が作成される。以上により、機器依存の色空間とインク量Ijとを対応付ける色補正LUT610が作成され、記録部12に記録される。以後、コンピューター10では、この色補正LUT610を使用して入力データに対応するインク量Ijを選択してプリンター20に出力する。また、本実施形態では、色補正LUT610の元となる色補正なしLUT510の作成において、最適なインク量Ijの評価項目としてブロンジング指標値BIIjが用いられる。そのため、色補正LUT610に記録されるインク量Ijは、ブロンジングを抑制するようその組合せが選択される。4.その他の実施形態:記録剤としてインクを用いることは一例である。例えば、プリンター20としてレザープリンターを用いる場合は、トナーのような記録剤に対して本指標値算出方法を適用するものであってもよい。また、ブロンジング指標値の適用例としてLUT作成方法を用いたことは一例であり、これ以外にも様々な適用例がある。なお、本発明は上記実施例に限られるものでないことは言うまでもない。即ち、上記実施例の中で開示した相互に置換可能な部材および構成等を適宜その組み合わせを変更して適用すること、上記実施例の中で開示されていないが、公知技術であって上記実施例の中で開示した部材および構成等と相互に置換可能な部材および構成等を適宜置換し、またその組み合わせを変更して適用すること、上記実施例の中で開示されていないが、公知技術等に基づいて当業者が上記実施例の中で開示した部材および構成等の代用として想定し得る部材および構成等と適宜置換し、またその組み合わせを変更して適用すること、は本発明の一実施例として開示されるものである。