PS 滤镜——旋转模糊

这里给出灰度图像的模糊算法，彩色图像只要分别对三个通道做模糊即可。

%%  spin blur

% 旋转模糊
clc;
clear all;
close all;

I=imread('4.jpg');
I=double(I);

% % % I_new=I;
% % % for kk=1:3
% % %     I_new(:,:,kk)=Spin_blur_Fun(I(:,:,kk), 30, 30);
% % % end

% % % imshow(I_new/255)

Image=I;
Image=0.2989 * I(:,:,1) + 0.5870 * I(:,:,2) + 0.1140 * I(:,:,3); 

[row, col]=size(Image);
Image_new=Image;
Center_X=(col+1)/2;
Center_Y=(row+1)/2;
validPoint=1;
angle=5;
radian=angle*pi/180;
radian2=radian*radian;
Num=30;
Num2=Num*Num;
for i=1:row
    for j=1:col
        validPoint=1;
        x0=j-Center_X;
        y0=Center_Y-i;
        x1=x0;
        y1=y0;
        Sum_Pixel=Image(i,j);
        for k=1:Num
            x0=x1;
            y0=y1;
            %%% 逆时针
            % x1=x0-radian*y0/Num-radian2*x0/Num2;
            % y1=y0+radian*x0/Num-radian2*y0/Num2;
            %%% 顺时针
            x1=x0+radian*y0/Num-radian2*x0/Num2;
            y1=y0-radian*x0/Num-radian2*y0/Num2;
            
            x=floor(x1+Center_X);
            y=floor(Center_Y-y1);
            
            if(x>1 && x<col && y>1 && y<row)
                validPoint=validPoint+1;
                Sum_Pixel=Sum_Pixel+Image(y,x);
            end   
        end
        Image_new(i,j)=Sum_Pixel/validPoint;
    end
end

 imshow(Image_new/255);      

原图 

效果图

效果图