地形图拼接方法初探

本人从事城市规划行业，用Python干过拼接地形图的事。发现python很好用，的确可以解决工作中的很多问题。

本人非大牛，代码能解决问题就OK，代码不一定最优，有优化的方法，欢迎去知乎。

代码仅介绍方法，非最终版本。

干过城市规划的人应该知道，地形图是我们进行踏勘、绘制用地现状图的基础文件。特别是做城市总体规划的时候，往往会获得整个城市的地形图，多的上万，少的上千个图幅文件。而我们知道，各个区域的行政边界（县、街道、镇、乡、村）是不规则的，跑现场的时候都是分区域跑的，因此就存在分区域打印地形图的需求。当然，你也可以整个区域一起打印，但问题是地物的精度就不够了。问题来了，几千个文件，我怎么知道哪个文件属于哪个区域呢？这涉及到空间上的对应关系，好在正规的地形图绘制单位都会做地形图的拼接表，上面有数字与地形图文件编号对应，这就成为我们后面筛选特定区域文件的关键。

我的解决思路是（1）利用ArcGIS实现空间上特定区域图幅编号的筛选，（2）然后根据编号分文件夹提取地形图文件，（3）最后进行拼接。