mysql 查询一个地点(经纬度) 附近N公里内的数据。(根据一个地点的经纬度查询这个地点方圆几公里内的数据)

1.创建测试表

CREATE TABLE `location` (

`id` int(10) unsigned NOT NULL AUTO_INCREMENT,

`name` varchar(50) NOT NULL,

`longitude` decimal(13,10) NOT NULL,

`latitude` decimal(13,10) NOT NULL,

PRIMARY KEY (`id`),

KEY `long_lat_index` (`longitude`,`latitude`)

) ENGINE=InnoDB DEFAULT CHARSET=utf8;

2.插入测试数据

insert into location(name,longitude,latitude) values

('广州东站',113.332264,23.156206),

('林和西',113.330611,23.147234),

('天平架',113.328095,23.165376);

1dcab1ec12aceb97a120a3a76a92dd54.png

3.搜寻1公里内的数据

搜寻点坐标:时代广场 113.323568, 23.146436

6370.996公里为地球的半径

计算球面两点坐标距离公式:

C = sin(MLatA)sin(MLatB)cos(MLonA-MLonB) + cos(MLatA)cos(MLatB)

Distance = RArccos(C)*Pi180

根据计算公式得到查询语句如下:

select * from `location` where (

acos(

sin(([#latitude#]*3.1415)/180) * sin((latitude*3.1415)/180) +

cos(([#latitude#]*3.1415)/180) * cos((latitude*3.1415)/180) * cos(([#longitude#]*3.1415)/180 - (longitude*3.1415)/180)

)*6370.996

)<=1;

执行查询:

SELECT * FROM `location` WHERE (ACOS(SIN((23.146436*3.1415)/180) * SIN((latitude*3.1415)/180) + COS((23.146436*3.1415)/180) * COS((latitude*3.1415)/180) * COS((113.323568*3.1415)/180 - (longitude*3.1415)/180))*6370.996)<=1;

58540a049b1b312d6f54544394fd723d.png

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐