由于要在内网开发地图项目,不能访问在线的地图服务了,就想把地图瓦片下载下来,网上找了一些下载器都是需要注册及收费的,否则下载到的图都是打水印的,如下:
因为地图瓦片就是按照层级、行、列规则组织的一张张图片(如下图,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加载离线瓦片参考例子
如果对您有用,请赞赏!谢谢!