[C#]地球坐标系 (WGS-84) 到火星坐标系 (GCJ-02) 的转换 帮助类

关键代码：

using System;

using YanZhiwei.DotNet2.Utilities.Models;



namespace YanZhiwei.DotNet2.Utilities.Common

{

    /// <summary>

    /// 地球坐标系 (WGS-84) 到火星坐标系 (GCJ-02) 的转换 帮助类 

    /// </summary>

    public class WGSGCJLatLonHelper

    {

        /*

         *参考：

         *1.http://blog.csdn.net/yorling/article/details/9175913

         *2.http://kongxz.com/2013/10/wgs-cgj/

         *3.https://on4wp7.codeplex.com/SourceControl/changeset/view/21483#EvilTransform.cs

         *4.https://github.com/googollee/eviltransform

         *5.https://github.com/shenqiliang/WGS2Mars

         *

         *WGS84:World Geodetic System 1984，是为GPS全球定位系统使用而建立的坐标系统。通过遍布世界的卫星观测站观测到

         *的坐标建立，其初次WGS84的精度为1-2m，在1994年1月2号，通过10个观测站在GPS测量方法上改正，得到了WGS84

         *（G730），G表示由GPS测量得到，730表示为GPS时间第730个周。1996年，National Imagery and Mapping Agency (NIMA) 

         *为美国国防部 (U.S.Departemt of Defense, DoD)做了一个新的坐标系统。这样实现了新的WGS版本：WGS（G873）。其因为

         *加入了USNO站和北京站的改正，其东部方向加入了31-39cm 的改正。所有的其他坐标都有在1分米之内的修正。

         *

         *GCJ-02:国家保密插件，也叫做加密插件或者加偏或者SM模组，其实就是对真实坐标系统进行人为的加偏处理，按照几行代码的算

         *法，将真实的坐标加密成虚假的坐标，而这个加偏并不是线性的加偏，所以各地的偏移情况都会有所不同。而加密后的坐标也常

         *被人称为火星坐标系统。GCJ-02 意味“国测局-2002”，也就是说，这是国家测绘局于2002年弄出的标准。

         */

        #region 常量以及私有方法

        const double pi = 3.14159265358979324;//Math.PI;//;

        const double a = 6378245.0;

        const double ee = 0.00669342162296594323;

        private static double TransformLat(double lonX, double latY)

        {

            double _ret = -100.0 + 2.0 * lonX + 3.0 * latY + 0.2 * latY * latY + 0.1 * lonX * latY + 0.2 * Math.Sqrt(Math.Abs(lonX));

            _ret += (20.0 * Math.Sin(6.0 * lonX * pi) + 20.0 * Math.Sin(2.0 * lonX * pi)) * 2.0 / 3.0;

            _ret += (20.0 * Math.Sin(latY * pi) + 40.0 * Math.Sin(latY / 3.0 * pi)) * 2.0 / 3.0;

            _ret += (160.0 * Math.Sin(latY / 12.0 * pi) + 320 * Math.Sin(latY * pi / 30.0)) * 2.0 / 3.0;

            return _ret;

        }

        private static double TransformLon(double lonX, double latY)

        {

            double _ret = 300.0 + lonX + 2.0 * latY + 0.1 * lonX * lonX + 0.1 * lonX * latY + 0.1 * Math.Sqrt(Math.Abs(lonX));

            _ret += (20.0 * Math.Sin(6.0 * lonX * pi) + 20.0 * Math.Sin(2.0 * lonX * pi)) * 2.0 / 3.0;

            _ret += (20.0 * Math.Sin(lonX * pi) + 40.0 * Math.Sin(lonX / 3.0 * pi)) * 2.0 / 3.0;

            _ret += (150.0 * Math.Sin(lonX / 12.0 * pi) + 300.0 * Math.Sin(lonX / 30.0 * pi)) * 2.0 / 3.0;

            return _ret;

        }

        private static LatLngPoint Transform(LatLngPoint point)

        {

            LatLngPoint _transPoint = new LatLngPoint();

            double _lat = TransformLat(point.LonX - 105.0, point.LatY - 35.0);

            double _lon = TransformLon(point.LonX - 105.0, point.LatY - 35.0);

            double _radLat = point.LatY / 180.0 * pi;

            double _magic = Math.Sin(_radLat);

            _magic = 1 - ee * _magic * _magic;

            double _sqrtMagic = Math.Sqrt(_magic);

            _lat = (_lat * 180.0) / ((a * (1 - ee)) / (_magic * _sqrtMagic) * pi);

            _lon = (_lon * 180.0) / (a / _sqrtMagic * Math.Cos(_radLat) * pi);

            _transPoint.LatY = _lat;

            _transPoint.LonX = _lon;



            return _transPoint;

        }

        #endregion

        #region 火星坐标转 (GCJ-02)地球坐标(WGS-84)

        /// <summary>

        /// 火星坐标转 (GCJ-02)地球坐标(WGS-84)

        /// </summary>

        /// <param name="gcjPoint">火星坐标转 (GCJ-02)</param>

        /// <returns>地球坐标(WGS-84)</returns>

        public static LatLngPoint GCJ02ToWGS84(LatLngPoint gcjPoint)

        {

            if (MapHelper.OutOfChina(gcjPoint))

            {

                return gcjPoint;

            }

            LatLngPoint _transPoint = Transform(gcjPoint);

            return new LatLngPoint(gcjPoint.LatY - _transPoint.LatY, gcjPoint.LonX - _transPoint.LonX);

        }

        #endregion

        #region 地球坐标(WGS-84)转火星坐标 (GCJ-02)

        /// <summary>

        /// 地球坐标(WGS-84)转火星坐标 (GCJ-02)

        /// </summary>

        /// <param name="wgsPoint">地球坐标(WGS-84)</param>

        /// <returns>火星坐标 (GCJ-02)</returns>

        public static LatLngPoint WGS84ToGCJ02(LatLngPoint wgsPoint)

        {

            if (MapHelper.OutOfChina(wgsPoint))

            {

                return wgsPoint;

            }

            LatLngPoint _transPoint = Transform(wgsPoint);

            return new LatLngPoint(wgsPoint.LatY + _transPoint.LatY, wgsPoint.LonX + _transPoint.LonX);

        }

        #endregion



    }

}