网格数据怎么提取山脊线

Posted 2023-03-27

参考技术A

山脊线的提取实质上也是分水线与汇水线的提取。因此，可以利用水文分析的方法进行提取。对于山脊线而言，由于它同时也是分水线，而分水线的性质即为水流的起源点。所以，通过地表径流模拟计算之后，这些栅格的水流方向都应该只具有流出方向而不存在流人方向，即概格的汇流累积量为零。因此，通过对零值的提取，就可得到分水线，即山脊线。

（1）填充洼地。点击【Spatial Analyst Tools】-【Hydrlogy】-【Fill】，在打开的Fill对话框中选择输入dem数据，输出栅格数据文件名为Fill_dem。

（2）计算水流方向。选择【Spatial Analyst Tools】-【Hydrology】-【Flow Direction】，在打开的Flow Direction对话框中选择输入Fill_dem数据，输出栅格数据文件名为flowdirfill。

（3）计算汇流累积量。点击【Spatial Analyst Tools】-【Hydrology】-【Flow Accumulation】，在弹出的Flow Accumulation对话框中选择输入flowdirfill数据，输出栅格数据文件名为flowacc1。

（4）汇流累积量为零值的提取。打开【Spatial Analyst Tools】-【Map Algebra】-【Raster Calculator】,在打开的Raster Calculator对话框中输入公式："flowacc1" == 0，输出栅格数据文件名为facc0。

（5）对facc0进行领域分析。点击【Spatial Analyst Tools】-【Neighborhood】-【Focal Statistics】，在打开的Focal Statistics对话框中选择输入facc0数据，输出栅格数据文件名为neiborfacc0。

（6）打开【Spatial Analyst Tools】-【Surface】-【Contour】和【Spatial Analyst Tools】-【Surface】-【Hillshade】，分别生成dem的等值线图和晕渲图hillshade。

（7）在neiborfacc0图层上右键选择属性，进行重新分级，分为2级，这时不断调整分级临界点，并以等值线图和晕渲图作为辅助判断，属性值越接近1的栅格越有可能是山脊线的位置，最终确定的分界阈值为0.5541。

（8）将进行过二值化的neiborfacc0进行重分类为Reclass_neiborfacc0，将属性值接近1的那一类的属性值赋值为1，其余赋值为0。

（9）打开栅格计算器，输入命令: "Reclass_neiborfacc0" * "zhengdixing"，以消除那些存在于负地形区域中的错误的山脊线。输出栅格数据文件名为：shanjixiannoworry。

（10）对shanjixiannoworry图层重分类，所有属性不为1的栅格属性赋值NO DATA，这样就得到了山脊线。