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

Краткая информация
Рисунок
Использование
Синтаксис
Пример кода
Параметры среды
Информация о лицензировании

Краткая информация

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

Рисунок

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

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

Код	Направление уклона	Диапазон угла экспозиции
-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 или быть таблицей базы геоданных). Имена полей не имеют принципиального значения, поскольку первое из них используется для разбиения на классы, а второе – для кодов экспозиции.

Граница	Код экспозиции
90.0	1
180.0	2
270.0	3
360.0	4

Синтаксис

SurfaceAspect_3d (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