一种植被总初级生产力的遥感反演方法

1.本发明涉及遥感技术领域，更具体的说是涉及一种植被总初级生产力的遥感反演方法。


背景技术：

2.总初级生产力是指在单位时间和单位面积上，绿色植物通过光合作用所固定的有机碳总量，包括自养呼吸和异养呼吸。总初级生产力是定量表征碳汇/碳源的关键参数，代表了特定区域或生态系统的生产力，同时也反映了区域植被或生态系统的生长特征和健康状况。对植被生产力研究根据不同精度的要求，主要的方法有生物量调査法、涡度相关通量观测法和模型估算法。目前，利用模型估算区域生态系统生产力已成为一种重要而广泛接受的研究方法。
3.gpp估算模型可分为统计模型、过程模型和光能利用率模型，其中光能利用率模型结构简单，数据较少且易获取，并且结果具有一定的可靠性，因此被广泛使用。光能利用率模型结构较为固定，在不同的模型中同一影响因子具有不同的表达形式，并且其估算效果也有差异，所以利用单一参数对植物总生产例进行估算时，其代表性不足，容易造成误差。


技术实现要素：

4.有鉴于此，本发明提供了一种植被总初级生产力的遥感反演方法，可以避免由单一参数代表性不足带来的误差，提高反演精度。
5.为了实现上述目的，本发明采用如下技术方案：
6.一种植被总初级生产力的遥感反演方法，包括以下步骤：
7.将不同的光能利用率模型中表示同一影响因子但表示形式不同的参数，采用多元线性方法组合为一个参数；
8.基于不同影响因子线性组合后的参数构建植被总初级生产力估算模型；
9.基于所述植被总初级生产力估算模型估算植被总初级生产力。
10.进一步的，所述影响因子包括：水分胁迫因子、温度胁迫因子和光合有效吸收比例。
11.进一步的，所述水分胁迫因子的不同表示形式包括：
[0012][0013][0014]
[0015][0016]
其中，w
vpm
、w
ec-lue
、w
casa
和w
mod17
分别表示vpm模型、ec-lue模型、casa模型和mod17模型中的水分胁迫因子；lswi表示地表水分指数，lswi
max
表示地表水分指数最大值；le表示潜热，h表示显热；e表示蒸散发，e
p
表示潜在蒸散发；vpd表示饱和水汽压差，vpd
min
表示最小饱和水汽压差，vpd
max
表示最大饱和水汽压差。
[0017]
进一步的，不同表示形式下的所述水分胁迫因子的线性组合为：
[0018]
f(w)＝a1w
vpm
+a2w
ec-lue
+a3w
casa
+a4w
mod17
+d1[0019]
其中，a1、a2、a3、a4和d1分别为经验参数。
[0020]
进一步的，所述温度胁迫因子的不同表示形式包括：
[0021]
t
casa
＝t
s1
×
t
s2
[0022]
t
s1
＝0.8+0.02
×
t
opt-0.0005
×
t
2opt
[0023][0024][0025][0026][0027][0028]
其中，t
casa
、t
tem
、t
tg
、t
c-fix
、t
mod17
分别表示casa模型、tem模型、tg模型、c-fix模型和mod17模型中的温度胁迫因子；t
s1
、t
s2
表示计算t
casa
的中间参数，t
opt
表示植物进行光合作用的最适温度，t表示近地表气温，t
max
表示植物进行光合作用的最大温度，t
min
表示植物进行光合作用的最小温度，δh
a，p
表示活化能，δh
d，p
表示去活化能，rg表示气体常数，c1表示常数，根据地形不同取不同值，δs表示co2变性平衡的熵，tmin
min
表示在各时间段内能进行光合作用的各个最低温度里的最小值，tmin
max
表示在各时间段内能进行光合作用的各个最低温度里的最大值。
[0029]
进一步的，不同表示形式下的所述温度胁迫因子的线性组合为：
[0030]
f(t)＝b1t
casa
+b2t
tem
+b3t
tg
+b4t
c-fix
+b5t
mod17
+d2[0031]
其中，b1、b2、b3、b4、b5和d2分别为经验参数。
[0032]
进一步的，所述光和有效辐射吸收比例的不同表示形式包括：
[0033]
fpar
c-fix
＝1.638
×
ndvi-0.1426
[0034][0035][0036][0037]
fpar
vpm
＝evi
[0038]
其中，fpar
c-fix
和fpar
vpm
分别表示c-fix模型和vpm模型中的光合有效辐射吸收比例，fpar
casa1
、fpar
casa2
和fpar
casa3
分别表示casa模型中的光合有效辐射吸收比例；ndvi表示归一化植被指数；ndvi
i,min
表示第i种植被的ndvi最小值，ndvi
i,max
表示第i种植被的ndvi最大值，fpar
max
表示光合有效辐射吸收比例的最大值，fpar
min
表示光合有效辐射吸收比例的最小值，sr表示比值植被指数，sr
i,min
表示第i种植被的sr最小值，sr
i,max
表示第i种植被的sr最大值，evi表示增强型植被指数。
[0039]
进一步的，不同表示形式下的所述光和有效辐射吸收比例线性组合为：
[0040]
fpar＝c1fpar
c-fix
+c2fpar
casa1
+c3fpar
casa2
+c4fpar
casa3
+c5fpar
evi
+d3[0041]
其中，c1、c2、c3、c4、c5和d3分别为经验参数。
[0042]
进一步的，所述植被总初级生产力估算模型的表达式为：
[0043]
gpp＝ε0×
f(w)
×
f(t)
×
par
×
fpar
[0044]
其中，ε0表示最大光能利用率，代表生态系统对太阳辐射的利用效率；f(w)表示水分胁迫因子，f(t)表示温度胁迫因子，fpar表示光合有效吸收比例，par表示光合有效辐射。
[0045]
经由上述的技术方案可知，与现有技术相比，本发明具有以下有益效果：
[0046]
由于不同光能利用率模型中对同一影响因子考虑的切入点不同，因此表示方式也不同，估算效果也有所不同。本发明将多个模型中表示同一影响因子但表示形式不同的参数使用多元线性方法将其组合为一个参数，从多个角度对同一影响因子进行表示，减少了由单一参数解释力不够带来的误差，提高了影响因子在估算模型中的精度，对于计算植被总初级生产力具有重要意义。
[0047]
同时，本发明实现方式简单、实用，模型使用的遥感数据易获取、处理技术成熟，所需的输入参数较少。
附图说明
[0048]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
[0049]
图1为本发明提供的植被总初级生产力的遥感反演方法的流程图；
[0050]
图2为本发明提供的制备总初级生产力的计算过程示意图。
具体实施方式
[0051]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0052]
如图1-图2所示，本发明实施例公开了一种植被总初级生产力的遥感反演方法，其特征在于，包括以下步骤：
[0053]
s1、将不同的光能利用率模型中表示同一影响因子但表示形式不同的参数，采用多元线性方法组合为一个参数；
[0054]
s2、基于不同影响因子线性组合后的参数构建植被总初级生产力估算模型；
[0055]
s3、基于所述植被总初级生产力估算模型估算植被总初级生产力。
[0056]
其中，所述影响因子包括：水分胁迫因子、温度胁迫因子和光合有效吸收比例。
[0057]
在一个具体实施例中，对水分胁迫因子、温度胁迫因子和光合有效吸收比例的构建过程和线性组合过程做了详细描述。
[0058]
(1)水分胁迫因子的不同表示形式包括：
[0059][0060][0061][0062][0063]
其中，w
vpm
、w
ec-lue
、w
casa
和w
mod17
分别表示vpm模型(vegetation photosynthesis model，植被光合作用模型)、ec-lue模型(eddy covarianc e-light use efficiency，涡度相关光能利用率模型)、casa模型(carnegie-ames-stanford approach模型)和mod17模型(modis gpp产品算法)中的水分胁迫因子；lswi表示地表水分指数，lswi
max
表示地表水分指数最大值；le表示潜热，h表示显热；e表示蒸散发，e
p
表示潜在蒸散发；vpd表示饱和水汽压差，vpd
min
表示最小饱和水汽压差，vpd
max
表示最大饱和水汽压差。
[0064]
上述不同模型中的水分胁迫因子因关注的生理过程不同，而使用了不同的参数和构造形式，而这些水分胁迫因子效果也各有优劣，因此考虑将其进行线性组合从而达到互补的效果。不同表示形式下的所述水分胁迫因子的线性组合为：
[0065]
f(w)＝a1w
vpm
+a2w
ec-lue
+a3w
casa
+a4w
mod17
+d1[0066]
其中，a1、a2、a3、a4和d1分别为经验参数。
[0067]
(2)温度胁迫因子的不同表示形式包括：
[0068]
t
casa
＝t
s1
×
t
s2
[0069]
t
s1
＝0.8+0.02
×
t
opt-0.0005
×
t
2opt
[0070][0071][0072][0073][0074][0075]
其中，t
casa
、t
tem
、t
tg
、t
c-fix
、t
mod17
分别表示casa模型、tem模型(terrestrial ecosystem model，陆地生态系统模型)、tg模型(temperature and greenness model，温度和绿度模型)、c-fix模型和mod17(modis gpp产品算法)模型中的温度胁迫因子；t
s1
、t
s2
表示计算t
casa
的中间参数，t
opt
表示植物进行光合作用的最适温度，t表示近地表气温，t
max
表示植物进行光合作用的最大温度，t
min
表示植物进行光合作用的最小温度，δh
a，p
表示活化能，δh
d，p
表示去活化能，rg表示气体常数，c1表示常数，根据地形不同而取不同值(无量纲，山地取21.9，平原或丘陵取21.6)，δs表示co2变性平衡的熵(j
·
k-1
mol-1
)，tmin
min
表示在各时间段内能进行光合作用的各个最低温度里的最小值，tmin
max
表示在各时间段内能进行光合作用的各个最低温度里的最大值。
[0076]
不同模型中的温度胁迫因子因关注的生理过程不同而使用了不同的参数和构造形式，而这些温度胁迫因子效果也各有优劣，因此同样考虑将其进行线性组合从而达到互补的效果。
[0077]
不同表示形式下的所述温度胁迫因子的线性组合为：
[0078]
f(t)＝b1t
casa
+b2t
tem
+b3t
tg
+b4t
c-fix
+b5t
mod17
+d2[0079]
其中，b1、b2、b3、b4、b5和d2分别为经验参数。
[0080]
(3)所述光和有效辐射吸收比例的不同表示形式包括：
[0081]
fpar
c-fix
＝1.638
×
ndvi-0.1426
[0082][0083][0084]
[0085][0086]
fpar
vpm
＝evi
[0087]
其中，fpar
c-fix
和fpar
vpm
分别表示c-fix模型和vpm模型(vegetation photosynthesis model，植被光合作用模型)中的光合有效辐射吸收比例，在casa模型中，ndvi和sr都可以用来表示光和有效辐射吸收比例，但是使用ndvi计算的fpar略高于实测值，使用sr计算的fpar略低于实测值，因此casa模型使用这两者的均值来表示光和有效辐射吸收比例，因此分别使用fpar
casa1
、fpar
casa2
和fpar
casa3
表示casa模型中的三种光合有效辐射吸收比例的计算方法；ndvi表示归一化植被指数；ndvi
i,min
表示第i种植被的ndvi最小值，ndvi
i,max
表示第i种植被的ndvi最大值，fpar
max
表示光合有效辐射吸收比例的最大值，为固定值0.95，fpar
min
表示光合有效辐射吸收比例的最小值，为固定值0.001，sr表示比值植被指数，sr
i,min
表示第i种植被的sr最小值，sr
i,max
表示第i种植被的sr最大值，evi表示增强型植被指数。
[0088]
同样，不同模型中的fpar因关注的生理过程不同而使用了不同的参数和构造形式，而这些fpar表示的植被冠层吸收比例效果也各有优劣，因此考虑将其进行线性组合从而达到互补的效果。
[0089]
不同表示形式下的所述光和有效辐射吸收比例线性组合为：
[0090]
fpar＝c1fpar
c-fix
+c2fpar
casa1
+c3fpar
casa2
+c4fpar
casa3
+c5fpar
evi
+d3[0091]
其中，c1、c2、c3、c4、c5和d3分别为经验参数。
[0092]
在一个具体实施例中，所述植被总初级生产力估算模型的表达式为：
[0093]
gpp＝ε0×
f(w)
×
f(t)
×
par
×
fpar
[0094]
其中，ε0表示最大光能利用率，代表生态系统对太阳辐射的利用效率；f(w)表示水分胁迫因子，f(t)表示温度胁迫因子，fpar表示光合有效吸收比例，par表示光合有效辐射。
[0095]
本实施例所涉及到的所有参数容易获取，其中，ε0通过观测站的nee和par数据经michaelis-menten函数拟合得到，ndvi、evi和lswi数据经陆地反射率产品mod09a1计算得出；地表蒸散发数据、地表潜在蒸散发数据和地表温度数据分别使用mod16a2产品和mod11a2产品；par数据使用glass par数据产品；潜热、显热和饱和水汽压差数据通过era5再分析数据计算得出；t
opt
、t
min
、t
max
和tmin
min
数据根据查找表结合植被类型图得到。
[0096]
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0097]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

技术特征：
1.一种植被总初级生产力的遥感反演方法，其特征在于，包括以下步骤：将不同的光能利用率模型中表示同一影响因子但表示形式不同的参数，采用多元线性方法组合为一个参数；基于不同影响因子线性组合后的参数构建植被总初级生产力估算模型；基于所述植被总初级生产力估算模型估算植被总初级生产力。2.根据权利要求1所述的一种植被总初级生产力的遥感反演方法，其特征在于，所述影响因子包括：水分胁迫因子、温度胁迫因子和光合有效吸收比例。3.根据权利要求2所述的一种植被总初级生产力的遥感反演方法，其特征在于，所述水分胁迫因子的不同表示形式包括：分胁迫因子的不同表示形式包括：分胁迫因子的不同表示形式包括：分胁迫因子的不同表示形式包括：其中，w
vpm
、w
ec-lue
、w
casa
和w
mod17
分别表示vpm模型、ec-lue模型、casa模型和mod17模型中的水分胁迫因子；lswi表示地表水分指数，lswi
max
表示地表水分指数最大值；le表示潜热，h表示显热；e表示蒸散发，e
p
表示潜在蒸散发；vpd表示饱和水汽压差，vpd
min
表示最小饱和水汽压差，vpd
max
表示最大饱和水汽压差。4.根据权利要求3所述的一种植被总初级生产力的遥感反演方法，其特征在于，不同表示形式下的所述水分胁迫因子的线性组合为：f(w)＝a1w
vpm
+a2w
ec-lue
+a3w
casa
+a4w
mod17
+d1其中，a1、a2、a3、a4和d1分别为经验参数。5.根据权利要求2所述的一种植被总初级生产力的遥感反演方法，其特征在于，所述温度胁迫因子的不同表示形式包括：t
casa
＝t
s1
×
t
s2
t
s1
＝0.8+0.02
×
t
opt-0.0005
×
t
2optopt
其中，c1、c2、c3、c4、c5和d3分别为经验参数。9.根据权利要求2所述的一种植被总初级生产力的遥感反演方法，其特征在于，所述植被总初级生产力估算模型的表达式为：gpp＝ε0×
f(w)
×
f(t)
×
par
×
fpar其中，ε0表示最大光能利用率，代表生态系统对太阳辐射的利用效率；f(w)表示水分胁迫因子，f(t)表示温度胁迫因子，fpar表示光合有效吸收比例，par表示光合有效辐射。

技术总结
本发明涉及遥感技术领域，具体涉及一种植被总初级生产力的遥感反演方法，包括以下步骤：将不同的光能利用率模型中表示同一影响因子但表示形式不同的参数，采用多元线性方法组合为一个参数；基于不同影响因子线性组合后的参数构建植被总初级生产力估算模型；基于所述植被总初级生产力估算模型估算植被总初级生产力。本发明可以避免由单一参数代表性不足带来的误差，提高反演精度。提高反演精度。提高反演精度。


技术研发人员：唐伯惠 冯思齐 李梦华 陈俊屹 葛中曦
受保护的技术使用者：昆明理工大学
技术研发日：2023.06.13
技术公布日：2023/9/7