CN103057292B - Printer calibration algorithm based on characteristics of human eye - Google Patents
Printer calibration algorithm based on characteristics of human eye Download PDFInfo
- Publication number
- CN103057292B CN103057292B CN201310029840.6A CN201310029840A CN103057292B CN 103057292 B CN103057292 B CN 103057292B CN 201310029840 A CN201310029840 A CN 201310029840A CN 103057292 B CN103057292 B CN 103057292B
- Authority
- CN
- China
- Prior art keywords
- calibration
- printer
- human eye
- spectrum
- wavelength
- 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.)
- Expired - Fee Related
Links
Landscapes
- Color Image Communication Systems (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种打印机的定标算法,具体是涉及一种基于人眼特性的打印机定标算法。The invention relates to a printer calibration algorithm, in particular to a printer calibration algorithm based on the characteristics of human eyes.
背景技术Background technique
随着工业技术的发展,打印机已经广泛地应用于人们的日常生活中。与此同时,大量科学工作者对打印机的墨水颜色管理进行了相应的研究。目前,相关的一个研究热点就是打印机的光谱特征化模型,该模型可以从墨水配方(即各墨水用量)预测其在某一纸张上打印之后的光谱数据。现有的打印机光谱特征化模型主要有三维查找表法、Yule-Clapper模型及其改进模型、Neugebauer模型及其改进模型等。With the development of industrial technology, printers have been widely used in people's daily life. At the same time, a large number of scientific workers have carried out corresponding research on ink color management of printers. At present, a related research hotspot is the spectral characterization model of the printer, which can predict the spectral data after printing on a certain paper from the ink formula (ie, the amount of each ink). The existing printer spectral characterization models mainly include three-dimensional look-up table method, Yule-Clapper model and its improved model, Neugebauer model and its improved model, etc.
就当前的技术而言,要建立一高精度的打印机光谱特征化模型至少需要打印并测量几千个样本,故建立一光谱模型是一个比较费时、费力的工程。但是,最终模型的预测精度仍受到多个因素的制约,如打印机的重复性、纸张的均匀性、环境温度、环境湿度、墨盒的更换等。为了解决这一方面的问题,就需要进行打印机的定标,即通过在现有打印条件下打印一些样本建立与已经建立的模型(下文简称“原模型”)之间的联系,从而实现原模型在现有打印条件下的高精度预测。虽然目前也有一些关于打印机定标方面的研究,但是所有的打印机定标算法均没有考虑人眼的特性,而且打印出的颜色是给人评判或欣赏的,因此,提出一基于人眼特性的打印机定标算法是非常有必要的。As far as the current technology is concerned, at least thousands of samples need to be printed and measured to establish a high-precision printer spectral characterization model, so establishing a spectral model is a relatively time-consuming and laborious project. However, the prediction accuracy of the final model is still constrained by multiple factors, such as printer repeatability, paper uniformity, ambient temperature, ambient humidity, ink cartridge replacement, etc. In order to solve this problem, it is necessary to calibrate the printer, that is, by printing some samples under the existing printing conditions to establish the connection with the established model (hereinafter referred to as "the original model"), so as to realize the original model. High-precision predictions under existing printing conditions. Although there are some studies on printer calibration, all printer calibration algorithms do not consider the characteristics of human eyes, and the printed colors are judged or appreciated by people. Therefore, a printer based on the characteristics of human eyes is proposed. Calibration algorithm is very necessary.
发明内容Contents of the invention
本发明为了解决背景技术中所述的问题,公开了一种基于人眼特性的打印机定标算法。本方法结合多项式和代表人眼特性的CIE(国际照明委员会)标准色度观察者的光谱三刺激值函数实现打印机的定标,其具体步骤如下:In order to solve the problems described in the background technology, the present invention discloses a printer calibration algorithm based on the characteristics of human eyes. This method realizes the calibration of the printer by combining the polynomial and the spectral tristimulus value function of the CIE (International Commission on Illumination) standard chromaticity observer representing the characteristics of the human eye. The specific steps are as follows:
1)选取n个基本覆盖各个色区的定标墨水配方,一般n为30-50个;1) Select n calibration ink formulas that basically cover each color area, generally n is 30-50;
2)利用建立原模型时所采用的打印机和纸张输出定标墨水配方的打印样本,即定标样本,并将定标样本在暗室中放置一段时间以确保墨水充分干燥;2) Utilize the printer and paper used to build the original model to output the print sample of the calibration ink formula, that is, the calibration sample, and place the calibration sample in the darkroom for a period of time to ensure that the ink is fully dry;
3)采用一分光光度计对所有定标样本进行测量,测得其对应的测量光谱Rm(λ),其中下标m表示测量,λ表示波长;3) A spectrophotometer is used to measure all calibration samples, and the corresponding measurement spectrum R m (λ) is obtained, wherein the subscript m indicates measurement, and λ indicates wavelength;
4)利用原模型对所有的定标墨水配方进行预测得预测光谱Rp(λ),其中,下标p表示预测;4) Use the original model to predict all calibration ink formulations to obtain the predicted spectrum R p (λ), where the subscript p represents prediction;
5)采用CIE1964标准色度观察者的三刺激值函数x(λ)、y(λ)和z(λ)在各波长之和S(λ)作为人眼特性函数;5) The sum S(λ) of the three-stimulus value functions x(λ), y(λ) and z(λ) of the CIE1964 standard chromaticity observer at each wavelength is used as the human eye characteristic function;
6)对于波长λ,假定Rm(λ)、Rp(λ)和S(λ)满足下列二次多项式公式:6) For wavelength λ, it is assumed that R m (λ), R p (λ) and S(λ) satisfy the following quadratic polynomial formula:
S(λ)×Rm(λ)=a(λ)×[S(λ)×Rp(λ)]2+b(λ)×S(λ)×Rp(λ)+c(λ)S(λ)×R m (λ)=a(λ)×[S(λ)×R p (λ)] 2 +b(λ)×S(λ)×R p (λ)+c(λ)
其中,a(λ)、b(λ)和c(λ)为定标系数;Among them, a(λ), b(λ) and c(λ) are calibration coefficients;
7)根据最小二乘法,计算最终的定标系数a(λ)、b(λ)和c(λ);7) Calculate the final calibration coefficients a(λ), b(λ) and c(λ) according to the least square method;
8)对于任一墨水配方,可先由原模型预测其光谱R’p(λ),然后根据步骤6)中的公式和步骤7)中计算的定标系数a(λ)、b(λ)和c(λ),即可求得最终定标后的光谱R’m(λ),其值为8) For any ink formula, its spectrum R' p (λ) can be predicted by the original model first, and then according to the formula in step 6) and the calibration coefficients a(λ) and b(λ) calculated in step 7) and c(λ), the final calibrated spectrum R' m (λ) can be obtained, and its value is
R’m(λ)={a(λ)×[S(λ)×R’p(λ)]2+b(λ)×S(λ)×R’p(λ)+c(λ)}/S(λ)。R' m (λ)={a(λ)×[S(λ)×R' p (λ)] 2 +b(λ)×S(λ)×R' p (λ)+c(λ)} /S(λ).
本发明通过打印和测量定标样本得到的测量光谱,以及原模型的预测光谱,结合多项式和人眼特性实现打印机的定标。该发明应用简单、精度高、具有很高的实用价值。The invention realizes the calibration of the printer by printing and measuring the measurement spectrum obtained from the calibration sample and the prediction spectrum of the original model, and combining the polynomial and the characteristics of the human eye. The invention has the advantages of simple application, high precision and high practical value.
附图说明Description of drawings
图1是基于人眼特性的打印机定标算法流程图;Figure 1 is a flow chart of the printer calibration algorithm based on the characteristics of the human eye;
图2是CIE1964标准色度观察者的三刺激值函数x(λ)、y(λ)和z(λ);Fig. 2 is the tristimulus value function x(λ), y(λ) and z(λ) of CIE1964 standard chromaticity observer;
图3是CIE1964标准色度观察者的三刺激值函数之和S(λ);Fig. 3 is the sum S(λ) of the tristimulus value function of CIE1964 standard chromaticity observer;
图4是实施例中最终计算的定标系数a(λ)、b(λ)和c(λ);Fig. 4 is the scaling factor a (λ), b (λ) and c (λ) of final calculation in the embodiment;
具体实施方式Detailed ways
以一台HP Designjet1055CM四墨CMYK(青品黄黑)四墨打印机和一基于该打印机的已经建立的Yule-Clapper光谱模型(下文简称“原YC模型”)为例,对上述基于人眼特性的打印机定标算法的具体实施方式进行阐述。如图1所示,其具体步骤如下:Taking an HP Designjet1055CM four-ink CMYK (green magenta yellow-black) four-ink printer and an established Yule-Clapper spectral model based on the printer (hereinafter referred to as "the original YC model") as an example, the above-mentioned analysis based on the characteristics of the human eye The specific implementation of the printer calibration algorithm will be described. As shown in Figure 1, the specific steps are as follows:
1)选取40个基本覆盖各个色区的定标墨水配方,即CMYK各墨水用量的组合;1) Select 40 calibration ink formulas that basically cover each color area, that is, the combination of CMYK ink consumption;
2)利用建立原模型时所采用的打印机和纸张输出定标墨水配方的打印样本,即定标样本,并将定标样本在暗室中放置1小时以确保墨水充分干燥;2) Use the printer and paper used to build the original model to output the printed sample of the calibration ink formula, that is, the calibration sample, and place the calibration sample in the dark room for 1 hour to ensure that the ink is fully dry;
3)采用Datacolor650分光光度计对所有40个定标样本进行测量,测得其对应的测量光谱Rm(λ),其中下标m表示测量,λ表示波长,测量范围为400nm到700nm以10nm为间隔;3) Use Datacolor650 spectrophotometer to measure all 40 calibration samples, and measure the corresponding measurement spectrum R m (λ), where the subscript m means measurement, λ means wavelength, and the measurement range is from 400nm to 700nm with 10nm as the interval;
4)利用原YC模型对所有的定标墨水配方进行预测得预测光谱Rp(λ),其中,下标p表示预测;4) Use the original YC model to predict all calibration ink formulations to obtain the predicted spectrum R p (λ), where the subscript p represents prediction;
5)采用CIE1964标准色度观察者的三刺激值函数x(λ)、y(λ)和z(λ)在各波长之和S(λ)作为人眼特性函数。其中,CIE1964标准色度观察者的三刺激值函数x(λ)、y(λ)和z(λ)如图2所示,CIE1964标准色度观察者的三刺激值函数之和S(λ)如图3所示;5) The sum S(λ) of the three-stimulus value functions x(λ), y(λ) and z(λ) at each wavelength of the CIE1964 standard chromaticity observer is used as the human eye characteristic function. Among them, the tristimulus value functions x(λ), y(λ) and z(λ) of the CIE1964 standard chromaticity observer are shown in Figure 2, and the sum of the tristimulus value functions of the CIE1964 standard chromaticity observer S(λ) As shown in Figure 3;
6)对于波长λ,假定Rm(λ)、Rp(λ)和S(λ)满足下列二次多项式公式:6) For wavelength λ, it is assumed that R m (λ), R p (λ) and S(λ) satisfy the following quadratic polynomial formula:
S(λ)×Rm(λ)=a(λ)×[S(λ)×Rp(λ)]2+b(λ)×S(λ)×Rp(λ)+c(λ)S(λ)×R m (λ)=a(λ)×[S(λ)×R p (λ)] 2 +b(λ)×S(λ)×R p (λ)+c(λ)
其中,a(λ)、b(λ)和c(λ)为定标系数;Among them, a(λ), b(λ) and c(λ) are calibration coefficients;
7)根据最小二乘法,计算最终的定标系数a(λ)、b(λ)和c(λ),如图4所示;7) Calculate the final calibration coefficients a (λ), b (λ) and c (λ) according to the least squares method, as shown in Figure 4;
8)对于任一墨水配方,可先由原模型预测其光谱R’p(λ),然后根据步骤6)中的公式和步骤7)中计算的定标系数a(λ)、b(λ)和c(λ),即可求得最终定标后的光谱R’m(λ),其值为8) For any ink formula, its spectrum R' p (λ) can be predicted by the original model first, and then according to the formula in step 6) and the calibration coefficients a(λ) and b(λ) calculated in step 7) and c(λ), the final calibrated spectrum R' m (λ) can be obtained, and its value is
R’m(λ)={a(λ)×[S(λ)×R’p(λ)]2+b(λ)×S(λ)×R’p(λ)+c(λ)}/S(λ)。R' m (λ)={a(λ)×[S(λ)×R' p (λ)] 2 +b(λ)×S(λ)×R' p (λ)+c(λ)} /S(λ).
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310029840.6A CN103057292B (en) | 2013-01-18 | 2013-01-18 | Printer calibration algorithm based on characteristics of human eye |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310029840.6A CN103057292B (en) | 2013-01-18 | 2013-01-18 | Printer calibration algorithm based on characteristics of human eye |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103057292A CN103057292A (en) | 2013-04-24 |
| CN103057292B true CN103057292B (en) | 2014-11-05 |
Family
ID=48100301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310029840.6A Expired - Fee Related CN103057292B (en) | 2013-01-18 | 2013-01-18 | Printer calibration algorithm based on characteristics of human eye |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103057292B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103528686A (en) * | 2013-09-24 | 2014-01-22 | 岑夏凤 | Formula optimization algorithm for LED (light emitting diode) illumination system with adjustable spectrum |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4479718A (en) * | 1982-06-17 | 1984-10-30 | E. I. Du Pont De Nemours And Company | Three direction measurements for characterization of a surface containing metallic particles |
| JP3615390B2 (en) * | 1998-05-12 | 2005-02-02 | 株式会社ニレコ | Measured value analysis method of on-line spectrometer |
| US7586656B2 (en) * | 2006-06-06 | 2009-09-08 | E. I. Du Pont De Nemours And Company | Process for producing optimised printing forms |
| JP4408915B2 (en) * | 2006-08-08 | 2010-02-03 | キヤノン株式会社 | Profile creation method, profile creation program, and printer |
| US8401289B2 (en) * | 2009-07-31 | 2013-03-19 | Eastman Kodak Company | Method for matching colors between two systems |
| CN102180003B (en) * | 2011-02-22 | 2013-04-24 | 浙江大学 | Method of correcting printer spectrum characterization model |
-
2013
- 2013-01-18 CN CN201310029840.6A patent/CN103057292B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN103057292A (en) | 2013-04-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104897374B (en) | A kind of color calibration method of camera module | |
| CN103954362B (en) | A kind of digital method for measuring color based on imaging device | |
| CN102082893B (en) | Digital proof color correcting method and digital proof color correcting device | |
| JP2005508125A (en) | Spectral color reproduction with 6-color output | |
| JP2001209788A (en) | Online model predicting method of dynamically fluctuating color reproducing device and color printer | |
| CN105522822B (en) | Method for calculating spot color database | |
| Li et al. | Comparison of real colour gamuts using a new reflectance database | |
| CN102180003B (en) | Method of correcting printer spectrum characterization model | |
| ES2469717T3 (en) | Method and apparatus for calculating metameric colors | |
| CN103057292B (en) | Printer calibration algorithm based on characteristics of human eye | |
| CN102238297B (en) | Method and system for generating international color consortium profile file | |
| CN105681624A (en) | Spot color overprinted color prediction method based on spectrum | |
| JP5736382B2 (en) | Optimized printer calibration method | |
| CN107389194A (en) | Color high-fidelity digital imaging system under opening luminous environment based on colour atla | |
| Chou et al. | Methods for designing characterisation targets for digital cameras | |
| US8253996B2 (en) | System, color image producing device, color measurement device and color measurement method | |
| US20070070367A1 (en) | Linearization and characterization process of an image producing device | |
| CN113119447B (en) | Method for color space conversion of color 3D printing | |
| CN106739484B (en) | A method of prediction printed matter drying process color change | |
| CN103540187B (en) | Method for improving formula of printer ink | |
| WO2022070773A1 (en) | Image analysis method, image analysis device, program, and recording medium | |
| JP6155766B2 (en) | Print reproduction color prediction method | |
| CN114013033B (en) | A method for controlling the color of photosensitive resin type full-color 3D printing | |
| Nussbaum et al. | Regression based characterization of color measurement instruments in printing applications | |
| CN103538385B (en) | Based on the method for the how black sample spectrum of the super black amount restriction of cell element formula Neugebauer spectral model prediction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: CIXI FUHAI TOWN SIWEI INDUSTRIAL PRODUCT DESIGN ST Free format text: FORMER OWNER: CEN XIAFENG Effective date: 20140924 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20140924 Address after: 315332 Cixi, China attached to the sea in the town of East China Sea Village view attached Highway No. 1285 Applicant after: Cixi Fuhai Siwei Industry Product Work Design Office Address before: Flower Village Fuhai Town Zhejiang city Cixi province 315332 triangle field Applicant before: Cen Xiafeng |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| CB02 | Change of applicant information |
Inventor after: Liao Dazhen Inventor before: Cen Xiafeng |
|
| COR | Change of bibliographic data | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20161214 Address after: 300400 8-110, Linglan Garden, Vanke Garden, Xinyibai Avenue, Beichen District, Tianjin City# Patentee after: Liao Dazhen Address before: 315332 Cixi, China attached to the sea in the town of East China Sea Village view attached Highway No. 1285 Patentee before: Cixi Fuhai Siwei Industry Product Work Design Office |
|
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20190118 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |