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.
Syntax
RasterDomain_3d (in_raster, out_feature_class, out_geometry_type)
Parameter | Explanation | Data Type |
in_raster | The raster to process. | Raster 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.
| 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
Licensing Information
- ArcGIS for Desktop Basic: Requires 3D Analyst
- ArcGIS for Desktop Standard: Requires 3D Analyst
- ArcGIS for Desktop Advanced: Requires 3D Analyst