六：大气效应校正和大气参数反演.pptVIP

大气的遥感探测 * * 第六章 大气效应及大气纠正和大气参数反演 大气效应纠正 ATMOSHPHERIC EFFECT CORRECTION 大气效应纠正 ATMOSHPHERIC EFFECT CORRECTION 利用遥感数据进行应用时，首要的一步是要从卫星对地观测所得上行信号中区分大气和地表各自的贡献，然后再从各自的贡献中分析反演大气和地表的物理参数。 从物理实质上看，这是地-气系统辐射传输问题，对地表遥感而言，即为大气效应校正问题，而对大气遥感而言，则是地表背景作用的扣除问题，确切的说，这是同一个问题的两个方面。 对同一波长而言，卫星对地观测在同一时刻只有一个观测值，而至少有两个或者两个以上的未知量（即大气光学厚度和地表反照率），因此问题的解事不确定的，必须要增加新信息，以解决反演求解的不确定性。 大气效应纠正 ATMOSHPHERIC EFFECT CORRECTION 回顾平面平行大气的辐射传输方程： 给定边界条件，可以求解辐射传输方程。 实际上，遥感图像大气校正的目的是求取地表反射率，从某种意义上来说，是要求解边界条件。从辐射传输方程来看，这是一个十分复杂的过程。 1．零反射地表假设 零反射地表是对地表—大气耦合系统辐射传输问题的简化，它假设地表反射率为零，传感器接收到的辐射全部来自大气。 通过零反射地表假设，我们可以求解大气程辐射，从而将遥感图像中来自地表的辐射和来自大气的辐射分离，使遥感大气订正模型得以大大简化。 （1）对于近红外与可见光波段，大气自身辐射可以忽略不计，大气路径辐射顶主要来源于大气对太阳辐射的多次散射。 对于多次散射的计算可以通过求解辐射传输方程（如离散坐标法），或者直接用蒙特卡洛方法模拟获得。在实际应用中，我们又往往通过已有的大气辐射传输模式软件来计算，如6S软件，Lowtran软件，Modtran软件等。 （2）对于热红外波段，多次散射一般可以忽略不计，但大气和地表自身发射必须考虑。 （3）对于中红外波段，则既需要考虑地表与大气自身的发射，同时又要考虑大气的多次散射作用，因此更加复杂，我们不展开讨论。 2．朗伯体地表大气订正方法 设N+代表离开大气上界的总光子数，则： 代表与地表碰撞i次以后穿越大气的光子数，若R(0)为地表发射率， “-”代表光子向下运动，n代表处在大气低层的光子数，令： g描述了与地表碰撞了i次以后，所有向上运动的光子数中，有多少个光子最终能穿出大气，所以g值只与大气状态有关，与地表状况无关，令： h描述了与地表碰撞i次以后，向上运动的光子数中究竟有多少光子会被大气弹回地面，所以h值只决定于大气状况 2．朗伯体地表大气订正方法 所以： 因为h与R均小于1，所以： 如果假定第一次碰撞的反射率为1，则： 2．朗伯体地表大气订正方法 仅讨论可见光/近红外波段 为了问题的简化，在地表朗伯体、大气水平均一假设条件下，我们可以得到： 其中 ， ； 分别为观测天顶角与太阳天顶角； 为传感器接受到的辐射亮度， 为观测方向的路径辐射项； 为地表反射率；S为大气下界的半球反射率； 为大气层顶与太阳光垂直方向的通量密度。 利用入射太阳辐射项 归一化上式可得： 式中 是大气顶部反射率， 是大气的路径辐射项等效反射率，将T可以分为直射和散射两部分： 是总体光学厚度， 是散射透射。 如何进行遥感图像的大气校正呢？即如何求得地表反射率？ （1）已知大气状况的校正方法 如果已知大气的垂直廓线（温、湿、压），大气水汽含量，大气光学厚度（气溶胶），以及气溶胶模式，我们可以通过大气辐射传输模式模拟，计算三个不同地表反射率条件下的大气层顶辐射亮度，进而求解得到 ，S和F值，而达到大气效应纠正的目的。对于可见光/近红外波段而言，大气状况最主要的影响因素是大气气溶胶的变化，即大气光学厚度的变化，因此关键是如何估算遥感图像的气溶胶分布。 （2）参考地物法(Invariant—Object Methods) 对于TM等高分辨率的图像，通常可以假设整幅图像的大气状况相同，图像的观测天顶角也可以近似看作相同（垂直观测），如果可以在一幅图像中找到苦干个（3个以上）地表反射率固定不变的地物，则可以通过求解方程组得到大气参数，进而达到整幅图像大气校正的目的。 ? 其中 ； 理论上，n≥3即可求解，但为避免方程的相关性，一般要寻找地表反射率分别为高、中、差的地物作为参考物。 该方法对于时间序列的多幅