CN113313671A - 一种结合卫星影像和modis数据的水稻估产方法 - Google Patents
一种结合卫星影像和modis数据的水稻估产方法 Download PDFInfo
- Publication number
- CN113313671A CN113313671A CN202110488775.8A CN202110488775A CN113313671A CN 113313671 A CN113313671 A CN 113313671A CN 202110488775 A CN202110488775 A CN 202110488775A CN 113313671 A CN113313671 A CN 113313671A
- Authority
- CN
- China
- Prior art keywords
- rice
- evi
- extracting
- yield
- image
- 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
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 167
- 235000009566 rice Nutrition 0.000 title claims abstract description 167
- 238000000034 method Methods 0.000 title claims abstract description 63
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 166
- 238000000605 extraction Methods 0.000 claims abstract description 39
- 238000003066 decision tree Methods 0.000 claims abstract description 18
- 230000000007 visual effect Effects 0.000 claims abstract description 12
- 230000001419 dependent effect Effects 0.000 claims abstract description 9
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 238000010276 construction Methods 0.000 claims abstract description 7
- 238000000611 regression analysis Methods 0.000 claims abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 5
- 230000003595 spectral effect Effects 0.000 claims description 8
- 230000004927 fusion Effects 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 7
- 238000007476 Maximum Likelihood Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 6
- 238000011835 investigation Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 3
- 235000013339 cereals Nutrition 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 239000002689 soil Substances 0.000 description 3
- 240000008042 Zea mays Species 0.000 description 2
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 2
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 235000005822 corn Nutrition 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000011362 coarse particle Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000035558 fertility Effects 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/24323—Tree-organised classifiers
-
- 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
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- 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/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种结合卫星影像和MODIS数据的水稻估产方法,对原始的卫星影像进行预处理操作;在预处理后的卫星影像上选取感兴趣区,分别采用基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种方法对感兴趣区进行分类和水稻提取,分类后的类型包括:林地、水稻、水体和建设用地;取三种分类方法得到的水稻提取结果的相交区域作为水稻区域,提取水稻区域的不同时期的MODIS影像中的EVI指数均值图像;将水稻EVI均值作为自变量,水稻产量统计数据作为因变量,进行回归分析,构建统计回归估算模型,既解决了特定时期内高分辨率影像用于多时序指数提取数据量不足的问题,又解决了低分辨率提取精度低的问题。
Description
技术领域
本发明涉及水稻估产领域,尤其涉及一种结合卫星影像和MODIS数据的水稻估产方法。
背景技术
水稻作为世界上最重要的粮食作物之一,水稻产量数据有着重要的意义。传统的水稻估产速度慢、工作量大、成本高。目前获取水稻产量数据的渠道大多来源于各级统计部门抽样调查统计的数据、农业部门调查汇总的数据、气象部门根据气候条件推算的数据,传统的水稻估产则基于以上数据,根据水稻的物理构建,通过作物面积抽样调查水稻穗禾的颗粒数,利用统计数学估算小范围的水稻产量,难以估测大面积的水稻产量。
与常规的统计方法相比,遥感技术因其具有的多光谱、多时相及高分辨率,全天候和宏观特性的特点,使其在水稻面积提取和估产方面有着独特的优势。通过遥感技术与地理信息系统和全球定位系统结合,实现更为快速、准确、动态的信息采集和分析,通过卫星传感器,获取更为丰富的地表,获取更为准确的水稻种植区域,提取不同时期水稻的相关植被指数信息。通常是使用基于时间序列NDVI(Normalized Difference Vegetation Index,归一化植被指数)与水稻产量拟合一种函数关系建立估产模型,能够在一定程度上满足估产的需求,但会受到数据源的限制,存在数据泛化的风险。另外还有一种基于无人机遥感归一化光谱植被指数来分析水稻不同时期的状态,建立一种时序模型来估测产量,但会受到地形的限制,增加了估产的难度。
发明内容
本发明针对现有技术中存在的技术问题,提供一种结合卫星影像和MODIS数据的水稻估产方法,解决现有技术中水稻估产数据低分辨率提取精度低的问题。
本发明解决上述技术问题的技术方案如下:一种结合卫星影像和MODIS数据的水稻估产方法,包括:
步骤1,对原始的卫星影像进行预处理操作;
步骤2,在预处理后的所述卫星影像上选取感兴趣区,分别采用基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种方法对所述感兴趣区进行分类和水稻提取,分类后的类型包括:林地、水稻、水体和建设用地;
步骤3,取三种分类方法得到的水稻提取结果的相交区域作为水稻区域,提取所述水稻区域的不同时期的MODIS影像中的EVI指数均值图像;
步骤4,将水稻EVI均值作为自变量,水稻产量统计数据作为因变量,进行回归分析,构建所述水稻EVI均值与所述水稻产量统计数据之间的统计回归估算模型。
本发明的有益效果是:本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法,使用高分辨率影像,通过不同提取特征、不同分类方法综合提取水稻区域,在此基础上,使用多时序影像来提取植被指数EVI,以描述植被在不同季节的差异。通过探究水稻植被指数与水稻产量之间的数学关系,建立估产的数学模型,能够更加快速、准确、有效地估测水稻的产量。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,所述步骤1中对所述卫星影像进行预处理操作包括:大气校正、影像融合、影像剪裁和影像拼接。
进一步,所述步骤2中选取所述感兴趣区包括:结合高清地图影像、实地考察水稻样本点和水稻种植标准农田数据,在所述卫星影像上使用人工标注的方式选取所述感兴趣区域。
进一步,所述步骤2中采用基于人工目视解译的监督分类进行分类和水稻提取的过程包括:采用最大似然法对所述感兴趣区进行分类和水稻提取。
进一步,所述步骤2中基于规则的面向对象分类进行分类和水稻提取的过程包括:
根据光谱范围的差异建立水稻的提取规则,使用基于所述提取规则的面向对象分类方法对所述感兴趣区进行分类;
所述提取规则包括:设置在可见光蓝色波段光谱的植被的提取阈值进行植被范围的提取,设置近红外波段光谱的水稻的提取阈值进行水稻范围的提取。
进一步,所述步骤2中基于专家知识的决策树分类进行分类和水稻提取的过程包括:
通过对植被物候期的分析得出所述感兴趣区所有植被的植被指数随时间变化的先验信息,选取一种差距显著的植被指数,确定基于所述植被指数进行分类的阈值,使用决策树分类方法对所述感兴趣区进行分类。
进一步,所述步骤3中提取EVI均值指数图像的过程包括:
以所述水稻区域作为掩膜,提取不同时期的MODIS影像中的EVI指数图像,对所述EVI指数图像做重投影、去除空值和取有效值处理,提取得到EVI有效值指数图像;
以统计区域的矢量边界为掩膜,基于所述EVI有效值指数图像提取得到所述统计区域的水稻EVI均值图像。
进一步,所述步骤4中采用逐步回归方法构建多元线性的所述统计回归估算模型,并根据输出参数选择最佳模型。
进一步,所述统计回归模型的公式为yield=135165.591-146486.295*EVId209-49164.16*EVId177+22347.074*EVId161,yield表示因变量水稻产量统计数据,EVId209、EVId177和EVId161分别表示全年第209、177和161天的EVI均值。
进一步,所述步骤4之后还包括:
步骤5,使用水稻产量统计数据验证所述统计回归估算模型的精度并做相应的评估。
采用上述进一步方案的有益效果是:
1、针对农作物面积统计和产量预估等传统农业应用需求,本发明提出了一种引入遥感技术手段、实现农业遥感,同时协调农村与自然资源领域统计数据的综合性行业需求解决技术方案,为农村产业升级融合提供了应用示范,推动农村相关多种类产业的升级和融合发展。
2、针对水稻面积提取过程中单一分类方法准确率低的问题,本发明提出了一种结合基于人工目视解译的监督分类,基于规则的面向对象分类和基于专家知识的决策树分类方式,通过3种方式结合取相交区域的方式,减少误分类带来的误差,提高准确率。
3、针对传统水稻估产速度慢、工作量大、成本高的问题,本发明提出了一种基于遥感影像提取水稻作物信息的方式,能够更加快速、准确地获取到水稻植被信息;提出了一种通过建立数理统计学模型来预估产量的方式,大大减少了工作量,降低了成本。
4、针对高分辨率影像重访周期较长,难以获取针对特定时期的多时序影像,低分辨率影像提取作物精度较低的问题,本发明提出了一种结合使用高低分辨率影像的方式,使用高分辨率影像提取作物基础面积,使用中低分辨率影像提取多时序作物指数信息,同时减少来自大气和土壤噪音的影响,既解决了高分辨率影像用于多时序指数提取数据量不足的问题,又解决了低分辨率提取精度低的问题。
附图说明
图1为本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法的实施例的流程图;
图2为本发明实施例提供的采用最大似然法进行分类后的分类结果示意图;
图3为本发明实施例提供的采用最大似然法进行分类后的水稻提取结果示意图;
图4为本发明实施例提供的面向对象分类的水稻提取结果图;
图5为本发明实施例提供的决策树进行分类的水稻提取结果图;
图6为本发明实施例提供的三种分类结果相交区域示意图;
图7为本发明实施例提供的水稻区域EVI有效值示意图;
图8为本发明实施例提供的水稻EVI均值示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示为本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法的实施例的流程图,结合图1可知,该方法包括:
步骤1,对原始的卫星影像进行预处理操作。
该卫星影像可以为高分二号(GF-2)卫星影像,GF-2卫星是空间分辨率优于1米的民用光学遥感卫星,搭载有两台高分辨率1米全色、4米多光谱相机,星下点空间分辨率可达0.8米,具有亚米级空间分辨率、高定位精度和快速姿态机动能力等特点,有效地提升了卫星综合观测效能,主要用户为国土资源部、住房和城乡建设部、交通运输部和国家林业局等部门,同时还将为其他用户部门和有关区域提供示范应用服务。
MODIS(Moderate-resolution Imaging Spectroradiometer,中分辨率成像光谱仪)是美国宇航局研制大型空间遥感仪器,以了解全球气候的变化情况以及人类活动对气候的影响。
步骤2,在预处理后的卫星影像上选取感兴趣区,分别采用基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种方法对感兴趣区进行分类和水稻提取,分类后的类型包括:林地、水稻、水体和建设用地。
为了提高水稻面积的提取精度,同时结合基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种影像分类方法的提取结果,以确定水稻估产模型基础面积。使用高分辨率GF-2数据,结合基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类3种方法进行农作物分类,分类精度较高,满足精度需求。
步骤3,取三种分类方法得到的水稻提取结果的相交区域作为水稻区域,提取水稻区域的不同时期的MODIS影像中的EVI(Enhanced VegetationIndex增强植被指数)指数均值图像。
使用多时序MODIS数据提取水稻相关植被指数,数据覆盖的周期更广,同时提取的植被指数满足估产需求;进行作物基础估产,估产成本较小,精度符合实际需要,模型可用性较强。
步骤4,将水稻EVI均值作为自变量,水稻产量统计数据作为因变量,进行回归分析,构建水稻EVI均值与水稻产量统计数据之间的统计回归估算模型。
为构建不同光谱值运算得到的植被指数与水稻产量之间的统计回归估算模型,使用重访周期短的MODIS数据获取多时序的影像数据。为减少来自大气和土壤噪音的影响,稳定反映地表植被的特征,使用EVI作为提取的特征指数,EVI的计算公式为其中G为增益系数,NIR为近红外波段反射率,R为红光波段发射率,B为蓝光波段发射率,C1、C2为气溶胶阻抗系数,L为冠层背景调整因子。
本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法,使用高分辨率影像,通过不同提取特征、不同分类方法综合提取水稻区域,在此基础上,使用多时序影像来提取植被指数EVI,以描述植被在不同季节的差异。通过探究水稻植被指数与水稻产量之间的数学关系,建立估产的数学模型,能够更加快速、准确、有效地估测水稻的产量;结合高分辨率的GF-2数据、中低分辨率的MODIS数据进行估产,既解决了特定时期内高分辨率影像用于多时序指数提取数据量不足的问题,又解决了低分辨率提取精度低的问题。
实施例1
本发明提供的实施例1为本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法的实施例,本发明提供的一种结合卫星影像和MODIS数据的水稻估产方法的实施例中,首先基于高分辨率GF-2卫星影像、高清地图影像、部分实地考察水稻样本点和水稻种植标准农田数据,结合人工目视解译的监督分类、规则的面向对象分类和专家知识的决策树分类3种方法,综合确定地区的水稻种植范围。再根据水稻生育物候期,选取多时相多光谱MODIS13Q1影像产品,通过探究影像经过不同波段运算得到的植被指数与水稻产量之间的数学关系,实现利用遥感数据的作物产量监测和估算。最后在SPSS中对提取出的水稻区域EVI均值和县域水稻年产量进行相关性分析,获得显著相关性结果。
结合图1可知,该实施例包括:
步骤1,对原始的卫星影像进行预处理操作。
优选的,对卫星影像进行预处理操作包括:大气校正、影像融合、影像剪裁和影像拼接等。
在使用GF-2数据时,需要对原始GF-2影像做相应的影像预处理操作,对GF-2影像进行融合和大气校正操作,为后续提取壁纸植被指数做准备。
如图1所示,本发明以湖北省恩施土家族苗族自治州宣恩县为例,由当地水稻种植的物候期可知,7月上旬、10月初分别为早稻孕穗抽穗期和晚稻收获期,水稻生长茂盛。选择对应时期的GF-2影像作为原始数据,并对影像做大气校正和重采样后,使用宣恩县行政区矢量边界作为掩膜,对影像进行裁剪,拼接等操作。
步骤2,在预处理后的高分辨率(2m)GF-2影像上选取感兴趣区,分别采用基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种方法对感兴趣区进行分类和水稻提取,分类后的类型包括:林地、水稻、水体和建设用地。
优选的,选取感兴趣区包括:结合高清地图影像、实地考察水稻样本点和水稻种植标准农田数据,结合谷歌地球的显示影像确定地类的解译标志,在卫星影像上使用人工标注的方式选取感兴趣区域。
基于实地考察结果、利用遥感技术手段实现农作物面积统计和产量预估等传统农业应用的同时,协调农村农业与自然资源部门及领域统计数据,形成农村行业综合体应用,能够推动农村相关多种类产业的升级和融合发展。
具体的,采用基于人工目视解译的监督分类进行分类和水稻提取的过程包括:采用最大似然法对感兴趣区进行分类和水稻提取。
如图2和图3所示分别为本发明实施例提供的采用最大似然法进行分类后的分类结果示意图和水稻提取结果示意图。
基于规则的面向对象分类进行分类和水稻提取的过程包括:
根据光谱范围的差异建立水稻的提取规则,使用基于提取规则的面向对象分类方法对感兴趣区进行分类。
进一步的,通过对各类感兴趣区样本进行光谱统计分析可以得出,植被在可见光蓝色波段与山体、建设用地的光谱范围差异明显,水稻等农作物在近红外波段光谱范围整体大于山区植被。
对应的,该提取规则包括:设置在可见光蓝色波段光谱的植被的提取阈值进行植被范围的提取,设置近红外波段光谱的水稻的提取阈值进行水稻范围的提取。如下表1所示为基于光谱范围的水稻提取规则表:
表1:水稻提取规则表
如图4所示为本发明实施例提供的面向对象分类的水稻提取结果图。
基于专家知识的决策树分类进行分类和水稻提取的过程包括:
通过对植被物候期的分析得出感兴趣区所有植被的植被指数随时间变化的先验信息,选取一种差距显著的植被指数,例如RVI(Ratio Vegetation Index,比值植被指数),确定基于植被指数进行分类的阈值,使用决策树分类方法对感兴趣区进行分类。
具体实施中,基于对水稻物候期的分析可得出,7月上旬、10月初当地种植的农作物分别为水稻和玉米,且10月初玉米已完成抽穗,植被指数迅速下降,同时RVI对高植被地区覆盖敏感,且山地植被的RVI指数均值较农作物区域更大,RVI的计算公式为其中NIR为近红外波段反射率,R为红光波段发射率,在对影像数据进行大气校正的基础上,通过阈值划分RVI的取值范围。如下表2所示为基于决策树分类规则表:
表2:水稻分类规则表
如图5所示为本发明实施例提供的决策树进行分类的水稻提取结果图。
如图6所示为本发明实施例提供的三种分类结果相交区域示意图,如图6所示,计算区域水稻面积,与水稻面积统计数据的误差处于合理误差范围内,可作为水稻估产的基础水稻面积。
步骤3,取三种分类方法得到的水稻提取结果的相交区域作为水稻区域,提取水稻区域的不同时期的MODIS影像中的EVI指数均值图像。
MODIS影像使用MODIS13Q1(250m16天合成)产品,相比原粗颗粒估产更具可行性。因为村县级地区分散稻田面积通常范围不大,粗颗粒估产在村县级地区可行性较低;同时,MODIS09A1获取EVI值需要通过多波段修改处理后进行组合和公式计算,MODIS13Q1直接获取第2波段并修改处理即为EVI有效值。
优选的,提取EVI均值指数图像的过程包括:
以水稻区域作为掩膜,提取不同时期的MODIS影像中的EVI指数图像,对EVI指数图像做重投影、去除空值和取有效值处理,提取得到EVI有效值指数图像。
以统计区域的矢量边界为掩膜,基于EVI有效值指数图像提取得到统计区域的水稻EVI均值图像。
如图7和图8所示为本发明实施例提供的水稻区域EVI有效值示意图和水稻EVI均值示意图。
步骤4,将水稻EVI均值作为自变量,水稻产量统计数据作为因变量,进行回归分析,构建水稻EVI均值与水稻产量统计数据之间的统计回归估算模型。
优选的,采用逐步回归方法构建多元线性的统计回归估算模型,并根据输出参数选择最佳模型。如下表3所示为统计回归估算模型的回归参数取值表:
表3:统计回归估算模型的回归参数取值表
由表3可以得到,统计回归模型的公式为yield=135165.591-146486.295*EVId209-49164.16*EVId177+22347.074*EVId161,yield表示因变量水稻产量统计数据,EVId209、EVId177和EVId161分别表示全年第209、177和161天的EVI均值。
对于具有因变量yield的最终模型,拟合效果较好,方差膨胀因子VIF值均接近于1,自变量之间基本不存在共线性问题,同时以2012年-2017年全年第161天(6月11日左右)、第177天(6月27日左右)、第209天(7月29日左右)作为自变量,分别对应早稻孕穗期、早稻黄熟期和晚稻黄熟期,植被指数较高,符合水稻估产模型的现实理论意义。
进一步优选的,步骤4之后还包括:
步骤5,使用水稻产量统计数据验证统计回归估算模型的精度并做相应的评估。
本发明提供的实施例中,对生长期间每连续16天的EVI值进行逐步回归的模型拟合,以数理统计规律确定最具影响因素的自变量,最终确定仅以3天EVI值为自变量的统计模型;具体的,使用2018-2020年全年第161天(6月11日左右)、第177天(6月27日左右)、第209天(7月29日左右)共9期MODIS影像数据提取EVI均值并代入统计回归估算模型的公式,与2018-2020年的水稻产量统计数据对比,对比结果如表4所示,模型精度得以验证。
表4:水稻估产量与实际产量对比表
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种结合卫星影像和MODIS数据的水稻估产方法,其特征在于,所述方法包括:
步骤1,对原始的卫星影像进行预处理操作;
步骤2,在预处理后的所述卫星影像上选取感兴趣区,分别采用基于人工目视解译的监督分类、基于规则的面向对象分类和基于专家知识的决策树分类三种方法对所述感兴趣区进行分类和水稻提取,分类后的类型包括:林地、水稻、水体和建设用地;
步骤3,取三种分类方法得到的水稻提取结果的相交区域作为水稻区域,提取所述水稻区域的不同时期的MODIS影像中的EVI指数均值图像;
步骤4,将水稻EVI均值作为自变量,水稻产量统计数据作为因变量,进行回归分析,构建所述水稻EVI均值与所述水稻产量统计数据之间的统计回归估算模型。
2.根据权利要求1所述的方法,其特征在于,所述步骤1中对所述卫星影像进行预处理操作包括:大气校正、影像融合、影像剪裁和影像拼接。
3.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤2中选取所述感兴趣区包括:结合高清地图影像、实地考察水稻样本点和水稻种植标准农田数据,在所述卫星影像上使用人工标注的方式选取所述感兴趣区域。
4.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤2中采用基于人工目视解译的监督分类进行分类和水稻提取的过程包括:采用最大似然法对所述感兴趣区进行分类和水稻提取。
5.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤2中基于规则的面向对象分类进行分类和水稻提取的过程包括:
根据光谱范围的差异建立水稻的提取规则,使用基于所述提取规则的面向对象分类方法对所述感兴趣区进行分类;
所述提取规则包括:设置在可见光蓝色波段光谱的植被的提取阈值进行植被范围的提取,设置近红外波段光谱的水稻的提取阈值进行水稻范围的提取。
6.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤2中基于专家知识的决策树分类进行分类和水稻提取的过程包括:
通过对植被物候期的分析得出所述感兴趣区所有植被的植被指数随时间变化的先验信息,选取一种差距显著的植被指数,确定基于所述植被指数进行分类的阈值,使用决策树分类方法对所述感兴趣区进行分类。
7.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤3中提取EVI均值指数图像的过程包括:
以所述水稻区域作为掩膜,提取不同时期的MODIS影像中的EVI指数图像,对所述EVI指数图像做重投影、去除空值和取有效值处理,提取得到EVI有效值指数图像;
以统计区域的矢量边界为掩膜,基于所述EVI有效值指数图像提取得到所述统计区域的水稻EVI均值图像。
8.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤4中采用逐步回归方法构建多元线性的所述统计回归估算模型,并根据输出参数选择最佳模型。
9.根据权利要求1所述的水稻估产方法,其特征在于,所述统计回归模型的公式为yield=135165.591-146486.295*EVId209-49164.16*EVId177+22347.074*EVId161,yield表示因变量水稻产量统计数据,EVId209、EVId177和EVId161分别表示全年第209、177和161天的EVI均值。
10.根据权利要求1所述的水稻估产方法,其特征在于,所述步骤4之后还包括:
步骤5,使用水稻产量统计数据验证所述统计回归估算模型的精度并做相应的评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110488775.8A CN113313671A (zh) | 2021-04-30 | 2021-04-30 | 一种结合卫星影像和modis数据的水稻估产方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110488775.8A CN113313671A (zh) | 2021-04-30 | 2021-04-30 | 一种结合卫星影像和modis数据的水稻估产方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113313671A true CN113313671A (zh) | 2021-08-27 |
Family
ID=77371560
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110488775.8A Pending CN113313671A (zh) | 2021-04-30 | 2021-04-30 | 一种结合卫星影像和modis数据的水稻估产方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113313671A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115578644A (zh) * | 2022-08-08 | 2023-01-06 | 江苏省水利科学研究院 | 基于多特征机器学习的平原水网水葫芦提取方法及装置 |
CN117197651A (zh) * | 2023-07-24 | 2023-12-08 | 移动广播与信息服务产业创新研究院(武汉)有限公司 | 一种结合边缘检测与语义分割的田块提取方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858971A (zh) * | 2010-06-02 | 2010-10-13 | 浙江大学 | 一种基于modis数据的水稻单产遥感估算方法 |
CN103761447A (zh) * | 2014-02-18 | 2014-04-30 | 中国科学院遥感与数字地球研究所 | 用于地市及县两级粮食补贴核算的种植面积遥感估算方法 |
CN105740759A (zh) * | 2016-01-15 | 2016-07-06 | 武汉珈和科技有限公司 | 基于多时相数据中特征提取的中稻信息决策树分类方法 |
CN105787457A (zh) * | 2016-03-08 | 2016-07-20 | 浙江工商大学 | 一种modis卫星集成dem提高植被分类遥感精度的估算方法 |
CA3007202A1 (en) * | 2015-12-02 | 2017-06-08 | The Climate Corporation | Forecasting field level crop yield during a growing season |
-
2021
- 2021-04-30 CN CN202110488775.8A patent/CN113313671A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101858971A (zh) * | 2010-06-02 | 2010-10-13 | 浙江大学 | 一种基于modis数据的水稻单产遥感估算方法 |
CN103761447A (zh) * | 2014-02-18 | 2014-04-30 | 中国科学院遥感与数字地球研究所 | 用于地市及县两级粮食补贴核算的种植面积遥感估算方法 |
CA3007202A1 (en) * | 2015-12-02 | 2017-06-08 | The Climate Corporation | Forecasting field level crop yield during a growing season |
CN105740759A (zh) * | 2016-01-15 | 2016-07-06 | 武汉珈和科技有限公司 | 基于多时相数据中特征提取的中稻信息决策树分类方法 |
CN105787457A (zh) * | 2016-03-08 | 2016-07-20 | 浙江工商大学 | 一种modis卫星集成dem提高植被分类遥感精度的估算方法 |
Non-Patent Citations (1)
Title |
---|
刘珊: "基于农业遥感技术的水稻种植信息方法的研究——以湘潭市为例", 《中国优秀硕士学位论文全文数据库 农业科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115578644A (zh) * | 2022-08-08 | 2023-01-06 | 江苏省水利科学研究院 | 基于多特征机器学习的平原水网水葫芦提取方法及装置 |
CN115578644B (zh) * | 2022-08-08 | 2023-04-18 | 江苏省水利科学研究院 | 基于多特征机器学习的平原水网水葫芦提取方法及装置 |
CN117197651A (zh) * | 2023-07-24 | 2023-12-08 | 移动广播与信息服务产业创新研究院(武汉)有限公司 | 一种结合边缘检测与语义分割的田块提取方法及系统 |
CN117197651B (zh) * | 2023-07-24 | 2024-03-29 | 移动广播与信息服务产业创新研究院(武汉)有限公司 | 一种结合边缘检测与语义分割的田块提取方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Improved estimation of rice aboveground biomass combining textural and spectral analysis of UAV imagery | |
Gašparović et al. | An automatic method for weed mapping in oat fields based on UAV imagery | |
Darvishzadeh et al. | Mapping leaf chlorophyll content from Sentinel-2 and RapidEye data in spruce stands using the invertible forest reflectance model | |
CN109345555B (zh) | 基于多时相多源遥感数据进行水稻识别的方法 | |
CN110472184B (zh) | 一种基于Landsat遥感数据的多云雨雾地区水稻识别方法 | |
Johnson | Temporal stability of an NDVI‐LAI relationship in a Napa Valley vineyard | |
Onojeghuo et al. | Applications of satellite ‘hyper-sensing’in Chinese agriculture: Challenges and opportunities | |
Raj et al. | Precision agriculture and unmanned aerial Vehicles (UAVs) | |
Nuarsa et al. | Spectral characteristics and mapping of rice plants using multi-temporal Landsat data | |
CN109685081B (zh) | 一种遥感提取撂荒地的联合变化检测方法 | |
Sharifi | Estimation of biophysical parameters in wheat crops in Golestan province using ultra-high resolution images | |
CN114821362B (zh) | 一种基于多源数据的水稻种植面积提取方法 | |
CN111028096A (zh) | 一种天、空、地一体化数据融合的系统和方法 | |
White et al. | Determining a robust indirect measurement of leaf area index in California vineyards for validating remote sensing-based retrievals | |
Ahmad | A review of remote sensing data change detection: Comparison of Faisalabad and Multan Districts, Punjab Province, Pakistan | |
CN113313671A (zh) | 一种结合卫星影像和modis数据的水稻估产方法 | |
Dong et al. | Using RapidEye imagery to identify within-field variability of crop growth and yield in Ontario, Canada | |
Saura et al. | Mapping multispectral Digital Images using a Cloud Computing software: applications from UAV images | |
Zhou et al. | Remote sensing of regional-scale maize lodging using multitemporal GF-1 images | |
Amankulova et al. | Comparison of PlanetScope, Sentinel-2, and landsat 8 data in soybean yield estimation within-field variability with random forest regression | |
Qu et al. | Monitoring lodging extents of maize crop using multitemporal GF-1 images | |
Naqvi et al. | Remote estimation of wheat yield based on vegetation indices derived from time series data of Landsat 8 imagery. | |
Ouyang et al. | UAV and ground-based imagery analysis detects canopy structure changes after canopy management applications | |
Mashonganyika et al. | Mapping of winter wheat using sentinel-2 NDVI data. a case of Mashonaland central province in Zimbabwe | |
Sahoo et al. | Drone remote sensing of wheat N using hyperspectral sensor and machine learning |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210827 |