function [final] =uniform(b,initial)

%  uniform() 	calculates the b-bit quantization of an image
%               according to the Uniform Quantization Method

%  usage  --->  QUANT_LM = uniform(b,initial,ss)
%  inputs --->	b: 	   number of bit quantization desired
%		initial:   initial image to be quantized
%  output --->	final:	   the final quantized image
%****************************************************************


  [M,N]=size(initial);		% M,N: number of rows,columns in initial
  L=2^b;			% L:  number of intensity levels in final
  Rmin=min(min(initial));	% Rmin: minimum intensity level in initial
  Rmax=max(max(initial));	% Rmax: maximum intensity level in initial
  A=Rmax-Rmin;  		% A: dynamic range of image in int. levels
  Q=A/(L-1);		        % Q: quantization step size
 
 
  if (A == 0)
    final = zeros(M,N);
  else
    for m=1:M
     for n=1:N
      final(m,n)=round((initial(m,n)-Rmin)./Q);
     end;
    end;
  end;