• <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抑制活性研究
    會說話的蘋果
    精品久久久久久久久久久久久| 春色校园在线视频观看| 一个人观看的视频www高清免费观看| 久久人人爽人人爽人人片va| 国产亚洲午夜精品一区二区久久 | 人体艺术视频欧美日本| 三级国产精品欧美在线观看| 中文字幕av成人在线电影| 中文字幕av成人在线电影| 亚洲成人精品中文字幕电影| 99久久精品热视频| 九九热线精品视视频播放| 免费观看人在逋| 青春草亚洲视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美成人一区二区免费高清观看| 人妻少妇偷人精品九色| 亚洲一区高清亚洲精品| 亚洲国产精品久久男人天堂| 成人欧美大片| 一个人看视频在线观看www免费| 久久久久久久久久久免费av| 两个人的视频大全免费| 哪个播放器可以免费观看大片| 丰满人妻一区二区三区视频av| 色综合站精品国产| 久久6这里有精品| 在线免费十八禁| 桃色一区二区三区在线观看| 99久久人妻综合| 秋霞伦理黄片| 看十八女毛片水多多多| 久久热精品热| 看非洲黑人一级黄片| 长腿黑丝高跟| 午夜亚洲福利在线播放| 亚洲av免费高清在线观看| 久热久热在线精品观看| 国产成人a区在线观看| 国产成人一区二区在线| 亚洲成人av在线免费| 国产亚洲5aaaaa淫片| 我的女老师完整版在线观看| 国产精品日韩av在线免费观看| 欧美激情在线99| 欧美色视频一区免费| 成人鲁丝片一二三区免费| 十八禁国产超污无遮挡网站| 欧美+日韩+精品| 亚洲av日韩在线播放| 波多野结衣高清无吗| 九九久久精品国产亚洲av麻豆| 国产精品爽爽va在线观看网站| 一区二区三区免费毛片| 久久精品夜色国产| 国产乱人偷精品视频| 日韩在线高清观看一区二区三区| 熟女电影av网| 亚洲精品日韩av片在线观看| 国产成年人精品一区二区| 亚洲av一区综合| 亚洲av免费高清在线观看| 中文字幕免费在线视频6| 精品久久久久久久人妻蜜臀av| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩在线中文字幕| 免费av毛片视频| 99九九线精品视频在线观看视频| 2021少妇久久久久久久久久久| 欧美丝袜亚洲另类| 国产大屁股一区二区在线视频| 欧美+日韩+精品| 欧美一区二区精品小视频在线| 精品国产三级普通话版| 淫秽高清视频在线观看| 午夜免费激情av| 极品教师在线视频| 久久精品熟女亚洲av麻豆精品 | 国产午夜精品一二区理论片| 最近中文字幕2019免费版| 国产午夜精品久久久久久一区二区三区| 国产精品国产三级国产av玫瑰| 嫩草影院新地址| 精品人妻熟女av久视频| 日韩欧美国产在线观看| 欧美人与善性xxx| 成人毛片60女人毛片免费| 久久99热这里只有精品18| 国产精品女同一区二区软件| 欧美色视频一区免费| 精品99又大又爽又粗少妇毛片| 青春草视频在线免费观看| 看十八女毛片水多多多| 亚洲成人av在线免费| 亚洲电影在线观看av| 亚洲成色77777| 国产一级毛片七仙女欲春2| 国产精品一区二区在线观看99 | 亚洲av免费高清在线观看| 内射极品少妇av片p| 99久久精品一区二区三区| 久久鲁丝午夜福利片| 女人久久www免费人成看片 | 国产淫片久久久久久久久| 国产精品久久电影中文字幕| 亚洲av中文字字幕乱码综合| 色网站视频免费| 级片在线观看| 欧美变态另类bdsm刘玥| 一级黄色大片毛片| 亚洲成人精品中文字幕电影| 国产精品日韩av在线免费观看| 午夜老司机福利剧场| 国产v大片淫在线免费观看| 国产激情偷乱视频一区二区| 免费av不卡在线播放| 日韩av不卡免费在线播放| 最近最新中文字幕免费大全7| 国产高清国产精品国产三级 | 夜夜爽夜夜爽视频| 成年版毛片免费区| 亚洲自拍偷在线| 国产亚洲午夜精品一区二区久久 | 热99在线观看视频| 久久久精品欧美日韩精品| 亚洲四区av| 亚洲精品乱码久久久v下载方式| 久久亚洲国产成人精品v| 国产极品精品免费视频能看的| 热99在线观看视频| 一夜夜www| 精品无人区乱码1区二区| 国产av不卡久久| 欧美日韩一区二区视频在线观看视频在线 | 久久精品久久精品一区二区三区| 97超碰精品成人国产| 亚洲精品一区蜜桃| h日本视频在线播放| 男人舔女人下体高潮全视频| 国产人妻一区二区三区在| 99热全是精品| 中文天堂在线官网| 中文字幕精品亚洲无线码一区| 免费电影在线观看免费观看| 精品久久久久久电影网 | 免费一级毛片在线播放高清视频| 成年女人看的毛片在线观看| 26uuu在线亚洲综合色| 99九九线精品视频在线观看视频| 最近2019中文字幕mv第一页| 十八禁国产超污无遮挡网站| 在线观看66精品国产| 精华霜和精华液先用哪个| 免费大片18禁| 午夜激情福利司机影院| 大又大粗又爽又黄少妇毛片口| 人体艺术视频欧美日本| 久久久久精品久久久久真实原创| 国产精品精品国产色婷婷| 欧美色视频一区免费| 一级黄片播放器| 看黄色毛片网站| 欧美日韩综合久久久久久| 波多野结衣高清无吗| 亚洲精品乱码久久久久久按摩| 亚洲经典国产精华液单| 亚洲成人中文字幕在线播放| 欧美bdsm另类| 又粗又硬又长又爽又黄的视频| 99久久精品国产国产毛片| 免费观看性生交大片5| 成人二区视频| 欧美色视频一区免费| 春色校园在线视频观看| av又黄又爽大尺度在线免费看 | 午夜精品一区二区三区免费看| 国产91av在线免费观看| 伊人久久精品亚洲午夜| 国内精品宾馆在线| 亚洲av成人精品一区久久| 成人亚洲欧美一区二区av| 国产精品野战在线观看| 身体一侧抽搐| 插逼视频在线观看| 精品一区二区免费观看| 国产一区二区亚洲精品在线观看| 久久人妻av系列| 麻豆乱淫一区二区| 国产一区有黄有色的免费视频 | 在现免费观看毛片| 超碰97精品在线观看| 国产91av在线免费观看| 成人美女网站在线观看视频| 国产精品电影一区二区三区| 久久精品夜色国产| 能在线免费看毛片的网站| 免费黄网站久久成人精品| 久久久久国产网址| 国产av码专区亚洲av| 亚洲国产精品sss在线观看| 成人国产麻豆网| 成人漫画全彩无遮挡| 免费黄色在线免费观看| 99久久精品热视频| 久久久久久久久久黄片| 午夜a级毛片| 久久久久久久久久黄片| 久久久午夜欧美精品| 色综合色国产| 嫩草影院新地址| 亚洲欧美精品专区久久| 亚洲成人精品中文字幕电影| 久久久欧美国产精品| 一区二区三区高清视频在线| 国产精品久久视频播放| 免费看美女性在线毛片视频| 国产精品,欧美在线| 国产高潮美女av| 亚洲一区高清亚洲精品| 熟妇人妻久久中文字幕3abv| 国产午夜精品久久久久久一区二区三区| 美女xxoo啪啪120秒动态图| 国产精品一二三区在线看| 中国美白少妇内射xxxbb| 国产av在哪里看| 欧美性猛交╳xxx乱大交人| 国产精品一二三区在线看| 麻豆成人av视频| 黄片无遮挡物在线观看| 校园人妻丝袜中文字幕| 国产精品蜜桃在线观看| 伦理电影大哥的女人| 哪个播放器可以免费观看大片| 一区二区三区四区激情视频| 精品久久久久久成人av| 欧美日韩综合久久久久久| 免费不卡的大黄色大毛片视频在线观看 | av免费在线看不卡| 一本一本综合久久| 在线播放国产精品三级| 国产成人a区在线观看| 亚洲欧美精品自产自拍| 亚洲国产欧美人成| 91aial.com中文字幕在线观看| 国产成人91sexporn| 青青草视频在线视频观看| 国产真实乱freesex| 欧美成人精品欧美一级黄| 欧美区成人在线视频| 国产欧美另类精品又又久久亚洲欧美| 国产精品.久久久| 久久99热这里只有精品18| 亚洲国产精品成人综合色| 成人午夜高清在线视频| 免费黄网站久久成人精品| 夜夜爽夜夜爽视频| 老司机影院毛片| 边亲边吃奶的免费视频| 中文字幕熟女人妻在线| or卡值多少钱| 午夜福利成人在线免费观看| 亚洲,欧美,日韩| 国产在线男女| 亚洲自拍偷在线| 我要看日韩黄色一级片| 老司机福利观看| 老司机影院毛片| 高清日韩中文字幕在线| 日韩成人伦理影院| 亚洲乱码一区二区免费版| 欧美一区二区亚洲| 免费看日本二区| 99热这里只有是精品50| 日韩一本色道免费dvd| 亚洲激情五月婷婷啪啪| 国产男人的电影天堂91| 少妇熟女欧美另类| 亚洲欧美清纯卡通| 免费一级毛片在线播放高清视频| 欧美日韩综合久久久久久| 97在线视频观看| 国产精品,欧美在线| 最近2019中文字幕mv第一页| 午夜福利视频1000在线观看| 中文字幕av在线有码专区| 成年女人看的毛片在线观看| 久久精品综合一区二区三区| 国产精品.久久久| 天堂√8在线中文| 久久久久网色| 久久精品熟女亚洲av麻豆精品 | 亚洲欧美清纯卡通| 中文字幕人妻熟人妻熟丝袜美| 麻豆久久精品国产亚洲av| 草草在线视频免费看| 99热网站在线观看| 国产成人精品婷婷| 三级国产精品片| 中文字幕制服av| 91精品一卡2卡3卡4卡| 亚洲精品日韩在线中文字幕| 99国产精品一区二区蜜桃av| 国产黄色视频一区二区在线观看 | 美女xxoo啪啪120秒动态图| 欧美成人精品欧美一级黄| 少妇的逼水好多| 国产女主播在线喷水免费视频网站 | 日韩一区二区视频免费看| 中文欧美无线码| 久久午夜福利片| 老司机福利观看| 亚洲av二区三区四区| 91aial.com中文字幕在线观看| 插阴视频在线观看视频| 亚洲av成人精品一区久久| 女人久久www免费人成看片 | 久热久热在线精品观看| 日日摸夜夜添夜夜添av毛片| 欧美bdsm另类| 久久国内精品自在自线图片| 亚洲av电影不卡..在线观看| 国产一级毛片七仙女欲春2| 一级二级三级毛片免费看| 高清毛片免费看| 精品少妇黑人巨大在线播放 | 日韩 亚洲 欧美在线| 观看免费一级毛片| 91精品一卡2卡3卡4卡| 九九久久精品国产亚洲av麻豆| 久99久视频精品免费| 美女被艹到高潮喷水动态| 免费观看a级毛片全部| 18禁裸乳无遮挡免费网站照片| 韩国av在线不卡| av福利片在线观看| 国产精品爽爽va在线观看网站| 麻豆乱淫一区二区| 中文乱码字字幕精品一区二区三区 | 51国产日韩欧美| 亚洲图色成人| 熟女电影av网| 免费大片18禁| 狂野欧美白嫩少妇大欣赏| av天堂中文字幕网| 一级毛片电影观看 | 欧美成人午夜免费资源| 1000部很黄的大片| 一个人免费在线观看电影| 国产69精品久久久久777片| av.在线天堂| 成人漫画全彩无遮挡| 99热全是精品| 人体艺术视频欧美日本| 少妇裸体淫交视频免费看高清| 色视频www国产| 国产黄片美女视频| 亚州av有码| 91精品国产九色| 亚洲经典国产精华液单| 青春草亚洲视频在线观看| 内地一区二区视频在线| 欧美又色又爽又黄视频| 男女边吃奶边做爰视频| 亚洲美女搞黄在线观看| 美女脱内裤让男人舔精品视频| 最近视频中文字幕2019在线8| 国产精华一区二区三区| 日韩国内少妇激情av| 国产精品一及| 日日干狠狠操夜夜爽| 欧美变态另类bdsm刘玥| 亚洲真实伦在线观看| 国产在视频线精品| 国产一级毛片在线| 永久免费av网站大全| 97在线视频观看| 国产乱人偷精品视频| 国产一区有黄有色的免费视频 | 国产午夜福利久久久久久| 青春草亚洲视频在线观看| 亚洲人成网站高清观看| www日本黄色视频网| 狂野欧美激情性xxxx在线观看| 国产单亲对白刺激| 寂寞人妻少妇视频99o| 狂野欧美激情性xxxx在线观看| 久久久久九九精品影院| 欧美人与善性xxx| 91午夜精品亚洲一区二区三区| 九九热线精品视视频播放| 桃色一区二区三区在线观看| 中文亚洲av片在线观看爽| 精品午夜福利在线看| 中文字幕亚洲精品专区| 日韩一本色道免费dvd| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 精品熟女少妇av免费看| 欧美激情在线99| 日韩欧美在线乱码| 精华霜和精华液先用哪个| 国产乱来视频区| 国产精品.久久久| 噜噜噜噜噜久久久久久91| 在线免费观看不下载黄p国产| 高清午夜精品一区二区三区| 又粗又硬又长又爽又黄的视频| 别揉我奶头 嗯啊视频| av福利片在线观看| 免费观看性生交大片5| 精品一区二区三区视频在线| 亚洲人成网站在线观看播放| 黑人高潮一二区| 久久久久久久久久久丰满| 日韩欧美精品v在线| 日本爱情动作片www.在线观看| a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 女的被弄到高潮叫床怎么办| 永久免费av网站大全| av免费在线看不卡| 亚洲美女搞黄在线观看| 亚洲精品日韩av片在线观看| 听说在线观看完整版免费高清| 成人欧美大片| 国产真实伦视频高清在线观看| 午夜激情福利司机影院| 国产乱人偷精品视频| 亚洲乱码一区二区免费版| 亚洲av熟女| 中文字幕亚洲精品专区| 搞女人的毛片| АⅤ资源中文在线天堂| 精品人妻一区二区三区麻豆| 欧美成人精品欧美一级黄| 三级国产精品欧美在线观看| 少妇猛男粗大的猛烈进出视频 | 丰满乱子伦码专区| 99久国产av精品| 老司机影院毛片| 嫩草影院新地址| 久久久久免费精品人妻一区二区| 日韩在线高清观看一区二区三区| 久久综合国产亚洲精品| 日韩精品青青久久久久久| 老司机影院成人| 色网站视频免费| 免费大片18禁| 日韩一区二区视频免费看| 国产欧美另类精品又又久久亚洲欧美| 日本黄色片子视频| 午夜福利在线在线| 日本-黄色视频高清免费观看| 一边摸一边抽搐一进一小说| www.av在线官网国产| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 亚洲欧美成人精品一区二区| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 99久久精品热视频| 国产亚洲精品久久久com| 最近中文字幕2019免费版| 国产精品日韩av在线免费观看| 九九爱精品视频在线观看| 伊人久久精品亚洲午夜| 少妇裸体淫交视频免费看高清| 中文字幕久久专区| 长腿黑丝高跟| 欧美成人精品欧美一级黄| 一级黄色大片毛片| 精品久久久久久久末码| 一级黄色大片毛片| АⅤ资源中文在线天堂| 免费观看精品视频网站| 亚洲av福利一区| a级毛色黄片| 日产精品乱码卡一卡2卡三| 99久久无色码亚洲精品果冻| 少妇熟女aⅴ在线视频| 校园人妻丝袜中文字幕| 国产成人一区二区在线| 热99在线观看视频| 国产精品嫩草影院av在线观看| 国产探花极品一区二区| 一个人免费在线观看电影| 成人一区二区视频在线观看| ponron亚洲| 国产精品美女特级片免费视频播放器| av在线蜜桃| 好男人视频免费观看在线| 欧美最新免费一区二区三区| 久久精品人妻少妇| 精品无人区乱码1区二区| 汤姆久久久久久久影院中文字幕 | 久久99热这里只频精品6学生 | 欧美不卡视频在线免费观看| 国产私拍福利视频在线观看| 91精品国产九色| 国产一区二区在线观看日韩| 精品午夜福利在线看| 欧美97在线视频| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽人人片av| 大香蕉久久网| 国产亚洲5aaaaa淫片| 亚洲一区高清亚洲精品| 国产 一区精品| 亚洲人与动物交配视频| 大话2 男鬼变身卡| 亚洲人与动物交配视频| 波野结衣二区三区在线| 国产精品久久久久久av不卡| 亚洲精品一区蜜桃| 国内少妇人妻偷人精品xxx网站| 午夜激情欧美在线| 你懂的网址亚洲精品在线观看 | kizo精华| 听说在线观看完整版免费高清| 伦精品一区二区三区| 午夜亚洲福利在线播放| 十八禁国产超污无遮挡网站| 久久人人爽人人爽人人片va| 啦啦啦韩国在线观看视频| av免费在线看不卡| 亚洲在线自拍视频| 欧美一区二区精品小视频在线| 亚洲最大成人手机在线| av在线蜜桃| 精品欧美国产一区二区三| 亚洲国产高清在线一区二区三| 男女国产视频网站| 美女被艹到高潮喷水动态| 国国产精品蜜臀av免费| 中文资源天堂在线| 狠狠狠狠99中文字幕| 99热这里只有精品一区| 亚洲国产精品专区欧美| 国产亚洲av片在线观看秒播厂 | 久久精品91蜜桃| 好男人在线观看高清免费视频| 国产一级毛片七仙女欲春2| 一本一本综合久久| 亚洲精品亚洲一区二区| 久久精品国产99精品国产亚洲性色| 久久精品国产鲁丝片午夜精品| 国产精品av视频在线免费观看| 乱码一卡2卡4卡精品| 男人舔女人下体高潮全视频| 亚洲av不卡在线观看| 久久久久国产网址| 欧美xxxx黑人xx丫x性爽| 国产午夜福利久久久久久| 亚洲av熟女| 国产熟女欧美一区二区| 色播亚洲综合网| 熟女电影av网| 亚洲成人久久爱视频| 午夜日本视频在线| 最近2019中文字幕mv第一页| 欧美日本亚洲视频在线播放| 美女国产视频在线观看| 中文字幕久久专区| 99九九线精品视频在线观看视频| 少妇熟女aⅴ在线视频| 综合色av麻豆| 欧美激情久久久久久爽电影| 春色校园在线视频观看| 在线a可以看的网站| 久久久国产成人免费| 亚洲成人久久爱视频| 欧美一区二区精品小视频在线| 美女cb高潮喷水在线观看| 午夜视频国产福利| 国产精品女同一区二区软件| 久久这里有精品视频免费| 日韩欧美精品v在线| 国产69精品久久久久777片| 看片在线看免费视频| 免费大片18禁| 三级男女做爰猛烈吃奶摸视频| 一个人观看的视频www高清免费观看| 成人高潮视频无遮挡免费网站| 国产精品久久久久久av不卡| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 国产伦理片在线播放av一区| 亚洲国产精品成人综合色| 精品熟女少妇av免费看| 一级黄片播放器| 午夜福利在线观看免费完整高清在| 午夜激情福利司机影院| 久久精品国产自在天天线| 亚洲av电影不卡..在线观看| 99久久中文字幕三级久久日本| 乱码一卡2卡4卡精品| 成人特级av手机在线观看| 国产在线男女| 2021天堂中文幕一二区在线观| 国产精品永久免费网站| 久久人人爽人人爽人人片va| 秋霞伦理黄片| 69人妻影院| 高清视频免费观看一区二区 | 少妇猛男粗大的猛烈进出视频 | 天天躁日日操中文字幕| 国产免费福利视频在线观看| 乱系列少妇在线播放| 中文亚洲av片在线观看爽| 在线免费观看不下载黄p国产| 国产精品无大码| 免费黄网站久久成人精品| 国产成人精品婷婷| 国产免费男女视频| 日本免费a在线|