您好,欢迎来到微智科技网。
搜索
您的当前位置:首页一个基于中心差分技术的实时高度场可视化方法

一个基于中心差分技术的实时高度场可视化方法

来源:微智科技网
自动化技术与应用》2009年第28卷第l 1期 计算机应用 Computer Applications 一个基于中心差分技术的实时高度场可视化方法★ 徐林 (北华大学计算机科学技术学院,吉林吉林13201 1) 摘 要:本文提出了一种基于中心差分技术进行高度场镶嵌的系统方法,结合视点相关技术实现了多分辨率模型的生成,以适应实时 应用的需要。此外,结合数据的预处理使本文的方法也能适用于真实采样数据的图形生成。本文的方法可用来实现对数据高 度场进行实时图形显示的目的。与同类方法相比,本方法灵活性高,能够适于多种应用,其视觉质量好并且执行效率高。 关键词:图形学;科学计算可视化;高度场;细节层次 中图分类号:TN9I1.73 文献标识码:B 文章编号:1003 7241(2009)11 0038-05 A Real・-time Scientific Visualization Method of Height Fields Based on Central Diference Technology XULin (Computer Science&Technology,Beihua University,Jilin 13201 1 China) Abstract:In this paper,a systemic method using central difference technology is presented tO tessellate height fields and generate view—dependent multiresolution models tO meet the necessity of real—time applications.In addition,through data pre— processing,this method adapts to figure generation of real data sampled.Moreover,this method also adapts to realize the purpose of real time figure displaying to data height field.Comparing to previous methods,the method presented in this paper provides good flexibility in various applications.It works well in realtime display application while has good visual quality. Key words:computer graphics;visualization of scientific calculation;height field;level of detail 1 引言 科学计算可视化技术(Visualization in Scien- tific C omputing)是对计算及数据进行探索,以获得 对数据的理解与洞察。它为抽象的观测数据赋予了 图形化解释。 来表现。针对如何由采样数据恢复空间上连续存在的 数据场并以曲面形式进行表现的问题,本文利用一种适 于利于标准化OpenGL图形接口实现的可视化方案,并 采用计算机图形学中可视效果较好的双三次曲面表现 可视化图元。 相比之前多采用的线性插值曲面拟合技术,本文 在一定的数字化描述形式下,使用参数曲线曲面进 行3 D建模,从而可达到物体的虚拟重构或可视化映射。 但如何解决真实感、观测精度与生成速度的矛盾成为 自由曲面普及应用的重要课题。 二维网格采样数据场的可视化表现形式有多种,这 里选择其中较常用的最具直观表现力的三维映射曲面 提出一种从 2域正交网格分布数据集自动重建平滑 具自适应细分特性表面的实时渲染的系统方法。此方 法克服了线性插值图形视觉质量表现较差的缺点。通 过提出的网格简化的高效算法,结合从真实采样数据 构建表面平滑的高度场模型的预处理方法,避免求解 多项式微分方程组而得到拟合曲面以降低计算复杂度 ・基金项目:北华大学科研发展基金北华大学科合字(2 0 0 8)第2 8号 收稿日期:2 0 09—06—1 2 【ll。较之此类技术,方法以动态方式产生LOD模型,利 用了快速中心差分操作,可支持帧率和建模策略的灵 计算机应,El 《自动化技术与应用》2009年第28卷第11期 Computer Applications 活调整,并通过支持表面细分特性进一步提高了视觉 即.Il((以+b)/2),遵循以下所述进行。 质量。方法执行效率高,完全可满足应用于实时采样 首先,求 和h(b)的均值,通过此项操作,就可得 数据的高质量可视化。因此,能更好适应实时高度场 到二者之间连线上的一点。 的图形表现需求。 接下来,对于 (日)和^ ( ),即函数JIl( 在a和b点 的二阶导数,求其均值可得到中点的二阶导数 2基于中心差分法的单面片生成方法 ’( + /2)。继而求得 ,它代表 和 (易)之间参数 以本章讨论的中心差分法为基础,本文最后将给出 距离的一半,即 ~a)/2。然后,计算一 ^ ((a+b)/2)/2。 由Bezier面片构成高度场(地形)系统的方法。因此,从 再将其与上面求得的均值点加和,就得到了欲求的函数 本节开始,将讨论如何找到一个面片镶嵌方法,它应该 中点,虽Ⅱ|Il((口+易)/2)。 较之以常规曲面生成原理为基础的方法(例如OpenGL 至此,可以确定,中心差分是一种镶嵌面片的更好 求值器)更好。 方式。这里,将在推导中利用泰勒多项式的性质。首先, 直接依据Bezier曲面定义的表达式对每个顶点进 对一些概念作一说明。这里,将中点称为矗(“),du为中 行计算,这是一种非常低效的方式。因此,应该研究如 点与端点之间的参数距离。这就得到第一个端点,将其 何设计一种更好的镶嵌方法。以此为目标,以下开始探 表示为h(u—du),以及另一个端点h(u+du)。 讨如何得到一种真正实用的执行方式。 以此为基础,若给定中点,对泰勒多项式求值, 下面将就一种更好的特殊镶嵌算法进行讨论,实际 便可以此方式等价的得到端点。这样就得到以下两 上,应用这种算法囊中心差分,将通常的Bezier面 个等式: 片表达形式。在计算机真实感图形的算法基础I21中提及 + : )+dtt(h。㈤+等矗 )+等^ (“) (2) 了一种形式的中心差分算法,但相比其中提及方法,本 h(u-du):|Il(1‘)+-du(h。(“))+(--du)2 h (“)+(-d u)3 h (文作了较大的修改。 。“) (3) 在中心差分中使用的基本数学结构是泰勒多项式。 在(3)中可将右侧的含有一du的几项中的负号抽取 在一个参数U附近的函数h(x)的泰勒多项式的构造公 出来,得到以下两个等式: 式如下: +du): )+du(h ( + h (卅华h…(“) (4) 在本文所述情况下,处理对象是立方函数。因此, 任何超过三阶的导数都将为零 因此,如果三次函数是 +du) -du(h ㈣+譬^” )一譬|ll…(“) (5) h(x)。其泰勒多项式形式如下: (4)、(5)式自身并未提供求得h(u)的足够信息,但仅通 (“+du):^(卅du(h ( + h”(“)+ h (“) (1) 过将二者加和简化,便能得到一个形式很好的等式。 二 叶 h(u+du)+h(u—du)=2h(u)+du h”(U) (6) 至此,可能会注意到一个问题,就是这些所起的作 正的dU与负的dU互相消去,这样就得到了一个 用就是给定一个三次多项式,结果又得到一个三次多项 很好的形式。然而,端点并非求值目标,所需要求得的 式。因此,这就可能会产生疑问,即式(1)的变换有何助 是中点。 益。但中心差分所起的作用正是以此为基础。中心差分 对得到的等式移项重新组织得到: 的思想在于,给定一曲线段的端点,以此为基础求其中 2h(u)=h(u+du)+h(u—du)一du ”(“) (7) 点。然后,在后继过程中,继续以这种方式,利用两边得 等式两边以2除之并作分解,变为: 到的中点作为新的端点,进行递归并求其中点,如此不 矗 )= h(u+d“)+h(u—du)du .1l (“) 断进行,直至得到所需的曲线。 2 2 下面来看一下这个过程是如何进行的,然后再来探 等式右侧第一个较大分项部分恰为两端点的平均, 讨求导问题。如果已知函数h(x)上的两点,比如说 求到它很容易,因为假设拥有任何所需的端点信息。但 和h(b),而且想找到h(a)和h(b)之间的函数中点,或者 求得右侧的第二个部分相对来说要困难一些。它依赖 《自动化技术与应用》2009年第28卷第ll期 计算机应用 Computer Applications 于中点的二阶导数。问题在于如何从端点出发得到它。 注意到二阶导数是一个线性函数。因而从端点出发,推 导它将会很容易;仅需将端点处h的二阶导数作平均, 就已经得到了结果: ”( ): 二 (9) 现在,将其代回式(8),得到: = h(u+du)+h(u一 2 (矗’ + +∥ 一 ) 4 (10) 以上过程较为复杂,但通过分析它可以得到一些有 用的结论。它表明通过对端点作平均,并加上端点二阶 导数的一个带权均值,就可以得到中点。至此,知道了 如何获得中点的二阶导数(刚刚在(9)中实现)。因此,当 欲对曲线的其余部分求值时,将有足够的信息进行递归 求值操作。 图1 中心差分镶嵌的一级递归 v向中点可使用大致相同的方式产生。这时,可忽 略u,并利用端点及其关于v的二阶导数来求得其v向 尽管如此,以上所述仅仅是一条曲线的情况,原因 在于 仅为一单变量函数。将这种思想应用于面片 中点及中点关于v的二阶导数。 为使这个过程看起来一目了然,在求值过程进行 中,明确当前阶段都得到了哪些点的什么信息是颇为有 效的方法。下面来看一下现阶段哪些点的何种信息是 已知的。 才是目的,而曲面是双变量的函数。在开始就曲面的 情况讨论之前,下一节先来讨论中心差分的一个有用 的性质。 3适用于高度场建模的中心差分技术 使用中心差分的方式,利用至此所讨论的内容,已 u,v)角点:,( ),百32f(u,v)32f(dIU向中点:,( ), du ,百u dv ; V向中点:,( , ), ;中心点:暂无。 经足以执行一条曲线的镶嵌。现在要讨论的是如何应 用中心差分来镶嵌曲面。较之一个变量的情况,执行二 维的中心差分更为困难一些。这里将问题重新作一描 在考虑如何求中心点之前。应该保证U向中点和v 向中点的完整性。即,若进行递归,需保证结果之中新 点的信息种类与从角点获得的信息种类完全相同。否 则,对于四个较小方形区域上的递归,其信息完整性将 无法满足要求。所以,目标应确定为寻找f在u向中点 对于v的二阶偏导。可是,如何做到这一点,目前却还 不是很清楚。对如何沿u向对变量v插值的问题,一条 线索来自Bezier曲面的表达式: 述。这次,期望能由一个面片之上的四个角点信息出 发,找到面片各个角点之间的中点及面片中心点的信 息,所要求的信息种类要求与角点的信息种类相同。其 图解见图1。 观察图1,可以明显地看出,给定各角,在此一层 次上,若能得到中点,则可在较小的四个方形的层次 上以同样的方式递归,等等,直至得到理想的面片模 型。当前,既然还不能肯定需要角点的何种信息,可 以首先来探讨一下需要什么来求得已知点之外的另 外五个点。 u向中点可通过一维的中心差分方式产生。这时可 ∑∑Po口 (“) ;(V) (11) 不妨先沿u向寻找关于v的二阶偏导。这时, 必须得到曲面的关于v的二阶导数,由上面的等 式,得到: 忽略v,因为它在此方向上保持不变,并利用端点和其 关于u的二阶导数来得到其u向中点及中点关于u的 二阶导数。 砉 ㈨ 现在,关于u的基函数,B(u),为一三阶函数。既然 它保持原函数形式未变,那么得到的曲面关于v的二阶 计算机应用 《自动化技术与应用》2009年第28卷第1 1期 Computer Applications 导数是一个u的三阶函数。对于任意u的三阶函数,利 中点的对应复合偏导值实现。 用一维中心差分,即可在U向对其插值。 最终得到了所有点的全部信息: 在此情况下,被插值的函数是关于v的二阶导数, 角点、u向中点、v向中点、中心点: 将其记为g(u)。 f(u,v , , , 。 g(“): (13) ay 因此,给定各自包含五个值的四个角点。通过使用 因而,若知道端点的g(u)和g关于u的二阶导数值, 中心差分,可以引出u向中点、v向中点以及中心点。以 就可对g(u)进行插值。这意味着在u向中点,欲得到关 此为基础,可在四个由新点构成的较小方形的层次继续 于v的二阶偏导需要以下数据: 进行递归。整个过程从显式地计算四个面片角点的所 a f(H,v)a f(u,v) — 一’ 需信息开始,继而递归进行。迄今已经明确了递归进行 在各个角点,已经有了第一个数据的值,即关于v 的步骤,但还有一个问题需要解决,即递归应进行至何 的二阶导数。还必须添上第二个值,即v的二阶偏导数 种程度才应终止。 的关于u的二阶偏导。 在式(1 0)中,通过求端点平均值并加上另一项—— 至此,得到了在U向中点关于v的二阶偏导。以同 二阶导数项,便得到了中点。值得注意之处在于,如果 样的方式,通过将变量互换,容易得到在v向中点关于u 二阶导数项为0,就意味着中点处于端点之间的连线上。 的二阶偏导。所以,关于u的二阶偏导是v的三阶函数, 因此,将该带权的二阶导数项称为中心差分的非线性项 因而将以下内容增添到角点即可求得它: 它决定着中点距离端点之间的连线有多远。 : : : 中心差分的非线性项是对于正在镶嵌点的表面弯 a“。 ’av 2a“ 曲程度的一种估计。可以使用中点生成过程中的非线 同样,第一项在各个角点已知。此处复合偏导可以 性项来决定正在处理的子片的曲率。进而可以指定某 轮换,即: 个阈值,若所有中点的非线性项的大小比阈值小,保持 a: 子片由各角决定的两个三角形构成的状态并返回。由 H av av u 2 O 2 (1、 4)最后,注意到复合偏导是u和v的线性函数,所以, 此,若阈值很大,最终将得到一个包含仅两个三角形的 通过将角点的复合偏导值作平均,即可得到u向中点和 粗糙面片。若阈值很小。则面片会极度细化,这样视觉 v向中点的复合偏导值。 质量很好但运行很慢。因此,选择一个恰当的阈值只是 一通过以上分析,明确了各角点所需的全部信息,进 个试验的问题,并且可将其作为系统设置功能的一个 而,也可在u向中点和v向中点得到相同种类的信息。出 有用的候选项。此外,需要使用一种恰当的方式,在镶 于保留线索的考虑,以下列出当前已有的信息: 嵌中使距离和视角的因素也能够被考虑进去,并且期望 角点、u向中点、v向中点: 执行要快。理想情况下,应能够知道非线性向量以像素 为计算单位的大小。 f(u,v , , , :暂无。 以前述基于中心差分法的单面片生成方法为基础, 现在,仅需指出如何得到中心点信息即可达到目的。 在对面片拼接构成地形的相关方法选择时,这里考虑到 回顾图l,正如每一v向中点可由其两个角点推得 图形真实感表现要求,对于面片间的连续性的问题,即 一样,中心点可从u向中点推得。所以可以利用u向中 共享邻接边两侧的面片的切向量应该相同。这样,即便 点和它关于v的二阶偏导来得到中心点和它关于v的二 在面片间,地形的观察效果也将是平滑的,采用【 ]中的方 阶偏导。然后,通过求v向中点关于u的二阶偏导的平 法,通过Cardinal插值生成曲面两个参数方向上的切向 均值,就可得到中心点的关于u的二阶偏导。最后,可得 量,然后进行Hermite曲面至Bezier曲面转换的方式解 到中心点信息中的最后一个值,关于u的二阶偏导的关 决了这一问题。此方法可避免使用一些计算复杂度更 于v的二阶偏导(从式(1 4)得到的值),这可通过平均v向 高的方法,如求解微分方程组【引,而显示质量相当,仍可 《自动化技术与应用》2009年第28卷第11期 很好的满足显示应用的需要。 计算机应用 Computer Applications 多边彤量臣0 ^O 0  O0 0 1 ^0 O 1 ¨ 与高度场本身起伏程度相关的LOD算法来决定地形 应有的细节程度,使需要绘制网格模型的三角形数大 大减少。 为了说明多边形裁剪和动态视点相关LOD建模的 有效性。这里给出对多边性数量、帧率等的定量分析结 果。距离衰减参数部分未作变换及保持线性变化。对于 以指定差分递归次数处理得到的高度场数据(测试中得 到的此项数据规模为1 X 10 个网格多边性)。以像素为 单位,计算每一正交网格中非线性项最大值的屏幕投影 大小,并取所有网格中此项数据的最大值。阈值的测试 用值以此为上限变化,并以投影像素值计量,记为 。 对于面片间的拼接,由于即相邻面片控制点定义 之间存在关联性,无法进行完全局部性的控制。对此, 结合应用了一种对数据进行预处理的方法,它通过 Cardinal插值生成曲面两个参数方向上的切向量,然 后进行Hermite曲面至Bezier曲面转换的方式解决了 这一问题。 首先将生成模型多边形数量作为阈值 的函数,检测其 变化。结果如图2所示。 在面片生成过程中,基于Bezier曲面的中心差分方 法为得到视觉质量良好的模型打下了基础。利用非线 性项决定镶嵌精度,不仅使判断和建模计算变得方便高 效,而且,结合考虑视点的方位等因素对视觉的影响,这 里的动态建模方法同时具有可随视点接近不断细分的 特性【 】,在视觉的细节度表现方面可以不受原始模型规 模的。并且执行效率得以保证。 0 2 —4 ——6  参考文献: 【l】IN—HUGH CHOI,MIN—YANG YANG,WON—PY0 HONG and TAE-SUNG JUNG.Curve interpolation with 譬线性l,I授髟大小 俏謇) -.●_一variable feedrate for surface requirement[M】.Publisher: Springer-Verlag London Ltd,Issue,2005,25(3-4):325-333. 【2】彭群生,鲍虎军等.计算机真实感图形的算法基础【M】. 图2 的变化对多边性简化数量的影响 北京:科学出版社,1999. 【3】徐林,孙长嵩.一个基于映射曲面形式的二维网格采样 数据场可视化方法【M】.自动化技术与应用.中国自动化学会, 2006,l1:16-l7. 当 值在2以内时,其值非常小,画面在此时没有明 显的跳变。在加入纹理与光照时,画面质量几乎与未经 L O D简化的画面没有区别。而此时多边性的数量明显 减少。取 =2,并将以随机方式确定视点的运动轨迹。 对应于图2的测试集,应用本文方法,在l20秒时间内, 获得了每秒不少于3 0帧的结果。 【4】【美】金凯德(KINCAID,D.),切尼(CHENEY,w.).数值 分析(英文版第3版)[M】.北京:机械工业出版社,2003.4. 【5】LINDSTROM,P.,K0LLER,D.,RIBARSKY,W., HODGES,L.,FAUST,N.,and TURNER,G.Real-time, continuous level of detail rendering of height fields[C】. Computer Graphics(SIGGRAPH’96 Proceedings),1996, 4 结束语 解决真实感、观测精度与生成速度的矛盾是自由 曲面广泛应用的重要课题。本文主要解决的问题是数 据高度场如何应用参数曲面进行三维模型的生成,同时 以保证实时显示作为要求。 本文讨论单面片的生成方法,在分析了参数曲线 的泰勒多项式表示形式的基础上,引出了中心差分的 概念。提出了结合基于Bezier曲面的中心差分技术 109-l1 8. 作者简介:徐林(1 9 7 9一),男,讲师,工学硕士,研究方向:计 算机图形学、数字图像处理与计算机视觉等。 和视点相关技术的递归镶嵌生成单面片的方法。结合 

因篇幅问题不能全部显示,请点此查看更多更全内容

Copyright © 2019- 7swz.com 版权所有 赣ICP备2024042798号-8

违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务