本发明涉及文物保护与鉴定技术领域，尤其涉及一种彩绘类文物混合颜料的分析方法。
背景技术：
目前，在文物保护与鉴定技术领域，混合颜料分析是挖掘古书画内在价值，为其正确的评估分析提供科学依据基础的技术。然而纸寿千年，绢寿八百，即使没有环境、人为等因素的影响，古书画出现发霉、虫蛀、污迹、破损等各种情况也非常普遍，为了更好地留存和展示古书画作品的价值，往往需要进行对其进行颜色修复和复制，这就需要确定颜料的种类。
现有技术中有很多科学技术应用于彩绘类文物的颜料类别鉴定，如x射线衍射分析、荧光x射线分析、拉曼光谱分析、电子显微分析、近红外光谱分析等技术，通过测定颜料的元素、结构以及外在形貌等特征确定物质的类别，但现有技术的方法在大多数情况下需要制样才能进行精确测量，而制样本身就是对文物的一个二次损害的过程，由于文物不可再生的特殊性，对于文物分析所采用的技术应当是无损的，而现有技术中也缺乏对混合颜料分析的解决方案。
技术实现要素：
本发明的目的是提供一种彩绘类文物混合颜料的分析方法，该方法能够永久留存彩绘类文物的颜料数字化信息，有助于对文物的数字化虚拟修复，并有效缩短文物实体的修复时间，对彩绘类文物的数字化保护与修复具有较大的经济效益和社会效益。
本发明的目的是通过以下技术方案实现的：
一种彩绘类文物混合颜料的分析方法，所述方法包括：
步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；
步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；
步骤3、将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；
步骤4、再利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得待处理彩绘类文物表面色彩组成颜料的种类和比例。
由上述本发明提供的技术方案可以看出，上述方法能够永久留存彩绘类文物的颜料数字化信息，有助于对文物的数字化虚拟修复，并有效缩短文物实体的修复时间，对彩绘类文物的数字化保护与修复具有较大的经济效益和社会效益。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。
图1为本发明实施例提供的彩绘类文物混合颜料的分析方法流程示意图。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
下面将结合附图对本发明实施例作进一步地详细描述，如图1所示为本发明实施例提供的彩绘类文物混合颜料的分析方法流程示意图，所述方法包括：
步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；
在该步骤中，宽波段连续反射光谱的覆盖范围为350nm-2500nm，采集环境为暗室，利用仪器自带的光源，测量多次取平均值，对同一种颜色采集至少二处以上浓淡不一区域的反射光谱；该宽波段连续反射光谱是由合成颜色的多种纯净颜料的光谱反射混合而成。
步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；
这里，非负矩阵分解算法要求矩阵中所有元素均为非负，在不考虑误差的情况下，将矩阵表示为两个非负矩阵完全相乘的结果，正是由于分解的结果是非负的特点，因而在实际应用中具有明确的物理意义，由于高光谱数据是对地物反射率数据的一种记载，满足非负性的要求，从而可以将非负矩阵分解算法用来高光谱数据的端元提取中来。
具体实现中，该非负矩阵分解算法nmf采用如下式子来表示：
xn×m＝an×rsr×m+en×m
其中，
在不考虑误差的影响，对于由一个已知量x求解出两个未知量a和s的问题，往往转化为下式来求解，即：
min||x-as||；
进一步，通过定义目标函数和选择优化准则来求解两个端元矩阵a和s，所定义的目标函数包括最小化欧式距离和k-l散度，其中：
所述最小化欧式距离表示为：
所述k-l散度表示为：
其中i(i＝1，2，…n)为混合光谱矩阵x的行号(即光谱通道)，其最大值就是n，n是所测量的一条光谱曲线的总的点数，与仪器性能有关，是仪器的总的波段数；j(j＝1，2，…m)是混合光谱矩阵x的列号，j的最大值是m，m代表选择几条混合反射光谱参与运算。
结合所述目标函数，不断交替迭代计算端元矩阵a和s的值，直到达到最优解，其迭代公式表示为：
从非负矩阵分解的公式来看，该算法在于寻找到一组基向量，使得目标函数最小化，对于混合颜料光谱分离来说，各个纯净颜料组分光谱就是这样的一组基向量，因而求解出这组基向量，就相当于得到了构成混合光谱的端元光谱。
步骤3、将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；
步骤4、再利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得该待处理彩绘类文物表面色彩组成颜料的种类和比例。
这里，基于众数的比值导数光谱法是建立在原有比值导数法的基础上，又考虑到存在非线性波段的因素，通过改进比值导数光谱法中特征波段的选取方法，引进众数这一统计概念，用来表征其总体特征，进而确定丰度。
具体实现中，上述利用基于众数的比值导数法来获取各个种类颜料的丰度的过程为：
首先根据色差值与颜料丰度之间的关系确定统计区间；
然后将基于比值导数法求解的各个种类颜料的丰度值按照非负且大于1的约束条件剔除异常值之后，计算落在统计区间内的个数，并将其所在区间认定为初步的丰度值范围；
再以步长0.01为比例间隔，构建一系列模拟混合光谱，计算模拟光谱与实测光谱的相关系数，相关系数最高的一组端元比例即为各个种类颜料最终的丰度值。
具体实现中，比值导数解混算法是建立在混合光谱满足线性混合模型的基础上，并满足端元丰度非负并且端元丰度之和为1的约束条件。当混合物仅包含两种矿物组分时，且不考虑误差限项时，该比值导数法的模型表示为：
r(λi)＝f1×r1(λi)+f2×r2(λi)
其中，i为光谱通道，r(λi)为混合光谱在λi波长位置的反射率，f1、f2为端元的丰度，r1(λi)、r2(λi)分别为端元f1、f2所对应的在λi波长位置的反射率；
将端元f2视为干扰组分，上述表达式两侧同时除以f2端元光谱的反射率，得到
再对r(λi)一阶求导，得到
通过该式可以看出，此时导数光谱已经和f2端元的含量无关，然后再两边同时除以
再将f1端元经过二阶求导之后的丰度结果作为各个种类颜料的初始丰度值，具体表达式为：
可以看出，比值导数光谱法求解过程，简单清晰，不需要复杂的迭代运算，能够快速、准确的进行端元光谱丰度的反演。
值得注意的是，本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
综上所述，本发明实施例所提供的方法能够永久留存文物的颜料数字化信息，有助于对文物的数字化虚拟修复，提高文物的艺术表达效果，有利于文物的线上线下数字化展示；另外，颜料的丰度反演对文物的实体修复，特别是对其全色过程能够提供科学的依据，能够有效缩短文物实体的修复时间，对我国具有大量馆藏彩绘类文物数字化保护与修复具有较大的经济效益和社会效益。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。