一种纠正卫星遥感数据辐射误差的方法技术

本发明专利技术公开了一种纠正卫星遥感数据辐射误差的方法，包括以下步骤：(A)对目标地区700米以上高山区各季节的大气散射光进行光谱测量，对目标地区的山区植被和清洁深水体进行光谱测量；(B)从目标地区的原始卫星图像数据提取高山区植被阴影像元和清洁深水体像元；(C)消除可见光波段的辐射误差；(D)消除红外波段的辐射误差。该方法可对光学卫星遥感数据的辐射误差进行有效纠正，使得各波段之间的相对辐射误差的纠正精度更高。

【技术实现步骤摘要】
一种纠正卫星遥感数据辐射误差的方法
本专利技术属于水质遥感数据处理方法
，特别涉及一种纠正卫星遥感数据辐射误差的方法。
技术介绍
遥感卫星发射上天后，由于温度、压力等环境因素改变等原因，传感器的响应函数常发生一定的变化，造成了卫星遥感得到的辐射值与实际辐射值之间的误差。目前解决该问题的方法主要有以下两种：一是：在卫星上测量太阳直射光和标准光源作为定标参数来对遥感数据进行辐射纠正。但是，由于多种原因，如卫星上物理化学条件与地面相差甚大，标准灯的辐射能量不可避免地存在一定的波动等，造成星上辐射纠正后仍然存在一定误差。二是：通过在地面与卫星同步测量地物光谱，然后根据测得的地物光谱对遥感数据进行辐射纠正。但是，由于地面测得的地面目标反射率与卫星测得的反射率相差大气的影响，且由于大气气溶胶的变化，大气的影响难以精确测量，故地面定标也常存在明显的误差。目前，国外卫星遥感原始数据提供的定标参数进行辐射纠正后，反射率误差仍达1％～2％；国产卫星的辐射误差更达2％～4％。该误差对水质遥感等方面的定量计算产生了极严重的影响，成为提高定量遥感精度的一大障碍。同时，该误差限制了水质遥感等定量遥感方法的应用。
技术实现思路
本专利技术的目的是提供一种纠正卫星遥感数据辐射误差的方法，该方法可对光学卫星遥感数据的辐射误差进行有效纠正，使得各波段之间的相对辐射误差的纠正精度更高。为了达到上述目的，本专利技术采用如下技术方案：一种纠正卫星遥感数据辐射误差的方法，包括以下步骤：(A)对目标地区700米以上高山区各季节的大气散射光进行光谱测量，得到大气散射光反射率ρ，对目标地区的山区植被和清洁深水体进行光谱测量，得到植被可见光波段反射率Rv和清洁水体红外波段出水反射率Rwp；(B)从目标地区的原始卫星图像数据提取高山区植被阴影像元和清洁深水体像元，得到原始图像植被阴影像元表观反射率R′1和原始图像清洁深水体像元表观反射率R′2；(C)根据高山区相应季节的大气散射光反射率ρ、植被可见光波段反射率Rv和原始图像植被阴影像元表观反射率R′1，计算消除可见光波段的辐射误差；(D)根据高山区相应季节的大气散射光反射率ρ、清洁水体红外波段出水反射率Rwp和原始图像清洁深水体像元表观反射率R′2，计算消除红外波段的辐射误差。作为一种具体的实施例，步骤(C)具体为：(C1)根据高山区相应季节的大气散射光反射率ρ和植被可见光波段反射率Rv，计算出植被阴影像元的可见光波段真实表观反射率Rvd；(C2)根据所述植被阴影像元的可见光波段真实表观反射率Rvd和所述原始图像植被阴影像元表观反射率R′1，计算出可见光波段的辐射误差ε1；(C3)根据所述可见光波段的辐射误差ε1和所述原始图像植被阴影像元表观反射率R′1，消除可见光波段的辐射误差。进一步地，步骤(C1)中植被阴影像元的可见光波段真实表观反射率Rvd的计算公式为：式中：ρ为测量得到的大气散射光反射率；为大气散射系数，P(θ)为大气散射相函数，θ为散射角，T＝e-ω为大气透过率，Rv为植被反射率。进一步地，步骤(C2)中可见光波段的辐射误差ε1的计算公式为：ε1＝R′1-Rvd，式中：R′1为像元原始表观反射率。进一步地，步骤(C3)中消除可见光波段的辐射误差的计算公式为：R1＝R′1-ε1。作为一种具体的实施例，步骤(D)具体为：(D1)根据高山区相应季节的大气散射光反射率ρ和清洁水体红外波段出水反射率Rwp，计算出清洁深水体像元的红外波段真实表观反射率Rw；(D2)根据所述清洁深水体像元的红外波段真实表观反射率Rw和原始图像清洁深水体像元表观反射率R′2，计算出红外波段的辐射误差ε2；(D3)根据所述红外波段的辐射误差ε2和原始图像清洁深水体像元表观反射率R′2，消除红外波段的辐射误差。进一步地，步骤(D1)中计算出清洁深水体像元的红外波段真实表观反射率Rw的计算公式为：Rw＝Rwp+ρ。进一步地，步骤(D2)中计算出红外波段的辐射误差ε2的计算公式为：ε2＝R′2-Rwp进一步地，步骤(D3)中消除红外波段的辐射误差的计算公式为：R2＝R′2-ε2。与现有技术相比，本专利技术的有益效果在于：本专利技术一种纠正卫星遥感数据辐射误差的方法，通过首先获取目标地区的高山区各个季节的大气散射光反射率、目标地区的山区植被和清洁深水体的植被可见光波段反射率、清洁水体红外波段出水反射率，以及目标地区的原始图像植被阴影像元表观反射率和原始图像清洁深水体像元表观反射率；然后根据获取到的数据，计算消除可见光波段的辐射误差，进一步地，计算消除红外波段的辐射误差。对光学卫星遥感数据的辐射误差进行有效纠正，使得各波段之间的相对辐射误差的纠正效果，尤其是低值区的辐射误差的纠正效果尤为显著。可将低值区的辐射误差减小近一个数量级。为提高定量遥感的精度打下了基础。本专利技术独立于传感器，可对各种光学传感器的遥感数据进行纠正。对国产卫星数据的使用意义尤为突出。经本方法进行辐射纠正后，环境卫星、中巴卫星和资源卫星等国产卫星数据均可用于水质遥感，并取得很好的结果。附图说明图1是本专利技术的纠正卫星遥感数据辐射误差的方法流程图。图2是步骤(C)的具体方法流程图。图3是步骤(D)的具体方法流程图。图4是珠江三角洲地区的原始中巴卫星的图像。图5是利用本专利技术方法纠正后的珠江三角洲地区的中巴卫星的图像。图6(a)—图6(c)是存在严重辐射误差的中巴卫星遥感数据的第二、三和四波段的纠正前图像。图6(d)—图6(f)是存在严重辐射误差的中巴卫星遥感数据的第二、三和四波段的纠正后图像。具体实施方式下面将结合附图以及具体实施方法来详细说明本专利技术，在本专利技术的示意性实施及说明用来解释本专利技术，但并不作为对本专利技术的限定。在有数字高程模型(DigitalElevationModel，简称：DEM)的情况下，为了使辐射误差纠正的精度更加精准，本专利技术定量考虑目标地区的地形对辐射误差纠正的影响。为了进一步使辐射误差纠正的精度更加精准，本专利技术优选以天空光为主的特征地物像元进行辐射纠正。因此，本专利技术分别定量考虑目标地区的山区植被和清洁深水体对辐射误差纠正的影响。同时，根据相同季节、晴天的情况下，高山区天空散射光稳定的特点，本专利技术对目标地区高山区各季节的大气散射光进行光谱测量，优选地，对目标地区700米以上高山区进行大气散射光光谱测量。如图1所示，本专利技术一种纠正卫星遥感数据辐射误差的方法，包括以下步骤：步骤(A)对目标地区700米以上高山区各季节的大气散射光进行光谱测量，得到大气散射光反射率ρ，对目标地区的山区植被和清洁深水体进行光谱测量，得到植被可见光波段反射率Rv和清洁水体红外波段出水反射率Rwp。具体地，用ASD(便携式地物光谱仪)等连续光谱仪测量目标地区700米以上高山区各季节的大气散射光、测量目标地区的植被可见光波段反射率和清洁水体红外波段出水反射率。步骤(B)从目标地区的原始卫星图像数据提取高山区植被阴影像元和清洁深水体像元，得到原始图像植被阴影像元反射率R′1和原始图像清洁深水体像元表观反射率R′2。优选地，此处需要提取的高山区植被阴影像元和清洁深水体像元为典型像元，也即纯粹的植被阴影像元和纯粹的清洁水体像元。步骤(本文档来自技高网...

【技术保护点】
一种纠正卫星遥感数据辐射误差的方法，其特征在于，包括以下步骤：(A)对目标地区700米以上高山区各季节的大气散射光进行光谱测量，得到大气散射光反射率ρ，对目标地区的山区植被和清洁深水体进行光谱测量，得到植被可见光波段反射率Rv和清洁水体红外波段出水反射率Rwp；(B)从目标地区的原始卫星图像数据提取高山区植被阴影像元和清洁深水体像元，得到原始图像植被阴影像元表观反射率R′1和原始图像清洁深水体像元表观反射率R′2；(C)根据高山区相应季节的大气散射光反射率ρ、植被可见光波段反射率Rv和原始图像植被阴影像元表观反射率R′1，计算消除可见光波段的辐射误差；(D)根据高山区相应季节的大气散射光反射率ρ、清洁水体红外波段出水反射率Rwp和原始图像清洁深水体像元表观反射率R′2，计算消除红外波段的辐射误差。

【技术特征摘要】
1.一种纠正卫星遥感数据辐射误差的方法，其特征在于，包括以下步骤：(A)对目标地区700米以上高山区各季节的大气散射光进行光谱测量，得到大气散射光反射率ρ，对目标地区的山区植被和清洁深水体进行光谱测量，得到植被可见光波段反射率Rv和清洁水体红外波段出水反射率Rwp；(B)从目标地区的原始卫星图像数据提取高山区植被阴影像元和清洁深水体像元，得到原始图像植被阴影像元表观反射率R′1和原始图像清洁深水体像元表观反射率R′2；(C1)根据高山区相应季节的大气散射光反射率ρ和植被可见光波段反射率Rv，计算出植被阴影像元的可见光波段真实表观反射率Rvd，其计算公式为：式中：为大气散射系数，P(θ)为大气散射相函数，θ为散射角，T＝e-ω为大气透过率；(C2)根据所述植被阴影像元的可见光波段真实表观反射率Rvd和所述原始图像植被阴影像元表观反射率R′1，计算出可见光波段的辐射误差ε1；(C3)根据所述可见光波段的辐射误差ε1和所述原始图像植被阴影像元表观反射率R′1，消除可见光波段的辐射误差；(D)根据高山区相应季节的大气散射光反射率ρ、清洁水体红外波段出水反射率Rwp和原始图像清洁深水体像元表观反射率R′2，计算消除红外波段的辐射误差。2.根据权利要求1所述的纠正卫星遥感数据辐射误差的方法，其特征在于，步...

【专利技术属性】
技术研发人员：秦雁，邓孺孺，梁业恒，何颖清，陈启东，熊龙海，刘旭拢，刘英飞，卢世军，刘永明，林梨，
申请(专利权)人：中山大学，
类型：发明
国别省市：广东;44