Classer le sol LAS

Résumé

Classe des points terrestres dans des données lidar.

Utilisation

• Seuls les derniers points LAS renvoyés avec des valeurs de code de classe de 0, 1 ou 2 peuvent être reclassés comme terrestres. Si vos fichiers LAS utilisent d'autres valeurs de code de classe pour représenter des mesures non classées ou terrestres, vous pouvez utiliser l'outil Modifier les codes de classes LAS pour les réattribuer correctement. Le processus de classification ignore également les points qui sont attribués avec l'indicateur de classification de superposition ou de bruit.

• Appliquez une étendue ou une limite de traitement pour analyser les résultats de l'algorithme de classification sur une petite surface des données lidar.

• Les emplacements qui comportent des ponts et des bretelles d'autoroute doivent être examinés, car ils risquent d'être classés par erreur comme terrestres.

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

  Remarque :

  Bien que les nouvelles attributions de codes de classes qui ont été introduites avec LAS 1.4 étaient auparavant désignées comme Réservé, les codes de classes 8 et 12 sont passés de Clé de modèle et Superposition à Réservé.

Syntaxe

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

Paramètre	Explication	Type de données
in_las_dataset	Jeu de données LAS à traiter. Seuls les derniers points LAS renvoyés avec les valeurs de code de classe 0, 1 et 2 sont évalués.	LAS Dataset Layer
method	Méthode utilisée pour détecter les points terrestres. STANDARD —Cette méthode utilise une tolérance de variation de pente qui permet de capturer les ondulations graduelles de la topographie du sol généralement non détectées par l'option conservatrice. En revanche, elle ne permet pas de capturer le type de relief accidenté détecté par l'option agressive. Il s'agit de la méthode par défaut. CONSERVATIVE — Comparée aux autres options, cette méthode utilise une restriction plus sévère sur la variation de la pente du sol pour différencier le sol de la végétation de basse altitude, comme l'herbe et les arbustes. Elle convient mieux aux topographies présentant une courbure minimale. AGGRESSIVE — Détecte les zones terrestres aux reliefs accentués, tels que crêtes et sommets, qui sont généralement ignorées par la méthode standard. Cette option convient mieux la deuxième fois que vous exécutez cet outil, lorsque l'option reuse_ground est définie sur REUSE_GROUND. N'utilisez pas cette méthode dans les zones urbaines ou rurales planes, car elle risquerait de classer comme terrestres les objets plus élevés, comme des pylônes, la végétation et certains bâtiments.	String
reuse_ground (Facultatif)	Indiquez si les points terrestres actuels doivent être reclassés ou réutilisés. Désactivé : les points terrestres actuels sont reclassés. Les points qui ne sont pas considérés comme terrestres se voient réattribuer la valeur de code de classe 1 représentant des points non classés. Il s'agit de l'option par défaut. Activé : les points terrestres actuels sont acceptés sans aucun test et aident à déterminer les points non classés. Indiquez si les points terrestres actuels doivent être acceptés ou réévalués. RECLASSIFY_GROUND — Les points terrestres actuels sont reclassés. Les points qui ne sont pas considérés comme terrestres se voient réattribuer la valeur de code de classe 1 représentant des points non classés. Il s'agit de l'option par défaut. REUSE_GROUND — Les points terrestres actuels sont acceptés sans aucun test et aident à déterminer les points non classés.	Boolean
dem_resolution (Facultatif)	En spécifiant une distance, vous obtenez un sous-ensemble de points qui peuvent être classés comme terrestres, ce qui accélère le processus. Envisagez cette option si vous avez besoin d'une méthode plus rapide pour générer une surface de MNA. La distance minimale est de 0,30 mètre, mais la distance spécifiée doit être au moins égale à 1,5 fois l'espacement moyen des points des données lidar pour que ce processus puisse avoir lieu.	Linear Unit
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
extent (Facultatif)	Spécifiez l'étendue des données qui seront évaluées par cet outil.	Extent
boundary	Entité surfacique définissant la zone d'intérêt qui sera traitée par cet outil.	Feature Layer
process_entire_files (Facultatif)	Indiquez comment appliquer l'étendue de traitement. PROCESS_EXTENT —Seuls les points LAS qui intersectent la zone d'intérêt sont traités. Il s'agit de l'option par défaut. PROCESS_ENTIRE_FILES —Si une partie d'un fichier LAS intersecte la zone d'intérêt, tous les points qu'il contient, y compris ceux situés à l'extérieur de la zone d'intérêt, sont traités.	Boolean

Exemple de code

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

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

arcpy.env.workspace = 'C:/data'
arcpy.ClassifyLasGround_3d('metro.lasd', 'CONSERVATIVE', 
                           boundary='study_area.shp', 
                           process_entire_files='PROCESS_ENTIRE_FILES')

Exemple 2 d'utilisation de l'outil ClassifyLasGround (script autonome)

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

'''****************************************************************************
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())

Environnements

• Espace de travail courant
• Espace de travail temporaire
• Etendue

Informations de licence

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