给定经纬度计算距离_通过经纬度坐标计算距离的方法（经纬度距离计算）ZZ

通过经纬度坐标计算距离的方法(经纬度距离计算)

最近在网上搜索“通过经纬度坐标计算距离的方法”，发现网上大部分都是如下的代码：

#define PI 3.14159265

static double Rc = 6378137;  // 赤道半径

static double Rj = 6356725;  // 极半径

class JWD

{

public:

double m_Longitude, m_Latitude;

double m_RadLo, m_RadLa;

double Ec;

double Ed;

public:

JWD(double longitude, double latitude)

{

m_Longitude = longitude;

m_Latitude = latitude;

m_RadLo = longitude * PI/180.;

m_RadLa = latitude * PI/180.;

Ec = Rj + (Rc - Rj) * (90.-m_Latitude) / 90.;

Ed = Ec * cos(m_RadLa);

}

~JWD() {};

};

static JWD GetJWDB(JWD A, double x,double y)

{

double dx=x;

double dy=y;

double BJD = (dx/A.Ed + A.m_RadLo) * 180./PI;

double BWD = (dy/A.Ec + A.m_RadLa) * 180./PI;

JWD B(BJD, BWD);

return B;

}

void main()

{

double referla=30.0;

double referlo=60.0;

double dx=500.0;

double dy=60.0;

JWD A(referla,referlo),B(0.0,0.0);

B=GetJWDB(A,dx,dy);

cout < < "  LA = " < < B.m_Latitude < < "  LO= " < < B.m_Longitude < < endl;

}

上面这段与之类似的代码是最容易通过搜索引擎找到的。大部分都是相互的转载，从来没有说明和注释。给初学者造成了很大的困扰。特别是：

Ec = Rj + (Rc - Rj) * (90.-m_Latitude) / 90.;

Ed = Ec * cos(m_RadLa);

Ec、Ed这2个参数，有人还在CSDN上发帖询问“函数中Ec，Ed是什么意思”。但从来没有见到有人回答。

提问的链接：

我刚开始接触这个问题，在中文世界中也只能搜到这些Ctrl+C 到Ctrl+V的东西。全球最大的中文互联网上也没有任何解答。已经明白的人知道很简单，但初学者肯定很疑惑。更何况现在LBS应用已经非常多了。肯定有非常多的人已经找到了更好的解答方式。 但对于我来说，用最常用的关键词，最容易的找到了这些复制的答案，但疑点重重。本着好奇的心，经过一番了解，我把结果给大家分析一下。下面是C#的代码：

public const double Ea = 6378137; //赤道半径

public const double Eb = 6356725; //极半径

private static void GetJWDB(double GLAT, double GLON, double distance, double angle, out double BJD, out doubleBWD)

{

double dx = distance * 1000 * Math.Sin(angle * Math.PI / 180.0);

double dy = distance * 1000 * Math.Cos(angle * Math.PI / 180.0);

//double ec = 6356725 + 21412 * (90.0 - GLAT) / 90.0;

//21412 是赤道半径与极半径的差

double ec = Eb + (Ea-Eb) * (90.0 - GLAT) / 90.0;

double ed = ec * Math.Cos(GLAT * Math.PI / 180);

BJD = (dx / ed + GLON * Math.PI /