地球化学化探数据处理与成图

用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．RE计算RE=abs（C1-C2）/（C1+C2）×200要求：小于3倍检出限时，RE≤85%，大于3倍检出限时，RE≤66%为合格。

总合格率一般要求大于70%。

需要说明的是：这种方法适用于简单对比，就是一个点取2个样时使用。

如果用三重套合分析，就不能用这种方法了。

所以写设计时，一定要用简单对比方法。

还有，这些年的化探中的RE还没有超标的，如果真的超标了，我也不知道怎样处理。

二、数据处理：1．剔除一级及二级标准样；2．剔除重复样；3．剔除0值及化验室输入错误的值，或小于检出限的值。

主要是指≤0.3这样数据的≤号；4．做原始数据图；5．计算异常下限，做单元素异常图，圈定单元素异常。

单点异常，只有外带的不圈，有中、内带的圈出。

外带用黄色，中带用浅红色，内带用深红色。

单元素异常编号为Au-1，Ag-1等。

需要指出的是，圈定异常时应该形成数据异常图，但交报告时，必须分开，就是形成一张数据图，再形成一张异常图；需要说明的是，如果面积较大（这个没有标准），总体说是水系面积超过一个5万图幅，就要分子区计算下限。

分子区的原则是不同年代、不同地质体都要划分成子区。

如果各个子区的异常下限接近，就采用总的，否则分别确定子区的异常下限，然后分别圈定子区的异常。

6．填单元素异常评序表。

异常点数、面积、平均值、极大值、标准离差、衬度、规模、浓度分带等。

评序有5参数和多参数两个评序，5参数不参与的参数有标准离差、浓度分带及异常点数。

需要说明的是，排序时，单项值高的给1，以下类推2、3等，一样的值给一样的排序。

不产生空的值。

浓度分带有内带的给1，中带的给2，外带的给3；7．做组合异常图，先把所有的元素做成1张组合异常图，只要外带，圈定组合异常。

内蒙古扎赉特旗东芒合矿和哈拉街吐矿化探数据处理及图件编制方法1 化探数据质量评价的数据处理(分矿区)⑴统计重采样和重分析抽查样所占样品总数的比例比例 = （重采样和重分析抽查样数/工作样总数）100%⑵作出SSPS数据文件将重采样和重分析样分别作成SSPS数据文件。

文件中列出项目为：①重采抽查样重采样号元素含量相应的工作样号元素含量②重分析抽查样重分析样号元素含量相应的工作样号元素含量⑶计算各元素相对误差重采样和重分析抽查样相对误差均按RE(%) = |C1-C2|/0.5×(C1+C2)×100%计算。

C1为重采样或重分析抽查样的分析含量C2为重采样或重分析抽查样的相应的工作样的分析含量| |为绝对值RE(%)≤30%为合格，＞30为超差（不合格）；（Au：RE(%)≤50%为合格，＞50为超差）⑷计算各元素的合格率η= （抽查样品中合格的样品数/抽查样品的总数）100%合格率（η）应＞80%，即这批样品的分析结果是可信的。

⑸列表表示检查或分析质量结果表××化探重采样抽查各元素的合格率(%)Cu Pb Zn Cr Ni Co Sn V Ag Ti2 矿区地球化学特征研究的数据处理(以哈拉街吐为例)⑴作出SSPS数据文件作出下列SSPS数据文件：①文件1：整个矿区数据文件；②文件2：矿区地层数据文件；③文件3：矿区岩浆岩数据文件；④文件4 ：下二叠统大石寨组（P1d）数据文件；⑤文件5 ：下白垩统大磨拐河含煤组（K1d）数据文件；⑥文件6 ：华力西晚期侵入岩数据文件；⑦文件7 ：燕山期早期侵入岩数据文件；⑧文件8 ：燕山期晚期侵入岩数据文件；⑨文件9：已知矿附近一定范围数据文件每一数据文件的内容项目包括：序号野外号 X坐标 Y坐标各元素的含量⑵整个矿区和各地质单元（各地层、各岩浆岩）样品各元素含量特征统计统计的参数包括：①元素含量平均值；②最大值；③最小值；④标准离差；⑤变化系数（标准离差/含量平均值）；⑥浓度克拉克值（元素含量平均值/该元素的克拉克值）整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。

1数据准备（1）成图前逐点采集样点的平面坐标。

数据排列顺序为x坐标（米）、y坐标（米）、元素含量或采样点线号，图面坐标与实际坐标要进行转换。

分别建立各元素含量数据及点位点线号的电子表格文件，在Excel电子表格中另存文件，文件类型为.csv,逗号分隔（*.csv）。

（2）用记事本软件将*.csv文件打开，在首行加入notgrid,将文件类型"*.csv"改为"*.det"保存。

2取样点位图成图（1）进入MapGIS6.5主菜单，点击空间分析—DTM分析。

（2）进入文件菜单打开三角剖分文件…，即上述的*.det文件。

（3）进入模型应用菜单下选择高程点标注制图，弹出下列菜单。

（4）高程点标注显示框中的“符号”是指采样点位置符号，可以点击缺省符号区域自定义为“ｏ”符号；标注是指元素点、线号。

（5）在标注位置框中设置数据的“对齐方式”选择“右边”，“角度”为“0”。

也可以自定义设置x、y位移量。

（6）在符号尺寸框中一般选择固定尺寸，为0.5～1。

（7）在“标注字体…”框中设置字体大小、颜色、类型。

（8）在“标注格式…”框中设置小数位数。

（9）设置完成点击“确认”。

（10）保存点文件、线文件。

3单元素异常图绘制（1）打开单元素高程数据三角剖分文件…，即单元素*.det文件。

（2）在TIM模型菜单下，进行快速生成三角剖分网，对于弹出的“*.tin”数据已修改是否先保存，点击否。

利用删除三角剖分网删除不必要的三角网。

（3）在TIN模型菜单下追踪剖分等值线。

（4）在设置等值线参数对话框中以前面计算的各元素异常下限为等值线外层，选择等值线套区、保留边界线、等值线光滑处理—高程度、轴向标尺。

（5）在等值线定层中设置所确定的等值线间隔、线参数类型、注记参数类型及修改区参数。

（6）在所有参数设置完毕，则点击确定，异常等值线自动生成。

（7）保存所生成的三角剖分文件、点文件、线文件、区文件。

化探数据处理及图件编制第二节分析方法及质量评述一、分析方法本次扫面和异常查证的全部样品均交由四川省地矿局华阳地矿检测中心测试，根据任务书要求共分析测试元素14种。

样品从加工到测试到质量监控均按中华人民共和国地质矿产行业标准DZ0130-1994《地质矿产实验室测试质量管理规范》、ISSN-1870《1?5万区域地质调查及地球化学样品分析方法及质量管理指导性规程》和2002年新疆地勘局试验管理科《1?5万化探样品分析质量过程管理规则报告》进行。

14种元素的分析方法见表3,3。

二、技术要求1、报出率十四种元素的总报出率应大于95%。

2、外检样对已测试样品，测试单位按照3%的比率进行外检。

3、分析质量检查及质量监控方案为了有重点地监控元素的分析质量，实验室在送样单位确定的分析元素中，要再选择若干种主要监控元素并根据这些元素在本省制备的全部GRS二级标样中选择四个在元素含量范围及基体组成均为合适的GRD二级标样作为本图幅质量检查监控之用。

主要监控元素和二级标样的选择均应和送样单位协商进行。

每一大批样品测定完毕后，应将数据交给质量管理人员，对每一小批中插入的四个二级标样及四个重复分析(内部检查)样进行统计计算，并及时绘制日常质量监控图，在日常金的分析工作中，必须进行不小于10%的内检抽查。

为满足在一个较大范围的成矿远景区带内的1?5万图幅的拼接，应对分析的准确度进行检查和考核，为此实验室应在每一个1?5万普查化探项目完成后，分析8个GSD一级标样一次，痕金分析也应用金标样作准确度检查。

准确度和精密度计算结果应符合表3,2的要求。

4、微量金由于金元素在自然界中的均匀度和赋存状态对分析检测影响比较大，为确保金元素的分析质量，化验室特采用两种监控措施:第一，在每一分析批次的50个样品中插入两个国家?级标准物质GBW系列，用以计算实测值与推荐值之间的对数偏差:ΔlgC,lgC,lgC; 定值实测值第二，该地区样品分析结果结束或阶段性结束后，再对高、低异常点进行随机抽样检查约20%.5、?级标样为严格监控各元素的分析质量，实验室选取了四个不同含量的GRD系列监控样，每批次50个样品密码插入一组，与样品同时分析。

化探数据处理成图的过程毕武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、>等字符，如果有这样的情况要找实验室给予纠正。

最新最详细化探数据处理与编图流程⼀、指导思想成矿地质背景地球化学研究就是从地球化学特征出发，借助已建⽴的地球化学信息提取技术，充分利⽤地球化学调查所获得的海量数据信息，提取有关反应成矿地质背景条件的地球化学信息，并编制相应地球化学图及相应的推断解释图件，为资源潜⼒评价有关成矿地质背景的研究提供地球化学⽀撑。

⼆、⼯作内容（⼀）基础图件成矿地质背景条件的地球化学信息提取⾸先是要编制有关基础地球化学图件。

主要有：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。

首先把EXCEL表Y坐标放第一列，另存为CSV格式。

1：50000地化图绘制：1、用记事本打开Ag.csv格式数据，将x（坐标符号），y（坐标符号），Ag （元素符号）改为notgrid，另存为Ag.det格式。

对于1：5万化探，x、y单位为m，成图比例尺为1：50000万，那么将x、y同除以50，就变成mm。

2、依次打开MAPGIS程序主菜单\空间分析\DTM分析（MAPGIS数字地面模型自系统），打开Ag.det文件。

3、打开Tin模型\快速生成三角剖分网。

4、打开Tin模型\删除三角剖分网边。

5、打开Tin模型\三角网内插网格化：左下角坐标设为最小整公里坐标，右上角坐标设为最大整公里坐标；1：50000地化图网格化间距实验证明设为X：125，Y：125最为适当，另存栅格化结果文件Ag.Grd.6、打开Grd模型\平面等值线图绘制：选择：等值线套区、绘制色阶、保留边界线、等值线光滑处理（高程度）、示坡线、线型。

线参数设定：默认。

注记参数设定：标注字体尺寸：2mm，小数部分位数：与该元素线值位数相同，其它默认。

(.VEL)。

以后重新绘图时，可直接装入该文件。

7、确定成图，保存点、线、面文件。

8、极值点标注：等值线图极值点标注，只能在“三角剖分-追踪剖分等值线后”，在“模型应用-高程点标注制图”子菜单下标注。

在标注输出设置中，选择“等值区内部极值点标注”，然后按自己的要求设网上收集资料一、建立数据表1、用电子表格建立成图元素表（不要表头，建议将平面坐标按Y、X、数据…列次序移到表前1、2…列，Y、X坐标分别为7位、6位整数）；2、另存为文本文件（制表符分隔）类型的新文件并取名(如文件名为茨竹箐原始数据表)。

二、Mapgiss投影转换1、打开Mapgiss程序主菜单→点击″实用服务″→点击″投影变换″→点击″P投影转换″→点击″U用户文件投影转换″→打开文本文件(茨竹箐原始数据表)；2、选″指定分隔符″→指定坐标X、Y列（注意X、Y列可能倒置）；3、点击″设置分隔符″→选″Tab″及″空格″→点击″数据类型″→将所有的字符串点为″双精度″→在″属性名称所在行″选框中选择数据（注意数据不能选″无″）→在″加入″列选定所需的坐标及元素值→确定；4、点击″点图元参数″→选择子图参数（建议选编号35的圆点的子图号，高度宽度一般均为1，颜色为1）→确定；5、选″不需要投影″→点击″数据生成″→确定→关闭；6、另取点文件名（如茨竹箐数据投影)→保存→退出。

地球化学实习主要流程步骤实习四概率分布型式的检验及背景值和异常下限的确定1.将实习用的数据新建一个文件夹，以自己的学号命名2.复制表4-1的数据，选择性粘贴→文本，然后保存。

文件→另存为→一共3个要点：保存的位置：自己建的文件夹里，必须和运行程序在同一级目录；文件名：推荐4-1；保存的文件类型：csv逗号分隔 (后缀必须是csv格式)；然后以同样的操作方式复制表4-2的数据进行相同操作。

3.运行BASIC程序，得出结果4.以写字板方式（EXCEL）打开运行结果.5.简单解释一下相关参数的意思6.以相同的方式运行4-2的数据，得出运行结果.7.演示做作业。

复制实习4的作业到另外一个空白的WORD文档里，单独以自己的学号命名。

具体方式为：姓名+学号+实习内容的顺序。

比如：200701010326朱国勋No.18.开始填表：首先求出均值和方差，用各组组上限值（Xi），计算，以4-1为例：均值为1.46，方差为0.26，用EXCEL编辑公式=ROUND((A2-1.46)/0.26,2)9. 据值查标准正态概率分布函数表，得到各Xi点的累积概率F（）。

10.计算Dn值，在F（）的基础上加减13.6分组间隔频数f频率f累计频率∑fμ值概率F(μ)%临界曲线（上）临界曲线（下）0.77-0.97 3 3 3 -1.88 3.01 16.61 -10.590.97-1.17 9 9 12 -1.12 13.14 26.74 -0.461.17-1.37 22 22 34 -0.35 36.32 49.92 22.721.37-1.57 32 32 66 0.42 66.28 79.88 52.681.57-1.77 20 20 86 1.19 88.3 101.9 74.71.77-1.97 12 12 98 1.96 97.5 111.1 83.91.97-2.17 2 2 100 2.73 99.68 113.28 86.0811.建立画图的文件。

地球化学探矿工作方法及注意事项异常下限计算公式：都采用Excle表格计算(原始数据参照表三)平均值：=A VERAGE(C2:C36)方差：STDEV(C2:C36)异常下限：平均值+K值*方差计算过程中自动剔除大于平均值加3倍离差的高值和小于平均值减3倍离差的低值。

K值自选（1-3）剔除次数的多少取决于没有大于和小于这些数据。

表一表二表三等值线图生成方法步骤：利用Excel先做一个三角部分文件，格式如图：A1：要打上notgrid.A2为生成点文件的X坐标 B2为Y坐标 C2为高程代表化学分析结果。

文件需要另存格式如图：点击保存，关闭表格。

打开MAPGIS6.7空间分析-DTM分析-点击文件打开三角部分文件。

如图：文件打开后在菜单栏中找到Grd模型打开里面有一栏离散数据网格化打开。

如图：网格化方法用：Kring泛克立格法网格化，文件换名保存文件。

“选择打开”，选三阶漂移如图：其它内容不用改动，改完后按确定。

“搜索打开”搜索类型改为八方向，其它内容不用改动，改完后按确定。

如图：全部改完后按确定。

再一次在Grd模型里找到平面等值线图绘制点击后打开离散数据网格化后保存的文件，如图：把文件打开如图：设置：等值线套区、绘制色阶、保留边界线全部选上。

光滑：等值线光滑处理选上，光滑度选择高程度。

等值线定层一般定四层，第一层为背景值，第二层为异常下限值，第三层为中带值（异常下限值的2倍），第四层为内带值（异常下限值的4倍）。

注记参数都改为YES。

区、线参数参照图例规范设定。

制图幅面：改为数据投影变换如图：原数据投影参数与结果数据投影参数都要改变，设置结果与投影点位图一致。

全部设置完后点击确定。

生成图像如图：生成图像后保存点、线、区。

保存好以后再一次打开MAPGIS6.7图形处理-输入编辑打开刚才已保存好的点、线、区文件进行图式修饰。

点位图制作方法步骤：利用Excel表格把所有化探数据先处理一下如图：注：ID代表总序号、YH代表样号、Y代表X轴坐标、X代表Y轴坐标、Au、Ag 代表元素化学分析结果。

地球化学+地质填图常用教程总结一、1:5万土壤采样点位设计(水产沟幅为例)1、标准图框生成首先利用Mapgis—实用服务—投影变换—系列标准图框—生成1:5万图框.图框模式：高斯坐标实线公里网.投影参数：起点纬度491000. 起点经度1210000. 网间间距1KM.图框文件名：水产沟图框.图框参数输入：图框内容：图幅名称：水产沟幅.图框参数选择：输入并绘制接图表+绘制图幅比例尺+图框外图廓线.2、图形校正图片JPG转换成msi文件： Mapgis—图像处理—图像分析—文件—数据输入—添加文件(图片JPG)—转换—转化成msi文件.校正到标准图框内：Mapgis—图像处理—图像分析—文件—打开影像—msi文件—镶嵌融合--打开参照文件—参照线文件(标准框). 镶嵌融合—添加控制点—(先点击带校正图片某点，然后点击对应的标准框点，最后按两次空格键，选择是，以此类推，尽量多的找几处控制点，通常不少于13个控制点，要求残差<1为准)—镶嵌融合—校正预览—透明度调节--回车--镶嵌融合—影像几何校正—几何校正--重命名保存文件.3、图片数字化对地形线、河流、公路、村庄等地形要素进行数字化过程. 同时，准备大格号、图例、网格线、采样位置、重复样和监控样分析表格等文件.4、采样点位设计四幅图框：据前述方法，生成5万四幅合体的图框（注意修饰）。

将四幅图之间的界线用红线标出，根据阵列复制线，画出网格线，在网格线的交汇处定出大格号标注，然后建立采样点位文件，在上面点击采样点位置，辅助工具1—导入导出功能—点位置转属性. 辅助工具2，水系沉积物自动编号，0001大格号的左上角坐标为5488000,355000，一行上的大格号从0001—0004，点击确定，重新显示文件，出来小格号标注文件.5、重复样、监控样在每个单幅的采样点位图中，设置重复样和监控样，按从上到下，从左到右的原则，即a—c—b —d，查出每50个样品，其中1个是重复样，4个是监控样，45个是普通样品，导出excel表格之后排序，对其重复样和监控样设置，然后在单幅采样点位图中标出来.二、1:1万土壤采样点位设计1、基点基线准备在地形图基础上，可以在图框文件里输入一个点，作为基点（待设计完毕可以删除，例如基点为200/100，点号为200，线号为100），图框外部边缘作为基线。