%*************** EECS_556  IMAGE PROCESSING *******************
%***************    Computer Assignment #1  *******************
%***************                            *******************
%***************       James J. Huth        *******************
%***************           2/6/96           *******************
%**************************************************************

%*********** Part d


%***** Initialize 
  close all;
   clear;

%***** Set resolution for original transform quantization
  resol = 12;

%***** Load Image
   load house_dct;
   load house_dft;

%**************************************************************
% Part(d1): Requantize the image to 8 bits using both a Lloyd-Max Quantizer

  house_dct_r = real(house_dct);
  house_dct_i = imag(house_dct);

  house_dft_r = real(house_dft);
  house_dft_i = imag(house_dft);

  house_dct_r = uniform(resol,house_dct_r);
  house_dct_i = uniform(resol,house_dct_i);
  house_dft_r = uniform(resol,house_dft_r);
  house_dft_i = uniform(resol,house_dft_i);

  house_dct_2bl = (lloyd_max(2,house_dct_r,resol) + j*lloyd_max(2,house_dct_i,resol));
  save house_dct_2bl house_dct_2bl
  house_dft_2bl = (lloyd_max(2,house_dft_r,resol) + j*lloyd_max(2,house_dft_i,resol));
  save house_dft_2bl house_dft_2bl

  house_dct_4bl = (lloyd_max(4,house_dct_r,resol) + j*lloyd_max(4,house_dct_i,resol));
  save house_dct_4bl house_dct_4bl
  house_dft_4bl = (lloyd_max(4,house_dft_r,resol) + j*lloyd_max(4,house_dft_i,resol));
  save house_dft_4bl house_dft_4bl

  house_dct_6bl = (lloyd_max(6,house_dct_r,resol) + j*lloyd_max(6,house_dct_i,resol));
  save house_dct_6bl house_dct_6bl
  house_dft_6bl = (lloyd_max(6,house_dft_r,resol) + j*lloyd_max(6,house_dft_i,resol));
  save house_dft_6bl house_dft_6bl

  house_dct_8bl = (lloyd_max(8,house_dct_r,resol) + j*lloyd_max(8,house_dct_i,resol));
  save house_dct_8bl house_dct_8bl
  house_dft_8bl = (lloyd_max(8,house_dft_r,resol) + j*lloyd_max(8,house_dft_i,resol));
  save house_dft_8bl house_dft_8bl

  figure;
  image(abs(house_dct_2bl));
  title('2-Bit L-M Quantization of DCT of Original Image');
  print -dgif8 partd/house_dct_2bl
  figure;
  image(abs(house_dft_2bl));
  title('2-Bit L-M Quantization of DFT of Original Image');
  print -dgif8 partd/house_dft_2bl

  figure;
  image(abs(house_dct_4bl));
  title('4-Bit L-M Quantization of DCT of Original Image');
  print -dgif8 partd/house_dct_4bl
  figure;
  image(abs(house_dft_4bl));
  title('4-Bit L-M Quantization of DFT of Original Image');
  print -dgif8 partd/house_dft_4bl

  figure;
  image(abs(house_dct_6bl));
  title('6-Bit L-M Quantization of DCT of Original Image');
  print -dgif8 partd/house_dct_6bl
  figure;
  image(abs(house_dft_6bl));
  title('6-Bit L-M Quantization of DFT of Original Image');
  print -dgif8 partd/house_dft_6bl

  figure;
  image(abs(house_dct_8bl));
  title('8-Bit L-M Quantization of DCT of Original Image');
  print -dgif8 partd/house_dct_8bl
  figure;
  image(abs(house_dft_8bl));
  title('8-Bit L-M Quantization of DFT of Original Image');
  print -dgif8 partd/house_dft_8bl