CN113592201B - 一种拔节后受冻冬小麦产量损失早期预测方法 - Google Patents
一种拔节后受冻冬小麦产量损失早期预测方法 Download PDFInfo
- Publication number
- CN113592201B CN113592201B CN202111013108.0A CN202111013108A CN113592201B CN 113592201 B CN113592201 B CN 113592201B CN 202111013108 A CN202111013108 A CN 202111013108A CN 113592201 B CN113592201 B CN 113592201B
- Authority
- CN
- China
- Prior art keywords
- winter wheat
- yield
- jointing
- frozen
- index
- 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.)
- Active
Links
- 241000209140 Triticum Species 0.000 title claims abstract description 131
- 235000021307 Triticum Nutrition 0.000 title claims abstract description 131
- 238000000034 method Methods 0.000 title claims abstract description 89
- 230000006378 damage Effects 0.000 claims abstract description 53
- 238000007710 freezing Methods 0.000 claims abstract description 40
- 230000008014 freezing Effects 0.000 claims abstract description 40
- 208000027418 Wounds and injury Diseases 0.000 claims abstract description 35
- 208000014674 injury Diseases 0.000 claims abstract description 35
- 230000009467 reduction Effects 0.000 claims abstract description 30
- 244000025254 Cannabis sativa Species 0.000 claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 21
- 238000011161 development Methods 0.000 claims abstract description 14
- 238000005070 sampling Methods 0.000 claims description 23
- 235000013339 cereals Nutrition 0.000 claims description 21
- 230000018109 developmental process Effects 0.000 claims description 12
- 238000001816 cooling Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 9
- 239000003814 drug Substances 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000004069 differentiation Effects 0.000 claims description 5
- 238000011835 investigation Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 5
- 210000005069 ears Anatomy 0.000 claims description 4
- 230000017260 vegetative to reproductive phase transition of meristem Effects 0.000 claims description 4
- 241000196324 Embryophyta Species 0.000 claims description 3
- 238000013135 deep learning Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000003306 harvesting Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 2
- 230000001174 ascending effect Effects 0.000 abstract 1
- 230000009746 freeze damage Effects 0.000 description 18
- 238000010586 diagram Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 208000024891 symptom Diseases 0.000 description 3
- 238000005259 measurement Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000008571 general function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000000246 remedial effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Marine Sciences & Fisheries (AREA)
- Primary Health Care (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- General Health & Medical Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种拔节后受冻冬小麦产量损失早期预测方法，主要步骤为：确定冬小麦拔节开始日期dj；在所述冬小麦拔节开始日期开始观测麦田逐小时草面温度Tt；根据冻害对冬小麦不同发育阶段的影响构建脆弱性系数VC曲线，所述脆弱性系数VC随冬小麦发育进程整体上呈上升趋势；根据以往观测数据，预选冻害温度阈值Th；构建综合霜冻胁迫指数IFS；通过设定不同所述温度阈值Th以计算所述IFS指数，利用所述IFS指数和受冻冬小麦的减产率进行拟合，依据二者间模型的决定系数和均方根误差来确定合适的冻害温度阈值Th。该方法不仅适用于拔节后受冻冬小麦产量损失的早期预测目的，而且预测精度显著提升，适用性强、推广应用价值高。
Description
技术领域
本发明属于小麦产量预测技术领域，特别涉及一种拔节后受冻冬小麦产量损失早期预测方法。
背景技术
黄淮冬麦区是中国小麦重要产区，大部分属黄淮平原，地势平坦，春季寒流入侵频繁，可对拔节后冬小麦幼穗造成显著伤害，进而导致产量损失。大田小麦遭受晚霜冻害的症状通常在冻后5至10天内显现，随冻害程度加重，叶片枯黄、幼穗萎蔫趋势愈发明显。鉴于冬小麦在冻后仍须较长时间(至少1个多月)才成熟，因而为实施各种补救措施或管理战略(如粮食宏观调控)等提供了契机，这就要求在短期内对冬小麦减产风险及空间分布做出早期预测，这是确保各类措施精准和高效实施的前提。因此，早期预测应将关注点放在冻害症状显现之前，这对冻后管理具有重要指导价值。
当前，可用于对拔节后受冻冬小麦产量损失进行早期预测的技术手段主要是田间调查、气象观测和遥感手段。田间调查手段通过考察幼穗冻害来评估产量损失，精度较高，但是仅限于田块内的产量损失预测，针对大范围空间区域的产量损失早期预测而言，该手段耗时且工作量巨大。遥感手段主要依靠逐日时相的地表温度(LST)或植被指数(如NDVI)数据。由于冻害通常发生于0至7点之间，与卫星在当地过境时间并不一致，且LST的遥感反演精度受下垫面影响而不确定性较大，因而基于LST数据的方法适用性较差。常见遥感方法是依据冻后植被指数随时间变化规律以构建冻害遥感指标，进而评估冬小麦实际冻害和预测产量损失。此类方法未从冬小麦冠层光谱响应冻害机制的角度来挖掘敏感的遥感指数，而是直接将植被指数作为可响应冻害的遥感指数，由于植被指数响应冻害具有时间滞后性，在冻害症状显现前的敏感性相对较差，会影响其预测精度。与遥感手段相比，气象手段评估冻害的时效性更强、数据精度和稳定性更高，因而在晚霜冻害评估中应用广泛。气象手段主要依靠气象观测台站记录的百叶箱温度数据，通过集成低温强度和其持续时长，来构建霜冻害胁迫指标。实际上，除了低温强度和持续时长外，冬小麦发育进程是另一个影响冻害的关键因子，但现有指标未能将发育进程因子集成进来。而且，由于受下垫面环境等的影响，百叶箱温度在判定冻害时的稳定性相对较差。可见，当前方法由于指标建设方面的不足，存在着预测精度较低的问题。为此，本发明立足于拔节后冻害影响机制、以及产量损失机制等的研究基础，考虑加入发育进程和草面温度，全新构建对冻害最敏感的指数，进而解决产量损失早期预测精度较低的问题。
发明内容
本发明提供一种受冻冬小麦产量损失预测方法，用于解决受冻冬小麦产量损失早期预测精度较低的技术问题。
本发明的实施例内容如下：
一种拔节后受冻冬小麦产量损失早期预测方法，包括以下步骤：
确定冬小麦拔节开始日期dj；在冬小麦拔节开始日期开始观测麦田逐小时草面温度Tt；根据冻害对冬小麦不同发育阶段的影响构建脆弱性系数VC曲线，脆弱性系数VC随冬小麦发育进程整体上呈上升趋势；根据以往观测数据，预选冻害温度阈值Th；
构建综合霜冻胁迫指数IFS：
其中，FDt为逐小时霜冻度日(单位：h℃)，FDDi为第i个降温过程中逐小时霜冻度日FDt的全天平均值(单位：d℃)，Tt为第i个降温过程中第t小时的最低草面温度(℃)，Th为冻害温度阈值(℃)，VCi为第i个降温过程中冬小麦的脆弱性系数(无量纲，介于0.1─1.0之间)，dj和da分别表示冬小麦拔节开始日期和开花期；
通过设定不同温度阈值Th以计算IFS指数，利用IFS指数和受冻冬小麦的减产率进行拟合，依据二者间模型的决定系数和均方根误差来确定合适的冻害温度阈值Th；基于IFS指数预测的减产率与实测减产率对比进行精度验证。
进一步地，冬小麦拔节开始日期dj由如下方法之一确定：
田间抽样调查，定期观察并记录植株生长情况，一旦发现节间露出地面1.5―2.0cm，则认为拔节开始；在较大范围空间区域内，如果超过半数的抽样点小麦均进入拔节，则认为该地区冬小麦拔节开始，对应日期为拔节开始日期；
基于积温的推算，将冬小麦返青至拔节所需有效积温的多年观测平均或经验值，作为判断冬小麦是否进入拔节的依据，利用气象观测台站记录的逐日平均气温，计算返青至拔节的有效积温，如果该有效积温达到或超过经验值，则认为拔节开始，对应日期为拔节开始日期；
基于植被指数检测，在冬小麦整个生育期内植被指数时序曲线上，将植被指数在春季由低值平缓区快速升高时对应的时序特征点，作为判断冬小麦是否进入拔节的依据，利用时序曲线函数模拟方法计算每个时序点上的曲率变化率，如果该点曲率变化率达到极大值，则认为拔节开始，对应日期为所述拔节开始日期。。
进一步地，麦田逐小时草面温度Tt为气象站点观测到的逐小时草面温度数据，或将气象站点观测到的逐小时草面温度数据进行空间插值而生成逐小时草面温度空间分布数据。空间插值方法可采用样条函数方法。
进一步地，脆弱性系数VC曲线呈4段式上升；第1时段从小花原基分化开始直至药隔形成期前1日；第2时段为药隔形成前1日内；第3时段从药隔形成期至四分体形成期；第4时段在四分体形成期后。
进一步地，脆弱性系数VC曲线呈线性上升，其中，第1时段中所述脆弱性系数从0.1升至0.4；第2时段中所述脆弱性系数从0.4升至0.7；第3时段中所述脆弱性系数从0.7升至1.0；第4时段中所述脆弱性系数为1.0。
进一步地，合适的冻害温度阈值Th确定方法为如下步骤：
测定所述的受冻冬小麦产量减产率；
设定所述冻害温度阈值Th的若干候选值为-6℃、-5℃、-4℃、-3℃、-2℃、-1℃和0℃，分别计算每个采样点上所述候选值对应的所述IFS指数；
将所述IFS指数和所述减产率进行线性拟合，计算不同所述冻害温度阈值Th条件下二者线性模型的决定系数R2和均方根误差RMSE，依据所述R2最大和所述RMSE最小的原则来确定所述冻害温度阈值Th。
进一步地，受冻冬小麦产量减产率测定方法为：
对照产量Yn测定，通过田间抽样调查的办法，在研究区晚霜冻害事件发生年份，逐样点框定1米×1米的采样范围，并将全部茎蘖收割、装袋且带回实验室凉晒干；针对每个样点，统计全部茎蘖数量，以作为正常生长条件下冬小麦的穗数ENn；挑选出全部正常穗，统计其粒数和粒重，进而计算穗粒数和千粒重，以作为正常生长条件下冬小麦的穗粒数GNPEn和千粒重TGWn；根据公式ENn×GNPEn×TGWn，计算得出正常生长条件下冬小麦的产量Yn，作为所述对照产量；实际产量Yf测定，统计全部茎蘖的籽粒重，作为所述实际产量Yf；利用公式(Yn–Yf)/Yn×100％，计算得出所述受冻冬小麦产量减产率。
进一步地，受冻冬小麦产量减产率测定方法也可以为：
基于多年测定的冬小麦产量数据，将实际产量分解为趋势产量、气象产量和随机产量3个部分，随机产量因其影响小而忽略不计，将趋势产量作为正常产量，计算实际产量和趋势产量之差，可以得到受冻冬小麦产量减产率。
进一步地，针对空间大范围区域内的受冻冬小麦，采用基于面向对象分类、深度学习等方法的遥感分类技术，可提取出冬小麦空间分布信息，屏蔽非小麦信息。
本发明提出拔节后受冻冬小麦产量损失早期预测方法，不仅适用于拔节后受冻冬小麦产量损失的早期预测目的，而且预测精度显著提升，仅需输入逐小时草面温度数据即可，指标和模型的区域适用性和推广应用价值较高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1拔节后受冻冬小麦产量损失早期预测技术体系示意图；
图2综合霜冻胁迫指数构建示意图；
图3脆弱性系数随幼穗分化进程的变化；
图4利用不同冻害温度阈值计算的累积霜冻度日AFDD1、AFDD2和IFS指数及其与减产率线性拟合模型的精度；
图5基于IFS指数预测的减产率与实测减产率的散点关系。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
根据本发明的实施例，提供了一种受冻冬小麦产量损失早期预测方法，图1为所述拔节后受冻冬小麦产量损失早期预测技术体系示意图，包括以下步骤：
确定冬小麦拔节开始日期dj；在冬小麦拔节开始日期开始观测麦田逐小时草面温度Tt；根据冻害对冬小麦不同发育阶段的影响构建脆弱性系数VC曲线，脆弱性系数VC随冬小麦发育进程整体上呈上升趋势；根据以往观测数据，预选冻害温度阈值Th；
构建综合霜冻胁迫指数IFS：本发明提出构建基于低温、低温持续时长以及发育进程综合影响的晚霜冻害胁迫指标。该指标的通用表达式如下，
IFS＝f(LT,DU,DP)
其中，IFS表示综合霜冻胁迫(integratedfroststress,IFS)指数，LT表示晚霜冻害事件中低于冻害温度阈值的温度，简称“低温”；DU表示低温持续时长；DP表示冬小麦发育进程；f表示一种通用函数式。
IFS的量化公式为：
其中，FDt为逐小时霜冻度日(单位：h℃)，FDDi为第i个降温过程中逐小时霜冻度日FDt的全天平均值(单位：d℃)，Tt为第i个降温过程中第t小时的最低草面温度(℃)，Th为冻害温度阈值(℃)，VCi为第i个降温过程中冬小麦的脆弱性系数(无量纲，介于0.1─1.0之间)，dj和da分别表示冬小麦拔节开始日期和开花期；
通过设定不同温度阈值Th以计算IFS指数，利用IFS指数和受冻冬小麦的减产率进行拟合，依据二者间模型的决定系数和均方根误差来确定合适的冻害温度阈值Th；
将各参数带入模型，对基于IFS指数预测的减产率与实测减产率对比进行精度验证。
图2是综合霜冻胁迫指数构建示意图。晚霜冻害事件的定义须考虑三个条件：一是冬小麦生育时期，二是晚霜冻害发生时段，三是温度阈值。在此，将晚霜冻害定义为发生在冬小麦拔节至开花期之间(大致处于春季3-5月份)的一种气象灾害。既然每日降温过程总是介于13:00至次日12:00之间，且晚霜冻害事件通常发生在春季夜晚，因此将每日降温过程作为一次晚霜冻害事件的发生时段。在这个时段内，一旦温度降至冻害温度阈值以下(即低温出现)，则一次晚霜冻害事件被认为开始发生；直至该时段内低温不再被检测到为止，则认为本次晚霜冻害事件结束。横坐标表示冬小麦拔节后日序(主刻度表示日，副刻度表示小时)；红色曲线表示逐小时草面温度，其值(单位：℃)如左侧纵坐标所示；水平方向的黄色虚线表示冻害温度阈值；低温(低于冻害温度阈值)的出现，意味着晚霜冻害事件发生，灰色柱状图即代表了一次晚霜冻害事件。
在本发明的一优选实施例中，冬小麦拔节开始日期dj由如下方法之一确定：
田间抽样调查，定期观察并记录植株生长情况，一旦发现节间露出地面1.5―2.0cm，则认为拔节开始；在较大范围空间区域内，如果超过半数的抽样点小麦均进入拔节，则认为该地区冬小麦拔节开始，对应日期为拔节开始日期；
基于积温的推算，将冬小麦返青至拔节所需有效积温的多年观测平均或经验值，作为判断冬小麦是否进入拔节的依据，利用气象观测台站记录的逐日平均气温，计算返青至拔节的有效积温，如果该有效积温达到或超过经验值，则认为拔节开始，对应日期为拔节开始日期；
基于植被指数检测，在冬小麦整个生育期内植被指数时序曲线上，将植被指数在春季由低值平缓区快速升高时对应的时序特征点，作为判断冬小麦是否进入拔节的依据，利用时序曲线函数模拟方法计算每个时序点上的曲率变化率，如果该点曲率变化率达到极大值，则认为拔节开始，对应日期为所述拔节开始日期。
在本发明的一优选实施例中，麦田逐小时草面温度Tt为气象站点观测到的逐小时草面温度数据，或将气象站点观测到的逐小时草面温度数据进行空间插值而生成逐小时草面温度空间分布数据。空间插值方法可采用样条函数方法。
在本发明的一优选实施例中，脆弱性系数VC曲线呈4段式上升；图3是脆弱性系数随幼穗分化进程的变化图，第1时段从小花原基分化(图中原点符号)开始直至药隔形成期(图中方形符号)前1日；第2时段为药隔形成前1日内；第3时段从药隔形成期至四分体形成期(图中三角形符号)；第4时段在四分体形成期后。图中对脆弱性系数VC变化及取值进行进一步优选，曲线呈线性上升，其中，第1时段中所述脆弱性系数从0.1升至0.4；第2时段中所述脆弱性系数从0.4升至0.7；第3时段中所述脆弱性系数从0.7升至1.0；第4时段中所述脆弱性系数为1.0。
在本发明的一优选实施例中，合适的冻害温度阈值Th确定方法为如下步骤：
测定所述的受冻冬小麦产量减产率；
设定所述冻害温度阈值Th的若干候选值为-6℃、-5℃、-4℃、-3℃、-2℃、-1℃和0℃，分别计算每个采样点上所述候选值对应的所述IFS指数；
将所述IFS指数和所述减产率进行线性拟合，计算不同所述冻害温度阈值Th条件下二者线性模型的决定系数R2和均方根误差RMSE，依据所述R2最大和所述RMSE最小的原则来确定所述冻害温度阈值Th。
在本发明的一优选实施例中，受冻冬小麦产量减产率测定方法为：对照产量Yn测定，通过田间抽样调查的办法，在研究区晚霜冻害事件发生年份，逐样点框定1米×1米的采样范围，并将全部茎蘖收割、装袋且带回实验室凉晒干；针对每个样点，统计全部茎蘖数量，以作为正常生长条件下冬小麦的穗数ENn；挑选出全部正常穗，统计其粒数和粒重，进而计算穗粒数和千粒重，以作为正常生长条件下冬小麦的穗粒数GNPEn和千粒重TGWn；根据公式ENn×GNPEn×TGWn，计算得出正常生长条件下冬小麦的产量Yn，作为所述对照产量；实际产量Yf测定，统计全部茎蘖的籽粒重，作为所述实际产量Yf；利用公式(Yn–Yf)/Yn×100％，计算得出所述受冻冬小麦产量减产率。
在本发明的一优选实施例中，受冻冬小麦产量减产率测定方法也可以为：基于多年测定的冬小麦产量数据，将实际产量分解为趋势产量、气象产量和随机产量3个部分，随机产量因其影响小而忽略不计，将趋势产量作为正常产量，计算实际产量和趋势产量之差，可以得到受冻冬小麦产量减产率。
在本发明的一优选实施例中，针对空间大范围区域内的受冻冬小麦，采用基于面向对象分类、深度学习等方法的遥感分类技术，可提取出冬小麦空间分布信息，屏蔽非小麦信息。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
下面结合一具体实例来说明本发明的应用。以河南省商丘地区2015至2019年冬小麦拔节后霜冻害事件为例，各年份均于冬小麦成熟后，在研究区内开展随机采样、产量测定和减产率计算工作，所采样点数量各年分别为42、27、32、43和19个，总计达163个。
经田间定期抽样调查和记录，得知2015至2019年商丘地区冬小麦分别于3月24日、3月18日、3月19日、3月20日和3月20日进入拔节。
获得商丘及其周边地区27个气象观测台站记录的逐小时草面温度数据，对拔节开始后的台站草面温度数据按逐小时进行空间插值，生成时空连续(时间分辨率：逐小时；空间分辨率：1km)的草面温度图像。在该图像上提取所有采样点对应的草面温度值，并将草面温度值输入综合霜冻胁迫(IFS)指数、累积霜冻度日AFDD1指数和AFDD2指数的计算公式中，分别计算-6、-5、-4、-3、-2、-1和0℃冻害温度阈值条件下的IFS、AFDD1和AFDD2指数值。
AFDD1和AFDD2的计算公式如下，
分别构建2015–2018年IFS、AFDD1、AFDD2与减产率之间的线性拟合模型，并计算不同冻害温度阈值下的线性模型决定系数(R2)和均方根误差(RMSE)，如图4所示。根据R2最大和RMSE最小原则，最终确定IFS指数的冻害温度阈值为-3℃。不论冻害温度阈值候选值大小，IFS指数总是比AFDD1、AFDD2有更高的R2和更低的RMSE。可见，IFS指数表现出比AFDD1和AFDD2指数更优的产量损失早期预测精度。AFDD1和AFDD2指数使用的是逐日最低百叶箱气温，缺乏发育进程对冻害影响的考虑，基于百叶箱气温的冻害阈值设定并不准确(在平流型霜冻、辐射型霜冻和混合型霜冻条件下，同一阈值对应的冻害程度明显不同)，预测精度相对较低。
如图5所示为基于IFS指数预测的减产率与实测减产率的散点关系，图中利用2019年减产率数据，对基于-3℃阈值的IFS指数的预测精度进行验证，发现其早期预测减产率的决定系数为0.463，达到0.015显著性水平，拟合线非常接近1:1线。
Claims (8)
1.一种拔节后受冻冬小麦产量损失早期预测方法，其特征在于，包括以下步骤：
确定冬小麦拔节开始日期dj；
在所述冬小麦拔节开始日期开始观测麦田逐小时草面温度Tt；
根据冻害对冬小麦不同发育阶段的影响构建脆弱性系数VC曲线，所述脆弱性系数VC随冬小麦发育进程整体上呈分段式线性上升，第1时段从小花原基分化开始直至药隔形成期前1日，所述脆弱性系数从0.1升至0.4；第2时段为药隔形成前1日内，所述脆弱性系数从0.4升至0.7；第3时段从药隔形成期至四分体形成期，所述脆弱性系数从0.7升至1.0；第4时段在四分体形成期后，所述脆弱性系数为1.0；
根据以往观测数据，预选冻害温度阈值Th；
构建综合霜冻胁迫指数IFS：
其中，FDt为逐小时霜冻度日(单位：h℃)，FDDi为第i个降温过程中逐小时霜冻度日FDt的全天平均值(单位：d℃)，Tt为第i个降温过程中第t小时的最低草面温度(℃)，Th为冻害温度阈值(℃)，VCi为第i个降温过程中冬小麦的脆弱性系数(无量纲，介于0.1-1.0之间)，dj和da分别表示冬小麦拔节开始日期和开花期；
通过设定不同所述温度阈值Th以计算所述IFS指数，利用所述IFS指数和受冻冬小麦的减产率进行拟合，依据二者间模型的决定系数和均方根误差来确定合适的冻害温度阈值Th；
基于IFS指数预测的减产率与实测减产率对比进行精度验证。
2.根据权利要求1所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述冬小麦拔节开始日期dj由如下方法之一确定：
田间抽样调查，定期观察并记录植株生长情况，一旦发现节间露出地面1.5-2.0cm，则认为拔节开始；在某一地区空间区域内，如果超过半数的抽样点小麦均进入拔节，则认为该地区冬小麦拔节开始，对应日期为所述拔节开始日期；
基于积温的推算，将冬小麦返青至拔节所需有效积温的多年观测平均或经验值，作为判断冬小麦是否进入拔节的依据，利用气象观测台站记录的逐日平均气温，计算返青至拔节的有效积温，如果该有效积温达到或超过经验值，则认为拔节开始，对应日期为所述拔节开始日期；
基于植被指数检测，在冬小麦整个生育期内植被指数时序曲线上，将植被指数在春季由低值平缓区快速升高时对应的时序特征点，作为判断冬小麦是否进入拔节的依据，利用时序曲线函数模拟方法计算每个时序点上的曲率变化率，如果该点曲率变化率达到极大值，则认为拔节开始，对应日期为所述拔节开始日期。
3.根据权利要求1所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述麦田逐小时草面温度Tt为气象站点观测到的逐小时草面温度数据，或将所述气象站点观测到的逐小时草面温度数据进行空间插值而生成逐小时草面温度空间分布数据。
4.根据权利要求3所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述空间插值方法采用样条函数方法。
5.根据权利要求1所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述合适的冻害温度阈值Th确定方法为：
测定所述的受冻冬小麦产量减产率；
设定所述冻害温度闽值Th的若干候选值为-6℃、-5℃、-4℃、-3℃、-2℃、-1℃和0℃，分别计算每个采样点上所述候选值对应的所述IFS指数；
将所述IFS指数和所述减产率进行线性拟合，计算不同所述冻害温度阈值Th条件下二者线性模型的决定系数R2和均方根误差RMSE，依据所述R2最大和所述RMSE最小的原则来确定所述冻害温度阈值Th。
6.根据权利要求5所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述受冻冬小麦产量减产率测定方法为：
对照产量Yn测定，通过田间抽样调查的办法，在研究区晚霜冻害事件发生年份，逐样点框定1米×1米的采样范围，并将全部茎蘖收割、装袋且带回实验室凉晒干；针对每个样点，统计全部所述茎蘖数量，以作为正常生长条件下冬小麦的穗数ENn；挑选出全部正常穗，统计其粒数和粒重，进而计算穗粒数和千粒重，以作为正常生长条件下冬小麦的穗粒数GNPEn和千粒重TGWn；根据公式ENn×GNPEn×TGWn，计算得出正常生长条件下冬小麦的产量Yn，作为所述对照产量；
实际产量Yf测定，统计全部茎蘖的籽粒重，作为所述实际产量Yf；
利用公式(Yn-Yf)/Yn×100％，计算得出所述受冻冬小麦产量减产率。
7.根据权利要求5所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，所述受冻冬小麦产量减产率测定方法为：
基于多年测定的冬小麦产量数据，将实际产量分解为趋势产量、气象产量和随机产量3个部分，所述随机产量因其影响小而忽略不计，将所述趋势产量作为正常产量，所述实际产量和所述趋势产量之差计算所述受冻冬小麦产量减产率。
8.根据权利要求1所述的拔节后受冻冬小麦产量损失早期预测方法，其特征在于，针对某一地区区域内的所述受冻冬小麦，采用基于面向对象分类、深度学习等方法的遥感分类技术，提取出冬小麦空间分布信息，屏蔽非小麦信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111013108.0A CN113592201B (zh) | 2021-08-31 | 2021-08-31 | 一种拔节后受冻冬小麦产量损失早期预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111013108.0A CN113592201B (zh) | 2021-08-31 | 2021-08-31 | 一种拔节后受冻冬小麦产量损失早期预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113592201A CN113592201A (zh) | 2021-11-02 |
| CN113592201B true CN113592201B (zh) | 2023-06-20 |
Family
ID=78240668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111013108.0A Active CN113592201B (zh) | 2021-08-31 | 2021-08-31 | 一种拔节后受冻冬小麦产量损失早期预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113592201B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108898499A (zh) * | 2018-05-31 | 2018-11-27 | 河南省气象科学研究所 | 一种冬小麦干热风保险天气指数计算方法 |
| CN110063158A (zh) * | 2019-04-25 | 2019-07-30 | 山西省农业科学院小麦研究所 | 一种利用可移动日光温室鉴定小麦耐春霜冻的方法 |
| CN110909933A (zh) * | 2019-11-20 | 2020-03-24 | 北京师范大学 | 一种耦合作物模型与机器学习语言的农业干旱快速诊断和评估方法 |
| CN113218887A (zh) * | 2021-06-23 | 2021-08-06 | 河南农业大学 | 一种小麦冻害严重程度快速监测评估方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010149504A2 (en) * | 2009-06-22 | 2010-12-29 | Basf Se | Method for enhancing harvest security of crops requiring vernalization |
-
2021
- 2021-08-31 CN CN202111013108.0A patent/CN113592201B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108898499A (zh) * | 2018-05-31 | 2018-11-27 | 河南省气象科学研究所 | 一种冬小麦干热风保险天气指数计算方法 |
| CN110063158A (zh) * | 2019-04-25 | 2019-07-30 | 山西省农业科学院小麦研究所 | 一种利用可移动日光温室鉴定小麦耐春霜冻的方法 |
| CN110909933A (zh) * | 2019-11-20 | 2020-03-24 | 北京师范大学 | 一种耦合作物模型与机器学习语言的农业干旱快速诊断和评估方法 |
| CN113218887A (zh) * | 2021-06-23 | 2021-08-06 | 河南农业大学 | 一种小麦冻害严重程度快速监测评估方法 |
Non-Patent Citations (3)
| Title |
|---|
| 农业旱灾脆弱性研究综述;别得进;《北京师范大学学报》;全文 * |
| 基于多因子关联的冬小麦晚霜冻害分析—以河南省为例;朱虹晖;《中国农业气象》;全文 * |
| 晚霜冻胁迫后冬小麦株高降低及其与籽粒产量关系;武永峰;《中国农业科学》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113592201A (zh) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103955606B (zh) | 一种基于遥感技术的草原蝗灾渐进式预测方法 | |
| CN109800921B (zh) | 基于遥感物候同化和粒子群优化的区域冬小麦估产方法 | |
| CN110751094A (zh) | 一种基于gee综合遥感影像和深度学习方法的作物估产技术 | |
| JP4719893B2 (ja) | 制御装置、制御方法、及びそのプログラム | |
| CN115099453B (zh) | 一种多变量栅格化卫星遥感综合干旱风险评估方法 | |
| Robertson et al. | Downscaling of seasonal precipitation for crop simulation | |
| CN108376265A (zh) | 一种冬小麦晚霜冻害多致灾因子权重的确定方法 | |
| CN110246051A (zh) | 作物生长周期模型及不同种类农产品建立生长周期的方法 | |
| CN109472283A (zh) | 一种基于多重增量回归树模型的危险天气事件预测方法和装置 | |
| CN113779113B (zh) | 基于雨洪时空过程相似性挖掘的洪水动态预估方法及系统 | |
| Wu et al. | Predicting yield loss in winter wheat due to frost damage during stem elongation in the central area of Huang-huai plain in China | |
| CN113592201B (zh) | 一种拔节后受冻冬小麦产量损失早期预测方法 | |
| CN111579565B (zh) | 农业干旱监测方法、系统及存储介质 | |
| Wang et al. | Exploring the feasibility of winter wheat freeze injury by integrating grey system model with RS and GIS | |
| CN116629453A (zh) | 一种适用于作物全生育时期的遥感估产方法 | |
| Spennemann et al. | Assessment of seasonal soil moisture forecasts over Southern South America with emphasis on dry and wet events | |
| CN109615149A (zh) | 一种确定甜菜气象产量的方法和系统 | |
| CN114936765A (zh) | 一种考虑空间异质性的农业干旱指数构建方法 | |
| Bagis et al. | An adaptive spatiotemporal agricultural cropland temperature prediction system based on ground and satellite measurements | |
| CN114994799A (zh) | 基于叶片尺度温度的棉花霜冻预报方法及系统 | |
| CN112561140A (zh) | 基于东亚副热带急流和极锋急流协同变化的中国四季降水预测方法 | |
| CN115965245B (zh) | 一种基于机器学习的山地果园冻害风险预测方法 | |
| Gill et al. | Pre-harvest wheat yield prediction using CERES-wheat model for Ludhiana district, Punjab, India | |
| CN116756591B (zh) | 一种基于水胁迫条件估产模型的油茶遥感估产方法 | |
| Pompei et al. | Validation of retrieval algorithms for the infrared remote sensing of precipitation with the Sardinian rain gauge network data |
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 |