add beaches, places, ports as observation points

[veenstrajelmer]

import osmnx as ox
import dfm_tools as dfmt
import matplotlib.pyplot as plt
plt.close("all")

lon_min, lon_max, lat_min, lat_max = -6, 5, 48, 50.5 # france
lon_min, lon_max, lat_min, lat_max = -2, 2, 48, 49.75 # france smaller
# lon_min, lon_max, lat_min, lat_max = 123, 148, 23, 47 # japan

beaches = ox.features.features_from_bbox(lat_max, lat_min, lon_max, lon_min, tags={"natural":["beach"]})
#cities = ox.features.features_from_bbox(lat_max, lat_min, lon_max, lon_min, tags={"place":["city"]})#,"town","village"]})
harbours = ox.features.features_from_bbox(lat_max, lat_min, lon_max, lon_min, tags={"landuse":["harbour"]})

fig, ax = plt.subplots()
beaches.geometry.plot(ax=ax, color='y')
#cities.geometry.plot(ax=ax, color='r')
harbours.geometry.plot(ax=ax, color='b')
dfmt.plot_coastlines(ax=ax)

Results in:

Can also be done for:

• waterways: https://wiki.openstreetmap.org/wiki/Waterways
• coastlines: https://osmdata.openstreetmap.de/data/
• probably way more

But is slow and warns:
"# UserWarning: This area is 22 times your configured Overpass max query area size. It will automatically be divided up into multiple sub-queries accordingly. This may take a long time."