Trouver les voisins les plus proches entre deux tables avec des emplacements de points dans SpatiaLite?

J'ai commencé à jouer avec SpatiaLite aujourd'hui et j'ai déjà rencontré un problème.

Pour chaque emplacement de point stocké dans tableOne, je voudrais sélectionner un point (distance linéaire) le plus proche de tableTwo.

Jusqu'à présent, j'ai trouvé une solution maladroite qui utilise 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;

Et alors:

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

semble faire le travail.

Deux questions:

Existe-t-il un moyen d'effectuer une telle requête sans créer de VUE?

Existe-t-il un autre moyen d'optimiser cette requête pour de meilleures performances? Dans un scénario réel, tableOne aura des enregistrements de plusieurs centaines de milliers de couples, et tableTwo - 1,3 million.

Je viens de tester ce SQL et ça marche:

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

Comme vous pouvez le lire ici "La façon naïve d'effectuer une requête de voisin le plus proche est d'ordonner la table candidate par distance de la géométrie de la requête, puis de prendre l'enregistrement avec la plus petite distance".

Meilleures salutations,

Andrea

Si vous ne souhaitez pas calculer les distances entre toutes les combinaisons de points, vous pouvez utiliser un index spatial sur l'une des tables:

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

Depuis la version 4.4.0, SpatiaLite prend en charge un index de table virtuelle KNN pour les problèmes de voisin le plus proche. Voici une requête qui recherche la ligne la plus proche dans une table de chaîne de lignes de chaque point d'une table de points.

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

Vous pouvez simplifier votre requête comme ceci.

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

Pour une solution plus générique, il pourrait être utile d'essayer de convertir cette fonction PostGIS Nearest Neighbour: http://blog.mackerron.com/2011/03/postgis-nearest-neighbour/