次表面散射指的是光线射入半透明材质，在内部发生散射后再透射出来的光线传播过程，考虑到有些项目会需要使用次表面散射，下面就给大家介绍下在Unity3D中次表面散射的简单实现，希望可以帮到大家。
一、前言
本文旨在与大家一起探讨学习新知识，如有疏漏或者谬误，请大家不吝指出。
以下内容参考了GPU精粹1中关于次表面散射一节的内容。
二、概述
次表面散射，英文全称为Subsurface Scattering，简称SSS。指的是光线射入半透明材质，在内部发生散射后再透射出来的光线传播过程。真实世界中拥有次表面散射的材质有蜡烛、大理石、玉石以及人的皮肤等。要模拟物理上真实的次表面散射是很复杂的一件事，比较经典的次表面散射模型有BSSRDF，全称叫双向次表面散射反射分布函数。本文中并未使用BSSRDF模型，而是简单的使用了exp指数函数对深度进行计算以此模拟散射。
以下为实现后的最终效果：
三、原理
首先我们来对次表面散射物体的光照做一个分解，也算是一个简单的建模过程。
Color = Diffuse * Scattering + SpecularColor;
Diffuse指的是物体表面的漫反射颜色，Scattering是散射的颜色，SpecularColor是高光颜色。
1、漫反射的计算
漫反射的计算公式，在这里我们使用的是环绕光照计算公式，即：
?
| 1 |
diff = ( dot(normal, lightDir) + wrap )/( 1+wrap )
|
其中normal为顶点的法线，lightDir是光照方向，wrap为环绕参数。传统的Lambert光照中，当物体表面的法线与光源方向垂直的时候，其产生的光照结果为0，但是在次表面散射物体中，由于内部散射光线的传播，导致其在上述情形下光照结果不会完全为0。所以为了减少传统Lambert光照中的黑暗区域，我们使用环绕光照公式，当dot(normal, lightDir)结果为0的时候，我们强制其至少有一点点亮度，即wrap/(1+wrap)。
2、散射的计算
我们假设物体表面的散射是均匀分布的，并且无视光源位置以及光照方向对散射的影响，取物体在视线方向上的深度值作为参数，带入exp指数函数中进行计算。当然上述假设并不符合物理规律，但是考虑到效果以及效率的问题，我们只好先这么干了。
为了获取物体在视线方向上的深度值，我们需要先以cull front模式渲染一遍物体，保存物体背面的顶点的深度值信息，然后再回到正常的cull back模式下渲染物体，使用(backDepth-frontDepth)来求出深度值，最后带入公式exp(-C*depth)中。C为外部传入的参数，用于调节物体的透光率。
3、高光的计算
高光的计算我们使用经典的BlinnPhong光照公式，即：
?
| 1 2 |
Specular = pow((dot(normal, half)), shiness);
Half = normalize(lightDir + viewDir);
|
其中normal为顶点法线向量，half为半角向量，是入射光向量与视点向量的角平分线向量，shiness为高光指数。
BlinnPhong光照模型相比较于Phong光照模型，其高光区域更平滑柔和，这也是为什么我们使用它。
4、半透明
由于在散射计算中，需要使用到物体表面顶点的深度值信息，导致我们在渲染时不 能关闭ZWrite，这就使得我们不能通过Unity3D中设置RenderType=Transparent、Queue=Transparent来实现半透明混合效果。在Unity3D中，要实现半透明，一般的做法是：
?| 1 2 3 |
Tags { "RenderType"="Transparent" "Queue"="Transparent"}
Blend SrcAlpha OneMinusSrcAlpha
ZWrite off
|
所以，我们需要寻找另外一种方法来实现半透明，通过GrabPass操作来获取除物体本身以外的屏幕渲染结果，然后我们在片段着色器中手动进行混合计算，以此达到半透明效果。当然需要注意的是，GrabPass本身的操作比Alpha混合要昂贵的多，需要牺牲更多的计算性能，另外GrabPass在某些手机平台上可能不被支持。
四、实现
下面给出顶点着色器以及片段着色器的实现代码：
?| 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 |
struct appdata
{
float4 vertex : POSITION;
float2 uv : TEXCOORD0;
float3 normal:NORMAL;
};
struct v2f
{
float2 uv : TEXCOORD0;
float4 screenUV:TEXCOORD1;
float3 lightDir:TEXCOORD2;
float3 viewDir:TEXCOORD3;
float3 normal:TEXCOORD4;
float4 grabUV:TEXCOORD5;
UNITY_FOG_COORDS(1)
float4 vertex : SV_POSITION;
};
sampler2D _MainTex;
float4 _MainTex_ST;
sampler2D _GrabTexture;
sampler2D _BackDepthTex;
float4 _AttenuationC;
float4 _Color;
float _Shininess;
float _ScatteringFactor;
float _Wrap;
struct LightingInput {
float3 Albedo;
float3 Normal;
float Gloss;
float Specular;
float Alpha;
};
float4 CalculateLighting (LightingInput i, float3 lightDir, float3 viewDir, float atten, float3 scattering)
{
float3 h = normalize (lightDir + viewDir);
float diff = (dot (i.Normal, lightDir)+_Wrap)/(1+_Wrap);
diff = saturate (diff);
float nh = (dot (i.Normal, h)+_Wrap)/(1+_Wrap);
nh = saturate(nh);
float spec = pow (nh, i.Specular*128.0) * i.Gloss;
float4 c;
c.rgb = (i.Albedo * _LightColor0.rgb * diff *scattering + _LightColor0.rgb * _SpecColor.rgb * spec) * (atten * 2);
c.a = i.Alpha + _SpecColor.a * spec * atten;
return c;
}
v2f vert (appdata v)
{
v2f o;
o.vertex = mul(UNITY_MATRIX_MVP, v.vertex);
o.uv = TRANSFORM_TEX(v.uv, _MainTex);
o.screenUV = ComputeScreenPos(o.vertex);
o.lightDir = ObjSpaceLightDir(v.vertex);
o.viewDir = ObjSpaceViewDir(v.vertex);
o.normal = v.normal;
o.grabUV = ComputeGrabScreenPos(o.vertex);
UNITY_TRANSFER_FOG(o,o.vertex);
return o;
}
float4 frag (v2f i) : SV_Target
{
// sample the texture
float4 col = tex2D(_MainTex, i.uv);
//
float frontDepth = LinearEyeDepth( i.screenUV.z/i.screenUV.w );
//
float2 backDepthUV = i.screenUV.xy/i.screenUV.w;
float4 backDepthColor = tex2D(_BackDepthTex, backDepthUV);
float backDepth = LinearEyeDepth(DecodeFloatRGBA(backDepthColor));
//do scattering
float depth = backDepth-frontDepth;
float3 scattering = exp(-_AttenuationC.xyz*depth);
//do lighting
LightingInput lightVar;
lightVar.Albedo = col.rgb * _Color.rgb;
lightVar.Gloss = col.a;
lightVar.Alpha = col.a * _Color.a;
lightVar.Specular = _Shininess;
lightVar.Normal = i.normal;
col = CalculateLighting (lightVar, i.lightDir, i.viewDir, _LightColor0.a, scattering);
//blend
//col.xyz = col.a*col.rgb + (1-col.a)*tex2D(_GrabTexture, i.grabUV.xy/i.grabUV.w);
col.xyz = lerp(tex2D(_GrabTexture, i.grabUV.xy/i.grabUV.w), col.rgb, col.a);
// apply fog
UNITY_APPLY_FOG(i.fogCoord, col);
return col;
}
|
照例，我们对上述代码中某些函数进行说明：
1、ComputeScreenPos函数是将经过透视投影的顶点变换到屏幕坐标系中，然后就可以使用xy/w的值作为UV取屏幕坐标系下的深度图的值。
具体细节可以参看UnityCG.cginc文件，这里也将代码贴出来：
?| 1 2 3 4 5 6 7 8 9 10 11 |
inline float4 ComputeScreenPos (float4 pos) {
float4 o = pos * 0.5f;
#if defined(UNITY_HALF_TEXEL_OFFSET)
o.xy = float2(o.x, o.y*_ProjectionParams.x) + o.w * _ScreenParams.zw;
#else
o.xy = float2(o.x, o.y*_ProjectionParams.x) + o.w;
#endif
o.zw = pos.zw;
return o;
}
|
2、ComputeGrabScreenPos函数做的事情跟上述ComputeScreenPos函数是一样的，只不过对于GrabPass取到的渲染结果与屏幕空间不太一致，这里也列出代码：
?| 1 2 3 4 5 6 7 8 9 10 11 |
inline float4 ComputeGrabScreenPos (float4 pos) {
#if UNITY_UV_STARTS_AT_TOP
float scale = -1.0;
#else
float scale = 1.0;
#endif
float4 o = pos * 0.5f;
o.xy = float2(o.x, o.y*scale) + o.w;
o.zw = pos.zw;
return o;
}
|
3、LinearEyeDepth函数是将经过透视投影变换的深度值还原成其在View坐标系中的值。具体细节读者可以参考此链接：
http://blog.sina.com.cn/s/blog_70f96aa90102v0wd.html
这篇博客是本人早期写的，大体意思有说明到，如果要穷究细节，则需要对透视投影做深入了解了。
4、DecodeFloatRGBA与EncodeFloatRGBA是一对函数。EncodeFloatRGBA用于将float值编码到RGBA四个通道上，Decode则是相应的解码过程。这两个函数是为了提高深度值的精度，以便于进行深度值计算时不会产生太大误差。
5、_BackDepthTex是物体背面的深度图，由Camera.RenderWithShader（）产生，这部分的代码在脚本中实现，读者可以参考完整示例。
下面给出完整示例程序：