WMTS：特定のズームレベルで、位置情報（緯度、経度）をタイルインデックスに変換しますか？

特定の地理位置（緯度、経度）およびズームレベルのWMTSタイルのインデックス（x、y）を取得する方法を知りたいと思いました。

たとえば、（48.675、2.7）にPOIがあります。ズーム10に対応するオープンストリートマップタイルを取得します。

数学はできますか？Webサービスが必要ですか？精度：これをプログラムで行う必要があります。

OSM wikiページは完璧です：http : //wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Lon..2Flat._to_tile_numbers_2

クイックリファレンス用に抽出された回答（擬似コード）を以下に示します。

経度/緯度/ズームがタイル番号に与えられた場合：

n = 2 ^ zoom
xtile = n * ((lon_deg + 180) / 360)
ytile = n * (1 - (log(tan(lat_rad) + sec(lat_rad)) / π)) / 2

log()この疑似コードでは、自然対数（一般的な数学構文ではしばしば「ln（）」ですが、多くのプログラミング言語ではしばしば「log（）」）を指すことに注意してください。

タイル番号を経度/緯度に与えます：

n = 2 ^ zoom
lon_deg = xtile / n * 360.0 - 180.0
lat_rad = arctan(sinh(π * (1 - 2 * ytile / n)))
lat_deg = lat_rad * 180.0 / π