第三次国土调查县级坡度图制作关键技术研究.docxVIP

第三次国土调查县级坡度图制作关键技术研究 　　摘要：针对安徽省分县坡度分级图生产过程中dem坡度计算模型分析、分幅坡度栅格数据矢量化与图斑合并、分县数据拼接和裁切、分幅和分县坡度分级图小面融合、图斑边界平滑等关键技术进行研究。其中构面线拓扑构建和海量图斑融合到邻接非同值面的处理是最关键环节和技术难点。研究成果为类似坡度分级图生产技术处理方案提供参考，在坡耕地确定、地灾、地貌模拟、环保等有广泛应用，使用的内存过渡、分段循环、递归判断等技术处理手段对类似的海量数据处理提供借鉴。 　　关键词：第三次国土调查坡度分级图horm拟合构面线构建融合 　　中图分类号：p208文献标识碼：a文章编号：1674-098x（2019）09（c）-0038-02 　　作为全国国土调查的重要基础资料和成果，坡度分级图（简称坡度图）是查清耕地和基本农田的数量、质量和空间分布的重要指标，在林地规划利用管理、水文模型建立、土壤侵蚀分析、水土流失监测、滑坡等地质灾害预警监测、地貌形态模拟、生态环境研究等地学分析领域也有着广泛的应用。二调坡度图由于原始dem数据生产时间距今太过久远、现势性不好;其测量手段落后、植被高扣除混乱等导致高程精度很差;生产坡度图方法不尽合理、图斑融合随意致使坡度准确性有欠缺。因此安徽省第三次国土调查考虑重新制作坡度图。 　　从生产流程顺序梳理需要解决的技术问题和重要环节如图1所示。 　　本文重点探讨微小图斑非同值融合到邻接面这一关键技术环节，研究基于arcgis进行数据处理，使用面向对象编程语言调用相关功能控件，进行相应图形对象获取、拓扑、更新等操作。 　　1准备工作 　　1.1dem数据准备 　　本文选择安徽省机载lidar地表高程数据获取制作的dem，该数据点间距2m，通过滤除噪声点、去除架空人工构筑物顶点、鉴别消除高程突变区域等技术手段，生产的dem高程中误差为平地优于0.3m，丘陵地优于0.8m，山地优于1.7m。 　　1.2不同dem计算坡度模型比选 　　本文选择三阶反距离平方差分模型的horm拟合曲面算法，计算的结果，不仅能够很好地反应坡度变化趋势，也可以较好地体现微型地貌和细部表达，是比较理想的计算坡度模型。 　　1.3坡度级栅格矢量化 　　以图幅为单位，将生成的分幅坡度级栅格图进行矢量化，并对相同坡度值的相邻图斑进行合并，得到简化的折线构成的矢量图斑。由于各矢量图斑因原始栅格导致的等距特性，图斑边界点很多在同一直线上，我们以0为容值对其进行抽稀处理，简化图斑边界点。 　　2微小图斑非同值融合到邻接面关键技术探讨 　　非同值融合是坡度图生产的重点和难点，目的是将微小图斑融合到邻接图斑，逐步扩大构成超过3000km2的面，核心是保持微型地貌和细部特征不变形，融合的原则是“逐渐就低”，融合次序（也是其难点）是每次都从该时刻的最小面积图斑（包括多个小面融合的面）入手进行判断，需要时刻保持所有原始的以及融合的面的动态排序和检索。 　　2.1构面线拓扑构建 　　这是进行图斑邻接判断的前提，相邻面拓扑查询判断，是通过构面边界线节的“左侧面id”和“右侧面id”判断面的相邻关系的。 　　arcgis编程环境提供的面拓扑操作组件不尽完善，调用过程往往出现异常，甚至死机，所以，我们从底层直接读取、判断、建立构面线节（数组）。程序读取各面的边界线的各点坐标，记录点与对应的面id，通过统计指向面次数进行判断：如果指向面次数只有1次，说明是外围面上点;超过2次判断为结点;其它正好2次的为线节中间点。结点之间的各点串起线节，线节指向左右两个面，记录线节以及对应的左右两面id，得到构面线节数组。为提高查询效率，程序进行了分块处理，较大面对应多个分块。 　　程序对标准图号的文件名计算内图廓线，将边界线与之重合的面认定为外围边界面，赋予边界不容标志。 　　2.2分段建立过渡数据 　　为了减少判断的范围和数目，提高效率，将生成的过渡数据进行分段管理，排序和融合处理只针对局部的某个特定数量段的数据进行，进行检索变得快捷。分段后可以大大减少排序和判断的数据对象数量，使处理海量数据成为可能。 　　2.3坡度级不同的邻面融合处理 　　首先通过递归找出当前处理面（也可能是融合面）所拥有的所有最终叶子（原始面）;然后对所属的原始面逐个查找相交面（排除自己所属的原始面）;通过相交面的其他面，逐个进行判断，取属性值符合规定的面（按照“就低不就高”原则）id纳入到融合数组中，供迭代判断使用（其个数按照需要动态调整）;最后将融合内容写入到涉及的相关面对象对应的数组中（原始面或融合面）;对融合结