返回笔记列表
01 · NOTES / 遥感

气溶胶反演

2026 年 6 月 12 日

气溶胶反演

目的

反演气溶胶厚度(Aerosol Optical Depth,AOD),表征大气污染时空格局特征

意义

气溶胶(aerosol)是固体或液体微粒均匀地悬浮于气体介质中形成的分散体系。由于其中的微粒比气态分子大而比粗尘颗粒小,因而它们不像气态分子那样服从气体分子运动规律,但也不会像粗尘颗粒那样受地心引力作用而沉降,使整个体系具有胶体性质,故称为气溶胶。气溶胶是大气污染和气候变化的重要因素。此外,当气溶胶的浓度达到足够高时,将对人类健康造成威胁,尤其是对哮喘病人和其他有呼吸性疾病的人群。空气中的气溶胶还能传播细菌和病毒,这可能会导致一些地区疾病的流行和暴发。

气溶胶光学厚度(AOD)是气溶胶的光学属性之一,表征光在单位截面面积垂直气柱上的透过率,也称为大气浑浊度。AOD是一个无量纲的正值,该值越大,光的大气透过率越低。由于人类生活和生产活动产生的大气污染物可以生成人为气溶胶,因此AOD在一定程度上可反映大气污染状况

对大气污染时空格局特征的大面积监测。MODIS-L1B是MODIS系列2级数据产品中的一种,产品编号为MOD02/MYD02。该数据产品已经过仪器定标,尚未进行大气校正;包含地理坐标,但是地理坐标与科学数据没有连接,当直接显示时,数据边缘会存在bow-tie效应(又名蝴蝶结效应)。因此,当对该数据产品进行相应处理后,即可完成AOD遥感反演,实现对大气污染时空格局特征的动态监测。

  1. 算对“全球变暖”:帮助科学家判断气溶胶(比如雾霾和黑碳)到底是在给地球降温还是升温,减少气候预测的误差。

  2. 看清“空气污染”:把看不见的PM2.5变成一张张污染地图,告诉我们哪里空气差、污染源在哪,以及出门要不要戴口罩。

  3. 校正“卫星照片”:就像PS去掉照片里的雾气,让卫星拍下的地面真实情况(农田、森林、城市)变得清晰准确。

数据

MODIS-L1B(MOD02)

步骤

  1. 数据预处理(辐射定标、几何校正)

  2. 选择算法:

    1. 暗目标(像元)法:适用于浓密植被区,利用浓密植被区蓝红波段反射率低,那么地表反射即为气溶胶浓度

    2. 深蓝算法:适用于城市、干旱区,利用深蓝波段,地表反射代表气溶胶浓度

  3. 云检测

  4. 气溶胶反演

  5. 后处理操作

打开数据

image.png

发射率数据、反射率数据处理

  • 发射率数据几何校正,输出GCP,坐标选择WGS84

image.png

输出GCP

image.png

image.png

结果如图

image.png

  • 反射率数据几何校正(步骤相同)

image.png

image.png

image.png

结果如图

image.png

反射率、发射率合成(波段组合)与裁剪

  • 波段合成(两张单波段图像合成一张双波段图像,注意反射率在发射率上面,选择研究区矢量裁剪)

image.png

image.png

结果如图

image.png

角度数据处理

  • 角度数据集几何校正(HDF4格式打开头文件,依次加入卫星天顶角、卫星方位角、太阳天顶角和太阳方位角)

MODIS数据的角度数据集的行列号和科学数据集的行列号不同,所以要使用校正发射率的GCP控制点校正角度数据集,在几何校正前,要对角度数据集进行波段合成和重采样。角度数据集包括卫星天顶角、卫星方位角、太阳天顶角和太阳方位角等数据集。

image.png

image.png

image.png

  • 对角度数据进行重采样

先查看角度数据集的合成结果与发射率数据集的行列数。view metadata,然后重采样

image.png

image.png

image.png

image.png

  • 最后进行角度数据集几何校正

利用反射率或者发射率几何校正时的GCP,校正角度数据集,注意更改几何校正方法(Triangulation)和重采样方法

image.png

image.png

image.png

  • 角度数据波段运输

HDF文件中的角度数据值为实际值的100倍,所以在气溶胶反演之前要把角度好吃数据乘0.01,利用波段计算器,点Map Variable to Input File可以选择一个图像所有波段,同时选择研究区裁剪

image.png

image.png

云检测

对反射率和发射率数据合成文件进行去云处理。

image.png

image.png

气溶胶反演

气溶胶反演原理是利用暗像元法反演地表真实反射率,并使用查找表进行气溶胶浓度查找,查找表是3-9月MODIS数据通用的

image.png

image.png

  • 结果如图

image.png

  • 可视化

image.png