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

    時序影像在冬小麥種植區(qū)提取中的應(yīng)用分析

    2021-12-12 12:09:30明艷芳劉春秀隋淞蔓
    無線電工程 2021年12期
    關(guān)鍵詞:冬小麥時序作物

    石 嫻,明艷芳,劉春秀,瞿 渝,隋淞蔓

    (山東科技大學(xué) 測繪與空間信息學(xué)院,山東 青島 266590)

    0 引言

    冬小麥為頭年秋季播種次年初夏收獲的一種農(nóng)作物[1],在全世界范圍內(nèi)都有著廣闊的種植面積[2]。它的種植面積對經(jīng)濟(jì)和社會發(fā)展有重要的影響[3-4]。精確、高效的冬小麥空間分布監(jiān)測具有十分重要的意義[5]。遙感技術(shù)能夠快速實(shí)現(xiàn)大面積作物種植區(qū)域的連續(xù)監(jiān)測,可以在冬小麥監(jiān)測中發(fā)揮重要作用[6]。

    農(nóng)作物的遙感識別方法主要包括2大類:單時相遙感影像識別和多時相遙感變化檢測。然而由于不同植被之間具有較為相似的光譜特征,單一的光譜信息難以達(dá)到較高的識別精度,而物候期的差異通常被作為農(nóng)作物識別的重要信息。此外,區(qū)域灌溉不均勻、土壤鹽漬化、病蟲害、植被密度差異甚至植物的含水量不同等因素都會導(dǎo)致相同作物在生長過程中產(chǎn)生不同程度的差異[7],因此很難利用單一時相的遙感影像進(jìn)行區(qū)分。多時相遙感影像能夠充分反映不同作物的物候特性及變化規(guī)律,增大光譜特征相似作物之間的可分離性,可以很好地解決“異物同譜”和“同物異譜”問題,大大提高識別精度[8]。隨著高時空分辨率遙感的快速發(fā)展,基于時間序列的農(nóng)作物遙感識別已廣泛應(yīng)用于作物種植面積提取、農(nóng)業(yè)遙感監(jiān)測等方面。

    時間序列遙感作物識別方法主要利用作物整個生育期固定時間間隔的時序數(shù)據(jù)獲取農(nóng)作物特定的時序曲線,然后分析其物候或數(shù)學(xué)特征并進(jìn)行農(nóng)作物識別。楊小喚等[9]繪制并分析了作物的歸一化植被指數(shù)(NDVI)時序曲線,制定了不同作物的提取依據(jù),獲取了北京市的冬小麥、玉米和大豆的種植空間分布圖;陳健等[10]以整個作物生育期內(nèi)35個時相的合成地表反射率數(shù)據(jù)生成MODIS數(shù)據(jù)的增強(qiáng)植被指數(shù)(EVI)時序數(shù)據(jù),結(jié)合作物物候歷和種植結(jié)構(gòu),構(gòu)建分類模型,最終達(dá)到了95.7%的整體分類精度。近年來隨著遙感衛(wèi)星研發(fā)技術(shù)的不斷進(jìn)步,云平臺計算能力和計算機(jī)圖像處理技術(shù)的不斷提高,作物識別逐漸向大空間尺度、長時間序列方向發(fā)展,分類算法也在不斷改進(jìn)。周珂等[11]根據(jù)每旬NDVI最大值合成Landsat8影像,提取了河南省的冬小麥,其平均分類精度達(dá)到了95.3%,且有效降低了與統(tǒng)計數(shù)據(jù)的相對誤差;You等[12]將3個波段、7個光譜指數(shù)作為特征,使用10 d中值合成Sentinel-2時間序列影像,利用隨機(jī)森林分類器進(jìn)行分類,首次在10 m分辨率下制作了中國年度區(qū)域作物圖。

    然而,目前大多數(shù)識別方法都是使用某一種時間序列數(shù)據(jù),很少有專門的研究通過試驗(yàn)比較分析不同間隔的時序影像用來做冬小麥提取的適宜性和準(zhǔn)確性。根據(jù)Blaes等[13]的研究,時間序列影像中過長的時間間隔或關(guān)鍵期一景影像的缺失都會對分類精度造成影響。Wardlow等[14]表示,每種作物有其獨(dú)特的時序光譜曲線,不同的構(gòu)造方法會造成時序數(shù)據(jù)的變化,使得作物的提取效果產(chǎn)生差異。由此,時序數(shù)據(jù)的構(gòu)造方法對于時序影像的遙感識別分類和基于物候特征的遙感監(jiān)測產(chǎn)生著至關(guān)重要的影響,此類研究也十分有必要開展。

    另外,由于中等分辨率衛(wèi)星自身的局限性以及小規(guī)模農(nóng)業(yè)種植區(qū)域的分類易受到背景因素的干擾,像素尺度下的冬小麥提取效果可能難以滿足要求。面向地塊的分類方法在以往的研究中被證明可以很好地克服上述問題,尤其是在小規(guī)模農(nóng)業(yè)盛行的地區(qū)要優(yōu)于面向像素的分類方法[15-16]。但是目前還未有在不同時間序列數(shù)據(jù)下關(guān)于地塊和像素提取方法的對比,也未有關(guān)于提取效率方面的評價。

    綜上所述,冬小麥種植區(qū)提取目前主要存在時間序列單一的問題,另外沒有考慮到面向地塊、像素方法與不同時序數(shù)據(jù)結(jié)合的適用性。因此本文利用GEE平臺和Sentinel-2A/B數(shù)據(jù)評估了不同時間間隔(10、15和30 d)合成的時序數(shù)據(jù)對冬小麥提取精度的影響,明確了更適合提取冬小麥種植區(qū)的Sentinel-2A/B時序數(shù)據(jù)。同時對比了使用不同間隔的時序數(shù)據(jù)時面向地塊和面向像素分類方法的優(yōu)劣性,為擴(kuò)大Sentinel-2A/B數(shù)據(jù)在作物分類中的應(yīng)用提供了依據(jù)。

    1 研究區(qū)和數(shù)據(jù)預(yù)處理

    1.1 研究區(qū)概況

    本文研究區(qū)膠州市(36°00′N~36°30′N,119°37′E~120°12′E)位于山東省青島市,占地面積為1 324 km2,總耕地面積82萬畝,地貌類型有丘陵、平原、洼地和濱海低地4類。膠州市年無霜期210 d,平均氣溫12 ℃,有效積溫4 599.2 ℃,日照時數(shù)2 170.5 h,降水量724.8 mm,為典型的暖溫帶半濕潤季風(fēng)區(qū)大陸性氣候,四季分明,雨熱同期,非常適宜小麥、馬鈴薯和花生等農(nóng)作物的種植。

    1.2 作物物候期

    通過實(shí)地調(diào)查了解到研究區(qū)內(nèi)農(nóng)作物多種多樣,主要種植的農(nóng)作物有冬小麥、馬鈴薯、花生和玉米等。該地玉米年耕種2次,冬小麥、馬鈴薯和花生則耕種一次,典型地表類型包括冬小麥-夏玉米雙熟制耕地、馬鈴薯-夏玉米雙熟制耕地、蔬菜-春玉米雙熟制耕地和花生單熟制耕地。冬小麥一般在10月初播種,12月下旬進(jìn)入越冬期,翌年3月開始返青生長,4月進(jìn)入生長旺期,6月中旬以前基本收割完畢,生長期大概8個月。本研究將2020年10月—2021年6月(共270 d)作為研究期,使用DOW(Day of Wheat)表示研究期內(nèi)的時間,表1為研究期內(nèi)膠州市主要農(nóng)作物的關(guān)鍵物候期。

    表1 膠州市主要農(nóng)作物物候期

    1.3 衛(wèi)星數(shù)據(jù)及預(yù)處理

    Sentinel-2是一種高分辨率多光譜成像衛(wèi)星,攜帶有一枚多光譜成像儀(MSI),包括Sentinel-2A和Sentinel-2B兩顆衛(wèi)星,一起提供5 d間隔的影像。GEE平臺提供了2種級別的Sentinel-2A/B數(shù)據(jù):Level-1C和Level-2A。Level-1C級產(chǎn)品是指經(jīng)過了輻射定標(biāo)、幾何校正(主要包括正射校正、空間配準(zhǔn)操作)的大氣表層(Top-of-Atmosphere,TOA)反射率數(shù)據(jù),Level-2A級產(chǎn)品則是在Level-1C級產(chǎn)品的基礎(chǔ)上又經(jīng)過了大氣校正的地表反射率(Surface Reflectance,SR)數(shù)據(jù)[17]。本次研究使用Sentinel-2 SR數(shù)據(jù)集,已經(jīng)消除了大氣吸收以及各項散射作用造成的誤差,能夠反演地物真實(shí)的表面反射率,使得作物分類更加準(zhǔn)確。數(shù)據(jù)共含有12個光譜波段,分別具有不同的空間分辨率(10,20,60 m),具體的數(shù)據(jù)介紹如表2所示。

    表2 Sentinel-2A/B Level-2A數(shù)據(jù)介紹

    研究期內(nèi)共218景影像,由于原始數(shù)據(jù)是按10 000倍縮放的SR數(shù)據(jù),將每景影像除以10 000獲得真實(shí)單位的地表反射率數(shù)據(jù)。另外,數(shù)據(jù)還有3個QA(Quality Assessment)波段,其中的QA60波段存儲有云掩膜信息,利用它消除Sentinel-2A/B數(shù)據(jù)中被云和陰影污染的觀測,又使用線性插值的方法填補(bǔ)去云留下的缺失,以實(shí)現(xiàn)整個時間、空間域的完全覆蓋。

    1.4 地面調(diào)查數(shù)據(jù)

    樣本數(shù)據(jù)來源于2021年3月和5月膠州市野外實(shí)地調(diào)查數(shù)據(jù),共獲得436個樣本點(diǎn),其中冬小麥212個,非冬小麥224個。采樣時記錄樣本的經(jīng)緯度及作物類別,并拍攝照片方便后續(xù)核查。

    2 材料和方法

    本研究主要包括以下2個部分。首先,利用圖像分割技術(shù)從高分辨率影像中提取膠州市的地塊,為了避免非農(nóng)作物區(qū)域的干擾,將非耕地地塊剔除。然后,在GEE平臺中使用實(shí)地調(diào)查得到的樣本結(jié)合構(gòu)造的作物時序特征作為輸入,訓(xùn)練基于RF算法的冬小麥分類器,將分類器分別應(yīng)用于耕地地塊或像素,以實(shí)現(xiàn)膠州市冬小麥的提取。研究的技術(shù)路線如圖1所示。

    圖1 技術(shù)路線Fig.1 Technical route

    2.1 耕地地塊獲取

    在面向地塊的分類方法中地塊的獲取過程至關(guān)重要,地塊與真實(shí)地物的符合程度直接影響信息提取的正確與否。圖像分割的過程其實(shí)是相鄰?fù)|(zhì)像素的結(jié)合和異質(zhì)像素的分離過程,在盡量減少圖像信息損失的基礎(chǔ)上將圖像分割成有意義的多邊形對象[18]。每個多邊形代表一個實(shí)地地物,內(nèi)部具有相同光譜、紋理等信息。

    本文使用易康(eCognition)軟件的多尺度分割方法獲得耕地地塊,影像為下載得到的2021年4月膠州市GF2 PMS數(shù)據(jù),經(jīng)過正射校正、波段融合后,空間分辨率達(dá)到了1 m。為了使分割結(jié)果與影像更加吻合,經(jīng)過多次試驗(yàn)達(dá)到了最佳效果,分割尺度設(shè)置為50,形狀因子為 0.3,緊致度為0.4。采用目視解譯的方法將非耕地地塊(主要包括建筑物、道路等)刪除,只在耕地范圍內(nèi)作圖,人工修正和細(xì)化邊界后得到了99 733個農(nóng)作物地塊,局部效果如圖2所示。近年來國家耕地保護(hù)政策十分嚴(yán)格,獲得的耕地塊矢量具有較長期的穩(wěn)定性,可以重復(fù)利用。

    圖2 地塊分割Fig.2 Parcel segmentation picture

    2.2 隨機(jī)森林分類

    隨機(jī)森林(RF)算法是Breiman在2001年提出的一種強(qiáng)大的機(jī)器學(xué)習(xí)算法[19],已經(jīng)在不同領(lǐng)域有了廣泛的應(yīng)用,例如風(fēng)險評估及預(yù)測、醫(yī)學(xué)檢驗(yàn)和金融投資等[20-22]。它是基于決策樹弱分類器的集成學(xué)習(xí)算法,在利用集成算法優(yōu)勢的同時規(guī)避了單個分類器的缺點(diǎn),比許多傳統(tǒng)的分類器(如單層神經(jīng)網(wǎng)絡(luò)、最大似然法等)更具有準(zhǔn)確性和魯棒性[23]。RF的隨機(jī)主要有2個方面:首先是對原始數(shù)據(jù)集采用bootstrap方法(隨機(jī)且有放回的抽樣方法)抽取多個訓(xùn)練集和其對應(yīng)的測試集;其次是在生成決策樹時隨機(jī)抽取一部分特征用于樹節(jié)點(diǎn)的分裂。最后每個決策樹進(jìn)行投票來確定最終分類結(jié)果。

    RF算法的模型表示為:

    (1)

    式中,I表示示性函數(shù);Y表示預(yù)測結(jié)果;hi(Si)表示單顆決策樹。

    RF分類步驟如下:

    ① 從原始數(shù)據(jù)集N中隨機(jī)且有放回地抽取n個樣本,重復(fù)n次,而未抽中的樣本組成袋外數(shù)據(jù)集(OOB),作為測試數(shù)據(jù);

    ② 從特征集合M中隨機(jī)抽取m個作為特征子集,且m

    ③ 基尼指數(shù)能夠選出數(shù)據(jù)集的最佳分界點(diǎn),找出使每棵決策樹的基尼指數(shù)最小的特征作為最優(yōu)特征,其對應(yīng)的切分點(diǎn)作為最優(yōu)切分點(diǎn),然后將剩余訓(xùn)練數(shù)據(jù)應(yīng)用到下面的2個子節(jié)點(diǎn)中;

    ④ 在每個子節(jié)點(diǎn)上重復(fù)進(jìn)行步驟③,直到分割完所有節(jié)點(diǎn);

    ⑤ 重復(fù)進(jìn)行步驟②~步驟④,得到需要數(shù)量的決策樹;

    ⑥ 每一顆決策樹都會對測試數(shù)據(jù)集中的每一條數(shù)據(jù)進(jìn)行投票,統(tǒng)計得票最多的類別即為該樣本的最終類別。

    RF算法有著較高的準(zhǔn)確性,GEE平臺中的RF算法已成功用于檢測土地利用變化和作物分類等研究中,本研究將它運(yùn)用在冬小麥提取中。

    2.2.1 光譜特征計算

    NDVI是監(jiān)測作物生長狀態(tài)的最佳指示因子,還能反映出植物冠層的背景(如土壤、潮濕地面、雪和枯葉等)影響,當(dāng)土地還未耕作時裸露地表的NDVI指數(shù)值接近0。NDVI計算公式為:

    (2)

    式中,ρNIR為近紅外波段的反射率;ρR為紅波段的反射率。

    EVI設(shè)置了更窄的紅邊波段,可以減少水汽的影響,加強(qiáng)了對植被稀疏地區(qū)的監(jiān)測能力。另外,它引入的藍(lán)光波段能夠降低土壤背景和大氣噪聲的干擾,可穩(wěn)定地反映地表植被特征。EVI計算公式為:

    (3)

    式中,ρNIR為近紅外波段的反射率;ρR為紅波段的反射率;ρB為藍(lán)波段的反射率。

    2.2.2 時序數(shù)據(jù)分析及合成

    為了對比不同時間間隔的合成對時序曲線的影響,本文結(jié)合Sentinel-2A/B NDVI,EVI時間序列數(shù)據(jù)和實(shí)地調(diào)查的樣本數(shù)據(jù),繪制不同作物原始的NDVI,EVI時間序列曲線。由于受到云、水汽、氣溶膠和傳感器等因素的干擾,時間序列數(shù)據(jù)容易出現(xiàn)異常值,時間序列曲線出現(xiàn)不規(guī)則波動現(xiàn)象,對波譜分析和特征構(gòu)造造成不可忽略的影響[24]。平滑和去噪操作可以更加真實(shí)地反映各種作物的生長曲線,本文采用S-G(Savitzky-Golay)濾波方法重構(gòu)Sentinel-2A/B的時間序列數(shù)據(jù),設(shè)置移動窗口大小為7(70 d觀測值),多項式次數(shù)為3。圖3為經(jīng)過平滑后的作物NDVI,EVI時間序列曲線,橫坐標(biāo)表示日期,縱軸表示NDVI,EVI的值。

    (a) NDVI時序曲線

    由NDVI時間序列曲線可以看出,冬小麥曲線整體呈現(xiàn)“兩峰一谷”的趨勢。夏玉米在9月下旬—10月上旬收獲,然后冬小麥開始播種,7 d左右出苗,NDVI逐漸增大;12月份出現(xiàn)一個小的波峰,然后開始冬眠,光合作用衰減,NDVI逐漸減??;波谷在1月份出現(xiàn);在2—3月份開始返青、起身,NDVI迅速升高;4月中下旬—5月上旬進(jìn)入抽穗期,達(dá)到第2個波峰;6月成熟后收獲,植被指數(shù)值明顯下降。花生、春玉米和馬鈴薯都是春季作物,一般去年的夏季作物收獲之后都會為其保留耕地,NDVI值較低,在11月中旬—次年4月都與冬小麥曲線有著明顯差異。EVI曲線與NDVI較為類似,冬小麥曲線也有著“兩峰一谷”的形態(tài),在11月中旬—次年4月都是區(qū)分冬小麥的關(guān)鍵時期。

    由于影像眾多,且直接使用單景影像易受到異常值的干擾,雖然進(jìn)行了濾波處理但是作物曲線依然出現(xiàn)了振蕩現(xiàn)象。本文分別取DOW 1~270中的每10,15,30 d的均值合成一景新的影像,得到NDVI,EVI的不同時間間隔合成的時序影像。面向地塊的分類方法最后還要計算每個地塊內(nèi)的均值作為地塊的特征。

    2.2.3 特征選擇

    經(jīng)過以上計算,得到了不同時間間隔合成的時序特征,如表3所示。

    表3 時序特征

    多光譜時序數(shù)據(jù)在提供更多光譜信息的同時,也帶來了高維輸入和輸出的新挑戰(zhàn)。由于大量光譜特征可能會攜帶高度相關(guān)的信息,增加模型的計算負(fù)擔(dān),出現(xiàn)過擬合現(xiàn)象。因此特征選擇顯得十分重要,它在很大程度上決定機(jī)器學(xué)習(xí)算法的效率。RF允許使用基尼指數(shù)在分類時對特征重要性進(jìn)行評價。因此使用GEE中的評價方法對特征波段進(jìn)行排序,最后分別保留14~24個最佳特征進(jìn)行冬小麥的提取。

    2.2.4 模型訓(xùn)練和分類

    在進(jìn)行地塊級別的分類時,由于膠州市面積較大,地塊數(shù)量較多,將它分為3個區(qū)分別計算和導(dǎo)出結(jié)果。面向地塊的分類方法將冬小麥樣本所屬的地塊標(biāo)記為1,非冬小麥地塊標(biāo)記為0;面向像素的分類方法將冬小麥樣本點(diǎn)標(biāo)記為1,非冬小麥點(diǎn)標(biāo)記為0。最后將樣本分別按照7∶3的比例隨機(jī)分為訓(xùn)練樣本和驗(yàn)證樣本。將樣本和特征都導(dǎo)入RF分類器中進(jìn)行分類,設(shè)置RF的參數(shù)如下:① numberOfTrees,樹數(shù)代表建立RF模型的決策樹的數(shù)量,在一定范圍內(nèi)隨著樹數(shù)的增多,精度可能會略有提高,設(shè)置為500;② seed,隨機(jī)種子設(shè)置為999;③ 其他幾個參數(shù),maxNodes(每棵樹中葉節(jié)點(diǎn)的最大數(shù)目)、minLeafPopulation(葉節(jié)點(diǎn)所需的最小樣本數(shù))、bagFraction(每棵樹輸入到bag的分?jǐn)?shù))、variablesPerSplit(每次分割的變量數(shù))保持默認(rèn)設(shè)置。

    為了方便與面向像素的分類方法進(jìn)行對比,基于像素的冬小麥提取也使用多尺度分割得到的耕地地塊進(jìn)行掩膜,時序特征及樣本也與面向地塊的分類方法保持一致。

    2.3 結(jié)果評價指標(biāo)

    評價指標(biāo)用來定量的衡量一個算法性能的優(yōu)劣水平?;煜仃囀亲畛S玫姆绞?,能一目了然地展示出分類結(jié)果和地表真實(shí)信息的差異,如表4所示。

    表4 混淆矩陣

    將真實(shí)類別為非小麥的樣本稱為負(fù)類,真實(shí)類別為小麥的樣本稱為正類,TN(True Negative)代表負(fù)類的正確預(yù)測,F(xiàn)N(False Negative)代表正類的錯誤預(yù)測,F(xiàn)P(False Positive)代表負(fù)類的錯誤預(yù)測,TP(True Positive)代表正類的正確預(yù)測。樣本總和為N,則:

    N=TN+FN+FP+TP。

    (4)

    常用的幾種精度評價方式:總體分類精度(OA)、Kappa系數(shù)、生產(chǎn)者精度(PA)和用戶精度(UA)都是基于混淆矩陣計算的。OA表示正確分類的像元(地塊)總和除以總像元(地塊)數(shù),計算公式如下:

    (5)

    但是OA只考慮分類正確的比例,不能對模型性能的好壞做出評價。Kappa系數(shù)不僅考慮了被正確分類的像素(地塊),又綜合了各種錯分和漏分誤差,是更為全面的評價指標(biāo)。Kappa系數(shù)計算公式為:

    (6)

    式中,

    (7)

    (8)

    生產(chǎn)者精度是指分類器將所有像素(地塊)正確分為某類和某類真實(shí)總數(shù)的比例,例如冬小麥類別的生產(chǎn)者精度計算公式為:

    (9)

    用戶精度是指正確分到某類的像素(地塊)數(shù)與分類器分為某類的總數(shù)的比例,例如冬小麥類別的用戶精度計算公式為:

    1.2.2 培訓(xùn)方法 我院臨床技能培訓(xùn)中心采用兒童仿真模型模擬臨床環(huán)境對學(xué)員進(jìn)行培訓(xùn),內(nèi)容包括理論及技能操作培訓(xùn),形式為教師講解、視頻演示、邊看邊做、實(shí)戰(zhàn)演練等。

    (10)

    3 結(jié)果與評價

    3.1 時序曲線對比

    圖4為冬小麥生長期內(nèi)(DOW 0~270)不同時間間隔合成的作物生長曲線圖,橫軸表示時間,縱軸表示NDVI,EVI的值。

    (a) 10 d_NDVI

    經(jīng)過時間合成以后,曲線變得更加平滑也消除了一些異常波動,并且合成的時間間隔越長,作物曲線越平滑,但同時也失去了更多的細(xì)節(jié),比如30 d合成的冬小麥NDVI,EVI時序曲線均不再呈現(xiàn)“兩峰一谷”的形態(tài)。時間間隔為10 d和15 d的曲線既抵消了一些干擾也保留了作物生長的細(xì)節(jié)變化,在不同特征的曲線中冬小麥與其他作物的差異都比較明顯。

    3.2 定性評價

    為了檢驗(yàn)本文方法的有效性,對比不同時序數(shù)據(jù)以及地塊和像素級別分類的優(yōu)劣性,分別使用定性和定量的方式評價最終分類結(jié)果。定量方法即使用第2節(jié)中的評價指標(biāo)進(jìn)行對比分析,定性方法主要將提取結(jié)果與遙感影像進(jìn)行對比,目視評判冬小麥提取效果,主要包括下面幾個方面:① 判斷冬小麥提取區(qū)域是否全面;② 判斷是否存在誤提現(xiàn)象;③ 基于像素的提取方法孤立像素現(xiàn)象是否嚴(yán)重。

    本文使用多尺度分割時利用的GF2 1m分辨率影像與各種方法得到的提取結(jié)果進(jìn)行局部對比,圖5為典型區(qū)域的原始影像,冬小麥表現(xiàn)為深綠色,圖6、圖7和圖8分別為這些區(qū)域?qū)?yīng)的提取結(jié)果。

    (a) 規(guī)整區(qū)域

    (a) 30 d_像元

    (a) 30 d_像元

    (a) 30 d_像元

    在圖6中,30 d和15 d合成的地塊、像素級別的提取均存在一定程度的誤提現(xiàn)象;10 d合成的像素級別的提取有少部分的漏提;10 d合成的地塊級別的提取結(jié)果符合影像顯示的冬小麥區(qū)域。在圖7中,地塊的不規(guī)整增大了提取的難度,基于像素的提取結(jié)果漏提現(xiàn)象嚴(yán)重;30 d和15 d合成的地塊級提取也出現(xiàn)了漏提現(xiàn)象;10 d合成的地塊級提取結(jié)果較為全面準(zhǔn)確。圖8中,在復(fù)雜種植區(qū)域,基于像素的提取結(jié)果漏提和孤立像素現(xiàn)象都比較嚴(yán)重,提取效果不佳;30 d合成的地塊級提取有小面積的漏提,15 d和10 d合成的地塊級別的提取結(jié)果一致,在復(fù)雜種植區(qū)域的提取結(jié)果也較為理想。整體來看,地塊和像素級別的提取效果都會隨著合成時間序列數(shù)據(jù)時間間隔的縮小變得更好;當(dāng)時間間隔相同時,地塊級別的提取并不全都優(yōu)于像素級別的提取,但隨著時間間隔的縮小地塊級別的提取進(jìn)步明顯,10 d合成的面向地塊的提取方法總是最優(yōu)的,幾乎沒有出現(xiàn)錯提、漏提現(xiàn)象。

    3.3 定量評價

    表5為最終得到的一系列量化結(jié)果,包括不同提取方法得到的OA、Kappa系數(shù)、冬小麥的PA和UA。

    表5 冬小麥提取精度

    以上方法均得到了較高的精度,證明了利用規(guī)律時間間隔合成的時間序列影像進(jìn)行分類的有效性。整體而言,與基于像素的分類結(jié)果相比,面向地塊的分類方法都得到了更高的OA和Kappa系數(shù),但PA和UA并不全都高于基于像素的分類。并且總體趨勢顯示,當(dāng)合成Sentinel-2A/B影像的時間間隔越短時,冬小麥提取的精度越高。15 d和10 d合成時間間隔相差不大,提取精度相差不大。

    3.4 效率評價

    表6為各種冬小麥提取方法所需要的計算時間和結(jié)果導(dǎo)出的時間。

    表6 冬小麥提取時間對比

    無論是計算時間還是導(dǎo)出時間,不同時間序列下地塊級別的提取之間或者像素級別的提取之間都相差不大,真正對提取時間產(chǎn)生影響的不是不同時間序列的合成方法,而是面向不同元素的分類方法。地塊與像素級別的提取在計算時間上最多相差3倍,但都在1 min之內(nèi),因此對計算效率的對比評價意義不大;然而在結(jié)果的導(dǎo)出上,二者相差了23 h,將近一天的時間,由于面向像素的分類方法精度也在比較高的水平,當(dāng)在時間緊急的情況下,面向地塊的分類方法優(yōu)勢不大。

    4 冬小麥分布

    通常在種植結(jié)構(gòu)單一、農(nóng)田連片區(qū)域的冬小麥更容易識別,但在地塊破碎區(qū)域由于內(nèi)部的光譜變異和邊界的光譜混合導(dǎo)致識別較為困難。本文的研究區(qū)為整個膠州市,地塊數(shù)量眾多、類型多樣,方法具有普適性。由以上對比結(jié)果得知,10 d合成的面向地塊的分類方法精度更高,提取效果更好,因此利用此方法做出2021年的膠州市冬小麥種植分布圖,如圖9所示,可以明顯看出,冬小麥種植區(qū)主要分布在膠州市的北部和西南部。

    圖9 2021年膠州市冬小麥空間分布Fig.9 Spatial distribution map of winter wheat in Jiaozhou in 2021

    5 結(jié)束語

    本文使用不同時間間隔合成的Sentinel-2A/B時間序列數(shù)據(jù),提取了膠州市的冬小麥種植區(qū),主要為了探索更適合提取冬小麥的時序數(shù)據(jù)與方法。為了對比在不同時序數(shù)據(jù)下面向像素和面向地塊的分類結(jié)果,2種方法在完全相同的條件下進(jìn)行。結(jié)果表明:① 當(dāng)Sentinel-2A/B合成影像的時間間隔縮短時,可以為冬小麥與其他作物的區(qū)分提供更多信息,地塊和像素級別的提取準(zhǔn)確度都會更高;② 在相同的時序數(shù)據(jù)下,面向地塊的分類方法并不總是優(yōu)于面向地塊的分類方法,但可以明顯改善面向像素分類結(jié)果中常見的“椒鹽現(xiàn)象”,且隨著合成時間序列數(shù)據(jù)時間間隔的縮小進(jìn)步更明顯。同時,更短時間間隔(10 d)合成的Sentinel-2時序數(shù)據(jù)結(jié)合面向地塊的分類方法,對冬小麥的提取幫助更大,在一些不規(guī)整、復(fù)雜種植區(qū)域也獲得了很好的提取效果;③ 由于地塊級別的結(jié)果導(dǎo)出時間較長,當(dāng)對精度和提取效果要求較高而時間比較充裕時,選擇地塊級別的分類更加合適;當(dāng)對時間的要求較高而對精度的需求較低時,面向像素的分類更加符合條件。分類是不斷優(yōu)化的過程,隨著遙感衛(wèi)星研發(fā)技術(shù)、計算機(jī)和云平臺計算能力的不斷提高,現(xiàn)有的分類算法也在逐漸改進(jìn),未來更高級的圖像分割技術(shù)與更強(qiáng)大的云平臺計算能力的結(jié)合將獲得更加準(zhǔn)確的冬小麥種植區(qū)提取結(jié)果。

    猜你喜歡
    冬小麥時序作物
    時序坐標(biāo)
    基于Sentinel-2時序NDVI的麥冬識別研究
    作物遭受霜凍該如何補(bǔ)救
    四種作物 北方種植有前景
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    甘肅冬小麥田
    一種毫米波放大器時序直流電源的設(shè)計
    電子制作(2016年15期)2017-01-15 13:39:08
    無人機(jī)遙感在作物監(jiān)測中的應(yīng)用與展望
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    冬小麥——新冬18號
    黑人猛操日本美女一级片| 成人亚洲精品一区在线观看| 免费人妻精品一区二区三区视频| 欧美精品av麻豆av| 午夜两性在线视频| av在线老鸭窝| 99精国产麻豆久久婷婷| 女人精品久久久久毛片| 脱女人内裤的视频| 一级黄色大片毛片| 另类亚洲欧美激情| 亚洲精品久久久久久婷婷小说| 久久久久国产一级毛片高清牌| 91精品伊人久久大香线蕉| 天堂俺去俺来也www色官网| 久久亚洲精品不卡| 久久人妻熟女aⅴ| 人人妻人人澡人人看| 欧美大码av| 丝袜脚勾引网站| 国产日韩欧美视频二区| 国产成人系列免费观看| 日本黄色日本黄色录像| 亚洲国产精品一区二区三区在线| 免费久久久久久久精品成人欧美视频| 日韩 亚洲 欧美在线| 亚洲第一av免费看| cao死你这个sao货| 国产亚洲午夜精品一区二区久久| 老司机亚洲免费影院| 精品亚洲成国产av| 色网站视频免费| 视频区图区小说| 亚洲精品久久午夜乱码| 欧美日韩亚洲高清精品| 国产在线视频一区二区| 香蕉国产在线看| 制服人妻中文乱码| 久久人人爽人人片av| 两个人免费观看高清视频| 国产福利在线免费观看视频| 亚洲精品国产一区二区精华液| 不卡av一区二区三区| 91成人精品电影| 90打野战视频偷拍视频| 国产成人av教育| 高清不卡的av网站| 中文字幕亚洲精品专区| 色网站视频免费| 99久久精品国产亚洲精品| 少妇猛男粗大的猛烈进出视频| 黄色 视频免费看| 亚洲男人天堂网一区| 丝袜人妻中文字幕| 50天的宝宝边吃奶边哭怎么回事| 免费看av在线观看网站| 看免费成人av毛片| 色婷婷av一区二区三区视频| 色婷婷久久久亚洲欧美| 午夜免费男女啪啪视频观看| 两性夫妻黄色片| 欧美久久黑人一区二区| 午夜两性在线视频| 深夜精品福利| 午夜老司机福利片| 99国产精品99久久久久| 无遮挡黄片免费观看| 一级毛片黄色毛片免费观看视频| 国产成人免费观看mmmm| 下体分泌物呈黄色| 深夜精品福利| 99热网站在线观看| 亚洲欧美中文字幕日韩二区| 亚洲,欧美精品.| 一区二区三区精品91| 青春草亚洲视频在线观看| 一二三四在线观看免费中文在| avwww免费| 欧美精品一区二区大全| 欧美久久黑人一区二区| 97精品久久久久久久久久精品| av在线老鸭窝| 最新的欧美精品一区二区| 男女边吃奶边做爰视频| 日韩一区二区三区影片| 国产精品九九99| 国产国语露脸激情在线看| 91精品伊人久久大香线蕉| 中文字幕制服av| 日韩,欧美,国产一区二区三区| 99精品久久久久人妻精品| 男男h啪啪无遮挡| 乱人伦中国视频| 丝袜美腿诱惑在线| 多毛熟女@视频| 亚洲精品第二区| a级片在线免费高清观看视频| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 晚上一个人看的免费电影| 九草在线视频观看| av又黄又爽大尺度在线免费看| 亚洲欧美中文字幕日韩二区| 日本黄色日本黄色录像| av又黄又爽大尺度在线免费看| 下体分泌物呈黄色| 亚洲,欧美精品.| 69精品国产乱码久久久| 又黄又粗又硬又大视频| 亚洲视频免费观看视频| 人人妻人人澡人人看| 91九色精品人成在线观看| 波多野结衣av一区二区av| 国产福利在线免费观看视频| 男女国产视频网站| 成年av动漫网址| 欧美日韩亚洲综合一区二区三区_| 一本—道久久a久久精品蜜桃钙片| 日本欧美国产在线视频| 欧美黄色片欧美黄色片| 婷婷色综合大香蕉| 欧美性长视频在线观看| 午夜福利乱码中文字幕| 99精品久久久久人妻精品| 精品国产国语对白av| 波多野结衣av一区二区av| 满18在线观看网站| 欧美日韩视频高清一区二区三区二| 国产高清视频在线播放一区 | 亚洲精品日韩在线中文字幕| 少妇猛男粗大的猛烈进出视频| 国产亚洲av片在线观看秒播厂| 国产成人精品无人区| 男女国产视频网站| 自线自在国产av| 久久人妻福利社区极品人妻图片 | 欧美黄色片欧美黄色片| 操美女的视频在线观看| 国产视频一区二区在线看| 91精品三级在线观看| 国产极品粉嫩免费观看在线| 男人爽女人下面视频在线观看| kizo精华| e午夜精品久久久久久久| 免费在线观看视频国产中文字幕亚洲 | 日韩电影二区| 亚洲熟女精品中文字幕| 天天躁日日躁夜夜躁夜夜| 一区福利在线观看| 久久精品久久久久久噜噜老黄| 一本综合久久免费| 人妻 亚洲 视频| 欧美成狂野欧美在线观看| 丝袜美足系列| 久久狼人影院| 免费人妻精品一区二区三区视频| av电影中文网址| 久久精品熟女亚洲av麻豆精品| 丝袜人妻中文字幕| 无限看片的www在线观看| 日韩欧美一区视频在线观看| 国产成人av教育| 欧美国产精品va在线观看不卡| 考比视频在线观看| bbb黄色大片| 久久久久网色| 性高湖久久久久久久久免费观看| 亚洲国产中文字幕在线视频| 国产日韩欧美视频二区| 少妇裸体淫交视频免费看高清 | 我的亚洲天堂| 人人妻,人人澡人人爽秒播 | 中文乱码字字幕精品一区二区三区| 午夜免费男女啪啪视频观看| 成人国产av品久久久| 操出白浆在线播放| 久久精品国产综合久久久| 精品视频人人做人人爽| 在线观看免费午夜福利视频| 午夜激情av网站| 91精品三级在线观看| 成人国产av品久久久| 国产人伦9x9x在线观看| 精品亚洲成a人片在线观看| 伊人亚洲综合成人网| 黄色a级毛片大全视频| 国产亚洲av片在线观看秒播厂| 婷婷色麻豆天堂久久| 黄色毛片三级朝国网站| 午夜福利乱码中文字幕| 亚洲人成电影免费在线| 97人妻天天添夜夜摸| 十八禁网站网址无遮挡| 午夜久久久在线观看| 亚洲精品中文字幕在线视频| 最近中文字幕2019免费版| 国产男人的电影天堂91| 亚洲国产欧美日韩在线播放| 母亲3免费完整高清在线观看| 人人澡人人妻人| 涩涩av久久男人的天堂| 亚洲熟女精品中文字幕| 久热爱精品视频在线9| 一二三四社区在线视频社区8| 国产97色在线日韩免费| 丝袜在线中文字幕| 亚洲激情五月婷婷啪啪| 一本综合久久免费| 大香蕉久久网| 男女国产视频网站| 巨乳人妻的诱惑在线观看| 国产精品成人在线| 午夜激情av网站| 国产成人精品在线电影| 久久久久久久精品精品| 亚洲国产av新网站| 97精品久久久久久久久久精品| 视频区欧美日本亚洲| 国产激情久久老熟女| 中文字幕亚洲精品专区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 90打野战视频偷拍视频| 国产高清国产精品国产三级| 国产日韩欧美在线精品| 少妇人妻久久综合中文| 女人被躁到高潮嗷嗷叫费观| 亚洲人成电影观看| av不卡在线播放| 一二三四社区在线视频社区8| 自拍欧美九色日韩亚洲蝌蚪91| 晚上一个人看的免费电影| 大型av网站在线播放| 久久国产亚洲av麻豆专区| 另类亚洲欧美激情| 中文精品一卡2卡3卡4更新| 在线亚洲精品国产二区图片欧美| av网站免费在线观看视频| 91精品国产国语对白视频| 热99久久久久精品小说推荐| 高清av免费在线| 中文字幕制服av| 一级,二级,三级黄色视频| 亚洲av片天天在线观看| 亚洲国产日韩一区二区| 国产精品一区二区在线不卡| tube8黄色片| 亚洲精品日韩在线中文字幕| 欧美性长视频在线观看| 国产av精品麻豆| 色视频在线一区二区三区| 国产又爽黄色视频| 黄色片一级片一级黄色片| 久热爱精品视频在线9| 久久人妻福利社区极品人妻图片 | 免费观看人在逋| 久久精品久久精品一区二区三区| 91麻豆精品激情在线观看国产 | 欧美国产精品va在线观看不卡| 国产一区二区激情短视频 | 国产精品一二三区在线看| 黄频高清免费视频| 午夜日韩欧美国产| 国产在线免费精品| 人人澡人人妻人| 在线观看免费视频网站a站| 老司机午夜十八禁免费视频| 久久久久久人人人人人| 久久九九热精品免费| www日本在线高清视频| 国产亚洲欧美在线一区二区| 人人妻人人添人人爽欧美一区卜| 久久毛片免费看一区二区三区| 国产97色在线日韩免费| 精品亚洲乱码少妇综合久久| 亚洲欧美成人综合另类久久久| 亚洲中文字幕日韩| 美女大奶头黄色视频| 欧美日韩国产mv在线观看视频| 男女床上黄色一级片免费看| 一级片'在线观看视频| 欧美+亚洲+日韩+国产| 最新的欧美精品一区二区| 五月天丁香电影| 亚洲人成77777在线视频| 亚洲欧美日韩高清在线视频 | 亚洲欧美色中文字幕在线| 99热国产这里只有精品6| 亚洲精品成人av观看孕妇| 亚洲av日韩在线播放| 午夜福利视频精品| 亚洲精品美女久久av网站| 亚洲欧洲精品一区二区精品久久久| 久久久国产欧美日韩av| 十分钟在线观看高清视频www| 飞空精品影院首页| 成年人黄色毛片网站| 久久久久国产精品人妻一区二区| 男女下面插进去视频免费观看| 丰满少妇做爰视频| 人人妻人人澡人人爽人人夜夜| 精品一区在线观看国产| 一区福利在线观看| 国产精品免费大片| avwww免费| 爱豆传媒免费全集在线观看| 欧美av亚洲av综合av国产av| 亚洲成国产人片在线观看| 18禁观看日本| 亚洲人成电影免费在线| 在线亚洲精品国产二区图片欧美| 最黄视频免费看| 秋霞在线观看毛片| 久久精品成人免费网站| 国产精品一区二区精品视频观看| 一级片免费观看大全| 欧美精品啪啪一区二区三区 | 久久天躁狠狠躁夜夜2o2o | 成人手机av| 精品少妇一区二区三区视频日本电影| 免费看十八禁软件| 久热爱精品视频在线9| 永久免费av网站大全| 精品久久久精品久久久| 香蕉丝袜av| 色播在线永久视频| 成人亚洲欧美一区二区av| 久久国产亚洲av麻豆专区| 精品人妻熟女毛片av久久网站| 亚洲av电影在线观看一区二区三区| 80岁老熟妇乱子伦牲交| 欧美中文综合在线视频| 国产精品麻豆人妻色哟哟久久| 亚洲欧美成人综合另类久久久| 欧美 日韩 精品 国产| 香蕉国产在线看| 建设人人有责人人尽责人人享有的| 男的添女的下面高潮视频| 在线 av 中文字幕| 永久免费av网站大全| 免费日韩欧美在线观看| 男女免费视频国产| 色婷婷久久久亚洲欧美| 国产精品久久久久久人妻精品电影 | 亚洲av男天堂| 波多野结衣av一区二区av| 七月丁香在线播放| 91成人精品电影| 女人精品久久久久毛片| 丁香六月天网| 婷婷丁香在线五月| 午夜免费观看性视频| 啦啦啦在线观看免费高清www| 精品亚洲成国产av| 国产男人的电影天堂91| 在线观看免费视频网站a站| 高清av免费在线| 国产精品香港三级国产av潘金莲 | 久久精品国产亚洲av涩爱| 黄片播放在线免费| 啦啦啦 在线观看视频| av国产精品久久久久影院| 国产一区二区三区综合在线观看| 一区二区三区四区激情视频| 成人亚洲精品一区在线观看| 人妻人人澡人人爽人人| 日本a在线网址| 亚洲欧美清纯卡通| 99国产精品一区二区蜜桃av | 亚洲熟女精品中文字幕| 午夜精品国产一区二区电影| 菩萨蛮人人尽说江南好唐韦庄| 伊人久久大香线蕉亚洲五| 欧美日韩福利视频一区二区| 亚洲 国产 在线| 久久国产精品影院| 亚洲一区二区三区欧美精品| 一级黄片播放器| 汤姆久久久久久久影院中文字幕| 久久久精品区二区三区| 亚洲天堂av无毛| 老司机靠b影院| 欧美日韩一级在线毛片| 伊人久久大香线蕉亚洲五| 午夜老司机福利片| 无遮挡黄片免费观看| 爱豆传媒免费全集在线观看| 亚洲精品久久成人aⅴ小说| 夫妻午夜视频| 国产91精品成人一区二区三区 | 中国美女看黄片| 老司机在亚洲福利影院| 亚洲精品自拍成人| 国产精品一二三区在线看| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 国产一区二区在线观看av| 亚洲av在线观看美女高潮| 国产精品秋霞免费鲁丝片| 亚洲少妇的诱惑av| 国产一区亚洲一区在线观看| 考比视频在线观看| 国产精品免费大片| 国产欧美日韩一区二区三 | av国产久精品久网站免费入址| 人人妻人人澡人人爽人人夜夜| 99九九在线精品视频| 少妇裸体淫交视频免费看高清 | 国产亚洲av高清不卡| 亚洲精品久久久久久婷婷小说| 久久精品久久久久久久性| 日本av免费视频播放| 国产高清视频在线播放一区 | 一区二区三区乱码不卡18| 久久中文字幕一级| 9热在线视频观看99| 一级片免费观看大全| 日韩一区二区三区影片| 欧美日韩福利视频一区二区| 精品亚洲成国产av| 夫妻午夜视频| www.999成人在线观看| 久久青草综合色| 五月开心婷婷网| 亚洲精品av麻豆狂野| 亚洲久久久国产精品| 欧美少妇被猛烈插入视频| 精品一区二区三卡| 精品久久久久久电影网| 欧美亚洲日本最大视频资源| 18在线观看网站| 九色亚洲精品在线播放| 在线看a的网站| 一本久久精品| 成人亚洲精品一区在线观看| 久久久久久久久免费视频了| 青春草视频在线免费观看| 国产麻豆69| 国产亚洲一区二区精品| 国产熟女欧美一区二区| 日日爽夜夜爽网站| 欧美老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| 亚洲第一青青草原| 国产精品三级大全| 国产成人系列免费观看| 日本vs欧美在线观看视频| 在线亚洲精品国产二区图片欧美| 青草久久国产| 777久久人妻少妇嫩草av网站| 欧美精品人与动牲交sv欧美| 亚洲精品国产区一区二| 国产成人av教育| 亚洲国产毛片av蜜桃av| www.熟女人妻精品国产| 在线天堂中文资源库| 99九九在线精品视频| 国产成人精品久久二区二区91| 久久人人爽人人片av| 观看av在线不卡| 男女午夜视频在线观看| 欧美国产精品一级二级三级| 成年人免费黄色播放视频| 亚洲国产欧美网| 最近最新中文字幕大全免费视频 | 视频区图区小说| 丝袜喷水一区| 国产片内射在线| 国产欧美亚洲国产| 欧美人与性动交α欧美软件| 久久免费观看电影| 少妇粗大呻吟视频| 成人国产av品久久久| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 亚洲成色77777| 一区二区三区四区激情视频| 在线天堂中文资源库| 99九九在线精品视频| 亚洲一码二码三码区别大吗| 超碰成人久久| 男女免费视频国产| 欧美成人午夜精品| 菩萨蛮人人尽说江南好唐韦庄| 国产淫语在线视频| 91老司机精品| 亚洲九九香蕉| 一级黄色大片毛片| 午夜老司机福利片| 久久久精品免费免费高清| 国产99久久九九免费精品| 免费观看人在逋| 日韩一区二区三区影片| 菩萨蛮人人尽说江南好唐韦庄| av欧美777| 国产成人av激情在线播放| 精品国产一区二区三区久久久樱花| 亚洲精品一卡2卡三卡4卡5卡 | 自拍欧美九色日韩亚洲蝌蚪91| 尾随美女入室| 一区二区三区乱码不卡18| 中国美女看黄片| 亚洲 国产 在线| 美女扒开内裤让男人捅视频| 欧美日本中文国产一区发布| 在现免费观看毛片| 波多野结衣一区麻豆| 亚洲国产毛片av蜜桃av| 青青草视频在线视频观看| 99国产综合亚洲精品| 侵犯人妻中文字幕一二三四区| 国产国语露脸激情在线看| 成年动漫av网址| 午夜免费鲁丝| 久久人妻熟女aⅴ| 激情视频va一区二区三区| 精品久久久久久电影网| 国产精品久久久人人做人人爽| 91精品伊人久久大香线蕉| 亚洲av欧美aⅴ国产| 亚洲精品中文字幕在线视频| 91麻豆av在线| 精品少妇黑人巨大在线播放| 免费日韩欧美在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区有黄有色的免费视频| 9色porny在线观看| 免费女性裸体啪啪无遮挡网站| 韩国高清视频一区二区三区| 国产精品亚洲av一区麻豆| 久久av网站| 一级毛片我不卡| 国产黄色免费在线视频| 乱人伦中国视频| 国产人伦9x9x在线观看| 精品国产一区二区三区久久久樱花| 99热全是精品| 成人国语在线视频| 欧美黄色片欧美黄色片| 欧美 亚洲 国产 日韩一| 久久精品亚洲av国产电影网| av在线app专区| 亚洲国产精品一区三区| 免费观看人在逋| 国产精品久久久久久精品古装| 十分钟在线观看高清视频www| 亚洲av美国av| 日本欧美国产在线视频| 最黄视频免费看| 一区二区日韩欧美中文字幕| 在线观看免费高清a一片| 日韩一本色道免费dvd| 亚洲伊人色综图| 在线观看国产h片| 欧美激情高清一区二区三区| 婷婷色av中文字幕| av又黄又爽大尺度在线免费看| 老汉色av国产亚洲站长工具| 亚洲国产精品成人久久小说| 人人妻人人澡人人看| 午夜免费男女啪啪视频观看| 美国免费a级毛片| 国产一区二区三区av在线| 国产真人三级小视频在线观看| 高清不卡的av网站| 男女边摸边吃奶| 亚洲伊人色综图| 亚洲情色 制服丝袜| 大话2 男鬼变身卡| 国产爽快片一区二区三区| 黑人猛操日本美女一级片| 黑丝袜美女国产一区| 精品国产超薄肉色丝袜足j| 亚洲精品乱久久久久久| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合一区二区三区| 99热网站在线观看| 欧美黑人欧美精品刺激| 欧美精品啪啪一区二区三区 | 国产成人啪精品午夜网站| 亚洲中文av在线| 少妇 在线观看| 亚洲成国产人片在线观看| 日韩,欧美,国产一区二区三区| 日韩 亚洲 欧美在线| 女人精品久久久久毛片| 国产在线视频一区二区| 亚洲欧美一区二区三区国产| 亚洲伊人久久精品综合| 亚洲,欧美精品.| 久久久久网色| 极品人妻少妇av视频| 黄片小视频在线播放| 五月天丁香电影| 亚洲国产av新网站| 高清视频免费观看一区二区| 老司机影院成人| 亚洲国产av影院在线观看| 男女床上黄色一级片免费看| a 毛片基地| 女人精品久久久久毛片| 看免费成人av毛片| 秋霞在线观看毛片| 亚洲国产欧美在线一区| 晚上一个人看的免费电影| 久久久精品94久久精品| 久久人妻福利社区极品人妻图片 | 国产男人的电影天堂91| 90打野战视频偷拍视频| 免费看av在线观看网站| 夫妻午夜视频| 免费观看av网站的网址| 久久人人97超碰香蕉20202|