ó
§u,Xc           @   sC   d  d l  m Z d  d l Z d  d l m Z d  d l Td   Z d S(   i˙˙˙˙(   t   fitsN(   t   interpolate(   t   *c         C   s_  t  d  } i  } d | d <d | d <d | d <d | d	 <d
 | d <| d } t j |   } | d j } | d j d } | d j d } | d j d } | d j d }	 | d j d }
 | j GHt j t j	 d | j d d  t j	 d | j d d  | d d } t
 |	 d  } t
 |
 d  } t j | | f d t j } |	 j GH|
 j GH| j GH| |	 d |
 d  } | d j } t j | | f d t j } | | (| | | } | j GHt j |  } | | } | d j } i  } xb | D]Z } | d | d k | d | d k @| d | d k @| d | d k @rö| } PqöqöW| d | } | d | } | Gd G| GHt j | | |  } | | } t j | j   d d } t j | j   d d } | Gd G| Gd G| | GHt j d |  } t j | g  } | j d d t t j d |  } t j | g  } | j d  d t d  S(!   Nt
   frame_infoi    t   ui   t   gi   t   ri   t   ii   t   zt   photo_fieldt   FILTERt   ALLSKYt   XINTERPt   YINTERPt   kindt   lineart   dtypet   FIELDt   CAMCOLt   RERUNt   RUNt   DARK_VARIANCEt   GAINs      t   axiss     t   datas   err.fitst   clobbers   sky_stretched.fits(   t   read_fiR    t   openR   t   headert   shapeR   t   interp2dt   npt   aranget   lent   zerost   floatt   sqrtt   stdt   ravelt   meant
   PrimaryHDUt   HDUListt   writetot   True(   t	   framenamet   infot   filters_numbert   photofieldnamet   hdulistt   scidatat   filtt	   fits_infot   skyt   x_arrt   y_arrt	   interpolet   heightt   widtht   sky_stretchedt	   calib_vect   calibrationt   dnt   hdulist_photoFieldt   filt_numt   datt   photoField_rowt   rowt   dark_variancet   gaint   dn_errt   img_errt   stddivR'   t   hdu(    (    s0   /db2/Data/gris/Fitting/super_objects/calc_err.pyt   calc_err_fr   sf    





G
J
(   t
   astropy.ioR    t   numpyR   t   scipyR   t   read_frameinfoRI   (    (    (    s0   /db2/Data/gris/Fitting/super_objects/calc_err.pyt   <module>   s   
