# Nearest Neighbor Search by Distance in Large Datasets with MySQL

Nearest neighbor search is common in geographical databases. Here we consider searching for nearby locations in a large geographical database. Many commercial databases have geographical indexes supporting such queries; however, in MySQL there is no simple implementation. We show a simple way below.

Consider the following table definition.

CREATE TABLE locations (
id INT NOT NULL AUTO_INCREMENT PRIMARY KEY,
lat DOUBLE NOT NULL,
lng DOUBLE NOT NULL,
address VARCHAR(64) NOT NULL,
KEY lat (lat),
KEY lng (lng)
);


Suppose we have a center at ($lat0,$lng0), and we wish to find nearby locations by distance. The following query contains a formula calculating the distance from the center. It returns the closet 20 locations within 10 km from the center. Here 6371 is the earth’s radius in km.

SELECT *,
( 6371 * acos(
cos(radians($lat0)) * cos(radians(lat)) * cos(radians(lng) - radians($lng0)) +
sin(radians($lat0)) * sin(radians(lat)) ) ) AS distance FROM locations HAVING distance < 10 ORDER BY distance LIMIT 20;  The query works well on small datasets; however, it is not efficient when the table is large (say, with 1 million records). This query scans all rows in the table, calculates the distance to the center, and finally sorts the result by distance. It runs in linear time of the number of rows. None of the indexes is effective. A quick fix is to specify a bounding box for latitude and longitude. The following query specifies a bounding box by simply adding/subtracting 1 degree from latitude/longitude of the center. SELECT * FROM locations WHERE lat <$lat0 + 1 AND lat > $lat0 - 1 AND lng <$lng0 + 1 AND lng > $lng0 - 1;  And the next query sorts the result by the distance to the center. SELECT *, ( 6371 * acos( cos(radians($lat0)) * cos(radians(lat)) * cos(radians(lng) - radians($lng0)) + sin(radians($lat0)) * sin(radians(lat))
) ) AS distance
FROM locations
WHERE lat < $lat0 + 1 AND lat >$lat0 - 1 AND
lng < $lng0 + 1 AND lng >$lng0 - 1
HAVING distance < 10
ORDER BY distance LIMIT 20;


However, this simple “degree manipulation” may not be good for applications requiring precise results. A better (but more complicated) way is to calculate the extremal degrees of points within certain distance from the center (the north/south/east/west point of 10 km to the center). Suppose we have parameters ($lat0,$lng0, $dist) for the center and the query distance. SELECT * FROM locations WHERE lat < degrees( asin( sin(radians($lat0)) * cos($dist / 6371) + cos(radians($lat0)) * sin($dist / 6371) * cos(radians(0)) )) AND lat > degrees( asin( sin(radians($lat0)) * cos($dist / 6371) + cos(radians($lat0)) * sin($dist / 6371) * cos(radians(180)) )) AND lng <$lng0 - degrees( atan2(sin(radians(90)) * sin(radians($dist / 6371)) * cos(radians($lat0)),                       cos(radians($dist / 6371)) - sin(radians($lat0)) * sin(radians($lat0))) ) AND lng >$lng0 + degrees( atan2(sin(radians(90)) * sin(radians($dist / 6371)) * cos(radians($lat0)),
cos(radians($dist / 6371)) - sin(radians($lat0)) * sin(radians($lat0))) )  The queries with bounding boxes will benefit from the indexes on “lat” and “lng”; the queries will be run with index scan on “lat” and “lng”. Note again that, we assume the dataset is large and scattered while the querying box is relatively small. The final query below searches the nearest 20 locations within distance$dist (in kilometers) from the center ($lat0,$lng0), with the result sorted by the distance.

SELECT *,
( 6371 * acos(
cos(radians($lat0)) * cos(radians(lat)) * cos(radians(lng) - radians($lng0)) +
sin(radians($lat0)) * sin(radians(lat)) ) ) AS distance FROM locations WHERE lat < degrees( asin( sin(radians($lat0)) * cos($dist / 6371) + cos(radians($lat0)) * sin($dist / 6371) * cos(radians(0)) )) AND lat > degrees( asin( sin(radians($lat0)) * cos($dist / 6371) + cos(radians($lat0)) * sin($dist / 6371) * cos(radians(180)) )) AND lng <$lng0 - degrees( atan2(sin(radians(90)) * sin(radians($dist / 6371)) * cos(radians($lat0)),
cos(radians($dist / 6371)) - sin(radians($lat0)) * sin(radians($lat0))) ) AND lng >$lng0 + degrees( atan2(sin(radians(90)) * sin(radians($dist / 6371)) * cos(radians($lat0)),
cos(radians($dist / 6371)) - sin(radians($lat0)) * sin(radians(\$lat0))) )
ORDER BY distance LIMIT 20;


One should note that the range of latitude is from -90 to 90 and the range of longitude is from -180 to 180; if the bounding box overlaps with the boundaries of the ranges, the WHERE condition has to be adapted!