C# 火星坐标系与WGS84坐标系转换方法

话不多说，直接上代码：

public static class CoordinateTransformation
{
    private static double pi = 3.1415926535897932384626; //π
    private static double a = 6378245.0; // 长半轴
    private static double ee = 0.00669342162296594323; // 扁率

    public static (double lng, double lat) Wgs84toGcj02(double lng, double lat)
    {
        var dlat = Transformlat(lng - 105.0, lat - 35.0);
        var dlng = Transformlng(lng - 105.0, lat - 35.0);
        var radlat = (lat / 180.0) * 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)) * pi);
        dlng = (dlng * 180.0) / ((a / sqrtmagic) * Math.Cos(radlat) * pi);
        var mglat = lat + dlat;
        var mglng = lng + dlng;
        return (mglng, mglat);
    }

    public static (double lng, double lat) Gcj02toWgs84(double lng, double lat)
    {
        // GCJ02(火星坐标系) 转GPS84: param lng: 火星坐标系的经度: param lat: 火星坐标系纬度: return :
        var dlat = Transformlat(lng - 105.0, lat - 35.0);
        var dlng = Transformlng(lng - 105.0, lat - 35.0);
        var radlat = (lat / 180.0) * 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)) * pi);
        dlng = (dlng * 180.0) / ((a / sqrtmagic) * Math.Cos(radlat) * pi);
        var mglat = lat + dlat;
        var mglng = lng + dlng;
        return (lng * 2 - mglng, lat * 2 - mglat);
    }

    private static double Transformlat(double lng, double lat)
    {
        var ret =
          -100.0 +
          2.0 * lng +
          3.0 * lat +
          0.2 * lat * lat +
          0.1 * lng * lat +
          0.2 * Math.Sqrt(Math.Abs(lng));
        ret +=
          ((20.0 * Math.Sin(6.0 * lng * pi) +
            20.0 * Math.Sin(2.0 * lng * pi)) *
            2.0) /
          3.0;
        ret +=
          ((20.0 * Math.Sin(lat * pi) +
            40.0 * Math.Sin((lat / 3.0) * pi)) *
            2.0) /
          3.0;
        ret +=
          ((160.0 * Math.Sin((lat / 12.0) * pi) +
            320 * Math.Sin((lat * pi) / 30.0)) *
            2.0) /
          3.0;
        return ret;
    }

    private static double Transformlng(double lng, double lat)
    {
        var ret =
          300.0 +
          lng +
          2.0 * lat +
          0.1 * lng * lng +
          0.1 * lng * lat +
          0.1 * Math.Sqrt(Math.Abs(lng));
        ret +=
          ((20.0 * Math.Sin(6.0 * lng * pi) +
            20.0 * Math.Sin(2.0 * lng * pi)) *
            2.0) /
          3.0;
        ret +=
          ((20.0 * Math.Sin(lng * pi) +
            40.0 * Math.Sin((lng / 3.0) * pi)) *
            2.0) /
          3.0;
        ret +=
          ((150.0 * Math.Sin((lng / 12.0) * pi) +
            300.0 * Math.Sin((lng / 30.0) * pi)) *
            2.0) /
          3.0;
        return ret;
    }
}