def make_dark(path=None, bias=None): """ Make master dark from raw data files """ if not path: path = set_path() if length(bias) == 1: bias = make_bias(path=path) darks,dct = tp.get_files(dir = path,tag='Dark') masterdark = tp.master_dark(darks,bias=bias,outdir='./') return masterdark
y0 = pix0[0,1] plt.figure(1) plt.plot([x0],[y0],'y+',markersize=30) plt.xlim(0,2048) plt.ylim(0,2048) darks,dct = tp.get_files(tag='Dark') biases,bct = tp.get_files(tag='Bias') flats,fct = tp.get_files(tag='SkyFlat') targetfiles,tct = tp.get_files(tag='KIC10935310') bias = tp.master_bias(biases,outdir='./Photometry/',readnoise=False) dark = tp.master_dark(darks,bias=bias,outdir='./Photometry/') flat = tp.master_flat(flats,bias=bias,dark=dark,outdir='./Photometry/') cal = (image - dark - bias)/flat qi.display_image(cal,siglo=2,sighi=2) plt.plot([x0],[y0],'y+',markersize=30,linewidth=1e6) plt.xlim(0,2048) plt.ylim(0,2048) # define the aperture radius,ynew,xnew,fwhm,aspect,snrmax,totflux,totap,chisq = \ tp.optimal_aperture(x0,y0,image,skyrad=[15,20]) pref = np.array([[x0,y0],[836,1154],[645,1152], [1234,1002], [1131,1338], [377,426]]) coords0 = w.wcs_pix2world(pref,1)