Topo vers raster (Spatial Analyst)

Niveau de licence :BasicStandardAdvanced

Récapitulatif

Interpole une surface raster hydrologiquement correcte à partir de données ponctuelles, linéaires et surfaciques.

Pour en savoir plus sur la fonction Topo vers raster

Utilisation

Syntaxe

TopoToRaster (in_topo_features, {cell_size}, {extent}, {Margin}, {minimum_z_value}, {maximum_z_value}, {enforce}, {data_type}, {maximum_iterations}, {roughness_penalty}, {discrete_error_factor}, {vertical_standard_error}, {tolerance_1}, {tolerance_2}, {out_stream_features}, {out_sink_features}, {out_diagnostic_file}, {out_parameter_file}, {profile_penalty}, {out_residual_feature}, {out_stream_cliff_error_feature}, {out_contour_error_feature})
ParamètreExplicationType de données
in_topo_features
topo_input

La classe Topo précise les entités en entrée contenant les valeurs z à interpoler dans un raster de surface.

La classe Topo accepte neuf types d'entrée de données : TopoPointElevation, TopoContour, TopoStream, TopoSink, TopoBoundary, TopoLake, TopoCliff, TopoExclusion et TopoCoast.

  • TopoPointElevation ([[inFeatures,{champ}],...])

    Classe d'entités points représentant l'altitude des surfaces.

    Le champ stocke l'altitude des points.

  • TopoContour ([[inFeatures,{champ}],...])

    Classe d'entités linéaires représentant des courbes de niveau.

    Le champ stocke l'altitude des isolignes.

  • TopoStream ([inFeatures,...])

    Classe d'entités linéaires d'emplacements de cours d'eau. Tous les arcs doivent être orientés pour pointer vers l'aval. La classe d'entités doit contenir uniquement les cours d'eau constitués d'un arc unique.

  • TopoSink ([[inFeatures,{champ}],...])

    Classe d'entités points représentant les dépressions topographiques connues. La fonction Topo vers raster ne tente pas de supprimer de l'analyse les points explicitement identifiés comme des cuvettes.

    Le champ utilisé doit être un champ qui stocke l'altitude de la cuvette en question. Si l'option NONE est sélectionnée, seul l'emplacement de la cuvette est utilisé.

  • TopoBoundary ([inFeatures,...])

    Une limite est une classe d'entités contenant un seul polygone représentant la limite extérieure du raster en sortie. Les cellules du raster en sortie situées hors de cette limite ont la valeur NoData. Cette option permet de découper les zones aquatiques le long des lignes de côtes avant de réaliser le raster en sortie final.

  • TopoLake ([inFeatures,...])

    Classe d'entités surfaciques spécifiant l'emplacement des lacs. Toutes les cellules d'un raster en sortie situées à l'intérieur d'un lac ont la valeur d'altitude minimale de toutes les cellules situées le long de la ligne de berge.

  • TopoCliff ([inFeatures,...])

    Classe d'entités linéaires des falaises. Les entités linéaires de falaises doivent être orientées de telle sorte que le côté gauche de la ligne se trouve en bas de la falaise et le côté droit en haut de la falaise.

  • TopoExclusion ([inFeatures,...])

    Classe d'entités surfaciques des zones dans lesquelles les données en entrée doivent être ignorées. Ces polygones autorisent la suppression de données d'altitude du processus d'interpolation. Ce paramètre est généralement utilisé pour supprimer des données d'altitude associées à des murs de barrage et des ponts. La vallée sous-jacente peut ainsi être interpolée avec la structure de drainage connectée.

  • TopoCoast ([inFeatures,...])

    Classe d'entités surfaciques contenant le contour d'une zone côtière. Les cellules du raster en sortie final qui se trouvent à l'extérieur de ces polygones sont définies sur une valeur qui est inférieure à la limite de hauteur minimale spécifiée par l'utilisateur.

Un champ contenant les valeurs z peut être spécifié pour les types PointElevation, Contour et Sink de l'entité en entrée. Il n'y a pas d'option de champ pour les types d'entrée Boundary, Lake, Cliff, Coast, Exclusion ou Stream.

TopoInput
cell_size
(Facultatif)

Taille des cellules qui sera utilisée pour la création du raster en sortie.

Il s'agit de la valeur dans l'environnement s'il est explicitement défini, ou de la plus petite largeur ou hauteur de l'étendue des entités ponctuelles en entrée, dans la référence spatiale en entrée, divisée par 250.

Analysis Cell Size
extent
(Facultatif)

La classe Etendue détermine l'étendue du jeu de données raster en sortie.

L'interpolation se produit en dehors des limites x et y, et les cellules situées hors de cette étendue ont la valeur NoData. Pour optimiser les résultats de l'interpolation le long des tronçons du raster en sortie, les limites x et y doivent être inférieures à l'étendue des données en entrée et ce, d'au moins 10 cellules de chaque côté.

La classe Etendue a le format suivant :

  • Etendue (X min., Y min., X max., Y max.)

    où :

    • X min. - La valeur par défaut est la plus petite coordonnée x de toutes les entrées.
    • Y min. - La valeur par défaut est la plus petite coordonnée y de toutes les entrées.
    • X max. - La valeur par défaut est la plus grande coordonnée x de toutes les entrées.
    • Y max. - La valeur par défaut est la plus grande coordonnée y de toutes les entrées.

L'étendue par défaut est la plus grande de toutes les étendues parmi les données d'entité en entrée.

Extent
Margin
margin
(Facultatif)

Distance dans les cellules à interpoler au-delà de l'étendue et de la limite en sortie spécifiées.

La valeur doit être supérieure ou égale à 0 (zéro). La valeur par défaut est 20.

Si l'Etendue et la TopoBoundary du jeu de données d'entité sont identiques à la limite des données en entrée (la valeur par défaut), les valeurs interpolées le long de la limite du MNT ne correspondent pas parfaitement aux données MNT adjacentes. Elles ont en effet été interpolées en utilisant la moitié des données en tant que points à l'intérieur du raster, lesquels sont entourés de tous les côtés par des données en entrée. L'option Marge permet d'utiliser, dans l'interpolation, les données en entrée supérieures à ces limites.

Long
minimum_z_value
(Facultatif)

Valeur z minimale à utiliser dans l'interpolation.

La valeur par défaut est fixée à 20% en dessous de la valeur en entrée la plus faible.

Double
maximum_z_value
(Facultatif)

Valeur z maximale à utiliser dans l'interpolation.

La valeur par défaut est fixée à 20% au-dessus de la valeur en entrée la plus élevée.

Double
enforce
(Facultatif)

Type de drainage à appliquer.

Vous pouvez définir l'option de drainage pour tenter de supprimer toutes les cuvettes ou dépressions, afin de générer un MNT hydrologiquement correct. Si les points de cuvettes ont été explicitement identifiées en tant que données d'entité en entrée, ces dépressions ne sont pas comblées.

  • ENFORCE L'algorithme tente de supprimer toutes les cuvettes rencontrées, qu'elles soient réelles ou fictives. Il s'agit de l'option par défaut.
  • NO_ENFORCE Aucune cuvette n'est comblée.
  • ENFORCE_WITH_SINK Les points identifiés comme cuvettes dans les données d'entité en entrée représentent des dépressions topographiques connues et ne seront pas modifiés. Les cuvettes non identifiées dans les données d'entité en entrée sont considérées fictives et l'algorithme tentera de les combler.L'existence de plus de 8 000 cuvettes fictives entraîne l'échec de l'outil.
String
data_type
(Facultatif)

Type de données d'altitude dominant dans les données d'entité en entrée.

  • CONTOUR Le type dominant de données en entrée est la courbe de niveau. Il s'agit de l'option par défaut.
  • SPOT Le type dominant de données en entrée est le point.

Sélectionnez l'option adéquate pour optimiser la méthode de recherche utilisée lors de la génération des cours d'eau et des crêtes.

String
maximum_iterations
(Facultatif)

Nombre maximal d'itérations d'interpolation.

Le nombre d'itérations doit être supérieur à zéro. Une valeur par défaut de 20 convient généralement aux données d'isoligne et aux données linéaires.

Une valeur de 30 supprime moins de cuvettes. Très rarement, des valeurs supérieures (de 45 à 50) peuvent s'avérer utiles pour supprimer davantage de cuvettes ou pour définir plus de cours d'eau et de crêtes. L'itération cesse pour chaque résolution de grille lorsque le nombre maximal d'itérations est atteint.

Long
roughness_penalty
(Facultatif)

Second terme de la dérivée au carré intégré en tant que mesure de la rugosité.

La pénalité de rugosité doit être égale ou supérieure à zéro. Si le type de données en entrée principal est CONTOUR, la valeur par défaut est zéro. Si le type de données principal est SPOT, la valeur par défaut est 0,5. Les valeurs plus élevées ne sont normalement pas recommandées.

Double
discrete_error_factor
(Facultatif)

Le facteur d'erreur de discrétisation permet d'ajuster le degré de lissage lors de la conversion des données en entrée en raster.

La valeur doit être supérieure à zéro. L'ajustement standard est compris entre 0,25 et 4, et la valeur par défaut est fixée à 1. Une valeur inférieure entraîne un lissage moins important, et une valeur supérieure un lissage plus important.

Double
vertical_standard_error
(Facultatif)

Quantité d'erreurs aléatoires dans les valeurs z des données en entrée.

La valeur doit être égale ou supérieure à zéro. La valeur par défaut est zéro.

L'erreur standard verticale peut être définie sur une petite valeur positive si les données présentent des erreurs verticales aléatoires (non systématiques) significatives, avec une variance uniforme. Dans ce cas, définissez l'erreur standard verticale sur l'écart type de ces erreurs. Pour la plupart des jeux de données altimétriques, l'erreur verticale doit être fixée à zéro, mais vous pouvez lui attribuer une petite valeur positive pour stabiliser la convergence lorsque vous rastérisez des données ponctuelles avec des données linéaires de cours d'eau.

Double
tolerance_1
(Facultatif)

Cette tolérance reflète la précision et la densité des points d'altitude par rapport au drainage de surface.

Pour les jeux de données ponctuelles, définissez l'erreur standard des hauteurs comme tolérance. Pour les jeux de données isolignes, utilisez la moitié de l'intervalle moyen des isolignes.

La valeur doit être égale ou supérieure à zéro. La valeur par défaut est 2,5 si le type de données est CONTOUR et zéro si le type de données est SPOT.

Double
tolerance_2
(Facultatif)

Cette valeur de tolérance évite tout problème lié à des hauteurs d'interruptions inappropriées.

La valeur doit être supérieure à zéro. La valeur par défaut est 100 si le type de données est CONTOUR et 200 si le type de données est SPOT.

Double
out_stream_features
(Facultatif)

Classe d'entités linéaires en sortie des entités polylignes d'écoulement et des entités linéaires de crête.

Les entités linéaires sont créées au début du processus d'interpolation. Cela fournit la morphologie générale du terrain pour l'interpolation. Vous pouvez l'utiliser pour vérifier l'exactitude du drainage et de la morphologie par comparaison avec des données connues de cours d'eau et de crêtes.

Les entités polylignes sont codées comme suit :

  1. Ligne d'écoulement en entrée ne dépassant pas la falaise.
  2. Ligne d'écoulement en entrée dépassant la falaise (cascade).
  3. Application du drainage pour supprimer une cuvette fictive.
  4. Ligne d'écoulement déterminée à partir de l'angle de l'isoligne.
  5. Ligne de dorsale déterminée à partir de l'angle de l'isoligne.
  6. Code inutilisé.
  7. Data stream line side conditions.
  8. Code inutilisé.
  9. Ligne indiquant des marges élevées de données d'altitude.
Feature Class
out_sink_features
(Facultatif)

Classe d'entités points en sortie des entités ponctuelles de cuvettes résiduelles.

Il s'agit des cuvettes qui n'ont pas été spécifiées dans les données d'entité en entrée de cuvette et n'ont pas été supprimées durant le drainage. L'ajustement des valeurs de tolérance, tolerance_1 et tolerance_2, peut réduire le nombre de cuvettes résiduelles. Les cuvettes résiduelles indiquent souvent des erreurs dans les données en entrée, que l'algorithme de drainage n'a pas pu résoudre. Cela permet de détecter efficacement les erreurs d'altitude minimes.

Feature Class
out_diagnostic_file
(Facultatif)

Fichier de diagnostic en sortie répertoriant tous les paramètres et entrées utilisés, et le nombre de cuvettes supprimées à chaque résolution et itération.

File
out_parameter_file
(Facultatif)

Fichier de paramètres en sortie répertoriant l'ensemble des paramètres et entrées utilisés, exploitable avec la fonction Topo vers raster - fichier pour réexécuter l'interpolation.

File
profile_penalty
(Facultatif)

La pénalité de rugosité de la courbure longitudinale est une pénalité adaptative applicable localement qui peut être utilisée pour remplacer partiellement la courbure totale.

Elle peut donner de bons résultats avec des données d'isolignes de qualité mais peut générer une situation de convergence instable si ce n'est pas le cas. Cette entrée est définie sur 0,0 en cas d'absence de courbure longitudinale (valeur par défaut), sur 0,5 pour une courbure longitudinale modérée et sur 0,8 pour une courbure longitudinale maximale. N'utilisez pas de valeurs supérieures à 0,8.

Double
out_residual_feature
(Facultatif)

Classe d'entités ponctuelles en sortie de toutes les valeurs résiduelles d'altitude élevées telles que mises à l'échelle par l'erreur de discrétisation locale.

Vous devez examiner toutes les valeurs résiduelles supérieures à 10 à la recherche d'éventuelles erreurs dans les données de cours d'eau et d'altitude en entrée. Les valeurs résiduelles élevées mises à l'échelle indiquent des conflits entre les données d'altitude et les données linéaires de cours d'eau en entrée. Cela peut également être dû à de mauvaises drainages automatiques. Pour résoudre ces conflits, commencez par rechercher des erreurs dans les données en entrée et, le cas échéant, corrigez-les, puis ajoutez des données linéaires de cours d'eau et/ou d'altitude de point. Les valeurs résiduelles élevées non mises à l'échelle signalent généralement des erreurs d'altitude en entrée.

Feature Class
out_stream_cliff_error_feature
(Facultatif)

Classe d'entités ponctuelles en sortie des emplacements où peuvent se produire des erreurs de cours d'eau et de falaise.

Il est possible d'identifier dans la classe d'entités ponctuelles les emplacements où les cours d'eau présentent des boucles fermées, des effluents et des cours d'eau au-dessus de falaises. Les falaises dont les cellules voisines sont incohérentes avec les côtés haut et bas de la falaise sont également identifiées. Cela permet d'identifier avec précision les falaises dont la direction est incorrecte.

Les points sont codés comme suit :

  1. Circuit réel dans le réseau des données linéaires de cours d'eau des données.
  2. Circuit dans le réseau de cours d'eau tel qu'il est codé dans le raster en sortie.
  3. Circuit dans le réseau hydrographique via des lacs connectés.
  4. Points des effluents.
  5. Cours d'eau au-dessus d'une falaise (cascade).
  6. Points indiquant plusieurs débordements de cours d'eau provenant de lacs.
  7. Code inutilisé.
  8. Points à proximité de falaises dont la hauteur ne correspond pas à la direction de la falaise.
  9. Code inutilisé.
  10. Effluent circulaire supprimé.
  11. Effluent sans arrivée d'eau.
  12. Effluent tramé de la cellule en sortie différent de l'endroit où commence l'effluent .
  13. Erreur de traitement des conditions - indicateur de données de lignes d'écoulement très complexes.
Feature Class
out_contour_error_feature
(Facultatif)

Classe d'entités ponctuelles en sortie des erreurs possibles concernant les données d'isolignes en entrée.

Les isolignes présentant une hauteur incorrecte cinq fois supérieure à l'écart type des valeurs d'isolignes représentées dans le raster en sortie sont signalées dans cette classe d'entités. Les isolignes qui rejoignent d'autres isolignes dont l'altitude est différente sont signalées dans cette classe d'entités par le code 1. Cela indique clairement une erreur d'étiquette d'isoligne.

Feature Class

Valeur renvoyée

NomExplicationType de données
out_surface_raster

Raster de surface interpolé en sortie.

Raster

Exemple de code

1er exemple de la commande TopoToRaster (fenêtre Python)

Cet exemple crée un raster de surface TIFF hydrologiquement correct à partir de données ponctuelles, linéaires et surfaciques.

import arcpy
from arcpy import env  
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outTTR = TopoToRaster([TopoPointElevation([['spots', 'spot_meter']]), 
                       TopoContour([['contours', 'spot_meter']]),
                       TopoCliff(['cliff'])], 60, 
                       "#", "#", "#", "#", "NO_ENFORCE")
outTTR.save("C:/sapyexamples/output/ttrout.tif")
2è exemple de la commande TopoToRaster (script autonome)

Cet exemple crée un raster de surface Grid hydrologiquement correct à partir de données ponctuelles, linéaires et surfaciques.

# Name: TopoToRaster_Ex_02.py
# Description: Interpolates a hydrologically correct surface 
#    from point, line, and polygon data.
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "C:/sapyexamples/data"

# Set local variables
inPointElevations = TopoPointElevation([['spots.shp', 'spot_meter'], 
                                        ['spots2.shp', 'elev']])
inBoundary = TopoBoundary(['boundary.shp'])
inContours = TopoContour([['contours.shp', 'spot_meter']])
inLake = TopoLake(['lakes.shp'])
inSinks = TopoSink([['sink1.shp', 'elevation'], ['sink2.shp', 'none']])
inStream = TopoStream(['streams.shp'])
inCliff = TopoCliff(['cliff.shp'])
inCoast = TopoCoast(['coast.shp'])
inExclusion = TopoExclusion(['ignore.shp'])
                    
inFeatures = ([inPointElevations, inContours, inLake, inStream, inBoundary, inSinks, inCliff, inCoast, inExclusion])

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute TopoToRaster
outTTR = TopoToRaster(inFeatures)

# Save the output 
outTTR.save("C:/sapyexamples/output/ttrout03")

Environnements

Thèmes connexes

Informations de licence

ArcGIS for Desktop Basic: Requis ArcGIS Spatial Analyst ou 3D Analyst
ArcGIS for Desktop Standard: Requis ArcGIS Spatial Analyst ou 3D Analyst
ArcGIS for Desktop Advanced: Requis ArcGIS Spatial Analyst ou 3D Analyst
5/10/2014