PS 滤镜— —挤压效果

    clc;
    clear all;
    close all;
    addpath('E:\PhotoShop Algortihm\Image Processing\PS Algorithm');
    I=imread('4.jpg');
    I=double(I);
    Image=I/255;
    angle = pi/4;
    centreX = 0.5;
    centreY = 0.5;
    radius=200;
    amount=0.75;
    [height, width, depth]=size(Image);
    Image_new=Image;
    icentreX = width * centreX;
    icentreY = height * centreY;
    if ( radius == 0 )
        radius=min(icentreX, icentreY);
    end
    radius2=radius*radius;
    for ii=1:height
        for jj=1:width
            dx = jj-icentreX;
            dy = ii-icentreY;
            distance = dx*dx + dy*dy;
            if (distance>radius2 || distance==0)
                x=jj;
                y=ii;
            else
                 d = sqrt( distance / radius2 );
                 t = sin( pi*0.5 * d ).^(-amount);
                dx =dx* t;
                dy =dy* t;
                e = 1 - d;
                a = angle * e * e;
                s = sin( a );
                c = cos( a );
                x = icentreX + c*dx - s*dy;
                y = icentreY + s*dx + c*dy;
            end
            if (x<=1)     x=1;  end
            if (x>=width)   x=width-1; end;
            if (y>=height)  y=height-1; end;
            if (y<1)  y=1;     end;
    % %         if (x<=1)     continue;  end
    % %         if (x>=width)   continue; end;
    % %         if (y>=height)  continue; end;
    % %         if (y<1)  continue;     end;
            x1=floor(x);
            y1=floor(y);
            p=x-x1;
            q=y-y1;
            Image_new(ii,jj,:)=(1-p)*(1-q)*Image(y1,x1,:)+p*(1-q)*Image(y1,x1+1,:)...
                +q*(1-p)*Image(y1+1,x1,:)+p*q*Image(y1+1,x1+1,:); 
        end
    end
    imshow(Image_new)
    imwrite(Image_new, 'out.jpg');

参考来源：http://www.jhlabs.com/index.html

原图：

效果图：