C# 高斯投影正算程序：将大地坐标转换为高斯平面直角坐标

下面是一个使用 C# 语言编写的高斯投影正算程序的示例，该程序通过输入点的大地坐标（经度、纬度、大地高），利用高斯投影正算公式计算出对应的高斯平面直角坐标 (X、Y)。其中，使用的参数为国测局 2000 大地坐标系 (CGCS2000) 的参数。

using System;

namespace GaussProjection
{
    class Program
    {
        static void Main(string[] args)
        {
            Console.WriteLine('请输入点的大地坐标（经度，纬度，大地高）：');
            Console.Write('经度：');
            double longitude = double.Parse(Console.ReadLine());
            Console.Write('纬度：');
            double latitude = double.Parse(Console.ReadLine());
            Console.Write('大地高：');
            double height = double.Parse(Console.ReadLine());

            double a = 6378137; // 长半轴
            double f = 1 / 298.257223563; // 扁率
            double b = a * (1 - f); // 短半轴
            double e2 = (a * a - b * b) / (a * a); // 第一偏心率的平方
            double e2_ = (a * a - b * b) / (b * b); // 第二偏心率的平方
            double n = (a - b) / (a + b); // 第三偏心率

            double L = longitude * Math.PI / 180; // 经度转换为弧度
            double B = latitude * Math.PI / 180; // 纬度转换为弧度

            double N = a / Math.Sqrt(1 - e2 * Math.Sin(B) * Math.Sin(B)); // 卯酉圈曲率半径
            double t = Math.Tan(B); // 正切值
            double η2 = e2_ * Math.Cos(B) * Math.Cos(B); // η^2

            double M = a * (1 - e2) / Math.Pow(1 - e2 * Math.Sin(B) * Math.Sin(B), 1.5); // 子午圈曲率半径
            double ρ = a * (1 - e2) / Math.Sqrt(1 - e2 * Math.Sin(B) * Math.Sin(B)); // 曲率半径

            double L0 = Math.Floor(longitude / 3) * 3; // 中央子午线经度
            double l = L - L0; // 经差

            double X = 1000000 + l * N * (1 + l * l * N * (1 - t * t + η2)) / 2; // 高斯平面坐标X
            double Y = l * N * (1 + l * l * N * (5 - t * t + 9 * η2 + 4 * η2 * η2)) / 6; // 高斯平面坐标Y

            Console.WriteLine('高斯平面直角坐标：');
            Console.WriteLine($'X: {X}米');
            Console.WriteLine($'Y: {Y}米');

            Console.ReadLine();
        }
    }
}

注意：该程序没有进行输入验证和异常处理，请确保输入的数据格式正确。