Gauss-Boaga to WGS84

Question

我正在尝试将一些坐标从Gauss-Boaga转换为WGS84。下面是我所做的(PHP):

function gaussToLatLng($Eutm, $Nutm){
    // parametri che dipendono dal fuso ovest
    $l0 = 9;
    $fraz = $Nutm/111092.0821;
    echo "$fraz ";
    $A= $fraz + 
        (0.1449300705 * sin(deg2rad(2*$fraz))) + 
        (0.0002138508 * sin(deg2rad(4*$fraz))) + 
        (0.0000004322 * sin(deg2rad(6*$fraz)));
    $v = sqrt(1 + (0.0067681702 * cos(deg2rad($A)) * cos(deg2rad($A)) ) );
    $y = $Eutm - 500000;
    $B = rad2deg(atan( ($v * sinh(deg2rad($y/6397376.633)) ) / cos(deg2rad($A))));
    $lng = rad2deg(atan(tan(deg2rad($A)) * cos(deg2rad($v * $B))));
    $lat = $B + $l0;
    echo "A=$A, B=$B n";
    return array(
        'lat' => $lat,
        'lng' => $lng
    );
}

使用(1517140,5036970)作为输入，我得到(9.2271558768758 45.485183518206)，而使用这个在线工具，我计算出我应该有(9.2189597,45.4859253)(大约1公里的差异)。

你能帮我找出代码中的错误吗?

Answer 1

使用PHP库解决:proj4php，使用此处找到的Gauss Boaga fuso oest的定义。

那些家伙太棒了!