C#转84坐标

public static class CoordinateTransformation
{
    private static double pi = 3.1415926535897932384626; //π
    private static double a = 6378245.0; // semi-major axis
    private static double ee = 0.00669342162296594323; // oblateness

 
    public static (double lng, double lat) Gcj02toWgs84(double lng, double lat)
    {
        // GCJ02->GPS84
	// Param lng: longitude in GCJ02
        // Param lat: latitude in GCJ02 : 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;
    }
}

 

posted on 2025-05-16 17:32  大大的橙子  阅读(19)  评论(0)    收藏  举报

导航