![[PukiWiki]](image/picc.png "[PukiWiki]")
地域メッシュコード
4分の1地域メッシュのコードから緯度経度を算出することになった.
簡単にExcelで出せるように関数がネット上によく転がっていた.
セル(A1)に10桁のコードを置いて、グリッドの境界部分は下記の緯度(南北)経度(東西)になる
南
=(LEFT(A1,2)/1.5*3600+MID(A1,5,1)*5*60+MID(A1,7,1)*30 +
IF(MID(A1,9,1)*1<>0,IF(OR(MID(A1,9,1)*1=3,MID(A1,9,1)*1=4),15,0)) +
IF(MID(A1,10,1)*1<>0,IF(OR(MID(A1,10,1)*1=3,MID(A1,10,1)*1=4),7.5,0))) /3600
北
=(LEFT(A1,2)/1.5*3600+MID(A1,5,1)*5*60+MID(A1,7,1)*30+
IF(MID(A1,9,1)*1<>0,IF(OR(MID(A1,9,1)*1=3,MID(A1,9,1)*1=4),15,0)) +
IF(MID(A1,10,1)*1<>0,IF(OR(MID(A1,10,1)*1=3,MID(A1,10,1)*1=4),15,7.5))) /3600
西
=((MID(A1,3,2)+100)*3600+MID(A1,6,1)*7.5*60+MID(A1,8,1)*45 +
IF(MID(A1,9,1)*1<>0,IF(OR(MID(A1,9,1)*1=2,MID(A1,9,1)*1=4),22.5,0)) +
IF(MID(A1,10,1)*1<>0,IF(OR(MID(A1,10,1)*1=2,MID(A1,10,1)*1=4),11.25,0))) /3600
東
=((MID(A1,3,2)+100)*3600+MID(A1,6,1)*7.5*60+MID(A1,8,1)*45 +
IF(MID(A1,9,1)*1<>0,IF(OR(MID(A1,9,1)*1=2,MID(A1,9,1)*1=4),22.5,0)) +
IF(MID(A1,10,1)*1<>0,IF(OR(MID(A1,10,1)*1=2,MID(A1,10,1)*1=4),22.5,11.25))) /3600
Google先生に確認をしてみた
<?xml version='1.0' encoding='UTF-8'?>
<kml xmlns='http://earth.google.com/kml/2.0'> <Document>
<Placemark>
<Polygon> <outerBoundaryIs> <LinearRing>
<coordinates>
139.7625,35.68125,0.
139.765625,35.68125,0.
139.765625,35.67916667,0.
139.7625,35.67916667,0.
</coordinates>
</LinearRing> </outerBoundaryIs> </Polygon>
<Style>
<PolyStyle>
<color>#900000ff</color>
<outline>0</outline>
</PolyStyle>
</Style>
</Placemark>
</Document> </kml>
残念なことに縦横250m位が本当なのだが、南北230m、東西280mほどな長方形である....
