Topography of the 2018 lava flows of the Sierra Negra volcano, Ecuador

Topography of the 2018 lava flows of the Sierra Negra volcano, Ecuador#

This is a structure-from-motion point cloud of the 2018 lava flows of the Sierra Negra volcano, located on the Galápagos islands, Ecuador. The survey covers a small region at the flank of the volcano and shows different structures and terrain roughness on the lava flows.

Original source: Carr, B. (2020). Sierra Negra Volcano (TIR Flight 3): Galápagos, Ecuador, October 22 2018. Distributed by OpenTopography.

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

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/v1/sierra-negra-topography.csv.xz

Load the CSV formatted data with pandas

longitude latitude elevation_m
0 -91.115651 -0.783062 930.1
1 -91.115658 -0.783056 930.7
2 -91.115649 -0.783060 930.3
3 -91.115656 -0.783063 929.7
4 -91.115655 -0.783068 929.2
... ... ... ...
1731379 -91.118421 -0.781943 990.7
1731380 -91.118303 -0.781933 990.2
1731381 -91.118357 -0.781971 992.4
1731382 -91.118354 -0.781940 991.2
1731383 -91.118374 -0.781945 991.4

1731384 rows × 3 columns



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

fig = pygmt.Figure()
fig.basemap(
    region=[
        data.longitude.min(),
        data.longitude.max(),
        data.latitude.min(),
        data.latitude.max(),
    ],
    projection="M15c",
    frame=True,
)
pygmt.makecpt(cmap="viridis", series=[data.elevation_m.min(), data.elevation_m.max()])
fig.plot(
    x=data.longitude, y=data.latitude, fill=data.elevation_m, cmap=True, style="c0.01c"
)
fig.colorbar(frame='af+l"elevation [m]"')
fig.show()
sierra negra topography

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

Gallery generated by Sphinx-Gallery