DEM分析与可视化-第6章.ppt

通视分析 1）点对点通视 基于格网DEM的通视问题，为了简化问题，可以将格网点作为计算单位。这样点对点的通视问题简化为离散空间直线与某一地形剖面线的相交问题。 已知视点V的坐标为（x0,y0,z0），以及P点的坐标（x1,y1,z1）。DEM为二维数组Z[M][N]，则V为（m0,n0,Z[m0,n0]），P为（m1,n1,Z[m1,n1]）。 1》使用Bresenham直线算法，生成V到P的投影直线点集{x , y}，K=||{x , y}||, 并得到直线点集{x , y}对应的高程数据{Z[k], ( k=1,...K-1 )}，这样形成V到P的DEM剖面曲线。 2》以V到P的投影直线为X轴，V的投影点为原点，求出视线在X-Z坐标系的直线方程： （0kK） 3》比较数组H[k]与数组Z[k]中对应元素的值，如果存在Z[k]H[k]，则V与P不可见，否则可见。 2）点对线通视 点对线的通视，实际上就是求点的视野。应该注意的是，对于视野线之外的任何一个地形表面上的点都是不可见的，但在视野线内的点有可能可见，也可能不可见。基于格网DEM点对线的通视算法如下： 1》设P点为一沿着DEM数据边缘顺时针移动的点，与计算点对点的通视相仿，求出视点到P点投影直线上点集{x, y}，并求出相应的地形剖面{x, y, Z(x, y)}。 2》计算视点至每个与Z轴的夹角： 3》求得 。对应的点就为视点视野线的一个点。 4》移动P点，重复以上过程，直至P点回到初始位置，算法结束。 3）点对区域通视 点对区域的通视算法是点对点算法的扩展。与点到线通视问题相同，P点沿数据边缘顺时针移动。逐点检查视点至P点的直线上的点是否通视。 一个改进的算法思想是，视点到P点的视线遮挡点，最有可能是地形剖面线上高程最大的点。因此，可以将剖面线上的点按高程值进行排序，按降序依次检查排序后每个点是否通视，只要有一个点不满足通视条件，其余点不再检查。 点对区域的通视实质仍是点对点的通视，只是增加了排序过程。 十一、等高线追踪 1、正方形网格DEM的等高线追踪 2、三角网DEM的等高线追踪 1）三角网内插 格网点的检索 给定一点的平面坐标P（x，y），要基于TIN内插该点的高程Z，首先要确定点P落在TIN的哪个三角形中。一般的做法是通过计算距离，得到据P点最近的点，设为Q1。 然后就要确定P所在的三角形。依次取出Q1为顶点的三角形，判断P是否位于该三角形内。可利用P是否与该三角形每一顶点均在该顶点所对边的同侧（点的坐标分别代人该边直线方程所得的值符号相同）加以判断。若P不在以Q1为顶点的任意一个三角形中，则取离P次最近的格网点。 重复上述处理，直至取出P所在的三角形，即检索到用于内插P点高程的三个格网点。 高程内插 若P（x，y）所在的三角形为ΔQ1Q2Q3，三顶点坐标为（x1，y1，z1），（x2，y2，z2）与（x3，y3，z3），则由Q1，Q2与Q3确定的平面方程为 或 令 则P点高程为 2）等高线追踪 基于三角形搜索的等高线绘制算法如下 对于记录了三角形表的TIN，按记录的三角形顺序搜索。其基本过程如下： 1）对给定的等高线高程h，与所有网点高程zi（i=1,2,…，n），进行比较，若zi=h，则将zi加上（或减）一个微小正数ε 0（如ε=10-4），以使程序设计简单而又不影响等高线的精度。 2）设立三角形标志数组，其初始值为零，每一元素与一个三角形对应，凡处理过的三角形将标志置为1，以后不再处理，直至等高线高程改变。 3）按顺序判断每一个三角形的三边中的两条边是否有等高线穿过。若三角形一边的两端点为P1（x1，y1，z1），P2（x2，y2，z2）则 (z1-h)(z2-h)0表明该边有等高线点； (z1-h)(z2-h)0表明该边无等高线点。 直至搜索到等高线与网边的第一个交点，称该点为搜索起点，也是当前三角形的等高线进入边、线性内插该点的平面坐标（x，y）： 4）搜索该等高线在该三角形的离去边，也就是相邻三角形的进人边，并内插其平面坐标。搜索与内插方法与上面的搜索起点相同，不同的只是仅对该三角形的另两边作处理。 5）进入相邻三角形，重复第（4）步，直至离去边没有相邻三角形（此时等高线为开曲线）或相邻三角形即搜索起点所在的三角形（此时等高线为闭曲线）时为止。 6）对于开曲线，将已搜索到的等高线点顺序倒过来，并回到搜索起点向另一方向搜索，直至到达边界（即离去边没有相邻三角形）。 7）当一条等高线全部跟