Calculando la distancia mínima entre puntos en QGIS

Question

Espero calcular la distancia entre cada propiedad en el área metropolitana de Washington DC y la estación de metro más cercana, y estoy haciendo esto para cada mes desde 2000 hasta 2017. Utilizo Distance Matrix en QGIS 3, establezco la capa de propiedades como la Capa de Entrada y la capa de estaciones de metro como la Capa de Destino, y utilizo solo el punto de destino más cercano. Desde 2000 hasta 2017, varias estaciones de metro nuevas se abrieron en el área metropolitana de DC mientras que ninguna estación fue cerrada.

Los resultados son en su mayoría buenos. La mayoría de las propiedades tienen las mismas estaciones de metro más cercanas a lo largo del período y las distancias permanecen iguales. Algunas propiedades tienen estaciones de metro más cercanas diferentes a medida que pasa el tiempo debido a la introducción de nuevas estaciones, y las distancias disminuyen a medida que se abren nuevas estaciones.

Sin embargo, después de la apertura de nuevas estaciones, una pequeña parte de las propiedades (aunque aún en gran número en términos absolutos) son asignadas a una nueva estación que está más lejos que la estación asignada a ellas antes de la apertura de nuevas estaciones. Por ejemplo, la Estación McLean de la Línea Plateada abrió en julio de 2014. Para un grupo de propiedades en mis datos, QGIS identifica la Estación Rockville de la Línea Roja como su estación de metro más cercana antes de julio de 2014. Después de julio de 2014, QGIS identifica la Estación McLean como su estación más cercana, pero las distancias entre estas propiedades y la Estación McLean son mayores que las distancias entre ellas y la Estación Rockville.

Revisé manualmente alrededor de 20-30 de estos casos y no encontré ningún error en las distancias calculadas por QGIS, pero ¿por qué QGIS no asigna estas propiedades a la Estación Rockville que sigue siendo la estación más cercana incluso después de la apertura de la Estación McLean?

Answer 1

El algoritmo de Matriz de distancias actualmente utiliza un filtro espacial cartesiano basado en el SRC de la capa de destino para encontrar los k puntos más cercanos, y calcula la distancia elipsoidal solo para esos puntos. [La Matriz de Distancias devuelve un registro/resultato incorrecto al usar Matriz de Distancias Lineal con k=1]

Esto suele ser problemático. Y es especialmente problemático cuando k=1 (el punto más cercano en coordenadas planas no necesariamente es el más cercano en la superficie elipsoidal) y aún más cuando el SRC de la capa de destino es geográfico (el filtro espacial asume las coordenadas geográficas como si fueran planas).

Por lo tanto, si tu intención es encontrar el punto más cercano en la superficie elipsoidal, asegúrate de reproyectar las capas a un SRC que no distorsione demasiado las distancias en tu área de trabajo, permitir que el filtro calcule de manera más amplia las distancias a varios de los puntos más cercanos, y finalmente elegir el más cercano en base a los valores resultantes en el campo de Distancia.