MySQL按坐标排序和查找指定范围的坐标

做的项目多了，总会遇到一些需要涉及需要用到地图坐标的。
既然有坐标，那肯定又得涉及位置距离。
例如我们平时使用美团，想搜索附近的美食店，这功能很方便是不是？
希望下面的经验能给你帮助，让你少走些弯路，同时自己做个笔记，方便自己。
毕竟老了，有些东西很快就自己也忘记了，也需要笔记来回忆~

数据库表结构

CREATE TABLE `markers` (
  `id` INT NOT NULL AUTO_INCREMENT PRIMARY KEY COMMENT '主键ID' ,
  `name` VARCHAR( 60 ) NOT NULL COMMENT '标记名称' ,
  `address` VARCHAR( 80 ) NOT NULL COMMENT '标记地址' ,
  `lat` FLOAT( 10, 6 ) NOT NULL COMMENT '纬度' ,
  `lng` FLOAT( 10, 6 ) NOT NULL COMMENT '经度'
) ENGINE = InnoDB COMMENT = '标记表' ;

数据库表数据

INSERT INTO `markers` (`name`, `address`, `lat`, `lng`) VALUES ('北京市天安门','北京市东城区东长安街','39.915599','116.402687');
INSERT INTO `markers` (`name`, `address`, `lat`, `lng`) VALUES ('广州华立科技职业学院','广东省广州市增城广州华立科技园华立路7号','23.248335','113.871302');
INSERT INTO `markers` (`name`, `address`, `lat`, `lng`) VALUES ('韶关市风采楼','广东省韶关市浈江区风采路34号','24.813028','113.606039');

SQL语句

以下是基于半正矢公式 a = sin²(Δφ/2) + cos φ1 ⋅ cos φ2 ⋅ sin²(Δλ/2)的SQL查询语句。

6371是地球的半径，单位：公里。如果想以英里搜索，将6371换成3959即可。
39.915599是搜索点中心纬度（例如想搜索北京天安门附近的标记点，则这里就是北京天安门的纬度）
116.402687是搜索点中心经度（例如想搜索北京天安门附近的标记点，则这里就是北京天安门的经度）
distance字段是标记点与搜索点中心的距离，单位：公里（如果地球半径是英里，则这里也是英里）
25是范围，表示搜索出搜索中心点25公里以内的标记点

SELECT `id` , `name`
    , 6371 * acos(cos(radians(39.915599)) * cos(radians(`lat`)) * cos(radians(`lng`) - radians(116.402687)) + sin(radians(39.915599)) * sin(radians(`lat`))) AS `distance`
FROM `markers`
HAVING `distance` < 25
ORDER BY `distance`
LIMIT 0, 20;

查询结果

从以下结果能看出，还是有些误差的，两个点坐标一样，距离应该是0的。但误差较小，问题不大吧，嘿嘿。

id	name	distance
1	北京市天安门	0.00009493529796600342

来自：https://developers.google.com/maps/articles/phpsqlsearch_v3#findnearsql

已有 4 条评论

水木易安

2020-03-25 10:47:48

这个也不错 https://www.cnblogs.com/longhao/p/4071275.html

回复
1. jwj
  
  2020-03-30 11:11:46
  
  用MySQL自带的空间索引的确会好很多
  
  回复
鄙人李晋虎

2019-01-02 14:53:40

你好，零零博客的主人，我叫李晋虎。是这样，今天看了你这个算法感觉不错，但是关于这个计算公式这一块我想与你探讨一番，鄙人有两个问题想请教一下，鄙人也是实数算尽数据，通过你的算法也没有得到想要的数据，不知能否赐教一二，鄙人感激不尽。我的微信号就是Email的那个手机号，忘能百忙之中回复。感谢。

回复
1. jwj
  
  2019-01-18 13:47:42
  
  高看我了，以上算法来自谷歌地图里的文档。我只是将那部分摘录下来总结下，算法只是照搬而已。
  
  回复