立地类型划分

立地类型划分

（二）立地类型划分

1、坡向图制作

（1）坡向初级分布图生成

打开ArcMap程序，并加载调查区（马连滩）DEM数据。打开ArcToolbox，选择Spatial Analyst工具→表面分析→坡向（图23），生成马连滩坡向初步分级分布图（图24）。

图23 坡向提取工具图24 马连滩坡向提取结果

（2）坡向分类

利用中重分类的功能（重分类）功能对上一步生成的坡向分布结果进行分类。选择Spatial Analyst工具→重分类→重分类工具，弹出对话框（图25）；点击对话框中“分类”按钮，弹出对话框（图26），根据“坡向方位界定标准”（表

3），对调查区坡地的坡向进行分类，先重分类为3类：0-112.5°、112.5-292.5°、292.5-360°，点击确定。右键重分类后的图层，点击“打开属性表”→“添加字段”，添加字段（如坡向），字段类型选择为“文本”。根据表3，0-112.5°和292.5-360°为阴坡，112.5-292.5°为阳坡，在新添加的字段中输入阴坡或阳坡（图27），得到重分类结果（图28）。

表3 坡向方位界定标准

坡向Aspect

方位Orientation 方位角azimuth angle（°）

阳坡Sunny Slope

南South 157.5 - 202.5 东南Southeast 112.5 - 157.5 西南Southwest 202.5 – 247.5 西west 247.5 - 292.5

阴坡Shady Slope

北North 337.5 - 22.5 西北Northwest 292.5 - 337.5 东北Northeast 22.5-67.5 东East 67.5- 112.5

图25打开重新分类工具

图26分类方法确定

图27坡向属性编辑图28坡向重新分类结果（3）去除小斑块

划分结果中存在大量面积很小，在实际生产中没有意义的斑块，可合并到邻近的大斑块中。

在ArcToolbox中，按照路径：Spatial Analyst 工具→栅格综合→众数滤波（图29），打开“众数滤波”对话框（图30）。将上一步重分类的结果进行输入，“要使用的相邻要素数”选项中选

择“FOUR”，“替换阈值”选项中选择“MAJORITY”，点击确定。查看属性表（打开属性表）中“坡向”字段，若没有，则需重新添加（添加字段）。

图29 图30

（4）矢量化

将去除小斑块的坡向图矢量化，在ArcToolbox中按照路径：转换工具→由栅格转出→栅格转面，打开栅格转换矢量（栅格转面）工具。指定输入输出路径；选中“简化面”单选框；“字段”选项中选择“poxiang” （图31）。确

定输出矢量（图32）

图31 栅格转换矢量工具图32 转换成矢量结果显示

2、坡度图制作

（1）坡度因子初步提取

在ArcMap中，加载调查区（马连滩）DEM 数据。打开ArcToolbox，选择Spatial Analyst 工具→表面分析→坡度工具（图33），生成马连滩坡度初步分级分布图（图34）。

图33 坡度提取工具图34 坡度提取结果

（2）坡度分级分级

利用重分类的功能（重分类）对上一步生成的坡度分布结果进行分级。选择Spatial Analyst 工具→→工具，弹出对话框（图35）；点击对话框中“分类”按钮，弹出对话框（图36），按照“坡度因子等级表”（表4），对上步初步提取的坡向进行分级，此处分为2类。右键重分类后的图层，点击“打开属性表”→“添加字段”，添加字段（如坡度），字段类型选择为“文本”（图37），得到分级后的结果（图38）。

图35 坡度分级工具

图36 坡度重新分级设置

表4 立地类型划分坡度因子等级

坡度等级Slope

坡度Slope Degree (°) Ranking

平缓坡Gentle 0 - 15

斜坡Incline 16 - 25

陡坡Steep 26 - 35

暂不可造林地

Temporarily

> 35

unavailable

afforestation land

图37 坡度属性编辑图38 坡度分级结果

（3）去除小斑块

将重分类的结果去除小斑块，在ArcToolbox 中按照路径：Spatial Analyst工具→栅格综合→众数滤波，打开“众数滤波”对话框（图39）。选择重分类的结果进行输入“要使用的相邻要素数”选项中选择“FOUR”，“替换阈值”选项中选择“MAJORITY”，点击确定。查看属性表（打开属性表）中“坡度”字段，若没有，则需重新添加（添加字段）。

图39

（4）矢量化

在ArcToolbox中按照路径：转换工具→由栅格转出→栅格转面，打开栅格转换矢量（栅格转面）工具（图40）。指定输入输出路径；选中“简化面”单选框；“字段”选项中选择“podu”。确定，输出矢量（图41）。

图40 栅格转换矢量工具图41 坡度图层栅格转矢量结果

3、坡位图制作

（1）山脊沟底线的提取

①提取坡向数据

点击DEM数据，使用表面分析工具中的坡向（Aspect）工具，提取DEM的坡向数据层，命名为A。

②提取坡度数据

点击数据层A，使用表面分析中的坡度（Slope）工具，提取数据A的坡度数据，命名为SOA1。

③生成反地形DEM

使用空间分析工具中的栅格计算器，由路径Spatial Analyst工具→地图代数→栅格计算器，用最大高程值H（此处H=1376.93）减去DEM，公式为:1335.57-“DEM”，得到与原来地形相反的数据层，即反地形DEM，命名为R-DEM。

④基于反地形DEM数据求算坡向值

利用反地形DEM数据，并使用空间分析工具的Aspect工具，求算R-DEM的坡向值。

⑤求算反地形坡向变率

利用SOA法求算反地形的坡向变率，对上一步求出的R-DEM坡向值再求算坡度值Slope，即得到反地形的坡向变率，记为SOA2。

⑥消除坡向变率误差

使用空间分析工具中的栅格计算器，Spatial Analyst工具→地图代数→栅格计算器，公式为SOA=((“SOA1”+“SOA2”)-Abs(“SOA1”-“SOA2”))/2，即可求出没有误差的坡向变率SOA。

⑦求算坡向平均值

再次点击初始DEM数据，使用空间分析工