0 引言
航磁异常是航空磁测的重要成果之一,它是以航磁原始测量数据为基础,辅以位场转化数据,结合区域地质构造及其他物化探资料,综合分析后进行拾取的局部异常。通过对航磁异常的定性分析及定量解释,可直接或间接地寻找矿床、区分和圈定各类磁性地质体、划分某些成矿远景区、研究地质构造,以及为国民经济各部门提供基础性地球物理资料[1]。
我国自开始进行航空磁法勘测至今已有60余年,测量比例尺由1∶100万至1∶1万不等,已覆盖了中国大部分陆地和海域[2]。以航磁异常为主的航空物探异常已成为地质找矿的重要信息之一,它不但对实时的找矿具有指导作用,更加对研究区域性地质构造具有重要意义,对地质找矿战略性部署提供重要参考。然而长期以来,由于严重缺乏对航磁异常信息全面、科学、有效的管理,虽然建立了基于Oracle的航空物探数据库,对航空物探数据进行存储管理[3],但对于航磁异常数据而言仅仅将数据保存至数据库中并不能充分发挥其作用,航磁异常的图面展示受到很多技术上的制约,一定程度上限制了其服务提供。基于此,笔者在对航磁异常数据进行深入研究的基础上,结合大量图件及相关规范,采用二次开发的形式,实现了对存储在Oracle数据库中的航磁异常数据进行自动渲染的功能。
1 渲染方法研究
1.1 管理平台研究
图1给出了航磁数据的生产流程,航空物探工作以项目为单位通过野外飞行作业获取到的航磁ΔT数据,经数据处理解释专家对数据进行处理分析,形成本地数据,主要包括本地数据库形式的基础数据、图层数据形式的解释评价成果数据、网格成果数据;本地数据经人工数字制图工作后,形成最终成果图件。
图1
根据信息化建设对全部数据的管理需要,通过数据采集软件对成果数据及图件进行采集录入,存储至以ArcSDE为引擎进行管理的Oracle数据库中。为形成可发布的图件,需要再次对数据进行渲染,因此需大量的人工干预进行数字制图。
笔者通过对数据库管理引擎的研究,明确了基于Oracle数据库的制图渲染的基本机制,即通过ArcSDE连接数据库,使用ArcGIS平台进行操作,依据数据属性列手动进行子类划分,确定渲染方式(唯一值渲染或多字段唯一值渲染),使用*.style文件形式的符号库与属性内容进行关联匹配,依照相关行业标准规范为数据要素选择对应的制图符号,若无符合要求的符号,则需打开符号管理器,新建或者在已有符号基础上制作完成所需符号,匹配完成后,刷新图面并将符号信息进行保存[3,4,5,6]。在渲染过程中,主要难点有二:其一,项目工作中,是按照解释工作逐步推进,逐一对异常数据进行分类及渲染,此种方法明显不适用于大数据量的数据库管理;其二,数据库平台的渲染方式是对数据进行分类,而后对每一类别进行渲染,最后将渲染信息保存至工程文件中,即数据库每更新一次,就要对数据重新进行一次渲染工作,工作量巨大且重复。基于此,亟需一种能够满足数据表达需求的自动渲染工具。
1.2 渲染方案制定及符号制作
航磁异常数据作为点类型数据存储于数据库中,在制图表达过程中所需的信息有五种,其一为坐标信息,以双精度浮点型存储,记录异常点的经纬度信息;其二,异常轴向为双精度浮点型,记录异常的轴向角度,以正东方向为零,遵循右手螺旋法则,单位为度;其三,异常性质为无符号整型,记录解释专家对磁异常的定性划分,通过数据字典可查找原始语义,即引起磁异常的岩体性质,如04表示由辉长岩引起的局部异常;其四,异常分类为无符号整型,记录解释专家对异常的半定量级别划分,如01表示甲一类异常;其五,异常编号为可变长度的字符串类型,用于记录解释专家为方便描述对异常进性的编号,如新C-2014-126,表示该异常地处新疆,为航磁异常,测量于2014年,在测量区内排序编号126。
由于异常点图层的图面表达信息非常丰富,与多个属性相关,每一属性又有各自的数据字典,各数据字典项对应一符号,若想表达完全,则需从三个符号类别中分别选出一个再进行叠加,而若将三类符号全部排列组合在一起形成一个新的符号,则会使符号库指数级增长,特别是异常轴向信息属性,其并非数据字典,是数值型数据,因其无固定值,在图面表达时是根据其轴向对符号角度值进行设置,也就无法成为固定的符号,因此在项目制图工作对航磁异常点进行表达时,一般使用叠加图层的方式来进行渲染,即使用三个相同的点图层进行不同的符号化,叠加显示在图面上以表达其相关信息。如表1所示,其中一数据层用于表达异常分类,选择中空的符号来表达,且不遮挡异常性质的花纹及颜色,一数据层用于表达异常性质,用颜色及花纹表达。同理,异常轴向也用一数据选用中空符号,外加轴向示意表达,对于异常编号,则对任一图层进行注释标注,完成整个异常的图面表达。
表1 航磁异常数据渲染方式
Table 1
| 属性类 | 属性信息 | 符号形式 |
|---|---|---|
| 异常类别 | 乙3 |  |
| 异常性质 | 辉长岩引起的局部异常 |  |
| 异常轴向 | 0 |  |
| 异常编号 | 新C-2014-192 | 新C-2014-192 |
| 异常位置 | (78.000,39.000) |  |
以上制图工作中所需的符号形式由《航空物探数字制图与印刷技术规定》明确,然而在实际工作中,规定所列举的符号信息有所欠缺,重要信息异常类别,体现了引起异常的岩体性质,其多样性基本涵盖了大部分岩石种类,在项目工作中,由于是手工制图,常常选取就近花纹的已有符号,在图例中加以标注区分。此种方法对小片区具有特定找矿方向的图件较为实用,但对于大范围、地质构造多样、岩性种类各异、异常数量多的综合研究图件而言,此种方法会导致在图面上使用相同符号表现不同性质异常的现象。因此,笔者结合数据库数据及相关图件,参照国家标准《GBT958-2015区域地质图图例》及其他相关资料,对航空物探工作中所用到的航磁异常类别符号进行了整理归纳,完善异常类别符号。在符号制作时,使用矢量字体文件,通过符号管理器,在模板的基础上制作完成相应符号,形成航空物探异常符号库[8]。
航空物探数据库中的异常数据,由于数据属性齐全,要表达在图面上的信息,已划分为有数据字典的属性字段,而通过数据采集软件录入至数据库中的数据,属性完备,可直接根据属性字段内容与符号库内符号进行匹配并进行其他属性设置,完成数据渲染。
2 自动渲染实现
为简化人工操作,实现数据库数据的自动渲染,主要需解决以下问题:第一,航空物探数据库的建设是为了实现信息化管理,但先前工作中多图层叠加显示的方法要在数据库中对异常数据图层存储三遍,有违信息化管理的初衷[6,7,8,9,10],本软件要对仅在数据库中存储一次的数据进行自动制图;第二,航空物探数据库并非针对某一项目,而是为将所有项目数据集中管理,因此,数据库内数据是在不断增长,以往针对项目的图件制作不存在此问题,若要对新增数据进行渲染,只能将数据拷贝至制图界面重新手工制作,工作量大,易出错,而航空物探数据库内信息量巨大,难以实现[11,12,13,14,15]。因此,笔者在航空物探数据库管理引擎的基础上,设计了对航磁异常单图层数据进行制图渲染的流程,并使用原引擎接口进行二次开发,实现航磁异常自动渲染。
2.1 渲染流程设计
图2
2.2 接口设计实现
航空物探数据库通过ArcSDE进行管理,在授权终端机上可通过ArcGIS平台进行属性查看及相关制图、空间分析等操作,因此笔者通过二次开发,集成原ArcEngine相关接口及虚拟类,重新定义实现自动渲染功能所需的类对象,并进行插件式开发,将其集成为ArcGIS平台的拓展工具,在数据进入数据库中后,对其进行自动渲染,并将符号信息记录至属性库中。实现自动渲染所需类主要包括三部分:实现数据库读取及写入等操作的DataBaseop类;生成符号及标注的MagNorSymbol类及Label类,其中MagNorSymbol类为PointSymbol类的派生类,PointSynbol类为点类符号的基类,实现点符号的属性控制,MagNorSymbol类则是针对航磁异常所设计的符号类,通过数据库接口读取到的属性合成最终所需的航磁异常符号;最后是实现图面渲染的MagRender类,定义渲染器并完成渲染。
图3
3 应用实例
笔者依据上述设计,构建了航磁异常自动渲染工具,并依照属性对航磁异常进行制图渲染。如图4所示,笔者选取部分数据进行还原,通过对比,可见二者基本一致,均符合相关标准,达到根据数据对图件进行恢复的目标。自动渲染功能的研发大大拓展了航磁异常数据的应用范围,不但能够将所有项目的航磁异常分布实时更新至一张图,还能够与其他解释评价数据相互叠加进行分析应用,拓展了航磁数据的应用领域。
图4
4 结束语
数据渲染成图是航空物探工作的重要成果之一,目前数据难以实现跨项目、跨平台渲染,因此笔者以C#为编程语言,构建航磁异常数据自动渲染工具,实现了工具的研发,最终成图效果满足航空物探工作图件的要求,符合相关国家及行业规范。作为公共函数,提供了标准的API接口,满足不同开发语言的调用,也方便对其进行统一维护。此种方法思路及对象类,不仅可用于航磁异常的自动渲染,更可加以应用,对其他点类数据进行数据渲染,推动航空物探数据的信息化应用。
参考文献
中国航磁异常管理信息系统开发与建设
[J].
The Development and construction of aeromagnetic anomaly management information system of China
[J].
航空物探测量数据共享技术研究与应用
[J].
The study and application of sharing technology to aeromagnetic data
[J].
海洋航空物探数据库建设进展
[J].
A marine aeromagnetic database for efficient data management
[J].
制图新技术与传统地图编制的比较和影响
[J].
Comparison and influence between new cartographic technology and traditional map compilation
[J].
地球物理数据处理解释系统中坐标转换的设计及实现
[J].
The design and software realization of coordinate transormation for the geophysical data processing and interpretation system
[J].
基于ArcGIS的电子地图快速制图方法研究
[J].
Research on quick mapping method of electronic map based on ArcGIS
[J].
全矢量地质符号库的设计与实现
[D].
Design and implementation of full vector geological symbol library
[D].
基于ArcSDE的空间数据库的创建及应用
[J].
Creation and application of spatial database based on ArcSDE
[J].
应急场景快速制图及地图信息传输
[J].
Rapid mapping of emergency scenario and cartographic information transmission
[J].
基于数据库驱动的地质图快速制图技术研究
[J].
Study on the technology of rapid geological mapping based on the database drive
[J].
利用GIS数据进行地图制图表达与更新
[J].
Use GIS data to express and update maps
[J].
The establishment of the SAR images database system based on Oracle and ArcSDE
[J].
一体化地图辅助设计系统的设计与实现
[J].
Design and Implementation of aided design system of integrated maps
[J].
浅谈基于Arcgis自动化专题制图的关键技术
[J].
Introduction to the key technology of automation thematic mapping based on Arcgis
[J].
航空物探工作程度图批量定制技术研究
[J].
Research on technology of batch customizing aero geophysical survey regions maps
[J].
基于ArcSDE的空间数据库技术的运用分析
[J].
Application analysis of spatial database technology based on ArcSDE
[J].
基于地理信息框架数据库出图的关键技术
[J].
Mapping key technology of based on geographical information frame database
[J].
浅谈新疆地质灾害快速制图系统
[J].
A brief discussion on the rapid mapping system of geological Hazards in Xinjiang
[J].
