WO1997031247A1

WO1997031247A1 - Computer color matching method and apparatus

Info

Publication number: WO1997031247A1
Application number: PCT/JP1996/000738
Authority: WO
Inventors: Hiroshi Kumamoto; Hideharu Imoto; Hiroshi Tamae; Nobuo Adachi
Original assignee: Toto Ltd; Hiroshi Kumamoto; Hideharu Imoto; Hiroshi Tamae; Nobuo Adachi
Priority date: 1996-02-22
Filing date: 1996-03-21
Publication date: 1997-08-28
Also published as: CN1214240C; CN1207808A; JP3870421B2; JPH09229773A

Abstract

A scattering coefficient (SW') for a mixture of an object of coloring and a white pigment is expressed by a function of the white pigment content (CW). Computer color matching is executed for a verification sample, and a value which gives high prediction accuracy is determined as coefficients (a, b, c and d) (or SE) contained in the function representing the scattering coefficient (SW'). An absorption coefficient of (the object + white pigment), absorption coefficients of other coloring pigments and scattering coefficients are determined by using the scattering coefficient (SW') given by the function so determined. Computer color matching is executed by using physical constants of each component obtained in this way, and the color and the matching proportion of the mixture are predicted. In computer color matching, a prediction error can be reduced by utilizing a neural network.

Description

明細害コンピュータカラ一マッチング方法および装置 Specification damage Computer color matching method and apparatus

[技術分野] [Technical field]

本発明は、コンピュータカラーマッチングによつて着色剤の調合割合の予測または混合物の色予測を行なう方法およびそのための装置に関するものである。 The present invention relates to a method and apparatus for predicting the mixing ratio of a colorant or the color of a mixture by computer color matching.

[背景技術] [Background technology]

顔料や染料などの着色剤を被着色物に混合した混合物の色を予測するために、いわゆるコンピュータカラ一マツチングが利用されている。コンピュータカラーマッチングでは、被着色物と着色剤の吸収係数 Ki (λ ) と散乱係数 ) とを用い、ダンカン（Duncan) の式（数式 1 ) と、クベルカ一ムンク（ Kubelka-Munk) の混色理論による式（数式 2) に基づいて、住意の混合物の分光反射率 R (/) を求めることができる。 In order to predict the color of a mixture in which a coloring agent such as a pigment or a dye is mixed with a material to be colored, so-called computer color matching is used. In computer color matching, the absorption coefficient Ki (λ) and scattering coefficient) of the coloring matter and the colorant are used, and the Duncan equation (Equation 1) and the Kubelka-Munk color mixing theory are used. Based on the equation (Equation 2), the spectral reflectance R (/) of the mixture of resident can be obtained.

(la) (la)

i=] i =]

S\i： j j .-.(lb) S \ i: j j .-. (Lb)

i=l i = l

(s)_M~~ R~ *··⁽²⁾ ここで、 K_M , S_M は混合物の吸収係数と散乱係数、 Ki , S_± は i番目の成分の吸収係数と散乱係数、 Ci は i番目の成分の調合率である。但し、この明細書の数式においては、波長 Iに依存していることを示す「（）」は省略されている。混合物の成分は、被着色物と着色剤である。 (s) _M ~~ R ~ * ⁽²⁾ where K _M and S _M are the absorption and scattering coefficients of the mixture, Ki and S _± are the absorption and scattering coefficients of the i-th component, and Ci is This is the mixing ratio of the i-th component. However, in the mathematical expressions in this specification, “()” indicating that the wavelength depends on the wavelength I is omitted. The components of the mixture are an object to be colored and a colorant.

混合物の分光反射率 R (λ ) が解れば、その混合物の三刺激値 X， Υ, Ζが計算できるので、混合物の色を算出することができる。 Knowing the spectral reflectance R (λ) of a mixture, the tristimulus values X, Υ, and の of the mixture are calculated. Therefore, the color of the mixture can be calculated.

コンピュータカラ一マッチングを行なう場合には、被着色物と種々の着色剤の吸収係数 (λ ) と散乱係数 S_± (λ ) を予め求めておく必要がある。ところカ^ 吸収係数 Ki (； J と散乱係数 Si (λ ) を求めるのは必ずしも容易ではない場合が多い。 When computer color matching is performed, the absorption coefficient (λ) and scattering coefficient S _± (λ) of the material to be colored and various colorants must be determined in advance. However, it is often not always easy to obtain the absorption coefficient Ki (; J) and the scattering coefficient Si (λ).

被着色物と種々の着色剤の吸収係数 Ki (；. ) と散乱係数 Si ( . ) を求める方法としては、絶対法と相対法がある。絶対法は、各物質の吸収係数 Κ_± (λ ) と散乱係数 (；.) の絶対値を求める方法である。一方、相対法は、基準となる顔料（通常は白色顔料）の散乱係数 s_w を 1と仮定して、各物質の吸収係数There are an absolute method and a relative method for obtaining the absorption coefficient Ki (;.) And the scattering coefficient Si (.) Of the material to be colored and various colorants. The absolute method is a method for calculating the absolute values of the absorption coefficient Κ _± (λ) and the scattering coefficient (;.) Of each substance. On the other hand, the relative method, a primary pigment (usually white pigment) assuming scattering coefficient s _w 1 and the absorption coefficient of each substance

(λ ) と散乱係数 Si (A) の相対値を求める方法である。数式 2に示すように、混合物の分光反射率 R ）は混合物の吸収係数 K_M と散乱係数 S_M の比で与えられる。従って、各物質の吸収係数 Ki (； .) と散乱係数 Si (λ) の絶対値が不明でも、各物質の吸収係数 (λ ) と散乱係数 Si (λ ) の相対値が解れば、数式 1および 2から正しい分光反射率 R (Λ) を求めることができる。絶対法に用いるサンプルの作成はかなり困難であり、煩雑な作業を要するので、相対法が用いられるのが普通である。 This is a method of calculating the relative value of (λ) and the scattering coefficient Si (A). As shown in Equation 2, the spectral reflectance R) of the mixture is given by the ratio of the absorption coefficient K _M and the scattering coefficient S _M of the mixture. Therefore, even if the absolute values of the absorption coefficient Ki (;.) And the scattering coefficient Si (λ) of each substance are unknown, if the relative values of the absorption coefficient (λ) and the scattering coefficient Si (λ) of each substance are known, From Equations 1 and 2, the correct spectral reflectance R (Λ) can be determined. Since the preparation of samples for the absolute method is quite difficult and requires complicated work, the relative method is usually used.

従来の相対法によって各物質の吸収係数 Ki (λ ) と散乱係数 Si (λ ) とを求める際には、以下のようにして、基準となる白色顔料の吸収係数 K_w (Λ) と散乱係数 S_w を決定する。まず、被着色物に白色顔料のみを混合した混合物を作成して、その分光反射率 R (λ ) を測定する。従来の相対法では、被着色物を無色透明とみなして、被着色物の吸収係数と散乱係数はともに 0であると仮定している。従って、被着色物と白色顔料のみを混合した混合物に関しては、上記の数式 1で与えられる吸収係数 Κ_Μ (λ ) と散乱係数 S„ (λ ) は、白色顔料の散乱係数 S_w と吸収係数 K_w (λ ) にそれぞれ等しい。また、上記の数式 2において S„ =S_W = 1と仮定すると、分光反射率 R (/) の測定値から白色顔料の吸収係数 K_w (λ ) を求めることができる。 When determining the absorption coefficient Ki (λ) and the scattering coefficient Si (λ) of each substance by the conventional relative method, the absorption coefficient K _w (Λ) of the reference white pigment is calculated as follows. to determine the scattering coefficient S _w. First, a mixture is prepared by mixing only a white pigment with an object to be colored, and its spectral reflectance R (λ) is measured. In the conventional relative method, the object to be colored is regarded as colorless and transparent, and it is assumed that both the absorption coefficient and the scattering coefficient of the object to be colored are zero. Thus, for a mixture obtained by mixing only the color bodies and the white pigment, the absorption coefficient given by Equation 1 above kappa _Micromax (lambda) and the scattering coefficient S "(lambda) is the scattering coefficient of the white pigment S _w and absorption coefficient K _w (λ) equal to, respectively. also, obtaining assuming in equation 2 above S "= S _W = 1 and the absorption coefficient of the white pigment from the measurements of the spectral reflectance R (/) K _w a (lambda) be able to.

上述のように、従来の相対法では、被着色物を無色透明とみなして、その吸収係数と散乱係数を 0と仮定していた。ところが、現実には陶磁器等をはじめ，被着色物（陶磁器の場合は釉薬層）が無色透明でない場合が多い。白色顔料の調合率が多い場合には、被着色物が無色透明であるとみなすことによる誤差は小さいが、白色顔料の調合率が少ない場合にはその誤差が無視できない程度に大きくなる。このような誤差を回避するためには、上記数式 1において、無色透明でない被着色物の吸収係数や散乱係数を考慮にいれる必要がある。しかし、被着色物の吸収係数や散乱係数を求めるには、被着色物単体の薄板を作成するなどの煩雑な作業を要する。特に、陶器の釉のように、単体で薄板を作成することが困難な被着色の場合には、無色透明でない被着色物の吸収係数と散乱係数を求めることは困難であった。 As described above, in the conventional relative method, an object to be colored is regarded as colorless and transparent, and its absorption coefficient and scattering coefficient are assumed to be zero. However, in reality, in many cases, the objects to be colored (glaze layers in the case of ceramics), such as ceramics, are not colorless and transparent. When the blending ratio of the white pigment is large, the error caused by assuming that the object to be colored is colorless and transparent is small, but when the blending ratio of the white pigment is small, the error becomes so large that it cannot be ignored. In order to avoid such errors, it is necessary to take into account the absorption coefficient and scattering coefficient of the colorless and non-transparent object in Equation 1 above. However, in order to obtain the absorption coefficient and the scattering coefficient of an object to be colored, a complicated operation such as creating a thin plate of the object to be colored is required. In particular, in the case of coloring, such as ceramic glaze, in which it is difficult to make a thin plate by itself, it was difficult to determine the absorption coefficient and scattering coefficient of a coloring object that is not colorless and transparent.

この発明の第 1の目的は、無色透明でない被着色物の吸収係数や散乱係数を考慮してコンピュータカラーマッチングを行なうことにある。 A first object of the present invention is to perform computer color matching in consideration of an absorption coefficient and a scattering coefficient of an object to be colored that is not colorless and transparent.

この発明の第 2の目的は、コンピュータカラ一マッチングにおける予測誤差を減少させることにある。 A second object of the present invention is to reduce a prediction error in computer color matching.

この発明の第 3の目的は、コンピュータカラ一マッチングを利用して，調合済み着色剤の有効利用を図りつつ、その再調合を簡略化することにある。 A third object of the present invention is to make effective use of the prepared colorant by using computer color matching and to simplify re-formulation.

[発明の開示] [Disclosure of the Invention]

この発明は，コンピュータカラ一マッチングによつて着色剤の調合割合の予測または混合物の色予測を行なう方法に関する。このコンピュータカラ一マツチング方法は， The present invention relates to a method for predicting a mixing ratio of a colorant or a color of a mixture by computer color matching. This computer color matching method is

無色透明でない被着色物に白色着色剤を混合した第 1の混合物の散乱係数 S _w_ ' を、前記白色着色剤の調合率 C_w に依存した関数 f ( C_K ) として準備するェ程と、 Preparing a scattering coefficient S _w ′ of a first mixture obtained by mixing a white colorant with a colorless and transparent material as a function f (C _K ) depending on the blending ratio C _w of the white colorant; ,

前記散乱係数 S _w' を基準として、前記第〗の混合物の吸収係数 K _w，を前記調合率 C_w に依存する形式で求める工程と、 Obtaining an absorption coefficient K _w , of the second mixture, based on the scattering coefficient S _w ′ in a form depending on the mixing ratio C _w ;

前記散乱係数 S _w' を基準として、白色でない有色着色剤の吸収係数 K _p と散乱係数 S _p とを前記有色着色剤の調合率 C_p に依存する形式で求める工程と、所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数 K_w' , _p および前記散乱係数 S_w' , S_p を用いたコンピュータカラーマッチングを行なうことによって求める工程と、 The basis of the scattering coefficient S _w ', a step of determining a format which depends the turbulent coefficient dispersion and absorption coefficient K _p of the colored colorant not white S _p in Formulation ratio C _p of the colored colorant, a desired color Mixing ratio of the colorant to adjust the mixture having The color of the mixture produced in formulation ratio, the absorption coefficient K _w ', _p, and the scattering coefficient S _w', a step of determining by performing computer color matching using S _p,

を備える。 Is provided.

上記方法では、無色透明でない被着色物に白色着色剤を混合した第】の混合物の散乱係数 s_w，を関数 f (C_w ) で表現し、この散乱係数 s_w，を基準として用いるので、散乱係数 s_w' に無色透明でない被着色物の影響を舍めることができる。 In the above method, the scattering coefficient s _w , of the second mixture obtained by mixing the white colorant with the non-colorless and transparent coloring object is expressed by a function f (C _w ), and this scattering coefficient s _w , is used as a reference. In addition, the influence of the non-colorless and transparent object on the scattering coefficient s _w 'can be suppressed.

好ましい実施例においては、 In a preferred embodiment,

前記第 1の混合物の散乱係数 s_w' を準備する工程は、 Preparing a scattering coefficient s _w ′ of the first mixture,

(a) 前記被着色物に前記白色着色剤を混合して、前記白色着色剤の調合率 C_w が異なる複数個の第 1の混合物を作成するとともに、前記複数個の第 1の混合物の分光反射率をそれぞれ測定する工程と、 (a) said mixing the white colorant to be colored object, the spectral of the with Formulation ratio C _w of the white colorant to create a plurality of different first mixture, the first mixture of said plurality Measuring the reflectivity,

( b ) 前記被着色物に前記有色着色物を混合した第 2の混合物を作成するとともに、前記第 2の混合物の分光反射率を測定する工程と、 (b) a step of preparing a second mixture in which the colored object is mixed with the colored object, and measuring a spectral reflectance of the second mixture;

(c) 前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記有色着色剤の調合率 C_p が異なる複数個の第 3の混合物を作成するとともに、前記複数個の第 3の混合物の分光反射率をそれぞれ測定する工程と、 (c) said mixing with the white colorant to be colored object and the colored colorant, along with Formulation ratio C _p of the colored colorant to create a plurality of different third mixture, several said double Measuring the spectral reflectances of the third mixture of

(d) 前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記第 3 の混合物とは調合割合が異なる第 4の混合物を作成するとともに、前記第 4の混合物の分光反射率を測定する工程と、 (d) mixing the white colorant and the colored colorant with the object to be colored to form a fourth mixture having a different mixing ratio from the third mixture, and the fourth mixture. Measuring the spectral reflectance of the

(e) 前記複数個の第】の混合物の散乱係数 S„' を前記白色着色剤の調合率 C_w の関数 f (C_w ) によって表わすとともに、前記関数 f (C_w ) に舍まれる係数の値を仮決定する工程と、 (e) The scattering coefficient S „′ of the plurality of the second mixtures is represented by a function f (C _w ) of the mixing ratio C _w of the white colorant, and the coefficient f (C _w ) Temporarily determining the value of

(f ) 前記複数個の第 1の混合物の分光反射率の測定値と前記関数 f (C_w ) とを用いて、前記調合率 C_w に依存した形式で前記第 1の混合物の吸収係数 K_w' を求める工程と、 (g) 前記第 2の混合物の分光反射率と、前記複数個の第 3の混合物の分光反射率と、前記関数 f (C_w ) と、前記吸収係数 K_w' とを用いて、前記調合率 C_p に依存した形式で前記有色着色剤の吸収係数と散乱係数とをそれぞれ求める工程と、 (f) Using the measured values of the spectral reflectances of the plurality of first mixtures and the function f (C _w ), the absorption coefficient K of the first mixture in a form dependent on the blending rate C _w _w ' (g) using the spectral reflectance of the second mixture, the spectral reflectances of the plurality of third mixtures, the function f (C _w ), and the absorption coefficient K _w ′, Determining an absorption coefficient and a scattering coefficient of the colored colorant in a format dependent on the rate C _p ;

(h) 前記吸収係数 K_w，， K_P および前記散乱係数 S_w，， S_P を用いて、前記第 4の混合物に関するコンピュータカラーマッチングを行なうとともに、前記コンピュータカラーマッチングによる予測精度を向上させるように前記関数 f (C _w ) に舍まれる係数を修正する工程と、 (h) using the absorption coefficients K _w , K _p and the scattering coefficients S _w , S _p to perform computer color matching on the fourth mixture and improve prediction accuracy by the computer color matching Correcting the coefficient included in the function f (C _w ),

( i ) 前記工程（f) ないし（h) を繰り返すことによって前記関数 ί (C_w ) に舍まれる係数を決定する 1:程と、 (i) determining the coefficients involved in the function ί (C _w ) by repeating the steps (f) to (h);

を備える。 Is provided.

こうして関数 f (C_w ) の係数を決定するようにすれば、コンピュータカラーマッチングによる予測精度が高くなるように関数 f (C_w ) の係数を求めることができる。 By determining the coefficients of the function f (C _w ) in this way, the coefficients of the function f (C _w ) can be obtained so that the prediction accuracy by computer color matching is improved.

前記関数 f (c_w ) が定数 _{a i} と係数 S_B とを含む次の式で与えられるようにすることが好ましい。： Wherein the function f (c _w) it is preferable to be given by the following equation including the constants _ai and coefficient S _B. :

f (C_w ) = (C_K + a, · S_B ) Z (C_w + _3l ) 。 f (C _w ) = (C _K + a, · S _B ) Z (C _w + _3l ).

上記の関数 f (C_H ) の形式は、白色着色剤単独の散乱係数を 1と仮定した場合に理論的に得られるものなので、散乱係数 S_w' を表わす関数として適切なものである。また、 1つの係数 S_B を舍むだけの簡単な形式を有しているので、係数 S_B を決定するのが容易である。 The form of the above function f (C _H ) is theoretically obtained assuming that the scattering coefficient of the white colorant alone is 1, and is therefore appropriate as a function representing the scattering coefficient S _w '. Also, since one coefficient S _B have a simple form only舍Mu, it is easy to determine the engagement number S _B.

この発明の他の局面によれば，コンピュータカラーマッチング方法は、（a) 複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備するェ稈と、（b) 前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、（c) 前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する工程と、 (d) 前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、（e ) 前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラ一マッチングにより行なう工程と、を備える。 According to another aspect of the present invention, a computer color matching method includes: (a) mixing a plurality of colorants to prepare a plurality of samples having different mixing ratios; and (b) the plurality of samples. Measuring the spectral reflectances of the plurality of samples, and obtaining actual measured values of coordinate values of a predetermined color system representing respective colors of the plurality of samples from the measured values of the spectral reflectances; Calculating, for each of the samples, a prediction error of the coordinate values of the color system; and Analyzing the relationship between the coordinate values of the color system with respect to the sample and the prediction error by a predetermined error correction method; and (e) calculating a target value or a prediction value of computer color matching using the error correction method. Predicting the blending ratio of the colorant in the new mixture or predicting the color of the mixture by computer color matching while making corrections.

複数のサンプルに関する所定の表色系の座標値とその予測誤差とを所定の誤差補正法で分析し、その誤差補正法によってコンピュータカラーマッチングの目標値または予測値を補正しつつ予測を行なうので、各成分の吸収係数 Ki と散乱係数 S _Α とを補正することなく予測誤差を減少させることができる。 Since the coordinate values of the predetermined color system and the prediction errors thereof for a plurality of samples are analyzed by a predetermined error correction method, and the prediction is performed while correcting the target value or the predicted value of the computer color matching by the error correction method. it is possible to reduce the prediction error without correcting the absorption coefficient Ki and scattering coefficient number S _Alpha of each component.

好ましい実施例においては、前記工程（d ) は、前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を二ユーラルネットワークに学習させる工程、を含み、前記工程（e ) は、学習済みのニューラルネットワークを用いてコンピュータカラ一マッチングによる予測を行なう工程を舍む。 In a preferred embodiment, the step (d) includes a step of causing a dual neural network to learn a relationship between the coordinate values of the color system and the prediction error for the plurality of samples, and the step (e). ) Shows a process of performing prediction by computer color matching using a trained neural network.

ニューラルネットワークを用いて予測誤差を補正するようにすれば、数多くのサンプルを学習させることによって、予測誤差を減少させることができる。また、前記ニューラルネットワークは、 3つのニューロンで構成される入力層と、複数のニューロンを舍む中間層と、 3つのニューロンで構成される出力層と、で構成される三層の階餍構造を有することが好ましい。 If the prediction error is corrected using a neural network, the prediction error can be reduced by learning a large number of samples. The neural network has a three-layer hierarchical structure including an input layer including three neurons, an intermediate layer including a plurality of neurons, and an output layer including three neurons. It is preferred to have.

この発明のさらに他の局面によれば，コンピュータカラーマッチング方法は、所望の色を有する目標混合物の着色剤の調合割合を予測する方法であって、 According to yet another aspect of the present invention, a computer color matching method is a method for predicting a mixing ratio of a colorant in a target mixture having a desired color,

( a ) 着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める工程と、 (a) obtaining an actual measurement value of coordinates of a predetermined color system for an adjacent color sample having a colorant formulation known and having a color close to the desired color;

( b ) 前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算使から計算誤差を求める工程と、 (b) obtaining a calculated value of the coordinates of the color system representing the color of the adjacent color sample from the known blending ratio of the adjacent color sample, and obtaining a calculation error from the actual measurement value and the calculator.

( c ) ¾UH目標混合物の色に対する ii己表色系の座標の目標値を設定する工程と、 (c) setting a target value of ii self-color system coordinates for the color of the ¾UH target mixture;

( d ) 前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマツチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する工程と、 (d) correcting the target value by using the calculation error, and executing computer color matching using the corrected target value, whereby the colorant of the target mixture is Estimating the blending ratio;

を備える。 Is provided.

近接色サンプルの表色系の座標値の計算誤差によって目標混合物の表色系の座標の目標彼を補正し、補正後の目標値を用いてコンピュータカラーマッチングを行なうので、表色系の座標値の計算誤差を低減した条件下でコンピュータカラーマッチングを実行できる。 The target colorimetric coordinate of the target mixture is corrected by the calculation error of the coordinate value of the colorimetric system of the adjacent color sample, and computer color matching is performed using the corrected target value. Computer color matching can be performed under conditions that reduce value calculation errors.

好ましい実施例によれば、前記工程（a ) は、複数のサンプルに関して、着色剤の調合割合と、前記表色系の座標の実測値とを含むデータベースから、前記目標混合物との色差が最小となるサンプルを選択することによつて前記近接色サンプルを検索する工程を舍む。 According to a preferred embodiment, in the step (a), for a plurality of samples, a color difference from the target mixture is minimized from a database including a mixing ratio of a coloring agent and an actual measurement value of coordinates of the color system. The step of searching for the adjacent color sample by selecting the sample to be performed is provided.

データベースから色差が最小のサンプルを近接色サンプルとして検索するようにすれば、任意の色を有する目標混合物に関して、コンピュータカラーマツチングの予測誤差を低減できる。 By retrieving the sample having the smallest color difference from the database as the closest color sample, it is possible to reduce the prediction error of computer color matching for a target mixture having an arbitrary color.

この発明のさらに他の局面によれば，コンピュータカラ一マッチング方法は、複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求める方法であって、 According to yet another aspect of the present invention, a computer color matching method includes: adjusting a compounding ratio of a plurality of colorants so that the mixture obtained by mixing the plurality of colorants exhibits a color similar to a desired target color. Is the method of seeking

( a ) 前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める工程と、 (a) obtaining a measured value of a color evaluation value in a predetermined color system for a preparation sample exhibiting the target color;

( b ) 既知の調合割合で前記着色剤が調合された〗次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、 (b) a step of obtaining an actual measurement value of a color evaluation value in the predetermined color system, for a secondary preparation in which the colorant is prepared at a known preparation ratio;

( c ) 前記 1次調合物についての前記既知の調合割合に基づいて、前記 1次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める工程と、 (c) obtaining a calculated value of a color evaluation value of the color represented by the primary formulation in the predetermined color system, based on the known blending ratio of the primary formulation;

( d ) 前記】次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記 1次調合物から前記着色剤增量調合物への前記色評価値の計算値の変化量を求める工程と、(d) Calculating the color evaluation value in the predetermined color system for the colorant-enhanced formulation assuming that the colorant in the next formulation has been corrected for the increased amount of the colorant, and obtaining the calculated value from the primary formulation Determining the amount of change in the calculated value of the color evaluation value to the colorant amount formulation,

( e ) 前記調合物見本と前記 1次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する丁-程と、を備える。 (e) the colorant based on the amount of change in the calculated value of the color evaluation value so that the difference between the measured value of the color evaluation value of the sample of the formulation and the primary formulation is equal within a predetermined range. And a step for calculating each amount of increase correction.

ヒ記のコンピュータカラーマッチング方法では、既知の調合割合の 1次調合物に着色剤を仮に増量補正した場合、その際の色評価値の計算値の変化量を求める。そして、この色評価値の計算ィ ¾の変化量に基づいてそれぞれの着色剤についての増量補正量を算出し、調合物見本と 1次調合物との色評価値の実測値の差を所定範囲で一致させる。このため、求めた增暈補正量だけそれぞれの着色剤を 1次調合物に実際に追加調合すれば、調合物見本の目標色若しくはこれに近似した色を呈する調合物に、 1次調合物を再調合できる。従って、本発明のコンピュータ力ラーマッチング方法によれば、着色剤の除去を意味する負の補正量を求めることがないので、着色剤が調合済みの調合物（ 1次調合物）の廃棄が不要となり、既存の調合物を有効に利用することができる。また、この際に 1次調合物を既存の調合物ではなく新たに試験的に調合するにしても、その調合に技術者が一回限り関与すればよく、しかもその際に、技術者の長年の勘や経験を要しないので、調合物の再調合を簡略化することができる。 According to the computer color matching method described in Hi, when the colorant is temporarily increased and corrected in a primary formulation having a known blending ratio, the amount of change in the calculated color evaluation value at that time is determined. Then, based on the amount of change of the color evaluation value calculation amount, the amount of increase correction for each colorant is calculated, and the difference between the actually measured value of the color evaluation value of the preparation sample and that of the primary preparation is determined within a predetermined range. To match. For this reason, if each colorant is actually added to the primary formulation by the calculated halo correction amount, the primary formulation will be added to the target color of the sample of the formulation or a formulation exhibiting a color similar to the target color. Can be reconstituted. Therefore, according to the computer color matching method of the present invention, there is no need to obtain a negative correction amount meaning removal of the colorant, so that it is not necessary to dispose of the preparation in which the colorant has been prepared (primary preparation). The existing formulation can be used effectively. In this case, even if the primary formulation is newly prepared on a trial basis instead of the existing formulation, the technician need only be involved once in the formulation, and at that time, Since no intuition or experience is required, re-mixing of the composition can be simplified.

好ましい実施例においては、 In a preferred embodiment,

前記工程（d ) は、前記〗次調合物における前記着色剤の調合量に比べて微量の量の前記着色剤を前記 1次調合物に増量補正したと仮定した場合について、前記色評価値の計算値の変化量を求める工程を含む。 In the step (d), the color evaluation value is calculated based on the assumption that a small amount of the coloring agent is corrected to be added to the primary formulation compared to the blending amount of the coloring agent in the secondary formulation. Calculating the amount of change in the calculated value.

このコンピュータカラーマツチング方法では、色評価値の計算値の変化量を求めるために行なう着色剤の仮の増量補正を微量単位で行なって、微小単位の変化量を求めることができる。よって、それぞれの着色剤の増量補正量を高い精度で求めることを通して、目標色によく一致した色を呈することができる調合物を得ることができる。 In this computer color matching method, a temporary increase correction of the colorant, which is performed in order to obtain a change amount of the calculated value of the color evaluation value, is performed in a minute unit, and a change amount in a minute unit can be obtained. Thus, by obtaining the amount of correction for the amount of each colorant with high accuracy, it is possible to obtain a composition that can exhibit a color that is in good agreement with the target color.

また、前記工程（e ) は、前記着色剤の増量に伴う派生費用を表わす費用関数を用いた線形計 ί面法にて、前記着色剤のそれぞれについての最小の増量補正量を算出する工程を舍むことが好ましい。 Further, the step (e) includes a step of calculating the minimum amount of increase correction for each of the colorants by a linear measurement surface method using a cost function representing a derivation cost associated with the increase in the amount of the colorant. It is preferable to settle.

このコンピュータカラーマツチング方法では、それぞれの着色剤を実際に追加増量する際の最小の増量補正量を算出できるので、着色剤の使用量の低減を通して派生費用の最小化を図ることができ、コストを低減することができる。 In this computer color matching method, since the minimum amount of increase correction when each colorant is actually added is calculated, the amount of colorant used can be reduced. As a result, derivation costs can be minimized and costs can be reduced.

この発明は，また，コンピュータカラ一マッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置に向けられている。この装置は、無色透明でない被着色物に白色着色剤を混合した第〗の混合物の散乱係数 S _w_The present invention is also directed to an apparatus for predicting a mixing ratio of a colorant or a color of a mixture by computer color matching. This device has a scattering coefficient S _w _ of a second mixture in which a white colorant is mixed with a colorless and transparent object to be colored.

' を、前記白色着色剤の調合率 C_w に依存した関数 f ( C_w ) として作成する手段と、 'As a function f (C _w ) depending on the blending ratio C _w of the white colorant,

前記散乱係数 s _w，を基準として、前記第〗の混合物の吸収係数 κ„，を前記調合率 c_w に依存する形式で求める手段と、 Means for determining an absorption coefficient κ „, of the second mixture in a form depending on the mixing rate c _w , based on the scattering coefficient s _w ,

前記散乱係数 s _w' を基準として、白色でない有色着色剤の吸収係数 K _P と散乱係数 s _P とを前記有色着色剤の調合率 c_p に依存する形式で求める孚段と、所望の色を有する混合物を調整するための若色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数 K_w' ， K _P および前記散乱係数 s _w' ， S _P を用いたコンピュータカラーマッチングを行なうことによって求める手段と、 With reference to the scattering coefficient s _w ′, the absorption coefficient K _P and the scattering coefficient s _P of the non-white colored colorant are determined in a form depending on the mixing ratio c _p of the colored colorant, and a desired color is obtained. The blending ratio of a young agent for preparing a mixture having the following formulas, or the color of the mixture produced at a predetermined blending ratio, is calculated using the absorption coefficients K _w ′, K _P and the scattering coefficients s _w ′, S _P Means to obtain by performing the computer color matching used;

を備える。 Is provided.

この発明の他の局面によれば，コンピュータカラ一マツチング装置は、複数の着色剤を混合することによつて作成された調合率が互いに異なる複数のサンプルの分光反射率をそれぞれ測定する手段と， According to another aspect of the present invention, a computer color matching device comprises: means for measuring spectral reflectances of a plurality of samples having different mixing ratios formed by mixing a plurality of colorants;

前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める手段と、 Means for obtaining, from the measured values of the spectral reflectance, actual measured values of coordinate values of a predetermined color system representing respective colors of the plurality of samples;

前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出する手段と、 Means for calculating a prediction error of the coordinate value of the color system for each of the plurality of samples;

前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する手段と、 Means for analyzing the relationship between the coordinate values of the color system for the plurality of samples and the prediction error, using a predetermined error correction method,

前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう手段と、を備える。 Means for correcting the target value or the predicted value of the computer color matching using the error correction method, and performing the prediction of the mixing ratio of the colorant in the new mixture or the color prediction of the mixture by computer color matching, Is provided.

この発明のさらに他の局面によれば，コンピュータカラ一マッチング装置は、着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める手段と、 According to still another aspect of the present invention, a computer color matching apparatus includes: an actual measurement value of coordinates of a predetermined color system for a proximity color sample having a colorant mixing ratio known and having a color close to the desired color; Means to determine

前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、 ii記実測値と前記計算値から計算誤差を求める丰段と、 Calculating a calculated value of the coordinates of the color system representing the color of the adjacent color sample from the known blending ratio of the adjacent color sample, and ii obtaining a calculation error from the actually measured value and the calculated value;

前記目標混合物の色に対する前記表色系の座標の目標値を設定する手段と、前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する手段と、 Means for setting a target value of the coordinates of the color system for the color of the target mixture; correcting the target value using the calculation error; and executing computer color matching using the corrected target value. Means for predicting the blending ratio of the colorant in the target mixture,

を備える。 Is provided.

この発明の他の局面によれば，コンピュータカラーマッチング装置は，複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるものである。このコンピュータカラーマツチング装置は、 According to another aspect of the present invention, a computer color matching device calculates a blending ratio of a plurality of coloring agents so that the formulation obtained by blending the plurality of coloring agents exhibits a color similar to a desired target color. Things. This computer color matching device

前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める手段と、 Means for obtaining an actual measurement value of a color evaluation value in a predetermined color system for a preparation sample exhibiting the target color;

既知の調合割合で前記着色剤が調合された】次調合物について、前記所定の表色系での色評価値の実測俯を求める手段と、 Means for obtaining the actually measured depression of the color evaluation value in the predetermined color system for the next preparation, wherein the colorant is prepared at a known preparation ratio;

前記〗次調合物についての前記既知の調合割合に基づいて、前記 1次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める丰段と、 A step of obtaining a calculated value of a color evaluation value of the color represented by the primary formulation in the predetermined color system based on the known blending ratio of the primary formulation;

前記 1次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記 1次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める手段と、 Calculating the color evaluation value in the predetermined color system for the colorant-enhanced formulation assuming that the colorant has been increased in the primary formulation and calculating the coloration from the primary formulation Means for determining the amount of change in the calculated value of the color evaluation value to the agent extension formulation,

前記調合物見本と前記 1次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの増量補正量を算出する手段と、を傰える。 Each of the colorants based on the amount of change in the calculated value of the color evaluation value so that the difference between the measured value of the color evaluation value of the preparation sample and that of the primary preparation matches within a predetermined range. Means for calculating the increase correction amount.

この発明は，さらに，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法に向けられている。この方法は、 The invention is further directed to a method for determining the absorption and scattering coefficients of a colorant used in computer color matching. This method

無色透明でない被着色物に白色着色剤を混合した第 1の混合物の散乱係数 s_w_ ' と吸収係数 κ„' を決定する工程と、 Determining a scattering coefficient s _w _ 'and an absorption coefficient κ „' of a first mixture obtained by mixing a white colorant with a colorless and transparent object;

前記散乱係数 s_w' と吸収係数 K_w，の少なくとも一方を基準として、白色でない有色着色剤の吸収係数 κ_ρ と散乱係数 Sp の相対値をそれぞれ決定する工程と、を備える。 Determining the relative values of the absorption coefficient κ _ρ and the scattering coefficient Sp of the non-white colored colorant on the basis of at least one of the scattering coefficient s _w ′ and the absorption coefficient K _w , respectively.

好ましい実施例によれば，前記第〗の混合物の散乱係数 s_w' 力前記白色着色剤の調合率 c_wと定数 _{3 l} と係数 S_B とを舍む次の関係で与えられる： According to a preferred embodiment, the scattering coefficient s _w ′ of the second mixture is given by the following relationship involving the mixing ratio c _w of the white colorant, a constant _{3 l} and a coefficient S _B :

S_w' =(C_W +a, -S_B )/(C„ +_a】 ). S _w '= (C _W + a, -S _B ) / (C „+ _a ]).

他の好ましい実施例によれば，前記第〗の混合物の散乱係数 S_w' 力，前記白色着色剤の調合率 C_wと定数 a， b， d， e ,C_WUを舍む次の関係で与えられる：According to another preferred embodiment, the scattering coefficient S _w 'of the second mixture, the mixing ratio C _w of the white colorant and the following relations a, a, b, d, e, C _WU Given by:

S_w， - a-C_w + b (C_w。 ≤ C_w); S _w ,-aC _w + b (C _w . ≤ C _w );

S_w， = d-C_w + e (C_w < C_w0). S _w , = dC _w + e (C _w <C _w0 ).

この発明は，また，コンピュータカラーマッチングを利用して作成された混合釉に向けられている。この混合釉は， The present invention is also directed to mixed glazes made using computer color matching. This mixed glaze

無色透明でないベース釉に白色着色剤を混合した第】の混合釉の散乱係数 S_w_ ' と吸収係数 K_w' を決定し、 Determine the scattering coefficient S _w _ 'and the absorption coefficient K _w ' of the mixed glaze in which a white colorant is mixed with a colorless base glaze,

前記散乱係数 S_w，と吸収係数 K_w' の少なくとも一方を基準として、白色でない有色着色剤の吸収係数 K_p と散乱係数 S_p の相対値をそれぞれ決定し、前記吸収係数 K_w，， K_P および前記散乱係数 s_w，， S_P を用いたコンビユータカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、 Based on at least one of the scattering coefficient S _w , and the absorption coefficient K _w ′, the relative values of the absorption coefficient K _p and the scattering coefficient S _p of the non-white colored colorant are determined, and the absorption coefficient K _w , , K _P and the scattering coefficient s _w,, by performing Konbiyu Takara over matching using S _P, formulation ratio of the colorant to adjust the mixed glaze having a desired color, or, at a predetermined compounding ratio Predict the color of the resulting mixed glaze,

前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，ことによって生成される。 Mixing the base glaze, the white colorant, and the colored colorant to have the predicted mix proportion or the predicted color; Generated by

この発明は，さらに，コンピュータカラーマッチングを利用して作成された混合釉を用いて製造された陶磁器にも向けられている。この陶磁器は、 The present invention is further directed to ceramics manufactured using a mixed glaze made using computer color matching. This ceramic is

無色透明でないベ一ス釉に白色着色剤を混合した第 1の混合釉の散乱係数 S _w— ' と吸収係数 K_w' を決定し、 Determine the scattering coefficient S _w — 'and the absorption coefficient K _w ' of the first mixed glaze in which a white colorant is mixed with a colorless and transparent base glaze,

前記散乱係数 s _w' と吸収係数 K _w' の少なくとも一方を基準として、白色でない有色着色剤の吸収係数 κ _ρ と散乱係数 S _p の相対値をそれぞれ決定し、前記吸収係数 K_w，， K _p および前記散乱係数 s _w，， S _p を用いたコンビユータカラーマツチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、 Based on at least one of the scattering coefficient s _w ′ and the absorption coefficient K _w ′, the relative values of the absorption coefficient κ _ρ and the scattering coefficient S _p of the non-white colored colorant are determined, and the absorption coefficient K _w , , K _p and the scattering coefficient s _w, by performing Konbiyu Takara over pine quenching with S _p, formulation ratio of the colorant to adjust the mixed glaze having a desired color or a predetermined compounding ratio Predict the color of the mixed glaze generated by

前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，ことによって生成された混合釉で陶磁器の素地の少なくとも一部が覆われる。 Mixing the base glaze, the white colorant, and the colored colorant to have the predicted blending ratio or the predicted color, wherein at least a portion of the ceramic body is Covered.

この発明の他の局面においては，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法は、 In another aspect of the present invention, a method for determining an absorption coefficient and a scattering coefficient of a colorant used in computer color matching comprises:

無色透明でない被着色物に基準となる顔料を混合した第〗の混合物の散乱係数 s と吸収係数 K_w' を決定する工程と、 A step of determining a scattering coefficient s and an absorption coefficient K _w ′ of a second mixture obtained by mixing a reference pigment with a colorless and non-colored object;

前記散乱係数 S _w，と吸収係数 K _w，の少なくとも一方を基準として、有色着色剤の吸収係数 κ _ρ と散乱係数 S _P の相対値をそれぞれ決定する工程と、を備える。 The scattering coefficient S _w, the absorption coefficient K _w, based on the at least one, and a step of determining the relative value of the absorption coefficient kappa _[rho and scattering coefficient S _P output colored colorant, respectively, the.

この発明のさらに他の局面によれば，コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する装置は、 According to yet another aspect of the present invention, an apparatus for determining an absorption coefficient and a scattering coefficient of a colorant used in computer color matching comprises:

無色透明でない被着色物に基準となる顔料を混合した第】の混合物の散乱係数 s と吸収係数 K_w' を決定する手段と、 Means for determining the scattering coefficient s and the absorption coefficient K _w ′ of the second mixture of the non-colorless and transparent coloring matter mixed with the reference pigment;

前記散乱係数 S _w' と吸収係数 K_w' の少なくとも一方を基準として、有色着色剤の吸収係数 K _P と散乱係数 S _p の相対値をそれぞれ決定する:に程と、を備える。 On the basis of at least one of the scattering coefficient S _w ′ and the absorption coefficient K _w ′, determine the relative values of the absorption coefficient K _P and the scattering coefficient S _p of the colorant, respectively: Is provided.

また，この発明は. コンピュータカラーマッチングを利用して作成された混合釉を用いて製造された便器にも向けられている。 The present invention is also directed to a toilet bowl manufactured using mixed glaze made using computer color matching.

[図面の鲔. な説明] [A description of the drawing]

図 1は，第 1実施例における処理の全体手順を示すフローチヤ一ト。 FIG. 1 is a flowchart showing the overall processing procedure in the first embodiment.

図 2は，数式 5で与えられる散乱係数 S _w' を示すグラフ。 Fig. 2 is a graph showing the scattering coefficient S _w 'given by Eq.

図 3は，ステップ S 1の詳細手順を示すフローチャート。 Fig. 3 is a flowchart showing the detailed procedure of step S1.

図 4は，（ベース釉 +白色顔料）の物性彼を求めるためのサンプルの調合率を示す説明図。 Fig. 4 is an explanatory diagram showing the mixing ratio of a sample to determine the physical properties of (base glaze + white pigment).

図 5は，サンプルの分光反射率 R (；. ) を示す概念図。 Figure 5 is a conceptual diagram showing the spectral reflectance R (;.) Of the sample.

図 6は，ステップ S 1 4の詳細手順を示すフローチャート。 FIG. 6 is a flowchart showing the detailed procedure of step S14.

図 7は，顔料物性値決定用サンプルの調合率を示す説明図。 Fig. 7 is an explanatory diagram showing the mixing ratio of the pigment physical property value determination sample.

図 8は，第 1実施例で用いた検証用サンプルの調合率を示す説明図。 FIG. 8 is an explanatory diagram showing the mixing ratio of the verification sample used in the first embodiment.

図 9は，吸収係数 Κ _Ρ ( Λ ) の調合率 C_pに対する依存性の一例を示すグラフ。図 1 0は，第 1実施例における調合率の予測結果を示す説明図。 Figure 9 is a graph showing an example of the dependence of the absorption coefficient Κ _Ρ (Λ) on the blending rate C _p . FIG. 10 is an explanatory diagram showing a prediction result of a mixing ratio in the first embodiment.

図 1 1は，第 2実施例における処理の全体手順を示すフローチャート。 FIG. 11 is a flowchart illustrating an overall procedure of a process according to the second embodiment.

図 1 2は，ニューラルネットワークの構成を示す説明図。 Figure 12 is an explanatory diagram showing the configuration of a neural network.

図 1 3は，ステップ S 3 1の詳細手順を示すフローチャート。 FIG. 13 is a flowchart showing the detailed procedure of step S31.

図 1 4は，コンピュータカラーマッチングで予測対象範囲 P Aと複数のサンプル M 1〜M 7の三刺激値の分布を示す概念図。 Fig. 14 is a conceptual diagram showing the distribution of tristimulus values in the prediction target range PA and multiple samples M1 to M7 in computer color matching.

図 1 5は，第 1実施例における 7つのサンプル M 1〜M 7の予測誤差を示す概念図。 FIG. 15 is a conceptual diagram showing prediction errors of seven samples M1 to M7 in the first embodiment.

図 1 6は，各サンプルの三刺激値の目標値（ステップ S 4 1で決定された値）と予測誤差 AMi (Δ Χ, Δ Υ, Δ Ζ) を示す説明図。 Figure 16 is an explanatory diagram showing the target tristimulus values (the values determined in step S41) and the prediction error AMi (ΔΧ, ΔΥ, ΔΖ) for each sample.

図 1 7は，第 2実施例におけるニューラルネットワークの学習の実証結果を示す説明図。図 1 8は，第 2実施例におけるコンピュータカラ一マッチングの予測の実証結果を示す説明図。 Figure 17 is an explanatory diagram showing the results of the neural network learning in the second embodiment. Fig. 18 is an explanatory diagram showing the verification results of the prediction of computer color matching in the second embodiment.

図 1 9は，設計された色（標準色）と実際に製造される陶器の色のばらつきを示す X— y色度図。 Fig. 19 is an X-Y chromaticity diagram showing the variation between the designed colors (standard colors) and the colors of the ceramics actually manufactured.

図 2 0は，第 3実施例における処理の全体手順を示すフ π—チャート。 FIG. 20 is a π-chart showing the overall procedure of the process in the third embodiment.

図 2 1は，ステップ S 5 7におけるコンピュータカラ一マッチングの詳細手順を示すフローチヤ一ト， Figure 21 is a flowchart showing the detailed procedure of computer color matching in step S57,

図 2 2は，第 3実施例において用いた檫準色見本の三刺激値と溢淡限度見本の三刺激値を示すテーブル。 FIG. 22 is a table showing the tristimulus values of the reference color sample and the tristimulus values of the overflow limit sample used in the third embodiment.

図 2 3は，第 3実施例による渔淡限度見本の三刺激値の予測結果と比较例と実調合率とを示すテーブル。 Fig. 23 is a table showing the results of the prediction of tristimulus values of the light limit sample, comparative examples, and the actual mixing ratio according to the third embodiment.

図 2 4は，各実施例のコンピュータカラ一マッチング方法を実施するための装置を示すブ口ック図。 Figure 24 is a block diagram showing the equipment for implementing the computer color matching method of each embodiment.

図 2 5は，第 4実施例のコンビュ一タカラーマッチング方法における処理の全体手順を示すフローチヤ一ト， FIG. 25 is a flowchart showing the overall procedure of the process in the computer color matching method of the fourth embodiment.

図 2 6は，図 2 5のステップ S 7 6の誶細処理を示すフローチヤ一ト。 FIG. 26 is a flowchart showing the detailed processing of step S76 in FIG.

図 2 7は、ステップ S 7 2で取得した目標色見につレヽての U激値（色値）とステップ S 7 3で取得した第 1回目試作釉についての三剌瀲値の対比、並びに目標色見本釉，第 1回目試作釉における各顔料の調合率（調合率）を表わすテ一ブル。 Figure 27 shows the comparison between the extreme U value (color value) for the target color cast obtained in step S72 and the value of the tristimulus for the first prototype glaze obtained in step S73, and Table showing the mixing ratio (mixing ratio) of each pigment in the target color sample glaze and the first prototype glaze.

図 2 8は、図 2 7に掲げる第 1回目試作釉に各顔料をそれぞれ微量だけずつ追力 U調合した場合の三刺激値の変化率（微係数）を示すテーブル。 Fig. 28 is a table showing the rate of change (derivative coefficient) of tristimulus values when each pigment is added to the first trial glaze shown in Fig. 27 by a small amount of each pigment.

図 2 9は、第 4実施例によるコンピュータカラ一マッチング方法による結果を示すテーブル。 FIG. 29 is a table showing the result of the computer color matching method according to the fourth embodiment.

図 3 0は、第 4実施例による他の結果を示すテーブル。 FIG. 30 is a table showing another result according to the fourth embodiment.

[発明を実施するための最良の形態] [Best Mode for Carrying Out the Invention]

A. 第 1実施例：図 1は、第】実施例における処理の全体手順を示すフローチャートである。なお、この第 1実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉 (ゆう）である。すなわち、顔料を入れないベ一ス釉（基礎釉）が被着色物であり、このべ一ス釉に顔料を添加した釉がコンピュータカラーマツチングの対象となる混合物である。 A. First embodiment: FIG. 1 is a flowchart showing the overall procedure of the process in the embodiment. The mixture of interest in the first embodiment is a glaze for covering the surface of the ceramic body. In other words, the base glaze (base glaze) containing no pigment is the object to be colored, and the glaze obtained by adding a pigment to this base glaze is a mixture to be subjected to computer color matching.

なお，陶磁器において，白色顔料としては乳濁剤が有効である。乳濁剤としては，ジルコン等のジルコニウム化合物や，リン酸カルシウム等のリン化合物が用いられる。 In ceramics, an emulsifier is effective as a white pigment. Emulsifiers include zirconium compounds such as zircon and phosphorus compounds such as calcium phosphate.

ステップ S 〗では、被着色物（ベース釉）と白色顔料の混合物の吸収係数 K _w_ ' と散乱係数 S _w' とを求める。ここで、従来の相対法と異なるのは、被着色物を無色透明と仮定していない点、および、被着色物や白色顔料単独の吸収係数や散乱係数を求めずに、その混合物の吸収係数 K_w' と散乱係数 s _w' を求める点にある。なお、散乱係数 S _w' は、白色顔料の調合率 C_w に依存した関数 f ( C_w In step S 吸収, the absorption coefficient K _w _ 'and the scattering coefficient S _w ' of the mixture of the object to be colored (base glaze) and the white pigment are determined. Here, the difference from the conventional relative method is that the material to be colored is not assumed to be colorless and transparent, and the absorption coefficient of the mixture is determined without calculating the absorption coefficient or scattering coefficient of the material to be colored or the white pigment alone. The point is to find K _w 'and the scattering coefficient s _w '. Incidentally, the scattering coefficient S _w 'is a function depending on the compounding ratio C _w of the white pigment f (C _w

) として求める。また、吸収係数 κ„' も、白色顔料の調合率 c_w に依存した形式で求める。 ). Also, the absorption coefficient κ „'is determined in a form that depends on the white pigment blending ratio c _w .

なお、ステップ S 1において，白色顔料以外の顔料（以下、「有色顔料」と呼ぶ）の吸収係数 K _p と散乱係数 S _p も決定される（ステップ S 2 ) 。なお、有色顔料の吸収係数と散乱係数は、その調合率 c_p に依存した形式で求められる。 Note that, in step S 1, pigment other than white pigment (hereinafter, "colored pigments" and hump) absorption coefficient K _p and scattering coefficient S _p of is also determined (step S 2). The absorption and scattering coefficients of the chromatic color pigment is determined Me in a format depending on the compounding ratio c _p.

ステップ S 3では、ステップ S l， S 2で求めた吸収係数と散乱係数に基づいて、上述した数式 1および数式 2を用いてコンピュータカラーマッチングを実行し、混合物の色予測や調合割合の予測を行なう。 In step S3, computer color matching is performed using the above-described equations 1 and 2 based on the absorption coefficient and the scattering coefficient obtained in steps S1 and S2, and the color prediction of the mixture and the prediction of the blending ratio are performed. Do.

ここで、ステップ S 1の詳細手順を説明する前に、ベース釉と白色顔料の混合物の散乱係数 S _w' を表現するための関数 f ( C_w ) について考察する。ベース釉と白色顔料の混合物については、数式 1を次の数式 3のように書き換えることができる。 Here, before describing the detailed procedure of step S1, a function f (C _w ) for expressing the scattering coefficient S _w ′ of the mixture of the base glaze and the white pigment will be considered. For a mixture of base glaze and white pigment, Equation 1 can be rewritten as Equation 3 below.

(C _W + C_B)K_W ' = C_wKw + RKB . . . ( 3a ) (C _W + C _B ) K _W '= C _w Kw + RKB... (3a)

(C、_V + C_B)S_W ' = C_wSw +し BSR — ( 3b ) ここで、 K_w，， S_w' は混合物の吸収係数と散乱係数、 K_w ， S_w は白色顔料単独の吸収係数と散乱係数、 K_B , S_B はベース釉単独の吸収係数と散乱係数、 c_w は白色顔料の調合率、 c_B はべ一ス釉の調合率である。なお、調合率 c_w， c_B は正確には体積率で表現されるべきであるが、重量率で表現してもその誤差は通常無視できる程度である。 (C, _V + C _B ) S _W '= C _w Sw + BSR — (3b) Here, K _{_w,,} S _w 'is the absorption coefficient and scattering coefficient of the mixture, K _{_w,} S _w is the absorption coefficient and scattering coefficient of the white pigment alone, K _{_B,} S _B is the scattering and absorption coefficient of the base glaze alone coefficient , C _w is the mixing ratio of the white pigment, and c _B is the mixing ratio of the base glaze. It should be noted that the mixing ratios c _w and c _B should be expressed accurately by volume ratio, but even if expressed by weight ratio, the error is usually negligible.

数式 3を変形すると、混合物の散乱係数 S_w' は次の数式 4で与えられる。 , _ Cw S_W+レ RSR . _Λ、 By transforming Equation 3, the scattering coefficient S _w 'of the mixture is given by Equation 4 below. , _ Cw S _W + Les _RSR. Λ,

^Sw- c_w+c_B ·'·^{(4 )} この発明では相対法を用いるので、白色顔料の調合率 c_wによらず散乱係数 s_w ^Sw -c _{w +} c _B · '· ^{(4) Since the} present invention uses the relative method, the scattering coefficient s _w is independent of the mixing ratio c _w of the white pigment.

= 1 と仮定する。また、通常はベース釉の調合率 C_B を一定とし、白色顔料の調合率 c_wを変化させて混合物を作成するので、 C_B =〗 0 0 (—定）と仮定する。この時、白色顔料の調合率 „ は、いわゆる外添加調合率（ベ一ス釉の重量を 1 0 0 %とした時の顔料の重量％) となる。これらの仮定をおくと、上記数式 4は次の数式 5のように害き換えられる。 C_w+ 100S_B = 1 Also, usually, the blending ratio C _B of the base glaze is constant and the blending ratio c _w of the white pigment is changed to form a mixture. Therefore, it is assumed that C _B =〗 0 0 (—constant). At this time, the mixing ratio の of the white pigment is the so-called external addition mixing ratio (weight% of the pigment when the weight of the base glaze is 100%). Is harmed by the following equation 5. C _w + 100S _B

^Sw'⁼ Cw+ 100 · ' · ^{( 5 )} ベース釉は無色透明に近いが、完全な無色透明ではない。すなわち、数式 5において、ベ一ス釉の散乱係数 S_B は 0でない正の小さな値である。従って、ベース釉と白色顔料の混合物の散乱係数 S_w' は、図 2 (A) に示すようなグラフとなる。このグラフは、白色顔料の調合率 C_w が大きくなると S_w' =〗に漸近する双曲線である。 ^Sw ' ⁼ Cw + 100 ·' · ^{(5) The} base glaze is almost colorless and transparent, but not completely colorless and transparent. That is, in Equation 5, the scattering coefficient S _B of the base glaze is a small positive non-zero value. Therefore, the scattering coefficient S _w 'of the mixture of the base glaze and the white pigment becomes a graph as shown in Fig. 2 (A). This graph is a hyperbola that gradually approaches S _w ′ =〗 as the blending ratio C _w of the white pigment increases.

ところで、図 2 (A) に示すような散乱係数 S_w' を実測することは容易ではない。そこで、この第 1実施例では、白色顔料の調合率 C_wに対する散乱係数 S_w— ，のグラフの形状を、図 2 (B) に示すように高濃度側の第 1の直線 L 1と低濃度側の第 2の直 L 2で近似できるものと仮定した。直線し L 2を表わす関数 S = f , (C_w ) , S = f ₂ (C_w ) は、それぞれ次の数式 6および 7で表わすことができる。 S_w' = aC_w + b -..(6) (C„。≤ C_w) By the way, it is not easy to measure the scattering coefficient S _w ′ as shown in Fig. 2 (A). Therefore, in the first embodiment, the shape of the graph of the scattering coefficient S _w —, with respect to the mixing ratio C _w of the white pigment is changed to the first straight line L 1 on the high concentration side as shown in FIG. It is assumed that it can be approximated by the second straight line L 2 on the concentration side. The functions S = f, (C _w ) and S = f ₂ (C _w ), which are linear and represent L 2, can be expressed by the following equations 6 and 7, respectively. S _w '= aC _w + b-.. (6) (C „.≤ C _w )

S_w' = d C_w + e ... (7) (C_K < C_wo) S _w '= d C _w + e ... (7) (C _K <C _wo )

ここで，は直線 L I, L 2の交点における白色顔料の調合率である。 Here, is the blending ratio of the white pigment at the intersection of the straight lines L I and L 2.

数式 6， 7における係数 a， b, d, eを決定すれば、（ベ一ス釉 +白色顔料）の混合物に対する散乱係数 S_w，が求まる。そして、散乱係数 S_w，が求まれば、分光反射率 R (A ) の測定値を用いて上記数式 2から吸収係数 K_w' を求めることができる。 By determining the coefficients a, b, d, and e in Equations 6 and 7, the scattering coefficient S _w , for the mixture of (base glaze + white pigment) is determined. Then, if the scattering coefficient S _w , is obtained, the absorption coefficient K _w ′ can be obtained from Equation 2 using the measured value of the spectral reflectance R (A).

図 3は、図 1のステップ S 1の詳細手順を示すフローチャートである。ステツプ S 1 〗では、ベース釉と白色顔料のみを混合した混合物（第 1の混合物）を作成した。図 4は、第〗実施例において準備した 1 3個のサンプル W】 2〜W0の調合率を示している。図 4から解るように、 1 3個のサンプル Wl 2〜W0は、白色顔料の調合率 C_w を 1 2%〜0%の範囲で 1 %ずつ変化させて作成したものである。なお、第 1実施例におけるサンプルは、通常の陶磁器の素地に、調合した釉をかけて焼成することによって作成した。また、後述する他のサンプルも同じ条件で作成した。 FIG. 3 is a flowchart showing a detailed procedure of step S1 in FIG. In step S 1。, a mixture (first mixture) was prepared in which only the base glaze and the white pigment were mixed. FIG. 4 shows the mixing ratios of 13 samples W] 2 to W0 prepared in the first embodiment. As can be seen from FIG. 4, 1 3 samples Wl 2～W0 are those created by the compounding ratio C _w of the white pigment is changed by 1% in the range of 1 2% 0%. The sample in the first example was prepared by baking a normal ceramic body with a blended glaze. In addition, other samples described later were prepared under the same conditions.

ステップ S 1 2では、各サンプル W〗 2〜W0について、分光光度計で分光反射率 R' (A ) を測定した。ステップ S 1 3では、分光反射率の測定値 R， In step S12, the spectral reflectance R '(A) of each sample W2 to W0 was measured with a spectrophotometer. In step S13, the measured values of spectral reflectance R,

(λ ) に基づいて、数式 6に舍まれる係数 a， bを以下のように実験的に決定した。 Based on (λ), the coefficients a and b in Equation 6 were experimentally determined as follows.

白色顔料の調合率 C_w が比較的大きなサンプルでは、釉の隠ぺぃ力が大きいので、 S_w' = 1と仮定することができる。具体的には、白色顔料の調合率 C_w が最も大きなサンプル W1 2に対する散乱係数 S_w' を 1に等しいと仮定した。この仮定を用いると、白色顔料の調合率 C_w が i %のサンプル W iの散乱係数 S_wi_ ' は、次の数式 8で与えられる。 In a sample having a relatively large blending ratio C _w of the white pigment, the hiding power of the glaze is large, so it can be assumed that S _w ′ = 1. Specifically, it was assumed that the scattering coefficient S _w ′ for the sample W 12 having the largest white pigment preparation ratio C _w was equal to 1. Using this assumption, the scattering coefficient S _wi ′ of the sample W i having the white pigment preparation rate C _w of i% is given by the following equation (8).

。w i = Swn ' ...(8). w i = Swn '... (8)

CT/JP96/00738

CT / JP96 / 00738

18 数式 8において、（KZS) _W12 は白色顔料の調合率が 1 2 %のサンプル W 1 2の吸収係数と散乱係数の比であり、（KZS) _wiは白色顔料の調合率が i %のサンプル W iの吸収係数と散乱係数の比である。また、数式 8では、吸収係数の変化は散乱係数の変化に比べて小さく、 K_{wl 2} ==K_{w i}であると仮定している。数式 8の右辺の（KZS) _wl2 ， (K/S) _wiはサンプル W iの分光反射率の測定 :R' ( J から、上述した数式 2に従って求めることができる。似し、混合物である釉の厚みは小さいので、分光光度計で測定される分光反射率 R， (λ ) を、次の数式 9 (サンダーソン（Saunderson) の式）に従って理想状態18 In Equation 8, (KZS) _W12 is the ratio of the absorption coefficient to the scattering coefficient of the sample _W12 in which the blending ratio of the white pigment is ₁₂ %, and (KZS) _wi is the sample in which the blending ratio of the white pigment is i%. It is the ratio between the absorption coefficient and the scattering coefficient of Wi. Equation 8 assumes that the change in the absorption coefficient is smaller than the change in the scattering coefficient, and that K _{wl 2} == K _wi . (KZS) _wl2 , (K / S) _{wi on} the right side of Equation 8 can be obtained from the measurement of the spectral reflectance of the sample W i: R ′ (J according to Equation 2 described above. Because the thickness of a certain glaze is small, the spectral reflectance R, (λ) measured by a spectrophotometer is calculated according to the following equation 9 (Saunderson's equation).

(混合物の厚みが無限大の状態）の分光反射率 R (； - ) に変換する。 p - k - k】 Q (The thickness of the mixture is infinite) and is converted into the spectral reflectance R (;-). p-k-k] Q

(1 -ki)(l -k₂) + k₂R'- kjk₂ · ここで、係数， k ₂ は、被着色物（ベース釉）の光学的性質に依存する値である。係数， k ₂ としては、被着色物の屈折率 nから次の数式 1 0に従って決定することができる。

(1 -ki) (l -k ₂ ) + k ₂ R'-kjk ₂ · Here, the coefficient, k ₂ is a value that depends on the optical properties of the object to be colored (base glaze). The coefficients, k ₂ can be determined from the refractive index n of the object to be colored according to the following equation 10:

k₂=0.68n-0.56 ... (10b) k ₂ = 0.68n-0.56 ... (10b)

なお、第 1実施例において用いたベース釉は、屈折率 πが約 1. 4である。こうして理想状態の分光反射率 R (λ ) が求まると、上述の数式 2から（KZS) _wiが得られる。 The base glaze used in the first embodiment has a refractive index π of about 1.4. When the spectral reflectance R (λ) in the ideal state is obtained in this way, (KZS) _wi is obtained from the above-described equation (2).

数式 9に従って分光反射率の測定値 R' U ) から理想状態の分光反射率 R {1 ) を求める演算は、後述する他の工程においてサンプルの分光反射率から (K/S) を求める際にも同様に実施される。 The calculation to determine the ideal state spectral reflectance R (1) from the spectral reflectance measured value R 'U) according to Equation 9 is performed when calculating (K / S) from the spectral reflectance of the sample in another process described later. Is similarly performed.

各サンプル W iについて得られた（K/S) _wiの値を上記数式 8に代入することによって、各調合率 C_w において混合物の散乱係数 S„' の値を求めることができ、係数 a， bを決定することができる。係数 a, bを含む数式 6は、白色顔料の調合率 C„ が比較的髙ぃ範囲にのみ適用される。従って、係数 a， bの決定においても、調合率 C_w が比較的高い範囲の（例えば 1 0%〜 1 2%の）数個のサンプルの（K/S) _wiのみを用いる。 The obtained values of (K / S) _wi by the child substituted into the equation 8 for each sample W i, it is possible to determine the value of the scattering coefficient S "'of the mixture in each of the mixing ratio C _w, the coefficient a , B can be determined Equation 6 including the coefficients a and b is applied only to the range where the mixing ratio C „of the white pigment is relatively high. Therefore, determination of coefficients a and b In (2), only (K / S) _{wi of} several samples in which the blending ratio _Cw is relatively high (for example, 10% to 12%) is used.

図 5は、 C_w= l 0%〜〗 2%のサンプルの分光反射率 R U ) を示す概念図である。このように、ベース釉に高濃度の白色顔料を混合したものは、可視光線の波長の全範囲（約 400 n m〜約 700 n m) においてほとんど一定の反射率を示す。そこで、この分光反射率 R (λ ) の値から数式 2と数式 8を用いて散乱係数 S_wi' を求めることができる。なお、この散乱係数 S は、可視光線の全波長範两において同一の値を有する。第 1実施例では、 a = 0. 0005, b =0. 94が得られたので、数式 6が次の数式 1 1に書き換えられた。 Figure 5 is a conceptual diagram showing a spectral reflectance RU) of C _w = l 0% ~〗 2% sample. Thus, the base glaze mixed with a high concentration of white pigment exhibits almost constant reflectance over the entire visible wavelength range (about 400 nm to about 700 nm). Then, the scattering coefficient S _wi ′ can be obtained from the value of the spectral reflectance R (λ) using Expressions 2 and 8. The scattering coefficient S has the same value in the entire wavelength range of visible light. In the first embodiment, since a = 0.0005 and b = 0.94 were obtained, Expression 6 was rewritten into Expression 11 below.

S_w ' = 0.005 C_w + 0. 4 ...(11) (髙濃度側） S _w '= 0.005 C _w + 0.4 ... (11) (髙 concentration side)

なお、数式 1 1は、 C_w = 1 0%付近において S_w' = 1に極めて近い値を与えるので、数式 1 1の代わりに S_w， = 1を用いてもよい。 Since Equation 11 gives a value very close to S _w ′ = 1 near C _w = 10%, S _w , = 1 may be used instead of Equation 11.

図 3のステップ S 1 4では、低濃度側の数式 7の係数 d, eの値を決定する。白色顔料の調合率 C_w が比較的小さい場合には下地の影響が大きくなるので、分光反射率 R (λ ) の測定値が真の釉の層のみの情報ではなくなる。従って、分光反射率 R (λ ) の測定値に基づいて係数 d， eを決定するのは容易ではない。そこで、係数 d， eの値は、図 6に示すような手順に従って決定する。 In step S14 of FIG. 3, the values of the coefficients d and e of Equation 7 on the low concentration side are determined. When the blending ratio C _w of the white pigment is relatively small, the influence of the ground becomes large, so that the measured value of the spectroscopic reflectance R (λ) is not information of only the true glaze layer. Therefore, it is not easy to determine the coefficients d and e based on the measured values of the spectral reflectance R (λ). Therefore, the values of the coefficients d and e are determined according to the procedure shown in Fig. 6.

ステップ S 2 1では、顔料物性値決定用サンプルを作成する。図 7は、顔料決定用サンプルの調合率を示す説明図である。顔料物性値決定用サンプルは、ベース袖と白色顔料と他の有色顔料を混合して作成したものである。第 1実施例では、白色顔料と有色顔料の合計の調合率 (Pigment Volume Concentration, P VC) を 1 2%—定とし、有色顔料の調合率 C_p を 1 %〜 1 2%の範囲で 1 %ずつ変えた 1 2個のサンプル M】〜M 1 2を作成した。なお、有色顔料としては、陶磁器の彩色に良く用いられる青色、赤色、黄色、等を選択し、各顔料について図 7に示すような 1 2個のサンプル M〗〜M 1 2をそれぞれ作成した。 In step S21, a sample for determining the physical properties of a pigment is prepared. FIG. 7 is an explanatory diagram showing the mixing ratio of the pigment determination sample. The sample for determining the physical properties of the pigment was prepared by mixing a base sleeve, a white pigment and another colored pigment. In the first embodiment, the total compounding ratio of the white pigment and colored pigment (Pigment Volume Concentration, P VC) 1 2% - a constant, 1 formulatory ratio C _p of colored pigment in the range of 1% to 1 2% 12 samples M] to M12 were prepared. As colored pigments, blue, red, yellow, etc., which are often used for coloring ceramics, were selected, and for each pigment, 12 samples M〗 to M12 as shown in FIG. 7 were prepared.

なお、サンプル M l 2は白色顔料を舍まない混合物であり、本願発明の第 2の混合物に相当する。また、サンプル M 1〜M 1 1は本願発明の第 3の混合物に相当する。 Note that the sample Ml2 is a mixture containing no white pigment, and corresponds to the second mixture of the present invention. Samples M1 to M11 correspond to the third mixture of the present invention. Hit.

ステップ S 2 2では、検証用サンプルを作成する。 |、|8は、第 1実施例で用いた検証用サンプルの調合率を示す説明図である。 4つの検証用サンプル D 1〜D 4は、 4色の顔料をベース釉に混合したものである。なお、ベース釉の調合率 C _B は省略しているが、いずれのサンプルも C_B = 1 00%である。 In step S22, a verification sample is created. |, | 8 are explanatory diagrams showing the mixing ratios of the verification samples used in the first example. The four verification samples D1 to D4 are four-color pigments mixed in the base glaze. Incidentally, Formulation ratio C _B of the base glaze is omitted, all of the samples are also C _B = 1 00%.

検証用サンプル D〗〜D 4は、本願発明における第 4の混合物に相当する。ステップ S 2 3では、係数 d, eに仮の値を割り当てる。図 2 (A) ， (B) を比校すれば解るように、係数 eの値は、ベース釉の散乱係数 S_B の値に相当しており、かなり小さな値であることが解っている。例えば、 e = 0と仮定することも可能である。また、図 2 (B) に示す第 1 と第 2の直線 L 1， L 2の交点は、横轴の調合率 C_w が数％のところに位置するので、係数 dも概略の値を仮決めすることが可能である。 The verification samples D〗 to D4 correspond to the fourth mixture in the present invention. In step S23, temporary values are assigned to the coefficients d and e. As can be seen by comparing Figs. 2 (A) and 2 (B), the value of the coefficient e corresponds to the value of the scattering coefficient S _B of the base glaze, and is a very small value. For example, it is possible to assume e = 0. Also, since the intersection of the first and second straight lines L 1 and L 2 shown in FIG. 2 (B) is located at a position where the horizontal mixing ratio C _w is several percent, the coefficient d is also a rough value. It is possible to decide.

係数 d， eを仮決めすると、調合率 C_w に依存する散乱係数 S_wi' が調合率 C_w の全範囲において決定される。散乱係数 S_Wi' が決まると、次の数式 1 2に従つて調合率 C_w に依存する吸収係数 K_wi' も得られる。 W_i ' = ( ■) xSw,' ...(12) Coefficient d, the temporarily decided e, scattering coefficient depends on the compounding ratio C _w S _wi 'is determined in the entire range of the compounding ratio C _w. When the scattering coefficient S _Wi 'is determined, an absorption coefficient K _wi ' depending on the mixing ratio C _w is obtained according to the following equation (12). W _i '= (■) xSw,' ... (12)

、へ；なお、散乱係数 S_wi' の値は波長スに依存しないが、吸収係数 K_w,' の値は波長に依存する。この理由は、この発明の相対法においては散乱係数 s_wi' を波長スに依存しない基準値として用い、他の物性値の相対値を求めるからである。分光反射率 R U) の値はもちろん波長スに依存するので、数式 2に従って得られる（KZS) _wiの愤も波長スに依存する。従って、数式 1 2に従って得られる吸収係数 K_wi' も波長スに依存する。換言すれば、吸収係数 K_wi' は白色顔料の調合率 C_w と波長とに依存する形式で求められる。 The value of the scattering coefficient S _wi 'does not depend on the wavelength, but the value of the absorption coefficient K _w ,' depends on the wavelength. The reason for this is that in the relative method of the present invention, the scattering coefficient s _wi ′ is used as a reference value that does not depend on the wavelength, and the relative values of other physical properties are obtained. Since the value of the spectral reflectance RU) depends of course on the wavelength, the value of (KZS) _wi obtained according to Equation 2 also depends on the wavelength. Therefore, the absorption coefficient K _wi 'obtained according to Equation 12 also depends on the wavelength. In other words, the absorption coefficient K _wi 'is obtained in a form that depends on the blending ratio C _w and the wavelength of the white pigment.

ステップ S 25では、顔料物性値決定用サンプルを用いて、各有色顔料の吸収係数 Κ_Ρ (λ ) と散乱係数 S _p (λ ) とを求める。なお、ステップ S 25は、図 1におけるステップ S 2に相当する。 ( . ) と S _P (ス）は、次のような手順で算出される。まず、顔料物性値決定用サンプルについては、ヒ述の数式 1から次の数式 1 3 が得られる。 In step S25, the absorption coefficient Κ _Ρ (λ) and the scattering coefficient S _p (λ) of each colored pigment are determined using the pigment property value determination sample. Step S25 corresponds to step S2 in FIG. (.) And _SP (s) are calculated by the following procedure. First, for a pigment physical property value determination sample, the following expression 13 is obtained from expression 1 described above.

(CR + ρ)Κ\ι _ ^_ΒΚ_Β + v^ K\v + Cprvp (1 (CR + ρ) Κ \ ι _ ^ _Β Κ _Β + v ^ K \ v + Cprvp (1

(C B + C\v + Cp)a.\i CB¾ W W + Cp ところで、上記の数式 3は、次の数式 1 4のように害き換えられるので、数式 1 4を数式〗 3に代入すると、数式 1 5が得られる。 (CB + C \ v + Cp) a. \ I CB¾ WW + Cp By the way, the above formula 3 can be replaced by the following formula 14.So, when formula 14 is substituted into formula〗 3, Equation 15 is obtained.

Cw T V ' = CB¾ + Cw \v · · · 4a ) Cw T V '= CB¾ + Cw \ v 4a)

C_w 'S w ¹ = C_BS_R + C_WS_W . - - (14b) C _w 'S w ¹ = C _B S _R + C _W S _W .--(14b)

Cw * = C_B + C_w ... (14c) Cw * = C _B + C _w ... (14c)

KM _ C\y ' Kv ' + C K · ( ) KM _ C \ y 'Kv' + C K · ()

SM C\y ' Sy ' + CpSp 数式〗 5を変形すると、次の数式 1 6が得られる。 SM C \ y 'Sy' + CpSp By transforming equation〗 5, the following equation 16 is obtained.

ここで、（KZS) _p は有色顔料単独での吸収係数と散乱係数の比であり、 (K/S) は顔料物性値決定用サンプルの吸収係数と散乱係数の比、（KZ S) は（ベース釉 +白色顔料）のサンプルの吸収係数と散乱係数の比である。また、 c_w，はベース釉と白色顔料の調合率の合計値、 s は（ベース釉 +白色顔料）の混合物の散乱係数であり、数式 6および 7で与えられる値である。数式 1 6の右辺の各項は、以下のように求めることができる。

Here, (KZS) _p is the ratio between the absorption coefficient and the scattering coefficient of the colored pigment alone, (K / S) is the ratio between the absorption coefficient and the scattering coefficient of the pigment physical property determination sample, and (KZ S) is ( This is the ratio between the absorption coefficient and the scattering coefficient of the sample (base glaze + white pigment). C _w , is the sum of the blending ratios of the base glaze and the white pigment, and s is the scattering coefficient of the mixture of (base glaze + white pigment), which is given by Equations 6 and 7. The terms on the right side of Equation 16 can be obtained as follows.

(K/S) _p の倘は、白色顔料を舍まない顔料物性値決定用サンプル（図 7のサンプル M l 2) の分光反射率 R (λ ) から上記数式 2によって求めることができる。但し、厳密に言えば、数式 1 6における（KZS) の値はベース釉の影響を含まない有色顔料^独に対する値である。一方、上述のようにして実測で得られる（KZS) _ρ は（ベース釉 +有色顔料）のサンプルに対する他なので、ベース釉の影響が舍まれている。し力し、 17のサンプル M 1 2の顔料の調合率は 1 2%と卨ぃ値なので、（KZS) p に対するベース釉の寄与は極めて小さい。従って、数式 1 6における（K/S) _p の値としては、サンプル M 1 2の分光反射率 R (λ ) から数式 2によって求めた値を用いても、その誤差は無視できる程度である。 K of (K / S) _p can be obtained by the above equation 2 from the spectral reflectance R (λ) of the pigment physical property determination sample (sample Ml 2 in FIG. 7) which does not contain white pigment. However, strictly speaking, the value of (KZS) in Equation 16 is a value for a colored pigment that does not include the effect of the base glaze. On the other hand, (KZS) _ρ obtained by actual measurement as described above is for the sample of (base glaze + colored pigment). The influence of the glaze is impaired. However, since the mixing ratio of the pigment of 17 samples M12 is 12%, which is a 卨ぃ value, the contribution of the base glaze to (KZS) p is extremely small. Therefore, as the value of (K / S) _p in Equation 16, even if the value obtained from Equation 2 from the spectral reflectance R (λ) of sample M 12 is used, the error is negligible. is there.

(KZS) _Μ の値は、白色顔料と有色顔料の両方を舍む顔料物性値決定用サンプル（図 7のサンプル Μ 1〜！ VI 1 1 ) の分光反射率 R (；. ) から数式 2に従って算出できる。従って、（KZS) _Μ は有色顔料の調合率 C_p に依存した形式で求められる。 (KZS) value of _Micromax is the spectral reflectance of the white pigment and Complex both colored pigments unpigmented physical properties determined for samples (samples of Figure 7 Μ 1~ VI 1 1!) R (;.) From the formula 2 Can be calculated according to Therefore, (KZS) _Μ is determined in a form depending on the mixing ratio C _p of the colored pigment.

(K/S) _w' の値は、 4に示す（ベース釉 +白色顔料）のサンプルの分光反射率 R (ス）力 ^数式 2に従って算出された値を用いることができる。 As the value of (K / S) _w ′, the value calculated according to the spectral reflectance R (s) power of the sample of (base glaze + white pigment) shown in 4 ^ Equation 2 can be used.

以上のようにして数式 1 6の右辺の各項の値を求めることができるので、顔料物性値決定用サンプル Ml〜M 1 1のそれぞれについて、散乱係数 S_p (λ ) を求めることができる。なお、散乱係数 S_p (Λ) は、有色顔料の調合率 C_p に依存する形式で求められる。 It is possible to determine the value of each term of the to the right side of Equation 1 6 as described above, for each of the pigment physical properties determined for the samples Ml～M 1 1, can be obtained scattering coefficient S _p a (lambda). The scattering coefficient S _p (Λ) is obtained in a form depending on the mixing ratio C _p of the colored pigment.

こうして得られた散乱係数 S_p (λ ) と、分光反射率 R ( .) の測定値から数式 2で算出された（KZS) ρ の値とに基づいて、次の数式 1 7に従って吸収係数 Κ_Ρ {?. ) も求められる。 Based on the scattering coefficient S _p (λ) thus obtained and the value of (KZS) ρ calculated from the measured value of the spectral reflectance R (. The number Κ _Ρ {?.) Is also required.

(17) (17)

. なお、吸収係数 Κ_ρ (λ ) は、散乱係数 S _p (λ ) と同様に、有色顔料の調合率 C_p に依存する形式で求められる。 . The absorption coefficient Κ ρ _(λ), as well as the scattering coefficient S _p (λ), obtained in a form which depends on the compounding ratio C _p of colored pigments.

図 9は、吸収係数 Κ_Ρ (Λ) の調合率 C_p に対する依存性の一例を示すグラフである。散乱係数 S _p (； J もこれと同様な依存性を示す。なお、図 9では、代表的な波長についてのグラフのみを示しているが、実際には可視光線の全波長範囲（約 400 nm〜約 7 00 n m) で〗 0 n m毎に吸収係数 Κ _ρ (λ ) が求められる。以上のようにして、（ベース釉 +白色顔料）の混合物の吸収係数 K_w，（ス）および散乱係数 S_w' と、他の有色顔料の吸収係数 K_P ）および散乱係数 S _p (； J が求まると、図 6のステップ S 26において、検証用サンプル（図 8) に対してコンピュータカラーマッチングによるシミュレーションを行なう。 FIG. 9 is a graph showing an example of the dependence of the absorption coefficient Κ _Ρ (Λ) on the mixing ratio C _p . The scattering coefficient S _p (; J shows a similar dependence. Note that FIG. 9 shows only a graph of typical wavelengths, but in actuality, the actual wavelength range of visible light (about 400 nm~ about 7 00 nm) absorption coefficient for each〗 0 nm in Κ ρ _(λ) is obtained. As described above, the absorption coefficient K _w , (s) and the scattering coefficient S _w ′ of the mixture of (base glaze + white pigment) and the absorption coefficient K _p of other colored pigments and the scattering coefficient S _p (; J When is obtained, in step S26 in FIG. 6, a simulation is performed on the verification sample (FIG. 8) by computer color matching.

ここで、コンピュータカラーマッチングによる色予測と調合割合の予測について簡単に説明する。数式 2を変形すると、混合物の分光反射率 R ）は、次の数式 1 8で与えられる。 Here, the color prediction by computer color matching and the prediction of the mixing ratio will be briefly described. By transforming Equation 2, the spectral reflectance R) of the mixture is given by Equation 18 below.

(18) (18)

' MM

混合物の吸収係数と散乱係数の比（KZS) _M は、数式 1に従って算出できるので、混合物の分光反射率 R (λ ) も数式 1 8から求めることができる。この分光反射率 R (λ ) は理想状態の分光反射率なので、分光光度計で測定できる分光反射率 (λ) を、数式 9 (サンダーソンの式）を変形した次式 19に従つて求める。 '' MM

Since the ratio (KZS) _M of the absorption coefficient and the scattering coefficient of the mixture can be calculated according to Equation 1, the spectral reflectance R (λ) of the mixture can also be determined from Equation 18. Since this spectroscopic reflectance R (λ) is an ideal state spectral reflectance, the spectral reflectance (λ) that can be measured by a spectrophotometer is calculated according to the following Equation 19, which is a modification of Equation 9 (Sanderson's equation). Ask.

R' = ki + (l-ki)(l-k₂)_] ...(19) 係数， k₂ は上記数式〗 0で与えられる。 R ′ = ki + (l-ki) (lk ₂ ) _] ... (19) The coefficient, k ₂ is given by the above equation〗 0.

分光反射率 R' (λ ) が求まると、混合物の三刺激値 X， Υ, Ζが次の数式 2 0によって求められる。 Once the spectral reflectance R '(λ) has been determined, the tristimulus values X, Υ, Ζ of the mixture are determined by the following equation 20.

• .. (20a) • .. (20a)

Y = k(S( )R )yWd .. - (20b) Z=kjS( )R'( )z( )dX ••■(20c) k =， 100 ...(20d) Y = k (S () R) yWd ..-(20b) Z = kjS () R '() z () dX •• ■ (20c) k =, 100 ... (20d)

ここで、 S ）は標準光の分光分布、 X (え）， y (λ ) , ζ (λ ) (数式中ではバ一付きである）は等色関数である。コンピュータカラーマッチングによって混合物の色を予測する場合には、上述のように、混合物の各成分の吸収係数と散乱係数に基づき、数式数式 1 8〜 20に従って、その混合物の三刺激値 X， Υ, Zを算出する。三刺激値は混合物の色を表わすので、任意の混合物の色を予測することができることになる。また、所望の色を有する混合物の調合割合を予測する場合には、混合物の調合割合を仮定して、上述した手順でその三刺激値を算出し、ニュートン一ラプソン法などの逐次近似法によって所望の色に所定の誤差内で一致するような調合割合を求める。 Here, S) is the spectral distribution of the standard light, and X (e), y (λ), ζ (λ) (in the formulas, with a dash) are color matching functions. When predicting the color of a mixture by computer color matching, as described above, based on the absorption coefficient and the scattering coefficient of each component of the mixture, the tristimulus values X, Υ, Calculate Z. Since the tristimulus values represent the color of the mixture, the color of any mixture can be predicted. Further, when predicting the blending ratio of a mixture having a desired color, the tristimulus value is calculated by the above-described procedure assuming the blending ratio of the mixture, and the desired value is calculated by a successive approximation method such as the Newton-Raphson method. Find a blending ratio that matches the color of the product within a specified error.

なお、丫∑表色系の代ゎりにし* 3*1)* 表色系のような他の表色系を用いてコンピュータカラ一マッチングを行なうことも可能である。 Note that it is also possible to perform computer color matching using another color system such as * 3 * 1) * color system instead of the color system.

図 6のステップ S 26におけるシミュレーションでは、コンピュータカラーマツチングによって、図 8に示す各検証用サンプル D 1〜D 4の調合率を予測する。検証用サンプルの調合率は既知なので、予測した調合率と実調合率との一致度は簡単に算出できる。例えば、調合率の一致度の指標として、次の数式 2 1で与えられる自乗平均誤差厶を使用する。 In the simulation in step S26 in FIG. 6, the compounding ratio of each of the verification samples D1 to D4 shown in FIG. 8 is predicted by computer color matching. Since the mixing ratio of the verification sample is known, the degree of coincidence between the predicted mixing ratio and the actual mixing ratio can be easily calculated. For example, the root mean square error given by the following equation 21 is used as an index of the degree of coincidence of the mixing ratio.

> (CR.-C,)² > (CR.-C,) ²

^Δ=-^ · '·(²¹) ^{Δ = - ^ · '· (} 21)

ここでは成分 iの実調合率、 Ci は予測された調合率である。ステップ S 26において、各検証用サンプル D 1〜D 4について得られた自乗平均誤差八が所定の許容値以下でない場合にはステップ S 2 3に戻り、係数 d， eの値を修正する。こうして、各検証用サンプル D 1〜D 4について得られる自乗平均誤差△ が所定の許容値以下になるまでステップ S 23〜S 26を繰り返し実行し、この結果、係数 d， eを決定することができる（ステップ S 2 7) 。第 1実施例では、 d = 0. 1 225, e = 0となり、数式 7が次の数式 22書き換えられた。 Where is the actual blending ratio of component i and Ci is the predicted blending ratio. In step S26, if the root mean square error 8 obtained for each of the verification samples D1 to D4 is not equal to or smaller than the predetermined allowable value, the process returns to step S23 to correct the values of the coefficients d and e. Thus, steps S23 to S26 are repeatedly executed until the root mean square error について obtained for each of the verification samples D1 to D4 is equal to or smaller than a predetermined allowable value. As a result, the coefficients d and e can be determined. Yes (step S27). In the first embodiment, d = 0.125, e = 0, and Expression 7 is rewritten as Expression 22 below.

S_w ' = 0.1225 C_w . - . (22) (低浪度側） _{_{S w '= 0.1225 C w -}} .. (22) ( low Sina degree side)

なお、数式 1 〗で与えられる高渙度側の散乱係数 S_w' のグラフ（図 2 (B) の直線 L I) と、数式 22で与えられる低濃度側の散乱係数 S_w' のグラフ（直線 L 2) との交点における調合率 C_w は 8. 0%である。 Note that the graph of the scattering coefficient S _w 'on the high lysis level given by Equation 1 (Fig. 2 (B) The mixing ratio C _{w at} the intersection of the straight line LI) of this and the graph of the scattering coefficient S _w ′ on the low concentration side given by Equation 22 (the straight line L 2) is 8.0%.

以上のように、上記第 1実施例では検証用サンプルの調合率を精度良く予測できるように数式 7における係数 d， eを決定したので、 A色顔料の調合率 C_w が比較的低い領域においても、（ベース釉 +白色顔料）の混合物の散乱係数 S_w' と吸収係数 K_w' (λ) とを精度良く決定することができた。 As described above, coefficients in Equation 7 as cut with accurately predicting the formulation ratio of the verification samples in the first embodiment d, since the determined e, Formulation ratio C _w of the A color pigment is relatively low area In, the scattering coefficient S _w 'and the absorption coefficient K _w ' (λ) of the mixture of (base glaze + white pigment) could be determined with high accuracy.

なお、図 6に示す手順で（ベ一ス釉 +白色顔料）の混合物の散乱係数 S_w' と吸収係数 K (/) とが求められる際には、ステップ S 25において他の有色顔料の吸収係数 Κ_Ρ (λ ) と散乱係数 S_p (λ ) も同時に求められる。 When the scattering coefficient S _w ′ and the absorption coefficient K (/) of the mixture of (base glaze + white pigment) are determined by the procedure shown in FIG. 6, the absorption of other colored pigments is determined in step S25. coefficient Κ Ρ _(λ) and scattering coefficient S _p (λ) is also required at the same time.

こうして各成分の吸収係数と散乱係数が求まると、図 1のステップ S 3において、コンピュータ力ラーマッチングによる色予測や調合割合の予測が実行される。図 1 0は、第 1実施例における調合率の予測結果を示す説明図である。図 1 0における第 1実施例の予測結果は、数式 1 1および数式 22を用いた場合の結果である。また、比較例の予測結果は、 S„' = 0 (一定）の仮定（図 2 (B) において破線で示すグラフ）を用いた場合の結果である。なお、図 1 0の第 1実施例と比較例では、 L*a*b* 表色系を用いて色差が最小となる調合率を予測した。 When the absorption coefficient and the scattering coefficient of each component are obtained in this way, in step S3 in FIG. 1, color prediction and blend ratio prediction by computer color matching are executed. FIG. 10 is an explanatory diagram showing a prediction result of the mixing ratio in the first embodiment. The prediction result of the first embodiment in FIG. 10 is a result obtained when Expressions 11 and 22 are used. In addition, the prediction result of the comparative example is a result when the assumption of S „′ = 0 (constant) (a graph shown by a broken line in FIG. 2B) is used. In the examples and comparative examples, the blending ratio that minimizes the color difference was predicted using the L * a * b * color system.

図 1 0に示されるように、第]実施例による予測結果は調合率の自乗平均誤差が比較例に比べてかなり小さく、より高精度で調合率を予測できたことが解る。特に、サンプル P 2， P 3のように白色顔料の調合率が小さな場合には、比較例の予測精度がかなり悪化しているのに対して、第 1実施例の予測精度は良好である。この理由は、図 2 (A) ， (B) に示すように、白色顔料の調合率 C_w が小さくなるに従って、実際の散乱係数 S_w，の値が 1から離れていくので、 S_w， = 1 (一定）と仮定した比較例では、各成分の吸収係数と散乱係数の値の誤差が大きくなるからであると考えられる。 As shown in FIG. 10, the prediction result of the Example] has a significantly smaller root mean square error of the mixing ratio than that of the comparative example, indicating that the mixing ratio can be predicted with higher accuracy. In particular, when the mixing ratio of the white pigment is small as in the samples P2 and P3, the prediction accuracy of the comparative example is considerably deteriorated, whereas the prediction accuracy of the first embodiment is good. . This is because, as shown in FIG. 2 (A), (B) , in accordance with Formulation ratio C _w of the white pigment is small fence, the actual scattering coefficient S _w, the value of moves away from 1, S _w , = 1 (constant), it is considered that the error between the absorption coefficient and the scattering coefficient of each component is large.

以上のように、この第 1実施例では、無色透明でない被着色物であるべ一ス釉の光学的影響を考慮して各成分の吸収係数と散乱係数とを決定しているので、コンピュータカラ一マッチングによつて高精度の予測を行なうことが可能である。また、サンプルとしては図 4、図 7、図 8に示すような混合物を通常の素地の上に塗布したものを準備すればよいので、サンプルの作成も容易である。 As described above, in the first embodiment, since the absorption coefficient and the scattering coefficient of each component are determined in consideration of the optical effect of the base glaze, which is not colorless and transparent, the computer is used. High-precision prediction can be performed by color matching. In addition, it is only necessary to prepare a sample as shown in Figs. 4, 7 and 8 on which a mixture is applied on a normal substrate, so that it is easy to prepare a sample.

なお、上記第 1実施例に関しては、例えば次のような変形も可能である。 Note that, for the first embodiment, for example, the following modifications are possible.

( 1) 散乱係数 S_w' の値を図 2 (B) に示される 2本の直線 L 1， L 2を数式 6， 7で近似する代わりに、数式 5をそのまま用いて散乱係数 S_w' を決定してもよい。この場合には、図 6のステップ S 2 3においてベース釉の散乱係数 S_B の値を仮決定し、ステップ S 2 3〜S 26を繰り返し実行することによって散乱係数 S_B の値を求めるようにすればよい。こうすれば、未知数が散乱係数 S_B の(1) scattering coefficient S _w 'linear L 1 2 pieces of indicated values in FIG. 2 (B) of the, L 2 instead be approximated by Equation 6, 7, the scattering coefficient using Equation 5 as S _w' May be determined. In this case, temporarily determines the values of the scattering coefficient S _B of the base glaze in step S 2 3 6, to determine the value of the scattering coefficient S _B by repeatedly executing step S 2 3~S 26 do it. In this way, the unknown becomes the scattering coefficient S _B

1つだけなので、ステップ S 2 3において散乱係数 S_B を修正するための逐次近似法の計算の収束が早くなるという利点がある。 Since there is only one, there is an advantage that the convergence of the calculation of the successive approximation method for correcting the scattering coefficient S _B in step S 23 becomes faster.

ただし、ステップ S 2 3では、逐次近似法によって係数を求める必要はなく、種々の係数の値を仮定して、ステップ S 26におけるシミュレーション結果が最適となる係数を選択するようにしてもよレ、。例えば、数式 5をそのまま用いる場合には、 S_s の値を約 0. 00 1〜約0. 0 1 0の範囲で数点設定しておき、ステツプ S 26において実調合率に最も近い調合率を予測できる S_B の値を決定するようにしてもよい。 However, in step S23, it is not necessary to find the coefficients by the successive approximation method, and it is possible to assume various coefficient values and select the coefficient that optimizes the simulation result in step S26. ,. For example, when using Equation 5 as it is, set several values of S _{s in} the range of about 0.001 to about 0.010, and in Step S26, the value closest to the actual mixing rate is set. the compounding ratio may be so that to determine the value of S _B predictable.

散乱係数 S を表現するための関数は数式 5〜 7に限られるわけではなく、一般に、散乱係数 S_w' を調合率 C_w に依存した関数 f (C_w ) で表現するようにすればよい。 The function for expressing the scattering coefficient S is not limited to Equations 5 to 7.In general, the scattering coefficient S _w ′ may be expressed by a function f (C _w ) that depends on the mixing ratio C _w .

なお、上記第 1実施例では、着色剤の調合率を外添加調合率（被着色物の調合率 S_B を 1 00 %とした時の着色剤の調合率）で定義したが、調合率を他の規定の仕方で定義することも可能である。従って、数式 5は、定数 _3l を含む次の数式 23によって表現することができる。 In the above-described first embodiment, defining the formulation ratio of the colorant in the outer additive formulation ratio (compounding ratio of the colorant when the 1 100% The formulation ratio S _B of the coloring material), the preparation rate It can be defined in other prescribed ways. Therefore, Equation 5 can be expressed by the following Equation 23 including the constant _3l .

Sw' = ^i¾ ...(23) Sw '= ^ i¾ ... (23)

し、 ν + ί Then ν + ί

(2) ΧΥΖ表色系や L*a*b* 表色系の値を一致させるコンピュータカラーマッチングは、一般にメタメリックマッチ法と呼ばれている。一方、混合物の分光反射率 R ( λ ) の曲線を一致させるァイソメリックマッチ法と呼ばれる方法もある。この発明は、メタメリックマッチ法のみでなく、ァイソメリックマッチ法によるコンピュータカラーマッチングにも適用可能である。 (2) コンピュータ Computer color matching that matches the values of the color system and the L * a * b * color system is generally called the metameric matching method. On the other hand, There is also a method called an isomeric matching method for matching the curves of the light reflectance R (λ). The present invention is applicable not only to the metameric matching method but also to computer color matching by the isometric matching method.

( 3 ) 上記第〗実施例では陶磁器の釉を对象とするコンピュータカラ一マツチングについて説明したが、本発明はこれに限らず、他の種類の混合物を対象とするコンピュータカラーマッチングにも適用することが可能である。但し、陶磁器の素地は、完全な白色や完全な黒色のものを作成することは困難なので、絶対法によって被着色物（ベース釉）単独の物性値を決定することは難しい。従って、釉に関するコンピュータカラ一マッチングに本発明を適用すれば、予測精度を向上させる上で特に効果が大きい。 (3) In the above-described first embodiment, computer color matching using ceramic glaze has been described. However, the present invention is not limited to this, and is also applicable to computer color matching for other types of mixtures. It is possible. However, since it is difficult to produce a completely white or completely black ceramic body, it is difficult to determine the physical properties of the object to be colored (base glaze) by the absolute method. Therefore, applying the present invention to computer color matching for glaze is particularly effective in improving prediction accuracy.

以下では，第 1実施例で求められた散乱係数と吸収係数とを用いたコンビユータカラーマッチングの種々の適用例としての各種の実施例を説明する。但し，以下の各実施例は，第 1実施例で決定された散乱係数と吸収係数とを用いる場合に限られるわけではなく，他の方法で決定された散乱係数と吸収係数とを用いることも可能であることに注意すべきである。 In the following, various embodiments will be described as various application examples of the computer color matching using the scattering coefficient and the absorption coefficient obtained in the first embodiment. However, each of the following embodiments is not limited to the case where the scattering coefficient and the absorption coefficient determined in the first embodiment are used, and the scattering coefficient and the absorption coefficient determined by other methods are used. It should be noted that both are possible.

Β . 第 2実施例： Β. Second embodiment:

一般に、コンピュータカラ一マッチングによる予測には誤差があるので、その予測誤差を小さくする工夫がである。従来は、予測誤差を減少させるために、各着色剤の吸収係数 Ki と反射係数 S i とを正確な値に近づけるように補正していた。 In general, there is an error in the prediction by computer color matching, so a method to reduce the prediction error is proposed. Conventionally, in order to reduce the prediction error, the absorption coefficient Ki and the reflection coefficient S i of each colorant have been corrected so as to be close to accurate values.

し力し、各成分の吸収係数と反射係数 S i とを正確な値に近づけるためには、被着色物に各着色剤を単独で舍む数多くの混合物サンプルを作成して、その分光反射率を測定しなければならず、膨大な作業を要していた。また、天然の顔料を着色剤として用いる場合には、その吸収係数 κ_± や散乱係数 S i が必ずしも一定の値にはならないので、吸収係数 Ki と散乱係数 S i を正確な値に近づけることが困難であった。 In order to make the absorption coefficient and the reflection coefficient S i of each component close to the exact values, a large number of mixture samples were prepared, in which each colorant was added to the object to be colored, and the spectral reflection The rate had to be measured and required a lot of work. When a natural pigment is used as a coloring agent, the absorption coefficient κ _± and the scattering coefficient S i do not always become constant, so that the absorption coefficient Ki and the scattering coefficient S i are close to accurate values. It was difficult to find out.

以下に示す第 2実施例は、各成分の吸収係数 Ki と散乱係数 S i とを補正することなく、予測誤差を減少させることを目的としたものである。 The second embodiment described below corrects the absorption coefficient Ki and the scattering coefficient S i of each component. The aim is to reduce the prediction error without any change.

図】 1は、第 2実施例における処理の全体手順を示すフローチャートである。なお、この第 2実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉（ゆう）である。すなわち、顔料を入れないベース釉（基礎釉）が被着色物であり、このべ一ス釉に顔料を添加した釉がコンピュータカラ一マッチングの対象となる混合物である。 FIG. 1 is a flowchart illustrating an overall procedure of a process according to a second embodiment. The mixture targeted in the second embodiment is a glaze for covering the surface of the ceramic body. In other words, the base glaze (base glaze) without pigment is the object to be colored, and the glaze obtained by adding pigment to this base glaze is the mixture to be subjected to computer color matching.

ステップ S 3 1では、コンピュータカラ一マッチングの予測結果（三刺激値）を補正するためのニューラルネットワークの学習を行なう。ステップ S 3 2〜S 3 5では、学習済みのニューラルネットワークを用いてコンピュータカラーマツチングの目標値を補正して、正確な予測結果を求めている。以下ではまず、ステツプ S 3 2〜S 3 5の内容を説明する前に、ニューラルネットワークの構成とステツプ S 3 1の詳細手順について説明する。 In step S31, learning of a neural network for correcting the prediction result (tristimulus value) of computer color matching is performed. In steps S32 to S35, the target value of computer color matching is corrected using the trained neural network to obtain an accurate prediction result. Hereinafter, before describing the contents of steps S32 to S35, the configuration of the neural network and the detailed procedure of step S31 will be described first.

図 1 2は、ニューラルネットワークの構成を示す説明図である。この二ユーラルネットワークは、入力層 1 0と中閒肩 2 0と出力層 3 0とで構成される三層の階層構造を有している。人力層 I 0は、 3つのニューロン N 11〜N 13で構成されており、中間層 2 0は 5つのニューロン N 21〜N 25で、出力肩 3 0は 3つのニューロン N31〜N33でそれぞれ構成されている。 FIG. 12 is an explanatory diagram showing the configuration of the neural network. This dual-purpose network has a three-layer hierarchical structure composed of an input layer 10, a middle shoulder 20 and an output layer 30. The human stratum I0 is composed of three neurons N11 to N13, the middle layer 20 is composed of five neurons N21 to N25, and the output shoulder 30 is composed of three neurons N31 to N33. Have been.

入力層 1 0の 3つのニューロン N 11〜N 13には、二:刺激値 X, Υ , Zがそれぞれ入力される。入力層 1 0のニューロン N ijから中間層 2 0のニューロン Nk に伝達される信号は、それぞれの入力信号に重み W , _k を乗じたものである。ここで、 iは注目している階屑を示す番号、 jは注目している階層内でのニューロンの順番を示す番号、 kは次の階層のニューロンの順番を示す番号である。例えば、入力餍 1 0の第 1のニューロン N ilから中聞層 2 0の第 1 のニューロン N 21に伝達される信号は W_u. であり、入力屑 1 0の第 1のニューロン N ilから中間層 2 0の第 2のニューロン N 22に伝達される信号は W_u, ₂Xである。 Two: stimulus values X, ,, and Z are input to the three neurons N11 to N13 of the input layer 10 respectively. The signal transmitted from the neuron Nij of the input layer 10 to the neuron Nk of the hidden layer 20 is _obtained by multiplying each input signal by weights W and _k . Here, i is a number indicating the debris of interest, j is a number indicating the order of neurons in the layer of interest, and k is a number indicating the order of neurons in the next layer. For example, the signal transmitted from the first neuron N il of the input 餍 10 to the first neuron N 21 of the middle layer 20 is W _u ., And the first neuron N il of the input dust 10 is signal transmitted to the second neuron N 22 of pressurized et intermediate layer 2 0 is W _u, ₂ X.

中聞層 2 0の各ニューロン N ijの入力と出力 Q _Ajとの関係は、次の数式 2 4に示す情報伝達関数 f ( _{U ij}) で与えられる。

The relationship between the input and output Q _Aj of each neuron N ij in the middle hearing layer 20 is given by an information transfer function f ( _{U ij} ) shown in the following Expression 24.

ここで、は（ i 一〗）番目の階層（すなわち入力層 1 0) の〗番目のニューロン N (i- l)jの出力であり、闵 1 2の例では Q_U=X， 0₁₂=Y, 0₁₃ = Ζである。また、

_k は（ i— 1) 番目の階層の j番目のニューロン N (i-1 ) jから、注目しているニューロン Nijに伝達される信号に掛かる重みである。 tはしきい値であり、一定の他が割当てられる。なお、数式 24の情報伝達関数 f (u , はシグモイド関数と呼ばれている。 Here, is the output of the〗 th neuron N (i-l) j in the (i〗 1) th hierarchy (that is, the input layer 10). In the example of 闵₁₂ , Q _U = X, 0 ₁₂ = Y, 0 ₁₃ = Ζ. Also,

_k is the weight applied to the signal transmitted from the j-th neuron N (i-1) j in the (i-1) th hierarchy to the neuron Nij of interest. t is a threshold, some other being assigned. Note that the information transfer function f (u, in Equation 24 is called a sigmoid function.

例えば、数式 24を中閒層 20の第〗のニューロン N21の入出力関係に適用すると、次の数式 25が得られる。

u₂i = W_n,iX + W,2jY + W₁₃₁Z -..(25b) 出力層 30の各ニューロンの入出力関係も上記数式 24で与えられる。図 1 2 に示す第 2実施例では、出力層 30の 3つのニューロン N31〜N33の出力For example, when Equation 24 is applied to the input / output relationship of the second neuron N21 of the middle layer 20, the following Equation 25 is obtained.

u ₂ i = W _n , iX + W, 2jY + W ₁₃₁ Z-.. (25b) The input / output relationship of each neuron of the output layer 30 is also given by the above equation (24). In the second embodiment shown in FIG. 12, the output of the three neurons N31 to N33 in the output layer 30 is

〜Q₁₃を、コンピュータカラ一マッチング（CCM) による三刺激値の予測誤差 ΔΧ, ΔΥ, ΔΖとしている。 The to Q _13, predicts errors ΔΧ tristimulus values by computer color one matching (CCM), ΔΥ, and Derutazeta.

ニューラルネットワークの学習は、入力屑 1 0への入力（X, Y, Z) と出力層 30からの出力（ΔΧ, ΔΥ, ΔΖ) との間の関係を数多く与えて、正しい入出力関係を与えるような重み W , _k の値を決定する作業である。 Neural network learning gives a number of relationships between the inputs (X, Y, Z) to the input debris 10 and the outputs (ΔΧ, ΔΥ, ΔΖ) from the output layer 30 to give the correct input / output relationships. This is the task of determining the values of the weights W and _k .

図 1 3は、図 1 1のステップ S 3 1の詳細手順を示すフローチャートである。ステップ S 4 1では、コンピュータカラ一マッチングにおいて予測の対象とする混合物の色の範囲をカバーするような複数の三刺激値（Xit. Yit, Zit) 定めて、これらの複数の三刺激値を有するような複数のサンプルの調合割合をコンピュータカラーマッチングで決定する。図 1 4は、コンピュータカラーマッチングでの色の予測対象範囲 P Aと、この予測対象範囲 P Aをカバーする複数の三刺激値の分布を示す概念図である。この第 2実施例では C I E— XYZ表色系で色を表現するものとしており、色の予測対象範囲 P Aは X Y Z座標系の 3次元的な範囲として^えられる。なお、予測対象範囲 P Aは、予測対象とする混合物が取り得る色の範囲を示すものであり、任意に設定し得る範囲である。 FIG. 13 is a flowchart showing the detailed procedure of step S31 in FIG. In step S41, a plurality of tristimulus values (Xit. Yit, Zit) that cover the color range of the mixture to be predicted in computer color matching Then, the blending ratio of a plurality of samples having the plurality of tristimulus values is determined by computer color matching. FIG. 14 is a conceptual diagram illustrating a color prediction target range PA in computer color matching and a distribution of a plurality of tristimulus values covering the prediction target range PA. In the second embodiment, colors are expressed in the CIE-XYZ color system, and the color prediction target range PA can be obtained as a three-dimensional range in the XYZ coordinate system. Note that the prediction target range PA indicates the range of colors that the mixture to be predicted can take, and is a range that can be set arbitrarily.

この第 2実施例においては、予測対象範囲 P Aをカバーするために、図〗 4に白丸で示す 7組の三刺激値を決定した。ステップ S 4 1では、さらに、これらの 7組の三刺激値（Xit, Y it, Z it) を有するような 7種類のサンプル M 1〜 M 7の調合割合を、コンピュータカラーマッチングによって予測した。 In the second example, seven sets of tristimulus values indicated by white circles in FIG. 4 were determined in order to cover the prediction target range PA. In step S41, the compounding ratio of the seven types of samples M1 to M7 having these seven sets of tristimulus values (Xit, Yit, Zit) was further predicted by computer color matching.

ここで、コンピュータカラーマッチングによる色予測と調合割合の予測について簡単に説明する。第 1実施例で説明した数式 2を変形すると、混合物の分光反射率 R (ス）は、次の数式 2 6で与えられる。 Here, the color prediction by computer color matching and the prediction of the mixing ratio will be briefly described. By transforming Equation 2 described in the first embodiment, the spectral reflectance R (s) of the mixture is given by Equation 26 below.

R= 1 + fK] ffi IK} ...(26) R = 1 + fK] ffi IK} ... (26)

M 混合物の吸収係数と散乱係数の比（KZS) _κは、各成分の吸収係数 Ki (1) と散乱係数 S_± ）と調合率 Ci から，第 1実施例で説明した数式 1に従って算出できるので、混合物の分光反射率 R (λ ) を上記数式 2 6から求めることができる。この分光反射率 R (λ ) は理想状態（被着色物の厚みが無限大の場合）の分光反射率なので、分光光度計で測定できる分光反射率 R' (； -) を、次の数式 2 7 (サンダーソンの式）に従って求める。 The ratio (KZS) _κ between the absorption coefficient and the scattering coefficient of the M mixture can be calculated from the absorption coefficient Ki (1) and the scattering coefficient S _± ) of each component and the mixing ratio Ci according to Equation 1 described in the first embodiment. The spectral reflectance R (λ) of the mixture can be obtained from the above equation 26. Since this spectral reflectance R (λ) is the spectral reflectance in an ideal state (when the thickness of the object to be colored is infinite), the spectral reflectance R '(;-) that can be measured by a spectrophotometer is expressed by the following equation. 2 7 (Sanderson's equation)

R' = k₁-Kl -ki)(l -k₂)_l ..- (27) ここで、係数， k ₂ は、被着色物（ベース釉）の光学的性質に依存する値である。係数， k ₂ としては、被着色物の屈折率 nから次の数式 2 8に従って決定することができる。

R '= k ₁ -Kl -ki) (l -k ₂ ) _l ..- (27) where the coefficient, k ₂ is a value that depends on the optical properties of the object to be colored (base glaze). The coefficients, k ₂ can be determined from the refractive index n of the object to be colored according to the following equation 28.

k₂=0.68n-0.56 . - . (28b) k ₂ = 0.68n-0.56 .-. (28b)

なお、第 2実施例において用いたベース釉は、屈折率 nが約 1. 4である。数式 27によって分光反射率 R' ）が求まると、混^の三刺激値 X， Y, Zが次の数式 29によって求められる。 The base glaze used in the second example had a refractive index n of about 1.4. When the spectral reflectance R ') is determined by Equation 27, the tristimulus values X, Y, and Z of the mixture are determined by Equation 29 below.

X = k]S( )R'( )x( )dX ...(29a) X = k] S () R '() x () dX ... (29a)

Y = kj S( )R'( )y( )dX - - . (29b) Y = kj S () R '() y () dX--. (29b)

Z = k)S(X)R )z^)dX ... (29c) k =,—— ... (29d) Z = k) S (X) R) z ^) dX ... (29c) k =, ---- ... (29d)

ここで、 S (λ ) は標準光の分光分布、 X (/) , y (λ ) , ζ (λ ) (数式中ではバー付きである）は等色関数である。 Here, S (λ) is the spectral distribution of the standard light, and X (/), y (λ), and ζ (λ) (with a bar in the formula) are color matching functions.

コンピュータカラ一マッチングによって混合物の色を予測する場合には、上述のように、混合物の各成分の吸収係数と散乱係数に基づき、数式 26〜29 に従って、その混合物の三刺激値 X, Υ, Ζを算出する。三刺激値は混合物の色を表わすので、任意の混合物の色を予測することができることになる。 When predicting the color of a mixture by computer color matching, as described above, based on the absorption coefficient and scattering coefficient of each component of the mixture, the tristimulus values X, Υ, Ζ Is calculated. Since the tristimulus values represent the color of the mixture, the color of any mixture can be predicted.

また、所望の色を有する混合物の調合割合を予測する場合には、混合物の調合割合を仮定して、上述した手順でその三刺激値を算出し、ニュートン一ラプソン法などの逐次近似法によって所望の色に所定の誤差内で一致するような調合割 1 を求める。 Further, when predicting the blending ratio of a mixture having a desired color, the tristimulus value is calculated by the above-described procedure assuming the blending ratio of the mixture, and the desired value is calculated by a successive approximation method such as the Newton-Raphson method. Formula 1 is determined so that the color matches the color within a specified error.

図 1 3のステップ S 42では、ステップ S 4 〗で予測された調合割合 C iを有する複数のサンプルを作成する。この第 2実施例では、図 1 4の 7組の三刺激値 (Xit, Yit, Zit) に対応する 7つのサンプル M 1〜M 7を作成した。ステツプ S 43では、各サンプル M iの分光反射率を分光光度計で測定し、上記数式 2 9に従ってその三刺激値（Xim Yim, Zin を求める。ステップ S 44では、各サンプル M iに関して、ステップ S 43で得られた実測値（Xim， Yim, Zim) とステップ S 4 1で決定した目標値（Xit, Yit, Zit) との差をとることによって、予測誤差 AMi (Xim-Xit, Yim-Yit- ， Zim— Zit) を求める。図 1 5は、第 2実施例において得られた 7つのサンプル M 1〜M 7の予測誤差を示す概念図である。また、図 1 6は、各サンプルの三刺激値の目標値（ステップ S 4 1で決定された値）と予測誤差 ΔΜΪ (ΔΧ, ΔΥ, ΔΖ) を示す説明「¾1である。 In step S42 of FIG. 13, a plurality of samples having the compounding ratio C i predicted in step S4〗 are created. In the second embodiment, seven samples M1 to M7 corresponding to the seven sets of tristimulus values (Xit, Yit, Zit) in FIG. 14 were created. In step S43, the spectral reflectance of each sample Mi is measured with a spectrophotometer, and its tristimulus values (Xim Yim, Zin) are obtained according to the above equation (29). In step S44, for each sample Mi, the difference between the measured value (Xim, Yim, Zim) obtained in step S43 and the target value (Xit, Yit, Zit) determined in step S41 is calculated. By doing so, the prediction error AMi (Xim-Xit, Yim-Yit-, Zim—Zit) is obtained. FIG. 15 is a conceptual diagram illustrating prediction errors of seven samples M1 to M7 obtained in the second example. FIG. 16 is a description “説明 1” showing the target value of the tristimulus value (the value determined in step S41) of each sample and the prediction error ΔΜΪ (ΔΧ, ΔΥ, ΔΖ).

なお、ステップ S 4 1で得られた調合割合 C iからコンピュータカラ一マッチングによって各サンプルの三刺激俩の予測値 (Xic, Yic, Zic) を求め、実測値 (Xim, Yim, Zim) と予測愤 (Xic, Yic, Zic) との差（Xim—X ic, Y im-Y ic, Zim— Zic) を予測誤差 AMiと定義してもよい。コンピュータカラーマッチングで調合割合 C iを決定する際には、三刺激値の冃標値 The predicted values (Xic, Yic, Zic) of the tristimulus の of each sample were obtained from the mixture ratio C i obtained in step S41 by computer color matching, and the measured values (Xim, Yim, Zim) were obtained. The difference (Xim—Xic, Yim-Yic, Zim—Zic) between the prediction error and the prediction 愤 (Xic, Yic, Zic) may be defined as the prediction error AMi. When determining the mixing ratio C i by computer color matching, the target value of tristimulus values

(Xit, Yit, Zit) と予測値 (Xic, Yic, Zic) との差が所定の許容誤差以下になるように調合割合 C iを決定するので、目標値（Xit, Yit, Zit- ) と予測値（Xic， Yic, Zic) は実質的にほぼ等しい値を有している。従つて、実測値 (Xim, Yim, Zim) と予測愤 (Xic, Yic, Zic) との差（X im-Xic, Yim— Yic, Zim— Zic) を予測誤差 ΔΜίと定義しても、実測値Since the mixing ratio C i is determined so that the difference between (Xit, Yit, Zit) and the predicted value (Xic, Yic, Zic) is equal to or less than a predetermined allowable error, the target value (Xit, Yit, Zit-) And the predicted values (Xic, Yic, Zic) have substantially the same value. Therefore, even if the difference (X im-Xic, Yim—Yic, Zim—Zic) between the measured value (Xim, Yim, Zim) and the prediction 愤 (Xic, Yic, Zic) is defined as the prediction error ΔΜί, Measured value

(Xim, Yim, Zim) と目標値（Xit， Yit, Zit) との差（Xim-Xit, Yim- Yit, Zim- Zit) を予測誤差 AMiと定義しても実質的にはほぼ同じである。 Even if the difference (Xim-Xit, Yim-Yit, Zim-Zit) between (Xim, Yim, Zim) and the target value (Xit, Yit, Zit) is defined as the prediction error AMi, it is virtually the same. .

図 1 3のステップ S 45では、図 1 6に示す各サンプル M iの三刺激値の目標値（Xit, Yit, Zit) と予測誤差 ΔΜΐとを用いてニューラルネットワークの学習を行ない、上記数式 24における重み Wij， kを決定する。ニューラルネットワークの学習方法としては、例えば逆誤差伝搬学習方式を用いる。 In step S45 in FIG. 13, learning of the neural network is performed using the target values (Xit, Yit, Zit) and the prediction error Δ サンプル of the tristimulus values of each sample Mi shown in FIG. Determine the weights Wij, k at. As a learning method of the neural network, for example, a back error propagation learning method is used.

図 1 7は、第 2実施例におけるニューラルネットワークの学習の実証結果を示す説明図である。ここでは、図 1 4， 5に示す予測対象範囲 P Aに含まれるもう 1組の三刺激値をコンピュータカラ一マッチングの目標値として設定し、この目標値を有する第 8のサンプル M 8を作成した。そして、このサンプル M 8の三刺激俯を実測した。図 I 7において「CCM目標値」とあるのは第 8のサンプル M 8のコンピュータカラーマッチングに用いた目標値を意味している。また、「誤差（真植）」とあるのは三刺激値の目標値と実測値との差である。「ニューロ予測誤差」は、 CCM目標値を学習済みのニューラルネットワーク（図 1 2) に入力した場合に得られる予測誤差である。図〗 7の結果から、学習済みの二ユーラルネットワークは、三刺激値の誤差を精度良く予測できることが解る。 FIG. 17 is an explanatory diagram showing the results of the verification of the neural network learning in the second embodiment. Here, another set of tristimulus values included in the prediction target range PA shown in Figs. 14 and 5 is set as the target value for computer color matching. An eighth sample M8 having a standard value was prepared. Then, the sagittal depression of this sample M8 was actually measured. In FIG. I7, “CCM target value” means the target value used for the computer color matching of the eighth sample M8. In addition, "error (true planting)" is the difference between the target value and the measured value of the tristimulus value. “Neuro prediction error” is the prediction error obtained when the CCM target value is input to a trained neural network (Fig. 12). From the results in Fig. 7, it can be seen that the trained bi-ural network can accurately predict the error of tristimulus values.

こうしてニューラルネットワークの学習が終了すると、図 1 1のステップ S 3 2〜S 35を実行して、調合率未知の色見本の調合率を予測する。ステップ S 3 2では、調合率未知の色見本の分光反射率を測定してその三刺激値 (X s, Y s, Z s ) を求める。ステップ S 3 3では、図 1 2に示すニューラルネットワークに色見本の三刺激値（X s， Y s, Z s) を入力して、予測誤差△ (AX s， AY s， ΔΖ s) を求める。テツプ S 34では、三刺激値 (X s, Y s, Z s ) を予測誤差厶で補正して、コンピュータカラーマッチングの調合割合予測に用いる三刺激値の目標储（X s— ΔΧ s， Y s -ΔΥ s , Zs—AZs) を求める。ステツプ S 35では、こうして補正された目標値を用いてコンピュータカラ一マツチングを実行し、色見本の調合率を予測する。 When learning of the neural network is completed in this way, steps S32 to S35 in FIG. 11 are executed to predict the mixing ratio of a color sample whose mixing ratio is unknown. In step S32, the spectral reflectance of the color sample whose mixing ratio is unknown is measured, and its tristimulus value (Xs, Ys, Zs) is obtained. In step S33, the tristimulus values (Xs, Ys, Zs) of the color samples are input to the neural network shown in Fig. 12, and the prediction error △ (AXs, AYs, ΔΖs) is obtained. . In step S34, the tristimulus values (Xs, Ys, Zs) are corrected by the prediction error, and the target of the tristimulus values 储 (Xs—ΔΧs, Y s -ΔΥ s, Zs—AZs). In step S35, computer color matching is performed using the corrected target value to predict the mixing ratio of the color sample.

図 1 8は、第 2実施例におけるコンピュータカラーマッチングの予測精度を検証するために行なった実験結果を示す説明図である。ここでは、予測精度を実証することを目的としたので、調合率既知の色見本について三刺激値を実測し、その実測値（X s， Y s, Z s) を実現する調合率をコンピュータカラーマツチングで予測した。図 1 7の結果から、調合率の真値とコンピュータカラーマツチングによる予測値とは極めて良く一致していることが解る。 FIG. 18 is an explanatory diagram showing the results of an experiment performed to verify the prediction accuracy of computer color matching in the second embodiment. Here, we aimed to verify the prediction accuracy, so we measured the tristimulus values for a color sample with a known blending ratio, and calculated the blending ratio that realized the measured values (Xs, Ys, Zs) using a computer. Predicted by color matching. From the results in Fig. 17, it can be seen that the true value of the blending ratio and the predicted value by computer color matching are in very good agreement.

なお、図 1 2に示す学習済みのニューラルネットワークを利用すれば、調合率既知の混合物の三刺激値も精度良く予測することが可能である。すなわち、予測対象の混合物の調合率 C iから上記数式】， 2 6〜 29を用いて三刺激値を求め、この三撒値を図 1 2のニューラルネッ卜ワークに入力して予測誤差△を求める。そして、数式 29で得られた三刺激値を予測誤差厶で補正すれば、実際の値に極めて近い :刺激値が得られる。 If the learned neural network shown in FIG. 12 is used, it is possible to accurately predict the tristimulus values of the mixture with a known mixing ratio. That is, tristimulus values are obtained from the mixture ratio C i of the mixture to be predicted using the above formulas, 26 to 29, and the tristimulus values are input to the neural network shown in FIG. Ask. Then, if the tristimulus value obtained by Equation 29 is corrected by the prediction error, it will be extremely close to the actual value. Very close: Stimulus value is obtained.

プラスチック等と異なり，陶磁器においては，焼成条件，熔融中の原料の化学反応等により，混色理論が成立しない場合がある。この場合，ニューラルネットワークに混合物の三刺激値と C C Mによる予測誤差の関係を学習させる本実施例の方法は有効である。 Unlike plastics, etc., in ceramics, the theory of color mixing may not be established due to firing conditions, chemical reaction of raw materials during melting, and the like. In this case, the method of this embodiment in which the neural network learns the relationship between the tristimulus value of the mixture and the prediction error due to CCM is effective.

なお、上記第 2実施例に関しては、例えば次のような変形も可能である。 Note that, for the second embodiment, for example, the following modifications are also possible.

( 1 ) 上記第 2実施例ではコンピュータカラーマッチングにおける三刺激値の目標値をニューラルネットワークの入力としていたが、ニューラルネットワークの人力としては、各サンプル M iの三刺激値の実測値（Xim， Yim, Zim) を用いても良く、また、コンピュータカラーマッチングによる予測値（Xic， Yic- ， Zic) を用いても良い。すなわち、ニューラルネットワークには、複数のサンプルの三刺激値（XYZ表色系の座標値）とその予測誤差との関係を学習させるようにすればよい。 (1) In the second embodiment, the target values of the tristimulus values in the computer color matching were input to the neural network. However, as the human power of the neural network, the measured tristimulus values (Xim, Yim, Zim) may be used, or predicted values (Xic, Yic-, Zic) by computer color matching may be used. In other words, the neural network may be made to learn the relationship between the tristimulus values (coordinate values in the XYZ color system) of a plurality of samples and their prediction errors.

(2) t記第 2実施例では、ニューラルネットワークを用いて三刺激値を補正していたが、回帰分析やニューロフアジィ技術等の他の誤差補正法を用いて補正することも可能である。 (2) In the second embodiment, the tristimulus values are corrected using a neural network, but may be corrected using other error correction methods such as regression analysis or neuro-fuzzy technology. is there.

(3) 上記第 2実施例では XYZ表色系を用いていたが、本発明は、 L*a*b* 表色系などの他の任意の表色系を用いる場合に適用できる。 (3) Although the XYZ color system is used in the second embodiment, the present invention can be applied to the case where any other color system such as the L * a * b * color system is used.

(4) XYZ表色系や L*a*b* 表色系の座標値を一致させるコンピュータ力ラーマッチングは、一般にメタメリックマッチ法と呼ばれている。一方、混合物の分光反射率 R (λ ) の曲線を一致させるァイソメリックマッチ法と呼ばれる方法もある。この発明は、メタメリックマッチ法のみでなく、ァイソメリックマツチ法によるコンピュータカラーマッチングにも適用可能である。 (4) Computer color matching that matches the coordinate values of the XYZ color system and the L * a * b * color system is generally called the metameric matching method. On the other hand, there is also a method called an isomeric matching method in which a curve of the spectral reflectance R (λ) of a mixture is matched. The present invention is applicable not only to the metameric matching method but also to the computer color matching by the isometric matching method.

( 5 ) 上記第 2実施例では陶磁器の釉を対象とするコンピュータカラーマツチングについて説明したが、本発明はこれに限らず、他の種類の混合物を対象とするコンピュータカラーマッチングにも適用することが可能である。 C. 第 3実施例： (5) In the second embodiment, computer color matching for ceramic glaze was described, but the present invention is not limited to this, and is also applicable to computer color matching for other types of mixtures. It is possible. C. Third Embodiment:

衛生陶器などのような工業陶器には、ユーザの様々な好みに応じて種々の色が設定されている。陶器の色は、釉に混合する顔料の調合率によって決定されるので、釉は、コンピュータカラーマッチングにおける混合物に相当する。陶器の色を設定する際には、まず、設計者が紙などを彩色したり、予め焼成されたサンプルから所望の色を有するサンプルを選択したりすることによって色を決定する。そして、その色の分光反射率を分光光度計で測定し、分光反射率の実測値からコンピュータマツチングを用いて顔料や染料等の着色剤の調合割合を予測する。ところで、同じ調合割合の釉薬を用いても、製造される陶器の実際の色は同一にはならず、色にばらつきがあるのが普通である。図〗 9は、設計された色と実際に製造される陶器の色のばらつきを示する X— y色度図である。図 1 9において、二重丸で示される色が設計された所望の色（標準色） L 0であり、小さな白丸が実際に製造される陶器の色の分布である。標準色 L 0から大きく異なる色を有する陶器は不良品として認識する必要がある。そこで、良品の色の濃淡の限度を示すものとして、図 1 9に黒丸で示す色を有する濃淡限度色 L 1， L 2が設定される。標準色 L 0の色度座標値は設計者によって決定されているので、この標準色 L 0の色度座標値に製造誤差等を考慮して、 2つの'濃淡限度色 L 1， L 2の色度座標値を決定する。 Industrial porcelain, such as sanitary ware, has various colors set according to various user preferences. Since the color of the pottery is determined by the blending ratio of the pigment mixed with the glaze, the glaze corresponds to the mixture in computer color matching. When setting the color of the pottery, the designer first determines the color by coloring the paper or the like or selecting a sample having a desired color from pre-baked samples. Then, the spectral reflectance of the color is measured with a spectrophotometer, and the mixing ratio of a coloring agent such as a pigment or a dye is predicted from the measured value of the spectral reflectance using computer matching. By the way, even if glazes with the same blending ratio are used, the actual colors of the pottery produced will not be the same, and the colors will usually vary. Figure 9 is an X-Y chromaticity diagram showing the variation between the designed colors and the colors of the pottery actually manufactured. In Fig. 19, the color indicated by the double circle is the designed color (standard color) L0, and the small white circle is the color distribution of the ceramics actually manufactured. Pottery with a color that is significantly different from the standard color L0 must be recognized as defective. Therefore, the shade limit colors L1 and L2, which have the colors indicated by the black circles in Fig. 19, are set to indicate the shade limits of non-defective products. Since the chromaticity coordinate values of the standard color L0 are determined by the designer, the chromaticity coordinate values of the standard color L0 are taken into consideration by considering the manufacturing error and the like, and the two 'shade limit colors L1, L2' Determine the chromaticity coordinate value of.

実際の陶器の製造においては、標準色 L 0を有する陶器を標準色見本として作成し、また、濃淡限度色 L I , L 2を有する 2つの陶器を濃淡限度見本として作成する。そして、検査工程において、標準色見本および渙淡限度見本の色と製造された陶器の色とを比較して、渙淡限度見本の色の範囲にある陶器のみを良品とする。 In the actual production of pottery, a pottery having the standard color L0 is created as a standard color sample, and two potteries having shading limit colors L I and L2 are created as shading limit samples. Then, in the inspection process, the color of the standard color sample and the color limit sample is compared with the color of the manufactured ceramic, and only the ceramics in the color range of the color limit sample are regarded as good products.

ところで、コンピュータカラ一マッチングによる予測にはかなりの誤差が伴うので、予測された調合割合から標準色見本ゃ澹淡限度見本を作成するにはかなりの手間を要するのが普通である。例えば、標準色見本を作成する際には、コンピュータカラーマッチングで予測された調合割合で調合した釉薬を用いて多数のサンプルを焼成し、その中で所望の標準色を有するサンプルがあればそれを標準色見本として採用する。しかし、標準色を有するサンプルが無ければ、勘と経験に従って調合割合を変更して再度多数のサンプルを焼成し、標準色を有するサンプルが得られるまでこの作業を繰り返す。このような作業は、濃淡限度色 L 1, L 2を有する濃淡限度見本についても行なわれる。 By the way, since the prediction by the computer color matching involves a considerable error, it usually takes a considerable amount of time to create a standard color sample and a shading limit sample from the predicted mixing ratio. For example, when preparing a standard color sample, a large number of samples are fired using glaze prepared at the mixing ratio predicted by computer color matching, and if there is a sample having the desired standard color, The standard Adopt as a color sample. However, if there is no sample having the standard color, the mixing ratio is changed according to intuition and experience, and a large number of samples are baked again, and this operation is repeated until a sample having the standard color is obtained. Such an operation is also performed on the gray scale sample having the gray scale colors L1 and L2.

このように、従来のコンピュータカラーマッチングでは、所望の色を有する混合物を作成するための調合割合を精度良く予測するのが困難であった。 As described above, in the conventional computer color matching, it has been difficult to accurately predict the mixing ratio for producing a mixture having a desired color.

以下に示す第 3実施例は、所望の色を有する混合物の調合割合を精度良く予測することを目的としたものである。 The third embodiment described below aims at accurately estimating the blending ratio of a mixture having a desired color.

図 2 0は、第 3実施例における処理の全体手順を示すフローチャートである。なお、この第 3実施例で対象とする混合物は、陶磁器の素地の表面を覆うための釉（ゆう）である。すなわち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラーマッチングの対象となる混合物である。 FIG. 20 is a flowchart showing the overall procedure of the process in the third embodiment. The mixture targeted in the third embodiment is a glaze for covering the surface of a ceramic body. In other words, the base glaze without pigment (base glaze) is the object to be colored, and the glaze obtained by adding pigment to this base glaze is the mixture to be subjected to computer color matching.

ステップ S 5 1では、まず標準色見本を作成する。従って、標準色見本の着色剤の混合割合は既知である。ステップ S 5 2では、作成された標準色見本の分光反射率 R' を分光光度計で測定し、この分光反射率 R' から次の数式 30に従つて三刺激値 X。， Υ。， Ζ₍₎ を算出する。 In step S51, first, a standard color sample is created. Therefore, the mixing ratio of the colorant of the standard color sample is known. In step S52, the spectral reflectance R 'of the created standard color sample is measured with a spectrophotometer, and the tristimulus value X is calculated from the spectral reflectance R' according to the following equation (30). , Υ. , Ζ Calculate ₍₎ .

Xo = klS(X)R*( )x( )dX - .. (30a) Xo = klS (X) R * () x () dX-.. (30a)

Y₀ = kjS(X)R*(X)y(X)dX ... (30b) Zo=k(S( )R'( )z( )dX ... (30c) k = 100 . - . (30d) Y ₀ = kjS (X) R * (X) y (X) dX ... (30b) Zo = k (S () R '() z () dX ... (30c) k = 100 .-. (30d)

ここで、 S (λ ) は標準-光の分光分布、 X (λ) , y ）， ζ ( Ο (数式中ではバー付きである）は等色関数である。 Here, S (λ) is the standard-light spectral distribution, and X (λ), y) and ζ (Ο (with a bar in the formula) are color matching functions.

ステップ S 5 3では、ステップ S 5 2で得られた標準色の三刺激値の実測値 X 。， Υ。， Z。に基づいて、製造誤差を考慮して、設計者が図〗 9に示す 2つの濃淡限度色 L I， L 2に対する三刺激値を設定する。なお、以下では説明の便宜、第〗の濃淡限度色 L〗を実現するための調合割合を予測する場合について説明する。第 2の濃淡限度色 L 2に対しても、同様な処理によってその調合割合を精度良く予測することができる。 In step S53, the measured value X of the tristimulus value of the standard color obtained in step S52. , Υ. , Z. Based on the above, considering the manufacturing error, the designer Set the tristimulus values for the shade limit colors LI, L2 of. In the following, for convenience of description, a case will be described in which a blending ratio for realizing the first shade limit L 限度 is predicted. With respect to the second shade limit color L2, the blending ratio can be accurately predicted by the same processing.

ステップ S 5 4では、ステップ S 5 1で作成した標準色見本の各着色剤の調合率 Ci ( iは着色剤の番号を示す）と、散乱係数 Si と、吸収係数 Ki とから、標準色見本の三刺激値の計算値 X_M , Y_M , Z_N を以下の手順で算出する。まず、第 1実施例で説明した数式 2を変形すると、混合物の分光反射率 R (λ ) は、次の数式 3 1で与えられる。

混合物の吸収係数と散乱係数の比（KZS) „は、各着色剤の吸収係数 (λ ) と散乱係数 Si ( ) と調合率 Ci から、第 1実施例で説明した数式 1に従って算出できる。従って、混合物の分光反射率 R (λ ) は、この比（KZS) _Μ から、数式 3 1に従って求めることができる。この分光反射率 R (/) は理想状態（被着色物の厚みが無限大の場合）の分光反射率なので、分光光度計で測定できる分光反射率 R' (λ ) を、次の数式 32 (サンダーソンの式）に従って求める。 In step S54, the standard color sample is obtained from the mixing ratio Ci (i indicates the number of the colorant), the scattering coefficient Si, and the absorption coefficient Ki of each colorant of the standard color sample created in step S51. The calculated values X _M , Y _M , and Z _N of the tristimulus values are calculated by the following procedure. First, by transforming Equation 2 described in the first embodiment, the spectral reflectance R (λ) of the mixture is given by the following Equation 31.

The ratio (KZS) 散乱 between the absorption coefficient and the scattering coefficient of the mixture can be calculated from the absorption coefficient (λ), the scattering coefficient Si (), and the blending ratio Ci of each colorant according to Equation 1 described in the first embodiment. Therefore, the spectral reflectance R of the mixture (lambda) from the ratio (KZS) _Micromax, can be determined according to equation 3 1. the thickness of the spectral reflectance R (/) is an ideal state (the colorings infinite ), The spectral reflectance R '(λ) that can be measured with a spectrophotometer is calculated according to the following equation 32 (Sanderson's equation).

R' = k₁ + (l-k₁)(l-k₂)₁ ...(32) ここで、係数 , k₂ は、被着色物（ベース釉）の光学的性質に依存する値である。係数 1^ ， k , の値は、被着色物の屈折率 nから次の数式 3 3に従つて決定することができる _c k_{1 =} Ml-J- (33a) R '= k ₁ + (lk ₁ ) (lk ₂ ) ₁ ... (32) where the coefficient, k ₂ is a value that depends on the optical properties of the object to be colored (base glaze). The values of the coefficients 1 ^, k, can be determined from the refractive index n of the object to be colored according to the following equation 33 _: _c k _{1 =} Ml-J- (33a)

\n+ 1 \ n + 1

k =0.68n-0.56 .. - (33b) k = 0.68n-0.56 ..-(33b)

なお、ベース釉の屈折率 nは、例えば約 1. 4である。 The refractive index n of the base glaze is, for example, about 1.4.

数式 3 2によって分光反射率 R' {?, ) が求まると、混合物の三刺激値 X- ， Y_M ， Z_H は、前述した数式 3 0と同様な式に従って算出できる。 When the spectral reflectance R '{?,) Is obtained from Equation 32, the tristimulus values X-, Y _M and Z _H can be calculated according to the same formula as the above-described formula 30.

ところで、各着色剤の吸収係数と散乱係数 S iは誤差を含んでおり、また、数式 3 2〜 3 3は理論式ではなく実験式なので、ステップ S 5 4で得られた三刺激値の計算値 X_M， Y_M ， Z_Mは誤差を有している。この計算誤差は、ステツプ S 5 2で得られた実測値 X。， Y。， Ζ。からの差である。ところで、コンピュータカラーマッチングによる調合割合の予測も、数式 1 , 2， 3 0〜 3 3に従って混合物の三刺激値を求める工程を含んでいる。そこで、濃淡限度見本の調合割合をコンピュータカラ一マッチングで予測する際に、標準色見本における三刺激値の計算誤差を考慮しておけば、その予測精度を向上させることができる。図 2 0のステップ S 5 5では、ステップ S 5 2で得られた標準色の三刺激値の実測俩 X。 , Υ₀ , Ζ₀ と、ステップ S 5 4でられた計算値 Χ_Μ ， Υ_Μ , Ζ_Μ との差分 ΔΧ， ΔΥ, ΔΖを次の数式 3 4に従って求める。 By the way, the absorption coefficient and the scattering coefficient S i of each colorant include an error, and the equations 32 to 33 are empirical rather than theoretical, so the tristimulus values obtained in step S 54 are calculated. The values X _M , Y _M , and Z _M have errors. This calculation error is the measured value X obtained in step S52. , Y. , Ζ. Is the difference from By the way, the prediction of the blending ratio by computer color matching also includes the step of calculating the tristimulus values of the mixture according to the formulas 1, 2, 30 to 33. Therefore, when predicting the mixing ratio of the shading limit sample by computer color matching, if the calculation error of the tristimulus value in the standard color sample is considered, the prediction accuracy can be improved. In step S55 of FIG. 20, an actual measurement 俩 X of the tristimulus value of the standard color obtained in step S52. , Upsilon _0, and Zeta _0, Step S 5 4 out was calculated values Χ _{_Μ,} Υ _Μ, difference ΔΧ with Ζ _Μ, ΔΥ, seek ΔΖ following Equation 3 4.

ΔΧ = Χ₀-ΧΜ . . - (34a) ΔΧ = Χ ₀ -ΧΜ..-(34a)

ΔΥ = Υ₀-Υ_Μ .. - (34b) ΔΥ = Υ ₀ -Υ _Μ ..-(34b)

ΔΖ = Zo - Z_M ... (34c) ΔΖ = Zo-Z _M ... (34c)

ステップ S 5 6では、ステップ S 5 3で設定された濃淡限度見木の三刺激値 X _τ ， Υ.,. , Ζ_τ を上言 tiの差分 ΔΧ， ΔΥ, 厶 Ζで補正することによって、コンピュータカラーマッチングにおける三刺激値の標値 X_c , Y_c , Z_c を求める。すなわち、渙淡限度見本に対する三刺激値の目標値 X_c , Y_c , Z_c は次の数式 3 5で与えられる。 In step S56, the tristimulus values X _τ , Υ.,., Ζ _τ set in step S53 are corrected by the differences ΔΧ, ΔΥ, Ζ of the above ti. target value X _c tristimulus values in computer color matching, Y _c, Ru seek Z _c. That is, the target values X _c , Y _c , and Z _c of the tristimulus values for the lysis limit sample are given by the following Expression 35.

Χτ = Χ_τ-ΔΧ ... (35a) Χτ = _{τ τ} -ΔΧ ... (35a)

Υ_Γ = Υ_Τ-ΔΥ ... (35b) Υ _Γ = Υ _Τ -ΔΥ ... (35b)

ΖΓ = Ζ_Τ-ΔΖ ... (35C) ΖΓ = Ζ _Τ -ΔΖ ... (35C)

ステップ S 5 7では、コンピュータカラーマッチングによって濃淡限度見本の調合割合を予測する。このコンピュータカラーマッチングでは、数式 3 5で与えられる目標値 X_c , γハ， Z_c が得られるような調合割合を求めている。図 2 1は、ステップ S 5 7におけるコンピュータカラーマッチングの詳細丰順を示すフローチャートである。なお、図 2 1の手順は、ニュートン一ラプソン法を用いたコンピュータカラーマッチングを適用したものである。 In step S57, the mixing ratio of the density limit sample is predicted by computer color matching. In this computer color matching, the mixing ratio is determined so that the target values _Xc , γc, and Zc given by Expression ₃₅ are obtained. FIG. 21 is a flowchart showing the detailed procedure of the computer color matching in step S57. The procedure in Fig. 21 applies computer color matching using the Newton-Raphson method.

ステップ S 6 】では、各着色剤（顔料）の調合率 Ci を微小変化させた時の三刺激値の変化を以下の手順で計算する。まず、図 20のステップ S 5 2で用いた標準色 ¾本の調合率 Ci から、〗つの着色剤の調合率のみを AC だけ変えて、その混合物の刺激値の計算値 Χ_Μ'， Υ_Μ', Ζ_Μ 'を求める。これは、前述した数式 1， 2, 3 0〜3 3を用いて、ステップ S 5 2と同様の手順で行なわれる。そして、ステップ S 5 2で得られていた標準色の三刺激値の計算値 Χ„ ， Υ_Μ ， Ζ _Μ との差分 AX_Ci, AY_ci, AZ_Ciを以下の数式 3 6に従って算出する。

In step S6], the change in tristimulus value when the mixing ratio Ci of each colorant (pigment) is slightly changed is calculated in the following procedure. First, from the mixing ratio Ci of the standard color used in step S52 in FIG. 20, only the mixing ratio of one colorant is changed by AC, and the calculated stimulus values of the mixture Χ _Μ ', Υ _Μ ', Ζ _Μ '. This is performed in the same manner as in step S52, using the above-described equations 1, 2, 30 to 33. Then, the differences AX _Ci , AY _ci , and AZ _Ci from the calculated values 三, Υ _，, Ζ _の of the tristimulus values of the standard color obtained in step S ₅₂ are calculated according to the following equation _36. .

ΔΥΟ, = ΥΜ-ΥΜ ...(36b) ΔΥΟ, = ΥΜ-ΥΜ ... (36b)

ΔΖΟ; = Ζ -Ζ_Μ . - . (36c) _{ΔΖΟ;.. = Ζ -Ζ Μ} - (36c)

数式 36を用いると、各着色剤の調合率 Ci のみを微小変化させた時の三刺激値の変化率が、次の数式 3 7で与えられる。

） Using Equation 36, the change rate of the tristimulus value when only the mixing ratio Ci of each colorant is slightly changed is given by the following Equation 37.

)

^ΔΥο, ...₍37b)

濃淡限度色は、標準色の製造誤差を示すものなので、両者の三刺激値は互いに近接した値である。従って、数式 3 7の変化率を用いると、図 20のステップ S 5 6で得られた濃淡限度見本の三刺激値の目標値 X_c ， Y_c , Z_c と、ステツプ S 5 4で得られた標準色見本の三刺激値の計算値 X_M， Y_H， Z_Mとの関係を、次の数式 38で表わすことができる。

^ ΔΥο, ... ₍ 37b)

Since the shading limit color indicates a manufacturing error of the standard color, the tristimulus values of both are close to each other. Thus, with the rate of change of Equation 3 7, obtained in step S 5 the target value X _c tristimulus values of the resulting shade limit sample at 6, Y _c, and Z _c, Sutetsu flop S 5 4 in FIG. 20 The relationship between the calculated tristimulus values X _M , Y _H , and Z _M of the obtained standard color samples can be expressed by the following equation (38).

... (38b)

△C , + AC₂ +AC₃ +ΔΟ) = APVC ... (38d) △ C, + AC ₂ + AC ₃ + ΔΟ) = APVC ... (38d)

なお、数式 38では、着色剤が 4種類であると仮定している。なお、各着色剤の調合率の変化量の合計値 ΔΡ CVには一定値が代入される。例えば、 4 種類の着色剤の調合率 Ci の和∑Ci を一定に保った場合には、？は0である。 Note that Equation 38 assumes that there are four types of colorants. Note that a fixed value is substituted for the total value ΔΡCV of the change amount of the mixing ratio of each colorant. For example, if the sum ∑Ci of the mixing ratios Ci of the four colorants is kept constant,? Is 0.

数式 3 8は、 4つの未知数 ( i = 1〜4) を舍む 4元連立一次方程式なので、これを解くことによって、各未知数の値を求めることができる（ステツプ S 62) 。ステップ S 63では、ステップ S 62で求められた彼を用いて、濃淡限度見本の各着色剤の調合率 C_iTが次の数式 3 9に従って算出される。 Since Equation 38 is a four-dimensional system of linear equations involving four unknowns (i = 1 to 4), the value of each unknown can be obtained by solving this equation (Step S62). In step S63, using him obtained in step _S62 , the mixing ratio C _iT of each colorant of the shading limit sample is calculated according to the following Expression ₃₉ .

C_lT = Ci + ACi .-.(39) C _lT = Ci + ACi .-. (39)

ステップ S 64では、ステップ S 63で得られた調合率 C_iTから、前述した数式 1， 2， 30〜 3 3に従って三刺激値の計算値 , Υ_χ , Z₂ を求める。ステツプ S 65では、こうして得られた三刺激値 , Υ_: , Z_: と、図 20のステツプ S 5 6で求めた三刺激値の目標値 X_c ， Y_c ， Z_c との色差が所定の許容誤差 <5以下であるか否かが判断される。なお、色差は、 L a b表色系を用いて次の数式 40で与えられる。 In step S 64, from the blending ratio C _iT obtained in step S 63, the number mentioned above Equation 1, 2, calculated values of tristimulus values according to. 30 to 3 3, Upsilon _chi, seek Z _2. In scan Tetsupu S 65, thus resulting tristimulus values, _Upsilon:, _Z: the target value X _c tristimulus values obtained in scan Tetsupu S 5 6 of FIG. 20, Y _c, the color difference between Z _c It is determined whether or not the predetermined tolerance is less than or equal to <5. The color difference is given by the following equation 40 using the Lab color system.

ΔΕ = ( - L,)² + (a_c - a !i)² + (b_c - b,!)² ...(40) ΔΕ = (-L,) ² + (a _c -a! I) ² + (b _c -b ,!) ² ... (40)

ステップ S 65において、色差 ΔΕが許容誤差 όよりも小さな場合にはステツプ S 63で得られた調合率 C_iTの値を予測値として採用し、コンピュータカラーマッチングを終了する。一方、色差厶 Eが許容誤差 ό以上の場合には、ステップ S 66において、ステップ S 64で得られた三刺激値の計算値 Χ₁ ， Υ₁ ， Ζ, を、 Χ_Μ , Υ_Μ ， Ζ_Μ に代入して、ステップ S 62〜S 65の処理を繰り返す。こうしてステップ S 62〜S 6 6を繰り返し実行することによって、色差 ΔΕの値が許容誤差 <5よりも小さくなるような調合率 C_iTを求めることができる。 In step S65, if the color difference Δ 許容 is smaller than the allowable error S, The value of the mixing ratio _CiT obtained in step _S63 is adopted as the predicted value, and the computer color matching is terminated. On the other hand, if the color difference E is equal to or larger than the allowable error ό, the calculated values 三₁ , 三₁ , Ζ, of the tristimulus values obtained in step S 64 are replaced with Χ _Μ , Υ _，, Ζ in step S 66. _Μ and the processing of steps S62 to S65 is repeated. By repeatedly executing steps S 62~S 6 6 Thus, the value of the color difference ΔΕ can be obtained compounding ratio C _iT as smaller than the allowable error <5.

図 22は、本発明の第 3実施例において用いた標準色見本の三刺激値と、濃淡限度見本の三刺激値を示すテ一ブルである。サンプル T l， Τ 2, Τ 3は、それぞれ異なる色を有する標準色見本である。図 2 2には、各標準色見本について、図 20のステップ S 5 2で実測された三刺激値 X。， Υ₀ , Ζ,, と、ステップ S 5 4で計算された三刺激値 Χ_Μ， Υ_Μ ， Ζ_Μ とが示されており、また、ステツプ S 5 3で設定された濃淡限度見本の三刺激値の設定値 Χ_τ , Υ_τ ， Ζ_τ も示されている。コンピュータカラーマッチングにおける濃淡限度見本の三刺激値の目標値 X_c , Y_c , .Z_c は、これらの値から、前述した数式 34および 35に従って算出される。 FIG. 22 is a table showing tristimulus values of a standard color sample and tristimulus values of a shading limit sample used in the third embodiment of the present invention. Samples T l, Τ 2, Τ 3 are standard color samples having different colors. Fig. 22 shows the tristimulus values X actually measured in step S52 in Fig. 20 for each standard color sample. , Upsilon _0, and Zeta ,,, Step S 5 4 tristimulus values chi _Micromax calculated in, Upsilon _Micromax, are shown and Zeta _Micromax is also shading boundary samples set in Sutetsu flop S 5 3 The set values of the tristimulus values Χ _τ , Υ _τ , and Ζ _τ are also shown. The target values X _c , Y _c , .Z _c of the tristimulus values of the gray scale limit sample in the computer color matching are calculated from these values according to the above-mentioned equations 34 and 35.

図 23は、第 3実施例における濃淡限度見本の三刺激値の予測結果と、比較例の予測結果と、実調合率とを示したものである。比較例の調合率は、濃淡限度見本の三刺激値の設定値 Χ_τ ， Υ_τ ， Ζ_τ をそのままコンピュータカラーマッチングの目標値とした時に得られる予測値である。また、第 3実施例の調合率は、標準色見本の計算誤差で補正した目標値 X_c , Y_c , Z_c を用いて得られる予測値である。なお、この第 3実施例では、予測精度を検証するために、図 2 3の右端に示す調合率を用いて濃淡限度見本を実際に製作しておき、その三刺激値の実測値を図 22に示す設定値 Χ_τ ， Υ_τ ， Ζ_τ として用いている。 FIG. 23 shows the prediction results of the tristimulus values of the density limit sample in the third embodiment, the prediction results of the comparative example, and the actual blending ratio. The blending ratio of the comparative example is a predicted value obtained when the set values 三_τ , Υ _τ , Ζ _τ of the tristimulus values of the gray scale limit sample are directly used as the target values of the computer color matching. Also, compounding ratio of the third embodiment is the predicted value obtained using the target value X _c were corrected in the calculation error of the standard color swatch, Y _c, the Z _c. In this third example, in order to verify the prediction accuracy, a gray scale limit sample was actually manufactured using the mixing ratio shown on the right end of FIG. 23, and the measured tristimulus values were compared with those in FIG. Are used as the set values Χ _τ , Υ _τ , Ζ _τ shown in Fig.

図 2 3から解るように、第 3実施例は比較例に比べてより高い精度で調合率を予測することが可能であった。また、 3つのサンプル丁】〜 Τ 3では 4つの顔料の成分がかなり異なるが、そのすベての場合において、第 3実施例の予測精度が比較例に比べて高いことが解る。以上のように、上記第 3実施例では、瀘淡限度見本に近い三刺激値を有する標準色見本に関して三刺激摘の計算誤差を求め、これを用いて濃淡限度見本の三刺激値の目標値を補正している。この結果、濃淡限度見本の調合率を予測する際の予測精度を向上させることが可能である。 As can be seen from FIG. 23, the third example was able to predict the mixing ratio with higher accuracy than the comparative example. In addition, although the components of the four pigments are significantly different in the three samples, the prediction accuracy of the third embodiment is higher than that of the comparative example in all cases. As described above, in the third embodiment, the calculation error of tristimulation was calculated for a standard color sample having a tristimulus value close to the filtering limit sample, and this was used to calculate the tristimulus value of the shading limit sample. The target value has been corrected. As a result, it is possible to improve the prediction accuracy when predicting the mixing ratio of the shading limit sample.

と記第 3実施例に関しては、例えば次のような変形も可能である。 With respect to the third embodiment, for example, the following modifications are also possible.

( 1 ) 上記第 3実施例では、濃淡限度見本に近い三刺激値を有するサンプル（近接色サンプル）として標準-色見本を使用し、標準色見本に関する三刺激値の計算誤差を用いて浪淡限度見本の三刺激値の目標値を補正していた。しかし、標準色見本以外の近接色サンプルを選択し、その近接色サンプルに関する計算誤差を用いることも可能である。例えば、種々の色を有するサンプルを収集したデータべ一スを用いて、濃淡限度見本に近い色を有するサンプルを、近接色サンプルとして選択するようにしてもよい。このデータベースは、各着色剤の調合率 C i と、三刺激値の実測値（または、反射率 R ( λ ) ) を少なくとも舍むようにするのが好ましい。なお、データベースを用いて近接色サンプルを検索する場合には、上述の数式 4 0で与えられる色差が、目標とする混合物（濃淡限度見本）に最も近いものを近接色サンプルとすることが好ましい。 (1) In the third embodiment, the standard-color sample is used as a sample having a tristimulus value close to the density limit sample (near-color sample), and the calculation error of the tristimulus value for the standard color sample is used. The target value of the tristimulus value of the light-limit sample was corrected. However, it is also possible to select a nearby color sample other than the standard color sample and use the calculation error for that nearby color sample. For example, using a database in which samples having various colors are collected, a sample having a color close to the shading limit sample may be selected as a nearby color sample. This database should preferably include at least the formula C i of each colorant and the measured tristimulus value (or the reflectance R (λ)). When a close color sample is searched using a database, it is preferable that the closest color sample whose color difference given by Expression 40 is closest to the target mixture (shade limit sample) is used.

なお、このようなデータベースを用いれば、漉淡限度見本の調合率をコンピュータカラーマッチングによって予測する場合に限らず、任意の三刺激値を有する目標混合物の調合率を予測する場合に、その予測精度を向上させることができるという利点がある。 It should be noted that using such a database is not limited to the case where the blending rate of the sample of the strain limit is predicted by computer color matching, but is used to predict the blending rate of a target mixture having arbitrary tristimulus values. There is an advantage that the prediction accuracy can be improved.

( 2 ) 上記第 3実施例では、 Χ Υ Ζ表色系を用いていたが、表色系としては Χ Υ Ζ表色系以外の任意の表色系を使用することができる。例えば、 L * a * b * 表色系などを使用することが可能である。 D . 第 4実施例： (2) Although the third embodiment uses the Χ / Ζ color system, any color system other than the Χ / Ζ color system can be used as the color system. For example, an L * a * b * color system can be used. D. Fourth embodiment:

顔料や染料などの着色剤の着色対象となるもの（被着色物）を所望の色とするには、その所望の色となるように着色剤を種々調合して調合割合を規定する必要がある。この場合、技術者が試行錯誤して所望の色を呈する着色剤の混合割合を求めていたが、近年では、個々の着色剤の濃度別の光学的データを基礎として、着色剤の混合割合を予測するいわゆるコンピュータカラーマッチング（以下、これを適宜 C C Mと略称する）が提案されている（特開平 4一 1 8 1 1 2 9，特公平 6— 9 8 8 8 0 ) 。 In order to obtain a desired color for an object to be colored by a coloring agent such as a pigment or a dye (object to be colored), it is necessary to formulate the coloring agent in various ways so as to obtain the desired color and to define a mixing ratio. . In this case, the technician determines the mixing ratio of the colorant that exhibits the desired color through trial and error. In recent years, so-called computer color matching (hereinafter abbreviated as CCM) has been proposed in recent years, which predicts the mixing ratio of colorants based on optical data for each colorant concentration. (Japanese Patent Application Laid-Open No. Hei 4-181, 1989, Japanese Patent Publication No. Hei 6-98880).

ところで、これら公報で提案された C CMでは、ニュートン 'ラプソン法等の逐次近似法を用いているため、その予測した調合割合が負となることがある。このため、負となった予測調合割合をゼロ若しくは正とするための Π夫がなされている。 By the way, in the CCMs proposed in these publications, since the successive approximation method such as the Newton-Raphson method is used, the predicted mixing ratio may be negative. For this reason, efforts have been made to make the negative forecast compounding ratio zero or positive.

上記した従来の C C Mでは、その被着色物を繊維とし、一旦着色（染色）した繊維の色が所望のものではない場合に、この染色済みの繊維を追加染色してその色を所望のものとする。しかしながら、陶器やタイルでは、その色は、釉に混合する顔料の調合率（調合割合）によって決定され、顔料が調合された釉（以ド、釉薬という）の焼成を経て当該色を呈する。このため、陶器やタイル等のように追加染色（着色）できない被着色物には、上記の従来の C C Mは適用できない。また、追加着色が可能な繊維や追加着色が不可能な陶器，タイルであっても、共に、着色剤（繊維にあっては染料，陶器やタイルにあっては釉薬）は、繊維，陶器等の工業的生産に合わせて多量に繰り返し調合される。この場合、着色剤の調合割合は維持されるが、調合工程に変動、例えば温度や調合のタイミング等のズレがあったりすると、それ以前に調合した着色剤で得られる色や色見本の色を再晃できない場合がある。特に、陶器やタイルでは、天然の顔料を用いる都合上、色の再現の信頼性にやや欠ける。 In the above-mentioned conventional CCM, the object to be colored is a fiber, and if the color of the dyed (dyed) fiber is not the desired one, the dyed fiber is additionally dyed to change the color to the desired one. I do. However, in ceramics and tiles, the color is determined by the mixing ratio (mixing ratio) of the pigment mixed with the glaze, and the color is obtained through firing of the glaze (hereinafter referred to as glaze) in which the pigment is mixed. For this reason, the above-mentioned conventional CCM cannot be applied to an object to be colored that cannot be additionally dyed (colored), such as pottery and tiles. In addition, the colorant (dye for fiber, glaze for pottery and tile) is used for fiber, pottery, etc. Is repeatedly formulated in large quantities according to the industrial production of In this case, the mixing ratio of the colorant is maintained, but if the mixing process fluctuates, for example, if there is a deviation in the temperature, the timing of the mixing, or the like, the color obtained from the colorant previously mixed or the color of the color sample is changed. There is a case that cannot be redone. Especially in ceramics and tiles, the reliability of color reproduction is somewhat lacking due to the use of natural pigments.

このような場合、追加着色が可能な繊維では、上記した C C Mにより追加染色して色修正ができるが、追加染色するために調合した着色剤と既に調合済みの着色剤を併用しなければならず、煩雑である。この煩雑さを回避するためには、色見本の色を得られる着色剤を新たに調合してこの着色剤だけで染色すればよいが、既に調合済みの着色剤は不要として廃棄するか、この調合済みの着色剤に新たな着色剤を追加して再調合する必要がある。また、陶器やタイルでは、追加着色が不可能な都合上、調合済みの着色剤の廃棄か再調合を採ることになる。しかしながら、調合済み着色剤の廃棄は無駄であり、その一方、着色剤の再調合は技術者の長年の勘と経験により着色剤を徐々に追加しながら行なうのでやはり煩雑であつた。 In such a case, for fibers that can be additionally colored, the color can be corrected by additional dyeing using the above-mentioned CCM.However, in order to perform additional dyeing, the colorant that has been prepared must be used in combination with the already prepared colorant. It is complicated. In order to avoid this complication, a colorant that can obtain a color sample color may be newly prepared and dyed with only this colorant. However, the already prepared colorant is discarded as unnecessary or discarded. It is necessary to add a new colorant to the prepared colorant and re-formulate. Also, for pottery and tiles, it is not possible to add additional color, so the prepared colorant must be discarded or remixed. But However, disposing of the prepared colorant is useless, while re-formulation of the colorant is also cumbersome because the colorant is gradually added based on the intuition and experience of engineers.

以下に説明する第 4実施例は、調合済み着色剤の有効利用を図りつつ、再調合を簡略化することを冃的としたものである。 The fourth embodiment described below aims at simplifying re-mixing while effectively utilizing the mixed colorant.

図 2 4は、実施例のコンピュータカラーマッチング方法を実施するための装置を示すブロック図である。この装置は，上述した第 1〜第 3実施例を実現する装置としても使用することが可能である。演算装置 4 0は汎用コンピュータであり，図示しない C P Uがソフトゥヱァプログラムを実行することによって，この発明によるコンピュータカラーマツチングの各工程および各手段を実現する。 FIG. 24 is a block diagram illustrating an apparatus for performing the computer color matching method according to the embodiment. This device can also be used as a device that implements the first to third embodiments described above. The arithmetic unit 40 is a general-purpose computer, and realizes each step and each means of computer color matching according to the present invention by causing a CPU (not shown) to execute a software program.

図 2 5は、この第 4実施例における処理の全体手順を示すフ口一チヤートである。なお、この第 4実施例で対象とする調合物は、陶磁器の素地の表面を覆うための釉である。即ち、顔料を入れないベース釉（基礎釉）が被着色物であり、このベース釉に顔料を添加した釉がコンピュータカラ一マツチングの対象となる調合物である。また、その表色系は、 X Y Z表色系とするが、これ以外の表色系、例えば L * a * b *表色系を採用してもよいことは勿論である。 FIG. 25 is a flowchart showing the entire procedure of the process in the fourth embodiment. The composition targeted in the fourth embodiment is a glaze for covering the surface of a ceramic body. In other words, the base glaze without pigment (base glaze) is the material to be colored, and the glaze with the pigment added to this base glaze is the target for computer color matching. Although the color system is the XYZ color system, it goes without saying that other color systems, for example, the L * a * b * color system may be adopted.

図 2 4に示すように、第 4実施例では、 C CMに関する後述する稀々の演算を実行する演算装置 4 0と、データ入力を行なうキーボード，マウス等の入力機器 4 2と、演算の状態や後述する合否判定の結果を表示する表示機器 4 4と、この合否判定の結果や種々の演算式等を記憶する記憶装置 4 6と、 C CMに必要なデータとしての分光反射率を取得するための分光光度計 4 8とを備える。そして、この演算装置 4 0では、以下の調合処理が行なわれる。 As shown in FIG. 24, in the fourth embodiment, an arithmetic unit 40 for executing a rare operation, which will be described later, relating to the CCM, an input device 42 such as a keyboard and a mouse for inputting data, and a state of the operation And a display device 44 for displaying the results of pass / fail judgment described later, a storage device 46 for storing the results of the pass / fail judgment and various arithmetic expressions, etc., and acquiring a spectral reflectance as data necessary for CCM. And a spectrophotometer 48. Then, in the arithmetic unit 40, the following blending process is performed.

まず、ステップ S 7 1では、目標色を呈する釉（見本釉）を準備する。この場合の H標色は、調合済みの釉で呈される色であるので、この釉（見本釉）における顔料の調合割合は既知である。続くステップ S 7 2では、この冃標色見本釉を分光光度計 4 8で測色し、 X Y Z表色系での色評価値である三刺激値（実測値） Χ_τ ， Υ _τ ， Ζ_τ を求める。この三刺激値は、分光光度計 4 8での測色により得られた目標色見本釉の分光反射率 R ' ( - ) から、次の数式 4 2に従って算出される。なお、算出された三刺激値は、表示機器 4 4に目標見本色とともに表示され、後述の処理に用いるために記憶装置 4 6に記憶される。また、以下に記す三刺激値等の演算結果は、その都度に、記憶装置 4 6に記憶される。

First, in step S71, a glaze (sample glaze) exhibiting a target color is prepared. Since the H standard color in this case is the color presented by the prepared glaze, the blending ratio of the pigment in this glaze (sample glaze) is known. In the following step S72, the 冃 standard color sample glaze is measured with a spectrophotometer 48, and tristimulus values (actually measured values) _{τ τ} , ， _τ , Ζ _τ which are color evaluation values in the XYZ color system Ask for. The tristimulus value is calculated from the spectral reflectance R ′ (-) of the target color sample glaze obtained by colorimetry with a spectrophotometer 48 according to the following equation 42. Is done. The calculated tristimulus value is displayed on the display device 44 together with the target sample color, and is stored in the storage device 46 for use in processing described later. The calculation results of the tristimulus values and the like described below are stored in the storage device 46 each time.

Y_T = K j S( )y( ) R'( )d ... (42b) Y _T = K j S () y () R '() d ... (42b)

Z_T = K |S( )z( ) R'(X)dX ... (42c) Z _T = K | S () z () R '(X) dX ... (42c)

K= 100/ jS(X)y( )dX ... (42d) ここで、 S (；. ) は標準光の分光分布、 X (λ ) , y U) ， z (/) (数式中ではバー付きである。以下同じ）は等色隱であり、い" mも既知の^ ある。なお、式中に記した（ ) は、分光反射率，分光分布，等色関数がいずれも波長 λに依存していることを表わす。 K = 100 / jS (X) y () dX ... (42d) where S (;.) Is the spectral distribution of standard light, X (λ), y U), z (/) (where (The same applies to the bars.) The same is true for color matching, and “m is also known. The parentheses in the formula indicate that the spectral reflectance, spectral distribution, and color matching function are all wavelengths λ. Indicates that it is dependent on

ステップ S 7 3では、上記の見本釉が呈する目標色と近似する色を呈すると想定した調合割合でそれぞれの顔料を添加し、第】回目の試作釉を調合し、この試作釉についてステップ S 7 2と同様にして測色する。これにより、この第 1回目試作釉についての三刺激値（実測値），， Z, を求める。この第】回目試作において、上記の調合割合での顔料の調合は技術者によってなされるものの、その値は任意性のある既知の値であり、従来のように試行錯誤して繰り返し調合する必要はなく、調合割合を定めるに当たって特段の経験や勘も要しない。この場合の第 1回目試作釉についての三刺激値並びに顔料調合量（調合割合）も、表示機器 4 4に表示されると共に、記憶装置 4 6に記憶される。また、この際の顔料調合量は、入力機器 4 2から入力される。 In step S73, each pigment is added at a blending ratio assumed to exhibit a color similar to the target color exhibited by the sample glaze, and the second trial glaze is blended. The color is measured in the same manner as in S72. The tristimulus values (actually measured values),, and Z for the first prototype glaze are obtained. In this second trial, the pigments are blended at the above blending ratios by a technician, but the values are arbitrary values that are known and need to be repeatedly blended by trial and error as in the past. No special experience or intuition is required to determine the mix ratio. In this case, the tristimulus value and the pigment preparation amount (preparation ratio) for the first prototype glaze are also displayed on the display device 44 and stored in the storage device 46. In addition, the amount of the mixture at this time is input from the input device 42.

続くステップ S 7 4では、ステップ S 7 1で作成した目標色見本釉の呈する色とステップ S 7 3における第 1回目試作釉の呈する色との色差 ΔΕ* (JIS Z 8"730)が所定範囲内に納まるか否か、即ち色差 ΔΕ* の合否判定を下す。この場合の色差 ΔΕ* の許容値は、目標色見本の呈する色と第 1回目試作釉の呈する色との違いが一見しては判別できない程度の値、例えば 0. 3〜0. 5程度の値が予め人力機器 4 2から設定される。なお、色差 ΔΕ* の許容値を 0. 3〜 0. 5以外の値とすることもできることは勿論である。 In the following step S74, the color difference ΔΕ * (JIS Z 8 "730) between the color of the target color sample glaze created in step S71 and the color of the first prototype glaze in step S73 is set within a predetermined range. In this case, the allowable value of the color difference ΔΕ * is determined by the difference between the color of the target color sample and the color of the first prototype glaze. Value that cannot be distinguished by The value is set in advance from the manpower device 42. Of course, the allowable value of the color difference ΔΕ * can be set to a value other than 0.3 to 0.5.

このステップ S 7 4で色差 ΔΕ* が 0. 3〜0. 5以内であると合格判定すれば、ステップ S 7 3での第 1回目試作の釉で目標色見本の色を再できるので、それ以上の調合処理は不要であるとして総ての処理を終了する。つまり、第 1回目試作の際の調合割合で調合した新たな釉は、目標色見本釉とほぼ同じ色を呈する。 If it is judged that the color difference ΔΕ * is within 0.3 to 0.5 in step S74, the color of the target color sample can be re-used with the glaze of the first prototype in step S73. It is determined that no further compounding processing is necessary, and all the processing ends. In other words, the new glaze blended at the blending ratio at the time of the first trial production has almost the same color as the target color sample glaze.

一方、ステップ S 7 4で色差厶 Ε* が 0. 3〜0. 5以内に納まらないとして不合格判定した場合には、続くステップ S 7 5で、ステップ S 7 2で求めた目標色見木の三刺激値 Χ_τ ， Υ_τ ， Ζ_τ とステップ S 7 3で求めた第 1回目試作釉の三刺激値 X Υ, , とから、次の数式 4 3に従って三刺激値の差 ΔΧ， ΔΥ， ΔΖを算出する。この三刺激値の差 ΔΧ, ΔΥ, ΔΖは、目標色見木の呈する色と第 1回 Η試作釉の呈する色との色差△ Ε * を反映した値となる。 On the other hand, if the color difference Ε * does not fall within the range of 0.3 to 0.5 in step S74, the rejection is determined, and in step S75, the target color obtained in step S72 is determined. From the tristimulus values 見_τ , Υ _τ , Ζ _τ of the miki and the tristimulus values X Υ,, of the first prototype glaze obtained in step S 73, the difference between the tristimulus values according to the following equation 43 is ΔΧ , ΔΥ, ΔΖ. The differences ΔΧ, ΔΥ, ΔΖ of the tristimulus values reflect the color difference △ * between the color presented by the target swatch and the color exhibited by the first prototype glaze.

ΔΧ = Χ_Τ-Χ) ... (43a) ΔΧ = Χ _Τ -Χ) ... (43a)

ΔΥ =Υ_Τ- Υ, ... (43b) ΔΥ = Υ _Τ -Υ, ... (43b)

ΔΖ = Ζ_Τ-Ζ】，"（43 ΔΖ = Ζ _Τ -Ζ], "(43

ステップ S 7 6では、ステップ S 7 3で行なった第 1回目試作釉の調合時の各顔料（調合割合既知）に微量だけその顔料を追加調合した場合の三刺激値の変化率（撖係数）を、 C CMの手法を用いて、以下の手順に従って求める。なお、各顔料を追加調合する際の微量増加調合量も入力機器 4 2から入力される。 In step S76, the rate of change in tristimulus value (撖 coefficient) when a small amount of each pigment (mixing ratio is known) is added to each pigment (mixing ratio is known) at the time of the first trial glaze blending performed in step S73 Is calculated using the CCM method according to the following procedure. It should be noted that a slight increase in the amount of each additional pigment to be prepared is also input from the input device 42.

まず、このステップ S 7 6の最初のステップ S 9 1では、図 2 6のフ口一チャートに示すように、第 1回目試作釉の呈する色の三刺激値を第 1回目試作時の各顔料の既知の調合割合から C C Mの手法で演算する。この場合の三刺激値 (演算値） X_1/E ， Y_1/E ， Z_1/E は、次の数式 44で表わされる。

First, in the first step S91 of this step S76, as shown in Fig. 26, the tristimulus value of the color exhibited by the first prototype glaze was calculated for each pigment at the time of the first prototype. Calculate using the CCM method from the known blending ratio. The tristimulus values (calculated values) X _{1 / E} , Y _{1 / E} , and Z _{1 / E} in this case are represented by the following Equation 44.

Y_1;E = K ) S( ) y( )R( )d ... (44b) Z_1;F: = K j ζ(λ) R (λ) d λ ... (44c) Y _1; E = K) S () y () R () d ... (44b) Z _{1; F} : = K j ζ (λ) R (λ) d λ ... (44c)

K= 100/ (S( )y( )dX ... (44d) この数式 4 4にあっても、上記した数式 4 2と同様、 S (A ) は標準光の分光分布であり、 X ( ， y (λ) ， ζ (λ) は等色関数（既知）である。し力し、数式 44における分光反射率 R U) は、三刺激値 Χ_1/Ε . Υ_1/Ε , Ζ_1/Ε を顔料の特性から算出する都合上、やはり顔料の特性を用いて以下のようにして算出する。そして、この算出した分光反射率 R ( J を用いて、上記の数式 4 4 から三刺激値 X_1/E ， Y_1/E ， Z_1/E を算出する。 K = 100 / (S () y () dX ... (44d) Even in the equation (44), as in the above equation (42), S (A) is the spectral distribution of the standard light, and X ( , Y (λ) and ζ (λ) are color-matching functions (known), and the spectral reflectance RU in Equation 44 is the tristimulus value Χ _{1 / Ε} . _Ε _{1 / Ε} , Ζ _{1 / For} the sake of convenience in calculating _Ε from the pigment properties, it is also calculated using the properties of the pigment as follows: Then, using the calculated spectral reflectance R (J, the tristimulus value is obtained from the above equation 44. Calculate X1 _{/ E} , Y1 _{/ E} , Z1 _{/ E.}

被着色物と着色剤の吸収係数 {?.) と散乱係数 Si ( .) とは、以下の数式 4 5で表わされるダンカン（Duncan) の式と、数式 4 6で表わされるクベルカ —ムンク（Kubelka-Munk) の混色理論による式の関係にあり、これら数式に基づいて、任意の調合物の分光反射率 R U) を C CMにより求めることができる。

The absorption coefficient {?.) And the scattering coefficient Si (.) Of the coloring matter and the coloring agent are expressed by the following Duncan's equation represented by the following equation 45 and Kubelka-munk represented by the equation 46. -Munk), the spectral reflectance RU) of any formulation can be determined by CCM based on these equations.

SM= ∑ Q S, ...(45b) SM = ∑ Q S, ... (45b)

i=l ここで、 K_M， S_M は調合物（釉）の吸収係数と散乱係数、 Ki ， Si は i番目の成分（顔料）の吸収係数と散乱係数、は i番目の成分（顔料）の調合割合である。 … この数式 4 6を変形すると、分光反射率 R {1) は、次の数式 4 7で与えられる。 i = l where, K _M, S _M is the absorption coefficient and scattering coefficient of the formulation (glazed), Ki, Si is the absorption coefficient and scattering coefficient of the i th component (pigment), the i th component (pigment ). … By transforming this formula 46, the spectral reflectance R {1) is given by the following formula 47.

(47a)(47a)

(47b)

ここで、調合物の吸収係数と散乱係数の比（KZS) は、この数式 47で示されているように、各顔料の吸収係数 Ki (A) と散乱係数 Si (； .) と調合率 Ci で規定される。従って、分光反射率 R (λ ) は、この比（KZS) から算出され、求めることができる。なお、 K_w ， S_w は白色成分（白顔料）の吸収係数と散乱係数、 C_wは白色成分の調合割合である。

Here, the ratio between the absorption coefficient and the scattering coefficient (KZS) of the preparation is represented by the absorption coefficient Ki (A), the scattering coefficient Si (;.) And the mixing ratio Ci of each pigment, as shown in Equation 47. Is defined by Therefore, the spectral reflectance R (λ) can be calculated and calculated from this ratio (KZS). K _w and S _w are the absorption coefficient and scattering coefficient of the white component (white pigment), and C _w is the blending ratio of the white component.

続くステップ S 92では、各顔料（顔料 1, 2, 3) について、第 1回目試作釉の調合割合をそれぞれ微量の調合量だけ侗別に増加させた釉の呈する色についての三刺激値（演算値） X_1/1/E ， Y_1/1/E , Z_1/1/E を、第 1回目試作時の各顔料の既知の調合割合並びに増加させた顔料の既知の調合割合から C CMの手法で演算する。より詳しく説明すると、まず、第〗回目試作釉に顔料 1を微量の調合量 (0. 1 *C, ) だけ加え他の顔料 2 3および白顔料は第 1回目試作釉と同一調合率の釉の呈する色についての三刺激値を演算する。この場合にあっても、上記した数式 44 数式 47が用いられ、数式 47におけるは +0. 1In the following step S92, for each pigment (pigments 1, 2, and 3), the tristimulus value (calculation) for the color of the glaze, in which the proportion of the first prototype glaze was separately increased by a small amount, was calculated. Value) X1 _{/ 1 / E} , Y1 _{/ 1 / E} , Z1 _{/ 1 / E} were calculated from the known blending ratio of each pigment at the time of the first trial production and the known blending ratio of the increased pigment. Calculate using the method described above. To explain in more detail, first, add a small amount of pigment 1 (0.1 * C,) to the first trial glaze, and add other pigments 23 and white pigment to the glaze with the same blending ratio as the first trial glaze. The tristimulus value for the color presented by is calculated. Even in this case, the above Equation 44 Equation 47 is used, and in Equation 47, +0.1

* ) となり、 C₂ ， C₃ ， C_wはそのままである。そして、このように計算された三刺激値 X_1/1/E Y_1/1/E , z_1/1/E は、顔料 1の微量増加調合に起因した三刺激値である。 *), And C ₂ , C ₃ , and C _w remain as they are. Then, thus calculated tristimulus values _{_{X 1/1 / E Y 1}} /1 / E, z 1/1 / E is a tristimulus value due to trace amounts increase preparation of pigment 1.

同様に、顔料 2を微量の調合量（0. 1 *C₂ ) だけ加え他の顔料は同一調合率の釉の呈する色についての三刺激値 X_1/2/E Y_1/2/E Z_1/2/E と、顔料 3を微量の調合量（0. 】 *C₃ ) だけ加え他の顔料は同一調合率の釉の呈する色についての三刺激値 X_1/3/E Y_1/3/E Z_1/3/E とを演算する。なお、白顔料についてもその三刺激値 X_1/w/E Y_1/w/E z_1/w/E が同様に演算される。そして、各顔料 (顔料 1 2 3および白顔料）をそれぞれ微量増加して調合したことによる三刺激愤の変化量と各顔料の微量増加量とから、次の数式 48に従って三刺激俄の変化率（微係数）を算出する。その後は図 25のステップ S 77に進む。 Similarly, Pigment 2 is added in only a small amount (0.1 * C ₂ ), and the other pigments are tristimulus values X _{1/2 / E} Y _{1/2 / E} Z for the color of the glaze of the same mixing ratio. _{1/2 / E} , and pigment 3 with only a small amount (0.) * C ₃ ), and the other pigments have the same stimulus tristimulus value for the color of the glaze X _{1/3 / E} Calculate Y _{1/3 / E} Z _{1/3 / E.} The tristimulus value X _{1 / w / E} Y _{1 / w / E} z _{1 / w / E} is similarly calculated for the white pigment. Then, based on the amount of change in tristimulus による due to each pigment (pigment 123 and white pigment) being slightly increased and mixed and the amount of slight increase in each pigment, the rate of change of tristimulus is calculated according to the following equation 48. (Derivative coefficient) is calculated. Thereafter, the process proceeds to step S77 in FIG.

顔料 1 微量増加調合： d6fcvJ Pigment 1 Minor increase formulation: d6fcvJ

ステップ S 7 7では、目標色見本釉と第 1回試作釉の測色を経てステップ S 7 5で求めた三刺激値の差 ΔΧ, ΔΥ, ΔΖ (実測値の差）を補正するために必要な各顔料の追加増量調合量（顔料 1は , 顔料 2は AC₂ , 顔料 3は AC ₃ ，白顔料は AC_W) を、次の数式 4 9と数式 5 0を用いて算出する。 In step S77, it is necessary to correct the differences Δ 三, ΔΥ, ΔΖ (actual measurement values) of the tristimulus values obtained in step S75 through colorimetry of the target color sample glaze and the first prototype glaze. The required additional compounding amount of each pigment (Pigment 1, AC ₂ for Pigment ₂ , AC ₃ for Pigment ₃ , AC _{W for} white pigment) is calculated using the following Equations 49 and 50.

ΔΧ AC_W ...(49a)ΔΧ AC _W ... (49a)

ΔΥ =—— ACi +—— ACつ +—— AC +—— AC_W ... (49b) ΔΥ = —— ACi + —— AC ++ — AC + —— AC _W ... (49b)

dCi ec₂ ac₃ dc_w dCi ec ₂ ac ₃ dc _w

ΔΖ = AC_W ... (49c)ΔΖ = AC _W ... (49c)

F= AiAC,+ A₂AC₂ + A₃AC₃ + A_VAC_W ... (50) F = AiAC, + A ₂ AC ₂ + A ₃ AC ₃ + A _V AC _W ... (50)

数式 4 9は、三刺激値の差 ΔΧ， ΔΥ, 厶 Ζ (実測値の差）を補正するためのものであり、各顔料の追加増量調合量（ACi , AC₂ ， AC₃ ， AC_W) を変数とする。また、数式 5 0は、各顔料を変動させた場合の費用関数であり、各顔料の追加増量調合量（Δ(^ ， AC₂ ，厶 C₃ ， AC_W) をやはり変数とする。従って、数式 4 9， 5 0を 4つの上記変数について解くことでこれら変数、即ち各顔料の追加増量調合量（厶(^ ， AC₂ ， AC₃ ， AC_W) が定まる。なお、数式 5 0におけるは、顔料 i を単位量変動させるのに必要な费用である。この場合、ステップ S 7 7では、厶が厶 ≥ 0を満たし、 AC_Wが AC_W ≥ 0を満たし、且つ、数式 5 0で表わされる費用関数 Fを最小とする Δ<^ 並びに AC_Wを線形計画法を用いて解く。これにより、並びに AC_Wが負の値として算出されることはない。 Equation 4-9, the difference between the tristimulus values ΔΧ, ΔΥ, provided for correcting厶Ζ (difference measured value), additional bulking formulation amount of each pigment _{_{(ACi, AC 2, AC 3}} , AC W) Is a variable. Equation 50 is a cost function when each pigment is varied, and the amount of addition of each pigment (Δ (^, AC ₂ , m C ₃ , AC _W ) is also a variable. , equation 4 9, 5 0 these variables by solving the four above variables, i.e. additional increase formulation amount of each pigment (厶_{_{(^, AC 2, AC 3}} , AC W) is determined. Here, in equation 5 0 is an expense required to a unit variation amount of the pigment i. in this case, in step S 7 7,厶satisfies厶≥ 0, AC _W satisfies AC _W ≥ 0, and, in equation 5 0 the AC _W a cost function F to delta <^ sequence that minimizes represented solved using linear programming. Accordingly, and never AC _W is calculated as a negative value.

続くステップ S 7 8では、各顔料（白顔料を含む）の調合割合を、求めた追加増量調合量を加味して算出し、追加増量調合した釉についての物性値を求め更新する。より具体的には、この調合割合と各顔料についての吸収係数および散乱係数等を用い、上記の数式 4 5〜数式 47に従って、分光反射率 R (； -) を求める。そして、その後のステップ S 7 9では、求めた分光反射率 R (λ ) を数式 4 4に代入して、各顔抖を追加増量調合した釉の呈する色についての三刺激値（演算値）In the following step S78, the blending ratio of each pigment (including white pigment) is calculated by taking into account the determined amount of added bulk compounding, and the physical property values of the additionally blended glaze are calculated and updated. I do. More specifically, the spectral reflectance R (;-) is obtained from the formulas 45 to 47 using the blending ratio, the absorption coefficient and the scattering coefficient of each pigment, and the like. Then, in the subsequent step S79, the calculated spectral reflectance R (λ) is substituted into Equation 44, and the tristimulus value (calculated value) for the color of the glaze that is obtained by adding and increasing each face bridge is added.

XEND/E ^B D/E' ZBOTJ/E¾ "演算る。 XEND / E ^ B D / E 'ZBOTJ / E¾ "Calculate.

ステップ S 8 0では、この追加増量調合後の三刺激値 X_END/E, Y_END/E, Z_END/E. と目標色見木釉についての三刺激値 Χ_τ , Υ_τ Ζ_τ を用い（JIS Z 8730)、その色差 ΔΕ* が所定範囲内に納まるか否かの合否判定を再度下す。このステツプ S 8 0で色差 ΔΕ* が 0. 3 0. 5以内であると合格判定すれば、それまでのステップで求めた補正量で調合した釉であれば目標色見本の色を再現できることになる。よって、この場合は、調合処理が完了したとして処理を終了する。つまり，求めた各顔料についての補正量（追加増量調合量）を加味した調合割合力調合済みの釉を再調合する際の最終的な調合割合として採用される。より具体的に説明すると、調合工程の変動等により色に変調を来たした調合済み釉は、その調合割合が既知であるので第 4実施例における第 1回目試作釉に相当するので、この色に変調を来たした調合済み釉に上記の補正量に従った配合割合で各顔料を追加して再調合すれば、目標色を呈する釉を再調合できることになる。 In step S80, the tristimulus values X _{END / E} , Y _{END / E} , Z _{END / E.} And the tristimulus values 目標_τ , Υ _τ Ζ _τ for the target swatches are used. (JIS Z 8730), the pass / fail judgment is again made as to whether or not the color difference Δ 納 * falls within a predetermined range. If it is judged that the color difference ΔΕ * is within 0.30.5 in this step S80, if the glaze prepared with the correction amount obtained in the previous steps reproduces the color of the target color sample You can do it. Therefore, in this case, the process is terminated assuming that the blending process is completed. In other words, the blending ratio that takes into account the correction amount (additional blending amount) obtained for each pigment is used as the final blending ratio when the glaze that has been blended is reblended. More specifically, the blended glaze whose color has been modulated due to variations in the blending process, etc., corresponds to the first prototype glaze in the fourth embodiment because the blending ratio is known, If each pigment is added to the blended glaze whose color has been modulated and added at the blending ratio according to the correction amount described above, the glaze exhibiting the target color can be blended again.

一方、ステップ S 8 0で色差 ΔΕ* が 0. 3 0. 5以内に納まらないとして不合格判定した場合には、微係数計算，補正量計算のための三刺激値をそれまでの値からシフトするとともに、各顔料の微量増加量を変更する。具体的に説明すると、数式 4 3における三刺激値 , Υ, , Z_X に、ステップ S 7 9で求めた三刺激値 X_BND/E Y_BND/E, をシフトする。これより、このシフト後の三刺激値（X _/E, Y_END/E, Z _/E) から、ステップ S 7 5では目標色見本との三刺激偾の差（厶 X ΔΥ, ΔΖ) が新たに算出される。 On the other hand, if the color difference ΔΕ * is determined not to be within 0.30.5 in step S80, and the rejection is determined, the tristimulus values for differential coefficient calculation and correction amount calculation are set to the previous values. And change the slight increase of each pigment. More specifically, the tristimulus value X _{BND / E} Y _{BND / E} , which is obtained in step _S79 , is shifted to the tristimulus values, Υ,, and Z _X in Expression 43. From the tristimulus values (X _{/ E} , YEND _{/ E} , Z _{/ E} ) after this shift, the difference (mmXΔΥ, ΔΖ) between the tristimulus and the target color sample is newly obtained in step S75. Is calculated.

また、各顔料（顔料 1 2, 3) について、ステップ S 7 6で考慮した微量增加量を、それまでの調合量（この場合には、顔料 1は Ci + 0. 1 氺）に 0. 1を乗じた微量増加量（（ ^ + 0. 1 ) * 0. 1 ) に変更する。これにより、各顔料は、ステップ S 8 0で不合格とされた場合よりそれぞれ 0. 〗だけ微量増量されるので、ステップ S 7 6で、この微 ftif量をカ喊した調合割合から、各顔料カ^量 ±tftされたことによる≡ί【撤値が上記した数式 4 4〜数式 4 7に従つて求められる。そして、各顔料が微量増加されたことによるこの新たな三刺激値が数式 4 8における三刺激値 X_1/1/E， Υ _1/1/Ε , Z_1/1/E 等に替わって用いられて改めて微係数が算出される。その後は既述したように各顔料の追加増量調合量 (補正量）が求められ、ステップ S 8 0で合格判定されるまで、上記の処理が繰り返される。 In addition, for each pigment (pigments 1 2 and 3), add the microaddition amount considered in step S76 to the previous formulation amount (in this case, pigment 1 is Ci + 0.1 氺). Change to a small increase ((^ + 0.1) * 0.1) multiplied by 1. This allows each pigment to be 0. 0 each less than if it were rejected in step S80. Since the amount is increased by a small amount, in step S76, the amount of each pigment is ± tft, based on the proportion of the fine ftif that has been reduced. It is required according to. Then, this new tristimulus value due to the slight increase of each pigment is used in place of the tristimulus values X 1/1 _{/ E} , Υ _{1/1 / Ε} , Z _1/1 _{/ E,} etc. in Equation ₄₈ . Then, the derivative is calculated again. After that, as described above, the amount of additional compounding (correction amount) of each pigment is determined, and the above-described processing is repeated until a pass is determined in step S80.

以上説明したように第 4実施例のコンピュータカラーマッチング方法では、それまで存在していた釉（目標色見本釉）に対して、その呈する色がある程度近似した色を呈するような調合割合で第 1回目試作釉を調合し、その後、この第 1问目試作釉に、標色見本釉の呈する色と所定範囲で一致するような調合割合で各顔料を追加させた場合の追加増量調合量を求める。このため、顔料の除去を意味する負の調合量を C CMにより求めることがない。よって、第 4実施例のコンピュータカラーマッチング方法によれば、顔料が調合済みの釉（着色剤）の廃棄が不要となるために既存の釉を有効に利用することができる。また、既存の釉の利用に当たり、技術者が関与する工程はステップ S 7 3における一回限りの試作釉の調合であり、この際に、技術者の長年の勘や経験を要しないので、釉の再調合を簡略化することができる。 As described above, in the computer color matching method of the fourth embodiment, the glaze (target color sample glaze) that has existed up to that point has a blending ratio such that the color presents a color that is somewhat similar. After blending the first trial glaze, and then adding each pigment to this first trial glaze at a blending ratio that matches the color of the standard sample glaze within a predetermined range, Ask. For this reason, the amount of negative compounding that means the removal of pigment is not determined by CCM. Therefore, according to the computer color matching method of the fourth embodiment, it is not necessary to dispose of the glaze (colorant) in which the pigment has been prepared, so that the existing glaze can be effectively used. In addition, when using the existing glaze, the process involving the engineer is the one-time trial production of the glaze in step S73. At this time, since the engineer's many years of intuition and experience are not required, the glaze is used. Can be simplified.

更に、第 4実施例のコンピュータカラ一マッチング方法では、釉を再調合するための各顔料の補正量の算出の際に、各顔料の補正量（A C_±並びに A C_W) を、数式 5 0で表わされる費用関数 Fを用いた線形計画法の丰法で求め、各顔料の補正に要する費用が最小となるようにした。このため、第 4実施例のコンピュータカラーマツチング方法によれば、上記した既存の袖の有効利用と再調合に簡略化に加え、コスト低減をも図ることができる。 Further, in the computer color matching method of the fourth embodiment, when calculating the correction amount of each pigment for re-mixing the glaze, the correction amount (AC _{± and} AC _W ) of each pigment is calculated by Equation 50. The cost was calculated by the linear programming method using the expressed cost function F so that the cost required to correct each pigment was minimized. For this reason, according to the computer color matching method of the fourth embodiment, in addition to simplifying the above-mentioned effective use and re-mixing of the existing sleeves, the cost can be reduced.

次に、上記した第 4実施例のコンピュータカラーマッチング方法を行なった際の各ステップでの処理により得られる三刺激値ゃ微係数について、図 2 7を用いて説明する。 Next, the tristimulus value ゃ differential coefficient obtained by the processing in each step when the computer color matching method of the fourth embodiment is performed will be described with reference to FIG.

図 2 7は、ステップ S 7 2で取得した目標色見本釉についての三刺激値（色値）とステップ S 7 3で取得した第 1回目試作釉についての三刺激値の対比、並びに目標色見本釉，第 1回目試作釉における各顔料の調合率（調合率）を表わす。そして、図中の Δ Ε * は、目標色見本釉と第 1回目試作釉との間の色差であり、この値に基づいてステップ S 7 4での合否判定が下される。 Figure 27 shows tristimulus values (color values) for the target color sample glaze obtained in step S72. And the tristimulus values of the first trial glaze obtained in step S73 and the blending ratio (mixing ratio) of each pigment in the target color sample glaze and the first trial glaze. Δ Ε * in the figure is a color difference between the target color sample glaze and the first prototype glaze, and a pass / fail judgment is made in step S74 based on this value.

図 2 8は、図 2 7に掲げる第 1回目試作釉に各顔料をそれぞれ微量だけずつ追加調合した場合の三刺激値の変化率（微係数）を示しており、上記の数式 4 8から演算される。この際、顔料 1 (赤顔料）の微量増量の場合における微係数算出には、演算した三刺激値 X _1/1/E， Y _1/1/E， Z_1/1/Eが用いられ、顔料 2 (黄顔料）では三刺激値 X _1/2/E ， Y _1/2/E ， Z_{1/2 /E} 力、顔料 3 (青顔料）では三刺激値 X _1/3/E ， Y "_3/E ， z_1/3/E 力;、白顔料では三刺激値 X _1/w/E ， Y _1/H/E ， ζ _{1 ί/Ε} が用いられる。 Figure 28 shows the rate of change (differential coefficient) of tristimulus values when each pigment was added to the first prototype glaze shown in Figure 27 in a very small amount. It is calculated from At this time, the calculated tristimulus values X 1/1 _{/ E} , Y 1/1 _{/ E} , and Z 1/1 _{/ E} are used to calculate the differential coefficient when the amount of pigment 1 (red pigment) is slightly increased. In pigment 2 (yellow pigment) tristimulus values _{X 1/2 / E, Y} 1/2 / E, Z 1/2 / E power, the pigment 3 (blue pigment) the tristimulus values X _{1/3 / E,} Y _{3 / E} , z _{1/3 / E} force; For white pigments, tristimulus values X _{1 / w / E} , Y _{1 / H / E} , _{{1} / ί} are used.

この図 2 9は、第 4実施例によるコンピュータカラーマッチング方法による結果を示しており、各顔料について定めた最終的な調合割合で調合した釉と目標色見本釉とを、その調合率と三刺激値について対比して表わす。また、この両釉についての色差 Δ Ε * は 0 . 4 7であり、ステップ S 8 0で合格判定されたことが判る。 FIG. 29 shows the result of the computer color matching method according to the fourth embodiment. The glaze prepared at the final blending ratio determined for each pigment and the target color sample glaze are shown in FIG. The stimulus values are shown in comparison. The color difference ΔΕ * for both glazes was 0.47, indicating that the pass was judged in step S80.

以上の第 4実施例に関しては，次のような変形も可能である。 Regarding the fourth embodiment described above, the following modifications are also possible.

例えば、数式 5 0の費用関数 Fに替わって、以下の数式 5 1で表わされる费用関数 Fを用いてもよい。

For example, instead of the cost function F of Equation 50, a utility function F represented by the following Equation 51 may be used.

更には、ステップ S 7 6にて顔料をそれぞれ微量ずつ追加増量した微係数を計算する際に、各顔料の追加増量調合量の最低値を最低追加増量調合量△ C step- として予め規定しておき、この△ C stepの整数倍ずつ各顔料を追加増量するよう構成することもできる。そして、このように各顔料を追加増量していき、最も目標色見本釉との色差 Δ Ε * が小さくなるときの追加増量量を、微係数算出の際の補正量とし、この補正量を加味した調合割合を調合済みの釉を再調合する際の最終的な調合割合とすればよい。なお、この場合の最低追加増量調合量 A C stepは、各顔料を追加したときにその呈する色が僅かに変わる最小単位の追加増量量として規定される。また、この場合には、各顔料の追加増量量の最大許容量 A Craax を予め定めておき、各顔料の追加増量の総量をこの A Cmax で規定すればよい。そして、この最大許容量 A Cmax は、その顔料を追加したときにその呈する色が大きく変わり、他の顔料の追加堉量では元の色への戻りができないと思われる追加増量として規定すればよレ、。 In addition, when calculating the differential coefficient of each additional pigment by a small amount in step S76, the minimum value of the additional additional compounding amount of each pigment is defined in advance as the minimum additional additional compounding amount △ C step-. In addition, it is also possible to configure so that each pigment is additionally increased by an integral multiple of this △ C step. Each pigment is additionally added in this way, and the additional amount when the color difference Δ Ε * from the target color sample glaze becomes the smallest is used as a correction amount for calculating the differential coefficient, and this correction amount is taken into account. The blending ratio determined may be used as the final blending ratio when re-mixing the blended glaze. In this case, the minimum additional compounding amount AC The step is defined as the minimum increment of the minimum unit in which the color changes slightly when each pigment is added. In this case, the maximum allowable amount A Craax of the additional amount of each pigment may be determined in advance, and the total amount of the additional amount of each pigment may be defined by this A Cmax. The maximum allowable amount A Cmax is defined as an additional increase in the color that the color changes when the pigment is added, and it is considered that the original color cannot be returned with the additional amount of the other pigment. Yeah.

また、上記の费用関数 Fを用いた線形計画法に限られるわけではなく、他の手法を採ることもできる。例えば、上記の第 4実施例で示したように、色差 Δ Ε * で合否判定する場合には、その合否判定にある程度の幅が許容される。従って、この合否判定にある程度の瞹眛さを導入したいわゆるファジィ線形計画法の手法を採ることもできる。このフアジィ«計西肤の手法を採ったコンピュータカラーマッチング方法での結果を以下の図 3 0に示す。この図 3 0に示すように、各顔料について定めた最終的な調合割合で調合した釉と目標色見本釉とは、その色差 Δ Ε * が 0. 2 0であり、両者の色はよく一致していることが判る。 Further, the present invention is not limited to the linear programming using the above-mentioned function F, and other methods can be adopted. For example, as shown in the fourth embodiment, when the pass / fail judgment is made based on the color difference ΔΕ *, a certain width is allowed for the pass / fail judgment. Therefore, a so-called fuzzy linear programming method that introduces a certain degree of ambiguity in this pass / fail decision can be adopted. The result of the computer color matching method using the method of Fujii «Keisei 肤 is shown in Figure 30 below. As shown in Fig. 30, the glaze blended at the final blending ratio determined for each pigment and the target color sample glaze have a color difference ΔΕ * of 0.20, and both colors are well It turns out that they match.

また、上記の第 4実施例では、陶器やタイルを色付けする釉を例に採り説明したが、繊維を染める染色剤にも適用できることは勿論である。 Further, in the fourth embodiment, the glaze for coloring pottery and tiles has been described as an example, but it is needless to say that the present invention can be applied to a dye for dyeing fibers.

[産業上の利用可能性] [Industrial applicability]

この発明にかかるコンピュータカラーマツチング方法および装置は、陶器ゃタィルを色付けするための釉に混合される着色剤の調合割合の予測とその混合物 A computer color matching method and apparatus according to the present invention include a method for predicting a mixing ratio of a colorant to be mixed with a glaze for coloring pottery tiles and a mixture thereof.

(すなわち釉）の色予測に適用できる他に，繊維を染めるための染色剤等の各種の混合物における着色剤の調合割合の予測とその混合物の色予測等に適用可能である。 In addition to being applicable to the color prediction of (glaze), it can be applied to the prediction of the blending ratio of the colorant in various mixtures such as dyes for dyeing fibers and the color prediction of the mixture.

Claims

請求の範囲 The scope of the claims

1 . コンピュータカラ一マッチングによって着色剤の調合割合の予測または混合物の色予測を行なう方法であって、 1. A method of estimating the mixing ratio of a colorant or the color of a mixture by computer color matching,

無色透明でない被着色物に白色着色剤を混合した第 1の混合物の散乱係数 s _w. ，を、前記白色着色剤の調合率 c_w に依存した関数 f ( C_w ) として準備するェ程と、 Preparing a scattering coefficient s _w ., Of a first mixture obtained by mixing a white colorant with a colorless and transparent material as a function f (C _w ) depending on the mixing ratio c _w of the white colorant; ,

前記散乱係数 s _w，を基準として、前記第 1の混合物の吸収係数 K _w，を前記調合率 c„ に依存する形式で求める工程と、 Obtaining an absorption coefficient K _w , of the first mixture in a form depending on the mixing rate c „, based on the scattering coefficient s _w ,

前記散乱係数 s _w' を基準として、白色でない有色着色剤の吸収係数 κ ,, と散乱係数 s _p とを前記有色着色剤の調合率 c_p に依存する形式で求める工程と、所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数 K_w' ， K _p および前記散乱係数 s_w' ， S _P を用いたコンピュータカラーマッチングを行なうことによって求める工程と、 The basis of the scattering coefficient s _w ', a step of determining a format which depends the absorption coefficient kappa ,, and the scattering coefficient s _p colored colorant not white in formulation ratio c _p of the colored colorant, desired color formulation proportion of the coloring agent in order to adjust the mixture with use or the color of the mixture produced in a predetermined compounding ratio, the absorption coefficient K _w ', K _p and the scattering coefficient s _w', the S _P, The process required by performing computer color matching,

を備えるコンピュータカラ一マッチング方法。 A computer color matching method comprising:

2. 請求項 1記載のコンピュータカラーマッチング方法であって、 2. The computer color matching method according to claim 1, wherein

( a ) 前記被着色物に前記白色着色剤を混合して、前記白色着色剤の調合率 C _w が異なる複数個の第 1の混合物を作成するとともに、前記複数個の第 1の混合物の分光反射率をそれぞれ測定する工程と、 (A) said mixing the white colorant to be colored object, the spectral of the with Formulation ratio C _w of the white colorant to create a plurality of different first mixture, the first mixture of said plurality Measuring the reflectivity,

( c ) 前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記有色着色剤の調合率 C_p が異なる複数個の第 3の混合物を作成するとともに、前記複数個の第 3の混合物の分光反射率をそれぞれ測定する工程と、 (d) 前記被着色物に前記白色着色剤と前記有色着色剤とを混合して、前記第 3 の混合物とは調合割合が異なる第 4の混合物を作成するとともに、前記第 4の混合物の分光反射率を測定する工程と、 (C) said mixing with the white colorant to be colored object and the colored colorant, along with Formulation ratio C _p of the colored colorant to create a plurality of different third mixture, several said double Measuring the spectral reflectances of the third mixture of (d) mixing the white colorant and the colored colorant with the object to be colored to form a fourth mixture having a different mixing ratio from the third mixture, and the fourth mixture. Measuring the spectral reflectance of the

(e) 前記複数個の第 1の混合物の散乱係数 S_w' を前記白色着色剤の調合率 C_w の関数 f (C_w ) によって表わすとともに、前記関数 f (C_w ) に含まれる係数の値を仮決定する工程と、 (e) The scattering coefficient S _w ′ of the plurality of first mixtures is represented by a function f (C _w ) of the mixing ratio C _w of the white colorant, and a coefficient f included in the function f (C _w ) Tentatively determining a value;

(f ) 前記複数個の第 1の混合物の分光反射率の測定値と前記関数 f (C_w ) とを用いて、前記調合率 c_w に依存した形式で前記第 1の混合物の吸収係数 K_w' を求める工程と、 (f) Using the measured values of the spectral reflectances of the plurality of first mixtures and the function f (C _w ), the absorption coefficient K of the first mixture in a form dependent on the blending rate c _w _w '

(g) 前記第 2の混合物の分光反射率と、前記複数個の第 3の混合物の分光反射率と、前記関数 f (C_w ) と、前記吸収係数 K_w' とを用いて、前記調合率 C_p に依存した形式で前記有色着色剤の吸収係数 K_p と散乱係数 S _Ρ とをそれぞれ求める工程と、 (g) using the spectral reflectance of the second mixture, the spectral reflectances of the plurality of third mixtures, the function f (C _w ), and the absorption coefficient K _w ′, Determining the absorption coefficient K _p and the scattering coefficient S _{の of the} colored colorant in a form dependent on the rate C _p ;

(h) 前記吸収係数 K_w，， Κ_Ρ および前記散乱係数 S_w，， S_p を用いて、前記第 4の混合物に関するコンピュータカラーマッチングを行なうとともに、前記コンピュータカラーマッチングによる予測精度を向上させるように前記関数 f (C _B ) に舍まれる係数を修正する工程と、 (h) the absorption coefficient K _w,, kappa _[rho and the scattering coefficient S _w,, using S _p, performs a computer color matching for said fourth mixture, to improve the prediction accuracy by the computer color matching Correcting the coefficient included in the function f (C _B ),

( i ) 前記工程（f ) ないし（h) を繰り返すことによって前記関数 f (C_w ) に舍まれる係数を決定する工程と、 (i) determining a coefficient included in the function f (C _w ) by repeating the steps (f) to (h);

を備えるコンピュータカラーマッチング方法。 Computer color matching method comprising:

3. 前記関数 f (C_H ) が定数 _ai と係数 S_B とを舍む次の式で与えられる、請求項 1または 2記載のコンピュータカラーマッチング方法： 3. The computer color matching method of claim 1 or 2, wherein the function f (C _H ) is given by the following equation that includes a constant _ai and a coefficient S _B :

f(C_w ) = (C_H +a, -S_B )/(C_w +a, ). f (C _w ) = (C _H + a, -S _B ) / (C _w + a,).

4. 請求項 1記載のコンピュータカラ一マッチング方法であって、前記コンピュータカラーマッチングを行なう T-程は， (a) 複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備する工程と、 4. The computer color matching method according to claim 1, wherein the T-step of performing the computer color matching comprises: (a) mixing a plurality of coloring agents to prepare a plurality of samples having different mixing ratios;

(b) 前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、 (b) measuring the spectral reflectances of the plurality of samples, respectively, and obtaining actual measured values of coordinate values of a predetermined color system representing respective colors of the plurality of samples from the measured values of the spectral reflectances; When,

(c) 前記複数のサンプルのそれぞれに関して、前記吸収係数 K_w' ， K_P および前記散乱係数 S_w' ， S_p を用いて前記表色系の座標値の予測値を求めるとともに，前記表色系の座標値の予測誤差を算出する工程と、 (c) for each of said plurality of samples, the absorption coefficient K _w ', K _P and the scattering coefficient S _w', also when using the S _p obtain the prediction value of the coordinate values of the color system Calculating a prediction error of the coordinate values of the color system;

( d ) 前記複数のサンプルに関する前記表色系の座標储と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、 (d) analyzing the relationship between the coordinate 储 of the color system and the prediction error for the plurality of samples by a predetermined error correction method,

(e) 前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測を，前記吸収係数 K_w' ， _P および前記散乱係数 S_w' ， S_p を用いて求める工程と、 (e) While correcting the target value or the predicted value of the computer color matching using the error correction method, the prediction of the mixing ratio of the colorant of the new mixture or the color prediction of the mixture is performed by the absorption coefficient K _w ′. , _P and the scattering coefficients S _w ′, S _p ,

5. 請求項 1記載のコンピュータカラーマッチング方法であって、前記コンピュータカラーマツチングを行なう工程は， 5. The computer color matching method according to claim 1, wherein the step of performing the computer color matching comprises:

(b) 前記近接色サンプルの既知の調合割合から、前記吸収係数 K_w' ， K_p および前記散乱係数 S_w' ， S_p を用いて前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める工程と、 (b) from a known compounding ratio of the proximity color sample, the absorption coefficient K _w ', K _p Contact and the scattering coefficient S _w', of the color system representing the color of the adjacent color sample using the S _p Obtaining a calculated value of the coordinates, obtaining a calculation error from the measured value and the calculated value;

(c) 目標混合物の色に対する前記表色系の座標の目標値を設定する工程と、 (c) setting a target value of the coordinates of the color system for the color of the target mixture;

(d) 前記計算誤差を用いて前記目標値を補正し、補正後の目標値と前記吸収係数 Κ„' ， Κ_Ρおよび前記散乱係数 S„' ， S_pとを用いたコンピュータカラーマツチングによって、前記目標混合物の着色剤の調合割合を予測する工程と、を備えるコンピュータカラーマッチング方法。 ; (d) using the calculated error by correcting the target value, the corrected target value and the absorption coefficient number kappa of "', kappa _[rho and the scattering coefficient S"', computer color pine using the S _p Estimating the blending ratio of the colorant in the target mixture by ching.

6 . 請求項 1記載のコンピュータカラーマッチング方法であって，前記コンピュータカラーマツチングを行なう工程は， 6. The computer color matching method according to claim 1, wherein the step of performing the computer color matching comprises:

複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求める工程を舍み， Providing a process of determining a blending ratio of the plurality of colorants so that the blended product of the plurality of colorants exhibits a color similar to a desired target color;

前記調合割合を求める工程は， The step of obtaining the blending ratio comprises:

( a ) 前記目標色を呈する調合物見本について、所定の表色系での色評価値の実測値を求める丁-程と、 (a) obtaining a measured value of a color evaluation value in a predetermined color system for a preparation sample exhibiting the target color;

( b ) 既知の調合割合で前記着色剤が調合された 1次調合物について、前記所定の表色系での色評価値の実測値を求める工程と、 (b) a step of obtaining an actual measurement value of a color evaluation value in the predetermined color system, for a primary preparation in which the colorant is prepared at a known preparation ratio,

( c ) 前記〗次調合物についての前記既知の調合割合に基づいて、前記 1次調合物の呈する色の前記所定の表色系での色評価値の計算値を，前記吸収係数 K_w' ， K _p および前記散乱係数 S _w' ， S _p を用いて求める工程と、（d ) 前記 1次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を，前記吸収係数 K_w' ， K _p および前記散乱係数 S _w' ， S _p を用いて求めるとともに、前記 1次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、 (c) calculating the calculated color evaluation value of the color of the primary formulation in the predetermined color system based on the known formulation ratio of the〗 -order formulation, using the absorption coefficient K _w ′ , K _p and the scattering coefficients S _w ′, S _p , and (d) the above-mentioned specifications for the colorant-extended formulation assuming that the colorant was increased in the primary formulation. of the calculated values of the color evaluation value in a color system, the absorption coefficient K _w ', K _p and the scattering coefficient S _w', with obtained using S _p, the colorant bulking formulation from the primary formulation Determining the amount of change in the calculated value of the color evaluation value for the object;

( e ) 前記調合物見木と前記 1次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれそれの増量補正量を算出する工程と、 (e) The coloring based on the change in the calculated value of the color evaluation value so that the difference between the measured value of the color evaluation value of the mixture look-ahead and the primary formulation is equal within a predetermined range. Calculating the respective amount-increasing correction amounts of the agent;

7 . コンピュータカラーマッチングによつて着色剤の調合割合の予測または混合物の色予測を行なう方法であって、 7. A method for predicting the mixing ratio of a colorant or the color of a mixture by computer color matching,

( a ) 複数の着色剤を混合して、調合率が互いに異なる複数のサンプルを準備する工程と、 (a) Prepare multiple samples with different blending ratios by mixing multiple colorants Process,

( b ) 前記複数のサンプルの分光反射率をそれぞれ測定するとともに、前記分光反射率の測定値から、前記複数のサンプルのそれぞれの色を表わす所定の表色系の座標値の実測値を求める工程と、 (b) a step of measuring the spectral reflectances of the plurality of samples, respectively, and obtaining an actual measurement value of a coordinate value of a predetermined color system representing each color of the plurality of samples from the measured values of the spectral reflectances; When,

( c ) 前記複数のサンプルのそれぞれに関して、前記表色系の座標値の予測誤差を算出するて程と、 (c) calculating a prediction error of coordinate values of the color system for each of the plurality of samples;

( d ) 前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する工程と、 (d) analyzing the relationship between the coordinate value of the color system and the prediction error for the plurality of samples by a predetermined error correction method,

( e ) 前記誤差補正法を用いてコンピュータカラーマッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう T-程と、 (e) While correcting the target value or the predicted value of the computer color matching using the error correction method, the prediction of the mixing ratio of the colorant in the new mixture or the color prediction of the mixture is performed by the computer color matching. About

8 . 請求項 7記載のコンピュータカラーマッチング方法であって、 8. The computer color matching method according to claim 7, wherein

前記工程（d ) は、前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係をニューラルネッ卜ワークに学習させる工程、を舍み、 The step (d) comprises a step of causing a neural network to learn a relationship between the coordinate values of the color system and the prediction error for the plurality of samples.

前記工程（e ) は、学習済みのニューラルネットワークを用いてコンピュータカラ一マツチングによる予測を行なう工程を舍む、 The step (e) includes a step of performing prediction by computer matching using a trained neural network.

コンピュータカラーマッチング方法。 Computer color matching method.

9 . 請求項 8記載のコンピュータカラーマツチング方法であって、 9. The computer color matching method according to claim 8, wherein

前記ニューラルネットワークは、 3つのニューロンで構成される入力層と、複数のニューロンを舍む中間層と、 3つのニューロンで構成される出力層と、で構成される三層の階層構造を有する、コンピュータカラーマッチング方法。 The neural network has a three-layer hierarchical structure composed of an input layer composed of three neurons, an intermediate layer composed of a plurality of neurons, and an output layer composed of three neurons. , Computer color matching method.

1 0 . 所望の色を有する目標混合物の着色剤の調合割合をコンピュータカラ一マツチングによって予測する方法であって、 ( ) 着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める工程と、 10. A method for predicting the blending ratio of a colorant in a target mixture having a desired color by computer color matching, (D) obtaining the actual measurement value of the coordinates of a predetermined color system for a nearby color sample having a colorant blending ratio known and having a color close to the desired color;

( b ) 前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める工程と、 (b) obtaining a calculated value of the coordinates of the color system representing the color of the adjacent color sample from a known blending ratio of the adjacent color sample, and obtaining a calculation error from the measured value and the calculated value;

( c ) ifS目合物の色に文る IWti表色系の座標の目標値を設定する工程と、 ( d ) 前記計算誤差を用いて前記目標値を補正し、補正後の目標値を用いてコンピュータカラーマツチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する工程と、 (c) setting a target value of the coordinates of the IWti color system written in the color of the ifS compound; and (d) correcting the target value using the calculation error, and using the corrected target value. Estimating the mixing ratio of the colorant of the target mixture by performing computer color matching with

1 1 . 請求項 1 0記載のコンピュータカラ一マッチング方法であって、 11. The computer color matching method according to claim 10, wherein

前記工程（a ) は、複数のサンプルに関して、着色剤の調合割合と、前記表色系の座標の実測値とを舍むデータベースから、前記目標混合物との色差が最小となるサンプル ¾1択することによって前記近接色サンプルを検索する工程を含む、コンピュータカラーマッチング方法。 In the step (a), for a plurality of samples, a sample that minimizes the color difference from the target mixture is selected from a database that stores the mixing ratio of the colorant and the actually measured values of the coordinates of the color system. Computer color matching method, comprising: searching for the close color sample according to:

1 2 . 複数の着色剤を調合した調合物が所望の目標色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラ一マッチング方法であって、 12. A computer color matching method for determining a blending ratio of a plurality of colorants so that a blended product of a plurality of colorants exhibits a color similar to a desired target color,

( c ) 前記 1次調合物についての前記既知の調合割合に基づいて、前記 1次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める工程と、（ d ) 前記 1次調合物に前記着色剤を増量補正したと仮定した着色剤増量調合物についての前記所定の表色系での色評価値の計算値を求め、前記 1次調合物から前記着色剤増量調合物への前記色評価値の計算値の変化量を求める工程と、 (c) a step of obtaining a calculated value of a color evaluation value of the color of the primary formulation in the predetermined color system based on the known blending ratio of the primary formulation; (d) Calculate the color evaluation value in the predetermined color system for the colorant-extended formulation assuming that the colorant has been increased in the primary formulation and calculate the color evaluation value from the primary formulation. Determining the amount of change in the calculated value of the color evaluation value to the colorant increase formulation,

( e ) 前記調合物見木と前記 1次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの增量補正量を算出する工程と、を備える (e) The coloring based on the change in the calculated value of the color evaluation value so that the difference between the measured value of the color evaluation value of the mixture look-ahead and the primary formulation is equal within a predetermined range. Calculating a correction amount for each of the agents.

1 3 . 請求項 1 2記載のコンピュータカラ一マッチング方法であって、 13. The computer color matching method according to claim 12, wherein

前記工程（d ) は、前記 1次調合物における前記着色剤の調合量に比べて微量の量の前記着色剤を前記 1次調合物に増量補正したと仮定した場合について、前記色評価値の計算値の変化量を求める工程を舍む。 In the step (d), the color evaluation value is based on a case where it is assumed that a small amount of the coloring agent is corrected to be added to the primary formulation compared to the blending amount of the coloring agent in the primary formulation. A process for obtaining the amount of change in the calculated value of is established.

1 4 . 請求項 1 2又は請求項 1 3記載のコンピュータカラーマッチング方法であつて、 14. The computer color matching method according to claim 12 or claim 13, wherein:

前記工程（e ) は、前記着色剤の増量に伴う派生費用を表わす費用関数を用いた線形計画法にて、前記着色剤のそれぞれについての最小の増量補正量を算出する工程を舍む。 The step (e) includes a step of calculating the minimum amount of increase correction for each of the colorants by a linear programming method using a cost function representing a derivation cost accompanying the increase in the amount of the colorant.

1 5 . コンピュータカラーマッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置であって、 15. An apparatus for predicting the mixing ratio of colorants or predicting the color of a mixture by computer color matching,

無色透明でない被着色物に白色着色剤を混合した第】の混合物の散乱係数 S _w_ ' を、前記白色着色剤の調合率 C_w に依存した関数 f ( C_w ) として作成する手段と、 Means for preparing the scattering coefficient S _w _ 'of the second mixture obtained by mixing a white colorant with a colorless and non-transparent substance as a function f (C _w ) depending on the blending ratio C _w of the white colorant; ,

前記散乱係数 S を基準として、前記第 1の混合物の吸収係数 K _w，を前記調合率 c_w に依存する形式で求める手段と、 Means for determining an absorption coefficient K _w , of the first mixture in a form depending on the mixing rate c _w , based on the scattering coefficient S,

前記散乱係数 S を基として、白色でない有色着色剤の吸収係数 K _P と散乱係数 S _p とを前記有色着色剤の調合率 C_p に依存する形式で求める手段と、所望の色を有する混合物を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合物の色を、前記吸収係数 K_w' ， K _p および前記散乱係数 S _w' ， S _P を用いたコンピュータカラーマッチングを行なうことによって求める手段と、 Means for determining, on the basis of the scattering coefficient S, an absorption coefficient K _P and a scattering coefficient S _p of a non-white colored colorant in a form depending on the mixing rate C _p of the colored colorant; and a mixture having a desired color. The proportion of the colorant to adjust the color, or The color of the mixture produced in formulation ratio, the absorption coefficient K _w ', K _p and the scattering coefficient S _w', and means for obtaining by performing a computer color matching using S _P,

を備えるコンピュータカラ一マッチング装置。 A computer color matching device comprising:

1 6 . コンピュータカラ一マッチングによって着色剤の調合割合の予測または混合物の色予測を行なう装置であって、 16. An apparatus for predicting the mixing ratio of colorants or predicting the color of a mixture by computer color matching,

複数の着色剤を混合することによって作成された調合率が互いに異なる複数のサンプルの分光反射率をそれぞれ測定する手段と， Means for measuring the spectral reflectances of a plurality of samples having different mixing ratios formed by mixing a plurality of colorants, respectively;

前記複数のサンプルに関する前記表色系の座標値と前記予測誤差との関係を、所定の誤差補正法で分析する丰段と、 Analyzing the relationship between the coordinate values of the color system for the plurality of samples and the prediction error by a predetermined error correction method,

前記誤差補正法を用いてコンピュータカラ一マッチングの目標値または予測値を補正しつつ、新たな混合物の着色剤の調合割合の予測または混合物の色予測をコンピュータカラーマッチングにより行なう手段と、 Means for correcting the target value or the predicted value of the computer color matching using the error correction method, and performing the prediction of the mixing ratio of the colorant of the new mixture or the color prediction of the mixture by computer color matching,

1 7 . 所望の色を有する目標混合物の着色剤の調合割合をコンピュータカラーマツチングによって予測する装置であって、 17. An apparatus for predicting, by computer color matching, the blending ratio of a colorant of a target mixture having a desired color,

着色剤の調合率が既知で前記所望の色に近い色を有する近接色サンプルについて所定の表色系の座標の実測値を求める手段と、 Means for determining an actual measurement value of coordinates of a predetermined color system for a proximity color sample having a colorant mixing ratio known and having a color close to the desired color;

前記近接色サンプルの既知の調合割合から、前記近接色サンプルの色を表わす前記表色系の座標の計算値を求め、前記実測値と前記計算値から計算誤差を求める手段と、前記目標混合物の色に対する前記表色系の座標の目標値を設定する手段と、前記計算誤差を用いて前記標値を補正し、補正後の目標値を用いてコンピュータカラーマッチングを実行することによって、前記目標混合物の着色剤の調合割合を予測する手段と、 Means for calculating a calculated value of the coordinates of the color system representing the color of the close color sample from the known blending ratio of the close color sample, and calculating a calculation error from the measured value and the calculated value; Means for setting a target value of the coordinates of the color system for the color of the target mixture; correcting the target value using the calculation error; and performing computer color matching using the corrected target value. Means for predicting the blending ratio of the colorant in the target mixture,

を備えるコンピュータカラーマッチング装置。 A computer color matching device comprising:

1 8 . 複数の着色剤を調合した調合物が所望の目檫色に近似した色を呈するように該複数の着色剤についての調合割合を求めるコンピュータカラーマッチング装置であって、 18. A computer color matching apparatus for determining a blending ratio of a plurality of colorants such that a blended product of a plurality of colorants exhibits a color similar to a desired target color,

既知の調合割合で前記着色剤が調合された〗次調合物について、前記所定の表色系での色評価値の実測値を求める手段と、 Means for obtaining an actual measurement value of a color evaluation value in the predetermined color system, for a secondary formulation in which the colorant is formulated at a known formulation ratio,

前記〗次調合物についての前記既知の調合割合に基づいて、前記 1次調合物の呈する色の前記所定の表色系での色評価値の計算値を求める手段と、 Means for obtaining a calculated value of a color evaluation value of the color represented by the primary formulation in the predetermined color system, based on the known formulation ratio of the primary formulation;

前記調合物見本と前記 1次調合物との前記色評価値の実測値の差が所定範囲で一致するように、前記色評価値の計算値の変化量に基づいて前記着色剤のそれぞれの增量補正量を算出する手段と、を備える Each of the coloring agents based on the amount of change in the calculated value of the color evaluation value so that the difference between the measured value of the color evaluation value of the preparation sample and that of the primary preparation matches within a predetermined range. Means for calculating the mass correction amount of

コンピュータカラ一マッチング装置。 Computer color matching device.

1 . コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法であって、 1. A method for determining the absorption and scattering coefficients of colorants used in computer color matching,

無色透明でない被着色物に白色着色剤を混合した第 1の混合物の散乱係数 S _w_ ' と吸収係数 K_w' を決定する工程と、前記散乱係数 S_w，と吸収係数 K_w，の少なくとも一方を基準として、 A色でない有色着色剤の吸収係数 Κ_μ と散乱係数 S の相対値をそれぞれ決定する I：程と、 A step of determining a scattering coefficient S _w _ 'and an absorption coefficient K _w ' of a first mixture obtained by mixing a white colorant with a colorless and transparent object; Based on at least one of the scattering coefficient S _w , and the absorption coefficient K _w , the relative values of the absorption coefficient Κ _μ and the scattering coefficient S of the colorant other than A are determined.

を備える方法。 A method comprising:

20. 前記第〗の混合物の散乱係数 S_w' 力，前記白色着色剤の調合率 ^.と定数 a, と係数 S_B とを舍む次の関係で与えられる、請求項 1 9記載の方法： 20. scattering coefficient S _w 'force of the mixture of the first〗, the preparation of the white colorant ^. And constants a, and given a coefficient S _B in Complex free following relationship The method of claim 1 9, wherein :

S_w， =(C_W +a, -S_B )/(C_w +a】 ). S _w , = (C _W + a, -S _B ) / (C _w + a]).

2 1. 前記第 1の混合物の散乱係数 S_w' 力前記白色着色剤の調合率 C_wと定数 a, b， d, e ,C_W。を含む次の関係で与えられる、請求項 1 9記載の方法： s_w， = a'C_w + b (C_wo ≤ C_w); 2 1. Scattering coefficient S _w 'force of the first mixture The mixing ratio C _w of the white colorant and constants a, b, d, e, and C _W. 10. The method of claim 19, given by the following relation: s _w , = a'C _w + b (C _wo ≤ C _w );

S_w， = d-C_w + e (C_w < C_wo). S _w , = dC _w + e (C _w <C _wo ).

22. コンピュータカラーマッチングを利用して作成された混合釉であって、無色透明でないベース釉に白色着色剤を混合した第 1の混合釉の散乱係数 S_w_22. Scattering coefficient S _w _ of the first mixed glaze, which is a mixed glaze made using computer color matching and is a colorless and transparent base glaze mixed with a white colorant

' と吸収係数 K_w' を決定し、 'And the absorption coefficient K _w '

前記散乱係数 S_w' と吸収係数 K_w' の少なくとも一方を基準として、白色でない有色着色剤の吸収係数 K_P と散乱係数 S_P の相対値をそれぞれ決定し、前記吸収係数 K_w，， K_p および前記散乱係数 S_w，， S_p を用いたコンビュ一タカラーマッチングを行なうことによって，所望の色を有する混合釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色を予測し、 Based on at least one of the scattering coefficient S _w ′ and the absorption coefficient K _w ′, a relative value of the absorption coefficient K _P and the scattering coefficient S _P of the non-white colorant is determined, and the absorption coefficient K _w , , K _p and the scattering coefficient S _w, by performing Konbyu one Takara over matching using S _p, formulation ratio of the colorant to adjust the mixed glaze having a desired color or a predetermined compounding ratio Predict the color of the mixed glaze generated by

前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する， Mixing the base glaze, the white colorant, and the colored colorant to have the predicted blending ratio or the predicted color;

ことによって生成された混合釉。 Mixed glaze produced by.

2 3 . コンピュータカラーマッチングを利用して作成された混合釉を用いて製造された陶磁器であって、 2 3. Ceramics manufactured using mixed glaze made using computer color matching.

無色透明でないベース釉に白色着色剤を混合した第 1の混合釉の散乱係数 S _w_ ' と吸収係数 K_w，を決定し、 The scattering coefficient S _w _ 'and the absorption coefficient K _w , of the first mixed glaze in which the white colorant is mixed with the colorless and transparent base glaze are determined,

前記散乱係数 S _w' と吸収係数 K _w' の少なくとも一方を基準として、白色でない有色着色剤の吸収係数 K _P と散乱係数 S _P の相対値をそれぞれ決定し、前記吸収係数 K ， K _p および前記散乱係数 S _w，，を用いたコンビユータカラーマッチングを行なうことによって，所望の色を有する混台釉を調整するための着色剤の調合割合、または、所定の調合割合で生成される混合釉の色をザ' 測し、 Based on at least one of the scattering coefficient S _w ′ and the absorption coefficient K _w ′, the relative values of the absorption coefficient K _P and the scattering coefficient S _P of the colorant that is not white are determined, and the absorption coefficients K and K are determined. By performing a combi- ter color matching using _p and the scattering coefficients S _w ,, a mixture ratio of a colorant for adjusting a mixed-plate glaze having a desired color or a predetermined mixture ratio is generated. Measure the color of the mixed glaze,

前記予測された調合割合または前記予測された色を有するように前記ベース釉と前記白色着色剤と前記有色着色剤とを混合する，ことによって生成された混合釉で陶磁器の素地の少なくとも一部が覆われた，陶磁器。 Mixing the base glaze, the white colorant and the colored colorant so as to have the predicted blending ratio or the predicted color, wherein at least a portion of the ceramic body is Covered ceramic.

2 4 . コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する方法であって、 24. A method for determining the absorption and scattering coefficients of colorants used in computer color matching,

無色透明でない被着色物に基準となる顔料を混合した第 1の混合物の散乱係数 s と吸収係数 K_w' を決定する工程と、 A step of determining a scattering coefficient s and an absorption coefficient K _w ′ of a first mixture obtained by mixing a reference pigment with a colorless and non-transparent object

前記散乱係数 S と吸収係数 K の少なくとも一方を基準として、有色着色剤の吸収係数 K _P と散乱係数 S _P の相対値をそれぞれ決定する工程と、を備える方法。 The method comprising, based on the at least one of the scattering coefficient S and the absorption coefficient K, and determining the relative value of the absorption coefficient K _P and the scattering coefficient S _P output colored colorant, respectively, the.

2 5 . コンピュータカラーマッチングに使用される着色剤の吸収係数と散乱係数を決定する装置であって、 25. A device for determining the absorption coefficient and scattering coefficient of a colorant used in computer color matching,

無色透明でない被着色物に基となる顔料を混合した第 1の混合物の散乱係数 S _H' と吸収係数 K_w' を決定する手段と、 Means for determining a scattering coefficient S _H ′ and an absorption coefficient K _w ′ of a first mixture obtained by mixing a base pigment that is not colorless and transparent;

前記散乱係数 s _w' と吸収係数 K _w' の少なくとも一方を基準として、有色着色剤の吸収係数 Κ,, と散乱係数 S,> の相対値をそれぞれ決定する上程と、を備える方法。 Colored coloring based on at least one of the scattering coefficient s _w ′ and the absorption coefficient K _w ′ Determining the relative values of the absorption coefficient Κ ,, and the scattering coefficient S,> of the agent, respectively.

2 6. コンピュータカラーマッチングを利用して作成された混合釉を用いて製造された便器。 2 6. A toilet made using mixed glaze made using computer color matching.