1.一种国土资源遥测数据管理方法,其特征在于,包括:S1,获取国土空间地理信息的地形数据点集合;
S2,基于所述地形数据点集合,拟合每个地形数据点在预设局部区域 的高程趋势面,得到每个地形数据点对应的第一局部高程趋势面方程;
S3,计算每个地形数据点的实际高程值与第一局部高程趋势面方程计算的对应点高程值之差,得到每个地形数据点的局部高差值;
S4,根据所有地形数据点的局部高差值,计算每个地形数据点的局部高差异常值;
S5,利用局部高差异常值,检测并标记地形异常区域;
S6,提取每个地形异常区域的边界点,拟合每个地形异常区域边界点的局部高程趋势面,得到每个地形异常区域对应的第二局部高程趋势面方程;
S7,利用第二局部高程趋势面方程,计算每个地形异常区域内部地形数据点修正后的高程值,利用修正后的高程值替换原有的高程值。
2.根据权利要求1所述的国土资源遥测数据管理方法,其特征在于:S2,基于所述地形数据点集合,拟合每个地形数据点在预设局部区域 的高程趋势面,得到每个地形数据点对应的第一局部高程趋势面方程,包括:S21,计算每个地形数据点 与对应邻域内点 的高程差异性指标 ;
S22,根据高程差异性指标 ,以地形数据点 为中心,计算局部区域 的半径:;
其中,k为预设的比例系数;
S23,将地形数据点集合划分为多个子区域,每个子区域包含多个局部区域 ;
S24,对每个子区域的局部区域 ,利用kriging算法拟合对应区域的局部高程趋势面方程;
S25,对每个局部区域 ,计算局部高程趋势面方程计算的地形数据点的高程值与实际高程值的残差;
S26,计算每个子区域内所有局部区域 的残差均值 和标准差 ;判断是否满足且 ,如果是,则输出局部高程趋势面方程,如果否,则返回步骤S22减小k值,继续迭代。
3.根据权利要求2所述的国土资源遥测数据管理方法,其特征在于:通过如下公式计算高程差异性指标 :;其中, 和 分别为 和其邻域内 的高程值,n为邻域内点的数量。
4.根据权利要求2所述的国土资源遥测数据管理方法,其特征在于:利用kriging算法拟合对应区域的局部高程趋势面方程:;
其中,(x,y)为估计点的坐标; 为估计点的高程值;m为局部区域内的地形数据点数量; 为第 个已知地形数据点的坐标; 为第 个已知地形数据点的高程值;
为趋势项,用于描述局部区域内高程的总体变化趋势;
为kriging插值系数,满足:
;
;
其中, 为半变异函数,用于描述空间自相关性:;
其中,h为数据点之间的距离,称为空间滞后距离; 为位置x处的高程值,是一个随机变量;E[*]表示数学期望; 为第i个已知点和第 个已知点之间的空间距离; 为估计点和第i个已知点之间的空间距离; 为Lagrange乘子,为了使kriging方程组满足约束条件 而引入的附加变量,在求解方程组时可根据约束条件求得。
5.根据权利要求4所述的国土资源遥测数据管理方法,其特征在于:趋势项 为:
;
其中, 为多项式系数,x和y为 点的坐标。
6.根据权利要求4或5所述的国土资源遥测数据管理方法,其特征在于:S3,计算每个地形数据点的实际高程值与第一局部高程趋势面方程计算的对应点高程值之差,得到每个地形数据点的局部高差值,包括:S31,将每个地形数据点 的坐标 作为输入,利用第一局部高程趋势面方程计算所述点 的高程估计值 :;
S32,计算所述点 的实际高程值 与高程估计值 的差值,得到所述点 的局部高差值 : 。
7.根据权利要求6所述的国土资源遥测数据管理方法,其特征在于:S4,根据所有地形数据点的局部高差值,计算每个地形数据点的局部高差异常值,包括:S41,计算所有地形数据点的局部高差值 的均值 和标准差 :;
;其中,N为地形数据点的总数;
S42,对每个地形数据点 ,计算其局部高差异常值 :。
8.根据权利要求7所述的国土资源遥测数据管理方法,其特征在于:S5,利用局部高差异常值,检测并标记地形异常区域,包括:S51,根据每个地形数据点 的局部高差异常值 ,判断所述地形数据点 是否为地形异常点,若 ,则认为 点的局部高程异常偏高,标记为正异常点;若 ,则认为点的局部高程异常偏低,标记为负异常点;其中, 和 为预设的阈值,用于判断局部高差异常的严重程度;
S52,对所有标记为正异常点或负异常点的地形数据点,根据其空间邻接关系进行聚类,得到一个或多个连通的地形异常区域R;
S53,对每个地形异常区域R,提取其边界点集合 ,其中k为边界点的数量;计算地形异常区域R的几何属性,几何属性包括面积、周长和形心;
S54,根据预设的异常区域筛选条件,对地形异常区域R进行筛选和合并;
S55,对筛选和合并后每个地形异常区域R,计算对应的异常指标 :;
其中, 为异常区域R内部的地形数据点数量;
S56,将异常指标 大于阈值 的地形异常区域R标记为最终的地形异常区域。
9.根据权利要求8所述的国土资源遥测数据管理方法,其特征在于:S54,根据预设的异常区域筛选条件,对地形异常区域进行筛选和合并,包括:将面积小于 或大于 的区域,判断为虚拟异常区域,予以删除;
将质心距离小于 的区域,如果均为正异常区域或负异常区域则合并为一个区域;
其中, 、 和 为预设阈值。
10.一种国土资源遥测数据管理系统,其特征在于,包括:
至少一个处理单元;用于执行指令以实现权利要求1至9任一项所述的国土资源遥测数据管理方法。