古树名木旅游最优路线设计及实现范文

《旅游》2018年第1期

摘要：为了提高公众对古树名木的研究、管理及保护政策的认知与响应度，更好地提高社会保护森林资源的意识，设计与实现综合评估价值最优的古树名木旅游路线。对研究区浙江省金华市各区县的人口密度、服务业发展状况、公路密度和古树名木区域分布等作相关性分析，为古树名木旅游线路设计提供可行性分析；通过层次分析法（AHP）对古树名木综合价值进行量化，并进行邻接点矩阵存储；利用Dijkstra算法遍历加权无向图，确定综合价值最优的古树名木节点集（区域内规划和手动预设2种形式）；古树名木综合价值最优节点集与百度地图JavascriptAPI进行集成，最终实现最优古树名木旅游路线的规划，从而吸引更多的社会资源投入到古树名木的相关研究中，在一定程度上解决了古树名木前期研究的可操作性不强、社会认知度不够等问题。

关键词：森林资源信息系统；古树名木；最优旅游线路；综合价值评价；邻接点矩阵

古树是指100年生以上的树木；名木是指具有社会影响力，闻名于世且树龄往往也是在100年生以上的树木［1］。古树名木具有科研价值、历史价值、社会价值、文化价值和景观价值等多元价值特征。对古树名木的前期研究主要在数据标准化采集与建档以及数据库建设，其中在数据库的建设方面国外文献相关研究较早，如PMRD（植物微数据库）［2］、MEDDB（药用植物数据库）［3］、NIASGBdb（尼亚斯资源库遗传资源和植物病害信息数据库）［4］的建设；国内研究起步相对较晚，但是发展很快，例如植物种子信息图像数据库［5］、药用植物基因数据库［6］、以及地方植物数据库［7］等的建设，实现了对古树名木的各种性状及基因特征永久性记录。古树馆的建设是从实践应用的层次上提高社会公众保护古树名木意识的成功尝试［8］。随着乡村生态游、民俗文化游和自驾游等旅游形式快速发展［9－11］，以古树名木为主题特色的森林生态旅游也逐渐成为一种旅游形式，促使古树名木价值研究逐步成为热点。CIAPALA等［12］对树木年轮进行分析，探究了影响游客评估的一个潜在方法；贾永生等［13］基于层次分析法（AHP）建立了古树名木的部分价值的评价指标体系；徐炜［14］通过选取古树名木中简单而且可操作性强、指标代表性广的9个因子，对古树名木的综合价值进行了定量分析；另外，对古树名木的景观价值［15］、等级价值［16］和损失额价值［17］等方面也有相关研究。古树名木大多立地于偏远地区、风景旅游胜地以及乡村，社会关注度较低；再加上近几年中国新农村建设不断推进、生态旅游快速发展、环境污染日益加剧等构成了古树名木管理与保护的诸多不利因素。本研究为了提高公众对古树名木的认知度以及社会参与度，以浙江省金华市古树名木为研究对象，设计综合价值最优的古树名木旅游路线———适应当今生态游、自驾游等旅游形式的发展潮流，使人们在游乐中感受古树名木的文化内涵。

1数据准备与可行性分析

1.1研究地区概况

金华市位于浙江省中部，28°32′~29°41′N，119°14′~120°46′30″E，东西跨度151km，南北跨度129km，土地面积为10941km2［18］。截至2014年底，全市公路总里程12269km（浙江省排名第4名），高速公路里程354km，一级公路580km，二级公路1201km，二级以上公路占公路总里程的17.41%［19］。

1.2古树名木信息概况

金华市9区县（兰溪市、永康市、义乌市和东阳市，浦江县、武义县和磐安县，金东区和婺城区。金东区和婺城区同归属于市区），截至2014年统计到的古树名木为10614株。古树名木隶属于51科98属181种，一级保护古树名木881株，二级保护古树名木2128株，三级保护古树名木7605株。其中武义县古树名木最多，为3016株，占据统计数据的28.41%（一级保护古树名木1.15%，二级保护古树名木4.36%，三级保护古树名木22.90%），占到各个级别的百分数为一级保护古树名木13.85%，二级保护古树名木21.76%，三级保护古树名木31.97%；浦江县最少，为575株，占据统计数据的5.42%（一级保护古树名木0.42%，二级保护古树名木1.10%，三级保护古树名木3.89%），占到各级别的百分比为一级保护古树名木5.11%，二级保护古树名木5.50%，三级保护古树名木5.43%。

1.3MATLAB程序最优旅游路线可行分析

Pearson相关系数是用来衡量定距变量间的线性关系，一般按r值的相关性强度进行划分。|r|＜0.4为相关不显著，0.4≤|r|≤0.7为显著相关，0.7＜|r|＜1.0为极显著相关［20］。设定金华市各区县人口密度变量x1，金华市各区县服务业发展状况变量x2和金华市各区县公路密度变量x3。分别为：x1＝［145.7047540；221.8769816；411.8295371；450.0978474；462.5646924；477.9234973；689.7044805；1116.7420810；497.0554839］；x2＝［37.8700000；82.3900000；223.5000000；910.6900000；352.7800000；70.7900000；474.5700000；1028.7300000］；x3＝［0.9700000；0.8700000；0.9100000；1.3400000；1.2800000；1.0400000；1.1000000；1.2900000］。数据经方差齐次性检验，利用Pearson相关性分析显示，金华市各区县人口密度与服务业发展状况极显著正相关（r=0.764，P＜0.01），金华市各区县人口密度与公路密度呈显著正相关（r=0.588，P＜0.05），金华市各区县服务业发展状况与公路密度呈极显著正相关（r=0.824，P＜0.01）。另外，金华市各区县古树名木分布变量设置为x4，x4＝［1866；3016；575；812；1830；912；863；740］，与金华市各区县服务业发展状况变量x2做相关性分析，显示金华市各区县古树名木分布与服务业发展状况呈显著负相关（r=-0.500，P＜0.05）。依据MATLAB程序运行结果分析，金华市各区县人口密度、服务发展状况和公路密度之间具有较强的正线性相关，各个要素之间具有一定平衡度，金华市各区县古树名木特征分布与金华市服务业发展状况呈比较明显的负线性相关关系，补缺了经济服务发展相对薄弱的区县。这样，古树名木旅游路线的设计可以引导人们选择偏远地区的民俗文化和乡村生态游产品，减小热点地区的旅游压力，从而达到双赢的目的。

2古树名木综合价值评价

2.1古树名木综合评估价值体系框架

本研究参照文献［13］和［14］从古树名木的生态价值（生物量价值如胸径、树高等，生物质量价值如生长势），科研价值（保护价值如保护级别，稀有价值依据中国稀有植物库和统计的株数为标准），景观价值（树木立地地段价值如立地小地名，观赏价值如树容树貌），历史价值（历史价值如历史典故）和社会价值（民俗文化如和当地民族文化是否相融合，文化信仰价值如民众供奉或者参拜祈福）等5个方面建立古树名木综合价值评价体系框架。

2.2MALAB程序一致性检验

在MATLAB窗体的CommandWindow中输入一致性检验程序，根据运行结果可得：方案层的初始加权值通过了一致性检验，同时得到各方案层的指标权重WPi。生物量价值P1的权重WP1＝0.0303，生物质量价值P2权重WP2＝0.0303，树龄调整价值P3的权重WP3＝0.0909，保护价值P4的权重WP4＝0.0606，稀有性价值P5的权重WP5＝0.1818，树木立地地段价值P6的权重WP6＝0.0303，观赏价值P7的权重WP7＝0.0909，历史史实价值P8的权重WP8＝0.2121，民俗文化价值P9的权重WP9＝0.1818，文化信仰价值P10的权重WP10=0.0909。依据方案层P1~P10的权重值和古树名木方案层的指标得分（TD1~TD10）。

2.3古树名木综合评估价指标调整价值设定

利用模糊综合评价（FSE法），对古树名木综合价值体系框架中的10个方案层（表2）进行等级量化。

3古树名木邻接点的存储

本研究采用邻接矩阵的方式存储古树名木节点信息，这样可以方便判断2个古树名木节点是否通路且很快计算出通路的综合评估价值总和，不需要开辟多余空间来存储古树名木之间的关系和通路信息。存储对象为：①一级保护古树名木；②二级保护古树名木和三级保护古树名木中的综合评估价值高于一级保护古树名木评估价值中值的部分。通过百度地图对选定的古树名木节点进行两两互联并用直线进行手绘划分，百度地图上标注的古树名木节点为随机选取的浦江县9株古树名木（Ti），地图缩放到乡镇/街道显示时，可以发现相离很近或同立地小地名的古树名木的标注重合在一起，如T4标注是位于浦江县的塔山公园内的3株古树名木的重合，T7标注是位于浦江县江南第一家的2株古树名木的重合。另外，对于不在同一个立地的小地名，且密度相对很高的古树名木，利用百度地图JavascriptAPI提供的点聚合思想，以古树名木为圆点，1000m为半径，添加圆形覆盖物。

4最优旅游路线的实现

4.1最优古树名木旅游路线的设计

起止点获取的方法：①金华市各区县内综合价值排名前2位的古树名木节点；②手动预选取的古树名木节点。利用Dijkstra算法遍历无向图（古树名木邻接点矩阵），得到综合价值最优的古树名木队列。具体实现为：①根据古树名木邻节点矩阵建立加权无向图。②初始化，即S＝{v}，v的权值为0。③设置U包含除v外的其他古树名木节点，即：U＝{其余古树名木节点}，若v与U中古树名木节点ui有边，则＜ui，v＞正常有权值；若u不是vi的出边古树名木邻接点，则＜ui，v＞权值为∞。④从U中选取一个权值v最小的顶点ui，把ui加入S中。⑤以ui为新中间点，更新U中各古树名木节点的权值；若从源点v到古树名木节点ui的权值（经过古树名木节点v）比原来权值（不经过古树名木节点v）短，则修改顶点ui的权值，修改后的权值的顶点v的权值加上边上的权值。⑥重复步骤④和⑤，直到所有古树名木节点都包含在S中。最优古树名木队列与百度地图API（JavascriptAPI）进行兼容性设计，通过百度地图渲染出古树名木旅游线路。

4.2代码实现

4.3实现效果

区域内古树名木最优旅游路线的实现，经过预提取区域内的古树名木综合评估价值最高的3个节点的提取，经过技术兼容性设计，调取百度地图javascriptAPI进行地图渲染出最优行车路线。

5结论与讨论

本研究以浙江省金华市古树名木数据为研究对象，分别对研究区各区县的人口密度、公路密度、服务业发展状况与古树名木区域分布特征进行相关性分析，论证了古树名木旅游最优路线设计的可行性；利用AHP法和FSE法实现了古树名木综合价值的评价，完成古树名木邻接点矩阵的存储；利用Dijkstra算法遍历古树名木的邻接点矩阵得到综合评价值最优的古树名木节点集，经过与百度地图JavascriptAPI进行兼容性设计，最终实现了古树名木最优旅游路线。实现的古树名木最优旅游路线，综合考虑了乡村生态游、民俗文化游和自驾体验游等因素，引导人们选择以古树名木为特色的旅游形式，在游乐中丰富古树名木知识。研究的最终目的是提高社会公众对古树名木的关注度，为古树名木相关研究及政策制定与实施提供一个良好的社会氛围；在此基础上，吸引更多的社会资源投入到古树名木的相关研究中去，一定程度上解决古树名木前期研究的可操作性不强、社会认知度不够等问题。

6参考文献

［1］中华人民共和国建设部.城市古树名木保护管理办法［Z/OL］.（2000-09-01）［2017-04-13］.

［5］宁静，杜中赫，刘雅慧，等.植物种子信息图像数据库的建立［J］.种子，2016，35（4）：126-129.

［6］朱英杰.药用植物基因资源的生物信息学研究［D］.北京：北京协和医学院，2014.

［7］盛耀，罗云波，谢子鑫，等.药用/工业用转基因植物数据库的构建［J］.食品安全质量检测学报，2014，5（8）：2544-2549.

［9］卢小丽，成宇行，王立伟.国内外乡村旅游研究热点：近20年文献回顾［J］.资源科学，2014，36（1）：200-205.

［10］翁钢民，李凌雁.中国旅游与文化产业融合发展的耦合协调度及空间相关分析［J］.经济地理，2016，36（1）：178-185.

［11］程哲，蔡建明，崔莉，等.乡村转型发展产业驱动机制：以盘锦乡村旅游为例［J］.农业现代化研究，2016，37（1）：143-150.

［13］贾永生，乌志颜，于海蛟，等.基于AHP的赤峰市古树名木综合价值评价［J］.内蒙古林业科技，2014，40（1）：43-46.

［14］徐炜.古树名木价值评估标准的探讨［J］.华南热带农业大学学报，2005，11（1）：66-69.

［15］孙超，车生泉.古树名木景观价值评价：程式专家法研究［J］.上海交通大学学报（农业科学版），2010，28（3）：209-217.

［16］杨韫嘉，王晓辉，乐也，等.古树名木价值等级的评估研究［J］.中国农学通报，2014，30（10）：28-34.

［17］米锋，李吉跃，张大红，等.北京地区林木损失额的价值计量研究；有关古树名木科学文化价值损失额计量方法的探讨［J］.北京林业大学学报，2006，28（增刊2）：141-148.

［18］金华市人民政府网.金华市服务业发展状况［EB/OL］.（2015-07-20）［2017-04-13］.

［20］贾平俊，何晓群，金勇俊.统计学［M］.4版.北京:中国人民大学出版社,2009:299-306.

作者：李记1,2；徐爱俊1,2 单位：1.浙江农林大学，2.浙江农林大学