cmap = gen_cmap([plt.cm.terrain(i) for i in np.int32(255*np.linspace(0.0, 0.97, 256))])