masking

[Rick Thompson]

Hi Robin,
  You can use any raster command-line command to produce a MASK. For example:
r.reclass in=township1 out=MASK
r.buffer in=area1 out=MASK dist=600
r.mapcalc MASK = "if(township$count)"

You can pass variables in the script to help the automation. Hope this helps.
Rick

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

> 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.
------------------------------------------------------------------------
Rick Thompson-   Research Specialist          E-mail: rick at cast.uark.edu
Center for Advanced Spatial Technologies      Telephone: (501) 575-8498
Ozark Hall Rm. 12                             Fax: (501) 575-5218
University of Arkansas              
Fayetteville, AR 72701
DISCLAIMER: Any opinions or comments about any product, organization,
            or individual expressed in the above text are my own.