#stackoverflow.com/questions/4913349/haversine-formula-in-python-bearing-and-distance-between-two-gps-points def haversine(lon1, lat1, lon2, lat2): lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, 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)) km = 6367 * c #convert from km to miles miles = km * 0.621371 return round(miles) #Takes an initial latitude, longitude, and specified radius, and returns domain of latitude def find_latitude_range(lat1, lon1, radius): lon1, lat1, lon2 = map(radians, [lon1, lat1, lon1]) #latitude is changing, not longitude, so dlon = 0 dlon = 0 c = radius / 6367 / 0.621363 a = (sin(c/2))**2 #dlon is 0, so the right side evaluates to 0 #a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2 #a = sin(dlat/2) ** 2 #sqrt(a) = sin(dlat/2) #asin(sqrt(a)) = (lat2- lat1) / 2 #( 2 * asin(sqrt(a))) + lat1 = lat2
Run
Reset
Share
Import
Link
Embed
Language▼
English
中文
Python Fiddle
Python Cloud IDE
Follow @python_fiddle
Browser Version Not Supported
Due to Python Fiddle's reliance on advanced JavaScript techniques, older browsers might have problems running it correctly. Please download the latest version of your favourite browser.
Chrome 10+
Firefox 4+
Safari 5+
IE 10+
Let me try anyway!
url:
Go
Python Snippet
Stackoverflow Question
