大气遥感中的方位角和散射角

学习大气辐射传输方程，会遇到各种方向问题，直接关系到表达式里的正负号。
对于大气遥感领域，角度的正负会影响到散射角的计算，进而会影响粒子散射相函数的正确性。
本博文将梳理气溶胶遥感中所涉及到的相对方位角、散射角的问题，以及在实际6s查找表反演中的运用情况。
注：以下代码仅示意

参数定义

太阳天顶角：SZA
太阳方位角：phi0
观测天顶角：VZA
观测方位角：phi
散射角：SC

根据球面几何学（“大气辐射导论K.N.LIOU”,附录C），公式1如下。

raa = abs(phi - phi0)
if raa > 180 then raa = 360 - raa
A = acos(cos(VZA*dtr)*cos(SZA*dtr) + sin(VZA*dtr)*sin(SZA*dtr)*cos(raa*dtr))*rtd

上式的定义与6S是一致的，这里附上6S的示意图。

这里的A为卫星和太阳之间的夹角。

大气卫星遥感

大气遥感领域的散射角通常定义为入射光线和散射光线的夹角，也就是太阳方向的延长线与卫星的夹角。这里附上NASA官网的示意图link。
这里$\theta$就是散射角。公式2如下。

raa = abs(phi - phi0)
if raa > 180 then raa = 360 - raa
raa = 180 - raa //对raa取余角
SC = acos(cos((180-vza)*dtr)*cos(sza*dtr) + sin(vza*dtr)*sin(sza*dtr)*cos(raa*dtr))*rtd  
   = acos(-cos(vza*dtr)*cos(sza*dtr) + sin(vza*dtr)*sin(sza*dtr)*cos(raa*dtr))*rtd

这里的入射光线不同于公式1，是太阳光的延长线，因此与天顶方向的夹角应为180 - vza

举例说明

来源于（Levy et al., 2007）

这个表完美的契合了公式2。
如果大家选用6s模型建立查找表进行气溶胶的反演的话，你会发现用table的值作为输入，结果并不同。

说明6s的相对方位角还是公式1中的相对方位角的定义。

所以这样设置就对啦。
告诉你们一个验证策略：假设卫星和太阳位置几乎相同，此时的散射角接近180°。

你们看明白了吗？？？
如果对你有帮助，可以👍👍👍支持，欢迎交流讨论。