C# 大地坐标系 转平面坐标系

以下是C#中大地坐标系转平面坐标系的代码示例：

using System;

namespace CoordinateConversion
{
    class Program
    {
        static void Main(string[] args)
        {
            double a = 6378137; // 长半轴
            double b = 6356752.3142; // 短半轴
            double f = (a - b) / a; // 扁率
            double e = Math.Sqrt(2 * f - f * f); // 第一偏心率
            double e2 = Math.Sqrt(f * (2 - f)) / (1 - f); // 第二偏心率
            double L = 116.404; // 经度
            double B = 39.915; // 纬度
            double L0 = 115; // 中央经线
            double k0 = 1; // 比例因子
            double X, Y;
            double N = a / Math.Sqrt(1 - e * e * Math.Sin(B * Math.PI / 180) * Math.Sin(B * Math.PI / 180)); // 卯酉圈曲率半径
            double T = Math.Tan(B * Math.PI / 180) * Math.Tan(B * Math.PI / 180);
            double C = e2 * Math.Cos(B * Math.PI / 180) * Math.Cos(B * Math.PI / 180);
            double A = (L - L0) * Math.PI / 180 * Math.Cos(B * Math.PI / 180);
            double M = a * ((1 - e * e / 4 - 3 * e * e * e * e / 64 - 5 * e * e * e * e * e * e / 256) * B * Math.PI / 180
                - (3 * e * e / 8 + 3 * e * e * e * e / 32 + 45 * e * e * e * e * e * e / 1024) * Math.Sin(2 * B * Math.PI / 180)
                + (15 * e * e * e / 256 + 45 * e * e * e * e * e / 1024) * Math.Sin(4 * B * Math.PI / 180)
                - (35 * e * e * e * e / 3072) * Math.Sin(6 * B * Math.PI / 180));
            X = k0 * N * (A + (1 - T + C) * A * A * A / 6 + (5 - 18 * T + T * T + 72 * C - 58 * e2) * A * A * A * A * A / 120) + 500000;
            Y = k0 * (M + N * Math.Tan(B * Math.PI / 180) * (A * A / 2 + (5 - T + 9 * C + 4 * C * C) * A * A * A * A / 24 + (61 - 58 * T + T * T + 600 * C - 330 * e2) * A * A * A * A * A * A / 720));
            Console.WriteLine("X坐标为：" + X);
            Console.WriteLine("Y坐标为：" + Y);
        }
    }
}

以上代码中，我们首先定义了长半轴a、短半轴b、扁率f、第一偏心率e、第二偏心率e2等参数，然后输入经度L、纬度B、中央经线L0、比例因子k0等参数，最后根据公式计算出平面坐标系的X、Y坐标。