CN104517299B - 视频流体物理驱动模型恢复及重新仿真的方法 - Google Patents
视频流体物理驱动模型恢复及重新仿真的方法 Download PDFInfo
- Publication number
- CN104517299B CN104517299B CN201410805222.0A CN201410805222A CN104517299B CN 104517299 B CN104517299 B CN 104517299B CN 201410805222 A CN201410805222 A CN 201410805222A CN 104517299 B CN104517299 B CN 104517299B
- Authority
- CN
- China
- Prior art keywords
- fluid
- density
- normalized
- geometry
- rho
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Abstract
本发明公开了一种视频流体物理驱动模型恢复及重新仿真的方法,该方法旨在利用视频流的样例帧,得到反向工程的流体仿真结果。该方法首先采用HS(Horn‑Schunck)方法,利用多尺度框架,在每层对第一帧流体粒子的二维运动速度进行计算,然后在每层中利用浅水波方程做约束,得到了较为准确的流体运动速度及非归一化几何。在流体仿真过程中,利用样例帧表面的几何作为约束,调整流体粒子运动的密度属性,从而得到了流体仿真的物理驱动模型。最后利用增强的粒子密度进行重仿真,得到了具有真实感细节的重仿真结果。本发明可有效地应用于流体的反向工程的研究及应用中,可以定制任意场景,克服了基于物理的传统流体仿真中计算量大以及实时性差的问题。
Description
技术领域
本发明涉及视频流体反向工程中物理属性及驱动模型的恢复。利用运动矢量计算的多尺度框架及浅水波方程优化得到流体表面的非归一化几何和初始帧的二维运动速度,恢复流体运动属性,得到反向工程中的驱动模型,获得具有真实感细节的重仿真结果。
背景技术
真实感流体仿真目前已经广泛地被应用于计算机游戏、动画、电影特效、军事仿真等领域,关于流体仿真计算的研究已经取得了很多的成果。限于流体物理仿真中的大量物理属性计算的问题,近些年人们开始从事于基于视频的流体仿真技术的研究。怎样充分利用流体视频中的丰富真实感信息及流体的运动信息,进一步生成实时的、具有真实感特性的流体仿真场景,这就是视频流体反向工程中亟待解决的问题。
近年来,一些研究工作已经开始从事基于视频的建模及绘制的研究,一些研究中,采用能量优化的策略,并权衡计算量及外观的特性,从而使得生成的流体仿真场景不具有流体真实的物理特性;在另外的相关研究中,人们仅仅恢复流体表面的归一化几何,由于这些几何信息仅能反映流体表面的外观特点,还不能用于对于流体表面的真实物理属性的恢复。
由于流体运动的特殊性,对于流体反向工程的研究存在着一定难度。目前人们在研究中利用流体运动速度跟踪的方法来实现流体重仿真的过程,存在的主要问题是流体运动的物理驱动模型并没有真正得到恢复;在另外的反向工程相关研究中,人们试图采用在低维空间中进行流体仿真,为了保持流体仿真的细节,如何从低维空间建立简化的仿真模型,这已经成为研究中瓶颈,尽管人们已经做了一些努力探索这个问题,目前如何在流体方向工程的低维空间中,快速生成与原流体场景相似结果以及长序列的流体重仿真结果,仍然是一个亟待解决的实际问题。
发明内容
本发明的目的是针对现有技术的不足而提供的一种视频流体反向工程中物理属性及驱动模型恢复的方法,利用该方法能够得到视频流体反向工程中流体物理属性及驱动模型,利用该模型可以得到流体重仿真的结果。
实现本发明目的的具体技术方案是:
(1)利用多尺度的框架,计算视频流体第一帧的表面几何和流体二维运动速度
为了计算视频流体的表面几何和流体粒子二维运动速度,采用以下两个条件约束:流体表面的外观的约束和流体运动的动态行为约束。
首先,利用多尺度光流计算的框架,采用5层尺度规模,初始化流体表面粒子的二维运动速度。利用视频流体初始两帧的灰度图像I0和I1产生5层不同尺度的图像其中s={0,1,2,3,4},及对应着较粗糙的模型,而及对应着较精细的模型。将每层流体粒子运动速度记为t表示帧的序号,s表示尺度规模,s={0,1,2,3,4}。从第0层开始,每一层得到的二维运动速度传给下一层。每层的计算按照以下三个步骤进行:
(a)二维运动速度的初始化
利用现有多尺度光流计算方法,计算流体第一帧的二维速度
(b)非归一化几何的计算
在每层中,利用现有方法计算流体表面的归一化几何,表示为z;非归一化几何表示为Z,比例系数(我们也称之为非归一化因子)假设为S,那么,
Z=Sz (1)
利用浅水波方程,根据现有的研究结果中速度与高度之间的计算关系
将二维运动速度的初始化的结果看作是较准确的,因此可以得到:
其中,
其中,ux、uy表示运动速度第一个分量u在x轴和y轴方向的导数,vx、vy表示运动速度第二个分量v在x轴和y轴方向的导数,uxy表示运动速度第一个分量u在x轴和y轴方向的混合偏导数,vxy表示运动速度第二个分量v在x轴和y轴方向的混合偏导数。Zx表示Z在x轴方向的导数,Zy表示Z在y轴方向的导数,Zt表示Z在时序上的导数。h表示相邻像素之间的坐标差,取1值,α是平滑系数,取7.0。
将式(1)带入到方程式(2)中,可以得到有关S的四次方程:
P1S4+P2S2+P3=0 (13)
利用t=S2对上式进行换元,得到有关t的一元二次方程:
P1t2+P2t+P3=0 (14)
其中,
P1=CDu+DFv-ABu-BEv (15)
根据式(14)求得t,进一步能够得到S;在求取t和S时,一元二次方程有两个解,舍去不合理的一个解;
(c)流体表面几何和速度的优化
按照以下步骤循环两次进行优化:步骤1,将步骤(b)得到的第k次非归一化几何看作是较准确的(第一次迭代假设k=0),利用以下公式(24)和(25)可以得到第k+1次更准确的速度(uk+1,vk+1);步骤2,将步骤1得到的速度看作是较准确的结果,再利用方程(13)计算得到非归一化因子S,从而得到更准确的表面几何;
每个粒子利用上述第L-1层的步骤(c)计算得到流体表面几何的非归一化因子,以下步骤是利用粒子的非归一化因子来计算整个场景的S:如果任意两个粒子i和j计算得到的非归一化因子分别为Si和Sj,聚类时,如果Si和Sj满足以下条件,那么Si和Sj就标记为同一类:
|Si-Sj|<ST i,j∈{0,1,2,3…N-1} (26)
其中N为粒子总数目,ST为S的阈值,设其值为1.0.
聚类时类的数量是自动得到的,如果聚得最大类的样本数为n,最大类中所有粒子求得的非归一化因子的和为Ss,那么,
S=Ss/n (27)
(2)利用步骤(1)得到流体表面几何的非归一化因子S,再计算视频样例帧非归一化几何取视频中前100帧作为样例帧,按照步骤(1)得到流体表面几何的非归一化因子S,再利用式(1)计算这些样例帧表面非归一化几何;
(3)采用LBM(Lattice Boltzmann Method),恢复流体运动的属性及驱动模型
对于所有样例帧,利用LBM进行仿真时,利用式(28)和(29)计算并调整粒子的粒子密度和分布函数;
其中,hL是粒子演化的高度,fi(t+1)和fi(t)分别是粒子在i方向的t+1时刻和t时刻的分布函数,i表示LBM的D2Q9模型中第i个方向,i=0,1,2…8,ρ(t+1)和ρ(t)分别是粒子在t+1时刻和t时刻的密度;
(4)对流体粒子仿真密度属性进行增强
将步骤(3)中得到的粒子密度我们称其为标准密度,记为ρs;LBM自动推演产生的密度,称之为自动演化密度,记为ρa;分析ρs的直方图的分布,并对ρa进行增强,增强后的密度称之为增强密度,记为ρe;
对于ρs,求其最大值和最小值将区间等分成256个子区间:
[iss,(i+1)ss] i={0,1...255} (30)
其中,ss为区间尺度:
标准密度ρs表示为ρs=Ss×ss+Os,将ρs记为[Ss,Os],其中Ss为区间号,Os为区间内密度位移量,其满足:
对ρa,求其最大值和最小值将区间等分成256个子区间:
[isa,(i+1)sa] i={0,1...255} (33)
其中,sa为区间尺度:
对于任何一个自动演化密度ρa,都可以表示为ρa=Sa×sa+Oa,将ρa记为[Sa,Oa]其中Sa为区间号,Oa为区间内密度位移量,其满足:
根据ρs的直方图的统计分布,利用图像增强理论中直方图规定化的方法,实现Ss和Sa的对应,如果对于任意一个自动演化密度它的区间号为如果和它匹配的区间号为它被增强后的密度ρe可以计算为:
(5)重仿真过程
视频流体的重仿真过程包括以下步骤:
(a)利用步骤(1)初始化流体初始帧的运动速度及非归一化几何
(b)利用步骤(2)计算样例帧表面非归一化几何
(c)利用初始帧运动速度,推演LBM,同时利用步骤(3)和(4)增强流体粒子的密度
利用增强的密度演化仿真,得到具有细节特征的真实感流体重仿真结果。
本发明具有简单、实用的特点,能够实现流体物理属性及驱动模型的恢复,有效地实现具有细节特征的真实感流体的重仿真过程。克服了现有流体仿真的反向工程中,缺少物理模型驱动的问题。进一步实例证明了该发明方法的有效性。可以实现流体场景的放大、缩小、定制功能,并且所恢复的物理属性及驱动模型能够满足边界条件。
附图说明
图1为本发明实施例流体粒子前两帧之间运动速度初始化结果图;
图2为本发明实施例流体粒子表面非归一化几何的计算结果图;
图3为本发明实施例流体粒子的密度计算可视化的结果图;
图4为本发明实施例流体重仿真的结果图;
图5为本发明实施例中恢复流体物理属性的边界条件图;
图6为本发明实施例中烟的重仿真结果,及与现有方法的比较图;
图7为本发明实施例反向工程中场景定制尺度的应用实例图;
图8为本发明实施例反向工程中定制场景的应用实例图。
具体实施方式
实施例
下面结合附图对于发明进一步说明。
本实施例采用DynTex动态纹理库中的流体视频进行物理属性及驱动模型的恢复。在PC机上Windows7操作系统下进行实施,其硬件配置是2.66GHz Intel Core(TM)2 DuoCPU、4GB RAM。(1)利用多尺度的框架,计算视频流体前两帧的表面几何和流体粒子二维运动速度
(a)二维运动速度的初始化
利用现有多尺度光流计算方法,采用5层尺度规模,计算流体第一帧的二维速度(u,v).
(b)非归一化几何的计算
在每层中,利用现有方法计算流体表面的归一化几何,对归一化几何我们表示为z;非归一化几何我们表示为Z,比例系数假设为S,那么,
Z=Sz(1)
利用浅水波方程,根据现有的研究结果中速度与高度之间的计算关系,并将二维运动速度的初始化的结果看作是较准确的,因此可以得到速度的计算关系:
其中,
其中,ux、uy表示运动速度第一个分量u在x轴和y轴方向的导数,vx、vy表示运动速度第二个分量v在x轴和y轴方向的导数,uxy表示运动速度第一个分量u在x轴和y轴方向的混合偏导数,vxy表示运动速度第二个分量v在x轴和y轴方向的混合偏导数;Zx表示Z在x轴方向的导数,Zy表示Z在y轴方向的导数,Zt在时序上的导数;h表示相邻像素之间的坐标差,取1值,α是平滑系数,取7.0;
将式(1)带入到浅水波方程式(2)中,可以得到有关S的四次方程:
P1S4+P2S2+P3=0 (13)
我们利用t=S2对上式进行换元,得到有关t的一元二次方程:
Plt2+P2t+P3=0 (14)
其中,
P1=CDu+DFv-ABu-BEv (15)
根据式(14)可以求得t,进一步可以得到S。
(c)流体表面几何和速度的优化
按照以下步骤循环两次进行优化:步骤1,将步骤(b)得到的第k次非归一化几何看作是较准确的(第一次迭代假设k=0),利用以下公式(24)和(25)可以得到第k+1次更准确的速度(uk+1,vk+1);步骤2,将步骤1得到的速度看作是较准确的结果,再利用方程(13)计算得到非归一化因子S,从而得到更准确的表面几何。
每个粒子利用上述第L-1层的步骤(c)计算得到流体表面几何的非归一化因子S,利用以下步骤是利用粒子的非归一化因子来计算整个场景的S;如果任意两个粒子i和j计算得到的S表示为Si和Sj,聚类时,如果Si和Sj满足以下条件,那么Si和Sj就标记为同一类:
|Si-Sj|<ST i,j∈{0,1,2,3…N-1} (26)
其中N为粒子总数目,ST为S的阈值,我们设其值为1.0.
聚类时类的数量是自动得到的,如果聚得最大类的样本数为n,最大类中所有粒子求得的非归一化因子的和为Ss,那么,
S=Ss/n (27)
表1是不同视频计算得到的S的结果。
表1 不同视频计算得到S的结果
(2)利用步骤(1)得到流体表面几何的非归一化因子S,再计算视频样例帧非归一化几何
取视频中前100帧作为样例帧,按照步骤(1)得到流体表面几何的非归一化因子S,再利用式(1)计算这些样例帧表面非归一化几何;
图2是计算得到的流体表面非归一化几何信息。
(3)采用LBM(Lattice Boltzmann Method),恢复流体运动的属性及驱动模型
对于所有样例帧,利用LBM进行仿真时,利用式(28)和(29)计算并调整粒子的粒子密度和分布函数;
其中,hL是粒子演化的高度,fi(t+1)和fi(t)分别是粒子在i方向的t+1时刻和t时刻的分布函数,i表示LBM的D2Q9模型中第i个方向,i=0,1,2…8,ρ(t+1)和ρ(t)分别是粒子在t+1时刻和t时刻的密度;
(4)对流体粒子仿真密度属性进行增强
将步骤(3)中得到的粒子密度我们称其为标准密度,记为ρs;LBM自动推演产生的密度,称之为自动演化密度,记为ρa;分析ρs的直方图的分布,并对ρa进行增强,增强后的密度称之为增强密度,记为ρe;
对于ρs,求其最大值和最小值将区间等分成256个子区间:
[iss,(i+1)ss] i={0,1...255} (30)
其中,ss为区间尺度:
标准密度ρs表示为ρs=Ss×ss+Os,将ρs记为[Ss,Os],其中Ss为区间号,Os为区间内密度位移量,其满足:
对ρa,求其最大值和最小值将区间等分成256个子区间:
[isa,(i+1)sa] i={0,1...255} (33)
其中,sa为区间尺度:
对于任何一个自动演化密度ρa,都可以表示为ρa=Sa×sa+Oa,将ρa记为[Sa,Oa]其中Sa为区间号,Oa为区间内密度位移量,其满足:
根据ρs的直方图的统计分布,利用图像增强理论中直方图规定化的方法,实现Ss和Sa的对应,如果对于任意一个自动演化密度它的区间号为如果和它匹配的区间号它被增强后的密度ρe可以计算为:
图3是视频流体标准密度ρs计算的可视化结果,从该结果可以看出该发明专利能够有效性地恢复物理驱动模型。
(5)重仿真过程
视频流体的重仿真过程包括以下四个步骤:
(a)利用步骤(1)初始化流体初始帧的运动速度及非归一化几何
(b)利用步骤(2)计算样例帧表面非归一化几何
(c)利用初始帧运动速度,推演LBM,同时利用步骤(3)和(4)增强流体粒子的密度
(d)利用增强的密度演化仿真,得到具有细节特征的真实感流体重仿真结果。
图4是重仿真的结果,每组中左上结果是原流体视频,右下结果是重仿真的结果。
利用上述的流体重仿真实施方法,进一步验证了边界条件。我们仍利用采用DynTex动态纹理库中的流体视频进行物理属性及驱动模型的恢复以后,再进行物理属性的边界条件测试,采用现有LBM的边界条件的控制方法,改变边界处的分布函数,得到了实施结果,证明边界条件的有效性。图5是边界条件的实施结果。
将本发明与现有方法进行比较:
为了验证本发明中,流体运动属性及物理驱动模型恢复的有效性,采用JamesGregson2014研究中烟的视频进行重仿真验证,并与现有的研究结果进行比较,图6是对烟视频进行重仿真的结果,左侧是原烟的视频,中间的是现有方法的结果,右面的结果是本发明的结果,从其可以证明流体物理属性及物理驱动模型的恢复也适用于烟流体的反向工程。与现有的基于运动跟踪的方法相比,本发明的仿真结果具有真实感强的特点,并且反向工程的结果与原视频更加相似。
为了验证该发明中物理属性恢复的准确性,与现有的方法进行比较,分别采用现有方法和本发明对流体运动速度进行初始化,然后进行驱动LBM仿真,计算演化高度和计算得到的表面几何之间的差异:
其中,n是粒子的数目,hi是粒子i的演化高度,Zi是粒子i的非归一化的表面几何。
表2列出了与现有方法进行比较的误差结果。
表2 与现有方法进行比较的误差结果
从表中的误差结果可以看出,本发明中恢复的运动速度更加准确一些。
在反向工程中应用的实施。
应用1:利用本发明定制场景规模。可以任意尺度的场景,实施方法:利用图像放大及缩小的理论,生成大场景(或小场景)的密度场,从而可以驱动流体物理模型进行驱动。具体地,如果场景尺度放大的倍数为2×3,采用图像放大理论中的后向映射法生成密度场,如果目标位置的坐标为(x,y),那么可以计算得到其对应的原场景的对应位置(x0,y0):
x0=x/2 y0=y/3 (38)
那么可以利用对应位置(x0,y0)处的密度生成新的坐标处(x,y)的密度。如果计算得到的结果非整型数值,可以采用最近邻插值法或双线性插值方法计算得到(x,y)点处的密度值。图7是利用该方法得到的不同放大倍数的场景结果,每组中,左上是原来尺度重仿真结果,右上是0.7×0.7的缩放比例,即水平方向尺度缩小到0.7倍,垂直方向也缩小到0.7倍。下面的放大比例为1.5×1.5。
应用2:利用本发明的方法,定制密度分布直方图,从而可以得到定制的流体场景。在实施中,我们改变密度直方图的分布,得到了新的流体定制场景。为了可视地表现出新的流通场景高度几何与外观的对应关系,我们生成了颜色与灰度之间对应的映射关系,如果灰度值为g,那么利用原来场景,生成对应关系:CR[g]、CG[g]和CB[g]。在流体仿真中可以得到定制场景的演化高度,再根据现有的基于SFS(Shape from Shading)的重建方法,可以计算得到粒子的灰度值,再利用所建立的颜色与灰度之间的映射关系,可以得到表面的颜色映射结果。图8是定制流体重仿真的结果。每组中左侧是流体重仿真的结果,右侧是定制的流体仿真结果。
Claims (1)
1.视频流体物理驱动模型恢复及重新仿真的方法,其特征在于该方法包括以下步骤:
(1)利用多尺度的框架,计算视频流体第一帧的表面几何和流体二维运动速度
按强度将图像共分L层,粗糙的层记为0层,精细的层为第L-1层,每层的计算按照以下三个步骤进行;从第0层开始,每一层得到的二维运动速度传给下一层;
(a)二维运动速度的初始化
利用现有多尺度光流计算方法,计算流体第一帧的二维速度(u,v);
(b)非归一化几何的计算
利用现有方法计算流体表面的归一化几何,表示为z;非归一化几何表示为Z,比例系数也称之为非归一化因子,假设为S,那么,
Z=Sz (1)
利用浅水波方程,根据现有的研究结果,得到速度与高度之间的计算关系:
其中
其中,ux、uy表示运动速度第一个分量u在x轴和y轴方向的导数,vx、vy表示运动速度第二个分量v在x轴和y轴方向的导数,uxy表示运动速度第一个分量u在x轴和y轴方向的混合偏导数,vxy表示运动速度第二个分量v在x轴和y轴方向的混合偏导数;Zx表示Z在x轴方向的导数,Zy表示Z在y轴方向的导数,Zt表示Z在时序上的导数;h表示相邻像素之间的坐标差,取1值,α是平滑系数,取7.0;
将式(1)带入到浅水波方程式(2)中,得到有关S的四次方程:
P1S4+P2S2+P3=0 (13)
利用t=S2对上式进行换元,得到有关t的一元二次方程:
P1t2+P2t+P3=0 (14)
其中,
P1=CDu+DFv-ABu-BEv (15)
根据式(14)求得t,进一步能够得到S;在求取t和S时,一元二次方程有两个解,舍去不合理的一个解;
(c)流体表面几何和速度的优化
按照以下步骤循环两次进行优化:步骤1,将步骤(b)得到的第k次非归一化几何看作是较准确的,利用以下公式(24)和(25)得到第k+1次更准确的速度(uk+1,vk+1);步骤2,将步骤1得到的速度看作是较准确的结果,再利用方程(13)计算得到非归一化因子S,从而得到更准确的表面几何;
每个粒子利用上述第L-1层的步骤(c)计算得到流体表面几何的非归一化因子,以下步骤是利用粒子的非归一化因子来计算整个场景的非归一化因子Sw;如果任意两个粒子i和j计算得到的非归一化因子分别为Si和Sj,聚类时,如果Si和Sj满足以下条件,那么Si和Sj就标记为同一类:
|Si-Sj|<ST i,j∈{0,1,2,3…N-1} (26)
其中N为粒子总数目,ST为非归一化因子差异的阈值,设其值为1.0;
聚类时类的数量是自动得到的,如果聚得最大类的样本数为n,最大类中所有粒子求得的非归一化因子的和为Ss,那么,
Sw=Ss/n (27)
(2)利用步骤(1)得到的整个场景的非归一化因子Sw,再计算视频样例帧非归一化几何
取视频中前100帧作为样例帧,按照步骤(1)得到的整个场景的非归一化因子Sw,再利用式(1)计算这些样例帧表面非归一化几何;
(3)采用LBM(Lattice Boltzmann Method),恢复流体运动的属性及驱动模型
对于所有样例帧,利用LBM进行仿真时,利用式(28)和(29)计算并调整粒子的粒子密度和分布函数;
其中,hL是粒子演化的高度,fi(t+1)和fi(t)分别是粒子在i方向的t+1时刻和t时刻的分布函数,i表示LBM的D2Q9模型中第i个方向,i=0,1,2…8,ρ(t+1)和ρ(t)分别是粒子在t+1时刻和t时刻的密度;
(4)对流体粒子仿真密度属性进行增强
将步骤(3)中得到的粒子密度我们称其为标准密度,记为ρs;LBM自动推演产生的密度,称之为自动演化密度,记为ρa;分析ρs的直方图的分布,并对ρa进行增强,增强后的密度称之为增强密度,记为ρe;
对于ρs,求其最大值和最小值将区间等分成256个子区间:
[iss,(i+1)ss] i={0,1…255} (30)
其中,ss为区间尺度:
标准密度ρs表示为ρs=Ss×ss+Os,将ρs记为[Ss,Os],其中Ss为区间号,Os为区间内密度位移量,其满足:
Os=ρs-Ss×ss (32)
对ρa,求其最大值和最小值将区间等分成256个子区间:
[isa,(i+1)sa] i={0,1…255} (33)
其中,sa为区间尺度:
对于任何一个自动演化密度ρa,都可以表示为ρa=Sa×sa+Oa,将ρa记为[Sa,Oa]其中Sa为区间号,Oa为区间内密度位移量,其满足:
Oa=ρa-Sa×sa (35)
根据ρs的直方图的统计分布,利用图像增强理论中直方图规定化的方法,实现Ss和Sa的对应,如果对于任意一个自动演化密度它的区间号为如果和它匹配的区间号为它被增强后的密度ρe可以计算为:
(5)重仿真过程
视频流体的重仿真过程包括以下步骤:
(a)利用步骤(1)初始化流体初始帧的运动速度及非归一化几何;
(b)利用步骤(2)计算样例帧表面非归一化几何;
(c)利用初始帧运动速度,推演LBM,同时利用步骤(3)和(4)增强流体粒子的密度;
(d)利用增强的密度演化仿真,得到具有细节特征的真实感流体重仿真结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410805222.0A CN104517299B (zh) | 2014-12-19 | 2014-12-19 | 视频流体物理驱动模型恢复及重新仿真的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410805222.0A CN104517299B (zh) | 2014-12-19 | 2014-12-19 | 视频流体物理驱动模型恢复及重新仿真的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104517299A CN104517299A (zh) | 2015-04-15 |
CN104517299B true CN104517299B (zh) | 2017-05-24 |
Family
ID=52792544
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410805222.0A Expired - Fee Related CN104517299B (zh) | 2014-12-19 | 2014-12-19 | 视频流体物理驱动模型恢复及重新仿真的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104517299B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204714B (zh) * | 2016-08-01 | 2019-02-01 | 华东师范大学 | 基于Phong模型的视频流体光照计算方法 |
CN106126865B (zh) * | 2016-08-01 | 2019-04-26 | 华东师范大学 | 真实感流体重仿真方法 |
CN106294955B (zh) * | 2016-08-01 | 2019-05-17 | 华东师范大学 | 一种物理数据驱动的真实感流体重仿真方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009155274A1 (en) * | 2008-06-16 | 2009-12-23 | Schlumberger Canada Limited | Streamline flow simulation of a model that provides a representation of fracture corridors |
CN101908086A (zh) * | 2010-07-09 | 2010-12-08 | 上海理工大学 | 一种对建筑室内风环境进行数字化动态模拟的分析方法 |
CN102819662A (zh) * | 2012-07-10 | 2012-12-12 | 华东师范大学 | 一种视频流体高度的计算方法 |
CN103294850A (zh) * | 2013-04-27 | 2013-09-11 | 苏州市数字城市工程研究中心有限公司 | 一种三维动态流体仿真算法智能匹配方法 |
CN103413346A (zh) * | 2013-04-09 | 2013-11-27 | 华东师范大学 | 一种真实感流体实时重建方法及其系统 |
-
2014
- 2014-12-19 CN CN201410805222.0A patent/CN104517299B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009155274A1 (en) * | 2008-06-16 | 2009-12-23 | Schlumberger Canada Limited | Streamline flow simulation of a model that provides a representation of fracture corridors |
CN101908086A (zh) * | 2010-07-09 | 2010-12-08 | 上海理工大学 | 一种对建筑室内风环境进行数字化动态模拟的分析方法 |
CN102819662A (zh) * | 2012-07-10 | 2012-12-12 | 华东师范大学 | 一种视频流体高度的计算方法 |
CN103413346A (zh) * | 2013-04-09 | 2013-11-27 | 华东师范大学 | 一种真实感流体实时重建方法及其系统 |
CN103294850A (zh) * | 2013-04-27 | 2013-09-11 | 苏州市数字城市工程研究中心有限公司 | 一种三维动态流体仿真算法智能匹配方法 |
Non-Patent Citations (3)
Title |
---|
computation of fluid and particle motion from a time-sequenced image pair: A global outlier identification approach;Nilanjan R 等;《IEEE Transations on Image Processing》;20111231;第20卷(第10期);2925-2936 * |
一种流体运动矢量计算的有效方法;全红艳 等;《计算机学报》;20130903;第36卷(第9期);1889-1897 * |
真实感流体实时重建;俞铭琪 等;《计算机辅助设计与图形学学报》;20130331;第25卷(第5期);622-630 * |
Also Published As
Publication number | Publication date |
---|---|
CN104517299A (zh) | 2015-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101404091B (zh) | 基于两步形状建模的三维人脸重建方法和系统 | |
Vanegas et al. | Modelling the appearance and behaviour of urban spaces | |
CN109147048A (zh) | 一种利用单张彩色图的三维网格重建方法 | |
CN100562895C (zh) | 一种基于区域分割和分段学习的三维人脸动画制作的方法 | |
Zhang et al. | A deformable surface model for real-time water drop animation | |
Halatsch et al. | Using shape grammars for master planning | |
Raveendran et al. | Blending liquids | |
Iizuka et al. | An interactive design system for pop-up cards with a physical simulation | |
Chao et al. | Realistic data-driven traffic flow animation using texture synthesis | |
CN102930583B (zh) | 互动式生成液滴效果的方法 | |
CN102708227A (zh) | 基于sph算法的洪水溃坝过程仿真方法与仿真系统 | |
Hu et al. | Sketch2VF: Sketch‐based flow design with conditional generative adversarial network | |
CN106027962A (zh) | 视频监控的覆盖率计算方法及装置、布点方法及系统 | |
CN104517299B (zh) | 视频流体物理驱动模型恢复及重新仿真的方法 | |
CN106127847B (zh) | 一种生成实时山水画特征线条以及适当渲染的方法 | |
CN107452058A (zh) | 基于特征保留的仿真模型分程度简化方法生成及系统 | |
CN106934192A (zh) | 一种参数优化的浅水方程模型水体建模方法 | |
CN107273617B (zh) | 一种利用浅水波方程获得表面流流体运动的实时模拟方法及系统 | |
Yan et al. | Interactive liquid splash modeling by user sketches | |
CN105006015A (zh) | 基于位置的流体模拟控制方法及系统 | |
CN105279781A (zh) | 基于多精度融合的流体动画生成方法 | |
CN104463934A (zh) | 一种“质点-弹簧”系统驱动的点集模型动画自动生成方法 | |
CN105321205B (zh) | 一种基于稀疏关键点的参数化人体模型重建方法 | |
CN110378023A (zh) | 屋面降雨径流模拟方法、装置、计算机设备和存储介质 | |
CN102819630B (zh) | 八叉树粗糙网格投影的快速的烟模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170524 Termination date: 20191219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |