《艳云双鹿》古画数字化修复研究

1、    艳云双鹿古画数字化修复研究    周杰段炼孺摘 要：古画作为华夏文明的重要载体，从根本上体现了我国传统艺术的活力和内涵。但许多古画由于年代久远等原因遭到损坏，古画的复原对历史文化、造型艺术具有极大的影响。为了更好地对古画进行保护，提升修复的效率和缩短修复耗时，文章将分类对古画形式和破坏情况进行研究，基于criminisi算法进行改进，插入了分散度项，确保信息最复杂的部分被先行修复，然后在相似度计算中引入协方差因子，以达到正确填充率。对艳云双鹿古画进行数字化修复的实验结果表明，该算法可以相对好地实现数字化修复在古画修复方面的研究。关键词：古画修复;

2、数字化修复;计算机图像处理;艳云双鹿古画不仅是中华灿烂文明的重要体现，而且绘画技法种类多样，历史悠久，具有极高的艺术价值，在历史、考古、美术、设计等领域有很重要的研究价值。如今古画修复对延长文物的“寿命”起到积极作用。但是由于时隔久远，保护力度不够，古代画作受到的破坏也是有差异的，导致价值损失严重。在流传的过程中有许多古代绘画作品出现各种损坏，如艳云双鹿由于在历史长河中保护不当，出现了糟朽断裂、脱色、褪色等问题。本文修复作品在破损程度上也是具有一定的典型性，所以选择艳云双鹿这幅画为代表进行修复。传统古画修复技术多为人为手工修复，效率低下，耗时较长，而数字化修复具有效率高、耗时短等优点，是未来的

3、修复趋势。从古画修复视角来看，古画修复本质为图像的数字化修复。目前国内外对于数字化图像的修复处在初步研究阶段。闫丽霞虚拟化地对古代画作运用多光谱技术进行修复研究。周平平等对古代画作中污渍部分运用高分光谱影像分类线性回归的方法进行修复。根据以上研究，数字化修复技术已经在古建、壁画、岩画等修复领域有了一定的应用，主要分为变分偏微分法、纹理合成法等方式。目前流行的criminisi算法仅仅参照了纹理信息获取图像样本，但研究发现大量的古画存在夸张、抽象、主观性强和色彩、线条等细节丰富多样的特殊性，导致出现结构性偏差的严重问题。因此考虑从实践出发，对 criminisi算法进行改进，运用图样数字化修复方

4、法针对艳云双鹿古画进行修复研究。1  criminisi 算法criminisi 算法在边缘上获取适当的点作为修复区域中央，用正常画作部分中的最相近部位代替被修复块，一直达到所有部位被修补完整。此算法主要包括3步：step1：以破损区域的边缘?上获取一点o为中点，把m×m 左右的面积o按顺序得到所有o的先行权，选择先行权值较大者为当前破损修复块。以点o为中点的图像块优先函数为：破损修复块中有用信息用置信度h（o）表示：o 表示破损块中总像素数;点o处结构的强弱反应由数据项  表示。表示点o的等照度矢量;no为在点o处的单位法向量;为归一化因子。step2：在已知区

5、域中通过相似度计算搜索最良匹配块r=arg min d（o，r ），其中，d（o，r ）表示o 于r中对应像素点 rgb 3 个通道的欧氏距离平方和;用r填充o，完成一次填充。step3：刷新置信度。等到一次图像破损块修复后，修复部位于其边缘同时变化，因此每当修复操作完成一步后，应当再次用公式（2）算出置信度项k（o）。重复step1step3，直至待修复区域完整为止。2  本文修复算法本文将从两个部分对 criminisi 算法进行改进，分别为相似性度量、优先权计算。第一，量区域中现存像素的疏散程度，用分散度项对其进行权衡。第二，优先权计算函数以达到填充正确的目的。第三，斟酌结构和

6、纹理信息相关性在原算法中匮乏的问题，首先在填充纹理时思考了其结构信息，其次在相似性计算式中插入协方差。2.1 优先权函数算法改进如圖1-a所示的修补区域不仅颜色变化程度较复杂，且须修复形状不规则。而图1-b所示修补区域不仅颜色变化程度较简单，且须修复形状较规则。在分布特点不足的条件下，仅仅考虑了已知像素的数目，用criminisi 算法用公式（1）计算优先权，有非常大的概率会使其进行错误填充。为了确保信息量最多的修复部分先被修复。本文增加了一个分散度项 w（o），古画中待修复的部分用o点作为中心点表示，作为现有图案中的已知像素，并且分散度于w值成正比，具体算法如下。输入：m×m待修复

7、块o。输出：分散度项 w（o）。step1：将o列向量化为zo。step2：zo中第x个像素的坐标值记为cx，计算 m×m个坐标值的平均值。step3：取o中像素r'的相对坐标为c'x。step4：计算o '的坐标方差=。step5：归一化坐标方差w（o）=，max、min分别为坐标方差的最大值和最小值。criminisi 算法中优先函数可以用置信度项乘以数据项获得，当修复到数据项 k（o）值为 0 则停止，导致结果出现错误。本文将优先函数改进为 ：在权值系数中、可调节权重，且+=1，破损区域修复在最佳效果时权值系数可在实验中得到。2.2 相似度函数改进cri

8、minisi 算法在计算最优匹配区域时，运用分化的图片修复计算方法可以补充对图像数据准备的不足;并且忽略了整体的结构和块间的结构一致性，只考虑块间相应像素的间隔。根据以上原因，本文算法将把协方差计算引入块间相似度计算公式，对破损古画进行整体修复。单通道的协方差计算公式为：其中，两个图案间相应的点在布局方向上的一致性，用均值的差值乘以像素点来表示，将破损区域的像素均值用y代替，匹配块的像素均值用z 代替。方向越一致，乘积结果为正，绝对值就越大，反之相反。这两个像素块的整体布局的一致性可由求和后的结果来体现;为归一化参数，灰度图像取值为 2442。首先归一化rgb各个通道，然后求协方差因子yz。本

9、文相似度计算改进后的公式为：2.3 算法步骤本文对古画的修复算法步 如下。输入：破损古画 y。输出：修复结果 z。step1：用公式（3）算出以破损边缘上o点为中点 m×m区域的优先权。step2：根据优先权值获取首批恢复区go。step3：通过公式（4）（5）的计算，在已知区域寻找最佳匹配块gr。step4：用gr替换go。step5：若存在未修复区域，转step1，否则，算法结束。3 实验结果及分析本文实验环境为：英特尔core i5-9400f  8 gb内存，64位操作系统，matlab r2014a 平台。实践中，像素块大小均为 9*9， 优先函数计算公式中的权值

10、系数、分别取为0.3、0.5、0.2。3.1 模拟破损古画的修复为分析本文算法的显著性，此次修复选不同古画中完整部分，对完好图样进行破损模拟处理后形成的破损图样样进行修复，表1为criminisi算法与本文算法通过修复结果图样的峰值信噪比（peak signal to noise ratio，psnr）对比，其中共有5类图样的修复结果，图2是其中的2幅样例的修复结果，绿色遮挡部分为将要修复部分。从而得出，本文算法比criminisi算法更优。3.2 真实破损古画的修复图3所示为截取艳云双鹿真实破损部分。图3中，左边古画图案修复部位主要在背景色彩丰富的位置，本文算法虽然也有少量缺失，但基本保持了色彩的流畅。中间古画的破损部位出现在鹿头边缘部分周边色差较大，从修复结果看，而本文算法的修复结果较好，但由于黑色区域将全部有效信息遮挡导致的，所以在部分区域出现了修复不足的现象。右边古画修复中既有色彩丰富部分，又有结构部分，结果显而易见，本文算法修复结果中背景及边缘纹理部分都有较好的修复结果。4 结语通过分析艳云双鹿古画纹理特點，本文提出一种古画修复算法。实验证明，相比现有的修复古画的方法，本文算法更适用于纸绘古画的修复。但是该算法对古画较小区域的脱色、脱落等修复问题适用，对于古画较大须修复区域，还有待提高。今后应对纸绘古画修