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) — 符合预期！? 注意事项：
- winY 计算中使用 1.0 - ndcY 是因为 OpenGL 的 NDC 坐标系中 Y=1 在顶部，而窗口坐标系中 Y=0 在顶部，需做垂直翻转；
- 务必检查 w ≠ 0，避免除零崩溃（尤其当点位于摄像机后方或投影平面时）；
- 该问题已存在于 mathgl 仓库长期未修复（截至 v0.2.0），建议在生产项目中弃用 mgl32.Project，统一采用自定义安全实现；
- 若需深度值（winZ），同样需经 w 归一化，否则 z 值将失真，影响深度测试一致性。
总结：mgl32.Project 的设计缺陷暴露了底层数学实现与 OpenGL 标准流程的偏差。开发者不应依赖有缺陷的封装，而应理解透视投影的本质——齐次坐标变换 + 透视除法 + 视口映射。通过手动补全透视除法，不仅能获得精确坐标，还能增强对图形管线的理解与可控性。
