申请/专利权人：西北农林科技大学

申请日：2019-12-26

公开（公告）日：2021-11-02

公开（公告）号：CN111177917B

主分类号：G06F30/20(20200101)

分类号：G06F30/20(20200101)

优先权：

专利状态码：有效-授权

法律状态：2021.11.02#授权;2020.06.26#实质审查的生效;2020.05.19#公开

摘要：本发明提出了一种基于SRTM的坡长提取方法，通过建立地理坐标下栅格模型，从该模型中直接推导出SRTM栅格边长的计算方式。计算坡长时考虑坡度截断、河网截断影响，使结果尽可能符合实际情况，最后使用正反遍历计算SRTM的累计坡长。与传统方法相比，本方法不再经过投影坐标变换，从而提高坡长的提取效率；同时也能保证原始数据的精度不下降，坡长提取结果相对精确，坡长结果与现有的DEM计算结果对比显示，坡长差值范围集中，具有规律性。

主权项：1.一种基于SRTM的坡长提取方法，其特征在于，该方法按照以下步骤进行：步骤1：存储SRTM数据头信息，同时获取SRTM数据中每一个栅格的坡度值和流向；步骤2：计算栅格单元坡长值；步骤2.1：遍历所有栅格，判断当前栅格是否是无值点：若是无值点：则设置该栅格的单元坡长值为0，进入下一个栅格的判断；若是非无值点，执行步骤2.2；步骤2.2：根据步骤1获取的栅格流向，若当前栅格的流向为东西方向，该当前栅格的单元坡长为LONG；若流向为南北方向，该当前栅格的单元坡长为LAT；若流向为东南、西南、东北、西北中的任意一种，则该当前栅格的单元坡长为diagCellSize；其中，LAT＝2π×6371000·cellsize360，LONG＝2π×6371000·cellsize360·cosθ，θ为当前栅格中任意一点与地球球心连线与赤道的夹角； 其中，cellsize为SRTM数据中每个栅格对应的头信息中cellsize这一项的取值；步骤2.3：遍历所有栅格，得到每一个栅格的单元坡长值；步骤3：计算栅格初始坡长值；步骤3.1：判断当前栅格是否截断，若为不截断，则当前格栅初始坡长值等于单元坡长值；若为截断，则当前栅格的初始坡长值为单元坡长值的一半；步骤3.2：遍历所有栅格，得到每一个栅格的初始坡长值；步骤4：计算栅格累积坡长：步骤4.1：将当前栅格作为初始累积栅格，初始累积栅格初始值设置为0，任意与初始累积栅格相邻的栅格中流向指向初始累积栅格的定义为相邻指向栅格，若其中任意一个相邻指向栅格截断，则取初始值与该相邻指向栅格的初始坡长值的一半之和作为初始累积栅格对该相邻指向栅格的累积坡长；若该相邻指向栅格不截断，则取初始值与该相邻指向栅格的初始坡长值之和作为初始累积栅格对该相邻指向栅格的累积坡长；依次计算初始累积栅格对其他相邻指向栅格的累积坡长，比较初始累积栅格与所有相邻指向栅格的累积坡长的总和SUM与初始累积栅格的初始坡长值的大小，若SUM＞初始累积栅格的初始坡长值，则将SUM赋值给初始累积栅格作为初始累积栅格的初始坡长值，否则初始累积栅格的初始坡长值不变；步骤4.2：将所有剩余栅格中的每个栅格均作为初始累积栅格，按照步骤4.1的方法，遍历所有剩余栅格，计算每个栅格对应的SUM并判断是否赋值，至不再出现赋值情况时遍历结束，最后一次赋值时获得的初始坡长值即为提取的对应栅格的累积坡长；所述步骤1中，存储SRTM数据的头信息，同时获取SRTM数据中每一个栅格的坡度值和流向，包括：步骤1.1：创建结构体DemData，用于存放SRTM的头信息cellsize；步骤1.2：按行读取SRTM的数据，SRTM的头信息的格式为：前边是名称，接着是空格，后边是头信息中cellsize这一项；将读取的每行数据存到对应的一个字符串中，然后对该字符串进行空格分割，最后将头信息转换类型并保存到数据结构DemData对应的属性中；步骤1.3：采用最陡坡降法计算每个单元栅格的坡度值和流向。

全文数据：

权利要求：

百度查询： 西北农林科技大学 一种基于SRTM的坡长提取方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD