利用Python自动化生成地质解释图件

张伟，邱崇涛，谢明宏，赵丛

(1.核工业航测遥感中心，河北 石家庄 050002；2.中核集团 铀资源地球物理勘查技术中心(重点实验室)，河北 石家庄 050002)

0 引言

在地质、物探工作中，自始至终都会应用不同形式的专题地图来直观地表现其研究成果[1]。为获取最佳的地质效果，技术人员需不断地对地质情况深入研究，反复斟酌，经多次易稿后才能获取最终满意的成果，而伴随其中的专题图件也在不断的修改中完善。

在物探工作中，地质推断解释图是反映物探成果的必备专题图件之一。总体而言，地质推断解释图中的图元要素并非十分复杂；但由于修改次数较多，又受到专业制图软件功能的限制，工作效率难以满足快节奏项目需求。以目前国内应用范围最广的国产软件MapGIS为例，虽有强大编辑功能，但凭手工绘图方式，其效率依然难有明显提升。如电法勘探中的反演电阻率断面的地质推断解释图，不整合线是一种常见的地质现象，图中使用波浪线表示，而在MapGIS软件中波浪线则由一个特殊线型呈现，实际坐标位置处于曲线的中轴线上；若直接成区，线与区间的分界面无法满足制图要求，仍需徒手绘制或借助于第三方软件生成，操作繁琐且费时费力。笔者针对该类项目制图特点和需求，利用Python语言及其强大的第三方库，设计了自动制图程序，主要功能包括波浪线(不整合界线)、自动剪断线、微短线识别与删除以及多边形化成区等功能，目的是提高工作效率，减少制图过程中不必要的错误发生。……