위도 경도에서 지리원 타일을 다운로드하다
7842 단어 Python
좌표 좌표 좌표의 위도 경도가 아닌 좌표 값을 얻다.
네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 * n)
return (xtile, ytile)
#今回は使わないけど,タイル座標から緯度経度
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)
지리원 타일에서 다운 받으면 이런 느낌이에요.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))
사용법은 이런 느낌이다. 이 예는 하늘나무의 좌표의 항공 사진이다.#zはズームレベル.18が一番拡大した状態
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)
결과는 이런 이미지를 얻었다.건설 도중 사진인가?
def_url지리원 타일 일람의 url을 참고하세요.
(주의사항은 jpg일 때도 있고 pg일 때도 있다.)
Reference
이 문제에 관하여(위도 경도에서 지리원 타일을 다운로드하다), 우리는 이곳에서 더 많은 자료를 발견하고 링크를 클릭하여 보았다 https://qiita.com/J-Taniguchi/items/e6d037ff938355c8d69e텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
우수한 개발자 콘텐츠 발견에 전념 (Collection and Share based on the CC Protocol.)