PHP 坐标转换

增加了GCJ02转换为WGS84(高德转北斗)

class cord_convert
{
	// π圆周率*3000/180
	private $x_PI = 52.3598775598;
	// π圆周率
	private $PI = 3.1415926535897932384626;
	// 地球长半轴
	private $a = 6378245.0;
	// 地球扁率
	private $ee = 0.00669342162296594323;

	/**
	 * 判断是否在国内，不在国内则不做偏移
	 * @param   float  $lng  经度
	 * @param   float  $lat  纬度
	 * @return  boolean
	 */
	private function out_of_china($lng, $lat)
	{
		return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
	}

	/**
	 * WGS84转GCJ02(北斗转高德)
	 * @param   float  $lng  经度
	 * @param   float  $lat  纬度
	 * @return  array        处理之后的GCJ02坐标
	 */
	public function wgs84togcj02($lng, $lat)
	{
		if ($this->out_of_china($lng, $lat)) {
			return array($lng, $lat);
		} else {
			$dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
			$dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
			$radlat = $lat / 180.0 * $this->PI;
			$magic = sin($radlat);
			$magic = 1 - $this->ee * $magic * $magic;
			$sqrtmagic = sqrt($magic);
			$dlat = ($dlat * 180.0) / (($this->a * (1 - $this->ee)) / ($magic * $sqrtmagic) * $this->PI);
			$dlng = ($dlng * 180.0) / ($this->a / $sqrtmagic * cos($radlat) * $this->PI);
			$mglng = $lng + $dlng;
			$mglat = $lat + $dlat;
			return array($mglng, $mglat);
		}
	}

	/**
	 * GCJ02转换为WGS84(高德转北斗)
	 * @param   float  $lng  经度
	 * @param   float  $lat  纬度
	 * @return  array        处理之后的WGS84坐标
	 */
	public function gcj02towgs84($lng, $lat)
	{
		if ($this->out_of_china($lng, $lat)) {
			return array($lng, $lat);
		} else {
			$dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
			$dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
			$radlat = $lat / 180.0 * $this->PI;
			$magic = sin($radlat);
			$magic = 1 - $this->ee * $magic * $magic;
			$sqrtmagic = sqrt($magic);
			$dlat = ($dlat * 180.0) / (($this->a * (1 - $this->ee)) / ($magic * $sqrtmagic) * $this->PI);
			$dlng = ($dlng * 180.0) / ($this->a / $sqrtmagic * cos($radlat) * $this->PI);
			$mglng = $lng + $dlng;
			$mglat = $lat + $dlat;
			return array($lng * 2 - $mglng, $lat * 2 - $mglat);
		}
	}

	private function transformlat($lng, $lat)
	{
		$ret = -100.0 + 2.0 * $lng + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lng * $lat + 0.2 * sqrt(abs($lng));
		$ret += (20.0 * sin(6.0 * $lng * $this->PI) + 20.0 * sin(2.0 * $lng * $this->PI)) * 2.0 / 3.0;
		$ret += (20.0 * sin($lat * $this->PI) + 40.0 * sin($lat / 3.0 * $this->PI)) * 2.0 / 3.0;
		$ret += (160.0 * sin($lat / 12.0 * $this->PI) + 320 * sin($lat * $this->PI / 30.0)) * 2.0 / 3.0;
		return $ret;
	}

	private function transformlng($lng, $lat)
	{
		$ret = 300.0 + $lng + 2.0 * $lat + 0.1 * $lng * $lng + 0.1 * $lng * $lat + 0.1 * sqrt(abs($lng));
		$ret += (20.0 * sin(6.0 * $lng * $this->PI) + 20.0 * sin(2.0 * $lng * $this->PI)) * 2.0 / 3.0;
		$ret += (20.0 * sin($lng * $this->PI) + 40.0 * sin($lng / 3.0 * $this->PI)) * 2.0 / 3.0;
		$ret += (150.0 * sin($lng / 12.0 * $this->PI) + 300.0 * sin($lng / 30.0 * $this->PI)) * 2.0 / 3.0;
		return $ret;
	}

	/**
	 * 火星坐标系（GCJ-02）与百度坐标系（BD-09）的转换
	 * 即谷歌、高德 转 百度
	 * @param   float  $lng  经度
	 * @param   float  $lat  纬度
	 * @return  array        处理之后的百度坐标
	 */
	public function gcj02tobd09($lng, $lat)
	{
		if ($this->out_of_china($lng, $lat)) {
			return array($lng, $lat);
		} else {
			$z = sqrt($lng * $lng + $lat * $lat) + 0.00002 * sin($lat * $this->x_PI);
			$theta = atan2($lat, $lng) + 0.000003 * cos($lng * $this->x_PI);
			$mglng = $z * cos($theta) + 0.0065;
			$mglat = $z * sin($theta) + 0.006;
			return array($mglng, $mglat);
		}
	}

	/**
	 * 百度坐标系（BD-09）与火星坐标系（GCJ-02）的转换
	 * 即百度 转 谷歌、高德
	 * @param   float  $lng  经度
	 * @param   float  $lat  纬度
	 * @return  array        处理之后的GCJ02坐标
	 */
	public function bd09togcj02($lng, $lat)
	{
		if ($this->out_of_china($lng, $lat)) {
			return array($lng, $lat);
		} else {
			$x = (float)$lng - 0.0065;
			$y = (float)$lat - 0.006;
			$z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $this->x_PI);
			$theta = atan2($y, $x) - 0.000003 * cos($x * $this->x_PI);
			$mglng = $z * cos($theta);
			$mglat = $z * sin($theta);
			return array($mglng, $mglat);
		}
	}
}