我定义了一个城市模型,用于保存城市的geoname_id和位置(如GeoPt).我想要实现两件事.
>我希望距离某个城市500公里半径范围内的所有城市.
>我想计算两个城市之间以km为单位的距离.
实现这一目标的最佳方式是什么,请记住,我有一个非常庞大的城市数据库,我不想在性能因素上牺牲很多.任何帮助或建议表示赞赏.
最佳答案
这很完美,但速度很慢:
计算距离的功能.传递给此函数的参数是位置的纬度和经度元组或Geopt():
def haversinedistance(location1,location2):
"""Method to calculate distance between two sets of Lat/Lon."""
lat1,lon1 = location1
lat2,lon2 = location2
earth = 6371 #Earth's Radius in Kms.
#Calculate distance based in haversine Formula
dlat = math.radians(lat2-lat1)
dlon = math.radians(lon2-lon1)
a = math.sin(dlat/2) * math.sin(dlat/2) + math.cos(math.radians(lat1)) * math.cos(math.radians(lat2)) * math.sin(dlon/2) * math.sin(dlon/2)
c = 2 * math.atan2(math.sqrt(a),math.sqrt(1-a))
d = earth * c
return d
用于计算半径内的周边城市的函数.这是City模型下存储所有城市的方法:
def get_closest_cities(self,kms):
cities = []
#Find surrounding Cities of a given city within a given radius
allcities = self.country.city_set
for city in allcities:
distance = haversinedistance((self.location.lat,self.location.lon),(city.location.lat,city.location.lon))
if not distance >= kms:
cities.append((city.name,int(distance)))
cities.remove(cities[0])
return cities
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。