WO2003029766A2 - Apparatus and method for measuring colour - Google Patents

Apparatus and method for measuring colour Download PDF

Info

Publication number
WO2003029766A2
WO2003029766A2 PCT/GB2002/004521 GB0204521W WO03029766A2 WO 2003029766 A2 WO2003029766 A2 WO 2003029766A2 GB 0204521 W GB0204521 W GB 0204521W WO 03029766 A2 WO03029766 A2 WO 03029766A2
Authority
WO
WIPO (PCT)
Prior art keywords
colour
values
enclosure
image
reflectance
Prior art date
Application number
PCT/GB2002/004521
Other languages
French (fr)
Other versions
WO2003029766A3 (en
Inventor
Ming Ronnier Luo
Chuangjun Li
Guihua Cui
Original Assignee
Digieye Plc
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
Priority claimed from GB0123810A external-priority patent/GB0123810D0/en
Application filed by Digieye Plc filed Critical Digieye Plc
Priority to US10/491,706 priority Critical patent/US20050018191A1/en
Priority to AU2002329470A priority patent/AU2002329470A1/en
Priority to EP02765098A priority patent/EP1436577A2/en
Publication of WO2003029766A2 publication Critical patent/WO2003029766A2/en
Publication of WO2003029766A3 publication Critical patent/WO2003029766A3/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/10Arrangements of light sources specially adapted for spectrometry or colorimetry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46Measurement of colour; Colour measuring devices, e.g. colorimeters

Definitions

  • the present invention relates to an apparatus and method for measuring colours, using a digital camera.
  • colour physics systems are widely used for colour quality control and recipe formulation purposes.
  • These systems generally include a computer and a colour measuring instrument, typically a spectrophotometer, which defines and measures colour in terms of its colorimetric values and spectral reflectance.
  • spectrophotometers are expensive and can only measure one colour at a time.
  • spectrophotometers are unable to measure the colours of curved surfaces or of very small areas.
  • a second area in which accurate colour characterisation is very important is in the field of graphic arts, where an original image must be reproduced on to a hard copy via a printing process.
  • colour management systems are frequently used for predicting the amounts of inks required to match the colours of the original image.
  • These systems require the measurement of a number of printed colour patches on a particular paper media via a colour measurement instrument, this process being called printer characterisation.
  • the colour measuring instruments can only measure one colour at a time.
  • the invention relates to the use of an apparatus including a digital camera for measuring colour.
  • a digital camera represents the colour of an object at each pixel within an image of the object in terms of red (R), green (G) and blue (B) signals, which may be expressed as follows:
  • S( ⁇ ) is the spectral power distribution of the illuminant
  • R( ⁇ ) is the reflectance function of a physical object captured by a camera at a pixel within the image (and is between 0 and 1)
  • the k factor is a normalising factor to make G equal to 100 for a reference white.
  • the colour of the object at each pixel may alternatively be expressed in terms of standard tristimulus (X, Y, Z) values, as defined by the CIE (International Commission on Illumination).
  • X, Y, Z standard tristimulus values
  • CIE International Commission on Illumination
  • the x,y,z are the CIE 1931 or 1964 standard colorimetric observer functions, also known as colour matching functions (CMF), which define the amounts of reference red, green and blue lights in order to match a monochromatic light in the visible range.
  • CMF colour matching functions
  • the k factor is a normalising factor to make Y equal to 100 for a reference white.
  • the reflectance function defines the extent to which light at each visible wavelength is reflected by the object and therefore provides an accurate characterisation of the colour.
  • any particular set of R, G, B or X, Y, Z values could define any of a large number of different reflectance functions.
  • the corresponding colours of these reflectance functions will produce the same colour under a reference light source, such as daylight.
  • a reference light source such as daylight.
  • an inappropriate reflectance function is derived from the camera R, G, B values, the colour of the object at the pixel in question may be defined in such a way that, for example, it appears to be a very different colour under a different light source, for example a tungsten light.
  • an apparatus for measuring colours of an object including: an enclosure for receiving the object; illumination means for illuminating the object within the enclosure; a digital camera for capturing an image of the object; a computer connected to the digital camera, for processing information relating to the image of the object; and display means for displaying information relating to the image of the object.
  • digital camera it should be taken to be interchangeable with or to include other digital imaging means such as a colour scanner.
  • the enclosure may include means for mounting an object therein such that its position may be altered.
  • These means may include a tiltable table for receiving the , object.
  • the tiltable table is controllable by the computer.
  • the illumination means are located within the enclosure.
  • the illumination means may include diffusing means for providing a diffuse light throughout the enclosure.
  • the illumination means includes a plurality of different light sources for providing respectively different illuminations for the object.
  • One or more of the light sources may be adjustable to adjust the level of the illumination or the direction of the illumination.
  • the light sources may be controllable by the computer.
  • the digital camera is mounted on the enclosure and is directed into the enclosure for taking an image of the object within the enclosure.
  • the camera is mounted such that its position relative to the enclosure may be varied.
  • the location and/or the angle of the digital camera may be varied.
  • the camera may be adjusted by the computer.
  • the display means may include a video display unit, which may include a cathode ray tube (CRT).
  • a video display unit which may include a cathode ray tube (CRT).
  • CRT cathode ray tube
  • a method for measuring colours of an object including the steps of: locating the object in an enclosure; illuminating the object within the enclosure; using a digital camera to capture an image of the object within the enclosure; using a computer to process information relating to the image of the object; and displaying selected information relating to the image of the object.
  • the method may include the step of illuminating the object with a number of respectively different light sources.
  • the light may be diffuse.
  • the light sources may be controlled by the computer.
  • the digital camera may also be controlled by the computer.
  • the method preferably includes the step of calibrating the digital camera, to transform its red, green, blue (R, G, B) signals into standard X, Y, Z values.
  • the calibration step may include taking an image of a reference chart under one or more of the light sources and comparing the camera responses for each known colour within the reference chart with the standard X, Y, Z responses for that colour.
  • the relationship between the measured R, G, B values and the predicted X, Y, Z values is preferably represented as follows:
  • X , Y , Z are the measured tristimulus values and X , Y , Z are the predicted tristimulus values.
  • the method may include the step of predicting a reflectance function for a pixel or group of pixels within the image of the object.
  • P Wr, where P includes known X , Y , Z values, W is a known weight matrix derived from the product of an illuminant function and the CIE x, y, z colour matching functions, W is the transposition of the matrix W and r is an unknown n component column vector representing reflectance function defined by:
  • the camera is initially calibrated so that measured R, G, B values can be trans ormed to predicted X , Y , Z values.
  • the X , Y , Z values may then be used to predict the reflectance functions.
  • R, G, B values may be used to predict reflectance functions directly using the following steps:
  • the weighting factors may be predetermined and are preferably calculated empirically.
  • n is at least 10. Most preferably n is at least 16, and n may be
  • the smoothness is defined by determining the following:
  • G is an (n-l) x (n) matrix defined by the following:
  • o is an ⁇ component zero vector and e is an n component column vector where all the elements are unity (equal one).
  • the colour constancy of the reflectance vector is calculated as follows:- compute tristimulus X, Y, Z values (denoted P R ) using the reflectance vector, under a reference illuminant; compute tristimulus X, Y, Z values (denoted P_) using the reflectance vector, under a test illuminant; using a chromatic adaptation transform, transfer P to a corresponding colour denoted by P TC under the reference illuminant; compute the difference ⁇ E between P ⁇ c and P R ; and define the colour inconstancy index (CON) as ⁇ E .
  • a plurality / of test illuminants may be used such that the colour j inconstancy index is defined as ⁇ ⁇ .AE j where ⁇ ⁇ is a weighting factor
  • the reference illuminant is preferably D65, which represents daylight.
  • the smoothness weighting function ⁇ may be set to zero, such that the reflectance is generated with the least colour inconstancy.
  • the colour constancy weighting factors ⁇ ⁇ . may alternatively be set to zero, such that the reflectance vector has smoothness only.
  • ⁇ and ⁇ . are set such that the method generates a reflectance function having a high degree of smoothness and colour constancy.
  • the values of ⁇ and ⁇ . may be determined by trial and error.
  • the method further includes the step of providing an indication of an appearance of texture within a selected area of the object.
  • the method may include the steps of: determining an average colour value for the whole of the selected area; and determining a difference value at each pixel within the image of the selected area, the difference value representing the difference between the measured colour at that pixel and the average colour value for the selected area.
  • the selected area has a substantially uniform colour.
  • the difference value may be a value A Y which represents the difference between the tristimulus value Y at that pixel and the average Y for the selected area.
  • the difference value may also include a value ⁇ X which represents the difference between the tristimulus value X at that pixel and the average X for the selected area and/or a value ⁇ Z which represents the difference between the tristimulus value Z at that pixel and the average Z for the selected area.
  • the texture of the selected area may be represented by an image comprising the difference values for all the respective pixels within the selected area.
  • the method may further include the step of simulating the texture of a selected area of an object, for example in an alternative, selected colour.
  • the method may include the step of:
  • X, Y, Z values for each pixel l,m may be transformed to:
  • Fig. 1 is a diagrammatic overview of an apparatus according to the invention
  • Fig. 2 is a diagrammatic sectional view of an illumination box for use with the apparatus of Fig. 1.
  • an apparatus includes an illumination box 10 in which an object 18 to be observed may be placed.
  • a digital camera 12 is located towards the top of the illumination box 10 so that the digital camera 12 may take a picture of the object 18 enclosed in the illumination box 10.
  • the digital camera 12 is connected to a computer 14 provided with a video display unit (VDU) 16, which includes a colour sensor 30.
  • VDU video display unit
  • the illumination box 10 is provided with light sources 20 which are able to provide a very carefully controlled illumination within the box 10.
  • Each light source includes a lamp 21 and a diffuser 22, through which the light passes in order to provide uniform, diffuse light within the illumination box 10.
  • the inner surfaces of the illumination box are of a highly diffusive material coated with a matt paint for ensuring that the light within the box is diffused and uniform.
  • the light sources are able to provide a variety of different illuminations within the illumination box 10, including: D65, which represents daylight; tungsten light; and lights equivalent to those used in various department stores, etc.
  • D65 represents daylight
  • tungsten light and lights equivalent to those used in various department stores, etc.
  • the illumination is fully characterised, i.e., the amounts of the various different wavelengths of light are known.
  • the illumination box 10 includes a tiltable table 24 on which the object 18 may be placed. This allows the angle of the object to be adjusted, allowing different parts of the object to be viewed by the camera.
  • the camera 12 is mounted on a slider 26, which allows the camera to move up and down as viewed in Fig. 2. This allows the lens of the camera to be brought closer to and further away from the object, as desired.
  • the orientation of the camera may also be adjusted.
  • the light sources 20, the digital camera 12 and its slider 26 and the tiltable table 24 may all be controllable automatically from the computer 14. Alternatively, control may be effected from control buttons on the illumination box or directly by manual manipulation.
  • the digital camera 12 is connected to the computer 14 which is in turn connected to the VDU 16.
  • the image taken by the camera 12 is processed by the computer 14 and all or selected parts of that image or colours or textures within that image may be displayed on the VDU and analysed in various ways. This is described in more detail hereinafter.
  • the digital camera describes the colour of the object at each pixel in terms of red (R), green (G) and blue (B) signals, which are expressed in the following equations
  • S( ⁇ ) is the spectral power distribution of the illuminant. Given that the object is illuminated within the illumination box 10 by the light sources 20, the spectral power distribution of any illuminant used is known.
  • R( ⁇ ) is the reflectance function of the object at the pixel in question (which is unknown) and r,g,b are the spectral sensitivities of the digital camera, i.e., the responses of the charge coupled device (CCD) sensors used by the camera.
  • CCD charge coupled device
  • CIE 1931 or 1964 standard colorimetric observer functions also known as colour matching functions (CMF), which define the amounts of reference red, green and blue lights in order to match a monochromatic light in the visible range.
  • the k factor in equation (2) is a normalising factor to make Y equal to 100 for a reference white.
  • R, G, B values captured by the digital camera may be transformed into X, Y, Z values, it is desirable to calibrate the digital camera before the apparatus is used to measure colours of the object 18. This is done each time the camera is switched on or whenever the light source or camera setting is altered.
  • the camera is calibrated by using a standard colour chart, such as a GretagMacbeth ColorChecker Chart or Digital Chart.
  • the chart is placed in the illumination box 10 and the camera 12 takes an image of the chart.
  • the X, Y, Z values are known.
  • the values are obtained either from the suppliers of the chart or by measuring the colours in the chart by using a colour measuring instrument.
  • a polynomial modelling technique may be used to transform from the camera R, G, B values to X, Y, Z values.
  • each pixel represented by R, G, B values is transformed using the following equation to predict X , Y , Z values, these being the X, Y, Z values at a particular pixel:
  • the coefficients in the 3 by 11 matrix M may be obtained via an optimisation method based on a least squares technique.
  • the digital camera may be calibrated such that its R, G, B readings for any particular colour may be accurately transformed into standard X, Y, Z values.
  • VDU 16 It is also necessary to characterise the VDU 16. This may be carried out using known techniques, such as are described in Berns R.S. et al, CRT colorimetry, Part I and II at Col, Res Appn, 1993.
  • a sample object may be placed into the illumination box 10.
  • the digital camera is controlled directly or via the computer 14, to take an image of the object 18.
  • the image may be displayed on the VDU 16.
  • the apparatus preferably predicts the' reflectance function of the object at each pixel. . This ensures that the colour of the object is realistically characterised and can be displayed accurably on the VDU, and reproduced on other objects if required.
  • W is a n x 3 matrix called the weight matrix, derived from the illuminant function and the sensors of the camera for equation (1), or from the illuminant used, and the colour matching functions for equation (2), W is the transposition of the matrix W and r is the unknown n component column vector (the reflectance vector) representing unknown reflectance function given by:
  • the 3-component column vector p consists of either the camera responses R, G and B for the equation (1), or the CIE tristimulus values X, Y and Z for the equation (2).
  • o is a ⁇ -component zero vector and e is a rz-component vector where all the elements are unity (equal one).
  • Some fluorescent materials have reflectances of more than 1, but this method is not generally applicable to characterising the colours of such materials.
  • the preferred method used with the present invention recovers the reflectance vector r satisfying equation (3) by knowing all the other parameters or functions in equations (1) and (2).
  • the method uses a numerical approach and generates a reflectance vector r defined by equation (4) that is smooth and has a high degree of colour constancy.
  • colour constant products i.e., the colour appearance of the goods will not be changed when viewed under a wide range of light sources such as daylight, store lighting, tungsten.
  • a smoothness constraint condition is defined as follows:
  • G is an (n-1) x n matrix referred to as the "smooth operator”, and defined by the following:
  • the chromatic transform CMCCAT97 is described in the following paper: “M R Luo and R W G Hunt, A chromatic adaptation transform and a colour inconstancy Index, Color Res Appn, 1998".
  • the colour difference formula is described in "M R Luo, G Cui and B Rigg, The development of the CIE 2000 colour difference Formula: CIEDE2000, Color Res Appn, 2001”.
  • the reference and test illuminants are provided by the illumination box 10 and are thus fully characterised, allowing the above calculations to be carried out accurately.
  • the method may be summarised as follows:
  • the above method If the smoothness weighting factor is set to 0, then the above method generates the reflectance with the least colour inconstancy. However, the reflectance vector r could be too fluctuated to be realistic. At the other extreme, if the weighting factors ⁇ s are all set to be zero, then the above method produces a reflectance vector r with smoothness only. By choosing appropriate weighting factors, ⁇ and ⁇ , the above method generates reflectances with smoothness and a high degree of colour constancy.
  • the weight matrix W should be known from the camera characterisation carried out before the apparatus is used to measure the colours of the object 18.
  • the above described method for predicting a reflectance function from the digital camera's red, green and blue signals results in a reflectance function which is smooth and colour constant across a number of illuminants.
  • the apparatus is able to characterise and reproduce a colour of the object 18 very realistically and in such a way that the colour is relatively uniform in appearance under various different illuminants.
  • An image of the existing object 18 is taken using the digital camera 12 and a particular area of uniform colour to be analysed is isolated from the background using known software.
  • the R, G, B values are transformed to standardised X, Y, Z values.
  • Average colour values X, Y, Z are calculated, these being the mean X, Y, Z values for the whole selected area of colour.
  • a difference value ⁇ Y is calculated, ⁇ Y being equal to the difference between the Y value at the pixel in question and the average Y value,
  • ⁇ Y l m equals Y lm -Y, where l,m represents a particular pixel.
  • the computer calculates ⁇ Y values at each pixel within the selected area of colour in the image. Because the colour of the area is uniform, the variations in the measured Y values from the average Y value must represent textural effects. Thus the computer can create a "texture profile" for the area of colour, the profile being substantially independent of the colour of the area.
  • ⁇ Y values are stored for each pixel in the selected area, providing the texture profile, this may be used to simulate a similar texture in a different colour. This is carried out as follows.
  • the new colour is measured or theoretical colour values provided.
  • the X, Y, Z values are transformed to x, y, Y, where
  • the X, Y, Z colour space is not very uniform, including very small areas of blue and very large areas of green.
  • the above transform transfers the colour to x, y, Y space in which the various colours are more uniformly represented.
  • the x, y values remain the same and the 7 value is replaced with a value Y l m , for a pixel l,m.
  • X, Y, Z values for each pixel l,m may be transformed to:
  • t varies with Y but there are different functions of t against Y for different materials, with the relationship between f and Y depending upon the coarseness of the material.
  • the appropriate values of t may be calculated empirically.
  • the illumination box 10 allows objects to be viewed in controlled conditions under a variety of accurately characterised lights. This, preferably together with the novel method for predicting reflectance functions, enables colours to be characterised in such a way that they are predictable and realistically characterised under all lights.
  • the apparatus and method also provide additional functions such as the ability to superimpose a texture of one fabric on to a different coloured fabric.

Abstract

An apparatus and method for measuring colours of an object includes: an enclosure for receiving the object; illumination means for illuminating the object within the enclosure; a digital camera for capturing an image of the object; a computer connected to the digital camera, for processing information relating to the image of the object; and display means for displaying information relating to the image of the object. The enclosure may include means for mounting an object therein such that its position may be altered. These means may include a tiltable table for receiving the object, the tiltable table being controllable by the computer. The illumination means are preferably located within the enclosure, and may include diffusing means for providing a diffuse light throughout the enclosure. The illumination means may include a plurality of different light sources for providing respectively different illuminations for the object, one or more of the light sources may be adjustable to adjust the level of the illumination or the direction of the illumination. The light sources may be controllable by the computer.

Description

APPARATUS AND METHOD FOR MEASURING COLOUR
The present invention relates to an apparatus and method for measuring colours, using a digital camera.
There are many applications in which the accurate measurement of colour is very important. Firstly, in the surface colour industries such as textiles, leather, paint, plastics, packaging, printing, paper and food, colour physics systems are widely used for colour quality control and recipe formulation purposes. These systems generally include a computer and a colour measuring instrument, typically a spectrophotometer, which defines and measures colour in terms of its colorimetric values and spectral reflectance. However, spectrophotometers are expensive and can only measure one colour at a time. In addition, spectrophotometers are unable to measure the colours of curved surfaces or of very small areas.
A second area in which accurate colour characterisation is very important is in the field of graphic arts, where an original image must be reproduced on to a hard copy via a printing process. Presently, colour management systems are frequently used for predicting the amounts of inks required to match the colours of the original image. These systems require the measurement of a number of printed colour patches on a particular paper media via a colour measurement instrument, this process being called printer characterisation. As mentioned above, the colour measuring instruments can only measure one colour at a time.
Finally, the accurate measurement of colour is very important in the area of professional photography, for example for mail order catalogues, internet shopping, etc. There is a need to quickly capture images with high colour fidelity and high image quality over time.
The invention relates to the use of an apparatus including a digital camera for measuring colour. A digital camera represents the colour of an object at each pixel within an image of the object in terms of red (R), green (G) and blue (B) signals, which may be expressed as follows:
b _
R = k,\S(λ)r(λ)R(λ) λ a
Figure imgf000004_0001
a where S(λ) is the spectral power distribution of the illuminant, R(λ) is the reflectance function of a physical object captured by a camera at a pixel within the image (and is between 0 and 1) and r , g, b are the responses of the CCD sensors used by the camera. All the above functions are defined within the visible range, typically between a=400 and b=700 nm. The k factor is a normalising factor to make G equal to 100 for a reference white.
The colour of the object at each pixel may alternatively be expressed in terms of standard tristimulus (X, Y, Z) values, as defined by the CIE (International Commission on Illumination). The tristimulus values are defined as follows:
b _
X = k\S{λ)x{λ)R{λ)dλ b b _
Y = k\S{λ)y{λ)R{λ)dλ b
Figure imgf000004_0002
where all the other functions were defined. The x,y,z are the CIE 1931 or 1964 standard colorimetric observer functions, also known as colour matching functions (CMF), which define the amounts of reference red, green and blue lights in order to match a monochromatic light in the visible range. The k factor is a normalising factor to make Y equal to 100 for a reference white. In order to provide full colour information about the object, it is desirable to predict the colorimetric values or reflectance function of the object at each pixel, from the RGB or X, Y, Z values. The reflectance function defines the extent to which light at each visible wavelength is reflected by the object and therefore provides an accurate characterisation of the colour. However, any particular set of R, G, B or X, Y, Z values could define any of a large number of different reflectance functions. The corresponding colours of these reflectance functions will produce the same colour under a reference light source, such as daylight. However, if an inappropriate reflectance function is derived from the camera R, G, B values, the colour of the object at the pixel in question may be defined in such a way that, for example, it appears to be a very different colour under a different light source, for example a tungsten light.
The apparatus according to a preferred embodiment of the invention allows:
• the colour of an object at a pixel or group of pixels to be measured in terms of tristimulus values;
• the colour of an object at a pixel or group of pixels to be measured in terms of reflectance values via spectral sensitivities of a camera (from the RGB equations above);
• the colour of an object at a pixel or group of pixels to be measured in terms of reflectance values via standard colour matching functions (from the X, Y, Z equations above).
According to the invention there is provided an apparatus for measuring colours of an object, the apparatus including: an enclosure for receiving the object; illumination means for illuminating the object within the enclosure; a digital camera for capturing an image of the object; a computer connected to the digital camera, for processing information relating to the image of the object; and display means for displaying information relating to the image of the object.
Where the term "digital camera" is used, it should be taken to be interchangeable with or to include other digital imaging means such as a colour scanner.
The enclosure may include means for mounting an object therein such that its position may be altered. These means may include a tiltable table for receiving the , object. Preferably the tiltable table is controllable by the computer.
Preferably the illumination means are located within the enclosure. The illumination means may include diffusing means for providing a diffuse light throughout the enclosure. Preferably the illumination means includes a plurality of different light sources for providing respectively different illuminations for the object. One or more of the light sources may be adjustable to adjust the level of the illumination or the direction of the illumination. The light sources may be controllable by the computer.
Preferably the digital camera is mounted on the enclosure and is directed into the enclosure for taking an image of the object within the enclosure. Preferably the camera is mounted such that its position relative to the enclosure may be varied. Preferably the location and/or the angle of the digital camera may be varied. The camera may be adjusted by the computer.
The display means may include a video display unit, which may include a cathode ray tube (CRT).
According to the invention there is further provided a method for measuring colours of an object, the method including the steps of: locating the object in an enclosure; illuminating the object within the enclosure; using a digital camera to capture an image of the object within the enclosure; using a computer to process information relating to the image of the object; and displaying selected information relating to the image of the object.
The method may include the step of illuminating the object with a number of respectively different light sources. The light may be diffuse. The light sources may be controlled by the computer.
The digital camera may also be controlled by the computer.
The method preferably includes the step of calibrating the digital camera, to transform its red, green, blue (R, G, B) signals into standard X, Y, Z values. The calibration step may include taking an image of a reference chart under one or more of the light sources and comparing the camera responses for each known colour within the reference chart with the standard X, Y, Z responses for that colour.
For each pixel, the relationship between the measured R, G, B values and the predicted X, Y, Z values is preferably represented as follows:
Figure imgf000008_0001
which can be expressed in the matrix form: X = MR , and hence M = XR '
The coefficients in the 3 by 11 matrix M are pre erably obtained via an optimisation method based on the least square technique, the measure used (Error) being as follows, where n=240 colours in a calibration chart:
Error = ± [ {XM ~XP Y + {YM ~YP 1 + {Zx ZP Y ] i=l
where X , Y , Z are the measured tristimulus values and X , Y , Z are the predicted tristimulus values.
The method may include the step of predicting a reflectance function for a pixel or group of pixels within the image of the object. The method may include the following steps: uniformly sampling the visible range of wavelengths (λ=a to λ = b) by choosing an integer n and specifying that
λ;
Figure imgf000008_0002
defining a relationship between camera output and reflectance function, using the following equation: P = Wr, where P includes known X , Y , Z values, W is a known weight matrix derived from the product of an illuminant function and the CIE x, y, z colour matching functions, W is the transposition of the matrix W and r is an unknown n component column vector representing reflectance function defined by:
*OD
where R(^.)\o R(^n) are the unknown reflectances of the observed object at each of the n different wavelengths; and finding a solution for P = Wr which includes a measure of both the smoothness and the colour constancy of the reflectance function, the relative importance of smoothness and of colour constancy being defined by respective weighting factors.
Using the above method, the camera is initially calibrated so that measured R, G, B values can be trans ormed to predicted X , Y , Z values. The X , Y , Z values may then be used to predict the reflectance functions.
Alternatively the R, G, B values may be used to predict reflectance functions directly using the following steps:
uniformly sampling the visible range of wavelengths (λ = to λ = b) by choosing an integer n and specifying that b -a
X. =a+(i - l)Aλ, i-1,2,. .n. with Δλ=- n-r defining a relationship between camera output and reflectance function, using the following equation: P = Wr, where P includes known camera R, G, B values, W is a known weight matrix derived from the product of an illuminant function and the CIE x, y, z colour matching functions, " is the transposition of the matrix W and r is an unknown n component column vector representing reflectance function defined by:
r = R _
R .)
where R ^^to E(^n) are the unknown reflectances of the observed object at each of the n different wavelengths; and finding a solution for P = Wr which includes a measure of both the smoothness and the colour constancy of the reflectance function, the relative importance of smoothness and of colour constancy being defined by respective weighting factors.
The weighting factors may be predetermined and are preferably calculated empirically.
Preferably n is at least 10. Most preferably n is at least 16, and n may be
31.
Preferably the smoothness is defined by determining the following:
Min Gr where G is an (n-l) x (n) matrix defined by the following:
Figure imgf000011_0001
where r is an unknown n component column vector representing reflectance function (referred to as the "reflectance vector") and |y|| is the 2- norm of the vector y, defined by
IMI = ∑. y~ ^ y 1S a vector with N components). l £=l
Preferably o ≤ r≤ e where o is an π component zero vector and e is an n component column vector where all the elements are unity (equal one).
Preferably the colour constancy of the reflectance vector is calculated as follows:- compute tristimulus X, Y, Z values (denoted PR) using the reflectance vector, under a reference illuminant; compute tristimulus X, Y, Z values (denoted P_) using the reflectance vector, under a test illuminant; using a chromatic adaptation transform, transfer P to a corresponding colour denoted by P TC under the reference illuminant; compute the difference ΔE between Pτc and PR; and define the colour inconstancy index (CON) as ΔE .
A plurality / of test illuminants may be used such that the colour j inconstancy index is defined as ^ β .AEj where β} is a weighting factor
defining the importance of colour constancy under a particular illuminant J. The reference illuminant is preferably D65, which represents daylight.
The preferred method for predicting the reflectance function may thus be defined as follows:- choose a reference illuminant and /test illuminants; choose a smoothness weighting factor α and weighting factors βf , j=l, 2,
•••J for CON; and for a given colour vector P and weight matrix W solve the following constrained non-linear problem:
Min aprf +] .βJAEJ '
T subject to o ≤ r ≤ e and P = J/fZ r for the reflectance vector
The smoothness weighting function α may be set to zero, such that the reflectance is generated with the least colour inconstancy.
The colour constancy weighting factors β}. may alternatively be set to zero, such that the reflectance vector has smoothness only.
Preferably α and β . are set such that the method generates a reflectance function having a high degree of smoothness and colour constancy. The values of α and β . may be determined by trial and error.
Preferably the method further includes the step of providing an indication of an appearance of texture within a selected area of the object. The method may include the steps of: determining an average colour value for the whole of the selected area; and determining a difference value at each pixel within the image of the selected area, the difference value representing the difference between the measured colour at that pixel and the average colour value for the selected area. Preferably the selected area has a substantially uniform colour.
The difference value may be a value A Y which represents the difference between the tristimulus value Y at that pixel and the average Y for the selected area.
Alternatively, the difference value may also include a value ΔX which represents the difference between the tristimulus value X at that pixel and the average X for the selected area and/or a value ΔZ which represents the difference between the tristimulus value Z at that pixel and the average Z for the selected area.
The texture of the selected area may be represented by an image comprising the difference values for all the respective pixels within the selected area.
The method may further include the step of simulating the texture of a selected area of an object, for example in an alternative, selected colour. The method may include the step of:
obtaining X, Y, Z values for the selected colour; converting these to x, y, Y values, where:
X Y x=- !«' = - Z =
X +Y + Z X + Y + Z X + Y +Z where x + y + z = I;
transforming the Y value for each pixel l,m to Y! m =Y + tΔY, m ,
where t is a function of Y. The x, ' y ' , ' and Y '.,« „', values for each p rixel may ' be converted to X l,m ,' Y l,m , Z l,,m values. The X, Y, Z values may then be transformed to monitor R, G, B values, for displaying the selected colour with the simulated texture on the display means.
Alternatively, the X, Y, Z values for each pixel l,m may be transformed to:
X l.m X + t A X l,,m
F l,m Y+ t A Y, z l.,m Z + t A Z l,,m
An embodiment of the invention will be described for the purpose of illustration only with reference to the accompanying drawings in which:
Fig. 1 is a diagrammatic overview of an apparatus according to the invention;
Fig. 2 is a diagrammatic sectional view of an illumination box for use with the apparatus of Fig. 1.
Referring to Fig. 1, an apparatus according to the invention includes an illumination box 10 in which an object 18 to be observed may be placed. A digital camera 12 is located towards the top of the illumination box 10 so that the digital camera 12 may take a picture of the object 18 enclosed in the illumination box 10. The digital camera 12 is connected to a computer 14 provided with a video display unit (VDU) 16, which includes a colour sensor 30.
Referring to Fig. 2, the illumination box 10 is provided with light sources 20 which are able to provide a very carefully controlled illumination within the box 10. Each light source includes a lamp 21 and a diffuser 22, through which the light passes in order to provide uniform, diffuse light within the illumination box 10. The inner surfaces of the illumination box are of a highly diffusive material coated with a matt paint for ensuring that the light within the box is diffused and uniform.
The light sources are able to provide a variety of different illuminations within the illumination box 10, including: D65, which represents daylight; tungsten light; and lights equivalent to those used in various department stores, etc. In each case the illumination is fully characterised, i.e., the amounts of the various different wavelengths of light are known.
The illumination box 10 includes a tiltable table 24 on which the object 18 may be placed. This allows the angle of the object to be adjusted, allowing different parts of the object to be viewed by the camera.
The camera 12 is mounted on a slider 26, which allows the camera to move up and down as viewed in Fig. 2. This allows the lens of the camera to be brought closer to and further away from the object, as desired. The orientation of the camera may also be adjusted.
Referring again to Fig. 1, the light sources 20, the digital camera 12 and its slider 26 and the tiltable table 24 may all be controllable automatically from the computer 14. Alternatively, control may be effected from control buttons on the illumination box or directly by manual manipulation.
The digital camera 12 is connected to the computer 14 which is in turn connected to the VDU 16. The image taken by the camera 12 is processed by the computer 14 and all or selected parts of that image or colours or textures within that image may be displayed on the VDU and analysed in various ways. This is described in more detail hereinafter.
The digital camera describes the colour of the object at each pixel in terms of red (R), green (G) and blue (B) signals, which are expressed in the following equations
Figure imgf000016_0001
G = k' )s{λ)g(λ)R(λ)dλ Equation 1
Figure imgf000016_0002
a
S(λ) is the spectral power distribution of the illuminant. Given that the object is illuminated within the illumination box 10 by the light sources 20, the spectral power distribution of any illuminant used is known. R(λ) is the reflectance function of the object at the pixel in question (which is unknown) and r,g,b are the spectral sensitivities of the digital camera, i.e., the responses of the charge coupled device (CCD) sensors used by the camera.
All the above functions are defined within the visible range, typically between a=400 and b=700 nm.
There are known calibration methods for converting a digital camera's R, G, B signals in the above equation into the CIE tristimulus values (X, Y, Z). The tristimulus values are defined in the following equations:
Figure imgf000016_0003
b _
Y = k S(λ)y{λ)R(λ)dλ Equation 2 b
Figure imgf000016_0004
where all the other functions in equation (1) were defined. The x,y,z are the
CIE 1931 or 1964 standard colorimetric observer functions, also known as colour matching functions (CMF), which define the amounts of reference red, green and blue lights in order to match a monochromatic light in the visible range. The k factor in equation (2) is a normalising factor to make Y equal to 100 for a reference white. In order that the R, G, B values captured by the digital camera may be transformed into X, Y, Z values, it is desirable to calibrate the digital camera before the apparatus is used to measure colours of the object 18. This is done each time the camera is switched on or whenever the light source or camera setting is altered. Preferably the camera is calibrated by using a standard colour chart, such as a GretagMacbeth ColorChecker Chart or Digital Chart.
The chart is placed in the illumination box 10 and the camera 12 takes an image of the chart. For each colour in the chart, the X, Y, Z values are known. The values are obtained either from the suppliers of the chart or by measuring the colours in the chart by using a colour measuring instrument. A polynomial modelling technique may be used to transform from the camera R, G, B values to X, Y, Z values. For a captured image from the camera, each pixel represented by R, G, B values is transformed using the following equation to predict X , Y , Z values, these being the X, Y, Z values at a particular pixel:
Figure imgf000017_0001
which can be expressed in the matrix form: X =MR , and hence, M = XR~
The coefficients in the 3 by 11 matrix M may be obtained via an optimisation method based on a least squares technique. The measure used (Error) is as follows, where n=240 colours in a standard calibration chart:
Figure imgf000018_0001
[ ( - + (Y -YP + (zM-zP f ]
where X 7 > ZM are the measured tristimulus values and X „_ YA Z^ are the predicted tristimulus values.
Using the above technique, the digital camera may be calibrated such that its R, G, B readings for any particular colour may be accurately transformed into standard X, Y, Z values.
It is also necessary to characterise the VDU 16. This may be carried out using known techniques, such as are described in Berns R.S. et al, CRT colorimetry, Part I and II at Col, Res Appn, 1993.
Once the camera 12 and VDU 16 have been calibrated, a sample object may be placed into the illumination box 10. The digital camera is controlled directly or via the computer 14, to take an image of the object 18. The image may be displayed on the VDU 16. In analysing and displaying the image, the apparatus preferably predicts the' reflectance function of the object at each pixel. . This ensures that the colour of the object is realistically characterised and can be displayed accurably on the VDU, and reproduced on other objects if required.
One method of predicting reflectance functions from R, G, B or X, Y, Z values is as follows.
If we uniformly sample the visible range (a, b) by choosing an integer n and
= a + (i-l)Δλ, i = l,2,- -,n with Δλ = ^X^- / ' n - 1 then the equations (1) and (2) can be rewritten as the following matrix vector form to define a relationship between camera output and reflectance function: P = WΥ Equation 3
Here, p is a 3 -component column vector consisting of the camera response, W is a n x 3 matrix called the weight matrix, derived from the illuminant function and the sensors of the camera for equation (1), or from the illuminant used, and the colour matching functions for equation (2), W is the transposition of the matrix W and r is the unknown n component column vector (the reflectance vector) representing unknown reflectance function given by:
R(λ r = R(λ_
R(λ.) Equation 4
The 3-component column vector p consists of either the camera responses R, G and B for the equation (1), or the CIE tristimulus values X, Y and Z for the equation (2).
Note also that the reflectance function R(λ) should satisfy:
0 ≤ R{λ) ≤ l Thus, the reflectance vector r defined by equation (4) should satisfy:
° - r - e Equation 5
Here o is a π-component zero vector and e is a rz-component vector where all the elements are unity (equal one).
Some fluorescent materials have reflectances of more than 1, but this method is not generally applicable to characterising the colours of such materials.
The preferred method used with the present invention recovers the reflectance vector r satisfying equation (3) by knowing all the other parameters or functions in equations (1) and (2).
The method uses a numerical approach and generates a reflectance vector r defined by equation (4) that is smooth and has a high degree of colour constancy. In the surface industries, it is highly desirable to produce colour constant products, i.e., the colour appearance of the goods will not be changed when viewed under a wide range of light sources such as daylight, store lighting, tungsten.
Firstly, a smoothness constraint condition is defined as follows:
Min « ,,2 \\Or\\ r Here G is an (n-1) x n matrix referred to as the "smooth operator", and defined by the following:
Figure imgf000020_0001
where r is the unknown reflectance vector defined by equation (4) and ||>>|| is the 2-norm of the vector y, defined by:
Figure imgf000020_0002
if y is a vector with N components. Since the vector r should satisfy equations (3) and (5), therefore, the smoothness vector r is the solution of the following constrained least squares problem: Min
\\Gι \\2 o ≤ r ≤ e
T subject to p = W r r is always between 0 and 1, i.e., within the defined boundary.
It is assumed that the reflectance vector r generated by the above smoothness approach has a high degree of colour constancy. However, it has been realised by the inventors that the colour constancy of such reflectance vector may be improved as follows.
A procedure for calculating a colour inconstancy index CON of the reflectance vector r is described below.
1. Compute tristimulus values denoted by PR. using the reflectance vector under a reference illuminant.
2. Compute tristimulus values denoted by P_, using the reflectance vector under a test illuminant.
3. Using a reliable chromatic adaptation transform such as CMCCAT97, transfer P_. to a corresponding colour denoted by Pτc under the reference illuminant.
4. Using a reliable colour difference formula such as CIEDE2000, compute the difference ΔE between P R_ and P T_C„ under the reference illuminant.
5. Define CON as ΔE .
The chromatic transform CMCCAT97 is described in the following paper: "M R Luo and R W G Hunt, A chromatic adaptation transform and a colour inconstancy Index, Color Res Appn, 1998". The colour difference formula is described in "M R Luo, G Cui and B Rigg, The development of the CIE 2000 colour difference Formula: CIEDE2000, Color Res Appn, 2001". The reference and test illuminants are provided by the illumination box 10 and are thus fully characterised, allowing the above calculations to be carried out accurately. The method may be summarised as follows:
Choose the reference illuminant (say D65) and / test illuminants (A, Fll, etc).
Choose the smoothness weighting factor α and the weighting factors β , j = 1,2, -,J for CON.
For a given colour vector p and using a known weight matrix W in equation (3), solve the following constrained non-linear problem:
— [4Grf + Σ β jAEj ] Equation 6
7 = 1
T subject to o ≤ r ≤ e, and p = Jf/~ r for the reflectance vector r.
If the smoothness weighting factor is set to 0, then the above method generates the reflectance with the least colour inconstancy. However, the reflectance vector r could be too fluctuated to be realistic. At the other extreme, if the weighting factors β s are all set to be zero, then the above method produces a reflectance vector r with smoothness only. By choosing appropriate weighting factors, α and β , the above method generates reflectances with smoothness and a high degree of colour constancy.
The weight matrix W should be known from the camera characterisation carried out before the apparatus is used to measure the colours of the object 18.
The above described method for predicting a reflectance function from the digital camera's red, green and blue signals results in a reflectance function which is smooth and colour constant across a number of illuminants. Using the above method, the apparatus is able to characterise and reproduce a colour of the object 18 very realistically and in such a way that the colour is relatively uniform in appearance under various different illuminants.
In industrial design, it is frequently also desired to simulate products in different colours. For example, a fabric of a particular texture might be available in green and the designer may wish to view an equivalent fabric in red. The apparatus according to the invention allows this to be done as follows.
An image of the existing object 18 is taken using the digital camera 12 and a particular area of uniform colour to be analysed is isolated from the background using known software.
Within the above selected area of colour, the R, G, B values are transformed to standardised X, Y, Z values.
Average colour values X, Y, Z are calculated, these being the mean X, Y, Z values for the whole selected area of colour.
At each pixel, a difference value ΔY is calculated, ΔY being equal to the difference between the Y value at the pixel in question and the average Y value,
7 , such that: ΔYl m equals Ylm -Y, where l,m represents a particular pixel.
The computer calculates ΔY values at each pixel within the selected area of colour in the image. Because the colour of the area is uniform, the variations in the measured Y values from the average Y value must represent textural effects. Thus the computer can create a "texture profile" for the area of colour, the profile being substantially independent of the colour of the area.
According to the above method only ΔY values (and not -INand ΔZ values) are used. The applicants have found that the perceived lightness of an area within an image has much to do with the green response and that the ΔY values give a very good indication of lightness and therefore of texture.
Once ΔY values are stored for each pixel in the selected area, providing the texture profile, this may be used to simulate a similar texture in a different colour. This is carried out as follows.
Firstly the new colour is measured or theoretical colour values provided. The X, Y, Z values are transformed to x, y, Y, where
X 7 Z y=-
X + Y + Z X + Y + Z X + Y + Z and x + y + z = 1
The X, Y, Z colour space is not very uniform, including very small areas of blue and very large areas of green. The above transform transfers the colour to x, y, Y space in which the various colours are more uniformly represented.
To retain the chosen colour but to superimpose the texture profile of the previously characterised colour, the x, y values remain the same and the 7 value is replaced with a value Yl m , for a pixel l,m.
Y U, rn Y +tΔY l,,m
Alternatively, the X, Y, Z values for each pixel l,m may be transformed to:
_= x,m X + t A X,
= Y + t yi« y Δ Y l,,m
=
Z,m Z + t z A Z l.m
This takes into account the lightness of the red and blue response as well as the green response. Thus, the lightness values and thus the texture profile of the previous material have been transferred to the new colour.
The term t varies with Y but there are different functions of t against Y for different materials, with the relationship between f and Y depending upon the coarseness of the material. The appropriate values of t may be calculated empirically.
There is thus provided an apparatus and method for providing accurate and versatile information about colours of objects, for capturing high colour fidelity and repeatable images and for simulating different colours of a product having the same texture. The illumination box 10 allows objects to be viewed in controlled conditions under a variety of accurately characterised lights. This, preferably together with the novel method for predicting reflectance functions, enables colours to be characterised in such a way that they are predictable and realistically characterised under all lights. The apparatus and method also provide additional functions such as the ability to superimpose a texture of one fabric on to a different coloured fabric.
Whilst endeavouring in the foregoing specification to draw attention to those features of the invention believed to be of particular importance it should be understood that the Applicants claim protection in respect of any patentable feature or combination of features hereinbefore referred to and/or shown in the drawings whether or not particular emphasis has been placed thereon.

Claims

1. Apparatus for measuring colours of an object, the apparatus including: an enclosure for receiving the object; illumination means for illuminating the object within the enclosure; a digital camera for capturing an image of the object; a computer connected to the digital camera, for processing information relating to the image of the object; and display means for displaying information relating to the image of the object.
2. Apparatus according to claim 1 wherein the enclosure includes means for mounting an object therein such that its position may be altered.
3. Apparatus according to claim 2 wherein the mounting means includes a tiltable table for receiving the object, the tiltable table being controllable by the computer.
4. Apparatus according to any preceding claim wherein the illumination means are located within the enclosure, and include diffusing means for providing a diffuse light throughout the enclosure.
5. Apparatus according to any preceding claim wherein the illumination means includes a plurality of different light sources for providing respectively different illuminations for the object, one or more of the light sources being adjustable to adjust the level of the illumination or the direction of the illumination, and the light sources being controllable by the computer.
6. Apparatus according to any preceding claim wherein the digital camera is mounted on the enclosure and is directed into the enclosure for taking an image of the object within the enclosure.
7. Apparatus according to claim 6 wherein the camera is mounted such that its position relative to the enclosure may be varied, and the location and/or the angle of the digital camera may be varied.
8. Apparatus according to claim 7 wherein the camera may be adjusted by the computer.
9. Apparatus according to any preceding claim wherein the display means includes a video display unit including a cathode ray tube (CRT).
10. A method for measuring colours of an object, the method including the steps of: locating the object in an enclosure; illuminating the object within the enclosure; using a digital camera to capture an image of the object within the enclosure; using a computer to process information relating to the image of the object; and displaying selected information relating to the image of the object.
11. A method according to claim 10 wherein the step of illuminating the object with a number of respectively different light sources.
12. A method according to claim 10 or claim 11, the method including the step of calibrating the digital camera, to transform its red, green, blue (R, G, B) signals into standard X, Y, Z values, the calibration step includes taking an image of a reference chart under one or more of the light sources and comparing the camera responses for each known colour within the reference chart with the standard X, Y, Z responses for that colour.
13. A method according to claim 11 or claim 12, the method including the following steps: uniformly sampling the visible range of wavelengths (λ = a to λ = b) by choosing an integer n and specifying that a λi =a+(i - l)Δλ, i = l,2,. with Δλ=- n - 1 defining a relationship between camera output and reflectance function, using the following equation: P = Wr, where P includes known X P , Y P , ' Z P values, ' W is a known weig °ht matrix derived from the product of an illuminant function and the CIE x, y, z colour matching functions, W is the transposition of the matrix W and r is an unknown n component column vector representing reflectance function defined by:
RW r = R(λ_)
(λ,)
where Rf^^Xo E(^() are the unknown reflectances of the observed object at each of the n different wavelengths; and finding a solution for P = Wr which includes a measure of both the smoothness and the colour constancy of the reflectance function, the relative importance of smoothness and of colour constancy being defined by respective weighting factors.
14. A method according to claim 11 or claim 12, the method including the following steps: uniformly sampling the visible range of wavelengths (λ=a to λ = b) by choosing an integer n and specifying that b — a λi -a+(i - l)Δλ, i=l,2, ...ti, with Δλ=- n- 1 defining a relationship between camera output and reflectance function, using the following equation: P = Wr, where P includes known camera R, G, B values, W is a known weight matrix derived from the product of an illuminant function and the CIE x, y, z colour matching functions, Wτ is the transposition of the matrix W and r is an unknown n component column vector representing reflectance function defined by:
R(λ r = _)
*ou
where
Figure imgf000029_0001
are the unknown reflectances of the observed object at each of the n different wavelengths; and finding a solution for P = Wr which includes a measure of both the smoothness and the colour constancy of the reflectance function, the relative importance of smoothness and of colour constancy being defined by respective weighting factors.
15. A method according to any of claims 12 to 14 wherein the weighting factors are predetermined, being calculated empirically.
16. A method according to any of claims 11 or claim 15 wherein n is at least 16.
17. A method according to claims 11 to 16, wherein the smoothness is defined by determining the following:
Min Gr where G is an (n-ϊ) x (ri) matrix defined by the following:
Figure imgf000030_0001
where r is an unknown n component column vector representing reflectance function (referred to as the "reflectance vector") and y| is the 2- norm of the vector y, defined by
_ _ y (if y is a vector with N components).
18. A method according to claims 17, wherein o ≤ r≤ e where o is an n component zero vector and e is an n component column vector where all the elements are unity (equal 1).
19. A method according to any of claims 11 to 18, wherein the colour constancy of the reflectance vector is calculated as follows:- compute tristimulus X, Y, Z values (denoted P ) using the reflectance vector, under a reference illuminant; compute tristimulus X, Y, Z values (denoted P_) using the reflectance vector, under a test illuminant; using a chromatic adaptation transform, transfer P,_ to a corresponding colour denoted by Pτc under the reference illuminant; compute the difference ΔE between Pτc and PR, and define the colour inconstancy index (CON) as ΔE .
20. A method according to claim 19 wherein a plurality /of test illuminants is
J used such that the colour inconstancy index is defined as β ΔEj where β }. is j=l a weighting factor defining the importance of colour constancy under a particular illuminant j.
21. A method according to any of claims 11 to 20 wherein the method further includes the step of providing an indication of an appearance of texture within a selected area of the object, the method including the steps of: determining an average colour value for the whole of the selected area; and determining a difference value at each pixel within the image of the selected area, the difference value representing the difference between the measured colour at that pixel and the average colour value for the selected area.
22. A method according to claim 21 wherein the selected area has a substantially uniform colour.
23. A method according to claim 21 or claim 22 wherein difference value is a value Δ Y which represents the difference between the tristimulus value 7 at that pixel and the average 7 for the selected area.
24. A method according to claim 23 wherein the difference value also includes a value AN which represents the difference between the tristimulus value X at that pixel and the average X for the selected area and/or a value ΔZ which represents the difference between the tristimulus value Z at that pixel and the average Z for the selected area.
25. A method according to claims 21 to claim 24 wherein texture of the selected area may be represented by an image comprising the difference values for all the respective pixels within the selected area.
26. A method according to any of claims 11 to 25, the method further including the step of simulating the texture of a selected area of an object, for example in an alternative, selected colour, by:
obtaining X, Y, Z values or the selected colour; converting these to x, y, 7 values, where:
X 7 x=- y- z =
X + Y + Z X + Y + Z X + Y + Z where x + y + z = 1;
transforming the 7value for each pixel l,m to Yl m =Y + t.ΔY ,
where t is a function of 7.
27. A method according to claim 26 wherein the x, y, and Yl m values for each p rixel are converted to X l,m , Y l,m , Z l,,ιn values and the X, ' Y, ' Z values are then transformed to monitor R, G, B values, for displaying the selected colour with the simulated texture on the display means.
28. A method according to claim 26, wherein the X, Y, Z values for each pixel l,m are be transformed to:
x, = X + t x A X l,,m
=
Y,,m Y+ t y A Y ι,,m
=
,m Z + t z A Z \,,m
29. Apparatus substantially as herein described with reference to the drawings.
30. A method substantially as herein described with reference to the drawings.
31. Any novel subject matter or combination including novel subject matter disclosed herein, whether or not within the scope of or relating to the same invention as any of the preceding claims.
PCT/GB2002/004521 2001-10-04 2002-10-04 Apparatus and method for measuring colour WO2003029766A2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US10/491,706 US20050018191A1 (en) 2001-10-04 2002-10-04 Apparatus and method for measuring colour
AU2002329470A AU2002329470A1 (en) 2001-10-04 2002-10-04 Apparatus and method for measuring colour
EP02765098A EP1436577A2 (en) 2001-10-04 2002-10-04 Apparatus and method for measuring colour

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
GB0123810A GB0123810D0 (en) 2001-10-04 2001-10-04 Method of predicting reflectance functions
GB0123810.4 2001-10-04
GB0124683.4 2001-10-15
GB0124683A GB0124683D0 (en) 2001-10-04 2001-10-15 Apparatus and method for measuring colour

Publications (2)

Publication Number Publication Date
WO2003029766A2 true WO2003029766A2 (en) 2003-04-10
WO2003029766A3 WO2003029766A3 (en) 2003-07-24

Family

ID=26246608

Family Applications (3)

Application Number Title Priority Date Filing Date
PCT/GB2002/004528 WO2003029811A1 (en) 2001-10-04 2002-10-04 Assessing colour fastiness
PCT/GB2002/004521 WO2003029766A2 (en) 2001-10-04 2002-10-04 Apparatus and method for measuring colour
PCT/GB2002/004500 WO2003030524A2 (en) 2001-10-04 2002-10-04 A method of predicting reflectance functions

Family Applications Before (1)

Application Number Title Priority Date Filing Date
PCT/GB2002/004528 WO2003029811A1 (en) 2001-10-04 2002-10-04 Assessing colour fastiness

Family Applications After (1)

Application Number Title Priority Date Filing Date
PCT/GB2002/004500 WO2003030524A2 (en) 2001-10-04 2002-10-04 A method of predicting reflectance functions

Country Status (3)

Country Link
US (1) US20050018191A1 (en)
EP (1) EP1436577A2 (en)
WO (3) WO2003029811A1 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006020833A2 (en) 2004-08-11 2006-02-23 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
WO2006064010A1 (en) 2004-12-14 2006-06-22 Akzo Nobel Coatings International B.V. Method and device for analysing visual properties of a surface
WO2006064011A1 (en) 2004-12-14 2006-06-22 Akzo Nobel Coatings International B.V. Method and device for measuring coarseness of a paint film
WO2006134429A2 (en) * 2005-02-08 2006-12-21 Alessandro Occelli Device for analysing the colour of a unhomogeneous material, like hair, and method thereof
US7599559B2 (en) 2004-05-13 2009-10-06 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US7751653B2 (en) 2004-05-13 2010-07-06 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US7974466B2 (en) 2004-11-23 2011-07-05 Color Savvy Systems Limited Method for deriving consistent, repeatable color measurements from data provided by a digital imaging device
US8044969B2 (en) 2005-03-04 2011-10-25 Chrometrics Limited Reflectance spectra estimation and colour space conversion using reference reflectance spectra
EP2526684B1 (en) 2010-01-19 2019-04-10 Akzo Nobel Coatings International B.V. Method and system for determining colour from an image

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2331196T3 (en) * 2004-09-17 2009-12-23 Akzo Nobel Coatings International Bv METHOD TO EQUAL A PAINT.
US8345252B2 (en) * 2005-04-25 2013-01-01 X-Rite, Inc. Method and system for enhanced formulation and visualization rendering
CN101346082B (en) * 2005-12-23 2012-03-21 皇家飞利浦电子股份有限公司 Color matching for display system for shops
FR2908427B1 (en) * 2006-11-15 2009-12-25 Skin Up PROCESS FOR IMPREGNATING FIBERS AND / OR TEXTILES WITH A COMPOUND OF INTEREST AND / OR AN ACTIVE INGREDIENT IN THE FORM OF NANOPARTICLES
DE102008013387B4 (en) * 2008-03-10 2020-02-13 Byk-Gardner Gmbh Device for determining the optical surface properties of workpieces
CN102236008B (en) * 2011-02-22 2014-03-12 晋江市龙兴隆染织实业有限公司 Method for detecting color fastness of fabric products to water
CN102359819B (en) * 2011-09-21 2013-10-02 温州佳易仪器有限公司 Color detection method of multi-light-source colorful image and color collection box used by color detection method
ES2955491T3 (en) 2012-11-02 2023-12-01 Variable Inc Computer-implemented system and method for color detection, storage and comparison
CN103925992B (en) * 2013-01-16 2016-03-16 光宝电子(广州)有限公司 There are the brightness measurement method and system of the device of backlight
CN103063310A (en) * 2013-01-18 2013-04-24 岑夏凤 Non-contact type color measurement method and non-contact type color measurement device based on digital technology
DE102014201124A1 (en) * 2014-01-22 2015-07-23 Zumtobel Lighting Gmbh Method for controlling a lighting arrangement
US20160034944A1 (en) * 2014-08-04 2016-02-04 Oren Raab Integrated mobile listing service
EP3218682B1 (en) * 2014-11-13 2022-08-10 BASF Coatings GmbH Reference number for determining a colour quality
US10156477B2 (en) 2015-05-01 2018-12-18 Variable, Inc. Intelligent alignment system and method for color sensing devices
CN104897374B (en) * 2015-06-16 2016-04-13 常州千明智能照明科技有限公司 A kind of color calibration method of camera module
CN105445182B (en) * 2015-11-17 2017-12-05 陕西科技大学 Color fastness detection sampler on the inside of a kind of footwear
CN105445271B (en) * 2015-12-02 2018-10-19 陕西科技大学 A kind of device and its detection method of real-time detection colour fastness to rubbing
US11002676B2 (en) 2018-04-09 2021-05-11 Hunter Associates Laboratory, Inc. UV-VIS spectroscopy instrument and methods for color appearance and difference measurement
US10746599B2 (en) 2018-10-30 2020-08-18 Variable, Inc. System and method for spectral interpolation using multiple illumination sources
CN109632647A (en) * 2018-11-29 2019-04-16 上海烟草集团有限责任公司 The binding strength detection method of printed matter, system, storage medium, electronic equipment
CN110286048B (en) * 2019-06-13 2021-09-17 杭州中服科创研究院有限公司 Textile fabric color fastness detection equipment
CN112362578A (en) * 2020-11-10 2021-02-12 云南中烟工业有限责任公司 Method for measuring cigarette tipping paper lip adhesion according to color fastness
EP4187217A1 (en) * 2021-11-26 2023-05-31 Kuraray Europe GmbH Mobile computing device for performing color fastness measurements

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4648051A (en) * 1984-10-15 1987-03-03 The Board Of Trustees Of The Leland Stanford Junior University Color imaging process
US4812904A (en) * 1986-08-11 1989-03-14 Megatronics, Incorporated Optical color analysis process
EP0570003A2 (en) * 1992-05-15 1993-11-18 Toyota Jidosha Kabushiki Kaisha Rendering apparatus, multispectral image scanner, and three-dimensional automatic gonio-spectrophotometer
US5526285A (en) * 1993-10-04 1996-06-11 General Electric Company Imaging color sensor
EP0780671A1 (en) * 1995-12-21 1997-06-25 Shimadzu Corporation Spectrophotometers and colorimeters
US5844680A (en) * 1994-09-24 1998-12-01 Byk-Gardner Gmbh Device and process for measuring and analysing spectral radiation, in particular for measuring and analysing color characteristics
US5850472A (en) * 1995-09-22 1998-12-15 Color And Appearance Technology, Inc. Colorimetric imaging system for measuring color and appearance

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5041328A (en) * 1986-12-29 1991-08-20 Canon Kabushiki Kaisha Recording medium and ink jet recording method by use thereof
JPH02258345A (en) * 1989-03-31 1990-10-19 Toppan Printing Co Ltd Decoloration tester for printed matter
JPH04199969A (en) * 1990-11-29 1992-07-21 Canon Inc Image reader
JPH05119672A (en) * 1991-10-25 1993-05-18 Mita Ind Co Ltd Decolorizing machine
JP3577503B2 (en) * 1992-09-28 2004-10-13 大日本インキ化学工業株式会社 Color code
JP3310786B2 (en) * 1994-08-22 2002-08-05 富士写真フイルム株式会社 Color thermal recording paper package and color thermal printer
US5633722A (en) * 1995-06-08 1997-05-27 Wasinger; Eric M. System for color and shade monitoring of fabrics or garments during processing
US5740078A (en) * 1995-12-18 1998-04-14 General Electric Company Method and system for determining optimum colorant loading using merit functions
JPH09327945A (en) * 1996-04-11 1997-12-22 Fuji Photo Film Co Ltd Recording material and image recording method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4648051A (en) * 1984-10-15 1987-03-03 The Board Of Trustees Of The Leland Stanford Junior University Color imaging process
US4812904A (en) * 1986-08-11 1989-03-14 Megatronics, Incorporated Optical color analysis process
EP0570003A2 (en) * 1992-05-15 1993-11-18 Toyota Jidosha Kabushiki Kaisha Rendering apparatus, multispectral image scanner, and three-dimensional automatic gonio-spectrophotometer
US5526285A (en) * 1993-10-04 1996-06-11 General Electric Company Imaging color sensor
US5844680A (en) * 1994-09-24 1998-12-01 Byk-Gardner Gmbh Device and process for measuring and analysing spectral radiation, in particular for measuring and analysing color characteristics
US5850472A (en) * 1995-09-22 1998-12-15 Color And Appearance Technology, Inc. Colorimetric imaging system for measuring color and appearance
EP0780671A1 (en) * 1995-12-21 1997-06-25 Shimadzu Corporation Spectrophotometers and colorimeters

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LUO M R ET AL: "A CHROMATIC ADAPTATION TRANSFORM AND A COLOUR INCONSTANCY INDEX" COLOR RESEARCH AND APPLICATIONS, JOHN WILEY AND SONS, NEW YORK, NY, US, June 1998 (1998-06), pages 154-158, XP000984144 ISSN: 0361-2317 *
See also references of EP1436577A2 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7599559B2 (en) 2004-05-13 2009-10-06 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US8320663B2 (en) 2004-05-13 2012-11-27 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US7907780B2 (en) 2004-05-13 2011-03-15 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US7751653B2 (en) 2004-05-13 2010-07-06 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
WO2006020833A3 (en) * 2004-08-11 2006-04-13 Color Savvy Systems Ltd Method for collecting data for color measurements from a digital electronic image capturing device or system
WO2006020833A2 (en) 2004-08-11 2006-02-23 Color Savvy Systems Limited Method for collecting data for color measurements from a digital electronic image capturing device or system
US7974466B2 (en) 2004-11-23 2011-07-05 Color Savvy Systems Limited Method for deriving consistent, repeatable color measurements from data provided by a digital imaging device
WO2006064011A1 (en) 2004-12-14 2006-06-22 Akzo Nobel Coatings International B.V. Method and device for measuring coarseness of a paint film
WO2006064010A1 (en) 2004-12-14 2006-06-22 Akzo Nobel Coatings International B.V. Method and device for analysing visual properties of a surface
WO2006134429A3 (en) * 2005-02-08 2007-06-14 Alessandro Occelli Device for analysing the colour of a unhomogeneous material, like hair, and method thereof
WO2006134429A2 (en) * 2005-02-08 2006-12-21 Alessandro Occelli Device for analysing the colour of a unhomogeneous material, like hair, and method thereof
US8044969B2 (en) 2005-03-04 2011-10-25 Chrometrics Limited Reflectance spectra estimation and colour space conversion using reference reflectance spectra
EP2526684B1 (en) 2010-01-19 2019-04-10 Akzo Nobel Coatings International B.V. Method and system for determining colour from an image

Also Published As

Publication number Publication date
WO2003030524A2 (en) 2003-04-10
WO2003030524A3 (en) 2003-05-15
US20050018191A1 (en) 2005-01-27
WO2003029766A3 (en) 2003-07-24
WO2003029811A1 (en) 2003-04-10
EP1436577A2 (en) 2004-07-14

Similar Documents

Publication Publication Date Title
WO2003029766A2 (en) Apparatus and method for measuring colour
CN109417586B (en) Color conversion system, color conversion device, and color conversion method
US4884130A (en) Method of describing a color in a triaxial planar vector color space
US5680327A (en) Apparatus and process for a digital swatchbook
JP5967441B2 (en) Color processing method, color processing apparatus, and color processing system
US5694227A (en) Method and apparatus for calibrating and adjusting a color imaging system
Luo Applying colour science in colour design
KR100238960B1 (en) Apparatus for chromatic vision measurement
EP1218706A1 (en) Methods for colour matching by means of an electronic imaging device
US7187797B2 (en) Color machine vision system for colorimetry
WO2005104659A2 (en) Method and system for approximating the spectrum of a plurality of color samples
JP2000184223A (en) Calibration method for scanner, image forming device and calibration processor
US6614530B1 (en) Method and device for the colorimetric measurement of a colored surface
MacDonald et al. Colour characterisation of a high-resolution digital camera
Zhu et al. Color calibration for colorized vision system with digital sensor and LED array illuminator
Connolly et al. Colour measurement by video camera
JP4088016B2 (en) Color management system
Hirschler Electronic colour communication in the textile and apparel industry
Pawlik et al. Color Formation in Virtual Reality 3D 360 Cameras
Connolly et al. Industrial colour inspection by video camera
JPH11132849A (en) Method for correcting colorimetry value and colorimetery machine
van Aken Portable spectrophotometer for electronic prepress
Alstergren Development of a spectral and goniophotometric imagingmeasurement system for optical characterization of materials
Rich Critical parameters in the measurement of the color of nonimpact printing
Pointer et al. Perceived colour differences in displayed colours Part 1: hard copy to soft copy matching

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BY BZ CA CH CN CO CR CU CZ DE DM DZ EC EE ES FI GB GD GE GH HR HU ID IL IN IS JP KE KG KP KR LC LK LR LS LT LU LV MA MD MG MN MW MX MZ NO NZ OM PH PL PT RU SD SE SG SI SK SL TJ TM TN TR TZ UA UG US UZ VC VN YU ZA ZM

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GH GM KE LS MW MZ SD SL SZ UG ZM ZW AM AZ BY KG KZ RU TJ TM AT BE BG CH CY CZ DK EE ES FI FR GB GR IE IT LU MC PT SE SK TR BF BJ CF CG CI GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2002765098

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2002765098

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 10491706

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP

WWW Wipo information: withdrawn in national office

Ref document number: 2002765098

Country of ref document: EP