ArcMap

Exposition de surface

Résumé

Crée des entités surfaciques qui représentent les mesures d'exposition découlant d'une surface du jeu de données TIN, de MNT ou LAS.

Illustration

Illustration de l'outil Exposition TIN

Utilisation

  • L'exposition représente l'orientation horizontale d'une surface et est définie en unités de degrés. Une valeur de code représentant la direction cardinale ou ordinale de sa pente est attribuée à chaque facette de la surface. Les surfaces contiguës portant le même code sont combinées en une seule entité. Le jeu de classification par défaut est défini comme suit :

    CodeDirection de pentePlage d'angle de pente

    -1

    Plat standard

    Pas de pente

    1

    Nord

    0° – 22.5°

    2

    Nord-Est

    22.5° – 67.5°

    3

    Est

    67.5° – 112.5°

    4

    Sud-Est

    112.5° – 157.5°

    5

    Sud

    157.5° – 202.5°

    6

    Sud-Ouest

    202.5° – 247.5°

    7

    Ouest

    247.5° – 292.5°

    8

    Nord-Ouest

    292.5° – 337.5°

    9

    Nord

    337.5° – 360°

  • Il est possible de spécifier des définitions de classe personnalisées via une Table de bornes de classes. La table doit comporter deux colonnes dont la première indique le point de rupture d'exposition en degrés et la seconde définit sa valeur de code. Prenons l'exemple suivant :

    InterruptionAspect_Code

    90.0

    1

    180.0

    2

    270.0

    3

    360.0

    4

    La table peut être enregistrée dans n'importe quel format pris en charge (.dbf, .txt ou une table de géodatabase). Les noms des champs n'ont pas d'importance, étant donné que le premier sera toujours utilisé pour les bornes de classe et le second pour les codes d'exposition.

Syntaxe

arcpy.ddd.SurfaceAspect(in_surface, out_feature_class, {class_breaks_table}, {aspect_field}, {pyramid_level_resolution})
ParamètreExplicationType de données
in_surface

Surface du jeu de données TIN, de MNT ou LAS à traiter.

LAS Dataset Layer; Terrain Layer; TIN Layer
out_feature_class

Classe d'entités qui va être produite par cet outil.

Feature Class
class_breaks_table
(Facultatif)

Table contenant les bornes utilisées pour définir les plages d'exposition dans la classe d'entités en sortie.

Table
aspect_field
(Facultatif)

Champ contenant les valeurs de code d'exposition.

String
pyramid_level_resolution
(Facultatif)

Résolution de taille de fenêtre ou de tolérance z du niveau de pyramide de MNT qui sera utilisée par cet outil. La valeur par défaut est 0 (ou résolution maximale).

Double

Exemple de code

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

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

import arcpy
from arcpy import env

arcpy.CheckOutExtension("3D")
env.workspace = "C:/data"
arcpy.SurfaceAspect_3d("sample.gdb/featuredataset/terrain", "sample.gdb/featuredataset/terrain_aspect")
2e exemple d'utilisation de l'outil SurfaceAspect (script autonome)

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

'''****************************************************************************
Name: SurfaceAspect Example
Description: This script demonstrates how to use the 
             SurfaceAspect and SurfaceSlope tools to generate a polygon
             that contains the intersection of both 
****************************************************************************'''

# Import system modules
import arcpy
from arcpy import env

# Obtain a license for the ArcGIS 3D Analyst extension
arcpy.CheckOutExtension("3D")

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

try:
    # List all TINs in workspace
    listTINs = arcpy.ListDatasets("","TIN")
    # Determine whether the list contains any TINs
    if len(listTINs) > 0:
        for dataset in listTINs:
            print dataset
            # Set Local Variables
            aspect = arcpy.CreateUniqueName("Aspect.shp")
            slope = arcpy.CreateUniqueName("Slope.shp")
            outFC = dataset + "_Aspect_Slope.shp"
            #Execute SurfaceAspect
            arcpy.SurfaceAspect_3d(dataset, aspect)
            #Execute SurfaceSlope
            arcpy.SurfaceSlope_3d(dataset, slope)
            #Execute SurfaceSlope
            print "Starting Intersect"
            arcpy.Intersect_analysis(aspect + " #;" + slope + " #", outFC, "ALL")
            print "Completed intersect for " + dataset
            del aspect, slope, outFC
    else:
        print "There are no TINs in the " + env.workspace + " directory."
except:
    # Returns any other error messages
    print arcpy.GetMessages(2)

del arcpy, listTINs

Environnements

Informations de licence

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

Rubriques connexes