基于MapGIS组件技术实现化探异常参数的计算
VBA和MapGIS在化探异常解释评价中的应用
中 图分 类 号 ： ６２ Ｐ ３
文献 标 识 码 ： Ａ
文章 编 号 ： ００— ９ ８ ２０ ）５— ３ ２— ５ １０ ８ １ （０ ８ ０ ０ ９ ０
化 探 是 地 球 化 学 找 矿 或 勘 查 地 球 化 学 的 简 称 ， 为地质 找矿 工作 的一种手 段 ， 作 在成 矿远 景
２ １ 工作 思 路 ．
均值 加 ２倍 离 差 计算 异 常下 限 ， 即单元 素 异 常
下 限 ＝均值 （ ）＋ ２倍 离 差 （ ６ 。 为单 元 素 的数 ２ ）
计及 编译 专业方 丽 的软 件 工具 包 等 ， Ｇ ｏ ｉ地 球 如 ｅｋｔ
其中， 均值加 ２ 倍离差法在有色金属矿产地质勘查
化探异 常解 释评价 工作 中应用 较多 。这种 方法数 据 处理量 大 ， 基本数 据 上 千 个 ， 算 过 程较 为繁 琐 ， 计 在 没有 专 门计 算 软件 的情 况下 ， 通过 编译 Ｖ Ａ程 序计 Ｂ 算处 理 ， 化 了计 算过 程 ， 高 了工 作效 率及成 果资 简 提 料 的精确 性 。其 程序 编译思 路 、 方法 如下 。
差 方 法 为 例 ， 讨 了运 用 Ｖ Ａ 编程 计 算 确 定 化 探 异 常 下 限 的基 本 思 路 和 一 般 方 法 ， 细 介 绍 了离 散 点 位 数 据 图 在 探 Ｂ 详
MAPGIS二次开发在化探元素异常提取中的应用
Ｇｅ ｏ ｃ ｈｅ ｍｉ ｃ ａ ｌ Ａｎｏ ｍａ ｌ ｙ Ｅｘｔ ｒ ａ ｃ ｔ ｉ ｏ ｎ Ｂａ ｓ ｅ ｄ ｏ ｎ ｔ ｈｅ Ｆｒ ａ ｃ ｔ ａ ｌ Ｃｏ ｎ ｔ ｅ ｎｔ Ａｒ ｅ ａ ａ ｎｄ Ｓ ｏ ｆ ｔ ｗａ ｒ ｅ Ｉ ｍｐ ｌ ｅ ｍｅ ｎｔ ａ ｔ ｉ ｏ ｎ
Ｚ ＨＡＮＧ Ｌｉ ｌ ， ２ Ｌ Ｉ Ｕ Ｂｉ ｎ Ｗ ＡＮ Ｇ Ｙｕａ ｎ ｍｉ ｎ
，
１ （ Ｋ ｅ ｙ Ｌ ａ ｂ ｆＧ ｏ ｅ ｏ ｍ ａ ｔ ｈ ｅ ｍ ａ ｔ ｉ ｃ ｓ ｏ ｆ Ｓ ｉ ｃ ｈ ｕ ａ ｎ Ｐ ｒ ｏ ｖ ｉ ｎ ｃ ｅ ， Ｃ ｈ ｅ ｎ ｇ ｄ ｕ ６ １ ０ ０ ５ ９ ， Ｃ ｈ ｉ ｎ ２（ Ｃ ｏ ｌ ｌ ｅ ｇ ｅ ｆＡ ｏ ｄ ｍ ｉ ｎ ｉ ｓ ｔ ｒ ａ ｔ ｉ ｖ ｅ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ｆＣ ｏ ｈ ｅ ｎ ｇ ｄ ｕ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ｆＴ ｏ ｅ ｃ ｈ ｎ ｏ ｌ ｏ ｇ ｙ ， Ｃ ｈ ｅ ｎ ｇ ｄ ｕ ６ １ ０ ０ ５ ９ ， ｃ ｈ ｉ ｎ ａ ） ３ （ Ｃ ｈ ｅ ｎ ｇ ｄ ｕＵ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｙｏ ｔ ｆＴ ｅ ｃ ｈ ｎ ｏ ｌ ｏ ｇ ｙ ， Ｃ ｈ ｅ ｎ ｇ ｄ ｕ ６ １ ０ ０ ５ ９ ， ｃ ｈ ｉ ｎ ａ ）
常下 限与传 统方法对 比， 得到的异常下 限值更加精确。同时， 通过 ＭＡ Ｐ ＧＩ Ｓ二次开发使地球化 学元素异常固定更加 自动化和智 能化。
本文 以湖北省铜绿 山铜 铁矿为例 ，取得 了较好 的效果。
关键 词：ＭＡ Ｐ ＧＩ Ｓ二次开发；分形含量面积；异常下 限；地球化学
中图分类号 ：Ｔ Ｐ ３ １ １ ． ５ 文献标识码 ：Ａ ＤＯＩ ：１ ０ ． ３ ９ ６ ９ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ３ — ６ ９ ７ ０ ． ２ ０ １ ３ ． ０ ５ ． ０ １ ６ 本文著录格式 ：［ １张立 １ ］ ， 刘斌 ， 王渊 民 ． ＭＡ Ｐ Ｇ Ｉ Ｓ二次开发在化探 元素异常提 取中的应用 Ⅱ １ ． 软件 ， ２ ０ １ ３ ， ３ ４（ ５ ）： ５ ６ — ５ ９
MAPGIS技术技巧-----化探数据处理步骤
化探数据处理步骤：
1、数据整理：
①将样品点位坐标和元素分析结果倒入Excel表
②将每个元素对应的坐标分别整理，另存为“。
Det“，CSV(逗号分隔)格式的文件，表中第一行表头为notgrid。
将文件保存到工作盘内。
3、打开MAPGIS的“空间分析”DTM分析”工作窗口，
4、打开三角剖分文件，即“Au.det”文件，
5点击Tin模型工具栏，选择“快速生成三角剖分网”
6删除三角剖分边，整理三角剖分网，删除不相关联线
7、追踪剖分等值线, 选择等值线套区等地设置后金星等值线层的设置后点击确定。
等值线层的设置是按照异常下限、2s、4s、8s三级异常进行设置。
形成下图。
8、分别另存点、线、面文件即可。
基于MAPGIS的分形方法确定化探异常
本文由中国地质调查局综合研究项目(编号:20010200015)资助。
改回日期:2003211218;责任编辑:宫月萱。
第一作者:李随民,男,1971年生,在读博士生,主要从事地学信息处理的研究;E -mail :smli71@ 。
基于MAPGIS 的分形方法确定化探异常李随民 姚书振石家庄经济学院,河北石家庄,050031摘 要 地球化学元素的异常下限值确定是地球化学中重要的问题之一,目前还没有一个令人满意的具有科学依据的计算方法。
传统的化探异常下限值计算是基于元素的地球化学分布呈正态分布或元素含量在空间上呈连续的变化这一假设为基础的,而事实上地球化学元素含量的空间分布是极其复杂的,研究表明,地球化学景观可能是一个具有低维吸引子的混沌系统,元素的地球化学背景值和异常具有各自独立的幂指数关系,由此导致了一种多重分形分布,因而可以利用元素的分形分布求出其异常下限。
利用元素的分形分布求异常下限的几种常用方法均是以求取不同尺度r 下,对应的N (r )数。
如果用手工统计方法计算,计算过程简单但繁琐,地理信息系统(GIS )可以实现图形单元的动态查询和属性统计,这种特性很适合于统计不同尺度r 下对应的N (r ),因此,可以将分形方法与地理信息系统加以结合,在GIS 平台上对地学图形和属性信息进行分析及统计处理,将原来比较繁杂的分形计算操作变得方便简捷,易于实现。
本文以武汉中地信息公司开发的MAPGIS 地理信息系统为例,对河北某地1∶5万Cu 元素化探数据进行了处理,具体说明了以MAPGIS 为工具,利用分形方法求取元素化探异常下限的过程,并将计算结果与传统异常下限计算方法所得的结果进行了对比,所得的异常下限值相近或一致,但由于传统的计算方法不仅要求对原始数据进行预处理,而且在下限值的确定上需要结合本地区的实际地质情况确定,因此人为干扰因素较大,而利用基于MAPGIS 的分形方法计算元素异常下限值不仅操作简单,而且计算时由于不受元素特高值的影响,因而不用对原始数据进行预处理,人为干扰因素较少。
基于MAPGIS平台的化探数据处理方法及异常图绘制
基于MAPGIS平台的化探数据处理方法及异常图绘制[摘要] 近些年来，随着经济的不断发展，我国地质勘测工作也有了更大的进步，尤其是地质绘图软件的广泛应用，为勘查工作带来了更有效的技术职称。
MAPGIS是一种广泛应用的制图技术，其具有较强的GIS制图功能，而且操作简单，在各种地质勘测中都较为实用，所以也受到了越来越多的关注。
本文就主要以MAPGIS6.7为研究对象，对其化探数据处理方法以及异常图绘制的相关问题进行简单的探讨。
[关键字] MAPGIS6.7 化探数据处理方法绘制异常图制图是地质工作中一项重要的内容，其不仅是对实际地质情况的展示，同时也是地质勘查人员工作的成果。
地质图也是开展各项地质工作的基本依据，对于对地质图的准确性和真实性也有着较高的要求，科学技术的不断进步为地质制图工作带来了更多的技术支撑，然而传统的地质图绘制技术在制度的精度和美观性等方面都存在着不同程度的缺陷，所以必须要对其进行不断的改进与完善，MAPGIS6.7 技术便是在这种情况被广泛的采纳，MAPGIS技术的运用，使得地质绘图更加准确和美观。
由于实际的地质情况十分复杂，所以在利用MAPGIS技术进行绘制时，难免会存在圈定不合理的现象，所以需要按照规定形成相应的异常图，以此来保证地质绘图的真实性。
1 MAPGIS6.7的化探数据处理方法在地址绘图的过程中，主要是通过选择坐标系和数字化绘图两个主要的环节来完成。
对于坐标系的选择，一般可以利用已知的坐标进行对比和计算来获得。
数字化绘图中，则包含了数据的导入、数据的采集和处理、数据的误差校对等几个关键的环节，只有确保各个环节的有效性和连续性，才能确保地质图的清晰和准确。
其中，数据处理是制图的核心环节。
1.1 原始数据的处理在进行实际的勘测过程中，对于地质沉积物的测量一般需要通过GPS技术对采样点进行定位人，获得相应的采样坐标X和Y，并且将这些原始的数据导入到数据处理表格中进行分析，进而获得不同元素的分析结果，利用制表符对不同的文本文件之间进行间隔，形成名称为“化探结果”的原始数据记录。
化探数据处理步骤
所有数据均输入计算机、以MAPGIS 制图系统为平台，以原始数据筛选替换特高值后，转计算成对数值后，用 X +2S 求出异常下限，分别以X －2S 、X －0.5S 、X +0.5S 、X +2S 分出色区，绘制各元素地球化学图，以上做图过程均在计算机上用MAPGIS 软件完成。
对化验室的样品分析结果取对数分组作直方图，证明所有元素均符合对数正态分布。
元素异常参数的确定：首先对原始数据进行假设正态检验，再作X －
+3S 特高值逐步剔除，然后进行各参数统计。
Au 元素含量为W×10－9，其它元素含量为W×10－6。
（1）背景平均值：f
fxc X ∑∑= （2）对数标准离差：1)(22-∑-∑=
n n fxc fxL S （3）对数异常下限：T0=X －L+2S
（4）变异系数：%1001%2230285.2⨯-=⋅S e Cv
e －自然对数，2.30285为常用对数与自然对数模数的倒数
S －对数标准离差，1为常数
（5）衬度：To
Xa K =（Xa 为异常平均值） （6）异常规模：P=k×km 2（km 2为异常面积）
各类系数计算和所利用的公式均符合规范要求。
MapGIS在制作地球化学异常图中的应用
MapGIS在制作地球化学异常图中的应用
王璨;张艳霞
【期刊名称】《黑龙江国土资源》
【年(卷),期】2009(000)005
【摘要】@@ 地球化学勘查始于20世纪50年代,随着找矿难度的加大,地球化学勘查已经不是单一的学科,目前以勘查地球化学为中心.综合地球化学物理空间、地形空间的数据处理方案,形成了一套计算机图像处理系统(MapGIs).
【总页数】1页(P61)
【作者】王璨;张艳霞
【作者单位】黑龙江省地质矿产测试研究所;黑龙江省地质矿产测试研究所
【正文语种】中文
【中图分类】P2
【相关文献】
1.MAPGIS在地球化学制图中的应用 [J], 潘洪儒;秦振宇
2.应用MAPGIS数字高程模型提取区域地球化学异常信息的方法探讨 [J], 陈希清;杨晓君;陈富文;高立宝
3.应用MAPGIS制作地球化学图单元素异常图及综合异常图 [J], 袁义生;刘应忠;罗明学;何彦南
4.Mapgis在地球化学找矿中异常特征参数提取的应用 [J], 丁浩;申开洪;曲锦
5.MAPGIS在地球化学元素异常图制作中的应用 [J], 窦继忠
因版权原因，仅展示原文概要，查看原文内容请购买。
基于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。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
图层 管理组 件文件 ： ｐ ｓ ａｅＣ ｍ．ｌ ＭａＧｉ ｙｒ ｏ ｄｌ Ｌ
地 图、 图例管理组 件文件 ： ｐ ｉ ｐ ｏ ｄｌ Ｍａ Ｇ ｓ Ｃ ｍ． ｌ ｍａ
图库 管理组 件文件 ： ｐ ｓ Ｂ Ｃ ｍ． ｌ ＭａＧｉ Ｓ ｏ ｄｌ Ｄ 图像分析管理组件 文件 ： ｐ ｉ ａｔｒ ｏ ． ｌ Ｍａ Ｇ ｓ ｓｅ ｍ ｄ Ｒ Ｃ ｌ 图形 编辑控 件文 件 ：ｄｔ ｅ ｏｘ Ｅ ｉ ｗ． ｅ Ｖｉ 工程管 理控 件文件 ： ｐ ｅ．Ｃ ＭａＴｒｅＯＸ
摘 要 ：在 Ｖ ｓａ Ｂ ｓ 环境中编写一个求化探异常参数的程序 ， ｉ ｌ ａｉ ｕ ｃ 程序首先运 用 ＭａＧ Ｓ ｐ Ｉ 组件 （ Ｏ 技术 ， Ｃ Ｍ） 构造 了含有元素分析数据的点工作 区， 然后读入 Ｍａ Ｇ Ｓ异常 区文件 （ｗ ） ｐＩ ． ｐ 。通 过点对 区的相交分析 ， 提取 出异常数据 ， 进而对数据进行计算， 求得异常参数。程序结果输 出到 Ｅ Ｃ Ｌ工作表 中， ＸＥ 非常实用。 关键词 ：Ｖ ｓａ Ｂ ｓ ； ｐ Ｉ ； ｉ ｌ ａｉ ＭａＧ Ｓ 组件 ； ｕ ｃ 化探异常参数 ； 计算
式 而研发 的一 种新 的软 件 开 发技 术 。组 件 有 自己 的属 性和 方法 ， 属性 是 组 件对 象数 据 的 特性 ， 方法 则是 组件 对象 的行 为 。在 Ｃ ＯＭ 构 架下 ， 们 可 以 人
琐 的工作 ， 内、 目前 都 是从 数 据 图上利 用 手 工 国 外 逐 点输 入元 素 的异 常数 据 ， 异 常 图上 得 到 面 积 ， 从 然 后求 取 各 项 参 数 。手 工 方 法 费 时费 力 ， 易 出 且 错 。作 者 在 本 文 中 利 用 ＭａＧＩ ｐ Ｓ组 件 技 术 ， 在
属性 编辑控 件文件 ： ｔ ｔｕ ｄｔＯ Ｘ ＡｔＳｒ Ｅ ｉ Ｃ ．
１ 期
刘俊长等 ： 于 ＭａＧ Ｓ 基 ｐ Ｉ 组件技术实现化探异常参数的计算
１３ １
属 性 结构 编辑 控件 文件 ： ｓ ｔ ｄｔＯＸ Ｇｉ Ｅ ｉ Ｃ Ａｔ ．
１ ２ 在 Ｖｉ ａ Ｂ ｓ ． ｓ ｌ ａｉ ｕ ｃ中使 用 Ｍａ ＧＩ 件 ｐ Ｓ组
Ｖｉ ａ ａｉ 环境 下开 发 出一 个化探 异 常参数 的计 ｓ ｌ ｓ ｕ Ｂ ｃ
开发 出各 种各样 的功 能专 一的组 件 ， 后将 它们按 然
照需 要组 合起来 ， 构成 复杂 的应用 系统 。 ＭＡ ＧＩ Ｐ Ｓ以组件 的方式 提 供 了完善 的 ＧＩ Ｓ功
算程序 ， 该程序能够 自动获得异 常数据 和异常面 积 。程 序只 需输入 Ｅ Ｅ ＸＣ Ｌ原 始 数据 文件 和 Ｍａ － ｐ ＧＩ 常 区文 件 ， 可 输 出各 元 素 的化 探 异 常 参 Ｓ异 即
数 ， 出结果 保存在 Ｅ Ｅ 输 ＸＣ Ｌ中。
Байду номын сангаас
能， 包含 在若干 动态链 接库 （ Ｌ ） Ａｃｉｅ Ｄ Ｌ和 ｔ Ｘ控件 ｖ （ Ｘ） 件 中＿ ］具体 如下 ： ＯＣ 文 ｌ ， ４
基础 组件 文件 ： ｐ ｓ ａＣ ｍ１ ｄｌ Ｍａ Ｇｉ ｓ ｏ ． ｌ Ｂ 基 础 控 件 文 件 ： ｐ ｉ ａＣ １ ｄｌ Ｍａ Ｇｉ Ｍａ Ｇ ｓ ｓ ｔ ｌ Ｂ ． 、 ｐ ｓ —
Ｓ ｔＭ ｙ ｈ ｅ ｅ ｓ ｅ ｔ— Ｍ ｙ ｏ ｋ Ｗ ｏ ｋ ｈ ｅｓ １ Ｂ ｏ ． ｒ ｓ ｅｔ（ ）
（） 行 ＭａＧＩ １运 ｐ Ｓ二次 开发 包 中的 ｓｔｐ ｅｅ ｅｕ ．ｘ ， 或 是对 需 要 的组 件 进 行 注 册 ： ｉｄ ｗｓ 始 菜 单 Ｗ ｎｏ 开 ＞ 运 ｒｇ ｖ３ ｔ ｐ ｉ ７｝ ｒｇ ａ ｔ ｅ ｓｒ ２ Ｃ ｍａ ｇｓ ６ ｐ ｏ ｒｍ
区域化 探 中的大量 异常 ， 此项 工作是 一项 细致 而繁
技 术 的二 次 开 发 可 以很 好 地 引 入 Ｖｉ ａ ａｉ ｓ ｌ ｓ ｕ Ｂ ｃ开
发 环境 中 。
１ １ Ｍａ Ｇ Ｓ组件介 绍 ． ｐＩ
组件 技 术 ， 是 组 件 对 象 模 型 （ ｏ ｏ ｅ ｔ 即 Ｃ ｍｐ ｎｎ ＯｂｅｔＭｏ ｅ， 写 Ｃ Ｍ ）是 对 对 象 及 其 数 据 和 ｊｃ ｄｌ缩 Ｏ ， 方 法 的封 装 ， 使软件 生产 更加符 合人 类 的行为方 是
ＢａＣｔ１ ｄｌＭ ａ ＧｉＢ ｓ ｌ． ｌ ｓ ｌ． ｌ、 ｐ ｓ ａ Ｃｔ２ ｄｌ
图形 显示组 件文 件 ： ｐ ｓ ｓ Ｃ ｍ． ｌ Ｍａ Ｇｉ ｐ ｏ ｄｌ Ｄ
１ ＭａＧＩ ｐ Ｓ组件技术
ＭａＧ Ｓ是地 质行 业 使 用最 广 泛 的 Ｇ Ｓ平 ｐＩ Ｉ 台［ 作为一款优秀的 Ｇ Ｓ软件 ， 翻， Ｉ 提供 了完备 的 数字制图、 数据库管理及空间分析等功能 。为 了 使 用 户更 方便 地使 用这 一 空 间信息 管理 系 统，
中图分 类号 ：Ｔ １ ． Ｐ３ ７４ 文献 标识 码 ：Ａ
ＭａＧ Ｓ还 提供 了二 次 开 发 功 能 ， 中 ， 于组 件 ｐＩ 其 基
０ 前 言
新 一轮 的区域 化探 数据 更新及 各省 矿调工 作 ， 获得 了大量 的地 球化 学异常 ， 中一项 重要 的工作 其 就是化 探异 常登 记 。异 常 的登 记 工 作 包括 计 算 一 系列 的异常 参数 ， 主要 有 各元 素 的异 常 面积 、 高 最 值 、 均值 、 平 衬度 、 模 、 度 分 带数 目等 。相 对 于 规 浓
第３ 卷 第１ ４ 期
物 探 化 探 计 算 技 术
２１ 年１ ０２ 月
文 章编号 ：１ ０一 ｌ ４ （０ ２ Ｏ — ０ １一 Ｏ Ｏ１ ７９２ １） １ １２ ４
基 于 Ｍａ ＧＩ ｐ Ｓ组件 技 术 实现 化 探 异 常参 数 的计算
刘俊长， 龚红蕾，陈军威， 刘军恒
（ 河北省 地球 物理勘 查 院 ， 北 廊 坊 河 Ｏ ５Ｏ ） ６ Ｏ Ｏ