CN101163190A - 一种基于线性插值运算中误差估计的图像放大方法 - Google Patents
一种基于线性插值运算中误差估计的图像放大方法 Download PDFInfo
- Publication number
- CN101163190A CN101163190A CNA2007102023544A CN200710202354A CN101163190A CN 101163190 A CN101163190 A CN 101163190A CN A2007102023544 A CNA2007102023544 A CN A2007102023544A CN 200710202354 A CN200710202354 A CN 200710202354A CN 101163190 A CN101163190 A CN 101163190A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- interpolation
- mfrac
- math
- 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
Links
- 230000003321 amplification Effects 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000003199 nucleic acid amplification method Methods 0.000 title claims abstract description 17
- 230000000694 effects Effects 0.000 abstract description 5
- 230000007850 degeneration Effects 0.000 abstract 1
- 238000012216 screening Methods 0.000 abstract 1
- 230000007704 transition Effects 0.000 description 5
- 230000014509 gene expression Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 101100127891 Caenorhabditis elegans let-4 gene Proteins 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明属于数字图像处理与视频显示技术领域,具体涉及基于线性插值运算的图像放大方法。本发明提出一种基于线性插值运算的图像放大方法,能较好的保留图像高频信息,具体包括以下步骤:a.计算插值点P的位置;b.获取插值点P所处位置的相邻4个像素点的像素值;c.通过相邻4个像素点的像素值计算线性插值的误差估计,将所述误差估计作为线性插值结果的补偿;d.计算线性插值结果与误差估计之和,其结果为插值点P的像素值。本发明克服了原线性插值运算导致高频退化、图像模糊的缺陷。通过本发明处理后的图像放大效果基本达到了双三次插值算法处理后的图像效果,且运算量小,易于硬件实现。
Description
技术领域
本发明属于数字图像处理与视频显示技术领域,具体涉及基于线性插值运算的图像放大方法。
背景技术
图像放大是数字图像处理中最重要的技术之一,常用的放大方法有最邻近插值、双线性插值、双三次插值等。最邻近插值简单、速度快,但放大后的图像有严重的马赛克现象;双线性插值会退化图像的高频部分,造成图像的模糊;双三次插值算法较好的保留图像高频信息,边缘更锐化、细节更清晰,但运算量和硬件实现难度都很大。
发明内容
本发明所要解决的技术问题是,提出一种基于线性插值运的图像放大方法,能较好的保留图像高频信息,使边缘更锐化、细节更清晰。
本发明为解决上述技术问题所采用的技术方案是,一种基于线性插值运算中误差估计的图像放大方法,其特征在于,具体包括以下步骤:
a、计算插值点P的位置;
b、获取插值点P所处位置的相邻4个像素点的像素值;
c、通过相邻4个像素点的像素值计算线性插值的误差估计,将所述误差估计作为线性插值结果的补偿;
d、计算线性插值结果与误差估计之和,其结果为插值点P的像素值。
所述相邻4个像素点,具体的是,在做水平方向的插值时为水平位置的4个相邻像素点;做垂直方向的插值时为垂直位置的4个相邻像素点;
所述误差估计具体为,在做水平方向的插值时误差估计为:
其中插值点P水平方向相邻4点的像素值分别为f(Pi-1)、f(Pi)、f(Pi+1)、f(Pi+2);x为插值点P与相邻像素点Pi的水平距离,x∈(0,1);那么,水平方向做插值时插值点P的像素值f(P)为:
在做垂直方向的插值时误差估计为: 其中插值点P垂直方向相邻4点的像素值分别为f(Pj-1)、f(Pj)、f(Pj+1)、f(Pj+2);y为插值点P与相邻像素点Pi的垂直距离,y∈(0,1);那么,垂直方向做插值时插值点P的像素值f(P)为:
本发明的有益效果是,原双线性插值方法的插值点由原图像周围4个像素点的值确定,通过本发明引入误差估计参变量的改进后,插值点像素值的确定参考了该点周围16个像素点的值。本发明克服了原线性插值运算导致高频退化、图像模糊的缺陷。通过本发明处理后的图像放大效果基本达到了双三次插值算法处理后的图像效果,且运算量小,易于硬件实现。
附图说明
图1为本发明原理图;
图2为误差估计的推导示意图。
具体实施方式
线性插值中相对于三次插值的误差估计E的计算原理如下:
线性插值如图1所示,Pi与Pi+1为插值点P在水平方向上的相邻2像素点,线性插值的公式f′(P)=f(Pi)+[f(Pi+1)-f(Pi)]·x,x∈(0,1);线性插值不光滑,会造成图像高频部分的退化。要保留图像高频信息,需要插值点P与像素点形成的曲线尽量光滑,所以引入误差估计E,完善插值点P的像素值的计算方法:
f(P)=f(Pi)+[f(Pi+1)-f(Pi)]·x+E,x∈(0,1);
误差估计E的推导方法如下:
牛顿一次插值多项式为:
其中, 为牛顿一次插值多项式的误差估计;
如图2所示,-1、0、1、2分别表示插值点x在水平方向上相邻4点的位置,若:x0=-1,h=3,代入(1)式,得到:
若:x0=0,h=1,代入(1)式,得到:
在本发明中作误差估计时,认为f″(x)在x∈(-1,2)上变换很小,即认为:f″(ξ1)≈f″(ξ2),并统一用f″(σ)表示,并由式(2)、式(3)两式相减得到:
令x=1/2,得到:
当x=1/2时,如图2所示,A点的值为:
f(0)+[f(1)-f(0)]·(1/2);B点的值为:
当x=1/2时,得出结论,误差估计E=AB/8。
再用二次函数并引入误差估计值E,来近似估计error曲线ferror(x),如图2所示得到:f1(x)=-4·E·x2+(4·E-a)·x;f2(x)=-ax,error曲线为曲线f1(x)与直线f2(x)之差:因此:ferror(x)=4·E·x·(1-x);
将E=AB/8带入上式,得到
E为ferror(x)的估计值,令插值点水平方向相邻4点分别为Pi-1、Pi、Pi+1、Pi+2,那么有水平方向插值时的误差估计值E为:
同理,令插值点水平方向相邻4点分别为Pj-1、Pj、Pj+1、Pj+2,由上面的推导可知垂直方向上插值时的误差估计值E为:
实施例
在实现图像放大过程中,本发明先对图像在水平方向作放大,得到一个过渡图像,然后再对过渡图像作垂直方向的放大,得到最终的放大图像。当然也可以先做垂直方向的放大再做水平方向的放大。放大的实质是将二维信号作一维方向上的插值运算。
令原始图像大小为M×N,放大后的图像大小为X×Y,并且第一行(列)以零行(列)开始计数,最后一行(列)为M-1行(N-1列)。
首先作水平方向的插值运算,具体步骤如下:
1、对于原始图像第r行(r=0,1,2……M-1),依次计算出插值点P的水平位置s:s=-0.5+N/(2*Y)+C*N/Y;其中,C=0,1,2……Y-1为图像放大后的列序数。
由此可以找到插值点P所处位置的前后相邻四个已知点:Pi-1=floor(s)-1、Pi=floor(s)、Pi+1=floor(s)+1、Pi+2=floor(s)+2;其中,函数floor(s)为s向下取整。
2、计算出P到Pi的距离x:x=s-floor(s);
3、根据本发明提出的插值公式,计算出插值点P的像素值:
4、依次完成0至M-1行的水平方向插值运算,并由此完成整个图像在水平方向的放大,得到一个过渡图像T,其大小为M×Y。
再进行垂直方向的插值,具体步骤如下:
5、对于过渡图像T的第C列,C=0,1,2……Y-1,在垂直方向上,依次计算出插值点P的垂直位置t为:t=-0.5+M/(2*X)+R*M/X,其中,R=0,1,2……X-1,为图像放大后的行序数;
得到待插点P所处位置的上下相邻四个已知点:Pj-1=floor(t)-1、Pj=floor(t)、Pj+1=floor(t)+1、Pj+2=floor(t)+2;
6、计算出P到Pj的距离y:y=t-floor(t);
7、根据本发明提出的插值公式,计算出待插点P的像素值:
8、依次完成0至Y-1列的垂直方向插值运算,并由此完成整个过渡图像T在垂直方向的放大,得到最终放大的图像,其大小为X×Y。
通过实验验证:经本发明处理过的图像比经双线性插值处理后的图像有明显改善,总体放大效果与双三次插值相当,但运算量远低于双三次插值运算。
Claims (3)
1.一种基于线性插值运算中误差估计的图像放大方法,其特征在于,具体包括以下步骤:
a、计算插值点P的位置;
b、获取插值点P所处位置的相邻4个像素点的像素值;
c、通过相邻4个像素点的像素值计算线性插值的误差估计,将所述误差估计作为线性插值结果的补偿;
d、计算线性插值结果与误差估计之和,其结果为插值点P的像素值。
2.如权利要求1所述一种基于线性插值运算中误差估计的图像放大方法,其特征在于,步骤b具体为:在做水平方向的插值时,获取插值点P水平位置的相邻4个像素点的像素值;做垂直方向的插值时,获取插值点P垂直位置的相邻4个像素点的像素值。
3.如权利要求2所述一种基于线性插值运算中误差估计的图像放大方法,其特征在于,步骤c所述误差估计具体为,在做水平方向的插值时误差估计为:
其中插值点P水平方向的相邻4个像素点的像素值分别为f(Pi-1)、f(Pi)、f(Pi+1)、f(Pi+2);x为插值点P与相邻像素点Pi的水平距离,x∈(0,1);步骤d中,水平方向做插值时插值点P的像素值f(P)为:
其中插值点P垂直方向相邻4个像素点的像素值分别为f(Pj-1)、f(Pj)、f(Pj+1)、f(Pj+2);y为插值点与相邻像素点Pj的垂直距离,y∈(0,1);步骤d中,垂直方向做插值时插值点P的像素值f(P)为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007102023544A CN100563295C (zh) | 2007-11-01 | 2007-11-01 | 一种基于线性插值运算中误差估计的图像放大方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007102023544A CN100563295C (zh) | 2007-11-01 | 2007-11-01 | 一种基于线性插值运算中误差估计的图像放大方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101163190A true CN101163190A (zh) | 2008-04-16 |
CN100563295C CN100563295C (zh) | 2009-11-25 |
Family
ID=39298010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007102023544A Expired - Fee Related CN100563295C (zh) | 2007-11-01 | 2007-11-01 | 一种基于线性插值运算中误差估计的图像放大方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100563295C (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101437137B (zh) * | 2008-12-19 | 2010-08-25 | 四川虹微技术有限公司 | 场内插值方法 |
CN102263924A (zh) * | 2010-05-29 | 2011-11-30 | 比亚迪股份有限公司 | 一种基于双三次插值的图像处理方法及其图像显示方法 |
CN102760281A (zh) * | 2011-04-26 | 2012-10-31 | 撖龙 | 图像缩放方法 |
CN101815157B (zh) * | 2009-02-24 | 2013-01-23 | 虹软(杭州)科技有限公司 | 图像及视频的放大方法与相关的图像处理装置 |
CN103366342A (zh) * | 2013-07-02 | 2013-10-23 | 天津大学 | 应用于视频图像放大的分段线性插值方法 |
CN104184981A (zh) * | 2014-08-27 | 2014-12-03 | 深圳市华星光电技术有限公司 | 一种基于缩减像素采样处理的低解析度显示方法及设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR950006776B1 (ko) * | 1993-01-14 | 1995-06-22 | 삼성전자주식회사 | 디지탈 영상 데이타의 보간방법 및 회로 |
CN100365660C (zh) * | 2004-12-13 | 2008-01-30 | 北京中星微电子有限公司 | 一种图像放大插值的方法 |
-
2007
- 2007-11-01 CN CNB2007102023544A patent/CN100563295C/zh not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101437137B (zh) * | 2008-12-19 | 2010-08-25 | 四川虹微技术有限公司 | 场内插值方法 |
CN101815157B (zh) * | 2009-02-24 | 2013-01-23 | 虹软(杭州)科技有限公司 | 图像及视频的放大方法与相关的图像处理装置 |
CN102263924A (zh) * | 2010-05-29 | 2011-11-30 | 比亚迪股份有限公司 | 一种基于双三次插值的图像处理方法及其图像显示方法 |
CN102263924B (zh) * | 2010-05-29 | 2014-05-28 | 比亚迪股份有限公司 | 一种基于双三次插值的图像处理方法及其图像显示方法 |
CN102760281A (zh) * | 2011-04-26 | 2012-10-31 | 撖龙 | 图像缩放方法 |
CN103366342A (zh) * | 2013-07-02 | 2013-10-23 | 天津大学 | 应用于视频图像放大的分段线性插值方法 |
CN103366342B (zh) * | 2013-07-02 | 2015-12-23 | 天津大学 | 应用于视频图像放大的分段线性插值方法 |
CN104184981A (zh) * | 2014-08-27 | 2014-12-03 | 深圳市华星光电技术有限公司 | 一种基于缩减像素采样处理的低解析度显示方法及设备 |
CN104184981B (zh) * | 2014-08-27 | 2017-12-15 | 深圳市华星光电技术有限公司 | 一种基于缩减像素采样处理的低解析度显示方法及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN100563295C (zh) | 2009-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100563295C (zh) | 一种基于线性插值运算中误差估计的图像放大方法 | |
CN102254301B (zh) | 基于边缘方向插值的cfa图像去马赛克方法 | |
JP4874904B2 (ja) | 画像処理装置及びその方法 | |
CN106169173B (zh) | 一种图像插值方法 | |
JP2011504682A (ja) | 画像シーケンスのサイズ変更 | |
CN107341765A (zh) | 一种基于卡通纹理分解的图像超分辨率重建方法 | |
US9105106B2 (en) | Two-dimensional super resolution scaling | |
JP2009212969A (ja) | 画像処理装置、画像処理方法、及び画像処理プログラム | |
CN103595981A (zh) | 基于非局部低秩的色彩滤波阵列图像去马赛克方法 | |
CN101499164A (zh) | 一种基于单张低分辨率图像的图像插值重构方法 | |
Mahajan et al. | Adaptive and non-adaptive image interpolation techniques | |
CN107392854A (zh) | 一种基于局部自适应增益因子的联合上采样方法 | |
CN104134189B (zh) | 一种图像放大的方法及装置 | |
CN105279742B (zh) | 一种快速的基于分块噪声能量估计的图像去噪方法 | |
JP2005217532A (ja) | 解像度変換方法及び解像度変換装置 | |
JP4868249B2 (ja) | 映像信号処理装置 | |
US20040086201A1 (en) | Fast edge directed polynomial interpolation | |
Ousguine et al. | A new image interpolation using gradient-orientation and cubic spline interpolation | |
JP4688753B2 (ja) | 補間処理装置および補間処理方法 | |
TWI406187B (zh) | 快速高畫質視頻影像內插方法及裝置 | |
US8346021B2 (en) | Content adaptive scaler based on a farrow structure | |
CN102930503B (zh) | 基于cfa图像多方向边缘插值的方法 | |
CN101459811B (zh) | 视频画面格式转换方法和对应设备 | |
CN102930504B (zh) | 一种保持边缘特征的数字图像放大方法 | |
JP2005057529A (ja) | 画像処理方法、画像処理装置およびプログラム |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091125 Termination date: 20161101 |