clear all;
%% 读取图片
img=imread('ground.jpg');
%img=rgb2gray(img); %转为灰度图
img=im2double(img);
%% 去模糊比较
subplot(2,2,1),imshow(img),title('initial');
len=8; %模糊距离 这个需要根据模糊的图像估算
theta=9; %模糊角度
PSF=fspecial('motion',len,theta);
sVar=var(img(:));
wnr1=deconvwnr(img,PSF,0);
subplot(2,2,2),imshow(wnr1),title('without using NSR')
wnr2=deconvwnr(img,PSF,0.0001/sVar);
subplot(2,2,3),imshow(wnr2),title('with using NSR=0.0001');
wnr3=deconvwnr(img,PSF,0.001/sVar);
subplot(2,2,4),imshow(wnr3),title('with using NSR=0.001');
%figure(2),imshow(wnr2)
%h=fspecial('average');
%f=filter2(h,wnr2);
%figure(3),imshow(f)
%% 频谱 估计运动模糊方向(光斑的角度与真实图像运动模糊距离垂直)
I2=imfilter(img,[0 0 0 0 0 0 0 0 5/25 4/25 4/25 3/25 3/25 2/25 2/25 1/25 1/25]);
I_hsv=rgb2hsv(I2);
I=I_hsv(:,:,3)*255;
img_fft=fftshift(fft2(I));
N=abs(img_fft);
P=(N-min(min(N)))/(max(max(N))-min(min(N)))*225;
figure(2);imshow(P);
%% 估计运动模糊长度(两个最小点之间的距离为近似运动模糊距离)
h=fspecial('sobel');
img_double=double(I);
J=conv2(img_double,h,'same');
IP=abs(fft2(J));
S=fftshift(real(ifft2(IP)));
figure(3);plot(S)
sport_M.m.zip
版权申诉
171 浏览量
2024-04-16
00:10:48
上传
评论
收藏 1KB ZIP 举报
手把手教你学AI
- 粉丝: 7972
- 资源: 4768
最新资源
- Python 手写实现 iD3 决策树算法-根据信息增益公式.zip
- 411675952289057车联助手-小窗版(三星)3.5.1.apk
- 三种快速排序方法合并在一个文件中以便直接运行的Python代码示例
- 937712277954201实习5.word
- 2程序语言基础知识pdf1_1716337722703.jpeg
- 简单的Python示例,演示了如何使用TCP/IP协议进行基本的客户端和服务器通信
- 考试.sql
- keil2 + proteus + 8051.exe
- 1961ee27df03bd4595d28e24b00dde4e_744c805f7e4fb4d40fa3f695bfbab035_8(1).c
- mediapipe-0.9.0.1-cp37-cp37m-win-amd64.whl.zip
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈