def superimpose_bound_vtx(self, vertices, label, color='red',
coord_in='C'):
'''Superimpose the footprint of a survey on the background image
by giving the ra/dec corners of the image. The enclosed survey
footprint is generated by calling healpy.query_polygon.
Parameters
----------
vertices : array-like with shape (n,2)
The n corners of the survey footprint in degrees.
label : string
The label to put on the colorbar for this survey
color : string or array-like with shape (3,)
The color to use when overlaying the survey footprint. Either a
string or rgb triplet. Default : 'red'
coord_in : 'C', 'E', or 'G'
The coordinate system of the input vertices. Default : 'C'
'''
hpx_map = util.gen_map_polygon(vertices, self.nside)
self.superimpose_hpxmap(hpx_map, label, color=color,
coord_in=coord_in)
def superimpose_bound_vtx(self,
vertices,
label,
color='red',
coord_in='C'):
'''Superimpose the footprint of a survey on the background image
by giving the ra/dec corners of the image. The enclosed survey
footprint is generated by calling healpy.query_polygon.
Parameters
----------
vertices : array-like with shape (n,2)
The n corners of the survey footprint in degrees.
label : string
The label to put on the colorbar for this survey
color : string or array-like with shape (3,)
The color to use when overlaying the survey footprint. Either a
string or rgb triplet. Default : 'red'
coord_in : 'C', 'E', or 'G'
The coordinate system of the input vertices. Default : 'C'
'''
hpx_map = util.gen_map_polygon(vertices, self.nside)
self.superimpose_hpxmap(hpx_map, label, color=color, coord_in=coord_in)
def get_polygon(self, survey_name):
'''Generates a healpix map for a given survey given vertices of
a polygon.
Handler in .cfg file: "polygon"
Parameters
----------
survey_name : string
Configuration file block specifying survey parameters
Returns
-------
hpx_map : array-like
The healpix map associated with the survey
Notes
-----
The configuration file must contain 'vertex1', ..., 'vertexXX'
entries containing the ra,dec locations of each vertex. Nside of
the output map is the same as the nside specified in the
initialization of this class.
'''
lons = []
lats = []
i = 1
while True:
radec_point = 'vertex'+str(i)
try:
radec_val = self.config.get(survey_name, radec_point)
except ConfigParser.NoOptionError:
break
lonlat_val = radec_val.split(',')
lonlat_val = [tmp.strip() for tmp in lonlat_val]
lonlat_val = SkyCoord(lonlat_val[0], lonlat_val[1])
lons.append(lonlat_val.ra.deg)
lats.append(lonlat_val.dec.deg)
i += 1
vtxs = np.transpose([lons, lats])
coord = self.config.get(survey_name, 'coord')
hpx_map = util.gen_map_polygon(vtxs, self.nside)
return [hpx_map], coord
