Générer un MNT (3D Analyst)

Résumé

Effectue des tâches requises pour l'analyse et l'affichage d'un jeu de données de MNT.

Utilisation

  • Cet outil doit être utilisé après avoir initialement défini les entités qui participent à un MNT, avoir fait des modifications ultérieures aux entités impliquées dans le MNT ou avoir modifié la définition d'erreur systématique.

  • Il est plus efficace d'apporter toute une série de modifications et de faire ensuite une construction plutôt que d'exécuter une opération de construction après chaque modification individuelle.

Syntaxe

BuildTerrain(in_terrain, {update_extent})
ParamètreExplicationType de données
in_terrain

Jeu de données de MNT à traiter.

Terrain Layer
update_extent
(Facultatif)

Permet de recalculer l'étendue des données d'un jeu de données de MNT basé sur une taille de fenêtre lorsque la zone de données a été réduite suite à une mise à jour. Cette opération n'est pas nécessaire si l'étendue des données a augmenté ou si le jeu de données de MNT est basé sur une tolérance z. Cette option parcourt toutes les données du MNT pour déterminer la nouvelle étendue.

  • NO_UPDATE_EXTENT L'étendue du jeu de données de MNT ne sera pas recalculée. Il s’agit de l’option par défaut.
  • UPDATE_EXTENT L'étendue du jeu de données de MNT sera recalculée.
String

Sortie dérivée

NomExplicationType de données
derived_out_terrain

MNT mis à jour.

Couche de MNT

Exemple de code

1er exemple d'utilisation de l'outil BuildTerrain (fenêtre Python)

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

arcpy.env.workspace = 'C:/data'
arcpy.BuildTerrain_3d('test.gdb/featuredataset/terrain')
2e exemple d'utilisation de l'outil BuildTerrain (script autonome)

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

"""****************************************************************************
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

# 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:
    # 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 Exception as err:
    print(err)

Informations de licence

  • Basic: Requiert 3D Analyst
  • Standard: Requiert 3D Analyst
  • Advanced: Requiert 3D Analyst

Rubriques connexes