Learning Forum

[Mu]

您好：

最近我在参考BSDF例程计算双向反射分布函数（BRDF）时遇到一些问题。例程的simple.lsf脚本中首先使用up面接收上方散射强度，然后用gratingpolar获得不同角度的散射的远场E^2，再将镜面反射去掉，接着对各个方位角使用farfield3dintegrate获得一定角度内的E^2积分，并按照下面的方式求BRDF，最后将结果归一化到生成粗糙面的理论值。

for(i=1:(length(phi)-1)) {

result = result + farfield3dintegrate(E2*cos_theta,ux,uy,5,theta,phi(i)) / (farfield3dintegrate(0*E2+1,ux,uy,5,theta,phi(i))+1e-20); #求BRDF

}

 

我的问题是：

1. 在这个例程的脚本代码中，有一行获取gratingpolar之后的E2 = (BSDF_up.E2) *(1+transmission("BSDF::up"));为什么还需要多加一个透射的系数?，这里导致最后的total_power_R加起来大于1
2. 按照常见的BRDF定义，BRDF = dP_s / (dΩ_s * Pi * cosθ_s)，如何理解代码中BRDF的计算方式，它是如何定义的？分母的(farfield3dintegrate(0*E2+1,ux,uy,5,theta,phi(i))+1e-20），它的意思应该是将E^2=1的强度沿ux,uy在5°内积分，得到的是5°对应的立体角。我理解的代码的定义是BRDF（θ_s） = Es^2 *cosθ_s / （Ω_s*1^2)，这与常见的定义不同，另外Ps/Pi 与Es^2 / Ei^2是否相等？
3. 我如何准确获取这一个5°半角积分区域的功率Ps？可以使用P=η E^2 *A直接获得吗，这样的结果与使用波印廷矢量计算是否相同？
4. 如果我不使用例程中的BSDF分析组（分析组中使用gratingpolar获取特定级次方向的E^2），而直接在上方用监视器获取近场分布，再远场投射到1m球面，这样就能获得更多的方向和数据点，这样获得的E2是否是和gratingpolar获得的相同？

谢谢