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

    基于MODIS-LAI 數(shù)據(jù)的廣西甘蔗物候期提取

    2021-03-01 09:11:36謝鑫昌楊云川田憶廖麗萍韋鈞培周津羽陳立華
    關(guān)鍵詞:物候時序甘蔗

    謝鑫昌 ,楊云川,2,3*,田憶,廖麗萍,2,3,韋鈞培,周津羽,陳立華,2,3

    (1.廣西大學(xué)土木建筑工程學(xué)院,廣西 南寧 530004;2.工程防災(zāi)與結(jié)構(gòu)安全教育部重點實驗室,廣西 南寧 530004;3.廣西防災(zāi)減災(zāi)與工程安全重點實驗室,廣西 南寧 530004)

    物候是指作物受氣候變化及各種環(huán)境因素的影響而呈現(xiàn)出以年為周期的自然現(xiàn)象,當(dāng)作物達(dá)到關(guān)鍵的生長發(fā)育狀態(tài)所對應(yīng)的時期即為作物物候期[1]。區(qū)域作物物候的動態(tài)變化是響應(yīng)氣候環(huán)境變化的重要指標(biāo)[2],即反映了作物的生長和發(fā)育規(guī)律,也影響著作物的產(chǎn)量和質(zhì)量[3]。因此,及時、準(zhǔn)確獲取作物物候信息,對作物農(nóng)情實時監(jiān)測、田間精準(zhǔn)管理、作物產(chǎn)量保障和氣候響應(yīng)研究等均具有重要意義[4-5]。

    目前,地面監(jiān)測和衛(wèi)星遙感是及時獲取作物物候信息的主要技術(shù)手段[6-7]。地面監(jiān)測主要是在田間設(shè)立觀測站點或手持光譜儀定時獲取作物物候信息,該方式較為費時費力且局限于小范圍田間管理[7];相比而言,衛(wèi)星遙感方式具有低成本、高時效、觀測區(qū)域廣及省時省力等多方面的顯著優(yōu)勢[8-9],已廣泛應(yīng)用于農(nóng)業(yè)監(jiān)測;遙感監(jiān)測常通過植被指數(shù)或作物理化參數(shù)來反映大空間尺度上作物的光合作用狀態(tài)及長勢季節(jié)變化,最終實現(xiàn)作物物候期的過程識別[10]。

    常用的遙感物候監(jiān)測指標(biāo)有歸一化植被指數(shù)(NDVI,Normalized differential vegetation index)、增強植被指數(shù)(EVI,Enhanced vegetation index)及葉面積指數(shù)(LAI,Leaf area index)等。近年代表性的研究有:辛景峰等[11]、鄭小波等[12]采用的條件時間內(nèi)差法重構(gòu)AVHRR、MODIS-NDVI 時序數(shù)據(jù)集,實現(xiàn)了冬小麥、水稻的關(guān)鍵物候期估測,其物候監(jiān)測誤差均小于1 旬。王堯等[13]、張喜旺等[7]、Chu等[14]、Liu 等[15]及Xu 等[5]在上述研究基礎(chǔ)上采用S-G 濾波法重構(gòu)MODIS-NDVI、HJ-NDVI 時序數(shù)據(jù)集,并通過動態(tài)閾值法提取了單雙季稻、冬小麥、夏玉米、甜菜等多種作物的物候期,結(jié)果表明作物在空間上的物候分布與實地研究具有較好的一致性。上述方法不僅減少云霧對遙感數(shù)據(jù)質(zhì)量的影響,而且降低了數(shù)據(jù)噪聲,能有效提高作物物候信息監(jiān)測的準(zhǔn)確性。楊琳等[1]進(jìn)一步采用非高斯對稱、雙Logistic 擬合法對MODIS-NDVI 時序曲線進(jìn)行平滑降噪處理,實現(xiàn)了冬小麥關(guān)鍵物候期的提取,并通過對比論證了這兩種方法與S-G 濾波方法的優(yōu)劣。李艷等[3]則采用改進(jìn)的最大值合成法重構(gòu)MODIS-NDVI 時序數(shù)據(jù),并進(jìn)一步通過S-G 濾波降噪處理、Logistic 函數(shù)擬合,最后通過曲率法、動態(tài)閾值法實現(xiàn)夏玉米關(guān)鍵物候期提取。凌洋等[16]則采用HANTS、TSF 濾波法重構(gòu)MODIS-EVI,提取了水稻的關(guān)鍵物候期,結(jié)果表明TSF 濾波法準(zhǔn)確度高于HANTS 濾波,并將平滑和濾波處理綜合,充分考慮作物生長的周期性規(guī)律和遙感數(shù)據(jù)自身特點。黃健熙等[4]基于S-G 濾波法重構(gòu)MODISLAI,并結(jié)合雙Logistic 函數(shù)擬合、歷史積溫法實現(xiàn)了冬小麥返青、拔節(jié)及開花期的提取,研究表明雙Logistic 函數(shù)擬合精度較高。Huang 等[17]則基于重構(gòu)的MODIS-NDVI 時序數(shù)據(jù),采用有效積溫法,實現(xiàn)了華北平原冬小麥抽穗期和開花期的提取,其中RMSE 均保持在4.5~6.5 d。上述研究表明,通過融入有效積溫、太陽輻射等氣象因素,能進(jìn)一步提高作物不同關(guān)鍵物候期的識別精度。

    綜上所述,可通過多種時序重構(gòu)法的對比分析來篩選具體區(qū)域、特定作物物候信息提取的最優(yōu)方法;但已有研究主要應(yīng)用于小麥、玉米、水稻等糧食作物,針對甘蔗等經(jīng)濟(jì)作物的物候期識別研究相對較少。廣西屬亞熱帶季風(fēng)氣候,水熱條件充沛,是我國最大的甘蔗產(chǎn)區(qū)和蔗糖產(chǎn)業(yè)基地[18]。然而,廣西多山地丘陵、云霧覆蓋度高、農(nóng)田地塊破碎、插花種植普遍[19],導(dǎo)致區(qū)域甘蔗種植面積雖廣但平均單位面積產(chǎn)量不高、歷年災(zāi)害損失嚴(yán)重。因此,及時、準(zhǔn)確獲取廣西甘蔗的物候信息,是提升甘蔗品質(zhì)和產(chǎn)量、優(yōu)化甘蔗種植結(jié)構(gòu)、統(tǒng)籌機(jī)械化收割等精準(zhǔn)管理的迫切需求。本文即采用覆蓋尺度廣、時間分辨率高的MODIS-LAI 數(shù)據(jù)集[4,8],選擇S-G濾波函數(shù)、非對稱高斯函數(shù)及雙Logistic 函數(shù)等方法進(jìn)行廣西甘蔗物候的LAI 時序擬合,并通過動態(tài)閾值方法實現(xiàn)廣西甘蔗關(guān)鍵物候期提取,對比論證各方法的識別精度差異,為廣西甘蔗種植的智慧管理及蔗糖產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化提供科學(xué)支撐。

    1 數(shù)據(jù)來源與研究方法

    1.1 數(shù)據(jù)來源

    本研究采用覆蓋廣西區(qū)(圖1(a))的MODISMOD15A2H 葉面積指數(shù)(LAI)數(shù)據(jù)產(chǎn)品,來源于美國航天航空局(NASA,https://ladsweb.modaps.eosdis.nasa.gov/),空間分辨率500 m,時間分辨率8 d,時段為2014—2018 年1—12 月,共230 期影像數(shù)據(jù)。文中使用的甘蔗種植空間分布數(shù)據(jù)來源于本課題組采用隨機(jī)森林分類方法提取2014—2018 年的LANDSAT8 OIL(30 m 空間分辨率)遙感影像獲得(圖1(b)—(f)),該數(shù)據(jù)總體識別精度達(dá)到90%以上,Kappa 系數(shù)均超過0.8。此外,采用廣西6 個農(nóng)業(yè)氣象站統(tǒng)計的2014—2018 年甘蔗物候期數(shù)據(jù)作為精度驗證(表1)。

    表1 廣西甘蔗物候期統(tǒng)計結(jié)果(d)Table 1 Statistical results of sugarcane phenology in Guangxi (d)

    1.2 數(shù)據(jù)處理與函數(shù)擬合方法

    1.2.1 MODIS-LAI 數(shù)據(jù)預(yù)處理 MOD15A2H 數(shù)據(jù)的預(yù)處理主要采用NASA 官網(wǎng)軟件MRT(MODIS Reprojection Tool)工具;首先對該數(shù)據(jù)進(jìn)行重投影、拼接及鑲嵌、數(shù)據(jù)重采樣等處理,然后通過ENVIIDL 軟件實現(xiàn)廣西研究區(qū)及甘蔗種植區(qū)數(shù)據(jù)的裁剪。

    1.2.2 函數(shù)擬合方法 MODIS-LAI 數(shù)據(jù)產(chǎn)品雖然經(jīng)過8 d最大值合成,減弱了受大氣中云和氣溶膠的影響,但是由于衛(wèi)星傳感器角度、太陽高度角、水蒸氣等因素的變化產(chǎn)生影像噪聲[1],致使LAI 數(shù)據(jù)產(chǎn)品的時序特征存在顯著的不規(guī)則波動,難以真實地反映作物的時序特征。為此,文中對MODIS-LAI 時序數(shù)據(jù)進(jìn)行必要的濾波處理,通過平滑降噪去空值處理的方法來削弱這種不規(guī)則波動的影響;具體選擇了S-G 濾波函數(shù)法(Savitzky-Golay)、非對稱高斯擬合法(Asymmetric Gauss)及雙Logistic 函數(shù)擬合法(Double Logistic)等。

    1)S-G 濾波函數(shù)法。該方法是由Savitzky[20]等提出的一種移動窗口式加權(quán)平均算法,其加權(quán)系數(shù)主要通過最小二乘法擬合得出,計算公式如下:

    式中:Y指像元對應(yīng)的LAI 原始值;Y*指對應(yīng)像元通過濾波處理后的LAI 值;Ci為第i期LAI 值濾波過程的系數(shù);j指原始LAI 數(shù)組的系數(shù)。

    2)非對稱高斯擬合法。Jonsson 和Eklundh[21]提出了一種通過局部時序曲線擬合拓展到整體擬合的平滑函數(shù)處理法,即非對稱高斯擬合法,主要是通過不同的分段高斯函數(shù)模擬植被或作物的生長情況,最后采用平滑連接的方法實現(xiàn)整個時序曲線擬合,完成時間序列重建。局部擬合函數(shù)如下式(2)和(3):

    式中,g(t,a1,…a5)為高斯函數(shù);a1為決定最大值和最小值所在整條時序曲線的位置;c1和c2為控制整條曲線的基準(zhǔn)及振幅;a4、a5和a2、a3分別為控制曲線左、右部分的寬度及斜率。

    式中:[tL,tR]表示時序曲線中待擬合部分的取值區(qū)間范圍;fL(t)、fC(t)和fR(t)分別為該區(qū)間[tL,tR]內(nèi)左邊最小值、中間最大值及右邊最小值所對應(yīng)的局部擬合函數(shù);α(t)和β(t)分別為介于0 到1 之間的剪切系數(shù)。

    3)雙Logistic 函數(shù)擬合法。雙Logistic 函數(shù)擬合法是由Beck 等[22]提出,其擬合過程與上述非對稱高斯擬合類似,都需要進(jìn)行局部擬合再到整體擬合,其相比基于傅里葉變換的濾波函數(shù)擬合法,對植被或作物生長過程的擬合更為準(zhǔn)確,對于高緯度地區(qū)的作物生長模擬也普遍適用,其擬合結(jié)果與非對稱高斯擬合相似,不同在于所采用的局部擬合函數(shù)為雙邏輯形式,且公式中缺少一個參數(shù),僅有在大量統(tǒng)計對比中非對稱高斯擬合對比該方法才會產(chǎn)生較小優(yōu)勢。其局部擬合公式如下:

    式中:g(t,a1,…,a5)為高斯函數(shù),a1為決定最大值和最小值所在整條時序曲線的位置,a2、a3和a4、a5分別為控制整條曲線左右半部分的寬度及斜率,最后通過整體擬合函數(shù)實現(xiàn)MODIS-LAI 時間序列重構(gòu)。

    1.3 物候期提取與精度檢驗方法

    1.3.1 物候期提取 TIMESAT 是一款由Jonsson 和Eklundh[23]基于Matlab 程序開發(fā)的專門用來處理時間序列衛(wèi)星遙感數(shù)據(jù)的軟件,能很好地實現(xiàn)植被或作物的物候期動態(tài)特性分析,并獲取短期甚至長期的植被或作物生長發(fā)育信息。該軟件中能夠?qū)崿F(xiàn)S-G濾波函數(shù)法(Savitzky-Golay)、非對稱高斯擬合法(Asymmetric Gauss)及雙Logistic函數(shù)擬合法(Double Logistic)等多種時間序列降噪重構(gòu)的平滑處理,本文主要采用動態(tài)閾值方法[23]并結(jié)合重構(gòu)的作物生長曲線來實現(xiàn)廣西甘蔗播種-萌芽期、莖伸期及成熟期的提取。其中本研究主要將20%閾值位置處確定為甘蔗的播種-萌芽期,也就是LAI 曲線上升階段最大值的20%位置處所對應(yīng)的時間即為甘蔗的播種-萌芽期開始時間;下降段最大值的20%處所對應(yīng)的時間即為甘蔗成熟期開始的時間[23];根據(jù)FAO作物生育期的規(guī)劃方法[18],其中作物生長中期對應(yīng)甘蔗生育的莖伸期,即甘蔗LAI 曲線達(dá)到峰值所對應(yīng)的時間即為莖伸期開始時間。

    1.3.2 精度檢驗 本文主要通過圖1(a)中6 個主要農(nóng)業(yè)氣象站點的2014—2018 年甘蔗物候期地面觀測數(shù)據(jù)對廣西甘蔗生育期遙感識別結(jié)果進(jìn)行精度驗證,分別采用最大誤差、最小誤差、平均誤差及均方根誤差等指標(biāo)驗證甘蔗物候期識別的準(zhǔn)確度。均方根誤差計算公式如下:

    式中:di為廣西甘蔗物候期提取的反演日期,為對應(yīng)的廣西甘蔗物候期的觀測日期,d;n為樣本總數(shù)。

    2 結(jié)果與分析

    2.1 擬合結(jié)果比較

    選取廣西崇左市扶綏縣甘蔗種植區(qū)域內(nèi)的某甘蔗像元,并查看該像元2014—2018 年LAI 時序重構(gòu)曲線結(jié)果;由圖2 可知,主要受到衛(wèi)星傳感器角度及水蒸氣等不確定性因素影響,不論是哪年的LAI 原始時序曲線均呈現(xiàn)出顯著的鋸齒狀波動,不能準(zhǔn)確地反映整個甘蔗物候期的長勢變化特征。因此,文中通過平滑函數(shù)處理方法可以很好地保持原始時序曲線的時間變化形態(tài),并校正極低的離異值,使整個時序曲線擬合更接近甘蔗生長發(fā)育的真實情況。從圖2 可以看出,S-G 濾波擬合效果較低,雖然消除了曲線不規(guī)則波動的影響,但局部曲線仍存在不規(guī)則波動變化,主要由于S-G 濾波處理的局限性導(dǎo)致的;根據(jù)楊琳等[1]、邊金虎等[24]研究可知,S-G 濾波窗口大小及函數(shù)階數(shù)的選取需要根據(jù)不同的條件進(jìn)行主觀人為地設(shè)定,如果選取的濾波窗口過大則容易遺漏掉較多的作物生長發(fā)育信息;反之偏小,則導(dǎo)致數(shù)據(jù)冗余,不能及時、準(zhǔn)確地獲取到整個物候期的長期長勢變化趨勢;S-G 函數(shù)的階數(shù)使用過高,則容易出現(xiàn)曲線過擬合,且局部曲線數(shù)值就會產(chǎn)生不規(guī)則波動的情況。而非高斯對稱及雙Logistic 函數(shù)擬合的效果較好,消除了S-G 擬合存在的鋸齒狀波動情況,較好地模擬了整個甘蔗生育期的生長發(fā)育,并保持了S-G 濾波及原始數(shù)據(jù)長勢季節(jié)性變化的趨勢。

    由圖2(a)—(e)結(jié)合廣西甘蔗的實際生長發(fā)育概況可知,廣西甘蔗的播種-萌芽期一般在2—6 月份,因此前期甘蔗的LAI 值增長緩慢;6—8 月份廣西甘蔗進(jìn)入分蘗期-莖伸期后,其生長速度加快,生長力旺盛,LAI 指數(shù)變化加快,曲線呈現(xiàn)出陡峭上升趨勢,直到8 月中旬左右達(dá)到最高值;9月下旬后進(jìn)入甘蔗莖伸期末,生長緩慢,曲線呈緩慢下降趨勢;11 月份前后進(jìn)入成熟期,曲線下降速度加快,該甘蔗物候曲線變化趨勢與丁美花[25-26]等在廣西甘蔗長勢監(jiān)測及長勢信息提取研究中的結(jié)論相一致,表明采用平滑函數(shù)降噪取空值的擬合方法可以較好地模擬廣西甘蔗長勢變化情況。

    綜上可知,S-G 濾波窗口及函數(shù)階數(shù)是S-G 函數(shù)準(zhǔn)確擬合及獲取作物長勢變化信息的關(guān)鍵因素;而非高斯對稱及雙Logistic 函數(shù)擬合的效果較好地彌補了S-G濾波處理中局部曲線值動蕩變化的缺點,并保持了原始數(shù)據(jù)中作物長勢季節(jié)性變化的趨勢,能更準(zhǔn)確地模擬廣西甘蔗生長過程。

    2.2 廣西甘蔗關(guān)鍵物候期提取結(jié)果精度驗證

    文中采用6 個主要農(nóng)業(yè)氣象站點的2014—2018年甘蔗生育期地面觀測數(shù)據(jù)對上述三種方法提取的廣西甘蔗物候期結(jié)果進(jìn)行了對比和精度驗證(圖3和表2)。

    由圖3 可知,2014—2018 年的甘蔗播種-萌芽期、莖伸期及成熟期的提取誤差普遍都在15 d 以內(nèi),只有少數(shù)是超過15 d;其中2014—2018 年,基于S-G濾波擬合提取的甘蔗物候期誤差超過15 d 的樣本數(shù)為播種-萌芽期3 個(10%),莖伸期7 個(23.33%),成熟期9 個(30%);2014—2018 年基于A-G 函數(shù)擬合提取的甘蔗物候期誤差超過15 d 的樣本數(shù)為播種-萌芽期0 個,莖伸期2 個(6.67%),成熟期6個(20%);2014—2018 年基于D-L 函數(shù)擬合提取的甘蔗物候期誤差超過15 d 的樣本數(shù)為播種-萌芽期2 個(6.67%),莖伸期3 個(10%),成熟期8 個(26.67%)。

    由表2 可知,基于S-G 濾波擬合提取的播種-萌芽期、莖伸期、成熟期開始時間的均方根誤差分別為10.43 d、11.67 d、12.46 d,而平均誤差分別為9.16 d、10.17 d、11.28 d;基于A-G 函數(shù)擬合提取的三個關(guān)鍵物候期開始時間的均方根誤差分別為7.63 d、8.43 d、10.48 d,其中播種-萌芽期及成熟期精度較高,平均誤差分別為6.69 d、6.91 d、9.13 d;基于D-L 函數(shù)擬合提取結(jié)果的均方根誤差分別為9.02 d、10.01 d、11.76 d,平均誤差為7.75 d、8.42 d、10.28 d。

    綜合圖3 和表2 結(jié)果表明,S-G、A-G 和D-L對廣西甘蔗的三個關(guān)鍵物候期提取結(jié)果的均方根誤差及平均誤差均在15 d 內(nèi),其中A-G、D-L 方法的提取結(jié)果優(yōu)于S-G 方法結(jié)果,而A-G 結(jié)果則較D-L結(jié)果優(yōu),也即A-G 函數(shù)擬合方法對廣西甘蔗的物候期識別效果最好,在廣西甘蔗物候監(jiān)測識別方面具有良好的應(yīng)用價值。

    2.3 廣西甘蔗關(guān)鍵物候期空間分布格局

    由3.2 節(jié)的甘蔗物候期提取精度分析可知,非對稱高斯函數(shù)(A-G)擬合法對廣西甘蔗物候期的提取效果最優(yōu),其中圖5 為基于該方法對MODISLAI 時序曲線重構(gòu),并采用動態(tài)閾值法提取的2014—2018 年廣西甘蔗播種-萌芽期、莖伸期和成熟期開始時間的空間差異分布圖。

    由圖4 可知:2014 年廣西甘蔗播種-萌芽期的開始時間呈“桂東向桂西推遲”的空間特征,桂東開始時間主要集中在第121~153 d(即4 月1 日到5 月2 日),桂西的開始時間則推遲至第153 d 后;2015—2018 年廣西甘蔗播種-萌芽期開始時間的空間格局相似,整體趨于5 月2 日后。

    2014、2016 年廣西甘蔗莖伸期開始時間呈“桂西南向桂東北提前”的空間分布格局,由第281 d后向217 d前過渡(即10月8日后提前至8月5日前);2015 年甘蔗莖伸期開始時間的空間格局則整體推遲至10 月8 日后;2017 年甘蔗莖伸期開始時間呈“桂東向桂西提前”的空間分布特征,即時間提前至8月5 日前;2018 年莖伸期開始時間則整體集中在8月5 日到10 月8 日。

    2014 年廣西甘蔗成熟期開始時間的空間特征與莖伸期相似,即由2014 年12 月11 日—2015 年1月12 日提前至2014 年11 月9 日—12 月10 日(即第346~377 d 提前至第313~345 d);2015 年成熟期開始時間則由11 月9 日前(桂東北)推遲至12 月10 日后(桂西南);2016、2018 年成熟期開始時間則呈“桂西向桂東提前”的空間變化特征,由2016年12 月10 日后提前至11 月9 日前。

    綜上可知,通過遙感技術(shù)提取的廣西甘蔗物候信息能夠很好地反映甘蔗生長發(fā)育的時空變化過程,對于廣西甘蔗物候監(jiān)測的年際、年內(nèi)變化分析具有較好的應(yīng)用價值。

    2.4 物候提取結(jié)果的影響因素分析

    表2 2014—2018 年廣西甘蔗物候期提取精度評價(d)Table 2 Accuracy evaluation of sugarcane phenology inversion in Guangxi from 2014 to 2018 (d)

    綜上可知,基于三種方法重構(gòu)MODIS-LAI 時序數(shù)據(jù),并通過動態(tài)閾值法提取的廣西甘蔗物候期存在顯著差異,這與三種序列重構(gòu)法本身的特點密切相關(guān)。首先,A-G 和D-L 擬合方法均屬于函數(shù)擬合法,通過局部擬合擴(kuò)展至整體,受到擬合函數(shù)自身形式的影響較大,可獲得平滑的重構(gòu)曲線;而S-G濾波屬于時域濾波法,屬于局部處理方法,時序曲線的平滑程度與濾波窗口及函數(shù)階數(shù)的選取有關(guān),人為主觀性較大,受原始數(shù)據(jù)質(zhì)量的影響大。其次,A-G 函數(shù)和D-L 函數(shù)擬合的基本原理相同,僅是局部擬合函數(shù)不同,導(dǎo)致兩種方法無論是時間還是空間上得到的甘蔗物候期模擬結(jié)果非常相似;S-G 濾波法的原理和前者完全不同,因此,造成其與A-G和D-L 函數(shù)擬合法提取甘蔗物候期的時空信息存在顯著差異,而該曲線上升及下降階段的左右偏移也是造成物候期出現(xiàn)提前或推遲的主要原因。

    雖然通過AG 函數(shù)提取的廣西甘蔗物候期效果最優(yōu),但其估測結(jié)果仍與地面觀測值存在誤差,而導(dǎo)致該誤差的主要原因有:

    1)數(shù)據(jù)方面,MODIS-LAI 數(shù)據(jù)受空間分辨率限制(空間分辨率500 m),其值僅代表像元平均值,而地面觀測數(shù)據(jù)來自農(nóng)業(yè)站點的物候觀測值,二者空間坐標(biāo)匹配存在明顯差異,影響了遙感提取的物候期與觀測值空間分辨率的一致性;但農(nóng)業(yè)氣象站點并不是在甘蔗田間的單點設(shè)置,而是在當(dāng)?shù)馗收岱N植區(qū)范圍內(nèi)布置多個觀測站點,通過多站點的共同協(xié)調(diào)作用實現(xiàn)當(dāng)?shù)馗收嵛锖虮O(jiān)測,所以站點布置的觀測范圍與MODIS-LAI 的空間分辨率接近,采用該遙感數(shù)據(jù)進(jìn)行廣西甘蔗物候信息監(jiān)測具有一定的合理性。

    2)廣西屬于亞熱帶季風(fēng)區(qū),降雨充沛,在甘蔗等作物生育季期間云雨較多,對于受云霧影響的區(qū)域,甘蔗物候期內(nèi)LAI 值普遍偏低,時序曲線波動變化顯著,不能反映真實的甘蔗生長發(fā)育狀況;而時序曲線重構(gòu)法具有較好的魯棒性,可更好地提高影像數(shù)據(jù)質(zhì)量,并獲得連續(xù)完整的甘蔗區(qū)LAI 時序曲線,使其更接近甘蔗的物候生長過程,具有一定的合理性;但其LAI 模擬值仍與甘蔗的實際值存在一定誤差,因此,導(dǎo)致遙感提取的甘蔗物候期與真實情況對比存在高估或低估現(xiàn)象。

    3)廣西山地丘陵偏多,地勢復(fù)雜,農(nóng)田地塊破碎,插花種植現(xiàn)象嚴(yán)重,導(dǎo)致甘蔗種植區(qū)解譯結(jié)果偏少或偏多,且影像獲取的蔗區(qū)位置相對偏移,與真實的甘蔗空間種植結(jié)構(gòu)分布存在一定偏差,進(jìn)而影響甘蔗的物候提取結(jié)果;且山地丘陵區(qū)也屬于云雨及大霧天氣易發(fā)區(qū),嚴(yán)重影響遙感數(shù)據(jù)質(zhì)量,與2)中影響效應(yīng)相同。

    隨著遙感技術(shù)的發(fā)展,GF-1、GF-6、Sentinel-2、ZY-3 等中高分辨率影像的日益豐富,有效積溫、太陽輻射等有關(guān)作物生長因素的融入,有助于實現(xiàn)更加精準(zhǔn)的作物物候信息提取。因此,采用遙感提取的作物物候信息與地面觀測相結(jié)合是實現(xiàn)作物物候信息監(jiān)測的有效途徑,也是未來農(nóng)業(yè)資源遙感監(jiān)測的研究熱點之一。

    3 結(jié)論

    本文基于MODIS-LAI 數(shù)據(jù),通過三種時序重構(gòu)方法實現(xiàn)了廣西甘蔗關(guān)鍵物候期提取,得到以下主要結(jié)論:

    1)針對MODIS-LAI 數(shù)據(jù)時序曲線存在鋸齒狀波動而難以準(zhǔn)確反映廣西甘蔗生育期長勢變化的問題,S-G 濾波函數(shù)法、非對稱高斯擬合法和雙Logistic 函數(shù)擬合法等均可有效消除該序列的不穩(wěn)定波動及其奇異值,但后兩種方法的擬合結(jié)果較好地彌補了S-G 濾波處理中局部曲線值動蕩變化的缺點,能更準(zhǔn)確地獲取廣西甘蔗生長過程信息。

    2)在廣西甘蔗物候期提取精度方面,三種方法提取結(jié)果的均方根誤差及平均誤差普遍在15 d內(nèi);其中,A-G 方法的精度最高,其次為D-L 方法,S-G 方法的提取結(jié)果精度相對最差。

    3)由于S-G 函數(shù)法的濾波窗口及函數(shù)階數(shù)選取的主觀性和A-G、D-L 擬合法自身函數(shù)形式的限制,導(dǎo)致三者在時空上獲取的廣西甘蔗物候期模擬結(jié)果存在明顯的差異性;地面觀測數(shù)據(jù)與遙感數(shù)據(jù)的空間一致性,云雨及大霧天氣導(dǎo)致遙感數(shù)據(jù)質(zhì)量的下降及廣西復(fù)雜地勢條件對甘蔗種植分布解譯結(jié)果的干擾等均是造成廣西甘蔗物候期估測結(jié)果與地面觀測值存在誤差的主要原因。

    猜你喜歡
    物候時序甘蔗
    時序坐標(biāo)
    “蔗”里時光
    花式賣甘蔗
    海南橡膠林生態(tài)系統(tǒng)凈碳交換物候特征
    清明甘蔗“毒過蛇”
    基于Sentinel-2時序NDVI的麥冬識別研究
    愛咬甘蔗的百歲爺爺
    特別健康(2018年3期)2018-07-04 00:40:08
    ‘灰棗’及其芽變品系的物候和生育特性研究
    一種毫米波放大器時序直流電源的設(shè)計
    電子制作(2016年15期)2017-01-15 13:39:08
    5種忍冬科植物物候期觀察和比較
    麻豆国产97在线/欧美| 美女高潮的动态| 日日摸夜夜添夜夜爱| 中文字幕熟女人妻在线| 最近中文字幕高清免费大全6| 亚洲熟妇中文字幕五十中出| 男的添女的下面高潮视频| 国内精品久久久久精免费| 美女内射精品一级片tv| 91av网一区二区| 亚洲性久久影院| 狂野欧美白嫩少妇大欣赏| 亚洲综合色惰| 国产极品精品免费视频能看的| 国产精品久久电影中文字幕| 欧美zozozo另类| avwww免费| 九草在线视频观看| a级毛色黄片| 18禁在线播放成人免费| 亚洲精品成人久久久久久| 在线观看一区二区三区| 色哟哟哟哟哟哟| 欧美+日韩+精品| 国产 一区 欧美 日韩| 青春草国产在线视频 | 哪里可以看免费的av片| 国产真实伦视频高清在线观看| 免费观看在线日韩| 国产伦精品一区二区三区四那| 欧美极品一区二区三区四区| 日韩亚洲欧美综合| 国国产精品蜜臀av免费| 国产一区亚洲一区在线观看| 日本五十路高清| 99视频精品全部免费 在线| 91狼人影院| 高清在线视频一区二区三区 | 亚洲av免费在线观看| 日韩一本色道免费dvd| 国产伦一二天堂av在线观看| 久久精品久久久久久噜噜老黄 | 国产一级毛片七仙女欲春2| 国产精品人妻久久久影院| 最好的美女福利视频网| av.在线天堂| 国产探花极品一区二区| 久久这里有精品视频免费| 日韩欧美国产在线观看| 免费看av在线观看网站| 国产老妇伦熟女老妇高清| 免费观看人在逋| 51国产日韩欧美| 久久精品影院6| 男人舔女人下体高潮全视频| 人妻夜夜爽99麻豆av| 久久久久久久久久成人| 久久人妻av系列| 日本黄色片子视频| 午夜精品在线福利| 国产高清有码在线观看视频| 深夜精品福利| 黄片wwwwww| 亚洲精品亚洲一区二区| 插逼视频在线观看| 亚洲国产日韩欧美精品在线观看| 久久人妻av系列| 成年av动漫网址| 免费在线观看成人毛片| 午夜精品在线福利| 日韩av在线大香蕉| 亚洲国产欧美人成| 一区福利在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲久久久久久中文字幕| 久久久久性生活片| 男插女下体视频免费在线播放| 日本av手机在线免费观看| 久久久久久大精品| 亚州av有码| 久久精品夜夜夜夜夜久久蜜豆| 亚洲乱码一区二区免费版| 12—13女人毛片做爰片一| 亚洲国产精品合色在线| 亚洲国产日韩欧美精品在线观看| 亚洲在线观看片| 久久精品国产鲁丝片午夜精品| 极品教师在线视频| 在线免费观看不下载黄p国产| 给我免费播放毛片高清在线观看| 99久久无色码亚洲精品果冻| 深爱激情五月婷婷| 久久久久久久午夜电影| 免费电影在线观看免费观看| 男女下面进入的视频免费午夜| 亚洲精品自拍成人| av在线蜜桃| 熟女电影av网| 久久久久久久久久久免费av| 精品少妇黑人巨大在线播放 | 成人欧美大片| 国产精品麻豆人妻色哟哟久久 | 久久久久久久久久黄片| 亚洲欧美成人综合另类久久久 | 中文字幕av成人在线电影| 精品一区二区三区人妻视频| 变态另类成人亚洲欧美熟女| 成人高潮视频无遮挡免费网站| 国产精品不卡视频一区二区| 青春草视频在线免费观看| 最近最新中文字幕大全电影3| 亚洲精品乱码久久久v下载方式| 成人鲁丝片一二三区免费| 国产三级中文精品| 国产成人精品久久久久久| 国产老妇女一区| 久久久久久久久久久免费av| 一个人看视频在线观看www免费| 欧美日本视频| 国产探花在线观看一区二区| 精品久久久久久久久久免费视频| 国内久久婷婷六月综合欲色啪| 欧美高清性xxxxhd video| 特大巨黑吊av在线直播| 亚洲一区高清亚洲精品| 中文字幕精品亚洲无线码一区| 成人av在线播放网站| 国产真实乱freesex| 国产淫片久久久久久久久| 1024手机看黄色片| 好男人在线观看高清免费视频| 国产成人a∨麻豆精品| 亚洲av一区综合| av在线天堂中文字幕| 狂野欧美激情性xxxx在线观看| 国产成人福利小说| 国产av不卡久久| 国产精品,欧美在线| 久久99热这里只有精品18| 国产精品福利在线免费观看| 日本免费一区二区三区高清不卡| 美女脱内裤让男人舔精品视频 | 欧美bdsm另类| 国产av不卡久久| 69人妻影院| 久久久国产成人免费| 高清毛片免费看| 一级毛片aaaaaa免费看小| 老熟妇乱子伦视频在线观看| 又粗又爽又猛毛片免费看| 欧美+亚洲+日韩+国产| 亚洲精品自拍成人| 国产极品天堂在线| 中文欧美无线码| 五月伊人婷婷丁香| 日韩欧美精品v在线| av免费观看日本| 久久精品国产自在天天线| 亚洲欧洲国产日韩| 久久久成人免费电影| 高清日韩中文字幕在线| 亚洲一级一片aⅴ在线观看| 久久中文看片网| 亚洲熟妇中文字幕五十中出| 中文精品一卡2卡3卡4更新| 哪里可以看免费的av片| 在线观看美女被高潮喷水网站| 日韩强制内射视频| 日韩三级伦理在线观看| avwww免费| www日本黄色视频网| 少妇人妻一区二区三区视频| 午夜福利在线在线| 亚洲成人久久性| 欧美日韩综合久久久久久| 亚洲欧美日韩东京热| 亚洲精华国产精华液的使用体验 | 91aial.com中文字幕在线观看| 91久久精品电影网| 欧美又色又爽又黄视频| 蜜臀久久99精品久久宅男| 免费观看在线日韩| 黄色日韩在线| 人人妻人人看人人澡| 级片在线观看| 国内久久婷婷六月综合欲色啪| 麻豆av噜噜一区二区三区| 国产精品美女特级片免费视频播放器| 99久久精品国产国产毛片| av卡一久久| 久久6这里有精品| 最近中文字幕高清免费大全6| 18禁在线播放成人免费| 我要看日韩黄色一级片| 在线观看美女被高潮喷水网站| 综合色丁香网| 一区二区三区免费毛片| 1024手机看黄色片| 日本撒尿小便嘘嘘汇集6| 国产精品蜜桃在线观看 | 久久鲁丝午夜福利片| 夫妻性生交免费视频一级片| 国产免费一级a男人的天堂| 人人妻人人澡欧美一区二区| 午夜a级毛片| 亚洲av成人av| 99九九线精品视频在线观看视频| 亚洲av熟女| 亚洲欧美清纯卡通| 国产黄色视频一区二区在线观看 | 啦啦啦啦在线视频资源| 国产精品久久电影中文字幕| 久久精品国产鲁丝片午夜精品| 免费人成在线观看视频色| 一个人看视频在线观看www免费| 成人亚洲欧美一区二区av| 亚洲中文字幕一区二区三区有码在线看| 日日啪夜夜撸| 婷婷六月久久综合丁香| 亚洲精品粉嫩美女一区| 亚洲av一区综合| 一本久久中文字幕| a级毛色黄片| 日韩欧美国产在线观看| 亚洲人成网站在线观看播放| 久久久久国产网址| 精品熟女少妇av免费看| 国产三级在线视频| 我的女老师完整版在线观看| or卡值多少钱| 精品久久久久久久末码| 能在线免费观看的黄片| 欧美另类亚洲清纯唯美| 中文字幕久久专区| 成人高潮视频无遮挡免费网站| 亚洲一区二区三区色噜噜| 国产亚洲5aaaaa淫片| 国产私拍福利视频在线观看| 中文字幕熟女人妻在线| 六月丁香七月| 久久人人爽人人片av| 高清毛片免费看| 日韩一区二区视频免费看| 久久鲁丝午夜福利片| 日本黄色视频三级网站网址| av天堂在线播放| 日韩欧美一区二区三区在线观看| 久久久久久国产a免费观看| 国产黄片视频在线免费观看| 免费观看的影片在线观看| 99热全是精品| 亚洲国产精品国产精品| 国产精品野战在线观看| 黄色一级大片看看| 午夜福利在线在线| 色综合亚洲欧美另类图片| 中国国产av一级| 美女cb高潮喷水在线观看| 嘟嘟电影网在线观看| 一级毛片aaaaaa免费看小| 大香蕉久久网| 亚洲精品亚洲一区二区| 91麻豆精品激情在线观看国产| 久久综合国产亚洲精品| 一本一本综合久久| 中文字幕免费在线视频6| 在线国产一区二区在线| 日本黄色视频三级网站网址| 亚洲精品日韩av片在线观看| 国产精品三级大全| 欧美xxxx黑人xx丫x性爽| 国产69精品久久久久777片| 亚洲三级黄色毛片| 亚洲色图av天堂| 欧美+亚洲+日韩+国产| 国产精品,欧美在线| av视频在线观看入口| 午夜激情福利司机影院| av黄色大香蕉| 国产精品不卡视频一区二区| 国产极品精品免费视频能看的| 国产高清激情床上av| 亚洲精品国产成人久久av| 美女cb高潮喷水在线观看| 成熟少妇高潮喷水视频| 国产在线男女| 日本-黄色视频高清免费观看| 26uuu在线亚洲综合色| a级毛片a级免费在线| 久久九九热精品免费| 内地一区二区视频在线| 中文资源天堂在线| 边亲边吃奶的免费视频| 日韩制服骚丝袜av| 精品免费久久久久久久清纯| 成人性生交大片免费视频hd| 精品人妻一区二区三区麻豆| 少妇熟女aⅴ在线视频| 中文字幕av成人在线电影| 亚洲精品色激情综合| 天天躁日日操中文字幕| 国产单亲对白刺激| 久久久久久久久中文| 三级国产精品欧美在线观看| 色综合亚洲欧美另类图片| 狠狠狠狠99中文字幕| 在线免费观看不下载黄p国产| 国产色爽女视频免费观看| www.av在线官网国产| 亚洲av中文av极速乱| 日韩三级伦理在线观看| 好男人视频免费观看在线| 国产精品伦人一区二区| 成人一区二区视频在线观看| 小说图片视频综合网站| 亚洲精品成人久久久久久| 狂野欧美激情性xxxx在线观看| 日本欧美国产在线视频| 久久精品国产亚洲av天美| 国产免费男女视频| 一进一出抽搐gif免费好疼| 亚洲av男天堂| 色播亚洲综合网| 久久精品国产99精品国产亚洲性色| a级一级毛片免费在线观看| 欧美日韩一区二区视频在线观看视频在线 | 少妇人妻精品综合一区二区 | 免费av不卡在线播放| 男人舔女人下体高潮全视频| 日韩成人av中文字幕在线观看| 麻豆精品久久久久久蜜桃| 国产成人午夜福利电影在线观看| 欧美一区二区精品小视频在线| 亚洲18禁久久av| 午夜免费激情av| 天堂影院成人在线观看| 免费av不卡在线播放| 国产高清激情床上av| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 欧美人与善性xxx| 久久99蜜桃精品久久| 青春草亚洲视频在线观看| 99热全是精品| 在线播放无遮挡| 99国产极品粉嫩在线观看| 亚洲无线在线观看| 少妇人妻精品综合一区二区 | 亚洲成人久久性| 国产国拍精品亚洲av在线观看| 亚洲四区av| 国产成人a区在线观看| 国产免费一级a男人的天堂| 国产亚洲91精品色在线| 欧美色视频一区免费| 99热精品在线国产| 成人毛片a级毛片在线播放| 国产男人的电影天堂91| 毛片女人毛片| 午夜精品一区二区三区免费看| 精品熟女少妇av免费看| 深夜精品福利| 久久久国产成人免费| 中文字幕av在线有码专区| 色综合亚洲欧美另类图片| 欧美丝袜亚洲另类| 亚洲国产日韩欧美精品在线观看| 2022亚洲国产成人精品| 一级毛片电影观看 | 国产精品,欧美在线| 美女内射精品一级片tv| 中文字幕制服av| 大香蕉久久网| 欧美日韩精品成人综合77777| 乱人视频在线观看| 能在线免费看毛片的网站| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久久成人| 99国产极品粉嫩在线观看| 成人特级av手机在线观看| 日韩亚洲欧美综合| 国产精品人妻久久久久久| 久久99蜜桃精品久久| 嫩草影院新地址| 亚洲国产精品久久男人天堂| 伊人久久精品亚洲午夜| 2021天堂中文幕一二区在线观| 国产一区二区在线av高清观看| 国产亚洲av嫩草精品影院| 又爽又黄a免费视频| 三级经典国产精品| 亚洲av.av天堂| 97热精品久久久久久| 国产乱人偷精品视频| 少妇高潮的动态图| eeuss影院久久| 91在线精品国自产拍蜜月| 久久草成人影院| 少妇熟女aⅴ在线视频| 精华霜和精华液先用哪个| 日本免费a在线| 特大巨黑吊av在线直播| 婷婷精品国产亚洲av| 爱豆传媒免费全集在线观看| 最近2019中文字幕mv第一页| 特级一级黄色大片| www.色视频.com| 人妻少妇偷人精品九色| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 久久久久网色| 亚洲性久久影院| 热99在线观看视频| 99久国产av精品国产电影| 午夜免费男女啪啪视频观看| 久久国产乱子免费精品| 能在线免费看毛片的网站| 看片在线看免费视频| eeuss影院久久| 精品欧美国产一区二区三| 亚洲无线观看免费| 色综合亚洲欧美另类图片| 亚洲成人久久爱视频| 九草在线视频观看| 国产高清视频在线观看网站| 麻豆av噜噜一区二区三区| 日本成人三级电影网站| 欧美性猛交黑人性爽| 精品99又大又爽又粗少妇毛片| 精品熟女少妇av免费看| 久久精品综合一区二区三区| 午夜福利成人在线免费观看| 久久精品91蜜桃| 国产在线精品亚洲第一网站| 网址你懂的国产日韩在线| 少妇的逼水好多| 大香蕉久久网| 亚洲欧美成人精品一区二区| 国产女主播在线喷水免费视频网站 | 99热这里只有精品一区| 一夜夜www| 99热这里只有是精品在线观看| 啦啦啦观看免费观看视频高清| 欧美+亚洲+日韩+国产| 中文精品一卡2卡3卡4更新| 精品欧美国产一区二区三| 国产av麻豆久久久久久久| 国产免费一级a男人的天堂| 国产精品精品国产色婷婷| 国产视频内射| 搞女人的毛片| 免费不卡的大黄色大毛片视频在线观看 | av在线播放精品| 国产又黄又爽又无遮挡在线| av福利片在线观看| АⅤ资源中文在线天堂| 99热6这里只有精品| 日日摸夜夜添夜夜添av毛片| 在线观看午夜福利视频| 亚洲国产精品国产精品| 国产成人精品久久久久久| 观看美女的网站| 国产精品蜜桃在线观看 | eeuss影院久久| 色播亚洲综合网| 成人综合一区亚洲| 亚洲自拍偷在线| 少妇猛男粗大的猛烈进出视频 | 极品教师在线视频| 成年版毛片免费区| 久久久久久久午夜电影| 免费在线观看成人毛片| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 最近的中文字幕免费完整| 婷婷精品国产亚洲av| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 亚洲国产精品合色在线| 国产精品久久久久久精品电影| 人妻系列 视频| 欧美不卡视频在线免费观看| 99在线人妻在线中文字幕| 日本黄色片子视频| 国产老妇女一区| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 亚洲精品自拍成人| 日韩av在线大香蕉| 看黄色毛片网站| 男人舔女人下体高潮全视频| 久久九九热精品免费| 麻豆成人av视频| 国产精品免费一区二区三区在线| 久久草成人影院| 欧美日韩综合久久久久久| 听说在线观看完整版免费高清| 特大巨黑吊av在线直播| 欧美区成人在线视频| 一级二级三级毛片免费看| 级片在线观看| 精品人妻视频免费看| 一边摸一边抽搐一进一小说| 天天一区二区日本电影三级| 美女黄网站色视频| 欧美成人免费av一区二区三区| 精品一区二区三区人妻视频| 一级二级三级毛片免费看| 青春草视频在线免费观看| 日韩欧美精品v在线| 三级毛片av免费| 欧美激情久久久久久爽电影| 国国产精品蜜臀av免费| 亚洲欧美精品综合久久99| 国产免费一级a男人的天堂| av卡一久久| 久久这里有精品视频免费| 高清日韩中文字幕在线| 一级黄片播放器| 国产精品综合久久久久久久免费| 一夜夜www| 亚洲国产日韩欧美精品在线观看| 久久久久久久久大av| 久久精品国产亚洲网站| 亚洲一区二区三区色噜噜| 一个人免费在线观看电影| 中文字幕av在线有码专区| 九九久久精品国产亚洲av麻豆| 极品教师在线视频| 99久久久亚洲精品蜜臀av| 日本一二三区视频观看| 欧美日韩国产亚洲二区| 国产黄a三级三级三级人| 精品日产1卡2卡| 免费看a级黄色片| 一本一本综合久久| 成年版毛片免费区| 免费看美女性在线毛片视频| 老师上课跳d突然被开到最大视频| 日韩精品有码人妻一区| 国产精品美女特级片免费视频播放器| 久久精品91蜜桃| а√天堂www在线а√下载| 伦精品一区二区三区| 国产男人的电影天堂91| 久久久久久久亚洲中文字幕| 深夜a级毛片| 三级毛片av免费| 精品熟女少妇av免费看| 99久久无色码亚洲精品果冻| 欧美高清成人免费视频www| 日日撸夜夜添| 久久久久久久久久久丰满| 国产三级在线视频| 99热这里只有是精品在线观看| 22中文网久久字幕| 99riav亚洲国产免费| 精品少妇黑人巨大在线播放 | 亚洲精品日韩在线中文字幕 | 欧美不卡视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色视频三级网站网址| 国产高清有码在线观看视频| 赤兔流量卡办理| 亚州av有码| 蜜桃久久精品国产亚洲av| www日本黄色视频网| 最近2019中文字幕mv第一页| 99久久久亚洲精品蜜臀av| 一边摸一边抽搐一进一小说| 免费av不卡在线播放| 天美传媒精品一区二区| 内射极品少妇av片p| 国产成人福利小说| 可以在线观看毛片的网站| 男女做爰动态图高潮gif福利片| 亚洲不卡免费看| 欧美日韩在线观看h| 嘟嘟电影网在线观看| 亚洲精品乱码久久久久久按摩| 中文字幕精品亚洲无线码一区| 寂寞人妻少妇视频99o| 久久久久网色| 国产熟女欧美一区二区| 一级毛片我不卡| 欧美激情久久久久久爽电影| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看| 一级二级三级毛片免费看| 99久久精品一区二区三区| 综合色av麻豆| 国产日韩欧美在线精品| 久久精品影院6| 丝袜喷水一区| 欧美精品国产亚洲| 国产在视频线在精品| 亚洲一级一片aⅴ在线观看| 亚洲人成网站在线播| 男女下面进入的视频免费午夜| 成年免费大片在线观看| 精华霜和精华液先用哪个| 夜夜夜夜夜久久久久| 免费观看a级毛片全部| 桃色一区二区三区在线观看| 免费看a级黄色片| 99热全是精品| 亚洲欧美日韩无卡精品| 99热6这里只有精品| av在线观看视频网站免费| 国产黄色视频一区二区在线观看 | 国产白丝娇喘喷水9色精品| 校园人妻丝袜中文字幕| 天天一区二区日本电影三级| 一级毛片久久久久久久久女| 久久久久国产网址|