Найти ближайших соседей между двумя таблицами с точечными местоположениями в SpatiaLite?

Я начал играть со SpatiaLite сегодня и уже наткнулся на проблему.

Для каждой точки, сохраненной в tableOne, я бы хотел выбрать одну, ближайшую (линейное расстояние) точку из tableTwo.

До сих пор я придумал неуклюжее решение, которое использует 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;

А потом:

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

кажется, делает работу.

Два вопроса:

Есть ли способ выполнить такой запрос без создания VIEW?

Есть ли другой способ оптимизировать этот запрос для повышения производительности? В реальном сценарии tableOne будет иметь сотни-пару тысяч записей, а tableTwo - 1,3 миллиона.

Я только что проверил этот SQL, и он работает:

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

Как вы можете прочитать здесь: «Наивный способ выполнить запрос ближайшего соседа - упорядочить таблицу кандидатов по расстоянию от геометрии запроса, а затем взять запись с наименьшим расстоянием».

С наилучшими пожеланиями,

Andrea

Если вы не хотите вычислять расстояния между всеми комбинациями точек, вы можете использовать пространственный индекс для одной из таблиц:

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

Начиная с версии 4.4.0 SpatiaLite поддерживает индекс виртуальной таблицы KNN для проблем ближайшего соседа. Вот запрос, который находит ближайшую строку в таблице линейных линий для каждой точки в таблице точек.

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

Вы можете упростить ваш запрос следующим образом.

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

Для более общего решения, возможно, стоит попытаться преобразовать эту функцию ближайшего соседа PostGIS: http://blog.mackerron.com/2011/03/postgis-nearest-neighbour/