29 """!Construct the indexer object 31 @param[in] depth depth of the hierarchy to construct 33 self.
htm = esutil.htm.HTM(depth)
36 """!Get all shards that touch a circular aperture 38 @param[in] ctrCoord afwGeom.SpherePoint ICRS center of the aperture 39 @param[in] radius afwGeom.Angle object of the aperture radius 40 @param[out] A pipeBase.Struct with the list of shards, shards, and a boolean arry, boundary_mask, 41 indicating whether the shard touches the boundary (True) or is fully contained (False). 43 pixel_id_list = self.
htm.intersect(ctrCoord.getLongitude().asDegrees(),
44 ctrCoord.getLatitude().asDegrees(),
45 radius.asDegrees(), inclusive=
True)
46 covered_pixel_id_list = self.
htm.intersect(ctrCoord.getLongitude().asDegrees(),
47 ctrCoord.getLatitude().asDegrees(),
48 radius.asDegrees(), inclusive=
False)
49 is_on_boundary = (pixel_id
not in covered_pixel_id_list
for pixel_id
in pixel_id_list)
50 return pixel_id_list, is_on_boundary
53 """!Generate trixel ids for each row in an input file 55 @param[in] ra_list List of RA coordinate in degrees 56 @param[in] dec_list List of Dec coordinate in degrees 57 @param[out] A list of pixel ids 59 return self.
htm.lookup_id(ra_list, dec_list)
63 """!Make a data id. Meant to be overridden. 64 @param[in] pixel_id An identifier for the pixel in question. 65 @param[in] dataset_name Name of the dataset to use. 66 @param[out] dataId (dictionary) 71 return {
'pixel_id': pixel_id,
'name': dataset_name}
def get_pixel_ids(self, ctrCoord, radius)
Get all shards that touch a circular aperture.
def make_data_id(pixel_id, dataset_name)
Make a data id.
def __init__(self, depth=8)
Construct the indexer object.
def index_points(self, ra_list, dec_list)
Generate trixel ids for each row in an input file.