<div dir="ltr"><div><div><div>Hi devs, <br></div><div><br></div><div>I found that r.mask does not behave as I would expect and I would really appreciate some clarification. <br></div><div><br></div><div>The issue is the following: I apply a mask, then I do some processing, I remove the mask and when I display the resulting map, all the pixels were processed. I would expect that those pixels that I masked were treated as NULL and not considered exactly as the message says: <br></div><div><br></div><div>"All subsequent raster operations will be limited to the MASK area. Removing<br>or renaming raster map named 'MASK' will restore raster operations to<br>normal."<br></div><div></div><div><br></div><div>Otherwise, what's the point of applying a mask?</div><div><br></div><div>This is an example with NC:<br></div><div><br></div><div>g.region raster=landcover_1m -p<br>r.category landcover_1m<br>r.mask raster=landcover_1m maskcats="1 thru 6"<br>d.mon wx0<br>d.rast landcover_1m<br>r.mapcalc "applied_mask = 1"<br>r.mask -r<br>d.erase<br>d.rast applied_mask <<-- everything is there, no holes... <br></div><div><br></div><div>The only operation that gave me what I would expect when I apply a mask was: <br></div><div></div><div>r.mapcalc expression="applied_mask=1-<wbr>MASK"</div><div><br></div><div>Can someone, please, test/clarify this?</div><br></div>thanks much in advance!<br></div>Vero<br></div>