De LAS a multipunto—Ayuda | ArcGIS for Desktop

Resumen
Ilustración
Uso
Sintaxis
Ejemplo de código
Entornos
Información sobre licencias

Resumen

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

Ilustración

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	Base
3	Vegetación baja
4	Vegetación media
5	Vegetación alta
6	Edificio
7	Ruido bajo
8	Clave de modelo/Reservada
9	Agua
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

LASToMultipoint_3d (input, out_feature_class, average_point_spacing, {class_code}, {return}, {attribute}, {input_coordinate_system}, {file_suffix}, {z_factor}, {folder_recursion})

Parámetro	Explicación	Tipo de datos
input [input,...]	Los archivos LAS 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.	Folder or 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. El rango de valores válidos varían de 1 a 32. Por defecto, no se aplica el 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 1 2 3 4 5 6 7 8 LAST_RETURNS	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 RETURN_NUMBER NUMBER_OF_RETURNS SCAN_DIRECTION_FLAG EDGE_OF_FLIGHTLINE CLASSIFICATION SCAN_ANGLE_RANK FILE_MARKER USER_BIT_FIELD GPS_TIME COLOR_RED COLOR_GREEN COLOR_BLUE	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.	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

Ejemplo 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 (secuencia de comandos independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en una secuencia de comandos 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)