[发明专利]一种沉陷裂缝区地形DEM修正方法

权利要求书

1.一种沉陷裂缝区地形DEM修正方法，其特征在于，包括以下步骤：

利用无人机获取沉陷裂缝区的厘米级高分辨率光学遥感影像，在遥感影像上勾绘出裂缝的中心线，并量算出裂缝的宽度；

基于裂缝的中心线和宽度数据，利用三角形相似定理计算裂缝中心线两侧每个栅格单元的高程修正值；

将待修正的沉陷裂缝区米级或分米级地形DEM重采样至厘米级分辨率，然后利用计算得到的裂缝中心线两侧栅格单元的高程修正值，修正沉陷裂缝区地形DEM，得到修正后的沉陷裂缝区地形DEM。

2.根据权利要求1所述的一种沉陷裂缝区地形DEM修正方法，其特征在于，所述裂缝中心线两侧栅格单元的高程修正值的计算公式为：

式中，M_(x,y)为坐标为(x,y)栅格单元的高程修正值；W为裂缝的宽度；l为坐标为(x,y)栅格单元到裂缝中心线的垂直距离；α为系数，取值为5～15；β为常数项，取值为0.5～0.9。

3.根据权利要求2所述的一种沉陷裂缝区地形DEM修正方法，其特征在于，当裂缝区域地表为沙时，α取值为5；当裂缝区域地表为土壤时，α取值为13，当裂缝区域地表为岩石时，α取值为15。

4.根据权利要求2或3所述的一种沉陷裂缝区地形DEM修正方法，其特征在于，当裂缝区域地表为沙时，β取值为0.5；当裂缝区域地表为土壤时，β取值为0.7；当裂缝区域地表为岩石时，β取值为0.9。

5.根据权利要求2所述的一种沉陷裂缝区地形DEM修正方法，其特征在于，沉陷裂缝区地形DEM的修正公式为：

Z′_(x,y)＝Z_(x,y)-M_(x,y)

式中，Z′_(x,y)为坐标为(x,y)栅格单元修正后的高程值；Z_(x,y)为坐标为(x,y)栅格单元在修正前地形DEM中的高程值。