system: [SYSTEM]

user: 3. 相关工作 目前，流可视化的主要研究机构有：美国国家宇航局Ames研究所和计算机科学与工程应用研究所主要针对计算流体力学数据进行研究；奥地利维也纳科技大学、德国柏林大学和意大利航空研究中心主要研究动力学系统产生的流场数据可视化；美国佛罗里达州立大学、美国加州大学、荷兰埃因霍温大学、加拿大西蒙非沙等大学主要从事基础理论研究；国内主要从事理论研究的单位有清华大学、上海交通大学、浙江大学、中科院CAD开放实验室、国防科技大学等；中国工程物理研究院、中国船舶科学研究中心、中国空气动力研究与发展中心等单位结合具体应用背景进行流场可视化的应用研究。 流可视化的算法很多，依据矢量数据的映射方式不同，一般可以分为三大类。 3.1. 几何方法 第一大类是几何方法，主要包括点图标法和矢量线法。 ¨ 点图标法：这是一种最简单直观的矢量数据映射方法，直接显示流中各点的速度矢量构成的“速度场”，优点是易实现、速度快，最常用的显示元素是箭头、椎体、有向线段等，所有这些点图标常被称作刺状体(hedgehog， glyph) [ 3 ] 。在二维矢量场中，箭头可以较好地反映出矢量大小和方向信息，但是需注意放大后使箭头不能相互覆盖；在三维矢量场中，普通箭头方法容易产生二义性，为克服这一点人们设计了一些特别的箭头，如有阴影的箭头、不透明有一定粗细的三维多边形等。这类方法的缺点是箭头疏密难以把握，而且难以揭示数据的内在连续性和某些结构特征。 ¨ 矢量线法：源于实验型流场可视化技术，通常被认为是对矢量场局部细节的中等程度抽象，包括流线、时线、脉线、迹线 [ 3 ] 。流线是用于稳定流场的最主要方法，其他方法一般用于非稳定流场，流线的计算可以采用基于流函数的构造方法或者像Euler方法、Euler修正法、二阶Runge-Kutta法、四阶Runge-Kutta法等数值积分方法 [ 4 ] [ 5 ] ；而迹线、时线和脉线均可通过求解迹线方程求解。与点图标法相比，矢量线能表示场的连续性，绘制效果较好，但是绘制质量依赖于初始点的选取。为充分利用图形学中的光强、色调、消隐和透视投影所提供的三维信息，还可将矢量线扩展为矢量管、矢量带或者矢量面等，而且这些方法可以可直观地显示流的一些特性，如可以用不同颜色、不同宽窄/粗细的色带/管子描述流场的速度的大小，而用色带/管子的走向描述速度的方向。近年来出现了一种被称为“流球”的流可视化技术。这一技术通过计算每一点对周围一个球形区域的贡献来获得一个流场的曲面表示，其中的可视化元素既可以是一小段曲面，也可以是一段圆管。脉线也可以用来描述流场的粒子运动情况，可以用脉线描述流场中粒子的发射起点和运动轨迹。还可以用动画来演示粒子在流体中的运动，称为“粒子动画” [ 6 ] [ 7 ] 。 几何方法 [ 6 ] [ 8 ] 生成三维场的可视化图像是离散的，而且使用者不得不自己决定箭头、流线或粒子的起点位置。颜色编码法在几何方法中起着非常重要的作用，对于标量场数据，比如温度、压力、密度等，可以通过在标量值和亚瑟之间建立一一映射关系，用颜色值的变化显示标量值的变化。 3.2. 纹理法 第二大类是纹理法 [ 8 ] ，纹理是色彩按一定方式排列组成的图案，兼具形状和色彩属性，纹理映射在图像空间进行，可连续生成具有高分辨率的细致图形。纹理映射技术是计算机图形学最常用的技术之一，利用纹理的灰度相关性来反映流体的运动特性。此类方法提供了流场的全局密集表示，可以很好地描绘流场的细节。典型方法有点噪声、线性卷积和IBFV等硬件加速算法。 ¨ 点噪声(spot noise)方法：Van Wijk将纹理引入流可视化 [ 2 ] ，他提出的点噪声(spot noise)方法中根据局部流的方向对周围数据点的取向进行调整，并通过在流的方向上绘制一些椭圆形的纹理“斑点”来实现流的可视化，并通过变换和分布控制来反映矢量场的特性。 ¨ 线积分卷积法(Line Integral Convolution，简称LIC)：在基于纹理的方法中，以Cabral和Leedom提出的线积分卷积法(Line Integral Convolution，简称LIC) [ 8 ] 最为典型，该方法对噪声纹理滤波卷积运算生成具有一定方向特征的纹理，即使在矢量方向变化很大的区域也能揭示出矢量的方向，可以较好地表达出矢量场的细节，卷积后的图像具有像素分辨率的连续性，描述流数据场的信息非常丰富，但该方法只考虑了沿速度方向那条折线段上的像素点对该点的作用，没有考虑与速度垂直方向上邻近点的影响，因而高频噪声较大。 这两种方法都可以在一个动态系统中产生对流体动力学特性的整体显示。其他还有少数与纹理相关但又不能简单归为以上的技术，如纹理溅射(Texture Splats)。有很多研究者对LIC算法进行了改进和扩展，提出了像加速LIC、并行计算LIC、曲线网格上的LIC (LIC on curvilinear grids)、三角化表面的LIC (LIC on triangulated surfaces)、多重频LIC (multi-frequency LIC)、有向LIC (oriented LIC)、增强 LIC、染色LIC (dyed LIC)、体素LIC (volume LIC)、超LIC (Hyper LIC)等 [ 8 ] 。另外，还有人引入间接像素纹理选择、加/减纹理混合技术等加速其积分和卷积速度；引入非线性扩散滤波来增强原始的低通滤波的效果；加入纹理映射来处理一些时变的、时变的流，用迹线代替流线作为积分线等。由于普通线性卷积生成的纹理偏暗，对比度太差，于是出现了许多纹理增强算法，如直方图均一化、双重卷积、二维高通滤波、偏微分、特征增强、属性映射、特殊区域增加权重等。 ¨ 硬件加速算法：具体来说就是在通用计算图形处理硬件的基础上进行编程。由于纹理法与GPU (Graphics Processing Units)所支持的纹理操作紧密结合，许多技术都应用GPU加速实现实时绘制。Jobard 等采用间接像素纹理选择技术设计了一种有效的显示管，Weiskopf等开发了像素纹理单元表征三维时变矢量场的潜力；Van Wijk在纹理平流技术和基于GPU的硬件加速的基础上提出了一种 IBFV (Image Based Flow Visualization) [ 2 ] 流场绘制技术，该方法通过向前的投影网格映射实现噪声纹理的流动，并形成新的帧图像，再结合连续帧图像的比例混合实现了流场体绘制；IBFV已被扩展运用到三维曲面矢量场和三维矢量场的可视化；UFLIC (Unsteady Flow LIC)是Shen和Kao提出的一种基于粒子模型的硬件加速算法，并采用二维高通滤波来增强纹理质量，并有效避免粒子连续采样卷积导致的值扩散；随后Liu和Moorhead提出了利用迹线提出了一种新的粒子管理机制AUFLIC (Accelerated UFLIC)，极大的加速了算法性能；他们的方法可以顺利过渡到VAUFLIC，G. S. Li给出了这一算法的GPU实现 [ 8 ] 。 这里还涉及动态纹理，动态纹理可以看作是二维静态图像纹理在三维时空维的扩展.动态纹理现象并不陌生，我们常发现这些特殊的图像特殊的图像序列如水波纹、气体和烟雾的运动，涡流等。计算机图形学和机器视觉都涉及动态纹理的研究，前者侧重于视频编缉和人造视频方法的研究，后者侧重于运动特征的识别和视频理解等。 相对于其他流场可视化方法，纹理法是一种有效的、多样化的矢量可视化方法，运用非常广泛，而且有许多拓展和优化，它可以致密地表征整个流场，可以表征特定的纹理特征和动态流场，而且没有初始点选取的问题，因而获得了更多的关注，但是这类方法计算复杂度高，因而研究者们在提高算法速度方面开展了许多研究工作，尤其是将纹理法从二维向三维扩展之后，这一问题更为突出。此外，研究者们关注比较多的问题还有：从稳定的矢量场到时变场的扩展、矢量场的大小和其他标量的显示、以及与矢量线技术的结合等。 3.3. 特征法 特征可视化一直是矢量可视化研究的一个热点，所谓特征具有两种含义：一种是指流场中有意义的形状、(拓扑或特征)结构、变化和现象，如涡流、激波等；另一种指分离出来的用户感兴趣的区域。特征可视化就是通过对场中这些特征重点可视化，减少可视化数据映射的数据量，却保持一定的准确性，有助于忽略许多冗余或不感兴趣的数据。 特征法主要有以下求解思路 [ 9 ] ： ¨ 直接映射，并不抽取特征，而是在直接显示原始数据的过程中通过选择合适的映射技术将特征表现出来。 ¨ 特征识别，提取特征数据并识别出用户感兴趣的部分。 ¨ 拓扑结构分析法，流场拓扑结构分析理论最早由Helman和Hesselink提出，以临界点理论为基础，认为一个矢量场的拓扑由临界点和连接临界点的积分曲线和积分曲面组成，更一般的矢量场结构分类法由Chong等归纳，临界点在矢量场中即为矢量为零的点，需要计算其位置并分类，具体算法由Globus等设计；最近著名的科学可视化软件Tecplot360已加入了三维数据网格的结合特征抽取功能，标志着流场拓扑结构特征的抽取盒可视化已经达到实用化阶段。 ¨ 时变流场的跟踪技术，特征结构可视化主要用于时变流场，一般是寻找、跟踪涡流、涡度等值面、激波等结构特征，然后进行事件检测和结果的可视化显示。Silver和Zabusky等提出了识别、量化和跟踪的方法，采用图像处理、形态学方法从标量、矢量场中提取不定型区域，用质量、重心、极值、体积、温度、压力等对特征量化，并将特征状态划分为产生、延续、消亡、分解和合并，并依据一定的匹配准则在相邻时间片之间建立特征对应关系，实现特征的跟踪与显示，K. L. Ma等针对非结构化网格上三维激波的探测提出了基于Mach数、基于密度梯度及基于方向导数的方法并比较了他们的优缺点和适用范围。 另外，还可选取用户感兴趣的作为特征，Walsum和Post提出了一种自动抽取用户感兴趣区域的技术，先用逻辑表达方式将用户感兴趣的区域定义出来，然后选择出那些满足用户要求的网格点。 自90年代以来，特征可视化在三个方面取得了进展：第一，Helman和Hesselink提出的一阶奇点(简单临界点)识别方法被扩展到了高阶奇点，实现了非线性向量场拓扑的可视化。这一进展把自动抽取特征的范围扩大到了场中较精细的结构，如缓慢弱小的漩涡等，并可定量地描述特征属性，为发展基于选择的特征可视化创造了条件。第二，提出了通过对标量场做向量化处理，增强及抽取标量场拓扑结构的方法。这一进展改变了此前多利用构建等值面提取标量场特征所带来的不够精确和难以选择的缺点，将标量场的特征抽取纳入到向量场特征可视化的统一处理构架。第三，提出了时变可视化的概念和处理时变场可视化的方法，其中Silver和Zabusky给出了时变场可视化的通用框架和算法，并将特征抽取和量化作为跟踪显示对象演变过程的主要部分。这表明带有属性量化功能的特征可视化已可以应用于四维的时变场，并成为当前处理时变场可视化的主要方法。 特征法从大量数据集中抽取仅仅包含特定信息的可视化方法，利用高度抽象的实体描述数据的特征，典型方法以拓扑结构分析法和针对时变场的特定结构抽取为主。单纯的拓扑结构分析法虽然可以把握流场中的本质信息，不过由于高度抽象，往往需要和其他映射方法结合使用 [ 8 ] 。