Estadísticas de punto de LAS por área (3D Analyst)

Resumen

Evalúa las estadísticas de los puntos LAS que se solapan con el área definida por entidades de polígono.

Uso

  • Es posible hacer que la capa del dataset LAS limite los puntos LAS a mostrar y procesar seleccionando cualquier combinación de códigos de clasificación, marcadores de clasificación y valores de retorno en la configuración del filtro de la capa. Los filtros se pueden definir en el cuadro de diálogo Propiedades de capa o con la herramienta Crear capa de dataset LAS.

  • Plantéese la posibilidad de utilizar esta herramienta para determinar las características de los puntos LIDAR sobre regiones de interés definidas por los límites de los polígonos. Solo se evaluará la elevación obtenida a partir de los puntos LAS. Se ignorará cualquier valor Z aportado por las entidades de restricción de superficie a las que se hace referencia en el dataset LAS.

  • Esta herramienta agrega varios campos a la tabla de atributos de las entidades de entrada. Si los campos ya existen, los valores se sobrescribirán. Los valores de los campos actualizados reflejarán las estadísticas indicadas sobre los puntos LAS que se superponen a la entidad de entrada.

    • Z_MIN: valor Z más bajo.
    • Z_MAX: valor Z más alto.
    • Z_MEAN: valor Z medio.
    • POINT_COUNT: la suma de los puntos LAS que se han evaluado.
    • STD_DEV: la desviación estándar de los valores z.

Sintaxis

arcpy.3d.LasPointStatsByArea(in_las_dataset, in_features, out_property)
ParámetroExplicaciónTipo de datos
in_las_dataset

El dataset LAS que se va a procesar.

LAS Dataset Layer
in_features

El polígono que define el área cuyas estadísticas se van a notificar.

Feature Layer
out_property
[out_property,...]

Las propiedades que se calcularán.

  • Z_MINEl valor Z más bajo de los puntos LAS que se superponen a la entidad.
  • Z_MAXEl valor Z más alto de los puntos LAS que se superponen a la entidad.
  • Z_MEANEl valor Z medio de los puntos LAS que se superponen a la entidad.
  • POINT_COUNTLa suma de los puntos LAS que se han evaluado.
  • STANDARD_DEVIATIONLa desviación estándar de los valores Z.
String

Salida derivada

NombreExplicaciónTipo de datos
output_polygons

Las entidades poligonales de entrada actualizadas.

Clase de entidad

Muestra de código

Ejemplo 1 de LasPointStatsByArea (ventana de Python)

En el siguiente ejemplo se muestra cómo usar esta herramienta en la ventana de Python.

import arcpy
from arcpy import env

env.workspace = 'C:/data'
arcpy.ddd.LasPointStatsByArea('city_lidar.lasd', 'study_area.shp', 
                              ['Z_MIN', 'Z_MAX'])
Ejemplo 2 de LasPointStatsByArea (script independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en un script independiente de Python.

'''****************************************************************************
       Name: Extract Building Footprints
Description: Extract footprint from lidar points classified as buildings, 
             regularize its geometry, and calculate the building height.

****************************************************************************'''
import arcpy

lasd = arcpy.GetParameterAsText(0)
footprint = arcpy.GetParameterAsText(1)

try:
    lasd_layer = 'building points'
    arcpy.management.MakeLasDatasetLayer(lasd, lasd_layer, class_code=6)
    temp_raster = 'in_memory/bldg_raster'
    arcpy.management.LasPointStatsAsRaster(lasd_layer, temp_raster,
                                           'PREDOMINANT_CLASS', 'CELLSIZE', 2.5)
    temp_footprint = 'in_memory/footprint'
    arcpy.conversion.RasterToPolygon(temp_raster, temp_footprint)
    arcpy.ddd.RegularizeBuildingFootprint(temp_footprint, footprint, 
                                          method='RIGHT_ANGLES')
    arcpy.ddd.LasPointStatsByArea(lasd_layer, footprint, ['MIN_Z', 'MAX_Z'])
    arcpy.management.AddField(footprint, 'Height', 'Double')
    arcpy.management.CalculateField(footprint, 'Height', 
                                    "round('!Z_Max! - !Z_Min!', 2)", 
                                    'PYTHON_9.3')


except arcpy.ExecuteError:
    print(arcpy.GetMessages())

Información de licenciamiento

  • Basic: Requiere 3D Analyst
  • Standard: Requiere 3D Analyst
  • Advanced: Requiere 3D Analyst

Temas relacionados