matlab_图像几何变换附带程序

实验四：图像几何变换（编程报告）
一、实验目的
(1)学习几种常见的图像几何变换，并通过实验体会几何变换的效果；
(2)掌握图像平移、剪切、缩放、旋转、镜像、错切等几何变换的算法原理及编
程实现
(3)掌握matlab编程环境中基本的图像处理函数
(4)掌握图像的复合变换
二、涉及知识点
(1)图像几何变换不改变图像像素的值，只改变像素所在的几何位置
(2)图像裁剪imcrop函数，语法格式为：
B=imcrop(A);交互式用鼠标选取区域进行剪切
B=imcrop(A,[left top right bottom]);针对指定的区域[left top right bottom]进行剪切
(3)图像缩放imresize函数，语法格式为：
B = imresize(A,m,method)
这里参数method用于指定插值的方法，可选用的值为'nearest'（最邻近法），'bilinear'（双线性插值），'bicubic'（双三次插值），默认为'nearest'。

B = imresize(A,m,method)返回原图A的m倍放大的图像（m小于1时效果是
缩小）。

(4)图像旋转imrotate函数，语法格式为：
B = imrotate(A,angle,’crop’)，参数crop用于指定裁剪旋转后超出图像的部分。

三、实验内容
(1)将图像hehua.bmp裁剪成200X200大小
(2)制作动画，将一幅图像逐渐向左上角平移移出图像区域，空白的地方用白色
填充
(3)利用剪切图像函数制作动画
(4)将图像分别放大1.5倍和缩小0.8倍，插值方法使用双线性插值法，分别显
示图像。

(5)将图像水平镜像，再顺时针旋转45度，显示旋转后的图。

(6)将图像分别进行水平方向30度错切，垂直方向45度错切，分别显示结果
四、实验环境
Windows下matlab编程环境
五、实验源代码及结果
1.f=imread('hehua.bmp');
figure;
imshow(f);
title('原图');
f2=imcrop(f,[50,50,250,250]);
figure;
imshow(uint8(f2));
title('裁剪后');
imwrite(f2,'d:/5/hehua1.bmp');
2.f=imread('hehua1.bmp');
[m,n,x]=size(f);
f=double(f);
for i=1:10
mx=10*i;
my=10*i;
g=zeros(m,n,x)+255;
%g(mx+1:m,my+1:n,1:x)=f(1:m-mx,1:n-my ,1:x);
g(1:m-mx,1:n-my ,1:x)=f(mx+1:m,my+1:n,1:x);
figure;
imshow(uint8(g));
end
3.f=imread('hehua1.bmp');
[m,n]=size(f);
for i=50:10:200
m=i;
n=i;
f2=imcrop(f,[n,n,m,m]);
figure;
imshow(uint8(f2));
end
4.f=imread('hehua1.bmp');
figure;
imshow(f);
title('原图');
f=double(f);
f1=imresize(f,1.5,'bilinear');
figure;
imshow(uint8(f1));
title('放大1.5倍');
f2=imresize(f,0.8,'bilinear');
figure;
imshow(uint8(f2));
title('缩小0.8倍');
5.f=imread('hehua1.bmp');
subplot(131);
imshow(f);
title('原图');
[m,n,x]=size(f);
g=zeros(m,n,x);
for i=1:m
for j=1:n
for k=1:x
g(i,j,k)=f(i,n-j+1,k);
end
end
end
subplot(132);
imshow(uint8(g));
title('水平镜像');
f2=imrotate(g,45,'crop');
subplot(133);
imshow(uint8(f2));
title('顺时针旋转45度');
6.f=imread('hehua1.bmp');
subplot(131);
imshow(f);
title('原图');
h=size(f);
f1=zeros(h(1)+round(h(2)*tan(pi/6)),h(2),h(3));
for m=1:h(1)
for n=1:h(2)
f1(m+round(n*tan(pi/6)),n,1:h(3))=f(m,n,1:h(3));
end
end
subplot(132);
imshow(uint8(f1));
title('水平30度');
f2=zeros(h(1),h(2)+round(h(2)*tan(pi/4)),h(3));
for m=1:h(1)
for n=1:h(2)
f2(m,n+round(m*tan(pi/4)),1:h(3))=f(m,n,1:h(3));
end
end
subplot(133);
imshow(uint8(f2));
title('垂直45度');
六、心得体会。