In [1]:
import opsimsummary
In [3]:
from opsimsummary import AllSkyMap,
In [4]:
%matplotlib inline
import matplotlib.pyplot as plt
In [5]:
import numpy as np
In [16]:
fig, ax = plt.subplots()
m = AllSkyMap(projection='moll', ax=ax, lon_0=0, lat_0=0)
_ = m.drawparallels(np.arange(-91.,91.,20.))
_ = m.drawmeridians(np.arange(-180., 181., 60.))#, labels=[1, 0, 0, 1], labelstyle='+/-')
ax.set_title('{:.6f}'.format(23.4))
Out[16]:
In [7]:
pts = m.tissot(lon_0=0, lat_0=0., radius_deg=4., npts=100, ax=ax,
**dict(fill=False, edgecolor='r', lw=2))
In [8]:
pts = m.tissot(lon_0=180, lat_0=0., radius_deg=4., npts=100, ax=ax,
**dict(fill=False, edgecolor='r', lw=2))
In [9]:
pts = m.tissot(lon_0=178, lat_0=20., radius_deg=4., npts=100, ax=ax,
**dict(fill=False, edgecolor='r', lw=2))
In [10]:
pts = m.tissot(lon_0=-178, lat_0=20., radius_deg=4., npts=100, ax=ax,
**dict(fill=False, edgecolor='r', lw=2))
In [11]:
pts = m.tissot(lon_0=182, lat_0=-20., radius_deg=4., npts=100, ax=ax,
**dict(fill=False, edgecolor='r', lw=2))
In [12]:
fig
Out[12]:
In [ ]:
In [ ]:
In [ ]: