def test_geometry(self): expected = ( utils.Extent((136700, 455795, 136705, 455800), utils.get_sr("EPSG:28992")) .as_geometry() .ExportToWkt() ) self.assertEqual(expected, self.source.geometry.ExportToWkt())
def test_extent(self): expected = ( utils.Extent((136700, 455795, 136705, 455800), utils.get_sr("EPSG:28992")) .transformed(utils.get_sr("EPSG:4326")) .bbox ) assert_allclose(self.source.extent, expected, atol=1e-10)
def test_extent_has_repr(self): sr = utils.get_sr("EPSG:4326") extent = utils.Extent(sr=sr, bbox=(0, 0, 1, 1)) self.assertTrue(repr(extent))
def test_extent(self): sr = utils.get_sr("EPSG:4326") extent = utils.Extent(sr=sr, bbox=(0, 0, 1, 1)) geometry = extent.as_geometry() self.assertEqual(str(geometry), "POLYGON ((0 0,1 0,1 1,0 1,0 0))") self.assertEqual(str(geometry.GetSpatialReference()), str(sr))
def _get_extent(self): bbox = self.geo_transform.get_bbox( (0, 0), (self.gdal_dataset.RasterYSize, self.gdal_dataset.RasterXSize) ) return utils.Extent(bbox, utils.get_sr(self.projection))
def _get_extent(self): if not self.data.size: return bbox = self.geo_transform.get_bbox((0, 0), self.data.shape[1:]) return utils.Extent(bbox, utils.get_sr(self.projection))