3D seismic velocity model with topography

Akroma · January 31, 2020, 10:44pm

Dear GMT users,

I want to plot 3d model like in the example figure. I have got separate columns for longitude-latitude-value-depth

What is the proper approach and which gmt commands should i use in order to get this model ?

Thanks in advance

example

Joaquim · February 3, 2020, 8:56pm

I’m afraid this is not doable in GMT.

Luojiahong · January 31, 2024, 9:18am

here is my solution：

github.com

Luojiahong/PyGMT_3D_mapview/blob/main/mapview_3d.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd \n",
    "import xarray as xr \n",
    "import pygmt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {

This file has been truncated. show original

Akroma · January 31, 2024, 3:10pm

@Luojiahong

It looks beautiful. After almost four years I wasn’t expecting an answer. Thank you so much : )