假设我们知道两点p1(-104.673178,39.861656),p2(-87.904842,41.978603)的经度/纬度,我们想要找到距离由p1和p2确定的大圆100Km的P3的坐标。还要求大圆中与p3最近的点是p2。
这是我找到p3的方法:
首先使用 finalBearing 函数在p2找到大圆的最终方位。
fbearing=finalBearing(p1,p2,sphere=TRUE)
fbearing = 139.4564
因为大圆到p3的最近点是p2,角度p3-p2-p1应该等于90度(是真的吗?)。然后我们可以使用 destPoint 来查找p3
p3=destPoint(p2,-(180-fbearing)-90,100*1000,r=6378137,sphere=TRUE)
P 3 =( - 105.555,39.27437)
然而,当我使用 dist2gc 来计算一个点和一个大圆之间的距离时
dist2gc(p1, p2, p3, r=R)
答案是-100106.4,而不是-100000。
我的方法是找错了p3吗?或者106米是否可以接受错误?
地圈包的链接: https://cran.r-project.org/web/packages/geosphere/geosphere.pdf