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

    基于RGB與高光譜圖像的小麥低溫脅迫評估模型

    2025-04-07 00:00:00余德炤江曉東楊瑩穎張建取忻樂張艷秦思容楊再強(qiáng)
    南方農(nóng)業(yè)學(xué)報 2025年1期

    DOI:10.3969/j.issn.2095-1191.2025.01.009

    摘要:【目的】探究低溫脅迫對小麥葉片葉綠素?zé)晒鈪?shù)、RGB圖像參數(shù)和高光譜指數(shù)的影響,建立小麥低溫脅迫評估模型,為小麥生產(chǎn)防災(zāi)減災(zāi)提供參考?!痉椒ā恳詽?jì)麥22為研究對象,在小麥拔節(jié)期開展低溫脅迫控制試驗(yàn),設(shè)白天(8:00—20:00)/夜間(20:00—次日8:00)平均溫度分別為8℃/0℃(T1)、6℃/-2℃(T2)和4℃/-4℃(T3)3個處理,持續(xù)時間3 d,以大田自然環(huán)境的盆栽小麥(23℃/8℃)為對照(CK),研究低溫脅迫處理結(jié)束后1、3和6 d小麥葉片葉綠素?zé)晒鈪?shù)、RGB圖像參數(shù)及高光譜指數(shù)的變化規(guī)律;使用一元線性回歸、隨機(jī)森林(RF)和人工神經(jīng)網(wǎng)絡(luò)(ANN)建立小麥低溫脅迫評估模型?!窘Y(jié)果】葉綠素?zé)晒鈪?shù)DIo/RC可作為評估小麥低溫脅迫的指標(biāo)。在一元線性回歸模型中,使用增強(qiáng)型植被指數(shù)(EVI)建立的一元線性回歸模型效果最佳,回歸方程為y=-1.261x+1.401,決定系數(shù)(R2)、均方根誤差(RMSE)、平均絕對誤差(MAE)、平均相對誤差(MRE)分別為0.536、0.058、0.045和11.31%。在RF和ANN模型中,基于RGB圖像參數(shù)建立的模型精度高于基于高光譜指數(shù)建立的模型,RF模型測試集R2、RMSE、MAE、MRE分別為0.771、0.042、0.033、8.57%,R2相比一元線性回歸模型提高43.78%,RMSE、MAE、MRE分別降低28.31%、28.06%、24.21%;ANN模型測試集R2、RMSE、MAE、MRE分別為0.742、0.046、0.037、9.01%,測試集R2相比一元線性回歸模型提高38.34%,RMSE、MAE、MRE分別降低20.33%、18.06%、20.32%?!窘Y(jié)論】基于RGB圖像參數(shù)的RF模型效果最好、精度最高,可用于小麥的低溫脅迫評估。

    關(guān)鍵詞:小麥;葉綠素?zé)晒鈪?shù);RGB圖像;高光譜圖像;低溫脅迫評估模型

    中圖分類號:S512.1文獻(xiàn)標(biāo)志碼:A文章編號:2095-1191(2025)01-0097-14

    Assessment model for low temperature stress in wheat based on RGB and hyperspectral images

    YU De-zhao1,JIANG Xiao-dong1*,YANG Ying-ying1,ZHANG Jian-qu1,XIN Le1,ZHANG Yan2,QIN Si-rong1,YANG Zai-qiang1

    (1School of Ecology and Applied Meteorology,Nanjing University of Information Scienceamp;Technology/Jiangsu Provin‐cial University Key laboratory of Agricultural and Ecological Meteorology,Nanjing,Jiangsu 210044,China;2School ofAtmospheric Sciences,Nanjing University of Information Scienceamp;Technology,Nanjing,Jiangsu 210044,China)

    Abstract:【Objective】This study investigated the effects of low temperature stress on chlorophyll fluorescence para-meters,RGB image parameters and hyperspectral indexes of wheat leaves.It aimed to establish an assessment model for wheat under low temperature stress,providing reference for disaster prevention and mitigation in wheat production.【Method】Using winter wheat variety Jimai 22 as the research material,a controlled low temperature stress experiment was conducted during the jointing stage.Three treatments were applied:daytime(8:00–20:00)/nighttime(20:00–8:00 on next day)mean temperatures of 8°C/0℃(T1),6°C/-2℃(T2)and 4°C/-4℃(T3),each lasting for three days.Potted wheat grown under natural field conditions(23°C/8℃)served as the control(CK).The changes in chlorophyllfluorescence parameters,RGB image parameters,and hyperspectral indexes of wheat leaves were analyzed on 1,3 and 6 d after the low temperature stress treatments.Wheat low temperature stress assessment models were developed using uni‐variate linear regression,random forest(RF)and artificial neural networks(ANN).【Result】The chlorophyll fluores‐cence parameter DIo/RC was identified as an effective indicator for assessing wheat low temperature stress.In the univaria-te linear regression model,the model using the enhanced vegetation index(EVI)performed the best,with a regression equation of y=-1.261x+1.401,yielding a coefficient of determination(R2),root mean square error(RMSE),mean abso-lute error(MAE)and mean relative error(MRE)of 0.536,0.058,0.045,and 11.31%respectively.For RF and ANN models,models based on RGB image parameters outperformed those based on hyperspectral indexes.The RF model achieved R2,RMSE,MAE and MRE values of 0.771,0.042,0.033 and 8.57%in the test set,with an R2 improvement of 43.78%and reductions in RMSE,MAE and MRE by 28.31%,28.06%and 24.21%respectively compared to the univaria-te linear regression model.The ANN model achieved R2,RMSE,MAE and MRE values of 0.742,0.046,0.037 and 9.01%,with R2 improving by 38.34%and RMSE,MAE and MRE decreasing by 20.33%,18.06%and 20.32%respec-tively compared to the univariate linear regression model.【Conclusion】The RF model based on RGB image parameters demonstrates the best performance and the highest accuracy,making it suitable for evaluating low temperature stress in wheat.

    Key words:wheat;chlorophyll fluorescence parameters;RGB image;hyperspectral image;low temperature stress assessment model

    Foundation items:National Key Research and Development Program of China(2022YFD2300202)

    0引言

    【研究意義】小麥?zhǔn)俏覈蠹Z食作物之一,對保障國家糧食安全具有重要意義。氣候變化背景下,極端天氣事件頻發(fā)(周波濤和錢進(jìn),2021),低溫脅迫風(fēng)險增大,嚴(yán)重影響我國的小麥生產(chǎn)(陳翔等,2021;代存虎等,2024)。調(diào)查發(fā)現(xiàn),我國85%的冬小麥種植區(qū)是低溫災(zāi)害的風(fēng)險區(qū)域(Yue et al.,2016),其中黃淮麥區(qū)是低溫災(zāi)害最嚴(yán)重的區(qū)域(Xiao et al.,2018),部分省份在低溫災(zāi)害嚴(yán)重年份可減產(chǎn)30%~50%(Liu et al.,2019)。低溫災(zāi)害已成為限制我國小麥產(chǎn)量的重要因素之一,開展小麥低溫脅迫監(jiān)測與評估研究刻不容緩。隨著信息技術(shù)和圖像處理技術(shù)的發(fā)展,基于RGB、高光譜和熒光圖像的計算機(jī)視覺技術(shù)為農(nóng)作物的低溫脅迫監(jiān)測提供了一種無接觸式的可行方法(Shammi etal.,2023)。相比傳統(tǒng)的人工方法,計算機(jī)視覺技術(shù)具有檢測效率高、時效性強(qiáng)、覆蓋范圍廣的優(yōu)點(diǎn)(Tian et al.,2020)。此外,結(jié)合計算機(jī)視覺技術(shù),可深入挖掘作物光學(xué)特征與生理指標(biāo)之間的關(guān)系,為低溫脅迫評估提供更多維度的數(shù)據(jù)支持與模型優(yōu)化的可能性(Li etal.,2020)。因此,開展基于計算機(jī)視覺技術(shù)的小麥低溫脅迫評估研究,對提升小麥低溫脅迫監(jiān)測能力、降低災(zāi)害損失以及保障國家糧食安全具有重要意義?!厩叭搜芯窟M(jìn)展】近年來,國內(nèi)外學(xué)者已廣泛基于葉綠素?zé)晒狻GB和高光譜圖像開展了小麥低溫脅迫的評估與監(jiān)測研究??焖偃~綠素?zé)晒庹T導(dǎo)動力學(xué)曲線蘊(yùn)含著豐富的光合信息,常被用于研究植物逆境脅迫(楊文強(qiáng)等,2024;Pessarakli,2024)。小麥?zhǔn)艿降蜏孛{迫后,光合作用受到影響,葉綠素?zé)晒鈪?shù)值發(fā)生改變(Dong et al.,2020;岳俊芹等,2021)。Li等(2015)研究發(fā)現(xiàn),拔節(jié)期低溫脅迫后,小麥葉綠素?zé)晒鈪?shù)ETo/RC和TRo/RC均顯著降低。張磊等(2017)、劉蕾蕾等(2018)、王瑞霞等(2018)研究均發(fā)現(xiàn),拔節(jié)期低溫脅迫處理后小麥葉綠素?zé)晒鈪?shù)Fv/Fm顯著降低,是評估低溫脅迫對小麥光合作用損傷的重要指標(biāo)。楊程等(2022)通過分析小麥越冬期34個葉綠素?zé)晒庹T導(dǎo)動力學(xué)參數(shù)與33個小麥品種凍害發(fā)生程度的相關(guān)性,發(fā)現(xiàn)熒光參數(shù)TRo/CSm比Fv/Fm更適合用于評價小麥凍害情況。RGB圖像技術(shù)憑借其簡單、廉價且高效的特點(diǎn),近年來逐漸成為作物低溫脅迫監(jiān)測的重要工具(Yuan and Choi,2021;Ku et al.,2023;樊小雪等,2024)。RGB圖像能直觀反映小麥的外觀特征,可準(zhǔn)確捕捉低溫脅迫下葉片顏色和形態(tài)變化(薛亞光等,2018),為小麥低溫災(zāi)害的監(jiān)測與評估提供了重要依據(jù)。Macedo-Cruz等(2011)通過獲取健康和凍害地區(qū)燕麥作物的RGB圖像,結(jié)合3種閾值處理方法,實(shí)現(xiàn)了燕麥作物凍害情況的無監(jiān)督分類。Fu等(2021)采集小麥RGB圖像后提取影響作物生長的各種重要指標(biāo),發(fā)現(xiàn)利用小麥苗期植被覆蓋度作為評價指標(biāo)對小麥凍害程度進(jìn)行分級具有良好效果。Zhang等(2024a)利用小麥冠層RGB圖像提取顏色特征,量化了小麥葉片在低溫脅迫下的黃化程度,驗(yàn)證了其對低溫脅迫監(jiān)測的可行性。小麥?zhǔn)艿降蜏孛{迫后,葉片對特定光譜波段的吸收強(qiáng)度會減弱,光譜特征發(fā)生變化(Zhang et al.,2024b)。Wang等(2016)、Zhao等(2020)結(jié)合歸一化植被指數(shù)(NDVI)評估了低溫脅迫對冬小麥生產(chǎn)的潛在損失,為低溫災(zāi)害監(jiān)測提供了技術(shù)支持。Feng等(2018)、Xie等(2020)篩選出低溫脅迫下小麥高光譜的敏感波段,揭示了小麥葉片在低溫脅迫條件下的光譜響應(yīng)特征。Wang等(2020)提出了基于遙感的春霜凍損害指數(shù),快速有效地量化了大范圍春霜凍對小麥的影響。Yao等(2024)構(gòu)建了基于高光譜圖像的小麥幼苗期凍害評估模型,模型精度超過86%。上述研究在一定程度上實(shí)現(xiàn)了小麥的低溫脅迫評估與監(jiān)測,快速葉綠素?zé)晒狻GB圖像和高光譜圖像3種技術(shù)各有優(yōu)勢,葉綠素?zé)晒鈪?shù)能表征小麥光合效率的變化,RGB圖像能直觀反映低溫脅迫對小麥外觀形態(tài)的影響,高光譜圖像的多維信息能深入揭示低溫脅迫對小麥的生理響應(yīng)?!颈狙芯壳腥朦c(diǎn)】不同圖像在特征選擇、模型構(gòu)建和應(yīng)用場景上存在一定局限性,前人研究在評估小麥低溫脅迫時多集中于單一數(shù)據(jù)源的分析,對不同技術(shù)融合應(yīng)用的研究較少,且在測量快速葉綠素?zé)晒庹T導(dǎo)動力學(xué)曲線時需進(jìn)行15~20 min的遮光處理,在進(jìn)行大面積監(jiān)測時不能及時獲取熒光數(shù)據(jù)。目前,鮮見使用RGB和高光譜圖像特征估算葉綠素?zé)晒鈪?shù)進(jìn)行小麥低溫脅迫評估的研究報道。【擬解決的關(guān)鍵問題】在小麥拔節(jié)期開展低溫脅迫控制試驗(yàn),測定并分析低溫脅迫處理結(jié)束后不同時間小麥葉片葉綠素?zé)晒鈪?shù)、RGB圖像參數(shù)及高光譜指數(shù)的變化規(guī)律,基于小麥葉片RGB圖像參數(shù)和高光譜指數(shù),利用一元線性回歸、隨機(jī)森林(RF)和人工神經(jīng)網(wǎng)絡(luò)(ANN)建立小麥低溫脅迫評估模型,并對模型的可行性和精準(zhǔn)度進(jìn)行評價,以期為小麥生產(chǎn)防災(zāi)減災(zāi)提供參考。

    1材料與方法

    1.1試驗(yàn)材料

    供試小麥品種為濟(jì)麥22,供試土壤采集自南京信息工程大學(xué)農(nóng)業(yè)氣象實(shí)驗(yàn)站小麥試驗(yàn)地耕層,土壤質(zhì)地為壤質(zhì)黏土,黏粒含量26.5%,土壤pH 6.4,有機(jī)碳、全氮含量分別為20.5和10.8 g/kg。

    1.2試驗(yàn)方法

    1.2.1試驗(yàn)設(shè)計試驗(yàn)于2023年11月—2024年5月在南京信息工程大學(xué)農(nóng)業(yè)氣象試驗(yàn)站進(jìn)行。小麥種植于直徑20 cm、高30 cm的聚乙烯盆中,每盆裝試驗(yàn)田耕層土5kg。小麥于2023年11月21日播種,2023年12月27日3葉期定苗,5株/盆。播種前施用基肥,每盆肥料用量為三元復(fù)合肥(N-P2O5-K2O=15-15-15)1.0 g、尿素1.5 g,在拔節(jié)期追肥,每盆施尿素1.5 g,其他管理措施同當(dāng)?shù)剞r(nóng)田。

    于拔節(jié)期(2024年3月15日)將小麥移入人工氣候箱中進(jìn)行為期3d(3月15日20:00至3月18日20:00)的低溫脅迫處理(圖1-A),設(shè)白天(8:00—20:00)/夜間(20:00—次日8:00)平均溫度分別為8℃/0℃(T1)、6℃/-2℃(T2)和4℃/-4℃(T3)3個處理,以大田自然環(huán)境的盆栽小麥(23℃/8℃)為對照(CK)。人工氣候箱相對濕度設(shè)置為(70±5)%,白天光照強(qiáng)度設(shè)置為800μmol/(m2?s),夜間光照強(qiáng)度設(shè)置為0μmol/(m2?s),低溫脅迫期間人工氣候箱內(nèi)溫度變化如圖1-B。每處理設(shè)3個重復(fù),每重復(fù)種植10盆,在低溫脅迫處理結(jié)束后當(dāng)天將各處理盆栽重新移回大田自然環(huán)境,直至成熟。

    1.2.2熒光參數(shù)測定在低溫脅迫處理結(jié)束后的1、3和6 d,以小麥主莖頂數(shù)第2~4片葉為觀測葉片,利用Pocket PEA植物效率分析儀(PP Systems公司,美國)測定葉片快速葉綠素?zé)晒庹T導(dǎo)動力學(xué)曲線,計算葉綠素?zé)晒鈪?shù)(Goltsev etal.,2017;曹雄軍等,2024)(表1)。

    1.2.3 RGB圖像采集及特征提取小麥葉片RGB圖像采集及特征提取過程見圖2,使用realme GT手機(jī)(深圳市銳爾覓移動通信有限公司)拍攝小麥葉片RGB圖像,共獲取有效RGB圖像229張。使用Python 3.12.1編程,借助OpenCV 4.6.0工具包對小麥葉片RGB圖像進(jìn)行預(yù)處理,部分圖像見圖3,并提取RGB圖像的顏色特征和紋理特征(楊俊等,2019;Chen et al.,2020;Lee et al.,2021)(表2)。

    1.2.4高光譜圖像采集及特征提取小麥葉片高光譜圖像采集及處理過程見圖4,采用SOC710便攜式高光譜成像儀(Surface Optics Corporation公司,美國)測定小麥葉片400~1000 nm高光譜圖像,共獲取有效高光譜圖像206張。使用SRAnal710e進(jìn)行光譜輻射定標(biāo),顯示小麥葉片RGB格式下的高光譜合成圖像,如圖5所示,使用ENVI 5.3進(jìn)行反射光譜提取,并計算相關(guān)高光譜特征參數(shù)和植被指數(shù)(表3和表4)(Qiaoetal.,2022;楊熙來等,2023)。

    1.3數(shù)據(jù)處理

    1.3.1建立模型將RGB圖像參數(shù)、高光譜指數(shù)與葉綠素?zé)晒鈪?shù)進(jìn)行線性回歸分析,建立一元線性回歸模型。使用RF和ANN建立回歸模型。RF決策樹數(shù)量為50,最大深度為10,節(jié)點(diǎn)最小分裂樣本數(shù)為5,葉子節(jié)點(diǎn)最小樣本數(shù)為1。ANN包含3層全連接層,用于預(yù)測目標(biāo)變量,其中,輸入層接收歸一化處理后的輸入特征(RGB圖像參數(shù)或高光譜指數(shù)),隱藏層分別包含32和16個神經(jīng)元,輸出層為單一節(jié)點(diǎn)用于生成預(yù)測值,模型的激活函數(shù)選用ReLU,以引入非線性特性。訓(xùn)練過程中,定義均方誤差(MSE)作為損失函數(shù),采用Adam優(yōu)化器更新網(wǎng)絡(luò)參數(shù),初始學(xué)習(xí)率設(shè)為0.01。為了提高收斂效果,引入學(xué)習(xí)率調(diào)度策略,每10次迭代后將學(xué)習(xí)率調(diào)整為原來的0.9倍。訓(xùn)練過程設(shè)置最大迭代次數(shù)為200,同時引入早停機(jī)制,若測試集的損失連續(xù)50次迭代無改善,則提前停止訓(xùn)練以防止過擬合。在使用RF和ANN模型訓(xùn)練時,將樣本數(shù)據(jù)進(jìn)行隨機(jī)劃分,其中70%樣本量用于訓(xùn)練,30%樣本量用于驗(yàn)證。

    1.3.2模型檢驗(yàn)使用決定系數(shù)(R2)、均方根誤差(RMSE)、平均絕對誤差(MAE)、平均相對誤差(MRE)作為評價指標(biāo),R2越大,RMSE、MAE和MRE越小代表模型精度越高,公式如下:

    式中,xi為實(shí)測值,x(ˉ)為實(shí)測值的平均值,yi為預(yù)測值,n為樣本個數(shù)。

    1.3.3統(tǒng)計分析使用Python 3.12.1、Excel 2021進(jìn)行數(shù)據(jù)處理,借助sklearn、PyTorch工具包進(jìn)行模型構(gòu)建,使用SPSS 21.0進(jìn)行方差分析及相關(guān)分析,使用Origin 2021繪圖。

    2結(jié)果與分析

    2.1不同低溫脅迫處理小麥葉片葉綠素?zé)晒鈪?shù)響應(yīng)特征

    以CK為標(biāo)準(zhǔn),將各處理數(shù)據(jù)歸一化,繪制不同處理葉綠素?zé)晒鈪?shù)雷達(dá)圖(圖6)。通過方差分析可知,在低溫脅迫處理結(jié)束后1 d,各處理TRo/RC、ETo/RC較CK顯著降低(Plt;0.05,下同),DIo/RC較CK顯著升高;在低溫脅迫處理結(jié)束后3 d,各處理Fv/Fm較CK顯著降低,DIo/RC較CK顯著升高;在低溫脅迫處理結(jié)束后6 d,各處理Fv/Fm、PIABS較CK顯著降低,DIo/RC較CK顯著升高,說明葉綠素?zé)晒鈪?shù)TRo/RC、ETo/RC、Fv/Fm、PIABS、DIo/RC可有效反映小麥的低溫脅迫狀況。

    將小麥葉綠素?zé)晒鈪?shù)與處理最低溫度進(jìn)行相關(guān)系數(shù)分析(圖7),結(jié)果顯示,F(xiàn)v/Fm、TRo/CSm和DIo/RC 3個參數(shù)相關(guān)性最強(qiáng),在低溫脅迫處理結(jié)束后1 d,3個參數(shù)與處理最低溫度的相關(guān)系數(shù)分別為0.70、0.63和-0.86;在低溫脅迫處理結(jié)束后3 d,相關(guān)系數(shù)分別為0.86、0.68和-0.86;在低溫脅迫處理結(jié)束后6 d,相關(guān)系數(shù)分別為0.99、0.85和-0.88,且均達(dá)到顯著水平。DIo/RC在各處理低溫脅迫處理結(jié)束后不同時間均較CK顯著升高,且與處理最低溫度均顯著相關(guān)(|r|gt;0.85),因此,本研究選取葉綠素?zé)晒鈪?shù)DIo/RC作為評估小麥低溫脅迫的指標(biāo)。

    2.2不同低溫脅迫處理小麥葉片RGB圖像特征

    不同低溫脅迫處理小麥葉片RGB圖像主要顏色特征及紋理特征見表5。在低溫脅迫處理結(jié)束后1 d,T2和T3處理小麥葉片RGB圖像顏色特征RRsd、SRsd、IRsd較CK顯著升高,且T3處理RRsd顯著高于T2處理。在低溫脅迫處理結(jié)束后3 d,各處理小麥葉片RGB圖像顏色特征Rmean、Gmean較CK均顯著降低,且在不同處理間呈現(xiàn)出T2gt;T1gt;T3。在低溫脅迫處理結(jié)束后6 d,各處理小麥葉片RGB圖像特征RRsd、GRsd、SRsd、IRsd較CK均顯著升高,Gmean、Correlation較CK均顯著降低。

    T3處理顏色特征g和ExG在低溫脅迫處理結(jié)束后1、3、6 d均與其他處理差異顯著,其中,在低溫脅迫處理結(jié)束后1 d時T3處理顯著低于其他處理,而3和6 d時T3處理顯著高于其他處理,表明T3處理對葉片綠色分量的影響隨時間發(fā)生變化。在低溫脅迫處理結(jié)束后1 d時,CK和T3處理的紋理特征Cor-relation顯著高于T1和T2處理,在低溫脅迫處理結(jié)束后3 d時,T3處理的紋理特征Correlation顯著低于其他處理,而在低溫脅迫處理結(jié)束后6 d時,各處理紋理特征Correlation均顯著低于CK,說明紋理特征Correlation對不同程度的低溫脅迫具有一定的區(qū)分性。

    2.3不同低溫脅迫處理小麥葉片高光譜圖像響應(yīng)特征

    低溫脅迫處理結(jié)束后1、3和6 d,不同處理與CK小麥葉片光譜反射率曲線(圖8)的形狀和變化趨勢基本相似。在可見光波段,低溫脅迫處理結(jié)束后1 d(圖8-A)和3 d(圖8-B),各處理反射率與CK差異較小。在低溫脅迫處理結(jié)束后6 d(圖8-C),T3處理綠峰附近反射率明顯高于CK,且伴隨有紅移現(xiàn)象,即綠峰中心位置向長波方向偏移,說明低溫脅迫下小麥葉片對自然光的吸收處于低值,小麥正常的生理活動受到抑制。在波長800~1000 nm內(nèi)的近紅外波段,低溫脅迫處理結(jié)束后1 d,小麥葉片平均反射率大小順序?yàn)镃Kgt;T1gt;T3gt;T2,分別為0.466、0.464、0.460和0.456;低溫脅迫處理結(jié)束后3 d,小麥葉片平均反射率大小順序?yàn)镃Kgt;T1gt;T2gt;T3,分別為0.560、0.558、0.546和0.529;低溫脅迫處理結(jié)束后6 d,小麥葉片平均反射率大小順序?yàn)镃Kgt;T2gt;T1gt;T3,分別為0.535、0.532、0.519和0.511。低溫脅迫處理結(jié)束后3 d各處理小麥葉片平均反射率較CK降低,但在低溫脅迫處理結(jié)束后6 d降幅減小,可能是低溫脅迫后小麥葉片細(xì)胞結(jié)構(gòu)受損,導(dǎo)致反射率降低,處理結(jié)束后6d小麥細(xì)胞結(jié)構(gòu)又有所恢復(fù)。

    從小麥葉片一階微分光譜曲線(圖9)可知,低溫脅迫處理結(jié)束后1 d(圖9-A)和3 d(圖9-B),各處理與CK差異較小,在低溫脅迫處理結(jié)束后6 d(圖9-C),T1、T2處理與CK差異較小,T3處理藍(lán)邊處一階微分光譜反射率為0.0023 nm-1,相比CK(0.0019 nm-1)明顯增大;T3處理紅邊處一階微分光譜反射率為0.0081 nm-1,較CK(0.0090 nm-1)明顯減小,且紅邊位置向短波方向移動,出現(xiàn)藍(lán)移現(xiàn)象,可見,低溫脅迫一定程度上導(dǎo)致小麥葉片紅邊藍(lán)移現(xiàn)象。

    2.4基于RGB圖像參數(shù)和高光譜指數(shù)的小麥低溫脅迫評估模型的構(gòu)建與檢驗(yàn)

    2.4.1一元線性回歸模型將葉綠素?zé)晒鈪?shù)DIo/RC與RGB圖像參數(shù)、高光譜指數(shù)分別進(jìn)行相關(guān)分析(圖10),選取顯著相關(guān)(Plt;0.05,|r|gt;0.6)的特征構(gòu)建一元線性回歸模型(表6)。結(jié)果顯示,RGB圖像參數(shù)中,使用紋理特征Correlation構(gòu)建的模型效果最好,回歸方程為y=-8.443x+8.722,R2、RMSE、MAE、MRE分別為0.470、0.062、0.052、12.60%。使用ExG和g構(gòu)建的模型R2、RMSE、MAE、MRE均分別為0.442、0.064、0.053、12.88%,原因是2個特征參數(shù)存在線性關(guān)系ExG=3g-1(楊俊等,2019)。高光譜指數(shù)中,使用植被指數(shù)EVI構(gòu)建的模型效果最好,回歸方程為y=-1.261x+1.401,R2、RMSE、MAE、MRE分別為0.536、0.058、0.045、11.31%。

    2.4.2 RF和ANN回歸模型利用RF和ANN分別構(gòu)建葉綠素?zé)晒鈪?shù)DIo/RC與RGB圖像參數(shù)、高光譜指數(shù)的回歸模型,結(jié)果見表7。

    2.4.3模型的評估與比較在模型的評估指標(biāo)中,R2反映實(shí)測值與預(yù)測值間的相關(guān)強(qiáng)度,RMSE、MAE、MRE則反映實(shí)測值與預(yù)測值間的偏差,最佳的估算模型應(yīng)選擇R2較大,而RMSE、MAE、MRE較小的模型。在使用一元線性回歸時,基于高光譜指數(shù)的估算模型相比基于RGB圖像參數(shù)的效果更好。而RF和ANN模型中,基于RGB圖像參數(shù)的估算模型效果均優(yōu)于基于高光譜指數(shù)的模型,說明不同模型選擇的最優(yōu)參數(shù)不同。

    RF和ANN模型測試集擬合結(jié)果見圖11,RF和ANN模型均取得較好效果,其中基于RGB圖像參數(shù)的RF模型效果最好,精度最高,測試集R2相比最優(yōu)一元線性回歸模型提高43.78%,RMSE、MAE、MRE分別降低28.31%、28.06%、24.21%?;赗GB圖像參數(shù)的ANN模型精度高于一元線性回歸模型,測試集R2相比最優(yōu)一元線性回歸模型提高38.34%,RMSE、MAE、MRE分別降低20.33%、18.06%、20.32%。綜合比較認(rèn)為基于RGB圖像參數(shù)的RF模型為最優(yōu)DIo/RC估算模型。

    3討論

    3.1低溫脅迫對小麥葉片葉綠素?zé)晒鈪?shù)的影響

    葉綠素?zé)晒庹T導(dǎo)動力學(xué)曲線蘊(yùn)含豐富的光合作用相關(guān)信息,對低溫脅迫響應(yīng)敏感(劉蕾蕾等,2018;Dong et al.,2020;方宇輝等,2023)。岳俊芹等(2021)、付延松等(2022)的研究表明,拔節(jié)期低溫脅迫后,小麥葉片葉綠素?zé)晒鈪?shù)Fv/Fm隨溫度的降低而下降。楊程等(2022)研究自然大田環(huán)境下33個小麥品種凍害情況,發(fā)現(xiàn)葉綠素?zé)晒鈪?shù)Fv/Fm、TRo/CSm和DIo/RC等均可作為評價小麥抗寒性的指標(biāo),其中TRo/CSm評價效果最好,說明利用葉綠素?zé)晒鈪?shù)評價小麥低溫脅迫的可行性。本研究結(jié)果亦表明,在受到低溫脅迫后,F(xiàn)v/Fm和TRo/CSm顯著降低,說明小麥?zhǔn)艿蜏孛{迫后葉片的光能轉(zhuǎn)化效率和捕獲光能的能力降低,而以熱能形式耗散的能量增多,光合活性受到抑制。Fv/Fm、TRo/CSm和DIo/RC 3個參數(shù)中,DIo/RC在低溫脅迫處理結(jié)束后1、3和6 d的數(shù)值均較CK顯著升高,說明DIo/RC能更靈敏地表征小麥低溫脅迫的受害情況。

    3.2低溫脅迫對小麥葉片RGB圖像和反射光譜特征的影響

    小麥葉片RGB圖像是在可見光范圍內(nèi)捕獲的彩色圖像,包含紅、綠、藍(lán)3個通道,可提供葉片顏色、紋理等視覺特征參數(shù)。小麥?zhǔn)艿降蜏孛{迫后,葉片會出現(xiàn)枯萎、斑點(diǎn)、發(fā)黃等現(xiàn)象導(dǎo)致外觀形態(tài)改變(薛亞光等,2018),從而使利用葉片圖像特征檢測低溫脅迫具有可能性。Fu等(2021)研究表明,低溫脅迫下小麥RGB圖像特征ExGR變化明顯,可用于評估小麥的受害狀況。Zhang等(2024a)通過分析小麥冠層圖像的顏色分布變化,發(fā)現(xiàn)隨著低溫脅迫的加重,葉片顏色從亮綠色轉(zhuǎn)變?yōu)辄S綠色和黃色,且顏色直方圖的扁平化和偏斜程度與低溫脅迫嚴(yán)重程度顯著相關(guān)。在本研究中,低溫脅迫后,小麥葉片RGB圖像顏色特征Gmean較CK明顯降低、GRsd較CK明顯升高,圖像綠色通道均值減小,相對標(biāo)準(zhǔn)偏差增大,反映了小麥葉片RGB圖像特征的變化。

    小麥反射光譜主要受其葉片內(nèi)部細(xì)胞結(jié)構(gòu)、色素含量等對光波的反射和吸收影響,與作物的生理生化成分相關(guān)(Silva-Perez et al.,2018;Zhang et al.,2024b)。Murphy等(2020)設(shè)置小麥低溫控制試驗(yàn),發(fā)現(xiàn)小麥旗葉在藍(lán)(415~494 nm)、紅(670~687 nm)和近紅外(727~889 nm)區(qū)域的反射特性均受到低溫脅迫的顯著影響,并在后續(xù)的5 d內(nèi)表現(xiàn)出低溫脅迫處理與時間的顯著交互作用。Xie等(2020)的研究結(jié)果表明,低溫脅迫會導(dǎo)致小麥冠層光譜的紅邊區(qū)域發(fā)生藍(lán)移,近紅外光譜反射率降低,說明小麥冠層光譜對低溫脅迫的敏感性。在本研究中,低溫脅迫下的小麥葉片顏色變黃,在可見光波段吸收減少,導(dǎo)致反射率增大,綠峰向長波方向偏移,而在800~1000 nm內(nèi)的近紅外波段,葉片反射率降低,紅邊位置向短波方向移動,出現(xiàn)藍(lán)移現(xiàn)象,同樣說明低溫脅迫對葉片細(xì)胞結(jié)構(gòu)造成了損傷。

    3.3基于RGB圖像參數(shù)和高光譜指數(shù)的小麥低溫脅迫評估模型

    一元線性回歸、RF和ANN 3種不同的建模方法比較,使用RF建模因其集成了多個決策樹來提高模型的穩(wěn)定性和準(zhǔn)確性(沈潤平等,2017),降低了數(shù)據(jù)誤差對模型的影響(Mellor et al.,2015),同時相較于一元線性回歸模型,其具有處理復(fù)雜非線性關(guān)系的能力,而其相較于ANN建模更不容易出現(xiàn)過擬合現(xiàn)象(Belgiu and Dr?gu?,2016),從而取得了最高的模型精度。Shammi等(2024)利用RF等機(jī)器學(xué)習(xí)方法基于紅外熱成像數(shù)據(jù)構(gòu)建了小麥霜凍早期檢測模型,模型準(zhǔn)確率達(dá)98%,驗(yàn)證了使用RF模型檢測小麥霜凍的合理性。

    本研究中,將RGB圖像參數(shù)、高光譜指數(shù)代入不同的模型進(jìn)行了比較,在RF和ANN模型中,基于RGB圖像參數(shù)建立的模型精度均高于基于高光譜指數(shù)建立的模型,主要是由于RGB圖像能直接反映小麥葉片的生理狀態(tài),故基于RGB圖像提取的特征能更準(zhǔn)確估算小麥葉片葉綠素?zé)晒鈪?shù)。相較于高光譜技術(shù),RGB圖像更容易獲取,成本較低,方便大范圍推廣使用。本研究使用RGB圖像構(gòu)建的模型取得了較好的效果,為高效、無損監(jiān)測小麥低溫脅迫提供了一定參考,但由于試驗(yàn)地點(diǎn)和供試品種的限制,模型能否用于其他區(qū)域及其他小麥品種還需進(jìn)一步驗(yàn)證。

    4結(jié)論

    低溫脅迫下小麥葉綠素?zé)晒鈪?shù)的變化及各葉綠素?zé)晒鈪?shù)與處理最低溫度的相關(guān)性分析結(jié)果表明,DIo/RC可作為評估小麥低溫脅迫的指標(biāo)。在一元線性回歸模型中,使用高光譜指數(shù)EVI建立的一元線性回歸模型效果最佳。在RF和ANN模型中,基于RGB圖像參數(shù)建立的模型精度高于基于高光譜指數(shù)建立的模型,效果均優(yōu)于一元線性回歸模型?;赗GB圖像參數(shù)的RF模型效果最好、精度最高,可用于小麥的低溫脅迫評估。

    參考文獻(xiàn)(References):

    曹雄軍,韓佳宇,謝蜀豫,黃秋秘,鄧海燕,盤豐平,王博,江春分,時曉芳,白先進(jìn).2024.葡萄葉綠素?zé)晒鈪?shù)測量方法的優(yōu)化與日變化特征分析[J].南方農(nóng)業(yè)學(xué)報,55(8):2248-2261.[Cao X J,Han J Y,Xie S Y,Huang Q M,Deng H Y,Pan F P,Wang B,Jiang C F,Shi X F,Bai X J.2024.Optimization of measurement methods for chloro‐phyll fluorescence parameters and analysis of diurnal varia‐tion characteristics in grape leaves[J].Journal of Southern Agriculture,55(8):2248-2261.]doi:10.3969/j.issn.2095-1191.2024.08.005.

    陳翔,于敏,蔡洪梅,吳宇,張樂樂,柯媛媛,許輝,李金才.2021.小麥倒春寒研究現(xiàn)狀與進(jìn)展[J].應(yīng)用生態(tài)學(xué)報,32(8):2999-3009.[Chen X,Yu M,Cai H M,Wu Y,Zhang L L,Ke Y Y,Xu H,Li J C.2021.Current status andresearch advances of late spring coldness in wheat[J].Chi‐nese Journal of Applied Ecology,32(8):2999-3009.]doi:10.13287/j.1001-9332.202108.021.

    代存虎,劉新,戶淑君,朱洋陽,楊舒然,朱敏,李春燕,丁錦峰,朱新開,郭文善.2024.拔節(jié)期低溫與施氮防御對小麥根系酶活性及產(chǎn)量形成的影響[J].江蘇農(nóng)業(yè)科學(xué),52(12):69-75.[Dai C H,Liu X,Hu S J,Zhu YY,Yang S R,Zhu M,Li C Y,Ding J F,Zhu X K,Guo W S.2024.Effects of low temperature and nitrogen application before low temperature at jointing stage on root enzyme activity and yield formation of wheat[J].Jiangsu Agricul‐tural Sciences,52(12):69-75.]doi:10.15889/j.issn.1002-1302.2024.12.009.

    樊小雪,李德翠,李遠(yuǎn),任妮.2024.基于RGB模型的草莓葉片光合作用指標(biāo)估測[J].江蘇農(nóng)業(yè)學(xué)報,40(4):675-681.[Fan X X,Li D C,Li Y,Ren N.2024.Estimation of photo‐synthetic indexes in strawberry leaves based on RGB model[J].Jiangsu Journal of Agricultural Sciences,40(4):675-681.]doi:10.3969/j.issn.1000-4440.2024.04.011.

    方宇輝,華夏,韓留鵬,趙明忠,齊學(xué)禮,董海濱,胡琳.2023.非生物脅迫因素對小麥光合作用的影響研究進(jìn)展[J].河南農(nóng)業(yè)科學(xué),52(10):1-13.[Fang Y H,Hua X,Han L P,Zhao M Z,Qi X L,Dong H B,Hu L.2023.Research progress on the effects of abiotic stress factors on wheat photosynthesis[J].Journal of Henan Agricultural Scien-ces,52(10):1-13.]doi:10.15933/j.cnki.1004-3268.2023.10.001.

    付延松,譚植,李文陽.2022.拔節(jié)期低溫對小麥籽粒產(chǎn)量與灌漿期旗葉熒光參數(shù)的影響[J].生態(tài)科學(xué),41(6):33-40.[Fu Y S,Tan Z,Li W Y.2022.Effects of low tempera‐ture at jointing stage on grain yield and fluorescence parameters of flag leaves in wheat grain filling[J].Eco‐logical Science,41(6):33-40.]doi:10.14108/j.cnki.1008-8873.2022.06.005.

    劉蕾蕾,紀(jì)洪亭,劉兵,馬吉鋒,肖瀏駿,湯亮,曹衛(wèi)星,朱艷.2018.拔節(jié)期和孕穗期低溫處理對小麥葉片光合及葉綠素?zé)晒馓匦缘挠绊懀跩].中國農(nóng)業(yè)科學(xué),51(23):4434-4448.[Liu L L,Ji H T,Liu B,Ma J F,Xiao L J,Tang L,Cao W X,Zhu Y.2018.Effects of jointing and booting low temperature treatments on photosynthetic and chloro‐phyll fluorescence characteristics in wheat leaf[J].Scien‐tia Agricultura Sinica,51(23):4434-4448.]doi:10.3864/j.issn.0578-1752.2018.23.004.

    沈潤平,郭佳,張婧嫻,李洛晞.2017.基于隨機(jī)森林的遙感干旱監(jiān)測模型的構(gòu)建[J].地球信息科學(xué)學(xué)報,19(1):125-133.[Shen R P,Guo J,Zhang J X,Li L X.2017.Construc‐tion of a drought monitoring model using the random fo-rest based remote sensing[J].Journal of Geo-information Science,19(1):125-133.]doi:10.3724/sp.j.1047.2017.00 125.

    王瑞霞,閆長生,張秀英,孫果忠,錢兆國,亓?xí)岳?,牟秋煥,肖世?2018.春季低溫對小麥產(chǎn)量和光合特性的影響[J].作物學(xué)報,44(2):288-296.[Wang R X,Yan C S,Zhang X Y,Sun G Z,Qian Z G,Qi X L,Mu Q H,Xiao S H.2018.Effect of low temperature in spring on yield and photosynthetic characteristics of wheat[J].Acta Agrono-mica Sinica,44(2):288-296.]doi:10.3724/sp.j.1006.2018.00288.

    薛亞光,石呂,魏亞鳳,李波,劉建.2018.稻秸全量還田方式下小麥苗期低溫的生理響應(yīng)及其抗寒性分析[J].南方農(nóng)業(yè)學(xué)報,49(9):1722-1728.[Xue Y G,Shi L,Wei Y F,Li B,Liu J.2018.Physiological responses of wheat under dif‐ferent all rice straw returning patterns to low temperature at seedling stage and cold resistance analysis[J].Journal of Southern Agriculture,49(9):1722-1728.]doi:10.3969/j.issn.2095-1191.2018.09.06.

    楊程,杜思夢,張德奇,時艷華,李向東,邵運(yùn)輝,方保停,王漢芳.2022.基于葉綠素?zé)晒獾男←溤蕉趦龊υu價[J].作物雜志,(1):154-160.[Yang C,Du S M,Zhang D Q,Shi Y H,Li X D,Shao Y H,F(xiàn)ang B T,Wang H F.2022.Evaluation of wheat freezing damage during overwintering period based on chlorophyll fluorescence[J].Crops,(1):154-160.]doi:10.16035/j.issn.1001-7283.2022.01.023.

    楊俊,丁峰,陳晨,劉濤,孫成明,丁大偉,霍中洋.2019.小麥生物量及產(chǎn)量與無人機(jī)圖像特征參數(shù)的相關(guān)性[J].農(nóng)業(yè)工程學(xué)報,35(23):104-110.[Yang J,Ding F,Chen C,Liu T,Sun C M,Ding D W,Huo Z Y.2019.Correlation of wheat biomass and yield with UAV image characteristic parameters[J].Transactions of the Chinese Society of Agri-cultural Engineering,35(23):104-110.]doi:10.11975/j.issn.1002-6819.2019.23.013.

    楊文強(qiáng),林榮呈,端木德強(qiáng),李磊,盧從明,歐陽敏,王文達(dá),孫韜,田利金,王柏臣,朱新廣.2024.近10年光合作用領(lǐng)域若干重要研究進(jìn)展[J].植物生理學(xué)報,60(2):211-247.[Yang W Q,Lin R C,Duanmu D Q,Li L,Lu C M,Ou-yang M,Wang W D,Sun T,Tian L J,Wang B C,Zhu X G.2024.Multiple of important progress on photosynthesis inthe last 10 years[J].Plant Physiology Journal,60(2):211-247.]doi:10.13592/j.cnki.ppj.600007.

    楊熙來,朱榴駿,馮兆忠.2023.臭氧脅迫冬小麥葉片高光譜特征和葉綠素含量估算[J].生態(tài)學(xué)報,43(8):3213-3223.[Yang X L,Zhu L J,F(xiàn)eng Z Z.2023.Hyperspectral charac-teristics and chlorophyll content estimation of winter wheat under ozone stress[J].Chinese Journal of Applied Ecology,43(8):3213-3223.]doi:10.5846/stxb20220120 0196.

    岳俊芹,張素,李向東,邵運(yùn)輝,方保停,葛勝修,王漢芳,張德奇,楊程,時艷華.2021.低溫脅迫對小麥葉綠素?zé)晒鈪?shù)及產(chǎn)量的響應(yīng)[J].麥類作物學(xué)報,41(1):105-110.[Yue J Q,Zhang S,Li X D,Shao Y H,F(xiàn)ang B T,Ge S X,Wang H F,Zhang D Q,Yang C,Shi Y H.2021.Effect of lowtemperature stress on chlorophyll fluorescence parameters and yield of wheat[J].Journal of Triticeae Crops,41(1):105-110.]doi:10.7606/j.issn.1009-1041.2021.01.13.

    張磊,李國領(lǐng),張建周,張德奇.2017.拔節(jié)期低溫脅迫對小麥生理特性和產(chǎn)量的影響[J].河南農(nóng)業(yè)科學(xué),46(11):13-18.[Zhang L,Li G L,Zhang J Z,Zhang D Q.2017.Effects oflow temperature stress on physiological characters and yield of two wheat varieties at jointing stage[J].Journal of Henan Agricultural Sciences,46(11):13-18.]doi:10.15933/j.cnki.1004-3268.2017.11.003.

    周波濤,錢進(jìn).2021.IPCCAR6報告解讀:極端天氣氣候事件變化[J].氣候變化研究進(jìn)展,17(6):713-718.[Zhou B T,Qian J.2021.Changes of weather and climate extremes in the IPCC AR6[J].Climate Change Research,17(6):713-718.]doi:10.12006/j.issn.1673-1719.2021.167.

    Belgiu M,Dr?gu?L.2016.Random forest in remote sensing:A review of applications and future directions[J].ISPRS Journal of Photogrammetry and Remote Sensing,114:24-31.doi:10.1016/j.isprsjprs.2016.01.011.

    Chen Z M,Wang F Z,Zhang P,Ke C D,Zhu Y,Cao W X,Jiang H D.2020.Skewed distribution of leaf color RGB model and application of skewed parameters in leaf color description model[J].Plant Methods,16:1-8.doi:10.1186/s13007-020-0561-2.

    Dong Z F,Men Y H,Liu Z Z,Li J P,Ji J W.2020.Application of chlorophyll fluorescence imaging technique in analysis and detection of chilling injury of tomato seedlings[J].Computers and Electronics in Agriculture,168:105109.doi:10.1016/j.compag.2019.105109.

    Feng M C,Guo X L,Wang C,Yang W D,Shi C C,Ding G W,Zhang X R,Xiao L J,Zhang M J,Song X Y.2018.Moni‐toring and evaluation in freeze stress of winter wheat(Triticum aestivum L.)through canopy hyperspectrumreflectance and multiple statistical analysis[J].Ecological Indicators,84:290-297.doi:10.1016/j.ecolind.2017.08.059.

    Fu X Q,Bai Y,Zhou J,Zhang H W,Xian J Y.2021.A method for obtaining field wheat freezing injury phenotype based on RGB camera and software control[J].Plant Methods,17:1-10.doi:10.1186/s 13007-021-00821-7.

    Goltsev V N,Kalaji H M,Paunov M,B?ba W,Horaczek T,Mojski J,Kociel H,Allakhverdiev S I.2016.Variable chlo‐rophyll fluorescence and its use for assessing physiological condition of plant photosynthetic apparatus[J].Russian Journal of Plant Physiology,63:869-893.doi:10.1134/s1021443716050058.

    Ku K B,Mansoor S,Han G D,Chung Y S,Tuan T T.2023.Identification of new cold tolerant Zoysia grass species using high-resolution RGB and multi-spectral imaging[J].Scientific Reports,13(1):13209.doi:10.1038/s41598-023-40128-2.

    Lee M K,Golzarian M R,Kim I.2021.A new color index forvegetation segmentation and classification[J].PrecisionAgriculture,22:179-204.doi:10.1007/s 11119-020-09735-1.

    Li X N,Cai J,Liu F L,Zhou Q,Dai T B,Cao W X,Jiang D.2015.Wheat plants exposed to winter warming are more susceptible to low temperature stress in the spring[J].Plant Growth Regulation,77:11-19.doi:10.1007/s 10725-015-0029-y.

    Li Z B,Guo R H,Li M,Chen Y R,Li G Y.2020.A review of computer vision technologies for plant phenotyping[J].Computers and Electronics in Agriculture,176:105672.doi:10.1016/j.compag.2020.105672.

    Liu L L,Ji H T,An J P,Shi K J,Ma J F,Liu B,Tang L,Cao W X,Zhu Y.2019.Response of biomass accumulation in wheat to low-temperature stress at jointing and booting stages[J].Environmental and Experimental Botany,157:46-57.doi:10.1016/j.envexpbot.2018.09.026.

    Macedo-Cruz A,Pajares G,Santos M,Villegas-Romero I.2011.Digital image sensor-based assessment of the status of oat(Avena sativa L.)crops after frost damage[J].Sen-sors,11(6):6015-6036.doi:10.3390/s 110606015.

    Mellor A,Boukir S,Haywood A,Jones S.2015.Exploring issues of training data imbalance and mislabelling on ran-dom forest performance for large area land cover classifica-tion using the ensemble margin[J].ISPRS Journal of Pho-togrammetry and Remote Sensing,105:155-168.doi:10.1016/j.isprsjprs.2015.03.014.

    Murphy M E,Boruff B,Callow J N,F(xiàn)lower K C.2020.Detec-ting frost stress in wheat:A controlled environment hyper-spectral study on wheat plant components and implications for multispectral field sensing[J].Remote Sensing,12(3):477.doi:10.3390/rs 12030477.

    Pessarakli M.2024.Handbook of photosynthesis[M].The 4th Edition.Boca Raton:CRC Press:714-722.doi:10.1201/b22922.

    Qiao L,Tang W J,Gao D H,Zhao R M,An L L,Li M Z,Sun H,Song D.2022.UAV-based chlorophyll content estima-tion by evaluating vegetation index responses under diffe-rent crop coverages[J].Computers and Electronics in Agri-culture,196:106775.doi:10.1016/j.compag.2022.106775.

    Shammi S,Diepeveen D,Zander S,Jones M G K,Sohel F.2024.Early frost detection in wheat using machine lear-ning from vertical temperature distributions[J].Computers and Electronics in Agriculture,221:108950.doi:10.1016/j.compag.2024.108950.

    Shammi S,Sohel F,Diepeveen D,Zander S,Jones M G K.2023.A survey of image-based computational learning techniques for frost detection in plants[J].Information Pro-cessingin Agriculture,10(2):164-191.doi:10.1016/j.inpa.2022.02.003.

    Silva-Perez V,Molero G,Serbin S P,Condon A G,Reynolds M P,F(xiàn)urbank R T,Evans J R.2018.Hyperspectral reflec-tance as a tool to measure biochemical and physiological traits in wheat[J].Journal of Experimental Botany,69(3):483-496.doi:10.1093/jxb/erx421.

    Tian H K,Wang T H,Liu Y D,Qiao X,Li Y Z.2020.Com-puter vision technology in agricultural automation—Areview[J].Information Processing in Agriculture,7(1):1-19.doi:10.1016/j.inpa.2019.09.006.

    Wang H F,Huo Z G,Zhou G S,Wu L,F(xiàn)eng H K.2016.Moni-toring and forecasting winter wheat freeze injury and yield from multi-temporal remotely sensed data[J].Intelligent Automationamp;Soft Computing,22(2):255-260.doi:10.1080/10798587.2015.1095475.

    Wang S,Chen J,Rao Y H,Liu L C,Wang W Q,Dong Q.2020.Response of winter wheat to spring frost from a remote sensing perspective:Damage estimation and influential fac-tors[J].ISPRS Journal of Photogrammetry and Remote Sensing,168:221-235.doi:10.1016/j.isprsjprs.2020.08.014.

    Xiao L J,Liu L L,Asseng S,Xia Y,Tang L,Liu B,Cao W X,Zhu Y.2018.Estimating spring frost and its impact on yield across winter wheat in China[J].Agricultural and Forest Meteorology,260-261:154-164.doi:10.1016/j.agrformet.2018.06.006.

    Xie Y K,Wang C,Yang W D,F(xiàn)eng M C,Qiao X X,Song J Y.2020.Canopy hyperspectral characteristics and yield esti-mation of winter wheat(Triticum aestivum)under low temperature injury[J].Scientific Reports,10(1):244.doi:10.1038/s41598-019-57100-8.

    Yao Z S,Shao R,Zain M,Zhao Y Y,Tian T,Wang J L,Zhang D S,Liu T,Song X X,Sun C M.2024.Achieving wheat seedling freezing injury assessment during the seedling stage using Unmanned Ground Vehicle(UGV)and hyper-spectral imaging technology[J].European Journal of Agro-nomy,161:127375.doi:10.1016/j.eja.2024.127375.

    Yuan W A,Choi D.2021.UAV-based heating requirement determination for frost management in apple orchard[J].Remote Sensing,13(2):273.doi:10.3390/rs 13020273.

    Yue Y J,Zhou Y,Wang J A,Ye X Y.2016.Assessing wheat frost risk with the support of GIS:An approach coupling a growing season meteorological index and a hybrid fuzzy neural network model[J].Sustainability,8(12):1308.doi:10.3390/su8121308.

    Zhang J B,Huan H J,Qiu C,Chen Q,Yi C X,Zhang P.2024a.Accurate classification of wheat freeze injury severity from the color information in digital canopy images[J].PLoS One,19(8):e0306649.doi:10.1371/journal.pone.03 06649.

    Zhang K,Yan F F,Liu P.2024b.The application of hyperspec-tral imaging for wheat biotic and abiotic stress analysis:A review[J].Computers and Electronics in Agriculture,221:109008.doi:10.1016/j.compag.2024.109008.

    Zhao L C,Li Q Z,Zhang Y,Wang H Y,Du X.2020.Normali-zed NDVI valley area index(NNVAI)-based framework for quantitative and timely monitoring of winter wheat frost damage on the Huang-Huai-Hai Plain,China[J].Agri-culture,Ecosystems and Environment,292:106793.doi:10.1016/j.agee.2019.106793.

    (責(zé)任編輯:麻小燕)

    欧美zozozo另类| 两性夫妻黄色片| 亚洲成a人片在线一区二区| 久久久精品大字幕| 免费在线观看影片大全网站| 免费在线观看视频国产中文字幕亚洲| 人人妻,人人澡人人爽秒播| 国产精品,欧美在线| а√天堂www在线а√下载| bbb黄色大片| 久久久精品大字幕| 999久久久国产精品视频| 国产av一区二区精品久久| 久久久精品欧美日韩精品| 欧美日韩一级在线毛片| 黄色a级毛片大全视频| 精品不卡国产一区二区三区| 亚洲av美国av| 亚洲精华国产精华精| 夜夜夜夜夜久久久久| 真人做人爱边吃奶动态| 久久欧美精品欧美久久欧美| or卡值多少钱| 国产69精品久久久久777片 | 亚洲乱码一区二区免费版| 亚洲精品国产精品久久久不卡| 欧美成狂野欧美在线观看| 国产人伦9x9x在线观看| 一区二区三区激情视频| 国产欧美日韩精品亚洲av| 国产精华一区二区三区| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 亚洲片人在线观看| 一边摸一边抽搐一进一小说| 国产精品一区二区免费欧美| 999久久久精品免费观看国产| 久久天堂一区二区三区四区| 成人手机av| 国产人伦9x9x在线观看| 少妇粗大呻吟视频| 欧美日韩中文字幕国产精品一区二区三区| 在线观看免费日韩欧美大片| 99精品在免费线老司机午夜| 欧美一区二区精品小视频在线| 长腿黑丝高跟| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 999久久久精品免费观看国产| 亚洲精品色激情综合| 精品久久久久久久末码| 99热这里只有是精品50| 一进一出抽搐gif免费好疼| 国产成年人精品一区二区| 99riav亚洲国产免费| 国产高清videossex| 国产精品1区2区在线观看.| avwww免费| 日韩av在线大香蕉| 中文字幕最新亚洲高清| 午夜视频精品福利| 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| 又紧又爽又黄一区二区| 丰满的人妻完整版| www.www免费av| 久久久久精品国产欧美久久久| 一二三四在线观看免费中文在| 久久亚洲精品不卡| 国产精品久久视频播放| 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 在线视频色国产色| av中文乱码字幕在线| 亚洲专区中文字幕在线| 亚洲成人精品中文字幕电影| 舔av片在线| 国产亚洲欧美98| 欧美乱色亚洲激情| 一a级毛片在线观看| 亚洲七黄色美女视频| 国产视频内射| 成人永久免费在线观看视频| 亚洲成人中文字幕在线播放| 最近最新免费中文字幕在线| 午夜成年电影在线免费观看| 变态另类成人亚洲欧美熟女| 国产精品av视频在线免费观看| 国产精品av久久久久免费| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久成人av| 老司机在亚洲福利影院| 久久精品aⅴ一区二区三区四区| 亚洲性夜色夜夜综合| 欧美日韩福利视频一区二区| 亚洲精品av麻豆狂野| 人妻夜夜爽99麻豆av| 精品乱码久久久久久99久播| 亚洲欧美日韩高清专用| 久久久久久久久免费视频了| 大型黄色视频在线免费观看| 精品国产乱码久久久久久男人| 久久久精品国产亚洲av高清涩受| 日日摸夜夜添夜夜添小说| 少妇粗大呻吟视频| 亚洲人成77777在线视频| 日本一区二区免费在线视频| 亚洲av成人一区二区三| 美女大奶头视频| 欧美日韩亚洲综合一区二区三区_| 黄色片一级片一级黄色片| 午夜久久久久精精品| 欧美精品亚洲一区二区| 国产精品免费一区二区三区在线| 99在线视频只有这里精品首页| 久久午夜综合久久蜜桃| 亚洲电影在线观看av| 黄色女人牲交| tocl精华| 精品久久久久久久久久免费视频| 精品国产亚洲在线| 女生性感内裤真人,穿戴方法视频| 国产熟女午夜一区二区三区| 两个人看的免费小视频| 成人av一区二区三区在线看| 国产精品自产拍在线观看55亚洲| 两人在一起打扑克的视频| 丁香六月欧美| 午夜免费成人在线视频| 大型av网站在线播放| 久久久久国产一级毛片高清牌| 熟妇人妻久久中文字幕3abv| 99re在线观看精品视频| 18美女黄网站色大片免费观看| 国产视频内射| 欧美激情久久久久久爽电影| 国产精品,欧美在线| 亚洲一码二码三码区别大吗| 成熟少妇高潮喷水视频| 精品少妇一区二区三区视频日本电影| 欧美日韩中文字幕国产精品一区二区三区| 亚洲免费av在线视频| 午夜日韩欧美国产| www日本黄色视频网| 两个人视频免费观看高清| 国产单亲对白刺激| 亚洲成人免费电影在线观看| 免费观看人在逋| 久久伊人香网站| 妹子高潮喷水视频| 一夜夜www| 国产人伦9x9x在线观看| 欧美日本亚洲视频在线播放| 亚洲精华国产精华精| 两人在一起打扑克的视频| 香蕉久久夜色| 99re在线观看精品视频| 禁无遮挡网站| 精品久久久久久久久久久久久| 成人高潮视频无遮挡免费网站| 男女午夜视频在线观看| 亚洲九九香蕉| 婷婷丁香在线五月| 国产爱豆传媒在线观看 | www日本在线高清视频| 日本 av在线| 久久精品成人免费网站| 黑人操中国人逼视频| 巨乳人妻的诱惑在线观看| 中国美女看黄片| 亚洲五月婷婷丁香| 日日摸夜夜添夜夜添小说| 国产野战对白在线观看| 亚洲国产精品成人综合色| 亚洲一区高清亚洲精品| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 好男人在线观看高清免费视频| 色综合欧美亚洲国产小说| 日日摸夜夜添夜夜添小说| 国产精品综合久久久久久久免费| 岛国视频午夜一区免费看| tocl精华| 国内精品久久久久久久电影| 午夜久久久久精精品| 亚洲国产精品合色在线| 亚洲中文字幕一区二区三区有码在线看 | 国产精品精品国产色婷婷| 欧美乱妇无乱码| 中出人妻视频一区二区| 日韩欧美在线二视频| 国产激情偷乱视频一区二区| www日本在线高清视频| 精品电影一区二区在线| 午夜a级毛片| 麻豆国产97在线/欧美 | 久久精品国产综合久久久| 国产又色又爽无遮挡免费看| 久久久精品国产亚洲av高清涩受| 精品国产超薄肉色丝袜足j| 99国产综合亚洲精品| 精品久久久久久久久久久久久| 嫩草影视91久久| 日本 欧美在线| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 1024手机看黄色片| 在线观看66精品国产| 色综合站精品国产| 国产高清videossex| 此物有八面人人有两片| 1024视频免费在线观看| 亚洲国产日韩欧美精品在线观看 | 亚洲av成人精品一区久久| 免费搜索国产男女视频| 亚洲色图 男人天堂 中文字幕| 中文在线观看免费www的网站 | 两个人免费观看高清视频| 精品免费久久久久久久清纯| 欧美丝袜亚洲另类 | av在线天堂中文字幕| 久9热在线精品视频| 精品久久久久久久久久久久久| 超碰成人久久| 欧美日本亚洲视频在线播放| videosex国产| 国产精品美女特级片免费视频播放器 | 丝袜人妻中文字幕| 国产单亲对白刺激| 一本综合久久免费| 午夜免费成人在线视频| 亚洲无线在线观看| 久久久精品国产亚洲av高清涩受| 日日摸夜夜添夜夜添小说| 亚洲熟女毛片儿| 国产精品av久久久久免费| 欧美久久黑人一区二区| 成人三级黄色视频| 亚洲色图 男人天堂 中文字幕| 欧美日韩瑟瑟在线播放| 最近视频中文字幕2019在线8| 黄色丝袜av网址大全| 午夜影院日韩av| tocl精华| 久久久久国产精品人妻aⅴ院| 久久久水蜜桃国产精品网| 午夜福利成人在线免费观看| 日本免费一区二区三区高清不卡| 黄片大片在线免费观看| 免费无遮挡裸体视频| 亚洲无线在线观看| 精品国内亚洲2022精品成人| 99热只有精品国产| 高清在线国产一区| 亚洲国产看品久久| 国产69精品久久久久777片 | 听说在线观看完整版免费高清| 国产高清videossex| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 人成视频在线观看免费观看| 日本熟妇午夜| 91av网站免费观看| 日日摸夜夜添夜夜添小说| 国产精品久久久久久精品电影| 国产男靠女视频免费网站| 国产成人影院久久av| 国产免费男女视频| 欧美一区二区国产精品久久精品 | 午夜福利欧美成人| 精品久久蜜臀av无| 久久国产乱子伦精品免费另类| 在线观看日韩欧美| 日本黄大片高清| 欧美黑人巨大hd| 欧美乱码精品一区二区三区| 亚洲激情在线av| 色综合欧美亚洲国产小说| 亚洲全国av大片| 日韩高清综合在线| av福利片在线观看| 最近最新免费中文字幕在线| 手机成人av网站| 舔av片在线| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三| 亚洲一区高清亚洲精品| 欧美日本亚洲视频在线播放| 国产精品美女特级片免费视频播放器 | 日韩 欧美 亚洲 中文字幕| 久久这里只有精品19| 黄色 视频免费看| 又大又爽又粗| 曰老女人黄片| 亚洲七黄色美女视频| 国产亚洲精品综合一区在线观看 | 免费在线观看视频国产中文字幕亚洲| 欧美一区二区国产精品久久精品 | 免费av毛片视频| 精品午夜福利视频在线观看一区| 午夜影院日韩av| 无人区码免费观看不卡| 大型黄色视频在线免费观看| 亚洲成人久久爱视频| 久久久久久免费高清国产稀缺| 亚洲国产欧美一区二区综合| 日韩精品中文字幕看吧| 日日干狠狠操夜夜爽| 国产亚洲精品一区二区www| 国产亚洲欧美在线一区二区| 日日夜夜操网爽| 美女高潮喷水抽搐中文字幕| 黄色丝袜av网址大全| 在线观看一区二区三区| 亚洲男人的天堂狠狠| 美女大奶头视频| e午夜精品久久久久久久| 伊人久久大香线蕉亚洲五| 最近在线观看免费完整版| 一边摸一边做爽爽视频免费| 一级片免费观看大全| 老司机深夜福利视频在线观看| 久久精品综合一区二区三区| 国产激情偷乱视频一区二区| 琪琪午夜伦伦电影理论片6080| 久久久精品大字幕| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| 亚洲,欧美精品.| 正在播放国产对白刺激| 精品国产美女av久久久久小说| 成熟少妇高潮喷水视频| 亚洲精品久久成人aⅴ小说| 首页视频小说图片口味搜索| 欧美av亚洲av综合av国产av| 久久精品国产99精品国产亚洲性色| 午夜免费成人在线视频| 国产区一区二久久| 村上凉子中文字幕在线| 成年女人毛片免费观看观看9| 99国产极品粉嫩在线观看| 三级毛片av免费| 天天躁夜夜躁狠狠躁躁| 99热6这里只有精品| 国产高清视频在线播放一区| 可以在线观看的亚洲视频| 三级国产精品欧美在线观看 | 人妻丰满熟妇av一区二区三区| 91麻豆精品激情在线观看国产| 国内久久婷婷六月综合欲色啪| 日本 av在线| 国产精品久久电影中文字幕| av在线播放免费不卡| 婷婷丁香在线五月| 亚洲精品久久成人aⅴ小说| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 一区二区三区国产精品乱码| 悠悠久久av| 国产精品日韩av在线免费观看| 色综合欧美亚洲国产小说| 12—13女人毛片做爰片一| 99国产综合亚洲精品| 丰满人妻熟妇乱又伦精品不卡| 在线观看免费视频日本深夜| 中文字幕人妻丝袜一区二区| 免费看日本二区| 日本黄大片高清| 免费观看人在逋| 中文字幕久久专区| 国产精品久久久av美女十八| 麻豆成人午夜福利视频| www日本在线高清视频| 亚洲精品国产精品久久久不卡| 美女扒开内裤让男人捅视频| 久久婷婷成人综合色麻豆| 性色av乱码一区二区三区2| 亚洲在线自拍视频| avwww免费| 国内精品久久久久久久电影| 亚洲电影在线观看av| 18禁美女被吸乳视频| 亚洲国产欧美网| 十八禁网站免费在线| 精品一区二区三区av网在线观看| 变态另类丝袜制服| 午夜日韩欧美国产| 搡老妇女老女人老熟妇| 免费无遮挡裸体视频| 最近最新免费中文字幕在线| 精品国产亚洲在线| 久99久视频精品免费| 91字幕亚洲| 国产日本99.免费观看| 色在线成人网| 色av中文字幕| 宅男免费午夜| 最新美女视频免费是黄的| 亚洲精品色激情综合| 脱女人内裤的视频| 亚洲人成77777在线视频| 国产午夜精品论理片| 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 18美女黄网站色大片免费观看| 亚洲熟女毛片儿| 制服诱惑二区| 久久精品人妻少妇| www国产在线视频色| 国产真人三级小视频在线观看| 欧美成人免费av一区二区三区| 日本 av在线| 亚洲天堂国产精品一区在线| 日韩三级视频一区二区三区| 国产熟女xx| 久久精品国产亚洲av高清一级| 老司机在亚洲福利影院| 青草久久国产| 日本免费一区二区三区高清不卡| 欧美黄色淫秽网站| 久久精品aⅴ一区二区三区四区| 午夜久久久久精精品| www.999成人在线观看| 国产1区2区3区精品| 舔av片在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 精品少妇一区二区三区视频日本电影| 久久久国产成人免费| 好男人电影高清在线观看| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| www.精华液| 免费在线观看完整版高清| 国产探花在线观看一区二区| 舔av片在线| 757午夜福利合集在线观看| 18美女黄网站色大片免费观看| 在线观看日韩欧美| 三级国产精品欧美在线观看 | 成人18禁高潮啪啪吃奶动态图| 亚洲熟妇熟女久久| 欧美日本视频| 欧美黑人巨大hd| 精品国产亚洲在线| 日本一区二区免费在线视频| 一级a爱片免费观看的视频| 在线播放国产精品三级| 亚洲精品国产精品久久久不卡| 久久久久免费精品人妻一区二区| 免费人成视频x8x8入口观看| 成人特级黄色片久久久久久久| 操出白浆在线播放| 18禁黄网站禁片免费观看直播| www.熟女人妻精品国产| 九色成人免费人妻av| 国产熟女xx| 男女视频在线观看网站免费 | xxxwww97欧美| 两个人免费观看高清视频| 欧美一级a爱片免费观看看 | 一边摸一边做爽爽视频免费| 女警被强在线播放| 中文字幕人成人乱码亚洲影| 国产精品一区二区免费欧美| 听说在线观看完整版免费高清| 亚洲午夜精品一区,二区,三区| 又大又爽又粗| 男人舔奶头视频| 成熟少妇高潮喷水视频| 国产高清视频在线观看网站| 国产午夜精品论理片| a级毛片在线看网站| 人妻夜夜爽99麻豆av| 两人在一起打扑克的视频| 久久亚洲精品不卡| 成人一区二区视频在线观看| 亚洲欧美一区二区三区黑人| 亚洲 欧美 日韩 在线 免费| 免费在线观看成人毛片| 99久久精品热视频| 国产激情欧美一区二区| 午夜亚洲福利在线播放| 久久精品91蜜桃| 老司机午夜十八禁免费视频| 两个人免费观看高清视频| 脱女人内裤的视频| 亚洲七黄色美女视频| 久久精品91蜜桃| 国产日本99.免费观看| 久久精品aⅴ一区二区三区四区| 一级毛片女人18水好多| 亚洲av第一区精品v没综合| 国产精品久久久久久久电影 | 久久午夜亚洲精品久久| 男女午夜视频在线观看| 一本精品99久久精品77| 天天添夜夜摸| 黑人巨大精品欧美一区二区mp4| 欧美日韩福利视频一区二区| 国产不卡一卡二| 人妻丰满熟妇av一区二区三区| 我要搜黄色片| 动漫黄色视频在线观看| xxxwww97欧美| 美女黄网站色视频| 在线观看舔阴道视频| 亚洲精华国产精华精| 亚洲专区中文字幕在线| 日韩欧美在线二视频| 亚洲精品粉嫩美女一区| 曰老女人黄片| 日本一区二区免费在线视频| 欧美黑人精品巨大| 成人精品一区二区免费| 亚洲精华国产精华精| 久久精品人妻少妇| 色尼玛亚洲综合影院| 国产亚洲欧美在线一区二区| 国产乱人伦免费视频| 亚洲成人久久爱视频| av欧美777| 国产激情偷乱视频一区二区| 母亲3免费完整高清在线观看| 两个人的视频大全免费| 婷婷丁香在线五月| 99在线人妻在线中文字幕| www日本在线高清视频| 最近在线观看免费完整版| 亚洲成人中文字幕在线播放| 免费在线观看日本一区| 69av精品久久久久久| 久久久久久人人人人人| 久久精品国产99精品国产亚洲性色| 精品乱码久久久久久99久播| 又黄又粗又硬又大视频| tocl精华| 亚洲国产看品久久| 天堂av国产一区二区熟女人妻 | 无限看片的www在线观看| 久久性视频一级片| 精品久久久久久,| 久久久久久久久免费视频了| 国产精品自产拍在线观看55亚洲| 欧美绝顶高潮抽搐喷水| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 欧美性猛交黑人性爽| 免费一级毛片在线播放高清视频| 欧美黄色片欧美黄色片| 亚洲欧美日韩高清在线视频| 亚洲熟妇中文字幕五十中出| 99精品在免费线老司机午夜| 欧美性猛交黑人性爽| 香蕉久久夜色| 日本a在线网址| 老司机在亚洲福利影院| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 无人区码免费观看不卡| 亚洲自拍偷在线| 熟女少妇亚洲综合色aaa.| 国产免费av片在线观看野外av| 欧美大码av| 麻豆国产97在线/欧美 | 国产精品久久久av美女十八| 国产久久久一区二区三区| 久久精品国产清高在天天线| 啦啦啦观看免费观看视频高清| videosex国产| 免费在线观看视频国产中文字幕亚洲| 亚洲国产欧美一区二区综合| 999久久久精品免费观看国产| 狂野欧美白嫩少妇大欣赏| 天天躁狠狠躁夜夜躁狠狠躁| 99精品在免费线老司机午夜| 日本 欧美在线| 人人妻人人澡欧美一区二区| 夜夜躁狠狠躁天天躁| 69av精品久久久久久| 国产三级中文精品| 级片在线观看| 免费看美女性在线毛片视频| 免费在线观看影片大全网站| 日本在线视频免费播放| 欧美日韩亚洲综合一区二区三区_| 色哟哟哟哟哟哟| 国产成人啪精品午夜网站| 成在线人永久免费视频| 日本黄大片高清| 欧洲精品卡2卡3卡4卡5卡区| 99国产综合亚洲精品| 黄色 视频免费看| 国产精品永久免费网站| 欧美日韩乱码在线| 日本免费一区二区三区高清不卡| 天堂影院成人在线观看| 51午夜福利影视在线观看| 91字幕亚洲| 男女那种视频在线观看| 狠狠狠狠99中文字幕| 午夜福利高清视频| 高清在线国产一区| 黄频高清免费视频| 亚洲全国av大片| 精品第一国产精品| 日韩av在线大香蕉| 亚洲精品一区av在线观看| 看免费av毛片| 在线观看午夜福利视频| av国产免费在线观看| 亚洲国产欧美网| 麻豆国产av国片精品|