根據圓周率和地球半徑係數以及搜尋點的經緯度，搜尋資料表中與搜尋點之間的距離為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);

mysql> select * from `location`;
+----+--------------+----------------+---------------+
| id | name         | longitude      | latitude      |
+----+--------------+----------------+---------------+
|  1 | 廣州東站      | 113.3322640000 | 23.1562060000 |
|  2 | 林和西        | 113.3306110000 | 23.1472340000 |
|  3 | 天平架        | 113.3280950000 | 23.1653760000 |
+----+--------------+----------------+---------------+
 

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;
 

執行查詢：

mysql> 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;
+----+-----------+----------------+---------------+
| id | name      | longitude      | latitude      |
+----+-----------+----------------+---------------+
|  2 | 林和西     | 113.3306110000 | 23.1472340000 |
+----+-----------+----------------+---------------+