from matplotlib import pyplot as plt
from ligo.skymap.plot import bluemarble, reproject_interp_rgb

obstime = '2017-08-17 12:41:04'
ax = plt.axes(projection='geo degrees aitoff', obstime=obstime)
ax.imshow(reproject_interp_rgb(bluemarble(obstime), ax.header))