这是我基本上想要实现的图片:
正如标题所说的那样,我想合并它们是半径(例如25Kms)的长/纬度点在长/纬点的边界框内触摸。
这是我非常简单的数据库结构:
+-------+-------------+------------+
| id | long | lat |
+-------+-------------+------------+
| 1 | -90.27137 | 50.00702 |
| 2 | -92.27137 | 52.00702 |
| 3 | -87.27137 | 48.00702 |
| 4 | -91.27137 | 51.00702 |
+-------+-------------+------------+
到目前为止,这是我的查询:
set @bottom_lat = -100.27137;
set @bottom_lon = 40.00702;
set @top_lat = -80.27137 ;
set @top_lon = 60.00702 ;
;
SELECT AVG(latitude), AVG(longitude)
FROM destination
WHERE latitude > @bottom_lat AND longitude > @bottom_lon AND latitude < @top_lat AND longitude < @top_lon
所以我的查询只是在假想的边界框内合并所有点而不考虑半径。
我知道我肯定会使用Haversine公式,但是我对数学和MySQL很苛刻,这让事情变得有点困难。事实上,如果我只有一个半径但是每个点都有自己的半径并且我正在努力,我最终可以合并点。
这是针对学生项目的,任何帮助都会得到很大的帮助。
参考文献:
- 我对SQL小提琴的查询:http://sqlfiddle.com/#!2/3a42b/2 (在评论中包含一个关于Haversine公式的SQL Fiddle例子)
- MySQL查询中的Haversine公式:(用于检查给定半径内的所有点)
SELECT*, ( 6371* acos( cos( radians(
@my_lat) ) * cos( radians(
destination.latitude ) ) * cos( radians(
destination.longitude ) - radians(
@my_lon) ) + sin( radians(
@my_lat) ) * sin( radians(
destination.latitude ) ) ) ) AS distance
FROM destination
ORDER BY distance limit 1
;
答案 0 :(得分:0)
如果没有PHP或其他编程语言的帮助,此操作可能太复杂而无法执行。以下是您在PHP中的表现:
<?
$link = mysqli_connect("host", "user", "pass", "database");
// Grab all the points from the db and push them into an array
$sql = "SELECT * FROM data";
$res = $link->query($sql);
$arr = array();
for($i = 0; $i < mysqli_num_rows($res); $i++){
array_push($arr, mysqli_fetch_assoc($res));
}
// Cycle through the point array, eliminating those points that "touch"
$rad = 1000; //radius in KM
for($i = 0; $i < count($arr); ++$i){
$lat1 = $arr[$i]['lat'];
$lon1 = $arr[$i]['long'];
for($j = 0; $j<count($arr); ++$j){
if($i != $j && isset($arr[$i]) && isset($arr[$j])){ // do not compare a point to itself
$lat2 = $arr[$j]['lat'];
$lon2 = $arr[$j]['long'];
// get the distance between each pair of points using the haversine formula
$dist = acos( sin($lat1*pi()/180)*sin($lat2*pi()/180) + cos($lat1*pi()/180)*cos($lat2*pi()/180)*cos($lon2*PI()/180-$lon1*pi()/180) ) * 6371;
if($dist < $rad){
echo "Removing point id:".$arr[$i]['id']."<br>";
unset($arr[$i]);
}
}
}
}
//display results
echo "Remaining points:<br>";
foreach($arr as $val){
echo "id=".$val['id']."<br>";
}
?>
此代码对您提供的数据的输出是:
Removing point id:1
Removing point id:2
Remaining points:
id=3
id=4
请注意,这只是删除重叠点,它不会对位置进行任何平均。你可以轻松添加它。希望这会有所帮助。