模型图.docxVIP

地貌类型划分模型问题重述：已知：地形起伏度是单位面积内最大相对高程差，反映地面相对高差，是描述地貌形态的定量指标。地形起伏度在土地利用评价、土壤侵蚀敏感性评价、生态环境评价、人居环境适宜性评价、地貌制图、地质环境评价等领域有广泛应用。使用不同尺度的地形起伏度会影响相关研究的基本结论。利用京津冀地区1:25万DEM数据按照地貌类型，计算获取京津冀地貌类型的最优分类方案，并回答如下问题：1、如图1，在正方形滤波、圆形滤波、环形滤波、扇形滤波中，你觉得应该选取哪种形态的滤波因子计算地形起伏度？2、选定滤波因子形态后，你觉得应该选取多大尺度的滤波因子（如应该选取圆形滤波因子的半径是多少为最佳？），才能获得地形起伏度的最佳计算结果，最佳尺度是多少？利用计算所得的最佳尺度地形起伏度编制专题地图。图1、不同形态的滤波因子3、利用计算得到的最佳尺度的地形起伏度数据集和原始提供的DEM海拔高度数据集，依据表1的地貌类型分类标准，测算京津冀地貌类型，并编制地貌类型专题地图。问题分析：从表1的地貌类型划分标准，我们可以知道在海拔高度一定时，地形起伏度可以作为划分地貌形态的重要指标。根据地形起伏度的定义可以知道计算地形起伏度的关键在于搜索区域的内的最高点和最低点，随着分析区域范围的增大，其高差也会发生一定变化，进而地形起伏度的值也随之发生变化，因此确定一个最佳分析区域是地形起伏度提取算法中的核心步骤。我们知道，按标准划分时，地貌类型与实际情况最相符时，即为最优分类。此时所对应的滤波因子的尺度即为最适尺度，形状即为最适形状。通过比较不同滤波因子下地形起伏度专题地图，得出与实际情况最相符的地图（所对应的滤波因子即为最佳尺度r）。2.利用计算所得的最佳尺度地形起伏度数据集和原始提供的DEM海拔高度数据集，依据表一的地貌类型分类标准，测算京津冀地貌类型，并编制地貌类型专题地图。基本假设假设一：题中所给的数据真实可靠，能反映京津冀地区每点的高程值。假设二：当统计单元中的中心点为无效点时，其所在的统计单元无效。（即不考虑中心点的地势起伏度）假设三：电脑等设备造成的误差忽略不计。假设四：滤波因子尺度仅为奇数不影响实验结果。假设五：当统计单元中包含无效点时（中心点除外），其对求统计单元的地势起伏度的值的影响忽略不计。假设六：符号说明模型的建立与求解（1）（模型一）不同尺度的滤波因子下的面积比（方形邻域）1.不考虑以无效数据为中心的邻域if MyDem(i+a-1,j+a-1)~=-9999M=MyDem(i:i+r-1,j:j+r-1);end2.去除不是以无效数据为中心的邻域中的无效值temp=-9999;%寻找出邻域中第一个有效值for x1=1:r-1for x2=1:r-1if M(x1,x2)~=-9999 temp=M(x1,x2); break;end end if temp~=-9999 break;end end%将邻域中的无效值改为有效值temp for x1=1:r-1 for x2=1:r-1 if M(x1,x2)==-9999 M(x1,x2)=temp;end end end3.判断邻域的地貌类型if big1200count_yuan=count_yuan+1;elseif R=30if big=100count_ping=count_ping+1;elsecount_pen=count_pen+1;endelseif R=200count_qiu=count_qiu+1;elseif R=500count_di=count_di+1;elseif R=1000count_zhong=count_zhong+1;elsecount_gao=count_gao+1;endendendendend（2）（模型二）不同尺度的滤波因子下的地势起伏度（圆形邻域）圆形邻域模型实在方形邻域模型的基础上更改的，如图圆形模型只是将方形模型的四个角出去，所以只需在方形算法的基础上加以下四条语句即可。 M(1,1)=M(1,2);M(1,r-1)=M(1,r-2); M(r-1,1)=M(r-1,2); M(r-1,r-1)=M(r-1,r-2);在最佳尺度下的地势起伏度：R=Hmax-Hmin通过对比，确定最佳滤波因子形态确定最佳尺度确定最佳尺度方形滤波：不同滤波因子下的平均地势起伏度圆形滤波：不同滤波因子下的平均地势起伏度不同形态的滤波因子研究区概况：京津冀位于处于环渤海地区和东北亚的核心重