Raster Domain (3D Analyst)

Summary

Creates a polygon or polyline footprint of the data portions of a raster dataset.

Illustration

Usage

• Single-band or multiband rasters can be provided as input.

• The domain feature will extend to the center of the perimeter cells in the contiguous data blocks of the raster. The cell center defines interpolation zone of a raster surface. NoData cells are ignored in the output feature.

• 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

RasterDomain_3d (in_raster, out_feature_class, out_geometry_type)

Parameter	Explanation	Data Type
in_raster	The input raster.	Raster Layer
out_feature_class	The output feature class.	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, Output Spatial Grid 1, Output Spatial Grid 2, Output Spatial Grid 3

Licensing Information

ArcGIS for Desktop Basic: Requires 3D Analyst

ArcGIS for Desktop Standard: Requires 3D Analyst

ArcGIS for Desktop Advanced: Requires 3D Analyst