轻松显示适用于Python的OpenStreetMap磁贴

Question

我想在我的python代码中包含开放式街道地图（OSM）。

我已经阅读了很多有关OSM的网页。但是不幸的是，我对哪种包装最满意感到迷茫。

我正在寻找一种在我的应用中获取OSM映像的简便方法。当我开始时，我想到的是：

import matplotlib.pyplot as plt

# Pseudo - Code for required function 'GetOSMImage'
Map = GetOSMImage(lat,long,delta_lat,delta_long)

imgplot = plt.imshow(Map)

稍后，我想在此plt中添加绘图我的其他数据。（我知道我需要处理预测等）

我不需要/想要的：

• 在我自己的网站上显示
• 要将我的数据上传到某些Internet服务器
• 互动功能，例如缩放，滚动（首先）
• 从OSM手动处理和呈现.xml数据
• 首先，我不想定义渲染样式的每个细节。我希望/期望存在一些默认样式。

你对我有一个很好的起点吗？还是我低估了这个主题的复杂性？

Answer 1

根据您的输入，我能够实现目标。这是我给其他人的代码，它们正在寻找OSM的起点。（当然，还有很多改进的余地）。

更新资料

请遵守开放街道地图的使用政策！

OpenStreetMap数据可供所有人免费使用。 我们的磁贴服务器不是。

要求

• 未经操作工作组的事先许可，禁止大量使用（例如，分发使用openstreetmap.org中的图块的应用程序）。请参阅下面的替代方法。
• 清楚显示许可证归属。
• 不要主动或被动地鼓励侵犯版权。
• 只能直接通过最终用户操作来触发对/ cgi-bin / export的调用。（例如：“单击此处进行导出”。）导出调用是运行时非常昂贵的（CPU
• RAM）函数，在服务器处于高负载时会经常拒绝。
• 推荐：请勿在tile.openstreetmap.org上对任何URL进行硬编码，因为这样做会限制您在服务中断或阻止时快速做出反应的能力。
  * 推荐：添加一个指向https://www.openstreetmap.org/fixthemap的链接，以允许您的用户报告并解决我们数据中的问题。
  
  
  

技术使用要求

• 有效的HTTP User-Agent识别应用程序。伪造另一个应用程序的User-Agent会让您被阻止。
• 如果已知，则为有效的HTTP引用。
• 不要发送无缓存头。（“缓存控制：不缓存”，“编译指示：不缓存”等）
• 缓存切片根据HTTP Expiry标头在本地下载，或者至少7天。
• 最多2个下载线程。（可接受未修改的Web浏览器的下载线程限制。）
  
  

有关更多详细信息，请参见：https :
//operations.osmfoundation.org/policies/tiles/

这是代码：

import matplotlib.pyplot as plt
import numpy as np

import math
import urllib2
import StringIO
from PIL import Image



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.log(math.tan(lat_rad) + (1 / math.cos(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)



def getImageCluster(lat_deg, lon_deg, delta_lat,  delta_long, zoom):
    smurl = r"http://a.tile.openstreetmap.org/{0}/{1}/{2}.png"
    xmin, ymax =deg2num(lat_deg, lon_deg, zoom)
    xmax, ymin =deg2num(lat_deg + delta_lat, lon_deg + delta_long, zoom)

    Cluster = Image.new('RGB',((xmax-xmin+1)*256-1,(ymax-ymin+1)*256-1) ) 
    for xtile in range(xmin, xmax+1):
        for ytile in range(ymin,  ymax+1):
            try:
                imgurl=smurl.format(zoom, xtile, ytile)
                print("Opening: " + imgurl)
                imgstr = urllib2.urlopen(imgurl).read()
                tile = Image.open(StringIO.StringIO(imgstr))
                Cluster.paste(tile, box=((xtile-xmin)*256 ,  (ytile-ymin)*255))
            except: 
                print("Couldn't download image")
                tile = None

    return Cluster



if __name__ == '__main__':

    a = getImageCluster(38.5, -77.04, 0.02,  0.05, 13)
    fig = plt.figure()
    fig.patch.set_facecolor('white')
    plt.imshow(np.asarray(a))
    plt.show()