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加載離線瓦片參考例子

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM