两点经纬度计算距离

东方忘忧

已于 2023-04-20 13:08:43 修改

阅读量133

点赞数

分类专栏： c/c++ 文章标签：算法

于 2023-04-20 13:08:26 首次发布

本文链接：https://blog.csdn.net/weixin_43246170/article/details/130264117

版权

c/c++ 专栏收录该内容

86 篇文章 1 订阅

订阅专栏

#include "math.h"
#include "assert.h"

namespace {
    #define _TRUE   1
    #define _FALSE  0
    static const int Earth_Radius = 6378137;	//地球半径
    static const double e = 0.081813369f;

    double radian(double d)
    {
        return d*M_PI/180.0;//角度1°=M_PI/180
    }
}

double LL_dst::Get_distance(double lat1,double lng1,double lat2,double lng2)
{
    double radLat1=radian(lat1);
    double radLat2=radian(lat2);
    double a=radLat1-radLat2;
    double b=radian(lng1)-radian(lng2);
    double dst=2*asin((sqrt(pow(sin(a/2),2)+cos(radLat1)*cos(radLat2)*pow(sin(b/2),2))));
    dst=dst*Earth_Radius;
    //	dst=round(dst*10000)/10000;
    return dst;
}