修复 Go-GL 中 mgl32.Project 方法的坐标投影错误-Golang-PHP中文网

修复 Go-GL 中 mgl32.Project 方法的坐标投影错误

Go-GL 的 mgl32.Project 方法因缺失透视除法（perspective divide）导致屏幕坐标计算严重失真，本文详解问题成因、验证方式及可靠替代方案。

go-gl 的 `mgl32.project` 方法因缺失透视除法（perspective divide）导致屏幕坐标计算严重失真，本文详解问题成因、验证方式及可靠替代方案。

在使用 Go-GL 进行 OpenGL 坐标变换时，mgl32.Project 是将世界坐标（或裁剪空间前坐标）映射到窗口坐标的常用工具。然而，如实际测试所示——当对点 (5, 0, 0) 执行投影时，得到 X ≈ 4021.32（远超 800px 宽度），这明显违背几何直觉，也与实际渲染结果（如位于 (5,0,0) 的立方体完全可见）矛盾。

根本原因在于 mgl32.Project 的实现缺陷：它直接对齐次坐标 vpp = projection × modelview × [x,y,z,1] 的前三分量进行视口变换，却遗漏了关键的透视除法步骤（即除以 vpp[3]，也就是 w 分量）。标准 OpenGL 投影流程必须包含：

顶点乘以 MVP 矩阵 → 得到齐次裁剪坐标 (xₚ, yₚ, zₚ, wₚ)
执行透视除法：(xₚ/wₚ, yₚ/wₚ, zₚ/wₚ) → 归一化设备坐标（NDC，范围 [-1,1]）
映射至窗口坐标：x_win = (x_ndc + 1) × width/2 + x₀，同理 y_win

而当前 mgl32.Project 跳过了第 2 步，直接用未归一化的 xₚ, yₚ 计算窗口坐标，导致结果被 wₚ 放大或扭曲。

以下为可复现问题的最小验证代码：

// 错误示例：原生 Project 方法（输出异常）
origin := mgl32.Vec3{0, 0, 0}
screenOrigin := mgl32.Project(origin, modelView, projection, 0, 0, 800, 600)
fmt.Printf("Origin (broken): (%.2f, %.2f)\n", screenOrigin[0], screenOrigin[1]) // → (400.00, 300.00)

test := mgl32.Vec3{5, 0, 0}
screenTest := mgl32.Project(test, modelView, projection, 0, 0, 800, 600)
fmt.Printf("Test (broken):   (%.2f, %.2f)\n", screenTest[0], screenTest[1]) // → (4021.32, 300.00)

✅ 正确解决方案：手动实现带透视除法的投影

func ProjectFixed(obj mgl32.Vec3, modelView, projection mgl32.Mat4, x, y, width, height float32) [3]float32 {
    // 1. 构造齐次向量并应用 MVP 变换
    obj4 := obj.Vec4(1.0)
    vpp := projection.Mul4(modelView).Mul4x1(obj4) // vpp = [x,y,z,w]

    // 2. ✅ 关键：透视除法 → 转换为 NDC
    if vpp[3] == 0 {
        panic("w component is zero — cannot project")
    }
    ndcX := vpp[0] / vpp[3]
    ndcY := vpp[1] / vpp[3]
    ndcZ := vpp[2] / vpp[3]

    // 3. NDC → 窗口坐标（OpenGL 默认 Y 向上，窗口 Y 向下需翻转）
    winX := float32(x) + (ndcX+1.0)*float32(width)/2.0
    winY := float32(y) + (1.0-ndcY)*float32(height)/2.0 // 注意：OpenGL NDC Y↑，窗口 Y↓，故用 (1-ndcY)
    winZ := (ndcZ + 1.0) / 2.0

    return [3]float32{winX, winY, winZ}
}

// 使用示例
screenOriginFixed := ProjectFixed(origin, modelView, projection, 0, 0, 800, 600)
screenTestFixed := ProjectFixed(test, modelView, projection, 0, 0, 800, 600)
fmt.Printf("Origin (fixed): (%.2f, %.2f)\n", screenOriginFixed[0], screenOriginFixed[1]) // ≈ (400.00, 300.00)
fmt.Printf("Test (fixed):   (%.2f, %.2f)\n", screenTestFixed[0], screenTestFixed[1])   // ≈ (625.00, 300.00) — 符合预期！

? 注意事项：