Clasificar el terreno LAS

Resumen

Clasifica puntos de suelo en datos LIDAR.

Uso

• Para la reclasificación como suelo, solo se tendrán en cuenta los últimos puntos LAS devueltos con valores de código de clase 0, 1 o 2. Si sus archivos LAS utilizan valores de código de clase diferentes para representar mediciones sin clasificar o de suelo, plantéese la posibilidad de utilizar la herramienta Cambiar códigos de clase de LAS para reasignarlos. En el proceso de clasificación también se ignorarán los puntos a los que se haya asignado el marcador de clasificación de superposición o ruido.

• Considere la posibilidad de aplicar una extensión o límite de procesamiento para revisar los resultados del algoritmo de clasificación en un área pequeña de los datos LIDAR.

• Se deben revisar las ubicaciones con puentes y vías de acceso de autovía, ya que pueden haberse clasificado erróneamente como suelo.

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

  Nota:

  Aunque el volumen de las nuevas asignaciones de código de clase que se introdujo con la versión de LAS 1.4 en principio se designó como Reservado, los códigos de clase 8 y 12 cambiaron de Clave de modelo y Superposición a Reservados.

Sintaxis

ClassifyLasGround_3d (in_las_dataset, method, {reuse_ground}, {dem_resolution}, {compute_stats}, {extent}, boundary, {process_entire_files})

Parámetro	Explicación	Tipo de datos
in_las_dataset	El dataset LAS que se va a procesar. Solo se evaluarán los últimos puntos LAS que se hayan devuelto con valores de código de clase 0, 1 y 2.	LAS Dataset Layer
method	El método usado para detectar puntos de suelo. STANDARD —Este método tiene una tolerancia para la variación de la pendiente que permite capturar ondulaciones graduales en la topografía del terreno que normalmente se omitiría con la opción conservadora, pero no se captura el tipo de relieves acusados que sí se haría mediante la opción agresiva. Este método es el predeterminado. CONSERVATIVE — Al compararlo con otras opciones, este método emplea una restricción más severa en cuanto a la variación de la pendiente del terreno, lo que le permite diferenciarlo de vegetación baja, como hierba y matorrales. Es más adecuado para topografías con curvaturas mínimas. AGGRESSIVE — Detecta áreas de terreno con relieves más acusados, como crestas y cerros que se podrían ignorar con el método estándar. Esta opción es más adecuada en una segunda pasada de esta herramienta con la opción reuse_ground establecida en REUSE_GROUND. Evite utilizar este método en zonas urbanas o en zonas rurales planas, dado que podría dar lugar a una clasificación incorrecta de objetos más altos, como torres de servicios, vegetación y partes de edificios como suelo.	String
reuse_ground (Opcional)	Especifica si se deben reclasificar o reutilizar los puntos de suelo existentes. Desactivado: se reclasificarán los puntos de suelo existentes. A los puntos que se haya determinado que no forman parte del suelo se les volverá a asignar el valor de código de clase 1, que representa puntos sin clasificar. Esta es la opción predeterminada. Activado: se aceptarán los puntos de suelo existentes, sin escrutinio, y contribuirán a la determinación de puntos sin clasificar. Especifica si se deben aceptar o reevaluar los puntos de suelo existentes. RECLASSIFY_GROUND — Se reclasificarán los puntos de suelo existentes. A los puntos que se haya determinado que no forman parte del suelo se les volverá a asignar el valor de código de clase 1, que representa puntos sin clasificar. Esta es la opción predeterminada. REUSE_GROUND — Se aceptarán los puntos de suelo existentes, sin escrutinio, y contribuirán a la determinación de puntos sin clasificar.	Boolean
dem_resolution (Opcional)	Especificar una distancia dará como resultado un único subconjunto de puntos que se evaluarán para su clasificación como suelo, agilizando así el proceso. Plantéese utilizar esta opción cuando necesite utilizar un método más rápido para generar una superficie de DEM. Aunque la distancia mínima es de 0,3 metros, para que este proceso surta efecto, la distancia especificada debe ser al menos 1,5 veces el espaciado de punto promedio de los datos LIDAR.	Linear Unit
compute_stats (Opcional)	Especifica si las estadísticas se deben calcular para los archivos LAS a los que hace referencia el dataset LAS. La existencia de las estadísticas permite las opciones de simbología y filtrado de la capa del dataset LAS para mostrar solo los valores de atributos LAS que existen en los archivos LAS. COMPUTE_STATS —Se calcularán las estadísticas. NO_COMPUTE_STATS —No se calcularán las estadísticas. Esta es la opción predeterminada.	Boolean
extent (Opcional)	Especifique la extensión de los datos que se evaluará con esta herramienta.	Extent
boundary	Entidad poligonal que define el área de interés que se va a procesar con esta herramienta.	Feature Layer
process_entire_files (Opcional)	Especifique cómo se va a aplicar la extensión de procesamiento. PROCESS_EXTENT —Solo se procesarán los puntos LAS que intersequen el área de interés. Esta es la opción predeterminada. PROCESS_ENTIRE_FILES —Si alguna parte del archivo LAS interseca el área de interés, se procesarán todos los puntos de ese archivo LAS, incluidos los que queden fuera del área de interés.	Boolean

Muestra de código

Ejemplo 1 de ClassifyLasGround (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.ClassifyLasGround_3d('metro.lasd', 'CONSERVATIVE', 
                           boundary='study_area.shp', 
                           process_entire_files='PROCESS_ENTIRE_FILES')

Ejemplo 2 de ClassifyLasGround (secuencia de comandos independiente)

En el siguiente ejemplo se muestra cómo usar esta herramienta en una secuencia de comandos independiente de Python.

'''****************************************************************************
Name:        Classify Ground & Vegetation in Forest Environment
Description: Classify points representing vegetation with LAS class code values
             of 3, 4, and 5. The code is designed for use as a script tool.
****************************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
lasd = arcpy.GetParameterAsText(2)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasd, folder_recursion=recursion)
    # Make an initial pass of ground classifier
    arcpy.ddd.ClassifyLasGround(lasd, method="Conservative")
    # Make a secondary pass to capture ridges
    arcpy.ddd.ClassifyLasGround(lasd, method="Aggressive", 
                                reuse_ground="REUSE_GROUND")
    # Classify vegetation
    arcpy.ddd.ClassifyLasByHeight(lasd, ground_source='GROUND', 
                                  height_classification=[[3, 5], 
                                                         [4, 17], 
                                                         [5, 120]], 
                                  noise='HIGH_NOISE', compute_stats="COMPUTE_STATS")
    arcpy.CheckInExtension('3D')

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

Entornos

• Espacio de trabajo actual
• Espacio de trabajo temporal
• Extensión

Información sobre licencias

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