Bathymetry single-beam surveys of the Caribbean

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

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

Load the CSV formatted data with pandas

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()
caribbean bathymetry

Total running time of the script: (0 minutes 8.261 seconds)

Gallery generated by Sphinx-Gallery