- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字图像处理第三专题报告
中国地质大学(武汉)
数字图像处理上机实习
(第三专题)
姓 名:
班级序号:
学 号:
指导老师:傅华明
一 、实验题目
因为我座次号为13,所以做第4小题
一、计算图象的傅氏变换频谱函数
4,设计图象f4(x,y) 为30*120/256*256,居中垂直排列;
二、计算显示图象的频谱函数
4,对p03-04的离散哈达玛变换,显示其频谱函数
三、设计图象的频域滤波
4,根据频率采样法设计一个带阻滤波器,对图象f4(x,y)进行带阻滤波,观察分析空域图象和频谱分布的变化(设计2个阻带)。
四、实际图象的频域滤波
4,根据频率采样法设计一个带阻滤波器,对图象p03-04进行带阻滤波,观察分析空域图象和频谱分布的变化(设计2个阻带)。
实验内容
一、计算图象的傅氏变换频谱函数:设计图象f4(x,y) 为30*120/256*256,居中垂直排列
算法设计
程序代码
f = zeros(256,256);
f(68:188,113:143) = 1;
subplot(121);
imshow(f);
title(原图);
F = fft2(f);
F2=fftshift(F);
F2 = log(abs(F2));
subplot(122);
imshow(F2);
title(频谱图);
运行结果
二、计算显示图象的频谱函数:对p03-04的离散哈达玛变换,显示其频谱函数
(1)程序代码
I=imread(E:\第03专题-频域增强\P03-04.tif);
H=hadamard(256);
I=double(I)/256;
hI=H*I*H;
hI=hI/256;
subplot(121);imshow(I);
title(原图);
subplot(122);imshow(hI);
title(图像的离散哈达玛变换);
运行结果
三、设计图象的频域滤波:根据频率采样法设计一个带阻滤波器,对图象f4(x,y)进行带阻滤波,观察分析空域图象和频谱分布的变化(设计2个阻带)
算法设计
程序代码
f = zeros(256,256);
f(68:188,113:143) = 1;
subplot(221);
imshow(f);
title(原图);
F = fft2(f);
F2=fftshift(F);
subplot(222);
imshow(log(abs(F2)));
title(频谱图);
[M,N]=size(F2);
D0=50;
D1=90;
m=fix(M/2);
n=fix(N/2);
for i=1:M
for j=1:N
D=sqrt((i-m)^2+(j-n)^2);
if(D=D0)
h=1;
else if(D=D1D=D0)
h=0;
else h=1;
end
end
Q(i,j)=h*F2(i,j);
end
end
subplot(223),imshow(Q);
title(理想带阻滤波器);
Q=ifftshift(Q);
J1=ifft2(Q);
J2=uint8(real(J1));
subplot(224),imshow(J2);
title(复原图);
(3)运行结果
实际图象的频域滤波:根据频率采样法设计一个带阻滤波器,对图象p03-04进行带阻滤波,观察分析空域图象和频谱分布的变化(设计2个阻带)
算法设计
同上一题
(2)程序代码
f = imread(E:\第03专题-频域增强\P03-04.tif);
subplot(221);
imshow(f);
title(原图);
f1=double(f);
F = fft2(f1);
F2=fftshift(F);
subplot(222);
imshow(log(abs(F2)),[-1 10]);
title(频谱图);
[M,N]=size(F2);
D0=50;
D1=90;
m=fix(M/2);
n=fix(N/2);
for i=1:M
for j=1:N
D=sqrt((i-m)^2+(j-n)^2);
if(D=D0)
h=1;
文档评论(0)