masking

[Steve Crabtree]

Try making a raster map for each area if you have not already done so. 
You can use r.resample to do this. Then use "g.copy rast=name,MASK"

This should solve your problem.

Steve Crabtree
GIS Coordinator
NRCS - Kentucky

On Mon, 8 Jan 1996, Robin L McNeely wrote:

> Thanks in advance for your help.
> 
> I am working with soils data in raster format at Iowa State University.  I
> have masked a township one section at a time and need to produce r.report
> files for each section.  I need to do this for seven maps per section and
> would like to automate the process since eventually I will need to make
> files for an entire county.  My problem is that r.mask is not interactive
> and can't be included in a shell script.  I tried to achieve a mask-like
> effect with g.region but it gave me a 90 degree box, not the skewed boundary
> of a section.  This affects the results listed in r.report such as acres and
> cell count.  
> 
> I am looking for a non-interactive mask command that I hope was written by
> someone intending to automate the masking of many files.  I would also be
> happy to hear if there is a way to achieve similar results with existing
> GRASS 4.1 commands.   Thank you.
> 
> 
> Robin McNeely
> Department of Landscape Architecture
> 
> 
>