緯度経度から地理院タイルをダウンロードする - Qiita

国土地理院は地図や航空写真等を地理院タイルという形で提供している 取得するにはタイル座標という緯度経度ではない座標値が必要なので，それを計算する． OpenStreetMapのwikiにそのコードがあるので拝借． #from #https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Python import math #緯度経度からタイル座標を計算 def deg2num(lat_deg, lon_deg, zoom): lat_rad = math.radians(lat_deg) n = 2.0 ** zoom xtile = int((lon_deg + 180.0) / 360.0 * n) ytile = int((1.0 - math.asinh(math.tan(lat_rad)) / math.pi) / 2.0