OpenGL实现球体绘制函数运行无显示 代码错误排查求助

存在的核心问题 #

• 数组越界：声明的顶点数组尺寸为30*30，合法索引范围是0~29，但绘制循环初始值设为i=30、j=30，还访问了a[i+1][j]、a[i][j+1]这类超出数组长度的位置，触发未定义行为，顶点数据完全错乱。
• 顶点数据填充不全：生成顶点的循环只覆盖了部分索引位置，数组内大量边缘位置、0索引位置的数值是未初始化的内存垃圾值，根本不是有效球面坐标。
• 球面坐标公式错误：极角、方位角对应的三角函数调用和遍历逻辑错误，镜像点生成规则有误，无法生成闭合的完整球面顶点。
• 三角面索引逻辑错误：拆分矩形为三角形时的顶点顺序和索引对应关系混乱，无法拼接成连续的球面网格。

修正后可正常运行的代码 #

#define TAU 6.28318530718    

// 球体绘制函数：r为球体半径，x/y/z为球心在场景中的坐标
void spherel(double r, double x=0, double y=0, double z=0)
{
    const int h_seg = 30; // 水平方向分段数
    const int v_seg = 30; // 垂直方向分段数
    // 数组多分配1位长度，避免相邻顶点访问越界
    double a[h_seg+1][v_seg+1][3];
    double h_step = TAU / h_seg;
    double v_step = TAU / (2 * v_seg);

    // 按标准球面坐标公式生成所有网格顶点
    for(int i = 0; i <= h_seg; i++){
        double phi = i * h_step;
        for(int j = 0; j <= v_seg; j++){
            double theta = -TAU/4 + j * v_step; // 极角从南极遍历到北极
            a[i][j][0] = r * cos(theta) * cos(phi);
            a[i][j][1] = r * sin(theta);
            a[i][j][2] = r * cos(theta) * sin(phi);
        }
    }

    // 按三角面图元绘制网格
    glBegin(GL_TRIANGLES);
    for(int i = 0; i < h_seg; i++){
        for(int j = 0; j < v_seg; j++){
            // 每个网格矩形拆为2个三角形
            glVertex3d(x+a[i][j][0], y+a[i][j][1], z+a[i][j][2]);
            glVertex3d(x+a[i+1][j][0], y+a[i+1][j][1], z+a[i+1][j][2]);
            glVertex3d(x+a[i][j+1][0], y+a[i][j+1][1], z+a[i][j+1][2]);

            glVertex3d(x+a[i+1][j][0], y+a[i+1][j][1], z+a[i+1][j][2]);
            glVertex3d(x+a[i+1][j+1][0], y+a[i+1][j+1][1], z+a[i+1][j+1][2]);
            glVertex3d(x+a[i][j+1][0], y+a[i][j+1][1], z+a[i][j+1][2]);
        }
    }
    glEnd();
}

调用注意事项 #

• 调用函数前需正确配置OpenGL视口、投影矩阵、模型视图矩阵，确保球体坐标落在视锥体裁剪范围内。
• 渲染前开启深度测试glEnable(GL_DEPTH_TEST)，避免三角面绘制顺序错误导致的显示异常。
• 如果需要开启光照，可直接将归一化后的顶点坐标作为法向量传入，无需额外计算。

内容的提问来源于stack exchange，提问作者S.M. Abdul Fattah