1++设计地形表面积测算的原理和方法研讨.doc

1 设计地形表面积测算的原理和方法 根据附件data5B.txt中所给的数据运用附录1中的MATLAB程序做出山体实体图和等高线图： 图2 实体图 图1 等高线图 根据山的实体图雏形为基础，由于附件data5B.txt中已给出了该区域的步长为5米的网格节点对应的海拔高度值，则数据已经网格化，因此引入数字地面模型中的TIN模型： 我们应用最大Z容差法来选择三角形： 最大Z容差法是将重要点的选取作为一个优化问题来处理，就是利用格网点原始高程和包含该点的三角形估算的高程差来动态选取重要点。 图3 最大容差法的思想 具体的做法是这样的： S1、连接格网DEM边界四个角点中任意对角的两个点,形成初始三角网。 S2、对每个格网点,找到包含它的三角形，内插该点在所在三角形面上的高程,求出内插高程与该点原始高程之差的绝对值（称为误差）。 S3、如果所有的格网点的误差都在最大容许的范围内则处理结束，输出TIN网；反之则进行下一步。 S4、将具有最大误差的格网点插入已存在的TIN中构成新的TIN，并返回S2。 这样我们就能构造出所有的三角形。若一个三角形的面积为则地表总面积为： ， 由于原数据为的坐标，现在改为的坐标，则总的数据点有个数据点。能反映出来的有效数据点总共有个,假设每一个有效点的坐标为 图4 有效点图 根据海伦公式： 公式里的p为半周长： ，