原创力文档- 1、本文档共18页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
图像处理上机业题
2、用Matlab语言编写直方图规定化程序
clc;
I=imread(tire.tif);
J=histeq(I,32);
[counts,x]=imhist(J);
Q=imread(eight.tif);
figure;
imshow(Q);
figure;
imhist(Q);
M=histeq(Q,counts);
figure;
imshow(M);
figure;
imhist(M);
3、编写程序实现半邻域选择平均法进行空间域图像平滑的程序。
clear;clc;
I=imread(pout.tif);
J=imnoise(I,salt pepper,0.02);
K=filter2(fspecial(average,3),J)/255;
subplot(131),imshow(I);
subplot(132),imshow(J);
subplot(133),imshow(K);
4、针对某种平滑方法，比较对高斯噪声，椒盐噪声和乘法性噪声的去噪效果
%%
% 使用二维中值滤波对噪声进行处理
clc;
%%
%原始图像以及加上噪声后的图像
[I,map]=imread(pout.tif);
imshow(I,map);
noisy1=imnoise(I , gaussian , 0.05);%加入0.05的高斯白噪声
noisy2=imnoise(I , salt pepper , 0.05);%加入0.05的椒盐噪声
noisy3=imnoise(I , speckle , 0.05);%加入0.05的乘法性噪声
figure;
imshow(noisy1,map);
figure;
imshow(noisy2,map);
figure;
imshow(noisy3,map);
%%
%二维窗口越大图像细节丢失越多
% 3*3窗口的二维中值滤波
hood=3;
filtered1=medfilt2(noisy1,[hood hood]);
filtered2=medfilt2(noisy2,[hood hood]);
filtered3=medfilt2(noisy3,[hood hood]);
figure;
imshow(filtered1);
figure;
imshow(filtered2);
figure;
imshow(filtered3);
原图像
加噪声后
去噪声后
5、编写掩模匹配法实现锐化的程序，对同一图像比较Robison，Prewitt及Kirsch模板的锐化效果
%%
% Robison?￡°??ú?￡?￥??·¨
clc;
[I,map]=imread(pout.tif);
imshow(I,map);
[M N]=size(I);
%%
myfilt1=[ 1 2 1;
0 0 0;
-1 -2 -1];
filtered1=filter2(myfilt1,I);
%%
myfilt2=[ 2 1 0;
1 0 -1;
0 -1 -2];
filtered2=filter2(myfilt2,I);
%%
myfilt3=[ 1 0 -1;
2 0 -2;
1 0 -1];
filtered3=filter2(myfilt3,I);
%%
myfilt4=[ 0 -1 -2;
1 0 -1;
2 1 0];
filtered4=filter2(myfilt4,I);
%%
myfilt5=[-1 -2 -1;
0 0 0;
1 2 1];
filtered5=filter2(myfilt5,I);
%%
myfilt6=[-2 -1 0;
-1 0 -1;
0 1 2];
filtered6=filter2(myfilt6,I);
%%
myfilt7=[-1 0 1;
-2 0 2;
-1 0 1];
filtered7=filter2(myfilt7,I);
%%
myfilt8=[ 0 1 2;
-1 0 1;
-2 -1 0];
filtered8=filter2(myfilt8,I);
%%
for i=1:M
for j=1:N
x=[filtered1(i,j),filtered2(i,j),filtered3(i,j),filtered4(i,j),...
文档评论（0）