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

I=imread('haifa1.bmp'); % load image (true color)
IR=double(I(:,:,1));
IG=double(I(:,:,2));
IB=double(I(:,:,3));

%parameters
theta = pi/1000;  % phase angle of diffusion coefficient

close all
figure(1); imshow(uint8(I)); 
figure; imshow(uint8(128+imag(double(I))));  % zero imaginary value
% run linear diffusion
J=double(I);
j=sqrt(-1);
C=exp(j*theta)
for iter=[1,10],
	J(:,:,1)=diffusion(IR,'lin',iter,C); 
	J(:,:,2)=diffusion(IG,'lin',iter,C); 
	J(:,:,3)=diffusion(IB,'lin',iter,C); 
	figure; imshow(uint8(real(J)));  % real
   figure; imshow(uint8(128+imag(J)/theta)); % imaginary
end % for iter

% imwrite(uint8(I),'cmplin_haifa_0.jpg')