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

    基于多光譜和氣象參數(shù)的菜心水分脅迫指數(shù)反演

    2022-06-27 02:28:28王衛(wèi)星楊明欣謝家興孫道宗曹亞芃駱潤玫藍(lán)于洋
    農(nóng)業(yè)工程學(xué)報 2022年6期
    關(guān)鍵詞:菜心植被指數(shù)冠層

    王衛(wèi)星,楊明欣,高 鵬,謝家興,孫道宗,曹亞芃,駱潤玫,藍(lán)于洋

    基于多光譜和氣象參數(shù)的菜心水分脅迫指數(shù)反演

    王衛(wèi)星1,2,楊明欣1,高 鵬1,謝家興1,2,孫道宗1,2,曹亞芃1,駱潤玫1,藍(lán)于洋1

    (1. 華南農(nóng)業(yè)大學(xué)電子工程學(xué)院(人工智能學(xué)院),廣州 510642;2. 廣東省農(nóng)情信息監(jiān)測工程技術(shù)研究中心,廣州 510642)

    作物水分脅迫指數(shù)(Crop Water Stress Index,CWSI)的監(jiān)測對掌握作物的水分狀況、指導(dǎo)灌溉具有重要意義。該研究以菜心為試驗對象,測量了不同土壤水分條件下的冠層溫度,采集了空氣溫度、相對濕度、風(fēng)速、光合有效輻射和4個波段(450、650、808、940 nm)的光譜反射圖像,并計算了歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)、差值植被指數(shù)(Difference Vegetation Index,DVI)、再歸一化差值植被指數(shù)(Re-Difference Vegetation Index,RDVI)和轉(zhuǎn)換型土壤調(diào)整指數(shù)(Optimized Soil-Adjusted Vegetation Index,OSAVI)等,通過支持向量回歸(Support Vector Regression,SVR)分別構(gòu)建了CWSI上基線、CWSI下基線和冠層溫度的反演模型。結(jié)果表明,菜心在450和650 nm的冠層光譜反射率在0~0.1之間,在808和940 nm的反射率較高,在0.4~0.6之間,當(dāng)菜心由營養(yǎng)生長階段進(jìn)入生殖生長階段,808和940 nm的反射率有所上升。植被指數(shù)能反映菜心的生長狀態(tài)和植被覆蓋度,隨著冠層溫度的升高,NDVI、DVI、RDVI上升,OSAVI下降;而同一個水分處理組在不同生長期的植被指數(shù)有明顯的差異,生殖生長期的植被指數(shù)變化范圍小于營養(yǎng)生長期。結(jié)果表明,使用空氣溫度、相對濕度、風(fēng)速、光合有效輻射反演CWSI上、下基線具有可行性,決定系數(shù)均大于0.75;使用植被指數(shù)反演菜心在兩個生長期的冠層溫度具有較好精度,決定系數(shù)均大于0.7。基于反演值計算的CWSI與基于測量值計算的CWSI有較好的相關(guān)性,決定系數(shù)為0.70;CWSI與氣孔導(dǎo)度是負(fù)相關(guān)的關(guān)系,決定系數(shù)為0.53。該研究應(yīng)用氣象參數(shù)反演CWSI上基線和CWSI下基線,利用植被指數(shù)反演冠層溫度,基于SVR的模型反演值達(dá)到了一定的擬合效果,為實現(xiàn)菜心水分脅迫指數(shù)的光譜監(jiān)測提供支持。

    水分;溫度;菜心;氣象參數(shù);植被指數(shù);支持向量回歸

    0 引 言

    菜心是十字花科蔬菜,其在廣東省內(nèi)分布較廣,是廣東省重要的蔬菜之一。研究表明,當(dāng)菜心長期遭受水分脅迫時,其產(chǎn)量隨之降低[1-3]。楊恒山等[4-5]指出土壤含水率反應(yīng)作物的水分脅迫程度具有局限性,應(yīng)結(jié)合作物本身的生理活動進(jìn)行分析。當(dāng)作物遭受水分脅迫時,其根部無法吸收充足水分,葉片氣孔開始關(guān)閉,蒸騰速率下降,最終表現(xiàn)為葉片溫度上升[6]。隨著紅外技術(shù)的發(fā)展,許多學(xué)者利用水分脅迫指數(shù)(Crop Water Stress Index,CWSI)反映和指示作物的水分脅迫程度。CWSI可實現(xiàn)對蔬菜等易損作物的無損監(jiān)測[7],又能根據(jù)與蒸散速率的關(guān)系建立方程計算灌溉量[8-9]。因此,研究水分脅迫指數(shù)有利于監(jiān)測菜心的需水狀況和制定對應(yīng)的灌溉計劃。CWSI需要測量作物在無水分脅迫脅迫狀態(tài)下的冠層溫度(記為下基線)和嚴(yán)重水分脅迫狀態(tài)下的冠層溫度(記為上基線)。由于作物的冠層溫度受環(huán)境影響,在溫度、濕度、輻射、風(fēng)速等環(huán)境因素變化時,上下基線隨之變化[10-11]。Kumar[12-14]等提出空氣溫度、相對濕度和冠層溫度反演基線的神經(jīng)網(wǎng)絡(luò)模型,均提高了CWSI的應(yīng)用范圍,減少上下基線的測量。

    隨著遙感技術(shù)迅速發(fā)展,光譜技術(shù)在農(nóng)業(yè)方面得到更多應(yīng)用,許多學(xué)者將光譜技術(shù)用于模擬和反演作物的生理參數(shù)[15-18]。研究表明,當(dāng)作物葉片溫度升高,光譜反射率隨之增加,歸一化植被值數(shù)則減少[19-21]。許改平等[22-23]利用光譜影像進(jìn)行高溫脅迫試驗時發(fā)現(xiàn),處于高溫狀態(tài)(40、39 ℃)的毛竹、水稻葉片,其光譜反射率大于溫度低的葉片(25、35 ℃)。Sagan等[24]獲取了700~1 400 nm的光譜影像,將歸一化植被指數(shù)與葉片溫度進(jìn)行擬合分析,發(fā)現(xiàn)歸一化植被指數(shù)與葉片溫度為負(fù)相關(guān)。Ihuoma等[25]以馬鈴薯為研究對象,建立了植被指數(shù)與冠層溫度的線性方程,其相關(guān)系數(shù)接近0.6。梁金晨等[26]獲取了水稻的葉片溫度及該溫度下的光譜數(shù)據(jù),構(gòu)建了多個植被指數(shù)與葉片溫度的多元線性方程,結(jié)果表明方程取得較好的反演效果。同時,許多學(xué)者將機(jī)器學(xué)習(xí)算法應(yīng)用在生理參數(shù)的分類和回歸方面,并獲得很好的效果[27-29]。譚丞軒[30]利用灰度關(guān)聯(lián)分析篩選出在玉米不同生育期下對土壤含水率敏感的植被指數(shù),使用支持向量機(jī)構(gòu)建了玉米土壤含水率的反演模型,相關(guān)系數(shù)均在0.68以上,表明利用機(jī)器學(xué)習(xí)與光譜數(shù)據(jù)相結(jié)合的方法估算土壤含水率是可行的。Li等[31]利用隨機(jī)森林對植被指數(shù)進(jìn)行分析并構(gòu)建模型,實現(xiàn)了竹子地表生物量的反演,模型的精度在0.7以上。

    根據(jù)King等[14]的研究,氣象因子反演CWSI上基線和下基線具有可行性,減少了充分灌溉和嚴(yán)重脅迫處理的試驗,相比于上下基線,氣象因子的采集更簡便;根據(jù)梁金晨等[26]的研究,植被指數(shù)反演冠層溫度具有可行性,相比于紅外測溫儀,多光譜可實現(xiàn)大范圍、快速的數(shù)據(jù)采集。目前,CWSI的研究主要基于氣象因子反演上基線和下基線,而結(jié)合光譜成像技術(shù)的研究仍較少?;谝陨戏治?,本研究以油青菜心為試驗對象,采集空氣溫度、相對濕度、風(fēng)速和光合有效輻射等氣象參數(shù),利用多光譜儀獲取菜心的光譜圖像,計算冠層反射率和植被指數(shù),使用支持回歸向量機(jī)構(gòu)建CWSI上基線、CWSI下基線和冠層溫度的反演模型,以期為菜心水分脅迫指數(shù)的光譜監(jiān)測提供支持。

    1 材料與方法

    1.1 試驗設(shè)計

    試驗于2021年10月26日—12月7日在廣東省廣州市華南農(nóng)業(yè)大學(xué)工程學(xué)院(113°28′06″E,23°12′51″N)室外進(jìn)行。試驗對象為油青菜心(L..()),物候期分別為營養(yǎng)生長期(10月26日—11月25日)和生殖生長期(11月26日—12月07日)。在四葉一心時期,將長勢相近的菜心移到相同的花盆中,開始對其進(jìn)行不同水分處理,當(dāng)植株菜苔高度與苔葉的先端齊平時結(jié)束試驗。花盆內(nèi)徑為40 cm×40 cm,高度為80 cm,土壤經(jīng)烘干、翻松后裝入花盆,每個花盆裝10 kg土壤。

    試驗共設(shè)置4個水分處理組,分別記為T1、T2、T3、T4,每個處理包含3株的油青菜心。T1為充分灌溉處理組,根據(jù)Yang[7]等的研究,充分灌溉條件下土壤最大田間持水量為32.2%,對應(yīng)的灌溉量為600 ml;T2和T3為控制水分處理組,土壤持水量分別為土壤最大田間持水量的70%和40%,對應(yīng)的灌溉量分別為400和200 ml。T4為嚴(yán)重脅迫處理組,作自然處理,試驗開始不進(jìn)行人工補(bǔ)充灌溉,直至試驗結(jié)束。在試驗期間對T1、T2、T3處理組進(jìn)行人工灌溉,灌溉量分別為600、400、200 ml,頻率為兩天1次。

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

    使用土壤水分傳感器(JXBS-J001-EC-RS,精迅,中國)獲取土壤含水率,測量時將該設(shè)備埋在菜心中間,與土壤表面的距離為10 cm。使用紅外測溫儀(T10S-B-HW,妙觀,中國)獲取菜心的葉片溫度,測量時將該設(shè)備固定在菜心頂部,距離約20 cm,保證視角與葉片呈45°。使用光合作用分析儀(SYS-GH30D,塞亞斯,中國)獲取氣孔導(dǎo)度,以5次測量的平均值為該對象的測量值。

    使用多光譜儀(k4 multi-spectrometer,PEAU,美國)獲取光譜影像,該儀器包含4個鏡頭。根據(jù)梁金晨等[26,32]的研究,近紅外波段和紅外波段的光譜反射率與作物葉片溫度相關(guān)性顯著,決定系數(shù)在0.5左右;而歸一化植被指數(shù)NDVI(由紅光波段與近紅外波段的反射率計算)與葉片溫度的決定系數(shù)在0.45。因此,本研究選用了650 nm(紅光波段)、808 nm(近紅外波段)和940 nm(紅外波段),并以450 nm(藍(lán)光波段)作為對比。影像數(shù)據(jù)分辨率為2 048像素×1 536像素,鏡頭焦距為9.60 mm。選擇晴朗無云的時刻進(jìn)行拍攝,拍攝時將配備的標(biāo)準(zhǔn)白板與菜心放在同一平面,多光譜儀置于正上方,固定距離為2 m,儀器通過信號線與計算機(jī)連接,由MAPIR軟件觸發(fā)拍照。

    使用距離試驗地點1 m的室外小型氣象站(WH2081,MISOL,中國)獲取環(huán)境參數(shù)(空氣溫度、相對濕度、風(fēng)速和降雨量),氣象站內(nèi)部包含空氣溫濕度傳感器、風(fēng)速傳感器、雨量傳感器和無線傳輸模塊,數(shù)據(jù)由無限傳輸模塊傳送到計算機(jī)。使用冠層分析儀(MH-G10,綠博,中國)獲取光合有效輻射,測量時將該設(shè)備置于在菜心頂部約20 cm的位置。

    所有的數(shù)據(jù)采集均在12:00—14:00進(jìn)行,氣孔導(dǎo)度和光譜影像的采集周期為每天1次,每次測量重復(fù)5次,其余數(shù)據(jù)的采集周期為10 min/次[13]。試驗期間的天氣情況如表1,各項數(shù)據(jù)的采集示意如圖1所示。在試驗期間,試驗區(qū)域的降雨量較少,降雨天數(shù)共5 d,相比于灌溉量,最大降雨量(僅1 d)對試驗的影響可忽略不計,室外處理可較大程度反映自然環(huán)境下CWSI對菜心水分脅迫的響應(yīng)。

    表1 試驗期間的天氣情況統(tǒng)計(12:00-14:00)

    1.3 水分脅迫指數(shù)計算

    水分脅迫指數(shù)CWSI的計算如下式所示[7]:

    式中Tc為作物冠層溫度,℃;Twet為下基線或無水分脅迫基準(zhǔn)線,是作物在無水分脅迫時或充分灌溉下的冠層溫度,℃;Tdry為上基線或無蒸騰作用基準(zhǔn)線,是作物在無蒸騰作用時、氣孔關(guān)閉狀態(tài)下的冠層溫度,℃。在本研究中,Twet為水分處理組T1的冠層溫度,Tdry為水分處理組T4的冠層溫度,Tc為水分處理組T2、T3的冠層溫度。因此,本研究只計算水分處理組T2、T3的CWSI。

    1.4 光譜數(shù)據(jù)處理

    使用MAPIR軟件對原始反射圖像進(jìn)行預(yù)處理,將MAPIR格式的圖片轉(zhuǎn)換為TIFF格式的灰度圖。圖2顯示4個不同波長的原始光譜圖像。

    圖2 不同波長下的原始光譜圖像

    使用ENVI5.1軟件的ROI工具提取目標(biāo)葉片的灰度值,每個波段的反射率i計算如下[15]:

    式中i為波段(450、650、808、940 nm)的灰度值,d為標(biāo)準(zhǔn)板中黑板的灰度值,w為標(biāo)準(zhǔn)板中白板的灰度值。

    歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)、差值植被指數(shù)(Difference Vegetation Index,DVI)、再歸一化差值植被指數(shù)(Re-Difference Vegetation Index,RDVI)和轉(zhuǎn)換型土壤調(diào)整指數(shù)(Optimized Soil-Adjusted Vegetation Index,OSAVI)的計算公式如下[13,26]:

    式中NIR是近紅外波段的反射率,RED是紅光波段的反射率,在本研究中,NIR、RED分別為808和650。

    1.5 模型構(gòu)建

    SVR最早于1995年提出[33],是一種有監(jiān)督的機(jī)器學(xué)習(xí)算法,在回歸和分類方面得到廣泛應(yīng)用。對于回歸問題,給定= {(1,1), (2,2), ..., (x,y)},SVR希望將每個樣本(x,y)盡量擬合到模型中()T·,=1,2,...,。將()記作最優(yōu)超平面,和分別決定了超平面的方向和位置,通過求到超平面的距離,引入松弛變量和懲罰系數(shù)求損失函數(shù),SVR將最優(yōu)超平面問題轉(zhuǎn)化為最優(yōu)化問題。

    約束條件如公式(6)~(9)所示,最終支持向量回歸模型如公式(11)所示。

    基于氣象因子構(gòu)建冠層溫度基線預(yù)測模型的數(shù)據(jù)集共有430個樣本,訓(xùn)練集有301個樣本,驗證集有129個樣本?;谥脖恢笖?shù)構(gòu)建冠層溫度反演模型的數(shù)據(jù)集共有267個樣本,在營養(yǎng)生長期,訓(xùn)練集有112個樣本,驗證集有48個樣本,在生殖生長期,訓(xùn)練集有75個樣本,驗證集有32個樣本。

    1.6 模型評估

    本研究采用誤差分析評估驗證集模型的精度,誤差分析包括決定系數(shù)(Coefficient of determination,2)、均方根誤差(Root Mean Square Error,RMSE)、平均絕對誤差(Mean Absolute Error,MAE)[13],計算如下所示,2越高,RMSE越小,MAE越小,則模型的預(yù)測精度越高。

    2 結(jié)果與分析

    2.1 冠層溫度的數(shù)據(jù)特征

    在試驗期間獲取了菜心在不同土壤含水率條件下的冠層溫度,4個水分處理組的冠層溫度日變化如圖3所示。由圖3可以看出,T1組的冠層溫度最低,T4組的冠層溫度最高,T2組和T3組的冠層溫度在T1組和T4組之間波動,這表明冠層溫度隨土壤含水量的減少而升高。當(dāng)菜心進(jìn)入生殖生長期,4個處理的冠層溫度保持較為明顯的差異,可能是在生殖生長期間,菜心的生理特征如光合作用速率、蒸騰作用速率變化趨于穩(wěn)定。

    圖3 冠層溫度日變化曲線

    4個水分處理組的冠層溫度分布特征如圖5所示。由圖4可知,T1組的冠層溫度低于34 ℃,T2組和T3組的冠層溫度范圍在22~36 ℃之間,T4組的冠層溫度大部分高于30 ℃。4個水分處理組的冠層溫度波動范圍較大,冠層溫度最大值與最小值的差異平均值約為10 ℃,而數(shù)據(jù)集的變異大有利于測試反演模型的精度[15]。

    圖4 冠層溫度的特征

    2.2 冠層溫度的光譜特征

    對T2和T3水分處理,每個處理3個采樣點的光譜數(shù)據(jù)取平均值,繪制不同生長期的菜心冠層光譜變化曲線,如圖5所示。

    圖5 光譜反射率變化曲線

    從圖5中可以看出,T2和T3在不同生長期的菜心冠層反射率呈現(xiàn)了先略微下降后快速上升的趨勢。同一生長期內(nèi),T3在4個波段的平均值比T2高,表明隨著冠層溫度的升高,葉片光譜反射率上升。隨著菜心成熟,進(jìn)入到生殖生長期后,葉片面積增大,在近紅外波段的反射率明顯上升。表3顯示了在菜心兩個生長期的植被指數(shù)變化。由表3可知,不同生長期的植被指數(shù)有差異,生殖生長期的植被指數(shù)變化范圍小于營養(yǎng)生長期,因為在營養(yǎng)生長期間,菜心的葉片面積不斷增加,生理活動較活躍,在生殖生長期間菜心趨于成熟狀態(tài),葉片面積達(dá)到最大值,表明植被指數(shù)能反映菜心的生長狀態(tài)和植被覆蓋度。

    表3 植被指數(shù)的特征

    注:NDVI為歸一化植被指數(shù);DVI為差值植被指數(shù);RDVI為再歸一化差值植被指數(shù);OSAVI轉(zhuǎn)換型土壤調(diào)整指數(shù)。

    Note: NDVI is normalized difference vegetation index; DVI is difference vegetation index; RDVI is re-difference vegetation index; OSAVI is optimized soil-adjusted vegetation index.

    2.3 基于氣象因子的Twet和Tdry反演

    以空氣溫度、相對濕度、風(fēng)速和光合有效輻射為自變量,使用SVR構(gòu)建wet和dry的反演模型,驗證集預(yù)測值和測量值的對比結(jié)果如圖6所示。結(jié)果表明使用氣象參數(shù)反演wet和dry具有較好精度,2均大于0.75。但wet的反演模型精度高于dry,文獻(xiàn)[7]指出dry與空氣溫度、光合有效輻射的相關(guān)性弱于wet,而且在生殖生長期,隨著氣象參數(shù)的變化,dry的變化趨于平緩,影響了反演的精度。由散點圖可知,SVR反演較高的冠層溫度時表現(xiàn)較差,這是因為該范圍的數(shù)據(jù)量要小于較低冠層溫度的數(shù)據(jù)量[14],后續(xù)的研究應(yīng)增加高溫天氣的數(shù)據(jù)集。

    注:Twet為作物在無水分脅迫時或充分灌溉下的冠層溫度,℃;Tdry為作物在無蒸騰作用、氣孔關(guān)閉狀態(tài)下的冠層溫度,℃。R2為決定系數(shù);RMSE為均方根誤差;MAE為平均絕對誤差。下同。

    2.4 基于植被指數(shù)的冠層溫度反演

    以植被指數(shù)為自變量,使用SVR構(gòu)建T2組和T3組冠層溫度的反演模型,驗證集預(yù)測值和測量值的對比結(jié)果如圖7所示。結(jié)果表明使用植被指數(shù)反演菜心在兩個生長期的冠層溫度具有較好精度,決定系數(shù)均大于0.7,但營養(yǎng)生長期的反演模型精度高于生殖生長期,表明生殖生長期的植被指數(shù)與冠層溫度相關(guān)性低于營養(yǎng)生長期。梁金晨[26]等分別利用RDVI和OSAVI建立水稻葉片溫度的線性反演方程,決定系數(shù)分別為0.524和0.537,低于本研究;RMSE分別為2.495和2.463℃,明顯高于本研究的RMSE。通過誤差分析,模型決定系數(shù)在0.7以上,表明使用植被指數(shù)反演菜心的冠層溫度具有可行性。

    圖7 冠層溫度估算

    2.5 基于反演值與預(yù)測值的CWSI計算

    根據(jù)反演的結(jié)果和公式(1)計算T2組和T3組的CWSI,其與基于測量值的CWSI的關(guān)系如圖8所示。

    基于反演值計算的CWSI精度較高,2為0.70。實現(xiàn)CWSI上基線和下基線的預(yù)測和冠層溫度的反演可以更快速地計算CWSI,縮短診斷灌溉所需的時長,因此兩個模型的結(jié)合對提高CWSI的適用性、提高灌溉效率具有重要意義。

    2.6 CWSI與氣孔導(dǎo)度的關(guān)系

    CWSI與氣孔導(dǎo)度的關(guān)系如圖9所示。T2組的CWSI變化范圍為0.1~0.4,T3組的CWSI變化范圍為0.4~0.8,CWSI越大,水分脅迫程度越深,而T3組的水分脅迫程度要明顯高于T2組。結(jié)果表明,基于反演結(jié)果計算的CWSI與氣孔導(dǎo)度是線性相關(guān)關(guān)系,2為0.53。

    圖8 基于反演值的CWSI與基于測量值的CWSI的對比

    圖9 CWSI與氣孔導(dǎo)度的關(guān)系

    由圖9可知,隨著CWSI增加,菜心的氣孔導(dǎo)度表現(xiàn)為下降,表明CWSI與氣孔導(dǎo)度為負(fù)相關(guān)的關(guān)系。氣孔導(dǎo)度的檢測屬于有損檢測,長期的有損檢測將對菜心產(chǎn)量產(chǎn)生影響,通過探討CWSI隨氣孔導(dǎo)度變化的規(guī)律,減少氣孔導(dǎo)度的有損檢測,并為下一步分析用于觸發(fā)灌溉的CWSI閾值、實現(xiàn)快速精準(zhǔn)灌溉提供參考依據(jù)[34]。

    3 結(jié) 論

    本文針對營養(yǎng)生長期和生殖生長期的菜心,測量了不同土壤含水率條件下的冠層溫度,采集了空氣溫度、相對濕度、風(fēng)速、光合有效輻射和光譜圖像,利用支持向量回歸(Support Vector Regression, SVR)建立了水分協(xié)迫指數(shù)(Crop Water Stress Index,CWSI)上、下基線和冠層溫度的反演模型,研究結(jié)果如下:

    1)同一生長期內(nèi),隨著冠層溫度的升高,葉片光譜反射率上升。隨著菜心成熟,進(jìn)入到生殖生長期后,葉片面積增大,在近紅外、紅外波段的反射率相較于營養(yǎng)生長期明顯上升。

    2)植被指數(shù)能反映菜心的生長狀態(tài)和植被覆蓋度。隨著冠層溫度的升高,歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)、差值植被指數(shù)(Difference Vegetation Index,DVI)、再歸一化差值植被指數(shù)(Re-Difference Vegetation Index,RDVI)上升,轉(zhuǎn)換型土壤調(diào)整指數(shù)(Optimized Soil-Adjusted Vegetation Index,OSAVI)下降;生殖生長期的植被指數(shù)變化范圍小于營養(yǎng)生長期。

    3)使用空氣溫度、相對濕度、風(fēng)速、光合有效輻射反演CWSI上、下基線具有可行性,決定系數(shù)均大于0.75。使用植被指數(shù)反演菜心在兩個生長期的冠層溫度具有較好精度,決定系數(shù)均大于0.7。

    4)基于反演值計算的CWSI具有較高的精度,決定系數(shù)為0.70;CWSI與氣孔導(dǎo)度是負(fù)相關(guān)的關(guān)系,決定系數(shù)為0.53。

    [1] 盧宇鵬,夏巖石,溫少波,等. 不同熟性菜心品質(zhì)性狀的多樣性分析[J]. 廣東農(nóng)業(yè)科學(xué),2020,47(5):161-164.

    Lu Yupeng, Xia Yanshi, Wen Shaobo, et al. Diversity of quality traits of Chinese flowering cabbage varieties with different maturity[J]. Guangdong Agricultural Sciences, 2020, 47(5): 161-164. (in Chinese with English abstract)

    [2] 徐燕. 土壤水分脅迫對菜心生理生化指標(biāo)及氣孔發(fā)育的影響[D]. 廣州:暨南大學(xué),2010.

    Xu Yan. Effects of Soil Water Stress on Physiology and Biochemistry and Stomatal Development in Brassica Chinensis[D]. Guangzhou: Jinan University, 2010. (in Chinese with English abstract)

    [3] 李中赫,占車生,胡實,等. 氣候變化條件下中國灌溉面積變化的產(chǎn)量效應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(19):94-104.

    Li Zhonghe, Zhan Chesheng, Hu Shi, et al. Yield effects of irrigated acreage change under climate change in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(19): 94-104. (in Chinese with English abstract)

    [4] 楊恒山,薛新偉,張瑞富,等. 灌溉方式對西遼河平原玉米產(chǎn)量及水分利用效率的影響[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(21):69-77.

    Yang Hengshan, Xue Xinwei, Zhang Ruifu, et al. Effects of irrigation methods on yield and water use efficiency of maize in the West Liaohe Plain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(21): 69-77. (in Chinese with English abstract)

    [5] Wu X, Shi J C, Zuo Q, et al. Parameterization of the water stress reduction function based on soil-plant water relations[J]. Irrigation Science, 2021, 39: 101-122.

    [6] Mukherjee S, Nandi R, Kundu A, et al. Soil water stress and physiological responses of chickpea (L. ) subject to tillage and irrigation management in lower Gangetic plain[J]. Agricultural Water Management, 2022, 263: 1-17.

    [7] Yang M X, Gao P, Zhou P, et al. Simulating Canopy Temperature using a random forest model to calculate the crop water stress index of Chinese Brassica[J]. Agronomy, 2021, 11: 2244.

    [8] Jamshidi S, Zand-Parsa S, Niyogi D. Assessing crop water stress index of citrus using in-situ measurements, landsat, and sentinel-2 data[J]. International Journal of Remote Sensing, 2021, 42(5): 1893-1916.

    [9] 毋海梅,閆浩芳,張川,等. 溫室滴灌黃瓜產(chǎn)量和水分利用效率對水分脅迫的響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(9):84-93.

    Wu Haimei, Yan Haofang, Zhang Chuan, et al. Responses of yield and water use efficiency of drip-irrigated cucumber in greenhouse to water stress[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(9): 84-93. (in Chinese with English abstract)

    [10] García L, Parra L, Jimenez J M, et al. Deployment strategies of soil monitoring WSN for precision agriculture irrigation scheduling in rural areas[J]. Sensors, 2021, 21(5): 1-27.

    [11] 張芮,王旺田,吳玉霞,等. 水分脅迫度及時期對設(shè)施延遲栽培葡萄耗水和產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(1):155-161.

    Zhang Rui, Wang Wangtian, Wu Yuxia, et al. Effect of moisture stress level and stage on evapotranspiration and yield of grape under protected and delayed cultivation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(1): 155-161. (in Chinese with English abstract)

    [12] Kumar N, Rustum R, Shankar V, et al. Self-organizing map estimator for the crop water stress index[J]. Computers and Electronics in Agriculture. 2021, 187, 106232.

    [13] Kumar N, Adeloye A J, Shankar V, et al. Neural computing modelling of the crop water stress index[J]. Agricultural Water Management. 2020, 239: 106259.

    [14] King B A, Tarkalson D D, Sharma V, et al. Thermal crop water stress index base line temperatures for sugarbeet in arid western U. S[J]. Agricultural Water Management, 2021, 243: 1-12.

    [15] 王偉,彭彥昆,馬偉,等. 冬小麥葉綠素含量高光譜檢測技術(shù)[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2010,41(5):172-177.

    Wang Wei, Peng Yankun, Ma Wei, et al. Prediction of chlorophyll content of winter wheat using leaf-level hyperspectral data[J]. Transactions of the Chinese Society for Agricultural Machinery, 2010, 41(5): 172-177. (in Chinese with English abstract)

    [16] 陳鵬,馮海寬,李長春,等. 無人機(jī)影像光譜和紋理融合信息估算馬鈴薯葉片葉綠素含量[J]. 農(nóng)業(yè)工程學(xué)報,2019,35(11):63-74.

    Chen Peng, Feng Haikuan, Li Changchun, et al. Estimation of chlorophyll content in potato using fusion of texture and spectral features derived from UAV multispectral image[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(11): 63-74. (in Chinese with English abstract)

    [17] 張鑫磊,劉連濤,孫紅春,等. 不同施氮水平下棉花葉片最大羧化速率的高光譜估測[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(11):166-173.

    Zhang Xinlei, Liu Liantao, Sun Hongchun, et al. Hyperspectral estimation of the maximum carboxylation rate of cotton leaves under different nitrogen levels[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(11): 166-173. (in Chinese with English abstract)

    [18] Gokhan C, Kursad D, Levent G. Use of infrared thermography and hyperspectral data to detect effects of water stress on pepper[J]. Quantitative InfraRed Thermography Journal, 2018, 15(1): 81-94.

    [19] Neiff N, Dhliwayo T, Suarez A E, et al. Using an airborne platform to measure canopy temperature and NDVI under heat stress in maize[J]. Journal of Crop Improvement, 2015, 29(6): 669-690.

    [20] Huang S, Tang L N, Joseph P H, et al. A commentary review on the use of normalized difference vegetation index (NDVI) in the era of popular remote sensing[J]. Journal of Forestry Research, 2021, 32: 1-6.

    [21] Mohammed A N, Raj K, Louis L, et al. Using NDVI to differentiate wheat genotypes productivity under dryland and irrigated conditions[J]. Remote sensing, 2020, 12(5): 1-17.

    [22] 許改平,吳興波,劉芳,等. 高溫脅迫下毛竹葉片色素含量與反射光譜的相關(guān)性[J]. 林業(yè)科學(xué),2014,50(5):41-48.

    Xu Gaiping, Wu Xingbo, Liu Fang, et al. The correlation between the pigment content and reflectance spectrum in phyllostachys edulis leaves subjected to high temperature[J]. Scientia Silvae Sinicae, 2014, 50(5): 41-48. (in Chinese with English abstract)

    [23] 程高峰. 高溫脅迫下水稻生理及光譜特性的研究[D]. 鎮(zhèn)江:江蘇大學(xué),2009.

    Cheng Gaofeng. Study on the Physiological and Hyperspectral Characteristics of Rice Under High Temperature Stress[D]. Zhenjiang: Jiangsu University, 2009. (in Chinese with English abstract)

    [24] Sagan V, Maimaitijiang M, Sidike P, et al. UAV-based high resolution thermal imaging for vegetation monitoring, and plant phenotyping using ICI 8640 P, FLIR Vue Pro R 640, and thermomap cameras[J]. Remote Sensing, 2019, 11(3): 1-29.

    [25] Ihuoma S O, Madramootoo C A. Sensitivity of spectral vegetation indices for monitoring water stress in tomato plants[J]. Computers and Electronics in Agriculture, 2019, 163: 1-10.

    [26] 梁金晨,江曉東,楊沈斌,等. 基于高光譜遙感數(shù)據(jù)的水稻葉溫反演[J]. 南方農(nóng)業(yè)學(xué)報,2020,51(1):230-236.

    Liang Jinchen, Jiang Xiaodong, Yang Shenbin, et al. Rice leaf temperature inversion based on hyperspectral remote sensing data[J]. Journal of Southern Agriculture, 2020, 51(1): 230-236. (in Chinese with English abstract)

    [27] Gyamerah S A, Ngare P, Ikpe D. Probabilistic forecasting of crop yields via quantile random forest and Epanechnikov Kernel function[J]. Agriculture and Forest Meteorology, 2020, 280: 107808.

    [28] Virnodkar S S, Pachghare V K, Patil V C, et al. Remote sensing and machine learning for crop water stress determination in various crops: a critical review[J]. Precision Agriculture, 2020, 21(5): 1121-1155.

    [29] 徐敏,趙艷霞,張顧,等. 基于機(jī)器學(xué)習(xí)算法的冬小麥?zhǔn)蓟ㄆ陬A(yù)報方法[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(11):162-171.

    Xu Min, Zhao Yanxia, Zhang Gu, et al. Method for forecasting winter wheat first flowering stage based on machine learning algorithm[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(11): 162-171. (in Chinese with English abstract)

    [30] 譚丞軒. 基于無人機(jī)多光譜遙感的大田玉米土壤含水率估算模型研究[D]. 楊凌:西北農(nóng)林科技大學(xué),2020.

    Tan Chengxuan. Research on the Estimation Model of Field Maize Soil Moisture Content Based on UAV Multispectral Remote Sensing[D]. Yangling: Northwest A&F University, 2020. (in Chinese with English abstract)

    [31] Li X J, Du H Q, Mao F J, et al. Estimating bamboo forest aboveground biomass using EnKF-assimilated MODIS LAI spatiotemporal data and machine learning algorithms[J]. Agriculture and Forest Meteorology, 2018, 256: 445-457.

    [32] 黃春燕,王登偉,肖莉娟,等. 不同水分條件下棉花光譜數(shù)據(jù)對冠層葉片溫度的響應(yīng)特征[J]. 棉花學(xué)報,2014,26(3):244-251.

    Huang Chunyan, Wang Dengwei, Xiao Lijuan, et al. The responsive characteristics between cotton canopy leaves temperature from infrared thermography and hyperspectral data under different water conditions[J]. Cotton Science, 2014, 26(3): 244-251. (in Chinese with English abstract)

    [33] Corinna C, Vladimir V. Support vector networks[J]. Machine Learning, 1995, 20(3):273-297.

    [34] Osroosh Y, Troy Peters R, Campbell C S, et al. Automatic irrigation scheduling of apple trees using theoretical crop water stress index with an innovative dynamic threshold[J]. Computer and Electronic in Agriculture, 2015, 118: 193-203.

    Inverting the water stress index of theusing multiple-spectral and meteorological parameters

    Wang Weixing1,2, Yang Mingxin1, Gao Peng1, Xie Jiaxing1,2, Sun Daozong1,2, Cao Yapeng1, Luo Runmei1, Lan Yuyang1

    (1.(),510642,; 2.,510642,)

    Monitoring the Crop Water Stress Index (CWSI) is of great significance for the water status and irrigation in crop production. Taking theas the test object, this study aims to measure the canopy temperature under different soil moisture conditions. Some meteorological parameters were collected, including the air temperature, relative humidity, wind speed, and photosynthetic active radiation. Meanwhile, the images of spectral reflectance were also collected for the four bands (450, 650, 808, and 940nm). Four vegetation indexes were then calculated by the canopy spectral reflectance, including the Normalized Difference Vegetation Index (NDVI), Difference Vegetation Index (DVI), Re-Difference Vegetation Index (RDVI), and Optimized Soil-Adjusted Vegetation Index (OSAVI). Support Vector Regression (SVR) was selected to construct the inversion models of the CWSI upper/lower baseline using the meteorological parameters, and the inversion models of the canopy temperature using the vegetation index. The results showed that the canopy spectral reflectance at 450 and 650 nm for theranged from 0 to 0.1, while the relatively higher one at 808 and 940 nm ranged from 0.4 to 0.6. The reflectance at 808 and 940 nm increased outstandingly, when thewas developed gradually from the vegetative to reproductive growth stage. The vegetation index reflected the growth state and vegetation coverage of the. There was a different response of vegetation indexes to the canopy temperature. The vegetation NDVI, DVI and RDVI increased, while the vegetation OSAVI decreased with the increase of the canopy temperature of the. The vegetation index under the same water treatment was slightly different in the various growth stages. Specifically, the range of the vegetation index in the reproductive growth stage was smaller than that in the vegetative growth stage. The error analysis showed that the inversion models were feasible to monitor the air temperature, relative humidity, wind speed, and photosynthetic radiation, further invert the upper/lower baseline of CWSI with the determination coefficient greater than 0.75. In the light of the error analysis of the inversion models, the vegetation index was inverted the canopy temperature of thein the vegetative and reproductive growth stage, indicating an excellent accuracy with the determination coefficient greater than 0.7. The calculated CWSI using the inversion models presented a significant correlation with the using the measurement, while the determination coefficient was equal to 0.70. And the CWSI showed the negative relationship with the stomatal conductance with the determination coefficient equal to 0.53. The meteorological parameters were used to invert the upper/lower baseline of CWSI, where the vegetation indexes were used to invert the canopy temperature. The inverted values using the SVR model shared the better fitting performance. The finding can provide a strong support for the spectral monitoring of the crop water stress index of the.

    water; temperature;; meteorological parameter; vegetation index; support vector regression

    10.11975/j.issn.1002-6819.2022.06.018

    S27

    A

    1002-6819(2022)-06-0157-08

    王衛(wèi)星,楊明欣,高鵬,等. 基于多光譜和氣象參數(shù)的菜心水分脅迫指數(shù)反演[J]. 農(nóng)業(yè)工程學(xué)報,2022,38(6):157-164.doi:10.11975/j.issn.1002-6819.2022.06.018 http://www.tcsae.org

    Wang Weixing, Yang Mingxin, Gao Peng, et al. Inverting the water stress index of theusing multiple-spectral and meteorological parameters[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(6): 157-164. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.06.018 http://www.tcsae.org

    2021-12-25

    2022-02-26

    廣東省重點領(lǐng)域研發(fā)計劃(2019B020214003);廣東省鄉(xiāng)村振興戰(zhàn)略專項省級組織實施項目(粵財農(nóng)〔2021〕37 號)“廣東省現(xiàn)代農(nóng)業(yè)關(guān)鍵技術(shù)模式集成與示范推廣”;廣東省科技專項資金(“大專項+任務(wù)清單”)項目(2020020103);廣東省教育廳特色創(chuàng)新類項目(2019KTSCX013);華南農(nóng)業(yè)大學(xué)新農(nóng)村發(fā)展研究院農(nóng)業(yè)科技合作共建項目(2021XNYNYKJHZGJ032)

    王衛(wèi)星,教授,博士生導(dǎo)師,研究方向為農(nóng)情信息獲取與智慧農(nóng)業(yè)。Email:weixin@scau.edu.cn

    猜你喜歡
    菜心植被指數(shù)冠層
    連州菜心秋種忙
    源流(2022年11期)2023-01-29 11:26:20
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達(dá)的樹形靶標(biāo)冠層葉面積探測模型研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    最好的老年人食譜——牛肝菌扒菜心
    一些菜心,一些故事
    西湖(2016年9期)2016-05-14 08:44:03
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    精品一区二区免费观看| 赤兔流量卡办理| 国产成人精品福利久久| 亚洲天堂av无毛| 亚洲内射少妇av| 最近手机中文字幕大全| 涩涩av久久男人的天堂| 夫妻性生交免费视频一级片| 久久久久精品久久久久真实原创| 日韩成人av中文字幕在线观看| 久久久久网色| 久久久久久久亚洲中文字幕| 大话2 男鬼变身卡| 大码成人一级视频| 在线 av 中文字幕| 日韩一区二区视频免费看| 亚洲国产高清在线一区二区三| 一区二区av电影网| 国产高清国产精品国产三级 | 99热这里只有是精品50| 日本黄大片高清| 中文天堂在线官网| 国产伦精品一区二区三区视频9| 国产精品人妻久久久影院| 亚洲精品国产色婷婷电影| 老司机影院成人| 国产大屁股一区二区在线视频| 人人妻人人澡人人爽人人夜夜| 少妇丰满av| 最近中文字幕高清免费大全6| 亚洲精品第二区| 亚洲丝袜综合中文字幕| 午夜日本视频在线| 偷拍熟女少妇极品色| 韩国av在线不卡| 成人午夜精彩视频在线观看| 亚洲精品亚洲一区二区| 久久综合国产亚洲精品| 在线亚洲精品国产二区图片欧美 | 一本久久精品| 九九爱精品视频在线观看| 欧美日本视频| 纵有疾风起免费观看全集完整版| 国产一区二区亚洲精品在线观看| 国产精品99久久99久久久不卡 | 一区二区三区四区激情视频| 亚洲美女搞黄在线观看| 九九在线视频观看精品| 少妇的逼好多水| 国产片特级美女逼逼视频| 在线观看一区二区三区激情| 国产乱人视频| a级毛片免费高清观看在线播放| 亚洲成人久久爱视频| 成人一区二区视频在线观看| 国产极品天堂在线| 亚洲色图综合在线观看| 亚洲精品乱码久久久v下载方式| 99久久九九国产精品国产免费| 26uuu在线亚洲综合色| 啦啦啦中文免费视频观看日本| 美女视频免费永久观看网站| 欧美激情久久久久久爽电影| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 天堂网av新在线| 免费观看在线日韩| 国产亚洲最大av| 最近的中文字幕免费完整| 免费看不卡的av| 神马国产精品三级电影在线观看| 中国国产av一级| 建设人人有责人人尽责人人享有的 | 成年免费大片在线观看| 欧美人与善性xxx| 国产精品无大码| 黄色怎么调成土黄色| 亚洲精品中文字幕在线视频 | 国产男女内射视频| 国产 精品1| 午夜视频国产福利| 国产 一区精品| 日日摸夜夜添夜夜爱| 中文字幕免费在线视频6| 伊人久久精品亚洲午夜| 视频区图区小说| 性插视频无遮挡在线免费观看| 久久久久久久久久久丰满| 美女脱内裤让男人舔精品视频| 国产精品国产三级国产av玫瑰| 免费观看的影片在线观看| 在线观看一区二区三区激情| 日韩一本色道免费dvd| av一本久久久久| 日韩一区二区三区影片| 国产探花极品一区二区| 最近最新中文字幕免费大全7| 毛片女人毛片| 91精品一卡2卡3卡4卡| 免费看a级黄色片| 超碰av人人做人人爽久久| 免费看日本二区| 国产成人精品福利久久| 午夜激情久久久久久久| 啦啦啦啦在线视频资源| 99精国产麻豆久久婷婷| 夫妻午夜视频| 成人二区视频| 久久女婷五月综合色啪小说 | 中文字幕制服av| 久久久久久久久久久免费av| 亚洲内射少妇av| 人妻系列 视频| 搞女人的毛片| 中国三级夫妇交换| 国产乱来视频区| 国产成人freesex在线| 日韩精品有码人妻一区| 久久久久网色| www.av在线官网国产| 99久久精品热视频| 国产美女午夜福利| 大香蕉久久网| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 亚洲经典国产精华液单| av又黄又爽大尺度在线免费看| 丝袜脚勾引网站| 免费看不卡的av| 亚洲精品成人久久久久久| 性插视频无遮挡在线免费观看| 国产色爽女视频免费观看| 亚洲精品aⅴ在线观看| 亚洲精品国产色婷婷电影| 嘟嘟电影网在线观看| 亚洲av二区三区四区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | .国产精品久久| 亚洲精品色激情综合| 欧美xxxx性猛交bbbb| 国产成人一区二区在线| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 麻豆成人av视频| 亚洲av成人精品一二三区| 一本久久精品| 国产探花极品一区二区| 国产黄a三级三级三级人| 韩国av在线不卡| 大香蕉久久网| 亚洲精品国产成人久久av| 国产精品蜜桃在线观看| 99视频精品全部免费 在线| 2018国产大陆天天弄谢| 欧美高清性xxxxhd video| 亚洲国产高清在线一区二区三| 精品一区在线观看国产| 晚上一个人看的免费电影| 全区人妻精品视频| 在线观看美女被高潮喷水网站| 丝瓜视频免费看黄片| 国产午夜福利久久久久久| a级毛片免费高清观看在线播放| 如何舔出高潮| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 欧美高清性xxxxhd video| 亚洲av日韩在线播放| 搡老乐熟女国产| 欧美极品一区二区三区四区| 日韩av不卡免费在线播放| 女人十人毛片免费观看3o分钟| 观看免费一级毛片| 男人舔奶头视频| 久久97久久精品| 国产高清三级在线| 3wmmmm亚洲av在线观看| 国产精品嫩草影院av在线观看| 久久亚洲国产成人精品v| 国产伦精品一区二区三区视频9| 亚洲精品日本国产第一区| 欧美+日韩+精品| 日韩大片免费观看网站| 18禁动态无遮挡网站| 你懂的网址亚洲精品在线观看| 久久久精品欧美日韩精品| 亚洲人成网站在线播| 国内精品美女久久久久久| 亚洲av一区综合| 最近最新中文字幕大全电影3| 2022亚洲国产成人精品| 99久久精品国产国产毛片| 亚洲精华国产精华液的使用体验| av免费在线看不卡| 男人狂女人下面高潮的视频| 成人毛片60女人毛片免费| 国产综合精华液| 中文字幕人妻熟人妻熟丝袜美| 精品人妻偷拍中文字幕| 精品久久久精品久久久| 欧美高清成人免费视频www| 精品国产三级普通话版| 99热全是精品| 免费观看在线日韩| 欧美激情在线99| 成人综合一区亚洲| 99久久精品热视频| 蜜桃亚洲精品一区二区三区| 国产永久视频网站| 国语对白做爰xxxⅹ性视频网站| 国产在线男女| 亚洲国产av新网站| 国产精品国产三级国产专区5o| 网址你懂的国产日韩在线| 久久久久久久午夜电影| 国产 一区 欧美 日韩| 日产精品乱码卡一卡2卡三| 美女内射精品一级片tv| 不卡视频在线观看欧美| 丰满乱子伦码专区| 国产免费又黄又爽又色| 亚洲天堂av无毛| 国产老妇伦熟女老妇高清| 国产探花极品一区二区| 国产高清国产精品国产三级 | 永久网站在线| 色播亚洲综合网| 免费看光身美女| 久久久国产一区二区| 香蕉精品网在线| 国内精品美女久久久久久| 日本爱情动作片www.在线观看| 欧美97在线视频| 日韩中字成人| 青春草国产在线视频| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 国产 一区 欧美 日韩| 欧美日韩视频精品一区| 国产中年淑女户外野战色| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 国产午夜精品一二区理论片| 自拍欧美九色日韩亚洲蝌蚪91 | 韩国高清视频一区二区三区| 精品人妻偷拍中文字幕| 国产亚洲av嫩草精品影院| 十八禁网站网址无遮挡 | 亚洲综合色惰| 欧美日韩亚洲高清精品| 亚洲精品,欧美精品| 亚洲精品国产av成人精品| 国产毛片在线视频| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| av免费观看日本| 日韩 亚洲 欧美在线| 99精国产麻豆久久婷婷| 午夜亚洲福利在线播放| 国产又色又爽无遮挡免| 国产精品三级大全| 最近最新中文字幕大全电影3| 哪个播放器可以免费观看大片| 中文乱码字字幕精品一区二区三区| 亚洲色图av天堂| 久久ye,这里只有精品| 日本三级黄在线观看| 全区人妻精品视频| 国产成年人精品一区二区| 插阴视频在线观看视频| 伊人久久国产一区二区| 两个人的视频大全免费| 99久久精品国产国产毛片| 日本熟妇午夜| 一级二级三级毛片免费看| 在线观看一区二区三区激情| 成人国产av品久久久| 哪个播放器可以免费观看大片| 国产成人aa在线观看| 熟女电影av网| 一级毛片我不卡| 国产一区二区三区综合在线观看 | 麻豆乱淫一区二区| 国内少妇人妻偷人精品xxx网站| 欧美日韩国产mv在线观看视频 | 久久久久精品久久久久真实原创| 嘟嘟电影网在线观看| 亚洲精品视频女| 免费观看的影片在线观看| 黄色一级大片看看| 免费观看性生交大片5| 久久久欧美国产精品| 国产成人精品婷婷| 午夜精品一区二区三区免费看| 精品国产三级普通话版| 国产91av在线免费观看| 26uuu在线亚洲综合色| 久久精品久久久久久久性| 边亲边吃奶的免费视频| 中文乱码字字幕精品一区二区三区| 国产精品一区二区在线观看99| 在线观看美女被高潮喷水网站| 80岁老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 亚洲av在线观看美女高潮| 精品人妻视频免费看| 亚洲人与动物交配视频| 日韩欧美 国产精品| 一区二区三区乱码不卡18| 欧美激情久久久久久爽电影| 久久6这里有精品| 亚洲无线观看免费| 青春草亚洲视频在线观看| 久久这里有精品视频免费| 国产一区亚洲一区在线观看| 亚洲人成网站高清观看| 国产成人免费观看mmmm| 极品少妇高潮喷水抽搐| 少妇人妻久久综合中文| 精品亚洲乱码少妇综合久久| 美女国产视频在线观看| 中文在线观看免费www的网站| 又大又黄又爽视频免费| 蜜桃久久精品国产亚洲av| 日韩欧美精品v在线| 少妇人妻精品综合一区二区| 神马国产精品三级电影在线观看| 韩国av在线不卡| 在线观看av片永久免费下载| 狂野欧美白嫩少妇大欣赏| 国产免费一级a男人的天堂| 男插女下体视频免费在线播放| 波野结衣二区三区在线| 成人无遮挡网站| 国产黄片美女视频| 大片免费播放器 马上看| 日本熟妇午夜| 69av精品久久久久久| 伦精品一区二区三区| 听说在线观看完整版免费高清| 赤兔流量卡办理| 国产精品三级大全| 亚洲人成网站在线观看播放| 综合色av麻豆| 欧美高清成人免费视频www| 国产男女内射视频| 精品久久久久久久人妻蜜臀av| 啦啦啦在线观看免费高清www| 国产精品麻豆人妻色哟哟久久| 国产精品久久久久久精品电影| 内地一区二区视频在线| 搡女人真爽免费视频火全软件| 国产精品久久久久久久久免| 精华霜和精华液先用哪个| 亚洲精品456在线播放app| 深夜a级毛片| 热re99久久精品国产66热6| 国产极品天堂在线| 性色av一级| 久久久午夜欧美精品| 日韩 亚洲 欧美在线| 久久久久久九九精品二区国产| 亚洲欧美中文字幕日韩二区| 22中文网久久字幕| 久久精品国产亚洲网站| 国产真实伦视频高清在线观看| 韩国高清视频一区二区三区| 久久国产乱子免费精品| 99热这里只有精品一区| 亚洲国产精品专区欧美| 中文字幕人妻熟人妻熟丝袜美| 日日撸夜夜添| 精品久久久噜噜| 亚洲精品乱久久久久久| 大又大粗又爽又黄少妇毛片口| 五月玫瑰六月丁香| av专区在线播放| 国国产精品蜜臀av免费| 色5月婷婷丁香| 欧美一区二区亚洲| 国产伦在线观看视频一区| 91久久精品国产一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 久热久热在线精品观看| 久久精品久久久久久噜噜老黄| 热re99久久精品国产66热6| 日本-黄色视频高清免费观看| 国产精品久久久久久av不卡| 国产精品久久久久久久电影| 在现免费观看毛片| 男女啪啪激烈高潮av片| 三级国产精品片| av在线观看视频网站免费| 日本黄色片子视频| 精华霜和精华液先用哪个| 中文字幕亚洲精品专区| 久久久久久伊人网av| 真实男女啪啪啪动态图| 一级片'在线观看视频| 嫩草影院入口| 欧美日韩综合久久久久久| 久久久久九九精品影院| 日韩欧美 国产精品| 欧美xxxx性猛交bbbb| 国产极品天堂在线| 日本黄色片子视频| 久久精品久久久久久噜噜老黄| 日本三级黄在线观看| 国产精品一及| 亚洲丝袜综合中文字幕| 亚洲av电影在线观看一区二区三区 | 日本wwww免费看| 中文字幕人妻熟人妻熟丝袜美| 黑人高潮一二区| 成年女人在线观看亚洲视频 | 三级经典国产精品| 美女内射精品一级片tv| 欧美xxⅹ黑人| 亚洲国产精品成人久久小说| freevideosex欧美| 美女视频免费永久观看网站| 精品一区二区免费观看| 精品一区二区三卡| 国产精品爽爽va在线观看网站| 国产伦在线观看视频一区| 一区二区三区四区激情视频| 精品国产露脸久久av麻豆| 国产精品嫩草影院av在线观看| 久久热精品热| 国产精品久久久久久精品古装| 26uuu在线亚洲综合色| 精品亚洲乱码少妇综合久久| 一级片'在线观看视频| 午夜免费男女啪啪视频观看| 91在线精品国自产拍蜜月| 波野结衣二区三区在线| 精品一区二区三区视频在线| 精品人妻偷拍中文字幕| 亚洲精品456在线播放app| 日本wwww免费看| 亚洲av成人精品一区久久| 国产一区二区亚洲精品在线观看| 亚洲在久久综合| 少妇猛男粗大的猛烈进出视频 | 国产乱来视频区| 亚洲成人中文字幕在线播放| 乱码一卡2卡4卡精品| 亚洲av福利一区| 亚洲av中文字字幕乱码综合| 综合色丁香网| 欧美日韩综合久久久久久| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区| 亚洲精品国产av蜜桃| 亚洲av男天堂| 日本熟妇午夜| av在线播放精品| 人人妻人人澡人人爽人人夜夜| 国产亚洲av嫩草精品影院| 久久久亚洲精品成人影院| 麻豆精品久久久久久蜜桃| 亚洲熟女精品中文字幕| 国产日韩欧美亚洲二区| 一级毛片电影观看| 欧美激情久久久久久爽电影| 2021少妇久久久久久久久久久| 久久精品人妻少妇| 中文字幕制服av| 一本色道久久久久久精品综合| 秋霞伦理黄片| 久久97久久精品| 欧美日韩一区二区视频在线观看视频在线 | 丝袜喷水一区| 男女国产视频网站| 大香蕉久久网| 国产精品嫩草影院av在线观看| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 大码成人一级视频| 麻豆成人av视频| 国模一区二区三区四区视频| 精品人妻偷拍中文字幕| 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 亚洲精华国产精华液的使用体验| 国产成人a∨麻豆精品| 亚洲av日韩在线播放| a级毛色黄片| 亚洲美女视频黄频| 夜夜爽夜夜爽视频| 美女xxoo啪啪120秒动态图| 久久女婷五月综合色啪小说 | 插逼视频在线观看| 久久人人爽人人片av| 久久精品久久久久久噜噜老黄| 国产成人a区在线观看| 中国国产av一级| 国产精品一区二区在线观看99| 亚洲国产色片| 亚洲熟女精品中文字幕| 久久久久久久亚洲中文字幕| 欧美一区二区亚洲| 熟女电影av网| a级毛色黄片| 少妇人妻久久综合中文| 另类亚洲欧美激情| 大片免费播放器 马上看| av网站免费在线观看视频| 特大巨黑吊av在线直播| 啦啦啦在线观看免费高清www| 欧美人与善性xxx| 久久久久九九精品影院| 亚洲精品,欧美精品| 97在线人人人人妻| a级毛片免费高清观看在线播放| 少妇高潮的动态图| 亚洲精品日本国产第一区| 亚洲欧美日韩无卡精品| 久久久a久久爽久久v久久| 青春草亚洲视频在线观看| 久久久久久久久久人人人人人人| 亚洲欧美一区二区三区国产| 日韩不卡一区二区三区视频在线| 久热这里只有精品99| 欧美zozozo另类| 欧美精品一区二区大全| 亚洲人成网站在线播| 精品久久久久久久久av| 99久久九九国产精品国产免费| av女优亚洲男人天堂| 3wmmmm亚洲av在线观看| 久久精品熟女亚洲av麻豆精品| 久久精品久久精品一区二区三区| 我的老师免费观看完整版| 女人被狂操c到高潮| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 18禁裸乳无遮挡动漫免费视频 | 能在线免费看毛片的网站| 日韩人妻高清精品专区| 成年av动漫网址| 国产黄a三级三级三级人| 国产精品一区二区三区四区免费观看| 国产老妇伦熟女老妇高清| 五月伊人婷婷丁香| 成人欧美大片| 18禁在线无遮挡免费观看视频| 亚洲一区二区三区欧美精品 | 丰满乱子伦码专区| .国产精品久久| 99热全是精品| 亚洲精品乱久久久久久| 十八禁网站网址无遮挡 | 国模一区二区三区四区视频| 亚洲精品久久久久久婷婷小说| 极品少妇高潮喷水抽搐| 偷拍熟女少妇极品色| 天堂中文最新版在线下载 | 亚洲精品日韩在线中文字幕| 中文字幕人妻熟人妻熟丝袜美| 深爱激情五月婷婷| 国产免费福利视频在线观看| 特大巨黑吊av在线直播| 男女边吃奶边做爰视频| 国产黄片美女视频| 人妻系列 视频| 亚州av有码| 午夜免费男女啪啪视频观看| 久久99热这里只频精品6学生| 三级经典国产精品| 日韩精品有码人妻一区| 毛片一级片免费看久久久久| av福利片在线观看| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 国产视频内射| 一区二区av电影网| 日韩电影二区| 成人漫画全彩无遮挡| 看免费成人av毛片| 一区二区三区免费毛片| 全区人妻精品视频| 日韩 亚洲 欧美在线| 精品人妻视频免费看| 欧美3d第一页| 日韩伦理黄色片| 国产色爽女视频免费观看| 五月玫瑰六月丁香| 亚洲欧美成人精品一区二区| av国产精品久久久久影院| 一个人看的www免费观看视频| 成年人午夜在线观看视频| 亚洲av欧美aⅴ国产| 秋霞在线观看毛片| 国产精品成人在线| 女人被狂操c到高潮| 免费播放大片免费观看视频在线观看| 最近最新中文字幕大全电影3| 亚洲精品国产色婷婷电影| 成人亚洲精品一区在线观看 | 搡女人真爽免费视频火全软件| 国内精品宾馆在线| 欧美精品人与动牲交sv欧美| 日本黄色片子视频| 免费观看a级毛片全部| 成年av动漫网址| 亚洲av免费在线观看| 国产精品不卡视频一区二区| 丝袜美腿在线中文| 一级毛片黄色毛片免费观看视频| 麻豆精品久久久久久蜜桃| 国产精品久久久久久精品电影| 精品国产三级普通话版| 国产精品久久久久久久久免| 免费电影在线观看免费观看|