相当于PostGIS在完美球体上的ST_Project

Question

我有一个PostGIS地理点，我想在理想球体（不是球体）上旋转一个设定的距离。我看到ST_Project（）函数，但它似乎在旋转球体上的点。例如，如果我向北旋转（0,10）1度，结果略大于1度

select ST_AsText(ST_Project(ST_Point(0,10), 111194.68229846345*1,radians(0.0)));

   st_astext         
---------------------------
 POINT(0 11.0052750273178)
(1 row)

旋转（0,60）相同的距离给出的结果略小于1度。

POINT(0 60.9979713953998)

我也可以用球体计算点之间的距离

select ST_Distance(ST_Project(ST_Point(60,89.9)::geography, 111194.68229846345*1,radians(0.0)),ST_Point(60,89.9),'t');
    st_distance    
------------------
 111194.682298475
(1 row)

还有完美的球体

select ST_Distance(ST_Project(ST_Point(60,89.9)::geography, 111194.68229846345*1,radians(0.0)),ST_Point(60,89.9),'f');
   st_distance   
-----------------
 110698.09475313
(1 row)

所以这并不是因为我把地球的周长搞错了。

一般来说，我的方位角将是0,90,180和270之一，它必须在极点周围正常工作。

Answer 1

您只需要创建一个新的空间投影，其中球体的半主轴和半次轴是相同的。所以对于半径为6370986 m的球体，这里有一个新的SRID=123456：

INSERT INTO spatial_ref_sys(srid, auth_name, srtext, proj4text)
VALUES (123456, NULL, '', '+proj=longlat +a=6370986 +b=6370986 +ellps=sphere +no_defs ');

对于地理类型，WGS84椭球体的默认SRID=4326。要使用新的球体椭球，您需要将SRID设置为非默认值。

SELECT ST_AsText(ST_Project(ST_SetSRID(ST_Point(0,10), 123456), 111194.68229846345, 0.0));

  st_astext
-------------
 POINT(0 11)
(1 row)