Encontrando vizinhos mais próximos entre duas tabelas com localizações de pontos no SpatiaLite?

Comecei a jogar com o SpatiaLite hoje e já me deparei com um problema.

Para cada local do ponto armazenado no tableOne, eu gostaria de selecionar um ponto mais próximo (distância linear) da tableTwo.

Até agora, eu vim com uma solução desajeitada que utiliza o 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;

E depois:

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

parece fazer o trabalho.

Duas questões:

Existe uma maneira de executar essa consulta sem criar uma VIEW?

Existe alguma outra maneira de otimizar essa consulta para obter melhor desempenho? Em um cenário do mundo real, o tableOne terá centenas de milhares de registros e o tableTwo - 1,3 milhão.

Acabei de testar este SQL e 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 você pode ler aqui "A maneira ingênua de realizar uma consulta vizinha mais próxima é ordenar a tabela candidata à distância da geometria da consulta e, em seguida, registrar o registro com a menor distância".

Cumprimentos,

Andrea

Se você não deseja calcular as distâncias entre todas as combinações de pontos, use um índice espacial em uma das tabelas:

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 a versão 4.4.0, o SpatiaLite suporta um índice de tabela virtual KNN para problemas de vizinhos mais próximos. Aqui está uma consulta que localiza a linha mais próxima em uma tabela de cadeia de linhas para cada ponto em uma tabela de pontos.

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

Você pode simplificar sua consulta assim.

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 uma solução mais genérica, talvez valha a pena tentar converter essa função de vizinho mais próximo do PostGIS: http://blog.mackerron.com/2011/03/postgis-nearest-neighbour/