[GRASS-user] [Solved] Re: cutting issue (again): r.mapcalc vs. r.resample

[Tim Michelsen]

> I while ago I asked how to cut a small map from a larger one using the
> region settings as boundary for the small map:
> http://n2.nabble.com/r-mapcalc-changes-color-table-td4072836.html#a4072836
> 
> I tried this again:
> 
> Approach 1: r.mapcalc
> 	
> 	r.mapcalc "cut=if(boundary_raster, large_map,
> null())"
> 
> Approach 2: r.resample
> 
> 	r.resample input=large_map output=cut
> 
> I notice the that using Approach 2, the areas that was outside the set
> region is white. The underlayung rasters are not visible (or only if I
> set the opacity of "cut" to a low value).
> If I use Approach 1 with r.mapcalc, the areas outside the region are
> really NULL() (=transparent) and underlaying larger rasters are visible.
> Only that the result "cut" from r.mapcalc has a new colortable.
> 
> Am I doing something wrong?
I did nothing wrong.

But when performing the resample in the wxGUI I activated the button
"add to layer manager" below in the dialog window.
The resampled map ("cut") is then automatically added to the map BUT
without the overlay option of d.rast!
I changed the settings in the layer manager and ticked the -o option and
everything worked fine!

Such a nitty detail!
Filed a ticket for this:
http://trac.osgeo.org/grass/ticket/892


Thanks for all suggestions!