地球化学涉及的基本图件制作及相关软件应用(点位图)

GeoPlot地球化学数据投图软件-使用说明GeoPlot地球化学数据投图软件-使用说明第一章：介绍1.1 软件背景与目的1.2 目标用户1.3 安装要求第二章：安装与配置2.1 软件包2.2 安装软件2.3 创建用户账户2.4 配置软件参数第三章：界面概览3.1 主界面3.2 工具栏3.3 菜单选项3.4 数据显示界面3.5 设置面板第四章：数据导入与处理4.1 导入数据文件4.2 数据清洗与预处理4.3 数据格式转换4.4 数据筛选与过滤4.5 数据计算与统计第五章：地球化学数据投图5.1 绘制散点图5.2 绘制地球化学图5.3 绘制剖面图5.4 绘制等值图5.5 绘制空间分布图第六章：图形设置与编辑6.1 修改图形样式6.2 添加图例与注释6.3 调整图像大小与比例6.4 图形导出与保存第七章：其他功能7.1 数据文件管理7.2 数据交互与共享7.3 扩展功能与插件第八章：常见问题解答8.1 如何导入多个数据文件？8.2 数据清洗过程中如何删除异常值？8.3 如何标准化图形？8.4 地球化学图如何进行符号设置？第九章：技术支持与反馈9.1 联系我们9.2 反馈建议9.3 常见问题咨询附件：附件一、安装包附件二、示例数据文件法律名词及注释：1、软件包：指GeoPlot地球化学数据投图软件的安装文件，包括执行文件和相关依赖文件。

2、安装软件：将软件包中的执行文件或脚本运行，并按照向导进行必要的安装配置的过程。

3、用户账户：为了保存用户个性化的设置和数据，软件需要创建用户账户来区分不同用户。

4、数据清洗：指对导入的数据进行去除异常值、填补缺失值等处理的过程。

5、数据格式转换：指将数据从一种格式转换为另一种格式的操作，例如从CSV格式转换为Excel格式。

6、数据筛选与过滤：指根据特定条件对数据进行过滤，例如筛选出特定元素浓度大于某个值的样点。

7、数据计算与统计：指对数据进行数学计算和统计分析的操作，例如计算元素比值、计算数据的平均值和标准差等。

利用Arcgis软件制作土壤等地球化学系列图
软件要求：Arcgis软件等。

数据要求
准备好删除重复样的原始分析测试数据。

第一步：把数据投点进入Arcgis
把准备好的数据，内容包含横坐标、纵坐标、样号及元素分析数据等，投点为点图层。

第二步：处理
1.打开工具箱，找到“Spatial Analyst”-“插值分析”-“反距离权重法”
2.在弹出的对话框中，第一行输入已经生成好的“点图层”；
3.第二行选择需要分析的地球化学元素；
4.输出栅格为生成shape的文件存储地；
5.搜索半径选择“固定”，半径大小根据需要自行确定；
6.点击下方的“环境”，弹出“环境设置”对话框，找到“处理范围”，在捕捉栅格中选择工作区的范围的面图层；
7.在栅格分析中输入指定的像元大小（图1中“6”对应的输出像元大小的值）；
8.“掩膜”仍然选择工作区范围的面图层；
9.点击“确定”。

图1
图2
第三步：设置
根据元素的需要设置类别（分级）数级分级间隔（即图3中右边的“中断值”），是需要可以勾选“4”中的显示标准差、平均值与否。

图4
说明：Arcgis软件首先需要在“自定义”-“扩展模块”中Spatial Analyst 勾选上，不然不能使用Spatial Analyst。

地球化学测量编图软件介绍目前，可用于地球化学测量等值线圈定的软件有GeoExpl、MAPGIS、SUFER、GEOIPAS等4种。

介绍如下：GeoExpl软件由中国地质调查局发展中心研制。

软件提供了数据预处理功能，提供了滑动平均、最近点、平均值、距离倒数加权、指数加权、累加、克里格插几种网格化的方法。

提供了等值线追踪及相关地球化学标注等一整套生成地球化学等值线图的功能。

但其仅提供了与MAPGIS的数据交换接口，没有矢量数据的编辑功能。

提供了大量地球化学相关的计算工具与地球化学其它图件的生成功能。

该软件在全国的使用较为广泛，软件性价比较高。

MAPGIS软件由中地软件公司研制生产，本次用于测试的版本为6．5。

该软件提供了数据预处理功能，距离幂函数反比加权、泛克里格、稠密数据中值选取、稠密数据高斯距离权等几种网格化的方法，提供了等值线追踪、相关地球化学标注等一整套生成地球化学等值线图的功能。

提供了离散点三角网等值线追踪功能：提供了与主要地理信息软件的数据文件交换接口，提供了强大的矢量数据交互编辑手段。

该软件作为地学矢量图件编辑软件广泛使用，拥有的用户众多。

该软件不提供生成等值线图以外的其它相关地球化学图件的生成功能。

MAPGIS是一个功能强大的制图软件，在数据极差相对较小时，用Kring泛克立格法网格化模型绘制的地球化学图效果较好；数据极差相对较大时，采用距离幂函数反比加权网格化模型绘制的地球化学图效果佳。

两种网格化模型都可以绘制地球化学图。

网格间距的确定一般可以根据工作区的地质背景，只要绘制的地球化学图能客观反映元素的地球化学特征即可。

由于目前地质调查项目要求提交的图形文件为MAPGIS格式，因此在MAPGIS软件下绘制地球化学图是十分必要的。

SUFER软件为GOLDENSOP-TWARE公司研制，本次用于测试的版本为8．0。

该软件提供了数据预处理功能；提供了滑动平均、最近点、最小曲率、改进谢别德、自然点临近、多项式回归、径向基函数、带线性插值三角剖分、移动平均值、数据度量、局部多项式、克里格等多种网格化方法；提供了等值线追踪的方法，提供了样条平滑的网格数据圆滑处理方法；提供了与主要地理信息软件的数据文件交换接口该软件的标注、图框等内容与我国现行地球化学等值线图的要求差别较大。

直接利用ExceI作地球化学微量元素蛛网图、稀土元素配分模式图以及地球化学散点图GeoKit，Geoplot都是很好的地球化学数据处理软件。

但是怎么样直接利用Excel作地球化学微量元素蛛网图、稀土元素配分模式图以及地球化学散点图等常用地球化学图解呢？原理如下：1. 准备好数据：La Ce Pr Hl Sm Em3T.5071 008.3531.50 5.70 1.5花岗闪畏斑岩33,3064 707. 73组30 5. E3 1.三化阳斑岩35.40fiT SO7.^730 10 5 &8 1.4花岗斑岩2B.QD52 30 6.2524. 40 4. &4I 1306.00322 0014 4026 90 3. 060 5夕卡岩8 7415. 70LTS s gg 1. 56 c a强强石□ 370.笑 C. 140 710.23 c.匸—ES-1-1102.1798BS-1-590. 丁3569BS-2-34577773. 297BS-S-5BS-6-S839. 2371BS-&-1O23. 81471,2. 对数据进行标准化，第一排：输入：=样品/球粒陨石，然后回车。

样品和球粒陨石的值用，鼠标去点击相应的格。

如第一个样品，=C2/C9，回车。

L ・ 3N氏IVIrT-nuVf SO T| » « » M 59 5T»H K 5 H s.n4 F 0 Si ? T ： D «lo • 屁阳•畑老 33 30 曲 > p T2 血30 S "El 1 31 5 ］： o. n 4 X. :E ： s2空□ 1» E.7. a 芒 昭帽"皆 3S 箱材.・1 " 10 IB S 世L «S 5汀_ o.?t 4 U *册±结□ M S.H. .4 此闻Bt«r z. MU yJ iv ■* Jfi -H 如嘴啊L JU * 1-4ii • J 4dV E4i 9SU 3a ■■ J v Ja $卡惜M DO 熾軽Mt*2 +析0眉 2 ■0利 I W0 H 1菇 4 |fe faM ______________ _B T4 谒” W1 W4 M_____ ■一9 M1.Xij 门L.护 ........... O 曲4 It _d,li _a ・ 4 2■呻■石 D 7? (fe in0 34辱71° t}0碑Q 11D m0 0Q Tl11 iM n 髀n n«102njROOKl MlMML 44 3A3t 2SH6M! iLStm 17.385«2 13.620B IL 7323 1&34078 LO L S^ IL.SltS 11.441.M 11. MSSfi 轨 THS 斂 皈 3$8新 騎■出扫輛 IL^Si 2爲934阳 14 4«23 值詢倨 13.103€11.33SEa SiS70?410.S6225 10.9^ ］0.»7«2 11.023^海時TH W.W1W St.lTSlt U33471 2CSS87I 1 離綁T 佩的讯 13.44323 1L&48» 9■阳82伸 10L 6<K7m 耶&K10.S2414 1O.T61^ 73-23T 5C 两66f ^6-3503*3fl,317ai20,08659： l£,643€B 13P 5^94||g p 344£39. 0026S 匸够06尅 S r I3ST3B3P J3S7M良 ETQ彌 宜 L 岳测丄ssazrri msm10a i (manr.sMoiis,24475&«JET 挺ingeoasBLMGBIT ^IOMCB E , MZies 1 氏爲nns c WSBHCBE ^SH翳.丄人些亚毗 muzi i 二亦工宜‘瓷巧「［:M 匚.的芒习 u 丄舵裂—工石IT 加 ：叮：rnr JJIOHD 久塑艶二B 見空1二“二刖芒3. 手动计算第一排的数，如上。

地球化学+地质填图常用教程总结一、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），图框外部边缘作为基线。

GeoPlot地球化学数据投图软件-使用说明一．软件介绍GeoPlot软件是用VBA开发并在Excel中使用的，主要用于地球化学数据计算和投图分析的专业宏程序软件；其实用性强、操作简便，可帮助研究人员快速分析地球化学数据。

功能如下：1)投图功能: (1)XY散点图、三角图;(2)岩石地球化学判别图;(3)蛛网图及REE配分图2)数据计算: (1)参数的公式计算;(2)CIPW标准矿物计算;(3)AFC及二元混合计算3)辅助功能: (1)样品分组投图;(2)数据统计结果;(3)数据排序及小数位数;(4)分子量计算;(5)Chart 2Picture;(6)系统设置及数据库管理4)二．安装删除运行自解压安装包程序文件“GeoPlot.exe”，选择一个安装目录（一般安装目录为：C:\Program Files\GeoPlot）就轻松安装成功，如图1。

启动GeoPlot软件的方法有两种：（1）运行安装目录下的宏程序GeoPlot.xla，自动启动Excel并出现GeoPlot的文件菜单和工具栏；打开数据文件后，就可以使用该软件的各项功能；（2）在Excel菜单中加载安装目录下的宏程序GeoPlot.xla，该方法可以让每次Excel运行后，都会自动启动GeoPlot。

如果需要不自动启动GeoPlot则可在Excel下取消宏GeoPlot.xla的加载。

该软件是绿色宏程序软件，所以反安装GeoPlot软件则可直接删除所安装的目录，比如“C:\Program Files\GeoPlot”。

图1 安装GeoPlot的目录选择三．软件使用1．数据准备使用GeoPlot软件进行计算和投图，地球化学数据存放的表格需要复合一个简单的、也是最常用的数据排列格式—元素名称排第一列，样品名称排第一行，其余数据对应元素名称和样品名称排列，如图2。

元素名称常用到软件中查询，所以元素名称中不能有非可见字符（网上复制的数据可能产生不可见字符，常导致不能正确计算和投图，需要用户手动重新输入元素名称一遍）。

地球化学重要图表1、太阳系元素丰度对数对应其原子序数作出曲线图1. H和He是丰度最高的两种元素。

这两种元素的原子几乎占了太阳中全部原子数目的98％。

2. 原子序数较低的范围内，元素丰度随原子序数增大呈指数递减，而在原子序数较大的范围内（Z＞45）各元素丰度值很相近。

3. 原子序数为偶数的元素其丰度大大高于相邻原子序数为奇数的元素。

具有偶数质子数（A）或偶数中子数（N）的核素丰度总是高于具有奇数A或N的核素。

这一规律称为奥多-哈根斯法则，亦即奇偶规律。

4. Li、Be和B具有很低的丰度，属于强亏损的元素，而O和Fe呈现明显的峰，它们是过剩元素。

5.质量数为4的倍数（即α粒子质量的倍数）的核素或同位素具有较高丰度。

例如，4He（Z=2，N＝2）、16O（Z=8，N=8）、40Ca（Z=20，N=20）和140Ce(Z=58,N=82)等都具有较高的丰度。

通过对上述规律的分析，人们认识到在元素丰度与原子结构及元素形成的整个过程有着一定的关系：1. 与元素的原子结构有关。

具有最稳定原子核的元素一般分布最广。

2. 与元素形成的整个过程有关。

H和He丰度占主导地位和Li、Be和B等元素的亏损适于元素的起源和形成的整个过程等方面来分析。

例如，根据恒星合成元素的假说，在恒星高温条件下（n×106），可以发生原子（H原子核）参加的热核反应，最初时刻H的“燃烧”产生He，另外在热核反应过程中Li、Be和B迅速转变为He的同位素42He，为此在太阳系中Li、Be和B等元素丰度偏低的原因可能是恒星热核反应过程中被消耗掉了。

2、金属元素电负性与亲和性的关系3、氧硫元素的性质对比在地球系统中，氧的克拉克值为46%，硫的克拉克值为0.05%，它们都是阴离子中分布最广的两种元素。

硫比氧是同族元素，硫的半径较大，电负性较小，形成离子键化合物的能力减小，而形成共价键化合物的能力增大。

4、格尔施密特的地球元素化学分类5、氧化还原电位对元素迁移的影响（铁离子的Eh-Ph图解）6、元素在水中的存在形式7、风化作用的类型8、元素的迁移能力。

使用Arcgis软件中的Geostatistical Analyst绘制地球化学图
软件要求：Arcgis软件等。

数据要求：准备好删除重复样的原始分析测试数据。

第一步：把数据投点进入Arcgis
把准备好的数据，内容包含横坐标、纵坐标、样号及元素分析数据等，投点为点图层。

第二步：处理
1.打开工具箱，找到“Geostatistical Analyst”-“地统计向导”，在确定性方法中选择“反距离权重法”。

2.在输入数据对话框中，“源数据集”选择已经投进图层的“点图层”，“数据字段”选择需要绘制地球化学图的元素，比如N元素，点击“下一步。

”
3.处理重合样本采用“使用平均值”。

4.在搜索邻域中，最大相邻要素数选择“8”，最小相邻要素数选择“4”，扇区类型选择“4个扇区”，见下图。

设置好后下一步-完成。

第三步：修改属性
右键选择“属性”，把显示质量和缩放优化勾选上，在分类中可以选择分格数，这里以15个分格为例，“方法”选择手动，“类”中选择15，在中断中输入中断的数值即可。

说明：Arcgis软件首先需要在“自定义”-“扩展模块”中Geostatistical Analyst勾选上，不然不能使用Geostatistical Analyst。

直接利用Excel作地球化学微量元素蛛网图、稀土元素配分模式图以及地球化学散点图GeoKit，Geoplot都是很好的地球化学数据处理软件。

但是怎么样直接利用Excel 作地球化学微量元素蛛网图、稀土元素配分模式图以及地球化学散点图等常用地球化学图解呢？原理如下：1.准备好数据：2.对数据进行标准化，第一排：输入：=样品/球粒陨石，然后回车。

样品和球粒陨石的值用，鼠标去点击相应的格。

如第一个样品，=C2/C9，回车。

3.手动计算第一排的数，如上。

然后拖动。

自动计算后面的数据。

如上。

现在开始画图。

如下，【图表向导】→【系列】4.分类（X）轴标志：用鼠标在表中拖动从La-到Lu，然后【确定】。

如上【完成】发现，Y轴是100，200，…这样的。

用鼠标放在Y轴上，右击→坐标轴格式→刻度，最大值1000。

对数刻度。

5.在图中空白出右击，图表选项，填上X,Y轴，用Geokit做的图对比右图是Geokit作的。

与左图一样。

上新的数据就会自动生成图。

按照上面的做法生成微量元素蛛网图：用Excel绘制地球化学散点图记下坐标轴x为35到79。

y为0到16（整数—描图时长点，变成整数）这是文献中常用到的火山岩TAS岩石分类图解。

遇到类似非常帅气的图片，我们可以先用CoreDraw或MapGIS描下来，形成底图，以便日后使用。

然后做XY散点图在散点图中任意右击，选【绘图区格式】→【填充效果】→【图片】→选择制作好的底图→【确定】。

设置坐标，将原图的坐标x为35到79，y为0到16，设置好。

图文均整理自网络，向原作者致敬……。

收稿日期:2007-04-12;修订日期:2007-04-28;作者E-mail:xjrs50@第一作者简介:韩小明(1957-),男,陕西蓝田人,高级工程师,1982年毕业于武汉地质学院,从事地质数据处理工作地球化学等值线图生成方法研究及相关软件使用对比韩小明,张晓梅,王瑞,李晓青,丁礼江,田玉川（新疆维吾尔自治区地质调查院第一区域地质调查所,新疆 乌鲁木齐 830011）摘 要:通过对目前地球化学等值线图的计算机成图方法进行分析,以及主要使用的相关软件对比,指出目前计算机自动生成地球化学等值线图存在的缺点.针对计算机自动生成地球化学等值线图存在的线条不够圆滑的问题进行了研究,提出了进行等值线自动圆滑处理的算法及其解决方案. 关键词:地球化学等值线;矩形网格化;三角网格化地球化学等值线图是地质工作者了解一个地区范围内地球化学情况的重要图件,目前多数是由计算机完成.人们对计算机自动完成的地球化学等值线图总存在一些不满,主要表现在既想保持地球化学分布的细微特性,又想要求地球化学等值线做的尽可能光滑.为此多要进行手工编辑,这项工作耗时又费力,笔者针对该问题进行了一些探索和研究.1 计算机自动生成等值线的主要算法计算机自动生成等值线的算法可以分为2类,即三角网格法和矩形网格法.其中三角网网格法的主要优势是能够保持采样点的确切位置,人为能干预它的结果的手段较少.矩形网格法有多种,每种都会形成略有差异的等值线.主要流程及可干预的手段如下：整理数据 在数据整理阶段一般说来只要能够将数据格式整理为使用软件要求的格式即可,但在特殊情况下可以在本阶段对数据进行干预.开窗口滑动平均处理,该处理也应当算做预处理的一个部分.开窗口的大小以及算法都有不少可选择的余地.总体而论滑动平均是缩小数据间局部差异的的方法.开的窗口越大,最后的等值线结果上显示局部线条越光滑.但是随着局部线条光滑程度的增加,能够反映平面数据的细部特征也被抹杀,至少在大比例尺地球化学等值线图中是这样,这不是我们要看到的结果.数据网格化 目前可使用的数据矩形网格化方法有十几种之多[1,2],主要有距离加权平均法、最近点、平均值、距离倒数加权、指数加权、累加、克里格插、泛克里格、趋势面等.一般在选择网格化方法的同时还要选择网格间距、网格起始坐标及其它的相关参数.应该说网格化的方法选择以及参数选择是在原始数据确定的情况下,改变平面等值线形态来选择最主要的环节.我们选择了同一块数据用不同的软件不同的方法对该数据进行了网格化及等值线追踪(图1).对已是网格数据进行的网格加密及光滑处理也会对形成的平面等值线形态产生影响.等值线追踪 一般只需选择所需要的等值线线值,在一些软件中如MAPGIS 等也会让你选择一些诸如是否光滑、如何标记等参数,其实质是先做完等值线追踪紧接着做后续工作.从算法的角度而言这是最能表现编程人员技巧的环节,但这个环节基本不能改变平面等值线的形态.光滑线 光滑线的效果是显而易见的,但其效果主要表现在较细微的部分.2 相关软件生成图件的对比我们选择GeoExpl 、MAPGIS 、SUFER 、GEOIPAS 4种软件进行了测试.GeoExpl 软件由中国地质调查局发展中心研制.软件提供了数据预处理功能,提供了滑动平均、最近点、平均值、距离倒数加权、指数加权、累加、克里格插几种网格化的方法.提供了等值线追踪及相关地球化学标注等一整套生成地球化学等值线图的功能.但其仅提供了与MAPGIS 的数据交换接口,没有328新疆地质图1不同方法网格化对比图Fig.1 Comparative diagram in different gridding矢量数据的编辑功能.提供了大量地球化学相关的计算工具与地球化学其它图件的生成功能.该软件在全国的使用较为广泛,软件性价比较高.MAPGIS软件由中地软件公司研制生产,本次用于测试的版本为 6.5.该软件提供了数据预处理功能,距离幂函数反比加权、泛克里格、稠密数据中值选取、稠密数据高斯距离权等几种网格化的方法,提供了等值线追踪、相关地球化学标注等一整套生成地球化学等值线图的功能.提供了离散点三角网等值线追踪功能;提供了与主要地理信息软件的数据文件交换接口,提供了强大的矢量数据交互编辑手段.该软件作为地学矢量图件编辑软件广泛使用,拥有的用户众多.该软件不提供生成等值线图以外的其它相关地球化学图件的生成功能.SUFER软件为GOLDENSOFTWARE公司研制,本次用于测试的版本为8.0.该软件提供了数据预处理功能;提供了滑动平均、最近点、最小曲率、改进谢别德、自然点临近、多项式回归、径向基函数、带线性插值三角剖分、移动平均值、数据度量、局部多项式、克里格等多种网格化方法;提供了等值线追踪的方法,提供了样条平滑的网格数据圆滑处理方法;提供了与主要地理信息软件的数据文件交换接口.该软件的标注、图框等内容与我国现行地球化学等值线图的要求差别较大.GEOIPAS软件为乌鲁木齐金维图文信息科技有限公司研制,本次用于测试的版本为 1.0.该软件提供了数据预处理,提供了滑动平均、按方位取点加权网格化、趋势面拟合等几种网格化的方法,提供了等值线追踪、相关地球化学标注等一整套生成地球化学等值线图的功能;提供了大量的地球化学相关计算工具与地球化学其它图件的生成功能.与其他通用软件的数据交换能力较弱.3 结论与对策3.1制定制作地球化学等值线图的规范地球化学等值线图是地球化学采样与分析后的一种表示图,由于可用于这种表示图的干预环节和干预手段较多,使人们对其认识及评价产生了很多的差异.因此需要对其评价和编制制定相应的规范,这样更利于对一个地区地球化学特性的把握.我们认为可以从以下几方面进行规范:①数据滑动平均窗口参数的选择和可以使用的范围;②采样点分析值可移位和被分解的程度;③推荐使用的数据网格化方法;④与地球化学等值线图同时提交的处理流程和处理参数的要求.3.2研制新的地球化学等值线算法研制非统一规则网格化算法对于同一块数据的不同部分使用不同的网格化算法.这种算法的基础是能判断已有算法产生网格局部缺陷的方法.研制对已产生的矢量线的自动修改算法分析地球化学工作者对自动生成地球化学等值线图的不满意主要为:对一些线条曲率过大的部分不能容忍.目前的算法中对一些局部的不光滑有较好的处理方法,本次主要针对范围较大、曲率较大线条提出了一些处理算法[3,4].模拟人工对生成地球化学等值线进行修改的方法,归纳其思想转化为算法思想.图2是一块数据直接生成曲线、人工修改曲线和计算机模拟人工修改方法生成的结果比较.可以看出这种算法能够在很多部位接近人工修改结果.该算法主要通过以下3个步骤实现目的:①以公式Q n=θ/(p n p n-1+p n p n+1)计算一条曲线上每一点上的曲图2对比图Fig.2 Comparative diagram1——数据直接生成曲线;2——人工修改曲线;3——计算机模拟修改曲线330新疆地质率,其中p n-1、p n、p n+1是曲线上任意连续的3个点,θ为以p n为顶点p n p n-1和p n p n+1为边组成的夹角,找出Q n>A的点,其中A是由经验确定的阀值;②以（Q n1+Q n2+…+ Q nm）/L<B确定处理曲线的范围,其中L是该处理弧段的长度,B是根据经验确定的阀值;③确定以C=θn1+θn2+…θnm的弧度和P n1P nm为弦长的圆,并以这段圆弧替代这段曲线.从初步测试结果看,该算法可以较大程度地减少人工交互编辑的工作量.参考文献[1] 成建梅,陈崇希,孙红林.三角网格等值线自动生成方法及程序实现[J].水利学报,1998,(10):23-26.[2] 刘幸,直接法自动等值线[J],武汉水利电力大学学报, 1995,28(6):588-593..[3] 黄淼云,张学忠,辛冬根,等.计算机图形技术[M].北京:人民交通出版社,1995.[4] 陆润民,杨惠英.计算机绘图[M].北京:清华大学出版社,1988.SDADY ABOUT THE MMETHOD OF MAKING GEOCHEMISTRY CONTOUR DIAGRAM AND COMPARISON ON THE MAINRELEV ANT SOFTWAREHAN Xiao-ming,ZHANG Xiao-mei,WANG Rui,LI Xiao-qing,DING Li-jiang,TIAN Yu-chuan (No.1 Regional Geological Survey Party,Xinjiang Bureau of Geology and Mineral Exploration and Development,Urumqi Xinjiang,830011, China)Abstract:By the analysis on the computing method of the geochemistry contour diagram, and the comparison on the main relevant software currently in use, point out the disadvantages in the geochemistry contour diagrams computed automatically now days. Aiming at the problem of unsmooth lines in the geochemical contour diagrams, conduct the research. Describe the calculation for auto-smooth of the contour. Also covers other problems to be solved and possible solutions thereof.Key words:Geochemistry contour; gridding; Triangulation。