Плотность (Kernel Density) (Spatial Analyst)
Резюме
Вычисляет величину (плотность) на единицу площади точеных или линейных объектов с использованием кернфункции для построения поверхности из сглаженных конусов для каждой точки или полилинии.
Более подробно о том, как работает инструмент Плотность ядер (Kernel Density)
Использование
-
Более высокие значения для параметра радиуса поиска приводят к построению более сглаженного и генерализованного растра плотности. Более низкие значения приводят к построению растра, на котором показано большее количество деталей.
-
При вычислении плотности учитываются только те точки или части линий, которые попадают в заданную область соседства. Если в область соседства конкретной ячейки не попадает ни одной точки или отрезка линии, такой ячейке будет присвоено значение NoData.
-
Если коэффициент преобразования для единиц измерения площади мал относительно объектов (расстояния между точками или длины сегментов линии, в зависимости от типа объектов), выходные значения могут быть очень маленькими. Чтобы получить более высокие значения, выберите коэффициент преобразования для более крупных единиц измерения площади (например, квадратные километры, а не квадратные метры).
Синтаксис
| Параметр | Объяснение | Тип данных |
in_features |
Входные объекты (точечные или линейные), для которых будет вычисляться плотность. | Feature Layer |
population_field |
Поле, в котором хранятся значения распределения (относительной величины показателя) для каждого объекта. Поле Population содержит значения количества или подсчета, которые должны быть распределены по ландшафту для построения непрерывной поверхности. Значения в поле численности могут быть целочисленными или с плавающей точкой. Опции и поведение по умолчанию для поля перечислены ниже.
| Field |
cell_size (дополнительно) | Размер ячейки для выходного набора растровых данных. Значение, устанавливаемое в параметрах среды. Если параметры среды не установлены, размер ячейки определяется, как более короткий размер экстента выходных объектов (ширины или высоты) в выходной пространственной привязке, деленный на 250. | Analysis Cell Size |
search_radius (дополнительно) |
Радиус поиска, в пределах которого будет вычислена плотность. Единицы измерения выбираются, основываясь на линейных единицах проекции пространственной привязки выходных данных. Например, если единицы в метрах, чтобы включить все пространственные объекты в окрестность в 1 милю, установите радиус поиска на 1609,344 (1 миля = 1609,344 метра). По умолчанию, используется значение меньшего из размеров ширины или высоты экстента выходных объектов на выходной пространственной привязке, деленное на 30. | Double |
area_unit_scale_factor (дополнительно) |
Желаемые единицы измерения площади выходных значений плотности. По умолчанию единицы измерения выбираются, основываясь на линейных единицах проекции пространственной привязки выходных данных. Можно изменить их на соответствующие единицы измерения, если вы хотите преобразовать выходные данные плотности. Значения для плотности линий преобразуют единицы измерения и длин, и площадей. Например, входные единицы измерения метры, по умолчанию единицы измерения плотности выходной области будут квадратные километры для точечных объектов или километр на квадратный километр для полилинейных объектов. Единицы плотности по умолчанию, на основе единиц измерения входных объектов:
| String |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_raster |
Выходной растр плотности. Это всегда растр с плавающей точкой. | Raster |
Пример кода
В этом примере вычисляется сглаженный растр плотности из точечного шейп-файла.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outKDens = KernelDensity("rec_sites.shp", "NONE", 45, 1200, "SQUARE_KILOMETERS")
outKDens.save("C:/sapyexamples/output/kdensout")
В этом примере вычисляется сглаженный растр плотности из точечного шейп-файла.
# Name: KernelDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from point or polyline
# features using a kernel function to fit a smoothly tapered
# surface to each point or polyline.
# 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
inFeatures = "rec_sites.shp"
populationField = "NONE"
cellSize = 60
searchRadius = 2500
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute KernelDensity
outKernelDensity = KernelDensity(inFeatures, populationField, cellSize,
searchRadius, "SQUARE_KILOMETERS")
# Save the output
outKernelDensity.save("C:/sapyexamples/output/kerneldout")