水面的简单渲染 – Gerstner波

发布于 2014年2月7日 作者： John Hany2,986次阅读

        渲染三维场景时经常会遇到需要渲染各种水体的情况，比如湖泊、河流、海洋等，不仅需要水体表面要有接近真实的随时间而变化的波动，还要有令人信服的颜色、反射光、透明度等细节。实时渲染水面的方法有很多，从简单的若干正弦波叠加，到《GPU Gems》中介绍的叠加Gerstner波的方法，再到如今GPU在线计算FFT得到表面高度，都是以追求效果更加逼真的同时保证计算的高效实时。我用OpenGL实现了通过合成Gestner波产生水波的方法，具体过程如下。

开发环境

Windows 7 x64，Visual Studio 2010，OpenGL版本3.0，GLSL版本1.3。

freeglut 2.8.1，GLM 0.9.5.1。GLM用于产生模型视图矩阵、透视投影矩阵和法线变换矩阵。

正弦函数波

        在一些数学书中介绍正弦函数时会提到“理想情况下的水波是正弦形状的”，但实际上，单独的水波应该是波峰尖、波谷宽的。如果用正弦波来表现这样的效果，可以选择如下变换：

equation-1
image-469
        由于正弦函数的值域是[-1,1]，缩放到[0,1]区间，再做幂运算，会使函数值减小，而且距离0越小的值减小得越多。这样就能产生波峰尖、波谷宽的形状。

 

        下面是一组k分别等于1.0，1.5和2.0时的情况，可见k越大（k>=1），形状就越明显。

sin-multi
image-470

 

        但是只有一个参数决定这种形状过于简单，而且在CG中希望在细节多的地方（波峰）网格点较密集，在细节少的地方（波谷）网格点较稀疏。用正弦函数绘制时，如果想提高细节，只能整体提高x的细分程度，也会在波谷处增加大量的多余计算。

Gerstner波

        Gerstner波的诞生早于计算机图形学（CG），它最初在物理中用于水波的模拟。由于它的形状比较真实，而且计算量不大，所以被广泛用于CG中水波的模拟。

        Gerstner波以参数方程的形式给出：

equation-2
image-471
        自变量为p，参数Q、D、A用来控制形状。Q控制波峰的尖锐度，D控制波长，A为振幅。

 

gerstner-single
image-472

        Q应为较小的值，若Q很大，会在波峰处产生环，破坏波的形状。比如：

 

gerstner-circle
image-473
        观察x(p)的表达式可以看出，与正弦波相比，Gerstner波在波峰处的点更紧凑，在波谷处更稀疏：

 

gerstner-and-sin
image-474

 

波的合成

        为了产生真实的水面，需要把若干不同方向、不同参数的Gerstner波合成为一个波。即：

equation-3
image-475
        在三维空间中绘制水波这样高度值频繁变化的面时，一般采用规则网格来绘制，即在x-y平面上画一张均匀的网格，对网格上的每一个点计算它的高度值（z值），这样就产生了一张高低起伏的面。随着时间的变化，每个点的高度也随之变化，就产生了动态的面。

 

        为了把这张网格与二维的Gerstner波结合起来，需要进行如下转换：

        假设二维Gerstner波表示为y=f(x)，三维网格表示为z=g(x,y)。则：

1564654
image-476

 

        (x0,y0)表示波的起点，theta角表示波传播的方向。

grid
image-477
        初始时，网格上每点的高度设为0，每叠加一个波，就根据上面的式子计算出一个高度，加在z上。计算完所有的波后，就实现了多个波的叠加。

 

        由于Gerstner参数方程也在改变x（即上图的d），直接应用原式计算会增加复杂度。同时，为了尽可能地减小计算量，我采用两种固定形状的Gerstner波，每种波用11对坐标表示，计算f(x)时只需要在这11个点中计算线性内插即可。

两种波形。第一个波峰较尖，用来绘制细小的水波，第二个波峰较宽，用来绘制波长较长的水波。

1
2
3
4
5
6
7
8
9
10
static const GLfloat gerstner_pt_a[22] = {
    0.0,0.0, 41.8,1.4, 77.5,5.2, 107.6,10.9,
    132.4,17.7, 152.3,25.0, 167.9,32.4, 179.8,39.2,
    188.6,44.8, 195.0,48.5, 200.0,50.0
};
static const GLfloat gerstner_pt_b[22] = {
    0.0,0.0, 27.7,1.4, 52.9,5.2, 75.9,10.8,
    97.2,17.6, 116.8,25.0, 135.1,32.4, 152.4,39.2,
    168.8,44.8, 184.6,48.5, 200.0,50.0
};

线性内插函数：

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
static float gerstnerZ(float w_length, float w_height, float x_in, const GLfloat gerstner[22])
{
    x_in = x_in * 400.0 / w_length;
 
    while(x_in < 0.0)
        x_in += 400.0;
    while(x_in > 400.0)
        x_in -= 400.0;
    if(x_in > 200.0)
        x_in = 400.0 - x_in;
 
    int i = 0;
    float yScale = w_height/50.0;
    while(i<18 && (x_in<gerstner[i] || x_in>=gerstner[i+2]))
        再分享一下我老师大神的人工智能教程吧。零基础！通俗易懂！风趣幽默！还带黄段子！希望你也加入到我们人工智能的队伍中来！https://blog.csdn.net/jiangjunshow

posted on 2019-01-31 21:53  jack船长大哥  阅读(1558)  评论(0)  编辑  收藏  举报

导航