OpenGL细分正二十面体出现T型接头的原因及调试方法

我尝试用OpenGL细分着色器对预细分的正二十面体做曲面细分，目标是让靠近相机的三角面获得更多细分，同时尽可能避免T型接头。我在视图空间计算每条边的中心和原始三角面的中心，用位置的Z分量作为距离，且认为同一条边在不同三角面的距离应该一致，但现在原始三角面的边上出现了T型接头——用equal_spacing时特别明显；用fractional_odd_spacing时接头变小，但叠加高度图后问题依然存在。

我怀疑的原因：

• 距离计算正确，但用错了边的细分层级（gl_TessLevelOuter索引错误）
• 对OpenGL细分的核心概念理解不到位
• 计算位置时的浮点精度限制
• 其他未考虑到的因素

如果我的方案逻辑没问题，该怎么调试这类问题？

顶点着色器 #

#version 430 core

layout (location = 0) in vec3 inPos;

uniform mat4 _modelMat;
uniform mat4 _viewMat;

out vec3 csPos;
out vec3 csVPos;

void main() {
    csPos = inPos;
    csVPos = vec3(_viewMat * _modelMat * vec4(inPos, 1.0));
}

细分控制着色器 #

#version 430 core

layout (vertices = 3) out;

in vec3 csPos[];
in vec3 csVPos[];

out vec3 esPos[];

uniform mat4 _modelMat;
uniform mat4 _viewMat;

void main() {
    esPos[gl_InvocationID] = csPos[gl_InvocationID];

    if(gl_InvocationID == 0) {
        float camDist;
        // 控制靠近相机时的细分速度
        const float tessFactor = 1.0;

        // 计算外部细分层级
        for(int vertexInd = 0; vertexInd < 3; vertexInd += 1) {
            vec3 edgePos = 0.5 * (csVPos[vertexInd] + csVPos[(vertexInd+1)%3]);
            camDist = -edgePos.z;

            gl_TessLevelOuter[vertexInd] = tessFactor / camDist;
        }

        // 计算内部细分层级
        vec3 triangleCenter = 0.33333 * (csVPos[0] + csVPos[1] + csVPos[2]);
        camDist = -triangleCenter.z;
        gl_TessLevelInner[0] = tessFactor / camDist;
    }
}

细分评估着色器 #

#version 430 core

layout (triangles, fractional_odd_spacing, ccw) in;

in vec3 esPos[];

void main() {
    vec3 vertexPos = vec3(0.0);
    vertexPos += esPos[0] * gl_TessCoord.x;
    vertexPos += esPos[1] * gl_TessCoord.y;
    vertexPos += esPos[2] * gl_TessCoord.z;
    vertexPos = normalize(vertexPos);

    gl_Position = vec4(vertexPos, 1.0);
}

几何着色器 #

#version 430 core

layout (triangles) in;
layout (triangle_strip, max_vertices = 3) out;

uniform mat4 _modelMat;
uniform mat4 _viewMat;
uniform mat4 _projMat;

void main() {
    for(int vertexInd = 0; vertexInd < 3; vertexInd += 1) {
        gl_Position = _projMat * _viewMat * _modelMat * gl_in[vertexInd].gl_Position;
        EmitVertex();
    }

    EndPrimitive();
}

调试步骤 #

1. 验证gl_TessLevelOuter索引匹配

   • 确保相邻三角面计算同一条边时，顶点顺序的反向不影响边中心的计算结果，且对应的gl_TessLevelOuter索引指向正确的边。比如三角面A的边是顶点0→1，相邻三角面B的同一条边是顶点1→0，此时两边计算出的边中心必须完全一致，细分层级值也必须相同。
   • 可以在细分控制着色器中添加输出变量，把每条边的细分层级值传递到后续着色器，最终渲染为颜色，直观查看相邻边的细分因子是否一致。

2. 强制统一共享边的细分层级

   • 临时将所有gl_TessLevelOuter设为固定值（比如4），如果T型接头消失，说明问题根源是相邻面细分层级不匹配。
   • 尝试改用共享顶点到相机的距离平均值计算边的细分层级，因为顶点是共享的，相邻面计算同一条边的两个顶点距离的平均值必然一致，避免边中心计算可能带来的差异。

3. 排查浮点精度问题

   • 在细分控制着色器中把计算出的camDist和细分层级值输出到CPU端，对比相邻三角面同一条边的数值是否存在微小浮点误差。
   • 尝试将计算用的float替换为double，提升计算精度，看是否能减少因浮点运算导致的细分层级差异。

4. 修正坐标变换流程

   • 当前细分评估着色器输出世界空间位置，几何着色器再做MVP变换，可能引入额外误差。建议在细分评估着色器中直接完成完整的MVP变换，确保所有顶点的变换流程一致。
   • 检查normalize(vertexPos)的时机：如果在细分后归一化，相邻面的共享边顶点可能因细分计算差异导致归一化后位置偏移。可以把归一化移到顶点着色器，确保原始顶点都是单位向量。

5. 可视化细分网格拓扑

   • 使用RenderDoc等调试工具捕获帧，查看细分后的网格结构，确认T型接头的位置和对应的细分层级差异。
   • 给每个原始三角面分配唯一颜色，直观查看相邻面的细分边界是否对齐。

内容的提问来源于stack exchange，提问作者Kacper Gruszka