根據經緯度,獲取'MBR'數據

/**

* 使用方法:

* 1.根據經緯度,獲取'MBR'數據

* @param float $lat 緯度

* @param float $lon 經度

* @param float $distance 距離(單位公里 如5就是5公里)

* @return array 'MBR'數據

* 使用實例:

* $gpsdis = (new Getmbr)->Main($lat,$lon,$distance)

*

* 2.根據經緯度,計算2個點之間的距離。

* @param float $lat 緯度1

* @param float $lon 經度1

* @param float $lat 緯度2

* @param float $lon 經度2

* @return float 距離(公里、千米)

* 使用實例:

* $distance = (new Getmbr)->Distance($lat1, $lon1, $lat2, $lon2)

*

* 說明:

* 這個根據一個經緯度坐標、距離然后求另外一個經緯度坐標的作用,

* 主要就是確定一個最小外包矩形(Minimum bounding rectangle,簡稱MBR)。

* 例如,我要找一個坐標點(lat,lon)的5公里范圍內的所有商戶信息、景點信息等。

* 這個MBR就是一個最大的范圍,

* 這個矩形是包含5公里范圍內所有這些有效信息的一個最小矩形。

* 利用公式,求出四個方向0度、90度、180度、270度方向上的四個坐標點就可以得到這個MBR。

*

* 如果有一個應用,表里存有100萬的數據,數據包含一個lat、lon的經緯度信息。

* 就可以先根據輸入的經緯度和距離得到一個MBR,然后通過類似

SELECT Id

FROM IdInfoTable

WHERE latitude >= minLat AND latitude < maxLat

AND longitude >= minLon AND longitude < maxLon

*

*/

namespace app\index\controller;

class Getmbr {

//最大緯度值

? ? public $MaxLatitude;

? ? //最小緯度值

? ? public $MinLatitude;

? ? public $MaxLongitude;

? ? public $MinLongitude;

? ? public $Ea = 6378137;? ? //? 赤道半徑

? ? public $Eb = 6356725;? ? //? 極半徑

? ? const EARTH_RADIUS = 6371.0;//km 地球半徑 平均值,千米

/**

* 獲取'MBR'數據

* @param float $lat 緯度

* @param float $lon 經度

* @param float $distance 距離(單位公里 如5就是5公里)

* @return array 'MBR'數據

* */

? ? public function Main($lat,$lon,$distance)

{

$this->GetMBR($lat, $lon, $distance);

? ? ? ? //返回'MBR'數據

? ? ? ? return [

'MaxLatitude'=>$this->MaxLatitude,

? ? ? ? ? ? 'MinLatitude'=>$this->MinLatitude,

? ? ? ? ? ? 'MaxLongitude'=>$this->MaxLongitude,

? ? ? ? ? ? 'MinLongitude'=>$this->MinLongitude

];

? ? }

private function GetMBR($centorlatitude, $centorLogitude, $distance)

{

//以下為核心代碼

? ? ? ? $range? = 180 / pi()* $distance / 6372.797; //里面的 $distance 就代表搜索 $distance 之內,單位km

? ? ? ? $lngR? = $range / cos($centorlatitude * pi()/ 180);

? ? ? ? $this->MaxLatitude= $centorlatitude + $range; //最大緯度

? ? ? ? $this->MinLatitude= $centorlatitude - $range; //最小緯度

? ? ? ? $this->MaxLongitude= $centorLogitude + $lngR; //最大經度

? ? ? ? $this->MinLongitude= $centorLogitude - $lngR; //最小經度

//得出這四個值以后,就可以根據你數據庫里存的經緯度信息查找記錄了~

? ? }

/**

* 將角度換算為弧度。

* @param float $degrees 角度

* @return float 弧度

* */

? ? private function ConvertDegreesToRadians($degrees)

{

return $degrees * M_PI / 180;

? ? }

/**

* 根據經緯度,計算2個點之間的距離。

* @param float $lat 緯度1

* @param float $lon 經度1

* @param float $lat 緯度2

* @param float $lon 經度2

* @return float 距離(公里、千米)

* */

? ? public function Distance($lat1, $lon1, $lat2, $lon2){

//用haversine公式計算球面兩點間的距離。

//經緯度轉換成弧度

? ? ? ? $lat1 = $this->ConvertDegreesToRadians($lat1);

? ? ? ? $lon1 = $this->ConvertDegreesToRadians($lon1);

? ? ? ? $lat2 = $this->ConvertDegreesToRadians($lat2);

? ? ? ? $lon2 = $this->ConvertDegreesToRadians($lon2);

? ? ? ? //差值

? ? ? ? $vLon = abs($lon1 - $lon2);

? ? ? ? $vLat = abs($lat1 - $lat2);

? ? ? ? //h is the great circle distance in radians, great circle就是一個球體上的切面,它的圓心即是球心的一個周長最大的圓。

? ? ? ? $h = $this->HaverSin($vLat)+ cos($lat1)* cos($lat2)* $this->HaverSin($vLon);

? ? ? ? $distance = 2 * self::EARTH_RADIUS * asin(sqrt($h));

? ? ? ? return $distance;

? ? }

private function HaverSin($theta)

{

$v = sin($theta / 2);

? ? ? ? return $v * $v;

? ? }

}

?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。

推薦閱讀更多精彩內容