Кривизна (Инструменты Spatial Analyst)
Краткая информация
Вычисляет кривизну растровой поверхности, дополнительно определяя профильную и плановую кривизну.
Использование
-
Основные выходные данные этого инструмента – кривизна поверхности, вычисленная по принципу "ячейка за ячейкой" и подобранная для конкретной ячейки и восьми соседних с ней ячеек. Кривизна – это вторая производная от поверхности или уклон уклона. Могут быть дополнительно получены два типа кривизны; профильная кривизна – это направление максимального уклона, и плановая кривизна – это кривизна в направлении, перпендикулярном направлению максимального уклона.
-
Положительное значение кривизны указывает на то, что поверхность в этой ячейке выпуклая. Отрицательное значение кривизны указывает на то, что поверхность в этой ячейке вогнутая. Значение 0 указывает, что поверхность плоская.
-
В выходных данных профильной кривизны отрицательное значение указывает на то, что поверхность в этой ячейке выпуклая. Положительное значение профильной кривизны указывает на то, что поверхность в этой ячейке вогнутая. Значение 0 указывает, что поверхность плоская.
-
В выходных данных плановой кривизны положительное значение указывает на то, что поверхность в этой ячейке выпуклая. Отрицательное значение плановой кривизны указывает на то, что поверхность в этой ячейке вогнутая. Значение 0 указывает, что поверхность плоская.
-
Единицы измерения выходного растра кривизны, так же, как и единицы измерения дополнительного выходного растра профильной кривизны и дополнительного выходного растра плановой кривизны, – одна сотая z-единицы (1/100). Приемлемые значения всех трех выходных растров для холмистого рельефа, как правило, находятся в диапазоне от -0.5 до 0.5; в то время, как для территорий с крутыми, скалистыми горами значения могут варьироваться от -4 до 4. Обратите внимание, что для некоторых растровых поверхностей значения могут выходить за рамки этого диапазона.
Когда возникает необходимость пересчитать входной растр, используется техника Билинейной интерполяции (Bilinear). Пример, когда входной растр может быть пересчитан, – когда выходные система координат, экстент или размер ячеек отличаются от входных.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_raster |
Входной растр поверхности. | Raster Layer |
z_factor (дополнительно) | Количество наземных единиц измерения координат x,y в одной единице измерения z-значения поверхности. Коэффициент z (z-factor) приводит в соответствие единицы измерения z-значений в том случае, если они отличаются от единиц измерения координат x,y входной поверхности. При вычислении результирующей выходной поверхности z-значения входной поверхности умножаются на коэффициент по z. Если координаты x,y и z-значения приведены в одной и той же системе координат, коэффициент z равен 1. Это значение используется по умолчанию. Если координаты x,y и z-значения приведены в отличающихся единицах измерения, для коэффициента по z должно быть задано соответствующее значение, или же результаты будут некорректными. Например, если единицы измерения для z-значений – футы, а координаты x,y приведены в метрах, для преобразования z-значений из футов в метры вы должны использовать коэффициент по z, равный 0,3048 (1 фут = 0,3048 метра). | Double |
out_profile_curve_raster (дополнительно) |
Выходной набор растровых данных профильной кривизны. Это кривизна поверхности в направлении уклона. | Raster Dataset |
out_plan_curve_raster (дополнительно) |
Выходной набор растровых данных плановой кривизны. Это кривизна поверхности в направлении, перпендикулярном направлению уклона. | Raster Dataset |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_curvature_raster |
Выходной растр кривизны. | Raster |
Пример кода
В этом примере создается растр кривизны из входного растра поверхности, и применяется коэффициент z.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outCurve = Curvature("elevation", 1.094)
outCurve.save("C:/sapyexamples/output/outcurv01")
В этом примере создается растр кривизны из входного растра поверхности, и применяется коэффициент z.
# Name: Curvature_Ex_02.py
# Description: Calculates the curvature of a raster surface,
# optionally including profile and plan curvature.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inRaster = "elevation"
zFactor = 1.094
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute Curvature
outCurve = Curvature(inRaster, 1.094)
# Save the output
outCurve.save("C:/sapyexamples/output/outcurv02")