CN102288956B - 一种遥感卫星多光谱数据的大气订正方法 - Google Patents
一种遥感卫星多光谱数据的大气订正方法 Download PDFInfo
- Publication number
- CN102288956B CN102288956B CN2011101196453A CN201110119645A CN102288956B CN 102288956 B CN102288956 B CN 102288956B CN 2011101196453 A CN2011101196453 A CN 2011101196453A CN 201110119645 A CN201110119645 A CN 201110119645A CN 102288956 B CN102288956 B CN 102288956B
- Authority
- CN
- China
- Prior art keywords
- reflectivity
- data
- ddv
- path
- optical depth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
Abstract
一种遥感卫星多光谱数据的大气订正方法,(1)读取数据,将数据转换为表观辐亮度数据和表观反射率数据;(2)依据表观反射率数据计算决策树法中的决策因子:比值植被指数、土壤可调植被指数和归一化水体指数,通过决策树法实现暗目标自动提取;(3)根据表观反射率数据获取上述暗目标区域的表观反射率数据,并计算大气路径反射率;(4)根据大气路径反射率和瑞利反射率比值与气溶胶光学厚度的关系,确定多光谱各波段的气溶胶光学厚度;并依据多光谱各波段的气溶胶光学厚度和瑞利散射光学厚度,确定出550nm处总的光学厚度;(5)依据上述550nm处总的光学厚度,利用查找表获取反演参数,通过辐射传输方程求解地表反射率,完成大气订正。
Description
技术领域
本发明涉及一种遥感卫星的大气订正方法,特别涉及一种与环境减灾卫星传感器波段设置相类似卫星的多光谱数据大气订正方法。
背景技术
环境减灾卫星自发射以来,在汶川地震连续检测、印度洋海啸灾情监测、澳大利亚火灾监测、日本地震和海啸灾情监测等灾害监测中发挥了独特的作用,并广泛应用于农业、林业、水利、土地利用、地况调查、城市规划等众多领域,取得了显著的效果。为了实现对不同地点、不同时间和不同的传感器获取的遥感数据进行比较和应用,及时为国家环境保护和灾害预报、抗灾、救灾提供准确、科学的决策依据,这必然要求遥感数据定量化。
对于遥感定量化,就要详细考虑大气对于成像波段传输过程的定量影响,对其影响做出准确的订正,实现地物目标辐射特性真实再现,达到提高遥感信息质量及定量化应用水平目的。
目前国外大气订正方面的成果很多,像LANDSAT、SPOT、MODIS、ASTER等其他国外遥感卫星都针对各自卫星的特点开展了相应的大气订正;国内学者也在CBERS-01卫星、CBERS-02、02B卫星等国产遥感卫星数据的大气订正研究产生了一些相关算法,但由于大气状况瞬息万变,当时、当地的大气对遥感数据影响程度和重要性不同,再加上各卫星的传感器波段设置不同,特性不同,研究目的和要求也不同,即使相同地区,大气状况也是复杂多变的,每景图像成像时的大气状况也存在明显差异,很难有一种算法是普遍适用的。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种普适性强的遥感卫星多光谱数据大气订正方法。
本发明的技术解决方案是:一种遥感卫星多光谱数据的大气订正方法,步骤如下:
(1)读取遥感卫星多光谱数据,对该多光谱数据进行预处理,将多光谱数据转换为表观辐亮度数据和表观反射率数据;
(2)依据表观反射率数据计算决策树法中的决策因子:比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI,通过决策树法实现DDV暗目标自动提取;
(3)根据步骤(1)中的表观反射率数据获取上述提取的DDV暗目标区域的DDV表观反射率数据;并根据该DDV表观反射率数据结合DDV地表反射率计算Path反射率;
(4)根据Path反射率和Rayleigh反射率比值与气溶胶光学厚度的关系,确定多光谱各波段的气溶胶光学厚度;并依据多光谱各波段的气溶胶光学厚度和Rayleigh散射光学厚度,确定出550nm处总的光学厚度τ;
(5)依据上述550nm处总的光学厚度,利用查找表获取反演参数ρpath(λ)、T(λ)、S(λ),通过辐射传输方程求解地表反射率,完成大气订正。
所述步骤(3)中的DDV地表反射率的计算公式如下:
ρDDV(λ)=Svegρveg(λ)+(1-Sveg)ρsoil(λ)
其中:ρDDV(λ)对应多光谱波段λ的DDV地表反射率;
ρveg(λ)、ρsoil(λ)是实测多光谱波段λ的植被反射率和土壤反射率;
所述步骤(5)中查找表的索引τ分别与ρpath(λ)、T(λ)、S(λ)之间的关系满足下式:
ρPath(λ)=c0(λ)+c1(λ)τ+c2(λ)τ2
S(λ)=d0(λ)+d1(λ)τ+d2(λ)τ2
T(λ)=e0(λ)+e1(λ)τ+e2(λ)τ2
其中,c0(λ)、c1(λ)、c2(λ)、d0(λ)、d1(λ)、d2(λ)、e0(λ)、e1(λ)、e2(λ)对应多光谱波段λ在确定的卫星成像几何参数太阳天顶角θS、卫星观测天顶角θV、太阳与卫星的方位角下的路径反射率ρpath(λ)、大气透过率T(λ)、半球反射率S(λ)与总光学厚度τ的系数。
本发明与现有技术相比有益效果为:
(1)本发明基于浓密植被暗目标区(DDV-density dark vegetation)反射率与土壤可调植被指数及波谱库实测植被和土壤地表反射率的相对稳定关系,利用辐射传输方程确定出路径反射率,然后利用为环境减灾卫星多光谱传感器建立的路径反射率和瑞利散射反射率比值与气溶胶光学厚度关系的查找表确立气溶胶光学厚,为此实现与环境减灾卫星传感器波段设置相似类卫星的大气订正算法,方法完整、合理可行、精度较高。
(2)本发明针对遥感卫星多光谱波段设置特点,依据不同地物特性,选取大量遥感数据,分析多种植被指数和水体指数,甄选出比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI作为决策树法的决策因子,实现了暗目标自动提取,使此类卫星多光谱数据的大气订正更具有普适性。
(3)本发明利用MODTRAN辐射传输模型,分析遥感卫星在不同几何状态下的各种情况,针对此类卫星多光谱数据波段设置特点,为大气订正各反演参数与主变量间建立二次线性查找表,有效地缩减了大气订正产品生产时间,使得业务化大气订正得以实现。
(4)本发明铜鼓大量实测地物波谱数据,结合遥感卫星地表反演量,针对多光谱数据的DDV暗目标区,分析实测地物波谱与多种植被指数的线性关系,确立了用土壤可调植被指数SAVI作为植被实测植被和土壤地表反射率的调整因子,使DDV暗目标区的地表反射率计算更加合理,确保了方法的合理可行性。
附图说明
图1为本发明总体流程图;
图2为本发明具体流程图;
图3为本发明大气订正结果对比曲线。
具体实施方式
如图1、2所示,本发明一种遥感卫星多光谱数据的大气订正方法,步骤如下:
(1)读取遥感卫星多光谱数据,对该多光谱数据进行预处理,将多光谱数据转换为表观辐亮度数据和表观反射率数据;
(1.1)表观辐亮度
读取多光谱CCD数据,依据数据附带文件中的绝对定标系数,按照公式①计算多光谱传感器入瞳处辐亮度即表观辐亮度:
其中:λ是多光谱的波段;
L(λ)对应多光谱波段λ的入瞳处辐亮度;
DN(λ)对应多光谱波段λ的数字DN值;
A(λ)、L0(λ)对应多光谱波段λ的定标系数。
(1.2)表观反射率数据
将上述计算得到的传感器入瞳处辐射亮度,按公式②计算得到传感器入瞳处的表观反射率数据:
其中:ρTOA(λ)对应多光谱波段λ的表观反射率;
d2是日-地距离校正因子可通过公式③计算
J是儒略历天数;
ES(λ)对应多光谱波段λ的大气外太阳光谱辐照度;
θS是卫星成像时的太阳天顶角,可从数据附带文件中成像几何参数中获取;
上述数据附带文件是指遥感卫星成像时的已知状态参数,根据卫星下传数据得到,对本领域技术人员来讲为已知量。
(2)依据表观反射率数据计算决策树法中的决策因子:比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI,通过决策树法实现DDV暗目标自动提取;
其中:ρTOA(λgreen)是多光谱数据绿光波段的表观反射率;
ρTOA(λred)是多光谱数据红光波段的表观反射率;
ρTOA(λnir)是多光谱数据近红外波段的表观反射率;
C是土壤调整因子,C∈[0,1]。
(3)针对此类卫星多光谱传感器的响应特性,确立大气组成中吸收、散射等衰减对多光谱数据的主要影响因子(Rayleigh散射衰减和气溶胶衰减),依此根据步骤(1)中的表观反射率数据获取上述提取的DDV暗目标区域的DDV表观反射率数据ρPath-DDV(λ);并根据该DDV表观反射率数据结合DDV地表反射率计算Path反射率;
(3.1)瑞利(Rayleigh)散射反射率
利用公式⑦,计算DDV暗目标选取区相对应DEM高程数据处多光谱波段λ的Rayleigh散射反射率:
ρRay-DDV(λ)=a0(λ)+a1(λ)H+a2(λ)H2⑦
其中:H是DEM高程数据获得海拔高度;
a0(λ)、a1(λ)、a2(λ)对应多光谱波段λ的Rayleigh散射反射率与海拔高度H关系的系数;
(3.2)DDV地表反射率
利用大量实测地物波谱库,建立DDV地表反射率与土壤可调植被指数SAVI的拟合关系,由公式⑧确定多光谱波段λ的DDV地表反射率:
ρDDV(λ)=Svegρveg(λ)+(1-Sveg)ρsoil(λ)⑧
其中:ρDDV(λ)对应多光谱波段λ的DDV地表反射率;
ρveg(λ)、ρsoil(λ)是实测地物波谱库中对应多光谱波段λ的植被反射率和土壤反射率;
上述实测地物波谱库是由不同时间、不同地区、不同类型的地物(植被、水体、土壤等)实地测量组成。
(3.3)Path反射率
在DDV暗目标选取区,利用公式⑩计算多光谱波段λ的大气路径反射率:
ρPath-DDV(λ)=ρTOA-DDV(λ)-ρDDV(λ)⑩
(4)根据Path反射率和Rayleigh反射率比值与气溶胶光学厚度的关系,确定多光谱各波段的气溶胶光学厚度;并依据多光谱各波段的气溶胶光学厚度和Rayleigh散射光学厚度,确定出550nm处总的光学厚度τ;
(4.1)多光谱波段λ的气溶胶光学厚度
利用DDV暗目标选取区Path路径反射率和Rayleigh散射反射率的比值与气溶胶光学厚度关系,通过公式求出对应多光谱波段λ的气溶胶光学厚度:
其中:τAer(λ)对应多光谱波段λ的气溶胶光学厚度;
b0(λ)、b1(λ)、b2(λ)对应多光谱波段λ的Path路径反射率和Rayleigh散射反射率的比值与气溶胶光学厚度关系的系数。
(4.2)550nm处的气溶胶光学厚度τ550
由上式确定出550nm处的气溶胶光学厚度τ550:
(4.3)550nm处的Rayleigh散射光学厚度τRay(550)
其中:f0(550)、f1(550)、f2(550)对应550nm处Rayleigh散射光学厚度与海拔高度H大气关系的系数。
(4.4)大气总的光学厚度τ
由上可知550nm处大气总的光学厚度τ为:
(5)依据上述550nm处总的光学厚度,利用查找表获取反演参数ρpath(λ)、T(λ)、S(λ),通过辐射传输方程求解地表反射率,完成大气订正。
上述查找表的建立为利用MODTRAN辐射传输模型,依据几何参数、大气模式、气溶胶模式等设置建立的;查找表的索引τ分别与ρpath(λ)、T(λ)、S(λ)之间的关系满足下式:
ρPath(λ)=c0(λ)+c1(λ)τ+c2(λ)τ2
S(λ)=d0(λ)+d1(λ)τ+d2(λ)τ2
T(λ)=e0(λ)+e1(λ)τ+e2(λ)τ2
c0(λ)、c1(λ)、c2(λ)、d0(λ)、d1(λ)、d2(λ)、e0(λ)、e1(λ)、e2(λ)对应多光谱波段λ在确定的卫星成像几何参数太阳天顶角θS、卫星观测天顶角θV、太阳与卫星的方位角下的大气路径反射率ρpath(λ)、大气透过率T(λ)、半球反射率S(λ)与总光学厚τ的系数。
(6)对步骤(5)中输出的大气订正产品进行精度检验
利用与其他同类卫星反演、其他方法反演和实地测量等多种方法进行对比验证,以确保大气订正的精度,如不符合精度,就从步骤(2)开始重新进行,符合精度输出大气订正产品。
附图3是植被的表观反射率、订正反射率与实测地表反射率曲线对比,从图比较可知,表观反射率虽然经过了绝对辐射校正,但它是地面反射率和大气反射率的总和,不能反映地物的特性,而经过本发明大气订正的反射率与实测地表反射率非常接近误差非常小,确保了订正地表反射率真实反映地物的辐射特性,这使大气订正后的图像更加适宜在地物识别、定量反演和模拟中应用。
本发明未详细说明部分属于本领域技术人员公知常识。
Claims (2)
1.一种遥感卫星多光谱数据的大气订正方法,其特征在于步骤如下:
(1)读取遥感卫星多光谱数据,对该多光谱数据进行预处理,将多光谱数据转换为表观辐亮度数据和表观反射率数据;
(2)依据表观反射率数据计算决策树法中的决策因子:比值植被指数RVI、土壤可调植被指数SAVI和归一化水体指数NDWI,通过决策树法实现DDV暗目标自动提取;
(3)根据步骤(1)中的表观反射率数据获取上述提取的DDV暗目标区域的DDV表观反射率数据;并根据该DDV表观反射率数据结合DDV地表反射率计算Path反射率;
所述的DDV地表反射率的计算公式如下:
ρDDV(λ)=Svegρveg(λ)+(1-Sveg)ρsoil(λ)
其中:ρDDV(λ)对应多光谱波段λ的DDV地表反射率;
ρveg(λ)、ρsoil(λ)是实测多光谱波段λ的植被反射率和土壤反射率;
Sveg是植被调整因子,
(4)根据Path反射率和Rayleigh反射率比值与气溶胶光学厚度的关系,确定多光谱各波段的气溶胶光学厚度;并依据多光谱各波段的气溶胶光学厚度和Rayleigh散射光学厚度,确定出550nm处总的光学厚度τ;
所述Path反射率和Rayleigh反射率比值与气溶胶光学厚度的关系为:
τAer(λ)对应多光谱波段λ的气溶胶光学厚度;
b0(λ)、b1(λ)、b2(λ)对应多光谱波段λ的Path路径反射率和Rayleigh散射反射率的比值与气溶胶光学厚度关系的系数;
(5)依据上述550nm处总的光学厚度,利用查找表获取反演参数ρpath(λ)、T(λ)、S(λ),通过辐射传输方程求解地表反射率,完成大气订正。
2.根据权利要求1所述的一种遥感卫星多光谱数据的大气订正方法,其特征在于所述步骤(5)中查找表的索引τ分别与ρpath(λ)、T(λ)、S(λ)之间的关系满足下式:
ρPath(λ)=c0(λ)+c1(λ)τ+c2(λ)τ2
S(λ)=d0(λ)+d1(λ)τ+d2(λ)τ2
T(λ)=e0(λ)+e1(λ)τ+e2(λ)τ2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101196453A CN102288956B (zh) | 2011-05-10 | 2011-05-10 | 一种遥感卫星多光谱数据的大气订正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011101196453A CN102288956B (zh) | 2011-05-10 | 2011-05-10 | 一种遥感卫星多光谱数据的大气订正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102288956A CN102288956A (zh) | 2011-12-21 |
CN102288956B true CN102288956B (zh) | 2013-04-03 |
Family
ID=45335524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011101196453A Active CN102288956B (zh) | 2011-05-10 | 2011-05-10 | 一种遥感卫星多光谱数据的大气订正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102288956B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650128A (zh) * | 2020-06-08 | 2020-09-11 | 电子科技大学 | 一种基于地表反射率库的高分辨率大气气溶胶反演方法 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267533B (zh) * | 2013-05-08 | 2015-09-30 | 中国科学院光电研究院 | 一种实用型高光谱遥感图像大气自校正方法 |
CN106446307B (zh) * | 2015-08-05 | 2020-01-14 | 中国科学院遥感与数字地球研究所 | 基于气溶胶地基数据的aod垂直订正效果评价方法及*** |
CN105953921B (zh) * | 2016-04-15 | 2018-11-06 | 北京航空航天大学 | 气溶胶参数差异条件下对地观测辐射图像的快速仿真方法 |
KR101877173B1 (ko) * | 2017-05-04 | 2018-07-10 | 경북대학교 산학협력단 | 위성영상을 이용한 해안선 검출 시스템 및 방법 |
CN108256186B (zh) * | 2018-01-04 | 2020-07-10 | 中国科学院遥感与数字地球研究所 | 一种在线计算查找表的逐像元大气校正方法 |
CN108763756B (zh) * | 2018-05-28 | 2021-06-25 | 河南工业大学 | 一种气溶胶光学厚度与pm2.5反演订正方法及其*** |
CN109029735B (zh) * | 2018-07-06 | 2020-06-12 | 湖南文理学院 | 一种地表比辐射率计算方法及装置 |
CN109472237B (zh) * | 2018-11-02 | 2020-12-22 | 中国资源卫星应用中心 | 一种可见光遥感卫星影像的大气订正方法和*** |
CN110779876B (zh) * | 2019-11-07 | 2022-12-02 | 长光禹辰信息技术与装备(青岛)有限公司 | 一种疫木识别方法、装置、设备及计算机可读存储介质 |
CN111521562A (zh) * | 2020-03-19 | 2020-08-11 | 航天信德智图(北京)科技有限公司 | 基于Sentinel-2卫星的棉花植被指数遥感检测方法 |
CN111795936B (zh) * | 2020-08-03 | 2021-11-12 | 长安大学 | 一种基于查找表的多光谱遥感影像大气校正***、方法及存储介质 |
CN112665829B (zh) * | 2020-12-16 | 2022-12-27 | 航天恒星科技有限公司 | 一种光学遥感卫星的波段间定标方法 |
CN113066057B (zh) * | 2021-03-17 | 2023-04-14 | 云南电网有限责任公司电力科学研究院 | 一种气溶胶光学厚度监测方法 |
CN114280694B (zh) * | 2021-12-17 | 2023-05-05 | 南京信息工程大学 | 一种基于气象卫星光谱成像仪的快速辐射传输方法及*** |
CN114943142B (zh) * | 2022-04-29 | 2023-11-28 | 中国科学院空天信息创新研究院 | 高光谱地表反射率和大气参数一体化反演方法及装置 |
CN117347282B (zh) * | 2023-08-22 | 2024-05-28 | 中南大学 | 星基气溶胶光学厚度反演方法、装置及***和存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE59910670D1 (de) * | 1999-07-16 | 2004-11-04 | Deutsch Zentr Luft & Raumfahrt | Verfahren zur Korrektur atmosphärischer Einflüsse in multispektralen optischen Fernerkundungsdaten |
CN100362318C (zh) * | 2005-07-05 | 2008-01-16 | 华东师范大学 | 航空高光谱遥感反演边界层气溶胶光学厚度的大气校正法 |
CN100394212C (zh) * | 2005-08-19 | 2008-06-11 | 广州地理研究所 | 一种大面积作物种植面积及其产量的遥感检测估算方法 |
JP4240132B2 (ja) * | 2007-04-18 | 2009-03-18 | 株式会社デンソー | 内燃機関の制御装置 |
CN101598543B (zh) * | 2009-07-29 | 2011-01-19 | 中国科学院对地观测与数字地球科学中心 | 一种实用的遥感影像大气校正方法 |
-
2011
- 2011-05-10 CN CN2011101196453A patent/CN102288956B/zh active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111650128A (zh) * | 2020-06-08 | 2020-09-11 | 电子科技大学 | 一种基于地表反射率库的高分辨率大气气溶胶反演方法 |
CN111650128B (zh) * | 2020-06-08 | 2021-07-13 | 电子科技大学 | 一种基于地表反射率库的高分辨率大气气溶胶反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102288956A (zh) | 2011-12-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102288956B (zh) | 一种遥感卫星多光谱数据的大气订正方法 | |
CN109581372B (zh) | 一种生态环境遥感监测方法 | |
CN101598543B (zh) | 一种实用的遥感影像大气校正方法 | |
CN106407656A (zh) | 一种基于高分辨率卫星影像数据的气溶胶光学厚度反演方法 | |
CN102435586B (zh) | 地表反照率产品的生成方法及*** | |
Rautiainen et al. | Coupling forest canopy and understory reflectance in the Arctic latitudes of Finland | |
CN109145494B (zh) | 一种海面温度反演方法及*** | |
CN113553907A (zh) | 一种基于遥感技术的森林生态环境状况评价方法 | |
CN102194127A (zh) | 一种多频率sar数据农作物遥感分类方法 | |
Patil et al. | Estimation of mangrove carbon stocks by applying remote sensing and GIS techniques | |
CN114218786B (zh) | 非偏振卫星传感器在轨偏振辐射特性反演方法 | |
CN114782834A (zh) | 基于比值法的农田土壤水分反演方法及*** | |
Jin et al. | Application of QuickBird imagery in fuel load estimation in the Daxinganling region, China | |
CN105403201B (zh) | 一种基于像元分解的遥感图像大气程辐射获取方法 | |
CN106198460B (zh) | 应用于海洋环境的随机起伏海表面反照率全谱段计算方法 | |
CN102636439A (zh) | 一种确定小麦植株吸氮量核心波段的方法 | |
He et al. | Retrieval of rugged mountainous areas Land Surface Temperature from high-spatial-resolution thermal infrared remote sensing data | |
Li | Dynamic monitoring algorithm of natural resources in scenic spots based on MODIS Remote Sensing technology | |
Shah et al. | Evaluation of land surface temperature and normalized difference vegetation index relationship using landsat 8 Satellite images in Mehar Taluka, Dadu | |
CN101319896B (zh) | 一种测量林窗立体结构的方法 | |
Daxiang et al. | A dynamic threshold cloud detecting approach based on the brightness temperature from fy-2 vissr data | |
Chen et al. | Forest NPP estimation based on MODIS data under cloudless condition | |
Mu et al. | An assessment of atmospheric correction methods for GOCI images in the Yellow River estuary | |
CN114354510B (zh) | 一种基于极轨气象卫星mersi的白天云相态反演方法 | |
Zheng et al. | A feasible atmospheric correction method to TM image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |