def test_puget_sound():
""" lots of islands """
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'f', *puget_sound)
gj = format_for_geojson(land, map_bounds)
open("puget_sound.geojson", 'w').write(gj)
def test_wa_coast():
""" simple straight shoreline """
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'i', *wa_coast)
gj = format_for_geojson(land, map_bounds)
open("wa_shore.geojson", 'w').write(gj)
def test_long_island():
""" a barrier island -- 3 polys """
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'i', *li)
gj = format_for_geojson(land)
open("long_island_shore.geojson", 'w').write(gj)
def test_south_florida():
""" Lake Okeechobee is in there! """
# wl, el, nl, sl,
south_florida = (-83.0, -79.5, 27.5, 25.0)
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'h', *south_florida)
gj = format_for_geojson(land, map_bounds)
open("south_florida.geojson", 'w').write(gj)
def test_lake_st_clair():
# wl, el, nl, sl,
bounds = (-83.166, -82.255, 42.849, 42.226)
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'f', *bounds)
gj = format_for_geojson(land,
# map_bounds=map_bounds
)
open("lake_st_clair.geojson", 'w').write(gj)
def test_west_lake_erie():
# wl, el, nl, sl,
bounds = (-83.859, -81.343, 42.255, 41.189)
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'f', *bounds)
gj = format_for_geojson(land,
# map_bounds=map_bounds
)
open("lake_erie.geojson", 'w').write(gj)
def test_great_lakes2():
land, map_bounds = load_and_clip_gshhs(GSHHS_dir, 'f', *gl2)
gj = format_for_geojson(land,
# map_bounds=map_bounds
)
open("great_lakes_2.geojson", 'w').write(gj)