html按经纬度查询距离,根据经纬度与距离计算另一个经纬度

/**

* 根据一个经纬度及距离角度,算出另外一个经纬度

* @param {*} lng 经度 113.3960698

* @param {*} lat 纬度 22.941386

* @param {*} brng 方位角 45 ---- 正北方:000°或360° 正东方:090° 正南方:180° 正西方:270°

* @param {*} dist 90000距离(米)

*

* 1、角度转换为弧度公式:弧度=角度×(π ÷度180 )

* 2、弧度转换为角度公式: 角度=弧度×(180÷π)

*/

mapHelper.prototype.getLonAndLat = function (lng, lat, brng, dist) {

//大地坐标系资料WGS-84 长半径a=6378137 短半径b=6356752.3142 扁率f=1/298.2572236

var a = 6378137;

var b = 6356752.3142;

var f = 1 / 298.257223563;

var lon1 = lng * 1;

var lat1 = lat * 1;

var s = dist;

var alpha1 = brng * (Math.PI / 180)//mapNumberUtil.rad(brng);

var sinAlpha1 = Math.sin(alpha1);

var cosAlpha1 = Math.cos(alpha1);

//var tanU1 = (1 - f) * Math.tan(mapNumberUtil.rad(lat1));

var tanU1 = (1 - f) * Math.tan(lat1 * (Math.PI / 180));

var cosU1 = 1 / Math.sqrt((1 + tanU1 * tanU1)), sinU1 = tanU1 * cosU1;

var sigma1 = Math.atan2(tanU1, cosAlpha1);

var sinAlpha = cosU1 * sinAlpha1;

var cosSqAlpha = 1 - sinAlpha * sinAlpha;

var uSq = cosSqAlpha * (a * a - b * b) / (b * b);

var A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));

var B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq)));

var sigma = s / (b * A), sigmaP = 2 * Math.PI;

while (Math.abs(sigma - sigmaP) > 1e-12) {

var cos2SigmaM = Math.cos(2 * sigma1 + sigma);

var sinSigma = Math.sin(sigma);

var cosSigma = Math.cos(sigma);

var deltaSigma = B * sinSigma * (cos2SigmaM + B / 4 * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) -

B / 6 * cos2SigmaM * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM)));

sigmaP = sigma;

sigma = s / (b * A) + deltaSigma;

}

var tmp = sinU1 * sinSigma - cosU1 * cosSigma * cosAlpha1;

var lat2 = Math.atan2(sinU1 * cosSigma + cosU1 * sinSigma * cosAlpha1,

(1 - f) * Math.sqrt(sinAlpha * sinAlpha + tmp * tmp));

var lambda = Math.atan2(sinSigma * sinAlpha1, cosU1 * cosSigma - sinU1 * sinSigma * cosAlpha1);

var C = f / 16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha));

var L = lambda - (1 - C) * f * sinAlpha *

(sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));

var revAz = Math.atan2(sinAlpha, -tmp); // final bearing

var lngLatObj = { lng: lon1 + L * (180 / Math.PI), lat: lat2 * (180 / Math.PI) }

return lngLatObj;

}

只是将弧度与角度直接写在了里面

原文:https://www.cnblogs.com/zhoushangwu/p/12922817.html

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值