Python爬取地图瓦片

由于要在内网开发地图项目,不能访问在线的地图服务了,就想把地图瓦片下载下来,网上找了一些下载器都是需要注册及收费的,否则下载到的图都是打水印的,如下:

因为地图瓦片就是按照层级、行、列规则组织的一张张图片(如下图,14是瓦片的级别,13478是行号,6202是列号),所以想到可以试试用python把这些图爬取下来。

主要可以分为三步:

1.确定要下载的瓦片的级别,经度和纬度范围(可以通过百度坐标拾取获得,http://api.map.baidu.com/lbsapi/getpoint/index.html)

2.计算出这个范围内瓦片的起始和终止行列号

1 # 经纬度反算切片行列号 3857坐标系
2 def deg2num(lat_deg, lon_deg, zoom):
3     lat_rad = math.radians(lat_deg)
4     n = 2.0 ** zoom
5     xtile = int((lon_deg + 180.0) / 360.0 * n)
6     ytile = int((1.0 - math.log(math.tan(lat_rad) + (1 / math.cos(lat_rad))) / math.pi) / 2.0 * n)
7     return (xtile, ytile)

3.根据行列号拼接出瓦片的url地址

1 #Google地图瓦片为例,zoom是地图层级,x,y为行列号
2 tilepath = 'http://www.google.cn/maps/vt/pb=!1m4!1m3!1i'+str(zoom)+'!2i'+str(x)+'!3i'+str(y)+'!2m3!1e0!2sm!3i345013117!3m8!2szh-CN!3scn!5e1105!12m4!1e68!2m2!1sset!2sRoadmap!4e0'

下载图片

 1 # 下载图片
 2 def getimg(Tpath, Spath, x, y):
 3     try:
 4         f = open(Spath, 'wb')
 5         req = urllib.request.Request(Tpath)
 6         req.add_header('User-Agent', random.choice(agents))  # 换用随机的请求头
 7         pic = urllib.request.urlopen(req, timeout=60)
 8         f.write(pic.read())
 9         f.close()
10         print(str(x) + '_' + str(y) + '下载成功')
11     except Exception:
12         print(str(x) + '_' + str(y) + '下载失败,重试')
13         getimg(Tpath, Spath, x, y)

保存图片时也需要按照层级,行创建好目录存储

1 path = "D:\\GoogleImg\\" + str(zoom) + "\\" + str(x)
2     if not os.path.exists(path):
3         os.makedirs(path)

结果如下:

用Openlayers加载出来的效果和在线地图是一样的。

Python获取地图瓦片的完整代码

Openlayers加载离线瓦片参考例子

Leaflet加载离线瓦片参考例子

原文地址:https://www.cnblogs.com/jyughynj/p/pythonGetMapTile.html