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

%*********** Part b

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

%***** Load Image
  load house;

%**************************************************************
% Reduce sampling rate...

[M,N] = size(house);

for m = 1:N
  for n = 1:N
    if ((rem((m-1),2) == 0) & (rem((n-1),2) == 0))
      house_r2(((m-1)/2)+1,((n-1)/2)+1) = house(m,n);
    end;
    if ((rem((m-1),3) == 0) & (rem((n-1),3) == 0))
      house_r3(((m-1)/3)+1,((n-1)/3)+1) = house(m,n);
    end;
    if ((rem((m-1),4) == 0) & (rem((n-1),4) == 0))
      house_r4(((m-1)/4)+1,((n-1)/4)+1) = house(m,n);
    end;
    if ((rem((m-1),5) == 0) & (rem((n-1),5) == 0))
      house_r5(((m-1)/5)+1,((n-1)/5)+1) = house(m,n);
    end;
    if ((rem((m-1),6) == 0) & (rem((n-1),6) == 0))
      house_r6(((m-1)/6)+1,((n-1)/6)+1) = house(m,n);
    end;
    if ((rem((m-1),7) == 0) & (rem((n-1),7) == 0))
      house_r7(((m-1)/7)+1,((n-1)/7)+1) = house(m,n);
    end;
    if ((rem((m-1),8) == 0) & (rem((n-1),8) == 0))
      house_r8(((m-1)/8)+1,((n-1)/8)+1) = house(m,n);
    end;
  end;
end;

figure;
image(log((abs(fftshift(fft2(house))).^1)))
title('House spectrum');
print -dgif8 house_spec;

figure;
image(log((abs(fftshift(fft2(house_r2))).^1)))
title('House resampled at every 2nd point spectrum');
print -dgif8 house_r2;

figure;
image(log((abs(fftshift(fft2(house_r3))).^1)))
title('House resampled at every 3rd point spectrum');
print -dgif8 house_r3;

figure;
image(log((abs(fftshift(fft2(house_r4))).^1)))
title('House resampled at every 4th point spectrum');
print -dgif8 house_r4;

figure;
image(log((abs(fftshift(fft2(house_r5))).^1)))
title('House resampled at every 5th point spectrum');
print -dgif8 house_r5;

figure;
image(log((abs(fftshift(fft2(house_r6))).^1)))
title('House resampled at every 6th point spectrum');
print -dgif8 house_r6;

figure;
image(log((abs(fftshift(fft2(house_r7))).^1)))
title('House resampled at every 7th point spectrum');
print -dgif8 house_r7;

figure;
image(log((abs(fftshift(fft2(house_r8))).^1)))
title('House resampled at every 8th point spectrum');
print -dgif8 house_r8;