下载文档原格式
用Excel及Mapgis处理化探数据1、软件要求：Microsoft Office Excel 2003及MAPGIS672、参考资料：如何利用Excel处理化探数据，应用MAPGIS制作地球化学图单元素异常图及综合异常图3、Excel准备工作：①在化探找矿工作中 ,大量数据的处理是一项重要的前期基础工作 ,数据处理的及时、准确与否 ,将直接涉及到下一步勘查工作的部署。
作者在本文中以 Microsoft Office XP Excel工具～数据分析功能为工具 ,简述了这些功能在化探数据处理当中的一些应用。
缺省的 Windows Professional XP一般不安装用来数据处理所需要的数理统计功能 ,故需重新加载 ,步骤如下:工具 -加载宏 -分析工具库 -确定。
之后 ,还需检查在工具菜单下有无数据分析选项。
在大量数据的录入过程中 ,可以设定工具 -语音选项。
在某个数据输入完毕之后 ,按 Enter键即可语音朗读 (需要音响或耳麦 ) ,这样可以实现数据录入的同步检查 ,确保录入数据的准确性。
②数据的处理正态性检验数据是否服从正态分布(或对数正态分布)是化探数据处理的一个重要前提条件，因此必须先进行数据正态性检验。
4、Excel操作步骤：1）、Excel2003准备：把测试数据及样号坐标整理为下图所示2）、打开：①工具——数据分析——描述统计（如下图）②再点确定如下图（Au的数据统计）③再点确定如下图（Au的数据统计）这样就可以得出一系列所需数据，但经验算该组数据为非正态分布，不能做进一步计算，但此为计算流程。
5、MAPGIS操作步骤：1）、首先将Excel数据复制粘贴为记事本格式（.TXT格式）2）打开Mapgis如下①、选择②、Grd模型——离散数据网格化③、打开（做Au图）选择数据列D（X、Y、Z）—网格化方法M（按不同数据选择不同网格化）—输出网格文件名F（存储数据的位置）—确定。
④、选择打开上面保存的⑤、打开⑥、确定（如下图）⑦、选择文件—数据存储于—点数据文件、线数据文件、区数据文件（以上点线区文件保存完毕，以下就是用Mapgis正常编辑）⑧修改要编辑的项目化探直方图制图软件的编程方法2008-12-21 23:36软件采用VB+MAPGIS SDK 开发数据载入格式： CSV文本文件，内容为包含分区编码的元素分析数据分区编码与元素代号均采用MAPGIS的书写方式如C1表示为C#-1，以便成图时正确显示。
最新最详细化探数据处理与编图流程⼀、指导思想成矿地质背景地球化学研究就是从地球化学特征出发，借助已建⽴的地球化学信息提取技术，充分利⽤地球化学调查所获得的海量数据信息，提取有关反应成矿地质背景条件的地球化学信息，并编制相应地球化学图及相应的推断解释图件，为资源潜⼒评价有关成矿地质背景的研究提供地球化学⽀撑。
⼆、⼯作内容（⼀）基础图件成矿地质背景条件的地球化学信息提取⾸先是要编制有关基础地球化学图件。
主要有：1. 单元素（化合物）地球化学图2. 地球化学组合异常图3. 地球化学综合异常图（⼆）解释推断图件地球化学解释推断图件，内容包括：1. 地球化学推断解译地质图2. 地球化学找矿预测图三、⼯作⽅法（⼀）数据校正处理1|数据检查的必要性，因为实验室的分析报告还是⼿⼯输⼊的，还是存在录⼊错误的，我们重点检查的是“>”，数据中间的空格等录⼊错误问题；另外还有畸变检查，数据的特⼤值,⽐如超过10倍变差，⼀般对这样的分析值实验室会很重视的,你也可以提出让他们再确认⼀下,做到⼼中有数。
另⼀类错误可能会是我们录⼊样号或者坐标时出现的错误，如：“56b” 写成“56 b”,程序是以空格分开数据的，数据如果写成这样就会产⽣错误结果，有时在完成处理后才可能发现，这样⼀来我们前⾯的⼯作就作废了。
所以数据检查是⾮常必要的。
2|异常下限值的确定采⽤逐渐剔除法：①计算全区各元素原始数据的均值(X)和标准偏差(S)；②按X1+3S1的条件剔除⼀批⾼值后获得⼀个新数据集,再计算此数据集的均值(X2)和标准偏差(S2)；③重复第⼆步，直⾄⽆特⾼值点存在，求出最终数据集的均值(X)和标准偏差(S),则X做为背景平均值，S为标准离差，T(异常下限值)= X (背景平均值)+2S(标准离差)求出理论异常下限值，再结合地球化学等量线、地质背景及圈定效果确定出实⽤异常下限值。
3|重复样样品合格率统计野外重采样品以密码样形式插⼊样品中进⾏了分析，结果(C2)与第⼀次分析结果(C1)进⾏了⽐对。
1:1万化探的数据处理成图方法探讨毕武1、2段新力1、2袁小龙1、2黄显义1、2彭仲秋1、21.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，8300912.新疆地矿局物化探大队信息中心，新疆，昌吉，831100摘要：对1:1万化探的数据处理成图方法，在实际工作中的一些应用经验。
关键词：1:1万化探数据处理成图方法0前言随着1:5万区域化探的实施，在1:5万区域化探异常范围内开展的1:1万化探勘查也越来越普遍，对1:1万化探目前没有具体的规范，一般是参照1:5万区域化探规范执行，可是成图效果却有时不尽理想，下面就两个例子举例说明1:1万化探的数据处理成图方法。
1编制的步骤及要求1.1图件编制的要求：a土壤地球化学测量图件分两部分：基础图与推断解释图。
b区域调查和普查工作的图件编制按化探区域调查有关规范执行。
c图件编制必须符合地球化学勘查图式图例及用色标准规定。
d成果报告需做交通测区位置图、实际材料图、等值线图、综合异常图及其它推断解释图。
1.2例子一：南北向测线，100×40网格距。
这样的数据本身符合网格数据的格式，只是以TXT文本格式保存，所以我们可以借助GeoIPAS软件的“数据预处理”→“数据格式转换”→“XYZ数据转为网格数据”功能，这里注意数据坐标的起始值、终止值和网格间距，网格间距就选择100×40，起始值和终止值要保证数据的原始点位都与网格点重合。
这里有一个问题，就是有些项目用的是GPS坐标，不一定是与网格点位完全重合，一般我们建议用规则网的坐标。
下面先看一下图1-1点位数据图和图1-2点位符号图；图1-3原始数据转为网格数据绘制的地球化学图以及其转换参数：X最小值320000，X最大值32500网格化0，间距100，列数51；Y最小值4121000，Y最大值4127000，间距40，行数151。
图1-4距离平方导数加权网格化绘制的地球化学图，GeoIPAS系统→数据预处理→离散数据网格化→距离平方导数；搜索范围：圆域，R=300米；网格化参数：X最小值319880，X最大值325120，间距40，列数132；Y最小值41209200，Y最大值4127080，间距40，行数155。
Mapgis勾绘地质图的基本操作一、扫描底图的装入1、打开Mapgis的主菜单，点击图形处理的输入编辑，在Mapgis编辑器中选择新建工程，确定后，进入Mapgis编辑子系统。
2、在“设置”菜单中，修改目录环境，使得以后所保存的项目都在某个预设的文件夹中。
3、在“v矢量化”菜单下，点击“装入光栅文件”，打开预先扫描好的地质底图（tif格式）。
二、制作坐标图框1、记下地图的比例尺，底图需要勾绘范围在工程图框中的对应坐标。
2、在Mapgis的主菜单，点击“实用服务”菜单下的“投影变换”。
3、在Mapgis投影变换系统中，点击自选图框进入“矩形图框参数输入”窗口。
在“矩形分幅方法”中选择“用户自定义坐标框”，在“网线类型中”选择“绘制实线坐标线”，根据勾图需要的坐标范围输入有关坐标起点、终点、间隔值、线宽、字体、文件名等参数，确定后计算机会自动生成需要的坐标图框。
然后保存。
4、如果所引入的坐标框与底图之间存在位置误差，可通过Mapgis 编辑子系统中“其他”菜单下的整体移动进行调整。
三、在工程图中引入坐标框在已打开地质底图的Mapgis编辑子系统中，于文件名底下的空白处右击鼠标，点击添加项目，将先前保存的坐标框点、线、面文件引入工程系统中。
四、点项目文件编辑所有的点图标、文字或数字标注都属于点文件。
1、在“文件名”下面的空栏中右击鼠标，后点击“新建点”，输入点文件名，确认。
2、在所要编辑的点文件前的小方框里打勾。
3、在“点编辑”菜单下或在图标栏中选择“输入点元”，进入“输入点图形”窗口，选择所要输入的类型，即注释或子图，确定后再输入相关的类型、大小、颜色等参数，确定后即完成一个点的点编辑。
4、点元输入后，如果位置、角度、大小、颜色和图案等参数不合适或需要改变，可在“点编辑”菜单下或在图标栏中选择相应的编辑工具进行调整。
5、某点元各项参数确定后，先点击需要保存的项目使其颜色变蓝，再点击鼠标右键，后点击“保存项目”保存该编辑项目。
化探数据处理成图的过程毕武1、2段新力1、2黄显义1、2袁小龙1、2彭仲秋1、2李永华1、21.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，8300912.新疆地矿局物化探大队计算中心，新疆，昌吉，8311000 前言GeoIPAS软件用户群不断扩大，由于各用户对系统的熟悉程度不同，对软件功能了解不够，有必要分专题将GeoIPAS处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。
1 处理步骤化探处理的成果包括：(1)参数统计表；(2)R型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。
在GeoIPAS系统中，化探数据处理分为以下几个主要步骤：1.1 数据检查数值检查，坐标检查，重复样坐标检查。
1.2 分析处理重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。
1.3 数据分析数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。
1.4 网格化离散数据网格化、XYZ数据转网格数据1.5 成图点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。
1.6 单元素异常参数统计1.7 综合异常登记卡图2 化探数据处理成图流程2 具体处理过程2.1 数据检查我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X坐标，不加带号，南北向纵坐标为Y坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、0、>等字符，如果有这样的情况要找实验室给予纠正。
数据准备好后，我们要把数据转换成TXT后缀的文本文件，这就做好了处理前的准备数据工作。
2.2 分析处理2.2.1 重复样三层套合方差分析一般是从分析样品的结果中挑出来重复样的分析值，每组四个样品，按如下顺序排列：11 第一次采样第一次分析12 第一次采样第二次分析21 第二次采样第一次分析22 第二次采样第二次分析结果：三层套合方差分析成果--元素:Cu三层套合方差分析成果--元素:Au查表：Fa0.05(18,19)=1.35 Fb0.05(19,38)=1.22Cu：11.152 > 1.35 1.482 > 1.22 即：F1>Fa F2>FbAu：10.781 > 1.35 0.483 < 1.22 即：F1>Fa F2<Fb说明：一般是F1＞临界值，F2＜临界值。
Fa>f0.05()理论值,说明采样与分析误差不显著，此次分析结果可靠，Fb>f0.05()理论值,说明采样与分析误差显著，且以采样误差为主。
例如报告中的描述：先列一个表格说明所有元素的F1和F2；重复样三层套合方差分析，39种元素F1值远超过临界值，说明采样与分析误差变化均低于区内地球化学自然变化，数据是可用的。
F2值表现为两种情况，Ag、Pb、P、U、MgO、La等6个元素F2值小于临界值，说明采样误差不大于分析误差，以分析误差为主；其余33种元素F2值均大于临界值，说明误差应以采样误差为主。
2.2.2 地层编码统计从样品相对应的地质图中自动读取地层的颜色号作为编码，实现每个样品的地层自动编码。
在GeoIPAS系统中的化探->化探特征参数统计->地层编码统计其处理结果是在输入的数据文件中增加一列，过程是将坐标投影到地质图相应的投影参数上，以相同地层颜色来套取数据，取地层的颜色号为落入此地层的采样点的地层编码，所以要想一次做好，要先对地质图做个预处理，最好是按编码顺序，将颜色号改成相应的顺序号，比如Q的颜色改为1，N的颜色改为2，E 的颜色改为3，等等。
相对应的出来的结果就是按顺序的地层编码了，可以直接用了，绘制的直方图也是这个顺序了。
2.2.3 化探特征参数统计化探参数统计成果--元素:Cu注：浓集克拉克值可以选用系统提供默认的，也可以直接输入用户所在不同地区的浓集克拉克值。
这个值可以保存为文件，以后导入使用。
2.2.4 R型聚类分析从谱系图中可以看到，Cr、Ni、Co相关性最好，我们在做组合异常图时可以参考谱系图来做元素的组合，但这个也不能完全作为组合参考，我们通常还是以常用组合方式来做组合异常图，如：Cu-Pb-Zn；W-Sn-Mo；Cr-Ni-Co；Au-Ag-As-Sb。
2.3 网格化2.3.1 工作比例尺与测网密度DZ／T 0145-94土壤地球化学测量规范从上表可以看到，化探在实际工作中一般和物探工作一起工作的常采用矩形网格，但在实际成图中却是用正方形网格的好，不用经过网格化处理，没有数据损失，也符合绘图的正方形网。
1:20万化探在全国范围已基本完成，只在新疆等高山区有少量的没有覆盖全，最常采用的是组合样，即：1点／4km2，少量的有单点样的，这主要用转网格数据的方式来将数据做网格化处理。
单点样的网格化一般采用500×500到1000×1000的网格距，搜索半径3000×3000到5000×5000。
组合样做转网格数据，注意起始坐标为奇数公里，间距为2000，然后要进行一次补空扩边处理，这主要是处理数据成图到图框边。
1:5万化探是现在化探普查的主要工作，一般在山区采用自由网的水系样，土壤测量一般用正方形网格采样，1:1万化探作为1:5万化探异常检查的主要工作手段，常采用100×40的矩形网格，少量有100×20、200×40、150×40等加密或者抽稀采样。
也有采用物探布设的测网非南北方向的测线。
2.3.2 下面用几个问题来具体说明。
问题1：1:5万土壤测量，网距330×300的怎么设置搜索半径和间距？ 5万化探，网距500*500的，用GeoIPAS做地球化学图的时候，怎么设置搜索半径?答：a、关于网格距，最好采用正方形网格，网格距为野外采样的点线距的一半，比如野外采样为500×500，网格化的网格距选250×250就可以了，100×40的选为50×40或者50×50。
b、网格化的搜索半径选择需要一定的经验，建议你首先看看点位数据图，选择一个既能将所有数据联系起来，又不至于网格后把未采样区填充的参数，一般是网格距的2-3倍。
问题2：在进行1:5万化探采样时，乙类异常区采样密度设计为12个点／km2，丙、丁类异常区及其它基岩区采样密度设计为个8个点／km2，网格化时X Y 坐标的间距、行列数如何选？答：一般1:5万化探点线距250×250米或者500×500米；搜索半径以满足丙、丁类异常区网格后只有个别空白区或者无空区，搜索半径一般1000×1000到2000×2000米；如果对网格化参数设置不太清楚，可以选择“计算网格化参数”；网格化后如果异常比较细碎零乱，可用数据预处理→二维数据滤波处理→周围点平均圆滑处理，X和Y方向点数选择3，原点倍数选择2-4倍。
问题3：请教如下问题． 1:2.5万土壤测量，线距250m ，点距40怎么设置搜索半径和间距？答：搜索半径500*500，间距40*40或者50*40，注意最好与原始点位重叠.问题4：1:1万化探，南北向测线，100×40网格距。
答：这样的数据本身符合网格数据的格式，只是以TXT文本格式保存，所以我们可以借助GeoIPAS软件的“数据预处理”→“数据格式转换”→“XYZ数据转为网格数据”功能，这里注意数据坐标的起始值、终止值和网格间距，网格间距就选择100×40，起始值和终止值要保证数据的原始点位都与网格点重合。
这里有一个问题，就是有些项目用的是GPS坐标，不一定是与网格点位完全重合，一般我们建议用规则网的坐标。
圆滑方法：周围点平均圆滑；圆滑参数：X方向点数：3，Y方向点数：3，原点倍数：3问题5：我的测网不是正南北向是45度方向，用咱们的软件进行网格化后变成正南北的，而且图幅形状不是原来测网的形状。
请问这样的问题怎样解决?1:1万化探，北西向测线，100×40网格距。
答：这样的数据虽然也是网格数据，但是经过旋转的测线方向与我们GRD 所定义的网格间距是有夹角的，最常采用的是离散数据网格化方法。
离散数据网格化选择幂指数加权，指数因子：10；搜索范围：圆域，300米，网格化间距：40×40。
要用XYZ数据转网格数据的方式，就要知道测线的方位角，采样的点线号，起始坐标；利用点线号转换为实际点线距，然后再把点线距转为网格数据，成图后，统一处理，包括旋转方位角，然后平移到图上的起始坐标位置套图框。
2.3.3网格化成图的比较：三角剖分绘制的地化图，有被测线方向拉长的趋势，虽然其能保留原占坐标位的原始值，但成图结果不是很美观实用。
距离平方导数加权和幂指数加权网格化，对数据的压制较大，异常线的形态较刻板。
Kriging网格化选用了两种搜索范围，一种是圆域：300×300；一种是椭圆域：300×100，并且旋转45度。
这两种方法整体上差异不大，只在局部的空区上表现出来不同，对数据的压制较小，异常的形态也保持较好，也是常选用的网格化方法，但其也有一个问题，就是对特高值边上的空区会处理出负值，这个不是我们希望得到的结果。
2.3.4经验：在处理1:1万化探数据时，可以根据测线方向和数据坐标点位的不同选择不同的数据处理方法，要遵循的一个原则是尽量使网格化后的坐标点位与原始数据坐标点位接近，网格化搜索的范围越小数据损失越小。
对于网格间距以小于等于最小的原始点距为准，如果要保留更多的细节，也可以将网格距减小一倍，如果图面太过零乱，可以用GeoIPAS的二维滤波处理的周围点平均圆滑处理，这里可以加大原点倍数，默认为2，可以选择3-5倍。
所有这些处理方法可以先选择一个元素用不同的方法做出结果来比较一下确定效果最好的方法，记住处理参数然后再统一批处理。
2.4 成图2.4.1图件编制的要求：a、土壤地球化学测量图件分两部分：基础图与推断解释图。
b、区域调查和普查工作的图件编制按化探区域调查有关规范执行。
c、图件编制必须符合地球化学勘查图式图例及用色标准规定。
d、成果报告需做交通测区位置图、实际材料图、等值线图、综合异常图及其它推断解释图。
