• <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白浆流出| 黄色成人免费大全| 男女高潮啪啪啪动态图| 欧美人与性动交α欧美精品济南到| 久久精品国产亚洲av香蕉五月| 国产精品98久久久久久宅男小说| 欧美成人免费av一区二区三区| 99riav亚洲国产免费| 久久久久久大精品| 精品人妻1区二区| 久久久久国产精品人妻aⅴ院| 欧美一级毛片孕妇| 日日干狠狠操夜夜爽| 国产av一区二区精品久久| 变态另类成人亚洲欧美熟女 | 国产xxxxx性猛交| 国产伦一二天堂av在线观看| 久久人妻av系列| 国产成人影院久久av| 欧美日本亚洲视频在线播放| 一级a爱视频在线免费观看| 久久久久久亚洲精品国产蜜桃av| 久久久水蜜桃国产精品网| 高清欧美精品videossex| 国产亚洲精品综合一区在线观看 | 国产亚洲精品久久久久久毛片| 夫妻午夜视频| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲综合一区二区三区_| 亚洲第一av免费看| 久久久久久免费高清国产稀缺| 成人永久免费在线观看视频| 搡老乐熟女国产| 日本vs欧美在线观看视频| 欧美日韩视频精品一区| 国产精品爽爽va在线观看网站 | 人人妻,人人澡人人爽秒播| 亚洲一区中文字幕在线| 欧美乱码精品一区二区三区| 国产精品一区二区在线不卡| 夜夜夜夜夜久久久久| 男人操女人黄网站| 日本a在线网址| 一级黄色大片毛片| 在线观看日韩欧美| 日韩高清综合在线| 一级,二级,三级黄色视频| 俄罗斯特黄特色一大片| 男女床上黄色一级片免费看| 亚洲一区二区三区色噜噜 | 国产精品久久久久成人av| 欧美最黄视频在线播放免费 | 国产精品一区二区免费欧美| 中文字幕色久视频| av超薄肉色丝袜交足视频| 成人免费观看视频高清| 国产午夜精品久久久久久| 男女床上黄色一级片免费看| 国产一卡二卡三卡精品| 黄色 视频免费看| 久久热在线av| 久久久精品国产亚洲av高清涩受| a在线观看视频网站| 99国产精品一区二区蜜桃av| 长腿黑丝高跟| 国产av一区二区精品久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品91蜜桃| 热99国产精品久久久久久7| 精品高清国产在线一区| 级片在线观看| 在线av久久热| 欧美丝袜亚洲另类 | 99精国产麻豆久久婷婷| 国产91精品成人一区二区三区| 亚洲少妇的诱惑av| 日本黄色视频三级网站网址| 热re99久久精品国产66热6| 亚洲va日本ⅴa欧美va伊人久久| av电影中文网址| 热99国产精品久久久久久7| 一级毛片女人18水好多| 久久精品91无色码中文字幕| 视频区欧美日本亚洲| 国产又爽黄色视频| 亚洲少妇的诱惑av| 国产深夜福利视频在线观看| 黄色 视频免费看| 9热在线视频观看99| 99久久99久久久精品蜜桃| 啦啦啦在线免费观看视频4| 夜夜躁狠狠躁天天躁| 亚洲欧美激情综合另类| 午夜精品久久久久久毛片777| 热99re8久久精品国产| 在线播放国产精品三级| 黄色a级毛片大全视频| 91在线观看av| 亚洲欧美日韩另类电影网站| 欧美日韩精品网址| 亚洲色图av天堂| 热re99久久国产66热| 欧美另类亚洲清纯唯美| 熟女少妇亚洲综合色aaa.| 欧美一级毛片孕妇| 老司机午夜福利在线观看视频| 一级片'在线观看视频| 国产亚洲欧美98| 日本撒尿小便嘘嘘汇集6| 亚洲人成电影免费在线| 真人一进一出gif抽搐免费| 视频在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 九色亚洲精品在线播放| 在线国产一区二区在线| 色综合婷婷激情| 看黄色毛片网站| 亚洲中文日韩欧美视频| 中文字幕色久视频| 久久久久久久久免费视频了| 亚洲av熟女| 欧美成人性av电影在线观看| 亚洲人成电影观看| 国产麻豆69| 中文字幕人妻熟女乱码| 欧美日韩中文字幕国产精品一区二区三区 | 高清欧美精品videossex| 又黄又粗又硬又大视频| 精品国产亚洲在线| 国产亚洲精品第一综合不卡| 欧美午夜高清在线| 国产亚洲精品久久久久5区| 超色免费av| 嫁个100分男人电影在线观看| 久久精品亚洲av国产电影网| 国产精品成人在线| 国产av又大| 久久九九热精品免费| 亚洲色图 男人天堂 中文字幕| 中文欧美无线码| 两性夫妻黄色片| 欧美精品一区二区免费开放| 亚洲成a人片在线一区二区| 欧美日韩亚洲国产一区二区在线观看| 亚洲成国产人片在线观看| 久久久精品国产亚洲av高清涩受| 老司机午夜十八禁免费视频| 久久草成人影院| 色综合欧美亚洲国产小说| √禁漫天堂资源中文www| 成人三级做爰电影| 母亲3免费完整高清在线观看| 一边摸一边抽搐一进一小说| 亚洲精华国产精华精| 国产成年人精品一区二区 | 19禁男女啪啪无遮挡网站| 欧美老熟妇乱子伦牲交| 美女高潮喷水抽搐中文字幕| 性色av乱码一区二区三区2| 国产成人影院久久av| 欧美一级毛片孕妇| 色综合欧美亚洲国产小说| 一级片'在线观看视频| 久久久久精品国产欧美久久久| 黑人猛操日本美女一级片| 夫妻午夜视频| 久久影院123| 人人妻人人添人人爽欧美一区卜| 人妻丰满熟妇av一区二区三区| 国产激情欧美一区二区| 看片在线看免费视频| 亚洲午夜理论影院| 老司机亚洲免费影院| 精品少妇一区二区三区视频日本电影| 手机成人av网站| 麻豆成人av在线观看| 在线观看www视频免费| 天天躁夜夜躁狠狠躁躁| 精品福利永久在线观看| 久久人妻福利社区极品人妻图片| 日韩av在线大香蕉| 亚洲一区二区三区欧美精品| 国产色视频综合| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 欧美中文综合在线视频| 精品国内亚洲2022精品成人| 久久久久久久久中文| netflix在线观看网站| 欧美日韩瑟瑟在线播放| 自线自在国产av| 男人舔女人下体高潮全视频| 在线永久观看黄色视频| 99精品欧美一区二区三区四区| 亚洲九九香蕉| 欧美精品一区二区免费开放| 精品国产一区二区久久| 视频区图区小说| 亚洲欧洲精品一区二区精品久久久| 最新美女视频免费是黄的| 真人做人爱边吃奶动态| 一级毛片高清免费大全| 老司机亚洲免费影院| 国产亚洲av高清不卡| 国产精品国产高清国产av| 精品福利观看| 色老头精品视频在线观看| 欧美成人午夜精品| 纯流量卡能插随身wifi吗| 女人高潮潮喷娇喘18禁视频| 美女午夜性视频免费| 黄色 视频免费看| 久久人人爽av亚洲精品天堂| 久久久久精品国产欧美久久久| 在线国产一区二区在线| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品久久成人aⅴ小说| 精品国产一区二区三区四区第35| 岛国视频午夜一区免费看| 露出奶头的视频| 满18在线观看网站| 中文欧美无线码| 国产精品成人在线| 青草久久国产| 日韩一卡2卡3卡4卡2021年| 亚洲免费av在线视频| 亚洲情色 制服丝袜| 好看av亚洲va欧美ⅴa在| 在线十欧美十亚洲十日本专区| 成年版毛片免费区| 每晚都被弄得嗷嗷叫到高潮| 国产又爽黄色视频| 99国产精品免费福利视频| 亚洲性夜色夜夜综合| 久久久久九九精品影院| 亚洲欧美激情综合另类| 国产亚洲欧美精品永久| www.熟女人妻精品国产| 50天的宝宝边吃奶边哭怎么回事| 一级毛片精品| 亚洲 欧美一区二区三区| 人人妻人人添人人爽欧美一区卜| 怎么达到女性高潮| 女人爽到高潮嗷嗷叫在线视频| 亚洲中文字幕日韩| 国产又爽黄色视频| 无遮挡黄片免费观看| 亚洲九九香蕉| 日本wwww免费看| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 国产av一区二区精品久久| 亚洲国产欧美一区二区综合| 夫妻午夜视频| 国产乱人伦免费视频| 国产精品久久久av美女十八| 国内毛片毛片毛片毛片毛片| 在线观看免费日韩欧美大片| 精品一区二区三区av网在线观看| 欧美精品亚洲一区二区| 国产又色又爽无遮挡免费看| 在线看a的网站| 亚洲人成伊人成综合网2020| 在线国产一区二区在线| 老汉色∧v一级毛片| 免费看a级黄色片| 亚洲精品美女久久av网站| 黄片大片在线免费观看| 欧美中文日本在线观看视频| 国产精品野战在线观看 | 精品国产一区二区三区四区第35| 一二三四社区在线视频社区8| 一夜夜www| 别揉我奶头~嗯~啊~动态视频| 国产aⅴ精品一区二区三区波| 国产精品乱码一区二三区的特点 | 99久久久亚洲精品蜜臀av| 午夜福利一区二区在线看| 精品免费久久久久久久清纯| 免费在线观看黄色视频的| 91老司机精品| 亚洲欧美精品综合一区二区三区| 欧美在线黄色| 国产精品九九99| 成人黄色视频免费在线看| 亚洲免费av在线视频| 欧美 亚洲 国产 日韩一| 午夜免费观看网址| 日本撒尿小便嘘嘘汇集6| 午夜激情av网站| 日本精品一区二区三区蜜桃| 满18在线观看网站| 极品人妻少妇av视频| 男男h啪啪无遮挡| 少妇粗大呻吟视频| www.999成人在线观看| 亚洲五月婷婷丁香| 淫妇啪啪啪对白视频| 成人国产一区最新在线观看| 女性生殖器流出的白浆| 欧美日韩福利视频一区二区| 最近最新中文字幕大全电影3 | 精品久久久久久,| 国产一区二区三区在线臀色熟女 | 在线观看66精品国产| 午夜a级毛片| 大型av网站在线播放| 少妇的丰满在线观看| 欧美日韩国产mv在线观看视频| 日本五十路高清| 午夜福利免费观看在线| 免费在线观看黄色视频的| 91精品国产国语对白视频| 水蜜桃什么品种好| 黄色成人免费大全| 国产xxxxx性猛交| 少妇裸体淫交视频免费看高清 | 夜夜看夜夜爽夜夜摸 | 久久精品影院6| 国产亚洲精品久久久久5区| 香蕉国产在线看| www.精华液| 亚洲午夜精品一区,二区,三区| 国产一区在线观看成人免费| 99精国产麻豆久久婷婷| 一进一出抽搐动态| 国产精品 欧美亚洲| 一本大道久久a久久精品| 99国产极品粉嫩在线观看| 在线免费观看的www视频| a在线观看视频网站| 欧美另类亚洲清纯唯美| 免费在线观看视频国产中文字幕亚洲| 免费看a级黄色片| 亚洲七黄色美女视频| 性欧美人与动物交配| 老司机亚洲免费影院| av欧美777| 亚洲一区二区三区欧美精品| 91大片在线观看| 动漫黄色视频在线观看| 人人妻人人爽人人添夜夜欢视频| 大香蕉久久成人网| 国产激情欧美一区二区| 日韩免费av在线播放| 视频区欧美日本亚洲| aaaaa片日本免费| 搡老熟女国产l中国老女人| 国产男靠女视频免费网站| av在线天堂中文字幕 | 成熟少妇高潮喷水视频| 又黄又粗又硬又大视频| 人人妻,人人澡人人爽秒播| 淫秽高清视频在线观看| 少妇的丰满在线观看| 97人妻天天添夜夜摸| 久久精品国产清高在天天线| 色精品久久人妻99蜜桃| 久久伊人香网站| 女同久久另类99精品国产91| 日韩精品免费视频一区二区三区| 日本a在线网址| 国产精品亚洲av一区麻豆| 国产精品乱码一区二三区的特点 | 天天躁夜夜躁狠狠躁躁| 极品教师在线免费播放| 女人被狂操c到高潮| 欧美日韩瑟瑟在线播放| 国产精品99久久99久久久不卡| 国产真人三级小视频在线观看| 99精国产麻豆久久婷婷| av网站在线播放免费| 九色亚洲精品在线播放| 80岁老熟妇乱子伦牲交| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 嫩草影视91久久| 国产精品av久久久久免费| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品合色在线| 久久狼人影院| 欧美成狂野欧美在线观看| 91在线观看av| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 激情在线观看视频在线高清| www.999成人在线观看| 国产欧美日韩综合在线一区二区| 午夜福利影视在线免费观看| 久久狼人影院| 桃红色精品国产亚洲av| 免费在线观看日本一区| 久久久久亚洲av毛片大全| 黄频高清免费视频| 欧美中文日本在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 久久久久久久久免费视频了| 午夜精品国产一区二区电影| 日韩视频一区二区在线观看| 97超级碰碰碰精品色视频在线观看| 午夜两性在线视频| 亚洲第一av免费看| 制服诱惑二区| 亚洲午夜精品一区,二区,三区| 久久久水蜜桃国产精品网| 美女扒开内裤让男人捅视频| 黑人巨大精品欧美一区二区蜜桃| 久久久久国内视频| 黑人操中国人逼视频| 88av欧美| 美女扒开内裤让男人捅视频| 黑人巨大精品欧美一区二区蜜桃| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| av在线播放免费不卡| 成人黄色视频免费在线看| 亚洲精品一卡2卡三卡4卡5卡| 国产av精品麻豆| 91老司机精品| 国产高清激情床上av| 亚洲免费av在线视频| 日本免费一区二区三区高清不卡 | 国产精品二区激情视频| 免费久久久久久久精品成人欧美视频| 亚洲片人在线观看| 欧美大码av| 亚洲av美国av| 日日摸夜夜添夜夜添小说| 欧美日韩中文字幕国产精品一区二区三区 | 日本vs欧美在线观看视频| 久久人人精品亚洲av| 色综合欧美亚洲国产小说| 80岁老熟妇乱子伦牲交| 久久人人精品亚洲av| 日日摸夜夜添夜夜添小说| 久久久久国产精品人妻aⅴ院| 亚洲欧美日韩另类电影网站| 大型黄色视频在线免费观看| 国产在线观看jvid| 亚洲欧美日韩无卡精品| 丁香欧美五月| 日韩三级视频一区二区三区| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 免费看十八禁软件| 老熟妇乱子伦视频在线观看| 久久亚洲真实| 亚洲专区国产一区二区| 亚洲国产精品合色在线| 老司机在亚洲福利影院| 巨乳人妻的诱惑在线观看| 国产成人影院久久av| 精品久久久久久电影网| 国产在线精品亚洲第一网站| 欧美性长视频在线观看| 日本 av在线| 日本a在线网址| 午夜精品在线福利| 亚洲 欧美一区二区三区| 国产又爽黄色视频| tocl精华| 级片在线观看| 欧美乱色亚洲激情| 亚洲三区欧美一区| 久久久久久久午夜电影 | 首页视频小说图片口味搜索| 国产精品香港三级国产av潘金莲| 国产日韩一区二区三区精品不卡| 美女高潮到喷水免费观看| 午夜91福利影院| 天堂中文最新版在线下载| 精品福利观看| 久久午夜亚洲精品久久| 国产av在哪里看| 老司机福利观看| 亚洲视频免费观看视频| 丁香六月欧美| 欧美日韩乱码在线| 黑人巨大精品欧美一区二区mp4| 欧美色视频一区免费| 人人妻,人人澡人人爽秒播| 成年女人毛片免费观看观看9| 亚洲精品久久成人aⅴ小说| 成人18禁高潮啪啪吃奶动态图| 天堂√8在线中文| 欧美乱码精品一区二区三区| 国产精品久久久久久人妻精品电影| 国产精品亚洲一级av第二区| 国产国语露脸激情在线看| 亚洲第一av免费看| 久久国产精品人妻蜜桃| 黑丝袜美女国产一区| 欧美日韩av久久| 成人手机av| 欧美+亚洲+日韩+国产| 精品欧美一区二区三区在线| 欧美老熟妇乱子伦牲交| 真人做人爱边吃奶动态| 日本a在线网址| 亚洲 欧美 日韩 在线 免费| 国产麻豆69| 精品一品国产午夜福利视频| 亚洲色图综合在线观看| 天堂中文最新版在线下载| 亚洲三区欧美一区| 两人在一起打扑克的视频| 国产高清激情床上av| 精品一区二区三区av网在线观看| 成在线人永久免费视频| 水蜜桃什么品种好| 99精品久久久久人妻精品| 国产成人系列免费观看| 久久人妻av系列| 国产成人精品在线电影| 国产乱人伦免费视频| 真人一进一出gif抽搐免费| 最好的美女福利视频网| svipshipincom国产片| 老汉色av国产亚洲站长工具| 免费在线观看视频国产中文字幕亚洲| 丝袜人妻中文字幕| 日韩欧美国产一区二区入口| 亚洲国产中文字幕在线视频| 99在线人妻在线中文字幕| 国产成+人综合+亚洲专区| a在线观看视频网站| 国产深夜福利视频在线观看| 国产激情久久老熟女| 亚洲av成人av| www日本在线高清视频| 久久性视频一级片| 国产三级在线视频| 伦理电影免费视频| 日韩视频一区二区在线观看| 国产精品免费视频内射| 免费看十八禁软件| 亚洲精品一区av在线观看| 免费av中文字幕在线| 黄色视频不卡| 国产乱人伦免费视频| 18禁黄网站禁片午夜丰满| 精品国产一区二区三区四区第35| 国产成人av激情在线播放| 欧美激情高清一区二区三区| 天堂俺去俺来也www色官网| 99久久精品国产亚洲精品| 国产精品国产av在线观看| 国产亚洲精品综合一区在线观看 | 色哟哟哟哟哟哟| 国产成人欧美在线观看| 99热只有精品国产| 一本综合久久免费| 91国产中文字幕| 欧美日本亚洲视频在线播放| 久久午夜综合久久蜜桃| 香蕉丝袜av| 国产99白浆流出| 精品欧美一区二区三区在线| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一出视频| 欧美成狂野欧美在线观看| 国产成人系列免费观看| 亚洲 欧美一区二区三区| 精品福利永久在线观看| 又紧又爽又黄一区二区| 国产成人影院久久av| 欧美一区二区精品小视频在线| 校园春色视频在线观看| 日本精品一区二区三区蜜桃| 黄色怎么调成土黄色| 亚洲久久久国产精品| 国产成+人综合+亚洲专区| 男人操女人黄网站| 国产又色又爽无遮挡免费看| 一级毛片精品| av网站免费在线观看视频| 欧美+亚洲+日韩+国产| 级片在线观看| 日本免费a在线| 不卡一级毛片| 人妻久久中文字幕网| 国产亚洲欧美在线一区二区| 一本综合久久免费| 欧美在线一区亚洲| 亚洲欧美精品综合一区二区三区| 丝袜人妻中文字幕| 国产成人一区二区三区免费视频网站| 涩涩av久久男人的天堂| 人妻久久中文字幕网| 亚洲激情在线av| 高潮久久久久久久久久久不卡| 亚洲av电影在线进入| 99国产精品一区二区三区| 亚洲一区二区三区欧美精品| 一区二区三区激情视频| 日韩有码中文字幕| 亚洲成人免费av在线播放| 国产精品 欧美亚洲| 久久香蕉国产精品| 看黄色毛片网站| 99国产综合亚洲精品| 九色亚洲精品在线播放| 一区二区三区精品91| 大型黄色视频在线免费观看| 成人国语在线视频| 久久狼人影院| 国产日韩一区二区三区精品不卡| 亚洲av日韩精品久久久久久密| 丝袜美足系列| 一夜夜www| 成年版毛片免费区| 一区二区三区精品91| 波多野结衣高清无吗|