基于齐次多项式数码相机参数可变的获取三刺激值方法技术

本发明专利技术公开了一种基于齐次多项式数码相机参数可变的获取三刺激值方法。该方法将数码相机拍摄设置参数引入作为输入变量，建立实测样本和训练样本不同参数设置下相机响应值关联，基于一阶齐次多项式模型，将等效RGB值转换为实测样本的XYZ三刺激值。本发明专利技术克服了传统方法中要求相机必须固定拍摄参数进行训练和实测的缺陷，使得数码相机可以在固定一种拍摄参数下训练，而在实际拍摄测量时则可以自由变化拍摄参数，从而为数码相机作为一种成像式测色仪器提供了便利。

【技术实现步骤摘要】

本专利技术涉及颜色管理系统中的基于数码相机成像式获取颜色三刺激值的方法，尤其是在数码相机的设置参数变化的情况下训练和实际测量的方法。
技术介绍
目前，相对传统的测色设备，使用商用的数码相机作为一种成像式的测色仪器具有成本低，操作简便等优点，因而在定量捕获场景色度信息的领域得到了越来越多的应用。然而，相机的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；LSi=LGreyRSiRGrey]]>(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的非负整数；ρki=rilgimbinl+m+n]]>(2)将数码相机应用于实际的成像式测色：由拍摄到的待测物体色的图像RGB值和拍摄时的拍摄参数，计算得到待测物体颜色的三刺激值XYZ；具体包括以下子步骤：(2.1)在实际成像式测色应用中，为了获取待测试物体的高质量图像，重新设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含测试物体，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；对于测试图像中任一待测点RGB值，记为RGBTj，角标Tj表示第j个测试点；记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TT、ST、NT；(2.2)将测试拍摄参数下得到的RGB响应值RGBTj转换至等效RGB值RGBEj，如下式：RGBEj=TSSSNT2TTSTNS2RGBTj]]>(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-120mmf/4GEDVR镜头为例，阐述使用这台数码相机成像式获取物体颜色三刺激值的方法。需说明的是，本专利技术不限于尼康数码单反相机，只要能获取到Raw文件格式的数码相机均适用于本专利技术。如图1所示，本专利技术一种基于齐次多项式数码相机参数可变的获取三刺激值方法，具体包括以下步骤：(1)对数码相机进行特征化训练：获得训练样本均匀性校正后的RGB值和CIE三刺激值XYZ，以及变换系数矩阵M，记录训练时的相机拍摄参数；具体包括以下子步骤：(1.1)拍摄训练样本。选择GretagMacbethColorChec本文档来自技高网...

【技术保护点】
一种基于齐次多项式数码相机参数可变的获取三刺激值方法，其特征在于，包括以下步骤：(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；LSi=LGreyRSiRGrey]]>(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的非负整数；ρki=rilgimbinl+m+n]]>(2)将数码相机应用于实际的成像式测色：由拍摄到的待测物体色的图像RGB值和拍摄时的拍摄参数，计算得到待测物体颜色的三刺激值XYZ；具体包括以下子步骤：(2.1)在实际成像式测色应用中，为了获取待测试物体的高质量图像，重新设置相机拍摄参数：ISO感光度、光圈、快门、焦距，使得拍摄图像包含测试物体，且曝光正确；同时设置相机储存参数，使得输出图像包含Raw文件格式；对于测试图像中任一待测点RGB值，记为RGBTj，角标Tj表示第j个测试点；记录相机的拍摄参数：快门时间、ISO感光度、光圈数，分别记为TT、ST、NT；(2.2)将测试拍摄参数下得到的RGB响应值RGBTj转换至等效RGB值RGBEj，如下式：RGBEj=TSSSNT2TTSTNS2RGBTj]]>(2.3)按照步骤1.3训练过程中的齐次多项式形式，扩展等效RGB值RGBEj成为形式上与之一致的一阶齐次多项式RE；(2.4)利用变换系数矩阵M，得到待测点的颜色三刺激值HT，计算方法如下式：HT＝MRE。...

【技术特征摘要】
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；
LSi=LGreyRSiRGrey]]>(1.3)由下式得到最小二乘的训练样本RGB值到XYZ值的变换系数矩阵M：
M＝HRT(RTR)-1式中，R为可以拓展项数的RGBSi齐次多项式矩阵，必须保证齐次多项式次数为一阶，R矩阵大
...

【专利技术属性】
技术研发人员：徐海松，吴晓玟，汪哲弘，方竞宇，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33