- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
图像复原程序
adaptive1.m
% adaptive filter
gimage=double(imread(fig513.jpg));
subplot(121);
imagesc(gimage);
colormap(gray);
fimage=gimage;
[m,n]=size(gimage);
for i=4:m-3
for j=4:n-3
plate=gimage(i-3:i+3,j-3:j+1);
ml=sum(sum(plate))/49;
el=sum(sum((plate-ml).*(plate-ml)))/49;
fimage(i,j)=gimage(i,j)-1000/el*(gimage(i,j)-ml);
end
end
subplot(122);
imagesc(fimage);
colormap(gray);
aver.m
% average filter
% filter 1
o_i=double(imread(fig57a.jpg));
subplot(221);
imagesc(o_i);
colormap(gray);
title(原始图像)
o_i_n=double(imread(fig57b.jpg));
subplot(222);
imagesc(o_i_n);
colormap(gray);
title(含噪图像)
[m,n]=size(o_i);
plate=zeros(3,3);
f=o_i_n;
for i=2:m-1
for j=2:n-1
plate=o_i_n(i-1:i+1,j-1:j+1);
f(i,j)=sum(sum(plate))/9;
end
end
subplot(223);
imagesc(f);
colormap(gray);
title(算术平均)
%filter 2
[m,n]=size(o_i);
plate=zeros(3,3);
f=o_i_n;
for i=2:m-1
for j=2:n-1
plate=o_i_n(i-1:i+1,j-1:j+1); f(i,j)=(plate(1,1)*plate(1,2)*plate(1,3)*plate(2,1)*plate(2,2)*plate(2,3)*plate(3,1)*plate(3,2)*plate(3,3))^(1/9);
end
end
subplot(224);
imagesc(f);
colormap(gray);
title(几何平均)
f43.m
woye
%image and noise
%their distogram
imageyin=double(imread(lena512.bmp));
[m,n]=size(imageyin);
imageline=imageyin(1,1:n);
for i=2:m
imageline=[imageline imageyin(i,1:n)];
end
hist(imageline,256);
f43abc.m
%image and noise
%their distogram
clear all
clc
imageyin=double(imread(fig53.jpg));
[m,n]=size(imageyin);
imageline=imageyin(1,1:n);
for i=2:m
imageline=[imageline imageyin(i,1:n)];
end
subplot(241);
imagesc(imageyin);
colormap(gray);
subplot(245);
hist(imageline,-1:1:255);
%高斯噪声
k=20;
imagegauss=k*randn(m,n);
imagey_g=imageyin+imagegauss;
imagey_g=(imagey_g-min(min(imagey_g)))/(max(max(imagey_g))-min(min(imagey_g)))*255;
subplot(242);
imagesc(imagey_g);
colormap(gray);
subplot(246);
imageline=imagey_g(1,1:n);
for i=2:m
imageline=[imageline imagey_g(i,1:n)];
end
hist(imageline,1:1:256);
% 随机噪声
k=40;
imagenorm=k*rand(m,n);
imag
文档评论(0)