GMT grdimage - doesn't create imagine

Q&A
1

Hi everyone,

I have a file that contains 3 columns: longitude, latitude and a height.
I would like to create a plot that is colored based on the height of the point that is contained in the file.
I have used the following commend:
PSFILE=“surface.ps”
PROJ=“-JM0/50/14.2i”
LIMS=“-R-50/50/30/70”

gmt xyz2grd coord_blh.longlat.txt -D+xdegree+ydegree+zGal -Gsoso.nc $LIMS $PROJ -I30s
gmt grd2cpt soso.nc $LIMS -Crainbow -T-0.01/2.1/0.02 -Z > soso.cpt
gmt grdimage soso.nc $LIMS $PROJ -Csoso.cpt -Q -K > $PSFILE
gmt grdcontour soso.nc -C30 $PROJ $LIMS -Wdarkgray -A- -K -O >> $PSFILE
gmt pscoast $PROJ $LIMS -Bx10 -By4 -W -Di -Ir -N1 -K -O -P >> $PSFILE
gmt psscale -D1.5i/1.5i+w2.5i/0.25i+h+jCT+m -Csoso.cpt -B0.5+l"Height Diff (cm)" -O >> $PSFILE

but unfortunately my plot is just grey - as if it doesn’t consider the grdimage command an the colors that should be in the soso.cpt file.
I’m doing something wrong but I don’t know what.
Thank you for the support in advance.

Sorin

coord_blh.longlat.txt (3.5 KB) .

2

No miracles, your data is not interpolated so xyz2grd cannot work. Must use one of the GMT interpolaters first (e.g. surface or triangulate or others) first to create one grid.

3

Hi Joaquim,

Thank you very much. Using both interpolates works - surface and triangulate.

Sorin

1 Like