import pylab as plt
import numpy as np
from astropy import units as u
from mw_plot import MWSkyMap

# setup a MWSkyMap instance with projection, other projection can be 'hammer', 'mollweide' etc
# grayscale: whether to turn the background image to grayscale
mw1 = MWSkyMap(projection="aitoff", grayscale=False)

# set up plot title
mw1.title = "LMC and SMC in red dots"

# LMC and SMC coordinates
lsmc_ra = [78.77, 16.26] * u.degree
lsmc_dec = [-69.01, -72.42] * u.degree

mw1.scatter(lsmc_ra, lsmc_dec, c="r", s=200)
plt.tight_layout()