<div dir="ltr"><div><div>I've narrowed the problem down to creating a masked copy using r.mapcalc. The LOS output works fine, but the LOS with elevation values does not. Is there another way to create this? Or can I apply a mask based on the area of one raster and the value of another?<br><br></div>Ta,<br></div>Mike<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 25, 2015 at 3:56 PM, Michael Spencer <span dir="ltr"><<a href="mailto:spencer.mike.r@gmail.com" target="_blank">spencer.mike.r@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>No problem.<br><br></div>The working mask is based on a DEM (<a href="http://www.ordnancesurvey.co.uk/business-and-government/products/terrain-50.html" target="_blank">http://www.ordnancesurvey.co.uk/business-and-government/products/terrain-50.html</a>), created with:<br><br></div>r.reclass input=DEM@Permanent output=DEM_reclass@Permanent rules=elev.txt<br><br></div><div>elev.txt<br>0 thru 150 = 0<br>150 thru 300 = 150<br>300 thru 450 = 300<br>450 thru 600 = 450<br>600 thru 750 = 600<br>750 thru 900 = 750<br>900 thru 1050 = 900<br>1050 thru 1200 = 1050<br>1200 thru 1500 = 1200<br><br></div><div>The none working mask uses a raster based on the above, created using:<br><br></div><div># in A mapset<br></div><div>r.mask -o input=viewshed maskcasts='*'<br></div><div>r.mapcalc 'subset=DEM_reclass@Permanent'<br><br></div><div>Strictly speaking, this step was done via spgrass6 as I have 145 subsets to apply:<br>lapply(st[,1], function(i){<br> # Which site?<br> print(i)<br> # Mask to LOS<br> system(paste0("r.mask -o input=", i, " maskcats='*'"))<br> # Reclass via mapcalc<br> system(paste0("r.mapcalc '", i, "_elev=Terrain50_SSGB@PERMANENT'"))<br>})<br></div><div><br></div><div>However, The rasters I want to mask exist in B mapset and I understand the mask and rasters to be applied to need to be in the same mapset, so I have:<br></div><div>g.copy rast='subset@A',subset<br></div><div>prior to applying the mask<br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 25, 2015 at 3:20 PM, Anna Petrášová <span dir="ltr"><<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, Feb 25, 2015 at 9:58 AM, Michael Spencer <span dir="ltr"><<a href="mailto:spencer.mike.r@gmail.com" target="_blank">spencer.mike.r@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div>I'm getting odd results from using r.mask. The expected behaviour is that when applied it masks off all rasters in the working mapset, however what I'm getting is it only masks off the raster used as a mask.<br><br></div>This doesn't happen when I create a mask using a different raster, this behaves as expected and masks all rasters in the current mapset.<br><br></div>Detail:<br><br></div>Working mask: reclass of another raster<br></div>Not working mask: subset of above, created by applying a mask and using mapcalc (newmap=oldmap).<br></div></div></div></blockquote><div><br></div></div></div><div>Could you post the actual commands you used?</div><div><br></div><div>Thanks,</div><div><br></div><div>Anna</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><div dir="ltr"><div><div><br></div>Any help greatly appreciated!<br></div>Mike<br></div>
<br></span>_______________________________________________<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="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote></div><br></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>