卫星遥感图像大气校正的光谱分析方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106339990A(43)申请公布日2017.01.18(21)申请号201610662270.8(22)申请日2016.08.13(71)申请人东北师范大学地址130024吉林省长春市人民大街5268号(72)发明人陈春张树清卢珊韩阳赵建军瞿瑛(74)专利代理机构长春市东师专利事务所22202代理人刘延军李荣武(51)Int.Cl.G06T5/00(2006.01)权利要求书2页说明书6页附图3页(54)发明名称卫星遥感图像大气校正的光谱分析方法(57)摘要本发明属于卫星遥感技术领域，具体涉及卫星遥感图像大气校正的光谱分析方法。本发明完全根据图像自身光谱信息进行大气校正，不需要DTM校正后的图像地形起伏区与平坦区精度相当,不需要同步大气测量数据或地物反射率测量数据，校正精度明显优于目前流行的大气校正软件校正的图像。输出的大气图像及地物图像均为反射率，可合成真彩色图像，用于定量分析。完整保留了原图像的地形信息及BRDF信息。对巨型幅面图像进行分区校正拼接成整幅无接痕，可研发成通用大气校正软件。本发明是中国国家自然科学基金(“遥感图像大气校正光谱矢量分析模型”，批准号:41371394)项目的科研成果。CN106339990ACN106339990A权利要求书1/2页1.卫星遥感图像大气校正的光谱分析方法，其特征在于包含以下步骤：(1).主要利用以下三项原理及步骤编写大气校正程序：①卫星遥感图像定标成反射率图像：ρ*:大气上界观测的反射率,L:卫星观测的辐射亮度[W/(m2·μm·sr)]，θ:太阳高度角，E0:大气上界太阳常数,d:成像时日地距离，以日地平均距离为单位,可表示为：d＝1.0-0.01673cos{π/180[0.9856(j-4)]}(2)j:从1月1日到成像日的天数,含成像日。若j＜4，则取j＝j+365.2422，QuickBird各波段太阳常数,单位:W/(m2·μm)波段Pan1234E01462.65721935.66871853.67661584.87111129.95292E0′(λ):日地平均距离时大气上界太阳辐照度光谱[W/(m·μm)]，Δλ:遥感器波段有效宽度(μm)，f(λ):QuickBird遥感器光谱响应效率函数，无量纲；②反射率图像分离成地物图像和大气图像：卫星遥感图像由大气图像和地物图像叠加形成，大气图像位于上层，地物图像位于下层，光谱分析法基于以下几点事实将反射率图像分离成地物图像和大气图像：a.大气反射率随波长的增加迅速减少；b.图像中的地物光谱不受地形影响，不受观测几何条件太阳高度角，太阳方位角，卫星观测天顶角，观测方位角的影响；c.大气图像没有明显边界，表现为逐渐过渡，地物一般有明显边界；d.阴影中的入射光(漫射光)与直射光比较，短波成分增加，卫星观测的波段i图像定标成反射率图像记作根据以上事实满足以下方程，用这些方程将分离成地物图像ρi和大气图像ρa(i)：波段i(地气总反射率，地物反射率，大气反射率)令：令：2CN106339990A权利要求书2/2页已知图像中的一个典型值且可求取(γ,K,F)的对应值(i，j)取(蓝:B，红:R)波段，根据上式计算蓝波段大气反射率ρa(B)，其它波段大气反射率应采用下式计算:ρa(j)是最终输出的大气图像，供大气遥感使用；③地物图像大气校正大气散射和吸收消减了地物的入射光，吸收导致的波段i大气透射率计作Ti，校正消减后波段i地物反射率图像记作e＝2.71828459045是最终输出的地物图像，供地物遥感使用；(2).对待校正的卫星遥感图像需要进行以下几项预处理：①图像辐射校正，误差小于5％；②图像几何校正；③图像波段间配准，配准误差小于0.5个像元；(3).利用步骤(1)所编写的卫星遥感图像大气校正程序进行卫星遥感图像大气校正。3CN106339990A说明书1/6页卫星遥感图像大气校正的光谱分析方法技术领域[0001]本发明属于卫星遥感技术领域，具体涉及卫星遥感图像大气校正的光谱分析方法。背景技术[0002]典型情况下，卫星遥感图像来自地物反射辐射贡献在450nm波段占50％，在850nm波段占80％，经过大气散射和吸收与大气图像叠置在一起，导致图像中地物模糊不清、光谱失真，构建模型校正图像大气影响以获得清晰的地物图像、获得精确的地物反射率数据，国内外学者进行了40多年的研究，取得了可喜的成果，成果纳入遥感理论体系，写入专著，优秀成果被研发成计算机程序收集到了大气校正软件中，得到了广泛应用,例如,目前应用最广的3套大气校正程序：[0003]FLAASH程序[0004]FLAASH(FastLine-of-sightAtmosphericAnalysis