GitLab

The script below uses the ocean class to load and interpolate gebco bathy data. The results look very weird:

from kadlu.geospatial.ocean import Ocean
import numpy as np
import matplotlib
import matplotlib.pyplot as plt

kwargs = dict(south=43.2, north=46.2, west=-61, east=-58)
o = Ocean(load_bathymetry='gebco', **kwargs)
lat = np.linspace(43.2, 46.2, 250)
lon = np.linspace(-61,-58,250)
bathy = o.bathy(lat=lat, lon=lon, grid=True)
fig, ax = plt.subplots()
img = ax.contourf(lon, lat, bathy, 100, cmap=matplotlib.cm.get_cmap('viridis_r'))
fig.colorbar(img, ax=ax, format='%2.0f')# colobar
plt.show()