GEOG 751: Digital Techniques of Remote Sensing

Exercise #10

Image Classification


Objectives



Image

Part I Training Site Selection


A. Signature Extraction

To begin, open a color infra-red composite of cola_tm7_2000-03-06.img in a viewer (RGB= bands 4,3,2) and fit to frame. The ERDAS Imagine Signature Editor allows you to create, manage, evaluate, edit, and classify signatures (.sig extension). Both parametric (statistical) and non-parametric (feature space) signatures can be defined. In this exercise, we will be defining signatures by collecting them from the image to be classified using the Signature Editor and Area of Interest (AOI) tools. The Signature Editor can be accessed through the Classifier icon in the Imagine icon panel. This device will enable you to select and save training sites and make them available for future use in a supervised classification. You may launch the Signature Editor without having obtained any previous signatures or you can retrieve a .sig file using Load under the File menu within the Signature Editor. The Signature Editor has many interesting and useful tools. The tools you should concern yourself with are the buttons directly beneath the menu bar, especially the three that have pluses and minuses on them. These will be used in conjunction with the AOI editor to enter training sites into a .sig file. The first button, it looks like an L with a plus next to it, is used to add a currently selected AOI site to the file. The next one to the right will replace the highlighted field with the current AOI site. The third button is used to merge training sites (signatures) once you feel they have similar spectral characteristics.

Create New Signature(s) from AOI
Replace Current Signature(s) with AOI
Merge Selected Signatures

To gather the spectral signature of the sites you would like to place in the signature editor as training sites, you will need to use the AOI (Area Of Interest) tools. The AOI menu can be accessed through the current viewer's menu bar. In the AOI pull down menu you will be presented with many choices (AOI Styles changes the way the cursor styles look). The Tools and Commands options are important because they allow you to select the type of polygon, modify the polygon, etc. with which you want to encompass your AOI. The Seed Properties option is also important because it allows you to modify the limits of seed area growth by area and/or distance in addition to letting you select the Neighborhood selection criteria. We will be using the Neighborhood default setting which specifies that four pixels are to be searched, then only those pixels above, below, to the left, and to the right of the seed or any accepted pixels are considered contiguous. Under Geographic Constraints, the Area check box should be turned on to constrain the region area in pixels. Enter 500 into the Area number field and press Return. This will be the maximum number of pixels that will be in the AOI. Enter 10.00 in the Spectral Euclidean Distance number field. The pixels that are accepted in the AOI will be within this spectral distance from the mean of the seed pixel. Before closing the Seed Properties window, click on Options and make sure that the Include Island Polygons box is turned on in order to include nonadjacent polygons within the logical growth region.

To begin the process, you must select an area on the image using one of the AOI tools, such as the polygon or rectangle tool, or you can place a seed and grow a region using the Region Grow tool (looks like a magnifying glass in the AOI menu). Use whatever you need in that particular instance, just make sure you think you know what the area represents in terms of ground cover. In the viewer, zoom into an area where you want to select an AOI using the viewer's magnifier tool and then select the AOI polygon tool and draw a polygon around your chosen area (or you may plant a seed to grow). After the AOI is created, a bounding box surrounds the polygon or region, indicating that it is currently selected. While the area is selected, use the Create New Signature button to add the selected area into the Signature Editor. Now click inside the Signature Name column for the signature you just added and give it a name (use names like urban1, urban2, etc. to define your individual AOIs). You may also want to change the color in the Color column. You can use the Image Alarm tool under View in the Signature Editor to get a preview of the extent that the classes you have chosen represent the rest of the image. If you select the Image Alarm option a pop-up box titled Signature Alarm will open. In this box you can choose to indicate classes that overlap and the color that represents overlap. This can be useful if you are considering merging classes. The signature alarm will also, as mentioned, let you see the extent of each of the classes (you will need to do this anyway before you can see the overlap). Do this by selecting (highlighting) a class or a set of classes in the signature editor using the cursor. You can select the color you would like to represent the class as by clicking on the colored square with the right mouse button. Once you have made your selections click on the OK button in the signature alarm and let Imagine do its work. Using this tool, you can see what areas are covered and which are not using the classes you have selected (according to the rules of parallellepiped classification logic).  For this lab, take at least three relatively distinct training sites for each of the following classes found in the Columbia scene:

  1. Urban
  2. Residential (mixture of vegetation, paved surfaces, etc.)
  3. Wetlands
  4. Forest
  5. Water

When you are done generating the training sites for these 5 classes and you feel they are representative of the whole scene based on your use of the signature alarm, save the signature editor file using the Save As menu item under file in the signature editor menu.

B. Feature Selection

The Signature Mean Plot button to the left of the histogram in the signature editor, allows you to view the mean plots of your training data on the screen and thus estimate which of the TM bands best discriminates between the different training sites that you have selected. Select this option as well as the histogram (if you feel this is more helpful, use the all selected signatures and all bands option and make sure you have all of the classes highlighted when you do this). The most precise way to accomplish the task of determining which bands to use is to through the use of the separability option under the Evaluate menu in the Signature Editor. Select 3 for the Layers Per Combination choice. Consult with Dr. Jensen's book to determine which Distance Measure to use and how to interpret the results that will appear in the cell array (choose Cell Array option). Use the Best Average listing method and click OK. The results of this operation will appear in a pop-up box titled Separability CellArray. Note which 3 bands seem to do the best job of spectrally separating your classes. Also note which classes overlap and which are spectrally separable. You will need this information for the next part of the exercise.

Display Signature Mean Plot Window
Display Signature Histogram Window


Part II. Supervised Classification


Now that you have specified your training sites, you are ready to proceed with the supervised classification. Under the Classify menu in the Signature Editor, choose the Supervised Classification option. Because you have already selected a signature file it will not ask for one. If you were to close the signature editor and access the supervised classification through the Imagine Classifier menu, you would be able to open a .sig file. In the Supervised Classification pop-up box that appears give a name for your output file. The Parametric Rule setting should be set to Minimum Distance (note: the textbook gives a good description of the differences, advantages and disadvantages of the various classification logic schemes) and everything else should be left as you find them. Select OK when everything is in place. Open a new viewer and display the results, then answer the following questions concerning signature extraction and supervised classification:

 

Part III. Unsupervised Classification (Clustering)


Unsupervised classification differs from a supervised classification in that the computer develops the signatures that will be used to classify the scene rather than the user. The classification process results in a number of spectral classes which the analyst must then assign (a posteriori) to information classes of interest. This requires a knowledge of the terrain present in the scene as well as its spectral characteristics.

The Unsupervised Classification option is selected in the Classification menu under the Imagine Classifier icon. You will notice that the Unsupervised Classification dialog box states that it is an ISODATA unsupervised classification. The Iterative Self-Organizing Data Analysis Technique (ISODATA) is a widely used clustering algorithm and is different from the formerly used chain method because it makes a large number of passes through the remote sensing dataset, not just two passes. It uses the minimum spectral distance formula to form clusters. It begins with either arbitrary cluster means or means of an existing signature set, and each time the clustering repeats, the means of these clusters are shifted. The new cluster means are used for the next iteration.

The ISODATA utility repeats the clustering of the image until either a maximum number of iterations has been performed, or a maximum percentage of unchanged pixels has been reached between two iterations. Performing an unsupervised classification is simpler than a supervised classification, because the signatures are automatically generated by the ISODATA algorithm. However, as stated before, the analyst must have ground reference information and knowledge of the terrain, or useful ancillary data if this approach is to be successful.

To begin the unsupervised classification, click on the Classification icon and then select Unsupervised Classification... Fill in the input and output information in the Unsupervised Classification dialog box. Give both the Output Cluster Layer and Output Signature Set a similar name. Make sure that under Clustering Options the Initialize from Statistics box is on and set Number of Classes to 30. Set Maximum Iterations to 20 and leave the Convergence Threshold set to 0.950. Maximum Iterations is the number of times that the ISODATA utility will recluster the data. It prevents the utility from running too long, or from getting stuck in a cycle without reaching the convergence threshold. The convergence threshold is the maximum percentage of pixels whose cluster assignments can go unchanged between iterations. This prevents the ISODATA utility from running indefinitely. Leave everything else in its default state. When you have entered all the relevant information click OK to begin the process.

Part IV. Cluster Identification

To aid in evaluation we will need to view the results of the clustering so that we may see how the clusters are arranged in feature space and thereby make informed decisions about the nature of the cluster. The first step that will allow us to do so is the creation of feature space images. The Feature Space Image button can be found on the Classification menu. When it has been selected a dialog box will appear saying Create Feature Space Images at the top. Select the original image (not the clustered one) as the Input Raster Layer and make sure the Output Root Name is similar to the raster layer and the directory path is correct. Leave the rest of the selections at their default settings and click OK. When the processing is complete open a new viewer and view the output images (i.e. the cola_tm7_2000-03-06.fsp.img file as the raster layer). Note that the 2_5 (and other options) represents the layers that are being shown in the image. In this case layer 2 (band 2) will be displayed on the x-axis and layer 5 (band 5) on the y-axis. Pay close attention when you look in the book for help in determining what clusters represent which ground elements.

The next step is to open the Signature Editor (under the Classification menu) with the *.sig file you created in the unsupervised classification. Select all the clusters (they should all be highlighted in yellow). In the Signature Editor main menu select Feature and then in that pull-down menu select Objects This will display a Signature Objects dialog box that allows you to tell Imagine which viewer you want to receive the signature editor information about the clusters. In this case we want the viewer in which you have displayed your chosen feature space graphic. Select that viewer # in the Signature Objects space provided that represents this viewer. Select Plot Ellipses and Plot Means (or you can try the others if you like). Leave everything else in it's default state and click OK. Only selected clusters in the Signature Editor window will be drawn. More than likely your ellipses and means are multi-tonal in nature. If you would like them all to be white, red, green, etc... select all the classes in the Signature Editor dialog box using the mouse and change the color to the one you desire. Save the information as an Annotation Layer.

To analyze the content of the clusters, you should use a combination of techniques. The first should be to use the mean scatter plot to make some educated guess about the information in each cluster. You might want to label each of the 30 clusters on the scatter plot using the Label option in the Signature Objects dialog box so you know what cluster is containing which class. You will more than likely have to zoom in to get a better look at some of the clusters given the close proximity of clusters to each other. You should also have a viewer open with the original scene displayed. This will further help you identify the land cover class. If you are feeling adventurous, you can overlay your classified image on the original image and set all the clustered image's colors to transparent using the Raster Attribute Editor found under the Viewer menu (Raster - Attributes). Once you have set all classes to transparent then you can individually color (by making them opaque) particular classes and see where they are on the image. Another method may be to use the Utility - Swipe or the Utility - Fade tools in the Viewer by opening the classified image on top of the raw data (do not Clear Display after opening the first raw image). Regardless of how you choose to proceed you should not rely on any one particular method but a combination of methods and some common sense to arrive at a sound classification.

When you have decided upon the class breakdowns, use the Raster Attribute editor to assign class names and colors to the classification image. Create the same five classes you used in the supervised classification and place each of the 30 clusters into one of the classes by giving it the same color and class name as every other cluster in that class. When you satisfied with your unsupervised classification, finish the lab by doing the following:






Back to Exercises



Last Modified: September 27, 2004