一种基于牛顿二阶插值的图像缩放方法技术

本发明专利技术专利提供了一种基于牛顿二阶插值的图像缩放方法，首先根据源图像和目标图像的分辨率大小计算出缩放倍率，将源图像和目标图像的像素位置进行坐标化；然后基于牛顿二阶插值算法对源图像进行水平方向上的缩放得到中间图像；再根据中间图像对垂直方向、45度方向和135度方向进行边缘位置判断，通过计算三个方向选取像素点的二阶差分来和给定阈值比较确定插值点是否处于边缘区域，最终选择对应的算法计算得到目标图像的像素值。本发明专利技术可以对源图像进行任意倍率的无极缩放，同时根据了插值点所在邻域内源像素的相关性对图像的边缘区域进行了保护，缩放后的图像视觉效果比传统的算法更佳，且算法复杂度低，利于硬件实现。

An image scaling method based on Newton second order interpolation

【技术实现步骤摘要】
一种基于牛顿二阶插值的图像缩放方法
本专利技术属于图像缩放
，特别涉及一种基于牛顿二阶插值的图像缩放方法，适用于FPGA实现对图像进行缩放。
技术介绍
图像是人类感知世界的视觉基础，是人类获取信息，表达信息的重要手段。图像缩放技术被广泛的应用于各类消费类电子产品、医学图像处理以及安防监控等领域中。图像是由离散的点阵构成，要对图像进行缩放处理，首先要提取原始图像中的像素数据，然后对像素数据进行处理得到目标图像。图像缩放的处理算法的种类很多，根据算法是否对边缘进行处理，可以大致分为两类。第一类是不考虑边缘信息的图像插值算法，如最邻近插值、双线性插值和双三次插值等一类经典的插值算法。这类算法在计算目标图像的像素时只通过其位置和邻域内的原始像素数据计算，不考虑边缘处的图像信息变化，缩放后效果不理想。其中，最邻近插值算法在实现过程中最为简单，但在处理后，目标图像会产生马赛克现象；双线性插值和双三次插值具有低通滤波效应，导致高频分量丢失，缩放后得到的目标图像会变得平滑。第二类是基于边缘的图像缩放算法，边缘检测算法有Sobel、Prewitt和Canny以及基于形态学的边缘检测算法。如江雯等提出基于Sobel算子的边缘插值算法，需要对图像边缘进行检测后存储，然后再根据边缘方向进行插值计算，其实现过于复杂。除了先采用边缘算子检测边缘再插值的方法，还有一类采用边缘自适应的算法，齐敏等人通过数据融合的方式对不同区域的插值点进行插值计算，算法通过多次归一化方法求系数，缩放效果并不理想，且计算过程繁杂。刘政林等人提出基于边缘的实时缩放算法，在边缘检测上考虑欠缺，导致边缘区域的缩放效果不佳。本领域需要提出缩放效果更好，计算复杂度低，更利于硬件实现的算法。
技术实现思路
基于上述分析，提出一种基于牛顿二阶插值的图像缩放方法，考虑到边缘的各个方向的相关性，处理结果相比传统的缩放算法，效果更好，计算复杂度低，且利于硬件实现。为了达到上述技术目的，本专利技术采用如下技术方案予以实现。一种基于牛顿二阶插值的图像缩放方法，包括以下步骤：步骤a，根据源图像和目标图像的分辨率大小确定输入图像和目标图像的缩放倍率，将源图像和目标图像的像素位置进行坐标化。步骤b，对坐标化后对原始图像进行水平缩放，根据目标像素点映射到原始图像中的位置选择邻域内相关性最大的一组源像素，通过牛顿二阶插值算法计算水平方向的目标像素值得到中间图像；步骤c，根据水平缩放后的中间图像，同样根据其映射到中间图像中的位置选择3x3邻域内的一组源像素，然后对垂直方向、45度方向和135度方向的各组源像素进行边缘判断，对处于边缘区域的插值点进行融合计算得到目标像素值，处于平缓区域根采用对应方向的源像素插值计算得到目标像素值，插值计算采用的是牛顿二阶插值算法，需要插值点邻域内相关性最大的三个源像素点参与计算。有益的技术效果由于传统缩放方法对图像的边缘区域缺少保护，现有的对边缘保护缩放算法存在一些算法复杂和效果不佳的现象；本专利技术方法计算量小且对边缘的保护更优，适合在硬件上实现。附图说明图1为本专利技术的一种基于牛顿二阶插值的图像缩放方法流程图；图2为目标像素水平方向插值示意图；图3为中间图像的插值示意图；图4为在FPGA平台下将512x512大小Lena原图放大到1024x1024大小输出的结果。具体实施方式参照图1，为本专利技术的一种基于牛顿二阶插值的图像缩放方法流程图；包括以下步骤：步骤a，根据源图像和目标图像的分辨率大小确定输入图像和目标图像的缩放倍率，将源图像和目标图像的像素位置进行坐标化。步骤b，对水平方向进行缩放只需考虑一维方向的缩放处理，假设对水平方向进行插值处理，将一行像素序列看成距离为1的等距节点，则目标像素点与原始图像在水平方向上的位置关系如图1所示。其中f是目标像素点在水平方向上映射到原始图像中的位置，f1、f2、f3和f4为f水平邻域内的四个原始像素点，Δx是f和f2之间的距离，大小为0≤Δx<1，且t＝1+Δx。选取二阶牛顿插值算法来计算f的像素值时，可以用f1、f2和f3三个源像素点作为一组像素或f2、f3和f4三个源像素点作为一组来计算目标像素f的像素值，根据插值点的位置选取一组像素计算，计算公式如式(1)所示。其中Δf1和Δf2分别是f1和f2的一阶差分，计算公式分别如式(2)和式(3)所示；Δ2f1和Δ2f2分别是f1和f2的二阶差分，计算公式如式(4)和式(5)所示。Δf1＝f2-f1(2)Δf2＝f3-f2(3)Δ2f1＝(f3-f2)-(f2-f1)(4)Δ2f2＝(f4-f3)-(f3-f2)(5)由以上公式可以对图像进行水平缩放，得到中间图像。步骤c，对中间图像45度、135度和垂直方向的进行缩放处理得到目标图像。其插值示意图如图2所示。图2中空心圆形和空心菱形分别表示插值点距离像素点f(i,j)小于0.5和大于0.5两种情况，实心圆点分别是插值点邻域内的各原始像素点。对水平缩放后的中间图像进行垂直方向上的缩放，假设图中空心圆为目标像素点所在位置，u的值小于等于0.5，根据水平方向的插值原理，这里在f(i,j)的3x3邻域内分别选取垂直方向、45度方向和135度方向各一组像素，计算其二阶差分分别为Δ2f1、Δ2f2和Δ2f3。设置阈值T分别与Δ2f1、Δ2f2和Δ2f3的绝对值进行比较，如果三者中某一个值小于T，则说明在该方向上为平缓区域，利用该方向上三个像素点进行二阶牛顿插值计算得到目标像素值，如果出现多个值都小于T，需给出三者顺序，先判断垂直方向，然后是45度方向，最后是135度方向。当目标像素在多个方向上处于平缓区域时，选择一组计算即可。对于三者都大于T时，则说明该目标像素在三个方向都不处于平缓区域，需利用三个方向的像素点进行计算得到目标像素值。方法如下：对于二阶差分越小的方向，其相关性越大所占权重应该越高，需对三个方向的二阶差分Δ2f1、Δ2f2和Δ2f3先归一化，然后处理得到各方向的权重系数，计算公式分别如式(6)、式(7)和式(8)所示。同时计算三个方向的二阶牛顿插值结果，为了简化计算，对于45度和135度方向上Δx的值直接选用垂直方向上Δx的值，即图3中u的值，且两源像素之间间距同样看作1处理，计算得到三个方向的结果分别为F1、F2和F3。再乘上各自的权重系数，因为三个方向权重系数之和为2，所以最后乘以1/2得到目标像素值，计算公式如式(9)所示。假设目标像素位置处于图2中空心菱形的位置，即u大于0.5，则选择f(i,j+1)像素点3x3邻域的像素点，采用同样的方法计算目标像素点的像素值。选用ALTERA公司StratixV系列的5SGXEA7N2F45C2NFPGA开发板实现上述算法，通过QuartusII15.0软件编译下载程序，将大小为512x512分辨率的Lena图像放大到1024x1024分辨率输出到显示器，如图4所示。为便于理解本专利技术技术效果起见，以下通过仿真实验来验证本专利技术的有效性：实验选用了House、Baboom、Items、Ruler和Bicycle图像，放大2倍再缩小还原，计算其峰值信噪比(PSNR)的值和边缘保护指数(EPI)的值。峰值信噪比(PSNR)是最常用的客观评价图像处理效果的量化指标，其计算本文档来自技高网...

【技术保护点】
1.一种基于牛顿二阶插值的图像缩放方法，其特征在于，包括以下步骤：步骤a，根据源图像和目标图像的分辨率大小确定输入图像和目标图像的缩放倍率，将源图像和目标图像的像素位置进行坐标化；步骤b，坐标化后，对原始图像进行水平缩放，根据目标像素点映射到原始图像中的位置选择邻域内相关性最大的一组源像素，通过牛顿二阶插值算法计算水平方向的目标像素值得到中间图像；步骤c，根据水平缩放后的中间图像，同样根据其映射到中间图像中的位置选择3x3邻域内的一组源像素，然后对垂直方向、45度方向和135度方向的各组源像素进行边缘判断，对处于边缘区域的插值点进行融合计算得到目标像素值，处于平缓区域根据对应方向的源像素插值计算得到目标像素值，插值计算采用的是牛顿二阶插值算法，需要插值点邻域内相关性最大的三个源像素点参与计算。

【技术特征摘要】
1.一种基于牛顿二阶插值的图像缩放方法，其特征在于，包括以下步骤：步骤a，根据源图像和目标图像的分辨率大小确定输入图像和目标图像的缩放倍率，将源图像和目标图像的像素位置进行坐标化；步骤b，坐标化后，对原始图像进行水平缩放，根据目标像素点映射到原始图像中的位置选择邻域内相关性最大的一组源像素，通过牛顿二阶插值算法计算水平方向的目标像素值得到中间图像；步骤c，根据水平缩放后的中间图像，同样根据其映射到中间图像中的位置选择3x3邻域内的一组源像素，然后对垂直方向、45度方向和135度方向的各组源像素进行边缘判断，对处于边缘区域的插值点进行融合计算得到目标像素值，处于平缓区域根据对应方向的源像素插值计算得到目标像素值，插值计算采用的是牛顿二阶插值算法，需要插值点邻域内相关性最大的三个源像素点参与计算。2.根据权利要求1所述的一种基于牛顿二阶插值的图像缩放方法，其特征在于，在步骤a中，坐标化的方式为，将源图像中每个像素点水平和垂直方向上的步...

【专利技术属性】
技术研发人员：王康，杨瑞琪，杨依忠，解光军，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：安徽,34