3D interpolation via gmt surface

Q&A
17

OK, so that’s what I said here (or in another similar thread). Not knowing why the greensplineis failing (it with many points it’s a very demanding method) the alternatives are.

  • With plain GMT: create individual grids, one by each depth layer, and stack them together with grdinterpolate
  • Use the GMT.jl Julia wrapper. With it it would be (hopefully) just do something like Cube = xyzw2cube("/path/to/yourfile"); Next, either save it with gmtwrite(Cube, "cube.nc") and go back to plain GMT. Or use G = cubeslice(Cube, newdepth) and get that desired (linear) interpolated layer.
1 Like