HR图像添加运动焦散模糊【Matlab】

直接上code和效果图：

I = imread('C:\Users\jizef\Desktop\flower.jpg');
figure,subplot(2,2,1);imshow(I);title('原图像');
%T = rgb2gray(M);figure;imshow(I);
PSF = fspecial('motion',20,15);
J = imfilter(I,PSF,'conv','circular');
subplot(2,2,2);imshow(J);title('模糊运动的图像');

%'计算中......'
% I=imread('C:\Users\jizef\Desktop\flower.jpg');
% figure,subplot(1,3,1);imshow(I);title('原图像');
r=7;%散焦半径r
PSF=fspecial('disk',r);   %得到点扩散函数
I1=imfilter(I,PSF,'symmetric','conv');  %实现散焦模糊
subplot(2,2,3);imshow(I1);title('散焦模糊的图像');

I2=imfilter(J,PSF,'symmetric','conv');  %实现运动散焦模糊
subplot(2,2,4);imshow(I2);title('运动散焦模糊的图像');

运行结果：