# Path: ugali/preprocess/maglims.py
# def inMangle(polyfile,ra,dec):
#     coords = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     logger.debug("Writing coordinates to %s"%coords.name)
#     np.savetxt(coords, np.array( [ra,dec] ).T, fmt='%.6g' )
#     coords.close()
# 
#     weights = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     cmd = "polyid -W %s %s %s"%(polyfile,coords.name,weights.name)
#     logger.debug(cmd)
#     subprocess.call(cmd,shell=True)
# 
#     tmp = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     cmd = """awk '{if($3==""){$3=0} print $1, $2, $3}' %s > %s"""%(weights.name,tmp.name)
#     logger.debug(cmd)
#     subprocess.call(cmd,shell=True)
# 
#     data = np.loadtxt(tmp.name,unpack=True,skiprows=1)[-1]
#     for f in [coords,weights,tmp]:
#         logger.debug("Removing %s"%f.name)
#         os.remove(f.name)
# 
#     return data > 0
# 
# Path: ugali/utils/projector.py
# class SphericalRotator:
# class Projector:
#     def __init__(self, lon_ref, lat_ref, zenithal=False):
#     def setReference(self, lon_ref, lat_ref, zenithal=False):
#     def cartesian(self,lon,lat):
#     def rotate(self, lon, lat, invert=False):
#     def __init__(self, lon_ref, lat_ref, proj_type = 'ait'):
#             def rotate(lon,lat,invert=False):
#     def sphereToImage(self, lon, lat):
#     def imageToSphere(self, x, y):
# def sphere2image(lon_ref,lat_ref,lon,lat):
# def image2sphere(lon_ref,lat_ref,x,y):
# def cartesianSphereToImage(lon, lat):
# def cartesianImageToSphere(x,y):
# def aitoffSphereToImage(lon, lat):
# def aitoffImageToSphere(x, y):
# def gnomonicSphereToImage(lon, lat):
# def gnomonicImageToSphere(x, y):
# def angsep2(lon_1, lat_1, lon_2, lat_2):
# def angsep(lon1,lat1,lon2,lat2):
# def galToCel(ll, bb):
# def celToGal(ra, dec):
# def estimate_angle(angle, origin, new_frame, offset=1e-7):
# def gal2cel_angle(glon,glat,angle,offset=1e-7):
# def cel2gal_angle(ra,dec,angle,offset=1e-7):
# def dec2hms(dec):
# def dec2dms(dec):
# def hms2dec(hms):
# def dms2dec(dms):
# def sr2deg(solid_angle):
# def deg2sr(solid_angle):
# def distanceToDistanceModulus(distance):
# def distanceModulusToDistance(distance_modulus):
# def ang2const(lon,lat,coord='gal'):
# def ang2iau(lon,lat,coord='gal'):
# def match(lon1, lat1, lon2, lat2, tol=None, nnearest=1):
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.


import healpy
import numpy as np
import pylab as plt
from ugali.preprocess.maglims import inMangle
from ugali.utils.projector import pixToAng
    from optparse import OptionParser
#!/usr/bin/env python
#import mangle


if __name__ == "__main__":
    usage = "Usage: %prog  [options] mangle1.ply [mangle2.ply ...]"
    description = "python script"
    parser = OptionParser(usage=usage,description=description)
    (opts, args) = parser.parse_args()

    nside = 2**8
    print("Creating ra, dec...")
    pix = np.arange(healpy.nside2npix(nside))

# Path: ugali/preprocess/maglims.py
# def inMangle(polyfile,ra,dec):
#     coords = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     logger.debug("Writing coordinates to %s"%coords.name)
#     np.savetxt(coords, np.array( [ra,dec] ).T, fmt='%.6g' )
#     coords.close()
# 
#     weights = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     cmd = "polyid -W %s %s %s"%(polyfile,coords.name,weights.name)
#     logger.debug(cmd)
#     subprocess.call(cmd,shell=True)
# 
#     tmp = tempfile.NamedTemporaryFile(suffix='.txt',delete=False)
#     cmd = """awk '{if($3==""){$3=0} print $1, $2, $3}' %s > %s"""%(weights.name,tmp.name)
#     logger.debug(cmd)
#     subprocess.call(cmd,shell=True)
# 
#     data = np.loadtxt(tmp.name,unpack=True,skiprows=1)[-1]
#     for f in [coords,weights,tmp]:
#         logger.debug("Removing %s"%f.name)
#         os.remove(f.name)
# 
#     return data > 0
# 
# Path: ugali/utils/projector.py
# class SphericalRotator:
# class Projector:
#     def __init__(self, lon_ref, lat_ref, zenithal=False):
#     def setReference(self, lon_ref, lat_ref, zenithal=False):
#     def cartesian(self,lon,lat):
#     def rotate(self, lon, lat, invert=False):
#     def __init__(self, lon_ref, lat_ref, proj_type = 'ait'):
#             def rotate(lon,lat,invert=False):
#     def sphereToImage(self, lon, lat):
#     def imageToSphere(self, x, y):
# def sphere2image(lon_ref,lat_ref,lon,lat):
# def image2sphere(lon_ref,lat_ref,x,y):
# def cartesianSphereToImage(lon, lat):
# def cartesianImageToSphere(x,y):
# def aitoffSphereToImage(lon, lat):
# def aitoffImageToSphere(x, y):
# def gnomonicSphereToImage(lon, lat):
# def gnomonicImageToSphere(x, y):
# def angsep2(lon_1, lat_1, lon_2, lat_2):
# def angsep(lon1,lat1,lon2,lat2):
# def galToCel(ll, bb):
# def celToGal(ra, dec):
# def estimate_angle(angle, origin, new_frame, offset=1e-7):
# def gal2cel_angle(glon,glat,angle,offset=1e-7):
# def cel2gal_angle(ra,dec,angle,offset=1e-7):
# def dec2hms(dec):
# def dec2dms(dec):
# def hms2dec(hms):
# def dms2dec(dms):
# def sr2deg(solid_angle):
# def deg2sr(solid_angle):
# def distanceToDistanceModulus(distance):
# def distanceModulusToDistance(distance_modulus):
# def ang2const(lon,lat,coord='gal'):
# def ang2iau(lon,lat,coord='gal'):
# def match(lon1, lat1, lon2, lat2, tol=None, nnearest=1):
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.
#     DEGREE = 360.
#     HOUR = 24.
#     MINUTE = 60.
#     SECOND = 3600.


# Path: ugali/scratch/PlotManglePolygon.py
import healpy
import numpy as np
import pylab as plt
from ugali.preprocess.maglims import inMangle
from ugali.utils.projector import pixToAng
    from optparse import OptionParser

#!/usr/bin/env python
#import mangle


if __name__ == "__main__":
    usage = "Usage: %prog  [options] mangle1.ply [mangle2.ply ...]"
    description = "python script"
    parser = OptionParser(usage=usage,description=description)
    (opts, args) = parser.parse_args()

    nside = 2**8
    print("Creating ra, dec...")
    pix = np.arange(healpy.nside2npix(nside))

    ra,dec = pixToAng(nside,pix)