CN105744267B - 基于齐次多项式数码相机参数可变的获取三刺激值方法 - Google Patents
基于齐次多项式数码相机参数可变的获取三刺激值方法 Download PDFInfo
- Publication number
- CN105744267B CN105744267B CN201610236691.4A CN201610236691A CN105744267B CN 105744267 B CN105744267 B CN 105744267B CN 201610236691 A CN201610236691 A CN 201610236691A CN 105744267 B CN105744267 B CN 105744267B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- rgb
- camera
- quantic
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N17/00—Diagnosis, testing or measuring for television systems or their details
- H04N17/002—Diagnosis, testing or measuring for television systems or their details for television cameras
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Spectrometry And Color Measurement (AREA)
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
- Color Image Communication Systems (AREA)
Abstract
本发明公开了一种基于齐次多项式数码相机参数可变的获取三刺激值方法。该方法将数码相机拍摄设置参数引入作为输入变量，建立实测样本和训练样本不同参数设置下相机响应值关联，基于一阶齐次多项式模型，将等效RGB值转换为实测样本的XYZ三刺激值。本发明克服了传统方法中要求相机必须固定拍摄参数进行训练和实测的缺陷，使得数码相机可以在固定一种拍摄参数下训练，而在实际拍摄测量时则可以自由变化拍摄参数，从而为数码相机作为一种成像式测色仪器提供了便利。
Description
技术领域
本发明涉及颜色管理系统中的基于数码相机成像式获取颜色三刺激值的方法，尤其是在数码相机的设置参数变化的情况下训练和实际测量的方法。
背景技术
目前，相对传统的测色设备，使用商用的数码相机作为一种成像式的测色仪器具有成本低，操作简便等优点，因而在定量捕获场景色度信息的领域得到了越来越多的应用。然而，相机的RGB空间是设备相关的空间，相同的场景使用不同的相机会得到不同的RGB值，因此需要建立相机RGB颜色空间到CIEXYZ标准颜色空间的映射关系，这就是相机的色度特征化。为保证得到原始相机捕获数据，一般基于数码相机的Raw文件进行数据处理。
ISO规定的数码相机特征化方法分为两种：基于光谱法和基于目标样本法。基于光谱法虽然精度较高，但需要使用专业设备如单色仪，搭建和调试系统比较复杂；而基于目标样本法由于其实验环境易于构建，得到广泛的应用。目前研究的基于目标样本的相机特征化方法有三维查找表法、多项式法等。其中，多项式特征化方法具有项数可扩展，方法简单易于使用，精度高等特点，故在相机特征化中优势明显。
基于目前的目标样本特征化方法将相机应用于场景颜色三刺激值时存在局限：相机是一种多个参数可以调节的成像式设备，其中主要影响Raw文件的可控参数有：ISO感光度、快门时间、光圈数、焦距(若使用变焦镜头)。但是，传统的特征化方法并没有将这些因素考虑在模型中，因此在训练和测试时相机的参数必须保持固定，不可变动。一旦测试场景需要相机改变参数，相机必须重新在训练环境中以改变后的参数来训练模型，这大大限制了传统特征化方法的使用，给实际使用带来了不便。
发明内容
为了克服现有相机测色方法在训练和实测时相机参数必须固定的缺陷,本发明提供一种基于齐次多项式数码相机参数可变的获取三刺激值方法。
本发明的目的是通过以下技术方案来实现的：一种基于齐次多项式数码相机参数可变的获取三刺激值方法，包括以下步骤：
(1)对数码相机进行特征化训练：获得训练样本均匀性校正后的RGB值和CIE三刺激值XYZ，以及变换系数矩阵M，记录训练时的相机拍摄参数；具体包括以下子步骤：
(1.1)拍摄训练样本：选择颜色信息丰富的标准色卡作为训练样本，在光谱稳定的光源照明条件下，设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含训练样本，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；拍摄获取色卡照片，得到色卡色块的RGB值；移除色卡，用一张均匀灰卡放置于色卡位置再次拍摄得到灰卡照片用于均匀性校正；均匀性校正后的色卡RGB值记为RGBSi,角标Si表示第i个训练样本，同时记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TS、SS、NS；
(1.2)获取训练样本的色度值，使用分光光度计测得色卡和灰卡的光谱反射比RSi和RGrey，使用分光辐射度计测得灰卡在步骤1.1中所在位置的光谱功率分布LGrey，由下式计算得到此光源照明条件下色卡的光谱功率分布LSi，进而得到其CIE三刺激值XYZ，记为XYZSi；
(1.3)由下式得到最小二乘的训练样本RGB值到XYZ值的变换系数矩阵M：
M＝HRT(RTR)-1
式中，R为可以拓展项数的RGBSi齐次多项式矩阵，必须保证齐次多项式次数为一阶，R矩阵大小为m*n，m为齐次多项式项数，n为训练样本的色块个数；H为色卡的三刺激值XYZSi构成的XYZ矩阵，大小为3*n；变换系数矩阵M大小为3*m。一阶齐次多项式矩阵R的第k行第i列的通式可以表示为下式，其中ri、gi、bi分别表示第i个训练样本的红、绿、蓝通道响应值，l、m、n为无公约数不全为0的非负整数；
(2)将数码相机应用于实际的成像式测色：由拍摄到的待测物体色的图像RGB值和拍摄时的拍摄参数，计算得到待测物体颜色的三刺激值XYZ；具体包括以下子步骤：
(2.1)在实际成像式测色应用中，为了获取待测试物体的高质量图像，重新设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含测试物体，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；对于测试图像中任一待测点RGB值，记为RGBTj，角标Tj表示第j个测试点；记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TT、ST、NT；
(2.2)将测试拍摄参数下得到的RGB响应值RGBTj转换至等效RGB值RGBEj，如下式：
(2.3)按照步骤1.3训练过程中的齐次多项式形式，扩展等效RGB值RGBEj成为形式上一致的一阶齐次多项式RE；
(2.4)利用变换系数矩阵M，得到待测点的颜色三刺激值HT，计算方法如下式：
HT＝MRE
本发明的有益效果是：本发明通过建立拍摄参数和响应值的关系，将相机拍摄参数引入特征化模型，可以方便地将不同参数的RGB响应值转换到训练参数对应的等效RGB响应值，通过拍摄一次训练样本就可以建立基于齐次多项式的颜色特征化模型，根据测试场景的需求自由调整相机参数以对场景进行色度数据采集和计算。本发明克服了传统数码相机特征化模型必须在固定相机参数下使用的局限，大大增加了数码相机作为测色设备的自由度，显著提高其适应性和实用性。针对目前普遍使用的多项式模型，一阶齐次多项式模型兼顾了等效RGB值与CIEXYZ色度值之间的全局线性特征和局部非线性特征，提高了预测的精度。相比于复杂的缩放搜索方法，齐次多项式模型的数学表达和程序实现都更简单，实际应用时运算速度更快，可以满足针对图像实时测量的要求。
附图说明
图1是基于齐次多项式数码相机参数可变的获取三刺激值方法的流程图；
图2是模型训练时实验装置图；
图3是实例中用作训练和检验的样本色卡；
图4是不同参数下的测试结果图。
具体实施方式
以尼康数码单反相机D3x配尼康AF-S尼克尔24-120mm f/4G ED VR镜头为例，阐述使用这台数码相机成像式获取物体颜色三刺激值的方法。需说明的是，本发明不限于尼康数码单反相机，只要能获取到Raw文件格式的数码相机均适用于本发明。
如图1所示，本发明一种基于齐次多项式数码相机参数可变的获取三刺激值方法，具体包括以下步骤：
(1)对数码相机进行特征化训练：获得训练样本均匀性校正后的RGB值和CIE三刺激值XYZ，以及变换系数矩阵M，记录训练时的相机拍摄参数；具体包括以下子步骤：
(1.1)拍摄训练样本。选择GretagMacbeth ColorChecker DC色卡为训练色卡，如图3所示。照明环境选择在LED灯箱中，使用D65色温50lx照度的照明条件。D3x的影像记录选项中选择NEF(RAW)格式输出，其他参数调整至构图包含训练样本，曝光正确(EV＝0)，ISO感光度为100，光圈数为5.6，快门为1s，焦距为30mm。由此，得到了色卡的RGB值。再用一张灰色均匀色卡放置于DC色卡前拍摄，用于色卡RGB值的均匀性校正。校正后的DC色卡RGB值记为RGBSi,角标Si表示第i个训练色块，除去外围一圈重复的灰阶，训练样本总共为180个色块。
(1.2)获取训练样本的色度值，使用SP64分光光度计测得色卡和灰卡的光谱反射比RSi和RGrey，使用CS-2000分光辐射度计测得灰色色卡在步骤1.1中色卡位置的光谱功率分布LGrey，如图2实验设置所示，CS-2000分光辐射度计放置在相机拍摄训练样本的位置，以保证相同的观察几何条件。由此，可计算得到此光源照明条件下色卡的CIE三刺激值XYZ，记为XYZSi,角标Si表示第i个训练色块。
(1.3)由下式得到最小二乘的RGB到XYZ的变换系数矩阵M
M＝HRT(RTR)-1
其中R选择使用7项多项式即分别为矩阵大小为7*180，变换系数矩阵为3*7大小。
(2)将数码相机应用于实际的成像式测色：由拍摄到的待测物体色的图像RGB值和拍摄时的拍摄参数，计算得到待测物体颜色的三刺激值XYZ；具体包括以下子步骤：
(2.1)获取测试图像，仍然使用DC色卡做测试，但是照明条件改变，分别在12lx，25lx，100lx，200lx，400lx，800lx共六种照明环境下使用合适的相机设置参数。在本次实验中，不改变ISO感光度、光圈、焦距，只改变快门，对应的快门分别为4s，2s，1/2s，1/4s，1/8s，1/15s，由此分别拍摄得到测试图像。
(2.2)将测试拍摄参数下得到的RGB响应值RGBTj转换至等效RGB值RGBEj，如下式：
式中NT、TT、ST表示拍摄测试样本的光圈数、曝光时间和ISO感光度，在此次实验中以12lx条件为例，分别为5.6、4s、100；NS，TS，SS表示拍摄训练样本的光圈数、曝光时间和ISO感光度，在此次实验中分别为5.6、1s、100。N、T、S这些参数均可以在Raw或者jpg的头文件中读取，可通过编写程序来自动读取，从而实现特征化的自动化过程。至此，完成了不同参数下的特征化映射。
(2.3)按照训练过程1.3步骤中的齐次多项式形式,扩展等效RGB值RGBEj成为形式上一致的一阶齐次多项式RE；
(2.4)利用变换系数矩阵M，得到待测点的颜色三刺激值HT，计算方法如下式：
HT＝MRE
实验测试结果如图4所示，不同参数设置条件下的特征化精度均在2个CIELAB色差单位以内，实现了数码相机在不同参数设置下对待测物体进行颜色三刺激值的测量。
Claims (1)
1.一种基于齐次多项式数码相机参数可变的获取三刺激值方法，其特征在于，包括以下步骤：
(1)对数码相机进行特征化训练：获得训练样本均匀性校正后的RGB值和CIE三刺激值XYZ，以及变换系数矩阵M，记录训练时的相机拍摄参数；具体包括以下子步骤：
(1.1)拍摄训练样本：选择颜色信息丰富的标准色卡作为训练样本，在光谱稳定的光源照明条件下，设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含训练样本，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；拍摄获取色卡照片，得到色卡色块的RGB值；移除色卡，用一张均匀灰卡放置于色卡位置再次拍摄得到灰卡照片用于均匀性校正；均匀性校正后的色卡RGB值记为RGBSi,角标Si表示第i个训练样本，同时记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TS、SS、NS；
(1.2)获取训练样本的色度值，使用分光光度计测得色卡和灰卡的光谱反射比RSi和RGrey，使用分光辐射度计测得灰卡在步骤(1.1)中所在位置的光谱功率分布LGrey，由下式计算得到此光源照明条件下色卡的光谱功率分布LSi，进而得到其CIE三刺激值XYZ，记为XYZSi；
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>G</mi>
<mi>r</mi>
<mi>e</mi>
<mi>y</mi>
</mrow>
</msub>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>S</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>R</mi>
<mrow>
<mi>G</mi>
<mi>r</mi>
<mi>e</mi>
<mi>y</mi>
</mrow>
</msub>
</mfrac>
</mrow>
(1.3)由下式得到最小二乘的训练样本RGB值到XYZ值的变换系数矩阵M：
M＝HRT(RTR)-1
式中，R为可以拓展项数的RGBSi齐次多项式矩阵，必须保证齐次多项式次数为一阶，R矩阵大小为m*n，m为齐次多项式项数，n为训练样本的色块个数；H为色卡的三刺激值XYZSi构成的XYZ矩阵，大小为3*n；变换系数矩阵M大小为3*m；一阶齐次多项式矩阵R的第k行第i列的通式可以表示为下式，其中ri、gi、bi分别表示第i个训练样本的红、绿、蓝通道响应值，l、m、n为无公约数不全为0的非负整数；
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>k</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mroot>
<mrow>
<msubsup>
<mi>r</mi>
<mi>i</mi>
<mi>l</mi>
</msubsup>
<msubsup>
<mi>g</mi>
<mi>i</mi>
<mi>m</mi>
</msubsup>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mi>n</mi>
</msubsup>
</mrow>
<mrow>
<mi>l</mi>
<mo>+</mo>
<mi>m</mi>
<mo>+</mo>
<mi>n</mi>
</mrow>
</mroot>
</mrow>
(2)将数码相机应用于实际的成像式测色：由拍摄到的待测物体色的图像RGB值和拍摄时的拍摄参数，计算得到待测物体颜色的三刺激值XYZ；具体包括以下子步骤：
(2.1)在实际成像式测色应用中，为了获取待测试物体的高质量图像，重新设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含测试物体，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；对于测试图像中任一待测点RGB值，记为RGBTj，角标Tj表示第j个测试点；记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TT、ST、NT；
(2.2)将测试拍摄参数下得到的RGB响应值RGBTj转换至等效RGB值RGBEj，如下式：
<mrow>
<msub>
<mi>RGB</mi>
<mrow>
<mi>E</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>S</mi>
</msub>
<msub>
<mi>S</mi>
<mi>S</mi>
</msub>
<msubsup>
<mi>N</mi>
<mi>T</mi>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mi>T</mi>
</msub>
<msub>
<mi>S</mi>
<mi>T</mi>
</msub>
<msubsup>
<mi>N</mi>
<mi>S</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<msub>
<mi>RGB</mi>
<mrow>
<mi>T</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
(2.3)按照步骤(1.3)训练过程中的齐次多项式形式，扩展等效RGB值RGBEj成为形式上与之一致的一阶齐次多项式RE；
(2.4)利用变换系数矩阵M，得到待测点的颜色三刺激值HT，计算方法如下式：
HT＝MRE。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610236691.4A CN105744267B (zh) | 2016-04-15 | 2016-04-15 | 基于齐次多项式数码相机参数可变的获取三刺激值方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610236691.4A CN105744267B (zh) | 2016-04-15 | 2016-04-15 | 基于齐次多项式数码相机参数可变的获取三刺激值方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105744267A CN105744267A (zh) | 2016-07-06 |
| CN105744267B true CN105744267B (zh) | 2017-08-25 |
Family
ID=56254499
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610236691.4A Expired - Fee Related CN105744267B (zh) | 2016-04-15 | 2016-04-15 | 基于齐次多项式数码相机参数可变的获取三刺激值方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105744267B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110738957B (zh) * | 2019-09-17 | 2023-04-25 | 苏州佳世达电通有限公司 | 显示系统及颜色特性测量方法 |
| CN115597726B (zh) * | 2022-09-30 | 2023-12-15 | 深圳市汇顶科技股份有限公司 | 色温检测的方法、装置和电子设备 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104168475A (zh) * | 2014-08-15 | 2014-11-26 | 浙江大学 | 一种数码相机参数可变的成像式获取颜色三刺激值的方法 |
-
2016
- 2016-04-15 CN CN201610236691.4A patent/CN105744267B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104168475A (zh) * | 2014-08-15 | 2014-11-26 | 浙江大学 | 一种数码相机参数可变的成像式获取颜色三刺激值的方法 |
Non-Patent Citations (2)
| Title |
|---|
| A study of digital camera colorimetric characterisation based on polynomial modeling;G Hong等;《Color Research &Application》;20011231;全文 * |
| Practical Spectral Characterization of Trichromatic Cameras;Rump M等;《ACM TRANSACTIONS ON GRAPHICS》;20111231;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105744267A (zh) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104168475B (zh) | 一种数码相机参数可变的成像式获取颜色三刺激值的方法 | |
| US9514535B1 (en) | Color calibration method of camera module | |
| JP6974397B2 (ja) | 測定装置 | |
| Berns | The science of digitizing paintings for color-accurate image archives: a review | |
| JP6455764B2 (ja) | 色補正パラメータ算出方法、色補正パラメータ算出装置および画像出力システム | |
| JP6257551B2 (ja) | 色忠実環境補正装置および色忠実環境補正方法 | |
| EP3888345B1 (en) | Method for generating image data for machine learning based imaging algorithms | |
| CN110017904B (zh) | 一种基于ccd相机的多光谱辐射测温方法 | |
| CN101933321A (zh) | 用于估测场景光源的图像传感器装置及方法 | |
| CN101939997A (zh) | 采用与光源相关的色彩校正矩阵进行色彩校正的图像传感器装置及方法 | |
| US7616314B2 (en) | Methods and apparatuses for determining a color calibration for different spectral light inputs in an imaging apparatus measurement | |
| CN1981512A (zh) | 用于对照相机和/或显示设备进行色彩校准和用于校正数字图像色彩缺陷的方法和设备 | |
| US20140267782A1 (en) | Apparatus And Method For Automated Self-Training Of White Balance By Electronic Cameras | |
| CN105812661B (zh) | 一种基于标准灯箱和灰卡的数码相机均匀性校正方法 | |
| CN115100082A (zh) | 一种基于高光谱相机的高精度颜色显示系统 | |
| Varghese et al. | Colorimetric calibration of high dynamic range images with a ColorChecker chart | |
| CN105744267B (zh) | 基于齐次多项式数码相机参数可变的获取三刺激值方法 | |
| Pointer et al. | Practical camera characterization for colour measurement | |
| Ejaz et al. | Development of a camera system for the acquisition of high-fidelity colors | |
| Kretkowski et al. | Development of an XYZ digital camera with embedded color calibration system for accurate color acquisition | |
| Guarini et al. | A practical procedure for obtaining calibrated material colors for CAD Systems | |
| CN106896068A (zh) | 解决网购颜色失真的商品颜色获取与再现方法及其系统 | |
| CN111896109A (zh) | 一种基于数码相机原始响应值的光谱重建方法 | |
| Guarnera et al. | Absolute colorimetric characterization of a dslr camera | |
| Kim et al. | Developing a multispectral HDR imaging module for a BRDF measurement system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170825 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |