java中两个地理位置之间的距离：NaN

Question

我有http://www.codereye.com/2010/10/how-to-calculate-distance-between-2.html的这段代码，它根据两个IP地址对应的纬度和经度计算它们之间的大致物理距离。

public static double distance(double lat1, double lon1, double lat2, double  lon2) 
{
double theta = lon1 - lon2;
double dist = Math.sin(deg2rad(lat1)) * Math.sin(deg2rad(lat2)) + 
Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2)) * Math.cos(deg2rad(theta));

dist = Math.acos(dist);
dist = rad2deg(dist);
dist = dist * 60 * 1.1515;
return (dist);
}

lat和long是有符号值。但是，此代码有时会将dist返回为NaN。

我错过什么了吗？θ需要是绝对值吗？

更新：返回南的测试值：-

lat1=-23.5477，lon1=-46.6358，lat2=-23.5477，lon2=-46.6358

lat1=53.3331，lon1=-6.2489，lat2=53.3331，lon2=-6.2489

Answer 1

似乎当两个输入点相同时，以下表达式返回略高于1 (1.0000000000000002)的结果。

double dist = Math.sin(deg2rad(lat1)) * Math.sin(deg2rad(lat2)) + 
              Math.cos(deg2rad(lat1)) * Math.cos(deg2rad(lat2)) * Math.cos(deg2rad(theta));

这必须是浮点运算精度有限（或三角函数精度有限）的结果。

数学。acos（dist）如果dist返回NaN

dist = Math.acos(dist);

到

dist = Math.acos(Math.min(dist,1));