def test_same_proj(self):
# this should work regardless of gdal or not:
p1 = pyproj.Proj('+proj=utm +zone=15 +datum=NAD83 '
'+ellps=GRS80 +towgs84=0,0,0 +units=m +no_defs')
p2 = pyproj.Proj('+proj=utm +zone=15 +datum=NAD83 +units=m +no_defs '
'+ellps=GRS80 +towgs84=0,0,0')
self.assertTrue(gis.proj_is_same(p1, p2))
# this needs gdal
p1 = gis.check_crs('+init=epsg:26915')
p2 = pyproj.Proj('+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 '
'+units=m +no_defs')
if gis.has_gdal:
self.assertTrue(gis.proj_is_same(p1, p2))
def test_same_proj(self):
# this should work regardless of gdal or not:
p1 = pyproj.Proj('+proj=utm +zone=15 +datum=NAD83 '
'+ellps=GRS80 +towgs84=0,0,0 +units=m +no_defs')
p2 = pyproj.Proj('+proj=utm +zone=15 +datum=NAD83 +units=m +no_defs '
'+ellps=GRS80 +towgs84=0,0,0')
self.assertFalse(p1.srs == p2.srs)
self.assertTrue(gis.proj_is_same(p1, p2))
# this needs gdal
p1 = pyproj.Proj(init='epsg:26915')
p2 = pyproj.Proj( '+proj=utm +zone=15 +datum=NAD83 +units=m +no_defs '
'+ellps=GRS80 +towgs84=0,0,0')
if gis.has_gdal:
self.assertTrue(gis.proj_is_same(p1, p2))