Распределение по стоимостному расстоянию (Spatial Analyst)
Краткая информация
Вычисляет для каждой ячейки ближайший к ней источник, основываясь на наименьшей совокупной стоимости перемещения по поверхности стоимости.
Более подробно о работе инструментов стоимостного расстояния
Рисунок
 |
Cost_Alloc = CostAllocation(Source_Ras, Cost_Ras) |
Использование
-
Входные данные источников могут быть представлены классом объектов или растром.
Когда входные данные источников – растр, набор ячеек источников состоит из всех ячеек в растре источников, которые имеют действительные значения. Ячейки, имеющие значение NoData, не включаются в набор источников. Нулевое значение рассматривается как истинный источник. Исходный растр можно легко создать с помощью инструментов извлечения.
-
Когда входные данные источников представлены классом пространственных объектов, местоположения источников внутренне конвертируются в растр до выполнения анализа. Разрешение растра можно настраивать с помощью параметра среды Размер ячейки (Cell Size). Выходное разрешение по умолчанию будет соответствовать разрешению входного растра стоимости.
При использовании данных полигонального объекта в качестве входных исходных данных следует особенно тщательно выбирать способ обработки размера выходной ячейки, если его можно охарактеризовать как грубый относительно подробных сведений во входных данных. Процесс внутренней растеризации задействует тот же метод присвоения значений ячейкам (Cell assignment type) по умолчанию, что и инструмент Полигон в растр (Polygon to Raster), который является CELL_CENTER. Это означает, что данные, не расположенные по центру ячейки, не будут включены в промежуточные растеризованные исходные выходные данные, поэтому не будут представлены в расчётах расстояния. Например, если источником является ряд небольших полигонов, например периметры знаний, которые невелики относительно размеру выходной ячейки, возможно, что только некоторые из них попадут по центру выходных растровых ячеек, и, видимо, большинство остальных объекты будут потеряны в анализе.
Во избежание этой ситуации в качестве промежуточного шага можно напрямую растеризовать входные объекты с помощью инструмента Полигон в растр (Polygon to Raster), задать Поле приоритета (Priority field) и использовать полученные выходные данные в качестве входных данных инструмента Расстояние (Distance). Кроме того, можно выбрать достаточно небольшой размер ячейки, чтобы получить достаточный объем данных из входных объектов.
Чтобы вычислить распределение, местоположения источника могут иметь связанное значение, которое может быть задано параметром Поле источника (Source field). Если входные исходные данные - целочисленный растр, поле по умолчанию - VALUE. Если это пространственный объект, он будет первым целочисленным полем в таблице атрибутов. Если входные исходные данные – это растр с плавающей точкой, должен быть задан параметр растра целочисленного значения.
Если входные данные источника являются объектом, по умолчанию используется первое доступное поле. Если доступных полей нет, используется поле ObjectID (например, OID или FID, в зависимости от типа входных данных объекта).
Ячейки со значением NoData во Входном растре стоимости ведут себя как барьеры в инструментах стоимостной поверхности. Любой ячейке, которой присвоено значение NoData на входной поверхности стоимости, будет присвоено значение NoData на всех выходных растрах (стоимость расстояния, распределение и направление).
Если входные данные источников и растр стоимости имеют различные экстенты, по умолчанию выходной экстент будет определяться как пересечение двух растров. Чтобы получить поверхность стоимостного расстояния для всего экстента, выберите опцию Объединение входных данных (Union of Inputs) на выходных параметров среды Экстент (Extent).
Если Маска (Mask) была задана в параметрах среды, любые ячейки, попадающие под маску, будут обрабатываться как значения NoData.
Если в окне Анализ растров (Raster Analysis) была определена маска, и маскируемые участки содержат источник, вычисления, выполняемые для поверхности стоимости, будут осуществляться для оставшихся ячеек источников. Ячейки источников, попадающие под маску, не будут рассматриваться в вычислениях. Этим ячейкам будет присвоено значение NoData на всех выходных растрах (расстояния, распределения и направления).
Входной растр значений (Input value raster) полезен при использовании альтернативных значений или зон, или если исходные данные были получены из операции, которая приводит к двоичному результату 0 или 1, с потерей исходных значений зоны, связанных с этими местоположениями. Растр значений может либо восстанавливать эти значения, либо предусматривать анализ дополнительных комбинаций значений зон в пределах местоположений источников.
Если используется растр значений, он может менять конфигурацию и результаты выходных данных инструмента распределения по стоимостному расстоянию; он не будет влиять на дополнительные результирующие растры стоимостного расстояния и стоимостного направления.
Максимальное расстояние (Maximum distance) задаётся в тех же единицах стоимости, что используются и на входном растре стоимости.
Для выходного растра расстояния, расстояние с наименьшей стоимостью (или минимальное суммарное стоимостное расстояние) ячейки до набора исходных местоположений - это нижняя граница расстояний с наименьшей стоимостью от ячейки до всех исходных местоположений.
Синтаксис
| Параметр | Объяснение | Тип данных | |
in_source_data |
Входные местоположения источников. Растр или набор пространственных данных, определяющий ячейки или местоположения, до которых вычисляются стоимостные расстояния с наименьшей стоимостью перемещения от всех ячеек. Для растров, входной тип может быть целочисленным или с плавающей точкой. Если входной исходный растр – значение с плавающей точкой, {in_value_raster} должен быть целочисленного типа. Растр значений будет иметь преимущество над любым параметром {source_field}. | Raster Layer | Feature Layer | |
in_cost_raster |
Растр, определяющий полное сопротивление или стоимость движения в плоскости через каждую ячейку. Значение в каждой ячейке представляет стоимость движения через ячейку по поверхности стоимости, выраженную на единицу расстояния. Для получения общей стоимости прохождения через ячейку, значение каждой ячейки умножается на разрешение ячейки, с учетом возможности движения по диагонали. Значения на входном растре стоимости могут быть представлены целыми числами или числами с плавающей точкой, но они не могут быть отрицательными (стоимость не может быть отрицательной). | Raster Layer | |
maximum_distance (дополнительно) | Определяет пороговое значение суммарной стоимости перемещения, которое не может быть превышено. Если значение суммарного стоимостного расстояния превышает это значение, выходной ячейке будет присвоено значение NoData. Максимальное расстояние определяет экстент, для которого вычисляются совокупные стоимостные расстояния. Значение по умолчанию равно расстоянию до границы выходного растра. | Double | |
in_value_raster (дополнительно) |
Входной целочисленный растр, определяющий значения зон, которые должны быть использованы для каждого местоположения источника. Для каждого исходного местоположения (ячейки или пространственного объекта), значение, определенное по {in_value_raster}, будет присвоено всем ячейкам, относящимся к источнику, при вычислении. Растр значений будет иметь преимущество над параметром поля источников {source_field}. | Raster Layer | |
source_field (дополнительно) | Поле, используемое для присвоения значений исходным местоположениям. Должен быть целочисленным. Если был установлено {in_value_raster}, значения для этих входных данных будут иметь приоритет над любым параметром для {source_field}. | Field | |
out_distance_raster (дополнительно) |
Выходной растр стоимостного расстояния. Растр стоимостного расстояния определяет, для каждой ячейки, наименьшее суммарное стоимостное расстояние по поверхности стоимости до установленных положений источников. Источник может быть представлен ячейкой, набором ячеек или векторным объектом. Выходной растр представлен числами с плавающей точкой. | Raster Dataset | |
out_backlink_raster (дополнительно) |
Выходной растр стоимостного направления. Растр направления содержит значения от 0 до 8, которые определяют направление или следующую соседнюю ячейку вдоль оптимального пути с наименьшей суммарной стоимостью перемещения из каждой ячейки до ближайшего (в терминах стоимости) источника. Если маршрут должен пройти через правую соседнюю ячейку, ячейке будет присвоено значение 1, значение 2 будет присвоено нижней правой диагональной ячейке, и так далее по часовой стрелке. Значение 0 резервируется для ячеек источников.
| Raster Dataset |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_allocation_raster |
Выходной растр распределения по стоимостному расстоянию. Этот растр определяет зону каждого источника (ячейки или пространственного объекта), который может быть достигнут с наименьшими суммарными затратами. Выходной растр будет целочисленным | Raster |
Пример кода
На следующем скрипте Python Window показано, как использовать инструмент Распределение по стоимостному расстоянию (Cost Allocation).
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
out = ()
costAllocOut = CostAllocation("observers.shp", "costraster", 25000,"elevation",
"FID", "c:/sapyexamples/output/distout",
"c:/sapyexamples/output/backlinkout")
costAllocOut.save("c:/sapyexamples/output/costalloc")
Этот скрипт использует стоимостный растр, векторный слой исходных точек и несколько дополнительных параметров для вычисления растра ячеек, которые содержат значение ближайшего источника.
# Name: CostAllocation_Ex_02.py
# Description: Calculates for each cell its nearest source
# based on the least accumulative cost over a
# cost surface.
# 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
inFeature = "observers.shp"
costRaster = "costraster"
maxDist = 25000
valRaster = "elevation"
featField = "FID"
outDistanceRaster = "c:/sapyexamples/output/distout"
outBacklink = "c:/sapyexamples/output/backlinkout"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute CostAllocation
costAllocOut = CostAllocation(inFeature, costRaster, maxDist,
valRaster, featField, outDistanceRaster,
outBacklink)
# Save the output
costAllocOut.save("c:/sapyexamples/output/costalloc01")