vegetation typing

Hello, GRASS users.

Could someone tell me if GRASS has the capability to type vegetation areas
using infrared reflectance. I have a key made from a small portion of the area that we are studying that has each vegetation type catagorized. I have been
told that this is possible but I have no idea how to do it. This was not my
idea and I did not start this task but it has been left for me to finish. Is
there software that will do such a thing? Any correspondance would be
appreciated. If I need to elaborate on this please let me know.

Thank you in advance,

Beau Bush
bbush@ticsys.tamu.edu

Beau Bush (bbush@ticsys.tamu.edu) writes on 27 Jul 94:

Could someone tell me if GRASS has the capability to type vegetation areas
using infrared reflectance. I have a key made from a small portion of the area that we are studying that has each vegetation type catagorized. I have been
told that this is possible but I have no idea how to do it. This was not my
idea and I did not start this task but it has been left for me to finish. Is
there software that will do such a thing? Any correspondance would be
appreciated. If I need to elaborate on this please let me know.

[my apologies if this sounds rudamentory. Hopefully this will be
of help to a few.]

Do you have several reflectance images, each representing a different
wavelength? If so, you want to do something called "supervised
multispectral classification."

There are two separate algorithms in GRASS for doing this:
maximum likelihood and sequential maximum a prior estimation.
(i.gensig/i.maxlik and i.gensigset/i.smap)

i.gensig/i.gensigset work on your training data. They statistically
define where in n-dimensional space lies each type of vegetation.
i.gensigset can generate "multimodal" training statistics
(a strength of i.smap)

i.maxlik/i.smap do the classification, which is the next step.

These two steps are fairly easy. It's the preparation that might
make things confusing to some. Particularly, importing the
data and somehow figuring out how to assign/delineate training data.
For your case, perhaps r.in.poly might be the best route since
you already know which sections of the image to do the training on.
Also, 'g.manual imagery' will help you understand the GRASS concept
of "groups," which is important for every step of the process.

[so, which do you use, maxlik or smap? depends on your "scene"
and your objective. there's a draft manuscript at
ftp://pasture.ecn.purdue.edu/pub/mccauley/papers/smap-echo-ml-compare.ps.gz
where Bernie Engel and I look at the performance of i.smap versus
i.maxlik and another algorithm. i.smap did the best for "homogeneous"
fields (usually found in cultivated scenes).]

Beau Bush
bbush@ticsys.tamu.edu

--
Gig 'em,
Darrell McCauley '90 '92

P.S. (for Aggies only) Robert Maggio over in Forest Science may be
a good person to get to know if you're going to be doing a lot
of remote sensing work.

We are getting very small DN values for TM raw data:
System: SunOS 4.1.3
GRASS version 4.1
Problem:
    When the bands 1-5 and 7 are displayed in the form of a
histogram, DN values never exceed the range 1-31. In fact,
the histogram display seems to be lumping these categories
into four greyscale levels. It appears that a great deal of
the information contained in the data set was lost when I
imported the files into GRASS. Furthermore, after
examining the statistics obtained from running an
unsupervised classification in GRASS, equally puzzling
numbers appeared. For example, mean DN values for clas ,
representing some type of forest on the ground, never exceed
10. Previous work with TM data sets tell me that highly
vegetated areas usually have a high mean DN value in Band 4.
This data set tells me that the mean DN for class 1 in band4
is 6.673981, and this number doesn't make sense. In
addition, class 15, representing clouds, has the highest
DNs: low 15.983657 (band1) to high 28.832482 (band2).

    The bands were imported by first converting them from
PGA files to sunraster format in xview. A flag in the xview
window signaled that the images "contained all 32
colors" . When I tried to save these bands to the
greyscale option the operation failed. In any case,
after all the initial parameters were set in GRASS, the
files were imported using the "r.in.sunrast" command.

   Any suggestions would be much appreciated.

Jonathan Deenik
c/o rsyost@Hawaii.Edu

-----
Russell Yost Office: 808-956-7066
University of Hawaii at Manoa FAX: 808-956-6539
1910 East West Road rsyost@uhunix.uhcc.Hawaii.Edu
Honolulu, HI, 96822