<?php
function LatLonPointUTMtoLL($f, $f1, $j=32) {
    // $f = Northing
    // $f1 = Easting
    // $j: UTM ZONE (Danmark = 32)
    $d = 0.99960000000000004; // scale along long0
    $d1 = 6378137; // Polar Radius
    $d2 = 0.0066943799999999998;

    $d4 = (1 - sqrt(1 - $d2)) / (1 + sqrt(1 - $d2));
    $d15 = $f1 - 500000;
    $d16 = $f;
    $d11 = (($j - 1) * 6 - 180) + 3;
    $d3 = $d2 / (1 - $d2);
    $d10 = $d16 / $d;
    $d12 = $d10 / ($d1 * (1 - $d2 / 4 - (3 * $d2 * $d2) / 64 - (5 * pow($d2,3) ) / 256));
    $d14 = $d12 + ((3 * $d4) / 2 - (27 * pow($d4,3) ) / 32) * sin(2 * $d12) + ((21 * $d4 * $d4) / 16 - (55 * pow($d4,4) ) / 32) * sin(4 * $d12) + ((151 * pow($d4,3) ) / 96) * sin(6 * $d12);
    $d13 = rad2deg($d14);
    $d5 = $d1 / sqrt(1 - $d2 * sin($d14) * sin($d14));
    $d6 = tan($d14) * tan($d14);    
    $d7 = $d3 * cos($d14) * cos($d14);
    $d8 = ($d1 * (1 - $d2)) / pow(1 - $d2 * sin($d14) * sin($d14), 1.5);    
    $d9 = $d15 / ($d5 * $d);
    $d17 = $d14 - (($d5 * tan($d14)) / $d8) * ((($d9 * $d9) / 2 - (((5 + 3 * $d6 + 10 * $d7) - 4 * $d7 * $d7 - 9 * $d3) * pow($d9,4) ) / 24) + (((61 + 90 * $d6 + 298 * $d7 + 45 * $d6 * $d6) - 252 * $d3 - 3 * $d7 * $d7) * pow($d9,6) ) / 720);
    $d17 = rad2deg($d17); // Breddegrad (N)
    $d18 = (($d9 - ((1 + 2 * $d6 + $d7) * pow($d9,3) ) / 6) + (((((5 - 2 * $d7) + 28 * $d6) - 3 * $d7 * $d7) + 8 * $d3 + 24 * $d6 * $d6) * pow($d9,5) ) / 120) / cos($d14);
    $d18 = $d11 + rad2deg($d18); // Længdegrad (Ø)
    return array('lat'=>$d17,'long'=>$d18);
}

include_once "../tlib_functions.php";
html_header("Kode - Omregn northing/easting til latitude/longtitude for danske koordinater (UTM, zone 32)");

echo "<div id=\"form\">\n";
echo "<h1>Omregn northing/easting til længdegrad/breddegrad</h1>";
echo "<p>Eks.<br/>Northing = 6175755<br/>Easting = 724435<br/>(Rådhuspladsen 1, Kbh.)</p>";
echo "<form action=\"".$_SERVER['PHP_SELF']."\" method=\"post\">\n";
echo "Northing: <input type=\"text\" name=\"northing\" /><br/>\n";
echo "Easting:  <input type=\"text\" name=\"easting\" /><br/>\n";
echo "<input type=\"submit\" name=\"submit\" value=\"Omregn til længdegrad/breddegrad\" />\n";
echo "</form>\n";
echo "</div>\n";

if ( (int)$_POST['northing']>0 && (int)$_POST['easting']>0 ) {
    $res = LatLonPointUTMtoLL((int)$_POST['northing'],(int)$_POST['easting']);
    echo "<p>&nbsp;</p>\n";
    echo "<hr/>\n";
    echo "<p>&nbsp;</p>\n";
    echo "<h3>Resultat</h3>\n";
    echo "<pre>";
    echo "Northing: ".(int)$_POST['northing']."\n";
    echo " Easting: ".(int)$_POST['easting']."</pre>\n";

    echo "<pre>";
    echo "Breddegrad: ".$res['lat']."\n";
    echo "Længdegrad: ".$res['long']."</pre>\n";
    
    echo "<h3>Kopier til Google-maps</h3>\n";
    echo "<pre>new GLatLng(".$res['lat'].", ".$res['long'].")</pre>";
}
echo "<p><a href=\"omregn_northing_easting.phps\" target=\"_blank\">Vis kilde</a></p>";
html_footer();
?>