GPS国际坐标(WGS-84)转火星坐标(GCJ-02)----golang版本

package main

 

import (

    "fmt"

    "math"

)

 

// GCJEncrypt ..

func GCJEncrypt(wgsLat, wgsLon float64) (float64, float64) {

    if wgsLon < 72.004 || wgsLon > 137.8347 || wgsLat < 0.8293 || wgsLat > 55.8271 {

        return wgsLat, wgsLon

    }

    PI := 3.14159265358979324

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

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

    var x = wgsLon - 105.0

    var y = wgsLat - 35.0

    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*PI) + 20.0*math.Sin(2.0*x*PI)) * 2.0 / 3.0

    ret += (20.0*math.Sin(y*PI) + 40.0*math.Sin(y/3.0*PI)) * 2.0 / 3.0

    ret += (160.0*math.Sin(y/12.0*PI) + 320*math.Sin(y