clc;
    clear all;
    close all;

    addpath('E:\PhotoShop Algortihm\Image Processing\PS Algorithm');

    I=imread('4.jpg');
    I=double(I);
    Image=I/255;

    xOffset=100;
    yOffset=200;
    warp=1;

    [height, width, depth]=size(Image);

    if (warp)
        while(xOffset<0)
            xOffset=xOffset+width;
        end
        while(yOffset<0)
            yOffset=yOffset+height;
        end

        xOffset=mod(xOffset, width);
        yOffset=mod(yOffset, height);
    end

    Image_new=Image;

    for i=1:height

        for j=1:width

            if (warp)
                x=mod(j+width-xOffset, width);
                y=mod(i+height-yOffset, height);
            else
                x=j-xOffset;
                y=i-yOffset;
            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(i,j,:)=(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

原图:
这里写图片描述

效果图:

这里写图片描述

posted on 2015-10-10 10:16  未雨愁眸  阅读(456)  评论(0编辑  收藏  举报