%matplotlib inline
import vaex as vx
from astropy import units as u
from astropy.coordinates import SkyCoord
from astropy.table import Column, Table, join
from astropy.cosmology import FlatLambdaCDM

import numpy as np

import ipyvolume as ipv

overview = Table.read("../dmu32/dr1_data_products_overview_20180627.csv")

overview

Table length=23

field	dr1_suffix	last_suffix	objects	area_sq_degrees	file_size_bytes	file_size_readable	xid_objects	photoz_objects	cigale_objects
str18	int64	int64	int64	float64	int64	str9	int64	int64	int64
AKARI-NEP	20180215	20180215	531746	9.19473235878	348848640	348.85 MB	0	0	0
AKARI-SEP	20180221	20180221	844172	8.71330647513	538606080	538.61 MB	0	0	0
Bootes	20180520	20180520	3481661	11.428152991	5002044480	5 GB	0	0	0
CDFS-SWIRE	20180613	20180613	2171051	12.9712464037	5705593920	5.71 GB	240489	139532	9308
COSMOS	20180619	20180619	2599374	5.0838634785	10517175360	10.52 GB	44310	694478	15747
EGS	20180501	20180501	1412613	3.56638327512	3808480320	3.81 GB	0	1183490	0
ELAIS-N1	20171016	20180216	4026292	13.5074845555	5242279680	5.24 GB	120282	2885116	50129
ELAIS-N2	20180218	20180218	1783240	9.16747990399	2007964800	2.01 GB	0	0	0
ELAIS-S1	20180416	20180416	1655564	9.00294064689	1970156160	1.97 GB	194276	1102319	25393
GAMA-09	20180601	20180601	12937982	62.0139341728	19911611520	19.91 GB	1386659	8834023	130293
...	...	...	...	...	...	...	...	...	...
Herschel-Stripe-82	20180307	20180307	50196455	363.227479306	127047306240	127.05 GB	0	21613821	0
Lockman-SWIRE	20180219	20180219	4366298	22.4138634029	5195931840	5.2 GB	242065	1382438	46719
NGP	20180219	20180501	6759591	177.698197108	4954806720	4.95 GB	0	3175339	0
SA13	20180501	20180501	9799	0.274163793284	2721600	2.72 MB	0	0	0
SGP	20180221	20180221	29790690	294.56845409	45192530880	45.19 GB	3511594	17057212	0
SPIRE-NEP	20180220	20180220	2674	0.128066046938	1203840	1.2 MB	0	0	0
SSDF	20180221	20180221	12661903	111.116249493	8078319360	8.08 GB	0	9268783	0
xFLS	20180501	20180501	977148	7.44043242141	1050474240	1.05 GB	0	0	0
XMM-13hr	20180501	20180501	38629	0.757024016658	10676160	10.68 MB	0	0	0
XMM-LSS	20180504	20180504	8705837	21.7545793933	28241838720	28.24 GB	0	6152920	0

overview['field', 
         'area_sq_degrees', 
         'objects', 
         'xid_objects', 
         'photoz_objects', 
         'cigale_objects'].write('overview_table.tex')

WARNING: AstropyDeprecationWarning: overview_table.tex already exists. Automatically overwriting ASCII files is deprecated. Use the argument 'overwrite=True' in the future. [astropy.io.ascii.ui]
WARNING:MainThread:astropy:AstropyDeprecationWarning: overview_table.tex already exists. Automatically overwriting ASCII files is deprecated. Use the argument 'overwrite=True' in the future.

np.sum()

from astropy.table import Table, Column
from astropy import units as u
from astropy.io import fits
from astropy.coordinates import SkyCoord

import numpy as np

import re

import yaml

import glob

from shutil import copyfile


import os
import time

import humanfriendly

def name_to_ra_dec(help_id):
    """Take a HELP id and return the ra dec in deg
    
    parameters
    ----------
    
    help_id: string
       The help ID
       
    returns
    -------
    
    ra: flaot
       Right Ascension in degrees
       
    dec: float
        Declination in degrees
    """
    if '-' in help_id:
        ra, dec = help_id.strip('HELP_J').split('-')
        dec = '-' + dec
    if '+' in help_id:
        ra, dec = help_id.strip('HELP_J').split('+')
        dec = '+' + dec
    
    ra = ra[:2] + 'h' + ra[2:]
    ra = ra[:5] + 'm' + ra[5:]
    ra = ra + 's' 
    
    
    
    dec = dec[:3] + 'd' + dec[3:]
    dec = dec[:6] + 'm' + dec[6:]
    dec = dec + 's'
    
    coords = SkyCoord('{} {}'.format(ra, dec))
    
    return coords.ra.value, coords.dec.value

ids = Table.read('cdfs_ids.fits')

ra_max = 0.0
for gal in ids:
    ra, dec = name_to_ra_dec(gal[0])
    if ra > ra_max:
        ra_max = ra
    if ra == np.nan:
        print(ra, dec, 'bad ra')
    if dec == np.nan:
        print(ra, dec, 'bad dec')
    

ra_max

51.701924999999996