Точка отсчёта за 50 миль (север, 45% к северо-востоку, 45% к юго-востоку) - PullRequest
1 голос
/ 21 января 2020

В PostGIS, есть ли способ вычислить еще одну точку в 50 милях в разных направлениях?

По заданной точке («Нью-Йорк», - 74,00,40,71), как рассчитать следующие точки? ?

1) 50 miles directly North
2) 50 miles 45% North East
4) 50 miles directly East
3) 50 miles 45% South West

Обновление: кажется, http://postgis.net/docs/ST_Project.html может быть решением.

ST_Project('POINT(-74.00 40.71)'::geography, 80467.2, radians(45.0))

Однако мне нужно обратиться к записи базы данных, чтобы сделать это. не жесткий код это.

1 Ответ

1 голос
/ 21 января 2020

Попробуйте объединить ST_Project с CTE - отрегулируйте значения radians до нужного азимута.

WITH j AS (
  SELECT poi::geography AS poi FROM t
)
SELECT 
  ST_AsText(ST_Project(j.poi, 80467.2, radians(90.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(45.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(180.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(135.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(270.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(225.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(360.0)),2),
  ST_AsText(ST_Project(j.poi, 80467.2, radians(315.0)),2)
FROM j;

      st_astext      |      st_astext      |    st_astext     |     st_astext      |      st_astext      |     st_astext      |    st_astext     |      st_astext      
---------------------+---------------------+------------------+--------------------+---------------------+--------------------+------------------+---------------------
 POINT(-73.05 40.71) | POINT(-73.32 41.22) | POINT(-74 39.99) | POINT(-73.33 40.2) | POINT(-74.95 40.71) | POINT(-74.67 40.2) | POINT(-74 41.43) | POINT(-74.68 41.22)
(1 Zeile)

enter image description here

Примечание : буфер на изображении только для иллюстрации.

...