首页 > 分享 > 一种快速鉴定小麦节水抗旱性的方法

一种快速鉴定小麦节水抗旱性的方法

花匠小妙招
2025-09-11 03:26

1.本发明涉及冬小麦品种节水抗旱性鉴定领域，尤其涉及一种利用遥感植被指数ndvi快速鉴定小麦节水抗旱性的方法。

背景技术：

2.水分在干旱和半干旱生态环境中，是限制粮食产量的主要因素，严重威胁着全球粮食安全和社会稳定。据有关统计，全球有大约三分之一的可耕种土地位于干旱和半干旱地区，涉及20多个国家和地区(sivakumar et al.2005)。在中国，干旱区域的总面积接近3百万平方公里，约占陆地面积的近4成，其中极度干旱的国土面积约有70万平方公里，毫无疑问，干旱已成为中国最主要自然灾害。小麦是世界分布范围最广、栽培面积最大的口粮作物，在口粮消费中占比超过60％，具有难以替换的优势地位(王一杰2018)。在中国，小麦是仅次于水稻和玉米的第三大粮食作物，种植区域主要集中于北方干旱半干旱地区。常年依靠灌溉维持干旱半干旱地区的作物产量，一段时间内，确实起到了保障了国家粮食安全，实现了由“南粮北调”向“北粮南运”的历史性转变。然而，深层土壤水的过度抽取使得北方地下水漏斗区不断加深和扩大，严重破坏自然生态，不可持续，地下水超采综合治理已刻不容缓，大幅度降低灌溉量甚至休耕都可能是各区域因地制宜的备选方案。因此，在未来有限灌溉甚至雨养条件下，确保有限土地生产更多粮食是唯一必走之路。研究小麦的节水抗旱性规律，快速挖掘和筛选优良种质资源并推广应用将是未来粮食生产的重要保障之一。
3.小麦的抗旱性可以从植株形态性状，生理生化性状以及产量性状这三个角度来鉴定(陈晓杰2013)。基于产量的抗旱性鉴定是一种综合性的鉴定方法，是小麦植株安全度过生命周期完成其生活史的总体性评价。目前，基于产量的抗旱系数(fischer and maurer 1978)和抗旱指数(sio
‑
se mardeh et al.2006)得到广泛采纳和应用并取得了较好的效果。但是在获取抗旱信息时，需要实地测量小麦的产量等数据，比较费时费力，这就限制了品种筛选的规模和速度。因此，寻求可以替代产量的抗旱性鉴定指标是一个较好的突破口。
4.ndvi作为遥感无损监测领域常用的植被指数，最早由rouse等人于1973年提出，在作物干物质(hansen and schjoerring 2003)，叶面积(qiu et al.2018)，叶绿素(gonzalez
‑
dugo et al.2015)，叶片氮素(shi et al.2015)等生理生化指标预测上，在作物营养诊断(lu et al.2017)，推荐施肥(maes and steppe 2019)，产量预估(nagy et al.2018)等实际应用中，都能找到ndvi的身影。但是利用ndvi指标来代替小麦产量来进行抗旱性鉴定还未见报道。

技术实现要素：

5.本发明的目的是通过获取小麦群体ndvi值来提供一种简单、快速、实用的小麦节水抗旱性鉴定的方法。本发明的方法与其它鉴定小麦节水抗旱性方法的显著区别在本发明是一种无损节水抗旱性鉴定。
6.本发明所提供的鉴定小麦节水抗旱性的方法，包括如下步骤：
7.1)在干旱或半干旱地区设置小麦处于雨养环境、限水灌溉环境以及水分充足灌溉环境下的品种节水抗旱性比较试验；
8.2)通过田间调查，记录不同小麦群体的具体开花日期，精确到某一自然日；
9.3)获取小麦群体开花后ndvi值时间序列动态；
10.4)构造ndvi相关指标
11.基于步骤3)构造ndvi相关指标mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
或者mndvi
x
；
12.其中，mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
分别表示花后0
‑
20天内的每隔1天、2天、4天、5天、10天测一次ndvi并求平均值，mndvi
x
表示花后籽粒形成期(花后0
‑
6天)某一天的ndvi值和灌浆中期(花后15
‑
20天)某一天的ndvi值的平均值；
13.5)基于ndvi相关指标计算评价小麦品种节水抗旱性的指标：抗旱系数drc(drought resistance coefficient)、抗旱指数di(drought resistance index)和节水指数wsi(water
‑
saving index)；
[0014][0015][0016][0017]
其中，n
a
是参试品种0水处理(雨养环境)的ndvi指标；n
m
是参试品种丰水处理(水分充足灌溉环境)的ndvi指标；
[0018]
n
a
是对照品种0水处理的ndvi指标；n
m
是对照品种丰水处理的ndvi指标；
[0019]
n
b
是参试品种节水处理的ndvi指标；n
b
是对照品种限水处理(限水灌溉环境)的ndvi指标；
[0020]
所述ndvi指标可为mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
或mndvi
x
；具体可为mndvi
x
；
[0021]
6)以某一节水抗旱性能已知的小麦品种(系)为对照，基于ndvi计算得到的抗旱系数drc、抗旱指数di、或节水指数wsc，将待测品种基于相同的ndvi指标计算得到的抗旱系数drc、抗旱指数di或节水指数wsi与之进行比较，确定待测品种的节水抗旱性能。
[0022]
上述方法步骤1)中，所述雨养环境、限水灌溉环境以及水分充足灌溉环境分别为：小麦足墒播种后，至收获前不再补充灌溉；小麦足墒播种后，春季灌拔节水；小麦足墒播种后，春季灌拔节水和开花水，每次补充灌溉水75mm；
[0023]
小麦品种由育种单位提供或者由试验单位通过商业途径收集；试验地土壤养分均匀，地势平坦，灌溉水盐碱度较低；
[0024]
上述方法步骤2)中，开花日期的确定标准为：开花期为麦田50％麦穗的中部花开放，花药开始露出的日期；
[0025]
上述方法步骤3)中，使用双波段光谱仪greenseeker(trimble，usa)测定冠层群体的ndvi值；
[0026]
ndvi值由红光波段的反射率和近红外光波段反射率计算得出，计算公式为：ndvi
＝(nir
‑
red)/(nir+red)，nir为近红外光波段反射率，red为红光波段反射率。
[0027]
ndvi相关指标在小麦节水抗旱性快速鉴定中的应用也属于本发明的保护范围。
[0028]
ndvi相关指标为：mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
或mndvi
x
；具体可为mndvi
x
；
[0029]
其中，mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
分别表示花后0
‑
20天内的每隔1天、2天、4天、5天、10天测一次ndvi并求平均值，mndvi
x
表示花后籽粒形成期(花后0
‑
6天)某一天的ndvi值和灌浆中期(花后15
‑
20天)某一天的ndvi值的平均值；
[0030]
利用ndvi相关指标进行小麦节水抗旱性快速鉴定在小麦品种筛选和育种中的应用也属于本发明的保护范围。
[0031]
ndvi相关指标为：mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
或者mndvi
x
；
[0032]
其中，mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
分别表示花后0
‑
20天内的每隔1天、2天、4天、5天、10天测一次ndvi并求平均值，mndvi
x
表示花后籽粒形成期(花后0
‑
6天)某一天的ndvi值和灌浆中期(花后15
‑
20天)某一天的ndvi值的平均值。
[0033]
本发明提供了一种基于ndvi相关指标对小麦节水抗旱性进行快速鉴定的方法，与现有的测定方法相比，具有以下明显的优点：
[0034]
1)本发明具有操作简单、快速、成本低、效率高等优点；无需实测小麦产量；
[0035]
2)本发明在操作过程中不使用有毒化学试剂，对环境无污染且无需破坏性取样，可以做到对小麦冠层无损害；
[0036]
3)本发明在获取小麦群体ndvi的时间上较为灵活，可以避免阴雨天等不利环境对小麦冠层测试的影响。
附图说明
[0037]
图1为本发明的利用遥感植被指数ndvi快速鉴定小麦节水抗旱性的示意图。
[0038]
图2为基于不同采样频率下ndvi平均值计算得到的抗旱系数和抗旱指数与实际抗旱系数和抗旱指数的线性关系。
[0039]
图3为基于花后20天内任意两天ndvi平均值计算得到的抗旱系数和抗旱指数/节水指数与实际抗旱系数和抗旱指数/节水指数的线性关系。
具体实施方式
[0040]
下面通过具体实施例对本发明进行说明，但本发明并不局限于此。
[0041]
下述实施例中所使用的实验方法如无特殊说明，均为常规方法；下述实施例中所用的试剂、材料等，如无特殊说明，均可从商业途径得到。
[0042]
实施例1、冬小麦节水抗旱性快速鉴定
[0043]
1、试验设置
[0044]
抗旱性品种筛选试验于2017
‑
2019年中国农业大学吴桥实验站进行。试验所用品种包括：衡观35、衡s29、衡136、衡4399、衡4444、衡6632、中麦23、中麦155、中麦175、中麦1062、农大139、农大399、农大3636、农大3753、农大5181、藁优2018、藁优5766、藁优5818、刑麦18、刑麦20、石麦22、济麦22、小偃81、临旱6号、山农32、周麦18、良星66、石农086、冀麦418和沧6005，供试品种材料由中国农业大学农学院提供。
[0045]
2、从小麦抽穗后开始，每日上午9点去试验地观察小麦开花情况，记录不同品种具体开花日期。
[0046]
3、小麦群体花后ndvi值时间序列动态的获取
[0047]
(1)使用双波段光谱仪greenseeker(trimble，usa)测定冠层群体的ndvi值(ndvi值由红光波段的反射率和近红外光波段反射率计算得出，计算公式为：ndvi＝(nir
‑
red)/(nir+red)，nir为近红外光波段反射率，red为红光波段反射)；
[0048]
(2)测试条件：晴朗无云或少云，无风或微分环境；测试人员身作深色衣帽；上午11点至下午1点；传感器探头距离冠层上方70
‑
90cm并保持水平；测试人员要避免身影遮挡被测冠层区域顶部；连续测试长度超过3米并往返各1次并求平均值作为每次测试的记录值；
[0049]
(3)阴雨天未测试的冠层ndvi数据使用数学插值法埃尔米特(hermite)插值法补齐。
[0050]
4、构建ndvi相关指标
[0051]
mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
和mndvi
y
；
[0052]
mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
分别表示花后0
‑
20天内的每隔1天、2天、4天、5天、10天测一次ndvi并求平均值，mndviy表示花后0
‑
20天内任意两天ndvi的平均值；
[0053]
5、获取产量
[0054]
小麦成熟后，在测试ndvi的相同区域选取2m2，单独收割脱粒测产(换算成含水量为12.5％的产量)。
[0055]
6、使用产量数据计算不同品种的节水抗旱性指标drc、di和wsi
[0056][0057][0058][0059]
上述各式中，ya是参试品种0水处理的产量；ym是参试品种丰水处理的产量；ya是对照品种0水处理的产量；ym是对照品种丰水处理的产量。yb是参试品种节水处理的产量；yb是对照品种节水处理的产量。本实施例对照数据采用的是当年试验中相同处理下所有品种产量的平均值。
[0060]
7、使用步骤4创建出的ndvi相关指标按照如下公式计算出不同品种的基于ndvi相关指标的节水抗旱性指标；
[0061][0062][0063][0064]
其中，n
a
是参试品种0水处理的ndvi指标；n
m
是参试品种丰水处理的ndvi指标；
[0065]
n
a
是对照品种0水处理的ndvi指标；n
m
是对照品种丰水处理的ndvi指标；
[0066]
n
b
是参试品种限水处理的ndvi指标；n
b
是对照品种限水处理的ndvi指标；
[0067]
所述ndvi指标为mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
或mndvi
y
；
[0068]
8、使用步骤7计算出的基于ndvi相关指标的节水抗旱性指标和步骤6计算处理的实际抗旱性指标进行线性回归分析；
[0069]
9、挑选出步骤8中参与线性拟合决定系数较高的ndvi指标作为快速节水抗旱性鉴定的备选指标；
[0070]
10、选取某一节水抗旱性能已知的小麦品种为对照。选用步骤9确定的ndvi备选指标，基于该ndvi备选指标计算得到对照品种的抗旱系数、抗旱指数di和节水指数wsi，获取待测品种的ndvi备选指标，基于待测品种的ndvi备选指标计算得到待测品种的抗旱系数、抗旱指数di和节水指数wsi，将待测品种计算得到的抗旱系数、抗旱指数di和节水指数wsi与对照品种的抗旱系数、抗旱指数di和节水指数wsi分别进行比较，确定待测品种的节水抗旱性能。
[0071]
线性回归分析结果如图2和图3，mndvi、mndvi
2d
、mndvi
4d
、mndvi
5d
、mndvi
10d
以及mndvix可以作为备选指标。
[0072]
mndvix是mndviy的一部分花后ndvi组合，表示籽粒形成期(花后0
‑
6天)某一天的ndvi值和灌浆中期(花后15
‑
20天)某一天的ndvi值的平均值。
[0073]
图3呈现的是以mndviy为指标得到的线性拟合决定系数r2热图。以红颜色为基准，只有以mndvix组合为指标得到的线性拟合决定系数r2达到0.8以上(各节水抗旱指标同时满足)。因此，可选定mndvix组合的ndvi指标来确定待测品种的节水抗旱性能。
[0074]
表1、表2和表3展现的是基于ndvi指标和基于产量计算出来的抗旱系数、抗旱指数和节水指数具有较好的同步性。
[0075]
表1基于产量和ndvi指标的抗旱系数比较
[0076][0077]
表2基于产量和ndvi指标的抗旱指数比较
[0078][0079]
表3基于产量和ndvi指标的节水指数比较
[0080][0081]
虽然，上文中已经用一般性说明、具体实施方式及试验，对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。