文档介绍：重磁勘探数据处理与解释实验一目录: 功能说明书一. 无限长薄板状体重力异常程序: clear; x=-15:1:15; y=- 10 :1: 10; [X,Y]=meshgrid(x,y); rd=; b=2; h=3; B=pi/2; G=; ra=((X+b).^2+h^2).^(1/2); rc=((X-b).^2+h^2).^(1/2); Qa=atan((X+b)./h); Qc=atan((X-b)./h); Vxz=20*G*rd*sin(B)*(sin(B)*log(rc./ra)+cos(B)*(Qa-Qc)) Vzz=20*G*rd*sin(B)*(sin(B)*(Qa-Qc)-cos(B)*log(rc./ra)) 图形等高线图图像图阴影地貌图一维向量图带线构架图图形等高线图图像图阴影地貌图一维向量图带线构架图功能说明书(一) 简介 Golden Software Surfer ( 以下简称 Surfer) 是一款画三维图(等高线,image map, 3d surface )的软件,该软件简单易学,可以在几分钟内学会主要内容,且其自带的英文帮助( help 菜单) 对如何使用 surfer 解释的很详细, 其中的 tutorial 教程更是清晰的介绍了 surfer 的简单应用,应该说 surfer 软件自带的帮助文件是相当完美且容易阅读的, 只要学过英语的人都可以很快上手。 Surfer 是具有插值功能的绘图软件, 因此, 即使你的数据是不等间距的, 依然可以用它作图。但依据作者的经验,最好不使用 Surfer 自带的插值功能,尤其是要精确确定等高线时。由于 surfer 是美国的一个软件,它不提供对中文的支持,这可以算的上一个小的遗憾。 Surfer 的主要功能是绘制等高线图( contour map ) ,此外它还可以绘制 post map , classed post map , vector map , image map , wireframe map , 3d surface map , 等形式的图形。其功能是比较强的, 但没有各种投影变化是它的一大缺点。尤其是在等高线领域,这不能不说是它的应用受到限制的地方。由于 surfer 软件没有中文说明书,对一些初学者来说可能会存在上手较难的问题,鉴于此种需求,编写了这一初学者参考手册,希望对大家有所帮助。(二) 等高线的绘制 Surfer 的最主要的功能是绘制等高线图, 但并不是我们具有了数据文件就可以直接绘制等高线, surfer 要求绘制等高线的数据有特殊的格式要求,即首先要将数据文件转换成 Surfe r 认识的 gr d 文件格式, 才能绘制等高线( 当然, 可以直接生成 surfe r 接受的 ascii 码的 grd 文件格式,这样就可以直接作图,此方法将在后面介绍,首先我们介绍常用的作图方法)。假设你有三列数据分别为 X ,Y ,Z , 其中 Z 为点(x ,y ) 处的值, 存在文件 t 中,其中第一列是 X 坐标,第二列是 Y 坐标,第三列是( x ,y )上的值 Z ,则绘制等高线的步骤如下: 步骤一:把数据文件转换成 grd 文件 1. 打开菜单 Grid | Data... , 在 open 对话框中选择数据文件 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 “中输入输出文件名 ,然后