CN115358095A - 山地高空间分辨率植被总初级生产力估算方法 - Google Patents
山地高空间分辨率植被总初级生产力估算方法 Download PDFInfo
- Publication number
- CN115358095A CN115358095A CN202211278846.2A CN202211278846A CN115358095A CN 115358095 A CN115358095 A CN 115358095A CN 202211278846 A CN202211278846 A CN 202211278846A CN 115358095 A CN115358095 A CN 115358095A
- Authority
- CN
- China
- Prior art keywords
- spatial resolution
- vegetation
- data
- primary productivity
- estimation result
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000029553 photosynthesis Effects 0.000 claims abstract description 43
- 238000010672 photosynthesis Methods 0.000 claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 14
- 239000002689 soil Substances 0.000 claims description 32
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 8
- 238000007637 random forest analysis Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 6
- 238000002310 reflectometry Methods 0.000 claims description 6
- 238000001556 precipitation Methods 0.000 claims description 5
- 230000002776 aggregation Effects 0.000 claims description 4
- 238000004220 aggregation Methods 0.000 claims description 4
- 230000004931 aggregating effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- GVVPGTZRZFNKDS-JXMROGBWSA-N geranyl diphosphate Chemical compound CC(C)=CCC\C(C)=C\CO[P@](O)(=O)OP(O)(O)=O GVVPGTZRZFNKDS-JXMROGBWSA-N 0.000 abstract description 22
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 abstract description 7
- 229910052799 carbon Inorganic materials 0.000 abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 7
- 230000004044 response Effects 0.000 abstract description 6
- 230000007613 environmental effect Effects 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000012876 topography Methods 0.000 description 7
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 230000000243 photosynthetic effect Effects 0.000 description 4
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 102000004190 Enzymes Human genes 0.000 description 2
- 108090000790 Enzymes Proteins 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 210000000473 mesophyll cell Anatomy 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006473 carboxylation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000027756 respiratory electron transport chain Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Image Processing (AREA)
Abstract
本发明属于遥感图像数据处理技术领域，涉及山地高空间分辨率植被总初级生产力估算方法，该方法包括：获取第一输入数据集；处理得到第二输入数据集与第三输入数据集；确定第一GPP估算结果与第二GPP估算结果；构建第一关系模型；计算第三GPP估算结果；计算差值；构建第二关系模型；估算第四GPP估算结果；计算第五植被总初级生产力估算结果；求和得到目标GPP估算结果。本发明采用生态水文模型，能够更好地表征山地生态系统的碳循环和水循环机制，通过光合作用过程对不同空间尺度环境变化的非线性响应过程来改进线性降尺度过程，减小了低分辨率估算空间尺度误差引起的不确定性。
Description
技术领域
本发明属于遥感图像数据处理技术领域，具体而言，涉及山地高空间分辨率植被总初级生产力估算方法。
背景技术
山地生态系统约占陆地表面的比例较大，在气候变化、碳储存、物质资源供应、生物多样性保护和水源涵养发挥着举足轻重的作用。植被作为监测陆地生态系统对气候变化响应的重要生物指标，准确估算山地植被总初级生产力(GPP，生态系统水平上叶片光合作用过程所固定的有机碳总量)对于理解山地生态系统的功能及其对全球气候系统的反馈至关重要。相较于涡度相关技术，遥感技术能够获取区域地表特征信息，可为植被生产力的模拟提供长时间序列大区域范围的地表信息。
基于多种卫星传感器提供不同空间分辨率的植被信息，通过生态系统模型能开展多种空间分辨率的植被GPP(gross primary productivity,植被生态系统总初级生产力)估算，然而，由于计算资源的限制使得绝大多数区域的GPP数据集只能利用LUE（Light useefficiency model, 光能利用效率模型）或VI（vegetation index，植被指数）方法进行生成，缺少像元之间生态过程的利用；同时，在山地生态系统中，地表地形对水循环具有明显的水文控制作用，但是传统山地植被GPP估算方法，无法应用在高分辨率下。
目前针对碳通量和水通量估算的空间降尺度算法通常假设在高空间分辨率中，低空间分辨率估算与地表异质性之间的关系是通用的（称为线性降尺度过程），实际上这种线性降尺度算法忽略了低分辨率估算的空间尺度误差，而采用低分辨率下得到的这些关系以获取高分辨率估算结果存在偏差，造成估算结果具有较高的不确定性。例如，专利文件CN110276304A中公开了一种基于降尺度的高分辨率植被生产力遥感估算方法，并具体公开了基于降尺度的高分辨率植被生产力遥感估算方法的实现包括以下步骤：(1)将植被生产力估算的因子LAI/FPAR（植被叶面积指数/光合有效辐射吸收比例）进行降尺度，得到时间序列的高分辨率的LAI/FPAR；(2)对温度进行空间插值和地形校正，得到高分辨率的温度因子；利用山地小气候模型模拟高分辨率的太阳短波辐射；将降尺度后的高分LAI/FPAR、经过地形校正后高分辨率的温度因子和太阳短波辐射数据输入植被生产力模型MuSyQ-NPP，得到连续时间序列的高分辨率植被生产力。
发明内容
针对山地环境高度时空异质性的地表水热条件和气候环境以及中低空间分辨率GPP估算中由于空间尺度误差造成的不确定性，本发明提供山地高空间分辨率植被总初级生产力估算方法，包括：
获取生态水文模型第一空间分辨率下的第一输入数据集；
对所述第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集；其中，第一空间分辨率为高空间分辨率；第二空间分辨率与第三空间分辨率为中低空间分辨率；
模拟植被光合作用过程，计算植被光合作用参数；
利用所述第二输入数据集、第三输入数据集与所述植被光合作用参数，采用生态水文模型进行植被总初级生产力估算，确定第二空间分辨率下的第一植被总初级生产力估算结果与第三空间分辨率下的第二植被总初级生产力估算结果；
根据所述第一输入数据集、所述第二输入数据集与所述第三输入数据集，获取第一空间分辨率下的第一地表异质性参数、第二空间分辨率下的第二地表异质性参数与第三空间分辨率下的第三地表异质性参数；
构建第二空间分辨率下的所述第一植被总初级生产力估算结果与所述第二地表异质性参数的第一关系模型；
将第三空间分辨率下的所述第三地表异质性参数作为所述第一关系模型的输入，得到第三空间分辨率下的第三植被总初级生产力估算结果；
计算第三空间分辨率下的所述第二植被总初级生产力估算结果与所述第三植被总初级生产力估算结果的差值；
在第二空间分辨率到第三空间分辨率范围内，构建所述差值与所述第二地表异质性参数、所述第三地表异质性参数的第二关系模型；
根据所述第一关系模型，将所述第一地表异质性参数作为所述第一关系模型的输入，估算第一空间分辨率下的第四植被总初级生产力估算结果；
将所述第一地表异质性参数、所述第二地表异质性参数作为所述第二关系模型的输入，计算第五植被总初级生产力估算结果；
计算所述第四植被总初级生产力估算结果与所述第五植被总初级生产力估算结果之和，得到目标植被总初级生产力估算结果。
本发明的有益效果是：本发明采用生态水文模型，替代常用于GPP数据集的LUE或VI模型，能够更好地表征山地生态系统的碳循环和水循环机制，通过光合作用过程对不同空间尺度环境变化的非线性响应过程来改进线性降尺度过程，减小了低分辨率估算空间尺度误差引起的不确定性。
在上述技术方案的基础上，本发明还可以做如下改进。
进一步，第一空间分辨率下的第一输入数据集包括站点气象数据、第一空间分辨率下的土地覆盖类型数据、MODIS与Landsat反射率数据、地形数据以及土壤类型数据。
采用上述进一步方案的有益效果是，
进一步，对所述第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集，包括：
获取第一空间分辨率下的第一时间序列LAI数据；
对所述第一时间序列LAI数据进行处理，确定第一空间分辨率下的第一土壤类型数据；
将第一空间分辨率下的第一输入数据集聚合为第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据：将第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据均由第一空间分辨率下的第一土壤类型数据组成；将第二空间分辨率下的第二土壤类型数据中的所有第一空间分辨率的第一时间序列LAI数据的均值作为第二空间分辨率下第一时间序列LAI数据的像元值，将第二空间分辨率下的第二土壤类型数据中所有第一空间分辨率的地形数据的均值作为第二空间分辨率下地形数据的像元值。进一步，获取第一空间分辨率下的第一时间序列LAI数据，包括：使用UofT LAI算法、STSG模型和IFSDAF时空重建模型从Landsat和MODIS反射率数据获得第一空间分辨率下的时间序列LAI数据。
进一步，模拟植被光合作用过程，计算植被光合作用参数，包括：
在光合作用过程中，将植被冠层叶片分为光照叶片和阴影叶片，分别模拟两种叶片的光合作用过程；
根据太阳天顶角、聚集指数和日长时间，确定冠层总叶片的叶面积指数、光照叶片的叶面积指数和阴影叶片的叶面积指数；
将Farquhar模型扩展到日尺度计算日尺度的光照叶片的光合作用速率和阴影叶片的光合作用速率。
进一步，所述第一地表异质性参数、所述第二地表异质性参数与所述第三地表异质性参数用于表征植被异质性指标与地形指标；
所述植被异质性指标根据土地覆盖类型数据、叶面积指数数据与气候指标进行表示；所述地形指标根据高程数据、坡度数据、坡向数据、土壤类型数据、地形湿度指数以及天空可视性因子进行表示；所述气候指标包括降水量数据、辐射数据、气温数据。
进一步，构建所述第一关系模型的方式为采用随机森林回归算法模拟第二空间分辨率下的所述第一植被总初级生产力估算结果与所述第二地表异质性参数的关系；在第二空间分辨率到第三空间分辨率范围内，构建所述第二关系模型的方式为采用随机森林回归算法模拟所述差值与所述第二地表异质性参数、所述第三地表异质性参数的关系。
附图说明
图1为本发明实施例1提供的山地高空间分辨率植被总初级生产力估算方法的流程图；
图2为本发明实施例1提供的山地高空间分辨率植被总初级生产力估算方法的原理图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
作为一个实施例，本发明实施例提供将生态水文模型与线性和非线性降尺度相结合的山地生态系统大尺度总初级生产力高分辨率估算方法，该方法可为区域或全球尺度的山地植被GPP遥感估算提供参考。如附图1所示，本实施例提供山地高空间分辨率植被总初级生产力估算方法，包括：
获取生态水文模型第一空间分辨率下的第一输入数据集；
对第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集；其中，第一空间分辨率为高空间分辨率；第二空间分辨率与第三空间分辨率为中低空间分辨率；
模拟植被光合作用过程，计算植被光合作用参数；
利用第二输入数据集、第三输入数据集与植被光合作用参数，采用生态水文模型进行植被总初级生产力估算，确定第二空间分辨率下的第一植被总初级生产力估算结果与第三空间分辨率下的第二植被总初级生产力估算结果；
根据第一输入数据集、第二输入数据集与第三输入数据集，获取第一空间分辨率下的第一地表异质性参数、第二空间分辨率下的第二地表异质性参数与第三空间分辨率下的第三地表异质性参数；
构建第二空间分辨率下的第一植被总初级生产力估算结果与第二地表异质性参数的第一关系模型；
将第三空间分辨率下的第三地表异质性参数作为第一关系模型的输入，得到第三空间分辨率下的第三植被总初级生产力估算结果；
计算第三空间分辨率下的第二植被总初级生产力估算结果与第三植被总初级生产力估算结果的差值；
在第二空间分辨率到第三空间分辨率范围内，构建差值与第二地表异质性参数、第三地表异质性参数的第二关系模型；
根据第一关系模型，将第一地表异质性参数作为第一关系模型的输入，估算第一空间分辨率下的第四植被总初级生产力估算结果；
将第一地表异质性参数、第二地表异质性参数作为第二关系模型的输入，计算第五植被总初级生产力估算结果；
计算第四植被总初级生产力估算结果与第五植被总初级生产力估算结果之和，得到目标植被总初级生产力估算结果。
可选的，第一空间分辨率下的第一输入数据集包括站点气象数据、第一空间分辨率下的土地覆盖类型数据、MODIS与Landsat反射率数据、地形数据以及土壤类型数据。
在实际应用过程中，生态水文模型能够更好地表征植被生理过程和水文过程之间存在复杂的交互作用，该生态水文模型输入数据包括：站点气象数据，涵盖逐日的太阳入射辐射、空气温度与降水数据；30m空间分辨率的土地覆盖类型数据；MODIS（moderate-resolution imaging spectroradiometer,中分辨率成像光谱仪）与Landsat反射率数据；30m空间分辨率的地形数据，涵盖高程、坡度、坡向数据以及土壤类型数据。
可选的，对第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集，包括：
获取第一空间分辨率下的第一时间序列LAI数据；
对所述第一时间序列LAI数据进行处理，确定第一空间分辨率下的第一土壤类型数据；
将第一空间分辨率下的第一输入数据集聚合为第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据：将第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据均由第一空间分辨率下的第一土壤类型数据组成；将第二空间分辨率下的第二土壤类型数据中的所有第一空间分辨率的第一时间序列LAI数据的均值作为第二空间分辨率下第一时间序列LAI数据的像元值，将第二空间分辨率下的第二土壤类型数据中所有第一空间分辨率的地形数据的均值作为第二空间分辨率下地形数据的像元值。
可选的，获取第一空间分辨率下的第一时间序列LAI数据，包括：使用UofT LAI算法、STSG模型和IFSDAF时空重建模型从Landsat和MODIS反射率数据获得第一空间分辨率下的时间序列LAI数据。
可选的，模拟植被光合作用过程，计算植被光合作用参数，包括：
在光合作用过程中，将植被冠层叶片分为光照叶片和阴影叶片，分别模拟两种叶片的光合作用过程；
根据太阳天顶角、聚集指数和日长时间，确定冠层总叶片的叶面积指数、光照叶片的叶面积指数和阴影叶片的叶面积指数；
将Farquhar模型扩展到日尺度计算日尺度的光照叶片的光合作用速率和阴影叶片的光合作用速率。
具体的，在BEPS-TerrainLab光合作用过程中，将植被冠层叶片分为光照叶片和阴
影叶片，分别模拟两种叶片的光合作用过程。设LAI为冠层总叶片的叶面积指数、LAI sun 为光
照叶片的叶面积指数和LAI shd 为阴影叶片的叶面积指数；Z s 为太阳天顶角，Ω为聚集指数，DAYL为日长时间，A sun 是光照叶片的光合作用速率，A shd 是阴影叶片的光合作用速率，
为植被总初级生产力估算结果，则：
其中，A sun,unsat 为不饱和水分状态下光照冠层的净光合作用速率，A sun,sat 为饱和水分状态下光照冠层的净光合作用速率，A shd,unsat 为不饱和水分状态下阴影冠层的净光合作用速率，A shd,sat 为饱和水分状态下阴影冠层的净光合作用速率；u是处于根系处于不饱和水分状态的比例，(1-u)则表示根系处于饱和水分状态的比例。上述四种净光合作用速率可基于Farquhar模型获得。设A为目标净光合作用速率，A c 为Rubiso酶限制下的光合作用速率，V cmax 为最大羧化反应速率，C i 为叶肉细胞的二氧化碳的浓度，O i 为叶肉细胞的氧气的浓度，K c 为二氧化碳的Michaelis-Menten常数，K o 为氧气的Michaelis-Menten常数，A j 为RuBP酶限制下的光合作用速率，J max 为最大电子传输速率，Γ为没有暗呼吸的二氧化碳补偿点，R d 为白天暗呼吸速率，则：
其中，R d 取值一般等于0.015V cmax ，则目标净光合作用速率：
可选的，第一地表异质性参数、第二地表异质性参数与第三地表异质性参数用于表征植被异质性指标与地形指标；
植被异质性指标根据土地覆盖类型数据、叶面积指数LAI数据与气候指标进行表示；地形指标根据高程数据、坡度数据、坡向数据、土壤类型数据、地形湿度指数（topographic wetness index,TWI）以及天空可视性因子（sky-view factor,SVF）进行表示；气候指标包括降水量数据、辐射数据、气温数据。
TWI与SVF可按下式进行计算：
其中，
为水流聚集面积，
为斜率，
为垂直高度的角度，
为表示方向
的总数。
可选的，构建第一关系模型的方式为采用随机森林回归算法模拟第二空间分辨率下的第一植被总初级生产力估算结果与第二地表异质性参数的关系；在第二空间分辨率到第三空间分辨率范围内，构建第二关系模型的方式为采用随机森林回归算法模拟差值与第二地表异质性参数、第三地表异质性参数的关系。
具体的，本发明实施例中，第一空间分辨率取30m分辨率作为高空间分辨率，第二空间分辨率取480m分辨率、第三空间分辨率取120m分辨率作为两个中低空间分辨率。对30m分辨率的数据集进行处理，聚合得到480m空间分辨率的数据集和120m空间分辨率的数据集。
如附图2所示，该山地高空间分辨率植被总初级生产力估算方法，包括：
（1）数据集获取与处理：
获取生态水文模型30m空间分辨率下的第一输入数据集；
对第一输入数据集进行处理，得到480m空间分辨率下的第二输入数据集与120m空间分辨率下的第三输入数据集；其中，30m空间分辨率为高空间分辨率；480m空间分辨率与120m空间分辨率为中低空间分辨率；
基于BEPS-TerrainLab模型，模拟植被光合作用过程，计算植被光合作用参数；
基于BEPS-TerrainLab生态水文模型（简称BTL）分别以480m空间分辨率和120m空
间分辨率估算山地植被GPP，得到480m空间分辨率的山地植被GPP：
，以及120m空间
分辨率的山地植被GPP：
；
利用第二输入数据集、第三输入数据集与植被光合作用参数，采用生态水文模型进行植被总初级生产力估算，确定480m空间分辨率下的第一植被总初级生产力估算结果与120m空间分辨率下的第二植被总初级生产力估算结果；
根据第一输入数据集、第二输入数据集与第三输入数据集，获取30m空间分辨率下的第一地表异质性参数、480m空间分辨率下的第二地表异质性参数与120m空间分辨率下的第三地表异质性参数；
（2）线性降尺度过程：
构建第二空间分辨率下的第一植被总初级生产力估算结果与第二地表异质性参数的第一关系模型；
在实际应用过程中，假设光合作用对环境变化的响应在空间尺度上是线性的，则
低分辨率上的GPP估算与地表异质性之间的关系可以用于高分辨率上的GPP估算，设
为480m分辨率下描述植被异质性和表面地形的指标，
为480m分辨率下描述表面地形
的指标，地形指标根据高程数据、坡度数据、坡向数据、土壤类型数据、地形湿度指数以及天
空可视性因子；设
为表示气候指标，气候指标包括降水量数据、辐射数据、气温数据；
第一关系模型用于表征采用随机森林回归算法模拟480m空间分辨率下GPP估算结果与地表
异质性参数的关系，用
表示；480m空间分辨率山地植被GPP与植被异质性参数之间的关
系表示为：
由于
为基于BEPS-TerrainLab生态水文模型得到的480m空间分辨率的山
地植被GPP；
（3）非线性降尺度过程：
将120m空间分辨率下的第三地表异质性参数作为第一关系模型的输入，估算120m
空间分辨率下的第三植被总初级生产力估算结果
；
设
为120m空间分辨率下描述植被异质性的指标，
为120m空间分辨率
下描述表面地形的指标；则120m空间分辨率下的第三植被总初级生产力估算结果为：
计算120m空间分辨率下的第二植被总初级生产力估算结果与第三植被总初级生产力估算结果的差值；
在实际过程中，由于植被光合作用对环境变化未表现为标准的线性响应，因此存
在空间尺度误差使得则
和
之间存在偏差，通过下式计算偏差，得到差值
：
在120m空间分辨率到480m空间分辨率范围内，采用随机森林回归算法构建差值
与第二地表异质性参数、第三地表异质性参数的第二关系模型，用
表示；
（4）采用线性和非线性降尺度过程相结合进行30m分辨率GPP估算：
将第一地表异质性参数作为第一关系模型的输入，估算30m空间分辨率下的第四
植被总初级生产力估算结果
；
设
为30m分辨率下描述植被异质性的指标，
为30m分辨率下描述表面
地形的指标，第四植被总初级生产力估算结果为：
在30m空间分辨率到480m空间分辨率范围内，将第一地表异质性参数、第二地表异
质性参数作为第二关系模型的输入，计算第五植被总初级生产力估算结果
：
第五植被总初级生产力估算结果计算第四植被总初级生产力估算结果
与第五植被总初级生产力估算结果
之和，得到30m空间分辨率下的目标植被总初
级生产力估算结果。
本发明采用生态水文模型，替代常用于GPP数据集的LUE或VI模型，能够更好地表征山地生态系统的碳循环和水循环机制，通过光合作用过程对不同空间尺度环境变化的非线性响应过程来改进线性降尺度过程，减小了低分辨率估算空间尺度误差引起的不确定性。
以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.山地高空间分辨率植被总初级生产力估算方法，其特征在于，包括：
获取生态水文模型第一空间分辨率下的第一输入数据集；
对所述第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集；其中，第一空间分辨率为高空间分辨率；第二空间分辨率与第三空间分辨率为中低空间分辨率；
模拟植被光合作用过程，计算植被光合作用参数；
利用所述第二输入数据集、第三输入数据集与所述植被光合作用参数，采用生态水文模型进行植被总初级生产力估算，确定第二空间分辨率下的第一植被总初级生产力估算结果与第三空间分辨率下的第二植被总初级生产力估算结果；
根据所述第一输入数据集、所述第二输入数据集与所述第三输入数据集，获取第一空间分辨率下的第一地表异质性参数、第二空间分辨率下的第二地表异质性参数与第三空间分辨率下的第三地表异质性参数；
构建第二空间分辨率下的所述第一植被总初级生产力估算结果与所述第二地表异质性参数的第一关系模型；
将第三空间分辨率下的所述第三地表异质性参数作为所述第一关系模型的输入，得到第三空间分辨率下的第三植被总初级生产力估算结果；
计算第三空间分辨率下的所述第二植被总初级生产力估算结果与所述第三植被总初级生产力估算结果的差值；
在第二空间分辨率到第三空间分辨率范围内，构建所述差值与所述第二地表异质性参数、所述第三地表异质性参数的第二关系模型；
根据所述第一关系模型，将所述第一地表异质性参数作为所述第一关系模型的输入，估算第一空间分辨率下的第四植被总初级生产力估算结果；
将所述第一地表异质性参数、所述第二地表异质性参数作为所述第二关系模型的输入，计算第五植被总初级生产力估算结果；
计算所述第四植被总初级生产力估算结果与所述第五植被总初级生产力估算结果之和，得到目标植被总初级生产力估算结果。
2.根据权利要求1所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，第一空间分辨率下的第一输入数据集包括站点气象数据、第一空间分辨率下的土地覆盖类型数据、MODIS与Landsat反射率数据、地形数据以及土壤类型数据。
3.根据权利要求1所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，对所述第一输入数据集进行处理，得到第二空间分辨率下的第二输入数据集与第三空间分辨率下的第三输入数据集，包括：
获取第一空间分辨率下的第一时间序列LAI数据；
对所述第一时间序列LAI数据进行处理，确定第一空间分辨率下的第一土壤类型数据；
将第一空间分辨率下的第一输入数据集聚合为第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据：将第二空间分辨率下的第二土壤类型数据与第三空间分辨率下的第三土壤类型数据均由第一空间分辨率下的第一土壤类型数据组成；将第二空间分辨率下的第二土壤类型数据中的所有第一空间分辨率的第一时间序列LAI数据的均值作为第二空间分辨率下第一时间序列LAI数据的像元值，将第二空间分辨率下的第二土壤类型数据中所有第一空间分辨率的地形数据的均值作为第二空间分辨率下地形数据的像元值。
4.根据权利要求3所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，获取第一空间分辨率下的第一时间序列LAI数据，包括：使用UofT LAI算法、STSG模型和IFSDAF时空重建模型从Landsat和MODIS反射率数据获得第一空间分辨率下的时间序列LAI数据。
5.根据权利要求1所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，模拟植被光合作用过程，计算植被光合作用参数，包括：
在光合作用过程中，将植被冠层叶片分为光照叶片和阴影叶片，分别模拟两种叶片的光合作用过程；
根据太阳天顶角、聚集指数和日长时间，确定冠层总叶片的叶面积指数、光照叶片的叶面积指数和阴影叶片的叶面积指数；
将Farquhar模型扩展到日尺度计算日尺度的光照叶片的光合作用速率和阴影叶片的光合作用速率。
6.根据权利要求1所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，所述第一地表异质性参数、所述第二地表异质性参数与所述第三地表异质性参数用于表征植被异质性指标与地形指标；
所述植被异质性指标根据土地覆盖类型数据、叶面积指数数据与气候指标进行表示；所述地形指标根据高程数据、坡度数据、坡向数据、土壤类型数据、地形湿度指数以及天空可视性因子进行表示；所述气候指标包括降水量数据、辐射数据、气温数据。
7.根据权利要求1所述山地高空间分辨率植被总初级生产力估算方法，其特征在于，构建所述第一关系模型的方式为采用随机森林回归算法模拟第二空间分辨率下的所述第一植被总初级生产力估算结果与所述第二地表异质性参数的关系；在第二空间分辨率到第三空间分辨率范围内，构建所述第二关系模型的方式为采用随机森林回归算法模拟所述差值与所述第二地表异质性参数、所述第三地表异质性参数的关系。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211278846.2A CN115358095B (zh) | 2022-10-19 | 2022-10-19 | 山地高空间分辨率植被总初级生产力估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211278846.2A CN115358095B (zh) | 2022-10-19 | 2022-10-19 | 山地高空间分辨率植被总初级生产力估算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115358095A true CN115358095A (zh) | 2022-11-18 |
| CN115358095B CN115358095B (zh) | 2023-03-17 |
Family
ID=84008549
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211278846.2A Active CN115358095B (zh) | 2022-10-19 | 2022-10-19 | 山地高空间分辨率植被总初级生产力估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115358095B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116227246A (zh) * | 2023-05-09 | 2023-06-06 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种植被总初级生产力遥感产品的地形校正方法 |
| CN116485201A (zh) * | 2023-04-24 | 2023-07-25 | 中国科学院大气物理研究所 | 一种基于par能量平衡的净生态系统生产力评估方法 |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110101239A1 (en) * | 2008-05-08 | 2011-05-05 | Iain Woodhouse | Remote sensing system |
| CN106446564A (zh) * | 2016-09-29 | 2017-02-22 | 华南农业大学 | 一种植被净初级生产力遥感估算方法 |
| CN110276304A (zh) * | 2019-06-25 | 2019-09-24 | 北京师范大学 | 基于降尺度的高分辨率植被生产力遥感估算方法 |
| CN110443504A (zh) * | 2019-07-30 | 2019-11-12 | 中国环境科学研究院 | 一种植被净初级生产力的遥感估算方法 |
| CN111242022A (zh) * | 2020-01-10 | 2020-06-05 | 西安科技大学 | 基于低分辨率遥感产品降尺度的高分辨率fapar估算方法 |
| CN111582703A (zh) * | 2020-04-30 | 2020-08-25 | 中国科学院空天信息创新研究院 | 全球中高分辨率植被总初级生产力产品的计算方法及装置 |
| CN112149295A (zh) * | 2020-09-17 | 2020-12-29 | 中国科学院空天信息创新研究院 | 一种全球通用植被总初级生产力遥感指数估算方法 |
| CN113139901A (zh) * | 2021-04-15 | 2021-07-20 | 青岛地质工程勘察院(青岛地质勘查开发局) | 一种流域尺度植被净初级生产力遥感精细反演方法 |
| CN113516362A (zh) * | 2021-05-17 | 2021-10-19 | 深圳大学 | 基于遥感与过程耦合模型的红树林净初级生产力估算方法 |
| CN113888386A (zh) * | 2021-12-03 | 2022-01-04 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种高分辨率的时空无缝地表土壤水分估算方法及系统 |
| WO2022032329A1 (en) * | 2020-08-14 | 2022-02-17 | Agriculture Victoria Services Pty Ltd | System and method for image-based remote sensing of crop plants |
| US20220061236A1 (en) * | 2020-08-25 | 2022-03-03 | The Board Of Trustees Of The University Of Illinois | Accessing agriculture productivity and sustainability |
| CN114331233A (zh) * | 2022-03-15 | 2022-04-12 | 航天宏图信息技术股份有限公司 | 植被总初级生产力估算方法、装置、电子设备及存储介质 |
| CN114511550A (zh) * | 2022-02-22 | 2022-05-17 | 江西财经大学 | 一种鄱阳湖湿地植被净初级生产力的遥感估算方法 |
-
2022
- 2022-10-19 CN CN202211278846.2A patent/CN115358095B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110101239A1 (en) * | 2008-05-08 | 2011-05-05 | Iain Woodhouse | Remote sensing system |
| CN106446564A (zh) * | 2016-09-29 | 2017-02-22 | 华南农业大学 | 一种植被净初级生产力遥感估算方法 |
| CN110276304A (zh) * | 2019-06-25 | 2019-09-24 | 北京师范大学 | 基于降尺度的高分辨率植被生产力遥感估算方法 |
| CN110443504A (zh) * | 2019-07-30 | 2019-11-12 | 中国环境科学研究院 | 一种植被净初级生产力的遥感估算方法 |
| CN111242022A (zh) * | 2020-01-10 | 2020-06-05 | 西安科技大学 | 基于低分辨率遥感产品降尺度的高分辨率fapar估算方法 |
| CN111582703A (zh) * | 2020-04-30 | 2020-08-25 | 中国科学院空天信息创新研究院 | 全球中高分辨率植被总初级生产力产品的计算方法及装置 |
| WO2022032329A1 (en) * | 2020-08-14 | 2022-02-17 | Agriculture Victoria Services Pty Ltd | System and method for image-based remote sensing of crop plants |
| US20220061236A1 (en) * | 2020-08-25 | 2022-03-03 | The Board Of Trustees Of The University Of Illinois | Accessing agriculture productivity and sustainability |
| CN112149295A (zh) * | 2020-09-17 | 2020-12-29 | 中国科学院空天信息创新研究院 | 一种全球通用植被总初级生产力遥感指数估算方法 |
| CN113139901A (zh) * | 2021-04-15 | 2021-07-20 | 青岛地质工程勘察院(青岛地质勘查开发局) | 一种流域尺度植被净初级生产力遥感精细反演方法 |
| CN113516362A (zh) * | 2021-05-17 | 2021-10-19 | 深圳大学 | 基于遥感与过程耦合模型的红树林净初级生产力估算方法 |
| CN113888386A (zh) * | 2021-12-03 | 2022-01-04 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种高分辨率的时空无缝地表土壤水分估算方法及系统 |
| CN114511550A (zh) * | 2022-02-22 | 2022-05-17 | 江西财经大学 | 一种鄱阳湖湿地植被净初级生产力的遥感估算方法 |
| CN114331233A (zh) * | 2022-03-15 | 2022-04-12 | 航天宏图信息技术股份有限公司 | 植被总初级生产力估算方法、装置、电子设备及存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 谢馨瑶: "基于邻域空间的山地植被GPP遥感估算建模及其空间尺度扩展研究", 《中国博士学位论文全文数据库 (基础科学辑)》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116485201A (zh) * | 2023-04-24 | 2023-07-25 | 中国科学院大气物理研究所 | 一种基于par能量平衡的净生态系统生产力评估方法 |
| CN116227246A (zh) * | 2023-05-09 | 2023-06-06 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种植被总初级生产力遥感产品的地形校正方法 |
| CN116227246B (zh) * | 2023-05-09 | 2023-08-11 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种植被总初级生产力遥感产品的地形校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115358095B (zh) | 2023-03-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115358095B (zh) | 山地高空间分辨率植被总初级生产力估算方法 | |
| Xu et al. | Assessing the impact of urbanization on regional net primary productivity in Jiangyin County, China | |
| CN110276304A (zh) | 基于降尺度的高分辨率植被生产力遥感估算方法 | |
| Ji et al. | Mapping land use/cover dynamics of the Yellow River Basin from 1986 to 2018 supported by Google Earth Engine | |
| Maeda et al. | Prospective changes in irrigation water requirements caused by agricultural expansion and climate changes in the eastern arc mountains of Kenya | |
| CN113139901A (zh) | 一种流域尺度植被净初级生产力遥感精细反演方法 | |
| Ju et al. | Prediction of summer grain crop yield with a process-based ecosystem model and remote sensing data for the northern area of the Jiangsu Province, China | |
| Wang et al. | Yield estimation of winter wheat in the North China Plain using the remote-sensing–photosynthesis–yield estimation for crops (RS–P–YEC) model | |
| CN112857583B (zh) | 一种基于卫星遥感数据的植被温度和裸地温度估算方法 | |
| CN114510824A (zh) | 同步考虑植被冠层和根系动态变化的区域尺度蒸散发模型的构建方法 | |
| Dudiak et al. | Soil and climatic bonitation of agricultural lands of the steppe zone of Ukraine | |
| CN114417646B (zh) | 一种高维异构降水数据融合方法及系统 | |
| CN110927120A (zh) | 一种植被覆盖度预警方法 | |
| CN115796659B (zh) | 一种基于区域尺度蒸发衰退估计城市蓄水能力的方法 | |
| CN113610421A (zh) | 一种基于地形和土地利用数据的流域生态资产价值评估方法 | |
| CN115239175B (zh) | 一种陆地生态系统生产总值的快速估算方法及其系统 | |
| Qin et al. | Reconstruction of 60-year (1961–2020) surface air temperature on the Tibetan Plateau by fusing MODIS and ERA5 temperatures | |
| Xian | Analysis of impacts of urban land use and land cover on air quality in the Las Vegas region using remote sensing information and ground observations | |
| Lee et al. | Combining flux variance similarity partitioning with artificial neural networks to gap-fill measurements of net ecosystem production of a Pacific Northwest Douglas-fir stand | |
| CN117496369A (zh) | 一种矿区碳汇数据集构建方法 | |
| Xie et al. | TAVIs: Topographically adjusted vegetation index for a reliable proxy of gross primary productivity in mountain ecosystems | |
| Ellenburg et al. | Evaluation of a regional crop model implementation for sub-national yield assessments in Kenya | |
| CN111078811A (zh) | 一种净初级生产力预警方法 | |
| Hazimeh et al. | Impact of ET and biomass model choices on economic irrigation water productivity in water-scarce basins | |
| CN116227246B (zh) | 一种植被总初级生产力遥感产品的地形校正方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |