def test_center(): ul = some_raster.corner('ul').get_shape(some_raster.crs) br = some_raster.corner('br').get_shape(some_raster.crs) expected_center = Point((ul.x + br.x) / 2, (ul.y + br.y) / 2) expected_center_vector = GeoVector(expected_center, some_raster.crs) assert expected_center_vector.almost_equals(some_raster.center())
def test_almost_equals(): some_crs = CRS({'init': 'epsg:32630'}) another_crs = CRS({'init': 'epsg:32631'}) pt = GeoVector(Point(0, 0), some_crs) similar = GeoVector(Point(0, 0.00000000001), some_crs) assert pt.almost_equals(similar.reproject(another_crs)) different = GeoVector(Point(0, 1), some_crs) assert pt != different