(57)【特許請求の範囲】
【請求項１】
印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出工程と、
観察光源での基準彩度と前記観察光源下における正反射での前記印刷物の測色値の彩度との差分を算出する彩度差算出工程と、
前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出工程と、を有することを特徴とするブロンジング指標値算出方法。
【請求項２】
前記指標値算出工程では、低彩度において前記ブロンジング指標値が低くなるよう重みを設定することを特徴とする請求項１に記載のブロンジング指標値算出方法。
【請求項３】
前記指標値算出工程では、前記色相角の差分と前記彩度の差分とを正規化することを特徴とする請求項１に記載のブロンジング指標値算出方法。
【請求項４】
前記ブロンジング指標値は、前記印刷物毎に取得されることを特徴とする請求項１〜請求項３のいずれか一項に記載のブロンジング指標値算出方法。
【請求項５】
印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出手段と、
観察光源下での基準彩度と前記観察光源下における正反射での前記印刷物の測色値の彩度との差分を算出する彩度差算出工程と、
前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出手段と、を有することを特徴とするブロンジング指標値算出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ブロンジングを評価するための指標化に関する。
【背景技術】
【０００２】
従来、印刷物の光沢性を強調するインク等の記録剤や、用紙等のメディアが知られている。これらの記録剤やメディアは、特定の用途（例えば、写真印刷等の写像性を高めたいもの。）において用いられ、印刷物の質感を際立たせる一因となっている（例えば、特許文献１−３参照）。
【０００３】
上記した光沢性に起因して発生する現象として、ブロンジング（又はブロンズ現象ともいう。）が知られている。ブロンジングは、印刷物が光源に対して正反射している状態から観察した際に、通常の拡散反射している状態の観察とは異なる色味がついて見え、あたかも金属光沢のように見える現象である。ブロンジングが生じている印刷物では、印刷物と観察者との相対的な位置関係で印刷物の色味が変化してしまうため、好ましくない。また、ブロンジングの要因は、インク等の記録剤の組合せや、記録剤とこの記録剤が記録されるメディアの組合せ等、様々な要因により発生する。
【０００４】
そのため、印刷物のブロンジング特性をその色味（厳密にはＬ＊ａ＊ｂ＊値）をもとに算出し、最適なインクの組合せを選択する際に、この算出値を用いる技術が開示されている（例えば、特許文献４参照。）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−１９４５００号公報
【特許文献２】特開２００５−２７４３４０号公報
【特許文献３】特開２０１０−５２２４７号公報
【特許文献４】特開２００７−２８８２７０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
印刷物のブロンジングを評価する値として、印刷物の色味を実際に測定して評価を行うことはあっても、印刷物のブロンジングの程度を定量的に示す指標値を算出することは開示されていない。
【０００７】
本発明は、上記課題にかんがみてなされたもので、印刷物に対してそのブロンジングの程度を定量的に示す指標値であるブロンジング指標値の、算出方法の提供を目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決するために、本発明では、印刷物に対して拡散反射での測色値の色相角と正反射での測色値の色相角の差分を算出する差分算出工程と、観察光源での基準彩度と前記観察光源下における正反射での前記記録剤の測色値の彩度との差分を算出する彩度差算出工程と、前記算出された色相角の差分と前記彩度の差分とを乗算した値をもとにブロンジング指標値を算出する指標値算出工程と、を有する構成としてある。
【０００９】
本発明にかかる具体的な構成及び作用は明細書中にて開示する。
【図面の簡単な説明】
【００１０】
【図１】指標値取得システム１００の構成を説明するためのブロック構成図である。
【図２】指標値取得システム１００により実行される処理を示すフローチャートである。
【図３】プリンター２０が形成するカラーチャートを示す図である。
【図４】測色機におけるカラーチャートの測色方法を示す図である。
【図５】ステップＳ５の処理を詳細に示すフローチャートである。
【図６】ブロンジング指標値を説明するための図である。
【図７】本実施形態に係る重みｗｅｉｇｈｔを示す図である。
【図８】指標値テーブルＴを説明する図である。
【図９】ＬＵＴ作成システム１１０を示すブロック構成図である。
【図１０】色補正なしＬＵＴ作成時の全体処理手順を示すフロー図である。
【図１１】ＦＭコンバーターとインバースモデル初期ＬＵＴの説明図である。
【図１２】スムージング処理における入力格子点とＬ＊ａ＊ｂ＊色空間の座標値との対応関係を示す図である。
【図１３】スムージング処理（図１０のステップＳ８００）の処理手順の一例を示すフロー図である。
【図１４】図１３のステップＳ８２０〜Ｓ８５０の処理内容を示す説明図である。
【図１５】色補正なしＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。
【発明を実施するための形態】
【００１１】
以下、下記の順序に従って本発明の実施形態を説明する。
１．ブロンジング指標値算出方法１：
２．ブロンジング指標値算出方法２：
３．ＬＵＴ作成方法：
４．その他の実施形態：
【００１２】
１．ブロンジング指標値算出方法１：
以下、図を参照して、発明に係るブロンジング指標値算出方法を具体化した実施の形態について説明する。図１は、指標値取得システム１００の構成を説明するためのブロック構成図である。
【００１３】
指標値取得システム１００は、コンピューター１０と、プリンター２０と、測色機２１,２２を備えて構成されている。指標値取得システム１００は、コンピューター１０の統合的な制御のもと、測色機２１,２２及びプリンター２０を駆動させて、プリンター（プリンター２０に限定されない）が使用するインク量Ｉｊに対してブロンジング指標（ブロンジング指標値ＢＩ、ブロンジング指標値ＢＩ（Ｉｊ）とも記載する。）を算出する。
【００１４】
コンピューター１０は、制御部１１、記憶部１２、ＵＳＢインターフェイス（Ｉ／Ｆ）２３、ビデオＩ／Ｆ３１、入力Ｉ／Ｆ４１などを備えている。そして、コンピューター１０が備える制御部１１が、記憶部１２に記憶されたプログラムを読み込み、このプログラムに沿った演算を実行することにより、所定の処理を実行する。さらに、コンピューター１０は、インク量Ｉｊにおける指標値を算出する処理を行なうブロンジング指標値算出装置としても機能する。
【００１５】
記憶部１２は、例えばハードディスクドライブ（ＨＤＤ）であり、ブロンジングに係る指標値が記録された指標値テーブルＴａ等が記憶される。また、記録部１２には、制御部１１に、対してプリンター２０の駆動制御、又は測色機２１に測色値を取得させるための各種ドライバープログラムが記録されている。
【００１６】
コンピューター１０は、プログラムに従った演算を実行し、プリンター２０をＵＳＢＩ／Ｆ２３等を介して制御することにより、印刷制御装置としても機能する。具体的にはコンピューター１０は、印刷対象の画像データを取得し、当該画像データを色変換プロファイルによって画素単位で色変換し、色変換後の画像データにハーフトーン処理やマイクロウェーブ処理を行って印刷データを生成する。そして、当該印刷データをＵＳＢＩ/Ｆ２３等を介してプリンター２０に出力する。これによりプリンター２０は、当該印刷データに基づく印刷を実行する。
また、コンピューター１０は、ビデオＩ／Ｆ３１を介してディスプレイ３０と接続されており、入力Ｉ／Ｆ４１を介してキーボードやマウス等の操作部４０と接続されている。
【００１７】
測色機２１は、測色対象物を測色してコンピューター１０に対して測色値を供給する。測色機２１は、キセノンランプ等により構成される光源（例えば、標準のＤ５０光源）や、測色対象物を反射した光源の反射強度（分光反射率とも記載する。）を測定する測色部、更には測色部により取得された分光反射率を用いて測色値（本実施形態では、Ｌ＊ａ＊ｂ＊値）を算出する測色値算出部を備えて構成される。ここで、測色機２１の測色部は、測色対象物を拡散反射した光（拡散反射光とも記載する。）を受光して、その反射強度を取得するよう構成されている。
【００１８】
測色機２２は、測色機２１同様、光源、測色部、更には測色値算出部を備えて構成され、コンピューター１０に対して測色値を供給する。また、本実施形態では、測色機２２の測色部は、測色対象物を正反射した光（正反射光とも記載する。）を受光して、その反射強度を取得するよう構成されている。
【００１９】
プリンター２０は、コンピューター１０による制御のもと、メディア（用紙）にインクを噴射して画像を形成する。本実施形態ではシアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、の４種類のインクを利用可能なカラープリンターを想定している。無論、プリンター２０が使用するインクの種類は４種類に限定されず、これに淡シアン（Ｌｃ）、淡マゼンダ（Ｌｍ）、淡ブラック（Ｌｋ）を加えるものであってもよい。
【００２０】
図２は、指標値取得システム１００により実行される処理を示すフローチャートである。図２に示す処理により、コンピューター１０は、各インク量Ｉｊに対してブロンジング指標値ＢＩを作成する。本処理により算出されるブロンジング指標値ＢＩにより、あるインク量Ｉｊに対してブロンジングの程度が指標化される。即ち、ブロンジング指標値ＢＩが大きいインク量Ｉｊほど、ブロンジングの程度が悪くなることを示し、値が小さいインク量Ｉｊほど、ブロンジングの程度が良いことを示す。ここで、ブロンジング指標値は、その値が大きいほど程度が悪いものと定義している。
【００２１】
利用者がコンピューター１０の操作部４０を操作して、図示しないＵＩ画面において入力を行うと、コンピューター１０は、ステップＳ１において、この入力を受付ける。また上記入力にて、ブロンジング指標値ＢＩが作成される用紙の種別が選択される。例えば、このＵＩ画面上において選択可能な用紙の一例として、「普通紙」、「厚紙（マット紙）」、「コート紙（光沢紙）」、「特殊紙」がある。
ここで、ブロンジング指標値ＢＩを作成する用紙を選択する理由は、発生するブロンジングの程度が用紙の種別によっても変化するからである。
【００２２】
ステップＳ２では、コンピューター１０は、プリンター２０にカラーチャートを形成させる。図３は、プリンター２０が形成するカラーチャートを示す図である。本実施形態に係るカラーチャートは、複数のパッチ（ｉ）が縦横に配列して構成される。また、各パッチ（ｉ）は、それぞれプリンター２０が噴射する１次色インク（Ｃ、Ｍ、Ｙ、Ｋ）を、そのインク量Ｉｊ化させて所定の階調値が形成されている。なお、各パッチ（ｉ）を形成する１次色インクの階調値は、このインクを搭載するプリンターが使用する階調値に応じて設定される。なお、本実施形態では、インク量Ｉｊを１次色インクにより構成するが、１次色インクの組合せ等とするものであってもよい。また、本明細書において「カラーチャート」という用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。
【００２３】
ステップＳ３では、コンピューター１０は、測色機２１にカラーチャートを測色させて、各パッチ（ｉ）の測色値を取得する。図４は、測色機におけるカラーチャートの測色方法を示す図である。ここで、図４Ａは、測色機２１における測色方法を示す。また、図４Ｂは、測色機２２における測色方法を示す。
【００２４】
図４Ａに示すように、測色機２１は、光源をカラーチャートに対して垂直（０度）に配置し、測色部をカラーチャートに対して４５度に配置し、各パッチ（ｉ）を測色する。そのため、測色機２１により、各パッチ（ｉ）に対して拡散反射での測色値が取得される。そして、測色機２１は、カラーチャートを形成するすべてのパッチ（ｉ）に対して測色値を取得する。
【００２５】
ステップＳ４では、コンピューター１０は、測色機２２にカラーチャートを測色させて各パッチ（ｉ）の測色値を取得する。図４Ｂに示すように、測色機２１は、光源をカラーチャートに対して４５度に配置し、測色部をカラーチャートに対して−４５度に配置し、各パッチ（ｉ）を測色する。そのため、測色機２１により、各パッチ（ｉ）に対して正反射での測色値が取得される。
【００２６】
ステップＳ５では、コンピューター１０は、ステップＳ３及びＳ４で取得した各測色値からブロンジング指標値ＢＩを取得する。
一般にブロンジングとして認識される現象として以下の２点が挙げられる。
（現象１）：観察光源に対して正反射である角度で観察した色と、拡散反射で観察した色との、色味（色相）が異なる。
（現象２）：観察光源に対して正反射である角度で観察した色と、観察光源の色との、色味（色相）が異なる。
そのため、本実施形態では、上記（現象１）及び（現象２）を定量化し、ブロンジング指標値ＢＩを算出する。
【００２７】
図５は、ステップＳ５の処理を詳細に示すフローチャートである。また、図６は、ブロンジング指標値を説明するための図である。
【００２８】
図５のステップＳ５１では、コンピューター１０は、拡散反射における色相角Ｈｕｅｓｔｄと、正反射における色相角Ｈｕｅｓｐｅｃｕｌａｒとの差分となるΔＨｕｅを算出する。即ち、本ステップでは、（現象１）として示した、正反射で観察した色味と、通常の観察（拡散反射）での色味の違いを、色相角Ｈｕｅの差として求める。図６Ａにおいて、色相角は、色相環上の角度として算出される。ここで、正反射における色味と拡散反射における色味が異なる場合、色相環上の色相角Ｈｕｅも異なるため、両色味の差が大きくなるほど、差となる角度も大きくなる。そのため、本ステップでは、（現象１）を色相角Ｈｕｅの差として定量化する。
【００２９】
そのため、コンピューター１０は、ステップＳ３により取得した拡散反射におけるパッチ（ｉ）の測色値（Ｌ＊ａ＊ｂ＊値）から色相角Ｈｕｅｓｔｄを算出する。また、コンピューター１０は、ステップＳ４により取得した正反射におけるパッチ（ｉ）の測色値（Ｌ＊ａ＊ｂ＊値）から色相角Ｈｕｅｓｐｅｃｕｌａｒを算出する。一般に、色相角Ｈｕｅは、下記式（１）により求めることができる。
【数１】
0005682759.tif 000002
ここで、ａ＊はｂ＊は、測色値（Ｌ＊ａ＊ｂ値）におけるａ＊、ｂ＊の値である。
【００３０】
そうして求めた角色相角Ｈｕｅｓｔｄ、Ｈｕｅｓｐｅｃｕｌａｒもとに下記式（２）を用いて、パッチ（ｉ）（インク階調値）毎の色相角の差ΔＨｕｅを算出する。
【数２】
0005682759.tif 000003
ただし、上記ΔＨｕｅの値が取りうる範囲は、０≦ΔＨｕｅ≦２π（Ｒａｄ）であるものとする。上記ステップＳ５１の処理により、本発明の差分算出工程、差分算出手段を実現する。
【００３１】
ステップＳ５２では、コンピューター１０は、正反射における基準彩度からの差ΔＣｈｒｏｍａを算出する。即ち、本ステップでは、（現象２）として示した、観察光源の色味に対して、正反射で異なる色味が見える現象を、正反射における彩度の差ΔＣｈｒｏｍａとして定量化する。図６Ｂに示すように、同一観察光源で正反射での色味が異なる場合、色相環上で観察光源下での彩度（基準彩度とも記載する。）に対して、正反射で対象物を観察した彩度Ｃｈｒｏｍａｓｐｅｃｕｌａｒが異なる。ここで、彩度Ｃｈｒｏｍａは、一般には、下記式（３）により算出することができる。
【数３】
0005682759.tif 000004
【００３２】
本実施形態では、測色機１及び測色機２において観察光源を白色光とするため、基準彩度は限りなく小さい値（≒０）となる。そのため、正反射における彩度を彩度の差ΔＣｈｒｏｍａとして使用することができる。上記ステップＳ５２の処理により、本発明の彩度差算出工程、彩度差算出手段を実現する。
【００３３】
ステップＳ５３では、コンピューター１０は、（現象１）（現象２）を定量化した各値（ΔＨｕｅ、ΔＣｈｒｏｍａ）を乗算してブロンジング指標値ＢＩ（Ｉｊ）を算出する。即ち、ブロンジング指標値ＢＩ（Ｉｊ）は、下記式（４）により算出することができる。
【数４】
0005682759.tif 000005
【００３４】
ここで、ｗｅｉｇｈｔは、低彩度においてブロンジング指標値ＢＩ（Ｉｊ）が低くなるよう調整する重みである。本実施形態では、重みｗｅｉｇｈｔは、下記式（５）に示すようにシグモイド関数として算出される。
【数５】
0005682759.tif 000006
【００３５】
図７は、本実施形態に係る重みｗｅｉｇｈｔを示す図である。重みｗｅｉｇｈｔは、シグモイド関数であるため、拡散反射での彩度（Ｃｈｒｏｍａｓｐｅｃｕｌａｒ）の値に応じて単調増加し、１つの変曲点をもつ。即ち、低彩度域においては、重みｗｅｉｇｈｔが限りなく小さく、ある彩度（Ｃｈｒｏｍａｍａｘ）を境に重みｗｅｉｇｈｔが増加し、やがて重みｗｅｉｇｈｔが「１」付近に飽和する。そのため、この重みｗｅｉｇｈｔをブロンジング指標値ＢＩ（Ｉｊ）の一部として付与することにより、低彩度域ではブロンジング指標値ＢＩ（Ｉｊ）の値が小さくなり（即ち、指標値として良い値となる。）、高彩度域ではブロンジング指標値ＢＩ（Ｉｊ）の値が大きくなる（即ち、指標値として悪い値となる。）。
【００３６】
ブロンジング指標値ＢＩ（Ｉｊ）を彩度に応じたシグモイド関数により重み付けする理由として、色相角の差ΔＨｕｅと彩度の差ΔＣｈｒｏｍａとはその取りうる値の範囲が異なるため、影響度合いを均一にする必要があるからである。即ち、色相角の差ΔＨｕｅは、０≦ΔＨｕｅ≦２π（なお、単位はｒａｄ）の範囲での値であるのに対して、彩度の差ΔＣｈｒｏｍａは、０≦Ｃｈｒｏｍａ≦５０の範囲での値となる。そのため、単純に両者（ΔＨｕｅ、ΔＣｈｒｏｍａ）を乗算すると、値が小さい色相角の差ΔＨｕｅが支配的となってしまう。また、ブロンジング指標値ＢＩ（Ｉｊ）は、彩度が低彩度であるほど、その値が良くなる（指標値が低くなる）よう設定される必要がある。即ち、ΔＣｈｒｏｍａの算出において正反射での彩度Ｃｈｒｏｍａｓｐｅｃｕｌａｒが小さくなるほどブロンジングが発生し難くなるよう値を調整する必要がある。そのため、低彩度ほど、ブロンジング指標値ＢＩ（Ｉｊ）が低い値を取るよう、重みｗｅｉｇｈｔを付与している。上記ステップＳ５３の処理により、本発明の指標値算出工程、指標値算出手段を実現する。
【００３７】
図２に戻り、ステップＳ６では、コンピューター１０は、ステップＳ５で取得した指標値を指標値テーブルＴａに記録する。図８は、指標値テーブルＴａを説明する図である。図８に示すように、指標値テーブルＴａには、プリンターが使用する用紙毎（本実施形態では、「普通紙」、「厚紙」、「コート紙」、「特殊紙」）に、ブロンジング指標値ＢＩ（Ｉｊ）が設定される。即ち、用紙毎に、インク量Ｉｊ（階調値）と、このインク量Ｉｊに対応するブロンジング指標値ＢＩ（Ｉｊ）が記録される。そのため、後述するプリンター２０におけるＬＵＴの作成において、この指標値テーブルＴａが参照され、最適なインク量Ｉｊの探索が行われる。無論、指標値テーブルＴａに記録される指標値は用紙毎の形式に限定されるものではなく、コンピューター１０における印刷制御情報（例えば、印刷時に噴射されるインク滴重量の情報）などを加えてもよい。
【００３８】
２．ブロンジング指標値算出方法２：
ブロンジング指標値ＢＩを算出する過程において、色相角の差ΔＨｕｅと彩度の差ΔＣｈｒｏｍａをそれぞれ正規化してもよい。上記のように、色相角の差ΔＨｕｅは、彩度の差ΔＣｈｒｏｍａに比べて値が小さいため、単純に乗算するとブロンジング指標値ＢＩにおいて、色相角の差ΔＨｕｅが支配的になってしまう。そのため、色相角の差ΔＨｕｅと彩度の差ΔＣｈｒｏｍａをそれぞれ０以上、１以下の値に正規化し、正規化後の値によりブロンジング指標値ＢＩ（Ｉｊ）を算出してもよい。
【００３９】
即ち、図５のステップＳ５３において、コンピューター１０は、正規化後の色相角の差ΔＨｕｅと彩度の差ΔＣｈｒｏｍａを用いてブロンジング指標値ＢＩ（Ｉｊ）を算出する。そして、図２のステップＳ６において、コンピューター１０は、算出されたブロンジング指標値ＢＩ（Ｉｊ）を指標値テーブルＴａに記録する。
【００４０】
以上説明したように、本ブロンジング指標値算出方法により、所定印刷条件下での印刷物におけるブロンジングの程度を指標値として算出するため、印刷物に発生しているブロンジングの程度の評価等を定量的に扱うことができる。
【００４１】
３．ＬＵＴ作成方法：
以下、本指標値作成方法により作成されたブロンジング指標値ＢＩを用いたＬＵＴ作成方法を説明する。図９は、ＬＵＴ作成システム１１０を示すブロック構成図である。本ＬＵＴ作成システム１１０は、コンピューター１０と、プリンター２０と、測色機（不図示）により実現される。
【００４２】
ＬＵＴ作成方法に係るコンピューター１０の記録部１２には、プログラムや色変換プロファイル（例えば、後述するインバースモデル初期ＬＵＴ４１０、色補正なしＬＵＴ５１０、色補正ＬＵＴ６１０）、ブロンジングに係る指標値が記録された指標値テーブルＴａ等が記憶される。なお、「ＬＵＴ」は、ルックアップテーブルの略語である。
【００４３】
色補正なしＬＵＴ５１０は、所定の入力色空間（例えばＲＧＢ色空間）の階調値をプリンター２０で使用される複数種類のインクの組み合わせ（インク量Ｉｊ）に変換するための色変換テーブルである。色補正なしＬＵＴ５１０の入力色空間であるＲＧＢ色空間は、いわゆる機器依存色空間では無く、特定のデバイスとは無関係に設定された仮想の色空間（あるいは抽象的な色空間）である。この色補正なしＬＵＴ５１０は、例えば色補正ＬＵＴ６１０を作成する際に使用される。色補正ＬＵＴ６１０は、標準的な機器依存の色空間（例えばｓＲＧＢ色空間やＪＡＰＡＮ  ＣＯＬＯＲ  ２００１色空間）を、特定のプリンターのインク量Ｉｊに変換するためのルックアップテーブルである。インバースモデル初期ＬＵＴ４１０については後述する。 
【００４４】
また、コンピューター１０のＬＵＴ作成モジュールは、初期値設定部１４と、スムージング処理部１５と、テーブル作成部１６と、ガマット予測部１７と、フォワードモデルコンバーター１８と、を有している。これらの各部の機能については後述する。
【００４５】
＜＜プロファイルの作成について＞＞
図１０は、色補正なしＬＵＴ作成時の全体処理手順を示すフロー図である。
ステップＳ１００では、インクデータの設定が行われる。ここで、インクデータとは、プリンター２０で使用されるインクの色や種類、デューティ制限値などである。
【００４６】
また、デューティ制限値とは、メディア（用紙）の単位面積の中に、打ち込むことのできるインク量Ｉｊを示すものである。このデューティ制限値は、メディアの種類や形成する画像の解像度に依存する。このため、デューティ制限値は、メディアの種類毎および印刷解像度毎に定められる。本実施形態では、写真用紙に１４４０×１４４０ｄｐｉの解像度で印刷を行うこととする。また、デューティ制限値として、インク単色によるデューティ制限値や、複数色の合計によるデューティ制限値（合計デューティ制限値）がある。例えば、ある紙のある解像度におけるデューティ制限値として、単色では８０％、各色の合計では１２０％のように定められる。
【００４７】
ステップＳ２００では、カラーチャートの印刷が行われる。印刷されるカラーチャートは、図３に示すものと同様であり、複数のパッチ（ｉ）が配列して構成されたものである。また、各パッチ（ｉ）は、所定のインク量Ｉｊを用いて形成されている。ここではデューティ制限値や色変換テーブルを用いずに、簡単なインク量Ｉｊの組み合わせによりカラーチャートの印刷を行なう。
【００４８】
ステップＳ３００では、例えば、不図示の測色機やスキャナなどを用いて、カラーチャートを構成する各パッチの測色が行われる。この測色結果は、コンピューター１０に取り込まれ、コンピューター１０においてインク量Ｉｊと測色値（Ｌ＊ａ＊ｂ＊値）とを対応させた仮想サンプルが形成される。
【００４９】
そして、ステップＳ４００では、カラーチャートの測色結果（仮想サンプル）に基づいてフォワードモデル（以下ＦＭともいう）の作成が行われる。なお、「フォワードモデル」とは、インク量Ｉｊを機器独立色空間の色彩値（測色値）に変換する変換モデルを意味する。逆に、「インバースモデル」とは、機器独立色空間の色彩値をインク量Ｉｊに変換する変換モデルを意味する。本実施形態では、機器独立色空間としてＣＩＥ-Ｌ＊ａ＊ｂ＊色空間を使用する。なお、以下では、ＣＩＥ-Ｌ＊ａ＊ｂ＊色空間の色彩値を、単に「Ｌ＊ａ＊ｂ＊値」又は「Ｌａｂ値」とも呼ぶ。
【００５０】
図１１は、ＦＭコンバーターとインバースモデル初期ＬＵＴの説明図である。図１１に示すように、ＦＭコンバーター１８の前段を構成する分光プリンティングコンバーターＲＣは、複数種類のインクのインク量Ｉｊを、そのインク量Ｉｊに応じて印刷されるパッチの分光反射率Ｒ（λ）に変換する。本実施形態では、前述した４色のインクを使用するプリンター２０を想定しており、分光プリンティングモデルコンバーターＲＣもこの４種類のインクのインク量Ｉｊを入力としている。但し、プリンター２０で使用する複数種類のインクとしては、任意のインク量Ｉｊを利用することが可能である。
【００５１】
色算出部ＣＣは、分光反射率Ｒ（λ）からＬ＊ａ＊ｂ＊色空間の色彩値を算出する。この色彩値の算出には、予め選択された光源（例えば標準の光Ｄ５０）がカラーチャートの観察条件として使用される。なお、分光プリンティングモデルコンバーターＲＣを作成する方法としては、例えば特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
【００５２】
インバースモデル初期ＬＵＴ４１０は、Ｌ＊ａ＊ｂ＊値を入力とし、インク量Ｉｊを出力とするルックアップテーブルである。この初期ＬＵＴ４１０は、例えば、Ｌ＊ａ＊ｂ＊色空間を複数の小セルに区分し、小セル毎に最適なインク量Ｉｊを選択して登録したものである。この選択は、例えば、そのインク量Ｉｊで印刷されるカラーチャートの画質を考慮して行われる。一般に、或る１つのＬ＊ａ＊ｂ＊値を再現するインク量Ｉｊの組み合わせは多数存在する。そこで、初期ＬＵＴ４１０では、ほぼ同じＬ＊ａ＊ｂ＊値を再現する多数のインク量Ｉｊの組み合わせの中から、画質等の所望の観点から最適なインク量Ｉｊを選択したものが登録されている。この初期ＬＵＴ４１０の入力値であるＬ＊ａ＊ｂ＊値は各小セルの代表値である。一方、出力値であるインク量Ｉｊはそのセル内のいずれかのＬ＊ａ＊ｂ＊値を再現するものである。従って、この初期ＬＵＴ４１０では、入力値であるＬ＊ａ＊ｂ＊と出力値であるインク量Ｉｊとが厳密に対応したものとなっておらず、出力値のインク量ＩｊをＦＭコンバーター１８でＬ＊ａ＊ｂ＊値に変換すると、初期ＬＵＴ４１０の入力値とは多少異なる値が得られる。但し、初期ＬＵＴ４１０として、入力値と出力値とが完全に対応するものを利用してもよい。また、初期ＬＵＴ４１０を用いずに色補正なしＬＵＴ５１０を作成することも可能である。なお、小セル毎に最適なインク量Ｉｊを選択して初期ＬＵＴ４１０を作成する方法としては、例えば上記特表２００７−５１１１７５号公報に記載された方法を採用することが可能である。
【００５３】
図１０のステップＳ５００では、フォワードモデル（ＦＭ）とインクのデューティ制限値を用いてガマットの予測が行われる。ステップＳ６００では、ステップＳ５００で予測されたガマット（以下、予測ガマットともいう）の評価が行われる。なお、この判断は、例えばユーザーが行う。予測ガマットの評価がＮＧであれば（Ｓ６００でＮＯ）、ステップＳ１００に戻る。
【００５４】
一方、予測ガマットの評価がＯＫであれば（Ｓ６００：ＹＥＳ）、ステップＳ７００で、色補正なしＬＵＴ５１０作成のための初期入力値がユーザーによって設定される。色補正なしＬＵＴ５１０の入力値としては、ＲＧＢの各値として予め定められたほぼ等間隔の値が設定される。１組のＲＧＢ値はＲＧＢ色空間内の点を表していると考えられるので、１組のＲＧＢ値を「入力格子点」とも呼ぶ。
【００５５】
ステップＳ７００においては、複数の入力格子点のうちから予め選択された幾つかの少数の入力格子点に対するインク量Ｉｊの初期値がユーザーによって入力される。この初期入力値が設定される入力格子点としては、ＲＧＢ色空間における３次元色立体の頂点に相当する入力格子点を少なくとも選択することが好ましい。この３次元色立体の頂点では、ＲＧＢの各値がその定義範囲の最小値又は最大値を取る。具体的には、ＲＧＢの各値を８ビットで表現した場合には、(R,G,B)=(0,0,0)、(0,0,255)、(0,255,0)、(255,0,0)、(0,255,255)、(255,0,255)、(255,255,0)、(255,255,255)である８つの入力格子点に関してインク量Ｉｊの初期入力値が設定される。
【００５６】
図１０のステップＳ８００では、スムージング処理部１５（図９）が、ステップＳ７００で設定された初期入力値に基づいてスムージング処理（平滑化処理）を実行する。
図１２は、スムージング処理における入力格子点とＬ＊ａ＊ｂ＊色空間の座標値との対応関係を示す図である。なお、図１２Ｂにはスムージング処理前の状態における複数の色点の分布が２重丸と白丸とで示されている。これらの座標値は、Ｌ＊ａ＊ｂ＊色空間における３次元色立体ＣＳを構成している。なお、各座標値のＬ＊ａ＊ｂ＊値は、色補正なしＬＵＴ５１０の複数の入力格子点におけるインク量Ｉｊを、ＦＭコンバーター１８を用いてＬ＊ａ＊ｂ＊値に変換した値である。なお、他の入力格子点に対するインク量Ｉｊの初期値は、初期入力値から初期値設定部１４（図９参照）によって設定される。
【００５７】
また、図１２Ｃには、スムージング処理後の色点の分布が示されている。スムージング処理は、Ｌ＊ａ＊ｂ＊色空間における複数の座標値を移動させて、それらの座標値の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各座標値のＬ＊ａ＊ｂ＊値を再現するために最適なインク量Ｉｊも決定される。このとき、入力格子点に割り当てられる最適なインク量Ｉｊを目的関数Ｅをもとに算出する。
【００５８】
Ｌ＊ａ＊ｂ＊色空間の３次元色立体ＣＳ（図１２Ｂ、図１２Ｃ）の上記各頂点は、色補正なしＬＵＴ５１０の入力色空間（ＲＧＢ色空間）の３次元色立体（図１２Ａ）の頂点と一対一に対応している。また、各頂点を結ぶ辺（稜線）も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前（図１２Ｂ）のＬ＊ａ＊ｂ＊色空間の各座標値は、色補正なしＬＵＴ５１０の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後（図１２Ｃ）のＬ＊ａ＊ｂ＊色空間の各色点も色補正なしＬＵＴ５１０の入力格子点にそれぞれ対応付けられる。なお、色補正なしＬＵＴ５１０の入力格子点はスムージング処理によって変化しない。
【００５９】
色補正なしＬＵＴ５１０を作成する際に、Ｌ＊ａ＊ｂ＊空間においてスムージング処理を行う理由は以下の通りである。色補正なしＬＵＴ５１０では、なるべく大きな色域を再現できるように出力色空間のインク量Ｉｊを設定したいという要望がある。特定のインクセットで再現可能な色域は、インクデューティ制限値（一定面積に吐出可能なインク量Ｉｊの制限値）などの所定の制限条件を考慮して決定される。一方、上述したフォワードモデルＦＭは、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にインクデューティ制限等の制限条件を考慮してＬ＊ａ＊ｂ＊空間内の座標値の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、座標値の移動を行うアルゴリズムとしては、例えば、力学モデルを使用したものが利用される（力学モデルについては、例えば特開２００６−１９７０８０号公報参照）。
【００６０】
図１３は、スムージング処理（図１０のステップＳ８００）の処理手順の一例を示すフロー図である。ステップＳ８１０では、初期値設定部１４（図９）が、スムージング処理の対象とする複数の座標値を初期設定する。具体的には、まず、インク量Ｉｊの初期入力値から、スムージング処理の対象となる各座標値の仮インク量Ｉｊが決定される。次に、ＦＭコンバーター１８を用いて、仮インク量Ｉｊに対応するＬ＊ａ＊ｂ＊値を求める。そして、得られたＬ＊ａ＊ｂ＊値を、インバースモデル初期ＬＵＴ４１０（図９）を用いてインク量Ｉｊに再度変換する。
【００６１】
上述のステップＳ８１０の処理の結果、スムージング処理の対象となる座標値について、以下の初期値が決定される。
（ｉ）色補正なしＬＵＴの入力格子点の値：(R,G,B)
（ii）各入力格子点に対応するＬ＊ａ＊ｂ＊色空間の座標値の初期座標値：L(R,G,B)
（iii）各入力格子点に対応する初期インク量Ｉｊ：I(R,G,B)
以上の説明から理解できるように、初期値設定部１４は、代表的な入力格子点に関する入力初期値から他の入力格子点に関する初期値を設定する機能を有している。なお、初期値設定部１４は、スムージング処理部１５に含まれるものとしてもよい。
【００６２】
ステップＳ８２０では、スムージング処理部１５は、Ｌ＊ａ＊ｂ＊空間内の座標値を移動させる。
図１４は、図１３のステップＳ８２０〜Ｓ８５０の処理内容を示す説明図である。図１４Ａに示すように、スムージング処理前には、座標値の分布にはかなりの偏りがある。図１４Ｂは、微少時間経過後の各色点の位置を示している。この移動後のＬ＊ａ＊ｂ＊値を「ターゲット値LABt」と呼ぶ。「ターゲット」という修飾語は、この値LABtが、以下で説明するインク量Ｉｊの最適値の探索処理の際の目標値として使用されるからである。
【００６３】
ステップＳ８３０では、スムージング処理部１５は、以下の式（６）で表される目的関数Ｅを用いて、ターゲット値LABtに対するインク量Ｉｊの最適値を探索する（図１４Ｃ参照）。
【数６】
0005682759.tif 000007
ここで、Ｉｊは目的関数Ｅを最小とするインク量Ｉｊ、LABFM(Ij)はインク量ＩｊをＦＭコンバーター１８で変換して得られるＬ＊ａ＊ｂ＊値、ＧＩ（Ｉｊ）はインク量Ｉｊで印刷されるカラーチャートの粒状性指数、ＢＩ（Ｉｊ）はブロンジング指標、α、βは定数である。（６）式の右辺第１項は、２つのベクトルLABt, LABFM(Ij)で表される２つの座標値の距離の２乗を求める演算を意味している。
【００６４】
上記式（６）は、ステップＳ８２０の処理において微少量だけ移動した後の色点の座標値LABtに近いＬ＊ａ＊ｂ＊値を再現するインク量Ｉｊの中で、粒状性指数ＧＩ及びブロンジング指標値ＢＩのより小さいインク量Ｉｊが最適なインク量Ｉｊとして決定されることを意味している。この最適なインク量Ｉｊの探索は、例えば準ニュートン法などの最適化手法を用いて実行される。また、インク量Ｉｊの探索は、ステップＳ８１０で設定された各入力格子点の初期インク値から開始される。従って、探索で得られるインク量Ｉｊは、この初期インク量Ｉｊを修正した値となる。
【００６５】
なお、一般に、最適化処理によるインク量Ｉｊの探索は以下の条件の下で実行される。
（条件ｉ）目的関数Ｅを最小とする
（条件ii）使用可能として予め指定されたインクのみを使用する
【００６６】
条件iの目的関数Ｅを与える式としては、上記（６）式以外の種々の式を利用可能である。また、条件iiにおける使用可能なインクとしては、プリンター２０で利用可能な複数種類のインクのうちの一部のインクのみが指定される場合と、全種類のインクが指定される場合とがある。使用可能なインクの種類は、スムージング処理の対象となる色点毎にユーザーによって予め設定されることが好ましい。これらの最適化条件の具体例については更に後述する。
【００６７】
ここで、ブロンジング指標値ＢＩは、図８に示す指標値テーブルＴａに記録された各値が取得される。このとき、カラーチャートが形成されたメディアの種別に応じて、インク量Ｉｊに対応するブロンジング指標値ＢＩ（Ｉｊ）が選択される。
【００６８】
また、粒状性指数ＧＩ（Ｉｊ）は、例えば以下の式（７）で与えられる。
【数７】
0005682759.tif 000008
ここで、ａＬは明度補正係数、ＷＳ（ｕ）はカラーチャートの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、ＶＴＦ（ｕ）は視覚の空間周波数特性、ｕは空間周波数である。粒状性指数ＧＩの計算方法としては、例えば、特表２００７−５１１１６１に記載されている。
【００６９】
図１３のステップＳ８４０では、ステップＳ８３０で探索されたインク量Ｉｊに対応するＬ＊ａ＊ｂ＊値が、ＦＭコンバーター１８で再計算される（図１４Ｄ）。ここで、Ｌ＊ａ＊ｂ＊値を再計算する理由は、探索されたインク量Ｉｊが目的関数Ｅを最小とするインク量Ｉｊなので、そのインク量Ｉｊで再現されるＬ＊ａ＊ｂ＊値は、最適化処理のターゲット値LABtから多少ずれているからである。こうして再計算されたＬ＊ａ＊ｂ＊値が、各色点の移動後の座標値として使用される。
【００７０】
ステップＳ８５０では、各色点の座標値の移動量の平均値(ΔLab)aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値(ΔLab)aveが閾値εよりも大きい場合には、ステップＳ８２０に戻りステップＳ８２０〜Ｓ８５０のスムージング処理が継続される。一方、移動量の平均値(ΔLab)aveが閾値ε以下の場合には、色点の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。
【００７１】
このように、本実施形態のスムージング処理（平滑化処理）では、各色点を微少時間間隔毎に移動させつつ、移動後の色点に対応する最適なインク量Ｉｊを最適化手法で探索する。そして、色点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図１４Ｃに示したように、平滑な色点分布を得ることが可能である。
【００７２】
図１０に戻り、ステップＳ９００では、スムージング処理の結果を用いて、テーブル作成部１６が色補正なしＬＵＴ５１０を作成する。すなわち、テーブル作成部１６は、各入力格子点に対応付けられたＬ＊ａ＊ｂ＊色空間の色点を再現するための最適なインク量Ｉｊを色補正なしＬＵＴ５１０の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、色補正なしＬＵＴ５１０の入力格子点の一部のみに対応する色点のみを処理対象として選択することも可能である。例えば、色補正なしＬＵＴ５１０の入力格子点におけるＲＧＢ値の間隔が１６である場合に、スムージング処理の対象となる入力格子点におけるＲＧＢ値の間隔を３２に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成部１６は、スムージング処理結果を補間することによって色補正なしＬＵＴ５１０のすべての入力格子点に対するインク量Ｉｊを決定して登録する。
【００７３】
そして、図１０のステップＳ１０００では、色補正なしＬＵＴ５１０の評価が行われる。評価がＮＧであれば（Ｓ１０００でＮＯ）、再度ステップＳ１００に戻る。一方、評価がＯＫであれば、色補正なしＬＵＴ５１０を記憶部１２に登録する。
なお、ＬＵＴ作成モジュールは、色補正なしＬＵＴ５１０を用いて、色補正ＬＵＴ６１０の作成も行う。
【００７４】
図１５は、色補正なしＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。上記のように色補正ＬＵＴ６１０は、標準的な機器依存の色空間を、特定のプリンターのインク量Ｉｊに変換するためのルックアップテーブルであり、本実施形態では、機器依存の色空間としてｓＲＧＢ色空間を採用する。以下、図１５に示すように、色補正なしＬＵＴが規定するＲＧＢ値とインク量Ｉｊの関係をもとに、色補正ＬＵＴ６１０の作成が行われる。
【００７５】
まず、色補正なしＬＵＴ５１０により、ＲＧＢ値がインク量Ｉｊに変換される。変換後のインク量Ｉｊは、ＦＭコンバーター１８によってＬ＊ａ＊ｂ＊値に変換される。一方、色補正ＬＵＴ６１０に対応付けられるｓＲＧＢ値は、既知の変換式に従ってＬ＊ａ＊ｂ＊値に変換される。そして、変換後のＬ＊ａ＊ｂ＊値は、その色域が、ＦＭコンバーター１８で変換されたＬ＊ａ＊ｂ＊値の色域と一致するようにガマットマッピングされる。
【００７６】
一方、色補正なしＬＵＴ５１０とＦＭコンバーター１８を通じて、ＲＧＢ値から変換したＬ＊ａ＊ｂ＊値を、逆方向ルックアップテーブルとして、逆変換ＬＵＴ５１１が作成される。即ち、逆変換ＬＵＴ５１１は、Ｌ＊ａ＊ｂ＊値とＲＧＢ値とを対応付けるＬＵＴである。さらに、ガマットマッピングされたＬ＊ａ＊ｂ＊値は、この逆変換ＬＵＴ５１１によってＲＧＢ値に変換される。このＲＧＢ値は、さらに、色補正なしＬＵＴ５１０によってインク量Ｉｊに再度変換される。そして、最後のインク量Ｉｊと最初のｓＲＧＢ値の対応関係をルックアップテーブルに登録することによって、色補正ＬＵＴ６１０が作成される。
【００７７】
以上により、機器依存の色空間とインク量Ｉｊとを対応付ける色補正ＬＵＴ６１０が作成され、記録部１２に記録される。以後、コンピューター１０では、この色補正ＬＵＴ６１０を使用して入力データに対応するインク量Ｉｊを選択してプリンター２０に出力する。また、本実施形態では、色補正ＬＵＴ６１０の元となる色補正なしＬＵＴ５１０の作成において、最適なインク量Ｉｊの評価項目としてブロンジング指標値ＢＩ（Ｉｊ）が用いられる。そのため、色補正ＬＵＴ６１０に記録されるインク量Ｉｊは、ブロンジングを抑制するようその組合せが選択される。
【００７８】
４．その他の実施形態：
記録剤としてインクを用いることは一例である。例えば、プリンター２０としてレザープリンターを用いる場合は、トナーのような記録剤に対して本指標値算出方法を適用するものであってもよい。
また、ブロンジング指標値の適用例としてＬＵＴ作成方法を用いたことは一例であり、これ以外にも様々な適用例がある。
【００７９】
なお、本発明は上記実施例に限られるものでないことは言うまでもない。即ち、上記実施例の中で開示した相互に置換可能な部材および構成等を適宜その組み合わせを変更して適用すること、上記実施例の中で開示されていないが、公知技術であって上記実施例の中で開示した部材および構成等と相互に置換可能な部材および構成等を適宜置換し、またその組み合わせを変更して適用すること、上記実施例の中で開示されていないが、公知技術等に基づいて当業者が上記実施例の中で開示した部材および構成等の代用として想定し得る部材および構成等と適宜置換し、またその組み合わせを変更して適用すること、は本発明の一実施例として開示されるものである。
【符号の説明】
【００８０】
１０…コンピューター、１１…制御部、１２…記憶部、２０…プリンター、２１,２２…測色機、２３…ＵＳＢＩ／Ｆ、３０…ディスプレイ、３１…ビデオＩ／Ｆ、４０…操作部、４１…入力Ｉ／Ｆ、１００…指標値取得システム、１１０…ＬＵＴ作成システム 
【図１】
0005682759.tif 000009
【図２】
0005682759.tif 000010
【図３】
0005682759.tif 000011
【図４】
0005682759.tif 000012
【図５】
0005682759.tif 000013
【図６】
0005682759.tif 000014
【図７】
0005682759.tif 000015
【図８】
0005682759.tif 000016
【図９】
0005682759.tif 000017
【図１０】
0005682759.tif 000018
【図１１】
0005682759.tif 000019
【図１２】
0005682759.tif 000020
【図１３】
0005682759.tif 000021
【図１４】
0005682759.tif 000022
【図１５】
0005682759.tif 000023