matlab之PS算法小程序（差值，溶解， 排除）

原理见http://blog.csdn.net/matrix_space/article/details/22488467

[filename,pathname] = uigetfile('*.jpg;*.bmp','选择图片','E:\pictures\For_Project\Matlab');
imgaepath = strcat(pathname,filename);
I1 = imread(imgaepath);

[filename,pathname] = uigetfile('*.jpg;*.bmp','选择图片','E:\pictures\For_Project\Matlab');
imgaepath = strcat(pathname,filename);
I2 = imread(imgaepath);

I1 = (imresize(I1,[600,800]));
I2 = (imresize(I2,[600,800]));

C = imabsdiff(I1,I2);

figure(1),set(gcf,'Name','差值');
subplot(1,3,1),imshow(I1),title('图片1');
subplot(1,3,2),imshow(I2),title('图片2');
subplot(1,3,3),imshow(C),title('图片3');

[filename,pathname] = uigetfile('*.jpg;*.bmp','选择图片','E:\pictures\For_Project\Matlab');
imgaepath = strcat(pathname,filename);
I1 = imread(imgaepath);

[filename,pathname] = uigetfile('*.jpg;*.bmp','选择图片','E:\pictures\For_Project\Matlab');
imgaepath = strcat(pathname,filename);
I2 = imread(imgaepath);

I1 = double(imresize(I1,[600,800]))/255;
I2 = double(imresize(I2,[600,800]))/255;

for i = 1:600
    for j = 1:800
        C(i,j) = I1(i,j) + I2(i,j) - 2*I1(i,j)*I2(i,j);
    end
end

figure(1),set(gcf,'Name','排除');
subplot(1,3,1),imshow(I1),title('图片1');
subplot(1,3,2),imshow(I2),title('图片2');
subplot(1,3,3),imshow(C),title('图片3');
figure(2),imshow(C);