Surfer软件应用.ppt

Surfer软件应用Surfer软件简介 各种等值线图和三维图形是地学、工程、矿业各种等值线图和三维图形是地学、工程、矿业等领域的研究人员最常用的图件，等领域的研究人员最常用的图件， SurferSurfer是美国是美国GoldenGolden软件公司的产品，主要用于绘制等值线图及软件公司的产品，主要用于绘制等值线图及相应的三维图形，因此成为上述各行业科技工作者相应的三维图形，因此成为上述各行业科技工作者的常备软件的常备软件GoldenGolden软件公司是由软件公司是由Dan SmithDan Smith和和Patrick Madison 1983Patrick Madison 1983年在美国创立的该公司一年在美国创立的该公司一直致立于图形软件的开发，于直致立于图形软件的开发，于19851985年推出年推出SurferSurfer，，后又于后又于19861986年推出年推出GrapherGrapher，，现在，现在，GoldenGolden软件公软件公司在科学图形软件领域已处于领先地位，用户包括司在科学图形软件领域已处于领先地位，用户包括矿业、工程、医药、地学、生物等领域的研究人员。

矿业、工程、医药、地学、生物等领域的研究人员 目前目前SurferSurfer的最高版本为的最高版本为8 8，， SurferSurfer是在是在WindowsWindows和和NTNT操作系统下最强大的、最灵活的和较操作系统下最强大的、最灵活的和较容易使用的绘制等值线图及三维立体图软件包容易使用的绘制等值线图及三维立体图软件包Surfer 8的主要功能 ØØ支持支持1212种内插方法对离散的种内插方法对离散的XYZXYZ数据格网化，生成数据格网化，生成规则的格网数据；规则的格网数据； ØØ强大的地学数据分析功能，支持多达强大的地学数据分析功能，支持多达1212种变异函种变异函数理论模型，可计算残差，进行地形分析和计算数理论模型，可计算残差，进行地形分析和计算体积、面积等，体积、面积等， ØØ绘制等值线图绘制等值线图(Contour Maps)(Contour Maps)；； ØØ绘制曲面图绘制曲面图(Surface(Surface）；）； ØØ从从Surfer8Surfer8开始有几种常用的开始有几种常用的GISGIS软件数据格式可软件数据格式可直接输人直接输人Surfer8Surfer8系统系统ØØSurfer8Surfer8具有简单的数学运算、数据统计、平滑、具有简单的数学运算、数据统计、平滑、滤波、微分、积分、傅立叶和谱分析等多种数据滤波、微分、积分、傅立叶和谱分析等多种数据分析功能，分析功能，ØØSurfer8Surfer8具有趋势面分析、体积、面积计算、地形具有趋势面分析、体积、面积计算、地形分析、剖面计算等三维空间分析功能。

分析、剖面计算等三维空间分析功能讲述主要内容pp数据文件的编辑数据文件的编辑ØØ转换转换( (TranformTranform) ) ØØ统计统计(Statistics)(Statistics) 常用图件绘制常用图件绘制ØØ绘制等值线图绘制等值线图ØØ绘制粘贴图（原始数据底图）绘制粘贴图（原始数据底图）ØØ绘制地图空白区域绘制地图空白区域ØØ绘制线框图绘制线框图ØØ绘制曲面图绘制曲面图ØØ绘制堆叠图绘制堆叠图ØØ绘制叠置图绘制叠置图pp三维空间分析ØØ趋势面分析体积、面积计算文件的输入与输出数据文件的基本操作ØØ数据文件至少要有三列数据值，前两列输入的数据为X、Y坐标，第三列以后的数据为Z值ØØ第一行的文字描述可帮助用户确定所在列数据的类型和含义ØØ可读入XLS、TXT、DAT等数据格式的文件转换(Tranform) ØØ计算从0到更多数据项的一列数据，计算结果取决于所定义的等式，这个等式能含有对应于数据表列的字母变量一般用来对某一列数据进行数据变换，如对数指数变换等，以修复某些等值线图形绘制中的X、Y轴比例差异过大的问题ØØ如果任何单元格含有文本或没有值，公式将不计算该数据行。

例如:使用公式C=A十B将计算列A的值和列B的值之和，并将结果放在列C统计(Statistics) ØØ在在工工作作表表中中选选定定要要统统计计的的数数据据，，选选所所需需的的统统计计特特征数ØØSurfer8Surfer8中中统统计计功功能能有有：："Number "Number of of Values(Values(数数据据点点 数数 ) )、、 "Sum("Sum(求求 和和 )")"、、 "Minimum("Minimum(最最 小小 值值 ) )、、"Maximum("Maximum(最最大大值值)")"、、"Range("Range(值值域域)")"、、"Mean("Mean(平平均均 值值 )")"、、 "Median("Median(中中 位位 数数 )")"、、 "First "First Quartile/25 Quartile/25 Percentile(Percentile(第第 2525百百 分分 位位 数数 )")"、、"Third "Third Quartile/75th Quartile/75th Percentile(Percentile(第第7575百百分分位位数数)")"、、"Standard "Standard Error Error of of the the Mean(Mean(平平均均值值的的标标 准准 误误 差差 )")"、、 "Variance("Variance(方方 差差 )")"、、 "Standard "Standard Deviation(Deviation(标标 准准 差差 )")"、、 "Skew("Skew(偏偏 斜斜 度度 )")"和和Kurtosis(Kurtosis(峰值峰值)")"等等。

等等等值线图绘制(Contour Map) 等值线图是指在一个表面上等值点的连线，在空间分析中常采用等值线图来表示数据的变化趋势从图形学的角度看，等值线图是指在X-Y平面上f(x,y)=z的轨迹分布这里的z为某一常数值，该值表示的物理意义可以是地形高程数据、化探数据场中的某一元素含量值、温度场中的温度数据和压数据等等在数据变化较明显的区域中，等值线密度较高，反之，在数据变化，等值线密度较小这样，数据的空间分布趋势就可通过等值线的密度和方向来分析等值线图绘制步骤1.1.建立或输入数据文件2.2.离散数据网格化，形成GRD文件 注意数据内插方法的选择注意数据内插方法的选择3.3.用菜单“Map｜Contour Map｜New Contour Map ”生成等值线图4.4.修改等值线图的各部分属性5.5.保存文件（SRF）绘制粘贴图(Post Map) 粘贴图的功能是用标志和符号把某些指定的点或原始数据标出来例如，在等值线图上标出观测站的位置、采样点位置、石油钻井位置、城镇村庄等地理位置信息 从包括X和y坐标的数据文件来新建粘贴图，坐标被用来确定符号在粘贴图上的位置。

数据文件可包含附加的信息如符号大小、颜色、标志、位置和角度等 新建的粘贴图独立于在绘图窗口的其它的地图例如，在绘图窗口中已经存在等值线图，再新建一个粘贴图，将产生2个没有联系的地图，它们各自具有自己的坐标轴和缩放比例参数在等值线图上粘贴数据点，在绘图窗口新建没有联系的等值线图和粘贴图，选中2个地图，选择菜单命令“Map（地图）Overlay Maps(叠置图），这将在一套坐标轴和缩放比例参数下将2个地图合并成一个地图可把任何数量的粘贴图和任何其它的地图合并绘制线框图（Wireframe） 线框图是格网文件的三维表达线框图是通过为XY格网(格网行、列)画线产生块状图表，在每个行和列(即在每个格网节点)的交叉点的格网Z值与表面的高度成比例格网文件里的行和列的数决定线框图里的绘制X和y线数线框图可用格网文件或美国地质勘察数字高程模型(USGS DEM)数据生成更改线框图的方向和倾斜度更改线框图的尺度绘制曲面图（Surface） 曲面图指空间坐标的函数，用以表示变量的分布情况绘制地图空白区域（Blank） 空白文件由定义一个封闭多边形的一系列了X、Y坐标值指定，多边形的边必须闭合，一个空白文件用来定义空白区域。

空白文件的格式为，第一行是标志行，它的第一个数字表示下边的点数，第二个数是边的标志符，如果边的标志是0，多边形外边的格网单元被空白，如果边的标志是1,多边形里的格网单元被空白整个数列表示多个空白的区域例：一个四边形区域被空白的数列实例5(点数) 1 边标志(四边形内部被空白)3 2 (第1点)8 6 (第2点)3 3 (第3点)4 2 (第4点)2 2 (第5点)绘制地图空白区域步骤I.I.创建一个空白文件ØØ在绘图窗口里选择菜单命令在绘图窗口里选择菜单命令"Map("Map(地图地图) )｜｜Contour Map(Contour Map(等值线图等值线图)|New Contour Map()|New Contour Map(新新建等值线图建等值线图)")"，在弹出的，在弹出的" "打开格网打开格网" "对话框里，对话框里，新建一个等值线图新建一个等值线图。

ØØ双击选中等值线图，选择菜单命令双击选中等值线图，选择菜单命令"Map("Map(地图地图) )｜｜Digitize(Digitize(数字化数字化)")"，光标变成十字细丝，光标变成十字细丝; ;ØØ在绘图窗口绘制空白所需的位置，系统自动在在绘图窗口绘制空白所需的位置，系统自动在“ “数字化编辑窗口数字化编辑窗口” ”里输入当前位置的坐标；里输入当前位置的坐标；ØØ保存为空白文件（保存为空白文件（* *. .blnbln））或数据文件或数据文件（（* *. .datdat）） II.II.空白等值线图绘制ØØGrid｜Blank创建新的格网文件（从原始格网文件中提取数据）ØØMap｜Contour Map|New Contour Map|选择新建的格网文件空白文件边标志的修改（1和0互换） 重新建立格网文件，绘制新建立的格网文件图空白文件的边标志＝1空白文件的边标志＝0绘制堆叠图（Stack Maps） 在我们的研究工作中，有时想在一个图件中上、下堆叠已经绘制好的平面图(如等值线图、粘贴图等)和线框图，这样既美观，又具有一目了然的效果绘制平面一立体图对ØØ绘制等高线图，设置地图倾斜角度，旋转角度ØØ绘制线框图ØØ全部选择后，选择Map｜Stack Map，添加垂直线，完成平面—立体图对的绘制绘制叠置图（Overlay Maps） 为了增强图件的可视化效果，可以框图或等值线图上叠置包含地理信息方向和大小的矢量地图，也可以框图或等值线图上叠置包含如采样点位置、石油钻井位置、城镇村庄等地理位置信息的粘贴图。

线框图上叠置粘贴图ØØ绘制线框图ØØ绘制粘贴图ØØ全部选择后，选择Map｜Overlay Maps三维空间分析 Surfer8中三维空间分析是指对XYZ数据的第三维变量Z的分析，变量Z可能是地形，也可能是降雨量、污染物浓度、土壤酸碱度、某一化探元素的值等变量p 趋势面分析趋势面分析p体积、面积计算体积、面积计算  趋势面分析 某种地理属性在空间的连续变化，可以用一个平某种地理属性在空间的连续变化，可以用一个平滑的数学平面加以描述思路是先用已知采样点数滑的数学平面加以描述思路是先用已知采样点数据拟合出一个平滑的数学平面方程，再根据该方程据拟合出一个平滑的数学平面方程，再根据该方程计算无测量值的点上的数据这种只根据采样点的计算无测量值的点上的数据这种只根据采样点的属性数据与地理坐标的关系，进行多元回归分析得属性数据与地理坐标的关系，进行多元回归分析得到平滑数学平面方程的方法，称为趋势面分析它到平滑数学平面方程的方法，称为趋势面分析它的理论假设是地理坐标（的理论假设是地理坐标（x x，，y y）是独立变量，属性）是独立变量，属性值值z z也是独立变量且是正态分布的，同样回归误差也也是独立变量且是正态分布的，同样回归误差也是与位置无关的独立变量。

是与位置无关的独立变量 多项式回归分析是描述长距离渐变特征最简单的多项式回归分析是描述长距离渐变特征最简单的方法多项式回归的基本思想是用多项式表示线、方法多项式回归的基本思想是用多项式表示线、面，按最小二乘法原理对数据点进行拟合线或面面，按最小二乘法原理对数据点进行拟合线或面多项式的选择取决于数据是一维的还是二维的多项式的选择取决于数据是一维的还是二维的 通常许多地质数据一般都包含着与空间地理坐标（x,y）相关的三部分信息： ①反映区域性变化的：数据中反映总体的规律性变化的部分，由地质区域构造、区域岩相、区域背景等大区域因素所决定 ②反映局部性变化的：反映局部范围的变化特征 ③反映随机性变化的：它是由各种随机因素造成的残差 趋势面分析是对地质特征的空间分布进行研究趋势面分析是对地质特征的空间分布进行研究和分析的一种方法，它是用某种形式的函数所代表和分析的一种方法，它是用某种形式的函数所代表的曲面来逼近该地质特征的空间分布这个函数从的曲面来逼近该地质特征的空间分布这个函数从总体上反映了采样数据的区域性变化趋势，称为趋总体上反映了采样数据的区域性变化趋势，称为趋势面部分；采样数据的实测值与这个函数对应值之势面部分；采样数据的实测值与这个函数对应值之差，称为残差部分。

差，称为残差部分 它反映了局部性的变化这就是说，把采样数它反映了局部性的变化这就是说，把采样数据的实测值分解成趋势面和残差两部分趋势面部据的实测值分解成趋势面和残差两部分趋势面部分用一个函数表示，它反映采样数据的总体变化，分用一个函数表示，它反映采样数据的总体变化，可以认为是由大范围的系统性因素引起的；残差部可以认为是由大范围的系统性因素引起的；残差部分反映了局部性的变化特点，可以认为由局部因素分反映了局部性的变化特点，可以认为由局部因素和随机因素引起的和随机因素引起的 趋势面分析，一方面是要研究某些地质参数的区趋势面分析，一方面是要研究某些地质参数的区域性变化规律，另一方面是要圈出该参数变化的局域性变化规律，另一方面是要圈出该参数变化的局部异常，发现通常难于发现的低缓异常带，这样就部异常，发现通常难于发现的低缓异常带，这样就需要求出残差值，作残差分析第需要求出残差值，作残差分析第i i点的残差值是点的残差值是观测值与趋势值之差可以按一定间隔作出残差等观测值与趋势值之差可以按一定间隔作出残差等值线图 引进残差以后，原始数据即被分解成为趋势值与引进残差以后，原始数据即被分解成为趋势值与残差值两部分，趋势值反应区域性大范围内的变化残差值两部分，趋势值反应区域性大范围内的变化情况，残差值反映局部变化特点。

二者结合起来可情况，残差值反映局部变化特点二者结合起来可以帮助人们深人地作地质分析利用这一方法，可以帮助人们深人地作地质分析利用这一方法，可以发现有意义的低缓异常带（区）用趋势面拟合以发现有意义的低缓异常带（区）用趋势面拟合物探和地球化学指标时，利用残差图可找出异常带，物探和地球化学指标时，利用残差图可找出异常带，为进一步深人研究矿产资源形成的地质条件提供不为进一步深人研究矿产资源形成的地质条件提供不可缺少的信息可缺少的信息分布趋势实现过程 SurfeSurfe软件中格网菜单下的数据命令（软件中格网菜单下的数据命令（GridGrid｜｜DataData），）， ““Polynomial RegressionPolynomial Regression（多元回归法）（多元回归法）””方法可根据方法可根据空间的采样数据拟合一个数学曲面，用该数学曲面来反映空间的采样数据拟合一个数学曲面，用该数学曲面来反映空间数据分布的变化情况，因此该插值方法实际上是一个空间数据分布的变化情况，因此该插值方法实际上是一个趋势面分析程序，可利用该程序完成从原始数据中拟合出趋势面分析程序，可利用该程序完成从原始数据中拟合出数据变化的趋势。

点击数据变化的趋势点击““Advanced OptionsAdvanced Options””按钮，则按钮，则弹出弹出““Regression Advanced OptionsRegression Advanced Options””对话框；可在对话框；可在Surface DefinitionSurface Definition选择框中选取多元回归类型，提供选择框中选取多元回归类型，提供选用的回归类型有选用的回归类型有““Simple Planar SurfSimple Planar Surf（简单平面）（简单平面）””、、““BiBi－－Linear Linear SadleSadle（双线性鞍）（双线性鞍）””、、““Quadratic Quadratic SurfaceSurface（二次曲面）（二次曲面）””、、““Cubic SurfaceCubic Surface（三次曲面）（三次曲面）””和和““User Defined PolynomialUser Defined Polynomial（用户自定义多项式）（用户自定义多项式）””选择不同的多项式类型，在。

选择不同的多项式类型，在““ParametersParameters””选择框中选择框中的值会相应改变，其中的值会相应改变，其中Max Total Max Total 框为指定多项式中框为指定多项式中X X和和Y Y最大的次数之和最大的次数之和“多元回归法高级选项”对话框剩余值（残参）计算 残参值是相应位置原数据的观测值与趋势值之残参值是相应位置原数据的观测值与趋势值之差计算剩余的命令在差计算剩余的命令在““GridGrid｜｜ResidualsResiduals””，首，首先在先在““Open GridOpen Grid””对话框中选择某一元素的格网对话框中选择某一元素的格网文件（文件（* *. .grdgrd），单击确定后出现），单击确定后出现““OpenOpen””对话框，对话框，在其中选择与该元素相对应的数据文件（在其中选择与该元素相对应的数据文件（* *. .datdat），），单击打开，出现单击打开，出现““Grid ResidualsGrid Residuals””对话框（图对话框（图2 2），分别指定数据文件中），分别指定数据文件中X X、、Y Y、、Z Z数据所在的列，数据所在的列，残参相应的存储在最后一列。

单击残参相应的存储在最后一列单击““OkOk””后，系统后，系统自动打开一个选择元素的工作表，并在表中最后一自动打开一个选择元素的工作表，并在表中最后一列显示新产生的残参值列显示新产生的残参值“Grid Residuals”对话框 元素异常圈定 ØØ剩余值是某一位置上原始观测值与经过该点的趋势值之差，如剩余值为正，说明剩余值为高出趋势的部分，如剩余值为负，说明剩余值为低出趋势的部分因此，求出各点的剩余值后，即可按一定的间隔作出剩余分布等值线图ØØ值得注意的是，剩余值并不完全等于局部异常，它只是局部异常分量与随机因素造成的随机分量之和 在地质找矿中，高于趋势部分的剩余意义较大，为元素富集区域，故实际工作中常用正剩余值圈定异常图随机分量常采用正剩余的平均值来代替式中： 是正剩余的个数； 表示第i点的正剩余不同方法圈定Ba元素异常分布对比图 不同方法圈定Cu元素异常分布对比图 趋势面问题讨论 趋势面的拟合精度过高会因数学曲面过于逼近实际分布而难于反映分布的主体特征，达不到描述空间趋势的目的因此一般用处理后的数据作二次或三次趋势面计算，趋势面即可代表区域元素的背景值。

体积、面积计算 体积计算方法通常用于确定两个格网文件之间的实体的体积所谓体积通常是指空间曲面与一基准平面之间的空间的体积，在绝大多数情况下，基准平面是一水平面基准平面的高度不同，尤其是当高度上升时，空间曲面的高度可能会低于基准平面，此时出现负的体积在对地形数据的处理中，当体积为正时，工程中称之为“挖方”，体积为负时称之为“填方” 通常采用近似的方法计算体积，由于空间曲面的表示方法的差异，近似的计算方法也不一样,Surferx8同时采用三种方法计算体积：Trapezoidal Rule(梯形法）、Simpson's Rule(普森法）和Simpson’s 3/8 Rule(辛普森3/8法）估计三种体积计算方法结果的偏差，可以对三种方法的计算结果进行对比选择菜单命令"Grid｜Volume”，在格网体积报告中显示近似体积值、正体积负体积、正负体积之差、正面积、负面积、总面积等，结果能被另存为一个ASCII格式的文本文件，或者粘贴在剪贴板上。