数字图像处理上机作业 2.docx
- 文档编号:15125693
- 上传时间:2023-07-01
- 格式:DOCX
- 页数:21
- 大小:1.25MB
数字图像处理上机作业 2.docx
《数字图像处理上机作业 2.docx》由会员分享,可在线阅读,更多相关《数字图像处理上机作业 2.docx(21页珍藏版)》请在冰点文库上搜索。
数字图像处理上机作业2
西安电子科技大学
数字图像处理
上机作业
教师:
刘清华
学院:
通信工程学院
班级:
021141
专业:
探测制导与控制技术
姓名:
苏二鹏
学号:
02114030
(1)
(2)
f1与f2的幅度谱相同,因为f1与f2振幅的绝对值是相同的。
(3)
f3的幅度谱与f2的幅度谱按顺时针旋转90度后相同
(4)
f4的幅度谱与f1的幅度谱按顺时针旋转90度后相同
(5)
f5的幅度谱是f1与f4的幅度谱相加,f6的幅度谱是f2与f3的幅度谱相加
程序:
f1=zeros(256,256);
f1(65:
192,113:
144)=100;
figure
subplot(121)
imshow(f1)
title('原图f1')
F=fft2(f1);
subplot(122)
imshow(abs(F))
title('f1幅度谱图')
form=(1:
256);
forn=(1:
256);
f2(m,n)=power(-1,m+n)*f1(m,n);
end
end
figure
subplot(121)
imshow(f2)
title('f2图')
F1=fft2(f2);
subplot(122)
imshow(abs(F1))
title('f2幅度谱图')
f3=imrotate(f2,-90);
figure
subplot(121)
imshow(f3)
title('f3图')
F2=fft2(f3);
subplot(122)
imshow(abs(F2))
title('f3幅度谱图')
f4=imrotate(f1,-90);
figure
subplot(121)
imshow(f4)
title('f4图')
F5=fft2(f4);
subplot(122)
imshow(abs(F5))
title('f4幅度谱图')
f5=f1+f4;
figure
subplot(121)
imshow(f5)
title('f5图')
F3=fft2(f5)
subplot(122)
imshow(abs(F3))
title('f5幅度谱图')
f6=f2+f3;
figure
subplot(121)
imshow(f6)
title('f6图')
F4=fft2(f6)
subplot(122)
imshow(abs(F4))
title('f6幅度谱图')
256*256图像如下:
从中值滤波和均值滤波后的图像可以看出,中值滤波把黑色方快完全完全隔离开,而白色方块则通过边角的两个多出来的像素点相连,而均值滤波则把图像变模糊了,边缘模糊的比较厉害。
中值滤波后图像的直方图没有变化,而均值滤波后的直方图变化了。
程序:
A=zeros(256,256);
fori=(1:
64:
256);
forj=(1:
64:
256);
A(j:
j+32,i:
i+32)=1;
end
end
form=(32:
64:
256);
forn=(32:
64:
256);
A(n:
n+32,m:
m+32)=1;
end
end
figure
imshow(A)
title('原图像')
figure
K=medfilt2(A,[3,3])
subplot(121)
imshow(K)
title('3*3中值滤波')
M=filter2(fspecial('average',3),A)
subplot(122)
imshow(M)
title('3*3均值值滤波')
figure
subplot(131)
imhist(A)
axis([0202])
title('原图直方图')
subplot(132)
imhist(K)
axis([0202])
title('中值滤波后的直方图')
subplot(133)
imhist(M)
axis([0202])
title('均值滤波后的直方图')
整体而言中值滤波的结果比均值滤波的结果要好,均值滤波模糊了轮廓边缘,但是中值滤波对椒盐噪声的滤波破坏了轮廓边缘。
中值滤波对高斯噪声的水平滤波不是很好。
程序:
I=zeros(256,256);
fori=(32:
24:
224);
I(23:
233,i:
i+7)=1;
end
figure
imshow(I)
title('原图像')
A=imnoise(I,'gaussian',0,0.025);
figure
subplot(121)
imshow(A)
title('高斯白噪声')
B=imnoise(I,'salt&pepper',0.025);
subplot(122)
imshow(B)
title('椒盐噪声')
h=[111;111;111];%%领域平均法消除噪声模板为1/9[111;111;111]
h=9./h;
Ja=conv2(A,h);
figure
subplot(121)
imshow(Ja,[])
title('对高斯白噪声的均值滤波')
Jb=conv2(B,h);
subplot(122)
imshow(Jb,[])
title('对椒盐噪声的均值滤波')
Ka=medfilt2(A,[3,3])
figure
subplot(121)
imshow(Ka,[])
title('对高斯白噪声的中值滤波')
Kb=medfilt2(B,[3,3]);
subplot(122)
imshow(Kb,[])
title('对椒盐噪声的中值波')
roberts自动选择的阈值为:
0.1305
sobel自动选择的阈值为:
0.1162
prewitt自动选择的阈值为:
0.1138
程序:
I=imread('rice.png');
figure
imshow(I)
title('原始图像')
[BW1,thresh1]=edge(I,'roberts');
disp('roberts自动选择的阈值为:
')
disp(thresh1)
figure
subplot(131)
imshow(BW1)
title('roberts算子')
[BW2,thresh2]=edge(I,'sobel');
disp('sobel自动选择的阈值为:
')
disp(thresh2)
subplot(132)
imshow(BW2)
title('sobel算子')
[BW3,thresh3]=edge(I,'prewitt');
disp('prewitt自动选择的阈值为:
')
disp(thresh3)
subplot(133)
imshow(BW3)
title('prwitte算子')
BW4=filter2(fspecial('Prewitt'),I);
BW5=filter2(fspecial('Prewitt'),I);
BW6=filter2(fspecial('Sobel'),I);
figure
subplot(2,2,1);
imshow(I);
title('灰度图像');
subplot(2,2,2);
imshow(BW4);
title('Roberts的锐化结果');
subplot(2,2,3);
imshow(BW5);
title('Prewitt的锐化结果');
subplot(2,2,4);
imshow(BW6);
title('Sobel的锐化结果');
迭代后的阈值:
131
程序:
clearall;
I=imread('rice.png');
ZMAX=max(max(I));
ZMIN=min(min(I));
TK=(ZMAX+ZMIN)/2;
bcal=1;
ISIZE=size(I);
while(bcal)
ifground=0;
ibground=0;
FgroundS=0;
BgroundS=0;
fori=1:
ISIZE
(1)
forj=1:
ISIZE
(2)
tmp=I(i,j);
if(tmp>=TK)
ifground=ifground+1;
FgroundS=FgroundS+double(tmp);
else
ibground=ibground+1;
BgroundS=BgroundS+double(tmp);
end
end
end
ZO=FgroundS/ifground;
ZB=BgroundS/ibground;
TKTmp=uint8((ZO+ZB)/2);
if(TKTmp==TK)
bcal=0;
else
TK=TKTmp;
end
end
disp(strcat('迭代后的阈值:
',num2str(TK)));
newI=im2bw(I,double(TK)/255);
subplot(121)
imshow(I)
title('原始图像')
subplot(122)
imshow(newI)
title('迭代法分割得到的二值图像')
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字图像处理上机作业 数字图像 处理 上机 作业