matlab绘制oc曲线,三角网方法绘制等值线图及其在重力资料成图中的应用

一目U舀物化探资抖成图过去常用的是网格化数据绘等值线图,然而实际工作中观测点一般是随机的,:主就需光转换成网格化数据,然后再绘图,这样做有很多不足取的地方,例如极值偏移、影响精度等等。最近随机点数据绘等值线图已有较多方法,其中笔者认为三角网方法颇具优点:,玛二角网方法绘制等值线,直接利用原始观测点插补等值点,这样可以提高等值线的精度;对待汪值部位任意小的等值线都能绘制,不会把很小面积的闭合等值线丢弃,极值也不会偏移;可以具有自由边界,即在有观测点的地区出现等值线,无观测点的地区不出现等值线,配_}_图框时切边简单,容易处理,易于接图。此外,三角网方法绘等值线图方法直观易懂,程序没计简单,占机器内存也较少,做特殊处理后运行速度也块,程序也便于移植。为了进汀1/50万,1/100万区域重力资料成图,我们编制了以三角网方法绘等值线图为基础的重力资抖成图系统,在M一1601机上采用FORTRAN语言。二三角网方法绘等值线图方法简介三角网方法绘等值线的步骤是:自动联结三角网;在三角形边上内插等值点,寻找等值线起始点和追踪等值点;找出等值线上的合适部位注记等值线的值;联结等值点绘制光滑曲线。自动联结三角网设有n个观测点(x;,y;,z;),(i=l,2…n),我们把z;认为是x一y平面上点坐标收稿日期:1日86年5月‘“12‘物化探计算技术,_客二(x;,yi)的函数,在x一y平面上,利用(x;,yi)建立起在可能条件下的最佳三角网,‘这些点(二,yi)(i“1,2.二n)就是三角形网的顶点,这个条件就是在使用周围邻近的观测点组成三角网时,应尽可能地确保每个三角形都是锐角三角形或者避免出现过大的钝角和过小的锐角。联结三角网的方法的过程是:先确定第一个三角形,然后对其三边进行扩展,又形成新的三角形,对各三角形的各边进行扩展,’这样不断下去,直到按要求联结完毕。任选一点作为第一个三角形的一个顶点,再选离该点最近的一点作为第二个顶点,第三点的选择是利用余弦定理:c“二a“+bZ一Zabcos(c),找出使C角为最大的点作为第三点(这是符合最佳三角形条件的,而且也为以后扩展三角形网时避免交叉提供了有效保证,因为这样确定的第三点,标志着它离第一、二点的距离之和为最短),这就形成了第一个三角形。从某个t三角形(例如第K个三角形)向外扩展是这样进行的:先考虑第一条边p;p:,显然,位于其对顶点p3同侧的点p.厂号应被排除,为此,建立直线方程:F(x,y)=丁ty一Ax一B;X一Xl,x:价xZ时x:=x:时其中:A二(y:一yl)/(x:一x,)图1扩边、找第三点形成B=(y:x:一yZx,)/(x:一x:)新三角形示意图(x,,v:),(xZ,了:)分别是p;点,pZ点的坐标。将各点代入F(x,y)中,可能出现的情况是F(x,y))0,F(x,y)=o,F(x,y)