def test_get_intersection(self):
pol1 = self.generate_rect_geojson([0,0,2,1])
pol2 = self.generate_rect_geojson([1,0,3,1])
shp1 = gis_util.geojson_to_shape(pol1)
shp2 = gis_util.geojson_to_shape(pol2)
intersection = gis_util.get_intersection(shp1, shp2)
gis_util.shape_to_wkt(intersection)
def test_polygon_to_multipolygon(self):
poly_geojson = self.generate_rect_geojson([1,-1,2,-2])
poly_shape = gis_util.geojson_to_shape(poly_geojson)
multipoly_shape = gis_util.polygon_to_multipolygon(poly_shape)
self.assertEquals(multipoly_shape.geom_type, 'MultiPolygon')
multipoly_wkt = gis_util.shape_to_wkt(multipoly_shape)
self.assertEquals(multipoly_wkt, 'MULTIPOLYGON (((1 -1, 1 -2, 2 -2, 2 -1, 1 -1)))')
def test_reproject_shape(self):
geojson = self.generate_rect_geojson([1, -1, 2, -2])
shape = gis_util.geojson_to_shape(geojson)
reprojected_shape = gis_util.reproject_shape(shape,
'EPSG:4326',
gis_util.get_mollweide_crs(),
)
rereprojected_shape = gis_util.reproject_shape(reprojected_shape,
gis_util.get_mollweide_crs(),
'EPSG:4326',
)
rereprojected_wkt = gis_util.shape_to_wkt(rereprojected_shape)
def test_wkt_conversion(self):
poly_geojson = self.generate_rect_geojson([1,-1,2,-2])
poly_shape = gis_util.geojson_to_shape(poly_geojson)
poly_wkt = gis_util.shape_to_wkt(poly_shape)
