Bands From Raster (Production Mapping)

Summary

Creates an elevation bands feature class from a Digital Terrain Elevation Data (DTED) raster dataset. The resulting feature class is used by the Elevation Guide Bar surround element to display the high and low elevation areas on a Topographic Line Map (TLM) using a series of hypsometric gray bands.

There can be two, three, or four bands per map sheet. The number of bands that display will depend on the difference between the high and low elevation and the contour interval in the map sheet, with the following designations:

Contour intervals can be 10, 20, or 40 meters by default. If no area of interest (AOI) features are selected, the bands will be calculated for the extent of the raster dataset.

Usage

Syntax

BandsFromRaster_production (in_raster, contour_interval, out_feature_class, {in_aoi_features}, {aoi_field}, {in_hydro_exclusion_features})
ParameterExplanationData Type
in_raster

The input raster dataset.

Raster Layer
contour_interval

Contour intervals with default value options of 10, 20, and 40.

Long
out_feature_class

The output elevation guide band feature class.

Feature Class
in_aoi_features
(Optional)

The optional area of interest feature class to be used in lieu of the current data frame extent to define the map extents.

Feature Layer
aoi_field
(Optional)

Field containing a unique identifier for each map sheet.

Field
in_hydro_exclusion_features
(Optional)

Used to define large water body areas that need to be excluded from the elevation band area calculations.

Feature Layer

Code Sample

BandsFromRaster example (Stand–alone Python script)

This stand-alone Python script uses the BandsFromRaster tool with a grid feature layer AOI.

# BandsFromRaster.py
# Description: Creates polygon bands from a Raster DTED
# Author: ESRI

# Import arcpy module
import arcpy

# check out a foundation license
arcpy.CheckOutExtension("Foundation")
arcpy.CheckOutExtension("Spatial")

# import the toolbox
arcpy.ImportToolbox("C:\\Program Files\\ArcGIS\\Desktop10.0\\ArcToolbox\\Toolboxes\\Production Mapping Tools.tbx")

# Local variables:
polbndA = "C:\\data\\Foundation_Data_Tutorials\\MPS-Atlas\\GDB\\SoCal.gdb\\SoCal\\PolbndA"
dted = "C:\\data\\Foundation_Data_Tutorials\\Reference\\DEM\\dted\\w117\\n32.dt0"
rasBnds = "C:\\data\\outExtents.gdb\\RasterBands"

# overwrite existing data
arcpy.env.overwriteOutput = True

# Process: Make Feature Layer
if arcpy.Exists("polBoundLyr") == False:
    arcpy.MakeFeatureLayer_management(polbndA, "polBoundLyr", "\"OBJECTID\" = 143")

# Process: Bands From Raster
arcpy.BandsFromRaster_production(dted, "10", rasBnds, "polBoundLyr", "gfid", "")

# check the extension in
arcpy.CheckInExtension("Foundation")
arcpy.CheckInExtension("Spatial")

Environments

Related Topics

Licensing Information

ArcView: No
ArcEditor: Requires Production Mapping and Spatial Analyst
ArcInfo: Requires Production Mapping and Spatial Analyst

5/4/2011