Disponible avec une licence Spatial Analyst.
Disponible avec une licence 3D Analyst.
Résumé
Interpole une surface raster hydrologiquement correcte à partir de données ponctuelles, linéaires et surfaciques.
Utilisation
Vous obtiendrez de meilleurs résultats si toutes les données en entrée sont stockées dans le même système de coordonnées planaires et ont les mêmes unités z. Des données non projetées (latitude-longitude) peuvent être utilisées ; il est possible toutefois que les résultats ne soient pas aussi précis, en particulier à des latitudes élevées.
Si plusieurs points en entrée sont compris dans une cellule en sortie,Topo vers raster utilise la valeur moyenne pour l'interpolation (seuls les 100 premiers points compris dans une cellule sont pris en compte, le reste est ignoré). Si l'algorithme rencontre trop de points, une erreur peut se produire, indiquant que le jeu de données ponctuelles contient trop de points. Le nombre maximal de points qui peuvent être utilisés équivaut à NRows * NCols, où NRows correspond au nombre de lignes dans le raster en sortie et NCols correspond au nombre de colonnes.
Si le type de données des entités en entrée est Isoligne, l'algorithme génère tout d'abord une morphologie généralisée de la surface en fonction de la courbure des isolignes. Il implémente ensuite les isolignes en tant que source d'informations sur l'altitude. Les isolignes conviennent mieux aux données à grande échelle, pour lesquelles isolignes et angles constituent des indicateurs fiables de cours d'eau et de crêtes. Elles s'avèrent aussi efficaces et moins coûteuses à plus petite échelle, pour numériser des points d'angle d'isolignes et les utiliser comme une classe d'entités points en entrée.
La taille de cellule en sortie (champ Output cell size) peut être définie par une valeur numérique ou obtenue à partir d’un jeu de données raster. Si la taille de cellule n’a pas été spécifiée de manière explicite comme étant la valeur du paramètre, elle est dérivée de l’environnement de taille de cellule si celui-ci a été spécifié. Lorsque ni le paramètre ni l’environnement de taille de cellule n’a été spécifié, mais que l’environnement Snap Raster (Raster de capture) est défini, c’est la taille de cellule du raster de capture qui est utilisée. Si aucun élément n’est spécifié, la taille de cellule est calculée d’après la largeur ou la hauteur la plus petite (selon celle qui est la plus petite des deux) de l’étendue spécifiée dans le système de coordonnées en sortie (option Output Coordinate System) de l’environnement, divisée par 250.
Si la taille de cellule est spécifiée à l’aide d’une valeur numérique, l’outil l’utilise directement pour le raster en sortie.
Si la taille de cellule est spécifiée à l’aide d’un jeu de données raster, le paramètre affiche le chemin du jeu de données raster au lieu de la valeur de la taille de cellule. La taille de cellule de ce jeu de données raster sera utilisée directement dans l’analyse, à condition que la référence spatiale du jeu de données soit identique à la référence spatiale en sortie. Si la référence spatiale du jeu de données diffère de la référence spatiale en sortie, elle sera projetée d’après la valeur du champ Cell Size Projection Method (Méthode de projection de la taille de cellule).
La représentation des cours d'eau enchevêtrés ou l'utilisation d'arcs pour représenter deux côtés d'un cours d'eau risquent de ne pas produire de résultats fiables. Les données de cours d'eau prévalent toujours sur les données de points ou d'isolignes ; c'est pourquoi les points de données d'altitude en conflit avec la pente vers chaque cours d'eau sont ignorés. Les données de cours d'eau constituent un moyen très efficace d'ajouter des informations topographiques à l'interpolation, afin de renforcer la qualité du MNT en sortie.
Pour les formats de données prenant en charge les valeurs nulles (par exemple, les classes d’entités de géodatabase fichier), une valeur nulle est ignorée lorsqu’elle est utilisée comme entrée.
Les valeurs types des paramètres Tolérance 1 et Tolérance 2 sont :
- Pour les données ponctuelles avec une échelle de 1:100 000, utilisez 5,0 et 200,0.
- Pour les données ponctuelles moins denses avec une échelle pouvant atteindre 1:500 000, utilisez 10,0 et 400,0.
- Pour les données d'isolignes avec un espacement d'isoligne de 10, utilisez 5,0 et 100,0.
La tolérance 2 dot être au moins 6 fois supérieure à la tolérance 1.
Afin de simplifier l'utilisation des entrées et des paramètres, utilisez la boîte de dialogue Topo vers raster pour créer un fichier de paramètres en sortie, qui peut être modifié dans un éditeur de texte quelconque et utilisé comme entrée pour l'outil Topo vers raster - fichier.
Cet outil étant une application qui consomme beaucoup de mémoire, il n'est pas possible de créer des rasters en sortie volumineux. Si vous devez générer une sortie de taille importante, utilisez le paramètre Marge pour produire des rasters en sortie plus petits. Pour savoir comment procéder, reportez-vous à la section Création et mosaïquage de rasters de la rubrique Fonctionnement de l'outil Topo vers raster.
Syntaxe
arcpy.3d.TopoToRaster(in_topo_features, out_surface_raster, {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ètre | Explication | Type de données |
in_topo_features topo_input | Entités en entrée contenant les valeurs z à interpoler dans un raster de surface. A chaque entité en entrée peut correspondre un champ spécifique qui contient les valeurs z et l'un des six types précisés.
Six types d'entrée sont acceptés :
Plusieurs jeux de données d'entité en entrée doivent être entourés de guillemets doubles. Les entrées individuelles sont séparées par un point-virgule, avec un espace entre chaque composant. Reportez-vous au premier exemple de code ci-dessous. | TopoInput |
out_surface_raster | Raster de surface interpolé en sortie. Il s'agit toujours d'un raster à virgule flottante. | Raster Dataset |
cell_size (Facultatif) | La taille de cellule du jeu de données raster en sortie sera créée. Ce paramètre peut être défini par une valeur numérique ou obtenu à partir d’un jeu de données raster. Si la taille de cellule n’est pas explicitement spécifiée en tant que valeur de paramètre, la valeur de la taille de cellule de l’environnement sera utilisée dans la mesure où elle est définie. Dans le cas contraire, le calcul fera appel à des règles supplémentaires pour la déterminer à partir d’autres données en entrée. Pour en savoir plus, consultez la section relative à l’utilisation. | Analysis Cell Size |
extent (Facultatif) | Etendue 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é.
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 les jeux de données d'entité {extent} et Boundary sont identiques à la limite des données en entrée (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 {margin} 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.
| String |
data_type (Facultatif) | Type de données d'altitude dominant dans les données d'entité en entrée.
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 Isoligne, la valeur par défaut est zéro. Si le type de données principal est Point, 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. La plage normale d'ajustement est comprise entre 0,25 et 4, la valeur par défaut est de 1. Une valeur plus petite diminue le lissage des données, une valeur plus grande l'augmente. | 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 Isoligne et zéro si le type de données est Point. | 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 Isoligne et 200 si le type de données est Point. | 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 |
Exemple de code
Cet exemple crée un raster de surface Grid hydrologiquement correct à partir de données ponctuelles, linéaires et surfaciques.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.TopoToRaster_3d("C:/data/spots.shp spot_meters PointElevation; C:/data/contours.shp spot_meters Contour; C:/data/cliff.shp # Cliff",
"c:/output/toporast", 32, "740825 4309275 748825 4317275",
"20","","","ENFORCE","CONTOUR","40","","1","0")Cet exemple crée un raster de surface Grid hydrologiquement correct à partir de données ponctuelles, linéaires et surfaciques.
# Name: TopoToRaster_3D_Ex_02.py
# Description: Interpolates a hydrologically correct surface
# from point, line, and polygon data.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inPointElevations = "C:/data/spots.shp spot_meters PointElevation"
inContours = "C:/data/contours.shp spot_meters Contour"
inFeatures = (inPointElevations + ";" + inContours)
outRaster = "C:/output/topoout"
# Execute TopoToRaster
arcpy.TopoToRaster_3d(inFeatures, outRaster)Environnements
Informations de licence
- Basic: Requiert 3D Analyst ou Spatial Analyst
- Standard: Requiert 3D Analyst ou Spatial Analyst
- Advanced: Requiert 3D Analyst ou Spatial Analyst
Rubriques connexes
Vous avez un commentaire à formuler concernant cette rubrique ?