• <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ì)天童常綠闊葉林凋落物量的影響
    窝窝影院91人妻| 男女做爰动态图高潮gif福利片| 1024香蕉在线观看| 久久精品91无色码中文字幕| www.熟女人妻精品国产| 搞女人的毛片| 男人舔女人的私密视频| 欧美黄色片欧美黄色片| 高清毛片免费观看视频网站| 精品久久久久久成人av| 婷婷精品国产亚洲av| www日本黄色视频网| 欧美一区二区国产精品久久精品| 亚洲国产色片| 亚洲欧美日韩高清在线视频| 日本黄色视频三级网站网址| 日韩欧美 国产精品| 丰满的人妻完整版| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 两人在一起打扑克的视频| 麻豆一二三区av精品| 亚洲熟女毛片儿| 久久99热这里只有精品18| 日本免费a在线| 非洲黑人性xxxx精品又粗又长| 成人欧美大片| 每晚都被弄得嗷嗷叫到高潮| 久久久久久国产a免费观看| 99在线视频只有这里精品首页| 久久人妻av系列| www日本在线高清视频| 精品一区二区三区视频在线观看免费| 99精品久久久久人妻精品| 这个男人来自地球电影免费观看| 男人和女人高潮做爰伦理| 法律面前人人平等表现在哪些方面| 757午夜福利合集在线观看| 亚洲欧美日韩无卡精品| 欧美激情在线99| 久久久国产成人精品二区| 国产亚洲精品av在线| 97碰自拍视频| 亚洲狠狠婷婷综合久久图片| 少妇的丰满在线观看| 无限看片的www在线观看| 午夜福利在线在线| 亚洲七黄色美女视频| 色视频www国产| 两个人看的免费小视频| 精品久久久久久久人妻蜜臀av| 又黄又爽又免费观看的视频| 亚洲专区中文字幕在线| 九九热线精品视视频播放| 国产精品美女特级片免费视频播放器 | 日韩欧美在线二视频| 国内精品一区二区在线观看| 一进一出好大好爽视频| 又大又爽又粗| 亚洲av成人精品一区久久| 精品国产超薄肉色丝袜足j| 99热只有精品国产| 禁无遮挡网站| 精品一区二区三区av网在线观看| 日韩免费av在线播放| 床上黄色一级片| 精品熟女少妇八av免费久了| 国产成人精品久久二区二区免费| 国产三级中文精品| 亚洲七黄色美女视频| 成人亚洲精品av一区二区| 国产精品一区二区三区四区免费观看 | 国产高清三级在线| 亚洲 欧美 日韩 在线 免费| 国产精品日韩av在线免费观看| 亚洲人与动物交配视频| 在线国产一区二区在线| 五月伊人婷婷丁香| 国产成人精品无人区| 成人亚洲精品av一区二区| 成年版毛片免费区| 中文字幕人妻丝袜一区二区| 午夜激情福利司机影院| 亚洲va日本ⅴa欧美va伊人久久| 97人妻精品一区二区三区麻豆| 18禁美女被吸乳视频| 婷婷丁香在线五月| 熟女电影av网| 色综合站精品国产| 亚洲精品456在线播放app | 欧美成狂野欧美在线观看| 高潮久久久久久久久久久不卡| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 无人区码免费观看不卡| 国产亚洲精品av在线| 成人三级黄色视频| 又大又爽又粗| 亚洲 欧美一区二区三区| 男人舔女人的私密视频| 观看美女的网站| 欧美极品一区二区三区四区| 亚洲国产高清在线一区二区三| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 不卡av一区二区三区| avwww免费| 国产精品 国内视频| 最新美女视频免费是黄的| 成年女人看的毛片在线观看| 精品久久久久久久毛片微露脸| 国产日本99.免费观看| 免费看光身美女| 国产主播在线观看一区二区| 国产成人aa在线观看| 亚洲国产日韩欧美精品在线观看 | 两个人视频免费观看高清| 国产成人啪精品午夜网站| 黑人巨大精品欧美一区二区mp4| 色吧在线观看| 国产精品美女特级片免费视频播放器 | 国产精品,欧美在线| 午夜精品在线福利| 亚洲专区字幕在线| 1024香蕉在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲最大成人中文| 免费高清视频大片| 欧美3d第一页| 日韩中文字幕欧美一区二区| 国产精品九九99| av片东京热男人的天堂| 精品不卡国产一区二区三区| 中文字幕熟女人妻在线| 日韩欧美在线二视频| 亚洲精品在线美女| 午夜亚洲福利在线播放| 国产在线精品亚洲第一网站| 少妇人妻一区二区三区视频| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 岛国视频午夜一区免费看| 黄色片一级片一级黄色片| 国产伦一二天堂av在线观看| 久久亚洲精品不卡| 国产美女午夜福利| 亚洲专区国产一区二区| 淫秽高清视频在线观看| 国产高清videossex| 麻豆久久精品国产亚洲av| 国产精品爽爽va在线观看网站| 国产69精品久久久久777片 | 成年免费大片在线观看| 欧美三级亚洲精品| 国产一区二区三区在线臀色熟女| 一卡2卡三卡四卡精品乱码亚洲| 精品不卡国产一区二区三区| 我要搜黄色片| 国产精品一区二区免费欧美| 成人三级黄色视频| 嫁个100分男人电影在线观看| 国内精品久久久久精免费| 欧美一区二区精品小视频在线| 88av欧美| 日本免费a在线| 国产成人影院久久av| 国产亚洲欧美98| 亚洲 欧美一区二区三区| 夜夜看夜夜爽夜夜摸| 日本 av在线| 一个人免费在线观看电影 | 国产男靠女视频免费网站| 久久精品综合一区二区三区| 1024手机看黄色片| 中文在线观看免费www的网站| 日韩精品中文字幕看吧| 国产精品久久久av美女十八| 日本a在线网址| 亚洲国产精品久久男人天堂| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| 亚洲午夜理论影院| 深夜精品福利| 国产97色在线日韩免费| 999久久久精品免费观看国产| 在线十欧美十亚洲十日本专区| 热99在线观看视频| 美女 人体艺术 gogo| 国产精品国产高清国产av| 国产精品,欧美在线| 久久久久性生活片| 网址你懂的国产日韩在线| 亚洲精品粉嫩美女一区| 久久久久精品国产欧美久久久| 国内少妇人妻偷人精品xxx网站 | 日本免费一区二区三区高清不卡| 嫩草影院入口| 成人精品一区二区免费| 美女黄网站色视频| 一卡2卡三卡四卡精品乱码亚洲| av在线蜜桃| 国产黄色小视频在线观看| 淫妇啪啪啪对白视频| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 小说图片视频综合网站| 一夜夜www| 日韩大尺度精品在线看网址| 久久久成人免费电影| 国产成人aa在线观看| 国产黄片美女视频| 日韩欧美免费精品| 中文在线观看免费www的网站| 婷婷精品国产亚洲av在线| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 一二三四社区在线视频社区8| 久久久国产成人精品二区| 国产午夜福利久久久久久| 观看美女的网站| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 久久久色成人| 免费av不卡在线播放| 欧美性猛交黑人性爽| svipshipincom国产片| 91字幕亚洲| 国产精品 国内视频| 中文字幕高清在线视频| 岛国视频午夜一区免费看| 亚洲精品456在线播放app | 国产单亲对白刺激| 欧美又色又爽又黄视频| 午夜精品久久久久久毛片777| 757午夜福利合集在线观看| 麻豆国产97在线/欧美| 欧美中文综合在线视频| 在线视频色国产色| 久久久久性生活片| 男人舔女人的私密视频| 成年人黄色毛片网站| 亚洲欧美日韩东京热| 国产欧美日韩精品亚洲av| 男女下面进入的视频免费午夜| 久久午夜亚洲精品久久| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 亚洲午夜理论影院| 国产精品日韩av在线免费观看| 国产精品 欧美亚洲| 神马国产精品三级电影在线观看| 99精品久久久久人妻精品| 黄色日韩在线| 成人午夜高清在线视频| 国产精品99久久久久久久久| 国产一区二区三区在线臀色熟女| 美女被艹到高潮喷水动态| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 九九在线视频观看精品| 中亚洲国语对白在线视频| 国产亚洲精品一区二区www| 天堂动漫精品| 级片在线观看| 国内揄拍国产精品人妻在线| 欧美一级毛片孕妇| 国产乱人视频| 亚洲国产色片| 色哟哟哟哟哟哟| 老鸭窝网址在线观看| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 成人三级做爰电影| 无遮挡黄片免费观看| 在线观看日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 国产精品99久久99久久久不卡| 叶爱在线成人免费视频播放| 看黄色毛片网站| 1024香蕉在线观看| 级片在线观看| 他把我摸到了高潮在线观看| svipshipincom国产片| 久久国产乱子伦精品免费另类| 成人性生交大片免费视频hd| av在线蜜桃| 丝袜人妻中文字幕| 亚洲专区国产一区二区| 国产精品 国内视频| 欧美绝顶高潮抽搐喷水| 宅男免费午夜| 亚洲黑人精品在线| 日本免费a在线| 日本撒尿小便嘘嘘汇集6| 国产av麻豆久久久久久久| netflix在线观看网站| 亚洲乱码一区二区免费版| 色精品久久人妻99蜜桃| 黄色丝袜av网址大全| 亚洲熟女毛片儿| 久久久久久大精品| 亚洲av电影不卡..在线观看| 欧美三级亚洲精品| 听说在线观看完整版免费高清| 色av中文字幕| 国产一区二区三区视频了| 亚洲性夜色夜夜综合| 蜜桃久久精品国产亚洲av| 久久久久久人人人人人| 舔av片在线| 精品一区二区三区四区五区乱码| 亚洲av片天天在线观看| 成人三级黄色视频| 免费看十八禁软件| 最新美女视频免费是黄的| 欧美三级亚洲精品| 听说在线观看完整版免费高清| 一本精品99久久精品77| 久久久久久九九精品二区国产| 亚洲欧美一区二区三区黑人| 级片在线观看| 伊人久久大香线蕉亚洲五| 99国产精品一区二区蜜桃av| 国产又黄又爽又无遮挡在线| 男人舔女人的私密视频| 亚洲七黄色美女视频| 国产在线精品亚洲第一网站| 亚洲精品一区av在线观看| 国产精品久久久久久人妻精品电影| 99热精品在线国产| 神马国产精品三级电影在线观看| 一本久久中文字幕| 久久久久免费精品人妻一区二区| 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器 | АⅤ资源中文在线天堂| 精品熟女少妇八av免费久了| 可以在线观看的亚洲视频| 特级一级黄色大片| 色吧在线观看| 色视频www国产| 亚洲 欧美 日韩 在线 免费| 国产精品爽爽va在线观看网站| 国产激情久久老熟女| 每晚都被弄得嗷嗷叫到高潮| 日本a在线网址| 精品国内亚洲2022精品成人| av欧美777| 美女 人体艺术 gogo| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| www国产在线视频色| 午夜福利在线观看免费完整高清在 | 成人国产一区最新在线观看| a级毛片在线看网站| 十八禁人妻一区二区| 国产97色在线日韩免费| 国产成人精品久久二区二区免费| 香蕉av资源在线| 97超级碰碰碰精品色视频在线观看| 亚洲色图av天堂| 亚洲av成人av| 国产日本99.免费观看| 老熟妇乱子伦视频在线观看| 亚洲精品国产精品久久久不卡| 91麻豆av在线| 人妻夜夜爽99麻豆av| 两个人视频免费观看高清| 脱女人内裤的视频| av欧美777| 精品无人区乱码1区二区| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 黑人欧美特级aaaaaa片| 亚洲成人精品中文字幕电影| 成年人黄色毛片网站| 欧美日韩瑟瑟在线播放| 久久香蕉国产精品| 日本在线视频免费播放| 免费看a级黄色片| 毛片女人毛片| 久久久久国产精品人妻aⅴ院| 国产一区二区激情短视频| 国产av一区在线观看免费| 人人妻人人澡欧美一区二区| 久久午夜亚洲精品久久| 最近最新免费中文字幕在线| 99久久精品一区二区三区| 少妇人妻一区二区三区视频| 两性夫妻黄色片| 国产精品乱码一区二三区的特点| 女生性感内裤真人,穿戴方法视频| 日本黄色片子视频| 12—13女人毛片做爰片一| 久久精品91无色码中文字幕| 两个人看的免费小视频| 18禁黄网站禁片免费观看直播| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 中亚洲国语对白在线视频| 在线a可以看的网站| 最近最新免费中文字幕在线| 国产99白浆流出| 最近在线观看免费完整版| 大型黄色视频在线免费观看| 欧美丝袜亚洲另类 | 99热6这里只有精品| 色在线成人网| 午夜福利成人在线免费观看| 国产精品久久久人人做人人爽| 真实男女啪啪啪动态图| 久久精品综合一区二区三区| 亚洲专区字幕在线| 夜夜爽天天搞| 亚洲一区二区三区色噜噜| 在线观看免费午夜福利视频| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 国产蜜桃级精品一区二区三区| 男女那种视频在线观看| 麻豆国产97在线/欧美| 亚洲av成人一区二区三| 最近最新中文字幕大全免费视频| 国产乱人伦免费视频| a在线观看视频网站| 欧美日韩瑟瑟在线播放| 精品国产乱码久久久久久男人| 久久精品亚洲精品国产色婷小说| 极品教师在线免费播放| 毛片女人毛片| 韩国av一区二区三区四区| 国产高清视频在线观看网站| 亚洲av电影不卡..在线观看| 人人妻,人人澡人人爽秒播| av女优亚洲男人天堂 | 国产成年人精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲自拍偷在线| 亚洲人与动物交配视频| 三级国产精品欧美在线观看 | 麻豆久久精品国产亚洲av| 国产伦在线观看视频一区| 成人精品一区二区免费| 黄片小视频在线播放| 欧美大码av| 一a级毛片在线观看| 成熟少妇高潮喷水视频| 久久欧美精品欧美久久欧美| 麻豆国产av国片精品| 99久久久亚洲精品蜜臀av| 欧美激情在线99| 亚洲午夜理论影院| 国产欧美日韩精品亚洲av| 国产欧美日韩精品一区二区| 欧美黄色片欧美黄色片| 观看免费一级毛片| 国产成年人精品一区二区| 黑人欧美特级aaaaaa片| 亚洲,欧美精品.| 亚洲欧美日韩高清在线视频| 亚洲熟妇熟女久久| 亚洲在线观看片| 一区福利在线观看| av女优亚洲男人天堂 | 淫秽高清视频在线观看| 亚洲avbb在线观看| 亚洲成人精品中文字幕电影| 久久这里只有精品19| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| 欧美av亚洲av综合av国产av| 国产v大片淫在线免费观看| avwww免费| 成年女人毛片免费观看观看9| 亚洲真实伦在线观看| 99精品久久久久人妻精品| 一区二区三区激情视频| 18美女黄网站色大片免费观看| 国产精品98久久久久久宅男小说| 丝袜人妻中文字幕| 999久久久精品免费观看国产| av女优亚洲男人天堂 | 国产成人欧美在线观看| 少妇人妻一区二区三区视频| 校园春色视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 很黄的视频免费| 又粗又爽又猛毛片免费看| 97超级碰碰碰精品色视频在线观看| 成人特级黄色片久久久久久久| 嫁个100分男人电影在线观看| 午夜两性在线视频| 国模一区二区三区四区视频 | 国产精品日韩av在线免费观看| 可以在线观看毛片的网站| 日本精品一区二区三区蜜桃| 丁香六月欧美| 叶爱在线成人免费视频播放| 久久精品亚洲精品国产色婷小说| 久久久久久大精品| 麻豆一二三区av精品| 亚洲精品美女久久久久99蜜臀| 一本精品99久久精品77| 国产精品免费一区二区三区在线| 天堂影院成人在线观看| 午夜免费观看网址| 制服人妻中文乱码| 久久中文看片网| 国产综合懂色| 国产私拍福利视频在线观看| 国产aⅴ精品一区二区三区波| 久久天堂一区二区三区四区| svipshipincom国产片| 国产伦精品一区二区三区视频9 | 午夜精品久久久久久毛片777| 99久久久亚洲精品蜜臀av| 精品一区二区三区视频在线 | 亚洲成人精品中文字幕电影| 在线免费观看不下载黄p国产 | 国产又黄又爽又无遮挡在线| 国产精品精品国产色婷婷| 小说图片视频综合网站| 在线看三级毛片| 久久精品国产清高在天天线| 欧美乱码精品一区二区三区| 中国美女看黄片| 久久久久亚洲av毛片大全| 欧美国产日韩亚洲一区| 午夜福利免费观看在线| 日本成人三级电影网站| 桃色一区二区三区在线观看| 亚洲,欧美精品.| 欧美乱色亚洲激情| 后天国语完整版免费观看| 美女免费视频网站| 亚洲 欧美一区二区三区| 99国产精品一区二区蜜桃av| 18禁观看日本| 在线播放国产精品三级| 极品教师在线免费播放| 精品一区二区三区视频在线观看免费| 国产单亲对白刺激| АⅤ资源中文在线天堂| 亚洲色图av天堂| 午夜激情欧美在线| 黄色女人牲交| netflix在线观看网站| 高清毛片免费观看视频网站| 日韩成人在线观看一区二区三区| 久久久精品大字幕| 国产成人精品无人区| 999久久久国产精品视频| 精品99又大又爽又粗少妇毛片 | 一个人看的www免费观看视频| 99久久国产精品久久久| 欧美午夜高清在线| 国产午夜精品论理片| 午夜精品在线福利| 在线观看一区二区三区| 久久久久久久久中文| 99riav亚洲国产免费| 精品久久久久久,| 国产aⅴ精品一区二区三区波| 国产午夜精品久久久久久| 黑人操中国人逼视频| 两个人视频免费观看高清| 国产视频一区二区在线看| 精品欧美国产一区二区三| 精品电影一区二区在线| 99热6这里只有精品| 黄色女人牲交| 亚洲专区字幕在线| 精品一区二区三区视频在线 | 又大又爽又粗| 成人三级黄色视频| 天堂av国产一区二区熟女人妻| 国产精品 国内视频| 人人妻人人澡欧美一区二区| 国内精品久久久久精免费| 免费搜索国产男女视频| 日本五十路高清| 午夜精品一区二区三区免费看| 日韩av在线大香蕉| 此物有八面人人有两片| 亚洲人成网站高清观看| 免费人成视频x8x8入口观看| 97碰自拍视频| 国产精品一区二区免费欧美| 欧美激情在线99| 大型黄色视频在线免费观看| 国产aⅴ精品一区二区三区波| а√天堂www在线а√下载| 成人国产一区最新在线观看| 久久热在线av| 国产成人精品久久二区二区91| 亚洲av成人不卡在线观看播放网| 久久热在线av| 免费人成视频x8x8入口观看| av中文乱码字幕在线| 亚洲在线观看片| 日本 欧美在线| 亚洲美女黄片视频| 国产欧美日韩精品一区二区| 性色avwww在线观看| 此物有八面人人有两片| 成熟少妇高潮喷水视频| 久久久久久久久久黄片| 婷婷丁香在线五月| 亚洲最大成人中文| 午夜福利在线观看吧| 精品一区二区三区视频在线观看免费| 精品一区二区三区视频在线 |