资源描述：
《数字图像插值算法分析与实现》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、一．问题的叙述插值是一种最基本、最常用的几何运算。尤其随着多媒体技术的发展以及大屏幕的图像显示，以及生物医学领域的发展，图像放大的重要性也越来越明显。它不仅应用广泛,而且插值的精度直接影响最终的图像处理结果,在图像处理软件中对图像进行缩放时,插值算法的好坏直接关系到图像的失真程序,插值函数的设计是插值算法的核心问题。常采用三种插值算法:最近邻点插值、双线性插值和双三次插值,其中双三次插值的效果最好,而且这一结论也得到了普遍的公认。本文对该三种算法进行了分析并通过matlab编程实现了算法。二．二维图像插值技术2．1最邻近插值法最邻近插值又称为像素重复放大法。它是一种最简单
2、也是最原始的图像插值放大，它的基本原理是用原始像素直接对插值点像素用原始像素填充。最近邻点插值又称零阶插值,它输出的像素值等于距离它映射到的位置最近的输入像素值。对于二维图像,该法是"取待采样点周围4个相邻像素点中距离最近的1个邻点的灰度值作为该点的灰度值。插值公式:(1)其中i、j均为非负整数,u、v为[0,1)区间的浮点数,一般可取0.5,f(i,j)表示源图像(i,j)处的的像素值。2.2双线性插值法二元函数的线性插值方法是一元函数线性插值方法的直接推广，将一维空间中的线段推广为二维空间中的矩形，由此产生二维线性插值。已知平面上一矩形域内四个定点P1,P2,P3,P
3、4处的函数值为：，，令，由此构造基函数：，，所以有差值函数为（2）2.3双三次差值法对于灰度变化复杂的图像，不能简单的采用线性插值。可以采用多项式插值技术对数据点做多项式插值，不仅考虑到四个直接邻点灰度值的影响,还考虑到各邻点间灰度值变化率的影响",利用了待采样点周围更大邻域内像素的灰度值作三次插值。多项式插值如下：。需用n+1个数据点来建立方程组，求出系数的值。用公式表达为：（4）定义 f(i+u,j+v)=[A]*[B]*[C] 三．算法实现：3.1最邻近插值法：clc;clear;%functiony=resize(a,K)%%%%%%%%%%%%%%%%%
4、%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I:输入图像%K:缩放倍数%type:1表示最邻近法，2表示双极性插值法%画出缩放后图像并返回其灰度值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I=imread('girl.bmp');%读入图像I=rgb2gray(I);K=3;[m,n]=size(I);m1=m*K;n1=n*K;%fori=1:m1forj=1:n1tic%计算插值所用的时间u=round(i/K);%取最邻近值v=round(j
5、/K);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%ifu<=0%避免取近似值时超过图像的边界u=1;endifv<=0v=1;endifu>=m1u=m1;endifv>=n1v=n1;end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%b(i,j)=I(u,v);%为新的图像赋值endendtocimshow(b);title('缩放后图像');3.2双线性插值关键部分程序：%双线性插值functionM=s1(x,y,nrow
6、s,ncols,K,I)%K表示放大倍数width=K*nrows;height=K*ncols;widthScale=nrows/width;heightScale=ncols/height;%xx，yy表示当前像素点在源图像中的坐标，x，y表示在插值图像中的坐标xx=x*widthScale;yy=y*heightScale;if(xx==double(uint16(xx)))&(yy==double(uint16(yy)))%如果xx，yy为整数，则将像素点在原图中的值I(x,y)直接赋给J(x,y)M=I(int16(xx),int16(yy));else%如果a
7、，b不是整数，则进行双线性插值a=double(uint16(xx));%xx,yy取最邻近的整数b=double(uint16(yy));ifa<=0%避免取整后超过图像大小a=1;endifb<=0b=1;endifa>=nrowsa=nrows-1;endifb>=ncolsb=ncols-1;endx11=double(I(a,b));%x11<-I(a,b)x12=double(I(a,b+1));%x12<-I(a,b+1)x21=double(I(a+1,b));%x21<-I(a+1,b)x22=double(
