最近在作小程序的时候遇到要计算两点之间的距离,但小程序没有相关的方法,只好本身来了.小程序
地球虽然是个椭圆,但也极像个正圆(相差大概20千米),因此就按正圆来算了,半径是6371.393千米.函数
L = 弧度 * R = 圆心角度数 × π × R / 180spa
2abcosC=a^2+b^2-c^2
c^2 = a^2+b^2-2abcosCcode
根据上面的图算 blog
假设A坐标(JA,WA),B坐标(JB,WB).
即∠AOC = WA, ∠BOD = WB. 图片
根据三角形函数 it
AC = sin(WA) R, BD = sin(WB) R
因此 BE = BD - ACio
由于∠COD = (JB - JA), 根据余弦定律
CO = cos(WA) R, DO = cos(WB) R
因此 AE = CD = 根号[CO^2 + DO^2 - 2 CO DO * cos(∠COD)]function
根据勾股定律
AB^2 = AE^2 + BE^2class
由于∆AOB是等腰三角形,三线合一,GO为垂直平分线
AB^2 = 2R^2 - 2R^2 * cosC;
cosC = [2R^2 - AB^2] / 2R^2;
简化的
∠AOB弧度 = acos[sin(WA)sin(WB) + cos(WA)cos(WB)cos(JB-JA)]
L = ∠AOB弧度 * R
function _Radian(num) { return num * Math.PI / 180; } function CalculateDistance(lata, lnga, latb, lngb) { var earthR = 6371.393; var WA, WB; WA = _Radian(lata); WB = _Radian(latb); var lngMinus = Math.abs(lngb - lnga) > 180? 360 - Math.abs(lngb - lnga): Math.abs(lngb - lnga); var lngRadian = _Radian(lngMinus); var ANGLE = Math.sin(WA) * Math.sin(WB) + Math.cos(WA) * Math.cos(WB) * Math.cos(lngRadian); var L = Math.acos(ANGLE) * earthR; return L; }
用_Radian发角度转化为弧度. 由于取两点的最小值,因此但经度之间差大于180时,再用360 - 经度差,取最小经度差.