摘要：针对分段迭代曲线拟合存在的重建区域轮廓不连续、重建区域尺寸有误差等问题，提出了一种基于融合细分的纹理图像重构模型。首先提取原始图像的分割区域，经过轮廓跟踪与下采样得到区域形状的特征向量;然后利用三重逼近与三重插值统一的融合细分方法，重建区域轮廓曲线;最后合成区域纹理，得到纹理图像重构结果。在多幅自然场景图像上进行实验验证，并给出相应的实验结果和分析。实验结果表明，所提模型正确有效，具有和人类视觉特性相符合的重构结果; 所提算法能够减少图像重建时的处理时间，并在图像质量主观评价指标上明显优于多区域图像重建算法。
关键词：纹理图像重构;融合细分;纹理合成;图像分割
一、引言
随着数字多媒体技术的发展，人们对高效的图像信息处理提出了更高的要求。现有的基于像素/块的图像处理技术忽略了图像的层级结构，无法直接用于内容分析;而基于对象的又难以满足图像处理在通用性方面的需求[1]。因此，如何找到一种更加有效的图像表征方法一直是图像处理领域的研究热点与难点问题。
纹理是表达图像内容的一个非常重要的属性，它广泛存在于各类图像中。纹理图像通常构成图像或视频的静止背景。图像中的纹理分为两大类：不重要主观细节纹理和重要主观细节纹理[2]。由于人类视觉系统固有的缺陷，人眼对平滑区域的敏感性远高于纹理密集区域，图像中的纹理通常是人眼不关注的那部分内容，因此纹理细节的变化不会影响对原始纹理的主观理解[3]。
近年来，基于样图的纹理合成在图像修复、压缩编码、纹理传输等方面有着广泛的应用[4-6]。Efros等[7]提出一种计算较简单的Image Quilting算法，通过计算纹理重叠区域的累积误差和最小误差路径进行纹理拼接。而旋转的Wang Tiles 纹理合成算法[8]能够克服Wang Tiles 存在的样图利用不完全、切割路径非最优、中心和拐角区域不匹配等缺点。
图像中的线结构是指用来定义目标形状的轮廓或划分区域的边界，是图像的形状特征表达[9-10]。早期的形状描述方法使用二进制图像，基于二进制边缘的方法有很多，例如多边形近似、曲率的频域表示等。Zhang等[11]系统阐述了两类图像形状的表示算法：一类是基于轮廓，另一类是基于区域。每一类形状表示算法可分为结构方法和全局方法，并适用于空域和变换域。
但是，利用分段迭代曲线拟合重建的区域轮廓曲线是不连续的，需使用膨胀算法进行修复，从而造成重建区域的尺寸存在误差，因此会引起重建质量差的问题。针对上述问题，提出利用三重逼近与三重插值统一的融合细分方法，重建区域轮廓曲线;将合成的区域纹理填充到重建的区域轮廓曲线中，从而得到纹理图像重构结果。实验结果表明本文算法重构的图像质量良好，图像质量的主观评价指标明显优于多区域图像重建算法。
图1是本文提出的基于融合细分的纹理图像重构模型结构流程，主要包括预处理、融合细分和重构纹理图像3个阶段。
预处理阶段首先利用图像分割与区域融合方法获得原始纹理图像的多个区域及其索引号;然后通过轮廓跟踪与下采样，获得代表区域轮廓结构特征的有序序列。
融合细分阶段利用三重逼近与三重插值统一的融合细分方法，重建区域轮廓曲线。
重构纹理图像阶段首先，得到包含区域纹理局部与全局特征的区域纹理样本;然后，利用基于样图的纹理合成算法合成区域纹理;最后，将合成的区域纹理填充到重建的区域轮廓曲线中，实现纹理图像重构。
1.1预处理
首先，原始图像被分割为许多包含同质颜色与纹理特征的小分割区域[13-14]，同时考虑到纹理样本选择对分割结果是敏感的，空间紧邻的小分割区域必须合并成较大的区域，每个区域获得一个唯一的索引号。
然后，利用轮廓跟踪算法[10-11]得到一个由区域边界曲线上的像素点组成的有序序列，按照随机间隔进行有序抽样即下采样提取，得到代表区域轮廓形状的特征向量(xs，ys)，其中s是向量元素在区域轮廓曲线上的位置序号。
为了得到插值与逼近统一的融合细分方法，本文用新的角度观察三重逼近细分。
然后定义位移i如下所示：
Pk+13i=Pki-4ki
根据式(4)定义的位移算子，得到三重插值细分新的表示：
Pk+13i=Pki
然后，利用权值参数(01)得到三重逼近与三重插值统一的融合细分。三重逼近细分规则和三重插值细分规则统一的表示：
Pk+13i=Pki-4ki
显然，式(3)和式(5)均为式(6)的特殊情况。当=1时，式(6)代表三重逼近细分;当=0时，式(6)代表三重插值细分;当01时，式(6)可生成介于三重逼近细分和三重插值细分之间的细分曲线。
最后，进一步整理式(6)，得到三重融合细分方法新的表示：
Pk+13i=a0Pki-1+a1Pki+a0Pki+1
1.3重构纹理图像
首先，利用自回归统计分析模型对区域内部纹理进行分析;然后，选择包含区域纹理局部与全局特征信息的区域纹理样本，合成区域纹理[7-8];最后，把合成的区域纹理填充到重建的区域轮廓曲线中，得到纹理图像重构结果。具体的过程示例如图3所示。
二、实验结果与分析
由实验结果(图4～8)可知，本文的纹理图像重构质量优于对比算法。影响纹理图像重构质量的因素主要有两个：
2) 重建区域尺寸。为了使重建的区域轮廓曲线是连续、封闭的，文献[12]使用了膨胀算法，造成重建区域的尺寸误差;原始区域尺寸越小，这种误差就越明显;本文中重建区域的尺寸是正确的。
本文提出的基于融合细分的纹理图像重构模型的目的并不是为了使重构纹理图像与原始纹理图像完全一致，而是使人眼尽可能看起来自然、不令人讨厌。本文采用主观质量评价方法――双刺激失真衡量阶梯(Double Stimulus Impairment Scale， DSIS)。挑选25名非图像专业的一年级学生作为观测者。表1是ITUR五分制评分等级。平均意见分(Mean Opinion Score，MOS)代表观测者评分的平均水平;标准差用于计算主观质量评分结果的可信范围即95%置信区间[18]。
三、结语
本文提出一种基于融合细分的纹理图像重构模型。该模型首先得到原始图像的分割区域，进而提取代表轮廓形状的特征向量;然后利用融合细分重建区域轮廓，最后将合成的区域纹理填充到重建的区域轮廓曲线中，实现纹理图像重建。
与其他基于区域形状与纹理信息的图像重建算法相比，用本文给出的三重逼近与三重插值统一的融合细分方法重建的区域轮廓，能更准确有效地表征区域轮廓形状。对于包含大量纹理的原始图像，本文算法能得到好的重建结果，但对于含有非纹理细节信息的原始图像，重建质量不理想。因此，如何能够依据图像自身特点，设计出更合理的重建方法，适用于更复杂的图像是今后进一步努力的方向。
参考文献：
[3]SUN K， YE L， YANG Y， et al. Image restoration using piecewise iterative curve fitting and texture synthesis[C]// Proceedings of the 4th International Conference on Intelligent Computing. Piscataway： IEEE， 2008： 1056-1063.
[5]BALLE J， STOJANOVIC A， OHM J R. Models for static and dynamic texture synthesis in image and video compression[J]. Journal of sel
ected Topics in Signal Processing， 2011， 5(7)： 1353-1365.
[6]ZHANG X， KIM Y J. Efficient texture synthesis using strict Wang Tiles[J]. Graphical Models， 2008， 70(3)： 43-56.
[7]EFROS A A， FREEMAN W T. Image quilting for texture synthesis and transfer[C]// Proceedings of the 28th Annual Conference on Computer Graphics and Interactive Techniques. New York： ACM， 2001： 341-346.
