def load_data(image,ra,dec,hsize): try: xcen,ycen = mi.get_center_coords(image,ra,dec,hsize) hdu=pyfits.open(image) image_stamp= hdu[0].data[ycen-hsize:ycen+hsize, xcen-hsize:xcen+hsize] hdu.close() except IndexError as err: if args.verbose: print(err) image_stamp=np.zeros([hsize*2,hsize*2]) return image_stamp
if Z[i]<6.0: color=MagRef1[i]-Mag[i] if MagRef1[i]==Mag[i]==-99: color=-99 else: color=MagRef2[i]-Mag[i] if MagRef2[i]==Mag[i]==-99: color=-99 image_stamp, Xgal, Ygal = load_data(imgname,RA[i],DEC[i],hsize) if args.zphot: sel_radec = (table_zphot['alpha']<RA[i]+args.size/3600.)*(table_zphot['alpha']>RA[i]-args.size/3600.)*(table_zphot['delta']>DEC[i]-args.size/3600.)*(table_zphot['delta']<DEC[i]+args.size/3600.) Positions=[] for ra,dec,z,zl,zu in zip(table_zphot['alpha'][sel_radec],table_zphot['delta'][sel_radec],z_phot[sel_radec],z_l68[sel_radec],z_u68[sel_radec]): X,Y=mi.get_center_coords(imgname,ra,dec,hsize,verify_limits=False) Positions.append([(X-Xgal)+hsize,(Y-Ygal)+hsize,z,zl,zu]) Positions = np.array(Positions) else: Positions=None FullCatalog[ID[i]] = find_pairs_and_clumps(image_stamp,Z[i],Mag[i],color,hsize,t,fractions,sblim,args.pixscale,args.zeropoint,ksky=args.ksky,Areamin=args.areamin,no_dilation=args.nodilation,degrade=args.degrade,Aperture=args.aperture,title="VUDS %i @ z=%.4f [%i]"%(ID[i],Z[i],Zflag[i]),plot_results=args.nosaving,erosion=eroders,verbose=args.verbose,ident=args.ident,zphot_sel=Positions) if not args.nosaving: with open(table_name, 'wb') as outfile: pickle.dump(FullCatalog, outfile, protocol=pickle.HIGHEST_PROTOCOL)