Raster Domain

Summary

Constructs a 3D polygon or polyline delineating the height along the boundary of a raster surface.

Illustration

Usage

• The tool is designed to operate on a raster surface and will only operate on a single-band raster or a specific band from a multi-band raster. The output feature captures the boundary of contiguous cells with data values, and its z value is interpolated along the perimeter of the surface. The perimeter is defined by connecting the center of the outermost cells in the raster.

• Output geometry is placed in one feature record and may comprise of multipart feature if the raster has discontinuous data cells separated by NoData cells.

  Note:

  3D polygons only contain elevation values along the perimeter of the features, as interior portions of the polygon will not contain any vertices. When drawn in 3D with an area fill, the boundary vertices are arbitrarily connected into triangles for rendering. Unless the polygon is planar, either sloped or horizontal, it's unlikely the fill will accurately represent the interior surface. For this reason, it is recommended that nonplanar 3D polygons are drawn without fill symbology.

Syntax

arcpy.ddd.RasterDomain(in_raster, out_feature_class, out_geometry_type)

Parameter	Explanation	Data Type
in_raster	The raster to process.	Raster Layer; Mosaic Layer
out_feature_class	The feature class that will be produced by this tool.	Feature Class
out_geometry_type	The geometry of the output feature class. LINE —The output will be a z-enabled line feature class. POLYGON —The output will be a z-enabled polygon feature class.	String

Code sample

RasterDomain example 1 (Python window)

The following sample demonstrates the use of this tool in the Python window.

import arcpy
from arcpy import env

arcpy.CheckOutExtension("3D")
env.workspace = "C:/data"
arcpy.RasterDomain_3d("dtm_grd", "raster_domain.shp", "POLYGON")

RasterDomain example 2 (stand-alone script)

The following sample demonstrates the use of this tool in a stand-alone Python script.

'''*********************************************************************
Name: RasterDomain Example
Description: This script demonstrates how to use the 
             Raster Domain tool to generate polygon footprints for all
             *.img rasters in a given workspace.
**********************************************************************'''

# Import system modules
import arcpy
from arcpy import env

# Obtain a license for the ArcGIS 3D Analyst extension
arcpy.CheckOutExtension("3D")

# Set environment settings
env.workspace = "C:/data"

try:
    # Create the list of IMG rasters
    rasterList = arcpy.ListRasters("*", "IMG")
    # Verify there are rasters in the list
    if rasterList:
        # Loop the process for each raster
        for raster in rasterList:
            # Set Local Variables
            outGeom = "POLYGON" # output geometry type
            # The [:-4] strips the .img from the raster name
            outPoly = "domain_" + raster[:-4] + ".shp"
            print "Creating footprint polygon for " + raster + "."
            #Execute RasterDomain
            arcpy.RasterDomain_3d(raster, outPoly, outGeom)
        print "Finished."
    else:
        "There are no IMG files in the " + env.workspace + " directory."
        
except Exception as e:
    # Returns any other error messages
    print e.message

Environments

• Current Workspace
• Scratch Workspace
• Output Coordinate System
• Geographic Transformations
• Extent
• XY Resolution
• XY Tolerance
• Output XY Domain
• Auto Commit
• Output CONFIG Keyword

Licensing information

• Basic: Requires 3D Analyst
• Standard: Requires 3D Analyst
• Advanced: Requires 3D Analyst