[GRASSLIST:2324] Re: Removing a section of a raster map

Hal Mueller hal at seanet.com
Wed Aug 8 17:50:06 EDT 2001

I haven't tried to do anything like this in a while, and might have 
the details a bit off, don't have time to pin it down right now.

I think the key is r.mask.  You can derive the mask by comparing 
layer with protrusion to layer without protrusion.  Or maybe just by 
looking at layer without protrusion.

(hip shot, untested!)

I would try

to make a binary raster, 1 where I want data and 0 where I don't

set the mask to be that binary I just made

r.mapcalc newrast=rastwithprotrusion

g.remove rast=MASK

d.rast newrast
should display a correct layer

make sure to g.remove the mask (or use r.mask to set no mask) before testing
your new layer; otherwise you might just _think_ you got it right


Hal Mueller                hal at seanet.com
Mobile Geographics LLC     http://www.mobilegeographics.com/
Seattle, Washington        (206) 297-9575
Get CelestNav 2.2.2! http://www.mobilegeographics.com/celestnav/ 


More information about the grass-user mailing list