基于MAPGIS6.7化探数据处理方法及异常图绘制

基于MAPGIS6.7的化探数据处理方法及异常图的绘制
[摘要] 本文以苏丹红海州捷拜拉吐哈度57区金矿区水系沉积物测量中au元素化探数据为例，介绍了mapgis6.7环境中化探数据的处理方法和流程，以及化探异常图的绘制。

[关键词] mapgis6.7 化探数据处理方法化探异常图
1. 前言
我队在苏丹红海州捷拜拉吐哈度57区金矿区开展了1∶5万～1∶20万水系沉积物测量，所有采样点位均使用gps精确定位，共采集样品数量2200个，并由化验单位完成对化探样品au元素含量的分析。

mapgis6.7是一个大型基础地理信息系统软件平台，它可以方便、快捷、有效的对经过化验分析的数据进行处理并自动圈定元素异常
[1]。

2. 数据处理
2.1 原始数据整理
水系沉积物测量通过gps对采样点进行定位，得到采样点坐标（x、y）。

在excel中导入原始数据字段，包括采样点编号、x坐标、y
坐标、au元素分析结果，然后另存为制表符分隔的文本文件，命名为“化探结果.txt”。

2.2 背景值和异常下限值的确定
异常下限主要根据剔除高值点后样品au元素含量的算术平均值和标准离差来确定。

其计算公式为：
其中，——au元素含量的异常下限值；
——剔除高值点剩余样品au元素含量的算术平均值；
——剔除高值点剩余样品au元素含量的标准离差值；
为经验系数，一般取1.65～3。

本文的取值为： =3。

统计计算时，利用excel相关函数将au元素含量所在的列进行从大到小排列，然后计算全部样品中au元素含量的算术平均值和标准离差值，将大于的含量数据剔除后，再次统计计算全部样品中au元素含量的算术平均值和标准离差，直至将大于的含量数据全部剔除后，用剩余样品中au元素含量的算术平均值来度量au元素的背景值，并以剩余样品的算术平均值加3倍的标准离差值来度量au元素的异常下限值。

表1为剔除高值点后au元素的相关统计数据。

表1 剔除高值点后au元素的相关统计数据
最大值最小值平均值标准离差异常下限
3.55 0.10 1.36 0.73 3.56
本文在绘制化探异常图时，将背景值设定为1.36，异常下限值设定为4.0。

3. 利用mapgis绘制化探异常图
3.1 数据的投影变换处理
投影变换处理就是将化探采样点坐标通过mapgis处理生成化探数据点位图。

投影变换在mapgis6.7中的操作步聚：
（1）在mapgis6.7下单击实用服务→投影变换→用户文件投影
变换，如图1；
（2）打开文本文件“化探结果.txt”，点击按指定分隔符→设置分隔符，将分隔符号中的“tab键”打勾，“空格”中的勾取消,结果显示如图2；
（3）设置x坐标、y坐标所对应的的列号；
（4）将用户投影参数、结果投影参数分别根据需求进行设置；（5）设置投影点图元参数，然后投影变换→确定，投影结果如图3；
（6）保存点文件，命名为“化探结果.wt”。

3.2 数据属性连接
数据属性连接是将投影点位图与化探数据相关联，实现了将外部数据库的属性数据输入到mapgis图形文件的属性数据中[1]。

属性连接在mapgis6.7中的操作步聚：
（1）在access中打开“化探结果.txt”→按步骤（默认选项）生成“化探结果.mdb”→在access表文件中右击“化探结果”→导出（保存类型：dbaseⅲ（*.dbf）），生成“化探结果.dbf”；（2）在mapgis6.7下单击库管理→属性库管理；
（3）在mapgis属性管理子系统中单击属性→连接属性；
（4）打开连接文件“化探结果.wt”，打开被连文件“化探结果.dbf”（注：均按序号连接），结果如图4；
3.3 离散网格数据化
化探数据都是离散数据，为了绘制平面等值线图，需要对原始数
据进行网格化处理，即生成grd模型。

mapgis6.7环境下生成grd 模型的步骤：
（1）在mapgis6.7下单击空间分析→dtm管理，即打开mapgis 数字地面模型子系统；
（2）单击“文件”→打开数据文件→点数据文件，打开“化探结果.txt”；
（3）单击“处理点线”→点数据高程点提取→选择au元素含量值为高程值（即字段4）；
（4）单击“grd模型”→离散网格数据化，如图5；
（5）根据需求设置相应参数，生成“化探数据.grd”。

3.4 绘制平面等值线图
化探采样就是为了形成该地区各元素的化学异常图，圈定工作区的矿产异常范围，为后期地质找矿提供依据[2]。

绘制平面等值线图就可以达到这个效果。

mapgis6.7环境下生成平面等值线图的步骤：
（1）单击“grd模型”→平面等值线图绘制，打开“化探数据.grd”；（2）设置相关参数：
a、选中等值线套区、绘制色阶、保留边界线、等值线光滑处理（中程度）、制图注记，在制图幅面中选择原始数据范围；
b、等值层值和区参数的设定方法如表2；
表2 等值层值和区参数表
色区着色及（区名）蓝
（低值区）浅蓝
（低背景区）黄色
（背景区）褐色
（高背景区）红色
（高值区）
元素含量值范围～
～
～
～
>
c、注记参数：
需要在等值线上标注数值的选择yes，不需要标注的选择no；（3）等值线参数设置结果如图6，再单击确定，等值线图会自动绘制完毕；
（4）单击“文件”→存数据于，分别存点数据文件为“化探异常图.wt”、线数据文件为“化探异常图.wl”、区数据文件为“化探异常图.wp”。

（5）在mapgis编辑子系统中打开上面的点、线、区文件，即可对等值线图进行编辑。

结果如图7。

4. 结束语
应用mapgis6.7来处理分析gps采集的坐标数据以及样品元素分析结果，具有成图快、精度高、准确、可操作性强的特点。

其绘制
的异常图能与地形地质图很好的套合，有利于异常的解释并进一步圈定异常区，为后期勘查工作提供依据。

但因参数设置的不同，生成的等值线图还会存在有线条生硬、圈定不合理等问题，需要不断的尝试，以使之更加合理、美观。

参考文献：
[1] 武汉中地信息工程有限公司.mapgis地理信息系统实用教程[m].湖北武汉,2002.
[2] 何明华.mapgis数字高程模拟系统在化探数据处理中的应用[j].矿山测量,2005,(2):13～15.
作者简介：
宋雪枫（1985—），男，河南省郑州市人，助理工程师，主要从事工程地质与矿产地质方面的研究。

李威（1989—），男，河南省驻马店市人，主要从事工程地质与矿产地质方面的研究。