acos(1)返回某些值的NaN,而不是其他值

Question

我有一个纬度和经度值列表,我试图找到它们之间的距离.使用标准的大圆方法,我需要找到:

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舍入的一个奇怪的错误.

Answer 1

如果没有看到您的数据(尝试dput),就无法确切地说出来,但这很可能是FAQ 7.31的结果.

(x1 <- 1)
## [1] 1
(x2 <- 1+1e-16)
## [1] 1
(x3 <- 1+1e-8)
## [1] 1
acos(x1)
## [1] 0
acos(x2)
## [1] 0
acos(x3)
## [1] NaN

也就是说,即使您的值非常相似以至于它们的打印表示相同,它们仍然可能会有所不同:一些将在内.Machine$double.eps,另一些则不会......

确保输入值受[-1,1]限制的一种方法是使用pmax和pmin:acos(pmin(pmax(x,-1.0),1.0))