double
hypot(
double
x,
double
y)
...
{
return sqrt(x * x + y * y);
}
double distance( double wd1, double jd1, double wd2, double jd2) ... {// 根据经纬度坐标计算实际距离
double x, y, out;
double PI = 3.1415926535898;
double R = 6.371229 * 1e6;
x = (jd2 - jd1) * PI * R * cos( ( (wd1 + wd2) / 2) * PI / 180) / 180;
y = (wd2 - wd1) * PI * R / 180;
out = hypot(x, y);
return out;
}
return sqrt(x * x + y * y);
}
double distance( double wd1, double jd1, double wd2, double jd2) ... {// 根据经纬度坐标计算实际距离
double x, y, out;
double PI = 3.1415926535898;
double R = 6.371229 * 1e6;
x = (jd2 - jd1) * PI * R * cos( ( (wd1 + wd2) / 2) * PI / 180) / 180;
y = (wd2 - wd1) * PI * R / 180;
out = hypot(x, y);
return out;
}