地球坐标,火星坐标,百度坐标转换的php实现

coordinate.class.php
<?php

class Coordinate
{
	public  $x = 0;//lat
	public  $y = 0; // lng
	public  $z = 0; // other
	
	function __construct($lat,$lng)
	{
		$this->x = $lat;
		$this->y = $lng;
	}
}
?>

GeoTransUtil.class.php

<?php

require 'Coordinate.class.php';
define('X_PI',3.14159265358979324*3000.0/180.0);

class GeoTransUtil
{
	private static $x_pi = X_PI;
	private static $pi = 3.14159265358979324;
	private static $a = 6378245.0;
	private static $ee = 0.00669342162296594323;
	
//火星坐标(GCJ02坐标,高德,谷歌,腾讯坐标)到百度坐标BD-09
	public static function gcjTObd($from){
		$to = new Coordinate();
		$x = $from->y;
		$y = $from->x;
		$z = sqrt($x * $x + $y * $y) + 0.00002 * sin($y * self::$x_pi);  
		$theta = atan2($y, $x) + 0.000003 * cos($x * self::$x_pi);
		$to->y = $z * cos($theta) + 0.0065; 
		$to->x = $z * sin($theta) + 0.006;
		
		return ($to);
	}
	
	//百度坐标BD-09到火星坐标GCJ02(高德,谷歌,腾讯坐标)
	public static function bdTOgcj($from){
		$to = new Coordinate();
		$x = $from->y - 0.0065;
		$y = $from->x - 0.006;
		$z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * self::$x_pi); 
		$theta = atan2($y, $x) - 0.000003 * cos($x * self::$x_pi);
		$to->y = $z * cos($theta);
		$to->x = $z * sin($theta);
		return ($to);
	}
	
	// WGS-84(GPS坐标,谷歌地球坐标) 到 GCJ-02(火星坐标) 的转换
    public static function wgsTOgcj($from){
    	//double wgLat, double wgLon, out double mgLat, out double mgLon
    		$wgLat = $from->x;
    		$wgLon = $from->y;
        if (self::outOfChina($wgLat, $wgLon)){
            return (new Coordinate($wgLat, $wgLon));
        }
        
        $dLat = GeoTransUtil::transformLat($wgLon - 105.0, $wgLat - 35.0);
        $dLon = GeoTransUtil::transformLon($wgLon - 105.0, $wgLat - 35.0);
        $radLat = $wgLat / 180.0 * self::$pi;
        $magic = sin($radLat);
        $magic = 1 - self::$ee * $magic * $magic;
        $sqrtMagic = sqrt($magic);
        $dLat = ($dLat * 180.0) / ((self::$a * (1 - self::$ee)) / ($magic * $sqrtMagic) * self::$pi);
        $dLon = ($dLon * 180.0) / (self::$a / $sqrtMagic * cos($radLat) * self::$pi);
        
        $mgLat = $wgLat + $dLat;
        $mgLon = $wgLon + $dLon;
        
        return (new Coordinate($mgLat, $mgLon));
    }

	// GCJ-02 到 WGS-84 的转换
    public static function gcjTOwgs($from){
    	$to = self::wgsTOgcj($from);
    	$lat = $from->x;
    	$lon = $from->y;
    	$g_lat = $to->x;
    	$g_lon = $to->y;
    	$d_lat = $g_lat - $lat;
    	$d_lon = $g_lon - $lon;
    	
    	return new Coordinate($lat - $d_lat, $lon - $d_lon);
    }
    
    / //把地球经纬度转换为lat和lng
    public static function earthCoordinate($t1, $t2) {
    	$tt1 = $t1->x + $t1->y/60.0 + $t1->z/3600.0;
    	$tt2 = $t2->x + $t2->y/60.0 + $t2->z/3600.0;
    	
    	return (new Coordinate($tt1, $tt2));
    }
    
    private static function outOfChina($lat,$lon)
    {
        if ($lon < 72.004 || $lon > 137.8347)
            return true;
        if ($lat < 0.8293 || $lat > 55.8271)
            return true;
        
        return false;
    }

    private static function transformLat($x,$y)
    {
        $ret = -100.0 + 2.0 * $x + 3.0 * $y + 0.2 * $y * $y + 0.1 * $x * $y + 0.2 * sqrt(abs($x));
        $ret += (20.0 * sin(6.0 * $x * self::$pi) + 20.0 * sin(2.0 * $x * self::$pi)) * 2.0 / 3.0;
        $ret += (20.0 * sin($y * self::$pi) + 40.0 * sin($y / 3.0 * self::$pi)) * 2.0 / 3.0;
        $ret += (160.0 * sin($y / 12.0 * self::$pi) + 320 * sin($y * self::$pi / 30.0)) * 2.0 / 3.0;
        return $ret;
    }

    private static function transformLon($x, $y)
    {
        $ret = 300.0 + $x + 2.0 * $y + 0.1 * $x * $x + 0.1 * $x * $y + 0.1 * sqrt(abs($x));
        $ret += (20.0 * sin(6.0 * $x * self::$pi) + 20.0 * sin(2.0 * $x * self::$pi)) * 2.0 / 3.0;
        $ret += (20.0 * sin($x * self::$pi) + 40.0 * sin($x / 3.0 * self::$pi)) * 2.0 / 3.0;
        $ret += (150.0 * sin($x / 12.0 * self::$pi) + 300.0 * sin($x / 30.0 * self::$pi)) * 2.0 / 3.0;
        return $ret;
    }
    
}




?>


cor_test.php

<?php

require 'GeoTransUtil.class.php';

$from = new Coordinate(55.8271, 137.8345);

$ret = GeoTransUtil::wgsTOgcj($from);
echo '</br>';
echo $ret->x.','.$ret->y;
?>



  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值