CN109884664B - 一种城市地上生物量光学微波协同反演方法及系统 - Google Patents

一种城市地上生物量光学微波协同反演方法及系统 Download PDF

Info

Publication number
CN109884664B
CN109884664B CN201910046622.0A CN201910046622A CN109884664B CN 109884664 B CN109884664 B CN 109884664B CN 201910046622 A CN201910046622 A CN 201910046622A CN 109884664 B CN109884664 B CN 109884664B
Authority
CN
China
Prior art keywords
microwave
optical
biomass
data
variable
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.)
Active
Application number
CN201910046622.0A
Other languages
English (en)
Other versions
CN109884664A (zh
Inventor
邵振峰
张邻晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201910046622.0A priority Critical patent/CN109884664B/zh
Publication of CN109884664A publication Critical patent/CN109884664A/zh
Priority to US16/741,778 priority patent/US11454534B2/en
Application granted granted Critical
Publication of CN109884664B publication Critical patent/CN109884664B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01G—WEIGHING
    • G01G9/00—Methods of, or apparatus for, the determination of weight, not provided for in groups G01G1/00 - G01G7/00
    • G01G9/005—Methods of, or apparatus for, the determination of weight, not provided for in groups G01G1/00 - G01G7/00 using radiations, e.g. radioactive
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
    • G01S13/865—Combination of radar systems with lidar systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
    • 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
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/25—Fusion techniques
    • G06F18/253—Fusion techniques of extracted features
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/90—Determination of colour characteristics
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
    • G06V10/80—Fusion, i.e. combining data from various sources at the sensor level, preprocessing level, feature extraction level or classification level
    • G06V10/806—Fusion, i.e. combining data from various sources at the sensor level, preprocessing level, feature extraction level or classification level of extracted features
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Abstract

一种城市地上生物量光学微波协同反演方法及系统，包括利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；进行预处理，得到冠层高度值CHM、地表反射率数据与后向散射系数；基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；提取LiDAR数据覆盖区域的生物量估测值；将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；采用变量筛选方法筛选出最佳光学与微波特征变量；构建AGB反演的光学模型、微波模型及光学微博协同的模型，选择最优的模型实现生物量反演。

Description

一种城市地上生物量光学微波协同反演方法及系统
技术领域
本发明涉及一种反演方法，具体来说，涉及一种城市地上生物量光学微波协同反演方法及系统。
背景技术
植被地上生物量对于生态系统固碳能力和碳储量有重要的指示作用，在国际地圈-生物圈研究计划(IGBP)中，碳循环被确定为全球变化和陆地生态系统(GCTE)等计划的重要研究内容。因此，植被地上生物量的准确提取对于全球变化及陆地生态系统监测具有重要意义。
遥感技术由于其宏观动态实时多源的特点，已在植被地上生物量以及其他植被参数研究中蓬勃开展。常见的方法是通过建立样地实测数据与遥感特征变量之间的关系模型来估测植被参数。但是，样地实测数据的获取面临众多挑战，一方面，用于建模及验证的样地数据的获取精度难以得到评估；另一方面，收集的样地数据难以覆盖整个研究区，无法充分反映研究区内植被结构及生长状态的变化；并且，由于时间、人力成本以及无法到达的区域等的限制，获取足够数量的样本来进行模型构建几乎是不可能的。近几年，越来越多的国内外研究表明，基于机载激光雷达(LiDAR)数据获取的植被参数，如冠层高度、冠层覆盖度以及地上生物量具有可与样地实测数据相媲美的精度。然而，相比卫星遥感数据，机载LiDAR数据在多时相应用以及大区域制图方面具有较大局限性。因此，结合机载LiDAR数据与卫星遥感数据进行区域尺度的植被参数估测具有高度的应用潜力，将机载LiDAR数据提取的植被参数作为训练和验证样本集，卫星遥感数据提供预测变量，是大区域植被参数反演的一个重要发展趋势。
在众多卫星遥感数据源中，多光谱卫星遥感数据，尤其是多光谱植被指数，仍然是植被地上生物量估测使用最为广泛的数据源之一。然而，被动光学遥感在获取生物量参数信息时，不但要受到云、雨、雪等天气现象的影响，而且当植被比较茂密时，光学遥感获得的植被反射波谱信息出现饱和。合成孔径雷达(SAR)卫星遥感数据具有无昼夜性、受天气影响小等优势，并且具有一定的穿透性，对植被结构信息较为敏感。许多研究表明，SAR后向散射系数与植被结构参数(如胸径、体积、树高、生物量等)存在显著的相关性。但是，利用SAR后向散射系数进行生物量估测仍然面临严峻的挑战，一方面，当植被覆盖度比较低时，微波遥感受地表土壤含水量、粗糙度等影响显著；另一方面，在浓密植被区域，SAR后向散射对植被信息的响应也会出现饱和。因此，针对目前生物量精准定量反演的迫切需求以及当前单源遥感数据存在的局限性，本领域亟待提供更有效的技术方案。
发明内容
对此，本发明提出了一种城市地上生物量光学微波协同反演方法及系统，有效结合光学与微波遥感提取生物量。。
本发明的技术方案提供一种城市地上生物量光学微波协同反演方法，包括以下步骤：
步骤a，利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；
步骤b，分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，相应得到冠层高度值CHM、地表反射率数据与后向散射系数；
步骤c，基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；
步骤d，以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
步骤e，将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；
步骤f，采用变量筛选方法筛选出最佳光学与微波特征变量；
步骤g，将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及光学微波协同的模型，根据步骤e选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
而且，所述步骤c中，选取的LiDAR变量包括最小值、最大值、平均值、高度分位数、标准差、变异系数、倾斜度、峰值以及冠层覆盖度。
而且，所述步骤c中，选取的光学特征植被指数包括归一化植被指数NDVI、比值植被指数SR、增强型植被指数EVI、土壤调节植被指数SAVI、修正土壤调节植被指数MSAVI，优化型土壤调节植被指数OSAVI、水分胁迫指数MSI、归一化水指数NDWI和叶绿素指数Clgreen，其计算公式如下：
Figure BDA0001949400640000031
Figure BDA0001949400640000032
Figure BDA0001949400640000033
Figure BDA0001949400640000034
Figure BDA0001949400640000035
Figure BDA0001949400640000036
Figure BDA0001949400640000037
Figure BDA0001949400640000038
Figure BDA0001949400640000039
式中，R、G、B、NIR和SWIR1分别为红光、绿光、蓝光、近红外及短波红外波段的反射率。
而且，所述步骤c中，选取的微波特征变量包括VV,HH,VH,HV,VV/HH,HH/HV,VV/HV以及RVI，其中RVI的计算公式如下：
Figure BDA00019494006400000310
其中，VV、HH、VH是极化方式，VV/HH、HH/HV、VV/HV是相应的比值。
而且，所述步骤f中，采用的变量选择方法为逐步筛选法。
而且，所述步骤f中，筛选出的最佳光学特征变量为SR,NDVI,OSAVI,MSI和NDWI，最佳微波特征变量为VV,HV,HH/HV和RVI。
而且，所述步骤g中，选用的预测模型包括多元逐步线性回归SLR、K最近邻KNN,支持向量机SVM,BP神经网络BPNN,随机森林RF以及深度学习DL，其中深度学习模型选用的是栈式稀疏自编码网络模型SSAE。
本发明提供一种城市地上生物量光学微波协同反演系统，包括以下模块：
第一模块，用于利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；
第二模块，用于分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，分别得到冠层高度值CHM、地表反射率数据与后向散射系数；
第三模块，用于基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；
第四模块，用于以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
第五模块，用于将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；
第六模块，用于采用变量筛选方法筛选出最佳光学与微波特征变量；
第七模块，用于将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及光学微波协同的模型，根据第五模块选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
综上所述，本发明具有如下特点：
(1)将基于LiDAR数据提取的生物量值作为训练及验证样本集，避免了使用费时费力的人工采样方法，对于充足且具代表性的野外调查数据难以获取的大面积区域具有尤为重要的意义。
(2)采用分层随机采样方法选择生物量样本，该数据集实现了对研究区生物量值在统计和地理位置上的全覆盖，充分反映了研究区植被群落的生长状况变化。
(3)构建了光学与微波数据协同反演模型，充分发挥两者估测生物量的优势，有效提高了地上生物量的定量反演精度。
附图说明
图1为本发明实施例的流程图。
具体实施方式
下面结合附图和实施例对本发明技术方案做进一步说明。
不同于以往的反演方法，本发明首先以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的地上生物量值；然后将LiDAR数据覆盖区域的生物量值作为生物量训练及验证样本集，利用分层随机采样的方法，选取一定数量的样本进行后续的建模及验证；最后采用变量筛选方法筛选出最佳光学与微波特征变量，并将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及二者协同的模型。
本发明的实施例是基于Landsat8Operational Land Imager(OLI)和RADARSAT-2，对广州市从化研究区森林地上生物量进行估算，参照图1，本发明实施例的具体步骤如下：
步骤a：研究区样地地上生物量计算，包括利用地面观测试验获取的样地单株参数数据集，基于地上生物量计算公式，获得样地生物量(AGB)观测值；
实施例中，通过地面观测获取了各样地内单木的树高和胸径参数的测量结果数据集。根据研究区主要树种的异速生长模型(J.Fang,G.Liu,and S.Xu,“Biomass and netproduction of forest vegetation in China,”ActaEcologicaSinica,vol.16,pp.497–508,1996.)，首先计算出样地内单木生物量，进而计算出样地总生物量，样地总生物量与样地面积之比即为样地生物量观测值，单位为Mg/ha。
步骤b：分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，分别得到冠层高度值(CHM)、地表反射率数据与后向散射系数；
实施例中进行LiDAR数据、Landsat8OLI以及RADARSAT-2数据预处理：
具体实施时可以利用Terrascan软件，首先，将高度大于其周围点高度中位数的点以及孤立点进行去除；其次，将剩余的LiDAR点云数据分为地面点和非地面点，再分别对地面点和所有首次回波插值生成数字高程模型(Digital Elevation Model，DEM)和数字表面模型(Digital Surface Model，DSM)，且分辨率均为1m；然后，将DSM减去DEM，得到1m分辨率的冠层高度模型(Crown Height Model，CHM)，即归一化的高度值，根据实地调查情况，保留值在2m-35m范围内的CHM像元，将其作为植被像元，以此排除研究区内的林下低矮植被和高于树木的地物；最后，以正射校正后的Landsat8影像为基准，对CHM影像进行几何校正。
具体实施时可以利用ENVI 5.1软件，首先，将数字信号值(Digital Number，DN)转换为表观辐亮度值，实现对影像的辐射定标；其次，采用FLAASH(Fast Line-of-sightAtmospheric Analysis of Spectral Hypercube)模型对定标后的影像进行大气校正，得到地表反射率影像；然后，利用全球公开的30m公开的DEM对大气校正后的影像进行正射校正；最后，基于预处理后的LiDAR CHM影像对正射校正后的Landsat8-OLI影像进行裁剪，得到覆盖从化研究区的光学影像。
具体实施时可以利用NEST(Next European Space Agency(ESA)SAR Toolbox)软件，首先对RADARSAT-2影像进行辐射定标得到后向散射系数；然后对SAR数据进行滤波来减小斑点噪声；利用预处理后的Landsat8OLI影像对SAR数据进行几何校正，并将其重采样到与Landsat8OLI一致的空间分辨率(30m)；最后将预处理后的Landsat8OLI和RADARSAT-2进行裁剪，只保留与研究区域一致的Landsat8OLI和RADARSAT-2影像。
步骤c：LiDAR变量、光学变量与微波变量提取：基于CHM数据，提取多种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量
基于CHM，计算的LiDAR变量包括最小值、最大值、平均值、高度分位数、标准差、变异系数、倾斜度、峰值以及冠层覆盖度。这些是LiDAR统计变量中常用的变量，本发明不予赘述。例如提取限定范围内的CHM值最小值。取值范围是与卫星影像相同的分辨率，即30m*30m内获取的这些lidar变量。
基于预处理后的Landsat8OLI地表反射率数据，计算归一化植被指数NDVI、比值植被指数SR、增强型植被指数EVI、土壤调节植被指数SAVI、修正土壤调节植被指数MSAVI，优化型土壤调节植被指数OSAVI、水分胁迫指数MSI、归一化水指数NDWI和叶绿素指数Clgreen，其数学表达式见公式(1)-(9)。
Figure BDA0001949400640000061
Figure BDA0001949400640000062
Figure BDA0001949400640000063
Figure BDA0001949400640000064
Figure BDA0001949400640000065
Figure BDA0001949400640000066
Figure BDA0001949400640000067
Figure BDA0001949400640000068
Figure BDA0001949400640000071
式中，R、G、B、NIR和SWIR1分别为红光、绿光、蓝光、近红外及短波红外波段的反射率。
基于预处理后的RADARSAT-2数据，计算的微波特征变量包括VV,HH,VH,HV,VV/HH,HH/HV,VV/HV以及RVI，这里的VV、HH、VH是RADARSAT-2数据的三种极化方式，VV/HH、HH/HV、VV/HV是相应的比值，其中RVI的计算公式如下：
Figure BDA0001949400640000072
步骤d：AGB基准图获取，包括以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
实施例中，以样地实测AGB为因变量，提取的LiDAR变量为自变量，采用多元逐步线性回归方法建立AGB与LiDAR变量之间的关系模型，得到AGB-LiDAR变量回归模型，进而得到LiDAR数据覆盖区域的生物量值，作为后续反演模型建立的生物量基准图。
步骤e：分层随机采样，包括将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取一定数量的样本进行后续的建模及验证；
将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，采用分层随机采样的方法，选取一定数量的样本进行后续的建模及验证。在实施分层随机采样时，选取平均值和标准差两个LiDAR高度变量作为先验知识来进行分层，这是因为这两个变量能够直接显示森林结构及其生长状态。首先，进行LiDAR条带覆盖区域的数据准备，包括生物量基准图中的AGB值及其对应的预测变量值(包括所有光学变量、SAR变量以及高度平均值和标准差)；其次，将以上数据集按照平均值进行升序排列，再平均分为10个大小相同的数据层；然后，将每个数据层按照标准差进行升序排列，再平均分为4等份，从而一共生成40个相同大小的数据层；最后，从每个数据层中随机选择20个样本，总共可获得800个样本数据，并按照3:1的比例随机分为训练和验证样本集。
步骤f：最佳光学与微波变量选择，包括采用变量筛选方法筛选出最佳光学与微波特征变量；
实施例优选利用逐步筛选法分别对光学与微波特征变量进行筛选，具体步骤为：(1)设n为预测变量的个数，每个变量均有一次排除在预测变量集之外，形成n组变量，每组包含n-1个变量，将n组变量分别输入预测模型中进行模型训练；(2)对(1)中所有模型进行精度验证，验证指标为RMSE和R2，具有最佳预测效果(即R2最大)的模型能够间接证明，该模型未采用的预测变量具有最差的解释能力；(3)将步骤(2)判别出的解释力最差的变量剔除出候选变量集；
(4)重复以上步骤，直到R2收敛为止。需要注意的是，模型的精度评价采用的是所有预测模型(即SLR、KNN、BPNN、SVM、RF和SSAE)验证精度的平均值。最终筛选出的最佳光学特征变量为SR,NDVI,OSAVI,MSI和NDWI，最佳微波特征变量为VV,HV,HH/HV和RVI。
步骤g：光学、微波及二者协同的生物量反演模型构建，包括将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及二者协同的模型，即生物量光学反演模型、生物量微波反演模型及生物量光学微波协同反演模型。根据步骤e选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
优选地，选用的预测模型包括多元逐步线性回归(SLR)、K最近邻(KNN),支持向量机(SVM),BP神经网络(BPNN),随机森林(RF)以及深度学习(DL)，其中深度学习模型选用的是栈式稀疏自编码网络模型(Stacked Sparse Autoencoder network，SSAE)。
相比传统的统计回归方法，K最近邻(KNN)、支持向量机(SVM)、误差反向传播神经网络(BPNN)、随机森林(RF)方法以及近几年关注度颇高的深度学习方法(DL)能够有效描述植被参数与预测变量之间复杂的非线性关系，在生物量以及其他植被参数的估测研究中展现出了更大的应用潜力，比参数化方法更易于工程化实施。
实施例中，将最佳光学变量(SR、NDVI、OSAVI、MSI、NDWI)、最佳微波变量(VV、HV、HH/HV、RVI)及二者的协同分别应用于6种预测模型(SLR、KNN、SVM、BPNN、RF、SSAE)中，协同策略是将所有的最佳光学变量和所有的最佳微波变量一起分别输入这6种预测模型，构建AGB的光学反演模型、微波反演模型及二者协同反演模型，其中光学反演模型的输入变量为最佳光学变量；微波反演模型的输入变量为最佳微波变量；协同反演模型的输入变量为所有的最佳光学变量和所有的最佳微波变量。
实验基于本发明提出的光学微波协同反演模型，利用Landsat8OLI和RADARSAT-2卫星数据，估算广州市从化研究区地上生物量的结果图。
以广州市从化研究区森林为实验对象，将光学反演模型和微波反演模型估测结果与协同反演结果进行对比验证，客观评价指标选取判定系数(R2)、均方根误差(RMSE)和相对均方根误差(RMSEr)。
表1生物量反演模型精度评价结果
Figure BDA0001949400640000091
验证结果如表1所示，发现在所有预测方法中，深度学习(SSAE)方法具有最高的估测精度，并且光学微波数据协同在所有预测方法中均好于单一数据源的预测结果，因此，光学微波数据协同并结合深度学习模型取得了最好的估测效果，判定系数为0.812，均方根误差为21.753，相对均方根误差为14.457％。表明光学与微波数据协同能够发挥两者在AGB反演中各自的优势，提高反演精度。
具体实施时，以上技术方案可采用计算机软件技术实现自动运行流程，也可以采用模块化方式提供相应系统。本发明实施例还提供一种城市地上生物量光学微波协同反演系统，包括以下模块：
第一模块，用于利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；
第二模块，用于分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，分别得到冠层高度值CHM、地表反射率数据与后向散射系数；
第三模块，用于基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；
第四模块，用于以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
第五模块，用于将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；
第六模块，用于采用变量筛选方法筛选出最佳光学与微波特征变量；
第七模块，用于将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及光学微波协同的模型，根据第五模块选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
具体各模块实现可参见相应步骤，本发明不予赘述。
本文中所描述的具体实施例仅仅是对本发明精神做举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或超越所附权利要求书所定义的范围。

Claims (8)

1.一种城市地上生物量光学微波协同反演方法，其特征在于，包括以下步骤：
步骤a，利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；
步骤b，分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，相应得到冠层高度值CHM、地表反射率数据与后向散射系数；
步骤c，基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；
步骤d，以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
步骤e，将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；
步骤f，采用变量筛选方法筛选出最佳光学与微波特征变量；
步骤g，将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及光学微波协同的模型，根据步骤e选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
2.根据权利要求1所述的城市地上生物量光学微波协同反演方法，其特征在于：所述步骤c中，选取的LiDAR变量包括最小值、最大值、平均值、高度分位数、标准差、变异系数、倾斜度、峰值以及冠层覆盖度。
3.根据权利要求1所述的城市地上生物量光学微波协同反演方法，其特征在于：所述步骤c中，选取的光学特征植被指数包括归一化植被指数NDVI、比值植被指数SR、增强型植被指数EVI、土壤调节植被指数SAVI、修正土壤调节植被指数MSAVI，优化型土壤调节植被指数OSAVI、水分胁迫指数MSI、归一化水指数NDWI和叶绿素指数Clgreen，其计算公式如下：
Figure FDA0001949400630000011
Figure FDA0001949400630000012
Figure FDA0001949400630000013
Figure FDA0001949400630000021
Figure FDA0001949400630000022
Figure FDA0001949400630000023
Figure FDA0001949400630000024
Figure FDA0001949400630000025
Figure FDA0001949400630000026
式中，R、G、B、NIR和SWIR1分别为红光、绿光、蓝光、近红外及短波红外波段的反射率。
4.根据权利要求1所述的城市地上生物量光学微波协同反演方法，其特征在于：所述步骤c中，选取的微波特征变量包括VV,HH,VH,HV,VV/HH,HH/HV,VV/HV以及RVI，其中RVI的计算公式如下：
Figure FDA0001949400630000027
其中，VV、HH、VH是极化方式，VV/HH、HH/HV、VV/HV是相应的比值。
5.根据权利要求1所述的城市地上生物量光学微波协同反演方法，其特征在于：所述步骤f中，采用的变量选择方法为逐步筛选法。
6.根据权利要求5所述的城市地上生物量光学微波协同反演方法，其特征在于：所述步骤f中，筛选出的最佳光学特征变量为SR,NDVI,OSAVI,MSI和NDWI，最佳微波特征变量为VV,HV,HH/HV和RVI。
7.根据权利要求1所述的城市地 上生物量光学微波协同反演方法，其特征在于：所述步骤g中，选用的预测模型包括多元逐步线性回归SLR、K最近邻KNN,支持向量机SVM,BP神经网络BPNN,随机森林RF以及深度学习DL，其中深度学习模型选用的是栈式稀疏自编码网络模型SSAE。
8.一种城市地上生物量光学微波协同反演系统，其特征在于，包括以下模块：
第一模块，用于利用地面观测试验获取的样地单株参数数据集，获得样地生物量AGB观测值；
第二模块，用于分别对覆盖研究区的LiDAR数据、光学遥感数据和微波遥感数据进行预处理，分别得到冠层高度值CHM、地表反射率数据与后向散射系数；
第三模块，用于基于CHM数据，提取各种LiDAR变量，并基于光学地表反射率数据，提取多种光学特征植被指数，同时基于微波后向散射系数数据，提取多种微波特征变量；
第四模块，用于以样地实测AGB为因变量，提取的LiDAR变量为自变量，建立生物量的多元逐步线性回归模型，得到LiDAR数据覆盖区域的生物量估测值；
第五模块，用于将LiDAR数据覆盖区域的生物量值作为训练及验证样本集，利用分层随机采样的方法，选取样本用于后续的建模及验证；
第六模块，用于采用变量筛选方法筛选出最佳光学与微波特征变量；
第七模块，用于将最佳光学变量、最佳微波变量及二者的结合分别应用于多种预测模型中，分别构建AGB反演的光学模型、微波模型及光学微波协同的模型，根据第五模块选取的样本实现建模及验证后，选择最优的模型实现生物量反演。
CN201910046622.0A 2019-01-14 2019-01-14 一种城市地上生物量光学微波协同反演方法及系统 Active CN109884664B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910046622.0A CN109884664B (zh) 2019-01-14 2019-01-14 一种城市地上生物量光学微波协同反演方法及系统
US16/741,778 US11454534B2 (en) 2019-01-14 2020-01-14 Method and system for optical and microwave synergistic retrieval of aboveground biomass

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910046622.0A CN109884664B (zh) 2019-01-14 2019-01-14 一种城市地上生物量光学微波协同反演方法及系统

Publications (2)

Publication Number Publication Date
CN109884664A CN109884664A (zh) 2019-06-14
CN109884664B true CN109884664B (zh) 2022-12-02

Family

ID=66926248

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910046622.0A Active CN109884664B (zh) 2019-01-14 2019-01-14 一种城市地上生物量光学微波协同反演方法及系统

Country Status (2)

Country Link
US (1) US11454534B2 (zh)
CN (1) CN109884664B (zh)

Families Citing this family (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109884664B (zh) * 2019-01-14 2022-12-02 武汉大学 一种城市地上生物量光学微波协同反演方法及系统
CN110287457B (zh) * 2019-07-02 2023-02-17 吉林大学 基于卫星雷达遥感数据的玉米生物量反演测算方法
CN110826764B (zh) * 2019-09-24 2020-11-24 广州地理研究所 基于随机森林算法的植被指数预测方法、系统及设备
CN113205475B (zh) * 2020-01-16 2022-07-12 吉林大学 基于多源卫星遥感数据的森林高度反演方法
CN111860328B (zh) * 2020-07-21 2021-04-06 杭州时光坐标影视传媒股份有限公司 一种基于双向反射函数和森林场景光照效果建模的生物量估算方法
CN113822107B (zh) * 2020-08-11 2024-01-19 山东理工大学 一种融合多种植被指数的土地荒漠化提取方法及装置
JP7217855B2 (ja) * 2020-10-14 2023-02-06 株式会社興和 森林資源情報生成構造
CN112270675B (zh) * 2020-11-11 2021-04-16 中山大学 一种基于极化雷达遥感影像的城市内涝区域检测方法
CN112434617B (zh) * 2020-11-26 2021-08-13 南京观微空间科技有限公司 一种基于多源遥感数据的森林生物量变化监测方法及系统
CN112785590B (zh) * 2021-02-09 2022-09-13 河北地质大学 一种基于双差值归一化的植被指数计算方法
WO2022173383A1 (en) * 2021-02-11 2022-08-18 Saengprachatanarug Khwantri Method for determining height of crop covered soil surface from sea level through use of aerial photographs
CN113205565A (zh) * 2021-04-13 2021-08-03 武汉大学 基于多源数据的森林生物量遥感制图方法
CN113297904B (zh) * 2021-04-16 2023-11-21 广州市华南自然资源科学技术研究院 一种基于卫星驱动模型的高寒草地生物量估算方法及系统
CN113393541B (zh) * 2021-06-17 2022-02-15 北京博创联动科技有限公司 一种地块面积计算方法、装置、计算机终端及存储介质
CN113128134B (zh) * 2021-06-17 2021-09-14 中国矿业大学(北京) 一种矿区生态环境演变驱动因子权重量化分析方法
CN113469586B (zh) * 2021-09-02 2021-12-07 中国农业大学 土地系统状态的评估方法、装置以及存储介质
CN113865523A (zh) * 2021-09-29 2021-12-31 中国科学院东北地理与农业生态研究所 一种基于雷达和光学遥感数据同时估算裸露土壤表面水分和粗糙度的方法
CN114076738B (zh) * 2021-11-22 2023-04-18 浙江大学 一种利用遥感构建指数估算秸秆残留农田土壤盐分的方法
WO2023095039A1 (en) * 2021-11-24 2023-06-01 Dendra Systems Ltd. Generating an above ground biomass density prediction model
CN114169165B (zh) * 2021-12-06 2024-04-26 湘潭大学 一种三波段植被指数估算镉胁迫下水稻叶绿素的模型方法
CN114241331B (zh) * 2021-12-16 2024-04-05 中国科学院南京地理与湖泊研究所 以UAV为地面和Sentinel-2中介的湿地芦苇地上生物量遥感建模方法
CN114324409A (zh) * 2021-12-28 2022-04-12 南京晓庄学院 一种基于微波遥感的土壤湿度监测系统
CN114414491B (zh) * 2021-12-30 2023-10-27 青海省草原总站 一种草原生态动态监测与分析系统
CN114490905B (zh) * 2022-01-24 2024-02-09 北京师范大学 一种晴空地表净长波辐射一体化反演方法及系统
CN114429591A (zh) * 2022-01-26 2022-05-03 中国农业科学院草原研究所 一种基于机器学习的植被生物量自动监测方法及系统
CN114821349B (zh) * 2022-03-01 2023-07-07 南京林业大学 顾及谐波模型系数和物候参数的森林生物量估算方法
CN114281871B (zh) * 2022-03-01 2022-10-18 国家林业和草原局调查规划设计院 航空激光雷达估算山地林区chm树高的精度验证方法
CN114689545B (zh) * 2022-03-02 2022-11-29 珠江水利委员会珠江水利科学研究院 一种基于dsm等高切片的植被覆盖度分层估算方法、介质
CN114694036A (zh) * 2022-03-18 2022-07-01 南京农业大学 一种基于高分影像和机器学习的高海拔地区作物分类识别的方法
CN114676636B (zh) * 2022-04-01 2023-08-29 自然资源部国土卫星遥感应用中心 一种综合植被和生境特征的草原区土壤水分快速反演方法
CN114965300B (zh) * 2022-05-23 2024-04-19 中国科学院东北地理与农业生态研究所 基于光学水体类型和bp神经网络算法构建bp-turb的湖泊浊度制图方法
CN114862941B (zh) * 2022-05-24 2024-05-03 中国农业大学 阈值修正的孔隙度模型计算棉花叶面积指数的方法及系统
CN114998728A (zh) * 2022-05-24 2022-09-02 中国农业大学 无人机多源遥感预测棉花叶面积指数的方法及系统
CN114819737B (zh) * 2022-05-26 2023-10-17 中交第二公路勘察设计研究院有限公司 公路路域植被的碳储量估算方法、系统及存储介质
CN115063437B (zh) * 2022-06-16 2023-04-28 广西壮族自治区气象科学研究所 一种红树林冠层可见光图像指数特征分析方法及系统
CN115115948B (zh) * 2022-07-26 2024-03-29 云南大学 一种基于随机森林和辅助因子的林地信息精细化提取方法
CN115512244B (zh) * 2022-08-19 2024-01-05 中国林业科学研究院资源信息研究所 一种单株树木碳储量确定方法及系统
CN115561181B (zh) * 2022-10-26 2024-05-07 中科合肥智慧农业协同创新研究院 一种基于无人机多光谱数据的水质反演方法
CN115840224B (zh) * 2022-12-07 2023-09-29 西南交通大学 一种用于反演土壤水分的多目标函数变化检测方法
CN115797760B (zh) * 2023-01-29 2023-05-12 水利部交通运输部国家能源局南京水利科学研究院 主被动融合的水质立体遥感反演方法、系统及存储介质
CN116385871B (zh) * 2023-02-28 2024-03-15 中国林业科学研究院资源信息研究所 一种森林生物量遥感估算方法、系统、设备及介质
CN116563718B (zh) * 2023-07-11 2023-09-05 成都垣景科技有限公司 一种基于遥感测绘的碳储量估算方法
CN117347282A (zh) * 2023-08-22 2024-01-05 中南大学 星基气溶胶光学厚度反演方法、装置及系统和存储介质
CN117494066A (zh) * 2023-11-15 2024-02-02 黑龙江省网络空间研究中心(黑龙江省信息安全测评中心、黑龙江省国防科学技术研究院) 长时间时序地上生物量反演量化评价方法
CN117292267B (zh) * 2023-11-27 2024-02-02 武汉大学 一种基于物候信息的水稻地上生物量分段估算方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101308215A (zh) * 2008-06-11 2008-11-19 武汉大学 一种ⅱ类水体离水辐亮度反演方法
CA2639677A1 (en) * 2008-09-16 2010-03-16 Gilles Roy Standoff determination of the size and concentration of low concentration aerosols
CN104483271A (zh) * 2014-12-19 2015-04-01 武汉大学 光学反射模型与微波散射模型协同的森林生物量反演方法
CN104656098A (zh) * 2015-02-03 2015-05-27 南京林业大学 一种遥感森林生物量反演的方法
CN105527616A (zh) * 2015-10-27 2016-04-27 中国科学院遥感与数字地球研究所 一种基于多极化RADARSAT-2和Landsat8数据的植被覆盖地表土壤含水量反演方法
CN105608293A (zh) * 2016-01-28 2016-05-25 武汉大学 融合光谱与纹理特征的森林地上生物量反演方法及系统
CN108009384A (zh) * 2017-12-26 2018-05-08 中国科学院东北地理与农业生态研究所 一种城市森林碳储量景观尺度推演方法
CN108921885A (zh) * 2018-08-03 2018-11-30 南京林业大学 一种综合三类数据源联合反演森林地上生物量的方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AUPR301401A0 (en) * 2001-02-09 2001-03-08 Commonwealth Scientific And Industrial Research Organisation Lidar system and method
WO2016025848A1 (en) * 2014-08-15 2016-02-18 Monsanto Technology Llc Apparatus and methods for in-field data collection and sampling
CA3084716A1 (en) * 2017-10-25 2019-05-02 Advanced Biological Marketing, Inc. Endophytic microbial seed treatment formulations and methods related thereto for improved plant performance
US11373399B2 (en) * 2018-07-26 2022-06-28 Hall Enterprise Llc Method of crop analysis using drone with flying and driving capability
CN109884664B (zh) * 2019-01-14 2022-12-02 武汉大学 一种城市地上生物量光学微波协同反演方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101308215A (zh) * 2008-06-11 2008-11-19 武汉大学 一种ⅱ类水体离水辐亮度反演方法
CA2639677A1 (en) * 2008-09-16 2010-03-16 Gilles Roy Standoff determination of the size and concentration of low concentration aerosols
CN104483271A (zh) * 2014-12-19 2015-04-01 武汉大学 光学反射模型与微波散射模型协同的森林生物量反演方法
CN104656098A (zh) * 2015-02-03 2015-05-27 南京林业大学 一种遥感森林生物量反演的方法
CN105527616A (zh) * 2015-10-27 2016-04-27 中国科学院遥感与数字地球研究所 一种基于多极化RADARSAT-2和Landsat8数据的植被覆盖地表土壤含水量反演方法
CN105608293A (zh) * 2016-01-28 2016-05-25 武汉大学 融合光谱与纹理特征的森林地上生物量反演方法及系统
CN108009384A (zh) * 2017-12-26 2018-05-08 中国科学院东北地理与农业生态研究所 一种城市森林碳储量景观尺度推演方法
CN108921885A (zh) * 2018-08-03 2018-11-30 南京林业大学 一种综合三类数据源联合反演森林地上生物量的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
基于多源数据的根河实验区生物量反演研究;李春梅等;《北京林业大学学报》;20160315(第03期);正文第64-71页 *
机载激光雷达森林叶面积指数反演研究;骆社周等;《地球物理学报》;20130515(第05期);正文第1467-1473页 *
森林地上生物量的极化相干层析估计方法;罗环敏等;《遥感学报》;20111125(第06期);正文第1147-1154页 *
湿地植被地上生物量遥感估算方法研究进展;赵天舸等;《生态学杂志》;20160715(第07期);正文第1936-1944页 *

Also Published As

Publication number Publication date
US11454534B2 (en) 2022-09-27
US20200225075A1 (en) 2020-07-16
CN109884664A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN109884664B (zh) 一种城市地上生物量光学微波协同反演方法及系统
Shao et al. Stacked sparse autoencoder modeling using the synergy of airborne LiDAR and satellite optical and SAR data to map forest above-ground biomass
Avitabile et al. Capabilities and limitations of Landsat and land cover data for aboveground woody biomass estimation of Uganda
CN114076738B (zh) 一种利用遥感构建指数估算秸秆残留农田土壤盐分的方法
CN108268735A (zh) 基于多源遥感卫星融合数据的地表土壤水分降尺度方法
CN111598045B (zh) 一种基于对象图谱和混合光谱的遥感耕地变化检测方法
CN112381013B (zh) 基于高分辨率遥感影像的城市植被反演方法及系统
CN106501186B (zh) 一种土壤含水量产品降尺度方法
Melin et al. Assessing the performance of aerial image point cloud and spectral metrics in predicting boreal forest canopy cover
CN113205475A (zh) 基于多源卫星遥感数据的森林高度反演方法
CN108981616B (zh) 一种由无人机激光雷达反演人工林有效叶面积指数的方法
CN114926748A (zh) Sentinel-1/2微波与光学多光谱影像结合的大豆遥感识别方法
CN102052918B (zh) 基于湖泊水生植被蔓延的沼泽化动态监测与预警技术
Lang et al. High carbon stock mapping at large scale with optical satellite imagery and spaceborne LIDAR
CN116229289A (zh) 一种基于无人机的冬小麦生物量快速高精准估测方法
CN113205014B (zh) 一种基于图像锐化的时序数据耕地提取方法
CN117035066A (zh) 一种耦合地理加权与随机森林的地表温度降尺度方法
CN115830464A (zh) 基于多源数据的高原山地农业大棚自动提取方法
CN115128013A (zh) 一种基于分区算法的土壤有机质含量空间预测评估方法
CN115527108A (zh) 一种基于多时相Sentinel-2快速识别水土流失人为扰动地块的方法
Lee et al. Two-step carbon storage estimation in human settlements in urban areas using airborne LiDAR and Sentinel-2 data based on machine learning
CN109359394B (zh) 土壤湿度降尺度因子模型构建方法及系统
Guo et al. Estimating aboveground biomass of alpine grassland during the wilting period using in situ hyperspectral, Sentinel-2 and Sentinel-1 data
Memduhoğlu Identifying impervious surfaces for rainwater harvesting feasibility using unmanned aerial vehicle imagery and machine learning classification
CN117423011A (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