化探数据处理成图的过程
毕武1、2 段新力1、2 黄显义1、2 袁小龙1、2 彭仲秋1、2 李永华1、2
1.乌鲁木齐金维图文信息科技有限公司，新疆，乌鲁木齐，830091
2.新疆地矿局物化探大队计算中心，新疆，昌吉，831100
0 前言
GeoIPAS软件用户群不断扩大，由于各用户对系统的熟悉程度不同，对软件功能了解不够，有必要分专题将GeoIPAS处理数据及成图过程做一系统总结，下面就化探数据处理成图的过程做一总结。
数据输入 数据预处理 数据检查 数据整理 原测数据 投影变换 数据网格化 格式转换 数据处理 重磁数据处理 电法数据处理 化探数据处理 数学地质数据处理 水文地质数据处理 图件绘制 点位数据图 点位符号图 剖面图 输出 打印输出 图像输出 文本文件 图形文件 …… 等值线图 场强图 三元图 三维立体图 图1 金维地学信息处理研究应用系统(GeoIPAS)图件制作流程 1 处理步骤
化探处理的成果包括：(1)参数统计表；(2)R型聚类分析-谱系图；(3)重复样三层套合方差分析或者重复样合格率计算结果；(4)点位数据图；(5)地球化学图；(6)直方图；(7)组合异常图；(8)综合异常图；(9)远景区划图；(10)单元素异常参数统计(附表册)；(11)异常剖析(附图册)；(12)综合异常登记卡(附表册)。
在GeoIPAS系统中，化探数据处理分为以下几个主要步骤： 1.1 数据检查
数值检查，坐标检查，重复样坐标检查。 1.2 分析处理
重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。
1.3 数据分析
数据变换；衬值、累加衬值；数据累加、累乘、比值；异常归一化。
1
1.4 网格化
离散数据网格化、XYZ数据转网格数据 1.5 成图
点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。
1.6 单元素异常参数统计 1.7 综合异常登记卡
开始 数据输入 分析样品与坐标对应 TXT 数据预处理 TXT 数值检查 坐标检查 重复样提取 采样编码统计 三层套合方差分析或TXT 数据分析处理 TXT 矢量 重复样合格率计算 化探特征参数统计 R型聚类分析 TXT 网格化 GRD TXT 离散数据网格化 XYZ数据转网格数据 点位数据图 直方图 成图 地球化学图 GRD 组合异常图 单元素异常图 综合异常图 剖析图 TXT 矢量 异常统计 CSV 结束 图2 化探数据处理成图流程
单元素异常参数统计 矢量
2
2 具体处理过程
2.1 数据检查
我们的数据处理工作从化验室提供的样品分析报告开始，项目要提供坐标和样品对应的分析数据，坐标我们一般取实际工作中的米单位，系统中默认东西向横坐标为X坐标，不加带号，南北向纵坐标为Y坐标，需要时还要提供样品对应的地质编码，我们拿到这个数据后首先进行数据检查，以确保数据中不出现写错、漏填、负数、0、>等字符，如果有这样的情况要找实验室给予纠正。数据准备好后，我们要把数据转换成TXT后缀的文本文件，这就做好了处理前的准备数据工作。
2.2 分析处理
2.2.1 重复样三层套合方差分析
一般是从分析样品的结果中挑出来重复样的分析值，每组四个样品，按如下顺序排列：
11 第一次采样第一次分析 12 第一次采样第二次分析 21 第二次采样第一次分析 22 第二次采样第二次分析 结果：
三层套合方差分析成果--元素:Cu
方差来源 地点间 样品间 分析间 离差平方和 13.224 1.252 1.689 自由度 18 19 38 平均离差平方和 0.735 0.066 0.044 F-检验值 11.152 1.482 三层套合方差分析成果--元素:Au
方差来源 地点间 样品间 分析间 离差平方和 0.573 0.056 0.232 自由度 18 19 38 平均离差平方和 0.032 0.003 0.006 F-检验值 10.781 0.483 查表：Fa0.05(18,19)=1.35 Fb0.05(19,38)=1.22 Cu：11.152 > 1.35 1.482 > 1.22 即：F1>Fa F2>Fb Au：10.781 > 1.35 0.483 < 1.22 即：F1>Fa F2 3
化探数据处理成图过程
