1. 从一个画不出来的困境说起如果你曾经尝试用代码去描述一个曲面比如一个简单的球体、一个扭曲的环面或者更复杂的分形结构你大概率经历过这样的场景用传统的三角网格去逼近它面片数量爆炸内存吃紧边缘还带着明显的多边形棱角。你想让曲面光滑就得加密网格你想让网格变密渲染时间就成倍增长。这个矛盾在实时渲染和离线渲染里都让人头疼。Sphere Tracing球体追踪就是在这种背景下进入我视野的。它属于隐式曲面渲染Implicit Surface Rendering的一种核心方法专门解决用数学公式描述曲面然后直接把它画出来这个问题。和传统的先建模成网格再渲染不同隐式曲面从一开始就不存在网格这个概念——曲面是由一个函数定义的比如f(x, y, z) x² y² z² - r²当这个函数值为零时对应的点就在曲面上。这篇文章我想聊的是隐式曲面和距离场Distance Fields到底怎么理解Sphere Tracing 为什么能高效地渲染它们以及在实际写代码时有哪些坑和技巧。适合对图形学有兴趣、写过一些光线追踪或者Shader、但还没系统接触过距离场渲染的开发者。如果你完全没写过渲染代码也没关系我会尽量用生活化的类比把原理讲清楚。先说结论Sphere Tracing 的核心思想非常优雅——它利用距离场的一个关键性质任意点到曲面的距离让光线以安全步长前进永远不会穿过曲面。这个性质让渲染隐式曲面变得既稳定又高效而且代码量少得惊人。我第一次实现它的时候核心循环不到二十行但渲染出来的效果让我盯着屏幕看了很久。2. 隐式曲面与距离场先搞清楚我们在画什么2.1 隐式曲面到底是什么为什么不用网格传统建模流程里一个物体通常被表示成三角网格顶点、边、面。GPU 渲染管线对三角形有原生支持所以这条路走得很顺。但网格有个根本问题它是对曲面的离散近似。你用一个二十面体去近似球体放大之后就能看到棱角你加密到一百万面内存和带宽又受不了。隐式曲面换了个思路不存储几何只存储一个函数。这个函数接受空间中的任意一点返回一个标量值。当这个值等于零时点就在曲面上大于零在外部小于零在内部。比如球体就是f(p) length(p) - r。这个函数叫符号距离函数Signed Distance FunctionSDF如果它返回的绝对值恰好等于点到曲面的最近距离那它就是一个距离场。距离场的好处在于它不只是告诉你在里还是在外还告诉你离表面有多远。这个多远的信息就是 Sphere Tracing 能工作的全部基础。2.2 距离场的核心性质Lipschitz 连续性距离场有一个非常重要的数学性质它的梯度长度不超过 1。用大白话说就是你在空间里移动一小段距离函数值的变化量不会超过你移动的距离。这个性质叫Lipschitz 连续性常数取 1。为什么这个性质重要因为它给了我们一个安全半径的保证。假设当前点p到曲面的距离是d那么以p为球心、d为半径的球体内一定不包含任何曲面。这个结论是严格的不是近似。你可以这样理解如果球内真的有曲面那p到曲面的距离就小于d和距离场的定义矛盾。这个安全球的概念就是 Sphere Tracing 名字的由来——每一步都沿着光线方向前进一个安全球的半径保证不会穿模。2.3 距离场和普通隐式函数的区别这里要区分两个概念。任意隐式函数f(p) 0都能定义曲面但只有满足|f(p)| ≈ dist(p, surface)的函数才是距离场。比如f(p) length(p)² - r²也能定义球体但它不是距离场——在离球很远的地方它的值增长得比真实距离快得多。如果你拿一个非距离场的函数去做 Sphere Tracing步长就会过大光线可能直接穿过曲面渲染结果出现空洞或者错误。所以实际使用中我们几乎总是构造近似距离场哪怕不是精确的也要保证它不会高估距离也就是|f(p)| ≤ dist(p, surface)。这个条件叫保守性是 Sphere Tracing 正确性的前提。注意很多教程直接拿length(p) - r当例子这确实是精确距离场。但一旦开始做布尔运算、变形、重复距离场就变成近似的了这时候保守性就成了必须检查的事情。3. Sphere Tracing 的算法逻辑为什么它能又快又稳3.1 从光线步进说起光线追踪的基本流程是从相机发射一条光线找到它和场景的第一个交点然后着色。对于网格我们通常用射线-三角形求交解析求解。对于隐式曲面没有解析解只能数值逼近。最朴素的数值方法是固定步长步进沿着光线每次前进一个固定的小距离检查当前点的函数值是否接近零。这个方法的问题很明显——步长太小效率极低步长太大直接穿过薄壁结构。而且对于不同尺度的场景步长很难统一。Sphere Tracing 的改进在于步长不是固定的而是由距离场动态决定的。每一步都查询当前点的距离值然后前进这么多。因为距离值保证是到曲面的最近距离所以这一步绝对不会跨过曲面。当距离值小于某个阈值时就认为命中了。3.2 算法伪代码与逐行解释我把核心循环写出来你可以对照着理解float sphereTrace(vec3 ro, vec3 rd) { float t 0.0; // 当前光线行进距离 for (int i 0; i MAX_STEPS; i) { vec3 p ro rd * t; // 当前采样点 float d map(p); // 查询距离场 if (d EPSILON) { // 足够接近曲面认为命中 return t; } t d; // 前进一个安全步长 if (t MAX_DIST) break; // 超出最远距离放弃 } return -1.0; // 未命中 }逐行看t是已经走过的距离p是当前点map(p)返回该点到最近曲面的距离。如果距离小于EPSILON比如 0.001就认为光线打到了曲面。否则前进d继续循环。MAX_STEPS防止死循环MAX_DIST防止光线无限远。这个循环的收敛速度取决于距离场的形状。对于光滑的球体几步就能命中对于复杂的场景可能需要几十步。但相比固定步长它的效率通常高出一个数量级。3.3 为什么不会穿模安全半径的严格保证我再强调一遍这个性质因为它是整个算法的基石。假设当前点p的距离值是d那么以p为球心、d为半径的球内没有曲面。光线从p沿方向rd前进d距离到达pp一定还在这个球内因为rd是单位向量。所以p不可能在曲面上也不可能穿过曲面。这个保证只在距离场保守时成立。如果你用的函数在某些区域高估了距离比如返回了2 * dist那步长就会过大光线可能直接跳到曲面另一侧导致漏检。这就是为什么构造距离场时宁可低估也不能高估。实操心得如果你不确定自己的距离场是否保守可以写一个测试在空间里随机采样大量点对每个点计算|f(p)|和真实最近距离可以用暴力遍历或者已知解析解检查是否始终|f(p)| ≤ dist。这个测试能帮你提前发现很多渲染 bug。4. 距离场的构造从简单图元到复杂场景4.1 基础图元的距离函数构造距离场的第一步是掌握基础图元。下面这张表是我常用的几个直接可以抄图元距离函数说明球体length(p) - r精确距离场立方体length(max(abs(p) - b, 0.0)) min(max(...), 0.0)精确距离场b是半边长平面dot(p, n) dn是法线d是偏移环面length(vec2(length(p.xz) - R, p.y)) - rR是大半径r是小半径圆柱max(length(p.xz) - r, abs(p.y) - h)近似需要修正立方体的公式看起来复杂但拆开就清楚了max(abs(p) - b, 0.0)处理外部min(max(...), 0.0)处理内部。这个公式在外部是精确的内部也是精确的所以整体是精确距离场。4.2 布尔运算并集、交集、差集有了基础图元下一步是把它们组合起来。距离场的布尔运算非常优雅并集min(d1, d2)交集max(d1, d2)差集max(d1, -d2)为什么并集是min因为并集的表面是离任意一个图元最近的地方取最小值正好对应这个逻辑。交集取最大值因为要同时在两个图元内部距离两个表面的最大值决定了边界。差集稍微绕一点d1减去d2相当于在 d1 内部但不在 d2 内部所以是max(d1, -d2)。这些运算的结果不一定是精确距离场但通常是保守的。比如两个球并集min(d1, d2)在交界处会低估距离但不会高估所以 Sphere Tracing 仍然安全。4.3 平滑混合与变形硬边布尔运算做出来的形状比较生硬实际项目里经常需要平滑过渡。最常用的平滑并集是smooth minfloat smin(float a, float b, float k) { float h clamp(0.5 0.5 * (b - a) / k, 0.0, 1.0); return mix(b, a, h) - k * h * (1.0 - h); }k控制混合的柔和程度越大越平滑。这个公式的推导涉及多项式插值但用起来很简单直接替换min就行。类似的还有smax用于平滑交集。变形方面常见的有扭曲twist、弯曲bend、重复repeat。重复特别有用比如做无限延伸的柱子或者阵列vec3 repeat(vec3 p, vec3 c) { return mod(p 0.5 * c, c) - 0.5 * c; }这个操作把空间切成周期性的格子每个格子里放同一个图元就能得到无限阵列。但要注意重复后的距离场不再是精确的因为格子边界处的距离被低估了。不过它仍然是保守的Sphere Tracing 可以正常工作。注意事项变形操作会破坏距离场的精确性尤其是非等距变换比如非均匀缩放。如果你发现渲染结果有瑕疵先检查变形后的距离场是否还保守。一个常用的修正方法是乘以变换的 Jacobian 范数的倒数但具体推导比较繁琐实际项目中更多是靠经验调整。5. 法线计算与着色让画面看起来对5.1 用梯度近似法线Sphere Tracing 找到交点后下一步是着色。着色需要法线而隐式曲面没有显式的法线信息只能通过距离场的梯度来求。梯度可以用中心差分近似vec3 calcNormal(vec3 p) { vec2 e vec2(0.001, 0.0); return normalize(vec3( map(p e.xyy) - map(p - e.xyy), map(p e.yxy) - map(p - e.yxy), map(p e.yyx) - map(p - e.yyx) )); }这里e是采样偏移通常取0.001左右。太小会因为浮点精度产生噪声太大则法线不够精确。这个方法的本质是用有限差分估计偏导数三个分量分别对应 x、y、z 方向的梯度。5.2 光照模型的选择有了法线就可以套用标准的光照模型。最常用的是Lambert 漫反射 Blinn-Phong 高光vec3 shade(vec3 p, vec3 n, vec3 rd) { vec3 lightDir normalize(vec3(1.0, 1.0, 1.0)); float diff max(dot(n, lightDir), 0.0); vec3 halfVec normalize(lightDir - rd); float spec pow(max(dot(n, halfVec), 0.0), 32.0); return vec3(0.2) vec3(0.8) * diff vec3(1.0) * spec; }这个模型简单但效果不错适合快速验证。如果想要更真实的效果可以加环境光遮蔽AO、软阴影、菲涅尔反射等。AO 在距离场渲染里特别容易实现——沿着法线方向采样几个点累加距离值就能得到近似遮蔽float calcAO(vec3 p, vec3 n) { float occ 0.0; float sca 1.0; for (int i 0; i 5; i) { float h 0.01 0.12 * float(i) / 4.0; float d map(p n * h); occ (h - d) * sca; sca * 0.95; } return clamp(1.0 - 3.0 * occ, 0.0, 1.0); }这段代码的逻辑是如果法线方向附近有几何体距离值会小于采样高度累加后得到遮蔽。sca是衰减系数让近处的贡献更大。5.3 软阴影的实现距离场渲染的另一个优势是软阴影很容易做。核心思路是从交点向光源方向步进记录最小距离和行进距离的比值比值越小阴影越硬float softShadow(vec3 ro, vec3 rd, float mint, float maxt, float k) { float res 1.0; float t mint; for (int i 0; i 32; i) { float h map(ro rd * t); res min(res, k * h / t); t clamp(h, 0.01, 0.2); if (res 0.005 || t maxt) break; } return clamp(res, 0.0, 1.0); }k控制阴影的柔和程度越大越硬。这个方法的原理是如果光线路径上某点的距离值很小说明它靠近遮挡物阴影就暗。h / t这个比值衡量了遮挡物张角比值越小阴影越深。实操心得软阴影的采样次数和步长需要平衡。步长太小性能差步长太大阴影会有条带。我通常用clamp(h, 0.01, 0.2)作为步长采样 32 次效果和性能比较均衡。如果场景简单可以降到 16 次。6. 性能优化与常见坑6.1 步数、精度与性能的三角关系Sphere Tracing 的性能主要取决于三个参数MAX_STEPS、EPSILON、MAX_DIST。它们之间是互相制约的MAX_STEPS越大能处理的复杂场景越多但每帧计算量越大。EPSILON越小命中判定越精确但需要的步数越多。MAX_DIST越大能看到越远的物体但远处物体的步进效率低。我的经验值是MAX_STEPS 128EPSILON 0.001MAX_DIST 100.0。这个配置能覆盖大多数场景。如果场景特别复杂可以提高到 256 步如果追求性能可以降到 64 步但可能出现漏检。6.2 距离场不保守导致的穿模这是最常见的 bug。表现是某些角度下物体表面出现空洞或者光线直接穿过薄壁。原因通常是某个变形操作破坏了保守性。排查方法是把距离场可视化用颜色表示距离值检查是否有区域的值大于真实距离。一个实用的技巧是在距离场函数里加一个全局缩放系数比如return 0.9 * d;。这会人为地低估距离保证保守性代价是步数增加。如果加了系数之后穿模消失就说明原来的距离场确实高估了。6.3 浮点精度问题在距离场渲染里浮点精度问题比网格渲染更明显。因为步进是累积的误差会放大。常见表现是远处物体表面出现噪点或者法线计算不稳定。解决方法有几个一是用highp精度在 GLSL 里显式声明二是把场景缩放到合理范围避免极大或极小的坐标三是在法线计算时用相对偏移比如e 0.001 * length(p)而不是固定值。6.4 常见问题速查表问题现象可能原因解决方法表面有空洞距离场不保守加缩放系数检查变形操作边缘有噪点浮点精度不足提高精度调整 EPSILON渲染太慢步数过多降低 MAX_STEPS优化距离场阴影有条带软阴影步长过大减小步长增加采样次数法线方向错误差分偏移过大减小 e检查距离场符号远处物体消失MAX_DIST 太小增大 MAX_DIST或加 LOD7. 从 Demo 到项目我的实操路线7.1 最小可运行示例如果你想快速上手我建议从 Shadertoy 风格的片段着色器开始。核心代码不到 50 行就能渲染一个带光照的球体。先把map函数写成最简单的球体跑通 Sphere Tracing 循环确认能出图。然后逐步加图元、加布尔运算、加光照。这个渐进路线能帮你快速建立直觉。7.2 进阶方向跑通基础之后可以往几个方向深入。一是分形比如 Mandelbulb用距离场渲染分形特别合适因为分形没有网格表示。二是体积渲染距离场可以扩展到体积密度场做云、烟、雾。三是实时性能优化比如用锥体追踪Cone Tracing做 LOD或者用时间重投影减少每帧计算量。7.3 我踩过的几个坑第一个坑是距离场符号搞反。差集运算max(d1, -d2)里的负号很容易漏导致结果完全错误。第二个坑是变形后忘记修正距离。比如非均匀缩放之后距离场不再精确需要乘以缩放因子的倒数。第三个坑是软阴影的 k 值调不好要么全黑要么全亮后来发现是h / t的比值范围没搞对。最后一个心得距离场渲染的调试可视化比断点更有效。把距离值映射成颜色把步数映射成亮度一眼就能看出问题在哪。这个习惯帮我省了大量时间。8. 一些个人体会Sphere Tracing 最吸引我的地方是它把渲染这件事还原成了纯粹的数学。没有网格、没有加速结构、没有复杂的求交代码只有一个函数和一个循环。这种简洁性让它在快速原型、艺术创作、教育演示里特别有优势。当然它也不是万能的。对于大规模静态场景网格加 BVH 仍然更快对于需要精确物理碰撞的场景距离场也有局限。但如果你要渲染的是数学定义的形状、分形、或者需要频繁变形的物体Sphere Tracing 几乎是最优解。我现在的习惯是任何可以用公式描述的形状优先考虑距离场。写起来快改起来也快调参直观效果可控。这个工具值得每个图形学爱好者花时间掌握。