一种自然地表特征表征方法.pdf

1、10申请公布号CN104318066A43申请公布日20150128CN104318066A21申请号201410508901122申请日20140926G06F19/0020110171申请人武汉大学地址430072湖北省武汉市武昌区珞珈山武汉大学72发明人刘耀林刘艳芳焦利民何建华74专利代理机构武汉科皓知识产权代理事务所特殊普通合伙42222代理人薛玲54发明名称一种自然地表特征表征方法57摘要本发明公开了一种自然地表特征表征方法，包括自然地形评价、自然地貌评价、自然地表覆盖评价、自然地理国情综合评价、自然地表格局特征表征；本发明使用综合评估模型对区域地理国情要素中的自然地表要素及其对区域

2、生产、生活和生态环境等活动的影响进行综合定量评估。同时使用景观格局分析方法，设计了包括斑块密度、多样性指数等面向自然地表要素空间格局特征表征的系列指数。克服了现有区域自然地表特征计算与表征技术在表征区域自然地表特征时缺乏综合评估模型、无法表征区域自然地表空间特征等不足。本发明提出的技术方案将服务于即将在全国范围内开展地理国情统计分析工作，同时也是面向地理国情时空数据挖掘的技术基础。51INTCL权利要求书5页说明书14页附图2页19中华人民共和国国家知识产权局12发明专利申请权利要求书5页说明书14页附图2页10申请公布号CN104318066ACN104318066A1/5页21一种自然地表

3、特征表征方法，其特征在于，包括以下步骤步骤1确定要进行自然地表特征评价与表征的空间范围和统计单元；步骤2从地理国情普查或监测数据集中读取计算范围内地理国情普查的10米空间分辨率数字高程模型DEM数据；其中，DEM以栅格像元的方式进行存储，每个栅格像元表示10米10米的一个空间位置，栅格像元的灰度值表示了该位置的高程海拔；遍历计算范围内的所有统计单元，依次计算各统计单元的地形指数TI；其中，地形指数用于综合刻画和反映统计单元内的自然地形条件和特征；步骤3遍历计算范围内的所有统计单元，依次计算各统计单元的地貌指数FI；步骤4遍历范围内的所有统计单元，依次计算各统计单元的地表覆盖指数CI；步骤5在步

4、骤2、步骤3、步骤4计算结果的基础上，遍历计算范围内的所有统计单元，计算所有统计单元的自然地表指数LI以综合评估区域自然地理国情综合条件对于生产、生活和生态环境保护的优势；步骤6遍历区域内的所有统计单元，从斑块格局、类型格局和区域总体格局三个不同层次，依次提取各统计单元的自然地理国情要素的空间格局分布特征。2根据权利要求1所述的自然地表特征表征方法，其特征在于步骤1中所述的空间范围是指县及县级以上的行政区或特定的流域、城市圈、城市群、经济带类在地域空间上连片的、封闭的一个区域；所述的统计单元包括标准网格单元和行政区划单元，是地理国情评价与表征的基本单位。3根据权利要求1所述的自然地表特征表征方

5、法，其特征在于步骤2中所述的地形指数TI的计算方法包括以下子步骤步骤21利用地理信息技术GIS，将各个统计单元和区域DEM空间叠置，获取落在各统计单元多边形内的DEM像元集；在此基础上统计出每个统计单元内栅格像元集的最大值EMAX，最小值EMIN，平均值EAVG；步骤22利用地理信息技术GIS，将各个统计单元和区域DEM叠置，计算各统计单元的平面投影面积AP及其在DEM上的表面面积AS；步骤23计算统计单元内的地形起伏度TA，宏观反映统计单元内的自然地形高低起伏程度，其中地形起伏度的计算公式为TAEMAXEMIN；步骤24计算统计单元内的地表粗糙度TR，宏观反映统计单元内的自然地表粗糙程度，其

6、中地表粗糙度的计算公式为TRAS/AP；步骤25计算统计单元内的地表切割深度TP，直观反映统计单元内地表被侵蚀切割的情况，并对这一地学现象进行了量化，其计算公式为TPEAVGEMIN；步骤26对区域内各个单元的EAVG、TA、TR和TP进行标准化，使各项指标的值的取值区间在1100之间；各项指标的标准化的方法分别为其中FU为统计单元U标准化后的指标值，TU为统计单元U标准化前的原始指标值，T为该项指标在统计范围内各单元的最大值；步骤27计算各统计单元内的地形指数TI，综合反映统计单元内的自然地形特征；对于任意统计单元U的地形指数TI计算公式为权利要求书CN104318066A2/5页3TIW1

7、FELEW2FTAW3FTRW4FTP；其中W1W4分别对应平均高程、地形起伏度、地表粗糙度和地表切割深度指标的权重，取值范围为0100；FELE为经过标准化后的统计单元U的平均高程指标值；FTA为经过标准化后的统计单元U的地形起伏度指标值；FTR为经过标准化后的统计单元U的地表粗糙度指标值；FTP为经过标准化后的统计单元U的地表切割深度指标值。4根据权利要求3所述的自然地表特征表征方法，其特征在于步骤27中所述的权重采用特尔斐法专家打分法或层次分析法AHP法进行确定。5根据权利要求3所述的自然地表特征表征方法，其特征在于步骤27中所述的权重W1取10，W2取30，W3取40，W4取20。6根

8、据权利要求1所述的自然地表特征表征方法，其特征在于步骤3中所述的地貌指数计算公式如下FI1F平原08F盆地06F丘陵F高原F山地150；其中，F平原、F盆地、F丘陵、F高原、F山地分别表示各种地貌在统计单元内所占的面积比例，取值为0100之间；上述面积比例值直接来自于地理国情普查工作的结果数据，计算方法为FIAI/A；其中，A指统计单元的面积，AI指该统计单元内某种地貌类型的面积，FI即为该地貌类型的面积比例。7根据权利要求1所述的自然地表特征表征方法，其特征在于步骤4中所述的地表覆盖指数计算公式如下CIF耕F园05F林F草F水F荒150；其中，F耕、F园、F林、F草、F水、F荒分别表示各种耕

9、地、园地、林地、草地、水域和荒地在统计单元内所占的面积比例，取值为0100之间；上述面积比例值直接来自于地理国情普查工作的结果数据，计算方法为FIAI/A；其中，A指统计单元的面积，AI指该统计单元内某种地表覆盖类型的面积，FI即为该地表覆盖类型的面积比例。8根据权利要求1所述的自然地表特征表征方法，其特征在于步骤5中所述的自然地表指数计算公式如下LIWTTIWFFIWCCI；其中WT、WF和WE分别对应地形指数、地貌指数和地表覆盖的权重，取值范围为0100。9根据权利要求8所述的自然地表特征表征方法，其特征在于步骤5中所述的权重采用特尔斐法专家打分法或层次分析法AHP法进行确定。10根据权利

10、要求8所述的自然地表特征表征方法，其特征在于步骤5中所述的权重WT取30，WF取30，WC取40。11根据权利要求1所述的自然地表特征表征方法，其特征在于步骤6中所述的区域总体格局表征包括地表斑块密度、地表边缘密度、地表多样性指数、地表优势度指数、地表均匀度指数和地表破碎度指数共6个定量指标；所述的地表斑块密度等于统计单元内地表覆盖图层中所有的斑块总数除以统计单元的总面积，反映统计单元内地理国情要素斑块在空间上的分布密度；其计算公式为式中PD表示统计单元内的地表斑块密度，M表示统计单元内的地表覆盖类型的个数，权利要求书CN104318066A3/5页4NI表示统计单元内第I种地表覆盖类型的斑块

11、个数在景观格局分析和空间分析中，斑块可以理解为图斑或地块，即用一个封闭多边形表示空间上连片的、具有相同地表覆盖类型的区域，对应地理国情普查数据集中地表覆盖要素图层中的一个图斑，A表示统计单元的地表总面积；所述的地表边缘密度反映统计单元内不同类型地理国情要素斑块间的边界长度，其计算公式为其中，PIJ表示景观中第I类地理国情要素斑块与相邻第J类地理国情要素斑块间的边界长度；所述的地表多样性指数是基于信息论基础上，用来度量一定范围内地理国情要素系统的空间结构组成复杂程度的指数，其计算公式为式中，H为多样性指数，M为地表覆盖类型总数；PK为地表覆盖类型面积比例，即某类地表覆盖类型斑块的总面积与统计单元

12、内所有地表覆盖类型斑块总面积的比值；当M1，说明一个地区只有一种地表覆盖类型，此时H0，多样性指数最小；如果某一个地区地表覆盖类型越多，则M越大；当M时，多样性指数H达到最大值；所述的地表优势度指数用于测定地表空间格局构成中一种或几种类型支配整个地表构成的程度，计算公式为式中，D为优势度指数；H是多样性指数，HMAX是其最大值，是给定区域内地表覆盖最大可能多样性指数，PK为地表覆盖类型面积比例，M为地表覆盖类型总数；所述的地表均匀度指数用以反映地表覆盖空间格局构成中各类型分布的均匀程度，均匀度可表达为所述的地表破碎度指数反映地表覆盖景观的破碎化程度，即地表覆盖被分割的破碎程度，计算公式为CN/

13、A；式中，N为区域地表覆盖图斑总数，A为整个区域土地总面积。12根据权利要求1所述的自然地表特征表征方法，其特征在于步骤6中所述的类型格局表征包括类型斑块密度、类型边缘密度、类型多样性指数、类型优势度指数、类型均匀度指数、类型破碎度指数共6个定量指标；所述的类型斑块密度等于类型斑块总数除以类型总面积，用于反映单元内某种地表覆盖类型的斑块分布状态，计算公式为权利要求书CN104318066A4/5页5式中PDI表示类型斑块密度，NI表示研究范围内类型斑块总数，AI表示研究范围某类型总面积；所述的类型边缘密度为统计单元内单位面积上不同类型景观要素斑块间的边缘长度，计算公式为式中PIJ表示景观中第I

14、类景观要素斑块与相邻第J类景观要素斑块间的边界长度；所述的类型多样性指数是基于信息论基础上，用来度量类型结构组成复杂程度的指数，计算公式为式中，H为多样性指数，M为某类型图斑总数，PK为斑块面积比例，即某类地表覆盖类型斑块面积与该类型斑块总面积的比值；所述的类型优势度指数用于测定某类型图斑构成中一个或几个图斑支配整个类型构成的程度，计算公式为式中，DT为类型T的优势度指数；N为类型T的图斑个数；HMAX是SHANNON多样性指数最大值，PI为某类型图斑I占该地类总面积的比例；所述的类型均匀度指数用以反映某类型构成中各图斑大小分布的均匀程度，均匀度可表达为所述的类型破碎度指数反映地表覆盖景观的破

15、碎化程度，即某类地表覆盖被分割的破碎程度，计算公式为CNI/AI；式中，NI为区域某类图斑总数，AI为某类地表总面积。13根据权利要求1所述的自然地表特征表征方法，其特征在于步骤6中所述的斑块格局表征包括斑块平均面积、最大和最小斑块面积、斑块面积标准差、斑块面积变动系数、斑块圆度、斑块矩度、斑块形状复杂度共7个定量指标；所述的斑块平均面积表示斑块平均大小，计算公式为式中，NI表示第I类斑块总数；AIJ表示第I类第I个斑块的面积；所述的最大和最小斑块面积表示统计单元内斑块面积的最大最小值；所述的斑块面积标准差表示统计单元内某类斑块面积的统计标准差；权利要求书CN104318066A5/5页6所述

16、的斑块面积变动系数表示统计单元内斑块规模的变异程度，计算公式为所述的斑块圆度表示斑块与圆的偏离程度，计算公式为所述的斑块矩度表示斑块与正方形的偏离程度，计算公式为所述的斑块形状复杂度采用斑块分维数表示斑块的形状复杂程度，计算公式为式中，P是斑块的周长，A是斑块的面积，D是分维数，K是常数。权利要求书CN104318066A1/14页7一种自然地表特征表征方法技术领域0001本发明属于地理国情统计分析与时空数据挖掘技术领域，特别是涉及一种自然地表特征表征方法。背景技术0002地理国情是地形地貌、地表覆盖、交通网络、城镇布局等多种自然和人文要素在宏观层面的综合表达，反映了国家自然资源、经济社会和人

17、文在地理空间上的时空特征与基本规律，是基本国情的重要组成部分。0003地理国情根据其表达的专题要素内容的不同，可分为自然地表、地理区位和生态覆被等多个专题内容以反映地理国情的不同方面。其中，自然地表主要反映区域地形、地貌、地表覆盖和地表格局等自然地理国情要素的时空格局及变化特征；地理区位则主要从区位的角度主要反映区域的自然资源环境、城镇布局和交通网络等区位特征的时空格局及变化特征；生态覆被则主要从生态环境的角度对区域的植被覆盖、水域湿地、生物丰度和土地退化等特征进行评价和表征。0004为全面掌握我国地理国情，满足社会经济和生态文明发展的决策需求，国务院提出了于2013年至2015年在全国范围内

18、开展第一次地理国情普查。地理国情信息的普查、监测和统计分析工作的开展对地理空间信息的采集、存储、可视化与统计分析提出了更高要求。遥感技术、地理信息技术、全球定位技术、并行计算技术和时空数据挖掘等技术被广泛用于地理国情要素信息的采集、存储和统计分析。通过开展地理国情普查和监测工作，能够准确、及时的掌握区域地理国情要素耕地、园地、公路、河流、湖泊等的空间位置、数量和形状等信息。地理国情普查与监测是地理国情评价、表征与应用的前提，为地理国情的表征和决策应用提供重要的数据基础。0005随着全国第一次地理国情普查试点工作的完成和后续工作的推进，如何基于地理国情普查和监测基础数据，结合必要的社会经济统计数

19、据，准确、客观、合理的表征区域地理国情，服务政府、部门和公众的多层次决策，已成为地理国情普查监测成果应用的迫切需求。0006基于我国地理普查监测工作的实践，地理国情评价与表征相关的技术主要可分为以下几种0007一基于点、线、面的地理国情矢量图形表征技术0008地理国情普查监测数据集中，地理国情相关的要素对象被抽象为点、线、面等几何形态进行表征。如，高速公路出入口、学校、医院等被抽象为带属性的点状要素；公路、铁路、河渠被抽象为带有属性的线状要素；湖泊、耕地、园地、城市建成区被抽象为带属性信息的面要素。以点、线、面表征地理国情要素能够非常详细、具体和直观的表达地理国情要素的空间位置、形状。其缺点主

20、要在于缺乏对区域地理国情信息的宏观整体把握，如无法直接表达区域内各类地理国情要素的面积、长度、个数、密度等信息。0009二基于汇总统计方法的地理国情评价与表征技术说明书CN104318066A2/14页80010汇总统计方法是以地理国情普查/监测数据为统计对象，针对不同的统计单元，根据地理国情实体的空间形态特征，将其划分为点状地理国情实体、线状地理国情实体和面状地理国情实体，以数量、密度、长度、面积、距离等作为统计指标，开展统计分析。以行政区划为单元，以地理国情要素为统计对象进行的地理国情汇总统计方法和步骤具体见国务院第一次全国地理国情普查领导小组办公室发布的系列技术规范。0011基于汇总统计

21、方法的地理国情表征技术一定程度能够从整体上反映区域地理国情要素的数量特征，如获得某个行政区范围内的耕地面积、园地面积、各等级道路的长度、学校设施的密度等。然而，地理国情是一种与空间位置相关的基本国情，以汇总统计的方法表征地理国情只能从数量上反映区域某一类地理国情要素的情况，无法实现对区域地理国情的综合评价，也不能反映区域地理国情的空间格局特征。0012自然地表特征是区域地形、地貌、地表覆盖和地表格局特征的综合体现，反映区域的基本自然地理客观条件和地表空间格局。自然地表特征是自然条件和长期人类活动共同作用的结果，也是未来人类社会经济发展物质基础和背景条件，是基本的地理国情。设计一种全面、客观、综

22、合的刻画区域自然地表特征、反映区域自然地理条件和地表空间格局的评价与表征技术，进而服务于社会经济发展决策，是当前地理国情普查与监测中亟待解决的关键技术问题。0013现有技术方案无法对区域自然地表特征和格局进行综合、全面的描述和表征。发明内容0014为了解决上述的技术问题，本发明面向区域自然地表特征综合表达与应用需求，提出一种自然地表特征表征方法，为地理国情普查监测的后续成果统计分析、数据挖掘与应用工作奠定技术基础，推进地理国情普查监测成果服务于区域发展决策水平。0015本发明所采用的技术方案是一种自然地表特征表征方法，其特征在于，包括以下步骤0016步骤1确定要进行自然地表特征评价与表征的空间

23、范围和统计单元；0017步骤2从地理国情普查或监测数据集中读取计算范围内地理国情普查的10米空间分辨率数字高程模型DEM数据；其中，DEM以栅格像元的方式进行存储，每个栅格像元表示10米10米的一个空间位置，栅格像元的灰度值表示了该位置的高程海拔；遍历计算范围内的所有统计单元，依次计算各统计单元的地形指数TI；其中，地形指数用于综合刻画和反映统计单元内的自然地形条件和特征；0018步骤3遍历计算范围内的所有统计单元，依次计算各统计单元的地貌指数FI；0019步骤4遍历范围内的所有统计单元，依次计算各统计单元的地表覆盖指数CI；0020步骤5在步骤2、步骤3、步骤4计算结果的基础上，遍历计算范围

24、内的所有统计单元，计算所有统计单元的自然地表指数LI以综合评估区域自然地理国情综合条件对于生产、生活和生态环境保护的优势；0021步骤6遍历区域内的所有统计单元，从斑块格局、类型格局和区域总体格局三个不同层次，依次提取各统计单元的自然地理国情要素的空间格局分布特征。0022作为优选，步骤1中所述的空间范围是指县及县级以上的行政区或特定的流域、城市圈、城市群、经济带类在地域空间上连片的、封闭的一个区域；所述的统计单元包括标说明书CN104318066A3/14页9准网格单元和行政区划单元，是地理国情评价与表征的基本单位。0023作为优选，步骤2中所述的地形指数TI的计算方法包括以下子步骤0024

25、步骤21利用地理信息技术GIS，将各个统计单元和区域DEM空间叠置，获取落在各统计单元多边形内的DEM像元集；在此基础上统计出每个统计单元内栅格像元集的最大值EMAX，最小值EMIN，平均值EAVG；0025步骤22利用地理信息技术GIS，将各个统计单元和区域DEM叠置，计算各统计单元的平面投影面积AP及其在DEM上的表面面积AS；0026步骤23计算统计单元内的地形起伏度TA，宏观反映统计单元内的自然地形高低起伏程度，其中地形起伏度的计算公式为TAEMAXEMIN；0027步骤24计算统计单元内的地表粗糙度TR，宏观反映统计单元内的自然地表粗糙程度，其中地表粗糙度的计算公式为TRAS/AP；

26、0028步骤25计算统计单元内的地表切割深度TP，直观反映统计单元内地表被侵蚀切割的情况，并对这一地学现象进行了量化，其计算公式为TPEAVGEMIN；0029步骤26对区域内各个单元的EAVG、TA、TR和TP进行标准化，使各项指标的值的取值区间在1100之间；各项指标的标准化的方法分别为00300031其中FU为统计单元U标准化后的指标值，TU为统计单元U标准化前的原始指标值，T为该项指标在统计范围内各单元的最大值；0032步骤27计算各统计单元内的地形指数TI，综合反映统计单元内的自然地形特征；对于任意统计单元U的地形指数TI计算公式为0033TIW1FELEW2FTAW3FTRW4TT

27、P；0034其中W1W4分别对应平均高程、地形起伏度、地表粗糙度和地表切割深度指标的权重，取值范围为0100；FELE为经过标准化后的统计单元U的平均高程指标值；FTA为经过标准化后的统计单元U的地形起伏度指标值；FTR为经过标准化后的统计单元U的地表粗糙度指标值；FTP为经过标准化后的统计单元U的地表切割深度指标值。0035作为优选，步骤27中所述的权重采用特尔斐法专家打分法或层次分析法AHP法进行确定。0036作为优选，步骤27中所述的权重W1取10，W2取30，W3取40，W4取20。0037作为优选，步骤3中所述的地貌指数计算公式如下0038FI1F平原08F盆地06F丘陵F高原F山地

28、150；0039其中，F平原、F盆地、F丘陵、F高原、F山地分别表示各种地貌在统计单元内所占的面积比例，取值为0100之间；上述面积比例值直接来自于地理国情普查工作的结果数据，计算方法为FIAI/A；其中，A指统计单元的面积，AI指该统计单元内某种地貌类型的面积，FI即为该地貌类型的面积比例。0040作为优选，步骤4中所述的地表覆盖指数计算公式如下0041CIF耕F园05F林F草F水F荒150；0042其中，F耕、F园、F林、F草、F水、F荒分别表示各种耕地、园地、林地、草地、水域和荒地在统计单元内所占的面积比例，取值为0100之间；上述面积比例值直接来自于地理国情说明书CN104318066

29、A4/14页10普查工作的结果数据，计算方法为FIAI/A；其中，A指统计单元的面积，AI指该统计单元内某种地表覆盖类型的面积，FI即为该地表覆盖类型的面积比例。0043作为优选，步骤5中所述的自然地表指数计算公式如下0044LIETTIWFFIWCCI；0045其中WT、WF和WC分别对应地形指数、地貌指数和地表覆盖的权重，取值范围为0100。0046作为优选，步骤5中所述的权重采用特尔斐法专家打分法或层次分析法AHP法进行确定。0047作为优选，步骤5中所述的权重WT取30，WF取30，WC取40。0048作为优选，步骤6中所述的区域总体格局表征包括地表斑块密度、地表边缘密度、地表多样性指

30、数、地表优势度指数、地表均匀度指数和地表破碎度指数共6个定量指标；0049所述的地表斑块密度等于统计单元内地表覆盖图层中所有的斑块总数除以统计单元的总面积，反映统计单元内地理国情要素斑块在空间上的分布密度；其计算公式为00500051式中PD表示统计单元内的地表斑块密度，M表示统计单元内的地表覆盖类型的个数，NI表示统计单元内第I种地表覆盖类型的斑块个数在景观格局分析和空间分析中，斑块可以理解为图斑或地块，即用一个封闭多边形表示空间上连片的、具有相同地表覆盖类型的区域，对应地理国情普查数据集中地表覆盖要素图层中的一个图斑，A表示统计单元的地表总面积；0052所述的地表边缘密度反映统计单元内不同

31、类型地理国情要素斑块间的边界长度，其计算公式为00530054其中，PIJ表示景观中第I类地理国情要素斑块与相邻第J类地理国情要素斑块间的边界长度；0055所述的地表多样性指数是基于信息论基础上，用来度量一定范围内地理国情要素系统的空间结构组成复杂程度的指数，其计算公式为00560057式中，H为多样性指数，M为地表覆盖类型总数；PK为地表覆盖类型面积比例，即某类地表覆盖类型斑块的总面积与统计单元内所有地表覆盖类型斑块总面积的比值；当M1，说明一个地区只有一种地表覆盖类型，此时H0，多样性指数最小；如果某一个地区地表覆盖类型越多，则M越大；当M时，多样性指数H达到最大值；0058所述的地表优势

32、度指数用于测定地表空间格局构成中一种或几种类型支配整个地表构成的程度，计算公式为0059说明书CN104318066A105/14页110060式中，D为优势度指数；H是多样性指数，HMAX是其最大值，是给定区域内地表覆盖最大可能多样性指数，PK为地表覆盖类型面积比例，M为地表覆盖类型总数；0061所述的地表均匀度指数用以反映地表覆盖空间格局构成中各类型分布的均匀程度，均匀度可表达为00620063所述的地表破碎度指数反映地表覆盖景观的破碎化程度，即地表覆盖被分割的破碎程度，计算公式为0064CN/A；0065式中，N为区域地表覆盖图斑总数，A为整个区域土地总面积。0066作为优选，步骤6中所

33、述的类型格局表征包括类型斑块密度、类型边缘密度、类型多样性指数、类型优势度指数、类型均匀度指数、类型破碎度指数共6个定量指标；0067所述的类型斑块密度等于类型斑块总数除以类型总面积，用于反映单元内某种地表覆盖类型的斑块分布状态，计算公式为00680069式中PDI表示类型斑块密度，NI表示研究范围内类型斑块总数，AI表示研究范围某类型总面积；0070所述的类型边缘密度为统计单元内单位面积上不同类型景观要素斑块间的边缘长度，计算公式为00710072式中PIJ表示景观中第I类景观要素斑块与相邻第J类景观要素斑块间的边界长度；0073所述的类型多样性指数是基于信息论基础上，用来度量类型结构组成复

34、杂程度的指数，计算公式为00740075式中，H为多样性指数，M为某类型图斑总数，PK为斑块面积比例，即某类地表覆盖类型斑块面积与该类型斑块总面积的比值；0076所述的类型优势度指数用于测定某类型图斑构成中一个或几个图斑支配整个类型构成的程度，计算公式为00770078式中，DT为类型T的优势度指数；N为类型T的图斑个数；HMAX是SHANNON多样性指数最大值，PI为某类型图斑I占该地类总面积的比例；0079所述的类型均匀度指数用以反映某类型构成中各图斑大小分布的均匀程度，均匀度可表达为说明书CN104318066A116/14页1200800081所述的类型破碎度指数反映地表覆盖景观的破碎

35、化程度，即某类地表覆盖被分割的破碎程度，计算公式为0082CNI/AI；0083式中，NI为区域某类图斑总数，AI为某类地表总面积。0084作为优选，步骤6中所述的斑块格局表征包括斑块平均面积、最大和最小斑块面积、斑块面积标准差、斑块面积变动系数、斑块圆度、斑块矩度、斑块形状复杂度共7个定量指标；0085所述的斑块平均面积表示斑块平均大小，计算公式为00860087式中，NI表示第I类斑块总数；AIJ表示第I类第I个斑块的面积；0088所述的最大和最小斑块面积表示统计单元内斑块面积的最大最小值；0089所述的斑块面积标准差表示统计单元内某类斑块面积的统计标准差；0090所述的斑块面积变动系数表

36、示统计单元内斑块规模的变异程度，计算公式为00910092所述的斑块圆度表示斑块与圆的偏离程度，计算公式为00930094所述的斑块矩度表示斑块与正方形的偏离程度，计算公式为00950096所述的斑块形状复杂度采用斑块分维数表示斑块的形状复杂程度，计算公式为00970098式中，P是斑块的周长，A是斑块的面积，D是分维数，K是常数。0099相对于已有的区域自然地表特征表征技术，本发明使用综合评估模型对区域地理国情要素中的自然地表要素及其对区域生产、生活和生态环境等活动的影响进行综合定量评估。同时使用景观格局分析方法，设计了包括斑块密度、多样性指数等面向自然地表要素空间格局特征表征的系列指数。克

37、服了现有区域自然地表特征计算与表征技术在表征区域自然地表特征时缺乏综合评估模型、无法表征区域自然地表空间特征等不足。本发明提出的技术方案将服务于即将在全国范围内开展地理国情统计分析工作，同时也是面向地理国情时空数据挖掘的技术基础。附图说明0100图1本发明实施例的流程图。0101图2本发明实施例的统计单元的高程信息提取原理图。说明书CN104318066A127/14页13具体实施方式0102为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。0103本发明的理论基础为01041多

38、因素综合评价模型0105多因素综合评价模型是通过系统、综合地的分析各类因素和因子对区域自然地表特征的影响和作用强度，进而基于定量的方法对评价区域的自然地表特征进行系统、综合的评价，反映区域自然地表对土地利用、经济发展、生态环境保护等人类生产生活活动的难易程度和优劣度。多因素综合评价模型的指导思想是从影响区域自然地表特征或人类活动难易程度的原因出发，采取由原因到结果的评价方法。采用多因素综合评价模型，综合系统的评价区域自然地表特征的关键步骤在于如何科学、合理的确定评价指标体系。0106一、自然地表特征评价指标体系设计原则0107基于多因素综合评价模型的基本原理，构建面向区域自然地表特征评价的自然

39、地表特征评价指标体系的基本原则如下01081以地理国情普查和监测数据为基础。区域自然地表特征评价是地理国情普查和监测工作的重要组成部分。地理国情普查和监测是区域自然地表特征评价的主要数据来源。因此，区域自然地表特征评价所设计的评价指标的数据来源也必须以地理国情普查和监测数据为基础，保证设计的指标可获取，具有可行性。01092具有时空动态变化特征。区域自然地表特征同时具有空间特征和时态特征，表现出时空动态特性。区域自然地表特征在年度或若干年度内具有可监测的变化，通过时序分析可反映区域自然地表特征的时空动态变化状况和变化趋势。01103综合分析和主导因素分析相结合。自然地表特征具有高度综合性。自然

40、地表特征的评价需要考虑多项指标进行综合计算。但是需要注意构成具体地理国情的指标可能不是同等重要的，在计算时需要考虑因素重要性进行选取或合理确定权重。0111二、自然地表特征评价指标体系设计方案0112基于上述原则，构建区域自然地表特征评价指标体系方案如表1；0113表1自然地表特征综合评价指标体系0114说明书CN104318066A138/14页1401152景观格局分析方法0116景观是在气候、地形、地貌、土壤、植被、水文等自然要素及人文干扰作用下形成的综合有机体。景观格局，一般指景观的空间格局，由一系列大小和形状各异，排列不同的景观斑块在景观空间的排列。景观格局分析是景观生态学的重要组成

41、部分，通过一系列景观格局指数反映区域景观要素在空间上的异质性。由此可见，景观格局分析方法是用于表征自然地表要素在空间上的分布特征和规律。0117一、景观格局分析的对象0118斑块。斑块是景观格局的基本组成单元，是指不同于周围背景的、相对均质区域。斑块反映了系统内部和系统间的相似性或相异性。斑块的大小、形状、边界以及斑块之间的距离等空间分布特征构成了不同的生态带，形成了生态系统的差异，对区域生态系统过程产生重要影响。地理国情普查与监测中，斑块对应地表覆盖要素图层中的某一个具体的多边形图形，如某一具体的耕地图斑、园地图斑等。0119类型。景观格局中要素类型相同的所有斑块的集合。例如，在地理国情普查

42、监测中所指的类型包括各类地表覆盖中的各类地理要素的集合，如耕地、园地、林地、草地、水体、裸地等要素。针对某一要素类型的统计即是针对某一种地理国情要素的所有斑块的统计。0120景观格局。区域内由大小、形状、类型各不相同的斑块构成了区域整体的景观格说明书CN104318066A149/14页15局。对于自然地表综合评价与表征问题，针对景观格局整体的计算就是对评价单元内所有斑块整体所呈现出的特征进行表征。0121二、景观格局分析的方法与指标0122斑块密度。反映一定范围内地理国情要素斑块在空间上的分布密度。计算公式为0123式中PD为统计单元内的地表斑块密度，M为统计单元内的地表覆盖类型个数，NI为

43、统计单元内的地表覆盖类型I的斑块个数在景观格局分析和空间分析中，斑块可以理解为图斑或地块，即用一个封闭多边形表示空间上连片的、具有相同地表覆盖类型的区域，A为统计单元内的地表总面积。0124边缘密度。边缘密度反映一定范围内单位面积上不同类型地理国情要素斑块间的边界长度。计算公式为0125式中PIJ为景观中第I类地理国情要素斑块与相邻第J类地理国情要素斑块间的边界长度。0126多样性指数。基于信息论基础上，用来度量一定范围内地理国情要素系统的空间结构组成复杂程度的指数。计算公式为0127式中，H为多样性指数，M为地表覆盖类型总数。PK为地表覆盖类型面积比例，即某类地表覆盖类型斑块的总面积与整个区

44、域所有地表覆盖类型斑块总面积的比值。当M1，说明一个地区只有一种地表覆盖类型，此时H0，多样性指数最小。如果某一个地区地表覆盖类型越多，则M越大。当M时，多样性指数H达到最大值。0128优势度指数。优势度指数用于测定地表空间格局构成中一种或几种类型支配整个地表构成的程度，计算公式为01290130式中，D为优势度指数；H是多样性指数，HMAX是其最大值，是给定区域内地表覆盖最大可能多样性指数，PK为地表覆盖类型面积比例，M为地表覆盖类型总数。0131均匀度指数。用以反映地表覆盖空间格局构成中各类型分布的均匀程度，均匀度可表达为0132破碎度指数。地表破碎度指数反映地表覆盖景观的破碎化程度，即地

45、表覆盖被分割的破碎程度。公式为CN/A；其中，N为区域地表覆盖图斑总数，A为整个区域土地总面积。0133基于自然地表特征的基本内涵和地理国情普查监测工作的数据特点，设计自然地表的景观格局指标体系如下表2；0134表2自然地表空间特征表征指标体系0135说明书CN104318066A1510/14页160136如上表所示，分别对统计单元从斑块、类型、景观整体三个不同尺度分别反映统计单元内自然地表的空间特征。0137基于上述理论基础，请见图1、图2，本发明所采用的技术方案是一种自然地表特征表征方法，具体包括以下步骤0138一、确定要进行自然地表特征评价与表征的空间范围和统计单元。其中，所指范围一般

46、是指县及县级以上的行政区或特定的流域、城市圈、城市群、经济带等在地域空间上连片的、封闭的一个区域；统计单元通常包括标准网格单元和行政区划单元等，是地理国情评价与表征的基本单位。计算范围和统计单元的类型、划分方法可参见技术规范地理国情普查基本统计技术规定GDPJ022013。设统计范围内共有N个统计单元。为了便于理解，本方案描述的范围可理解为县级行政区，统计单元可理解为乡镇级行政区，但具体实施时不限于以上内容，可按照地理国情普查基本统计技术规定GDPJ022013确定。0139二、从地理国情普查或监测数据集中读取计算范围内地理国情普查的10米空间说明书CN104318066A1611/14页17

47、分辨率数字高程模型DEM数据。其中，DEM以栅格像元的方式进行存储。每个栅格像元表示10米10米的一个空间位置，栅格像元的灰度值表示了该位置的高程海拔。遍历计算范围内的所有统计单元，依次计算各统计单元的地形指数TI。其中，地形指数用于综合刻画和反映统计单元内的自然地形条件和特征。请见图2，地形指数TI的计算方法包括以下子步骤0140利用地理信息技术GIS，将各个统计单元和区域DEM空间叠置，获取落在各统计单元多边形内的DEM像元集。在此基础上统计出每个统计单元内栅格像元集的最大值EMAX，最小值EMIN，平均值EAVG。0141利用地理信息技术GIS，将各个统计单元和区域DEM叠置，计算各统计

48、单元的平面投影面积AP及其在DEM上的表面面积AS。多边形的投影面积和表面面积计算方法可利用现有GIS软件实现，非本发明的关注重点，具体计算方法不做阐述。0142计算统计单元内的地形起伏度TA，宏观反映统计单元内的自然地形高低起伏程度。其中地形起伏度的计算公式为TAEMAXEMIN。0143计算统计单元内的地表粗糙度TR，宏观反映统计单元内的自然地表粗糙程度。其中地表粗糙度的计算公式为TRAS/AP。0144计算统计单元内的地表切割深度TP，直观反映统计单元内地表被侵蚀切割的情况，并对这一地学现象进行了量化，其计算公式为TPEAVGEMIN。0145对区域内各个单元的EAVG,TA,TR和TP

49、进行标准化，使各项指标的值的取值区间在1100之间。各项指标的标准化的方法分别为01460147其中FU为统计单元U标准化后的指标值，TU为统计单元U标准化前的原始指标值，T为该项指标在统计范围内各单元的最大值。0148计算各统计单元内的地形指数TI,综合反映统计单元内的自然地形特征。对于任意统计单元U的地形指数TI计算公式为0149TIW1FELEW2FTAW3FTRW4FTP0150其中W1W4分别对应平均高程、地形起伏度、地表粗糙度和地表切割深度指标的权重，取值范围为0100。FELE为经过标准化后的统计单元U的平均高程指标值；FTA为经过标准化后的统计单元U的地形起伏度指标值；FTR为经过标准化后的统计单元U的地表粗糙度指标值；FTP为经过标准化后的统计单元U的地表切割深度指标值；0151权重一般可采用特尔斐法专家打分法或层次分析法AHP法进行确定。也可采用本发明提供的参考权重值W1取10，W2取30，W3取40，W4取200152三、遍历计算范围内的所有统计单元，依次计算各统计单元的地貌指数FI。地貌指数综合反映了统计单元内的地貌类型条件对生产生活等人类活动的便利度。FI的值越高，则表明该单元内的地貌条件越有利于人类活动的开展，对该区域实施开发利用的成