CN107067444B - Optimized spectrum color gamut mapping method - Google Patents

Optimized spectrum color gamut mapping method Download PDF

Info

Publication number
CN107067444B
CN107067444B CN201710286928.4A CN201710286928A CN107067444B CN 107067444 B CN107067444 B CN 107067444B CN 201710286928 A CN201710286928 A CN 201710286928A CN 107067444 B CN107067444 B CN 107067444B
Authority
CN
China
Prior art keywords
spectral
color
point
lms
spectrum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710286928.4A
Other languages
Chinese (zh)
Other versions
CN107067444A (en
Inventor
孔玲君
孙叶维
杨晟炜
张建青
周颖梅
方恩印
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Publishing and Printing College
Original Assignee
Shanghai Publishing and Printing College
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 Shanghai Publishing and Printing College filed Critical Shanghai Publishing and Printing College
Priority to CN201710286928.4A priority Critical patent/CN107067444B/en
Publication of CN107067444A publication Critical patent/CN107067444A/en
Application granted granted Critical
Publication of CN107067444B publication Critical patent/CN107067444B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention provides an optimized spectral gamut mapping method, which comprises the following steps: firstly, the method comprises the following steps: according to the visual characteristics of human eyes, selecting a proper LMS cone response spectrum sensitivity curve as a weighting function to perform weighting processing on high-dimensional spectrum data; II, secondly: extracting the first three principal components of the spectrum after weighting processing by adopting a principal component analysis method, and constructing an LMS-PCA spectrum link space; thirdly, the method comprises the following steps: describing the maximum spectral boundaries of the output device in LMS-PCA space; fourthly, the method comprises the following steps: describing specific spectral boundaries of the output device in LMS-PCA space; fifthly: judging whether the color point of the spectrum of the image to be mapped needs to be mapped or not; sixthly, the method comprises the following steps: mapping the image spectral color points beyond the spectral color gamut of the output equipment by adopting an LSLINceLmax color gamut compression algorithm; seventhly, the method comprises the following steps: and performing spectral reconstruction based on the color data in the LMS-PCA space obtained by mapping processing. The spectral color gamut mapping method has high spectral precision and chromaticity precision, and has stable color difference precision under variable observation environments.

Description

Optimized spectrum color gamut mapping method
Technical Field
The invention relates to the fields of color display, reproduction and multispectral color management, in particular to a spectral gamut mapping method in the field of multispectral color reproduction and reproduction.
Background
At present, researches on spectral color gamut mapping at home and abroad are still few, and the mainstream method is to construct a low-dimensional ICS space based on Principal Component Analysis (PCA) or LabPQR dimension reduction method, then design a spectral color gamut description method in the low-dimensional ICS space, and finally complete the mapping process from an image to equipment by using a mapping algorithm. However, in the spectrum color gamut mapping model based on the PCA space, because the visual characteristics of human eyes are not considered, the spectrum error before and after mapping is larger, and the mapping effect is not ideal; the spectral gamut mapping model based on the LabPQR space has a relatively complex calculation method, so that the mapping efficiency is reduced. The invention can effectively solve the problems, improve the spectral accuracy and the chromaticity accuracy of mapping on the basis of ensuring the mapping efficiency, and obtain an image with better subjective visual perception after mapping.
Disclosure of Invention
The present invention overcomes the above-mentioned deficiencies of the prior art and provides an optimized spectral gamut mapping method.
The invention adopts the following technical scheme:
the invention provides an optimized spectral gamut mapping method, which is characterized by comprising the following steps: comprises the following steps:
the method comprises the following steps: according to the visual characteristics of human eyes, selecting a proper LMS cone response spectrum sensitivity curve as a weighting function to perform weighting processing on high-dimensional spectrum data;
step two: extracting the first three principal components of the spectrum after weighting treatment by adopting a Principal Component Analysis (PCA) method, and constructing an LMS-PCA spectrum link space (ICS);
step three: describing the maximum spectral boundary GBD of the output device in LMS-PCA space;
step four: describing the specific spectral boundaries LBD of the output device in LMS-PCA space;
step five: judging whether the color point of the spectrum of the image to be mapped needs to be mapped or not;
step six: mapping the image spectral color points beyond the spectral color gamut of the output equipment by adopting an LSLINceLmax color gamut compression algorithm;
step seven: and performing spectral reconstruction based on the color data in the LMS-PCA space obtained by mapping processing.
Further, an optimized spectral gamut mapping method according to the present invention may further have the following features, wherein the first step includes the steps of:
step 1-1: and according to the visual characteristics of human eyes, selecting a proper LMS cone response spectrum sensitivity curve as a weighting function to perform weighting processing on the high-dimensional spectrum data. When the spectrum tristimulus values are converted into LMS cone responses, selecting a Bradford transformation matrix as a transformation matrix:
Figure BDA0001280915630000021
step 1-2: the LMS weighting function is constructed as equation 2:
Figure BDA0001280915630000022
in the formula, p is a variable coefficient and ranges from-0.5 to 50, and k is 1-10.
Step 1-3: and carrying out weighting processing on the high-dimensional spectral data by utilizing an LMS weighting function.
Further, an optimized spectral gamut mapping method according to the present invention may further have the following features, in step two: and extracting the first three principal components of the spectrum after weighting processing by adopting a principal component analysis method, and constructing an LMS-PCA spectrum link space.
Further, the optimized spectral gamut mapping method of the present invention may further have the following characteristics, and the step two specifically includes the following steps:
step 2-1: performing singular value decomposition on the weighted spectrum data;
step 2-2: and reducing the dimension of the weighted spectral data through the eigenvector of singular value decomposition, and extracting the first three principal components to construct an LMS-PCA rectangular coordinate system as a spectral link space.
Further, an optimized spectral gamut mapping method according to the present invention may further have the following features, where the step three includes the steps of:
step 3-1: using a rectangular coordinate system P of LMS-PCA space1,P2,P3And converting the coordinate system into a spherical coordinate system α, theta and r, and calculating spherical coordinate values corresponding to the spectral color points of the output device.
And 3-2, uniformly dividing the LMS-PCA space into n multiplied by n equal parts along α and theta directions respectively, calculating the spherical coordinates of each output device spectral color sample point, determining which sub-region of the LMS-PCA space each spectral color sample point is located in according to the values of α and theta, comparing the distances r from all spectral color sample points in each sub-region to the central point, keeping the spectral color sample point with the maximum r value as the boundary point of the sub-region, and finally storing all the obtained boundary points into a matrix to obtain the maximum spectral boundary GBD.
Further, an optimized spectral gamut mapping method according to the present invention may further have the following features, where the step four includes the steps of:
step 4-1, calculating a spherical coordinate α corresponding to the hue angle of the spectrum color point of the image to be mapped, and calculating the color point according to the value α and the spectrum plane coordinate (C)o,p1) And a center point E, finding α planes that are equal to the spectral color point, wherein the coordinates
Figure BDA0001280915630000031
Step 4-2, finding out a pair of adjacent spectral color points from the GBD submatrix of the first horizontal division layer according to the hue angle α, wherein the hue angle of one point is more than α, and the hue angle of the other point is less than α, then connecting the two adjacent spectral color points into a straight line, and calculating the intersection points of the straight line and the equal α plane, because one horizontal division layer has n total, n intersection points can be finally obtained, the intersection points are the linear boundaries of the spectral color gamut of the output device on the equal α plane,
step 4-3: and obtaining the maximum and minimum linear boundary points in the spectral color sample points of the output device according to the spectral color sample points with the maximum lightness value and the minimum lightness value in the maximum spectral boundary GBD matrix, and combining the maximum and minimum linear boundary points with the linear boundary obtained in the step 4-2 to form a specific spectral boundary LBD of the output device.
Further, the optimized spectral gamut mapping method of the present invention may also have the following characteristics, and the specific steps in step five are as follows:
and judging the relation between the mapping color point and the spectral color gamut of the output device. If the mapped color point is outside the spectral color gamut of the output device, processing according to a sixth step; if the mapped color point is within the spectral gamut of the output device, no processing is performed and the seventh step is entered directly.
Further, an optimized spectral gamut mapping method according to the present invention may further have such a feature that the method for determining whether the mapped color point is located inside or outside the spectral gamut of the output device in step five is as follows:
coordinate point M (C) of spectral plane of point Mo,p1) And p1The center point E of the shaft is connected to and extends to a straight line intersecting the LBD linear boundary at point H (C'o,p'1) If C iso<C'oAnd p ismin<p'1<pmaxThen the M points are within the printer spectral gamut, otherwise they are outside the printer spectral gamut.
Further, an optimized spectral gamut mapping method according to the present invention may further have the following feature, in step six, including the steps of:
step 6-1: and mapping the color points of the spectral image to be mapped of the spectral color gamut of the super-output device into the spectral domain of the output device by utilizing the LSLINCELmax color gamut compression in the LMS-PCA-based space. Based on LSLINceLmax color gamut compression in LMS-PCA space, on the basis of a traditional LSLIN algorithm, the influence of a spectral color gamut median value of output equipment and the influence of a spectral image brightness maximum value to be mapped are added, and meanwhile, the spectral chroma and the spectral brightness of a spectral color point of a spectral image to be mapped are compressed, so that a mapping process with higher precision is completed. The spectrum mapping expression is as follows:
Figure BDA0001280915630000041
in the formula, CrFor the spectral chroma value, C, of the mapped output imager(max)Maximum spectral chroma value, C, of the spectral gamut of the output deviceoAs a spectral color point chroma value, C, of the spectral image to be mappedo(max)Is the maximum spectral chroma value of the spectral image gamut to be mapped. p is a radical ofrFor the spectral brightness value, p, of the mapped output image1The spectral image brightness value is the spectral color point brightness value of the spectral image to be mapped, the center is the brightness value of the central point of the equal hue surface spectral brightness axis corresponding to the spectral color point of the spectral image to be mapped, and the cemmax is the weighting function of the maximum brightness value Lmax and the center of the spectral image to be mapped.
Step 6-2: the cemmax expression is as follows:
Figure BDA0001280915630000051
in the formula, Lmax is the maximum brightness value of the spectral image to be mapped, cc, dd and hh are constants, the value range of cc is-15 to 15, the value range of dd is 0 to 10, and the value range of hh is 1 to 3.
Advantageous effects of the invention
According to the invention, the high-dimensional spectral data is weighted by using a proper LMS weighting function, and the first three principal components of the weighted spectrum are extracted by using a principal component analysis method so as to construct an LMS-PCA spectral link space. And introducing a partition mature SMGBD algorithm into an LMS-PCA space to describe the spectral color gamut of the output device, and mapping the spectral image color point beyond the spectral color gamut of the output device into the spectral color gamut of the output device by adopting LSLINceLmax color gamut compression. The invention fully considers the visual characteristics of human eyes, constructs a spectral color gamut mapping model on the basis of optimizing the traditional LSLIN algorithm, has higher spectral precision and chromaticity precision, has stable color difference precision under variable observation environment, has better color difference stability simultaneously, is suitable for various types of spectral images, and has better practicability.
Drawings
Fig. 1 is a schematic diagram of segmentation.
FIG. 2 is C of a linear boundaryoP1Plane dieDrawing.
Fig. 3 is a diagram illustrating mapping determination.
Detailed Description
The following describes embodiments of the present invention with reference to the drawings.
In this example, the source gamut of the spectral gamut mapping is the spectral gamut of the standard spectral image, the target gamut is the spectral gamut of the HPDesign Z3200 multi-channel printer, and the mapping light source is set to D50. And obtaining detailed spectral color information of the HPdesign Z3200 printer through a self-made training color target. The optimized spectral gamut mapping method comprises the following specific steps:
firstly, according to the visual characteristics of human eyes, a proper LMS cone response spectrum sensitivity curve is selected as a weighting function, and weighting processing is carried out on high-dimensional spectrum data. The first step specifically comprises the following steps:
step 1-1: the bradford transformation Matrix (MBFD) is selected as the transformation matrix, and the spectral tristimulus values are converted into LMS cone responses according to equation 1.
Figure BDA0001280915630000061
Step 1-2: the LMS weighting function is constructed as equation 2:
Figure BDA0001280915630000062
step 1-3: and respectively carrying out weighting processing on the standard spectral image to be mapped and the high-dimensional spectral data of the HP Design Z3200 multi-channel printer by utilizing the LMS weighting function. If the spectral information of the printer is R and the dimension is N, the weighted spectral expression is as follows:
Figure BDA0001280915630000063
in the formula (I), the compound is shown in the specification,
Figure BDA0001280915630000064
is the LMS weighting function
Figure BDA0001280915630000065
A diagonal matrix of transitions.
Step two: and performing dimensionality reduction on the spectrum after weighting treatment by adopting a Principal Component Analysis (PCA), extracting the first three principal components, and constructing an LMS-PCA space. The second step specifically comprises the following steps:
step 2-1: the weighted spectrum is subjected to a singular value decomposition,
Figure BDA0001280915630000066
obtaining a characteristic vector U ═ U1,U2,...UnExtracting the first three feature vectors, and recording the vectors as
Figure BDA0001280915630000067
Step 2-2: by using
Figure BDA0001280915630000068
Reducing the dimension of the weighted spectrum to obtain a three-dimensional weighted spectrum after dimension reduction
Figure BDA0001280915630000069
Figure BDA00012809156300000610
Constructing an LMS-PCA rectangular coordinate system space by using the first three principal components obtained by dimensionality reduction, and setting the first three principal components as P1,P2,P3
Step three: the maximum spectral boundary GBD of the output device is described in LMS-PCA space. The third step specifically comprises the following steps:
step 3-1: the LMS-PCA space is transformed from a rectangular coordinate system P1,P2,P3And converted to a spherical coordinate system α, theta, r.
r and θ are calculated as follows:
Figure BDA0001280915630000071
α, the calculation of the value is calculated according to different formulas according to different quadrants of the included angle:
Figure BDA0001280915630000072
and 3-2, uniformly dividing the LMS-PCA space into n multiplied by n equal parts along α and theta directions respectively, calculating the spherical coordinates of each printer spectral color sample point, determining which sub-region of the LMS-PCA space each spectral color sample point is located in according to the values of α and theta, comparing the distances r from all spectral color sample points in each sub-region to the central point, reserving the spectral color sample point with the maximum r value as the boundary point of the sub-region, and finally storing all the obtained boundary points into a matrix to obtain the maximum spectral boundary GBD.
And step four, describing the specific spectral boundary LBD of the output device in the LMS-PCA space. The fourth step specifically comprises the following steps:
step 4-1: a certain spectral color point M in the spectral image is taken as a checking object. Reducing the dimension of the color point of the test spectrum to the LMS-PCA space to obtain the coordinate value of (p)1,p2,p3)。
Step 4-2, calculating a spherical coordinate α corresponding to the hue angle of the M point, and calculating a α value and a spectrum plane coordinate (C)o,p1) And the center point E finds the equal α planes of the inspection spectral color point.
And 4-3, finding out a pair of adjacent spectral color points from the maximum spectral boundary GBD sub-matrix of the first horizontal segmentation layer according to the hue angle α of the M points, wherein the hue angle of one point is more than α, and the hue angle of the other point is less than α, connecting the two adjacent spectral color points into a straight line, and solving the intersection point of the straight line and the equal α plane, wherein n intersection points are finally obtained as n horizontal layers are total, and the intersection points are linear boundaries of the printer spectral color gamut on the equal α plane.
Step 4-4: calculating p1Maximum and minimum linear boundary points on the axis. Among all GBD sub-matrices, the one with the largest is selectedp1Spectral color point of value as maximum linear boundary point pmaxSimilarly, choose to have the minimum p1Spectral color point of value as minimum linear boundary point pmin. Combined with the linear boundaries obtained in step 4-3 to form the printer specific spectral boundaries LBD.
Step five: and judging whether the color point of the spectrum of the image to be mapped needs to be mapped or not.
And judging the relation between the mapping color point and the spectral color gamut of the output device. If the mapped color point is outside the spectral color gamut of the output device, processing according to a sixth step; if the mapped color point is within the spectral gamut of the output device, no processing is performed and the seventh step is entered directly.
The method for judging whether the mapping color point is positioned in the spectral color gamut of the output device or not in the step five is as follows: coordinate point M (C) of spectral plane of point Mo,p1) And p1The center points E of the axes are connected and extend to a line intersecting the linear boundary of LBD at a pointH(C'o,p'1) If C iso<C'oAnd p ismin<p'1<pmaxIf the M points are located within the printer spectral color gamut, otherwise, the M points are located outside the printer spectral color gamut, and the next lslincellmax color gamut compression is performed, where fig. 3 is a mapping determination diagram.
Step six: the image spectral color points that exceed the printer spectral gamut are mapped according to equations 3 and 4 using the lslincellmax gamut compression algorithm in the LMS-PCA based space. The sixth step specifically comprises the following steps:
step 6-1: and mapping the color points of the spectral image to be mapped of the spectral color gamut of the super-output device into the spectral domain of the output device by utilizing the LSLINCELmax color gamut compression in the LMS-PCA-based space. Based on LSLINceLmax color gamut compression in LMS-PCA space, on the basis of a traditional LSLIN algorithm, the influence of a spectral color gamut median value of output equipment and the influence of a spectral image brightness maximum value to be mapped are added, and meanwhile, the spectral chroma and the spectral brightness of a spectral color point of a spectral image to be mapped are compressed, so that a mapping process with higher precision is completed. The spectrum mapping expression is as follows:
Figure BDA0001280915630000091
in the formula, CrFor the spectral chroma value, C, of the mapped output imager(max)Maximum spectral chroma value, C, of the spectral gamut of the output deviceoAs a spectral color point chroma value, C, of the spectral image to be mappedo(max)Is the maximum spectral chroma value of the spectral image gamut to be mapped. p is a radical ofrFor the spectral brightness value, p, of the mapped output image1The spectral image brightness value is the spectral color point brightness value of the spectral image to be mapped, the center is the brightness value of the central point of the equal hue surface spectral brightness axis corresponding to the spectral color point of the spectral image to be mapped, and the cemmax is the weighting function of the maximum brightness value Lmax and the center of the spectral image to be mapped.
Step 6-2: the cemmax expression is as follows:
Figure BDA0001280915630000092
where Lmax is the maximum brightness value of the spectral image to be mapped, cc, dd, and hh are constants, cc ranges from-15 to 15, dd ranges from 0 to 10, and hh generally ranges from 1 to 3.
Step seven: and carrying out spectral reconstruction on the color data in the LMS-PCA space obtained after mapping processing to obtain the mapped image spectral information.
Step eight: and (6) ending.
Of course, this embodiment is only used for illustrating an optimized spectral gamut mapping method of the present invention, and is not used to limit the scope of the present invention.

Claims (1)

1. An optimized spectral gamut mapping method, characterized by: comprises the following steps:
the method comprises the following steps:
step 1-1: according to the visual characteristics of human eyes, an LMS cone response spectral sensitivity curve is selected as a weighting function, high-dimensional spectral data are weighted, and when a spectral tristimulus value is converted into an LMS cone response, a Bradford transformation matrix is selected as a transformation matrix:
Figure FDA0002175989060000011
step 1-2: the LMS weighting function is constructed as equation 2:
Figure FDA0002175989060000012
wherein p is a variable coefficient and has a value ranging from-0.5 to 50, k is 1 to 10,
step 1-3: the high-dimensional spectral data is weighted by using an LMS weighting function,
step two: extracting the first three principal components of the spectrum after weighting treatment by adopting a Principal Component Analysis (PCA) method, and constructing an LMS-PCA spectrum link space (ICS);
the second step specifically comprises the following steps:
step 2-1: performing singular value decomposition on the weighted spectrum data;
step 2-2: reducing the dimension of the weighted spectral data through the eigenvector of singular value decomposition, extracting the first three principal components to construct an LMS-PCA rectangular coordinate system as a spectral link space, constructing the LMS-PCA rectangular coordinate system space by using the first three principal components obtained by dimension reduction, and setting the first three principal components as P1,P2,P3,
Step three: describing the maximum spectral boundary GBD of the output device in LMS-PCA space;
step 3-1: using a rectangular coordinate system P of LMS-PCA space1,P2,P3Converting into a spherical coordinate system α, theta, r, calculating the spherical coordinate value corresponding to the spectral color point of the output device,
step 3-2, evenly dividing the LMS-PCA space into n multiplied by n equal parts along α and theta directions, respectively, calculating the spherical coordinates of each output device spectrum color sample point, determining which sub-region of the LMS-PCA space each spectrum color sample point is located in according to the values of α and theta, comparing the distances r from all spectrum color sample points in each sub-region to the central point, reserving the spectrum color sample point with the maximum r value as the boundary point of the sub-region, finally storing all the obtained boundary points into a matrix to obtain the maximum spectrum boundary GBD,
step four: describing the specific spectral boundaries LBD of the output device in LMS-PCA space;
the fourth step comprises the following steps:
step 4-1, calculating a spherical coordinate α corresponding to the hue angle of the spectrum color point of the image to be mapped, and calculating the color point according to the value α and the spectrum plane coordinate (C)o,p1) And a center point E, finding α planes that are equal to the spectral color point, wherein the coordinates are
Figure FDA0002175989060000021
p1、p2、p3Is a rectangular coordinate system P1、P2、P3The coordinate values of (a) and (b),
step 4-2, finding out a pair of adjacent spectrum color points from the GBD submatrix of the first horizontal division layer according to the hue angle α, wherein the hue angle of one point is more than α, the hue angle of the other point is less than α, then connecting the two adjacent spectrum color points into a straight line, solving the intersection points of the straight line and the equal α plane, and finally obtaining n intersection points which are linear boundaries of the output device spectrum color gamut on the equal α plane as a total of n horizontal layers,
step 4-3: according to maximum spectral margin GBD matrix with maximum p1Value and minimum p1Spectral color sample points of values, resulting in a maximum linear boundary point p among the spectral color sample points of the output devicemaxAnd a minimum linear boundary point pminCombined with the linear boundaries obtained in step 4-2 to form specific spectral boundaries LBD of the output device,
step five: judging whether the color point of the spectrum of the image to be mapped needs to be mapped or not;
judging the relation between the mapping color point and the spectral color gamut of the output device, and if the mapping color point is positioned outside the spectral color gamut of the output device, processing according to the sixth step; if the mapped color point is within the spectral gamut of the output device, then proceed directly to step seven without processing,
the method for judging whether the mapping color point is positioned in the spectral color gamut of the output device or not in the step five is as follows:
coordinate point M (C) of spectral plane of point Mo,p1) And p1The center point E of the shaft is connected to and extends to a straight line intersecting the LBD linear boundary at point H (C'o,p'1) If C iso<C'oAnd p ismin<p'1<pmaxThen the M points are within the printer spectral gamut, otherwise they are outside the printer spectral gamut,
step six: mapping the image spectral color points beyond the spectral color gamut of the output equipment by adopting an LSLINceLmax color gamut compression algorithm;
in the sixth step, the method comprises the following steps:
step 6-1: the method comprises the following steps of utilizing LSLINCELmax color gamut compression in an LMS-PCA space to map color points of a spectral image to be mapped of a spectral color gamut of super-output equipment into a spectral domain of the output equipment, adding the influence of a median value of spectral color gamut brightness of the output equipment and the influence of a maximum value of the brightness of the spectral image to be mapped on the basis of a traditional LSLIN algorithm and simultaneously compressing the spectral chroma and the spectral brightness of the spectral color points of the spectral image to be mapped, thereby completing the mapping process with higher precision, wherein the spectral mapping expression is as follows:
Figure FDA0002175989060000031
in the formula, CrFor the spectral chroma value, C, of the mapped output imager(max)Maximum spectral chroma value, C, of the spectral gamut of the output deviceoAs a spectral color point chroma value, C, of the spectral image to be mappedo(max)Is the maximum spectral chroma value, p, of the spectral image gamut to be mappedrFor the spectral brightness value, p, of the mapped output image1Spectral color point brightness for spectral images to be mappedThe center is the brightness value of the central point of the isochromatic plane spectral brightness axis corresponding to the spectral color point of the spectral image to be mapped, the cemmax is the weighting function of the maximum brightness value Lmax and the center of the spectral image to be mapped,
step 6-2: the cemmax expression is as follows:
Figure FDA0002175989060000041
wherein Lmax is the maximum brightness value of the spectral image to be mapped, cc, dd and hh are constants, cc ranges from-15 to 15, dd ranges from 0 to 10, hh ranges from 1 to 3,
step seven: and performing spectral reconstruction based on the color data in the LMS-PCA space obtained by mapping processing.
CN201710286928.4A 2017-04-27 2017-04-27 Optimized spectrum color gamut mapping method Active CN107067444B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710286928.4A CN107067444B (en) 2017-04-27 2017-04-27 Optimized spectrum color gamut mapping method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710286928.4A CN107067444B (en) 2017-04-27 2017-04-27 Optimized spectrum color gamut mapping method

Publications (2)

Publication Number Publication Date
CN107067444A CN107067444A (en) 2017-08-18
CN107067444B true CN107067444B (en) 2020-03-06

Family

ID=59605058

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710286928.4A Active CN107067444B (en) 2017-04-27 2017-04-27 Optimized spectrum color gamut mapping method

Country Status (1)

Country Link
CN (1) CN107067444B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107680142B (en) * 2017-10-23 2019-12-24 深圳市华星光电半导体显示技术有限公司 Method for improving out-of-gamut color overlay mapping
CN113916817B (en) * 2021-10-21 2024-02-02 深圳市中科云驰环境科技有限公司 Spectrum method chromaticity online measurement method for urban living drinking water
CN115546323B (en) * 2022-08-31 2023-06-09 广东工业大学 Image compression reconstruction method based on spherical coordinate principal component analysis

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6556932B1 (en) * 2000-05-01 2003-04-29 Xerox Corporation System and method for reconstruction of spectral curves using measurements from a color sensor and a spectral measurement system model
CN103679709A (en) * 2013-11-28 2014-03-26 深圳职业技术学院 Spectral image processing method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6556932B1 (en) * 2000-05-01 2003-04-29 Xerox Corporation System and method for reconstruction of spectral curves using measurements from a color sensor and a spectral measurement system model
CN103679709A (en) * 2013-11-28 2014-03-26 深圳职业技术学院 Spectral image processing method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Ja'n Morovic等.Calculating Medium and Image Gamut Boundaries for Gamut Mapping.《Color Research and Application》.2000,第25卷(第6期),第394-401页. *
一种基于LMS加权的残差补偿光谱降维模型研究;于海琦 等;《包装工程》;20150228;第36卷(第3期);第98-102页 *
加权视觉特性的PCA空间内光谱域映射模型;刘攀 等;《光学学报》;20150630;第35卷(第6期);第0633001-(1-8)页摘要,第1小节,第3小节 *
基于人眼视觉特性的光谱降维模型研究;何颂华 等;《光谱学与光谱分析》;20150630;第35卷(第6期);第1459-1463页 *
基于修改的LMS曲线加权残差补偿光谱降维法;孙叶维 等;《包装工程》;20160531;第37卷(第9期);第114-119页摘要,第1小节 *
彩色复制中的色域映射方法;王淑真 等;《北京印刷学院学报》;20060630;第14卷(第3期);第4-6页 *

Also Published As

Publication number Publication date
CN107067444A (en) 2017-08-18

Similar Documents

Publication Publication Date Title
CN106471567B (en) Color calibration system and method
CN107067444B (en) Optimized spectrum color gamut mapping method
CN104599636B (en) LED display bright chroma bearing calibration and bright chroma correction coefficient generating means
CN113222823B (en) Hyperspectral image super-resolution method based on mixed attention network fusion
US9947258B2 (en) Image processing method and image processing apparatus
CN103528968A (en) Reflectance spectrum reconstruction method based on iterative threshold method
CN105761234A (en) Structure sparse representation-based remote sensing image fusion method
CN103268596B (en) A kind of method for reducing picture noise and making color be near the mark
CN105303542A (en) Gradient weighted-based adaptive SFIM image fusion algorithm
CN106709888B (en) A kind of high dynamic range images production method based on human vision model
CN108038852A (en) A kind of Quality Measures for Image Fusion represented based on joint sparse
CN109859153B (en) Multispectral image fusion method based on adaptive spectrum-spatial gradient sparse regularization
CN106296749A (en) RGB D image eigen decomposition method based on L1 norm constraint
CN108680535A (en) Based on the spectral reflectance recovery method for improving R matrixes
CN111563866B (en) Multisource remote sensing image fusion method
CN111010557A (en) White point range determining method, white balance correcting method and device, and storage medium
CN112862906B (en) Color space conversion method based on neural network
CN104680150A (en) Face identification method on basis of duplex multi-kernel discriminant analysis of color features
US9384534B2 (en) Method and system for establishing model based on virtual algorithm
Zhang et al. A method for selecting training samples based on camera response
CN104794690A (en) Three-dimensional image color-to-gray correcting method based on SIFT (scale invariant feature transform) matching
Luo Classification precision analysis on different fusion algorithm for ETM+ remote sensing image
Tsagaris et al. A measure for evaluation of the information content in color images
CN106210554B (en) A method of synthesis high dynamic range images
Petrović et al. Study of objective evaluation of natural colour image fusion

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant