Major Description
Related to the usage of `org.apache.sis.coverage.grid.GridRoundingMode` :
At the moment, its usage is often considered as a global configuration for every axes treated (Example: GridDerivation.rounding()).
A possible improvement would be to specify the rounding strategy independently for each axis of the grid. One possible use-case would be to ease management of 2D resample over 3D+ datasets.
For example, if the source coverage of a resample defines a grid [ 0..99, 0..99, 0..0 ], with the axes being respectively longitude, latitude and time, the pipeline of operation over the grid could extend all axes, even time slice depending on chosen rounding mode. Such behavior can fail resampling operation, but would be easily avoidable if we can specify that time is a special axis to round using NEAREST strategy.