Matlab:盲反卷积还原失焦图像

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%       盲反卷积还原失焦图像             %
%                                       %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

I = imread('G:\BaiduNetdiskDownload\5.jpg');
I=rgb2gray(I)
PSF=fspecial('gaussian',7,30);
Blurred=imfilter(I,PSF,'symmetric','conv');  
figure;imshow(Blurred);title('失焦图像');   
 
 
UNDERPSF = ones(size(PSF)-4);
[J1 P1] = deconvblind(Blurred,UNDERPSF);
figure;imshow(J1);title('Deblurring with Undersized PSF'); 
   
OVERPSF = padarray(UNDERPSF,[4  4],'replicate','both');
[J2 P2] = deconvblind(Blurred,OVERPSF);
figure;imshow(J2);title('Deblurring with Oversized PSF');     
INITPSF = padarray(UNDERPSF,[2 2],'replicate','both');
[J3 P3] = deconvblind(Blurred,INITPSF);
figure;imshow(J3);title('Deblurring with INITPSF');  
   
figure;
subplot(221);imshow(P1,[],'InitialMagnification','fit');
title('重建 Undersized PSF');
subplot(222);imshow(P2,[],'InitialMagnification','fit');
title('重建 Oversized PSF');
subplot(223);imshow(P3,[],'InitialMagnification','fit');
title('重建  init PSF');  
 
   
WEIGHT = edge(I,'sobel',.3);  
se = strel('disk',2);
WEIGHT = 1-double(imdilate(WEIGHT,se));  
 
WEIGHT([1:3 end-[0:2]],:) = 0;
WEIGHT(:,[1:3 end-[0:2]]) = 0;
figure;imshow(WEIGHT);title('Weight array'); 
[J P] =                                  deconvblind(Blurred,INITPSF,30,[],WEIGHT);
figure;imshow(J);title('复原图像');  
    
P1 = 2;
P2 = 2;
FUN = @(PSF) padarray(PSF(P1+1:end-P1,P2+1:end-P2),[P1 P2]);  

[JF PF] = deconvblind(Blurred,OVERPSF,30,[],WEIGHT,FUN);
figure;imshow(JF);title('复原图像');   
displayEndOfDemoMessage(mfilename) 

阅读更多

扫码向博主提问

一个处女座的程序猿

非学,无以致疑;非问,无以广识
  • 擅长领域:
  • ML
  • DL
  • BlockChain
  • CV
  • Programmin
去开通我的Chat快问
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/qq_41185868/article/details/80345551
上一篇Matlab:绘制简单能量的接收机工作特性曲线(Energy_detection_simulation_ok)
下一篇MAT之GA:遗传算法(GA)解决M-TSP多旅行商问题
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭