【GIS实验】山脊山谷鞍部提取

*用于ArcGIS山脊山谷鞍部提取，速提备注，不做详细解释。

一、正负地形提取（用于校正）

1. Spatial Analyst——领域分析——焦点统计（输出：meandem）

   •   11*11 mean

2. Spatial Analyst——地图代数——栅格计算器（输出：jianfa）

   •   "DEM" - "meandem"

3. Spatial Analyst——重分类——重分类（输出：zhengdixing，fudixing）

   •   正地形：>0→1，<0→0

   •   负地形：>0→0，<0→1

二、山脊线

1. Spatial Analyst——水文分析——填洼（输出：filldem）

2. Spatial Analyst——水文分析——流向（输出：flowdirfill）

3. Spatial Analyst——水文分析——流量（输出：flowacc1）

4. Spatial Analyst——地图代数——栅格计算器（输出：facc0）

   •   "flowacc1" == 0

5. Spatial Analyst——领域分析——焦点统计（输出：neiborfacc0）

   •   3*3 mean

6. Spatial Analyst——表面分析——等值线（输出：denggao）

7. Spatial Analyst——表面分析——山体阴影（输出：yinying）

8. 属性——符号系统——间断点0.5541

9. Spatial Analyst——重分类——重分类（输出：shanji1）

   •   >0.5541→1，其他→0

10. Spatial Analyst——地图代数——栅格计算器（输出：shanji2）

    •   "zhengdixing" * "shanji1"

11. Spatial Analyst——重分类——重分类（输出：shanjixian）

    •   ≠1→NO DATA

三、山谷线

1. Spatial Analyst——地图代数——栅格计算器（输出：fandem）

   •   Abs("DEM" -  5000")

2. Spatial Analyst——水文分析——流向（输出：flowdirfan）

3. Spatial Analyst——水文分析——流量（输出：flowacc2）

4. Spatial Analyst——地图代数——栅格计算器（输出：facc0fan）

   •   "flowacc2" == 0

5. Spatial Analyst——领域分析——焦点统计（输出：nbfacc0fan）

   •   3*3 mean

6. 属性——符号系统——间断点0.65677

7. Spatial Analyst——重分类——重分类（输出：shangu1）

   •   >0.65677→1，其他→0

8. Spatial Analyst——地图代数——栅格计算器（输出：shangu2）

   •   "fudixing" * "shangu1"

9. Spatial Analyst——重分类——重分类（输出：shanguxian）

   •   ≠1→NO DATA

四、鞍部

1. Spatial Analyst——地图代数——栅格计算器（输出：anbuqu）

   •   "flowacc0" * "facc0fan"

2. Spatial Analyst——地图代数——栅格计算器（输出：rasteranbu）

   •   "zhengdixing" * "anbuqu"

3. Spatial Analyst——重分类——重分类（输出：rasteranbu2）

   1.   =0→NO DATA

4. 转换工具——由栅格转出——栅格转点（输出：anbudian）

5. 对照等高线和山体阴影人工剔除伪鞍部点