Ajouter un niveau de pyramide au MNT (3D Analyst)
Récapitulatif
Ajoute un ou plusieurs niveaux de pyramide à un jeu de données de MNT existant.
Pour en savoir plus sur la définition d'un niveau de pyramide de MNT
Utilisation
Chaque niveau de pyramide est entré comme une paire Résolution des niveaux de pyramide/Échelle de référence, délimitée par des espaces (par exemple, "20 24000" pour une taille de fenêtre de 20 et une échelle de référence de 1:24000, ou "1,5 10000" pour une tolérance Z de 1,5 et une échelle de référence de 1:10000).
La résolution des niveaux de pyramide désigne la valeur de tolérance Z ou de taille de fenêtre, utilisée pour la pyramide.
- La tolérance Z spécifie l'écart maximal possible vis-à-vis de l'altitude du MNT en résolution maximale.
- La taille de fenêtre définit la surface de tuile utilisée pour affiner les points d'altitude en sélectionnant un ou deux points de la surface selon la méthode de taille de fenêtre spécifiée lors de la création du MNT.
L'échelle de référence représente la plus grande échelle de la carte pour l'affichage d'un niveau de pyramide.
-
L'ajout d'un niveau de pyramide à un MNT l'invalide, car le MNT doit renseigner la pyramide avec les points d'altitude échantillonné à partir du niveau de pyramide précédent. Utilisez l'outil Générer un MNT après avoir ajouté un niveau de pyramide.
-
Lorsqu'il est utilisé dans une base de données ArcSDE, le MNT en entrée ne peut pas être inscrit comme versionné.
Syntaxe
| Paramètre | Explication | Type de données |
in_terrain |
The input terrain dataset. | Terrain Layer |
pyramid_level_definition ["<pyramid_level_resolution> <reference_scale>",...] |
Tolérance Z ou taille de fenêtre et échelle de référence associée pour chaque niveau de pyramide ajouté au MNT. Chaque niveau de pyramide est entré comme une paire Résolution des niveaux de pyramide/Échelle de référence, délimitée par des espaces (par exemple, "20 24000" pour une taille de fenêtre de 20 et une échelle de référence de 1:24000, ou "1,5 10000" pour une tolérance Z de 1,5 et une échelle de référence de 1:10000). La résolution des niveaux de pyramide peut être une valeur à virgule flottante, tandis que l'échelle de référence doit être entrée sous la forme d'un nombre entier. La valeur de tolérance Z représente l'écart maximal possible vis-à-vis de l'altitude du MNT en résolution maximale, tandis que la valeur de taille de fenêtre définit la surface de tuile de MNT utilisée pour affiner les points d'altitude en sélectionnant un ou deux points de la surface selon la méthode de taille de fenêtre spécifiée lors de la création du MNT. L'échelle de référence représente la plus grande échelle de carte d'application du niveau de pyramide. Lorsque le MNT s'affiche à une échelle plus grande que cette valeur, le niveau de pyramide supérieur s'affiche. | String |
pyramid_type (Facultatif) |
Type de pyramide utilisé par le jeu de données de MNT. Ce paramètre n'est pas utilisé dans ArcGIS 9.3 et versions ultérieures, car il garantit la compatibilité ascendante avec les scripts et les modèles écrits via ArcGIS 9.2. | String |
Exemple de code
The following sample demonstrates the use of this tool in the Python window:
import arcpy
from arcpy import env
arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.AddTerrainPyramidLevel_3d('test.gdb/featuredataset/terrain', 'WINDOWSIZE',
'2.5 10000; 5 25000; 10 50000')
The following sample demonstrates the use of this tool in a stand-alone Python script:
"""****************************************************************************
Name: Create Terrain from TIN
Description: This script demonstrates how to create a terrain dataset using
features extracted from a TIN. It is particularly useful in
situations where the source data used in the TIN is not available,
and the amount of data stored in the TIN proves to be too large
for the TIN. The terrain's scalability will allow improved
display performance and faster analysis. The script is designed
to work as a script tool with 5 input arguments.
****************************************************************************"""
# Import system modules
import arcpy
import exceptions, sys, traceback
from arcpy import env
# Set local variables
tin = arcpy.GetParameterAsText(0) # TIN used to create terrain
gdbLocation = arcpy.GetParameterAsText(1) # Folder that will store terran GDB
gdbName = arcpy.GetParameterAsText(2) # Name of terrain GDB
fdName = arcpy.GetParameterAsText(3) # Name of feature dataset
terrainName = arcpy.GetParameterAsText(4) # Name of terrain
try:
arcpy.CheckOutExtension("3D")
# Create the file gdb that will store the feature dataset
arcpy.management.CreateFileGDB(gdbLocation, gdbName)
gdb = '{0}/{1}'.format(gdbLocation, gdbName)
# Obtain spatial reference from TIN
SR = arcpy.Describe(tin).spatialReference
# Create the feature dataset that will store the terrain
arcpy.management.CreateFeatureDataset(gdb, fdName, SR)
fd = '{0}/{1}'.format(gdb, fdName)
# Export TIN elements to feature classes for terrain
arcpy.AddMessage("Exporting TIN footprint to define terrain boundary...")
boundary = "{0}/boundary".format(fd)
# Execute TinDomain
arcpy.ddd.TinDomain(tin, tinDomain, 'POLYGON')
arcpy.AddMessage("Exporting TIN breaklines...")
breaklines = "{0}/breaklines".format(fd)
# Execute TinLine
arcpy.ddd.TinLine(tin, breaklines, "Code")
arcpy.AddMessage("Exporting TIN nodes...")
masspoints = "{0}/masspoints".format(fd)
# Execute TinNode
arcpy.ddd.TinNode(sourceTIN, TIN_nodes)
arcpy.AddMessage("Creating terrain dataset...")
terrain = "terrain_from_tin"
# Execute CreateTerrain
arcpy.ddd.CreateTerrain(fd, terrainName, 10, 50000, "",
"WINDOWSIZE", "ZMEAN", "NONE", 1)
arcpy.AddMessage("Adding terrain pyramid levels...")
terrain = "{0}/{1}".format(fd, terrainName)
pyramids = ["20 5000", "25 10000", "35 25000", "50 50000"]
# Execute AddTerrainPyramidLevel
arcpy.ddd.AddTerrainPyramidLevel(terrain, "", pyramids)
arcpy.AddMessage("Adding features to terrain...")
inFeatures = "{0} Shape softclip 1 0 10 true false boundary_embed <None> "\
"false; {1} Shape masspoints 1 0 50 true false points_embed "\
"<None> false; {2} Shape softline 1 0 25 false false lines_embed "\
"<None> false".format(boundary, masspoints, breaklines)
# Execute AddFeatureClassToTerrain
arcpy.ddd.AddFeatureClassToTerrain(terrain, inFeatures)
arcpy.AddMessage("Building terrain...")
# Execute BuildTerrain
arcpy.ddd.BuildTerrain(terrain, "NO_UPDATE_EXTENT")
arcpy.GetMessages()
except arcpy.ExecuteError:
print arcpy.GetMessages()
except:
# Get the traceback object
tb = sys.exc_info()[2]
tbinfo = traceback.format_tb(tb)[0]
# Concatenate error information into message string
pymsg = "PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}"\
.format(tbinfo, str(sys.exc_info()[1]))
msgs = "ArcPy ERRORS:\n {0}\n".format(arcpy.GetMessages(2))
# Return python error messages for script tool or Python Window
arcpy.AddError(pymsg)
arcpy.AddError(msgs)
finally:
arcpy.CheckInExtension("3D")