一种彩绘类文物混合颜料的分析方法技术

本发明专利技术公开了一种彩绘类文物混合颜料的分析方法，首先利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；然后将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；再将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得待处理彩绘类文物表面色彩组成颜料的种类和比例。上述方法能够永久留存彩绘类文物的颜料数字化信息，有助于对文物的数字化虚拟修复和实体修复。

An analytical method for mixed pigments of coloured cultural relics

The present invention discloses an analytical method for mixed pigments of painted cultural relics. Firstly, the wide-band continuous reflection spectrum of the region of interest on the surface of the painted cultural relics to be processed is collected by hyperspectral imaging and geospectroscopy, and then the wide-band continuous reflection spectrum is regarded as a multi-source mixed information, which is realized by non-negative matrix decomposition algorithm. Now the pigments in the same color region are decomposed and the reflective spectrum curves of different pigments are obtained; then the reflective spectrum curves of the different pigments are matched with the spectrum of the pigments spectrum library to get the kinds of mixed pigments; the abundance of each kind of pigments is obtained by the ratio derivative method based on the mode number, and then the pigments to be processed are obtained. The surface color of painted cultural relic forms the proportion and proportion of pigments. The above methods can permanently preserve the digital information of painted cultural relics, which is conducive to the digital virtual restoration and physical restoration of cultural relics.

【技术实现步骤摘要】
一种彩绘类文物混合颜料的分析方法
本专利技术涉及文物保护与鉴定
，尤其涉及一种彩绘类文物混合颜料的分析方法。
技术介绍
目前，在文物保护与鉴定
，混合颜料分析是挖掘古书画内在价值，为其正确的评估分析提供科学依据基础的技术。然而纸寿千年，绢寿八百，即使没有环境、人为等因素的影响，古书画出现发霉、虫蛀、污迹、破损等各种情况也非常普遍，为了更好地留存和展示古书画作品的价值，往往需要进行对其进行颜色修复和复制，这就需要确定颜料的种类。现有技术中有很多科学技术应用于彩绘类文物的颜料类别鉴定，如X射线衍射分析、荧光X射线分析、拉曼光谱分析、电子显微分析、近红外光谱分析等技术，通过测定颜料的元素、结构以及外在形貌等特征确定物质的类别，但现有技术的方法在大多数情况下需要制样才能进行精确测量，而制样本身就是对文物的一个二次损害的过程，由于文物不可再生的特殊性，对于文物分析所采用的技术应当是无损的，而现有技术中也缺乏对混合颜料分析的解决方案。
技术实现思路
本专利技术的目的是提供一种彩绘类文物混合颜料的分析方法，该方法能够永久留存彩绘类文物的颜料数字化信息，有助于对文物的数字化虚拟修复，并有效缩短文物实体的修复时间，对彩绘类文物的数字化保护与修复具有较大的经济效益和社会效益。本专利技术的目的是通过以下技术方案实现的：一种彩绘类文物混合颜料的分析方法，所述方法包括：步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；步骤3、将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；步骤4、再利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得待处理彩绘类文物表面色彩组成颜料的种类和比例。由上述本专利技术提供的技术方案可以看出，上述方法能够永久留存彩绘类文物的颜料数字化信息，有助于对文物的数字化虚拟修复，并有效缩短文物实体的修复时间，对彩绘类文物的数字化保护与修复具有较大的经济效益和社会效益。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。图1为本专利技术实施例提供的彩绘类文物混合颜料的分析方法流程示意图。具体实施方式下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术的保护范围。下面将结合附图对本专利技术实施例作进一步地详细描述，如图1所示为本专利技术实施例提供的彩绘类文物混合颜料的分析方法流程示意图，所述方法包括：步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；在该步骤中，宽波段连续反射光谱的覆盖范围为350nm-2500nm，采集环境为暗室，利用仪器自带的光源，测量多次取平均值，对同一种颜色采集至少二处以上浓淡不一区域的反射光谱；该宽波段连续反射光谱是由合成颜色的多种纯净颜料的光谱反射混合而成。步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；这里，非负矩阵分解算法要求矩阵中所有元素均为非负，在不考虑误差的情况下，将矩阵表示为两个非负矩阵完全相乘的结果，正是由于分解的结果是非负的特点，因而在实际应用中具有明确的物理意义，由于高光谱数据是对地物反射率数据的一种记载，满足非负性的要求，从而可以将非负矩阵分解算法用来高光谱数据的端元提取中来。具体实现中，该非负矩阵分解算法NMF采用如下式子来表示：Xn×m＝An×rSr×m+En×m其中，为待分解的宽波段连续反射光谱的混合光谱矩阵；和为分解后的两个端元矩阵；En×m表示为误差矩阵；r值为端元个数；并且||E||尽可能小，这里x∈R+表示x≥0；在不考虑误差的影响，对于由一个已知量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端元的丰度；再将F1端元经过二阶求导之后的丰度结果作为各个种类颜料的初始丰度值，具体表达式为：可以看出，比值导数光谱法求解过程，简单清晰，不需要复杂的迭代运算，能够快速、准确的进行端元光谱丰度的反演。值得注意的是，本专利技术实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技本文档来自技高网...

【技术保护点】
1.一种彩绘类文物混合颜料的分析方法，其特征在于，所述方法包括：步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；步骤3、将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；步骤4、再利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得该待处理彩绘类文物表面色彩组成颜料的种类和比例。

【技术特征摘要】
1.一种彩绘类文物混合颜料的分析方法，其特征在于，所述方法包括：步骤1、利用高光谱成像和地物光谱仪采集待处理彩绘类文物表面感兴趣区域的宽波段连续反射光谱；步骤2、将该宽波段连续反射光谱视为多源混合信息，利用非负矩阵分解算法实现相同颜色区域的颜料分解，得到不同颜料的反射光谱曲线；步骤3、将所述不同颜料的反射光谱曲线与颜料光谱库中的光谱进行匹配，得到混合颜料的种类；步骤4、再利用基于众数的比值导数法来获取各个种类颜料的丰度，进而获得该待处理彩绘类文物表面色彩组成颜料的种类和比例。2.根据权利要求1所述彩绘类文物混合颜料的分析方法，其特征在于，在步骤1中，所述宽波段连续反射光谱的覆盖范围为350nm-2500nm；且该宽波段连续反射光谱是由合成颜色的多种纯净颜料的光谱反射混合而成。3.根据权利要求1所述彩绘类文物混合颜料的分析方法，其特征在于，在步骤2中，所述非负矩阵分解算法采用如下式子来表示：Xn×m＝An×rSr×m+En×m其中，为待分解的宽波段连续反射光谱的混合光谱矩阵；和为分解后的两个端元矩阵；En×m表示为误差矩阵；r值为端元个数；进一步，通过定义目标函数和选择优化准则来求解两个端元矩阵A和S，所定义的目标函数包括最小化欧式距离和K-L散度，其中：所述最小化欧式距离表示为：所述K-L散度表示为：其中i(i＝1，2，…n)为混合光谱矩阵X的行号，其最大值是n，n是所测量的一条光谱曲...

【专利技术属性】
技术研发人员：吕书强，侯妙乐，尹琴丽，刘依依，刘冰，董庆豪，
申请(专利权)人：北京建筑大学，
类型：发明
国别省市：北京,11