Skip to content

Instantly share code, notes, and snippets.

@epifanio
Created June 27, 2017 23:08
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save epifanio/5376f86a16daf46077d54d8938e25916 to your computer and use it in GitHub Desktop.
Save epifanio/5376f86a16daf46077d54d8938e25916 to your computer and use it in GitHub Desktop.
gdal_translate -of GMT grassdata/nc_basic_spm_grass7/PERMANENT/cellhd/elevation elevation.gmt
gdalwarp -overwrite elevation.gmt out_elevation.gmt -t_srs EPSG:4326
gmt grdgradient out_elevation.gmt -Ne0.8 -A100 -fg -Gus_i.nc
gmt makecpt -Chaxby -T50/160/10 -Z > topo.cpt
gmt grdimage out_elevation.gmt -JM6i -P -Ba -Ctopo.cpt -V -K -Ius_i.nc >> GMT_tut_12.ps
gmt grdcontour out_elevation.gmt -JM6i -O -BWSne -C20 -A20 >> GMT_tut_12.ps
gmt psscale -DjTC+w5i/0.25i+h+o0/-1i -Rout_elevation.gmt -JM6i -Ctopo.cpt -B+tcontinuous -I0.5 -O -Ba1000:Elevation:/:m: >> GMT_tut_12.ps
@epifanio
Copy link
Author

GMT commands for osgeolive quickstart

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment