Summary
Converts the U.S. Federal Aviation Administration (FAA) Digital Obstacle File (DOF) to obstruction points and obstruction buffer features.
Usage
The U.S. FAA DOF specification and dataset can be reviewed at https://www.faa.gov/air_traffic/flight_info/aeronav/digital_products/dof/. The table is a comma-separated values file (.csv).
Output obstacle buffers are created at a distance of 10 times the obstruction object height (meters) based on criteria from US Army Field Manual 3-21.38 Pathfinder Operations, Chapter 4 Helicopter Landing Zones, Section 4-16. For example, a 5 meter obstacle will have a 50 meter buffer.
The Input Table must contain numeric AMSL and AGL fields pursuant to FAA DOF standards.
Syntax
arcpy.intelligence.DOFToObstacleFeatures(in_table, out_obstacle_features, out_obstacle_buffers, {clip_features})
Parameter | Explanation | Data Type |
in_table | The input DOF table to convert into obstacle features. | Table View |
out_obstacle_features | The point obstacle features created from the Input Table. | Feature Class |
out_obstacle_buffers | The distance buffers created at 10 times the value of the AGL field in the Input Table. | Feature Class |
clip_features (Optional) | An area to clip from the Input Table. Only obstacles within this area will be created and buffered. | Feature Layer |
Code sample
The following Python window script demonstrates how to use the DOFToObstacles function in immediate mode.
import arcpy
arcpy.DOFToFeatures_intelligence("C:/data/dof.csv",
"C:/data/results.gdb/obstacles",
"C:/data/results.gdb/buffers",
"C:/data/tasking.gdb/aoi")
The following Python script demonstrates how to use the DOFToObstacleFeatures function in a stand-alone script.
# Name: DOFToObstacleFeatures_Example2.py
# Description: Convert DOF records to obstacle features.
# Import system modules
import os
import arcpy
# Set local variables
working_mty = r"d:\working\monterey"
in_dof_csv = os.path.join(working_mty, "dof.csv")
out_obx = os.path.join(working_mty, "results.gdb", "obstacles")
out_buf = os.path.join(working_mty, "results.gdb", "buffers")
aoi = os.path.join(working_mty, "tasking.gdb", "mtyAOI")
arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(32611) # UTM Zone 11
# Execute DOFToObstacleFeatures
arcpy.DOFToFeatures_intelligence(in_dof_csv, out_obx, out_buf, aoi)
Environments
Licensing information
- Basic: Requires 3D Analyst and Spatial Analyst
- Standard: Requires 3D Analyst and Spatial Analyst
- Advanced: Requires 3D Analyst and Spatial Analyst