CN109766871A - 一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 - Google Patents
一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 Download PDFInfo
- Publication number
- CN109766871A CN109766871A CN201910094028.9A CN201910094028A CN109766871A CN 109766871 A CN109766871 A CN 109766871A CN 201910094028 A CN201910094028 A CN 201910094028A CN 109766871 A CN109766871 A CN 109766871A
- Authority
- CN
- China
- Prior art keywords
- lai
- crop
- assimilation
- time point
- matrix
- 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.)
- Withdrawn
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法,包括:S1:提取各网格LAI数据时序特征和作物关键物候期数据;S2:获得各网格的作物生长参数和生物物理过程参数;S3:进行逐日作物生长模拟;S4:分别构建模拟LAI矩阵及遥感LAI矩阵;S5:对上述两个LAI矩阵进行归一化处理以提取空间差异特征,然后对空间差异进行同化运算以获得分析LAI矩阵;S6:在各网格对作物模型中的LAI模拟参数进行分别修正,以使模拟LAI空间差异与目标空间差异相符,完成修正后继续运行作物模型至下一同化时间点;S7:重复步骤S4~S6直到设定的同化终点;S8:将作物模型运行至成熟期,输出产量结果。
Description
技术领域
本发明涉及农业作物产量估算技术领域,更具体涉及一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法。
背景技术
作物模型在进行区域尺度的作物产量估算时常常面临输入数据不足、精度不够等问题。地表特征、近地表环境以及作物管理措施往往存在明显的空间差异,而一些必要的模型输入数据,例如作物物候数据,气象数据,作物管理信息等均在站点尺度记录,这导致作物模型在应用到区域尺度时难以获取足够的数据以代表初始条件、作物参数、生长过程等关键因素的空间异质性。卫星遥感数据提供了大范围地表信息的连续监测数据,可以反映地表信息的空间连续性以及时序变化特征。这一优势有效地补充了作物模型在区域尺度上应用中的弱点,因此作物模型与遥感数据同化技术成为了当前改进作物模型区域模拟精度的重要途径,也发展出了多种数据同化方法及其同化流程。
现有的作物模型与遥感数据同化技术在实际应用上仍面临很大的局限性。这些同化技术往往要求高精度的遥感反演数据,例如准确的地表叶面积指数(LAI)反演数据。然而,现有成熟的遥感反演产品,例如MODIS LAI,实际中对LAI的反演存在明显误差,直接同化这些LAI产品会导致模拟产量严重偏离实际情况。常用的解决办法包括利用地表观测数据修正低分辨率遥感LAI数据产品,或使用地表观测数据协助高分辨遥感数据反演发展LAI数据产品。这两种解决办法存在一个共同的不足之处就是要求高空间密度及时间连续性的地表观测数据,人力成本及时间成本非常高,这导致此类解决办法难以在大的空间范围及长的时间尺度上进行推广。
也即,现有的作物模型-遥感数据技术的主要缺点是要求大量的地表观测数据,数据收集的人力及时间成本高,而且点位数据在大面积范围的代表性存在较大的不确定性。使得现有作物模型-遥感数据同化技术难以在大的空间范围及长的时间尺度上进行推广应用。
因此,需要新的技术以至少部分解决现有技术中存在的局限。
发明内容
本发明的发明人在实践以及研究中发现,通过对遥感产品所包含的信息进行提炼,剔除误差较大的信息,保留关键的、可以相对准确地反映地表特征的信息,并使用这些信息与作物模型进行同化。该同化技术将有效避免遥感数据误差导致的模拟结果不切实际地偏离实际观测的问题。同时,该方法可以基于成熟的遥感反演产品,不要求使用地表观测数据对遥感数据进行二次修正或单独反演,可以克服数据收集成本带来的同化技术应用能力不足的问题,并且模拟结果与实际产量之间能够实现良好的吻合。
根据本发明的一方面,提供一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法,包括如下步骤:
S1:获取遥感数据,基于遥感数据提取各网格LAI数据时序特征,按照作物生长期特征识别作物种植区并提取作物关键物候期数据,包括返青期、抽穗期和成熟期数据;
S2:基于获得的物候期数据以及种植区的样点产量记录数据分别对作物模型的作物生长模块参数以及生物物理过程模块参数进行校准,获得各网格各异的作物生长参数以及各网格相同的生物物理过程参数;
S3:基于校准的参数,运行作物模型进行逐日作物生长模拟;
S4:在具备遥感观测的时间点(称为同化时间点)提取种植区内作物模型模拟的各网格LAI值及其对应的遥感反演LAI 值,分别构建模拟LAI矩阵(也称为预报LAI矩阵)及遥感LAI 矩阵(也称为观测LAI矩阵)。
S5:对上述两个LAI矩阵进行归一化处理以提取空间差异特征,然后对空间差异进行同化运算以获得目标空间差异矩阵,也即分析LAI矩阵;
S6:在各网格对作物模型中的LAI模拟参数进行分别修正,以使模拟LAI空间差异与目标空间差异相符,完成修正后继续运行作物模型至下一同化时间点;以及
S7:重复步骤S4~S6直到设定的同化终点;
S8:将作物模型运行至成熟期,输出产量结果。
根据本发明的一个实施方案,步骤S1中,所述遥感数据选自网格分辨率不大于1km×1km的遥感产品。应该理解的是,本发明的原理也可适用于那些高网格分辨率的遥感产品。
根据本发明的一个实施方案,所述遥感产品为GLASS LAI。
根据本发明的一个实施方案,步骤S2中,所述作物模型为 MCWLA作物模型,所述模型参数还包括叶面积生长模块参数。
根据本发明的一个实施方案,步骤S5中,所述归一化处理包括将LAI矩阵除以其矩阵元素的最大值,使原矩阵归一化为取值范围为(0,1]的归一化矩阵,由下式(1)和(2)表示:
其中,表示同化时间点k上的模拟LAI矩阵,Dk代表同一时间点k上对应的遥感LAI矩阵;和分别表示由和Dk转化而来的归一化模拟LAI矩阵和归一化遥感LAI矩阵; Sf和SD是用来对原矩阵进行归一化的归一化函数;
所述同化运算为固定增益卡尔曼滤波算法,由下式(3)表示:
其中是同化时间点k的归一化分析LAI矩阵,H是先验权重系数;通过归一化函数Sf的反函数进行反归一化转化以得到分析LAI矩阵由下式(4)表示:
根据本发明的一个实施方案,步骤S6中,所述LAI模拟参数包括返青期最大叶面积指数LAImax,r1,抽穗期最大叶面积指数 LAImax以及“生殖生长期叶面积衰老系数LAIdg。
根据本发明的一个实施方案,步骤S6中,对LAI模拟参数修正的修正包括:
将MCWLA模型中对LAI的模拟简化为下式(5):
式中是LAI在同化时间点k网格i的LAI模拟值,M代表 MCWLA模型中用来模拟LAI的方程组,代表在同化时间点k-1 到同化时间点k期间模型采用的LAI生长参数(叶面积生长参数),称为预报参数;代表同化时间点k-1网格i的分析 LAI;
在同化时间点k,当计算得到分析LAI,也即后,将基于进行反向优化,由下式(6)和(7)表示:
其中,w是同化时间点k上的模型模拟LAI值与分析LAI值之间的误差绝对值;ΔP代表对的修正量,通过调整ΔP使得 w最小化;如果模拟LAI值与分析LAI值的误差在0.05以下,则判断为达到优化目的,此时优化的调整量被标记为ΔPo;经过ΔPo修正之后得到同化时间点k网格i的分析参数
将该分析参数应用于同化时间点k-1到同化时间点k期间对网格i的LAI模拟;此时,在同化后的模型中,式(5)转化为下式(8):
同时,也作为下一个同化时段内,也即时间点k到时间点k+1期间的预报参数以模拟时间点k+1时的LAI模拟值
根据本发明的一个实施方案,步骤S8包括将网格尺度的模拟产量按照行政边界汇总,输出行政单元尺度的产量模拟结果。
根据本发明的一个实施方案,所述作物为小麦。
本发明避免了同化技术对于地表观测数据的需求,仅使用中低分辨率遥感数据与作物模型同化,即可在大面积、长时间开展较高精度的作物估产。这一新的作物模型-遥感数据同化技术不仅能够提高大面积作物估产的精度,而且极大地降低了应用成本。
附图说明
附图中相同的附图标记标示了相同或类似的部件或部分。本发明的目标及特征考虑到如下结合附图的描述将更加明显,附图中:
图1是根据本发明一个实施方案的基于空间差异同化遥感数据与作物模型的区域农作物估产方法中的同化流程示意图。
图2是根据本发明一个实施方案的实施本发明方法的研究区的示意图;
图3为图2所示研究区中2001-2008年数据同化前后平均模拟产量与实际产量空间分布对比图。
具体实施方式
为清楚的说明本发明中的方案,下面给出优选的实施例并结合附图详细说明。以下的说明本质上仅仅是示例性的而并不是为了限制本公开的应用或用途。
应该理解的是,本发明所引用的作物模型和遥感数据本身是已知的,例如模型的各个子模块、各种参数、运行机制等等,因此本发明重点阐述基于空间差异化的作物模型和遥感数据之间同化过程。
本发明提出一种新的基于LAI时空格局特征的同化方案,具体实施步骤如下:
首先,基于遥感数据提取各网格LAI数据时序特征,按照作物生长期特征识别作物种植区并提取作物关键物候期数据,包括返青期、抽穗期、和成熟期数据。所述遥感数据可以来源于成熟的遥感产品,特别是那些中低网格分辨率的产品,例如网格分辨率不大于1km×1km的那些。例如网格分辨率为1km× 1km的GLASS LAI遥感产品。当然,本发明亦可应用于那些高网格分辨率的产品。
其次,基于物候期数据对作物模型的作物生长模块参数进行校准,获得各网格各异的作物生长参数,并且基于种植区的样点产量记录数据对作物模型生物物理过程模块参数进行校准,获得各网格相同的生物物理过程参数。例如,所述作物模型可以为MCWLA作物模型,应该理解是,本发明的原理也可应用于其他适当的作物模型。模型模拟在网格尺度运行,即对各网格分别模拟作物生长。模型参数可以分为三类,包括作物生长参数,用于控制物候模拟;叶面积生长参数,用于控制叶面积生长模拟;生物物理过程参数,用于控制水平衡、光合作用、呼吸作用、干物质累积等过程。这些具体参数本身为本领域所熟知,例如生物物理过程参数可以包括与土壤水分平衡、光合、呼吸作用以及产量形成等相关的参数。
接着,基于校准后的参数,运行该作物模型进行逐日作物生长模拟。所述模型运行的初始参数为基于研究区内(也即上述识别的作物种植区)典型种植区校准的参数,此套参数作为初始参数应用于所有研究区内网格的模拟初始化。
图1是根据本发明一个实施方案的基于空间差异同化遥感数据与作物模型的区域农作物估产方法中的同化流程示意图。下面结合附图对本发明的同化步骤进行详细的说明,图中以 MCWLA作物模型为例。
如图所示,在具备遥感观测的时间点(称为同化时间点) 提取研究区内作物模型模拟的各网格LAI值及其对应的遥感反演LAI值,分别构建模拟LAI矩阵(也可称为预报LAI矩阵,也即模拟所得的各LAI值的集合)及遥感LAI矩阵(也可称为观测LAI 矩阵,也即遥感反演所得的各LAI值的集合)。图中示出了k+1 个同化时间点,其中k可以为大1的自然数,可以根据具体情况在作物的生长期内进行适当的选择,图中示出了在播种期至成熟期之间。
对两个LAI矩阵进行归一化处理以提取空间差异特征,对空间差异进行同化运算以获得目标空间差异矩阵(称为分析LAI 矩阵)。
所述归一化方法为将LAI矩阵除以其矩阵元素的最大值,使原矩阵归一化为取值范围为(0,1]的归一化矩阵:
其中,表示时间节点k上的LAI预报矩阵,Dk代表由同一时间点上对应的观测矩阵。和分别表示由和 Dk转化而来的归一化的预报LAI矩阵和观测LAI矩阵。Sf和SD是用来对原始矩阵进行归一化的归一化函数。
其中,S5所述同化方法为固定增益卡尔曼滤波算法:
其中是同化时间点k的归一化分析矩阵,和分别表示同时间的归一化模拟LAI矩阵与归一化遥感反演LAI矩阵。H 是一个固定的先验权重系数,取值范围可以为(0,1)。通过归一化函数Sf的反函数(记做)进行反归一化转化以得到分析LAI矩阵如下式(4)所示。
之后,在各网格对模型中的LAI生长参数进行分别修正以使模拟LAI空间差异与目标空间差异相符。
更具体地,所述LAI生长参数包括返青期最大叶面积指数 LAImax,r1,抽穗期最大叶面积指数LAImax以及“生殖生长期叶面积衰老系数LAIdg等。
所述LAI参数的修正并非更改模型输入参数的初始状态,而是对当前时间点上的参数值进行修正,该参数值仅应用于当前同化时间点到下一同化时间点之间的模型模拟。对参数的修正基于以下步骤:
将MCWLA模型中对LAI的模拟简化为如下式(5):
式中是LAI在同化时间点k网格i的LAI模拟值,M代表 MCWLA模型中用来模拟LAI的方程组,代表在同化时间点k-1 到同化时间点k期间模型采用的LAI生长参数,称为预报参数。代表同化时间点k-1网格i的分析LAI。在同化时间点k,当计算得到分析LAI之后,将基于进行反向优化,如下式(6)和(7):
其中,w是同化时间点k模型模拟LAI值与分析LAI值之间的误差绝对值。ΔP代表对的修正量。对的反向修正的目的是通过调整ΔP使得w最小化。如果模拟LAI值与分析LAI值的误差在0.05以下,则认为达到优化目的,满足精度的要求,此时优化的调整量被标记为ΔPo。经过ΔPo修正之后可以得到被称作同化时间点k网格i的分析参数。该参数将应用于同化时间点k-1到同化时间点k期间对网格i的LAI模拟。此时,在同化后的模型中,式(5)可以被转化为下式(8):
同时,也作为下一个同化时段内(时间点k到时间点k+1 期间)的预报参数模拟时间点k+1时的LAI值
也即,在完成修正后继续运行模型至下一同化时间点。如此循环运行,直到设定的同化终点。
最后,将模型运行至成熟期,输出模拟产量结果。例如可以将网格尺度的模拟产量按照行政边界汇总,输出行政单元尺度的产量模拟结果。
本发明的方法可以适用于多种作物例如小麦、大豆等等。
实施例
本案例以估测华北平原地区冬小麦产量为例进一步阐述本发明的技术方案。包括以下步骤:
选择华北平原中部作为研究区域,研究时间为2001-2008 年。基于1km×1kmGLASS LAI遥感数据,提取各网格LAI数据时序特征,按照作物生长期特征识别作物种植区并提取作物关键物候期数据,研究区如图2所示。
基于物候期数据对MCWLA作物模型的作物生长模块参数进行校准,获得各网格各异的作物生长参数。获取样点县的记录产量数据,并依据产量数据对生物物理过程参数进行校准。校准所得生物物理过程参数应用到整个研究区的模拟中。对于叶面积生长参数(“返青期最大叶面积指数LAImax,r1”,“抽穗期最大叶面积指数LAImax”以及“生殖生长期叶面积衰老系数LAIdg”) 采用其模型默认值。
基于本发明提出的技术方案的叶面积同化实施过程如下:
从播种时期开始使用MCWLA模型进行模拟,在具备遥感观测的时间点(称为同化时间点)提取研究区内作物模型模拟的各网格LAI值及其对应的遥感反演LAI值,分别构建模拟LAI矩阵 (也可称为预报LAI矩阵)及遥感LAI矩阵(也可称为观测LAI 矩阵)。
对两个LAI矩阵进行归一化处理以提取空间差异特征,对空间差异进行同化运算以获得目标空间差异矩阵(称为分析LAI 矩阵)。
所述归一化方法为将LAI矩阵除以其矩阵元素的最大值,使原矩阵归一化为取值范围为(0,1]的归一化矩阵:
其中,表示时间节点k上的LAI预报矩阵,Dk代表由同一时间点上对应的观测矩阵。和分别表示由和 Dk转化而来的归一化的预报LAI矩阵和观测LAI矩阵。Sf和SD是用来对原始矩阵进行归一化的归一化函数。
其中,S5所述同化方法为固定增益卡尔曼滤波算法:
其中是同化时间点k的归一化分析矩阵,和分别表示同时间的归一化模拟LAI矩阵与归一化遥感反演LAI矩阵。H 是一个固定的先验权重系数,本案例中取值0.53。通过归一化函数Sf的反函数(记做)进行反归一化转化以得到分析 LAI矩阵如下式(4)所示。
之后,在各网格对上述提及的模型中的LAI生长参数进行分别修正以使模拟LAI空间差异与目标空间差异相符。对参数的修正基于以下步骤:
将MCWLA模型中对LAI的模拟简化为如下式(5):
式中是LAI在同化时间点k网格i的LAI模拟值,M代表 MCWLA模型中用来模拟LAI的方程组,代表在同化时间点k-1 到同化时间点k期间模型采用的LAI生长参数,称为预报参数。代表同化时间点k-1网格i的分析LAI。在同化时间点k,当计算得到分析LAI之后,将基于进行反向优化,如下式(6)和(7):
其中,w是同化时间点k模型模拟LAI值与分析LAI值之间的误差绝对值。ΔP代表对的修正量。对的反向修正的目的是通过调整ΔP使得w最小化。如果模拟LAI值与分析LAI值的误差在0.05以下,则认为达到优化目的,满足精度的要求,此时优化的调整量被标记为ΔPo。经过ΔPo修正之后可以得到被称作同化时间点k网格i的分析参数。该参数将应用于同化时间点k-1到同化时间点k期间对网格i的LAI模拟。此时,在同化后的模型中,式(5)可以被转化为下式(8):
同时,也作为下一个同化时段内(时间点k到时间点k+1 期间)的预报参数模拟时间点k+1时的LAI值
也即,在完成修正后继续运行模型至下一同化时间点。如此循环运行,直到设定的同化终点。在本案例中,同化步长同 GLASS LAI数据步长,为8天。同化起点设置为返青期,终点设置为灌浆期。同化运行至灌浆期后结束,模型继续运行至成熟期,输出作物产量并聚合至县级尺度与产量记录数据进行比较。
本发明实施例所述的估产方法,融合了遥感数据和作物模型的优势,将LAI的空间差异作为同化变量,避免了遥感数据中LAI偏低的不利影响,实现了对区域产量估算精度的提升。与未同化相比,2001-2008年均方根误差(RMSE)平均减小了33.1%, 最大减小50.2%。R2平均提升0.13,最高提升0.27。估产精度明显提升,区域空间分布趋势也更为符合现实。图3为2001-2008 年研究区内同化前后平均模拟产量与实际产量空间分布对比图。
本发明融合了遥感数据和作物模型的优势,利用成熟的 LAI遥感反演产品以及MCWLA作物模型进行同化估产。使用空间差异同化的方法对模型模拟的LAI进行了调整,实现了遥感数据与作物模型的同化。同化后的模拟产量与未同化相比较,均方根误差(RMSE)减小而决定系数(R2)明显上升,作物模型产量估算的精度有显著提高,产量空间分布趋势与统计产量一致。同时,本发明不要求地表观测数据,不要求对遥感产品进行二次修正,显著降低了同化技术的成本,提升了同化技术的实际应用能力。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的装置及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法,包括如下步骤:
S1:获取遥感数据,基于遥感数据提取各网格LAI数据时序特征,按照作物生长期特征识别作物种植区并提取作物关键物候期数据,包括返青期、抽穗期、和成熟期数据;
S2:基于获得的物候期数据以及种植区的样点产量记录数据分别对作物模型的作物生长模块参数以及生物物理过程模块参数进行校准,获得各网格各异的作物生长参数以及各网格相同的生物物理过程参数;
S3:基于校准后的参数,运行作物模型进行逐日作物生长模拟;
S4:在具备遥感观测的时间点提取种植区内作物模型模拟的各网格LAI值及其对应的遥感反演LAI值,分别构建模拟LAI矩阵及遥感LAI矩阵;
S5:对上述两个LAI矩阵进行归一化处理以提取空间差异特征,然后对空间差异进行同化运算以获得目标空间差异矩阵,也即分析LAI矩阵;
S6:在各网格对作物模型中的LAI模拟参数进行分别修正,以使模拟LAI空间差异与目标空间差异相符,完成修正后继续运行作物模型至下一同化时间点;以及
S7:重复步骤S4~S6直到设定的同化终点;
S8:将作物模型运行至成熟期,输出产量结果。
2.根据权利要求1所述的区域农作物估产方法,其特征在于,步骤S1中,所述遥感数据选自网格分辨率不大于1km×1km的遥感产品。
3.根据权利要求2所述的区域农作物估产方法,其特征在于,所述遥感产品为GLASSLAI。
4.根据权利要求1所述的区域农作物估产方法,其特征在于,步骤S2中,所述作物模型为MCWLA作物模型,所述模型参数还包括叶面积生长模块参数。
5.根据权利要求1所述的区域农作物估产方法,其特征在于,步骤S5中,所述归一化处理包括将LAI矩阵除以其矩阵元素的最大值,使原矩阵归一化为取值范围为(0,1]的归一化矩阵,由下式(1)和(2)表示:
其中,表示同化时间点k上的模拟LAI矩阵,Dk代表同一时间点k上对应的遥感LAI矩阵;和分别表示由和Dk转化而来的归一化模拟LAI矩阵和归一化遥感LAI矩阵;Sf和SD是用来对原矩阵进行归一化的归一化函数;
所述同化运算为固定增益卡尔曼滤波算法,由下式(3)表示:
其中是同化时间点k的归一化分析LAI矩阵,H是先验权重系数;通过归一化函数Sf的反函数进行反归一化转化以得到分析LAI矩阵由下式(4)表示:
6.根据权利要求1所述的区域农作物估产方法,其特征在于,步骤S6中,所述LAI模拟参数包括返青期最大叶面积指数LAImax,r1,抽穗期最大叶面积指数LAImax以及生殖生长期叶面积衰老系数LAIdg。
7.根据权利要求4所述的区域农作物估产方法,其特征在于,步骤S6中,对LAI模拟参数的修正包括:
将MCWLA模型中对LAI的模拟简化为如下式(5):
式中是LAI在同化时间点k网格i的LAI模拟值,M代表MCWLA模型中用来模拟LAI的方程组,代表在同化时间点k-1到同化时间点k期间模型采用的LAI生长参数,称为预报参数;代表同化时间点k-1网格i的分析LAI;
在同化时间点k,当计算得到分析LAI,也即后,将基于进行反向优化,由下式(6)和(7)表示:
其中,w是同化时间点k上的模型模拟LAI值与分析LAI值之间的误差绝对值;ΔP代表对的修正量,通过调整ΔP使得w最小化;如果模拟LAI值与分析LAI值的误差在0.05以下,则判断为达到优化目的,此时优化的调整量被标记为ΔPo;经过ΔPo修正之后得到同化时间点k网格i的分析参数
将该分析参数应用于同化时间点k-1到同化时间点k期间对网格i的LAI模拟;此时,在同化后的模型中,式(5)转化为下式(8):
同时,也作为下一个同化时段内,也即时间点k到时间点k+1期间的预报参数以模拟时间点k+1时的LAI模拟值
8.根据权利要求1所述的区域农作物估产方法,其特征在于,步骤S8包括将网格尺度的模拟产量按照行政边界汇总,输出行政单元尺度的产量模拟结果。
9.根据权利要求1所述的区域农作物估产方法,其特征在于,所述作物为小麦。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910094028.9A CN109766871A (zh) | 2019-01-30 | 2019-01-30 | 一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910094028.9A CN109766871A (zh) | 2019-01-30 | 2019-01-30 | 一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109766871A true CN109766871A (zh) | 2019-05-17 |
Family
ID=66455732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910094028.9A Withdrawn CN109766871A (zh) | 2019-01-30 | 2019-01-30 | 一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109766871A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766308A (zh) * | 2019-10-17 | 2020-02-07 | 中国科学院地理科学与资源研究所 | 一种基于集合同化策略的区域农作物估产方法 |
CN112446155A (zh) * | 2020-12-09 | 2021-03-05 | 四川省农业科学院农业信息与农村经济研究所 | 一种目标农作物的空间格局模拟模型获取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050234691A1 (en) * | 2004-04-20 | 2005-10-20 | Singh Ramesh P | Crop yield prediction |
CN106845428A (zh) * | 2017-01-26 | 2017-06-13 | 中国科学院遥感与数字地球研究所 | 一种作物产量遥感估算方法及系统 |
US20170228743A1 (en) * | 2016-02-05 | 2017-08-10 | Weather Analytics, LLC | Crop forecasting with incremental feature selection and spectrum constrained scenario generation |
-
2019
- 2019-01-30 CN CN201910094028.9A patent/CN109766871A/zh not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050234691A1 (en) * | 2004-04-20 | 2005-10-20 | Singh Ramesh P | Crop yield prediction |
US20170228743A1 (en) * | 2016-02-05 | 2017-08-10 | Weather Analytics, LLC | Crop forecasting with incremental feature selection and spectrum constrained scenario generation |
CN106845428A (zh) * | 2017-01-26 | 2017-06-13 | 中国科学院遥感与数字地球研究所 | 一种作物产量遥感估算方法及系统 |
Non-Patent Citations (1)
Title |
---|
YI CHEN ETC.: ""Improving regional winter wheat yield estimation through assimilation of phenology and leaf area index from remote sensing data"", 《EUROPEAN JOURNAL OF AGRONOMY》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766308A (zh) * | 2019-10-17 | 2020-02-07 | 中国科学院地理科学与资源研究所 | 一种基于集合同化策略的区域农作物估产方法 |
CN112446155A (zh) * | 2020-12-09 | 2021-03-05 | 四川省农业科学院农业信息与农村经济研究所 | 一种目标农作物的空间格局模拟模型获取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112163375B (zh) | 一种基于神经网络的长时间序列近地面臭氧反演方法 | |
CN106485002B (zh) | 在复杂地形气候区域估算甘蔗潜在产量的方法 | |
CN102162850A (zh) | 基于模型的小麦产量遥感监测预报方法 | |
CN102651096A (zh) | 同化叶面积指数时序曲线特征的冬小麦估产方法 | |
CN112598277B (zh) | 一种评估跨区域缩小冬小麦产量差与提高氮肥效率的方法 | |
CN105740759A (zh) | 基于多时相数据中特征提取的中稻信息决策树分类方法 | |
CN106483147B (zh) | 基于多源数据的长时间序列被动微波土壤水分精度改进研究方法 | |
CN107941713A (zh) | 一种基于耦合作物模型同化光谱反射率的水稻估产方法 | |
CN110705182B (zh) | 耦合作物模型和机器学习的作物育种适应时间预测方法 | |
CN109766871A (zh) | 一种基于空间差异同化遥感数据与作物模型的区域农作物估产方法 | |
CN108537679A (zh) | 遥感与作物模型融合的区域尺度作物出苗日期估算方法 | |
CN115759524A (zh) | 一种基于遥感影像植被指数的土壤生产力等级识别方法 | |
CN110008621B (zh) | 基于双重流依赖集合平方根滤波同化算法的作物模型遥感同化估产方法 | |
CN104933699A (zh) | 基于高斯函数拟合方差自动提取地表植被物候信息的方法 | |
CN114967798A (zh) | 一种基于互联网的天麻种植管理控制系统 | |
CN109614763B (zh) | 一种基于多源信息分步校正作物模型的区域作物估产方法 | |
Cui et al. | Estimating maize yield in the black soil region of Northeast China using land surface data assimilation: integrating a crop model and remote sensing | |
Cheng et al. | Improving soil available nutrient estimation by integrating modified WOFOST model and time-series earth observations | |
CN116595333A (zh) | 土壤-气候智慧型的水稻目标产量与氮肥用量确定方法 | |
Ahuja et al. | A synthesis of current parameterization approaches and needs for further improvements | |
CN116579872A (zh) | 一种基于作物生长模型和天气预报的精准灌溉决策方法 | |
CN115392016A (zh) | 一种基于遥感数据同化的青贮玉米生长发育预测方法 | |
CN109272416A (zh) | 一种大棚农作物种植系统 | |
CN109359862A (zh) | 一种粮食作物实时估产方法及系统 | |
CN110766308B (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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20190517 |
|
WW01 | Invention patent application withdrawn after publication |