/**
* 使用方法:
* 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;
? ? }
}
