林业科学研究Forest Research 2008,21(4):429~435 文章编号:1001—1498(2008)04-0429-07 基于GIS建模技术的县级林业区划方法研究 张 超 , 黄清麟 , 张晓红 , 涂年旺 (1.中国林业科学研究院资源信息研究所,国家林业局林业遥感与信息技术重点实验室,北京2.福建省永安市林业局,福建永安366000) 100091; 摘要:提出了基于GIS的县级林业区划方法框架;建立了永安市林业区划因素体系,并运用空间分析方法对各因子 进行了提取;探讨了基于GIS建模技术的县级林业区划方法。通过主成分分析、建立综合模型、判别分析、确定合理 分区个数和清除细小多边形等过程,最终得到的区划结果能够将各相对优势区域合理分开,分区林业生产状况统计 数据能较清晰地反映各分区的资源状况和相对优势,林业建设发展方向较为明确。 关键词:地理信息系统;林业区划;县级;主成分分析;判别分析 中图分类号:s757.4 文献标识码:A Study on Forestry Division Method at County Level Based on GIS with Integrated Modeling ZHANG Chao ,HUANG Qing—lin ,ZHANG Xiao—hong ,TU Nian—wang (1.Research Institute of Forest Resource Information Techniques,CAF;Key Laboratory of Forestry Remote Sensing and Information Technology,State Forestry Administration,Beijing 100091,China; 2.Yong'n Forestary Bumau,F ̄in Proviance,Yong'n 366000,F ̄ian,Chiana) Abstract:Forestry division has been discussed rarely since 1995 in China.The growing interest in forestry provides an opportunity to improve the research tools and approaches that can help achieve a broad undemtanding of forestry division.Based on the research in Yong’an,Fujian Province of China,a framework that may be useful in correlative foresty dirvision study was put forward:(1)The factor—system of Yong’an foresty dirvision was established;(2)The factors were extracted wih tthe spatial analysis;(3)The foresty dirvision method of integrated modeling at county level based on GIS was processed.The method was taken through Principal Components Analysis,building the integrated model,Discriminate Analysis,seeking the ideal number of divisions and eliminating the scrappy polygons.The outcome showed that the method based on GIS with integrated modeling can reflect the present situation of the foresty development.The relative dominantr subdivisions were distinguished reasonably,and a specific direction in foresty develropment was given.According to the statistics of the subdivisions’state,the relative dominant subdivisions were showed obviously. Key words:GIS;forest division;county level;PCA;Discriminate Analysis 林业区划作为林业生产和森林经营的基础,在 1995年以后研究较少。近年来,林业的地位和作用 我国2O世纪7O、8O年代曾经全面研究和开展过,自 得到越来越多人的重视,随着由传统林业向现代林 收稿日期:2007—10-25 基金项目:国家林业局“948”项目“森林经营单位水平森林可持续经营技术引进(2006-4-23)” 作者简介:张超(198O一),男,河北丰南人,中国林业科学研究院在读博士,主要从事森林可持续经营研究. 通讯作者:黄清麟(1967一),男,福建莆田人,研究员,博士,博士生导师,主要从事森林经理学研究,主要研究方向为森林可持续经营 理论与技术. 维普资讯 http://www.cqvip.com 430 林业科学研究 第21卷 业的重大转变,原来的林业区划已难以适应全面建 2 940 km ,大体上是九山半水半分田。属亚热带海 洋性气候,年平均气温14.3~19.2℃,年平均降水 量1 490~2 060 mm;地貌以山地、丘陵和山间盆地 设现代林业的要求;林业的经营管理实现手段发生 了较大变化,随着计算机技术、地理信息技术和航空 航天遥感等新技术在林业上的广泛应用,信息的获 取、处理和管理技术大大提高,为实现更加客观、科 为主;根据福建省林地立地等级划分结果,永安市属 I、Ⅱ级立地林地面积占全市林业用地面积 的63% 。 学和准确的县级林业区划提供了条件。 近年来,地理信息系统技术在林业专业性区划 中的应用取得了较大的进展,用于辅助指导林业生 永安市现有罗坊、上坪、槐南和青水4个乡,小 陶、安砂、贡川、西洋、曹远、大湖和洪田7个镇,燕 产,并已取得了较好成果。2002年,洪军等 以东 北林业大学帽儿山实验林场为例,以森林资源调查 资料为数据源,结合帽儿山实验林场的地形与自身 特点,在GIS支持下,按照森林的功能类型,将研究 区域的森林划分为商品经济林和生态公益林2个经 营区,在此基础上又划分了4个一级林种和16个二 级林种;2004年,曾思齐等 以武冈林场2003年森 林资源调查资料为数据源,结合武冈林场的地形与 区位特征,利用Arc/Info的空间分析与数据管理功 能,按照森林的功能类型及森林区划技术标准,将武 冈林场的森林资源划分为商品林和生态公益林2个 类型,在此基础上,二级区划分为4个林种,三级区 划分为9个林种。这些研究为探讨GIS技术在林业 区划中的应用奠定了基础。 传统的县级林业区划方法有很多,例如多元统计 分析中的谱系图聚类法、正交函数排序法,模糊数学 中的直接模糊聚类法、模糊相似优先比法和模糊综合 评判法,灰色系统中的灰色聚类法、灰色局势决策法 等等 。大多数指标只能依靠直接获取的数据,限制 了区划因子的选择范围 ;数据计算过程较为复杂, 工作量较大,区划结果不直观,不能快速反复地调整 分区个数和分区界线。而利用GIS的空间分析方法 可以提取到传统方法不能获取的因子信息 ,生成 各种专题图,如地形地貌分布图、林种分布图、优势树 种分布图和社会经济因子分布图等,能直观反映林地 空间分布信息;结合数学建模原理与方法 ,可直接 在图上快速进行县级林业区划工作,并能根据实际情 况和需要快速调整分区个数和分区界线,为林业规 划、管理和决策提供科学依据。本文以福建省永安市 为例,以村为区划基本单元,研究和探讨了基于GIS 建模技术的县级林业区划方法。 1 研究区概况 永安市位于福建省中部偏西,地处116。56 ~ 117。47 E,25。33 ~26。12 N,全市总面积 东、燕西、燕南和燕北4个街道办事处,228个行政 村 。根据永安市林业局2003年森林资源基础地 理信息数据统计:全市林业用地面积25.48万hm , 占总面积的86.26%。林业用地中有林地23.69万 hm ,疏林地0.06万hm ,灌木林地0.20万hm ,未 成林造林地0.65万hm ,无林地0.88万hm 。有林 地中用材林面积12.94万hm ,蓄积1 541.99万 m ;防护林面积4.60万hm ,蓄积548.74万m ;特 用林面积0.37万hm ,蓄积50.42万m 。 2 研究方法 2.1 基于GIS的县级林业区划方法框架 空间分析是对分析空间数据有关技术的统 称 “ 。空间分析方法包括各种几何或逻辑运算、 数理统计分析和代数运算等数学手段。GIS的空间 分析方法主要有空间查询与量算、叠置分析、缓冲区 分析、表面分析、邻域分析和空间统计分类分 析等 。 本文将GIS的空间分析技术应用于县级林业区 划方法中,提出了基于GIS的县级林业区划方法框 架,如图1所示。 首先,遵循县级林业区划的原则和依据,并在对 区划影响因素进行综合分析的基础上,分别从当地 自然条件、社会经济和林业状况等方面构建县级林 业区划的指标体系,运用GIS的空间分析技术对各 指标的数据进行提取,建立包含各指标数据的综合 数据库。然后,应用空问统计分类分析方法对综合 数据库中各指标的数据进行综合分类(分区),同时 辅以专家评判分级和决策支持,最终形成县级林业 区划结果和林业区划图,同时建立了县级林业区划 的基础数据资源库,在此基础上可进一步开展各种 专业性区划工作,并得到各类区划专题图。通过对 各指标数据的不断更新,能够及时得到调整后的县 级林业区划结果,为林业生产和森林经营提供科学 的依据。 维普资讯 http://www.cqvip.com 第4期 张超等:基于GIS建模技术的县级林业区划方法研究 431 t 划 土壤分布 一 千 影。 区 篙 区 地形地貌 _1查询与量算卜 r 划 原 l专家评判分级 主成分分析卜 甄琚贫猓库 I d I层次分析法l l聚类分析l 林业区划图 则 篓 气候分布 行政区划 l叠置分析l 交通网络 河流水系 lt 区划指标体系 t小班数据 l人口分布 缓冲区分析l D,- 因子 ■ l社会经济I 遥感影像 l表面分析l 析 一判别分析卜 各类专题图 DEM数据 l自然条件I 社会经济因子 l林业状况l _1邻域分析 专家辅助决策卜 图1 基于GIS的县级林业区划方法框架 2.2分区依据 标,建立综合模型;对综合值I,按其一维分布划定若 干区间,对各分区结果进行判别分析得到修正后的 预测分区结果,确定合理分区个数;对各区间按综合 值I,融合,生成初步分类结果,在此基础上清除细小 多边形,最终生成县级林业区划图和区划结果。主 要步骤如下: (1)以整个区域生态、经济和社会可持续发展 为指导原则 ; (2)考虑社会、市场发展的需求。全面了解和 研究社会、市场对林业发展和布局的要求; (3)考虑自然环境因子¨ 。自然条件,如地形 地貌、水分、温度和土壤等,是县级林业区划确定分 (1)分别提取各指标的数据,进行主成分分析; (2)根据各主成分的贡献率确定主成分个数m, 按照因子负荷量对各主成分命名; (3)确定各主成分方程 , ,…, ; 区林业发展方向的限制因素; (4)考虑生物生长规律,遵循适地适树的 原则 ; (5)考虑地区经济发展水平、林业所占地位、林 业科技及应用水平¨。。等。 2.3永安市林业区划因素体系 (4)建立综合模型 : + + 赢 一+.一 ( ) 县级林业区划因素体系的建立需要充分考虑区 域林业发展条件与优势、发展方向和水平等的相似 性和差异性,同时考虑林业发展的主要矛盾,并依据 区域发展水平进行 。本文通过综合分析永安 市自然条件、社会经济条件和林业生产条件对全市 式中:Y为综合值;A。,A:,…,A 为各主分量的特征 值; , ,…, 为各主分量; (5)对I,按其一维分布划分 个区间,对各分 林业区划的影响,遵循县级林业区划因素体系的建 立原则 ,主要选取平均海拔、年均降水、年均温度、 ≥1O℃活动积温、农民人均纯收入、人均林业产值、 区结果进行判别分析,得到各区间的正判率和修正 后的预测分类结果; (6)统计各分区的林业生产现状数据,并进行 极差标准化处理: ‰: 【L 一 人均林地面积、林业产值占工业、农业、林业总产值 比重、林地面积比重、森林覆盖率、天然林面积比重、 生态公益林面积比重和竹林面积比重13项指标作 为永安市林业区划因子。 2.4基于GIS建模技术的县级林业区划方法 (2) m n』 m (7)以各分区林业生产现状数据及极差标准化 处理后的数据为参照,分别判断各分区的林业发展 方向,将方向相同或相近的分区合并,讨论合理的分 区个数; (8)对以上得到的合理分区结果进行融合; 基于GIS建模技术 的县级林业区划方法的 基本思路:利用数理统计分析中的主成分分析法,把 原来的多个指标化为少数几个相互独立的综合指 维普资讯 http://www.cqvip.com 432 林业科学研究 第21卷 (9)清除细小多边形:统计细小多边形的林业 生产现状,与各分区林业现状平均数据对比,将其融 合到发展方向相同或相近的分区内; (10)生成县级林业区划图及区划结果。 2.5数据及处理方法 的比例。 根据式(3)计算各村单元的综合值y。在计算 y1、y2过程中,代人的各因子的数值应是观测值标 准化处理后的数值。 得到综合值l,后,即可对其分区划段,从而得到 初步分区结果。分几个区较为合适?怎样分区较为 科学合理?这是比较关键的。根据y的一维分布, 本文使用的主要数据有:1:25万比例尺永安市 行政区划图,lay格式;1:1万比例尺永安市森林资 源基础地理信息数据库,lay格式;永安市各乡镇多 年平均温度、年均降水量、≥10℃活动积温,纸质; 永安市2003年总人口数、农业人口数、农民人均纯 收人等数据,纸质。这些数据主要来源于永安市林 业局、气象局和统计局等单位。 本文使用的基础空间数据库中,数据图层、专题 (Theme)多,数据来源不同,其投影方式、坐标系和 数据格式等不完全一致。建立规范化的标准空间数 据库是进行空间分析的前提和必要条件,因此,首先 对数据进行检查和规范化处理。 (1)地图投影转换:利用ArcView GIS的View Properties将各数据图层进行投影转换。 (2)数据格式转换:主要在ArcMap中AreTool— Box的数据转换向导中完成。 (3)空间数据的编辑、裁剪和合并等:空间数据 的编辑主要在ArcView GIS中完成;数据裁剪和合并 在扩展模块中的Geoprocessing Wizard中完成。 3 结果与分析 3.1建立综合模型 利用GIS的查询、量算和统计等功能分别提取 13项指标的数据,并对以上提取的数据进行主成分 分析 ,结果如表1、2所示。从得到的相关阵特征 值矩阵表可知,前2个主成分的累计贡献率为 84.39%,包含了原始数据的主要信息,若将其余主 成分舍去只损失总信息量的15.61%,因此选取前 2个主成分构造综合模型。 由表1可知,第一主成分由多因素综合影响,其 中平均海拔和天然林比重2因子负荷量为最大;第 二主成分主要受经济因素影响,农民人均纯收人和 人均林业产值2因子的负荷量为最大。 根据表2分别得到各主成分方程,建立综合 模型: l,:0.700 6 Y1+0.299 4 y2 (3) 式中:Y为综合值;Y1、y2为第一、二主分量,其前面 的系数为各自特征值所占一、二主分量特征值总和 利用系统聚类方法,将l,分割为/"L个区间,对每个区 间分别以1,2,…,/"L进行标记。为使区划结果更加 合理,应取/"L:3,4,…,分别对各分类结果进行判 别分析,对修正后的预测分类结果统计各分区的林 业现状数据,判断其发展方向。对/"L进行讨论,最终 得到合理的分区个数。根据永安实际情况,本文依 次取/"L:3,4,5,6,7,利用SPSS软件对l,进行系统 聚类,得到各初步分类结果。 表1特征向量 表2主成分得分系数矩阵 3.2确定合理分区个数 以得出的初步分类结果作为分组变量,以13项 指标数据标准化后的数值作为判别依据变量,采用 Fisher判别系数,对初步分类结果进行判别分析,得 到修正后的预测分类结果 。限于篇幅,下面以 :3为例: 维普资讯 http://www.cqvip.com 第4期 张超等:基于GIS建模技术的县级林业区划方法研究433 分区—一 用材 —— 经济_蚕 。 防护。j亘丛 善 特用 竹 —T— —用材 经济—_蚕 — 堕 防护笔善 — —_ 特用 竹 ……。—一发展方向 注:加下划线部分为相对优势较为明显的发展方向 图2基于GIS综合模型的林业区划图(n=3) 第1分区中,用材林和经济林相对优势较为明 显,分区发展方向确定为“用材林+经济林”;第2分 区中,用材林和竹林相对优势较为明显,分区发展方 相对优势均较明显,各分区发展方向明确;凡=6,7 时,随着凡逐渐增大,分区结果趋向于破碎化,出现 了较多相对优势相同或相近的重复分区组合,因此, 不再赘述凡>7的情况,本文选取凡=5,即认为将研 究区划分为5个分区较为合理。 3.3清除细小多边形 向确定为“用材林+竹林”;第3分区中,防护林、特 用林和竹林相对优势较为明显,分区发展方向确定 为“防护林+特用林+竹林”。 从表4和图2中可以看出,各分区多为若干个 相对优势区域的组合,分区结果较为笼统,区划界线 粗糙,在此基础上可进一步细分,因此凡=3不合适。 以上通过主成分分析和判别分析得到的分类结 果中,由于空间异质性,必然会存在一些与周围区域 相互独立的图斑,表现为较大区域中存在的较细碎、 独立的若干块图斑。 按同样方法依次取凡=4,5,6,7,得到各分区结 果。对得到的结果对比分析:凡=3,4时,分区界线 较为粗糙,局部区域需进一步细分;凡=5耐,各分区 县级林业区划的任务是确定研究区总体范围内的 林业生产战略布局及分区林业事业的发展方向,属于 维普资讯 http://www.cqvip.com 434 林业科学研究 第21卷 宏观决策,因此要将细碎的图斑融合到与其相邻的较 大区域内。判定细小多边形的标准应根据图斑在总体 以确定图斑的融合方向(具体过程略)。 3.4区划结果 范围内的实际分布以及研究区的面积大小、区划基本 单元大小等具体情况确定。综合考虑以上因素,在本 文中选取细小多边形的原则为:明显独立于周围较大 按照最终调整完成的分类结果进行融合,最终 生成区划专题图和区划结果,打印输出。区划图如 图3所示,区划结果的正判率情况如表5所示。 研究结果表明,应用基于GIS建模技术的区划 方法,以村为区划基本单元,得到的区划结果中,各 分区能够较为明确地反映分区林业发展方向,通过 区域的一块或两块空间上连续分布的图斑(村单元)。 要融合这些图斑,即将它们与周围区域合并,是 应有一定原则的,不能随意合并。本文的方法:将这 些图斑的森林资源状况与每个分区的森林资源平均 状况相比较,以及查看其在区划图中的具体位置,据 与永安现状统计情况对比,区划结果在局部表现为 对现状布局的调整。 图3基于GIS综合模型的永安市林业区划图 表5分区融合结果与l3项指标的符合情况分析 够将具有相对优势的区域合理分开,林业建设发展 方向较为明确。 (2)通过将综合值l,按其一维分布进行聚类划 分若干区间,分别讨论不同分区个数与现状的实际 符合情况,参照各分区各林种面积百分比及极差标 准化后的数据,比较、判断各分区的林业发展方向。 得到的较理想的区划结果表明分区个数的确定方法 在本文中区划方法的适用性。 (3)由于县级林业区划本身的复杂性,未能形 成较为科学的县级林业区划因素体系,本文着重探 讨了基于GIS建模技术的县级林业区划方法,对县 级林业区划的因素体系未做深入研究。如何科学确 定一个具体区域林业区划的因素体系有待进一步 研究。 4结论与讨论 (1)通过与永安市实际情况对比,参照区划结 果的正判率,基于GIS建模技术的县级林业区划方 法对于各分区划分得较为准确,各分区的融合结果 与选定的13项区划指标的符合情况均超过80%,能 (4)由于时间和数据的限制,本文未能探讨表 面分析、邻域分析等其它空间分析功能在县级林业 区划中的应用。如何充分利用GIS提取相关县级林 维普资讯 http://www.cqvip.com 第4期 张超等:基于GIS建模技术的县级林业区划方法研究 435 业区划指标的功能有待进一步研究和完善。 参考文献: [1]洪军,蔡体久.基于GIS的森林分类经营区划[J].东北林业 大学学报,2002,30(4):14一l8 [2]曾思齐,肖化顺,石军南.地理信息系统支持下的森林分类区划 界定[J].中南林学院学报,2004,24(5):ll—l5 [3]刘建国,袁嘉祖.林业区划原理与方法[M].北京:中国林业出 版社,1994:3—107 [4]张超,黄清麟.林业区划研究综述[J].林业资源管理,2005 (5):l6—20 [5]Birkin M,Clarke G,Clrake M,et a1.Intelligent GIS,Location De— cisions and Strategic Planning[M].Cambridge:Geoinformation In— ternationa1.1996 [6]赵宪文,张淑娟,陆显祥.遥感资料和聚类方法在县级林业区划 中的应用[J].林业科学研究,1992,5(2):239—242 [7]张超.基于GIS的县级林业区划方法研究——以福建省永安 市为例[D].北京:中国林业科学研究院,2006:69—70 [8]永安市林业区划办公室.福建省永安市林业区划文件汇编[c]. 福建永安,1984 [9]永安市统计局.永安统计年鉴[c].福建永安,2004 [10]邬伦,刘瑜,张晶,等.地理信息系统原理、方法和应用 [M].北京:科学出版社,2001:161—194 [11]张超,陈丙咸.地理信息系统[M].北京:高等教育出版社, 1995:10—45 [12]边馥苓.地理信息系统原理和方法[M].北京:测绘出版社, 1996:87—93 [13]黄清麟.浅谈德国的“近自然森林经营”[J].世界林业研究, 2005,18(3):73—77 [14]谢守鑫,姜姜,余树全.区域森林资源分类区划方法的研究 [J].林业资源管理,2006(4):6一ll [15]陈建新,王明怀,殷祚云,等.广东省秃杉引种栽培效果及栽培 区划分研究[J].林业科学研究,2002,15(4):399—405 [16]黄清麟.森林可持续经营综述 J].福建林学院学报,1999,19 (3):282—285 [17]Ian Fry.Twists and turns in thejungIe:Exploring the evolution of lnad use.1nad—use change and forestry decisions within the Kyoto Protocol[J].Review of European Community and International Environmentla Law,2002,ll(2):159—168 [18]Haber W. Using lnadscape ecology in plnaning and management [M]//Zonneveld I S,Forman R T T.Changing Landscapes:An Ecological Perspective.New Y0rk:Springer—Verl ̄.1990: 2l7—23l [19]FAO.FESLM.An Imemational Framework for Evaluating Sustain— able land management[R].Rome:World Soil Resources Report 73,1993 [20]陈述彭,鲁学军,周成虎.地理信息系统导论[M].北京:科学 出版社,2000:33—58 [21]朱秀珍.大型灌区运行状况综合评价研究[D].武汉:武汉大 学,2005:89—90 [22]白亚恒.洪湖市可持续发展评价及对策研究[D].北京:中国 科学院测量与地球物理研究所,2001:ll—l4 [23]李建兵.基于GIS的鞍山市投资环境评价系统的研究[D].鞍 山:鞍山科技大学,2005:58—59 [24]廖伟华.基于GIS甘肃省可持续发展及其空间差异研究[D]. 兰州:兰州大学,2003:17—20 [25]唐守正.多元统计分析方法[M].北京:中国林业出版社, 1984:20—22 [26]卢纹岱.SPSS for Windows统计分析[M].北京:电子工业出版 社,2000:379—404
因篇幅问题不能全部显示,请点此查看更多更全内容