Dataset LAS para ráster (Conversión)

Nivel de licencia:BasicStandardAdvanced

Resumen

Crea un ráster utilizando valores de elevación, intensidad o RGB almacenados en los archivos lidar (*.las) a los que hace referencia el dataset LAS.

Ilustración

LAS Dataset to Raster

Uso

Sintaxis

LasDatasetToRaster_conversion (in_las_dataset, out_raster, {value_field}, {interpolation_type}, {data_type}, {sampling_type}, {sampling_value}, {z_factor})
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
value_field
(Opcional)

Especifica la información lidar que se utilizará para generar la salida del ráster.

  • ELEVATIONLa elevación desde los archivos lidar se utilizará para crear el ráster. Esta es la opción predeterminada.
  • INTENSITYLa información de intensidad a partir de los archivos lidar se utilizará para crear el ráster.
  • RGBLas imágenes derivadas de los valores RGB incorporados con los puntos lidar se utilizarán para crear el ráster.
String
interpolation_type
"BINNING {cell_assignment_type} {void_fill_method}" or "TRIANGULATION {interpolation_method} {point_thinning_type} {point_selection_method} {resolution}"
(Opcional)

El método de interpolación que se usa para producir el ráster.

BINNING: Los valores de las celdas se obtienen utilizando los puntos que caen dentro de la extensión de la celda, excepto las celdas que no contienen puntos en su extensión. Las siguientes opciones están disponibles para esta técnica:

  • Tipo de asignación de celda: método utilizado para definir el valor de cualquiera de las celdas que contenga puntos dentro de su extensión.
    • PROMEDIO: asigna el valor promedio de todos los puntos dentro de la celda. Esta es la opción predeterminada.
    • MÍNIMO: asigna el valor mínimo de los puntos que se encuentran dentro de la celda.
    • MÍNIMO: asigna el valor máximo de los puntos que se encuentran dentro de la celda.
    • IDW: utiliza la interpolación ponderada por distancia inversa para determinar el valor de la celda.
    • NEAREST: utiliza la asignación de vecino más cercano para determinar el valor de celda.
  • Método de relleno de vacío: método de interpolación que se utiliza para definir los valores para las celdas que no tienen puntos dentro de su extensión.
    • NONE: se asigna NoData a la celda.
    • SIMPLE: calcula el promedio de los valores de los datos de las celdas que rodean a la celda NoData para eliminar pequeños vacíos.
    • LINEAR: realiza la triangulación en las áreas vacías y utiliza la interpolación lineal en el valor de triangulación para determinar el valor de celda. Esta es la opción predeterminada.
    • NATURAL_NEIGHBOR: utiliza la interpolación de vecinos naturales para determinar el valor de celda.

TRIANGULACIÓN: se obtienen los valores de las celdas interpolando las mediciones desde una representación de triángulos del dataset LAS. Las siguientes opciones están disponibles para esta técnica:

  • Método de interpolación: método de interpolación que define los valores de celda:
    • LINEAL: utiliza la interpolación lineal contra la superficie de dataset LAS triangulada para determinar el valor de celda.
    • NATURAL_NEIGHBOR: utiliza la interpolación de vecinos naturales para determinar el valor de celda.
  • Tipo de simplificación de puntos: determina si los puntos de datos LAS se simplifican:
    • NONE: los puntos LAS no se simplifican. Esta es la opción predeterminada.
    • WINDOW_SIZE: los puntos LAS se simplifican al identificar el punto que satisface los criterios de selección dentro del área definida por el tamaño de la ventana.
  • Método de selección de puntos: método de selección utilizado para simplificar puntos de datos LAS al utilizar simplificación WINDOW_SIZE:
    • MAXIMUM: se mantiene el punto con el valor más alto en cada tamaño de ventana. Esta es la opción predeterminada.
    • MINIMUM: se mantiene el punto con el valor más bajo en cada tamaño de ventana.
    • CLOSEST_TO_MEAN: se mantiene el punto cuyo valor es más cercano a la media de todos los valores de los puntos en el tamaño de la ventana.
  • Resolución: un valor numérico que define el área del tamaño de ventana que se utiliza para la simplificación de puntos.
String
data_type
(Opcional)

El tipo de datos del ráster de salida se puede definir mediante las siguientes palabras clave:

  • FLOATEl ráster de salida utilizará una coma flotante de 32 bits, que admite valores que varían de -3,402823466e+38 a 3,402823466e+38. Esta es la opción predeterminada.
  • INTEl ráster de salida utilizará una profundidad de bit de enteros apropiada. Esta opción aproximará los valores z al número entero más cercano y escribirá un entero para cada valor de celda ráster.
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
z_factor
(Opcional)

El factor por el que se multiplicarán los valores de elevación. Esto se suele usar para convertir unidades lineales Z que coinciden con las de las unidades lineales XY. El valor predeterminado es 1, que no altera los valores de elevación.

Double

Ejemplo de código

Ejemplo 1 de LasDatasetToRaster (ventana de Python)

El siguiente ejemplo muestra cómo usar de esta herramienta en la ventana Python:

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.LasDatasetToRaster_3d('baltimore.lasd', 'baltimore.tif', 'INTENSITY', 
                          'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT', 
                          'CELLSIZE', 10, 3.28)
Ejemplo 2 de LasDatasetToRaster (secuencia de comandos independiente)

El siguiente ejemplo muestra cómo usar esta herramienta en una secuencia de comandos independiente de Python:

'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
             ground measurements from LAS files to a raster using a 
             LAS dataset. This sample is designed to be used as a script
             tool.
*********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback

try:
    # Set Local Variables
    inLas = arcpy.GetParameterAsText(0)
    recursion = arcpy.GetParameterAsText(1)
    surfCons = arcpy.GetParameterAsText(2)
    classCode = arcpy.GetParameterAsText(3)
    returnValue = arcpy.GetParameterAsText(4)
    spatialRef = arcpy.GetParameterAsText(5)
    lasD = arcpy.GetParameterAsText(6)
    outRaster = arcpy.GetParameterAsText(7)
    cellSize = arcpy.GetParameter(8)
    zFactor = arcpy.GetParameter(9)
    if arcpy.ProductInfo == 'ArcView':
        arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
    # Execute MakeLasDatasetLayer
    lasLyr = arcpy.CreateUniqueName('Baltimore')
    arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
    # Execute LasDatasetToRaster
    arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
                              'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                              'CELLSIZE', cellSize, zFactor)
    arcpy.GetMessages()
    
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)
    
finally:
    arcpy.management.Delete(lasLyr)

Entornos

Temas relacionados

Información sobre licencias

ArcGIS for Desktop Basic: Requiere ArcGIS 10.1 Spatial Analyst o ArcGIS 10.1 3D Analyst
ArcGIS for Desktop Standard: Requiere ArcGIS 10.1 Spatial Analyst o ArcGIS 10.1 3D Analyst
ArcGIS for Desktop Advanced: Sí
9/11/2013