Plot terrain and bathymetry

[blaylockbk]

I have always wanted to implement a way to plot terrain and bathymetry data.

This tweet reminded I still want to do this: https://twitter.com/DanJonesOcean/status/1448597935508369411?s=20

I could implement the example from here: https://github.com/mattphysics/plot_bathymetry/blob/master/plot_bathymetry.ipynb

Any better methods and data sources for adding Bathymetry? Could get data from:

• https://ngdc.noaa.gov/mgg/global/
• https://topex.ucsd.edu/WWW_html/mar_topo.html
• https://www.ngdc.noaa.gov/mgg/topo/globe.html
• https://www.gebco.net/data_and_products/gridded_bathymetry_data/
• http://rabernat.github.io/intro_to_physical_oceanography/01_ocean_bathymetry.html
• http://iridl.ldeo.columbia.edu/SOURCES/.NOAA/.NGDC/

[blaylockbk]

from toolbox.cartopy_tools import common_features
ax = common_features('50m', dpi=200).BATHYMETRY(10, kind='contourf').TERRAIN(10, kind='contourf').ax

ax = common_features('50m', dpi=200).BATHYMETRY(10, kind='contourf').TERRAIN(10, kind='contourf', cmap='Greys').ax

# Full resolution
# This took a long time to run (>2 mins)
ax = common_features('10m', dpi=200).BATHYMETRY(1).TERRAIN(1, cmap='Greys').ax
ax.center_extent(state="CA", pad=5)

(Issue: notice that some areas that are "land" below 0 are shaded blue. Unfortunately the ETOPO1 dataset does not have a land mask.)

[blaylockbk]

The performance could be improved. This is not very good to use when you need to transform the data to a projection that is not PlateCarree (especially if you use fine resolution).

Perhaps I could do a set_extent method, get the map limits, trim off the data not used in the frame, and then plot TERRAIN and BATHYMETRY.