def test_clip_geometry(self):
"""Test that we can clip a geometry using another geometry."""
geometry = QgsGeometry.fromPolyline([
QgsPoint(10, 10),
QgsPoint(20, 20),
QgsPoint(30, 30),
QgsPoint(40, 40)
])
clip_polygon = QgsGeometry.fromPolygon([[
QgsPoint(20, 20),
QgsPoint(20, 30),
QgsPoint(30, 30),
QgsPoint(30, 20),
QgsPoint(20, 20)
]])
result = clip_geometry(clip_polygon, geometry)
expected_wkt = 'LINESTRING(20.0 20.0, 30.0 30.0)'
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
# Now poly on poly clip test
clip_polygon = QgsGeometry.fromWkt(
'POLYGON((106.8218 -6.1842,106.8232 -6.1842,'
'106.82304963 -6.18317148,106.82179736 -6.18314774,'
'106.8218 -6.1842))')
geometry = QgsGeometry.fromWkt(
'POLYGON((106.8216869 -6.1852067,106.8213233 -6.1843051,'
'106.8212891 -6.1835559,106.8222566 -6.1835184,'
'106.8227557 -6.1835076,106.8228588 -6.1851204,'
'106.8216869 -6.1852067))')
result = clip_geometry(clip_polygon, geometry)
expected_wkt = (
'POLYGON((106.82179833 -6.18353616,106.8222566 -6.1835184,'
'106.8227557 -6.1835076,106.82279996 -6.1842,'
'106.8218 -6.1842,106.82179833 -6.18353616))')
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
# Now point on poly test clip
geometry = QgsGeometry.fromWkt('POINT(106.82241 -6.18369)')
result = clip_geometry(clip_polygon, geometry)
if qgis_version() > 10800:
expected_wkt = 'POINT(106.82241 -6.18369)'
else:
expected_wkt = 'POINT(106.822410 -6.183690)'
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
def test_clip_geometry(self):
"""Test that we can clip a geometry using another geometry."""
geometry = QgsGeometry.fromPolyline([
QgsPoint(10, 10), QgsPoint(20, 20), QgsPoint(30, 30),
QgsPoint(40, 40)])
clip_polygon = QgsGeometry.fromPolygon([
[QgsPoint(20, 20),
QgsPoint(20, 30),
QgsPoint(30, 30),
QgsPoint(30, 20),
QgsPoint(20, 20)]])
result = clip_geometry(clip_polygon, geometry)
expected_wkt = 'LINESTRING(20.0 20.0, 30.0 30.0)'
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
# Now poly on poly clip test
clip_polygon = QgsGeometry.fromWkt(
'POLYGON((106.8218 -6.1842,106.8232 -6.1842,'
'106.82304963 -6.18317148,106.82179736 -6.18314774,'
'106.8218 -6.1842))')
geometry = QgsGeometry.fromWkt(
'POLYGON((106.8216869 -6.1852067,106.8213233 -6.1843051,'
'106.8212891 -6.1835559,106.8222566 -6.1835184,'
'106.8227557 -6.1835076,106.8228588 -6.1851204,'
'106.8216869 -6.1852067))')
result = clip_geometry(clip_polygon, geometry)
expected_wkt = (
'POLYGON((106.82179833 -6.18353616,106.8222566 -6.1835184,'
'106.8227557 -6.1835076,106.82279996 -6.1842,'
'106.8218 -6.1842,106.82179833 -6.18353616))')
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
# Now point on poly test clip
geometry = QgsGeometry.fromWkt('POINT(106.82241 -6.18369)')
result = clip_geometry(clip_polygon, geometry)
if qgis_version() > 10800:
expected_wkt = 'POINT(106.82241 -6.18369)'
else:
expected_wkt = 'POINT(106.822410 -6.183690)'
# There should only be one feature that intersects this clip
# poly so this assertion should work.
assert compare_wkt(expected_wkt, str(result.exportToWkt()))
def test_minimum_extent(self):
"""Test we can compute the minimum extent in the IF."""
# Without aggregation layer
hazard_layer = load_test_vector_layer('hazard',
'flood_multipart_polygons.shp')
exposure_layer = load_test_vector_layer('exposure', 'roads.shp')
impact_function = ImpactFunction()
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent without an aggregation layer is '
'failing.')
self.assertTrue(
compare_wkt(
'Polygon (('
'106.8080099999999959 -6.19531000000000009, '
'106.8080099999999959 -6.16752599999999962, '
'106.83456946836641066 -6.16752599999999962, '
'106.83456946836641066 -6.19531000000000009, '
'106.8080099999999959 -6.19531000000000009))',
impact_function.analysis_extent.exportToWkt()), message)
# Without aggregation layer but with a requested_extent
hazard_layer = load_test_vector_layer('hazard',
'flood_multipart_polygons.shp')
exposure_layer = load_test_vector_layer('exposure', 'roads.shp')
impact_function = ImpactFunction()
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
impact_function.requested_extent = wkt_to_rectangle(
'POLYGON (('
'106.772279 -6.237576, '
'106.772279 -6.165415, '
'106.885165 -6.165415, '
'106.885165 -6.237576, '
'106.772279 -6.237576'
'))')
impact_function.requested_extent_crs = QgsCoordinateReferenceSystem(
4326)
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent without an aggregation layer but '
'with a requested extent is failing.')
self.assertTrue(
compare_wkt(
'Polygon (('
'106.8080099999999959 -6.19531000000000009, '
'106.8080099999999959 -6.16752599999999962, '
'106.83456946836641066 -6.16752599999999962, '
'106.83456946836641066 -6.19531000000000009, '
'106.8080099999999959 -6.19531000000000009))',
impact_function.analysis_extent.exportToWkt()), message)
# With an aggregation layer, without selection
hazard_layer = load_test_vector_layer('gisv4', 'hazard',
'classified_vector.geojson')
exposure_layer = load_test_vector_layer('gisv4', 'exposure',
'building-points.geojson')
aggregation_layer = load_test_vector_layer('gisv4', 'aggregation',
'small_grid.geojson')
impact_function = ImpactFunction()
impact_function.aggregation = aggregation_layer
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
impact_function.use_selected_features_only = False
impact_function.aggregation.select(0)
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent with an aggregation layer is '
'failing.')
self.assertTrue(
compare_wkt(
'Polygon ((106.9033179652593617 -6.18324454090033182, '
'106.90331796525939012 -6.2725478115989306, '
'106.72365490843547775 -6.2725478115989306, '
'106.72365490843547775 -6.18324645462287137, '
'106.72365490843547775 -6.09392810187095257, '
'106.81348643684744104 -6.09392810187095257, '
'106.9033179652593617 -6.09392810187095257, '
'106.9033179652593617 -6.18324454090033182))',
impact_function.analysis_extent.exportToWkt()), message)
# With an aggregation layer, with selection
impact_function.use_selected_features_only = True
impact_function.aggregation = aggregation_layer
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent with an aggregation layer and '
'a selection is failing.')
self.assertTrue(
compare_wkt(
'Polygon ((106.72365490843547775 -6.09392810187095257, '
'106.81348643684744104 -6.09392810187095257, '
'106.81348643684744104 -6.18324645462287137, '
'106.72365490843547775 -6.18324645462287137, '
'106.72365490843547775 -6.09392810187095257))',
impact_function.analysis_extent.exportToWkt()), message)
def test_minimum_extent(self):
"""Test we can compute the minimum extent in the IF."""
# Without aggregation layer
hazard_layer = load_test_vector_layer(
'hazard', 'flood_multipart_polygons.shp')
exposure_layer = load_test_vector_layer('exposure', 'roads.shp')
impact_function = ImpactFunction()
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent without an aggregation layer is '
'failing.')
self.assertTrue(
compare_wkt(
'Polygon (('
'106.8080099999999959 -6.19531000000000009, '
'106.8080099999999959 -6.16752599999999962, '
'106.83456946836641066 -6.16752599999999962, '
'106.83456946836641066 -6.19531000000000009, '
'106.8080099999999959 -6.19531000000000009))',
impact_function.analysis_extent.exportToWkt()),
message
)
# Without aggregation layer but with a requested_extent
hazard_layer = load_test_vector_layer(
'hazard', 'flood_multipart_polygons.shp')
exposure_layer = load_test_vector_layer('exposure', 'roads.shp')
impact_function = ImpactFunction()
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
impact_function.requested_extent = wkt_to_rectangle(
'POLYGON (('
'106.772279 -6.237576, '
'106.772279 -6.165415, '
'106.885165 -6.165415, '
'106.885165 -6.237576, '
'106.772279 -6.237576'
'))')
impact_function.requested_extent_crs = QgsCoordinateReferenceSystem(
4326)
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent without an aggregation layer but '
'with a requested extent is failing.')
self.assertTrue(
compare_wkt(
'Polygon (('
'106.8080099999999959 -6.19531000000000009, '
'106.8080099999999959 -6.16752599999999962, '
'106.83456946836641066 -6.16752599999999962, '
'106.83456946836641066 -6.19531000000000009, '
'106.8080099999999959 -6.19531000000000009))',
impact_function.analysis_extent.exportToWkt()),
message
)
# With an aggregation layer, without selection
hazard_layer = load_test_vector_layer(
'gisv4', 'hazard', 'classified_vector.geojson')
exposure_layer = load_test_vector_layer(
'gisv4', 'exposure', 'building-points.geojson')
aggregation_layer = load_test_vector_layer(
'gisv4', 'aggregation', 'small_grid.geojson')
impact_function = ImpactFunction()
impact_function.aggregation = aggregation_layer
impact_function.exposure = exposure_layer
impact_function.hazard = hazard_layer
impact_function.use_selected_features_only = False
impact_function.aggregation.select(0)
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent with an aggregation layer is '
'failing.')
self.assertTrue(
compare_wkt(
'Polygon ((106.9033179652593617 -6.18324454090033182, '
'106.90331796525939012 -6.2725478115989306, '
'106.72365490843547775 -6.2725478115989306, '
'106.72365490843547775 -6.18324645462287137, '
'106.72365490843547775 -6.09392810187095257, '
'106.81348643684744104 -6.09392810187095257, '
'106.9033179652593617 -6.09392810187095257, '
'106.9033179652593617 -6.18324454090033182))',
impact_function.analysis_extent.exportToWkt()),
message
)
# With an aggregation layer, with selection
impact_function.use_selected_features_only = True
impact_function.aggregation = aggregation_layer
status, message = impact_function.prepare()
self.assertEqual(PREPARE_SUCCESS, status, message)
message = (
'Test about the minimum extent with an aggregation layer and '
'a selection is failing.')
self.assertTrue(
compare_wkt(
'Polygon ((106.72365490843547775 -6.09392810187095257, '
'106.81348643684744104 -6.09392810187095257, '
'106.81348643684744104 -6.18324645462287137, '
'106.72365490843547775 -6.18324645462287137, '
'106.72365490843547775 -6.09392810187095257))',
impact_function.analysis_extent.exportToWkt()),
message
)