%%% linear diffusion demo %%
%%% By Guy Gilboa

I=imread('haifa1.bmp'); % load image (true color)
%I=double(I(20:120,10:105,:)); % cut a piece, convert to double
IR=double(I(:,:,1));
IG=double(I(:,:,2));
IB=double(I(:,:,3));

close all
figure(1); imshow(uint8(I)); 
% run linear diffusion
J=I;
for iter=[1,10],
	J(:,:,1)=diffusion(IR,'lin',iter); 
	J(:,:,2)=diffusion(IG,'lin',iter); 
	J(:,:,3)=diffusion(IB,'lin',iter); 
   figure; imshow(uint8(J)); 
end % for iter