On Fri, Sep 11, 2015 at 5:22 PM, Helmut Kudrnovsky <[email protected]> wrote: > > > I define the computational area with 'g.region raster=dem zoom=dem' ('dem' > being my raster) before launching the script so it should match the > elevation map. > > try > > g.region -a raster=elevation align=elevation > > to match and align your dem > > https://grass.osgeo.org/grass70/manuals/g.region.html > > -a Align region to resolution > align=name Adjust region cells to cleanly align with this raster map
consider to also add -p in order to actually see what you have there. Ideally the resolution should match that of the map, also in order to avoid non-square pixels as you report above " nsres: 30.00043309 ewres: 29.99901663 " which are not ideal. best Markus _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
