Classer le bruit LAS (3D Analyst)

Synthèse

Classe comme bruit les points LAS avec des caractéristiques spatiales anormales.

Utilisation

• Les points de bruit dans les retours LAS ont généralement un impact néfaste sur la visualisation et l’analyse des données. Par exemple, les retours d’oiseaux volant à haute altitude et des impulsions dispersées qui ont rencontré une couverture nuageuse, du brouillard, des plans d’eau et des miroirs très réfléchissants, peuvent déformer la plage z des points qui entourent cet emplacement. L’identification de ces points en tant que bruit permet ensuite de les filtrer de l’affichage et de les éliminer lorsque des données dérivées (surface d’élévation, pente ou profil de d’exposition) sont produites.

• La méthode Isolement traite les données LAS dans des groupes tuilés de 3 sur 3 en fonction de la région définie par les paramètres Largeur du voisinage et Hauteur du voisinage. Si le nombre de points LAS dans le volume d’analyse est inférieur à la valeur du paramètre Limite des points de voisinage, les points LAS sont traités en tant que bruit. La limite des points doit refléter une approximation raisonnable en fonction de la densité des points lidar et du nombre de points LAS pouvant être anticipés dans le volume d'analyse.

• Si certains retours de points LAS possèdent des valeurs z anormalement élevées ou faibles pour la région capturée par la collection lidar, utilisez la méthode Hauteur absolue pour définir le seuil de valeurs z des données et identifier rapidement les points aberrants en tant que bruit.

• Si des points LAS possèdent des valeurs z anormalement élevées ou faibles pour certaines régions, mais que ces valeurs se trouvent dans la plage des mesures valides capturées dans la collection lidar, utilisez la méthode Hauteur relative pour définir le seuil de valeurs z des données valides en fonction d'un décalage à partir du sol. Pour générer une surface au sol, filtrez le jeu de données LAS pour les points classés comme terrestres, puis utilisez l'outil Jeu de données LAS vers raster.

• Seuls les points LAS avec un code de classe de valeur 0 ou 1 sont reclassés. Si les points non classés sont représentés par une autre valeur, utilisez l'outil Modifier les codes de classes LAS pour leur attribuer la valeur 1. Lorsque les points de bruit sont en cours de classification et que la méthode Isolation (Isolement) ou Absolute Height (Hauteur absolue) est utilisée, un code de classe de valeur 7 est attribué à tous les points de bruit. Si la méthode Hauteur relative est utilisée, la valeur 7 est attribuée aux points de bruit inférieurs au seuil de paramètre Hauteur minimum, ce qui représente un niveau sonore faible. La valeur 18 est attribuée aux points de bruit supérieurs au seuil de paramètre Hauteur maximum, ce qui représente un niveau sonore élevé.

• Si vous n’êtes pas certain des paramètres à utiliser pour déterminer les points de bruit, vous pouvez exporter les points LAS détectés comme bruit sous forme d’entité ponctuelle tout en désactivant l’option de mise à jour du code de classification LAS. Si les points en sortie reflètent les résultats attendus, vous pouvez reclasser les points LAS à l’aide de ces entités via l’outil Localiser les points LAS par proximité.

• La méthode d’isolement est une opération performante qui s’exécute plus rapidement avec des groupes de plus grande taille. Indiquez la taille de groupe la plus importante possible selon la nature et la répartition des données pour obtenir de meilleures performances.

Paramètres

Boîte de dialogue

Étiquette	Explication	Type de données
Jeu de données LAS en entrée	Jeu de données LAS qui va être traité.	LAS Dataset Layer
Méthode (Facultatif)	Indique la méthode de détection du bruit qui va être utilisée. Isolement—La proximité spatiale des points LAS sera analysée dans des volumes tuilés pour déterminer les mesures de bruit, ainsi que la détection du bruit en fonction de la hauteur. Il s’agit de l’option par défaut. Hauteur relative par rapport au sol—Tous les points inférieurs à la hauteur minimum spécifiée à partir de la surface au sol et supérieurs à la hauteur maximum à partir de la surface au sol seront identifiés comme bruit. Hauteur absolue—Tous les points inférieurs à la hauteur minimum spécifiée et supérieurs à la hauteur maximum par rapport au niveau moyen de la mer seront identifiés comme bruit.	String
Modifier la classification (Facultatif)	Spécifie si les points LAS identifiés comme bruit vont être reclassés. Activé : les points de bruit vont être reclassés. Il s’agit de l’option par défaut. Désactivé : les points de bruit ne seront pas reclassés.	Boolean
Attribuer un indicateur de retenue (Facultatif)	Spécifie si l'indicateur de classification Retenu sera attribué aux points de bruit. Ce paramètre est activé lorsque le paramètre Mettre à jour la classification est sélectionné. Activé - L'indicateur de classification Retenu sera attribué aux points de bruit. Désactivé - L'indicateur de classification Retenu ne sera pas attribué aux points de bruit. Il s’agit de l’option par défaut.	Boolean
Calculer les statistiques (Facultatif)	Spécifie si des statistiques seront calculées pour les fichiers .las référencés par le jeu de données LAS. Les statistiques de calcul fournissent un index spatial pour chaque fichier .las, ce qui contribue à améliorer les performances d’analyse et d’affichage. Les statistiques améliorent également l’expérience de filtrage et de symbologie en limitant l’affichage des attributs LAS (par exemple, les codes de classification et les informations renvoyées) aux valeurs présentes dans le fichier .las. Activé : les statistiques sont calculées. Il s’agit de l’option par défaut. Désactivé : les statistiques ne sont pas calculées.	Boolean
Surface raster au sol (Facultatif)	La surface raster au sol qui sera utilisée pour calculer les valeurs de hauteur relative pour chaque point. Ce paramètre est requis lorsque le paramètre Méthode est défini sur Hauteur relative par rapport au sol.	Raster Layer
Hauteur minimale (Facultatif)	La hauteur qui définira le seuil de valeur z le moins élevé pour l’identification des points de bruit. Tout point inférieur à la valeur fournie sera classé comme bruit. Si une surface au sol est spécifiée, ce seuil repose sur un décalage par rapport au sol. Par exemple, la valeur -3 pieds signifie que tous les points situés à 3 pieds sous la surface au sol sont classés comme bruit.	Linear Unit
Hauteur maximale (Facultatif)	La hauteur qui définira le seuil de valeur z le plus élevé pour l'identification des points de bruit. Tout point supérieur à la valeur fournie sera classé comme bruit. Si une surface au sol est fournie, ce seuil repose sur un décalage par rapport au sol. Par exemple, la valeur 250 mètres signifie que tous les points situés à 250 mètres au-dessus de la surface au sol sont classés comme bruit.	Linear Unit
Limite de points du voisinage (Facultatif)	Le nombre maximum de points dans le volume d'analyse pouvant être qualifiés comme bruit à l'aide de la méthode Isolement. Si le volume d'analyse contient des points LAS égaux ou inférieurs à cette valeur, ces points sont classés comme bruit.	Long
Largeur du voisinage (Facultatif)	La taille de chaque dimension dans l’espace x,y du volume d’analyse avec de la méthode Isolement.	Linear Unit
Hauteur du voisinage (Facultatif)	La hauteur du volume d’analyse à l’aide de la méthode Isolation (Isolement).	Linear Unit
Étendue de traitement (Facultatif)	Spécifie l’étendue des données qui seront évaluées. Current Display Extent (Étendue actuellement affichée)  : l’étendue repose sur la carte ou la scène active. Cette option est uniquement disponible lorsqu’une carte est active. Dessiner l’étendue  : l’étendue est basée sur un rectangle dessiné sur la carte ou dans la scène. Cette option crée une classe d’entités dans la géodatabase du projet et ajoute une couche à la carte. La classe d’entités possède le même système de coordonnées que la carte. Remarque : Cette option n’est pas disponible dans la boîte de dialogue Environnements. Elle est disponible uniquement dans un paramètre de l’outil correspondant à un type de données d’étendue ou dans l’onglet Environnements de la boîte de dialogue de l’outil. Remarque : Lorsque l’option de mise à jour Activer et désactiver la mise à jour de l’onglet Mettre à jour est activée, vous devez activer la mise à jour dans l’onglet du ruban Mettre à jour pour dessiner l’étendue. Extent of a Layer (Étendue d’une couche)  : l’étendue repose sur une couche de carte active. Utilisez la liste déroulante pour choisir une couche disponible ou utilisez l’option Extent of data in all layers (Étendue des données dans toutes les couches) pour obtenir l’étendue combinée de toutes les couches de carte actives, y compris le fond de carte. Cette option est uniquement disponible lorsqu’une carte avec des couches est active. Chaque couche de carte comporte les options suivantes : Toutes les entités  : étendue de toutes les entités dans la couche. Entités sélectionnées : étendue de toutes les entités sélectionnées dans la couche. Entités visibles  : étendue de toutes les entités visibles dans la couche. Remarque : Les étendues des options Entités sélectionnées  et Entités visibles  sont disponibles uniquement pour les couches d’entités. Browse (Parcourir)  : l’étendue repose sur un jeu de données existant. Intersection of Inputs (Intersection des entrées)  : l’étendue repose sur l’étendue minimale ou l’étendue d’intersection de toutes les entrées. S’il n’existe aucune intersection entre les entrées, il en résulte une étendue nulle avec uniquement des zéros. Union of Inputs (Union des entrées)  : l’étendue est basée sur l’étendue maximale ou combinée de toutes les entrées. Presse-papiers  : l’étendue peut être copiée depuis et vers le Presse-papiers. Copier l’étendue  : copie les coordonnées et le système de coordonnées de l’étendue dans le Presse-papiers. Coller l’étendue  : colle les coordonnées et, éventuellement, le système de coordonnées de l’étendue depuis le Presse-papiers. Si les valeurs du Presse-papiers n’incluent pas de système de coordonnées, l’étendue utilise celui de la carte. Remarque : Les coordonnées de l’étendue sont copiées vers et collées depuis le Presse-Papiers avec la même mise en forme et le même ordre que ceux de l’objet ArcPy Extent : X min, Y min, X max, Y max et la référence spatiale. Reset Extent (Réinitialiser l’étendue)  : la valeur par défaut de l’étendue est rétablie. Manually entered coordinates (Coordonnées saisies manuellement) : les coordonnées doivent être des valeurs numériques incluses dans le système de coordonnées de la carte active. Attention : La carte peut utiliser des unités d’affichage différentes de celles des coordonnées saisies. L’utilisation d’une direction cardinale (N, S, E, O) n’est pas prise en charge. Utilisez le signe moins pour les coordonnées sud et ouest.	Extent
Traiter l'intégralité des fichiers LAS qui intersectent l'étendue (Facultatif)	Spécifie comment la zone d’intérêt est utilisée pour déterminer le traitement des fichiers .las. La zone d’intérêt est définie par la valeur du paramètre Étendue de traitement, la valeur du paramètre Limite de traitement ou les deux à la fois. Désactivé : seuls les points LAS qui intersectent la zone d'intérêt sont traités. Il s’agit de l’option par défaut. Activé : si une partie d’un fichier .las intersecte la zone d’intérêt, tous les points que ce fichier contient, y compris ceux situés à l’extérieur de la zone d’intérêt, seront traités.	Boolean
Points de bruit en sortie (Facultatif)	Les entités ponctuelles en sortie qui représentent les points LAS identifiés comme bruit.	Feature Class
Mettre à jour la pyramide (Facultatif)	Spécifie si la pyramide du jeu de données LAS est mise à jour après la modification des codes de classe. Activé : la pyramide du jeu de données LAS est mise à jour. Il s’agit de l’option par défaut. Désactivé : la pyramide du jeu de données LAS n’est pas mise à jour.	Boolean

Sortie obtenue

Étiquette	Explication	Type de données
Jeu de données LAS en sortie	Le jeu de données LAS à modifier.	LAS Dataset Layer

Python

arcpy.ddd.ClassifyLasNoise(in_las_dataset, {method}, {edit_las}, {withheld}, {compute_stats}, {ground}, {low_z}, {high_z}, {max_neighbors}, {step_width}, {step_height}, {extent}, {process_entire_files}, {out_feature_class}, {update_pyramid})

Nom	Explication	Type de données
in_las_dataset	Jeu de données LAS qui va être traité.	LAS Dataset Layer
method (Facultatif)	Indique la méthode de détection du bruit qui va être utilisée. ISOLATION—La proximité spatiale des points LAS sera analysée dans des volumes tuilés pour déterminer les mesures de bruit, ainsi que la détection du bruit en fonction de la hauteur. Il s’agit de l’option par défaut. RELATIVE_HEIGHT—Tous les points inférieurs à la hauteur minimum spécifiée à partir de la surface au sol et supérieurs à la hauteur maximum à partir de la surface au sol seront identifiés comme bruit. ABSOLUTE_HEIGHT—Tous les points inférieurs à la hauteur minimum spécifiée et supérieurs à la hauteur maximum par rapport au niveau moyen de la mer seront identifiés comme bruit.	String
edit_las (Facultatif)	Spécifie si les points LAS identifiés comme bruit vont être reclassés. CLASSIFY—Les points de bruit vont être reclassés. Il s’agit de l’option par défaut. NO_CLASSIFY—Les points de bruit ne seront pas reclassés.	Boolean
withheld (Facultatif)	Spécifie si l'indicateur de classification Retenu sera attribué aux points de bruit. Ce paramètre est activé lorsque le paramètre edit_las est défini sur CLASSER. WITHHELD—L'indicateur de classification Retenu sera attribué aux points de bruit. NO_WITHHELD—L'indicateur de classification Retenu ne sera pas attribué aux points de bruit. Il s’agit de l’option par défaut.	Boolean
compute_stats (Facultatif)	Spécifie si des statistiques seront calculées pour les fichiers .las référencés par le jeu de données LAS. Les statistiques de calcul fournissent un index spatial pour chaque fichier .las, ce qui contribue à améliorer les performances d’analyse et d’affichage. Les statistiques améliorent également l’expérience de filtrage et de symbologie en limitant l’affichage des attributs LAS (par exemple, les codes de classification et les informations renvoyées) aux valeurs présentes dans le fichier .las. COMPUTE_STATS—Les statistiques sont calculées. Il s’agit de l’option par défaut. NO_COMPUTE_STATS—Les statistiques ne sont pas calculées.	Boolean
ground (Facultatif)	La surface raster au sol qui sera utilisée pour calculer les valeurs de hauteur relative pour chaque point. Ce paramètre est requis lorsque le paramètre method est défini sur RELATIVE_HEIGHT.	Raster Layer
low_z (Facultatif)	La hauteur qui définira le seuil de valeur z le moins élevé pour l’identification des points de bruit. Tout point inférieur à la valeur fournie sera classé comme bruit. Si une surface au sol est spécifiée, ce seuil repose sur un décalage par rapport au sol. Par exemple, la valeur -3 pieds signifie que tous les points situés à 3 pieds sous la surface au sol sont classés comme bruit.	Linear Unit
high_z (Facultatif)	La hauteur qui définira le seuil de valeur z le plus élevé pour l'identification des points de bruit. Tout point supérieur à la valeur fournie sera classé comme bruit. Si une surface au sol est fournie, ce seuil repose sur un décalage par rapport au sol. Par exemple, la valeur 250 mètres signifie que tous les points situés à 250 mètres au-dessus de la surface au sol sont classés comme bruit.	Linear Unit
max_neighbors (Facultatif)	Le nombre maximum de points dans le volume d'analyse pouvant être qualifiés comme bruit à l'aide de la méthode Isolement. Si le volume d'analyse contient des points LAS égaux ou inférieurs à cette valeur, ces points sont classés comme bruit.	Long
step_width (Facultatif)	La taille de chaque dimension dans l’espace x,y du volume d’analyse avec de la méthode Isolement.	Linear Unit
step_height (Facultatif)	La hauteur du volume d’analyse à l’aide de la méthode Isolation (Isolement).	Linear Unit
extent (Facultatif)	Spécifie l’étendue des données qui seront évaluées. MAXOF : l’étendue maximale de toutes les entrées est utilisée. MINOF : l’étendue minimale commune à 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 l’étendue de traitement est appliquée. 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 que le fichier contient, y compris ceux situés à l’extérieur de la zone d’intérêt, seront traités.	Boolean
out_feature_class (Facultatif)	Les entités ponctuelles en sortie qui représentent les points LAS identifiés comme bruit.	Feature Class
update_pyramid (Facultatif)	Spécifie si la pyramide du jeu de données LAS est mise à jour après la modification des codes de classe. UPDATE_PYRAMID—La pyramide du jeu de données LAS est mise à jour. Il s’agit de l’option par défaut. NO_UPDATE_PYRAMID—La pyramide du jeu de données LAS n’est pas mise à jour.	Boolean

Sortie obtenue

Nom	Explication	Type de données
out_las_dataset	Le jeu de données LAS à modifier.	LAS Dataset Layer

Exemple de code

Exemple 1 d’utilisation de la fonction ClassifyLasNoise (fenêtre Python)

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

arcpy.env.workspace = 'C:/data'

arcpy.ddd.ClassifyLasNoise('Denver_2.lasd', "ABSOLUTE_HEIGHT", 
                           edit_las='CLASSIFY', withheld='WITHHELD', 
                           high_z='450 Feet')

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

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

'''****************************************************************************
       Name: Classify Lidar & Extract Building Footprints
Description: Extract footprint from lidar points classified as buildings, 
             regularize its geometry, and calculate the building height.

****************************************************************************'''
import arcpy

lasd = arcpy.GetParameterAsText(0)
dem = arcpy.GetParameterAsText(1)
footprint = arcpy.GetParameterAsText(2)

try:
    desc = arcpy.Describe(lasd)
    if desc.spatialReference.linearUnitName in ['Foot_US', 'Foot']:
        unit = 'Feet'
    else:
        unit = 'Meters'
    ptSpacing = desc.pointSpacing * 2.25
    sampling = '{0} {1}'.format(ptSpacing, unit)
    # Classify overlap points
    arcpy.ddd.ClassifyLASOverlap(lasd, sampling)
    # Classify ground points
    arcpy.ddd.ClassifyLasGround(lasd)
    # Filter for ground points
    arcpy.management.MakeLasDatasetLayer(lasd, 'ground', class_code=[2])
    # Generate DEM
    arcpy.conversion.LasDatasetToRaster('ground', dem, 'ELEVATION', 
                                        'BINNING NEAREST NATURAL_NEIGHBOR', 
                                        sampling_type='CELLSIZE', 
                                        sampling_value=desc.pointSpacing)
    # Classify noise points
    arcpy.ddd.ClassifyLasNoise(lasd, method='ISOLATION', edit_las='CLASSIFY', 
                               withheld='WITHHELD', ground=dem, 
                               low_z='-2 feet', high_z='300 feet', 
                               max_neighbors=ptSpacing, step_width=ptSpacing, 
                               step_height='10 feet')
    # Classify buildings
    arcpy.ddd.ClassifyLasBuilding(lasd, '7.5 feet', '80 Square Feet')
    #Classify vegetation
    arcpy.ddd.ClassifyLasByHeight(lasd, 'GROUND', [8, 20, 55], 
                                  compute_stats='COMPUTE_STATS')
    # Filter LAS dataset for building points
    lasd_layer = 'building points'
    arcpy.management.MakeLasDatasetLayer(lasd, lasd_layer, class_code=[6])
    # Export raster from lidar using only building points
    temp_raster = 'in_memory/bldg_raster'
    arcpy.management.LasPointStatsAsRaster(lasd_layer, temp_raster,
                                           'PREDOMINANT_CLASS', 'CELLSIZE', 2.5)
    # Convert building raster to polygon
    temp_footprint = 'in_memory/footprint'
    arcpy.conversion.RasterToPolygon(temp_raster, temp_footprint)
    # Regularize building footprints
    arcpy.ddd.RegularizeBuildingFootprint(temp_footprint, footprint, 
                                          method='RIGHT_ANGLES')

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

Environnements

Espace de travail courant, Système de coordonnées en sortie, Transformations géographiques