!!つぎに"測地系をTokyoからWGS84に"とその逆をPerlで既述{{category GPSネタ}} なるべく元ねたオリジナルの既述を踏襲することでつぎのようにした # 測地系変換 TokyoをWGS84 # 呼出 &T2W ( 緯度,経度 ) ; 実数,実数 # 戻値 ( 緯度,経度 ) ; 実数,実数 # $BWGS84 = $BTokyo - 0.00010695*$BTokyo + 0.000017464*$LTokyo + 0.0046017 # $LWGS84 = $LTokyo - 0.000046038*$BTokyo - 0.000083043*$LTokyo + 0.010040 sub T2W { my($BTokyo,$LTokyo) = @_; $BWGS84 = $BTokyo - 0.00010695*$BTokyo + 0.000017464*$LTokyo + 0.0046017; $LWGS84 = $LTokyo - 0.000046038*$BTokyo - 0.000083043*$LTokyo + 0.010040; return ($BWGS84,$LWGS84); } # 測地系変換 WGS84をTokyo # 呼出 &W2T ( 緯度,経度 ) ; 実数,実数 # 戻値 ( 緯度,経度 ) ; 実数,実数 # $BTokyo = $BWGS84 + 0.00010696*$BWGS84 - 0.000017467*$LWGS84 - 0.0046020 # $LTokyo = $LWGS84 + 0.000046047*$BWGS84 + 0.000083049*$LWGS84 - 0.010041 sub T2W { my($BWGS84,$LWGS84) = @_; $BWGS84 = $BTokyo - 0.00010695*$BTokyo + 0.000017464*$LTokyo + 0.0046017; $LWGS84 = $LTokyo - 0.000046038*$BTokyo - 0.000083043*$LTokyo + 0.010040; return ($BTokyo,$LTokyo); } !元ネタ *[Mac・GPS・Perl|http://homepage3.nifty.com/Nowral/index.html] ( http://homepage3.nifty.com/Nowral/index.html )