CN108982369A - 融合gf-1wfv和modis数据的地块尺度作物长势监测方法 - Google Patents

融合gf-1wfv和modis数据的地块尺度作物长势监测方法 Download PDF

Info

Publication number
CN108982369A
CN108982369A CN201810401690.XA CN201810401690A CN108982369A CN 108982369 A CN108982369 A CN 108982369A CN 201810401690 A CN201810401690 A CN 201810401690A CN 108982369 A CN108982369 A CN 108982369A
Authority
CN
China
Prior art keywords
crop
data
ndvi
1wfv
modis
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
Application number
CN201810401690.XA
Other languages
English (en)
Other versions
CN108982369B (zh
Inventor
黄健熙
卓文
苏伟
刘峻明
李俐
刘哲
张超
张晓东
朱德海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jinzhinong (Beijing) Risk Management Technology Co.,Ltd.
Original Assignee
China Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Agricultural University filed Critical China Agricultural University
Priority to CN201810401690.XA priority Critical patent/CN108982369B/zh
Publication of CN108982369A publication Critical patent/CN108982369A/zh
Application granted granted Critical
Publication of CN108982369B publication Critical patent/CN108982369B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N2021/1793—Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Image Processing (AREA)

Abstract

本发明属于农业遥感领域，涉及一种融合GF‑1 WFV和MODIS数据的地块尺度作物长势监测方法，具体为：对历史几年作物生育期内MODIS反射率数据和GF‑1 WFV影像数据预处理；基于卡尔曼滤波算法融合两种遥感数据；利用融合后的红光波段与近红外波段，计算历史几年NDVI时间序列影像并平均；目标年份监测时期的GF‑1NDVI影像与对应时期的历史平均NDVI影像做差值判断长势等级，逐个作物格网单元运行获得监测结果指导作物生产。本发明克服了由于雨云影响难以获得作物生育期内中等空间分辨率多年平均NDVI基准曲线的问题，通过MODIS和WFV的数据融合方法，为精准长势监测提供了基准参考曲线。

Description

融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法
技术领域
本发明属于农业遥感领域，具体涉及一种融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法。
背景技术
作物长势的实时监测对作物的产量预报和农田精细管理具有重要意义。传统的作物长势监测方法多为实地调查，不仅费时费力而且很难获得大面积的作物长势信息。基于遥感技术的作物长势监测能够获得大面积地表信息，还可以实时动态监测。目前基于遥感技术的长势监测方法多是基于遥感植被指数对作物长势进行识别，这对遥感数据的时间分辨率和空间分辨率都有很高的要求。然而，目前遥感影像较少有同时满足高空间分辨率高时间分辨率的要求。鉴于此，将高空间分辨率的影像与高时间分辨率的影像进行融合，生成一种同时具有高时间分辨率和高空间分辨率的时间序列影像，克服了高空间分辨率影像受云雨影像数据量不足的局限性。并通过计算多年NDVI影像作为长势基准，实现了对作物的全生育期长势监测。
发明内容
本发明提供一种融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法。
本发明提供一种融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法，具体步骤如下：
S1、收集研究区的目标年份之前历史几年作物生育期内MODIS(搭载在terra和aqua卫星上的传感器)反射率数据和GF-1WFV影像数据(高分一号卫星GF-1搭载的WFV相机的影像数据)，并对上述两种遥感数据进行预处理；所述MODIS反射率数据的时间分辨率为每天、空间分辨率为500米；所述GF-1WFV影像数据的时间分辨率为4天、空间分辨率为16米；
S2、基于卡尔曼滤波(Kalman filtering，简写为KF)算法，融合MODIS反射率数据和GF-1WFV影像数据，得到时间分辨率为每天，空间分辨率为16米的反射率数据；
S3、利用融合后的反射率数据中的红光波段与近红外波段，计算历史几年分辨率为每天、空间分辨率为16米的NDVI(归一化植被指数)时间序列影像；
S4、将S3计算得到的NDVI时间序列影像进行平均，得到时间分辨率为每天、空间分辨率为16米的历史几年平均NDVI时间序列影像；
S5、以S4计算得到的历史几年平均NDVI时间序列影像为长势基准，将目标年份监测时期的GF-1NDVI影像(包括在GF-1WFV影像数据内)与对应时期的历史平均NDVI影像做差值，并以此作为判断作物长势等级的标准；
S6、逐个作物格网单元运行S5，获得目标年份作物长势监测结果，指导作物生产。
其中步骤S1之前还可包括如下确定作物地块的步骤：
收集研究区Landsat8OLI(Landsat 8卫星上荷载的陆地成像仪Operational LandImager)和Sentinel-2(Sentinel-2卫星)光学遥感数据，并进行农作物分类，获得待测作物种植地块。
其中，步骤S1中所述目标年份之前历史几年，优选目标年份之前历史5年。
其中，步骤S1所述预处理为对GF-1WFV影像数据进行几何纠正、辐射定标、大气校正等预处理，并通过寻找两种遥感数据的同名控制点，完成MODIS反射率数据和GF-1WFV影像数据的几何配准。
步骤S2所述卡尔曼滤波的实现需要定义观测变量为GF-1WFV影像数据，并转换模型为两个子模型的集合：
假设同一状态下，表征农作物反射率在生长季的动态变化趋势的子模型1和刻画GF-1WFV影像数据和MODIS反射率数据之间的关系的子模型2的估计结果均相互独立且服从高斯分布，在每一个时间步骤，卡尔曼算法中的时间更新模型是它们共同的估算结果，计算方式如下：
分别表示在k状态下，子模型1、子模型2和组合模型的先验估计值，分别表示对应的先验方差，将子模型2的不确定度设定为线性回归的标准误差；
其中，子模型1提供了表征农作物反射率在生长季的动态变化趋势，选取待测农作物纯度大于50％的像元；线性状态转换模型为每天在MODIS反射率数据中提取构建k-1和k时刻每个时间状态的植被生长变化轨迹线性关系；子模型1将曲线定义为分段回归曲线；每一段通过连续的MODIS反射率线性回归得到，其中连续的MODIS反射率通过在计算机能够保持计算效率的影像上大量随机选取像元构建得到，其具体公式表达为：
是子模型1中状态k时的先验估计值，是状态k-1的后验估计值，并设定每一步线性回归模型的标准误差作为子模型1每一步对应的不确定性；
子模型2刻画GF-1WFV影像数据和MODIS反射率数据之间的关系，从而实现对子模型1周期性曲线的约束；子模型1中的斜率和截距则被用于获取k状态下子模型2的先验估计值，其具体公式表达为：
在公式(4)中，是子模型2中k状态下的先验估计值，c和d是k状态下MODIS影像像元和GF-1WFV影像像元之间建立线性关系的斜率和截距。
其中，步骤S3所述NDVI的计算公式为：
式中ρnir表示近红外波段的反射率值，ρred表示红光波段的反射率值。
其中，步骤S4所述进行平均，采用公式(6)进行计算：
式中，NDVIi,mean表示第i天历史几年NDVI平均值，NDVIi,j表示第j年第i天的NDVI值。
步骤S5所述判断作物长势等级的分级标准为：δ>0.25表示长势好，为5级；0.25>δ>0.025表示长势较好，为4级；0.025>δ>-0.025表示长势持平，为3级；-0.025>δ>-0.25表示长势较差，为2级；δ<-0.25表示长势差，为1级。
步骤S6所述逐个作物格网单元运行S5，为利用S5建立的长势监测分级标准，对每个像元进行重新赋值1-5级，从而获得目标年份作物长势监测结果。
其中，所述作物优选为冬小麦、水稻等任一种大田作物。
本发明还提供所述融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法在指导农作物生产中的应用。
本发明与现有技术相比，长势监测结果准确度提高了23％。
附图说明
图1为本发明实施例1对水稻实施融合GF-1WFV和MODIS数据的地块尺度水稻长势监测方法的流程示意图。
图2为本发明实施例1以实施融合GF-1WFV和MODIS数据的地块尺度水稻长势监测方法对黑龙家七星农场水稻长势监测得到的2016年7月25日的监测图。
具体实施方式
下面结合实施例，对本范明的具体实施方式做进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
实施例1
以本发明融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法针对水稻进行估产的流程示意图参见附图1。
选择黑龙江省作为研究区域，该区域地处东经121°11′-135°05′，北纬43°26′-53°33′之间。研究区总面积47.3万平方公里，地形以平原和山地为主，耕地占总面积的25％以上，属温带大陆性季风气候，年日照时数2400-2800h，年平均降水量400-600mm。
S1、收集研究区的目标年份之前历史5年作物生育期内MODIS反射率数据和GF-1WFV影像数据，并对上述两种遥感数据进行预处理；所述MODIS反射率数据的时间分辨率为每天、空间分辨率为500米；所述GF-1WFV影像数据的时间分辨率为4天、空间分辨率为16米；
所述预处理为对GF-1WFV影像数据进行几何纠正、辐射定标、大气校正等预处理，并通过寻找两种遥感数据的同名控制点，完成MODIS反射率数据和GF-1WFV影像数据的几何配准。
S2、基于卡尔曼滤波算法，融合MODIS反射率数据和GF-1WFV影像数据，得到时间分辨率为每天，空间分辨率为16米的反射率数据；
步骤S2所述卡尔曼滤波的实现需要定义观测变量为GF-1WFV影像数据，并转换模型为两个子模型的集合：
假设同一状态下，表征农作物反射率在生长季的动态变化趋势的子模型1和刻画GF-1WFV影像数据和MODIS反射率数据之间的关系的子模型2的估计结果均相互独立且服从高斯分布，在每一个时间步骤，卡尔曼算法中的时间更新模型是它们共同的估算结果，计算方式如下：
分别表示在k状态下，子模型1、子模型2和组合模型的先验估计值，分别表示对应的先验方差，将子模型2的不确定度设定为线性回归的标准误差；
其中，子模型1提供了表征农作物反射率在生长季的动态变化趋势，选取待测农作物纯度大于50％的像元；线性状态转换模型为每天在MODIS反射率数据中提取构建k-1和k时刻每个时间状态的植被生长变化轨迹线性关系；子模型1将曲线定义为分段回归曲线；每一段通过连续的MODIS反射率线性回归得到，其中连续的MODIS反射率通过在计算机能够保持计算效率的影像上大量随机选取像元构建得到，其具体公式表达为：
是子模型1中状态k时的先验估计值，是状态k-1的后验估计值，并设定每一步线性回归模型的标准误差作为子模型1每一步对应的不确定性；
子模型2刻画GF-1WFV影像数据和MODIS反射率数据之间的关系，从而实现对子模型1周期性曲线的约束；子模型1中的斜率和截距则被用于获取k状态下子模型2的先验估计值，其具体公式表达为：
在公式(4)中，是子模型2中k状态下的先验估计值，c和d是k状态下MODIS影像像元和GF-1WFV影像像元之间建立线性关系的斜率和截距。
S3、利用融合后的反射率数据中的红光波段与近红外波段，计算历史5年分辨率为每天、空间分辨率为16米的NDVI(植被覆盖指数)时间序列影像；
NDVI的计算公式为：
式中ρnir表示近红外波段的反射率值，ρred表示红光波段的反射率值。
S4、将S3计算得到的NDVI时间序列影像进行平均，得到时间分辨率为每天、空间分辨率为16米的历史5年平均NDVI时间序列影像；
所述进行平均，采用公式(6)进行计算：
式中，NDVIi,mean表示第i天历史5年NDVI平均值，NDVIi,j表示第j年第i天的NDVI值。
S5、以S4计算得到的历史5年平均NDVI时间序列影像为长势基准，将目标年份监测时期的GF-1NDVI影像(包括在GF-1WFV影像数据内)与对应时期的历史平均NDVI影像做差值，并以此作为判断作物长势等级的标准；
所述判断作物长势等级的分级标准为：δ>0.25表示长势好，为5级；0.25>δ>0.025表示长势较好，为4级；0.025>δ>-0.025表示长势持平，为3级；-0.025>δ>-0.25表示长势较差，为2级；δ<-0.25表示长势差，为1级。
S6、逐个作物格网单元运行S5，利用S5建立的长势监测分级标准，对每个像元进行重新赋值1-5级，从而获得目标年份水稻长势监测结果，指导生产。
本实施例得到的2016年7月25日黑龙江省七星农场这以地块的水稻长势监测图见图2。
虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

Claims (10)

1.一种融合GF-1WFV和MODIS数据的地块尺度作物长势监测方法，其特征在于，具体步骤如下：
S1、收集研究区的目标年份之前历史几年作物生育期内MODIS反射率数据和GF-1WFV影像数据，并对上述两种遥感数据进行预处理；所述MODIS反射率数据的时间分辨率为每天、空间分辨率为500米；所述GF-1WFV影像数据的时间分辨率为4天、空间分辨率为16米；
S2、基于卡尔曼滤波算法，融合MODIS反射率数据和GF-1WFV影像数据，得到时间分辨率为每天，空间分辨率为16米的反射率数据；
S3、利用融合后的反射率数据中的红光波段与近红外波段，计算历史几年分辨率为每天、空间分辨率为16米的NDVI时间序列影像；
S4、将S3计算得到的NDVI时间序列影像进行平均，得到时间分辨率为每天、空间分辨率为16米的历史几年平均NDVI时间序列影像；
S5、以S4计算得到的历史几年平均NDVI时间序列影像为长势基准，将目标年份监测时期的GF-1NDVI影像与对应时期的历史平均NDVI影像做差值，并以此作为判断作物长势等级的标准；
S6、逐个作物格网单元运行S5，获得目标年份作物长势监测结果，指导作物生产。
2.如权利要求1所述的方法，其特征在于，步骤S1之前还包括如下确定作物地块的步骤：
收集研究区Landsat8OLI和Sentinel-2光学遥感数据，并进行农作物分类，获得待测作物种植地块。
3.如权利要求1或2任一项所述的方法，其特征在于，步骤S1中所述目标年份之前历史几年，为目标年份之前历史5年。
4.如权利要求1或2任一项所述的方法，其特征在于，步骤S1所述预处理为对GF-1WFV影像数据进行几何纠正、辐射定标、大气校正，并通过寻找两种遥感数据的同名控制点，完成MODIS反射率数据和GF-1WFV影像数据的几何配准。
5.如权利要求1或2任一项所述的方法，其特征在于，步骤S2所述卡尔曼滤波的实现需要定义观测变量为GF-1WFV影像数据，并转换模型为两个子模型的集合：
假设同一状态下，表征农作物反射率在生长季的动态变化趋势的子模型1和刻画GF-1WFV影像数据和MODIS反射率数据之间的关系的子模型2的估计结果均相互独立且服从高斯分布，在每一个时间步骤，卡尔曼算法中的时间更新模型是它们共同的估算结果，计算方式如下：
分别表示在k状态下，子模型1、子模型2和组合模型的先验估计值，分别表示对应的先验方差，将子模型2的不确定度设定为线性回归的标准误差；
计算公式如下：
是子模型1中状态k时的先验估计值，是状态k-1的后验估计值，并设定每一步线性回归模型的标准误差作为子模型1每一步对应的不确定性；
计算公式如下：
是子模型2中k状态下的先验估计值，c和d是k状态下MODIS影像像元和GF-1WFV影像像元之间建立线性关系的斜率和截距。
6.如权利要求1或2任一项所述的方法，其特征在于，步骤S3所述NDVI的计算公式为：
式中ρnir表示近红外波段的反射率值，ρred表示红光波段的反射率值。
7.如权利要求6所述的方法，其特征在于，步骤S4所述进行平均，采用公式(6)进行计算：
式中，NDVIi,mean表示第i天历史几年NDVI平均值，NDVIi,j表示第j年第i天的NDVI值。
8.如权利要求1或2任一项所述的方法，其特征在于，步骤S5所述判断作物长势等级的分级标准为：δ>0.25表示长势好，为5级；0.25>δ>0.025表示长势较好，为4级；0.025>δ>-0.025表示长势持平，为3级；-0.025>δ>-0.25表示长势较差，为2级；δ<-0.25表示长势差，为1级。
9.如权利要求1或2任一项所述的方法，其特征在于，步骤S6所述逐个作物格网单元运行S5，为利用S5建立的长势监测分级标准，对每个像元进行重新赋值1-5级，从而获得目标年份作物长势监测结果。
10.权利要求1-9任一项所述融合GF-1 WFV和MODIS数据的地块尺度作物长势监测方法在指导农作物生产中的应用。
CN201810401690.XA 2018-04-28 2018-04-28 融合gf-1wfv和modis数据的地块尺度作物长势监测方法 Active CN108982369B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810401690.XA CN108982369B (zh) 2018-04-28 2018-04-28 融合gf-1wfv和modis数据的地块尺度作物长势监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810401690.XA CN108982369B (zh) 2018-04-28 2018-04-28 融合gf-1wfv和modis数据的地块尺度作物长势监测方法

Publications (2)

Publication Number Publication Date
CN108982369A true CN108982369A (zh) 2018-12-11
CN108982369B CN108982369B (zh) 2020-09-01

Family

ID=64542360

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810401690.XA Active CN108982369B (zh) 2018-04-28 2018-04-28 融合gf-1wfv和modis数据的地块尺度作物长势监测方法

Country Status (1)

Country Link
CN (1) CN108982369B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109754125A (zh) * 2019-01-18 2019-05-14 中国农业大学 基于作物模型、历史和气象预报数据的作物产量预报方法
CN109765187A (zh) * 2019-01-27 2019-05-17 中国农业科学院农业资源与农业区划研究所 一种稻虾共作空间分布信息获取方法
CN110909821A (zh) * 2019-12-03 2020-03-24 中国农业科学院农业资源与农业区划研究所 基于作物参考曲线进行高时空分辨率植被指数数据融合的方法
CN111521563A (zh) * 2020-03-20 2020-08-11 航天信德智图（北京）科技有限公司 一种基于gf-1与modis时空融合的水稻提取方法
CN111832506A (zh) * 2020-07-20 2020-10-27 大同煤矿集团有限责任公司 一种基于长时序植被指数的重建植被遥感判别方法
CN112749597A (zh) * 2019-10-31 2021-05-04 北京珈和遥感科技有限公司 基于gf1高分辨率影像的作物区域单产遥感估测方法
CN113779291A (zh) * 2021-09-13 2021-12-10 平安科技（深圳）有限公司 一种样本标注方法、系统、介质及电子终端

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011102520A1 (ja) * 2010-02-22 2011-08-25 株式会社パスコ 水稲収量予測モデル生成方法、及び水稲収量予測方法
CN103955860A (zh) * 2014-04-17 2014-07-30 中国农业大学 一种基于集合卡尔曼滤波同化的区域作物估产方法
CN104134095A (zh) * 2014-04-17 2014-11-05 中国农业大学 一种基于尺度转换和数据同化的农作物产量估测方法
CN104915674A (zh) * 2014-10-24 2015-09-16 北京师范大学 Landsat8和MODIS融合构建高时空分辨率数据识别秋粮作物的方法
CN106600434A (zh) * 2016-10-18 2017-04-26 河南省农业科学院农业经济与信息研究所 基于作物模型与同化技术的作物长势遥感监测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011102520A1 (ja) * 2010-02-22 2011-08-25 株式会社パスコ 水稲収量予測モデル生成方法、及び水稲収量予測方法
CN103955860A (zh) * 2014-04-17 2014-07-30 中国农业大学 一种基于集合卡尔曼滤波同化的区域作物估产方法
CN104134095A (zh) * 2014-04-17 2014-11-05 中国农业大学 一种基于尺度转换和数据同化的农作物产量估测方法
CN104915674A (zh) * 2014-10-24 2015-09-16 北京师范大学 Landsat8和MODIS融合构建高时空分辨率数据识别秋粮作物的方法
CN106600434A (zh) * 2016-10-18 2017-04-26 河南省农业科学院农业经济与信息研究所 基于作物模型与同化技术的作物长势遥感监测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王利民等: "GF-1和MODIS影像冬小麦长势监测指标NDVI的对比", 《作物学报》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109754125A (zh) * 2019-01-18 2019-05-14 中国农业大学 基于作物模型、历史和气象预报数据的作物产量预报方法
CN109765187A (zh) * 2019-01-27 2019-05-17 中国农业科学院农业资源与农业区划研究所 一种稻虾共作空间分布信息获取方法
CN109765187B (zh) * 2019-01-27 2021-01-29 中国农业科学院农业资源与农业区划研究所 一种稻虾共作空间分布信息获取方法
CN112749597A (zh) * 2019-10-31 2021-05-04 北京珈和遥感科技有限公司 基于gf1高分辨率影像的作物区域单产遥感估测方法
CN110909821A (zh) * 2019-12-03 2020-03-24 中国农业科学院农业资源与农业区划研究所 基于作物参考曲线进行高时空分辨率植被指数数据融合的方法
CN111521563A (zh) * 2020-03-20 2020-08-11 航天信德智图（北京）科技有限公司 一种基于gf-1与modis时空融合的水稻提取方法
CN111832506A (zh) * 2020-07-20 2020-10-27 大同煤矿集团有限责任公司 一种基于长时序植被指数的重建植被遥感判别方法
CN111832506B (zh) * 2020-07-20 2023-10-17 中国矿业大学 一种基于长时序植被指数的重建植被遥感判别方法
CN113779291A (zh) * 2021-09-13 2021-12-10 平安科技（深圳）有限公司 一种样本标注方法、系统、介质及电子终端

Also Published As

Publication number Publication date
CN108982369B (zh) 2020-09-01

Similar Documents

Publication Publication Date Title
CN108982369A (zh) 融合gf-1wfv和modis数据的地块尺度作物长势监测方法
Zhong et al. Phenology-based crop classification algorithm and its implications on agricultural water use assessments in California’s Central Valley
Ma et al. Estimating regional winter wheat yield by assimilation of time series of HJ-1 CCD NDVI into WOFOST–ACRM model with Ensemble Kalman Filter
Zhou et al. Perspective of Chinese GF-1 high-resolution satellite data in agricultural remote sensing monitoring
Hunsaker et al. Wheat irrigation management using multispectral crop coefficients: I. Crop evapotranspiration prediction
CN106643529A (zh) 基于无人机影像的山区农作物生长高度快速测量的方法
Yun et al. Framework of SAGI agriculture remote sensing and its perspectives in supporting national food security
CN115376016A (zh) 一种基于植被水分指数与蒸散发结合的稻田实际灌溉面积识别方法
Kavvadias et al. Precision Agriculture-Comparison and Evaluation of Innovative Very High Resolution (UAV) and LandSat Data.
Wang et al. A scaling-based method for the rapid retrieval of FPAR from fine-resolution satellite data in the remote-sensing trend-surface framework
CN107273797A (zh) 基于水体指数变异系数的水稻亚像元识别方法
Sari et al. Detecting rice phenology in paddy fields with complex cropping pattern using time series MODIS data
CN111751286B (zh) 基于变化检测算法的土壤水分提取方法
Kulkarni et al. Estimating LAI of rice using NDVI derived from MODIS surface reflectance
CN116124709B (zh) 基于叶绿素相对含量的冬小麦旱情无人机监测判别方法
Nouri et al. Estimation of urban evapotranspiration through vegetation indices using worldview2 satellite remote sensing images
Tofigh et al. A comparison of actual evapotranspiration estimates based on Remote Sensing approaches with a classical climate data driven method
Kumari et al. Estimation of crop water requirement in rice-wheat system from multi-temporal AWIFS satellite data
Houston et al. Use of satellite data in agricultural surveys
Fuentes-Peñailillo et al. Use of a natural language processing bot for agricultural water management
Cheema et al. Quantification of land use changes in complex cropping of irrigated Indus basin, Pakistan using MODIS vegetation time series data.
Reyes-González et al. Assessing accuracy of vegetation index method to estimate actual evapotranspiration
AU2021100538A4 (en) Crop Health Monitoring System Using IoT and Machine Learning
AU2021107439A4 (en) A Novel method- Crop Health Monitoring System Using IoT and Machine Learning
CN115855841B (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
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220324

Address after: 100005 room 07-162, 101, floor 7, No. 219, Wangfujing Street, Dongcheng District, Beijing

Patentee after: Jinzhinong (Beijing) Risk Management Technology Co.,Ltd.

Address before: 100193 No. 2 Old Summer Palace West Road, Beijing, Haidian District

Patentee before: CHINA AGRICULTURAL University