GIS Coordinate transformation #2

GIS 座標変換の計算でハマるが自己解決

Unfortunately this text is in Japanese only.
座標からタイルのx,y座標を求めるがハマった件。
OpenStreetMapの [Wiki : Slippy map tilenames] を参照して真球メルカトル投影の座標変換を行ったつもりが。。。

OpenStreetMap Wikiから転写。

表計算ソフトでうまくいかないのでOpenStreetMap WikiのVB.NETのソースをそのまま使用して計算してみると．．．

== ======= ====== ZOOM {X} {Y} == ======= ====== 0 0 0 1 1 0 2 3 1 3 6 3 4 13 6 5 27 12 6 55 25 7 110 51 8 221 102 9 442 204 10 884 409 11 1,768 819 12 3,536 1,639 13 7,072 3,279 14 14,145 6,559 15 28,290 13,118 16 56,580 26,236 17 113,160 52,472 == ======= ======

何の事はない、ちゃんと計算できる(T_T
試しにZoomレベル１７の場合

OpenStreetMap
URLは
http://a.tile.openstreetmap.org/17/113160/52472.png

実際に取得するなら [ここをクリック]

電子国土
URLは
http://cyberjapandata.gsi.go.jp/xyz/std/17/113160/52472.png

実際に取得するなら [ここをクリック]

これでイロイロなソフトウェアが組めるゾ♪
でも、何で表計算ソフトではうまくいかない？？？








2014/01/05 追記
原因判明。
表計算ソフトのLOG()は基数省略すると底は１０だと．．．
e を底とする自然対数はLN()
うまく計算できました。