步骤一:把数据文件转换成grd文件
1. 打开菜单Grid | Data... ,在open对话框中选择数据文件test.dat
2. 这会打开“Grid Data”对话框。在“Data Columns”中选择要进行GRID的网格数据(X和Y坐标)以及格点上的值(Z列),这里我们不用选择,因只有3列数据且它们的排列顺利已经是XYZ了,如果是多列数据,则可以在下拉菜单中选择所需要的列数据。选择好坐标XY和Z值后,在“Griding Method"中选择一种插值方法(如果你需要比原始数据的网格X和Y更密的Z数据,或是你的网格是非均匀的),则在Grid的过程中,Surfer会自动进行插值计算,生成更密网格的数据。如果你只是想绘制原始数据的图,不想插值,则最好选择距离平方反比法(inverse distance to a power)或Kriging方法。因为这两种方法在插值点与取样点重合时,插值点的值就是样本点的值,而其它方法不能保证如此。在"Output GridFile“中输入输出文件名test.grd,然后在“Grid Line Geometry”中设置网格点数。这里需要注意的是,当X和Y的数值相差很大时,这里显示的最大最小值可能有错误(即与原始数据不同),这是Surfer软件本身的问题,遇到这种情况,必须手动改正这种错误,即输入正确的最大最小值。由于我们的数据没有此类问题,因此不必手动改正最大最小值。但必须手动改正X和Y的间距“spacing” 或数据个数“#f”,这二者是相关的,改动一个,另一个自动改正。如果你的原始数据是等间距的,这里的X和Y的间距“spacing"或”#f“ 最好与原数据一致,这样可以减少插值带来的误差。我们的数据是不均匀的,所以必须插值,这里可以不进行任何改动。好了,点“ok”,画图所需要的grd文件test.grd就生成了。
步骤二:画等高线打开菜单“Map|Countour Map|New Countour Map",在“Open Grid”对话框中选择刚才输出的grd文件,如test.grd,点“ok”,则一副等高线图就画完了。
步骤三:等高线图的设置在所画的等高线图中双击鼠标,或点击右键,选中“属性”,就会出现设置等高线的各种属性,在图的边缘双击鼠标则可以改动坐标轴的属性。
四、bln文件及blank作图大家看了第三部分的图有可能会奇怪,为什么你的等高线只在海洋上有,陆地上怎么没有呢?这就是此节要讲的内容。首先讲一下bln文件的
格式
pdf格式笔记格式下载页码格式下载公文格式下载简报格式下载
,bln文件具有以下格式: length,flag x1,y1 x2,y2 ... xn,yn length,flag x1,y1 x2,y2 ... xn,yn 其中,length是背景地图曲线中点坐标的个数,如果xn=x1,yn=y1,那么所有的点将围成一个圈。如果不是这样,则是一条不闭合的曲线。 flag的值为1或0,当flag=1时,闭合曲线里面的等高线可以用此bln文件blank掉(即将闭合曲线中的等高线消除),如果flag=0,则是曲线外面的等高线被blank。我所提供的bln文件为flag=1的格式,当然可以将其改成0。简单介绍了bln文件的格式,再来介绍怎样blank等高线。
近些日子,有人问surfer8中的插值哪种好? 我个人的观点,如果你的数据的网格是均匀的,那么,最好用我提供的这种ascii 码的grd文件直接作图。这样做的好处是显而易见的,即不会产生作图误差,且不会出现其它小问题,如不能
正确处理x和y坐标(当二者差异较大时等)。我不主张在作图过程中进行插值计算。如果你的网格是非均匀的,则可能需要进行插值,我在前面介绍过在哪里选择,我个人认为比较好的是前两种( inverse distance to a power 和Kriging 方法),至于插值方法的详细介绍请参考help|contents|index 中查找interpolation主题词。因为时间关系,不再介绍插值问题。
浙公网安备 33021202002483