Python proj_is_sameの例

コード例 #1

ファイル: test_gis.py プロジェクト: shawndegroot/salem

    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))

コード例 #2

ファイル: test_gis.py プロジェクト: fmaussion/salem

    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))