%%%%%  motion blur
clc;
clear all;
close all;
Image=imread('4.jpg');
Image=double(Image);
theta=pi/4;
len=20;
row=floor(len*sin(theta))+1;
col=floor(len*cos(theta))+1;
motion_f1(1:row,1:col)=0;
K=tan(theta);
center_x=col/2;
center_y=row/2;

for i=1:row
    for j=1:col
        x=j-center_x;
        y=center_y-i;
        Dis=abs(K*x-y)/sqrt(K*K+1);
        motion_f1(i,j)=max(1-Dis,0);
    end
end

motion_f1=motion_f1/sum(motion_f1(:));
Image1=imfilter(Image,motion_f1,'replicate','conv');
imshow(Image1/255);

motion_filter=fspecial('motion',20,0);
Image2=imfilter(Image,motion_filter,'replicate','conv');
figure,

imshow(Image2/255);


原图 


效果图 


效果图 


posted on 2014-10-08 13:51  未雨愁眸  阅读(212)  评论(0编辑  收藏  举报