一种利用多光谱相机虚拟响应值获取光谱反射比的方法技术

本发明专利技术公开了一种利用多光谱相机虚拟响应值获取光谱反射比的方法，首先将多光谱相机的所有通道划分为奇数通道和偶数通道，利用训练样本建立多光谱相机偶数和奇数通道响应值之间的函数关系，用奇数通道采集测试目标的相机响应值，根据建立的偶数和奇数通道响应值之间的函数关系，计算出偶数通道的虚拟相机响应值。利用测试目标所有奇数和偶数通道的响应值计算出其光谱反射比。本发明专利技术只需要利用多光谱相机一半的通道去拍摄就能获得全部通道的多光谱相机响应值，显著减少了通道数量的使用而又能充分利用所有通道的信息来获取物体的光谱反射比，从而在确保高精度获取物体光谱反射比的同时有效提高获取效率。

【技术实现步骤摘要】

本专利技术涉及成像式获取物体光谱及色度信息的方法，尤其是利用具有多个可见光波段通道的多光谱相机获取物体的反射光谱信息以进一步获取物体色度信息的方法。
技术介绍
由于同色异谱现象的存在，物体的色度信息难以可靠表征物体的本源信息，而物体的光谱反射比是设备无关且与捕获光源无关的物理量，可以表征物体的本源信息。因此，准确获取物体的光谱反射比可以实现高精度的颜色测量，且可以忠实复现物体在任意光源下的颜色外貌。成像式获取物体光谱反射比可以克服分光光度计逐点式测量和接触式测量的缺陷，能够快速获取成像范围内所有物体的光谱反射比。普通商用相机可以用于成像式获取物体光谱反射比，但由于只有R、G、B三个通道而无法高精度获取光谱反射比。多光谱相机用于成像式获取物体光谱反射比时，由于增加了更多的通道进行成像，故可大幅度提高光谱反射比获取的精度。多光谱相机通道数的增加也延长了成像的时间，增加了操作的复杂度，例如需要对每个通道调整曝光时间，对每个通道聚焦，多通道图像之间配准等，这一系列操作随着通道数的增多而变得更加复杂和耗时。同时，成像时间的延长所带来的光源和拍摄场景的时间不稳定性也会进一步影响物体光谱反射比的获取精度。
技术实现思路
为了降低多光谱相机的成像时间和操作复杂度，同时确保光谱反射比的获取精度，本专利技术提出了一种利用多光谱相机虚拟响应值获取光谱反射比的方法。本专利技术的目的是通过以下技术方案来实现的：一种利用多光谱相机虚拟响应值获取光谱反射比的方法，包括以下步骤：(1)利用训练样本建立偶数通道响应值和奇数通道相机响应值之间的函数关系。具体包括以下子步骤：(1.1)用多光谱相机拍摄训练样本。采用包含各种颜色块的标准色卡作为训练样本，选取光谱功率分布稳定的光源作为照明光源，将标准色卡以45/0照明与测量几何条件放置于照明光源下，用多光谱相机的所有通道拍摄标准色卡，提取标准色卡的通道响应值P为n行q列的矩阵，其中n为标准色卡上的色块数，q为多光谱相机通道的数量，每一行表示一个色块对应的q个通道的响应值，每一列表示一个通道输出的n个色块的响应值。(1.2)建立每一个偶数通道响应值与前后相邻奇数通道响应值之间的函数关系假设p1，p2和p3分别为标准色卡响应值矩阵P的第1，2和3列，各对应第1，2和3个通道的标准色卡上n个色块的响应值。通过二阶多项式拟合第2个通道响应值p2与第1个和第3个通道响应值p1和p3之间的关系，即p2＝a0+a1*p1+a2*p3+a3*p12+a4*p32+a5*p1*p3利用p2和p1，p3，通过最小二乘法求解上述多项式的系数A＝[a0,a1,a2,a3,a4,a5]，即A＝(CTC)-1CTp2其中，为含有n个元素都为1的列向量，为矩阵或向量的Hadamard积。用上述同样的方法，类似地可以求出第4,6,8等偶数通道响应值与其前后相邻通道响应值之间的二阶多项式函数。(2)计算光谱反射比与多通道响应值之间的转换矩阵。用分光光度计测量标准色卡的所有色块的光谱反射比R为n行m列的矩阵，其中n为标准色卡上的色块数，m为光谱反射比的维数。所有奇数通道的响应值为所有偶数通道的响应值为根据每一个偶数通道响应值和前后相邻奇数通道响应值的二阶多项式函数，计算出标准色卡中所有偶数通道的虚拟响应值则标准色卡的所有通道的响应值然后，由最小二乘法计算从多通道响应值到光谱反射比的转换矩阵M，即 M = R T P ^ ( P ^ T P ^ ) - 1 ]]>(3)用奇数通道拍摄目标样本，并计算目标样本的光谱反射比。在相同光源下，用多光谱相机的所有奇数通道拍摄目标样本，得到目标样本的奇数通道响应值Pto,根据步骤(1)中偶数通道与奇数通道响应值之间的函数关系，计算出目标样本所有偶数通道的响应值则目标样本的所有通道的响应值根据步骤(2)中多通道响应值到光谱反射比的转换矩阵M,计算目标样本的光谱反射比 R ^ t = ( M P ^ t T ) T . ]]>本专利技术的有益效果是：本专利技术预先建立偶数通道响应值和奇数通道响应值之间的函数关系，然后仅利用奇数通道对目标成像，根据偶数通道响应值和奇数通道响应值之间的函数关系，利用奇数通道响应值计算出偶数通道响应值，从而利用所有奇数通道和偶数通道的响应值来获取物体的光谱反射比，确保了高精度光谱反射比的获取，同时显著减少了实际拍摄时所用通道的数量，有效降低了多通道图像的采集时间和操作复杂度。附图说明图1是利用多光谱相机虚拟响应值获取光谱反射比的流程图；图2是实例中多光谱相机结构示意图；图3是实例中采用虚拟通道响应值前后分别获得的光谱反射比。具体实施方式以一台滤色片轮式的多光谱相机为例，阐述利用虚拟通道响应值获取光谱反射比的方法。多光谱相机的结构示意图如图2，由镜头、滤色片轮、滤色片轮控制单元和单色CCD传感器构成。该多光谱相机利用16个干涉滤色片形成16个通道，滤色片轮一次能安装8个滤色片，因此16个滤色片需要两次安装。16个滤色片的峰值透过率波长为从400nm到700nm,FWHM(峰值半高宽,Full Width of Half Maximum)为20nm,因此8个奇数通道对应的峰值透过率波长分别为400nm，440nm，480nm，520nm，560nm，600nm，640nm，680nm；8个偶数通道对应的峰值透过率波长分别为420nm，460nm，500nm，540nm，580nm，620nm，660nm，700nm。需要说明的是，本专利技术不仅限于滤色片轮式的多光谱相机，对于其他形式的多光谱相机，如基于LCTF(液晶可调滤色片，Liquid Crystal Tuneable Filter)的多光谱相机也适用；同时，本专利技术也不仅限于通道数量为16的多光谱相机，少于或多于该滤色片数量的多光谱相机也适用。此外，也可以预先建立奇数通道与偶数通道之间的二阶多项式函数关系，然后仅用偶数通道对目标样本成像，然后利用奇数通道和偶数通道之间的函数关系计算出目标样本的奇数通道响应值，最后则利用全部通道响应值来获取物体光谱反射比。如图1所示，本专利技术利用多光谱相机虚拟响应值获取光谱反射比的方法，具体包括以下步骤：(1)利用训练样本建立偶数通道响应值和奇数通道相机响应值之间的函数关系。具体包括以下子步骤：(1.1)用多本文档来自技高网...

【技术保护点】
一种利用多光谱相机虚拟响应值获取光谱反射比的方法，其特征在于，包括以下步骤：(1)利用训练样本建立偶数通道响应值和奇数通道相机响应值之间的函数关系。具体包括以下子步骤：(1.1)用多光谱相机拍摄训练样本。采用包含各种颜色块的标准色卡作为训练样本，选取光谱功率分布稳定的光源作为照明光源，将标准色卡以45/0照明与测量几何条件放置于照明光源下，用多光谱相机的所有通道拍摄标准色卡，提取标准色卡的通道响应值P为n行q列的矩阵，其中n为标准色卡上的色块数，q为多光谱相机通道的数量，每一行表示一个色块对应的q个通道的响应值，每一列表示一个通道输出的n个色块的响应值。(1.2)建立每一个偶数通道响应值与前后相邻奇数通道响应值之间的函数关系。假设p1，p2和p3分别为标准色卡响应值矩阵P的第1，2和3列，各对应第1，2和3个通道的标准色卡上n个色块的响应值。通过二阶多项式拟合第2个通道响应值p2与第1个和第3个通道响应值p1和p3之间的关系，即p2＝a0+a1*p1+a2*p3+a3*p12+a4*p32+a5*p1*p3利用p2和p1，p3，通过最小二乘法求解上述多项式的系数A＝[a0,a1,a2,a3,a4,a5]，即A＝(CTC)‑1CTp2其中，为含有n个元素都为1的列向量，为矩阵或向量的Hadamard积。用上述同样的方法，类似地可以求出第4,6,8等偶数通道响应值与其前后相邻通道响应值之间的二阶多项式函数。(2)计算光谱反射比与多通道响应值之间的转换矩阵。用分光光度计测量标准色卡中所有色块的光谱反射比R为n行m列的矩阵，其中n为标准色卡上的色块数，m为光谱反射比的维数。所有奇数通道的响应值为所有偶数通道的响应值为根据每一个偶数通道响应值和前后相邻奇数通道响应值的二阶多项式函数，计算出标准色卡的所有偶数通道的虚拟响应值则标准色卡的所有通道的响应值然后，由最小二乘法计算从多通道响应值到光谱反射比的转换矩阵M，即M=RTP^(P^TP^)-1]]>(3)用奇数通道拍摄目标样本，并计算目标样本的光谱反射比。在相同光源下，用多光谱相机的所有奇数通道拍摄目标样本，得到目标样本的奇数通道响应值Pto,根据步骤(1)中偶数通道与奇数通道响应值之间的函数关系，计算出目标样本所有偶数通道的响应值则目标样本的所有通道的响应值根据步骤(2)中多通道响应值到光谱反射比的转换矩阵M,计算目标样本的光谱反射比R^t=(MP^tT)T.]]>...

【技术特征摘要】
1.一种利用多光谱相机虚拟响应值获取光谱反射比的方法，其特征在于，包括以下步骤：(1)利用训练样本建立偶数通道响应值和奇数通道相机响应值之间的函数关系。具体包括以下子步骤：(1.1)用多光谱相机拍摄训练样本。采用包含各种颜色块的标准色卡作为训练样本，选取光谱功率分布稳定的光源作为照明光源，将标准色卡以45/0照明与测量几何条件放置于照明光源下，用多光谱相机的所有通道拍摄标准色卡，提取标准色卡的通道响应值P为n行q列的矩阵，其中n为标准色卡上的色块数，q为多光谱相机通道的数量，每一行表示一个色块对应的q个通道的响应值，每一列表示一个通道输出的n个色块的响应值。(1.2)建立每一个偶数通道响应值与前后相邻奇数通道响应值之间的函数关系。假设p1，p2和p3分别为标准色卡响应值矩阵P的第1，2和3列，各对应第1，2和3个通道的标准色卡上n个色块的响应值。通过二阶多项式拟合第2个通道响应值p2与第1个和第3个通道响应值p1和p3之间的关系，即p2＝a0+a1*p1+a2*p3+a3*p12+a4*p32+a5*p1*p3利用p2和p1，p3，通过最小二乘法求解上述多项式的系数A＝[a0,a1,a2,a3,a4,a5]，即A＝(CTC)-1CTp2其中，为含有n个元素都为1的列向量，为矩阵或向量的Hadamard积。用上述同样的方法，类似地可以求出第4,6,8等偶数通道响应值与其前后相邻通道响应值之间的二阶多项式函数。(2)计算光谱反射比与多通道响应值之间的转换矩阵。用分光光度计测量标准色卡中所有色块的光谱反射比R为n行m列的矩阵，其中n为标准色卡上的色块数，m为光谱反射比的维数。所有奇数通道的响应值为所有偶数通道的响应值为根据每一个偶数通道...

【专利技术属性】
技术研发人员：徐海松，徐鹏，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江;33