cmap = gen_cmap([colors.to_rgba('lemonchiffon')]\
                 + [plt.cm.YlOrRd(i) for i in np.int32(255*np.linspace(0.15, 0.8, 6))])