162--slightly-weather.zip_evil

I = imread('F:\work\rain\aa.jpg'); %figure;imshow(I);title('Original Image'); %text(size(I,2),size(I,1)+15, ... % 'Image courtesy of Massachusetts Institute of Technology', ... % 'FontSize',7,'HorizontalAlignment','right'); I1 = imread('F:\work\rain\211.jpg'); %figure;imshow(I1);title('midle Image'); I2 = imread('F:\work\rain\212.jpg'); % I3 = imread('F:\work\rain\25.jpg'); %figure;imshow(I2);title('third Image'); Igary=rgb2gray(I); % figure;imshow(Igary);title('Original Image'); Igary1=rgb2gray(I1); figure;imshow(Igary1);title('midle Image'); Igary2=rgb2gray(I2); % Igary4=rgb2gray(I3); % figure;imshow(Igary2);title('third Image'); Igary3=Igary1; for k=1:240 for l=1:320 % if (((Igary1(k,l)-Igary(k,l))>10)&&((Igary1(k,l)-Igary2(k,l))>10)) a=(Igary1(k,l)-Igary(k,l)); b=(Igary1(k,l)-Igary2(k,l)); c=(0.039*double(Igary1(k,l))-double(a))/double(Igary(k,l)); % if((a==b)&&(a>1)) if((abs(a-b)<3)&&a>3) if((c<0.039)) Igary3(k,l)=255; % Igary3(k,l)=Igary2(k,l); end else Igary3(k,l)=0; % Igary3(k,l)=Igary3(k,l); end end end figure;imshow(Igary3);title('procesing Image'); Ip=Igary1-Igary3; % Ip = imsubtract(Igary3,Igary1); figure;imshow(Ip);title('subimage Image');