Statistiques de points LAS sous forme de raster (Gestion des données)

Synthèse

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 Propriétés de la couche ou à l’aide de l’outil 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.

  • L’option Nombre de points du paramètre 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, car ce type de données peut produire des performances dégradées lors du traitement des fichiers LAS par d’autres outils. Si une grande densité de points chauds est détectée, l’outil LAS fin permet de produire un nouveau nuage de points avec une distribution spatiale plus significative.

  • L’option Code de classe le plus fréquent du paramètre 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. Créez ensuite le raster de statistiques à l’aide d’une taille de cellule suffisamment grande pour garantir que la connectivité des cellules représentant les points appartiennent au même objet et que les écarts entre les points des différents objets soient 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 une entité surfacique grâce à l’outil Raster vers polygone ou en attribuant une valeur unique à chaque agrégat de cellules à l’aide de l’outil Groupes par régions.

Paramètres

ÉtiquetteExplicationType de données
Jeu de données LAS en entrée

Jeu de données LAS qui va être traité.

LAS Dataset Layer
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 que Esri Grid, 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 lorsque 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é dans un 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 impulsionsLe nombre de derniers points renvoyés sera collecté.
  • Nombre de pointsLe nombre de tous les points renvoyés sera collecté.
  • Valeur retournée la plus fréquenteLa valeur renvoyée la plus fréquente sera collectée.
  • Code de classe le plus fréquentLe 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évationLa plage des valeurs d’élévation sera collectée.
String
Type d’échantillonnage
(Facultatif)

Détermine le modèle d’interprétation du paramètre Valeur d’échantillonnage pour définir la taille de cellule du raster en sortie.

  • ObservationsLe paramètre Valeur d’échantillonnage définit le nombre de colonnes ou de lignes dans le raster en sortie en fonction du plus long. La taille de cellule sera déduite en divisant le côté le plus long de l’étendue de la sortie par l’entrée du paramètre Valeur d’échantillonnage. Si la valeur d’observation 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 sera de 7,89 mètres. Cette méthode est utile pour générer rapidement une sortie avec une taille prévisible.
  • Taille de celluleLa taille de cellule sera directement définie par le paramètre Valeur d’échantillonnage. Il s’agit de l’option par défaut.
String
Valeur d’échantillonnage
(Facultatif)

Valeur utilisée conjointement avec le paramètre Type d’échantillonnage afin de définir la taille de cellule du raster en sortie.

Double

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
NomExplicationType de données
in_las_dataset

Jeu de données LAS qui va être traité.

LAS Dataset Layer
out_raster

Nom et emplacement du raster en sortie. Pour stocker un jeu de données raster dans une géodatabase ou dans un dossier tel que Esri Grid, 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 lorsque 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é dans un 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_COUNTLe nombre de derniers points renvoyés sera collecté.
  • POINT_COUNTLe nombre de tous les points renvoyés sera collecté.
  • PREDOMINANT_LAST_RETURNLa valeur renvoyée la plus fréquente sera collectée.
  • PREDOMINANT_CLASSLe code de classe le plus fréquent sera collecté.
  • INTENSITY_RANGELa plage de valeurs d’intensité sera collectée.
  • Z_RANGELa plage des valeurs d’élévation sera collectée.
String
sampling_type
(Facultatif)

Détermine le modèle d’interprétation du paramètre Valeur d’échantillonnage pour définir la taille de cellule du raster en sortie.

  • OBSERVATIONSLe paramètre Valeur d’échantillonnage définit le nombre de colonnes ou de lignes dans le raster en sortie en fonction du plus long. La taille de cellule sera déduite en divisant le côté le plus long de l’étendue de la sortie par l’entrée du paramètre Valeur d’échantillonnage. Si la valeur d’observation 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 sera de 7,89 mètres. Cette méthode est utile pour générer rapidement une sortie avec une taille prévisible.
  • CELLSIZELa taille de cellule sera directement définie par le paramètre Valeur d’échantillonnage. Il s’agit de l’option par défaut.
String
sampling_value
(Facultatif)

Valeur utilisée conjointement avec le paramètre Type d’échantillonnage afin de définir la taille de cellule du raster en sortie.

Double

Exemple de code

Exemple 1 d’utilisation de l’outil LasPointStatsAsRaster (fenêtre Python)

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

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
Exemple 2 d’utilisation de l’outil 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)

# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
                                       "PREDOMINANT_RETURNS_PER_PULSE",
                                       "CELLSIZE", 15)

Informations de licence

  • Basic: Nécessite 3D Analyst ou Spatial Analyst
  • Standard: Oui
  • Advanced: Oui

Rubriques connexes