XDRender_LightMode_Anisotropic(2) 各项异性着色(2) 切线
一、为何是切线
还是这个图, 需要更深的理解.
为何需要切线来做?
其实我们去看一个图和一个算法能更清楚, 这里就需要说明什么叫各项异性的异性来源. 可以理解为当不同的观察角度,会有不同的影响因子. (有点薛定谔了 哈哈).
观察角度
这里观察角度就是ViewDir
影响因子
得出的影响因子,其实还是法线. 这是这时候的法线由于我们的观察角度被改变了.
其实这样做的原因,还是由于精度和性能, 可以详细非常细非常细的丝线, 环绕的法线就有一圈. 这时候我们不可能完全物理来计算. 所以做了一个近似. 只需要有个切线(变化了一丢丢的)加上观察角度. 就能算出一个大概.
下面是一个还不够细的线, 想象我们再进入微观一点.
细看这个不完全符合物理的
float Aniso_StrandSpecular(float3 T, float3 V, float L, float exponent)
{
float3 H = normalize(L + V);
float dotTH = dot(T, H);
float sinTH = sqrt(1.0 - dotTH*dotTH);
float dirAtten = smoothstep(-1.0, 0.0, dot(T, H));
return dirAtten * pow(sinTH, exponent);
}
二、切线、法线到切线, 切线和方向影响到法线
1、法线到切线的一种步骤
法线纹理
在切线空间
法线
根据计算出的切线空间的法线,经过空间转换,得到世界空间
切线
这时候(世界空间下)根据法线和副切线 得到切线
(世界空间)扰动了切线
法线加重对切线的影响
inputData.tangentWS = ShiftTangent(tangentWS,input.normalWS,_CustomData.z);
影响计算
几种经验:
dot(N,H)^specularity
其中切线方向可以用外部纹理等给到
fixed HdotA = dot(normalize(Normal + AnisoDirection), halfVector);
float aniso = max(0, sin(radians((HdotA + _AnisoOffset) * 180f)));
float spec = saturate(pow(aniso, Gloss * 128) );
sin(T,H)^specularity
1、Kajiya-Kay
2、 circle各向异性, 对Kajiya-Kay一种近似
ward各向异性
这里主要是有粗糙度对两个方向的影响
第一张
ward各向异性分布函数的优点是可以使用comb贴图，也就是各向异性方向贴图，有作为tangent方向矫正各向异性高光形状及方向的功能
if (dotLN < 0.0) // light source on the wrong side?
{
specularReflection = vec3(0.0, 0.0, 0.0);
// no specular reflection
}
else // light source on the right side
{
float dotHN = dot(halfwayVector, normalDirection);
float dotVN = dot(viewDirection, normalDirection);
float dotHTAlphaX =
dot(halfwayVector, tangentDirection) / _AlphaX;
float dotHBAlphaY = dot(halfwayVector,
binormalDirection) / _AlphaY;
specularReflection = attenuation * vec3(_SpecColor)
* sqrt(max(0.0, dotLN / dotVN))
* exp(-2.0 * (dotHTAlphaX * dotHTAlphaX
+ dotHBAlphaY * dotHBAlphaY) / (1.0 + dotHN));
}
}
迪斯尼的G和D两项
// specular
float aspect = sqrt(1-anisotropic*.9);
float ax = max(.001, sqr(roughness)/aspect);
float ay = max(.001, sqr(roughness)*aspect);
float D_GTR2_aniso(float dotHX, float dotHY, float dotNH, float ax, float ay)
{
float deno = dotHX * dotHX / (ax * ax) + dotHY * dotHY / (ay * ay) + dotNH * dotNH;
return 1.0 / (PI * ax * ay * deno * deno);
}
/ Smith GGX G项，各项异性版本
// Derived G function for GGX
float smithG_GGX_aniso(float dotVN, float dotVX, float dotVY, float ax, float ay)
{
return 1.0 / (dotVN + sqrt(pow(dotVX * ax, 2.0) + pow(dotVY * ay, 2.0) + pow(dotVN, 2.0)));