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

Résumé

Crée un raster dont les valeurs de cellule reflètent les statistiques sur les mesures de fichiers LAS référencés par un jeu de données LAS.

Utilisation

• La couche de jeu de données LAS peut également servir à contrôler la mise en application des entités de contraintes de surface qui peuvent être référencées par le jeu de données LAS. Les contraintes sont appliquées lors de l’affichage ou du traitement du jeu de données LAS en tant que surface triangulée.

• 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é.

• Utilisez PREDOMINANT_LAST_RETURN avec le paramètre Méthode pour identifier des emplacements aux valeurs renvoyées supérieures qui pourraient indiquer la présence de végétation.

• Utilisez Z_RANGE avec le paramètre Méthode pour déterminer les emplacements contenant d'éventuels points aberrants.

Syntaxe

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})

Paramètre	Explication	Type de données
in_las_dataset	Jeu de données LAS à traiter.	LAS Dataset Layer
out_raster	Nom et emplacement du raster en sortie. Pour enregistrer un jeu de données raster dans une géodatabase ou dans un dossier, tel qu’une grille Esri, vous ne devez ajouter aucune 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 TIFF 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)	Type de statistiques collectées à propos des points LAS de chaque cellule du raster en sortie. PULSE_COUNT —Nombre de derniers points renvoyés. POINT_COUNT —Nombre de points provenant de toutes les valeurs renvoyées. PREDOMINANT_LAST_RETURN —Valeur retournée la plus fréquente. PREDOMINANT_CLASS —Code de classe le plus fréquent. INTENSITY_RANGE —Plage de valeurs d'intensité. Z_RANGE —Plage de valeurs d'altitude.	String
sampling_type (Facultatif)	Spécifie la méthode utilisée pour interpréter le paramètre Valeur d'échantillonnage afin de définir la résolution du raster en sortie. OBSERVATIONS —Définit le nombre de cellules qui divisent le côté le plus long de l'étendue du jeu de données LAS. CELLSIZE —Définit la taille des cellules du raster en sortie. Il s’agit de l’option par défaut.	String
sampling_value (Facultatif)	Spécifie la valeur utilisée conjointement avec le paramètre Type d'échantillonnage afin de définir la résolution 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)

2e exemple 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)

Environnements

Espace de travail courant, Etendue, Système de coordonnées en sortie, Transformations géographiques, Taille de cellule, Pyramide, Statistiques raster, Mot-clé de configuration en sortie, Validation automatique

Informations de licence

• Basic: Requiert 3D Analyst ou Spatial Analyst
• Standard: Oui
• Advanced: Oui