Estadística de puntos LAS como ráster (Administración de datos)

Nivel de licencia:BasicStandardAdvanced

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

Sintaxis

LasPointStatsAsRaster_management (in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
ParámetroExplicaciónTipo de datos
in_las_dataset

El dataset de LAS de entrada.

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. Cuando se almacena el ráster en una carpeta, se puede proporcionar una extensión de archivo para definir el formato del ráster:

  • .bil—Esri BIL
  • .bip—Esri BIP
  • .bsq—Esri BSQ
  • .dat: ENVI DAT
  • .img—ERDAS IMAGINE
  • .png: PNG
  • .tif: TIFF

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_COUNTLa cantidad de puntos de la última devolución.
  • POINT_COUNTLa cantidad de puntos de todas las devoluciones.
  • PREDOMINANT_LAST_RETURNEl último valor de devolución más frecuente.
  • PREDOMINANT_CLASSCódigo de clase más frecuente.
  • INTENSITY_RANGEEl rango de valores de intensidad.
  • Z_RANGEEl 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.

  • OBSERVATIONSDefine la cantidad de celdas que dividen el lado más largo de la extensión del dataset LAS.
  • CELLSIZEDefine 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

Ejemplo de código

Ejemplo 1 de LasPointStatsAsRaster (ventana de Python)

El siguiente ejemplo muestra cómo usar de esta herramienta en la ventana 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)

El siguiente ejemplo muestra cómo usar esta herramienta en una secuencia de comandos 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
import exceptions, sys, traceback

# 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)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
    # Execute LasPointStatsAsRaster
    arcpy.management.LasPointStatsAsRaster(lasD, statsRaster, 
                                           "PREDOMINANT_RETURNS_PER_PULSE", 
                                           "CELLSIZE", 15)
    arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
    print arcpy.GetMessages()
except:
    # Get the traceback object
    tb = sys.exc_info()[2]
    tbinfo = traceback.format_tb(tb)[0]
    # Concatenate error information into message string
    pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
          .format(tbinfo, str(sys.exc_info()[1]))
    msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
    # Return python error messages for script tool or Python Window
    arcpy.AddError(pymsg)
    arcpy.AddError(msgs)

Entornos

Temas relacionados

Información sobre licencias

ArcGIS for Desktop Basic: Requiere 3D Analyst o Spatial Analyst
ArcGIS for Desktop Standard: Sí
ArcGIS for Desktop Advanced: Sí
5/12/2014