function y=normalizacja(OR)
	[w , k ,N]=size(OR);
	clear D;
	clear HIST;
	D(w,k,1)=uint8(0);
	D=OR;
	HIST=zeros(256);
	for i=1:w
	   for j=1:k
		   for c=1:N
			   if(D(i,j,c)>=1 && D(i,j,c)<256)
				   HIST(D(i,j,c))=HIST(D(i,j,c))+1;
			   end
		   end
	   end
	end

	L2R=0;
	L1R=0;
	for i=1:256
	   if(HIST(i)>0)
		   L1R=i;
	   end
	   if(HIST(i)>0 && L2R==0)
		   L2R=i;
	   end
	end

	L2G=0;
	L1G=0;
	for i=1:256
	   if(HIST(i)>0)
		   L1G=i;
	   end
	   if(HIST(i)>0 && L2G==0)
		   L2G=i;
	   end
	end

	L2B=0;
	L1B=0;
	for i=1:256
	   if(HIST(i)>0)
		   L1B=i;
	   end
	   if(HIST(i)>0 && L2B==0)
		   L2B=i;
	   end
	end

	ob_w=OR;
	for i=1:w
	   for j=1:k
		 ob_w(i,j,1)=255*(double(OR(i,j,1))-L2R)/(L1R-L2R);
	   end
	end
	for i=1:w
	   for j=1:k
		 ob_w(i,j,2)=255*(double(OR(i,j,2))-L2G)/(L1G-L2G);
	   end
	end
	for i=1:w
	   for j=1:k
		 ob_w(i,j,3)=255*(double(OR(i,j,3))-L2B)/(L1B-L2B);
	   end
	end
	y=ob_w;
end