使用 Python 合并地图瓦片

港控/mmm° 2022-02-28 14:08 947阅读 0赞

2018 年 04 月 06 日 18:34 本文 1946 字，阅读全文约需 6 分钟  
前文提到了合并瓦片图，而瓦片图应用比较多的则是瓦片地图。对地图本就感兴趣的我，也想试试合并互联网地图的某个范围内的地图图层。

随着技术的发展，国内的地图服务商相继将地图瓦片更新为矢量瓦片\[1\]，这下想取到相应的瓦片图，我感觉没那么简单。那还是找个相对简单的，还使用栅格瓦片的地图服务网站来试试。

我先想到的是国字号的“天地图”，天地图·广东\[2\]还在使用栅格瓦片，其地图放大到一定级别时，可以看到水域及建筑物图层等。

广东的天地图比起江浙等省份的还是差一大截，精确度也有待提高。不过观察对比后发现，澄海这边的图层也还是要比高德、百度等地图要细腻一些，如图是莲下镇的一部分，有些大巷子都能体现。

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MDk3NjI2MQ_size_16_color_FFFFFF_t_70]  
莲下部分区域▲ 莲下部分区域  
查看请求可知其瓦片图的规则如下：

[http://services.tianditugd.com/gdvec2014/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=GDVEC2014&STYLE=GDVEC2014&TILEMATRIXSET=Matrix\_0&TILEMATRIX=17&TILEROW=\{row\}&TILECOL=\{col\}&FORMAT=image%2Ftile][http_services.tianditugd.com_gdvec2014_wmts_SERVICE_WMTS_REQUEST_GetTile_VERSION_1.0.0_LAYER_GDVEC2014_STYLE_GDVEC2014_TILEMATRIXSET_Matrix_0_TILEMATRIX_17_TILEROW_row_TILECOL_col_FORMAT_image_2Ftile]  
其中 \{row\} 和 \{col\}，分别为行和列，也就可以说对应的是 y 和 x。因按链接下载到的文件名不规则，故每张图需另存为有规则的文件名，aria2 很好解决，只需加个参数即可。下面取的是接近澄海区域内，放大级别为 17 的地图瓦片。

\#!/usr/bin/env python  
output =’’  
r = 24170  
while (r <= 24255):  
c = 108020  
while (c <= 108100 ):  
output += ‘[http://services.tianditugd.com/gdvec2014/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=GDVEC2014&STYLE=GDVEC2014&TILEMATRIXSET=Matrix\_0&TILEMATRIX=17&TILEROW=’+str(r)+’&TILECOL=’+str(c)+’&FORMAT=image%2Ftile\\n][http_services.tianditugd.com_gdvec2014_wmts_SERVICE_WMTS_REQUEST_GetTile_VERSION_1.0.0_LAYER_GDVEC2014_STYLE_GDVEC2014_TILEMATRIXSET_Matrix_0_TILEMATRIX_17_TILEROW_str_r_TILECOL_str_c_FORMAT_image_2Ftile_n] out=’+str©+’\_’+str®+’.jpg\\n’  
c += 1  
r += 1

file = open(‘tiles.txt’, ‘w’)  
file.write(output)  
file.close()  
合并时，和上篇差不多，或许是数字比较大，到后面排序乱了，再进行一次排序后正常。代码如下：

\#!/usr/bin/env python  
import glob  
import re  
from PIL import Image

files = glob.glob(’./*\_*.jpg’)  
files.sort(key=lambda x: tuple(int(i) for i in re.findall(’\\d+’, x)\[:2\]))

imagefiles = \{\}  
for item in files:  
match = re.search(r’\\d+’, item)  
pre = int(match.group())  
if not imagefiles.get(pre):  
imagefiles\[pre\] = \[\]  
imagefiles\[pre\].append(item)

imagefiles = sorted(zip(imagefiles.keys(),imagefiles.values()));

total\_width = len(imagefiles) \* 256  
total\_height = len(imagefiles\[0\]\[1\]) \* 256

new\_image = Image.new(‘RGB’, (total\_width, total\_height))

x\_offset = 0  
for item in imagefiles:  
y\_offset = 0  
images = map(Image.open, item\[1\])  
for subitem in images:  
new\_image.paste(subitem, (x\_offset, y\_offset))  
y\_offset += subitem.size\[0\]  
x\_offset += images\[0\].size\[0\]

new\_image.save(‘merge.jpg’, quality = 90)  
下面是经缩小后的地图，原图长宽各两万多像素， 40 多兆。

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MDk3NjI2MQ_size_16_color_FFFFFF_t_70 1]  
包含建筑物、水域等图层的澄海地图底图▲ 包含建筑物、水域等图层的澄海地图底图  
对了，它的底图跟地名图层是分开的，且并非所有底图瓦片都有对应的地名瓦片，等哪时闲了再折腾。

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MDk3NjI2MQ_size_16_color_FFFFFF_t_70]: /images/20220228/a6e5c335242b4a1e8c272343c92a20d6.png
[http_services.tianditugd.com_gdvec2014_wmts_SERVICE_WMTS_REQUEST_GetTile_VERSION_1.0.0_LAYER_GDVEC2014_STYLE_GDVEC2014_TILEMATRIXSET_Matrix_0_TILEMATRIX_17_TILEROW_row_TILECOL_col_FORMAT_image_2Ftile]: http://services.tianditugd.com/gdvec2014/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=GDVEC2014&STYLE=GDVEC2014&TILEMATRIXSET=Matrix_0&TILEMATRIX=17&TILEROW=%7Brow%7D&TILECOL=%7Bcol%7D&FORMAT=image%2Ftile
[http_services.tianditugd.com_gdvec2014_wmts_SERVICE_WMTS_REQUEST_GetTile_VERSION_1.0.0_LAYER_GDVEC2014_STYLE_GDVEC2014_TILEMATRIXSET_Matrix_0_TILEMATRIX_17_TILEROW_str_r_TILECOL_str_c_FORMAT_image_2Ftile_n]: http://services.tianditugd.com/gdvec2014/wmts?SERVICE=WMTS&REQUEST=GetTile&VERSION=1.0.0&LAYER=GDVEC2014&STYLE=GDVEC2014&TILEMATRIXSET=Matrix_0&TILEMATRIX=17&TILEROW='%20rel=
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MDk3NjI2MQ_size_16_color_FFFFFF_t_70 1]: /images/20220228/f4d3b78eb9294d1cafc7cd585a2fe4cb.png