[发明专利]宽视域卫星影像地表反射率反演方法

摘要

太阳照射天顶角和卫星观测天顶角是卫星遥感影像大气校正和地表反射率反演中很重要的输入参数。针对视场角较大的影像，本发明逐像元计算太阳照射天顶角和卫星观测天顶角，然后反演地表反射率，提出一种宽视域陆地卫星遥感影像地表反射率获取的方法。

主权项

一种宽视域卫星影像地表反射率反演的方法，包括如下步骤：第一步、对卫星影像进行几何校正；第二步、辐射定标，将影像像元亮度值转换为星上辐射亮度Lλ；Lλ＝Gain×DN+BiasLλ为波段λ的星上辐射亮度，DN为波段λ的像元值，Gain和Bias分别为波段λ的增益和偏置；第三步、获取大气程辐射LpLp＝Gain×DNmin+Bias‑0.01E0cosθzTzTv/πd2Lp是大气程辐射，DNmin是暗目标DN值，E0是大气层外相应波长的太阳光谱辐照度，θz是太阳天顶角，Tz是太阳照射方向的大气透过率，Tv是卫星观测方向的大气透过率，d是日地距离；第四步、逐像元计算太阳天顶角(4‑1)计算太阳时：t=ts+0.170sin(4π(J-80)373)-0.129sin(2π(J-8)355)+12(SM-L)π]]>其中t是太阳时，ts是标准时，SM是该时区标准经线的经度，L是该像元点的经度，J是儒略日；(4‑2)计算太阳赤纬：δ=0.4093sin(2π(J-81)368)]]>其中δ是太阳赤纬，J是儒略日；(4‑3)计算太阳天顶角θz=π2-arcsin(sinlsinδ+coslcosδcosπt12)]]>其中θz是太阳天顶角，l是该像元点的纬度，δ是太阳赤纬，t是太阳时；第五步、在陆地卫星影像星下点成像情况下，按照以下步骤逐像元计算卫星观测天顶角：(5‑1)获取任一像元点E与影像中心点O的行距Δx和列距Δy；(5‑2)由行距和列距计算像元点E与影像中心点O的距离ss=(Δx)2+(Δy)2;]]>(5‑3)由s和卫星高度h之间的三角关系获取该像元点的卫星观测天顶角θvθv＝arctan(s/h)；第六步、计算大气光学厚度τ(6‑1)根据陆地卫星遥感影像的经纬度和影像成像时间，获取同一地理范围、同步过境的MODIS 0.47μm和0.66μm的气溶胶光学厚度、总可降水汽、地表气压和臭氧浓度影像，并将各影像做几何校正，转换为和陆地卫星遥感影像一致的投影方式和分辨率，且影像行列数一致；(6‑2)由MODIS总可降水汽w计算大气水蒸汽光学厚度；τw=0.2385awλwM(1+20.07awλwM)0.45]]>其中，τw是大气水蒸汽光学厚度，awλ是水汽吸收系数，w是总可降水汽，M是相对大气量；(6‑3)计算臭氧吸收光学厚度τo＝Cozone×Aozone(λ)τo是臭氧吸收光学厚度，Cozone为MODIS臭氧浓度，Aozone(λ)为波段λ的臭氧吸收系数；(6‑4)由MODIS 0.47μm和0.66μm的气溶胶光学厚度，利用下式计算波长指数α及大气浑浊度系数βα=-lnτa(λ1)-lnτa(λ2)lnλ1-lnλ2]]>β＝τa(λ1)λ1α＝τa(λ2)λ2ατa为气溶胶光学厚度，λ1和λ2分别为0.47μm和0.66μm，τa(λ1)和τa(λ2)分别为0.47μm和0.66μm的气溶胶光学厚度；在得到α和β后，可以利用下式得到任意波长λ的气溶胶光学厚度τa(λ)：τa(λ)＝βλ‑α其中，λ为波长，τa(λ)是波长λ的气溶胶光学厚度；(6‑5)计算瑞利散射光学厚度；τr=p1013.250.008569λ-4(1+0.0133λ-2+0.00013λ-4)]]>其中，τr为瑞利散射光学厚度，λ为波长；p为MODIS的地表气压产品；(6‑6)τ＝τw+τo+τa+τrτ是大气光学厚度；第七步、计算太阳照射方向和卫星观测方向的大气透过率；Tz＝exp(‑τ/cosθz)＝exp{(‑τr‑τa‑τo‑τw)/cosθz}Tv＝exp(‑τ/cosθv)＝exp{(‑τr‑τa‑τo‑τw)/cosθv}其中，Tz和Tv分别是太阳照射方向和卫星观测方向的大气透过率，τ是大气光学厚度，θz是太阳天顶角，θv是卫星观测天顶角，τr、τa、τo和τw分别是瑞利散射光学厚度、气溶胶光学厚度、臭氧吸收光学厚度和大气水蒸汽光学厚度；第八步、计算地表反射率ρρ=π(Lλ-Lp)d2E0cosθzTzTv]]>d＝1+0.0167sin[2π(J‑93.5)/365]E0是大气层外相应波长的太阳光谱辐照度，d是日地距离，J是儒略日。