PyGMT v0.1.1 assumes that input
xarray.DataArray grids are gridline registered (
GMT_GRID_NODE_REG ), and of cartesian type (
GMT_GRID_IS_CARTESIAN ). Currrently we are working on allowing for either pixel/gridline registered and cartesian/geographic grids at #476, by changing the
virtualfile_from_grid mechanism. But we’ll need to decide on a good default going forward (i.e. for v0.2.0).
GMT defaults to using
gridline registration (aka data values represent top left corner of pixel), and PyGMT currently follows this convention. Should we change it to
pixel registration which is what
xarray follows? Place your vote below!
- Gridline (left image)
- Pixel (right image)
- Cartesian, i.e.
- Geographic, i.e.
Feel free to leave your comments below, or on Github at https://github.com/GenericMappingTools/pygmt/issues/487 where there will be some more context.