ArcMap

Définir des codes de classes LAS à l’aide d’entités

Résumé

Classe les points LAS qui intersectent les entités en entrée.

Utilisation

  • Les points, lignes et polygones sont des entrées valides. La distance de la zone tampon est particulièrement importante pour les points et les lignes, car il est moins probable qu'un point LAS donné se situe exactement sur la ligne ou qu'il ait la même coordonnée que le point en entrée.

  • Utilisez cet outil si des polygones délimitent des plans d'eau et que vous souhaitez attribuer à tous les points LAS situés au sein du polygone la valeur de code de classification 9, qui représente l'eau.

  • Le format LAS prend en charge la classification de chaque point en fonction des spécifications définies par l'ASPRS (American Society for Photogrammetry and Remote Sensing). La plateforme ArcGIS applique le schéma de classification spécifié pour la version 1.4 du fichier LAS :

    Valeur de classification Type de classification

    0

    Jamais classé

    1

    Non attribuée

    2

    Sol

    3

    Végétation basse

    4

    Moyenne végétation

    5

    Haute végétation

    6

    Bâtiment

    7

    Niveau sonore faible

    8

    Clé de modèle/Réservé

    9

    Eau

    10

    Ferroviaire

    11

    Surface routière

    12

    Superposition/Réservé

    13

    Fil métallique

    14

    Conducteur métallique

    15

    Tour de transmission

    16

    Connecteur métallique

    17

    Tablier de pont

    18

    Niveau sonore élevé

    19 – 63

    Réservé pour la définition ASPRS (prise en charge de LAS 1.1 à 1.3 jusqu'au code de classe 31)

    32 – 255

    Défini par l'utilisateur (uniquement pris en charge dans LAS 1.0 et certaines versions de 1.4)

Syntaxe

SetLasClassCodesUsingFeatures_3d (in_las_dataset, feature_class, {compute_stats})
ParamètreExplicationType de données
in_las_dataset

Jeu de données LAS à traiter.

LAS Dataset Layer
feature_class
[[features, buffer_distance, new_class, synthetic, key_point, withheld, overlap],...]

Entrez chaque entité et ses options associées qui serviront à définir l'opération de classification sous forme d'une liste de listes (par exemple [['feature1', 6, 9, 'NO CHANGE', 'SET', 'CLEAR', 'NO CHANGE'], ['feature 2', 0, 6, 'NO CHANGE', 'NO CHANGE', 'NO CHANGE', 'NO CHANGE']]. Chaque entité possède les options suivantes :

  • features : entités servant à la classification des points LAS.
  • buffer_distance : la distance de la zone tampon s'applique aux entités en entrée avant de déterminer les points LAS qui intersectent la zone bufférisée.
  • new_class : code de classe à attribuer.
  • synthetic : l'indicateur de classification Synthétique permet d'identifier les points qui n'ont pas été obtenus par un capteur lidar.
  • key_point : l'indicateur de classification Point clé de modèle représente un sous-ensemble de points pouvant servir à capturer un niveau de détail souhaité d'un code de classe donné.
  • withheld : l'indicateur de classification Retenu désigne les données erronées qui doivent être exclues de l'analyse et de la visualisation.
  • overlap : la désignation Superposition identifie les points des analyses de superposition et n'est prise en charge que dans les fichiers LAS 1.4.
Value Table
compute_stats
(Facultatif)

Spécifie si des statistiques doivent être calculées pour les fichiers LAS référencés par le jeu de données LAS. La présence de statistiques permet au filtrage de la couche du jeu de données LAS et aux options de symbologie de n'afficher que les valeurs d'attributs LAS qui existent dans les fichiers LAS.

  • COMPUTE_STATS —Les statistiques sont calculées.
  • NO_COMPUTE_STATS —Les statistiques ne sont pas calculées. Il s’agit de l’option par défaut.
Boolean

Sortie dérivée

NomExplicationType de données
derived_las_dataset

Exemple de code

Exemple 1 d'utilisation de l'outil SetLasClassCodesUsingFeatures (fenêtre Python)

L'exemple suivant illustre l'utilisation de cet outil dans la fenêtre Python.

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.SetLasClassCodesUsingFeatures_3d("test.lasd", [["lake.shp 0 9"],
                                       ["outliers.shp", 5, "NO_CHANGE", 
                                        "NO_CHANGE", "NO_CHANGE", "SET"]], 
                                       "COMPUTE_STATS")
Exemple 2 d'utilisation de l'outil SetLasClassCodesUsingFeatures (script autonome)

L'exemple suivant illustre l'utilisation de cet outil dans un script Python autonome.

'''**********************************************************************
Name: Assign Withheld Classification Flag to Outlier Points in LAS Files
Description: Uses Locate Outliers to identify points in LAS files that
             should be assigned the 'withheld' classification flag.
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy
import exceptions, sys, traceback

try:
    arcpy.CheckOutExtension('3D')
    # Set Local Variables
    lasD = arcpy.GetParameterAsText(0)
    outliers = 'in_memory/outliers'
    # Execute LocateOutliers
    arcpy.ddd.LocateOutliers(lasD, outliers, 'APPLY_HARD_LIMIT', -10, 
                             350, 'APPLY_COMPARISON_FILTER', 1.2, 120, 
                             0.8, 8000)
    # Execute SetLasClassCodeUsingFeatures
    arcpy.ddd.SetLasClassCodesUsingFeatures(lasd, [["outliers.shp", 5, 
                                                    "NO_CHANGE", "NO_CHANGE",
                                                    "NO_CHANGE", "SET"]])
    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)

Environnements

Informations de licence

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

Rubriques connexes