c# - Lambert 72 转换为 Lat/Lon 对

Question

我正在尝试合并两个数据库以合并两个客户的网站。但是，客户端 A 一直使用常规的纬度/经度对进行地理定位，而客户端 B 使用兰伯特 72 (X/Y) 坐标。

我已经构建了一个应该转换这些坐标的脚本(因为我不确定最终合并的数据库中将使用哪个坐标，我正在尝试以任何一种方式转换它们)。
我从这里拿了一些片段:http://zoologie.umh.ac.be/tc/algorithms.aspx

请注意，下面提到的所有坐标都指向比利时的位置。

我正在转换一些坐标以查看计算是否正确，但我得到的坐标似乎很遥远。作为引用，比利时的中心大致是(北 50.84323737103243，东 4.355735778808594)，所以我希望所有坐标都接近这些值。

我将 Lambert 72 值 (X: 151488250, Y: 170492909) 转换为 Lat/Lon 对，但结果是: (-87.538.... , -50.724....) 这与预期值相差甚远.
如果我转换整圈(Lambert->LatLon->Lambert，反之亦然)，我会得到与输入相同的结果值，所以我知道我的转换至少是一致的，并且转换是彼此完美的反转。

我也尝试了一些在线转换器工具，它们给了我相同的 (-87.538...., -50.724....) 结果。

由于多个来源产生相同的结果，并且我的转换是彼此的正确反转，我认为计算本身是正确的，但结果值仍然需要进一步转换/偏移？

我认为自己在代数方面已经足够了，但我完全无法使用制图投影。

有人可以对此有所了解吗？

额外信息

我希望我在正确的论坛上发布了这个。我不太确定把它放在哪里，因为这是地理、数学和编码/转换的混合...

提到的兰伯特坐标(X:151488250，Y:170492909)指向布鲁塞尔的一个位置，因此纬度/经度结果应该非常接近(北 50.84323737103243，东 4.355735778808594)。

请在下面找到我的转换函数:

public static Lambert72 LatLon_To_Lambert72(LatLon latlon)
{
    var lat = latlon.Lat;
    var lng = latlon.Lon;

    double LongRef = 0.076042943;
    //=4°21'24"983
    double bLamb = 6378388 * (1 - (1 / 297));
    double aCarre = Math.Pow(6378388, 2);
    double eCarre = (aCarre - Math.Pow(bLamb, 2)) / aCarre;
    double KLamb = 11565915.812935;
    double nLamb = 0.7716421928;

    double eLamb = Math.Sqrt(eCarre);
    double eSur2 = eLamb / 2;

    //conversion to radians
    lat = (Math.PI / 180) * lat;
    lng = (Math.PI / 180) * lng;

    double eSinLatitude = eLamb * Math.Sin(lat);
    double TanZDemi = (Math.Tan((Math.PI / 4) - (lat / 2))) * (Math.Pow(((1 + (eSinLatitude)) / (1 - (eSinLatitude))), (eSur2)));

    double RLamb = KLamb * (Math.Pow((TanZDemi), nLamb));

    double Teta = nLamb * (lng - LongRef);

    double x = 0;
    double y = 0;

    x = 150000 + 0.01256 + RLamb * Math.Sin(Teta - 0.000142043);
    y = 5400000 + 88.4378 - RLamb * Math.Cos(Teta - 0.000142043);

    return new Lambert72(x, y);
}

public static LatLon Lambert72_To_LatLon(Lambert72 lb72)
{
    double X = lb72.X;
    double Y = lb72.Y;

    double LongRef = 0.076042943;
    //=4°21'24"983
    double nLamb = 0.7716421928;
    double aCarre = Math.Pow(6378388, 2);
    double bLamb = 6378388 * (1 - (1 / 297));
    double eCarre = (aCarre - Math.Pow(bLamb, 2)) / aCarre;
    double KLamb = 11565915.812935;

    double eLamb = Math.Sqrt(eCarre);
    double eSur2 = eLamb / 2;

    double Tan1 = (X - 150000.01256) / (5400088.4378 - Y);
    double Lambda = LongRef + (1 / nLamb) * (0.000142043 + Math.Atan(Tan1));
    double RLamb = Math.Sqrt(Math.Pow((X - 150000.01256), 2) + Math.Pow((5400088.4378 - Y), 2));

    double TanZDemi = Math.Pow((RLamb / KLamb), (1 / nLamb));
    double Lati1 = 2 * Math.Atan(TanZDemi);

    double eSin = 0;
    double Mult1 = 0;
    double Mult2 = 0;
    double Mult = 0;
    double LatiN = 0;
    double Diff = 0;

    double lat = 0;
    double lng = 0;

    do {
        eSin = eLamb * Math.Sin(Lati1);
        Mult1 = 1 - eSin;
        Mult2 = 1 + eSin;
        Mult = Math.Pow((Mult1 / Mult2), (eLamb / 2));
        LatiN = (Math.PI / 2) - (2 * (Math.Atan(TanZDemi * Mult)));
        Diff = LatiN - Lati1;
        Lati1 = LatiN;
    } while (Math.Abs(Diff) > 2.77777E-08);

    lat = (LatiN * 180) / Math.PI;
    lng = (Lambda * 180) / Math.PI;

    return new LatLon(lat, lng);
}

Answer 1

我是您在帖子中提到的页面的作者。
我不知道您是否解决了您的问题，但您提供的兰伯特坐标不正确。我认为你必须将它们除以 1000。这给出了 x=151488.250 和 y=170492.909，它们是可能的坐标，对应于......布鲁塞尔的一条街道。
在与 lat/lng 值之间转换时要小心选择基准。