¿Consulta PostGIS para seleccionar características no conectadas al resto de la red de carreteras?

Tengo un conjunto de datos de carreteras en una ciudad almacenado en una base de datos habilitada para PostGIS. Me gustaría identificar y eliminar cualquier característica que no esté conectada al resto de la red de carreteras. Este comportamiento es común en los bordes de las redes, donde se recortaron los datos.

He encontrado otras respuestas en el SIG GIS que hacen esto, pero todas parecen usar ArcGIS, y me gustaría resolver el problema usando PostGIS, QGIS o un script de Python.

Mis conjuntos de datos pueden incluir más de un millón de funciones, por lo que el rendimiento es una preocupación.

Tuve un problema similar, así que desarrollé un complemento QGIS que ahora está disponible en el repositorio oficial de complementos QGIS y se puede encontrar en el menú QGIS: Complementos / Administrar e instalar complementos ... y buscar Islas desconectadas.

http://plugins.qgis.org/plugins/disconnected-islands/

Este complemento se ejecuta en una capa de polilínea, construyendo un gráfico de red de carretera (o ferrocarril, etc.) de enlaces conectados. Luego analiza subgrafías conectadas, unas que están conectadas entre sí, pero que no están conectadas a enlaces aislados o flotantes. Crea un atributo adicional que contiene el ID de grupo del subgrafo. Esto se puede usar para dar estilo a la capa con estilos categorizados, o Zoom a selección. Los enlaces desconectados pueden repararse o eliminarse.

El código fuente se puede bifurcar desde: https://github.com/AfriGIS-South-Africa/disconnected-islands

Mi conjunto de datos tiene alrededor de 2 millones de enlaces y se analizó en menos de 15 minutos, usando 55 GB de RAM.

Puede usar fácilmente PostGIS para seleccionar carreteras que no se crucen con ninguna otra carretera:

SELECT id, geom FROM roads a
WHERE NOT EXISTS 
     (SELECT 1 FROM roads b 
      WHERE a.id != b.id
      AND   ST_Intersects(a.geom, b.geom))

Puede usar este script QGIS Python para detectar líneas que no están conectadas a nada:

from qgis.utils import iface

layer = iface.mapCanvas().currentLayer() # Selected layer

featureList = list( layer.getFeatures() ) # Gets all features
allFeatures = { feature.id(): feature for feature in featureList }

# Spatial indexing
spatialIdx = QgsSpatialIndex()
map( spatialIdx.insertFeature, featureList )

resList = [] # ids of features not connected to anything

for f in featureList:

    # List of potentially connected features from spatial indexing
    ids = spatialIdx.intersects( f.geometry().boundingBox() )

    hasNeighbor = False

    for id in ids:
        ifeature = allFeatures[id]

        if ifeature.id() == f.id():
            continue

        # Checks if f is really intersecting with ifeature
        if f.geometry().intersects(ifeature.geometry()):
            hasNeighbor = True
            break # Exit current for loop

    if (not hasNeighbor) and (not f.id() in resList):
        resList.append( f.id() )

print resList

Tenga en cuenta que esto no funcionará en líneas de varias partes. No creo que se pueda hacer mucho más rápido ...

Puede eliminar las carreteras fáciles con el método de @ dbaston primero, luego use pgRouting para encontrar los casos más complicados, como cuando tiene una red de carreteras que no está conectada a otra red.

Elija un segmento de carretera que esté definitivamente dentro de la red principal, luego intente encontrar una ruta entre cada segmento a ese. Si no existe una ruta, elimine el segmento.