關於兩個經緯度之間的距離計算【更新至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)
    );
}