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
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@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.
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