一种结合多物理属性的动态纹理瀑布建模方法

1.一种结合多物理属性的动态纹理瀑布建模方法，其特征在于由以下步骤组成：
1)在不同物理条件下采集各种瀑布的实地物理属性以及相应纹理数据；
2)预处理视频纹理数据，提取瀑布纹理控制参数；
3)将海量瀑布物理属性、纹理控制参数用数据库组织存储；
4)分析各种物理属性和每个纹理控制参数的统计分析；
5)根据订制的瀑布场景从数据库中提取相应纹理，遵循分布规律计算生成正确动态纹理；
所述步骤2)对视频纹理的预处理和控制参数的提取方法如下：
2.1)对视频纹理进行去噪、校正、缩放至固定分辨率大小；
2.2)对视频进行目标分类计算，结合地形数据分析得出瀑布的流向线；
2.3)沿着流向线等距离采样得到一个纹理样本矩阵M，其元素mij代表第j帧的第i个位置对应的纹理元，所述纹理元为一个32*32的像素矩阵；
2.4)从颜色光滑和宽度两方面简化纹理样本矩阵，得到一个纹理样本向量 其元素ni表示第i个位置的纹理元；
2.5)以日照参数为标准将视频分为早上，中午，晚上三组分别进行以上的组织，每个瀑布场景分别得到三组纹理样本向量；
所述步骤4)分析各种物理属性和每个纹理控制参数的统计分析如下：
设多物理属性集合为P；其中日照情况为s，流量为f，流速为v，落差为d；设整体纹理为T，其中纹理元为m，动态随机描述r，流线为l，亮度为b，纹理元数目n；提取纹理元和落差、流速相关性描述函数：m＝fm(d，v)；亮度与日照、流量的相关性描述函数：bfb(s，f)，纹理单元数目n＝fn(f)；最终瀑布纹理描述函数T＝FT(m，b，r，l，n)；
对控制参数的统计映射即是求解得到T＝FT(m，b，r，l，n)的映射。
2.如权利要求1所述结合多物理属性的动态纹理瀑布建模方法，其特征在于：所述步骤1)对瀑布的实地采集具体为：
1.1)对不同形态的多个瀑布进行有针对性的数据采集工作；
1.2)单个瀑布实地测量时需要在一年中不同季节、一天中不同日照情况分别进行采集；
1.3)每一次独立的采集工作都要同时记录当时的时间、太阳高度角、瀑布上游水流量、流速，以及使用高速摄像机多角度拍摄瀑布视频纹理。
3.如权利要求1所述结合多物理属性的动态纹理瀑布建模方法，其特征在于：所述函数T＝FT(m，b，r，l，n)估计方法如下：
4.1)在T＝FT(m，b，r，l，n)函数中，由于r是给合成的纹理加入时变信息，l是纹理的一个走向，b用于调整纹理亮度，n是使用粒子系统渲染时粒子的个数；四者都是对其物理状态有明确表达式约束的；而核心在于纹理元m和物理属性的映射；
4.2)对一组真实瀑布物理属性和相应的纹理进行一一对应；构建一个高维的散点图，因变量为m，自变量为f，v和d；
4.3)在三维空间中基于统计的使用高斯拟合一个包含每个对应散点的连续函数，即为物理属性和纹理特征对应的映射函数；
4.4)评价步骤4.3)中物理属性和纹理特征的映射函数的各个区间的可信度；
4.5)对计算的可信度通过指定阈值将所有区间划分为可靠区间和非可靠区间。
4.如权利要求3中所述结合多物理属性的动态纹理瀑布建模方法，其特征在于：所述步骤5)对根据订制的瀑布场景从数据库中提取相应纹理，遵循分布规律计算生成正确动态纹理其特征在于：
5.1)将需要建立的瀑布模型物理属性数据作为查询输入，从已组织好的数据库中取出与之最相邻的真实瀑布采样数据；
5.2)使用步骤4.3)中通过散点图拟合物理属性和纹理特征得到的物理属性和纹理数据的映射关系计算出该瀑布模型中的动态纹理和相邻纹理的欧式距离；
5.3)使用与真实纹理数据的距离加权计算出给定情况下的纹理状态；
5.4)将计算出的动态纹理映射到瀑布模型的三维表面上，形成瀑布模型。
5.如权利要求4中所述结合多物理属性的动态纹理瀑布建模方法，其特征在于：所述步骤1)还包括对真实瀑布采样数据的选择具体如下：
根据所述步骤4.5)对物理属性和纹理特征拟合映射函数所划分的可靠区间和非可靠区间，确定订制瀑布场景的物理属性组合所属区间是否可靠；若处于可靠区间则使用欧氏距离最近的4个样本值进行差值计算，否则计算最邻近的8个样本的差值结果。

一种结合多物理属性的动态纹理瀑布建模方法\n技术领域\n[0001] 本发明属于虚拟现实技术领域，具体来说是虚拟现实技术中物体真实感仿真中动态纹理的一种建模方法。\n背景技术\n[0002] 虚拟现实所需要模拟千差万别的目标，通常要求输入和输出图像序列具有稳定的内容和一定程度的运动形态重复特征。对于瀑布场景这样动态性很高的仿真一般由于模型表示能力的局限性，在表达运动形态的目标时，动态纹理的真实效果和可信程度至关重要。\n[0003] 使用传统计算机图形学的三角面片的方法来模拟瀑布场景计算量很大而且效果往往不理想。目前对瀑布仿真的研究主要集中在以下两种途径：基于粒子系统的模拟方法；\n图像或视频纹理合成的瀑布仿真方法。\n[0004] 基于粒子系统的瀑布模拟：1983年，Reeves首先提出了粒子系统。粒子系统用非常简单的体素来构造复杂的物体，充分体现了不规则模糊物体的动态性和随机性。由于粒子系统的模拟是按照一定的约束方程对每个粒子分别进行位置、速度、性质以及生存状态的计算，这一特性使得其在瀑布建模上克服了传统三角面片对不规则物体仿真的运动失真。但是粒子系统中每个粒子都使用统一的纹理使得其表达瀑布的真实感上存在缺陷。\n[0005] 利用图像或视频纹理合成的瀑布模拟：1999年的ICCV会议上，Efros和Leung提出了基于MRF模型的点匹配纹理合成，简单有效，能很好的合成一些随机纹理。Arno Schodl[1]在2000年的SIGGRAPH年会上提出了一种新的媒体——视频纹理。视频纹理同时具有图像和视频的特征：时序稳定性和内容动态性。利用视频纹理，可以实现利用有限样本表达无限信息的目标。以上特性使得其在瀑布建模上能表现出瀑布表面的视觉真实感。\n但缺点在于由于这些合成的纹理仅仅具有二维属性，使得模拟的瀑布场景不能进行三维漫游。\n[0006] 可以把关于动态纹理的研究方法分为三类：第一类方法把动态纹理当作时间序列分析，应用ARMA模型或者ARMA模型的扩展形式进行研究。第二类方法是基于纹理合成的扩展方法，由于其原本只是针对纹理图像设计的。第三类方法如视频纹理方法，它介于基于物理模型的仿真方法和时间序列方法之间。\n[0007] 时间序列分析主要指采用参数模型(特指ARMA模型)对所观测的有序的随机数据进行分析与处理的一种数据处理、系统辨识和系统分析方法。Szummer和Picard提出的时间纹理模型采用时-空自回归(STAR)模型表示动态纹理，STAR模型是自回归(AR)模型的三维扩展形式。该方法对非线性、不稳定动态纹理和受控合成动态纹理的想过任然不够理想。\n[0008] 纹理合成是计算机图形学和计算机视觉中相当活跃的研究领域，应用纹理合成方法的扩展形式合成。Bar-Jaseph提出的时变纹理方法是一种基于统计学习的通用的纹理信号建模方法。该方法由于其原本只是针对纹理图像设计的，所以也局限于处理由纹理图像序列所组成动态纹理。\n[0009] Schodl等人提出的视频纹理方法的思想是：给定一段动态纹理视频，只要视频中存在足够的相似帧，在视频的播放过程中，就可以不断的从跳变点处往回切换，合成任意长度的视频。由于视频播放次序的重组总是要求切换帧之间相似，然而许多动态纹理不符合这一要求，这是第三类方法难以跨越的障碍。\n[0010] 近期研究者结合粒子系统和动态纹理的方法对瀑布进行了具有真实感、动态性的模拟，但是在纹理的合成上没有考虑瀑布自身各项物理属性的对其纹理的影响。所以存在一个缺点：对一个没有真实纹理采集的订制瀑布场景，无法合成符合其真实物理特性的动态纹理，以至于渲染结果的可信度较差。\n[0011] 总的来说，在瀑布的建模上粒子系统和动态纹理的结合能进行动态性和真实感较强的仿真。而动态纹理合成技术的研究工作尚未成熟，没有考虑到多种物理属性对纹理特性的影响，无法真实可信的估计出位置物理状态下的目标纹理形态特征。而在虚拟现实对瀑布的模拟仿真中可信度有着有至关重要的作用。\n发明内容\n[0012] 本发明解决的技术问题是：克服了目前已有的所有瀑布动态纹理生成方法中都忽略瀑布本身的多种物理属性对其形成的纹理特征有着一定相关影响这一缺点。为已有的瀑布纹理建模方法在仿真一个订制瀑布场景时使用的纹理可信度较低的问题提出了一个可行的解决方案。\n[0013] 本发明提出一种综合考虑水流流速、流量、落差、地形等物理因素对瀑布动态纹理的影响，进而建立与特定场景相对应的真实感瀑布模型的方法。包括以下步骤：\n[0014] 1)在不同物理条件下采集各种瀑布的实地物理属性以及相应纹理数据；\n[0015] 2)预处理视频纹理数据，提取瀑布纹理控制参数；\n[0016] 3)将海量瀑布物理属性、纹理控制参数用数据库组织存储；\n[0017] 4)分析各种物理属性和每个纹理控制参数的统计分析；\n[0018] 5)根据订制的瀑布场景从数据库中提取相应纹理，遵循分布规律计算生成正确动态纹理；\n[0019] 所述步骤2)对视频纹理的预处理和控制参数的提取方法如下：\n[0020] 2.1)对视频纹理进行去噪、校正、缩放至固定分辨率大小；\n[0021] 2.2)对视频进行目标分类计算，结合地形数据分析得出瀑布的流向线；\n[0022] 2.3)沿着流向线等距离采样得到一个纹理样本矩阵M，其元素mij代表第j帧的第i个位置对应的纹理元，所述纹理元为一个32*32的像素矩阵；\n[0023] 2.4)从颜色光滑和宽度两方面简化纹理样本矩阵，得到一个纹理样本向量N，其元素ni表示第i个位置的纹理元；\n[0024] 2.5)以日照参数为标准将视频分为早上，中午，晚上三组分别进行以上的组织，[0025] 每个瀑布场景分别得到三组纹理样本向量；\n[0026] 所述步骤4)分析各种物理属性和每个纹理控制参数的统计分析如下：\n[0027] 设多物理属性集合为P；其中日照情况为s，流量为f，流速为v，落差为d；设整体纹理为T，其中纹理元为m，动态随机描述r，流线为l，亮度为b，纹理元数目n；提取纹理元和落差、流速相关性描述函数：m＝fm(d，v)；亮度与日照、流量的相关性描述函数：b＝fb(s，f)，纹理单元数目n＝fn(f)；最终瀑布纹理描述函数T＝FT(m，b，r，l，n)；\n[0028] 对控制参数的统计映射即是求解得到T＝FT(m，b，r，l，n)的映射。\n[0029] 所述步骤1)对瀑布的实地采集具体为：\n[0030] 1.1)对不同形态的多个瀑布进行有针对性的数据采集工作；\n[0031] 1.2)单个瀑布实地测量时需要在一年中不同季节、一天中不同日照情况分别进行采集；\n[0032] 1.3)每一次独立的采集工作都要同时记录当时的时间、太阳高度角、瀑布上游水流量、流速，以及使用高速摄像机多角度拍摄瀑布视频纹理。\n[0033] 所述的T＝FT(m，b，r，l，n)映射的计算，其特征在于：所述函数T＝FT(m，b，r，l，n)估计方法如下：\n[0034] 4.1)在T＝FT(m，b，r，l，n)函数中，由于r是给合成的纹理加入时变信息，l是纹理的一个走向，b用于调整纹理亮度，n是使用粒子系统渲染时粒子的个数；四者都是对其物理状态有明确表达式约束的；而核心在于纹理元m和物理属性的映射；\n[0035] 4.2)对一组真实瀑布物理属性和相应的纹理进行一一对应；构建一个高维的散点图，因变量为m，自变量为f，v和d；\n[0036] 4.3)在三维空间中基于统计的使用高斯拟合一个包含每个对应散点的连续函数，即为物理属性和纹理特征对应的映射函数；\n[0037] 4.4)评价步骤4.3)中物理属性和纹理特征的映射函数的各个区间的可信度；\n[0038] 4.5)对计算的可信度通过指定阈值将所有区间划分为可靠区间和非可靠区间。\n[0039] 所述步骤5)对根据订制的瀑布场景从数据库中提取相应纹理，遵循分布规律计算生成正确动态纹理其特征在于：\n[0040] 5.1)将需要建立的瀑布模型物理属性数据作为查询输入，从已组织好的数据库中取出与之最相邻的真实瀑布采样数据；\n[0041] 5.2)使用步骤4.3)中通过散点图拟合物理属性和纹理特征得到的物理属性和纹理数据的映射关系计算出该瀑布模型中的动态纹理和相邻纹理的欧式距离；\n[0042] 5.3)使用与真实纹理数据的距离加权计算出给定情况下的纹理状态；\n[0043] 5.4)将计算出的动态纹理映射到瀑布模型的三维表面上，形成瀑布模型。\n[0044] 所述步骤1)中对真实瀑布采样数据的选择具体如下：\n[0045] 根据所述步骤4.5)对物理属性和纹理特征拟合映射函数所划分的可靠区间和非可靠区间，确定订制瀑布场景的物理属性组合所属区间是否可靠；若处于可靠区间则使用欧氏距离最近的4个样本值进行差值计算，否则计算最邻近的8个样本的差值结果。\n[0046] 本发明的有益效果是：\n[0047] 1)将瀑布的多种物理属性和其表现的纹理特征的相关性考虑入计算机仿真的建模方法中。克服了传统方法中忽略这一相关性而导致的对虚拟场景中瀑布外观纹理的生成失真和模型不可信问题。\n[0048] 2)在本发明中不是将纹理图片作为基本的映射单元，而是将纹理进一步分解为纹理元，同时提取出纹理控制参数，能更加有效的控制动态纹理的形成。\n[0049] 3)本发明将与瀑布纹理联系最紧密的流量、流速、落差、日照情况这几个典型物理属性作为影响纹理特征的改变量，并使用海量的各种条件下的真实瀑布参数作为样本来估计各物理属性和纹理的映射。使得对瀑布动态纹理的生成具有更加真实可信的效果。\n附图说明\n[0050] 图1是本发明中对整体瀑布动态纹理计算的方法流程图；\n[0051] 图2是本发明中对真实瀑布采集数据库建立示意图；\n[0052] 图3是本发明中纹理控制参数提取流程图；\n[0053] 图4是本发明中物理属性和纹理元映射关系提取示意图；\n[0054] 图5是本发明中使用粒子系统渲染瀑布场景流程图。\n具体实施方式\n[0055] 以下，参照附图，对本发明的具体过程作说明，但本发明不限于图事例。\n[0056] 如图1所示的完整动态纹理映射计算流程，通过参数入库，特征描述，统计映射，生成动态纹理，模拟仿真这个流程来使用瀑布多种物理属性计算出真实可信的动态纹理[0057] 步骤1为了捕获各种物理属性下的真实瀑布样本，需要在一年中各个季节，每天的不同日照情况下对大量的瀑布进行实地采集。从多个角度对一个瀑布的外观进行视频纹理采集。获取其落差、流速、流量、日照、地形情况等多种物理属性的途径可以使和视频纹理同时用相应的一起进行实地测量，也可以是从地质、气象部门采购以上瀑布状态属性数据。\n以图2所示的方式建立采集数据库。\n[0058] 步骤2对视频纹理进行预处理后，提取出纹理控制参数：亮度b、动态性r、纹理元m、流向l，单位纹理数目n。纹理元m由纹理样本向量来描述，提取过程如图3所示。\n[0059] 步骤2.1预处理视频纹理，进行图像去噪并缩放至同样的分辨率大小。\n[0060] 步骤2.2对视频逐帧使用目标提取算法得到视频中水的部分；再进行非监督分类，每一类即是一个瀑布的流向线。利用地形数据中的瀑布表面几何形状对流向线进行检验。\n[0061] 步骤2.3沿流线对纹理等距离采样10组32*32象素的纹理元。对每一帧进行类似处理，一个瀑布场景视频的得到一个纹理样本矩阵M，其元素mij代表第j帧的第i个位置对应的纹理元。\n[0062] 步骤2.4简化样本矩阵M，获得纹理样本向量 其元素ni表示第i个位置的纹理元特征。\n[0063] 步骤2.4.1计算纹理元中颜色期望：对矩阵M每行分别求纹理元中的颜色期望值记录在向量Ui中。用以描述各帧相应纹理元形状的变化。\n[0064] 步骤2.4.2计算纹理元中颜色的平滑度：引入光滑项V(i，j)：\n[0065] j＝1，2，...，N\n[0066] 其中(s，t)表示一个像素，N为矩阵的总帧数。Ui为步骤2.4.1中颜色期望向量，K为纹理像素的尺寸本例中为32。以此描述各帧相应纹理元颜色空间差异。\n[0067] 步骤2.4.3得到纹理样本向量：使用第m(i，j*)元素为纹理样本向量该位置的值，*\n其中j 求满足：\n[0068] \n[0069] 步骤2.4.4对每个M矩阵每行进行上面三步计算之后即得到描述瀑布纹理的 向量。\n[0070] 步骤3定义物理属性和其他纹理控制参数的约束函数。\n[0071] 设多物理属性集合为P；其中日照情况为s，流量为f，流速为v，落差为d。设整体纹理为T，动态随机描述r，流线为l，亮度为b，纹理数目n。亮度与日照、流量的相关性描述函数：b＝fb(s)，纹理单元数目n＝fn(f)。\n[0072] 其中对于亮度b＝fb(s)和纹理单元数目n＝fn(f)，由其定义可知各自均为与特定一项物理属性相关，为一个简单的二维约束关系。其映射可以通过将对应的样本点使用简单的线性插值的方式来获得新的物理属性值对应的控制参数值。\n[0073] 而对于动态性参数r的描述：由于动态性是瀑布内在的变化情况与外在物理属性无关。可以通过在时间维上合理的改变纹理的像素来给瀑布增加真实的动态效果。使用Hermite差值方法来实现，具体定义如下：\n[0074] \n[0075] \n[0076] t表示时间ti和ti+1分别表示第i和第i+1帧时刻。c为颜色ci和ci+1分别表示第i和i+1帧的像素颜色值，可以对RGB三通道分别进行同样的计算。d′i和d′i+1分别表示在第i和第i+1帧流线上两位置之间的切方向。N为差值矩阵。用fr(r，c)表示动态性和颜色的关联函数。\n[0077] 流线l与订制瀑布场景有关，可以是交互式的由用户完成流线的指定，也可以根据瀑布场景的几何表面凹凸情况可以使用简单流体运动方法来确定。\n[0078] 步骤4统计分析纹理特征与各种物理属性的映射关系。\n[0079] 用纹理元m参数来表示纹理特征，在步骤2中每个瀑布使用了一个纹理样本矩阵来组织纹理元m。物理属性流速和落差和纹理特征的映射关系使用函数：m＝fm(d，v)来表示。由该函数的定义可知其为一个二元函数，可以使用高斯逼近的方式来进行函数曲面拟合。示意如图4所示。具体方法是将数据库中已经存储的真实样本值的每组纹理元m和落差d、流速v像对应，在三维空间使用三维高斯拟合公式迭代的计算出分布函数，如图4曲面所示。\n[0080] 步骤5以各参数和相关物理属性的映射关系为基础，对订制的瀑布场景从数据库中提取相应纹理，遵循映射关系计算生成符合需求场景的瀑布动态纹理。\n[0081] 步骤5.1确定订制的场景的各项参数，以步骤3定义的方法通过日照和流量情况确定亮度、和纹理数目值。\n[0082] 步骤5.2进行区间可信度评价：根据映射中样本的分布情况来评价该映射的每个区间的可信度。将样本分布的密集程度作为可信度标准标准，密度越高该区间可信度越高，反之亦然。\n[0083] 步骤5.3对纹理元m和流速v，落差d的映射函数中，依照可信度分为可靠区间和不可靠区间。\n[0084] 在订制瀑布场景中新非样本物理属性(d，v的任意组合)所处区间若为可靠区间，即可找到离该物理属性组合最近的4个邻域样本点，选取符合该分布函数上特征的权值。\n最后按确定的权值以样本点差值生成新的纹理元。若是不可靠区间则找最邻近8个样本点进行上述计算，以提高真实性。\n[0085] 步骤6用粒子系统在瀑布三维场景中渲染，完成对瀑布的模拟仿真。渲染流程如图5所示。\n[0086] 步骤6.1用纹理数目n做为粒子系统的初始粒子数目。\n[0087] 步骤6.2粒子沿流线方向运动，每个粒子根据所处离线的位置选取纹理样本向量中相应位置上的纹理元做纹理映射。\n[0088] 步骤6.3使用步骤3中定义的动态函数fr(r，c)进行对瀑布纹理时变性质的渲染控制。\n[0089] 最后应说明的是，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的人员来说，在不脱离结合多种物理属性对纹理特征的影响来进行建模仿真的前提下，还可以做出若干改进或等同替换，这些改进和替换也应视为本发明的保护范围。