紧急求助，python处理表格数据

西河沿的风

知道两个表格中的GPS坐标点，求出每个任务坐标点的半径1千米范围内有几个会员。

要求记录该任务的坐标，价格，半径1千米内会员数，完成情况等数据，保存到第三个表格

拜托了，各位亲们，谢谢，谢谢，谢谢。





附件.zip (194.48 KB, 下载次数: 1)

2017-9-16 15:07 上传

点击文件名下载附件

西河沿的风

根据GPS坐标求距离的公式是
double rad(double d)

{

  const double PI = 3.1415926535898;

  return d * PI / 180.0;

}

// 从两个gps坐标点（经纬度）获得两点的直线距离，单位是米

int CalcDistance(float fLati1, float fLong1, float fLati2, float fLong2)

{

  const float EARTH_RADIUS = 6378.137;

  double radLat1 = rad(fLati1);

  double radLat2 = rad(fLati2);

  double a = radLat1 - radLat2;

  double b = rad(fLong1) - rad(fLong2);

  double s = 2 * asin(sqrt(pow(sin(a/2),2) + cos(radLat1)*cos(radLat2)*pow(sin(b/2),2)));

  s = s * EARTH_RADIUS;

  s = (int)(s * 10000000) / 10000;

  return s;

}

}

西河沿的风

from math import radians, cos, sin, asin, sqrt  
  
def haversine(lon1, lat1, lon2, lat2): # 经度1，纬度1，经度2，纬度2 （十进制度数）  
    """
    Calculate the great circle distance between two points  
    on the earth (specified in decimal degrees)
    """  
    # 将十进制度数转化为弧度  
    lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])  
  
    # haversine公式  
    dlon = lon2 - lon1   
    dlat = lat2 - lat1   
    a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2  
    c = 2 * asin(sqrt(a))   
    r = 6371 # 地球平均半径，单位为公里  
    return c * r * 1000  

西河沿的风

大佬们帮帮忙啊，感激不尽