neo1973-germany/pylgrim/osmdump.py

65 lines
1.6 KiB
Python
Raw Normal View History

#!/usr/bin/python
#coding=utf8
import os
import sys
import math
import time
z = 10
def download(x,y,z):
import urllib
try:
webFile = urllib.urlopen("http://a.tile.openstreetmap.org/%d/%d/%d.png"%(z,x,y))
if not os.path.exists("%d"%z):
os.mkdir("%d"%z)
if not os.path.exists("%d/%d"%(z,x)):
os.mkdir("%d/%d"%(z,x))
localFile = open("%d/%d/%d.png"%(z,x,y), 'w')
localFile.write(webFile.read())
webFile.close()
localFile.close()
except Exception, e:
print e
def lon2km(lat):
return math.cos(lat*math.pi/180)*2*math.pi*6378.137/360
def getxy(lat,lon,z):
x = (lon+180)/360 * 2**z
y = (1-math.log(math.tan(lat*math.pi/180) + 1/math.cos(lat*math.pi/180))/math.pi)/2 * 2**z
return int(x),int(y)
lat = 49.009051
lon = 8.402481
r = 10
lat1 = lat-r/111.32
lon1 = lon-r/lon2km(49.009051)
lat2 = lat+r/111.32
lon2 = lon+r/lon2km(49.009051)
tiles = 0
#do not download zoom 18
for z in range(5,18):
x1,y1 = getxy(lat1, lon1, z)
x2,y2 = getxy(lat2, lon2, z)
tiles += (x2+1-x1)*(y1+1-y2)
print "do you really want to download %d tiles? [Y/n]"%tiles,
data = sys.stdin.read(1)
if data in ("y", "Y"):
i = 1;
for z in range(5,18):
x1,y1 = getxy(lat1, lon1, z)
x2,y2 = getxy(lat2, lon2, z)
for x in xrange(x1,x2+1):
for y in xrange(y2,y1+1):
if not os.path.exists("%d/%d/%d.png"%(z,x,y)):
download(x,y,z)
print "\r%i"%i,
sys.stdout.flush()
i+=1