# Patent

## Title
Color calibration method for imaging color measurement device

## Abstract

A color calibration method for an imaging color measurement device utilizes a detector array, a plurality of optical elements, and multiple instances of irradiation of the detector array for a single measurement. A flat-fielding correction error correction matrix of the imaging color measurement device for each instance of irradiation of the detector array is obtained prior to color calibration. The response for each instance of irradiation of the detector array is flat-fielded with the corresponding error matrix to obtain a flat-fielded, spectrally weighted irradiance response for each instance of irradiation of the detector array. An illuminant light source with known spectral output or chromaticity coordinates is measured to obtain an irradiance response of the imaging color measurement device for each instance of irradiation of the detector array. A color correction coefficient is calculated using the known spectral output and chromaticity coordinates of the light source and the corresponding flat-fielded irradiance response.

## Background
CROSS REFERENCE OF RELATED APPLICATION
(Not Applicable)
STATEMENT RE: FEDERALLY SPONSORED RESEARCH/DEVELOPMENT
(Not Applicable)
BACKGROUND OF THE INVENTION
The present invention relates generally to a color calibration method for imaging color measurement devices. More particularly, the present invention relates to a calibrating means for a color measurement system for the measurement of color, in quantifiable units.
It is often desired to know the color of light emitted by an illumination system, a light source, or a display system in order to characterize its performance. Traditionally, this has been accomplished by a calorimeter or a spectroradiometer, which measures the average color over an angular extent.
Currently, the spectroradiometer is the most accurate instrument for determining the color of light in quantifiable units. The spectroradiometer captures the entire spectrum of a visible light emitted by a light source. The light source includes any surface or volume that emits, transmits or reflects optical radiation. The measured spectral data may be mathematically analyzed to determine the color of the light defined by a standard color space. Commonly used color spaces for illumination systems are CIE (Commission Internationale de l'Eclairage) 1931 (x,y) chromaticity coordinates, CIE 1964 (u′,v′) chromaticity coordinates, and the L*, a*, b* color space. In principle, an infinite number of color space are possible. However, only a few are routinely used in industrial practice.
A diagram of the experimental set-up for display measurements made with a spectroradiometer are shown in FIG. 1. The spectroradiometer 12 captures light from some small portion of the display surface 10 (the light is emitted in the direction of the spectroradiometer) and acquires data of the spectral content of the light. A sample of four spectroradiometer measurement data sets is shown in FIG. 2. The vertical axis of the graph is normally recorded in W/nm/sr. The horizontal axis of the graph shows the wavelength of light, commonly in nm. The spectral power distribution of light is commonly written as function S(λ). The spectroradiometer is able to calculate the CIE 1931 tristimulus values by employing the following formulas defined by CIE: 



  
    X
    =
    
      
        ∫
        
          380
          ⁢
          
            
          
          ⁢
          nm
        
        
          780
          ⁢
          
            
          
          ⁢
          nm
        
      
      ⁢
      
        
          S
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          x
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          d
          λ
        
      
    
  
  ,
  
    Y
    =
    
      
        ∫
        
          380
          ⁢
          
            
          
          ⁢
          nm
        
        
          780
          ⁢
          
            
          
          ⁢
          nm
        
      
      ⁢
      
        
          S
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          y
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          d
          λ
        
      
    
  
  ,
  
    

  
  ⁢
  
    Z
    =
    
      
        ∫
        
          380
          ⁢
          
            
          
          ⁢
          nm
        
        
          780
          ⁢
          
            
          
          ⁢
          nm
        
      
      ⁢
      
        
          S
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          z
          ⁡
          
            (
            λ
            )
          
        
        ⁢
        
          d
          λ
        
      
    
  




In the above formulas, x(λ), y(λ) and z(λ) are the CIE 1931 2 degree standard observer color matching functions specifically defined by the CIE to serve as international definitions for the CIE 1931 (x,y) color space. The chromaticity coordinates are then determined by the following formulas: 



  
    x
    =
    
      X
      
        X
        +
        Y
        +
        Z
      
    
  
  ,
  
    y
    =
    
      Y
      
        X
        +
        Y
        +
        Z
      
    
  




If the spectroradiometer has been calibrated to accurately measure the spectral content of light, there is no further calibration required of a spectroradiometer to accurately measure the CIE 1931 tristimulus values and calculate the 1931 chromaticity coordinates accurately.
A calorimeter is often used as a color measurement device instead of a spectroradiometer because the cost of a calorimeter is substantially less than a spectroradiometer. Most calorimeters are tristimulus calorimeters, which means that they attempt to measure the X, Y and Z CIE 1931 chromaticity matching functions. A calorimeter normally consists of a lens to focus the light onto a detector (sometimes multiple detectors) and color filters disposed between the lens and the detector(s). The spectral response of the traditional tristimulus calorimeter closely approximates the x(λ), y(λ), and z(λ) tristimulus matching functions shown as FIG. 3. As shown in FIG. 3, the spectral response of a real calorimeter system is not an exact match to the spectral curves defined by the CIE 1931 color space.
In order to account for this spectral mismatch, the prior art calorimeters must be calibrated to measure the color more accurately. In general, both the magnitude and the shape of the spectral response curves do not match the CIE 1931 tristimulus values. The most common means of performing a calorimeter calibration has been to use a single point color calibration technique. During the calibration procedure, the calorimeter first measures uncorrected tristimulus values (Xuncorrected, Yuncorrected, Zuncorrected) for the spectrum of the standard light source, SIlluminant A(λ) Illuminant A defined by CIE is the most commonly used light source spectrum for calibration of calorimeter. Since the spectrum of the standard light source SIlluminant A(λ) is known, the chromaticity coordinates of the standard light source can be calculated. The chromaticity coordinates are designated as (xIlluminantA, yIlluminantA). Given these values, the calibration coefficients of the calorimeter (k1, k2, k3) can be calculated by the following equations. 



  
    
      
        
          x
          corrected
        
        =
        
          
            
              k1
              *
            
            ⁢
            
              X
              uncorrected
            
          
          
            
              
                k1
                *
              
              ⁢
              
                X
                uncorrected
              
            
            +
            
              
                k2
                *
              
              ⁢
              
                Y
                uncorrected
              
            
            +
            
              
                k3
                *
              
              ⁢
              
                Z
                uncorrected
              
            
          
        
      
    
  
  
    
      
        
          y
          corrected
        
        =
        
          
            
              k2
              *
            
            ⁢
            
              Y
              uncorrected
            
          
          
            
              
                k1
                *
              
              ⁢
              
                X
                uncorrected
              
            
            +
            
              
                k2
                *
              
              ⁢
              
                Y
                uncorrected
              
            
            +
            
              
                k3
                *
              
              ⁢
              
                Z
                uncorrected
              
            
          
        
      
    
  
  
    
      
        
          
            
              k1
              *
            
            ⁢
            
              X
              uncorrected
            
          
          +
          
            
              k2
              *
            
            ⁢
            
              Y
              uncorrected
            
          
          +
          
            
              k3
              *
            
            ⁢
            
              Z
              uncorrected
            
          
        
        =
        1
      
    
  
  
    
      where
    
  
  
    
      
        
          X
          corrected
        
        =
        
          
            
              k1
              *
            
            ⁢
            
              X
              uncorrected
            
          
          =
          
            k1
            ⁢
            
              
                ∫
                
                  380
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
                
                  780
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
              
              ⁢
              
                
                  
                    
                      S
                      
                        illuminant
                        ⁢
                        
                          
                        
                        ⁢
                        A
                      
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                  ·
                  
                    
                      x
                      colorimeter
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                
                ⁢
                
                  d
                  λ
                
              
            
          
        
      
    
  
  
    
      
        
          Y
          corrected
        
        =
        
          
            
              k2
              *
            
            ⁢
            
              Y
              uncorrected
            
          
          =
          
            k2
            ⁢
            
              
                ∫
                
                  380
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
                
                  780
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
              
              ⁢
              
                
                  
                    
                      S
                      
                        illuminant
                        ⁢
                        
                          
                        
                        ⁢
                        A
                      
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                  ·
                  
                    
                      y
                      colorimeter
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                
                ⁢
                
                  d
                  λ
                
              
            
          
        
      
    
  
  
    
      
        
          Z
          corrected
        
        =
        
          
            
              k3
              *
            
            ⁢
            
              Z
              uncorrected
            
          
          =
          
            k3
            ⁢
            
              
                ∫
                
                  380
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
                
                  780
                  ⁢
                  
                    
                  
                  ⁢
                  nm
                
              
              ⁢
              
                
                  
                    
                      S
                      
                        illuminant
                        ⁢
                        
                          
                        
                        ⁢
                        A
                      
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                  ·
                  
                    
                      z
                      colorimeter
                    
                    ⁡
                    
                      (
                      λ
                      )
                    
                  
                
                ⁢
                
                  d
                  λ
                
              
            
          
        
      
    
  




By storing the values k1, k2 and k3 and using them to calculate Xcorrected, Ycorrected and Zcorrected, the colorimeter is able to make accurate measurement of the chromaticity coordinates for calibration light source.

In FIG. 1, the measurements made with a spectroradiometer or a colorimeter extend over only a small area of the display device under test. Therefore, the traditional measurement comprises an average color or luminance measurement taken over a small portion of the light emitting area of the device under test. When a color measurement is performed on a display at multiple test locations, one has to reposition the prior art colorimeter or spectroradiometer from one point to another to measure the color of light at each test location.
BRIEF SUMMARY OF THE INVENTION
The present invention provides an improved color calibration method to measure color of light from multiple locations at one time. An imaging calorimeter comprising multiple individual detectors able to capture multiple light emitting areas is employed in this calibration method. Therefore, one can acquire much more data without moving the color measurement instrument (the imaging calorimeter).
In the color calibration method of the present invention, a flat-field calibration of the color measurement device is performed prior to color calibration.
For color measurements, multiple instances of irradiance of the array of detectors are performed for a single measurement. Optical elements in the imaging calorimeter modify the spectral power distribution of incoming light for each instance of irradiation of the detector array. The detector array response is flat-fielded for each instance of irradiation within a single color measurement. Color calibration coefficients are obtained for each instance of irradiation by selecting a subset of pixels from the different instances of irradiation, applying a flat-field correction matrix to the matrix of the pixel values, obtaining the total signal within the subset of pixels, and then calculating the appropriate color calibration coefficients based on the known spectrum or chromaticity coordinates of the light from the device under test.



BRIEF DESCRIPTION OF THE DRAWINGS
These, as well as other features of the present invention, will become more apparent upon reference to the drawings wherein:
FIG. 1 shows a traditional spot measurement technique of measuring a display device light and color output;
FIG. 2 shows a sample of four spectroradiometer measurement data sets using the spectroradiometer as shown in FIG. 1;
FIG. 3 shows both the spectral response of a calorimeter and the CIE 1931 chromaticity matching functions;
FIG. 4 shows an example of an imaging light and color measurement system provided in the invention;
FIG. 5 shows a typical imaging light and color measurement system camera;
FIG. 6 shows a non-uniform irradiance formed on a charge coupled device (image plane) by an imaging system and uniform luminance light source (object plane);
FIG. 7 shows two CCD images, one with and one without applying a flat-field correction matrix; and
FIG. 8 shows the x 1931 CIE chromaticity function and filter response for each hump of the function.



DETAILED DESCRIPTION OF THE INVENTION
As disclosed above, measurements made with a spectroradiometer or calorimeter as shown in FIG. 1 are taken over only a small area of the display device under test. Therefore, the traditional measurement comprises an average color or luminance measurement taken over an area at one location on the display surface. To reduce the measurement time and to analyze the color characteristics of the source of light, an imaging calorimeter is provided as the imaging light and color measurement instrument. The imaging calorimeter comprises a plurality, preferably thousands or millions of individual detectors. In one embodiment, a charge-coupled device (CCD) is cited as an example of the detector contained within the imaging calorimeter. FIG. 4 shows an imaging light and color measurement device 40 able to collect data for many points within the field of view thereof. As shown in FIG. 4, the light source 42 includes a plurality of light emitting areas designated as numeral references 1–9. All of the light emitting areas 1–9 are within the field of view (the solid circle in FIG. 4) of the imaging light and color measurement device 40. That is, the imaging light and color measurement device 40 is able to measure all the light-emitting areas 1–9 at once. The primary advantage the imaging color measurement system is the speed with which thousands or millions of data points are acquired. Thereby a complete analysis of the color characteristics of the source of light is achieved so that relative and absolute color differences can be efficiently measured.
The imaging calorimeter comprises a lens which images light onto an array of detectors (pixels), optical elements designed to modify the spectral power distribution of the light incident on the array of detectors (optionally designed to allow the spectral response of the system to match a specific color space), and software for camera control, calibration and data analysis. To measure color, the imaging calorimeter will have one or more instances of irradiation of the array of detectors within a single color measurement. Each instance of irradiation comprises taking an image with the imaging calorimeter with the optical elements set to modify the spectral power distribution of incoming light in a different manner. Specific embodiments of the optical elements used to modify the spectral power distribution of light incident on the array of detectors include using a plurality of color filters, using a single tunable color filter, using a plurality of tunable color filters, using an interferometer, using an interference grating, using a plurality of arrays of detectors for different spectra of incident light, using a spectrally tunable mirror or micromirror device, using a dispersive element such as a prism or holographic grating, using a spectrograph, or using a liquid crystal on silicon (LCOS) reflective element that controls the color of light reflected in a particular direction. In the preferred embodiment of the imaging calorimeter, a plurality of color filters are used as the optical elements designed to modify the spectral power distribution of the light incident on the array of detectors and a CCD is used as the array of detectors. A view of the imaging light and color measurement system is shown in FIG. 5. As shown in FIG. 5, in the CCD camera 50, light beams 52 from the image to be captured are focused onto the CCD 56 by an imaging lens 54. According to specific requirements, a filter wheel 58 that comprises various filters are disposed between the imaging lens 54 and the CCD 56.
For all of the embodiments of the imaging calorimeter mentioned above, it is known that the irradiance on the image plane such as the CCD 56 in FIG. 5 is not linearly related to the luminance in the object plane (the image to be captured). As shown in FIG. 6, for a perfect imaging system, the irradiance in the image plane falls off as cosine function of the field angle of the system.
In addition, for a real imaging system, factors such as lens vignetting can further affect the irradiance distribution in the image plane. This necessitates the CCD camera to be flat-fielded to accurately measure the luminance, radiance or other spectrally weighted radiance quantities in the object plane.
To accurately measure the luminance, radiance or other spectrally weighted radiance quantities in the object plane, the CCD camera must be “flat-fielded” for each instance of irradiation of the array of detectors within a single color measurement. In the flat-fielding process, the CCD is aimed at a uniform radiance light source, of which the actual or relative radiance values are known as a matrix L. The acquired image is in a form of a voltage value on each CCD pixel. The voltage value is proportional to the irradiation incident on each CCD pixel. The acquired image of the light source can thus be represented by a matrix M with a matrix of voltage values, i.e., the irradiance response. For example, when the CCD comprises a 512×512 matrix of pixels, the resultant image is then measured with a 512×512 matrix of voltage values. The degree that the measured irradiance at each pixel [i,j] deviates from the actual or relative radiance measured is represented by a correction coefficient K[i,j]. The coefficient for the whole matrix of pixels can thus be obtained by the relationship of M*K=L, which means that [i,j] the matrix element of L is equal to the [i,j] matrix element of M multiplied by the [i,j] matrix element of K. This should not be confused as standard matrix multiplication.
FIG. 6 shows a flat-field correction method for a camera with a detector array. In FIG. 7, the icon 70 illustrates the image of a light source with uniform luminance captured by a charge coupled device (CCD) without applying the correction matrix. The icon 72 shows the image of a light source with uniform luminance captured by a charge coupled device with the correction matrix applied.
The above flat fielding process is performed before the color calibration is performed. In the color calibration, each instance of irradiation of the CCD within a single color measurement must be flat-fielded. In the preferred embodiment, various filters are used in front of the as the optical elements designed to modify the spectral power distribution of the light incident on the CCD. In the preferred embodiment, four filters are used to allow the spectral response of the system to match the specific color space of the CIE (Commission Internationale de l'Eclairage) 1931 (x,y) chromaticity coordinates. The filters are used for the X1, X2, −Y and Z measurements respectively, whereas X1 is calculated using the hump of the x(λ) curve, designated by the curve H1 in FIG. 8, and X2 is calculated using the red hump of the x(λ) curve, designated by the curve H2 in FIG. 8. The two filters X1 and X2 are used together to additively obtain the measurement for X. For this embodiment, each instance of irradiation of the CCD through the different filters can be referred to as a different color image. The above flat-fielding process for each of the four color images allows calculation of the four corresponding error correction matrices. That is, for each filter, the flat-fielding calibration is performed to obtain the error correction matrices KX1, KX2, KY and KZ. The relationship between the 3 actual irradiance (and therefore actual measured voltage response of the CCD pixels) MX1, MX2, MY and MZ and the flat-fielded corrected spectrally weighted irradiance (power per unit area) EX1, EX2, EY and EZ incident on the array of CCD pixels for each of the four filters can be represented as:

EX1 flat-fielded=Kx1*MX1,

EX2 flat-fielded=Kx2*MX2,

EY flat-fielded=KY*MY,

EZ flat-fielded=KZ*MZ,

where Kx1 Kx2, KY, and Kz are the error correction matrices required for color calibration which can be obtained using the procedure as disclosed above.

In order to calculate the tristimulus values for each pixel of the CCD, appropriate scaling factors have to be determined by a calibration process for each of the tristimulus values. The matrices of tristimulus measurements for the CCD array measurements may be written as:

X1corrected=Cx1*EX1flat-fielded,

X2corrected=Cx2*EX2flat-fielded,

Ycorrected=CY*EYflat-fielded,

Zcorrected=CZ*EZflat-fielded,

where CX1, Cx2, CY, Cz are constants, called color correction coefficients, to be determined by the calibration process.

To determine the color correction coefficients, a standard illuminant light source with a known spectral content S(λ) and therefor known chromaticity coordinates (x,y) is measured. Four images are obtained by measuring the standard illuminant light source using the imaging colorimeter incorporating four filters. The same portion of each of the four images is selected as the color calibration area for each image, which amounts to some subset of the pixels in the matrices EX1, EX2, EY, EZ. This subset of the images is chosen to correspond to that portion of the images containing the light of known spectral content, and therefore the known chromaticity coordinates. The spectral characteristic of the calibration light source may be known by virtue of the fact that it is a calibrated standard illuminant, a special calibrated light source, or by measurement with a spectroradiometer. When a spectroradoimeter is used to measure the light source, it has to be ensured that the light measured by the spectroradiometer is identical to the light captured in the subset of the image matrices. If the subset of the image matrices is selected as EX1, EX2, EY, EZ and the subset of CCD pixels contains the same pixels for each image, and if the total (or the average) optical power per unit area (spectrally weighted irradiance) on each subset of pixels is designated as the scalar quantities ĒX1, ĒX2, ĒY, ĒZ, the color correction coefficients CX1, Cx2, CY, Cz can be derived from the following four equations. 



  
    
      
        
          x
          CalibrationStandard
        
        =
        
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  
                    X1
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
          
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  flat
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                Y
              
              *
              
                
                  E
                  _
                
                
                  Yflat
                  -
                  fielded
                
              
            
            +
            
              
                C
                Z
              
              *
              
                
                  E
                  _
                
                
                  Zflat
                  -
                  fielded
                
              
            
          
        
      
    
  
  
    
      
        
          y
          CalibrationStandard
        
        =
        
          
            
              C
              Y
            
            *
            
              
                E
                _
              
              
                
                  Y
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  
                    X1
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                Y
              
              *
              
                
                  E
                  _
                
                
                  Yflat
                  -
                  fielded
                
              
            
            +
            
              
                C
                Z
              
              *
              
                
                  E
                  _
                
                
                  Zflat
                  -
                  fielded
                
              
            
          
        
      
    
  
  
    
      
        
          
            
              C
              X1
            
            *
            
              
                E
                _
              
              
                
                  X1
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          +
          
            
              C
              X2
            
            *
            
              
                E
                _
              
              
                
                  X2
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          +
          
            
              C
              Y
            
            *
            
              
                E
                _
              
              
                Yflat
                -
                fielded
              
            
          
          +
          
            
              C
              Z
            
            *
            
              
                E
                _
              
              
                Zflat
                -
                fielded
              
            
          
        
        =
        1
      
    
  
  
    
      
        
          
            
              C
              X1
            
            *
            
              
                E
                _
              
              
                
                  X1
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          
            
              C
              X2
            
            *
            
              
                E
                _
              
              
                
                  X2
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
        
        =
        
          known
          ⁢
          
            
          
          ⁢
          constant
        
      
    
  




In the above equation, the known constant is a known value based on the CIE definition or spectroradiometer measurements of the calibration standard light source. Given the known parameters of 



  
    
      x
      
        Calibration
        ⁢
        
          
        
        ⁢
        Standard
      
    
    ,
    
      y
      
        Calibration
        ⁢
        
          
        
        ⁢
        Standard
      
    
    ,
    
      
        
          C
          X1
        
        *
        
          
            E
            _
          
          
            
              X1
              ⁢
              
                
              
              ⁢
              flat
            
            -
            fielded
          
        
      
      
        
          C
          X2
        
        *
        
          
            E
            _
          
          
            
              X2
              ⁢
              
                
              
              ⁢
              flat
            
            -
            fielded
          
        
      
    
  
  ⁢
  
    
  




and the average measured spectrally weighted irradiance of ĒX1flat-fielded, ĒX2flat-fielded, ĒYflat-fielded, ĒZflat-fielded, the color correction coefficients CX1, Cx2, CY, Cz can be calculated from the above four equations.

Having calculated the correction coefficients CX1, Cx2, CY, Cz in the calibration process using the above method, all subsequent color measurements are calculated using the correction coefficients. In other words, for the calibration of the imaging calorimeter, the color coordinates x and y are known values, and CX1, Cx2, CY, Cz are calculated based on data measured with the color measurement device as described above; however, after calibration, when using the system to take new color measurements, CX1, Cx2, CY, Cz are known values and the color coordinates x and y are calculated based on data measured with the color measurement device as described above.
According to the above, the color calibration method provided by the invention has the advantage of allowing accurate evaluation of light prior to calculating the color correction coefficients. From empirical data, it has been shown that the method of color calibration yields chromaticity coordinate accuracy of the measurements with the illuminant A of about ±0.001. With the new calibration technique utilizing the flat-fielding calibration, an accuracy of about ±0.0003 for both the x and y chromaticity coordinates of illuminant A is obtained.
Indeed, each of the features and embodiments described herein can be used by itself, or in combination with one or more of other features and embodiment. Thus, the invention is not limited by the illustrated embodiment but is to be defined by the following claims when read in the broadest reasonable manner to preserve the validity of the claims.

## Claims

1. A color calibration method for an imaging color measurement device, comprising:
selecting an array of detectors of optical radiation;
selecting a plurality of optical elements designed to modify spectral power distribution of light incident on the array of detectors, such that a plurality of instances of irradiance of the array of detectors in a single measurement is resulted;
obtaining an error correction matrix for the imaging color measurement device for each instance of irradiance of the array of detectors in a single color measurement;
measuring a light source with known spectral output or chromaticity coordinates to obtain an irradiance response matrix of the imaging color measurement device for each of instance of irradiation of the array of the detectors in a single color measurement;
applying the error correction matrix to the irradiance response matrix obtained by measuring the light source, so as to obtain a corrected irradiance response matrix for each instance of irradiation of the array of detectors in a single color measurement; and
calculating a color correction coefficient according to the known spectral output or chromaticity coordinate of the light source and the corrected response for each instance of irradiation of the array of detectors in a single color measurement.

2. The color calibration method according to claim 1, wherein the step of obtaining an error correction matrix of the imaging color measurement device for each instance of irradiation of the array of detectors in a single color measurement further comprises calculation of a flat-fielding correction matrix by:
aiming the imaging color measurement device at a light source with a known radiance value matrix L to acquire an image for each instance of irradiation of the array of detectors in a single color measurement, wherein a response M for each image is obtained;
calculating the flat-fielding correction matrix K according to the equation of M*K=L, wherein each element of L is equal to each corresponding element of M multiplied by each corresponding element of L; wherein
each element of the flat-fielding correction matrix K represent an error correction factor of a corresponding light point of the illuminant light source.

3. The color calibration method according to claim 1, wherein the step of selecting a plurality of optical elements designed to modify the spectral power distribution of light incident on the array of detectors includes selecting a set of optical filters designated to be used in conjunction with spectral response of a detector and imaging optics of the color calibration method to define a color space.

4. The color calibration method according to claim 1, wherein the step of selecting a plurality of optical elements designed to modify the spectral power distribution of light incident on the array of detectors includes selecting four filters, designated to be used in conjunction with the spectral response of the detector and imaging optics to match the tristimulus values X1, X2, Y and Z.

5. The color calibration method according to claim 4, wherein color correction coefficients CX1, Cx2, CY, Cz for each of the filters X1, X2, Y and Z are calculated according to the equations of: 



  
    
      
        
          x
          CalibrationStandard
        
        =
        
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  
                    X1
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
          
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  flat
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                Y
              
              *
              
                
                  E
                  _
                
                
                  Yflat
                  -
                  fielded
                
              
            
            +
            
              
                C
                Z
              
              *
              
                
                  E
                  _
                
                
                  Zflat
                  -
                  fielded
                
              
            
          
        
      
    
  
  
    
      
        
          y
          CalibrationStandard
        
        =
        
          
            
              C
              Y
            
            *
            
              
                E
                _
              
              
                
                  Y
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  
                    X1
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            +
            
              
                C
                Y
              
              *
              
                
                  E
                  _
                
                
                  Yflat
                  -
                  fielded
                
              
            
            +
            
              
                C
                Z
              
              *
              
                
                  E
                  _
                
                
                  Zflat
                  -
                  fielded
                
              
            
          
        
      
    
  
  
    
      
        
          
            
              C
              X1
            
            *
            
              
                E
                _
              
              
                
                  X1
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          +
          
            
              C
              X2
            
            *
            
              
                E
                _
              
              
                
                  X2
                  ⁢
                  
                    
                  
                  ⁢
                  flat
                
                -
                fielded
              
            
          
          +
          
            
              C
              Y
            
            *
            
              
                E
                _
              
              
                Yflat
                -
                fielded
              
            
          
          +
          
            
              C
              Z
            
            *
            
              
                E
                _
              
              
                Zflat
                -
                fielded
              
            
          
        
        =
        1
      
    
  
  
    
      
        
          
            
              
                C
                X1
              
              *
              
                
                  E
                  _
                
                
                  
                    X1
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
            
              
                C
                X2
              
              *
              
                
                  E
                  _
                
                
                  
                    X2
                    ⁢
                    
                      
                    
                    ⁢
                    flat
                  
                  -
                  fielded
                
              
            
          
          =
          
            known
            ⁢
            
              
            
            ⁢
            constant
          
        
        ,
      
    
  




wherein xColibriationStandard and yCalibrationStandard are the known chromaticity coordinates of the illuminate light source and ĒX1flat-fielded, ĒX2flat-fielded, ĒYflat-fielded, ĒZflat-fielded are total or average spectrally weighted irradiance on a subset of pixels of the array of detectors.

6. The color calibration method according to claim 1, wherein the step of selecting an array of detectors of optical radiation further comprising the step of selecting a charge coupled device.

7. The color calibration method according to claim 1, wherein the step of selecting an array of detectors of optical radiation further comprises the step of selecting a CMOS device.

8. The color calibration method according to claim 4, further comprising a step of obtaining tristimulus values X1, X2, Y and Z for the imaging color measurement device by multiplying the total or average spectrally weighted irradiance on a subset of pixels of the array of detectors ĒX1flat-fielded, ĒX2flat-fielded, ĒYflat-fielded, ĒZflat-fielded for each of the color filters with the corresponding color correction coefficient CX1, CX2, CY and CZ, respectively.

9. The color calibration method according to claim 1, wherein the imaging color measurement device includes an array of pixels, such that:
the error correction matrix for each instance of irradiation of the array of detectors in a single color measurement is in a form of a flat-fielding correction matrix;
the irradiance response for each of the color filters is in a form of a irradiance response matrix; and
the flat-fielded spectrally weighted irradiance response for each of the color filters is in a form of a flat-fielded spectrally weighted irradiance response matrix.

10. A color calibration method, comprising:
providing an imaging color measurement device able to measure a plurality of light emitting spots at once;
flat fielding the imaging color measurement device;
color calibrating the imaging color measurement device, further comprising the following steps:
measuring an illuminant light source with known spectral response or chromaticity coordinates to obtain a flat-fielded irradiance response and an average flat-fielded irradiance response;
calculating one or more color correction coefficients for the imaging color measurement device according to the known chromaticity coordinates of the illuminant light source, the flat-fielded irradiance response, and the average flat-fielded irradiance response; and
obtaining trismulus values of the imaging color measurement device by multiplying the flat-fielded irradiance response with the color correction coefficient.

11. The color calibration method according to claim 10, wherein the flat-fielding step is performed prior to the step of calculating a color correction coefficient of the imaging color measurement.

12. The color calibration method according to claim 10, further comprising using a plurality of color filters to measure the illuminant light source.

13. The color calibration method according to claim 12, wherein the imaging color measurement device is flat fielded for each of the color filters.

14. The color calibration method according to claim 12, wherein the step of obtaining the color correction coefficient includes calculating the color correction coefficient for each of the color filters.