CN113642579A - 南荻生长区域确定、分割方法及系统 - Google Patents
南荻生长区域确定、分割方法及系统 Download PDFInfo
- Publication number
- CN113642579A CN113642579A CN202110772118.6A CN202110772118A CN113642579A CN 113642579 A CN113642579 A CN 113642579A CN 202110772118 A CN202110772118 A CN 202110772118A CN 113642579 A CN113642579 A CN 113642579A
- Authority
- CN
- China
- Prior art keywords
- area
- triarrhena
- vegetation
- sacchariflora
- spectral feature
- 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
- 241001326934 Triarrhena Species 0.000 title claims abstract description 160
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000003595 spectral effect Effects 0.000 claims abstract description 110
- 230000012010 growth Effects 0.000 claims abstract description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 27
- 230000011218 segmentation Effects 0.000 claims description 57
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 31
- 230000008859 change Effects 0.000 claims description 11
- 241000219000 Populus Species 0.000 claims description 10
- 241001290610 Abildgaardia Species 0.000 claims description 8
- 241000878007 Miscanthus Species 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 235000007164 Oryza sativa Nutrition 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 235000009566 rice Nutrition 0.000 claims description 5
- 241000226665 Anaphalis Species 0.000 claims 4
- 240000007594 Oryza sativa Species 0.000 claims 1
- 238000000605 extraction Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 239000002028 Biomass Substances 0.000 description 4
- 241000209094 Oryza Species 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003709 image segmentation Methods 0.000 description 3
- 244000025254 Cannabis sativa Species 0.000 description 2
- 244000273256 Phragmites communis Species 0.000 description 2
- 235000014676 Phragmites communis Nutrition 0.000 description 2
- 238000003066 decision tree Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241000722731 Carex Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229920002678 cellulose Polymers 0.000 description 1
- 239000001913 cellulose Substances 0.000 description 1
- 235000019646 color tone Nutrition 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000243 photosynthetic effect Effects 0.000 description 1
- 230000008636 plant growth process Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
Abstract
本发明公开了南荻生长区域确定、分割方法及系统，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。相比现有技术，本发明能利用南荻在不同生长阶段的光谱特征来区分待监测区域的南荻与其他植被，从而准确获取到待监测区域南荻的分布情况。
Description
技术领域
本发明涉及南荻生长监控技术领域，尤其涉及南荻生长区域确定、分割方法及系统。
背景技术
南荻是我国特有的芒属植物，其具有生物质产量高，光合效率强，纤维素含量高和耐贫瘠等诸多优点，是一种极具开发潜能的生态修复植物和生物质原料，近年来随着国人对环境保护越来越重视，南荻这一生态友好型植物也越来越受到研究者的青睐。
发明内容
本发明提供了南荻生长区域确定、分割方法及系统，用于解决现有技术无法准确获取南荻分布情况的技术问题。
为解决上述技术问题，本发明提出的技术方案为：
一种南荻生长区域确定方法，包括以下步骤：
根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；
获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。
优选的，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相，包括以下步骤：
分别从历史数据中获取南荻以及其他植被在不同时期的原始光谱特征数据，构建南荻以及其他植被生长过程的时序光谱特征序列；
对于其他植被的每一种植被A的时序光谱特征序列，均执行以下步骤：
分别计算植被A的时序光谱特征序列内不同时期的光谱特征数据与其对应时期的南荻的光谱特征数据之间的差值，选取植被A与南荻光谱特征数据之间的差值最大的时期或植被A 与南荻光谱特征数据之间的差值大于预设差阈值的时期作为植被A与南荻的区分时期；
基于南荻与其他植被中每一种植被的区分时期构建南荻在待监测区域内的最佳监测时相。
优选的，构建南荻以及其他植被生长过程的时序光谱特征序列，包括以下步骤：
基于南荻或其他植被在不同时期的原始光谱特征数据拟合南荻或其他植被的时序光谱特征曲线，判断时序光谱特征曲线上各个部分的光谱特征数据时序变化是否在对应的正常变化范围内，若存在任意部分的光谱特征数据的时序变化不在对应的正常变化范围内，则判断部分对应原始光谱特征数据有噪声干扰，修正部分的对应的原始光谱特征数据，使修正后的光谱特征数据拟合的时序光谱特征曲线上各个部分的光谱特征数据时序变化均在其对应的正常范围内；再使用修正后光谱特征数据构建时序光谱特征序列。
优选的，基于南荻与其他植被中每一种植被的区分时期构建南荻在待监测区域内的最佳监测时相，具体包括以下步骤：
若存在任意两种以上植被的区分时期存在重合，选取重合的区分时期作为两种以上植被的最优区分时期；
对于区分时期与其他植被不存在重合的植被，选取植被的区分时期中差值最大的时期作为最优时期；
基于其他植被中每一种植被的最优区分时期构建南荻在待监测区域内的最佳监测时相。
优选的，光谱特征数据包括：归一化植被指数NDVI、归一化水体指数NDWI、绿度Ratio Green中一种或任意几种的组合。
优选的，当待监测区域为大陆性亚热带季风湿润气候，且其他植被包括水稻、苔草、杨树、油菜时，光谱特征数据包括归一化植被指数NDVI、归一化水体指数NDWI、绿度Ratio Green 时，南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值如下：
优选的，还包括以下步骤：
获取待监测区域内各个子区域的植被的形状指数、长宽比指数；
将待监测区域内各个子区域的植被的形状指数分别与预设的南荻的形状阈值进行比较，将待监测区域内各个子区域植被的长宽比指数分别与预设的南荻的长宽比阈值进行比较，当存在任意子区域的光谱特征数据在南荻的光谱特征阈值内，子区域植被的形状指数在南荻的形状阈值内，且子区域植被的长宽比指数在南荻的长宽比阈值内，则判断子区域为南荻的生长区域。
优选的，形状阈值为＜10.2，长宽比阈值为＜8.7。
一种南荻生长区域分割方法，包括以下步骤：
确定待分割区域的水体区域和非水体区域，采用第一分割尺度分割待监测区域的水体区域和非水体区域；
确定非水体区域上的道路、建筑区域和植被区域，采用第二分割尺度从非水体区域内分割出的道路、建筑区域，得到植被区域；
在植被区域内，采用上述任意一项的南荻生长区域确定方法确定植被区域内的南荻生长区域，并采用第三分割尺度分割南荻的生长区域；其中，第一分割尺度大于第二分割尺度，且第二分割尺度大于第三分割尺度；第一分割尺度为200，第二分割尺度为150；第三分割尺度为50。
一种计算机系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述算法的步骤。
本发明具有以下有益效果：
1、本发明中的南荻生长区域确定、分割方法及系统，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。相比现有技术，本发明能利用南荻在不同生长阶段的光谱特征来区分待监测区域的南荻与其他植被，从而准确获取到待监测区域南荻的分布情况。
2、在优选方案中，本发明中通过设置分类参数归一化植被指数(NDVI)、归一化水体指数(NDWI)、绿度(Ratio Green)、亮度值(Brightness)、形状指数(Shape Index)、长宽比(Length/width)等的阈值来区分南荻与其他植被，能进一步提高待监测区域南荻的分布情况的准确性。
3、在优选方案中，本技术方案在从待监测区域的卫星图像中分割南荻的生长区域时，针对不同的分割对象，采用不同的分割尺度，在保证分割的精确性的同时，能有效提高分割的速度。
除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照附图，对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是本发明优选实施例中的洞庭湖典型植被NDVI时间序列曲线；
图2是本发明优选实施例中的不同分割尺度分割效果图，其中，(a)为分割尺度50的效果图，(b)为分割尺度150的效果图，(c)为分割尺度200的效果图；
图3是本发明优选实施例中的关键时相卫星影像和NDVI数据图，其中，(a)为4月16日的卫星影像，(b)为4月16日的NDVI数据图；(c)为8月4日的卫星影像；(d) 8月4日的NDVI数据图；(e)为12月16日的卫星影像；(f)为12月16日的NDVI数据图；(g)为NDVI组合影像图；
图4是本发明优选实施例中的南荻空间分布提取技术路线图；
图5是本发明优选实施例中的东洞庭湖实地2015-2019年南荻空间分布图，其中，(a) -(e)分别为2015-2019年的南荻空间分布图；
图6为本发明中的南荻生长区域确定方法的流程图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
实施例一：
如图6所示，本实施中公开了一种南荻生长区域确定方法，包括以下步骤：
根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；
获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。
此外，在本实施例中，还公开了一种南荻生长区域分割方法，包括以下步骤：
确定待分割区域的水体区域和非水体区域，采用第一分割尺度分割待监测区域的水体区域和非水体区域；
确定非水体区域上的道路、建筑区域和植被区域，采用第二分割尺度从非水体区域内分割出的道路、建筑区域，得到植被区域；
在植被区域内，采用上述的南荻生长区域确定方法确定植被区域内的南荻生长区域，并采用第三分割尺度分割南荻的生长区域；其中，第一分割尺度大于第二分割尺度，且第二分割尺度大于第三分割尺度。
此外，在本实施例中，还公开了一种计算机系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述算法的步骤。
本发明中的南荻生长区域确定、分割方法及系统，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。相比现有技术，本发明能利用南荻在不同生长阶段的光谱特征来区分待监测区域的南荻与其他植被，从而准确获取到待监测区域南荻的分布情况。
实施例二：
实施例二是实施例一的优选实施例，其与实施例一的不同之处在于，对本发明中的南荻生长区域确定、分割方法的具体步骤进行了细化：
在本实施例中，公开了一种南荻生长区域确定、分割方法，该方法综合时间序列数据重构、面向对象特征提取、光谱特征变量提取、阈值选取和决策树构建等技术，采用综合时空谱多维遥感特征的精细分类方法，提取南荻空间分布。
下面以洞庭湖作为待监测区域，来介绍南荻生长区域确定、分割方法的具体实施步骤：
一、时间序列数据重构
本次研究选用GF-1宽幅数据、Landsat-8、Sentinel-2等多源中高分辨率卫星影像构建洞庭湖湿地植被全生育期月度NDVI时间序列数据集，对比洞庭湖湿地主要植被时间序列曲线，确定南荻最佳监测时相。NDVI时间序列重建的实质就是降低序列本身的噪声水平，重建高质量的数据集。各种重建方法都基于以下两个前提假设：①NDVI的时序变化对应于植被的生长过程；②在NDVI时序数据中出现的陡升或陡降都认为是与植物生长过程不一致的噪声干扰。
为了寻找南荻与其他植被生长差异最大时期，确定南荻最佳监测时相，本研究将全年12 个月的卫星影像数据，分别计算NDVI，并按照时间顺序将所得的12个月NDVI数据组合，获得时间序列数据。洞庭湖湿地几种植被NDVI时间序列曲线见图1。
从图1中的时间序列曲线可以看出南荻在5月与水稻、油菜的NDVI差异最大，可将南荻与水稻和油菜区分；在8月，南荻与苔草NDVI差异最大，可将南荻与苔草区分；在4～11月，南荻与杨树NDVI时间序列曲线接近，但在1月、2月、12月，南荻NDVI值均低于杨树，因而选取该时间段的NDVI数据可以区分南荻和杨树。因此，本研究选取当年5月、8 月、12月(或1月)各一期卫星影像，作为南荻提取关键时相数据，其中，选取具体情况如表1所示，在5月、8月没有有效数据时，分别考虑4月、9月数据作为关键时相补充数据。
表1南荻遥感监测时相表
本阶段研究所用遥感数据详情见表2：
表2选用卫星影像数据详情表
二、面向对象特征提取
采用面向对象分类方法，选择合理分割尺度、分类特征和分类规则，提取南荻空间分布并统计面积。
①多尺度影像分割
遥感影像的多尺度分割是面向对象分类方法的基本方法与手段，分割的最终目的是将具有共同光谱特征和属性的临近像元归类划分为整体性的小区域(区域的大小由分割尺度及临近像元相似性而定)，从而为进一步细分提供信息构建基础。有效的分割尺度不仅可减少计算机运算冗沉、节省工作时间，而且可以更有效地结合对象的空间(如对象形状，长宽等) 与光谱属性差异。
不同的分割尺度产生的对象多边形的数量、大小及提取精度会有所不同，同时不同尺度提取的信息也能反映出不同层次地物的结构差异，因此针对不同类型的地物选择不同的尺度进行分割。洞庭湖湿地景观类型中南荻群落浓疏相间，同时还有水体、岛屿、道路、房屋、南荻垛，图像分割时，兼顾分类效率和精度，结合洞庭湖湿地特点，经选择不同分割尺度进行分割试验，最终选择50、150、200为最有分割尺度。当分割尺度为200时，主要依据色调将影像分为大范围混合对象，可将水体与非水体区分开；当分割尺度为150时，道路有较完整的多边形，能将植被与道路明显区分；当分割尺度为50时，影像分割图斑破碎，但芦苇和杨树、苔草等分布边界较为清晰完整。分割形状因子设置为0.2，紧致度设置为0.5。
②分类层次建立
分类层次体系是指在不同尺度分割后，不同尺度层的对象间的纵向关系和相同尺度分割产生的同层中对象间的横向关系。利用分割产生的各层次地类的特征信息，建立分类规则集，并分析其上下层关系，使之具有继承关系，将分类信息在层次间传递，并在最后一层继承所有地物。
本研究分三层设置不同分割尺度，最终获取不同地物空间分布和南荻空间分布信息。 Level 1为第一层，分割尺度为200，研究区水体、非水体在Level 1中进行提取；Level2在 Level 1的基础上分割，尺度为150，继续进行分割，主要提取中尺度的地物(道路、建筑等)；Level 3在Level 2的基础上以尺度为50继续分割，主要提取相对较小斑块的地物对象(芦苇、苔草、杨树等)，其中，研究区(待监测区域)的主要地物特征如表3所示。分割效果图如图2所示。
表3研究区主要地物特征(假彩色影像中)
③分类参数及阈值选取
NDVI反映植被生长情况，NDVI值介于-1～1之间，植被长势越好、生物量越大，NDVI值越大，而无植被覆盖区域NDVI值小。图3中水体NDVI值最小，在NDVI图上显示为棕黄色，而植被在NDVI图上显示为绿色或蓝色，且植被长势越好、生物量越大，图中颜色越趋于蓝色。图3(b)中可以看出，在4月16日，南荻和苔草比洞庭湖外的耕地NDVI值大；图3(d)中可以看出，在8月4日，南荻比周边的苔草NDVI值大；图3(f)中可以看出，在12月16日，草和杨树比南荻NDVI值大。将4月16日、8月4日和12月16日NDVI值各作为一个波段，组合得到RGB影像，即图3(g)所示影像。综合图3(b)、3(d)、3 (f)分析结果，判断图3(g)中黄色部分为南荻，故提取该黄色部分，即为南荻空间分布；其中，各时相南荻与其他植被NDVI差异情况如表4所示。
表4各时相南荻与其他植被NDVI差异
三、光谱特征变量提取
由于地物存在“同物异谱、同谱异物”现象，图3(g)中黄色部分除南荻外，可能存在其他植被，还需要考虑形状特征，区分其他植被。因而本研究综合光谱和形状特征，选取NDVI、NDWI、Ratio Green、Shape Index、Length/width作为分类参数，其中，光谱特征以及形状特征的规则描述如表5所示。
表5特征变量及规则描述
四、阈值选取
由于每年各分类参数阈值由于植被生长差异，因而具有上下浮动，以2017年为例，各分类参数阈值见表6：
表6南荻提取各分类参数阈值
五、决策树构建
如图4所示，依据洞庭湖湿地主要植被南荻、杨树、苔草、水稻和油菜因生育期不同所表现在影像上的颜色、纹理及光谱等的差异，选取多个影像分割尺度，通过设置分类参数归一化植被指数(NDVI)、归一化水体指数(NDWI)、绿度(Ratio Green)、亮度值(Brightness)、形状指数(Shape Index)、长宽比(Length/width)等的阈值，结合人工目视修正，提取南荻分布信息。
六、南荻空间分布提取结果的应用分析
采用上述南荻的生长区域分割方法对2015-2019年东洞庭湖南卫星图像上的南荻进行分割，得到如图5所示的南荻空间分布图，其中，图5(a)-(e)分别为2015-2019年的南荻空间分布图，从图5可以得出：
在2015-2019年这几年，东洞庭湖南荻面积大致呈下降趋势，2015年南荻面积最大。但 2019年南荻面积增加，经分析有两点原因：
①2019年所用影像数据均是Sentinel-2，空间分辨率为10米，而其他年份数据空间分辨率为15或16米，2019年影像空间分辨率提升，混合像元减少，南荻辨识度也相应提升，并且较小面积分布的南荻也能提取出来。
②洞庭湖杨树砍伐后，湿地生态恢复，可能导致部分区域生长南荻。
③2019年南荻所用影像为6月以前影像，未利用到8月(或9月)、12月影像，2019 年提取结果中可能混有草类，导致提取面积偏大。
综上所述，本发明中的南荻生长区域确定、分割方法及系统，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；不同时期对应南荻与其他植被不同的生长阶段；获取待监测区域内各个子区域在最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在光谱特征阈值内，则判断子区域为南荻的生长区域。相比现有技术，本发明能利用南荻在不同生长阶段的光谱特征来区分待监测区域的南荻与其他植被，从而准确获取到待监测区域南荻的分布情况。
在优选方案中，本发明中通过设置分类参数归一化植被指数(NDVI)、归一化水体指数 (NDWI)、绿度(Ratio Green)、亮度值(Brightness)、形状指数(Shape Index)、长宽比(Length/width)等的阈值来区分南荻与其他植被，能进一步提高待监测区域南荻的分布情况的准确性。
在优选方案中，本技术方案在从待监测区域的卫星图像中分割南荻的生长区域时，针对不同的分割对象，采用不同的分割尺度，在保证分割的精确性的同时，能有效提高分割的速度。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种南荻生长区域确定方法，其特征在于，包括以下步骤：
根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值；所述不同时期对应南荻与其他植被不同的生长阶段；
获取待监测区域内各个子区域在所述最佳监测时相的光谱特征数据，分别将各个子区域的光谱特征数据与所述光谱特征阈值进行比较，当存在任意子区域的光谱特征数据在所述光谱特征阈值内，则判断所述子区域为南荻的生长区域。
2.根据权利要求1所述的南荻生长区域确定方法，其特征在于，根据待监测区域内南荻与其他植被在不同时期的光谱特征数据之间的差异性确定南荻在待监测区域内的最佳监测时相，包括以下步骤：
分别从历史数据中获取所述南荻以及其他植被在不同时期的原始光谱特征数据，构建南荻以及其他植被生长过程的时序光谱特征序列；
对于其他植被的每一种植被A的时序光谱特征序列，均执行以下步骤：
分别计算所述植被A的时序光谱特征序列内不同时期的光谱特征数据与其对应时期的南荻的光谱特征数据之间的差值，选取所述植被A与南荻光谱特征数据之间的差值最大的时期或所述植被A与南荻光谱特征数据之间的差值大于预设差阈值的时期作为所述植被A与所述南荻的区分时期；
基于所述南荻与其他植被中每一种植被的区分时期构建南荻在待监测区域内的最佳监测时相。
3.根据权利要求1所述的南荻生长区域确定方法，其特征在于，构建南荻以及其他植被生长过程的时序光谱特征序列，包括以下步骤：
基于南荻或其他植被在不同时期的原始光谱特征数据拟合所述南荻或其他植被的时序光谱特征曲线，判断所述时序光谱特征曲线上各个部分的光谱特征数据时序变化是否在对应的正常变化范围内，若存在任意部分的光谱特征数据的时序变化不在对应的正常变化范围内，则判断所述部分对应原始光谱特征数据有噪声干扰，修正所述部分的对应的原始光谱特征数据，使修正后的光谱特征数据拟合的时序光谱特征曲线上各个部分的光谱特征数据时序变化均在其对应的正常范围内；再使用修正后光谱特征数据构建时序光谱特征序列。
4.根据权利要求2所述的南荻生长区域确定方法，其特征在于，基于所述南荻与其他植被中每一种植被的区分时期构建南荻在待监测区域内的最佳监测时相，具体包括以下步骤：
若存在任意两种以上植被的区分时期存在重合，选取重合的区分时期作为所述两种以上植被的最优区分时期；
对于区分时期与其他植被不存在重合的植被，选取所述植被的区分时期中差值最大的时期作为最优时期；
基于其他植被中每一种植被的最优区分时期构建南荻在待监测区域内的最佳监测时相。
5.根据权利要求1所述的南荻生长区域确定方法，其特征在于，所述光谱特征数据包括：归一化植被指数NDVI、归一化水体指数NDWI、绿度Ratio Green中一种或任意几种的组合。
6.根据权利要求5所述的南荻生长区域确定方法，其特征在于，当待监测区域为大陆性亚热带季风湿润气候，且所述其他植被包括水稻、苔草、杨树、油菜时，所述光谱特征数据包括归一化植被指数NDVI、归一化水体指数NDWI、绿度Ratio Green时，所述南荻在待监测区域内的最佳监测时相及其对应的光谱特征阈值如下：
7.根据权利要求1-6任意一项所述的南荻生长区域确定方法，其特征在于，还包括以下步骤：
获取待监测区域内各个子区域的植被的形状指数、长宽比指数；
将待监测区域内各个子区域的植被的形状指数分别与预设的南荻的形状阈值进行比较，将待监测区域内各个子区域植被的长宽比指数分别与预设的南荻的长宽比阈值进行比较，当存在任意子区域的光谱特征数据在南荻的光谱特征阈值内，所述子区域植被的形状指数在南荻的形状阈值内，且所述子区域植被的长宽比指数在南荻的长宽比阈值内，则判断所述子区域为南荻的生长区域。
8.根据权利要求7所述的南荻生长区域确定方法，其特征在于，所述形状阈值为＜10.2，所述长宽比阈值为＜8.7。
9.一种南荻生长区域分割方法，其特征在于，包括以下步骤：
确定待分割区域的水体区域和非水体区域，采用第一分割尺度分割所述待监测区域的水体区域和非水体区域；
确定所述非水体区域上的道路、建筑区域和植被区域，采用第二分割尺度从所述非水体区域内分割出的道路、建筑区域，得到植被区域；
在所述植被区域内，采用权利要求1-6任意一项所述的南荻生长区域确定方法确定所述植被区域内的南荻生长区域，并采用第三分割尺度分割所述南荻的生长区域；其中，第一分割尺度大于第二分割尺度，且第二分割尺度大于第三分割尺度，所述第一分割尺度为200，所述第二分割尺度为150；所述第三分割尺度为50。
10.一种计算机系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至9任一所述算法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110772118.6A CN113642579A (zh) | 2021-07-08 | 2021-07-08 | 南荻生长区域确定、分割方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110772118.6A CN113642579A (zh) | 2021-07-08 | 2021-07-08 | 南荻生长区域确定、分割方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113642579A true CN113642579A (zh) | 2021-11-12 |
Family
ID=78416861
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110772118.6A Pending CN113642579A (zh) | 2021-07-08 | 2021-07-08 | 南荻生长区域确定、分割方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113642579A (zh) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6167167A (en) * | 1996-07-05 | 2000-12-26 | Canon Kabushiki Kaisha | Image extractions apparatus and method |
| CN105372672A (zh) * | 2015-12-07 | 2016-03-02 | 广州地理研究所 | 基于时间序列数据的南方冬种作物种植面积提取方法 |
| JP2016111939A (ja) * | 2014-12-11 | 2016-06-23 | 富士通株式会社 | 耕作放棄地判定方法、耕作放棄地判定プログラム及び耕作放棄地判定装置 |
| CN106022224A (zh) * | 2016-05-12 | 2016-10-12 | 北京师范大学 | 一种冬小麦识别方法 |
| US20170161627A1 (en) * | 2015-12-02 | 2017-06-08 | The Climate Corporation | Forecasting field level crop yield during a growing season |
| CN108364284A (zh) * | 2018-01-18 | 2018-08-03 | 中国科学院遥感与数字地球研究所 | 一种基于遥感的罂粟提取方法及系统 |
| CN109635731A (zh) * | 2018-12-12 | 2019-04-16 | 中国科学院深圳先进技术研究院 | 一种识别有效耕地的方法及装置、存储介质及处理器 |
| CN109684929A (zh) * | 2018-11-23 | 2019-04-26 | 中国电建集团成都勘测设计研究院有限公司 | 基于多源遥感数据融合的陆生植物生态环境监测方法 |
| CN109816707A (zh) * | 2018-12-25 | 2019-05-28 | 中铁第四勘察设计院集团有限公司 | 一种基于高分辨率卫星影像的露天采矿场信息提取方法 |
| CN111079846A (zh) * | 2019-12-20 | 2020-04-28 | 中国科学院遥感与数字地球研究所 | 基于时间序列高分辨率遥感数据的苹果识别方法 |
| CN111345214A (zh) * | 2020-03-13 | 2020-06-30 | 武汉禾大科技有限公司 | 一种基于卫星影像数据的新疆棉花区域识别方法及系统 |
| CN111666815A (zh) * | 2020-05-06 | 2020-09-15 | 武汉大学 | 一种基于Sentinel-2遥感影像的大蒜种植信息自动化提取方法 |
| CN112085781A (zh) * | 2020-09-08 | 2020-12-15 | 中国农业科学院农业资源与农业区划研究所 | 一种基于光谱重构技术提取冬小麦种植面积的方法 |
| CN112800973A (zh) * | 2021-01-29 | 2021-05-14 | 宁波大学 | 一种基于植被物候特征决策的互花米草提取方法 |
| CN113009485A (zh) * | 2021-03-10 | 2021-06-22 | 安徽皖南烟叶有限责任公司 | 一种基于改进植被指数的遥感烟田识别方法 |
-
2021
- 2021-07-08 CN CN202110772118.6A patent/CN113642579A/zh active Pending
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6167167A (en) * | 1996-07-05 | 2000-12-26 | Canon Kabushiki Kaisha | Image extractions apparatus and method |
| JP2016111939A (ja) * | 2014-12-11 | 2016-06-23 | 富士通株式会社 | 耕作放棄地判定方法、耕作放棄地判定プログラム及び耕作放棄地判定装置 |
| US20170161627A1 (en) * | 2015-12-02 | 2017-06-08 | The Climate Corporation | Forecasting field level crop yield during a growing season |
| CN105372672A (zh) * | 2015-12-07 | 2016-03-02 | 广州地理研究所 | 基于时间序列数据的南方冬种作物种植面积提取方法 |
| CN106022224A (zh) * | 2016-05-12 | 2016-10-12 | 北京师范大学 | 一种冬小麦识别方法 |
| CN108364284A (zh) * | 2018-01-18 | 2018-08-03 | 中国科学院遥感与数字地球研究所 | 一种基于遥感的罂粟提取方法及系统 |
| CN109684929A (zh) * | 2018-11-23 | 2019-04-26 | 中国电建集团成都勘测设计研究院有限公司 | 基于多源遥感数据融合的陆生植物生态环境监测方法 |
| CN109635731A (zh) * | 2018-12-12 | 2019-04-16 | 中国科学院深圳先进技术研究院 | 一种识别有效耕地的方法及装置、存储介质及处理器 |
| CN109816707A (zh) * | 2018-12-25 | 2019-05-28 | 中铁第四勘察设计院集团有限公司 | 一种基于高分辨率卫星影像的露天采矿场信息提取方法 |
| CN111079846A (zh) * | 2019-12-20 | 2020-04-28 | 中国科学院遥感与数字地球研究所 | 基于时间序列高分辨率遥感数据的苹果识别方法 |
| CN111345214A (zh) * | 2020-03-13 | 2020-06-30 | 武汉禾大科技有限公司 | 一种基于卫星影像数据的新疆棉花区域识别方法及系统 |
| CN111666815A (zh) * | 2020-05-06 | 2020-09-15 | 武汉大学 | 一种基于Sentinel-2遥感影像的大蒜种植信息自动化提取方法 |
| CN112085781A (zh) * | 2020-09-08 | 2020-12-15 | 中国农业科学院农业资源与农业区划研究所 | 一种基于光谱重构技术提取冬小麦种植面积的方法 |
| CN112800973A (zh) * | 2021-01-29 | 2021-05-14 | 宁波大学 | 一种基于植被物候特征决策的互花米草提取方法 |
| CN113009485A (zh) * | 2021-03-10 | 2021-06-22 | 安徽皖南烟叶有限责任公司 | 一种基于改进植被指数的遥感烟田识别方法 |
Non-Patent Citations (3)
| Title |
|---|
| GUAN XUDONG，等: "Mapping Rice Cropping Systems in Vietnam Using an NDVI-Based Time-Series Similarity Measurement Based on DTW Distance", 《REMOTE SENSING》, vol. 8, no. 1, pages 19 * |
| 李中元，等: "利用物候差异与面向对象决策树提取油菜种植面积", 《地球信息科学学报》, vol. 21, no. 5, pages 720 - 730 * |
| 杜保佳，等: "应用Sentinel-2ANDVI时间序列和面向对象决策树 方法的农作物分类", 《地球信息科学学报》, vol. 21, no. 5, pages 740 - 751 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108363949B (zh) | 一种基于物候分析的棉花遥感监测方法 | |
| CN107609526A (zh) | 基于规则的精细尺度城市不透水面快速提取方法 | |
| CN111626269B (zh) | 一种实用的大空间范围滑坡提取方法 | |
| CN112800973B (zh) | 一种基于植被物候特征决策的互花米草提取方法 | |
| CN108830844B (zh) | 一种基于多时相高分辨率遥感影像的设施蔬菜提取方法 | |
| CN108710864B (zh) | 基于多维度识别及图像降噪处理的冬小麦遥感提取方法 | |
| CN108986116B (zh) | 一种基于遥感影像的红树林提取方法及系统 | |
| CN113033670A (zh) | 一种基于Sentinel-2A/B数据的水稻种植面积提取方法 | |
| CN105447274B (zh) | 一种利用面向对象分类技术对中等分辨率遥感图像进行滨海湿地制图的方法 | |
| CN114387528A (zh) | 松材线虫病监测天空地一体化监测方法 | |
| CN111563430A (zh) | 一种基于卷积神经的公益林地类图像智能识别方法及系统 | |
| CN110334623B (zh) | 一种基于Sentinel-2A卫星遥感影像提取崩岗信息的方法 | |
| CN113971769A (zh) | 基于多源大数据的海岸带地域功能长时序识别方法 | |
| CN116912578A (zh) | 一种农作物分类方法、系统及电子设备 | |
| CN114049564A (zh) | 基于高光谱遥感影像的松材线虫病等级预测模型构建方法 | |
| CN113642579A (zh) | 南荻生长区域确定、分割方法及系统 | |
| Roy et al. | Comparative analysis of object based and pixel based classification for mapping of mango orchards in Sitapur district of Uttar Pradesh | |
| CN115049920B (zh) | 基于gee结合生育期属性和面向对象的大区域作物提取方法 | |
| CN116129284A (zh) | 一种基于时间序列变化特征的撂荒地遥感提取方法 | |
| CN113887940B (zh) | 一种基于卷积神经网络的牧草优势度识别方法 | |
| Novković et al. | Towards UAV Assisted Monitoring Of An Aquatic Vegetation Within The Large Rivers–The Middle Danube | |
| CN113177441A (zh) | 面向对象与物候知识相融合的遥感互花米草制图方法 | |
| Shao et al. | Large-scale mapping of new mixed rice cropping patterns in southern China with phenology-based algorithm and MODIS dataset | |
| Kuria et al. | Multi-temporal land cover mapping of the Kakamega Forest utilising LANDSAT imagery and GIS | |
| CN117409313B (zh) | 一种基于Sentinel-2光学影像的互花米草物候衰退期指数构建方法 |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Li Inventor after: Liu Zhijie Inventor after: Hu Wen Inventor before: Chen Li Inventor before: Liu Zhijie |
|
| CB03 | Change of inventor or designer information |