专利名称：太阳直射光遥感数字图像的辐射变换虚拟仿真方法
太阳直射光遥感数字图像的辐射变换虚拟仿真方法技术领域
本发明属于遥感技术与测绘学领域，尤其是一种太阳直射光遥感数字图像的辐射 变换虚拟仿真方法。
背景技术：
随着卫星遥感技术的高速发展，卫星图像在众多领域有了广泛的应用。由于卫星 图像采集的周期性，影像采集的时间相对固定，很难根据实际需求限时采集图像。采用遥感 数字图像的地形变换方法，虽然可以消除地形起伏带来的地表光辐射差异与山体阴影的影 响，但是仍然没有从根本上解决卫星无法从任意方向、任意时间探测地面的难题。发明内容
本发明的目的在于克服现有技术的不足，提供一种太阳直射光遥感数字图像的辐 射变换虚拟仿真方法，本发明将太阳虚拟为天然的探照灯，模拟太阳从任意方向照射地面 生成的对应遥感图像，从而解决卫星从任意方向、任意时间情况下探测地面的问题。
本发明解决现有的技术问题是采取以下技术方案实现的
一种太阳直射光遥感数字图像的辐射变换虚拟仿真方法，包括以下步骤
(1)根据数字地形图条件和卫星遥感图像条件，计算太阳直射光改正系数、天空散 射光改正系数、散射辐射比，并根据坡面像元太阳直射光、水平地面太阳直射光与对应像元 遥感值之间的关系，获得像元的遥感值和大气程辐射遥感值；
(2)根据卫星遥感图像及太阳直射光改正系数、天空散射光改正系数、散射辐射比 生成水平地面太阳直射光遥感图像和水平地面天空散射光遥感图像；
(3)根据任意太阳位置和任意地形条件，逐点计算任意太阳位置和任意地形条件 下的遥感值，并生成任意太阳位置和任意地形条件下的虚拟遥感数字图像，从而实现数字 地形图上任意太阳方向的地形辐射变换功能。
而且，所述步骤(1)中的数字地形图条件包括地形高度、地形坡度和地形坡向，所 述步骤(1)中的卫星遥感图像条件包括太阳高度角、太阳方位角和像元遥感值。
而且，所述步骤O)中的太阳位置包括太阳高度角和方位角所述步骤O)中的地 形条件包括地形高度、地形坡度和地形坡向。
而且，所述的步骤(3)中计算任意太阳位置和任意地形条件下的遥感值的方法包 括如下步骤
(1)计算虚拟坡面上像元的太阳直射光遥感值Dfsu
DN、SiJ = (DNirDAiJ) / (Fij^Lij X Gij)
(2)计算虚拟坡面上像元的太阳直射光改正系数
ij = 1-tg α Jj Xctg θ ~ X cos (AL~ ^.-Aij)
其中，α ,^Aij分别为该像元的地形坡度和地形坡向； [σοι ] Θ -^arij分别虚拟的太阳高度角和太阳入射的方位角；
(3)计算虚拟的太阳直射光遥感数字图像值)(DNsiJ
XDNsij = DN、SiJ X F、= (DNirDAiJ) X F、υ/ (Fij^Lij X Gij)
(4)计算虚拟坡面上像元的天空散射光遥感值Df diJ
DN、diJ = (DNirDAiJ) Lij/ (Fij^Lij X Gij)
(5)计算虚拟坡面上像元的天空散射光改正系数G— ,j
G ij = Gij ；
(6)计算虚拟的天空散射光遥感数字图像值HMdij
XDNdij = DN、diJ X G、= (DNirDAiJ) Lij X G、(Fij^Lij X Gij)
(7)计算任意太阳位置和任意地形条件下的虚拟遥感数字图像值^(DNij
XDNij = XDNSiJ+XDNdiJ
= (DNirDAiJ) (F、^.+Lij X G、J / (Fij^Lij X Gij)
上述公式中
i = 1，2，……，m ;j = 1，2，……，η ;m，η是该遥感数字图像最大纵向元数和最大横向元数；
Fij为太阳直射光改正系数，Gij为天空散射光改正系数，Lij为散射辐射比；为虚拟的太阳直射光水平变换系数，为虚拟的天空散射光水平变换系 数；
DNij为像元的遥感值，DNAiJ为大气程辐射遥感值。
本发明的优点和积极效果是
本发明将太阳虚拟为天然的探照灯，通过模拟仿真方法实现了遥感数字图像的任 意地形和任意太阳位置之间的相互转换，解决了卫星无法从任意方向、任意时间探测地面 的难题，避免了卫星遥感图像上出现低亮度和永久阴影的现象，完成具有严格解析意义上 的遥感数字图像地形和辐射的虚拟与仿真功能。
图1是本发明处理流程图2是某地区数字高程模型示意图3是某地区太阳光从东南方向照射的遥感图像；
图4是图3利用图2进行虚拟该地区太阳从西北方向照射的遥感图像。
具体实施方式
以下结合附图对本发明实施例做进一步详述。
一种太阳直射光遥感数字图像的辐射变换虚拟仿真方法，如图1所示，包括如下 步骤
1、根据数字地形图条件(包括地形高度、地形坡度和地形坡向)和卫星遥感图像 条件(包括太阳高度角、太阳方位角和像元遥感值)，计算太阳直射光改正系数Fij、天空散 射光改正系数Gu、散射辐射比Lu，并根据坡面像元太阳直射光、水平地面太阳直射光与对 应像元遥感值之间的关系，获得像元的遥感值DNij和大气程辐射遥感值DNau
2、根据卫星遥感图像及太阳直射光改正系数Fij、天空散射光改正系数Gij、散射辐射比Lij生成水平地面太阳直射光遥感图像和水平地面天空散射光遥感图像；
3、根据任意太阳位置(包括太阳高度角和方位角)和任意地形条件(包括地形高 度、地形坡度和地形坡向)，逐点计算任意太阳位置和任意地形条件下的遥感值，并生成任 意太阳位置和任意地形条件下的虚拟遥感数字图像，从而实现数字地形图上任意太阳方向 的地形辐射变换功能。其具体实现方法包括如下步骤
(1)计算虚拟坡面上像元的太阳直射光遥感值DNsij
DN、SiJ = (DNirDAiJ) / (Fij^Lij X Gij)
其中山队」、DNAiJ分别为像元的遥感值和大气程辐射遥感值；
(2)计算虚拟坡面上像元的太阳直射光改正系数Γ ,j
= 1-tga ^ij Xctg θ ^ij X cos ω
由于 ω = AL^。·-Α~ ,
因此F~ ij = 1-tg a \jX ctg θ \j X cos ω ~
= 1-tg a \ j X ctg θ \j X cos (AL~ ΓΑ~
在本次发明方法中前后虚拟的地形模型不变，α = α⑴A^ij = Aij
所以F~ij = 1-tg a Xctg θ ~ X cos (AL~ ^.-Aij)
其中，a ^aij分别为该像元的地形坡度和地形坡向；
θ AIZij分别虚拟的太阳高度角和太阳入射的方位角；
α\_、θ ^ij分别为该像元虚拟的地形坡度和太阳高度角，Α\、ΑΙΛ」、ω ^ij分别为 该像元虚拟的地形坡向、太阳入射的方位角以及它们之间的夹角；
(3)计算虚拟的太阳直射光遥感数字图像值)(DNsiJ
XDNsij = DN、SiJ X F、= (DNirDAiJ) X F、υ/ (Fij^Lij X Gij)
(4)计算虚拟坡面上像元的天空散射光遥感值Df diJ
DN、diJ = (DNij-DAij) Lij/ (Fij^Lij X Gij)
(5)计算虚拟坡面上像元的天空散射光改正系数G— ,j
虚拟数字地面模型上像元点的天空散射光改正系数，在本次发明方法中前后 虚拟的地形模型不变，因此G— ,j = Gij ；
(6)计算虚拟的天空散射光遥感数字图像值)(DNdiJ
XDNdij = DN、diJ X G、= (DNirDAiJ) Lij X G、(Fij^Lij X Gij)
(7)计算任意太阳位置和任意地形条件下的虚拟遥感数字图像值^(DNij
XDNij = XDNSiJ+XDNdiJ = DN~ SiJ X F\j+DN~ diJ X G^ij = ( Ν Γ Α^.) XF^ij/ (FyLij X Gu)+ (DNu-Dau) Lij X G^y(FyLuXGu) = (DNu-Dau) (F^j+L.J X G^ij) / (WGij)
其中i = l，2，……，m;j = l，2，……，η ;m，η是该遥感数字图像最大纵向元数和最大横向元数。
在计算虚拟阴影像元的遥感值时，由于，XDNsij = 0，因此，将上述虚拟遥感数字图 像值计算公式进行推导得到
XDNij = XDNsiJ+XDNdiJ = XDNdij
= DN' diJXG\j
= DNdij X G、^./Gij = (DNirDAiJ) Lij X G、(Fij^LijGij)
本发明的工作原理数字地形图上的卫星太阳直射光遥感数字图像辐射变换虚拟 与仿真本质上是在严格遵守卫星成像规律的原则下，实现遥感数字图像水平地形变换的逆 过程-将水平地面的遥感数字图像转换为任何地形(非水平)上的遥感数字图像，相当于 在一个虚拟的数字地形图上叠加一幅经过水平地形变换的卫星遥感数字图像后，再利用这 个虚拟数字地面模型和太阳位置参数对水平地面上的太阳直射光遥感图像各点遥感值进 行逐点辐射变换的过程，最终生成任意太阳位置和任意地形条件下的虚拟遥感数字图像。 如图2至图4所示，图2为某地区数字高程模型示意图，图3是某地区太阳光从东南方向照 射的遥感图像，经模拟仿真后生成了太阳从西北方向照射的虚拟遥感数字图像。
本发明所述的实施例是说明性的，而不是限定性的，因此发明并不限于具体实施 方式中所述的实施例，凡是由本领域技术人员根据本发明的技术方案得出的其他实施方 式，同样属于本发明保护的范围。
权利要求
1.一种太阳直射光遥感数字图像的辐射变换虚拟仿真方法，其特征在于包括以下步骤(1)根据数字地形图条件和卫星遥感图像条件，计算太阳直射光改正系数、天空散射光 改正系数、散射辐射比，并根据坡面像元太阳直射光、水平地面太阳直射光与对应像元遥感 值之间的关系，获得像元的遥感值和大气程辐射遥感值；(2)根据卫星遥感图像及太阳直射光改正系数、天空散射光改正系数、散射辐射比生成 水平地面太阳直射光遥感图像和水平地面天空散射光遥感图像；(3)根据任意太阳位置和任意地形条件，逐点计算任意太阳位置和任意地形条件下的 遥感值，并生成任意太阳位置和任意地形条件下的虚拟遥感数字图像，从而实现数字地形 图上任意太阳方向的地形辐射变换功能。
2.根据权利要求1所述的太阳直射光遥感数字图像的辐射变换虚拟仿真方法，其特 征在于所述步骤(1)中的数字地形图条件包括地形高度、地形坡度和地形坡向，所述步骤 (1)中的卫星遥感图像条件包括太阳高度角、太阳方位角和像元遥感值。
3.根据权利要求1所述的太阳直射光遥感数字图像的辐射变换虚拟仿真方法，其特征 在于所述步骤O)中的太阳位置包括太阳高度角和方位角所述步骤O)中的地形条件包 括地形高度、地形坡度和地形坡向。
4.根据权利要求1所述的太阳直射光遥感数字图像的辐射变换虚拟仿真方法，其特征 在于所述的步骤(3)中计算任意太阳位置和任意地形条件下的遥感值的方法包括如下步 骤(1)计算虚拟坡面上像元的太阳直射光遥感值Dfsu DN~SiJ = (DNirD奶)/( +、X Gij)(2)计算虚拟坡面上像元的太阳直射光改正系数 ij = 1-tg α Jj Xctg θ ~ Jj X cos (AL~ ΓΑ~其中，α ^Aij分别为该像元的地形坡度和地形坡向； θ \、AL\.分别虚拟的太阳高度角和太阳入射的方位角；(3)计算虚拟的太阳直射光遥感数字图像值)(DNsU: XDNsij = DN^sijXF^ij = (DNirDAiJ) XF^ij/(Fij^LijXGij)(4)计算虚拟坡面上像元的天空散射光遥感值 DN^dij = (DNirDAij)IV(FiZLijXGij)(5)计算虚拟坡面上像元的天空散射光改正系数 G ij = Gij ；(6)计算虚拟的天空散射光遥感数字图像值)(DNdu: XDNdij = DN、diJXG\. = (DNirDAiJ) LijXG^ij/(Fij^LijXGij)(7)计算任意太阳位置和任意地形条件下的虚拟遥感数字图像值 XDNij = XDNSiJ+XDNdiJ=(DNirDAiJ) (F-^LijXG^j)/(Fi^LijXGij) 上述公式中i = 1,2,……，m;j = l，2，……，η ;m，η是该遥感数字图像最大纵向元数和最大横 向元数；Fij为太阳直射光改正系数，Gij为天空散射光改正系数，Lij为散射辐射比；Fij 为虚拟的太阳直射光水平变换系数，G^ij 为虚拟的天空散射光水平变换系数；DNij为像元的遥感值，DNAiJ为大气程辐射遥感值。
全文摘要
本发明涉及一种太阳直射光遥感数字图像的辐射变换虚拟仿真方法，包括以下步骤(1)计算太阳直射光改正系数、天空散射光改正系数、散射辐射比，获得像元的遥感值和大气程辐射遥感值；(2)生成水平地面太阳直射光遥感图像和水平地面天空散射光遥感图像；(3)逐点计算任意太阳位置和任意地形条件下的遥感值，并生成任意太阳位置和任意地形条件下的虚拟遥感数字图像，从而实现数字地形图上任意太阳方向的地形辐射变换功能。本发明将太阳虚拟为天然的探照灯，通过模拟仿真方法实现了遥感数字图像的任意地形和任意太阳位置之间的相互转换，解决了卫星无法从任意方向、任意时间探测地面的难题，完成具有严格解析意义上的遥感数字图像地形和辐射的虚拟与仿真。
文档编号G06T17/05GK102034269SQ20101058293
公开日2011年4月27日 申请日期2010年12月10日 优先权日2010年12月10日
发明者彭玲, 李先华, 池天河 申请人:天津中科遥感信息技术有限公司