• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于無人機多源影像數(shù)據(jù)的灌漿期人工合成小麥抗旱性評價

    2024-05-14 03:24:00燕雯金秀良李龍徐子涵蘇悅張躍強景蕊蓮毛新國孫黛珍
    中國農(nóng)業(yè)科學(xué) 2024年9期
    關(guān)鍵詞:抗旱性抗旱種質(zhì)

    燕雯,金秀良,李龍,徐子涵,蘇悅,張躍強,景蕊蓮,毛新國,孫黛珍

    基于無人機多源影像數(shù)據(jù)的灌漿期人工合成小麥抗旱性評價

    燕雯1,2,金秀良2,李龍2,徐子涵3,蘇悅1,2,張躍強3,景蕊蓮2,毛新國2,孫黛珍

    1山西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,山西太谷 030801;2中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所/作物基因資源與育種全國重點實驗室,北京 100081;3新疆農(nóng)業(yè)科學(xué)院核技術(shù)生物技術(shù)研究所/農(nóng)業(yè)農(nóng)村部荒漠綠洲區(qū)作物生理生態(tài)與耕作重點實驗室/新疆作物生物技術(shù)重點實驗室,烏魯木齊 830091

    【目的】基于無人機多源影像及產(chǎn)量數(shù)據(jù)評價人工合成小麥種質(zhì)的抗旱性,優(yōu)選高通量抗旱性鑒定指標(biāo),發(fā)掘抗旱人工合成小麥種質(zhì)資源,為加快拓展小麥抗旱遺傳資源、提升旱地小麥育種水平提供技術(shù)支撐和種質(zhì)材料?!痉椒ā恳?0份人工合成小麥種質(zhì)及對照小麥品種新春37為試驗材料,在田間進行小區(qū)播種,設(shè)置干旱和灌溉2種水分處理;利用無人機搭載多光譜及熱紅外相機采集試驗材料灌漿期多源影像進行拼接處理,通過閾值分割等方法提取各試驗材料的光譜指數(shù);利用相關(guān)性分析和主成分分析鑒選抗旱相關(guān)光譜指標(biāo),結(jié)合單指標(biāo)及綜合評價方法鑒定人工合成小麥種質(zhì)的抗旱性?!窘Y(jié)果】基于無人機多源影像數(shù)據(jù)提取了80份人工合成小麥種質(zhì)的19種光譜指數(shù)。不同光譜指數(shù)抗旱系數(shù)與小區(qū)產(chǎn)量抗旱指數(shù)的相關(guān)性分析結(jié)果表明,OSAVI的抗旱系數(shù)與抗旱指數(shù)的關(guān)聯(lián)度最高,NDVI、CIre和NDRE的抗旱系數(shù)與抗旱指數(shù)的關(guān)聯(lián)度較高。部分光譜指數(shù)的抗旱系數(shù)間相關(guān)性較高,存在冗余信息,通過主成分分析,將19個光譜指數(shù)的抗旱系數(shù)轉(zhuǎn)換為3個相互獨立的綜合指標(biāo),3個綜合指標(biāo)的貢獻度分別為59.6%、12.0%和9.6%。利用加權(quán)隸屬函數(shù)法聚合綜合指標(biāo),通過公式計算獲得各人工合成小麥種質(zhì)的綜合抗旱性度量值?;诳购抵笖?shù)鑒定出6份強抗旱人工合成小麥種質(zhì),基于綜合抗旱性度量值鑒定出5份強抗旱種質(zhì),其中,SW004和SW009在2種方法的評價結(jié)果中均被評為強抗旱種質(zhì)?;贠SAVI的抗旱系數(shù)對80份人工合成小麥種質(zhì)進行抗旱性分級,分級結(jié)果與基于綜合抗旱性度量值的分級結(jié)果基本一致。根據(jù)OSAVI的抗旱系數(shù)鑒定出的6份強抗旱種質(zhì)中,有5份在基于綜合抗旱性度量值分級中也被鑒定為強抗旱種質(zhì)?!窘Y(jié)論】基于無人機多源影像提取的光譜指數(shù)NDVI、OSAVI、CIre和NDRE,以及基于光譜指數(shù)的綜合抗旱性度量值均可用于輔助鑒定小麥種質(zhì)抗旱性。

    多源影像;光譜指數(shù);人工合成小麥;抗旱性;灌漿期

    0 引言

    【研究意義】小麥?zhǔn)俏覈狈街饕Z食作物,也是農(nóng)業(yè)用水大戶[1]。在水資源短缺的背景下,選育抗旱節(jié)水品種是應(yīng)對環(huán)境脅迫、維持小麥生產(chǎn)可持續(xù)發(fā)展、保障糧食安全的重要途徑[2]。受長期人工選擇的影響,現(xiàn)代小麥品種的遺傳本底逐漸趨于同質(zhì)化。由四倍體小麥和粗山羊草經(jīng)過人為雜交及染色體加倍得到的人工合成小麥,蘊含豐富的野生近緣種基因資源,遺傳變異范圍廣,對不良環(huán)境適應(yīng)性強,是改良現(xiàn)代小麥的重要基因資源,例如,利用人工合成小麥育成的川麥42、龍春1號等品種抗病、抗逆、適應(yīng)性廣。因此,發(fā)掘并利用人工合成小麥抗旱種質(zhì)資源,有助于拓寬小麥栽培品種的遺傳基礎(chǔ),提高小麥抗旱育種核心材料的多樣性[3],加快突破性抗旱小麥新品種的培育進程。表型性狀反映了作物對環(huán)境條件的響應(yīng)[4]。傳統(tǒng)的作物表型獲取方法大多需要人工操作,費時費力且準(zhǔn)確度較低,限制了優(yōu)異種質(zhì)資源的篩選效率。得益于現(xiàn)代農(nóng)業(yè)技術(shù)的發(fā)展,無人機高通量表型平臺以其靈活高效的特點正在被應(yīng)用于田間作物表型信息的采集。發(fā)掘抗旱相關(guān)高通量表型鑒定指標(biāo),用于描述和評估抗旱表型特征,是高效發(fā)掘和利用人工合成小麥種質(zhì)資源的有效途徑?!厩叭搜芯窟M展】作物的抗旱性是由多基因控制的數(shù)量性狀,評價較為困難[5]。國家標(biāo)準(zhǔn)《小麥抗旱性鑒定評價技術(shù)規(guī)范》(GB/T 21127—2007)明確了小麥抗旱性鑒定技術(shù)和評價標(biāo)準(zhǔn),例如種子萌發(fā)期采用高滲溶液法、苗期采用反復(fù)干旱法評價抗旱性,水分臨界期和全生育期的抗旱性則直接通過產(chǎn)量計算抗旱指數(shù)。灌漿期是小麥籽粒形成的關(guān)鍵時期,水分脅迫影響籽粒中光合產(chǎn)物的積累,導(dǎo)致減產(chǎn),是抗旱性研究中具有代表性的時期。白志英等[6]通過主成分分析和聚類分析方法,對干旱脅迫下21個小麥代換系進行了抗旱性鑒定評價,將其劃分為高度抗旱類型、中度抗旱類型和不抗旱類型;王士強等[7]采用灰色關(guān)聯(lián)度分析法分析了65個小麥品種的抗旱指數(shù)與其主要農(nóng)藝性狀和部分生理生化性狀,提出葉片氣孔相關(guān)的生理生化性狀與抗旱指數(shù)關(guān)聯(lián)度最大,可以選擇相關(guān)性狀進行小麥抗旱性的鑒定。張軍等[8]檢測了8個強筋小麥品種萌發(fā)期和幼苗期的抗旱相關(guān)生理指標(biāo),通過主成分分析和隸屬函數(shù)法評價其抗旱性,篩選出1個在萌發(fā)期和幼苗期均表現(xiàn)抗旱的品種;馬玉慧等[9]通過主成分分析及加權(quán)隸屬函數(shù)法對70份春小麥材料在萌發(fā)期的生理指標(biāo)進行了比較,提出利用綜合抗旱指數(shù)評價春小麥抗旱性。隨著作物表型組學(xué)的發(fā)展,高通量表型鑒定方法逐漸被應(yīng)用于小麥抗旱性鑒定評價。孟雨等[10]借助手持式熱紅外儀檢測技術(shù),明確了冠層溫度特征參數(shù)與抗旱指數(shù)間的定量關(guān)系,構(gòu)建了基于冠層溫度特征參數(shù)的抗旱指數(shù)回歸模型;Qin等[1]獲取了小麥在不同生育期的無人機采集圖像和生理性狀,結(jié)合層次聚類分析方法評價10個小麥品種的抗旱性,明確了在缺水條件下,小麥冠層溫度與葉面積指數(shù)、葉綠素含量、光合速率、蒸騰作用和氣孔導(dǎo)度呈負相關(guān)。DAS等[11]通過無人機遙感量化了基因型與環(huán)境的相互作用,評估不同基因型栽培小麥的耐旱性,篩選出在堿土上最耐旱的小麥基因型Bremer。傳統(tǒng)的小麥抗旱性鑒定方法準(zhǔn)確性較高,但鑒定時間較長,不適用于大田環(huán)境中對大量材料進行抗旱性鑒定。與之相比,高通量表型鑒定方法具有快速、高效、不受小麥生育時期限制的優(yōu)點,同時減少了傳統(tǒng)方法中耗時費力、效率較低,以及人為主觀影響的問題。迄今為止,研究者利用高通量檢測方法評價抗旱性的對象多為普通小麥,尚缺少針對人工合成小麥抗旱性評價的報道。【本研究切入點】與普通小麥相比,人工合成小麥遺傳變異范圍廣,表型多樣性豐富,現(xiàn)有的小麥抗旱性鑒定評價指標(biāo)體系是否適用于人工合成小麥,尚需實踐檢驗,因此,亟須研究適宜的高通量鑒定技術(shù)及評價指標(biāo),為從人工合成小麥中發(fā)掘抗旱優(yōu)異種質(zhì)資源提供技術(shù)支持?!緮M解決的關(guān)鍵問題】本研究基于無人機多源影像數(shù)據(jù),通過主成分分析和隸屬函數(shù)法建立從高通量光譜指數(shù)中優(yōu)選適用于人工合成小麥的抗旱性評價指標(biāo),發(fā)掘優(yōu)異抗旱種質(zhì)資源,為加快拓展小麥抗旱遺傳資源、提升旱地小麥育種水平提供技術(shù)支撐和種質(zhì)材料。

    1 材料與方法

    1.1 試驗材料與田間種植

    以80份從國際玉米小麥改良中心(CIMMYT)引進的人工合成小麥種質(zhì)(附表1)為試驗材料,以小麥品種新春37為對照品種。新春37由新疆農(nóng)業(yè)科學(xué)院核技術(shù)生物技術(shù)研究所選育,是新疆種植面積最大的春小麥品種,生育期110 d左右,具有抗倒伏、綜合抗逆性好等特征。

    試驗在新疆維吾爾自治區(qū)昌吉回族自治州昌吉市軍戶農(nóng)場(43°96′N,87°01′E)進行,試驗基地海拔高度為1 266 m,年平均溫度6.8 ℃,土壤質(zhì)地為砂壤土,有機質(zhì)含量1.2%,播種前0—10 cm土壤平均含水量為19.5%。分別于2022年4月和2023年3月播種,2022年7月和2023年8月收獲。2個小麥生長季內(nèi)的自然降水量分別為42.4和60.8 mm。試驗設(shè)置干旱脅迫和灌溉對照2組處理,各2次重復(fù),隨機區(qū)組設(shè)計。每個小區(qū)6行,行長2 m,行距0.2 m,參考大田生產(chǎn)播種量,每行播240粒種子。干旱脅迫組與灌溉對照組均為10 d左右澆一次水,共澆水7—8次。全生育期內(nèi),干旱脅迫組灌水總量為2 700 m3·hm-2,對照組灌水總量為4 500 m3·hm-2。施肥及病蟲害防治等田間管理措施與當(dāng)?shù)匦←溕a(chǎn)模式相同。

    1.2 數(shù)據(jù)采集

    于小麥灌漿期選擇晴朗無風(fēng)的天氣進行2次遙感影像采集。采集系統(tǒng)由DJI Matrice M300型無人機搭載MicaScence RedEdge-MX多光譜相機和FLIR Duo Pro R640熱紅外相機組成。MicaScence RedEdge-MX相機包括藍光(blue)、綠光(green)、紅光(red)、紅邊(red edge,Re)和近紅外(near infrared,NIR)5個波段,中心波長分別為475、560、668、717和840 nm。每次無人機飛行前,使用與多光譜相機配套的MicaScence校正白板獲取參照圖像。無人機飛行高度為50 m,飛行速度為4.5 m·s-1,影像航向和旁向重疊率設(shè)置為80%,拍照時間間隔1 s。試驗地內(nèi)放置6塊帶有溫度傳感器的黑色鋁板,用于標(biāo)定實際溫度并構(gòu)建熱紅外圖像定標(biāo)公式。小麥開花后測量株高,收獲后稱取小區(qū)產(chǎn)量,采用自動考種儀測量千粒重。

    1.3 圖像處理

    使用Agisoft PhotoScan軟件進行多光譜圖像及熱紅外圖像的拼接與校正,生成試驗區(qū)域的正射影像。為避免邊緣效應(yīng)的影響,使用ArcGIS10.2軟件,在各小區(qū)的中間部位選定目標(biāo)區(qū)域。對于多光譜圖像,采用閾值分割法分離出各小區(qū)的植株部分,使用重分類工具制作掩膜后再利用柵格計算器工具進行運算,得到各波段去除土壤背景后的影像數(shù)據(jù),最終提取目標(biāo)區(qū)域在不同波段的多光譜圖像反射率(圖1)。對于熱紅外圖像,采用閾值分割法提取各小區(qū)的平均灰度值(digital number,DN)。

    1.4 數(shù)據(jù)分析

    根據(jù)國家標(biāo)準(zhǔn)《小麥抗旱性鑒定評價技術(shù)規(guī)范》(GB/T 2127—2007)計算供試材料的抗旱指數(shù)(drought resistance index,),公式如下:

    =S.T2×S.W-1×CK.W×(CK.T2)-1(1)

    式中,S.T和S.W分別為試驗材料在干旱脅迫和灌溉條件下的籽粒產(chǎn)量,CK.W和CK.T分別為對照品種(新春37)在干旱脅迫和灌溉條件下的籽粒產(chǎn)量。使用IBM SPSS Statistics 19.0軟件進行基本統(tǒng)計學(xué)分析。利用表1所述公式計算各項光譜指數(shù)。

    利用前期獲得的熱紅外圖像定標(biāo)公式(2)計算小區(qū)冠層溫度(canopy temperature,CT):

    表1 光譜指數(shù)計算公式

    BNDVI:藍色歸一化植被指數(shù);CIgreen:綠邊葉綠素指數(shù);CIre:紅邊葉綠素指數(shù);DVI:差值植被指數(shù);GBNDVI:藍-綠歸一化植被指數(shù);GCC:綠色色度坐標(biāo);GNDVI:綠色歸一化植被指數(shù);MCARI:修正葉綠素調(diào)節(jié)植被指數(shù);mNDblue:修正藍光歸一化植被指數(shù);MTCI:MERIS陸地葉綠素指數(shù);NDRE:標(biāo)準(zhǔn)化差分紅邊指數(shù);NDVI:歸一化植被指數(shù);NGBDI:歸一化綠藍差異指數(shù);NPCI:歸一化色素葉綠素比值指數(shù);NRCT:歸一化相對冠層溫度;OSAVI:優(yōu)化土壤調(diào)整植被指數(shù);SIPI:結(jié)構(gòu)不敏感色素指數(shù);SR:簡單比值植被指數(shù);TVI:三角形植被指數(shù)。a:RBlue、RGreen、RRed、RRe、RNIR分別表示藍、綠、紅、紅邊、近紅外波段反射率。CTi表示第i個種質(zhì)的冠層溫度;CTmin表示80份種質(zhì)最低冠層溫度;CTmax表示80份種質(zhì)最高冠層溫度。下同

    BNDVI: Blue normalized difference vegetation index; CIgreen: Green chlorophyll index; CIre: Red-edge chlorophyll index; DVI: Difference vegetation index;GBNDVI: Normalized blue-green band difference vegetation index; GCC: Green chromatic coordinate; GNDVI: Green-band normalized vegetation index; MCARI: Modified chlorophyll absorption ratio index; mNDblue: Modified normalized difference blue index; MTCI: Meris terrestrial chlorophyll index; NDRE: Normalized difference red edge index; NDVI: Normalized difference vegetation index; NGBDI: Normalized green-blue difference index; NPCI: Normalized pigment chlorophyll vegetation index; NRCT: Normalized relative canopy temperature; OSAVI: Optimized soil-adjusted vegetation index; SIPI: Structure insensitive pigment index; SR: Simple ratio index; TVI: Triangle vegetation index.a:RBlue, RGreen, RRed, RRe, RNIRrepresent the spectral reflectance of blue, green, red, red edge, and near infrared bands, respectively. CTirepresents the canopy temperature of the i-th genotype, CTminrepresents the lowest canopy temperature among all genotypes, and CTmaxrepresents the highest canopy temperature among all genotypes. The same as below

    圖1 圖像處理步驟

    CT=0.0353×DN-233.5 (2)

    采用公式(3)計算各光譜指數(shù)的抗旱系數(shù)(drought resistance coefficient,):

    =X/CK(=1,2,…n) (3)

    式中,X、CK分別表示干旱脅迫、灌溉條件下的光譜指數(shù)值;為光譜指數(shù)指標(biāo);為光譜指數(shù)數(shù)量。

    使用SPSS 19.0軟件對數(shù)據(jù)進行主成分分析(principal component analysis,PCA),利用公式(4)計算各綜合指標(biāo)的隸屬函數(shù)值():

    (X)=(X-min)/(max-min)=1,2,3,…,n (4)

    式中,X表示第個綜合指標(biāo),min表示第個綜合指標(biāo)的最小值,max表示第個綜合指標(biāo)的最大值。

    采用公式(5)計算各綜合指標(biāo)的權(quán)重():

    式中,W為第個綜合指標(biāo)在所有指標(biāo)中的權(quán)重,P表示各人工合成小麥材料中第個綜合指標(biāo)的貢獻率。

    利用公式(6)計算各人工合成小麥材料的綜合抗旱性度量值():

    式中,為綜合指標(biāo)數(shù)量;為綜合指標(biāo)。

    2 結(jié)果

    2.1 干旱脅迫對人工合成小麥農(nóng)藝性狀的影響

    在干旱脅迫和灌溉條件下(圖2-a),2022和2023年的試驗結(jié)果相似,在干旱脅迫和灌溉條件下,對照品種新春37的平均產(chǎn)量都高于人工合成小麥(圖2-b)。2022年,對照品種在干旱脅迫下的產(chǎn)量為 6 250.5 kg·hm-2,比灌溉條件下減少916.5 kg·hm-2,減產(chǎn)12.8%;而人工合成小麥在干旱脅迫下的產(chǎn)量為4 722.2 kg·hm-2,較灌溉條件下減少337.3 kg·hm-2,減產(chǎn)6.7%。2023年,與灌溉條件相比,干旱脅迫條件下對照品種減產(chǎn)17.7%,人工合成小麥減產(chǎn)21.7%。干旱脅迫也導(dǎo)致人工合成小麥株高和千粒重降低(圖2-c—d),2年平均株高降低5.2%,千粒重減少3.4%。人工合成小麥群體中不同農(nóng)藝性狀呈現(xiàn)連續(xù)分布且分布范圍較廣。2023年,在灌溉條件下,株高、千粒重和產(chǎn)量的分布范圍分別為70.7—118.0 cm、29.2—50.5 g和1 959.3—8 754.4 kg·hm-2,變異系數(shù)依次為10.6%、11.2%和30.6%;而在干旱脅迫下,3個性狀的變異系數(shù)依次為11.3%、12.3%和29.4%,變異系數(shù)均較大(圖2-c—e)。結(jié)果表明,無論在灌溉條件下,還是干旱脅迫條件下,人工合成小麥的農(nóng)藝性狀都表現(xiàn)出較大差異,尤其是在干旱脅迫條件下不同種質(zhì)各性狀的離散程度仍然較大,反映了不同種質(zhì)適應(yīng)水分虧缺的能力差異。

    2.2 單指標(biāo)抗旱性評價

    根據(jù)干旱和灌溉條件下的產(chǎn)量計算抗旱指數(shù),結(jié)果表明,人工合成小麥的平均抗旱指數(shù)為0.73,分布范圍為0.34—1.30,變異系數(shù)為28.1%,進一步表明,不同種質(zhì)的抗旱性存在較大差異。采用平均值-標(biāo)準(zhǔn)差分類法[29]對不同材料抗旱性進行分級,即使用群體材料的抗旱指數(shù)平均值加或減1.5個標(biāo)準(zhǔn)差,將材料分為強抗、中間類型和敏感材料;再使用抗旱指數(shù)平均值加或減0.5個標(biāo)準(zhǔn)差,將中間類型材料劃分為抗旱、中抗和較敏感材料?;诖朔椒ǎ瑢?0份人工合成小麥的抗旱性分為強抗(≥1.04)、抗(1.04>≥0.83)、中抗(0.83>≥0.63)、較敏感(0.63>≥0.42)和敏感(<0.42)5個等級。各等級抗旱性材料的數(shù)目依次為6、17、29、24和4份,分別占總數(shù)的7.5%、21.3%、36.3%、30.0%和5.0%(圖3-a)。其中,6份強抗旱種質(zhì)SW004、SW009、SW019、SW034、SW037和SW051的抗旱指數(shù)依次為1.30、1.24、1.23、1.10、1.08和1.10(表2)。

    2.3 抗旱指數(shù)與光譜指數(shù)抗旱系數(shù)的相關(guān)性分析

    抗旱指數(shù)與NDVI、OSAVI等14個光譜指數(shù)的抗旱系數(shù)(光譜指數(shù)-)顯著正相關(guān),與mNDblue、NPCI和SIPI的抗旱系數(shù)呈顯著負相關(guān)關(guān)系(圖4)。其中,CIre、NDRE、NDVI和OSAVI的抗旱系數(shù)與抗旱指數(shù)的相關(guān)性較高,均大于0.68,說明基于無人機多源影像數(shù)據(jù)提取的光譜指數(shù)的抗旱系數(shù)能在較大程度上反映不同種質(zhì)的抗旱性。然而,不同光譜指數(shù)抗旱系數(shù)間普遍存在較高的相關(guān)性,如,BNDVI、CIre、GCC、NDRE、OSAVI的抗旱系數(shù)與大多數(shù)光譜指數(shù)抗旱系數(shù)顯著正相關(guān),CIre-與NDRE-的相關(guān)系數(shù)高達0.99;NPCI和SIPI的抗旱系數(shù)與多個光譜指數(shù)抗旱系數(shù)顯著負相關(guān),且SIPI的抗旱系數(shù)與大部分光譜指數(shù)抗旱系數(shù)的相關(guān)系數(shù)超過0.60。由此可見,通過不同波段光譜參數(shù)所計算的各類光譜指數(shù)間存在大量的信息重疊,需要對光譜指數(shù)指標(biāo)進行篩選與優(yōu)化。

    a:基于抗旱指數(shù)的評價結(jié)果;b:基于綜合抗旱性度量值的評價結(jié)果;c:不同指標(biāo)評價結(jié)果韋恩圖;d:強抗旱及抗旱人工合成小麥種質(zhì)信息。HR:強抗旱;R:抗旱;MR:中等抗旱;S:干旱較敏感;HS:干旱敏感。DRI-HR:由抗旱指數(shù)評級為強抗旱的種質(zhì);DRI-R:由抗旱指數(shù)評級為抗旱的種質(zhì);D-HR:由綜合抗旱性度量值評級為強抗旱的種質(zhì);D-R:由綜合抗旱性度量值評級為抗旱的種質(zhì);OSAVI-DC-HR/R:OSAVI-DC值評級為強抗旱及抗旱的種質(zhì)

    表2 強抗旱人工合成小麥種質(zhì)產(chǎn)量差異

    2.4 光譜指數(shù)抗旱系數(shù)主成分分析

    為了消除大量光譜指數(shù)抗旱系數(shù)中冗余的表型信息,優(yōu)選抗旱性評價指標(biāo),進一步對所提取的19個光譜指數(shù)抗旱系數(shù)進行主成分分析(表3)。結(jié)果表明,19個光譜指數(shù)抗旱系數(shù)可轉(zhuǎn)化為3個累積貢獻率大于80%且相互獨立的主成分,其特征值分別為11.31、2.28和1.81,貢獻率分別為59.6%、12.0%和9.6%,說明此3個主成分所包含的信息量可以反映19個原始光譜指數(shù)的大部分信息。第一主成分中,CIreNDRE、NDVI和OSAVI具有較高的載荷系數(shù)絕對值;第二主成分中,MTCI相對于其他光譜指數(shù)的載荷系數(shù)絕對值較高;第三主成分中,起主要作用的是GBNDVI和NPCI。該結(jié)果與相關(guān)性分析結(jié)果相似,即貢獻率最高的第一主成分中載荷系數(shù)最大的CIreNDRE、NDVI和OSAVI與抗旱指數(shù)的相關(guān)系數(shù)最高(圖4),可優(yōu)選為高通量抗旱性鑒定評價指標(biāo)。

    表3 光譜指數(shù)抗旱系數(shù)主成分分析

    2.5 抗旱性綜合評價

    使用公式(4)計算各綜合指標(biāo)的隸屬函數(shù)值,通過公式(5)分別計算各綜合指標(biāo)的權(quán)重,并利用公式(6)計算出各試驗材料的綜合抗旱性度量值()。80份人工合成小麥的平均值為0.42,分布范圍為0.12—0.72,變異系數(shù)為24.5%,小于單指標(biāo)評價的變異系數(shù)?;谥?,使用平均值-標(biāo)準(zhǔn)差分類法將80份人工合成小麥的抗旱性分為強抗(≥0.57)、抗(0.57>≥0.47)、中抗(0.47>≥0.37)、較敏感(0.37>≥0.26)和敏感(<0.26)5個等級。各等級種質(zhì)數(shù)量分別為5、17、35、18和5份,占總數(shù)的6.3%、21.3%、43.8%、22.5%和6.3%(圖3-b)。其中,強抗旱種質(zhì)為SW004、SW009、SW024、SW047和SW089,值依次為0.60、0.72、0.58、0.61、0.64。此外,8份種質(zhì)(SW013、SW017、SW040、SW042、SW68、SW100、SW106、SW123)在單指標(biāo)和綜合評價結(jié)果中均為抗旱,SW004和SW009均為強抗旱(圖4-d)。

    抗旱指數(shù)與綜合抗旱性度量值()呈顯著正相關(guān),相關(guān)系數(shù)為0.72(圖5),表明2種方法的評價結(jié)果具有較高的相似性。此外,各光譜指數(shù)的抗旱系數(shù)中,OSAVI-與抗旱指數(shù)的關(guān)聯(lián)度最高,進一步基于OSAVI-對人工合成小麥進行抗旱性分級。使用平均值-標(biāo)準(zhǔn)差分類法將80份人工合成小麥種質(zhì)抗旱性分為強抗(≥0.91)、抗(0.91>≥0.85)、中抗(0.85>≥0.78)、較敏感(0.78>≥0.72)和敏感(<0.72)5個等級。各等級種質(zhì)數(shù)量分別為6、19、30、22和3份,占總數(shù)的7.5%、23.8%、37.5%、27.5%和3.8%。結(jié)果表明,僅使用OSAVI-單指標(biāo)進行抗旱性分級的結(jié)果與基于綜合抗旱性度量值的分級結(jié)果基本一致(圖3-c)。根據(jù)OSAVI-鑒定出的6份強抗旱種質(zhì)中,有5份(SW004、SW009、SW024、SW047、SW089)在綜合評價體系中也被鑒定為強抗旱。因此,在大規(guī)模小麥種質(zhì)資源抗旱性鑒定中,也可使用OSAVI-作為輔助評判依據(jù)。

    圖5 綜合抗旱性度量值及OSAVI-DC與抗旱指數(shù)的相關(guān)性

    3 討論

    3.1 基于光譜指標(biāo)的灌漿期人工合成小麥抗旱性綜合評價

    小麥的灌漿期是決定籽粒產(chǎn)量的關(guān)鍵時期[30],灌漿期干旱脅迫對小麥產(chǎn)生的影響遠大于其他時期[31-32]。而在實際生產(chǎn)中,我國北方的小麥常在灌漿期遇到不同程度的干旱[33]。Qin等[1]研究結(jié)果表明,在虧缺灌溉處理下,灌漿期前無人機熱成像評價小麥抗旱性的效果較差。因此,本試驗選擇在灌漿期采集人工合成小麥田間圖像數(shù)據(jù)并進行分析。作物抗旱性遺傳改良的重要目標(biāo)是提高干旱脅迫條件下的產(chǎn)量,產(chǎn)量是評價作物抗旱性的直接指標(biāo)[34]。但作物產(chǎn)量只可在成熟期獲取,僅使用產(chǎn)量性狀無法在灌漿期對小麥的抗旱性進行評價。干旱脅迫條件下,作物與產(chǎn)量相關(guān)的生理生化指標(biāo),或者生長發(fā)育指標(biāo)也可用于評價其抗旱性。然而,這些指標(biāo)的采集往往費時費力,對環(huán)境條件和性狀采集時間要求嚴(yán)格,甚至是苛刻,難以滿足對大批量材料進行抗旱性鑒定的需求。近年來,隨著無人機的普及與圖像采集設(shè)備的不斷升級,高通量表型觀測平臺逐步完善,目前,利用無人機搭載的多光譜圖像采集設(shè)備已能夠?qū)Υ筇镒魑锏亩鄻有怨庾V指數(shù)進行快速、準(zhǔn)確、高效采集。Candiago等[35]使用無人機搭載的多光譜成像儀獲取了NDVI、GNDVI和SAVI等植被指數(shù),對作物的植被活力進行評價;Shao等[36]通過相同的方法獲取了GNDVI等6種植被指數(shù),用于估算玉米的Kc(crop coefficient,Kc)值。干旱脅迫影響作物的生理生化特征,這些特征進而影響作物對光的吸收、投射和反射[37-38]。光譜指數(shù)的差異能夠間接反映作物在干旱脅迫下的生理狀態(tài),因此,借助光譜指數(shù)高效評價小麥抗旱性的策略備受關(guān)注,如李升東等[38]研究了NDVI與小麥產(chǎn)量及抗旱性的關(guān)系,證明抽穗期小麥的NDVI值與小麥抗旱性呈正相關(guān)關(guān)系。本研究所使用的無人機平臺同時搭載熱紅外及多光譜成像儀,可同時獲取多源影像。植被指標(biāo)與植株的生理狀態(tài)及形態(tài)性狀密切相關(guān),本研究在小麥灌漿期選擇植株形態(tài)結(jié)構(gòu)、葉綠素含量和冠層溫度等與植株生理狀態(tài)相關(guān)的19種光譜指數(shù),從水分狀態(tài)、生長狀態(tài)、株型結(jié)構(gòu)等層次對人工合成小麥的抗旱性進行評價。通過結(jié)合相關(guān)性分析和主成分分析發(fā)現(xiàn)OSAVI-等多種光譜指數(shù)的抗旱系數(shù)與基于產(chǎn)量的抗旱指數(shù)顯著正相關(guān),可用于輔助評價小麥抗旱性,從而提高抗旱種質(zhì)鑒定效率。

    然而,作物抗旱性是多種性狀的綜合體現(xiàn),僅利用單一性狀或生理指標(biāo)難以準(zhǔn)確評價作物抗旱性,聚合多項主要指標(biāo)進行綜合評價已逐漸成為作物抗旱性鑒定的重要理念。例如,李春情等[39]通過灰色關(guān)聯(lián)度分析篩選出超氧化物歧化酶、葉綠素含量等6種與玉米抗旱性相關(guān)的指標(biāo),通過計算加權(quán)抗旱系數(shù)對42個玉米品種的抗旱性進行綜合評價,提高了鑒定結(jié)果的可靠性;王倩等[40]比較分析3種抗旱性評價方法,發(fā)現(xiàn)利用抗旱性度量值()的綜合評價方法最適宜對成株期糜子進行抗旱性評價。李龍等[29]研究也發(fā)現(xiàn)采用綜合抗旱性度量值()有利于區(qū)分干旱對不同小麥種質(zhì)產(chǎn)量的影響力。本研究通過主成分分析將19種光譜指數(shù)的抗旱系數(shù)轉(zhuǎn)化為3個相互獨立的綜合指標(biāo),去除其中的冗余信息,進一步通過加權(quán)隸屬函數(shù)法聚合綜合指標(biāo)獲得綜合抗旱性度量值。與單指標(biāo)評價指標(biāo)相比,多指標(biāo)綜合評價結(jié)果考慮了各光譜指數(shù)的重要性及其互作關(guān)系,更具客觀性和可靠性,有利于快速準(zhǔn)確地篩選優(yōu)異等位基因種類多、育種應(yīng)用潛力大的種質(zhì)資源。

    3.2 人工合成小麥抗旱種質(zhì)資源發(fā)掘

    種質(zhì)資源的創(chuàng)新與品種改良至關(guān)重要[41]。長期追求高產(chǎn)的定向選擇和骨干親本重復(fù)利用導(dǎo)致小麥栽培品種中的抗旱遺傳基礎(chǔ)收窄[42],制約了抗旱高產(chǎn)小麥新品種的培育,因此,亟需從更加廣泛的種質(zhì)資源中挖掘新的優(yōu)異抗旱種質(zhì)以拓寬抗旱育種的遺傳基礎(chǔ)[43]。本研究所分析的人工合成小麥(AABBDD)是利用野生或栽培二粒小麥(AABB)與粗山羊草(DD)雜交創(chuàng)制的新型小麥種質(zhì),攜帶了小麥野生近緣種特有的耐逆性,同時還能與普通小麥(AABBDD)正常雜交,是拓展小麥抗逆遺傳基礎(chǔ)的優(yōu)選基因庫,具有極大的優(yōu)良基因遺傳潛力[44]。張澤全等[45]采用性狀相對值和隸屬函數(shù)法對大量人工合成小麥進行耐低溫綜合評價,篩選出了一批低溫耐性強的人工合成小麥種質(zhì)并提供育種利用。王秋云等[46]對56份人工合成小麥進行白粉病抗性鑒定,發(fā)現(xiàn)80%以上的供試材料對白粉病具有良好的抗性。然而,目前對于人工合成小麥抗旱種質(zhì)資源的發(fā)掘工作相對薄弱,在本研究中,人工合成小麥也存在產(chǎn)量較低、性狀不穩(wěn)定的缺點。與2022年相比,2023年人工合成小麥產(chǎn)量降低,這可能與2023年開花期有降水相關(guān)。本研究結(jié)合不同評價方法發(fā)掘到多個強抗旱和抗旱人工合成小麥,將為小麥抗旱高產(chǎn)優(yōu)異親本資源圃注入“新鮮血液”,是種質(zhì)資源發(fā)掘工作的重要進展。

    4 結(jié)論

    利用無人機搭載多功能成像儀采集多源影像是高效獲取小麥抗旱表型信息的有效方法。NDVI、OSAVI、CIre、NDRE和綜合抗旱性度量值可用于輔助鑒定小麥種質(zhì)抗旱性。結(jié)合單指標(biāo)及綜合指標(biāo)評價方法鑒定出8份抗旱人工合成小麥種質(zhì)(SW013、SW017、SW040、SW042、SW68、SW100、SW106和SW123)和2份強抗旱人工合成小麥種質(zhì)(SW004、SW009)。

    [1] QIN W L, WANG J, MA L F, WANG F L, HU N Y, YANG X Y, XIAO Y Y, ZHANG Y H, SUN Z C, WANG Z M, YU K. UAV-based multi-temporal thermal imaging to evaluate wheat drought resistance in different deficit irrigation regimes. Remote Sensing, 2022, 14(21): 5608.

    [2] BAPELA T, SHIMELIS H, TSILO T J, MATHEW I. Genetic improvement of wheat for drought tolerance: Progress, challenges and opportunities. Plants, 2022, 11(10): 1331.

    [3] 廖祥政, 王瑾, 周榮華, 任正隆, 賈繼增. 發(fā)掘人工合成小麥中千粒重QTL的有利等位基因. 作物學(xué)報, 2008, 34(11): 1877-1884.

    LIAO X Z, WANG J, ZHOU R H, REN Z L, JIA J Z. Mining favorable alleles of QTLs conferring 1000-grain weight from synthetic wheat. Acta Agronomica Sinica, 2008, 34(11): 1877-1884. (in Chinese)

    [4] 王曉鳴, 邱麗娟, 景蕊蓮, 任貴興, 李英慧, 李春輝, 秦培友, 谷勇哲, 李龍. 作物種質(zhì)資源表型性狀鑒定評價: 現(xiàn)狀與趨勢. 植物遺傳資源學(xué)報, 2022, 23(1): 12-20.

    WANG X M, QIU L J, JING R L, REN G X, LI Y H, LI C H, QIN P Y, GU Y Z, LI L. Evaluation on phenotypic traits of crop germplasm: status and development. Journal of Plant Genetic Resources, 2022, 23(1): 12-20. (in Chinese)

    [5] 李龍, 毛新國, 王景一, 李超男, 柳玉平, 景蕊蓮. 小麥抗旱性研究進展與展望. 干旱地區(qū)農(nóng)業(yè)研究, 2023, 41(3): 11-20

    LI L, MAO X G, WANG J Y, LI C N, LIU Y P, JING R L. Progress and prospect of wheat research on drought resistance. Agricultural Research in the Arid Areas, 2023, 41(3): 11-20. (in Chinese)

    [6] 白志英, 李存東, 孫紅春, 趙金鋒. 小麥代換系抗旱生理指標(biāo)的主成分分析及綜合評價. 中國農(nóng)業(yè)科學(xué), 2008, 41(12): 4264-4272.doi: 10.3864/j.issn.0578-1752.2008.12.043.

    Bai Z Y, Li C D, Sun H C, Zhao J F. Principal component analysis and comprehensive evaluation on physiological indices of drought resistance in wheat substitution lines. Scientia Agricultura Sinica, 2008, 41(12): 4264-4272. doi: 10.3864/j.issn.0578-1752.2008.12.043. (in Chinese)

    [7] 王士強, 胡銀崗, 佘奎軍, 周琳璘, 孟凡磊. 小麥抗旱相關(guān)農(nóng)藝性狀和生理生化性狀的灰色關(guān)聯(lián)度分析. 中國農(nóng)業(yè)科學(xué), 2007, 40(11): 2452-2459. doi: 10.3864/j.issn.0578-1752.at-2006-8460.

    Wang S Q, Hu Y G, She K J, Zhou L L, Meng F L. Gray relational grade analysis of agronomical and physi-biochemical traits related to drought tolerance in wheat. Scientia Agricultura Sinica, 2007, 40(11): 2452-2459. doi: 10.3864/j.issn.0578-1752.at-2006-8460. (in Chinese)

    [8] 張軍, 魏國, 彭彥珉, 汪洪濤, 黃修利, 吳金芝, 黃明, 李友軍. 8份強筋小麥品種抗旱性評價. 麥類作物學(xué)報, 2023. https://link. cnki.net/urlid/61.1359.S.20231107.0924.006.

    ZHANG J, WEI G, PENG Y M, WANG H T, HUANG X L, WU JI Z, HUANG M, LI Y J. Drought resistance evaluation of eight strong gluten wheat varieties. Journal of Triticeae Crops, 2023. https://link. cnki.net/urlid/61.1359.S.20231107.0924.006. (in Chinese)

    [9] 馬玉慧, 張小虎, 馬小樂, 司二靜, 孟亞雄, 姚立蓉, 王化俊, 汪軍成, 李葆春. 干旱脅迫下春小麥品種(系)萌發(fā)期抗旱性鑒定與評價. 分子植物育種, 2022, 20(19): 6459-6473.

    MA Y H, ZHANG X H, MA X L, SI E J, MENG Y X, YAO L R, WANG H J, WANG J C, LI B C. Identification and evaluation of drought resistance of spring wheat varieties (lines) at germination under drought stress. Molecular Plant Breeding, 2022, 20(19): 6459-6473. (in Chinese)

    [10] 孟雨, 溫鵬飛, 丁志強, 田文仲, 張學(xué)品, 賀利, 段劍釗, 劉萬代, 馮偉. 基于熱紅外圖像的小麥品種抗旱性鑒定與評價. 中國農(nóng)業(yè)科學(xué), 2022, 55(13): 2538-2551. doi: 10.3864/j.issn.0578-1752.2022. 13.005.

    MENG Y, WEN P F, DING Z Q, TIAN W Z, ZHANG X P, HE L, DUAN J Z, LIU W D, FENG W. identification and evaluation of drought resistance of wheat varieties based on thermal infrared image. Scientia Agricultura Sinica, 2022, 55(13): 2538-2551. doi: 10.3864/ j.issn.0578-1752.2022.13.005. (in Chinese)

    [11] DAS S, CHRISTOPHER J, CHOUDHURY M R, APAN A, CHAPMAN S W, MENZIES N P, DANG Y P. Evaluation of drought tolerance of wheat genotypes in rain-fed sodic soil environments using high-resolution UAV remote sensing techniques. Biosystems Engineering, 2022, 217: 68-82.

    [12] WANG F M, HUANG J F, TANG Y L, WANG X Z. New vegetation index and its application in estimating leaf area index of rice. Rice Science, 2007, 14(3): 195-203.

    [13] GITELSON A A, KEYDAN G P, MERZLYAK M N. Three-band model for noninvasive estimation of chlorophyll, carotenoids, and anthocyanin contents in higher plant leaves. Geophysical Research Letters, 2006, 33(11): 11402.

    [14] Richardson A, Wiegand C. Distinguishing vegetation from soil background information. Photogrammetric Engineering and Remote Sensing, 1977, 43(12): 1541-1552.

    [15] GILLESPIE A R, KAHLE A B, WALKER R E. Color enhancement of highly correlated images: Ⅱchannel ratio and “chromaticity” transformation techniques. Remote Sensing of Environment, 1987, 22(3): 343-365.

    [16] GITELSON A A, KAUFMAN Y J, MERZLYAK M N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sensing of Environment, 1996, 58(3): 289-298.

    [17] DAUGHTRY C S T, WALTHALL C L, KIM M S, DE COLSTOUN E B, MCMURTREY J E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sensing of Environment, 2000, 74(2): 229-239.

    [18] JAY S, GORRETTA N, MOREL J, MAUPAS F, BENDOULA R, RABATEL G, DUTARTRE D, COMAR A, BARET F. Estimating leaf chlorophyll content in sugar beet canopies using millimeter to centimeter-scale reflectance imagery. Remote Sensing of Environment, 2017, 198: 173-186.

    [19] DASH J, CURRAN P J. Theterrestrial chlorophyll index. International Journal of Remote Sensing, 2004, 25(23): 5403-5413.

    [20] Fitzgerald G J, Rodriguez D, Christensen L K, Belford R, Sadras V O, Clarke T R. Spectral and thermal sensing for nitrogen and water status in rainfed and irrigated wheat environments. Precision agriculture, 2006, 7(4): 233-248.

    [21] ROUSE J W J, HAAS R H, SCHELL J A, DEERING D W. Monitoring vegetation systems in the great Plains with ERTS. National Aeronautics and Space Administration Special Publication, 1974, 351(1): 309.

    [22] XU F X, GAO Z Q, JIANG X P, SHANG W T, NING J C, SONG D B, AI J Q. A UAV and S2A data-based estimation of the initial biomass of green algae in the South Yellow Sea. Marine pollution bulletin, 2018, 128: 408-414.

    [23] Pe?uelas J, Gamon J A, Fredeen A L, Merino J, Field C B. Reflectance indices associated with physiological changes in nitrogen- and water-limited sunflower leaves. Remote Sensing of Environment, 1994, 48(2): 135-146.

    [24] Elsayed S, Rischbeck P, Schmidhalter U. Comparing the performance of active and passive reflectance sensors to assess the normalized relative canopy temperature and grain yield of drought- stressed barley cultivars. Field Crops Research, 2015, 177: 148-160.

    [25] RONDEAUX G, STEVEN M, BARET F. Optimization of soil-adjusted vegetation indices. Remote Sensing of Environment, 1996, 55(2): 95-107.

    [26] PENUELAS J, BARET F, FILELLA I. Semi-empirical indices to assess carotenoids/chlorophyll a ratio from leaf spectral reflectance. Photosynthetica, 1995, 31(2): 221-230.

    [27] Jordan C F. Derivation of leaf‐area index from quality of light on the forest floor. Ecology, 1969, 50(4): 663-666.

    [28] BROGE N, LEBLANC E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sensing of Environment, 2001, 76(2): 156-172.

    [29] 李龍, 毛新國, 王景一, 昌小平, 柳玉平, 景蕊蓮. 小麥種質(zhì)資源抗旱性鑒定評價. 作物學(xué)報, 2018, 44(7): 988-999.

    LI L, MAO X G, WANG J Y, CHANG X P, LIU Y P, JING R L. Drought tolerance evaluation of wheat germplasm resources. Acta Agronomica Sinica, 2018, 44(7): 988-999. (in Chinese)

    [30] 胡芬, 趙聚寶, 姜雁北, 梅旭榮, 鐘兆站. 小麥灌漿期土壤干旱對子粒發(fā)育影響的研究. 中國農(nóng)業(yè)氣象, 1993, 14(4): 8-11.

    HU F, ZHAO J B, JIANG Y B, MEI X R, ZHONG Z Z. Study on the effects of soil drought on grain development during wheat filling stage. Chinese Journal of Agrometeorology, 1993, 14(4): 8-11. (in Chinese)

    [31] YANG J C, ZHANG J H. Grain filling of cereals under soil drying. The New Phytologist, 2006, 169(2): 223-236.

    [32] FAROOQ M, HUSSAIN M, SIDDIQUE K H M. Drought stress in wheat during flowering and grain-filling periods. Critical Reviews in Plant Sciences, 2014, 33(4): 331-349.

    [33] 呂金印, 山侖, 高俊鳳, 覃鳳云, 楊淑慎. 干旱對小麥灌漿期旗葉光合等生理特性的影響. 干旱地區(qū)農(nóng)業(yè)研究, 2003, 21(2): 77-81.

    Lü J Y, SHAN L, GAO J F, QIN F Y, YANG S S. Effect of drought stress on photosynthesis and some other physiological characteristics in flag leaf during grain filling of wheat. Agricultural Research in the Arid Areas, 2003, 21(2): 77-81. (in Chinese)

    [34] 王淑榮, 費德友, 唐宗祥. 干旱脅迫對四川小麥品種產(chǎn)量及農(nóng)藝性狀的影響. 農(nóng)業(yè)與技術(shù), 2020, 40(13): 14-17.

    WANG S R, FEI D Y, TANG Z X. Effects of drought stress on yield and agronomic traits of wheat varieties in Sichuan. Agriculture and Technology, 2020, 40(13): 14-17. (in Chinese)

    [35] CANDIAGO S, REMONDINO F, DE GIGLIO M, DUBBINI M, GATTELLI M. Evaluating multispectral images and vegetation indices for precision farming applications from UAV images. Remote Sensing, 2015, 7(4): 4026-4047.

    [36] SHAO G M, HAN W T, ZHANG H H, LIU S Y, WANG Y, ZHANG L Y, CUI X. Mapping maize crop coefficient Kc using random forest algorithm based on leaf area index and UAV-based multispectral vegetation indices. Agricultural Water Management, 2021, 252: 106906.

    [37] 李映雪, 朱艷, 田永超, 姚霞, 秦曉東, 曹衛(wèi)星. 小麥葉片氮含量與冠層反射光譜指數(shù)的定量關(guān)系. 作物學(xué)報, 2006, 32(3): 358-362.

    LI Y X, ZHU Y, TIAN Y C, YAO X, QIN X D, CAO W X. Quantitative relationship between leaf nitrogen concentration and canopy reflectance spectra. Acta Agronomica Sinica, 2006, 32(3): 358-362. (in Chinese)

    [38] 李升東, 王法宏, 司紀(jì)升, 孔令安, 馮波, 張賓. 不同基因型小麥NDVI值與產(chǎn)量的關(guān)系. 干旱地區(qū)農(nóng)業(yè)研究, 2008, 26(6): 47-50.

    LI S D, WANG F H, SI J S, KONG L A, FENG B, ZHANG B. The relationship between NDVI values and yield of different genotype wheat varieties. Agricultural Research in the Arid Areas, 2008, 26(6): 47-50. (in Chinese)

    [39] 李春情, 劉翔宇, 閆鵬, 周六千, 盧霖, 董志強, 徐江. 不同玉米品種抗旱性的生理鑒定與綜合評價. 作物雜志, 2023. https://link. cnki.net/urlid/11.1808.S.20230926.1429.006.

    LI C Q, LIU X Y, YAN P, ZHOU L Q, LU L, DONG Z Q, XU J. Physiological identification and comprehensive evaluation of drought resistance of different maize varieties. Crop Journal, 2023. https://link. cnki.net/urlid/11.1808.S.20230926.1429.006. (in Chinese)

    [40] 王倩, 董孔軍, 薛亞鵬, 劉少雄, 王若楠, 楊佳琪, 陸平, 王瑞云, 楊天育, 劉敏軒. 糜子核心種質(zhì)成熟期抗旱性鑒定評價與抗旱種質(zhì)篩選. 中國農(nóng)業(yè)科學(xué), 2023, 56(21): 4163-4174.doi:10.3864/j.issn. 0578-1752.2023.21.003.

    Wang Q, Dong K J, Xue Y P, Liu S X, Wang R N, Yang J Q, Lu P, Wang R Y, Yang T Y, Liu M X. Drought resistance identification and evaluation of mature stage in foxtail millet core germplasm and selection of drought-resistant germplasm. Scientia Agricultura Sinica, 2023, 56(21): 4163-4174.doi:10.3864/j.issn.0578- 1752.2023.21.003. (in Chinese)

    [41] 張青顯. 野生遺傳基因種質(zhì)資源在小麥育種中的應(yīng)用及遺傳特征. 農(nóng)業(yè)科技通訊, 2022(7): 183-184, 188.

    ZHANG Q X. Application and genetic characteristics of wild genetic gene germplasm resources in wheat breeding. Bulletin of Agricultural Science and Technology, 2022(7): 183-184, 188. (in Chinese)

    [42] EVENSON R E, GOLLIN D. Assessing the impact of the green revolution, 1960 to 2000. Science, 2003, 300(5620): 758-762.

    [43] DRECCER M F, BORGOGNONE M G, OGBONNAYA F C, TRETHOWAN R M, Winter B. CIMMYT-selected derived synthetic bread wheats for rainfed environments: Yield evaluation in Mexico and Australia. Field Crops Research, 2007, 100(2/3): 218-228.

    [44] RATTEY A, SHORTER R. Evaluation of CIMMYT conventional and synthetic spring wheat germplasm in rainfed sub-tropical environments: I. Grain yield. Field Crops Research, 2010, 118(3): 273-281.

    [45] 張澤全, 舒長生, 賴運平, 董雪芳, 李俊, 魏會廷, 彭正松, 楊武云. 人工合成小麥幼苗耐低溫特性鑒定. 中國農(nóng)學(xué)通報, 2010, 26(12): 135-138.

    ZHANG Z Q, SHU C S, LAI Y P, DONG X F, LI J, WEI H T, PENG Z S, YANG W Y. Evaluation of low-temperature tolerant germplasm resource on seedling traits of synthetic hexaploid wheats. Chinese Agricultural Science Bulletin, 2010, 26(12): 135-138. (in Chinese)

    [46] 王秋云, 馮雯杰, 張巧玲, 趙凱, 張憲軍, 劉淑珍, 鮑印廣. 56份人工合成六倍體小麥的遺傳多樣性分析. 麥類作物學(xué)報, 2021, 41(8): 960-968.

    WANG Q Y, FENG W J, ZHANG Q L, ZHAO K, ZHANG X J, LIU S Z, BAO Y G. Genetic diversity analysis of 56 synthetic hexaploid wheats. Journal of Triticeae Crops, 2021, 41(8): 960-968. (in Chinese)

    Drought Resistance Evaluation of Synthetic Wheat at Grain Filling Using UAV-based Multi-source Imagery Data

    YAN Wen1,2, JIN XiuLiang2, LI Long2, XU ZiHan3, SU Yue1,2, ZHANG YueQiang3, JING RuiLian2, MAO XinGuo2, SUN DaiZhen1

    1College of Agronomy, Shanxi Agricultural University, Taigu 030801, Shanxi;2Institute of Crop Sciences, Chinese Academy of Agricultural Sciences/State Key Laboratory of Crop Gene Resources and Breeding, Beijing 100081;3Institute of Nuclear and Biological Technologies, Xinjiang Academy of Agricultural Sciences/Key Laboratory of Crop Ecophysiology and Farming System in Desert Oasis Region, Ministry of Agriculture and Rural Affairs/Xinjiang Key Laboratory of Crop Biotechnology, Urumqi 830091

    【Objective】To evaluate the drought resistance of synthetic wheat germplasm based on multi-source images collected by unmanned aerial vehicle (UAV) and yield data, explore high-throughput indices for drought resistance evaluation, and identify synthetic wheat germplasm resources with drought resistance.This provides technical support and germplasm materials for accelerating the expansion of drought-resistant genetic resources for wheat and enhancing the level of breeding for dryland wheat.【Method】Eighty synthetic wheat germplasm and the control variety Xin Chun 37 were used as plant materials, which were sown in the field and treated with a water regime of drought stress and irrigation. Multi-source images of test materials during filling stage were collected by multi-spectral and thermal infrared cameras equipped with unmanned aerial vehicle, and the spectral index of each test material was extracted by threshold segmentation. The analyses of Pearson’s correlation and principal component were performed to identify drought resistance-related spectral indices, and the drought resistance of each synthetic wheat germplasm was determined by single index and comprehensive evaluation methods. 【Result】The drought resistance coefficients of 19 spectral indices of 80 synthetic wheat germplasm were calculated based on multisource imagery data obtained from unmanned aerial vehicles. The correlation analysis between the spectral indices and the yield-based drought index () showed that among the drought resistance coefficients of the 19 spectral indices, OSAVI exhibited the highest correlation with the drought index, while NDVI, CIre, and NDRE demonstrated relatively strong associations with the drought index. However, the different drought indices showed a high correlation, resulting in redundant information. The drought resistance coefficients of the 19 spectral indices were transformed into three independent comprehensive indicators through principal component analysis, with contribution rates of 59.6%, 12.0% and 9.6%, respectively. The comprehensive drought resistance index () for each synthetic wheat germplasm were calculated by aggregating the three independent comprehensive indicators using the weighted membership function method. 6 and 5 synthetic wheat germplasms with strong drought resistance were identified based onand, respectively. Among them, 2 germplasms (SW004 and SW009) with high drought resistance were detected based on bothand. Furthermore, the drought resistance of the 80 synthetic wheat germplasms was graded based on the drought resistance coefficient of OSAVI, and the grading results were found to be consistent with that based on thevalue. Among the six strongly drought-resistant germplasms identified based on the drought resistance coefficient of OSAVI, five of them were also classified as strongly drought-resistant germplasms based on comprehensive drought resistance evaluation.【Conclusion】The spectral indices NDVI, OSAVI, CIreand NDRE extracted from UAV-based multi-source images, as well as the drought resistance comprehensive evaluation value can be used to assist in the identification of drought resistance of wheat germplasm.

    multi-source images; spectral index; synthetic wheat; drought resistance; grain filling

    10.3864/j.issn.0578-1752.2024.09.005

    2024-01-05;

    2024-03-11

    國家重點研發(fā)計劃(2023YFD1201003,2022YFD1200201)、國家小麥產(chǎn)業(yè)技術(shù)體系(CARS-03-5)

    燕雯,E-mail:yanwenzzl@163.com。通信作者孫黛珍,E-mail:sdz64@126.com。通信作者毛新國,E-mail:maoxinguo@caas.cn

    (責(zé)任編輯 李莉)

    猜你喜歡
    抗旱性抗旱種質(zhì)
    美國聯(lián)邦政府撥款8 400萬美元用于西部地區(qū)抗旱
    華南地區(qū)最大農(nóng)作物種質(zhì)資源保護庫建成
    種子穿新“外衣”鎖水抗旱長得好
    果園抗旱“24字方針”
    亞麻抗白粉病種質(zhì)資源的鑒定與篩選
    防汛抗旱
    不同光質(zhì)對黃瓜幼苗抗旱性的影響
    貴州玉米種質(zhì)資源遺傳多樣性及核心種質(zhì)庫構(gòu)建
    基于稱重法的煙草伸根期抗旱性鑒定
    紅錐種質(zhì)早期生長表現(xiàn)
    大片电影免费在线观看免费| 高清不卡的av网站| 免费黄频网站在线观看国产| 最近最新中文字幕免费大全7| 人妻一区二区av| 国产亚洲一区二区精品| 黄片播放在线免费| 91aial.com中文字幕在线观看| 中文字幕人妻丝袜一区二区 | 18禁观看日本| 丁香六月欧美| 国产97色在线日韩免费| 久久人妻熟女aⅴ| 成人午夜精彩视频在线观看| 日韩一区二区三区影片| 啦啦啦在线免费观看视频4| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久人妻精品电影 | 你懂的网址亚洲精品在线观看| a级片在线免费高清观看视频| 免费黄色在线免费观看| 日韩av免费高清视频| 欧美日韩视频精品一区| 飞空精品影院首页| 丝袜美足系列| 亚洲欧洲国产日韩| 这个男人来自地球电影免费观看 | 精品一区二区免费观看| 亚洲精品久久久久久婷婷小说| 人人妻人人爽人人添夜夜欢视频| 美女视频免费永久观看网站| 午夜激情av网站| 国精品久久久久久国模美| 国产精品国产三级专区第一集| 国产免费又黄又爽又色| 午夜福利影视在线免费观看| 国产亚洲午夜精品一区二区久久| 日韩一本色道免费dvd| 天天躁狠狠躁夜夜躁狠狠躁| 色婷婷久久久亚洲欧美| 免费av中文字幕在线| 久久久精品国产亚洲av高清涩受| 国产精品偷伦视频观看了| 美女脱内裤让男人舔精品视频| 另类亚洲欧美激情| 亚洲精品中文字幕在线视频| 久久综合国产亚洲精品| xxx大片免费视频| 国产极品粉嫩免费观看在线| 欧美激情 高清一区二区三区| 久久久久久免费高清国产稀缺| av电影中文网址| 欧美 亚洲 国产 日韩一| 少妇被粗大猛烈的视频| 大话2 男鬼变身卡| 性高湖久久久久久久久免费观看| 成年动漫av网址| 欧美 日韩 精品 国产| 午夜免费鲁丝| 国产精品亚洲av一区麻豆 | 国产精品女同一区二区软件| 菩萨蛮人人尽说江南好唐韦庄| 一级,二级,三级黄色视频| 女人久久www免费人成看片| 国产一区有黄有色的免费视频| 国产亚洲av高清不卡| 制服人妻中文乱码| 操出白浆在线播放| 人成视频在线观看免费观看| 香蕉丝袜av| 欧美激情高清一区二区三区 | 男女边摸边吃奶| 久久久国产精品麻豆| 99国产精品免费福利视频| 亚洲精品在线美女| 久久精品久久精品一区二区三区| 亚洲欧洲国产日韩| 欧美人与性动交α欧美精品济南到| 香蕉丝袜av| 国产精品国产av在线观看| 国产精品99久久99久久久不卡 | 亚洲精品久久成人aⅴ小说| 婷婷成人精品国产| 妹子高潮喷水视频| 啦啦啦啦在线视频资源| 久久韩国三级中文字幕| 中国国产av一级| 爱豆传媒免费全集在线观看| 美女主播在线视频| 欧美亚洲 丝袜 人妻 在线| 国产精品 国内视频| 亚洲成色77777| 久久久久久人妻| 最近2019中文字幕mv第一页| 国产精品熟女久久久久浪| 日韩中文字幕视频在线看片| 国产精品国产三级专区第一集| 男女之事视频高清在线观看 | av国产精品久久久久影院| av一本久久久久| 各种免费的搞黄视频| 色婷婷av一区二区三区视频| 中文欧美无线码| 99国产综合亚洲精品| 欧美日本中文国产一区发布| 日本91视频免费播放| 国产精品久久久久久久久免| 婷婷色av中文字幕| 国产欧美日韩一区二区三区在线| 青春草国产在线视频| 女人被躁到高潮嗷嗷叫费观| 精品视频人人做人人爽| 国产乱人偷精品视频| 黑人猛操日本美女一级片| 黄色毛片三级朝国网站| 狂野欧美激情性bbbbbb| 久久久国产欧美日韩av| 国产在线免费精品| 十八禁人妻一区二区| 大码成人一级视频| 国产成人a∨麻豆精品| 在线观看www视频免费| 精品一品国产午夜福利视频| 精品国产一区二区三区久久久樱花| 欧美精品一区二区免费开放| 美女扒开内裤让男人捅视频| 咕卡用的链子| 国产爽快片一区二区三区| 人妻一区二区av| 天天躁夜夜躁狠狠久久av| e午夜精品久久久久久久| 亚洲av福利一区| 国产日韩欧美亚洲二区| 久久久久久久大尺度免费视频| 亚洲欧美日韩另类电影网站| 老司机在亚洲福利影院| 中文字幕高清在线视频| 精品一区二区三区四区五区乱码 | 最新在线观看一区二区三区 | 男人舔女人的私密视频| 久久久久久久大尺度免费视频| 999久久久国产精品视频| 9热在线视频观看99| 亚洲婷婷狠狠爱综合网| 精品一区二区三区av网在线观看 | 国产xxxxx性猛交| 丁香六月天网| 97在线人人人人妻| 丁香六月天网| 伊人久久国产一区二区| 美女视频免费永久观看网站| 天天操日日干夜夜撸| 亚洲国产欧美一区二区综合| bbb黄色大片| 校园人妻丝袜中文字幕| 毛片一级片免费看久久久久| 高清欧美精品videossex| 大片免费播放器 马上看| 一区二区三区乱码不卡18| 777久久人妻少妇嫩草av网站| 一边摸一边做爽爽视频免费| 悠悠久久av| 亚洲成人一二三区av| 国产成人欧美| 晚上一个人看的免费电影| 热re99久久精品国产66热6| 伊人久久大香线蕉亚洲五| 老熟女久久久| 十八禁高潮呻吟视频| 亚洲图色成人| 国产日韩欧美亚洲二区| www.自偷自拍.com| 男人添女人高潮全过程视频| 亚洲人成网站在线观看播放| 日韩,欧美,国产一区二区三区| 宅男免费午夜| 久久久精品免费免费高清| 午夜福利乱码中文字幕| 亚洲国产毛片av蜜桃av| 国产午夜精品一二区理论片| 亚洲精品一区蜜桃| 免费黄网站久久成人精品| 欧美日韩亚洲高清精品| 999精品在线视频| 国产成人欧美| 丝瓜视频免费看黄片| 亚洲国产毛片av蜜桃av| 日韩欧美一区视频在线观看| 久久亚洲国产成人精品v| 尾随美女入室| 尾随美女入室| 90打野战视频偷拍视频| 亚洲中文av在线| av有码第一页| 国产熟女欧美一区二区| 久久人妻熟女aⅴ| 一级片免费观看大全| 在线观看免费日韩欧美大片| 一区在线观看完整版| 看非洲黑人一级黄片| 极品人妻少妇av视频| 成人国语在线视频| 亚洲国产中文字幕在线视频| 欧美少妇被猛烈插入视频| 久久国产亚洲av麻豆专区| 亚洲国产精品一区三区| 国产一区二区 视频在线| 久久久久人妻精品一区果冻| 久久精品久久精品一区二区三区| 久久青草综合色| 久久精品亚洲熟妇少妇任你| 亚洲精品成人av观看孕妇| 亚洲精品第二区| 侵犯人妻中文字幕一二三四区| 国产片特级美女逼逼视频| 飞空精品影院首页| 国产淫语在线视频| 欧美老熟妇乱子伦牲交| 国产探花极品一区二区| 狠狠婷婷综合久久久久久88av| 国产黄色免费在线视频| 欧美人与性动交α欧美精品济南到| 国产福利在线免费观看视频| 欧美激情高清一区二区三区 | 亚洲国产欧美在线一区| 久久ye,这里只有精品| 免费观看性生交大片5| 亚洲精品久久午夜乱码| av国产精品久久久久影院| 国产精品嫩草影院av在线观看| 日韩一区二区三区影片| 丰满饥渴人妻一区二区三| videos熟女内射| 中文字幕人妻丝袜一区二区 | 丰满迷人的少妇在线观看| 午夜福利视频精品| 亚洲精品日本国产第一区| 日韩一区二区三区影片| 18禁动态无遮挡网站| 少妇精品久久久久久久| 国产又色又爽无遮挡免| 国产淫语在线视频| 赤兔流量卡办理| 在线天堂中文资源库| 欧美日韩一区二区视频在线观看视频在线| 高清黄色对白视频在线免费看| 天天操日日干夜夜撸| 99精品久久久久人妻精品| 国产成人欧美在线观看 | 亚洲成人手机| 侵犯人妻中文字幕一二三四区| 亚洲精品久久午夜乱码| 咕卡用的链子| 777久久人妻少妇嫩草av网站| 亚洲精品在线美女| 一个人免费看片子| 久久精品aⅴ一区二区三区四区| 欧美黑人欧美精品刺激| 欧美日韩亚洲综合一区二区三区_| 久久免费观看电影| 久久久久精品久久久久真实原创| 成年女人毛片免费观看观看9 | 各种免费的搞黄视频| 青春草亚洲视频在线观看| 久久精品亚洲熟妇少妇任你| 蜜桃国产av成人99| 在线观看一区二区三区激情| 老汉色av国产亚洲站长工具| 18在线观看网站| 男人操女人黄网站| 日韩一卡2卡3卡4卡2021年| 久久久国产一区二区| 亚洲情色 制服丝袜| 午夜激情av网站| av在线老鸭窝| 黄色一级大片看看| 国产成人精品久久二区二区91 | 可以免费在线观看a视频的电影网站 | 久久久精品94久久精品| 成年女人毛片免费观看观看9 | 午夜久久久在线观看| 亚洲欧洲精品一区二区精品久久久 | 欧美精品一区二区免费开放| 十分钟在线观看高清视频www| 国产人伦9x9x在线观看| 最近手机中文字幕大全| 韩国精品一区二区三区| 精品人妻熟女毛片av久久网站| 女人久久www免费人成看片| 只有这里有精品99| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人av在线免费| 夫妻性生交免费视频一级片| 黄频高清免费视频| 亚洲七黄色美女视频| 国产成人欧美在线观看 | 亚洲五月色婷婷综合| 少妇 在线观看| 精品人妻一区二区三区麻豆| 国产一卡二卡三卡精品 | 免费高清在线观看视频在线观看| 成人国产av品久久久| 一级片免费观看大全| 又大又爽又粗| 一级毛片 在线播放| 久久久国产一区二区| 久久久久久久精品精品| 久久精品亚洲熟妇少妇任你| 亚洲专区中文字幕在线 | 国产国语露脸激情在线看| av一本久久久久| 我的亚洲天堂| 亚洲一区中文字幕在线| 老司机靠b影院| a级片在线免费高清观看视频| 男女之事视频高清在线观看 | 亚洲熟女毛片儿| 国产亚洲av片在线观看秒播厂| 国产一卡二卡三卡精品 | 亚洲av欧美aⅴ国产| 香蕉国产在线看| 欧美另类一区| 一级a爱视频在线免费观看| 久热这里只有精品99| 国产熟女欧美一区二区| 久久久亚洲精品成人影院| 99热国产这里只有精品6| 国产成人欧美在线观看 | 日韩 亚洲 欧美在线| 麻豆精品久久久久久蜜桃| 午夜福利乱码中文字幕| 一级毛片黄色毛片免费观看视频| 久久久精品国产亚洲av高清涩受| 国产精品女同一区二区软件| 麻豆乱淫一区二区| a级片在线免费高清观看视频| 亚洲av电影在线进入| 欧美日韩成人在线一区二区| 男女免费视频国产| 亚洲精品久久久久久婷婷小说| 日韩大片免费观看网站| 黄色毛片三级朝国网站| 侵犯人妻中文字幕一二三四区| 久热这里只有精品99| 亚洲欧美一区二区三区黑人| 哪个播放器可以免费观看大片| 亚洲国产av新网站| 亚洲欧洲日产国产| 久久久久精品性色| 亚洲,欧美精品.| 丝袜美腿诱惑在线| 亚洲一卡2卡3卡4卡5卡精品中文| 多毛熟女@视频| 一边摸一边抽搐一进一出视频| 亚洲欧美色中文字幕在线| 在线天堂最新版资源| 最新在线观看一区二区三区 | 岛国毛片在线播放| 欧美人与性动交α欧美精品济南到| 啦啦啦视频在线资源免费观看| 久久久精品区二区三区| 成人午夜精彩视频在线观看| 国产免费福利视频在线观看| 日本vs欧美在线观看视频| 91精品伊人久久大香线蕉| 丝袜在线中文字幕| 精品少妇黑人巨大在线播放| 国产黄频视频在线观看| 日本猛色少妇xxxxx猛交久久| 最近手机中文字幕大全| 三上悠亚av全集在线观看| www.av在线官网国产| 久久这里只有精品19| 制服人妻中文乱码| 国产精品久久久久久精品电影小说| 人人妻,人人澡人人爽秒播 | 国产成人a∨麻豆精品| 中文字幕亚洲精品专区| 日本黄色日本黄色录像| 久久天躁狠狠躁夜夜2o2o | 欧美人与性动交α欧美精品济南到| 青青草视频在线视频观看| 少妇人妻精品综合一区二区| 涩涩av久久男人的天堂| 国产成人精品在线电影| 亚洲精品,欧美精品| av天堂久久9| 韩国高清视频一区二区三区| 看免费成人av毛片| 人人妻人人澡人人看| 人人澡人人妻人| 丰满饥渴人妻一区二区三| 十八禁网站网址无遮挡| 日韩av在线免费看完整版不卡| 亚洲一级一片aⅴ在线观看| 精品久久久精品久久久| 一级爰片在线观看| 又大又爽又粗| 国产视频首页在线观看| 亚洲一区中文字幕在线| 国产淫语在线视频| 日韩av不卡免费在线播放| 久久久久精品国产欧美久久久 | a 毛片基地| 国产一卡二卡三卡精品 | 男女国产视频网站| 日韩一区二区三区影片| 国产探花极品一区二区| 欧美日韩福利视频一区二区| av国产久精品久网站免费入址| 国产亚洲精品第一综合不卡| 永久免费av网站大全| 国产视频首页在线观看| 久久人人97超碰香蕉20202| 午夜福利免费观看在线| 观看av在线不卡| 一级毛片我不卡| 日韩 欧美 亚洲 中文字幕| 国产极品天堂在线| 日本爱情动作片www.在线观看| 国产在线一区二区三区精| 满18在线观看网站| 国产日韩一区二区三区精品不卡| 国产精品蜜桃在线观看| 国产成人精品在线电影| 亚洲av成人精品一二三区| 精品久久久精品久久久| 一区福利在线观看| 国产极品天堂在线| 久久影院123| 国产精品蜜桃在线观看| 成人亚洲欧美一区二区av| xxx大片免费视频| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 纵有疾风起免费观看全集完整版| 成人18禁高潮啪啪吃奶动态图| 国产爽快片一区二区三区| 久久女婷五月综合色啪小说| 精品国产一区二区久久| 亚洲av电影在线进入| 日韩一本色道免费dvd| 日韩成人av中文字幕在线观看| 亚洲av成人不卡在线观看播放网 | av网站在线播放免费| 国产精品久久久久久精品古装| 精品国产乱码久久久久久男人| 桃花免费在线播放| 涩涩av久久男人的天堂| 亚洲第一区二区三区不卡| 久久久久精品国产欧美久久久 | 十八禁高潮呻吟视频| 一级,二级,三级黄色视频| 热re99久久国产66热| 欧美黑人精品巨大| av福利片在线| 女性被躁到高潮视频| 久久鲁丝午夜福利片| 少妇的丰满在线观看| 狂野欧美激情性bbbbbb| 19禁男女啪啪无遮挡网站| av网站在线播放免费| 久久久久久久大尺度免费视频| 51午夜福利影视在线观看| 91老司机精品| 巨乳人妻的诱惑在线观看| 亚洲天堂av无毛| 无限看片的www在线观看| 1024香蕉在线观看| 桃花免费在线播放| 丰满饥渴人妻一区二区三| 在线天堂最新版资源| 国产人伦9x9x在线观看| 亚洲图色成人| 麻豆av在线久日| 欧美乱码精品一区二区三区| 宅男免费午夜| 久久久久视频综合| 婷婷色av中文字幕| 国产福利在线免费观看视频| 亚洲 欧美一区二区三区| 成年动漫av网址| 国产一区二区在线观看av| 一级黄片播放器| 亚洲成色77777| 亚洲国产中文字幕在线视频| 自线自在国产av| 久久99精品国语久久久| 咕卡用的链子| 午夜免费男女啪啪视频观看| 国产精品秋霞免费鲁丝片| 少妇人妻精品综合一区二区| 久久久国产精品麻豆| 欧美少妇被猛烈插入视频| 亚洲图色成人| 这个男人来自地球电影免费观看 | 秋霞伦理黄片| 日本wwww免费看| 欧美成人午夜精品| 国产有黄有色有爽视频| 黄片播放在线免费| 午夜影院在线不卡| 日本午夜av视频| 777久久人妻少妇嫩草av网站| 久久精品人人爽人人爽视色| 男女高潮啪啪啪动态图| 女人高潮潮喷娇喘18禁视频| 国产精品偷伦视频观看了| 2021少妇久久久久久久久久久| 久久精品国产亚洲av涩爱| 免费久久久久久久精品成人欧美视频| 久久久久久久国产电影| 男女边摸边吃奶| 亚洲综合精品二区| 91精品三级在线观看| 亚洲国产中文字幕在线视频| 亚洲av在线观看美女高潮| 大香蕉久久成人网| 女性被躁到高潮视频| 国产免费现黄频在线看| 欧美久久黑人一区二区| av在线app专区| 亚洲熟女精品中文字幕| 欧美精品亚洲一区二区| 天堂中文最新版在线下载| 91精品伊人久久大香线蕉| 午夜福利,免费看| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| 大片电影免费在线观看免费| 亚洲成人免费av在线播放| 亚洲美女视频黄频| 婷婷成人精品国产| 一级,二级,三级黄色视频| 成年人午夜在线观看视频| 成人亚洲欧美一区二区av| 日本av手机在线免费观看| 日韩电影二区| 各种免费的搞黄视频| 如何舔出高潮| 一级片免费观看大全| 日本欧美国产在线视频| 免费观看人在逋| 欧美乱码精品一区二区三区| 看十八女毛片水多多多| 天美传媒精品一区二区| 久久毛片免费看一区二区三区| 中文精品一卡2卡3卡4更新| 人人妻,人人澡人人爽秒播 | 亚洲专区中文字幕在线 | av女优亚洲男人天堂| 精品国产露脸久久av麻豆| 久久97久久精品| 欧美黑人欧美精品刺激| 精品一区在线观看国产| 精品一品国产午夜福利视频| 少妇的丰满在线观看| 91成人精品电影| 国产成人91sexporn| 嫩草影视91久久| 在线观看免费午夜福利视频| 久久久精品区二区三区| 汤姆久久久久久久影院中文字幕| 九色亚洲精品在线播放| 男的添女的下面高潮视频| 亚洲一区中文字幕在线| 在线看a的网站| 不卡av一区二区三区| 亚洲人成77777在线视频| 国产99久久九九免费精品| 精品少妇黑人巨大在线播放| 久久热在线av| 欧美老熟妇乱子伦牲交| 久久ye,这里只有精品| 久久毛片免费看一区二区三区| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 少妇的丰满在线观看| 亚洲国产毛片av蜜桃av| 狂野欧美激情性xxxx| 波多野结衣av一区二区av| 七月丁香在线播放| 老熟女久久久| 亚洲美女视频黄频| www.自偷自拍.com| 一区二区av电影网| 亚洲图色成人| e午夜精品久久久久久久| 国产在线一区二区三区精| 国产麻豆69| 国产乱人偷精品视频| 极品人妻少妇av视频| 丝瓜视频免费看黄片| 国产国语露脸激情在线看| 自线自在国产av| 精品一区二区三区四区五区乱码 | xxx大片免费视频| 男人添女人高潮全过程视频| 大片电影免费在线观看免费| 久久精品熟女亚洲av麻豆精品| 高清av免费在线| 午夜免费观看性视频| 欧美另类一区| 亚洲精品一二三| 精品一区二区三区av网在线观看 | 王馨瑶露胸无遮挡在线观看| 大陆偷拍与自拍| 老司机深夜福利视频在线观看 | 精品少妇一区二区三区视频日本电影 | 赤兔流量卡办理| 免费av中文字幕在线| 久久婷婷青草| 一级毛片电影观看| 亚洲第一青青草原| 综合色丁香网|