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 zones that consist of a multipolygon, as it
is located at the 180deg/dateline.
"""
grid = LatLonGrid(0.0006)
kamchatka_geom = setup_geom_kamchatka()
kamchatka_geom_tiles = sorted(["GL600U_E342N144T6", "GL600U_E348N144T6", "GL600U_E354N144T6"])
tiles = sorted(grid.search_tiles_in_roi(kamchatka_geom, coverland=False))
assert sorted(tiles) == sorted(kamchatka_geom_tiles)
def test_search_tiles_kamchatka():
"""
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 zones that consist of a multipolygon, as it
is located at 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 = sorted(grid.search_tiles_in_roi(kamchatka_geom, coverland=False))
assert sorted(tiles) == sorted(kamchatka_geom_tiles)