import matplotlib.pyplot as plt
%matplotlib inline
I0 = plt.imread('SiebelCtr_01.jpg')
print("Dimension of I0: ", I0.shape)
plt.figure(figsize=(10,10))
plt.imshow(I0, interpolation='none')
I1 = I0[::2, ::2, :]
print("Dimension of I1: ", I1.shape)
plt.figure(figsize=(10,10))
plt.imshow(I1, interpolation='none')
I2 = I1[::2, ::2, :]
print("Dimension of I2: ", I2.shape)
plt.figure(figsize=(10,10))
plt.imshow(I2, interpolation='none')
I3 = I2[::2, ::2, :]
print("Dimension of I3: ", I3.shape)
plt.figure(figsize=(10,10))
plt.imshow(I3, interpolation='none')
I4 = I3[::2, ::2, :]
print("Dimension of I4: ", I4.shape)
plt.figure(figsize=(10,10))
plt.imshow(I4, interpolation='none')
I5 = I4[::2, ::2, :]
print("Dimension of I5: ", I5.shape)
plt.figure(figsize=(10,10))
plt.imshow(I5, interpolation='none')
I6 = I5[::2, ::2, :]
print("Dimension of I6: ", I6.shape)
plt.figure(figsize=(10,10))
plt.imshow(I6, interpolation='none')
I7 = I6[::2, ::2, :]
print("Dimension of I7: ", I7.shape)
plt.figure(figsize=(10,10))
plt.imshow(I7, interpolation='none')
J = plt.imread('somerock.jpg')
print("J is ", J.shape)
Jold = J.copy()
for j in range(9):
Jcoarse = Jold[::2, ::2, :]
Jold = Jcoarse.copy()
print("J is ", Jcoarse.shape)
plt.figure(figsize=(10,10))
plt.imshow(Jcoarse, interpolation='none')
plt.figure(figsize=(10,10))
plt.imshow(J, interpolation='none')