CN103530494A - Subspace search algorithm in reverse CYNSN model - Google Patents

Subspace search algorithm in reverse CYNSN model Download PDF

Info

Publication number
CN103530494A
CN103530494A CN201210595446.4A CN201210595446A CN103530494A CN 103530494 A CN103530494 A CN 103530494A CN 201210595446 A CN201210595446 A CN 201210595446A CN 103530494 A CN103530494 A CN 103530494A
Authority
CN
China
Prior art keywords
subspace
spectrum
ink
candidate
summit
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.)
Pending
Application number
CN201210595446.4A
Other languages
Chinese (zh)
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201210595446.4A priority Critical patent/CN103530494A/en
Publication of CN103530494A publication Critical patent/CN103530494A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a subspace search algorithm in a reverse CYNSN model. By calculating RRMS values between all top point spectrums and target spectrums, the top point spectrum with the minimum RRMS value is selected, all subspaces including the top point are found out afterwards and serve as candidate subspaces, and then candidate ink formulas corresponding to all the candidate subspaces are calculated through a linear iterative algorithm, spectrums corresponding to the ink formulas are predicted, and finally the candidate subspace corresponding to the minimum RRMS value between the top point spectrums and the target spectrums is selected as the optimum subspace. The algorithm is high in search speed, precision and practical value, and the optimum subspace can be acquired after calculation through few iterations.

Description

Subspace search algorithm in a kind of reverse CYNSN model
Technical field
The present invention relates to a kind of subspace search algorithm, specifically relate to the subspace search algorithm in a kind of reverse CYNSN model.
Background technology
Along with scientific and technical development, printer has obtained increasingly extensive application in daily life work, when being greatly convenient for people to Working Life, the series of problems such as cross-color have also been brought, as do not mated with the color on display through the color of printer output etc.
The characterization of equipment is to solve one of problem of the primary solution of cross-color problem.Apparatus characteristic refers to by specific mathematical model and training sample sets up the correspondent transform relation between device independent look and equipment related colour.Device independent look refers generally to the color characterizing in CIE standard color space, as CIE XYZ tristimulus values, CIELAB colorimetric parameter etc.; Equipment related colour is that the color characteristics by respective digital vision facilities determines, for printer, its equipment related colour is often referred to the consumption of each ink, i.e. ink set.Characterization model generally comprises forward and reverse model.Forward characterization model is set up from equipment related colour to the contact device independent look, and opposite feature model refers to by device independent look and predicts corresponding equipment related colour.In printer characteristic model investigation in the past, often select CIE XYZ tristimulus values or CIELAB value as device independent look, reason is to calculate relatively simple, but because CIEXYZ or CIELAB value do not comprise the information of all colours, can cause occurring the problem of metamerism, be that master sample and reproduction sample its color under a kind of lighting condition are mated mutually, but under another kind of lighting condition unmatched phenomenon.Because spectroscopic data has recorded all information of object color, therefore the characterization model based on spectrum can solve the problem of metamerism well.
CYNSN (Cellular Yule-Nielsen Spectral Neugebauer, the spectrum Neugebauer that subspace formula Yule-Nielsen revises) model is the printer spectrum characterization model of commonly using the most, it is the improved model of basic Neugebauer model, comprises forward and reverse model.The opposite feature model of this model (based target Forecast of Spectra ink set) always is the focus of correlative study.In CYNSN model, whole ink space is divided into a plurality of subspaces, so in reverse CYNSN model, for a certain target optical spectrum, first need to find the subspace at target optical spectrum place, then according to linear iterative algorithm, this target optical spectrum is recorded to ink set in advance.Up to the present, seldom there is simple, effective subspace search algorithm to find fast the affiliated subspace of optimum target optical spectrum.
Summary of the invention
The present invention, in order to solve the problem described in background technology, discloses the subspace search algorithm in a kind of reverse CYNSN model.Its concrete steps are as follows:
1) choose n ink of printer, each ink is divided into m section from 0 (without ink) to 1 (maximum quantity of ink), altogether whole ink space is divided into m nsub spaces, each subspace comprises 2 nindividual summit, (m+1) contained in whole ink space nindividual summit; In general, for obtaining higher model accuracy, the value of m should be more than or equal to 6;
2) print all (m+1) nindividual summit is as training sample, adopt afterwards spectrophotometer to measure to obtain corresponding spectrum to all training samples, according to the establishment of spectrum forward CYNSN model of all training samples that record, the foundation of forward CYNSN model can be referring to journal of writings (the Optimizationof the spectral Neugebauer model for printer characterization.Journal of ElectronicImaging of R.Balasubramanian, 1999, the 8th the 2nd phase of volume, 156-166);
3) for a certain target optical spectrum, the RRMS of calculating and each summit spectrum, computing formula is as follows:
RRMS = Σ λ = λ 1 λ 2 ( R 1 ( λ ) - R 2 ( λ ) ) 2 N
In formula, R 1(λ) and R 2(λ) be the spectrum on target sample and a certain summit, λ 1 and λ 2 are minimum value and the maximal value within the scope of measure spectrum, and N is the number of samples within the scope of measure spectrum.
4) select that summit minimum with target optical spectrum RRMS value, find out afterwards all subspaces that comprise this summit, be called candidate subspace, be assumed to t candidate subspace;
5) in each candidate subspace, with linear iterative algorithm, calculate the ink set of target optical spectrum, obtain altogether t candidate's ink set, linear iterative algorithm adopts UG algorithm (P.Urban and R.R.Grigat, Spectralbased color separation using linear regression iteration.Color Research Application, 2006, the 31 phase the 3rd volumes, 229-238);
6) for each candidate's ink set, utilize its corresponding spectrum of forward CYNSN model prediction, calculate the RRMS value of this prediction spectrum and target optical spectrum simultaneously;
7) finally selecting minimum corresponding that candidate subspace of RRMS value in all prediction spectrum is optimal subspace, the final ink set that the ink set calculating in this candidate subspace is target optical spectrum.
The present invention is by selecting to find all subspaces that comprise this summit with that summit of target optical spectrum RRMS minimum, thereby greatly dwindle the hunting zone of subspace, accelerate the speed of algorithm, then utilize linear iterative algorithm to calculate candidate's ink set corresponding to all candidates subspace, and predict the spectrum that it is corresponding, finally select corresponding that candidate subspace of RRMS value minimum between prediction spectrum and target optical spectrum as optimal subspace.The method search speed is fast, precision is high, only needs less iterations just can calculate optimal subspace, can obviously improve the solution efficiency of the more reverse CYNSN model in subspace.
Accompanying drawing explanation
Fig. 1 is the subspace search algorithm flow chart in a kind of reverse CYNSN model;
Fig. 2 is the target optical spectrum in embodiment;
Fig. 3 is the RRMS value of prediction spectrum and the target optical spectrum of 8 candidate's ink sets in embodiment;
Fig. 4 is prediction spectrum that in embodiment, final ink set is corresponding and the comparison of target optical spectrum.
Embodiment
Traditional three black CMY (green grass or young crops, product, the Huang) printer of take is example, and the embodiment of the subspace search algorithm in above-mentioned a kind of reverse CYNSN model is set forth.As shown in Figure 1, its concrete steps are as follows:
1) so choose 3 inks of printer, each ink is evenly divided into 7 sections from 0 (without ink) to 1 (maximum quantity of ink), altogether whole ink space is divided into 7 3=343 sub spaces, each subspace comprises 2 3=8 summits, whole ink space contains 8 3=512 summits;
2) print all 512 summits as training sample, adopt afterwards spectrophotometer to measure to obtain corresponding spectrum to all training samples, according to the establishment of spectrum forward CYNSN model (Yule-Nielsen modified value is made as 2.5) of all training samples that record; Wherein, spectrophotometric measurement range is made as 400nm to 700nm, take 10nm as interval;
3) for a certain target optical spectrum, see Fig. 2, the RRMS of calculating and each summit spectrum, computing formula is as follows:
RRMS = Σ λ = λ 1 λ 2 ( R 1 ( λ ) - R 2 ( λ ) ) 2 N
In formula, R 1(λ) and R 2(λ) be the spectrum on target sample and a certain summit, λ 1 and λ 2 are minimum value 400nm and the maximal value 700nm within the scope of measure spectrum, and N=31 is the number of samples within the scope of measure spectrum.
4) select that summit minimum with target optical spectrum RRMS value, find out afterwards all subspaces that comprise this summit, be called candidate subspace, in this embodiment t=8 candidate subspace; Wherein, the value of t is relevant with the position in whole ink space on this summit, t=1 or 4 or 8;
5) in each candidate subspace, with linear iterative algorithm, calculate the ink set of target optical spectrum, obtain altogether t candidate's ink set, linear iterative algorithm adopts UG algorithm;
6) for each candidate's ink set, utilize its corresponding spectrum of forward CYNSN model prediction, calculate the RRMS value of this prediction spectrum and target optical spectrum simultaneously, see Fig. 3;
7) finally selecting minimum corresponding that candidate subspace of RRMS value in all prediction spectrum is optimal subspace, the final ink set that the ink set calculating in this candidate subspace is target optical spectrum, this final ink set is C=0.369, M=0.098, Y=0.122, the prediction spectrum of its correspondence and target optical spectrum are as shown in Figure 4.

Claims (6)

1. a subspace search algorithm in reverse CYNSN model, is characterized in that comprising the following steps:
1) choose n ink of printer, each ink is divided into m section from 0 (without ink) to 1 (maximum quantity of ink), altogether whole ink space is divided into m nsub spaces, each subspace comprises 2 nindividual summit, (m+1) contained in whole ink space nindividual summit;
2) print all (m+1) nindividual summit, as training sample, adopts spectrophotometer to measure to obtain corresponding spectrum to all training samples, according to the establishment of spectrum forward CYNSN model of all training samples that record afterwards;
3) for a certain target optical spectrum, the RRMS of calculating and each summit spectrum, computing formula is as follows:
RRMS = Σ λ = λ 1 λ 2 ( R 1 ( λ ) - R 2 ( λ ) ) 2 N
In formula, R 1(λ) and R 2(λ) be the spectrum on target sample and a certain summit, λ 1 and λ 2 are minimum value and the maximal value within the scope of measure spectrum, and N is the number of samples within the scope of measure spectrum.
4) select that summit minimum with target optical spectrum RRMS value, find out afterwards all subspaces that comprise this summit, be called candidate subspace, be assumed to t candidate subspace;
5) in each candidate subspace, with linear iterative algorithm, calculate the ink set of target optical spectrum, obtain altogether t candidate's ink set, linear iterative algorithm adopts UG algorithm;
6) for each candidate's ink set, utilize its corresponding spectrum of forward CYNSN model prediction, calculate the RRMS value of this prediction spectrum and target optical spectrum simultaneously;
7) finally selecting minimum corresponding that candidate subspace of RRMS value in all prediction spectrum is optimal subspace, the final ink set that the ink set calculating in this candidate subspace is target optical spectrum.
2. the subspace search algorithm in a kind of reverse CYNSN model according to claim 1, it is characterized in that described step 1) in choose n ink of printer, each ink is divided into m section from 0 (without ink) to 1 (maximum quantity of ink), altogether whole ink space is divided into m nsub spaces.In general, for obtaining higher model accuracy, the value of m should be more than or equal to 6.
3. the subspace search algorithm in a kind of reverse CYNSN model according to claim 1,, it is characterized in that described step 2) in print all (m+1) nindividual summit, as training sample, adopts spectrophotometer to measure to obtain corresponding spectrum to all training samples, according to the establishment of spectrum forward CYNSN model of all training samples that record afterwards.The foundation of forward CYNSN model can adopt the method for R.Balasubramanian.
4. the subspace search algorithm in a kind of reverse CYNSN model according to claim 1, it is characterized in that described step 4) in select that summit with target optical spectrum RRMS value minimum, find out afterwards all subspaces that comprise this summit, be called candidate subspace.
5. the subspace search algorithm in a kind of reverse CYNSN model according to claim 1, it is characterized in that described step 5) in each candidate subspace, with linear iterative algorithm, calculate the ink set of target optical spectrum, linear iterative algorithm can adopt UG algorithm.
6. the subspace search algorithm in a kind of reverse CYNSN model according to claim 1, it is characterized in that described step 7) in final to select corresponding that candidate subspace of RRMS value minimum in all prediction spectrum be optimal subspace, the final ink set that the ink set calculating in this candidate subspace is target optical spectrum.
CN201210595446.4A 2012-12-10 2012-12-10 Subspace search algorithm in reverse CYNSN model Pending CN103530494A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210595446.4A CN103530494A (en) 2012-12-10 2012-12-10 Subspace search algorithm in reverse CYNSN model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210595446.4A CN103530494A (en) 2012-12-10 2012-12-10 Subspace search algorithm in reverse CYNSN model

Publications (1)

Publication Number Publication Date
CN103530494A true CN103530494A (en) 2014-01-22

Family

ID=49932499

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210595446.4A Pending CN103530494A (en) 2012-12-10 2012-12-10 Subspace search algorithm in reverse CYNSN model

Country Status (1)

Country Link
CN (1) CN103530494A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103862858A (en) * 2014-03-11 2014-06-18 武汉大学 Color separation method of spectrum of printing system with multicolor printing

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101042832A (en) * 2007-04-27 2007-09-26 浙江大学 Method for building LCD colors characterizing segment and space division model
CN102180003A (en) * 2011-02-22 2011-09-14 浙江大学 Method of correcting printer spectrum characterization model

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101042832A (en) * 2007-04-27 2007-09-26 浙江大学 Method for building LCD colors characterizing segment and space division model
CN102180003A (en) * 2011-02-22 2011-09-14 浙江大学 Method of correcting printer spectrum characterization model

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PHILIPP URBAN,ET AL.: "Spectral-Based Color Separation Using Linear Regression Iteration", 《COLOR RESEARCH AND APPLICATION》, vol. 31, no. 3, 30 June 2006 (2006-06-30), pages 229 - 238 *
RAJA BALASUBRAMANIAN.: "Optimization of the spectral Neugebauer model for printer characterization", 《JOURNAL OF ELECTRONIC IMAGING》, vol. 8, no. 2, 30 April 1999 (1999-04-30), pages 156 - 166, XP 000825806, DOI: doi:10.1117/1.482694 *
郭晋一.: "基于光谱匹配法Neugebauer方程的彩色打印系统特征化模型研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, no. 07, 15 July 2011 (2011-07-15), pages 138 - 646 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103862858A (en) * 2014-03-11 2014-06-18 武汉大学 Color separation method of spectrum of printing system with multicolor printing
CN103862858B (en) * 2014-03-11 2015-12-30 武汉大学 A kind of multi-color printing print system spectrum color separation method

Similar Documents

Publication Publication Date Title
CN104217409B (en) A kind of image color correction method based on simulated annealing optimization algorithm
CN100483483C (en) Method for building LCD colors characterizing segment and space division model
CN201464052U (en) Color CCD-based multi-spectral color radiation temperature measurement device
CN103747161B (en) The hard proofing method of a kind of high-fidelity spot color numeral
CN104410850A (en) Colorful digital image chrominance correction method and system
CN101476942B (en) Radiation temperature field measurement apparatus and method based on colorful CCD
US20150339552A1 (en) Ink set optimization
CN107424197B (en) Method for realizing cross-media color reproduction based on spectral domain mapping
CN102799895A (en) Offset printing ink color matching method based on least square support vector machine
JP2011071924A (en) Image processing apparatus and image processing method
CN105069234A (en) Spectrum dimensionality reduction method and system based on visual perception feature
CN112484856A (en) Method for obtaining high-precision chromaticity and spectrum image
CN108872156B (en) Method and device for predicting ink component proportion based on reciprocal of spectral reflectance
CN102180003B (en) Method of correcting printer spectrum characterization model
CN106153193B (en) A kind of method for obtaining spectral reflectance using the double light source responses of multispectral camera
CN103530494A (en) Subspace search algorithm in reverse CYNSN model
CN103973980A (en) Method for acquiring multispectral images based on digital camera
CN107452042A (en) A kind of halftoning apparatus spectral characteristic model optimization method and system
CN103862858A (en) Color separation method of spectrum of printing system with multicolor printing
Sun et al. Modified spectral Neugebauer model for printer characterization
CN103544334B (en) Many ink printer formulation selection algorithms based on CII and aberration
Hu et al. Study of color rendering evaluation method of light sources for printing matter
CN103057292B (en) Printer calibration algorithm based on characteristics of human eye
CN107766896A (en) A kind of spectrum dimension reduction method based on form and aspect cluster
CN105628190A (en) Method and device for measuring optical radiation on the basis of light-filtering units

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140122