ArcMap

De LAS a multipunto

Resumen

Crea entidades multipunto que utilizan uno o más archivos LIDAR.

Ilustración

De LAS a multipunto

Uso

  • Esta herramienta admite las versiones de archivos ZLAS y LAS de la 1.0 a la 1.4.

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

  • El formato LAS admite el almacenamiento de numerosos atributos predefinidos para cada punto LIDAR. Si no sabe qué atributos están disponibles en un determinado grupo de archivos LAS, plantéese la posibilidad de utilizar un dataset LAS para revisar las propiedades del archivo LAS.

  • Al cargar varios atributos LAS en una base de datos de Oracle, tendrá que asegurarse de que todas las palabras clave DBTUNE del parámetro attribute_binary estén establecidas para utilizar objetos binarios grandes (BLOB) y no LONGRAW. Esto se debe a que los atributos LAS se cargan como BLOBs y Oracle no admite varios BLOBs en las tablas LONGRAW. Consulte con el administrador de la base de datos de Oracle para obtener más ayuda.

Sintaxis

arcpy.ddd.LASToMultipoint(input, out_feature_class, average_point_spacing, {class_code}, {return}, {attribute}, {input_coordinate_system}, {file_suffix}, {z_factor}, {folder_recursion})
ParámetroExplicaciónTipo de datos
input
[input,...]

Los archivos LAS o ZLAS que se importarán a una clase de entidad multipunto. Si se especifica una carpeta, se importarán todos los archivos LAS incluidos en ella.

En la ventana del cuadro de diálogo de la herramienta, también se puede especificar una carpeta como entrada seleccionando la carpeta deseada en el Explorador de Windows y arrastrándola al cuadro de entrada del parámetro.

Folder; File
out_feature_class

La clase de entidad que generará esta herramienta.

Feature Class
average_point_spacing

La distancia 2D promedio entre puntos en el archivo o los archivos de entrada. Esto puede ser una aproximación. Si las áreas se han muestreado en diferentes densidades, especifique el espaciado más pequeño. El valor se debe proporcionar en las unidades de proyección del sistema de coordenadas de salida.

Double
class_code
[class_code,...]
(Opcional)

Los códigos de clasificación para utilizar como un filtro de consulta para los puntos de datos LAS. Los valores válidos varían de 1 a 32. De manera predeterminada, no se aplica ningún filtro.

Long
return
[return,...]
(Opcional)

Los valores de retorno que se utilizarán para filtrar los puntos LAS que se importan a entidades multipunto.

  • ANY_RETURNS —Cualquier valor de retorno
  • 1 —1
  • 2 —2
  • 3 —3
  • 4 —4
  • 5 —5
  • 6 —6
  • 7 —7
  • 8 —8
  • LAST_RETURNS —Retornos más recientes
String
attribute
[[keyword, name],...]
(Opcional)

Las propiedades de los puntos LAS cuyos valores se almacenarán en campos de objeto binario grande (BLOB) en la tabla de atributos de la salida. Si las entidades resultantes van a participar en un dataset de terreno, también se pueden utilizar los atributos almacenados para simbolizar el terreno. La columna Nombre indica el nombre del campo que se utilizará para almacenar los atributos especificados. Se admiten las siguientes propiedades de LAS:

  • INTENSITY —Intensidad
  • RETURN_NUMBER —Número de retorno
  • NUMBER_OF_RETURNS —Número de retornos por pulso
  • SCAN_DIRECTION_FLAG —Marcador de dirección de escaneo
  • EDGE_OF_FLIGHTLINE —Borde de la línea de vuelo
  • CLASSIFICATION —Clasificación
  • SCAN_ANGLE_RANK —Clasificación de ángulo de escaneo
  • FILE_MARKER —Marcador de archivos
  • USER_BIT_FIELD —Valor de datos de usuario
  • GPS_TIME —Tiempo del GPS
  • COLOR_RED —Banda roja
  • COLOR_GREEN —Banda verde
  • COLOR_BLUE —Banda azul
Value Table
input_coordinate_system
(Opcional)

El sistema de coordenadas del archivo LAS de entrada

Coordinate System
file_suffix
(Opcional)

El sufijo de los archivos para importar desde una carpeta de entrada. Este parámetro es necesario cuando se especifica una carpeta como entrada.

String
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
folder_recursion
(Opcional)

Escanea todas las subcarpetas cuando se selecciona una carpeta de entrada que contenga datos en un directorio de subcarpetas. La clase de entidad de salida se generará con una fila para cada archivo encontrado en la estructura del directorio.

  • NO_RECURSION —solo los archivos LAS que se encuentran en una carpeta de entrada se convertirán a entidades de multipunto. Esta es la opción predeterminada.
  • RECURSION —Todos los archivos LAS que residen en los subdirectorios de una carpeta de entrada se convertirán a entidades de multipunto.
Boolean

Muestra de código

Ejemplo 1 de LASToMultipoint (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

arcpy.CheckOutExtension("3D")
env.workspace = "C:/data"
arcpy.LASToMultipoint_3d("001.las", "Test.gdb/feature_dataset/sample_1", 1.5, 
                        "2", "ANY_RETURNS", "INTENSITY", "Coordinate Systems"\
                        "/Projected Coordinate Systems/UTM/NAD 1983/NAD 1983 "\
                        "UTM Zone 17N.prj", "las", 1)
Ejemplo 2 de LASToMultipoint (script independiente)

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

'''****************************************************************************
Name: Define Data Boundary of LAS File
Description: This script demonstrates how to delineate data boundaries of 
             LAS files with irregularly clustered points. It is intended for 
             use as a script tool with one input LAS file.
****************************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback

# Set local variables
inLas = arcpy.GetParameterAsText(0) #input LAS file
ptSpacing = arcpy.GetParameterAsText(1) # LAS point spacing
classCode = arcpy.GetParameterAsText(2) # List of integers
returnValue = arcpy.GetParameterAsText(3) # List of strings
outTin = arcpy.GetParameterAsText(4) # TIN created to delineate data area
outBoundary = arcpy.GetParameterAsText(5) # Polygon boundary file

try:
    arcpy.CheckOutExtension("3D")
    # Execute LASToMultipoint
    arcpy.AddMessage("Creating multipoint features from LAS...")
    lasMP = arcpy.CreateUniqueName('lasMultipoint', 'in_memory')
    arcpy.ddd.LASToMultipoint(inLas, LasMP, ptSpacing, class_code, 
                             "ANY_RETURNS", "", sr, inFormat, zfactor)
    # Execute CreateTin
    arcpy.AddMessage("Creating TIN dataset...")
    arcpy.ddd.CreateTin(outTin, sr, "{0} Shape.Z masspoints"\
                       .format(lasMP), "Delaunay")
    # Execute CopyTin
    arcpy.AddMessage("Copying TIN to delineate data boundary...")
    arcpy.ddd.CopyTin(outTin, "{0}_copy".format(outTin))
    # Execute DelineateTinDataArea
    arcpy.AddMessage("Delineating TIN boundary...")
    maxEdge = ptSpacing * 4
    arcpy.ddd.DelineateTinDataArea(outTin, maxEdge, "PERIMETER_ONLY")
    # Execute TinDomain
    arcpy.AddMessage("Exporting data area to polygon boundary...")
    arcpy.ddd.TinDomain(outTin, outBoundary, "POLYGON")
    arcpy.AddMessage("Finished")
    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

Información de licenciamiento

  • Basic: Requiere 3D Analyst
  • Standard: Requiere 3D Analyst
  • Advanced: Requiere 3D Analyst

Temas relacionados