《Surfer使用说明书》由会员分享，可在线阅读，更多相关《Surfer使用说明书（29页珍藏版）》请在金锄头文库上搜索。
1、11.网格文件及其建立一、简介网格文件规则矩形网格的三维数组 X，Y，Z 文件。用来画等值线图、网框图、矢量图、粘贴图、立体阴影图、色谱图，计算面积、体积和画剖面图。.GRD 文件格式：.GRD 文件：GS 二进制文件( 缺省网格文件格式 )、GS ASCII 文件.DAT 文件：ASCII XYZ 数据文件 .Surfer7 还可以直接使用 USGS .DEM, GTopo30 .HDR和 SDTS .DDF 文件。 二、由离散点的 XYZ 文件创建 .GRD 文件为了产生等值线图或线网图 , 或执行需要网格文件的操作如网格代数处理、计算体积和面积, 光滑或计算残差，必须预先对数据文件进行网
2、格化处理。由 GridData 打开所需的数据文件后，出现散点数据插值法对话框 Scattered Data Interpolationd。该对话框由 4 个 Tab 组成。1 Date数据列组（Data Columns）：设定 X、Y 和 Z 的值对应于数据文件内的哪一列。相同数据点处理方式选择。由下拉菜单 Duplicates 可以选择对有相同的 X、Y 坐标但有不同的 Z 值的两个或多个点 Z 值处理办法，包括取平均值（ Average）、保留一个并删除其它值(Delete)、对各点 Z 值加和(Sum) 等。有些网格化方法还可以选择忽略这类数据点(Ignore)。数据筛选（Filter
3、）：可以使用布尔表达式对原始数据点进筛选。原始 XYZ 数据对话框提供数据统计。显示在数据文件内指定列中，相应的 X, Y和 Z 值的数据点数、最小值、最大值、平均值、方差和标准差等。2 Genaral网格线参数（Grid Line Geometry）：选择 X,Y 方向网格的取值范围；设定产生网格文件时网格间隔和网格线的密度。改变网格线的密度，网格的间隔宽度会自动随着改变。网格线的密度高则较精密，但文件长度大，构筑网格费时间，网格化方法（Gridding Method）：设定采用网格化的方法。有下拉式菜单提供多种可用的网格化方法清单。9 种网格化的数学方法见本章第二节。用不同的网格化2方法处
4、理数据时会产生不同的结果。选择时，首先应使你的数据最接近实际情况。其次应使表现数据的图形为最好。Option 每种网格化方法都有相应的参数设定。输出网格文件（Output Grid File）：指定生成的网格文件的路径和文件名。缺省名与数据文件同名，扩展名.GRD，可点击浏览按钮以改变其路径或文件名并选择网格文件类型。3 Search搜寻功能决定在内插网格结点时，那些数据点参加插值计算。对于要用到全部数据文件的网格化方法，例如多项式回归法、最小曲率法和线性插值的三角形法等，搜寻功能不起作用，Search 键不能激活。有以下 4 种搜寻类型：全部数据（All Data）：每一个结点 Z 值的计算
5、都要用到全部数据点。对于较少的数据组，如少于 250 组数据时，最常用到全部数据搜寻。简单搜寻(Sample)：只用最靠近某结点的一些数据来计算该结点的 Z 值。数据点的数目由窗口右面 Search Rules 框内的 Data Per Sector 来设定。只有在搜寻椭圆（Search ellipse）以内的数据点才被用来计算。象限搜寻（Quadrant）：将一个结点周围分为 4 个 90的象限，软件自动分析每个象限的数据点数目，搜寻 4 个扇区内最近的数据点，而忽略掉该象限中这些最近点以外的数据点。卦限搜寻（Octant）：将一个结点周围分为 8 个 45的卦限，软件自动分析每个卦限的数据
6、点数目，搜寻 8 个扇区内最近的数据点。较象限搜寻慢，但构成比较平滑的网格。当观察点明显成簇状不均匀分布时，推荐使用象限或卦限搜寻方法，以消除或减少简单搜寻所造成的失真。搜寻椭圆（Search Ellipse）：在搜寻方式为简单搜寻、象限搜寻或卦限搜寻情况下，计算每个结点 Z 值时搜寻数据点的范围。搜寻椭圆以外的数据点，进行插值时不予考虑。缺省为圆形，表示在所有方向搜寻距离相同。搜寻椭圆的参数有：半径：半径 1 和半径 2 分别是用数据单位表示的 X 和 Y 方向的搜寻距离。角度：+X 轴和半径 1 之间的倾斜角度（-360+360）。结点 Z 值为 0 说明在特定位置，根据用户设定的搜寻规律
7、，没有足够的数据点来3计算结点 Z 值。搜寻规律（Search Rules）：在搜寻方式为简单搜寻、象限搜寻或卦限搜寻情况下，设置用来计算每个结点 Z 值的数据点数。与搜寻椭圆相关。每扇区点数：设置每扇区的数据点数，对于简单搜寻、象限搜寻和卦限搜寻，扇区数分别是 1、4 、8 ，缺省每扇区点数分别为 24、8、3 。最少总点数：对一个网格结点进行插值时，最少所需数据点数。最多空扇区：对一个网格结点进行插值时所容许的最多空扇区数。4 Faults and BreaklinesFaults 用.BLN 文件格式保存的 X、Y 二维多段线或多边形，用来限定网格化的范围。如对多段线一侧的原始数据点进行
8、网格化，忽略另一侧。或对一多边形之内（或之外）的原始数据点进行网格化。File Containing Fault Traces 框给出 Faults 的.BLN 文件所在的路径。有一些网格化模型可以使用 Faults。Breaklines 用.BLN 文件格式保存的三维数据集，X 、Y、Z 给出每个点的坐标。不同于 Faults，Breaklines 不阻断信息流。当 Grid 算法遇到 Breaklines 时，用线形内插的方法联合确定结点的 Z 值。File Containing Breaklines 框给出Breaklines 的 .BLN 文件所在的路径。Breaklines 用于定义
9、流线、山脊和其他斜率的突变。三、由函数关系生成.GRD 文件GridFunction 打开网格函数对话框。在函数编辑框中输入数学函数 Z= f(X,Y)。数学函数逐点的计算每一个 Z 值，并将其写入到输出的网格文件中。其中网格结点的数目取决于所选定的 X 和 Y 的最小值、 最大值和增量。所采用的数学函数中能使用系统的任何内部函数。例如, 按给定的函数关系，X、Y 的取值范围和增量产生一个网格文件。缺省的输出网格文件名为 OUT.GRD，可改变文件的路径和名字。四、网格代数处理由具有相同 X,Y 网格尺度的两个网格文件按照一定的函数关系 C= f(A,B) 生成一个新的输出网格文件。这里 A
10、和 B 是输入的网格文件，C 是输出网格文件。从输入网格文件的相应结点上(具有相同 X,Y 坐标的网格结点) 计算函数值，计算的结果送到具有相同坐标的输出网格文件内。网格数学也能用于单一的网格文件。例如 ,使用函数 C= log10(A)，则网格文件4C 的每一个结点的值将为网格 A 相应的每一网格结点值的以 10 为底的对数。也可以进行简单的数学操作,例如 C=A-100。这时输出网格文件 C 比输入网格文件 A 低100 单位。由 GridMath 打开网格代数对话框。输入网格文件 A 组指定执行网格数学操作的第一个网格文件。输入网格文件 B 组指定用于网格数学操作的第二个网格文件。如果只
11、想在一个网格文件上执行操作, 可使其为空白。第二个网格文件应包含与第一个网格文件相同的网格行、列数, 和相同的 X 和 Y 列。缺省的输出网格文件名为MATHT.GRD，可改变文件的路径和名字。底部的编辑框用于输入函数 C= f(A,B)。 五、网格算法（Grid Calculus）提供解析网格文件的工具，用来解释一些在等值线图和线网图中不明显的网格文件的特点。由 GridCalculus 打开网格算法对话框。选择网格算法进行处理，生成一个新的网格文件。由网格算法对话框可以选择以下 4 种算法之一：方向导数、地体模型、微分和积分算符以及 Fourier 和谱分析1）方向导数（Direction
12、al Derivative）方向导数给出沿着规定的方向网格表面的斜率或斜率变化率的信息。斜率是有方向性的，如果规定方向为东，但倾斜方向为北，则虽然向北有斜率，在那个点上方向导数仍然为零。操作时可分为一阶导数、二阶导数和方向曲率三种方式。一阶导数（First Derivative）: 一阶导数给出表面沿方向线的斜率. 由一阶导数的网格文件生成的等值线图绘制出沿着指定的方向具有相同斜率的等值线。在各网格结5点处, 如向上倾斜则斜率为正，如向下倾斜则斜率为负。若斜率接近正或负的无限大，表示在向上或向下的向上接近垂直。二阶导数（Second Derivative）：二阶导数给出沿着方向线斜率的变化率。
13、由二阶导数的网格文件生成的等值线图绘制出穿过表面斜率变化率相同的等值线。其值随表面上升下降也有正值和负值；数值的大小反映表面缓陡变化的快慢。方向曲率（Directional Curvature）：方向曲率是曲面沿着规定方向线，在切面倾角的变化率。方向曲率采用变化率的绝对值，因此它是一个正数。方向曲率类似二阶导数。2） 地形模型(Terrain Modeling) 在表面上每一网格结点处 ,Surfer 确定了该点最陡的斜坡和最陡斜率的方向(梯度方向) 。地形模型的值是从 90 到-90 度. 最陡斜率零度是水平线,90 度是垂直上升的，而 -90 是垂直下降的。5 种地形模型的操作如下；地形斜
14、率、地形方位、纵剖面曲率、平面曲率和切线曲率。可以选择地形模型的临界值，缺省为一极小值。地形斜率（Terrain Slope）：穿过表面，斜率梯度的方向可以改变。由地形斜率的网格文件生成的等值线图或线网图绘制出等最大斜率线。该功能与一阶方向导数相同，用来确定表面任何一点的斜率，在考虑自动确定图中每一个点的的梯度方向时更为有用。梯度方向可由地形方位功能得到。地形方位（Terrain Aspect）：给出每一网格结点最大斜率的方向。比如表面的水流方向严格垂直于等值线的方向。其值用方位角来表示，0 度为北；90 度为东等。可以利用地形方位功能创建一幅反映流动方向的粘贴图。纵剖面曲率（Profile
15、Curvature）：纵剖面曲率确定每一个结点延梯度方向斜率变化率的增减情况。由纵剖面曲率的网格文件生成的等值线图显示穿过表面，最大斜率的等变化率线。该功能可与二阶方向导数对比，但更加有用，因为它可以自动地判定表面每一点的下降方向，然后确定该点沿此方向斜率的变化率。平面曲率（Plan Curvature）：平面曲率反映了在水平面上地形方位角测量的变化率，描述了表面水流的会聚或发散。由平面曲率的网格文件生成的等值线图显示等会聚或等发散线。在网格文件中的负值表示发散，而正值表示会聚。切线曲率（Tangential Curvature）：切线曲率与平面曲率类似，也用来反映会聚或发散，但它测量的是在垂
16、直面上的情况。该功能描述表面水流的会聚或发散，可以想象为沿垂直于等值线方向的垂直面的表面曲率。3）微分和积分算符(Differential Integral Operator)：包括了 4 种功能：梯度算符、拉普拉斯算符、双谐算符和积分体积6梯度算符（Gradient Operator）：梯度算符在表面的任何点上产生最大斜率的网格，等于地形斜率操作, 但是梯度算符为一数值，而不是角度。水平表面梯度算符为零, 而接近于垂直面时梯度算符接近无限大。拉普拉斯算符（Laplacian）：拉普拉斯算符用于提供表面充电和放电的测量。地下水、热量和电荷是三个守恒物理量的例子，其局部流动速率与梯度成正比。拉普拉斯算符对这些物理量的局部体积进行量化，流入时或在充电区 0，而流出时或在放电区1000 数据点) ，最小曲率法最快，网格足以代表原始数据特点。线性内插三角形法网格化较慢，网格有足够的代表性。以下的建议仅仅是一般的推荐。距离反比法最快，但是围绕数据点，有产生 “牛眼”效应的趋
