Bathymetry single-beam surveys of the Caribbean

This dataset is a compilation of several public domain single-beam bathymetry surveys of the ocean in the Caribbean. The data display a wide range of tectonic activity and uneven distribution.

Original source: NOAA NCEI

Pre-processing: Source code for preparation of the original dataset for redistribution in Ensaio

Changes in version 2

In version 1, there were 1,938,095 data taking up a larger area. The data were depth_m and positive downward. Version 2, cropped the data to make it more manageable and converted the depths to bathymetric heights (negative downward).

import pandas as pd
import pygmt

import ensaio

Download and cache the data and return the path to it on disk

fname = ensaio.fetch_caribbean_bathymetry(version=2)
print(fname)

/home/runner/work/_temp/cache/ensaio/v2/caribbean-bathymetry.csv.xz

Load the CSV formatted data with pandas

data = pd.read_csv(fname)
data

	survey_id	latitude	longitude	bathymetry_m
0	86005311	16.09652	-61.52117	-187
1	86005311	16.09415	-61.52104	-177
2	86005311	16.09177	-61.52091	-185
3	86005311	16.08940	-61.52078	-188
4	86005311	16.08703	-61.52066	-192
...	...	...	...	...
294316	JR336	15.28529	-57.01258	-5276
294317	JR336	15.28705	-57.00994	-5277
294318	JR336	15.28883	-57.00732	-5278
294319	JR336	15.29057	-57.00467	-5277
294320	JR336	15.29234	-57.00203	-5276

294321 rows × 4 columns

Make a PyGMT map with the data points colored by the bathymetry.

fig = pygmt.Figure()
pygmt.makecpt(
    cmap="cmocean/topo+h",
    series=[data.bathymetry_m.min(), data.bathymetry_m.max()],
)
fig.plot(
    x=data.longitude,
    y=data.latitude,
    fill=data.bathymetry_m,
    cmap=True,
    style="c0.02c",
    projection="M15c",
    frame=True,
)
fig.colorbar(frame='af+l"bathymetry [m]"')
fig.coast(land="#666666")
fig.show()