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)