简介：面向计算机图形学初学者与OpenGL开发者的Bezier曲面演示程序，基于Visual C++工程完整实现了曲面建模、纹理映射与光照渲染的整合。源代码通过GLU函数库构建二维Bezier曲面，涵盖gluBeginSurface、gluEndSurface与gluBuild2DMultitexture等接口；纹理部分演示glBindTexture、glTexImage2D及glTexCoord2f的用法，并开启GL_TEXTURE_2D实现贴图；光照部分则设置环境光、漫反射与镜面高光，配合glLight、glMaterial以及glNormal3f定义表面法线，清晰展示各要素如何影响最终渲染效果。压缩包共13个文件，以.cpp源程序、.dsw/.dsp工程文件、可运行的exe可执行文件为主，同时包含pdb、idb等调试符号文件，便于单步跟踪和二次开发。资源整体仅1.47MB，结构紧凑。目前已有454人学习下载，适合希望快速理解OpenGL中Bezier曲面、纹理与光照结合原理的读者直接参考。
1. Bezier曲面不是“画出来”的：先有顶点、法线和UV，才有光照
很多第一次接触Bezier曲面的读者，照着教程启用了glMap2、开了GL_LIGHTING，看到屏幕上有彩色曲面出现便觉得已经完成。可是旋转视角时高光像贴纸一样固定在屏幕上，或者纹理在曲面接缝处错开半格，问题多半出在同一个根源上：Bezier曲面本身不携带法线和纹理坐标，网格生成阶段必须从控制点的偏导数单独算出法线，再从参数域映射出UV。以下按构造曲面、求法线、选光照模型、绑定纹理这条流水线逐层展开，每一步都和OpenGL状态机对应，手边有OpenGL环境配置基础的读者可以直接对照实验。
2. OpenGL中生成Bezier曲面网格：从数学控制点到细分三角片
很多资料会把Bezier曲面的公式摆出来，然后直接跳到glEvalMesh2。但对于现代OpenGL核心profile，求值器已不再可用，因此我更习惯先讲清控制点到网格的数学，再给出手动细分和求值器两条路线。
2.1 Bezier曲面的两个关键偏导数：如何用控制点直接算切线
Bezier曲面是控制点矩阵和两条Bernstein基函数基的张量积，写成：
P(u,v) = Σ_i Σ_j B_i^n(u) B_j^m(v) P_ij
在实际编码时，我不会直接展开这层求和，而是预先把控制点按行存成矩阵。3次Bezier曲面有4x4个控制点，对u求偏导，等价于把n次基函数降阶为n-1次，同时把控制点序列变成相邻控制点之差再乘以n，对v方向的偏导也同理。要得到曲面法线，先对u求导得到切向量du，对v求导得到切向量dv，再取叉积cross(du,dv)。
这个叉积顺序一旦写反，法线就会朝里，光照会呈现完全错误的暗面。所以我在细分循环里会先画一段临时法线验证朝向，再继续后续工作。
2.2 使用OpenGL求值器快速验证：glMap2/glEvalMesh2的参数含义
OpenGL 2.1固定管线提供glMap2和glEvalMesh2，适合快速验证控制点是否合理。代码：
GLfloat ctrlpoints[4][4][3] = { {{-1.5f, -1.5f, 4.0f}, {-0.5f, -1.5f, 2.0f}, {0.5f, -1.5f, -1.0f}, {1.5f, -1.5f, 2.0f}}, {{-1.5f, -0.5f, 1.0f}, {-0.5f, -0.5f, 3.0f}, {0.5f, -0.5f, 0.0f}, {1.5f, -0.5f, -1.0f}}, {{-1.5f, 0.5f, 4.0f}, {-0.5f, 0.5f, 0.0f}, {0.5f, 0.5f, 3.0f}, {1.5f, 0.5f, 2.0f}}, {{-1.5f, 1.5f, -2.0f},{-0.5f, 1.5f, 2.0f}, {0.5f, 1.5f, 0.0f}, {1.5f, 1.5f, -1.0f}} }; glEnable(GL_MAP2_VERTEX_3); glMap2f(GL_MAP2_VERTEX_3, 0.0f, 1.0f, 3, 4, 0.0f, 1.0f, 12, 4, &ctrlpoints[0][0][0]); glEnable(GL_AUTO_NORMAL); glEnable(GL_NORMALIZE); glMapGrid2f(20, 0.0f, 1.0f, 20, 0.0f, 1.0f); glEvalMesh2(GL_FILL, 0, 20, 0, 20);这里的glMap2f参数序列是u范围、v范围、u方向跨度、u方向阶数、v方向跨度、v方向阶数和数据指针。具体来说，3代表每个顶点占用3个GLfloat，4代表u和v方向各有4个控制点，12代表v方向每跨一行需要跳过12个浮点数。glEnable(GL_AUTO_NORMAL)让求值器自动把偏导数叉积作为法线，GL_NORMALIZE防止非均匀缩放导致法线长度偏离1。
这个方案只在固定管线里有效，在核心profile中会退化。我通常只用它做控制点飞线的快速取证，真正出图还是走手动细分。
固定管线求值器与现代VBO细分的差异如下表：
| 对比项 | 固定管线求值器 | 手动细分VBO |
|---|---|---|
| GPU控制权 | 无法在顶点里附加自定义属性 | 可以同时写入位置、法线、UV、切线 |
| 法线生成 | 依赖GL_AUTO_NORMAL自动计算 | 自己选择解析法线或数值法线 |
| 纹理坐标 | 需要单独用glEvalCoord2控制 | 参数域(u,v)直接映射为UV |
| 核心兼容性 | OpenGL 3.0后废弃 | 所有OpenGL版本都稳定 |
2.3 手动细分路线：固定步长生成顶点、法线与UV的代码模板
手动细分的核心是把(u,v)参数域均匀切分成step份，然后生成(step+1)*(step+1)个顶点。下面是一个可复用的C++函数模板：
#include <glm/glm.hpp> #include <vector> using glm::vec2; using glm::vec3; struct BezierVertex { vec3 position; vec3 normal; vec2 uv; }; using ControlGrid = vec3[4][4]; vec3 evalBezier(const ControlGrid& c, float u, float v) { // 用 Bernstein 基函数或 de Casteljau 实现 return c[0][0]; // 占位 } vec3 evalBezierDx(const ControlGrid& c, float u, float v) { // 返回 dP/du return vec3(0.0f); // 占位 } vec3 evalBezierDy(const ControlGrid& c, float u, float v) { // 返回 dP/dv return vec3(0.0f); // 占位 } void generateBezierMesh(const ControlGrid& ctrl, int step, std::vector<BezierVertex>& out) { out.clear(); for (int j = 0; j <= step; ++j) { float v = static_cast<float>(j) / static_cast<float>(step); for (int i = 0; i <= step; ++i) { float u = static_cast<float>(i) / static_cast<float>(step); BezierVertex vert; vert.position = evalBezier(ctrl, u, v); vec3 du = evalBezierDx(ctrl, u, v); vec3 dv = evalBezierDy(ctrl, u, v); vert.normal = glm::normalize(glm::cross(du, dv)); vert.uv = vec2(u, v); out.push_back(vert); } } }2.3.1 细分步长STEP与顶点数量控制
step参数是曲面质量和渲染性能之间的杠杆。step=16时每个面片有289个顶点和512个三角形；step=32时顶点变成1089个，三角形数量翻了一倍多。我一般先用step=8快速看形状，确认控制点没有飞点后再提到16。如果曲面在屏幕上的投影小于200像素，step=8和16的差异很难察觉。
2.3.2 为什么法线在GenerateBezier中就要算出来
法线如果不在这里提前算好，即使OpenGL在光栅化时对顶点法线插值，得到的也只是顶点位置附近法线的近似，在光滑的Bezier曲面上会形成明显马赫带。所以在网格生成阶段直接调用导数函数求cross，比在着色器里用ddx/ddy去估稳定得多。
接下来把顶点和索引上传到VBO，用glDrawElements绘制，就能看到纯几何线框。这一步不需要光照，只需一个简单的把顶点位置映射到颜色的shader即可。
3. Bezier曲面光照模型：逐像素Blinn-Phong的GLSL实现与参数调整
3.1 固定管线光照模型的局限：为什么Bezier曲面需要可编程管线
固定管线的GL_LIGHTING提供的是逐顶点Phong光照，也就是Gouraud着色。Bezier曲面本身很光滑，但顶点数量受限于step，顶点之间由三角形线性插值，高光很容易出现多边形棱边。我在step=16的曲面上试过固定管线，镜面高光看起来像被折叠过。所以后续一定要使用可编程管线，把光照计算放到片元着色器，让每个像素独立计算明暗。
3.2 解析法线 vs 数值法线：两种求法及性能边界
在2.3节中，我使用解析法线。解析法线直接来自Bernstein基函数的导数，只与当前(u,v)坐标有关，不依赖网格密度。数值法线则通过对相邻顶点点做差再叉积得到，在控制点间距不均匀时会引入额外误差。比如step较小、相邻顶点跨越较长曲面弧长时，数值法线会把弧线的“割线”方向当成切向，导致法线偏移。
我通常只在无法使用解析导数的情况下才用数值法线，例如从外部导入的无参数网格。对Bezier曲面这条链路，解析法线几乎是零成本。
3.3 GLSL渲染Bezier曲面的光照模型着色器：漫反射、镜面反射与指数设置
3.3.1 顶点着色器：传递法线、灯向量与半程向量
以下是一个可直接嵌入现代OpenGL管线的顶点着色器。它把位置、法线和UV直接传递给片元着色器，光照核心放在片元阶段：
#version 330 core layout(location = 0) in vec3 aPos; layout(location = 1) in vec3 aNormal; layout(location = 2) in vec2 aUV; uniform mat4 modelViewMatrix; uniform mat4 projectionMatrix; uniform mat3 normalMatrix; out vec3 vNormal; out vec3 vViewPos; out vec2 vUV; void main() { vec4 viewPos = modelViewMatrix * vec4(aPos, 1.0); gl_Position = projectionMatrix * viewPos; vNormal = normalMatrix * aNormal; vViewPos = viewPos.xyz; vUV = aUV; }normalMatrix必须由模型视图矩阵的逆转置矩阵计算，尤其在模型做了非均匀缩放时，否则法线方向会被拉斜，光照轮廓会变形。
3.3.2 片元着色器：把材质参数做成uniform或纹理
#version 330 core in vec3 vNormal; in vec3 vViewPos; in vec2 vUV; uniform vec3 lightDirection; // 在视图空间中归一化 uniform vec3 lightColor; uniform vec3 ambientColor; uniform float shininess; // 镜面指数 uniform sampler2D tex; out vec4 fragColor; void main() { vec3 N = normalize(vNormal); vec3 L = normalize(-lightDirection); vec3 V = normalize(-vViewPos); vec3 H = normalize(L + V); vec3 albedo = texture(tex, vUV).rgb; float diff = max(dot(N, L), 0.0); float spec = pow(max(dot(N, H), 0.0), shininess); vec3 color = ambientColor * albedo + diff * lightColor * albedo + spec * lightColor * 0.8; fragColor = vec4(color, 1.0); }shininess是Bezier曲面高光强度的主要旋钮。shininess=32会让高光区分散，适合哑光表面；shininess=128会产生很窄的高光，适合像汽车漆面那样高度抛光的曲面。如果纹理本身亮暗不均匀，可能需要将采样值先做sRGB解码，否则高光计算会在错误颜色空间里进行。
3.4 参数对标：Phong与Blinn-Phong在高光形状上的取舍
很多人担心半程向量H的归一化开销，但在现代GPU上这比反射向量计算更快。下表给出两种方法的对照：
| 光照模型 | 高光条件 | 计算核心 | 适用场景 |
|---|---|---|---|
| Phong | R·V | 反射向量 | 逐顶点Gouraud着色，老旧管线 |
| Blinn-Phong | N·H | 半程向量 | 逐像素着色，Bezier曲面等光滑表面更自然 |
对于Bezier曲面这样法线变化平缓的几何体，Blinn-Phong在高光边缘的过渡更自然，视觉上更容易接近真实镜面反射。因此大多数计算机图形学课程在Bezier曲面练习中都建议直接写Blinn-Phong计算。
4. Bezier曲面的纹理映射：UV生成、纹理过滤与压缩配置
4.1 参数域(u,v)直接当纹理坐标：从曲面公式导出UV映射
Bezier曲面最重要的属性就是参数(u,v)天然适合直接映射为纹理坐标。若纹理尺寸不超过OpenGL最大纹理范围，直接让UV=(u,v)，不需要额外变换。这样的映射在整个曲面上是连续的，不会因为细分步长变化而跳变。
但如果控制点间距相差悬殊，直接用(u,v)会导致纹理拉伸，比如同一段纹理在控制点密集区域被压扁，在稀疏区域被拉长。要解决需要先做弧长参数化，常见做法是把Bezier曲面按等弧长重新采样，再映射到UV。代码增加会比较多，初学者阶段直接用(u,v)即可。
4.2 OpenGL纹理状态配置：Mipmap、过滤模式和纹理压缩方案
纹理的最终观感不只取决于贴图本身，OpenGL采样参数会直接影响曲面边缘的抗闪烁能力。代码：
glGenTextures(1, &texId); glBindTexture(GL_TEXTURE_2D, texId); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR_MIPMAP_LINEAR); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_REPEAT); glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_REPEAT); glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA8, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE, pixels); glGenerateMipmap(GL_TEXTURE_2D);GL_LINEAR_MIPMAP_LINEAR意味着三线性过滤，曲面远离摄像机缩小时不会产生摩尔纹。如果做医学体积渲染之类的数据可视化，纹理坐标经常动态变化，可以把MAG_FILTER设为GL_LINEAR，MIN_FILTER设为GL_NEAREST_MIPMAP_LINEAR来减少模糊副作用。
纹理压缩也值得关注。桌面上可以用GL_COMPRESSED_RGBA_S3TC_DXT5_EXT或ASTC，但需要确认驱动程序支持。更稳妥的做法是加载DDS或KTX文件，用glCompressedTexImage2D上传。压缩纹理不仅省显存，还能在带宽紧张时缓解压力，只是采样返回的颜色已经被近似，高光细节可能变得粗糙。
4.3 在GLSL中把纹理颜色和光照结果组合：差异化权重
在实际渲染中，纹理颜色代表表面的漫反射率，光照漫反射亮度应该与纹理相乘，而不是相加。如果在片元着色器中写albedo + diff*lightColor，纹理的高光被冲淡，曲面看起来像没贴图。我在3.3节的代码中用的是ambient*albedo + diff*lightColor*albedo，这样纹理明暗会随光照变化。若要模拟金属质感，可以单独从另一个纹理通道读取高光贴图，或者用uniform控制镜面反射颜色。
4.4 多块Bezier曲面纹理接缝的通用解法：统一UV空间
单个Bezier曲面往往拼不成一个闭合几何体，多块Patch拼接时，各块使用自己的参数坐标，贴图会在接缝处断开。常见解法是让相邻Patch共享同一套UV：构造全局参数空间，每个Patch的上限正好是下一个Patch的下限。如果做不到，则需要把多个Patch的参数域偏移到一张纹理的不同区块，类似texture atlas。在Qt/OpenGL项目里，我通常把Patch的UV偏移量作为attribute传入顶点，而不在着色器里硬算。
5. Bezier曲面光照与纹理的验证技巧和优化坑位
5.1 用极短路径验证法线方向：绘制顶点法线而不是整个曲面
在调试光照前，先把法线可视化。方法是在VBO中额外生成一段从顶点位置出发、沿法线方向延长0.1个单位的线段，绘制模式改为GL_LINES。如果看到的朝向和曲面弯曲方向相反，就把叉积顺序换成cross(dv, du)。这个技巧比盯着高光判断快很多，一次渲染就能定位所有翻面区域。
5.2 光照+纹理的调试顺序：先单色光照，再叠纹理
调试Bezier曲面光照时不应立即启用真实纹理，而应先用单色棋盘格纹理验证UV连续性。棋盘格上同一条棋盘线在曲面上的间距应该平滑变化，如果出现突跳或扭曲，说明UV生成代码有bug。在光照阶段，把纹理替换成纯白色，只观察漫反射和镜面轮廓，确认高光形状符合Blinn-Phong预期后，再打开真实纹理。
5.3 从GLSL切换到固定管线作为对照基准
如果你还在用OpenGL 2.1环境做实验，可以在固定管线的glEvalMesh2流程里直接启用GL_AUTO_NORMAL，观察高光轮廓，快速验证控制点数据。这个对照意义在于：固定管线的法线计算结果是OpenGL标准实现，如果GLSL渲染与其不一致，优先检查normalMatrix是否传错，再用5.1的法线可视化定位差异。排除这两点后，光照参数和纹理采样顺序就成了最后的怀疑对象。
至此，一个带光照和纹理的Bezier曲面就能在OpenGL中正确显示。后续可以考虑追加置换贴图或动态LOD，但无论怎么扩展，核心都在于把Bezier参数域转换为稳定的顶点数据、法线和UV，这一小步永远没有捷径可走。
本文还有配套的精品资源，点击获取
