请教Pixar OpenSubdiv中OsdComputePostProjectionSphereExtent函数的计算逻辑

问题背景 #

我在浏览Pixar的OpenSubdiv源码时，遇到一个无法理解的简短函数：

float OsdComputePostProjectionSphereExtent(vec3 center, float diameter)
{
    vec4 p = OsdProjectionMatrix() * vec4(center, 1.0);
    return abs(diameter * OsdProjectionMatrix()[1][1] / p.w);
}

已知信息：

• OsdProjectionMatrix()仅返回当前帧的4x4投影矩阵
• 参数center是相机空间中球体的中心，diameter是球体直径
• 从函数名推测，它返回球体经投影矩阵投影到屏幕空间后的某种“范围”

我搜索相关关键词未找到有效信息，已知球体投影后通常是椭圆，想确认：

• 这个函数是计算投影后椭圆长轴的长度，还是估算值？
• 如果是估算值，是否为保守估算？
• 或者它计算的是完全不同的内容？

函数解析 #

核心作用 #

这个函数是对球体投影后在屏幕垂直方向（Y轴）上的范围的快速近似估算，并非精确计算椭圆的长轴长度。

逐行拆解逻辑 #

1. 透视投影变换：相机空间的球心center被投影矩阵转换为齐次坐标p，其中p.w是透视除法的关键分母——在标准透视投影下，p.w等于-center.z（投影矩阵第四行通常为(0,0,-1,0)）。
2. 垂直方向缩放系数：投影矩阵的[1][1]元素对应垂直方向的缩放比例，标准透视投影中这个值是2/(tan(fovy/2))，作用是把相机空间的Y轴范围映射到NDC（标准化设备坐标）的[-1,1]区间。
3. 投影范围计算：
   • 相机空间里球体的垂直方向总跨度是diameter，乘以投影矩阵的垂直缩放系数，得到未做透视除法前的垂直范围。
   • 除以p.w完成透视除法，将结果转换到NDC空间——这一步体现了透视投影“近大远小”的特性：球体离相机越远，p.w的绝对值越大，最终投影范围就越小。

估算性质与保守性 #

• 属于近似估算：它只聚焦于垂直方向的投影范围，完全没考虑水平方向的投影结果，也忽略了球心不在相机光轴上时，透视变形对椭圆长轴的影响。
• 是垂直方向的保守估算：该计算结果不会小于球体在屏幕垂直方向的实际投影范围，能保证后续基于此的逻辑（比如曲面细分密度）不会低估所需精度。

实际应用场景 #

在OpenSubdiv的曲面细分管线中，这个函数是用来快速计算曲面片包围球的屏幕投影大小，以此决定曲面细分的级别——投影范围越大，细分越精细，确保曲面在屏幕上的显示精度达标。

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