Monday, September 20, 2010

Isolating hot spots by reclassification and zone area calculations in SAGA GIS

Another way to summarize and visualize a continuous grid surface e.g. a smooth surface that modeled the density of occurrences of point data incidents, is to reclassify the continuous data into discrete groups. An example of a continuous data surface is shown below.

In SAGA GIS, the Reclassify Grid Values module can do the zoning. Once the continuous surface has been classified into groups, meaningful numerical information such as the area of the hot spots can be derived.

The following describes how I use SAGA GIS to reclassify and calculate the zone area values.

Reclassifying into discrete zones
  1. Start up SAGA GIS. Load in a continuous grid layer e.g. hotspot.sgrd.
  2. In the Workspace pane on the left, select the grid layer. In the Options pane on the right, click the Description pane.
  3. Examine the Value Range and note down the values, e.g. 0 - 10. Note: this will be used to determine how to group the data.
  4. Select Modules | Grid | Tools | Values | Reclassify Grid Values.

    The Reclassify Grid Values dialog box appears.
  5. In the Grid system drop down list, choose the grid system e.g. "30;354x420y...".
  6. In the Grid drop down list, choose the grid layer for reclassification, e.g. "01.hotspot".
  7. In the Method drop down list, choose simple table.

  8. Click the Method simple table's Lookup Table row. Click [...].

    The Lookup Table dialog box appears.
  9. Use the minimum, maximum, and new columns to define the input ranges to reclassify, e.g. 0~3 to 1, 3~6 to 2 and so on.

  10. Click Okay twice.

    The grid layer is reclassified.
Visualizing the discrete zones
  1. In the Workspace pane, double click the discrete zone grid layer e.g. "03.Reclassified Grid" and display in a new map window.

    Note: the default colors may not be optimal in showing off the zones as shown below.
  2. On the Options pane on the right, click the Settings tab. Change the Display:Color Classification Type to Lookup Table.

  3. Click the Lookup Table's Table row. Click [...].

    The Table dialog box appears.
  4. Enter the NAME, DESCRIPTION, MINIMUM, MAXIMUM values appropriate for the discrete zones grid layer and define unique colors for each range, as shown below.

  5. Click Okay.
  6. Click Apply.

    The map window is updated with the new colors which show off the zones better.
Calculating the zone areas
  1. Select Modules | Geostatistics | Grids | Zonal Grid Statistics.

    The Zonal Grid Statistics dialog box appears.
  2. Click the Grid system drop down list and choose the discrete grid layer's system.
  3. In the Zone Grid drop down list, choose the discrete zone grid layer, e.g. "03.Reclassified Grid".

  4. Click Okay.

    The zonal statistics table is created.
  5. In the Workspace pane, double click the newly created zonal statistics table e.g. "01.Zonal Statistics" to display it.

    Note: the Count column counts the number of grid cells within a zone.
  6. Select Modules | Table | Calculus | Table Calculator.

    The Table calculator dialog box appears.
  7. In the Table drop down list, choose the newly create table, e.g. "01.Zonal Statistics".
  8. In the Formula field, type in the formula to calculate the zone area, e.g. if the area of a grid cell is 2 square meter, then the zone area formula is b*2, where b refers to the second column.
  9. In the Field Name, type in the new field name, e.g. area.

  10. Click Okay.

    The new table is generated.
  11. In the Workspace pane, double click on the newly created table.

    The new table is displayed. Note the new area field showing the calculated zone area.

No comments: