Add or remove NaturalEarth features

cf. https://www.naturalearthdata.com/downloads

%matplotlib inline
# for figures in notebook

# import & initialize epygram
import epygram
epygram.init_env()

import os
INPUTS_DIR = os.path.join('..', 'inputs')

r = epygram.formats.resource(os.path.join(INPUTS_DIR, 'pgd.albachir-07km50.SAND.fa'), 'r')

f = r.readfield('SFX.SAND')

# no coastlines, no borders
fig, ax = f.cartoplot(natural_earth_features=[])

import cartopy.feature as cfeature
# no coastlines, no borders
fig, ax = f.cartoplot(natural_earth_features=[])
# coastlines
ax.add_feature(cfeature.NaturalEarthFeature(category='physical',
                                            name='coastline',
                                            scale='10m',
                                            facecolor='none',
                                            edgecolor='k'))
# roads
ax.add_feature(cfeature.NaturalEarthFeature(category='cultural',
                                            name='roads',
                                            scale='10m',
                                            facecolor='none',
                                            edgecolor='r'))

<cartopy.mpl.feature_artist.FeatureArtist at 0x74c6bd76ec50>