ó
˘(ËWc           @   so   d  d l  Td  d l Z d  d l m Z d  d l m Z d  d l m Z d  d l	 m
 Z d  d l Z d   Z d S(   i˙˙˙˙(   t   *N(   t   fits(   t   WCS(   t   SkyCoord(   t   unitsc         C   sQ  t  d d  } t  d d  } t j  |   } t | d j  } | d j } | d j } | j d } | j d } | j d d d  \ }	 }
 | j | | d  \ } } |	 GHt d |	 t j	 d |
 t j	 d	 d
  } t d | t j	 d | t j	 d	 d
  } | j
 |  } | j d } |	 | d } |
 | d } t | | |  } | d d } d } xŇ | D]Ę } | d } | d } | j | | d  \ } } | d k | d k @| | k  @| | k  @rfd t |  d t |  d } t |  d t |  } | j | d  | j | d  | d } qfqfW| j   | j   | GHd  S(   Ns   regions_stars_cfhtm.regt   ws   assoc_stars_list_cfhtm.txtt   SCIi    i   t   rat   dect   framet   icrsi   t   Rowss   circle(t   ,s   ,7.0")t   	s   
(   t   openR   R   t   headert   datat   shapet   wcs_pix2worldR   t   ut   degreet
   separationt   arcmint   search_sdss_start   wcs_world2pixt   strt   writet   close(   t   namet   f_regt   f_assoct   hdulistR   t   hdrt   datt   widtht   heightt   ra1t   dec1t   ra2t   dec2t   c1t   c2t   sept   radt   ra_ct   dec_ct
   stars_jsont   jsont   it   rowR   R   t   xt   yt   string_to_reg_filet   string_to_assoc_file(    (    s6   /db2/Data/gris/Fitting/super_objects/make_assoc_hst.pyt   create_assoc
   sB    ))

*"

(   t   request_sdsst   numpyt   npt
   astropy.ioR   t   astropy.wcsR   t   astropy.coordinatesR   t   astropyR   R   t   osR6   (    (    (    s6   /db2/Data/gris/Fitting/super_objects/make_assoc_hst.pyt   <module>   s   
