Color matching method and system based on osterwald full color gamut mixing model

Abstract

The application discloses a color matching method and system based on an Ostwald full-color gamut color mixing model, belongs to the field of color matching theory, color model construction, and color design methods and application development technology in the textile field, and discloses the following technical scheme: taking actually obtained red, yellow, green, blue, the highest brightness gray and the lowest brightness gray as the basis, an Ostwald color solid that can be applied to actual industrial design is constructed, and the traditional 24 color phases are broken through, and the Ostwald color solid can be flexibly applied to different color design scenes; in addition, in the grid model of the Ostwald color solid, the grid point position corresponding to each color is directly related to the mixing concentration ratio of black, white and pure color pigments required for matching the color, so that a designer can directly obtain the proportion of each pigment from the model to obtain a target color when matching colors, and therefore, the design efficiency is greatly improved.

Description

TECHNICAL FIELD

[0001] The present application relates to a color matching method and system based on Ostwald full color gamut color mixing model, belonging to the field of color matching theory, color model construction and color design method and application development technology in the field of textiles. BACKGROUND

[0002] Ostwald color solid, also known as Ostwald's color solid, is a three-dimensional biconical color solid composed of 24 isochromatic surfaces, 8 equal whiteness surfaces and 8 equal blackness surfaces, and is characterized by hue, whiteness and blackness.

[0003] Ostwald color solid is based on standard four primary colors to construct a hue ring containing 24 hues, and is based on standard white and black to construct an eight-grade gray axis from standard white to standard black. The pure and turbid changes and light and dark changes of various hues are processed with the same step to construct 24 isochromatic surfaces. Therefore, Ostwald color solid explains the quantitative relationship between color and mixing ratio with a color distribution rule. However, when used in industrial design, Ostwald color solid has the following problems:

[0004] First, the color values of the red, yellow, green, blue, black and white primary colors produced by industrial production are not the standard colors defined by color optics. The color values of the four primary colors are not actually uniformly distributed on the hue ring, and the black and white colors are actually high-lightness gray close to white and low-lightness gray close to black. Therefore, when designers perform color matching, if the actual six primary colors are used as the basis and the mixing concentrations of each color are calculated according to the color distribution rule of the standard Ostwald color solid, the ideal color matching effect cannot be obtained, resulting in a large color matching error.

[0005] Second, the existing Ostwald color solid only provides color cards, and even if the color distribution rule is disclosed, designers cannot directly obtain the mixing ratio of each color when performing color matching. In combination with the error caused by the actual color values of the six primary colors, the amount of each color can only be adjusted manually by referring to the color card, and the color matching efficiency is very low. Therefore, the existing Ostwald color solid cannot be applied to efficient and accurate industrial color design. SUMMARY

[0006] In order to provide an Ostwald color solid that can be applied to actual industrial color design, reduce color difference and improve color matching efficiency, the present application provides a color matching method and system based on Ostwald full color gamut color mixing model, which comprises:

[0007] The first object of the present application is to provide a color matching method, which comprises:

[0008] Step one: based on the actual color measurement, the color values of the six primary colors of the Ostwald color solid are obtained, and an Ostwald full-color domain gridding color mixing model is constructed;

[0009] Step two: the target color is obtained from the Ostwald full-color domain gridding color mixing model, and the mixing concentration values of the six primary colors required for adjusting the color are obtained according to the grid point coordinates of the target color in the model;

[0010] Step three: the six primary colors are mixed based on the mixing concentration obtained in step two, and the target color material is obtained;

[0011] The construction process of the Ostwald full-color domain gridding color mixing model includes:

[0012] Step 1: taking the color values of red, yellow, green and blue color materials as the reference, an Ostwald color solid color ring is constructed, including:

[0013] Four color mixing intervals are constructed, and binary coupling color mixing is performed based on the color values of the two end points of the color mixing interval to obtain the color value C(ξ) of each hue on the Ostwald color solid color ring, which is expressed as:

[0014]

[0015] The color mixing concentration of each hue on the Ostwald color solid color ring is:

[0016]

[0017] Wherein, ξ represents the grid point serial number of each hue on the Ostwald color solid color ring, ε m ≥2; δ=1, 2, 3, 4 respectively represent the serial numbers of the four color mixing intervals; [r α (δ),g α (δ),b α (δ)] and [r β (δ),g β (ξ),b β (δ)] respectively represent the color values of the two end points of the δth color mixing interval; And respectively represent the mixing concentration of the color mixing interval end point color C α (δ) and C β (δ) corresponding to the ξth grid point, and the increasing gradient of the mixing concentration is 1 / ε m ;

[0018] The grid point coordinate value of each hue on the Ostwald color solid color ring is:

[0019]

[0020] wherein θ ξ , ρ ξ and z ξ represent the polar angle, polar radius and height of the grid point, respectively;

[0021] Step 2: Construct the gray axis of the Ostwald color solid with the color values of white and black as the color values of the highest and lowest lightness nodes, and the polar coordinates of the 8 nodes on the gray axis are:

[0022]

[0023] The color values corresponding to the 8 nodes on the gray axis and the mixing concentrations are obtained by equal ratio increment, equal difference increment or non-linear increment;

[0024] The concentration value matrix of the 8 nodes on the gray axis is:

[0025]

[0026] The color value matrix of the 8 nodes on the gray axis is:

[0027]

[0028] wherein η = 1 corresponds to the gray axis constructed by common ratio coefficient increment, η = 2 corresponds to the gray axis constructed by equal difference coefficient increment, and η = 3 corresponds to the gray axis constructed by non-linear increment;

[0029] Step 3: Construct each isochromatic plane of the Ostwald color solid with a certain hue color value C ξ = (r ξ , g ξ , b ξ ), the highest lightness node color value C W = (r W , g W , b W ) and the lowest lightness node color value C K = (r K , g K , b K ), and the three-primary color mixing concentrations of the color blocks on each isochromatic plane are:

[0030]

[0031] The color value of the color block corresponding to the grid point P(τ, δ, ξ) is:

[0032]

[0033] ​Where δ = 1, 2, ..., 7, 8 are the iso-blackness indices, τ = 1, 2, ..., 7, 8 are the iso-whiteness indices, γ = 1, 2, ..., 7, 8 are the iso-purity indices, and Ω... ξ ,Ω W ,Ω K These represent the weights of the pigments for the three primary colors.

[0034] Optionally, step 2 constructs a grayscale axis using a method that increases brightness proportionally, and the color mixing density value matrix of the eight nodes of the grayscale axis is as follows:

[0035]

[0036] in,

[0037]

[0038] The color values ​​of each color node on the grayscale axis are:

[0039]

[0040] Among them, C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ] represent the color values ​​of the highest and lowest brightness nodes, respectively.

[0041] Optionally, step 2 uses an arithmetic progression of brightness to construct the Ostwald grayscale axis, and the color mixing density value matrix of the eight nodes of the grayscale axis is as follows:

[0042]

[0043] The color values ​​of each color node on the grayscale axis are:

[0044]

[0045] Among them, C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ] represent the color values ​​of the highest and lowest brightness nodes, respectively.

[0046] Optionally, step 2 uses a non-linear brightness increment method to construct the grayscale axis, and the color mixing density value matrix of the eight nodes of the grayscale axis is as follows:

[0047]

[0048] The color values of each color node of the gray scale axis are:

[0049]

[0050] where C W = [r W ,g W ,b W ] and C K = [r K ,g K ,b K ] represent the color values of the highest and lowest lightness nodes, respectively.

[0051] Optionally, when ξ = 1, 2,..., 23, 24, the equal white planes of the Ostwald solid color are respectively:

[0052] The color matrix of the equal white plane when τ = 1 is:

[0053]

[0054] The color matrix of the equal white plane when τ = 2 is:

[0055]

[0056] The color matrix of the equal white plane when τ = 3 is:

[0057]

[0058] The color matrix of the equal white plane when τ = 4 is:

[0059]

[0060] The color matrix of the equal white plane when τ = 5 is:

[0061]

[0062] The color matrix of the equal white plane when τ = 6 is:

[0063]

[0064] The color matrix of the equal white plane when τ = 7 is:

[0065]

[0066] The color matrix of the equal white plane when τ = 8 is:

[0067]

[0068] Optionally, when ξ = 1, 2,..., 23, 24, the isochromatic surfaces of the Ostwald solid of colors are respectively:

[0069]

[0070] The color matrix of the isochromatic surface when δ = 2 is:

[0071]

[0072] The color matrix of the isochromatic surface when δ = 3 is:

[0073]

[0074] The color matrix of the isochromatic surface when δ = 4 is:

[0075]

[0076] The color matrix of the isochromatic surface when δ = 5 is:

[0077]

[0078] The color matrix of the isochromatic surface when δ = 6 is:

[0079]

[0080] The color matrix of the isochromatic surface when δ = 7 is:

[0081]

[0082] The color matrix of the isochromatic surface when δ = 8 is:

[0083]

[0084] Optionally, when ξ = 1, 2,..., 23, 24, the isopurity surfaces of the Ostwald solid of colors are respectively:

[0085]

[0086] The color matrix of the isopurity surface when γ = 2 is:

[0087]

[0088] The color matrix of the isopurity surface when γ = 3 is:

[0089]

[0090] The color matrix of the equal purity surface when γ=4 is:

[0091]

[0092] The color matrix of the equal purity surface when γ=5 is:

[0093]

[0094] The color matrix of the equal purity surface when γ=6 is:

[0095]

[0096] The color matrix of the equal purity surface when γ=7 is:

[0097]

[0098] The color matrix of the equal purity surface when γ=8 is:

[0099]

[0100] A second object of the present application is to provide a color matching system for implementing the color matching method according to any one of the preceding items, the system comprising a color measurement device, an Ostwald full color gamut gridding color mixing model construction module and a visualization module;

[0101] The color measurement device is configured to obtain the color values of the Ostwald color solid six primary colors in the actual colorant;

[0102] The Ostwald full color gamut gridding color mixing model construction module comprises:

[0103] A colorant color value acquisition module configured to obtain the color values of the actual colorant from the color measurement device;

[0104] A hue ring construction module configured to construct an equatorial hue ring based on the actually obtained Ostwald six primary colors, and give the color values, position coordinate values and mixing concentrations of each hue;

[0105] A gray axis construction module configured to construct a gray axis based on the actually obtained black and white as the lowest lightness node and the highest lightness node respectively, and give the color values, position coordinate values and mixing concentrations of black and white of each lightness node;

[0106] An isochromatic surface color spectrum construction module configured to construct 4ε m isochromatic surface color spectrums of the Ostwald color solid;

[0107] The visualization module is configured to display the model color spectrum, and output the color mixing concentrations and color values corresponding to the six primary colors based on the grid point coordinate values of the target color.

[0108] Optionally, the Ostwald full color gamut gridding color mixing model construction module further comprises:

[0109] An isotri-chromaticity surface color spectrum construction module is configured to construct, based on the 4ε m equal-chromaticity surface color spectrums, a conical surface formed by all grid points with equal white content on each chromaticity surface as an isotri-chromaticity surface of the Ostwald color solid;

[0110] An isotetra-chromaticity surface color spectrum construction module is configured to construct, based on the 4ε m equal-chromaticity surface color spectrums, a conical surface formed by all grid points with equal black content on each chromaticity surface as an isotetra-chromaticity surface of the Ostwald color solid;

[0111] An isopurity surface color spectrum construction module is configured to construct, based on the 4ε m equal-chromaticity surface color spectrums, a cylindrical surface formed by all grid points with equal color content on each chromaticity surface as an isopurity surface of the Ostwald color solid.

[0112] A third object of the present application is to provide a computer-readable storage medium, the storage medium having stored thereon a computer program which, when executed by a processor, implements the method according to any one of the preceding aspects.

[0113] The present application has the following advantages:

[0114] 1. The obtained red, yellow, green and blue are used as the primary colors of the four color pigments, the position of the grid point is defined by the mixing concentration of the four primary color pigments, the color value of each grid point is obtained by the binary coupling color mixing equation, and a chromaticity ring matching the classical Ostwald color solid is constructed, thereby solving the problem of uneven distribution of the chromaticity ring based on the actual color obtained, and the number of chromaticity of the chromaticity ring can be arbitrarily set, thereby breaking through the traditional 24 chromaticity, and the industrial design scene with different color matching requirements can be flexibly applied.

[0115] 2. The highest brightness gray and the lowest brightness gray obtained by actual dyeing are used as the south and north poles of the gray axis, the common ratio coefficient of eight times of equal ratio increment or equal ratio decrement of the color value of the south and north poles is obtained, and the brightness color value of the eight levels on the gray axis is obtained, and the actual gray axis is constructed in this way, thereby solving the problem of mismatching between the theoretically constructed gray axis and the actual brightness distribution of the gray axis.

[0116] 3. The highest brightness gray and the lowest brightness gray obtained by actual dyeing are used as the south and north poles of the gray axis, a brightness value nonlinear increment function between the color value of the south and north poles and the brightness color value of the eight levels on the gray axis are designed, and the actual gray axis is constructed in this way, thereby solving the problem of mismatching between the theoretically constructed gray axis and the actual brightness distribution of the gray axis.

[0117] 4. The highest lightness gray and the lowest lightness gray obtained by actual dyeing are taken as the south and north poles of the gray scale axis, and an equal difference increasing function of lightness value is designed between the color values of the south and north poles, and 8 levels of lightness color values on the gray scale axis are designed, so as to construct an actual gray scale axis, and solve the problem of mismatching between the theoretically constructed gray scale axis and the actual lightness distribution of the gray scale axis.

[0118] 5. In an embodiment, the Ostwald color solid is divided into 24 levels of hue, 8 levels of lightness, and 8 levels of purity, so as to obtain an Ostwald color solid composed of 8 white lightness planes, 8 black lightness planes, and 24 hue planes, and obtain the spatial coordinates of each grid point and the corresponding mixing concentration, color value and other parameters, and the mixing concentration of the primary color pigments corresponding to the grid point and the color value corresponding to the grid point can be obtained through the grid point coordinate value.

[0119] 6. In an embodiment, the mixing concentration value matrix and the color spectrum matrix of the 24 hue planes, the 8 black lightness planes and the 8 white lightness planes are given respectively with the grid point coordinates as the independent variable.

[0120] 7. In an embodiment, based on the color values of the six primary color pigments actually obtained, the color values of 36 color blocks on the 24 hue planes are obtained by taking hue as the reference through the constructed mathematical model, the color values of 169, 145, 121, 97, 73, 49, 25 and 1 color blocks on the 8 black lightness planes are obtained by taking the mixing concentration of the black pigment as the reference, and the color values of 169, 145, 121, 97, 73, 49, 25 and 1 color blocks on the 8 white lightness planes are obtained by taking the mixing concentration of the white pigment as the reference.

[0121] In the Ostwald full color domain gridding color mixing model constructed by the present application, the grid point position corresponding to each color is directly related to the mixing concentration ratio of the black, white and pure color pigments required for matching the color, so that the designer can directly obtain the proportion of each pigment from the model to obtain the target color when matching colors, and therefore the design efficiency is greatly improved. BRIEF DESCRIPTION OF DRAWINGS

[0122] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed in the embodiment description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor.

[0123] Figure 1 is the construction flowchart of the Ostwald full color domain gridding color mixing model of the present application.

[0124] Figure 2is a mixed color concentration contrast curve diagram of three gray scale axis construction modes in embodiment one of the present application.

[0125] Figure 3 is a flow chart of embodiment two of the present application. DETAILED DESCRIPTION

[0126] In order to make the purpose, technical solutions and advantages of the present application more clear, the embodiments of the present application are further described in detail below.

[0127] Embodiment one:

[0128] This embodiment takes a 24-hue ring as an example to further describe the specific derivation and construction process of the Ostwald full-color gamut gridding mixed color model of the present application, mainly including the following contents:

[0129] (I) Digital construction of Ostwald hue ring

[0130] 1. Acquisition of four primary colors of hue ring

[0131] There are four primary colors of red, yellow, green and blue in the Ostwald model, and their color values are respectively: C 2R =(r 2R ,g 2R ,b 2R ), C 2Y =(r 2Y ,g 2Y ,b 2Y ), C 2SG =(r 2SG ,g 2SG ,b 2SG ), C 2UB =(r 2UB ,g 2UB ,b 2UB ). In practical operation, the color values of red, yellow, green and blue dyes obtained based on actual dyeing and measurement can be used as the color values of the four primary colors of the Ostwald model, and the weights of the four primary colors of dyes are taken as Ω 2R ,Ω 2Y ,Ω 2SG ,Ω 2UB .

[0132] 2. Construction of two-by-two coupling mixed color mode of four primary colors

[0133] Take the four primary colors red, yellow, green and blue (denoted as 2R, 2Y, 2SG and 2UB) of the Ostwald color solid four-primary-color hue ring as the reference, and form four mixed color intervals (2R-2Y), (2Y-2SG), (2SG-2UB) and (2UB-2R). Take δ=1, 2, 3, 4, and set C α (δ) and C β(δ) and Ω α (δ) and Ω β (δ) represent the color values of the two end points of the four color mixing intervals respectively, Ω

[0134]

[0135] The end points C α (δ) and C β (δ) are coupled to mix colors, and the increasing gradient of the mixing concentration is set as 1 / ε m (ε m ≥2) and ε = 1, 2, …, ε m -1, ε m , then (ε α +1) grid points are obtained in the color mixing interval of C β (δ)-C m (δ), and the mixing concentration corresponding to each grid point is set as The weight of the corresponding mixed sample is Ω αβ (ε), and the color value corresponding to each grid point is C α-β (ε), then:

[0136]

[0137] 3. Ostwald color wheel spectrum construction

[0138] In order to obtain the Ostwald color wheel spectrum, four independent color mixing intervals (2R~2Y), (2Y~2SG), (2SG~2UB), and (2UB~2R) with four primary colors as nodes can be combined into a continuous cycle color mixing process, and ξ is the grid point number of the continuous cycle color mixing process (2R~2Y~2SG~2UB~2R), and:

[0139]

[0140] From equation (6), we have:

[0141]

[0142] Substitute equation (7) into equations (3), (4), and (5), and set ξ = 1, 2, …, (2ε m -1), 2ε m , 2ε m +1, …, (3ε m -1), 3ε m , (3ε m +1), …, (4ε m -1), 4ε mBased on the change of grid point index ξ, a continuous cyclic color mixing process (2R~2Y~2SG~2UB~2R) can be achieved.

[0143] Let the color mixing concentration value corresponding to each grid point index be... Let the weight of the mixed sample corresponding to each grid point number be Ω(ξ), and the color value corresponding to each grid point number be C(ξ). Then:

[0144]

[0145] At this point, the grid point index ξ and the color mixing concentration of the mixed sample are obtained. The correspondence between the weight Ω(ξ) of the mixed sample and the color value C(ξ) of the mixed sample. When ξ=1,2,...,(2ε m -1), 2ε m ,2ε m +1,...,(3ε m -1),3ε m ,(3ε m +1),...,(4ε m -1), 4ε m At this time, the color value corresponding to each grid point can be obtained, thereby obtaining the mixed color spectrum matrix. as follows:

[0146]

[0147] From equation (11), it can be seen that through the grid points ξ=1,2,...,(2ε m -1), 2ε m ,2ε m +1,...,(3ε m -1),3ε m ,(3ε m +1),...,(4ε m -1), 4ε m The change causes the colors of each grid point, C(ξ)=[r(ξ),g(ξ),b(ξ)], to achieve a continuous cyclic color mixing of (2R~2Y~2SG~2UB~2R) in sequence.

[0148] 4. Position coordinates of each hue on the color wheel

[0149] Based on equation (11), the color is divided into 4ε m A color spectrum composed of all hues, with each level of hue. Mapping this onto the hue circle of the Ostwald color solid yields 4ε. m grid point P ξ (θ ξ ,ρ ξ )=[θ ξ,ρ ξ ,z ξ The polar coordinates of each grid point in the polar coordinate system are θ. ξ The polar radius coordinates are ρ ξ The height coordinate is z ξ The coordinates of each grid point on the color wheel are as follows:

[0150]

[0151] When ξ=1,2,...,(2ε m -1), 2ε m ,2ε m +1,...,(3ε m -1),3ε m ,(3ε m +1),...,(4ε m -1), 4ε m At that time, the coordinate values ​​corresponding to each grid point are combined into a grid point coordinate matrix. as follows:

[0152]

[0153] or

[0154]

[0155] 5. Polar coordinates and color values ​​corresponding to the 24-hue color wheel of the Oswald color solid.

[0156] Based on equations (11) to (14), when ε m When θ = 6, the hue angle of the 24-hue color wheel of the Ostwald color solid is obtained as θ. ξ = (ξ-1)×15, then the position coordinates and color values ​​of each hue are shown in Table 1.

[0157] Table 1. Polar coordinates and color values ​​corresponding to the 24-hue color wheel of the Oswald color solid.

[0158]

[0159]

[0160] (II) Digital Construction of Ostwald Gray Axis

[0161] 1. The composition of the grayscale axis and the position coordinates of each node

[0162] Ostwald's grayscale axis is divided into 8 brightness levels, with gray O being the highest brightness. a The color standard is white, and the lowest lightness gray is O. p Marked as black on the color scale, and white on the color scale. a and color mark black Op Insert 6 gray levels of gray between them, get 8 nodes on the Ostwald gray axis with different lightness from low to high: O p ,O n ,O l ,O i ,O g ,O e ,O c ,O a Then the polar coordinate values of the 8 nodes on the gray axis are:

[0163]

[0164] 2、Obtain the mixing concentration of each node on the Ostwald gray axis

[0165] This embodiment constructs the gray axis of Ostwald color solid in three ways: the first is to obtain the lightness value (color value) of the eight nodes of the gray axis with equal ratio increment of lightness, and to construct the gray axis; The second is to obtain the lightness value (color value) of the eight nodes of the gray axis with equal difference increment of lightness, and to construct the gray axis; The third is to obtain the lightness value (color value) of the eight nodes of the gray axis with non-linear increment of lightness, and to construct the gray axis.

[0166] (1) Ostwald gray axis constructed by equal ratio increment of lightness

[0167] ①Obtain the common ratio coefficient of the gray axis with equal ratio increment of lightness

[0168] Assume that the lightness values of the eight nodes of the Ostwald defined gray axis are increased in equal ratio series, and the color values of the highest lightness and the lowest lightness nodes of the Ostwald color system gray axis are C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ], respectively, and the fiber weight is Ω W ,Ω K . The fiber weight corresponding to each node a, c, e, g, i, l, n, p of the gray axis is Let Ω W ,Ω K be mixed with mixing concentration to obtain Then the mixing concentration corresponding to the eight nodes of the gray axis is The corresponding color value is Then:

[0169]

[0170] Let the common ratio coefficient λ1 of the eight nodes of the Ostwald grayscale axis be:

[0171]

[0172] but:

[0173]

[0174] If r W =g W =b W ;r K =g K =b K ,but but:

[0175]

[0176] ② Obtaining the mixed concentration values ​​of each node on the grayscale axis with proportionally increasing brightness

[0177] make τ Taking values ​​of 1 and 2 respectively, according to formula (19), the common ratio coefficient λ1 is:

[0178]

[0179] From equation (22), we can obtain:

[0180]

[0181] Based on formula (23), we can obtain for:

[0182]

[0183] make τ Taking values ​​of 1 and 2 respectively, the common ratio coefficient λ1 can also be written as:

[0184]

[0185] From equation (25), we can obtain:

[0186]

[0187] Then we can obtain for:

[0188]

[0189] Concentration value matrix of 8 nodes on the grayscale axis It can be written as:

[0190]

[0191] ③ Obtaining the color values ​​of each node on the grayscale axis with proportionally increasing brightness

[0192] Based on equation (20), the highest luminance C that can be obtained in practice can be derived. W and lowest brightness C K Find the common ratio coefficient λ1. Let:

[0193]

[0194] Then, the color value C with the highest brightness on the grayscale axis W =[r W ,g W ,b W ] and the lowest brightness color value C K =[r K ,g K ,b K The color values ​​of the remaining nodes on the grayscale axis are obtained as follows:

[0195]

[0196] It can also be written as:

[0197]

[0198] Color values ​​based on 8 nodes of the grayscale axis The chromatographic matrix of the grayscale axis can be obtained. for:

[0199]

[0200] If the highest brightness color value C on the grayscale axis is known... W =[r W ,g W ,b W ] and the lowest brightness color value C K =[r K ,g K ,b K ], and convert it into a concentration value matrix Color mixing yields the blended colors of the 8 nodes on the grayscale axis. Then we can obtain:

[0201]

[0202] As can be seen from equation (33), based on the highest brightness C of the pigment W =[r W ,g W ,b W and minimum brightness C K =[r K ,gK ,b K ], and the mixture concentration corresponding to the eight nodes on the gray axis, the color values corresponding to the eight nodes on the gray axis can be obtained. If C W = [247, 247, 247], C K = [11, 11, 11], τ = 1, 2,..., 7, 8 are substituted into equations (22), (24), (27), we have:

[0203]

[0204] Then equation (35) can be written as:

[0205]

[0206] (2) Osterwald gray axis constructed by increasing the brightness equidifference

[0207] ①Obtain the equidifference coefficient of the brightness equidifference gray axis

[0208] Similarly, equations (16)-(18), let the color values of the eight nodes of the Osterwald model gray axis increase by the equidifference coefficient λ2, then:

[0209]

[0210] From equation (37), we have:

[0211]

[0212] If r W = g W = b W , r K = g K = b K , then C W , C K can be written as:

[0213]

[0214] Then equation (38) can be written as:

[0215]

[0216] ②Obtain the mixture concentration values of each node of the brightness equidifference gray axis

[0217] Let τ take 1 and 2 respectively, according to equations (37)-(40), then the equidifference coefficient λ2 can be written as:

[0218]

[0219] From equation (41), we have is:

[0220]

[0221] Let τ respectively, according to formula (37)-(40), the arithmetic progression coefficient λ2 can also be written as:

[0222]

[0223] From formula (43), we can get is:

[0224]

[0225] The concentration value matrix of the 8 nodes of the gray axis can be written as:

[0226]

[0227] ③Obtaining the color value of each node of the lightness arithmetic progression increasing gray axis

[0228] Based on formula (38)-(40), the color value of each node on the gray axis can be obtained as follows:

[0229]

[0230] Based on the color value of the 8 nodes of the gray axis the color spectrum matrix of the gray axis can be obtained is:

[0231]

[0232] (3) Construction of Ostwald gray axis with non-linear increase of lightness

[0233] ① Construction of non-linear coupling color mixing mode of the highest lightness white and the lowest lightness black

[0234] In order to make the lightness of the Ostwald gray axis change in a non-linear increasing rule, this patent realizes the non-linear coupling color mixing of the highest lightness white color material and the lowest lightness black color material, that is, under the premise of keeping the total weight of the mixed sample unchanged, the mixed concentration of the white color material is increased non-linearly, and the mixed concentration of the black color material is decreased non-linearly, so that the lightness value corresponding to the eight nodes of the gray axis changes in a non-linear rule from the lowest lightness to the highest lightness. Let the weight of the highest lightness white color material be Ω w , and the weight of the lowest lightness black color material be Ω k The above fiber weights are designed as follows:

[0235]

[0236] White color material Ω based on the highest brightness w Black color material Ω with the lowest brightness k , the nonlinear coupling color mixing mode is constructed as follows:

[0237]

[0238] When τ = 1, 2,..., 7, 8, respectively corresponding to the weight of the mixed sample of each node of the gray axis.

[0239] ② Acquisition of mixed concentration values of each node of the gray axis

[0240] Based on formula (49), when τ = 1, 2,..., 7, 8, the mixed concentration of each node of the gray axis corresponding to is:

[0241]

[0242] Based on formula (50), when τ = 1, 2,..., 7, 8, let the concentration values of white and black fibers in the mixed sample corresponding to each node of the gray axis are:

[0243]

[0244] By comparing the color mixing concentration matrix equations based on the three gray axis construction methods of brightness equal ratio increment, equal difference increment and nonlinear increment, as shown in Figure 2 The blue, green and red respectively correspond to the gray axis fiber mixing concentration curves under the construction methods of equal ratio increment, equal difference increment and nonlinear increment. The comparison of the values of the three construction methods can be obtained: with the increase of the brightness value τ, the proportion of the decrease of the black fiber mixing concentration of the nonlinear mode gradually increases, which is significantly smaller than the black fiber mixing concentration under the equal ratio and equal difference two modes, which can improve the color recognition of the color mixing blocks in practical application.

[0245] ③ Acquisition of color values of each node of the brightness nonlinear increment gray axis

[0246] Assuming that the brightness values of the eight nodes of the Ostwald definition gray axis are increased in a nonlinear manner, the color values of the highest brightness and the lowest brightness of the Ostwald color system gray axis are C W = [r W , g W , b W ] and C K = [r K , g K , b K ], respectively, and the fiber weight is Ω W , ΩK The fiber weight corresponding to each node a, c, e, g, i, 1, n, p of the gray axis is Ω O The color value corresponding to each node of the gray axis is Let Ω W Ω K The mixed concentration is The mixture obtained is And its color value is Then:

[0247]

[0248] The color value of the eight nodes of the Ostwald gray axis is :

[0249]

[0250] The digital color matrix of the color value of the eight nodes of the gray axis is :

[0251]

[0252] (4) The digital gray axis constructed based on the coordinate value of each node, the mixed concentration, and the color value

[0253] Based on the foregoing derivation, the construction of the Ostwald gray axis can be divided into three modes. Let η = 1 correspond to the gray axis constructed with the common ratio coefficient increasing, η = 2 correspond to the gray axis constructed with the arithmetic coefficient increasing, and η = 3 correspond to the gray axis constructed with the non-linear increasing.

[0254] Let the mixed weight matrix of the eight nodes of the gray axis obtained by the three construction modes be Then:

[0255]

[0256] Let the concentration value matrix of the eight nodes of the gray axis obtained by the three construction modes be Then:

[0257]

[0258] Let the color value matrix of the eight nodes of the gray axis obtained by the three construction modes be Then:

[0259]

[0260] Based on formulas (14), (48), and (49), the position coordinates, mixed concentrations, and color values corresponding to each node of the gray axis can be obtained, as shown in Table 2. The digital gray axis is constructed therefrom.

[0261] Table 2 Color values and position coordinates of each node of the gray axis

[0262]

[0263] (III) Digital construction of Ostwald color solid

[0264] Ostwald color solid is a three-dimensional biconic color solid, mainly including a series of isochromatic phase planes, isogreyscale planes, and isowhiteness planes. The main task of constructing Ostwald color solid is to construct its 24 isochromatic phase planes, 8 isowhiteness planes, 8 isoblackness planes, and 8 isopurity planes. Among them, the construction of isochromatic phase planes is the key.

[0265] In constructing isochromatic phase planes, a color phase in 24 color phases is mixed with white and black to obtain the color phase plane. There are two mixing methods. One mixing method is to combine a certain proportion of black on the gray axis with another eight different proportions of white, eliminate combinations with a mixing concentration proportion greater than 100%, leave combinations with a mixing concentration proportion less than or equal to 100%, and fill the part with a mixing concentration proportion less than 100% with the color phase in the remaining proportion, thereby obtaining eight isoblackness lines on the isochromatic phase plane and further obtaining an isochromatic phase plane containing 36 color blocks. The other mixing method is to combine a certain proportion of white on the gray axis with another eight different proportions of black, eliminate combinations with a mixing concentration proportion greater than 100%, leave combinations with a mixing concentration proportion less than or equal to 100%, and fill the part with a mixing concentration proportion less than 100% with the color phase in the remaining proportion, thereby obtaining eight isowhiteness lines on the isochromatic phase plane and further obtaining an isochromatic phase plane containing 36 color blocks.

[0266] The construction of isowhiteness planes is to regard the conical surface formed by all grid points with equal white content on each color phase plane as an isowhiteness plane, thereby obtaining eight isowhiteness planes.

[0267] The construction of isoblackness planes is to regard the conical surface formed by all grid points with equal black content on each color phase plane as an isoblackness plane, thereby obtaining eight isoblackness planes.

[0268] The construction of isopurity planes is to regard the cylindrical surface formed by all grid points with equal (white + black) content, i.e., equal color content, on each color phase plane as an isopurity plane, thereby obtaining eight isopurity planes.

[0269] 1. Base color combination design of each isochromatic phase plane

[0270] According to equal hue, equal whiteness, equal blackness and equal purity, the Ostwald color solid is divided into grids, and the hue serial number of each grid point is ξ = 1, 2,..., 23, 24; the blackness serial number is δ = 1, 2,..., 7, 8; the whiteness serial number is τ = 1, 2,..., 7, 8; and the purity serial number is γ = 1, 2,..., 7, 8. The color block in the Ostwald color solid is formed by mixing a certain hue (denoted by C ξ ) with the highest brightness white (denoted by C W ) and the lowest brightness black (denoted by C K ). The base color combination of each equal hue plane is shown in Table 3.

[0271] Table 3 Base color combination of Ostwald color solid equal hue plane

[0272]

[0273] 2. Construction of three-base color grid mixed color mode based on equal hue plane

[0274] On the ξth equal hue plane of the Ostwald color solid, let the color values of the three base colors be C ξ = (r ξ , g ξ , b ξ ), C W = (r W , g W , b W ), and C K = (r K , g K , b K ), and the weights be Ω ξ , Ω W , Ω K . When the equal hue plane θ ξ = (ξ-1) × 15 (ξ = 1, 2,..., 23, 24), let Ω(τ, δ, ξ) = Ω ξ = Ω W = Ω K = Ω, and according to the mixing rule of each color block on the Ostwald equal hue plane, the coupling mixing mode corresponding to each grid point is constructed as follows:

[0275]

[0276] Based on formula (60), the three-base color mixing concentration in the mixed color sample Ω η (τ, δ, ξ) can be obtained by the following formula:

[0277]

[0278] or written as:​

[0279]

[0280] The color value corresponding to the grid point P(τ, δ, ξ) can be obtained as

[0281]

[0282] 3. Grid Osterberg color cube based on digitalized color mixing construction

[0283] Let where (τ, δ, ξ) is the grid point number. According to the construction method of Osterberg color cube, eliminate the color blocks with mixing concentration ratio greater than 100% in the Osterberg color cube constructed under the three modes of equal ratio, equal difference and non-linear when η = 1, 2, 3, and the mixing matrix of the remaining 680 color blocks can be expressed as follows:

[0284]

[0285] Equation (64) is the quality matrix of the 680 color mixing color blocks of the Osterberg color cube.

[0286]

[0287] Equation (65) is the color mixing concentration matrix of the 680 color mixing color blocks of the Osterberg color cube.

[0288]

[0289] Equation (66) is the color matrix of the 680 color mixing color blocks of the Osterberg color cube.

[0290] 4. Construction of Osterberg color cube isochromatic surface and its color spectrum matrix

[0291] In the Osterberg model, mixing one color with white and black in 24 hues can obtain isochromatic surfaces, a total of 24 isochromatic surfaces, each containing 36 color blocks.

[0292] Let ξ = 1, 2,..., 23, 24, when η = 1, 2, 3, the color matrix of the 24 isochromatic surfaces defined by the grid point number is as follows:

[0293]

[0294] 5. Construction of Osterberg color cube isochromatic surface and its color spectrum matrix

[0295] ​In Ostwald model, the conical surface composed of all grid points with equal white content on each hue plane is regarded as the equal-white plane, and 8 equal-white planes can be obtained, each of which contains 169, 145, 121, 97, 73, 49, 25, 1 color blocks respectively. Let ξ = 1, 2,..., 23, 24, the color matrix of the 8 equal-white planes is shown as follows.

[0296] The color matrix of the equal-white plane when τ = 1 is as follows:

[0297]

[0298] The color matrix of the equal-white plane when τ = 2 is as follows:

[0299]

[0300] The color matrix of the equal-white plane when τ = 3 is as follows:

[0301]

[0302] The color matrix of the equal-white plane when τ = 4 is as follows:

[0303]

[0304] The color matrix of the equal-white plane when τ = 5 is as follows:

[0305]

[0306] The color matrix of the equal-white plane when τ = 6 is as follows:

[0307]

[0308] The color matrix of the equal-white plane when τ = 7 is as follows:

[0309]

[0310] The color matrix of the equal-white plane when τ = 8 is as follows:

[0311]

[0312] 6、Ostwald color solid equal-black plane and construction of its color spectrum matrix

[0313] In Ostwald model, the conical surface composed of all grid points with equal black content on each hue plane is regarded as the equal-black plane, and 8 equal-black planes can be obtained, each of which contains 1, 25, 49, 73, 97, 121, 145, 169 color blocks respectively. Let ξ = 1, 2,..., 23, 24, the color matrix of the 8 equal-white planes is shown as follows.

[0314] The color matrix of the equal brightness surface when δ = 1 is:

[0315]

[0316] The color matrix of the equal brightness surface when δ = 2 is:

[0317]

[0318] The color matrix of the equal brightness surface when δ = 3 is:

[0319]

[0320] The color matrix of the equal brightness surface when δ = 4 is:

[0321]

[0322] The color matrix of the equal brightness surface when δ = 5 is:

[0323]

[0324] The color matrix of the equal brightness surface when δ = 6 is:

[0325]

[0326] The color matrix of the equal brightness surface when δ = 7 is:

[0327]

[0328] The color matrix of the equal brightness surface when δ = 8 is:

[0329]

[0330] 7. Ostwald color solid equal purity surface and construction of its color matrix

[0331] The construction of the equal purity surface is to take all the grid points on each hue surface with equal (white + black) content, i.e. equal color content, to form a cylindrical surface as the equal purity surface, from which eight equal purity surfaces can be obtained. Take the (white + black) content minimum (0%) as the first equal purity surface, at this time γ = 1, and so on. Let τ, δ = 1, 2,..., 8; ξ = 1, 2,..., 23, 24, then the color matrix of the eight equal purity surfaces is as follows.

[0332] The color matrix of the equal purity surface when γ = 1 is:

[0333]

[0334] The color matrix of the equal purity surface when γ = 2 is:

[0335]

[0336] The color matrix of the equal-purity surface when γ = 3 is:

[0337]

[0338] The color matrix of the equal-purity surface when γ = 4 is:

[0339]

[0340] The color matrix of the equal-purity surface when γ = 5 is:

[0341]

[0342] The color matrix of the equal-purity surface when γ = 6 is:

[0343]

[0344] The color matrix of the equal-purity surface when γ = 7 is:

[0345]

[0346] The color matrix of the equal-purity surface when γ = 8 is:

[0347]

[0348] Example 2:

[0349] In this example, the color values of the six primary color pigments obtained by dyeing and color measurement are: red (235, 12, 23), yellow (235, 225, 32), green (23, 224, 36), blue (28, 32, 238), black (11, 11, 11), and white (247, 247, 247). With the aid of the Ostwald color solid color spectrum matrix model and algorithm constructed independently, the design of the 24-hue color spectrum of the Ostwald color solid is completed; the design of the eight-grade gray axis color spectrum with equal-ratio increment of lightness; the design of the second Y isochromatic hue surface, the seventh equal-white surface, the second equal-black surface, and the third equal-purity surface color spectrum. The specific process is shown in Figure 3 , and specifically includes the following steps.

[0350] (1) Design of the hue ring color spectrum of the Ostwald model

[0351] Let the color values of the four primary colors of red, yellow, green, and blue in the Ostwald model be: C 2R = (235, 12, 23), C 2Y = (235, 225, 32), C 2SG = (23, 224, 36), and C 2UB= (28, 32, 238), the color values of 24 hues can be obtained as shown in Table 4.

[0352] Table 4 Color values of 24 hues of Ostwald color system

[0353]

[0354]

[0355] (2) Design of equi-interval gray scale chroma of Ostwald model

[0356] Let the color values of the highest lightness gray O a and the lowest lightness gray O p be C W = (247, 247, 247) and C K = (11, 11, 11) respectively, and the common ratio coefficient λ1 of the color values of 8 nodes on the gray scale can be obtained according to the algorithm of the gray scale model with equi-interval lightness:

[0357]

[0358] The position coordinates and color values of the 8 nodes on the gray scale are shown in Table 5.

[0359] Table 5 Color values and position coordinates of 8 nodes on the gray scale

[0360]

[0361] Meanwhile, the white-black mixing concentrations of the remaining 7 nodes on the gray scale can be obtained as shown in Table 6.

[0362] Table 6 White-black mixing ratios of 8 nodes on the gray scale

[0363] Reference Signs O p ]]> <![CDATA[O n ]]> O l ]]> O i ]]> O g ]]> O e ]] O c ]]> O a <!-- 28 -->]]> White Proportion 0 2.6％ 6.7％ 13％ 22.9％ 38.4％ 62.4％ 100％ Black Proportion 1 97.4％ 93.3％ 87％ 77.1％ 61.6％ 37.6％ 0

[0364] (3) Design of equihue surface chroma of Ostwald model

[0365] Take the 2Y yellow hue surface as an example, where ξ = 2, and the color value C 2Y = (235, 225, 32) is known, the color values of each color block on the 2Y hue surface can be obtained according to the mixing concentration formula (67) on the equihue surface as shown in Table 7, and the color values of each color block on the other 23 equihue surfaces can be obtained in the same way.

[0366] Table 7 Color values of each color block on the yellow 2Y hue surface

[0367]

[0368]

[0369] (4) Ostwald model isohyperchromatic surface color spectrum design

[0370] Taking the 7th isohyperchromatic surface as an example, τ=7, the color values of each color block on the isohyperchromatic surface can be obtained based on the mixing concentration of each color block, as shown in Table 8, and the color values of each color block on the other seven isohyperchromatic surfaces can be obtained in the same way.

[0371] Table 8 Color values of each color block on the seventh isohyperchromatic surface

[0372]

[0373]

[0374] (5) Ostwald model isohypochromatic surface color spectrum design

[0375] Taking the 2nd isohypochromatic surface as an example, δ=2, the color values of each color block on the isohypochromatic surface can be obtained based on the mixing concentration of each color block, as shown in Table 9, and the color values of each color block on the other seven isohypochromatic surfaces can be obtained in the same way.

[0376] Table 9 Color values of each color block on the seventh isohypochromatic surface

[0377]

[0378]

[0379] (6) Ostwald model isopurity surface color spectrum design

[0380] Taking the 3rd isopurity surface as an example, γ=3, the color values of each color block on the isopurity surface can be obtained based on the mixing concentration of each color block, as shown in Table 10, and the color values of each color block on the other seven isopurity surfaces can be obtained in the same way.

[0381] Table 10 Color values of each color block on the third isopurity surface

[0382]

[0383]

[0384] Some steps in the embodiments of the present application can be implemented by software, and the corresponding software program can be stored in a readable storage medium, such as an optical disc or a hard disk.

[0385] The above only describes the preferred embodiments of the present application, and is not intended to limit the present application, and any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application shall be included in the protection scope of the present application.

Claims

1. A color mixing method, characterized in that, The method includes: Step 1: Based on actual color measurement, obtain the color values ​​of the six primary color pigments of the Ostwald color solid, and construct the Ostwald full-gamut gridded color mixing model; Step 2: Obtain the target color from the Ostwald full-gamut meshed color mixing model, and obtain the mixing concentration values ​​of the six primary colors required to mix the target color based on the grid point coordinates of the target color in the model; Step 3: Based on the mixing concentration obtained in Step 2, mix the six primary color pigments to obtain the target pigment; The construction process of the Ostwald full-gamut meshed color mixing model includes: Step 1: Using the color values ​​of red, yellow, green, and blue pigments as a basis, construct the color wheel of the Ostwald color solid, including: Four color mixing intervals are constructed. Binary coupling color mixing is performed based on the color values ​​at the two endpoints of each interval to obtain the color value C(ξ) of each hue on the Oswald color solid hue wheel, expressed as: The mixing concentration of each hue on the Oswald color solid color wheel is as follows: Wherein, ξ represents the grid point number of each hue on the Ostwald color solid hue wheel, and ε m ≥2; δ=1, 2, 3, 4 represent the indices of the four color mixing intervals; [r α (δ),g α (δ),b α (δ)]、[r β (δ),g β (ξ),b β [(δ)] represent the color values ​​of the two endpoints of the δth color mixing interval, respectively; and Let C represent the endpoint colors of the color mixing interval corresponding to the ξ-th grid point. α (δ) and C β The mixing concentration of (δ), wherein the increasing gradient of the mixing concentration is 1 / ε m ; The coordinates of the grid points for each hue on the Ostwald color solid are: Where, θ ξ ρ ξ and z ξ These represent the polar angle, polar radius, and height of the grid point, respectively. Step 2: Using the color values ​​of white and black pigments as the color values ​​of the highest and lowest brightness nodes, construct the grayscale axis of the Ostwald color solid. The polar coordinates of the 8 nodes on the grayscale axis are: The color values ​​and mixing concentrations corresponding to the eight nodes on the grayscale axis are obtained by increasing proportionally, arithmetically, or nonlinearly. The concentration value matrix of the 8 nodes of the grayscale axis is as follows: The color value matrix of the 8 nodes of the grayscale axis is as follows: Among them, η=1 corresponds to a grayscale axis constructed with increasing common ratio coefficient, η=2 corresponds to a grayscale axis constructed with increasing arithmetic coefficient, and η=3 corresponds to a grayscale axis constructed with non-linear increasing coefficient. Step 3: Using a certain hue color value C ξ =(r ξ ,g ξ ,b ξ ), highest brightness node color value C W =(r W ,g W ,b W ) and the color value C of the lowest brightness node K =(r K ,g K ,b K Construct the isohue planes of the Ostwald color solid, and the three primary color mixing concentrations of the color patches on each isohue plane are: The color value of the color patch corresponding to grid point P(τ,δ,ξ) for: Where δ = 1, 2, ..., 7, 8 are the iso-blackness indices, τ = 1, 2, ..., 7, 8 are the iso-whiteness indices, γ = 1, 2, ..., 7, 8 are the iso-purity indices, and Ω... ξ ,Ω W ,Ω K These represent the weights of the pigments for the three primary colors.

2. The color matching method according to claim 1, characterized in that, Step 2 involves constructing a grayscale axis using a proportionally increasing brightness method. The color mixing density value matrix of the eight nodes of the grayscale axis is as follows: in, The color values ​​of each color node on the grayscale axis are: Among them, C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ] represent the color values ​​of the highest and lowest brightness nodes, respectively.

3. The color matching method according to claim 1, characterized in that, Step 2 uses an arithmetic progression of brightness to construct the Ostwald grayscale axis, and the color mixing density matrix of the eight nodes of the grayscale axis is as follows: The color values ​​of each color node on the grayscale axis are: Among them, C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ] represent the color values ​​of the highest and lowest brightness nodes, respectively.

4. The color matching method according to claim 1, characterized in that, Step 2 uses a non-linear brightness increment method to construct the grayscale axis. The color mixing density value matrix of the eight nodes of the grayscale axis is as follows: The color values ​​of each color node on the grayscale axis are: Among them, C W =[r W ,g W ,b W ] and C K =[r K ,g K ,b K ] represent the color values ​​of the highest and lowest brightness nodes, respectively.

5. The color matching method according to claim 1, characterized in that, When ξ = 1, 2, ..., 23, 24, the isowhiteness surfaces of the Ostwald color solid are as follows: The color matrix of the iso-whiteness surface when τ = 1 is: The color matrix of the iso-whiteness surface when τ = 2 is: The color matrix of the iso-whiteness surface when τ = 3 is: The color matrix of the iso-whiteness surface when τ=4 is: The color matrix of the iso-whiteness surface when τ=5 is: The color matrix of the iso-whiteness surface when τ = 6 is: The color matrix of the iso-whiteness surface when τ = 7 is: The color matrix of a homoluminous surface when τ = 8 is: 。 6. The color matching method according to claim 1, characterized in that, When ξ = 1, 2, ..., 23, 24, the iso-echo surfaces of the Ostwald color solid are as follows: The color matrix of the iso-echo surface when δ = 1 is: The color matrix of the iso-echo surface when δ = 2 is: The color matrix of the iso-echo surface when δ=3 is: The color matrix of the iso-echo surface when δ=4 is: The color matrix of the iso-echo surface when δ=5 is: The color matrix of the iso-echo surface when δ=6 is: The color matrix of the iso-echo surface when δ = 7 is: The color matrix of the iso-echo surface when δ = 8 is: 。 7. The color mixing method according to claim 1, characterized in that, When ξ = 1, 2, ..., 23, 24, the isopurity surfaces of the Ostwald color solid are as follows: The color matrix of the isopurity surface when γ = 1 is: The color matrix of the isopurity surface when γ = 2 is: The color matrix of the isopurity surface when γ = 3 is: The color matrix of the isopurity surface when γ = 4 is: The color matrix of the isopurity surface when γ = 5 is: The color matrix of the isopurity surface when γ = 6 is: The color matrix of the isopurity surface when γ = 7 is: The color matrix of the isopurity surface when γ = 8 is:

8. A color matching system, characterized in that, The system for implementing the color mixing method as described in any one of claims 1-7 includes: a color measuring device, an Ostwald full-gamut gridded color mixing model construction module, and a visualization module; The color measuring device is used to obtain the color values ​​of the six primary colors of Ostwald solid color in actual pigments; The Ostwald full-gamut meshed color mixing model construction module includes: The color value acquisition module is configured to acquire the actual color value of the colorant from the color measuring device; The hue wheel construction module is configured to construct an equatorial hue wheel based on the actual acquired Ostwald six primary colors, providing the color value, position coordinate value, and mixing concentration of each hue. The grayscale axis construction module is configured to construct a grayscale axis based on the actual acquired black and white as the lowest and highest brightness nodes, respectively, and to provide the color value, position coordinate value, and mixing density of black and white for each brightness node; The isochromatic surface chromatogram construction module is configured to construct the 4ε Ostwald chromatic solid. m A colorimetric spectrum; The visualization module is configured to display the model color spectrum and output the mixing concentration and color value corresponding to the six primary colors based on the grid point coordinates of the target color.

9. The color matching system according to claim 8, characterized in that, The Ostwald full-gamut meshed color mixing model construction module also includes: The isowhiteness surface chromatography construction module is configured to be based on the 4ε m A colorimetric surface spectrum, in which the conical surface formed by all grid points with equal white content on each colorimetric surface is used as the isowhiteness surface of the Ossian color solid; The isocyanate surface chromatography construction module is configured to be based on the 4ε m A colorimetric spectrum with equal black content on each colorimetric surface is used as the isoblackness surface of the austenitic solid, which is a conical surface composed of all grid points with equal black content on each colorimetric surface. The isopurity surface chromatography construction module is configured to be based on the 4ε m An isochromatic color plane chromatogram uses the cylindrical surface formed by all grid points with equal color content on each color plane as the isopurity surface of the austenitic color solid.

10. A computer-readable storage medium, characterized in that, The storage medium stores a computer program that, when executed by a processor, implements the method as described in any one of claims 1 to 7.

Images