CN103093418A - 一种改进的数字图像缩放方法 - Google Patents

一种改进的数字图像缩放方法 Download PDF

Info

Publication number
CN103093418A
CN103093418A CN2013100556107A CN201310055610A CN103093418A CN 103093418 A CN103093418 A CN 103093418A CN 2013100556107 A CN2013100556107 A CN 2013100556107A CN 201310055610 A CN201310055610 A CN 201310055610A CN 103093418 A CN103093418 A CN 103093418A
Authority
CN
China
Prior art keywords
sad
dir
diagonal
confidence
smooth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013100556107A
Other languages
English (en)
Other versions
CN103093418B (zh
Inventor
欧阳平
黄海东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Risheng Polytron Technologies Inc
Original Assignee
SHENZHEN GENIUSUN TECHNOLOGY CO LTD
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by SHENZHEN GENIUSUN TECHNOLOGY CO LTD filed Critical SHENZHEN GENIUSUN TECHNOLOGY CO LTD
Priority to CN201310055610.7A priority Critical patent/CN103093418B/zh
Publication of CN103093418A publication Critical patent/CN103093418A/zh
Priority to PCT/CN2013/089614 priority patent/WO2014127665A1/zh
Application granted granted Critical
Publication of CN103093418B publication Critical patent/CN103093418B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

一种改进的数字图像缩放方法,包括根据对目标图像方向和强度检测的结果确定一个插值滤波器对原始图像进行处理得到目标像素值,其特征在于所述插值滤波器由相对平滑的插值滤波器和相对锐利的插值滤波器分别根据其对应的权重修正后相加复合而成,所述方向和强度检测包括分别计算水平、垂直、对角线、反对角线四个方向的SAD值,对水平缩放和垂直缩放分别根据四个方向的SAD值计算水平边缘置信度和垂直边缘置信度,对水平缩放和垂直缩放分别按水平边缘置信度越高以及垂直边缘置信度越高,权重越大的原则设定相对平滑的插值滤波器的权重,本发明解决了插值法缩放长期以来无法克服的边缘锯齿和边缘模糊问题。

Description

一种改进的数字图像缩放方法
技术领域
本发明涉及一种改进的数字图像缩放方法,具体涉及到数字图像的缩放技术,其属于数字图像处理领域类。
背景技术
图像缩放(image scaling)是数字图像处理的基本技术之一。它是指对数字图像的宽高进行调整的过程。图像缩放可以通过两次缩放即水平缩放和垂直缩放来实现。水平缩放改变图像水平像素数量,垂直缩放改变图像垂直像素的数量。如图像从640(宽)×480(高)(宽高比为4:3)到1920×1080(宽高比为16:9),可以先进行垂直缩放将图像放大到640×1080,然后通过水平缩放将图像像素从640×1080调整到1920×1080。这两次缩放的算法原理是相同的。在图像缩放中有一个重要的问题就是缩放过程中图像边缘处理的问题,所述图像边缘指的是图像中物体与图像中其它物体或背景间的分界线,如果这分界线平行于水平轴线我们定义为水平边缘,如果这分界线垂直于水平轴线我们定义为垂直边缘。
在现有技术中,数字图像缩放方法有最近邻方法、双线性插值方法、双三次插值等。最近邻插值比较简单,容易实现,但该方法会在结果图像中产生有明显的边缘锯齿。双线性插值法具有平滑功能, 能有效地克服最近邻法的不足,但会使图像细节和边缘变模糊。双三次插值能提供了较好的折衷,但没有完全有效地解决边缘锯齿和边缘模糊问题,这是插值缩放最大的问题。为了最终解决数字图像缩放过程中产生的边缘锯齿和边缘模糊,亦有学者提出过基于边缘检测的缩放方法,其要点是先检测图像的亚像素边缘,然后对边缘像素和非边缘像素分别进行插值。对于边缘像素,进行基于曲线拟合的插值。而对非边缘像素,进行双线性插值。此方法的缩放效果依赖于亚像素边缘检测结果的好坏。而亚像素边缘检测方法(如Canny边缘检测方法)不仅非常复杂,而且往往难以得到稳定的检测结果。依以上所述,显然现有技术有进一步改进的必要。
发明内容
本发明的目的在于提供一种改进的数字图像缩放方法,以克服现有技术存在的问题。
本发明的一种改进的数字图像缩放方法,包括对目标图像对应的原始图像进行方向和强度检测,然后根据方向和强度检测的结果确定一个插值滤波器对原始图像进行处理得到目标像素值,其特征在于所述插值滤波器为复合插值滤波器,其由相对平滑的插值滤波器和相对锐利的插值滤波器分别根据其对应的权重修正后相加复合而成,所述方向和强度检测包括以与目标像素点对应的原始图像像素点为中心取n×n个像素点,分别计算水平、垂直、对角线、反对角线四个方向的SAD值(方向绝对差之和),对水平缩放根据四个方向的SAD值计算水平边缘置信度,对于垂直缩放根据四个方向的SAD值计算垂直边缘置信度,对水平缩放按水平边缘置信度越高,权重越大的原则基于水平边缘置信度设定相对平滑的插值滤波器的权重,对于垂直缩放,按垂直边缘置信度越高,权重越大的原则基于垂直边缘置信度设定相对平滑的插值滤波器的权重,相对平滑的插值滤波器定义为filter_smooth,相对锐利的插值滤波器定义为filter_sharp,复合滤波器定义为filter_intel,相对平滑的插值滤波器权重定义为w_smooth,则:
filter_intel=(w_smooth)×filter_smooth+(1- w_smooth)×filter_sharp
其中0≤w_smooth≤1。
作为优选,所述水平、垂直、对角线、反对角线四个方向的SAD值计算方向SAD的计算,水平方向的SAD定义为SAD_H,垂直方向的SAD定义为SAD_V,对角线方向的SAD定义为SAD_D,反对角线方向的SAD定义为SAD_A,所述SAD_H为n×n个像素点中行相邻像素点像素值绝对差之和,SAD_V为n×n个像素点中列相邻像素点像素值绝对差之和,SAD_D为n×n个像素点中对角线相邻像素点像素值绝对差之和,SAD_A为n×n个像素点中反对角线相邻像素点像素值绝对差之和。
作为前述优选方案的进一步优先之一,在原始图像是RGB图像的情况下,取RGB图像的Y分量(灰度值)分别计算所述水平、垂直、对角线、反对角线四个方向的SAD值;
作为前述优选方案的进一步优选之二,在原始图像是RGB图像的情况下,分别计算RGB三个通道的水平、垂直、对角线、反对角线四个方向的SAD值,分别取它们中的最大值作为当前目标像素点的水平方向、垂直方向、对角线方向以及反对角线方向的SAD值。
作为进一步优选,对于水平缩放,所述水平边缘置信度为w_dir;
所述w_dir = min(1,max(0,k_p×w_dir_p -k_n×w_dir_n));
其中:
w_dir为水平边缘置信度;
w_dir_p为水平边缘正置信度;
w_dir_n为水平边缘负置信度;
k_p>0为正因素调整系数;
k_n>0为负因素调整系数;
其中:
水平边缘正置信度w_dir_p为max(0,SAD_V-SAD_H);
水平边缘负置信度w_dir_n为w_dir_n0+w_dir_n1+w_dir_n2+w_dir_n3;
其中:
w_dir_n0 = min_SAD, 
w_dir_n1 = avg_SAD,
w_dir_n2 = max(0,SAD_H-min_SAD),
w_dir_n3 = min( max(0,SAD_A-SAD_V), max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值;
本优选方案,给出了通过水平、垂直、对角线、反对角线四个方向的SAD值计算水平边缘置信度的一种方法,对于任何一幅图像来说,一般会存在平行于水平方向的水平边缘以及垂直于水平方向的垂直边缘,在通过水平、垂直、对角线、反对角线四个方向的SAD值计算水平边缘置信度的函数关系式中,我们限定了当目标像素点对应原始图像中的水平边缘或者近水平方向的边缘时,水平边缘置信度较大,如果对应原始图像中的垂直边缘或者近垂直方向边缘时,水平边缘置信度较小;
作为优选,对于水平缩放,所述相对平滑的插值滤波器权重w_smooth依据下式得到:
w_smooth = min(1,k×( w_dir ^alpha)),
其中,k>=1,alpha>=1均为变换调整参数,所述w_dir ^alpha我们称之为水平变换函数,其意义在于当水平边缘置信度较小的时候,此时不大可能是水平边缘,可以将平滑滤波器权重进一步降低,使复合滤波器filter_intel更倾向于相对锐利滤波器;当水平边缘置信度较大的时候,此时是比较可靠的水平边缘,将相对平滑滤波器权重进一步提高,使复合滤波器filter_intel更倾向于相对平滑滤波器。
作为进一步优选,对于垂直缩放,定义所述垂直边缘置信度为v_dir;
所述v_dir = min(1,max(0,k_p×v_dir_p - k_n×v_dir_n));
其中:
v_dir为垂直边缘置信度;
v_dir_p为垂直边缘正置信度;
v_dir_n为垂直边缘负置信度; 
k_p>0为正因素调整系数;
k_n>0为负因素调整系数;
其中:
垂直边缘正置信度v_dir_p为max(0,SAD_H-SAD_V);
垂直边缘负置信度v_dir_n为v_dir_n0+v_dir_n1+v_dir_n2+v_dir_n3;
其中:
v_dir_n0 = min_SAD, 
v_dir_n1 = avg_SAD,
v_dir_n2 = max(0,SAD_V-min_SAD),
v_dir_n3 = min( max(0,SAD_A-SAD_V, max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值;
本优选方案,给出了通过水平、垂直、对角线、反对角线四个方向的SAD值计算垂直边缘置信度的一种方法,亦如前所述那样对于任何一幅图像来说,一般会存在平行于水平方向的水平边缘以及垂直于水平方向的垂直边缘,在通过水平、垂直、对角线、反对角线四个方向的SAD值计算垂直边缘置信度的函数关系式中,我们限定了当目标像素点对应原始图像中的垂直边缘或者近垂直方向的边缘时,垂直边缘置信度较大,如果目标像素点对应原始图像中的水平边缘或者近水平方向边缘时,垂直边缘置信度较小;
作为优选,对于垂直缩放,所述相对平滑的插值滤波器权重w_smooth依据下式得到:
w_smooth = min(1,k×( v_dir ^alpha)),
其中,k>=1,alpha>=1均为变换调整参数,所述“v_dir ^alpha”我们称之垂直变换函数,其意义在于当垂直边缘置信度较小的时候,此时不大可能是垂直边缘,可以将相对平滑的插值滤波器权重进一步降低,使复合滤波器filter_intel更倾向于相对锐利滤波器;当垂直边缘置信度较大的时候,此时是比较可靠的垂直边缘,将相对平滑的插值滤波器权重进一步放大,使复合滤波器filter_intel更倾向于相对平滑滤波器。
本发明的方法,本质上是水平缩放时对于水平边缘最终采用的复合插值滤波器相对平滑,垂直缩放时则对于垂直边缘最终采用的复合插值滤波器相对平滑,同时本发明的复合插值滤波器是由相对平滑的插值滤波器和相对锐利的插值滤波器分别根据其对应的权重修正后相加复合而成的,而相对平滑的插值滤波器和相对锐利的插值滤波器的权重则是以目标像素点对应的原始图像像素点为中心,选定n×n个像素点,并通过计算其水平、垂直、对角线、反对角线四个方向的SAD值,最终以SAD值为基础,在水平缩放中计算出相应的水平边缘置信度,在垂直缩放中计算出相应的垂直边缘置信度而最终确定的,本发明最终实现了在水平缩放时,对水平边缘或近水平方向的边缘的目标像素点采用相对平滑的复合插值滤波器来进行插值,对垂直边缘或近垂直方向的边缘的目标像素则采用相对锐利的复合插值滤波器来进行插值。对于垂直缩放,对垂直边缘或近垂直方向的边缘的目标像素采取相对平滑的复合插值滤波器来进行插值,对水平边缘或近水平方向的边缘的目标像素采用相对锐利的复合插值滤波器来进行插值。显然对于水平缩放时的水平边缘以及垂直缩放时的垂直边缘采用相对平滑的插值滤波器插值,不会带来锯齿问题,而于水平缩放时的垂直边缘以及垂直缩放时的水平边缘采用相对锐利的插值滤波器插值,可确保边缘的清晰度,本发明在水平、垂直、对角线、反对角线四个方向的SAD值的计算的基础上,只附加了水平强边缘和垂直强边缘的检测,就较好地解决了插值法缩放长期以来无法克服的边缘锯齿和边缘模糊问题,另外检测比较简单而且结果可靠鲁棒,同时对系统资源的需求,即对计算机系统的处理能力要求亦低,综上所述,显然本发明的目的得以实现。
附图说明
图1是拟缩放原始图像示意图;
图2是相对平滑的插值滤波器filter_smooth示意图,图中水平轴为像素位置,与目标像素点对应的原始图像像素点位置为0,垂直轴为滤波器取值;
图3是相对锐利的插值滤波器filter_sharp示意图,图中水平轴为像素位置,与目标像素点对应的原始图像像素点位置为0,垂直轴为滤波器取值;
图4是以原始图像像素点为中心形成的3×3矩阵各像素点邻域位置关系图;
图5是以原始图像像素点为中心形成的3×3矩阵水平方向各像素点邻域位置关系图;
图6是以原始图像像素点为中心形成的3×3矩阵垂直方向各像素点邻域位置关系图;
图7是以原始图像像素点为中心形成的3×3矩阵对角线方向各像素点邻域位置关系图;
图8是以原始图像像素点为中心形成的3×3矩阵反对角线方向各像素点邻域位置关系图;
图9是相对平滑的插值滤波器变换函数图图中水平轴为边缘置信度,垂直轴为平滑滤波器权重;
各图中的主要标识如下:
1为图像中物体;
101为垂直边缘;
102为水平边缘;
103为近水平方向边缘;
104为近垂直方向边缘;
2为水平轴线;
3为垂直轴线;
4为水平边缘置信度曲线;
5为经相对平滑的插值滤波器变换函数变换得来的相对平滑的插值滤波器权重曲线;
6为相对平滑的插值滤波器曲线;
7为相对锐利的插值滤波器曲线。
具体实施方式
以下将结合本发明较佳实施例所提供的一种改进的数字图像缩放方法及其附图对本发明作进一步说明。
附图1给出了拟进行缩放处理的原始图像示意图,如附图1所示,图像中包括物体1,包括明显与水平轴线2垂直的垂直缘101,以及与水平轴线2平行的水平缘102, 当然垂直边缘101与垂直轴线3平行,水平边缘102与垂直轴线3垂直,同时,所述物体1还可能存在近103,亦称水平强边缘和近垂直方向边缘10[tao1] 4,亦称垂直强边缘,本较佳实施例给出的这种改进的数字图像缩放方法,首先确定一个相对平滑的插值滤波器filter_smooth和一个相对锐利的插值滤波器filter_sharp,在本较佳实施例中,所述相对平滑的插值滤波器filter_smooth如附图2所示,所述相对锐利的插值滤波器filter_sharp如附图3所示,同时定义复合滤波器filter_intel等于(w_smooth)×filter_smooth+(1- w_smooth)×filter_sharp,其中w_smooth为相对平滑的插值滤波器权重,且0≤w_smooth≤1;
具体缩放时,先以与目标像素点对应的原始图像像素点22为中心取3×3个像素点,附图4同时标示出了以原始图像像素点22为中心形成的3×3矩阵各像素点邻域位置关系,如附图4所示,3×3矩阵中第一行从左到右各像素点分别标示为:11、12、13,第二行分别标示为21、22、23,第三行分别标示为:31、32、33;如附图5、附图6、附图7、附图8所示,所述3×3个像素点中,水平方向、垂直方向、对角线方向以及反对角线方向均为三行,定义P为像素点的像素值,则我们有SAD_H为∣p12-p11∣、∣p13-p12∣、∣p22-p21∣、∣p23-p22∣、∣p32-p31∣、∣p33-p32∣六个邻域像素点像素值绝对差之和,SAD_V为∣p21-p11∣、∣p31-p21∣、∣p22-p12∣、∣p32-p22∣、∣p23-p13∣、∣p33-p13∣六个邻域像素点像素点像素值绝对差之和,SAD_D为∣p32-p21∣、∣p22-p11∣、∣p33-p22∣、∣p23-p12∣四个邻域像素点像素点像素值绝对差之和,SAD_A为∣p21-p12∣、∣p22-p13∣、∣p31-p22∣、∣p32-p23∣四个邻域像素点像素值绝对差之和,在本较佳实施例中,所述原始图像1为RGB图像,为此取RGB图像的Y分量,即灰度值作为计算SAD的依据,
在本较佳实施例中,如果[tao2] ,可对水平边缘置信度进行计算,在本较佳实施例中定义水平边缘正置信度为w_dir_p,然后根据w_dir_p=max(0,SAD_V-SAD_H)求得水平边缘正置信度w_dir_p的值;再定义水平边缘负置信度为w_dir_n,然后根据w_dir_n=w_dir_n0+w_dir_n1+w_dir_n2+w_dir_n3取得水平边缘负置信度w_dir_n的值,其中:
w_dir_n0 = min_SAD, 
w_dir_n1 = avg_SAD,
w_dir_n2 = max(0,SAD_H-min_SAD),
w_dir_n3 = min( max(0,SAD_A-SAD_V), max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值;
最后根据水平边缘置信度为w_dir = min(1,max(0,k_p×w_dir_p -k_n×w_dir_n))取得水平边缘置信度为w_dir的值;
其中:
k_p>0为正因素调整系数,其取值范围为0-4,在本较佳实施例中取值2;
k_n>0为负因素调整系数,其取值范围为0-4,在本较佳实施例中取值2;
过大的k_p值和太小k_n值会有模糊图像细节的风险,反之会有导致水平边缘锯齿的风险。
 而所述相对平滑的插值滤波器权重w_smooth依据下式给出的相对平滑的插值滤波器变换函数得到:
w_smooth = min(1,k×( w_dir ^alpha)),
其中,k>=1,alpha>=1均为变换调整参数,在本较佳实施例中k值为1.5,alpha值为2;w_dir ^alpha为水平变换函数,附图9给出了在k值为1.5,alpha值为2的情况下,经水平变换函数处理后的相对平滑的插值滤波器权重w_smooth的变化情况,图中虚线3为经水平变换函数处理前的相对平滑的插值滤波器权重w_smooth,实线3为经水平变换函数处理后的相对平滑的插值滤波器权重w_smooth。
依前所述,相对平滑的插值滤波器定义为filter_smooth,相对锐利的插值滤波器定义为filter_sharp,复合插值滤波器定义为filter_intel,相对平滑的插值滤波器权重定义为w_smooth,则我们有:
filter_intel=(w_smooth)×filter_smooth+(1- w_smooth)×filter_sharp;
最后按我们得到的复合插值滤波器filter_intel进行水平缩放的插值处理得到目标图像;
如果进行垂直缩放,在本较佳实施例中定义垂直边缘正置信度为v_dir_p,然后根据v_dir_p=max(0,SAD_H-SAD_V)求得垂直边缘正置信度v_dir_p的值;再定义垂直边缘负置信度为v_dir_n,然后根据v_dir_n=v_dir_n0+v_dir_n1+v_dir_n2+v_dir_n3取得垂直边缘负置信度v_dir_n的值;
其中:
v_dir_n0 = min_SAD, 
v_dir_n1 = avg_SAD,
v_dir_n2 = max(0,SAD_V-min_SAD),
v_dir_n3 = min( max(0,SAD_A-SAD_V, max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值;
然后根据垂直边缘置信度v_dir = min(1,max(0,k_p×v_dir_p – k_n×v_dir_n))取得垂直边缘置信度v_dir的值,其中:
k_p>0为正因素调整系数,其取值范围为0-4,在本较佳实施例中取值2; 
k_n>0为负因素调整系数,其取值范围为0-4,在本较佳实施例中取值2; 
过大的k_p值和太小k_n值会有模糊图像细节的风险,反之会有导致垂直边缘锯齿的风险。
相对平滑的插值滤波器权重w_smooth依据下式给出的相对平滑的插值滤波器变换函数得到:
w_smooth = min(1,k×( v_dir ^alpha)),
其中,k>=1,alpha>=1均为变换调整参数,所述“v_dir ^alpha”我们称之垂直变换函数,意义同水平缩放,在此不再详细说明。
依前所述,相对平滑的插值滤波器定义为filter_smooth,相对锐利的插值滤波器定义为filter_sharp,复合插值滤波器定义为filter_intel,相对平滑的插值滤波器权重定义为w_smooth,则我们有:
filter_intel=(w_smooth)×filter_smooth+(1- w_smooth)×filter_sharp;
最后按我们得到的复合插值滤波器filter_intel进行垂直缩放的插值处理得到目标图像。
综上所述,本发明的一种改进的数字图像缩放方法,包括根据对目标图像方向和强度检测的结果确定一个插值滤波器对原始图像进行处理得到目标像素值,其特征在于所述插值滤波器由相对平滑的插值滤波器和相对锐利的插值滤波器分别根据其对应的权重修正后相加复合而成,所述方向和强度检测包括分别计算水平、垂直、对角线、反对角线四个方向的SAD值,对水平缩放和垂直缩放分别根据四个方向的SAD值计算水平边缘置信度和垂直边缘置信度,对水平缩放和垂直缩放分别按水平边缘置信度越高以及垂直边缘置信度越高,权重越大的原则设定相对平滑的插值滤波器的权重,本发明解决了插值法缩放长期以来无法克服的边缘锯齿和边缘模糊问题。
以上所述,仅为本发明的具体实施方式。但本发明的保护范围并不局限于此。任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内,因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种改进的数字图像缩放方法,包括对目标图像对应的原始图像进行方向和强度检测,然后根据方向和强度检测的结果确定一个插值滤波器对原始图像进行处理得到目标像素值,其特征在于所述插值滤波器为复合插值滤波器,其由相对平滑的插值滤波器和相对锐利的插值滤波器分别根据其对应的权重修正后相加复合而成,所述方向和强度检测包括以与目标像素点对应的原始图像像素点为中心取n×n个像素点,分别计算水平、垂直、对角线、反对角线四个方向的SAD值,对水平缩放根据四个方向的SAD值计算水平边缘置信度,对于垂直缩放根据四个方向的SAD值计算垂直边缘置信度,对水平缩放按水平边缘置信度越高,权重越大的原则基于水平边缘置信度设定相对平滑的插值滤波器的权重,对于垂直缩放,按垂直边缘置信度越高,权重越大的原则基于垂直边缘置信度设定相对平滑的插值滤波器的权重。
2.根据权利要求1所述的一种改进的数字图像缩放方法,其特征在于相对平滑的插值滤波器定义为filter_smooth,相对锐利的插值滤波器定义为filter_sharp,复合滤波器定义为filter_intel,相对平滑的插值滤波器权重定义为w_smooth,则:
filter_intel=(w_smooth)×filter_smooth+(1- w_smooth)×filter_sharp;
其中0≤w_smooth≤1。
3.根据权利要求2所述的一种改进的数字图像缩放方法,其特征在于所述水平、垂直、对角线、反对角线四个方向的SAD值计算方向SAD的计算,水平方向的SAD为n×n个像素点中行相邻像素点像素值绝对差之和,垂直方向的SAD为n×n个像素点中列相邻像素点像素值绝对差之和,对角线方向的SAD为n×n个像素点中对角线相邻像素点像素值绝对差之和,反对角线方向的SAD为n×n个像素点中反对角线相邻像素点像素值绝对差之和。
4.根据权利要求3所述的一种改进的数字图像缩放方法,其特征在于在原始图像是RGB图像的情况下,取RGB图像的Y分量(灰度值)分别计算所述水平、垂直、对角线、反对角线四个方向的SAD值。
5.根据权利要求3所述的一种改进的数字图像缩放方法,其特征在于在原始图像是RGB图像的情况下,分别计算RGB三个通道的水平、垂直、对角线、反对角线四个方向的SAD值,分别取它们中的最大值作为当前目标像素点的水平方向、垂直方向、对角线方向以及反对角线方向的SAD值。
6.根据权利要求3所述的一种改进的数字图像缩放方法,其特征在于对于水平缩放,所述水平边缘置信度为w_dir;
所述w_dir = min(1,max(0,k_p×w_dir_p – k_n×w_dir_n));
其中:
SAD_H为水平方向的SAD;
SAD_V为垂直方向的SAD;S
AD_D为对角线方向的SAD;
SAD_A为反对角线方向的SAD;
w_dir为水平边缘置信度;
w_dir_p为水平边缘正置信度;
w_dir_n为水平边缘负置信度;
k_p>0为正因素调整系数;
k_n>0为负因素调整系数;
其中:
水平边缘正置信度w_dir_p为max(0,SAD_V-SAD_H);
水平边缘负置信度w_dir_n为w_dir_n0+w_dir_n1+w_dir_n2+w_dir_n3;
其中:
w_dir_n0 = min_SAD, 
w_dir_n1 = avg_SAD,
w_dir_n2 = max(0,SAD_H-min_SAD),
w_dir_n3 = min( max(0,SAD_A-SAD_V), max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值。
7.根据权利要求6所述的一种改进的数字图像缩放方法,其特征在于所述相对平滑的插值滤波器权重w_smooth依据下式得到:
w_smooth = min(1,k×( w_dir ^alpha));
其中,k>=1,alpha>=1均为变换调整参数。
8.根据权利要求3所述的一种改进的数字图像缩放方法,其特征在于对于垂直缩放,所述垂直边缘置信度为v_dir;
所述v_dir = min(1,max(0,k_p×v_dir_p – k_n×v_dir_n));
其中:
SAD_H为水平方向的SAD;
SAD_V为垂直方向的SAD;
S AD_D为对角线方向的SAD;
SAD_A为反对角线方向的SAD;
v_dir为垂直边缘置信度;
v_dir_p为垂直边缘正置信度;
v_dir_n为垂直边缘负置信度; 
k_p>0为正因素调整系数;
k_n>0为负因素调整系数;
其中:
垂直边缘正置信度v_dir_p为max(0,SAD_H-SAD_V);
垂直边缘负置信度v_dir_n为v_dir_n0+v_dir_n1+v_dir_n2+v_dir_n3;
其中:
v_dir_n0 = min_SAD, 
v_dir_n1 = avg_SAD,
v_dir_n2 = max(0,SAD_V-min_SAD),
v_dir_n3 = min( max(0,SAD_A-SAD_V, max(0,SAD_D-SAD_V) )
其中:
所述min_SAD为水平、垂直、对角线、反对角线四个方向的SAD的最小值;
所述avg_SAD为水平、垂直、对角线、反对角线四个方向的SAD的均值。
9.根据权利要求8所述的一种改进的数字图像缩放方法,其特征在于对于垂直缩放,所述相对平滑的插值滤波器权重w_smooth依据下式得到:
w_smooth = min(1,k×( v_dir ^alpha));
其中,k>=1,alpha>=1均为变换调整参数。
CN201310055610.7A 2013-02-21 2013-02-21 一种改进的数字图像缩放方法 Expired - Fee Related CN103093418B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310055610.7A CN103093418B (zh) 2013-02-21 2013-02-21 一种改进的数字图像缩放方法
PCT/CN2013/089614 WO2014127665A1 (zh) 2013-02-21 2013-12-16 一种改进的数字图像缩放方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310055610.7A CN103093418B (zh) 2013-02-21 2013-02-21 一种改进的数字图像缩放方法

Publications (2)

Publication Number Publication Date
CN103093418A true CN103093418A (zh) 2013-05-08
CN103093418B CN103093418B (zh) 2015-08-26

Family

ID=48205954

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310055610.7A Expired - Fee Related CN103093418B (zh) 2013-02-21 2013-02-21 一种改进的数字图像缩放方法

Country Status (2)

Country Link
CN (1) CN103093418B (zh)
WO (1) WO2014127665A1 (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103530845A (zh) * 2013-10-19 2014-01-22 深圳市晶日盛科技有限公司 一种改进的图像缩放方法
WO2014127665A1 (zh) * 2013-02-21 2014-08-28 深圳市晶日盛科技有限公司 一种改进的数字图像缩放方法
CN104794692A (zh) * 2015-04-16 2015-07-22 中国科学院自动化研究所 一种图像去锯齿的系统
US9246990B2 (en) 2014-02-14 2016-01-26 Google Inc. Methods and systems for predicting conversion rates of content publisher and content provider pairs
US9461936B2 (en) 2014-02-14 2016-10-04 Google Inc. Methods and systems for providing an actionable object within a third-party content slot of an information resource of a content publisher
WO2016165097A1 (zh) * 2015-04-16 2016-10-20 中国科学院自动化研究所 一种图像去锯齿的系统
CN112534466A (zh) * 2018-08-02 2021-03-19 苹果公司 定向缩放系统和方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11854157B2 (en) * 2020-03-18 2023-12-26 Intel Corporation Edge-aware upscaling for improved screen content quality

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101188017A (zh) * 2007-12-18 2008-05-28 上海广电集成电路有限公司 数字图像的缩放方法以及系统
CN102568443A (zh) * 2010-12-29 2012-07-11 无锡华润矽科微电子有限公司 一种数字图像缩放方法
CN102596050A (zh) * 2009-10-27 2012-07-18 株式会社日立医疗器械 超声波成像装置、超声波成像方法及超声波成像用程序
CN102779335A (zh) * 2011-05-13 2012-11-14 华晶科技股份有限公司 数字图像处理装置及其处理方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050094899A1 (en) * 2003-10-29 2005-05-05 Changick Kim Adaptive image upscaling method and apparatus
JP2010157877A (ja) * 2008-12-26 2010-07-15 Sharp Corp 画像処理装置、画像処理方法及びコンピュータプログラム
CN103093418B (zh) * 2013-02-21 2015-08-26 深圳市晶日盛科技有限公司 一种改进的数字图像缩放方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101188017A (zh) * 2007-12-18 2008-05-28 上海广电集成电路有限公司 数字图像的缩放方法以及系统
CN102596050A (zh) * 2009-10-27 2012-07-18 株式会社日立医疗器械 超声波成像装置、超声波成像方法及超声波成像用程序
CN102568443A (zh) * 2010-12-29 2012-07-11 无锡华润矽科微电子有限公司 一种数字图像缩放方法
CN102779335A (zh) * 2011-05-13 2012-11-14 华晶科技股份有限公司 数字图像处理装置及其处理方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014127665A1 (zh) * 2013-02-21 2014-08-28 深圳市晶日盛科技有限公司 一种改进的数字图像缩放方法
CN103530845A (zh) * 2013-10-19 2014-01-22 深圳市晶日盛科技有限公司 一种改进的图像缩放方法
US10210140B2 (en) 2014-02-14 2019-02-19 Google Llc Methods and systems for providing an actionable object within a third-party content slot of an information resource of a content publisher
US9246990B2 (en) 2014-02-14 2016-01-26 Google Inc. Methods and systems for predicting conversion rates of content publisher and content provider pairs
US9461936B2 (en) 2014-02-14 2016-10-04 Google Inc. Methods and systems for providing an actionable object within a third-party content slot of an information resource of a content publisher
US10067916B2 (en) 2014-02-14 2018-09-04 Google Llc Methods and systems for providing an actionable object within a third-party content slot of an information resource of a content publisher
CN104794692B (zh) * 2015-04-16 2017-08-25 中国科学院自动化研究所 一种图像去锯齿的系统
WO2016165097A1 (zh) * 2015-04-16 2016-10-20 中国科学院自动化研究所 一种图像去锯齿的系统
CN104794692A (zh) * 2015-04-16 2015-07-22 中国科学院自动化研究所 一种图像去锯齿的系统
US10410326B2 (en) 2015-04-16 2019-09-10 Institute Of Automation Chinese Academy Of Sciences Image anti-aliasing system
CN112534466A (zh) * 2018-08-02 2021-03-19 苹果公司 定向缩放系统和方法
CN112534466B (zh) * 2018-08-02 2022-05-13 苹果公司 定向缩放系统和方法
US11941785B2 (en) 2018-08-02 2024-03-26 Apple Inc. Directional scaling systems and methods

Also Published As

Publication number Publication date
WO2014127665A1 (zh) 2014-08-28
CN103093418B (zh) 2015-08-26

Similar Documents

Publication Publication Date Title
CN103093418B (zh) 一种改进的数字图像缩放方法
CN110276734B (zh) 图像畸变校正方法和装置
CN110264426B (zh) 图像畸变校正方法和装置
CN103500435B (zh) 一种边缘导向的自适应图像插值方法及其vlsi实现装置
US10614551B2 (en) Image interpolation methods and related image interpolation devices thereof
CN109741356B (zh) 一种亚像素边缘检测方法及系统
CN102620673A (zh) 基于图像分析的隧道变形在线监测系统及其应用
US11062464B2 (en) Image processing apparatus, method, and storage medium to derive optical flow
JP2012191465A (ja) 画像処理装置、および画像処理方法、並びにプログラム
WO2010021009A1 (ja) 画像補正装置および画像補正方法
JP6465982B2 (ja) 画像処理装置、画像処理方法及びプログラム
CN104700360B (zh) 基于边缘自适应的图像缩放方法及系统
CN103167218A (zh) 一种基于非局部性的超分辨率重建方法和设备
CN114862725B (zh) 基于光流法实现运动感知模糊特效的方法及装置
CN104200426A (zh) 图像插值方法和装置
KR102404557B1 (ko) 영상 처리 장치 및 영상 처리 방법
US20150227779A1 (en) Image processing device, image processing method, and program
JP2012208553A5 (zh)
JP5012967B2 (ja) 画像処理装置及び方法、並びにプログラム
CN110349090A (zh) 一种基于牛顿二阶插值的图像缩放方法
CN108648205A (zh) 一种亚像素边缘检测方法
US20150110405A1 (en) Point spread function cost function with non-uniform weights
WO2009150696A1 (ja) 画像補正装置および画像補正方法
CN107977981B (zh) 一种运动目标跟踪方法及装置
CN104574276A (zh) 一种基于光流法的图像对齐的方法和装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C56 Change in the name or address of the patentee
CP03 Change of name, title or address

Address after: 518000 Guangdong city of Shenzhen province Baoan District manhole Street Maan Shan Industrial Zone, the first saddle Long Road No. 8 on the third.

Patentee after: Shenzhen Risheng Polytron Technologies Inc

Address before: 518000, Shenzhen, Guangdong province Baoan District Fuyong Avenue deep color building, 8 floor

Patentee before: Shenzhen Geniusun Technology Co.,Ltd.

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150826

Termination date: 20190221

CF01 Termination of patent right due to non-payment of annual fee