ArcMap

Estadísticas de puntos LAS por área

Resumen

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

Uso

  • Cuando un dataset LAS se especifica como entrada, todos los puntos de datos en los archivos LAS a los que hace referencia se procesarán. También es posible seleccionar un subconjunto de datos LIDAR por los códigos de clasificación, indicadores de clasificación y valores de retoro aplicando los filtros de punto LAS que desee mediante una capa de dataset LAS. Los filtros se pueden definir mediante el cuadro de diálogo de propiedades de la capa 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: Suma de los puntos LAS que se han evaluado.
    • STD_DEV: la desviación estándar de los valores Z.

Sintaxis

LasPointStatsByArea_3d (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

Las propiedades que se calcularán.

  • Z_MIN —El valor Z más bajo de los puntos LAS que se superponen a la entidad.
  • Z_MAX —El valor Z más alto de los puntos LAS que se superponen a la entidad.
  • Z_MEAN —El valor Z medio de los puntos LAS que se superponen a la entidad.
  • POINT_COUNT —La suma de los puntos LAS que se han evaluado.
  • STANDARD_DEVIATION —La desviación estándar de los valores Z.
String

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 (secuencia de comandos independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en una secuencia de comandos 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())

Entornos

Información sobre licencias

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

Temas relacionados