Estadísticas de punto de LAS como ráster (Administración de datos)

Resumen

Crea un ráster cuyos valores de celda reflejan información estadística sobre las mediciones de los archivos LAS a los que hace referencia un dataset LAS.

Uso

• La capa de dataset LAS también se puede utilizar para controlar la aplicación de las funciones de restricción de superficie a las que un dataset LAS puede hacer referencia. Las restricciones se aplican cuando el dataset LAS se muestra o procesa como una superficie triangulada.

• El formato LAS admite la clasificación de cada punto según las especificaciones definidas por la American Society for Photogrammetry and Remote Sensing (ASPRS, por sus siglas en inglés o Sociedad americana de fotogrametría y teledetección). La plataforma ArcGIS aplica el esquema de clasificación especificado para la versión 1.4 de los archivos LAS:

  Valor de clasificación	Tipo de clasificación
  0	Nunca clasificado
  1	No asignado
  2	Suelo
  3	Vegetación baja
  4	Vegetación media
  5	Vegetación alta
  6	Edificio
  7	Ruido bajo
  8	Clave de modelo/Reservada
  9	Water
  10	Ferrocarril
  11	Superficie de la carretera
  12	Superposición/Reservado
  13	Protector de cable
  14	Conductor de cable
  15	Torre de transmisión
  16	Cable - Conector
  17	Plataforma del puente
  18	Ruido alto
  19 – 63	Reservado para la definición de la ASPRS (las versiones de LAS de la 1.1 a la 1.3 admiten hasta el código de clase 31)
  32 – 255	Lo puede definir el usuario (solo se admite en la versión de LAS 1.0 y en algunas versiones de 1.4)

  Nota:

  Aunque el volumen de las nuevas asignaciones de código de clase que se introdujo con la versión de LAS 1.4 en principio se designó como Reservado, los códigos de clase 8 y 12 cambiaron de Clave de modelo y Superposición a Reservados.

• Considere utilizar PREDOMINANT_LAST_RETURN para el parámetro Método para identificar las ubicaciones con valores de devolución más altos que podrían indicar la presencia de vegetación.

• Considere utilizar Z_RANGE para el parámetro Método para determinar las ubicaciones con valores atípicos potenciales.

Sintaxis

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})

Parámetro	Explicación	Tipo de datos
in_las_dataset	El dataset LAS que se va a procesar.	LAS Dataset Layer
out_raster	Ubicación y nombre del ráster de salida. Al almacenar un dataset ráster en una geodatabase o en una carpeta como una cuadrícula de Esri, no se debe adicionar extensión de archivo al nombre del dataset ráster. Se puede indicar una extensión de archivo para definir el formato del ráster al almacenarlo en una carpeta, como .tif para generar un GeoTIFF o .img para generar un archivo con el formato ERDAS IMAGINE. Si el ráster se almacena como un archivo TIFF o en una geodatabase, su tipo de compresión de ráster y la calidad se pueden especificar utilizando la configuración del entorno de geoprocesamiento.	Raster Dataset
method (Opcional)	El tipo de estadísticas que se recopiló sobre los puntos LAS en cada celda del ráster de salida. PULSE_COUNT —La cantidad de puntos de la última devolución. POINT_COUNT —La cantidad de puntos de todas las devoluciones. PREDOMINANT_LAST_RETURN —El último valor de devolución más frecuente. PREDOMINANT_CLASS —Código de clase más frecuente. INTENSITY_RANGE —El rango de valores de intensidad. Z_RANGE —El rango de valores de elevación.	String
sampling_type (Opcional)	Especifica el método utilizado para interpretar el Valor de muestreo para definir la resolución del ráster de salida. OBSERVATIONS —Define la cantidad de celdas que dividen el lado más largo de la extensión del dataset LAS. CELLSIZE —Define el tamaño de celda del ráster de salida. Esta es la opción predeterminada.	String
sampling_value (Opcional)	Especifica el valor utilizado junto con el Tipo de muestreo para definir la resolución del ráster de salida.	Double

Muestra de código

Ejemplo 1 de LasPointStatsAsRaster (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.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)

Ejemplo 2 de LasPointStatsAsRaster (secuencia de comandos independiente)

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

'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for
             each pulse in a given set of LAS files.
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)

# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
                                       "PREDOMINANT_RETURNS_PER_PULSE",
                                       "CELLSIZE", 15)

Entornos

Espacio de trabajo actual, Extensión, Sistema de coordenadas de salida, Transformaciones geográficas, Tamaño de celda, Pirámide, Estadísticas de ráster, Palabra clave CONFIG de salida, Auto ejecución

Información de licenciamiento

• Basic: Requiere 3D Analyst o Spatial Analyst
• Standard: Sí
• Advanced: Sí