[PYTHON] Laden Sie die Kacheln des Geographical Institute von Längen- und Breitengrad herunter

Das National Land Research Institute stellt Karten, Luftbilder usw. in Form von [Kachel des Geography Institute] zur Verfügung (https://maps.gsi.go.jp/development/ichiran.html). Um es zu erhalten, benötigen wir einen Koordinatenwert namens Kachelkoordinaten, der nicht Längen- und Breitengrad ist, also berechnen wir ihn. Ich habe [den Code im OpenStreetMap-Wiki] ausgeliehen (https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Python).

#from
#https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#Python

import math

#Berechnen Sie die Kachelkoordinaten aus Breiten- und Längengrad
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 * n)
    return (xtile, ytile)

#Ich werde es diesmal nicht verwenden, sondern von Kachelkoordinaten zu Breiten- und Längengraden
def num2deg(xtile, ytile, zoom):
    n = 2.0 ** zoom
    lon_deg = xtile / n * 360.0 - 180.0
    lat_rad = math.atan(math.sinh(math.pi * (1 - 2 * ytile / n)))
    lat_deg = math.degrees(lat_rad)
    return (lat_deg, lon_deg)

Das Herunterladen von Kacheln vom Geographical Institute sieht folgendermaßen aus

import requests
def download_from_gsi(z, x, y, def_url, fname):
    url = def_url.format(z,x,y)
    response = requests.get(url)
    if response.status_code == 200:
        image = response.content
        with open(fname, "wb") as f:
            f.write(image)
    else:
        raise Exception("{} returned {}".format(response.url, response.status_code))

Wie man es benutzt ist so. Dieses Beispiel ist eine Luftaufnahme der Koordinaten des Himmelsbaums.

#z ist die Zoomstufe.18 ist der am weitesten verbreitete Staat
z = 17
lat = 35.710163
lon = 139.8105428

def_url = "https://cyberjapandata.gsi.go.jp/xyz/ort/{}/{}/{}.jpg "
fname ="test.png "

x,y = deg2num(lat, lon, z)

download_from_gsi(z, x, y, def_url, fname)

Als Ergebnis habe ich ein Bild wie dieses bekommen. Es ist ein Bild von der Zeit, als es noch im Bau war.

Informationen zu def_url finden Sie in der URL unter Kachelliste des Geografischen Instituts. (Beachten Sie, dass es je nach Objekt jpg oder png sein kann.)