没有合适的资源?快使用搜索试试~ 我知道了~
matlab去模糊的几种方法.doc
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 32 浏览量
2022-07-05
20:01:06
上传
评论
收藏 36KB DOC 举报
温馨提示
试读
3页
matlab去模糊的几种方法
资源推荐
资源详情
资源评论
I = imread('peppers.png');%¶ÁȡͼÏñ
I = I(10+[1:256],222+[1:256],:);%Ñ¡¶¨Í¼ÏñµÄ·¶Î§
figure(1); imshow(I);%ÏÔʾÔͼÏñ
title('Original Image');
LEN = 31; THETA = 11;
PSF = fspecial('motion',LEN,THETA);%²úÉúÔ˶¯Ä£ºýº¯Êý
blurred = imfilter(I,PSF,'circular','conv');%²úÉúÔ˶¯Ä£ºýͼÏñ
figure(2); subplot(221),
imshow(blurred);%ÏÔʾÔ˶¯Ä£ºýͼÏñ
title('Blurred');
wnr1 = deconvwnr(blurred,PSF);%ʹÓÃάÄÉÂ˲¨Æ÷½øÐÐÂ˲¨
subplot(222); imshow(wnr1);%ÏÔʾάÄÉÂ˲¨Æ÷ºóÂ˲¨µÄͼÏñ
title('Restored, True PSF');
wnr2 = deconvwnr(blurred,...
fspecial('motion',2*LEN,THETA));%ʹÓÃÁ½±¶³¤¶ÈµÄµãÀ©É¢º¯Êý½øÐÐάÄÉÂ˲¨
subplot(223); imshow(wnr2); %ÏÔʾ»Ö¸´µÄͼÏñ
title('Restored, "Long" PSF');
wnr3 = deconvwnr(blurred,...
fspecial('motion',LEN,2*THETA));%ʹÓÃÁ½±¶½Ç¶ÈµÄµãÀ©É¢º¯Êý½øÐÐάÄÉÂ˲¨
subplot(224); imshow(wnr3); %ÏÔʾ»Ö¸´µÄͼÏñ
title('Restored, Steep');
noise = 0.1*randn(size(I));%Éú³ÉËæ»úÔëÉù
blurredNoisy = imadd(blurred,im2uint8(noise));%Ä£ºýͼÏñ¼ÓÈëËæ»úÔëÉù
figure(3); subplot(221)
imshow(blurredNoisy);%ÏÔʾ¼ÓÈëÔëÉùµÄÄ£ºýͼÏñ
title('Blurred & Noisy');
wnr4 = deconvwnr(blurredNoisy,PSF);%ʹÓÃάÄÉÂ˲¨Æ÷½øÐÐÂ˲¨
subplot(222),
imshow(wnr4);%ÏÔʾ»Ö¸´ºóµÄͼÏñ
title('Inverse Filtering of Noisy Data');
NSR = sum(noise(:).^2)/sum(im2double(I(:)).^2);%¼ÆËãÔëÉùÐźŹ¦ÂʱÈ
wnr5 = deconvwnr(blurredNoisy,PSF,NSR);%¼ÓÈëNSR²ÎÊýʹÓÃάÄÉÂ˲¨Æ÷Â˲¨
subplot(223); imshow(wnr5);%ÏÔʾ»Ö¸´ºóµÄͼÏñ
title('Restored with NSR');
wnr6 =deconvwnr(blurredNoisy,PSF,NSR/2);%¸Ä±äNSR²ÎÊýʹÓÃάÄÉÂ˲¨Æ÷Â˲¨
subplot(224),
imshow(wnr6);%ÏÔʾ»Ö¸´ºóµÄͼÏñ
title('Restored with NSR/2');
NP = abs(fftn(noise)).^2;%¼ÆËãÔëÉù¹¦ÂÊÆ×
NPOW = sum(NP(:))/prod(size(noise)); % ÔëÉùÄÜÁ¿
NCORR = fftshift(real(ifftn(NP))); % ÔëÉùµÄ×ÔÏà¹Øº¯Êý
IP = abs(fftn(im2double(I))).^2;%¼ÆËãͼÏñ¹¦ÂÊÆ×
IPOW = sum(IP(:))/prod(size(I)); % ͼÏñÄÜÁ¿
ICORR = fftshift(real(ifftn(IP))); % ͼÏñ×ÔÏà¹Øº¯Êý
wnr7 = deconvwnr(blurredNoisy,PSF,NCORR,ICORR);%Ôö¼Ó²ÎÊý½øÐÐάÄÉÂ˲¨
figure(4); subplot(121)
imshow(wnr7);%ÏÔʾ»Ö¸´ºóµÄͼÏñ
title('Restored with ACF');
ICORR1 = ICORR(:,ceil(size(I,1)/2));
wnr8 = deconvwnr(blurredNoisy,PSF,NPOW,...
ICORR1);%ʹÓÃÔëÉùÄÜÁ¿ºÍͼÏñһά×ÔÏà¹Øº¯Êý½øÐÐάÄÉÂ˲¨
subplot(122); imshow(wnr8);%ÏÔʾ»Ö¸´ºóµÄͼÏñ
title('Restored with NP & 1D-ACF');
PSF2 = fspecial('gaussian',11,5);%¸ß˹ģºýº¯Êý
Blurred = imfilter(I,PSF2,'conv');%Ä£ºýºóµÄͼÏñ
figure(5); subplot(121)
imshow(Blurred);%ÏÔʾģºýºóµÄͼÏñ
资源评论
老帽爬新坡
- 粉丝: 82
- 资源: 2万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功