-
毛竹Phyllostachys edulis又名楠竹，它不仅具有较高的经济价值和社会价值，而且发挥着十分重要的生态功能[1]。科学经营管理毛竹林不仅能有效地保护生态环境，而且对中国毛竹生产应用等方面经济效益的发挥具有十分重要的现实意义。因此，如何提高毛竹林经营管理的技术水平一直是关注的焦点问题。可视化是利用图像处理技术或是计算中产生的数字信息变为直接的计算机图形，将已有的数据呈现在研究者面前，是研究数据表示、数据处理、决策分析等一系列问题的综合技术[2]。计算机可视化技术的发展为森林可视化管理提供了可能。森林可视化主要包括林木可视化和林地可视化。林木可视化技术按照建模方法的不同，可分为基于规则几何体建模[3]、交互式参数化建模[4-6]和基于图像纹理的建模[7-8]。基于规则几何体建模和基于图像纹理的建模具有建模方法简单，内存占用小的特点，但是林木表达比较粗略，真实感较差，近视点效果不够理想。交互式参数化建模是通过设定相关参数完成树木造型，能够快速地创建出逼真的树木模型，是林木可视化的实用方法，但是在大场景的可视化上，容易造成数据冗余，影响系统运行速度[9]。因此，有学者研究运用细节层次技术（levels of detail，LOD）提高森林场景的渲染效率，取得了不错的可视化效果[10-12]。林地可视化是地面起伏状态的可视化描述。数字高程模型（digital elevation model，DEM）是地形三维可视化中最常用的技术[13]。DEM有许多种表示方法，其中不规则三角网（triangulated irregular network，简称TIN）模型与其他DEM模型相比有着不可替代的优点。TIN模型可以根据地形的具体特点和复杂程度而确定采样点的密度和位置，这样可以避免数据冗余，并且可以相对准确地表达地形的结构和细部分析[14]。另外，还可以通过叠加图像纹理、分形纹理和地表地物的方式，提高地形的真实性[15]。森林可视化的方法有很多种，每种方法都有各自的优缺点，并且每一种方法都是针对具体情况而言。目前，毛竹林分可视化方面的研究还很少见报道。本研究以浙江省天目山国家级自然保护区的毛竹林为研究对象，在林分尺度上，探讨集成应用美国Onyx Computing公司的OnyxTREE BAMBOO软件[16]、Autodesk公司的3ds Max建模软件[17]和美国环境系统研究所的ArcGIS[18-19]等多种软件相结合实现毛竹林三维可视化的技术方法。采用专业竹子建模软件OnyxTREE BAMBOO进行毛竹基本模型的创建，具有逼真的视觉效果。以此为基础，在3ds Max中创建出活竹、死竹、断竹、裂竹、倒竹等多种状态的毛竹模型。最后，在ArcGIS中，采用TIN模型上叠加林地地表贴图的方法，对毛竹林地进行可视化表达，再叠加基于毛竹模型的符号化点要素图层，实现毛竹林可视化。本研究提出集成应用多种软件实现毛竹林可视化的技术方法，操作简单，用户容易掌握，可视化效果逼真，不需要复杂而费时的编程，为毛竹林可视化经营管理提供了新的技术方法。
-
天目山国家级自然保护区位于浙江省西北部临安市境内，地处30°18′30″～30°24′55″N，119°23′47″～119°28′27″ E，总面积为1 050 hm2，主峰仙人顶海拔1 506 m。年平均气温为8.8~14.8 ℃，≥10 ℃年积温2 500~5 100 ℃；年降水量为1 390~1 870 mm，相对湿度76%~81%。保护区地处中亚热带北缘向北亚热带过渡的地带，受亚热带季风性湿润气候影响，森林植被十分茂盛。受海拔的影响，山体植被分布有明显的垂直界限，海拔870 m以下为常绿阔叶林，870~1 100 m为常绿、落叶阔叶林，1 100~1 380 m为落叶阔叶林，1 380~1 500 m为落叶矮林。近自然毛竹林作为一种特殊的森林植被类型多分布在海拔350~900 m，且由于林立地条件和林内环境的相互作用，林下植被生物多样性较低[20-22]。
-
2009年7月，在浙江省天目山国家级自然保护区内，用全站仪建立了1块面积为100 m × 100 m的近自然毛竹林固定标准地，标准地中心海拔840 m，主坡向南偏东30°。利用相邻网格法[23]，将标准地划分100个10 m × 10 m单元网格。用全站仪对单元网格内每株毛竹进行精确定位，测定毛竹基部的三维坐标（x, y, z）。同时，调查每株毛竹的胸径、竹高、枝下高、冠幅等因子，以及死活、受损等生长状态因子（表 1）。此后，每年对毛竹固定标准地进行复查。本研究采用2014年的复查数据。
表 1 毛竹生长状态分类
Table 1. Classification of moso bamboo growth state
存活状态 受损状态 竹秆形态 状态名称 存活 完整 直立 活完直 弯曲 活完曲 倒伏 活损倒 受损 破裂 活损裂 折断 活损折 多种损伤 活损多 死亡 完整 直立 死完直 弯曲 死完曲 倒伏 死损倒 受损 破裂 死损裂 折断 死损折 多种损伤 死损多 -
为了满足毛竹单株建模的需要，以样地调查数据中毛竹的平均胸径为参考，选取胸径大小适中的5株毛竹，伐倒后详细测量各项参数，包括竹秆每竹节的长度、宽度，每个竹节点的扩展、宽度、枝下高、角度、簇级数等。选取最长一级分枝和二级分枝，按与竹秆相同的参数进行测量。此外，还测量竹叶的长度、宽度、底长和梢长等参数。
-
OnyxTREE BAMBOO是由Onyx Computing公司开发的一款专门用来创建竹子及竹林模型的建模软件，它可以创建单株竹子模型及在指定的区域内创建由多株竹子组成的竹林。通过调整各项参数的值，它可以创建不同种类竹子模型，还可以创建处于不同生长阶段的竹子模型[16]。
3ds Max是由美国Autodesk公司开发的一种基于矢量的三维造型和动画设计软件，它主要是运用计算机图形生成和图形处理技术模仿真实现象，建立具体的模型，可以制作出非常逼真的三维实体及动画[17]。
ArcGIS是美国环境系统研究所开发的软件，是世界上最广泛的地理信息系统（GIS）软件之一。ArcScene是ArcGIS三维分析模块3D Analyst的核心，它具有管理3D GIS数据，进行3D分析，编辑3D要素，创建3D图层以及把二维数据生成3D要素等功能，同时提供强大的三维可视化和空间分析能力[18-19]。
-
毛竹林的可视化包括毛竹单株建模和林地可视化2个环节。首先，通过OnyxTREE BAMBOO和3ds Max软件创建不同生长状态的毛竹模型。然后在ArcScene中与林地DEM模型进行叠加，实现毛竹林分可视化。
-
毛竹有多种不同生长状态，如死活、弯直、折断、开裂、倒伏（表 1）。为了更好地反映毛竹林的真实情况，需要对这些不同状态的毛竹分别建模。OnyxTREE BAMBOO官网上提供一些免费的竹子模型，可作为毛竹建模的基础模型。
毛竹建模的过程分3步。首先，下载与需要建模的毛竹类似的竹子模型。其次，将下载的模型加载到OnyxTREE BAMBOO中，根据毛竹单株参数调查数据，对模型进行修改、调整，以达到所需要的效果。最后，将模型保存为3DS格式文件，以便在ArcGIS中调用。
根据毛竹实际参数建立的模型，加载到ArcGIS中，仅显示枝干，并不显示竹叶。因为按毛竹实际大小建立的模型，产生大量的枝、叶，从而导致单株毛竹的数据量很大，不能正常显示毛竹的叶子。鉴于以上原因，在建模时将模型参数按比例调整为原来的50%，以达到良好的显示效果（图 1）。
图 1 毛竹单株建模
Figure 1. Single modeling of moso bamboo
目前，OnyxTREE BAMBOO软件只能创建出直立、弯曲和倒伏状态的竹子模型，无法创建折断、开裂的竹子模型。需要将它与3ds Max 软件相结合以解决此问题。方法是将OnyxTREE BAMBOO中创建的毛竹基本模型导入到3ds Max，在修改面板下，针对顶点、边、面等不同的子对象层级，综合运用移动、旋转、切片等工具，对模型进行修改，得到折断、开裂的毛竹模型。对于死竹的建模，只需在活竹模型基础上，去掉竹叶，调整竹秆、竹枝的颜色即可。最终创建出各种状态的毛竹模型，建立毛竹模型库（图 2）。
图 2 不同生长状态的毛竹模型
Figure 2. Moso bamboo models with different growth state
-
毛竹林地可视化的实现主要分为2个部分：一是数字高程模型DEM的建立；二是林地地表贴图的制作，并与DEM进行叠加。
-
将包含毛竹各种调查信息的Excel数据文件加载到ArcScene，根据表格中毛竹的坐标生成点要素（图 3）。图 3中每个点要素代表 1株毛竹。
图 3 毛竹三维点图
Figure 3. 3D points of moso bamboo
-
直接由毛竹点要素进行TIN的创建，得到的TIN表面会太过尖锐、粗糙，显示效果不理想。实践表明：先对点要素进行自然邻域插值处理（图 4），对得到的结果进行生成等值线的操作（图 5），再用生成的等值线创建TIN。这样创建出的TIN表面较为光滑，比较符合真实的地面起伏情况。
图 4 插值结果
Figure 4. Interpolation result
图 5 生成等值线
Figure 5. Create contour
-
为了提高林地地表的真实性，在标准地中选择有代表性的地段，拍摄标准地地表纹理照片，在Photoshop中对地表纹理照片进行处理，合成一张整体的标准地地表图片。在ArcMap中对其重采样，提高图像分辨率。再对其进行地理配准，使其空间范围与TIN保持一致。最后，加载到ArcScene中与TIN进行叠加。显示效果如图 6所示。
图 6 叠加地表贴图后效果
Figure 6. Superposition of surface map
-
在上述工作基础上，将已生成的毛竹点要素按毛竹的不同状态进行分类，生成各自对应的点要素图层，用相应状态的毛竹模型替换点要素符号，并进行相应调整，以实现整个林分的三维可视化。
以“活完曲”毛竹的可视化为例。具体实现过程如下：激活毛竹点要素图层，利用查询工具筛选出“活完曲”状态的毛竹点，生成新的图层。打开图层属性对话框，在符号系统下点击点符号，进入符号选择器，点击“编辑符号”按钮，在符号属性编辑器对话框下，导入“活完曲”模型替代原来点符号（图 7）。毛竹模型的大小，弯曲或倒伏的方向，可以通过点击“高级”按钮，对毛竹符号进行“大小”和“旋转”的调整来实现。这里选择毛竹属性中的“竹高”和“方位角”作为调整因子。其他状态毛竹可视化过程与上述相同，全部加载完毕，最终实现整个林分的可视化。
图 7 导入毛竹模型
Figure 7. Import moso bamboo model
从毛竹林分整体进行观察，由于死倒竹数量较少，显示不明显，林分基本呈现一片绿色（图 8）。如果将视野拉近，则进入竹林内部，可以很清楚地观察到毛竹的枝叶和不同生长状态的毛竹如活竹、死竹、倒竹、断竹、裂竹等，以及地表的起伏状况，有一种身临其境的感觉（图 9）。借助ArcScene的漫游功能，可以实现竹林内部实时漫游观察。通过设定漫游路径，可以实现固定路径漫游。单击任一株毛竹，还可以很方便地查看其相关信息，包括毛竹的竹高、胸径、冠幅、枝下高、坐标等信息（图 10）。
图 8 毛竹林分整体三维图
Figure 8. 3D map of whole moso bamboo stand
图 9 毛竹林内部情况
Figure 9. Internal situation of moso bamboo stand
图 10 查看毛竹信息
Figure 10. View information of moso bamboo
利用ArcGIS强大的空间统计功能，在属性表中可以很方便地根据毛竹地不同属性信息进行汇总统计，不仅可以很方便地统计出整个林分毛竹相关属性的极值、平均值、标准差等，也可以根据毛竹年份、所在的网格号、存活状态分别进行汇总统计，并且根据需要创建相应的图表（图 11，图 12）。
图 11 毛竹径阶分布
Figure 11. Distribution of moso bamboo diameter grade
图 12 毛竹林新竹株数
Figure 12. Number of new bamboo in moso bamboo stand
-
本研究集成应用OnyxTREE BAMBOO，3ds Max建模软件和ArcGIS中的ArcScene软件，主要将前者的三维建模功能和后者的三维显示功能相结合，实现毛竹林分的三维可视化，取得较好的可视化效果。观察者既可以俯视整片竹林，掌握竹林的整体生长情况，又可以近距离观察毛竹，精确了解任一毛竹相关信息。
研究探索出一套完整的毛竹林可视化技术方法，此方法可操作性强，将它应用于毛竹林可视化经营管理，不仅可以真实反映毛竹林分的当前生长状况，也可用于直观展现竹林采伐或自然干扰后的林分状况，为进一步研究毛竹林动态变化的可视化奠定了基础。
在毛竹林可视化的过程中还存在一些有待进一步研究的问题。交互式参数化建模所建立的模型精度越高，内存占有量越大，特别是在大场景森林可视化中，如何平衡模型精度和计算机硬件性能之间的关系是有待解决的重要问题。本研究聚焦于林分层次上的静态毛竹林可视化技术探索，中等配置计算机即可满足要求，可以很方便地实现竹林内的漫游，查看任意单株毛竹当前的生长状态和属性信息，也可以对竹林采伐等干扰后的状态进行模拟。在林地可视化中，用竹子坐标点生成的TIN模型表面与真实的地表起伏仍存在一定的差距。另外，由于叠加到TIN表面的林地地表贴图属于栅格图像，当对地表近距离观察时，可能会出现栅格地面，导致图像一定程度的失真，从而影响林地的可视化效果。特别当从林分拓展到包括景观等大尺度的动态场景可视化模拟时，则可以考虑引入注重不同尺度的细节层次技术（LOD），以提高可视化效率。
Visualization for a Phyllostachys edulis stands
-
摘要: 毛竹Phyllostachys edulis林是一种具有较高经济效益和生态效益的特殊森林类型，科学经营管理毛竹林不仅能有效地保护生态环境，而且对发挥毛竹林的经济效益具有十分重要的现实意义。如何提高毛竹林经营管理的技术水平一直是经营者关注的焦点。采用集成OnyxTREE BAMBOO，3ds Max和ArcGIS等种软件，重点将前者的三维建模功能和后者的三维显示功能相结合，实现毛竹林分的三维可视化，旨在为毛竹林科学经营管理提供新的技术方法。选择浙江省天目山国家级自然保护区内的毛竹林为研究对象，建立100 m × 100 m固定标准地，用全站仪测定每株毛竹的基部三维坐标（x，y，z），调查毛竹的胸径、竹高、枝下高和冠幅等因子。根据调查数据，基于OnyxTREE BAMBOO和3ds Max软件，对不同生长状态的毛竹进行三维建模，再导入ArcGIS，在ArcScene中与数字高程模型（DEM）叠加，最终实现毛竹林分的三维可视化。研究取得较好的可视化效果，观察者既可以俯瞰整片竹林，掌握竹林整体状况，也可以深入竹林内部，直观了解毛竹的具体生长状态，如毛竹的死活、弯直、折断、倒伏等，以及地表的起伏状况。单击任一株毛竹，还可以很方便地查看其相关信息包括毛竹的竹高、胸径、冠幅、枝下高、坐标等。此外，借助ArcScene的漫游功能，可以实现竹林内部实时漫游观察，可以设定漫游路径，实现固定路径漫游，以及将漫游动画保存为视频格式进行输出。研究探索出了一套完整的毛竹林可视化技术方法，此方法可操作性强，应用于毛竹林可视化经营管理，不仅可以真实反映毛竹林分的当前生长状况，还可以直观展现竹林采伐或自然干扰后的林分状况，为进一步研究毛竹林动态变化的可视化奠定了基础。Abstract: Moso bamboo (Phyllostachys edulis), a special type of forest with high economic value and ecological benefits, requires scientific management that will effectively protect the ecological environment as it promotes economic benefits. Therefore, to provide a new technology and methodology for improved scientific management of the moso bamboo stand, the 3D modeling function of BAMBOO OnyxTREE and 3ds Max software were combined with the 3D display function from ArcScene software in ArcGIS to realize a 3D visualization of the bamboo stand. Next, a moso bamboo stand with a standard sample plot size of 100 m × 100 m was chosen in the National Nature Reserve of Mount Tianmu, Zhejiang Province. Total Station was used to measure the three dimensional coordinates of each moso bamboo, as well as diameter at breast height (DBH), bamboo height, branch height, and crown width. Based on survey data, 3D models of different growth state of bamboo were made by OnyxTREE BAMBOO and 3ds Max. These were then imported into ArcGIS and overlaid with DEM in ArcScene to make a visualization of the moso bamboo stand. Results showed a favorable visual effect. The observer could not only overlook the entire bamboo stand, but could also virtually go into the stand to see the specific growth state of moso bamboo, such as living or dead, standing or curving, and broken or lodging, as well as the rolling woodland surface. By clicking on a bamboo, relevant information was easily obtained including height, DBH, crown width, and branch height. In addition, real-time roaming within the bamboo stand was realized in ArcScene, and fixed path roaming could be realized by setting the roaming path. Additionally, roaming animation could be saved as a video format for output. From this study, a set of complete visualization techniques for the moso bamboo stand with favorable operability was obtained that could not only reflect the current growth situation of the stand, but could also display the stand status after cutting or natural disturbance, and provided a foundation for further study on the visualization of a moso bamboo stand's dynamic change.
-
表 1 毛竹生长状态分类
Table 1. Classification of moso bamboo growth state
存活状态 受损状态 竹秆形态 状态名称 存活 完整 直立 活完直 弯曲 活完曲 倒伏 活损倒 受损 破裂 活损裂 折断 活损折 多种损伤 活损多 死亡 完整 直立 死完直 弯曲 死完曲 倒伏 死损倒 受损 破裂 死损裂 折断 死损折 多种损伤 死损多 
下载: 导出CSV
-
[1] 余林. 皖南毛竹林密度效应研究[D]. 北京: 中国林业科学研究院, 2011. YU Lin. Study on Density Effect of Phyllostachys edulis Stands in Southern Anhui Province [D]. Beijing: Chinese Academy of Forestry, 2011. [2] 常禹, 胡远满, 布仁仓, 等. 景观可视化及其应用[J]. 生态学杂志, 2008, 27(8): 1422-1429. CHANG Yu, HU Yuanman, BU Rencang, et al. Landscape visualization and its applications [J]. Chin J Ecol, 2008, 27(8): 1422-1429. [3] de REFFYE P, EDELIN C, FRANCON J, et al. Plant models faithful to botanical structure and development [J]. Comp Graph, 1988, 22(4): 151-158. [4] LINTERMAN B, DEUSSEN O. Interactive structural and geometrical modeling of plants [J]. Comp Graph Appl IEEE, 1999, 19(1): 56-65. [5] LINTERMAN B, DEUSSEN O. Interactive modeling of branching structures [G]//BOULIC R, HEGRON G. Computer Animation and Simulation 96. New York: Springer Vienna, 1996. [6] KRSZEWSKI P. An algorithm for sculpting trees [J]. Comp & Graph, 1999, 23(5): 739-749. [7] MAX N, OHSAKI K. Rendering trees from pre-computed Z-buffer views [C] //HANRAHAN P, PURGATHOFER W. Rendering Techniques'95. New York: Springer Vienna, 1995: 74-81. [8] 宋铁英. 一种基于图象的林分三维可视模型[J]. 北京林业大学学报, 1998, 20(4): 93-97. SONG Tieying. A visible stand model based on image of tree [J]. J Beijing For Univ, 1998, 20(4): 93-97. [9] 舒娱琴. 基于林分生长规律的虚拟森林环境的构建研究[D]. 武汉: 武汉大学, 2004. SHU Yuqin. Study on Generation of Virtual Forest Environments Based on the Law of Stand Growth [D]. Wuhan: Wuhan University, 2004. [10] 孙文全. 复杂森林场景的实时绘制技术研究[D]. 福州: 福州大学, 2011. SUN Wenquan. Research on Real-Time Rendering Technology for Complex Forest Scene [D]. Fuzhou: Fuzhou University, 2011. [11] 刘真余, 芮小平, 董承玮. 森林三维真实感建模与可视化LOD 技术研究[J]. 中国科学院研究生院学报, 2011, 28(3): 322-327. LIU Zhenyu, RUI Xiaoping, DONG Chengwei. Study on realistic modeling of three-dimensional forest and visualization based on LOD technology [J]. J Grad School Chin Acad Sci, 2011, 28(3): 322-327． [12] 侯爽, 陈世红. 真实感森林场景的可视化研究与实现[J]. 北京联合大学学报, 2015, 29(1): 58-62. HOU Shuang, CHEN Shihong. Visualization research and realization of realistic forest scene [J]. J Beijing Union Univ, 2015, 29(1): 58-62. [13] 袁玉珠. 地形三维可视化方法研究[D]. 西安: 西安科技大学, 2006. YUAN Yuzhu. Study on the Method of 3D Visualization of Terrain [D]. Xi'an: Xi'an University of Science and Technology, 2006. [14] 宋哲, 刘衍聪, 牛文杰. GIS中TIN模型的实现算法[J]. 计算机应用, 2007, 23(z2): 94-96. SONG Zhe, LIU Yancong, NIU Wenjie. Realization algorithm of TIN model in GIS [J]. Comp Appl, 2007, 23(z2): 94-96. [15] 王永明. 地形可视化[J]. 中国图象图形学报, 2000, 5(6): 449-456. WANG Yongming. Overview of terrain visualization [J]. J Image Graph, 2000, 5(6): 449-456. [16] 李积元. 基于3DMAX三维虚拟景观的设计技术探讨[J]. 大众科技, 2010(7): 36-37. LI Jiyuan. Study on design technology of three dimensional virtual landscape based on 3DMAX [J]. Popul Sci & Technol, 2010(7): 36-37. [17] 吴兆艳, 汤孟平. 基于SketchUp与GIS的森林景观可视化实现[J]. 浙江农林大学学报, 2012, 29(3): 352-358. WU Zhaoyan, TANG Mengping. Forest landscape visualization based on SketchUp and GIS [J]. J Zhejiang A & F Univ, 2012, 29(3): 352-358. [18] 陈军, 唐治锋, 周勇前. 关于城市GIS空间分析模型研究[J]. 经济地理, 1994, 14(5): 5-9. CHEN Jun, TANG Zhifeng, ZHOU Yongqian. Study on spatial analysis model of urban GIS [J]. Econ Geogr, 1994, 14(5): 5-9. [19] 宋仁波, 江南. 基于GIS的三维地质可视化系统设计与实现[J]. 工程地质计算机应用, 2011(4): 22-28. SONG Renbo, JIANG Nan. Design and implementation of 3D geological visualization system based on GIS [J]. Eng Geol Comp Appl, 2011(4): 22-28. [20] 汤孟平, 周国模, 施拥军, 等. 不同地形条件下群落物种多样性与胸高断面积的差异分析[J]. 林业科学, 2007, 43(6): 27-31． TANG Mengping, ZHOU Guomo, SHI Yongjun, et al. Difference analysis of community biodiversity and basal area under different terrain [J]. Sci Silv Sin, 2007, 43(6): 27-31． [21] 高志勤, 傅懋毅. 不同结构毛竹林下植被物种多样性比较[J]. 浙江林业科技, 2005, 25(4): 1-5. GAO Zhijun, FU Maoyi. Comparison of vegetation biodiversity under different structure of Phyllostachys edulis stands [J]. J Zhejiang For Sci Technol, 2005, 25(4): 1-5. [22] 仇建习, 汤孟平, 沈利芬, 等. 天目山近自然毛竹林空间结构与胸径的关系[J]. 生态学报, 2015, 35(12): 1-13. QIU Jianxi, TANG Mengping, SHEN Lifen, et al. 2015. The relationship between spatial structure and DBH of close-to-nature Phyllostachys edulis stands in Tianmu Mountain [J]. Acta Ecol Sin, 35(12): 1-13. [23] 汤孟平, 周国模, 施拥军, 等. 天目山常绿阔叶林优势种群及其空间分布格局[J]. 植物生态学报, 2006, 30(5): 743-752. TANG Mengping, ZHOU Guomo, SHI Yongjun, et al. Study of dominant plant populations and their spatial patterns in evergreen broadleaved forest in Tianmu Mountain, China [J]. Chin J Plant Ecol, 2006, 30(5): 743-752. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2016.05.014
点击查看大图
图(12) / 表(1)
计量
- 文章访问数: 2682
- HTML全文浏览量: 495
- PDF下载量: 380
- 被引次数: 0
