CN105678236B - 一种陆地植被冠层偏振反射建模方法 - Google Patents
一种陆地植被冠层偏振反射建模方法 Download PDFInfo
- Publication number
- CN105678236B CN105678236B CN201511021375.7A CN201511021375A CN105678236B CN 105678236 B CN105678236 B CN 105678236B CN 201511021375 A CN201511021375 A CN 201511021375A CN 105678236 B CN105678236 B CN 105678236B
- Authority
- CN
- China
- Prior art keywords
- canopy
- polarization
- polarized reflectance
- vegetation
- soil
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/188—Vegetation
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种陆地植被冠层偏振反射建模方法，其步骤如下：在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式；依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数；基于几何光学原理，分别计算单个叶片和下垫面土壤的偏振反射率；对冠层进行等间距分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率；分析植被冠层偏振反射率的分布特性。本发明能够实现不同类型植被、太阳位置及观测位置等条件下，陆地植被冠层的偏振反射率特性的精确描述，是分析植被偏振方向反射特性及其影响因素的有效技术方法，可以为地气系统偏振信息解耦提供有效的技术方案。
Description
(一)技术领域
本发明涉及一种陆地植被冠层偏振反射建模方法，属于光学遥感领域，在对植被偏振反射特性分析和地气偏振信息解耦研究方面具有重要的理论意义和应用价值。
(二)背景技术
偏振探测具有广泛的应用前景，目前大气偏振遥感已经开展了探索性应用，植被偏振的研究较少。植被偏振信息作为陆地地表偏振信息的重要组成部分，不仅能够直接用于地表参数反演，而且是实现地气偏振信息解耦的关键，然而，目前已有的植被偏振模型过于简化，并不能够充分描述植被偏振特性，尤其是忽略了入射天空光的偏振特性，这在很大程度上制约了植被偏振遥感应用的发展。
本发明首先建立上半球空间的天空光偏振度分布模式，再分别建立植被单个叶片的偏振反射率模型、下垫面土壤的偏振反射率模型、植被冠层的偏振反射率模型，实现大气——植被——土壤耦合系统的偏振反射建模方法。冠层的方向反射率分解为单次散射和多次散射之和，多次散射不具备偏振特性，因此只计算单次散射产生的偏振反射率，单次散射包括冠层单次散射和土壤单次散射，建模过程中将植被冠层沿竖直方向平均分为水平方向均一，竖直方向变化的若干层，土壤表面单独看作一层，采用累加法计算得到整个冠层的偏振反射率；陆地植被偏振反射建模，对于对地遥感中正确计算植被偏振反射特性，实现陆地植被的分类、参数反演与地气偏振信息解耦具有重要的科学意义与应用价值。
(三)发明内容
本发明涉及一种陆地植被冠层偏振反射建模方法，技术解决方案是：在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式；依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数；基于几何光学原理，代入偏振菲涅尔反射系数，分别构建单个叶片和下垫面土壤的偏振反射率模型；对植被冠层进行等间距分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率；分析植被冠层偏振反射率的分布特性。
1一种陆地植被冠层偏振反射建模方法，其特征在于包含以下步骤：
(1)在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式；
(2)依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数；
(3)基于几何光学原理，代入偏振菲涅尔反射系数，分别构建单个叶片和下垫面土壤的偏振反射率模型；
(4)对冠层进行等间距分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率；
(5)分析植被冠层偏振反射率的分布特性。
2步骤(1)中所述的“在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式”，具体实现过程如下：
①输入太阳入射方向在上半球空间对天空光的散射方向进行等立体角采样，得到采样点处大气分子的散射方向计算散射角γ
其中θs为太阳天顶角，为太阳方位角，θ为散射天顶角，为散射方位角；
②采用半经验瑞利散射模型计算天空光偏振度δ
其中δmax为特定太阳天顶角时的最大偏振度，计算方法为
3步骤(2)中所述的“依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数”，具体实现过程如下：
①依据输入的太阳方向和观测方向计算入射角α
其中θv为观测天顶角，为观测方位角；
②计算偏振菲涅尔反射系数F
F(α,n)＝κsrs 2-κprp 2
其中rs为反射光垂直于入射面的分量，rp为反射光平行于入射面的分量，计算方法为
其中α为入射角，n为被测目标和空气的折射率比；
κs和κp分别为垂直分量和平行分量占总能量的比例，且κs+κp＝1。
4步骤(3)中所述的“基于几何光学原理，代入偏振菲涅尔反射系数，分别构建单个叶片和下垫面土壤的偏振反射率模型”，具体实现过程如下：
①将叶片表面看作很多准镜面的微面元，计算微面元法线方向的天顶角θh
叶片镜面反射存在几何影响和面元之间的相互遮挡，计算衰减因子K
其中
②假设叶片表面的微面元法向服从高斯分布，计算单个叶片的偏振反射率ρleafp
其中nl为植被叶片和空气的折射率比，σ为叶片表面粗糙因子；
③计算下垫面土壤的偏振反射率ρsoilp
其中ns为土壤和空气的折射率比，A是与土壤表面粗糙度有关的自由参数。
5步骤(4)中所述的“对冠层进行等间距分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率”，具体实现过程如下：
①给定地面与冠层顶部间任意位置高度z，计算该位置处的二向孔隙率函数Q
其中H为冠层高度，ps(z)·pv(z)是冠层中高度z处能被太阳照亮，同时能由冠层外自由观测到的联合概率，G(Ωs)和G(Ωv)分别为单位面积叶片在太阳入射方向和观测方向上的投影，S(z)为高度z处的单位叶面积指数；
②将植被冠层沿竖直方向平均分为N层，每层介质的光学特性视为水平方向均一，
竖直方向变化，土壤表面单独看作一层，计算第i个植被单位层的高度z
其中N为层数，则第i个植被层对整个冠层的单次散射贡献为
其中为叶倾角分布概率密度函数，Q(z)为第i个植被层处的二向孔隙率函数；
土壤层对整个冠层的单次散射贡献为
其中Q(H)为冠层底部的二向孔隙率函数；
③分别计算N-1个植被层及下垫面土壤层的偏振反射率，采用累加法计算整个冠层的偏振反射率。
6步骤(5)中所述的“分析植被冠层偏振反射率的分布特性”，具体实现过程如下：对观测方向进行等立体角采样，根据设定的地表类型，输入叶片、土壤和冠层光学及结构特性参数，计算每个采样点的植被冠层偏振反射率，分析偏振反射率在上半球空间的分布特性。
本发明与现有技术相比的优点在于：
(1)目前学术界缺少有效的考虑入射光偏振特性的描述陆地植被偏振反射特性的技术方法，本发明采用累加法计算叶片和下垫面土壤对整个冠层偏振反射率的贡献，具有显著的创新性，丰富了植被偏振光学遥感领域的研究。
(2)本发明考虑了大气分子散射作用引起的天空光偏振特性、土壤的粗糙因子与遮挡效应、植被的空间结构特征等因素，使得模型模拟结果更加接近真实情况，能适用于不同植被覆盖度下的植被冠层或土壤的偏振特性模拟，具有良好的应用适用性。
(四)附图说明
图1为本发明的总体技术流程图。
(五)具体实施方式
为了更好地说明本发明涉及的一种陆地植被冠层偏振反射建模方法，利用本发明的建模方法进行植被冠层偏振反射特性进行分析。本发明一种陆地植被冠层偏振反射建模方法，实现流程如图1所示，具体实施步骤如下：
步骤(1)在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式：
①输入正午时刻太阳天顶角θs和太阳方位角在上半球空间对天空光的散射方向进行等立体角采样，得到1500个采样点处的散射天顶角θ和散射方位角计算散射角γ
②采用半经验瑞利散射模型计算天空光偏振度δ
其中δmax为特定太阳天顶角时的最大偏振度，当太阳天顶角为30°时，δmax取值63％，步骤(2)依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数：
①输入观测天顶角θv和观测方位角计算入射角α
②计算偏振菲涅尔反射系数F
F(α,n)＝κsrs 2-κprp 2
其中rs为反射光垂直于入射面的分量，rp为反射光平行于入射面的分量，计算方法为
其中α为入射角，n为被测目标和空气的折射率比；
κs和κp分别为垂直分量和平行分量占总能量的比例，以自然光入射为例，κs＝κp＝0.5。
步骤(3)基于几何光学原理，代入偏振菲涅尔反射系数，分别构建单个叶片和下垫面土壤的偏振反射率模型：
①将叶片表面看作很多准镜面的微面元，计算微面元法线方向的天顶角θh
叶片镜面反射存在几何影响和面元之间的相互遮挡，计算衰减因子K
其中
②假设叶片表面的微面元法向服从高斯分布，计算单个叶片的偏振反射率ρleafp
其中nl为植被叶片和空气的折射率比，σ为叶片表面粗糙因子；
③计算下垫面土壤的偏振反射率ρsoilp
其中ns为土壤和空气的折射率比，A是与土壤表面粗糙度有关的自由参数。
步骤(4)对冠层进行等间距分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率：
①给定地面与冠层顶部间任意位置高度z，计算该位置处的二向孔隙率函数Q
其中H为冠层高度，ps(z)·pv(z)是冠层中高度z处能被太阳照亮，同时能由冠层外自由观测到的联合概率，G(Ωs)和G(Ωv)分别为单位面积叶片在太阳入射方向和观测方向上的投影，S(z)为高度z处的单位叶面积指数；
②将植被冠层沿竖直方向平均分为N层，每层介质的光学特性视为水平方向均一，竖直方向变化，土壤表面单独看作一层，计算第i个植被单位层的高度z
其中N为层数，则第i个植被层对整个冠层的单次散射贡献为
其中为叶倾角分布概率密度函数，Q(z)为第i个植被层处的二向孔隙率函数；
土壤层对整个冠层的单次散射贡献为
其中Q(H)为冠层底部的二向孔隙率函数；
③分别计算N-1个植被层及下垫面土壤层的偏振反射率，采用累加法计算整个冠层的偏振反射率。
步骤(5)分析植被冠层偏振反射率的分布特性：对观测方向进行等立体角采样，根据设定的地表类型，输入叶片、土壤和冠层光学及结构特性参数，分别计算每个采样点处太阳光和天空光入射时的冠层偏振反射率，进行加权求和，最后分析偏振反射率在上半球空间的分布特性。
Claims (5)
1.一种陆地植被冠层偏振反射建模方法，其特征在于包含以下步骤：
(1)在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式；
(2)依据天空光偏振度的分布模式，计算部分偏振光入射的菲涅尔反射系数；
(3)基于几何光学原理，代入偏振菲涅尔反射系数，分别建立单个叶片和下垫面土壤的偏振反射率模型；
(4)对植被冠层进行等间距垂向分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率；
(5)分析植被冠层偏振反射率的分布特性。
2.根据权利要求1 所述的一种陆地植被冠层偏振反射建模方法，其特征在于：步骤(1)中所述的“在上半球空间对天空光的散射方向进行采样，构建天空光的偏振度分布模式”，具体实现过程如下：
①输入太阳入射方向在上半球空间对天空光的散射方向进行等立体角采样，得到采样点处大气分子的散射方向计算散射角γ
其中θs为太阳天顶角，为太阳方位角，θ为散射天顶角，为散射方位角；
②采用半经验瑞利散射模型计算天空光偏振度δ
其中δmax为特定太阳天顶角时的最大偏振度，计算方法如下：
3.根据权利要求1所述的一种陆地植被冠层偏振反射建模方法，其特征在于：步骤(3)中所述的“基于几何光学原理，代入偏振菲涅尔反射系数，分别建立单个叶片和下垫面土壤的偏振反射率模型”，具体实现过程如下：
①将叶片表面看作很多准镜面的微面元，计算微面元法线方向的天顶角θh
叶片镜面反射存在几何影响和面元之间的相互遮挡，计算衰减因子K
其中
②假设叶片表面的微面元法向的服从高斯分布，计算单个叶片的偏振反射率ρleafp
其中nl为植被叶片和空气的折射率比，σ为叶片表面粗糙因子；θs为太阳天顶角，为太阳方位角，θv为观测天顶角，为观测方位角，F(α,nl)为偏振菲涅尔反射系数，α为相对入射角；
③计算下垫面土壤的偏振反射率ρsoilp
其中ns为土壤和空气的折射率比，A是与土壤表面粗糙度有关的自由参数。
4.根据权利要求1所述的一种陆地植被冠层偏振反射建模方法，其特征在于：步骤(4)中所述的“对植被冠层进行等间距垂向分层，逐层计算叶片的偏振反射率，土壤单独看作一层，采用累加法计算整个冠层的偏振反射率”，具体实现过程如下：
①给定地面与冠层顶部间任意位置高度z，计算该位置处的二向孔隙率函数Q
其中H为冠层高度，ps(z)·pv(z)是冠层中高度z处能被太阳照亮，同时能由冠层外自由观测到的联合概率，G(Ωs)和G(Ωv)分别为单位面积叶片在太阳入射方向和观测方向上的投影，S(z)为高度z处的单位叶面积指数；
②将植被冠层沿竖直方向平均分为N层，每层介质的光学特性视为水平方向均一，竖直方向变化，土壤表面单独看作一层，计算第i个植被单位层的高度z
其中N为层数，则第i个植被层对整个冠层的单次散射贡献为
其中为叶倾角分布概率密度函数，Q(z)为第i个植被层处的二向孔隙率函数，θs为太阳天顶角，为太阳方位角，θv为观测天顶角，为观测方位角，ρleafp为叶片偏振反射率；
土壤层对整个冠层的单次散射贡献为
其中Q(H)为冠层底部的二向孔隙率函数，ρsoilp为土壤偏振反射率；
③分别计算N-1个植被层及下垫面土壤层的偏振反射率，采用累加法计算整个冠层的偏振反射率。
5.根据权利要求1所述的一种陆地植被冠层偏振反射建模方法，其特征在于：步骤(5)中所述的“分析植被冠层偏振反射率的分布特性”，具体实现过程如下：对观测方向进行等立体角采样，根据设定的地表类型，输入叶片、土壤和冠层光学及结构特性参数，计算每个采样点的植被冠层偏振反射率，分析偏振反射率在上半球空间的分布特性。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511021375.7A CN105678236B (zh) | 2015-12-31 | 2015-12-31 | 一种陆地植被冠层偏振反射建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511021375.7A CN105678236B (zh) | 2015-12-31 | 2015-12-31 | 一种陆地植被冠层偏振反射建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105678236A CN105678236A (zh) | 2016-06-15 |
| CN105678236B true CN105678236B (zh) | 2018-11-20 |
Family
ID=56298120
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511021375.7A Expired - Fee Related CN105678236B (zh) | 2015-12-31 | 2015-12-31 | 一种陆地植被冠层偏振反射建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105678236B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2740871C2 (ru) * | 2019-05-31 | 2021-01-21 | Российская Федерация, от имени которой выступает Государственная корпорация по космической деятельности "РОСКОСМОС" | Способ определения коэффициентов отражения исследуемой поверхности и автоматизированный комплекс для его осуществления |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108318458B (zh) * | 2017-01-16 | 2020-10-09 | 北京航空航天大学 | 一种适用于不同天气条件下的室外典型地物pBRDF的测量方法 |
| CN107688003B (zh) * | 2017-09-04 | 2020-06-30 | 南京大学 | 一种消除植被冠层结构和地表背景影响的叶片反射率卫星遥感提取方法 |
| CN108009392B (zh) * | 2017-10-19 | 2022-01-04 | 桂林航天工业学院 | 一种浓密植被地表的遥感反射率模型构建及标定应用方法 |
| CN109740180B (zh) * | 2018-11-30 | 2021-02-02 | 北京航空航天大学 | 一种基于几何光学模型的稀疏水生植被冠层方向反射率模拟方法 |
| CN110686776B (zh) * | 2019-10-09 | 2021-05-04 | 北京航空航天大学 | 一种室内高光谱遥感成像仿真装置与方法 |
| CN114676386A (zh) * | 2022-02-28 | 2022-06-28 | 中国人民解放军火箭军工程大学 | 一种基于光学微面元理论几何衰减因子修正方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102538783A (zh) * | 2012-02-14 | 2012-07-04 | 北京大学 | 基于遥感天空偏振模式图的仿生导航方法及导航定位系统 |
| CN103398957A (zh) * | 2013-08-12 | 2013-11-20 | 河海大学 | 基于高光谱与激光雷达提取叶面积垂直分布的方法 |
| CN103632040A (zh) * | 2013-11-14 | 2014-03-12 | 北京航空航天大学 | 一种通用的水生植被辐射传输模型 |
-
2015
- 2015-12-31 CN CN201511021375.7A patent/CN105678236B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102538783A (zh) * | 2012-02-14 | 2012-07-04 | 北京大学 | 基于遥感天空偏振模式图的仿生导航方法及导航定位系统 |
| CN103398957A (zh) * | 2013-08-12 | 2013-11-20 | 河海大学 | 基于高光谱与激光雷达提取叶面积垂直分布的方法 |
| CN103632040A (zh) * | 2013-11-14 | 2014-03-12 | 北京航空航天大学 | 一种通用的水生植被辐射传输模型 |
Non-Patent Citations (2)
| Title |
|---|
| Polarized reflectance of bare soils and vegetation: measurements and models;Breon F M, Tanre D, Lecomte P;《Polarized reflectance of bare soils and vegetation: measurements and models》;19951231;全文 * |
| 基于垂直观测的植被冠层高光谱偏振反射特性研究;吕云峰;《光谱学与光谱分析》;20130430;全文 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2740871C2 (ru) * | 2019-05-31 | 2021-01-21 | Российская Федерация, от имени которой выступает Государственная корпорация по космической деятельности "РОСКОСМОС" | Способ определения коэффициентов отражения исследуемой поверхности и автоматизированный комплекс для его осуществления |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105678236A (zh) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105678236B (zh) | 一种陆地植被冠层偏振反射建模方法 | |
| CN104699952B (zh) | 一种湿地水生植被冠层brdf蒙特卡洛模型 | |
| Kanda et al. | A simple energy balance model for regular building arrays | |
| CN106248003B (zh) | 一种三维激光点云提取植被冠层聚集度指数的方法 | |
| CN103630495B (zh) | 一种水生植被-大气耦合辐射传输模型 | |
| CN102914501A (zh) | 一种利用激光点云计算三维森林冠层消光系数的方法 | |
| Zhang et al. | Doubling–adding method for delta-four-stream spherical harmonic expansion approximation in radiative transfer parameterization | |
| CN103940515B (zh) | 一种波浪水面天空漫反射光偏振场模拟方法 | |
| CN106874621B (zh) | 一种针叶植被冠层反射率计算方法 | |
| CN106446307B (zh) | 基于气溶胶地基数据的aod垂直订正效果评价方法及系统 | |
| CN101915914A (zh) | 一种基于查找表的遥感影像逐像元大气校正方法 | |
| CN106909750B (zh) | 一种阔叶植被冠层反射率的计算方法 | |
| CN103400364B (zh) | 一种森林资源变化监测方法 | |
| CN103632040A (zh) | 一种通用的水生植被辐射传输模型 | |
| CN111595801B (zh) | 一种整层大气气溶胶棕色碳成分的遥感识别和估计方法 | |
| CN101936777A (zh) | 一种基于热红外遥感反演近地层气温的方法 | |
| CN106225693A (zh) | 一种细粒子气溶胶光学厚度和气溶胶类型同时反演方法 | |
| CN109740180A (zh) | 一种基于几何光学模型的稀疏水生植被冠层方向反射率模拟方法 | |
| CN116452055B (zh) | 融合表面网格模型和深度学习的城市光伏潜力计算方法 | |
| CN101876700B (zh) | 一种基于辐射度的复杂地形区域辐射传输模拟方法 | |
| CN116822141A (zh) | 利用卫星微光遥感反演夜间大气气溶胶光学厚度的方法 | |
| CN105302980B (zh) | 一种基于sar数据的城区空气动力学粗糙度反演方法 | |
| CN113870147B (zh) | 一种考虑阴影区域的遥感图像brdf校正方法 | |
| CN111709180B (zh) | 耦合吸收性气溶胶模型的大气矢量辐射传输模拟方法 | |
| Mao et al. | Dust particle size distribution inversion based on the multi population genetic algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181120 Termination date: 20191231 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |