Continuation of: SECTION II. MAP DEVELOPMENT
IIA. Upland Land Cover Classification from LANDSAT Thematic Mapper Data
From: Influences on Wetlands and Lakes in the Adirondack Park of New York State:
A Catalog of Existing and New GIS Data Layers for the
400,000 Hectare Oswegatchie/Black River Watershed,
1997

 

Old Forge Quadrangle Test Area OB2 Contents

Because file sizes and, therefore, processing times were quite large, the Old Forge 7 1/2' quadrangle was chosen to develop and test image processing techniques and methods. Old Forge is in the center of the southern scene and upland land cover types expected to be found in the study area are well-represented. An initial perusal of the air photos (1:58000 May 12, 1986 color infrared transparencies) produced field check sites. Detailed notes of upland land cover were made on a 7 1/2' X 15' 1:25000 USGS quadrangle during a field trip designed specifically for this phase of the project. A transparency was also made that fit over the air photos and indicated wetland areas. Upland cover types were sketched on this transparency to aid in signature and classification evaluation.

The Old Forge test area was subset from both the wetland GIS file and from the May and July LANDSAT imagery. In addition, a scanned NYS planimetric map of the Old Forge quadrangle in TIFF format (NYSDEC, 1992) was remade into an ERDAS image file with Image Alchemy and was subsequently brought into ArcView for registration. This latter file was used to help develop signatures in residential areas.

Because so much land cover information was known about the Old Forge quadrangle, a major concentration was directed on developing supervised classification techniques for the area. In total, 44 Old Forge GIS files, 29 image files (including the planimetric file), and 45 individual band calculations were created to evaluate signature and classification options.

Signature training polygons were digitized on the computer screen. Imagery was shown at 3X magnification and three to four band combinations were used with both the May and July imagery to highlight and accentuate various signatures. Signature statistics were created from training polygons with SIGEXT and viewed with SIGMAN (Table IIA.3).

 


Table IIA.3. Signature statistics produced by SIGMAN.

Signature Listing
Signature File: E:UPLALLA/NWSIGS
Signature Statistics Listing

Signature Name: CONIFER

Number of points = 240

  
Band 1 2 3 4 5 6 7 8 9 10 11 12
Minimum 58 23 18 53 34 9 49 17 11 50 27 7
Mean 61.93 25.30 20.59 70.26 45.26 14.15 52.82 19.80 15.00 68.30 39.54 11.43
Standard 1.55 1.18 1.45 7.83 6.52 2.67 1.61 1.07 1.24 10.41 6.77 2.32
Maximum 70 29 26 89 66 28 58 23 19 93 59 22
Covariance Matrix
1 2.41 0.29 0.62 0.13 1.26 1.18 0.62 0.34 0.40 -0.43 0.96 0.64
2 0.29 1.40 0.46 5.24 3.78 1.21 0.54 0.45 -0.07 8.34 4.72 1.10
3 0.62 0.46 2.11 -1.86 0.70 1.38 1.13 0.82 0.88 -1.30 0.67 0.78
4 0.13 5.24 -1.86 61.29 24.75 2.60 -0.28 0.26 -2.64 65.24 28.60 4.42
5 1.26 3.78 0.70 24.75 42.56 12.81 -0.14 0.77 -1.22 39.55 36.09 10.12
6 1.18 1.21 1.38 2.60 12.81 7.11 0.84 0.81 0.45 6.74 11.08 4.08
7 0.62 0.54 1.13 -0.28 -0.14 0.84 2.59 0.99 0.92 1.10 0.34 0.62
8 0.34 0.45 0.82 0.26 0.77 0.81 0.99 1.15 0.59 2.76 1.31 0.59
9 0.40 -0.07 0.88 -2.64 -1.22 0.45 0.92 0.59 1.53 -3.17 -0.98 0.29
10 -0.43 8.34 -1.30 65.24 39.55 6.74 1.10 2.76 -3.17 108.26 48.09 9.22
11 0.96 4.72 0.67 28.60 36.09 11.08 0.34 1.31 -0.98 48.09 45.88 12.07
12 0.64 1.10 0.78 4.42 10.12 4.08 0.62 0.59 0.29 9.22 12.07 5.37

The initial signature polygons were developed based upon field notes on the topographic map and examination of aerial photographs. Both Conifer and Deciduous forest cover types were taken from the center of relatively pure stands. Additional signature polygons were developed based upon comparisons of the scanned planimetric map, aerial photography, and field notes in an effort to discover the separability of Urban Open, Commercial, and Mixed Urban/Vegetation categories. After initial perusals of the signatures, additional Conifer, Cloud, and Cloud Shadow samples were taken.

Within-class signatures were examined with ELLIPSE (Figure IIA.2). Samples that were distinctly different from others of a similar category were removed as were samples that encompassed too broad a distribution of digital numbers. Widely different signature statistics must be considered to represent poor samples in within-class comparisons.

ERDAS SIGNATURE ELLIPSE PLOT

ob2fg2a2.gif (2643 bytes) Figure IIA.2. Within-class signatures shown with ELLIPSE. Note the strong overlap between the two signature samples. These two samples were eventually merged with other samples to form the training signature Deciduous. Axes represent the digital numbers of each band while the scatterplot shows the distribution of digital numbers.
ob2fg2aa.gif (1752 bytes)

Signatures were merged into classes, and between-class separability was evaluated with ELLIPSE (Figure IIA.3). While some overlap of signature statistics is expected, particularly with such a simple representation of distribution as ELLIPSE, excessive overlap suggests inseparable classes. An initial signature file with the following classes was created: Deciduous, Mixed, Conifer, Open, Barren, Urban Open (non-vegetated urban), Urban Mixed (residential and forest), Cloud, and Cloud Shadow. It seemed important to attempt a separation of urban land covers even though ELLIPSE did not support such separability. CMATRIX, which classifies using the signature samples as the data set (Table IIA.4), was used with the May imagery and the July imagery with each of the three primary classifiers: Maximum Likelihood, Minimum Distance, and Mahalanobis. The Maximum Likelihood classifier demonstrated the best overall discrimination between classes.

ERDAS SIGNATURE ELLIPSE PLOT

ob2fg2a3.gif (2763 bytes) Figure IIA.3. Between-class signatures shown with ELLIPSE. Deciduous and Conifer are distinct, illustrating excellent between-class separability. However, Shrub and Deciduous show strong overlap in bands 1 and 4. This overlap was consistent with all band comparisons, indicating that Shrub and Deciduous were not able to be separated into distinct classes.
ob2fg2ab.gif (1752 bytes)

 


Table IIA.4. CMATRIX output.
CMATRIX classifies using the signature sample pixels as the data set.

Image Classification Contingency Table Page 1- 1
Signature File: NWSIGS
Image File : NWXWTLRF
Maximum likelihood classifier

Signature
Name

Training polygon names

WETLAND OPENVEG BARREN CONIFER CLOUD07
WETLAND 186 100.0% 0 0.0% 0 0.0% 0 0.0% 0 0.0%
OPENVEG 0 0.0% 257 90.8% 9 1.5% 3 1.2% 0 0.0%
BARREN 0 0.0% 7 2.5% 584 96.8% 0 0.0% 3 0.2%
CONIFER 0 0.0% 0 0.0% 0 0.0% 247 95.4% 0 0.0%
DECIDUOUS 0 0.0% 3 1.1% 0 0.0% 0 0.0% 0 0.0%
MIX 0 0.0% 10 3.5% 0 0.0% 8 3.1% 0 0.0%
CLOUD07 0 0.0% 6 2.1% 10 1.7% 0 0.0% 1773 99.8%
CLDSHD07 0 0.0% 0 0.0% 0 0.0% 1 0.4% 0 0.0%
Total Points 186   283   603   259   1776  
 
Image Classification Contingency Table Page 2- 1
Signature File: NWSIGS
Image File : NWXWTLRF
Maximum likelihood classifier
Signature
Name

Training polygon names

CLDSHDO7 MIX DECIDUOUS
WETLAND 0 0.0% 1 0.1% 0 0.0%
OPENVEG 0 0.0% 42 3.5% 1 0.4%
BARREN 2 0.2% 1 0.1% 0 0.0%
CONIFER 11 0.8% 27 2.2% 0 0.0%
DECIDUOUS 0 0.0% 213 17.7% 209 87.1%
MIX 0 0.0% 921 76.4% 30 12.5%
CLOUD07 0 0.0% 1 0.1% 0 0.0%
CLDSHD07 1308 99.0% 0 0.0% 0 0.0%
Total Points 1321   1206   240  

As a result of ELLIPSE signature analysis, the Commercial class was dropped because all samples demonstrated huge statistical variations. All non-forested areas, such as residential open, residential vegetated, residential trailer park, commercial/industrial, open fields, airports, ski slopes, golf courses, gravel pits, and shrubby agricultural were combined into open categories and were further subdivided into Open and Open with Vegetation. While Barren (such as bare rock and gravel pits) was a desirable category, the available sample areas were very small and the class tended to strongly overlap urban non-vegetated.

A fourteen-band image file was created for Old Forge with all bands from both the May and July imagery and the Old Forge subset of the wetland GIS file was used to mask out wetlands. An initial classification was performed on the 14-band (all LANDSAT bands from the southern May and July scenes) and two 7-band (May all bands, July all bands) images. Signatures were evaluated both on-screen and with THRESH and additional signatures samples taken. Different a priori values were also tested. Some signatures in particular showed strong banding on the screen display. Consequently, individual bands in each scene were evaluated for striping (Table IIA.5). Although the thermal band (band 6) appeared to be very helpful in classifying urban areas and clouds, it created strong banding in all classifications and was removed from analyses. Numerous standard image processing band manipulations were tried and evaluated (Table IIA.6). Some band combinations demonstrated interesting and bizarre colors that did little to increase knowledge of the sample study area.


Table IIA.5. Notes on Image Banding (Striping)

 

Band North South
May 1992 July 1990 May 1992 July 1992
1 some subtle subtle slight banding banding with 1 standard deviation
2 very slight almost imperceptible banding ok heavy striping with digital numbers 14 - 30
3 some no pattern very streaky almost imperceptible banding heavy striping with 1 standard deviation
4 ok ok ok ok
5 ok ok ok ok
6 real streaky real streaky slight banding unusable
7 ok ok ok; lots of clouds (SE) slight banding
Very heavy banding noted in South July 1992 bands 1, 2, and 3 in the southern section of the Old Forge quadrangle when the image was spectrally stretched through the available digital numbers (as opposed to 0-255). By removing haze on the July image, banding is increased because the dynamic range is decreased. Also band ratioing decreases the digital number range, tending to produce images (e.g., 3/7, 3/4, and 5/2) that are very gray.
Continue reading next section of OB2 Report -- Land Cover continued OB2 Contents