Classer LAS par hauteur (3D Analyst)

Résumé

Reclasse les points Lidar en fonction de leur hauteur par rapport à la surface du sol.

Utilisation

  • Cet outil permet de reclasser les points LAS dont les valeurs de codes de classes sont égales à 0 ou 1 en fonction de leur hauteur par rapport à la surface du sol créée à l'aide de retours lidar auxquels les valeurs de codes de classes 2 ou 8 sont attribuées. Si les retours terrestres lidar utilisent des valeurs différentes de 2 ou 8, utilisez l'outil Modifier les codes de classes LAS pour que la définition du code de classe soit conforme aux spécifications de format de fichier LAS.

  • Vous pouvez utiliser cet outil pour classer la végétation dans les données lidar collectées sur des zones éloignées avec une présence de bâtiments minimale.

  • Appliquez une étendue de traitement pour juger de la pertinence des valeurs de plage Z désignées sur un sous-ensemble de points LAS. Lorsque vous êtes satisfait de la région reclassée, vous pouvez l'appliquer à une plus grande étendue.

  • 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

    Construction

    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

arcpy.3d.ClassifyLasByHeight(in_las_dataset, ground_source, height_classification, {noise}, {compute_stats}, {extent}, {process_entire_files}, boundary)
ParamètreExplicationType de données
in_las_dataset

Jeu de données LAS qui va être traité. Seuls les points LAS avec les valeurs de code de classe 0 et 1 sont évalués.

LAS Dataset Layer
ground_source

Source des mesures terrestres qui seront utilisées pour déterminer la hauteur au-dessus du sol.

  • GROUNDLes points LAS désignés avec le code de classification terrestre de valeur 2 et le code de clé de modèle de valeur 8 seront utilisés.
  • MODEL_KEYSeuls les points LAS désignés avec le code de classification de clé de modèle de valeur 8 seront utilisés.
String
height_classification
[[class_code, height_from_ground],...]

La valeur de code de classe qui sera attribuée aux points LAS compris dans la plage de valeurs découlant de la hauteur spécifiée par rapport au sol. L'ordre d'entrée influence les plages de hauteur qui sont utilisées pour définir la reclassification des points LAS. La plage Z de la première entrée s'étendra de la surface du sol à la valeur height_from_ground spécifiée. La plage Z des entrées suivantes s'étendra de la limite supérieure de l'entrée précédente jusqu'à sa propre valeur height_from_ground.

Value Table
noise
(Facultatif)

Indique si les points doivent être reclassés comme sonores en fonction de leur proximité avec le sol. Les artefacts sonores peuvent être introduits dans les données Lidar par des erreurs de capteur et l'interception accidentelle d'obstructions aériennes, par exemple des oiseaux sur la trajectoire de l'impulsion Lidar.

  • ALL_NOISELes niveaux sonores faibles et élevés peuvent tous les deux être classés.
  • HIGH_NOISESeuls les points situés au-dessus de la hauteur maximum dans la table de classification LAS seront reclassés comme niveau sonore élevé.
  • LOW_NOISESeuls les points situés au-dessous de la surface au sol seront reclassés comme niveau sonore. Cela est possible lorsque tous les points au sol servent à définir la surface au sol.
  • NONEAucun point ne sera reclassé comme niveau sonore.
String
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. Le calcul des statistiques fournit un index spatial pour chaque fichier LAS, ce qui améliore les performances d’analyse et d’affichage. Les statistiques permettent également d’optimiser l’utilisation des filtres et de la symbologie en limitant l’affichage des attributs LAS (tels que les codes de classification et les informations renvoyées) aux valeurs présentes dans le fichier LAS.

  • COMPUTE_STATSLes statistiques sont calculées.
  • NO_COMPUTE_STATSLes statistiques ne sont pas calculées. Il s’agit de l’option par défaut.
Boolean
extent
(Facultatif)

Spécifie l’étendue des données qui seront évaluées par cet outil.

  • MAXOF : l’étendue maximale de toutes les entrées est utilisée.
  • MINOF : l’étendue minimale de toutes les entrées est utilisée.
  • DISPLAY : l’étendue est égale à l’affichage visible.
  • Nom de la couche : l’étendue de la couche spécifiée est utilisée.
  • Objet Extent : l’étendue de l’objet spécifié est utilisée.
  • Chaîne délimitée par des espaces de coordonnées géographiques : l’étendue de la chaîne spécifiée est utilisée. Les coordonnées sont exprimées dans l’ordre x-min, y-min, x-max, y-max.
Extent
process_entire_files
(Facultatif)

Indique comment appliquer l’étendue de traitement.

  • PROCESS_ENTIRE_FILESSeuls les points LAS qui se trouvent dans l'étendue de traitement seront évalués. Il s’agit de l’option par défaut.
  • PROCESS_EXTENTTous les points des fichiers LAS qui intersectent l'étendue de traitement sont évalués.
Boolean
boundary

Entité surfacique qui définit la région pour laquelle des points terrestres LAS sont évalués.

Feature Layer

Sortie dérivée

NomExplicationType de données
out_las_dataset

Le jeu de données LAS qui a été modifié.

Couche du jeu de données LAS

Exemple de code

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

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

arcpy.env.workspace = 'C:/data'
arcpy.ClassifyLasByHeight_3d('lidar.lasd', 'Ground', 
                             [[3, 5], [4, 17], [5, 120]], 'HIGH_NOISE')
Exemple 2 d'utilisation de l'outil ClassifyLasByHeight (script autonome)

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

'''****************************************************************************
Name:        Classify Vegetation Points
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
import exceptions, sys, traceback

# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
lasd = arcpy.GetParameterAsText(2)
extent = arcpy.GetParameter(3)
calcStats = arcpy.GetParameter(4)

try:
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasd, folder_recursion=recursion)
    # Execute ChangeLasClassCodes
    arcpy.ddd.ClassifyLasByHeight(lasd, ground_source='GROUND', 
                                  height_classification=[[3, 5], 
                                                         [4, 17], 
                                                         [5, 120]], 
                                  noise='ALL_NOISE', compute_stats=calcStats, 
                                  extent=extent)

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

Informations de licence

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

Rubriques connexes