WO2012074014A1 - 画像色推定方法、画像色推定装置、及び画像色推定プログラム - Google Patents

画像色推定方法、画像色推定装置、及び画像色推定プログラム Download PDF

Info

Publication number
WO2012074014A1
WO2012074014A1 PCT/JP2011/077687 JP2011077687W WO2012074014A1 WO 2012074014 A1 WO2012074014 A1 WO 2012074014A1 JP 2011077687 W JP2011077687 W JP 2011077687W WO 2012074014 A1 WO2012074014 A1 WO 2012074014A1
Authority
WO
WIPO (PCT)
Prior art keywords
component
value
signal value
display
gradation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2011/077687
Other languages
English (en)
French (fr)
Inventor
美文 下平
大野　勝也
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shizuoka University NUC
Original Assignee
Shizuoka University NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shizuoka University NUC filed Critical Shizuoka University NUC
Priority to US13/989,978 priority Critical patent/US9311885B2/en
Priority to JP2012546917A priority patent/JP5825600B2/ja
Publication of WO2012074014A1 publication Critical patent/WO2012074014A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
    • G09G5/00—Control arrangements or circuits for visual indicators common to cathode-ray tube indicators and other visual indicators
    • G09G5/02—Control arrangements or circuits for visual indicators common to cathode-ray tube indicators and other visual indicators characterised by the way in which colour is displayed
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46—Colour picture communication systems
    • H04N1/56—Processing of colour picture signals
    • H04N1/60—Colour correction or control
    • H04N1/6027—Correction or control of colour gradation or colour contrast
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/46—Colour picture communication systems
    • H04N1/56—Processing of colour picture signals
    • H04N1/60—Colour correction or control
    • H04N1/603—Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00—Details of colour television systems
    • H04N9/64—Circuits for processing colour signals
    • H04N9/67—Circuits for processing colour signals for matrixing
    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
    • G09G2320/00—Control of display operating conditions
    • G09G2320/06—Adjustment of display parameters
    • G09G2320/0693—Calibration of display systems
    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
    • G09G2360/00—Aspects of the architecture of display systems
    • G09G2360/14—Detecting light within display terminals, e.g. using a single or a plurality of photosensors
    • G09G2360/145—Detecting light within display terminals, e.g. using a single or a plurality of photosensors the light originating from the display screen

Definitions

  • the present invention relates to an image color estimation method for estimating an image color displayed on an image display device such as a display, an image color estimation device for estimating an image color, and an image color estimation program for estimating an image color.
  • Patent Document 1 describes a technique related to a display characteristic calibration method for calibrating display characteristics of a color display device.
  • Patent Document 2 describes a technique related to a measurement method for measuring display characteristics of a display.
  • Patent Document 3 describes a technique related to a system for creating a monitor profile. According to this technique, a monitor profile for performing high-accuracy correction according to the display characteristics of each monitor can be easily created.
  • CMS color management system
  • an image color estimation method capable of accurately estimating the image color of an image displayed on a display in which additive color mixture is not established and color tracking occurs, an image color estimation device for estimating an image color, and an image color
  • An image color estimation program for estimating the image color is provided.
  • An image color estimation method is an XYZ representing a device-independent color of an image displayed on a display from RGB color system gradation values (R, G, B) input to the display.
  • the second display characteristic is to set two invariant components of the gradation values (R, G, B) to fixed values and change the remaining variable components.
  • each component of the image color is expressed as the sum of the output image signal value (Rs, Gs, Bs) and the offset value (Ro, Go, Bo).
  • the offset value is calculated based on the first display characteristic, the second display characteristic, and the gradation values (R, G, B). The output image signal values (Rs, Gs, Bs) can be corrected using this offset value. Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately estimated based on the gradation values (R, G, B) input to the display. can do.
  • the offset value (Ro, Go, Bo) is a function of the gradation value (R, G, B), and the following equation (1); The function shown in FIG.
  • C1 represents a first component selected from the components of gradation values (R, G, B).
  • C2 represents a second component different from the first component selected from the components of the gradation values (R, G, B).
  • C3 represents a third component different from the first component and the second component selected from the components of the gradation values (R, G, B).
  • ⁇ Cs indicates a numerical value calculated based on a first display characteristic when the first component is an invariable component and the second component is a variable component.
  • ⁇ C's indicates a numerical value calculated based on another first display characteristic when the first component is an invariable component and the third component is a variable component.
  • the coefficient ⁇ indicates a coefficient of an approximate function that approximates one second display characteristic.
  • the coefficient ⁇ represents a coefficient of an approximate function that approximates another second display characteristic.
  • the component corresponding to the first component of the output image signal value (Rs, Gs, Bs) when only the first component is turned on, the first component, and the second component And the component corresponding to the first component of the output image signal value (Rs, Gs, Bs) when the component is turned on can be calculated by the first term of the above equation (1).
  • the component corresponding to the first component of the output image signal values (Rs, Gs, Bs) when only the first component is turned on, and the first component and the third component are turned on.
  • the difference from the component corresponding to the first component of the output image signal value (Rs, Gs, Bs) can be calculated by the second term of the above equation (1).
  • the offset value of the 1st component which considered two differences can be computed. Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately estimated based on the gradation values (R, G, B) input to the display. can do.
  • the fourth step further includes correction using a zero bias value, and the zero bias value is input to the display as a gradation value (R, G, B) in which each component is zero.
  • the measured value (Xk, Yk, Zk) of the XYZ color system displayed on the display may be used. According to this, when the gradation value (R, G, B) is zero in the display, the image color displayed on the display in the portion where the offset of the component of the image color exists can be estimated with higher accuracy. .
  • An image color estimation apparatus represents a device-independent color of an image displayed on a display from RGB color system gradation values (R, G, B) input to the display.
  • An image color estimation device for estimating an XYZ color system signal value (XYZ), based on a first display characteristic and a gradation value (R, G, B), an RGB color system signal value ( Rs, Gs, Bs), a signal value component (Rs, Gs) based on the first display characteristic, the second display characteristic, and the gradation value (R, G, B).
  • the second display characteristic is a case where two invariant components are set to fixed values among gradation values (R, G, B) and the remaining variable components are changed.
  • each component of the image color is expressed as the sum of the output image signal value (Rs, Gs, Bs) and the offset value (Ro, Go, Bo).
  • the offset value is calculated based on the first display characteristic, the second display characteristic, and the gradation values (R, G, B). The output image signal values (Rs, Gs, Bs) can be corrected using this offset value. Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately estimated based on the gradation values (R, G, B) input to the display. can do.
  • An image color estimation program represents device-independent colors of an image displayed on a display from RGB color system gradation values (R, G, B) input to the display.
  • An image color estimation program for estimating a signal value (XYZ) of an XYZ color system, the computer using an RGB color system based on a first display characteristic and a gradation value (R, G, B) Based on the signal value calculation unit for calculating the signal value (Rs, Gs, Bs), the first display characteristic, the second display characteristic, and the gradation value (R, G, B), the component ( Rs, Gs, Bs), an offset value calculator for calculating an offset value (Ro, Go, Bo) corresponding to each of the values, and adding the offset value (Ro, Go, Bo) to the signal value (Rs, Gs, Bs) As a result, correction signal values (Rc, c, Bc) and a signal value conversion unit that converts the correction signal value (Rc, Gc, Bc) into
  • the first display characteristic is that when at least one of the gradation values (R, G, B) is a variable component, the variable component of the gradation value (R, G, B) and the signal value ( Rs, Gs, Bs) is a relationship with the component corresponding to the variable component
  • the second display characteristic is that the two invariant components of the gradation values (R, G, B) are set to fixed values and the rest Is a relationship between the variable component of the gradation value (R, G, B) and the component corresponding to the invariant component of the signal value (Rs, Gs, Bs).
  • each component of the image color is expressed as the sum of the output image signal value (Rs, Gs, Bs) and the offset value (Ro, Go, Bo).
  • the offset value is calculated based on the first display characteristic, the second display characteristic, and the gradation values (R, G, B). The output image signal values (Rs, Gs, Bs) can be corrected using this offset value. Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately estimated based on the gradation values (R, G, B) input to the display. can do.
  • additive color mixture is not established, and the image color of the image displayed on the display where color tracking occurs can be estimated with high accuracy.
  • FIG. 1 shows one Embodiment of the image color estimation apparatus for performing the image color estimation method which concerns on this invention.
  • It is a block diagram which shows the functional structure of the image processing apparatus for performing the image color estimation method of 1st Embodiment.
  • It is a flowchart for demonstrating the image color estimation method of 1st Embodiment.
  • It is a figure for demonstrating a 1st display characteristic.
  • It is a figure for demonstrating a 1st display characteristic.
  • It is a figure for demonstrating a 2nd display characteristic.
  • It is a figure for demonstrating a 2nd display characteristic.
  • It is a figure for demonstrating a 2nd display characteristic.
  • FIG. 1 is a diagram illustrating a configuration of an information processing apparatus (image color estimation apparatus) 1 that executes an image color estimation program that implements the image color estimation method according to the first embodiment.
  • the information processing apparatus 1 includes a personal computer 3 (hereinafter referred to as “PC”), a display 5, a printer 7, and a colorimeter 9.
  • the information processing apparatus 1 has a function of accurately estimating the image color displayed on the display 5.
  • the estimated image color data is output to a printer 7 capable of colorimetric printing. As a result, the information processing apparatus 1 can accurately match the image color of the image displayed on the display 5 with the color of the printed matter printed by the printer 7.
  • FIG. 2 is a block diagram showing a functional configuration of the information processing apparatus 1.
  • the information processing apparatus 1 executes the image color estimation method of the first embodiment.
  • the PC 3 includes a memory 10 in which image data is stored, a graphic processing unit 20, a hard disk drive 30 (hereinafter referred to as “HDD”), and an image color estimation unit 40.
  • the memory 10 has a function of storing image data.
  • the color of the image data is represented as gradation data of the RGB color system.
  • the gradation value T of this image data is assumed to be (R, G, B).
  • R represents a component of a red gradation value T.
  • G indicates a component of a green gradation value T.
  • B represents a component of a blue tone value T.
  • the graphic processing unit 20 is connected to the memory 10 and the image color estimation unit 40.
  • the graphic processing unit 20 is connected to the display 5.
  • the graphic processing unit 20 edits the image data on the memory 10 in response to the operation of the input device by the operator.
  • the graphic processing unit 20 outputs the image data to the display 5 to display an image.
  • the graphic processing unit 20 has a function of outputting the image data on the memory 10 to the printer 7 through the image color estimating unit 40 and printing the image.
  • the graphic processor 20 has a function of outputting reference image data.
  • the reference image data is displayed on the display 5 when a first display characteristic P1 and a second display characteristic P2 described later are acquired.
  • the graphic processing unit 20 has a function of calculating a coefficient ⁇ and a coefficient ⁇ described later.
  • the graphic processing unit 20 is a functional component that is partially or entirely realized by software.
  • the HDD 30 is configured to be referred to from the image color estimation unit 40.
  • the HDD 30 stores a first display characteristic P1, a second display characteristic P2, a coefficient ⁇ , a coefficient ⁇ , and a conversion matrix M.
  • the first display characteristic P1 and the second display characteristic P2 represent the correlation between the gradation value T (R, G, B) of the image data and the signal value SG (Rs, Gs, Bs). is there.
  • the gradation value T (R, G, B) is image data input to the display 5 from the graphic processing unit 20.
  • the signal value SG (Rs, Gs, Bs) is a signal value of an actual image displayed on the display 5 based on the input image data.
  • the coefficient ⁇ and the coefficient ⁇ are coefficients for calculating an approximate curve of the second display characteristic P2.
  • the conversion matrix M is a matrix for converting RGB values of the RGB color system into tristimulus values (XYZ) of the XYZ color system.
  • the image color estimation unit 40 is connected to the HDD 30, the graphic processing unit 20, and the printer 7.
  • the image color estimation unit 40 includes a signal value calculation unit 41, an offset value calculation unit 42, a signal value correction unit 43, and a signal value conversion unit 44.
  • the image color estimation unit 40 has a function of estimating the image color displayed on the display 5. The image color is estimated based on the image data input to the display 5 from the graphic processing unit 20. Further, the image color estimation unit 40 has a function of outputting a signal value (XYZ) corresponding to the estimated color to the printer 7.
  • the image color estimation unit 40 is a functional component that is partially or entirely realized by software.
  • the signal value calculation unit 41 is connected to the graphic processing unit 20, the signal value correction unit 43, and the HDD 30.
  • the signal value calculation unit 41 has a function of calculating a signal value SG (Rs, Gs, Bs).
  • the signal value SG (Rs, Gs, Bs) is calculated based on the gradation value T (R, G, B) that is image data and the first display characteristic P1. Further, the signal value calculation unit 41 has a function of outputting the signal value SG (Rs, Gs, Bs) to the signal value correction unit 43.
  • the offset value calculation unit 42 is connected to the graphic processing unit 20, the HDD 30, and the signal value correction unit 43.
  • the offset value calculation unit 42 has a function of calculating an offset value Co (Ro, Go, Bo).
  • the offset value Co (Ro, Go, Bo) includes a gradation value T (R, G, B) that is image data, a first display characteristic P1, a second display characteristic P2, a coefficient ⁇ , It is calculated based on the coefficient ⁇ . Further, the offset value calculation unit 42 has a function of outputting the offset value Co (Ro, Go, Bo) to the signal value correction unit 43.
  • the signal value correction unit 43 is connected to the signal value calculation unit 41, the offset value calculation unit 42, and the signal value conversion unit 44.
  • the signal value correction unit 43 has a function of calculating a correction signal value SC (Rc, Gc, Bc).
  • the correction signal value SC (Rc, Gc, Bc) is calculated based on the signal value SG (Rs, Gs, Bs) and the offset value Co (Ro, Go, Bo). Further, the signal value correction unit 43 has a function of outputting the correction signal value SC (Rc, Gc, Bc) to the signal value conversion unit 44.
  • the signal value conversion unit 44 is connected to the signal value correction unit 43, the HDD 30, and the printer 7.
  • the signal value conversion unit 44 has a function of converting the correction signal value SC (Rc, Gc, Bc) of the XYZ color system into the correction signal value ST (Xc, Yc, Zc).
  • the correction signal value SC (Rc, Gc, Bc) is converted into the correction signal value ST (Xc, Yc, Zc) using the conversion matrix M.
  • the signal value conversion unit 44 has a function of outputting the correction signal value ST (Xc, Yc, Zc) to the printer 7.
  • the display 5 is connected to the graphic processing unit 20 of the PC 3.
  • the display 5 is a device that displays image data.
  • the printer 7 is connected to the image color estimation unit 40 of the PC 3.
  • the printer 7 is a device that outputs image data as a printed matter.
  • the printer 7 has a function capable of realizing colorimetric color reproduction. In colorimetric color reproduction, the signal value (Xp, Yp, Zp) of the printed image and the signal value VT (Xs, Ys, Zs) of the image displayed on the display 5 are the same. Color reproduction.
  • the colorimeter 9 is connected to the PC 3.
  • the colorimeter 9 includes a probe 9 a that is applied to the display screen of the display 5.
  • the colorimeter 9 is a device that obtains measurement values (measurement signal values VT) of an XYZ color system of an image displayed on the screen by placing the probe 9a on the display 5.
  • measurement values VT measurement signal values
  • a colorimeter for a light source such as a display colorimeter can be suitably used.
  • FIG. 3 is a flowchart showing a process of acquiring the display characteristic P.
  • the display characteristic P is created by the following process.
  • the graphic processing unit 20 generates image data.
  • the graphic processing unit 20 inputs the gradation value T1 (R, G, B), which is the generated image data, to the display 5 and displays an image composed of one color on the display 5 (S101).
  • the gradation value is represented by a numerical range of 0 to 255 (8 bits), but may be represented by another numerical range such as 0 to 1023 (10 bits).
  • the probe 9a of the colorimeter 9 is applied to the screen of the display 5, and the measurement signal value VT1 (Xs, Ys, Zs) of the image displayed on the screen is acquired (S102). .
  • the measurement signal value VT1 (Xs, Ys, Zs) is input from the colorimeter 9 to the PC 3.
  • the average of the measurement signal values VT1 may be calculated by applying the probe 9a to a plurality of locations on the screen.
  • the measurement signal value VT1 (Xs, Ys, Zs) is converted into an output image signal value SG1 (Rs, Gs, Bs) using the following equation (2) (S103).
  • the graphic processing unit 20 associates the gradation value T1 (R, G, B) and the output image signal value SG1 (Rs, Gs, Bs) and stores them in the memory 10 (S104). Subsequently, the graphic processing unit 20 changes the image data input to the display 5 and repeats the processing from S101. That is, the graphic processing unit 20 outputs the image data of the next gradation value T2 (R, G, B) to the display 5 and displays the image on the display 5 (S101). Then, the measurement signal value VT2 (Xs, Ys, Zs) is acquired by the above-described processing (S102).
  • the acquired measurement signal value VT2 (Xs, Ys, Zs) is converted into an output image signal value SG2 (Rs, Gs, Bs) (S103). Then, the graphic processing unit 20 associates the gradation value T2 (R, G, B) and the output image signal value SG2 (Rs, Gs, Bs) and stores them in the memory 10 (S104).
  • the above-described processing is repeated n times, which is a predetermined number of times, from the gradation value T1 (R, G, B) to the gradation value Tn (R, G, B) (S105).
  • n times which is a predetermined number of times, from the gradation value T1 (R, G, B) to the gradation value Tn (R, G, B) (S105).
  • the graphic processing unit 20 inputs the obtained n sets of correlations as display characteristics P to the HDD 20 and stores them as a database file (S106).
  • the acquired display characteristics P include a first display characteristic and a second display characteristic.
  • the first display characteristic is a characteristic indicating a relationship between a variable component of the gradation value T (R, G, B) and a component corresponding to the variable component of the output image signal value SG (Rs, Gs, Bs).
  • the variable component of the gradation value T (R, G, B) is at least one component among the gradation values T (R, G, B).
  • the first display characteristic includes a first display characteristic P1R (representing a relationship between R of the gradation value T (R, G, B) and Rs of the output image signal value SG (Rs, Gs, Bs).
  • R of the gradation value T (R, G, B) is a variable component.
  • the first display characteristic includes a first display characteristic P1G indicating a relationship between G of the gradation value T (R, G, B) and Gs of the output image signal value SG (Rs, Gs, Bs). (See FIG. 5).
  • G of the gradation value T (R, G, B) is a variable component.
  • the first display characteristic includes a first display characteristic P1B indicating a relationship between B of the gradation value T (R, G, B) and Bs of the output image signal value SG (Rs, Gs, Bs). (See FIG. 6).
  • B of the gradation value T (R, G, B) is a variable component.
  • the first display characteristic P1R is a relationship between the component R of the gradation value T (R, G, B) and the component Rs of the output image signal value SG (Rs, Gs, Bs).
  • the gradation value T (R, G, B) is a gradation value input to the display 5.
  • the output image signal value SG (Rs, Gs, Bs) is an output image signal value displayed on the display 5.
  • the first display characteristic P1R includes characteristics under four input conditions indicating the relationship between the gradation value T (R, G, B) and the output image signal value SG (Rs, Gs, Bs). ing.
  • the first display characteristic P1R includes a characteristic in the case where R of the gradation value T (R, G, B) is a variable component and G and B are invariable components (FIG. 4). D1). R of the gradation value T (R, G, B) is changed in the range of 0 to 255. Further, the first display characteristic P1R includes a characteristic when R and G of the gradation value T (R, G, B) are variable components and B is an invariant component (FIG. 4). D2). R and G of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • the first display characteristic P1R includes a characteristic in the case where R and B of the gradation value T (R, G, B) are variable components and G is an invariant component (FIG. 4). D3). R and B of the gradation value T (R, G, B) are changed in the range of 0 to 255. Further, the first display characteristic P1R includes a characteristic in the case where R, G, and B of the gradation value T (R, G, B) are variable components (see D4 in FIG. 4). R, G, and B of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • a curve D2 in FIG. 4 shows a relationship between R of the gradation value T (R, G, B) and Rs of the output image signal value SG (Rs, Gs, Bs).
  • the first display characteristic P1G has a relationship between the component G of the gradation value T (R, G, B) and the component Gs of the output image signal value SG (Rs, Gs, Bs).
  • the gradation value T (R, G, B) is a gradation value input to the display 5.
  • the output image signal value SG (Rs, Gs, Bs) is an output image signal value displayed on the display 5.
  • the first display characteristic P1G includes characteristics under four input conditions indicating the relationship between the gradation value T (R, G, B) and the output image signal value SG (Rs, Gs, Bs). ing.
  • the first display characteristic P1G includes a characteristic in the case where G of the gradation value T (R, G, B) is a variable component and R and B are invariant components (FIG. 5). D5). G of the gradation value T (R, G, B) is changed in the range of 0 to 255. Further, the first display characteristic P1G includes a characteristic in the case where G and B of the gradation values T (R, G, B) are variable components and R is an invariant component (FIG. 5). D6). G and B of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • the first display characteristic P1G includes a characteristic in the case where R and G of the gradation value T (R, G, B) are variable components and B is an invariant component (FIG. 5). D7). R and G of the gradation value T (R, G, B) are changed in the range of 0 to 255. Further, the first display characteristic P1G includes a characteristic in the case where R, G, and B of the gradation value T (R, G, B) are variable components (see D8 in FIG. 5). R, G, and B of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • a curve D6 in FIG. 5 shows a relationship between G of the gradation value T (R, G, B) and Gs of the output image signal value SG (Rs, Gs, Bs).
  • the first display characteristic P1B is a relationship between the component B of the gradation value T (R, G, B) and the component Bs of the output image signal value SG (Rs, Gs, Bs).
  • the gradation value T (R, G, B) is a gradation value input to the display 5.
  • the output image signal value SG (Rs, Gs, Bs) is an output image signal value displayed on the display 5.
  • the first display characteristic P1B includes characteristics under four input conditions indicating the relationship between the gradation value T (R, G, B) and the output image signal value SG (Rs, Gs, Bs). ing.
  • the first display characteristic P1B includes a characteristic when B of the gradation value T (R, G, B) is a variable component and R and G are invariant components (FIG. 6). D9). B of the gradation value T (R, G, B) is changed in the range of 0 to 255. Further, the first display characteristic P1B includes a characteristic in the case where R and B of the gradation value T (R, G, B) are variable components and G is an invariant component (FIG. 6). D10). R and B of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • the first display characteristic P1B includes a characteristic when G and B of the gradation values T (R, G, B) are variable components and R is an invariant component (FIG. 6). D11). G and B of the gradation value T (R, G, B) are changed in the range of 0 to 255. Furthermore, the first display characteristic P1B includes a characteristic in the case where R, G, and B of the gradation value T (R, G, B) are variable components (see D12 in FIG. 6). R, G, and B of the gradation value T (R, G, B) are changed in the range of 0 to 255.
  • a curve D10 in FIG. 6 shows a relationship between B of the gradation value T (R, G, B) and Bs of the output image signal value SG (Rs, Gs, Bs).
  • the range of the variable component when acquiring the first display characteristics P1R, P1G, and P1B is set to 0 to 255, but this range is an example.
  • the range of the variable component is not limited to the above-described range of 0 to 255, and may be another range.
  • the range of the variable component may be 0 to 1023.
  • the second display characteristic is a characteristic indicating the relationship between the variable component of the gradation value T (R, G, B) and the component corresponding to the invariant component of the output image signal value SG (Rs, Gs, Bs). It is. This characteristic is acquired by setting two invariant components of the gradation value T (R, G, B) to fixed values and changing the remaining variable components.
  • the second display characteristic includes, for example, a second display characteristic P2R indicating the relationship between G of the gradation value T (R, G, B) and Rs of the output image signal value SG (Rs, Gs, Bs).
  • G of the gradation value T (R, G, B) is a variable component.
  • the second display characteristic includes a second display characteristic P2G indicating a relationship between B of the gradation value T (R, G, B) and Gs of the output image signal value SG (Rs, Gs, Bs). (See FIG. 8).
  • B of the gradation value T (R, G, B) is a variable component.
  • the second display characteristic includes a second display characteristic P2B indicating a relationship between R of the gradation value T (R, G, B) and Bs of the output image signal value SG (Rs, Gs, Bs). (See FIG. 9).
  • R of the gradation value T (R, G, B) is a variable component.
  • the second display characteristic P2R is a relationship between the component G of the gradation value T (R, G, B) and the component Rs of the output image signal value SG (Rs, Gs, Bs). (See D13 in FIG. 7).
  • R of the gradation value T (R, G, B) is fixed to 128, G is changed from 0 to 255, and B is fixed to 0. In this case, the value of the component G of the gradation value T (R, G, B).
  • the second display characteristic P2G is a relationship between the component B of the gradation value T (R, G, B) and the component Gs of the output image signal value SG (Rs, Gs, Bs). (See D16 in FIG. 8).
  • the component B of the gradation value T (R, G, B) is such that R of the gradation value T (R, G, B) is fixed to 0, G is fixed to 128, and B is changed from 0 to 255. Is the value of the component B of the gradation value T (R, G, B).
  • the second display characteristic P2B is a relationship between the component R of the gradation value T (R, G, B) and the component Bs of the output image signal value SG (Rs, Gs, Bs). (See D19 in FIG. 9).
  • the R of the gradation value T (R, G, B) is changed from 0 to 255, G is fixed to 0, and B is fixed to 128. This is the value of the component R of the gradation value T (R, G, B) in this case.
  • the above-described second display characteristics P2R, P2G, and P2B are examples of the second display characteristics P2.
  • the combination of the invariant component and the variable component selected from the gradation value T (R, G, B) may be a combination different from the combination shown in the second display characteristics P2R, P2G, P2B.
  • the fixed values are 128 and 0, but may be fixed values different from 128 and 0.
  • the graphic processing unit 20 calculates the coefficient ⁇ and the coefficient ⁇ based on the second display characteristic P2 when storing the second display characteristic P2 in the HDD 30.
  • the coefficients ⁇ and ⁇ are constants for determining an approximate curve of the second display characteristic P2 as shown in FIGS. 7 to 9, for example.
  • the coefficient ⁇ 1 is a constant for setting an approximate curve of the relationship between the component G of the gradation value T (R, G, B) and the component Rs of the output image signal value SG (Rs, Gs, Bs) ( (See D14 in FIG. 7).
  • the coefficient ⁇ 1 is set so that the difference between the approximate curve (D14) and the actual display characteristic (D13) is small.
  • the coefficient ⁇ 2 is a constant for setting an approximate curve of the relationship between the component B of the gradation value T (R, G, B) and the component Gs of the output image signal value SG (Rs, Gs, Bs) ( (See D17 in FIG. 8).
  • the coefficient ⁇ 2 is set so that the difference between the approximate curve (D17) and the actual display characteristic (D16) is small.
  • the coefficient ⁇ 3 is a constant for setting an approximate curve of the relationship between the component R of the gradation value T (R, G, B) and the component Bs of the output image signal value SG (Rs, Gs, Bs) ( (See D20 in FIG. 9).
  • the coefficient ⁇ 3 is set so that the difference between the approximate curve (D20) and the actual display characteristic (D19) is small.
  • the coefficient ⁇ 1 is a constant for setting an approximate curve of the relationship between the component B of the gradation value T (R, G, B) and the component Rs of the output image signal value SG (Rs, Gs, Bs). is there.
  • the coefficient ⁇ 2 is a constant for setting an approximate curve of the relationship between the component R of the gradation value T (R, G, B) and the component Gs of the output image signal value SG (Rs, Gs, Bs).
  • the coefficient ⁇ 3 is a constant for setting an approximate curve of the relationship between the component G of the gradation value T (R, G, B) and the component Bs of the output image signal value SG (Rs, Gs, Bs).
  • the first display characteristic P1, the second display characteristic P2, the coefficient ⁇ , and the coefficient ⁇ may be acquired by the operator using the image color estimation device 1 immediately after the display 5 is introduced. , May be acquired regularly. Moreover, you may acquire every time before starting work.
  • image data composed of one color of the gradation values T1 (R, G, B) to Tn (R, G, B) is sequentially displayed on the display 5 one by one.
  • a color chart combining image data of gradation values T1 (R, G, B) to Tn (R, G, B) is displayed on one screen, and a measurement signal value is displayed for each color chart.
  • Display characteristics P may be acquired by measuring VT (Xs, Ys, Zs).
  • first display characteristic P1 and the second display characteristic P2 may be represented by measured values obtained from input values of 0 to 255. Further, the first display characteristic P1 and the second display characteristic P2 may be represented by an approximate function created based on a measured value obtained by inputting only a predetermined condition from among input values of 0 to 255.
  • FIG. 10 is a flowchart showing the steps of the image color estimation method.
  • the output image signal value SG (Rs, Gs, Bs) is calculated (S121).
  • the output image signal value SG (Rs, Gs, Bs) corresponds to the gradation value T (Rt, Gt, Bt).
  • This process is executed by the signal value calculation unit 41 (see FIG. 2).
  • Rs corresponding to Rt is calculated.
  • Rs is a component of the output image signal value SG when only red is lit with the value of Rt.
  • a curve D1 that is the first display characteristic P1R regarding red is referred to.
  • Rs associated with the value of Rt is searched on the curve D1.
  • Gs corresponding to Gt is calculated.
  • Gs is a component of the output image signal value SG when only green is lit with the value of Gt.
  • Gs associated with the value of Gt is searched on the curve D5.
  • Bs corresponding to Bt is calculated.
  • Bs is a component of the output image signal value SG when only blue is lit with the value of Bt.
  • Bs associated with the value of Bt is searched on the curve D9.
  • ⁇ Second calculation step> Subsequently, based on the gradation value T (Rt, Gt, Bt), the offset value calculation unit 42 calculates an offset value Co (Ro, Go, Bo) (S122). This process is executed by the offset value calculation unit 42 (see FIG. 2).
  • This offset value Co (Ro, Go, Bo) is calculated by the following equation (3).
  • ⁇ Cs and ⁇ C ′s are functions of the gradation value T (Rt, Gt, Bt).
  • C1, C2, and C3 are components of the gradation value T selected from the gradation values T (Rt, Gt, Bt). That is, the offset value Co (Ro, Go, Bo) is calculated as a value expressed by a function of the gradation value T (Rt, Gt, Bt).
  • ⁇ Rs in the following equation (5) is a value when R is a first component that is an invariant component and G is a second component that is a variable component.
  • ⁇ R ′s in the following formula (5) is a numerical value when R is a first component that is an invariant component and B is a third component that is a variable component.
  • rg0 (Rt) is a component Rs of an output image signal value SG of an image displayed when image data is input to the display 5.
  • R and G are Rt values
  • B is zero.
  • a curve D2 showing the display characteristics in which R and G are variable components among the display characteristics P relating to red.
  • Rs. r00 (Rt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • R is a value of Rt.
  • a curve D1 showing the display characteristics with red as a variable component among the display characteristics P relating to red.
  • Rs. r00 (Rt) is searched.
  • Rs. rg0 (Rt) is calculated.
  • rgb (Rt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • R, G, and B are Rt values.
  • a curve D4 showing display characteristics with R, G, and B as variable components is referred to.
  • Rs. rgb (Rt) is searched.
  • ⁇ Gs and ⁇ G ′s are expressed by the following formula (7).
  • ⁇ Gs in the following formula (7) is a value when G is a first component that is an invariant component and B is a second component that is a variable component.
  • ⁇ G ′s in the following formula (5) is a numerical value when G is a first component that is an invariant component and R is a third component that is a variable component.
  • Gs are Gs.
  • 0 gb (Gt) is a component Gs of an output image signal value SG of an image displayed when image data is input to the display 5.
  • G and B are Gt values, and R is zero.
  • D6 indicating a display characteristic in which G and B are variable components among the display characteristic P relating to green (see FIG. 5).
  • Gs. 0 gb (Gt) is searched.
  • Gs. 0g0 (Gt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • G is a value of Gt
  • R and B are 0.
  • a curve D5 indicating display characteristics with G as a variable component is referred to.
  • Gs. 0g0 (Gt) is searched.
  • Gs. 0 gb (Gt) to Gs. ⁇ Gs (Gt) is calculated by subtracting 0g0 (Gt).
  • Gs. rgb (Gt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • R, G, and B are Gt values.
  • a curve D8 indicating a display characteristic having R, G, and B as variable components is referred to.
  • Gs. rgb (Gt) is searched.
  • Gs. 0 gb (Gt) is obtained by a process similar to the process described above. And Gs. rgb (Gt) to Gs.
  • ⁇ G ′s (Gt) is calculated.
  • ⁇ G's (Gt) indicates a deviation of the green component Gs generated by lighting the green component and the red component at the same Gt value at the same time.
  • ⁇ G ′s is 0.
  • the coefficient ⁇ 2 is referred to. Thereby, since all the variables of the above formula (6) have been determined, Go is calculated.
  • Bo is calculated.
  • the said Formula (3) is represented like the following formula (8).
  • ⁇ Bs and ⁇ B ′s are expressed by the following formula (9).
  • ⁇ Bs in the following equation (9) is a value when B is a first component that is an invariant component and R is a second component that is a variable component.
  • ⁇ B ′s in the following equation (5) is a numerical value when B is a first component that is an invariant component and G is a third component that is a variable component.
  • Bs. r0b (Bt) is the component Bs of the output image signal value SG of the image displayed when the image data is input to the display 5.
  • R and B are values of Bt
  • G is zero.
  • R and B are variable components among the display characteristics P relating to blue (see FIG. 6).
  • Bs. r0b (Bt) is searched.
  • Bs. 00b (Bt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • B is a value of Bt
  • R and G are zero.
  • a curve D9 indicating the display characteristic with B as a variable component is referred to.
  • Bs. 00b (Bt) is searched.
  • Bs. r0b (Bt) By subtracting 00b (Bt), ⁇ Bs (Bt) is calculated.
  • Bs. rgb (Bt) is a component of an output image signal value SG of an image displayed when image data is input to the display 5.
  • R, G, and B are Bt values.
  • P2B relating to blue see FIG. 6
  • a curve D12 indicating a display characteristic having R, G, and B as variable components is referred to.
  • Bs. rgb (Bt) is searched.
  • Bs. r0b (Bt) is obtained by the same process as described above.
  • ⁇ B's (Bt) is calculated by subtracting r0b (Bt).
  • ⁇ B ′s (Bt) indicates a deviation of the blue component Bs generated by lighting the blue component and the green component at the same value at the same time.
  • ⁇ B ′s is 0.
  • the coefficient ⁇ 3 is referred to.
  • step S122 the offset value Co (Ro, Go, Bo) is calculated based on the formula (1), but the present invention is not limited to this method.
  • the measured values D13, D16, and D19 shown in FIGS. 7 to 9 may be approximated by a regression curve, and the offset value Co (Ro, Go, Bo) may be calculated from the regression curve.
  • the approximate expression in this case is different from the above expression (1).
  • the image color displayed on the display 5 can be estimated with higher accuracy.
  • the color difference of the image color estimated using the regression curve is about 1/4 to 1/5 of the color difference of the image color estimated using Equation (1).
  • the signal value correction unit 43 calculates a correction value (correction signal value SC (Rc, Gc, Bc)) of the output image signal (S123).
  • the correction signal value SC (Rc, Gc, Bc) includes the output image signal value SG (Rs, Gs, Bs) calculated in the first calculation step and the offset value Co (Ro) calculated in the second calculation step. , Go, Bo). This process is executed by the signal value correction unit 43 (see FIG. 2).
  • the correction signal value SC (Rc, Gc, Bc) is calculated by the following equation (10).
  • Rs, Gs, and Bs correspond to each component of the output image signal value SG (Rs, Gs, Bs) calculated by the signal value calculation unit 41.
  • Ro, Go, and Bo correspond to the components of the offset value Co (Ro, Go, Bo) calculated by the offset value calculation unit 42.
  • the correction signal value SC (Rc, Gc, Bc) is calculated by adding each component of the output image signal value SG and each component of the offset value Co.
  • the correction signal value SC (Rc, Gc, Bc) is converted into a correction signal value ST (Xc, Yc, Zc) (S124).
  • the correction signal value ST (Xc, Yc, Zc) is an XYZ color system signal value representing a device-independent color.
  • This process is executed by the signal value converter 44 (see FIG. 2).
  • the following formula (11) is used for the conversion.
  • M is a transformation matrix and is represented by the following equation (12).
  • the image color estimation unit 40 performs the following processing. First, an output image signal value SG (Rs, Gs, Bs) of the RGB color system is generated by the signal value calculation unit 41 based on the first display characteristic P1 and the gradation value T (R, G, B). Calculated. Next, the offset value calculation unit 42 uses the first display characteristic P1, the second display characteristic P2, and the gradation value T (R, G, B) based on the components (Rs, The calculation of the offset value Co (Ro, Go, Bo) corresponding to each of Gs, Bs) is calculated.
  • the signal value correction unit 43 adds the offset value Co (Ro, Go, Bo) to the output image signal value SG (Rs, Gs, Bs), thereby correcting the RGB color system correction signal value SC (Rc, Gc, Bc) is calculated.
  • the signal value conversion unit 44 converts the correction signal value SC (Rc, Gc, Bc) into the correction signal value ST (Xc, Yc, Zc) of the XYZ color system using the conversion matrix M.
  • each component of the image color is expressed as the sum of the output image signal value SG (Rs, Gs, Bs) and the offset value Co (Ro, Go, Bo).
  • each component of the image color is expressed as the sum of the output image signal value SG (Rs, Gs, Bs) and the offset value Co (Ro, Go, Bo). Yes.
  • the output image signal value SG (Rs, Gs, Bs) when the RGB components are turned on simultaneously and the RGB components are turned on individually. This is different from the value obtained by synthesizing the output image signal value SG (Rs, Gs, Bs).
  • this different amount is calculated based on the first display characteristic, the second display characteristic, and the gradation value T (R, G, B) as the offset value Co.
  • the output image signal value SG (Rs, Gs, Bs) can be corrected using the offset value Co. Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately determined based on the gradation value T (R, G, B) input to the display. Can be estimated.
  • the component corresponding to the first component of the output image signal value SG (Rs, Gs, Bs) when only the first component is turned on and the first
  • the difference between the component and the component corresponding to the first component of the output image signal value SG (Rs, Gs, Bs) when the component and the second component are turned on is expressed by the first term of the above equation (1). Can be calculated. Further, the component corresponding to the first component of the output image signal value SG (Rs, Gs, Bs) when only the first component is turned on, the first component, and the third component are turned on.
  • the difference from the component corresponding to the first component of the output image signal value SG (Rs, Gs, Bs) in this case can be calculated by the second term of the above equation (1). Therefore, even in a display in which additive color mixture is not established and color tracking occurs, the image color displayed on the display is accurately determined based on the gradation value T (R, G, B) input to the display. Can be estimated.
  • the offset value Co (Ro, Go, Bo) is calculated using approximate equations (Equations (3) to (9)).
  • This approximate expression is a total of 9 graphs including 6 graphs shown in FIGS. 4 to 9 plus 3 graphs obtained from data obtained when the display 5 is turned on with all the RGB components at the same gradation. Based on two graphs. That is, this approximate expression is obtained based on the first value, the second value, and the third value.
  • the first value is when two gradation values of the RGB signal are equal and constant (may be zero), and one gradation value of the RGB signal changes from zero to the maximum value.
  • the second value is one of the three gradation levels when one gradation value of the RGB signal is a constant value (may be zero), and two gradation values of the RGB signal change from zero to the maximum value. This is the display color value for the combination of key values.
  • the third value is the value of the display color when the three gradation values of the RGB signal are equal and change from zero to the maximum value.
  • FIG. 11 shows an example of the configuration of an information processing apparatus 1A used in the image color estimation method according to the second embodiment.
  • An information processing apparatus 1A shown in FIG. 11 stores a correction matrix L for correcting a zero bias value in the HDD 20 in addition to the configuration of the information processing apparatus 1 of the first embodiment shown in FIG.
  • the processing from the first calculation step (S121) to the third calculation step (S123) shown in FIG. 10 is the same as the image processing method of the first embodiment described above. It is the same.
  • a correction matrix L is created.
  • the graphic processing unit 20 outputs the image data to the display 5.
  • each component of the gradation value T (R, G, B) is all zero.
  • the measurement signal value VT (Xk, Yk, Zk) is acquired by applying the probe 9a of the colorimeter 9 to the screen of the display 5. This measured value is input from the colorimeter 9 to the PC 3.
  • Each component of this tristimulus value corresponds to each component of the correction matrix L. That is, in the case of the measurement signal value VT (Xk, Yk, Zk), the correction matrix L is the following equation (13).
  • Mc is a conversion matrix including a correction component of zero bias value.
  • an accurate correction signal value ST (Xc, Yc, Zc) of a display image is obtained. Therefore, by including the correction of the zero bias value, the image color displayed on the display 5 can be estimated with higher accuracy.
  • Example 1 the image color estimated by using the method according to the present embodiment from the predetermined image data (R, G, B) and the image data (R, G, B) are output to the display 5 and are actually displayed.
  • the color difference from the image color was confirmed.
  • the color difference between the image color estimated from the predetermined image data (R, G, B) using the method according to the comparative example and the actually displayed image color was confirmed.
  • an estimation method using SMM was selected. This SMM is used as standard in the ICC profile.
  • four different liquid crystal displays were used.
  • the gamma value of the liquid crystal display 5 is 2.2.
  • the colorimeter 9 used was TOPCONSR-3AL1.
  • the gamma value is a numerical value indicating the response characteristic of the gradation value T (R, G, B) of the image data.
  • T gradation value
  • the relationship between an input value and an output value is often not a linear function but a relationship approximate to an exponential function.
  • the exponent of the power in this exponential function is called the gamma value.
  • the information processing apparatus 1 was installed in a dark room and left for about 1 hour after turning on the power to stabilize the equipment.
  • the color target was displayed on the entire screen of the liquid crystal display 5, and the output image signal value (Xs, Ys, Zs) at the center of the liquid crystal display 5 was measured using the colorimeter 9.
  • the acquired data is 35937 colors in the sRGB color gamut.
  • the measured output image signal values (Xs, Ys, Zs) were converted into RGB color system output image signal values SG (Rs, Gs, Bs) using the above equation (2).
  • RGB color system output image signal values SG were calculated by the method according to this embodiment. Further, the output image signal value D (RpGpBp) of the RGB color system was calculated by the method according to the comparative example. Next, the color difference between the output image signal value SG (Rs, Gs, Bs) acquired by the colorimeter 9 and the output image signal value SG (Rs, Gs, Bs) calculated by the method according to the present embodiment is calculated. did. Further, the color difference between the output image signal value SG (Rs, Gs, Bs) acquired by the colorimeter 9 and the output image signal value D (RpGpBp) calculated by the method according to the comparative example was calculated.
  • FIGS. 12 (a) to 12 (d) show the comparison results of the respective liquid crystal displays.
  • FIG. 12A shows the result of the display A.
  • FIG. 12B shows the result of the display B.
  • FIG. 12C shows the result of the display C.
  • FIG. 12D shows the result of the display D.
  • the average value of the color difference was 4.80.
  • the average value of color differences was 0.63.
  • the maximum value of the color difference was 9.99.
  • the maximum color difference is 1.47.
  • the standard deviation of the color difference was 1.83, and when the method according to the present embodiment was used, the standard deviation of the color difference was 0.26. From this result, it was found that the method according to this embodiment can estimate the image color of the display with higher accuracy than the method according to the comparative example.
  • the estimated value obtained by the method according to the present embodiment is the average of the color differences than the estimated value obtained by the method according to the comparative example. Good values were obtained for all of the values, maximum values, and standard deviations. It was also found that the average color difference between the image color estimated by the method according to this embodiment and the image color actually displayed on the display is 1.0 or less. Accordingly, it has been found that the method according to the present embodiment does not establish additive color mixing and can accurately estimate the image color displayed on the display where color tracking occurs.
  • additive color mixture is not established, and the image color of the image displayed on the display where color tracking occurs can be estimated with high accuracy.

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Color Image Communication Systems (AREA)
  • Image Processing (AREA)
  • Processing Of Color Television Signals (AREA)

Abstract

画像色推定方法は、第１の表示特性、及び階調値Ｔに基づいて、信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を算出し（Ｓ１２１）、第１の表示特性Ｐ１、第２の表示特性Ｐ２、及び階調値Ｔ（Ｒ，Ｇ，Ｂ）に基づいて、信号値ＳＧの成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出する（Ｓ１２２）。そして信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）にオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を加え、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を算出する（Ｓ１２３）。補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列Ｍを用いて、補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換する（Ｓ１２４）。

Description

画像色推定方法、画像色推定装置、及び画像色推定プログラム
　本発明は、ディスプレイ等の画像表示装置に表示される画像色を推定する画像色推定方法、画像色を推定する画像色推定装置、及び画像色を推定する画像色推定プログラムに関する。
　これまで、加法混色が成立しないディスプレイの表示特性を算出する方法が提案されている。特許文献１には、カラー表示装置の表示特性を較正する表示特性較正方法に関する技術が記載されている。また、特許文献２には、ディスプレイの表示特性を測定する測定方法に関する技術が記載されている。また、特許文献３には、モニタプロファイルを作成するシステムに関する技術が記載されている。この技術によれば、個々のモニタの表示特性に応じた高精度な補正を行うモニタプロファイルが容易に作成される。
特開２００５―１２８２５４号公報 国際公開０１／０１５１２９ 特開２００５―２０８９８２号公報
Dawn Wallner, "Building ICC profiles-the Mechanics andEngineering", [online],Thesight ofInternational Color Consortium, [平成２２年１１月２２日検索], インターネット, ＜URL;http://www.color.org/icc-book1.pdf＞
　従来から、ディスプレイに表示された画像を印刷した場合に両方の画像の色を一致させたいという強い要請がある。この要請に応えるため、異なったデバイス間で色再現を行う方法としてカラーマネジメントシステム（ｃｏｌｏｒ　ｍａｎａｇｅｍｅｎｔ　ｓｙｓｔｅｍ（ＣＭＳ））がある。このＣＭＳではデバイス依存の色空間とデバイス非依存の色空間（ＸＹＺ色空間やＬａｂ色空間）とを用いて色再現を行う。ディスプレイではＲＧＢ色空間とＸＹＺ色空間との間で色変換を行う。そして、この色変換ではＩＣＣプロファイル（Ｓｈａｐｅｒ　ｍａｔｒｉｘ　Ｍｏｄｅｌ（ＳＭＭ））（非特許文献１参照。）を用いて実施することが一般的である。しかし、ＳＭＭでは、加法混色が成立し、カラートラッキング現象が生じないディスプレイを想定している。この想定のため、一般的なディスプレイでは正確な色変換が行えていない。よって、上記の要請を満たすことができない。
　そこで、本発明では、加法混色が成立せず、カラートラッキングが発生するディスプレイに表示される画像の画像色を精度良く推定できる画像色推定方法、画像色を推定する画像色推定装置、及び画像色を推定する画像色推定プログラムを提供する。
　本発明の一側面に係る画像色推定方法は、ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定方法であって、第１の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する第１の算出ステップと、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出する第２の算出ステップと、信号値（Ｒｓ，Ｇｓ，Ｂｓ）にオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する第３の算出ステップと、補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する第４の算出ステップと、を含み、第１の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の可変成分に対応する成分との関係であり、第２の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の不変成分に対応する成分との関係である。
　この画像色推定方法によれば、画像色の各成分を、出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）との和として表現している。加法混色が成立せず、カラートラッキングが発生するディスプレイでは、ＲＧＢの各成分を同時に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、ＲＧＢの各成分を個々に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を合成した値とは、相違する。この画像推定装方法では、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいてオフセット値を算出する。そして、このオフセット値を用いて出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を補正することが可能である。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　上述した画像色推定方法において、オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）は、階調値（Ｒ，Ｇ，Ｂ）の関数であり、下記式（１）；
Figure JPOXMLDOC01-appb-M000002
 に示す関数であってもよい。
　上記式（１）において、Ｃ１は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分を示す。Ｃ２は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分とは別の第２の成分を示す。Ｃ３は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分及び第２の成分とは別の第３の成分を示す。更に、ΔＣｓは第１の成分を不変成分とし、第２の成分を可変成分としたときの一の第１の表示特性に基づいて算出される数値を示す。ΔＣ’ｓは第１の成分を不変成分とし、第３の成分を可変成分としたときの他の第１の表示特性に基づいて算出される数値を示す。そして、係数αは一の第２の表示特性を近似する近似関数の係数を示す。係数βは他の第２の表示特性を近似する近似関数の係数を示す。
　この上記式（１）によれば、第１の成分のみを点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、第１の成分と第２の成分とを点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、の差異を、上記式（１）の第１項により算出することができる。更に、第１の成分のみを点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、第１の成分と第３の成分とを点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、の差異を、上記式（１）の第２項により算出することができる。そして、２つの差異を加味した第１の成分のオフセット値を算出することができる。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　上記画像色推定方法において、第４のステップでは、更にゼロバイアス値を用いた補正を含み、ゼロバイアス値は、各成分がゼロである階調値（Ｒ，Ｇ，Ｂ）をディスプレイに入力した場合に、ディスプレイに表示されるＸＹＺ表色系の測定値（Ｘｋ，Ｙｋ，Ｚｋ）であってもよい。これによれば、ディスプレイにおいて階調値（Ｒ，Ｇ，Ｂ）がゼロに対して、画像色の成分のオフセットが存在する部分のディスプレイに表示される画像色を更に精度良く推定することができる。
　本発明の別の側面に係る画像色推定装置は、ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定装置であって、第１の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する信号値算出部と、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出するオフセット値算出部と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）にオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する信号値補正部と、補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する信号値変換部と、を含み、第１の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の可変成分に対応する成分との関係であり、第２の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の不変成分に対応する成分との関係である。
　この画像色推定装置によれば、画像色の各成分を、出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）との和として表現している。加法混色が成立せず、カラートラッキングが発生するディスプレイでは、ＲＧＢの各成分を同時に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、ＲＧＢの各成分を個々に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を合成した値とは、相違する。この画像推定装方法では、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいてオフセット値を算出する。そして、このオフセット値を用いて出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を補正することが可能である。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　本発明の別の側面に係る画像色推定プログラムは、ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定プログラムであって、コンピュータを、第１の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する信号値算出部と、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいて、信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出するオフセット値算出部と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）にオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する信号値補正部と、補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する信号値変換部として機能させ、第１の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の可変成分に対応する成分との関係であり、第２の表示特性は、階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、階調値（Ｒ，Ｇ，Ｂ）の可変成分と、信号値（Ｒｓ，Ｇｓ，Ｂｓ）の不変成分に対応する成分との関係である。
　この画像色推定プログラムによれば、画像色の各成分を、出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）との和として表現している。加法混色が成立せず、カラートラッキングが発生するディスプレイでは、ＲＧＢの各成分を同時に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）と、ＲＧＢの各成分を個々に点灯させた場合の出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を合成した値とは、相違する。この画像推定装方法では、第１の表示特性、第２の表示特性、及び階調値（Ｒ，Ｇ，Ｂ）に基づいてオフセット値を算出する。そして、このオフセット値を用いて出力画像信号値（Ｒｓ，Ｇｓ，Ｂｓ）を補正することが可能である。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　本発明による画像色推定方法、画像色推定装置、及び画像色推定プログラムによれば、加法混色が成立せず、カラートラッキングが発生するディスプレイに表示される画像の画像色を精度良く推定できる。
本発明に係る画像色推定方法を実行するための画像色推定装置の一実施形態を示す図である。 第１実施形態の画像色推定方法を実行するための画像処理装置の機能的な構成を示すブロック図である。 第１実施形態の画像色推定方法を説明するためのフローチャートである。 第１の表示特性を説明するための図である。 第１の表示特性を説明するための図である。 第１の表示特性を説明するための図である。 第２の表示特性を説明するための図である。 第２の表示特性を説明するための図である。 第２の表示特性を説明するための図である。 第１実施形態の画像色推定方法を説明するためのフローチャートである。 第２実施形態の画像色推定方法を実行するための画像処理装置の機能的な構成を示すブロック図である。 実施例１を説明するための図表である。
　以下、添付図面を参照しながら本発明による画像色推定方法、画像色推定装置、及び画像色推定プログラムの実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
＜第１実施形態＞
　まず、第１の実施形態について説明する。図１は、第１実施形態の画像色推定方法を実現する画像色推定プログラムを実行する情報処理装置（画像色推定装置）１の構成を示す図である。この情報処理装置１は、パーソナルコンピュータ３（以下「ＰＣ」という）、ディスプレイ５、プリンタ７、及び測色計９を備えている。この情報処理装置１は、ディスプレイ５に表示される画像色を精度良く推定する機能を有している。推定した画像色のデータは、測色的な印刷が可能なプリンタ７に出力される。これにより、情報処理装置１は、ディスプレイ５に表示された画像の画像色と、プリンタ７により印刷された印刷物の色とを精度良く一致させることができる。
　図２は、情報処理装置１の機能的な構成を示すブロック図である。情報処理装置１は、第１実施形態の画像色推定方法を実行する。ＰＣ３は、画像データが記憶されるメモリ１０、グラフィック処理部２０、ハードディスクドライブ３０（以下「ＨＤＤ」という）、及び画像色推定部４０を備えている。メモリ１０は、画像データを格納する機能を有する。この画像データの色は、ＲＧＢ表色系の階調データとして表される。この画像データの階調値Ｔを、（Ｒ，Ｇ，Ｂ）とする。Ｒは赤色の階調値Ｔの成分を示す。Ｇは緑色の階調値Ｔの成分を示す。Ｂは青色の階調値Ｔの成分を示す。
　グラフィック処理部２０は、メモリ１０及び画像色推定部４０に接続されている。また、グラフィック処理部２０は、ディスプレイ５と接続されている。このグラフィック処理部２０は、操作者による入力装置の操作に対応して、メモリ１０上の画像データを編集する。また、グラフィック処理部２０は、上記画像データをディスプレイ５に出力し、画像を表示させる。また、グラフィック処理部２０は、画像色推定部４０を通じてメモリ１０上の画像データをプリンタ７に出力し、画像を印刷させる機能を有する。また、グラフィック処理部２０は、基準画像データを出力する機能を有する。基準画像データは、後述する第１の表示特性Ｐ１、及び第２の表示特性Ｐ２を取得する際に、ディスプレイ５に表示される。また、グラフィック処理部２０は、後述する係数α及び係数βを算出する機能を有する。なお、このグラフィック処理部２０は、一部又は全部がソフトウェア的に実現される機能的な構成要素である。
　ＨＤＤ３０は、画像色推定部４０から参照可能に構成されている。このＨＤＤ３０には、第１の表示特性Ｐ１、第２の表示特性Ｐ２、係数α、係数β、及び変換行列Ｍが格納されている。第１の表示特性Ｐ１及び第２の表示特性Ｐ２は、画像データの階調値Ｔ（Ｒ，Ｇ，Ｂ）と、信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）との相関関係を表したものである。階調値Ｔ（Ｒ，Ｇ，Ｂ）は、グラフィック処理部２０からディスプレイ５に入力される画像データである。信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、入力された画像データに基づいて、ディスプレイ５に画面表示される実際の画像の信号値である。係数α及び係数βは、第２の表示特性Ｐ２の近似曲線を算出するための係数である。変換行列Ｍは、ＲＧＢ表色系のＲＧＢ値をＸＹＺ表色系の三刺激値（ＸＹＺ）に変換するための行列である。
　画像色推定部４０は、ＨＤＤ３０、グラフィック処理部２０及びプリンタ７と接続されている。この画像色推定部４０は、信号値算出部４１、オフセット値算出部４２、信号値補正部４３、及び信号値変換部４４を備えている。画像色推定部４０は、ディスプレイ５に表示される画像色を推定する機能を有する。この画像色は、グラフィック処理部２０からディスプレイ５に入力された画像データを基に推定される。更に、画像色推定部４０は、推定した色に相当する信号値（ＸＹＺ）をプリンタ７へ出力する機能を有する。なお、この画像色推定部４０は、一部又は全部がソフトウェア的に実現される機能的な構成要素である。
　信号値算出部４１は、グラフィック処理部２０、信号値補正部４３、及びＨＤＤ３０に接続されている。この信号値算出部４１は、信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を算出する機能を有する。この信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、画像データである階調値Ｔ（Ｒ，Ｇ，Ｂ）と、第１の表示特性Ｐ１とに基づいて算出される。更に、信号値算出部４１は、信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を信号値補正部４３に出力する機能を有する。
　オフセット値算出部４２は、グラフィック処理部２０、ＨＤＤ３０、及び信号値補正部４３に接続されている。このオフセット値算出部４２は、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出する機能を有する。このオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）は、画像データである階調値Ｔ（Ｒ，Ｇ，Ｂ）と、第１の表示特性Ｐ１と、第２の表示特性Ｐ２と、係数αと、係数βとに基づいて算出される。更に、オフセット値算出部４２は、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を信号値補正部４３に出力する機能を有する。
　信号値補正部４３は、信号値算出部４１、オフセット値算出部４２、及び信号値変換部４４に接続されている。この信号値補正部４３は、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を算出する機能を有する。この補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）は、信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）とに基づいて算出される。更に、信号値補正部４３は、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を信号値変換部４４に出力する機能を有する。
　信号値変換部４４は、信号値補正部４３、ＨＤＤ３０、及びプリンタ７に接続されている。この信号値変換部４４は、ＸＹＺ表色系の補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を、補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換する機能を有する。補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）は、変換行列Mを用いて補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換される。更に、信号値変換部４４は、補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）をプリンタ７に出力する機能を有する。
　ディスプレイ５は、ＰＣ３のグラフィック処理部２０に接続されている。このディスプレイ５は、画像データを表示する装置である。プリンタ７は、ＰＣ３の画像色推定部４０に接続されている。このプリンタ７は、画像データを印刷物として出力する装置である。プリンタ７は、測色的な色再現を実現できる機能を有する。測色的な色再現とは、印刷された画像の信号値（Ｘｐ，Ｙｐ，Ｚｐ）と、ディスプレイ５に表示された画像の信号値ＶＴ（Ｘｓ，Ｙｓ，Ｚｓ）とが同一となるような色再現である。測色計９は、ＰＣ３に接続されている。測色計９は、ディスプレイ５の表示画面に当てるプローブ９ａを備えている。この測色計９は、プローブ９ａをディスプレイ５に当てて、画面上に表示された画像のＸＹＺ表色系の測定値（測定信号値ＶＴ）を取得する装置である。測色計９としては、例えば、ディスプレイ測色計など光源用の測色計を好適に用いることができる。
　次に、図３～図１０を参照し、第１実施形態における画像色推定方法について説明する。図３は、表示特性Ｐを取得する工程を示すフローチャートである。表示特性Ｐは、以下のような工程により作成される。まず、グラフィック処理部２０は、画像データを生成する。次に、グラフィック処理部２０は、生成した画像データである階調値Ｔ１（Ｒ，Ｇ，Ｂ）をディスプレイ５に入力し、ディスプレイ５において１色からなる画像を表示させる（Ｓ１０１）。
　この階調値Ｔ１（Ｒ，Ｇ，Ｂ）としては、例えば階調値Ｔ１（Ｒ＝２５５、Ｇ＝０、Ｂ＝０）がある。階調値は０～２５５（８ビット）の数値範囲で表されるが、０～１０２３（１０ビット）等の他の数値範囲で表されてもよい。この画像が表示された状態において、ディスプレイ５の画面に測色計９のプローブ９ａを当てて、画面上に表示された画像の測定信号値ＶＴ１（Ｘｓ，Ｙｓ，Ｚｓ）を取得する（Ｓ１０２）。この測定信号値ＶＴ１（Ｘｓ，Ｙｓ，Ｚｓ）は、測色計９からＰＣ３に入力される。なお、プローブ９ａを画面の複数個所に当てて測定信号値ＶＴ１（Ｘｓ，Ｙｓ，Ｚｓ）の平均を算出してもよい。この測定信号値ＶＴ１（Ｘｓ，Ｙｓ，Ｚｓ）は、下記式（２）を用いて出力画像信号値ＳＧ１（Ｒｓ，Ｇｓ，Ｂｓ）に変換される（Ｓ１０３）。
Figure JPOXMLDOC01-appb-M000003
　次に、グラフィック処理部２０は、階調値Ｔ１（Ｒ，Ｇ，Ｂ）と、出力画像信号値ＳＧ１（Ｒｓ，Ｇｓ，Ｂｓ）とを関連付けてメモリ１０に記憶させる（Ｓ１０４）。続いて、グラフィック処理部２０は、ディスプレイ５に入力する画像データを変更して、Ｓ１０１からの処理を繰り返す。すなわち、グラフィック処理部２０は、次の階調値Ｔ２（Ｒ，Ｇ，Ｂ）の画像データをディスプレイ５に出力し、ディスプレイ５に画像を表示させる（Ｓ１０１）。そして、上述した処理により、測定信号値ＶＴ２（Ｘｓ，Ｙｓ，Ｚｓ）を取得する（Ｓ１０２）。取得した測定信号値ＶＴ２（Ｘｓ，Ｙｓ，Ｚｓ）を出力画像信号値ＳＧ２（Ｒｓ，Ｇｓ，Ｂｓ）に変換する（Ｓ１０３）。そして、グラフィック処理部２０は、階調値Ｔ２（Ｒ，Ｇ，Ｂ）と出力画像信号値ＳＧ２（Ｒｓ，Ｇｓ，Ｂｓ）とを関連付けてメモリ１０に記憶させる（Ｓ１０４）。
　上述した処理を、階調値Ｔ１（Ｒ，Ｇ，Ｂ）～階調値Ｔｎ（Ｒ，Ｇ，Ｂ）について所定の回数であるｎ回繰り返す（Ｓ１０５）。上述した処理の繰り返しにより、ディスプレイ５における階調値Ｔｉ（Ｒ，Ｇ，Ｂ）（ｉ＝１～ｎ）と、出力画像信号値ＳＧｉ（Ｒｓ，Ｇｓ，Ｂｓ）との相関関係を示すカーブが得られる。そして、グラフィック処理部２０は、得られた上記ｎ組の相関関係を、表示特性ＰとしてＨＤＤ２０に入力し、データベースファイルとして保存させる（Ｓ１０６）。
　ここで、表示特性Ｐについて詳細に説明する。取得される表示特性Ｐには、第１の表示特性と、第２の表示特性がある。第１の表示特性は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の可変成分と、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の可変成分に対応する成分との関係を示す特性である。階調値Ｔ（Ｒ，Ｇ，Ｂ）の可変成分は、階調値Ｔ（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分である。
　この第１の表示特性には、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＲｓとの関係を示す第１の表示特性Ｐ１Ｒ（図４参照）がある。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲは可変成分である。また、第１の表示特性には、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＧｓとの関係を示す第１の表示特性Ｐ１Ｇ（図５参照）がある。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧは可変成分である。更に、第１の表示特性には、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＢｓとの関係を示す第１の表示特性Ｐ１Ｂ（図６参照）がある。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢは可変成分である。
　第１の表示特性Ｐ１Ｒは、図４に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓとの関係を表す。
階調値Ｔ（Ｒ，Ｇ，Ｂ）は、ディスプレイ５に入力される階調値である。出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、ディスプレイ５に表示される出力画像信号値である。この第１の表示特性Ｐ１Ｒには、階調値Ｔ（Ｒ，Ｇ，Ｂ）と出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）との関係を示す、４つの入力条件下における特性が含まれている。すなわち、第１の表示特性Ｐ１Ｒには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲを可変成分とし、ＧとＢとを不変成分とした場合の特性が含まれている（図４のＤ１参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｒには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとを可変成分とし、Ｂを不変成分とした場合の特性が含まれている（図４のＤ２参照）。
階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｒには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＢとを可変成分とし、Ｇを不変成分とした場合の特性が含まれている（図４のＤ３参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＢとは０～２５５の範囲で変化される。更に、第１の表示特性Ｐ１Ｒには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとを可変成分とした場合の特性が含まれている（図４のＤ４参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとは０～２５５の範囲で変化される。
　また、可変成分には、全て同じ値が入力される。例えば、図４のカーブＤ２は、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＲｓとの関係を示す。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲは、階調値Ｔ（Ｒ，Ｇ，Ｂ）がＲ＝ａ、Ｇ＝ａ、Ｂ＝０（ａ＝０～２５５）である画像データをディスプレイ５に入力した場合におけるａの値である。
　第１の表示特性Ｐ１Ｇは、図５に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓとの関係を表す。階調値Ｔ（Ｒ，Ｇ，Ｂ）は、ディスプレイ５に入力される階調値である。出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、ディスプレイ５に表示される出力画像信号値である。この第１の表示特性Ｐ１Ｇには、階調値Ｔ（Ｒ，Ｇ，Ｂ）と出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）との関係を示す、４つの入力条件下における特性が含まれている。すなわち、第１の表示特性Ｐ１Ｇには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧを可変成分とし、ＲとＢとを不変成分とした場合の特性が含まれている（図５のＤ５参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｇには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧとＢとを可変成分とし、Ｒを不変成分とした場合の特性が含まれている（図５のＤ６参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧとＢとは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｇには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとを可変成分とし、Ｂを不変成分とした場合の特性が含まれている（図５のＤ７参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとは０～２５５の範囲で変化される。更に、第１の表示特性Ｐ１Ｇには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとを可変成分とした場合の特性が含まれている（図５のＤ８参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとは０～２５５の範囲で変化される。
　また、可変成分には、全て同じ値が入力される。例えば、図５のカーブＤ６は、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＧｓとの関係を示す。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧは、階調値Ｔ（Ｒ，Ｇ，Ｂ）がＲ＝０、Ｇ＝ａ、Ｂ＝ａ（ａ＝０～２５５）である画像データをディスプレイ５に入力した場合におけるａの値である。
　第１の表示特性Ｐ１Ｂは、図６に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓとの関係を表す。階調値Ｔ（Ｒ，Ｇ，Ｂ）は、ディスプレイ５に入力される階調値である。出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、ディスプレイ５に表示される出力画像信号値である。この第１の表示特性Ｐ１Ｂには、階調値Ｔ（Ｒ，Ｇ，Ｂ）と出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）との関係を示す、４つの入力条件下における特性が含まれている。すなわち、第１の表示特性Ｐ１Ｂには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢを可変成分とし、ＲとＧとを不変成分とした場合の特性が含まれている（図６のＤ９参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｂには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＢとを可変成分とし、Ｇを不変成分とした場合の特性が含まれている（図６のＤ１０参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＢとは０～２５５の範囲で変化される。また、第１の表示特性Ｐ１Ｂには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧとＢとを可変成分とし、Ｒを不変成分とした場合の特性が含まれている（図６のＤ１１参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧとＢとは０～２５５の範囲で変化される。更に、第１の表示特性Ｐ１Ｂには、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとを可変成分とした場合の特性が含まれている（図６のＤ１２参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲとＧとＢとは０～２５５の範囲で変化される。
　また、可変成分には、全て同じ値が入力される。例えば、図６のカーブＤ１０は、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＢｓとの関係を示す。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢは、階調値Ｔ（Ｒ，Ｇ，Ｂ）がＲ＝ａ、Ｇ＝０、Ｂ＝ａ（ａ＝０～２５５）である画像データをディスプレイ５に入力した場合におけるａの値である。
　なお、第１の表示特性Ｐ１Ｒ，Ｐ１Ｇ，Ｐ１Ｂを取得する場合の可変成分の範囲は、０～２５５に設定しているが、この範囲は例示である。この可変成分の範囲は、上述した０～２５５の範囲に限られず、他の範囲であっても良い。例えば、可変成分の範囲は、０～１０２３の範囲であってもよい。
　また、第２の表示特性は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の可変成分と、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の不変成分に対応する成分との関係を示す特性である。この特性は、階調値Ｔ（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させて取得される。
　この第２の表示特性には、例えば、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧと出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＲｓとの関係を示す第２の表示特性Ｐ２Ｒがある（図７参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧは可変成分である。また、第２の表示特性には、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＧｓとの関係を示す第２の表示特性Ｐ２Ｇがある（図８参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢは可変成分である。更に、第２の表示特性には、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＢｓとの関係を示す第２の表示特性Ｐ２Ｂ（図９参照）等がある。階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲは可変成分である。
　第２の表示特性Ｐ２Ｒは、図７に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓとの関係を表す（図７のＤ１３参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇは、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲを１２８に固定し、Ｇを０～２５５に変化させ、Ｂを０に固定した場合における、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇの値である。上述した入力条件では、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲが１２８に固定されているために出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＲｓは、理想的には一定のはずである（図７のＤ１５参照）。しかし、加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲが一定であっても、Ｇの成分を変化させると、図７のＤ１３に示されるように出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＲｓが変化する。この理想的な特性を示すＤ１５と、実際の特性を示すＤ１３との差異は、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓのオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）の成分Ｒｏに対応する。
　第２の表示特性Ｐ２Ｇは、図８に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓとの関係を表す（図８のＤ１６参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂは、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲを０に固定し、Ｇを１２８に固定し、Ｂを０～２５５に変化させた場合における、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂの値である。加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＧが一定であっても、Ｂの成分を変化させると、図８のＤ１６に示されるように出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＧｓが変化する。理想的な特性を示すＤ１８と、実際の特性を示すＤ１６との差異は、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓのオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）の成分Ｇｏに対応する。
　第２の表示特性Ｐ２Ｂは、図９に示されるように、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓとの関係を表す（図９のＤ１９参照）。階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒは、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＲを０～２５５に変化させ、Ｇを０に固定し、Ｂを１２８に固定した場合における階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒの値である。加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）のＢが一定であっても、Ｒの成分を変化させると、図９のＤ１９に示されるように出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）のＢｓが変化する。理想的な特性を示すＤ２１と、実際の特性を示すＤ１９との差異は、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓのオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）の成分Ｂｏに対応する。
　なお、上述した第２の表示特性Ｐ２Ｒ，Ｐ２Ｇ，Ｐ２Ｂは、第２の表示特性Ｐ２の一例である。階調値Ｔ（Ｒ，Ｇ，Ｂ）から選択される不変成分及び可変成分の組み合わせは、第２の表示特性Ｐ２Ｒ，Ｐ２Ｇ，Ｐ２Ｂにおいて示した組み合わせと異なる組み合わせであってもよい。また、上述した第２の表示特性Ｐ２Ｒ，Ｐ２Ｇ，Ｐ２Ｂにおいて固定値は１２８及び０であるが、１２８及び０と異なる固定値であってもよい。
　また、グラフィック処理部２０は、第２の表示特性Ｐ２をＨＤＤ３０に保存する際に、第２の表示特性Ｐ２に基づいて、係数α及び係数βを算出する。この係数α、βは、例えば、図７～図９に示されるような第２の表示特性Ｐ２の近似曲線を決定するための定数である。
　係数α１は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓとの関係の近似曲線を設定するための定数である（図７のＤ１４参照）。係数α１は、近似曲線（Ｄ１４）と、実際の表示特性（Ｄ１３）との差異が小さくなるように設定される。係数α２は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓとの関係の近似曲線を設定するための定数である（図８のＤ１７参照）。係数α２は、近似曲線（Ｄ１７）と、実際の表示特性（Ｄ１６）との差異が小さくなるように設定される。係数α３は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓとの関係の近似曲線を設定するための定数である（図９のＤ２０参照）。係数α３は、近似曲線（Ｄ２０）と、実際の表示特性（Ｄ１９）との差異が小さくなるように設定される。また、係数β１は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓとの関係の近似曲線を設定するための定数である。係数β２は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓとの関係の近似曲線を設定するための定数である。係数β３は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇと、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓとの関係の近似曲線を設定するための定数である。
　なお、このような第１の表示特性Ｐ１、第２の表示特性Ｐ２、係数α、及び係数βは、ディスプレイ５の導入直後に操作者が画像色推定装置１を用いて取得してもよいし、定期的に取得してもよい。また、作業を開始する前に毎回取得してもよい。また、上述した手法では、階調値Ｔ１（Ｒ，Ｇ，Ｂ）～Ｔｎ（Ｒ，Ｇ，Ｂ）の１色から成る画像データを１つづつ順次ディスプレイ５に表示させている。他の方法として、例えば、階調値Ｔ１（Ｒ，Ｇ，Ｂ）～Ｔｎ（Ｒ，Ｇ，Ｂ）の画像データを組み合わせたカラーチャートを１画面で表示させ、そのカラーチャート毎に測定信号値ＶＴ（Ｘｓ，Ｙｓ，Ｚｓ）を測定し、表示特性Ｐを取得してもよい。更に、第１の表示特性Ｐ１及び第２の表示特性Ｐ２は、０～２５５の入力値により取得される測定値で表されてもよい。また、第１の表示特性Ｐ１及び第２の表示特性Ｐ２は、０～２５５の入力値のうち所定の条件のみ入力して取得した測定値を基に作成した近似関数で表されてもよい。
　次に、メモリ１０に保存された画像データがディスプレイ５に表示された場合の画像色を推定する方法を説明する。ここでは、ディスプレイ５に入力される画像データの階調値Ｔ（Ｒ，Ｇ，Ｂ）を（Ｒｔ，Ｇｔ，Ｂｔ）とする。図１０は、画像色の推定方法の工程を示すフローチャートである。
＜第１の算出ステップ＞
　まず、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を算出する（Ｓ１２１）。出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）は、階調値Ｔ（Ｒｔ，Ｇｔ，Ｂｔ）に対応する。この処理は、信号値算出部４１（図２参照）により実行される。まず、Ｒｔに対応するＲｓを算出する。ここで、Ｒｓは赤色のみをＲｔの値で点灯させた場合の出力画像信号値ＳＧの成分である。赤色に関する第１の表示特性Ｐ１ＲであるカーブＤ１が参照される。次に、カーブＤ１上で、Ｒｔの値に関連付けられたＲｓが検索される。
　次に、Ｇｔに対応するＧｓを算出する。ここで、Ｇｓは緑色のみをＧｔの値で点灯させた場合の出力画像信号値ＳＧの成分である。緑色に関する第１の表示特性Ｐ１ＧであるカーブＤ５が参照される。次に、カーブＤ５上で、Ｇｔの値に関連付けられたＧｓが検索される。
　次に、Ｂｔに対応するＢｓを算出する。ここで、Ｂｓは青色のみをＢｔの値で点灯させた場合の出力画像信号値ＳＧの成分である。青色に関する第１の表示特性Ｐ１ＢであるカーブＤ７が参照される。次に、カーブＤ９上で、Ｂｔの値に関連付けられたＢｓが検索される。
＜第２の算出ステップ＞
　続いて、階調値Ｔ（Ｒｔ，Ｇｔ，Ｂｔ）を基に、オフセット値算出部４２は、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出する（Ｓ１２２）。この処理は、オフセット値算出部４２（図２参照）により実行される。このオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）は下記式（３）により算出される。ここで、ΔＣｓ、ΔＣ’ｓは、階調値Ｔ（Ｒｔ，Ｇｔ，Ｂｔ）の関数である。Ｃ１、Ｃ２、Ｃ３は階調値Ｔ（Ｒｔ，Ｇｔ，Ｂｔ）から選択される階調値Ｔの成分である。すなわち、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）は、階調値Ｔ（Ｒｔ，Ｇｔ，Ｂｔ）の関数で表現される値として算出される。
Figure JPOXMLDOC01-appb-M000004
　詳細には、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出する場合、上記式（３）は、下記式（４）のように表される。また、式（４）中のΔＲｓ、ΔＲ’ｓは下記式（５）より表される。下記式（５）のΔＲｓは、Ｒを不変成分である第１の成分とし、Ｇを可変成分である第２の成分とした場合の値である。また、下記式（５）のΔＲ’ｓは、Ｒを不変成分である第１の成分とし、Ｂを可変成分である第３の成分とした場合の数値である。
Figure JPOXMLDOC01-appb-M000006
　上記式（５）におけるＲｓ．ｒｇ０（Ｒｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分Ｒｓである。この画像データは、ＲとＧとがＲｔの値であり、Ｂがゼロである。ここでは、赤色に関する表示特性Ｐのうち、ＲとＧとを可変成分とした表示特性を示すカーブＤ２（図４参照）が参照される。次に、カーブＤ２上において、Ｒｔの値に関連付けられたＲｓ．ｒｇ０（Ｒｔ）が検索される。
　また、Ｒｓ．ｒ００（Ｒｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＲがＲｔの値である。ここでは、赤色に関する表示特性Ｐのうち、Ｒを可変成分とした表示特性を示すカーブＤ１（図４参照）が参照される。次に、カーブＤ１上において、Ｒｔの値に関連付けられたＲｓ．ｒ００（Ｒｔ）が検索される。そして、Ｒｓ．ｒｇ０（Ｒｔ）からＲｓ．ｒ００（Ｒｔ）を減算することにより、ΔＲｓ（Ｒｔ）が算出される。上述したように、加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｒが一定であっても、他の色の成分Ｇを点灯させることにより、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｒｓは変化する。すなわち、ΔＲｓ（Ｒｔ）は、赤成分と緑成分とを同時に、同じＲｔの値で点灯させたことにより発生する赤成分の信号値Ｒｓのずれを示す（図４のＷ１参照）。
　上記式（５）におけるＲｓ．ｒｇｂ（Ｒｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＲとＧとＢとがＲｔの値である。赤色に関する表示特性Ｐのうち、ＲとＧとＢとを可変成分とした表示特性を示すカーブＤ４（図４参照）が参照される。次に、カーブＤ４上において、Ｒｔの値に関連付けられたＲｓ．ｒｇｂ（Ｒｔ）が検索される。
　次に、Ｒｓ．ｒｇ０（Ｒｔ）は、上述した処理と同様の処理により得られる。そして、Ｒｓ．ｒｇｂ（Ｒｔ）からＲｓ．ｒｇ０（Ｒｔ）を減算することにより、ΔＲ’ｓ（Ｒｔ）が算出される。ΔＲ’ｓ（Ｒｔ）は、赤成分と青成分を同時に、同じＲｔの値で点灯させたことにより発生する赤成分の信号値Ｒｓのずれを示す。しかし、図４に示されるように、赤成分と青成分とを同時に点灯させても、赤成分の信号値Ｒｓのずれは発生しない（図４のＷ２参照）。従って、ΔＲ’ｓは０である。そして、係数α１が参照される。これにより、上記式（４）の全ての変数が決定された後に、Ｒｏが算出される。
　続いて、Ｇｏを算出する。Ｇｏを算出する場合、上記式（３）は、下記式（６）のように表される。また、ΔＧｓ、ΔＧ’ｓは下記式（７）により表される。下記式（７）のΔＧｓは、Ｇを不変成分である第１の成分とし、Ｂを可変成分である第２の成分とした場合の値である。また、下記式（５）のΔＧ’ｓは、Ｇを不変成分である第１の成分とし、Ｒを可変成分である第３の成分とした場合の数値である。
Figure JPOXMLDOC01-appb-M000007
Figure JPOXMLDOC01-appb-M000008
　上記式（７）におけるＧｓ．０ｇｂ（Ｇｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分Ｇｓである。この画像データは、ＧとＢとがＧｔの値であり、Ｒがゼロである。ここでは、緑色に関する表示特性Ｐ（図５参照）のうち、ＧとＢとを可変成分とした表示特性を示すカーブＤ６が参照される。
　次に、カーブＤ６上において、Ｇｔの値に関連付けられたＧｓ．０ｇｂ（Ｇｔ）が検索される。次に、Ｇｓ．０ｇ０（Ｇｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＧがＧｔの値であり、ＲとＢとが０である。緑色に関する表示特性Ｐ（図５参照）のうち、Ｇを可変成分とした表示特性を示すカーブＤ５が参照される。次に、カーブＤ５上において、Ｇｔの値に関連付けられたＧｓ．０ｇ０（Ｇｔ）が検索される。そして、Ｇｓ．０ｇｂ（Ｇｔ）からＧｓ．０ｇ０（Ｇｔ）を減算することにより、ΔＧｓ（Ｇｔ）が算出される。上述したように、加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｇが一定であっても、他の色の成分Ｂを点灯させることにより、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｇｓは変化する。すなわち、ΔＧｓ（Ｇｔ）は、緑成分と青成分とを同時に、同じＧｔの値で点灯させたことにより発生する緑成分Ｇｓのずれを示す（図５のＷ３参照）。
　上記式（７）における、Ｇｓ．ｒｇｂ（Ｇｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＲとＧとＢとがＧｔの値である。緑色に関する表示特性Ｐ（図５参照）のうち、ＲとＧとＢとを可変成分とした表示特性を示すカーブＤ８が参照される。次に、カーブＤ８上において、Ｇｔの値に関連付けられたＧｓ．ｒｇｂ（Ｇｔ）が検索される。
　次に、Ｇｓ．０ｇｂ（Ｇｔ）は、上述した処理と同様の処理により得られる。そして、Ｇｓ．ｒｇｂ（Ｇｔ）からＧｓ．０ｇｂ（Ｇｔ）を減算することにより、ΔＧ’ｓ（Ｇｔ）が算出される。ΔＧ’ｓ（Ｇｔ）は、緑成分と赤成分とを同時に、同じＧｔの値で点灯させたことにより発生する緑成分Ｇｓのずれを示す。しかし、図５に示されるように、緑成分と赤成分とを同時に点灯させても、緑成分Ｇｓのずれは発生しない（図５のＷ４参照）。従って、ΔＧ’ｓは０である。そして、係数α２が参照される。これにより、上記式（６）の全ての変数が決定されたため、Ｇｏが算出される。
　続いて、Ｂｏを算出する。Ｂｏを算出する場合、上記式（３）は、下記式（８）のように表される。ΔＢｓ及びΔＢ’ｓは下記式（９）により表される。下記式（９）のΔＢｓは、Ｂを不変成分である第１の成分とし、Ｒを可変成分である第２の成分とした場合の値である。また、下記式（５）のΔＢ’ｓは、Ｂを不変成分である第１の成分とし、Ｇを可変成分である第３の成分とした場合の数値である。
Figure JPOXMLDOC01-appb-M000009
Figure JPOXMLDOC01-appb-M000010
　上記式（９）におけるＢｓ．ｒ０ｂ（Ｂｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分Ｂｓである。この画像データは、ＲとＢとがＢｔの値であり、Ｇがゼロである。ここでは、青色に関する表示特性Ｐ（図６参照）のうち、ＲとＢとを可変成分とした表示特性を示すカーブＤ１０が参照される。
　次に、カーブＤ１０上において、Ｂｔの値に関連付けられたＢｓ．ｒ０ｂ（Ｂｔ）が検索される。次に、Ｂｓ．００ｂ（Ｂｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＢがＢｔの値であり、ＲとＧとがゼロである。青色に関する表示特性Ｐ（図６参照）のうち、Ｂを可変成分とした表示特性を示すカーブＤ９が参照される。次に、カーブＤ９上において、Ｂｔの値に関連付けられたＢｓ．００ｂ（Ｂｔ）が検索される。そして、Ｂｓ．ｒ０ｂ（Ｂｔ）からＢｓ．００ｂ（Ｂｔ）を減算することにより、ΔＢｓ（Ｂｔ）が算出される。上述したように、加法混色が成立せず、カラートラッキングが発生するディスプレイ５においては、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分Ｂが一定であっても、他の色の成分Ｒを点灯させることにより、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の成分Ｂｓは変化する。すなわち、ΔＢｓ（Ｂｔ）は、青成分と赤成分とを同時に、成分Ｂｔの値で点灯させたことにより発生する青成分の信号値Ｂｓのずれを示す（図６のＷ５参照）。
　上記式（９）におけるＢｓ．ｒｇｂ（Ｂｔ）は、画像データをディスプレイ５に入力したときに表示される画像の出力画像信号値ＳＧの成分である。この画像データは、ＲとＧとＢとがＢｔの値である。青色に関する第２の表示特性Ｐ２Ｂ（図６参照）のうち、ＲとＧとＢとを可変成分とした表示特性を示すカーブＤ１２が参照される。
　次に、カーブＤ１２上において、Ｂｔの値に関連付けられたＢｓ．ｒｇｂ（Ｂｔ）が検索される。Ｂｓ．ｒ０ｂ（Ｂｔ）は、上述した処理と同様の処理により得られる。そして、Ｂｓ．ｒｇｂ（Ｂｔ）からＧｓ．ｒ０ｂ（Ｂｔ）を減算することにより、ΔＢ’ｓ（Ｂｔ）が算出される。ΔＢ’ｓ（Ｂｔ）は、青成分と緑成分とを同時に、同じ値で点灯させたことにより発生する青成分Ｂｓのずれを示す。しかし、図６に示されるように、青成分と緑成分とを同時に点灯させても、青成分Ｂｓのずれは発生しない（図６のＷ６参照）。従って、ΔＢ’ｓは０である。そして、係数α３が参照される。これにより、式（８）の全ての変数が決定されたため、Ｂｏが算出される。
　なお、工程Ｓ１２２では、式（１）に基づいてオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出したが、この方法に限定されることはない。例えば、図７～図９に示された測定値Ｄ１３，Ｄ１６，Ｄ１９を回帰曲線で近似して、該回帰曲線からオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出してもよい。この場合の近似式は上記式（１）とは異なる。このような方法によれば、ディスプレイ５に表示された画像色を更に精度よく推定することができる。例えば、回帰曲線を用いて推定された画像色の色差は、式（１）を用いて推定された画像色の色差の１／４～１／５程度である。
＜第３の算出ステップ＞
　続いて、信号値補正部４３は、出力画像信号の補正値（補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ））を算出する（Ｓ１２３）。補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）は、第１の算出ステップで算出された出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、第２の算出ステップで算出されたオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）とを基に算出される。この処理は、信号値補正部４３（図２参照）により実行される。補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）は、下記式（１０）により算出される。
Figure JPOXMLDOC01-appb-M000011
　ここで、Ｒｓ，Ｇｓ，Ｂｓは、信号値算出部４１において算出された出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の各成分に対応する。また、Ｒｏ、Ｇｏ、Ｂｏは、オフセット値算出部４２において算出されたオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）の各成分に対応する。出力画像信号値ＳＧの各成分と、オフセット値Ｃｏの各成分を加算することにより、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）が算出される。
＜第４の算出ステップ＞
　続いて、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）を、補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換する（Ｓ１２４）。補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）は、デバイス非依存の色を表すＸＹＺ表色系の信号値である。この処理は、信号値変換部４４（図２参照）により実行される。変換には、下記式（１１）が用いられる。ここで、Ｍは変換行列であり、下記式（１２）により示される。この行列の要素は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分が全て最高輝度（Ｒ＝２５５、Ｇ＝２５５、Ｂ＝２５５）であるときの三刺激値（ＸＹＺ）である。上述の手順により、メモリ１０に保存された画像データがディスプレイ５に表示された場合の画像色が推定される。
Figure JPOXMLDOC01-appb-M000012
Figure JPOXMLDOC01-appb-M000013
　以上説明したような第１実施形態による画像色推定方法では、画像色推定部４０により、以下の処理が行われる。まず、信号値算出部４１により、第１の表示特性Ｐ１、及び階調値Ｔ（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）が算出される。次に、オフセット値算出部４２により、第１の表示特性Ｐ１、第２の表示特性Ｐ２、及び階調値Ｔ（Ｒ，Ｇ，Ｂ）に基づいて、出力画像信号値ＳＧの成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出が算出される。次に、信号値補正部４３により、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）にオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）が算出される。そして、信号値変換部４４により、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）が、変換行列Ｍを用いて、ＸＹＺ表色系の補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換される。第１実施形態では、画像色の各成分を、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）との和として表現している。
　第１実施形態に係る画像色推定方法では、画像色の各成分を、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、オフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）との和として表現している。加法混色が成立せず、カラートラッキングが発生するディスプレイでは、ＲＧＢの各成分を同時に点灯させた場合の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、ＲＧＢの各成分を個々に点灯させて出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を合成した値とは、相違する。第１実施形態の画像推定装置では、この相違する量を、オフセット値Ｃｏとして第１の表示特性、第２の表示特性、及び階調値Ｔ（Ｒ，Ｇ，Ｂ）に基づいて算出し、そのオフセット値Ｃｏを用いて出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を補正することが可能である。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値Ｔ（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　また、第１実施形態による画像色推定方法では、第１の成分のみを点灯させた場合の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、第１の成分と第２の成分とを点灯させた場合の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、の差異を、上記式（１）の第１項により算出することができる。更に、第１の成分のみを点灯させた場合の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、第１の成分と第３の成分とを点灯させた場合の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の第１の成分に相当する成分と、の差異を、上記式（１）の第２項により算出することができる。従って、加法混色が成立せず、カラートラッキングが発生するディスプレイであっても、ディスプレイに入力される階調値Ｔ（Ｒ，Ｇ，Ｂ）を基に、ディスプレイに表示される画像色を精度よく推定することができる。
　また、第１実施形態による画像色推定方法では、必要なデータ量を削減しつつ、所望の推定精度を得ることができる。所望の推定精度とは、例えば色差が１以下である。すなわち、第１実施形態による画像色推定方法では、近似式（式（３）～式（９））を用いてオフセット値Ｃｏ（Ｒｏ，Ｇｏ，Ｂｏ）を算出する。この近似式は、図４～図９に示された６つのグラフに加え、ＲＧＢの各成分を全て同じ階調でディスプレイ５を点灯させたときのデータから得られる３つのグラフを加えた合計９つのグラフを基にして得られる。すなわち、この近似式は、第１の値と、第２の値と、第３の値とを基にして得られる。ここで、第１の値とは、ＲＧＢ信号のうち２つの階調値が等しくかつ一定値（ゼロでもよい）であり、ＲＧＢ信号のうち１つの階調値がゼロから最大値まで変化するときにおける３種の階調値の組み合わせに対する表示色の値である。第２の値とは、ＲＧＢ信号のうち１つの階調値が一定値（ゼロでもよい）であり、ＲＧＢ信号のうち２つの階調値がゼロから最大値まで変化するときにおける３種の階調値の組み合わせに対する表示色の値である。第３の値とは、ＲＧＢ信号の３つの階調値が等しくかつゼロから最大値まで変化するときにおける表示色の値である。
＜第２実施形態＞
　次に、第２実施形態について説明する。第２実施形態と、上述した第１実施形態とで異なる点は、第４のステップにおいて、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）に対して、ＸＹＺ表色系の補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）に変換すると共に、ゼロバイアス値の補正を行う点である。
　液晶ディスプレイのような画像表示装置では、各色の成分が全てゼロの階調値Ｔ（Ｒ，Ｇ，Ｂ）を入力して画面上の画像色を測定した場合、出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）の各成分はゼロにならない。すなわち、画面上の画像色を測定すると、所定の測定信号値ＶＴ（Ｘｋ，Ｙｋ，Ｚｋ）が得られる。これをゼロバイアス値と呼ぶことにする。図１１は、第２実施形態に係る画像色推定方法に用いられる情報処理装置１Ａの構成の一例を示す。図１１に示される情報処理装置１Ａは、図２に示される第１実施形態の情報処理装置１の構成に加え、更にＨＤＤ２０にゼロバイアス値を補正するための補正行列Ｌを格納している。
　次に、第２実施形態に係る画像色推定方法を説明する。第２実施形態に係る画像色推定方法は、図１０に示される第１の算出ステップ（Ｓ１２１）から第３の算出ステップ（Ｓ１２３）までの処理は、上述した第１実施形態の画像処理方法と同様である。まず、補正行列Ｌを作成する。グラフィック処理部２０は、画像データを、ディスプレイ５に出力する。この画像データは、階調値Ｔ（Ｒ，Ｇ，Ｂ）の各成分が全てゼロである。次に、ディスプレイ５の画面に測色計９のプローブ９ａを当てて、測定信号値ＶＴ（Ｘｋ，Ｙｋ，Ｚｋ）を取得する。この測定値は、測色計９からＰＣ３に入力される。この三刺激値の各成分が、補正行列Ｌの各成分に対応する。すなわち、測定信号値ＶＴ（Ｘｋ，Ｙｋ，Ｚｋ）の場合、補正行列Ｌは、下記式（１３）である。
Figure JPOXMLDOC01-appb-M000014
＜第４の算出ステップ＞
　次に、補正信号値ＳＣ（Ｒｃ，Ｇｃ，Ｂｃ）をＸＹＺ値に変換する処理と併せて、ゼロバイアス値の補正を行う。この処理は、信号値変換部４４（図１１参照）により実行される。変換には、下記式（１４）が用いられる。ここで、Ｍｃはゼロバイアス値の補正成分を含んだ変換行列である。Ｍｃは下記式（１５）により示される。この行列の要素は、出力画像信号値（ＸＹＺ）から、補正行列Ｌの各成分を減算したものである。この出力画像信号値（ＸＹＺ）は、階調値Ｔ（Ｒ，Ｇ，Ｂ）の成分が全て最高輝度（Ｒ＝２５５、Ｇ＝２５５、Ｂ＝２５５）である。上述の手順により、ゼロバイアス補正がなされた画像色が算出される。
Figure JPOXMLDOC01-appb-M000015
Figure JPOXMLDOC01-appb-M000016
　第２実施形態の画像色推定方法では、表示画像の正確な補正信号値ＳＴ（Ｘｃ，Ｙｃ，Ｚｃ）が得られる。従って、ゼロバイアス値の補正を含んでいることで、ディスプレイ５に表示される画像色を更に精度良く推定することができる。
＜実施例１＞
　次に、所定の画像データ（Ｒ，Ｇ，Ｂ）から本実施形態に係る方法を用いて推定した画像色と、画像データ（Ｒ，Ｇ，Ｂ）をディスプレイ５に出力し、実際に表示された画像色との色差を確認した。また、所定の画像データ（Ｒ，Ｇ，Ｂ）から比較例に係る方法を用いて推定した画像色と、実際に表示された画像色との色差を確認した。比較例に係る方法として、ＳＭＭを用いて推定する方法を選択した。このＳＭＭは、ＩＣＣプロファイルにおいて標準的に使用されている。また、評価にはそれぞれ異なる４台の液晶ディスプレイを用いた。この液晶ディスプレイ５のガンマ値は２．２である。また、測色計９は、ＴＯＰＣＯＮＳＲ－３ＡＬ１を用いた。なお、ガンマ値とは、画像データの階調値Ｔ（Ｒ，Ｇ，Ｂ）の応答特性を示す数値をいう。画像表示装置では、入力値と出力値の関係は一次関数でなく、指数関数に近似した関係であることが多い。この指数関数におけるべき乗の指数をガンマ値という。
　まず、情報処理装置１を暗室に設置し、機器を安定させるために電源を投入後、約１時間放置した。次に、液晶ディスプレイ５に色標を全画面表示し、測色計９を用いて液晶ディスプレイ５の中心における出力画像信号値（Ｘｓ，Ｙｓ，Ｚｓ）を測定した。なお、取得したデータはｓＲＧＢ色域内の３５９３７色である。測定された出力画像信号値（Ｘｓ，Ｙｓ，Ｚｓ）は、上記式（２）を用いてＲＧＢ表色系の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）に変換した。
　次に、本実施形態に係る方法により、ＲＧＢ表色系の出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）を算出した。また、比較例に係る方法により、ＲＧＢ表色系の出力画像信号値Ｄ（ＲｐＧｐＢｐ）を算出した。次に、測色計９により取得した出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、本実施形態に係る方法により算出した出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）との色差を算出した。また、測色計９により取得した出力画像信号値ＳＧ（Ｒｓ，Ｇｓ，Ｂｓ）と、比較例に係る方法により算出した出力画像信号値Ｄ（ＲｐＧｐＢｐ）との色差を算出した。
　図１２（ａ）～（ｄ）は、それぞれの液晶ディスプレイにおける比較の結果を示す。図１２（ａ）はディスプレイＡの結果を示す。図１２（ｂ）はディスプレイＢの結果を示す。図１２（ｃ）はディスプレイＣの結果を示す。図１２（ｄ）はディスプレイＤの結果を示す。図１２（ａ）を参照すると、比較例に係る方法を用いた場合は色差の平均値が４．８０であった。一方、本実施形態に係る方法を用いた場合は色差の平均値が０．６３であった。また、比較例に係る方法を用いた場合は色差の最大値が９．９９であった。一方、本実施形態に係る方法を用いた場合は色差の最大値が１．４７であった。更に、比較例に係る方法を用いた場合は色差の標準偏差が１．８３であり、本実施形態に係る方法を用いた場合は色差の標準偏差が０．２６であった。この結果により、本実施形態に係る方法は、比較例に係る方法よりもディスプレイの画像色を精度良く推定できることがわかった。
　上記したディスプレイＡに関する結果と同様に、他のディスプレイＢ、Ｃ、Ｄについても、比較例に係る方法により得た推定値よりも、本実施形態に係る方法により得た推定値が、色差の平均値、最大値、及び標準偏差の全てにおいて良好な値を示した。また、本実施形態に係る方法により推定した画像色と、実際にディスプレイに表示された画像色との色差の平均は、１．０以下であることがわかった。従って、本実施形態に係る方法は、加法混色が成立せず、カラートラッキングが発生するディスプレイに表示されている画像色を精度良く推定できることがわかった。
　本発明による画像色推定方法、画像色推定装置、及び画像色推定プログラムによれば、加法混色が成立せず、カラートラッキングが発生するディスプレイに表示される画像の画像色を精度良く推定できる。
１…情報処理装置、３…パーソナルコンピュータ、５…ディスプレイ、７…プリンタ、９…測色計、９ａ…プローブ、１０…メモリ、２０…グラフィック処理部、３０…ハードディスクドライブ、４０…画像色推定部、４１…信号値算出部、４２…オフセット値算出部、４３…信号値補正部、４４…信号値変換部、Ｄ１～Ｄ２１…カーブ、P…表示特性、Ｐ１…第１の表示特性、Ｐ２…第２の表示特性、α，β…係数、Ｍ，Ｍｃ…変換行列、Ｌ…補正行列。

Claims (6)

  1. 　ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、前記ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定方法であって、
    　第１の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する第１の算出ステップと、
    　前記第１の表示特性、第２の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、前記信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出する第２の算出ステップと、
    　前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）に前記オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する第３の算出ステップと、
    　前記補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、前記信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する第４の算出ステップと、を含み、
    　前記第１の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記可変成分に対応する成分との関係であり、
    　前記第２の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記不変成分に対応する成分との関係である、画像色推定方法。
  2. 　前記オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）は、前記階調値（Ｒ，Ｇ，Ｂ）の関数である、請求項１に記載の画像色推定方法。
  3. 　前記オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）は、下記式（１）；
    Figure JPOXMLDOC01-appb-M000001
     　Ｃ１は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分、
    　Ｃ２は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分とは別の第２の成分、
    　Ｃ３は階調値（Ｒ，Ｇ，Ｂ）の成分から選択される第１の成分及び第２の成分とは別の第３の成分、
    　ΔＣｓは前記第１の成分を前記不変成分とし、前記第２の成分を前記可変成分としたときの一の前記第１の表示特性に基づいて算出される数値、
    　ΔＣ’ｓは前記第１の成分を前記不変成分とし、前記第３の成分を前記可変成分としたときの他の前記第１の表示特性に基づいて算出される数値、
    　αは一の前記第２の表示特性を近似する近似関数の係数、
    　βは他の前記第２の表示特性を近似する近似関数の係数、に示す関数である、請求項２に記載の画像色推定方法。
  4. 　前記第４のステップでは、更にゼロバイアス値を用いた補正を含み、
    　前記ゼロバイアス値は、各成分がゼロである前記階調値（Ｒ，Ｇ，Ｂ）を前記ディスプレイに入力した場合に、前記ディスプレイに表示されるＸＹＺ表色系の信号値（Ｘｋ，Ｙｋ，Ｚｋ）である、請求項１～請求項３のいずれか一項に記載の画像色推定方法。
  5. 　ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、前記ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定装置であって、
    　第１の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する信号値算出部と、
    　前記第１の表示特性、第２の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、前記信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出するオフセット値算出部と、
    　前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）に前記オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する信号値補正部と、
    　前記補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、前記信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する信号値変換部と、を含み、
    　前記第１の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記可変成分に対応する成分との関係であり、
    　前記第２の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記不変成分に対応する成分との関係である、画像色推定装置。
  6. 　ディスプレイに入力されるＲＧＢ表色系の階調値（Ｒ，Ｇ，Ｂ）から、前記ディスプレイに表示される画像のデバイス非依存の色を表すＸＹＺ表色系の信号値（ＸＹＺ）を推定する画像色推定プログラムであって、
    　コンピュータを、
    　第１の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、ＲＧＢ表色系の信号値（Ｒｓ，Ｇｓ，Ｂｓ）を算出する信号値算出部と、
    　前記第１の表示特性、第２の表示特性、及び前記階調値（Ｒ，Ｇ，Ｂ）に基づいて、前記信号値の成分（Ｒｓ，Ｇｓ，Ｂｓ）のそれぞれに対応するオフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を算出するオフセット値算出部と、
    　前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）に前記オフセット値（Ｒｏ，Ｇｏ，Ｂｏ）を加えることにより、ＲＧＢ表色系の補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を算出する信号値補正部と、
    　前記補正信号値（Ｒｃ，Ｇｃ，Ｂｃ）を、変換行列を用いて、前記信号値（Ｘｃ，Ｙｃ，Ｚｃ）に変換する信号値変換部として機能させ、
    　前記第１の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、少なくとも１つの成分を可変成分とした場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記可変成分に対応する成分との関係であり、
    　前記第２の表示特性は、前記階調値（Ｒ，Ｇ，Ｂ）うち、２つの不変成分を固定値に設定し、残りの可変成分を変化させた場合の、前記階調値（Ｒ，Ｇ，Ｂ）の前記可変成分と、前記信号値（Ｒｓ，Ｇｓ，Ｂｓ）の前記不変成分に対応する成分との関係である、画像色推定プログラム。
PCT/JP2011/077687 2010-11-30 2011-11-30 画像色推定方法、画像色推定装置、及び画像色推定プログラム Ceased WO2012074014A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US13/989,978 US9311885B2 (en) 2010-11-30 2011-11-30 Image color estimation method, image color estimation device, and image color estimation program
JP2012546917A JP5825600B2 (ja) 2010-11-30 2011-11-30 画像色推定方法、画像色推定装置、及び画像色推定プログラム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010267389 2010-11-30
JP2010-267389 2010-11-30

Publications (1)

Publication Number Publication Date
WO2012074014A1 true WO2012074014A1 (ja) 2012-06-07

Family

ID=46171941

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2011/077687 Ceased WO2012074014A1 (ja) 2010-11-30 2011-11-30 画像色推定方法、画像色推定装置、及び画像色推定プログラム

Country Status (3)

Country Link
US (1) US9311885B2 (ja)
JP (1) JP5825600B2 (ja)
WO (1) WO2012074014A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102349384B1 (ko) * 2017-03-15 2022-01-11 삼성전자주식회사 투명 디스플레이 장치, 그 제어 방법 및 컴퓨터 판독가능 기록 매체

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08294009A (ja) * 1995-04-20 1996-11-05 Canon Inc 画像処理装置及びその方法、及びシステム
JP2008515340A (ja) * 2004-10-01 2008-05-08 キヤノン株式会社 プロジェクタのカラーキャラクタリゼーション

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001015129A1 (en) 1999-08-25 2001-03-01 Fujitsu Limited Display measuring method and profile preparing method
US6844881B1 (en) * 2002-03-29 2005-01-18 Apple Computer, Inc. Method and apparatus for improved color correction
US20060050079A1 (en) * 2002-10-31 2006-03-09 Moshe Ben-Chorin System and method of selective adjustment of a color display
JP3751621B2 (ja) 2003-10-23 2006-03-01 株式会社ナナオ 表示特性較正方法、表示特性較正装置及びコンピュータプログラム
CN100533543C (zh) * 2003-12-12 2009-08-26 松下电器产业株式会社 彩色图像显示装置、色变换装置、色彩模拟装置及方法
JP2005208982A (ja) 2004-01-23 2005-08-04 Dainippon Printing Co Ltd モニタプロファイル作成システム、プログラム及び記録媒体
US7668365B2 (en) * 2004-03-08 2010-02-23 Seiko Epson Corporation Determination of main object on image and improvement of image quality according to main object
US9330630B2 (en) * 2008-08-30 2016-05-03 Sharp Laboratories Of America, Inc. Methods and systems for display source light management with rate change control
JP2010281888A (ja) * 2009-06-02 2010-12-16 Seiko Epson Corp 集積回路装置、電気光学装置及び電子機器

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08294009A (ja) * 1995-04-20 1996-11-05 Canon Inc 画像処理装置及びその方法、及びシステム
JP2008515340A (ja) * 2004-10-01 2008-05-08 キヤノン株式会社 プロジェクタのカラーキャラクタリゼーション

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NOBUHIKO TAMURA ET AL.: "Calibration of LCD Colorimetry Considering Channel Interaction", THE JOURNAL OF THE INSTITUTE OF IMAGE INFORMATION AND TELEVISION ENGINEERS, vol. 56, no. 8, 2002, pages 1315 - 1320 *
TOSHIYUKI WAKEBE ET AL.: "Color Reproduction Method Reducing Reproduction Error Induced by Color Tracking, Device Drift and Imperfect Additivity", IEICE TECHNICAL REPORT, vol. 105, no. 617, 17 February 2006 (2006-02-17), pages 37 - 40 *

Also Published As

Publication number Publication date
JP5825600B2 (ja) 2015-12-02
JPWO2012074014A1 (ja) 2014-05-19
US20130328911A1 (en) 2013-12-12
US9311885B2 (en) 2016-04-12

Similar Documents

Publication Publication Date Title
CN100583945C (zh) 图像处理装置、图像处理方法
JP6108933B2 (ja) キャリブレーション装置及びその制御方法
JP4623137B2 (ja) 色処理装置、方法及びプログラム
JP5037311B2 (ja) カラー再現システム及び方法
US20100309218A1 (en) Method for color calibration and device using the same
JP2012023627A (ja) 色変換係数生成装置、色再現特性修正装置及びプログラム
CN115035835B (zh) 一种屏幕色域的校准方法及其装置
US20070139678A1 (en) Profile creating apparatus, profile creating program storage medium, and image output apparatus
CN102511050B (zh) 用于优化显示器配置文件的方法
KR20070043535A (ko) 모니터의 색특성을 보정하는 방법 및 장치
WO1999056088A1 (en) Method to estimate the white point on a display device
JP4938669B2 (ja) 色処理方法およびその装置
US8687013B2 (en) Method, apparatus and system for providing a color device characterization with a quality evaluation
JP5825600B2 (ja) 画像色推定方法、画像色推定装置、及び画像色推定プログラム
JP2008177783A (ja) 色変換装置及びプログラム
US7847986B2 (en) Multiband color management
JP3815542B2 (ja) 色変換装置、色変換係数算出装置、色変換方法、及び色変換係数算出方法
JP2011205467A (ja) 色処理装置及びプログラム
JP2002247396A (ja) プリンターの色補正方法及びこの装置
WO2007083717A1 (ja) 色変換マトリクス作成方法及び色変換方法
JP2011197428A (ja) 階調補正装置、ディスプレイおよび階調補正方法
JPH11298746A (ja) 色信号変換テーブル作成方法および色信号変換方法
JP6753238B2 (ja) 指標表示装置および指標表示プログラム
JP5651663B2 (ja) ３ｄｌｕｔ生成方法、３ｄｌｕｔ生成装置及びコンピュータプログラム
JP2006217377A (ja) ディスプレイ較正方法及びディスプレイ較正用プログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11844987

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2012546917

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 13989978

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 11844987

Country of ref document: EP

Kind code of ref document: A1