GMT

Running grdfilter with -I < input increment produces odd behaviour. Increasing the spacing works fine in the x direction, but the y direction keeps the appearance of original input spacing. Upon closer inspection, the new points have been created, but not interpolated

I can achieve the desired result ( point –> Gaussian blob) if I run grdfilter rotate the grid 90 degrees, and run it again, but that seems wrong that I should have to do this manually. Am I missing something?

The following script generates a single point, and shows the result of grdfilter.

#!/bin/bash

echo "0.5 0.5 1" | xyz2grd -N0 -I0.1 -An -Gtemp.grd -R0/1/0/1

grdfilter temp.grd -Gtemp2.grd -R0/1/0/1 -I0.01 -D0 -Fg0.5

grd2cpt temp2.grd -E100 -D -Cjet > sn.cpt

grdimage temp2.grd -R0/1/0/1 -B0.1 -JX5i -K -Csn.cpt > filter.ps

grdimage temp.grd -R0/1/0/1 -B0.1 -JX5i -O -Csn.cpt -Y6i >> filter.ps

evince filter.ps