Пошук найближчих сусідів між двома таблицями з точковими місцями в SpatiaLite?

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

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

Поки що я придумав незграбне рішення, яке використовує ПОГЛЯД:

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

Як ви можете прочитати тут, "Наївний спосіб здійснити запит найближчого сусіда - це замовити таблицю кандидата на відстані від геометрії запиту, а потім взяти запис з найменшою відстані".

З повагою,

Андреа

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

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-nevable-neighbour/