Jeu de données LAS vers raster (Conversion)

Synthèse

Crée un raster à l’aide de valeurs d’altitude, d’intensité ou RVB stockées dans les points Lidar référencés par le jeu de données LAS.

Illustration

Jeu de données LAS vers raster

Utilisation

  • Il est possible que limiter les points LAS affichés et traités par la couche de jeu de données LAS en sélectionnant une combinaison de codes de classification, de balises de classification et de valeurs renvoyées dans les paramètres de filtrage de la couche. Les filtres peuvent être définis dans la boîte de dialogue Propriétés de la couche ou l’outil Créer une couche de jeu de données LAS.

  • 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 sous forme de surface triangulée.

  • Lors de l’exportation d’un grand raster, spécifiez la valeur du paramètre Type de données en sortie en tant qu’entier pour économiser de l’espace disque si les exigences en termes de précision des valeurs z peuvent être représentées par des données de type entier.

  • Nous vous conseillons vivement d'inclure une limite de zone d'étude sous la forme d'une contrainte de découpage dans la définition du jeu de données LAS en entrée. Tout d'abord, cela empêche toute interpolation en dehors de l'étendue des données réelles de l'étude. Ensuite, les performances risquent d'être considérablement réduites si vous utilisez des options de voisin naturel et que la zone de données n'est pas correctement définie.

    Remarque :

    Si vous utilisez BINNING, seules les contraintes de découpage, d'effacement et de suppression sont reconnues. Les lignes de fracture et les points d'ancrage ne le sont pas. L'option de triangulation reconnaît tous les types de contraintes, mais son exécution prend davantage de temps.

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
Champ de valeur
(Facultatif)

Données lidar qui seront utilisées pour générer la sortie raster.

  • ELEVATIONL'altitude des fichiers lidar sera utilisée pour créer le raster. Il s'agit de l'option par défaut.
  • INTENSITYLes informations sur l'intensité provenant des fichiers lidar seront utilisées pour créer le raster.
  • RGBLes valeurs RVB des points lidar seront utilisées pour créer une image 3 canaux.
String
Type d’interpolation
(Facultatif)

Technique d'interpolation qui servira à déterminer les valeurs de cellule du raster en sortie.

L'approche de classification fournit une méthode d'attribution de cellule qui permet de déterminer chaque cellule en sortie à l'aide des points compris dans son étendue, ainsi qu'une méthode de remplissage de vide servant à déterminer la valeur des cellules ne contenant aucun point LAS.

Méthodes d'attribution de cellule

  • AVERAGE - Attribue la valeur moyenne de tous les points de la cellule. Il s'agit de l'option par défaut.
  • MINIMUM - Attribue la valeur minimale trouvée dans les points de la cellule.
  • MAXIMUM - Attribue la valeur maximale trouvée dans les points de la cellule.
  • IDW - Utilise l'interpolation Pondération par l'inverse de la distance pour déterminer la valeur de la cellule.
  • NEAREST - Utilise le rapport voisin le plus proche pour déterminer la valeur de la cellule.

Méthodes de remplissage de vide

  • NONE - NoData est attribué à la cellule.
  • SIMPLE - Calcule les valeurs des cellules de données se trouvant dans l'entourage immédiat d'une cellule NoData pour supprimer les petits vides.
  • LINEAR - Effectue une triangulation sur les surfaces vides et utilise l'interpolation linéaire sur la valeur triangulée pour déterminer la valeur de la cellule. Il s'agit de l'option par défaut.
  • NATURAL_NEIGHBOR - Utilise l'interpolation par voisins naturels pour déterminer la valeur de la cellule.

Les méthodes d'interpolation de triangulation calculent les valeurs de cellule au moyen d'une approche TIN et permettent d'accélérer le traitement en affinant l'échantillonnage des données LAS à l'aide de la technique Taille de fenêtre.

Méthodes de triangulation

  • Linéaire - Utilise l'interpolation linéaire pour déterminer les valeur de cellule.
  • Voisins naturels - Utilise l'interpolation par voisins naturels pour déterminer la valeur de la cellule.

Méthodes de sélection de la taille de fenêtre

  • Maximum - Le point à la valeur la plus élevée dans chaque taille de fenêtre est conservé. Il s'agit de l'option par défaut.
  • Minimum - Le point à la valeur la plus faible dans chaque taille de fenêtre est conservé.
  • La plus proche de la moyenne - Le point dont la valeur est la plus proche de la moyenne de toutes les valeurs de points dans la taille de la fenêtre est conservé.
Interpolate
Type de données en sortie
(Facultatif)

Détermine le type de valeurs numériques qui seront stockées dans le raster en sortie.

  • Virgule flottanteLe raster en sortie utilisera la virgule flottante 32 bits, qui prend en charge les valeurs comprises entre -3,402823466e+38 et 3,402823466e+38. Il s’agit de l’option par défaut.
  • EntierLe raster en sortie utilisera une profondeur de bit d’entiers appropriée. Cette option arrondit les valeurs z au nombre entier le plus proche et écrit un entier dans chaque valeur de cellule raster.
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
Facteur Z
(Facultatif)

Facteur par lequel les valeurs z sont multipliées. Cette option permet généralement de convertir les unités linéaires z afin qu’elles correspondent aux unités linéaires x,y. La valeur par défaut est 1, qui laisse les valeurs d'altitude inchangées. Ce paramètre n’est pas disponible si la référence spatiale de la surface en entrée possède un datum z dont l’unité linéaire est spécifiée.

Double

arcpy.conversion.LasDatasetToRaster(in_las_dataset, out_raster, {value_field}, {interpolation_type}, {data_type}, {sampling_type}, {sampling_value}, {z_factor})
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
value_field
(Facultatif)

Données lidar qui seront utilisées pour générer la sortie raster.

  • ELEVATIONL'altitude des fichiers lidar sera utilisée pour créer le raster. Il s'agit de l'option par défaut.
  • INTENSITYLes informations sur l'intensité provenant des fichiers lidar seront utilisées pour créer le raster.
  • RGBLes valeurs RVB des points lidar seront utilisées pour créer une image 3 canaux.
String
interpolation_type
"BINNING {cell_assignment_type} {void_fill_method}" or "TRIANGULATION {interpolation_method} {point_thinning_type} {point_selection_method} {resolution}"
(Facultatif)

Technique d'interpolation qui servira à déterminer les valeurs de cellule du raster en sortie.

L'approche de classification fournit une méthode d'attribution de cellule qui permet de déterminer chaque cellule en sortie à l'aide des points compris dans son étendue, ainsi qu'une méthode de remplissage de vide servant à déterminer la valeur des cellules ne contenant aucun point LAS.

Méthodes d'attribution de cellule

  • AVERAGE - Attribue la valeur moyenne de tous les points de la cellule. Il s'agit de l'option par défaut.
  • MINIMUM - Attribue la valeur minimale trouvée dans les points de la cellule.
  • MAXIMUM - Attribue la valeur maximale trouvée dans les points de la cellule.
  • IDW - Utilise l'interpolation Pondération par l'inverse de la distance pour déterminer la valeur de la cellule.
  • NEAREST - Utilise le rapport voisin le plus proche pour déterminer la valeur de la cellule.

Méthodes de remplissage de vide

  • NONE - NoData est attribué à la cellule.
  • SIMPLE - Calcule les valeurs des cellules de données se trouvant dans l'entourage immédiat d'une cellule NoData pour supprimer les petits vides.
  • LINEAR - Effectue une triangulation sur les surfaces vides et utilise l'interpolation linéaire sur la valeur triangulée pour déterminer la valeur de la cellule. Il s'agit de l'option par défaut.
  • NATURAL_NEIGHBOR - Utilise l'interpolation par voisins naturels pour déterminer la valeur de la cellule.

Les méthodes d'interpolation de triangulation calculent les valeurs de cellule au moyen d'une approche TIN et permettent d'accélérer le traitement en affinant l'échantillonnage des données LAS à l'aide de la technique Taille de fenêtre.

Méthodes de triangulation

  • Linéaire - Utilise l'interpolation linéaire pour déterminer les valeur de cellule.
  • Voisins naturels - Utilise l'interpolation par voisins naturels pour déterminer la valeur de la cellule.

Méthodes de sélection de la taille de fenêtre

  • Maximum - Le point à la valeur la plus élevée dans chaque taille de fenêtre est conservé. Il s'agit de l'option par défaut.
  • Minimum - Le point à la valeur la plus faible dans chaque taille de fenêtre est conservé.
  • La plus proche de la moyenne - Le point dont la valeur est la plus proche de la moyenne de toutes les valeurs de points dans la taille de la fenêtre est conservé.
Interpolate
data_type
(Facultatif)

Détermine le type de valeurs numériques qui seront stockées dans le raster en sortie.

  • FLOATLe raster en sortie utilisera la virgule flottante 32 bits, qui prend en charge les valeurs comprises entre -3,402823466e+38 et 3,402823466e+38. Il s’agit de l’option par défaut.
  • INTLe raster en sortie utilisera une profondeur de bit d’entiers appropriée. Cette option arrondit les valeurs z au nombre entier le plus proche et écrit un entier dans chaque valeur de cellule raster.
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
z_factor
(Facultatif)

Facteur par lequel les valeurs z sont multipliées. Cette option permet généralement de convertir les unités linéaires z afin qu’elles correspondent aux unités linéaires x,y. La valeur par défaut est 1, qui laisse les valeurs d'altitude inchangées. Ce paramètre n’est pas disponible si la référence spatiale de la surface en entrée possède un datum z dont l’unité linéaire est spécifiée.

Double

Exemple de code

Exemple 1 d'utilisation de l'outil LasDatasetToRaster (fenêtre Python)

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

arcpy.env.workspace = 'C:/data'
arcpy.ddd.LasDatasetToRaster('baltimore.lasd', 'baltimore.tif', 'INTENSITY',
                          'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                          'CELLSIZE', 10, 3.28)
Exemple 2 d'utilisation de l'outil LasDatasetToRaster (script autonome)

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

'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
             ground measurements from LAS files to a raster using a
             LAS dataset. This sample is designed to be used as a script
             tool.
*********************************************************************'''
# Import system modules
import arcpy

try:
    # Set Local Variables
    inLas = arcpy.GetParameterAsText(0)
    recursion = arcpy.GetParameterAsText(1)
    surfCons = arcpy.GetParameterAsText(2)
    classCode = arcpy.GetParameterAsText(3)
    returnValue = arcpy.GetParameterAsText(4)
    spatialRef = arcpy.GetParameterAsText(5)
    lasD = arcpy.GetParameterAsText(6)
    outRaster = arcpy.GetParameterAsText(7)
    cellSize = arcpy.GetParameter(8)
    zFactor = arcpy.GetParameter(9)

    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
    # Execute MakeLasDatasetLayer
    lasLyr = arcpy.CreateUniqueName('Baltimore')
    arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
    # Execute LasDatasetToRaster
    arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
                              'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                              'CELLSIZE', cellSize, zFactor)
    print(arcpy.GetMessages())

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

except Exception as err:
    print(err.args[0])

finally:
    arcpy.management.Delete(lasLyr)

Rubriques connexes