Плотность ядер пространства-времени (Spatial Analyst)

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

Расширяет вычисления плотности ядер от анализа относительного положения и величины входных объектов до включения других измерений, таких как время и глубина (высота). Итоговые выходные данные идентифицируют величину на единицу площади, используя несколько функций ядра, чтобы вписать плавно сужающуюся поверхность в каждую входную точку.

Более подробно о том, как работает инструмент Плотность ядер пространства-времени

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

  • Необходимо предоставить значение параметра Поле высоты (elevation_field в Python) или Поле времени (time_field в Python), или обоих параметров. Хотя исходно инструмент показывает оба параметра как обязательные, как только предоставляется один из них, другой становится дополнительным.

    Функции ядра задаются на основе предоставленных значений параметра. Когда предоставлено значение параметра Поле высоты, ядро расширяется вдоль оси направления (z). Когда предоставлено значение параметра Поле времени, ядро расширяется вдоль оси времени (t). Если предоставлены оба значения, ядра расширяются вдоль z-оси и t-оси.

    Если предоставлен только один из этих параметров, выходные данные будут трехмерным растром. Если предоставлены оба значения параметров, выходные данные будут четырехмерным растром.

  • Очень большие или очень малые значения в параметре Поле генеральной совокупности (population_field в Python) могут привести к результатам, которые кажутся не интуитивными. Если среднее значение для поля генеральной совокупности намного больше 1 (например, для численности городского населения), радиус поиска по умолчанию может оказаться очень малым, в результате чего будут возникать небольшие кольца вокруг входных точек. Если среднее значение для поля генеральной совокупности намного меньше 1, то расчетный радиус поиска может показаться неоправданно большой. В таких случаях, вы указываете пользовательский радиус поиска.

  • Параметр Размер ячейки (cell_size в Python) может быть задан числовым значениям или получен из имеющегося набора растровых данных. Если размер ячейки не был задан как значение параметра, он извлекается из параметра среды Размер ячейки, если такой указан. Если же размер ячейки не указан ни как параметр, ни как параметр среды, но настроен параметр среды Растр привязки, то будет использоваться значение величины его ячейки. Если же ни один из параметров не задан, то размер ячейки будет вычислен как наименьшее значение ширины или высоты экстента растра, поделенное на 250, при этом экстент берется из параметра среды Выходная система координат.

  • Если размер ячейки задан с использованием числового значения, инструмент будет использовать его непосредственно для выходного растра.

    Если размер ячейки задан с использованием набора растровых данных, параметр покажет путь к набору растровых данных вместо значения размера ячейки. Размер ячейки этого набора растровых данных будет использоваться непосредственно в анализе при условии, что пространственная привязка набора данных совпадает с выходной пространственной привязкой. Если пространственная привязка набора данных отличается от выходной пространственной привязки, она будет спроецирована на основе значения Метод проецирования размера ячейки.

  • Параметр Результирующие значения (resultant_values в Python) указывает, что будут показывать значения выходного растра. Если задано Плотности (DENSITIES в Python), то значения будут показывать плотность ядер на единицу площади для каждой ячейки. Если задано Ожидаемое количество (EXPECTED_COUNTS в Python), то значения будут показывать плотность ядер на площадь ячейки.

  • Опция Плоскостной в параметре Метод (method в Python) может быть целесообразной, если анализ будет выполнен в локальном масштабе с проекцией, которая точно поддерживает правильные расстояния и площади. Опция Геодезический подходит, если анализ будет выполняться в масштабе региона или других больших площадей (например, с использованием проекции Web Mercator или любой географической системы координат). Этот метод учитывает кривизну сфероида и корректно обрабатывает данные вблизи полюсов и международной линии смены дат.

  • При вычислении плотности учитываются только те точки, которые попадают в заданную окрестность. Если в окрестность определенной ячейки не попадает никаких точек, такой ячейке будет присвоено значение NoData.

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

  • См. раздел Параметры среды анализа и Spatial Analyst для получения дополнительной информации о среде геообработки данного инструмента.

  • Ссылки:

    Härdle, W. K., Müller, M, Sperlich, S., and Werwatz, A. Nonparametric and semiparametric models (Vol. 1). Berlin: Springer, 2004.

    Hu, Y., Wang, F., Guin, C., and Zhu, H. "A spatio-temporal kernel density estimation framework for predictive crime hotspot mapping and evaluation." Applied geography, 99, 2018, 89-97.

    Nakaya, T., and Yano, K. "Visualising crime clusters in a space‐time cube: An exploratory data analysis approach using space time kernel density estimation and scan statistics." Transactions in GIS, 14(3), 2010, 223-239.

    Silverman, B. W. Density Estimation for Statistics and Data Analysis. New York: Chapman and Hall, 1986.

Параметры

ПодписьОписаниеТип данных
Входные точечные объекты

Входные точечные объекты, для которых будет рассчитываться плотность.

Feature Layer
Поле генеральной совокупности

Поле, в котором хранятся значения распределения (относительной величины показателя) для каждого объекта. Поле генеральной совокупности содержит значения количества или подсчета, которые должны быть распределены по ландшафту для построения непрерывной поверхности.

Значения в поле генеральной совокупности могут быть целочисленными или с плавающей точкой.

Используйте значение Нет, если не будет использовано какого-либо показателя или специального значения, и каждый объект будет подсчитан один раз.

Field
Поле высоты
(Дополнительный)

Поле, обозначающее значения высоты для каждого объекта.

Значения в поле высот могут быть целочисленными или с плавающей точкой.

Используйте Пусто, чтобы поддерживать плотность ядер 3D с временем.

Для 3D-объектов в список полей будет добавлено псевдо поле Shape.Z.

Field
Единицы измерения поля высот
(Дополнительный)

Единицы измерения, которые будут использоваться для входного значения поля высот. По умолчанию – Метры.

Используйте соответствующие единицы измерений для представления значений параметра Поле значений.

  • ДюймБудут использоваться дюймы.
  • ФутБудут использоваться футы.
  • ЯрдБудут использоваться ярды.
  • Миля (США)Будут использоваться мили США.
  • Морская миляБудут использоваться морские мили.
  • МиллиметрБудут использоваться миллиметры.
  • СантиметрБудут использоваться сантиметры.
  • МетрБудут использоваться метры.
  • КилометрыБудут использоваться километры.
  • ДециметрБудут использоваться дециметры.
String
Поле времени
(Дополнительный)

Поле, обозначающее значения времени для каждого объекта.

Field
Размер ячейки
(Дополнительный)

Размер ячейки создаваемого выходного многомерного растра.

Значение можно задать как числовым значением, так и получить из существующего набора растровых данных. Если размер ячейки не указан в качестве значения параметра, будет использоваться значение размера ячейки в параметрах среды, если оно указано; в противном случае для его вычисления на основе других входных данных будут использоваться дополнительные правила. Подробности смотрите в разделе использования инструмента.

Analysis Cell Size
Радиус поиска (x и y)
(Дополнительный)

Радиус поиска в плоскости x,y, в пределах которого будет вычислена плотность.

Задайте единицы измерения для использования. Например, чтобы включить все объекты в пределах окрестности в одну милю, когда единицами измерения являются метры, задайте радиус поиска равным 1609.344 (1 миля = 1609.344 метра).

Linear Unit
Радиус поиска (z)
(Дополнительный)

Вертикальное расстояние поиска в z- направлении, в пределах которого будет вычислена плотность. Это вертикальное расстояние будет использоваться при поиске объектов в направлении вверх и вниз по оси z.

Задайте единицы измерения для использования.

Linear Unit
Временное окно поиска (t)
(Дополнительный)

Диапазон поиска во времени, в пределах которого будет вычислена плотность.

Задайте единицы измерения для использования.

Time Unit
Результирующие значения
(Дополнительный)

Определяет, что представляют значения в выходном растре.

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

  • ПлотностьВыходные значения будут показывать вычисленную плотность на единицу площади для каждой ячейки. Это значение по умолчанию
  • Ожидаемое количествоВыходные значения представляют вычисленное значение плотности на площадь ячейки.
String
Метод
(Дополнительный)

Определяет, какой метод будет использоваться: на плоской поверхности земли (планарный) или кратчайший путь на сфероиде (геодезический).

  • ПлоскостнойБудет использоваться расстояние на плоскости между объектами. Это значение по умолчанию
  • ГеодезическийБудет использоваться геодезическое расстояние между объектами.
String
Минимальная высота
(Дополнительный)

Начальная высота, которая будет использоваться для выходного многомерного растра.

Double
Максимальная высота
(Дополнительный)

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

Double
Интервал высот
(Дополнительный)

Шаг высот между срезами в выходном многомерном растре.

Double
Единицы измерения высот
(Дополнительный)

Задает единицы измерения шага высот, которые будут использоваться для выходного многомерного растра. По умолчанию метры.

  • ДюймБудут использоваться дюймы.
  • ФутБудут использоваться футы.
  • ЯрдБудут использоваться ярды.
  • Миля (США)Будут использоваться мили США.
  • Морская миляБудут использоваться морские мили.
  • МиллиметрБудут использоваться миллиметры.
  • СантиметрБудут использоваться сантиметры.
  • МетрБудут использоваться метры.
  • КилометрыБудут использоваться километры.
  • ДециметрБудут использоваться дециметры.
String
Время начала
(Дополнительный)

Начальное время, которое будет использоваться для выходного многомерного растра.

Date
Время окончания
(Дополнительный)

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

Date
Временной интервал
(Дополнительный)

Временной интервал между срезами в выходном многомерном растре.

Double
Единицы измерения временного интервала
(Дополнительный)

Задает единицы измерения интервала времени, которые будут использоваться для выходного многомерного растра. По умолчанию - день.

  • СекундаЕдиницами временного интервала будут секунды.
  • МинутаЕдиницами временного интервала будут минуты.
  • ЧасЕдиницами временного интервала будут часы.
  • ДеньЕдиницами временного интервала будут дни.
  • НеделяЕдиницами временного интервала будут недели.
String

Возвращаемое значение

ПодписьОписаниеТип данных
Выходной растр

Выходной набор многомерных растровых данных плотности ядер в формате Cloud Raster Format (.crf). Другие форматы выходных данных в настоящее время не поддерживаются.

Это всегда растр с плавающей точкой.

Raster

SpaceTimeKernelDensity(in_features, population_field, {elevation_field}, {elevation_field_unit}, {time_field}, {cell_size}, {kernel_search_radius_xy}, {kernel_search_radius_z}, {kernel_search_time_window}, {resultant_values}, {method}, {min_elevation}, {max_elevation}, {elevation_interval}, {elevation_unit}, {start_time}, {end_time}, {time_interval}, {time_interval_unit})
ИмяОписаниеТип данных
in_features

Входные точечные объекты, для которых будет рассчитываться плотность.

Feature Layer
population_field

Поле, в котором хранятся значения распределения (относительной величины показателя) для каждого объекта. Поле генеральной совокупности содержит значения количества или подсчета, которые должны быть распределены по ландшафту для построения непрерывной поверхности.

Значения в поле генеральной совокупности могут быть целочисленными или с плавающей точкой.

Используйте значение '', если не будет использовано какого-либо показателя или специального значения, и каждый объект будет подсчитан один раз.

Field
elevation_field
(Дополнительный)

Поле, обозначающее значения высоты для каждого объекта.

Значения в поле высот могут быть целочисленными или с плавающей точкой.

Используйте '', чтобы поддерживать плотность ядер 3D с временем.

Для 3D-объектов в список полей будет добавлено псевдо поле Shape.Z.

Field
elevation_field_unit
(Дополнительный)

Единицы измерения, которые будут использоваться для входного значения поля высот. По умолчанию – Метры.

Используйте соответствующие единицы измерений для представления значений параметра elevation_field.

  • INCHБудут использоваться дюймы.
  • FOOTБудут использоваться футы.
  • YARDБудут использоваться ярды.
  • MILE_USБудут использоваться мили США.
  • NAUTICAL_MILEБудут использоваться морские мили.
  • MILLIMETERБудут использоваться миллиметры.
  • CENTIMETERБудут использоваться сантиметры.
  • METERБудут использоваться метры.
  • KILOMETERБудут использоваться километры.
  • DECIMETERБудут использоваться дециметры.
String
time_field
(Дополнительный)

Поле, обозначающее значения времени для каждого объекта.

Field
cell_size
(Дополнительный)

Размер ячейки создаваемого выходного многомерного растра.

Значение можно задать как числовым значением, так и получить из существующего набора растровых данных. Если размер ячейки не указан в качестве значения параметра, будет использоваться значение размера ячейки в параметрах среды, если оно указано; в противном случае для его вычисления на основе других входных данных будут использоваться дополнительные правила. Подробности смотрите в разделе использования инструмента.

Analysis Cell Size
kernel_search_radius_xy
(Дополнительный)

Радиус поиска в плоскости x,y, в пределах которого будет вычислена плотность.

Задайте единицы измерения для использования. Например, чтобы включить все объекты в пределах окрестности в одну милю, когда единицами измерения являются метры, задайте радиус поиска равным 1609.344 (1 миля = 1609.344 метра).

Linear Unit
kernel_search_radius_z
(Дополнительный)

Вертикальное расстояние поиска в z- направлении, в пределах которого будет вычислена плотность. Это вертикальное расстояние будет использоваться при поиске объектов в направлении вверх и вниз по оси z.

Задайте единицы измерения для использования.

Linear Unit
kernel_search_time_window
(Дополнительный)

Диапазон поиска во времени, в пределах которого будет вычислена плотность.

Задайте единицы измерения для использования.

Time Unit
resultant_values
(Дополнительный)

Определяет, что представляют значения в выходном растре.

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

  • DENSITIESВыходные значения будут показывать вычисленную плотность на единицу площади для каждой ячейки. Это значение по умолчанию
  • EXPECTED_COUNTSВыходные значения представляют вычисленное значение плотности на площадь ячейки.
String
method
(Дополнительный)

Определяет, какой метод будет использоваться: на плоской поверхности земли (планарный) или кратчайший путь на сфероиде (геодезический).

  • PLANARБудет использоваться расстояние на плоскости между объектами. Это значение по умолчанию
  • GEODESICБудет использоваться геодезическое расстояние между объектами.
String
min_elevation
(Дополнительный)

Начальная высота, которая будет использоваться для выходного многомерного растра.

Double
max_elevation
(Дополнительный)

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

Double
elevation_interval
(Дополнительный)

Шаг высот между срезами в выходном многомерном растре.

Double
elevation_unit
(Дополнительный)

Задает единицы измерения шага высот, которые будут использоваться для выходного многомерного растра. По умолчанию метры.

  • INCHБудут использоваться дюймы.
  • FOOTБудут использоваться футы.
  • YARDБудут использоваться ярды.
  • MILE_USБудут использоваться мили США.
  • NAUTICAL_MILEБудут использоваться морские мили.
  • MILLIMETERБудут использоваться миллиметры.
  • CENTIMETERБудут использоваться сантиметры.
  • METERБудут использоваться метры.
  • KILOMETERБудут использоваться километры.
  • DECIMETERБудут использоваться дециметры.
String
start_time
(Дополнительный)

Начальное время, которое будет использоваться для выходного многомерного растра.

Date
end_time
(Дополнительный)

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

Date
time_interval
(Дополнительный)

Временной интервал между срезами в выходном многомерном растре.

Double
time_interval_unit
(Дополнительный)

Задает единицы измерения интервала времени, которые будут использоваться для выходного многомерного растра. По умолчанию - день.

  • SECONDЕдиницами временного интервала будут секунды.
  • MINUTEЕдиницами временного интервала будут минуты.
  • HOURЕдиницами временного интервала будут часы.
  • DAYЕдиницами временного интервала будут дни.
  • WEEKЕдиницами временного интервала будут недели.
String

Возвращаемое значение

ИмяОписаниеТип данных
out_raster

Выходной набор многомерных растровых данных плотности ядер в формате Cloud Raster Format (.crf). Другие форматы выходных данных в настоящее время не поддерживаются.

Это всегда растр с плавающей точкой.

Raster

Пример кода

SpaceTimeKernelDensity, пример 1 (окно Python)

В этом примере вычисляется сглаженная временная плотность ядер с использованием шейп-файла.

from arcpy import env  
from arcpy.sa import * 

env.workspace = "C:/sapyexamples/data" 
STKD_out_raster = SpaceTimeKernelDensity("WOD_subset.shp", "Salinity", "Z", "Meter",
                                         "Time", "0.001", resultant_values="Densities",
                                         method="Planar", elevation_unit="Meter")  

STKD_out_raster.save("C:/sapyexamples/output/STKD_out.crf")
SpaceTimeKernelDensity, пример 2 (автономный скрипт)

В этом примере вычисляется сглаженная временная плотность ядер с использованием многомерного набора данных.

## Name: SpaceTimeKernelDensity_Ex_standalone.py  
## Description: Calculate spatial temporal salinity concentration using a multidimensional dataset 
## Requirements: Spatial Analyst Extension 
 
## Import system modules 
import arcpy  
from arcpy import env   
from arcpy.sa import *

## Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
 
## Set environment settings 
env.workspace = r" C:\STKD_Test"
# To allow overwriting outputs change overwriteOutput option to True. 
env.overwriteOutput = False 
  
## Set local variables 
in_features = "WOD_subset"  
Population_Field = "Salinity"  
Elevation_Field = "Z"  
Elevation_Field_Unit = "Meter"  
Time_Field = "Time"  
Cell_Size = "30"  
Resultant_values = "Densities"  
Method = "Planar" 
Elevation_Unit = "Meter"  
  
## Execute: Space Time Kernel Density  
STKD_out_raster = SpaceTimeKernelDensity(in_features, Population_Field,   
                                Elevation_Field, Elevation_Field_Unit,   
                                Time_Field, Cell_Size,   
                                resultant_values=Resultant_values,   
                                method=Method, 
                                elevation_unit=Elevation_Unit) 
  
## Save the output 
STKD_out_raster.save("STKD_test.crf")

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