发明名称 真实感流体重仿真方法
摘要 本发明公开了一种真实感流体的重仿真方法,该方法旨在利用视频流体的样例帧,得到真实感的重仿真结果。该方法首先采用计算样例帧流体的物理驱动数据,通过小波变换,在频域空间中获取驱动数据的高频成分,进而在流体重仿真的过程中,利用提取的高频细节成分对LBM仿真过程中的驱动数据进行增强,从而得到了具有真实感细节的重仿真结果。本发明可有效地应用于流体的反向工程的研究及应用中,可以实现流体场景真实感的重仿真,获得与样例流体外观接近的真实感效果,克服了现有的流体真实感中缺少细节的问题。
申请公布号 CN106126865A 申请公布日期 2016.11.16
申请号 CN201610620770.5 申请日期 2016.08.01
申请人 华东师范大学 发明人 全红艳
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 上海蓝迪专利商标事务所(普通合伙) 31215 代理人 徐筱梅;张翔
主权项 一种真实感流体重仿真方法,其特征在于,利用视频的前两个样例帧,记为F<sub>1</sub>和F<sub>2</sub>,利用密度的频域中高频成分,实现流体真实感重仿真,具体步骤如下:步骤1:预处理计算F<sub>1</sub>中任意一个粒子的高度H<sub>i</sub>,利用F<sub>1</sub>和F<sub>2</sub>计算F<sub>1</sub>中粒子的速度v<sub>i</sub>及密度P<sub>i</sub>,其中i是非负整数,且0≤i&lt;N,N是流体场景的粒子个数;(a)建立标准密度的频谱块集合S<sub>F</sub>,具体地:将F<sub>1</sub>中的所有粒子分块处理,并按照每块中粒子的运动方向,将所有的块确定为不同类别;先进行粒子块的划分与标记,将F<sub>1</sub>中的所有粒子按照指定大小b×b分成正方形的块,如果图像的分辨率为X×Y,那么,水平分割块数<img file="FDA0001066814170000011.GIF" wi="242" he="143" />垂直分割块数<img file="FDA0001066814170000012.GIF" wi="219" he="143" /><img file="FDA0001066814170000013.GIF" wi="76" he="75" />表示向下取整;利用这WA个频谱块构建频谱块集合S<sub>F</sub>;其中<img file="FDA0001066814170000014.GIF" wi="67" he="70" />为F<sub>1</sub>中任意一个块,它位于频谱空间中水平的第w块和垂直的第h块,w和h为非负整数,0≤w&lt;W,0≤h&lt;A;假设F<sub>1</sub>中任意一个块<img file="FDA0001066814170000015.GIF" wi="66" he="71" />中所有粒子2D速度的均值为<img file="FDA0001066814170000016.GIF" wi="64" he="53" />并且u是<img file="FDA0001066814170000017.GIF" wi="36" he="47" />的水平方向分量、v是<img file="FDA0001066814170000018.GIF" wi="37" he="46" />的垂直方向的分量,那么<img file="FDA0001066814170000019.GIF" wi="67" he="71" />的类型T利用(1)式确定:<maths num="0001"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>u</mi><mo>&gt;</mo><mn>0</mn><mo>,</mo><mi>v</mi><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>2</mn></mtd><mtd><mrow><mi>u</mi><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>v</mi><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>3</mn></mtd><mtd><mrow><mi>u</mi><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>v</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>4</mn></mtd><mtd><mrow><mi>u</mi><mo>&gt;</mo><mn>0</mn><mo>,</mo><mi>v</mi><mo>&lt;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010668141700000110.GIF" wi="1199" he="318" /></maths>计算F<sub>1</sub>中所有粒子密度,称其为标准密度,将标准密度按照指定层数进行离散小波变换,假设层数为L,那么变换后的高频区一共有3L个,对于这3L个高频区中任意一个区A<sub>k</sub>(0≤k&lt;3L),按照上述方法分块:水平分割为W块,垂直分割为A块;在每个高频区中,将WA个频谱块按如下方法构建四个集合<img file="FDA00010668141700000111.GIF" wi="322" he="71" />具体地:对于A<sub>k</sub>(0≤k&lt;3L)中任意一个块b<sub>j</sub>(0≤j&lt;W*H),如果A<sub>k</sub>在S<sub>F</sub>中对应块的类型假设为T,那么,将b<sub>j</sub>置入集合<img file="FDA00010668141700000112.GIF" wi="301" he="71" />中,按照这样方法,构建每个区域四个频谱的集合<img file="FDA00010668141700000113.GIF" wi="561" he="70" />称<img file="FDA00010668141700000114.GIF" wi="60" he="55" />为标准密度频谱块集;(b)对每个集合内部的频谱块排序,具体地:统计每个集合<img file="FDA0001066814170000021.GIF" wi="542" he="71" />内部所有粒子的频谱和,并将集合<img file="FDA0001066814170000022.GIF" wi="65" he="70" />的频谱块,按照统计结果从大到小排序,排序后的集合记为<img file="FDA0001066814170000023.GIF" wi="557" he="70" />称<img file="FDA0001066814170000024.GIF" wi="53" he="63" />为标准密度频谱块有序集;步骤2:按照以下步骤,共进行N<sub>R</sub>次仿真推演:(a)当前次是第I次仿真,设I=1;如果满足I≤N<sub>R</sub>,N<sub>R</sub>仿真一共进行的推演次数,继续进行以下步骤,否则,跳转到步骤3;(b)利用速度v<sub>i</sub>(0≤i&lt;Z),初始化LBM模型,使仿真开始推演,然后按照以下步骤对推演后的密度进行增强,Z是流体场景的粒子个数,如果流体视频F<sub>1</sub>的分辨率为g×m,那么Z可以计算为Z=gm;在推演过程中,获取第I次推演中流体场景粒子的密度p<sub>l</sub>(0≤l&lt;Z),p<sub>l</sub>为推演密度,Z是流体场景的粒子个数,其中I是自然数(1≤I≤N<sub>R</sub>);然后将p<sub>l</sub>按照步骤1的(a)和(b)建立p<sub>l</sub>的频谱块有序集<img file="FDA0001066814170000025.GIF" wi="579" he="71" />(c)推演密度的增强假设标准密度频谱块有序集<img file="FDA0001066814170000026.GIF" wi="530" he="71" />中有c个频谱块、推演密度的频谱块有序集<img file="FDA0001066814170000027.GIF" wi="554" he="70" />中有e个频谱块;对于指定的阈值θ(0&lt;0&lt;1),对于第m次推演中,<img file="FDA0001066814170000028.GIF" wi="557" he="71" />中前θc个频谱块进行增强处理;具体地:推演密度的频谱块b<sub>i</sub>块号i,计算标准密度频谱块B<sub>s</sub>块号s:s=ci/e     (2)选取标准密度频谱块B<sub>s</sub>,然后对于b<sub>i</sub>中粒子的密度频谱p<sub>t</sub>增强为:p<sub>t</sub>=K<sub>t</sub>Q<sub>b</sub>/X<sub>B</sub>       (3)其中K<sub>t</sub>为B<sub>s</sub>中序号为t的粒子的密度频谱;Q<sub>b</sub>为<img file="FDA0001066814170000029.GIF" wi="78" he="63" />频谱块中所有频谱的和,X<sub>B</sub>为<img file="FDA00010668141700000210.GIF" wi="61" he="70" />频谱块中所有频谱的和;(d)转步骤2,进行下一次推演。步骤3:重仿真结束。
地址 200241 上海市闵行区东川路500号