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 en 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, puede especificar el valor del parámetro Tipo de datos de salida como un entero para ahorrar espacio en disco si los requisitos de precisión de los valores z 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.

Parámetros

EtiquetaExplicaciónTipo de datos
Dataset LAS de entrada

El dataset LAS que se procesará.

LAS Dataset Layer
Ráster de salida

Ubicación y nombre del ráster de salida. Cuando almacene un dataset ráster en una geodatabase o en una carpeta como una cuadrícula de Esri, no agregue ninguna 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 .tif o en una geodatabase, el tipo de compresión de ráster y la calidad se pueden especificar utilizando la configuración del entorno de geoprocesamiento.

Raster Dataset
Campo de valor
(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
Tipo de interpolación
(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
Tipo de datos de salida
(Opcional)

Especifica el tipo de valores numéricos que se almacenarán en el ráster de salida.

  • FlotanteEl ráster de salida utilizará un punto flotante de 32 bits, que admite valores que oscilan de -3,402823466e+38 a 3,402823466e+38. Esta es la opción predeterminada.
  • EnteroEl 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
Tipo de muestra
(Opcional)

Especifica cómo se interpretará el parámetro Valor de muestreo para definir el tamaño de celda del ráster de salida.

  • ObservacionesEl Valor de muestreo definirá el número de columnas o filas en el ráster de salida en función del que sea más largo. El tamaño de celda se obtendrá dividiendo el lado más largo de la extensión de la salida con la entrada en el parámetro Valor de muestreo. Si se utiliza un valor de observación de 3000 en un dataset cuyo lado más largo es de 23,67 kilómetros, la resolución del ráster de salida será de 7,89 metros. Este método ofrece una forma útil de crear una salida con un tamaño predecible que pueda generarse rápidamente.
  • Tamaño de celdaEl tamaño de celda se definirá directamente mediante el parámetro Valor de muestreo. Esta es la opción predeterminada.
String
Valor de muestra
(Opcional)

El valor utilizado junto con el parámetro Tipo de muestreo para definir el tamaño de celda del ráster de salida.

Double
Factor Z
(Opcional)

El factor por el que se multiplicarán los valores z. Generalmente se utiliza para convertir unidades lineales para que coincidan con las unidades lineales x,y. El valor predeterminado es 1, que no altera los valores de elevación. Este parámetro no está disponible si la referencia espacial de la superficie de entrada tiene un datum z con una unidad lineal especificada.

Double

arcpy.conversion.LasDatasetToRaster(in_las_dataset, out_raster, {value_field}, {interpolation_type}, {data_type}, {sampling_type}, {sampling_value}, {z_factor})
NombreExplicaciónTipo de datos
in_las_dataset

El dataset LAS que se procesará.

LAS Dataset Layer
out_raster

Ubicación y nombre del ráster de salida. Cuando almacene un dataset ráster en una geodatabase o en una carpeta como una cuadrícula de Esri, no agregue ninguna 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 .tif o en una geodatabase, el 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 que se almacenarán en el ráster de salida.

  • FLOATEl ráster de salida utilizará un punto flotante de 32 bits, que admite valores que oscilan 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 cómo se interpretará el parámetro Valor de muestreo para definir el tamaño de celda del ráster de salida.

  • OBSERVATIONSEl Valor de muestreo definirá el número de columnas o filas en el ráster de salida en función del que sea más largo. El tamaño de celda se obtendrá dividiendo el lado más largo de la extensión de la salida con la entrada en el parámetro Valor de muestreo. Si se utiliza un valor de observación de 3000 en un dataset cuyo lado más largo es de 23,67 kilómetros, la resolución del ráster de salida será de 7,89 metros. Este método ofrece una forma útil de crear una salida con un tamaño predecible que pueda generarse rápidamente.
  • CELLSIZEEl tamaño de celda se definirá directamente mediante el parámetro Valor de muestreo. Esta es la opción predeterminada.
String
sampling_value
(Opcional)

El valor utilizado junto con el parámetro Tipo de muestreo para definir el tamaño de celda del ráster de salida.

Double
z_factor
(Opcional)

El factor por el que se multiplicarán los valores z. Generalmente se utiliza para convertir unidades lineales para que coincidan con las unidades lineales x,y. El valor predeterminado es 1, que no altera los valores de elevación. Este parámetro no está disponible 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.ddd.LasDatasetToRaster('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 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

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)

Temas relacionados