def make_bias(path=None): """ Make master bias from raw data files """ if not path: path = set_path() biases,bct = tp.get_files(dir=path, tag='Bias') masterbias = tp.master_bias(biases,outdir='./') return masterbias
pix0 = w.wcs_world2pix(world0,1) x0 = pix0[0,0] 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])
import thacherphot as tp import matplotlib as mpl from matplotlib.dates import DayLocator, HourLocator, DateFormatter, drange import matplotlib.gridspec as gridspec import datetime import jdutil from astropysics.coords import AngularCoordinate as angcor path = '/Users/jonswift/Astronomy/ThacherObservatory/data/22Dec2016/' #path = '/Users/ONeill/astronomy/Data/04Dec2016/' files,fct = tp.get_files(dir=path, prefix='WASP-50b',suffix='solved.fits') bfiles,bct = tp.get_files(dir=path,prefix='WASP-50b',suffix='bias.fit') master_bias = tp.master_bias(bfiles) data = tp.batch_total_flux(files,object='WASP-50b',mag=11.6,SpT='g8v',bias=master_bias, nearest=True,camera='iKON-L',filter='V',gain=3.0,network='swift') #tau = data['tau'] #secz = data['secz'] #plt.clf() #plt.ion() #plt.figure() #plt.plot(secz,tau,'o') #plt.ylabel('tau') #plt.xlabel('secz') #plt.xlim(1.85,1.95) #plt.ylim(0.64,0.67)