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

    基于多種算法的果樹果實生長模型研究—以云南昭通蘋果為例

    2021-09-18 06:22:54孫擎趙艷霞程晉昕曾廳余張祎
    中國農(nóng)業(yè)科學 2021年17期
    關鍵詞:莖干直徑蘋果

    孫擎,趙艷霞,程晉昕,曾廳余,張祎

    基于多種算法的果樹果實生長模型研究—以云南昭通蘋果為例

    1中國氣象科學研究院,北京 100081;2云南省氣候中心,昆明 650000;3云南昭通農(nóng)業(yè)氣象試驗站,云南昭通 657000

    【】針對果樹果實與生長過程中的氣象因子關聯(lián)密切,且生長過程多為非線性、非平穩(wěn)序列,直接對其連續(xù)測定難度較大的問題,對比多種模型對果實直徑的模擬能力,為果樹及其果實的生長發(fā)育監(jiān)測和預測、適時灌溉施肥、生長環(huán)境調控等提供科學參考。以云南昭通蘋果為例,分析2019和2020年果實生長期間直徑變化特征及其與環(huán)境氣候因子的關系。引入深度學習中的長短期記憶模型(LSTM),使用LSTM模型對蘋果果實直徑進行模擬及預測,與多元線性回歸模型(MLR)和機器學習模型中的決策樹(DT)及隨機森林(RF)模型的模擬結果進行對比分析,并使用3種采樣方法對不同模型模擬的結果進行評估。蘋果果實直徑有明顯日變化特征,呈夜間直徑增長而白天縮小為主的規(guī)律,一般早晨直徑達到最大,然后逐漸微縮,在日落前后直徑到達當日最小。蘋果果實直徑的增長速率在果實膨大初期較高,在果實生長后期降低。蘋果果實小時和日平均直徑與土壤溫度和土壤濕度呈中度或高度正相關,與紫外線指數(shù)(UVI)呈高度負相關。蘋果果實直徑的日平均增長量(FMDG)、日增長量(FDG)、日最大變化量(MDFS)與60 cm土壤溫度和20 cm、40 cm土壤濕度呈低負相關(-0.5≤<-0.3)。4個模型的模擬結果相比,LSTM模型的模擬精度高于MLR、DT和RF模型。LSTM模型比MLR模型在相關系數(shù)增加3%—20%的情況下,RMSE和MAE下降約50%—75%,而機器學習模型DT和RF對蘋果果實直徑的預測相對較差,可能存在過度擬合。對比統(tǒng)計學、機器學習和深度學習等方法,LSTM模型在蘋果果實直徑的模擬中表現(xiàn)出更高的精度和可靠性,能更好地解決果實生長過程中的復雜非線性問題。

    蘋果直徑;生長模型;氣象因子;深度學習;長短期記憶模型(LSTM)

    0 引言

    【研究意義】云南昭通市是云南蘋果種植的主要區(qū)域,也是當?shù)剞r(nóng)民主要的經(jīng)濟收入來源之一。2019年,昭通蘋果種植規(guī)模達4.51萬hm2,總產(chǎn)量65萬 t,總產(chǎn)值60億元[1]。云南地形復雜,呈立體氣候特征,氣象災害種類多、發(fā)生頻率高、危害強度大[2],對蘋果產(chǎn)量造成較大影響,亟需在模擬果實生長規(guī)律的同時,建立氣候因子對蘋果果實生長的影響模型。隨著物聯(lián)網(wǎng)在農(nóng)業(yè)上的廣泛應用,積累了大量的觀測數(shù)據(jù),為嘗試探索應用新的技術和方法來模擬果實生長變化帶來了機會。本研究嘗試提出模擬蘋果果實直徑變化的一種新方法,以彌補目前研究方法的不足。研究結果可精確模擬和預測蘋果果實直徑,對適時進行果樹災害防御、灌溉施肥、生長環(huán)境調控等提供重要參考和應用價值。【前人研究進展】果樹果實直徑每日膨脹和收縮量是果樹生長的重要指標,與作物生長環(huán)境如溫度、日照、土壤水分等環(huán)境因子密切相關[3]。大部分關于果樹生長的研究主要針對不同灌溉量[4]、施肥量[5]對果實生長發(fā)育、品質、直徑等的影響,或莖干微變的研究[6-7]。對蘋果果實生長的研究來說,主要是關于不同施肥量對蘋果果樹生長及果實產(chǎn)量的影響[8],或不同干旱程度對蘋果果實、莖干直徑和葉片水勢的影響[9],徑流與環(huán)境因子的關系[10]。傳統(tǒng)的果實生長發(fā)育模型包括作物生理模型、多元線性回歸(multi-linear regression,MLR)、移動平均模型(moving average model,MA)、自動回歸模型(autoregressive moving average model,ARMA),或基于機器學習方法如BP神經(jīng)網(wǎng)絡、人工神經(jīng)網(wǎng)絡(artificial neural network,ANN)、支持向量機(support vector machines,SVM)等。例如,Li等[11]使用作物生理模型模擬了蘋果果實的生長發(fā)育情況;夏桂敏等[10]建立了蘋果果樹徑流與環(huán)境因子的多元回歸模型;張海輝等[12]使用機器學習方法對蘋果霉心透射光譜進行研究,結果表明機器學習方法在修正果實直徑對透射光譜的影響方面表現(xiàn)出良好的精度;張彪等[13]使用BP神經(jīng)網(wǎng)絡評價了蘋果制干的適宜性。但傳統(tǒng)作物生理模型、統(tǒng)計模型或機器學習方法存在一定局限,如模型過于簡單導致模擬精度較低,或將動態(tài)時間建模問題轉化為靜態(tài)空間建模研究,忽略了過去時間的輸入對于預測的影響[14]。部分研究使用了深度學習中的遞歸神經(jīng)網(wǎng)絡(recurrent neural network,RNN),表現(xiàn)出良好的精度[15],但RNN在訓練過程中存在梯度消失和梯度爆炸等缺陷。長短期記憶模型(long short-term memory,LSTM)作為一種特殊的RNN模型可以有效規(guī)避上述問題,特別是LSTM在時間序列預測方面的優(yōu)良特性,促使越來越多的研究者將LSTM模型應用在各個方面,如面部識別[16]、質量控制[17]、文本識別[18]、病蟲害識別[19]?!颈狙芯壳腥朦c】對于果樹生長發(fā)育中果實或莖干直徑微變來說,僅有的研究主要涉及的是作物莖干微變;同時,目前關于深度學習應用在果樹方面的研究很少,且主要應用在機器視覺方面[20-21],對與經(jīng)濟價值直接關聯(lián)的果實直徑的監(jiān)測及模擬研究鮮見報道?!緮M解決的關鍵問題】本研究首先分析蘋果果實直徑在生長過程中的變化趨勢,進而分析蘋果果實直徑與生長過程中環(huán)境因子的相關性,最后使用多元回歸模型、機器學習模型中的決策樹和隨機森林以及深度學習中的LSTM模型分別模擬昭通蘋果果實直徑在生長過程中的變化情況,并對模擬及預測結果進行對比評估,以驗證模型在果實直徑監(jiān)測及模擬中的精度和可靠性。

    1 材料與方法

    1.1 試驗概況

    云南省昭通市夏秋季節(jié)受西南海洋季風控制,空氣濕度大、溫度高、雨熱同季;冬、春主要受極地大陸季風(干冷氣團)控制,空氣干燥、溫度低、干冷同季。年平均氣溫11.7℃,月平均氣溫≥10℃的天數(shù)為191.3 d,≥10℃的活動積溫3 217.4℃,全年降雨量750—1 100 mm,無霜期221 d,年日照1 902.2 h。蘋果開花期溫度為15—23℃,花芽分化期溫度為20—25℃,果實發(fā)育期溫度25℃左右,云南昭通自然條件適宜蘋果生長,肉質脆且風味濃。

    試驗于2019和2020年5—8月在昭通蘋果果實生長期開展。蘋果品種為‘嘎啦’,在2016年樹齡2年時移栽,果樹種植行距4 m,間距0.5 m,屬于自根砧矮化密植,種植地點為云南省昭通市昭陽區(qū)永豐鎮(zhèn)昭通超越有限公司蘋果種植基地(27.25°N,103.69°E),總種植面積約5萬畝,選用長勢一致、無病蟲害的一棵蘋果果樹及其上正常生長的一個蘋果果實進行觀測(圖1)。2019年2月14日蘋果果樹進入花芽膨大期,3月4日花蕾出現(xiàn),開花期為3月20日—4月3日,5月29日在蘋果果實直徑大于3 cm后開始監(jiān)測蘋果果實直徑,2019年8月16日收獲,其中7月1—4日,由于儀器故障,數(shù)據(jù)未能正常采集。2020年2月24日進入花芽膨大期,3月12日花蕾出現(xiàn),開花期為3月22日—4月8日,在5月18日在蘋果果實直徑大于3 cm后開始監(jiān)測蘋果果實直徑,8月20日收獲,其中7月15—16日,由于儀器故障,數(shù)據(jù)未能正常采集。依據(jù)果園常規(guī)管理方式開展施肥、灌溉和農(nóng)藥管理。

    圖1 蘋果果實和莖干直徑監(jiān)測儀示意圖

    1.2 觀測指標和方法

    本試驗采用摩爾多瓦Bio Instruments公司生產(chǎn)的線性位移傳感器連續(xù)自動監(jiān)測蘋果樹果實和莖干直徑(型號FI-LM),傳感器與物聯(lián)網(wǎng)自動監(jiān)測系統(tǒng)Caipos(北京科百宏業(yè)科技有限公司生產(chǎn))連接,每間隔1 h自動記錄作物果實直徑的變化量(μm)。果園內(nèi)物聯(lián)網(wǎng)監(jiān)測系統(tǒng)可以采集以下環(huán)境氣象要素:降水(mm)、氣壓(mbar)、太陽輻射(W·m-2)、空氣溫度(℃)、空氣相對濕度(%)、風速(m?s-1)、風向(度)、紫外線指數(shù)(UVI)、葉面溫度(℃)、土壤濕度(%)、土壤溫度(℃)等。其中土壤濕度采用英國Delta-T公司C200A型土壤水分傳感器,可測定20、40和60 cm的土壤濕度(體積含水量)和土壤溫度。2019年5月29日至8月16日與2020年5月18日至8月20日蘋果果實直徑監(jiān)測期間觀測的主要氣象數(shù)據(jù)如圖2所示。

    基于果實直徑傳感器測量的果實直徑小時數(shù)據(jù)(μm),本研究采用以下4種果實直徑變化指標來衡量果實直徑的變化程度:直徑日增長量(FDG)、日平均增長量(FMDG)、日最大變化量(MDFS)和小時增量(FHG)。直徑日增長量(FDG)是當日果實直徑最大值與前一天果實直徑最大值的差值;小時增量(FHG)是當前小時的果實直徑與前一小時的果實直徑的差值;日平均增長量(FMDG)是當日果實的平均直徑與前一日果實平均直徑的差值;日最大變化量(MDFS)是果實直徑日最大值與最小值之差。以上日平均值均為0—23時的平均值。

    1.3 模型簡介

    選取統(tǒng)計學、深度學習和機器學習3種方法中常用的4個模型來模擬蘋果果實直徑,包括多元回歸模型(MLR)、深度學習中的長短期記憶模型(LSTM)、機器學習中的決策樹(decision trees,DT)和隨機森林(random forests,RF)模型。

    1.3.1 多元回歸模型 多元回歸模型的公式如下所示:

    其中,y是因變量,即蘋果果實直徑,1,2,3….是自變量,即輸入的環(huán)境氣象因子,i為各個因子的系數(shù),為誤差項。本研究將使用Python 3.6的基于TensorFlow的Keras框架對多元回歸模型參數(shù)i和進行最優(yōu)擬合,擬合方法為隨機梯度下降法。

    1.3.2 LSTM模型 在深度學習中,RNN常用來處理時間序列數(shù)據(jù),但是RNN對長序列處理會出現(xiàn)記憶丟失的問題,從而導致梯度消失和梯度爆炸,這使得RNN模型無法深入訓練。長短期記憶(long short-term memory,LSTM)是在克服RNN這一缺點基礎上形成的一種RNN變形結構,對時間上有依賴的時序數(shù)據(jù)有著優(yōu)秀的預測能力。通過在LSTM內(nèi)部結構中引入門控機制,通過“門”(gate)來控制丟棄或者保留信息,使得時間序列上的記憶信息可控,具有記憶性,使不同時間點之間的數(shù)據(jù)網(wǎng)絡存在連接[22]。RNN只有一個傳輸狀態(tài)(h)相比,LSTM中有兩個傳輸狀態(tài):c(cell state)和h(hidden state),如圖3-a所示。

    圖3-b為LSTM模型的基本結構圖,其中c為LSTM的記憶狀態(tài)(即記憶神經(jīng)元),又稱為長期記憶,h為短期記憶。LSTM模型的網(wǎng)絡更新為:設為LSTM單元輸出,為LSTM模型的記憶模塊值為單元輸入,WW分別為輸入數(shù)據(jù)和上一時刻的輸出權重系數(shù),則時間的記憶單元值為:

    式中,為c的逐點乘積。輸入門i可控制當前數(shù)據(jù)輸入對記憶模塊狀態(tài)的影響,公式為:

    i=σ(Wx+Wh-1+Wc-1+b) (3)

    遺忘門t可控制歷史學習對當前記憶模塊狀態(tài)值的影響,決定單元狀態(tài)被遺忘的信息公式為:

    f=σ(Wx+Wh-1+Wc-1+b) (4)

    當前時刻記憶單元狀態(tài)則為:

    輸出門t可控制記憶模塊狀態(tài)制的輸出,公式為:

    o=σ(Wx+Wh-1+Wc-1+b) (6)

    LSTM的輸出單元h為:

    上述公式(2)—(7)的表示相關的權重向量,表示相關的偏差向量。

    本研究將使用Python3.6的基于TensorFlow的Keras深度學習框架,進行LSTM模型的擬合與訓練,利用dropout函數(shù)防止過度擬合,預測的數(shù)據(jù)為蘋果的果實直徑,訓練LSTM模型時采用Adam算法動態(tài)調整每個輸入?yún)?shù)的學習速率,使用MSE最小化為損失函數(shù)的優(yōu)化目標。

    1.3.3 決策樹模型 決策樹(DT)是一種逼近離散值目標函數(shù)的方法,DT學習到的函數(shù)被表示為一棵決策樹,代表的是對象屬性與對象值之間的一種映射關系。決策樹計算復雜度不高,輸出結果易于理解,對中間值缺失不敏感,可以處理不相關特征數(shù)據(jù)[23]。

    圖3 長短期記憶網(wǎng)絡LSTM模型結構圖

    1.3.4 隨機森林模型 隨機森林(RF)是由CART算法衍生出的分類與回歸算法[24]。RF是一種集成算法(ensemble learning algorithm),使用樣本數(shù)據(jù)中的不同子集建立多個分類回歸樹,其中每一個子集都在樣本數(shù)據(jù)中隨機選擇,且每個子集間相互獨立。每一個樹都是一個獨立的分類與回歸樹,樹的每一個節(jié)點內(nèi)可以建立一個簡單、準確的模型來模擬自變量與因變量之間的關系。由于該算法對每個節(jié)點進行拆分時運用的方案隨機,并從這些隨機方案中選出最優(yōu)點,所以隨機森林算法對異常值和噪聲有較高的容忍度[25]。

    1.4 模型評價指標

    果實直徑模擬結果評估使用的指標為:均方根誤差(root mean square error,RMSE)、相關系數(shù)(correlation coefficient,)和均方誤差(mean square error, MSE)。其中,一般認為0.9<≤1.0為極高相關,0.7<≤0.9為高度相關,0.5<≤0.7為中度相關,0.3<≤0.5為低相關,0.0<≤0.3為幾乎不相關[26-28]。

    2 結果

    2.1 蘋果果實直徑變化趨勢

    從2019年5月29日—8月16日(圖4-a)與2020年5月18日—8月20日(圖4-b)的蘋果果實直徑小時值和日平均值圖可以看出,2019和2020年蘋果果實直徑約從30 000 μm逐漸增大,成熟收獲時約為70000 μm。蘋果果實直徑具有明顯的日變化規(guī)律,在早晨06:00—08:00達到最大,然后果實直徑變化開始微縮,在傍晚日落前后18:00—20:00達到最小。蘋果果徑的日變化在果實生長前期波動較大,在生長后期波動較小。

    從圖5可以看出,蘋果的果徑在果實膨大初期增長較快,在整個6月都處在較高的增長階段,大約每日增長在600—800 μm。在果實生長中期以后,果徑的日增量逐漸減少,在果實生長后期接近成熟時果徑的日增量約200—400 μm。蘋果果徑的日增量和小時增量都有明顯的周期性波動。對于蘋果果徑小時增量FHG來說,從日落前后到第二天早晨,蘋果果徑處在增長的階段,在白天光合作用和蒸騰作用之后,蘋果果徑處于縮小階段,F(xiàn)HG為負值。蘋果果徑的MDFS、FMDG和FDG指標呈較為一致的變化趨勢,周期性地增大或減少,總體呈現(xiàn)出果實生長初期到后期逐漸減小的趨勢。

    圖4 蘋果直徑累積增量2019年(a)和2020年(b)小時變化和日平均變化

    圖5 蘋果直徑2019年(a)和2020年(b)小時增量(FHG)、日平均增長量(FMDG)、日增長量(FDG)、日最大變化量(MDFS)

    2.2 蘋果果實直徑與氣象因子的關系

    圖6為蘋果果實直徑指標(包括小時果徑、日平均果徑、FHG、MDFS、FMDG和FDG)與氣象及生長因子(包括降水,氣壓,氣溫,空氣相對濕度,風速,風向,UVI,太陽輻射,(20、40和60)cm土壤溫濕度,葉片溫度和莖干直徑等)的相關系數(shù)熱力圖,為2019年和2020年相關系數(shù)的平均值。圖中顏色偏藍說明相關系數(shù)呈正相關,顏色偏紅表明呈負相關,其中<0.05。從圖6可以看出,蘋果果實直徑的小時和日平均值與其他因子的相關性較為一致,與莖干直徑呈極高正相關(=0.97),與土壤溫度、土壤濕度、莖干直徑呈中等或高度正相關,與UVI呈高度負相關,與MDFS、DFG和FMDG呈中度負相關,與其他因子相關性較低。MDFS、FDG和FMDG與其他因子的相關性較為一致,與UVI呈中度正相關,與風向、60 cm土壤溫度、20和40 cm土壤濕度呈低負相關,與蘋果日平均直徑、莖干日平均直徑呈中度負相關。FHG與其他因子無明顯的相關性。

    圖6 蘋果果實直徑指標與生長氣候因子相關系數(shù)圖

    2.3 模型模擬結果

    使用MLR、LSTM、DT和RF 4種模型模擬蘋果直徑的生長變化時,考慮到更好的應用目的,使用自動氣象站觀測氣象要素的小時數(shù)據(jù)作為輸入變量(包括降水、氣壓、太陽輻射、空氣溫度、空氣相對濕度、風速、風向、紫外線指數(shù)、土壤濕度、土壤溫度等)。同時,訓練數(shù)據(jù)使用以下3種采樣方法對比模型在不同采樣方法下的適用性:1)將2019和2020年的數(shù)據(jù)按照時間順序放在一起,采樣方法為使用前3 003個數(shù)據(jù)進行訓練,剩余1 001個數(shù)據(jù)作為驗證;2)使用2019年數(shù)據(jù)作為訓練數(shù)據(jù),2020年數(shù)據(jù)作為驗證;3)將2019和2020年的數(shù)據(jù)按照時間順序放在一起,使用前3 003個數(shù)據(jù)中的80%進行隨機采樣訓練,剩余20%作為驗證,最后剩余的1 001個數(shù)據(jù)作為模型預測的驗證。

    從圖7中可以看出,LSTM模型驗證結果在4種模型中最好,模擬值與觀測值達到極高相關(3種采樣方法均大于0.9),其次是MLR模型,3種采樣方法下的觀測值與模擬值的相關系數(shù)在0.7以上,屬于高度相關;DT和RF模型訓練結果較差;所有相關系數(shù)的均小于0.01。使用不同采樣方法進行訓練和驗證的效果也不同,采樣方法1的驗證結果中LSTM模型的最高(>0.9),達到極高度相關,RMSE約為MLR模型的50%,約為DT和RF模型驗證結果的1/3,DT和RF模型的模擬結果相對較差(<0.3),RMSE在9 000 μm左右。采樣方法2使用了2020年的數(shù)據(jù)進行驗證,可以看出MLR和LSTM模型可以較好模擬出蘋果果實的變化規(guī)律,均達到極高度相關,LSTM模型驗證結果的RMSE約為MLR模型驗證結果的50%,為1 757.4 μm;DT和RF模型模擬結果相對較差,可以模擬出果實直徑隨時間上升的趨勢,但是RMSE和MAE較大(>15 000 μm),偏離1﹕1線較遠。對采樣方法3的驗證結果來說,進行隨機采樣的MLR、LSTM、DT和RF模型的驗證結果較好,大部分點在1﹕1線附近,模型驗證結果的均達到極高度相關,除MLR的RMSE為3 106.02 μm,其余模型的RMSE小于1 000 μm,MAE小于300 μm。從采樣方法3的預測結果可以看出LSTM模型模擬的精度最高,大于0.99,RMSE和MAE均在1 000 μm左右;MLR模型的RMSE和MAE在5 000 μm左右,達到高度相關;而DT和RF模型預測的驗證結果較差,小于0.3,RMSE在9 000 μm左右。

    圖8為使用3種采樣方法訓練的模型預測和實測值對比的時間序列圖。整體來看,MLR和LSTM能夠預測出蘋果果實直徑隨時間增大的趨勢,DT和RF模型預測的結果較差。MLR模型在果實生長中期模擬效果較好,在前期和后期模擬結果稍差,在果實生長后期的觀測值與模擬值的最大差值約為13 000 μm。LSTM模型的模擬結果最好,在2020年果實生長前期的預測結果好于MLR模型,使用采樣方法3的預測結果好于采樣方法1的結果,更接近觀測值,在果實生長后期的觀測值與模擬值的最大差值小于2 000 μm,使用采樣方法1和2在果實生長后期模擬結果相近,觀測值與模擬值的最大差值約為8 000 μm。DT和RF模型使用3種采樣方法的預測結果類似,只能預測部分時間的果實直徑,模擬效果較差。

    3 討論

    蘋果果實直徑有明顯的日變化特征,在早晨約6:00—8:00達到最大,日落前后約18:00—20:00達到當日最小,與蘋果樹莖干直徑的日變化較為一致,在LIU等[29]的研究結果中也得到證實。夜間果樹吸水導致果實和莖干的直徑膨大,日間果樹蒸騰致失水增加[30]進而導致果實和莖干直徑收縮,同時日間果樹形成的光合產(chǎn)物在晚上分配轉移至果實和莖干,是引起蘋果果實和莖干日變化的主要原因。蘋果果實直徑的FMDG、FDG和MDFS的波動較為一致,這些指標能夠反映植物體內(nèi)的水分狀況[6],蘋果果實直徑指標與土壤濕度和土壤溫度相關系數(shù)較高也從側面印證了果實中水分的變化。

    不同品種的蘋果樹對氣候環(huán)境因子的響應也不同[31],定量確定不同品種蘋果樹對氣候因子的響應是今后的研究方向。有研究表明果實產(chǎn)量與莖干直徑年增長量成反比,莖干直徑年增長量越低,蘋果果實產(chǎn)量就越高[32],本研究只使用了2年觀測數(shù)據(jù),在未來長時間序列的研究中可以進一步探究果實產(chǎn)量與莖干直徑增長的關系。果實和莖干直徑的大小除了與氣象條件相關外,與土壤養(yǎng)分、樹體營養(yǎng)、施肥管理措施等也密切相關,但由于較難測定或無法定量確定,因此,本研究僅考慮便于監(jiān)測的氣象因子作為果實監(jiān)測及模擬的依據(jù)。

    總體來看,所有采樣方法以LSTM模型模擬結果最好,MLR模型模擬結果較為穩(wěn)健,使用3種采樣方法對果實生長后期果實直徑的模擬結果相近,但均有一定程度低估。使用LSTM模型在果實生長后期預測的結果比采樣方法1和2的結果稍好,更接近觀測值。采樣方法2中MLR和LSTM模型對蘋果果實生長中期的模擬效果較好,采樣方法1和3對果實生長后期的果實直徑有一定程度的低估,這可能是由于使用了2年的果實生長前期和中期數(shù)據(jù)進行訓練,使用1年的果實生長后期數(shù)據(jù)進行訓練所導致。LSTM模型在果實生長前期和后期的模擬效果高于其他模型,說明LSTM模型可以自動學習出果實在不同生長階段的直徑與環(huán)境氣象因子的關系。有研究表明LSTM模型的RMSE比使用生理生長模型[10](physiological development time,PDT模型)會降低50%以上。

    在機器學習和深度學習訓練過程中可能會存在過度擬合現(xiàn)象,導致預測結果低于模型訓練結果。在使用采樣方法3時,使用隨機數(shù)據(jù)進行驗證時,RF和DT模型模擬結果較好,但在使用剩余1 001個蘋果果實直徑進行預測驗證時預測結果較差,對蘋果果實直徑的預測精度遠低于MLR和LSTM模型,這可能主要是由于過度擬合導致。本研究在LSTM模型訓練時采用dropout函數(shù)避免過度擬合,訓練和測試的損失曲線在經(jīng)過迭代后最終重合并未分散,未發(fā)生過度擬合或欠擬合。因此,在4種模型中,LSTM模型對蘋果果實直徑的模擬效果最好,精度最高,且能夠有效防止過度擬合的出現(xiàn)。

    本研究成果可用于無田間果實直徑監(jiān)測儀器情況下的果實直徑監(jiān)測及預測,是一種新的模擬思路和方法,一定程度上解決了傳統(tǒng)統(tǒng)計學方法如多元回歸以及機器學習方法模擬精度不高的問題,深度學習方法表現(xiàn)出更強的穩(wěn)定性和可靠性,在果實直徑監(jiān)測預測方面有廣闊的應用前景[33-35]。

    4 結論

    本文分析了蘋果果實在生長期間直徑的變化情況以及與氣候環(huán)境因子的相關性,使用統(tǒng)計模型MLR、深度學習模型LSTM和機器學習模型DT及RF共4種模型,建立了蘋果果實直徑模擬模型并進行驗證。結果表明,蘋果果實直徑有明顯日變化特征,一般在早晨直徑達到最大,在日落前后直徑為當日最小。蘋果果實直徑的FMDG、FDG和MDFS呈現(xiàn)出果實生長初期變化量增大,生長后期變化量減小的趨勢。蘋果果實直徑和莖干直徑呈極高度正相關,變化趨勢較為一致。蘋果果實直徑的小時和日平均值與土壤濕度和土壤溫度呈中高程度正相關,與40和60 cm土壤溫度及20 cm土壤濕度呈高度正相關,與UVI呈高度負相關。果實直徑的MDFS、FDG、FMDG與60 cm土壤溫度、20及40 cm土壤濕度呈低負相關。在蘋果果實直徑模擬方面,LSTM模型比MLR、DT和RF模型表現(xiàn)出更高的可靠性和精度。

    [1] 王順富, 張丹, 夏敏, 王成軼, 劉毅, 鄧芳, 楊艷群, 郎平勇, 全勇, 李云國. 昭通蘋果產(chǎn)業(yè)助推脫貧攻堅的經(jīng)驗及成效. 云南農(nóng)業(yè)大學學報(社會科學版), 2020, 14(6): 36-41.

    WANG S F, ZHANG D, XIA M, WANG C Y, LIU Y, DENG F, YANG Y Q, LANG P Y, QUAN Y, LI Y G. The experience and effectiveness of apple industry in promoting poverty alleviation program in Zhaotong. Journal of Yunnan Agricultural University (Social Science), 2020, 14(6): 36-41. (in Chinese)

    [2] 周桂華, 吳惠標. 2016年云南主要自然災害特點及減災對策研究. 災害學, 2017, 32(3): 112-119.

    ZHOU G H, WU H B. 2016 Yunnan major natural disasters and the characteristics of the mitigation measures research. Journal of Catastrophology, 2017, 32(3): 112-119. (in Chinese)

    [3] CONEJERO W, ALARCóN J J, GARCíA-ORELLANA Y, ABRISQUETA J M, TORRECILLAS A. Daily sap flow and maximum daily trunk shrinkage measurements for diagnosing water stress in early maturing peach trees during the post-harvest period. Tree Physiology, 2007, 27(1): 81-88.

    [4] 陳毓瑾, 陳立, 婁玉穗, 秦澤冠, 董肖, 馬超, 苗玉彬, 張才喜, 許文平, 王世平. 基于新梢直徑和果粒投影面積變化確定‘巨峰’葡萄果實發(fā)育期的灌溉閾值. 果樹學報, 2019, 36(5): 612-620.

    CHEN Y J, CHEN L, LOU Y S, QIN Z G, DONG X, MA C, MIAO Y B, ZHANG C X, XU W P, WANG S P. Determination of thresholds to trigger irrigation of ‘Kyoho’ grapevine during berry development periods based on variations of shoot diameter and berry projected area. Journal of Fruit Science, 2019, 36(5): 612-620. (in Chinese)

    [5] 王國棟, 肖元松, 彭福田, 張亞飛, 郜懷峰, 孫希武, 賀月. 尿素配施不同用量納米碳對桃幼樹生長及氮素吸收利用的影響. 中國農(nóng)業(yè)科學, 2018, 51(24): 4700-4709.

    WANG G D, XIAO Y S, PENG F T, ZHANG Y F, GAO H F, SUN X W, HE Y. Effects of urea application combined with different amounts of nano-carbon on plant growth along with nitrogen absorption and utilization in young peach trees. Scientia Agricultura Sinica, 2018, 51(24): 4700-4709. (in Chinese)

    [6] 王曉森, 孟兆江, 段愛旺, 劉祖貴. 基于莖直徑變化監(jiān)測番茄水分狀況的機理與方法. 農(nóng)業(yè)工程學報, 2010, 26(12): 107-113.

    WANG X S, MENG Z J, DUAN A W, LIU Z G. Mechanism and methodology for monitoring tomato water content based on stem diameter variation. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(12): 107-113. (in Chinese)

    [7] 林琭, 湯昀, 閆萬麗, 張紀濤, 籍增順. 基于莖直徑變化的精準灌溉技術研究進展. 農(nóng)業(yè)工程學報, 2016, 32(5): 1-13.

    LIN L, TANG Y, YAN W L, ZHANG J T, JI Z S. Overview of precise irrigation scheduling based on stem diameter variations. Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(5): 1-13. (in Chinese)

    [8] 賀琦琦, 郭向紅, 楊凱, 石美娟, 馬娟娟, 孫西歡, 雷濤, 劉艷武. 滴灌灌施磷鉀肥對矮砧蘋果樹生理生長和產(chǎn)量的影響. 節(jié)水灌溉, 2019(8): 24-27.

    HE Q Q, GUO X H, YANG K, SHI M J, MA J J, SUN X H, LEI T, LIU Y W. Effects of phosphorus and potassium fertilizer on physiological growth and yield of dwarf rootstock apple trees under drip irrigation. Water Saving Irrigation, 2019(8): 24-27. (in Chinese)

    [9] 高木旺, 聶佩顯, 孫瀅, 姜遠茂, 張金春, 張繼祥. 避雨環(huán)境下蘋果幼樹水分狀態(tài)指標對干旱脅迫的響應. 應用生態(tài)學報, 2016, 27(8): 2459-2466.

    GAO M W, NIE P X, SUN Y, JIANG Y M, ZHANG J C, ZHANG J X. Response of water status indicators in apple saplings to drought stress under shelter from rain. Chinese Journal of Applied Ecology, 2016, 27(8): 2459-2466. (in Chinese)

    [10] 夏桂敏, 孫媛媛, 王瑋志, 吳奇, 遲道才. ‘寒富’蘋果樹莖流特征及其對環(huán)境因子的響應. 中國農(nóng)業(yè)科學, 2019, 52(4): 701-714.

    XIA G M, SUN Y Y, WANG W Z, WU Q, CHI D C. The characteristics of sap flow of Hanfu apple trees and its response to environmental factors. Scientia Agricultura Sinica, 2019, 52(4): 701-714. (in Chinese)

    [11] LI M, CHEN M X, ZHANG Y, FU C X, XING B, LI W Y, QIAN J P, LI S, WANG H, FAN X D, YAN Y J, WANG Y A, YANG X T. Apple fruit diameter and length estimation by using the thermal and sunshine hours approach and its application to the digital orchard management information system. PLoS ONE, 2015, 10(4): e0120124.

    [12] 張海輝, 田世杰, 馬敏娟, 趙娟, 張軍華, 張佐經(jīng). 考慮直徑影響的蘋果霉心病透射光譜修正及檢測. 農(nóng)業(yè)機械學報, 2019, 50(1): 313-320.

    ZHANG H H, TIAN S J, MA M J, ZHAO J, ZHANG J H, ZHANG Z J. Detection method of moldy core in apples using modified transmission spectrum based on size of fruit. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(1): 313-320. (in Chinese)

    [13] 張彪, 劉璇, 畢金峰, 吳昕燁, 金鑫, 李旋, 李瀟. 基于BP人工神經(jīng)網(wǎng)絡算法的蘋果制干適宜性評價. 中國農(nóng)業(yè)科學, 2019, 52(1): 129-142.

    ZHANG B, LIU X, BI J F, WU X Y, JIN X, LI X, LI X. Suitability evaluation of apple for chips-processing based on BP artificial neural network. Scientia Agricultura Sinica, 2019, 52(1): 129-142. (in Chinese)

    [14] 陳舒蘊. 基于LSTM多類型機組風電場等值建模研究. 山東農(nóng)業(yè)大學學報(自然科學版), 2020, 51(2): 294-297.

    CHEN S Y. Study on the equivalent modeling of wind power farm based on LSTM multi-type wind turbine generators. Journal of Shandong Agricultural University (Natural Science Edition), 2020, 51(2): 294-297. (in Chinese)

    [15] 王子文. 基于深度學習和模糊決策的蘋果自動分揀方法研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2019.

    WANG Z W. Use of deep learning and fuzzy decision making in automated sorting of apple [D]. Harbin: Harbin Institute of Technology, 2019. (in Chinese)

    [16] SUN X D, WU P C, HOI S C H. Face detection using deep learning: an improved faster RCNN approach. Neurocomputing, 2018, 299: 42-50.

    [17] XUE W D, CHAI Y, LI Q G, HONG Y Q, ZHENG G F. Short-term relay quality prediction algorithm based on long and short-term memory. Instrumentation, 2018, 5(4): 46-54.

    [18] LI P, LI X, PAN H, KHYAM O, NOOR-A-RAHIM M. Text-based indoor place recognition with deep neural network. Neurocomputing, 2020, 390: 239-247.

    [19] 沈利言, 姜海燕, 胡濱, 謝元澄. 水稻病蟲草害與藥劑實體關系聯(lián)合抽取算法. 南京農(nóng)業(yè)大學學報, 2020, 43(6): 1151-1161.

    SHEN L Y, JIANG H Y, HU B, XIE Y C. A study on joint entity recognition and relation extraction for rice diseases pests weeds and drugs. Journal of Nanjing Agricultural University, 2020, 43(6): 1151-1161. (in Chinese)

    [20] 鄭艷梅, 黃光耀, 蘆碧波, 王永茂. 結合深度學習和引導濾波的蘋果葉片圖像分割. 中國農(nóng)業(yè)大學學報, 2020, 25(9): 112-120.

    ZHENG Y M, HUANG G Y, LU B B, WANG Y M. An apple leaf image segmentation algorithm using deep learning and guided filtering. Journal of China Agricultural University, 2020, 25(9): 112-120. (in Chinese)

    [21] 曲之琳. 深度學習網(wǎng)絡在蘋果表面病變圖像識別中的應用[D]. 南京: 南京郵電大學, 2017.

    QU Z L. Application of deep learning network in recognition of apple surface lesions [D]. Nanjing: Nanjing University of Posts and Telecommunications, 2017. (in Chinese)

    [22] MTIBAA F, NGUYEN K K, AZAM M, PAPACHRISTOU A, VENNE J S, CHERIET M. LSTM-based indoor air temperature prediction framework for HVAC systems in smart buildings. Neural Computing and Applications, 2020, 32(23): 17569-17585.

    [23] BLANQUERO R, CARRIZOSA E, MOLERO-RIO C, MORALES D R. Sparsity in optimal randomized classification trees. 2020, 284(1): 255-272.

    [24] Breiman L. Random Forests. Machine Learning, 2001, 45(1): 5-32.

    [25] GISLASON P O, BENEDIKTSSON J A, SVEINSSON J R. Random Forests for land cover classification. Pattern Recognition Letters, 2006, 27(4): 294-300.

    [26] DE ANDRADE C F, DELGADO R C, BARBOSA M L, TEODORO P E, DA SILVA J C A, WANDERLEY H S, CAPRISTO-SILVA G F. Fire regime in Southern Brazil driven by atmospheric variation and vegetation cover. Agricultural and Forest Meteorology, 2020, 295: 108194.

    [27] MUKAKA M M. Statics corner: A guide to appropriate use of correlation coefficient in medical research. Malawi Medical Journal, 2012, 24(3): 69-71.

    [28] VAJSOVá B, FASBENDER D, WIRNHARDT C, LEMAJIC S, DEVOS W. Assessing spatial limits of Sentinel-2 data on arable crops in the context of checks by monitoring. Remote Sensing, 2020, 12(14): 2195.

    [29] LIU C W, KANG S Z, LI F S, LI S E, DU T S, TONG L. Relationship between environmental factor and maximum daily stem shrinkage in apple tree in arid region of northwest China. Scientia Horticulturae, 2011, 130(1): 118-125.

    [30] ORTUNO M F, CONEJERO W, MORENO F, MORIANA A, INTRIGLIOLO D S, BIEL C, MELLISHO C D, PEREZ-PASTOR A, DOMINGO R, RUIZ-SANCHEZ M C, CASADESUS J, BONANY J, TORRECILLAS A. Could trunk diameter sensors be used in woody crops for irrigation scheduling? A review of current knowledge and future perspectives. Agricultural Water Management, 2010, 97(1): 1-11.

    [31] EGEA G, PAGáN E, BAILLE A, DOMINGO R, NORTES P A, PéREZ-PASTOR A. Usefulness of establishing trunk diameter based reference lines for irrigation scheduling in almond trees. Irrigation Science, 2009, 27(6): 431-441.

    [32] BERMAN M E, DEJONG T M. Seasonal patterns of vegetative growth and competition with reproductive sinks in peach (). The Journal of Horticultural Science and Biotechnology, 2003, 78(3): 303-309.

    [33] 孫紅, 李松, 李民贊, 劉豪杰, 喬浪, 張瑤. 農(nóng)業(yè)信息成像感知與深度學習應用研究進展. 農(nóng)業(yè)機械學報, 2020, 51(5): 1-17.

    SUN H, LI S, LI M Z, LIU H J, QIAO L, ZHANG Y. Research progress of image sensing and deep learning in agriculture. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(5): 1-17. (in Chinese)

    [34] 邵澤中, 姚青, 唐健, 李罕瓊, 楊保軍, 呂軍, 陳軼. 面向移動終端的農(nóng)業(yè)害蟲圖像智能識別系統(tǒng)的研究與開發(fā). 中國農(nóng)業(yè)科學, 2020, 53(16): 3257-3268.

    SHAO Z Z, YAO Q, TANG J, LI H Q, YANG B J, LV J, CHEN Y. Research and development of the intelligent identification system of agricultural pests for mobile terminals. Scientia Agricultura Sinica, 2020, 53(16): 3257-3268. (in Chinese)

    [35] 劉慧芳, 賀正, 賈彪, 劉志, 李振洲, 付江鵬, 慕瑞瑞, 康建宏. 基于機器學習的滴灌玉米光合響應特征. 中國農(nóng)業(yè)科學, 2019, 52(17): 2939-2950.

    LIU H F, HE Z, JIA B, LIU Z, LI Z Z, FU J P, MU R R, KANG J H. Photosynthetic response characteristics of maize under drip irrigation based on machine learning. Scientia Agricultura Sinica, 2019, 52(17): 2939-2950. (in Chinese)

    Fruit Growth Modelling Based on Multi-Methods - A Case Study of Apple in Zhaotong, Yunnan

    1Chinese Academy of Meteorological Sciences, Beijing 100081;2Yunnan Climate Center, Kunming 650000;3Agricultural Meteorological Experimental Station of Zhaotong, Zhaotong 657000, Yunnan

    【】Meteorological factors are closely related to fruit diameter during growth process, but this relationship between them tends to be non-linear and non-stationary, thus making it hard to monitor the fruit and trunk diameter continuously. Comparing the simulation capabilities of various growth models for fruit diameter could provide scientific support for fruit growth monitoring and predicting, timely irrigation and fertilization, and the regulation of growth environment. 【】Taking apples in Zhaotong, Yunan Province as an example, this study first analyzed the characteristics of diameter change during apple growth in 2019 and 2020 and its relationship with environmental and climate factors. Subsequently, a deep learning method of Long Short-Term Memory (LSTM) model was adopted to simulate and predict the fruit diameter by integrating these factors, which was evaluated with the multi-linear regression (MLR) model and machine learning methods including Decision Tree (DT) and Random Forests (RF) using three sampling methods. 【】The apple diameter had obvious diurnal cycle characteristics, which shrunk in the daytime and expanded in the nighttime. The maximum diameter was in the morning, while the minimum diameter was near the sunset. The growth rate of apple diameter was higher in the early growth period than near mature. The hourly and daily mean apple diameters were moderately or highly-positive correlated with soil temperature and soil moisture, while there was a highly-negative correlation with UVI. The daily mean increase (FMDG), daily increase (FDG), and maximum daily shrinkage (MDFS) of apple diameter had a weak negative correlation with 60 cm soil temperature as well as 20 and 40 cm soil moisture (-0.5≤<-0.3). The simulation accuracy of the LSTM model was significantly higher than that of MLR, DT and RF model. The correlation coefficient () of LSTM model increased (3% -20%) compared with MLR, and the RMSE and MAE were approximately decreased by 50%-75%.The machine learning methods showed relatively poor performance in apple diameter simulation and might have overfitting problems.【】 Compared to statistics and machine learning approaches, the LSTM model demonstrated higher accuracy and robust performance because of the incapability of considering the complex non-linear correlations in the fruit growth simulation.

    apple diameter; growth model; meteorological factors; deep learning; Long Short-Term Memory (LSTM)

    10.3864/j.issn.0578-1752.2021.17.015

    2020-11-25;

    2020-12-21

    國家重點研發(fā)計劃(2019YFD1002201)、云南省科技計劃(2018BC007)、云南省省部合作協(xié)議重點工程“高原特色農(nóng)業(yè)氣象服務系統(tǒng)建設專項”資助項目(2019.01—2021.12)、中國氣象科學研究院基本科研業(yè)務費(2020Y003)

    孫擎,E-mail:sunq@cma.gov.cn;sunqingmeteo@gmail.com。通信作者趙艷霞,E-mail:zhaoyanxia@cma.gov.cn

    (責任編輯 趙伶俐)

    猜你喜歡
    莖干直徑蘋果
    基于L1中值骨架提取的植物莖干補全研究
    各顯神通測直徑
    基于潛熱效應的活立木凍融檢測傳感器設計與實驗
    山水(直徑40cm)
    云南檔案(2019年7期)2019-08-06 03:40:50
    基于干型差異的云南松子代苗木生物量研究
    種子(2019年5期)2019-07-02 11:37:44
    收獲蘋果
    有毒的蘋果
    拿蘋果
    巖木瓜莖干化學成分及PTP1B抑制活性研究
    會說話的蘋果
    最新美女视频免费是黄的| 99久久综合精品五月天人人| 熟女少妇亚洲综合色aaa.| 久久久久久久久免费视频了| 女同久久另类99精品国产91| 国产不卡一卡二| 久久久国产精品麻豆| 后天国语完整版免费观看| 亚洲精品国产色婷婷电影| 一本大道久久a久久精品| 亚洲 国产 在线| 日本精品一区二区三区蜜桃| 99riav亚洲国产免费| 亚洲avbb在线观看| 变态另类成人亚洲欧美熟女 | www.熟女人妻精品国产| 欧美成人午夜精品| 亚洲激情在线av| 精品免费久久久久久久清纯| 夜夜夜夜夜久久久久| 国产精品野战在线观看| 嫩草影视91久久| 老汉色av国产亚洲站长工具| 国产麻豆成人av免费视频| 免费观看精品视频网站| 亚洲 国产 在线| 国产精品九九99| 两性夫妻黄色片| 中文字幕最新亚洲高清| 久久精品91蜜桃| 在线观看一区二区三区| 啪啪无遮挡十八禁网站| 亚洲精品国产色婷婷电影| 变态另类成人亚洲欧美熟女 | 亚洲精品av麻豆狂野| 妹子高潮喷水视频| 精品乱码久久久久久99久播| 午夜精品久久久久久毛片777| 亚洲在线自拍视频| 午夜福利高清视频| 久久久精品欧美日韩精品| 久久香蕉激情| 亚洲av成人一区二区三| 久久人人精品亚洲av| 日韩视频一区二区在线观看| 亚洲成国产人片在线观看| 91在线观看av| 啦啦啦免费观看视频1| 欧美成人性av电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 美女国产高潮福利片在线看| 亚洲激情在线av| 国产成人影院久久av| 伦理电影免费视频| 在线观看www视频免费| 亚洲,欧美精品.| 这个男人来自地球电影免费观看| www.精华液| 国产精品国产高清国产av| 一级a爱片免费观看的视频| 免费观看人在逋| 最新在线观看一区二区三区| 一边摸一边抽搐一进一出视频| 久久国产精品男人的天堂亚洲| 免费女性裸体啪啪无遮挡网站| 首页视频小说图片口味搜索| 激情在线观看视频在线高清| av欧美777| 久久久久久久精品吃奶| 级片在线观看| 老熟妇仑乱视频hdxx| 琪琪午夜伦伦电影理论片6080| 操出白浆在线播放| 无限看片的www在线观看| 村上凉子中文字幕在线| 久久久久久久精品吃奶| 国产片内射在线| 国产精品亚洲一级av第二区| 19禁男女啪啪无遮挡网站| 一级作爱视频免费观看| 真人做人爱边吃奶动态| 波多野结衣高清无吗| 国产成人精品在线电影| 最新在线观看一区二区三区| 久热爱精品视频在线9| 丰满的人妻完整版| 久9热在线精品视频| 亚洲伊人色综图| 亚洲精品美女久久久久99蜜臀| 后天国语完整版免费观看| 精品国内亚洲2022精品成人| 91成人精品电影| 欧美日韩乱码在线| 中出人妻视频一区二区| 侵犯人妻中文字幕一二三四区| 国产蜜桃级精品一区二区三区| 免费看十八禁软件| 怎么达到女性高潮| 男人舔女人的私密视频| 久久人妻熟女aⅴ| 欧美成人免费av一区二区三区| 岛国在线观看网站| 免费高清视频大片| 亚洲一区高清亚洲精品| 亚洲国产精品合色在线| av视频免费观看在线观看| 丝袜人妻中文字幕| 亚洲欧洲精品一区二区精品久久久| 99国产精品一区二区三区| 色哟哟哟哟哟哟| 精品国产美女av久久久久小说| 天天添夜夜摸| 亚洲精华国产精华精| 国产97色在线日韩免费| 久久久久精品国产欧美久久久| 免费在线观看视频国产中文字幕亚洲| 日韩精品青青久久久久久| 欧美黑人精品巨大| 久久久久久国产a免费观看| 18美女黄网站色大片免费观看| 亚洲免费av在线视频| 91大片在线观看| 国产精品亚洲美女久久久| 不卡一级毛片| 日韩 欧美 亚洲 中文字幕| 19禁男女啪啪无遮挡网站| 18美女黄网站色大片免费观看| av片东京热男人的天堂| 人妻丰满熟妇av一区二区三区| 岛国在线观看网站| 久久欧美精品欧美久久欧美| 18禁国产床啪视频网站| 欧美成人性av电影在线观看| 麻豆国产av国片精品| 多毛熟女@视频| 午夜激情av网站| 久久中文看片网| 国产熟女午夜一区二区三区| 成人永久免费在线观看视频| 午夜福利,免费看| 国产欧美日韩综合在线一区二区| 国产精品一区二区精品视频观看| www.999成人在线观看| 操出白浆在线播放| 热re99久久国产66热| 人人妻人人澡欧美一区二区 | 又紧又爽又黄一区二区| 国产精品二区激情视频| 国产精品秋霞免费鲁丝片| 中文字幕精品免费在线观看视频| 国产高清激情床上av| 男人的好看免费观看在线视频 | 久久国产精品人妻蜜桃| 午夜影院日韩av| 人成视频在线观看免费观看| svipshipincom国产片| 亚洲国产欧美日韩在线播放| 日日干狠狠操夜夜爽| 国产精品一区二区在线不卡| 自线自在国产av| 精品免费久久久久久久清纯| 可以在线观看毛片的网站| 欧美日本亚洲视频在线播放| 啪啪无遮挡十八禁网站| 91精品国产国语对白视频| 亚洲美女黄片视频| 日韩大码丰满熟妇| 热re99久久国产66热| 亚洲国产看品久久| 日本一区二区免费在线视频| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 一卡2卡三卡四卡精品乱码亚洲| 午夜两性在线视频| 欧美日韩乱码在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲av片天天在线观看| 免费不卡黄色视频| 伦理电影免费视频| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区三区在线| 91麻豆精品激情在线观看国产| 麻豆av在线久日| 正在播放国产对白刺激| 色老头精品视频在线观看| 男女下面进入的视频免费午夜 | 欧美性长视频在线观看| 大型黄色视频在线免费观看| 亚洲中文av在线| 久99久视频精品免费| 久久国产精品人妻蜜桃| 又黄又爽又免费观看的视频| 99国产极品粉嫩在线观看| 亚洲久久久国产精品| 18美女黄网站色大片免费观看| 成人亚洲精品av一区二区| 超碰成人久久| 亚洲精品在线观看二区| 国产一区二区激情短视频| 国产av又大| 国产人伦9x9x在线观看| 国产精品免费视频内射| 久久狼人影院| 51午夜福利影视在线观看| x7x7x7水蜜桃| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 最近最新中文字幕大全免费视频| 亚洲欧洲精品一区二区精品久久久| 精品国产亚洲在线| 亚洲精品美女久久av网站| 久久久国产成人免费| 久久伊人香网站| 亚洲欧美激情综合另类| 欧美日本视频| 不卡av一区二区三区| 国产又色又爽无遮挡免费看| 国产精品久久久久久人妻精品电影| 一进一出抽搐动态| 一区二区三区国产精品乱码| 国产成人欧美| 中文字幕人妻丝袜一区二区| 久久天躁狠狠躁夜夜2o2o| 啦啦啦观看免费观看视频高清 | 人成视频在线观看免费观看| 大码成人一级视频| 色婷婷久久久亚洲欧美| 天天躁狠狠躁夜夜躁狠狠躁| 精品国产乱码久久久久久男人| 99久久精品国产亚洲精品| 久久精品国产综合久久久| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品在线电影| 亚洲一区高清亚洲精品| 黄网站色视频无遮挡免费观看| 国内毛片毛片毛片毛片毛片| 一级毛片精品| 国产精品,欧美在线| 亚洲成a人片在线一区二区| 国产精品亚洲一级av第二区| 免费女性裸体啪啪无遮挡网站| 精品一区二区三区av网在线观看| 午夜福利欧美成人| 久久精品91无色码中文字幕| 国产成人欧美在线观看| 色播亚洲综合网| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 人人妻,人人澡人人爽秒播| 中文字幕另类日韩欧美亚洲嫩草| 又紧又爽又黄一区二区| 黄片小视频在线播放| av福利片在线| 精品高清国产在线一区| 不卡一级毛片| 淫秽高清视频在线观看| 精品无人区乱码1区二区| 午夜福利影视在线免费观看| 波多野结衣av一区二区av| 国产高清有码在线观看视频 | 激情视频va一区二区三区| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 日韩欧美一区视频在线观看| www国产在线视频色| 18禁观看日本| a在线观看视频网站| 成人永久免费在线观看视频| 亚洲性夜色夜夜综合| 两个人看的免费小视频| 1024香蕉在线观看| 亚洲男人的天堂狠狠| 久久久久亚洲av毛片大全| 午夜福利欧美成人| 大码成人一级视频| 国产成人系列免费观看| 性欧美人与动物交配| ponron亚洲| 99riav亚洲国产免费| 亚洲国产精品合色在线| 亚洲性夜色夜夜综合| 国产精品电影一区二区三区| 天堂动漫精品| 91老司机精品| 国产成年人精品一区二区| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品免费视频一区二区三区| 成年女人毛片免费观看观看9| 91字幕亚洲| 国产色视频综合| 亚洲最大成人中文| 久久久精品欧美日韩精品| 中文字幕久久专区| 免费在线观看完整版高清| 亚洲国产中文字幕在线视频| 人人妻人人澡欧美一区二区 | 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 国产三级在线视频| 精品无人区乱码1区二区| 悠悠久久av| 18禁裸乳无遮挡免费网站照片 | 亚洲精品中文字幕在线视频| 青草久久国产| 亚洲人成网站在线播放欧美日韩| 最近最新免费中文字幕在线| 国语自产精品视频在线第100页| 香蕉丝袜av| 精品久久久久久久人妻蜜臀av | 丰满人妻熟妇乱又伦精品不卡| 一级a爱片免费观看的视频| 日本vs欧美在线观看视频| 午夜视频精品福利| 可以在线观看的亚洲视频| 在线观看www视频免费| 中国美女看黄片| 老熟妇乱子伦视频在线观看| 最近最新中文字幕大全免费视频| 欧美乱色亚洲激情| 精品久久久久久久毛片微露脸| 精品国产乱子伦一区二区三区| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 久久精品国产99精品国产亚洲性色 | 国产成人免费无遮挡视频| 亚洲成人国产一区在线观看| 婷婷丁香在线五月| 一边摸一边抽搐一进一出视频| 欧美中文综合在线视频| 在线观看66精品国产| 精品久久久久久久毛片微露脸| 黄色a级毛片大全视频| 国产成人系列免费观看| 亚洲五月婷婷丁香| 男男h啪啪无遮挡| 国产又色又爽无遮挡免费看| 精品无人区乱码1区二区| 麻豆成人av在线观看| 视频区欧美日本亚洲| 欧美一级毛片孕妇| 国产日韩一区二区三区精品不卡| 成年版毛片免费区| 久热这里只有精品99| 日韩有码中文字幕| 免费观看精品视频网站| 老司机靠b影院| 日韩视频一区二区在线观看| av福利片在线| 亚洲自偷自拍图片 自拍| 中文字幕人妻丝袜一区二区| www国产在线视频色| 国产一区二区三区视频了| 91av网站免费观看| 亚洲在线自拍视频| 国产又爽黄色视频| 久久午夜亚洲精品久久| 亚洲专区字幕在线| 美女高潮喷水抽搐中文字幕| 国产单亲对白刺激| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 日韩成人在线观看一区二区三区| 精品欧美一区二区三区在线| 男女下面进入的视频免费午夜 | 成人国语在线视频| 亚洲国产精品sss在线观看| 女人被狂操c到高潮| 波多野结衣av一区二区av| 成人国语在线视频| 午夜福利欧美成人| 高清在线国产一区| 深夜精品福利| 久久 成人 亚洲| 成人免费观看视频高清| 国产熟女xx| 日韩欧美免费精品| 岛国在线观看网站| 1024视频免费在线观看| 精品国产超薄肉色丝袜足j| 在线观看免费日韩欧美大片| 国产99白浆流出| 真人一进一出gif抽搐免费| 妹子高潮喷水视频| 久久香蕉激情| av天堂在线播放| 欧美色欧美亚洲另类二区 | www.精华液| 男女做爰动态图高潮gif福利片 | 国产又爽黄色视频| 丰满人妻熟妇乱又伦精品不卡| 免费av毛片视频| 亚洲精品国产色婷婷电影| 1024香蕉在线观看| 亚洲 欧美一区二区三区| 国产三级黄色录像| 一本综合久久免费| 国产精品久久久久久人妻精品电影| 一本大道久久a久久精品| 99在线人妻在线中文字幕| 日韩一卡2卡3卡4卡2021年| 又大又爽又粗| 亚洲中文日韩欧美视频| 他把我摸到了高潮在线观看| 欧美色欧美亚洲另类二区 | videosex国产| 亚洲情色 制服丝袜| 欧美日本亚洲视频在线播放| 好男人在线观看高清免费视频 | 欧美日韩黄片免| 69av精品久久久久久| 日日爽夜夜爽网站| 免费av毛片视频| 老司机福利观看| 日韩欧美国产一区二区入口| 99在线视频只有这里精品首页| 亚洲欧美日韩另类电影网站| 久久中文字幕一级| 啦啦啦韩国在线观看视频| 欧美老熟妇乱子伦牲交| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 一级片免费观看大全| 757午夜福利合集在线观看| 一级片免费观看大全| 国产亚洲精品久久久久久毛片| 亚洲色图 男人天堂 中文字幕| 丝袜在线中文字幕| 夜夜夜夜夜久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 免费久久久久久久精品成人欧美视频| 久久精品国产清高在天天线| 中文字幕久久专区| 丝袜美腿诱惑在线| 手机成人av网站| 亚洲精品久久国产高清桃花| 久久婷婷成人综合色麻豆| 久久久精品国产亚洲av高清涩受| 啦啦啦免费观看视频1| 色尼玛亚洲综合影院| 亚洲男人天堂网一区| 亚洲,欧美精品.| 久久性视频一级片| www.自偷自拍.com| 在线观看免费午夜福利视频| 成人特级黄色片久久久久久久| 电影成人av| 激情在线观看视频在线高清| 18禁美女被吸乳视频| 欧美国产精品va在线观看不卡| 露出奶头的视频| 国内精品久久久久精免费| 一进一出好大好爽视频| 亚洲国产精品久久男人天堂| 久久亚洲精品不卡| 非洲黑人性xxxx精品又粗又长| 91字幕亚洲| 在线观看免费午夜福利视频| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 又紧又爽又黄一区二区| 最近最新免费中文字幕在线| 欧美精品亚洲一区二区| 精品免费久久久久久久清纯| 岛国在线观看网站| 亚洲国产精品成人综合色| 精品福利观看| 国产亚洲精品久久久久5区| 国产精品电影一区二区三区| 国产一区二区激情短视频| 国产1区2区3区精品| 19禁男女啪啪无遮挡网站| 中文字幕av电影在线播放| 啦啦啦 在线观看视频| 精品国产乱子伦一区二区三区| 99热只有精品国产| 国产高清videossex| 久久精品国产清高在天天线| 天堂影院成人在线观看| 亚洲成人精品中文字幕电影| 国内精品久久久久久久电影| ponron亚洲| 中文字幕人妻熟女乱码| 动漫黄色视频在线观看| 国产色视频综合| 国语自产精品视频在线第100页| 中文字幕另类日韩欧美亚洲嫩草| av电影中文网址| 国产一区二区三区在线臀色熟女| 丁香六月欧美| 亚洲精品美女久久久久99蜜臀| 脱女人内裤的视频| 黄色视频不卡| 亚洲色图综合在线观看| 黄色视频不卡| 亚洲中文日韩欧美视频| 国产av精品麻豆| 日韩国内少妇激情av| 免费人成视频x8x8入口观看| 亚洲午夜理论影院| 露出奶头的视频| 国产精品野战在线观看| 日本 av在线| aaaaa片日本免费| 婷婷精品国产亚洲av在线| 99久久精品国产亚洲精品| 纯流量卡能插随身wifi吗| 国产精华一区二区三区| 日本黄色视频三级网站网址| 桃红色精品国产亚洲av| 亚洲人成网站在线播放欧美日韩| 丝袜人妻中文字幕| 中文亚洲av片在线观看爽| 亚洲人成电影观看| 国产精品亚洲一级av第二区| a在线观看视频网站| 亚洲电影在线观看av| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女 | 国产免费男女视频| 欧美激情高清一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 大型av网站在线播放| 亚洲午夜精品一区,二区,三区| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| 一级,二级,三级黄色视频| 免费无遮挡裸体视频| 国产色视频综合| 天天躁狠狠躁夜夜躁狠狠躁| 50天的宝宝边吃奶边哭怎么回事| 精品免费久久久久久久清纯| 久久婷婷成人综合色麻豆| 欧美一级毛片孕妇| 久久久水蜜桃国产精品网| 两性午夜刺激爽爽歪歪视频在线观看 | 色在线成人网| 精品一区二区三区av网在线观看| 一夜夜www| 激情视频va一区二区三区| av片东京热男人的天堂| 国产精品免费一区二区三区在线| 欧美日本视频| 国产亚洲av嫩草精品影院| 欧美av亚洲av综合av国产av| 怎么达到女性高潮| 性色av乱码一区二区三区2| 身体一侧抽搐| 老司机福利观看| 国产精品九九99| av免费在线观看网站| 很黄的视频免费| 欧美老熟妇乱子伦牲交| 亚洲av成人一区二区三| 在线观看午夜福利视频| 精品人妻1区二区| xxx96com| 亚洲七黄色美女视频| 老熟妇仑乱视频hdxx| 美女扒开内裤让男人捅视频| 女性生殖器流出的白浆| 亚洲午夜精品一区,二区,三区| 国产国语露脸激情在线看| 两个人免费观看高清视频| 久久精品国产综合久久久| 禁无遮挡网站| 中文字幕人成人乱码亚洲影| 一级毛片精品| 欧美成人一区二区免费高清观看 | 中文字幕最新亚洲高清| 少妇裸体淫交视频免费看高清 | 一区二区三区精品91| 免费在线观看影片大全网站| 国产成人av教育| 女人高潮潮喷娇喘18禁视频| 夜夜夜夜夜久久久久| 日韩精品免费视频一区二区三区| 熟女少妇亚洲综合色aaa.| 黄色女人牲交| 夜夜躁狠狠躁天天躁| 久久精品人人爽人人爽视色| 一本大道久久a久久精品| 啪啪无遮挡十八禁网站| 亚洲电影在线观看av| 一进一出抽搐动态| 亚洲国产毛片av蜜桃av| 级片在线观看| 国产私拍福利视频在线观看| 国产精品野战在线观看| 久久人妻熟女aⅴ| 国产精品亚洲一级av第二区| 午夜精品久久久久久毛片777| 男女午夜视频在线观看| 午夜成年电影在线免费观看| 91精品三级在线观看| 91九色精品人成在线观看| 一个人观看的视频www高清免费观看 | 亚洲在线自拍视频| 成在线人永久免费视频| 精品福利观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲美女黄片视频| 日韩欧美国产在线观看| 美女高潮到喷水免费观看| 欧美绝顶高潮抽搐喷水| 久久久久国内视频| 欧美日本视频| 两个人免费观看高清视频| 亚洲欧美日韩高清在线视频| 亚洲人成77777在线视频| 麻豆一二三区av精品| 两个人视频免费观看高清| 亚洲熟女毛片儿| 一级作爱视频免费观看| 男女午夜视频在线观看|