根据经纬度坐标计算实际距离

 

 

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;
  }