地理位置SQL查询未找到确切位置

Question

我已经测试了一段时间的地理位置查询，到目前为止，我还没有发现任何问题。

我试图搜索给定半径内的所有城市，通常是使用城市坐标搜索城市周围的城市，但是最近我尝试在城市周围搜索，发现该城市本身未归还。

我在数据库中摘录了以下城市：

city            latitude    longitude  
Saint-Mathieu   45.316708   -73.516253  
Saint-脡douard   45.233374   -73.516254  
Saint-Michel    45.233374   -73.566256  
Saint-R茅mi      45.266708   -73.616257

但是，当我在Saint-R茅米市内运行查询时，出现以下查询…

SELECT tblcity.city, tblcity.latitude, tblcity.longitude, 
truncate((degrees(acos( sin(radians(tblcity.latitude)) 
* sin(radians(45.266708)) 
+ cos(radians(tblcity.latitude)) 
* cos(radians(45.266708)) 
* cos(radians(tblcity.longitude - -73.616257) ) ) ) 
* 69.09*1.6),1) as distance 
FROM tblcity HAVING distance < 10 ORDER BY distance desc

我得到这些结果：

city            latitude    longitude     distance  
Saint-Mathieu   45.316708   -73.516253    9.5  
Saint-脡douard   45.233374   -73.516254    8.6  
Saint-Michel    45.233374   -73.566256    5.3

搜索中缺少圣鲁茅米镇。

所以我尝试了一个修改后的查询，希望得到更好的结果：

SELECT tblcity.city, tblcity.latitude, tblcity.longitude, 
truncate(( 6371 * acos( cos( radians( 45.266708 ) ) 
* cos( radians( tblcity.latitude ) ) 
* cos( radians( tblcity.longitude ) 
- radians( -73.616257 ) ) 
+ sin( radians( 45.266708 ) ) 
* sin( radians( tblcity.latitude ) ) ) ),1) AS distance 
FROM tblcity HAVING distance < 10 ORDER BY distance desc

但我得到了相同的结果…

但是，如果我通过将经纬度的最后一位或长整数更改为1来稍微修改Saint-R茅mi的坐标，则两个查询都将返回Saint-
R茅mi。另外，如果我将查询集中在上述任何其他城市上，则搜索到的城市也会返回到结果中。

任何人都可以弄清是什么原因导致我的上述查询无法显示所搜索的圣鲁茅米市？我在下面添加了表格示例（已删除多余的字段）。

我使用的是MySQL 5.0.45，请先感谢。

CREATE TABLE `tblcity` ( 
`IDCity` int(1) NOT NULL auto_increment, 
`City` varchar(155) NOT NULL default '', 
`Latitude` decimal(9,6) NOT NULL default '0.000000', 
`Longitude` decimal(9,6) NOT NULL default '0.000000', 
PRIMARY KEY (`IDCity`) 
) ENGINE=MyISAM AUTO_INCREMENT=52743 DEFAULT CHARSET=latin1 AUTO_INCREMENT=52743;

INSERT INTO `tblcity` (`city`, `latitude`, `longitude`) VALUES 
('Saint-Mathieu', 45.316708, -73.516253), 
('Saint-脡douard', 45.233374, -73.516254), 
('Saint-Michel', 45.233374, -73.566256), 
('Saint-R茅mi', 45.266708, -73.616257);

Answer 1

在您的第一个查询中，我相信您已经颠倒了减法中的经度。余弦的球形定律是：

d = acos(sin(lat1)*sin(lat2) + cos(lat1)*cos(lat2)*cos(long2鈭抣ong1))*R

如果将lat1替换为tblcity.latitude，则必须将long1替换为tblcity.longitude。我认为您不小心在查询中替换了long2。这个更好吗？

SELECT tblcity.city, tblcity.latitude, tblcity.longitude, 
truncate((degrees(acos( sin(radians(tblcity.latitude)) 
* sin(radians(45.266708)) 
+ cos(radians(tblcity.latitude)) 
* cos(radians(45.266708)) 
* cos(radians(-73.616257 - tblcity.longitude) ) ) ) 
* 69.09*1.6),1) as distance 
FROM tblcity HAVING distance < 10 ORDER BY distance desc

我尚未调查您的第二个查询，但希望能对您有所帮助。