基于冠层内辐射传输机理的光学遥感影像地形校正方法技术

本发明专利技术涉及光学遥感影像地形校正技术领域，从光子‑冠层‑地形交互的物理机理出发，通过简化辐射传输方程，提出一种兼顾物理机理完备性和数学形式简洁性的光学遥感影像地形校正方法，包括以下步骤：收集目标的光学遥感影像和DEM数据；根据DEM数据计算坡度、坡向，根据影像头文件，计算逐像元的成像几何；对遥感影像进行大气校正，获得地表反射率；然后计算各像元对应的平地或坡地消光路径；将地表反射率和消光路径作为输入，将坡地反射率校正为平地反射率。本发明专利技术的地形校正方法普适性强，具备业务化运行能力，有助于提高卫星遥感技术在山区的应用潜力。

A topographic correction method for optical remote sensing images based on the mechanism of radiation transmission in the canopy

The present invention relates to the technical field of optical remote sensing image terrain correction, starting from the physical mechanism of the photon canopy terrain interaction, by simplifying the radiative transfer equation, a physical mechanism of both completeness and concise mathematical form correction method for optical remote sensing images of terrain, which comprises the following steps: Optical Remote Sensing images and DEM data were collected according to the target; DEM data to calculate the slope, according to the image header files, calculate the imaging geometry of each pixel of the remote sensing image; atmospheric correction for surface reflectance; flat or sloping path of each pixel corresponds to the extinction and calculation of surface reflectance and extinction; the path as input, the slope reflectance correction to the ground reflectance. The topographic correction method of the invention has the advantages of strong universality and operational capability, and helps to improve the application potential of satellite remote sensing technology in mountain areas.

【技术实现步骤摘要】
基于冠层内辐射传输机理的光学遥感影像地形校正方法
本专利技术涉及光学遥感影像地形校正
，具体涉及一种基于冠层内辐射传输机理的光学遥感影像地形校正方法。
技术介绍
遥感技术因其宏观、动态和客观的优点，在环境监测、防灾减灾、农业估产等众多领域取得了广泛应用。但在地形复杂的山区，地形显著影响光学遥感的成像过程，使遥感观测的辐射信号中包含了过多非地表实况的地形扰动信息。抑制光学遥感信号中的地形信息，获得影像的等效平地反射率(即地形校正)，是山地光学遥感应用的基本预处理技术。地形以多种途径使光学遥感影像产生畸变，如：①地形对太阳入射辐射的再分配；②地形对大气底层辐射传输过程的影响；③地形倾斜对冠层结构的调制。相应地，已有地形校正方法也可以大致分为三类：①入射辐射校正方法，如余弦校正；②地形-大气协同校正方法，如Kobayashi方法；③坡地冠层反射率地形校正方法，如SCS校正。SCS校正是一种常用的坡地冠层反射率校正方法，它认为冠层反射率主要有光照树冠决定，通过考虑地形倾斜引起的光照树冠面积的变化，将地形校正后的反射率表达为：式中，ρt、ρSCS分别是地形校正前后的反射率，α为坡度，θ1为太阳入射高度角，θ1’为坡地上的局地入射高度角，且θ1′＝cosθ1cosα+sinθ1sinαcosβ，式中，β为坡向。相对于余弦校正、C校正等，SCS方法考虑了地形对冠层结构的调制，物理机理更加完备。但是，此方法基于几何光学基本原理，将树冠视为不透光的刚体，并且仅考虑了地形对光照树冠面积的影响，而忽略了阴影树冠、光照背景、阴影背景三分量对冠层反射率的贡献。SCS方法在稀疏植被区的校正精度还存在不足。亟需发展对坡地冠层结构刻画更为准确的反射率校正方法。D-S校正也是一种常用的坡地冠层反射率校正方法，它认为平地和坡地冠层反射率模型均可以表达为以下形式：式中，为入射方向和观测方向之间夹角的函数，θ1、θ2分别为为入射天顶角和观测天顶角。基于以上假设，地形校正之后的反射率可表达为：式中，ρt、ρD-S分别是地形校正前后的反射率，θ’1、θ’2分别为坡地上的局地入射高度角和观测高度角，可由式θ′1＝cosθ1cosα+sinθ1sinαcosβ计算得到。D-S校正物理机理清晰，数学形式简洁，是一种有效的坡地冠层反射率校正方法。但是，坡地冠层内辐射传输是一个高度非线性的物理过程，其满足式的假设过于简单，使D-S校正存在较大的不确定性，限制了它的广泛应用。亟需发展对坡地冠层内辐射传输过程刻画更为细致的反射率校正方法。本专利技术将基于辐射传输基本理论，考虑由地表倾斜引起的冠层内光子传输自由路径长度的改变，据此构建坡地辐射传输方程；最后，通过简化辐射传输方程，得到形式简洁、机理完备的地形校正方法，普适性强，具备业务化运行能力，有助于提高光学卫星遥感技术在山区的应用潜力。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种基于冠层内辐射传输机理的光学遥感影像地形校正方法，从光子-冠层-地形交互的物理机理出发，通过简化辐射传输方程，提出一种兼顾物理机理完备性和数学形式简洁性的坡地冠层反射率归一化算法。本专利技术的目的是通过以下技术方案来实现的：一种基于冠层内辐射传输机理的光学遥感影像地形校正方法，包括以下步骤：S1.收集目标的光学遥感影像和DEM数据；S2.根据DEM数据计算坡度、坡向，根据影像头文件，计算逐像元的成像几何，然后根据下列公式计算各像元对应的平地或坡地消光路径，S(θ)＝1/cos(θ)式中，θ、φ分别为太阳/观测天顶角和方位角，α、β分别为坡度、坡向，S、St分别为相同成像几何条件下平地和坡地的归一化消光路径；同时对遥感影像进行大气校正，获得地表反射率；S3.将地表反射率和消光路径作为输入，利用下列公式将坡地反射率校正为平地反射率，式中，ρt、ρRT分别是地形校正前后的反射率，Ω1为太阳入射方向，Ω2为观测方向。进一步地，所述步骤S2中，利用ENVI、ArcGIS或ReliefVisualizationToolbox软件处理DEM数据计算坡度、坡向。进一步地，所述步骤S2中，影像头文件由USGS网站下载，并据此计算得到逐像元成像几何产品。进一步地，所述步骤S2中，利用6S、MORTRAN或ENVI软件的FLAASH模块完成对遥感影像进行大气校正，获得地表反射率。本专利技术的工作原理：假设坡地冠层是由倾斜的冠层顶和冠层底界定的浑浊介质。在重力影响下，植被具有向地生长特性，单株植被结构，如树冠形状、叶倾角等在水平面(而非倾斜坡面)的方位角上具有旋转不变性，因此适于在水平坐标系，而非局地坡面坐标系上描述冠层内的辐射传输过程。辐射传输方程描述了辐射量在冠层内的传输过程，可表示为：式中，I为特定高度和方向的分谱辐亮度(W/m2·nm·sr)，μ为叶面积体密度(m2/m3)，G(Ω2)为所有朝向的叶片向Ω2方向的平均投影面积，Γ(Ω1→Ω2)为散射相函数，即光子由Ω1方向碰撞叶片后被散射到Ω2方向的概率。z为由冠层顶(z＝0)至冠层底(z＝-1)的归一化距离；在以下假设下：假设I：天空散射光为0，即忽略天空散射光的影响；假设II：传感器接收的信号为植被的单次散射，即土壤背景的反射以及植被自身的多次散射可以忽略；基于以上假设，传感器接收到的辐亮度可以表达为：式中，L(z)为高度z处的累积叶面积指数，s(z,Ω2)为高度z处，沿Ω2方向逃逸出冠层所经过的消光路径长度；F(z,Ω2)为源函数，且式中，J为太阳入射辐照度。根据双向反射因子的定义，可知冠层反射率可表示为：式中，S(Ω1)和S(Ω2)分别为Ω1、Ω2方向相对于天顶方向的归一化消光路径。进一步假设：假设III：叶倾角类型为球型，即G≡0.5；则，由于植被的单体结构特征和叶片光学特性受地形的影响可以近似忽略，所以，式(9)中仅有1-exp[-LG(Ω1)S(Ω1)-LG(Ω2)S(Ω2)]与地形有关，对该项进行泰勒级数展开，并忽略二阶以上的高阶项，得：将平地和坡地反射率都写为式(9)的形式，两者取比值，并进一步整理可得本专利技术坡地反射率校正方法：式中，ρt、ρRT分别是地形校正前后的反射率，S、St分别为相同成像几何条件下平地和坡地的归一化消光路径；且S(θ)＝1/cos(θ)式中，θ、分别为太阳/观测天顶角和方位角，α、β分别为坡度、坡向；可见，山地情况下消光路径不再具有方位角上的旋转不变性，并且受坡度、坡向的影响。本专利技术的有益效果是：本专利技术基于冠层内辐射传输机理的光学遥感影像地形校正方法，考虑了坡地冠层内的辐射传输过程机理，克服了已有方法物理机理不完备的缺点；进一步通过化简辐射传输方程，得到了解析的坡地冠层反射率校正方法；由于本专利技术兼顾了物理机理的完备性和数学形式的简洁性，可以很方便地嵌入已有入射辐射校正方法和地形-大气协同校正方法，构造全链路地形校正方法，有助于显著提高光学遥感技术在山区的应用潜力。附图说明图1为本专利技术光学遥感影像地形校正方法的流程图；图2为本专利技术试验例中校正前的影像图；图3为本专利技术试验例中的DEM数据图；图4为本专利技术试验例中采用本专利技术方法校正后的影像图；图5为本专利技术试验例中采用SCS校正方法校正后的影像图；图6为本专利技术试验例中采用D-S校正方法校正后的影像图。本文档来自技高网...

【技术保护点】
一种基于冠层内辐射传输机理的光学遥感影像地形校正方法，其特征在于，包括以下步骤：S1.收集目标的遥感影像和DEM数据；S2.根据DEM数据计算坡度、坡向，根据影像头文件，计算逐像元的成像几何，然后根据下列公式计算各像元对应的平地或坡地消光路径，S(θ)＝1/cos(θ)

【技术特征摘要】
1.一种基于冠层内辐射传输机理的光学遥感影像地形校正方法，其特征在于，包括以下步骤：S1.收集目标的遥感影像和DEM数据；S2.根据DEM数据计算坡度、坡向，根据影像头文件，计算逐像元的成像几何，然后根据下列公式计算各像元对应的平地或坡地消光路径，S(θ)＝1/cos(θ)式中，θ、φ分别为太阳/观测天顶角和方位角，α、β分别为坡度、坡向，S、St分别为相同成像几何条件下平地和坡地的归一化消光路径；同时对遥感影像进行大气校正，获得地表反射率；S3.将地表反射率和消光路径作为输入，利用下列公式将坡地反射率校正为平地反射率，

【专利技术属性】
技术研发人员：尹高飞，李爱农，
申请(专利权)人：中国科学院，水利部成都山地灾害与环境研究所，
类型：发明
国别省市：四川,51