我有一个纬度和经度值列表,我试图找出它们之间的距离。使用标准的大圆方法,我需要找到:
acos(sin(lat1)*sin(lat2) + cos(lat1)*cos(lat2) * cos(long2-long1))
并将其乘以地球的半径,以我使用的单位表示。只要我们取 acos 的值在 [-1,1] 范围内,这就是有效的。如果它们稍微超出这个范围,它会返回NaN
,即使差异是由于四舍五入造成的。
我遇到的问题是,有时,当两个纬度/经度值相同时,这会给我一个NaN
错误。并非总是如此,即使对于同一对数字,但总是列表中的相同数字。例如,我让一个人停在沙漠的路上:
Time |lat |long
1:00PM|35.08646|-117.5023
1:01PM|35.08646|-117.5023
1:02PM|35.08646|-117.5023
1:03PM|35.08646|-117.5023
1:04PM|35.08646|-117.5023
例如,当我计算连续点之间的距离时,第三个值将始终为NaN
,即使其他值不是。这似乎是 R 舍入的一个奇怪错误。