一种基于非局部几何相似性的图像插值方法
【技术领域】
[0001] 本发明设及一种图像计算方法，尤其设及一种基于非局部几何相似性的图像插值 方法。
【背景技术】
[0002] 图像的分辨率代表着图像所包含细节内容的多少，分辨率越高，图像提供的细节 就越多，反之，分辨率越低，图像提供的细节就越少。因此，为了获取更多的细节内容，让图 像变得更清楚，就需要增大图像的尺寸，即提高图像的分辨率，而提高图像分辨率的一种有 效途径就是对图像进行插值。图像的插值技术已经在军事、医学、安防等领域得到了广泛的 应用。图像插值的实质是通过未知像素点周围若干已知像素点灰度值的组合来预测未知像 素点的灰，传统的插值方法包括双S次插值，参见参考文献"Cubic convolution inte;rpolation for digital image processing";基于边缘指导的插值算法，参见参考文 献"New edge-directed interpolation";基于自回归模型的插值算法，参见参考文献 ('Image interpolation via regularized local linear regression"，；基于自回归模型 白勺插值算'法，参见参考文南犬（'Image Interpolation Via Regularized Local Linear Regression"。在运些传统的图像插值方法中，双=次插值算法可W使生成的像素灰度值延 续原图像灰度值变化的趋势，从而使放大图像呈现出自然的平滑，特别是图像内部的像素 点呈现连续变化时，运种算法更为有效。但当图像内部的像素灰度值出现突变时，使用运种 算法在放大后的图像中则无法有效控制和体现出像素灰度值的突变。基于边缘指导的插值 方法，主要是利用低分辨率图像块和对应的高分辨率图像块之间存在的几何相似性，从而 利用低分辨率图像的统计量来代替插值模型中的高分辨率统计量，进而求得插值因子。基 于图像局部几何相似性的插值算法在图像的局部区域存在几何相似性时，可W有效地产生 未知像素点的灰度值，但当图像的局部区域内不存在几何相似性时，运种算法就会失效。在 基于自回归模型的算法，插值模型利用局部统计特征，通过求解有约束的最小二成问题，获 取均方误差最小条件小的插值因子。但运种算法因为仅仅使用了图像内部的局部统计特 征，而造成了算法不够高效。
【发明内容】
[0003] 本发明的目的就在于为了解决上述问题而提供一种基于非局部几何相似性的图 像插值方法。
[0004] 本发明通过W下技术方案来实现上述目的：
[0005] 本发明包括W下步骤：
[0006] (1)图像预处理:首先，取待插值图像的左边界和右边界分别向左和向右重复地扩 充两列;然后，取扩充后图像的上边界和下边界分别向上和向下重复地扩充两行；
[0007] (2)产生参考列向量:将步骤1产生的图像内每一个大小为7X7的正方形图像块内 的所有像素点按列依次取出，组成一个大小为49 X 1的参考列向量；
[0008] (3)初始插值:用传统的双=次插值法对待插值的图像进行初始化插值，得到一个 放大两倍的初始插值图像；
[0009] (4)第一类插值:对步骤3中得到的初始插值图像，在位于偶数行和偶数列的位置 上进行第一类插值:首先，W该位置为中屯、，取出一个大小为7X7的待插值正方形图像块； 接着，将该图像块内的全部像素点按列依次取出，组成一个大小为49X1待插值列向量;其 次，将待插值列向量和每一个参考列向量一起代入公式（1)，计算待插值图像块与每一个参 考图像块之间的几何相似度，并找出其中具有最小几何相似度的前100个参考列向量;然 后，将运100个参考列向量中每个列向量的第17、19、31、33和25个元素^及待插值列向量的 第17、19、31、33和25个元素一起代入公式(2)，进一步计算几何相似度，并找出其中具有最 小几何相似度的前60个参考列向量;最后，将运60个参考列向量中每个列向量的第17、19、 31、33和25个元素都代入公式(3)中，计算出4个插值加权因子，并将得到的加权因子和待插 值列向量的第17、19、31、和33个元素一起代入公式(6)，计算出插值灰度值；
[0010] 公式(1)几何相似g
[OOW 式中：Pi表示由待插值图像块组成的列向量中的元素，riW表示第k个参考列向量 中的元素；
[001^ 公式(2)几何相似度
[OOU]式中：义1、枯、义3、义4和(3表示从待插值列向量中取出的5个像素点的灰度值，曰1^、曰2 W、a3^、a4W和表示从第k个参考列向量中取出的5个像素点的灰度值，0由公式(3)计算 得到；
[0014]公式（3
[001 引式中：戈.=[々.,抱,從JW=[巧("，媒。尸'皆，，6脚].，.X1、X2、X 3、X^C 表示从待插值列向量中取出的5个像素点的灰度值，U,"，、皆\ :冷> 和bW表示从 第k个参考列向量中取出的5个像素点的灰度值，符号T表示矩阵的转置操作；
[0016] 公式(-
[0017] 式中：O由匹配的参考列向量依次按列摆放组成，B是由所有匹配列向量的第25个 元素组成的列向量，X= [XI，X2，X3，X4] ,Qi = X-Gi,
[001引 Gi二[幻,W-A…，却'-々w，為护'、，幻f知')]?'，对、又2、村和勒表示从待插值列向 量中取出的4个像素点的灰度值，af、af、af、af和表示从第k个参考列向量中 取出的5个像素点的灰度值，A由用户自己定义，符号T表示矩阵的转置操作;Sk由公式(5)计 算得到;化由公式(3)计算得到；
[0019]公式（5):
[0020] 式中
表示从待插值列向量中取出的5个像素点的灰度值，af、皆\媒>和b?表示从第 k个参考列向量中取出的5个像素点的灰度值;Pk由公式（3)计算得到；A是正则化因子，
O由用户自己定义，6 = 2.7183;
[0021] 公式（6) :Xi,j = [Wl W2 W3 W4] ? [XI X2 X3 X4]T;
[002^ 式中：W1、W2、W3和W4是插值因子，X1、X2、X3和X4表示从待插值列向量中取出的与待插 值位置相邻的4个已知像素点的灰度值，符号T表示矩阵的转置操作。
[0023] (5)第二类插值:对步骤4中得到的插值图像，在位于偶数行和奇数列上的位置上 进行第二类插值:首先，W该位置为中屯、，取出一个大小为7X7的待插值正方形图像块;接 着，将该图像块内的全部像素点按列依次取出，组成一个大小为49X1待插值列向量;其次， 将待插值列向量和每一个参考列向量一起代入公式（1)，计算待插值图像块与每一个参考 图像块之间的几何相似度，并找出其中具有最小几何相似度的前100个参考列向量;然后， 将运100个参考列向量中每个列向量的第18、24、26、32和25个元素 W及待插值列向量的第 18、24、26、32和25个元素一起代入公式(2)，进一步计算几何相似度，并找出其中具有最小 几何相似度的前60个参考列向量;最后，将运60个参考列向量中每个列向量的第18、24、26、 32和25个元素都代入公式(3)中，计算出4个插值加权因子，并将得到的加权因子和待插值 列向量的第18、24、26、和32个元素一起代入公式(6)，计算出插值灰度值；
[0024] (6)第=类插值:对步骤5中得到的插值图像，在位于奇数行和偶数列上的像素点 进行第二类插值:首先，W该位置为中屯、，取出一个大小为7X7的待插值正方形图像块;接 着，将该图像块内的全部像素点按列依次取出，组成一个大小为49X1待插值列向量;其次， 将待插值列向量和每一个参考列向量一起代入公式（1)，计算待插值图像块与每一个参考 图像块之间的几何相似度，并找出其中具有最小几何相似度的前100个参考列向量;然后， 将运100个参考列向量中每个列向量的第18、24、26、32和25个元素 W及待插值列向量的第 18、24、26、32和25个元素一起代入公式(2)，进一步计算几何相似度，并找出其中具有最小 几何相似度的前60个参考列向量;最后，将运60个参考列向量中每个列向量的第18、24、26、 32和25个元素都代入公式(3)中，计算出4个插值加权因子，并将得到的加权因子和待插值 列向量的第18、24、26和32个元素一起代入公式(6)，计算出插值灰度值。
[0025] 本发明的有益效果在于：
[0026] 本发明是一种基于非局部几何相似性的图像插值方法，与现有技术相比，本发明 首先将插值模型中所要使用的统计特征从图像的局部区域扩展到了整个图像内部，并引入 了一种改进的欧式距离来衡量图像内部所有图像块之间的几何相似性，通过两步捜索法来 寻找相似的图像块W构建插值模型，同时，所有用来构建插值模型的相似图像块的中屯、点 被应用于插值模型来优化所产生的插值因子，最终所构建出的插值算法，效果更优，稳定性 更强。
【具体实施方式】
[0027]下面对本发明作进一步说明：
[002引本发明包括W下步骤：
[0029] (1)图像预处理:首先，取待插值图像的左边界和右边界分别向左和向右重复地扩 充两列;然后，取扩充后图像的上边界和下边界分别向上和向下重复地扩充两行；
[0030] (2)产生参考列向量:将步骤1产生的图像内每一个大小为7X7的正方形图像块内 的所有像素点按列依次取出，组成一个大小为49 X 1的参考列向量；
[0031] (3)初始插值:用传统的双=次插值法对待插值的图像进行初始化插值，得到一个 放大两倍的初始插值图像；
[0032] (4)第一类插值:对步骤3中得到的初始插值图像，在位于偶数行和偶数列的位置 上进行第一类插值:首先，W该位置为中屯、，取出一个大小为7X7的待插值正方形图像块； 接着，将该图像块内的全部像素点按列依次取出，组成一个大小为49X1待插值列向量;其 次，将待插值列向量和每一个参考列向量一起代入公式（1)，计算待插值图像块与每一个参 考图像块之间的几何相似度，并找出其中具有最小几何相似度的前100个参考列向量;然 后，将运100个参考列向量中每个列向量的第17、19、31、33和25个元素^及待插值列向量的 第17、19、31、33和25个元素一起代入公式(2)，进一步计算几何相似度，并找出其中具有最 小几何相似度的前60个参考列向量;最后，将运60个参考列向量中每个列向量的第17、19、 31、33和25个元素都代入公式(3)中，计算出4个插值加权因子，并将得到的加权因子和待插 值列向量的第17、19、31、和33个元素一起代入公式(6)，计算出插值灰度值；
[0033] 公式（1)几何相似度