
数据库几十万万条数据中查询出符合“周围1公里范围内”条件的坐标点
Mysql本身是支持空间索引的,但是在5.X版本中取消了Distance()和Related(),无法使用空间的距离函数去直接查询距离在一定范围内的点。所以,我首先想到的是,对每条数据去进行遍历,跟数据库中的每个点进行距离计算,当距离小于1公里时候,认为匹配成功。经测试,这样做确实能得到结果,但是效率极其低下,因为每条数据都得去和数据库中的几十万条数据进行比对,其耗费的时间可想而知
后来经过自己的仔细考虑,以及查询各种资料,想到一种利用正方形将方圆1公里这个圆包围起来。利用正方形的四个点,和大厦经纬度进行比较。问题转向了,如何计算正方形四个点经纬度

假设已知点的经纬度分别为lng,lat
先实现经度范围的查询,
在haversin公式中令φ1 = φ2
/**
* 计算某个经纬度的周围某段距离的正方形的四个点
*
* @param
* radius 地球半径 平均6371km
* @param
* lng float 经度
* @param
* lat float 纬度
* @param
* distance float 该点所在圆的半径,该圆与此正方形内切,默认值为1千米
* @return array 正方形的四个点的经纬度坐标
*/
public function returnSquarePoint($lng, $lat, $distance = 1, $radius = 6371)
{
$dlng = 2 * asin(sin($distance / (2 * $radius)) / cos(deg2rad($lat)));
$dlng = rad2deg($dlng);
$dlat = $distance / $radius;
$dlat = rad2deg($dlat);
return array(
'left-top' => array(
'lat' => $lat + $dlat,
'lng' => $lng - $dlng
),
'right-top' => array(
'lat' => $lat + $dlat,
'lng' => $lng + $dlng
),
'left-bottom' => array(
'lat' => $lat - $dlat,
'lng' => $lng - $dlng
),
'right-bottom' => array(
'lat' => $lat - $dlat,
'lng' => $lng + $dlng
)
);
}
$array[0]就是起点终点坐标数组
$room = $this->returnSquarePoint($array[0]['lng'], $array[0]['lat']);
->andwhere(['>','start_lat',$start['right-bottom']['lat']])
->andWhere(['<','start_lat',$start['left-top']['lat']])
->andWhere(['>','start_lng',$start['left-top']['lng']])
->andWhere(['<','start_lng',$start['right-bottom']['lng']]);
2.-----第二种
$range = 180 / pi() * 1 / 6372.797; //里面的 1 就代表搜索 1km 之内,单位km
$lngR = $range / cos($myLat * pi() / 180);
$maxLat = $myLat + $range; //最大纬度
$minLat = $myLat - $range; //最小纬度
$maxLng = $myLng + $lngR; //最大经度
$minLng = $myLng - $lngR; //最小经度
3.*----第三种
mysql两点间距,
公式如下,单位米:
第一点经纬度:lng1 lat1
第二点经纬度:lng2 lat2
round(6378.138*2*asin(sqrt(pow(sin(
(lat1*pi()/180-lat2*pi()/180)/2),2)+cos(lat1*pi()/180)*cos(lat2*pi()/180)*
pow(sin( (lng1*pi()/180-lng2*pi()/180)/2),2)))*1000)
例如:
SELECT store_id,lng,lat,
ROUND(6378.138*2*ASIN(SQRT(POW(SIN((22.299439*PI()/180-lat*PI()/180)/2),2)+COS(22.299439*PI()/180)*COS(lat*PI()/180)*POW(SIN((114.173881*PI()/180-lng*PI()/180)/2),2)))*1000)
AS
juli
FROM store_info
ORDER BY juli DESC
LIMIT 316
4.
select ROUND(6378.138*2*ASIN(SQRT(POW(SIN(($latitude*PI()/180-latitude*PI()/180)/2),2)+COS($latitude*PI()/180)*COS(latitude*PI()/180)*POW(SIN(($longitude*PI()/180-longitude*PI()/180)/2),2)))*1000) AS distance FROM shop having distance <= 5000 order by distance asc
6378.138这个是地球的直径,单位千米.latitude是用户位置的纬度,longitude是用户位置经度.latitude 为商户的纬度字段,longitude为商户的经度字段。上面一段sql计算得出根据用户经纬度计算与商户的距离。