|
according to the algorithm of Shi and Dozier (1996). For
this project, theoutput of this algorithm which is based on
a decision tree classificationis an image that is partitioned
into four classes:water, bareground/short vegetation,forest,
and snow (figure 1).
For validation, a Landsat Thematic Mapper (TM) image
acquired on April 14, 1994 over the same area was used.
We performed two types of snow mapping on this image.
For our primary validation image, we used the algorithm
of Rosenthal andDozier (1996) to map sub-pixel snow frac-
tion for the scene (figure 4) [note: to facilitate comparison,
the figures are not in order]. This algorithm has demon-
strated accuracies on par with those obtainable from air-
photo interpretation (Rosenthal and Dozier, 1996). Adding
to our confidencein this method is knowledge that its ac-
curacy assessment was performed over Mammoth Moun-
tain, eliminatingquestionsabout the geographicportability
of the algorithm. As a secondary validation image, we
classified theTM sceneinto the above four categoriesusing
a decision tree classifier (figure 2). This classification
remains unverified so any conclusions drawn from
information derived through it must be tempered. To
facilitate intercomparison, wecoregistered andwarpedthe
TM products to the SIR-C ground range geometry.
The primary validation required a method of compar-
ing the output of the snow fraction image (quantitative or
ratio data) to the output of the SIR-C algorithm (qualitative
or nominal data). These are two fundamentally different
data types. In order to compare we needed to transform
one of the data types to the other. In other words, we could
either synthesize a fractional snow cover image from the
SIR-C classification, or we could generatea nominalimage
of snow-covered/snow-freefrom the TM snowcoveredarea
|
map. The transformation from fractional SCA to a binary
SCA image would require setting some snow fraction
threshold above which pixels are classified as snow and
below which pixels are classified as snow-free. Sincethis
process wouldinvolve a significant loss of information and
the choice of a threshold seemed arbitrary, we pursued the
synthesized fractional SCA SIR-C image option.
One method for synthesizing a fractional SCA image
from a nominal snow/snow-free image is to select a win-
dow size in which to calculate the percentage of pixels
that are classified as snow. This window size is bounded,
at one extreme, by a window size of one in which all out-
put pixels are either 100% or 0% snow and, at the other
extreme, by a window size equal to the entire image in
which the output is simply the fraction of the image that is
classified as snow cover. Thus wehavea situation in which
we must weigh the value placedon cell resolution against
the value of fractional precision. We selected a window
size of 10 by 10 pixels. This choiceseemedreasonablesince
it maintained a moderate resolution with a large number
of samples in the scene and the conversion to percent is
conceptually simple. Since the square of 10 is 100, the
percentis simply the number of pixels in the window that
were classified as snow (figure 3)..
Withtwo coregisteredimages of fractional SCA, all that
remained wassimply to createa differenceimage inwhich
every pixel in the image is the SIR-C fraction minus the
TM fraction. Positive numbers on this image, therefore,
represent areas where SIR-C finds more snow than TM
while negative numbers on this image are areas where TM
finds more snow than SIR-C (figure 5).
In addition to this image, we deemedit useful to perform a
|
 |