Python OpenGL绘制一场烟花盛会

忙碌了一年，今天终于放假了。原本打算好好休息一下，没成想只过了半天就觉得有点无聊。看家人和朋友们都在忙年，那我就用opengl导演一场烟花盛会，献给即将到来的新年吧。

一说到opengl，很多人都会觉得复杂，其实不然。只要掌握了几个基本的概念，借助于工具软件，任何人都可以很轻松地上手。在制作烟花之前，我先介绍一下wxgl这个三维数据快速可视化工具。

1. 安装wxgl

wxgl是一个基于pyopengl的三维数据可视化库，以wx为显示后端，提供matplotlib风格的交互式应用方式。wxgl也可以和wxpython无缝结合，在wx的窗体上绘制三维模型。使用pip命令即可快速安装wxgl及其所依赖的其他模块。

pip install wxgl

2. 快速体验

下面这几行代码，绘制了一个中心在坐标原点半径为1的纯色圆球。忽略模块名的话，这些代码和matplotlib的风格是完全一致的。

>>> import wxgl.wxplot as plt
>>> plt.uvsphere((0,0,0), 1, color='cyan')
>>> plt.title('快速体验：$x^2+y^2=1$')
>>> plt.show()

生成一个地球模型是如此简单。

>>> plt.uvsphere((0,0,0), 1, texture='res/earth.jpg', xflip=true, yflip=false)
>>> plt.show()

让地球自转，更是易如反掌。

>>> plt.uvsphere((0,0,0), 1, 
    texture='res/earth.jpg', 
    xflip=true, 
    yflip=false,
    transform = lambda tn,gms,tms : ((0, 1, 0, (0.01*tms)%360),)
)
>>> plt.show()

勾选“屏幕录制”，点击“播放”按钮，即可保存为gif文件或mp4/avi格式的视频文件。

这是代码中用的的地球纹理图片，可以直接下载使用。

3. 编写自己的着色器

wxgl不仅提供了线段、散点、曲面、三维等值面等一系列绘图函数，还支持用户定制着色器程序，以实现更复杂的功能。下面这个例子，用粒子技术模拟了烟花升空的过程。

# -*- coding: utf-8 -*-

import numpy as np
import wxgl
import wxgl.wxplot as plt

def rise(n, pos, h, v, a, cycle):
    """烟花升空模型
    
    n       - 粒子数量
    pos     - 初始位置
    h       - 上升行程
    v       - 初始速度
    a       - 上升加速度
    cycle   - 循环周期
    """
    
    vshader_src = """
        #version 330 core
        in vec4 a_position;
        in vec4 a_color;
        in float a_delay; // 粒子发射延迟时间（s）
        uniform float u_ts; // 持续时间（s）
        uniform float u_v; // 初始速度
        uniform float u_a; // 上升加速度
        uniform mat4 u_mvpmatrix;
        out vec4 v_color;
        out float v_ts;
        void main() { 
            float t = u_ts - a_delay;
            if (t < 0) t = 0;
            
            float s = u_v * t + 0.5 * u_a * t * t;
            gl_position = u_mvpmatrix * vec4(a_position.x, a_position.y+s, a_position.z, a_position.w); 
            gl_pointsize = 1;
            v_color = a_color;
            v_ts = u_ts;
        }
    """
    
    fshader_src = """
        #version 330 core
        in vec4 v_color;
        uniform float u_tmax;
        in float v_ts;
        void main() { 
            if(v_ts > u_tmax) discard;
            
            vec2 temp = gl_pointcoord - vec2(0.5);
            float f = dot(temp, temp);
            if(f > 0.25) discard;
            
            gl_fragcolor = vec4(v_color.rgb, 1); 
        } 
    """
    
    vs = np.array(pos) + (np.random.random((n,3)) - 0.5) * h/100
    color = np.tile(np.array((1.0,1.0,0.8)), (n,1))
    delay = np.float32(np.absolute(np.random.randn(n))) / 10
    tmax = (pow(v*v+2*a*h, 0.5)-v)/a + delay.max()
    
    m = wxgl.model(wxgl.points, vshader_src, fshader_src, sprite=true)
    m.set_vertex('a_position', vs)
    m.set_color('a_color', color)
    m.set_argument('a_delay', delay)
    m.set_argument('u_ts', lambda tn,gms,tms:(tms/1000)%cycle)
    m.set_argument('u_v', v)
    m.set_argument('u_a', a)
    m.set_argument('u_tmax', tmax)
    m.set_mvp_matrix('u_mvpmatrix') # 设置模型矩阵、视点矩阵和投影矩阵
    
    return m

vs = np.array([
    [-1.5,2,1], [-1.5,0,1], [1.5,2,1], [1.5,0,1], 
    [-1.5,2,-1], [-1.5,0,-1], [1.5,2,-1], [1.5,0,-1]])
vs = vs[[0,1,2,3,0,2,1,3,4,5,6,7,4,6,5,7,0,4,1,5,2,6,3,7]]
m = rise(n=500, pos=(0,0,0), h=1.5, v=2, a=-1.2, cycle=5)

plt.figure(zoom=0.7, elev=10)
plt.line(vs, color=(0,1,1), method='isolate') # 六面体线框，表示烟花燃放的空间
plt.model(m)
plt.show()

4. 绽放的烟花

只要理解了烟花升空的代码，很容易写出烟花在空中爆炸的着色器程序。下面的代码除了烟花升空的着色器，还提供了两种烟花爆炸的着色器，其中用到了一个纹理图片，可直接下载下面这张图使用。

如果将上面的纹理图片替换成文字，就可以在烟花爆炸的瞬间显示出文字了。wxgl提供了一个文本转pil图形对象的函数，可以直接作为纹理使用。

# -*- coding: utf-8 -*-

import numpy as np
import wxgl
import wxgl.wxplot as plt

def rise(n, pos, h, v, a, cycle):
    """烟花升空模型
    
    n       - 粒子数量
    pos     - 初始位置
    h       - 上升行程
    v       - 初始速度
    a       - 上升加速度
    cycle   - 循环周期
    """
    
    vshader_src = """
        #version 330 core
        in vec4 a_position;
        in vec4 a_color;
        in float a_delay; // 粒子发射延迟时间（s）
        uniform float u_ts; // 持续时间（s）
        uniform float u_v; // 初始速度
        uniform float u_a; // 上升加速度
        uniform mat4 u_mvpmatrix;
        out vec4 v_color;
        out float v_ts;
        void main() { 
            float t = u_ts - a_delay;
            if (t < 0) t = 0;
            
            float s = u_v * t + 0.5 * u_a * t * t;
            gl_position = u_mvpmatrix * vec4(a_position.x, a_position.y+s, a_position.z, a_position.w); 
            gl_pointsize = 1;
            v_color = a_color;
            v_ts = u_ts;
        }
    """
    
    fshader_src = """
        #version 330 core
        in vec4 v_color;
        uniform float u_tmax;
        in float v_ts;
        void main() { 
            if(v_ts > u_tmax) discard;
            
            vec2 temp = gl_pointcoord - vec2(0.5);
            float f = dot(temp, temp);
            if(f > 0.25) discard;
            
            gl_fragcolor = vec4(v_color.rgb, 1); 
        } 
    """
    
    vs = np.array(pos) + (np.random.random((n,3)) - 0.5) * h/100
    color = np.tile(np.array((1.0,1.0,0.8)), (n,1))
    delay = np.float32(np.absolute(np.random.randn(n))) / 10
    tmax = (pow(v*v+2*a*h, 0.5)-v)/a + delay.max()
    
    m = wxgl.model(wxgl.points, vshader_src, fshader_src, sprite=true)
    m.set_vertex('a_position', vs)
    m.set_color('a_color', color)
    m.set_argument('a_delay', delay)
    m.set_argument('u_ts', lambda tn,gms,tms:(tms/1000)%cycle)
    m.set_argument('u_v', v)
    m.set_argument('u_a', a)
    m.set_argument('u_tmax', tmax)
    m.set_mvp_matrix('u_mvpmatrix') # 设置模型矩阵、视点矩阵和投影矩阵
    
    return m, tmax

def bomb_1(n, pos, start, a, cycle):
    """烟花爆炸模型
    
    n       - 粒子数量
    pos     - 位置
    start   - 时间
    a       - 下降加速度
    cycle   - 循环周期
    """
    
    vshader_src = """
        #version 330 core
        in vec4 a_position;
        in vec3 a_data;
        uniform float u_ts;
        uniform float u_start;
        uniform float u_a;
        uniform mat4 u_mvpmatrix;
        out vec4 v_color;
        out float v_ts;
        void main() { 
            float t = u_ts - u_start;
            if (t < 0) t = 0;
            
            float lat = radians((a_data.x - 0.5) * 90);
            float lon = radians(a_data.y * 360);
            float r = (a_data.z * 0.3 + 0.7) * 0.3 * t * (1 + 0.3 * a_position.z);
            float y = r * sin(lat) + a_position.y - 0.5*u_a*t*t;
            float xz = r * cos(lat);
            float x = xz * cos(lon) + a_position.x;
            float z = xz * sin(lon) + a_position.z;
            
            gl_position = u_mvpmatrix * vec4(x,y,z,a_position.w); 
            gl_pointsize = 3 * t;
            v_ts = t;
            
            int i = gl_vertexid % 6;
            if (i == 0) v_color = vec4(1,0,0,1);
            else if (i == 1) v_color = vec4(0,1,0,1);
            else if (i == 2) v_color = vec4(0,0,1,1);
            else if (i == 3) v_color = vec4(1,1,0,1);
            else if (i == 4) v_color = vec4(0,1,1,1);
            else v_color = vec4(1,0,1,1);
        }
    """
    
    fshader_src = """
        #version 330 core
        in vec4 v_color;
        in float v_ts;
        void main() { 
            if(v_ts <= 0 || v_ts > 2) discard;
            
            vec2 temp = gl_pointcoord - vec2(0.5);
            float f = dot(temp, temp);
            if(f > 0.25)  discard;
            
            //float alpha = v_color.a * exp(1-30*f) * (4-v_ts*v_ts)/2;
            float alpha = v_color.a * (1-4*f) * (4-v_ts*v_ts)/2;
            gl_fragcolor = vec4(v_color.rgb, alpha); 
        } 
    """
    
    vs = np.tile(np.array(pos), (n,1))
    data = np.float32(np.random.random((n,3)))
    
    m = wxgl.model(wxgl.points, vshader_src, fshader_src, sprite=true, opacity=false)
    m.set_vertex('a_position', vs)
    m.set_argument('a_data', data)
    m.set_argument('u_start', start)
    m.set_argument('u_a', a)
    m.set_argument('u_ts', lambda tn,gms,tms:(tms/1000)%cycle)
    m.set_mvp_matrix('u_mvpmatrix') # 设置模型矩阵、视点矩阵和投影矩阵
    
    return m

def bomb_2(pos, start, texture, a, size, cycle):
    """烟花爆炸模型
    
    pos     - 位置
    start   - 时间
    texture - 纹理
    a       - 下降加速度
    cycle   - 循环周期
    """
    
    vshader_src = """
        #version 330 core
        in vec4 a_position;
        uniform float u_ts;
        uniform float u_start;
        uniform float u_a;
        uniform float u_size;
        uniform mat4 u_mvpmatrix;
        out float v_ts;
        void main() { 
            float t = u_ts - u_start;
            if (t < 0) t = 0;
            
            if (t < 2) gl_pointsize = t * u_size/2 * (1 + 0.3 * a_position.z);
            else gl_pointsize = u_size * (1 + 0.3 * a_position.z);
            
            gl_position = u_mvpmatrix * vec4(a_position.x, a_position.y-0.5*u_a*t*t, a_position.z, a_position.w); 
            v_ts = t;
        }
    """
    
    fshader_src = """
        #version 330 core
        uniform sampler2d u_fireworks;
        in float v_ts;
        void main() { 
            if(v_ts <= 0 || v_ts > 2) discard;
            
            vec4 color = texture2d(u_fireworks, gl_pointcoord);
            gl_fragcolor = vec4(color.rgb, color.a*(4-v_ts*v_ts)/2);
        } 
    """
    
    vs = np.array(pos).reshape(-1,3)
    
    m = wxgl.model(wxgl.points, vshader_src, fshader_src, sprite=true)
    m.set_vertex('a_position', vs)
    m.set_argument('u_a', a)
    m.set_argument('u_size',size)
    m.set_argument('u_start', start)
    m.set_argument('u_ts', lambda tn,gms,tms:(tms/1000)%cycle)
    m.add_texture('u_fireworks', texture, wxgl.texture_2d, yflip=false)
    m.set_mvp_matrix('u_mvpmatrix') # 设置模型矩阵、视点矩阵和投影矩阵
    
    return m

if __name__ == '__main__':
    vs = np.array([[-1.5,2,1], [-1.5,0,1], [1.5,2,1], [1.5,0,1], [-1.5,2,-1], [-1.5,0,-1], [1.5,2,-1], [1.5,0,-1]])
    vs = vs[[0,1,2,3,0,2,1,3,4,5,6,7,4,6,5,7,0,4,1,5,2,6,3,7]]
    
    plt.figure(zoom=0.5, elev=10)
    plt.line(vs, color=(0,1,1,0), method='isolate') # 六面体线框，表示烟花燃放的空间
    
    # ------------------------------
    h, v, a, cycle = 1.7, 2.2, -1.2, 4
    for i, ch in enumerate('新春快乐'):
        x = -1.5 + i
        m1, start = rise(n=300, pos=(x,0,1), h=h, v=v, a=a, cycle=cycle)
        m2 = bomb_1(200, (x,h,1), start, a=0.1, cycle=cycle)
        m3 = bomb_2((x,h,1), start, wxgl.text2image(ch, 96, (1,0,0)), a=0.1, size=100, cycle=cycle)
        plt.model(m1)
        plt.model(m2)
        plt.model(m3)
    
    # -------------------------------
    for i in range(20):
        x, z = (np.random.random()-0.5)*4, (np.random.random()-0.5)*2
        h, v, a = 1.5+(np.random.random()-0.5)*0.4, 2.2, -1.2
        cycle = np.random.randint(4, 7)
        m1, start = rise(n=300, pos=(x,0,z), h=h, v=v, a=a, cycle=cycle)
        m2 = bomb_1(200, (x,h,z), start, a=0.1, cycle=cycle)
        plt.model(m1)
        plt.model(m2)
    
    # -------------------------------
    for i in range(20):
        x, z = (np.random.random()-0.5)*4, (np.random.random()-0.5)*2
        h, v, a = 1.5+(np.random.random()-0.5)*0.4, 2.3, -1.2
        cycle = np.random.randint(4, 7)
        m1, start = rise(n=300, pos=(x,0,z), h=h, v=v, a=a, cycle=cycle)
        m2 = bomb_2((x,h,z), start, 'res/fw.png', a=0.1, size=300, cycle=cycle)
        plt.model(m1)
        plt.model(m2)
    
    plt.show()    

最终的效果如下面的gif所示。

到此这篇关于python opengl绘制一场烟花盛会的文章就介绍到这了,更多相关python opengl烟花内容请搜索以前的文章或继续浏览下面的相关文章希望大家以后多多支持！