def test_zonal_statistics_raster(self):
"""Test we can do zonal statistics."""
# Same projection
raster = load_test_raster_layer(
'exposure', 'pop_binary_raster_20_20.asc')
raster.keywords['inasafe_default_values'] = {}
vector = load_test_vector_layer(
'aggregation', 'grid_jakarta_4326.geojson')
vector.keywords['hazard_keywords'] = {}
vector.keywords['aggregation_keywords'] = {}
number_fields = vector.fields().count()
vector = zonal_stats(raster, vector)
self.assertEqual(vector.fields().count(), number_fields + 1)
self.assertEqual(vector.geometryType(), QgsWkbTypes.PolygonGeometry)
# With different projections
raster = load_test_raster_layer(
'exposure', 'pop_binary_raster_20_20.asc')
raster.keywords['inasafe_default_values'] = {}
vector_b = load_test_vector_layer(
'aggregation', 'grid_jakarta_4326.geojson')
vector_b.keywords['hazard_keywords'] = {}
vector_b.keywords['aggregation_keywords'] = {}
number_fields = vector_b.fields().count()
vector_b = reproject(vector, QgsCoordinateReferenceSystem(3857))
vector_b = zonal_stats(raster, vector_b)
self.assertEqual(vector_b.fields().count(), number_fields + 1)
self.assertEqual(vector_b.geometryType(), QgsWkbTypes.PolygonGeometry)
# We compare the results between these 2 zonal stats.
for feature_a, feature_b in zip(
vector.getFeatures(), vector_b.getFeatures()):
self.assertEqual(feature_a.attributes(), feature_b.attributes())
def test_zonal_statistics_raster(self):
"""Test we can do zonal statistics."""
raster = load_test_raster_layer('exposure',
'pop_binary_raster_20_20.asc')
raster.keywords['inasafe_default_values'] = {}
vector = load_test_vector_layer('aggregation',
'grid_jakarta_4326.geojson')
vector.keywords['hazard_keywords'] = {}
vector.keywords['aggregation_keywords'] = {}
number_fields = vector.fields().count()
vector = zonal_stats(raster, vector)
self.assertEqual(vector.fields().count(), number_fields + 1)
self.assertEqual(vector.geometryType(), QGis.Polygon)
def test_zonal_statistics_raster(self):
"""Test we can do zonal statistics."""
raster = load_test_raster_layer(
'exposure', 'pop_binary_raster_20_20.asc')
raster.keywords['inasafe_default_values'] = {}
vector = load_test_vector_layer(
'aggregation', 'grid_jakarta_4326.geojson')
vector.keywords['hazard_keywords'] = {}
vector.keywords['aggregation_keywords'] = {}
number_fields = vector.fields().count()
vector = zonal_stats(raster, vector)
self.assertEqual(vector.fields().count(), number_fields + 1)
self.assertEqual(vector.geometryType(), QGis.Polygon)