CN112446522A - 面对多尺度分割的产草量估算方法、装置及存储介质 - Google Patents
面对多尺度分割的产草量估算方法、装置及存储介质 Download PDFInfo
- Publication number
- CN112446522A CN112446522A CN201910822474.7A CN201910822474A CN112446522A CN 112446522 A CN112446522 A CN 112446522A CN 201910822474 A CN201910822474 A CN 201910822474A CN 112446522 A CN112446522 A CN 112446522A
- Authority
- CN
- China
- Prior art keywords
- data
- random forest
- remote sensing
- sample
- grass
- 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.)
- Pending
Links
- 244000025254 Cannabis sativa Species 0.000 title claims abstract description 87
- 230000011218 segmentation Effects 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000003860 storage Methods 0.000 title claims abstract description 9
- 238000007637 random forest analysis Methods 0.000 claims abstract description 81
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 43
- 238000004519 manufacturing process Methods 0.000 claims abstract description 29
- 238000007781 pre-processing Methods 0.000 claims abstract description 27
- 239000013598 vector Substances 0.000 claims abstract description 16
- 238000012549 training Methods 0.000 claims description 46
- 238000012360 testing method Methods 0.000 claims description 25
- 238000003066 decision tree Methods 0.000 claims description 21
- 238000012937 correction Methods 0.000 claims description 19
- 238000011156 evaluation Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000010801 machine learning Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 239000012535 impurity Substances 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 description 12
- 238000009826 distribution Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 7
- 230000003595 spectral effect Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000004927 fusion Effects 0.000 description 4
- 238000002310 reflectometry Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 101100284396 Drosophila melanogaster Hayan gene Proteins 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000000699 topical effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Life Sciences & Earth Sciences (AREA)
- Economics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法,包括以下步骤:获取遥感数据并进行预处理;对预处理后的数据进行多尺度分割获得分割结果;获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;使用分割结果作为随机森林估算模型的输入向量预测产草量;还公开了一种面对多尺度分割遥感数据和随机森林算法的产草量估算装置及可读存储介质。运算高效、结果准确。
Description
技术领域
本发明涉及卫星测算技术领域,尤其涉及一种面对多尺度分割遥感数据和随机森林的产草量估算方法、装置及存储介质。
背景技术
草地资源是全球陆地生态系统的重要组成部分,在生态环境中起着举足轻重的作用。对草地资源监测有助于了解掌握草地实际情况,合理开发利用草地资源,达到维护草原生态系统平衡的目的。草地遥感估产是通过卫星传感器或地面光谱仪获取地球表面信息,在地物光谱理论指导下,将获取的地球表面信息经过复杂的综合处理,以识别草地及其长势,实现草地面积和长势监测及单产和总产的预报。使用遥感技术进行草地估产最早起源于国外,国内的研究起步较晚,但是发展很快。经研究,最初人们基于AVHRR-NDVI数据进行草原遥感估产,并且通过建立回归模型探究植被指数NDVI和产草量之间的关系。随着技术的发展,利用地面光谱实验数据或高光谱数据、地面监测样本数据结合Landsat、MODIS等数据产品,通过建立不同的反演模型估测草地产草量、探究草地生长状况。
哨兵2号(Sentinel-2)是欧洲太空局于2015年6月发射的多光谱成像卫星,携带一枚多光谱成像仪(MSI),具有最高可达10m分辨率、可覆盖13个光谱波段的遥感数据,用于陆地监测,可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像,还可用于紧急救援服务。目前已发送Sentinel-2A和Sentinel-2B两颗卫星,一颗卫星的重访周期为10天,两颗互补,重访周期为5天。但是由于Sentinel-2是新发射的卫星,目前将Sentinel-2数据用于草原地区遥感估产的研究较少。
高分一号(GF-1)卫星是中国高分辨率对地观测卫星系统重大专项的第1颗卫星,于2013-04-26由长征二号丁运载火箭在酒泉卫星发射基地成功发射入轨,搭载了4台空间分辨率为16m的宽视场(Wide Field View,WFV)传感器,2台全色和多光谱传感器(Pan andMultispectral Sensor,PMS),其中全色影像的空间分辨率为2m,多光谱影像的空间分辨率为8m。GF-1号卫星为草原资源监测提供了新的数据源。
合成孔径雷达(synthetic aperture radar,SAR)对地表覆盖区具有一定的穿透性,具有全天时全天候工作的能力,能穿透地表浅覆盖层,对覆盖层以下的地质特征进行有效探测,穿透深度与波长、地表介电常数、地表粗糙度有关。高分三号(GF-3)卫星是我国首颗C波段全极化合成孔径雷达成像卫星,于2016年8月10日在太原卫星发射中心成功发射升空,可以全天时全天候实现海洋和陆地信息的监测。
高分六号(GF-6)是我国2018年6月2日发射的一颗低轨光学遥感卫星,具有高分辨率、宽覆盖、高质量和高效成像等特点,能有力支撑农业资源监测、林业资源调查、防灾减灾救灾等工作。GF-6可以与GF-1组网运行,将使遥感数据获取的时间分辨率从4天缩短到2天。
发明内容
鉴于目前存在的上述不足,本发明提供一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法,能够估算获得结果较准确的产草量。
为达到上述目的,本发明的实施例采用如下技术方案:
一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法,所述产草量估算方法包括以下步骤:
获取遥感数据并进行预处理;
对预处理后的数据进行多尺度分割获得分割结果;
获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;
根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;
根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;
使用分割结果作为随机森林估算模型的输入向量预测产草量。
依照本发明的一个方面,所述获取遥感数据并进行预处理包括:
Sentinel数据获取及预处理;
GF数据获取及预处理。
依照本发明的一个方面,所述Sentinel数据获取及预处理包括:对Sentinel原始数据进行辐射定标以及大气校正;采用双三次卷积的方式将影像波段重采样至质量最好的10m分辨率。
依照本发明的一个方面,所述GF数据获取及预处理包括:GF-1数据、GF-3数据和GF-6数据的获取及预处理,其中:
GF-1数据获取及预处理包括:对于GF1-1 PMS影像数据,需要进行多光谱影像辐射定标、大气校正和正射校正,再对全色影像进行辐射定标与正射校正,之后将校正后的多光谱影像与全色影像融合,得到具有多光谱信息和高空间分辨率的遥感融合影像;
GF-3数据获取及预处理包括:对GF-3影像进行辐射定标、复数据转换、多视处理、滤波以及转DB影像的预处理操作;
GF-6数据获取及预处理包括:对GF-6号数据进行辐射定标、大气校正和几何校正。
依照本发明的一个方面,在对预处理后的数据进行多尺度分割获得分割结果之前,需对预处理后的数据进行影像镶嵌和影像裁剪步骤。
依照本发明的一个方面,包括以下步骤:
将需要计算的实测点总产草量鲜重/可食产草量鲜重作为Y值存入样本之中,将从分割结果提取出的波段值以及EVI、NDVI值作为X值;
使用机器学习库进行随机森林回归模型建模,构建的随机森林回归模型用{h(X,Θk),k=1,…}表示,其中X为输入向量,{Θk}为独立同分布随机向量。
依照本发明的一个方面,所述随机森林分类器和随机森林估算模型均需要进行随机森林回归模型建模,具体包括以下步骤:
对输入样本数据集做归一化处理;
将样本数据集划分为训练集和测试集;
对于训练集D={(x1,y1),(x2,y2),…,(xn,yn)},利用Bootstrap方法重采样,随机产生T个训练集S1,S2,…,ST;
对每个训练集生成对应的决策树C1,C2,…,CT;在对每个非叶节点上选择属性前,从全部M个属性中随机抽取m个(m<<M)作为当前节点的分裂属性集,并从中选择出一个最佳分割属性作为节点进行分裂;
将生成的多棵决策树构成随机森林,对于测试集样本X,利用每棵决策树进行测试,得到预测结果C1(x),C2(x),…,CT(x);
作为随机森林估算模型时,测试集样本X的预测值是这些树的结果的平均值;
作为随机森林分类器时,预测的最终类别为该样本点所到叶节点中投票数最多的类别。
依照本发明的一个方面,所述作为随机森林估算模型时,测试集样本X的预测值是这些树的结果的平均值包括:
对任意划分特征A,对应的任意划分点S两边划分成的数据集D1,D2,求出使D1,D2各自集合均方差之和最小,同时D1,D2的均方差之和最小所对应的特征和特征值划分点,表达式为:
其中,c1为D1数据集的样本输出均值,c2为D2数据集的样本输出均值。
依照本发明的一个方面,作为随机森林分类器时,预测的最终类别为该样本点所到叶节点中投票数最多的类别包括:
RF中的每一棵分类树都为二叉树,其生成遵循自顶向下的递归分裂原则,即从根节点开始依次对训练集进行划分;在二叉树中,根节点包含全部训练数据,按照节点不纯度最小原则分裂为左右节点,它们分别包含训练数据的一个子集,按照同样的规则节点继续分裂,直到满足分支停止规则而停止生长;若节点n上的分类数据全部来自于同一类别,则此节点的不纯度为0,不纯度度量方法是Gini准则:
对于一个给定的训练集D,假设训练集中由K类,C1,C2,…Ck,则Gini系数为:
如果训练样本被特征A的是否取某一个值a而被分成D1和D2两部分,那么在特征A的条件下的基尼指数为:
相当于对根据特征A划分的两个小集合分别求基尼指数然后根据经验概率取期望得到特征A下的D的基尼指数。
依照本发明的一个方面,所述进行随机森林回归模型建模包括判定变量重要性,具体包括:
在每一棵决策树的变量中加入随机噪声,然后检验袋外误差的增减,如果误差增加,则改变量比较重要,反之则不重要;
计算方法为:
依照本发明的一个方面,包括进行模型评价,所述模型评价选择选择决定系数(R2)和均方误差(MSE)两个指标,具体包括:
依照本发明的一个方面,所述对预处理后的数据进行多尺度分割获得分割结果包括:
将遥感图像分割成若干个互不重叠同时没有孔洞的子集区域,每个区域内部具有相同或相似的特征,表征了一种地物;选择同质异质法使对象类同质性尽量大,对象间同质性尽量小;同质性由影像对像内部的面积和标准差来表示,具体公式为:
其中,ai和vi分别表示第i个区域的面积和标准差,n为该尺度下分割区域的个数,wVar表示同质性指数;
异质性通过影像对象间的空间相关性指标Moran’sI(MI)表示,具体公式为:
一种面对多尺度分割遥感数据和随机森林算法的产草量估算装置,所述产草量估算装置包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如上所述的任一面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被执行时实现如上所述的任一面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
本发明实施的优点:本发明所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,包括以下步骤:获取遥感数据并进行预处理;对预处理后的数据进行多尺度分割获得分割结果;获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;使用分割结果作为随机森林估算模型的输入向量预测产草量;无需对变量的正态性和独立性等假设条件进行检验,同时也不需要考虑多变量的共线问题,且运算高效、结果准确。有较高的准确率,对异常值和噪声具有较好的容忍度,对高光谱遥感等高维度数据训练和学习效果较佳。同时,机器学习中的一个重大问题是过拟合,对于随机森林而言,只要森林中有足够多的树,分类器就不会过度拟合模型,泛化能力强。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所述的一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法流程示意图;
图2为本发明所述的海晏研究区总产草量分布图;
图3为本发明所述的海晏研究区可食产草量分布图;
图4为本发明所述的祁连研究区总产草量分布图;
图5为本发明所述的祁连研究区可食产草量分布图;
图6为本发明所述的产草量估算装置示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一
如图1、图2、图3、图4和图5所示,一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法,所述产草量估算方法包括以下步骤:
步骤S1:获取遥感数据并进行预处理;
所述获取遥感数据并进行预处理包括:
1)Sentinel数据获取及预处理;
2)GF数据获取及预处理;
其中,Sentinel数据获取及预处理包括以下具体实施:
根据实际项目需要,以及实地采样区与采样时间,选择所需Sentinel-2A或Sentinel-2B数据,数据要求基本无云,影像质量高。通过欧洲太空局(ESA)的数据共享网站(https://scihub.copernicus.eu/dhus/)直接获取的Sentinel-2数据是Level-1C级多光谱数据,是经过几何精校正的正射影像,投影坐标系为UTM/WGS-84。ESA还对Sentinel-2L2A级数据进行了定义,L2A级数据主要包含经过辐射定标和大气校正的大气底层反射率数据,这个L2A级数据需要用户自行处理。Sentinel-2数据波段参考如下:
表1 Sentinel-2卫星数据部分参数信息
在实际应用中,需要对Sentinel原始数据进行辐射定标以及大气校正。从表1可以看出,Sentinel-2波段的空间分辨率并不统一,在这里使用双三次卷积的方式将影像波段重采样至质量最好的10m分辨率,并且根据研究区域需要,对数据进行波段组合、影像拼接和裁剪。
在实际应用中,地物辐射定标原理如下:
L=gain×DN+bias
其中,L是地物在大气顶部的辐射亮度,DN是象元值,gain是定标增益,gain和bias可从头文件中得到,ρ是地物反射率,ds是日地天文单位距离,E0太阳辐照度,θ是太阳天顶角。
Flaash大气校正原理如下:
其中,等式右边前2项分别为辐射亮度的直射、散射分量,是太阳辐射传输中地表与大气共同作用的结果;La为辐射亮度的大气程辐射分量,是大气分子和气溶胶作用的结果;ρ为像元表面反射率;ρe为像元周围平均表面反射率;S为大气球面反照率;A、B为取决于大气条件和几何条件的2个独立系数。
GF数据获取及预处理具体可包括如下实施:
GF-1、GF-3和GF-6数据均都由中国资源卫星应用中心提供。
GF-1数据将高空间分辨率、高时间分辨率等优势集合于一体,在这里选择GF-1PMS数据作为数据源,其参数设置如表2所示。两台PMS相机组合使传感器在空间分辨率高的同时具备较大的幅宽,另外传感器具有良好的侧摆成像能力,可以在应急状态下对特定的区域进行高时间频次观测。对于GF1-1 PMS影像数据,需要进行多光谱影像辐射定标、大气校正、正射校正,再对全色影像进行辐射定标与正射校正,之后将校正后的多光谱影像与全色影像融合,得到具有多光谱信息和高空间分辨率的遥感融合影像。根据具体需求,对影像进行镶嵌拼接。
表2 GF-1PMS参数设置
在这里使用超分辨率贝叶斯算法(Pansharp)对GF-1数据进行融合。该算法基于最小二乘法对参与融合的原始多光谱影像、全色影像与融合后的多光谱影像、全色影像之间的灰度值关系进行最佳匹配,并采用统计的方法来解决融合过程标准化和自动化的问题,使得融合影像在获得高空间分辨率的同时,还具备良好的光谱保真度。该算法因不受波段数限制,光谱和细节特征保真度高,是目前高分辨率遥感影像融合中公认比较好的算法之一。
GF-3号卫星聚束、条带、扫描、全极化、波成像、全球观测、高低入射角12种成像模式,最高分辨率可达1m。这里选择成像方式如表3所示,并对GF-3影像进行了辐射定标、复数据转换、多视处理、滤波以及转DB影像(输出TIFF)等一系列预处理操作。
表3 可用于产草量估算的GF-3号卫星成像模式
对GF-3数据辐射定标一般是把后向散射能量强度转换为后向散射系数,定标公式为:
其中,σ0为后向散射系数,Pt为雷达发射器发射功率,Gt为雷达信号接收通道增益,Gθ为双程天线方向图,λ为发射波波长,δa为方位向分辨率,δt为距离向分辨率,M为像处理器的功率增强,R为斜距,β为角度(入射),I′为热噪声功率。
SAR影像最主要的特征之一是强度或者幅度,利用强度/幅度特征,可以提取地物目标信息,因此需要将SAR复数数据转换为SAR强度/幅度数据。SAR数据是由实部和虚部构成的复数据,也称为同相(In-phase)和正交通道(Quadrature channels)。强度公式为:
其中:A=Amplitude,I=In-phase,Q=Quadrature channels。强度和幅度之间的转换公式为:Intensity=Amplitude2,其中:Intensity为强度,Amplitude为振幅(幅度)。
多视处理,可以提高图像的可视效果,并优化每个像元后向散射的估计精度,具体是对邻域像素相干和的平均,公式如下:
使用Non-Local Means滤波对GF-3进行滤波处理,Non-Local Means滤波是一种极化SAR影像滤波方法,对斑点噪声有高抑制性和细节信息有高保真性。其原理是通过寻找整幅图像中与待恢复区域相似得块,再将找出得所有相似块通过某种相似性准则(如K均值聚类)进行聚类后联合滤波,从而减少图像噪声。而对于图像相似性的计算,可以通过比较一个窗口内的各个像素值来得到。Non-local Means算法公式为:
其中,为去噪后的图像,v(y)为含噪声图像,w(x,y)为权值,Ωx为像素x的邻域。其中w(x,y)的确定是该算法的关键。w(x,y)表示在原始图像v(y)中,像素x和像素y的相似度,w(x,y)的值由两个窗口v(x)与v(y)之间的欧式距离的平方决定,如下:
其中:
Z(x)表示归一化系数,h为平滑参数h为平滑参数。图像块之间的相似性由欧氏距离决定,从w(x,y)计算式可以看出,欧氏距离越大,该算法去噪能力对平滑参数的选择依赖性越高。
GF-6配置2米全色/8米多光谱高分辨率相机、16米多光谱中分辨率宽幅相机,2米全色/8米多光谱相机观测幅宽90公里,16米多光谱相机观测幅宽800公里。GF-6还实现了8谱段CMOS探测器的国产化研制,国内首次增加了能够有效反映作物特有光谱特性的“红边”波段。对GF-6号数据同样进行辐射定标、大气校正和几何校正。
步骤S2:对预处理后的数据进行多尺度分割获得分割结果;
所述对预处理后的数据进行多尺度分割获得分割结果包括:
遥感影像的分割是将其分割成若干个互不重叠同时没有孔洞的子集区域,每个区域内部具有相同或相似的特征,表征了一种地物。在这里选择同质异质法使对象类同质性尽量大,对象间同质性尽量小。同质性由影像对像内部的面积和标准差来表示,具体公式为:
其中,ai和vi分别表示第i个区域的面积和标准差,n为该尺度下分割区域的个数。引入面积权重有利于优先合并面积较小的区域,防止面积较小的区域产生的不稳定性。同质性指数wVar越小,则表示影像对象内部同质性越高。
异质性通过影像对象间的空间相关性指标Moran’sI(MI)表示,具体公式为:
其中,n为该尺度下分割区域的个数,wij为空间邻接关系权重,yi为区域i的光谱均值,为整个影像的光谱均值。如果区域i,j相邻则wi,j=1,否则wi,j=0。异质性指数MI越小,则表示影像对象的空间相关性越低,异质性越高。
所述分割结果包括:Sentinel-2特征波段,GF-1PMS特征波段,GF-3特征波段,GF-6特征波段。
步骤S3:获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;
采用的实测数据主要包括草本样方编号、经纬度、植被盖度、主要植物种类、总产草量鲜重及风干重、可食产草量鲜重及风干重、调查时间等信息。根据牧草生长规律,调查时间选择集中在7-8月为佳,这是牧草生长最旺盛的时期。并且考虑产草量估算精度,实测数据分布应尽可能均匀,并且为研究区内的草地类型选择数量足够的样方进行实测。
根据实测数据分别创建草地类别分类样本和产草量估算样本。具体包括:
植被指数是将遥感影像上不同波段的光谱反射率通过线性或非线性方式组合,得到反应地表植被状况的度量。根据过往研究,选择NDVI、EVI两种能较好的反映地表植被覆盖情况的植被指数作为特征值。根据样本点坐标,提取样本点对应波段值以及增强植被指数EVI、归一化植被指数NDVI作为样本数据,存储格式为.csv。
其中,ρnir为近红外波段值,ρred为红光波段值,ρblue为蓝光波段值。
在本实施例中,需要对植被指数NDVI进行判断,判断NDVI是否>0.15,若大于0.15则说明有植被,该数据可参与后续步骤;若小于0.15则说明无植被,数据不采用。
步骤S4:根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;
根据草地类别分类样本,结合预处理并多尺度分割后的遥感数据,通过建立随机森林模型进行草地类型分类。在这里使用python的机器学习库scikit-learn(sklearn)进行随机森林模型建模,sklearn支持包括分类、回归、降维和聚类四大机器学习算法,还包含了特征提取、数据处理和模型评估三大模块。所述随机森林模型的建立包括以下步骤:
数据归一化
由于选取的评价值之间有不同得量纲和量纲单位,这种情况会影响最终结果,为了消除这种影响,需要对输入样本数据集做归一化处理,使样本特征分布接近标准正态分布,以达到更好得效果。在这里使用MinMacScaler将特征缩放到给定得最小值和最大值之间,通常取值(0,1)。
其中,xi为该点实测样本值,xmin(axis=0)为最小样本,xmax(axis=0)为最大样本,max,min是给定缩放范围的最大值最小值。
样本训练集与测试集选取
在进行建模之前,需要将样本数据集划分为训练集和测试集,在划分训练集和测试集时,需要进行权衡,如果测试集数据越小,对模型得泛化误差得估计越不准确。一般在实际应用中,基于整个数据集的大小考虑,训练集数据和测试集数据的划分比例是6:4、7:3或8:2。对于庞大的数据可以使用9:1,甚至是99:1。
RF模型建立
RF是在Bagging基础上做了修改演化而来的。Bagging算法是在原始的数据集上采用有放回的随机取样的方式来抽取m个子样本,从而利用这m个子样本训练m个学习器,降低了模型的方差,然后将数据放在这m个分类器上,最后根据这个m个分类器的投票结果,决定数据属于哪一类。RF在此基础上改动了两处:一是在训练每个学习器时,不是从所有特征中选择最优特征来进行节点切分,而是随机的选取k个特征,从这k个特征中选择最优特征来切分节点;二是使用CART决策树。
RF回归模型用{h(X,Θk),k=1,…}表示,其中X为输入向量,{Θk}为独立同分布随机向量。
具体算法步骤如下:
对于训练数据集D={(x1,y1),(x2,y2),…,(xn,yn)},
(1)利用Bootstrap方法重采样,随机产生T个训练集S1,S2,…,ST;
(2)对每个训练集生成对应的决策树C1,C2,…,CT;在对每个非叶节点上选择属性前,从全部M个属性中随机抽取m个(m<<M)作为当前节点的分裂属性集,并从中选择出一个最佳分割属性作为节点进行分裂。这里m控制了随机性的引入程度,若m=M则基决策树的构建与传统决策树相同;若令m=1,则是随机选择一个属性进行划分,一般推荐m=log2M;
(3)将生成的多棵决策树构成随机森林,对于测试集样本X,利用每棵决策树进行测试,得到预测结果C1(x),C2(x),…,CT(x);
(4)预测的最终类别为该样本点所到叶节点中投票数最多的类别。
CART分类树创建
RF中的每一棵分类树都为二叉树,其生成遵循自顶向下的递归分裂原则,即从根节点开始依次对训练集进行划分;在二叉树中,根节点包含全部训练数据,按照节点不纯度最小原则分裂为左右节点,它们分别包含训练数据的一个子集,按照同样的规则节点继续分裂,直到满足分支停止规则而停止生长;若节点n上的分类数据全部来自于同一类别,则此节点的不纯度为0,不纯度度量方法是Gini准则:
对于一个给定的训练集D,假设训练集中由K类,C1,C2,…Ck,则Gini系数为:
如果训练样本被特征A的是否取某一个值a而被分成D1和D2两部分,那么在特征A的条件下的基尼指数为:
相当于对根据特征A划分的两个小集合分别求基尼指数然后根据经验概率取期望得到特征A下的D的基尼指数。
RF模型判定变量重要性的方法使在每一棵决策树的变量中加入随机噪声,然后检验袋外误差的增减,如果误差增加,则改变量比较重要,反之则不重要。计算方法为:
模型评价
模型评价选择决定系数(R2)和均方误差(MSE)两个指标。其中,R2表征回归方程在多大程度上解释了因变量的变化,或者说模型对观测值的拟合程度如何。
步骤S5:根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;
所述步骤S5的具体实现过程可如下:
根据产草量估算样本。将需要计算的实测点总产草量鲜重/可食产草量鲜重作为Y值存入样本之中,将提取出的波段值以及EVI、NDVI值作为X值,在这里使用python的机器学习库scikit-learn(sklearn)进行随机森林回归模型建模,sklearn支持包括分类、回归、降维和聚类四大机器学习算法,还包含了特征提取、数据处理和模型评估三大模块。
所述随机森林模型的构建如下:
数据归一化
由于选取的评价值之间有不同得量纲和量纲单位,这种情况会影响最终结果,为了消除这种影响,需要对输入样本数据集做归一化处理,使样本特征分布接近标准正态分布,以达到更好得效果。在这里使用MinMacScaler将特征缩放到给定得最小值和最大值之间,通常取值(0,1)。
其中,xi为该点实测样本值,xmin(axis=0)为最小样本,xmax(axis=0)为最大样本,max,min是给定缩放范围的最大值最小值。
样本训练集与测试集选取
在进行建模之前,需要将样本数据集划分为训练集和测试集,在划分训练集和测试集时,需要进行权衡,如果测试集数据越小,对模型得泛化误差得估计越不准确。一般在实际应用中,基于整个数据集的大小考虑,训练集数据和测试集数据的划分比例是6:4、7:3或8:2。对于庞大的数据可以使用9:1,甚至是99:1。
RF模型建立
RF是在Bagging基础上做了修改演化而来的。Bagging算法是在原始的数据集上采用有放回的随机取样的方式来抽取m个子样本,从而利用这m个子样本训练m个学习器,降低了模型的方差,然后将数据放在这m个分类器上,最后根据这个m个分类器的投票结果,决定数据属于哪一类。RF在此基础上改动了两处:一是在训练每个学习器时,不是从所有特征中选择最优特征来进行节点切分,而是随机的选取k个特征,从这k个特征中选择最优特征来切分节点;二是使用CART决策树。
RF回归模型用{h(X,Θk),k=1,…}表示,其中X为输入向量,{Θk}为独立同分布随机向量。
具体算法步骤如下:
对于训练数据集D={(x1,y1),(x2,y2),…,(xn,yn)},
(1)利用Bootstrap方法重采样,随机产生T个训练集S1,S2,…,ST;
(2)对每个训练集生成对应的决策树C1,C2,…,CT;在对每个非叶节点上选择属性前,从全部M个属性中随机抽取m个(m<<M)作为当前节点的分裂属性集,并从中选择出一个最佳分割属性作为节点进行分裂。这里m控制了随机性的引入程度,若m=M则基决策树的构建与传统决策树相同;若令m=1,则是随机选择一个属性进行划分,一般推荐m=log2M;
(3)将生成的多棵决策树构成随机森林,对于测试集样本X,利用每棵决策树进行测试,得到预测结果C1(x),C2(x),…,CT(x);
(4)对于回归问题,测试集样本X的预测值是这些树的结果的平均值。
CART回归树创建
使用CART构建回归树时,采用的原则是最小方差。即对任意划分特征A,对应的任意划分点S两边划分成的数据集D1,D2,求出使D1,D2各自集合和均方差最小,同时D1,D2的均方差之和最小所对应的特征和特征值划分点。表达式为:
其中,c1为D1数据集的样本输出均值,c2为D2数据集的样本输出均值。CART树的预测是根据叶子结点的均值,因此RF的预测是所有树的预测值的平均值。
模型评价
模型评价选择决定系数(R2)和均方误差(MSE)两个指标。其中,R2表征回归方程在多大程度上解释了因变量的变化,或者说模型对观测值的拟合程度如何。
步骤S6:使用分割结果作为随机森林估算模型的输入向量预测产草量。
所述步骤S6使用分割结果作为随机森林估算模型的输入向量预测产草量的具体实施步骤如下:
采用的遥感数据为进行了包含分割的预处理后的遥感数据,使用gdal类库读取预处理之后的遥感数据,将影像波段值按样本的格式存入数组作为输入向量。进行影像多尺度分割、植被覆盖判断、草地类型分类、分类别草地产草量估算,并将预测出的总产草量结果输出为TIFF影像。
实施例二
如图6所示,一种面对多尺度分割遥感数据和随机森林算法的产草量估算装置,所述产草量估算装置包括:
存储器100,用于存储计算机程序;
处理器200,用于执行所述计算机程序时实现如上实施例所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
实施例三
一种可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序被执行时实现如上实施例所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
在实际应用中,包括以下实施数据:
设定有两个研究区,分别为海晏研究区和祁连研究区;
其中,海晏研究区地处100.708769~101.136171°E,36.897002~37.188647°N,主要分布草地类型是高寒草甸和温性草原;祁连研究区地处100.495069~100.860690°E,37.588238~37.801124°N,主要分布草地类型是高寒草甸和温性草原。
进行RF模型建立,并且对总产草量和可食产草量分别建模,按照实施例一和实施例二的步骤及流程图方法估算出产草量,并将结果保存为TIFF影像。如图2所示,为海晏研究区的总产草量专题图集;如图3所示,为海晏研究区的可食产草量专题图集。如图4所示,为祁连研究区的总产草量专题图集;如图5所示,为祁连研究区的可食产草量专题图集。
无需对变量的正态性和独立性等假设条件进行检验,同时也不需要考虑多变量的共线问题,且运算高效、结果准确。有较高的准确率,对异常值和噪声具有较好的容忍度,对高光谱遥感等高维度数据训练和学习效果较佳。同时,机器学习中的一个重大问题是过拟合,对于随机森林而言,只要森林中有足够多的树,分类器就不会过度拟合模型,泛化能力强。
本发明实施的优点:本发明所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,包括以下步骤:获取遥感数据并进行预处理;对预处理后的数据进行多尺度分割获得分割结果;获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;使用分割结果作为随机森林估算模型的输入向量预测产草量;无需对变量的正态性和独立性等假设条件进行检验,同时也不需要考虑多变量的共线问题,且运算高效、结果准确。有较高的准确率,对异常值和噪声具有较好的容忍度,对高光谱遥感等高维度数据训练和学习效果较佳。同时,机器学习中的一个重大问题是过拟合,对于随机森林而言,只要森林中有足够多的树,分类器就不会过度拟合模型,泛化能力强。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域技术的技术人员在本发明公开的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (14)
1.一种面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述产草量估算方法包括以下步骤:
获取遥感数据并进行预处理;
对预处理后的数据进行多尺度分割获得分割结果;
获取产草区样本点实测数据并根据实测数据分别创建草地类别分类样本和产草量估算样本;
根据分割结果和草地类别分类样本构建随机森林分类器并获得草地类型分类结果;
根据草地类型分类结果、分割结果和产草量估算样本建立随机森林估算模型;
使用分割结果作为随机森林估算模型的输入向量预测产草量。
2.根据权利要求1所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述获取遥感数据并进行预处理包括:
Sentinel数据获取及预处理;
GF数据获取及预处理。
3.根据权利要求2所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述Sentinel数据获取及预处理包括:对Sentinel原始数据进行辐射定标以及大气校正;采用双三次卷积的方式将影像波段重采样至质量最好的10m分辨率。
4.根据权利要求2所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述GF数据获取及预处理包括:GF-1数据、GF-3数据和GF-6数据的获取及预处理,其中:
GF-1数据获取及预处理包括:对于GF1-1 PMS影像数据,需要进行多光谱影像辐射定标、大气校正和正射校正,再对全色影像进行辐射定标与正射校正,之后将校正后的多光谱影像与全色影像融合,得到具有多光谱信息和高空间分辨率的遥感融合影像;
GF-3数据获取及预处理包括:对GF-3影像进行辐射定标、复数据转换、多视处理、滤波以及转DB影像的预处理操作;
GF-6数据获取及预处理包括:对GF-6号数据进行辐射定标、大气校正和几何校正。
5.根据权利要求2所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,在对预处理后的数据进行多尺度分割获得分割结果之前,需对预处理后的数据进行影像镶嵌和影像裁剪步骤。
6.根据权利要求1所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,包括以下步骤:
将需要计算的实测点总产草量鲜重/可食产草量鲜重作为Y值存入样本之中,将从分割结果提取出的波段值以及EVI、NDVI值作为X值;
使用机器学习库进行随机森林回归模型建模,构建的随机森林回归模型用{h(X,Θk),k=1,…}表示,其中X为输入向量,{Θk}为独立同分布随机向量。
7.根据权利要求1所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述随机森林分类器和随机森林估算模型均需要进行随机森林回归模型建模,具体包括以下步骤:
对输入样本数据集做归一化处理;
将样本数据集划分为训练集和测试集;
对于训练集D={(x1,y1),(x2,y2),…,(xn,yn)},利用Bootstrap方法重采样,随机产生T个训练集S1,S2,…,ST;
对每个训练集生成对应的决策树C1,C2,…,CT;在对每个非叶节点上选择属性前,从全部M个属性中随机抽取m个(m<<M)作为当前节点的分裂属性集,并从中选择出一个最佳分割属性作为节点进行分裂;
将生成的多棵决策树构成随机森林,对于测试集样本X,利用每棵决策树进行测试,得到预测结果C1(x),C2(x),…,CT(x);
作为随机森林估算模型时,测试集样本X的预测值是这些树的结果的平均值;
作为随机森林分类器时,预测的最终类别为该样本点所到叶节点中投票数最多的类别。
9.根据权利要求7所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,作为随机森林分类器时,预测的最终类别为该样本点所到叶节点中投票数最多的类别包括:
RF中的每一棵分类树都为二叉树,其生成遵循自顶向下的递归分裂原则,即从根节点开始依次对训练集进行划分;在二叉树中,根节点包含全部训练数据,按照节点不纯度最小原则分裂为左右节点,它们分别包含训练数据的一个子集,按照同样的规则节点继续分裂,直到满足分支停止规则而停止生长;若节点n上的分类数据全部来自于同一类别,则此节点的不纯度为0,不纯度度量方法是Gini准则:
对于一个给定的训练集D,假设训练集中由K类,C1,C2,…Ck,则Gini系数为:
如果训练样本被特征A的是否取某一个值a而被分成D1和D2两部分,那么在特征A的条件下的基尼指数为:
相当于对根据特征A划分的两个小集合分别求基尼指数然后根据经验概率取期望得到特征A下的D的基尼指数。
12.根据权利要求1至11之一所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法,其特征在于,所述对预处理后的数据进行多尺度分割获得分割结果包括:
将遥感图像分割成若干个互不重叠同时没有孔洞的子集区域,每个区域内部具有相同或相似的特征,表征了一种地物;选择同质异质法使对象类同质性尽量大,对象间同质性尽量小;同质性由影像对像内部的面积和标准差来表示,具体公式为:
其中,ai和vi分别表示第i个区域的面积和标准差,n为该尺度下分割区域的个数,wVar表示同质性指数;
异质性通过影像对象间的空间相关性指标Moran’sI(MI)表示,具体公式为:
13.一种面对多尺度分割遥感数据和随机森林算法的产草量估算装置,其特征在于,所述产草量估算装置包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至12任一项所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
14.一种可读存储介质,其特征在于,所述可读存储介质上存储有计算机程序,所述计算机程序被执行时实现如权利要求1至12任一项所述的面对多尺度分割遥感数据和随机森林算法的产草量估算方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910822474.7A CN112446522A (zh) | 2019-09-02 | 2019-09-02 | 面对多尺度分割的产草量估算方法、装置及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910822474.7A CN112446522A (zh) | 2019-09-02 | 2019-09-02 | 面对多尺度分割的产草量估算方法、装置及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112446522A true CN112446522A (zh) | 2021-03-05 |
Family
ID=74734967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910822474.7A Pending CN112446522A (zh) | 2019-09-02 | 2019-09-02 | 面对多尺度分割的产草量估算方法、装置及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112446522A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113049530A (zh) * | 2021-03-17 | 2021-06-29 | 北京工商大学 | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 |
CN113160262A (zh) * | 2021-03-23 | 2021-07-23 | 珠海欧比特宇航科技股份有限公司 | 基于高分辨率卫星遥感影像的蚝排提取方法、系统及介质 |
CN113761790A (zh) * | 2021-07-27 | 2021-12-07 | 河海大学 | 一种基于Stacking集成学习的果树叶片氮含量估算方法 |
CN113780076A (zh) * | 2021-08-05 | 2021-12-10 | 北京市测绘设计研究院 | 建筑垃圾的图像识别方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101728137B1 (ko) * | 2016-02-04 | 2017-04-19 | (주)한라지리정보 | 위성 영상과 gis를 사용한 토지피복 항목별 영상분류 방법 |
CN108846338A (zh) * | 2018-05-29 | 2018-11-20 | 南京林业大学 | 基于面向对象随机森林的极化特征选择及分类方法 |
US20180366225A1 (en) * | 2017-06-19 | 2018-12-20 | Viz.ai, Inc. | Method and system for computer-aided triage |
-
2019
- 2019-09-02 CN CN201910822474.7A patent/CN112446522A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101728137B1 (ko) * | 2016-02-04 | 2017-04-19 | (주)한라지리정보 | 위성 영상과 gis를 사용한 토지피복 항목별 영상분류 방법 |
US20180366225A1 (en) * | 2017-06-19 | 2018-12-20 | Viz.ai, Inc. | Method and system for computer-aided triage |
CN108846338A (zh) * | 2018-05-29 | 2018-11-20 | 南京林业大学 | 基于面向对象随机森林的极化特征选择及分类方法 |
Non-Patent Citations (6)
Title |
---|
LI ET AL: "A method for scale parameter selection and segments refinement for multi-resolution image segmentation", 《MULTISPECTRAL, HYPERSPECTRAL, AND ULTRASPECTRAL REMOTE SENSING TECHNOLOGY, TECHNIQUES AND APPLICATIONS V》, vol. 9263, 19 March 2015 (2015-03-19), pages 926301 * |
LIU ET AL: "MODELLING METHOD OF PREDICTION MODEL FOR SALT FIELD ION CONCENTRATION UNDER SOLAR THERMAL SYSTEM USING RANDOM FOREST", 《THERMAL SCIENCE》, vol. 23, no. 5, 18 May 2019 (2019-05-18), pages 2623 - 2630 * |
WU ET AL: "Surface Roughness Prediction in Additive Manufacturing Using Machine Learning", 《PROCEEDINGS OF THE ASME 2018 13TH INTERNATIONAL MANUFACTURING SCIENCE AND ENGINEERING CONFERENCE》, 18 June 2018 (2018-06-18), pages 1 - 6 * |
李晓红;陈尔学;李增元;李世明;: "综合应用多源遥感数据的面向对象土地覆盖分类方法", 林业科学, no. 02, 15 February 2018 (2018-02-15) * |
欧阳玲;毛德华;王宗明;李慧颖;满卫东;贾明明;刘明月;张淼;刘焕军;: "基于GF-1与Landsat8 OLI影像的作物种植结构与产量分析", 农业工程学报, no. 11, 8 June 2017 (2017-06-08), pages 147 - 156 * |
欧阳玲等: "基于GF-1与Landsat8 OLI影像的作物种植结构与产量分析", 《农业工程学报》, vol. 33, no. 11, 8 June 2017 (2017-06-08), pages 147 - 156 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113049530A (zh) * | 2021-03-17 | 2021-06-29 | 北京工商大学 | 一种基于近红外高光谱的单粒玉米种子水分含量检测方法 |
CN113160262A (zh) * | 2021-03-23 | 2021-07-23 | 珠海欧比特宇航科技股份有限公司 | 基于高分辨率卫星遥感影像的蚝排提取方法、系统及介质 |
CN113761790A (zh) * | 2021-07-27 | 2021-12-07 | 河海大学 | 一种基于Stacking集成学习的果树叶片氮含量估算方法 |
CN113761790B (zh) * | 2021-07-27 | 2024-04-23 | 河海大学 | 一种基于Stacking集成学习的果树叶片氮含量估算方法 |
CN113780076A (zh) * | 2021-08-05 | 2021-12-10 | 北京市测绘设计研究院 | 建筑垃圾的图像识别方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112446522A (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 | |
Hame et al. | Improved mapping of tropical forests with optical and SAR imagery, Part II: Above ground biomass estimation | |
Bagan et al. | Combination of AVNIR-2, PALSAR, and polarimetric parameters for land cover classification | |
Hou et al. | Use of ALS, Airborne CIR and ALOS AVNIR-2 data for estimating tropical forest attributes in Lao PDR | |
Persson et al. | Estimating forest biomass and height using optical stereo satellite data and a DTM from laser scanning data | |
CN112446397A (zh) | 基于遥感和随机森林的产草量估算方法、装置及存储介质 | |
US20240193938A1 (en) | System and method for vegetation detection from aerial photogrammetric multispectral data | |
CN112396019A (zh) | 一种基于无人机的植被分布识别方法、系统及可读存储介质 | |
Koch | Remote sensing supporting national forest inventories NFA | |
Nkeumoe Numbisi et al. | Multi-date sentinel1 SAR image textures discriminate perennial agroforests in a tropical forest-savannah transition landscape | |
Nuthammachot et al. | Fusion of Sentinel-1A and Landsat-8 images for improving land use/land cover classification in Songkla Province, Thailand. | |
Kaushik et al. | Crop health assessment using Sentinel-1 SAR time series data in a part of central India | |
KR102476565B1 (ko) | Sar 레이더 영상의 통계적 특성을 이용한 호소의 부유성 폐기물 탐지시스템 및 그 방법 | |
Safari et al. | Integration of synthetic aperture radar and multispectral data for aboveground biomass retrieval in Zagros oak forests, Iran: An attempt on Sentinel imagery | |
Li et al. | A new method for surface water extraction using multi-temporal Landsat 8 images based on maximum entropy model | |
CN117075138A (zh) | 一种区域30米森林冠层高度遥感测算方法、系统及介质 | |
Tian et al. | A new method for estimating signal-to-noise ratio in UAV hyperspectral images based on pure pixel extraction | |
Wang et al. | [Retracted] Remote Sensing Satellite Image‐Based Monitoring of Agricultural Ecosystem | |
CN114545410B (zh) | 基于合成孔径雷达双极化数据相干性的作物倒伏监测方法 | |
Fikriyah | Mapping land cover based on time series synthetic aperture radar (SAR) data in Klaten, Indonesia | |
Guo et al. | Estimating aboveground biomass of alpine grassland during the wilting period using in situ hyperspectral, Sentinel-2 and Sentinel-1 data | |
Jackson et al. | A machine learning approach to mapping canopy gaps in an indigenous tropical submontane forest using WorldView-3 multispectral satellite imagery | |
Sexton et al. | Earth science data records of global forest cover and change | |
Martinez et al. | Multi-sensor approach to leaf area index estimation using statistical machine learning models: A case on mangrove forests |
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 |