WO2007061117A1 - 色変換マトリクス作成方法及び色変換方法 - Google Patents

色変換マトリクス作成方法及び色変換方法 Download PDF

Info

Publication number
WO2007061117A1
WO2007061117A1 PCT/JP2006/323741 JP2006323741W WO2007061117A1 WO 2007061117 A1 WO2007061117 A1 WO 2007061117A1 JP 2006323741 W JP2006323741 W JP 2006323741W WO 2007061117 A1 WO2007061117 A1 WO 2007061117A1
Authority
WO
WIPO (PCT)
Prior art keywords
color
color conversion
image data
conversion matrix
color space
Prior art date
Application number
PCT/JP2006/323741
Other languages
English (en)
French (fr)
Inventor
Yoshifumi Shimodaira
Ken Itou
Original Assignee
National University Corporation Shizuoka University
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 National University Corporation Shizuoka University filed Critical National University Corporation Shizuoka University
Priority to JP2007546535A priority Critical patent/JP4806778B2/ja
Publication of WO2007061117A1 publication Critical patent/WO2007061117A1/ja

Links

Classifications

    • 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
    • H04N1/6033—Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer using test pattern analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/001—Texturing; Colouring; Generation of texture or colour
    • 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
    • G—PHYSICS
    • G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09G—ARRANGEMENTS OR CIRCUITS FOR CONTROL OF INDICATING DEVICES USING STATIC MEANS TO PRESENT VARIABLE INFORMATION
    • G09G2340/00—Aspects of display data processing
    • G09G2340/06—Colour space transformation

Definitions

  • the present invention relates to a color conversion matrix creation method and a color conversion method, and in particular, a color conversion matrix creation method when performing color conversion by matrix calculation, and a color conversion created by the color conversion matrix creation method Background art on color conversion method for color conversion of input image using matrix
  • printer color conversion is based on the principle of subtractive color mixing, paper, and ink. It is difficult to formulate a theoretical formula because the characteristics such as are non-linear.
  • ICC profile International Color Consortium
  • Non-Patent Document 2 a method using an ICC profile (ICC: International Color Consortium) that describes device color information is standardized as color conversion of input / output devices (see Non-Patent Document 2, for example).
  • An ICC profile applied to a printer is based on a LUT (Look-uptable) in which a plurality of input / output relationships are associated, and performs color conversion by combining the LUT and a table interpolation method (for example, non-display). (See Patent Document 3).
  • LUT Look-uptable
  • table interpolation method for example, non-display
  • Patent Documents 1 and 2 disclose a technique for creating a color conversion table and a color conversion matrix by overlapping each partial region.
  • Patent Document 1 Japanese Patent Laid-Open No. 2002-112053
  • Patent Document 2 JP 2000-188695 A
  • Patent Document 3 Japanese Patent Application Laid-Open No. 2004-266590
  • Non-Patent Document 1 Koji Tashiro: “To think about using an inkjet printer to produce the right colors”, Journal of the Japan Society of Photography, vol.68, No.l, pp.70-72 (2005)
  • Non-Patent Document 3 Dawn Wallner: "Building ICC profiles- the Mechanics and Engineering ⁇ (2000)
  • Non-Patent Document 4 Hiroshi Kodera, Tomomi Ishige, Ryoichi Saito: High-precision color matching of input / output devices by dividing amber space ", Journal of the Institute of Image Electronics Engineers of Japan, vol.29, pp.128-139 (2000)
  • Non-Patent Document 5 Tomomi Ishige, Ryoichi Saito, Hiroshi Kodera: ⁇ High-precision color management by polar coordinate division of L * a * b * space ", Proceedings of the IEICE General Conference, P.188 (1998) Disclosure of the invention
  • Patent Document 1 divides the color space into three parts corresponding to predetermined primary colors C (cyan), M (magenta), and Y (yellow).
  • C cyan
  • M magenta
  • Y yellow
  • the above-mentioned problem of gradation discontinuity cannot always be effectively solved.
  • Patent Document 2 describes the point of overlapping partial areas! However, it does not specifically describe how much the partial areas overlap, and the above-described gradation characteristics are not described. It is not always possible to effectively resolve the discontinuity problem. As described above, in the conventional method, the color conversion method for an output device such as an ink jet printer has not been sufficiently studied to realize high-accuracy color conversion.
  • the present invention has been made to solve the above problems, and at a boundary portion of a partial region. It is an object of the present invention to provide a color conversion matrix creation method and a color conversion method for creating a color conversion matrix capable of suppressing discontinuity of gradation and performing highly accurate color conversion.
  • the color conversion matrix creation method includes a dividing step of dividing the first color space into a plurality of partial areas, and adjacent partial areas overlap each other at a predetermined ratio.
  • a step for setting an enlarged partial area in which each partial area is enlarged, and enlargement of the colorimetric data obtained by measuring the color of the test pattern on which color patches of a plurality of colors in the color gamut of the predetermined output device are printed Creating a color conversion matrix for each partial area, based on the colorimetric data corresponding to the partial area, for converting the image data of the second color space into the image data of the second color space. It is characterized by that.
  • the color conversion matrix creation method converts image data in a first color space that does not depend on a device into image data in a second color space corresponding to a predetermined output device.
  • a division method including the number of divisions is set for each of the three axes constituting the first color space, and the set division method is created. Dividing the first color space into a plurality of partial areas, and (2) for each of the plurality of partial areas, so that the adjacent partial areas overlap each other at a predetermined ratio.
  • Color measurement obtained by measuring the color test pattern Create a color conversion matrix for each partial area to convert the image data of the first color space to the image data of the second color space based on the colorimetric data included in the enlarged partial area of the data And a step.
  • a division method including the number of divisions is set for each of the three axes constituting the color space, and the first color space is divided into a plurality of partial regions. It is divided. Thereby, the first color space can be preferably divided into a plurality of partial regions.
  • the division method can be set only by the number of divisions. Alternatively, if necessary, additional conditions such as a division width may be set in addition to the number of divisions.
  • an enlarged partial region is set by enlarging each partial region so that adjacent partial regions overlap each other at a predetermined ratio with respect to the plurality of partial regions obtained by dividing the first color space.
  • a partial area is enlarged at a predetermined enlargement ratio, and an overlapping area overlapping with an adjacent partial area is set.
  • the enlarged partial area is an area obtained by combining the overlapping area and the original partial area.
  • the predetermined ratio can be set such that the gradation does not become discontinuous at the boundary portion of the partial area of the image after color conversion.
  • a test pattern for a plurality of colors in the color gamut of the output device determined as a conversion target in advance (for example, a test pattern in which color patches of a plurality of colors in the color gamut of the output device are printed) is measured.
  • the color conversion matrix for converting the image data of the first color space into the image data of the second color space based on the colorimetric data corresponding to the enlarged partial area in the measured color measurement data for each partial area create.
  • the color conversion matrix creating method in the setting step, for each of the three axes constituting the first color space, only a predetermined ratio r to the length I of the partial area in the axial direction is provided. It is preferable that the enlarged partial area is set by enlarging the area in both directions of the axis and that the predetermined ratio r satisfies 0.1.l ⁇ r ⁇ 0.3. As described above, by setting the enlarged partial area at the ratio r of the value selected from the range of 0.1 to 0.3, it is possible to suitably create a color conversion matrix using the enlarged partial area.
  • the predetermined ratio for enlarging the partial area to the enlarged partial area may be set separately for each direction in which the partial area is enlarged. Thereby, the expansion range of the area can be suitably set for each direction in which the partial area is expanded.
  • the first color space is a CIELAB space
  • the second color space is an RGB space.
  • the colorimetric data is the correspondence between RGB data and Lab data when each color of the test pattern is measured (for example, RGB data and Lab data when each color patch is measured) It is assumed that the data represents a correspondence relationship between
  • the division number D is set to 3 or more for the L axis of the three axes constituting the first color space. In a range
  • the division step it is preferable to set the division number D to an even number for the a axis of the three axes constituting the first color space. Also, split steps
  • the division number D is 3 or more.
  • a color conversion method is a color conversion method using a color conversion matrix for each partial area created by the above-described color conversion matrix creation method, and is an input in the input first color space.
  • the input image data is input so as to approach an achromatic color at a predetermined conversion ratio.
  • the process of converting the image data and executing the conversion step again may be repeated until the image data converted by the conversion step falls within a predetermined range. As a result, color conversion can be appropriately performed even when image data outside the color gamut is input.
  • a target point is set in the first color space, and the input image data is set to the target point by a predetermined moving method. It is preferable to convert the converted image data into image data that falls within a predetermined range by moving in the direction. Thus, even when image data outside the color gamut is input, color conversion can be performed appropriately.
  • 0 0 is an integer greater than or equal to 2
  • the process of moving the input image data so that the input image data approaches the target point at the set movement rate lZx 0 and executing the conversion step again with the force is the image after conversion by the conversion step. It is preferable to repeat until the data falls within a predetermined range. At this time, it is particularly preferable to set the movement ratio lZx as 10 ⁇ 50.
  • the target point for the input image data when there is no color conversion matrix corresponding to the input image data, or the image data after conversion to the image data of the second color space is within a predetermined range. In any case, it is preferable to set the achromatic point in the first color space as the target point.
  • FIG. 1 is a schematic diagram showing a color management system.
  • FIG. 2 is a schematic block diagram of a color conversion device.
  • FIG. 3 is a flowchart of a color conversion matrix creation process executed by a color conversion matrix creation unit.
  • FIG. 4 is a diagram showing division of Lab space.
  • FIG. 5 is a diagram showing an enlargement of a partial region.
  • FIG. 6 is a diagram showing the color gamut of the printer and the chromaticity of the SHIPP color chart.
  • FIG. 7 is a diagram showing an image (a) after color conversion by the conventional color conversion method and an image (b) after color conversion by the color conversion method of the first embodiment.
  • FIG. 8 is a diagram showing images (a) to (d) after color conversion by the color conversion method of the first embodiment.
  • FIG. 9 is a diagram showing an average color difference, a maximum color difference, and a standard deviation when the number of divisions on the L axis is changed.
  • FIG. 10 is a diagram showing an average color difference, a maximum color difference, and a standard deviation when the number of a-axis divisions is changed.
  • FIG. 11 is a diagram showing an average color difference, a maximum color difference, and a standard deviation when the number of divisions on the b-axis is changed.
  • FIG. 12 is a flowchart of color conversion processing according to the second embodiment.
  • FIG. 13 is a schematic diagram showing an example of a color conversion method when there is no color conversion matrix corresponding to input image data.
  • FIG. 14 is a schematic diagram showing an example of a color conversion method when the image data after being converted into image data in the second color space is out of the color gamut.
  • FIG. 15 shows an image (a) after color conversion by the color conversion method of the first embodiment when data outside the color gamut is input, and the second implementation when data outside the color gamut is input. It is a figure which shows the image (b) after color conversion by the color conversion method of form.
  • a color management system as a mechanism for transmitting accurate color information from an input device such as a camera to an output device such as a printer or display.
  • Fig. 1 is a schematic diagram schematically showing the color management system.
  • PCS Profile Connection Space
  • the color signal from the input device is converted into color information in PCS independent of the device by conversion C1. Further, by converting this color information by conversion C2, a color signal for the output device is generated.
  • CIELAB space, CI EXYZ space, etc. are used as such PCS.
  • the device-independent PCS color information (first color space image data) is converted into a color signal (second color space) corresponding to the output device.
  • the latter conversion C2 to be converted into image data) is provided, and a color conversion matrix creation method used for this conversion C2 and a color conversion method using a color conversion matrix obtained thereby are provided.
  • a color conversion matrix creation method and a color conversion method according to the present invention will be described using specific embodiments and examples.
  • the CIELAB color space is a color system that uses L * for lightness and a * and b * for color, and is widely used as a PCS.
  • L *, a *, and b * in the L * a * b * color space are defined by cubes in which L * is 0 to 100, a * and b * are in the range of -127 to 127, respectively.
  • the three axes (L * axis, a * axis, and b * axis) that make up such an L * a * b * color space can be A color space with a normal range of 0 to 1 shall be used.
  • a color space in which the gradation values given in 0 to 255 are similarly normalized to 0 to 1 is used.
  • FIG. 2 shows a schematic block diagram of the color conversion apparatus 10 to which the present invention is applied.
  • the color conversion device 10 includes a color conversion matrix creation unit 12, a colorimetric data storage unit 14, and a color conversion unit 16.
  • the color conversion matrix creation unit 12 calculates the coefficients of the color conversion matrix based on the color measurement data stored in the color measurement data storage unit 14.
  • the color converter 16 Using the color conversion matrix created by the conversion matrix creation unit 12, image data in the first color space that does not depend on the input device is different from the first color space and corresponds to a predetermined output device. It is converted into image data in the second color space, which is the color space to be output, and output.
  • the color converter 16 converts the input L * a * b * image data (hereinafter referred to as Lab data) into RGB image data (hereinafter referred to as RGB data).
  • the colorimetric data storage unit 14 the correspondence between RGB data and Lab data when measuring a test pattern (color chart) on which many color patches are printed is learned data (colorimetric data). Data) in advance.
  • the color of the color patch is the color gamut (color gamut) of the output device used (not shown), that is, the output device (printer or the like) that outputs an image based on the image data color-converted by the color conversion device 10.
  • Overall power is chosen to include a variety of colors.
  • a test pattern for a plurality of colors in the color gamut of the output device that is previously determined as a conversion target can be used as the test pattern.
  • the color conversion matrix creating unit 12 divides the Lab color space into a plurality of partial areas (division step). As shown in Fig. 4, a division method including the number of divisions is set for each of the three axes (L axis, a axis, and b axis) that make up the Lab color space that is the first color space. To divide Lab color space into multiple subregions. Thus, the first color space such as the Lab color space can be preferably divided into a plurality of partial areas.
  • the division width is equally divided, and the L axis is D for the number of divisions of each axis.
  • a axis is divided into D parts
  • b axis is divided into D parts (D, D, and D are natural numbers) and divided into (DXDXD) ab Lab Lab subregions 20.
  • D, D, and D are natural numbers
  • DXDXD ab Lab Lab subregions 20.
  • the method of force division showing an example in which the number of divisions is divided into (3 ⁇ 4 ⁇ 2) is not limited to this.
  • the values of D, D and D may be determined in advance or can be set by the user
  • a value set by the user may be used.
  • the division method can be set only by the number of divisions.
  • additional conditions such as a division width may be set in addition to the number of divisions.
  • the processing power of the process is as constant as the division width, as in the example above, and the color space is divided by setting only the number of divisions for each axis. ,.
  • the color conversion matrix creation unit 12 outputs division information indicating how the Lab color space is divided, for example, information such as the division numbers D 1, D 2 and D, to the color conversion unit 16. Split like this
  • the information is used when the color conversion unit 16 performs color conversion of the image data.
  • step 102 the color conversion matrix creation unit 12 sets each partial area 20 that overlaps each partial area 20 to a predetermined ratio for each of the divided partial areas.
  • the central partial area 20E has a predetermined ratio (enlargement ratio) r with respect to the length I in the a-axis direction of the partial area 20E.
  • a predetermined ratio (enlargement ratio) r to the length I of the partial area 20E in the L-axis direction r
  • the partial region 20E has a width rl on both the left and right sides in FIG.
  • the length of the enlarged partial region 21E in the a-axis direction is (l + 2r) I.
  • the partial area 20E is enlarged with a width rl on both the top and bottom sides.
  • the enlarged partial region 21E is a region including a part of each of the partial regions 20A to 20D and 20F to 20I adjacent to the partial region 20E on the L-axis a-axis plane.
  • the length I of the partial area in the axial direction is set.
  • the predetermined ratio r in the enlargement of the area is set to a value that makes the discontinuity of gradation near the boundary between the partial areas inconspicuous in reproducing a natural image. It is preferable to set the value of the ratio r to a value that satisfies, for example, 0.1 ⁇ r ⁇ 0.3. In this way, it is possible to suitably create a color conversion matrix using an enlarged partial area by setting the enlarged partial area at a ratio of values selected from the range of 0.1 to 0.3! Can do.
  • the ratio r is not limited to such a range.
  • the value of the force r may be made different between the L axis and the a axis by increasing the same ratio r in the L axis and a axis directions. Furthermore, the value of r may be varied depending on the direction of expansion of each axis, for example, the ratio of expanding the right direction of the a axis in FIG. 5 is different from the ratio of expanding the left direction of the a axis. In other words, the value of r may be different for each enlargement direction.
  • the enlarged partial areas that overlap each other in the respective axial directions are set by enlarging the partial areas in the same manner as described above.
  • step 104 a color conversion matrix is created among the plurality of partial areas. Select one partial area, and measure the colorimetric data contained in the enlarged partial area corresponding to that partial area, that is, the Lab data contained in the enlarged partial area and the RGB data corresponding to this Lab data. Read from data storage unit 14.
  • the black spot shown in FIG. 5 is the colorimetric data
  • the selected partial area is the partial area 20E shown in FIG.
  • the enlarged partial region 21E corresponding to this is a region including a part of the partial regions 20A to 20D and 20F to 20I adjacent to the partial region 20E.
  • the color measurement data 24 included in the area overlapping with the adjacent partial area is also read simultaneously.
  • step 106 by executing a multiple regression analysis using the colorimetric data in the read enlarged partial area, Lab data (image data in the first color space) in the selected partial area is executed.
  • Each coefficient of the color conversion matrix shown in the following equation (1) for converting the image data into RGB data (image data in the second color space) is calculated (creation step).
  • the second-order term was used because the color conversion matrix with high accuracy cannot be obtained with the regression coefficient based on the first-order term, and the regression analysis is performed if the third-order or higher term is used. This is because a very large number of colorimetric data (the number of learning data) is required.
  • step 108 whether or not the color conversion matrix has been created for all of the partial areas obtained by dividing the Lab color space, that is, the coefficients for calculating the coefficients of the color conversion matrix. Judgment is made. If there is a partial area for which no color conversion matrix has been created, the process returns to step 104 and repeats the same processing as described above. If a color conversion matrix has been created for all partial areas, the process proceeds to step 110, and the coefficients of the created color conversion matrix for all partial areas are output to the color conversion unit 16 and this routine is terminated. To do.
  • the color conversion unit 16 in the color conversion information storage unit provided in the color conversion unit 16, the coefficient of the color conversion matrix of each partial area output from the color conversion matrix creation unit 12 is stored.
  • the color conversion information storage unit stores information necessary for color conversion processing such as division information indicating how the Lab color space is divided. If the color conversion matrix creation has been executed in advance, the color conversion unit 16 can read the information stored in the color conversion information storage unit and execute the color conversion.
  • the color conversion data included in the adjacent partial areas is used in part, so that the adjacent partial areas are created. Therefore, it is possible to perform highly accurate color conversion in consideration of part of the color information. Therefore, it is possible to suppress the occurrence of discontinuity in gradation near the boundary of the partial area, and accurate color reproduction is possible.
  • the inventor has verified three points of the conversion accuracy, processing time, and color reproducibility of a natural image of the present invention. To examine these, we used five SHIPP CIELAB standard images. For comparison purposes, color conversion was also performed for the ICC profile created with the profile creation software. In this embodiment, an image obtained by capturing a printed matter to which the present invention is applied with a scanner is shown.
  • the enlargement ratio r of the partial area is a value within the range of 0.1 to 0.3. Two.
  • the average color difference due to the color conversion according to the present invention is 2.32, which is smaller than the prism interpolation calculation. It is the same level as when using Photoshop, and the third class color difference (1.2 to 2.5) ( Non-patent document 8 below), and high-precision color conversion is possible.
  • the conversion processing time was measured for four natural images (3072 X 4096) of SHIPP. We also calculated the conversion processing time when prism interpolation was performed on ICC profile data. As a result, the average processing time of the present invention was 67.2% of the prism interpolation calculation, and the calculation processing time could be greatly shortened.
  • the present invention was applied to four SHIPP natural images, and the created images were printed.
  • the image (b) has no pseudo contour due to the overlap of the partial areas! / I can confirm that.
  • image (a) is a female image, ⁇ Bride ⁇ for human skin reproducibility evaluation
  • image (b) is an image with a fine geometric structure.
  • image (c) is an image containing abundant color information of high saturation, ⁇ Wool ⁇ for color reproducibility evaluation
  • image (d) is ⁇ B for metal luster and neutral color reproducibility evaluation. In all images, good images without pseudo contours were obtained, and the effectiveness of the present invention was confirmed.
  • Non-Patent Document 6 G. Sharma, W. Wu, ENDalal: "The CIEDE2000 Color- Difference r ormula: Implementation Notes, supplementary 1 est Data ⁇ and Mathematical Obs ervations, Color Research and Application ⁇ vol.30 , No. l (2005)
  • Non-Patent Document 7 Triquebs: ⁇ Practical digital color image design and evaluation ⁇ (2000)
  • Non-Patent Document 8 The Japan Society of Color Science: New Color Science Handbook [2nd edition] Tsuji, University of Tokyo Press P.290 (1998)
  • Example 1-2 five SHIPP CIELAB standard images were used as in Example 1-1.
  • EPSON PX-G900 was used as the printer. Photo paper was used for printing.
  • the PX-G900 can print highly saturated colors with a wider color gamut in the cyan, yellow, and red directions than sRGB. In particular, the cyan direction has a wider color gamut than Adobe RGB.
  • the color material is pigment ink.
  • the measurement point force using the measured values of 11,284 color charts is evenly distributed in the above printer gamut. It was made to distribute.
  • L Measures the number conversion dependency of color conversion accuracy when the other 2 axes are not divided, and averages 161 colors in the color gamut The color difference, the maximum color difference, and the standard deviation were compared, and suitable division numbers for each axis were examined.
  • FIG. 9 shows the average color difference, maximum color difference, and standard deviation when the number of divisions D on the L axis is changed.
  • the change in average color difference has become smaller in the above range. Therefore, it is preferable to set the division number D in the range of 3 or more for the L axis. In addition, when the number of divisions is greater than D force ⁇ , the change in average color difference is even smaller.
  • the standard deviation is the smallest in 7 and 8 divisions. In addition, the maximum color difference is significantly smaller at 7 divisions compared to others.
  • Fig. 10 shows the average color difference, maximum color difference, and standard deviation when the division number D of the a-axis is changed.
  • Figure 11 shows the average color difference, maximum color difference, and standard deviation when the number of divisions D on the b-axis is changed.
  • Graph (a) shows the D dependence of average color difference and maximum color difference
  • graph (b) shows the standard b.
  • the color conversion apparatus 10 described in the first embodiment uses the color conversion apparatus 10 described in the first embodiment. Can obtain a good image. However, there is a possibility that the printer cannot faithfully reproduce the partial area near the boundary of the color gamut! /, And the color is also subjected to matrix calculation and accurate color conversion cannot be performed! /. In other words, since the color conversion described in the first embodiment deals with normal RGB, the solution should be 0 or more and 1 or less, and may be negative or greater than 1. In the color conversion described in the first embodiment, there is no colorimetric data (learning data)! Since the color conversion matrix is not created for the area (outside the color gamut), color conversion is not possible. Nao.
  • color conversion processing executed by the color conversion unit 16 and corresponding to the case where data outside the color gamut is input will be described with reference to the flowchart shown in FIG.
  • the color conversion process shown in FIG. 12 is executed every time Lab data is input.
  • step 200 the input Lab data is color-converted by the same method as in the first embodiment. That is, it is determined which partial area in the input Lab data force Lab space belongs to, and using the color conversion matrix corresponding to the partial area to which the input Lab data belongs, it is converted into RGB data by the above equation (1). Convert.
  • step 202 it is determined whether or not all values of R, G, and B data after conversion are 0 or more and 1 or less. If all values are 0 or more and 1 or less, step 206 is performed. To output RGB data after conversion. On the other hand, if at least one of the R, G, and B data is 0 or more and 1 or less, the process proceeds to step 204.
  • a color conversion matrix that exists closest to the direction of the target point (for example, the achromatic color direction) is applied to the partial area where no color conversion matrix exists.
  • the image data converted into the image data of the second color space is out of the predetermined range, for example, any one of the converted RGB as in the above example is other than 0 to 1.
  • the target point is set in the first color space, and the input image data is moved in the direction of the target point by a predetermined movement method, so that the converted image data falls within the predetermined range. It is preferable to convert to data.
  • color conversion can be appropriately performed even when there is no color conversion matrix or when image data outside the color gamut is input. This also enables color conversion that does not become unnatural even when there is no color conversion matrix or when image data outside the color gamut is input.
  • the target point set for the input image data is the achromatic point (gray point) in the first color space, even in the case of! / ) Is preferably set as the target point. Further, the target point may be set by position coordinates in the color space or may be set by a target partial area.
  • this target point is a point of (0.5, 0.5, 0.5) when considered as a normal value.
  • other points for example, predetermined points in the vicinity of the achromatic point may be set as the target points.
  • FIG. 13 is a schematic diagram showing an example of a color conversion method when there is no color conversion matrix corresponding to input image data. In Fig.
  • the a-axis-b-axis plane is divided into 7 by each axis, and among the obtained 7 X 7 partial areas, the partial area 50A where the color conversion matrix a to o exists is shown.
  • ⁇ 50O is indicated by a solid line
  • the ⁇ partial area where no color conversion matrix exists outside the color gamut is indicated by a broken line.
  • the center point of the achromatic partial region 50G located at the center of the a-axis and b-axis plane is set as the target point.
  • the color conversion matrix b of the nearest partial area 50B (indicated by a thick line in the figure) where the color conversion matrix exists is assigned as the color conversion matrix of the area 50P, that is, the color conversion matrix for the input image data. .
  • the color conversion process with respect to the input image data outside the color gamut can be suitably executed.
  • all partial areas including areas outside the color gamut are displayed. .
  • Color conversion matrix is set
  • FIG. 14 is a schematic diagram showing an example of a color conversion method when the image data after being converted into image data in the second color space is out of a predetermined range (out of the color gamut).
  • input image data outside the color gamut is moved in the direction of the target point (moving target point) by a predetermined moving method.
  • the input image data and the movement target point are divided into X (X is an integer of 2 or more). And input image data and target point
  • the process of moving the input image data outside the color gamut into the color gamut is suitably executed. It becomes possible.
  • 0 and X is the number of divisions between the input image data and the movement target point.
  • the initial value of ⁇ corresponding to the input image data is 0, and 1 is added to X for each movement process.
  • the inventor of the present invention applies the color according to the present embodiment to an XYZ image of a wide color gamut photographed by a camera (see Non-Patent Document 9 below) that can acquire accurate XYZ tristimulus values of a subject having a wide color gamut. A conversion method was applied.
  • the image (a) in Fig. 15 shows the case where all the R, G, B data after color conversion is not between 0 and 1, or there is no learning data in the partial area corresponding to the input Lab data. And an image printed with black. As shown in this image (a), colors outside the color gamut cannot be reproduced.
  • an image (b) in FIG. 15 shows an image printed by applying the color conversion method according to the present embodiment. In this image (b), the black part of the image (a) is also converted into a natural color, and a good image can be obtained even for a wide color gamut image.
  • Non-Patent Document 9 Tomohiro Horiuchi, Ijietztersen, Gosuke Ohashi, Mifumi Shimohira: “Development of an imaging device for accurate color information acquisition", 11th Image Sensing Symposium, pp5 41-542 (2005)
  • Patent Document 10 Caro Fujiya, "Corresponding Color Reproduction from Softcopy Images to Hardcopy Images", Chiba University Doctoral Letter, January 2002
  • the present invention relates to a color conversion matrix creation method and a color conversion method for creating a color conversion matrix capable of suppressing color discontinuity at a boundary portion of a partial region and capable of highly accurate color conversion. Is available.

Landscapes

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

Abstract

　Ｌａｂ色空間を複数の部分領域に分割し（ステップ１００）、隣接する部分領域が所定割合で相互に重複するように各部分領域を拡大した拡大部分領域を設定し（ステップ１０２）、出力装置の色域内の複数の色についてのテストパターンを測色した測色データのうち拡大部分領域に応じた測色データに基づいて、第１の色空間の画像データを第２の色空間の画像データに変換するための色変換マトリクスを部分領域毎に作成する（ステップ１０４～１０８）。これにより、高精度な色変換が可能な色変換マトリクスを作成する色変換マトリクス作成方法及び色変換方法が実現される。

Description

色変換マトリクス作成方法及び色変換方法
技術分野
[0001] 本発明は、色変換マトリクス作成方法及び色変換方法に係り、特に、行列演算によ り色変換を行う際の色変換マトリクス作成方法、及び当該色変換マトリクス作成方法 により作成した色変換マトリクスを用いて入力画像を色変換する色変換方法に関する 背景技術
[0002] 近年、プラズマテレビや液晶ディスプレイに代表される表示装置にぉ ヽて、高解像 度、広色域の表示装置が発売されている。また、カラー印刷においてもインクジェット プリンタの普及により、忠実な色再現に関する人々の関心が高まっている（例えば非 特許文献 1参照)。インクジェットプリンタは、機構が比較的単純で小型化が可能、プ リンタ本体が低価格である、などの特徴を有し、インターネットで検索した画像や、デ ジタルカメラで撮影した画像、プレゼン資料などの印刷に多く利用されて 、る。
[0003] プリンタの色変換は、 R (赤）、 G (緑)、 B (青）の加法混色に基づいて行うことが可能 なディスプレイの色変換とは異なり、減法混色の原理や紙、インクなどの特性が非線 形であるため、理論式を立てるのは困難である。これに対し、入出力機器の色変換と して、機器の色情報を記述した ICCプロファイル（ICC :InternationalColor Conso rtium)を用いた手法が標準化されて 、る（例えば非特許文献 2参照)。
[0004] プリンタに対して適用される ICCプロファイルは、複数の入出力関係を対応づけた L UT(Look-uptable)によるものであり、 LUTとテーブル補間法を組み合わせて色 変換を行う（例えば非特許文献 3参照)。し力しながら、画像に対する変換では 1画素 ずつ補間演算を行うため、高解像度の画像では非常に多くの変換時間を必要とする
[0005] また、インクジェットプリンタにおける高精度な色変換手法として、色空間を複数の 部分領域に分割し、その部分領域ごと行列演算を行う手法が知られている (例えば 非特許文献 4、 5参照)。し力しながら、これらの手法では、自然画像の再現において 階調性の不連続が生じ、擬似輪郭が見られるなどの問題が発生する。
[0006] また、特許文献 1及び特許文献 2には、各部分領域を相互に重複させて色変換テ 一ブルや色変換マトリクスを作成する技術が開示されている。
特許文献 1 :特開 2002— 112053号公報
特許文献 2 :特開 2000— 188695号公報
特許文献 3：特開 2004— 266590号公報
非特許文献 1：田代耕二： "インクジェットプリンタで思 、通りの色を出すには"、 日本写 真学会誌、 vol.68, No.l、 pp.70- 72 (2005)
特干文献 2： International し olor Consortium： iし C Profile Specincation Version4.2 .0.0〃 (2004)
非特許文献 3 : Dawn Wallner： "Building ICC profiles- the Mechanics and Engineering 〃 (2000)
非特許文献 4:小寺宏曄、石毛淳美、斎藤了一:〃色空間の分割による入出力デバイ スの高精度カラーマッチング"、画像電子学会誌、 vol.29, pp.128-139 (2000) 非特許文献 5 :石毛淳美、斎藤了一、小寺宏曄:〃 L*a*b*空間の極座標分割による高 精度カラーマネージメント"、電子情報通信学会総合大会講演論文集、 P.188 (1998) 発明の開示
発明が解決しょうとする課題
[0007] し力しながら、上記特許文献 1に記載された技術は、予め定めた一次色 C (シアン） 、 M (マゼンダ）、 Y (イェロー）に対応して色空間を 3分割するものであり、分割方法 や分割数の点で、上記の階調性の不連続の問題を効果的に解消できるとは限らな い。
[0008] また、上記特許文献 2には、部分領域を重複させる点は記載されて!、るものの、ど の程度重複させるのか等については具体的に記載されておらず、上記の階調性の 不連続の問題を効果的に解消できるとは限らない。以上のように、従来の方法では、 インクジェットプリンタ等の出力装置に対する色変換手法について、高精度の色変換 を実現するための検討が充分になされて 、るとは言えな 、。
[0009] 本発明は上記課題を解決するために成されたものであり、部分領域の境界部分で 階調が不連続になるのを抑制し高精度な色変換が可能な色変換マトリクスを作成す る色変換マトリクス作成方法及び色変換方法を提供することを目的とする。
課題を解決するための手段
[0010] 上記目的を達成するため、本発明による色変換マトリクス作成方法は、第 1の色空 間を複数の部分領域に分割する分割ステップと、隣接する部分領域が所定割合で 相互に重複するように各部分領域を拡大した拡大部分領域を設定する設定ステップ と、予め定めた出力装置の色域内の複数の色のカラーパッチが印刷されたテストバタ 一ンを測色した測色データのうち拡大部分領域に応じた測色データに基づいて、第 ェの色空間の画像データを第 2の色空間の画像データに変換するための色変換マト リクスを部分領域毎に作成する作成ステップとを含むことを特徴とする。
[0011] より具体的には、色変換マトリクス作成方法は、デバイスに依存しない第 1の色空間 の画像データを、予め定めた出力装置に対応する第 2の色空間の画像データに変 換するための色変換マトリクスを作成する作成方法であって、（1)第 1の色空間を構 成する 3軸のそれぞれにつ ヽて分割数を含む分割方法を設定し、設定された分割方 法によって第 1の色空間を複数の部分領域に分割する分割ステップと、（2)複数の部 分領域のそれぞれにつ!ヽて、隣接する部分領域が所定割合で相互に重複するよう に、第 1の色空間を構成する 3軸のそれぞれの方向に部分領域を拡大した拡大部分 領域を設定する設定ステップと、（3)複数の部分領域のそれぞれに対し、上記出力 装置の色域内の複数の色についてのテストパターンを測色して得られた測色データ のうち拡大部分領域に含まれる測色データに基づいて、第 1の色空間の画像データ を第 2の色空間の画像データに変換するための色変換マトリクスを部分領域毎に作 成する作成ステップとを含むことを特徴とする。
[0012] この発明によれば、第 1の色空間に対し、色空間を構成する 3軸のそれぞれについ て分割数を含む分割方法を設定して、第 1の色空間を複数の部分領域へと分割して いる。これにより、第 1の色空間を好適に複数の部分領域へと分割することができる。 分割方法の設定については、分割数のみによって分割方法を設定することができる 。あるいは、必要に応じて、分割数に加えて分割幅などの追加の条件を設定しても良 い。 [0013] また、第 1の色空間を分割した複数の部分領域に対し、隣接する部分領域が所定 割合で相互に重複するように各部分領域を拡大した拡大部分領域を設定する。すな わち、部分領域を所定の拡大率で拡大させ、隣接する部分領域と重複する重複領 域を設定する。拡大部分領域は、この重複領域と元の部分領域とを併せた領域であ る。なお、所定割合は、色変換後の画像の部分領域の境界部分において階調が不 連続にならないような割合に設定することができる。
[0014] そして、予め変換対象として定められた出力装置の色域内の複数の色についての テストパターン (例えば、出力装置の色域内の複数の色のカラーパッチが印刷された テストパターン)を測色した測色データのうち拡大部分領域に応じた測色データに基 づいて、第 1の色空間の画像データを第 2の色空間の画像データに変換するための 色変換マトリクスを部分領域毎に作成する。
[0015] これにより、部分領域毎の色変換マトリクスを作成する際に、部分領域内の測色デ ータだけでなぐ隣接する部分領域と重複する領域の測色データも考慮されて色変 換マトリクスが作成されるため、部分領域同士の境界部分において階調が不連続に なるのを抑制することができる。なお、予め定められる出力装置については、上記方 法は、出力装置をプリンタとして適用することができる。このようなプリンタとしては、例 えば、インクジェットプリンタがある。
[0016] また、色変換マトリクス作成方法は、設定ステップにおいて、第 1の色空間を構成す る 3軸のそれぞれについて、その軸方向の部分領域の長さ Iに対して所定割合 r分だ け軸の両方向に領域を拡大することによって拡大部分領域を設定するとともに、所定 割合 rが 0. l≤r≤0. 3を満たすことが好ましい。このように、 0. 1〜0. 3の範囲から 選択された値の割合 rで拡大部分領域を設定することにより、拡大部分領域を用いた 色変換マトリクスの作成を好適に行うことができる。
[0017] また、設定ステップにお!/、て、部分領域を拡大部分領域へと拡大するための上記 所定割合は、部分領域を拡大する方向毎に別個に設定可能であるものとしても良い 。これにより、部分領域を拡大するそれぞれの方向について、領域の拡大範囲を好 適に設定することができる。
[0018] また、第 1の色空間は、 CIELAB空間であると共に、第 2の色空間は、 RGB空間で あり、測色データは、テストパターンにおける複数の色のそれぞれを測色した時の RG Bデータと Labデータとの対応関係（例えば、各カラーパッチを測色した時の RGBデ ータと Labデータとの対応関係）を表すデータであるものとすることができる。
[0019] また、第 1の色空間の複数の部分領域への分割については、分割ステップにおい て、第 1の色空間を構成する 3軸のうちの L軸について、分割数 Dを 3以上の範囲で
し
設定することが好ましい。また、分割ステップにおいて、第 1の色空間を構成する 3軸 のうちの a軸について、分割数 Dを偶数に設定することが好ましい。また、分割ステツ
a
プにおいて、第 1の色空間を構成する 3軸のうちの b軸について、分割数 Dを 3以上
b
の範囲で設定することが好ま 、。
[0020] 本発明による色変換方法は、上記の色変換マトリクス作成方法により作成された部 分領域毎の色変換マトリクスを用いた色変換方法であって、入力された第 1の色空間 における入力画像データが属する部分領域を判定する判定ステップと、上記色変換 マトリクス作成方法により作成された部分領域毎の色変換マトリクスのうち、入力画像 データが属する部分領域に対応した色変換マトリクスにより、入力画像データを、出 力装置に対応する第 2の色空間の画像データに変換する変換ステップとを含むこと を特徴とする。
[0021] このように、部分領域内の測色データだけでなぐ隣接する部分領域と重複する領 域の測色データも考慮されて作成された色変換マトリクスにより入力画像データを変 換することにより、変換後の画像の部分領域の境界部分において階調が不連続にな るのを抑 ff¾することができる。
[0022] また、色変換方法にお!、て、第 2の色空間の画像データに変換後の画像データが 所定範囲外の場合、入力画像データが所定変換割合で無彩色に近づくように入力 画像データを変換して力も再度変換ステップを実行する処理を、変換ステップによる 変換後の画像データが所定範囲内になるまで繰り返すようにしてもよい。これにより、 色域外の画像データが入力された場合でも適切に色変換することができる。
[0023] また、色変換方法において、入力画像データに対応する色変換マトリクスが存在し ない場合、第 1の色空間において目標点を設定し、目標点の方向に向力つて進み、 色変換マトリクスが存在する最も近 、部分領域の色変換マトリクスを、入力画像デー タに割り当てることが好ましい。これにより、色変換マトリクスが存在しない場合であつ ても、適切に色変換することができる。
[0024] あるいは、第 2の色空間の画像データに変換後の画像データが所定範囲外の場合 、第 1の色空間において目標点を設定し、入力画像データを所定の移動方法で目標 点の方向に移動することによって、変換後の画像データが所定範囲内となる画像デ ータに変換することが好ましい。これにより、色域外の画像データが入力された場合 であっても、適切に色変換することができる。
[0025] この場合、具体的な画像データの移動方法につ!ヽては、入力画像データと目標点 との間の色空間上の距離に対して 1回当たりの移動割合を lZx (X
0 0は 2以上の整数
)と設定し、設定された移動割合 lZx 0で入力画像データが目標点に近づくように入 力画像データを移動して力も再度変換ステップを実行する処理を、変換ステップによ る変換後の画像データが所定範囲内になるまで繰り返すことが好ましい。また、このと き、移動割合 lZxを 10≤χ≤ 50として設定することが特に好ましい。
0 0
[0026] また、入力画像データに対する目標点の設定については、入力画像データに対応 する色変換マトリクスが存在しない場合、あるいは、第 2の色空間の画像データに変 換後の画像データが所定範囲外の場合のいずれにおいても、第 1の色空間におけ る無彩色の点を、目標点に設定することが好ましい。
発明の効果
[0027] 本発明によれば、部分領域の境界部分で階調が不連続になるのを抑制し高精度 な色変換が可能な色変換マトリクスを作成する、 、う効果がある。
図面の簡単な説明
[0028] [図 1]図 1は、カラーマネージメントシステムを示す模式図である。
[図 2]図 2は、色変換装置の概略ブロック図である。
[図 3]図 3は、色変換マトリクス作成部で実行される色変換マトリクス作成処理のフロー チャートである。
[図 4]図 4は、 Lab空間の分割について示す図である。
[図 5]図 5は、部分領域の拡大について示す図である。
[図 6]図 6は、プリンタの色域と SHIPPカラーチャートの色度を示す図である。 [図 7]図 7は、従来の色変換方法による色変換後の画像 (a)、及び第 1実施形態の色 変換方法による色変換後の画像 (b)を示す図である。
[図 8]図 8は、第 1実施形態の色変換方法による色変換後の画像 (a)〜 (d)を示す図 である。
[図 9]図 9は、 L軸の分割数を変化させたときの平均色差、最大色差、標準偏差を示 す図である。
[図 10]図 10は、 a軸の分割数を変化させたときの平均色差、最大色差、標準偏差を 示す図である。
[図 11]図 11は、 b軸の分割数を変化させたときの平均色差、最大色差、標準偏差を 示す図である。
[図 12]図 12は、第 2実施形態に係る色変換処理のフローチャートである。
[図 13]図 13は、入力画像データに対応する色変換マトリクスが存在しな!ヽ場合の色 変換方法の一例を示す模式図である。
[図 14]図 14は、第 2の色空間の画像データに変換後の画像データが色域外の場合 の色変換方法の一例を示す模式図である。
[図 15]図 15は、色域外のデータが入力された場合における第 1実施形態の色変換 方法による色変換後の画像 (a)、及び色域外のデータが入力された場合における第 2実施形態の色変換方法による色変換後の画像 (b)を示す図である。
符号の説明
[0029] 10· · ·色変換装置、 12· · ·色変換マトリクス作成部、 14· · ·測色データ記憶部、 16· · · 色変換部、 20· · ·部分領域、 21 · · ·拡大部分領域、 50· · ·部分領域。
発明を実施するための最良の形態
[0030] 以下、本発明の好ましい実施の形態について図面を参照しながら詳細に説明する
[0031] まず、本発明の概要について説明する。カメラなどの入力装置から、プリンタやディ スプレイなどの出力装置へと正確な色情報を伝達するための仕組みとしてカラーマ ネージメントシステムがある。図 1は、カラーマネージメントシステムを概略的に示す模 式図である。 [0032] デバイス間での色情報の伝達には、図 1に示すように、デバイスに依存しない色情 報である PCS (Profile Connection Space)が用いられている。すなわち、図 1に 示した例では、入力装置からの色信号は、変換 C1によってデバイスに依存しない P CSでの色情報に変換される。さらに、この色情報を変換 C2で変換することによって、 出力装置に対する色信号が生成される。このような PCSとしては、 CIELAB空間、 CI EXYZ空間などが用いられる。
[0033] 本発明は、これらの変換 Cl、 C2のうちで、デバイスに依存しない PCSの色情報（ 第 1の色空間の画像データ)を、出力装置に対応する色信号 (第 2の色空間の画像 データ）に変換する後者の変換 C2を対象とし、この変換 C2に用いられる色変換マト リクス作成方法、及びそれによつて得られる色変換マトリクスを用いた色変換方法を 提供するものである。以下、本発明による色変換マトリクス作成方法、及び色変換方 法について、具体的な実施形態及び実施例を用いて説明する。
[0034] なお、以下の説明においては、上記した第 1の色空間から第 2の色空間への色変 換について、 CIELAB空間力 RGB空間への色変換を中心に説明する。 CIELAB 色空間は、明度を表す L*と、色合いを表す a*、 b*とによる表色系であり、 PCSとして 広く用いられている。
[0035] また、 L*a*b*色空間における L*、 a*、 b*についてそれぞれ、 L*が 0〜100、 a*、 b*が— 127〜127の範囲の立方体で定義する。また、このような L*a*b*色空間（第 ェの色空間）を構成する 3軸 (L*軸、 a*軸、 b*軸）について、必要に応じて、上記した 各軸の範囲を 0〜1に正規ィ匕した色空間を用いることとする。また、これに合わせて、 RGB色空間についても、必要に応じて、 0〜255で与えられている階調値を同様に 0 〜1に正規ィ匕した色空間を用いることとする。
[0036] (第 1実施形態）
図 2には、本発明が適用された色変換装置 10の概略ブロック図を示した。図 2に示 すように、色変換装置 10は、色変換マトリクス作成部 12、測色データ記憶部 14、及 び色変換部 16を含んで構成されて 、る。
[0037] 詳細は後述するが、色変換マトリクス作成部 12は、測色データ記憶部 14に記憶さ れた測色データに基づいて色変換マトリクスの係数を算出する。色変換部 16は、色 変換マトリクス作成部 12で作成された色変換マトリクスを用いて、入力されたデバイス に依存しない第 1の色空間における画像データを、第 1の色空間とは異なり、予め定 めた出力装置に対応する色空間である第 2の色空間における画像データに変換して 出力する。
[0038] 本実施形態では、上記したように、第 1の色空間を CIELAB表色系の色空間とし、 また、第 2の色空間を RGB表色系の色空間とした場合について説明する。すなわち 、色変換部 16は、入力された L*a*b*の画像データ（以下、 Labデータという）を、 R GBの画像データ（以下、 RGBデータという）に変換する。
[0039] 測色データ記憶部 14には、多数の色のカラーパッチが印刷されたテストパターン（ 色票）を測色した時の RGBデータと Labデータとの対応関係が学習データ（測色デ ータ）として予め記憶されている。なお、カラーパッチの色は、使用する出力装置（図 示省略)、すなわち色変換装置 10により色変換された画像データに基づいて画像を 出力する出力装置 (プリンタ等)の色域 (色再現域)全体力 まんべんなく様々な色を 含むように選択される。また、テストパターンとしては、一般には、予め変換対象として 定められた出力装置の色域内の複数の色についてのテストパターンを用いることが できる。
[0040] 次に、色変換マトリクス作成部 12において実行される色変換マトリクス作成処理に ついて、図 3に示すフローチャートを参照して説明する。
[0041] まずステップ 100では、色変換マトリクス作成部 12は、 Lab色空間を複数の部分領 域に分割する（分割ステップ)。図 4に示すように、第 1の色空間である Lab色空間を 構成する 3軸 (L軸、 a軸、 b軸）のそれぞれについて分割数を含む分割方法を設定し 、設定された分割方法によって Lab色空間を複数の部分領域に分割する。これにより 、 Lab色空間などの第 1の色空間を、好適に複数の部分領域へと分割することができ る。
[0042] 図 4に示す例では、分割幅については等分とし、各軸の分割数について、 L軸を D
し 等分、 a軸を D等分、 b軸を D等分して (D 、 D、 Dは自然数）、（D X D X D )個の a b L a b L a b 部分領域 20に分割する。具体的には、図 4では、分割数の設定について（3 X 4 X 2 )個に分割した例を示している力 分割の仕方はこれに限られるものではない。なお、 D 、 D 、 Dの値は、予め定めておいてもよいし、ユーザーによって設定可能な構成
L a b
として、ユーザーが設定した値を用いてもよい。
[0043] ここで、第 1の色空間の各軸に対する分割方法の設定については、分割数のみに よって分割方法を設定することができる。あるいは、必要に応じて、分割数に加えて 分割幅などの追加の条件を設定しても良い。ただし、処理の簡単ィ匕の観点力もは、 上記した例のように分割幅等につ!、ては一定とし、分割数のみを各軸で設定して色 空間の分割を行うことが好まし 、。
[0044] また、色変換マトリクス作成部 12は、 Lab色空間をどのように分割したかを示す分割 情報、例えば分割数 D 、 D 、 D等の情報を色変換部 16に出力する。このような分割
L a b
情報は、色変換部 16において画像データの色変換を行う際に用いられる。
[0045] 次に、ステップ 102において、色変換マトリクス作成部 12は、分割された複数の部 分領域のそれぞれについて、各部分領域 20を相互に重複させるベぐ各部分領域 2 0を所定の割合で Labの各軸の方向（第 1の色空間を構成する 3軸のそれぞれの方 向)へ拡大させた部分領域 (以下、拡大部分領域と!/、う）を設定する (設定ステップ)。
[0046] 例えば図 5に示すように、 L軸 a軸平面が 3 X 3に分割されて、部分領域 20A〜2 01に分割された場合について説明する。この場合、中央の部分領域 20Eについては 、部分領域 20Eの a軸方向の長さ Iに対して所定割合 (拡大率) r分 a軸方向両側に a
拡大させるとともに、部分領域 20Eの L軸方向の長さ Iに対して所定割合 (拡大率) r
し
分 L軸方向両側に拡大させ、この領域を拡大部分領域 21Eとして設定する。
[0047] このとき、 a軸方向については、図 5中の左右両側において幅 rlで部分領域 20Eが a
拡大されている。これにより、拡大部分領域 21Eの a軸方向の長さは（l + 2r) Iとなつ a ている。また、 L軸方向については、上下両側において幅 rlで部分領域 20Eが拡大
し
されている。これにより、拡大部分領域 21Eの L軸方向の長さは（l + 2r) Iとなってい
し
る。以上により、拡大部分領域 21Eは、 L軸 a軸平面において、部分領域 20Eに隣 接する部分領域 20A〜20D、 20F〜20Iのそれぞれの一部を含む領域となる。
[0048] 一般には、この拡大部分領域の設定については、上記した例のように、第 1の色空 間を構成する 3軸のそれぞれについて、その軸方向の部分領域の長さ Iに対して所 定割合 r分だけ軸の両方向に領域を拡大することによって拡大部分領域を設定する ことが好ましい。
[0049] また、この領域の拡大における所定割合 rは、自然画像の再現において部分領域 間の境界付近の階調性の不連続が目立たなくなる値に設定される。このような割合 r の値としては、例えば 0. l≤r≤0. 3を満たす範囲の値に設定することが好ましい。こ のように、 0. 1〜0. 3の範囲から選択された値の割合!:で拡大部分領域を設定するこ とにより、拡大部分領域を用いた色変換マトリクスの作成を好適に行うことができる。 ただし、割合 rの値については、このような範囲の値に限られるものではない。
[0050] また、本実施形態では、 L軸及び a軸方向にそれぞれ同じ割合 rずつ拡大させて 、 る力 rの値を L軸、 a軸のそれぞれで異ならせてもよい。さらに、例えば図 5において a 軸の右方向に拡大する割合を a軸の左方向に拡大する割合と異ならせる等、各軸の 拡大方向に応じて rの値を異ならせても良い。すなわち、拡大する方向毎に rの値を 異ならせても良い。
[0051] そして、 L軸 b軸平面、 a軸 b軸平面でも上記と同様に部分領域を拡大させるこ とにより、各軸方向に相互に重複する拡大部分領域が設定される。
[0052] ステップ 104では、複数の部分領域のうちで色変換マトリクスが作成されて 、な！/、部 分領域を一つ選択し、その部分領域に対応した拡大部分領域に含まれる測色デー タ、すなわち、拡大部分領域に含まれる Labデータとこの Labデータに対応する RGB データを測色データ記憶部 14から読み込む。
[0053] 例えば図 5に示す黒点を測色データとし、選択された部分領域が図 5に示す部分 領域 20Eである場合について考える。このとき、これに対応する拡大部分領域 21E は部分領域 20Eに隣接する部分領域 20A〜20D、 20F〜20Iの一部を含む領域と なっている。このため、部分領域 20E内の測色データ 22に加えて、隣接する部分領 域と重複した領域に含まれる測色データ 24も同時に読み込まれる。
[0054] ステップ 106では、読み込んだ拡大部分領域内の測色データを用いて重回帰分析 を実行することにより、選択された部分領域内の Labデータ (第 1の色空間の画像デ ータ）を RGBデータ (第 2の色空間の画像データ）に変換するための次式（1)で示す 色変換マトリクスの各係数を算出する（作成ステップ)。
[0055] すなわち、この例では、次式（1)で示すように、 L*a*b*の 2次の項を用いて、色変 換マトリクスの係数 A 〜A を重回帰分析により求める。
11 39
[数 1]
Figure imgf000014_0001
ここで、 2次の項を用いたのは、 1次の項による回帰係数では精度の良い色変換マト リクスを得ることができず、また 3次以上の項を用いたのでは回帰分析を行うのに非常 に多くの測色データ数 (学習データ数)を必要とするためである。
[0056] ステップ 108では、 Lab色空間を分割した複数の部分領域に対し、全ての部分領 域にっ 、て色変換マトリクスを作成したか否力、すなわち色変換マトリクスの各係数を 算出したカゝ否かを判断する。そして、色変換マトリクスを作成していない部分領域が 存在する場合には、ステップ 104へ戻って上記と同様の処理を繰り返す。また、全て の部分領域にっ ヽて色変換マトリクスを作成した場合には、ステップ 110へ移行し、 作成した全部分領域の色変換マトリクスの係数を色変換部 16へ出力して本ルーチン を終了する。
[0057] 色変換部 16では、色変換部 16に設けられた色変換情報記憶部において、色変換 マトリクス作成部 12から出力された各部分領域の色変換マトリクスの係数を記憶して おく。また、この色変換情報記憶部には、 Lab色空間をどのように分割したかを示す 分割情報などの色変換処理に必要な情報が格納される。色変換マトリクス作成が予 め実行済みの場合には、色変換部 16は、色変換情報記憶部に記憶された情報を読 み出して、色変換を実行することができる。
[0058] そして、画像データの変換時には、入力された Labデータが Lab空間内のどの部分 領域に属するかを判定する (判定ステップ)。そして、入力された Labデータが属する 部分領域に対応する色変換マトリクスを用いて、上記（1)式により、プリンタなどの指 定された出力装置に対応する RGBデータに変換する (変換ステップ)。
[0059] このように、本実施形態では、各部分領域の色変換マトリクスを作成する際に、隣接 する部分領域に含まれる測色データを一部利用して作成するため、隣接する部分領 域の色情報も一部考慮した高精度な色変換が可能となる。従って、部分領域の境界 付近にお 、て階調性の不連続が発生するのを抑えることができ、正確な色再現が可 能となる。
[0060] (実施例 1 1)
以下、第 1実施形態の実施例について説明する。なお、本発明は下記実施例によ り限定されるものではない。
[0061] 本発明者は、本発明の変換精度、処理時間、及び自然画像に対する色再現性の 3 点について検証した。これらを調べるために、 SHIPPの CIELAB標準画像 5枚を使 用した。今回は比較のために、プロファイル作成ソフトにより作成した ICCプロフアイ ルでも色変換を行った。なお、本実施例では、本発明を適用した印刷物をスキャナで 取り込んだ画像を示す。
[0062] 検証には、広色域の色が再現可能な 8色顔料インクタイプのプリンタを用いた。また 、印刷には写真用紙を使用した。この時の色域を図 6に示す。
[0063] 色変換マトリクスを作成するために、学習データ（測色データ）は 11, 281色の色票 の測定値を用いた。また、測定器は GretagMacbethSpectroScanを用いた。
[0064] Lab色空間の分割数は D =6、 D =4、 D = 11に設定し、得られた 6 X 4 X I I = 2
L a b
64個の部分領域に対応して、 264個の色変換マトリクスを作成した。また、拡大部分 領域への領域の拡大については、部分領域の拡大率 rは、 0. 1〜0. 3の範囲内の 値である。. 2とした。
[0065] SHIPPのカラーチャートにおいて、使用したプリンタにおける色域内の 161色に対 する入力値との差を CIE2000色差式（下記非特許文献 6参照）で求めた。 161色の 色度を図 6に、結果を以下の表 1に示す。また、表 1には比較として、 Adobe Photosho p 6.0で ICCプロファイルを用いたときの色差の他に、 ICCプロファイルデータをプリズ ム補間処理 (下記非特許文献 7参照)したときの結果も示す。 [表 1]
Figure imgf000016_0001
表 1より、本発明に係る色変換による平均色差は 2. 32となり、プリズム補間演算よりも 色差は小さぐ Photoshopを用いたときと同程度で 3級色差（1. 2〜2. 5) (下記非特 許文献 8参照）となり、高精度な色変換が可能となった。
[0066] また、 SHIPPの 4枚の自然画像（3072 X 4096)に対して変換処理時間を測定し た。また、 ICCプロファイルデータをプリズム補間処理したときの変換処理時間も求め た。その結果、本発明の平均処理時間は、プリズム補間演算の 67. 2%となり、演算 処理時間を大幅に短縮することができた。
[0067] また、本発明を SHIPPの自然画像 4枚に適用し、作成した画像を印刷した。部分 領域の重複による影響を確認するために、 Woolの画像を部分領域の拡大率 r=0で 変換した時の画像を図 7の画像 (a)に、本発明のように部分領域の拡大率 r=0. 2で 変換した時の画像を図 7の画像 (b)に示す。この図 7の画像 (a)、画像 (b)の楕円内 を見比べればわ力るように、画像 (b)では、部分領域の重複によって擬似輪郭がなく なって！/、ることが確認できる。
[0068] その他の結果を図 8の画像 (a)〜（d)に示す。ここで、画像 (a)は女性の画像で、人 肌の再現性評価用の〃Bride〃、画像 (b)は細かな幾何学的構造を持つ画像で、画像 処理手法の評価用の "Harb0r"、画像 (c)は高彩度の色情報を豊富に含む画像で、 色再現性評価用の〃Wool〃、画像 (d)は金属光沢感と中性色の再現性評価用の〃B ottles"である。全ての画像において擬似輪郭もなく良好な画像が得られ、本発明の 有効性を確認できた。
[0069] (非特許文献 6) G.Sharma、 W.Wu、 E.N.Dalal : "The CIEDE2000 Color- Difference r ormula : Implementation Notes、 supplementary 1 est Data^ and Mathematical Obs ervations 、 Color Research and Application^ vol.30、 No. l (2005)
(非特許文献 7) トリケッブス：〃実践ディジタルカラー画像の設計と評価〃 (2000) (非特許文献 8) 日本色彩学会：〃新編色彩科学ハンドブック [第 2版]〃、東京大学 出版会 P.290 (1998)
[0070] (実施例 1 2)
第 1実施形態の第 2の実施例について説明する。
[0071] 本実施例 1—2では、上記実施例 1—1と同様に、 SHIPPの CIELAB標準画像 5枚 を使用した。また、検証には、プリンタとして EPSON PX— G900を用いた。また、 印刷には写真用紙を使用した。 PX—G900は、 sRGBよりもシアン、イェロー、レッド 方向で色域が広ぐ彩度の高い色を印刷することができる。特に、シアン方向に対し ては、 Adobe RGBよりも広い色域を持っている。また、色材は顔料インクである。
[0072] 色変換マトリクスを作成するために用いられる学習データ (測色データ）については 、 11, 284色の色票の測定値を用いた測定点力 上記プリンタの色域にまんべんな く分布するようにした。
[0073] このような条件において、第 1の色空間である Lab色空間を複数の部分領域に分割 する際の分割数 D 、 D 、 Dの好適な値について検討した。具体的には、 L軸、 a軸、
L a b
b軸の 3軸について、 1つの軸のみを 1〜： L 1分割し、他の 2軸を分割しないときの色変 換精度の分割数依存性を測定し、色域内の色 161色の平均色差、最大色差、標準 偏差を比較して、それぞれの軸に対して好適な分割数を検討した。
[0074] 図 9は、 L軸の分割数 Dを変化させたときの平均色差、最大色差、標準偏差を示す
し
図であり、グラフ (a)は平均色差及び最大色差の D依存性を示し、グラフ (b)は標準
し
偏差の D依存性を示している。これらのグラフより、 L軸については、分割数 Dが 3
し し 以上の範囲で平均色差の変化が小さくなつている。したがって、 L軸については、分 割数 Dを 3以上の範囲で設定することが好ましい。さらに、分割数 D力 ^以上になる し し と平均色差の変化はさらに小さくなる。また、標準偏差は 7、 8分割で最も小さい。また 、最大色差については、 7分割のときに他と比較して著しく小さくなつている。
[0075] 図 10は、 a軸の分割数 Dを変化させたときの平均色差、最大色差、標準偏差を示 a
す図であり、グラフ (a)は平均色差及び最大色差の D依存性を示し、グラフ (b)は標 a
準偏差の D依存性を示している。これらのグラフより、 a軸については、平均色差、最 a
大色差、標準偏差のすべてが、偶数個で分割する方が奇数個で分割するよりも値が 小さくなつている。したがって、 a軸については、分割数 Dを偶数に設定することが好 a
ましい。これは、偶数個で分割する場合、 a軸の中心である a* =0 (正規ィ匕した値で考 えると a* =0. 5)を通る面で領域が分割されているためだと考えられる。また、分割数 Dを大きくすると、色差等の値は小さくなる。
a
[0076] 図 11は、 b軸の分割数 Dを変化させたときの平均色差、最大色差、標準偏差を示 b
す図であり、グラフ (a)は平均色差及び最大色差の D依存性を示し、グラフ (b)は標 b
準偏差の D依存性を示している。これらのグラフより、 b軸については、分割数 Dが 3 b b 以上の範囲で平均色差の変化が小さくなつている。したがって、 b軸については、 L 軸と同様に、分割数 Dを 3以上の範囲で設定することが好ましい。また、標準偏差に b
つ!、ては、 11分割のときに最も小さ 、値となって!/、る。
[0077] 以上の図 9〜図 11の結果に基づいて、暫定的な分割数を D = 7、 D = 10、 D = 1
L a b
1に設定した。そして、 L軸、 a軸、 b軸の 3軸について、 1つの軸のみで分割数を変化 させ、他の 2軸を上記の暫定分割数としたときの色変換精度の分割数依存性を測定 した。このような測定の結果により、好適な分割数の設定の一例として、分割数を D
し
= 3、 D = 10、 D = 11に設定した。
a b
[0078] このような分割条件において、表 1に示した場合と同様に、本発明を適用した場合 の色差、 ICC + Photoshopでの色差、及び ICC +プリズム補間での色差を求めた結 果を表 2に示す。また、ここでは、参考としてプリンタドライバによる色補正を行ったと きの精度も示している。
[表 2]
Figure imgf000018_0001
表 2より、本発明を適用した場合の色変換精度は、すべてにおいてプリズム補間より も値が小さくなつた。また、 Photoshopと比較すると、最大色差は本発明の方が小さく 、平均色差についても Photoshopとほぼ同程度となった。 [0079] (第 2実施形態）
次に、本発明の第 2実施形態について説明する。
[0080] 上記の第 1実施形態の実施例 1 1で示したように、 SHIPP自然画像はすべて使 用したプリンタの色域内の色であるため、第 1実施形態で説明した色変換装置 10で は良好な画像を得ることができる。し力しながら、色域の境界付近の部分領域では， プリンタが忠実に再現できな!/、色も行列演算し、正確な色変換ができな!/、可能性が ある。つまり、上記第 1実施形態で説明した色変換では正規ィ匕した RGBを扱っている ために解は 0以上 1以下となるはず力 負または 1より大きくなることがある。また、上 記第 1実施形態で説明した色変換では、測色データ（学習データ）が存在しな!、領 域 (色域外）に対しては色変換マトリクスが作成されな ヽため色変換できな ヽ。
[0081] そこで、本実施形態では、色変換部 16で実行される、色域外のデータが入力され た場合に対応した色変換処理について、図 12に示すフローチャートを参照して説明 する。図 12に示す色変換処理は、 Labデータが入力される毎に実行される。
[0082] まず、ステップ 200では、入力された Labデータを第 1実施形態と同様の方法で色 変換する。すなわち、入力された Labデータ力Lab空間内のどの部分領域に属する かを判定し、入力された Labデータが属する部分領域に対応する色変換マトリクスを 用いて、上記（1)式により RGBデータに変換する。
[0083] ステップ 202では、変換後の R, G, Bの各データ全ての値が 0以上 1以下であるか 否かを判断し、全ての値が 0以上 1以下である場合にはステップ 206へ移行して変換 後の RGBデータを出力する。一方、 R, G, Bの各データのうち一つでも 0以上 1以下 でな 、値のデータが存在する場合は、ステップ 204へ移行する。
[0084] ステップ 204では、入力された Labデータを無彩色（正規ィ匕した値で考えて L*， a* , b* = 0. 5)の方向に所定割合 (例えば 20分の 1)だけ近づけた値に変換し、ステツ プ 200へ戻って再度色変換を行う。
[0085] このように、色変換後の R, G, Bの各データが全て 0以上 1以下であればその値を 出力し、そうでなければ R, G, Bの各データが全て 0以上 1以下になるまでステップ 2 00〜204の処理を繰り返す。これにより、出力装置の色域外のデータが入力された 場合でも適切に色変換することができる。 [0086] 本実施形態における色変換方法について、さらに説明する。上記したように、第 1 実施形態の色変換方法においても、色域の境界付近の色、及び色域外の色につい ては、高精度に色変換することができない場合がある。これに対して、本実施形態で は、（1)色変換マトリクスが存在しない部分領域には、目標点の方向（例えば、無彩 色方向）に向力つて最も近くに存在する色変換マトリクスを割り当てる、及び（2)部分 領域に対応する色変換マトリクスにより得られた RGBのうち、どれか 1つでも 0〜1以 外であったときに色域外処理を行う、 t 、う処理を行う。
[0087] すなわち、入力画像データに対応する色変換マトリクスが存在しない場合、第 1の 色空間において目標点を設定し、目標点の方向に向かって進み、色変換マトリクス が存在する最も近い部分領域の色変換マトリクスを、入力画像データに割り当てるこ とが好ましい。
[0088] また、第 2の色空間の画像データに変換後の画像データが所定範囲外、例えば上 記した例のように変換後の RGBのうちでどれか 1つでも 0〜1以外であった場合、第 1 の色空間にお!、て目標点を設定し、入力画像データを所定の移動方法で目標点の 方向に移動することによって、変換後の画像データが所定範囲内となる画像データ に変換することが好ましい。
[0089] このような方法を用いることにより、色変換マトリクスが存在しない場合、あるいは色 域外の画像データが入力された場合であっても、適切に色変換することができる。ま た、これにより、色変換マトリクスが存在しない場合、あるいは色域外の画像データが 入力された場合であっても、不自然とならない色変換が可能となる。
[0090] 入力画像データに対して設定される目標点については、図 12に関して上述したよ うに、上記した!/、ずれの場合でも、第 1の色空間での無彩色の点（グレーの点）を目 標点に設定することが好ましい。また、この目標点は、色空間での位置座標で設定し ても良ぐあるいは目標とする部分領域で設定しても良い。
[0091] 上記した Lab色空間の場合、目標点として好適に用いられる無彩色の点は、（L* , a*, b*) = (50, 0, 0)の点となる。また、この目標点は、正規ィ匕した値で考えると (0. 5, 0. 5, 0. 5)の点となる。また、必要に応じてこれ以外の点、例えば無彩色の点の 近傍にある所定の点、を目標点としても良い。 [0092] 図 13は、入力画像データに対応する色変換マトリクスが存在しない場合の色変換 方法の一例を示す模式図である。図 13 (a)においては、 a軸—b軸平面を各軸で 7分 割したとし、得られた 7 X 7の部分領域のうちで、色変換マトリクス a〜oが存在する部 分領域 50A〜50Oを実線で示し、色域外で色変換マトリクスが存在しな ヽ部分領域 を破線で示している。また、ここでは、 a軸 b軸平面の中心に位置する無彩色の部 分領域 50Gの中心点を目標点とする。
[0093] このような入力色空間において、入力画像データが図中の左上にある色域外の領 域 50Pにあったとする。この場合、領域 50Pには色変換マトリクスが存在しないので、 図 13 (a)に矢印で示すように、入力画像データから、グレーに対応する部分領域 50 Gの目標点に向力つて進む方向を想定する。
[0094] そして、色変換マトリクスが存在する最も近い部分領域 50B (図中、太線で示す)の 色変換マトリクス bを、領域 50Pの色変換マトリクス、すなわち入力画像データに対す る色変換マトリクスとして割り当てる。これにより、色域外の入力画像データに対する 色変換処理を、好適に実行することができる。さらに、このような割り当て処理を色変 換マトリクスが存在しないすベての領域に対して行うことにより、図 13 (b)に示すよう に、色域外の領域を含むすべての部分領域に対して、色変換マトリクスが設定される
[0095] 図 14は、第 2の色空間の画像データに変換後の画像データが所定範囲外 (色域 外)の場合の色変換方法の一例を示す模式図である。このような場合、色域外の入 力画像データを目標点 (移動目標点）の方向に所定の移動方法で移動する処理を 行う。
[0096] 画像データの移動方法については、図 14に示す例では、入力画像データと移動 目標点との間を X分割 (Xは 2以上の整数)する。そして、入力画像データと目標点と
0 0
の間の色空間上の距離に対して i回当たりの移動割合を iZxと設定するとともに、
0
設定された移動割合 lZx
0で入力画像データが移動目標点に近づくように入力画像 データを移動して力も再度変換ステップを実行する処理を、変換ステップによる変換 後の画像データが色域内になるまで繰り返す方法を用いて 、る。
[0097] これにより、色域外の入力画像データの色域内への移動処理を、好適に実行する ことが可能となる。また、上記方法での移動割合 ΐΖχ
0については、移動割合 ΐΖχ
0 を 10≤x ≤ 50として設定することが好ましぐ特に、図 12に関して上述したように、 X
0 0
= 20とすることが好ましい。
[0098] ここで、正規化した L*， a* , b*の値を Cとし、設定された移動目標点での値を Cと
0 し、入力画像データと移動目標点との間の分割数を Xとする。このとき、上記した移
0
動処理を X回実行した後の画像データは、下記の式（2)で表される。
[数 2]
Figure imgf000022_0001
ここで、入力画像データに対応する χの初期値は 0で、移動処理毎に Xに 1を足してい
<ο
[0099] また、無彩色の点を移動目標点 (C =0. 5)とし、分割数を X = 20と設定すると、上
0 0
記の式（2)は下記の式 (3)のようになる。
[数 3]
C = C -— (C - 0.5) C ^ L a ,b* ' · · (3)
20 入力画像データが色域外の場合、上記の式によってデータを移動して再度、色変換 を行う。そして、 RGBがすべて 0〜1の値であれば得られた値を出力し、そうでなけれ ば Xに 1を足してさらにデータを移動する。このような処理を行うことにより、色域外の 入力画像データを色域内の画像データへと好適に変換することができる。
[0100] (実施例 2)
以下、第 2実施形態の実施例について説明する。なお、本発明は下記実施例によ り限定されるものではない。
[0101] 本発明者は、広 ヽ色域を持つ被写体の正確な XYZ三刺激値が取得できるカメラ（ 下記非特許文献 9参照)で撮影した広色域 XYZ画像に対して本実施形態による色 変換方法を適用した。
[0102] 図 15の画像 (a)は、色変換後の R, G, B各データ全てが 0以上 1以下でない場合、 若しくは入力された Labデータに対応する部分領域に学習データが存在しない場合 に、黒色を印刷した画像である。この画像 (a)に示すように、色域外の色は再現でき ていない。一方、図 15の画像 (b)は、本実施形態に係る色変換方法を適用して印刷 した画像を示す。この画像 (b)では、画像 (a)の黒の部分も自然な色に変換されてお り、広色域画像に対しても良好な画像を得ることができた。
[0103] (非特許文献 9) 堀内智博、イジエッツターセン、大橋剛介、下平美文: "正確な色 情報の取得を目指した撮像装置の開発"、第 11回画像センシングシンポジウム、 pp5 41-542 (2005)
( 特許文献 10) カロ藤 哉、 "Corresponding Color Reproduction from Softcopy I mages to Hardcopy Images",千葉大学博士餘文、 2002年 1月
産業上の利用可能性
[0104] 本発明は、部分領域の境界部分で階調が不連続になるのを抑制し、高精度な色変 換が可能な色変換マトリクスを作成する色変換マトリクス作成方法及び色変換方法と して利用可能である。

Claims

請求の範囲
[1] デバイスに依存しない第 1の色空間の画像データを、予め定めた出力装置に対応 する第 2の色空間の画像データに変換するための色変換マトリクスを作成する作成 方法であって、
前記第 1の色空間を構成する 3軸のそれぞれについて分割数を含む分割方法を設 定し、設定された前記分割方法によって前記第 1の色空間を複数の部分領域に分割 する分割ステップと、
前記複数の部分領域のそれぞれにつ ヽて、隣接する部分領域が所定割合で相互 に重複するように、前記第 1の色空間を構成する 3軸のそれぞれの方向に前記部分 領域を拡大した拡大部分領域を設定する設定ステップと、
前記複数の部分領域のそれぞれに対し、前記出力装置の色域内の複数の色につ いてのテストパターンを測色して得られた測色データのうちで前記拡大部分領域に 含まれる測色データに基づいて、前記第 1の色空間の画像データを前記第 2の色空 間の画像データに変換するための前記色変換マトリクスを前記部分領域毎に作成す る作成ステップと
を含むことを特徴とする色変換マトリクス作成方法。
[2] 前記設定ステップにおいて、前記第 1の色空間を構成する 3軸のそれぞれについ て、その軸方向の前記部分領域の長さ Iに対して所定割合 r分だけ軸の両方向に領 域を拡大することによって前記拡大部分領域を設定するとともに、前記所定割合 rが 0. l≤r≤0. 3を満たすことを特徴とする請求項 1記載の色変換マトリクス作成方法。
[3] 前記設定ステップにおいて、前記部分領域を前記拡大部分領域へと拡大するため の前記所定割合は、前記部分領域を拡大する方向毎に別個に設定可能であること を特徴とする請求項 1または 2記載の色変換マトリクス作成方法。
[4] 前記第 1の色空間は、 CIELAB空間であると共に、前記第 2の色空間は、 RGB空 間であり、
前記測色データは、前記テストパターンにおける前記複数の色のそれぞれを測色 したときの RGBデータと Labデータとの対応関係を表すデータであることを特徴とす る請求項 1〜3のいずれか一項記載の色変換マトリクス作成方法。
[5] 前記分割ステップにおいて、前記第 1の色空間を構成する 3軸のうちの L軸につい て、前記分割数 Dを 3以上の範囲で設定することを特徴とする請求項 4記載の色変
し
換マトリクス作成方法。
[6] 前記分割ステップにおいて、前記第 1の色空間を構成する 3軸のうちの a軸につい て、前記分割数 Dを偶数に設定することを特徴とする請求項 4または 5記載の色変
a
換マトリクス作成方法。
[7] 前記分割ステップにおいて、前記第 1の色空間を構成する 3軸のうちの b軸につい て、前記分割数 Dを 3以上の範囲で設定することを特徴とする請求項 4〜6のいずれ
b
か一項記載の色変換マトリクス作成方法。
[8] 請求項 1〜7のいずれか一項記載の色変換マトリクス作成方法により作成された前 記部分領域毎の前記色変換マトリクスを用いた色変換方法であって、
入力された前記第 1の色空間における入力画像データが属する前記部分領域を判 定する判定ステップと、
前記色変換マトリクス作成方法により作成された前記部分領域毎の前記色変換マト リクスのうち、前記入力画像データが属する部分領域に対応した色変換マトリクスによ り、前記入力画像データを、前記出力装置に対応する前記第 2の色空間の画像デー タに変換する変換ステップと
を含むことを特徴とする色変換方法。
[9] 前記入力画像データに対応する前記色変換マトリクスが存在しない場合、前記第 1 の色空間において目標点を設定し、前記目標点の方向に向力つて進み、色変換マト リクスが存在する最も近!、部分領域の色変換マトリクスを、前記入力画像データに割 り当てることを特徴とする請求項 8記載の色変換方法。
[10] 前記第 2の色空間の画像データに変換後の画像データが所定範囲外の場合、前 記第 1の色空間において目標点を設定し、前記入力画像データを所定の移動方法 で前記目標点の方向に移動することによって、前記変換後の画像データが前記所定 範囲内となる画像データに変換することを特徴とする請求項 8または 9記載の色変換 方法。
[11] 前記入力画像データと前記目標点との間の色空間上の距離に対して 1回当たりの 移動割合を lZx (χ
0 0は 2以上の整数)と設定し、設定された前記移動割合 lZx
0で 前記入力画像データが前記目標点に近づくように前記入力画像データを移動してか ら再度前記変換ステップを実行する処理を、前記変換ステップによる前記変換後の 画像データが前記所定範囲内になるまで繰り返すことを特徴とする請求項 10記載の 色変換方法。
[12] 前記移動割合 lZxを 10≤x≤ 50として設定することを特徴とする請求項 11記載
0 0
の色変換方法。
[13] 前記第 1の色空間における無彩色の点を、前記目標点に設定することを特徴とする 請求項 9〜 12のいずれか一項記載の色変換方法。
PCT/JP2006/323741 2005-11-28 2006-11-28 色変換マトリクス作成方法及び色変換方法 WO2007061117A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007546535A JP4806778B2 (ja) 2005-11-28 2006-11-28 色変換マトリクス作成方法を用いた色変換方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2005342855 2005-11-28
JP2005-342855 2005-11-28

Publications (1)

Publication Number Publication Date
WO2007061117A1 true WO2007061117A1 (ja) 2007-05-31

Family

ID=38067337

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2006/323741 WO2007061117A1 (ja) 2005-11-28 2006-11-28 色変換マトリクス作成方法及び色変換方法

Country Status (2)

Country Link
JP (1) JP4806778B2 (ja)
WO (1) WO2007061117A1 (ja)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09130626A (ja) * 1995-11-01 1997-05-16 Ricoh Co Ltd 色信号変換方法
JP2000188695A (ja) * 1998-12-22 2000-07-04 Toshiba Corp 色変換対応テーブル構築方法および色変換装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09130626A (ja) * 1995-11-01 1997-05-16 Ricoh Co Ltd 色信号変換方法
JP2000188695A (ja) * 1998-12-22 2000-07-04 Toshiba Corp 色変換対応テーブル構築方法および色変換装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KODERA K. ET AL: "IRo Kukan no Bunkatsu ni yoru Nyushutsuryoku Device no Koseido Color Muching", JOURNAL OF THE INSTITUTE OF IMAGE ELECTRONICS ENGINEERS OF JAPAN, vol. 29, no. 2, 25 March 2000 (2000-03-25), pages 128 - 139, XP003013630 *

Also Published As

Publication number Publication date
JPWO2007061117A1 (ja) 2009-05-07
JP4806778B2 (ja) 2011-11-02

Similar Documents

Publication Publication Date Title
US7206100B2 (en) Image processing method and apparatus
US6897988B1 (en) Image processing apparatus, image processing method, and storage medium for color matching
JP3403157B2 (ja) 画像処理装置およびその方法
US7312891B2 (en) Image processing method and apparatus
US7154634B2 (en) Calibration method, calibration apparatus, and memory medium storing the calibration method of printer
JP2006197457A (ja) カラーチャート作成方法、カラーチャート作成装置、およびカラーチャート作成プログラム、ならびに色変換定義修正方法
JP3990859B2 (ja) 色処理方法およびその装置
EP1441505A2 (en) A smooth gray component replacement strategy that utilizes the full device gamut
JP6275058B2 (ja) カラープロファイル作成方法、及びカラープロファイル作成プログラム
US7229146B2 (en) System and method for characterizing a printing device
US7515300B2 (en) Image processing apparatus and image processing method
JP2008042711A (ja) 色変換装置、色変換方法、色変換プログラム、色変換係数作成装置、色変換係数作成方法、及び色変換係数作成プログラム
JPH10164381A (ja) 色変換方法
JP2007312313A (ja) 画像処理装置、画像処理方法及びプログラム
JP3749102B2 (ja) 画像処理装置、画像形成装置、並びに画像処理方法
JP4806778B2 (ja) 色変換マトリクス作成方法を用いた色変換方法
JP3776326B2 (ja) 色変換テーブル作成方法
JP4466037B2 (ja) 画像処理装置及び画像処理方法
US20200031136A1 (en) Color prediction method
US20100253960A1 (en) Sensitivity matrix determination via chain rule of differentiation
JPH1169176A (ja) 色再現方法
JPH10164352A (ja) 画像信号処理方法
JP4603665B2 (ja) 色処理装置およびその方法
JP2005210225A (ja) 画像処理装置及び画像処理方法
JP3292519B2 (ja) 画像処理装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application
DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2007546535

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 06833545

Country of ref document: EP

Kind code of ref document: A1

DPE1 Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101)