Resumen
Calcula estadísticas sobre la distribución de mediciones de elevación de puntos de vegetación capturados en datos LAS.
Uso
Plantéese usar esta herramienta para calcular la biomasa en capturas de LIDAR de áreas boscosas.
Esta herramienta requiere la presencia de puntos clasificados del suelo y procesa un punto LAS con un valor de código de clase de 0, 1, 3, 4 y 5.
Sintaxis
LasHeightMetrics_3d (in_las_dataset, out_location, {base_name}, {statistics}, {height_percentiles}, min_height, min_points, cell_size)| Parámetro | Explicación | Tipo de datos |
in_las_dataset | El dataset LAS que se va a procesar. | LAS Dataset Layer |
out_location | La carpeta o geodatabase en la que residirán los datasets ráster de salida. | Workspace |
base_name (Opcional) | El nombre de base de los datasets ráster de salida. | String |
statistics [statistics,...] (Opcional) | Las estadísticas calculadas para los puntos no clasificados y de vegetación sobre el suelo que se encuentran dentro del área de cada celda del ráster de salida.
| String |
height_percentiles [height_percentiles,...] (Opcional) | Los percentiles de altura sobre el suelo para los que se calculará el histograma. | Long |
min_height | La altura más baja sobre el suelo de los puntos que se evaluarán. | Linear Unit |
min_points | El número mínimo de puntos que debe haber en una celda para calcular métricas de altura LAS. Las celdas con un número de puntos menor que el mínimo especificado no tendrán datos en la salida. | Long |
cell_size | El tamaño de celda para los datasets ráster de salida. | Linear Unit |
Salida derivada
| Nombre | Explicación | Tipo de datos |
| derived_out_location | La ubicación de los datasets ráster de salida. | Espacio de trabajo |
| output_rasters | Los datasets ráster de salida. |
Muestra de código
Ejemplo 1 de LasHeightMetrics (ventana de Python)
En el siguiente ejemplo se muestra cómo usar esta herramienta en la ventana de Python.
arcpy.env.workspace = 'C:/data'
arcpy.ClassifyLasBuilding_3d('Highland.lasd', minHeight='9 feet',
minArea='30 Square Feet', compute_stats=True)
Ejemplo 2 de LasHeightMetrics (script independiente)
En el siguiente ejemplo se muestra cómo usar esta herramienta en un script independiente de Python.
'''****************************************************************************
Name: Tile & Classify LAS Files
Description: Creates & classifies tiled LAS files.
****************************************************************************'''
# Import system modules
import arcpy
import tempfile
import math
in_las = arcpy.GetParameterAsText(1) # The LAS files that need to be tiled
out_folder = arcpy.GetParameterAsText(2) # folder for LAS files
basename = arcpy.GetParameterAsText(3) # basename for output files
out_lasd = arcpy.GetParameterAsText(4) # output LAS dataset
try:
# Create temp LAS dataset to reference LAS files that will be tiled
temp_lasd = arcpy.CreateUniqueName('temp.lasd', tempfile.gettempdir())
arcpy.management.CreateLasDataset(in_las, temp_lasd)
arcpy.ddd.TileLas(temp_lasd, out_folder, basename, out_lasd, las_version=1.4,
point_format=7, file_size=300)
arcpy.management.Delete(temp_lasd)
arcpy.ddd.ClassifyLasGround(out_lasd, method='AGGRESSIVE')
arcpy.ddd.ClassifyLasBuilding(out_lasd, min_height='3 Meters', min_area='4 Meters')
arcpy.ddd.ClassifyLasByHeight(out_lasd, height_classification=[(3, 6), (4,20), (5,70)],
noise='All Noise', compute_stats='COMPUTE_STATS')
except arcpy.ExecuteError:
print(arcpy.GetMessages())
Entornos
Información sobre licencias
- ArcGIS Desktop Basic: Requiere 3D Analyst
- ArcGIS Desktop Standard: Requiere 3D Analyst
- ArcGIS Desktop Advanced: Requiere 3D Analyst