Прежде всего, запустите это в режиме консоли. Это создаст новую функцию MySQL (см. Ниже для phpMyAdmin).
DELIMITER |
DROP FUNCTION IF EXISTS GPS_Distance_Meters|
CREATE FUNCTION GPS_Distance_Meters (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
BEGIN
DECLARE rlo1 DOUBLE;
DECLARE rla1 DOUBLE;
DECLARE rlo2 DOUBLE;
DECLARE rla2 DOUBLE;
DECLARE dlo DOUBLE;
DECLARE dla DOUBLE;
DECLARE a DOUBLE;
SET rlo1 = RADIANS(lng1);
SET rla1 = RADIANS(lat1);
SET rlo2 = RADIANS(lng2);
SET rla2 = RADIANS(lat2);
SET dlo = (rlo2 - rlo1) / 2;
SET dla = (rla2 - rla1) / 2;
SET a = SIN(dla) * SIN(dla) + COS(rla1) * COS(rla2) * SIN(dlo) * SIN(dlo);
RETURN (6378137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END|
DELIMITER ;
Или вы можете запустить это в phpMyAdmin, но не забудьте изменить разделитель на | вместо;
DROP FUNCTION IF EXISTS GPS_Distance_Meters|
CREATE FUNCTION GPS_Distance_Meters (lat1 DOUBLE, lng1 DOUBLE, lat2 DOUBLE, lng2 DOUBLE) RETURNS DOUBLE
BEGIN
DECLARE rlo1 DOUBLE;
DECLARE rla1 DOUBLE;
DECLARE rlo2 DOUBLE;
DECLARE rla2 DOUBLE;
DECLARE dlo DOUBLE;
DECLARE dla DOUBLE;
DECLARE a DOUBLE;
SET rlo1 = RADIANS(lng1);
SET rla1 = RADIANS(lat1);
SET rlo2 = RADIANS(lng2);
SET rla2 = RADIANS(lat2);
SET dlo = (rlo2 - rlo1) / 2;
SET dla = (rla2 - rla1) / 2;
SET a = SIN(dla) * SIN(dla) + COS(rla1) * COS(rla2) * SIN(dlo) * SIN(dlo);
RETURN (6378137 * 2 * ATAN2(SQRT(a), SQRT(1 - a)));
END|
Теперь для использования функции. Пример:
SELECT *
FROM restaurants
WHERE GPS_Distance_Meters('48.858205', '2.294359', lat, lng) < 1000 ORDER BY proximite ASC
LIMIT 10
Здесь мы ищем limit 10 результаты 48.858205, 2.294359 (по сравнению с полем lat и lng) в пределах 1000 метров (1KM).