How to create a heatmap

[PlanetGus]

If you have a good data density, you can play with the bin size, and you need a more continuous colorbar (makecpt).
Ultimately you can apply grdfilter after gmtbinstats

[Raisha]

just to make sure, the gmtbinstats replace the xyz2grd function in your example, right?

[PlanetGus]

Yep.
Have a look at the thread, there’s a second script

[mkononets]

@PlanetGus

gmt xyz2grd gps.txt -Ggridded.nc -I0.625/0.5+e $region $projection -An -di0

Did you try adding -rp option to produce pixel registered grid (instead of default grid line registration)? My idea is to avoid shifting grid lines when plotting, but IDK whether it actually works the way I imagine it.

Cannot try myself as the data haven’t been posted.

[mkononets]

@mkononets

adding -rp option to produce pixel registered grid (instead of default grid line registration)? My idea is to avoid shifting grid lines when plotting, but IDK whether it actually works the way I imagine it.

Cannot try myself as the data haven’t been posted.

Yes it does, changing grid registration to pixel with -rp makes the grid to match the grid lines without manual shifting, checked with some of my own data.

UPD xyz2grd -A... -rp produces essentially a different grid in a way that many data points fall into neighboring grid cells when grid registration is changed from gridline to pixel.