Cómo funciona la densidad kernel

Disponible con una licencia de Spatial Analyst.

La herramienta Densidad kernel calcula la densidad de las entidades en la vecindad de esas entidades. Puede calcularse para las entidades de punto y de línea.

Los posibles usos incluyen analizar la densidad de las casas o delitos para la planificación social, así como descubrir cómo las líneas de carreteras o de servicios públicos influyen en el hábitat natural. El campo de población se puede utilizar para ponderar algunas entidades más que otras o para permitir que un punto represente varias observaciones. Por ejemplo, una dirección puede representar un condominio con seis unidades, o algunos crímenes se pueden ponderar más que otros al determinar los niveles de crímenes en general. En el caso de las entidades de línea, es posible que una autovía tenga más impacto que un camino estrecho sin asfaltar.

Cómo se calcula la densidad kernel

Para entidades de puntos

La Densidad kernel calcula la densidad de las entidades de punto de alrededor de cada celda ráster de salida.

Conceptualmente, se ajusta una superficie curva uniforme sobre cada punto. El valor de superficie es más alto en la ubicación del punto y disminuye a medida que aumenta la distancia desde el punto, alcanza cero en la distancia Radio de búsqueda desde el punto. Sólo es posible un vecindario circular. El volumen bajo la superficie es igual al valor Campo de población del punto, o 1 si se especifica NONE. Para calcular la densidad de cada celda ráster de salida, se agregan los valores de todas las superficies de kernel en donde se superponen con el centro de la celda ráster. La función kernel se basa en la función kernel cuártica que se describió en Silverman (1986, p. 76, ecuación 4.5).

Si se utiliza una configuración de campo de población que no sea NONE, el valor de cada elemento determina la cantidad de veces que se va a contar el punto. Por ejemplo, un valor de 3 haría que el punto se cuente como tres puntos. Los valores pueden ser de tipo entero o punto flotante.

De forma predeterminada, una unidad se selecciona según la unidad lineal de la definición de proyección de los datos de la entidad de punto de entrada o como se haya especificado en la configuración del entorno del Sistema de coordenadas de salida. Si se selecciona una unidad de área, la densidad que se calculó para la celda se multiplica por el factor apropiado antes de que se escriba en el ráster de salida.

Por ejemplo, si las unidades de entrada son metros, las unidades de área de salida serán kilómetros cuadrados, de manera predeterminada. Comparar un factor de escala de unidades de metros a kilómetros tendrá como resultado que los valores serán diferentes por un multiplicador de 1.000.000 (1.000 metros x 1.000 metros).

Para entidades de línea

La Densidad kernel también puede calcular la densidad de las entidades lineales en la vecindad de cada celda ráster de salida.

Conceptualmente, se ajusta una superficie curva uniforme sobre cada línea. Su valor es más grande en la línea y disminuye a medida que se aleja de dicha línea, alcanza cero en la distancia Radio de búsqueda que se especificó desde la línea. La superficie se define de modo que el volumen bajo la superficie sea igual al producto de la longitud de línea y el valor Campo de población. Para calcular la densidad de cada celda ráster de salida, se agregan los valores de todas las superficies de kernel en donde se superponen con el centro de la celda ráster. El uso de la función kernel para las líneas se adapta desde la función de kernel de cuarto grado para las densidades de punto como se describe en Silverman.

Un segmento de línea y la superficie kernel que se ajusta sobre éste.

La ilustración anterior muestra un segmento de línea y la superficie kernel que se ajusta sobre el mismo. La contribución del segmento de línea a la densidad es igual al valor de la superficie de kernel en el centro de la celda ráster.

De forma predeterminada, una unidad se selecciona según la unidad lineal de la definición de proyección de los datos de la entidad de polilínea de entrada o como se haya especificado en la configuración del entorno del Sistema de coordenadas de salida.

Cuando se especifica un factor de Unidades de área de salida, este convierte las unidades de longitud y de área. Por ejemplo, si la unidad lineal es metros, las unidades de área de salida estarán, de manera predeterminada, en SQUARE_KILOMETERS y las unidades de densidad de línea resultantes se convertirán a kilómetros por kilómetro cuadrado. El resultado final, al comparar un factor de escala de área de metros a kilómetros, será que los valores de densidad serán diferentes por un multiplicador de 1.000.

Puede controlar las unidades de densidad al seleccionar manualmente el factor apropiado. Para establecer que la densidad esté en metros por metro cuadrado (en lugar del valor predeterminado de kilómetros por kilómetro cuadrado), establezca las unidades de área en SQUARE_METERS. Del mismo modo, para que las unidades de densidad de la salida estén en millas por milla cuadrada, establezca las unidades de área en SQUARE_MILES.

Consulte el siguiente tema para obtener más información sobre unidades de distancia específicas.

• Más información sobre unidades y constantes de conversión en ArcGIS

Si se utiliza un campo de población que no sea NONE, se considera que la longitud de la línea es su longitud real multiplicada por el valor del campo de población para esa línea.

Fórmulas para calcular la densidad kernel

Estas fórmulas definen cómo se calcula la densidad kernel para puntos y cómo se determina el radio de búsqueda predeterminado en la fórmula de densidad kernel.

Predecir la densidad de puntos

La densidad prevista en una nueva ubicación (x, y) se determina con esta fórmula:

donde:

• i = 1,…,n son los puntos de entrada. Solo incluya puntos en la suma si están dentro de la distancia de radio de la ubicación (x, y).
• pop_i es el valor de campo de población del punto I, que es un parámetro opcional.
• dist_i es la distancia entre el punto i y la ubicación (x, y).

Después, la densidad calculada se multiplica por el número de puntos o la suma del campo de población, si se proporcionó alguno. Esta corrección hace que el integral espacial sea igual que el número de puntos (o la suma o el campo de población) en lugar de que siempre sea igual que 1. Esta implementación usa una función kernel cuártica (Silverman, 1986). Será necesario calcular la fórmula para cada ubicación en la que desea estimar la densidad. Puesto que se está creando un ráster, los cálculos se aplican al centro de cada celda en el ráster de salida.

Radio de búsqueda predeterminado (ancho de banda)

En ArcGIS 10.2.1 se ha mejorado la forma de cálculo del radio predeterminado de búsqueda.

El algoritmo que se usa para determinar el radio predeterminado de búsqueda, también conocido como ancho de banda, es el siguiente:

1. Calcular el centro medio de los puntos de entrada. Si se proporcionó un campo de Población, este cálculo y los siguientes se ponderarán con los valores de ese campo.
2. Calcular la distancia desde el centro medio (ponderado) para todos los puntos.
3. Calcular la mediana (ponderada) de estas distancias, D_m.
4. Calcular el valor de la Distancia estándar (ponderada), SD.

   Consulte la herramienta de Estadísticas espaciales Distancia estándar para obtener más detalles sobre el tema.

5. Aplique esta fórmula para calcular el ancho de banda.

donde:

• D_m es la distancia mediana (ponderada) desde el centro medio (ponderado).
• n es el número de puntos si no se usa ningún campo de población, o si se proporciona un campo de población, la suma de los valores de campo de población.
• SD es la distancia estándar.

Tenga en cuenta que la parte min de la ecuación significa que se utilizará la opción cuyo resultado sea el menor valor de cualquiera de las dos, ya sea SD u .

Existen dos métodos para calcular la distancia estándar, Sin ponderar y Ponderado.

Distancia sin ponderar

donde:

• x _i , y _i and z _i son las coordenadas de la entidad i
• {x̄, ȳ, z̄} representa el centro medio de las entidades
• n es igual al número total de entidades.

Distancia ponderada

donde:

• w_i es la ponderación de la entidad i
• {x _w, y _w, z _w} representa el centro medio ponderado

Metodología

Esta metodología de seleccionar el radio de búsqueda se basa en la fórmula de estimación de ancho de banda de la regla general de Silverman, pero se ha adaptado a las dos dimensiones. Por lo general, esta forma de calcular un radio predeterminado evita el fenómeno de generar "anillos alrededor de puntos" que a menudo se producía con datasets de puntos dispersos, y resiste a valores atípicos espaciales –algunos puntos que están alejados del resto de puntos–.

Recreación de los resultados de versiones anteriores de ArcGIS

En versiones anteriores de ArcGIS a la 10.2.1, el radio de búsqueda predeterminado se calculaba a partir del ancho o la altura menores de la extensión de entrada, dividido entre 30. Si desea reproducir esos resultados, multiplique el valor predeterminado del parámetro Tamaño de celda de salida por 25/3 o 8,333333.

Referencias

Silverman, B. W. Estimación de densidad para las estadísticas y el análisis de datos. New York: Chapman and Hall, 1986.