ó
TûTYc           @   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   sƒ  t  d d ƒ } t  d d ƒ } y3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  @r)d t | ƒ d t | ƒ d } t | ƒ d t | ƒ } | j | d ƒ | j | d ƒ n  | d } qiW| j ƒ  | j ƒ  | GHWn+ t k
 rid GHn t k
 r~d GHn Xd  S(   Ns   regions_stars_cfhtm.regt   ws   assoc_stars_list_cfhtm.txti    i   t   rat   dect   framet   icrsi   t   Rowss   circle(t   ,s   ,7.0")t   	s   
t   file_errt   valuefile_err(   t   openR   R   t   headert   datat   shapet   wcs_pix2worldR   t   ut   degreet
   separationt   arcmint   search_sdss_start   wcs_world2pixt   strt   writet   closet   IOErrort
   ValueError(   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(    (    s.   /db2/Data/gris/Fitting/garik7/make_assoc_cm.pyt   create_assoc
   sL    ))

*"

	(   t   request_sdsst   numpyt   npt
   astropy.ioR   t   astropy.wcsR   t   astropy.coordinatesR   t   astropyR   R   t   osR9   (    (    (    s.   /db2/Data/gris/Fitting/garik7/make_assoc_cm.pyt   <module>   s   
