from pylab import *
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.cbook as cbook
import random
import time
from scipy.misc import imread
from scipy.misc import imresize
import matplotlib.image as mpimg
import os
from scipy.ndimage import filters

def downsample(im):
    N = im.shape[0]
    M = im.shape[1]
    m = meshgrid(range(0, N, 2), range(0, M, 2))
    return im[m]

gray()
os.chdir('c:/Users/Guerzhoy/Desktop/csc320/fourier/')
im = imread("grumpy.jpg")[:,:,0].astype(float)



close('all')
im = filters.gaussian_filter(im, 2)
im = downsample(im)

im = filters.gaussian_filter(im, 2)
im = downsample(im)

im = filters.gaussian_filter(im, 2)

imshow(im)
show()

im = imread("grumpy.jpg")[:,:,1].astype(float)
im = filters.gaussian_filter(im, sqrt(4+16+64))
im = downsample(im)
im = downsample(im)
figure(2)
imshow(im)
show()