matlab图像双线性插值,最近邻插值与几何变换

1、实验二图像几何变换与插值、实验目的巩固图像处理编程的步骤格式，理解数据插值及图像几何变换的原理，掌握图像几何变换的实现方法。二、实验内容1、理解反向变换的实现思路2、图像缩放及插值Matlab取整命令：floor,ceil,round分别实验最近邻插值和双线性插值f=imread(lena.bmp);beishu=0.5;row,col=size(f);r1=round(row*beishu);c1=round(col*beishu);b=zeros(r1,c1);fori=1:r1forj=1:c1i1=round(i/beishu);j1=round(j/beishu);Difi11i1=

2、1;endifj11j1=1;endb(i,j)=f(i1,j1);endendb=uint8(b);figure;imshow(f);figureimshow(b);3、图像旋转及插值以图像中心为轴实现任意角度(逆时针为正)的图像旋转，分别实验两种插值算法f=imread(lena.bmp);B=imrotate(f,45,nearest,crop);C=imrotate(f,45,bilinear,crop);figure;subplot(121);imshow(f);title(原图像);subplot(122);imshow(B);title(最近邻插值);figure;subplot

3、(121);imshow(f);title(原图像);subplot(122);imshow(C);title(双线性插值);原图像最近邻插值处理原图像双线性插值处理7%.7%，最邻近插值）%matlab练习程序F魁ib%9（图像放大/缩小%变换矩阵x=h*u,y=w*vccccl;w=0.6;h=1.4;%放大或缩小的宽度%放大或缩小的高度img=imread(Corner.png);imshow(img);mn=size(img);imgn=zeros(h*m,w*n);rot=h00;0w0;001;inv_rot=inv(rot);forx=1:h*mfory=1:w*npix=xy1

4、*inv_rot;imgn(x,y)=img(round(pix(1),round(pix(2);endendfigure,imshow(uint8(imgn)%matlab练习程序(图像放大/纟宿小，双线性插值)cccm=1.8;%放大或宿小的高度n=2.3;%放大或宿小的宽度img=imread(lena.jpg);imshow(img);hw=size(img);imgn=zeros(h*m,w*n);rot=m00;0n0;001;%变换矩阵fori=1:h*mforj=1:w*npix=ij1/rot;float_Y=pix(1)-floor(pix(1);float_X=pix(2

5、)-floor(pix(2);ifpix(1)hpix(1)=h;endifpix(2)wpix(2)=w;endpix_up_left=floor(pix(1)floor(pix(2);%四个相邻的点pix_up_right=floor(pix(1)ceil(pix(2);pix_down_left=ceil(pix(1)floor(pix(2);%计算临近四个点的权重pix_down_right=ceil(pix(1)ceil(pix(2);value_up_left=(1-float_X)*(1-float_Y);value_up_right=float_X*(1-float_Y);va

6、lue_down_left=(1-float_X)*float_Y;value_down_right=float_X*float_Y;%按权重进行双线性插值imgn(i,j)=value_up_left*img(pix_up_left(1),pix_up_left(2)+.value_up_right*img(pix_up_right(1),pix_up_right(2)+.value_down_left*img(pix_down_left(1),pix_down_left(2)+.value_down_right*img(pix_down_right(1),pix_down_right(2)

7、;endendfigure,imshow(uint8(imgn)%matlab练习程序（图像旋转，ccc最邻近插值）cl;H=1;W=2;jiaodu=45;img=imread(Corner.png);imshow(img);vu=size(img);%索引pix中第一个兀素，即咼度%索引pix中第二个元素，即宽度%要旋转的角度，旋转方向为顺时针%这里V为原图像的咼度，u为原图像的宽度%这里y为变换后图像的高度，x为变换后图像的宽度theta=jiaodu/180*pi;rot=cos(theta)-sin(theta)0;sin(theta)cos(theta)0;001;inV_rot=

8、inV(rot);pix1=111*rot;pix2=1u1*rot;pix3=V11*rot;pix4=Vu1*rot;%变换后图像左上点的坐标%变换后图像右上点的坐标%变换后图像左下点的坐标%变换后图像右下点的坐标height=round(max(abs(pix1(H)-pix4(H)+0.5abs(pix2(H)-pix3(H)+0.5);%变换后图像的高度width=round(max(abs(pix1(W)-pix4(W)+0.5abs(pix2(W)-pix3(W)+0.5);%变换后图像的宽度%取得y方向的负轴超出的偏移量%取得x方向的负轴超出的偏移量imgn=zeros(hei

9、ght,width);delta_y=abs(min(pix1(H)pix2(H)pix3(H)pix4(H);delta_x=abs(min(pix1(W)pix2(W)pix3(W)pix4(W);fory=1-delta_y:height-delta_yforx=1-delta_x:width-delta_xpix=yx1*inV_rot;%用变换后图像的点的坐标去寻找原图像点的坐标，%否则有些变换后的图像的像素点无法完全填充ifpix(H)=0.5&pix(W)=0.5&pix(H)=v&pix(W)=l&pix(2)=l&pix(l)=h&pix(2)=wpix_up_left=fl

10、oor(pix(l)floor(pix(2);%四个相邻的点pix_up_right=floor(pix(l)ceil(pix(2);pix_down_left=ceil(pix(l)floor(pix(2);pix_down_right=ceil(pix(l)ceil(pix(2);value_up_left=(l-float_X)*(l-float_Y);%计算临近四个点的权重value_up_right=float_X*(l-float_Y);value_down_left=(l-float_X)*float_Y;value_down_right=float_X*float_Y;imgn

11、(i+delta_y,j+delta_x)=value_up_left*img(pix_up_left(l),pix_up_left(2)+.value_up_right*img(pix_up_right(l),pix_up_right(2)+.value_down_left*img(pix_down_left(l),pix_down_left(2)+.value_down_right*img(pix_down_right(l),pix_down_right(2);endendendfigure,imshow(uint8(imgn)%实验一Matlab图像基本操作、实验目的熟悉利用Matlab进行图像处理的基本操作，了解图像数据的存储形式及进行图像处理编程的步骤方法。二、实验内容1、图像读写与显示重点函数：imread,imwrite,imshow2、彩色图像灰度化计算公式：Gray=R*0.299+G*0.587+B*0.1143、图像马赛克局部平均，改变窗口大小比较处理结果，如取2X2、4X4或更大尺寸的窗口4、图像平移分别完成图像水平方向、竖直方向和两个方向的平移三、实验要求1、编写代码，完成各项实验内容2、总结实验中遇到问题及解决方案，书写实验报告实验二图像几何变换与插值一、实验目的巩固图像处理编程的步骤格式，理解数