c#编写高斯投影正反算法代码

在 C# 中，编写高斯投影（也称为 UTM 或 Universal Transverse Mercator）的正向和反向转换算法需要利用地理空间库，如.NET Spatial 或是通过.NET Core的 GeoJSON.NET 等。这里是一个简单的示例，展示了如何使用这些库来进行转换：

using System;
using NetTopologySuite.Geometries; // 使用.NET Spatial 库

public static class GaussProjection
{
    public static Point ToProjectedCoordinates(LatLng latLng, int zoneNumber, double falseEasting)
    {
        // 正向转换：从经纬度到UTM坐标
        var point = new GeographyPoint(latLng.Longitude, latLng.Latitude);
        var utmZone = UtmZone.FromZoneNumber(zoneNumber);
        var projectedPoint = point.ToProjected(utmZone.Wgs84ToUtm);

        return new Point(
            projectedPoint.X + falseEasting, // 假定的东部偏移值
            projectedPoint.Y);
    }

    public static LatLng FromProjectedCoordinates(Point projectedPoint, int zoneNumber, double falseNorthing)
    {
        // 反向转换：从UTM坐标到经纬度
        var utmZone = UtmZone.FromZoneNumber(zoneNumber);
        var wgs84Point = projectedPoint.Clone() as Geometry;
        wgs84Point.Transform(utmZone.UtmToWgs84);

        var latLng = new LatLng(wgs84Point.GetY(), wgs84Point.GetX() - falseNorthing); // 假定的北部偏移值

        return latLng;
    }
}

// 示例用法
var latLng = new LatLng(51.5074, -0.1278); // 假设伦敦坐标
int zoneNumber = 30; // 欧洲西部的一个UTM区号
double falseEasting = 300000; // 假定的东部偏移值
Point projectedCoord = GaussProjection.ToProjectedCoordinates(latLng, zoneNumber, falseEasting);

// 反向转换
LatLng unprojectedLatLng = GaussProjection.FromProjectedCoordinates(projectedCoord, zoneNumber, falseEasting);