基于对称样本块的图像修复方法技术

一种基于对称样本块的图像修复方法；过程包括对原破损图像的预处理、破损区域的分割和破损区域修复。对原破损图像进行预处理包括：对图像进行高斯平滑、图像灰度化处理；而后进入破损区域分割与破损区域修复；破损区域修复是一个反复循环的过程，包括获得破损区域的边界点、获取破损区域边界上优先权最高的点和最先修复的破损块，以及选择搜索最佳对称样本块方法；选择搜索最佳对称样本块方法包括任意方向搜索对称样本块或者八方向寻找对称样本块两种方法，两种方法都包括更新破损块的像素、更新破损区域的面积和判断破损区域的面积是否为零，是，修复结束；否，回到获得破损区域的边界点重复循环。修复实验说明了本发明专利技术方法的有效性。

【技术实现步骤摘要】
基于对称样本块的图像修复方法
本专利技术属于数字图像处理与数字图像修复领域，具体涉及基于对称样本块的图像修复方法。
技术介绍
自从2000年“图像修复”概念引入数字图像处理后，短短十多年，数字图像修复得到广泛的应用，主要包括：古画数字复原、恢复破损的照片和影片、图像中的文字和对象移除等。所谓数字图像修复就是对图像上缺损区域进行信息填充的过程，其目的就是为了对有缺损的图像进行恢复，并且视觉上无法察觉到图像曾经缺损或已被修复的效果。图像修复技术中，一类是基于偏微分方程的修复算法或者基于几何图像模型的变分修复技术，利用待修补区域的边缘信息，修复小尺度缺损问题。另一类是用于填充图像中的大块丢失信息的图像补全技术，包含两种方法：一种是基于分解的方法，将图像分解为结构和纹理部分，其中结构部分用图像修复算法修复，纹理部分用纹理合成的方法填充。还有一种是基于块的纹理合成技术填充丢失的信息，它是从待修补区域边界上选取一个像素点，并以该点为中心，根据图像的纹理特征，选取大小合适的纹理块，然后在待修复块的周围寻找与之最相似的纹理匹配块来代替。然而图像修复仍然是一件困难的事情。现有的数字图像修复技术还不能有效修复大量的对称图像，如唐卡图像、壁画等，很多图像有一定的对称结构，图像中一些关键块的缺损，如一幅唐卡中释迦摩尼佛的一边耳部丢失，若仅利用该图像中耳部周围信息修复，则无法得到所需要的结果。
技术实现思路
为了利用被修复图像中所具备的对称信息，来修复破损区域，本专利技术提供一种基于对称样本块的图像修复方法，它是利用破损区域在图像上某一对称方向的信息填充破损区域的一种有效的数字化图像修复方法。为了实现本专利技术目的，本专利技术采取的技术方案为：一种基于对称样本块的图像修复方法，包括对原破损图像进行预处理、破损区域分割与破损区域修复，具体步骤如下：a.对原破损图像进行预处理a.1对图像进行高斯平滑因为高斯平滑对边缘保持较好，通常用如下5×5或3×3的高斯平滑模板掩模图像：a.2图像灰度化处理采用平均值法将彩色图像灰度化，即求出图像每个像素点的R、G、B三个分量像素值的平均值，并生成其灰度图像；b.破损区域分割常用的图像分割方法包括阈值分割法、边界分割法和区域增长法等分割方法；本专利技术采用区域增长的分割方法，其过程为：b.1.采用人机交互的方法，在被修复图像的破损区内任意选取一点作为种子点坐标(xr,yr)，并在灰度图像上记录下其像素值pr(xr,yr)；b.2.创建与原图同样大小的一个二值模板图像mask，将模板图像mask中像素值pr(xr,yr)设为1，其余像素设置为0；b.3.在灰度图像上通过算法1，得出模板图像mask里像素值为1的点集合，就对应原图像的破损区域，也就是分割出了原图像的破损区域；算法1包括的步骤：(1)创建一个堆栈；(2)获取种子点(xr,yr)及其像素值pr(xr,yr)；(3)以种子点(xr,yr)为中心，计算其像素值pr(xr,yr)与八邻域像素值pi(xi,yi)之差，如果|pr(xr,yr)-pi(xi,yi)|＜M，(i＝1,2,3,4,5,6,7,8)，M的取值根据实验确定，一般取10，并且mask相对应位置上像素值为0时，则将点(xi,yi)压入堆栈，mask相应位置上的像素值设为1；(4)判断堆栈是否为空，如果不为空则从堆栈中取出一个像素，把它当作点(xr,yr)，转回步骤(3)，否则转向步骤(5)；(5)结束；c.破损区域修复c.1获得破损区域的边界点c.1.1创建一个队列，用来存储破损区域边界点；c.1.2获得破损区域边界点：对二值模板图像mask，其中像素值为1的区域为原图破损区域，遍历整幅mask图像，如果像素点的像素值为1，并且它八个相邻像素值至少有一个是0，就将像素值为1的像素点存入队列中，最后得到的队列就是破损区域的边界点；c.2获取破损区域边界上优先权最高的点和最先修复的破损块图2所示为边界上优先权最高点和最先修复的破损块计算示意图，Ι代表全图，Ω表示破损区域，δΩ表示破损区域Ω的边界，Φ表示图像未破损的部分，p表示破损区域边界δΩ上的点，np是点p处的单位法向量，表示p点的等照线方向，Ψp是边界上以p点为中心、包含了Ω内的点和Φ内点的破损块；破损区域边界上优先权最高的点和以该点为中心的最先修复破损块计算过程如下：c.2.1计算边界上点p的优先权，由式(1)计算，P(p)＝C(p)D(p)(1)式(1)中：C(p)表示p点的置信因子，D(p)表示p点的数据因子，C(p)由式(2)计算，式(2)中：|Ψp|是Ψp的面积，图像中所有点的置信因子按照式(3)初始化，式(1)中D(p)通过式(4)计算，其中α是归一化因子，是图像灰度级的最大值；c.2.2计算边界上所有点的优先权，并比较值的大小从而获得优先权最高的点，记为p0；c.2.3获取以p0为中心的最先修复的块这个块的大小根据破损区域周围的纹理结构，选择与纹理结构大小相当的块，采用人机交互的方式，块的边长可以是3到99之间的奇数，即3×3到99×99大小的块，一般选择9×9、11×11到33×33；而后进入步骤c.3；如果图像具有左右、上下、右上与左下、左上与右下的八方向对称的情况之一，选择基于八方向的对称样本块修复方法，否则选择在任意方向搜索对称样本块的修复方法；实际上，任意方向搜索对称样本块修复的方法具有普遍性，完全可以取代基于八方向对称样本块的修复；但是如果一幅图像左右对称，左侧有破损区域，基于八方向对称样本块的修复速度更快，而任意方向搜索最佳对称样本块修复的方法则相对耗时；所以进入步骤c.3后，可以选择在任意方向搜索对称样本块或基于八方向的对称样本块的图像修复方式，一旦选定一种方法，则按照选定的这种方法循环修复直到修复结束；c.3在任意方向搜索最佳对称样本块或者在八方向寻找最佳对称样本块的图像修复c.3.1在任意方向寻找最佳对称样本块c.3.1.1在I-Ω中以任意点p1为中心的样本块的选取如图3所示，p0是最先修复块的中心点，p1是待寻找样本块的中心点，获得p1和p0连线与水平方向夹角θ，利用式(5)获得p1旋转到与p0同一水平线位置上p2的坐标：式(5)中：x0、x1和x2分别是p0、p1和p2的横坐标，y0、y1和y2分别是p0、p1和p2的纵坐标；获得以p2为中心、大小与相同并与成水平对称的样本块以3×3为例，破损块与样本块的对应关系如图4与图5所示，图4是破损块图5表示的水平对称样本块其中数字表示两个块对称位置；再利用式(6)获得样本块反方向旋转θ角的位置：式(6)中：(x,y)表示反向旋转之前中任意点的坐标，(x’,y’)表示旋转θ角后的中对应点的坐标；旋转后，与点的对应关系和与的对应关系相同，即得到以p1为中心的对称样本块；c.3.1.2在I-Ω搜索的最佳对称样本块，由式(7)计算，式(7)中：表示最佳对称样本块，是以p1为中心的样本块，表示最先待修复块与对称样本块之间的相似性度量，用式(8)计算，式(8)中：xij表示中点的像素值，yij表示中点的像素值，块的大小决定m的值，若选择块的大小为3×3，则m的值为3；c.3.1.3破损块的更新c.3.1.3.1按照图4与图5中的对应关系，将破损块中的破损像素值用最佳本文档来自技高网...

【技术保护点】
一种基于对称样本块的图像修复方法，其特征在于包括对原破损图像进行预处理、破损区域分割与破损区域修复，具体步骤如下：a.对原破损图像进行预处理a.1对图像进行高斯平滑用如下5×5或3×3的高斯平滑模板掩模图像：1273×1474141626164726412674162616414741116×121242121]]>a.2图像灰度化处理采用平均值法将彩色图像灰度化，即求出图像每个像素点的R、G、B三个分量像素值的平均值，并生成其灰度图像；b.破损区域分割采用区域增长的分割方法，其过程为：b.1.采用人机交互的方法，在被修复图像的破损区内任意选取一点作为种子点坐标(xr,yr)，并在灰度图像上记录下其像素值pr(xr,yr)；b.2.创建与原图同样大小的一个二值模板图像mask，将模板图像mask中像素值pr(xr,yr)设为1，其余像素设置为0；b.3.在灰度图像上通过算法1，得出模板图像mask里像素值为1的点集合，就对应原图像的破损区域，也就是分割出了原图像的破损区域；算法1包括步骤：(1)创建一个堆栈；(2)获取种子点(xr,yr)及其像素值pr(xr,yr)；(3)以种子点(xr,yr)为中心，计算其像素值pr(xr,yr)与八邻域像素值pi(xi,yi)之差，如果|pr(xr,yr)‑pi(xi,yi)|＜M，(i＝1,2,3,4,5,6,7,8)，M的取值根据实验确定，一般取10，并且mask相对应位置上像素值为0时，则将点(xi,yi)压入堆栈，mask相应位置上的像素值设为1；(4)判断堆栈是否为空，如果不为空则从堆栈中取出一个像素，把它当作点(xr,yr)，转回步骤(3)，否则转向步骤(5)；(5)结束；c.破损区域修复c.1获得破损区域的边界点c.1.1创建一个队列，用来存储破损区域边界点；c.1.2获得破损区域边界点：对二值模板图像mask，其中像素值为1的区域为原图破损区域，遍历整幅mask图像，如果像素点的像素值为1，并且它八个相邻像素值至少有一个是0，就将像素值为1的像素点存入队列中，最后得到的队列就是破损区域的边界点；c.2获取破损区域边界上优先权最高的点和最先修复的破损块Ι代表全图，Ω表示破损区域，δΩ表示破损区域Ω的边界，Φ表示图像未破损的部分，p表示破损区域边界δΩ上的点，np是点p处的单位法向量，表示p点的等照线方向，Ψp是边界上以p点为中心、包含了Ω内的点和Φ内点的破损块；破损区域边界上优先权最高的点和以该点为中心的最先修复破损块计算过程如下：c.2.1计算边界上点p的优先权，由式(1)计算，P(p)＝C(p)D(p)  (1)式(1)中：C(p)表示p点的置信因子，D(p)表示p点的数据因子，C(p)由式(2)计算，C(p)=Σq∈Ψp∩(I-Ω)C(q)|Ψp|---(2)]]>式(2)中：|Ψp|是Ψp的面积，图像中所有点的置信因子按照式(3)初始化，C(k)=0∀k∈Ω1∀k∈I-Ω---(3)]]>式(1)中D(p)通过式(4)计算，D(p)=|▿Ip⊥·np|α---(4)]]>其中α是归一化因子，是图像灰度级的最大值；c.2.2计算边界上所有点的优先权，并比较值的大小从而获得优先权最高的点，记为p0；c.2.3获取以p0为中心的最先修复的块这个块的大小根据破损区域周围的纹理结构，选择与纹理结构大小相当的块，采用人机交互的方式，块的边长是3到99之间的奇数，即3×3到99×99大小的块；而后进入步骤c.3；如果图像具有左右、上下、右上与左下、左上与右下的八方向对称的情况之一，选择基于八方向的对称样本块修复方法，否则选择在任意方向搜索对称样本块的修复方法；c.3在任意方向搜索最佳对称样本块或者在八方向寻找最佳对称样本块的图像修复c.3.1在任意方向寻找最佳对称样本块c.3.1.1在I‑Ω中以任意点p1为中心的样本块的选取p0是最先修复块的中心点，p1是待寻找样本块的中心点，获得p1和p0连线与水平方向夹角θ，利用式(5)获得p1旋转到与p0同一水平线位置上p2的坐标：x2=(x1-x0)cosθ+(y1-y0)sinθ+x0y2=-(x1-x0)sinθ+(y1-y0)cosθ+y0---(5)]]>式(5)中：x0、x1和x2分别是p0、p1和p2的横坐标，y0、y1和y2分别是p0、p1和p2的纵坐标；获得以p2为中心、大小与相同并与成水平对称的样本块以3×3为例，破损块与样本块的对应关系如图4与图5所示，...

【技术特征摘要】
1.一种基于对称样本块的图像修复方法，其特征在于包括对原破损图像进行预处理、破损区域分割与破损区域修复，具体步骤如下：a.对原破损图像进行预处理a.1对图像进行高斯平滑用如下5×5或3×3的高斯平滑模板掩模图像：a.2图像灰度化处理采用平均值法将彩色图像灰度化，即求出图像每个像素点的R、G、B三个分量像素值的平均值，并生成其灰度图像；b.破损区域分割采用区域增长的分割方法，其过程为：b.1.采用人机交互的方法，在被修复图像的破损区内任意选取一点作为种子点坐标(xr,yr)，并在灰度图像上记录下其像素值pr(xr,yr)；b.2.创建与原图同样大小的一个二值模板图像mask，将模板图像mask中像素值pr(xr,yr)设为1，其余像素设置为0；b.3.在灰度图像上通过算法1，得出模板图像mask里像素值为1的点集合，就对应原图像的破损区域，也就是分割出了原图像的破损区域；算法1包括步骤：(1)创建一个堆栈；(2)获取种子点(xr,yr)及其像素值pr(xr,yr)；(3)以种子点(xr,yr)为中心，计算其像素值pr(xr,yr)与八邻域像素值pi(xi,yi)之差，如果|pr(xr,yr)-pi(xi,yi)|＜M，(i＝1,2,3,4,5,6,7,8)，M＝10，并且mask相对应位置上像素值为0时，则将点(xi,yi)压入堆栈，mask相应位置上的像素值设为1；(4)判断堆栈是否为空，如果不为空则从堆栈中取出一个像素，把它当作点(xr,yr)，转回步骤(3)，否则转向步骤(5)；(5)结束；c.破损区域修复c.1获得破损区域的边界点c.1.1创建一个队列，用来存储破损区域边界点；c.1.2获得破损区域边界点：对二值模板图像mask，其中像素值为1的区域为原图破损区域，遍历整幅mask图像，如果像素点的像素值为1，并且它八个相邻像素值至少有一个是0，就将像素值为1的像素点存入队列中，最后得到的队列就是破损区域的边界点；c.2获取破损区域边界上优先权最高的点和最先修复的破损块Ι代表全图，Ω表示破损区域，δΩ表示破损区域Ω的边界，Φ表示图像未破损的部分，p表示破损区域边界δΩ上的点，np是点p处的单位法向量，表示p点的等照线方向，Ψp是边界上以p点为中心、包含了Ω内的点和Φ内点的破损块；破损区域边界上优先权最高的点和以该点为中心的最先修复破损块计算过程如下：c.2.1计算边界上点p的优先权，由式(1)计算，P(p)＝C(p)D(p)(1)式(1)中：C(p)表示p点的置信因子，D(p)表示p点的数据因子，C(p)由式(2)计算，式(2)中：|Ψp|是Ψp的面积，图像中所有点的置信因子按照式(3)初始化，式(1)中D(p)通过式(4)计算，其中α是归一化因子，是图像灰度级的最大值；c.2.2计算边界上所有点的优先权，并比较值的大小从而获得优先权最高的点，记为p0；c.2.3获取以p0为中心的最先修复的块这个块的大小根据破损区域周围的纹理结构，选择与纹理结构大小相当的块，采用人机交互的方式，块的边长是3到99之间的奇数，即3×3到99×99大小的块；而后进入步骤c.3；如果图像具有左右、上下、右上与左下、左上与右下的八方向对称的情况之一，选择基于八方向的对称样本块修复方法，否则选择在任意方向搜索对称样本块的修复方法；c.3在任意方向搜索最佳对称样本块或者在八方向寻找最佳对称样本块的图像修复c.3.1在任意方向寻找最佳对称样本块c.3.1.1在I-Ω中以任意点p1为中心的样本块的选取p0是最先修复块的中心点...

【专利技术属性】
技术研发人员：王维兰，贾艳军，
申请(专利权)人：西北民族大学，
类型：发明
国别省市：甘肃;62