Расстояние между 2 точками в Постгисе в сриде 4326 в метрах

Это, вероятно, простой вопрос, но я не очень хорош в PostGIS и не до конца понимаю все это.

В основном у меня есть стол (nodes) со столбцом ТОЧКА (point). Я создал индекс для этого столбца

create index nodes__points on nodes using gist (point)

Столбец был создан с

select addgeometrycolumn('nodes', 'points', 'POINT', 4326, 2)

Я использую srid 4326, потому что я добавляю данные в форме (широта, долгота). (то есть система координат, где позиция Дублина, Ирландия широта =53,353 лон =-6,264 (который я добавил с GeomFromText('POINT(-6.264 53.535)'))).

Для каждой точки я хочу найти все точки, которые находятся примерно в пределах 1-километровой рамки с центром в этой точке (поэтому selcet a.id, count(*) from nodes as a, nodes as b where SOME_DISTANCE_FUNCTION_HERE(a.point, b.point, 1000) group by a.id;, Это не должно быть точным, просто грубая фигура. 1 км bbox в порядке, 1 км круг в порядке. Это не должно быть точно 1 км, только этот порядок величины.

ST_Distance/ST_DWithin/так далее. все используют единицы измерения SRID, который для 4326/WGS64 равен градусам (поэтому 1 = 1 градус широты / долготы). Но я хочу использовать метры.

Я старался ST_distance_sphere а также st_dwithin которые могут использовать метры, но если я сделаю это, то explain показывает, что индекс не используется.

Как я могу получить примерно то, что я хочу, и использовать географический индекс?

ОБНОВЛЕНИЕ: это в PostgreSQL 9.1 и PostGIS 2.0 svn build.

2 ответа

Решение

С тех пор, как я написал это, я обнаружил "географический" в отличие от типа "геометрия" в PostGIS, который может делать именно то, что я хочу.

Вы можете использовать ST_Transform для использования счетчиков, также помните, что не все функции доступны с типами географии, но если вам действительно нужна скорость, используйте ST_DWithin, это самый быстрый способ. Вот примерное соотношение между градусами и градусами:

| places | degrees    | distance |
| ------ | ---------- | -------- |
| 0      | 1.0        | 111 km   |
| 1      | 0.1        | 11.1 km  |
| 2      | 0.01       | 1.11 km  |
| 3      | 0.001      | 111 m    |
| 4      | 0.0001     | 11.1 m   |
| 5      | 0.00001    | 1.11 m   |
| 6      | 0.000001   | 0.111 m  |
| 7      | 0.0000001  | 1.11 cm  |
| 8      | 0.00000001 | 1.11 mm  |

Старый вопрос, но вы просто должны бросить так: ST_Distance(p.geom::geography, u.geom::geography) где p и u - псевдоним для двух таблиц

Другие вопросы по тегам