plot_rougher_digits.py

pca = PCA(2)
PC = pca.fit_transform(digits.data)
inversed = pca.inverse_transform(PC)
_, axes = plt.subplots(4, 10, figsize=(10, 4), 
                         subplot_kw={'xticks':[], 'yticks':[]},
                         gridspec_kw=dict(hspace=0.1, wspace=0.1))
for i, ax in enumerate(axes.flat):
    ax.imshow(inversed[i].reshape(8, 8),
              cmap='binary', interpolation='nearest', clim=(0, 16))