[發明專利]一種基于等高線拉伸的數字高程模型自動調整方法及系統

有效

• 申請號: 202111386931.6
• 申請日: 2021-11-22
• 公開（公告）號: CN113836751B
• 公開（公告）日: 2022-02-08
• 發明（設計）人: 段玲娟;吳佳盈
• 申請（專利權）人: 武漢峰嶺科技有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;G01C11/00
• 代理公司: 武漢謙源知識產權代理事務所(普通合伙) 42251
• 代理人: 王力

1.一種基于等高線拉伸的數字高程模型自動調整方法, 其特征在于，包括如下步驟：

S1：加載航測立體像對和包含有等高線的數字高程模型，在所述數字高程模型中選定目標調整區域，并根據所述目標調整區域選定拉伸矢量；其中，所述目標調整區域內包含有所述等高線，所述等高線由多個節點順次連接形成；

S2：根據所述拉伸矢量計算所述目標調整區域內的所述等高線的形變量，并根據所述形變量對所述目標調整區域內的所述等高線進行形變，得到新的等高線；

S3：根據形變后的所述等高線調整數字高程模型中等高線上特征點的高程；

其中，所述拉伸矢量由所述目標調整區域內所有等高線的平均方向和所述目標調整區域內等高線上的節點投影至航測立體像對中的左影像和右影像上抽取的兩個像素塊的相關系數確定。

2.根據權利要求1所述的基于等高線拉伸的數字高程模型自動調整方法, 其特征在于，所述根據所述目標調整區域選定拉伸矢量包括：

S11：在所述目標調整區域內或目標調整區域外選定拉伸起點A，確定拉伸方向；

S12：沿著所述拉伸方向匹配拉伸終點B，所述拉伸起點A和拉伸終點B形成所述拉伸矢量；

其中，所述拉伸起點與拉伸終點之間的連線或二者連線的反向延長線貫穿所述目標調整區域。

3.根據權利要求2所述的基于等高線拉伸的數字高程模型自動調整方法, 其特征在于，所述確定拉伸方向具體包括如下步驟：

S111：在所述目標調整區域每條所述等高線上等間隔選取所述節點，依次計算每條所述等高線上相鄰兩個所述節點形成的線段的方向；

S112：對所有線段的方向取平均值，得到所述等高線的平均方向；

S113：將所有等高線的方向取平均值，得到所述目標調整區域內所有等高線的平均方向；

S114：將所述目標調整區域內所有等高線的平均方向旋轉90度，得到所述拉伸方向。

4.根據權利要求2所述的基于等高線拉伸的數字高程模型自動調整方法, 其特征在于，所述沿著所述拉伸方向匹配拉伸終點B具體包括如下步驟：

S121：給定初始拉伸長度len 0，并計算所述初始拉伸 長度len 0下所有所述等高線的變形匹配值f₀；

S122：給定初始步距s₀，進行循環迭代計算當前拉伸長度len i，并計算當前拉伸長度下的變形匹配值f_i，當前拉伸長度len i的計算公式為：

其中，i表示初次迭代次數，且0＜i＜n，n為迭代總次數；

S123：取所述變形匹配值f_i中的最大值f_m，給定迭代步距s_j= s_j-1/2，并進行迭代計算當前拉伸長度len i，計算當前拉伸長度len i長度下的變形匹配值f_i，當前拉伸長度len i的計算公式為：

其中，j表示二次迭代次數，且0＜j＜n，n為迭代總次數，len m表示變形匹配值f_m對應的拉伸長度；

S124：重復所述S123，直至所述迭代步距s_j＜toler，結束迭代結算，確定當前變形匹配值f_i中的最大值f_m對應的拉伸長度len m為最佳拉伸長度，并確定所述拉伸終點B，其中，toler為預設迭代步距閾值。