课件-遥感物理

第一章基本概念遥感物理第二节辐射传输(radiancetransfer)√§1.2.1传输方程

§1.2.2源函数中散射的表达Maxwell方程组与辐射传输方程麦克斯韦方程组描述了电磁场的基本规律。一般而言，波长较长的电磁波波动性较为突出。短波部分干涉与衍射等波动现象则不明显，而更多地表现为粒子性。在光学和热红外遥感领域，为方便和直观起见，常用辐射传输方程描述电磁波与介质的相互作用。部分辐射传输方程加入了反映波动性的修正因子。麦克斯韦方程组与辐射传输方程是不矛盾的，可以相互转换，不存在难易和优劣之分，只不过形式和求解方法有所区别，在不同的领域，有各自的优势。1/14消光截面在光散射和辐射传输领域中，通常用“截面”这一术语，它与几何面积类似，用来表示粒子由初始光束中所移除的能量大小。当对粒子而言时，截面的单位是面积（厘米2），因此，以面积计的消光截面等于散射截面与吸收截面之和。但当对单位质量而言时，截面的单位是每单位质量的面积（厘米2·克-1），这时，在传输研究中用术语质量消光截面，因而，质量消光截面等于质量散射截面与质量吸收截面之和。此外，当消光截面乘以粒子数密度（厘米-3）或当质量消光截面乘以密度（克·厘米-3）时，该量称为“消光系数”，它具有长度倒数（厘米-1）的单位。2/14传输方程在介质中传输的一束辐射，将因它与物质的相互作用而减弱。如果辐射强度Iλ，在它传播方向上通过ds厚度后变为Iλ+dIλ，则有：

dIλ=-kλρIλds式中ρ是物质密度，kλ表示对辐射波长λ的质量消光截面。辐射强度的减弱是由物质中的吸收以及物质对辐射的散射所引起。设σe为粒子消光截面，N为单位体积的总粒子数，上式如何表达？消光系数=？3/144/14I(0)I(s1)I+dII0dsS1另一方面，辐射强度也可以由于相同波长上物质的发射以及多次散射而增强，多次散射使所有其它方向的一部分辐射进入所研究的辐射方向。我们如下定义源函数系数，使由于发射和多次散射造成的强度增大为：

dIλ=jλρds式中源函数系数jλ具有和质量消光截面类似的物理意义。联合上述两个方程得到辐射强度总的变化为：

dIλ=-kλρIλds+jλρdsjλ的单位与kλ的单位不同：前者带有强度概念。5/14进一步为方便起见，定义源函数Jλ如下：

Jλ≡jλ/kλ这样一来，源函数则具有辐射强度的单位。因此有：

dIλ=-kλρIλds+kλJλρds即：6/14这就是不加任何座标系的普遍传输方程，它是讨论任何辐射传输过程的基础。求解辐射传输方程时，最难解决的是Jλ。比尔-布格-朗伯(Beer-Bouguer-Lambert)定律当忽略多次散射和发射的增量贡献时，辐射传输方程可以简化为：7/14如果在s=0处的入射强度为Iλ(0)，则在s1处，其射出强度可以通过对上式的积分获得：假定介质消光截面均一不变，即kλ不依赖于距离s，并定义路径长度：8/14这就是著名的比尔定律，或称布格定律，也可称朗伯定律。它叙述了忽略多次散射和发射影响时，通过均匀介质传播的辐射强度按简单的指数函数减弱，该指数函数的自变量是质量吸收截面和路径长度的乘积。由于该定律不涉及方向关系，所以它不仅适用于强度量，而且也适用于通量密度。介质完全均一（ρ也不依赖s），出射强度？则此时出射强度为：光学厚度(opticalthickness,opticaldepth)定义点s1和s2之间的介质的光学厚度为：9/14并有：

dτλ(s)=-kλρds因此传输方程可以写为：在实际应用中，τ的定义使τ永远是正数。而且I与τ的关系一般为exp(-τ)。平面平行(planeparallel)介质在遥感定量分析过程中，为简化起见，我们通常假设电磁波穿过的介质（如大气与植被冠层）是平面平行的，或称水平均一(horizontallyuniform)的。即介质可以分成若干或无穷多相互平行的层，各层内部（对辐射影响）的性质一样，各层之间的性质不同。10/14θθ为辐射方向与分层方向法线的夹角。z上述传输方程用z、θ替换s后，具体表达式？对于平面平行介质，辐射传输方程可以写为：11/14或其中μ=cosθ，τ是光学厚度(此时已是垂直计量)

。注意μ

，多数情况下，它会代替θ在辐射传输中出现对于平面平行大气，τ的定义为由大气上界向下测量的垂直光学厚度（省略下标λ）：12/14对于水平均一植被，τ的定义为由冠层表面向下测量到z处的垂直光学厚度：其中uL为叶面积密度。大气植被冠层0zz在植被中，dτ与dz关系如何？以平面平行大气为例，比尔定律具体表达式？对于平面平行大气，且忽略大气中的多次散射和发射，则传输方程为：13/14上式的解为：定义τ0=τ(0)为大气整层光学厚度，注意到τ(∞)=0，因此有：请注意指数形式在辐射传输中的作用。总结两个概念：光学厚度、平面平行介质14/14一组不同表达形式的传输方程：传输方程的简单解（比尔定律）：e的指数形式第一章基本概念遥感物理第二节辐射传输(radiancetransfer)

§1.2.1传输方程√

§1.2.2源函数中散射的表达散射电磁波通过介质时，会发生散射，即电磁波有可能改变方向。因此使某一方向的电磁波强度发生变化，可能减弱，也可能增强。1/13散射相函数（scatteringphasefunction）为描述电磁波被介质散射后在各个方向上的强度分布比例，引入散射相函数P(Ω,Ω’)，并且P(Ω,Ω’)/4π是归一化的，即：根据互易原理：因此同样有：2/13物理意义：P(Ω,Ω’)/4π为方向Ω’的电磁波被散射到方向Ω的比例。作业：对于在4π空间内各向均一的散射（散射辐射强度不随散射方向变化），散射相函数的表达式是什么？对于散射光只在入射方向Ω’存在，其它方向均为0的情况下，散射相函数的表达式是什么？3/13通常散射相函数P(Ω,Ω’)只与方向Ω’和方向Ω之间的夹角Θ有关，可以写为P(cosΘ)。散射角Θ定义为入射光束和散射光束之间的夹角。散射角的余弦可以表示为：请注意P与两个方向的天顶角，以及相对方位角有关。散射角余弦推导公式参见重要参考资料中的学生来信。4/13单次散射反射率（singlescatteringalbedo）实际上辐射被介质散射的同时，也被介质吸收，即消光过程既包括散射，也包括吸收。单次散射反射率ω定义为辐射发生每一次消光（或简称散射）过程中，遭受散射的百分比。入射为1，散射后各个方向的总和（积分）即为ω5/13由于位置τ处散射产生的源函数为来自所有方向、并经散射，到方向Ω的辐射总和。即上式对方向Ω’在4π空间的积分：6/13我们假设位于τ处、传播方向为Ω’的辐射强度为I(τ,Ω’)，则它散射到方向Ω的辐射强度为：源函数中散射的表达Ω’Ωτ因此，含有散射源函数的辐射传输方程展开为：通常情况下，这个方程没有解析解，只能靠数值解法或简化求解。在很多应用中，我们经常把上式等号右面第二项拆分成单次散射和多次散射两项。回忆上一小节中提到的平面平行介质（大气）中的传输方程为：7/13当电磁波由方向Ω0前进时，它被介质散射到方向Ω的散射过程可以分解为单（一）次散射和多次散射过程。多次散射是为了区别单次散射而定义的，凡是辐射被介质散射超过1次，均称为多次散射。区分单次散射和多次散射是为了方便于求解辐射传输方程。8/13Ω0Ω单次散射多次散射τ下面我们把前式9/13在τ处发生单次散射后，散射到方向Ω的辐射强度即为：中的I(τ,Ω’)包含的单次散射单独抽出来，推导其表达式。对于单次散射，我们假设入射辐射强度的初始值为I0，传播方向为Ω’，则它到达τ处的辐射强度为：对上式中入射方向在4π空间积分，并考虑只有一个入射方向：比如由于太阳光是平行光入射，所以只在Ω0一个方向上，即I0=δ(Ω’,Ω0)F0，其中F0为太阳的通量密度。则单次散射产生的源函数，即散射到方向Ω的辐射强度为

：10/13太阳为什么是单次散射的源?11/13注意此处的多次散射为二次以上散射源函数中的散射可以表达为单次散射与多次散射之和，其中多次散射项表达式与前述相同，即：J(τ,Ω)=上式结果肯定是强度单位12/13其中B(T)为普朗克函数，是物体亮温为T时发射的出射辐