实验二 图像几何变换与傅里叶变换.doc

实验二图像几何变换与傅里叶变换
课程名称数字图像处理实验名称图像几何变换与傅里叶变换
姓名__吴征宇学号____ 3101110002 专业班级_
实验日期__2013 年_ _4月_24日成绩_____ __ 指导教师___ _
一、实验目的与要求
、缩放、旋转与镜像变换;
。
,观察图像的频域表示的特点
二、实验例程
1、几何变换
图像平移,
function outImage=immove(inImage,Tx,Ty,i)
[m, n,i] = size(inImage);%计算图象大小
outImage = uint8(zeros(m+abs(Ty), n+abs(Tx),i));
if(Tx>0&&Ty>0) %往右下角平移的情况
outImage(1+Ty:m+Ty,1+Tx:n+Tx,1:i) = inImage;
else if(Tx<0&&Ty<0) %往左上角平移的情况
outImage(1:m,1:n,1:i) = inImage;
else if (Tx<0&&Ty>0) %往左下角平移的情况
outImage(1+Ty:m+Ty,1:n, 1:i) = inImage;
else % 往右上角平移的情况
outImage(1:m,1+Tx:n+Tx, 1:i) = inImage;
end
end
end
傅里叶变换
matlab计算和形象化2DDFT
傅立叶变换函数和逆向傅立叶变换函数
F=fft2(f) f=real(ifft2(F))
例子:I=imread('');
f=real(ifft2(F));
F=fft2(I);
imshow(F,[])
f1=real(ifft2(F));
imshow(f,[])
傅立叶频谱计算函数
s=abs(F)
例子:频谱可视化
I=imread('');
F=fft2(I);
S=abs(F);
imshow(S,[])
中心频谱函数: 把变换的原点移到频域区域中心
Fc=fftshift(F)
把变换的频域区域中心移到原点:
F=ifftshift(F)
例子:中心频谱可视化
I=imread('');
F=fft2(I);
Fc=fftshift(F);
S=abs(Fc);
imshow(S,[])
f2=ifftshift(Fc);
imshow(abs(f),[])
频谱动态范围压缩(频谱增强):
例子: I=imread('');
F=fft2(I);
S=abs(Fc);
S2=log(1+abs(Fc));
imshow(S2,[])
三、实验内容及步骤
***像,实现图像的平移、缩放、旋转和平移
I=imread('d:\');
H=immove(I,80,-80);
figure;imshow(H);title('平移后图像');
B=immove(I,80,-50);
figure;imshow(B);title('平移后图像');
G=imresize(I,);
figure;imshow(G);title('缩放后图像');
C=imresize(I,2);
figure;imshow(C);title('缩放后图像');
D=imrotate(I,30);
figure;imshow(D);title('旋转后图像');
E=imrotate(I,54);
figure;imshow(E);title('旋转后图像');
F = flipdim(I,1);
figure;imshow(F);title('垂直镜像图像');
G = flipdim(I,2);
figure;imshow(G);title('水平镜像图像')
imcrop(I,[80,80,80,80])
imcrop(I,[45,75,65,61])
2、读入一幅灰度图像,实现图像的平移、缩放、旋转和平移
J=rgb2gray(I);
imshow(J);title('原图像');
H=immove(J,80,-80);
figure;imshow(H);title('平移后图像');
B=immove(J,50,-50);
figure;imshow(B);title('平移后图像');
G=imresize(J,);
figure;imsho