# Apply flatfield img1 = flat.ImgNormalize(image, Method="mean", source="image") #img1 = flat.ImgNormalize(img1, Method="mean", source="image") img2 = sm.colmeansub(image) # Find stars centers = centroid.findstars(img1, zreject=3, zthresh=3.0, zpeakthresh=4, min_pix_per_star=6, max_pix_per_star=60, oblongness=5,debug=False) C = centroid.imgcentroid(image,centers) #img3 = sm.colmeansub(img2) #img4 = sm.colmeansub(image2) #img3 = np.delete(img3, np.r_[1080:1080+32], 0) #imgutil.dispimg(image,1) #imgutil.dispimg(img1,1) #imgutil.dispimg(img2,6) #imgutil.dispimg(sm.bgsub(img1), 6) #imgutil.dispimg(sm.colmeansub(img1), 6) imgutil.circstars(sm.colmeansub(img1), C, viewfactor=3)
print 'Opening: ' + fname image = imgutil.loadimg(fname, from_database=True) tic = time.clock() # Clean-up the image: #image = flatfield.NormalizeColumnGains(image,Plot=0,Wiener=0) image = flatfield.ImgNormalize(image, Method="mean") toc = time.clock() # Find star centroids: centroids = [] (centers) = centroid.findstars(image, zreject=3, zthresh=3.0, zpeakthresh=4, min_pix_per_star=6, max_pix_per_star=60, oblongness=2, debug=True) centroids = centroid.imgcentroid(image, centers) # Display image with stars circled: vf = 3 if centroids: imgutil.circstars(image, centroids, 20, viewfactor=vf) else: print 'Oopsies... we found ZERO centroids. Now how do you feel?' imgutil.dispimg(image, viewfactor=vf) print 'Total time: ', toc - tic, ' s'
################################################################################### # Load the image: db = database.Connect() fnames = db.select('select raw_fn from rawdata where burst_num = 185 limit 2').raw_fn.tolist() fnames.sort() fname = fnames[0] print 'Opening: ' + fname image = imgutil.loadimg(fname,from_database=True) tic = time.clock() # Clean-up the image: #image = flatfield.NormalizeColumnGains(image,Plot=0,Wiener=0) image = flatfield.ImgNormalize(image, Method="mean") toc = time.clock() # Find star centroids: centroids = [] (centers) = centroid.findstars(image,zreject=3, zthresh=3.0, zpeakthresh=4, min_pix_per_star=6, max_pix_per_star=60, oblongness=2,debug=True) centroids = centroid.imgcentroid(image,centers) # Display image with stars circled: vf = 3 if centroids: imgutil.circstars(image,centroids,20,viewfactor=vf) else: print 'Oopsies... we found ZERO centroids. Now how do you feel?' imgutil.dispimg(image,viewfactor=vf) print 'Total time: ',toc - tic,' s'