关于两个经纬度之间的距离计算【更新至2021.01.04】

一般来说 我们通常见到的经纬度计算时两点之间的直线距离 随着app对接地图，公众号小程序等对接地图插件，直接显示的是驾车或者骑车的曲线距离，所以两者无法等同需要注意

一下提供两种直线距离的写法

调取前请先注意一下问题

经度 longitude
纬度 latitude

 数据调用不出来或者距离不对 先考虑是不是经纬度写反了喔

1 公共方法调用

/**
 * 计算两点地理坐标之间的距离
 * @param Decimal $longitude1 起点经度
 * @param Decimal $latitude1 起点纬度
 * @param Decimal $longitude2 终点经度
 * @param Decimal $latitude2 终点纬度
 * @param Int   $unit    单位 1:米 2:公里
 * @param Int   $decimal  精度 保留小数位数
 * @return Decimal
 */

function getDistance($longitude1, $latitude1, $longitude2, $latitude2, $unit=1, $decimal=2){
    $EARTH_RADIUS = 6370.996; // 地球半径系数
    $PI = 3.1415926;
    $radLat1 = $latitude1 * $PI / 180.0;
    $radLat2 = $latitude2 * $PI / 180.0;
    $radLng1 = $longitude1 * $PI / 180.0;
    $radLng2 = $longitude2 * $PI /180.0;
    $a = $radLat1 - $radLat2;
    $b = $radLng1 - $radLng2;
    $distance = 2 * asin(sqrt(pow(sin($a/2),2) + cos($radLat1) * cos($radLat2) * pow(sin($b/2),2)));
    $distance = $distance * $EARTH_RADIUS * 1000;
    if($unit==2){
        $distance = $distance / 1000;
    }
    return round($distance, $decimal);
}

　　公共方法调用2  

function getdistance2($lng1,$lat1,$lng2,$lat2){
    //将角度转为狐度
    $radLat1=deg2rad($lat1);//deg2rad()函数将角度转换为弧度
    $radLat2=deg2rad($lat2);
    $radLng1=deg2rad($lng1);
    $radLng2=deg2rad($lng2);
    $a=$radLat1-$radLat2;
    $b=$radLng1-$radLng2;
    $s=2*asin(sqrt(pow(sin($a/2),2)+cos($radLat1)*cos($radLat2)*pow(sin($b/2),2)))*6378.137*1000;
    return $s;
}

2 分页显示距离远近数据

public function list(){
        $data = $this->request->param();
        $page=$data['page'];
        $where['status']= 1;

        $lat=$data['lat'];
        $lng=$data['lng'];

        $order='mi asc';
        $list= StoreApply::field("*,ROUND(6378.138*2*ASIN(SQRT(POW(SIN(({$lat}*PI()/180-lng*PI()/180)/2),2)+COS({$lat}*PI()/180)*COS(lng*PI()/180)*POW(SIN(({$lng}*PI()/180-lat*PI()/180)/2),2)))*1000) as mi")
            ->where($where)
            ->order($order)
            ->page($page,16)
            ->select();

        $this->success('请求成功',compact('list'));
    }     

关于sql语句，如果看不懂的 也可以直接以下方sql带入

ROUND(6370.996 * 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 mi

有兴趣的可以直接研究下自定义的函数与上面这个sql语句的带入关系

以上两种形式的算法有轻微差别 大差不差

另外涉及方法

/*
 * 获取当前位置附近N范围
 * */

function SquarePoint($lng, $lat,$distance = 5){
    $dlng =  2 * asin(sin($distance / (2 * 6371)) / cos(deg2rad($lat)));
    $dlng = rad2deg($dlng);
    $dlat = $distance/6371;
    $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)
    );
}