Introducing the Spatial Analyst Supplemental tools

With the introduction of the capability for sharing custom tools in ArcGIS 10.1, we here in the Spatial Analyst team have endeavored to provide some new functionality as a collection of script tools. These tools can be used as-is for the applications they are identified for, or you can modify them to handle specific use cases that you might have.

The following sections describe what each of these sample tools do and how they fit in with the rest of the available tools.

Create Dendrogram tool

In  the process of classifying multiband raster input using the Class Probability or Maximum Likelihood tools, a signature file is a required input.  In order to have a good classification, ideally there should be little overlap between the individual clusters in the multidimensional attribute space of the input bands.

By analyzing the statistics contained in a signature file, the Dendrogram tool in Spatial Analyst depicts the relationships between the classes graphically in a specially-formatted ASCII text file.  The Create Dendrogram tool gives you the option to output the dendrogram to a PDF file, which gives a better graphic result. By analyzing at the dendrogram, you can determine if the clusters identified in the signature file are distinct enough to give a good classification result.

Draw Signatures tool

When classifying raster images, better results are obtained when there is minimal overlap between the classes. Similar to the Dendrogram tool described above, the Draw Signatures tool provides a way to visualize relationships between the classes, but instead of representing them as a line diagram, this approach displays the class signatures as a series of ellipses.  The ellipses outline the boundary of the sample distribution defined by a specified number of standard deviations.

(Long time users may remember the nature of this output as being similar to that of the DrawSig command in ArcINFO Workstation)

Erase Raster Values tool

Have you ever wanted to just replace certain cells in a raster without having to go through t a sequence of Con, SetNull and other operations?  With the Erase Raster Values tool, you can simply provide an input (feature or raster) that defines the areas of cells you want to change, and then specify if you want those cells to become NoData or some other value in the output raster.

Filled Contours tool

Contours are lines that connect locations of equal value on a raster surface dataset, also known as isolines. They are useful for making areas of rapid change in the phenomena being mapped more apparent and easier to see.

The Filled Contours tool creates a polygon feature class of filled contours from the raster surface.  As compared to regular contour lines, filled contour polygons can be used to easily find the area of elevation values between two particular elevation contours, or quickly summarize all the vegetation types within each contour range.

Maximum Upstream Elevation tool

One type of analysis that agencies like the USGS may perform is to model the potential debris flow an extreme rain event might cause over a mountainous terrain after a wildfire has burned through a portion of an old growth forest.  The Maximum Upstream Elevation tool creates a result that is useful to this process. It uses an elevation raster and a flow direction raster to generate an output raster of the maximum upstream elevation raster along the longest upstream flow path for each cell.

Peak tool

A well-known tool used in hydrological analysis is the Sink tool, which identifies the locations of sinks (localized depressions where all neighboring cells are higher than a particular cell, or group of cells), used in the process of creating a depressionless DEM for further analysis.  The mechanism used for that tool can be applied in such a way to provide essentially the inverse of the result, which would be the locations of all of the localized high points, or peaks, on the surface raster.  The Peak tool takes an input surface raster and outputs a raster with unique cell values for each of the peaks, as well as an optional point feature output.

Tabulate Area 2 tool

The original Tabulate Area tool calculates cross-tabulated areas between two datasets.  However, if the input zones have any overlap, those features will be flattened down into one layer before performing the operation.  This new Tabulate Area 2 tool offers an improvement over the original tool by handling overlapping polygons in feature zone input.

Viewshed Along Path tool

Ordinarily, visibility is determined based on a full 360 degree field of view for each of the specified observer locations.  However, sometimes it would be helpful to constrain the field of view to smaller angle for an observer that moves along a path through a landscape.  For example, imagine being the driver of a parade float or an armored vehicle with a very limited viewport to look out of.  The new Viewshed Along Path tool creates an output that can give an indication of areas of concern along a planned route where the visibility may be of some concern.

Zonal Statistics as Table 2 tool

The original Zonal Statistics as Table tool summarizes into a table the values of a raster within the zones of another dataset.  However, if the input zones have any overlap, those will be flattened down into one layer before performing the operation.  This new Zonal Statistics as Table 2 tool processes the input in a different way, preserving the overlapping nature of zonal inputs.

Downloading and installing

The Spatial Analyst Supplemental Toolbox is available for download here:
http://esriurl.com/SpatialAnalystSupplementalTools

Once you’ve downloaded the file to your machine, extract the contents of the zip file.  Then, in Catalog, navigate to the location you extracted the files to. You should now see the Python toolbox containing the Script tools .

(Note that the samples are geared for ArcGIS version 10.1 and beyond.)

Update for 10.1 SP1

In order to support the use of keywords for parameter datatypes that was introduced in ArcGIS 10.1 SP1, the toolbox has been updated to version 1.1.  One of the benefits of this change is better support for different locales.

The tools can still be used with ArcGIS 10.1 without the SP1 service pack, but some small alterations need to be made to the python scripts.  Simply follow the instructions contained in the ReleaseNotes.txt file in the download.

Version history

  • v1.3 – Added Tabulate Area 2 tool, added Zonal Statistics as Table 2 tool (to replace Zonal Statistics as Table for Overlapping Polygons), and added a new parameter to the Maximum Upstream Elevation tool.
  • v1.2 – Added Erase Raster Values and Zonal Statistics as Table for Overlapping Polygons tools.
  •  v1.1 – Updates for ArcGIS 10.1 SP1.
  • v1.0 – Create Dendrogram, Draw Signatures, Filled Contours, Maximum Upstream Elevation, Viewshed Along Path, Peak, and Viewshed tools made available.
This entry was posted in Analysis & Geoprocessing and tagged , . Bookmark the permalink.

Leave a Reply

13 Comments

  1. jeffrey.vanetten says:

    Juan
    Great addition, however it appears that some of the tools require not only Spatial Analyst but a higher version of ArcGIS to successfully run. We are running ArcGIS Basic and was looking forward to using the Filled Contours tool only to get a license error when the tool got to the Feature to Line tool. Any chance of a work around?

    Thanks
    Jeff

    • Steve Lynch says:

      Jeff

      The Filled Contours tool makes use of FeatureToPoint, FeatureToLine and FeatureToPolygon which required an Advanced license. It also makes use of Contour and IsNull which require a Spatial Analyst license. We have updated the license requirements on the tool box download page.
      As far as a work around, you could look at developing your own tool. A very basic approach would be to use Reclass followed by RasterToPolygon.

      Regards
      Steve

  2. niwayu says:

    Hi
    It may not work on some localized versions (e.g. Japanese) of ArcGIS for Desktop. You need to replace the values of datatype (e.g. “File”) in each “arcpy.Parameter()” functions to the keywords (e.g. “DEFile”) which were introduced in 10.1 SP1.
    Thanks

  3. Steve Lynch says:

    Jeff

    The other tools will work with a Basic license. Filled Contours is the only one that requires Advanced.

    Regards
    Steve

  4. jyh says:

    I tried the Peak script tool but get ” ERROR 010123: Could not get field name” from the Zonal Statistics step. I can see from the script file the steps that are used, and I can run through them manually, but this tool would save me time. Any suggestions?

  5. debbow73 says:

    I want to use the draw signatures tool – but I only have ArcGIS 10.0. Will these tools work for me or do I need 10.1??

    Thanks

  6. Xuguang Wang says:

    The tool would require ArcGIS 10.1 SP1 to work. It does not support 10.0.
    Xuguang

  7. hymiller says:

    hi! the Esera raster value tool is unavailable,the others can be used,i dont kown why? wish your reply!

  8. fstony says:

    Zonal statistics as table 2 is still giving me no output for overlapping polygons.