In [1]:

    

import numpy as np
from matplotlib import pyplot as plt
%matplotlib inline
plt.rcParams['figure.figsize'] = 20,15
from myPyLib.Mollweide import moll


    

In [2]:

    

xao = np.loadtxt('field.xao.out.txt',skiprows=1)
bok = np.loadtxt('field.bok.out.txt',skiprows=1)


    

In [3]:

    

xra0 = [f[1] for f in xao if f[5] == 0]
xde0 = [f[2] for f in xao if f[5] == 0]
xgl0 = [f[3] for f in xao if f[5] == 0]
xgb0 = [f[4] for f in xao if f[5] == 0]
xra1 = [f[1] for f in xao if f[5] == 1]
xde1 = [f[2] for f in xao if f[5] == 1]
xgl1 = [f[3] for f in xao if f[5] == 1]
xgb1 = [f[4] for f in xao if f[5] == 1]
xra2 = [f[1] for f in xao if f[5] == 2]
xde2 = [f[2] for f in xao if f[5] == 2]
xgl2 = [f[3] for f in xao if f[5] == 2]
xgb2 = [f[4] for f in xao if f[5] == 2]

bra0 = [f[1] for f in bok if f[5] == 0]
bde0 = [f[2] for f in bok if f[5] == 0]
bgl0 = [f[3] for f in bok if f[5] == 0]
bgb0 = [f[4] for f in bok if f[5] == 0]
bra1 = [f[1] for f in bok if f[5] == 1]
bde1 = [f[2] for f in bok if f[5] == 1]
bgl1 = [f[3] for f in bok if f[5] == 1]
bgb1 = [f[4] for f in bok if f[5] == 1]


    

In [4]:

    

mxe = moll(lat_range=(-25,80))


    

In [12]:

    

plt.rcParams['figure.figsize'] = 16,12
mxe.grid(lat_lab_lon=0,lon_lab_lat=-22)
mxe.scatter(xra0, xde0, 'bs', label='Abondoned')
mxe.scatter(xra1, xde1, 'rs', label='SDSS')
mxe.scatter(xra2, xde2, 'gs', label='To Do')
plt.legend()
plt.title('XAO Obs Plan in Equatorial System')
plt.savefig('xao.1.eps')
plt.savefig('xao.1.png')


    

In [6]:

    

mxg = moll(lat_range=(-80,80))


    

In [13]:

    

plt.rcParams['figure.figsize'] = 16,12
mxg.grid(lat_lab_lon=0,lon_lab_lat=0)
mxg.scatter(xgl0, xgb0, 'bs', label='Abondoned')
mxg.scatter(xgl1, xgb1, 'rs', label='SDSS')
mxg.scatter(xgl2, xgb2, 'gs', label='To Do')
plt.legend()
plt.title('XAO Obs Plan in Galactic System')
plt.savefig('xao.2.eps')
plt.savefig('xao.2.png')


    

In [9]:

    

mbe = moll(lat_range=(-5,80))


    

In [14]:

    

plt.rcParams['figure.figsize'] = 16,12
mbe.grid(lat_lab_lon=0,lon_lab_lat=-3)
mbe.scatter(bra0, bde0, 'bs', label='Abondoned')
mbe.scatter(bra1, bde1, 'gs', label='To Do')
plt.legend()
plt.title('Bok Obs Plan in Equatorial System')
plt.savefig('bok.1.eps')
plt.savefig('bok.1.png')


    

In [11]:

    

mbg = moll(lat_range=(-80,80))


    

In [15]:

    

plt.rcParams['figure.figsize'] = 16,12
mbg.grid(lat_lab_lon=0,lon_lab_lat=-3)
mbg.scatter(bgl0, bgb0, 'bs', label='Abondoned')
mbg.scatter(bgl1, bgb1, 'gs', label='To Do')
plt.legend()
plt.title('Bok Obs Plan in Galactic System')
plt.savefig('bok.1.eps')
plt.savefig('bok.1.png')


    

In [ ]: