一种保持边缘特征的数字图像放大方法技术

本发明专利技术具体公开了一种保持边缘特征的数字图像放大方法，包括以下步骤：步骤1确定每个待放大像素的待放大单元以及放大单元；步骤2根据像素与其待放大单元中像素的对应位置，确定其曲面拟合的区域，并重新采样放大单元最外围一圈的像素；步骤3利用空间自适应方法确定放大单元里面中间像素的系数；步骤4利用带约束的最小二乘方法确定放大单元的中间像素；步骤5在放大单元中提取高分辨率像素。本发明专利技术的有益效果：该方法不仅能够保持图像的边缘特征，而且操作比较方便。

【技术实现步骤摘要】

本专利技术属于数字图像处理的图像放缩领域，具体涉及。
技术介绍
空间分辨率是衡量数字图像最重要的参照标准之一。尽管由数字摄像机或扫描仪等获得的原始图像分辨率稳步增长，总有新的技术需要更高的空间分辨率。比如为了计算机医学辅助诊断中观察图像中一个小的病理变化区域需要将病理区域放大；为了适应在高清电视或者显示器上的播放需要标清视频帧的向上转换。除此之外，图像放大在一些领域如计算机视觉、监视、遥感等领域起着有益甚至不可替代的作用。图像放缩最基本的方法是像素复制，这种方法的插值精度是常数，因此放大后的图像具有比较明显的锯齿。双线性插值和双三次插值分别使用一次和三次多项式进行插值，具有较高的精度，且放大效果比像素复制好得多。此外，其它的一些内插方法最小二乘、 分型、和B样条改进了传统的线性和三次插值方法，但是效果并不明显。总之，以图像数据点作为插值数据构造插值曲面的方法，图像数据本质上是区域的采样，区域上的采样是对原场景的平滑，因此，直接应用图像数据构造曲面将丢失一些图像细节。针对传统插值方法的问题，空间自适应的插值算法被提出来。这些算法根据图像局部结构自适应调整插值系数，以更好地保持边缘特征。根据边缘信息的获取方式，自适应算法分为显示和隐式两种。显式方法的基本思想是，将图像的边缘信息显式地表达出来，并引导插值过程沿着边缘方向调整插值系数，但该方法容易受到噪声、模糊等因素的影响。隐式方法则可以较好地解决显式方法的不足，其基本思想是，通过对局部图像邻域内统计量进行估计，从而根据隐式包含有边缘特征信息的统计量对插值过程进行自适应调整。
技术实现思路
针对现有技术存在的缺点，本专利技术所提出了。本专利技术采用的技术方案如下，包括以下步骤步骤I确定各待放大像素以及以该待放大像素为中心的相应的待放大单元以及放大单元；步骤2确定放大单元最外围像素中每个像素的曲面拟合的区域，并重新采样放大单元最外围一圈的像素；步骤3利用空间自适应方法确定放大单元里面中间像素的系数；步骤4利用带约束的最小二乘方法确定出放大单元的中间像素；步骤5在放大单元中提取高分辨率像素。所述的步骤I每个待放大像素的待放大单元以及放大单元的确定方法如下定义4每一个待放大的像素及其周围的3X3邻域为待放大单元，将待放大单元放大两倍之后的 6X6像素单元表示为放大单元。所述的步骤2曲面拟合并重新采样放大单元最外围一圈的像素包括如下步骤(I)根据放大单元像素与其待放大单元中像素的对应位置，确定放大单元最外围一圈像素中每个像素的曲面拟合的区域首先待放大单元里面的每一个像素对应放大单元的里面的每一个2*2像素；其次利用上述的对应关系，确定最外围一圈像素在待放大单元中的对应像素，此对应像素的3*3邻域单元即为曲面拟合的区域。(2)利用二次多项式对拟合的区域进行曲面拟合；首先，步骤(I)中的曲面拟合的区域看成是原场景采样得到，并假设该场景能用分片二次多项式表示；其次，利用式(I)对该拟合区域进行曲面拟合j (X, y) =a1u2+a2uv+a3v2+a4u+a5v+a6u, V e [-1. 5，I. 5] (I)其中u= x-i, V = y-j ；x, y, u, v分别表示曲面参数变量；A,」(x，y)表示对以Pi^1为中心的3*3像素进行拟合的曲面函数；系数a” a2、a3、a4、a5、a6通过带约束的最小二乘求解。(3)利用两倍采样频率的重新采样方法获得最外围一圈的像素；根据步骤(I)中放大单元与待放大单元中像素的对应关系和曲面函数fi, j (X，y)，设定放大单元的最外围一圈的每一个像素的积分范围为待放大单元中对应像素的积分范围的四分之一，即扩大两倍的采样频率进行重采样。所述的步骤3中中间像素系数的确定步骤如下(I)将自回归模型运用于待放大单元的像素待放大单元的每一个像素分别与其相邻的水平_垂直方向及对角线方向的四个像素满足线性一阶关系；(2)通过带约束的最小二乘法把待放大单元像素的系数求解出来首先，将步骤(I)中的一阶线性关系设为约束条件；其次，用带约束的最小二乘法求解出待放大单元中像素的系数；(3)通过待放大单元与放大单元的对应关系求解放大单元中间像素的系数。所述的步骤4中的中间像素的求解过程如下首先，放大单元的中间像素与相邻像素之间满足线性关系，即中间像素分别表示成对角方向上像素以及水平_垂直方向上像素的一阶线性形式；其次，构造一个包含上述线性关系的函数，并根据函数对所有的中间像素求微分来构建方程组；最后，求解方程组得到中间像素。本专利技术的有益效果该方法不仅能够保持图像的边缘特征，而且操作比较方便，另外该方法在曲面拟合求解的过程中，增加了像素与邻域像素之间的相关性，使得图像的边缘结构能够保持，放大图像的边缘特征保留下来；而且在曲面拟合的过程中，我们增加插值条件，使得放大像素的求解更加准确，放大图像更加清晰，而不模糊。实验证明，本方法优于已有的曲面拟合放大方法，而且计算量小，操作比较方便。附图说明图I本专利技术具体操作流程；图2(a)待放大像素对应的待放大单元；图2 (b)放大单元；图3放大单元中的最外围一圈像素对应的拟合区域；图4(a)放大单元里中间像素对角线方向的系数示意图；图4(b)放大单元里中间像素水平_垂直方向的系数示意图；图5(a)待放大单元中像素对角线方向的系数示意图；图5(b)待放大单元中像素水平_垂直方向的系数示意图；图6放大单元与待放大单元的系数对应关系示意图；图7(a)放大单元里中间像素对角线方向的相关性示意图；图7(b)放大单元里中间像素水平_垂直方向的相关性示意图；图8 (a)待放大图像a;图8(b)待放大图像b;图9对应于图8 (a)的高分辨率图像；图10对应于图8 (b)的高分辨率图像；具体实施例方式下面结合附图和实施例作进一步说明如图1所示，基于边缘保持的图像放大，它的步骤描述为步骤1确定各待放大像素以及以该待放大像素为中心的相应的待放大单元以及 放大单元；步骤2确定放大单元最外围像素中每个像素的曲面拟合的区域，并重新采样放大 单元最外围一圈的像素；步骤3利用空间自适应方法确定放大单元里面中间像素的系数；步骤4利用带约束的最小二乘方法确定出放大单元的中间像素；步骤5在放大单元中提取高分辨率像素。所述的步骤1中每个待放大像素的待放大单元以及放大单元的确定设定待放大像素为Pi，」，(i表示横坐标，j表示纵坐标)则以Pi，」为中心的3X3邻 域称为对应Pu的待放大单元，如图2 (a)，图2 (b)中的6X6单元为对应的放大单元，标 记为X”其中i的取值范围0彡i彡35的整数。所述的步骤2中曲面拟合并重采样放大单元最外围一圈的像素如图3所示，以图2 (b)中的Xc V1； X6为例说明最外围像素的求解过程均为最外围像素，a.确定最外围像素的拟合曲面X。，X1； X6三个像素对应待放大单元的Ph, 点。对原图像中以Ph, 为中心的 3X3像素，利用一种已有的二次多项式拟合方法进行曲面拟合，方法如下fi,j(x y) =a1u2+a2uv+a3v2+a4u+a5v+a6u, v G [-1. 5，1. 5](2)其中u= x-i, v = y-j ；x, y, u, v分别表示曲面参数变量；y)表示对以PiH1为中心的3X3像素进行拟合的曲本文档来自技高网...

【技术保护点】
一种保持边缘特征的数字图像放大方法，其特征在于，包括以下步骤：步骤1确定各待放大像素以及以该待放大像素为中心的相应的待放大单元以及放大单元；步骤2确定放大单元最外围像素中每个像素的曲面拟合的区域，并重新采样放大单元最外围一圈的像素；步骤3利用空间自适应方法确定放大单元里面中间像素的系数；步骤4利用带约束的最小二乘方法确定出放大单元的中间像素；步骤5在放大单元中提取高分辨率像素。

【技术特征摘要】
1.一种保持边缘特征的数字图像放大方法，其特征在于，包括以下步骤 步骤I确定各待放大像素以及以该待放大像素为中心的相应的待放大单元以及放大单元； 步骤2确定放大单元最外围像素中每个像素的曲面拟合的区域，并重新采样放大单元最外围一圈的像素； 步骤3利用空间自适应方法确定放大单元里面中间像素的系数； 步骤4利用带约束的最小二乘方法确定出放大单元的中间像素； 步骤5在放大单元中提取高分辨率像素。2.如权利要求I所述的图像放大方法，其特征在于所述的步骤I每个待放大像素的待放大单元以及放大单元的确定方法如下定义每一个待放大的像素及其周围的3*3邻域为待放大单元，将待放大单元放大两倍之后的6*6像素单元表示为放大单元。3.如权利要求I所述的图像放大方法，其特征在于所述的步骤2曲面拟合并重新采样放大单元最外围一圈的像素包括如下步骤 (1)根据放大单元像素与其待放大单元中像素的对应位置，确定放大单元最外围一圈像素中每个像素的曲面拟合的区域首先待放大单元里面的每一个像素对应放大单元的里面的每一个2*2像素；其次利用上述的对应关系，确定最外围一圈像素在待放大单元中的对应像素，此对应像素的3*3邻域单元即为曲面拟合的区域； (2)利用二次多项式对拟合的区域进行曲面拟合； 首先，步骤(I)中的曲面拟合的区域看成是原场景采样得到，并假设该场景能用分片二次多项式表示；其次，利用式(I)对该拟合区域进行曲面拟合f i, j (X，y) =a1u2+a2uv+a3v2+a4u+a5v+a6...

【专利技术属性】
技术研发人员：周元峰，郭莹莹，张彩明，
申请(专利权)人：山东大学，
类型：发明
国别省市：