根据经纬度求距离(算法)

public function howLang(Request $request){
        $navi_id = $request->get('id');//分类id,或其他查询条件
        $storex = $request->get('longitude');//经度
        $storey = $request->get('latitude');//纬度
        $distance = $request->get('distance');//搜索距离范围
        $result = DB::select("SELECT
                        *,(
                        6371 * acos (
                        cos ( radians({$storey}) )
                        * cos( radians( latitude ) ) 			//表字段中纬度字段
                        * cos( radians( longitude ) - radians({$storex}) )		//表字段中经度字段
                        + sin ( radians({$storey}) )
                        * sin( radians( latitude ) )			//表字段中纬度字段
                        )) AS distance
                        FROM storelist
                        WHERE navi_id = {$navi_id}		//条件，这里是分类id
                        HAVING distance < {$distance}		//一定范围内的地点
                        ORDER BY distance				//通过距离大小排序
                        LIMIT 0 , 20; ");					//分段查询，偏移量设为20
        foreach($result as $k=>$v){         
            $result[$k] = (array)$v;            		      //laravel框架，对象转数组
        }
        $result = $this->distance($result,$storex,$storey);
        return json_encode($result,JSON_UNESCAPED_UNICODE);
    }
    /**
     * @param $result 数组
     * @param $storex 经度
     * @param $storey 纬度
     * @return mixed 计算出距离赋值到数组里
     */
    public function distance($result,$storex,$storey){
        foreach($result as $k=>$v){
            $km = $this->km(round($this->getdistance($storex,$storey,$v['longitude'],$v['latitude'])));
            $result[$k]['km'] =$km ;
        }

        return $result;
    }
    /**
    *显示方式，km 或 m
    */
    public function km($km){
        if($km < 1000){
            return $km.'m';
        }else{
            return ($km/1000).'km';
        }
    }
    /**
     * 求两个已知经纬度之间的距离,单位为米
     * @param lng1 $ ,lng2 经度
     * @param lat1 $ ,lat2 纬度
     * @return float 距离，单位米
     */
    function getdistance($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))) * 6371 * 1000;
        return $s;
    }