haversine

I'm trying to implement a geospacial search in a php application. Currently, I'm using the following query to find points within 10km of a given latitude and longitude: SELECT * FROM ( SELECT *, (6378.1 * ACOS( COS(RADIANS(48.856614)) * COS(RADIANS(latitude)) * COS(RADIANS(2.3522219) - RADIANS(longitude)) + SIN(RADIANS(48.856614)) * SIN(RADIANS(latitude))) ) AS distance FROM `destinations_info` WHERE latitude BETWEEN 48.856614 - (10 / 69) AND 48.856614 + (10 / 69) AND longitude BETWEEN 2.3522219 - (10 / (69 * COS(RADIANS(48.856614)))) AND 2.3522219 + (10 / (69 * COS(RADIANS(48.856614)))) ) d

I have two arrays with lat and long. I want to calculate distance between every pair of lat and long with every other pair of lat and long in the array. Here are my two arrays. lat_array array([ 0.33356456, 0.33355585, 0.33355585, 0.33401788, 0.33370132, 0.33370132, 0.33370132, 0.33371075, 0.33371075, 0.33370132, 0.33370132, 0.33370132, 0.33356488, 0.33356488, 0.33370132, 0.33370132, 0.33370132, 0.33401788, 0.33362632, 0.33362632, 0.33364007, 0.33370132, 0.33401788, 0.33401788, 0.33358399, 0.33358399, 0.33358399, 0.33370132, 0.33370132, 0.33362632, 0.33370132, 0.33370132, 0.33370132, 0

I am implementing Haversine formula in PHP as the follows $result=mysqli_query($mysqli,"SELECT *,( 6371 * acos( cos( radians({$lat}) ) * cos( radians( `latitude` ) ) * cos( radians( `longitude` ) -radians({$lon}) ) +sin( radians({$lat}) ) * sin( radians( `latitude` ) ) ) ) AS distance FROM `places` HAVING distance <= {$radius} ORDER BY distance ASC") or die(mysqli_error($mysqli)); And Inside the Haversine fetch loop, I have the query that iterates through the results of haversine to select records that matches the IDs returned by the haversine formula. The Query is as follows. while($row =

I have the following model. class Location(models.Model): name = models.CharField(max_length = 128, blank = True) address =models.CharField(max_length = 200, blank= True) latitude = models.DecimalField(max_digits=6, decimal_places=3) longitude = models.DecimalField(max_digits=6, decimal_places=3) def __unicode__(self): return self.name If my current latitude & longitude is: current_lat = 43.648 current_long = 79.404 I did some research and came across the Haversine Equation which calculates the distance between two location coordinates. Below is the equation I found: import math def distance

I have an input of 36,742 points which means if I wanted to calculate the lower triangle of a distance matrix (using the vincenty approximation) I would need to generate 36,742*36,741*0.5 = 1,349,974,563 distances. I want to keep the pair combinations which are within 50km of each other. My current set-up is as follows shops= [[id,lat,lon]...] def lower_triangle_mat(points): for i in range(len(shops)-1): for j in range(i+1,len(shops)): yield [shops[i],shops[j]] def return_stores_cutoff(points,cutoff_km=0): below_cut = [] counter = 0 for x in lower_triangle_mat(points): dist_km = vincenty(x[0]

I currently have a MySQL table that is structured as follows: id name lon lat ----- ----- ----------- ----------- 1 Mark -76.316528 40.036027 2 John -95.995102 41.25716 3 Paul -82.337036 29.645095 4 Dave -82.337036 29.645095 5 Chris -76.316528 40.036027 The way I currently query the DB to see if a user's location is within a certain mile radius of a given location is by doing this. function Haversine($lat_from, $lon_from, $lat_to, $lon_to) { $radius = 6371000; $delta_lat = deg2rad($lat_to-$lat_from); $delta_lon = deg2rad($lon_to-$lon_from); $a = sin($delta_lat/2) * sin($delta_lat/2) + cos