如何根据离散点自动绘制等值线(等高线)之 三角形法

自动绘制等值线的方法从技术方向上看可以分为两大类,插值和曲线拟合.其中曲线拟合总的来说效果不如插值算法经典和应用广泛,效果也较逊色.这里着重介绍插值算法.

其中插值算法中,按照方式不同分为离散点客观化和三角网方式.两者区别在于三角网计算主要在生成三角网过程,省去了插值到格点的过程.而客观分析过程则是将离散点分析到格点后再内插到细网格,然后大多利用追踪法生成等值线,也有在这里再使用曲线拟合.

1、三角形算法

a、首先生成delaunay三角形，这一点在我的帖子"delaunay triangulation之丰衣足食“内有源程序，大家可以参考。

b、随后需要在三角形的边上插补等值点。要确定某个三角形的边上是否有等值点，需要进行判断和处理。注意：如果某原始数据点和等值线值相同，将该点改变一个微量。如果一个三角形三顶点的值相同则各边无等值点。如果一个三角形的任意边两端点（A、B〕的Z值（Za、Zb）满足满足（Zd-Za）*（Zd-Zb）<0,其中Zd代表等值线的值，则该边必有等值点，其平面位置是Xd=Xa+(Xb-Xa)*(Zd-Za)/(Zd-Za) ,  Yd=Ya+(Yb-Ya)*(Zd-Za)/(Zb-Za)。每个三角形上不可能三边都有同值的等值点，另一边上必定有同值的等值点。

c、等值点的追踪。

为了能将内插的等值点顺序追踪排列，绘出等值线，还必须找出相互重叠的环形网内所计算的等值点间的平面位置关系。因每个环形网都是由多个三角形组成的，我们先简单分析一下单个三角形中存在等值点的情况。由于不必考虑等值线穿过端点，如果一个三角形的边上存在等值点的话，只可能在某两条边上存在等值点，而不可能三条边上同时都有。也就是说，只要三角形一边上存在等值点，则其余的两条边中必有一边存在等值点