Resumen
Clasifica puntos de tejados de edificios en datos LIDAR aéreos.
Ilustración
Uso
Los datos LAS deben tener puntos de suelo clasificados antes de clasificar los puntos de tejado de los edificios. Considere el uso de la herramienta Clasificar el terreno LAS si los puntos de suelo no se han clasificado. Los puntos de suelo deben tener un valor de código de clase de 2, de modo que, si los puntos de suelo tienen un valor de código de clase diferente, utilice la herramienta Cambiar los códigos de clase LAS para reasignar los códigos de clase según corresponda.
Los puntos LAS con valores de código de clase de 0 y 1 se evaluarán para determinar si encajan en la altura y área mínimas especificadas para identificar tejados de edificios. Los puntos de edificio existentes clasificados con un valor de código de clase de 6 también se utilizarán para determinar si los puntos adyacentes encajan en un plano que conecte la misma estructura. Este proceso utiliza un algoritmo RANSAC optimizado que comienza con un punto de semilla seleccionado al azar, que podría hacer que los resultados de la clasificación variasen ligeramente entre ejecuciones.
Si se clasifican por error pequeños clústeres de puntos que no pertenecen a los edificios, plantéese aumentar el parámetro Área mínima para eliminar la inclusión de tales puntos. Por el contrario, si parece que la herramienta ha olvidado los edificios, plantéese utilizar un área más pequeña para capturar las regiones que faltan. Se puede proporcionar una extensión de procesamiento o un polígono de límite para limitar aún más la región que se está procesando, lo que puede ser útil al probar los resultados de los valores experimentales.
La herramienta es adecuada para el procesamiento de levantamientos topográficos LIDAR aéreos. Es posible que los archivos LAS creados a partir de orígenes como nubes de puntos fotogramétricas produzcan resultados menos fiables debido a que la probabilidad de puntos de ruido es mayor.
Esta herramienta funciona mejor si los puntos LAS tienen una densidad de puntos uniforme. Por lo general, los levantamientos topográficos LIDAR aéreos tendrán algo de superposición en su cobertura, lo que producirá áreas con clústeres de puntos más densos. La inclusión de esos puntos distorsiona el espaciado nominal de los puntos, además de producir errores de comisión en las regiones superpuestas y errores de omisión en el área no superpuesta. Para obtener mejores resultados, utilice la herramienta Clasificar solapamiento de LAS para clasificar puntos superpuestos y modifique el valor Espaciado de punto promedio para ofrecer una estimación del espaciado de puntos nominal fuera de las regiones de superposición.
Sintaxis
ClassifyLasBuilding_3d (in_las_dataset, {min_height}, min_area, {compute_stats}, {extent}, boundary, {process_entire_files}, point_spacing)| Parámetro | Explicación | Tipo de datos |
in_las_dataset | El dataset LAS que se va a clasificar. | LAS Dataset Layer |
min_height (Opcional) | La altura desde el suelo que define el punto más bajo a partir del cual se identificarán los puntos del tejado. | Linear Unit |
min_area | El área más pequeña de puntos coplanares que se usará para establecer la presencia del tejado de un edificio. | Areal Unit |
compute_stats (Opcional) | Especifica si las estadísticas se deben calcular para los archivos LAS a los que hace referencia el dataset LAS. La existencia de las estadísticas permite las opciones de simbología y filtrado de la capa del dataset LAS para mostrar solo los valores de atributos LAS que existen en los archivos LAS.
| Boolean |
extent (Opcional) | Especifica la extensión de los datos que se evaluarán con esta herramienta. | Extent |
boundary | Entidad poligonal que define el área de interés que se va a procesar con esta herramienta. | Feature Layer |
process_entire_files (Opcional) | Especifica cómo se va a aplicar la extensión de procesamiento.
| Boolean |
point_spacing | El espaciado promedio de los puntos LAS que se usará para determinar la presencia de tejados de edificios. Este valor será de forma predeterminada el espacio promedio calculado para el dataset LAS. | Linear Unit |
Salida derivada
| Nombre | Explicación | Tipo de datos |
| derived_las_dataset | El dataset LAS clasificado para los tejados de los edificios. | Capa de dataset LAS |
Muestra de código
Ejemplo 1 de ClassifyLasBuilding (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 ClassifyLasBuilding (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