由于要在内网开发地图项目,不能访问在线的地图服务了,就想把地图瓦片下载下来,网上找了一些下载器都是需要注册及收费的,否则下载到的图都是打水印的,如下:
因为地图瓦片就是按照层级、行、列规则组织的一张张图片(如下图,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加载出来的效果和在线地图是一样的。