Diary/2012-4-21の変更点

• 追加された行はこのように表示されます。
• 削除された行はこのように表示されます。

!!つぎに"測地系を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･Perlhttp://homepage3.nifty.com/Nowral/index.html] ( http://homepage3.nifty.com/Nowral/index.html )
*[Mac･GPS･Perl|http://homepage3.nifty.com/Nowral/index.html] ( http://homepage3.nifty.com/Nowral/index.html )