I have a land use/land cover classification of a SPOT XS image
for a region in northwestern Venezuela. This area includes ten
villages. I would like to define a catchment around each village
(an irregular polygon would be nice but a circle will do) and assign
a name to each of these catchments. Then, I would like a list of the
total area in each LU/LC category for each of the named catchments.
Can this be done in a single step, e.g., by using a script? Has
anyone already written such a script? Perhaps some use of GRASS mask
layers might be used for this?
Thanks in advance of your help.
Concerning your request for script materials for 'site catchment analysis',
I have a script which 1) puts up a raster map
2) puts up sites and labels
3) allows the user to draw a polygon with mouse
4) takes the polygon and makes a vector file out of it
5) converts the vector file into a raster file
6) uses the raster file as a psuedo mask for finding sites
within the polygon ( mask )
This is probably different from what you want to do, but maybe some of the
techniques, hacked together though they be, would be of use. You could view
my stuff as a prototype in sh ( with calls to little FORTRAN executables ) of
what might be better accomplished wioth some judicious C programming!
Anyway, drop me a line if it sounds of interest.
Joe Vaughan
Joseph K. Vaughan
University of Arizona
Institute of Atmospheric
Physics
PAS Bldg, Rm 470
Tucson, AZ 85721
I can also be reached by return mail on the
internet or as
jkv@air.atmo.arizona.edu and by phone at
602/621-9663.