Почему Postgres (PostGis) не верно вычисляет вхождения точек в радиус?

Question

[Данил Чекалин]

Существует таблица с таким DDL:

CREATE TABLE buildings
(
    id INTEGER PRIMARY KEY NOT NULL,
    city VARCHAR(255) NOT NULL,
    street VARCHAR(255) NOT NULL,
    build_number DOUBLE PRECISION NOT NULL,
    location GEOMETRY(POINT,4326) // пробовал и GEOGRAPHY(POINT,4326)
);

В таблице 100 записей, каждое поле location имеет значение:

sprintf('SRID=4326;POINT(%s %s)', $faker->longitude(55, 55.08), $faker->latitude(82.93, 83.01));

То есть точки сосредоточены в одном месте не далеко от друг друга.

Теперь отображение точека на карте:


Левая карта, это вывод всех точек. Права карта имеет вывод красной точки в радиусе которой будет происходить поиск.
Красная точка находится в координате

55.0502, 82.9615

Производим поиск с радиусом 500м. Запрос:

SELECT *
FROM "buildings"
WHERE ST_DWithin(location, ST_GeomFromText('POINT(55.0502 82.9615)'), 500, TRUE)

Результат:


Производим поиск с радиусом 1500м. Запрос:

SELECT *
FROM "buildings"
WHERE ST_DWithin(location, ST_GeomFromText('POINT(55.0502 82.9615)'), 1500, TRUE)

Результат:


Результат с обеими картами в меньшем маштабе.

Answer 1

[Rsa97]

Вы используете ST_DWithin с четырьмя параметрами, это Geography-режим. Для него объекты должны иметь одинаковый SRID. Попробуйте убрать четвёртый параметр или создавать точку с явным указанием SRID
ST_GeomFromText('POINT(55.0502 82.9615)', 4326)

Comments

[Данил Чекалин]

Удаление 4 агрумента и явное указание второго помогло

[Данил Чекалин]

Возможно подскажите, в 3й аргумент метры уже не пойду как единицы изменериня. Как конвертировать degrees в meters, чтоб получать корректную выборку?

[Rsa97]

Данил Чекалин: Напрямую никак. Если один градус широты это примерно 111 км, то градус долготы на экваторе те же 111 км, а на полюсе 0.

[Данил Чекалин]

Какое решение то этой проблемы?

[Rsa97]

Данил Чекалин: Попробуйте преобразование в Geography

ST_DWithin(location::geography , ST_GeomFromText('POINT(120.9982 24.788)',4326)::geography , 100)

[Данил Чекалин]

Это я изначально пробовал делать и как я уже писал в DDL, что GEOGRAPHY(POINT,4326) я пробовал. Результат такой же как на скриншотах

[Rsa97]

Данил Чекалин: Сейчас проверил, если колонка имеет тип GEOGRAPHY(POINT), тогда сравнение идёт в метрах. Запрос

ST_DWithin(location, ST_GeomFromText('POINT(120.9982 24.788)', 100)

dв этом случае вроде работает нормально.

[Данил Чекалин]

увы ни в какую не хочет нормально работать с geography

[Алексей Мартьянов]

Доброго времени суток! Столкнулся с такой же проблемой.. Подскажите, получилось победить сей недуг?

[Данил Чекалин]

Добрый день, не получилось, да и задача была тестовая и не требовала 100% решения. Это я больше для себя решал.