1.一种格状水系识别方法，其特征在于该方法包括：

(1)获取shp格式的待识别水系线图层数据，读入至数据集L＝{li|i＝0,1,2,…,m}，其中，li为表示河流的线要素，m为线要素的个数；

(2)根据数据集L计算得到待识别水系的总长度RL和水系平直特征因子SL；

(3)根据数据集L统计得到各方位角区间内的河流累积长度，统计结果存入列表DL；

(4)根据列表DL和水系总长度RL计算得到水系主要走向占比因子Rratio和水系主要走向夹角因子Rangle；

(5)基于水系平直特征因子SL、水系主要走向占比因子Rratio和水系主要走向夹角因子Rangle进行格状水系类型识别；

(6)对识别出的格状水系绘制走向玫瑰花图。

2.根据权利要求1所述的格状水系识别方法，其特征在于：步骤(2)具体包括：(2-1)针对数据集L中任一要素li，采用以下公式计算其支流长度rli和支流首尾点长度hti：式中，hti表示要素li所代表的支流首尾点长度，rli表示线要素li所代表的支流长度，形如pi,*为线要素li上第*个点，形如 表示两点间的距离，i∈[0,m-1],j∈[0,n-1]，n代表线要素li上点的个数；

(2-2)循环执行步骤(2-1)，直至完成所有线要素的计算；

(2-3)按照以下公式计算整个水系的总长度RL和水系平直特征因子SL：式中，RL表示整个水系的总长，HT表示所有支流的首尾点长度，SL为水系的平直特征因子。

3.根据权利要求1所述的格状水系识别方法，其特征在于：步骤(3)具体包括：(3-1)针对数据集L中任一线要素li，采用以下公式计算得到线要素上各河流分段的方位角angleij：式中，j∈[0,n-1]，形如(xi,*,yi,*)为线要素li上第*个点的坐标，angleij的阈值范围为[0,180)；

(3-2)循环执行步骤(3-1)，直至完成所有线要素的计算；

(3-3)对于计算得到的所有方位角angleij，按10度间隔，分类统计各方位角区间上的河流累积长度，统计结果存入列表DL＝{dlk|k＝0,1,..,17},其中，k为方位角区间序号，其区间范围为[k*10,k*10+10)，dlk为相应区间中各线要素河流分段的累积长度。

4.根据权利要求1所述的格状水系识别方法，其特征在于：步骤(4)具体包括：(4-1)从列表DL中抽取出极大值元素并存入队列MAXL中，当队列MAXL中极大值元素个数小于2时，判定待识别水系为非格状水系，执行步骤(6)；否则，执行步骤(4-2)；

(4-2)从队列MAXL中抽取最大、第二大的两个极大值元素，其在DL中对应方位角区间序号分别记为a、b，其值分别为dla、dlb；

(4-3)针对dla、dlb两个极大值所对应的两个波峰，计算各波峰所包含元素的累加值，分别记为Sa、Sb；其中，计算时不包含波峰的端点元素；

(4-4)根据以下公式分别计算水系主要走向占比Rratio和水系主支流夹角Rangle：Rangle＝|a-b|*10。

5.根据权利要求1所述的格状水系识别方法，其特征在于：步骤(5)具体包括：(5-1)判断基于水系平直特征因子SL、水系主要走向占比因子Rratio和水系主要走向夹角因子Rangle是否满足以下条件：(5-2)若满足条件，则判定待识别水系为格状水系，否则判定为非格状水系。