De dataset LAS a ráster (Conversión)

Resumen

Crea un ráster utilizando valores de elevación, intensidad o RGB almacenados en los puntos LIDAR a los que hace referencia el dataset LAS.

Ilustración

De dataset LAS a ráster

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 con la herramienta Crear capa de dataset LAS.

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

  • Al exportar un ráster grande, plantéese especificar el Tipo de datos de salida como un entero para ahorrar espacio en disco si los requisitos de precisión de sus valores z son tales que se pueden representar con datos enteros.

  • La inclusión de un límite del área de estudio como una restricción de recorte en la definición de la entrada dataset LAS es muy recomendable. Una razón es para evitar que se ocurra interpolación fuera de la extensión de datos reales de la topografía. En segundo lugar, puede haber una grave penalización de rendimiento al utilizar las opciones de vecinos naturales si el área de datos no está definida correctamente.

    Nota:

    Cuando se utiliza BINNING, solo se respetan las restricciones de recorte, borrado y reemplazo. Las líneas de corte y los puntos de anclaje no. La opción de triangulación respeta todos los tipos de restricciones, pero tarda más tiempo en ejecutarse.

Sintaxis

arcpy.conversion.LasDatasetToRaster(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 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
value_field
(Opcional)

Los datos LIDAR que se utilizarán para generar la salida de 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.
  • RGBLos valores RGB de los puntos LIDAR se utilizarán para crear imágenes de tres bandas.
String
interpolation_type
"BINNING {cell_assignment_type} {void_fill_method}" or "TRIANGULATION {interpolation_method} {point_thinning_type} {point_selection_method} {resolution}"
(Opcional)

Técnica de interpolación que se usará para determinar los valores de celda del ráster de salida.

El método de bin proporciona un Método de asignación de celdas para determinar cada celda de salida usando los puntos contenidos en su extensión, junto con un Método de relleno de vacíos para determinar el valor de las celdas que no contienen ningún punto LAS.

Métodos de asignación de celdas

  • AVERAGE: asigna el valor promedio de todos los puntos dentro de la celda. Esta es la opción predeterminada.
  • MINIMUM: asigna el valor mínimo de los puntos que se encuentran dentro de la celda.
  • MAXIMUM: 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étodos de relleno de vacíos

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

Los métodos de interpolación por triangulación derivan valores de celdas usando un método basado en TIN y ofrecen además la oportunidad de acortar el tiempo de procesamiento reduciendo el muestreo de datos LAS por medio de la técnica Tamaño de la ventana.

Métodos de triangulación

  • Linear: utiliza la interpolación lineal para determinar los valores de celda.
  • Natural Neighbors: utiliza la interpolación de vecinos naturales para determinar el valor de celda.

Métodos de selección de tamaño de ventana

  • 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.
Interpolate
data_type
(Opcional)

Especifica el tipo de valores numéricos almacenados en el ráster de salida.

  • 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 redondea los valores z al número entero más cercano y escribe el número entero en 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 z. Esto se utiliza generalmente para convertir las unidades lineales Z para que coincidan con las unidades lineales XY. El valor predeterminado es 1, que no altera los valores de elevación. Este parámetro está deshabilitado si la referencia espacial de la superficie de entrada tiene un datum Z con una unidad lineal especificada.

Double

Muestra de código

Ejemplo 1 de LasDatasetToRaster (ventana de Python)

En el siguiente ejemplo se muestra cómo usar esta herramienta en la ventana de Python.

arcpy.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)

En el siguiente ejemplo se muestra cómo usar esta herramienta en un script 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

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)

    # 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)
    print(arcpy.GetMessages())

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

except Exception as err:
    print(err.args[0])

finally:
    arcpy.management.Delete(lasLyr)

Información de licenciamiento

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

Temas relacionados