地理信息系统 第五章 地理信息系统的空间分析功能1.docxVIP

147

5．2．4坡度、坡向分析

据空间解析几何的原理，基本单位长度为Δx,Δy,原点坐标为

→Pi,j(x0,y0)的格网模型(图5-12(a))任一格网交点的标准向量为

→

Pi,j

（图

根

(a)(b)图5-12格网模型空间向量分析示意图

(i=1,2,...,M)(5-11)

(j=1,2,...,N)

→→

对于由四个相邻格网点确定的地表基本单元，其基本向量a,b的计算公式如下：

→→→

ai,j=Pi+1,j+1—Pi,j=(Δx,Δy,zi+1,j+1—zi,j)

→→→bi,j=Pi,j+1—Pi+1,j=(—Δx,Δy,zi,j+1—zi+1,j)

→→→

(5-12)

→

→→a,b

→

→→a,b

向量的

→

向量积就是基本单元的法向量ni,j。

148

(5-13)—Δx(zi,j+1+zi+1,j+1—zi+1,j—zi,j),2ΔxΔy)(i=1,2,...,M)

(5-13)

(j=1,2,...,N)

→

利用ni,j，就可以进行地表单元坡度、坡向等的分析和计算。

1．坡度分析

→

格网模型中地表基本单元的坡度等于其法向量ni,j与z轴之夹

角（图5-13中的slopei,j），而两向量的夹角余弦等于两向量的数量积与模的乘积之商，即：

图5-13地表基本单元坡度、坡向示意图图5-14格网点配置示意图

→→

slopei,j=arccos()

=arccos(2ΔxΔy/((Δy(zi,j+1+zi,j

(5-14)

—zi+1,j+1—zi+1,j))2

1

+(Δx(zi,j+1+zi+1,j+1—zi+1,j—zi,j))2+4Δx2Δy2)2)

根据(5-14)式计算格网模型的坡度的过程较为复杂。一般情况下，如果格网为正方形，那么可以采用下面介绍的简化公式（参见图5-14）进行计算。

假设za,zb,zc,zd为地表基本单元四个格网点的高程数据，Δxy为格

149

网的基本单位长度，那么P点的微分公式为：

(5-15)

(5-15)

对(5-14)、(5-15)而言，0o≤slopei,j≤90o。

具体应用时，可根据需要对度数进行分级，以形成坡度分析的分级标准。

当需要时，也可以把度数转化为百分比。

2．坡向分析

在图5-13中，x轴的方向向南，因此，地表基本单元的坡向即

→n

→

np

为其法向量ni,j在oxy平面上的投影

与x轴的夹角（结合图5-15）。

计算公式为：

(5-16)

(5-16)

=tan一1(Δx●A●B(j))

在(5-16)中，A(j)=zi,j+zi+1,j一zi,j+1一zi+1,j+1

在(5-16)中，

B(j)=zi+1,j+zi+1,j+1一zi,j一zi,j+1

图5-15坡向变量示意图

按式(5-16)计算出的值是与正南方向的夹角。在实际应用中，如果所建立的DEM数据是按从南到北获取的，那么在采用上述公式时，求出的DIR是与正北方向的夹角。

另外，计算出的坡向DIR有与x轴正向和负向夹角之分。坡向的正负号可以通过A(j),B(j)的符号加以确定（见表5-1）。

150

从表5-1可以看出，采用上述方法得到的坡向分级比较详细。在实际应用过程中，需要给予综合，得到四种坡向（图5-17）:平缓坡，阳坡，半阳坡，阴坡。在图5-17中，这四种坡向分别用1、2、3、4加以表示。

图5-17坡向的综合

与计算坡度类似，在进行具体的坡向分析时，同样可以采用如式(5-15)所示的简化公式，即计算出U,V后，可以按照下列规则直接确定不同的坡向（图5-17(b)）。

(1)当U0，在tan一1中，

151

如果V=0，则θ=0o，DIR为E；

如果V0，则