发明名称 一种数码相机参数可变的成像式获取颜色三刺激值的方法
摘要 本发明公开了一种数码相机参数可变的成像式获取颜色三刺激值的方法。该方法将数码相机拍摄设置参数引入作为输入变量,利用缩放和还原,将实测样本和训练样本不同参数设置建立关联,在训练样本的尺度上使用多项式模型,将RGB空间转化至XYZ空间,最后还原实测样本的真实XYZ三刺激值。本发明克服了传统方法中要求相机必须固定拍摄参数进行训练和实测的缺陷,使得数码相机可以在固定一种拍摄参数下训练,自由变化拍摄参数实际拍摄测量,为数码相机作为一种成像式测色仪器提供了便利。
申请公布号 CN104168475A 申请公布日期 2014.11.26
申请号 CN201410403603.6 申请日期 2014.08.15
申请人 浙江大学 发明人 徐海松;方竞宇;汪哲弘;吴晓玟
分类号 H04N17/00(2006.01)I;H04N17/02(2006.01)I 主分类号 H04N17/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种数码相机参数可变的成像式获取颜色三刺激值的方法,其特征在于,包括以下步骤:(1)对数码相机进行特征化训练:获得训练样本均匀性校正后的RGB值和XYZ值,以及变换系数矩阵M,记录训练时的相机拍摄参数;具体包括以下子步骤:(1.1)拍摄训练样本:选择颜色信息丰富的标准色卡作为训练样本,在光谱稳定的光源照明条件下,设置相机拍摄参数:ISO感光度、光圈、快门、焦距,使得拍摄图像包含训练样本,且曝光正确;同时设置相机储存参数,使得输出图像包含Raw文件格式;拍摄获取色卡照片,得到色卡色块的RGB值;移除色卡,用一张均匀灰卡放置于色卡位置再次拍摄得到灰卡照片用于均匀性校正;均匀性校正后的色卡RGB值记为R<sub>Si</sub>,G<sub>Si</sub>,B<sub>Si</sub>,角标Si表示第i个训练样本;通过下式将RGB<sub>Si</sub>从RGB空间转换到rg空间,记为r<sub>Si</sub>,g<sub>Si</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mi>Si</mi></msub><mo>=</mo><mfrac><msub><mi>R</mi><mi>Si</mi></msub><mrow><msub><mi>R</mi><mi>Si</mi></msub><mo>+</mo><msub><mi>G</mi><mi>Si</mi></msub><mo>+</mo><msub><mi>B</mi><mi>Si</mi></msub></mrow></mfrac><mo>,</mo><msub><mi>g</mi><mi>Si</mi></msub><mo>=</mo><mfrac><msub><mi>G</mi><mi>Si</mi></msub><mrow><msub><mi>R</mi><mi>Si</mi></msub><mo>+</mo><msub><mi>G</mi><mi>Si</mi></msub><mo>+</mo><msub><mi>B</mi><mi>Si</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000554550920000011.GIF" wi="867" he="140" /></maths>(1.2)获取训练样本的色度值,使用分光光度计测得色卡和灰卡的光谱反射比R<sub>Si</sub>和R<sub>Grey</sub>,使用分光辐射度计测得灰卡在步骤1.1中所在位置的光谱功率分布L<sub>Grey</sub>,计算得到此光源照明条件下色卡的CIE三刺激值XYZ,记为X<sub>Si</sub>,Y<sub>Si</sub>,Z<sub>Si</sub>,角标Si表示第i个训练样本;(1.3)由下式得到最小二乘的训练样本RGB值到XYZ值的变换系数矩阵M:M=(R<sup>T</sup>R)<sup>‑1</sup>R<sup>T</sup>H式中,R为可以拓展项数的R<sub>Si</sub>,G<sub>Si</sub>,B<sub>Si</sub>多项式矩阵,矩阵大小为m*n,m为多项式项数,n为训练样本的色块个数;H为色卡的三刺激值构成的XYZ矩阵,大小为3*n;变换系数矩阵M大小为3*m;(2)将数码相机应用于实际的成像式测色:由拍摄到的待测色物体的照片RGB值和拍摄时的拍摄参数,计算得到待测物体的颜色三刺激值XYZ;具体包括以下子步骤:(2.1)在实际成像式测色应用中,获取待测试物体的图像,重新设置相机拍摄参数:ISO感光度、光圈、快门、焦距,使得拍摄图像包含测试物体,且曝光正确;同时设置相机储存参数,使得输出图像包含Raw文件格式;对于测试图像中任一待测点RGB值,记为RGB<sub>Tj</sub>,角标Tj表示第j个测试点,将RGB<sub>Tj</sub>从RGB空间转换到rg空间,记为rg<sub>Tj</sub>;(2.2)缩放测试点的RGB值:在rg空间上,寻找与rg<sub>Tj</sub>距离最近的训练样本,记为rg<sub>Si’</sub>,缩放比例因子sl计算如下式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>sl</mi><mo>=</mo><mfrac><mrow><msub><mi>R</mi><msup><mi>Si</mi><mo>&prime;</mo></msup></msub><mo>+</mo><msub><mi>G</mi><msup><mi>Si</mi><mo>&prime;</mo></msup></msub><mo>+</mo><msub><mi>B</mi><msup><mi>Si</mi><mo>&prime;</mo></msup></msub></mrow><mrow><msub><mi>R</mi><mi>Tj</mi></msub><mo>+</mo><msub><mi>G</mi><mi>Tj</mi></msub><mo>+</mo><msub><mi>B</mi><mi>Tj</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000554550920000012.GIF" wi="411" he="145" /></maths>然后对测试点进行缩放:<img file="FDA0000554550920000013.GIF" wi="307" he="84" />其中X=R,G,B,<img file="FDA0000554550920000014.GIF" wi="86" he="84" />表示测试点缩放后的RGB值;(2.3)计算由测试点缩放的XYZ值构成的矩阵<img file="FDA0000554550920000021.GIF" wi="111" he="78" />由多项式变换系数矩阵计算<img file="FDA0000554550920000022.GIF" wi="266" he="78" /><img file="FDA0000554550920000023.GIF" wi="77" he="78" />为步骤2.2测试点缩放后的RGB值构成的多项式矩阵;(2.4)将步骤2.3测试点缩放的XYZ值还原,得到待测点的颜色三刺激值H<sub>T</sub>,计算如下式:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>H</mi><mi>T</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>sl</mi></mfrac><mo>&CenterDot;</mo><mfrac><mrow><msub><mi>N</mi><mi>T</mi></msub><msub><mi>T</mi><mi>S</mi></msub><msub><mi>S</mi><mi>S</mi></msub></mrow><mrow><msub><mi>N</mi><mi>S</mi></msub><msub><mi>T</mi><mi>T</mi></msub><msub><mi>S</mi><mi>T</mi></msub></mrow></mfrac><msubsup><mi>H</mi><mi>T</mi><mi>sl</mi></msubsup></mrow>]]></math><img file="FDA0000554550920000024.GIF" wi="457" he="140" /></maths>式中N<sub>T</sub>,T<sub>T</sub>,S<sub>T</sub>表示拍摄测试图像的光圈数、曝光时间和ISO感光度,N<sub>S</sub>,T<sub>S</sub>,S<sub>S</sub>表示拍摄训练样本的光圈数、曝光时间和ISO感光度。
地址 310058 浙江省杭州市西湖区余杭塘路866号