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

    融合無人機光譜信息與紋理信息的冬小麥生物量估測

    2018-09-03 12:42:44劉暢楊貴軍李振海湯伏全王建雯張春蘭張麗妍
    中國農(nóng)業(yè)科學(xué) 2018年16期
    關(guān)鍵詞:植被指數(shù)冬小麥紋理

    劉暢,楊貴軍,李振海,湯伏全,王建雯,張春蘭,張麗妍

    ?

    融合無人機光譜信息與紋理信息的冬小麥生物量估測

    劉暢1,2,3,4,楊貴軍2,3,4,李振海2,3,4,湯伏全1,王建雯2,3,4,張春蘭1,2,3,4,張麗妍2,3,4

    (1西安科技大學(xué)測繪科學(xué)與技術(shù)學(xué)院,西安 710054;2國家農(nóng)業(yè)信息化工程技術(shù)研究中心,北京 100097;3農(nóng)業(yè)部農(nóng)業(yè)遙感機理與定量遙感 重點實驗室,北京 100097;4北京市農(nóng)業(yè)物聯(lián)網(wǎng)工程技術(shù)研究中心,北京 100097)

    【目的】生物量是表征植被生命活動的重要參數(shù),對植被長勢監(jiān)測、產(chǎn)量預(yù)測有重要意義。以無人機為平臺的高光譜遙感技術(shù),具有機動靈活、成本低、空間覆蓋廣的優(yōu)勢,能夠及時準確地估測植被生物量,已成為遙感估算研究的熱點之一。由于光譜特征反演生物量存在飽和問題,因此,本研究嘗試結(jié)合紋理特征與植被指數(shù)構(gòu)建一種“圖-譜”融合指標,探究“圖-譜”融合指標的抗飽和能力及生物量估測能力。【方法】首先,利用無人機高光譜影像,提取其光譜信息和紋理信息,分別基于植被指數(shù)和紋理特征構(gòu)建生物量模型。其次,針對光譜特征存在的飽和問題,將植被指數(shù)與對生物量敏感的紋理指標相乘或相除兩種形式構(gòu)建“圖-譜”融合指標,分析“圖-譜”融合指標的飽和性,并基于“圖-譜”融合指標構(gòu)建生物量估算模型。最后,對比不同指標構(gòu)建的生物量模型的估測效果,來分析“圖-譜”融合指標估測生物量的能力。【結(jié)果】(1)植被指數(shù)多在LAI=5時出現(xiàn)飽和現(xiàn)象,而“圖-譜”融合指標VI×658, VI/658, VI/658, VI/658, VI/514, VI/514, VI/514, VI×802, VI×802均在LAI>5時才出現(xiàn)飽和現(xiàn)象,相比之下,這些“圖-譜”融合指標一定程度上改善了飽和問題;(2)與植被指數(shù)相比(除了GNDVI、NDVI之外),抗飽和能力提高的“圖-譜”融合指標VI×658、VI/658、VI/658、VI/658、VI/514、VI/514、VI/514、VI×802、VI×802,其與生物量的相關(guān)性也相對提高,所構(gòu)建的生物量模型精度較高(R=0.81,=826.02 kg·hm-2)。(3)對比單一植被指數(shù)、紋理特征,將紋理特征與光譜特征相結(jié)合的“圖-譜”融合指標估算小麥生物量的能力相對最強,模型精度明顯高于單一植被指數(shù)(R=0.69)和單一紋理特征(R=0.71)構(gòu)建的生物量模型?!窘Y(jié)論】“圖-譜”融合指標的抗飽和能力明顯提高,其構(gòu)建的生物量模型精度也有效提高,實現(xiàn)了結(jié)合光譜信息和紋理信息的冬小麥生物量遙感估測,為生物量定量反演提供一種新思路。

    生物量;“圖-譜”融合指標;紋理特征;飽和性;冬小麥

    0 引言

    【研究意義】生物量是生態(tài)系統(tǒng)中重要的理化參數(shù)之一[1],是評估植被生命活動的一個重要指標[2],并且對長勢監(jiān)測和作物估產(chǎn)有重要意義[3]。傳統(tǒng)的生物量估算方法不僅費時費力、破壞性大,且難以實現(xiàn)生物量的大范圍監(jiān)測[4]。近年來,隨著遙感技術(shù)的發(fā)展,高光譜遙感技術(shù)以其準確、快速、大尺度、非破壞性等優(yōu)勢已廣泛應(yīng)用于植被生物量的遙感估測中[2]?!厩叭搜芯窟M展】已有許多研究者[5-7]基于植被光譜特征,利用植被指數(shù)及其衍生光譜指數(shù)等對植株生物量進行了估測[8],并取得了一定的進展。SHIBAYAMA等[9]利用差值和比值植被指數(shù)構(gòu)建了水稻生物量的估算模型,精度R高達0.93。侯學(xué)會等[10]基于多種植被指數(shù)構(gòu)建了生物量估算模型,研究發(fā)現(xiàn),基于紅邊位置的指數(shù)模型估算冬小麥生物量的精度最高,決定系數(shù)R=0.83。劉瓊閣等[11]以TM遙感影像為數(shù)據(jù)源,基于波段及波段組合構(gòu)建了森林生物量的偏最小二乘回歸模型,估算精度較高,決定系數(shù)為0.72。然而,僅利用光譜特征進行生物量估測,在高植被密度條件下,植被指數(shù)會表現(xiàn)出對冠層生物量變化不敏感,即出現(xiàn)所謂的“飽和”現(xiàn)象[12],在識別地物及理化參數(shù)反演中容易出現(xiàn)偏差,影響反演的精度和效果。除光譜信息外,紋理信息作為另一重要遙感信息,是物體的一種固有屬性,不容易受到外界環(huán)境的影響,反映的是圖像灰度性質(zhì)及其空間關(guān)系,能夠提供地物的幾何信息,擴大了基于原始影像亮度的空間信息辨識度,能夠改善單一光譜信息源反演時存在的飽和問題,一定程度上提升理化參數(shù)反演的潛力[13]。對于紋理特征在定量反演中的應(yīng)用,許多學(xué)者進行了相關(guān)研究。GU等[14]對于植被覆蓋度進行了研究,發(fā)現(xiàn)植被指數(shù)和紋理特征的結(jié)合可以有效的反演植被覆蓋度。SARKER等[15]對森林生物量進行了研究,發(fā)現(xiàn)植被指數(shù)和紋理特征共同反演生物量的效果較好。曹慶先等[16]基于TM影像的紋理和光譜特征構(gòu)建的生物量多元回歸模型,可以高效的進行紅樹林生物量的估算,為紅樹林濕地生物量的估算提供了一種新方法。牧其爾等[17]以RapidEye多光譜影像為數(shù)據(jù)源,選取4種典型植被指數(shù),并提取了不同植被指數(shù)的灰度共生矩陣紋理特征值,構(gòu)建了植被指數(shù)、紋理特征各自與生物量的多元回歸模型,研究發(fā)現(xiàn),紋理特征估測生物量的能力高于單一的植被指數(shù)。【本研究切入點】目前,研究者多基于紋理特征與光譜特征的多元回歸模型進行生物量估算研究,并發(fā)現(xiàn)結(jié)合紋理與光譜信息進行生物量反演的效果更好。然而,從飽和問題出發(fā),聯(lián)系飽和現(xiàn)象與估算生物量的能力之間的關(guān)系,將光譜信息與紋理信息有機結(jié)合進行生物量反演的研究較為少見?!緮M解決的關(guān)鍵問題】為進一步探究植被指數(shù)的飽和問題,本文嘗試將對生物量敏感的紋理特征加入到植被指數(shù)當中,考慮各紋理指標隨著葉面積指數(shù)(LAI)增大的變化情況不同,采用將植被指數(shù)與紋理指標相乘或相除兩種形式構(gòu)建一種“圖-譜”融合指標,探究“圖-譜”融合指標的抗飽和能力及估測生物量的能力,以實現(xiàn)結(jié)合光譜信息和紋理信息的冬小麥生物量遙感估測,同時發(fā)揮光譜信息與紋理信息的優(yōu)勢,從而更好的反演植被生物量,為生物量定量反演提供一種新思路。

    1 材料與方法

    1.1 試驗區(qū)概況

    試驗于2014—2015年冬小麥生長季,在北京市昌平區(qū)小湯山國家精準農(nóng)業(yè)研究示范基地(40°00′—40°21′N,116°34′—117°00′E,海拔高度36 m)進行。試驗區(qū)氣候溫和,四季分明,屬于典型的暖溫帶半濕潤大陸季風(fēng)氣候。前茬作物為玉米,試驗土壤類型為潮土,土壤中0—30 cm土層中硝態(tài)氮含量為3.16—14.82 mg·kg-1,全氮含量為1.0—1.2 g·kg-1,有機質(zhì)含量為15.8—20.0 g·kg-1,有效磷含量為3.14—21.18 mg·kg-1,速效鉀含量為86.83—120.62 mg·kg-1。

    試驗設(shè)計為品種、氮肥、水分的正交試驗(圖1)。冬小麥品種分別為京9843(J9843)和中麥175(ZM175);氮肥(尿素)施用量包括0(N1)、195 kg·hm-2(N2)、390 kg·hm-2(N3)和585 kg·hm-2(N4)4個水平;水分包括雨養(yǎng)(W1,0)、正常灌水(W2,192 mm)和2倍正常灌水(W3,384 mm)3個水平。試驗共有48個小區(qū),16個小區(qū)為一組,重復(fù)3次。

    圖1 研究區(qū)概況及試驗小區(qū)分布

    1.2 數(shù)據(jù)獲取

    1.2.1 地面數(shù)據(jù)獲取 試驗于冬小麥挑旗期(2015年4月26日)進行,試驗數(shù)據(jù)主要為地上部生物量數(shù)據(jù)和葉面積指數(shù)。

    在室外調(diào)查固定樣方2行×1m(行距15 cm),進行群體密度調(diào)查,并在各小區(qū)選取具有代表性的20個小麥植株作為樣本。獲取的樣本分器官處理,測定鮮重后放入紙袋中,置于烘箱105℃殺青30 min,然后設(shè)置溫度到75℃烘干至恒重(約24—48 h),測定各器官干重,最終通過群體密度和樣本干重,計算得到單位面積上的冬小麥生物量。

    總生物量=(葉干重+莖干重+穗干重)/取樣莖數(shù)×畝莖數(shù) (1)

    在試驗小區(qū)樣方內(nèi)選取20個具有代表性的小麥作為樣本,將獲取的樣本莖葉分離,選取樣本中具有代表性的葉片,采用手持式葉面積激光儀CI-203(CID Inc.,USA)測定其葉片面積,求和計算出樣本葉片總面積,最后通過測定的群體密度,計算出冬小麥的葉面積指數(shù)。

    葉面積指數(shù)(LAI)= 樣方葉片總面積/樣方面積 (2)

    1.2.2 高光譜數(shù)據(jù)獲取 本研究以八旋翼電動無人機為近地飛行平臺搭載UHD185成像光譜儀(Cubert,德國)(圖2),選擇晴朗無云天氣,在中午12:00—13:00對小麥試驗田進行航拍,飛行高度為50 m。獲取空間分辨率為0.012 m的小湯山冬小麥挑旗期(2015年4月26日)一景高光譜影像。

    為消除成像過程中傳感器性能、大氣傳輸?shù)纫蛩氐挠绊?,集成高光譜分辨率和高空間分辨率為一體的遙感影像,本文對原始影像進行數(shù)據(jù)預(yù)處理。主要包括:利用無人機自帶Cubert-Pilot軟件(Cubert,德國)進行影像融合,利用Agisoft PhotoScan軟件(Agisoft LLC,俄羅斯)對融合后影像進行圖像拼接[18],利用ENVI軟件(Exelis Visual Information Solutions,美國)進行輻射校正,大氣校正、裁剪[19]。

    圖2 搭載高光譜傳感器(UHD185)的無人機平臺

    1.3 高光譜影像光譜分析

    根據(jù)前人研究及相關(guān)參考文獻,選取了10個與生物量相關(guān)的植被指數(shù)對冬小麥生物量進行估算(表1)。

    表1 本文用到的植被指數(shù)

    公式中i代表波段i處的反射率i denotes reflectance at band i

    1.4 高光譜影像紋理分析

    1.4.1 波段篩選 通過原始波段反射率與生物量的相關(guān)性分析,篩選出3個與生物量相關(guān)性較大的波段,并將其作為紋理特征分析波段。

    1.4.2 紋理特征提取 灰度共生矩陣法[20]是1973 年由HARALICK提出的目前應(yīng)用最廣的紋理提取方法之一,因其具有旋轉(zhuǎn)不變性和多尺度特性,且計算復(fù)雜度小,受到了較多的關(guān)注[21]。在ENVI 5.3中,利用灰度共生矩陣方法(GLCM)對所篩選的3個波段,進行0°、45°、90°、135° 4個方向的8個紋理特征的提取(表2),對不同方向的紋理特征進行平均,得到各波段的8個紋理特征。然后在ArcGIS 10.3.1(Esri,美國)中,對于各波段的紋理特征影像分別進行感興趣區(qū)域劃定,提取出所劃區(qū)域的紋理值,并將該值作為這個小區(qū)的紋理特征值。

    1.4.3 紋理特征與生物量的相關(guān)性分析 利用SPSS 22.0軟件(IBM,USA),采用皮爾遜相關(guān)分析法分析紋理特征與冬小麥地上生物量之間的相關(guān)性,篩選出與生物量相關(guān)性較大的紋理特征構(gòu)建生物量估算模型。采用多元逐步回歸法,基于所篩選紋理特征構(gòu)建生物量估算模型。

    1.5 飽和性分析及“圖-譜”融合指標的建立

    植被指數(shù)在植被覆蓋度較大時易過早飽和,嘗試將與生物量相關(guān)的紋理特征加入到植被指數(shù)當中,考慮各紋理指標隨著LAI增大的變化情況不同,采用將植被指數(shù)與紋理指標相乘或相除兩種形式,構(gòu)建“圖-譜”融合指標(表3),并通過飽和點位置的變化來探究“圖-譜”融合指標的抗飽和能力。

    表2 研究用到的紋理特征及其計算公式

    表3 本研究用到的“圖-譜”融合指標

    表中658代表波段658 nm的紋理特征;其他指標代表的含義類似。VI分別代表植被指數(shù)RVI、MTVI2、MSAVI、EVI、DVI、SAVI、RDVI、OSAVI、NDVI、GNDVI,將它們分別與紋理特征相乘或相除,構(gòu)成“圖-譜”融合指標。例如,表中VI/658分別代表“圖-譜”融合指標RVI/658,MTVI2/658,MSAVI/658,EVI/658,DVI/658,SAVI/658,RDVI/658,OSAVI/658,NDVI/658,GNDVI/658,其他指標代表的含義類似

    658denotes textureat band 658 nm, and the same as the others index; VI denotes vegetation index of RVI, MTVI2, MSAVI, EVI, DVI, SAVI, RDVI, OSAVI, NDVI, GNDVI, and establishing 'image and spectrum' fusion index by multiplying or dividing each of vegetation index with the texture features. For example, VI/658denotes RVI/658, MTVI2/658, MSAVI/658, EVI/658, DVI/658, SAVI/658, RDVI/658, OSAVI/658, NDVI/658, GNDVI/658, and the same as the other index

    1.6 模型構(gòu)建及數(shù)據(jù)分析

    本研究中,選取挑旗期2個重復(fù)試驗的數(shù)據(jù)進行模型構(gòu)建(樣本量n=32),該生育期的另一重復(fù)試驗數(shù)據(jù)進行模型驗證(n=16)。最后根據(jù)建模精度(R作為評判標準)和驗證精度(作為評判標準)來評判模型的擬合效果。一般來說,R越高,越低,建模與驗證精度越高,依此選出最佳生物量模型。

    2 結(jié)果

    2.1 基于光譜指數(shù)的生物量估算

    2.1.1 植被指數(shù)與生物量的相關(guān)性分析 對選取的10個植被指數(shù)和生物量的相關(guān)性進行分析(表4)。由表中可以得到,植被指數(shù)與生物量之間的相關(guān)性均達到極顯著水平(<0.01),其中與生物量相關(guān)性最大的是GNDVI,相關(guān)系數(shù)為0.776,最小是DVI,相關(guān)系數(shù)為0.547;相關(guān)性由大到小依次為:GNDVI、NDVI、RVI、OSAVI、MTVI2、MSAVI、RDVI、SAVI、EVI、DVI。

    表4 植被指數(shù)和生物量的相關(guān)性分析

    ** 表示相關(guān)性在0.01水平下達到顯著。下同

    ** means the correlation is significant at 0.01 level. The same as below

    2.1.2 基于植被指數(shù)的生物量估算模型 采用多元逐步回歸法,基于10個植被指數(shù)建立生物量回歸模型(圖3),同時基于植被指數(shù)的回歸模型得到冬小麥生物量的結(jié)果分布(圖4)。圖3中可以看出,生物量模型精度為R=0.69,=874.25 kg·hm-2,擬合效果較好,表明利用植被指數(shù)估算冬小麥生物量是可行的。同時發(fā)現(xiàn),在生物量大于7 000 kg·hm-2時,部分預(yù)測值略低于實測值,可能是植被初期生長的過程中,隨著葉片逐漸增大,出現(xiàn)了輕微的飽和現(xiàn)象。由圖4中可以看出,不施氮處理(N1)和0.5倍常規(guī)施氮處理(N2)的冬小麥生物量較小,多分布在5 100 kg·hm-2以下,部分小區(qū)的冬小麥生物量分布在5 100—6 300 kg·hm-2,分布范圍較大;而常規(guī)施氮處理(N3)和1.5倍常規(guī)施氮處理(N4)的冬小麥生物量相對較大,多分布在5 100—7 500 kg·hm-2。雨養(yǎng)處理下(W1),冬小麥生物量分布在5 100 kg·hm-2以下;在正常水(W2)及2倍正常水(W3)處理下,冬小麥的生物量較大,生物量多集中在5 100—7 500 kg·hm-2;但是W2與W3處理之間的生物量差異并不很明顯。能夠看出,冬小麥生物量隨著施氮水平的提高,以及水分處理水平的提升,整體呈逐步增加的趨勢,但是,生物量的分布范圍較大,具備的規(guī)律性較差,可能是植被指數(shù)進行生物量反演的過程中受到了飽和現(xiàn)象的影響,因此,僅利用光譜特征進行生物量反演并不能十分準確的反映出實際生物量的分布情況,存在一定的局限。

    圖3 多個植被指數(shù)的生物量模型預(yù)測值與實際值對比

    2.2 基于紋理特征的生物量估算

    2.2.1 紋理特征波段篩選 原始波段反射率與生物量的相關(guān)性分析如圖5所示,由圖中看出,454—726 nm波段光譜反射率與生物量呈負相關(guān),730—882 nm波段區(qū)間呈正相關(guān),其中分別在514 nm(= -0.697),658 nm(= -0.775)處形成兩個波谷,在802 nm(=0.565)波段處形成高反射平臺,因此最終選擇與生物量相關(guān)性較大的3個波段,綠波段(514 nm)、紅波段(658 nm)以及近紅外波段 (802 nm)作為紋理特征分析波段。

    2.2.2 基于紋理特征的生物量估算模型 經(jīng)過對比不同窗口的紋理特征后發(fā)現(xiàn),本文所采用的高光譜影像在3×3窗口下的分辨率最適合進行紋理特征分析,因此,選取3×3窗口下,與生物量相關(guān)性較好的紋理特征,即514(514 nm波段處的特征)、514、514、514、514、658、658、658、658、658、658、802、802、802這些紋理特征作為生物量的估測指標,對于冬小麥生物量的估算有積極意義。

    圖4 基于植被指數(shù)的冬小麥生物量反演結(jié)果

    圖5 光譜反射率與生物量的相關(guān)性曲線圖

    采用多元逐步回歸法,考慮前文所篩選的全部紋理特征(514、514、514、514、514、658、658、658、658、658、658、802、802、802),構(gòu)建生物量回歸模型(圖6),同時基于紋理特征的回歸模型得到冬小麥生物量的結(jié)果分布(圖7)。圖6中可以看出,生物量模型精度為R=0.71,=828.87 kg·hm-2,與光譜特征構(gòu)建的生物量模型(R=0.69,=874.25 kg·hm-2)相比,估算模型精度R有所提高,均方根誤差減小,擬合效果較好,可能原因是紋理特征的模型中含有豐富的紋理信息,一定程度上彌補了光譜信息的不足,緩解了光譜特征反演時存在的飽和以及“同物異譜,同譜異物”現(xiàn)象,從而提高了模型的擬合精度。由圖7中可以看出,不施氮處理(N1)和0.5倍常規(guī)施氮處理(N2)的冬小麥生物量較小,多分布在3 900—6 300 kg·hm-2,部分小區(qū)的冬小麥生物量分布在3 900 kg·hm-2以下;而常規(guī)施氮處理(N3)和1.5倍常規(guī)施氮處理(N4)的冬小麥生物量相對較大,多分布在5 100—7 500 kg·hm-2。雨養(yǎng)處理下(W1),冬小麥生物量分布在5 100 kg·hm-2以下;在正常水(W2)及2倍正常水(W3)處理下,冬小麥的生物量較大,生物量多集中在5 100— 7 500 kg·hm-2,相比之下,W3處理比W2處理下的生物量更大,冬小麥長勢更好一些。整體上,不同處理之間的生物量差異及規(guī)律能夠較好的呈現(xiàn)出來,相比基于植被指數(shù)的冬小麥生物量反演結(jié)果要更加具有規(guī)律性,更加貼近實際情況。因此,基于紋理特征的生物量模型的擬合效果較好,能較準確地反演冬小麥地上生物量,具有一定的應(yīng)用價值。

    2.3 基于“圖-譜”融合指標的生物量估算

    2.3.1“圖-譜”融合指標的飽和性分析 在利用植被指數(shù)估測生物量時,發(fā)現(xiàn)當生物量超過一定量時,出現(xiàn)部分預(yù)測值低于實測值的現(xiàn)象。主要原因可能是隨著植被生長,葉片逐漸增大,植被覆蓋度逐漸增加,使得植被指數(shù)對于生物量變化不再敏感,造成了光譜的飽和現(xiàn)象。因此,可通過分析光譜指標隨LAI增大的變化情況來分析各光譜指標的飽和性。不同光譜指數(shù)隨LAI的變化情況如圖所示(圖8—圖11)。本研究對14種“圖-譜”融合指標都進行了飽和性探究,圖8—圖11僅列出了部分“圖-譜”融合指標(VI,VI×658,VI×802,VI×802)隨著LAI 的變化情況,但是分析時仍對其他的“圖-譜”融合指標飽和性進行分析,最終的結(jié)果分析包括所有的“圖-譜”融合指標。

    對比植被指數(shù)與“圖-譜”融合指標的飽和點位置可以發(fā)現(xiàn)(圖8—圖11),大部分植被指數(shù)在LAI=5的時候出現(xiàn)飽和現(xiàn)象;“圖-譜”融合指標VI×658、VI/658、VI/658、VI/658、VI/514、VI/514、VI/514、VI×802、VI×802隨LAI增大幾乎呈線性趨勢遞增,或是先逐步增大后趨于穩(wěn)定,均在LAI>5后出現(xiàn)飽和現(xiàn)象,延后了飽和點位置,其抗飽和能力有所提高;但是“圖-譜”融合指標VI/514、VI×514、VI/658、VI×658、VI×802隨著LAI的增大,先逐步增大后略微下降,LAI在5以內(nèi)或大約等于5 時出現(xiàn)了飽和點,其飽和現(xiàn)象并沒有得到改善。

    2.3.2“圖-譜”融合指標與生物量的相關(guān)性分析 將植被指數(shù)、“圖-譜”融合指標分別與生物量進行相關(guān)性分析(表5)。結(jié)果表明,所有指標與生物量的相關(guān)性都達到了極顯著水平(<0.01),與生物量相關(guān)性相對較大的2種“圖-譜”融合指標為VI×802和VI×802,其中,相關(guān)性最大的“圖-譜”融合指標是RVI×802(=0.815)。對比植被指數(shù)和生物量之間的相關(guān)性,VI×658,VI/658,VI×658,VI/658,VI/658,VI/514,VI×514,VI/514,VI/514,VI×802,VI×802中大部分指標與生物量的相關(guān)性都有所提高;而VI/658,VI/514,VI×802中大多指標的相關(guān)性沒有提高,反而有所下降。

    圖6 多個紋理特征的生物量模型預(yù)測值與實際值對比

    圖7 基于紋理特征的冬小麥生物量反演結(jié)果

    圖8 各植被指數(shù)隨LAI的變化情況

    圖9 “圖-譜”融合指標VI×sm658隨LAI的變化情況

    圖10 “圖-譜”融合指標VI×con802隨LAI的變化情況

    圖11 “圖-譜”融合指標VI×mean802隨LAI的變化情況

    因此,最終篩選抗飽和能力強,且與生物量相關(guān)性相對較大的“圖-譜”融合指標作為生物量的估測指標,即VI×658,VI/658, VI/658, VI/658,VI/514, VI/514,VI/514,VI×802,VI×802。

    表5 植被指數(shù)和“圖-譜”融合指標與生物量的相關(guān)性分析

    2.3.3 基于“圖-譜”融合指標的生物量估算模型 采用多元逐步回歸法,基于所篩選的“圖-譜”融合指標構(gòu)建生物量估算模型(圖12),同時基于“圖-譜”融合指標的回歸模型得到冬小麥生物量的結(jié)果分布(圖13)。圖12可以看出,生物量模型精度為R=0.81,=826.02 kg·hm-2,與紋理特征構(gòu)建的生物量估算模型(R=0.71,=828.87 kg·hm-2)和植被指數(shù)構(gòu)建的生物量估算模型(R=0.69,=874.25 kg·hm-2)相比,“圖-譜”融合指標構(gòu)建的生物量估算模型精度R明顯提高,均方根誤差較小,擬合效果最好。由圖13中可以看出,不施氮處理(N1)和0.5倍常規(guī)施氮處理(N2)的冬小麥生物量較小,多分布在2 000—5 900 kg·hm-2;而常規(guī)施氮處理(N3)和1.5倍常規(guī)施氮處理(N4)的冬小麥生物量相對較大,多分布在4 600—7 200 kg·hm-2。雨養(yǎng)處理下(W1),冬小麥生物量分布在5 900 kg·hm-2以下;在正常水(W2)及2倍正常水(W3)處理下,冬小麥的生物量較大,生物量多集中在4 600—7 200 kg·hm-2。不同水氮處理間的冬小麥生物量值分布存在明顯的階梯性,而每個處理水平內(nèi)的冬小麥生物量分布較集中,能夠很好的表征冬小麥生物量。

    結(jié)果表明,融合光譜信息與紋理信息的“圖-譜”融合指標,綜合考慮了紋理特征及光譜特征對生物量的貢獻性。敏感的植被指數(shù)對于不同生物量下的光譜特征有很好的表征,而豐富的紋理信息減弱了光譜特征存在的飽和問題,抗飽和能力得到明顯改善,因此,其生物量回歸模型估算精度提高,相比單一紋理特征、單一光譜特征的生物量模型,紋理與光譜共同反演冬小麥地上生物量的效果最好。

    圖12 多個“圖-譜”融合指標的生物量模型預(yù)測值與實際值對比

    3 討論

    前人多利用光譜信息對生物量進行估算,取得了一定的進展,但是基于光譜信息進行生物量反演,在植被覆蓋度較大時易出現(xiàn)飽和問題。紋理特征應(yīng)用于分類研究中,研究者[22-24]發(fā)現(xiàn)加入紋理特征后分類精度有了明顯的提升;而在森林結(jié)構(gòu)參數(shù)反演的研究中[25-27],紋理特征能夠反映樹木的冠層結(jié)構(gòu)及異質(zhì)性,引入紋理特征估算結(jié)構(gòu)參數(shù),能夠取得較好的估算效果?;诩y理特征進行生物量反演的研究,李明詩等[28]發(fā)現(xiàn),少數(shù)紋理特征ME(均值)、VA(方差)在森林生物量估算中是有效的。前人多是直接選取對生物量敏感的紋理及光譜特征,同時基于所選紋理和光譜特征構(gòu)建多元回歸模型進行生物量估算,而本研究從飽和問題出發(fā),利用無人機獲取的高光譜影像,提取其光譜特征與紋理特征,主要探究了光譜與紋理信息相結(jié)合的“圖-譜”融合指標的抗飽和能力和估測生物量的能力。研究發(fā)現(xiàn),相比植被指數(shù),本文所構(gòu)建的“圖-譜”融合指標VI×658,VI/658,VI/658,VI/658,VI/514,VI/514,VI/514,VI×802,VI×802延后了對LAI響應(yīng)的飽和點位置,其抗飽和能力有所提升,同時,對比植被指數(shù)與生物量的相關(guān)性發(fā)現(xiàn),這些指標與生物量的相關(guān)性也有所提高,因此選取具有抗飽和特性并且對生物量敏感的指標,探究其反演生物量的效果。最終結(jié)果表明,基于“圖-譜”融合指標構(gòu)建的生物量估算模型精度明顯提高,主要是因為“圖-譜”融合指標考慮了光譜特征與紋理特征兩者共同對生物量的貢獻性,敏感的植被指數(shù)對于不同生物量下的光譜特征有了很好的表征,而豐富的紋理信息提供了作物的空間特征,減弱了單一植被指數(shù)反演時存在的飽和現(xiàn)象,一定程度上彌補了光譜信息的不足,“圖-譜”融合指標結(jié)合了兩者的優(yōu)勢,提高了模型的擬合精度。本文的研究結(jié)果發(fā)現(xiàn),相比單一光譜特征、單一紋理特征,紋理與光譜共同反演冬小麥地上生物量的效果最佳,與前人[16-17, 28]的相關(guān)研究結(jié)果一致。然而,前人基于紋理特征的生物量反演,主要應(yīng)用在紅樹林[14]、綠洲沙漠植被[15]、森林[29-30]等方面,對于冬小麥的生物量反演研究較為少見。同時,本研究對于飽和問題具有一定的針對性。本文嘗試從提升光譜特征的抗飽和能力出發(fā),聯(lián)系抗飽和能力與估算生物量能力之間的關(guān)系,探究了不同指標的生物量估算能力,為冬小麥生物量反演提供了一種新思路。同時,拓寬了紋理特征在長勢監(jiān)測、產(chǎn)量估測等方面的應(yīng)用前景。

    圖13 基于“圖-譜”融合指標的冬小麥生物量反演結(jié)果

    本研究基于“圖-譜”融合指標構(gòu)建的生物量模型精度提高,但是均方根誤差并沒有明顯改善,可能的原因是,研究僅利用了冬小麥挑旗期的高光譜影像及地面實測數(shù)據(jù),所使用的樣本量較少,對所構(gòu)建模型的通用性可能會有一定影響,未來的相關(guān)研究應(yīng)增加試驗生育期、品種,拓展模型應(yīng)用范圍,從而提高結(jié)果的普適性和穩(wěn)定性。另外,目前存在的紋理特征提取方法多是針對單波段影像,不適用于包含多波段的高光譜影像,而本文是通過原始波段反射率與生物量的相關(guān)性分析,篩選得到紋理特征分析波段,考慮因素不夠全面,對于最終估算模型的構(gòu)建可能會產(chǎn)生一定的影響,因此,還需要進一步的研究,探究出一種更適用高光譜影像波段篩選的方法,從而更加合理有效地利用高光譜影像的波段信息。

    4 結(jié)論

    本文基于無人機高光譜影像,對光譜信息和紋理信息估測冬小麥地上部生物量的能力及其飽和性進行了分析,主要獲得以下結(jié)論:

    (1)相比植被指數(shù),“圖-譜”融合指標VI×658,VI/658, VI/658, VI/658,VI/514, VI/514,VI/514,VI×802,VI×802延后了LAI的飽和點位置,其抗飽和能力明顯提高,而“圖-譜”融合指標VI/658,VI/514,VI×802,VI×658,VI×514的抗飽和性并沒有得到提升。

    (2)對比植被指數(shù)與生物量的相關(guān)性(除NDVI,GNDVI外)發(fā)現(xiàn),“圖-譜”融合指標VI×658,VI/658,VI×658,VI/658, VI/658,VI/514,VI×514,VI/514,VI/514,VI×802,VI×802中大部分指標與生物量的相關(guān)性都明顯提高,而VI/658,VI/514,VI×802中大多指標與生物量的相關(guān)性沒有提高,反而有所下降。

    (3)與植被指數(shù)相比,所有抗飽和能力提高的“圖-譜”融合指標,其與生物量的相關(guān)性也明顯增強;抗飽和性沒有得到改善的“圖-譜”融合指標(除VI×658,VI×514外),其與生物量的相關(guān)性并沒有增大,反而有所減小。同時,各指標與生物量的相關(guān)性大小和抗飽和能力的高低是成正比的。

    (4)相比單一光譜特征、紋理特征,將紋理特征與光譜特征相結(jié)合的“圖-譜”融合指標估算小麥生物量的能力最強,模型精度(R=0.81)明顯高于植被指數(shù)(R=0.69)、紋理特征(R=0.71)構(gòu)建的生物量模型,表明融合光譜與紋理信息共同反演生物量的效果較好,具有一定的優(yōu)勢。

    [1] 李佳佳. 小麥生物物理與生物化學(xué)參數(shù)的高光譜遙感監(jiān)測[D]. 南京: 南京信息工程大學(xué), 2015.

    LI J J. Monitoring bio-physical and bio-chemical parameters of wheat by hyper-spectral remote sensing[D]. Nanjing: Nanjing University of Information Science and Technology, 2015. (in Chinese)

    [2] 張凱, 王潤元, 王小平, 趙鴻, 韓海濤. 黃土高原春小麥地上鮮生物量高光譜遙感估算模型. 生態(tài)學(xué)雜志, 2009, 28(6): 1155-1161.

    ZHANG K, WANG R Y, WANG X P, ZHAO H, HAN H T. Hyperspectral remote sensing estimation models for aboveground fresh biomass of spring wheat on Loess Plateau., 2009, 28(6): 1155-1161. (in Chinese)

    [3] 尚艷. 不同氮水平下小麥冠層光譜特征及其與農(nóng)學(xué)參數(shù)關(guān)系研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2015.

    SHANG Y. Wheat canopy spectral features and its research relationship with agronomy parameter under different nitrogen levels[D]. Yangling: Northwest Agriculture and Forestry University, 2015. (in Chinese)

    [4] 范云豹, 宮兆寧, 趙文吉, 張敏. 基于高光譜遙感的植被生物量反演方法研究. 河北師范大學(xué)學(xué)報(自然科學(xué)版), 2016, 40(3): 267-271.

    FAN Y B, GONG Z N, ZHAO W J, ZHANG M. Study on vegetation biomass inversion method based on hyperspectral remote sensing., 2016, 40(3): 267-271. (in Chinese)

    [5] 譚昌偉, 楊昕, 羅明, 馬昌, 嚴翔, 陳亭亭. 以HJ-CCD影像為基礎(chǔ)的冬小麥孕穗期關(guān)鍵苗情參數(shù)遙感定量反演. 中國農(nóng)業(yè)科學(xué), 2015, 48(13): 2518-2527.

    TAN C W, YANG X, LUO M, MA C, YAN X, CHEN T T. Quantitative inversion of key seedling condition parameters in winter wheat at booting stage using remote sensing based on HJ-CCD images.

    , 2015, 48(13): 2518-2527. (in Chinese)

    [6] 陳鵬飛, 王卷樂, 廖秀英, 尹芳, 陳寶瑞, 劉睿. 基于環(huán)境減災(zāi)衛(wèi)星遙感數(shù)據(jù)的呼倫貝爾草地地上生物量反演研究. 自然資源學(xué)報, 2010, 25(7): 1122-1131.

    CHEN P F, WANG J L, LIAO X Y, YIN F, CHEN B R, LIU R. Using data of HJ-1A/B for hulunbeier grassland aboveground biomass estimation., 2010, 25(7): 1122-1131. (in Chinese)

    [7] 高明亮, 宮兆寧, 趙文吉, 高陽, 胡東. 基于植被指數(shù)的北京軍都山荊條灌叢生物量反演研究. 生態(tài)學(xué)報, 2014, 34(5): 1178-1188.

    GAO M L, GONG Z N, ZHAO W J, GAO Y, HU D. The study ofshrubs canopy biomass inversion in Beijing Jundu mountains area based on vegetation indices.2014, 34(5): 1178-1188. (in Chinese)

    [8] 趙天舸, 于瑞宏, 張志磊, 白雪松, 曾慶奧. 濕地植被地上生物量遙感估算方法研究進展. 生態(tài)學(xué)雜志, 2016, 35(7): 1936-1946.

    ZHAO T G, YU R H, ZHANG Z L, BAI X S, ZENG Q A. Estimation of wetland vegetation aboveground biomass based on remote sensing data: A review., 2016, 35(7): 1936-1946. (in Chinese)

    [9] SHIBAYAMA M, AKIYAMA T. Seasonal visible, near-infrared and mid-infrared spectra of rice canopies in relation to LAI and above-ground dry phytomass., 1989, 27(2): 119-127.

    [10] 侯學(xué)會, 牛錚, 黃妮, 許時光. 小麥生物量和真實葉面積指數(shù)的高光譜遙感估算模型. 國土資源遙感, 2012, 24(4): 30-35.

    HOU X H, NIU Z, HUANG N, XU S G. The hyperspectral remote sensing estimation models of total biomass and true LAI of wheat., 2012, 24(4): 30-35. (in Chinese)

    [11] 劉瓊閣, 彭道黎, 涂云燕, 李艷麗, 高東啟. 基于偏最小二乘的森林生物量遙感估測. 東北林業(yè)大學(xué)學(xué)報, 2014, 42(7): 44-47.

    LIU Q G, PENG D L, TU Y Y, LI Y Y, GAO D Q. Estimating forest biomass by partial least squares regression., 2014, 42(7): 44-47. (in Chinese)

    [12] 陳鵬飛, Nicolas T, 王紀華, PHILIPPE V, 黃文江, 李保國. 估測作物冠層生物量的新植被指數(shù)的研究. 光譜學(xué)與光譜分析, 2010, 30(2): 512-517.

    CHEN P F, NICOLAS T, WANG J H, PHILIPPE V, HUANG W J, LI B G. New index for crop canopy fresh biomass estimation., 2010, 30(2): 512-517. (in Chinese)

    [13] 劉俊, 畢華興, 朱沛林, 孫菁, 朱金兆, 陳濤. 基于 ALOS 遙感數(shù)據(jù)紋理及紋理指數(shù)的柞樹蓄積量估測. 農(nóng)業(yè)機械學(xué)報, 2014, 45(7): 245-254.

    LIU J, BI H X, ZHU P L, SUN J, ZHU J Z, CHEN T. Estimating stand volume offorest based on texture parameters and derivative texture indices of ALOS imagery., 2014, 45(7): 245-254. (in Chinese)

    [14] GU Z J, JU W M, LI L, LI D Q, LIU Y B, FAN W L. Using vegetation indices and texture measures to estimate vegetation fractional coverage (VFC) of planted and natural forests in Nanjing city, China., 2013, 51(7): 1186-1194.

    [15] SARKER L R, NICHOL J E. Improved forest biomass estimates using ALOS AVNIR-2 texture indices., 2011, 115(4): 968-977.

    [16] 曹慶先, 徐大平, 鞠洪波. 基于TM影像紋理與光譜特征的紅樹林生物量估算. 林業(yè)資源管理, 2010, 12(6): 102-107.

    CAO Q X, XU D P, JU H B. The biomass estimation of mangrove community based on the textural features and spectral information of TM images., 2010, 12(6): 102-107. (in Chinese)

    [17] 牧其爾, 高志海, 包玉海, 王琫瑜, 白黎娜. 植被指數(shù)紋理特征信息估測稀疏植被生物量. 遙感信息, 2016, 31(1): 58-63.

    MU Q E, GAO Z H, BAO Y H, WANG B Y, BAI L N. Estimation of sparse vegetation biomass based on grey-level co-occurrence matrix of vegetation indices., 2016, 31(1): 58-63. (in Chinese)

    [18] YUE J B, YANG G J, LI C C, LI Z H, WANG Y J, FENG H K, XU B. Estimation of winter wheat above-ground biomass using unmanned aerial vehicle-based snapshot hyperspectral sensor and crop height improved models., 2017, 9(70): 801-819.

    [19] 鄧書斌. ENVI遙感圖像處理方法. 北京: 科學(xué)出版社, 2010.

    DENG S B.Beijing: Science Press, 2010. (in Chinese)

    [20] HARALICK R M, SHANMUGAM K, DINSTEIN I. Textural features for image classification.1973, 3(6): 768-780.

    [21] 劉廣東. 基于TM影像植被指數(shù)和紋理特征的對比關(guān)系研究[D]. 重慶: 重慶師范大學(xué), 2010.

    LIU G D. Study on the contrast relationships between NDVI and texture features based on TM image [D]. Chongqing:Chongqing Normal University, 2010. (in Chinese)

    [22] FRANKLIN S E, HALL R J, MOSKAL L M. Incorporating texture into classification of forest species composition from airborne multispectral images., 2000, 21(1): 61-79.

    [23] TREITZ P, HOWARTH P. Integrating spectral, spatial, and terrain variables for forest ecosystem classification., 2000, 66(3): 305-318.

    [24] JOHANSEN K, COOPS N C, GERGEL S E. Application of high spatial resolution satellite imagery for riparian and forest ecosystem classification., 2007, 110(1): 29-44.

    [25] FRANKLIN S E, WULDER M A, GERYLO G R. Texture analysis of IKONOS panchromatic data for Douglas-fir forest age class separability in British Columbia., 2001, 22(13): 2627-2632.

    [26] KAYITAKIRE F, HAMEL C, DEFOURNY P. Retrieving forest structure variables based on image texture analysis and IKONOS-2 imagery., 2006, 102(3/4): 390-401.

    [27] OZDEMIR I, KARNIELI A. Predicting forest structural parameters using the image texture derived from WorldView-2 multispectral imagery in a dryland forest, Israel., 2011, 13(5): 701-710.

    [28] 李明詩, 譚瑩, 潘潔, 彭世揆. 結(jié)合光譜、紋理及地形特征的森林生物量建模研究. 遙感信息, 2006(6): 6-9.

    LI M S, TAN Y, PAN J, PENG S K. Modeling forest aboveground biomass by combining the spectrum, textures with topographic features., 2006, 2006(6): 6-9. (in Chinese)

    [29] SARKER M L, NICHOL J, IZ H B, AHMAD B B, RAHMAN A A. Forest biomass estimation using texture measurements of high- resolution dual-polarization C-band SAR data., 2013, 51(6): 3371-3384.

    [30] MIGUEL A C, MARTIN R, BERNARDUS H J. Estimation of tropical forest structure from SPOT-5 satellite images., 2010, 31(10): 2767-2782.

    (責(zé)任編輯 楊鑫浩)

    Biomass Estimation in Winter Wheat by UAV Spectral Information and Texture Information Fusion

    LIU Chang1,2,3,4, YANG GuiJun2,3,4, LI ZhenHai2,3,4, TANG FuQuan1, WANG JianWen2,3,4, ZHANG ChunLan1,2,3,4, ZHANG LiYan2,3,4

    (1College of Geomatics, Xi’an University of Science and Technology, Xi’an 710054;2National Engineering Research Center for Information Technology in Agriculture, Beijing 100097;3Key Laboratory of Quantitative Remote Sensing in Agriculture, Ministry of Agriculture, Beijing 100097;4Beijing Engineering Research Center for Agriculture Internet of Things, Beijing 100097)

    【Objective】Biomass, an important parameter to characterize vegetation activities, is of great significance for plant growth monitoring and yield forecasting. Hyperspectral remote sensing technology based on the unmanned aerial vehicle (UAV) has the advantages of flexibility, non-destructive and wide coverage, and could also timely and accurately estimate vegetation biomass, so it has become one attention topic in remote sensing application. Since saturation problem existed in the inversion of biomass by spectral features, the objective of this study was to propose a 'image and spectrum' fusion index by integrating the biomass-related texture feature into vegetation index.【Method】In this study, the extracted spectral indices and texture features from UAV hyperspectral imagery were used to first construct biomass models, respectively. Secondly, the 'image and spectrum' fusion indices by combining (multiplying or dividing) the biomass-sensitive vegetation index and texture feature were established to solve the saturation problem by spectral information and texture information fusion and to construct biomass model. Finally, the estimation effect of the biomass model constructed by different indices was compared, and then analyze the ability of the 'image and spectrum' fusion indices to estimate biomass.【Result】(1) The vegetation index was almost saturated when LAI was no larger than 5, while these 'image and spectrum' fusion indices, VI×658, VI/658, VI/658, VI/658, VI/514, VI/514, VI/514, VI×802, VI×802, began to perform saturation when at LAI was larger than 5. Compared with the vegetation index, the anti-saturation ability of the 'image and spectrum' fusion index was improved obviously. (2) Compared with the vegetation index (excepting for GNDVI、NDVI), the anti-saturation ability of the 'image and spectrum' fusion indices (VI×658, VI/658, VI/658, VI/658, VI/514, VI/514, VI/514, VI×802, VI×802) improved effectively, and their correlations with biomass improved as well. Meanwhile the biomass model based on the 'image and spectrum' fusion indices performed well, withRandvalues of 0.81 and 826.02 kg·hm-2, respectively. (3) Compared with spectral index and texture feature, biomass model accuracy by 'image and spectrum' fusion index (R=0.81) was significantly higher than that of the vegetation index (R= 0.69) and texture feature (R= 0.71).【Conclusion】Results showed that both the anti-saturation ability and the accuracy of biomass model constructed by the 'image and spectrum' fusion index were effectively improved, which indicated that spectral information and texture information fusion could achieve a great estimation of winter wheat biomass. The research provided a new way for quantitative inversion of biomass.

    biomass; 'image and spectrum' fusion index; texture feature; saturation; winter wheat

    2018-02-05;

    2018-05-29

    國家自然科學(xué)基金(61661136003,41471285)、國家重點研究計劃(2016YFD0300602-04,2016YFD0300603-05)

    劉暢,Tel:18292488506;E-mail:1224129134@qq.com。

    李振海,E-mail:lizh323@126.com

    10.3864/j.issn.0578-1752.2018.16.003

    猜你喜歡
    植被指數(shù)冬小麥紋理
    基于BM3D的復(fù)雜紋理區(qū)域圖像去噪
    軟件(2020年3期)2020-04-20 01:45:18
    使用紋理疊加添加藝術(shù)畫特效
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    TEXTURE ON TEXTURE質(zhì)地上的紋理
    Coco薇(2017年8期)2017-08-03 15:23:38
    甘肅冬小麥田
    植物保護(2017年1期)2017-02-13 06:44:34
    消除凹凸紋理有妙招!
    Coco薇(2015年5期)2016-03-29 23:22:15
    主要植被指數(shù)在生態(tài)環(huán)評中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
    日韩熟女老妇一区二区性免费视频| 婷婷色综合www| 大香蕉久久成人网| 欧美日韩综合久久久久久| 免费看不卡的av| 国产野战对白在线观看| 一级黄片播放器| 亚洲精品第二区| 精品99又大又爽又粗少妇毛片| 亚洲国产色片| 久久精品国产综合久久久| 国精品久久久久久国模美| 精品少妇久久久久久888优播| 中国国产av一级| 亚洲精品乱久久久久久| 26uuu在线亚洲综合色| 久久免费观看电影| 亚洲av电影在线进入| 国产精品麻豆人妻色哟哟久久| 国产在线一区二区三区精| 亚洲欧美色中文字幕在线| 久久国产精品男人的天堂亚洲| 高清在线视频一区二区三区| 国产黄色视频一区二区在线观看| 亚洲精品av麻豆狂野| 久久久久精品久久久久真实原创| 国产一区二区三区av在线| 成人毛片a级毛片在线播放| 大话2 男鬼变身卡| 国产一区二区三区综合在线观看| 久久精品熟女亚洲av麻豆精品| 乱人伦中国视频| 新久久久久国产一级毛片| 亚洲第一av免费看| 国产xxxxx性猛交| 久久av网站| 女人高潮潮喷娇喘18禁视频| 久久国产亚洲av麻豆专区| 午夜激情久久久久久久| 欧美日韩视频精品一区| 免费高清在线观看日韩| 精品人妻在线不人妻| 香蕉精品网在线| 99热全是精品| 日韩中文字幕欧美一区二区 | 丰满饥渴人妻一区二区三| 欧美少妇被猛烈插入视频| 久久99蜜桃精品久久| 叶爱在线成人免费视频播放| 蜜桃在线观看..| www.自偷自拍.com| 少妇人妻精品综合一区二区| 亚洲五月色婷婷综合| 热re99久久国产66热| 亚洲欧洲精品一区二区精品久久久 | 欧美日韩一级在线毛片| 麻豆av在线久日| 只有这里有精品99| 中文字幕另类日韩欧美亚洲嫩草| 国产毛片在线视频| 赤兔流量卡办理| 国产又色又爽无遮挡免| 国产又色又爽无遮挡免| 99re6热这里在线精品视频| 国产亚洲av片在线观看秒播厂| 精品99又大又爽又粗少妇毛片| 久久久久久久大尺度免费视频| 天天躁日日躁夜夜躁夜夜| 毛片一级片免费看久久久久| 黄色怎么调成土黄色| 久久综合国产亚洲精品| 国产黄色免费在线视频| 国产欧美亚洲国产| 午夜福利乱码中文字幕| 亚洲国产精品国产精品| 亚洲精品一区蜜桃| 婷婷成人精品国产| 国产精品国产三级国产专区5o| av女优亚洲男人天堂| 亚洲成人手机| 亚洲国产欧美网| videossex国产| 亚洲国产最新在线播放| 亚洲 欧美一区二区三区| 99热网站在线观看| 免费观看a级毛片全部| 日本免费在线观看一区| 大陆偷拍与自拍| 人人妻人人爽人人添夜夜欢视频| 精品亚洲成a人片在线观看| 国产精品国产三级国产专区5o| 亚洲欧美清纯卡通| 黄色怎么调成土黄色| √禁漫天堂资源中文www| 亚洲伊人久久精品综合| 免费人妻精品一区二区三区视频| 国产有黄有色有爽视频| 国产欧美亚洲国产| 1024香蕉在线观看| 欧美av亚洲av综合av国产av | 男男h啪啪无遮挡| av天堂久久9| 国产精品人妻久久久影院| 成人国语在线视频| 十分钟在线观看高清视频www| 日韩制服骚丝袜av| 久久99蜜桃精品久久| 一区二区三区精品91| 一级,二级,三级黄色视频| 精品人妻偷拍中文字幕| 日本色播在线视频| 日本-黄色视频高清免费观看| 在线精品无人区一区二区三| 在线精品无人区一区二区三| 欧美av亚洲av综合av国产av | 精品一区二区三区四区五区乱码 | 成人黄色视频免费在线看| 欧美日韩视频精品一区| 亚洲av成人精品一二三区| 国产成人91sexporn| 成人国产麻豆网| 极品人妻少妇av视频| 欧美日韩精品成人综合77777| 婷婷色综合大香蕉| 欧美日韩视频精品一区| 少妇被粗大猛烈的视频| 性少妇av在线| 亚洲经典国产精华液单| av线在线观看网站| 纯流量卡能插随身wifi吗| 婷婷色av中文字幕| 天天躁夜夜躁狠狠躁躁| 老汉色∧v一级毛片| 亚洲av综合色区一区| 久久精品夜色国产| 建设人人有责人人尽责人人享有的| 久久人人爽人人片av| 性色av一级| 亚洲av在线观看美女高潮| 久久狼人影院| 90打野战视频偷拍视频| 国产福利在线免费观看视频| 一区二区三区精品91| 成人国语在线视频| 我的亚洲天堂| 精品少妇一区二区三区视频日本电影 | 成年人免费黄色播放视频| 午夜福利在线免费观看网站| 成年美女黄网站色视频大全免费| 观看av在线不卡| 国产精品 欧美亚洲| av卡一久久| 黄网站色视频无遮挡免费观看| 亚洲av电影在线观看一区二区三区| 午夜激情久久久久久久| 成人毛片60女人毛片免费| 丝袜美腿诱惑在线| 亚洲国产看品久久| 一区二区日韩欧美中文字幕| 午夜免费观看性视频| 国产免费又黄又爽又色| 久久久久久久亚洲中文字幕| 亚洲国产精品国产精品| videos熟女内射| 18禁观看日本| 涩涩av久久男人的天堂| 丁香六月天网| 日韩电影二区| 亚洲人成网站在线观看播放| 18+在线观看网站| 欧美人与善性xxx| 国产av一区二区精品久久| 国产精品国产av在线观看| 80岁老熟妇乱子伦牲交| 亚洲人成网站在线观看播放| 国产毛片在线视频| 女人久久www免费人成看片| 久久久亚洲精品成人影院| 亚洲精品久久久久久婷婷小说| 亚洲一码二码三码区别大吗| 天天操日日干夜夜撸| 日韩中字成人| 激情五月婷婷亚洲| 热99国产精品久久久久久7| 国产乱来视频区| 久久精品国产亚洲av天美| 久久婷婷青草| 黄色配什么色好看| 国产熟女欧美一区二区| 男女下面插进去视频免费观看| 69精品国产乱码久久久| 精品久久蜜臀av无| 青青草视频在线视频观看| 美国免费a级毛片| 激情视频va一区二区三区| 一个人免费看片子| 亚洲国产精品999| 日本免费在线观看一区| 国产黄频视频在线观看| 极品人妻少妇av视频| 伦理电影免费视频| 亚洲成人手机| 蜜桃国产av成人99| 日韩视频在线欧美| 两个人免费观看高清视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 熟妇人妻不卡中文字幕| 精品视频人人做人人爽| 欧美变态另类bdsm刘玥| 亚洲成色77777| 91午夜精品亚洲一区二区三区| 日韩伦理黄色片| 有码 亚洲区| 免费观看无遮挡的男女| 免费在线观看视频国产中文字幕亚洲 | 中国三级夫妇交换| 日韩中文字幕视频在线看片| 亚洲av.av天堂| 少妇人妻久久综合中文| 夫妻午夜视频| 欧美精品国产亚洲| 美国免费a级毛片| 9191精品国产免费久久| 亚洲av在线观看美女高潮| 久久毛片免费看一区二区三区| 精品国产露脸久久av麻豆| 久久精品久久久久久噜噜老黄| 少妇猛男粗大的猛烈进出视频| 成年av动漫网址| 国产在线视频一区二区| 久久精品亚洲av国产电影网| 一级毛片电影观看| 人人澡人人妻人| 2018国产大陆天天弄谢| av电影中文网址| 美国免费a级毛片| 国产欧美日韩一区二区三区在线| 激情视频va一区二区三区| 妹子高潮喷水视频| 丝袜在线中文字幕| 国产精品蜜桃在线观看| 少妇精品久久久久久久| 涩涩av久久男人的天堂| 欧美av亚洲av综合av国产av | 一二三四在线观看免费中文在| 在线免费观看不下载黄p国产| 亚洲精品中文字幕在线视频| 青春草国产在线视频| 女人高潮潮喷娇喘18禁视频| 亚洲av在线观看美女高潮| 久久狼人影院| h视频一区二区三区| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 国产成人精品无人区| 成人手机av| 午夜av观看不卡| 五月伊人婷婷丁香| 熟女少妇亚洲综合色aaa.| 亚洲国产欧美在线一区| 搡老乐熟女国产| 国产精品一区二区在线不卡| 黄色怎么调成土黄色| 2021少妇久久久久久久久久久| 一本大道久久a久久精品| 国产 精品1| 侵犯人妻中文字幕一二三四区| 亚洲五月色婷婷综合| 精品亚洲成a人片在线观看| 波多野结衣av一区二区av| 深夜精品福利| 精品人妻熟女毛片av久久网站| 成人免费观看视频高清| 亚洲av日韩在线播放| a级毛片在线看网站| 视频区图区小说| 国产精品秋霞免费鲁丝片| 少妇的丰满在线观看| 女性被躁到高潮视频| 999久久久国产精品视频| 黄色一级大片看看| 狠狠婷婷综合久久久久久88av| www.精华液| 国产日韩欧美视频二区| 黄色配什么色好看| 国产黄色免费在线视频| 日本vs欧美在线观看视频| 中文字幕色久视频| 国产亚洲最大av| 国产一区二区 视频在线| 午夜老司机福利剧场| 另类亚洲欧美激情| 2022亚洲国产成人精品| 久久影院123| 女的被弄到高潮叫床怎么办| 久久女婷五月综合色啪小说| av免费观看日本| 熟女少妇亚洲综合色aaa.| av国产精品久久久久影院| 国产av国产精品国产| 国产精品久久久久成人av| 亚洲欧美一区二区三区国产| 在线观看www视频免费| 看十八女毛片水多多多| 日韩制服骚丝袜av| 男女高潮啪啪啪动态图| 免费观看无遮挡的男女| 中文字幕精品免费在线观看视频| 亚洲国产精品成人久久小说| 韩国精品一区二区三区| 一级爰片在线观看| 成人毛片a级毛片在线播放| 精品一区二区三区四区五区乱码 | 校园人妻丝袜中文字幕| 亚洲欧美一区二区三区国产| 999精品在线视频| 午夜福利视频在线观看免费| 成年女人在线观看亚洲视频| 亚洲成av片中文字幕在线观看 | 亚洲人成电影观看| 久久久久视频综合| 国产片内射在线| 亚洲国产色片| 亚洲第一青青草原| 最黄视频免费看| 美女福利国产在线| 国产乱人偷精品视频| 美女脱内裤让男人舔精品视频| 日日撸夜夜添| 日本欧美国产在线视频| 亚洲国产精品国产精品| 黑丝袜美女国产一区| 日韩大片免费观看网站| 日本欧美国产在线视频| 久久精品熟女亚洲av麻豆精品| 久久免费观看电影| 91久久精品国产一区二区三区| 另类精品久久| 久久精品国产亚洲av涩爱| 国产欧美日韩综合在线一区二区| 国语对白做爰xxxⅹ性视频网站| 一级爰片在线观看| 久久99一区二区三区| av又黄又爽大尺度在线免费看| 丝袜美足系列| 三上悠亚av全集在线观看| 亚洲成色77777| 少妇的丰满在线观看| 国产成人aa在线观看| 欧美日本中文国产一区发布| xxx大片免费视频| 韩国高清视频一区二区三区| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 啦啦啦中文免费视频观看日本| 在线亚洲精品国产二区图片欧美| 国产精品.久久久| av福利片在线| 久久热在线av| 国产亚洲一区二区精品| 亚洲,欧美,日韩| 美女主播在线视频| 夜夜骑夜夜射夜夜干| 国产精品国产av在线观看| 菩萨蛮人人尽说江南好唐韦庄| 免费日韩欧美在线观看| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| 两个人看的免费小视频| 老汉色∧v一级毛片| 精品酒店卫生间| 国产国语露脸激情在线看| 久久久久精品久久久久真实原创| 国产视频首页在线观看| 少妇熟女欧美另类| 另类亚洲欧美激情| 亚洲av国产av综合av卡| 国产视频首页在线观看| 99热全是精品| 亚洲,欧美精品.| 婷婷色麻豆天堂久久| 久久ye,这里只有精品| 亚洲精品成人av观看孕妇| 18+在线观看网站| 久久99精品国语久久久| 久久青草综合色| 秋霞伦理黄片| av网站免费在线观看视频| 午夜福利视频精品| 久久毛片免费看一区二区三区| 成人黄色视频免费在线看| 精品久久久精品久久久| 欧美黄色片欧美黄色片| 五月伊人婷婷丁香| 啦啦啦在线免费观看视频4| 韩国精品一区二区三区| 99精国产麻豆久久婷婷| 熟女少妇亚洲综合色aaa.| 999精品在线视频| 男女免费视频国产| xxxhd国产人妻xxx| 日本欧美视频一区| a级毛片在线看网站| 午夜日韩欧美国产| 两个人免费观看高清视频| 日韩三级伦理在线观看| 亚洲国产欧美日韩在线播放| 久久久久久久久久人人人人人人| 成年女人毛片免费观看观看9 | 国产精品 欧美亚洲| 美女视频免费永久观看网站| 香蕉丝袜av| 观看美女的网站| 欧美av亚洲av综合av国产av | 亚洲av欧美aⅴ国产| 色哟哟·www| 亚洲精品中文字幕在线视频| 国产男女超爽视频在线观看| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久久久免| 99精国产麻豆久久婷婷| 精品一区在线观看国产| 欧美成人精品欧美一级黄| 久久久久久久久久久免费av| 一本色道久久久久久精品综合| 亚洲精品第二区| 捣出白浆h1v1| 爱豆传媒免费全集在线观看| 韩国av在线不卡| 久久久久国产网址| 亚洲av.av天堂| 在线观看www视频免费| 久久久国产精品麻豆| 黑人欧美特级aaaaaa片| 国产免费又黄又爽又色| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| 九九爱精品视频在线观看| 精品国产乱码久久久久久男人| 高清av免费在线| 天堂俺去俺来也www色官网| 九草在线视频观看| 色婷婷av一区二区三区视频| 老司机影院成人| kizo精华| 免费播放大片免费观看视频在线观看| 日韩免费高清中文字幕av| 人妻 亚洲 视频| 五月伊人婷婷丁香| 婷婷色综合www| 另类亚洲欧美激情| 中文字幕人妻熟女乱码| 春色校园在线视频观看| 少妇的逼水好多| 日韩伦理黄色片| 欧美人与性动交α欧美软件| 久久久精品免费免费高清| 国产精品成人在线| 精品国产国语对白av| 一区二区日韩欧美中文字幕| 国产熟女午夜一区二区三区| 亚洲国产最新在线播放| 少妇的逼水好多| 最新中文字幕久久久久| 中文字幕制服av| 国产精品不卡视频一区二区| 精品国产一区二区久久| 国产精品嫩草影院av在线观看| 久久鲁丝午夜福利片| 欧美少妇被猛烈插入视频| 成年av动漫网址| 麻豆av在线久日| 亚洲国产最新在线播放| 热99久久久久精品小说推荐| 国产麻豆69| 国产成人aa在线观看| 久久精品人人爽人人爽视色| 丰满乱子伦码专区| 女的被弄到高潮叫床怎么办| 男女下面插进去视频免费观看| 国精品久久久久久国模美| 人妻 亚洲 视频| 成人国产av品久久久| 国产精品麻豆人妻色哟哟久久| 亚洲av欧美aⅴ国产| videossex国产| 免费在线观看完整版高清| 国产精品亚洲av一区麻豆 | 青春草国产在线视频| 人人妻人人澡人人爽人人夜夜| 可以免费在线观看a视频的电影网站 | 七月丁香在线播放| 在线看a的网站| 亚洲色图综合在线观看| 成人午夜精彩视频在线观看| 一区二区av电影网| 国产日韩欧美视频二区| 欧美日韩一区二区视频在线观看视频在线| av国产精品久久久久影院| 一区福利在线观看| 久久精品久久久久久噜噜老黄| 秋霞在线观看毛片| 成人毛片60女人毛片免费| 国产在视频线精品| 夫妻性生交免费视频一级片| 欧美精品av麻豆av| 一级爰片在线观看| 久久这里只有精品19| 婷婷色av中文字幕| 国产极品天堂在线| 黄色视频在线播放观看不卡| 最近中文字幕2019免费版| 毛片一级片免费看久久久久| 亚洲第一青青草原| 国产精品免费大片| 国产日韩欧美在线精品| 久久久久久伊人网av| 国产欧美日韩一区二区三区在线| 亚洲精品,欧美精品| 午夜激情av网站| 亚洲精品aⅴ在线观看| 最近手机中文字幕大全| 日韩,欧美,国产一区二区三区| 欧美精品av麻豆av| h视频一区二区三区| 中文精品一卡2卡3卡4更新| 在线观看一区二区三区激情| 国产一区亚洲一区在线观看| 国产xxxxx性猛交| 国产激情久久老熟女| 欧美97在线视频| 色吧在线观看| 999久久久国产精品视频| 国产精品成人在线| 亚洲国产色片| 大陆偷拍与自拍| 日日爽夜夜爽网站| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 女人高潮潮喷娇喘18禁视频| 婷婷色综合大香蕉| 久久人人爽人人片av| 男女下面插进去视频免费观看| 9191精品国产免费久久| 亚洲欧美成人精品一区二区| 欧美日韩亚洲国产一区二区在线观看 | 色94色欧美一区二区| 国产精品女同一区二区软件| 久热这里只有精品99| 一边亲一边摸免费视频| 亚洲国产精品国产精品| 日韩精品有码人妻一区| 美女视频免费永久观看网站| 免费高清在线观看视频在线观看| 国产成人aa在线观看| 国产免费福利视频在线观看| 亚洲国产最新在线播放| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 另类亚洲欧美激情| 日韩av免费高清视频| 一级爰片在线观看| 少妇人妻久久综合中文| 国产精品蜜桃在线观看| 中文欧美无线码| 欧美成人午夜精品| 婷婷色麻豆天堂久久| 妹子高潮喷水视频| 免费观看a级毛片全部| 欧美 日韩 精品 国产| 亚洲欧洲国产日韩| 国产精品久久久久久久久免| 午夜福利网站1000一区二区三区| 一区二区三区乱码不卡18| 人体艺术视频欧美日本| 国产免费一区二区三区四区乱码| 90打野战视频偷拍视频| 日韩欧美精品免费久久| 国产男女超爽视频在线观看| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 老司机亚洲免费影院| 亚洲精品aⅴ在线观看| 成年人午夜在线观看视频| 欧美精品一区二区大全| 丝袜脚勾引网站| 亚洲精品视频女| 亚洲av综合色区一区| 少妇人妻久久综合中文| 亚洲欧美色中文字幕在线| 国产成人免费观看mmmm| 国产一区有黄有色的免费视频| 中文精品一卡2卡3卡4更新| 国产精品香港三级国产av潘金莲 | 久久久久久人妻| 精品国产一区二区久久| 97在线人人人人妻| 亚洲综合色惰| 精品亚洲成a人片在线观看| 国产男人的电影天堂91| 两个人免费观看高清视频| 国产亚洲最大av| kizo精华| 亚洲精品久久午夜乱码| 五月开心婷婷网| 国产精品久久久久久av不卡| 色94色欧美一区二区| 午夜91福利影院| 香蕉丝袜av| 男人爽女人下面视频在线观看| 母亲3免费完整高清在线观看 | 日本91视频免费播放| 可以免费在线观看a视频的电影网站 | 久久人人爽av亚洲精品天堂| 欧美日本中文国产一区发布|