Finden Sie die nächsten Nachbarn zwischen zwei Tabellen mit Punktpositionen in SpatiaLite?

Ich habe heute angefangen, mit SpatiaLite zu spielen und bin bereits auf ein Problem gestoßen.

Für jede in tableOne gespeicherte Punktposition möchte ich einen nächstgelegenen (linearen Abstand) Punkt aus tableTwo auswählen.

Bisher habe ich eine ungeschickte Lösung gefunden, die VIEW verwendet:

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;

Und dann:

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

scheint den Job zu machen.

Zwei Fragen:

Gibt es eine Möglichkeit, eine solche Abfrage durchzuführen, ohne eine ANSICHT zu erstellen?

Gibt es eine andere Möglichkeit, diese Abfrage für eine bessere Leistung zu optimieren? In einem realen Szenario wird tableOne Hunderttausende von Datensätzen und tableTwo - 1,3 Millionen haben.

Ich habe gerade dieses SQL getestet und es funktioniert:

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

Wie Sie hier lesen können "Die naive Möglichkeit, eine Abfrage zum nächsten Nachbarn durchzuführen, besteht darin, die Kandidatentabelle nach Entfernung von der Abfragegeometrie zu ordnen und dann den Datensatz mit der geringsten Entfernung zu nehmen".

Freundliche Grüße,

Andrea

Wenn Sie die Abstände zwischen allen Punktkombinationen nicht berechnen möchten, können Sie einen räumlichen Index für eine der Tabellen verwenden:

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

Seit Version 4.4.0 unterstützt SpatiaLite einen virtuellen KNN-Tabellenindex für Probleme mit dem nächsten Nachbarn. Hier ist eine Abfrage, die die nächstgelegene Linie in einer Linestring-Tabelle zu jedem Punkt in einer Punktetabelle findet.

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

Sie können Ihre Abfrage so vereinfachen.

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

Für eine allgemeinere Lösung kann es sich lohnen, diese PostGIS Nearest Neighbor-Funktion zu konvertieren: http://blog.mackerron.com/2011/03/postgis-nearest-neighbour/