自然水系(或分水岭)自动提取的matlab实现方案3页.doc

水系提取算法（D8）的matlab实现D8算法是当今非常成熟的提取水系（或分水岭）的计算机实现方法。这里笔者结合自身的编程实践介绍D8算法的matlab实现过程，谨与大家分享。D8算法分为两部分实现：每个栅格点水流方向的计算；每个栅格点汇水面积的计算。1、水流方向的计算在水流方向的计算中，本人结合matlab内嵌的滤波工具，利用一个简单的非线性滤波来实现，具体内容包括一个自己编辑的确定水流方向的函数flowdirection和若干脚本命令，如下：以下是自己编辑的flowdirection函数代码：147258369function output=flowdirection(a)n=size(a);for i=1:n(2) k=0; b=-inf; for j=1:4,6:9 if rem(j,2)=0 r(j)=a(5,i)-a(j,i); else r(j)=(a(5,i)-a(j,i)/