def plot_aitoff_wcs_outline(wcsfn, plotfn, W=256, zoom=True):
#anutil.log_init(3)
H = int(W // 2)
# Create Hammer-Aitoff WCS of the appropriate size.
wcs = anutil.anwcs_create_allsky_hammer_aitoff(0., 0., W, H)
plot = ps.Plotstuff(outformat='png', size=(W, H))
plot.wcs = wcs
plot.linestep = 1.
plot.color = 'verydarkblue'
plot.apply_settings()
plot.line_constant_ra(180, -90, 90)
plot.line_constant_ra(-180, 90, -90)
plot.fill()
#plot.plot_grid(60, 30, 60, 30)
plot.fontsize = 12
ras = [-180, -120, -60, 0, 60, 120, 180]
decs = [-60, -30, 0, 30, 60]
# dark gray
plot.rgb = (0.3, 0.3, 0.3)
plot.apply_settings()
for ra in ras:
plot.line_constant_ra(ra, -90, 90)
plot.stroke()
for dec in decs:
plot.line_constant_dec(dec, -180, 180)
plot.stroke()
plot.color = 'gray'
plot.apply_settings()
for ra in ras:
plot.move_to_radec(ra, 0)
plot.text_radec(ra, 0, '%i' % ((ra + 360) % 360))
plot.stroke()
for dec in decs:
if dec != 0:
plot.move_to_radec(0, dec)
plot.text_radec(0, dec, '%+i' % dec)
plot.stroke()
plot.color = 'white'
plot.lw = 3
out = plot.outline
#out.fill = 1
out.wcs_file = wcsfn
anutil.anwcs_print_stdout(out.wcs)
plot.plot('outline')
# Not helpful to add constellations in this view
#ann = plot.annotations
#ann.NGC = ann.bright = ann.HD = 0
#ann.constellations = 1
#plot.plot('annotations')
if zoom:
owcs = anutil.Tan(wcsfn, 0)
# MAGIC 15 degrees radius
#if owcs.radius() < 15.:
if True:
ra, dec = owcs.radec_center()
# MAGIC 36-degree width zoom-in
# MAGIC width, height are arbitrary
zoomwcs = anutil.anwcs_create_box(ra, dec, 36, 1000, 1000)
out.wcs = zoomwcs
#plot.color = 'gray'
plot.lw = 1
plot.dashed(3)
plot.plot('outline')
plot.write(plotfn)
def debug(self):
from astrometry.util.util import anwcs_print_stdout
print 'WcslibWcs:'
anwcs_print_stdout(self.wcs)
def debug(self):
from astrometry.util.util import anwcs_print_stdout
print('WcslibWcs:')
anwcs_print_stdout(self.wcs)
def plot_aitoff_wcs_outline(wcsfn, plotfn, W=256, zoom=True):
#anutil.log_init(3)
H = W/2
# Create Hammer-Aitoff WCS of the appropriate size.
wcs = anutil.anwcs_create_allsky_hammer_aitoff(0., 0., W, H)
plot = ps.Plotstuff(outformat='png', size=(W, H))
plot.wcs = wcs
plot.linestep = 1.
plot.color = 'verydarkblue'
plot.apply_settings()
plot.line_constant_ra(180, -90, 90)
plot.line_constant_ra(-180, 90, -90)
plot.fill()
#plot.plot_grid(60, 30, 60, 30)
plot.fontsize = 12
ras = [-180, -120, -60, 0, 60, 120, 180]
decs = [-60, -30, 0, 30, 60]
# dark gray
plot.rgb = (0.3,0.3,0.3)
plot.apply_settings()
for ra in ras:
plot.line_constant_ra(ra, -90, 90)
plot.stroke()
for dec in decs:
plot.line_constant_dec(dec, -180, 180)
plot.stroke()
plot.color = 'gray'
plot.apply_settings()
for ra in ras:
plot.move_to_radec(ra, 0)
plot.text_radec(ra, 0, '%i'%((ra+360)%360))
plot.stroke()
for dec in decs:
if dec != 0:
plot.move_to_radec(0, dec)
plot.text_radec(0, dec, '%+i'%dec)
plot.stroke()
plot.color = 'white'
plot.lw = 3
out = plot.outline
#out.fill = 1
out.wcs_file = wcsfn
anutil.anwcs_print_stdout(out.wcs)
plot.plot('outline')
# Not helpful to add constellations in this view
#ann = plot.annotations
#ann.NGC = ann.bright = ann.HD = 0
#ann.constellations = 1
#plot.plot('annotations')
if zoom:
owcs = anutil.Tan(wcsfn, 0)
# MAGIC 15 degrees radius
#if owcs.radius() < 15.:
if True:
ra,dec = owcs.radec_center()
# MAGIC 36-degree width zoom-in
# MAGIC width, height are arbitrary
zoomwcs = anutil.anwcs_create_box(ra, dec, 36, 1000,1000)
out.wcs = zoomwcs
#plot.color = 'gray'
plot.lw = 1
plot.dashed(3)
plot.plot('outline')
plot.write(plotfn)
