# coding: utf-8
 
def LatLng2FirstMeshCode( lat, lng ):
if ( ( lat < 20.0 ) or ( lat >= 46.0 ) ):
print "エラー：緯度が範囲を超えています。"; return;
if ( ( lng < 122.0 ) or ( lng >= 154.0 ) ):
print "エラー：経度が範囲を超えています。"; return;
latd = int( lat * 1.5 );
lngd = int( lng - 100.0 );
print "%d%d" % ( latd, lngd )
 
def FirstMeshCode2LatLng( code ):
latd = int( code[0:2] );
lngd = int( code[2:4] );
if ( ( latd < 30 ) or ( latd >= 69 ) ):
print "エラー：緯度が範囲を超えています。"; return;
if ( ( lngd < 22 ) or ( lngd >= 54 ) ):
print "エラー：経度が範囲を超えています。"; return;
lat = latd / 1.5;
lng = lngd + 100.0;
print "(%f,%f)" % ( lat, lng )
 
LatLng2FirstMeshCode( 34.1, 136.1 );
FirstMeshCode2LatLng( "5136" )
 
LatLng2FirstMeshCode( 19.1, 136.1 );
FirstMeshCode2LatLng( "5160" )

# coding: utf-8
 
def LatLng2SecondMeshCode( lat, lng ):
# 範囲検査
if ( ( lat < 20.0 ) or ( lat >= 46.0 ) ):
print "エラー：緯度が範囲を超えています。"; return;
if ( ( lng < 122.0 ) or ( lng >= 154.0 ) ):
print "エラー：経度が範囲を超えています。"; return;
# 緯度・経度を秒単位の整数で表す。
latd = int( lat );
tmp  = ( lat - latd ) * 60.0;
latm = int( tmp );
lats = int( int( ( tmp - latm ) * 60.0 * 1000.0 ) / 1000.0 );
lngd = int( lng );
tmp  = ( lng - lngd ) * 60.0;
lngm = int( tmp );
lngs = int( int( ( tmp - lngm ) * 60.0 * 1000.0 ) / 1000.0 );
I = latd * 3600 + latm * 60 + lats;
K = lngd * 3600 + lngm * 60 + lngs;
# 2次メッシュコード表示
print "%d%d%d%d" % ( (I/2400)%100, (K/3600)%100, \
(I%2400)/300, (K%3600)/450 )
 
def SecondMeshCode2LatLng( code ):
# 2次メッシュコード文字列分解
lat1 = int( code[0:2] );
lng1 = int( code[2:4] );
lat2 = int( code[4:5] );
lng2 = int( code[5:6] );
# 範囲検査
if ( ( lat1 < 30 ) or ( lat1 >= 69 ) ):
print "エラー：緯度が範囲を超えています。"; return;
if ( ( lng1 < 22 ) or ( lng1 >= 54 ) ):
print "エラー：経度が範囲を超えています。"; return;
# 緯度・経度を秒単位の整数で表す。
I = lat1 * 2400 + lat2 * 300;
K = ( lng1 + 100 ) * 3600 + lng2 * 450;
# 度単位に変換し、表示する。
lat = I / 3600.0;
lng = K / 3600.0;
print "(%f,%f)" % ( lat, lng )
 
LatLng2SecondMeshCode( 34.345, 136.234 )
SecondMeshCode2LatLng( "513641" )