How to use the pytileproj.geometry.setup_geom_kamchatka function in pytileproj

To help you get started, we’ve selected a few pytileproj examples, based on popular ways it is used in public projects.

Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.

github TUW-GEO / Equi7Grid / tests / test_equi7grid.py View on Github external
def test_search_tiles_kamchatka(self):
        """
        Tests the tile searching over Kamchatka in far east Sibiria;
    
        This test is especially nice, as it contains also a tile that covers both,
        the ROI and the continental zone, but the intersection of the tile and
        the ROI is outside of the zone.
    
        Furthermore, it also covers Equi7Grid subgrids that consist of a multipolygon,
        as they overspan the 180deg/dateline.
        """
    
        grid = Equi7Grid(500)
    
        kamchatka_geom = setup_geom_kamchatka()
        kamchatka_geom_tiles = sorted(['AS500M_E072N078T6', 'AS500M_E078N078T6',
                                       'AS500M_E078N084T6', 'NA500M_E036N078T6',
                                       'NA500M_E036N084T6', 'NA500M_E042N078T6',
                                       'NA500M_E042N084T6'])

        tiles = grid.search_tiles_in_roi(kamchatka_geom, coverland=False)
    
        assert sorted(tiles) == sorted(kamchatka_geom_tiles)