Как я могу расширить этот SQL-запрос, чтобы найти k ближайших соседей?

у меня есть база данных, полная двумерных точек данных на карте. Каждая запись имеет поле типа "геометрия". Что мне нужно сделать, так это передать точку в хранимую процедуру, которая возвращает k ближайшие точки (k также будет передан sproc, но это легко). Я нашел запрос вhttp://blogs.msdn.com/isaac/archive/2008/10/23/nearest-neighbors.aspx который получает единственного ближайшего соседа, но я не могу понять, как расширить его, чтобы найти k ближайших соседей.

это текущий запрос - T стол, g поле геометрии, @x это точка для поиска вокруг,Numbers таблица с числами от 1 до n:

DECLARE @start FLOAT = 1000; 
WITH NearestPoints AS
(
     SELECT TOP(1) WITH TIES *,  T.g.STDistance(@x) AS dist
     FROM Numbers JOIN T WITH(INDEX(spatial_index)) 
     ON T.g.STDistance(@x) < @start*POWER(2,Numbers.n)
     ORDER BY n
)
SELECT TOP(1) * FROM NearestPoints
ORDER BY n, dist

внутренний запрос выбирает ближайшую непустую область, а внешний запрос затем выбирает верхний результат из этой области; внешний запрос можно легко изменить на (например) SELECT TOP(20), но если ближайший регион содержит только один результат, ты застрял с этим.

Я полагаю, что мне, вероятно, нужно рекурсивно искать первую область, содержащую k записи, но без использования переменной таблицы (что вызовет проблемы с обслуживанием, поскольку вам нужно создать структуру таблицы, и она может измениться - есть много полей), я не вижу, как.

9
geospatial nearest-neighbor sql sql-server-2008

2 ответов

что произойдет, если убрать TOP (1) WITH TIES из внутреннего запроса и установите внешний запрос для возврата верхнего k строк?

мне также было бы интересно узнать, помогает ли эта поправка вообще. Он должен быть более эффективным, чем использование TOP:

DECLARE @start FLOAT = 1000
        ,@k INT = 20
        ,@p FLOAT = 2;

WITH NearestPoints AS
(
     SELECT *
            ,T.g.STDistance(@x) AS dist
            ,ROW_NUMBER() OVER (ORDER BY T.g.STDistance(@x)) AS rn
     FROM Numbers 
     JOIN T WITH(INDEX(spatial_index)) 
     ON   T.g.STDistance(@x) <  @start*POWER(@p,Numbers.n)
     AND (Numbers.n - 1 = 0 
          OR T.g.STDistance(@x) >= @start*POWER(@p,Numbers.n - 1)
         )
)
SELECT * 
FROM NearestPoints
WHERE rn <= @k;

NB-непроверенный-у меня нет доступа к SQL 2008 здесь.

2

цитата из внутри Microsoft® SQL Server® 2008: Программирование T-SQL. Раздел 14.8.4.

следующий запрос вернет 10 достопримечательности, ближайшие к @input:

DECLARE @input GEOGRAPHY = 'POINT (-147 61)';
DECLARE @start FLOAT = 1000;
WITH NearestNeighbor AS(
  SELECT TOP 10 WITH TIES
    *, b.GEOG.STDistance(@input) AS dist
  FROM Nums n JOIN GeoNames b WITH(INDEX(geog_hhhh_16_sidx)) -- index hint
  ON b.GEOG.STDistance(@input) < @start*POWER(CAST(2 AS FLOAT),n.n)
  AND b.GEOG.STDistance(@input) >=
    CASE WHEN n = 1 THEN 0 ELSE @start*POWER(CAST(2 AS FLOAT),n.n-1) END
  WHERE n <= 20
  ORDER BY n
)
  SELECT TOP 10 geonameid, name, feature_code, admin1_code, dist
  FROM NearestNeighbor
  ORDER BY n, dist;

Примечание: только часть этого запроса, где предложение поддерживается пространственным индекс. Однако оптимизатор запросов правильно оценивает поддерживаемую деталь (сравнение" = " необходимо проверить, и запрос выполняет хорошо. Изменение значение @start может иногда ускорить запрос, если он медленнее чем хотелось.

листинг 2-1. Создание и заполнение вспомогательной таблицы чисел

SET NOCOUNT ON;
USE InsideTSQL2008;

IF OBJECT_ID('dbo.Nums', 'U') IS NOT NULL DROP TABLE dbo.Nums;

CREATE TABLE dbo.Nums(n INT NOT NULL PRIMARY KEY);
DECLARE @max AS INT, @rc AS INT;
SET @max = 1000000;
SET @rc = 1;

INSERT INTO Nums VALUES(1);
WHILE @rc * 2 <= @max
BEGIN
  INSERT INTO dbo.Nums SELECT n + @rc FROM dbo.Nums;
  SET @rc = @rc * 2;
END

INSERT INTO dbo.Nums
  SELECT n + @rc FROM dbo.Nums WHERE n + @rc <= @max;
2