JPH04274716A - Photoelectric chromaticity meter - Google Patents

Photoelectric chromaticity meter

Info

Publication number
JPH04274716A
JPH04274716A JP3625291A JP3625291A JPH04274716A JP H04274716 A JPH04274716 A JP H04274716A JP 3625291 A JP3625291 A JP 3625291A JP 3625291 A JP3625291 A JP 3625291A JP H04274716 A JPH04274716 A JP H04274716A
Authority
JP
Japan
Prior art keywords
correction
correction coefficient
measurement
photoelectric
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP3625291A
Other languages
Japanese (ja)
Other versions
JP3213952B2 (en
Inventor
Masahito Inaba
稲葉 政仁
Hitoshi Osugi
大椙 均
Masahiro Fukuda
福田 正浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Minolta Co Ltd
Original Assignee
Minolta Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Minolta Co Ltd filed Critical Minolta Co Ltd
Priority to JP03625291A priority Critical patent/JP3213952B2/en
Publication of JPH04274716A publication Critical patent/JPH04274716A/en
Application granted granted Critical
Publication of JP3213952B2 publication Critical patent/JP3213952B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Abstract

PURPOSE:To improve the correcting accuracy and to reduce the instrumental error by dividing the color space into a plurality of areas and performing correction for each area. CONSTITUTION:This photoelectric chromaticity meter is provided with a plurality of spectral sensors comprised of optical filters F1-F6 and photoelectric converting elements P1-P6. A plurality of reference samples in each of a plurality of areas divided in the color space are measured by the spectral sensors. A correcting coefficient is obtained so that the sum of the color difference between the measured value and the reference value of the reference sample is minimum for every area, which is in turn stored in a correcting coefficient memory means 8. Taking the correcting coefficient into consideration for tone measured value of an object to be measured, the measuring result is operated by an operating means 7.

Description

【発明の詳細な説明】[Detailed description of the invention]

【0001】0001

【産業上の利用分野】本発明は、光学フィルタと光電変
換素子からなる複数の分光センサを有し、測定対象から
の入射光を分光して該試料の色彩を測定する光電色彩計
に係り、特に分光感度に補正を施して測定精度の高い色
彩測定を行う光電色彩計に関する。
[Field of Industrial Application] The present invention relates to a photoelectric colorimeter that has a plurality of spectral sensors each consisting of an optical filter and a photoelectric conversion element, and measures the color of a sample by spectrally incoming light from an object to be measured. In particular, the present invention relates to a photoelectric colorimeter that performs color measurement with high measurement accuracy by correcting spectral sensitivity.

【0002】0002

【従来の技術】今日、試料の色彩値を容易、迅速に測定
できる色彩計が望まれている。ところで、色彩計の分光
感度については国際照明委員会(CIE)がその等色関
数を規定しているが、該等色関数に正確に一致した分光
感度を有する色彩計を製造することは極めて困難であり
、この不一致が測定時の絶対値誤差発生の主要因となっ
ている。また、同一機種間でも分光感度にばらつきが存
在するため、機器間誤差(以下、器差という)が生じて
いる。
BACKGROUND OF THE INVENTION Today, there is a need for a colorimeter that can easily and quickly measure the color value of a sample. By the way, the Commission Internationale de l'Eclairage (CIE) specifies color matching functions for the spectral sensitivity of colorimeters, but it is extremely difficult to manufacture colorimeters with spectral sensitivities that exactly match the color matching functions. This discrepancy is the main cause of absolute value errors during measurement. Furthermore, since there are variations in spectral sensitivity even between devices of the same model, inter-device errors (hereinafter referred to as instrumental errors) occur.

【0003】従来、試料の色彩値をより正確に測定する
ために、種々の方法が提案されている。例えば公知のT
Vカラーアナライザーでは、RGBの各分光分布に対し
その裾引き部分の干渉成分を除去すべく、RGBの各螢
光体毎に補正値を求めるようにしている。
Conventionally, various methods have been proposed to more accurately measure the color value of a sample. For example, the known T
In the V color analyzer, a correction value is obtained for each RGB phosphor in order to remove the interference component of the trailing portion of each RGB spectral distribution.

【0004】また、分光センサ出力x1,x2,x3,
……において、例えば出力x1の裾引き量を0にするた
めに他のセンサ出力x2,x3,……を利用するように
したもの、また三刺激値型光電色彩計のセンサ出力X=
x1+x2において、x1の代わりにkzを用い、X=
kz+x2として求めるようにしたものが知られている
[0004] Also, the spectral sensor outputs x1, x2, x3,
In ..., for example, in order to make the skirting amount of output x1 0, other sensor outputs x2, x3, ... are used, and the sensor output of tristimulus value type photoelectric colorimeter
In x1+x2, use kz instead of x1, and X=
It is known that the value is calculated as kz+x2.

【0005】更に、色度領域を複数に分割し、校正時に
領域毎の補正係数を求め、測定値に補正を加えるように
したものや、x,y,zそれぞれのセンサの分光感度に
他の2つのセンサの分光感度の定数倍を加減して補正を
行うものが提案されている(特開昭62−142239
号公報、特開平2−45718号公報)。
Furthermore, there are methods in which the chromaticity region is divided into a plurality of regions, correction coefficients are determined for each region during calibration, and correction is applied to the measured values, and spectral sensitivities of the x, y, and z sensors are adjusted by other factors. It has been proposed to perform correction by multiplying the spectral sensitivities of two sensors by a constant (Japanese Patent Laid-Open No. 142239/1983).
(Japanese Patent Application Laid-open No. 2-45718).

【0006】また、複数種類の試料について、ビデオプ
リンターの像の色度とCRT上の像の色度との色差の和
が最小になるように係数を決定するものが提供されてい
る(COLOR CORRECTION AND IN
TERPOLATION OF TELEVISION
 PICTURES FOR COMPACT PRI
NTING SYSTEMS)。
[0006] Furthermore, there has been proposed a method for determining coefficients for a plurality of types of samples so that the sum of the color differences between the chromaticity of an image on a video printer and the chromaticity of an image on a CRT is minimized (COLOR CORRECTION). AND IN
TERPOLATION OF TELEVISION
PICTURES FOR COMPACT PRI
NTING SYSTEMS).

【0007】[0007]

【発明が解決しようとする課題】上記の補正方法は一般
的に使用され、ある程度の測定精度が期待できるもので
あるが、色空間全体に亘っての平均的な効果を得るため
の補正であり、校正点も少ないため補正が粗く、特に特
定の領域の色のみを測定する場合等には十分な効果が得
られない。この場合、校正色を、例えば数千色程度用い
ると、所望の精度を期待することが出来るが、実際の生
産工程では、コスト、工数の点から現実的ではない。ま
た、上記従来の方法では前述の器差を減少する点は考慮
されていない。
[Problem to be Solved by the Invention] The above correction method is commonly used and can be expected to have a certain degree of measurement accuracy, but it is a correction to obtain an average effect over the entire color space. Since the number of calibration points is small, the correction is rough, and a sufficient effect cannot be obtained, especially when measuring only the color of a specific area. In this case, if several thousand proof colors are used, for example, desired accuracy can be expected, but this is not practical in the actual production process from the viewpoint of cost and man-hours. Furthermore, the above-mentioned conventional method does not take into consideration the reduction of the above-mentioned instrumental error.

【0008】更に、特開昭62−142239号公報記
載の色彩計では、分光感度自身を補正するものではない
ため、充分な効果が得られていない。また、特開平2−
45718号公報記載の色彩計では、分光感度自身を補
正するものではあるが、リニア系の色空間での補正であ
るため、色彩計で最も汎用的に利用されるノンリニア系
の色空間である等色差空間では必ずしもどの領域でも同
様に器差が低減されるとは限らないという問題がある。
Furthermore, the colorimeter described in JP-A-62-142239 does not correct the spectral sensitivity itself, and therefore does not provide a sufficient effect. Also, JP-A-2-
Although the colorimeter described in Publication No. 45718 corrects the spectral sensitivity itself, it is corrected in a linear color space, so it is a non-linear color space that is most commonly used in colorimeters. In the color difference space, there is a problem in that the instrumental error is not necessarily reduced in the same way in all regions.

【0009】本発明は、上記に鑑みてなされたもので、
色空間を複数の領域に分割し、各領域単位で補正を施す
ようにして、補正精度の向上及び器差の減少を図る光電
色彩計を提供することを目的とする。
The present invention has been made in view of the above, and
It is an object of the present invention to provide a photoelectric colorimeter that divides a color space into a plurality of regions and performs correction for each region to improve correction accuracy and reduce instrumental error.

【0010】0010

【課題を解決するための手段】本発明は、光学フィルタ
と光電変換素子からなる複数の分光センサを有する光電
色彩計において、複数の領域に分割された色空間の各領
域における複数の基準試料を上記分光センサで測定し、
得られた測定値と上記基準試料に付けられた基準値との
色差の和が領域毎に最小となるようにして求めた補正係
数を記憶する補正係数記憶手段と、測定対象を測定して
得られた測定値に上記補正係数を加味して測定結果を算
出する演算手段とを備えたものである。
[Means for Solving the Problems] The present invention provides a photoelectric colorimeter having a plurality of spectral sensors each consisting of an optical filter and a photoelectric conversion element, in which a plurality of reference samples are detected in each region of a color space divided into a plurality of regions. Measured with the above spectroscopic sensor,
a correction coefficient storage means for storing a correction coefficient obtained by minimizing the sum of color differences between the obtained measurement value and the reference value attached to the reference sample for each region; and calculating means for calculating a measurement result by adding the correction coefficient to the measured value.

【0011】また、各領域毎の前記基準試料の色度点を
代表する代表点を記憶する代表点記憶手段と、上記各代
表点と測定対象の測定点との距離をパラメータとする補
間関数を記憶する関数記憶手段とを備え、得られた前記
補正係数を上記補間関数を用いて補正し、この補正後の
補正係数を前記演算手段に導くようにしてもよい(請求
項2)。この場合、距離の2乗に反比例するパラメータ
を用いて補正係数の補間を行うことが好ましい(請求項
3)。
[0011] Also, representative point storage means for storing representative points representative of the chromaticity points of the reference sample for each area, and an interpolation function whose parameter is the distance between each of the representative points and the measurement point of the measurement object. The correction coefficient may be corrected using the interpolation function, and the corrected correction coefficient may be guided to the calculation means (Claim 2). In this case, it is preferable to interpolate the correction coefficient using a parameter that is inversely proportional to the square of the distance (claim 3).

【0012】0012

【作用】本発明によれば、複数の領域に分割された色空
間の各領域における複数の基準試料に予め基準値が値付
けされている。そして、補正を施すべき光電色彩計で上
記基準試料について測定を行い、得られた測定値と上記
基準値との差の和が領域毎に最小となるようにしたとき
の補正係数が求められ、この補正係数が補正係数記憶手
段に記憶される。この後、測定対象を測定して得られた
測定値に対して上記補正係数が加味され、この加味され
た後の測定値が測定結果とされる。
According to the present invention, reference values are assigned in advance to a plurality of reference samples in each region of a color space divided into a plurality of regions. Then, the photoelectric colorimeter to which correction is to be applied measures the reference sample, and a correction coefficient is determined so that the sum of the differences between the obtained measurement value and the reference value is minimized for each region, This correction coefficient is stored in the correction coefficient storage means. Thereafter, the correction coefficient is added to the measurement value obtained by measuring the measurement object, and the measurement value after this addition is taken as the measurement result.

【0013】また、請求項2記載の発明によれば、各領
域に対して代表点が決定されており、測定対象の測定点
と上記代表点とが一致しない場合でも、上記各代表点と
測定点間の距離をパラメータとする補間関数を用いるこ
とにより補正係数を補正すべく補間処理が施される。こ
のため、補正係数は連続的となる。そして、補正された
補正係数を測定値に加味することにより測定結果として
の測定値が得られる。更に、請求項3記載の発明によれ
ば、各代表点と測定点との距離の2乗に反比例するパラ
メータによって、補正係数が補間補正される。
Further, according to the invention as claimed in claim 2, a representative point is determined for each region, and even if the measurement point of the measurement object and the representative point do not match, the representative point and the measurement point are Interpolation processing is performed to correct the correction coefficient by using an interpolation function that uses the distance between points as a parameter. Therefore, the correction coefficient becomes continuous. Then, by adding the corrected correction coefficient to the measured value, a measured value as a measurement result is obtained. Furthermore, according to the third aspect of the invention, the correction coefficient is interpolated and corrected using a parameter that is inversely proportional to the square of the distance between each representative point and the measurement point.

【0014】[0014]

【実施例】図1は本発明にかかる光電色彩計の一実施例
を示すブロック図である。
DESCRIPTION OF THE PREFERRED EMBODIMENTS FIG. 1 is a block diagram showing an embodiment of a photoelectric colorimeter according to the present invention.

【0015】F1〜F6はフィルタで、F1〜F3は試
料1(測定対象)からの反射光のみが入射する位置ある
いはそのための構成を備えて配設され、F4〜F6は上
記試料1を照射する自然光に近い発光を行うキセノンラ
ンプ等の光源2からの光のみが入射する位置あるいはそ
のための構成を備えて配設されている。上記フィルタF
1〜F6の内、フィルタF1とF4、フィルタF2とF
5、フィルタF3とF6がそれぞれ対応する分光感度を
有するものである。P1〜P6はフィルタF1〜F6に
対応し、各フィルタ透過光の光強度に比例した電流を発
生する、例えばフォトダイオード等の光電変換素子で、
それぞれ対応するフィルタF1〜F6の後面に臨むよう
に配設されている。上記フィルタF1〜F6及び光電変
換素子P1〜P6との組み合せによって、前記CIE規
定の等色関数に出来るだけ近似した(ずれが小さい)分
光感度特性が得られる3個の分光センサが形成されてい
る。なお、光源2の光自身を測定する構成としているの
は、光源の発光色が試料の色彩判定に影響を及ぼしてい
るからであり、この影響分補正のための処理は後述する
第1演算部5により公知の方法を用いて行われる。
[0015] F1 to F6 are filters, F1 to F3 are arranged at positions where only the reflected light from the sample 1 (measurement object) enters, or have a configuration for that purpose, and F4 to F6 irradiate the sample 1. It is arranged at a position where only light from a light source 2 such as a xenon lamp that emits light similar to natural light is incident, or with a configuration for that purpose. Above filter F
Among 1 to F6, filters F1 and F4, filters F2 and F
5. Filters F3 and F6 each have corresponding spectral sensitivities. P1 to P6 correspond to filters F1 to F6, and are photoelectric conversion elements such as photodiodes, which generate a current proportional to the light intensity of light transmitted through each filter,
They are arranged so as to face the rear surfaces of the corresponding filters F1 to F6, respectively. The combination of the filters F1 to F6 and the photoelectric conversion elements P1 to P6 form three spectral sensors that can obtain spectral sensitivity characteristics that are as close as possible to the color matching function specified by the CIE (with a small deviation). . The reason why the light itself from the light source 2 is measured is because the emitted light color of the light source affects the color judgment of the sample, and the process for correcting this influence is performed by the first calculation unit, which will be described later. 5 using a known method.

【0016】3は光電変換素子P1〜P6からの各出力
電流値をそれぞれ電圧に変換する電流電圧変換回路であ
る。4は該電流電圧変換回路3の出力電圧をデジタル値
に変換するA/D変換回路である。5は入力される上記
デジタル値から上記の処理及び三刺激値としての測定値
X,Y,Zに変換する第1演算部で、該測定値X,Y,
Zは色情報等格納部6に格納されている色空間変換式等
に基づいて処理される。7は上記第1演算部5で得られ
た測定値X,Y,Zに対して、後述する補正及び色空間
の変換を施す第2演算部で、補正係数等格納部8に格納
されている前以て求められた補正係数等を用いて補正、
あるいは更に補間が行われるとともに所要の色空間への
変換を行うものである。9は表示部で、上記第2演算部
7で求めた測定結果としての測定値を表示するものであ
る。
Reference numeral 3 denotes a current-voltage conversion circuit that converts each output current value from the photoelectric conversion elements P1 to P6 into voltages. 4 is an A/D conversion circuit that converts the output voltage of the current-voltage conversion circuit 3 into a digital value. 5 is a first calculation unit that performs the above processing and converts the input digital values into measured values X, Y, and Z as tristimulus values;
Z is processed based on the color space conversion formula etc. stored in the color information storage section 6. Reference numeral 7 denotes a second calculation unit that performs correction and color space conversion, which will be described later, on the measured values X, Y, and Z obtained by the first calculation unit 5, which are stored in the correction coefficient storage unit 8. Correction using correction coefficients etc. determined in advance,
Alternatively, interpolation is further performed and conversion to a required color space is performed. Reference numeral 9 denotes a display section for displaying the measured value as the measurement result obtained by the second calculation section 7.

【0017】10は本光電色彩計の動作を統括的に制御
する制御部で、例えばマイクロコンピュータ等から構成
されている。この制御部10には測定に必要な各種デー
タをキー入力するキーボード11及び上記第1、2演算
部5,7での色空間の変換、その他必要なプログラムを
記憶している色空間変換システム等プログラム12が接
続されている。また、上記制御部10は上記キーボード
11等から入力された基準色を有する基準試料の基準値
データ等を用いて、後述する補正係数を算出するための
プログラム及び補正係数を補間するための補間関数及び
測定値に補正を加味するためのプログラムが内蔵ROM
等(不図示)に記憶されているものである。13は上記
制御部10により駆動制御され、前記電源2を点灯、消
灯させるものである。なお、補正係数を算出するプログ
ラムの換わりに、算出された補正係数を記憶しておいて
もよい。
Reference numeral 10 denotes a control section that comprehensively controls the operation of the photoelectric colorimeter, and is composed of, for example, a microcomputer. This control unit 10 includes a keyboard 11 for key-inputting various data necessary for measurement, a color space conversion system that stores color space conversion in the first and second calculation units 5 and 7, and other necessary programs, etc. Program 12 is connected. The control unit 10 also uses a reference value data of a reference sample having a reference color inputted from the keyboard 11 etc., and a program for calculating a correction coefficient, which will be described later, and an interpolation function for interpolating the correction coefficient. and a built-in ROM with programs for adding corrections to measured values.
etc. (not shown). Reference numeral 13 is driven and controlled by the control section 10 to turn on and off the power source 2. Note that instead of the program for calculating the correction coefficients, the calculated correction coefficients may be stored.

【0018】次に、色彩計における三刺激値X,Y,Z
、すなわち等色関数x(λ),y(λ),z(λ)に補
正を施す方法を以下に説明する。なお、説明は任意の等
色差空間で適用されるが、ここでは一例として、L*a
*b*空間を用いて説明する。
Next, the tristimulus values X, Y, Z on the colorimeter
, that is, a method of correcting the color matching functions x(λ), y(λ), and z(λ) will be described below. Note that the explanation applies to any isochromatic space, but here, as an example, L*a
This will be explained using *b* space.

【0019】今、ある色彩計について、最終的に得られ
る三刺激値を、Xc,Yc,Zcとするとき、Xc=X
+k11X+k12Y+k13Z    …(1)Yc
=Y+k21X+k22Y+k23Z    …(2)
Zc=Z+k31X+k32Y+k33Z    …(
3)となるような補正係数k11〜k33を求めること
にする。
Now, if the tristimulus values finally obtained for a certain colorimeter are Xc, Yc, and Zc, then Xc=X
+k11X+k12Y+k13Z...(1) Yc
=Y+k21X+k22Y+k23Z...(2)
Zc=Z+k31X+k32Y+k33Z...(
3) The correction coefficients k11 to k33 will be found.

【0020】先ず、基準試料についての基準値の値付け
が行われる。すなわち、a*b*平面をM個の領域、例
えば図2のように領域I〜VIIIのように8分割し、
このa*b*平面の領域Iに示すように、好ましくは満
遍なく点在するような色度点(図中、円形マークで示す
)を有するN個の基準試料を用意する。他の領域につい
ても同様に各領域内に満遍なく点在するような色度点を
有する所要個数の基準試料を用意する。そして、今、領
域Iについて説明すれば、該領域I内の各基準試料をマ
スター色彩計で測定し、それぞれ、 (L*1R,a*1R,b*1R),…………,(L*
NR,a*NR,b*NR) の値を得る。この値は基準値となり、各基準試料にそれ
ぞれ値付けされる。
First, a reference value for a reference sample is determined. That is, the a*b* plane is divided into M regions, for example, 8 regions I to VIII as shown in FIG.
As shown in region I of this a*b* plane, N reference samples having chromaticity points (indicated by circular marks in the figure) that are preferably evenly scattered are prepared. Similarly, for other regions, a required number of reference samples having chromaticity points evenly scattered within each region are prepared. Now, to explain region I, each reference sample in region I is measured with a master colorimeter, and (L*1R, a*1R, b*1R), ……, (L *
Obtain the values of NR, a*NR, b*NR). This value becomes a standard value, and is assigned a value to each standard sample.

【0021】以下、図3のフローチャートにより、例え
ば生産工程時において行われる上記補正係数を求める手
順を説明する。先ず、上記M個の領域をIの領域から順
次指定すべく,i=1にして(ステップ#1,#2)、
補正係数を求める手順が開始される。次に、補正が施さ
れるべき等色関数を持つ任意の色彩計で、領域I内のN
個の基準試料について測定を実行し、この測定により、
(L*1,a*1,b*1),…………,(L*N,a
*N,b*N) の測定値が得られる。(ステップ#3)。
Hereinafter, a procedure for determining the above-mentioned correction coefficient, which is carried out, for example, during the production process, will be explained with reference to the flowchart shown in FIG. First, in order to sequentially designate the M areas starting from area I, set i=1 (steps #1 and #2),
A procedure for determining a correction coefficient is started. Next, with any colorimeter whose color matching function is to be corrected, N
Measurements are carried out on reference samples of
(L*1,a*1,b*1),…………,(L*N,a
*N, b*N) measurement values are obtained. (Step #3).

【0022】測定が終了すると、ステップ#4で、以下
の補正係数算出のための演算が実行される。すなわち、
前記式(1)〜(3)内の補正係数k11〜k33に適
当な初期値を設定し、続いて測定の結果得られた測定値
(L*1,a*1,b*1),………,(L*N,a*
N,b*N)を(X1,Y1,Z1),………,(XN
,YN,ZN)で表される色空間の値に変換する。次に
、上記補正係数k11〜k33を用いて前記(1)〜(
3)式の補正演算を実行して(X1C,Y1C,Z1C
),………,(XNC,YNC,ZNC)を求め、更に
これらの値を(L*1C,a*1C,b*1C),……
…,(L*NC,a*NC,b*NC)に変換して元の
a*b*平面に戻す。
[0022] When the measurement is completed, the following computation for calculating the correction coefficient is executed in step #4. That is,
Appropriate initial values are set for the correction coefficients k11 to k33 in the above formulas (1) to (3), and then the measured values (L*1, a*1, b*1), . . . ..., (L*N,a*
N, b*N) as (X1, Y1, Z1), ......, (XN
, YN, ZN). Next, using the above correction coefficients k11 to k33, the above (1) to (
3) Execute the correction calculation of equation (X1C, Y1C, Z1C
), ......, (XNC, YNC, ZNC) and further convert these values into (L*1C, a*1C, b*1C),...
..., (L*NC, a*NC, b*NC) and return to the original a*b* plane.

【0023】上述の値付けされたN個の基準試料の色度
点と上記補正後の色度点との色差ΔE*abの和を、
The sum of the color differences ΔE*ab between the chromaticity points of the N reference samples valued above and the chromaticity points after the above correction is expressed as:


0024】
[
0024

【数1】[Math 1]

【0025】の(4)式から求める。It is determined from equation (4).

【0026】続いて、上記補正係数k11〜k33に他
の適当な値を入れて、上記と同様にして上記(4)式に
より色差ΔE*abの和を求める。以下、係る手順を複
数回繰り返し、求めた色差ΔE*abの和が最小となる
ときの補正係数k11〜k33を算出し、この補正係数
k11〜k33を領域Iの補正係数(k11)I〜(k
33)Iとして補正係数等格納部8に格納する(ステッ
プ#5)。iがMでなければ(ステップ#6でNO)、
iを1だけインクリメントして、上記領域Iに対して実
行した一連の手順(ステップ#3〜#5)を領域IIに
ついて実行する。このようにして順次同様な処理を領域
III,IV,…についても繰り返し、M個の領域全て
について行うことにより(ステップ#6でYES)、補
正係数(k11)I〜(k33)I,………,(k11
)VIII〜(k33)VIIIの補正係数等格納部8
への格納処理が終了する(ステップ#7)。
Next, other appropriate values are entered into the correction coefficients k11 to k33, and the sum of the color differences ΔE*ab is determined by the above equation (4) in the same manner as above. Hereinafter, such a procedure is repeated multiple times to calculate correction coefficients k11 to k33 when the sum of the obtained color differences ΔE*ab becomes the minimum, and these correction coefficients k11 to k33 are used as correction coefficients (k11)I to ( k
33) Store it in the correction coefficient etc. storage unit 8 as I (step #5). If i is not M (NO in step #6),
Increment i by 1, and execute the series of procedures (steps #3 to #5) executed for area I above for area II. In this way, the same process is repeated for areas III, IV, . . . , and by performing it for all M areas (YES in step #6), correction coefficients (k11)I to (k33)I, . ..., (k11
) VIII to (k33) VIII correction coefficient storage unit 8
The storage process ends (step #7).

【0027】かかる補正係数を求める一連の作業により
上記補正係数(k11)I〜(k33)I,………,(
k11)VIII〜(k33)VIIIが上述したよう
に色彩計本体内の補正係数等格納部8に予め格納される
こととなる。
Through a series of operations to obtain such correction coefficients, the above correction coefficients (k11)I to (k33)I, ......, (
k11) VIII to (k33) VIII are stored in advance in the correction coefficient storage section 8 within the colorimeter body as described above.

【0028】続いて、この色彩計を用いて実際の測定時
における測定値の補正について図4のフローチャートを
用いて説明する。先ず、測定対象に対して測定を行い、
三刺激値X,Y,Zを算出する(ステップ#11,#1
2)。次に、上記算出した三刺激値X,Y,ZをL*a
*b*の等色差空間の値に変換する(ステップ#13)
。第2演算部7は変換されたa*,b*の値より、その
試料が領域I〜VIIIのどの領域に属するかを判定す
る(ステップ#14)。属する領域が予め判断できると
きは、マニュアル入力も可能である。この第2演算部7
は判定された領域の補正係数(k11)i〜(k33)
iを前記補正係数等格納部8から読み込み(ステップ#
15)、上記の測定により得られた三刺激値X,Y,Z
に対して前記(1)〜(3)式による補正計算を施し、
補正後の三刺激値としてXC,YC,ZCを算出する(
ステップ#16)。続いて、この三刺激値XC,YC,
ZCはL*a*b*の等色差空間の値L*C,a*C,
b*Cに変換され(ステップ#17)、変換後測定結果
としてそのまま、あるいは他の表示形態で表示部9に表
示される(ステップ#18)。このような補正を各色彩
計で施すようにすることによって、複数台の色彩計で所
要の試料を測定したときの器差が等色差空間の全ての領
域で同様に縮小されることになる。
Next, correction of measured values during actual measurements using this colorimeter will be explained with reference to the flowchart of FIG. First, measure the object to be measured,
Calculate tristimulus values X, Y, Z (steps #11, #1
2). Next, the tristimulus values X, Y, and Z calculated above are calculated as L*a
Convert to a value in the isochromatic space of *b* (Step #13)
. The second calculation unit 7 determines to which region of regions I to VIII the sample belongs based on the converted values of a* and b* (step #14). Manual input is also possible when the area to which it belongs can be determined in advance. This second calculation unit 7
is the correction coefficient (k11)i~(k33) of the determined area
i is read from the correction coefficient etc. storage section 8 (step #
15), tristimulus values X, Y, Z obtained by the above measurements
Perform correction calculations according to formulas (1) to (3) above,
Calculate XC, YC, and ZC as tristimulus values after correction (
Step #16). Next, these tristimulus values XC, YC,
ZC is the value L*C, a*C, in the isochromatic space of L*a*b*
b*C (step #17), and displayed on the display unit 9 as the converted measurement result as it is or in another display form (step #18). By performing such correction in each colorimeter, the instrumental error when a required sample is measured with a plurality of colorimeters is similarly reduced in all regions of the isochromatic space.

【0029】次に、叙上の補正を色空間全体に亘つて連
続的に行う方法について説明する。
Next, a method of performing the above correction continuously over the entire color space will be explained.

【0030】この補正計算は色空間の各領域毎に求めら
れた前記補正係数を用いて行う。これは、領域の境界に
発生する不連続をなくし、かかる境界付近に色度点を持
つ試料の測定誤差をより低減させるようにするもので、
その手法は色度点間の距離の関数を用いて各領域間の補
正係数を補間するものである。
This correction calculation is performed using the correction coefficients determined for each region of the color space. This eliminates discontinuities that occur at the boundaries of regions and further reduces measurement errors for samples with chromaticity points near such boundaries.
This method interpolates correction coefficients between regions using a function of the distance between chromaticity points.

【0031】上記補正方法と同様、予め領域I〜VII
Iの補正係数(k11)I〜(k33)I,………,(
k11)VIII〜(k33)VIIIを求め、色彩計
本体内の補正係数等格納部8に格納しておく。また、同
時に、領域I〜VIIIを代表する色度点としての代表
点(a*ti,b*ti)(i=1〜8)も予め補正係
数等格納部8に格納しておく。この代表点は、特定の1
つの色度点の他、一領域内に存する基準点の位置平均で
もよい。図5はかかる場合の代表点を説明するための図
で、代表点は三角形マークで示している。
Similar to the above correction method, regions I to VII are
Correction coefficient of I (k11) I ~ (k33) I, ......, (
k11) VIII to (k33) VIII are determined and stored in the correction coefficient storage section 8 within the colorimeter body. At the same time, representative points (a*ti, b*ti) (i=1 to 8) as chromaticity points representing regions I to VIII are also stored in the correction coefficient storage section 8 in advance. This representative point is a specific point.
In addition to one chromaticity point, the positional average of reference points within one area may be used. FIG. 5 is a diagram for explaining the representative points in such a case, and the representative points are indicated by triangular marks.

【0032】図6はかかる補正係数の補間の手順を示す
フローチャートである。なお、この実施例では、領域数
Mを8として説明する。
FIG. 6 is a flowchart showing the procedure for interpolating the correction coefficients. In addition, in this example, the number M of regions will be described as eight.

【0033】先ず、測定対象に対して測定を行い(ステ
ップ#21)、三刺激値Xm,Ym,Zmを算出する(
ステップ#22)。次に、算出した三刺激値Xm,Ym
,ZmをL*ma*mb*mに色空間変換する(ステッ
プ#23)。そして、得られたa*b*平面における色
度点(a*m,b*m)と前述の各領域I〜VIIIの
代表点(a*ti,b*ti)(i=1〜8)との距離
diを(5)式より、
First, the measurement object is measured (step #21), and the tristimulus values Xm, Ym, and Zm are calculated (
Step #22). Next, the calculated tristimulus values Xm, Ym
, Zm to L*ma*mb*m (step #23). Then, the chromaticity points (a*m, b*m) on the obtained a*b* plane and the representative points (a*ti, b*ti) of each of the aforementioned regions I to VIII (i=1 to 8) From equation (5), the distance di from

【0034】[0034]

【数2】[Math 2]

【0035】と求める(ステップ#24)。この距離d
iと、格納されている各領域の補正係数(k11)i〜
(k33)i(i=1〜8)を用いて、下記(6)式に
より補間計算を行う(ステップ#25)。
[0035] is determined (step #24). This distance d
i and the correction coefficient (k11) of each stored area i~
(k33) Interpolation calculation is performed using the following equation (6) using i (i=1 to 8) (step #25).

【0036】[0036]

【数3】[Math 3]

【0037】上式で求めたk’lmが試料(a*m,b
*m)に最適な補正係数k’11〜k’33となる。こ
の補正係数を用いて、補正後の三刺激値Xm,Ym,Z
mを前記(1)〜(3)式から求める(ステップ#26
)。続いて、この三刺激値Xm,Ym,ZmはL*a*
b*の等色差空間の値L*m,a*m,b*mに変換さ
れ(ステップ#27)、該変換後測定結果として表示部
9に表示される(ステップ#28)。
k'lm obtained by the above formula is the sample (a*m, b
*m) are the optimum correction coefficients k'11 to k'33. Using this correction coefficient, the corrected tristimulus values Xm, Ym, Z
Find m from equations (1) to (3) above (step #26
). Next, these tristimulus values Xm, Ym, Zm are L*a*
b* is converted into values L*m, a*m, b*m in the isochromatic space (step #27), and displayed on the display unit 9 as the converted measurement result (step #28).

【0038】以上の補間手法は、一例としてL*a*b
*空間で領域数M=8としたものであるが、適用する等
色差空間、分割領域数及び分割方法は本実施例に限定さ
れるものではなく、目的等に応じて所要の形態を採用す
ることが出来る。また、補正係数を補間する際に用いる
距離の関数は、関数(1/x)のn乗(nは整数)に限
らず、例えば距離dをパラメータとする他の関数に拡張
することができる。なお、本実施例の関数において、n
=2とした場合、他の関数を用いた場合に比してより滑
らかな補間が行え、補正効果は大きい。
The above interpolation method uses, for example, L*a*b
*Although the number of regions M = 8 in the space, the applied isochromatic space, the number of divided regions, and the division method are not limited to this example, and a required form may be adopted depending on the purpose etc. I can do it. Furthermore, the distance function used when interpolating the correction coefficient is not limited to the function (1/x) to the nth power (n is an integer), but can be extended to other functions that use the distance d as a parameter, for example. Note that in the function of this example, n
When =2, smoother interpolation can be performed compared to when other functions are used, and the correction effect is large.

【0039】また、マスター色彩計での基準試料の測定
に際して行われる値付けは、通常付される値の他、マス
ター色彩計1台に限らず複数台の平均値でもよいし、ま
たCIE規定の2゜視野、10゜視野、等色関数と標準
の光、あるいは何等かの方法で測定された該基準試料の
分光反射率から計算上で求められた色度点でもよい。
In addition, the value assigned when measuring a reference sample with a master colorimeter may be the average value of not only one master colorimeter but multiple colorimeters, in addition to the normally assigned value, or the value determined by the CIE regulations. It may be a chromaticity point calculated from a 2° visual field, a 10° visual field, a color matching function and standard light, or the spectral reflectance of the reference sample measured by some method.

【0040】更に、補正係数の補間をする際の色空間の
各領域の代表点は、前記各領域の基準点の平均値に限ら
ず、該平均値に近い色度点を持つ基準試料の測定値その
ものでもよい。このようにすれば、代表値を求める処理
が容易となる。また、本発明が適用可能な色彩計につい
ては、三刺激値タイプに限定されず、複数のセンサを持
つ色測定装置に拡大することができる。
Furthermore, the representative point of each region of the color space when interpolating correction coefficients is not limited to the average value of the reference points of each region, but also the representative point of the reference sample having a chromaticity point close to the average value. It can also be the value itself. In this way, the process of determining the representative value becomes easy. Moreover, the colorimeter to which the present invention is applicable is not limited to the tristimulus value type, but can be expanded to color measuring devices having a plurality of sensors.

【0041】図7〜図12は測定精度を説明するための
図で、図7は補正なしの場合、図8は従来方法の補正が
施された場合、図9は本発明にかかる方法による補正が
施された場合を示している。
FIGS. 7 to 12 are diagrams for explaining measurement accuracy; FIG. 7 shows the case without correction, FIG. 8 shows the case when the conventional method has been corrected, and FIG. This shows the case where

【0042】図7〜図9において、100は分光反射率
R(λ)のBCRAタイルRED等の試料からの反射光
、200は等色関数の内の、例えばy(λ)のマスター
色彩計の分光感度曲線を示している。200aは各色彩
計の分光感度曲線のばらつきを示している。なお、図中
、分光感度曲線200はばらつきの分光感度曲線200
aの中の略中央に含まれており、見えていない。
In FIGS. 7 to 9, 100 is the reflected light from a sample such as BCRA tile RED with spectral reflectance R(λ), and 200 is the color matching function, for example, the master colorimeter's y(λ). The spectral sensitivity curve is shown. 200a shows variations in the spectral sensitivity curves of each colorimeter. In addition, in the figure, the spectral sensitivity curve 200 is the spectral sensitivity curve 200 with variations.
It is included approximately in the center of a and is not visible.

【0043】上記試料を測定した場合の測定結果は、図
7(2)〜図9(2)に示されており、測定値のばらつ
きを見るために、(y(λ)「ある色彩計」−y(λ)
「マスター色彩計」)×R(λ)という量を定義して、
マスター色彩計での測定値を基準(0%)とした場合の
各色彩計のばらつきを%表示で表している。上記図7(
2)〜図9(2)より、本発明に係る補正方法による場
合が最もばらつきが低減されていることが分かる。
The measurement results obtained when the above sample was measured are shown in FIGS. 7(2) to 9(2). In order to see the variation in the measured values, (y(λ) −y(λ)
Define the quantity "master colorimeter") × R (λ),
The variation of each colorimeter is expressed in % when the measurement value with the master colorimeter is used as a reference (0%). Figure 7 above (
2) to FIG. 9(2), it can be seen that the variation is reduced most when using the correction method according to the present invention.

【0044】図10〜図12は上記図7〜図9をL*a
*b*空間で色度分布したもので、それぞれ図7〜図9
の分光感度曲線200aのばらつきに対応している。本
発明に係る補正方法による図12は、図11に比して器
差が縮小されており、かつマスター色彩計による色度点
(図中、横線1と縦線2の交点)に接近している。すな
わち、器差が縮小されるとともに、絶対値精度も向上し
ていることが分かる。
FIGS. 10 to 12 are L*a
Chromaticity distribution in *b* space, shown in Figures 7 to 9, respectively.
This corresponds to variations in the spectral sensitivity curve 200a. In FIG. 12 obtained by the correction method according to the present invention, the instrumental error is reduced compared to FIG. There is. That is, it can be seen that the instrumental error is reduced and the absolute value accuracy is also improved.

【0045】[0045]

【発明の効果】以上説明したように、本発明によれば、
複数の領域に分割された色空間の各領域における複数の
基準試料を上記分光センサで測定し、得られた測定値と
上記基準試料に付けられた基準値との色差の和が領域毎
に最小となるようにして求めた補正係数を記憶する補正
係数記憶手段と、測定対象を測定して得られた測定値に
上記補正係数を加味して測定結果を算出する演算手段と
を備えたので、分光感度曲線自身が色域毎に補正され、
全領域に亘つて補正を行う従来の方法に比してより高精
度の補正が可能となり、特に測定する領域が限られてい
る場合には有効となるとともに、等色差空間での誤差を
縮小する補正のため視感に合った補正を施すことが出来
る。
[Effects of the Invention] As explained above, according to the present invention,
Multiple reference samples in each region of a color space divided into multiple regions are measured using the above spectroscopic sensor, and the sum of the color differences between the obtained measurement value and the reference value attached to the above reference sample is the minimum for each region. The present invention includes a correction coefficient storage means for storing the correction coefficient obtained as follows, and a calculation means for calculating the measurement result by adding the correction coefficient to the measurement value obtained by measuring the measurement target. The spectral sensitivity curve itself is corrected for each color gamut,
Compared to the conventional method that performs correction over the entire area, this method enables more accurate correction, is particularly effective when the area to be measured is limited, and reduces errors in isochromatic space. For correction purposes, it is possible to perform corrections that match visual perception.

【0046】また、補正係数の補間を行って連続にする
ため、色空間全体に亘つて不連続点のない補正が行え器
差の縮小がより一層図れる。更に、補正関数の補間を代
表点と測定点間の距離の2乗に反比例させることにより
、他の関数を用いた場合に比して滑らかな補間が行え、
補正効果を上げることが出来る。
Furthermore, since the correction coefficients are interpolated to make them continuous, correction without discontinuities can be performed over the entire color space, and instrumental differences can be further reduced. Furthermore, by making the interpolation of the correction function inversely proportional to the square of the distance between the representative point and the measurement point, smoother interpolation can be performed compared to when other functions are used.
It is possible to increase the correction effect.

【図面の簡単な説明】[Brief explanation of the drawing]

【図1】  本発明にかかる光電色彩計の一実施例を示
すブロック図である。
FIG. 1 is a block diagram showing an embodiment of a photoelectric colorimeter according to the present invention.

【図2】  a*b*平面の領域内に色度点を有するN
個の基準試料が示されるM個の領域、例えば領域I〜V
IIIのように8分割されたa*b*平面の図である。
[Figure 2] N with chromaticity points within the area of the a*b* plane
M regions in which reference samples are shown, e.g. regions I to V
It is a figure of the a*b* plane divided into 8 like III.

【図3】  例えば生産工程時において行われる補正係
数を求める手順を説明するフローチャートである。
FIG. 3 is a flowchart illustrating a procedure for determining a correction coefficient, which is performed, for example, during a production process.

【図4】  実際の測定時における測定値の補正につい
て説明するフローチャートである。
FIG. 4 is a flowchart illustrating correction of measured values during actual measurement.

【図5】  各領域における代表点を説明するためのa
*b*平面の図である。
[Figure 5] a for explaining representative points in each area
It is a diagram of the *b* plane.

【図6】  補正係数を補間して測定値を算出する手順
を説明するフローチャートである。
FIG. 6 is a flowchart illustrating a procedure for calculating a measured value by interpolating correction coefficients.

【図7】  測定精度を説明するための図で、補正なし
の場合を示すものである。
FIG. 7 is a diagram for explaining measurement accuracy, and shows a case without correction.

【図8】  測定精度を説明するための図で、従来方法
の補正が施された場合のものである。
FIG. 8 is a diagram for explaining measurement accuracy, and is a diagram after correction of the conventional method is applied.

【図9】  測定精度を説明するための図で、本発明に
かかる方法による補正が施された場合のものである。
FIG. 9 is a diagram for explaining measurement accuracy when correction is performed by the method according to the present invention.

【図10】  図7で示される測定値をL*a*b*空
間で色度分布した状態を示す図である。
10 is a diagram showing the chromaticity distribution of the measured values shown in FIG. 7 in L*a*b* space.

【図11】  図8で示される測定値をL*a*b*空
間で色度分布した状態を示す図である。
11 is a diagram showing a chromaticity distribution of the measured values shown in FIG. 8 in L*a*b* space.

【図12】  図9で示される測定値をL*a*b*空
間で色度分布した状態を示す図である。
12 is a diagram showing the chromaticity distribution of the measured values shown in FIG. 9 in L*a*b* space.

【符号の説明】[Explanation of symbols]

1  試料(測定対象) 2  光源 3  電流電圧変換回路 4  A/D変換回路 5  第1演算部 6  色情報等格納部 7  第2演算部 8  補正係数等格納部 9  表示部 10  制御部 11  キーボード 12  色空間変換システム等プログラム13  照明
回路 F1〜F6  フィルタ P1〜P6  光電変換素子
1 Sample (measurement target) 2 Light source 3 Current-voltage conversion circuit 4 A/D conversion circuit 5 First calculation unit 6 Color information storage unit 7 Second calculation unit 8 Correction coefficient storage unit 9 Display unit 10 Control unit 11 Keyboard 12 Color space conversion system, etc. program 13 Illumination circuits F1 to F6 Filters P1 to P6 Photoelectric conversion elements

Claims (3)

【特許請求の範囲】[Claims] 【請求項1】  光学フィルタと光電変換素子からなる
複数の分光センサを有する光電色彩計において、複数の
領域に分割された色空間の各領域における複数の基準試
料を上記分光センサで測定し、得られた測定値と基準試
料の上記基準値との色差の和が領域毎に最小となるよう
にして求めた補正係数を記憶する補正係数記憶手段と、
測定対象を測定して得られた測定値に上記補正係数を加
味して測定結果を算出する演算手段とを備えたことを特
徴とする光電色彩計。
Claim 1: A photoelectric colorimeter having a plurality of spectral sensors consisting of an optical filter and a photoelectric conversion element, in which a plurality of reference samples in each region of a color space divided into a plurality of regions are measured with the spectral sensor, correction coefficient storage means for storing a correction coefficient obtained by minimizing the sum of color differences between the measured value and the reference value of the reference sample for each region;
A photoelectric colorimeter comprising a calculation means for calculating a measurement result by adding the correction coefficient to a measurement value obtained by measuring a measurement object.
【請求項2】  請求項1記載の光電色彩計において、
各領域毎の前記基準試料の色度点を代表する代表点を記
憶する代表点記憶手段と、上記各代表点と測定対象の測
定点との距離をパラメータとする補間関数を記憶する関
数記憶手段とを備え、得られた前記補正係数を上記補間
関数を用いて補正し、この補正後の補正係数を前記演算
手段に導くようにしたことを特徴とする光電色彩計。
2. The photoelectric colorimeter according to claim 1, comprising:
representative point storage means for storing representative points representing the chromaticity points of the reference sample for each region; and function storage means for storing an interpolation function whose parameter is the distance between each of the representative points and the measurement point of the measurement object. A photoelectric colorimeter comprising the following: the obtained correction coefficient is corrected using the interpolation function, and the corrected correction coefficient is guided to the calculation means.
【請求項3】  前記補間関数は、前記代表点と測定対
象の測定点との距離をdとするとき、(1/d)2で表
されることを特徴とする請求項2記載の光電色彩計。
3. The photoelectric coloring device according to claim 2, wherein the interpolation function is expressed by (1/d)2, where d is a distance between the representative point and a measurement point of the measurement target. Total.
JP03625291A 1991-03-01 1991-03-01 Photoelectric colorimeter Expired - Lifetime JP3213952B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP03625291A JP3213952B2 (en) 1991-03-01 1991-03-01 Photoelectric colorimeter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP03625291A JP3213952B2 (en) 1991-03-01 1991-03-01 Photoelectric colorimeter

Publications (2)

Publication Number Publication Date
JPH04274716A true JPH04274716A (en) 1992-09-30
JP3213952B2 JP3213952B2 (en) 2001-10-02

Family

ID=12464584

Family Applications (1)

Application Number Title Priority Date Filing Date
JP03625291A Expired - Lifetime JP3213952B2 (en) 1991-03-01 1991-03-01 Photoelectric colorimeter

Country Status (1)

Country Link
JP (1) JP3213952B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011107114A (en) * 2009-10-20 2011-06-02 Hioki Ee Corp Photometric device

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011107114A (en) * 2009-10-20 2011-06-02 Hioki Ee Corp Photometric device

Also Published As

Publication number Publication date
JP3213952B2 (en) 2001-10-02

Similar Documents

Publication Publication Date Title
Bessell UBVRI passbands.
US4773761A (en) Photoelectric colorimeter
JP2007010364A (en) Spectroscope wavelength calibration method and spectroscope
JP5412694B2 (en) Light source measurement method and light source measurement system
US20040135750A1 (en) Transmission type display device and a method for controlling its display Colors
US4989982A (en) Spectral sensitivity correcting device in a photoelectric tristimulus colorimeter
JP6733667B2 (en) Spectral colorimeter and method for calculating spectral reflectance
US5363197A (en) Colorimeter for color displays
JP2007093477A (en) Method and apparatus of calibrating color measuring device, and color measuring method and device
WO2021118556A1 (en) Color calibration of display modules using a reduced number of display characteristic measurements
JPH04274716A (en) Photoelectric chromaticity meter
JPH04301530A (en) Dental color meter
JP2008157782A (en) Calibration method of tristimulus value direct reading type measuring device, calibration system, color measuring method and color measuring device
JPH0526731A (en) Photoelectric tint meter
JP3389678B2 (en) Color measuring device
WO2016181750A1 (en) Spectroscopic color measurement device and method for setting conversion rule
JPH0660846B2 (en) Photoelectric colorimeter with spectral sensitivity correction function
JP6787411B2 (en) How to calibrate a color filter type color luminance meter and a color filter type color luminance meter
JPH07111380B2 (en) Photoelectric colorimeter with spectral sensitivity correction function
US9311885B2 (en) Image color estimation method, image color estimation device, and image color estimation program
JPH0756461B2 (en) Color separation method
JP3208907B2 (en) Color value correction device
JPH0949765A (en) Device for measuring color
JPH08294009A (en) Image processing unit, its method and system
JPH09186892A (en) Read characteristic setting method of color image reader

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080727

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080727

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090727

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090727

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100727

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110727

Year of fee payment: 10

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110727

Year of fee payment: 10