专利名称：一种遥感卫星多光谱数据的大气订正方法
技术领域：
本发明涉及一种遥感卫星的大气订正方法，特别涉及一种与环境减灾卫星传感器波段设置相类似卫星的多光谱数据大气订正方法。
背景技术：
环境减灾卫星自发射以来，在汶川地震连续检测、印度洋海啸灾情监测、澳大利亚火灾监测、日本地震和海啸灾情监测等灾害监测中发挥了独特的作用，并广泛应用于农业、 林业、水利、土地利用、地况调查、城市规划等众多领域，取得了显著的效果。为了实现对不同地点、不同时间和不同的传感器获取的遥感数据进行比较和应用，及时为国家环境保护和灾害预报、抗灾、救灾提供准确、科学的决策依据，这必然要求遥感数据定量化。对于遥感定量化，就要详细考虑大气对于成像波段传输过程的定量影响，对其影响做出准确的订正，实现地物目标辐射特性真实再现，达到提高遥感信息质量及定量化应用水平目的。目前国外大气订正方面的成果很多，像LANDSAT、SPOT、MODIS、ASTER等其他国外遥感卫星都针对各自卫星的特点开展了相应的大气订正；国内学者也在CBERS-Ol卫星、 CBERS-02、02B卫星等国产遥感卫星数据的大气订正研究产生了一些相关算法，但由于大气状况瞬息万变，当时、当地的大气对遥感数据影响程度和重要性不同，再加上各卫星的传感器波段设置不同，特性不同，研究目的和要求也不同，即使相同地区，大气状况也是复杂多变的，每景图像成像时的大气状况也存在明显差异，很难有一种算法是普遍适用的。
发明内容
本发明的技术解决问题是克服现有技术的不足，提供一种普适性强的遥感卫星多光谱数据大气订正方法。本发明的技术解决方案是一种遥感卫星多光谱数据的大气订正方法，步骤如下(1)读取遥感卫星多光谱数据，对该多光谱数据进行预处理，将多光谱数据转换为表观辐亮度数据和表观反射率数据；(2)依据表观反射率数据计算决策树法中的决策因子比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI，通过决策树法实现DDV暗目标自动提取；(3)根据步骤(1)中的表观反射率数据获取上述提取的DDV暗目标区域的DDV表观反射率数据；并根据该DDV表观反射率数据结合DDV地表反射率计算I^ath反射率；(4)根据I^ath反射率和Rayleigh反射率比值与气溶胶光学厚度的关系，确定多光谱各波段的气溶胶光学厚度；并依据多光谱各波段的气溶胶光学厚度和feyleigh散射光学厚度，确定出550nm处总的光学厚度τ ；(5)依据上述550nm处总的光学厚度，利用查找表获取反演参数P path( λ )、Τ( λ )、 S(A),通过辐射传输方程求解地表反射率，完成大气订正。
所述步骤(3)中的DDV地表反射率的计算公式如下P 丽(入)=SvegPveg(A) + (I-Sveg) P soil ( λ )其中ρ ■ ( λ )对应多光谱波段λ的DDV地表反射率；Pveg(A), Psoil(A)是实测多光谱波段λ的植被反射率和土壤反射率；
\SAVI SAVI >0Sreg是植被调整因子Aeg =I SAVI<0所述步骤(5)中查找表的索引τ分别与Ρρ_(λ)、Τ(λ)、S(X)之间的关系满足下式Phth(X) = C0U)+C1U) τ+C2(X) τ2S(A) = d0(A)+d1(A) T+d2(A) τ2T ( λ ) = e0(A )+θι(λ ) τ +e2 ( λ ) τ 2其中，C。( λ )、C1 ( λ )、C2 ( λ )、d。( λ )、Cl1 ( λ )、d2 ( λ )、％ ( λ )、Θι ( λ )、％ ( λ )对应多光谱波段λ在确定的卫星成像几何参数太阳天顶角θ s、卫星观测天顶角θ v、太阳与卫星的方位角炉下的路径反射率0_(入)、大气透过率1~(入)、半球反射率3(入)与总光学厚度τ的系数。本发明与现有技术相比有益效果为(1)本发明基于浓密植被暗目标区(DDV-density dark vegetation)反射率与土壤可调植被指数及波谱库实测植被和土壤地表反射率的相对稳定关系，利用辐射传输方程确定出路径反射率，然后利用为环境减灾卫星多光谱传感器建立的路径反射率和瑞利散射反射率比值与气溶胶光学厚度关系的查找表确立气溶胶光学厚，为此实现与环境减灾卫星传感器波段设置相似类卫星的大气订正算法，方法完整、合理可行、精度较高。(2)本发明针对遥感卫星多光谱波段设置特点，依据不同地物特性，选取大量遥感数据，分析多种植被指数和水体指数，甄选出比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI作为决策树法的决策因子，实现了暗目标自动提取，使此类卫星多光谱数据的大气订正更具有普适性。(3)本发明利用M0DTRAN辐射传输模型，分析遥感卫星在不同几何状态下的各种情况，针对此类卫星多光谱数据波段设置特点，为大气订正各反演参数与主变量间建立二次线性查找表，有效地缩减了大气订正产品生产时间，使得业务化大气订正得以实现。(4)本发明铜鼓大量实测地物波谱数据，结合遥感卫星地表反演量，针对多光谱数据的DDV暗目标区，分析实测地物波谱与多种植被指数的线性关系，确立了用土壤可调植被指数SAVI作为植被实测植被和土壤地表反射率的调整因子，使DDV暗目标区的地表反射率计算更加合理，确保了方法的合理可行性。
图1为本发明总体流程图；图2为本发明具体流程图；图3为本发明大气订正结果对比曲线。
具体实施方式
如图1、2所示，本发明ー种遥感卫星多光谱数据的大气订正方法，步骤如下(1)读取遥感卫星多光谱数据，对该多光谱数据进行预处理，将多光谱数据转换为 表观辐亮度数据和表观反射率数据；(1. 1)表观辐亮度读取多光谱CCD数据，依据数据附帯文件中的绝对定标系数，按照公式①计算多 光谱传感器入瞳处辐亮度即表观辐亮度
权利要求
1.一种遥感卫星多光谱数据的大气订正方法，其特征在于步骤如下(1)读取遥感卫星多光谱数据，对该多光谱数据进行预处理，将多光谱数据转换为表观辐亮度数据和表观反射率数据；(2)依据表观反射率数据计算决策树法中的决策因子比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI，通过决策树法实现DDV暗目标自动提取；(3)根据步骤(1)中的表观反射率数据获取上述提取的DDV暗目标区域的DDV表观反射率数据；并根据该DDV表观反射率数据结合DDV地表反射率计算I^ath反射率；(4)根据I^ath反射率和Rayleigh反射率比值与气溶胶光学厚度的关系，确定多光谱各波段的气溶胶光学厚度；并依据多光谱各波段的气溶胶光学厚度和feyleigh散射光学厚度，确定出550nm处总的光学厚度τ ；(5)依据上述550nm处总的光学厚度，利用查找表获取反演参数Ppath( λ )、T ( λ )、 S(A),通过辐射传输方程求解地表反射率，完成大气订正。
2.根据权利要求1所述的一种遥感卫星多光谱数据的大气订正方法，其特征在于所述步骤(3)中的DDV地表反射率的计算公式如下P DDV( λ ) = SvegPveg(A)+ (I-Sveg) Psoil(X)其中P DDV( λ )对应多光谱波段λ的DDV地表反射率；Pveg(A), Psoil(A)是实测多光谱波段λ的植被反射率和土壤反射率；\SAVI SAVI >0Sveg是植被调整因子，&g = 10 SAVKOο
3.根据权利要求1所述的一种遥感卫星多光谱数据的大气订正方法，其特征在于所述步骤(5)中查找表的索引τ分别与ppathU)、T(X)、S(X)之间的关系满足下式Ppath(A) = C0(A)+C1(A) τ+C2(X) τ2S(A) = d0 ( λ ) +Cl1 ( λ ) τ +d2 ( λ ) τ 2Τ(λ) = e(1U)+eiU) τ+θ2(λ) τ2其中，C。( λ )、Cl ( λ )、c2 ( λ )、d。( λ )、Cl1 ( λ )、d2 ( λ )、％ ( λ )、ei ( λ )、％ ( λ )对应多光谱波段λ在确定的卫星成像几何参数太阳天顶角θ s、卫星观测天顶角θ v、太阳与卫星的方位角P下的大气路径反射率0_11(入)、大气透过率仪入)、半球反射率3(入)与总光学厚度τ的系数。
全文摘要
一种遥感卫星多光谱数据的大气订正方法，(1)读取数据，将数据转换为表观辐亮度数据和表观反射率数据；(2)依据表观反射率数据计算决策树法中的决策因子比值植被指数、土壤可调植被指数和归一化水体指数，通过决策树法实现暗目标自动提取；(3)根据表观反射率数据获取上述暗目标区域的表观反射率数据，并计算大气路径反射率；(4)根据大气路径反射率和瑞利反射率比值与气溶胶光学厚度的关系，确定多光谱各波段的气溶胶光学厚度；并依据多光谱各波段的气溶胶光学厚度和瑞利散射光学厚度，确定出550nm处总的光学厚度；(5)依据上述550nm处总的光学厚度，利用查找表获取反演参数，通过辐射传输方程求解地表反射率，完成大气订正。
文档编号G01S7/497GK102288956SQ20111011964
公开日2011年12月21日 申请日期2011年5月10日 优先权日2011年5月10日
发明者傅俏燕, 李杏朝, 潘志强, 王爱春, 闵祥军, 韩启金 申请人:中国资源卫星应用中心