import matplotlib.pyplot as plt
import astropy.units as u

import sunpy.map
from sunpy.data.sample import HMI_LOS_IMAGE

smap = sunpy.map.Map(HMI_LOS_IMAGE)

fig = plt.figure()
ax = fig.add_subplot(111, projection=smap)

smap.plot(axes=ax)

plt.show()