error = False
if not os.path.isfile(simage):
try:
show_overlay(lhdu,
shdu,
ra,
dec,
size,
firsthdu=None,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=simage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='magenta',
title=title,
lw_ellipse=3,
ellipse_style=ls,
ellipse_color=cs,
peak=peak,
ppsize=350)
except (Exception, RuntimeError) as e:
error = True
if error:
try:
show_overlay(lhdu,
lhdu,
ra,
dec,
size,
firsthdu=None,
overlay_cat=ots,
show_grid=False,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=cimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title,
lw_ellipse=3,
ellipse_style=ls,
ellipse_color=cs,
peak=peak,
ppsize=350,
cmap='inferno',
lofar_colorscale=True,
show_lofar=False)
show_overlay(lhdu,
peak = None
plist = []
if not os.path.isfile(iimage):
show_overlay(lhdu,
ihdu,
ra,
dec,
size,
firsthdu=None,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=iimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title,
lw_ellipse=3,
ellipse_style=ls,
peak=peak,
noisethresh=1)
plist.append(mogrify(iimage))
if not os.path.isfile(ipimage):
show_overlay(lhdu,
ihdu,
del(gals)
'''
ots=ot[separation(ra,dec,ot['RA'],ot['DEC'])<(size*2)]
#pshdu=extract_subim(imagedir+'/downloads/'+psmaps[i],ra,dec,size*2,hduid=1)
print 'Lofarfile is',lofarfile
lhdu=extract_subim(lofarfile,ra,dec,size*2)
firsthdu=extract_subim(imagedir+'/downloads/'+firstmaps[i],ra,dec,size*2)
whdu=extract_subim(imagedir+'/downloads/'+wisemaps[i],ra,dec,size)
try:
peak==r['Peak_flux']/1000.0
except:
peak=None
f=show_overlay(lhdu,whdu,ra,dec,size,firsthdu=firsthdu,coords_color='red',coords_ra=r['RA'],coords_dec=r['DEC'],coords_lw=3,lw=2,no_labels=True,marker_ra=marker_ra,marker_dec=marker_dec,marker_lw=3,marker_color='cyan',title=title,block=False,interactive=False,drlimit=8000,peak=peak)
ora,odec=ss.odict[sourcename]
components=ss.get_comps(sourcename)
notcomponents=ss.get_ncomps(sourcename)
print 'components is',components
I=Interactive(f,ora,odec,ots,components,notcomponents)
fig=plt.gcf()
fig.canvas.mpl_connect('button_press_event', I.onclick)
plt.ion()
plt.show(block=False)
try:
print 'Blend probability',r['Blend_prob']
except:
pass
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=None,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_lw=3,
lw=1,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title,
peak=peak,
plot_coords=False,
show_grid=False,
lw_ellipse=3,
ellipse_style=ls,
ellipse_color=colors,
noisethresh=1.5,
drlimit=1000,
rms_use=rms,
lofarlevel=2.5,
logfile=logfile,
sourcename=sourcename,
vmax_cap=0.5)
stop = False
while not (stop):
if display_mode == 'Spitzer':
ohdu = extract_subim(spitzerfile, ra, dec, size)
else:
ohdu = extract_subim(ibandfile, ra, dec, size)
f = show_overlay(lhdu,
ohdu,
ra,
dec,
size,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
no_labels=True,
title=title,
block=False,
interactive=False,
drlimit=8000,
peak=peak,
plotpos=[(pwg, '+'), (mlt, 'h', 'red')])
I = Interactive(f, ga)
fig = plt.gcf()
fig.canvas.mpl_connect('button_press_event', I.onclick)
plt.ion()
plt.show(block=False)
mode='wise'
whdu=extract_subim(wisefile,ra,dec,size*scalefactor)
else:
whdu=extract_subim(legacyfile,ra,dec,size*scalefactor)
if np.all(np.isnan(whdu[0].data)):
print('Legacy image is blanked, using WISE')
mode='wise'
whdu=extract_subim(wisefile,ra,dec,size*scalefactor)
try:
peak==r['Peak_flux']/1000.0
except:
peak=None
try:
f=show_overlay(lhdu,whdu,ra,dec,size*scalefactor/initfactor,coords_color='red',coords_ra=s.sd[sourcename]['RA'],coords_dec=s.sd[sourcename]['DEC'],coords_lw=3,lw=2,no_labels=True,marker_ra=marker_ra,marker_dec=marker_dec,marker_lw=3,marker_color='cyan',title=title,block=False,interactive=False,drlimit=8000,peak=peak,plotpos=[(pwg,'+')])
except Exception as e:
print('*** Error making image. Press RETURN to skip and continue ***')
print(e)
raw_input()
stop=True
break
ora,odec=s.optid(sourcename)
components=s.get_comps(sourcename)
notcomponents=[]
for r in ots:
name=r['Source_Name']
if name in s.cd and s.cd['Parent']!=sourcename:
notcomponents.append(name)
size,
hduid=1)
lhdu = extract_subim(lofarfile, ra, dec, size)
firsthdu = extract_subim(imagedir + '/downloads/' + firstmaps[i], ra, dec,
size)
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=firsthdu,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title)
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
dec,
size,
hduid=1)
lhdu = extract_subim(lofarfile, ra, dec, size)
firsthdu = extract_subim(imagedir + '/downloads/' + firstmaps[i], ra, dec,
size)
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=firsthdu,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='magenta',
title=title,
peak=r['Peak_flux'] / 1000.0)
os.system('mogrify -quality 90 -trim ' + sourcename + '*.png')
size)
try:
peak == r['Peak_flux'] / 1000.0
except:
peak = None
print ots, cols
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=firsthdu,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title,
ellipse_color=cols,
peak=peak)
pshdu = extract_subim(imagedir + '/downloads/' + psmaps[name],
ra,
dec,
size,
hduid=1)
lhdu = extract_subim(lofarfile, ra, dec, size)
firsthdu = extract_subim(imagedir + '/downloads/' + firstmaps[name],
ra, dec, size)
plot_positions = check_positions[name]
print 'Plotting positions', plot_positions
plot_ra = []
plot_dec = []
for p in plot_positions:
plot_ra.append(p[0])
plot_dec.append(p[1])
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=firsthdu,
overlay_cat=ot,
overlay_scale=scale,
lw=2,
save_name=outname,
no_labels=True,
coords_ra=plot_ra,
coords_dec=plot_dec,
coords_lw=2)
color = bits[2]
plotpos = [(tpp, shape)]
print '*** PLOTPOS is ', plotpos
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=firsthdu,
vlasshdu=vlasshdu,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_ra=r['RA'],
coords_dec=r['DEC'],
coords_lw=3,
lw=2,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='magenta',
title=title,
ellipse_color=cols,
peak=peak,
plotpos=plotpos)
whdu = extract_subim(imagedir + '/downloads/' + wisemaps[i], ra, dec,
size)
show_overlay(lhdu,
whdu,
peak == r['Peak_flux'] / 1000.0
except:
peak = None
show_overlay(lhdu,
pshdu,
ra,
dec,
size,
firsthdu=None,
overlay_cat=ots,
overlay_scale=scale,
coords_color='red',
coords_lw=3,
lw=1,
save_name=psimage,
no_labels=True,
marker_ra=marker_ra,
marker_dec=marker_dec,
marker_lw=3,
marker_color='cyan',
title=title,
peak=peak,
plot_coords=False,
show_grid=False,
lw_ellipse=3,
ellipse_style=ls,
noisethresh=1.5)
show_overlay(lhdu,
pshdu,
ra,
