使用未投影SpatialPointsDataFrame时spDists函数的返回单位是什么

你使用spDists计算CRS为+proj=longlat +datum=WGS84 +no_defs的未投影空间点大圆距离时，返回结果的单位确实是千米（km），你的猜测正确。

相关补充说明：

• 当输入空间对象为经纬度坐标系时，spDists会自动将longlat参数设为TRUE，调用球面大圆距离算法计算，输出单位固定为千米；如果输入为平面投影坐标系，输出单位才会和投影单位保持一致。
• 你可以用已知坐标点做简单验证：例如北京（116.40°E, 39.90°N）到上海（121.47°E, 31.23°N）的实测大圆距离约为1068km，代入函数计算得到的结果和该数值基本吻合。
• 如果你需要以米为单位的结果，直接将计算得到的距离值乘以1000即可。

内容的提问来源于stack exchange，提问作者Leah Bevis