[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    }}