Densidad de kernel espaciotemporal (Spatial Analyst)

Disponible con una licencia de Spatial Analyst.

Resumen

Amplía los cálculos de densidad de kernel a partir del análisis de la posición relativa y la magnitud de las entidades de entrada para incluir otras dimensiones como el tiempo y la profundidad (elevación). El resultado identifica el área de magnitud por unidad utilizando múltiples funciones kernel para ajustar una superficie suavemente cónica a cada punto de entrada.

Más información sobre cómo funciona Densidad de kernel espaciotemporal

Uso

  • Debe proporcionar un valor para el parámetro Campo de elevación (elevation_field en Python) o el parámetro Campo de tiempo (time_field en Python), o ambos. Aunque la herramienta mostrará inicialmente ambos como obligatorios, una vez que se ha proporcionado un valor para uno de ellos, el otro pasa a ser opcional.

    Las funciones kernel se definen a partir de los valores de los parámetros proporcionados. Cuando se proporciona el valor del parámetro Campo de elevación, el kernel se extiende a lo largo del eje de dirección (z). Cuando se proporciona el valor del parámetro Campo de tiempo, el kernel se extiende a lo largo del eje temporal (t). Cuando se proporcionan ambos valores, los kernels se extienden a lo largo del eje z y del eje t.

    Si solo se proporciona uno de los valores de estos parámetros, la salida será un ráster tridimensional. Si se proporcionan los dos valores de los parámetros, la salida será un ráster de cuatro dimensiones.

  • Valores muy grandes o muy pequeños en el parámetro Campo de población (population_field en Python) pueden producir resultados que parezcan poco intuitivos. Si el valor medio del campo de población es mucho mayor que 1 (por ejemplo, como sucede con las poblaciones en las ciudades), el radio predeterminado de búsqueda puede ser muy pequeño, lo cual generaría pequeños anillos alrededor de los puntos de entrada. Si el valor medio del campo de población es mucho menor que 1, el radio de búsqueda calculado puede ser injustificadamente grande. En estos casos, puede proporcionar un radio de búsqueda personalizado.

  • El parámetro Tamaño de celda (cell_size en Python) puede definirse mediante un valor numérico u obtenerse a partir de un dataset ráster existente. Si no se ha proporcionado el tamaño de celda como valor del parámetro, se deriva del entorno Tamaño de celda si se ha especificado. Si no se ha proporcionado el tamaño de celda del parámetro o el tamaño de celda del entorno, pero se ha configurado el entorno Ráster de alineación, se utilizará el tamaño de celda del ráster de alineación. Si no se especifica nada, el tamaño de celda se calculará a partir de la anchura o la altura de la extensión dividida entre 250, donde la extensión se encuentra en el sistema de coordenadas de salida especificado en el entorno.

  • Si el tamaño de celda se especifica mediante un valor numérico, la herramienta lo utilizará directamente para el ráster de salida.

    Si el tamaño de celda se especifica utilizando un dataset ráster, el parámetro mostrará la ruta del dataset ráster en lugar del valor del tamaño de celda. El tamaño de celda de ese dataset ráster se utilizará directamente en el análisis, siempre que la referencia espacial del dataset sea la misma que la referencia espacial de salida. Si la referencia espacial del dataset es diferente de la referencia espacial de salida, se proyectará en función del valor de Método de proyección de tamaño de celda especificado.

  • El parámetro Valores resultantes (resultant_values en Python) especifica lo que representan los valores ráster de salida. Si se especifica Densidades (DENSITIES en Python), los valores representan el valor de la densidad del kernel por área de unidad para cada celda. Si se especifica Recuentos esperados (EXPECTED_COUNTS en Python), los valores representan la densidad del kernel por área de celda.

  • La opción Planar del parámetro Método (method en Python) resulta adecuada si el análisis se va a realizar a escala local con una proyección que mantiene el área y la distancia correctas de manera precisa. La opción Geodésica es adecuada si el análisis se va a realizar a escala regional o a gran escala (por ejemplo, con Web Mercator o cualquier sistema de coordenadas geográficas). Este método tiene en cuenta la curvatura del esferoide y trata correctamente los datos cercanos a los polos y a la línea internacional de cambio de fecha.

  • Solo los puntos que se encuentran dentro de una vecindad se tienen en cuenta en el cálculo de la densidad. Si ningún punto cae dentro de la vecindad de una celda en particular, NoData se asigna a esa celda.

  • Para los formatos de datos que admiten valores nulos, como las clases de entidades de geodatabase de archivos, se ignorará un valor nulo cuando se utilice como entrada.

  • Consulte Entornos de análisis y Spatial Analyst para obtener detalles adicionales sobre los entornos de geoprocesamiento válidos para esta herramienta.

  • Referencias:

    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.

Parámetros

EtiquetaExplicaciónTipo de datos
Entidades de puntos de entrada

Las entidades de punto de entrada cuya densidad se calculará.

Feature Layer
Campo de población

Campo que denota los valores de población para cada entidad. La población es el recuento o la cantidad que se expandirá a lo largo del paisaje para crear una superficie continua.

Los valores del campo de población pueden ser enteros o de punto flotante.

Use Ninguno si no se utilizará un elemento o valor especial y cada entidad se contará una sola vez.

Field
Campo de elevación
(Opcional)

El campo que denota los valores de elevación para cada entidad.

Los valores del campo de elevación pueden ser enteros o de punto flotante.

Utilice Vacío para admitir la densidad de kernel 3D con tiempo.

Para las entidades 3D, se agregará un pseudocampo, Shape.Z, a la lista de campos.

Field
Unidad del campo de elevación
(Opcional)

La unidad de medición que se utilizará para el valor del campo de elevación de entrada. La opción predeterminada es metros.

Utilice la unidad apropiada para representar los valores del valor del parámetro Campo de elevación.

  • PulgadaSe utilizarán las pulgadas.
  • PieSe utilizarán los pies.
  • YardaSe utilizarán las yardas.
  • Millas (EE. UU.)Se utilizarán las millas de EE. UU.
  • Milla náuticaSe utilizarán las millas náuticas.
  • MilímetroSe utilizarán los milímetros.
  • CentímetroSe utilizarán los centímetros.
  • ContadorSe utilizarán los metros.
  • KilómetrosSe utilizarán los kilómetros.
  • DecímetroSe utilizarán los decímetros.
String
Campo de tiempo
(Opcional)

El campo que denota los valores de tiempo de cada entidad.

Field
Tamaño de celda
(Opcional)

El tamaño de celda del ráster multidimensional de salida que se creará.

El valor se puede definir mediante un valor numérico u obtenerse a partir de un dataset ráster existente. Si el tamaño de celda no se ha proporcionado como valor de parámetro, se usa el valor del tamaño de celda del entorno si se ha especificado; de lo contrario, se utilizan reglas adicionales para calcularlo a partir de las otras entradas. Consulte el uso de la herramienta para más detalles.

Analysis Cell Size
Radio de búsqueda (x e y)
(Opcional)

El radio de búsqueda en el plano x,y dentro del cual se calculará la densidad.

Defina las unidades que se van a utilizar. Por ejemplo, para incluir todas las entidades dentro de una vecindad de 1 milla cuando las unidades son metros, fije el radio de búsqueda en 1609,344 (1 milla = 1609,344 metros).

Linear Unit
Radio de búsqueda (z)
(Opcional)

La distancia de búsqueda vertical en la dirección z dentro de la cual se calculará la densidad. Esta distancia vertical se utilizará para buscar entidades en las direcciones ascendente y descendente a lo largo del eje z.

Defina las unidades que se van a utilizar.

Linear Unit
Buscar en la ventana de tiempo (t)
(Opcional)

El rango de búsqueda de tiempo dentro del cual se calculará la densidad.

Defina las unidades que se van a utilizar.

Time Unit
Valores resultantes
(Opcional)

Especifica qué representan los valores del ráster de salida.

Dado que el valor de la celda de salida está vinculado al tamaño de celda especificado, el ráster resultante no se puede remuestrear con otro tamaño de celda.

  • DensidadesLos valores de salida representan el valor de densidad calculado por área de unidad para cada celda. Esta es la opción predeterminada.
  • Recuentos esperadosLos valores de salida representan el valor de densidad calculado por área de celda.
String
Método
(Opcional)

Especifica si se va a utilizar la tierra plana (planar) o la ruta más corta en un método de esferoide (geodésico).

  • PlanarSe utilizará la distancia planar entre entidades. Esta es la opción predeterminada.
  • GeodésicaSe utilizará la distancia geodésica entre entidades.
String
Elevación mínima
(Opcional)

La elevación inicial que se utilizará para la salida del ráster multidimensional.

Double
Elevación máxima
(Opcional)

La elevación final que se utilizará para la salida del ráster multidimensional.

Double
Intervalo de elevación
(Opcional)

El intervalo de elevación entre las divisiones de la salida del ráster multidimensional.

Double
Unidad de elevación
(Opcional)

Especifica la unidad del intervalo de elevación que se utilizará para la salida del ráster multidimensional. El valor predeterminado es metros.

  • PulgadaSe utilizarán las pulgadas.
  • PieSe utilizarán los pies.
  • YardaSe utilizarán las yardas.
  • Millas (EE. UU.)Se utilizarán las millas de EE. UU.
  • Milla náuticaSe utilizarán las millas náuticas.
  • MilímetroSe utilizarán los milímetros.
  • CentímetroSe utilizarán los centímetros.
  • ContadorSe utilizarán los metros.
  • KilómetrosSe utilizarán los kilómetros.
  • DecímetroSe utilizarán los decímetros.
String
Hora de inicio
(Opcional)

El tiempo inicial que se utilizará para la salida del ráster multidimensional.

Date
Hora de finalización
(Opcional)

El tiempo final que se utilizará para la salida del ráster multidimensional.

Date
Intervalo de tiempo
(Opcional)

El intervalo de tiempo entre las divisiones de la salida del ráster multidimensional.

Double
Unidad de intervalo de tiempo
(Opcional)

Especifica la unidad del intervalo de tiempo que se utilizará para la salida del ráster multidimensional. El valor predeterminado es día.

  • SegundoLa unidad de intervalo de tiempo será segundos.
  • MinutoLa unidad de intervalo de tiempo será minutos.
  • HoraLa unidad de intervalo de tiempo será horas.
  • DíaLa unidad de intervalo de tiempo será días.
  • SemanaLa unidad de intervalo de tiempo será semanas.
String

Valor de retorno

EtiquetaExplicaciónTipo de datos
Ráster de salida

El dataset ráster multidimensional de densidad del kernel de salida en formato de ráster de nube (.crf). En este momento, no se admite ningún otro formato de salida.

Siempre es un ráster de punto flotante.

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})
NombreExplicaciónTipo de datos
in_features

Las entidades de punto de entrada cuya densidad se calculará.

Feature Layer
population_field

Campo que denota los valores de población para cada entidad. La población es el recuento o la cantidad que se expandirá a lo largo del paisaje para crear una superficie continua.

Los valores del campo de población pueden ser enteros o de punto flotante.

Use '' si no se utilizará ningún elemento ni valor especial y cada entidad se contará una sola vez.

Field
elevation_field
(Opcional)

El campo que denota los valores de elevación para cada entidad.

Los valores del campo de elevación pueden ser enteros o de punto flotante.

Utilice '' para admitir la densidad de kernel 3D con tiempo.

Para las entidades 3D, se agregará un pseudocampo, Shape.Z, a la lista de campos.

Field
elevation_field_unit
(Opcional)

La unidad de medición que se utilizará para el valor del campo de elevación de entrada. La opción predeterminada es metros.

Utilice la unidad apropiada para representar los valores del valor del parámetro elevation_field.

  • INCHSe utilizarán las pulgadas.
  • FOOTSe utilizarán los pies.
  • YARDSe utilizarán las yardas.
  • MILE_USSe utilizarán las millas de EE. UU.
  • NAUTICAL_MILESe utilizarán las millas náuticas.
  • MILLIMETERSe utilizarán los milímetros.
  • CENTIMETERSe utilizarán los centímetros.
  • METERSe utilizarán los metros.
  • KILOMETERSe utilizarán los kilómetros.
  • DECIMETERSe utilizarán los decímetros.
String
time_field
(Opcional)

El campo que denota los valores de tiempo de cada entidad.

Field
cell_size
(Opcional)

El tamaño de celda del ráster multidimensional de salida que se creará.

El valor se puede definir mediante un valor numérico u obtenerse a partir de un dataset ráster existente. Si el tamaño de celda no se ha proporcionado como valor de parámetro, se usa el valor del tamaño de celda del entorno si se ha especificado; de lo contrario, se utilizan reglas adicionales para calcularlo a partir de las otras entradas. Consulte el uso de la herramienta para más detalles.

Analysis Cell Size
kernel_search_radius_xy
(Opcional)

El radio de búsqueda en el plano x,y dentro del cual se calculará la densidad.

Defina las unidades que se van a utilizar. Por ejemplo, para incluir todas las entidades dentro de una vecindad de 1 milla cuando las unidades son metros, fije el radio de búsqueda en 1609,344 (1 milla = 1609,344 metros).

Linear Unit
kernel_search_radius_z
(Opcional)

La distancia de búsqueda vertical en la dirección z dentro de la cual se calculará la densidad. Esta distancia vertical se utilizará para buscar entidades en las direcciones ascendente y descendente a lo largo del eje z.

Defina las unidades que se van a utilizar.

Linear Unit
kernel_search_time_window
(Opcional)

El rango de búsqueda de tiempo dentro del cual se calculará la densidad.

Defina las unidades que se van a utilizar.

Time Unit
resultant_values
(Opcional)

Especifica qué representan los valores del ráster de salida.

Dado que el valor de la celda de salida está vinculado al tamaño de celda especificado, el ráster resultante no se puede remuestrear con otro tamaño de celda.

  • DENSITIESLos valores de salida representan el valor de densidad calculado por área de unidad para cada celda. Esta es la opción predeterminada.
  • EXPECTED_COUNTSLos valores de salida representan el valor de densidad calculado por área de celda.
String
method
(Opcional)

Especifica si se va a utilizar la tierra plana (planar) o la ruta más corta en un método de esferoide (geodésico).

  • PLANARSe utilizará la distancia planar entre entidades. Esta es la opción predeterminada.
  • GEODESICSe utilizará la distancia geodésica entre entidades.
String
min_elevation
(Opcional)

La elevación inicial que se utilizará para la salida del ráster multidimensional.

Double
max_elevation
(Opcional)

La elevación final que se utilizará para la salida del ráster multidimensional.

Double
elevation_interval
(Opcional)

El intervalo de elevación entre las divisiones de la salida del ráster multidimensional.

Double
elevation_unit
(Opcional)

Especifica la unidad del intervalo de elevación que se utilizará para la salida del ráster multidimensional. El valor predeterminado es metros.

  • INCHSe utilizarán las pulgadas.
  • FOOTSe utilizarán los pies.
  • YARDSe utilizarán las yardas.
  • MILE_USSe utilizarán las millas de EE. UU.
  • NAUTICAL_MILESe utilizarán las millas náuticas.
  • MILLIMETERSe utilizarán los milímetros.
  • CENTIMETERSe utilizarán los centímetros.
  • METERSe utilizarán los metros.
  • KILOMETERSe utilizarán los kilómetros.
  • DECIMETERSe utilizarán los decímetros.
String
start_time
(Opcional)

El tiempo inicial que se utilizará para la salida del ráster multidimensional.

Date
end_time
(Opcional)

El tiempo final que se utilizará para la salida del ráster multidimensional.

Date
time_interval
(Opcional)

El intervalo de tiempo entre las divisiones de la salida del ráster multidimensional.

Double
time_interval_unit
(Opcional)

Especifica la unidad del intervalo de tiempo que se utilizará para la salida del ráster multidimensional. El valor predeterminado es día.

  • SECONDLa unidad de intervalo de tiempo será segundos.
  • MINUTELa unidad de intervalo de tiempo será minutos.
  • HOURLa unidad de intervalo de tiempo será horas.
  • DAYLa unidad de intervalo de tiempo será días.
  • WEEKLa unidad de intervalo de tiempo será semanas.
String

Valor de retorno

NombreExplicaciónTipo de datos
out_raster

El dataset ráster multidimensional de densidad del kernel de salida en formato de ráster de nube (.crf). En este momento, no se admite ningún otro formato de salida.

Siempre es un ráster de punto flotante.

Raster

Muestra de código

Ejemplo 1 de SpaceTimeKernelDensity (ventana de Python)

Este ejemplo calcula una densidad kernel temporal suavizada utilizando un shapefile.

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")
Ejemplo 2 de SpaceTimeKernelDensity (script independiente)

Este ejemplo calcula una densidad kernel temporal suavizada utilizando un dataset multidimensional.

## 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")

Información de licenciamiento

  • Basic: Requiere Spatial Analyst
  • Standard: Requiere Spatial Analyst
  • Advanced: Requiere Spatial Analyst

Temas relacionados