根据经度和纬度计算两点之间的距离
//
// 测地导航公式
//
function geoDistance(lat1, lng1, lat2, lng2, precision) {
// 自变量 precision 精度是小数点后的位数(距离精度)
var distance = 0;
if ((Math.abs(lat1 - lat2) < 0.00001) && (Math.abs(lng1 - lng2) < 0.00001))
{
distance = 0;
}
else {
lat1 = lat1 * Math.PI / 180;
lng1 = lng1 * Math.PI / 180;
lat2 = lat2 * Math.PI / 180;
lng2 = lng2 * Math.PI / 180;
var A = 6378140;
var B = 6356755;
var F = (A - B) / A;
var P1 = Math.atan((B / A) * Math.tan(lat1));
var P2 = Math.atan((B / A) * Math.tan(lat2));
var X = Math.acos(Math.sin(P1) * Math.sin(P2) + Math.cos(P1) * Math.cos(P2) * Math.cos(lng1 - lng2));
var L = (F / 8) * ((Math.sin(X) - X) * Math.pow((Math.sin(P1) + Math.sin(P2)), 2) / Math.pow(Math.cos(X / 2), 2) - (Math.sin(X) - X) * Math.pow(Math.sin(P1) - Math.sin(P2), 2) / Math.pow(Math.sin(X), 2));
distance = A * (X + L);
var decimal_no = Math.pow(10, precision);
distance = Math.round(decimal_no * distance / 1) / decimal_no;
// km转换为时(除以1000)
}
return distance;
}
从纬度/经度计算方向(方向)
function geoDirection(lat1, lng1, lat2, lng2) {
// 从纬度和经度lat1,lng1开始,一直到纬度和经度lat2,lng2
// 北方0度,顺时针角度0-360度
var Y = Math.cos(lng2 * Math.PI / 180) * Math.sin(lat2 * Math.PI / 180 - lat1 * Math.PI / 180);
var X = Math.cos(lng1 * Math.PI / 180) * Math.sin(lng2 * Math.PI / 180) - Math.sin(lng1 * Math.PI / 180) * Math.cos(lng2 * Math.PI / 180) * Math.cos(lat2 * Math.PI / 180 - lat1 * Math.PI / 180);
var dirE0 = 180 * Math.atan2(Y, X) / Math.PI; // 东0度
if (dirE0 < 0) {
dirE0 = dirE0 + 360; //0~360
}
var dirN0 = (dirE0 + 90) % 360; //(dirE0+90)÷360输出北向的其余部分为0度
return dirN0;
}