Research on technology of forest subcompartments management plan assistant decision

系统的研建目标是实现一个森林小班经营方案编制辅助决策支持系统。经营者可根据经营任务，通过空间分析，从森林资源数据库中选择符合条件的小班，并保存为一个单独的图层；用户加载图层后，对每个小班进行立地质量评价、树种适宜性评价和收获预估，并作出经营类型推理，辅助用户编制每个小班的经营方案并输出。所有经营小班设计完成后，生成经营方案设计报告。

根据系统设计的目标，对小班经营方案编制的整个过程进行分解。根据经营任务，合理选择小班并保存为一个新的图层，用户加载经营小班图层以后，可对每个小班进行立地质量评价、生长收获预估、适宜性评价和经营方案辅助编制。系统的主要功能模块包括经营小班选择模块、立地质量评价模块、小班适宜性评价模块、生长收获预估模块和小班经营方案辅助编制模块。

系统采用C/S结构，基于. NET平台，使用C#语言，结合ArcGIS Engine（AE）9.3组件开发包进行组件式开发，实现可扩展、易维护的决策支持系统^[2]。系统自底向上分为数据层、服务层、人机交互层和用户。如图 1所示：数据层使用Personal Geodatabase空间数据模型，包括空间数据库、属性数据库、知识库、模型库等。服务层包括基本的信息服务、地图服务和小班经营方案编制辅助决策等。

Figure 1.  System structure diagram of decision support system (DSS)

森林资源数据库存储了所有小班的空间数据和属性数据，小班选择主要通过GIS空间分析实现。选择流程如图 2所示。

Figure 2.  Flow chart of subcompartments selection

对于空间定位选择，AE提供了ISpatialFilter接口，设置Geometry属性为空间二维目标几何体（圆、矩形或多边形等），设置WhereClause为基本属性过滤条件，基本属性条件通过XML文件进行读写。对于空间属性选择，则忽略Geometry属性，直接设置WhereClause属性。关键代码如下所示。

ISpatialFilter conditionFilter = new SpatialFilterClass（）; //定义ISpatialFilter

IGeometry myGeometry = this._mapControl.TrackCircle（）; //定义并获取几何体

conditionFilter.Geometry = myGeometry; //赋值空间属性条件

conditionFilter.WhereClause = CurLayerTypeOper.GetXMLConditions（）; //读XML文件并赋值

对于智能选择，采用面积或蓄积进行控制。先使用空间定位选择或点缓冲区分析方式选择出满足基本条件的小班，构成初始小班集合，然后计算这些符合条件的小班总面积或总蓄积，并以指定的目标总面积或总蓄积指标值和上限值进行控制；若累加的小班总面积或总蓄积大于指标值且小于上限值，则选择为目标小班，否则自动调整空间大小，重新选择生成初始小班集合，计算所有符合条件的小班总面积或总蓄积，再次与目标面积或蓄积进行比较，如此自动进行运算，直到选择出符合目标条件的小班集合为止。对于空间几何体或缓冲区调整算法如下：以面积控制为例，先计算出林场单位面积小班数和小班平均面积，再计算初始小班集合总面积与指标值的差值，然后通过差值与单位面积小班数量和小班平均面积间的换算即可确定空间几何体或缓冲区调整的大小。蓄积控制方法类似。

选择的小班集合保存为一个单独的经营小班图层，保存时使用IFeatureDataConverter2接口的ConvertFeatureClass方法进行；复制原小班的所有字段和属性值到新的图层，并增加类型为esriFieldTypeInteger的字段IsManageEdit，默认赋值属性值0，用于标识是否已对小班进行经营设计；新图层以经营任务年份+任务编号命名。

小班经营方案主要内容包括小班基本情况、地类、林种、立地类型、经营目标、经营措施、经营投资与效益分析、地点与责任单位、经营措施作业类型与作业时间、设计单位与设计人等内容^{[1, 13]}。小班经营方案辅助编制主要有2种方式：① 智能设计。根据小班林分状态、立地情况和经营目标等条件，利用专家知识进行推理，推荐多个合适的经营方案，用户根据经营目标和林分实际条件从推荐的方案中选择；然后由用户通过收获预估、适宜性评价等过程对选择的方案进行修改完善后直接生成为小班的经营方案。智能设计的核心是利用专家知识构建知识库，通过推理机进行方案推荐。② 手动辅助设计。用户通过用户界面（UI），从经营知识库中选择经营目标、经营技术措施等内容，系统利用经营方案模板自动根据用户的选择内容生成小班经营方案初稿；用户在小班经营方案初稿的基础上，根据立地质量评价、收获预估、适宜性评价等结果，修改完善小班经营方案初稿形成正式方案。小班经营方案推理流程如图 3所示。

Figure 3.  Flow chart of management plan inference rules

系统的知识库由规则表和事实表组成，采用关系数据库构建^{[2, 14]}。以森林经营专家多年经营实践为根据预先编制若干经营方案，利用森林经营方案的知识储备，根据小班的条件确定经营类型，得到根据该经营类型编制的经营方案。一种经营类型，可以有多个经营方案。本系统选择产生式规则表示专家进行经营方案推理规则。产生式规则表示法适合于表示因果之间的关系，经营者根据林分状态、立地条件和经营目标等条件判断而确定采用的经营类型的推理思维就属于因果关系。产生式由前件和后件构成，如“IF（华北落叶松∧中龄林∧郁闭度∈[0.5, 0.7] ∧缓坡∧坡下∧……）THEN华北落叶松大径材”，括号中内容为前件，“华北落叶松大径材”为后件。根据前文所述，同一后件可能具有不同的前件，故本研究集成了相同后件的所有前件，采用星号分隔存储不用的规则条件，有效地减少了规则条数^[13]。这种存储方式符合经营专家的知识构成，同时也方便对知识库的更新。经营类型推理规则表的关系模式采用关系表进行存储，经营类型推理规则的关系模式实例如表 1所示。

推理机按图 3给出的经营方案推理流程以及经营类型推理规则，使用正向推理机来搜索匹配经营作业知识^[13]。系统通过读取小班属性，获取小班当前状态，结合经营目标等条件，自动匹配知识库并给出经营类型推荐结果。用户还可以手动推理，与知识库中的知识进行匹配，最终给出经营类型推理结果。推理机工作流程如图 4所示。

Figure 4.  Working process of the inference engine

通过微软公司提供的Microsoft Word X Object Library操作Word，进而生成经营方案文档，其中X为版本号^[15]。为了操作方便，设计了专门的Word操作类WordOperate.cs，通过实例化该类，可以实现对Word的所有操作。对于手动辅助设计模式，用户可以通过设置的经营目标和选择的经营技术措施，从知识库中读取相关数据，自动写入文档，生成方案；其中方案标题由“立地类型+树种+林种”构成，如中山华北落叶松人工防风固沙林，然后通过代码控制写入到Word对应位置。以方案标题写入Word为例，主要代码如下。

using Word = Microsoft.Office.Interop.Word; //缩写引用

Word.Application wordApp = new Word.ApplicationClass（）; //初始化

Word.Document wordDoc = wordApp.Documents.Add（ref Nothing, ref Nothing, ref Nothing, ref Nothing）; //实例化文档

wordDoc.Paragraphs.Last.Range.Font.Name = "黑体"; //设置字体

wordDoc.Paragraphs.Last.Range.Font.Bold = 1; //加粗

wordDoc.Paragraphs.Last.Alignment = WdParagraphAlignment.wdAlignParagraphCenter; // //居中

wordDoc.Paragraphs.Last.Range.Text ="中山华北落叶松人工防风固沙林"; //写入标题

其他Word操作方式类似。完成后用户根据需要可适当修改立地质量评价、收获预估等结果并保存方案。对于智能设计模式，用户通过推理决策，得到适合当前小班、符合经营目标的一个或一组经营类型，结合立地质量评价、收获预估等过程，选择合适的方案并做适当修改，生成小班经营方案。小班经营方案设计完成以后，设置方案名称为“小班号+经营方案名.doc”，对应小班IsManageEdit字段属性值置为1，并用“√”号对系统左侧经营小班树列表进行已经营标识，同时经营小班图层中对应小班颜色修改为蓝色。修改小班颜色使用IUniqueValueRendererhe和ISymbol接口实现，根据IsManageEdit字段属性值，分别定义不同的Color对象并实例化ISymbol，最后使用IUniqueValueRendererhe对象的AddValue方法修改小班颜色。

系统作为森林经营基础平台的子系统之一，以内蒙古赤峰市桦木沟林场小班数据为对象进行了测试。林场总经营面积为68 667 hm²，有林地面积为45 113 hm²，活立木总蓄积为186.42万m³，灌木林面积为1 922 hm²。用材林面积为8 673 hm²，以华北落叶松Larix principis-rupprechtii为主。林场有小班2 612个，其中有林地小班2 251个；中龄林小班1 445个，幼龄林小班125个，近熟林小班143个，成熟林小班5个，过熟林小班63个，其他未标识龄组属性的有林地小班为470个。林场所有小班的属性数据和空间数据存储在森林资源数据库中。

经营任务指定面积为300 hm²。通过面积控制智能选择得到小班40个。保存为新的图层并加载以后，新图层叠加到底图上，默认填充绿色。用户通过对小班进行立地质量评价、适宜性评价和生长收获预估，进而进行经营方案辅助编制。系统主界面如图 5所示。

Figure 5.  Main interface of system

系统提供了2种编案方式。在手动辅助设计编案界面，用户根据系统读取到的小班信息和经营要求选择经营目标、经营技术措施等选项并提交，系统据此生成小班的基本经营方案；用户可再通过立地质量评价和生长收获预估等辅助功能进一步修改完善小班经营方案。智能设计编案方式则是在系统自动获取小班基本信息后，通过推理机推理得到多个合适的经营方案，用户选择其中之一并完善方案。生成小班经营方案功能如图 6所示。

Figure 6.  Example for management plan of subcompartments

小班是森林经营措施实施的最小单元，也是森林经营方案可操作性的最终体现。精细化经营方案要求利用现代信息技术以及地理信息系统技术，针对具体的森林类型和小班微环境定制林分作业法，制定全周期森林经营方案并实施。本研究以地理信息系统（GIS）、专家系统、决策支持系统等信息技术为基础，基于.NET平台，使用C#语言，结合AE 9.3组件开发包，研究设计了经营方案编制推理决策流程，设计实现了面向小班的经营方案编制的辅助决策支持系统，为林场森林经营管理人员提供了森林经营管理智能决策支持平台。系统根据经营任务选择符合经营要求的小班集合，通过对小班的立地质量评价、适宜性评价、生长收获预估等，结合经营方案辅助编制模块，智能辅助生成小班的经营方案。系统的研究与应用，填补了小班级经营方案编制辅助系统的空白，有效地提高了森林经营管理的信息化和智能化水平。