CN114612794A - 一种细碎化农业区土地覆被与种植结构的遥感识别方法 - Google Patents
一种细碎化农业区土地覆被与种植结构的遥感识别方法 Download PDFInfo
- Publication number
- CN114612794A CN114612794A CN202210193796.1A CN202210193796A CN114612794A CN 114612794 A CN114612794 A CN 114612794A CN 202210193796 A CN202210193796 A CN 202210193796A CN 114612794 A CN114612794 A CN 114612794A
- Authority
- CN
- China
- Prior art keywords
- agricultural
- classification
- land
- agricultural area
- growth period
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 238000012706 support-vector machine Methods 0.000 claims abstract description 20
- 238000003066 decision tree Methods 0.000 claims abstract description 13
- 238000012549 training Methods 0.000 claims abstract description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 31
- 238000012952 Resampling Methods 0.000 claims description 12
- 238000002310 reflectometry Methods 0.000 claims description 12
- 230000001174 ascending effect Effects 0.000 claims description 8
- 239000003086 colorant Substances 0.000 claims description 7
- 238000007635 classification algorithm Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 5
- 241000196324 Embryophyta Species 0.000 claims description 4
- 235000003222 Helianthus annuus Nutrition 0.000 claims description 4
- 241000209140 Triticum Species 0.000 claims description 4
- 235000021307 Triticum Nutrition 0.000 claims description 4
- 240000008042 Zea mays Species 0.000 claims description 4
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 claims description 4
- 235000002017 Zea mays subsp mays Nutrition 0.000 claims description 4
- 235000005822 corn Nutrition 0.000 claims description 4
- 235000012055 fruits and vegetables Nutrition 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000013139 quantization Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000011160 research Methods 0.000 claims description 2
- 244000020551 Helianthus annuus Species 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000010287 polarization Effects 0.000 description 4
- 241000208818 Helianthus Species 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 230000003595 spectral 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
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Image Processing (AREA)
Abstract
本发明属于农业遥感领域,涉及一种细碎化农业区土地覆被与种植结构的遥感识别方法,采用机器学习算法提取土地利用,不依赖实测样本点,自动训练后得到高精度农业区土地利用栅格化数据。本发明通过监督算法支持向量机SVM先对农业区及非农业区进行分类,而后基于已得到的农业区分类栅格数据,采用决策树算法进一步对作物类型进行分类判别。该方法不依赖于实测数据,大大降低了分类时间及经济成本,且有效避免了传统机器算法样本训练集质量对分类结果造成的负面影响。同时该方法在保证分类高效性的前提下,满足农业遥感应用所需的土地利用制图精度。
Description
技术领域
本发明属于农业遥感领域,涉及一种细碎化农业区土地覆被与种植结构的遥感识别方法,采用机器学习算法提取土地利用,不依赖实测样本点,自动训练后得到高精度农业区土地利用栅格化数据。
背景技术
准确及时地获取农业区土地覆被信息对监测农情,调整农业结构,估产和制定粮食政策等均具有重要意义。传统的田间调查方法费时费力,且难以获取大尺度区域土地利用信息。遥感技术能克服单点观测在复杂地表条件下空间代表性差的缺点,为区域土地利用识别提供多时相、多光谱及多角度的地表信息。
目前土地利用遥感反演主要分为基于多源数据物候和极化特征及基于实测样本点分类两种方式。基于样本点的识别算法只需输入样本训练集,利用卫星遥感图像,即可快速实现大尺度区域土地利用反演。其中支持向量机(SVM)和随机森林等算法因其训练效率高,精度满足应用要求,被广泛应用于土地利用遥感识别,但该方法对样本精度要求高,且应用时需要大量实测样本,时间及金钱成本较高。同时对于光谱及纹理特征较为类似的作物来说,该类算法在识别时会出现作物类型混分的情况。基于多源物候及极化特征算法只需建立特定的分类规则,不需实测样本训练集,即可利用决策树算法实现土地利用反演。这种方法只针对于物候/极化特征差异十分明显的地物识别,且需对原始卫星图像光谱数据进行处理分析,计算得到物候/极化特征用于分类。总而言之,现有的分类算法在应用于卫星遥感图像时,单一分类算法运算时间长,对计算机配置要求高,需要大量实测样本,费时费力,因此难以实现大尺度土地利用识别。
在实际应用中,SVM、随机森林、决策树等各种分类器都具有其优缺点,因此本发明结合各种算法的优势,针对细碎化农业区不同地物的特点应用不同算法,取长补短,同时提高分类的效率以及精度。
发明内容
针对上述技术问题,本发明的目的是提供一种细碎化农业区土地覆被与种植结构的遥感识别方法,通过卫星遥感图像光谱及纹理特征,以及识别农业区和非农业区效果较好的LSWI指数,利用机器学习算法实现农业区和非农业区的识别。在此基础上,针对农业区不同作物所具有的物候特征建立分类规则,利用NDVI时序数据计算得到的物候指标作为输入数据,基于决策树算法实现农业区不同作物的进一步划分,最终实现细碎化农业区高精度土地利用遥感识别。
为了实现上述目的,本发明提供了如下技术方案:
一种细碎化农业区土地覆被与种植结构的遥感识别方法,包括如下步骤:
S1、利用支持向量机SVM监督学习算法对原始卫星图像进行农业区及非农业区识别:
S1.1、采用双线性插值法将多张不同时期不同气象条件的原始卫星图像的每个波段重采样为统一的空间分辨率;选择一张八月初作物生长旺盛的原始卫星图像,并选取该图像中能较好分辨地物的可见光波段以及能判断植物生长状况的红边波段、近红外波段或短波红外波段作为分类识别波段;利用二阶概率统计方法用一个灰色调空间相关性矩阵,基于3×3的窗口和64灰度量化级别计算分类识别波段的纹理值;
S1.2、利用生育期内多期晴空无云条件下的原始卫星图像,通过公式一得到陆地水体指数LSWI时序数据集,采用线性插值的方法将陆地水体指数LSWI时序数据集重采样为统一的时间分辨率,并以一预设阈值对数据集中的陆地水体指数LSWI时序数据进行重新赋值,大于预设阈值的像元赋值为0,小于预设阈值的像元赋值为1,将重新赋值的陆地水体指数LSWI时序数据依次相加求和,获得一个陆地水体指数LSWI重分类数据;
式中:LSWI为陆地水体指数;ρNIR表示近红外波段的反射率;ρSWIR表示短波红外波段的反射率;
S1.3、从Google Earth中选取土地利用样本集作为训练数据,利用支持向量机SVM监督分类算法对步骤S1.1得到的分类识别波段及纹理值和步骤S1.2获得的陆地水体指数LSWI重分类数据进行训练学习,计算得到农业区及非农业区的土地利用分类结果;
S2、应用决策树算法对农业区进行进一步种植结构划分,得到最终农业区高精度土地利用数据:
S2.1、将步骤S1得到的土地利用分类结果转变为栅格数据并对其进行重新赋值,将农业区赋值为1,非农业区赋值为NoData;将重新赋值的栅格数据转变为矢量数据,获得只包含农业区的shp矢量数据;
S2.2、基于步骤S2.1得到的只包含农业区的shp矢量数据对生育期内多期晴空无云条件下的原始卫星图像进行裁剪,得到只包含农业区的原始卫星图像,并通过公式二计算得到归一化植被指数NDVI时序数据集,采用线性插值的方法将归一化植被指数NDVI时序数据集重采样为统一的时间分辨率;
式中:NDVI为归一化植被指数;ρNIR表示近红外波段的反射率;ρRed表示红光波段的反射率;
S2.3、利用Savitzky-Golay滤波器对步骤S2.2重采样后的归一化植被指数NDVI时序数据集进行平滑处理;基于平滑后的NDVI生育期过程曲线建立物候指标判断标准,对农业区每一个像元进行判断,最终得到物候指标栅格数据;所述物候指标栅格数据中包含生育期起始时间、生育期终止时间和生育期长度三个物候指标;
所述物候指标判断标准的建立过程如下:
确定在平滑后的NDVI生育期过程曲线的上升阶段、下降阶段以及NDVI设定值,当NDVI=设定值时,上升阶段所对应的时刻为生育期起始时间,下降阶段所对应的时刻为生育期终止时间,生育期起始时间与生育期终止时间内的时间长度为生育期长度;
S2.4、根据不同作物的三个物候指标的值域范围,建立作物分类规则;基于作物分类规则,利用决策树算法对步骤S2.3获得的物候指标栅格数据中的每一个栅格进行判断,获得农业种植结构分类结果,最终使得每个栅格都被划分为某一类具体的作物类型,从而实现对农业区的作物识别;
S3、将步骤S1得到的土地利用分类结果中非农业区与步骤S2获得的农业种植结构分类结果镶嵌到一起,得到一张带有分类标签色彩的农业种植区土地利用及种植结构分类图。
所述步骤S1.2中,预设阈值为0.2。
所述步骤S1.3中,土地利用分类结果中,农业区为农田,非农业区包括水体、沙丘、居民区和天然地。
所述步骤S2.3中,设定值为0.3。
所述步骤S2.4得到的农业种植结构分类结果包括小麦、玉米、葵花和果蔬。
所述步骤S2.4中,不同作物的三个物候指标的值域范围基于当地实际经验和实地调研分别确定。
所述步骤S3中,农业种植区土地利用及种植结构分类图中的每一个像元都有一个分类,不同的值和颜色代表着不同的分类。
与现有技术相比,本发明的有益效果在于:
1、本发明解决了实地样本获取数量有限,且实测样本质量要求较高,获取困难的问题,在无需实测样本的情况下,保证分类精度及稳定性。
2、本发明与传统算法相比,降低了时间成本及金钱成本,且可操作性较强,在应用于大区域大面积高精度遥感图像时,也可用普通硬件设备实现高精度土地利用识别,使得大区域的农业区土地利用自动获取成为可能。
附图说明
图1为本发明的细碎化农业区土地覆被与种植结构的遥感识别方法的流程图;
图2为本发明实施例的原始卫星遥感图像(在本描述中举例使用Sentinel-2历史图像,其他卫星原始图像均可);
图3为本发明实施例的原始卫星遥感图像纹理值计算过程示意图;
图4为本发明实施例监督算法支持向量机SVM分类的计算过程及经过分类后得到的土地利用图像;
图5为本发明选用的物候指标示意图;
图6为本发明实施例用于作物分类的决策树分类规则;
图7为本发明实施例对农业区进行具体种植结构分类后的图像;
图8为本发明实施例最终得到的土地利用及种植结构分类图;
图9为分类结果中各个类别的图例。
具体实施方式
下面结合附图和实施例对本发明进行进一步说明。
如图1所示,一种细碎化农业区土地覆被与种植结构的遥感识别方法,包括如下步骤:
S1、利用支持向量机SVM监督学习算法对原始卫星图像进行农业区及非农业区识别:
S1.1、采用双线性插值法将多张不同时期不同气象条件的原始卫星图像的每个波段重采样为统一的空间分辨率;选择一张八月初作物生长旺盛的原始卫星图像,并选取该图像中能较好分辨地物的可见光波段以及能判断植物生长状况的红边波段、近红外波段或短波红外波段作为分类识别波段;利用二阶概率统计方法用一个灰色调空间相关性矩阵,基于3×3的窗口和64灰度量化级别计算分类识别波段的纹理值;
S1.2、利用生育期内多期晴空无云条件下的原始卫星图像,通过公式一得到陆地水体指数LSWI时序数据集,采用线性插值的方法将陆地水体指数LSWI时序数据集重采样为统一的时间分辨率,并以一预设阈值对数据集中的陆地水体指数LSWI时序数据进行重新赋值,大于预设阈值的像元赋值为0,小于预设阈值的像元赋值为1,将重新赋值的陆地水体指数LSWI时序数据依次相加求和,获得一个陆地水体指数LSWI重分类数据;
式中:LSWI为陆地水体指数;ρNIR表示近红外波段的反射率;ρSWIR表示短波红外波段的反射率。
所述预设阈值为0.2。
S1.3、从Google Earth中选取土地利用样本集作为训练数据,利用支持向量机SVM监督分类算法对步骤S1.1得到的分类识别波段及纹理值和步骤S1.2获得的陆地水体指数LSWI重分类数据进行训练学习,计算得到农业区及非农业区的土地利用分类结果;
所述土地利用分类结果中,农业区为农田,非农业区包括水体、沙丘、居民区和天然地。
S2、应用决策树算法对农业区进行进一步种植结构划分,得到最终农业区高精度土地利用数据:
S2.1、将步骤S1得到的土地利用分类结果转变为栅格数据并对其进行重新赋值,将农业区赋值为1,非农业区赋值为NoData;将重新赋值的栅格数据转变为矢量数据,获得只包含农业区的shp矢量数据;
S2.2、基于步骤S2.1得到的只包含农业区的shp矢量数据对生育期内多期晴空无云条件下的原始卫星图像进行裁剪,得到只包含农业区的原始卫星图像,并通过公式二计算得到归一化植被指数NDVI时序数据集,采用线性插值的方法将归一化植被指数NDVI时序数据集重采样为统一的时间分辨率;
式中:NDVI为归一化植被指数;ρNIR表示近红外波段的反射率;ρRed表示红光波段的反射率。
S2.3、利用Savitzky-Golay(S-G)滤波器对步骤S2.2重采样后的归一化植被指数NDVI时序数据集进行平滑处理;基于平滑后的NDVI生育期过程曲线建立物候指标判断标准,对农业区每一个像元进行判断,最终得到物候指标栅格数据;所述物候指标栅格数据中包含生育期起始时间、生育期终止时间和生育期长度三个物候指标;
所述物候指标判断标准的建立过程如下:
确定在平滑后的NDVI生育期过程曲线的上升阶段、下降阶段以及NDVI设定值,当NDVI=设定值时,上升阶段所对应的时刻为生育期起始时间,下降阶段所对应的时刻为生育期终止时间,生育期起始时间与生育期终止时间内的时间长度为生育期长度。
所述设定值为0.3。
S2.4、根据不同作物的三个物候指标的值域范围,建立作物分类规则;基于作物分类规则,利用决策树算法对步骤S2.3获得的物候指标栅格数据中的每一个栅格进行判断,获得农业种植结构分类结果,最终使得每个栅格都被划分为某一类具体的作物类型,从而实现对农业区的作物识别。
所述步骤S2得到的农业种植结构分类结果包括小麦、玉米、葵花和果蔬。
不同作物的三个物候指标的值域范围基于当地实际经验和实地调研分别确定。
S3、将步骤S1得到的土地利用分类结果中非农业区与步骤S2获得的农业种植结构分类结果镶嵌到一起,得到一张带有分类标签色彩的农业种植区土地利用及种植结构分类图。
所述农业种植区土地利用及种植结构分类图中的每一个像元都有一个分类,不同的值和颜色代表着不同的分类。
实施例
S1、利用支持向量机(SVM)监督学习算法对原始卫星图像进行农业区及非农业区识别;
将Sentinel-2原始卫星图像每个波段采用双线性插值法重采样为统一的10m空间分辨率,如图2所示;选择一张八月初作物生长旺盛的原始卫星图像,选取该图像中能较好分辨地物的可见光(红绿蓝)波段,以及能判断植物生长状况的红边波段、近红外波段以及短波红外波段作为分类识别波段。利用二阶概率统计方法用一个灰色调空间相关性矩阵,基于3×3的窗口和64灰度量化级别对最优波段的纹理值进行计算,最终得到基于八个纹理滤波(均值、方差、协同性、对比度、相异性、信息熵、二阶矩和相关性)依次计算得到的纹理值,如图3所示;
利用生育期内多期晴空无云条件下的Sentinel-2原始卫星图像,经过公式一计算得到陆地水体指数LSWI时序数据集,采用线性插值的方法将陆地水体指数LSWI数据集重采样为统一的5天时间分辨率,并以0.2为阈值对其进行重新赋值,大于0.2的像元赋值为0,小于0.2的像元赋值为1,将重新赋值的陆地水体指数LSWI时序数据依次相加求和,获得一个陆地水体指数LSWI重分类数据;
从Google Earth中人工选取土地利用样本集(居民地、沙丘、荒地、水体和农田)作为训练数据,利用支持向量机SVM监督分类算法对最优波段及纹理值以及LSWI重分类数据进行训练学习,计算得到农业区及非农业区土地利用分类结果,如图4所示。
S2、应用决策树算法对农业区进行进一步种植结构划分,得到最终农业区高精度土地利用数据;
将土地利用分类结果转变为栅格数据并对其进行重新赋值,将农业区赋值为1,非农业区统一赋值为NoData,将重新赋值的栅格数据转变为矢量数据,得到农业区的shp矢量数据;
基于农业区矢量数据对生育期内多期晴空无云条件下的原始卫星图像进行裁剪,得到只包含农业区的原始卫星图像,通过公式二计算得到归一化植被指数NDVI时序数据集,将时序数据集采用线性插值的方法重采样为统一的5天时间分辨率;
利用Savitzky-Golay(S-G)滤波器对生育期NDVI时序数据集进行平滑处理。如图5所示,基于平滑后的NDVI生育期过程曲线,确定在曲线的上升阶段以及下降阶段,当NDVI=0.3时,上升阶段以及下降阶段所对应的时刻分别为生育期起始时间及生育期终止时间,生育期起始时间与生育期终止时间内的时间长度为生育期长度,基于上述物候指标判断标准对农业区每一个像元进行判断,最终得到包含生育期起始时间、生育期终止时间以及生育期长度三个物候指标的物候指标栅格数据,用于进一步的作物分类。
基于当地实际经验和实地调研,分别确定不同作物的三个物候指标的值域范围,建立作物分类规则,如图6所示。基于分类规则,利用决策树算法对物候指标栅格数据中的每一个栅格进行判断,最终使得每个栅格都被划分为某一类具体的作物类型,从而实现对农业区的作物识别,得到如图7所示的农业区种植结构分类结果。
将土地利用结果中的非农业区与种植结构识别结果镶嵌到一起,如图8所示,进而得到农业区高精度土地利用分类结果,分类的结果是得到带有分类标签色彩的一张图像,包括水体、沙丘、居民区、天然地非农业区土地利用类型以及小麦、玉米、葵花和果蔬种植结构分类结果,各类地物及作物的图例如图9所示。
Claims (7)
1.一种细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,包括如下步骤:
S1、利用支持向量机SVM监督学习算法对原始卫星图像进行农业区及非农业区识别:
S1.1、采用双线性插值法将多张不同时期不同气象条件的原始卫星图像的每个波段重采样为统一的空间分辨率;选择一张八月初作物生长旺盛的原始卫星图像,并选取该图像中能较好分辨地物的可见光波段以及能判断植物生长状况的红边波段、近红外波段或短波红外波段作为分类识别波段;利用二阶概率统计方法用一个灰色调空间相关性矩阵,基于3×3的窗口和64灰度量化级别计算分类识别波段的纹理值;
S1.2、利用生育期内多期晴空无云条件下的原始卫星图像,通过公式一得到陆地水体指数LSWI时序数据集,采用线性插值的方法将陆地水体指数LSWI时序数据集重采样为统一的时间分辨率,并以一预设阈值对数据集中的陆地水体指数LSWI时序数据进行重新赋值,大于预设阈值的像元赋值为0,小于预设阈值的像元赋值为1,将重新赋值的陆地水体指数LSWI时序数据依次相加求和,获得一个陆地水体指数LSWI重分类数据;
式中:LSWI为陆地水体指数;ρNIR表示近红外波段的反射率;ρSWIR表示短波红外波段的反射率;
S1.3、从Google Earth中选取土地利用样本集作为训练数据,利用支持向量机SVM监督分类算法对步骤S1.1得到的分类识别波段及纹理值和步骤S1.2获得的陆地水体指数LSWI重分类数据进行训练学习,计算得到农业区及非农业区的土地利用分类结果;
S2、应用决策树算法对农业区进行进一步种植结构划分,得到最终农业区高精度土地利用数据:
S2.1、将步骤S1得到的土地利用分类结果转变为栅格数据并对其进行重新赋值,将农业区赋值为1,非农业区赋值为NoData;将重新赋值的栅格数据转变为矢量数据,获得只包含农业区的shp矢量数据;
S2.2、基于步骤S2.1得到的只包含农业区的shp矢量数据对生育期内多期晴空无云条件下的原始卫星图像进行裁剪,得到只包含农业区的原始卫星图像,并通过公式二计算得到归一化植被指数NDVI时序数据集,采用线性插值的方法将归一化植被指数NDVI时序数据集重采样为统一的时间分辨率;
式中:NDVI为归一化植被指数;ρNIR表示近红外波段的反射率;ρRed表示红光波段的反射率;
S2.3、利用Savitzky-Golay滤波器对步骤S2.2重采样后的归一化植被指数NDVI时序数据集进行平滑处理;基于平滑后的NDVI生育期过程曲线建立物候指标判断标准,对农业区每一个像元进行判断,最终得到物候指标栅格数据;所述物候指标栅格数据中包含生育期起始时间、生育期终止时间和生育期长度三个物候指标;
所述物候指标判断标准的建立过程如下:
确定在平滑后的NDVI生育期过程曲线的上升阶段、下降阶段以及NDVI设定值,当NDVI=设定值时,上升阶段所对应的时刻为生育期起始时间,下降阶段所对应的时刻为生育期终止时间,生育期起始时间与生育期终止时间内的时间长度为生育期长度;
S2.4、根据不同作物的三个物候指标的值域范围,建立作物分类规则;基于作物分类规则,利用决策树算法对步骤S2.3获得的物候指标栅格数据中的每一个栅格进行判断,获得农业种植结构分类结果,最终使得每个栅格都被划分为某一类具体的作物类型,从而实现对农业区的作物识别;
S3、将步骤S1得到的土地利用分类结果中非农业区与步骤S2获得的农业种植结构分类结果镶嵌到一起,得到一张带有分类标签色彩的农业种植区土地利用及种植结构分类图。
2.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S1.2中,预设阈值为0.2。
3.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S1.3中,土地利用分类结果中,农业区为农田,非农业区包括水体、沙丘、居民区和天然地。
4.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S2.3中,设定值为0.3。
5.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S2.4得到的农业种植结构分类结果包括小麦、玉米、葵花和果蔬。
6.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S2.4中,不同作物的三个物候指标的值域范围基于当地实际经验和实地调研分别确定。
7.根据权利要求1所述的细碎化农业区土地覆被与种植结构的遥感识别方法,其特征在于,所述步骤S3中,农业种植区土地利用及种植结构分类图中的每一个像元都有一个分类,不同的值和颜色代表着不同的分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210193796.1A CN114612794B (zh) | 2022-03-01 | 2022-03-01 | 一种细碎化农业区土地覆被与种植结构的遥感识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210193796.1A CN114612794B (zh) | 2022-03-01 | 2022-03-01 | 一种细碎化农业区土地覆被与种植结构的遥感识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114612794A true CN114612794A (zh) | 2022-06-10 |
CN114612794B CN114612794B (zh) | 2024-06-07 |
Family
ID=81859076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210193796.1A Active CN114612794B (zh) | 2022-03-01 | 2022-03-01 | 一种细碎化农业区土地覆被与种植结构的遥感识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114612794B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114782838A (zh) * | 2022-06-17 | 2022-07-22 | 中化现代农业有限公司 | 水稻识别方法、装置、电子设备和存储介质 |
CN117437475A (zh) * | 2023-11-02 | 2024-01-23 | 清华大学 | 种植结构分类方法、装置、计算机设备和存储介质 |
CN118279431A (zh) * | 2024-06-04 | 2024-07-02 | 中国农业科学院农业资源与农业区划研究所 | 大区域低样本依赖度的作物制图方法与*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101728137B1 (ko) * | 2016-02-04 | 2017-04-19 | (주)한라지리정보 | 위성 영상과 gis를 사용한 토지피복 항목별 영상분류 방법 |
CN110796001A (zh) * | 2019-09-23 | 2020-02-14 | 武汉珈和科技有限公司 | 一种卫星影像覆膜农田识别与提取方法及*** |
CN110852262A (zh) * | 2019-11-11 | 2020-02-28 | 南京大学 | 基于时间序列高分一号遥感影像的农业用地提取方法 |
WO2021007665A1 (en) * | 2019-07-17 | 2021-01-21 | Farmers Edge Inc. | Automatic crop classification system and method |
CN112395914A (zh) * | 2019-08-15 | 2021-02-23 | 中国科学院遥感与数字地球研究所 | 一种融合遥感影像时序与纹理特征的地块作物识别方法 |
CN112906666A (zh) * | 2021-04-07 | 2021-06-04 | 中国农业大学 | 一种农业种植结构的遥感识别方法 |
CN113657158A (zh) * | 2021-07-13 | 2021-11-16 | 西安电子科技大学 | 基于Google Earth Engine的大尺度大豆种植区域提取算法 |
-
2022
- 2022-03-01 CN CN202210193796.1A patent/CN114612794B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101728137B1 (ko) * | 2016-02-04 | 2017-04-19 | (주)한라지리정보 | 위성 영상과 gis를 사용한 토지피복 항목별 영상분류 방법 |
WO2021007665A1 (en) * | 2019-07-17 | 2021-01-21 | Farmers Edge Inc. | Automatic crop classification system and method |
CN112395914A (zh) * | 2019-08-15 | 2021-02-23 | 中国科学院遥感与数字地球研究所 | 一种融合遥感影像时序与纹理特征的地块作物识别方法 |
CN110796001A (zh) * | 2019-09-23 | 2020-02-14 | 武汉珈和科技有限公司 | 一种卫星影像覆膜农田识别与提取方法及*** |
CN110852262A (zh) * | 2019-11-11 | 2020-02-28 | 南京大学 | 基于时间序列高分一号遥感影像的农业用地提取方法 |
CN112906666A (zh) * | 2021-04-07 | 2021-06-04 | 中国农业大学 | 一种农业种植结构的遥感识别方法 |
CN113657158A (zh) * | 2021-07-13 | 2021-11-16 | 西安电子科技大学 | 基于Google Earth Engine的大尺度大豆种植区域提取算法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114782838A (zh) * | 2022-06-17 | 2022-07-22 | 中化现代农业有限公司 | 水稻识别方法、装置、电子设备和存储介质 |
CN117437475A (zh) * | 2023-11-02 | 2024-01-23 | 清华大学 | 种植结构分类方法、装置、计算机设备和存储介质 |
CN118279431A (zh) * | 2024-06-04 | 2024-07-02 | 中国农业科学院农业资源与农业区划研究所 | 大区域低样本依赖度的作物制图方法与*** |
Also Published As
Publication number | Publication date |
---|---|
CN114612794B (zh) | 2024-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114612794B (zh) | 一种细碎化农业区土地覆被与种植结构的遥感识别方法 | |
CN110751019B (zh) | 基于深度学习的高分辨率影像农作物自动提取方法及装置 | |
US7058197B1 (en) | Multi-variable model for identifying crop response zones in a field | |
CN102829739B (zh) | 一种面向对象的作物叶面积指数遥感反演方法 | |
CN111104858A (zh) | 一种基于形态模型法的大范围作物物候提取方法 | |
CN113033279A (zh) | 一种基于多源遥感影像的作物精细分类方法及其*** | |
CN112766155A (zh) | 一种基于深度学习的海水养殖区提取方法 | |
CN112434569B (zh) | 一种无人机热成像*** | |
CN111199192A (zh) | 一种采用平行线取样的田间红提葡萄整体成熟度检测方法 | |
CN111832480B (zh) | 一种基于光谱特征的油菜种植区遥感识别方法 | |
CN112949607A (zh) | 一种基于JM Relief F湿地植被特征优选与融合方法 | |
CN116721385A (zh) | 一种基于机器学习的rgb摄像头数据蓝藻水华监测方法 | |
CN115063690B (zh) | 一种基于ndvi时序特征的植被分类方法 | |
CN115424006A (zh) | 应用于作物表型参数反演的多源多层次数据融合方法 | |
CN115471753A (zh) | 一种结合专家知识和机器学习的树种精细分层分类方法 | |
CN111860038B (zh) | 一种农作物前端识别装置及方法 | |
CN112179853A (zh) | 基于影像阴影去除的果树冠层氮素含量遥感反演方法及*** | |
Hall et al. | A method for extracting detailed information from high resolution multispectral images of vineyards | |
CN116052141B (zh) | 作物生育期的识别方法、装置、设备及介质 | |
CN114119536B (zh) | 基于云平台技术的作物氮含量快速监测与诊断方法 | |
Wang et al. | Scaling effects of chlorophyll content in walnut leaves estimations with coupling Boruta algorithm and machine learning model | |
CN115266649A (zh) | 一种结合无人机影像纹理特征的果树叶片氮含量估算方法 | |
CN117371603A (zh) | 基于多模态数据和集成学习的小麦倒伏参数预测方法 | |
CN113723195A (zh) | 基于图像识别技术的农作物旱情识别***及管理平台 | |
CN117011702A (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 |