ArcMap

Экспозиция поверхности

Сводка

Создает полигональные объекты, представляющие измерения экспозиции, полученные из поверхности TIN, terrain или набора данных LAS.

Иллюстрация

Иллюстрация инструмента Экспозиция TIN

Использование

  • Экспозиция представляет собой горизонтальную ориентацию поверхности и определяется в единицах градусов. Каждой грани поверхности присваивается значение, в котором закодировано направление наклона поверхности по сторонам света; соприкасающиеся области с одинаковым кодовым значением объединяются в один объект. По умолчанию определена следующая схема классификации:

    КодНаправление уклонаДиапазон угла экспозиции

    -1

    Плоский

    Нет уклона

    1

    Север

    0° – 22.5°

    2

    Северо-восток

    22.5° – 67.5°

    3

    Восток

    67.5° – 112.5°

    4

    Юго-восток

    112.5° – 157.5°

    5

    Юг

    157.5° – 202.5°

    6

    Юго-запад

    202.5° – 247.5°

    7

    Запад

    247.5° – 292.5°

    8

    Северо-запад

    292.5° – 337.5°

    9

    Север

    337.5° – 360°

  • Возможно применение пользовательской классификации, определенной при помощи Таблицы разбиения на классы. Таблица должна иметь два столбца, первый из которых указывает границу класса экспозиции в градусах, а второй – соответствующее значение кода. Рассмотрите следующий пример:

    ГраницаКод экспозиции

    90.0

    1

    180.0

    2

    270.0

    3

    360.0

    4

    Таблица может быть в любом поддерживаемом формате (.dbf, .txt или быть таблицей базы геоданных). Имена полей не имеют принципиального значения, поскольку первое из них используется для разбиения на классы, а второе – для кодов экспозиции.

Синтаксис

arcpy.ddd.SurfaceAspect(in_surface, out_feature_class, {class_breaks_table}, {aspect_field}, {pyramid_level_resolution})
ПараметрОбъяснениеТип данных
in_surface

Поверхность TIN, набора данных terrain или LAS, которые будут обрабатываться.

LAS Dataset Layer; Terrain Layer; TIN Layer
out_feature_class

Класс объектов, который будет получен в результате работы инструмента.

Feature Class
class_breaks_table
(Дополнительный)

Таблица, содержащая границы классов, которые будут использоваться для определения диапазонов экспозиции в выходном классе пространственных объектов.

Table
aspect_field
(Дополнительный)

Поле, содержащее значения кодов экспозиции.

String
pyramid_level_resolution
(Дополнительный)

z-допуск или разрешение размера окна для уровня пирамидных слоев Terrain, которое будет использоваться этим инструментом. Значением по умолчанию является 0, полное разрешение.

Double

Пример кода

SurfaceAspect – пример 1 (окно Python)

В следующем примере показано использование этого инструмента в окне 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")
SurfaceAspect, пример 2 (автономный скрипт)

В следующем примере показано использование этого инструмента в автономном скрипте Python.

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

Параметры среды

Информация о лицензиях

  • Basic: Требуется 3D Analyst
  • Standard: Требуется 3D Analyst
  • Advanced: Требуется 3D Analyst

Связанные разделы