josch/mapbender
1
0
Fork 0
Code Issues Pull requests Projects Releases Wiki Activity

disable plotting and add hint how to generate tiles

Browse source
This commit is contained in:
josch 2014-06-21 14:08:12 +02:00
parent 6220549097
commit ddf598c8d9
1 changed files with 17 additions and 15 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

32
mapbender.py
View file

@ -280,6 +280,8 @@ def main(x,y,width,smoothing,subdiv):
# v = off+t[0]*h # v = off+t[0]*h
# sx.append(u) # sx.append(u)
# sy.append(v) # sy.append(v)
# create map with
# python -c 'import logging; logging.basicConfig(level=logging.DEBUG); from landez import ImageExporter; ie = ImageExporter(tiles_url="http://{s}.tile.opencyclemap.org/cycle/{z}/{x}/{y}.png"); ie.export_image(bbox=(8.0419921875,51.25160146817652,10.074462890625,54.03681240523652), zoomlevel=14, imagepath="image.png")'
im = Image.open("map.png") im = Image.open("map.png")
bbox = [8.0419921875,51.25160146817652,10.074462890625,54.03681240523652] bbox = [8.0419921875,51.25160146817652,10.074462890625,54.03681240523652]
# apply mercator projection # apply mercator projection
@ -314,21 +316,21 @@ def main(x,y,width,smoothing,subdiv):
data.append((box,quad)) data.append((box,quad))
im_out = im.transform((int(iw*width/(bbox[2]-bbox[0])),int(ih*height/(bbox[3]-bbox[1]))),Image.MESH,data,Image.BICUBIC) im_out = im.transform((int(iw*width/(bbox[2]-bbox[0])),int(ih*height/(bbox[3]-bbox[1]))),Image.MESH,data,Image.BICUBIC)
im_out.save("out.png") im_out.save("out.png")
np.random.seed(seed=0) #np.random.seed(seed=0)
colors = 100*np.random.rand(len(patches)/2)+100*np.random.rand(len(patches)/2) #colors = 100*np.random.rand(len(patches)/2)+100*np.random.rand(len(patches)/2)
p = PatchCollection(patches, cmap=matplotlib.cm.jet, alpha=0.4) #p = PatchCollection(patches, cmap=matplotlib.cm.jet, alpha=0.4)
p.set_array(np.array(colors)) #p.set_array(np.array(colors))
plt.figure() #plt.figure()
plt.axes().set_aspect('equal') #plt.axes().set_aspect('equal')
#plt.axhspan(0, height, xmin=0, xmax=width) ##plt.axhspan(0, height, xmin=0, xmax=width)
fig, ax = plt.subplots() #fig, ax = plt.subplots()
#ax.add_collection(p) ##ax.add_collection(p)
ax.set_aspect('equal') #ax.set_aspect('equal')
plt.axis((0,width,0,height)) #plt.axis((0,width,0,height))
plt.imshow(np.asarray(im_out),extent=[0,width,0,height]) #plt.imshow(np.asarray(im_out),extent=[0,width,0,height])
plt.imshow(np.asarray(im),extent=[bbox[0],bbox[2],bbox[1],bbox[3]]) #plt.imshow(np.asarray(im),extent=[bbox[0],bbox[2],bbox[1],bbox[3]])
plt.plot(x,y,out[0],out[1],px,py,qx,qy,tx,ty) #plt.plot(x,y,out[0],out[1],px,py,qx,qy,tx,ty)
plt.show() #plt.show()
return True return True
if __name__ == '__main__': if __name__ == '__main__':