Statistiques de points LAS sous forme de raster (Outils de gestion des données)
Résumé
Permet de créer un raster dont les valeurs de cellule reflètent les informations statistiques sur les points LAS.
Utilisation
Vous pouvez filtrer les points traités par cet outil à l’aide d’une combinaison de codes de classification, d’indicateurs de classification et de valeurs renvoyées dans les filtres de points de la couche de jeu de données LAS. Les filtres peuvent être définis dans la boîte de dialogue Layer Properties (Propriétés de la couche) ou à l’aide de l’outil Make LAS Dataset Layer (Créer une couche de jeu de données LAS). Les filtres de points LAS se reflètent dans les résultats obtenus pour toutes les sorties statistiques sauf la méthode Total des impulsions qui évalue les derniers points renvoyés.
En plus de filtrer des points selon certains attributs, vous pouvez filtrer les points LAS par élévation à l’aide du contrôle Elevation Range (Plage d’élévation) disponible dans le groupe Filters (Filtres) du ruban pour une couche de jeu de données LAS. Cet outil tiendra compte de la plage d’élévation définie avec ce contrôle.
L’option Point Count (Nombre de points) du paramètre Method (Méthode) évalue la densité des points du jeu de données LAS. Ces informations sont utiles en cas de détection de points chauds irréguliers présentant une forte densité de points, car ce type de données peut entraîner une dégradation des performances du traitement des fichiers LAS à l’aide d’autres outils. Si des points chauds à forte densité sont détectés, l’outil Thin LAS (LAS fin) peut être utilisé pour produire un nouveau nuage de points avec une distribution spatiale plus cohérente.
L’option Most Frequent Class Code (Code de classe le plus fréquent) du paramètre Method (Méthode) permet d’identifier le nombre d’objets uniques définis par un code de classe spécifique, en particulier si ces objets sont séparés spatialement dans le sens horizontal. Pour ce faire, appliquez d’abord un filtre de points pour isoler le code de classe spécifique. Ensuite, créez le raster de statistiques avec une taille de cellule suffisamment grande pour garantir que la connectivité des cellules représentant les points appartient au même objet et que les écarts entre les points des différents objets sont préservés. Par exemple, pour déterminer la taille de cellule adaptée à la connexion de points appartenant à un éclairage de rue, considérez l’espacement général des points et la distance entre les lampadaires adjacents. Le raster obtenu permet d’évaluer l’emprise de chaque objet. Le nombre d’objets uniques peut être obtenu en convertissant le raster en entité surfacique à l’aide de l’outil Raster To Polygon (Raster vers polygone) ou en attribuant une valeur unique à chaque agrégat de cellules à l’aide de l’outil Region Group (Groupes par régions).
Nom et emplacement du raster en sortie. Pour stocker un jeu de données raster dans une géodatabase ou dans un dossier tel qu’une grille Esri, n’ajoutez pas d’extension de fichier au nom du jeu de données raster. Une extension de fichier peut être fournie afin de définir le format du raster si vous le stockez dans un dossier (par exemple, .tif pour générer un GeoTIFF ou .img pour générer un fichier au format ERDAS IMAGINE).
Si le raster est stocké en tant que fichier .tif ou dans une géodatabase, la qualité et le type de compression raster peuvent être spécifiés à l’aide de paramètres d’environnement de géotraitement.
Raster Dataset
Méthode
(Facultatif)
Spécifie le type de statistiques à collecter pour les points LAS dans chaque cellule du raster en sortie.
Total des impulsions—Le nombre de derniers points renvoyés sera collecté.
Nombre de points—Le nombre de tous les points renvoyés sera collecté.
Valeur de dernier retour la plus fréquente—La valeur de dernier retour la plus fréquente est collectée.
Code de classe le plus fréquent—Le code de classe le plus fréquent sera collecté.
Plage de valeurs d’intensité—La plage de valeurs d’intensité sera collectée.
Plage de valeurs d’élévation—La plage des valeurs d’élévation sera collectée.
String
Type d’échantillonnage
(Facultatif)
Indique la manière dont le paramètre Sampling Value (Valeur d’échantillonnage) est interprété pour définir la taille de cellule du raster en sortie.
Observations—La valeur indiquée dans Sampling Value (Valeur d’échantillonnage) définit le nombre de colonnes ou de lignes du raster en sortie, la valeur la plus grande étant retenue. La taille de cellule est déduite en divisant le côté le plus long de l’étendue de la sortie par l’entrée du paramètre Sampling Value (Valeur d’échantillonnage). Si une valeur d’observation de 3 000 est utilisée sur un jeu de données dont le côté le plus long est de 23,67 kilomètres, la résolution du raster en sortie est de 7,89 mètres. Cette méthode offre un moyen utile de créer une sortie de taille prévisible qui peut être générée rapidement.
Taille de cellule—La taille de cellule est directement définie par le paramètre Sampling Value (Valeur d’échantillonnage). Il s’agit de l’option par défaut.
String
Valeur d’échantillonnage
(Facultatif)
Valeur utilisée en conjonction avec le paramètre Sampling Type (Type d’échantillonnage) pour définir la taille de cellule du raster en sortie.
Nom et emplacement du raster en sortie. Pour stocker un jeu de données raster dans une géodatabase ou dans un dossier tel qu’une grille Esri, n’ajoutez pas d’extension de fichier au nom du jeu de données raster. Une extension de fichier peut être fournie afin de définir le format du raster si vous le stockez dans un dossier (par exemple, .tif pour générer un GeoTIFF ou .img pour générer un fichier au format ERDAS IMAGINE).
Si le raster est stocké en tant que fichier .tif ou dans une géodatabase, la qualité et le type de compression raster peuvent être spécifiés à l’aide de paramètres d’environnement de géotraitement.
Raster Dataset
method
(Facultatif)
Spécifie le type de statistiques à collecter pour les points LAS dans chaque cellule du raster en sortie.
PULSE_COUNT—Le nombre de derniers points renvoyés sera collecté.
POINT_COUNT—Le nombre de tous les points renvoyés sera collecté.
PREDOMINANT_LAST_RETURN—La valeur de dernier retour la plus fréquente est collectée.
PREDOMINANT_CLASS—Le code de classe le plus fréquent sera collecté.
INTENSITY_RANGE—La plage de valeurs d’intensité sera collectée.
Z_RANGE—La plage des valeurs d’élévation sera collectée.
String
sampling_type
(Facultatif)
Indique la manière dont le paramètre Sampling Value (Valeur d’échantillonnage) est interprété pour définir la taille de cellule du raster en sortie.
OBSERVATIONS—La valeur indiquée dans Sampling Value (Valeur d’échantillonnage) définit le nombre de colonnes ou de lignes du raster en sortie, la valeur la plus grande étant retenue. La taille de cellule est déduite en divisant le côté le plus long de l’étendue de la sortie par l’entrée du paramètre Sampling Value (Valeur d’échantillonnage). Si une valeur d’observation de 3 000 est utilisée sur un jeu de données dont le côté le plus long est de 23,67 kilomètres, la résolution du raster en sortie est de 7,89 mètres. Cette méthode offre un moyen utile de créer une sortie de taille prévisible qui peut être générée rapidement.
CELLSIZE—La taille de cellule est directement définie par le paramètre Sampling Value (Valeur d’échantillonnage). Il s’agit de l’option par défaut.
String
sampling_value
(Facultatif)
Valeur utilisée en conjonction avec le paramètre Sampling Type (Type d’échantillonnage) pour définir la taille de cellule du raster en sortie.
Double
Exemple de code
Exemple 1 d’utilisation de la fonction LasPointStatsAsRaster (fenêtre Python)
L’exemple suivant illustre l’utilisation de cet outil dans la fenêtre Python :
Exemple 2 d’utilisation de la fonction LasPointStatsAsRaster (script autonome)
L’exemple suivant illustre l’utilisation de cet outil dans un script Python autonome :
'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for
each pulse in a given set of LAS files.
Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy
# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)
# Run CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Run LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
"PREDOMINANT_RETURNS_PER_PULSE",
"CELLSIZE", 15)