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

    基于ICESat-GLAS數(shù)據(jù)和模糊模式識(shí)別算法識(shí)別森林類型

    2021-07-28 04:06:56蔡龍濤邢艷秋丁建華黃佳鵬
    關(guān)鍵詞:針葉林闊葉林混交林

    蔡龍濤,邢 濤,邢艷秋,丁建華,黃佳鵬,崔 陽(yáng),秦 磊

    (東北林業(yè)大學(xué)工程技術(shù)學(xué)院, 黑龍江 哈爾濱 150040)

    精確估測(cè)森林生物量,可以為陸地生態(tài)系統(tǒng)碳儲(chǔ)量及碳循環(huán)研究提供重要的數(shù)據(jù)支撐[1]。為提高森林生物量估測(cè)精度,現(xiàn)有的研究多采用不同森林類型生物量估測(cè)法,估算整體森林生物量[2-4]。因此,在森林生物量估測(cè)之前有必要對(duì)森林類型進(jìn)行分類研究。

    為提高森林類型分類精度,本研究通過(guò)對(duì)不同森林類型冠層分布特征以及GLAS冠層回波波形進(jìn)行研究分析,提取了波形特征參數(shù)并結(jié)合不同森林類型冠層分布特征,選取并聯(lián)合其他波形特征參數(shù)[16,18],建立波形特征參數(shù)組合,最后基于模糊模式識(shí)別算法,實(shí)現(xiàn)了針葉林和闊葉林以及針葉林、闊葉林和混交林森林類型高精度分類研究。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于吉林省汪清林業(yè)局,屬于長(zhǎng)白山系中低山區(qū)(129°56′~131°04′E,43°05′~43°40′N),地處寒溫帶,總面積30.4萬(wàn) hm2,南北長(zhǎng)約60 km,東西長(zhǎng)約85 km。其中地面高程為360~1 477 m,坡度變化范圍0°~45°。林區(qū)內(nèi)森林覆蓋率達(dá)到95.95%,深山區(qū)林相以針闊混交林為主,呈帶狀分布于海拔500~1 100 m。針葉樹主要有紅松(Pinuskoraiensis)、云杉(Piceaasperata)、臭松(Abiesnephrolepis)和落葉松(Larixgmelinii),闊葉樹多為椴樹(Tiliatuan)、蒙古櫟(Quercusmomglica)、楓樺(Betulacosta)、色木槭(Acermono)和白樺(Betulaplatyphylla)等[19]。

    1.2 GLAS波形數(shù)據(jù)及野外調(diào)查數(shù)據(jù)

    GLAS回波波形數(shù)據(jù)可從美國(guó)國(guó)家冰雪數(shù)據(jù)中心(http://nsidc.org/data/ice-sat/)免費(fèi)下載,該中心擁有ICESat衛(wèi)星2003年至2009年采集的所有GLAS回波波形數(shù)據(jù)。本研究所用數(shù)據(jù)為2006年10月、2009年10月和2010年10月的GLA01和GLA14數(shù)據(jù)產(chǎn)品,可通過(guò)經(jīng)緯度查找樣地所對(duì)應(yīng)的GLAS回波波形[19]。

    野外調(diào)查數(shù)據(jù)分別于2006年9月、2007年9月和2010年9月在長(zhǎng)白山汪清林區(qū)采集獲取。具體采集方法參見文獻(xiàn)[19]。3年野外實(shí)測(cè)調(diào)查共收集了286組樣地?cái)?shù)據(jù)。其中闊葉林樣地137組,針葉林樣地50組,混交林樣地為60組,其他類型樣地如水地、裸地和草地等共計(jì)39組。在對(duì)森林類型分類時(shí),結(jié)合我國(guó)森林資源調(diào)查主要技術(shù)規(guī)定把針葉林蓄積量比例占總蓄積量65%以上的林分定義為針葉林,闊葉林蓄積量比例占總蓄積量65%以上的林分定義為闊葉林,任何一個(gè)樹種蓄積量占總蓄積量不到65%的林分定義為混交林[19]。

    1.3 波形數(shù)據(jù)處理

    2012年11月1日美國(guó)冰雪數(shù)據(jù)中心提供了HDF5格式的GLAS波形數(shù)據(jù),HDF5文件中GLAS波形數(shù)據(jù)以十進(jìn)制形式進(jìn)行存儲(chǔ),省去了最初的GLAS波形數(shù)據(jù)二進(jìn)制轉(zhuǎn)換等操作。然而GLAS波形數(shù)據(jù)中仍存在噪聲數(shù)據(jù),需要對(duì)波形數(shù)據(jù)進(jìn)行去噪處理。常見的波形去噪算法有小波去噪算法[20-21]和高斯去噪算法[22-23]。但是小波去噪和高斯去噪算法不可避免地會(huì)濾除GLAS回波波形中部分信號(hào)數(shù)據(jù),為減少GLAS回波波形數(shù)據(jù)中信號(hào)數(shù)據(jù)被過(guò)多濾除,本研究直接利用高斯擬合方程對(duì)GLAS回波波形進(jìn)行去噪處理。

    另外,在定義GLAS回波波形背景噪聲閾值時(shí),本研究將其定義為波形中前100幀能量值均值與其4倍標(biāo)準(zhǔn)差之和;然而,研究時(shí)發(fā)現(xiàn)GLAS回波波形中存在有效回波波形(攜帶地物信息的GLAS回波波形)集中于波形前半部分的情況,針對(duì)這類波形數(shù)據(jù),把背景噪聲閾值定義為波形中后100幀能量值均值與其4倍標(biāo)準(zhǔn)差之和。確定背景噪聲閾值后利用高斯擬合方程對(duì)GLAS回波波形進(jìn)行去噪處理,然后基于去噪后的擬合波形提取波形特征參數(shù)。

    地形坡度較大時(shí),GLAS回波波形會(huì)出現(xiàn)波形展寬現(xiàn)象[24-27],影響波形特征參數(shù)提取精度。為削弱地形坡度對(duì)森林類型識(shí)別結(jié)果的影響,本研究把有效回波波形中右半部分回波波形能量值最大處定義為地面回波波形,并把該地面回波波形作為GLAS光斑內(nèi)實(shí)際地面回波波形[19]。

    1.4 波形特征參數(shù)提取

    針葉樹葉片如針,呈扁平尖細(xì)狀,冠形多為圓錐體形;闊葉樹葉片平展寬大,冠形多為橢球體、上半球體和下半球體形。受冠形形態(tài)影響,針葉林冠層頂部枝葉較為稀疏,反映在GLAS回波波形上為靠近波形起波點(diǎn)位置處的能量值偏低;闊葉林冠層多集中于林層頂部,反映在GLAS回波波形上為靠近波形起波點(diǎn)位置處的能量值偏高;而混交林林分內(nèi)闊葉樹和針葉樹個(gè)數(shù)相當(dāng),使得混交林冠層對(duì)應(yīng)的GLAS冠層回波波形能量值呈均勻分布。

    依據(jù)針葉林、闊葉林和混交林冠層頂部對(duì)應(yīng)在GLAS回波波形上的能量分布差異提取波形特征參數(shù),波形特征參數(shù)提取示意圖如圖1所示。

    圖1 波形特征參數(shù)提取示意圖

    參數(shù)提取公式如下所示。

    (1)

    (2)

    (3)

    (4)

    (5)

    據(jù)此,將Zhang等[16]提出波形特征參數(shù)AGS(回波波形分解后高斯擬合分量斜率值)、SGS(回波波形分解后高斯擬合分量斜率標(biāo)準(zhǔn)差)和MSGS(回波波形分解后高斯擬合分量斜率標(biāo)準(zhǔn)差校正值),與本研究提出的波形特征參數(shù)相結(jié)合,建立新的波形特征參數(shù)組合,用于森林類型分類研究。其中,波形特征參數(shù)AGS[記為σ(AGS)]、SGS[記為σ(SGS)]和MSGS[記為σ(MSGS)]計(jì)算方法如下所示。

    (6)

    (7)

    (8)

    (9)

    (10)

    1.5 模糊模式識(shí)別算法

    1.5.1 指標(biāo)歸一化

    因提取方式不同,波形特征參數(shù)之間存在量綱差異。為避免波形特征參數(shù)因量綱不同而影響森林類型分類精度,在對(duì)森林類型分類之前,先對(duì)波形特征參數(shù)進(jìn)行了指標(biāo)歸一化處理。其中,指標(biāo)歸一化公式為:

    (11)

    式中:xi為歸一化后的波形特征參數(shù)值;x′i為歸一化之前波形特征參數(shù)值;x′min為某一類波形特征參數(shù)最小值;x′max為某一類波形特征參數(shù)最大值。

    1.5.2 奇異點(diǎn)檢測(cè)

    受地形、噪聲以及森林自身等因素影響,提取的波形特征參數(shù)中存在奇異點(diǎn)樣本數(shù)據(jù),使得本研究建立的森林類型模式與實(shí)際森林類型模式存在一定差異,影響森林類型分類精度。因此,在建立不同森林類型模式之前,需剔除樣本數(shù)據(jù)中的奇異點(diǎn)樣本。而本研究在對(duì)波形特征參數(shù)進(jìn)行奇異點(diǎn)檢測(cè)時(shí),選用Jihyun等[28]提出的重心不穩(wěn)定系數(shù)法對(duì)樣本數(shù)據(jù)進(jìn)行奇異點(diǎn)檢測(cè),然后予以剔除。其中,重心不穩(wěn)定系數(shù)計(jì)算公式為:

    (12)

    Δθi(p)=|θi(p)-θi+1(p)|,

    i=1,2,…,k-2;

    (13)

    θi(p)=d(mi,mi+1),i=1,2,…,k-1;

    (14)

    (15)

    xq={xq1,xq2,…,xqn}。

    (16)

    式中:INSF為重心不穩(wěn)定系數(shù);Δθi為隨著k值增加兩個(gè)相鄰θi的差值;θi(p)為p的i鄰域與i+1鄰域幾何中心距離值,隨著k值增加,相鄰幾何中心距離值會(huì)發(fā)生變化;d(mi,mi+1)為mi與mi+1之間距離;mk為p的k鄰域k+1個(gè)點(diǎn)的幾何中心;p為某測(cè)試集樣本;k為測(cè)試集樣本p周圍的測(cè)試集樣本個(gè)數(shù);xq為某測(cè)試樣本特征參數(shù)數(shù)據(jù)集合;n為測(cè)試樣本q特征參數(shù)個(gè)數(shù)值。

    1.5.3 計(jì)算樣本均值

    奇異點(diǎn)檢測(cè)之后,需選取不同森林類型訓(xùn)練集樣本數(shù)據(jù)和測(cè)試集樣本數(shù)據(jù),利用選取的訓(xùn)練集樣本數(shù)據(jù)計(jì)算不同森林類型樣本均值。本研究在對(duì)訓(xùn)練集樣本數(shù)據(jù)和測(cè)試集樣本數(shù)據(jù)進(jìn)行選取時(shí),按照2∶1的比例進(jìn)行隨機(jī)選取,這樣選取的樣本中每個(gè)森林類型約有1/3的樣本數(shù)據(jù)用于森林類型分類測(cè)試。其中,樣本均值cij計(jì)算公式為:

    (17)

    式中:cijl為訓(xùn)練集森林類型i波形特征參數(shù)j中第l個(gè)樣本波形特征參數(shù)值;ki為訓(xùn)練集森林類型i樣本個(gè)數(shù)。

    1.5.4 計(jì)算模糊模式隸屬度

    在計(jì)算測(cè)試集樣本針葉林、闊葉林和混交林森林類型隸屬度時(shí),首先計(jì)算測(cè)試集樣本與針葉林、闊葉林和混交林模式之間的距離,然后通過(guò)計(jì)算得到距離值進(jìn)而求測(cè)試集樣本針葉林、闊葉林和混交林隸屬度。若某測(cè)試集樣本針葉林隸屬度大于闊葉林和混交林隸屬度,該樣地則被識(shí)別為針葉林,否則被識(shí)別為闊葉林或混交林。另外,在計(jì)算測(cè)試集樣本與針葉林、闊葉林和混交林模式之間的距離時(shí),本研究采用了歐幾里德(Euclid)距離算法,計(jì)算公式為:

    (18)

    式中:di(u,ci)為測(cè)試集樣本u與不同森林類型模式ci之間的距離;uj為第j個(gè)樣本特征參數(shù)組合。

    利用計(jì)算得到的距離值,可求測(cè)試集樣本針葉林、闊葉林和混交林隸屬度,計(jì)算公式為:

    (19)

    式中:Ci(u)為樣本u隸屬于森林類型i的隸屬度;p為森林類型分類總個(gè)數(shù)。

    2 結(jié)果與分析

    2.1 不同林型奇異點(diǎn)檢測(cè)比較

    奇異點(diǎn)波形特征參數(shù)值與正常林分波形特征參數(shù)值相差較大,故需對(duì)針葉林、闊葉林和混交林波形特征參數(shù)值統(tǒng)一進(jìn)行奇異點(diǎn)檢測(cè)處理。奇異點(diǎn)檢測(cè)時(shí),按照闊葉林、混交林和針葉林順序排序,闊葉林、針葉林和混交林樣本個(gè)數(shù)分別為137、50和60,檢測(cè)結(jié)果如圖2所示。

    圖2 不穩(wěn)定系數(shù)(INSF)隨樣本點(diǎn)位變化示意圖

    從圖2可以看出,247個(gè)樣本中7個(gè)樣本不穩(wěn)定系數(shù)(INSF)值明顯高于其他樣本,樣本位數(shù)分別為121、149、153、194、221、230和235。按照奇異點(diǎn)檢測(cè)時(shí)森林類型排列順序,發(fā)現(xiàn)針葉林、闊葉林和混交林奇異點(diǎn)個(gè)數(shù)值分別為3、1和3,結(jié)合針葉林、闊葉林和混交林樣本個(gè)數(shù)值,得到不同森林類型奇異點(diǎn)個(gè)數(shù)占比大小排序?yàn)椋横樔~林>混交林>闊葉林。分析其原因:通過(guò)對(duì)不同森林類型冠層特征研究,發(fā)現(xiàn)針葉林林分冠層枝葉多集中于冠層底部,而本研究所用樣地為天然林地,地形坡度變化較大,受地形坡度影響,同一高程水平面針葉林冠層有效投影面積與地面有效投影面積重疊部分較多,反映在GLAS回波波形上為冠層回波波形與地面回波波形有較大面積重疊,以至于GLAS回波波形數(shù)據(jù)處理時(shí)難以高效提取地面回波波形,使得GLAS波形特征參數(shù)受地面回波波形影響較大,最終表現(xiàn)為針葉林波形特征參數(shù)中奇異點(diǎn)個(gè)數(shù)占比較高;相對(duì)于針葉林,闊葉林林分冠層枝葉多集中于冠層頂部,與地面在垂直方向上距離較大,表現(xiàn)在同一高程水平面上為闊葉林冠層有效投影面積與地面有效投影面積重疊部分較少,在GLAS回波波形上顯示為冠層回波波形與地面回波波形重疊面積較小,這種情況使得在對(duì)GLAS回波波形進(jìn)行數(shù)據(jù)處理時(shí)更容易識(shí)別并提取地面回波波形,最終表現(xiàn)為闊葉林波形特征參數(shù)中奇異點(diǎn)個(gè)數(shù)占比較低;而混交林林分中針葉樹和闊葉樹個(gè)數(shù)相當(dāng),冠層枝葉在林分冠層中呈均勻分布,使得混交林冠層回波波形與地面回波波形重疊面積高于針葉林,低于闊葉林,反映在GLAS回波波形上為混交林冠層回波能量值受地面回波波形影響大于針葉林,小于闊葉林,最終體現(xiàn)在波形特征參數(shù)奇異點(diǎn)個(gè)數(shù)占比上為混交林低于針葉林,高于闊葉林。

    2.2 不同波形特征參數(shù)對(duì)針葉林和闊葉林分類的影響

    結(jié)合奇異點(diǎn)檢測(cè)結(jié)果,從針葉林和闊葉林樣本中分別剔除3個(gè)和1個(gè)奇異點(diǎn)樣本,針葉林和闊葉林剩余樣本數(shù)分別為47和136;然后按照2∶1的比例隨機(jī)選取訓(xùn)練集和測(cè)試集樣本,選取之后針葉林和闊葉林訓(xùn)練集樣本數(shù)分別為33和96,測(cè)試集樣本數(shù)分別為14和40;最后利用選取的訓(xùn)練集和測(cè)試集樣本數(shù)據(jù)進(jìn)行針葉林和闊葉林森林類型分類研究。

    在對(duì)針葉林和闊葉林進(jìn)行森林類型識(shí)別研究時(shí),首先利用訓(xùn)練集樣本數(shù)據(jù),計(jì)算針葉林和闊葉林樣木均值;然后通過(guò)模糊識(shí)別算法最大隸屬度原則,計(jì)算測(cè)試集樣本針葉林和闊葉林隸屬度。當(dāng)某樣本針葉林隸屬度大于闊葉林隸屬度時(shí),該樣本被識(shí)別為針葉林,否則該樣本被識(shí)別為闊葉林。不同參數(shù)組合針葉林和闊葉林森林類型分類結(jié)果見表1,其中最優(yōu)特征參數(shù)組合分類結(jié)果見表2。

    表1 不同參數(shù)組合針葉林和闊葉林森林類型分類結(jié)果

    表2 基于波形特征參數(shù)組合MSGS的森林類型分類結(jié)果

    從表1可以看出,3種波形特征參數(shù)組合森林類型分類結(jié)果Kappa系數(shù)值都大于0.75時(shí),說(shuō)明森林類型分類結(jié)果與實(shí)際結(jié)果一致性較高,即通過(guò)模糊識(shí)別算法能夠有效識(shí)別針葉林和闊葉林。

    從表2可以得出,針葉林和闊葉林正確分類精度分別為92.86%和97.50%,針葉林森林類型正確分類精度低于闊葉林。分析其原因:首先,針葉林冠層枝葉多集中于林分冠層中部和底部,而且本研究所用樣地為天然林地,地形坡度變化較大,受地形坡度影響,針葉林冠層回波波形與地面回波波形重疊面積較大,在對(duì)GLAS回波波形進(jìn)行數(shù)據(jù)處理時(shí)難以高效提取地面回波波形,容易把地面回波波形誤分為針葉林回波波形;其次,針葉樹和闊葉樹樹高相同時(shí),針葉樹地面覆蓋面積小于闊葉樹,反映在GLAS回波波形上為針葉樹冠層回波波形能量值低于闊葉樹冠層回波波形能量值。相對(duì)于針葉林,闊葉林冠層枝葉多集中于林分冠層頂部,與地面具有相對(duì)較高的高程差,反映在GLAS回波波形上為冠層回波波形與地面回波波形具有較少的重疊面積;而且,當(dāng)闊葉樹和針葉樹樹高值相同時(shí),闊葉樹在地面覆蓋面積大于針葉樹,反映在GLAS回波波形上為闊葉樹冠層回波波形能量值高于針葉樹冠層回波波形能量值。兩種原因使得闊葉林冠層回波波形能量值高于針葉林,且冠層回波波形受地面回波波形影響低于針葉林,最終表現(xiàn)在森林類型分類精度上為闊葉林正確分類精度高于針葉林。

    2.3 不同林型分類結(jié)果比較

    結(jié)合奇異點(diǎn)檢測(cè)結(jié)果,從混交林樣地中剔除3個(gè)奇異點(diǎn)樣本,之后混交林剩余樣本個(gè)數(shù)為57。另外,在對(duì)混交林樣本進(jìn)行訓(xùn)練集和測(cè)試集選取時(shí),同樣采用2∶1選取原則,從混交林樣本中隨機(jī)選取40組樣本作為訓(xùn)練集樣本,17組樣本作為測(cè)試集樣本,然后結(jié)合選取的針葉林和闊葉林訓(xùn)練集以及測(cè)試集樣本數(shù)據(jù),對(duì)針葉林、闊葉林和混交林進(jìn)行森林類型分類研究。

    對(duì)針葉林、闊葉林和混交林進(jìn)行森林類型分類時(shí),首先計(jì)算針葉林、闊葉林和混交林樣本均值;然后采用模糊識(shí)別算法最大隸屬度原則,計(jì)算測(cè)試集樣本針葉林、闊葉林和混交林的隸屬度。當(dāng)某樣本針葉林隸屬度大于闊葉林和混交林隸屬度時(shí),該樣本被識(shí)別為針葉林;當(dāng)該樣本闊葉林隸屬度大于針葉林和混交林隸屬度時(shí),該樣本被識(shí)別為闊葉林;否則該樣本被識(shí)別為混交林。針葉林、闊葉林和混交林森林類型分類結(jié)果如表3所示。

    表3 針葉林、闊葉林和混交林分類結(jié)果比較

    從表3可以得出,混交林森林類型分類精度為52.94%,與闊葉林森林類型分類精度(97.50%)和針葉林森林類型分類精度(85.71%)相差較大。分析其原因:按照混交林定義(針葉樹和闊葉樹蓄積量占總蓄積量不到65%時(shí)林分定義為混交林),混交林林分內(nèi)針葉樹和闊葉樹個(gè)數(shù)相近,當(dāng)混交林林分內(nèi)存在優(yōu)勢(shì)樹種時(shí),優(yōu)勢(shì)樹種冠形特征在GLAS回波波形上表現(xiàn)更為明顯,提取的波形特征參數(shù)值與優(yōu)勢(shì)樹種波形特征參數(shù)值更為相近,使得在森林類型分類時(shí),容易將混交林誤分為針葉林或闊葉林。與混交林相比,針葉林和闊葉林林分內(nèi)優(yōu)勢(shì)樹種為針葉樹或闊葉樹,兩種類型的樹種冠形差別明顯,反映在GLAS回波波形上為針葉林和闊葉林波形特征參數(shù)值相差較大,利用相差較大的波形特征參數(shù)值識(shí)別針葉林和闊葉林,能夠提高針葉林和闊葉林森林類型分類精度。

    另外,從表3可以看出,混交林被誤分為闊葉林的個(gè)數(shù)多于針葉林。分析其原因:相對(duì)于針葉林,闊葉林枝葉多集中于林分冠層頂部,反映在GLAS回波波形上其波形特征更為明顯,所以對(duì)混交林進(jìn)行森林類型分類時(shí),未正確分類的樣本被識(shí)別為闊葉林的個(gè)數(shù)多于針葉林。

    加入混交林之后,森林類型分類總精度從96.30%下降到84.51%,Kappa系數(shù)從0.903 6下降到0.725 7。分析其原因:加入混交林之后,受混交林冠形特征影響,混交林正確分類個(gè)數(shù)與針葉林和闊葉林分類正確分類個(gè)數(shù)相差較大,且存在針葉林和闊葉林誤分為混交林的情況,使得森林類型分類總精度下降;另外,Kappa系數(shù)反映的是分類結(jié)果與實(shí)際結(jié)果一致性情況,而加入混交林之后,混交林分類精度較低,且存在針葉林和闊葉林誤分為混交林的情況,使得Kappa系數(shù)值較混交林加入之前明顯下降。

    3 討 論

    模糊識(shí)別算法能夠在波形數(shù)據(jù)處理過(guò)程中剔除波形特征參數(shù)中的奇異點(diǎn)數(shù)據(jù),并利用最大隸屬度的原則對(duì)森林類型進(jìn)行分類研究,提高森林類型識(shí)別精度。本研究通過(guò)模糊識(shí)別方法,利用Zhang等[16]提取的波形特征參數(shù)組合AGS、SGS、MSGS對(duì)針葉林和闊葉林進(jìn)行森林類型識(shí)別研究。研究表明,通過(guò)模糊識(shí)別方法針葉林和闊葉林森林類型識(shí)別精度為90.74%,優(yōu)于Zhang等[16]支持向量機(jī)分類精度。研究結(jié)果表明模糊模式識(shí)別方法在森林類型識(shí)別方面具有一定優(yōu)勢(shì)。

    然而,與針葉林和闊葉林森林類型識(shí)別精度相比,加入混交林后,森林類型識(shí)別精度較低,其中混交林森林類型識(shí)別精度為52.94%。為提高混交林森林類型識(shí)別精度,將嘗試從波形數(shù)據(jù)處理方面著手對(duì)混交林森林類型進(jìn)行分類研究,以便為GEDI波形數(shù)據(jù)和陸地生態(tài)碳衛(wèi)星波形數(shù)據(jù)森林類型分類研究提供技術(shù)支持。

    參考文獻(xiàn)(reference):

    [1]FIELD C B, RAUPACH M R, EBRARY I,et al. The global carbon cycle: integrating humans, climate, and the natural world[M]. Chicago: Island Press,2004.

    [2]湯旭光.基于激光雷達(dá)與多光譜遙感數(shù)據(jù)的森林地上生物量反演研究[D].北京:中國(guó)科學(xué)院研究生院(東北地理與農(nóng)業(yè)生態(tài)研究所),2013.TANG X G.Estimation of forest aboveground biomass by integrating ICESat/GLAS waveform and TM data[D].Beijing: Graduate School of Chinese Academy of Sciences (Northeast Institute of Geography and Agricultural Ecology),2013.

    [3]邱賽.ICESat-GLAS波形與HJ-1A高光譜影像聯(lián)合反演森林地上生物量的研究[D].哈爾濱:東北林業(yè)大學(xué),2016.QIU S.The research of regional forest above ground biomass inversion combining ICESat-GLAS waveform and HJ-1A hyperspectral imageries[D].Harbin:Northeast Forestry University,2016.

    [4]吳文躍,姚順彬,徐志揚(yáng).基于森林資源清查數(shù)據(jù)的江西省主要森林類型凈生產(chǎn)力研究[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,43(5):193-198.WU W Y,YAO S B,XU Z Y, et al.A study on the net productivity of main forest types in Jiangxi Province based on forest resources inventory data[J]. J Nanjing For Univ (Nat Sci Ed),2019,43(5):193-198. DOI: 10.3969/j.issn.1000-2006.201903045.

    [5]MOHOLDT G,NUTH C,HAGEN J O,et al.Recent elevation changes ofSvalbardglaciersderived from ICESat laser altimetry[J].Remote Sens Environ,2010,114(11):2756-2767.DOI:10.1016/j.rse.2010.06.008.

    [6]CARABAJAL C C,HARDING D J,SUCHDEO V P,et al.Development of AN ICESat geodetic control database and evaluation of global topographic assets[J]. American Geophysical Union, Fall Meeting, 2010(6):1907-1910. DOI:10.1038/nmeth.2835.

    [7]LEFSKY M A.A global forest canopy height map from the moderate resolution imaging spectroradiometer and the geoscience laser altimeter system[J].Geophys Res Lett,2010,37(15):L15401.DOI:10.1029/2010GL043622.

    [8]KHALEFA E,SMIT I P J,NICKLESS A,et al.Retrieval of savanna vegetation canopy height from ICESat-GLAS spaceborne LiDAR with terrain correction[J].IEEE Geosci Remote Sens Lett,2013,10(6):1439-1443.DOI:10.1109/LGRS.2013.2259793.

    [9]XING Y Q,DE GIER A,ZHANG J J,et al.An improved method for estimating forest canopy height using ICESat-GLAS full waveform data over sloping terrain:a case study in Changbai Mountains,China[J].Int J Appl Earth Obs Geoinformation,2010,12(5):385-392.DOI:10.1016/j.jag.2010.04.010.

    [10]WANG C,TANG F X,LI L W,et al.Wavelet analysis for ICESat/GLAS waveform decomposition and its application in average tree height estimation[J].IEEE Geosci Remote Sens Lett,2013,10(1):115-119.DOI:10.1109/LGRS.2012.2194692.

    [11]邱賽,邢艷秋,田靜,等.基于ICESat-GLAS波形數(shù)據(jù)估測(cè)森林郁閉度[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,40(5):99-106.QIU S,XING Y Q,TIAN J,et al.Estimation of forest canopy density based on ICESat-GLAS waveform data[J].J Nanjing For Univ (Nat Sci Ed),2016,40(5):99-106.DOI:10.3969/j.issn.1000-2006.2016.05.016.

    [12]XI X H,HAN T T,WANG C,et al.Forest above ground biomass inversion by fusing GLAS with optical remote sensing data[J].ISPRS Int J Geo-Inf,2016,5(4):45.DOI:10.3390/ijgi5040045.

    [13]WANG Y,NI W J,SUN G Q,et al.Slope-adaptive waveform metrics of large footprint lidar for estimation of forest aboveground biomass[J].Remote Sens Environ,2019,224:386-400.DOI:10.1016/j.rse.2019.02.017.

    [14]FATOYINBO T E,SIMARD M.Height and biomass of mangroves in Africa from ICESat/GLAS and SRTM[J].Int J Remote Sens,2013,34(2):668-681.DOI:10.1080/01431161.2012.712224.

    [15]RANSON K J,SUN G,KOVACS K,et al.Landcover attributes from ICESat GLAS data in Central Siberia[C]//IGARSS 2004.2004 IEEE International Geoscience and Remote Sensing Symposium.September 20-24,2004,Anchorage,AK,USA.IEEE,2004:753-756.DOI:10.1109/IGARSS.2004.1368511.

    [16]ZHANG J J,DE GIER A,XING Y Q,et al.Full waveform-based analysis for forest type information derivation from large footprint spaceborne lidar data[J].Photogramm Eng Remote Sensing,2011,77(3):281-290.DOI:10.14358/pers.77.3.281.

    [17]LI L C,XING Y Q.Study on forest type classification based on ICESat-GLAS lidar data[C]//2011 2nd International Conference on Artificial Intelligence,Management Science and Electronic Commerce (AIMSEC).August 8-10,2011,Deng Feng,China.IEEE,2011:3379-3382.DOI:10.1109/AIMSEC.2011.6011388.

    [18]劉美爽,邢艷秋,李立存,等.基于星載激光雷達(dá)數(shù)據(jù)和支持向量分類機(jī)方法的森林類型識(shí)別[J].東北林業(yè)大學(xué)學(xué)報(bào),2014,42(2):124-128.LIU M S,XING Y Q,LI L C,et al.Forest type identification with spaceborne LIDAR data and C-support vector classification[J].J Northeast For Univ,2014,42(2):124-128.DOI:10.13759/j.cnki.dlxb.2014.02.030.

    [19]蔡龍濤,邢艷秋 ,黃佳鵬,等. 基于 ICESat-GLAS 數(shù)據(jù)和回波仿真原理識(shí)別森林類型[J]. 中 南 林 業(yè) 科 技 大 學(xué) 學(xué) 報(bào),2021,41(1): 60-68. CAI L T, XING Y Q, HUANG J P, et al. Forest types identification based on ICESat-GLAS data and echo simulation principles[J]. J Central South Univ For Technol, 2021,41(1):60-68.

    [20]WU H B,XING Y Q.Wavelet transform and its application to ICESat-GLAS full waveform data[C]//2010 International Symposium on Intelligence Information Processing and Trusted Computing.October 28-29,2010,Huanggang,China.IEEE,2010:373-377.DOI:10.1109/IPTC.2010.38.

    [21]邱賽,邢艷秋,李立存,等.基于小波變換的ICESAT-GlAS波形處理[J].森林工程,2012,28(5):33-35,59.QIU S,XING Y Q,LI L C,et al.ICESAT-GLAS data processing based on wavelet transform[J].For Eng,2012,28(5):33-35,59.DOI:10.16270/j.cnki.slgc.2012.05.022.

    [22]LI D,XU L J,LI X L,et al.Full-waveform LiDAR signal filtering based on Empirical Mode Decomposition method[C]//2013 IEEE International Geoscience and Remote Sensing Symposium-IGARSS.July 21-26,2013,Melbourne,VIC,Australia.IEEE,2013:3399-3402.DOI:10.1109/IGARSS.2013.6723558.

    [23]YANG T,WANG C,LI G C,et al.Forest canopy height mapping over China using GLAS and MODIS data[J].Sci China Earth Sci,2015,58(1):96-105.DOI:10.1007/s11430-014-4905-5.

    [24]MAHONEY C,KLJUN N,LOS S,et al.Slope estimation from ICESat/GLAS[J].Remote Sens,2014,6(10):10051-10069.DOI:10.3390/rs61010051.

    [25]PANG Y,LEFSKY M,SUN G Q,et al.Impact of footprint diameter and off-nadir pointing on the precision of canopy height estimates from spaceborne lidar[J].Remote Sens Environ,2011,115(11):2798-2809.DOI:10.1016/j.rse.2010.08.025.

    [26]汪泉,宋文龍,張怡卓,等. 基于改進(jìn) VGG16 網(wǎng)絡(luò)的機(jī)載高光譜針葉樹種分類研究[J]. 森林工程,2021,37(3): 79-87.WANG Q,SONG W L,ZHANG Y Z,et al. Study on hyperspectral conifer species classification based on improved VGG16 network[J]. For Eng,2021,37(3): 79-87.

    [27]DUNCANSON L I,NIEMANN K O,WULDER M A.Estimating forest canopy height and terrain relief from GLAS waveform metrics[J].Remote Sens Environ,2010,114(1):138-154.DOI:10.1016/j.rse.2009.08.018.

    [28]JIHYUN H A,SEOK S,LEE J S.Robust outlier detection using the instability factor[J].Knowl-Based Syst,2014,63:15-23.DOI:10.1016/j.knosys.2014.03.001.

    猜你喜歡
    針葉林闊葉林混交林
    “China tendrá nuestro apoyo para ingresar al CPTPP”
    亞熱帶常綠闊葉林的世界之窗 錢江源,探路國(guó)家公園
    大興安嶺主要針葉樹種移植容器苗的培育技術(shù)與造林效果的研究
    施肥對(duì)油松人工林針葉養(yǎng)分變化的影響
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    營(yíng)造混交林的必要性和技術(shù)分析
    闊葉林培育及保護(hù)對(duì)策
    馬尾松、巨尾桉及其混交林土壤微量元素調(diào)查
    桉樹純林和桉-珍混交林水土流失規(guī)律的分析
    臺(tái)風(fēng)干擾對(duì)天童常綠闊葉林凋落物量的影響
    精品福利永久在线观看| 少妇的丰满在线观看| 亚洲av电影在线进入| 亚洲久久久国产精品| www.自偷自拍.com| 一二三四在线观看免费中文在| 精品欧美一区二区三区在线| 午夜福利在线观看吧| 丝袜美腿诱惑在线| 一夜夜www| 久久狼人影院| 亚洲精品在线观看二区| 黑人巨大精品欧美一区二区mp4| 99国产精品99久久久久| 久久久久久亚洲精品国产蜜桃av| 日韩av在线大香蕉| 俄罗斯特黄特色一大片| 亚洲第一青青草原| 亚洲国产毛片av蜜桃av| 每晚都被弄得嗷嗷叫到高潮| 欧美 亚洲 国产 日韩一| 窝窝影院91人妻| 成人亚洲精品av一区二区 | 两人在一起打扑克的视频| 一区二区三区精品91| 午夜福利在线免费观看网站| 欧美老熟妇乱子伦牲交| 国产av在哪里看| 亚洲三区欧美一区| 中文字幕高清在线视频| 亚洲欧洲精品一区二区精品久久久| 国产极品粉嫩免费观看在线| 久久99一区二区三区| 久久精品国产清高在天天线| 男女午夜视频在线观看| 国产亚洲精品久久久久久毛片| 99久久99久久久精品蜜桃| 免费久久久久久久精品成人欧美视频| 欧美精品一区二区免费开放| 99精国产麻豆久久婷婷| 亚洲成av片中文字幕在线观看| 亚洲第一av免费看| 每晚都被弄得嗷嗷叫到高潮| 国产精品1区2区在线观看.| 搡老熟女国产l中国老女人| xxxhd国产人妻xxx| 欧美中文综合在线视频| 老司机午夜十八禁免费视频| 十八禁人妻一区二区| 精品久久久久久久毛片微露脸| 亚洲激情在线av| 亚洲自拍偷在线| 老汉色∧v一级毛片| 老司机深夜福利视频在线观看| 国产在线精品亚洲第一网站| 国产成人欧美在线观看| 高清毛片免费观看视频网站 | 久久精品成人免费网站| 国产有黄有色有爽视频| 欧美激情极品国产一区二区三区| 两个人免费观看高清视频| 亚洲欧美日韩无卡精品| 欧美日韩亚洲高清精品| 国产精品综合久久久久久久免费 | 欧美黄色淫秽网站| 亚洲精品成人av观看孕妇| 免费在线观看视频国产中文字幕亚洲| 精品福利观看| 亚洲五月色婷婷综合| 少妇 在线观看| 99久久久亚洲精品蜜臀av| 自线自在国产av| 成人三级黄色视频| 午夜免费观看网址| 美女高潮到喷水免费观看| 精品国产乱码久久久久久男人| 中国美女看黄片| 中文字幕人妻熟女乱码| 亚洲成国产人片在线观看| 国内久久婷婷六月综合欲色啪| 可以在线观看毛片的网站| 欧美 亚洲 国产 日韩一| 乱人伦中国视频| 国产xxxxx性猛交| 久久中文字幕人妻熟女| 欧美中文日本在线观看视频| 欧美日韩精品网址| 岛国视频午夜一区免费看| 后天国语完整版免费观看| 亚洲精品久久午夜乱码| 亚洲美女黄片视频| 久久99一区二区三区| 成人18禁在线播放| 999久久久国产精品视频| 国产成人精品在线电影| 国产激情久久老熟女| 午夜激情av网站| av视频免费观看在线观看| 丝袜人妻中文字幕| 国产精品野战在线观看 | 色综合欧美亚洲国产小说| 欧美激情久久久久久爽电影 | 亚洲国产欧美一区二区综合| 怎么达到女性高潮| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 好看av亚洲va欧美ⅴa在| 精品久久久久久,| 满18在线观看网站| av天堂久久9| 热re99久久精品国产66热6| 国产黄a三级三级三级人| 久久久精品国产亚洲av高清涩受| 午夜激情av网站| 欧美激情 高清一区二区三区| 国产av又大| 91精品国产国语对白视频| 国产野战对白在线观看| 嫁个100分男人电影在线观看| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 一个人免费在线观看的高清视频| 可以免费在线观看a视频的电影网站| 亚洲一区中文字幕在线| 岛国视频午夜一区免费看| 欧美日韩福利视频一区二区| 9热在线视频观看99| 国产91精品成人一区二区三区| 午夜免费激情av| 啦啦啦在线免费观看视频4| 国产精品影院久久| 高清黄色对白视频在线免费看| www.熟女人妻精品国产| 欧美一区二区精品小视频在线| 亚洲狠狠婷婷综合久久图片| 免费av中文字幕在线| 久久久久久免费高清国产稀缺| 国产三级黄色录像| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 交换朋友夫妻互换小说| 精品国产乱子伦一区二区三区| 成人av一区二区三区在线看| 一夜夜www| 国产人伦9x9x在线观看| 国产男靠女视频免费网站| 亚洲av熟女| 久久午夜综合久久蜜桃| 伊人久久大香线蕉亚洲五| 少妇裸体淫交视频免费看高清 | 亚洲自拍偷在线| 人成视频在线观看免费观看| 亚洲情色 制服丝袜| www.自偷自拍.com| 亚洲精品中文字幕一二三四区| 在线播放国产精品三级| www日本在线高清视频| 91麻豆精品激情在线观看国产 | 亚洲中文字幕日韩| 日本精品一区二区三区蜜桃| 国产精品一区二区三区四区久久 | 757午夜福利合集在线观看| 大型黄色视频在线免费观看| 俄罗斯特黄特色一大片| 熟女少妇亚洲综合色aaa.| 激情视频va一区二区三区| 黄色女人牲交| 亚洲欧美激情综合另类| 午夜两性在线视频| 免费女性裸体啪啪无遮挡网站| 国产97色在线日韩免费| 一区二区三区国产精品乱码| 一区二区三区激情视频| 99精品欧美一区二区三区四区| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看 | 国产成人欧美| 日韩欧美一区视频在线观看| 国产麻豆69| 黑人操中国人逼视频| av在线播放免费不卡| 久久99一区二区三区| 久久久久久久久免费视频了| 身体一侧抽搐| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 久久久久久久久久久久大奶| 欧美最黄视频在线播放免费 | 性色av乱码一区二区三区2| 黄色片一级片一级黄色片| 久久影院123| 在线观看舔阴道视频| 老司机靠b影院| 午夜免费观看网址| 十八禁网站免费在线| 叶爱在线成人免费视频播放| 亚洲欧美一区二区三区黑人| 中文字幕色久视频| 怎么达到女性高潮| 欧美性长视频在线观看| 精品国产亚洲在线| 天天躁夜夜躁狠狠躁躁| 国产免费av片在线观看野外av| 丝袜美腿诱惑在线| 欧美成人性av电影在线观看| 国产极品粉嫩免费观看在线| 99久久精品国产亚洲精品| 最新在线观看一区二区三区| 可以免费在线观看a视频的电影网站| 在线av久久热| 淫妇啪啪啪对白视频| 午夜福利影视在线免费观看| 精品少妇一区二区三区视频日本电影| 亚洲精品一区av在线观看| 亚洲久久久国产精品| a级毛片在线看网站| 国产欧美日韩综合在线一区二区| 成人精品一区二区免费| 老司机靠b影院| 岛国视频午夜一区免费看| 亚洲欧美一区二区三区黑人| 国产亚洲欧美精品永久| 少妇 在线观看| 深夜精品福利| 在线观看日韩欧美| 水蜜桃什么品种好| 国产蜜桃级精品一区二区三区| 中文欧美无线码| 两个人看的免费小视频| 亚洲专区中文字幕在线| 大陆偷拍与自拍| av片东京热男人的天堂| 悠悠久久av| 国产亚洲欧美在线一区二区| 国产精品美女特级片免费视频播放器 | 又黄又粗又硬又大视频| 日韩欧美一区二区三区在线观看| 免费在线观看黄色视频的| a级毛片黄视频| 久久精品aⅴ一区二区三区四区| a级毛片在线看网站| 日韩欧美在线二视频| 侵犯人妻中文字幕一二三四区| 午夜福利免费观看在线| 怎么达到女性高潮| 黄色女人牲交| 久久草成人影院| 国产精品一区二区精品视频观看| 日本精品一区二区三区蜜桃| 婷婷丁香在线五月| 国产乱人伦免费视频| 另类亚洲欧美激情| 中文字幕人妻熟女乱码| 亚洲va日本ⅴa欧美va伊人久久| 长腿黑丝高跟| а√天堂www在线а√下载| 国内久久婷婷六月综合欲色啪| 免费看a级黄色片| 男人操女人黄网站| 亚洲人成网站在线播放欧美日韩| 一夜夜www| 久久影院123| 亚洲精品中文字幕在线视频| 亚洲一区二区三区欧美精品| 精品国内亚洲2022精品成人| 在线视频色国产色| 一个人观看的视频www高清免费观看 | 大陆偷拍与自拍| 国产精品一区二区免费欧美| 99久久久亚洲精品蜜臀av| av有码第一页| 久久久久久久久免费视频了| 日韩成人在线观看一区二区三区| av超薄肉色丝袜交足视频| 制服人妻中文乱码| 精品一区二区三区视频在线观看免费 | 久久久久久久久久久久大奶| 久久青草综合色| 波多野结衣高清无吗| 日本a在线网址| 国产激情久久老熟女| 老汉色∧v一级毛片| 久久天躁狠狠躁夜夜2o2o| 50天的宝宝边吃奶边哭怎么回事| 国产亚洲精品第一综合不卡| av超薄肉色丝袜交足视频| 黄色a级毛片大全视频| 欧美黑人精品巨大| 成人黄色视频免费在线看| 久久精品亚洲精品国产色婷小说| 欧美黄色片欧美黄色片| aaaaa片日本免费| 女人被躁到高潮嗷嗷叫费观| 久久99一区二区三区| 国产日韩一区二区三区精品不卡| 丁香欧美五月| 欧美黄色淫秽网站| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久毛片微露脸| 老司机福利观看| 一区二区三区激情视频| 在线观看免费日韩欧美大片| 久久久久久久久中文| 在线视频色国产色| 欧美色视频一区免费| 国产精品野战在线观看 | 欧美一级毛片孕妇| 国产99白浆流出| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 村上凉子中文字幕在线| 制服诱惑二区| 岛国在线观看网站| 男女之事视频高清在线观看| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 日韩人妻精品一区2区三区| 国产激情欧美一区二区| 一区二区日韩欧美中文字幕| 午夜免费激情av| 麻豆一二三区av精品| 黄色女人牲交| 美女扒开内裤让男人捅视频| 身体一侧抽搐| tocl精华| 国产区一区二久久| 高潮久久久久久久久久久不卡| 欧美亚洲日本最大视频资源| 日本a在线网址| 精品久久久久久久久久免费视频 | 国产成人一区二区三区免费视频网站| 丰满迷人的少妇在线观看| 激情视频va一区二区三区| 国产精品免费一区二区三区在线| 视频在线观看一区二区三区| 中文字幕色久视频| 欧美日韩瑟瑟在线播放| 美女高潮到喷水免费观看| 在线播放国产精品三级| 亚洲全国av大片| 嫩草影院精品99| 色婷婷av一区二区三区视频| 窝窝影院91人妻| 精品一区二区三卡| 国产精品美女特级片免费视频播放器 | 成人av一区二区三区在线看| 美女高潮喷水抽搐中文字幕| 中文字幕人妻丝袜一区二区| 精品久久久久久电影网| 亚洲全国av大片| a级毛片在线看网站| 一级毛片精品| 成人免费观看视频高清| 午夜成年电影在线免费观看| 久久久国产欧美日韩av| а√天堂www在线а√下载| 夜夜躁狠狠躁天天躁| 波多野结衣av一区二区av| 国产成人影院久久av| 亚洲自偷自拍图片 自拍| 欧美日韩乱码在线| 黄片播放在线免费| 亚洲片人在线观看| 真人一进一出gif抽搐免费| 亚洲精品国产区一区二| a级毛片在线看网站| 国产99久久九九免费精品| 精品人妻在线不人妻| 国产精品 欧美亚洲| 波多野结衣高清无吗| 99riav亚洲国产免费| 99国产精品一区二区三区| 亚洲欧美一区二区三区黑人| 黑人操中国人逼视频| 天堂动漫精品| 99国产精品99久久久久| 长腿黑丝高跟| 黄色丝袜av网址大全| 久久天堂一区二区三区四区| 亚洲avbb在线观看| 久久精品91蜜桃| 免费观看人在逋| 国产精品久久久久成人av| 国产精品一区二区三区四区久久 | 免费人成视频x8x8入口观看| 亚洲色图av天堂| 黄色视频不卡| 久久久久久久久中文| 88av欧美| 精品日产1卡2卡| 久久久国产欧美日韩av| bbb黄色大片| 午夜激情av网站| 少妇裸体淫交视频免费看高清 | av电影中文网址| 叶爱在线成人免费视频播放| 淫妇啪啪啪对白视频| avwww免费| 亚洲免费av在线视频| 亚洲成人精品中文字幕电影 | 久久热在线av| 久久精品aⅴ一区二区三区四区| 欧美色视频一区免费| 侵犯人妻中文字幕一二三四区| 超色免费av| 亚洲精品中文字幕在线视频| 精品福利永久在线观看| 女同久久另类99精品国产91| 高清黄色对白视频在线免费看| 在线观看舔阴道视频| 91av网站免费观看| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| 可以免费在线观看a视频的电影网站| 麻豆久久精品国产亚洲av | 久久草成人影院| 国产亚洲av高清不卡| 日本免费a在线| 亚洲色图综合在线观看| 国产成年人精品一区二区 | 国产又色又爽无遮挡免费看| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| av福利片在线| 在线观看免费视频日本深夜| 亚洲av美国av| 成人三级做爰电影| www.自偷自拍.com| 精品日产1卡2卡| 免费在线观看影片大全网站| av网站免费在线观看视频| 国产亚洲精品一区二区www| 高清黄色对白视频在线免费看| 丁香六月欧美| 成在线人永久免费视频| 美女国产高潮福利片在线看| a在线观看视频网站| 12—13女人毛片做爰片一| 后天国语完整版免费观看| 成人国产一区最新在线观看| 在线观看午夜福利视频| 国产精品国产高清国产av| 久久精品91无色码中文字幕| 国产精品久久久久成人av| 亚洲精品美女久久av网站| 日本撒尿小便嘘嘘汇集6| 啦啦啦在线免费观看视频4| 黄色 视频免费看| 国产精品 国内视频| 一级毛片女人18水好多| 1024香蕉在线观看| 日韩人妻精品一区2区三区| 老司机福利观看| 少妇 在线观看| 午夜免费激情av| 天天躁夜夜躁狠狠躁躁| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产美女av久久久久小说| xxxhd国产人妻xxx| 极品教师在线免费播放| 中文字幕另类日韩欧美亚洲嫩草| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 亚洲久久久国产精品| 人妻久久中文字幕网| 久久人人精品亚洲av| 大陆偷拍与自拍| 中出人妻视频一区二区| 男人舔女人的私密视频| 香蕉国产在线看| 成年版毛片免费区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美成人免费av一区二区三区| 中文字幕精品免费在线观看视频| 亚洲精品久久成人aⅴ小说| 在线观看日韩欧美| 欧美日韩亚洲国产一区二区在线观看| 亚洲av成人一区二区三| 日韩精品青青久久久久久| 亚洲精品美女久久av网站| 色综合欧美亚洲国产小说| 999久久久国产精品视频| 午夜福利免费观看在线| 色尼玛亚洲综合影院| 在线观看66精品国产| 久久久久久久午夜电影 | av天堂在线播放| 精品欧美一区二区三区在线| 亚洲第一青青草原| netflix在线观看网站| 国产黄a三级三级三级人| 久久影院123| 午夜福利在线免费观看网站| 亚洲欧美一区二区三区久久| 色综合站精品国产| 亚洲一区中文字幕在线| 国产一区二区三区在线臀色熟女 | 啦啦啦 在线观看视频| 香蕉久久夜色| 淫秽高清视频在线观看| 香蕉久久夜色| 高清毛片免费观看视频网站 | 国产精品免费一区二区三区在线| 中出人妻视频一区二区| 人人妻人人爽人人添夜夜欢视频| 亚洲在线自拍视频| 欧美日韩亚洲综合一区二区三区_| 欧美日韩精品网址| www.www免费av| 波多野结衣高清无吗| 国产成人一区二区三区免费视频网站| 精品福利观看| 国产精华一区二区三区| 国产亚洲精品久久久久久毛片| 91国产中文字幕| 久久人人精品亚洲av| 中出人妻视频一区二区| 十八禁网站免费在线| 中文字幕精品免费在线观看视频| 大香蕉久久成人网| 一级片'在线观看视频| 搡老岳熟女国产| 桃色一区二区三区在线观看| av在线播放免费不卡| av国产精品久久久久影院| 在线免费观看的www视频| 色播在线永久视频| 免费高清在线观看日韩| 看免费av毛片| 亚洲精品久久成人aⅴ小说| 法律面前人人平等表现在哪些方面| 久久精品影院6| 窝窝影院91人妻| 欧美成人性av电影在线观看| 一级黄色大片毛片| 久久99一区二区三区| 久久久久久免费高清国产稀缺| 一边摸一边抽搐一进一小说| 国产精品一区二区精品视频观看| 99久久国产精品久久久| 黑人巨大精品欧美一区二区蜜桃| 欧美黄色淫秽网站| 亚洲男人天堂网一区| 国产精品一区二区免费欧美| 99国产精品免费福利视频| 欧美人与性动交α欧美软件| 一个人观看的视频www高清免费观看 | av有码第一页| 在线观看免费高清a一片| 亚洲精品粉嫩美女一区| 美女国产高潮福利片在线看| 十八禁人妻一区二区| 久久久久久免费高清国产稀缺| 十分钟在线观看高清视频www| 搡老岳熟女国产| 99精品在免费线老司机午夜| 成年人免费黄色播放视频| 99国产精品99久久久久| 久久久水蜜桃国产精品网| 久久久精品国产亚洲av高清涩受| 日本五十路高清| 亚洲成国产人片在线观看| 亚洲精品美女久久久久99蜜臀| 怎么达到女性高潮| 亚洲人成77777在线视频| 一二三四在线观看免费中文在| 一区二区三区精品91| 动漫黄色视频在线观看| 精品日产1卡2卡| 久久精品国产亚洲av香蕉五月| 精品国产国语对白av| 90打野战视频偷拍视频| 91在线观看av| 久久九九热精品免费| 色哟哟哟哟哟哟| 国产区一区二久久| av天堂在线播放| 黄片播放在线免费| 90打野战视频偷拍视频| av电影中文网址| 99在线视频只有这里精品首页| 999久久久国产精品视频| 正在播放国产对白刺激| 亚洲熟妇中文字幕五十中出 | 黄片小视频在线播放| 搡老熟女国产l中国老女人| 免费人成视频x8x8入口观看| 一区二区三区激情视频| 国产av又大| 日韩欧美国产一区二区入口| 9色porny在线观看| 国产精品爽爽va在线观看网站 | 亚洲avbb在线观看| 在线观看一区二区三区| 91在线观看av| 一个人免费在线观看的高清视频| 黄色毛片三级朝国网站| 日日干狠狠操夜夜爽| 国产欧美日韩精品亚洲av| 狂野欧美激情性xxxx| 久久狼人影院| 久久久久久久久中文| 国产精品久久久久久人妻精品电影| 亚洲伊人色综图| 在线av久久热| 在线天堂中文资源库| 国产精品乱码一区二三区的特点 | 久久久久久亚洲精品国产蜜桃av| 亚洲伊人色综图| 亚洲欧美精品综合久久99| 国产99白浆流出| 久久久精品国产亚洲av高清涩受| 欧美黄色片欧美黄色片| 手机成人av网站| 亚洲人成伊人成综合网2020|