(57)【要約】      （修正有）【課題】  各成分の物性値を用い、無色透明でない被着物質の吸収係数や散乱係数を考慮してコンピュータカラーマッチングにより混合物の色や調合割合を予測する。【解決手段】  被着色物と白色顔料の混合物に対する散乱係数ＳW'を、白色顔料の調合率ＣW の関数で表現する。そして、検証用サンプルに対してコンピュータカラーマッチングを実行し、散乱係数ＳW'を表わす関数に含まれる係数ａ，ｂ，ｃ，ｄ（またはＳB ）として、予測精度が高くなる値を決定する。そして、決定された関数で与えられる散乱係数ＳW'を用いて、（被着色物＋白色顔料）の吸収係数と、他の有色顔料の吸収係数および散乱係数とを求める。こうして求められた各成分の物性値を用いて、コンピュータカラーマッチングを実行し、混合物の色や調合割合を予測する。このコンピュータカラーマッチングにおいては，ニューラルネットワークを利用することによって、予測誤差を減少させることができる。【選択図】  図１
2006030210.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう方法であって、
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'を、前記白色着色剤の調合率ＣW に依存した関数ｆ（ＣW ）として準備する工程と、
前記散乱係数ＳW'を基準として、前記第１の混合物の吸収係数ＫW'を前記調合率ＣW に依存する形式で求める工程と、
前記散乱係数ＳW'を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP とを前記有色着色剤の調合率ＣP に依存する形式で求める工程と、
所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって求める工程と、
を備えるコンピュータカラーマッチング方法。
【請求項２】
請求項１記載の コンピュータカラーマッチング方法であって、
前記第１の混合物の散乱係数ＳW'を準備する工程は、
（ａ）前記被着色物に前記白色着色剤を混合して、前記白色着色剤の調合率ＣW が異なる複数個の第１の混合物を作製するとともに、前記複数個の第１の混合物の分光反射率をそれぞれ測定する工程と、
（ｂ）前記被着色物に前記有色着色物を混合した第２の混合物を作製するとともに、前記第２の混合物の分光反射率を測定する工程と、
（ｃ）前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記有色着色剤の調合率ＣP が異なる複数個の第３の混合物を作製するとともに、前記複数個の第３の混合物の分光反射率をそれぞれ測定する工程と、
（ｄ）前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記第３の混合物とは調合割合が異なる第４の混合物を作製するとともに、前記第４の混合物の分光反射率を測定する工程と、
（ｅ）前記複数個の第１の混合物の散乱係数ＳW'を前記白色着色剤の調合率ＣW の関数ｆ（ＣW ）によって表わすとともに、前記関数ｆ（ＣW ）に含まれる係数の値を仮決定する工程と、
（ｆ）前記複数個の第１の混合物の分光反射率の測定値と前記関数ｆ（ＣW ）とを用いて、前記調合率ＣW に依存した形式で前記第１の混合物の吸収係数ＫW'を求める工程と、
（ｇ）前記第２の混合物の分光反射率と、前記複数個の第３の混合物の分光反射率と、前記関数ｆ（ＣW ）と、前記吸収係数ＫW'とを用いて、前記調合率ＣP に依存した形式で前記有色着色剤の吸収係数ＫP と散乱係数ＳP とをそれぞれ求める工程と、
（ｈ）前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて、前記第４の混合物に関するコンピュータカラーマッチングを行なうとともに、前記コンピュータカラーマッチングによる予測精度を向上させるように前記関数ｆ（ＣW ）に含まれる係数を修正する工程と、
（ｉ）前記工程（ｆ）ないし（ｈ）を繰り返すことによって前記関数ｆ（ＣW ）に含まれる係数を決定する工程と、
を備えるコンピュータカラーマッチング方法。
【請求項３】
前記関数ｆ（ＣW ）が定数ａ1 と係数ＳB とを含む次の式で与えられる、請求項１または２記載のコンピュータカラーマッチング方法：
f(CW )=(CW +a1 ・SB )/(CW +a1 )。
【請求項４】
請求項１記載のコンピュータカラーマッチング方法であって、前記コンピュータカラーマッチングを行なう工程は，
（ａ）複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備する工程と、
（ｂ）前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、
（ｃ）前記複数のサンプルのそれぞれに関して、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて前記表色系の座標値の予測値を求めるとともに，前記表色系の座標値の予測誤差を算出する工程と、
（ｄ）前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、
（ｅ）前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測を，前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて求める工程と、
を備えるコンピュータカラーマッチング方法。
【請求項５】
請求項１記載のコンピュータカラーマッチング方法であって、前記コンピュータカラーマッチングを行なう工程は，
（ａ）着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める工程と、
（ｂ）前記近接色サンプルの既知の調合割合から、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める工程と、
（ｃ）目標混合物の色に対する前記表色系の座標の目標値を設定する工程と、
（ｄ）前記計算誤差を用いて前記目標値を補正し、補正後の目標値と前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP とを用いたコンピュータカラーマッチングによって、前記目標混合物の着色剤の調合割合を予測する工程と、
を備えるコンピュータカラーマッチング方法。
【請求項６】
請求項１記載のコンピュータカラーマッチング方法であって，前記コンピュータカラーマッチングを行なう工程は，
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求める工程を含み，
前記調合割合を求める工程は，
（ａ）前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める工程と、
（ｂ）既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、
（ｃ）前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を，前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて求める工程と、
（ｄ）前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を，前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて求めるとともに、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、
（ｅ）前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する工程と、
を備えるコンピュータカラーマッチング方法。
【請求項７】
コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう方法であって、
（ａ）複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備する工程と、
（ｂ）前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、
（ｃ）前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する工程と、
（ｄ）前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、
（ｅ）前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう工程と、
を備えるコンピュータカラーマッチング方法。
【請求項８】
請求項７記載のコンピュータカラーマッチング方法であって、
前記工程（ｄ）は、前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係をニューラルネットワークに学習させる工程、を含み、
前記工程（ｅ）は、学習済みのニューラルネットワークを用いてコンピュータカラーマッチングによる予測を行なう工程を含む、
コンピュータカラーマッチング方法。
【請求項９】
請求項８記載のコンピュータカラーマッチング方法であって、
前記ニューラルネットワークは、３つのニューロンで構成される入力層と、複数のニューロンを含む中間層と、３つのニューロンで構成される出力層と、で構成される三層の階層構造を有する、コンピュータカラーマッチング方法。
【請求項１０】
所望の色を有する目標混合物の着色剤の調合割合をコンピュータカラーマッチングによって予測する方法であって、
（ａ）着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める工程と、
（ｂ）前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める工程と、
（ｃ）前記目標混合物の色に対する前記表色系の座標の目標値を設定する工程と、
（ｄ）前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する工程と、
を備えるコンピュータカラーマッチング方法。
【請求項１１】
請求項１０記載のコンピュータカラーマッチング方法であって、
前記工程（ａ）は、複数のサンプルに関して、着色剤の調合割合と、前記表色系の座標の実測値とを含むデータベースから、前記目標混合物との色差が最小となるサンプルを選択することによって前記近接色サンプルを検索する工程を含む、
コンピュータカラーマッチング方法。
【請求項１２】
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラーマッチング方法であって、
（ａ）前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める工程と、
（ｂ）既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、
（ｃ）前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める工程と、（ｄ）前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、
（ｅ）前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する工程と、を備える
コンピュータカラーマッチング方法。
【請求項１３】
請求項１２記載のコンピュータカラーマッチング方法であって、
前記工程（ｄ）は、前記１次調合物における前記着色剤の調合量に比べて微量の量の前記着色剤を前記１次調合物に増量補正したと仮定した場合について、前記色評価値の計算値の変化量を求める工程を含む。
【請求項１４】
請求項１２又は請求項１３記載のコンピュータカラーマッチング方法であって、
前記工程（ｅ）は、前記着色剤の増量に伴う派生費用を表わす費用関数を用いた線形計画法にて、前記着色剤のそれぞれについての最小の増量補正量を算出する工程を含む。
【請求項１５】
コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置であって、
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'を、前記白色着色剤の調合率ＣW に依存した関数ｆ（ＣW ）として作製する手段と、
前記散乱係数ＳW'を基準として、前記第１の混合物の吸収係数ＫW'を前記調合率ＣW に依存する形式で求める手段と、
前記散乱係数ＳW'を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP とを前記有色着色剤の調合率ＣP に依存する形式で求める手段と、
所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって求める手段と、
を備えるコンピュータカラーマッチング装置。
【請求項１６】
コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置であって、
複数の着色剤を混合することによって作製された調合率が互いに異なる複数のサンプルの分光反射率をそれぞれ測定する手段と，
前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める手段と、
前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する手段と、
前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する手段と、
前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう手段と、
を備えるコンピュータカラーマッチング装置。
【請求項１７】
所望の色を有する目標混合物の着色剤の調合割合をコンピュータカラーマッチングによって予測する装置であって、
着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める手段と、
前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める手段と、
前記目標混合物の色に対する前記表色系の座標の目標値を設定する手段と、
前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する手段と、
を備えるコンピュータカラーマッチング装置。
【請求項１８】
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラーマッチング装置であって、
前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める手段と、
既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める手段と、
前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める手段と、
前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める手段と、
前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する手段と、を備える
コンピュータカラーマッチング装置。
【請求項１９】
コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法であって、
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する工程と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、
を備える方法。
【請求項２０】
前記第１の混合物の散乱係数ＳW'が，前記白色着色剤の調合率ＣＷと定数ａ1 と係数ＳB とを含む次の関係で与えられる、請求項１９記載の方法：
SW'=(CW +a1 ・SB )/(CW +a1 )。
【請求項２１】
前記第１の混合物の散乱係数ＳW'が，前記白色着色剤の調合率ＣWと定数ａ，ｂ，ｄ，ｅ,ＣW0を含む次の関係で与えられる、請求項１９記載の方法：
SW'= a・CW + b          (CW0 ≦ CW);
SW'= d・CW + e          (CW < CW0)。
【請求項２２】
コンピュータカラーマッチングを利用して作製された混合釉であって、
無色透明でないベース釉に白色着色剤を混合した第１の混合釉の散乱係数ＳW'と吸収係数ＫW'を決定し、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定し、
前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、
前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，
ことによって生成された混合釉。
【請求項２３】
コンピュータカラーマッチングを利用して作製された混合釉を用いて製造された陶磁器であって、
無色透明でないベース釉に白色着色剤を混合した第１の混合釉の散乱係数ＳW'と吸収係数ＫW'を決定し、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定し、
前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、
前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，ことによって生成された混合釉で陶磁器の素地の少なくとも一部が覆われた，陶磁器。
【請求項２４】
コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法であって、
無色透明でない被着色物に基準となる顔料を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する工程と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、
を備える方法。
【請求項２５】
コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する装置であって、
無色透明でない被着色物に基準となる顔料を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する手段と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、
を備える方法。
【請求項２６】
コンピュータカラーマッチングを利用して作製された混合釉を用いて製造された便器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう方法およびそのための装置に関するものである。
【背景技術】
【０００２】
顔料や染料などの着色剤を被着色物に混合した混合物の色を予測するために、いわゆるコンピュータカラーマッチングが利用されている。コンピュータカラーマッチングでは、被着色物と着色剤の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）とを用い、ダンカン（Duncan）の式（数式１）と、クベルカ−ムンク（Kubelka-Munk）の混色理論による式（数式２）に基づいて、任意の混合物の分光反射率Ｒ（λ）を求めることができる。
【０００３】
【数１】
2006030210.tif 000003
【数２】
2006030210.tif 000004
【０００４】
ここで、ＫM ，ＳM は混合物の吸収係数と散乱係数、Ｋi ，Ｓi はｉ番目の成分の吸収係数と散乱係数、Ｃi はｉ番目の成分の調合率である。但し、この明細書の数式においては、波長λに依存していることを示す「（λ）」は省略されている。混合物の成分は、被着色物と着色剤である。
【０００５】
混合物の分光反射率Ｒ（λ）が解れば、その混合物の三刺激値Ｘ，Ｙ，Ｚが計算できるので、混合物の色を算出することができる。
【０００６】
コンピュータカラーマッチングを行なう場合には、被着色物と種々の着色剤の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）を予め求めておく必要がある。ところが、吸収係数Ｋi（λ）と散乱係数Ｓi（λ）を求めるのは必ずしも容易ではない場合が多い。
【０００７】
被着色物と種々の着色剤の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）を求める方法としては、絶対法と相対法がある。絶対法は、各物質の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）の絶対値を求める方法である。一方、相対法は、基準となる顔料（通常は白色顔料）の散乱係数ＳW を１と仮定して、各物質の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）の相対値を求める方法である。数式２に示すように、混合物の分光反射率Ｒ（λ）は混合物の吸収係数ＫM と散乱係数ＳM の比で与えられる。従って、各物質の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）の絶対値が不明でも、各物質の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）の相対値が解れば、数式１および２から正しい分光反射率Ｒ（λ）を求めることができる。絶対法に用いるサンプルの作製はかなり困難であり、煩雑な作業を要するので、相対法が用いられるのが普通である。
【０００８】
従来の相対法によって各物質の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）とを求める際には、以下のようにして、基準となる白色顔料の吸収係数ＫW（λ）と散乱係数ＳW を決定する。まず、被着色物に白色顔料のみを混合した混合物を作製して、その分光反射率Ｒ（λ）を測定する。従来の相対法では、被着色物を無色透明とみなして、被着色物の吸収係数と散乱係数はともに０であると仮定している。従って、被着色物と白色顔料のみを混合した混合物に関しては、上記の数式１で与えられる吸収係数ＫM（λ）と散乱係数ＳM（λ）は、白色顔料の散乱係数ＳW と吸収係数ＫW（λ）にそれぞれ等しい。また、上記の数式２においてＳM ＝ＳW ＝１と仮定すると、分光反射率Ｒ（λ）の測定値から白色顔料の吸収係数ＫW（λ）を求めることができる。
【０００９】
上述のように、従来の相対法では、被着色物を無色透明とみなして、その吸収係数と散乱係数を０と仮定していた。
【００１０】
ところが、現実には陶磁器等をはじめ，被着色物（陶磁器の場合は釉薬層）が無色透明でない場合が多い。白色顔料の調合率が多い場合には、被着色物が無色透明であるとみなすことによる誤差は小さいが、白色顔料の調合率が少ない場合にはその誤差が無視できない程度に大きくなる。このような誤差を回避するためには、上記数式１において、無色透明でない被着色物の吸収係数や散乱係数を考慮にいれる必要がある。しかし、被着色物の吸収係数や散乱係数を求めるには、被着色物単体の薄板を作製するなどの煩雑な作業を要する。特に、陶器の釉のように、単体で薄板を作製することが困難な被着色の場合には、無色透明でない被着色物の吸収係数と散乱係数を求めることは困難であった。
【発明の開示】
【発明が解決しようとする課題】
【００１１】
この発明の第１の目的は、無色透明でない被着色物の吸収係数や散乱係数を考慮してコンピュータカラーマッチングを行なうことにある。
この発明の第２の目的は、コンピュータカラーマッチングにおける予測誤差を減少させることにある。
この発明の第３の目的は、コンピュータカラーマッチングを利用して，調合済み着色剤の有効利用を図りつつ、その再調合を簡略化することにある。
【課題を解決するための手段】
【００１２】
この発明は，コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう方法に関する。このコンピュータカラーマッチング方法は，
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'を、前記白色着色剤の調合率ＣW に依存した関数ｆ（ＣW ）として準備する工程と、
前記散乱係数ＳW'を基準として、前記第１の混合物の吸収係数ＫW'を前記調合率ＣW に依存する形式で求める工程と、
前記散乱係数ＳW'を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP とを前記有色着色剤の調合率ＣP に依存する形式で求める工程と、
所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって求める工程と、
を備える。
【００１３】
上記方法では、無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'を関数ｆ（ＣW ）で表現し、この散乱係数ＳW'を基準として用いるので、散乱係数ＳW'に無色透明でない被着色物の影響を含めることができる。
【００１４】
好ましい実施例においては、
前記第１の混合物の散乱係数ＳW'を準備する工程は、
（ａ）前記被着色物に前記白色着色剤を混合して、前記白色着色剤の調合率ＣW が異なる複数個の第１の混合物を作製するとともに、前記複数個の第１の混合物の分光反射率をそれぞれ測定する工程と、
（ｂ）前記被着色物に前記有色着色物を混合した第２の混合物を作製するとともに、前記第２の混合物の分光反射率を測定する工程と、
（ｃ）前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記有色着色剤の調合率ＣP が異なる複数個の第３の混合物を作製するとともに、前記複数個の第３の混合物の分光反射率をそれぞれ測定する工程と、
（ｄ）前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記第３の混合物とは調合割合が異なる第４の混合物を作製するとともに、前記第４の混合物の分光反射率を測定する工程と、
（ｅ）前記複数個の第１の混合物の散乱係数ＳW'を前記白色着色剤の調合率ＣW の関数ｆ（ＣW ）によって表わすとともに、前記関数ｆ（ＣW ）に含まれる係数の値を仮決定する工程と、
（ｆ）前記複数個の第１の混合物の分光反射率の測定値と前記関数ｆ（ＣW ）とを用いて、前記調合率ＣW に依存した形式で前記第１の混合物の吸収係数ＫW'を求める工程と、
（ｇ）前記第２の混合物の分光反射率と、前記複数個の第３の混合物の分光反射率と、前記関数ｆ（ＣW ）と、前記吸収係数ＫW'とを用いて、前記調合率ＣP に依存した形式で前記有色着色剤の吸収係数ＫP と散乱係数ＳP とをそれぞれ求める工程と、
（ｈ）前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いて、前記第４の混合物に関するコンピュータカラーマッチングを行なうとともに、前記コンピュータカラーマッチングによる予測精度を向上させるように前記関数ｆ（ＣW ）に含まれる係数を修正する工程と、
（ｉ）前記工程（ｆ）ないし（ｈ）を繰り返すことによって前記関数ｆ（ＣW ）に含まれる係数を決定する工程と、
を備える。
【００１５】
こうして関数ｆ（ＣW ）の係数を決定するようにすれば、コンピュータカラーマッチングによる予測精度が高くなるように関数ｆ（ＣW ）の係数を求めることができる。
【００１６】
前記関数ｆ（ＣW ）が定数ａ1 と係数ＳB とを含む次の式で与えられるようにすることが好ましい。
ｆ（ＣW ）＝（ＣW ＋ａ1 ・ＳB ）／（ＣW ＋ａ1 ）
【００１７】
上記の関数ｆ（ＣW ）の形式は、白色着色剤単独の散乱係数を１と仮定した場合に理論的に得られるものなので、散乱係数ＳW'を表わす関数として適切なものである。また、１つの係数ＳB を含むだけの簡単な形式を有しているので、係数ＳB を決定するのが容易である。
【００１８】
この発明の他の局面によれば，コンピュータカラーマッチング方法は、（ａ）複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備する工程と、（ｂ）前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、（ｃ）前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する工程と、（ｄ）前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、（ｅ）前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう工程と、を備える。
【００１９】
複数のサンプルに関する所定の表色系の座標値とその予測誤差とを所定の誤差補正法で分析し、その誤差補正法によってコンピュータカラーマッチングの目標値または予測値を補正しつつ予測を行なうので、各成分の吸収係数Ｋi と散乱係数Ｓi とを補正することなく予測誤差を減少させることができる。
【００２０】
好ましい実施例においては、前記工程（ｄ）は、前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係をニューラルネットワークに学習させる工程、を含み、前記工程（ｅ）は、学習済みのニューラルネットワークを用いてコンピュータカラーマッチングによる予測を行なう工程を含む。
【００２１】
ニューラルネットワークを用いて予測誤差を補正するようにすれば、数多くのサンプルを学習させることによって、予測誤差を減少させることができる。
【００２２】
また、前記ニューラルネットワークは、３つのニューロンで構成される入力層と、複数のニューロンを含む中間層と、３つのニューロンで構成される出力層と、で構成される三層の階層構造を有することが好ましい。
【００２３】
この発明のさらに他の局面によれば，コンピュータカラーマッチング方法は、所望の色を有する目標混合物の着色剤の調合割合を予測する方法であって、
（ａ）着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める工程と、
（ｂ）前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める工程と、
（ｃ）前記目標混合物の色に対する前記表色系の座標の目標値を設定する工程と、
（ｄ）前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する工程と、
を備える。
【００２４】
近接色サンプルの表色系の座標値の計算誤差によって目標混合物の表色系の座標の目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを行なうので、表色系の座標値の計算誤差を低減した条件下でコンピュータカラーマッチングを実行できる。
【００２５】
好ましい実施例によれば、前記工程（ａ）は、複数のサンプルに関して、着色剤の調合割合と、前記表色系の座標の実測値とを含むデータベースから、前記目標混合物との色差が最小となるサンプルを選択することによって前記近接色サンプルを検索する工程を含む。
【００２６】
データベースから色差が最小のサンプルを近接色サンプルとして検索するようにすれば、任意の色を有する目標混合物に関して、コンピュータカラーマッチングの予測誤差を低減できる。
【００２７】
この発明のさらに他の局面によれば，コンピュータカラーマッチング方法は、
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求める方法であって、
（ａ）前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める工程と、
（ｂ）既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、
（ｃ）前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める工程と、
（ｄ）前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、
（ｅ）前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する工程と、を備える。
【００２８】
上記のコンピュータカラーマッチング方法では、既知の調合割合の１次調合物に着色剤を仮に増量補正した場合、その際の色評価値の計算値の変化量を求める。そして、この色評価値の計算値の変化量に基づいてそれぞれの着色剤についての増量補正量を算出し、調合物見本と１次調合物との色評価値の実測値の差を所定範囲で一致させる。このため、求めた増量補正量だけそれぞれの着色剤を１次調合物に実際に追加調合すれば、調合物見本の目標色若しくはこれに近似した色を呈する調合物に、１次調合物を再調合できる。従って、本発明のコンピュータカラーマッチング方法によれば、着色剤の除去を意味する負の補正量を求めることがないので、着色剤が調合済みの調合物（１次調合物）の廃棄が不要となり、既存の調合物を有効に利用することができる。また、この際に１次調合物を既存の調合物ではなく新たに試験的に調合するにしても、その調合に技術者が一回限り関与すればよく、しかもその際に、技術者の長年の勘や経験を要しないので、調合物の再調合を簡略化することができる。
【００２９】
好ましい実施例においては、
前記工程（ｄ）は、前記１次調合物における前記着色剤の調合量に比べて微量の量の前記着色剤を前記１次調合物に増量補正したと仮定した場合について、前記色評価値の計算値の変化量を求める工程を含む。
【００３０】
このコンピュータカラーマッチング方法では、色評価値の計算値の変化量を求めるために行なう着色剤の仮の増量補正を微量単位で行なって、微小単位の変化量を求めることができる。よって、それぞれの着色剤の増量補正量を高い精度で求めることを通して、目標色によく一致した色を呈することができる調合物を得ることができる。
【００３１】
また、前記工程（ｅ）は、前記着色剤の増量に伴う派生費用を表わす費用関数を用いた線形計画法にて、前記着色剤のそれぞれについての最小の増量補正量を算出する工程を含むことが好ましい。
【００３２】
このコンピュータカラーマッチング方法では、それぞれの着色剤を実際に追加増量する際の最小の増量補正量を算出できるので、着色剤の使用量の低減を通して派生費用の最小化を図ることができ、コストを低減することができる。
【００３３】
この発明は，また，コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置に向けられている。この装置は、
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'を、前記白色着色剤の調合率ＣW に依存した関数ｆ（ＣW ）として作製する手段と、
前記散乱係数ＳW'を基準として、前記第１の混合物の吸収係数ＫW'を前記調合率ＣW に依存する形式で求める手段と、
前記散乱係数ＳW'を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP とを前記有色着色剤の調合率ＣP に依存する形式で求める手段と、
所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって求める手段と、
を備える。
【００３４】
この発明の他の局面によれば，コンピュータカラーマッチング装置は、
複数の着色剤を混合することによって作製された調合率が互いに異なる複数のサンプルの分光反射率をそれぞれ測定する手段と，
前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める手段と、
前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する手段と、
前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する手段と、
前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう手段と、
を備える。
【００３５】
この発明のさらに他の局面によれば，コンピュータカラーマッチング装置は、
着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める手段と、
前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める手段と、
前記目標混合物の色に対する前記表色系の座標の目標値を設定する手段と、
前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する手段と、
を備える。
【００３６】
この発明の他の局面によれば，コンピュータカラーマッチング装置は，複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるものである。このコンピュータカラーマッチング装置は、
前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める手段と、
既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める手段と、
前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める手段と、
前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める手段と、
前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する手段と、を備える。
【００３７】
この発明は，さらに，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法に向けられている。この方法は、
無色透明でない被着色物に白色着色剤を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する工程と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、を備える。
【００３８】
好ましい実施例によれば，前記第１の混合物の散乱係数ＳW'が，前記白色着色剤の調合率ＣＷと定数ａ1 と係数ＳB とを含む次の関係で与えられる。
SW'=(CW +a1・SB )/(CW +a1 )
【００３９】
他の好ましい実施例によれば，前記第１の混合物の散乱係数ＳW'が，前記白色着色剤の調合率ＣWと定数ａ，ｂ，ｄ，ｅ，ＣW0を含む次の関係で与えられる。
SW'= a・CW + b          (CW0 ≦ CW)
SW'= d・CW + e          (CW < CW0)
【００４０】
この発明は，また，コンピュータカラーマッチングを利用して作製された混合釉に向けられている。この混合釉は，
無色透明でないベース釉に白色着色剤を混合した第１の混合釉の散乱係数ＳW'と吸収係数ＫW'を決定し、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定し、
前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、
前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，
ことによって生成される。
【００４１】
この発明は，さらに，コンピュータカラーマッチングを利用して作製された混合釉を用いて製造された陶磁器にも向けられている。この陶磁器は、
無色透明でないベース釉に白色着色剤を混合した第１の混合釉の散乱係数ＳW'と吸収係数ＫW'を決定し、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、白色でない有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定し、
前記吸収係数ＫW'，ＫP および前記散乱係数ＳW'，ＳP を用いたコンピュータカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、
前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，ことによって生成された混合釉で陶磁器の素地の少なくとも一部が覆われる。
【００４２】
この発明の他の局面においては，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法は、
無色透明でない被着色物に基準となる顔料を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する工程と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、
を備える。
【００４３】
この発明のさらに他の局面によれば，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する装置は、
無色透明でない被着色物に基準となる顔料を混合した第１の混合物の散乱係数ＳW'と吸収係数ＫW'を決定する手段と、
前記散乱係数ＳW'と吸収係数ＫW'の少なくとも一方を基準として、有色着色剤の吸収係数ＫP と散乱係数ＳP の相対値をそれぞれ決定する工程と、
を備える。
【００４４】
また，この発明は，コンピュータカラーマッチングを利用して作製された混合釉を用いて製造された便器にも向けられている。
【発明を実施するための最良の形態】
【００４５】
Ａ．第１実施例：
図１は、第１実施例における処理の全体手順を示すフローチャートである。なお、この第１実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉（ゆう）である。すなわち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラーマッチングの対象となる混合物である。
【００４６】
なお，陶磁器において，白色顔料としては乳濁剤が有効である。乳濁剤としては，ジルコン等のジルコニウム化合物や，リン酸カルシウム等のリン化合物が用いられる。
【００４７】
ステップＳ１では、被着色物（ベース釉）と白色顔料の混合物の吸収係数ＫW'と散乱係数ＳW'とを求める。ここで、従来の相対法と異なるのは、被着色物を無色透明と仮定していない点、および、被着色物や白色顔料単独の吸収係数や散乱係数を求めずに、その混合物の吸収係数ＫW'と散乱係数ＳW'を求める点にある。なお、散乱係数ＳW'は、白色顔料の調合率ＣW に依存した関数ｆ（ＣW ）として求める。また、吸収係数ＫW'も、白色顔料の調合率ＣW に依存した形式で求める。
【００４８】
なお、ステップＳ１において、白色顔料以外の顔料（以下、「有色顔料」と呼ぶ）の吸収係数ＫP と散乱係数ＳP も決定される（ステップＳ２）。なお、有色顔料の吸収係数ＫP と散乱係数ＳP は、その調合率ＣP に依存した形式で求められる。
【００４９】
ステップＳ３では、ステップＳ１，Ｓ２で求めた吸収係数と散乱係数に基づいて、上述した数式１および数式２を用いてコンピュータカラーマッチングを実行し、混合物の色予測や調合割合の予測を行なう。
【００５０】
ここで、ステップＳ１の詳細手順を説明する前に、ベース釉と白色顔料の混合物の散乱係数ＳW'を表現するための関数ｆ（ＣW ）について考察する。ベース釉と白色顔料の混合物については、数式１を次の数式３のように書き換えることができる。
【００５１】
【数３】
2006030210.tif 000005
【００５２】
ここで、ＫW'，ＳW'は混合物の吸収係数と散乱係数、ＫW ，ＳW は白色顔料単独の吸収係数と散乱係数、ＫB ，ＳB はベース釉単独の吸収係数と散乱係数、ＣW は白色顔料の調合率、ＣB はベース釉の調合率である。なお、調合率ＣW，ＣB は正確には体積率で表現されるべきであるが、重量率で表現してもその誤差は通常無視できる程度である。
【００５３】
数式３を変形すると、混合物の散乱係数ＳW'は次の数式４で与えられる。
【００５４】
【数４】
2006030210.tif 000006
【００５５】
この発明では相対法を用いるので、白色顔料の調合率ＣW によらず散乱係数ＳW ＝１と仮定する。また、通常はベース釉の調合率ＣB を一定とし、白色顔料の調合率ＣW を変化させて混合物を作製するので、ＣB ＝１００（一定）と仮定する。この時、白色顔料の調合率ＣW は、いわゆる外添加調合率（ベース釉の重量を１００％とした時の顔料の重量％）となる。これらの仮定をおくと、上記数式４は次の数式５のように書き換えられる。
【００５６】
【数５】
2006030210.tif 000007
【００５７】
ベース釉は無色透明に近いが、完全な無色透明ではない。すなわち、数式５において、ベース釉の散乱係数ＳB は０でない正の小さな値である。従って、ベース釉と白色顔料の混合物の散乱係数ＳW'は、図２（Ａ）に示すようなグラフとなる。このグラフは、白色顔料の調合率ＣW が大きくなるとＳW'＝１に漸近する双曲線である。
【００５８】
ところで、図２（Ａ）に示すような散乱係数ＳW'を実測することは容易ではない。そこで、この第１実施例では、白色顔料の調合率ＣW に対する散乱係数ＳW'のグラフの形状を、図２（Ｂ）に示すように高濃度側の第１の直線Ｌ１と低濃度側の第２の直Ｌ２で近似できるものと仮定した。直線Ｌ１，Ｌ２を表わす関数ＳW'＝ｆ1（ＣW ），ＳW'＝ｆ2（ＣW ）は、それぞれ次の数式６および７で表わすことができる。
【００５９】
【数６】
2006030210.tif 000008
【数７】
2006030210.tif 000009
【００６０】
ここで，CW0 は直線Ｌ１，Ｌ２の交点における白色顔料の調合率である。
【００６１】
数式６，７における係数ａ，ｂ，ｄ，ｅを決定すれば、（ベース釉＋白色顔料）の混合物に対する散乱係数ＳW'が求まる。そして、散乱係数ＳW'が求まれば、分光反射率Ｒ（λ）の測定値を用いて上記数式２から吸収係数ＫW'を求めることができる。
【００６２】
図３は、図１のステップＳ１の詳細手順を示すフローチャートである。ステップＳ１１では、ベース釉と白色顔料のみを混合した混合物（第１の混合物）を作製した。図４は、第１実施例において準備した１３個のサンプルＷ１２〜Ｗ０の調合率を示している。図４から解るように、１３個のサンプルＷ１２〜Ｗ０は、白色顔料の調合率ＣW を１２％〜０％の範囲で１％ずつ変化させて作製したものである。なお、第１実施例におけるサンプルは、通常の陶磁器の素地に、調合した釉をかけて焼成することによって作製した。また、後述する他のサンプルも同じ条件で作製した。
【００６３】
ステップＳ１２では、各サンプルＷ１２〜Ｗ０について、分光光度計で分光反射率Ｒ'（λ）を測定した。ステップＳ１３では、分光反射率の測定値Ｒ'（λ）に基づいて、数式６に含まれる係数ａ，ｂを以下のように実験的に決定した。
【００６４】
白色顔料の調合率ＣW が比較的大きなサンプルでは、釉の隠ぺい力が大きいので、ＳW'＝１と仮定することができる。具体的には、白色顔料の調合率ＣW が最も大きなサンプルＷ１２に対する散乱係数ＳW'を１に等しいと仮定した。この仮定を用いると、白色顔料の調合率ＣW がｉ％のサンプルＷｉの散乱係数ＳWi'は、次の数式８で与えられる。
【００６５】
【数８】
2006030210.tif 000010
【００６６】
数式８において、（Ｋ／Ｓ）W12 は白色顔料の調合率が１２％のサンプルＷ１２の吸収係数と散乱係数の比であり、（Ｋ／Ｓ）Wiは白色顔料の調合率がｉ％のサンプルＷｉの吸収係数と散乱係数の比である。また、数式８では、吸収係数の変化は散乱係数の変化に比べて小さく、ＫW12 ＝ＫWiであると仮定している。
【００６７】
数式８の右辺の（Ｋ／Ｓ）W12 ，（Ｋ／Ｓ）WiはサンプルＷｉの分光反射率の測定値Ｒ'（λ）から、上述した数式２に従って求めることができる。但し、混合物である釉の厚みは小さいので、分光光度計で測定される分光反射率Ｒ'（λ）を、次の数式９（サンダーソン（Saunderson）の式）に従って理想状態（混合物の厚みが無限大の状態）の分光反射率Ｒ（λ）に変換する。
【００６８】
【数９】
2006030210.tif 000011
【００６９】
ここで、係数ｋ1 ，ｋ2 は、被着色物（ベース釉）の光学的性質に依存する値である。係数ｋ1 ，ｋ2 としては、被着色物の屈折率ｎから次の数式１０に従って決定することができる。
【００７０】
【数１０】
2006030210.tif 000012
【００７１】
なお、第１実施例において用いたベース釉は、屈折率ｎが約１．４である。こうして理想状態の分光反射率Ｒ（λ）が求まると、上述の数式２から（Ｋ／Ｓ）Wiが得られる。
【００７２】
数式９に従って分光反射率の測定値Ｒ'（λ）から理想状態の分光反射率Ｒ（λ）を求める演算は、後述する他の工程においてサンプルの分光反射率から（Ｋ／Ｓ）を求める際にも同様に実施される。
【００７３】
各サンプルＷｉについて得られた（Ｋ／Ｓ）Wiの値を上記数式８に代入することによって、各調合率ＣW において混合物の散乱係数ＳW'の値を求めることができ、係数ａ，ｂを決定することができる。係数ａ，ｂを含む数式６は、白色顔料の調合率ＣW が比較的高い範囲にのみ適用される。従って、係数ａ，ｂの決定においても、調合率ＣW が比較的高い範囲の（例えば１０％〜１２％の）数個のサンプルの（Ｋ／Ｓ）Wiのみを用いる。
【００７４】
図５は、ＣW＝１０％〜１２％のサンプルの分光反射率Ｒ（λ）を示す概念図である。このように、ベース釉に高濃度の白色顔料を混合したものは、可視光線の波長λの全範囲（約４００ｎｍ〜約７００ｎｍ）においてほとんど一定の反射率を示す。そこで、この分光反射率Ｒ（λ）の値から数式２と数式８を用いて散乱係数ＳWi'を求めることができる。なお、この散乱係数ＳWi'は、可視光線の全波長範囲において同一の値を有する。第１実施例では、ａ＝０．０００５，ｂ＝０．９４が得られたので、数式６が次の数式１１に書き換えられた。
【００７５】
【数１１】
2006030210.tif 000013
【００７６】
なお、数式１１は、ＣW ＝１０％付近においてＳW'＝１に極めて近い値を与えるので、数式１１の代わりにＳW'＝１を用いてもよい。
【００７７】
図３のステップＳ１４では、低濃度側の数式７の係数ｄ，ｅの値を決定する。白色顔料の調合率ＣW が比較的小さい場合には下地の影響が大きくなるので、分光反射率Ｒ（λ）の測定値が真の釉の層のみの情報ではなくなる。従って、分光反射率Ｒ（λ）の測定値に基づいて係数ｄ，ｅを決定するのは容易ではない。そこで、係数ｄ，ｅの値は、図６に示すような手順に従って決定する。
【００７８】
ステップＳ２１では、顔料物性値決定用サンプルを作製する。図７は、顔料決定用サンプルの調合率を示す説明図である。顔料物性値決定用サンプルは、ベース釉と白色顔料と他の有色顔料を混合して作製したものである。第１実施例では、白色顔料と有色顔料の合計の調合率（Pigment Volume Concentration，ＰＶＣ）を１２％一定とし、有色顔料の調合率ＣP を１％〜１２％の範囲で１％ずつ変えた１２個のサンプルＭ１〜Ｍ１２を作製した。なお、有色顔料としては、陶磁器の彩色に良く用いられる青色、赤色、黄色、等を選択し、各顔料について図７に示すような１２個のサンプルＭ１〜Ｍ１２をそれぞれ作製した。
【００７９】
なお、サンプルＭ１２は白色顔料を含まない混合物であり、本願発明の第２の混合物に相当する。また、サンプルＭ１〜Ｍ１１は本願発明の第３の混合物に相当する。
【００８０】
ステップＳ２２では、検証用サンプルを作製する。図８は、第１実施例で用いた検証用サンプルの調合率を示す説明図である。４つの検証用サンプルＤ１〜Ｄ４は、４色の顔料をベース釉に混合したものである。なお、ベース釉の調合率ＣB は省略しているが、いずれのサンプルもＣB ＝１００％である。
【００８１】
検証用サンプルＤ１〜Ｄ４は、本願発明における第４の混合物に相当する。
【００８２】
ステップＳ２３では、係数ｄ，ｅに仮の値を割り当てる。図２（Ａ），（Ｂ）を比較すれば解るように、係数ｅの値は、ベース釉の散乱係数ＳB の値に相当しており、かなり小さな値であることが解っている。例えば、ｅ＝０と仮定することも可能である。また、図２（Ｂ）に示す第１と第２の直線Ｌ１，Ｌ２の交点は、横軸の調合率ＣW が数％のところに位置するので、係数ｄも概略の値を仮決めすることが可能である。
【００８３】
係数ｄ，ｅを仮決めすると、調合率ＣW に依存する散乱係数ＳWi'が調合率ＣW の全範囲において決定される。散乱係数ＳWi'が決まると、次の数式１２に従って調合率ＣW に依存する吸収係数ＫWi'も得られる。
【００８４】
【数１２】
2006030210.tif 000014
【００８５】
なお、散乱係数ＳWi'の値は波長λに依存しないが、吸収係数ＫWi'の値は波長λに依存する。この理由は、この発明の相対法においては散乱係数ＳWi'を波長λに依存しない基準値として用い、他の物性値の相対値を求めるからである。分光反射率Ｒ（λ）の値はもちろん波長λに依存するので、数式２に従って得られる（Ｋ／Ｓ）Wiの値も波長λに依存する。従って、数式１２に従って得られる吸収係数ＫWi'も波長λに依存する。換言すれば、吸収係数ＫWi'は白色顔料の調合率ＣW と波長とに依存する形式で求められる。
【００８６】
ステップＳ２５では、顔料物性値決定用サンプルを用いて、各有色顔料の吸収係数ＫP（λ）と散乱係数ＳP（λ）とを求める。なお、ステップＳ２５は、図１におけるステップＳ２に相当する。ＫP（λ）とＳP（λ）は、次のような手順で算出される。
【００８７】
まず、顔料物性値決定用サンプルについては、上述の数式１から次の数式１３が得られる。
【００８８】
【数１３】
2006030210.tif 000015
【００８９】
ところで、上記の数式３は、次の数式１４のように書き換えられるので、数式１４を数式１３に代入すると、数式１５が得られる。
【００９０】
【数１４】
2006030210.tif 000016
【数１５】
2006030210.tif 000017
【００９１】
数式１５を変形すると、次の数式１６が得られる。
【００９２】
【数１６】
2006030210.tif 000018
【００９３】
ここで、（Ｋ／Ｓ）P は有色顔料単独での吸収係数と散乱係数の比であり、（Ｋ／Ｓ）M は顔料物性値決定用サンプルの吸収係数と散乱係数の比、（Ｋ／Ｓ）W'は（ベース釉＋白色顔料）のサンプルの吸収係数と散乱係数の比である。また、ＣW'はベース釉と白色顔料の調合率の合計値、ＳW'は（ベース釉＋白色顔料）の混合物の散乱係数であり、数式６および７で与えられる値である。数式１６の右辺の各項は、以下のように求めることができる。
【００９４】
（Ｋ／Ｓ）P の値は、白色顔料を含まない顔料物性値決定用サンプル（図７のサンプルＭ１２）の分光反射率Ｒ（λ）から上記数式２によって求めることができる。但し、厳密に言えば、数式１６における（Ｋ／Ｓ）P の値はベース釉の影響を含まない有色顔料単独に対する値である。一方、上述のようにして実測で得られる（Ｋ／Ｓ）P は（ベース釉＋有色顔料）のサンプルに対する値なので、ベース釉の影響が含まれている。しかし、図７のサンプルＭ１２の顔料の調合率は１２％と高い値なので、（Ｋ／Ｓ）P に対するベース釉の寄与は極めて小さい。従って、数式１６における（Ｋ／Ｓ）P の値としては、サンプルＭ１２の分光反射率Ｒ（λ）から数式２によって求めた値を用いても、その誤差は無視できる程度である。
【００９５】
（Ｋ／Ｓ）M の値は、白色顔料と有色顔料の両方を含む顔料物性値決定用サンプル（図７のサンプルＭ１〜Ｍ１１）の分光反射率Ｒ（λ）から数式２に従って算出できる。従って、（Ｋ／Ｓ）M は有色顔料の調合率ＣP に依存した形式で求められる。
【００９６】
（Ｋ／Ｓ）W' の値は、図４に示す（ベース釉＋白色顔料）のサンプルの分光反射率Ｒ（λ）から数式２に従って算出された値を用いることができる。
【００９７】
以上のようにして数式１６の右辺の各項の値を求めることができるので、顔料物性値決定用サンプルＭ１〜Ｍ１１のそれぞれについて、散乱係数ＳP（λ）を求めることができる。なお、散乱係数ＳP（λ）は、有色顔料の調合率ＣP に依存する形式で求められる。
【００９８】
こうして得られた散乱係数ＳP（λ）と、分光反射率Ｒ（λ）の測定値から数式２で算出された（Ｋ／Ｓ）P の値とに基づいて、次の数式１７に従って吸収係数ＫP（λ）も求められる。
【００９９】
【数１７】
2006030210.tif 000019
【０１００】
なお、吸収係数ＫP（λ）は、散乱係数ＳP（λ）と同様に、有色顔料の調合率ＣP に依存する形式で求められる。
【０１０１】
図９は、吸収係数ＫP（λ）の調合率ＣP に対する依存性の一例を示すグラフである。散乱係数ＳP（λ）もこれと同様な依存性を示す。なお、図９では、代表的な波長λについてのグラフのみを示しているが、実際には可視光線の全波長範囲（約４００ｎｍ〜約７００ｎｍ）で１０ｎｍ毎に吸収係数ＫP（λ）が求められる。
【０１０２】
以上のようにして、（ベース釉＋白色顔料）の混合物の吸収係数ＫW'（λ）および散乱係数ＳW'と、他の有色顔料の吸収係数ＫP（λ）および散乱係数ＳP（λ）が求まると、図６のステップＳ２６において、検証用サンプル（図８）に対してコンピュータカラーマッチングによるシミュレーションを行なう。
【０１０３】
ここで、コンピュータカラーマッチングによる色予測と調合割合の予測について簡単に説明する。数式２を変形すると、混合物の分光反射率Ｒ（λ）は、次の数式１８で与えられる。
【０１０４】
【数１８】
2006030210.tif 000020
【０１０５】
混合物の吸収係数と散乱係数の比（Ｋ／Ｓ）M は、数式１に従って算出できるので、混合物の分光反射率Ｒ（λ）も数式１８から求めることができる。この分光反射率Ｒ（λ）は理想状態の分光反射率なので、分光光度計で測定できる分光反射率Ｒ'（λ）を、数式９（サンダーソンの式）を変形した次の数式１９に従って求める。
【０１０６】
【数１９】
2006030210.tif 000021
【０１０７】
係数ｋ1 ，ｋ2 は上記数式１０で与えられる。
【０１０８】
分光反射率Ｒ'（λ）が求まると、混合物の三刺激値Ｘ，Ｙ，Ｚが次の数式２０によって求められる。
【０１０９】
【数２０】
2006030210.tif 000022
【０１１０】
ここで、Ｓ（λ）は標準光の分光分布、ｘ（λ），ｙ（λ），ｚ（λ）（数式中ではバー付きである）は等色関数である。
【０１１１】
コンピュータカラーマッチングによって混合物の色を予測する場合には、上述のように、混合物の各成分の吸収係数と散乱係数に基づき、数式１，数式１８〜２０に従って、その混合物の三刺激値Ｘ，Ｙ，Ｚを算出する。三刺激値は混合物の色を表わすので、任意の混合物の色を予測することができることになる。
【０１１２】
また、所望の色を有する混合物の調合割合を予測する場合には、混合物の調合割合を仮定して、上述した手順でその三刺激値を算出し、ニュートン−ラプソン法などの逐次近似法によって所望の色に所定の誤差内で一致するような調合割合を求める。
【０１１３】
なお、ＸＹＺ表色系の代わりにＬ*ａ*ｂ* 表色系のような他の表色系を用いてコンピュータカラーマッチングを行なうことも可能である。
【０１１４】
図６のステップＳ２６におけるシミュレーションでは、コンピュータカラーマッチングによって、図８に示す各検証用サンプルＤ１〜Ｄ４の調合率を予測する。検証用サンプルの調合率は既知なので、予測した調合率と実調合率との一致度は簡単に算出できる。例えば、調合率の一致度の指標として、次の数式２１で与えられる自乗平均誤差Δを使用する。
【０１１５】
【数２１】
2006030210.tif 000023
【０１１６】
ここでＣRiは成分ｉの実調合率、Ｃi は予測された調合率である。ステップＳ２６において、各検証用サンプルＤ１〜Ｄ４について得られた自乗平均誤差Δが所定の許容値以下でない場合にはステップＳ２３に戻り、係数ｄ，ｅの値を修正する。こうして、各検証用サンプルＤ１〜Ｄ４について得られる自乗平均誤差Δが所定の許容値以下になるまでステップＳ２３〜Ｓ２６を繰り返し実行し、この結果、係数ｄ，ｅを決定することができる（ステップＳ２７）。第１実施例では、ｄ＝０．１２２５，ｅ＝０となり、数式７が次の数式２２書き換えられた。
【０１１７】
【数２２】
2006030210.tif 000024
【０１１８】
なお、数式１１で与えられる高濃度側の散乱係数ＳW'のグラフ（図２（Ｂ）の直線Ｌ１）と、数式２２で与えられる低濃度側の散乱係数ＳW'のグラフ（直線Ｌ２）との交点における調合率ＣW は８．０％である。
【０１１９】
以上のように、上記第１実施例では検証用サンプルの調合率を精度良く予測できるように数式７における係数ｄ，ｅを決定したので、白色顔料の調合率ＣW が比較的低い領域においても、（ベース釉＋白色顔料）の混合物の散乱係数ＳW'と吸収係数ＫW'（λ）とを精度良く決定することができた。
【０１２０】
なお、図６に示す手順で（ベース釉＋白色顔料）の混合物の散乱係数ＳW'と吸収係数ＫW'（λ）とが求められる際には、ステップＳ２５において他の有色顔料の吸収係数ＫP（λ）と散乱係数ＳP（λ）も同時に求められる。
【０１２１】
こうして各成分の吸収係数と散乱係数が求まると、図１のステップＳ３において、コンピュータカラーマッチングによる色予測や調合割合の予測が実行される。図１０は、第１実施例における調合率の予測結果を示す説明図である。図１０における第１実施例の予測結果は、数式１１および数式２２を用いた場合の結果である。また、比較例の予測結果は、ＳW'＝１．０（一定）の仮定（図２（Ｂ）において破線で示すグラフ）を用いた場合の結果である。なお、図１０の第１実施例と比較例では、Ｌ*ａ*ｂ* 表色系を用いて色差が最小となる調合率を予測した。
【０１２２】
図１０に示されるように、第１実施例による予測結果は調合率の自乗平均誤差が比較例に比べてかなり小さく、より高精度で調合率を予測できたことが解る。特に、サンプルＰ２，Ｐ３のように白色顔料の調合率が小さな場合には、比較例の予測精度がかなり悪化しているのに対して、第１実施例の予測精度は良好である。この理由は、図２（Ａ），（Ｂ）に示すように、白色顔料の調合率ＣW が小さくなるに従って、実際の散乱係数ＳW'の値が１から離れていくので、ＳW'＝１（一定）と仮定した比較例では、各成分の吸収係数と散乱係数の値の誤差が大きくなるからであると考えられる。
【０１２３】
以上のように、この第１実施例では、無色透明でない被着色物であるベース釉の光学的影響を考慮して各成分の吸収係数と散乱係数とを決定しているので、コンピュータカラーマッチングによって高精度の予測を行なうことが可能である。また、サンプルとしては図４、図７、図８に示すような混合物を通常の素地の上に塗布したものを準備すればよいので、サンプルの作製も容易である。
【０１２４】
なお、上記第１実施例に関しては、例えば次のような変形も可能である。
【０１２５】
（１）散乱係数ＳW'の値を図２（Ｂ）に示される２本の直線Ｌ１，Ｌ２を数式６，７で近似する代わりに、数式５をそのまま用いて散乱係数ＳW'を決定してもよい。この場合には、図６のステップＳ２３においてベース釉の散乱係数ＳB の値を仮決定し、ステップＳ２３〜Ｓ２６を繰り返し実行することによって散乱係数ＳB の値を求めるようにすればよい。こうすれば、未知数が散乱係数ＳB の１つだけなので、ステップＳ２３において散乱係数ＳB を修正するための逐次近似法の計算の収束が早くなるという利点がある。
【０１２６】
ただし、ステップＳ２３では、逐次近似法によって係数を求める必要はなく、種々の係数の値を仮定して、ステップＳ２６におけるシミュレーション結果が最適となる係数を選択するようにしてもよい。例えば、数式５をそのまま用いる場合には、ＳB の値を約０．００１〜約０．０１０の範囲で数点設定しておき、ステップＳ２６において実調合率に最も近い調合率を予測できるＳB の値を決定するようにしてもよい。
【０１２７】
散乱係数ＳW'を表現するための関数は数式５〜７に限られるわけではなく、一般に、散乱係数ＳW'を調合率ＣW に依存した関数ｆ（ＣW ）で表現するようにすればよい。
【０１２８】
なお、上記第１実施例では、着色剤の調合率を外添加調合率（被着色物の調合率ＳB を１００％とした時の着色剤の調合率）で定義したが、調合率を他の規定の仕方で定義することも可能である。従って、数式５は、定数ａ1 を含む次の数式２３によって表現することができる。
【０１２９】
【数２３】
2006030210.tif 000025
【０１３０】
（２）ＸＹＺ表色系やＬ*ａ*ｂ* 表色系の値を一致させるコンピュータカラーマッチングは、一般にメタメリックマッチ法と呼ばれている。一方、混合物の分光反射率Ｒ（λ）の曲線を一致させるアイソメリックマッチ法と呼ばれる方法もある。この発明は、メタメリックマッチ法のみでなく、アイソメリックマッチ法によるコンピュータカラーマッチングにも適用可能である。
【０１３１】
（３）上記第１実施例では陶磁器の釉を対象とするコンピュータカラーマッチングについて説明したが、本発明はこれに限らず、他の種類の混合物を対象とするコンピュータカラーマッチングにも適用することが可能である。但し、陶磁器の素地は、完全な白色や完全な黒色のものを作製することは困難なので、絶対法によって被着色物（ベース釉）単独の物性値を決定することは難しい。従って、釉に関するコンピュータカラーマッチングに本発明を適用すれば、予測精度を向上させる上で特に効果が大きい。
【０１３２】
以下では，第１実施例で求められた散乱係数と吸収係数とを用いたコンピュータカラーマッチングの種々の適用例としての各種の実施例を説明する。但し，以下の各実施例は，第１実施例で決定された散乱係数と吸収係数とを用いる場合に限られるわけではなく，他の方法で決定された散乱係数と吸収係数とを用いることも可能であることに注意すべきである。
【０１３３】
Ｂ．第２実施例：
一般に、コンピュータカラーマッチングによる予測には誤差があるので、その予測誤差を小さくする工夫が必要である。従来は、予測誤差を減少させるために、各着色剤の吸収係数Ｋi と反射係数Ｓi とを正確な値に近づけるように補正していた。
【０１３４】
しかし、各成分の吸収係数Ｋi と反射係数Ｓi とを正確な値に近づけるためには、被着色物に各着色剤を単独で含む数多くの混合物サンプルを作製して、その分光反射率を測定しなければならず、膨大な作業を要していた。また、天然の顔料を着色剤として用いる場合には、その吸収係数Ｋi や散乱係数Ｓi が必ずしも一定の値にはならないので、吸収係数Ｋi と散乱係数Ｓi を正確な値に近づけることが困難であった。
【０１３５】
以下に示す第２実施例は、各成分の吸収係数Ｋi と散乱係数Ｓi とを補正することなく、予測誤差を減少させることを目的としたものである。
【０１３６】
図１１は、第２実施例における処理の全体手順を示すフローチャートである。なお、この第２実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉（ゆう）である。すなわち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラーマッチングの対象となる混合物である。
【０１３７】
ステップＳ３１では、コンピュータカラーマッチングの予測結果（三刺激値）を補正するためのニューラルネットワークの学習を行なう。ステップＳ３２〜Ｓ３５では、学習済みのニューラルネットワークを用いてコンピュータカラーマッチングの目標値を補正して、正確な予測結果を求めている。以下ではまず、ステップＳ３２〜Ｓ３５の内容を説明する前に、ニューラルネットワークの構成とステップＳ３１の詳細手順について説明する。
【０１３８】
図１２は、ニューラルネットワークの構成を示す説明図である。このニューラルネットワークは、入力層１０と中間層２０と出力層３０とで構成される三層の階層構造を有している。入力層１０は、３つのニューロンＮ11〜Ｎ13で構成されており、中間層２０は５つのニューロンＮ21〜Ｎ25で、出力層３０は３つのニューロンＮ31〜Ｎ33でそれぞれ構成されている。
【０１３９】
入力層１０の３つのニューロンＮ11〜Ｎ13には、三刺激値Ｘ，Ｙ，Ｚがそれぞれ入力される。入力層１０のニューロンＮijから中間層２０のニューロンＮk に伝達される信号は、それぞれの入力信号に重みＷij,k を乗じたものである。ここで、ｉは注目している階層を示す番号、ｊは注目している階層内でのニューロンの順番を示す番号、ｋは次の階層のニューロンの順番を示す番号である。例えば、入力層１０の第１のニューロンＮ11から中間層２０の第１のニューロンＮ21に伝達される信号はＷ11,1Ｘであり、入力層１０の第１のニューロンＮ11から中間層２０の第２のニューロンＮ22に伝達される信号はＷ11,2Ｘである。
【０１４０】
中間層２０の各ニューロンＮijの入力ｕijと出力Ｑijとの関係は、次の数式２４に示す情報伝達関数ｆ（ｕij）で与えられる。
【０１４１】
【数２４】
2006030210.tif 000026
【０１４２】
ここで、Ｑ(i-1)jは（ｉ−１）番目の階層（すなわち入力層１０）のｊ番目のニューロンＮ(i-1)jの出力であり、図１２の例ではＱ11＝Ｘ，Ｑ12＝Ｙ，Ｑ13＝Ｚである。また、Ｗ(i-1)j,k は（ｉ−１）番目の階層のｊ番目のニューロンＮ(i-1)jから、注目しているニューロンＮijに伝達される信号に掛かる重みである。ｔはしきい値であり、一定の値が割当てられる。なお、数式２４の情報伝達関数ｆ（ｕij）はシグモイド関数と呼ばれている。
【０１４３】
例えば、数式２４を中間層２０の第１のニューロンＮ21の入出力関係に適用すると、次の数式２５が得られる。
【０１４４】
【数２５】
2006030210.tif 000027
【０１４５】
出力層３０の各ニューロンの入出力関係も上記数式２４で与えられる。図１２に示す第２実施例では、出力層３０の３つのニューロンＮ31〜Ｎ33の出力Ｑ11〜Ｑ13を、コンピュータカラーマッチング（ＣＣＭ）による三刺激値の予測誤差ΔＸ，ΔＹ，ΔＺとしている。
【０１４６】
ニューラルネットワークの学習は、入力層１０への入力（Ｘ，Ｙ，Ｚ）と出力層３０からの出力（ΔＸ，ΔＹ，ΔＺ）との間の関係を数多く与えて、正しい入出力関係を与えるような重みＷij,k の値を決定する作業である。
【０１４７】
図１３は、図１１のステップＳ３１の詳細手順を示すフローチャートである。ステップＳ４１では、コンピュータカラーマッチングにおいて予測の対象とする混合物の色の範囲をカバーするような複数の三刺激値（Ｘit，Ｙit，Ｚit）を定めて、これらの複数の三刺激値を有するような複数のサンプルの調合割合をコンピュータカラーマッチングで決定する。図１４は、コンピュータカラーマッチングでの色の予測対象範囲ＰＡと、この予測対象範囲ＰＡをカバーする複数の三刺激値の分布を示す概念図である。この第２実施例ではＣＩＥ−ＸＹＺ表色系で色を表現するものとしており、色の予測対象範囲ＰＡはＸＹＺ座標系の３次元的な範囲として与えられる。なお、予測対象範囲ＰＡは、予測対象とする混合物が取り得る色の範囲を示すものであり、任意に設定し得る範囲である。
【０１４８】
この第２実施例においては、予測対象範囲ＰＡをカバーするために、図１４に白丸で示す７組の三刺激値を決定した。ステップＳ４１では、さらに、これらの７組の三刺激値（Ｘit，Ｙit，Ｚit）を有するような７種類のサンプルＭ１〜Ｍ７の調合割合を、コンピュータカラーマッチングによって予測した。
【０１４９】
ここで、コンピュータカラーマッチングによる色予測と調合割合の予測について簡単に説明する。第１実施例で説明した数式２を変形すると、混合物の分光反射率Ｒ（λ）は、次の数式２６で与えられる。
【０１５０】
【数２６】
2006030210.tif 000028
【０１５１】
混合物の吸収係数と散乱係数の比（Ｋ／Ｓ）Mは、各成分の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）と調合率Ｃi から，第１実施例で説明した数式１に従って算出できるので、混合物の分光反射率Ｒ（λ）を上記数式２６から求めることができる。この分光反射率Ｒ（λ）は理想状態（被着色物の厚みが無限大の場合）の分光反射率なので、分光光度計で測定できる分光反射率Ｒ'（λ）を、次の数式２７（サンダーソンの式）に従って求める。
【０１５２】
【数２７】
2006030210.tif 000029
【０１５３】
ここで、係数ｋ1 ，ｋ2 は、被着色物（ベース釉）の光学的性質に依存する値である。係数ｋ1 ，ｋ2 としては、被着色物の屈折率ｎから次の数式２８に従って決定することができる。
【０１５４】
【数２８】
2006030210.tif 000030
【０１５５】
なお、第２実施例において用いたベース釉は、屈折率ｎが約１．４である。
【０１５６】
数式２７によって分光反射率Ｒ'（λ）が求まると、混合物の三刺激値Ｘ，Ｙ，Ｚが次の数式２９によって求められる。
【０１５７】
【数２９】
2006030210.tif 000031
【０１５８】
ここで、Ｓ（λ）は標準光の分光分布、ｘ（λ），ｙ（λ），ｚ（λ）（数式中ではバー付きである）は等色関数である。
【０１５９】
コンピュータカラーマッチングによって混合物の色を予測する場合には、上述のように、混合物の各成分の吸収係数と散乱係数に基づき、数式１，２６〜２９に従って、その混合物の三刺激値Ｘ，Ｙ，Ｚを算出する。三刺激値は混合物の色を表わすので、任意の混合物の色を予測することができることになる。
【０１６０】
また、所望の色を有する混合物の調合割合を予測する場合には、混合物の調合割合を仮定して、上述した手順でその三刺激値を算出し、ニュートン−ラプソン法などの逐次近似法によって所望の色に所定の誤差内で一致するような調合割合を求める。
【０１６１】
図１３のステップＳ４２では、ステップＳ４１で予測された調合割合Ｃｉを有する複数のサンプルを作製する。この第２実施例では、図１４の７組の三刺激値（Ｘit，Ｙit，Ｚit）に対応する７つのサンプルＭ１〜Ｍ７を作製した。ステップＳ４３では、各サンプルＭｉの分光反射率を分光光度計で測定し、上記数式２９に従ってその三刺激値（Ｘim，Ｙim，Ｚim）を求める。
【０１６２】
ステップＳ４４では、各サンプルＭｉに関して、ステップＳ４３で得られた実測値（Ｘim，Ｙim，Ｚim）とステップＳ４１で決定した目標値（Ｘit，Ｙit，Ｚit）との差をとることによって、予測誤差ΔMi（Ｘim−Ｘit，Ｙim−Ｙit，Ｚim−Ｚit）を求める。図１５は、第２実施例において得られた７つのサンプルＭ１〜Ｍ７の予測誤差を示す概念図である。また、図１６は、各サンプルの三刺激値の目標値（ステップＳ４１で決定された値）と予測誤差ΔMi（ΔＸ，ΔＹ，ΔＺ）を示す説明図である。
【０１６３】
なお、ステップＳ４１で得られた調合割合Ｃiからコンピュータカラーマッチングによって各サンプルの三刺激値の予測値（Ｘic，Ｙic，Ｚic）を求め、実測値（Ｘim，Ｙim，Ｚim）と予測値（Ｘic，Ｙic，Ｚic）との差（Ｘim−Ｘic，Ｙim−Ｙic，Ｚim−Ｚic）を予測誤差ΔMiと定義してもよい。コンピュータカラーマッチングで調合割合Ｃｉを決定する際には、三刺激値の目標値（Ｘit，Ｙit，Ｚit）と予測値（Ｘic，Ｙic，Ｚic）との差が所定の許容誤差以下になるように調合割合Ｃiを決定するので、目標値（Ｘit，Ｙit，Ｚit）と予測値（Ｘic，Ｙic，Ｚic）は実質的にほぼ等しい値を有している。従って、実測値（Ｘim，Ｙim，Ｚim）と予測値（Ｘic，Ｙic，Ｚic）との差（Ｘim−Ｘic，Ｙim−Ｙic，Ｚim−Ｚic）を予測誤差ΔMiと定義しても、実測値（Ｘim，Ｙim，Ｚim）と目標値（Ｘit，Ｙit，Ｚit）との差（Ｘim−Ｘit，Ｙim−Ｙit，Ｚim−Ｚit）を予測誤差ΔMiと定義しても実質的にはほぼ同じである。
【０１６４】
図１３のステップＳ４５では、図１６に示す各サンプルＭｉの三刺激値の目標値（Ｘit，Ｙit，Ｚit）と予測誤差ΔMiとを用いてニューラルネットワークの学習を行ない、上記数式２４における重みＷij,k を決定する。ニューラルネットワークの学習方法としては、例えば逆誤差伝搬学習方式を用いる。
【０１６５】
図１７は、第２実施例におけるニューラルネットワークの学習の実証結果を示す説明図である。ここでは、図１４，１５に示す予測対象範囲ＰＡに含まれるもう１組の三刺激値をコンピュータカラーマッチングの目標値として設定し、この目標値を有する第８のサンプルＭ８を作製した。そして、このサンプルＭ８の三刺激値を実測した。図１７において「ＣＣＭ目標値」とあるのは第８のサンプルＭ８のコンピュータカラーマッチングに用いた目標値を意味している。また、「誤差（真値）」とあるのは三刺激値の目標値と実測値との差である。「ニューロ予測誤差」は、ＣＣＭ目標値を学習済みのニューラルネットワーク（図１２）に入力した場合に得られる予測誤差である。図１７の結果から、学習済みのニューラルネットワークは、三刺激値の誤差を精度良く予測できることが解る。
【０１６６】
こうしてニューラルネットワークの学習が終了すると、図１１のステップＳ３２〜Ｓ３５を実行して、調合率未知の色見本の調合率を予測する。ステップＳ３２では、調合率未知の色見本の分光反射率を測定してその三刺激値（Ｘｓ，Ｙｓ，Ｚｓ）を求める。ステップＳ３３では、図１２に示すニューラルネットワークに色見本の三刺激値（Ｘｓ，Ｙｓ，Ｚｓ）を入力して、予測誤差Δ（ΔＸｓ，ΔＹｓ，ΔＺｓ）を求める。テップＳ３４では、三刺激値（Ｘｓ，Ｙｓ，Ｚｓ）を予測誤差Δで補正して、コンピュータカラーマッチングの調合割合予測に用いる三刺激値の目標値（Ｘｓ−ΔＸｓ，Ｙｓ−ΔＹｓ，Ｚｓ−ΔＺｓ）を求める。ステップＳ３５では、こうして補正された目標値を用いてコンピュータカラーマッチングを実行し、色見本の調合率を予測する。
【０１６７】
図１８は、第２実施例におけるコンピュータカラーマッチングの予測精度を検証するために行なった実験結果を示す説明図である。ここでは、予測精度を実証することを目的としたので、調合率既知の色見本について三刺激値を実測し、その実測値（Ｘｓ，Ｙｓ，Ｚｓ）を実現する調合率をコンピュータカラーマッチングで予測した。図１７の結果から、調合率の真値とコンピュータカラーマッチングによる予測値とは極めて良く一致していることが解る。
【０１６８】
なお、図１２に示す学習済みのニューラルネットワークを利用すれば、調合率既知の混合物の三刺激値も精度良く予測することが可能である。すなわち、予測対象の混合物の調合率Ｃiから上記数式１，２６〜２９を用いて三刺激値を求め、この三刺激値を図１２のニューラルネットワークに入力して予測誤差Δを求める。そして、数式２９で得られた三刺激値を予測誤差Δで補正すれば、実際の値に極めて近い三刺激値が得られる。
【０１６９】
プラスチック等と異なり，陶磁器においては，焼成条件，熔融中の原料の化学反応等により，混色理論が成立しない場合がある。この場合，ニューラルネットワークに混合物の三刺激値とＣＣＭによる予測誤差の関係を学習させる本実施例の方法は有効である。
【０１７０】
なお、上記第２実施例に関しては、例えば次のような変形も可能である。
【０１７１】
（１）上記第２実施例ではコンピュータカラーマッチングにおける三刺激値の目標値をニューラルネットワークの入力としていたが、ニューラルネットワークの入力としては、各サンプルＭｉの三刺激値の実測値（Ｘim，Ｙim，Ｚim）を用いても良く、また、コンピュータカラーマッチングによる予測値（Ｘic，Ｙic，Ｚic）を用いても良い。すなわち、ニューラルネットワークには、複数のサンプルの三刺激値（ＸＹＺ表色系の座標値）とその予測誤差との関係を学習させるようにすればよい。
【０１７２】
（２）上記第２実施例では、ニューラルネットワークを用いて三刺激値を補正していたが、回帰分析やニューロファジイ技術等の他の誤差補正法を用いて補正することも可能である。
【０１７３】
（３）上記第２実施例ではＸＹＺ表色系を用いていたが、本発明は、Ｌ*ａ*ｂ* 表色系などの他の任意の表色系を用いる場合に適用できる。
【０１７４】
（４）ＸＹＺ表色系やＬ*ａ*ｂ* 表色系の座標値を一致させるコンピュータカラーマッチングは、一般にメタメリックマッチ法と呼ばれている。一方、混合物の分光反射率Ｒ（λ）の曲線を一致させるアイソメリックマッチ法と呼ばれる方法もある。この発明は、メタメリックマッチ法のみでなく、アイソメリックマッチ法によるコンピュータカラーマッチングにも適用可能である。
【０１７５】
（５）上記第２実施例では陶磁器の釉を対象とするコンピュータカラーマッチングについて説明したが、本発明はこれに限らず、他の種類の混合物を対象とするコンピュータカラーマッチングにも適用することが可能である。
【０１７６】
Ｃ．第３実施例：
衛生陶器などのような工業陶器には、ユーザの様々な好みに応じて種々の色が設定されている。陶器の色は、釉に混合する顔料の調合率によって決定されるので、釉は、コンピュータカラーマッチングにおける混合物に相当する。陶器の色を設定する際には、まず、設計者が紙などを彩色したり、予め焼成されたサンプルから所望の色を有するサンプルを選択したりすることによって色を決定する。そして、その色の分光反射率を分光光度計で測定し、分光反射率の実測値からコンピュータマッチングを用いて顔料や染料等の着色剤の調合割合を予測する。
【０１７７】
ところで、同じ調合割合の釉薬を用いても、製造される陶器の実際の色は同一にはならず、色にばらつきがあるのが普通である。図１９は、設計された色と実際に製造される陶器の色のばらつきを示するｘ−ｙ色度図である。図１９において、二重丸で示される色が設計された所望の色（標準色）Ｌ０であり、小さな白丸が実際に製造される陶器の色の分布である。標準色Ｌ０から大きく異なる色を有する陶器は不良品として認識する必要がある。そこで、良品の色の濃淡の限度を示すものとして、図１９に黒丸で示す色を有する濃淡限度色Ｌ１，Ｌ２が設定される。標準色Ｌ０の色度座標値は設計者によって決定されているので、この標準色Ｌ０の色度座標値に製造誤差等を考慮して、２つの濃淡限度色Ｌ１，Ｌ２の色度座標値を決定する。
【０１７８】
実際の陶器の製造においては、標準色Ｌ０を有する陶器を標準色見本として作製し、また、濃淡限度色Ｌ１，Ｌ２を有する２つの陶器を濃淡限度見本として作製する。そして、検査工程において、標準色見本および濃淡限度見本の色と製造された陶器の色とを比較して、濃淡限度見本の色の範囲にある陶器のみを良品とする。
【０１７９】
ところで、コンピュータカラーマッチングによる予測にはかなりの誤差が伴うので、予測された調合割合から標準色見本や濃淡限度見本を作製するにはかなりの手間を要するのが普通である。例えば、標準色見本を作製する際には、コンピュータカラーマッチングで予測された調合割合で調合した釉薬を用いて多数のサンプルを焼成し、その中で所望の標準色を有するサンプルがあればそれを標準色見本として採用する。しかし、標準色を有するサンプルが無ければ、勘と経験に従って調合割合を変更して再度多数のサンプルを焼成し、標準色を有するサンプルが得られるまでこの作業を繰り返す。このような作業は、濃淡限度色Ｌ１，Ｌ２を有する濃淡限度見本についても行なわれる。
【０１８０】
このように、従来のコンピュータカラーマッチングでは、所望の色を有する混合物を作製するための調合割合を精度良く予測するのが困難であった。
【０１８１】
以下に示す第３実施例は、所望の色を有する混合物の調合割合を精度良く予測することを目的としたものである。
【０１８２】
図２０は、第３実施例における処理の全体手順を示すフローチャートである。なお、この第３実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉（ゆう）である。すなわち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラーマッチングの対象となる混合物である。
【０１８３】
ステップＳ５１では、まず標準色見本を作製する。従って、標準色見本の着色剤の混合割合は既知である。ステップＳ５２では、作製された標準色見本の分光反射率Ｒ'を分光光度計で測定し、この分光反射率Ｒ'から次の数式３０に従って三刺激値Ｘ0 ，Ｙ0 ，Ｚ0 を算出する。
【０１８４】
【数３０】
2006030210.tif 000032
【０１８５】
ここで、Ｓ（λ）は標準光の分光分布、ｘ（λ），ｙ（λ），ｚ（λ）（数式中ではバー付きである）は等色関数である。
【０１８６】
ステップＳ５３では、ステップＳ５２で得られた標準色の三刺激値の実測値Ｘ0 ，Ｙ0 ，Ｚ0 に基づいて、製造誤差を考慮して、設計者が図１９に示す２つの濃淡限度色Ｌ１，Ｌ２に対する三刺激値を設定する。なお、以下では説明の便宜上、第１の濃淡限度色Ｌ１を実現するための調合割合を予測する場合について説明する。第２の濃淡限度色Ｌ２に対しても、同様な処理によってその調合割合を精度良く予測することができる。
【０１８７】
ステップＳ５４では、ステップＳ５１で作製した標準色見本の各着色剤の調合率Ｃi（ｉは着色剤の番号を示す）と、散乱係数Ｓi と、吸収係数Ｋi とから、標準色見本の三刺激値の計算値ＸM ，ＹM ，ＺM を以下の手順で算出する。
【０１８８】
まず、第１実施例で説明した数式２を変形すると、混合物の分光反射率Ｒ（λ）は、次の数式３１で与えられる。
【０１８９】
【数３１】
2006030210.tif 000033
【０１９０】
混合物の吸収係数と散乱係数の比（Ｋ／Ｓ）Mは、各着色剤の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）と調合率Ｃi から、第１実施例で説明した数式１に従って算出できる。従って、混合物の分光反射率Ｒ（λ）は、この比（Ｋ／Ｓ）M から、数式３１に従って求めることができる。この分光反射率Ｒ（λ）は理想状態（被着色物の厚みが無限大の場合）の分光反射率なので、分光光度計で測定できる分光反射率Ｒ'（λ）を、次の数式３２（サンダーソンの式）に従って求める。
【０１９１】
【数３２】
2006030210.tif 000034
【０１９２】
ここで、係数ｋ1 ，ｋ2 は、被着色物（ベース釉）の光学的性質に依存する値である。係数ｋ1 ，ｋ2 の値は、被着色物の屈折率ｎから次の数式３３に従って決定することができる。
【０１９３】
【数３３】
2006030210.tif 000035
【０１９４】
なお、ベース釉の屈折率ｎは、例えば約１．４である。
【０１９５】
数式３２によって分光反射率Ｒ'（λ）が求まると、混合物の三刺激値ＸM ，ＹM ，ＺM は、前述した数式３０と同様な式に従って算出できる。
【０１９６】
ところで、各着色剤の吸収係数Ｋi と散乱係数Ｓi は誤差を含んでおり、また、数式１，３２〜３３は理論式ではなく実験式なので、ステップＳ５４で得られた三刺激値の計算値ＸM ，ＹM ，ＺM は誤差を有している。この計算誤差は、ステップＳ５２で得られた実測値Ｘ0 ，Ｙ0 ，Ｚ0 からの差である。ところで、コンピュータカラーマッチングによる調合割合の予測も、数式１，２，３０〜３３に従って混合物の三刺激値を求める工程を含んでいる。そこで、濃淡限度見本の調合割合をコンピュータカラーマッチングで予測する際に、標準色見本における三刺激値の計算誤差を考慮しておけば、その予測精度を向上させることができる。
【０１９７】
図２０のステップＳ５５では、ステップＳ５２で得られた標準色の三刺激値の実測値Ｘ0 ，Ｙ0 ，Ｚ0 と、ステップＳ５４で得られた計算値ＸM ，ＹM ，ＺM との差分ΔＸ，ΔＹ，ΔＺを次の数式３４に従って求める。
【０１９８】
【数３４】
2006030210.tif 000036
【０１９９】
ステップＳ５６では、ステップＳ５３で設定された濃淡限度見本の三刺激値ＸT ，ＹT ，ＺT を上記の差分ΔＸ，ΔＹ，ΔＺで補正することによって、コンピュータカラーマッチングにおける三刺激値の目標値ＸC ，ＹC ，ＺC を求める。すなわち、濃淡限度見本に対する三刺激値の目標値ＸC ，ＹC ，ＺC は次の数式３５で与えられる。
【０２００】
【数３５】
2006030210.tif 000037
【０２０１】
ステップＳ５７では、コンピュータカラーマッチングによって濃淡限度見本の調合割合を予測する。このコンピュータカラーマッチングでは、数式３５で与えられる目標値ＸC ，ＹC ，ＺC が得られるような調合割合を求めている。
【０２０２】
図２１は、ステップＳ５７におけるコンピュータカラーマッチングの詳細手順を示すフローチャートである。なお、図２１の手順は、ニュートン−ラプソン法を用いたコンピュータカラーマッチングを適用したものである。
【０２０３】
ステップＳ６１では、各着色剤（顔料）の調合率Ｃi を微小変化させた時の三刺激値の変化を以下の手順で計算する。まず、図２０のステップＳ５２で用いた標準色見本の調合率Ｃi から、１つの着色剤の調合率のみをΔＣi だけ変えて、その混合物の三刺激値の計算値ＸM'，ＹM'，ＺM'を求める。これは、前述した数式１，２，３０〜３３を用いて、ステップＳ５２と同様の手順で行なわれる。そして、ステップＳ５２で得られていた標準色の三刺激値の計算値ＸM ，ＹM ，ＺM との差分ΔＸCi，ΔＹCi，ΔＺCiを以下の数式３６に従って算出する。
【０２０４】
【数３６】
2006030210.tif 000038
【０２０５】
数式３６を用いると、各着色剤の調合率Ｃi のみを微小変化させた時の三刺激値の変化率が、次の数式３７で与えられる。
【０２０６】
【数３７】
2006030210.tif 000039
【０２０７】
濃淡限度色は、標準色の製造誤差を示すものなので、両者の三刺激値は互いに近接した値である。従って、数式３７の変化率を用いると、図２０のステップＳ５６で得られた濃淡限度見本の三刺激値の目標値ＸC ，ＹC ，ＺC と、ステップＳ５４で得られた標準色見本の三刺激値の計算値ＸM ，ＹM ，ＺM との関係を、次の数式３８で表わすことができる。
【０２０８】
【数３８】
2006030210.tif 000040
【０２０９】
なお、数式３８では、着色剤が４種類であると仮定している。なお、各着色剤の調合率の変化量ΔＣi の合計値ΔＰＣＶには一定値が代入される。例えば、４種類の着色剤の調合率Ｃi の和ΣＣi を一定に保った場合には、ΔＰＣＶは０である。
【０２１０】
数式３８は、４つの未知数ΔＣi（ｉ＝１〜４）を含む４元連立一次方程式なので、これを解くことによって、各未知数ΔＣi の値を求めることができる（ステップＳ６２）。ステップＳ６３では、ステップＳ６２で求められた値ΔＣi を用いて、濃淡限度見本の各着色剤の調合率ＣiTが次の数式３９に従って算出される。
【０２１１】
【数３９】
2006030210.tif 000041
【０２１２】
ステップＳ６４では、ステップＳ６３で得られた調合率ＣiTから、前述した数式１，２，３０〜３３に従って三刺激値の計算値Ｘ1 ，Ｙ1 ，Ｚ1 を求める。ステップＳ６５では、こうして得られた三刺激値Ｘ1 ，Ｙ1 ，Ｚ1 と、図２０のステップＳ５６で求めた三刺激値の目標値ＸC ，ＹC ，ＺC との色差ΔＥが所定の許容誤差δ以下であるか否かが判断される。なお、色差ΔＥは、Ｌａｂ表色系を用いて次の数式４０で与えられる。
【０２１３】
【数４０】
2006030210.tif 000042
【０２１４】
ステップＳ６５において、色差ΔＥが許容誤差δよりも小さな場合にはステップＳ６３で得られた調合率ＣiTの値を予測値として採用し、コンピュータカラーマッチングを終了する。一方、色差ΔＥが許容誤差δ以上の場合には、ステップＳ６６において、ステップＳ６４で得られた三刺激値の計算値Ｘ1 ，Ｙ1 ，Ｚ1 を、ＸM ，ＹM ，ＺM に代入して、ステップＳ６２〜Ｓ６５の処理を繰り返す。こうしてステップＳ６２〜Ｓ６６を繰り返し実行することによって、色差ΔＥの値が許容誤差δよりも小さくなるような調合率ＣiTを求めることができる。
【０２１５】
図２２は、本発明の第３実施例において用いた標準色見本の三刺激値と、濃淡限度見本の三刺激値を示すテーブルである。サンプルＴ１，Ｔ２，Ｔ３は、それぞれ異なる色を有する標準色見本である。図２２には、各標準色見本について、図２０のステップＳ５２で実測された三刺激値Ｘ0 ，Ｙ0 ，Ｚ0 と、ステップＳ５４で計算された三刺激値ＸM ，ＹM ，ＺM とが示されており、また、ステップＳ５３で設定された濃淡限度見本の三刺激値の設定値ＸT ，ＹT ，ＺT も示されている。コンピュータカラーマッチングにおける濃淡限度見本の三刺激値の目標値ＸC ，ＹC ，ＺC は、これらの値から、前述した数式３４および３５に従って算出される。
【０２１６】
図２３は、第３実施例における濃淡限度見本の三刺激値の予測結果と、比較例の予測結果と、実調合率とを示したものである。比較例の調合率は、濃淡限度見本の三刺激値の設定値ＸT ，ＹT ，ＺT をそのままコンピュータカラーマッチングの目標値とした時に得られる予測値である。また、第３実施例の調合率は、標準色見本の計算誤差で補正した目標値ＸC ，ＹC ，ＺC を用いて得られる予測値である。なお、この第３実施例では、予測精度を検証するために、図２３の右端に示す調合率を用いて濃淡限度見本を実際に製作しておき、その三刺激値の実測値を図２２に示す設定値ＸT ，ＹT ，ＺT として用いている。
【０２１７】
図２３から解るように、第３実施例は比較例に比べてより高い精度で調合率を予測することが可能であった。また、３つのサンプルＴ１〜Ｔ３では４つの顔料の成分がかなり異なるが、そのすべての場合において、第３実施例の予測精度が比較例に比べて高いことが解る。
【０２１８】
以上のように、上記第３実施例では、濃淡限度見本に近い三刺激値を有する標準色見本に関して三刺激値の計算誤差を求め、これを用いて濃淡限度見本の三刺激値の目標値を補正している。この結果、濃淡限度見本の調合率を予測する際の予測精度を向上させることが可能である。
【０２１９】
上記第３実施例に関しては、例えば次のような変形も可能である。
【０２２０】
（１）上記第３実施例では、濃淡限度見本に近い三刺激値を有するサンプル（近接色サンプル）として標準色見本を使用し、標準色見本に関する三刺激値の計算誤差を用いて濃淡限度見本の三刺激値の目標値を補正していた。しかし、標準色見本以外の近接色サンプルを選択し、その近接色サンプルに関する計算誤差を用いることも可能である。例えば、種々の色を有するサンプルを収集したデータべースを用いて、濃淡限度見本に近い色を有するサンプルを、近接色サンプルとして選択するようにしてもよい。このデータベースは、各着色剤の調合率Ｃi と、三刺激値の実測値（または、反射率Ｒ（λ））を少なくとも含むようにするのが好ましい。なお、データベースを用いて近接色サンプルを検索する場合には、上述の数式４０で与えられる色差ΔＥが、目標とする混合物（濃淡限度見本）に最も近いものを近接色サンプルとすることが好ましい。
【０２２１】
なお、このようなデータベースを用いれば、濃淡限度見本の調合率をコンピュータカラーマッチングによって予測する場合に限らず、任意の三刺激値を有する目標混合物の調合率を予測する場合に、その予測精度を向上させることができるという利点がある。
【０２２２】
（２）上記第３実施例では、ＸＹＺ表色系を用いていたが、表色系としてはＸＹＺ表色系以外の任意の表色系を使用することができる。例えば、Ｌ*ａ*ｂ* 表色系などを使用することが可能である。
【０２２３】
Ｄ．第４実施例：
顔料や染料などの着色剤の着色対象となるもの（被着色物）を所望の色とするには、その所望の色となるように着色剤を種々調合して調合割合を規定する必要がある。この場合、技術者が試行錯誤して所望の色を呈する着色剤の混合割合を求めていたが、近年では、個々の着色剤の濃度別の光学的データを基礎として、着色剤の混合割合を予測するいわゆるコンピュータカラーマッチング（以下、これを適宜ＣＣＭと略称する）が提案されている（特開平４−１８１１２９，特公平６−９８８８０）。
【０２２４】
ところで、これら公報で提案されたＣＣＭでは、ニュートン・ラプソン法等の逐次近似法を用いているため、その予測した調合割合が負となることがある。このため、負となった予測調合割合をゼロ若しくは正とするための工夫がなされている。
【０２２５】
上記した従来のＣＣＭでは、その被着色物を繊維とし、一旦着色（染色）した繊維の色が所望のものではない場合に、この染色済みの繊維を追加染色してその色を所望のものとする。しかしながら、陶器やタイルでは、その色は、釉に混合する顔料の調合率（調合割合）によって決定され、顔料が調合された釉（以下、釉薬という）の焼成を経て当該色を呈する。このため、陶器やタイル等のように追加染色（着色）できない被着色物には、上記の従来のＣＣＭは適用できない。
【０２２６】
また、追加着色が可能な繊維や追加着色が不可能な陶器，タイルであっても、共に、着色剤（繊維にあっては染料，陶器やタイルにあっては釉薬）は、繊維，陶器等の工業的生産に合わせて多量に繰り返し調合される。この場合、着色剤の調合割合は維持されるが、調合工程に変動、例えば温度や調合のタイミング等のズレがあったりすると、それ以前に調合した着色剤で得られる色や色見本の色を再現できない場合がある。特に、陶器やタイルでは、天然の顔料を用いる都合上、色の再現の信頼性にやや欠ける。
【０２２７】
このような場合、追加着色が可能な繊維では、上記したＣＣＭにより追加染色して色修正ができるが、追加染色するために調合した着色剤と既に調合済みの着色剤を併用しなければならず、煩雑である。この煩雑さを回避するためには、色見本の色を得られる着色剤を新たに調合してこの着色剤だけで染色すればよいが、既に調合済みの着色剤は不要として廃棄するか、この調合済みの着色剤に新たな着色剤を追加して再調合する必要がある。また、陶器やタイルでは、追加着色が不可能な都合上、調合済みの着色剤の廃棄か再調合を採ることになる。しかしながら、調合済み着色剤の廃棄は無駄であり、その一方、着色剤の再調合は技術者の長年の勘と経験により着色剤を徐々に追加しながら行なうのでやはり煩雑であった。
【０２２８】
以下に説明する第４実施例は、調合済み着色剤の有効利用を図りつつ、再調合を簡略化することを目的としたものである。
【０２２９】
図２４は、実施例のコンピュータカラーマッチング方法を実施するための装置を示すブロック図である。この装置は，上述した第１〜第３実施例を実現する装置としても使用することが可能である。演算装置４０は汎用コンピュータであり，図示しないＣＰＵがソフトウェアプログラムを実行することによって，この発明によるコンピュータカラーマッチングの各工程および各手段を実現する。
【０２３０】
図２５は、この第４実施例における処理の全体手順を示すフローチャートである。なお、この第４実施例で対象とする調合物は、陶磁器の素地の表面を覆うための釉である。即ち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラーマッチングの対象となる調合物である。また、その表色系は、ＸＹＺ表色系とするが、これ以外の表色系、例えばＬ*ａ*ｂ* 表色系を採用してもよいことは勿論である。
【０２３１】
図２４に示すように、第４実施例では、ＣＣＭに関する後述する種々の演算を実行する演算装置４０と、データ入力を行なうキーボード，マウス等の入力機器４２と、演算の状態や後述する合否判定の結果を表示する表示機器４４と、この合否判定の結果や種々の演算式等を記憶する記憶装置４６と、ＣＣＭに必要なデータとしての分光反射率を取得するための分光光度計４８とを備える。そして、この演算装置４０では、以下の調合処理が行なわれる。
【０２３２】
まず、ステップＳ７１では、目標色を呈する釉（見本釉）を準備する。この場合の目標色は、調合済みの釉で呈される色であるので、この釉（見本釉）における顔料の調合割合は既知である。続くステップＳ７２では、この目標色見本釉を分光光度計４８で測色し、ＸＹＺ表色系での色評価値である三刺激値（実測値）ＸT ，ＹT ，ＺT を求める。この三刺激値は、分光光度計４８での測色により得られた目標色見本釉の分光反射率Ｒ'（λ）から、次の数式４１，４２に従って算出される。なお、算出された三刺激値は、表示機器４４に目標見本色とともに表示され、後述の処理に用いるために記憶装置４６に記憶される。また、以下に記す三刺激値等の演算結果は、その都度に、記憶装置４６に記憶される。
【０２３３】
【数４１】
2006030210.tif 000043
【数４２】
2006030210.tif 000044
【０２３４】
ここで、Ｓ（λ）は標準光の分光分布、ｘ（λ），ｙ（λ），ｚ（λ）（数式中ではバー付きである。以下同じ）は等色関数であり、いずれも既知の値である。なお、式中に記した（λ）は、分光反射率，分光分布，等色関数がいずれも波長λに依存していることを表わす。
【０２３５】
ステップＳ７３では、上記の見本釉が呈する目標色と近似する色を呈すると想定した調合割合でそれぞれの顔料を添加し、第１回目の試作釉を調合し、この試作釉についてステップＳ７２と同様にして測色する。これにより、この第１回目試作釉についての三刺激値（実測値）Ｘ1 ，Ｙ1 ，Ｚ1 を求める。この第１回目試作において、上記の調合割合での顔料の調合は技術者によってなされるものの、その値は任意性のある既知の値であり、従来のように試行錯誤して繰り返し調合する必要はなく、調合割合を定めるに当たって特段の経験や勘も要しない。この場合の第１回目試作釉についての三刺激値並びに顔料調合量（調合割合）も、表示機器４４に表示されると共に、記憶装置４６に記憶される。また、この際の顔料調合量は、入力機器４２から入力される。
【０２３６】
続くステップＳ７４では、ステップＳ７１で作製した目標色見本釉の呈する色とステップＳ７３における第１回目試作釉の呈する色との色差ΔＥ* (JIS Z 8730)が所定範囲内に納まるか否か、即ち色差ΔＥ* の合否判定を下す。この場合の色差ΔＥ* の許容値は、目標色見本の呈する色と第１回目試作釉の呈する色との違いが一見しては判別できない程度の値、例えば０．３〜０．５程度の値が予め入力機器４２から設定される。なお、色差ΔＥ* の許容値を０．３〜０．５以外の値とすることもできることは勿論である。
【０２３７】
このステップＳ７４で色差ΔＥ* が０．３〜０．５以内であると合格判定すれば、ステップＳ７３での第１回目試作の釉で目標色見本の色を再現できるので、それ以上の調合処理は不要であるとして総ての処理を終了する。つまり、第１回目試作の際の調合割合で調合した新たな釉は、目標色見本釉とほぼ同じ色を呈する。
【０２３８】
一方、ステップＳ７４で色差ΔＥ* が０．３〜０．５以内に納まらないとして不合格判定した場合には、続くステップＳ７５で、ステップＳ７２で求めた目標色見本の三刺激値ＸT ，ＹT ，ＺT とステップＳ７３で求めた第１回目試作釉の三刺激値Ｘ1 ，Ｙ1 ，Ｚ1 とから、次の数式４３に従って三刺激値の差ΔＸ，ΔＹ，ΔＺを算出する。この三刺激値の差ΔＸ，ΔＹ，ΔＺは、目標色見本の呈する色と第１回目試作釉の呈する色との色差ΔＥ* を反映した値となる。
【０２３９】
【数４３】
2006030210.tif 000045
【０２４０】
ステップＳ７６では、ステップＳ７３で行なった第１回目試作釉の調合時の各顔料（調合割合既知）に微量だけその顔料を追加調合した場合の三刺激値の変化率（微係数）を、ＣＣＭの手法を用いて、以下の手順に従って求める。なお、各顔料を追加調合する際の微量増加調合量も入力機器４２から入力される。
【０２４１】
まず、このステップＳ７６の最初のステップＳ９１では、図２６のフローチャートに示すように、第１回目試作釉の呈する色の三刺激値を第１回目試作時の各顔料の既知の調合割合からＣＣＭの手法で演算する。この場合の三刺激値（演算値）Ｘ1/E ，Ｙ1/E ，Ｚ1/E は、次の数式４４で表わされる。
【０２４２】
【数４４】
2006030210.tif 000046
【０２４３】
この数式４４にあっても、上記した数式４１，４２と同様、Ｓ（λ）は標準光の分光分布（既知）であり、ｘ（λ），ｙ（λ），ｚ（λ）は等色関数（既知）である。しかし、数式４４における分光反射率Ｒ（λ）は、三刺激値Ｘ1/E ，Ｙ1/E ，Ｚ1/E を顔料の特性から算出する都合上、やはり顔料の特性を用いて以下のようにして算出する。そして、この算出した分光反射率Ｒ（λ）を用いて、上記の数式４４から三刺激値Ｘ1/E ，Ｙ1/E ，Ｚ1/E を算出する。
【０２４４】
被着色物と着色剤の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）とは、以下の数式４５で表わされるダンカン（Duncan）の式と、数式４６で表わされるクベルカ−ムンク（Kubelka-Munk）の混色理論による式の関係にあり、これら数式に基づいて、任意の調合物の分光反射率Ｒ（λ）をＣＣＭにより求めることができる。
【０２４５】
【数４５】
2006030210.tif 000047
【０２４６】
ここで、ＫM ，ＳM は調合物（釉）の吸収係数と散乱係数、Ｋi ，Ｓi はｉ番目の成分（顔料）の吸収係数と散乱係数、Ｃi はｉ番目の成分（顔料）の調合割合である。
【０２４７】
【数４６】
2006030210.tif 000048
【０２４８】
この数式４６を変形すると、分光反射率Ｒ（λ）は、次の数式４７で与えられる。
【０２４９】
【数４７】
2006030210.tif 000049
【０２５０】
ここで、調合物の吸収係数と散乱係数の比（Ｋ／Ｓ）は、この数式４７で示されているように、各顔料の吸収係数Ｋi（λ）と散乱係数Ｓi（λ）と調合率Ｃi で規定される。従って、分光反射率Ｒ（λ）は、この比（Ｋ／Ｓ）から算出され、求めることができる。なお、ＫW ，ＳW は白色成分（白顔料）の吸収係数と散乱係数、ＣWは白色成分の調合割合である。
【０２５１】
続くステップＳ９２では、各顔料（顔料１，２，３）について、第１回目試作釉の調合割合をそれぞれ微量の調合量だけ個別に増加させた釉の呈する色についての三刺激値（演算値）Ｘ1/1/E ，Ｙ1/1/E ，Ｚ1/1/E を、第１回目試作時の各顔料の既知の調合割合並びに増加させた顔料の既知の調合割合からＣＣＭの手法で演算する。より詳しく説明すると、まず、第１回目試作釉に顔料１を微量の調合量（０．１×Ｃ1 ）だけ加え他の顔料２，３および白顔料は第１回目試作釉と同一調合率の釉の呈する色についての三刺激値を演算する。この場合にあっても、上記した数式４４〜数式４７が用いられ、数式４７におけるＣ1 は（Ｃ1 ＋０．１×Ｃ1 ）となり、Ｃ2 ，Ｃ3 ，ＣWはそのままである。そして、このように計算された三刺激値Ｘ1/1/E ，Ｙ1/1/E ，Ｚ1/1/E は、顔料１の微量増加調合に起因した三刺激値である。
【０２５２】
同様に、顔料２を微量の調合量（０．１×Ｃ2 ）だけ加え他の顔料は同一調合率の釉の呈する色についての三刺激値Ｘ1/2/E ，Ｙ1/2/E ，Ｚ1/2/E と、顔料３を微量の調合量（０．１×Ｃ3 ）だけ加え他の顔料は同一調合率の釉の呈する色についての三刺激値Ｘ1/3/E ，Ｙ1/3/E ，Ｚ1/3/E とを演算する。なお、白顔料についてもその三刺激値Ｘ1/W/E ，Ｙ1/W/E ，Ｚ1/W/E が同様に演算される。そして、各顔料（顔料１，２，３および白顔料）をそれぞれ微量増加して調合したことによる三刺激値の変化量と各顔料の微量増加量とから、次の数式４８に従って三刺激値の変化率（微係数）を算出する。その後は図２５のステップＳ７７に進む。
【０２５３】
【数４８】
2006030210.tif 000050
【０２５４】
ステップＳ７７では、目標色見本釉と第１回目試作釉の測色を経てステップＳ７５で求めた三刺激値の差ΔＸ，ΔＹ，ΔＺ（実測値の差）を補正するために必要な各顔料の追加増量調合量（顔料１はΔＣ1 ，顔料２はΔＣ2 ，顔料３はΔＣ3 ，白顔料はΔＣW）を、次の数式４９と数式５０を用いて算出する。
【０２５５】
【数４９】
2006030210.tif 000051
【数５０】
2006030210.tif 000052
【０２５６】
数式４９は、三刺激値の差ΔＸ，ΔＹ，ΔＺ（実測値の差）を補正するためのものであり、各顔料の追加増量調合量（ΔＣ1 ，ΔＣ2 ，ΔＣ3 ，ΔＣW）を変数とする。また、数式５０は、各顔料を変動させた場合の費用関数であり、各顔料の追加増量調合量（ΔＣ1 ，ΔＣ2 ，ΔＣ3 ，ΔＣW）をやはり変数とする。従って、数式４９，５０を４つの上記変数について解くことでこれら変数、即ち各顔料の追加増量調合量（ΔＣ1 ，ΔＣ2 ，ΔＣ3 ，ΔＣW）が定まる。なお、数式５０におけるＡi は、顔料i を単位量変動させるのに必要な費用である。
【０２５７】
この場合、ステップＳ７７では、ΔＣi がΔＣi ≧０を満たし、ΔＣWがΔＣW≧０を満たし、且つ、数式５０で表わされる費用関数Ｆを最小とするΔＣi 並びにΔＣWを線形計画法を用いて解く。これにより、ΔＣi 並びにΔＣWが負の値として算出されることはない。
【０２５８】
続くステップＳ７８では、各顔料（白顔料を含む）の調合割合を、求めた追加増量調合量を加味して算出し、追加増量調合した釉についての物性値を求め更新する。より具体的には、この調合割合と各顔料についての吸収係数および散乱係数等を用い、上記の数式４５〜数式４７に従って、分光反射率Ｒ（λ）を求める。そして、その後のステップＳ７９では、求めた分光反射率Ｒ（λ）を数式４４に代入して、各顔料を追加増量調合した釉の呈する色についての三刺激値（演算値）ＸEND/E，ＹEND/E，ＺEND/Eを演算する。
【０２５９】
ステップＳ８０では、この追加増量調合後の三刺激値ＸEND/E，ＹEND/E，ＺEND/Eと目標色見本釉についての三刺激値ＸT ，ＹT ，ＺT を用い(JIS Z 8730)、その色差ΔＥ* が所定範囲内に納まるか否かの合否判定を再度下す。このステップＳ８０で色差ΔＥ* が０．３〜０．５以内であると合格判定すれば、それまでのステップで求めた補正量で調合した釉であれば目標色見本の色を再現できることになる。よって、この場合は、調合処理が完了したとして処理を終了する。つまり、求めた各顔料についての補正量（追加増量調合量）を加味した調合割合が、調合済みの釉を再調合する際の最終的な調合割合として採用される。より具体的に説明すると、調合工程の変動等により色に変調を来たした調合済み釉は、その調合割合が既知であるので第４実施例における第１回目試作釉に相当するので、この色に変調を来たした調合済み釉に上記の補正量に従った配合割合で各顔料を追加して再調合すれば、目標色を呈する釉を再調合できることになる。
【０２６０】
一方、ステップＳ８０で色差ΔＥ* が０．３〜０．５以内に納まらないとして不合格判定した場合には、微係数計算，補正量計算のための三刺激値をそれまでの値からシフトするとともに、各顔料の微量増加量を変更する。具体的に説明すると、数式４３における三刺激値Ｘ1 ，Ｙ1 ，Ｚ1 に、ステップＳ７９で求めた三刺激値ＸEND/E，ＹEND/E，ＺEND/Eをシフトする。これより、このシフト後の三刺激値（ＸEND/E，ＹEND/E，ＺEND/E）から、ステップＳ７５では目標色見本との三刺激値の差（ΔＸ，ΔＹ，ΔＺ）が新たに算出される。
【０２６１】
また、各顔料（顔料１，２，３）について、ステップＳ７６で考慮した微量増加量を、それまでの調合量（この場合には、顔料１はＣ1 ＋０．１×Ｃ1 ）に０．１を乗じた微量増加量（（Ｃ1 ＋０．１×Ｃ1 ）×０．１）に変更する。これにより、各顔料は、ステップＳ８０で不合格とされた場合よりそれぞれ０．１だけ微量増量されるので、ステップＳ７６で、この微量増量を加味した調合割合から、各顔料が微量増量されたことによる三刺激値が上記した数式４４〜数式４７に従って求められる。そして、各顔料が微量増加されたことによるこの新たな三刺激値が数式４８における三刺激値Ｘ1/1/E ，Ｙ1/1/E ，Ｚ1/1/E 等に替わって用いられて改めて微係数が算出される。その後は既述したように各顔料の追加増量調合量（補正量）が求められ、ステップＳ８０で合格判定されるまで、上記の処理が繰り返される。
【０２６２】
以上説明したように第４実施例のコンピュータカラーマッチング方法では、それまで存在していた釉（目標色見本釉）に対して、その呈する色がある程度近似した色を呈するような調合割合で第１回目試作釉を調合し、その後、この第１回目試作釉に、目標色見本釉の呈する色と所定範囲で一致するような調合割合で各顔料を追加させた場合の追加増量調合量を求める。このため、顔料の除去を意味する負の調合量をＣＣＭにより求めることがない。よって、第４実施例のコンピュータカラーマッチング方法によれば、顔料が調合済みの釉（着色剤）の廃棄が不要となるために既存の釉を有効に利用することができる。また、既存の釉の利用に当たり、技術者が関与する工程はステップＳ７３における一回限りの試作釉の調合であり、この際に、技術者の長年の勘や経験を要しないので、釉の再調合を簡略化することができる。
【０２６３】
更に、第４実施例のコンピュータカラーマッチング方法では、釉を再調合するための各顔料の補正量の算出の際に、各顔料の補正量（ΔＣi 並びにΔＣW）を、数式５０で表わされる費用関数Ｆを用いた線形計画法の手法で求め、各顔料の補正に要する費用が最小となるようにした。このため、第４実施例のコンピュータカラーマッチング方法によれば、上記した既存の釉の有効利用と再調合に簡略化に加え、コスト低減をも図ることができる。
【０２６４】
次に、上記した第４実施例のコンピュータカラーマッチング方法を行なった際の各ステップでの処理により得られる三刺激値や微係数について、図２７を用いて説明する。
【０２６５】
図２７は、ステップＳ７２で取得した目標色見本釉についての三刺激値（色値）とステップＳ７３で取得した第１回目試作釉についての三刺激値の対比、並びに目標色見本釉，第１回目試作釉における各顔料の調合率（調合率）を表わす。そして、図中のΔＥ* は、目標色見本釉と第１回目試作釉との間の色差であり、この値に基づいてステップＳ７４での合否判定が下される。
【０２６６】
図２８は、図２７に掲げる第１回目試作釉に各顔料をそれぞれ微量だけずつ追加調合した場合の三刺激値の変化率（微係数）を示しており、上記の数式４８から演算される。この際、顔料１（赤顔料）の微量増量の場合における微係数算出には、演算した三刺激値Ｘ1/1/E ，Ｙ1/1/E ，Ｚ1/1/E が用いられ、顔料２（黄顔料）では三刺激値Ｘ1/2/E ，Ｙ1/2/E ，Ｚ1/2/E が、顔料３（青顔料）では三刺激値Ｘ1/3/E ，Ｙ1/3/E ，Ｚ1/3/E が、白顔料では三刺激値Ｘ1/W/E ，Ｙ1/W/E ，Ｚ1/W/E が用いられる。
【０２６７】
この図２９は、第４実施例によるコンピュータカラーマッチング方法による結果を示しており、各顔料について定めた最終的な調合割合で調合した釉と目標色見本釉とを、その調合率と三刺激値について対比して表わす。また、この両釉についての色差ΔＥ* は０．４７であり、ステップＳ８０で合格判定されたことが判る。
【０２６８】
以上の第４実施例に関しては，次のような変形も可能である。
【０２６９】
例えば、数式５０の費用関数Ｆに替わって、以下の数式５１で表わされる費用関数Ｆを用いてもよい。
【０２７０】
【数５１】
2006030210.tif 000053
【０２７１】
更には、ステップＳ７６にて顔料をそれぞれ微量ずつ追加増量した微係数を計算する際に、各顔料の追加増量調合量の最低値を最低追加増量調合量ΔＣstepとして予め規定しておき、このΔＣstepの整数倍ずつ各顔料を追加増量するよう構成することもできる。そして、このように各顔料を追加増量していき、最も目標色見本釉との色差ΔＥ* が小さくなるときの追加増量量を、微係数算出の際の補正量とし、この補正量を加味した調合割合を調合済みの釉を再調合する際の最終的な調合割合とすればよい。なお、この場合の最低追加増量調合量ΔＣstepは、各顔料を追加したときにその呈する色が僅かに変わる最小単位の追加増量量として規定される。また、この場合には、各顔料の追加増量量の最大許容量ΔＣmax を予め定めておき、各顔料の追加増量の総量をこのΔＣmax で規定すればよい。そして、この最大許容量ΔＣmax は、その顔料を追加したときにその呈する色が大きく変わり、他の顔料の追加増量では元の色への戻りができないと思われる追加増量として規定すればよい。
【０２７２】
また、上記の費用関数Ｆを用いた線形計画法に限られるわけではなく、他の手法を採ることもできる。例えば、上記の第４実施例で示したように、色差ΔＥ* で合否判定する場合には、その合否判定にある程度の幅が許容される。従って、この合否判定にある程度の曖昧さを導入したいわゆるファジィ線形計画法の手法を採ることもできる。このファジィ線形計画法の手法を採ったコンピュータカラーマッチング方法での結果を以下の図３０に示す。この図３０に示すように、各顔料について定めた最終的な調合割合で調合した釉と目標色見本釉とは、その色差ΔＥ* が０．２０であり、両者の色はよく一致していることが判る。
【０２７３】
また、上記の第４実施例では、陶器やタイルを色付けする釉を例に採り説明したが、繊維を染める染色剤にも適用できることは勿論である。
【０２７４】
この発明にかかるコンピュータカラーマッチング方法および装置は、陶器やタイルを色付けするための釉に混合される着色剤の調合割合の予測とその混合物（すなわち釉）の色予測に適用できる他に，繊維を染めるための染色剤等の各種の混合物における着色剤の調合割合の予測とその混合物の色予測等に適用可能である。
【図面の簡単な説明】
【０２７５】
【図１】第１実施例における処理の全体手順を示すフローチャート。
【図２】数式５で与えられる散乱係数ＳW'を示すグラフ。
【図３】ステップＳ１の詳細手順を示すフローチャート。
【図４】（ベース釉＋白色顔料）の物性値を求めるためのサンプルの調合率を示す説明図。
【図５】サンプルの分光反射率Ｒ（λ）を示す概念図。
【図６】ステップＳ１４の詳細手順を示すフローチャート。
【図７】顔料物性値決定用サンプルの調合率を示す説明図。
【図８】第１実施例で用いた検証用サンプルの調合率を示す説明図。
【図９】吸収係数ＫP（λ）の調合率ＣP に対する依存性の一例を示すグラフ。
【図１０】第１実施例における調合率の予測結果を示す説明図。
【図１１】第２実施例における処理の全体手順を示すフローチャート。
【図１２】ニューラルネットワークの構成を示す説明図。
【図１３】ステップＳ３１の詳細手順を示すフローチャート。
【図１４】コンピュータカラーマッチングで予測対象範囲ＰＡと複数のサンプルＭ１〜Ｍ７の三刺激値の分布を示す概念図。
【図１５】第２実施例における７つのサンプルＭ１〜Ｍ７の予測誤差を示す概念図。
【図１６】各サンプルの三刺激値の目標値（ステップＳ４１で決定された値）と予測誤差ΔMi（ΔＸ，ΔＹ，ΔＺ）を示す説明図。
【図１７】第２実施例におけるニューラルネットワークの学習の実証結果を示す説明図。
【図１８】第２実施例におけるコンピュータカラーマッチングの予測の実証結果を示す説明図。
【図１９】設計された色（標準色）と実際に製造される陶器の色のばらつきを示すｘ−ｙ色度図。
【図２０】第３実施例における処理の全体手順を示すフローチャート。
【図２１】ステップＳ５７におけるコンピュータカラーマッチングの詳細手順を示すフローチャート。
【図２２】第３実施例において用いた標準色見本の三刺激値と濃淡限度見本の三刺激値を示すテーブル。
【図２３】第３実施例による濃淡限度見本の三刺激値の予測結果と比較例と実調合率とを示すテーブル。
【図２４】各実施例のコンピュータカラーマッチング方法を実施するための装置を示すブロック図。
【図２５】第４実施例のコンピュータカラーマッチング方法における処理の全体手順を示すフローチャート。
【図２６】図２５のステップＳ７６の詳細処理を示すフローチャート。
【図２７】ステップＳ７２で取得した目標色見本釉についての三刺激値（色値）とステップＳ７３で取得した第１回目試作釉についての三刺激値の対比、並びに目標色見本釉，第１回目試作釉における各顔料の調合率（調合率）を表わすテーブル。
【図２８】図２７に掲げる第１回目試作釉に各顔料をそれぞれ微量だけずつ追加調合した場合の三刺激値の変化率（微係数）を示すテーブル。
【図２９】第４実施例によるコンピュータカラーマッチング方法による結果を示すテーブル。
【図３０】第４実施例による他の結果を示すテーブル。
【図１】
2006030210.tif 000054
【図２】
2006030210.tif 000055
【図３】
2006030210.tif 000056
【図４】
2006030210.tif 000057
【図５】
2006030210.tif 000058
【図６】
2006030210.tif 000059
【図７】
2006030210.tif 000060
【図８】
2006030210.tif 000061
【図９】
2006030210.tif 000062
【図１０】
2006030210.tif 000063
【図１１】
2006030210.tif 000064
【図１２】
2006030210.tif 000065
【図１３】
2006030210.tif 000066
【図１４】
2006030210.tif 000067
【図１５】
2006030210.tif 000068
【図１６】
2006030210.tif 000069
【図１７】
2006030210.tif 000070
【図１８】
2006030210.tif 000071
【図１９】
2006030210.tif 000072
【図２０】
2006030210.tif 000073
【図２１】
2006030210.tif 000074
【図２２】
2006030210.tif 000075
【図２３】
2006030210.tif 000076
【図２４】
2006030210.tif 000077
【図２５】
2006030210.tif 000078
【図２６】
2006030210.tif 000079
【図２７】
2006030210.tif 000080
【図２８】
2006030210.tif 000081
【図２９】
2006030210.tif 000082
【図３０】
2006030210.tif 000083
【手続補正書】
【提出日】平成17年8月30日(2005.8.30)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラーマッチング方法であって、
（ａ）前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める工程と、
（ｂ）既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、
（ｃ）前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める工程と、（ｄ）前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、
（ｅ）前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する工程と、を備える
コンピュータカラーマッチング方法。
【請求項２】
請求項１記載のコンピュータカラーマッチング方法であって、
前記工程（ｄ）は、前記１次調合物における前記着色剤の調合量に比べて微量の量の前記着色剤を前記１次調合物に増量補正したと仮定した場合について、前記色評価値の計算値の変化量を求める工程を含む、コンピュータカラーマッチング方法。
【請求項３】
請求項１又は請求項２記載のコンピュータカラーマッチング方法であって、
前記工程（ｅ）は、前記着色剤の増量に伴う派生費用を表わす費用関数を用いた線形計画法にて、前記着色剤のそれぞれについての最小の増量補正量を算出する工程を含む、コンピュータカラーマッチング方法。
【請求項４】
複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラーマッチング装置であって、
前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める手段と、
既知の調合割合で前記着色剤が調合された１次調合物について、前記所定の表色系での色評価値の実測値を求める手段と、
前記１次調合物についての前記既知の調合割合に基づいて、前記１次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める手段と、
前記１次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記１次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める手段と、
前記調合物見本と前記１次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する手段と、を備える
コンピュータカラーマッチング装置。