<div dir="ltr"><div><div>For cropping, see also a new addon module called r.clip which clips according the computational region while respecting mask and the *original* raster resolution:<br><br><a href="https://grass.osgeo.org/grass72/manuals/addons/r.clip.html">https://grass.osgeo.org/grass72/manuals/addons/r.clip.html</a><br><br></div>Testing and feedback welcome.<br><br></div>Vaclav<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 26, 2017 at 10:35 AM, Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 26/09/17 16:17, Moritz Lennert wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On 26/09/17 15:57, Marco Alicera wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear list,<br>
I am trying to crop a raster map into a smaller area, but I am failing.<br>
<br>
This wiki page link <<a href="https://grasswiki.osgeo.org/wiki/Cropping_maps" rel="noreferrer" target="_blank">https://grasswiki.osgeo.org/w<wbr>iki/Cropping_maps</a>> is<br>
quite clear:<br>
"Make boundary into mask using r.mask<br>
<<a href="https://grass.osgeo.org/grass72/manuals/r.mask.html" rel="noreferrer" target="_blank">https://grass.osgeo.org/grass<wbr>72/manuals/r.mask.html</a>> then copy raster<br>
source to new map (g.copy<br>
<<a href="https://grass.osgeo.org/grass72/manuals/g.copy.html" rel="noreferrer" target="_blank">https://grass.osgeo.org/grass<wbr>72/manuals/g.copy.html</a>>). All cells<br>
outside of the mask will be omitted from the new copy."<br>
<br>
r.mask vector=myarea<br>
g.copy --overwrite raster=OldRasterWithBigArea,Ne<wbr>wRasterWithSmallArea<br>
r.mask -r<br>
<br>
I expected all cell outside were not going to be omitted, but replaced<br>
by NULL.<br>
</blockquote>
<br>
I'm not sure I understand your last sentence, but AFAIK, the wiki info<br>
is wrong. g.copy works directly on the files and ignores mask and region<br>
settings. In your example, OldRasterWithBigArea and<br>
NewRasterWithSmallArea should be strictly identical.<br>
<br>
You need to use r.mapcalc "new = old" while the mask is set to get a<br>
raster with everything outside the mask set to NULL.<br>
</blockquote>
<br></span>
I've rewritten the wiki page to try to make things clearer.<div class="HOEnZb"><div class="h5"><br>
<br>
Moritz<br>
______________________________<wbr>_________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/grass-user</a></div></div></blockquote></div><br></div>