【python学习笔记】Python对经纬度处理

说明

因为地球是球面，所以地球平面间的距离也得根据球面来计算。
连个问题：
- 两经纬度点之间的距离
- 根据一个给定经纬度的点，进行附近若干距离地点查询

两点之间的距离

# 根据经纬度计算距离def __distance(lon1, lat1, lon2, lat2):  # 经度1，纬度1，经度2，纬度2 （十进制度数）    """    根据经纬度计算距离    :param lon1: 点1经度    :param lat1: 点1纬度    :param lon2: 点2经度    :param lat2: 点2纬度    :return:distance    """    # 将十进制度数转化为弧度    lon1, lat1, lon2, lat2 = map(radians, [float(lon1), float(lat1), float(lon2), float(lat2)])    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.137  # 地球平均半径，单位为公里    return float('%.2f' % (c * r))

Mysql-根据一个给定经纬度的点，进行附近若干距离地点查询

需求：已知当前位置，查找附近10km内的所有站点
效果应该是地图上有一个圆形的范围，但是这种思路只想到对每条记录，去进行遍历，跟数据库中的每一个点进行距离计算，当距离小于500米时，认为匹配。这样做确实能够得到结果，但是效率极其低下，因为每条记录都要去循环匹配n条数据，其消耗的时间可想而知。

看到网上有个思路:利用我们想要得到圆的外接正方形，然后拿正方形的经纬度范围去搜索数据库。

# 确定查询经纬度范围def __get_area(latitude, longitude, dis):    """    确定查询经纬度范围    :param latitude:中心纬度    :param longitude:中心经度    :param dis:半径    :return:(minlat, maxlat, minlng, maxlng)    """    r = 6371.137    dlng = 2 * math.asin(math.sin(dis / (2 * r)) / math.cos(latitude * math.pi / 180))    dlng = dlng * 180 / math.pi    dlat = dis / r    dlat = dlat * 180 / math.pi    minlat = latitude - dlat    maxlat = latitude + dlat    minlng = longitude - dlng    maxlng = longitude + dlng    return minlat, maxlat, minlng, maxlng

leason | 博客