|
@@ -4,6 +4,7 @@ from grass.gunittest.main import test
|
|
|
from unittest import skip
|
|
|
|
|
|
from grass.pygrass.raster import RasterRow
|
|
|
+from grass.pygrass.raster import raster2numpy
|
|
|
from grass.pygrass.gis.region import Region
|
|
|
|
|
|
class RasterRowRegionTestCase(TestCase):
|
|
@@ -78,5 +79,35 @@ class RasterRowRegionTestCase(TestCase):
|
|
|
|
|
|
rast.close()
|
|
|
|
|
|
+ def test_resampling_to_numpy(self):
|
|
|
+
|
|
|
+ region = Region()
|
|
|
+ region.ewres = 0.1
|
|
|
+ region.nsres = 0.1
|
|
|
+ region.adjust()
|
|
|
+ region.set_raster_region()
|
|
|
+
|
|
|
+ a = raster2numpy(self.name)
|
|
|
+
|
|
|
+ self.assertEqual(len(a), 400)
|
|
|
+
|
|
|
+ region.ewres = 1
|
|
|
+ region.nsres = 1
|
|
|
+ region.adjust()
|
|
|
+ region.set_raster_region()
|
|
|
+
|
|
|
+ a = raster2numpy(self.name)
|
|
|
+
|
|
|
+ self.assertEqual(len(a), 40)
|
|
|
+
|
|
|
+ region.ewres = 5
|
|
|
+ region.nsres = 5
|
|
|
+ region.adjust()
|
|
|
+ region.set_raster_region()
|
|
|
+
|
|
|
+ a = raster2numpy(self.name)
|
|
|
+
|
|
|
+ self.assertEqual(len(a), 8)
|
|
|
+
|
|
|
if __name__ == '__main__':
|
|
|
test()
|