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 en 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.

Parámetros

EtiquetaExplicaciónTipo de datos
Dataset LAS de entrada

El dataset LAS que se va a procesar.

LAS Dataset Layer
Polígonos de entrada

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

Feature Layer
Propiedad de salida

Las propiedades que se calcularán.

  • Z mínimaEl valor Z más bajo de los puntos LAS que se superponen a la entidad.
  • Z máximaEl valor Z más alto de los puntos LAS que se superponen a la entidad.
  • Z promedioEl valor Z medio de los puntos LAS que se superponen a la entidad.
  • Recuento de puntos LASLa suma de los puntos LAS que se han evaluado.
  • Desviación estándarLa desviación estándar de los valores Z.
String

Salida derivada

EtiquetaExplicaciónTipo de datos
Polígonos de entrada actualizados

Las entidades poligonales de entrada actualizadas.

Clase de entidad

arcpy.ddd.LasPointStatsByArea(in_las_dataset, in_features, out_property)
NombreExplicació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