using System;
using System.Collections.Generic;
using System.Text;

namespace Shared.Phone.UserCenter
{
    /// <summary>
    /// GPS的坐标转换逻辑
    /// </summary>
    public class HdlGpsLocationConvertLogic
    {
        #region ■ 变量声明___________________________

        /// <summary>
        /// GPS的坐标转换逻辑
        /// </summary>
        private static HdlGpsLocationConvertLogic m_Current = null;
        /// <summary>
        /// GPS的坐标转换逻辑
        /// </summary>
        public static HdlGpsLocationConvertLogic Current
        {
            get
            {
                if (m_Current == null)
                {
                    m_Current = new HdlGpsLocationConvertLogic();
                }
                return m_Current;
            }
        }

        // 保留的小数位数
        private const string decimals = "f6";

        // 椭球参数-圆周率
        private const double pi = 3.14159265358979324;

        private const double x_pi = pi * 3000.0 / 180.0;

        // (北京54)椭球长半轴，卫星椭球坐标投影到平面地图坐标系的投影因子
        private const double a = 6378245.0;
        /*
            * Krasovsky 1940 (北京54)椭球长半轴第一偏心率平方
            * 计算方式：
            * 长半轴：
            * a = 6378245.0
            * 扁率：
            * 1/f = 298.3（变量相关计算为：(a-b)/a）
            * 短半轴：
            * b = 6356863.0188 (变量相关计算方法为：b = a * (1 - f))
            * 第一偏心率平方:
            * e2 = (a^2 - b^2) / a^2;
        */
        private const double ee = 0.00669342162296594323;
        // 地球半径
        private const double earthR = 6371000.0;

        #endregion

        /// <summary>
        /// 计算偏差
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="dLng"></param>
        /// <param name="dLat"></param>
        private void CalculateDev(double in_lng, double in_lat, out double dLng, out double dLat)
        {
            dLat = TransformLat(in_lng - 105.0, in_lat - 35.0);
            dLng = TransformLng(in_lng - 105.0, in_lat - 35.0);

            double radLat = in_lat / 180.0 * pi;
            double magic = Math.Sin(radLat);
            magic = 1 - ee * magic * magic;
            double sqrtMagic = Math.Sqrt(magic);
            dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
            dLng = (dLng * 180.0) / (a / sqrtMagic * Math.Cos(radLat) * pi);
        }

        #region 高德坐标转WGS84坐标
        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void GCJ02_to_WGS84(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            if (OutOfChina(in_lat, in_lng))
            {
                out_lng = in_lng;
                out_lat = in_lat;
                return;
            }

            CalculateDev(in_lng, in_lat, out out_lng, out out_lat);
            out_lng = in_lng - out_lng;
            out_lat = in_lat - out_lat;
        }

        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] GCJ02_to_WGS84(double lng, double lat)
        {
            double[] lnglat = new double[2];
            GCJ02_to_WGS84(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] GCJ02_to_WGS84(double[] lnglat)
        {
            return GCJ02_to_WGS84(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] GCJ02_to_WGS84(string lng, string lat)
        {
            double _lng, _lat;
            GCJ02_to_WGS84(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] GCJ02_to_WGS84(string[] _lnglat)
        {
            return GCJ02_to_WGS84(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// 高德坐标转WGS84坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void GCJ02_to_WGS84(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            GCJ02_to_WGS84(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }

        /// <summary>
        /// 高德坐标转WGS84坐标 精确(二分极限法)
        /// 默认设置的是精确到小数点后9位，这个值越小，越精确
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        /// <param name="threshold"></param>
        public void GCJ02_to_WGS84_Exact(double in_lng, double in_lat, out double out_lng, out double out_lat, double threshold = 0.000000001)
        {
            double dLng = 0.01, dLat = 0.01;
            double mLng = in_lng - dLng, mLat = in_lat - dLat;
            double pLng = in_lng + dLng, pLat = in_lat + dLat;
            double wgsLng, wgsLat;
            int i = 0;

            while (true)
            {
                wgsLng = (mLng + pLng) / 2;
                wgsLat = (mLat + pLat) / 2;

                GCJ02_to_WGS84(wgsLat, wgsLng, out out_lng, out out_lat);
                dLng = out_lng - in_lng;
                dLat = out_lat - in_lat;

                if ((Math.Abs(dLat) < threshold) && (Math.Abs(dLng) < threshold))
                    break;

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

                if (dLng > 0)
                    pLng = wgsLng;
                else
                    mLng = wgsLng;

                if (++i > 10000)
                    break;
            }
        }
        #endregion

        #region 高德坐标转百度坐标
        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void GCJ02_to_BD09(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            double x = in_lng, y = in_lat;
            double z = Math.Sqrt(x * x + y * y) + 0.00002 * Math.Sin(y * x_pi);
            double theta = Math.Atan2(y, x) + 0.000003 * Math.Cos(x * x_pi);
            out_lng = z * Math.Cos(theta) + 0.0065;
            out_lat = z * Math.Sin(theta) + 0.006;
        }

        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] GCJ02_to_BD09(double lng, double lat)
        {
            double[] lnglat = new double[2];
            GCJ02_to_BD09(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] GCJ02_to_BD09(double[] lnglat)
        {
            return GCJ02_to_BD09(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] GCJ02_to_BD09(string lng, string lat)
        {
            double _lng, _lat;
            GCJ02_to_BD09(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] GCJ02_to_BD09(string[] _lnglat)
        {
            return GCJ02_to_BD09(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// 高德坐标转百度坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void GCJ02_to_BD09(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            GCJ02_to_BD09(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }

        #endregion

        #region 百度坐标转高德坐标
        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void BD09_to_GCJ02(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            double x = in_lng - 0.0065, y = in_lat - 0.006;
            double z = Math.Sqrt(x * x + y * y) - 0.00002 * Math.Sin(y * x_pi);
            double theta = Math.Atan2(y, x) - 0.000003 * Math.Cos(x * x_pi);
            out_lng = z * Math.Cos(theta);
            out_lat = z * Math.Sin(theta);
        }

        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] BD09_to_GCJ02(double lng, double lat)
        {
            double[] lnglat = new double[2];
            BD09_to_GCJ02(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] BD09_to_GCJ02(double[] lnglat)
        {
            return BD09_to_GCJ02(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] BD09_to_GCJ02(string lng, string lat)
        {
            double _lng, _lat;
            BD09_to_GCJ02(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] BD09_to_GCJ02(string[] _lnglat)
        {
            return BD09_to_GCJ02(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// 百度坐标转高德坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void BD09_to_GCJ02(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            BD09_to_GCJ02(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }
        #endregion

        #region 百度坐标转WGS84坐标
        /// <summary>
        /// 百度坐标转WGS84坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void BD09_to_WGS84(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            BD09_to_GCJ02(in_lng, in_lat, out out_lng, out out_lat);
            GCJ02_to_WGS84(out_lng, out_lat, out out_lng, out out_lat);
        }
        #endregion

        #region WGS-84 to Web mercator
        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void WGS84_to_Mercator(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            out_lng = in_lng * 20037508.34 / 180.0;
            double y = Math.Log(Math.Tan((90.0 + in_lat) * pi / 360.0)) / (pi / 180.0);
            out_lat = y * 20037508.34 / 180.0;
        }

        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] WGS84_to_Mercator(double lng, double lat)
        {
            double[] lnglat = new double[2];
            WGS84_to_Mercator(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] WGS84_to_Mercator(double[] lnglat)
        {
            return WGS84_to_Mercator(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] WGS84_to_Mercator(string lng, string lat)
        {
            double _lng, _lat;
            WGS84_to_Mercator(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] WGS84_to_Mercator(string[] _lnglat)
        {
            return WGS84_to_Mercator(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// WGS-84 to Web mercator
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void WGS84_to_Mercator(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            WGS84_to_Mercator(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }
        #endregion

        #region Web mercator to WGS-84
        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void Mercator_to_WGS84(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            out_lng = in_lng / 20037508.34 * 180.0;
            double y = in_lat / 20037508.34 * 180.0;
            out_lat = 180.0 / pi * (2 * Math.Atan(Math.Exp(y * pi / 180.0)) - pi / 2);
        }

        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] Mercator_to_WGS84(double lng, double lat)
        {
            double[] lnglat = new double[2];
            Mercator_to_WGS84(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] Mercator_to_WGS84(double[] lnglat)
        {
            return Mercator_to_WGS84(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] Mercator_to_WGS84(string lng, string lat)
        {
            double _lng, _lat;
            Mercator_to_WGS84(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] Mercator_to_WGS84(string[] _lnglat)
        {
            return Mercator_to_WGS84(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// Web mercator to WGS-84
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void Mercator_to_WGS84(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            Mercator_to_WGS84(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }
        #endregion

        #region WGS84坐标转高德坐标
        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void WGS84_to_GCJ02(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            if (OutOfChina(in_lat, in_lng))
            {
                out_lng = in_lng;
                out_lat = in_lat;
                return;
            }

            CalculateDev(in_lng, in_lat, out out_lng, out out_lat);
            out_lng = in_lng + out_lng;
            out_lat = in_lat + out_lat;
        }

        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        public double[] WGS84_to_GCJ02(double lng, double lat)
        {
            double[] lnglat = new double[2];
            WGS84_to_GCJ02(lng, lat, out lnglat[0], out lnglat[1]);
            return lnglat;
        }

        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="lnglat"></param>
        /// <returns></returns>
        public double[] WGS84_to_GCJ02(double[] lnglat)
        {
            return WGS84_to_GCJ02(lnglat[0], lnglat[1]);
        }

        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="lng"></param>
        /// <param name="lat"></param>
        /// <returns></returns>
        public string[] WGS84_to_GCJ02(string lng, string lat)
        {
            double _lng, _lat;
            WGS84_to_GCJ02(ToDouble(lng), ToDouble(lat), out _lng, out _lat);
            return new string[] { _lng.ToString(decimals), _lat.ToString(decimals) };
        }

        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="_lnglat"></param>
        /// <returns></returns>
        public string[] WGS84_to_GCJ02(string[] _lnglat)
        {
            return WGS84_to_GCJ02(_lnglat[0], _lnglat[1]);
        }

        /// <summary>
        /// WGS84坐标转高德坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void WGS84_to_GCJ02(string in_lng, string in_lat, out string out_lng, out string out_lat)
        {
            double _lng, _lat;
            WGS84_to_GCJ02(ToDouble(in_lng), ToDouble(in_lat), out _lng, out _lat);
            out_lng = _lng.ToString(decimals);
            out_lat = _lat.ToString(decimals);
        }
        #endregion

        #region WGS84坐标转百度坐标
        /// <summary>
        /// WGS84坐标转百度坐标
        /// </summary>
        /// <param name="in_lng"></param>
        /// <param name="in_lat"></param>
        /// <param name="out_lng"></param>
        /// <param name="out_lat"></param>
        public void WGS84_to_BD09(double in_lng, double in_lat, out double out_lng, out double out_lat)
        {
            WGS84_to_GCJ02(in_lng, in_lat, out out_lng, out out_lat);
            GCJ02_to_BD09(out_lng, out_lat, out out_lng, out out_lat);
        }
        #endregion

        /// <summary>
        /// 两点直接的距离
        /// </summary>
        public double GetDistance(double latA, double lonA, double latB, double lonB)
        {
            double x = Math.Cos(latA * pi / 180.0) * Math.Cos(latB * pi / 180.0) * Math.Cos((lonA - lonB) * pi / 180.0);
            double y = Math.Sin(latA * pi / 180.0) * Math.Sin(latB * pi / 180.0);
            double s = x + y;
            if (s > 1) s = 1;
            if (s < -1) s = -1;
            double alpha = Math.Acos(s);
            double distance = alpha * earthR;
            return distance;
        }

        /// <summary>
        /// 坐标是否在中国境内
        /// </summary>
        /// <param name="lat"></param>
        /// <param name="lng"></param>
        /// <returns></returns>
        public bool OutOfChina(double lat, double lng)
        {
            if (lng < 72.004 || lng > 137.8347)
                return true;
            if (lat < 0.8293 || lat > 55.8271)
                return true;
            return false;
        }

        private double TransformLat(double x, double y)
        {
            double 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 * pi / 30.0)) * 2.0 / 3.0;
            return ret;
        }

        private double TransformLng(double x, double y)
        {
            double 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 * pi) + 20.0 * Math.Sin(2.0 * x * pi)) * 2.0 / 3.0;
            ret += (20.0 * Math.Sin(x * pi) + 40.0 * Math.Sin(x / 3.0 * pi)) * 2.0 / 3.0;
            ret += (150.0 * Math.Sin(x / 12.0 * pi) + 300.0 * Math.Sin(x / 30.0 * pi)) * 2.0 / 3.0;
            return ret;
        }

        /// <summary>
        /// 将string转换为十进制双精度数，如果格式错误将转换失败并返回默认值
        /// </summary>
        /// <param name="str"></param>
        /// <returns></returns>
        private double ToDouble(string str, double defaultValue = default)
        {
            if (double.TryParse(str, out double result))
            {
                return result;
            }

            return defaultValue;
        }
    }
}