CN115205688A - 茶树种植区域提取方法及系统 - Google Patents
茶树种植区域提取方法及系统 Download PDFInfo
- Publication number
- CN115205688A CN115205688A CN202211089743.1A CN202211089743A CN115205688A CN 115205688 A CN115205688 A CN 115205688A CN 202211089743 A CN202211089743 A CN 202211089743A CN 115205688 A CN115205688 A CN 115205688A
- Authority
- CN
- China
- Prior art keywords
- tea tree
- feature
- remote sensing
- growth
- target area
- 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
- 235000009024 Ceanothus sanguineus Nutrition 0.000 title claims abstract description 245
- 235000015459 Lycium barbarum Nutrition 0.000 title claims abstract description 245
- 238000000605 extraction Methods 0.000 title claims abstract description 42
- 241000366182 Melaleuca alternifolia Species 0.000 title 1
- 240000003553 Leptospermum scoparium Species 0.000 claims abstract description 244
- 230000012010 growth Effects 0.000 claims abstract description 75
- 238000013145 classification model Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 29
- 230000010287 polarization Effects 0.000 claims abstract description 26
- 238000012549 training Methods 0.000 claims description 25
- 238000011156 evaluation Methods 0.000 claims description 15
- 238000012795 verification Methods 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004891 communication Methods 0.000 claims description 8
- 238000001579 optical reflectometry Methods 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000007637 random forest analysis Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 8
- 244000269722 Thea sinensis Species 0.000 description 17
- 238000010586 diagram Methods 0.000 description 11
- 235000006468 Thea sinensis Nutrition 0.000 description 9
- 230000008635 plant growth Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000002310 reflectometry Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000013136 deep learning model Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000013138 pruning Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 102000006463 Talin Human genes 0.000 description 2
- 108010083809 Talin Proteins 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000006266 hibernation Effects 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 230000001932 seasonal effect Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 101100148729 Caenorhabditis elegans sar-1 gene Proteins 0.000 description 1
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000002420 orchard Substances 0.000 description 1
- 230000020477 pH reduction Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000007704 transition 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
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
-
- 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/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- 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/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Image Processing (AREA)
Abstract
本申请涉及一种茶树种植区域提取方法及系统,通过引入卫星多光谱遥感影像、SAR卫星极化遥感影像和DEM数据三种不同的遥感数据源,增加了遥感数据的多样性,通过同时获取与茶树生长相关的初级特征、二级特征和三级特征,构建了多层级遥感特征集,此外在获取特征时引入了目标区域的地理环境信息和茶树生长特性信息,充分考虑了茶树在生长过程中的生长环境,生长特性等重要因素,使得最终茶树种植区域分类模型的分类结果精度大大提升,更符合茶树的实际生长情况。
Description
技术领域
本申请涉及数据识别技术领域,特别是涉及一种茶树种植区域提取方法及系统。
背景技术
茶树是我国重要的经济树种之一,茶树的种植面积和茶叶的产量常年位居世界第一,近年来我国茶园面积和产量也出现大幅增长的现象。茶树的扩张种植能够促进区域的经济发展,但是单一的树种会降低生态系统的生物多样性,加速水土流失、土壤酸化等生态环境问题。目前,茶园面积仅有统计数据,且多为茶农自主上报统计,具有很大的不确定性,且很难在空间上获取茶园种植区域。传统的获取空间分布的方法需要测绘人员挨家挨户使用专业测绘仪器测量、核对,需要投入大量的人力和物力,时间周期长,而且成本高,因此茶园种植区的数据很难获取或缺失。
传统的茶树种植区域提取方法一般使用基于遥感数据的二分类法。这种方法一般将茶树种植区域提取转化为二分类问题,这里采用遥感影像的特征分级方法,将遥感影像的色调或颜色纳入初级特征,将遥感影像的大小,形状和纹理特征纳入二级特征。基于遥感数据的二分类法具体是先采集遥感影像,再提取遥感影像的二级特征输入至二分类器中,输出分类结果。
传统的基于遥感数据的二分类法有如下问题:
1)分类精度差,偏离实际情况。由于二分类法只能单纯的区分正样本(茶树种植区域)和负样本(非茶树种植区域),而正样本(非茶树种植区域)实际上样本组成非常复杂,正样本中包含很多与茶树特征十分相似的果园、林地等,不可避免的造成正样本数量高于实际值的现象,导致提取到的茶树种植区域的面积大于实际茶树种植区域的面积。
2)遥感影像特征选取单一。在利用遥感影像方面,基于遥感数据的二分类法计算或使用的是二级特征,没有使用更适用于茶树种植提取的更高级别的遥感图像特征,例如植被指数的纹理特征、植被指数的时序变化特征等。
3)此外,没有考虑茶树的地理环境和茶树生长特性的影响,例如某省茶树多种植于丘陵山地,平原地区也有少量分布。
发明内容
基于此,有必要针对传统茶树种植区域提取方法分类精度差,遥感影像特征选取单一且没有考虑茶树的地理环境和茶树生长特性的影响的问题,提供一种茶树种植区域提取方法及系统。
本申请提供一种茶树种植区域提取方法,所述茶树种植区域提取方法包括:
获取目标区域的遥感影像数据,构成遥感影像数据集;所述遥感影像数据包括目标区域的卫星多光谱遥感影像、目标区域的SAR卫星极化遥感影像、以及目标区域的DEM数据;
采集目标区域内的多个地物样本,构成地物样本数据集;
依据目标区域的地理环境信息和茶树生长特性信息,依据遥感影像数据集获取与茶树生长相关的初级特征、二级特征和三级特征,将所有与茶树生长相关的初级特征、二级特征和三级特征均作为待评估特征纳入特征集;
对特征集合中的每一个待评估特征均进行重要性评估,依照重要性评估结果从特征集合中筛选出至少一个优选特征纳入优选特征集;
将地物样本数据集中的至少一个地物样本纳入训练集,利用训练集和优选特征集对茶树种植区域分类模型进行训练;
将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,输出目标区域的茶树种植区域分布图。
本申请还提供一种茶树种植区域提取系统,包括:
服务器,用于执行如前述内容提及的茶树种植区域提取方法;
卫星多光谱遥感影像数据库,与所述服务器通信连接,用于存储卫星多光谱遥感影像;
SAR卫星极化遥感影像数据库,与所述服务器通信连接,用于存储SAR卫星极化遥感影像;
DEM数据库,与所述服务器通信连接,用于存储DEM数据。
本申请涉及一种茶树种植区域提取方法及系统,通过引入卫星多光谱遥感影像、SAR卫星极化遥感影像和DEM数据三种不同的遥感数据源,增加了遥感数据的多样性,通过同时获取与茶树生长相关的初级特征、二级特征和三级特征,构建了多层级遥感特征集,此外在获取特征时引入了目标区域的地理环境信息和茶树生长特性信息,充分考虑了茶树在生长过程中的生长环境,生长特性等重要因素,使得最终茶树种植区域分类模型的分类结果精度大大提升,更符合茶树的实际生长情况。
附图说明
图1为本申请一实施例提供的茶树种植区域提取方法的流程示意图。
图2为本申请一实施例提供的茶树种植区域提取系统的结构示意图。
图3为本发明一实施例提供的茶树种植区域提取方法中目标区域示意图。
图4为本发明一实施例提供的茶树种植区域提取方法中茶树种植区域地物样本分布示意图。
图5为本发明一实施例提供的茶树种植区域提取方法中重要性排名前四的遥感特征示意图。
图6为本发明一实施例提供的茶树种植区域提取方法中目标区域的茶树种植区域分布图。
附图标记:
100-服务器;200-卫星多光谱遥感影像数据库;300-SAR卫星极化遥感影像;400-DEM数据库。
具体实施方式
为了使本申请的目的.技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
一方面,本申请提供一种茶树种植区域提取方法。需要说明的是,本申请提供的茶树种植区域提取方法应用于大到市县,小到乡镇村,甚至是一个具体尺寸土地的茶树种植区域提取工作中。
此外,本申请提供的茶树种植区域提取方法不限制其执行主体。可选地,本申请提供的茶树种植区域提取方法的执行主体可以为一种茶树种植区域提取系统。具体地,本申请提供的茶树种植区域提取方法的执行主体可以为所述茶树种植区域提取系统中的服务器。
在本申请的一实施例中,所述茶树种植区域提取方法包括:
S100,获取目标区域的遥感影像数据,构成遥感影像数据集。所述遥感影像数据包括目标区域的卫星多光谱遥感影像、目标区域的SAR卫星极化遥感影像、以及目标区域的DEM数据。
具体地,在S100之前,所述茶树种植区域提取方法还包括:
选取目标区域。
目标区域可以大到市县,小到乡镇村,甚至是一个具体尺寸土地。
如图3所示,图3示出的就是Z省Q市K县的行政区示意图。星星所示的位置为县政府所在地。图3左下角为比例尺。
步骤S100中,卫星多光谱遥感影像可以采用哨兵2号卫星采集的多光谱遥感影像,其具有多个波长,是由哨兵2号卫星搭载相机拍摄得到,选择覆盖目标区域影像的云量小于20%。
SAR卫星极化遥感影像可以采用哨兵1号采集的卫星SAR极化影像,也具有多个波长,其来源是由哨兵1号卫星发射雷达波至地面的目标区域,反射回哨兵1号卫星。哨兵1号卫星自带极化功能,极化方式为VV+VH双极化方式,VV为垂直单极化方式,VH为垂直水平双极化方式。
DEM数据则可以采用SRTM-DEM-v3数据,即使用SRTM数据库的V3精度的DEM数据。
SRTM英文全称为Shuttle Radar Topography Mission,中文名称为航天飞机雷达地形测绘任务。DEM的英文全称为Digital Elevation Model,中文名称为数字高程模型。V3精度即30米精度。
S200,采集目标区域内的多个地物样本,构成地物样本数据集。
具体地,地物样本的采集方式多种多样。可选地,可以采用高精度图室内采样与野外实地采样方法结合的方式采集目标区域内的多个地物样本,作为室内地物样本。
室内采样使用计算机软件对目标区域在目标时间的历史高精度地图上采集所需要的目标地物类型的地物样本。野外实地采样采用移动设备(例如手机)在目标区域内实地采集地物样本,作为实地地物样本。
无论是室内地物样本还是实地地物样本,都需要覆盖整个目标区域。
在获得室内地物样本和实地地物样本后,依据实地地物样本检查并纠正室内地物样本,最终形成地物样本数据集。
采集目标区域内的多个地物样本后,还可以生成茶树种植区域地物样本分布示意图,如图4所示。
S300,依据目标区域的地理环境信息和茶树生长特性信息,依据遥感影像数据集获取与茶树生长相关的初级特征、二级特征和三级特征,将所有与茶树生长相关的初级特征、二级特征和三级特征均作为待评估特征纳入特征集。
具体地,与茶树生长相关的初级特征依据遥感影像数据集得到。二级特征和三级特征由初级特征计算而来,并不能依据遥感影像数据集直接得到。二级特征和三级特征首先需要经过选取,选取的过程依据目标区域的地理环境信息和茶树生长特性信息。在选取好合适的二级特征和三级特征后,再通过初级特征计算得到。
二级特征和三级特征可以包括与茶树垄状种植相关的纹理特征,与茶树季节变化相关的特征,与海拔和雷达极化相关的特征等等,充分考虑了目标区域的地理环境信息和茶树生长特性信息。
S400,对特征集合中的每一个待评估特征均进行重要性评估,依照重要性评估结果从特征集合中筛选出至少一个优选特征纳入优选特征集。
具体地,为了提供重要性评估效率,可以采取人工智能方法进行评估。例如,创建并训练机器学习模型或深度学习模型,使得机器学习模型或深度学习模型具有对一级特征、二级特征和三级特征进行重要性评估的能力,从而在一级特征、二级特征和三级特征的总量较大时,大大重要性评估的效率。
S500,将地物样本数据集中的至少一个地物样本纳入训练集,利用训练集和优选特征集对茶树种植区域分类模型进行训练。
具体地,本步骤中,创建一个茶树种植区域分类模型,所述茶树种植区域分类模型为深度学习模型。将在S100中采集的所有地物样本中的一部分样本做为训练数据,对茶树种植区域分类模型进行训练,同时在训练的过程中还同步加入优选特征集,从而使得茶树种植区域分类模型具有以下能力:在输入待测区域的遥感图像的特征集后,可以自动将地物样本进行分类,辨别哪些地物样本是茶树,哪些地物样本不是茶树,并将茶树样本注在遥感特征图像的不同位置。
S600,将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,输出目标区域的茶树种植区域分布图。
具体地,目标区域的茶树种植区域分布图如图6所示。茶树种植区域分布图中,可以看出茶树种植区域在目标区域的分布情况。
本实施例中,通过引入卫星多光谱遥感影像、SAR卫星极化遥感影像和DEM数据三种不同的遥感数据源,增加了遥感数据的多样性,通过同时获取与茶树生长相关的初级特征、二级特征和三级特征,构建了多层级遥感特征集,此外在获取特征时引入了目标区域的地理环境信息和茶树生长特性信息,充分考虑了茶树在生长过程中的生长环境,生长特性等重要因素,使得最终茶树种植区域分类模型的分类结果精度大大提升,更符合茶树的实际生长情况。
在本申请的一实施例中,当采集目标区域内的多个地物样本,构成地物样本数据集时,控制目标区域内地物样本的地物类别大于或等于6类,控制目标区域内每平方公里内采集的地物样本数量大于或等于5个。
具体地,控制目标区域内每平方公里内采集的地物样本数量大于或等于5个的目的是为了控制所有地物样本均匀覆盖整个目标区域。
在本申请的一实施例中,在所述S100之后,所述茶树种植区域提取方法还包括:
S110,对目标区域的遥感影像数据进行预处理。
具体地,对目标区域的卫星多光谱遥感影像进行预处理包括辐射定标、去云去雾、大气校正、几何校正、正射校正、镶嵌和裁剪中的一种或多种。
对目标区域的SAR卫星极化遥感影像进行预处理包括标定轨道参数、边界噪音去除、去除热噪声、辐射定标、地形校正、镶嵌和裁剪中的一种或多种。
对目标区域的DEM数据进行预处理包括镶嵌、裁剪和调整空间分辨率的操作。DEM数据采用SRTM-DEM-v3数据,即精度为30米精度,而卫星多光谱遥感影像和SAR卫星极化遥感影像都是10米精度,精度即空间分辨率,那么需要将SRTM-DEM-v3数据调整为10米精度的DEM数据。
三种遥感数据都涉及到了镶嵌,这是因为在目标时间采集的遥感数据不一定完整,即产生缺失,此时需要在目标时间的临近时间再次采集一次遥感数据,镶嵌到在目标时间采集的遥感数据中,这样目标区域的遥感数据方可采集完整。
三种遥感数据都涉及到了裁剪,和镶嵌同理,在目标时间采集的遥感数据可能覆盖面积过大,超出了目标区域的原有面积,此时就需要对遥感数据进行裁剪。
此外,在S200之后,还可以包括:S210,对目标区域内的多个地物样本进行预处理。
对目标区域内的多个地物样本进行预处理包括样本错误成分核查、样本纯净度计算、样本提纯、样本空间分布均匀化和样本数量分布均衡化中的一种或多种。
在本申请的一实施例中,所述S300包括:
S310,依据遥感影像数据集获取与茶树生长相关的初级特征。
具体地,获取到的与茶树生长相关的初级特征可以如表1所示。
表1-与茶树生长相关的初级特征表
如表1所示,包括12个波段的哨兵2号卫星多光谱遥感影像,3个波段的哨兵1号SAR极化遥感影像,以及1个波段的SRTM-DEM-v3数据。
S320,依据目标区域的地理环境信息和茶树生长特性信息,以及与茶树生长相关的初级特征,计算与茶树生长相关的二级特征。
具体地,二级特征是先选取后计算的。在选取二级特征时引入了目标区域的地理环境信息和茶树生长特性信息。
S330,依据目标区域的地理环境信息和茶树生长特性信息,以及与茶树生长相关的二级特征,计算与茶树生长相关的三级特征。
具体地,三级特征是先选取后计算的。在选取三级特征时引入了目标区域的地理环境信息和茶树生长特性信息。与茶树生长相关的二级特征和三级特征可以如表2所示。我们把二级特征和三级特征统称为高级特征。
表2-与茶树生长相关的高级特征表
如图表2所示,表2列出了三类二级特征和两类三级特征。
表2中各个特征都是使用英文简称来表示,具体释义请参见表3和表4。
表3-表2中的英文简称解释表
表4-表2中GLCM的18个特征分量的介绍表
表2中的GLCM(B8)、GLCM(VV)、GLCM(VH)和GLCM(NDTI)均各自具有18个特征分量。
本实施例中,在获取二级特征和三级特征时引入了目标区域的地理环境信息和茶树生长特性信息,充分考虑了茶树在生长过程中的生长环境,生长特性等重要因素,使得最终茶树种植区域分类模型的分类结果精度大大提升,更符合茶树的实际生长情况。
在本申请的一实施例中,所述S320包括:
S321,分析目标区域的地理环境信息和茶树生长特性信息,依据分析结果获取12月的近红外光反射率、12月的红色光反射率、5月的近红外光反射率和5月的红色光反射率。所述近红外光反射率和所述红色光反射率均为与茶树生长相关的初级特征。
S322,依据公式1分别计算12月的NDVI值和5月的NDVI值。所述NDVI值为二级特征。
其中,NDVI为NDVI值,band8为近红外光反射率,band4为红色光反射率。
具体地,本实施例具体列出了在获取二级特征时引入目标区域的地理环境信息和茶树生长特性信息的实施例。
茶树作为常绿作物,虽然全年均呈现绿色,但是其绿度在全年是随着季节变化和修剪与否密切相关的,有一定的周期性规律,这也是茶树与其他常绿乔木的区别之处。通过获取目标区域的地理环境信息知晓了目标区域归属于长江中下游区域。这个地理位置的江南茶区每年都会在春茶采摘后进行一次定型修剪,时间一般为4月下旬至5月上旬,修剪后NDVI值(即归一化差异植被指数,此处只解释一次,后续不再解释,查表2和表3可以直接得到)会降低。修剪后的茶树随着营养积累与生殖生长,经过多次生长与休止期,茶树的冠层在冬季冬眠之前达到绿度最高的时期,那么此时NDVI值最高。
因此本实施例分别采集了12月的NDVI值和5月的NDVI值,着重考虑了茶树修剪之后与冬眠之前两个时期的NDVI值植被特征及其变化,充分考虑了茶树在生长过程中的生长环境,生长特性等重要因素,使得后续计算的三级特征也兼顾这些要素。
在本申请的一实施例中,所述S330包括:
S331,分析目标区域的地理环境信息和茶树生长特性信息,依据分析结果采用公式2计算NDVI_DVI值。所述NDVI_DVI值为三级特征。
NDVI_DVI=NDVI(12)-NDVI(5) 公式2。
其中,NDVI_DVI为NDVI_DVI值,NDVI(12)为12月的NDVI值,NDVI(5)为5月的NDVI值。
具体地,上述实施例已经提及到过,三级特征是依据二级特征计算而来,那么承接上述S321至S322的实施例,本实施例,将12月的NDVI值和5月的NDVI值做差值,得到三级特征NDVI_DVI值,因此在更深层次的植被指数的差值充分体现了茶树生长过程中茶树植株绿度变化之一重要的时序信息,使得后续的茶树种植区域提取过程更符合实际情况。
在本申请的一实施例中,所述S320还包括:
S323,分析目标区域的地理环境信息和茶树生长特性信息,依据分析结果获取短波红外1反射率和短波红外2反射率。短波红外1反射率即表1中的Band11,短波红外2反射率即表1中的Band12,均为与茶树生长相关的初级特征。
S322,依据公式4计算NDTI值。所述NDTI值为二级特征。
在本申请的一实施例中,所述S330还包括:
S332,分析目标区域的地理环境信息和茶树生长特性信息,依据分析结果和上述公式4得到的NDTI计算GLCM的18个特征分量,GLCM的18个特征分量如表4所示。所述GLCM为三级特征。
具体地,由于当前绝大部分茶树均是人工按照垄状种植,山地地区和丘陵地区沿垂直于山脊或山谷的方向成片分布,种植于山地地区和丘陵地区的茶树的垄的形状一般为平行的不闭合的曲线形状。而部分种植于平原耕地上的茶树的垄形结构一般呈现平行的短直线形状。此外,还有部分茶树种植在整个小丘陵上,茶树的垄的形状呈现闭合的环状结构,与等高线分布形状极为相似,均具有十分明显且区别于其他地物的纹理特征信息。具体表现为二级特征和三级特征的纹理特征部分。
在本申请的一实施例中,在所述S400之前,所述茶树种植区域提取方法还包括:
S350,依据公式3对特征集合中的每一个待评估特征进行归一化处理;
其中,k为待评估特征的波段序号,i为行序号,j为列序号,为波段序号为k的待评估特征的第i行第j列像素的像素值,为波段序号为k的待评估特征的所有像素的平均像素值,为波段序号为k的待评估特征的所有像素的像素值标准差,为波段序号为k的待评估特征的第i行第j列像素的像素值归一化后的值。
具体地,本申请提取的波长较多,特征较多,本实施例的归一化的目的是避免各个特征的数量级差异导致出现误差。
在本申请的一实施例中,所述S400包括:
S410,采用随机森林平均精确度下降算法对特征集合中的每一个待评估特征进行重要性评估,得到每一个待评估特征的重要性指数。
具体地,所述S410包括如下步骤:
S411a,创建一个遥感特征分布模型。将地物样本数据集中的所有地物样本标注在目标区域的特征集上,生成遥感特征图像(图中未显示),将这张遥感特征图像作为标准图像C,每一个地物样本被标注后生成一个标签。
S411b,将特征集和所有标签作为训练数据输入至遥感特征分布模型,对遥感特征分布模型进行训练,使得遥感特征分布模型具有自动在特征集上标注的功能。
S411c,将特征集和所有标签输入至训练后的遥感特征分布模型,运行训练后的遥感特征分布模型,获取训练后的遥感特征分布模型输出的初始预测图像C’。
S412a,在特征集中,随机抽出一个待评估特征置入待评估特征队列的随机一个位置,以改变待评估特征的排列顺序,将改变待评估特征的排列顺序的特征集合和所有标签输入至训练后的遥感特征分布模型,运行训练后的遥感特征分布模型,获取训练后的遥感特征分布模型输出的第一预测图像C1’。
S412c,在特征集中,再抽出S412a中同一个待评估特征置入待评估特征队列的随机一个位置,但是需要注意的是随机位置不可以是此待评估特征处于待评估特征队列的初始位置和S412a放置的位置,避免重复,这样又改变了一次待评估特征的排列顺序,将改变待评估特征的排列顺序的特征集输入至训练后的遥感特征分布模型,运行训练后的遥感特征分布模型,获取训练后的遥感特征分布模型输出的第二预测图像C2’。
S413,那么通过第二次抽出待评估特征改变待评估特征的排列顺序后,此待评估特征预测正确率的精度就可以通过公式8计算。
S414,那么我们后续继续依照上述的原理继续对其余每个待评估特征进行抽出并置入待评估特征队列的随机位置,改变待评估特征的排列顺序,计算预测正确率,进行再次抽出待评估特征置入待评估特征队列另一随机位置,改变待评估特征的排列顺序,计算预测正确率,从而计算此待评估特征预测正确率的精度,重复直至将所有待评估特征均计算出预测正确率的精度J。
正确率的精度即正确率的变化率,体现了抽出一个待评估特征后改变待评估特征的排列顺序后,预测结果的正确率的波动情况。预测正确率的精度越高,代表这个待评估特征的抽出对预测结果的正确率的波动越大,那么代表这个待评估特征的重要性就越高。因此,采用随机森林平均精确度下降算法评估各个待评估特征的重要性,实际上就是得出各个待评估特征对应的预测正确率的精度,我们将预测正确率的精度作为重要性指数。
S420,将重要性指数按从大到小的顺序进行排序。
具体地,预测正确率的精度越高,重要性指数越高,代表重要性越高。
S430,选取重要性指数前N%对应的待评估特征作为优选特征。
具体地,例如,N可以取大于等于20且小于等于40的数值范围中的一个数字。可选地,N为40。
S440,创建优选特征集,将所有优选特征纳入所述优选特征集。
具体地,例如我们通过S430选取了前40%的优选特征,其中重要性指数前六的待评估特征为二级特征NDVI_RE2(归一化差异红边2植被指数,见表2和表3),二级特征GLCM(VV)-savg(依据一级特征VV计算的GLCM纹理特征的savg波段,见表2),一级特征Band2(蓝色光反射率,见表1),一级特征ELevation(地面高程,见表1),三级特征NDVI_DVI(见表2),一级特征Band12(短波红外光2反射率,见表1)。那么我们将前40%的优选特征作为优选特征纳入所述优选特征集。
可选地,我们可以生成重要性排名前四的遥感特征示意图,如图5所示。图5包括a,b,c,d四张图,图5中的a对应二级特征NDVI_RE2,图5中的b对应二级特征GLCM(VV)-savg,图5中的c对应一级特征Band2,图5中的d对应一级特征ELevation。图5中的a中的左下角的图示为比例尺。图5中的a中的右上角的图示为指北方向标,图5中的a中的右下角的渐变条的含义是:颜色越浅,代表特征值越大。颜色越深,代表特征值越小。渐变条的最浅色端代表特征值的最大值,渐变条的最深色端代表特征值的最小值。图5中的b的图示,图5中的c的图示和图5中的d的图示和图5中的a的图示相同,此处不再赘述。
在本申请的一实施例中,所述S500包括:
S510,将地物样本数据集内占地物样本总量X%的地物样本纳入训练集。
S520,利用训练集和优选特征集对茶树种植区域分类模型进行训练。
S530,将地物样本数据集内占地物样本总量Y%的地物样本纳入验证集。
S540,将验证集输入至茶树种植区域分类模型,运行茶树种植区域分类模型得出分类结果。所述分类结果包括目标区域内归属于茶树地物类别的地物样本数量,以及归属于非茶树地物类别的地物样本数量。
S550,将分类结果和验证集的茶树分类实际情况进行比对,得到茶树分类精度。
S560,判断茶树分类精度是否大于或等于预设茶树分类精度。
S570,若茶树分类精度大于或等于预设茶树分类精度,则执行后续将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,输出目标区域的茶树种植区域分布图的步骤。
S580,若茶树分类精度小于预设茶树分类精度,则将地物样本数据集内占地物样本总量Z%的地物样本纳入测试集。X,Y和Z的和为100,且X大于Y且X大于Z。
S590,将测试集输入茶树种植区域分类模型,利用测试集对预设茶树分类精度进行调整,返回所述将验证集输入至茶树种植区域分类模型,运行茶树种植区域分类模型得出分类结果的步骤。
具体地,本实施例是将地物样本数据集内占地物样本总量按X:Y:Z的比例分配,X%作为训练集,Y%作为验证集,Z%作为测试集,三个不同的集合来分别完成不同的功能。
X:Y:Z可以为6:2:2,即X=60,Y=20,Z=20。X:Y:Z也可以为8:1:1,即X=80,Y=10,Z=10。X:Y:Z也可以为7:3:1,即X=70,Y=30,Z=10。训练集占比是最高的,X大于Y且X大于Z,可选地,X还要大于或等于60%,这样的目的是避免训练出来的茶树种植区域分类模型的泛化能力特别弱,避免茶树种植区域分类模型无法用于非训练集的预测与分类。
训练集用于训练茶树种植区域分类模型。验证集用于验证茶树种植区域分类模型的分类结果是否正确,从而可以计算出茶树种植区域分类模型的分类精度,若茶树分类精度大于或等于预设茶树分类精度,则茶树种植区域分类模型是合格的,可以执行后续步骤投入使用。若茶树分类精度小于预设茶树分类精度,则茶树种植区域分类模型是不合格的,需要使用测试集去调整预设茶树分类精度的值。预设茶树分类精度可以为95%。
S540至S550通过验证集对茶树种植区域分类模型的分类结果进行评价的过程中,使用了混淆矩阵的原理,例如,目标区域的茶树分布的正确状态参见验证集的实际分布状态,例如实际为50个茶树地物样本,50个非茶树地物样本。而样验证集输入至茶树种植区域分类模型,运行茶树种植区域分类模型,茶树种植区域分类模型输出的分类结果是45个茶树地物样本,55个非茶树地物样本。
将原来的50个茶树地物样本中的30个茶树地物样本成功分类出来,但是20个茶树地物样本错认为非茶树地物样本。
此外茶树种植区域分类模型输出的分类结果是将原来的50个非茶树地物样本中的35个认定为非茶树地物样本,其中的15个认定为茶树样本。
那么茶树分类精度就是30除以50乘以100%为60%。
在本申请的一实施例中,还可以计算下整体分类精度,计算整体分类精度的结果就是[(30+35)]/100=75%。
在本申请的一实施例中,S560中,若整体分类精度大于等于预设整体分类精度阈值且茶树分类精度高于预设茶树分类精度,才执行S570,否则执行S580。
可选地,预设茶树分类精度为95%,预设整体分类精度为90%。
测试集用于调整预设茶树分类精度的值,找到一个适合茶树种植区域分类模型当前分类精度(换言之是符合茶树种植区域分类模型的分类能力的预设茶树分类精度)。
S600中,将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,茶树种植区域分类模型可以输出(0,1)二值图像,1表示茶树,0表示非茶树,得到最终精确的茶树提取结果,可以具象化为目标区域的茶树种植区域分布图输出,如图6所示。
另一方面,本申请提供一种茶树种植区域提取系统。
如图2所示,在本申请的一实施例中,所述茶树种植区域提取系统包括服务器100、卫星多光谱遥感影像数据库200、SAR卫星极化遥感影像300和DEM数据库400。
所述服务器100用于执行如前述任意一个实施例所提及的茶树种植区域提取方法。所述卫星多光谱遥感影像数据库200与所述服务器100通信连接。所述卫星多光谱遥感影像数据库200用于存储卫星多光谱遥感影像。
SAR卫星极化遥感影像数据库300与所述服务器100通信连接。SAR卫星极化遥感影像数据库300用于存储SAR卫星极化遥感影像。
所述DEM数据库400与所述服务器100通信连接。所述DEM数据库400用于存储DEM数据。
具体地,为了行文简洁,本实施例中提及的服务器100、卫星多光谱遥感影像数据库、SAR卫星极化遥感影像和DEM数据库仅在本实施例中进行标号,前述茶树种植区域提取方法中出现相同名称的设备或数据库不进行标号。
以上所述实施例的各技术特征可以进行任意的组合,各方法步骤也并不做执行顺序的限制,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本申请专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利要求为准。
Claims (10)
1.一种茶树种植区域提取方法,其特征在于,所述茶树种植区域提取方法包括:
获取目标区域的遥感影像数据,构成遥感影像数据集;所述遥感影像数据包括目标区域的卫星多光谱遥感影像、目标区域的SAR卫星极化遥感影像、以及目标区域的DEM数据;
采集目标区域内的多个地物样本,构成地物样本数据集;
依据目标区域的地理环境信息和茶树生长特性信息,依据遥感影像数据集获取与茶树生长相关的初级特征、二级特征和三级特征,将所有与茶树生长相关的初级特征、二级特征和三级特征均作为待评估特征纳入特征集;
对特征集合中的每一个待评估特征均进行重要性评估,依照重要性评估结果从特征集合中筛选出至少一个优选特征纳入优选特征集;
将地物样本数据集中的至少一个地物样本纳入训练集,利用训练集和优选特征集对茶树种植区域分类模型进行训练;
将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,输出目标区域的茶树种植区域分布图。
2.根据权利要求1所述的茶树种植区域提取方法,其特征在于,当采集目标区域内的多个地物样本,构成地物样本数据集时,控制目标区域内地物样本的地物类别大于或等于6类,控制目标区域内每平方公里内采集的地物样本数量大于或等于5个。
3.根据权利要求1所述的茶树种植区域提取方法,其特征在于,在所述获取目标区域的遥感影像数据,构成遥感影像数据集之后,所述茶树种植区域提取方法还包括:
对目标区域的遥感影像数据进行预处理。
4.根据权利要求1所述的茶树种植区域提取方法,其特征在于,所述依据目标区域的地理环境信息和茶树生长特性信息,依据遥感影像数据集获取与茶树生长相关的初级特征、二级特征和三级特征,将所有与茶树生长相关的初级特征、二级特征和三级特征均作为待评估特征纳入特征集,包括:
依据遥感影像数据集获取与茶树生长相关的初级特征;
依据目标区域的地理环境信息和茶树生长特性信息,以及与茶树生长相关的初级特征,计算与茶树生长相关的二级特征;
依据目标区域的地理环境信息和茶树生长特性信息,以及与茶树生长相关的二级特征,计算与茶树生长相关的三级特征。
6.根据权利要求5所述的茶树种植区域提取方法,其特征在于,所述依据目标区域的地理环境信息和茶树生长特性信息,以及与茶树生长相关的二级特征,计算与茶树生长相关的三级特征,包括:
分析目标区域的地理环境信息和茶树生长特性信息,依据分析结果采用公式2计算NDVI_DVI值;所述NDVI_DVI值为三级特征;
NDVI_DVI=NDVI(12)-NDVI(5) 公式2;
其中,NDVI_DVI为NDVI_DVI值,NDVI(12)为12月的NDVI值,NDVI(5)为5月的NDVI值。
8.根据权利要求1所述的茶树种植区域提取方法,其特征在于,所述对特征集合中的每一个待评估特征均进行重要性评估,依照重要性评估结果从特征集合中筛选出至少一个优选特征纳入优选特征集,包括:
采用随机森林平均精确度下降算法对特征集合中的每一个待评估特征进行重要性评估,得到每一个待评估特征的重要性指数;
将重要性指数按从大到小的顺序进行排序;
选取重要性指数前N%对应的待评估特征作为优选特征;
创建优选特征集,将所有优选特征纳入所述优选特征集。
9.根据权利要求1-8中任意一项所述的茶树种植区域提取方法,其特征在于,所述将地物样本数据集中的至少一个地物样本纳入训练集,利用训练集和优选特征集对茶树种植区域分类模型进行训练,包括:
将地物样本数据集内占地物样本总量X%的地物样本纳入训练集;
利用训练集和优选特征集对茶树种植区域分类模型进行训练;
将地物样本数据集内占地物样本总量Y%的地物样本纳入验证集;
将验证集输入至茶树种植区域分类模型,运行茶树种植区域分类模型得出分类结果;所述分类结果包括目标区域内归属于茶树地物类别的地物样本数量,以及归属于非茶树地物类别的地物样本数量;
将分类结果和验证集的茶树分类实际情况进行比对,得到茶树分类精度;
判断茶树分类精度是否大于或等于预设茶树分类精度;
若茶树分类精度大于或等于预设茶树分类精度,则执行后续将目标区域的地物样本数据集和优选特征集输入至茶树种植区域分类模型,输出目标区域的茶树种植区域分布图的步骤;
若茶树分类精度小于预设茶树分类精度,则将地物样本数据集内占地物样本总量Z%的地物样本纳入测试集;X,Y和Z的和为100,且X大于Y且X大于Z;
将测试集输入茶树种植区域分类模型,利用测试集对预设茶树分类精度进行调整,返回所述将验证集输入至茶树种植区域分类模型,运行茶树种植区域分类模型得出分类结果的步骤。
10.一种茶树种植区域提取系统,其特征在于,包括:
服务器,用于执行如权利要求1-9中任意一项所述的茶树种植区域提取方法;
卫星多光谱遥感影像数据库,与所述服务器通信连接,用于存储卫星多光谱遥感影像;
SAR卫星极化遥感影像数据库,与所述服务器通信连接,用于存储SAR卫星极化遥感影像;
DEM数据库,与所述服务器通信连接,用于存储DEM数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211089743.1A CN115205688A (zh) | 2022-09-07 | 2022-09-07 | 茶树种植区域提取方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211089743.1A CN115205688A (zh) | 2022-09-07 | 2022-09-07 | 茶树种植区域提取方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115205688A true CN115205688A (zh) | 2022-10-18 |
Family
ID=83573175
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211089743.1A Pending CN115205688A (zh) | 2022-09-07 | 2022-09-07 | 茶树种植区域提取方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115205688A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117456378A (zh) * | 2023-12-20 | 2024-01-26 | 山东锋士信息技术有限公司 | 基于卫星遥感的水利数字孪生底座要素实现方法和系统 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108898070A (zh) * | 2018-06-07 | 2018-11-27 | 湖南中科星图信息技术有限公司 | 一种基于无人机平台的高光谱遥感提取薇甘菊装置及方法 |
CN108960300A (zh) * | 2018-06-20 | 2018-12-07 | 北京工业大学 | 一种基于深度神经网络的城市土地利用信息分析方法 |
CN111242224A (zh) * | 2020-01-16 | 2020-06-05 | 贵州省草业研究所 | 一种基于无人机提取分类样本点的多源遥感数据分类方法 |
AU2020101054A4 (en) * | 2020-06-19 | 2020-07-30 | Guizhou Institute Of Pratacultural | A Multi-source Remote Sensing Data Classification Method Based On the Classification Sample Points Extracted By the UAV |
CN111881953A (zh) * | 2020-07-14 | 2020-11-03 | 安徽大学 | 基于局部二值模式和knn分类器的遥感高光谱图像分类方法 |
CN112949612A (zh) * | 2021-04-22 | 2021-06-11 | 辽宁工程技术大学 | 一种基于无人机高分辨率遥感影像海岸带地物分类方法 |
CN113435254A (zh) * | 2021-05-27 | 2021-09-24 | 云南师范大学 | 一种基于哨兵二号影像的耕地深度学习提取方法 |
CN113657158A (zh) * | 2021-07-13 | 2021-11-16 | 西安电子科技大学 | 基于Google Earth Engine的大尺度大豆种植区域提取算法 |
CN114219795A (zh) * | 2021-12-17 | 2022-03-22 | 青岛农业大学 | 一种基于高光谱成像系统的茶树干旱诱导成分及干旱程度评估的预测方法和系统 |
CN114694038A (zh) * | 2022-04-11 | 2022-07-01 | 中国矿业大学 | 基于深度学习的高分辨率遥感影像分类方法及系统 |
CN114926748A (zh) * | 2022-06-15 | 2022-08-19 | 安徽理工大学 | Sentinel-1/2微波与光学多光谱影像结合的大豆遥感识别方法 |
-
2022
- 2022-09-07 CN CN202211089743.1A patent/CN115205688A/zh active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108898070A (zh) * | 2018-06-07 | 2018-11-27 | 湖南中科星图信息技术有限公司 | 一种基于无人机平台的高光谱遥感提取薇甘菊装置及方法 |
CN108960300A (zh) * | 2018-06-20 | 2018-12-07 | 北京工业大学 | 一种基于深度神经网络的城市土地利用信息分析方法 |
CN111242224A (zh) * | 2020-01-16 | 2020-06-05 | 贵州省草业研究所 | 一种基于无人机提取分类样本点的多源遥感数据分类方法 |
AU2020101054A4 (en) * | 2020-06-19 | 2020-07-30 | Guizhou Institute Of Pratacultural | A Multi-source Remote Sensing Data Classification Method Based On the Classification Sample Points Extracted By the UAV |
CN111881953A (zh) * | 2020-07-14 | 2020-11-03 | 安徽大学 | 基于局部二值模式和knn分类器的遥感高光谱图像分类方法 |
CN112949612A (zh) * | 2021-04-22 | 2021-06-11 | 辽宁工程技术大学 | 一种基于无人机高分辨率遥感影像海岸带地物分类方法 |
CN113435254A (zh) * | 2021-05-27 | 2021-09-24 | 云南师范大学 | 一种基于哨兵二号影像的耕地深度学习提取方法 |
CN113657158A (zh) * | 2021-07-13 | 2021-11-16 | 西安电子科技大学 | 基于Google Earth Engine的大尺度大豆种植区域提取算法 |
CN114219795A (zh) * | 2021-12-17 | 2022-03-22 | 青岛农业大学 | 一种基于高光谱成像系统的茶树干旱诱导成分及干旱程度评估的预测方法和系统 |
CN114694038A (zh) * | 2022-04-11 | 2022-07-01 | 中国矿业大学 | 基于深度学习的高分辨率遥感影像分类方法及系统 |
CN114926748A (zh) * | 2022-06-15 | 2022-08-19 | 安徽理工大学 | Sentinel-1/2微波与光学多光谱影像结合的大豆遥感识别方法 |
Non-Patent Citations (6)
Title |
---|
YONG XIE 等: "Multielement Classification of a Small Fragmented Planting Farm Using Hyperspectral Unmanned Aerial Vehicle Image", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
刘雪萍 等: "基于多特征优选的平和县蜜柚园遥感提取与扩张分析", 《贵州大学学报(自然科学版)》 * |
徐伟燕 等: "基于资源三号卫星影像的茶树种植区提取", 《农业工程学报》 * |
熊皓丽 等: "基于GEE云平台的福建省10m分辨率茶园专题空间分布制图", 《地球信息科学学报》 * |
秦泉等: "面向对象的GF-1卫星影像苹果树种植面积遥感提取研究――以山地丘陵地区的栖霞市为例", 《沙漠与绿洲气象》 * |
赵晓晴 等: "时序Sentinel-2A影像光谱特征的茶园提取应用", 《测绘科学》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117456378A (zh) * | 2023-12-20 | 2024-01-26 | 山东锋士信息技术有限公司 | 基于卫星遥感的水利数字孪生底座要素实现方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108764255B (zh) | 一种冬小麦种植信息的提取方法 | |
Tomppo et al. | The multi-source national forest inventory of Finland-methods and results 2011 | |
CN107527014A (zh) | 县级作物种植面积遥感统计抽样调查方案设计方法 | |
CN106918816B (zh) | 农作物长势监测系统及其数据处理方法和模块 | |
CN113205014B (zh) | 一种基于图像锐化的时序数据耕地提取方法 | |
CN111798327A (zh) | 基于高光谱影像的小麦产量计算模型的构建方法及其应用 | |
CN108458978B (zh) | 基于敏感波段和波段组合最优的树种多光谱遥感识别方法 | |
CN113033453A (zh) | 一种适用于景观破碎区作物类型遥感识别的方法及系统 | |
CN114120132A (zh) | 一种气象和红边波段遥感结合的农作物估产方法和装置 | |
CN116543316B (zh) | 一种利用多时相高分辨率卫星影像识别稻田内草皮的方法 | |
Kang et al. | Support vector machine classification of crop lands using sentinel-2 imagery | |
CN112861435A (zh) | 一种红树林质量遥感反演方法及智能终端 | |
CN112036313A (zh) | 烟草种植面积检测方法、装置、设备及可读存储介质 | |
Koranteng et al. | Remote sensing study of land use/cover change in West Africa | |
CN115205688A (zh) | 茶树种植区域提取方法及系统 | |
Liu et al. | Evaluating how lodging affects maize yield estimation based on UAV observations | |
Aguirre-Salado et al. | Modelling site selection for tree plantation establishment under different decision scenarios | |
Erdenee et al. | Particular agricultural land cover classification case study of Tsagaannuur, Mongolia | |
WO2023131949A1 (en) | A versatile crop yield estimator | |
CN115115948B (zh) | 一种基于随机森林和辅助因子的林地信息精细化提取方法 | |
CN115861836A (zh) | 基于多时相Sentinel-2数据的大豆种植区提取方法 | |
CN115861793A (zh) | 一种基于最小累积阻力模型的区域生态安全格局构建方法 | |
CN115527108A (zh) | 一种基于多时相Sentinel-2快速识别水土流失人为扰动地块的方法 | |
CN114997730A (zh) | 一种基于多维特征的城乡规划设计区域数据智能监测分析评价系统 | |
CN106682756B (zh) | 一种基于rs/gis的冬虫夏草产量预测模型 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20221018 |
|
RJ01 | Rejection of invention patent application after publication |