[发明专利]基于“窗口法”的流域地貌分形维数量化方法

摘要

本发明公开了一种基于“窗口法”的流域地貌分形维数量化方法，首先输入所需研究的流域地貌的数字高程，确定所需地貌的像元尺寸；然后确定窗口尺寸，将整个流域地貌的数字高程分成若干个相等窗口；提取窗口内中心像元点的高程值与窗口内其它各像元点高程值的均值；然后计算各像元点高程值的均值与中心像元点高程值之间的差值的绝对值，差值称为高程增量值，逐个计算流域地貌的数字高程所分成的若干个相等窗口内的高程增量值，直至整个流域全部计算完成；计算不同窗口尺寸下的流域高程增量绝对值的数学期望；计算得到流域地貌特征分形维数F_D。本发明解决了现有技术中存在的地貌分形计算准确度较差且覆盖不全面的问题。

主权项

1.一种基于“窗口法”的流域地貌分形维数量化方法，其特征在于，具体按照以下步骤实施：/n步骤1、输入所需研究的流域地貌的数字高程，确定所需地貌的像元尺寸；然后依据像元尺寸的奇数倍，确定窗口尺寸R的大小，R大于1，且窗口为正方形；然后将整个流域地貌的数字高程分成若干个相等窗口；/n步骤2、利用ArcGIS软件里的aggregate-tools工具提取窗口内中心像元点(x₀，y₀)的高程值与窗口内其它各像元点高程值的均值；然后计算各像元点高程值的均值与中心像元点(x₀,y₀)高程值之间的差值的绝对值，差值称为高程增量值ΔI(r)，r为各像元到中心像元点(x₀,y₀)的距离；/n步骤3、按照步骤2的方法，逐个计算流域地貌的数字高程所分成的若干个相等窗口内的高程增量值ΔI(r)，直至整个流域全部计算完成；/n步骤4、计算流域所有窗口内高程增量值ΔI(r)绝对值的数学期望E{|ΔI(r)|}，以满足流域地貌分形布朗运动的计算的要求；/n步骤5、改变覆盖窗口R的大小，重复以上步骤，计算不同窗口尺寸R下的流域高程增量绝对值的数学期望E{|ΔI(r)|}；/n步骤6、以高程增量绝对值的数学期望E{|ΔI(r)|}的对数Lg(E{|ΔI(r)|})为纵坐标，窗口内点(x,y)距中心点(x₀,y₀)之间的距离r的对数Lg(r)为横坐标，将相关计算结果点绘在双对数坐标上，利用最小二乘法将数据点拟合出一条直线，这条直线斜率的估计值为豪斯特指数H；/n步骤7，计算得到流域地貌特征分形维数F_D。/n