function Demo_TW() { $p = array(282745, 2641869); print_r($p); echo "cs2cs法:\n"; $p1 = cs2cs_t67to97($p[0], $p[1], 0); print_r($p1); echo "cs2cs法 invert:\n"; $pp = cs2cs_t97to67($p1[0], $p1[1], 0); print_r($pp); echo "平面四參數轉換:\n"; $p1 = t67to97($p[0], $p[1]); print_r($p1); echo "平面四參數轉換 Invert:\n"; $p2 = t97to67($p1[0], $p1[1]); print_r($p2); echo "使用平面四參數轉換:\n"; $q = proj_67toge($p); print_r($q); echo "使用 cs2cs:\n"; $q1 = proj_67toge2($p); print_r($q); }
function proj_geto67($p) { $pp = array(deg2dms($p[0]) . "E", deg2dms($p[1]) . "N"); // print_r($pp); $proj = "proj +proj=tmerc +ellps=aust_SA +lon_0=121 +x_0=250000 +k=0.9999"; $k = addslashes("{$pp['0']} {$pp['1']}"); $cmd = "echo {$k} | {$proj}"; // echo $cmd ."\n"; $ret = shell_exec("echo {$k} | {$proj}"); list($x, $y) = preg_split("/\\s+/", $ret); return t97to67($x, $y); }