• <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在线人妻在线中文字幕| 女生性感内裤真人,穿戴方法视频| 老熟妇仑乱视频hdxx| 99久久99久久久精品蜜桃| 欧美黄色淫秽网站| 麻豆av在线久日| av在线天堂中文字幕| av欧美777| 最近在线观看免费完整版| 亚洲中文av在线| 男女下面进入的视频免费午夜| 国产欧美日韩一区二区三| cao死你这个sao货| 特大巨黑吊av在线直播| 老司机深夜福利视频在线观看| 好男人电影高清在线观看| 亚洲熟妇中文字幕五十中出| bbb黄色大片| 国产探花在线观看一区二区| 日韩欧美 国产精品| 免费在线观看完整版高清| 亚洲五月天丁香| 熟女少妇亚洲综合色aaa.| 国产精品亚洲美女久久久| 亚洲真实伦在线观看| a级毛片a级免费在线| av福利片在线| 在线国产一区二区在线| 亚洲成人精品中文字幕电影| 99国产精品一区二区蜜桃av| 精品久久久久久久久久久久久| 欧美激情久久久久久爽电影| www日本黄色视频网| 制服丝袜大香蕉在线| 97超级碰碰碰精品色视频在线观看| 在线国产一区二区在线| 美女大奶头视频| 香蕉丝袜av| 成年免费大片在线观看| 日本熟妇午夜| 国产精品久久久久久精品电影| 99久久无色码亚洲精品果冻| 国产又黄又爽又无遮挡在线| 两个人的视频大全免费| 亚洲狠狠婷婷综合久久图片| 精品欧美一区二区三区在线| 一二三四在线观看免费中文在| 久久久国产欧美日韩av| 一进一出好大好爽视频| 久久婷婷人人爽人人干人人爱| 国产精品,欧美在线| 俄罗斯特黄特色一大片| 熟女少妇亚洲综合色aaa.| 欧美不卡视频在线免费观看 | 18美女黄网站色大片免费观看| 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 午夜福利18| 叶爱在线成人免费视频播放| 欧美极品一区二区三区四区| av免费在线观看网站| 午夜a级毛片| 亚洲aⅴ乱码一区二区在线播放 | 欧美性猛交黑人性爽| 国产成人aa在线观看| 国产精品电影一区二区三区| 国产三级中文精品| 嫩草影院精品99| 夜夜夜夜夜久久久久| 国产1区2区3区精品| 日韩免费av在线播放| 欧美色欧美亚洲另类二区| 国产在线观看jvid| 欧美 亚洲 国产 日韩一| 久久久久国产精品人妻aⅴ院| 国内少妇人妻偷人精品xxx网站 | 日韩欧美免费精品| www.精华液| 少妇熟女aⅴ在线视频| 免费看a级黄色片| 我要搜黄色片| 级片在线观看| 大型黄色视频在线免费观看| 欧美性猛交╳xxx乱大交人| 国产视频内射| 亚洲欧美精品综合一区二区三区| 国产精品一及| 国产亚洲精品av在线| 亚洲全国av大片| 成在线人永久免费视频| 成人亚洲精品av一区二区| 99国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 午夜激情福利司机影院| 三级毛片av免费| 亚洲av片天天在线观看| 欧美黑人精品巨大| av福利片在线| 国内精品一区二区在线观看| 免费av毛片视频| 国产亚洲av嫩草精品影院| 观看免费一级毛片| av天堂在线播放| 一进一出好大好爽视频| 亚洲av第一区精品v没综合| 在线国产一区二区在线| 男人舔女人的私密视频| 少妇被粗大的猛进出69影院| 亚洲精品国产精品久久久不卡| 国产一区二区激情短视频| 中国美女看黄片| 中文字幕最新亚洲高清| 国产精品亚洲一级av第二区| 两个人免费观看高清视频| 欧美性猛交╳xxx乱大交人| 俺也久久电影网| 黄色丝袜av网址大全| 99久久99久久久精品蜜桃| 久久香蕉精品热| av国产免费在线观看| 亚洲av第一区精品v没综合| 亚洲人成电影免费在线| 蜜桃久久精品国产亚洲av| 免费观看精品视频网站| 在线观看66精品国产| 久久久久亚洲av毛片大全| 美女高潮喷水抽搐中文字幕| 欧美成人一区二区免费高清观看 | 国产成人精品久久二区二区91| 成年女人毛片免费观看观看9| 天天一区二区日本电影三级| 国产精品影院久久| 国产精品一区二区三区四区久久| 午夜精品在线福利| 怎么达到女性高潮| 欧美乱码精品一区二区三区| 小说图片视频综合网站| 色在线成人网| 亚洲av第一区精品v没综合| 欧美一级毛片孕妇| 国产成人av激情在线播放| 亚洲 欧美一区二区三区| 日本黄色视频三级网站网址| 久久久久久久久久黄片| 最近最新中文字幕大全免费视频| АⅤ资源中文在线天堂| 国产精品 国内视频| 狂野欧美白嫩少妇大欣赏| 日韩大码丰满熟妇| 欧美人与性动交α欧美精品济南到| 在线播放国产精品三级| 99热这里只有是精品50| 特级一级黄色大片| 日本免费一区二区三区高清不卡| 婷婷六月久久综合丁香| a级毛片在线看网站| 两人在一起打扑克的视频| 国产成人啪精品午夜网站| 岛国在线观看网站| 少妇被粗大的猛进出69影院| 久久精品91无色码中文字幕| 国产成人精品无人区| 嫁个100分男人电影在线观看| 久99久视频精品免费| 欧美一级毛片孕妇| 伦理电影免费视频| 亚洲av中文字字幕乱码综合| 亚洲自偷自拍图片 自拍| 欧美成人午夜精品| 天天躁夜夜躁狠狠躁躁| 欧美在线黄色| 亚洲熟妇熟女久久| 制服丝袜大香蕉在线| 嫁个100分男人电影在线观看| 脱女人内裤的视频| 在线观看一区二区三区| 99精品久久久久人妻精品| 国产午夜精品论理片| av欧美777| 国产熟女午夜一区二区三区| 亚洲午夜精品一区,二区,三区| 丰满人妻一区二区三区视频av | 美女午夜性视频免费| 中文字幕熟女人妻在线| 丝袜人妻中文字幕| 国产av一区在线观看免费| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美精品v在线| 久久久久国内视频| 日韩国内少妇激情av| 国产高清激情床上av| 久久久国产成人免费| 无限看片的www在线观看| 久久天堂一区二区三区四区| 国产亚洲欧美98| 国产黄片美女视频| 久久久精品欧美日韩精品| 色尼玛亚洲综合影院| 欧美一区二区精品小视频在线| 亚洲av第一区精品v没综合| 亚洲av成人精品一区久久| 天天躁夜夜躁狠狠躁躁| 男人的好看免费观看在线视频 | 免费人成视频x8x8入口观看| 99国产精品99久久久久| 黄片小视频在线播放| 免费在线观看亚洲国产| 身体一侧抽搐| 老鸭窝网址在线观看| 窝窝影院91人妻| 国产又色又爽无遮挡免费看| 精品一区二区三区四区五区乱码| 亚洲中文字幕一区二区三区有码在线看 | 免费av毛片视频| 精品一区二区三区四区五区乱码| 身体一侧抽搐| 这个男人来自地球电影免费观看| 欧美绝顶高潮抽搐喷水| 十八禁网站免费在线| 久久精品国产99精品国产亚洲性色| 一边摸一边做爽爽视频免费| 精品免费久久久久久久清纯| 麻豆一二三区av精品| 久久精品aⅴ一区二区三区四区| 在线观看免费午夜福利视频| 亚洲电影在线观看av| 一个人免费在线观看电影 | 999精品在线视频| 国产黄色小视频在线观看| 黑人欧美特级aaaaaa片| 男人舔女人下体高潮全视频| 国产麻豆成人av免费视频| 欧美精品亚洲一区二区| 免费看美女性在线毛片视频| 免费观看精品视频网站| 国产91精品成人一区二区三区| 国内精品久久久久久久电影| 三级国产精品欧美在线观看 | 久9热在线精品视频| 欧美日本视频| 国产精品一区二区三区四区久久| 欧美极品一区二区三区四区| 欧美+亚洲+日韩+国产| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清在线视频| 欧美在线黄色| 黄色片一级片一级黄色片| 曰老女人黄片| 最近视频中文字幕2019在线8| 国产精品久久电影中文字幕| 国产精品久久久久久亚洲av鲁大| 黄片大片在线免费观看| 欧美性长视频在线观看| 久久久久久国产a免费观看| 欧美国产日韩亚洲一区| 色在线成人网| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| 在线永久观看黄色视频| av片东京热男人的天堂| 亚洲成人久久性| 亚洲七黄色美女视频| 男女那种视频在线观看| 久久婷婷人人爽人人干人人爱| 很黄的视频免费| 国产精品免费视频内射| 国产亚洲精品久久久久5区| 一a级毛片在线观看| 成年人黄色毛片网站| 欧美绝顶高潮抽搐喷水| 在线免费观看的www视频| 欧美成狂野欧美在线观看| 九色国产91popny在线| 亚洲精品中文字幕在线视频| 一个人免费在线观看的高清视频| 少妇的丰满在线观看| 少妇熟女aⅴ在线视频| 色在线成人网| 一个人免费在线观看的高清视频| 香蕉国产在线看| 99热这里只有是精品50| 国产一区二区三区视频了| 变态另类丝袜制服| 后天国语完整版免费观看| 国产99久久九九免费精品| 精品一区二区三区视频在线观看免费| 美女黄网站色视频| 十八禁人妻一区二区| 亚洲精品av麻豆狂野| 香蕉av资源在线| 成年女人毛片免费观看观看9| 18禁观看日本| 99在线视频只有这里精品首页| 黄色成人免费大全| 日本免费a在线| 国产精品永久免费网站| 亚洲国产高清在线一区二区三| 美女午夜性视频免费| 国产av一区二区精品久久| 国产精品一及| 欧美另类亚洲清纯唯美| 午夜福利欧美成人| 国产亚洲欧美98| 亚洲 国产 在线| 后天国语完整版免费观看| 在线观看免费视频日本深夜| 免费看十八禁软件| 又紧又爽又黄一区二区| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 非洲黑人性xxxx精品又粗又长| 一区二区三区激情视频| 一级片免费观看大全| 欧美乱色亚洲激情| 日本免费a在线| 亚洲av成人一区二区三| 国产三级黄色录像| 一级毛片精品| 久久香蕉国产精品| 国产激情欧美一区二区| 国产一区二区激情短视频| 久久久精品大字幕| 99精品欧美一区二区三区四区| 校园春色视频在线观看| av国产免费在线观看| 国产区一区二久久| 久久精品人妻少妇| 老鸭窝网址在线观看| 成人精品一区二区免费| 久久香蕉国产精品| 成人18禁在线播放| 亚洲专区中文字幕在线| 亚洲人与动物交配视频| 黑人巨大精品欧美一区二区mp4| 最近视频中文字幕2019在线8| 精品久久久久久久久久久久久| 亚洲av成人不卡在线观看播放网| 国产精品 国内视频| 91老司机精品| 欧美日本亚洲视频在线播放| 一个人观看的视频www高清免费观看 | 精品久久久久久,| 99riav亚洲国产免费| 国产av麻豆久久久久久久| 国产av在哪里看| 丝袜人妻中文字幕| 成人一区二区视频在线观看| 日本五十路高清| 日本 欧美在线| 日日爽夜夜爽网站| 欧美在线黄色| 99精品久久久久人妻精品| 又黄又爽又免费观看的视频| АⅤ资源中文在线天堂| 亚洲精品国产一区二区精华液| 神马国产精品三级电影在线观看 | 搡老岳熟女国产| 桃色一区二区三区在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 最新在线观看一区二区三区| av中文乱码字幕在线| 久久久精品国产亚洲av高清涩受| 午夜激情av网站| 久久这里只有精品中国| 久久国产精品影院| 丰满的人妻完整版| aaaaa片日本免费| 午夜两性在线视频| 国产成人精品久久二区二区91| 麻豆一二三区av精品| 久久中文看片网| 美女扒开内裤让男人捅视频| 18禁观看日本| 一本一本综合久久| 免费无遮挡裸体视频| 香蕉丝袜av| 日日摸夜夜添夜夜添小说| 国产蜜桃级精品一区二区三区| 国产区一区二久久| 亚洲成人久久性| 久久精品成人免费网站| 日韩国内少妇激情av| 久久久久久国产a免费观看| 99热这里只有是精品50| 两个人免费观看高清视频| 狂野欧美激情性xxxx| 18禁观看日本| 香蕉国产在线看| 午夜两性在线视频| 岛国在线免费视频观看| 欧美日韩国产亚洲二区| 亚洲国产欧美人成| 日本三级黄在线观看| 欧美又色又爽又黄视频| 欧美日韩一级在线毛片| 欧美三级亚洲精品| cao死你这个sao货| 熟妇人妻久久中文字幕3abv| 亚洲中文av在线| 黑人欧美特级aaaaaa片| 国产真人三级小视频在线观看| 久久久久性生活片| 久久婷婷成人综合色麻豆| 欧美成人性av电影在线观看| 日韩大尺度精品在线看网址| 国产精品av视频在线免费观看| 欧美黄色淫秽网站| 伊人久久大香线蕉亚洲五| 亚洲欧美日韩高清专用| 麻豆成人av在线观看| 日日干狠狠操夜夜爽| 久久久久免费精品人妻一区二区| 亚洲精品美女久久久久99蜜臀| 一本精品99久久精品77| 久久婷婷人人爽人人干人人爱| 亚洲第一欧美日韩一区二区三区| 999久久久国产精品视频| 一本久久中文字幕| 在线a可以看的网站| 日本 欧美在线| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| 中文字幕av在线有码专区| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 成人三级黄色视频| 欧美色欧美亚洲另类二区| 日韩免费av在线播放| 国产在线观看jvid| 午夜免费激情av| 亚洲欧美精品综合一区二区三区| 精品久久久久久久毛片微露脸| 欧美国产日韩亚洲一区| 97人妻精品一区二区三区麻豆| 欧美色欧美亚洲另类二区| www.自偷自拍.com| 午夜精品久久久久久毛片777| 日韩中文字幕欧美一区二区| 一级片免费观看大全| 亚洲国产日韩欧美精品在线观看 | 国产精品久久电影中文字幕| 啦啦啦免费观看视频1| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 长腿黑丝高跟| 亚洲第一电影网av| 搡老岳熟女国产| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 最近最新免费中文字幕在线| 又黄又粗又硬又大视频| 久久精品综合一区二区三区| 亚洲av成人精品一区久久| 丰满的人妻完整版| 亚洲第一欧美日韩一区二区三区| 91九色精品人成在线观看| 人妻夜夜爽99麻豆av| 久久天躁狠狠躁夜夜2o2o| 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕人成人乱码亚洲影| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 亚洲中文av在线| 熟妇人妻久久中文字幕3abv| 久久久久国内视频| 一进一出抽搐gif免费好疼| 最新在线观看一区二区三区| 国产三级中文精品| 久99久视频精品免费| 久久精品91蜜桃| 18美女黄网站色大片免费观看| 黑人操中国人逼视频| 长腿黑丝高跟| 午夜久久久久精精品| 国产精品香港三级国产av潘金莲| 欧美黄色淫秽网站| 又粗又爽又猛毛片免费看| 丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲欧美在线一区二区| √禁漫天堂资源中文www| 无遮挡黄片免费观看| 一本一本综合久久| 天天一区二区日本电影三级| 变态另类丝袜制服| 欧美极品一区二区三区四区| 久久久久久人人人人人| 夜夜爽天天搞| 国产精品亚洲美女久久久| 久久天堂一区二区三区四区| 国产日本99.免费观看| 精品不卡国产一区二区三区| 99国产综合亚洲精品| 国产精品一区二区免费欧美| 丁香欧美五月| 1024手机看黄色片| 午夜激情av网站| 18禁美女被吸乳视频| 国产亚洲精品一区二区www| 精品乱码久久久久久99久播| 老司机午夜福利在线观看视频| 神马国产精品三级电影在线观看 | 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 国产高清videossex| 亚洲自拍偷在线| 成年女人毛片免费观看观看9| 日韩 欧美 亚洲 中文字幕| 色综合欧美亚洲国产小说| 久久精品夜夜夜夜夜久久蜜豆 | 成人欧美大片| 色综合亚洲欧美另类图片| 最新在线观看一区二区三区| 亚洲精品中文字幕在线视频| 亚洲av熟女| 国产亚洲精品一区二区www| 波多野结衣高清作品| 午夜成年电影在线免费观看| 日韩欧美精品v在线| 精品熟女少妇八av免费久了| 熟女电影av网| 在线观看午夜福利视频| 国产黄片美女视频| 男女视频在线观看网站免费 | 国产黄a三级三级三级人| 伦理电影免费视频| 国产成人系列免费观看| 操出白浆在线播放| 精品一区二区三区视频在线观看免费| a级毛片在线看网站| 亚洲专区字幕在线| 老汉色av国产亚洲站长工具| www.精华液| 男女做爰动态图高潮gif福利片| 成年女人毛片免费观看观看9| 国产亚洲精品第一综合不卡| 国产成人一区二区三区免费视频网站| 久久久久久九九精品二区国产 | 日韩免费av在线播放| 久久午夜亚洲精品久久| 成人av一区二区三区在线看| 久久九九热精品免费| 国产亚洲精品久久久久久毛片| 一边摸一边抽搐一进一小说| 国产亚洲av高清不卡| 桃红色精品国产亚洲av| 亚洲一区二区三区色噜噜| 嫁个100分男人电影在线观看| 亚洲电影在线观看av| 久久久久国内视频| 精品一区二区三区视频在线观看免费| 最近在线观看免费完整版| 97碰自拍视频| 国产成人一区二区三区免费视频网站| 亚洲精品在线观看二区| 国内精品久久久久久久电影| 嫩草影视91久久| 亚洲一卡2卡3卡4卡5卡精品中文| 999精品在线视频| 俺也久久电影网| 一级黄色大片毛片| 国产高清有码在线观看视频 | 不卡一级毛片| 国产av又大| 两个人的视频大全免费| 香蕉国产在线看| 国产午夜福利久久久久久| 麻豆av在线久日| 亚洲精品一卡2卡三卡4卡5卡| 啪啪无遮挡十八禁网站| 欧美日韩瑟瑟在线播放| 最近最新中文字幕大全免费视频| 国产午夜精品论理片| 听说在线观看完整版免费高清| 黑人操中国人逼视频| 91成年电影在线观看| 一级毛片高清免费大全| 九九热线精品视视频播放| 免费一级毛片在线播放高清视频| 国产一区二区在线av高清观看| 欧美最黄视频在线播放免费| 深夜精品福利| 国产在线精品亚洲第一网站| 久久久久久久久久黄片| 欧美久久黑人一区二区| 黄色视频,在线免费观看| 国产熟女xx| 久久中文字幕一级| 国产伦在线观看视频一区| 搡老妇女老女人老熟妇| 欧美精品啪啪一区二区三区| 欧美成人性av电影在线观看| 国产精品一区二区三区四区免费观看 | 中文字幕人成人乱码亚洲影| www.精华液| 亚洲激情在线av| 中文在线观看免费www的网站 | 欧美一级毛片孕妇| 午夜a级毛片| 中文字幕最新亚洲高清| 久久精品国产清高在天天线| 欧美久久黑人一区二区| 十八禁网站免费在线| 国产激情欧美一区二区| 露出奶头的视频| 亚洲国产日韩欧美精品在线观看 | netflix在线观看网站| 最近视频中文字幕2019在线8| 欧美乱色亚洲激情| 日日干狠狠操夜夜爽| 黑人巨大精品欧美一区二区mp4| 1024手机看黄色片| 国产三级在线视频| 久久草成人影院| 亚洲精品久久成人aⅴ小说|