基于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采集的坐标数据以及样品元素分析结果，具有成图快、精度高、准确、可操作性强的特点。其绘制