¿Encontrar vecinos más cercanos entre dos tablas con ubicaciones de puntos en SpatiaLite?

Hoy comencé a jugar con SpatiaLite y ya me topé con un problema.

Para cada ubicación de punto almacenada en tableOne, me gustaría seleccionar un punto más cercano (distancia lineal) de tableTwo.

Hasta ahora se me ocurrió una solución torpe que utiliza VIEW:

CREATE VIEW testview AS 
SELECT 
A.id , 
B.myValue, 
Distance(A.Geometry, B.Geometry) AS distance
FROM tableOne AS A, tableTwo AS B
WHERE distance < 10000
ORDER BY A.Id, distance;

Y entonces:

SELECT * FROM testview
WHERE distance = (SELECT MIN(distance) FROM testview AS t WHERE t.id = testview.id)

Parece hacer el trabajo.

Dos preguntas:

¿Hay alguna manera de realizar dicha consulta sin crear una VISTA?

¿Hay alguna otra forma de optimizar esta consulta para un mejor rendimiento? En un escenario del mundo real, tableOne tendrá cientos de miles de registros, y tableTwo - 1.3 millones.

Acabo de probar este SQL y funciona:

SELECT g1.OGC_FID As id1, g2.OGC_FID As id2, MIN(ST_Distance(g1.GEOMETRY,g2.GEOMETRY)) AS DIST
FROM table_01 As g1, table_02 As g2   
WHERE g1.OGC_FID <> g2.OGC_FID
AND ST_Contains(ST_Expand(g1.geometry,50),g2.geometry)
GROUP BY id1
ORDER BY id1

Como puede leer aquí "La manera ingenua de llevar a cabo una consulta vecina más cercana es ordenar la tabla de candidatos por distancia desde la geometría de la consulta, y luego tomar el registro con la distancia más pequeña".

Atentamente,

Andrea

Si no desea calcular las distancias entre todas las combinaciones de puntos, puede usar un índice espacial en una de las tablas:

SELECT 
  A.id , 
  B.myValue, 
  MIN(Distance(A.Geometry, B.Geometry)) AS distance
FROM tableOne AS A, tableTwo AS B
WHERE A.ROWID IN (
  SELECT ROWID
  FROM SpatialIndex WHERE
    f_table_name = 'A' 
    AND search_frame = BuildCircleMbr(ST_X(B.Geometry), ST_Y(B.Geometry), 10000))
GROUP BY A.id, B.myValue

Desde la versión 4.4.0, SpatiaLite admite un índice de tabla virtual KNN para los problemas de vecinos más cercanos. Aquí hay una consulta que encuentra la línea más cercana en una tabla de cadenas lineales a cada punto en una tabla de puntos.

SELECT k.* FROM knn k, points p
WHERE f_table_name = 'linestrings' 
AND ref_geometry = p.geometry
AND max_items = 1;

Puede simplificar su consulta de esta manera.

SELECT 
   A.id , 
   B.myValue, 
   MIN(Distance(A.Geometry, B.Geometry)) AS distance
FROM tableOne AS A, tableTwo AS B
GROUP BY A.id, B.myValue

Para una solución más genérica, podría valer la pena intentar convertir esta función de vecino más cercano de PostGIS: http://blog.mackerron.com/2011/03/postgis-nearest-neighbour/