CN108616701A - 一种基于Sinc滤波器的图像缩放方法 - Google Patents

一种基于Sinc滤波器的图像缩放方法 Download PDF

Info

Publication number
CN108616701A
CN108616701A CN201611140258.7A CN201611140258A CN108616701A CN 108616701 A CN108616701 A CN 108616701A CN 201611140258 A CN201611140258 A CN 201611140258A CN 108616701 A CN108616701 A CN 108616701A
Authority
CN
China
Prior art keywords
image
horizontal direction
scaling
phase
filter
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.)
Pending
Application number
CN201611140258.7A
Other languages
English (en)
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.)
Xian Aeronautics Computing Technique Research Institute of AVIC
Original Assignee
Xian Aeronautics Computing Technique Research Institute of AVIC
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 Xian Aeronautics Computing Technique Research Institute of AVIC filed Critical Xian Aeronautics Computing Technique Research Institute of AVIC
Priority to CN201611140258.7A priority Critical patent/CN108616701A/zh
Publication of CN108616701A publication Critical patent/CN108616701A/zh
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N5/00Details of television systems
    • H04N5/222Studio circuitry; Studio devices; Studio equipment
    • H04N5/262Studio circuits, e.g. for mixing, switching-over, change of character of image, other special effects ; Cameras specially adapted for the electronic generation of special effects
    • H04N5/2628Alteration of picture size, shape, position or orientation, e.g. zooming, rotation, rolling, perspective, translation
    • 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
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/80Details of filtering operations specially adapted for video compression, e.g. for pixel interpolation

Landscapes

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

Abstract

本发明属于嵌入式计算机图像处理技术领域,特别是涉及一种图像缩放方法。图像缩放一般是采用近邻取样法,双线性插值,双三次插值等方法。本发明目的是克服已有的缩放方法的不足,设计一种使用Sinc滤波器进行图像缩放的方法。本发明的实现方法灵活多样,既可以使用软件实现也可以使用硬件比如DSP和可编程逻辑实现。

Description

一种基于Sinc滤波器的图像缩放方法
技术领域
本发明属于嵌入式计算机图像处理技术领域,特别是涉及一种图像缩放方法。
背景技术
近年来,高分辨率显示器被广泛应用于各个行业。但传统的低分辨率图像采样装置无法进行及时的升级和更换,这就造成低分辨率图像无法在高分辨率显示器上显示的问题。而在某些特殊行业还存在高分辨率图像无法在较低分辨率显示设备上正确显示的问题。这都要用到图像缩放技术。
图像缩放一般是采用近邻取样法,双线性插值,双三次插值等方法。这些方法要么效果较差,要么仅能利用相邻像素信息,要么计算复杂度较高,不适合嵌入式设备实现。所以就需要一种能够利用较多临近像素信息兼顾效果并适合嵌入式设备实现的方法来实现图像缩放功能。
发明内容
本发明目的是克服上述已有技术的不足,设计一种能够兼顾效果的更适应嵌入式设备的图像缩放方法。
本发明的技术方案是利用Sinc滤波器的滤波方法来实现图像的缩放。归一化Sinc函数表达式为为理想低通滤波函数,应用于图像的缩放具有很好的效果,使用此滤波器进行缩放操作时可以使用任意数量的抽头,抽头数量越大滤波器越接近理想滤波器,缩放的效果越好,但消耗的计算资源越多。所以可以根据嵌入式设备的计算资源量来指定抽头数量。图像缩放时,水平和竖直方向可以独立进行,所以可以以相似的结构实现水平和竖直方向的缩放。
本发明包括以下步骤:
一种基于Sinc滤波器的图像缩放方法,其特征是,本方法包括以下步骤:
第一步,指定滤波器的抽头数量,抽头数量大于等于2;
第二步,根据缩放比例确定水平方向和竖直方向滤波相位数量,由此确定水平方向和竖直方向的滤波相位偏移;
第三步,根据相位偏移和抽头数量在归一化Sinc函数上进行采样;
第四步,将水平方向的Sinc函数采样值作为水平方向滤波器的系数,将竖直方向Sinc函数采样值作为竖直方向滤波器的系数;
第五步,使用水平方向滤波器对图像的行数据进行滤波,得到水平方向缩放后的行数据;
第六步,将水平方向缩放后的行数据进行缓存;
第七步,对水平方向缩放后的图像数据进行竖直方向的滤波,得到最终缩放后的图像数据。
本发明具有如下优点:
1、可以利用多于2个相邻像素信息;
2、经典滤波器结构,更适用于嵌入式设备;
3、此方法不仅适用于通用平台的软件实现,还尤其适用于DSP和FPGA实现。
附图说明
图1:本发明流程示意图;
图2:Sinc函数曲线图。
具体实施方式
下面结合附图对本发明做进一步的描述:
以1024x768分辨率图像放大到1280x1024分辨率图像为例:
第一步,指定滤波器的抽头数量为4,即TAP_NUM=4,此处在计算资源足够情况下可以选择更多抽头,为更易于硬件实现和适应图像特点,建议使用如2,4,8,12,16等值,抽头数量越大滤波器越接近理想滤波器,缩放的效果越好,但消耗的计算资源越多;
第二步,水平方向上缩放比例为1024/1280=0.8,所以水平方向的相位数量为5可以满足放大需求,并得到每个滤波相位的偏移为(0:0.8:0.6:0.4:0.2);竖直方向上缩放比例为768/1024=0.75,所以竖直方向上相位数量为4可以满足放大需求,并得到每个滤波相位的偏移为(0:0.75:0.5:0.25);此相位偏移即为相对于原像素位置新生成的新像素位置,这里以行为例,假定新旧图像的第一个像素是重合的,所以偏移为0,行上的第二个新像素生成位置为原像素位置的第0.8,第三个像素为第1.6(所以相对相邻像素偏移为0.6)……以此类推,并在第6个像素时循环;
第三步,在不同相位偏移上对归一化Sinc函数(如图2)进行采样,以行为例,当取PHASE-1时,初始偏移为0,所以将x=(-2:-1:0:1)分别代入表达式得到PHASE-1的抽头系数;当取PHASE-2时,初始偏移为0.8,将x=(-2+0.8:-1+0.8:0+0.8:1+0.8)即x=(-1.2:-0.2:0.8:1.8)分别代入表达式得到PHASE-2的抽头系数……以此类推可以计算出如下两表的采样值:
水平采样 TAP-1 TAP-2 TAP-3 TAP-4
PHASE-1 -0.000 0.000 1.000 0.000
PHASE-2 -0.104 0.234 0.935 -0.156
PHASE-3 -0.189 0.505 0.757 -0.216
PHASE-4 -0.216 0.757 0.505 -0.189
PHASE-5 -0.156 0.935 0.234 -0.104
竖直采样 TAP-1 TAP-2 TAP-3 TAP-4
PHASE-1 -0.000 0.000 1.000 0.000
PHASE-2 -0.129 0.300 0.900 -0.180
PHASE-3 -0.212 0.636 0.636 -0.212
PHASE-4 -0.180 0.900 0.300 -0.129
第四步,将水平方向采样值作为水平滤波器的系数,竖直方向采样值作为竖直滤波器的系数。
第五步,对图像行进行水平方向滤波,由于前述相位偏移去除了所参考的原图像像素的位置信息,这里需要使用之前计算出的缩放比例对每个新像素所参考的原像素位置进行计算,如行方向上缩放比例为1024/1280=0.8,设新像素位置为n(n从0开始),所参考的TAP-1位置为n×0.8向上取整再减2,如第一个新像素(n=0)所参考的TAP-1位置为(0×0.8向上取整再减2)=-2,第二个新像素参考TAP-1位置-1,第三个TAP-1位置参考0……以此类推,可以得到如下计算结果:
第一个像素使用PHASE-1:
xout0=xin-2·(-0.000)+xin-1·(0.000)+xin0·(1.000)+xin1·(0.000);
第二个像素使用PHASE-2:
xout1=xin-1·(-0.104)+xin0·(0.234)+xin1·(0.935)+xin2·(-0.156)
第三个像素使用PHASE-3:
xout2=xin0·(-0.189)+xin1·(0.505)+xin2·(0.757)+xin3·(-0.216)
第四个像素使用PHASE-4:
xout3=xin1·(-0.216)+xin2·(0.757)+xin3·(0.505)+xin4·(-0.189)
第五个像素使用PHASE-5:
xout4=xin2·(-0.156)+xin3·(0.935)+xin4·(0.234)+xin5·(-0.104)
第六个像素使用PHASE-1:
xout5=xin2·(-0.000)+xin3·(0.000)+xin4·(1.000)+xin5·(0.000)
……
依次循环,其中xin-2、xin-1为像素向左侧的扩展,数值等于xin0
第六步,将水平方向缩放后的行数据进行缓存;
第七步,对水平方向缩放后的图像数据进行竖直方向的滤波,原理同第五步,得到最终缩放后的图像数据:
第一个像素使用PHASE-1:
pout(0,0)=x(out0,-2)·(-0.000)+x(out0,-1)·(0.000)+x(out0,0)·(1.000)+x(out0,1)·(0.000);
第二个像素使用PHASE-2:
pout(0,1)=x(out0,-1)·(-0.129)+x(out0,0)·(0.300)+x(out0,1)·(0.900)+x(out0,2)·(-0.180);
第三个像素使用PHASE-3:
pout(0,2)=x(out0,0)·(-0.212)+x(out0,1)·(0.636)+x(out0,2)·(0.636)+x(out0,3)·(-0.212);
第四个像素使用PHASE-4:
pout(0,3)=x(out0,1)·(-0.180)+x(out0,2)·(0.900)+x(out0,3)·(0.300)+x(out0,3)·(-0.129);
第五个像素使用PHASE-1:
pout(0,4)=x(out0,1)·(-0.000)+x(out0,2)·(0.000)+x(out0,3)·(1.000)+x(out0,4)·(0.000);
……
依次循环,其中x(out0,-2)、x(out0,-1)为像素向上侧的扩展,数值等于x(out0,0)

Claims (1)

1.一种基于Sinc滤波器的图像缩放方法,其特征是,本方法包括以下步骤:
第一步,指定滤波器的抽头数量,抽头数量大于等于2;
第二步,根据缩放比例确定水平方向和竖直方向滤波相位数量,由此确定水平方向和竖直方向的滤波相位偏移;
第三步,根据相位偏移和抽头数量在归一化Sinc函数上进行采样;
第四步,将水平方向的Sinc函数采样值作为水平方向滤波器的系数,将竖直方向Sinc函数采样值作为竖直方向滤波器的系数;
第五步,使用水平方向滤波器对图像的行数据进行滤波,得到水平方向缩放后的行数据;
第六步,将水平方向缩放后的行数据进行缓存;
第七步,对水平方向缩放后的图像数据进行竖直方向的滤波,得到最终缩放后的图像数据。
CN201611140258.7A 2016-12-12 2016-12-12 一种基于Sinc滤波器的图像缩放方法 Pending CN108616701A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611140258.7A CN108616701A (zh) 2016-12-12 2016-12-12 一种基于Sinc滤波器的图像缩放方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611140258.7A CN108616701A (zh) 2016-12-12 2016-12-12 一种基于Sinc滤波器的图像缩放方法

Publications (1)

Publication Number Publication Date
CN108616701A true CN108616701A (zh) 2018-10-02

Family

ID=63643266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611140258.7A Pending CN108616701A (zh) 2016-12-12 2016-12-12 一种基于Sinc滤波器的图像缩放方法

Country Status (1)

Country Link
CN (1) CN108616701A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112258396A (zh) * 2020-12-17 2021-01-22 恒银金融科技股份有限公司 一种用于字符图像缩放的方法
CN112862673A (zh) * 2019-11-12 2021-05-28 上海途擎微电子有限公司 自适应图像缩放方法及自适应图像缩放装置、存储装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454003A (zh) * 2003-05-29 2003-11-05 上海交通大学 用二维多相插值滤波器实现视频图象格式转换的方法
CN102194214A (zh) * 2010-03-05 2011-09-21 英特尔公司 用于图片缩放的系统、方法和计算机程序产品
CN103327218A (zh) * 2012-03-01 2013-09-25 苹果公司 图像处理系统和方法
WO2015138311A1 (en) * 2014-03-11 2015-09-17 Qualcomm Incorporated Phase control multi-tap downscale filter
US20160217553A1 (en) * 2015-01-22 2016-07-28 Samsung Display Co., Ltd. Image processing device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454003A (zh) * 2003-05-29 2003-11-05 上海交通大学 用二维多相插值滤波器实现视频图象格式转换的方法
CN102194214A (zh) * 2010-03-05 2011-09-21 英特尔公司 用于图片缩放的系统、方法和计算机程序产品
CN103327218A (zh) * 2012-03-01 2013-09-25 苹果公司 图像处理系统和方法
WO2015138311A1 (en) * 2014-03-11 2015-09-17 Qualcomm Incorporated Phase control multi-tap downscale filter
US20160217553A1 (en) * 2015-01-22 2016-07-28 Samsung Display Co., Ltd. Image processing device

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112862673A (zh) * 2019-11-12 2021-05-28 上海途擎微电子有限公司 自适应图像缩放方法及自适应图像缩放装置、存储装置
CN112258396A (zh) * 2020-12-17 2021-01-22 恒银金融科技股份有限公司 一种用于字符图像缩放的方法
CN112258396B (zh) * 2020-12-17 2021-04-06 恒银金融科技股份有限公司 一种用于字符图像缩放的方法

Similar Documents

Publication Publication Date Title
CN108154474B (zh) 一种超分辨率图像重构方法、装置、介质及设备
TWI714397B (zh) 視頻處理方法、裝置以及電腦儲存介質
US8989509B2 (en) Streaming wavelet transform
JP2016524423A (ja) ユーザ選択関心領域を超解像するためのシステムおよび方法
CN103996170A (zh) 一种具有超高解析度的图像边缘锯齿消除方法
CN110428382B (zh) 一种用于移动终端的高效视频增强方法、装置和存储介质
KR20200067920A (ko) 이미지 프로세싱 파이프라인 내의 이미지 데이터의 다중 레이트 프로세싱
US20130108186A1 (en) Image processing method for boundary resolution enhancement
US11887350B2 (en) Some automated and semi-automated tools for linear feature extraction in two and three dimensions
US11107191B2 (en) Apparatus and method for detail enhancement in super-resolution imaging using mobile electronic device
CN103927169A (zh) 一种实现gdi+绘制电子地图放缩的方法
CN108616701A (zh) 一种基于Sinc滤波器的图像缩放方法
CN102222317A (zh) 一种图像缩放方法和系统
JP2020042774A (ja) 人工知能推論演算装置
JP2017528765A (ja) ディスプレイ装置、システム及びその制御方法
WO2021013139A1 (zh) 图像处理的方法和装置
CN103747269B (zh) 一种滤波器插值方法及滤波器
CN106558021B (zh) 基于超分辨率技术的视频增强方法
TW202240535A (zh) 用於影像處理的方法及系統
CN103390266A (zh) 一种图像超分辨率方法和装置
Yuan et al. Depth map super-resolution via low-resolution depth guided joint trilateral up-sampling
Zou et al. Edgeformer: Edge-Enhanced Transformer for High-Quality Image Deblurring
JP7316130B2 (ja) 画像処理装置、画像処理装置の制御方法、プログラム
US9626565B2 (en) Automated conversion of two-dimensional hydrology vector models into valid three-dimensional hydrology vector models
CN101888557A (zh) 基于梯度域变换的卡通图像及视频上采样方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20181002

RJ01 Rejection of invention patent application after publication