def __getitem__(self, proj_point): """Get value for given projection point""" # XXX why isn't this next line using the self._point_transformer matrix # that's areadly been created? tm = geom.get_transform_matrix(self, self._spec) # blech (x, y) = proj_point (tx, ty) = geom.transform_coord(proj_point, tm, src="map", dst="image") return self._data[tx, ty]
def setUp(self): params = { 'band': 'LANDSAT_8/OLI_TIRS/sr_band2', 'x': -1850865, 'y': 2956785, 't1': '2013-01-01', 't2': '2015-01-01' } # Load the saved data that was originally returned by running the above # query against the LCMAP REST service. data = json.load(open("tests/data/tiles-result-payload.json")) fake_resp = FakeLCMAPRESTResponse(data) fake_http = FakeLCMAPHTTP(fake_resp) sr = surface_reflectance.SurfaceReflectance(fake_http) (self.spec, self.tiles) = sr.tiles(**params) self.tile = self.tiles[0] self.xform_matrix = geom.get_transform_matrix(self.tile, self.spec)
def __init__(self, tile, spec, mask=True, shape=True, unscale=True): self._tile = tile self._spec = spec self._data = decode(spec, tile, mask, shape, unscale) self._point_transformer = geom.get_transform_matrix(self, spec) pass