CN103236035B - 基于无偏移双边二次b-样条插值的图像放大方法 - Google Patents

基于无偏移双边二次b-样条插值的图像放大方法 Download PDF

Info

Publication number
CN103236035B
CN103236035B CN201310103996.4A CN201310103996A CN103236035B CN 103236035 B CN103236035 B CN 103236035B CN 201310103996 A CN201310103996 A CN 201310103996A CN 103236035 B CN103236035 B CN 103236035B
Authority
CN
China
Prior art keywords
target
pixel points
interpolation
image
pixel
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.)
Expired - Fee Related
Application number
CN201310103996.4A
Other languages
English (en)
Other versions
CN103236035A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201310103996.4A priority Critical patent/CN103236035B/zh
Publication of CN103236035A publication Critical patent/CN103236035A/zh
Application granted granted Critical
Publication of CN103236035B publication Critical patent/CN103236035B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

本发明公开了一种基于无偏移双边二次B-样条插值的图像放大方法,包括:根据原始图像的尺寸和放大的倍数,确定目标图像的像素数量;针对目标图像的一个目标像素点,确定其在原始图像中的位置;选择原图像中距离该位置最近的N×N个像素作为插值的采样点;把每个采样点与目标像素在同一方向上的距离带入核函数计算权重,每一行或列的N个采样点得到一个插值结果;再次把N个插值结果作为采样点,把各采样点与目标像素位置的竖直距离带入核函数计算权重,插值出最终结果,即为目标像素点的像素值,完成图像的放大。本发明能够在获得理想的插值放大效果的同时提高计算效率,并且减少图像信号的震荡。

Description

基于无偏移双边二次B-样条插值的图像放大方法
技术领域
本发明涉及数字图像处理及计算机视觉技术领域,具体涉及基于无偏移双边二次B-样条插值的图像放大方法。
背景技术
随着数字化应用的不断推广,数字图像信息在社会生产和生活中起到越来越大的作用。在可视化技术研究、医疗、航天、数字娱乐等领域对高质量画面的需求越来越大。图像显示设备的多样化以及尺寸的不断变大,对高分辨率图像的需求变得更加急切。而目前高分辨率数字图像采集的设备以及网络传输支持还远远跟不上需求,因此,对低分辨率图像放大的方法变得越来越重要,一直是计算机视觉领域研究的热点。
图像放大是图像处理的基本操作之一,它广泛应用于医学图像、遥感图像、网页制作以及一些商用图像处理软件中。图像放大即将一幅低分辨率的图像转换为高分辨率图像的一种图像处理技术,对一幅图像进行放大,实质上是对图像插值的过程。
图像放大算法的选择直接影响到放大图像的质量,所以寻找合适的算法是提高放大图像质量的关键。目前主要的图像放大方法大致可以分为两类:第一类是常规插值,包括最临近点插值、双线性插值、拉格朗日插值及三次样条插值等,这类方法是根据离散的点建立一个连续函数,用这个重建的函数求出任意位置处的函数值。第二类是利用图像中包含不同的高、低频成分的特点,经过对图像的数学统计特征的分析,采用不同的方式对图像不同部分进行插值的非线性的、移变的插值方法。
数字图像放大算法最经典和常用的算法是基于线性插值的方法。例如申请号为CN201010196205.3的专利文献公开了一种图像放大方法及装置,以双线性插值方法为基础,对于放大后得到目标图像的每一个要插入像素点的灰度值计算,都转换为要插入方向上两相邻像素点灰度值的加法运算及2的幂次项移位运算,然后根据该运算对相邻两个像素点灰度值进行加法运算及移位处理就得到每一个要插入像素点的灰度值。
基于线性插值的图像放大算法比较简单,计算效率高,是被广泛应用的一类算法。但是目前几种常用的该类算法也有比较明显的缺陷。双线性插值算法造成放大图像模糊,三次的核函数插值算法会造成图像信号的大幅震荡,且计算量比较大。
发明内容
本发明提供了一种基于无偏移双边二次B-样条插值的图像放大方法,该方法采用基于B-样条的二次线性核函数,插值效果与最好的三次核函数插值十分接近,同时减少了插值信号的震荡,并且减少了计算量。
一种基于无偏移双边二次B-样条插值的图像放大方法,包括以下实施步骤:
1)选择原始图像中的部分区域作为目标图像,根据原始图像的尺寸和放大倍数,确定目标图像的像素数量;
2)选取目标图像内的一个目标像素点,确定目标像素点在原始图像中的坐标,依据该坐标计算得到目标浮点坐标;
3)选择原始图像中距离目标像素点最近的N×N个像素作为插值的采样像素点;
4)建立核函数,计算每个采样像素点与目标浮点坐标在同一方向上的距离,并代入核函数计算得到权重,然后根据权重得到每行或列采样像素点的插值结果;
5)以步骤4)中的插值结果作为采样点,计算每个采样点和目标浮点坐标之间的距离,并代入所述的核函数计算得到第二权重,由第二权重计算得到目标像素点的值。
在所述的步骤2)中,目标像素点在原始图像中的坐标乘以原始图像的放大倍数得到目标浮点坐标。
作为优选的,所述步骤3)中的N×N个像素由目标浮点坐标的左边、右边N/2列像素以及上边、下边N/2行像素组成。
进一步优选的,所述的N×N个像素为4×4个像素,4×4个像素表示为4×4的像素矩阵,该矩阵由目标浮点坐标的左边、右边两列像素以及上边、下边两列像素组成,4列像素和4行像素的像素组合即为要选择的4×4的像素矩阵,所采用4×4的像素矩阵是根据本发明的核函数的定义域来确定,当核函数的定义域改变时,像素矩阵的选取与核函数相对应。
在所述的步骤4)中,每个采样像素点与目标浮点坐标在同一方向上的距离优选为每个采样像素点与目标浮点坐标的水平距离。
所述核函数的方程为:
H ( t ) = ( - 5 | x | 2 - | x | + 4 ) / 4 , 0 < = | x | < = 1 / 2 ; ( 3 | x | 2 - 9 | x | + 6 ) / 4 , 1 / 2 < | x | < = 1 ; ( 7 | x | 2 - 19 | x | + 12 ) / 12 , 1 < | x | < = 2 / 3 ; ( - | x | 2 + 5 | x | - 6 ) / 12 , 2 / 3 < | x | < = 2 ;
其中,在计算步骤4)中的权重时,x为步骤4)中每个采样像素点与目标像素点在同一方向上的距离;在计算步骤5)中的第二权重时,x为步骤5)中采样点和目标像素点之间的距离。
当N×N个像素为4×4个像素时,步骤4)中对每行采样像素点进行插值时,作如下处理:计算每行采样像素点与步骤2)得到的目标浮点坐标的水平距离,把该距离作为参数带入核函数,得到的值作为该采样像素点的权重,该行四个采样像素点的像素值乘以各自得到的权重然后相加得到该行采样像素点的插值结果。
还可以在步骤4)中对4×4个像素的每列采样像素点进行插值计算,处理过程如下:计算每列采样像素点与步骤2)得到的目标浮点坐标的竖直距离,把该距离作为参数带入相应的核函数,得到的值作为该采样像素点的权重,该列四个采样像素点的像素值乘以各自得到的权重然后相加得到该列采样像素点的插值结果。
采样像素点为4×4个像素时,通过所述的步骤4)得到四个插值结果,四个插值结果的竖直坐标就是步骤4)插值时该行的竖直坐标,四个插值结果的水平坐标就是步骤4)插值时该列的水平坐标,计算当前每个采样点与目标浮点坐标的距离,带入核函数计算第二权重,然后采样点的像素值乘以权重,最后相加得到目标像素点的值,即为目标像素点的像素值。
与现有技术相比,本发明的有益效果体现在:
(1)本发明使用二次的核函数达到了目前公认最好的三次核函数插值放大的效果,同时相对于三次的核函数插值,节省了计算量。
(2)本发明的插值结果减少了图像信号的震荡,放大后目标图像的像素取值在更加合理的范围。
附图说明
图1为本发明基于无偏移双边二次B-样条插值的图像放大方法的流程示意图。
图2为本发明实施例中4×4个像素采样点插值的计算示意图。
具体实施方式
如图1所示,本实施例基于无偏移双边二次B-样条插值的图像放大方法的实施步骤如下:
(1)选择原始图像中的部分区域作为目标图像,根据原始图像的尺寸和放大的倍数,确定目标图像的像素数量。
(2)针对目标图像的一个目标像素点,确定其在原始图像中的坐标,然后用目标像素点的坐标除以原始图像放大的倍数得到目标浮点坐标。
(3)如图2所示,选择原图像中水平和竖直坐标距离该位置最近的4×4个像素作为插值的采样点。根据步骤(2)中得到的目标浮点坐标,选择该坐标的左边、右边各两列像素以及上边、下边各两行像素,4列像素和4行像素的交点像素即为要选择的4×4的像素矩阵。
(4)在水平方向上,计算每个采样像素点与目标像素点的水平距离并把该距离代入核函数:
H ( t ) = ( - 5 | x | 2 - | x | + 4 ) / 4 , 0 < = | x | < = 1 / 2 ; ( 3 | x | 2 - 9 | x | + 6 ) / 4 , 1 / 2 < | x | < = 1 ; ( 7 | x | 2 - 19 | x | + 12 ) / 12 , 1 < | x | < = 2 / 3 ; ( - | x | 2 + 5 | x | - 6 ) / 12 , 2 / 3 < | x | < = 2 ;
其中,x为采样像素点与目标浮点坐标在水平方向上的距离。
由核函数计算得到的值作为每个采样像素点的权重系数,每行的四个采样像素点的像素值乘以各自得到的权重系数然后相加得到该行采样像素点的插值结果。
在本实施例中,先计算每个采样像素点与目标像素点的水平距离,经核函数计算,得到一行每个采样像素点的权重系数,由该行每个采样像素点的像素值乘以各自得到的权重系数然后相加得到该行采样像素点的插值结果,在4×4的像素矩阵中,产生四个插值结果。
本实施例还可以先计算每个采样像素点与目标像素点的竖直距离,经核函数计算,得到一列每个采样像素点的权重系数,由该列每个采样像素点的像素值乘以各自得到的权重系数然后相加得到该列采样像素点的插值结果,在4×4的像素矩阵中,也产生四个插值结果。
综上所述,在计算插值结果时,可以根据每列采样像素点与目标像素点的竖直距离获得,也还可以根据每行采样像素点与目标像素点的水平距离获得。
(5)重新以步骤4)中算得的四个插值结果作为采样点,计算每个采样点和目标浮点坐标之间的竖直距离,再代入核函数:
H ( t ) = ( - 5 | x | 2 - | x | + 4 ) / 4 , 0 < = | x | < = 1 / 2 ; ( 3 | x | 2 - 9 | x | + 6 ) / 4 , 1 / 2 < | x | < = 1 ; ( 7 | x | 2 - 19 | x | + 12 ) / 12 , 1 < | x | < = 2 / 3 ; ( - | x | 2 + 5 | x | - 6 ) / 12 , 2 / 3 < | x | < = 2 ;
其中,x为采样点与目标浮点坐标在竖直方向上的距离。
由核函数计算得到的值作为每个采样点的权重,再把四个采样点的像素值乘以各自得到的权重,然后相加即可得到目标像素点的像素值,将该像素值输入目标像素点,完成图像的放大。

Claims (6)

1.一种基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,包括以下实施步骤:
1)选择原始图像中的部分区域作为目标图像,根据原始图像的尺寸和放大倍数,确定目标图像的像素数量;
2)选取目标图像内的一个目标像素点,确定目标像素点在原始图像中的坐标,依据该坐标计算得到目标浮点坐标;
3)选择原始图像中距离目标像素点最近的N×N个像素作为插值的采样像素点;
4)建立核函数,计算每个采样像素点与目标浮点坐标在同一方向上的距离,并代入核函数计算得到权重,然后根据权重得到每行或列采样像素点的插值结果;
5)以步骤4)中的插值结果作为采样点,计算每个采样点和目标浮点坐标之间的距离,并代入所述的核函数计算得到第二权重,由第二权重计算得到目标像素点的值;
所述核函数的方程为:
H ( t ) = ( - 5 | x | 2 - | x | + 4 ) / 4 , 0 < = | x | < = 1 / 2 ; ( 3 | x | 2 - 9 | x | + 6 ) / 4 , 1 / 2 < | x | < = 1 ; ( 7 | x | 2 - 19 | x | + 12 ) / 12 , 1 < | x | < = 2 / 3 ; ( - | x | 2 + 5 | x | - 6 ) / 12 , 2 / 3 < | x | < = 2 ;
其中,x为步骤4)中每个采样像素点与目标像素点在同一方向上的距离或为步骤5)中采样点和目标像素点之间的距离。
2.如权利要求1所述的基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,利用所述步骤2)中的坐标除以原始图像的放大倍数,得到所述的目标浮点坐标。
3.如权利要求2所述的基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,所述步骤3)中的N×N个像素由目标浮点坐标的左边、右边N/2列像素以及上边、下边N/2行像素组成。
4.如权利要求3所述的基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,所述的N×N个像素为4×4个像素。
5.如权利要求1所述的基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,在步骤4)中,每行或列采样像素点的像素值乘以各自得到的权重,然后相加得到所述的每行或列采样像素点的插值结果。
6.如权利要求5所述的基于无偏移双边二次B-样条插值的图像放大方法,其特征在于,在步骤5)中,每个采样点的像素值乘以各自对应的第二权重,并相加得到所述的目标像素点的值。
CN201310103996.4A 2013-03-28 2013-03-28 基于无偏移双边二次b-样条插值的图像放大方法 Expired - Fee Related CN103236035B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310103996.4A CN103236035B (zh) 2013-03-28 2013-03-28 基于无偏移双边二次b-样条插值的图像放大方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310103996.4A CN103236035B (zh) 2013-03-28 2013-03-28 基于无偏移双边二次b-样条插值的图像放大方法

Publications (2)

Publication Number Publication Date
CN103236035A CN103236035A (zh) 2013-08-07
CN103236035B true CN103236035B (zh) 2015-12-02

Family

ID=48884074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310103996.4A Expired - Fee Related CN103236035B (zh) 2013-03-28 2013-03-28 基于无偏移双边二次b-样条插值的图像放大方法

Country Status (1)

Country Link
CN (1) CN103236035B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107240064B (zh) * 2016-03-28 2021-04-06 展讯通信(天津)有限公司 图像校正方法、装置及终端
CN106204444B (zh) * 2016-06-29 2019-06-25 青岛海信移动通信技术股份有限公司 一种图像放大的方法和装置
CN109658332B (zh) * 2017-10-10 2023-03-28 凌云光技术股份有限公司 一种基于流式数据处理的图像放大方法及终端
CN108596835B (zh) * 2018-05-04 2022-04-26 嘉兴南湖学院 一种圆形区域的图像放大方法
CN111161145B (zh) * 2019-12-18 2023-12-22 兰州恒达彩印包装有限责任公司 一种改进的双三次图像插值模型的参数优化方法
CN111179169B (zh) * 2019-12-18 2023-10-27 兰州恒达彩印包装有限责任公司 一种双三次图像插值模型的参数优化方法
CN113160051B (zh) * 2021-03-30 2023-12-22 珠海全志科技股份有限公司 基于边缘方向的图像插值采样方法及装置
GB2617220A (en) * 2022-09-30 2023-10-04 Imagination Tech Ltd Upsampling blocks of pixels

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102098509A (zh) * 2010-11-19 2011-06-15 浙江大学 基于Farrow结构的可重构插值滤波器
CN101163252B (zh) * 2007-11-27 2011-10-26 中国科学院计算技术研究所 一种多媒体视频图像的缩放方法
CN102750668A (zh) * 2012-07-04 2012-10-24 西南交通大学 结合局部方向特征的数字图像三倍插值放大方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100611179B1 (ko) * 2004-06-23 2006-08-10 삼성전자주식회사 영상처리장치

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101163252B (zh) * 2007-11-27 2011-10-26 中国科学院计算技术研究所 一种多媒体视频图像的缩放方法
CN102098509A (zh) * 2010-11-19 2011-06-15 浙江大学 基于Farrow结构的可重构插值滤波器
CN102750668A (zh) * 2012-07-04 2012-10-24 西南交通大学 结合局部方向特征的数字图像三倍插值放大方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《图像缩放及其GPU实现》;夏海宏;《中国优秀硕士学位论文全文数据库信息科技辑(月刊)》;20110415(第04期);论文第37-39页 *
《边缘自适应图像缩放算法研究及VLSI实现》;叶森;《中国优秀硕士学位论文全文数据库信息科技辑(月刊)》;20120615(第06期);论文第6-11页 *

Also Published As

Publication number Publication date
CN103236035A (zh) 2013-08-07

Similar Documents

Publication Publication Date Title
CN103236035B (zh) 基于无偏移双边二次b-样条插值的图像放大方法
Wu et al. An effective method to compute the box-counting dimension based on the mathematical definition and intervals
CN107644398B (zh) 图像插补方法及其相关图像插补装置
CN103020897B (zh) 基于多区块的单帧图像的超分辨率重建的装置、系统和方法
CN102881000B (zh) 一种视频图像的超分辨率方法、装置和设备
CN108986033A (zh) 基于双线性插值算法处理的图像缩放方法
TWI539406B (zh) 影像插補方法以及應用影像插補方法之影像插補裝置與影像裝置
CN103366342B (zh) 应用于视频图像放大的分段线性插值方法
CN103854268A (zh) 基于多核高斯过程回归的图像超分辨重建方法
CN103136752B (zh) 基于边缘提取的图像放大方法
CN102186044A (zh) 一种边缘相关性图像无极缩放算法及其硬件实现装置
EP4207051A1 (en) Image super-resolution method and electronic device
CN107680043A (zh) 基于图模型的单幅图像超分辨率输出方法
CN106780336B (zh) 一种图像缩小方法及装置
CN104166992A (zh) 基于网格变形的内容感知双目图像缩放方法
CN107038199B (zh) 一种绘图方法和装置
CN100563295C (zh) 一种基于线性插值运算中误差估计的图像放大方法
Huang et al. A visual detection method for PCI blockage and coke size distribution in tuyere raceway
CN110895790B (zh) 基于后验降质信息估计的场景图像超分辨方法
KR101668716B1 (ko) 사분위수 그래프에 기반하여 데이터를 표시하는 가시화 방법 및 시스템
CN109325909A (zh) 一种图像放大方法和图像放大装置
CN104318236B (zh) 一种获取图像局部特征的方法及系统
CN105590294A (zh) 一种图像处理方法和电子设备
TWI406187B (zh) 快速高畫質視頻影像內插方法及裝置
CN108537810B (zh) 一种改进的Zernike矩亚像素边缘检测方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151202

Termination date: 20210328

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