用MATLAB做简单的运动模糊

clear all
I=imread('cameraman.tif');
subplot(121);imshow(I);
xlabel('(a)原始图像');%x轴标签
len=35;%设置运动位移为35个像素
theta=45;%设置运动角度为45度
psf=fspecial('motion',len,theta);%建立二维运动仿真滤波器psf
mf=imfilter(I,psf,'circular','conv');%用psf产生退化图像
subplot(122);imshow(mf);
xlabel('(b)运动模糊图像');
imwrite(mf,'cameraman1.tif');%将运动得到的模糊图像进行保存