B
    %d_>                 @   sP   d dl Zd dlmZ d dlmZ d dlmZ ddd	Zd
d ZedkrLe  dS )    N)shift)bin_2d)rebin
   皙?   c             C   s   t | |}t|jd |jd f}tj|dd}tj|dd}	||| |	 k  ||| |	 k < t| j}
t|||
d||jd  d||jd  f< |
S )Nr   r   )axis)r   npzerosshapemedianstdr   )Z
imcube_inpZbgr_inpmaskZbinsizeZnegsigZn_chbbgr_bZm_valZs_valbgr r   "/gal/db1/Data/HiPS/clean_negbgr.pyclean_negbgr   s    
 .r   c              C   sl   ddl m}  | d}| d}t|d j|d jd |d jd }|d j| |d _|jddd d S )	Nr   )fitszi/data1/rcsed_v2/HECTOSPEC/MILES/2015.0414/skysub_cosmos15a2_1/spSpec279.cosmos15a2_1_1527.ms_results.fitsz5images/UKIDSS_J125957.20+280300.00/swarp/K_COADD.fitsg        ztest_bgr.fitsT)	overwrite)
astropy.ior   openr   datawriteto)r   hhhr   r   r   r   test_cnegbg   s    

&r   __main__)r   r   r   )	numpyr
   Zscipy.ndimage.interpolationr   r   r   r   r   __name__r   r   r   r   <module>   s   
	