%*************** EECS_556  IMAGE PROCESSING *******************
%***************    Computer Assignment #1  *******************
%**********                                     ***************
%**********     James J. Huth   Joe Mancuso     ***************
%**********               2/6/96                ***************
%**************************************************************
% Explore the application of Lloyd-Max and uniform quantizers,
% 2D:DCT, 2D:DFT, etc.  to an image.  Investigate sampling rates
% and aliasing, SNR.  


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

%***** Load Image
  load house;


%**************************************************************
%*********** Part c *******************************************
%*** Find the 2D DCT of Image *********************************

  w=0:1/255:1;

  house_dct=dct2(house);
   save house_dct.mat house_dct;
   dcterm1=house_dct(1,1);
   house_dct(1,1)=house_dct(2,2);
%   DC1= (log(dcterm1));



   house_dct_un8=uniform(8,log(1000.*abs(house_dct)));
  
  figure;
   mesh(w,w,round(house_dct_un8));view(0,90);colormap(gray);
  title(['2D DCT of Original Image - DC Term Excluded']);