GPS坐标互转：WGS-84(GPS)、GCJ-02(Google地图)、BD-09(百度地图) 筋斗云网络

<?php
class GPS {
    private $PI = 3.14159265358979324;
    private $x_pi = 0;
 
    public function __construct()
    {
        $this->x_pi = 3.14159265358979324 * 3000.0 / 180.0;
    }
    //WGS-84 to GCJ-02
    public function gcj_encrypt($wgsLat, $wgsLon) {
        if ($this->outOfChina($wgsLat, $wgsLon))
            return array('lat' => $wgsLat, 'lon' => $wgsLon);
 
        $d = $this->delta($wgsLat, $wgsLon);
        return array('lat' => $wgsLat + $d['lat'],'lon' => $wgsLon + $d['lon']);
    }
    //GCJ-02 to WGS-84
    public function gcj_decrypt($gcjLat, $gcjLon) {
        if ($this->outOfChina($gcjLat, $gcjLon))
            return array('lat' => $gcjLat, 'lon' => $gcjLon);
         
        $d = $this->delta($gcjLat, $gcjLon);
        return array('lat' => $gcjLat - $d['lat'], 'lon' => $gcjLon - $d['lon']);
    }
    //GCJ-02 to WGS-84 exactly
    public function gcj_decrypt_exact($gcjLat, $gcjLon) {
        $initDelta = 0.01;
        $threshold = 0.000000001;
        $dLat = $initDelta; $dLon = $initDelta;
        $mLat = $gcjLat - $dLat; $mLon = $gcjLon - $dLon;
        $pLat = $gcjLat + $dLat; $pLon = $gcjLon + $dLon;
        $wgsLat = 0; $wgsLon = 0; $i = 0;
        while (TRUE) {
            $wgsLat = ($mLat + $pLat) / 2;
            $wgsLon = ($mLon + $pLon) / 2;
            $tmp = $this->gcj_encrypt($wgsLat, $wgsLon);
            $dLat = $tmp['lat'] - $gcjLat;
            $dLon = $tmp['lon'] - $gcjLon;
            if ((abs($dLat) < $threshold) && (abs($dLon) < $threshold))
                break;
 
            if ($dLat > 0) $pLat = $wgsLat; else $mLat = $wgsLat;
            if ($dLon > 0) $pLon = $wgsLon; else $mLon = $wgsLon;
 
            if (++$i > 10000) break;
        }
        //console.log(i);
        return array('lat' => $wgsLat, 'lon'=> $wgsLon);
    }
    //GCJ-02 to BD-09
    public function bd_encrypt($gcjLat, $gcjLon) {
        $x = $gcjLon; $y = $gcjLat;  
        $z = sqrt($x * $x + $y * $y) + 0.00002 * sin($y * $this->x_pi);  
        $theta = atan2($y, $x) + 0.000003 * cos($x * $this->x_pi);  
        $bdLon = $z * cos($theta) + 0.0065;  
        $bdLat = $z * sin($theta) + 0.006; 
        return array('lat' => $bdLat,'lon' => $bdLon);
    }
    //BD-09 to GCJ-02
    public function bd_decrypt($bdLat, $bdLon)
    {
        $x = $bdLon - 0.0065; $y = $bdLat - 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);  
        $gcjLon = $z * cos($theta);  
        $gcjLat = $z * sin($theta);
        return array('lat' => $gcjLat, 'lon' => $gcjLon);
    }
    //WGS-84 to Web mercator
    //$mercatorLat -> y $mercatorLon -> x
    public function mercator_encrypt($wgsLat, $wgsLon)
    {
        $x = $wgsLon * 20037508.34 / 180.;
        $y = log(tan((90. + $wgsLat) * $this->PI / 360.)) / ($this->PI / 180.);
        $y = $y * 20037508.34 / 180.;
        return array('lat' => $y, 'lon' => $x);
        /*
        if ((abs($wgsLon) > 180 || abs($wgsLat) > 90))
            return NULL;
        $x = 6378137.0 * $wgsLon * 0.017453292519943295;
        $a = $wgsLat * 0.017453292519943295;
        $y = 3189068.5 * log((1.0 + sin($a)) / (1.0 - sin($a)));
        return array('lat' => $y, 'lon' => $x);
        //*/
    }
    // Web mercator to WGS-84
    // $mercatorLat -> y $mercatorLon -> x
    public function mercator_decrypt($mercatorLat, $mercatorLon)
    {
        $x = $mercatorLon / 20037508.34 * 180.;
        $y = $mercatorLat / 20037508.34 * 180.;
        $y = 180 / $this->PI * (2 * atan(exp($y * $this->PI / 180.)) - $this->PI / 2);
        return array('lat' => $y, 'lon' => $x);
        /*
        if (abs($mercatorLon) < 180 && abs($mercatorLat) < 90)
            return NULL;
        if ((abs($mercatorLon) > 20037508.3427892) || (abs($mercatorLat) > 20037508.3427892))
            return NULL;
        $a = $mercatorLon / 6378137.0 * 57.295779513082323;
        $x = $a - (floor((($a + 180.0) / 360.0)) * 360.0);
        $y = (1.5707963267948966 - (2.0 * atan(exp((-1.0 * $mercatorLat) / 6378137.0)))) * 57.295779513082323;
        return array('lat' => $y, 'lon' => $x);
        //*/
    }
    // two point's distance
    public function distance($latA, $lonA, $latB, $lonB)
    {
        $earthR = 6371000.;
        $x = cos($latA * $this->PI / 180.) * cos($latB * $this->PI / 180.) * cos(($lonA - $lonB) * $this->PI / 180);
        $y = sin($latA * $this->PI / 180.) * sin($latB * $this->PI / 180.);
        $s = $x + $y;
        if ($s > 1) $s = 1;
        if ($s < -1) $s = -1;
        $alpha = acos($s);
        $distance = $alpha * $earthR;
        return $distance;
    }
 
    private function delta($lat, $lon)
    {
        // Krasovsky 1940
        //
        // a = 6378245.0, 1/f = 298.3
        // b = a * (1 - f)
        // ee = (a^2 - b^2) / a^2;
        $a = 6378245.0;//  a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
        $ee = 0.00669342162296594323;//  ee: 椭球的偏心率。
        $dLat = $this->transformLat($lon - 105.0, $lat - 35.0);
        $dLon = $this->transformLon($lon - 105.0, $lat - 35.0);
        $radLat = $lat / 180.0 * $this->PI;
        $magic = sin($radLat);
        $magic = 1 - $ee * $magic * $magic;
        $sqrtMagic = sqrt($magic);
        $dLat = ($dLat * 180.0) / (($a * (1 - $ee)) / ($magic * $sqrtMagic) * $this->PI);
        $dLon = ($dLon * 180.0) / ($a / $sqrtMagic * cos($radLat) * $this->PI);
        return array('lat' => $dLat, 'lon' => $dLon);
    }
 
    private 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 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 * $this->PI) + 20.0 * sin(2.0 * $x * $this->PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($y * $this->PI) + 40.0 * sin($y / 3.0 * $this->PI)) * 2.0 / 3.0;
        $ret += (160.0 * sin($y / 12.0 * $this->PI) + 320 * sin($y * $this->PI / 30.0)) * 2.0 / 3.0;
        return $ret;
    }
 
    private 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 * $this->PI) + 20.0 * sin(2.0 * $x * $this->PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($x * $this->PI) + 40.0 * sin($x / 3.0 * $this->PI)) * 2.0 / 3.0;
        $ret += (150.0 * sin($x / 12.0 * $this->PI) + 300.0 * sin($x / 30.0 * $this->PI)) * 2.0 / 3.0;
        return $ret;
    }
}

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

147

148

149

150

151

152

153

154

155

156

157

158

159

160

161

162

163

164

<?php

class GPS {

private $PI = 3.14159265358979324;

private $x_pi = 0;

public function __construct()

{

$this->x_pi = 3.14159265358979324 * 3000.0 / 180.0;

}

//WGS-84 to GCJ-02

public function gcj_encrypt($wgsLat, $wgsLon) {

if ($this->outOfChina($wgsLat, $wgsLon))

return array('lat' => $wgsLat, 'lon' => $wgsLon);

$d = $this->delta($wgsLat, $wgsLon);

return array('lat' => $wgsLat + $d['lat'],'lon' => $wgsLon + $d['lon']);

}

//GCJ-02 to WGS-84

public function gcj_decrypt($gcjLat, $gcjLon) {

if ($this->outOfChina($gcjLat, $gcjLon))

return array('lat' => $gcjLat, 'lon' => $gcjLon);

$d = $this->delta($gcjLat, $gcjLon);

return array('lat' => $gcjLat - $d['lat'], 'lon' => $gcjLon - $d['lon']);

}

//GCJ-02 to WGS-84 exactly

public function gcj_decrypt_exact($gcjLat, $gcjLon) {

$initDelta = 0.01;

$threshold = 0.000000001;

$dLat = $initDelta; $dLon = $initDelta;

$mLat = $gcjLat - $dLat; $mLon = $gcjLon - $dLon;

$pLat = $gcjLat + $dLat; $pLon = $gcjLon + $dLon;

$wgsLat = 0; $wgsLon = 0; $i = 0;

while (TRUE) {

$wgsLat = ($mLat + $pLat) / 2;

$wgsLon = ($mLon + $pLon) / 2;

$tmp = $this->gcj_encrypt($wgsLat, $wgsLon);

$dLat = $tmp['lat'] - $gcjLat;

$dLon = $tmp['lon'] - $gcjLon;

if ((abs($dLat) < $threshold) && (abs($dLon) < $threshold))

break;

if ($dLat > 0) $pLat = $wgsLat; else $mLat = $wgsLat;

if ($dLon > 0) $pLon = $wgsLon; else $mLon = $wgsLon;

if (++$i > 10000) break;

}

//console.log(i);

return array('lat' => $wgsLat, 'lon'=> $wgsLon);

}

//GCJ-02 to BD-09

public function bd_encrypt($gcjLat, $gcjLon) {

$x = $gcjLon; $y = $gcjLat;

$z = sqrt($x * $x + $y * $y) + 0.00002 * sin($y * $this->x_pi);

$theta = atan2($y, $x) + 0.000003 * cos($x * $this->x_pi);

$bdLon = $z * cos($theta) + 0.0065;

$bdLat = $z * sin($theta) + 0.006;

return array('lat' => $bdLat,'lon' => $bdLon);

}

//BD-09 to GCJ-02

public function bd_decrypt($bdLat, $bdLon)

{

$x = $bdLon - 0.0065; $y = $bdLat - 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);

$gcjLon = $z * cos($theta);

$gcjLat = $z * sin($theta);

return array('lat' => $gcjLat, 'lon' => $gcjLon);

}

//WGS-84 to Web mercator

//$mercatorLat -> y $mercatorLon -> x

public function mercator_encrypt($wgsLat, $wgsLon)

{

$x = $wgsLon * 20037508.34 / 180.;

$y = log(tan((90. + $wgsLat) * $this->PI / 360.)) / ($this->PI / 180.);

$y = $y * 20037508.34 / 180.;

return array('lat' => $y, 'lon' => $x);

if ((abs($wgsLon) > 180 || abs($wgsLat) > 90))

return NULL;

$x = 6378137.0 * $wgsLon * 0.017453292519943295;

$a = $wgsLat * 0.017453292519943295;

$y = 3189068.5 * log((1.0 + sin($a)) / (1.0 - sin($a)));

return array('lat' => $y, 'lon' => $x);

//*/

}

// Web mercator to WGS-84

// $mercatorLat -> y $mercatorLon -> x

public function mercator_decrypt($mercatorLat, $mercatorLon)

{

$x = $mercatorLon / 20037508.34 * 180.;

$y = $mercatorLat / 20037508.34 * 180.;

$y = 180 / $this->PI * (2 * atan(exp($y * $this->PI / 180.)) - $this->PI / 2);

return array('lat' => $y, 'lon' => $x);

if (abs($mercatorLon) < 180 && abs($mercatorLat) < 90)

return NULL;

if ((abs($mercatorLon) > 20037508.3427892) || (abs($mercatorLat) > 20037508.3427892))

return NULL;

$a = $mercatorLon / 6378137.0 * 57.295779513082323;

$x = $a - (floor((($a + 180.0) / 360.0)) * 360.0);

$y = (1.5707963267948966 - (2.0 * atan(exp((-1.0 * $mercatorLat) / 6378137.0)))) * 57.295779513082323;

return array('lat' => $y, 'lon' => $x);

//*/

}

// two point's distance

public function distance($latA, $lonA, $latB, $lonB)

{

$earthR = 6371000.;

$x = cos($latA * $this->PI / 180.) * cos($latB * $this->PI / 180.) * cos(($lonA - $lonB) * $this->PI / 180);

$y = sin($latA * $this->PI / 180.) * sin($latB * $this->PI / 180.);

$s = $x + $y;

if ($s > 1) $s = 1;

if ($s < -1) $s = -1;

$alpha = acos($s);

$distance = $alpha * $earthR;

return $distance;

}

private function delta($lat, $lon)

{

// Krasovsky 1940

// a = 6378245.0, 1/f = 298.3

// b = a * (1 - f)

// ee = (a^2 - b^2) / a^2;

$a = 6378245.0;// a: 卫星椭球坐标投影到平面地图坐标系的投影因子。

$ee = 0.00669342162296594323;// ee: 椭球的偏心率。

$dLat = $this->transformLat($lon - 105.0, $lat - 35.0);

$dLon = $this->transformLon($lon - 105.0, $lat - 35.0);

$radLat = $lat / 180.0 * $this->PI;

$magic = sin($radLat);

$magic = 1 - $ee * $magic * $magic;

$sqrtMagic = sqrt($magic);

$dLat = ($dLat * 180.0) / (($a * (1 - $ee)) / ($magic * $sqrtMagic) * $this->PI);

$dLon = ($dLon * 180.0) / ($a / $sqrtMagic * cos($radLat) * $this->PI);

return array('lat' => $dLat, 'lon' => $dLon);

}

private 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 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 * $this->PI) + 20.0 * sin(2.0 * $x * $this->PI)) * 2.0 / 3.0;

$ret += (20.0 * sin($y * $this->PI) + 40.0 * sin($y / 3.0 * $this->PI)) * 2.0 / 3.0;

$ret += (160.0 * sin($y / 12.0 * $this->PI) + 320 * sin($y * $this->PI / 30.0)) * 2.0 / 3.0;

return $ret;

}

private 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 * $this->PI) + 20.0 * sin(2.0 * $x * $this->PI)) * 2.0 / 3.0;

$ret += (20.0 * sin($x * $this->PI) + 40.0 * sin($x / 3.0 * $this->PI)) * 2.0 / 3.0;

$ret += (150.0 * sin($x / 12.0 * $this->PI) + 300.0 * sin($x / 30.0 * $this->PI)) * 2.0 / 3.0;

return $ret;

}

var GPS = {
    PI : 3.14159265358979324,
    x_pi : 3.14159265358979324 * 3000.0 / 180.0,
    delta : function (lat, lon) {
        // Krasovsky 1940
        //
        // a = 6378245.0, 1/f = 298.3
        // b = a * (1 - f)
        // ee = (a^2 - b^2) / a^2;
        var a = 6378245.0; //  a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
        var ee = 0.00669342162296594323; //  ee: 椭球的偏心率。
        var dLat = this.transformLat(lon - 105.0, lat - 35.0);
        var dLon = this.transformLon(lon - 105.0, lat - 35.0);
        var radLat = lat / 180.0 * this.PI;
        var magic = Math.sin(radLat);
        magic = 1 - ee * magic * magic;
        var sqrtMagic = Math.sqrt(magic);
        dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * this.PI);
        dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * this.PI);
        return {'lat': dLat, 'lon': dLon};
    },
     
    //WGS-84 to GCJ-02
    gcj_encrypt : function (wgsLat, wgsLon) {
        if (this.outOfChina(wgsLat, wgsLon))
            return {'lat': wgsLat, 'lon': wgsLon};
 
        var d = this.delta(wgsLat, wgsLon);
        return {'lat' : wgsLat + d.lat,'lon' : wgsLon + d.lon};
    },
    //GCJ-02 to WGS-84
    gcj_decrypt : function (gcjLat, gcjLon) {
        if (this.outOfChina(gcjLat, gcjLon))
            return {'lat': gcjLat, 'lon': gcjLon};
         
        var d = this.delta(gcjLat, gcjLon);
        return {'lat': gcjLat - d.lat, 'lon': gcjLon - d.lon};
    },
    //GCJ-02 to WGS-84 exactly
    gcj_decrypt_exact : function (gcjLat, gcjLon) {
        var initDelta = 0.01;
        var threshold = 0.000000001;
        var dLat = initDelta, dLon = initDelta;
        var mLat = gcjLat - dLat, mLon = gcjLon - dLon;
        var pLat = gcjLat + dLat, pLon = gcjLon + dLon;
        var wgsLat, wgsLon, i = 0;
        while (1) {
            wgsLat = (mLat + pLat) / 2;
            wgsLon = (mLon + pLon) / 2;
            var tmp = this.gcj_encrypt(wgsLat, wgsLon)
            dLat = tmp.lat - gcjLat;
            dLon = tmp.lon - gcjLon;
            if ((Math.abs(dLat) < threshold) && (Math.abs(dLon) < threshold))
                break;
 
            if (dLat > 0) pLat = wgsLat; else mLat = wgsLat;
            if (dLon > 0) pLon = wgsLon; else mLon = wgsLon;
 
            if (++i > 10000) break;
        }
        //console.log(i);
        return {'lat': wgsLat, 'lon': wgsLon};
    },
    //GCJ-02 to BD-09
    bd_encrypt : function (gcjLat, gcjLon) {
        var x = gcjLon, y = gcjLat;  
        var z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * this.x_pi);  
        var theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * this.x_pi);  
        bdLon = z * Math.cos(theta) + 0.0065;  
        bdLat = z * Math.sin(theta) + 0.006; 
        return {'lat' : bdLat,'lon' : bdLon};
    },
    //BD-09 to GCJ-02
    bd_decrypt : function (bdLat, bdLon) {
        var x = bdLon - 0.0065, y = bdLat - 0.006;  
        var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * this.x_pi);  
        var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * this.x_pi);  
        var gcjLon = z * Math.cos(theta);  
        var gcjLat = z * Math.sin(theta);
        return {'lat' : gcjLat, 'lon' : gcjLon};
    },
    //WGS-84 to Web mercator
    //mercatorLat -> y mercatorLon -> x
    mercator_encrypt : function(wgsLat, wgsLon) {
        var x = wgsLon * 20037508.34 / 180.;
        var y = Math.log(Math.tan((90. + wgsLat) * this.PI / 360.)) / (this.PI / 180.);
        y = y * 20037508.34 / 180.;
        return {'lat' : y, 'lon' : x};
        /*
        if ((Math.abs(wgsLon) > 180 || Math.abs(wgsLat) > 90))
            return null;
        var x = 6378137.0 * wgsLon * 0.017453292519943295;
        var a = wgsLat * 0.017453292519943295;
        var y = 3189068.5 * Math.log((1.0 + Math.sin(a)) / (1.0 - Math.sin(a)));
        return {'lat' : y, 'lon' : x};
        //*/
    },
    // Web mercator to WGS-84
    // mercatorLat -> y mercatorLon -> x
    mercator_decrypt : function(mercatorLat, mercatorLon) {
        var x = mercatorLon / 20037508.34 * 180.;
        var y = mercatorLat / 20037508.34 * 180.;
        y = 180 / this.PI * (2 * Math.atan(Math.exp(y * this.PI / 180.)) - this.PI / 2);
        return {'lat' : y, 'lon' : x};
        /*
        if (Math.abs(mercatorLon) < 180 && Math.abs(mercatorLat) < 90)
            return null;
        if ((Math.abs(mercatorLon) > 20037508.3427892) || (Math.abs(mercatorLat) > 20037508.3427892))
            return null;
        var a = mercatorLon / 6378137.0 * 57.295779513082323;
        var x = a - (Math.floor(((a + 180.0) / 360.0)) * 360.0);
        var y = (1.5707963267948966 - (2.0 * Math.atan(Math.exp((-1.0 * mercatorLat) / 6378137.0)))) * 57.295779513082323;
        return {'lat' : y, 'lon' : x};
        //*/
    },
    // two point's distance
    distance : function (latA, lonA, latB, lonB) {
        var earthR = 6371000.;
        var x = Math.cos(latA * this.PI / 180.) * Math.cos(latB * this.PI / 180.) * Math.cos((lonA - lonB) * this.PI / 180);
        var y = Math.sin(latA * this.PI / 180.) * Math.sin(latB * this.PI / 180.);
        var s = x + y;
        if (s > 1) s = 1;
        if (s < -1) s = -1;
        var alpha = Math.acos(s);
        var distance = alpha * earthR;
        return distance;
    },
    outOfChina : function (lat, lon) {
        if (lon < 72.004 || lon > 137.8347)
            return true;
        if (lat < 0.8293 || lat > 55.8271)
            return true;
        return false;
    },
    transformLat : function (x, y) {
        var ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(y * this.PI) + 40.0 * Math.sin(y / 3.0 * this.PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(y / 12.0 * this.PI) + 320 * Math.sin(y * this.PI / 30.0)) * 2.0 / 3.0;
        return ret;
    },
    transformLon : function (x, y) {
        var ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(x * this.PI) + 40.0 * Math.sin(x / 3.0 * this.PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(x / 12.0 * this.PI) + 300.0 * Math.sin(x / 30.0 * this.PI)) * 2.0 / 3.0;
        return ret;
    }
};

100

101

102

103

104

105

106

107

108

109

110

111

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

130

131

132

133

134

135

136

137

138

139

140

141

142

143

144

145

146

147

148

149

var GPS = {

PI : 3.14159265358979324,

x_pi : 3.14159265358979324 * 3000.0 / 180.0,

delta : function (lat, lon) {

// Krasovsky 1940

// a = 6378245.0, 1/f = 298.3

// b = a * (1 - f)

// ee = (a^2 - b^2) / a^2;

var a = 6378245.0; // a: 卫星椭球坐标投影到平面地图坐标系的投影因子。

var ee = 0.00669342162296594323; // ee: 椭球的偏心率。

var dLat = this.transformLat(lon - 105.0, lat - 35.0);

var dLon = this.transformLon(lon - 105.0, lat - 35.0);

var radLat = lat / 180.0 * this.PI;

var magic = Math.sin(radLat);

magic = 1 - ee * magic * magic;

var sqrtMagic = Math.sqrt(magic);

dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * this.PI);

dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * this.PI);

return {'lat': dLat, 'lon': dLon};

//WGS-84 to GCJ-02

gcj_encrypt : function (wgsLat, wgsLon) {

if (this.outOfChina(wgsLat, wgsLon))

return {'lat': wgsLat, 'lon': wgsLon};

var d = this.delta(wgsLat, wgsLon);

return {'lat' : wgsLat + d.lat,'lon' : wgsLon + d.lon};

//GCJ-02 to WGS-84

gcj_decrypt : function (gcjLat, gcjLon) {

if (this.outOfChina(gcjLat, gcjLon))

return {'lat': gcjLat, 'lon': gcjLon};

var d = this.delta(gcjLat, gcjLon);

return {'lat': gcjLat - d.lat, 'lon': gcjLon - d.lon};

//GCJ-02 to WGS-84 exactly

gcj_decrypt_exact : function (gcjLat, gcjLon) {

var initDelta = 0.01;

var threshold = 0.000000001;

var dLat = initDelta, dLon = initDelta;

var mLat = gcjLat - dLat, mLon = gcjLon - dLon;

var pLat = gcjLat + dLat, pLon = gcjLon + dLon;

var wgsLat, wgsLon, i = 0;

while (1) {

wgsLat = (mLat + pLat) / 2;

wgsLon = (mLon + pLon) / 2;

var tmp = this.gcj_encrypt(wgsLat, wgsLon)

dLat = tmp.lat - gcjLat;

dLon = tmp.lon - gcjLon;

if ((Math.abs(dLat) < threshold) && (Math.abs(dLon) < threshold))

break;

if (dLat > 0) pLat = wgsLat; else mLat = wgsLat;

if (dLon > 0) pLon = wgsLon; else mLon = wgsLon;

if (++i > 10000) break;

}

//console.log(i);

return {'lat': wgsLat, 'lon': wgsLon};

//GCJ-02 to BD-09

bd_encrypt : function (gcjLat, gcjLon) {

var x = gcjLon, y = gcjLat;

var z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * this.x_pi);

var theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * this.x_pi);

bdLon = z * Math.cos(theta) + 0.0065;

bdLat = z * Math.sin(theta) + 0.006;

return {'lat' : bdLat,'lon' : bdLon};

//BD-09 to GCJ-02

bd_decrypt : function (bdLat, bdLon) {

var x = bdLon - 0.0065, y = bdLat - 0.006;

var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * this.x_pi);

var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * this.x_pi);

var gcjLon = z * Math.cos(theta);

var gcjLat = z * Math.sin(theta);

return {'lat' : gcjLat, 'lon' : gcjLon};

//WGS-84 to Web mercator

//mercatorLat -> y mercatorLon -> x

mercator_encrypt : function(wgsLat, wgsLon) {

var x = wgsLon * 20037508.34 / 180.;

var y = Math.log(Math.tan((90. + wgsLat) * this.PI / 360.)) / (this.PI / 180.);

y = y * 20037508.34 / 180.;

return {'lat' : y, 'lon' : x};

if ((Math.abs(wgsLon) > 180 || Math.abs(wgsLat) > 90))

return null;

var x = 6378137.0 * wgsLon * 0.017453292519943295;

var a = wgsLat * 0.017453292519943295;

var y = 3189068.5 * Math.log((1.0 + Math.sin(a)) / (1.0 - Math.sin(a)));

return {'lat' : y, 'lon' : x};

//*/

// Web mercator to WGS-84

// mercatorLat -> y mercatorLon -> x

mercator_decrypt : function(mercatorLat, mercatorLon) {

var x = mercatorLon / 20037508.34 * 180.;

var y = mercatorLat / 20037508.34 * 180.;

y = 180 / this.PI * (2 * Math.atan(Math.exp(y * this.PI / 180.)) - this.PI / 2);

return {'lat' : y, 'lon' : x};

if (Math.abs(mercatorLon) < 180 && Math.abs(mercatorLat) < 90)

return null;

if ((Math.abs(mercatorLon) > 20037508.3427892) || (Math.abs(mercatorLat) > 20037508.3427892))

return null;

var a = mercatorLon / 6378137.0 * 57.295779513082323;

var x = a - (Math.floor(((a + 180.0) / 360.0)) * 360.0);

var y = (1.5707963267948966 - (2.0 * Math.atan(Math.exp((-1.0 * mercatorLat) / 6378137.0)))) * 57.295779513082323;

return {'lat' : y, 'lon' : x};

//*/

// two point's distance

distance : function (latA, lonA, latB, lonB) {

var earthR = 6371000.;

var x = Math.cos(latA * this.PI / 180.) * Math.cos(latB * this.PI / 180.) * Math.cos((lonA - lonB) * this.PI / 180);

var y = Math.sin(latA * this.PI / 180.) * Math.sin(latB * this.PI / 180.);

var s = x + y;

if (s > 1) s = 1;

if (s < -1) s = -1;

var alpha = Math.acos(s);

var distance = alpha * earthR;

return distance;

outOfChina : function (lat, lon) {

if (lon < 72.004 || lon > 137.8347)

return true;

if (lat < 0.8293 || lat > 55.8271)

return true;

return false;

transformLat : function (x, y) {

var ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));

ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;

ret += (20.0 * Math.sin(y * this.PI) + 40.0 * Math.sin(y / 3.0 * this.PI)) * 2.0 / 3.0;

ret += (160.0 * Math.sin(y / 12.0 * this.PI) + 320 * Math.sin(y * this.PI / 30.0)) * 2.0 / 3.0;

return ret;

transformLon : function (x, y) {

var ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));

ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;

ret += (20.0 * Math.sin(x * this.PI) + 40.0 * Math.sin(x / 3.0 * this.PI)) * 2.0 / 3.0;

ret += (150.0 * Math.sin(x / 12.0 * this.PI) + 300.0 * Math.sin(x / 30.0 * this.PI)) * 2.0 / 3.0;

return ret;

}

};

gcj_encrypt中最后的
return {'lat' : wgsLat + d.lat,'lon' : wgsLon + d.lon};
似乎需要替换为
return {'lat' : parseFloat(wgsLat) + parseFloat(d.lat), 'lon' : parseFloat(wgsLon) + parseFloat(d.lon) };以避免字符串拼接.

gps转换谷歌地图接口：http://ipao.150643.com/aaa.php?a=36.71713&b=119.10078

谷歌定位接口：http://www.google.cn/maps/api/geocode/json?latlng=36.717650500549,119.10667902067&sensor=true

我的微信

把最实用的经验，分享给最需要的读者，希望每一位来访的朋友都能有所收获！

CommentCancel reply

微信