专利名称：一种河网密度空间化方法
技术领域：
本发明公开了一种河网密度空间化方法，属于气象灾害风险综合区划技术领域。
背景技术：
在热带气旋、暴雨洪涝的孕灾环境分析中，河网水系的区域性分布特征是至关重要的影响因素，在很大程度上决定了评价区域遭受洪涝侵袭的难易程度。河网密集程度和河流级别是河网水系的区域性分布特征的两个重要因素，例如在河槽汇流超过河道及水库的蓄洪和排泄能力以后，水流就会向河道周边蔓延、泛滥，假若该地区河网密布且均发生漫堤现象时，就会形成严重且难以抵御的洪涝灾害；传统的方法假定不同级别的河系的影响范围和致灾效能相同，对区域河网密度进行分析和计算，从而实现河网密度的空间化，而实际情况中干流要较一级支流、一级支流较二级支流具有更强的影响范围和致灾效能，所以传统方法对河网密度的空间化具有局限性。
发明内容
为了实现对区域河网密度的空间化，避免以往方法的局限性，提高方法的合理性和适用性，本发明的目的是提供一种利用层次分析法，根据不同等级河系影响范围和致灾效能的不同，确定各等级河系的权重系数，然后通过面积加权法获得单位网格内的河网缓冲区面积，即河网密度，来实现对河网密度空间化方法。本发明为实现上述发明目的采用如下技术方案，一种河网密度空间化方法，其特征在于，包括如下步骤:步骤(I)利用GIS软件中的网格划分功能对河系所在区域进行网格划分；步骤(2)引用层次分析法确定不同等级河系的权重系数Wi ；步骤(3)利用GIS软件中的缓冲区分析功能，设置缓冲区半径Ri,分别对不同等级的河系进行缓冲区分析，生成代表该河系的影响范围和致灾效能的河系缓冲区；步骤(4)利用GIS软件中的叠加分析功能，分别对不同等级河系的缓冲区和网格区域进行求交分析，利用网格单元打断河系缓冲区，即连续的一条河系缓冲区被不同的单位网格分割成若干段；步骤(5)在进行过求交分析后，从生成的结果数据集中提取相应的网格单元中被打断的河系缓冲区段的面积字段信息area，不同等级河系求交分析后生成的面积字段分别记作areai，通过GIS软件中的数据集追加列功能把area字段分别赋予给网格区域数据集；步骤(6)在网格区域数据集中新建河网密度字段，结合GIS软件中的更新列功能，利用面积加权法，计算出单位网格内的河网缓冲面积，即河网密度M ；上述步骤中，i代表河系等级，i的取值范围是自然数。进一步的，步骤(2)中所述不同等级河系的权重系数Wi的具体计算方式为:
(201)构造判断矩阵:判断矩阵中的元素axy采用1_9及其倒数的标度方法，其中x、y均表示河系等级，x、y的取值范围是自然数，axy表示第X级河流相对于第I级河流的比较结果，且axy=l/ayx，相互比较不同等级河流的重要性，用数字化的比值来定义比较结果；(202)根据判断矩阵，计算判断矩阵的特征根和特征向量:对判断矩阵A，计算矩阵满足AW= λ Maxff的最大特征根与特征向量，式中λ Max为判断矩阵A的最大特征根，W为对应于λ Max的正规化特征向量，最大特征根λ Max对应的特征向量归一化为=W=WnW2,…，Wi,其中Wi即为各等级河系的权重系数，i代表河系等级，i的取值范围是自然数；(203) 一致性检验:一致性检验涉及一致性指标C.1.，平均随机一致性指标R.1.，一致性比例c.R.，一致性检验的步骤如下:(2031)计算一致性指标C.1.，
权利要求
1.一种河网密度空间化方法，其特征在于，包括如下步骤: 步骤(I)利用GIS软件中的网格划分功能对河系所在区域进行网格划分； 步骤(2)引用层次分析法确定不同等级河系的权重系数Wi ； 步骤(3)利用GIS软件中的缓冲区分析功能，设置缓冲区半径Ri,分别对不同等级的河系进行缓冲区分析，生成代表该河系的影响范围和致灾效能的河系缓冲区； 步骤(4)利用GIS软件中的叠加分析功能，分别对不同等级河系的缓冲区和网格区域进行求交分析，利用网格单元打断河系缓冲区，即连续的一条河系缓冲区被不同的单位网格分割成若干段； 步骤(5)在进行过求交分析后，从生成的结果数据集中提取相应的网格单元中被打断的河系缓冲区段的面积字段信息area，不同等级河系求交分析后生成的面积字段分别记作areap通过GIS软件中的数据集追加列功能把area字段分别赋予给网格区域数据集； 步骤(6)在网格区域数据集中新建河网密度字段，结合GIS软件中的更新列功能，利用面积加权法，计算出单位网格内的河网缓冲面积，即河网密度M ； 上述步骤中，i代表河系等级，i的取值范围是自然数。
2.如权利要求1所述的一种河网密度空间化方法，其特征在于，步骤(2)中所述不同等级河系的权重系数Wi的具体计算方式为: (201)构造判断矩阵:判断矩阵中的元素axy采用1-9及其倒数的标度方法，其中x、y均表示河系等级，x、y的取值 范围是自然数，axy表示第X级河流相对于第y级河流的比较结果，且axy=l/ayx，相互比较不同等级河流的重要性，用数字化的比值来定义比较结果； (202)根据判断矩阵，计算判断矩阵的特征根和特征向量:对判断矩阵A，计算矩阵满足AW= λ Maxff的最大特征根与特征向量，式中λ Max为判断矩阵A的最大特征根，W为对应于入Max的正规化特征向量，最大特征根λ Max对应的特征向量归一化为=W=WpW2,…，Wi,其中Wi即为各等级河系的权重系数，i代表河系等级，i的取值范围是自然数； (203)一致性检验:一致性检验涉及一致性指标C.1.，平均随机一致性指标R.1.，一致性比例C.R.，一致性检验的步骤如下: (2031)计算一致性指标C.1.C.1.=,其中η为矩阵阶数； n~l (2032)确定相应的平均随机一致性指标R.1.; (2033)计算一致性比例C.R.并进行判断 (1=(上 RJ.当C.R.〈0.1时，认为判断矩阵的一致性可以接受的， 当C.R.>0.1时，认为判断矩阵不符合一致性要求，需要对该判断矩阵进行重新修正。
3.如权利要求1所述的一种河网密度空间化方法，其特征在于，步骤(6冲所述河网密度M，其具体计算方式为: 通过公式M=Srea1 X Wjarea2 X ff2+area3 X W3+…+Sreai X Wi,计算出单位网格内的河网缓冲面积，即河网密度[M=O^M2,…Mt)]，其中，t为网格单元编号。
4.如权利要求1所述的一种河网密度空间化方法，其特征在于:所述步骤(I)中对河系所在区域进行网格划 分，每个网格单元的划分标准为IOOmX 100m。
全文摘要
本发明公开了一种河网密度空间化方法，它首先对区域进行网格划分，其次利用层次分析法，根据不同等级河系的影响范围和致灾效能的不同，确定不同等级河系的影响权重系数，结合GIS(地理信息系统)软件功能获取单位格网内不同等级河网的缓冲面积，然后通过面积加权法计算出空间网格河网密度。该技术方案避免以往方法的局限性，能很好地反映不同等级河系的影响范围和致灾效能，具有较高的学性和合理性。
文档编号G06F19/00GK103226662SQ201310146468
公开日2013年7月31日 申请日期2013年4月24日 优先权日2013年4月24日
发明者薛丰昌, 钱洪亮, 耿焕同, 盛洁如 申请人:南京信息工程大学