全 文 :林业科学研究 2012,25(2):207 211
ForestResearch
文章编号:10011498(2012)02020705
杉木单木生长可视化模拟系统设计与实现
卢康宁1,张怀清1,刘 闽1,欧阳国良2
(1.中国林业科学研究院资源信息研究所,北京 100091;2.湖南省攸县黄丰桥国有林场,湖南 攸县 412300)
收稿日期:20110520
基金项目:林业公益性行业科研专项(201104028)、中央级公益性科研院所项目(IFRIT200904)、国家林业局948引进项目(2008461)
作者简介:卢康宁(1983—),男,河北任丘人,博士研究生,主要从事林业虚拟现实技术研究.
通讯作者: 张怀清(1973—),男,湖南宁乡人,研究员,硕士生导师,主要从事林业可视化模拟技术研究.
摘要:以杉木为例,采用从可视化模拟角度设计的调查参数调查杉木形态结构,根据实际调查数据对杉木形态结构
特征进行了分析。结合生长曲线控制下的迭代函数系统方法,采用面向对象的设计方法,根据DIRECTX9.0提供的
API及模型渲染技术,在.NET和DIRECTX9.0平台上,使用C#语言建立了杉木单木生长可视化模拟系统,实现了单
株杉木静态和动态的三维可视化模拟。
关键词:杉木;三维可视化模拟;形态结构;迭代函数系统
中图分类号:S791.27 文献标识码:A
DesignandImplementationofIndividualTreeGrowth
VisualizationSystemofCunninghamialanceolata
LUKangning1,ZHANGHuaiqing1,LIUMin1,OUYANGGuoliang2
(1.ResearchInstituteofForestryResourceInformationTechniques,ChineseAcademyofForestry,Beijing 100091,China;
2.HuangfengqiaoForestryFarmofYouCounty,HunanProvince,YouCounty 412300,Hunan,China)
Abstract:Alongwiththedevelopmentofcomputertechniques,theapplicationofcomputerinthreedimensionalvi
sualizationsimulationofplantgrowthprocesshasbeenpaidmoreandmoreatention.TakingCunninghamialanceo
lataasanexample,themorphologicalparametersfromtheperspectiveofthevisualsimulationweredesignedand
thetrees’morphologicalcharacteristicswereanalyzedbasedonthesurveydata.CombinedwithIFS(IteratedFunc
tionSystem)underthecontrolofthegrowthcurve,byusingObjectorienteddesignmethods,andaccordingtoAPI
andmodelrenderingtechnologyprovidedbyDIRECTX9.0,athreedimensionalvisualizationsystemonthegrowth
anddevelopmentofC.lanceolataindividualtreeweredevelopedonthe.NETandDIRECTX9.0platforms,which
realizedthestaticanddynamicthreedimensionalvisualizationsimulation.
Keywords:Cunninghamialanceolata;threedimensionalvisualizationsimulation;morphology;IteratedFunction
System(IFS)
随着计算机技术的不断进步,对植物生长的三
维可视化模拟正受到越来越多的学者重视,近年来
更是取得了长足的发展。植物的可视化模型按照建
模方法和目的的不同分为静态结构模型和动态结构
模型。静态结构模型是利用三维可视化的方法对植
物的形态结构进行重建。如:郭焱等[2]建立了玉米
(ZeamaysL.)的静态结构模型,分析了玉米冠层的
三维结构特征。动态结构模型主要研究植物形态结
构的变化规律,是目前比较流行的建模方法。在当
前研究中最具代表性的是法国农业发展国际合作中
心(CIRAD)建立的 AMAP方法[3]以及加拿大的 L
系统方法[4]。目前国外已经建立了许多虚拟植物软
林 业 科 学 研 究 第25卷
件,如法国的 AMAP、芬兰的 LIGNUM[5]、加拿大的
VLAB以及美国的 SpeedTree等软件。中国科学院
自动化研究所中法实验室与中国农业大学在 AMAP
的基础上建立的 GreenLab模型[6]在植物模拟中也
有广泛的应用。
依据林业可视化技术应用范围及其复杂程度,
可以将林业可视化技术分为单木可视化技术、森林
场景可视化技术和森林景观系统可视化技术3个层
次[7]。其中,单木可视化是指利用科学的计算方法,
将自然界中千变万化的树木形态转换为直观的几何
图形,并将描述其生理结构、生长过程的复杂数据在
计算机中进行直接的计算和模拟,以达到研究树木
生长趋势和木材的收获评价等目的[8]。树木模拟理
论是和树木结构理论、树木三维可视化研究并行交
替发展的[9]。本文以杉木(Cunninghamialanceolata
(Lamb.)Hook.)为例,结合实际测量的杉木形态结
构数据,采用面向对象的设计方法,根据 DI
RECTX9.0提供的 API及模型渲染技术在.NET和
DIRECTX9.0平台上,使用 C#语言建立了杉木单木
生长可视化模拟系统,实现了单株杉木静态和动态
的三维可视化模拟。
1 杉木形态结构
从杉木可视化模拟的角度,结合模拟所需要的参
数设计了杉木形态结构调查因子,主要包括:年龄、分
枝个数、根径、胸径、树高、枝下高、枝下径、东西冠幅、
南北冠幅、枝根直径、枝下距、枝长、仰角、方位角等。
调查因子的确定及具体调查方法详见文献[10]。
对实测得来的杉木一级枝仰角、方位角,二级枝
仰角、方位角应用统计软件进行分析得出,一级枝仰
角基本符合正态分布,主要集中在50° 80°之间,方
位角垂直分布,基本符合均匀分布;二级枝仰角基本
符合正态分布,主要集中在50° 100°之间,方位角
主要集中在(135°,225°]和(315°,45°]区间内,分布
在水平方向上。具体分析详见文献[11]。
2 杉木形态结构模型的建立
本研究采用IFS和冠形曲线相结合的方法来建
立杉木的形态结构模型。设有 N个三维仿射变
换为:
Wn(X)=AnX+Bn =
an bn cn
dn en fn
gn hn i
n
x
y
x
+
un
vn
w
n
其中:X为三维欧氏空间的一个点集,n=1,2,3
…N;An为压缩矩阵,表示绕Y轴旋转γn角,绕Z轴
旋转αn角的两个旋转矩阵和一个收缩矩阵(X、Y、Z
3个方向的收缩率为Sn,取值范围为0<Sn<1)的乘
积;Bn为位移矩阵,un、vn、wn为X、Y、Z上的位移。
通过胸径、根径和树高,应用削度方程来建立主
干的形状;枝下高(枝下距)、仰角和方位角来确定
分枝的位置;枝根枝径和枝长来确定分枝的大小。
根据IFS本身的特点以及调查参数的意义,应用分
枝个数来控制IFS的迭代次数;仰角、方位角来控制
旋转矩阵;枝下高(枝下距)来控制位移矩阵。通过
树高、胸径和冠幅的生长曲线数据来反演不同年龄
的冠形,如图1所示,每个一级枝的年生长量就是代
表一级枝的斜线在与冠形的交点所形成的各个线段
长度。加入影响树木生长的环境因素,通过设定环
境影响率来调整生长模型。树木建模过程如图 2
所示。
图1 冠形及分枝结构图
图2 杉木建模流程图
3 杉木单木生长可视化模拟系统设计
3.1 杉木形态结构模块设计
根据调查参数以及杉木形态结构模型建立的方
法来设计杉木形态结构模块。基于树木形态特征的
分析,采用面向对象的建模方法,将树木主要构件主
干和分枝作为矢量模型,树叶、树枝纹理、树干纹理
802
第2期 卢康宁等:杉木单木生长可视化模拟系统设计与实现
作为栅格模型处理,如图3所示。
图3 面向对象的树木模型
模型中的基本对象是一个基类,由基类派生出
矢量类和栅格类,其中,树叶、树枝纹理、树干纹理继
承于栅格类,主要利用图像的形式来描述其相应的
外观特征。树木模型的主要定义如下:
(1)树模型:{[主干,分枝],[分枝级数、分枝类
型,树高],[纹理]}
(2)主干:{[根径、胸径、高度]}
(3)分枝:{[仰角、方位角、枝径、枝长、枝下距
(枝下高)]}
(4)纹理:{[树叶]、[分枝纹理]、[主干纹理]}
3.2 杉木动态生长模块设计
根据前面提到的杉木形态结构模型建立的方
法,选择树高、胸径和冠幅生长模型作为驱动杉木生
长的模型。本模块采用开放性设计,用户可以选择
模块提供的常见树木生长模型,也可根据实际情况
输入生长模型。选取 Gompertz、Mitscherlich、Logis
tic、Richards、Korf和 Weibul[12]作为供用户选择模
型。其形式如表1所示。
表1 常见的树木生长模型
名称 Gompertz Mitscherlich Logistic
形式
y=k·exp(-ea-b·x)
k,a,b为参数,k>0,b>0
y=A·(1-l·e-m·x)
A,m,l为参数,A>0,m>0,0<l≤1
y= C
1+ep-q·x
C,p,q为参数,C>0,q>0
名称 Richards Korf Schumacher
形式 y=A·(1-b·e
-k·x)
1
1-m
A,b,k,m为参数,A>0,k>0,m>0
y=A·e-b·
1
xc
A,b,c为参数,A>0,b>0,c>0
y=k·exp -b
x( )c
k,b,c为参数,k>0,b>0,c>0
3.3 杉木单木生长可视化模拟系统总体框架设计
杉木单木生长可视化模拟系统通过给定的参数
和模型应用计算机模拟出杉木的静态以及动态的三
维形态。其系统框架如图4所示。系统层是指系统
开发与运行的硬件和操作系统;开发层是系统的开
发环境以及开发工具,建立在系统层之上,主要负责
模型的建立与程序的开发;应用层是系统的运行层,
也就是提供给用户使用的层次。在应用层用户可以
观察和控制运行结果以及场景。
4 杉木单木生长可视化模拟系统实现
4.1 关键技术
计算机图形学在本系统的开发中起着举足轻重
的作用。MicrosoftDirectX旨在使基于 Windows的
计算机成为运行和显示具有丰富多媒体元素(例如
全色图形、视频、3D动画和丰富音频)的应用程序的
理想平台。其中Direct3D是负责三维图形的一个组
件,集成开发了 DirectX程序三维部分的所有 API
(ApplicationProgrammingInterface)函数,也是 Di
rectX的最重要、最复杂的组成部分。本系统中的三
维图形显示都是借助于Direct3D的功能来实现的。
图4 杉木单木生长可视化模拟系统框架
902
林 业 科 学 研 究 第25卷
4.1.1 Direct3D Direct3D是基于微软的通用对
象模式COM(CommonObjectMode)的3D图形API。
可绕过图形显示接口(GDI)直接进行支持该 API的
各种硬件的底层操作,大大提高了图形处理速度。
D3DDevice分为硬件抽象层设备(HALDevice)和参
考设备(REFDevice),主要提供矩阵变换、三维图像
渲染和三维图像光栅化显示等功能。Direct3D在
完成坐标变换、光照处理和纹理映射后,将绘制结
构的值存储于颜色缓冲区中,最后通过图形刷新,
颜色缓冲区的内容显示在屏幕上。本系统的三维
显示部分正是利用 Direct3D提供的图形操作技术
来实现的。
4.1.2 纹理贴图技术 用真实的自然界中的纹理
来表现模型的表面,可以使三维图形显得更逼真、更
自然。DirectX提供相应的函数对纹理进行控制,如
贴图的方式,纹理坐标的转换等。在DirectX中使用
关键色透明技术,通过设置渲染状态 ALPHATSTE
ABLE和参考值ALPHAREF进行 Alpha测试。程序
根据 Alpha测试函数比较当前像素的 Alpha值与
ALPHAREF的值,由比较结果决定是否绘制该像
素,如果返回 TRUE,则通过测试并绘制像素,反之
则不予绘制[13]。本系统中主干、分枝以及树叶的纹
理都是基于这种技术实现的。
4.2 系统功能简介
系统主界面如图5所示。主要分成菜单栏、工
具栏、控制面板和显示区4个部分。控制面板包括4
个选项卡,分别用来控制树木的整体信息,如胸径、
树高、年龄、冠幅等;形态结构信息,如主干、分枝、叶
子以及相应的纹理设置等;生长模型选择以及模拟
结果的环境状况,如天空、土地等。工具栏包括最常
见的操作命令,如新建、保存、放大、缩小、平移、旋
转、节点显示、网格显示等。生成结果可以通过图片
的形式或者点文件的形式输出保存。
4.3 模拟结果
在湖南省黄丰桥国有林场广黄分场选取具有代
表性的杉木进行形态结构测量。通过对1株13年
生的杉木调查得出其形态结构参数如表2所示。
表2 13年生杉木形态结构调查
调查指标 取值 调查指标 取值
树高/m 9.2 一级枝个数/个 64
胸径/cm 18.4 枝下高/m 1.78
根径/cm 23.2 仰角/° 45° 80°
东西冠幅/m 4.3 南北冠幅/m 4.2
图5 杉木单木生长可视化模拟系统主界面
其中,顶端很小的一级枝并没有计算在内,一级
枝方位角基本上符合正态分布。将调查的参数输入
到系统,模拟得到杉木形态结构和实际杉木形态结
构对比如图6所示。
图6 杉木形态结构模拟图和真实图的对比
采用生长模型来控制杉木树高、胸径和冠幅的
生长情况,应用环境因子来简化环境因素对树木生
长的影响,通过调整环境因子大小来改变树木的生
长率和生长曲线的形状,如图7所示。
图7 树高、胸径和冠幅的生长曲线图
012
第2期 卢康宁等:杉木单木生长可视化模拟系统设计与实现
图8为10年生和15年生杉木的模拟图和真实 图的对比。
图8 模拟生成的杉木与真实杉木的对比
5 结论与讨论
在本研究中应用调查得到基础数据做为系统模
拟的入口参数,使用树木生长曲线改善传统的 IFS
方法进行形态结构建模,通过计算机技术完成了杉
木单木生长可视化模拟系统的设计与开发,实现了
单株杉木静态和动态的三维可视化模拟。纹理贴图
技术,环境渲染技术的使用大大提高了本系统的模
拟效果。选取调查的杉木为例,应用建立的系统进
行了树木三维形态结构模拟。通过设定生长曲线的
参数以及环境影响因子来实现杉木动态生长变化。
环境因素和树木生理因素对树木生长的影响是
目前树木研究的重点和难点。在本研究中需要进一
步的改进。从生物量生产和分配的角度来研究树木
生长与环境交互的作用是本研究的发展方向。另
外,本系统的输入输出功能有待进一步加强,在以后
的工作中将实现本系统输入输出 obj等格式的文
件,提高与成熟三维软件的交互性,把系统产生的树
木模型导出到3DMAX等行业软件中进一步操作。
由于树木形态结构的差异性、生长的随机性以及计
算机运算的有限性,建立通用的树木可视化模拟系
统也是本方向进一步研究的内容。
参考文献:
[1]郭 焱,李保国.虚拟植物的研究进展[J].科学通报,2001,46
(4):273-280
[2]郭 焱,李保国.玉米冠层三维结构研究[J].作物学报,1998,
24(6):1006-1009
[3]DeRefyeP,EdelinC,FranonJ,etal.Plantmodelsfaithfultobo
tanicalstructureanddevelopment[C].SIGGRAPH'88Proceedings
ofthe15thAnnualConferenceonComputerGraphicsandInteractive
Techniques,1988,22(4):151-158
[4]LindenmayerA.Mathematicalmodelsforcelularinteractionsinde
velopmentI.Simpleandbranchingfilamentswithtwosidedinputs
[J].JournalofTheoreticalBiology,1968,18(3):300-315
[5]PertunenJ,SievaEnenR,NikinmaaE.LIGNUM:amodelcombi
ningthestructureandfunctioningoftrees[J].EcologicalModel
ling,1998,108:189-198
[6]HuBG,DeRefyeP,ZhaoX,etal.GreenLab:aNewmethodolo
gytowardsplantfunctionalstructuralmodelstructuralaspect[C].
ProceedingsofPlantInternationalSymposiumonPlantGrowthMod
eling,Simulation,VisualizationandTheirApplications.Beijing,
China,2003:21-35
[7]侯军岐,杨 静.林业可视化技术推广路径选择研究[J].林业
经济问题,2010,30(3):211-214
[8]龙 洁,苏喜友.国内树木三维可视化研究进展[J].林业调查
规划,2007,32(6):44-47
[9]景向欣.樟子松人工林单木动态生长三维可视化模型的研究
[D].哈尔滨:东北林业大学,2007
[10]卢康宁,张怀清,刘 闽.杉木形态结构可视化模拟调查方法
研究[J].中南林业科技大学学报:自然科学版,2010,30(1):
34-40
[11]卢康宁,张怀清,刘 闽.基于实测数据的杉木构筑型研究
[J].林业科学研究,2011,24(1):132-136
[12]张建国,段爱国.理论生长方程直径结构模型的研究[M].北
京:科学出版社,2004
[13]吴 谦,张怀清,陈永富.基于生长规律与图像相结合的树木
叶片动态模拟[J].林业科学研究,2008,21(增刊):122-125
112