如何绘制带坐标的png-osm-map

时间:2012-05-16 12:36:05

标签: python matplotlib gis openstreetmap

我想在Python中创建一个包含几个给定点的地图。为此,我想使用matplotlib中的Basemap。它效果很好,但我不知道如何获得合适的背景图。

如何导入OSM地图?或者我应该使用不同的映射包?我只是想创建一个栅格地图并将其保存为png。

2 个答案:

答案 0 :(得分:3)

这不是我的解决方案;我已经从问题中粘贴了它,因为提问者没有足够的声誉来回答他自己的问题。

  

我找到了解决方案:

     

在底图中使用imshow包含一个png到图中   背景图片。为了获得正确的背景图像,我使用了   export feature of OSM的边界取自Basemap   构造:

m = Basemap(llcrnrlon=7.4319, urcrnrlat=52.0632, urcrnrlon=7.848, llcrnrlat=51.8495,
        resolution='h', projection='merc')

im = plt.imread('background.png')
m.imshow(im, interpolation='lanczos', origin='upper')

答案 1 :(得分:2)

我从NASA GIBS tileserver中找到了一些可访问的底图图像。您可以对其他tileservers使用相同的方法。

http://earthdata.nasa.gov/wiki/main/index.php/GIBS_Supported_Clients#Script-level_access_to_imagery

在python子shell中使用GDAL的gdal_translate:

import subprocess
import matplotlib.pyplot
import mpl_toolkits.basemap

l,u,r,d=(7.4319,52.0632,7.848,51.8495)

subprocess.call ('gdal_translate -of GTiff -outsize 400 400 -projwin {l} {u} {r} {d} TERRA.xml Background.tif'.format(l=l,u=u,r=r,d=d),shell=True )

im=matplotlib.pyplot.imread('Background.tif')

m = mpl_toolkits.basemap.Basemap(llcrnrlon=l, urcrnrlat=u, urcrnrlon=r, llcrnrlat=d,
    resolution='h', projection='merc')
m.imshow(im, interpolation='lanczos', origin='upper')

matplotlib.pyplot.show()

这需要上面链接中的TERRA.xml文件,尽管你也可以内联XML。