ArcMap

Statistiques de points LAS sous forme de raster

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)

  • 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ètreExplicationType 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)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
    # Execute LasPointStatsAsRaster
    arcpy.management.LasPointStatsAsRaster(lasD, statsRaster, 
                                           "PREDOMINANT_RETURNS_PER_PULSE", 
                                           "CELLSIZE", 15)
    arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
    print(arcpy.GetMessages())
except Exception as err:
    print(err.args[0])

Environnements

Informations de licence

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

Rubriques connexes