ฉันมีสถานที่ 2 แห่งที่กำหนดโดย gps coords, lat / long เช่นที่ส่งคืนโดย google maps:
http://gmaps-samples.googlecode.com/svn/trunk/geocoder/singlegeocode.html
ฉันต้องคำนวณระยะทางระหว่างพวกเขา ฉันรู้ว่าฉันสามารถใช้ google API ได้ แต่ฉันจะประมวลผลข้อความค้นหาเป็นจำนวนมากดังนั้นฉันควรทำกับเซิร์ฟเวอร์ของฉันเอง
ฉันใช้เวลาสองสามชั่วโมงกับเอกสารติดตั้ง geodjango แล้ว แต่ฉันไม่สามารถหาตัวอย่างได้ ทุกอย่างในวรรณคดีมีความซับซ้อนกว่าที่ฉันต้องการ
คำตอบ:
ดูเหมือนว่าคำตอบจะอยู่ในกลุ่ม Google Groupsนี้:
from django.contrib.gis.geos import GEOSGeometry
pnt = GEOSGeometry('SRID=4326;POINT(40.396764 -3.68042)')
pnt2 = GEOSGeometry('SRID=4326;POINT( 48.835797 2.329102 )')
pnt.distance(pnt2) * 100
* 100อะไร?
Pointจะไม่ได้ใช้
ฉันคิดว่ามันใช้ดีกว่าpyproj:
geod = pyproj.Geod(ellps='WGS84')
angle1,angle2,distance = geod.inv(long1, lat1, long2, lat2)
ดูเพิ่มเติม: http://blog.tremily.us/posts/pyproj/
คุณสามารถใช้คะแนนได้เช่นกัน
from django.contrib.gis.geos import Point
p1 = Point(37.2676483,-6.9273579)
p2 = Point(37.2653293,-6.9249401)
distance = p1.distance(p2)
distance_in_km = distance * 100
คุณสามารถใช้รหัส Python ของ Sven Marnachเพื่อรับผลลัพธ์ที่คุณต้องการ ฉันได้เพิ่มบรรทัดของรหัสในการรับผลเมตร
รหัส:
from math import sin, cos, radians, degrees, acos
def calc_dist(lat_a, long_a, lat_b, long_b):
lat_a = radians(lat_a)
lat_b = radians(lat_b)
long_diff = radians(long_a - long_b)
distance = (sin(lat_a) * sin(lat_b) +
cos(lat_a) * cos(lat_b) * cos(long_diff))
resToMile = degrees(acos(distance)) * 69.09
resToMt = resToMile / 0.00062137119223733
return resToMt
หากคุณต้องการคำตอบที่ไม่ได้ใช้ไลบรารีหรือฟังก์ชัน geodjango มองหาคำถามและคำตอบด้วยแท็กระยะทาง พวกเขาให้สูตรที่จะทำงานกับภาษาหรือกรอบงานใด ๆ คำถามหนึ่งคือระยะห่างระหว่างพิกัด GPS
ฉันชอบวิธีการแก้ปัญหาที่ฉันเคยเห็นครั้งหนึ่งด้วยความช่วยเหลือของ django และ geopy อย่างไรก็ตามฉันเปลี่ยนรหัสเล็กน้อยเพื่อให้มีอิสระในการป้อนมากกว่าสองจุด
from django.contrib.gis.geos import Point
from geopy.distance import distance as geopy_distance
from itertools import tee, izip
def pairwise(iterable):
a, b= tee(iterable)
next(b, None)
return izip(a,b)
chicago = Point(41.50, 87.37)
san_francisco = Point(37.47, 122.26)
st_louis = Point(38.62, 90.19)
washington = Point(38.53, 77.02)
points = (washington, st_louis, chicago, san_francisco)
d = sum(geopy_distance(a,b).meters for (a,b) in pairwise(points))
distance_km = d/1000
distance_miles = distance_km*0.621371
print "Distance in kilometre: ",distance_km
print "Distance in miles: ",distance_miles