Diary/2012-4-21

 つぎに"測地系をTokyoからWGS84に"とその逆をPerlで既述[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 )