grid_example.py

import matplotlib.pyplot as plt
import matplotlib.gridspec as gridspec
import numpy as np
import cv2

def jemma():
    img=cv2.imread('jemma.png')
    #img = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)
    img = img[100:300,100:300]
    return img


def squiggle_xy(a, b, c, d):
    i = np.arange(0.0, 2*np.pi, 0.05)
    return np.sin(i*a)*np.cos(i*b), np.sin(i*c)*np.cos(i*d)


fig = plt.figure(figsize=(8, 8))

# gridspec inside gridspec
outer_grid = gridspec.GridSpec(8, 8, wspace=0.0, hspace=0.0)

for i in range(64):
    ax=plt.Subplot(fig, outer_grid[i])
    a = i // 4 + 1
    b = i % 4 + 1
    c = b
    d = a
    ax.imshow(jemma(), cmap='Greys')
    #ax.plot(*squiggle_xy(a, b, c, d))
    ax.set_xticks([])
    ax.set_yticks([])
    fig.add_subplot(ax)

plt.show()
"""
for i in range(16):
    inner_grid = gridspec.GridSpecFromSubplotSpec(3, 3,
            subplot_spec=outer_grid[i], wspace=0.0, hspace=0.0)
    a = i // 4 + 1
    b = i % 4 + 1
    for j, (c, d) in enumerate(product(range(1, 4), repeat=2)):
        ax = plt.Subplot(fig, inner_grid[j])
        ax.plot(*squiggle_xy(a, b, c, d))
        ax.set_xticks([])
        ax.set_yticks([])
        fig.add_subplot(ax)

all_axes = fig.get_axes()

#show only the outside spines
for ax in all_axes:
    for sp in ax.spines.values():
        sp.set_visible(False)
    if ax.is_first_row():
        ax.spines['top'].set_visible(True)
    if ax.is_last_row():
        ax.spines['bottom'].set_visible(True)
    if ax.is_first_col():
        ax.spines['left'].set_visible(True)
    if ax.is_last_col():
        ax.spines['right'].set_visible(True)
"""

plt.show()