我正在尝试使用Haversine formula
在维基百科上给出的这个实现进行试验,但这个公式没有给出预期的结果。
public class Haversine {
public static final double R = 6372.8; // In kilometers
public static double haversine(double lat1, double lon1, double lat2, double lon2) {
double dLat = Math.toRadians(lat2 - lat1);
double dLon = Math.toRadians(lon2 - lon1);
lat1 = Math.toRadians(lat1);
lat2 = Math.toRadians(lat2);
double a = Math.sin(dLat / 2) * Math.sin(dLat / 2) + Math.sin(dLon / 2) * Math.sin(dLon / 2) * Math.cos(lat1) * Math.cos(lat2);
double c = 2 * Math.asin(Math.sqrt(a));
return R * c;
}
public static void main(String[] args) {
System.out.println(haversine(36.12, -86.67, 33.94, -118.40));
}
}
Input GPS latitude format : ddmm.mmmm
Input GPS longitude format : dddmm.mmmm
上述经纬度格式在需求文档中指定。
示例输入坐标如下:
lat1 = 3359.64868, lon1 = 8356.178
lat2 = 3359.649, lon2 = 8356.178
在将这些值传递给 Haversine 方法之前,我将这些值转换为度数格式。如果不需要此步骤,请纠正我。
我正在使用下面的公式将度数格式转换为十进制度数格式:
Decimal Degree = degree + (minute / 60)
所以新坐标变成
lat1 = 33 + (59.64868 / 60) = 33.994144666666664
lon1 = 83 + (56.178 / 60) = 83.9363
lat2 = 33 + (59.649 / 60) = 33.99415
lon2 = 83 + (56.178 / 60) = 83.9363
调用 hasrsine 方法变得像
haversine(33.994144666666664, 83.9363, 33.99415, 83.9363)
这是返回值5.932071604620887E-4
为了验证这些值,我向本网站(33.994144666666664, 83.9363, 33.99415, 83.9363)
上的转换器提供了相同的输入,但它给出的结果为km。0.001
我试图在不转换为十进制度的情况下提供输入值,但是两种方法的输出也不匹配。
谁能告诉我我在这里犯了什么错误?