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

    GF-3全極化SAR數(shù)據(jù)極化分解估算人工林冠層生物量*

    2020-10-17 04:58:24魏晶昱范文義毛學(xué)剛
    林業(yè)科學(xué) 2020年9期
    關(guān)鍵詞:冠層樣地極化

    魏晶昱 范文義 于 穎 毛學(xué)剛

    (東北林業(yè)大學(xué)林學(xué)院 森林生態(tài)系統(tǒng)可持續(xù)經(jīng)營(yíng)教育部重點(diǎn)實(shí)驗(yàn)室 哈爾濱 150040)

    森林生物量指單位面積森林群落在一定時(shí)間內(nèi)積累的有機(jī)質(zhì)總量(尹惠妍等, 2014),其可直接反映森林植被固碳現(xiàn)狀,是評(píng)價(jià)森林生態(tài)系統(tǒng)結(jié)構(gòu)、功能和生產(chǎn)力的重要指標(biāo)之一,也是進(jìn)行森林資源清查、規(guī)劃和管理的依據(jù)。樹(shù)冠是樹(shù)木光合作用、呼吸作用、蒸騰作用的主要發(fā)生場(chǎng)所,影響樹(shù)木生長(zhǎng)和變化(沈浩等, 2017; Ozanneetal., 2003),冠層生物量是森林生物量研究的重要內(nèi)容,其對(duì)環(huán)境變化十分敏感,在一定程度上影響著森林生態(tài)系統(tǒng)生產(chǎn)力(李鳳日, 2004; 劉兆剛等, 2005; Wangetal., 1990; Kramer, 1996; Tahvanainenetal., 2008)。

    傳統(tǒng)的冠層生物量估算通常采用平均枝法或標(biāo)準(zhǔn)枝法(Ibrahim, 1995),存在測(cè)定困難、耗時(shí)耗力等問(wèn)題,且因葉片水分損失較快也會(huì)帶來(lái)難以估計(jì)的誤差(劉琪璟, 2009; Sonetal., 2001)。隨著遙感技術(shù)不斷發(fā)展,快速、準(zhǔn)確、大范圍地估算冠層生物量成為可能,傳統(tǒng)光學(xué)遙感估算冠層生物量的研究已十分成熟(Curranetal., 1992; 馬澤清等, 2008; 閔志強(qiáng)等, 2010),但光學(xué)遙感只與葉生物量發(fā)生反應(yīng),且在生物量較高時(shí)具有局限性(Sinhaetal., 2015),無(wú)法實(shí)現(xiàn)對(duì)森林冠層生物量的準(zhǔn)確測(cè)定。

    與傳統(tǒng)光學(xué)遙感相比,合成孔徑雷達(dá)(synthetic aperture radar,SAR)波長(zhǎng)較長(zhǎng),對(duì)樹(shù)冠具有一定穿透能力,尤其是高頻率SAR(X波段或C波段)能夠獲得豐富的冠層信息(主要為冠層中枝、葉的散射信息)(Santoroetal., 2007); 同時(shí),SAR還可以穿透云層和一定程度的雨區(qū),具有全天時(shí)全天候觀測(cè)能力(陳爾學(xué), 1999)。目前,SAR估算森林生物量主要利用后向散射和極化干涉技術(shù)(Toanetal., 1992; Dobsonetal., 1992; Metteetal., 2004; 宋茜等, 2011; 李蘭等, 2017),但后向散射存在飽和點(diǎn)問(wèn)題,極化干涉技術(shù)需要多基線或多頻段的全相干極化SAR數(shù)據(jù),對(duì)數(shù)據(jù)質(zhì)量要求較高(Cloudeetal., 1999),而SAR估算冠層生物量的研究多集中于后向散射(Beaudoinetal., 1994; Kasischkeetal., 1995; Harrelletal., 1997)。極化目標(biāo)分解是從極化SAR數(shù)據(jù)中提取地物極化特征的主要方法(Jawaketal., 2015),目前已有許多提取地物信息的極化分解方法(Krogager, 1990; Cloudeetal., 1996; Freemanetal., 1998; Yamaguchietal., 2005),且一些分解方法已用于生物量估算(Goncalvesetal., 2011; Kobayashietal., 2012; Chowdhuryetal., 2013)。

    綜上可以看出,利用SAR數(shù)據(jù)提取極化分解參數(shù)在森林冠層生物量估算中具有一定潛力和優(yōu)勢(shì)。本研究基于內(nèi)蒙古赤峰市旺業(yè)甸林場(chǎng)2017年GF-3全極化SAR數(shù)據(jù),采用Freeman三分量分解、Freeman二分量分解、Yamaguchi三分量分解提取SAR參數(shù),并結(jié)合同一時(shí)期旺業(yè)甸林場(chǎng)外業(yè)調(diào)查數(shù)據(jù),確定極化分解分量與森林冠層生物量的相關(guān)關(guān)系,構(gòu)建冠層-地面散射比參數(shù),應(yīng)用多元逐步回歸方法建立森林冠層生物量與SAR提取參數(shù)回歸模型,旨在探討GF-3全極化SAR數(shù)據(jù)在人工林冠層生物量估算中的潛力,提出一種準(zhǔn)確估算森林冠層生物量的方法。

    1 研究區(qū)概況、數(shù)據(jù)來(lái)源與研究方法

    1.1 研究區(qū)概況 內(nèi)蒙古赤峰市喀喇沁旗旺業(yè)甸林場(chǎng)(41°20′—41°40′N,118°10′—118°29′E)地處燕山山脈北麓七老圖山支脈,平均海拔1 120 m。土壤以典型棕壤為主。屬溫帶季風(fēng)氣候,年均降水量522.6 mm,年均氣溫3.9 ℃,年日照時(shí)間2 700 h以上。林場(chǎng)土地總面積25 958 hm2,有林地面積23 118 hm2,以人工林和通過(guò)封山育林形成的次生林為主,主要樹(shù)種包括油松(Pinustabulaeformis)、華北落葉松(Larixprincipis-rupprechtii)、白樺(Betulaplatyphylla)、黑樺(Betuladahurica)等。

    1.2 數(shù)據(jù)來(lái)源 1) 樣地?cái)?shù)據(jù) 采用2017年10月外業(yè)調(diào)查數(shù)據(jù)。選擇林分整齊、受干擾小、具有代表性的區(qū)域設(shè)置樣地,樣地大小25 m×25 m。使用RTK對(duì)樣地中心點(diǎn)和4個(gè)角點(diǎn)進(jìn)行定位,共22塊樣地分布在SAR數(shù)據(jù)范圍內(nèi),記錄樣地內(nèi)胸徑大于5 cm的樹(shù)木胸徑、樹(shù)高和冠幅,其中1—10號(hào)樣地為油松純林,11—22號(hào)樣地為華北落葉松純林(表1)。

    表1 樣地?cái)?shù)據(jù)統(tǒng)計(jì)Tab.1 Statistics data of plots

    2) 全極化SAR數(shù)據(jù) 采用GF-3合成孔徑雷達(dá)獲取的C波段全極化單視復(fù)數(shù)數(shù)據(jù),數(shù)據(jù)為L(zhǎng)1級(jí)輻射定標(biāo)產(chǎn)品,SAR數(shù)據(jù)獲取于2017年8月5日,記錄方式為16 bit復(fù)數(shù)形式,中心入射角29.7°,中心點(diǎn)坐標(biāo)41°42′N、118°12′E,方位向分辨率5.3 m,距離向分辨率4.5 m,幅寬30 km×30 km(圖1)。

    1.3 研究方法

    1.3.1 樣地生物量計(jì)算 采用本研究區(qū)域油松單木生物量模型(馬欽彥, 1989)和華北落葉松單木生物量模型(羅云建等, 2009)(表2),將樣地每木檢尺數(shù)據(jù)代入模型,分別計(jì)算單木枝、葉生物量。油松枝、葉生物量模型為二元模型,落葉松枝、葉生物量模型為一元模型,將枝、葉生物量之和作為單木冠層生物量,計(jì)算得到樣地內(nèi)所有單木冠層生物量,求和后除以樣地面積(0.062 5 hm2),得到各樣地冠層生物量(表3)。

    圖1 研究區(qū)域及樣地分布Fig.1 Study area and sample plots distribution

    1.3.2 SAR數(shù)據(jù)預(yù)處理 單視復(fù)數(shù)數(shù)據(jù)存在輻射誤差,為了精確反映地物回波特性,對(duì)SAR數(shù)據(jù)進(jìn)行輻射定標(biāo)處理(Shimadaetal., 2009),定標(biāo)系數(shù)根據(jù)數(shù)據(jù)頭文件設(shè)置,獲得SAR的后向散射矩陣[S],并隨之生成極化相干矩陣[T]。

    表2 枝、葉生物量模型①Tab.2 Biomass model of branches and leaves

    表3 樣地冠層生物量Tab.3 Canopy biomass of sample plots

    通常采用多視處理和濾波壓制SAR數(shù)據(jù)固有的斑點(diǎn)噪聲(Leeetal., 2009)。在生成極化相干矩陣[T]的過(guò)程中,方位向和距離向上分別設(shè)置視數(shù)為2,對(duì)SAR數(shù)據(jù)進(jìn)行多視處理,該視數(shù)可保持處理后像元接近方形且盡量與樣地大小相匹配,之后采用Refined Lee方法進(jìn)行濾波,濾波窗口設(shè)置為7(Leeetal., 2002)。

    為了消除地形對(duì)SAR數(shù)據(jù)的影響,應(yīng)用AW3D 30(ALOS World 3D)數(shù)據(jù)對(duì)SAR數(shù)據(jù)進(jìn)行正射校正。AW3D是日本宇宙航空研究開(kāi)發(fā)機(jī)構(gòu)(http:∥www.eorc.jaxa.jp)和日本遙感技術(shù)中心聯(lián)合研制的全新一代高分辨率DEM,數(shù)據(jù)集水平分辨率約30 m,高程精度優(yōu)于5 m,正射校正后SAR數(shù)據(jù)空間分辨率10 m,采用WGS84(world geodetic system-1984)地理坐標(biāo)系和UTM(universal transverse mercator projection)的50 N帶投影坐標(biāo)系。

    1.3.3 目標(biāo)極化分解 極化分解可對(duì)目標(biāo)的散射機(jī)制和特性進(jìn)行合理物理解釋。Cloude等(1996)根據(jù)目標(biāo)散射特性將極化分解方法分為相干分解和非相干分解,本研究區(qū)域?yàn)樯?,地物基本為非相干目?biāo),適合采用非相干分解方法。非相干分解通常以二階統(tǒng)計(jì)的極化協(xié)方差矩陣[C]和極化相干矩陣[T]為基礎(chǔ),一般對(duì)圖像上每個(gè)像素鄰域內(nèi)協(xié)方差或相干矩陣進(jìn)行平均來(lái)估算,得到空間統(tǒng)計(jì)平均后的極化協(xié)方差矩陣〈[C]〉和極化相干矩陣〈[T]〉,且二者可通過(guò)酉矩陣相互轉(zhuǎn)換。極化相干矩陣更易于解釋散射的物理過(guò)程,本研究應(yīng)用極化相干矩陣進(jìn)行極化分解,由于互易性的限制(Cloude, 2010),極化相干矩陣最終簡(jiǎn)化為三維〈[T3]〉。

    采用Freeman三分量分解、Freeman二分量分解、Yamaguchi三分量分解3種極化分解方法,其中后二者對(duì)Freeman三分量分解模型進(jìn)行了不同改進(jìn)。

    1) Freeman三分量分解 Freeman三分量分解(Freemanetal., 1998)是一種基于物理實(shí)際的三分量散射機(jī)制模型,為非相干極化分解方法,對(duì)森林植被后向散射具有良好解譯能力,將森林植被后向散射分解為體散射、二次散射和表面散射3個(gè)在統(tǒng)計(jì)上獨(dú)立不相關(guān)的分量,體散射代表冠層的直接散射,二次散射代表地面與樹(shù)干間的二面角反射,表面散射代表地面的直接單次后向散射。本研究將3個(gè)分解分量功率分別表示為PV、PD、PS:

    〈[T3]〉=fS〈[T3]〉S+fD〈[T3]〉D+fV〈[T3]〉V。

    (1)

    式中:fS、fD、fV為展開(kāi)系數(shù); 〈[T3]〉S、〈[T3]〉D、〈[T3]〉V分別為表面散射、二次散射、體散射的散射模型。

    2) Freeman二分量分解 Freeman(2007)以Freeman三分量分解為基礎(chǔ),提出應(yīng)用于森林的二分量分解模型,該模型的2種散射機(jī)制分別來(lái)自于冠層和地面,體散射分量代表來(lái)自冠層的直接散射,地面散射分量代表地面與樹(shù)干間的二次散射或地面的直接單次后向散射。相比Freeman三分量分解固定的體散射模型,F(xiàn)reeman二分量分解在模型中加入了可變動(dòng)參數(shù),對(duì)體散射模型的限制更小。本研究將體散射、地面散射分量功率分別表示為PC、PG:

    〈[T3]〉=fG〈[T3]〉G+fC〈[T3]〉C。

    (2)

    式中:fG、fC為展開(kāi)系數(shù); 〈[T3]〉G、〈[T3]〉C分別為地面散射、體散射的散射模型。

    3) Yamaguchi三分量分解 Cui等(2012)在Freeman三分量分解的基礎(chǔ)上提出了Yamaguchi三分量分解,該模型與Freeman二分量分解類似,同樣在體散射模型中加入了可變動(dòng)系數(shù),但對(duì)系數(shù)范圍進(jìn)行了限制,且該分解方法首先對(duì)〈[T3]〉進(jìn)行正交變換和酉變換(Anetal., 2010; Singhetal., 2011; Yamaguchietal., 2011),之后應(yīng)用自適應(yīng)體散射模型進(jìn)行極化分解,得到體散射、二次散射、表面散射3個(gè)分量,功率分別表示為PV(θ)、PD(θ)、PS(θ):

    〈[T3(θ)]〉=fS(θ)〈[T3(θ)]〉S+fD(θ)〈[T3(θ)]〉D+

    fV(θ)〈[T3(θ)]〉V。

    (3)

    式中:fS(θ)、fD(θ)、fV(θ)為展開(kāi)系數(shù);〈[T3(θ)]〉S、〈[T3(θ)]〉D、〈[T3(θ)]〉V分別為表面散射、二次散射、體散射的散射模型。

    1.3.4 參數(shù)構(gòu)建 極化分解中,體散射分量代表森林冠層的后向散射信息,二次散射分量、表面散射分量代表SAR信號(hào)透過(guò)冠層后所獲得的樹(shù)干、地表后向散射信息。冠層散射與地面散射比值對(duì)森林冠層結(jié)構(gòu)具有一定敏感性(Freeman, 2007),因此本研究分別將Freeman三分量分解、Yamaguchi三分量分解的體散射分量與二次散射分量和表面散射分量之積進(jìn)行比值運(yùn)算,將Freeman二分量分解體散射分量與地面散射分量進(jìn)行比值運(yùn)算,從而獲得3個(gè)新的參數(shù),各分量均在對(duì)數(shù)變換后參與運(yùn)算,新參數(shù)可以看作各分解方法下的冠-地散射比。參數(shù)構(gòu)建方法如下:

    (4)

    (5)

    (6)

    式中:PV、PD、PS為Freeman三分量分解的體散射、二次散射、表面散射功率;PC、PG為Freeman二分量分解的體散射、地面散射功率;PV(θ)、PD(θ)、PS(θ)為Yamaguchi三分量分解的體散射、二次散射、表面散射功率。

    1.3.5 像元值提取 在提取固定樣地對(duì)應(yīng)的遙感信息時(shí),根據(jù)樣地中心點(diǎn)坐標(biāo)將樣地冠層生物量數(shù)據(jù)與SAR數(shù)據(jù)一一對(duì)應(yīng)。由于樣地面積為0.062 5 hm2,而SAR數(shù)據(jù)像元大小為10 m,所以在提取樣地對(duì)應(yīng)的SAR數(shù)據(jù)信息時(shí),采用八鄰域平均法提取,即樣地中心點(diǎn)對(duì)應(yīng)像元,相鄰8個(gè)像元的像元值求平均,即為該樣地對(duì)應(yīng)的SAR數(shù)據(jù)。

    1.3.6 相關(guān)性分析 相關(guān)性分析是通過(guò)相關(guān)系數(shù)來(lái)表示2個(gè)變量之間相關(guān)關(guān)系密切程度的統(tǒng)計(jì)方法。分別對(duì)冠層生物量、各極化分解分量進(jìn)行以10為底的對(duì)數(shù)變換(Austinetal., 2003; Gamaetal., 2010),將變換后的冠層生物量數(shù)據(jù)與分解分量參數(shù)進(jìn)行相關(guān)性分析,繪制散點(diǎn)圖進(jìn)行直觀評(píng)價(jià),并計(jì)算相關(guān)系數(shù),從而確定變換的冠層生物量與分解分量之間的線性關(guān)系。

    1.3.7 模型構(gòu)建與檢驗(yàn) 應(yīng)用多元逐步回歸方法構(gòu)建模型,其基本思想是將自變量逐個(gè)引入模型中,對(duì)每個(gè)引入的解釋變量進(jìn)行顯著性檢驗(yàn),保留顯著變量,當(dāng)引入的新變量使得原變量變得不再顯著時(shí)剔除原變量,逐步回歸獲得最優(yōu)模型。

    為了評(píng)價(jià)模型的預(yù)測(cè)能力,需對(duì)模型進(jìn)行評(píng)價(jià)檢驗(yàn),通常將數(shù)據(jù)分為建模樣本和檢驗(yàn)樣本。本研究樣地為22塊,若保留部分樣地?cái)?shù)據(jù)用于驗(yàn)證將導(dǎo)致用于建模的樣地?cái)?shù)量過(guò)少,因此運(yùn)用留一法交叉驗(yàn)證(leave-one-outcross-validation,LOOCV)對(duì)模型進(jìn)行評(píng)價(jià)(Geisser, 1974; Stone, 1974; Woldetal., 1984),即每次將1塊樣地作為驗(yàn)證樣本,其余21塊樣地作為訓(xùn)練樣本進(jìn)行建模,不斷重復(fù)該過(guò)程并記錄所有交叉驗(yàn)證結(jié)果,最終獲得22組實(shí)測(cè)值和預(yù)測(cè)值對(duì)原始模型進(jìn)行評(píng)價(jià)。該方法能夠提供模型真實(shí)擬合能力的無(wú)偏估計(jì)(Cawleyetal., 2004),不會(huì)造成數(shù)據(jù)浪費(fèi),且檢驗(yàn)過(guò)程中所得模型與使用全部數(shù)據(jù)所得模型基本一致。采用決定系數(shù)(R2)、均方根誤差(RMSE)、預(yù)測(cè)誤差平方和(PRESS)、平均偏差(ME)、平均絕對(duì)偏差(MAE)、平均相對(duì)偏差(MRE)、平均相對(duì)偏差絕對(duì)值(AMRE)和預(yù)測(cè)精度(P)對(duì)模型擬合結(jié)果和檢驗(yàn)結(jié)果進(jìn)行評(píng)價(jià)(李鳳日, 2004; 董利虎, 2015),計(jì)算公式如下:

    (7)

    (8)

    (9)

    (10)

    (11)

    (12)

    (13)

    其中,

    (14)

    2 結(jié)果與分析

    2.1 極化分解分量與冠層生物量的相關(guān)性 不同極化分解方法所得分量與冠層生物量均有較存在顯著的負(fù)相關(guān)關(guān)系,F(xiàn)reeman三分量分解表面散射分量(PS)與冠層生物量在0.05水平上顯著相關(guān),其余分量均在0.01水平上與冠層生物量顯著相關(guān)(表4),各極化分解方法提取分量對(duì)冠層生物量的變化均較為敏感。Freeman三分量分解體散射分量(PV)與冠層生物量的相關(guān)性高于二次散射分量(PD)和表面散射分量(PS)(圖2a),F(xiàn)reeman二分量分解有著相似的特征(圖2b),但Freeman二分量分解體散射分量(PC)與冠層生物量的相關(guān)性稍好于Freeman三分量分解體散射分量,說(shuō)明在Freeman三分量和二分量分解中代表冠層散射特征的體散射分量對(duì)冠層生物量更為敏感,且Freeman二分量分解比Freeman三分量分解效果更好。Yamaguchi三分量分解體散射分量[PV(θ)]與冠層生物量的相關(guān)性和Freeman三分量分解體散射分量(PV)相似,而二次散射分量[PD(θ)]、表面散射分量[PS(θ)]與冠層生物量的相關(guān)性均高于體散射分量[PV(θ)]。在所有極化分解分量中,Yamaguchi三分量分解二次散射分量[PD(θ)]與冠層生物量的相關(guān)性最高(圖2c)。

    2.2 模型構(gòu)建 本研究將對(duì)數(shù)變換后的冠層生物量(lgW)作為因變量,8個(gè)SAR極化分解分量(表4)和3個(gè)新構(gòu)建冠-地散射比參數(shù)[R1、R2和R3,見(jiàn)式(4)、(5)、(6)]共11個(gè)參數(shù)作為自變量,進(jìn)行逐步回歸,最終冠-地散射比參數(shù)R2、R3進(jìn)入回歸模型,且具有較強(qiáng)的顯著性,共線性容差大于0.5,二者間無(wú)明顯多重共線性(表5)。模型R2為0.658,RMSE為4.943 t·hm-2。

    2.3 模型檢驗(yàn) 將應(yīng)用多元逐步回歸方法獲得的最優(yōu)模型進(jìn)行留一法交叉驗(yàn)證,記錄每次驗(yàn)證預(yù)測(cè)結(jié)果,以交叉檢驗(yàn)所得冠層生物量預(yù)測(cè)值為自變量,以冠層生物量實(shí)測(cè)值為因變量,建立一元線性回歸方程,并繪制散點(diǎn)圖(圖3)和殘差圖(圖4),對(duì)模型進(jìn)行評(píng)價(jià)(徐小軍等, 2011)。結(jié)果表明,模型預(yù)測(cè)誤差均較低,平均偏差(ME)為-0.665 t·hm-2,預(yù)測(cè)模型有輕微高估現(xiàn)象,平均絕對(duì)偏差(MAE)為4.845 t·hm-2,平均相對(duì)偏差(MRE)和平均相對(duì)偏差絕對(duì)值(AMRE)分別為3.33%和23.233%,預(yù)測(cè)精度(P)達(dá)91.5%,說(shuō)明模型存在輕微偏差,但預(yù)測(cè)精度良好。對(duì)冠層生物量實(shí)測(cè)值與預(yù)測(cè)值之間的回歸方程進(jìn)行置信橢圓F檢驗(yàn),其F為0.145,小于F0.05(2,20)=3.49,表明實(shí)測(cè)值與預(yù)測(cè)值之間差異不明顯(胡珍珠等, 2016; 王樹(shù)文等, 2016),預(yù)測(cè)值大致分布在1∶1直線接近,模型未出現(xiàn)飽和點(diǎn)。模型標(biāo)準(zhǔn)化殘差全部在2倍標(biāo)準(zhǔn)殘差內(nèi),且呈隨機(jī)分布,說(shuō)明模型估算效果受生物量變化影響較小,具有良好的適用性和預(yù)測(cè)性。

    表4 極化分解分量與冠層生物量的相關(guān)系數(shù)①Tab.4 Correlation coefficient of polarization decomposition component and canopy biomass

    表5 模型回歸系數(shù)Tab.5 Regression coefficients of the model

    圖3 實(shí)測(cè)值與模型預(yù)測(cè)值的相關(guān)性Fig.3 The correlation between measured value and model predicted value

    圖4 實(shí)測(cè)值-模型預(yù)測(cè)值標(biāo)準(zhǔn)殘差Fig.4 Standardized residuals of measured value-model predicted value

    3 討論

    3.1 極化分解分量與冠層生物量的關(guān)系 Kobayashi等(2012)采用Yamaguchi四分量分解估算熱帶森林蓄積,極化相干矩陣經(jīng)正交變換后,表面散射、體散射與蓄積量的相關(guān)性變得稍好,表面散射與蓄積量為負(fù)相關(guān)關(guān)系,體散射與蓄積量為正相關(guān)關(guān)系,二次散射與蓄積量無(wú)相關(guān)性。Chowdhury等(2013)同樣采用Yamaguchi四分量分解估算寒帶森林蓄積,極化相干矩陣經(jīng)正交變換后,二次散射與蓄積量的相關(guān)性顯著提高,但對(duì)表面散射和體散射幾乎沒(méi)有影響。在Goncalves等(2011)的研究中,采用Freeman三分量分解估算熱帶森林蓄積,蓄積量與3個(gè)分解分量之間均呈顯著正相關(guān)關(guān)系。

    本研究采用3種極化分解方法提取分解分量,極化分解分量與冠層生物量均存在較為顯著的相關(guān)關(guān)系,與Freeman三分量分解相比,F(xiàn)reeman二分量分解、Yamaguchi三分量分解的極化分解分量與冠層生物量相關(guān)性更好,二者對(duì)體散射模型進(jìn)行了改進(jìn),均在體散射模型中加入可變系數(shù),說(shuō)明可變系數(shù)加入可提高森林區(qū)域極化分解效果,能夠更準(zhǔn)確地獲取森林垂直結(jié)構(gòu)信息。同時(shí),Yamaguchi三分量分解效果稍好于Freeman二分量分解,Yamaguchi三分量分解在極化分解前對(duì)極化相干矩陣進(jìn)行正交變換,改善了極化分解效果,從而提高了與冠層生物量的相關(guān)性,與Kobayashi等(2012)和Chowdhury等(2013)的研究結(jié)果在理論上一致。不同的是,本研究所得極化分解分量與冠層生物量均為負(fù)相關(guān)關(guān)系,這可能是因?yàn)槎紊⑸浜捅砻嫔⑸錇镾AR信號(hào)透過(guò)冠層的后向散射功率,冠層生物量越大,枝葉分布越密集,二次散射和表面散射功率越低; 體散射為SAR信號(hào)在冠層中的散射過(guò)程,C波段SAR穿透性較弱,隨著冠層生物量增加,枝葉對(duì)SAR信號(hào)衰減作用增強(qiáng),使得體散射功率降低。此外,分解分量特征受含水率、介電常數(shù)、灌木層、地形等多種因素影響, 森林結(jié)構(gòu)、樹(shù)葉及樹(shù)干的形狀和大小對(duì)生物量估算也存在一定影響(Santoroetal., 2009; Hoekmanetal., 2002; Rowlandetal., 2008)。

    3.2 冠層生物量估算效果 Goncalves等(2011)應(yīng)用多元線性回歸建立了包括Freeman三分量分解在內(nèi)的SAR參數(shù)與森林蓄積間的關(guān)系; Chowdhury等(2013)構(gòu)建體散射分量、二次散射分量之積與表面散射分量的比值,與單獨(dú)使用分解分量相比顯著提高了與森林蓄積間的相關(guān)性。類似地,本研究構(gòu)建各極化分解方法體散射分量與二次散射分量、表面散射分量之積的比值,應(yīng)用多元線性回歸建立SAR提取參數(shù)與森林冠層生物量回歸模型,得到了較好結(jié)果(R2=0.658,RMSE=4.943 t·hm-2),相比于單一分解分量,冠-地散射比參數(shù)定量化描述了森林冠層散射與冠層下部散射間的關(guān)系,對(duì)森林的極化分解信息有效結(jié)合,當(dāng)比值較大時(shí),冠層生物量較高; 當(dāng)比值較小時(shí),冠層生物量較低。

    由于數(shù)據(jù)限制,本研究?jī)H采用22塊樣地?cái)?shù)據(jù)進(jìn)行建模,在統(tǒng)計(jì)學(xué)中屬于小樣本,雖然具有一定說(shuō)服力,但還需要更多的樣地觀測(cè)數(shù)據(jù)進(jìn)行驗(yàn)證(Garestieretal., 2009; Neumannetal., 2012; Gamaetal., 2010; Goncalvesetal., 2011)。同時(shí),本研究?jī)H對(duì)3種極化分解方法進(jìn)行分析,更多分解方法對(duì)冠層生物量的估算效果有待進(jìn)一步研究。

    4 結(jié)論

    1) 3種極化分解方法獲得的極化分解分量均與冠層生物量具有顯著相關(guān)關(guān)系,采用極化分解手段估算冠層生物量具有可行性。

    2) Freeman二分量分解、Yamaguchi三分量分解的極化分解分量與冠層生物量相關(guān)性更高,說(shuō)明極化相干矩陣旋轉(zhuǎn)變換、體散射模型優(yōu)化可有效提高極化分解效果,能夠獲取更豐富的森林垂直結(jié)構(gòu)信息。

    3) 相比于各極化分解分量,冠-地散射比參數(shù)對(duì)冠層生物量的敏感性更高,多種SAR極化分解參數(shù)共同使用能夠較好估算冠層生物量,且未出現(xiàn)明顯的飽和點(diǎn)。

    應(yīng)該指出的是,由于C波段SAR數(shù)據(jù)穿透性較弱,且本研究?jī)H對(duì)Freeman三分量分解、Freeman二分量分解和Yamaguchi三分量分解3種極化分解方法估算冠層生物量的效果進(jìn)行了探討,應(yīng)進(jìn)一步研究其他分解方法估算冠層生物量是否與本研究獲得相同的結(jié)論,驗(yàn)證本研究方法對(duì)估算地上生物量是否具有適用性。

    猜你喜歡
    冠層樣地極化
    認(rèn)知能力、技術(shù)進(jìn)步與就業(yè)極化
    基于低空遙感的果樹(shù)冠層信息提取方法研究
    基于激光雷達(dá)的樹(shù)形靶標(biāo)冠層葉面積探測(cè)模型研究
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    安徽省淮南森林冠層輻射傳輸過(guò)程的特征
    昆明市主要綠化樹(shù)種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    雙頻帶隔板極化器
    施氮水平對(duì)冬小麥冠層氨揮發(fā)的影響
    基于PWM控制的新型極化電源設(shè)計(jì)與實(shí)現(xiàn)
    精品国内亚洲2022精品成人| 美女免费视频网站| 亚洲av成人av| 久久九九热精品免费| 亚洲男人的天堂狠狠| 国产激情欧美一区二区| 麻豆久久精品国产亚洲av| 高清日韩中文字幕在线| 丝袜美腿在线中文| 亚洲人成网站在线播| 久久精品亚洲精品国产色婷小说| 成人三级黄色视频| 欧美激情在线99| 首页视频小说图片口味搜索| 久久精品亚洲精品国产色婷小说| 午夜福利欧美成人| 日日干狠狠操夜夜爽| 夜夜夜夜夜久久久久| 国产激情偷乱视频一区二区| 91字幕亚洲| 国产国拍精品亚洲av在线观看 | 窝窝影院91人妻| 精品一区二区三区视频在线观看免费| 深夜精品福利| 麻豆国产av国片精品| 制服人妻中文乱码| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 老司机福利观看| 日本熟妇午夜| 亚洲成人精品中文字幕电影| 天美传媒精品一区二区| 国产精品久久久久久久久免 | 成年版毛片免费区| 久久婷婷人人爽人人干人人爱| 两个人视频免费观看高清| 性欧美人与动物交配| 中文资源天堂在线| av国产免费在线观看| 叶爱在线成人免费视频播放| 人妻久久中文字幕网| 国产探花极品一区二区| 色噜噜av男人的天堂激情| 国产高清有码在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 午夜两性在线视频| 亚洲精品成人久久久久久| 国内毛片毛片毛片毛片毛片| 成人亚洲精品av一区二区| 伊人久久精品亚洲午夜| 久久国产乱子伦精品免费另类| 亚洲狠狠婷婷综合久久图片| 性欧美人与动物交配| 国产精品99久久久久久久久| 日本与韩国留学比较| 麻豆国产97在线/欧美| 最近在线观看免费完整版| 国产成人av激情在线播放| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 免费一级毛片在线播放高清视频| 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 欧美bdsm另类| 一本久久中文字幕| a级一级毛片免费在线观看| 国产免费一级a男人的天堂| 老汉色av国产亚洲站长工具| 免费在线观看成人毛片| 久久精品亚洲精品国产色婷小说| 欧美av亚洲av综合av国产av| 亚洲精华国产精华精| 岛国视频午夜一区免费看| 嫁个100分男人电影在线观看| 国产 一区 欧美 日韩| 亚洲精品在线美女| 国产精品一区二区三区四区久久| 夜夜夜夜夜久久久久| 特级一级黄色大片| 免费av毛片视频| 俺也久久电影网| 亚洲人成网站高清观看| 国产爱豆传媒在线观看| 99国产综合亚洲精品| 国产97色在线日韩免费| 一个人看的www免费观看视频| 99久久精品热视频| 日韩亚洲欧美综合| 在线国产一区二区在线| 男女下面进入的视频免费午夜| 国产精品 国内视频| 国产综合懂色| 久久欧美精品欧美久久欧美| 偷拍熟女少妇极品色| 欧美中文日本在线观看视频| 麻豆成人av在线观看| 午夜免费观看网址| 99精品久久久久人妻精品| 亚洲av中文字字幕乱码综合| 大型黄色视频在线免费观看| 欧美激情在线99| 国内久久婷婷六月综合欲色啪| 亚洲人成网站在线播放欧美日韩| 麻豆国产97在线/欧美| 免费在线观看亚洲国产| 欧美日韩精品网址| 哪里可以看免费的av片| 国产97色在线日韩免费| 一个人看的www免费观看视频| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 亚洲午夜理论影院| 日本 av在线| 十八禁人妻一区二区| 高潮久久久久久久久久久不卡| 国产精品综合久久久久久久免费| 天美传媒精品一区二区| 国产精品一区二区免费欧美| 亚洲专区中文字幕在线| 欧美zozozo另类| 日韩欧美国产在线观看| 午夜福利高清视频| 久久6这里有精品| www.熟女人妻精品国产| 免费在线观看影片大全网站| 久久久精品欧美日韩精品| 99久久综合精品五月天人人| 亚洲一区二区三区色噜噜| 久久精品综合一区二区三区| 色视频www国产| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 一本一本综合久久| 亚洲成av人片在线播放无| 啦啦啦免费观看视频1| 欧美成人性av电影在线观看| 怎么达到女性高潮| 在线视频色国产色| 午夜精品在线福利| 国产精品99久久久久久久久| 国产私拍福利视频在线观看| 首页视频小说图片口味搜索| 精品国产超薄肉色丝袜足j| 变态另类丝袜制服| 成熟少妇高潮喷水视频| 99久久久亚洲精品蜜臀av| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 亚洲av二区三区四区| 日本精品一区二区三区蜜桃| 久久草成人影院| 99视频精品全部免费 在线| 国产精品三级大全| 首页视频小说图片口味搜索| 国内精品美女久久久久久| 乱人视频在线观看| 黄色丝袜av网址大全| 日韩欧美 国产精品| 日本免费a在线| 国产精品 欧美亚洲| 身体一侧抽搐| 久久久久久久精品吃奶| 国产伦精品一区二区三区四那| 成人三级黄色视频| 国产真实乱freesex| 欧美乱码精品一区二区三区| www.色视频.com| 成人av在线播放网站| 欧美+日韩+精品| 精品国产美女av久久久久小说| 国产精华一区二区三区| 日韩国内少妇激情av| 午夜免费观看网址| 成年女人毛片免费观看观看9| 女同久久另类99精品国产91| 欧美一区二区亚洲| 久久人妻av系列| 亚洲精品亚洲一区二区| 桃红色精品国产亚洲av| 欧美+日韩+精品| 国产91精品成人一区二区三区| 一进一出抽搐动态| 国产精品亚洲一级av第二区| 90打野战视频偷拍视频| www.熟女人妻精品国产| 在线视频色国产色| 亚洲av免费高清在线观看| 在线播放无遮挡| 国产主播在线观看一区二区| 一进一出抽搐动态| 亚洲中文字幕日韩| 久久6这里有精品| 一级毛片女人18水好多| 婷婷六月久久综合丁香| 久99久视频精品免费| 久久久久国产精品人妻aⅴ院| 深夜精品福利| 日本三级黄在线观看| 欧美一区二区精品小视频在线| 一夜夜www| 亚洲人成网站高清观看| 亚洲美女视频黄频| 国产成人欧美在线观看| 国产乱人视频| 亚洲五月婷婷丁香| 手机成人av网站| 亚洲精品国产精品久久久不卡| 熟女少妇亚洲综合色aaa.| 国产高潮美女av| 69人妻影院| 亚洲男人的天堂狠狠| 欧美区成人在线视频| 国产精品久久久人人做人人爽| 听说在线观看完整版免费高清| 精品一区二区三区视频在线 | 午夜福利在线观看免费完整高清在 | 色av中文字幕| 色哟哟哟哟哟哟| 一夜夜www| 脱女人内裤的视频| 天堂影院成人在线观看| 中文字幕人成人乱码亚洲影| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 欧美性猛交╳xxx乱大交人| 中文在线观看免费www的网站| 免费看a级黄色片| 午夜福利视频1000在线观看| 免费av观看视频| 色综合亚洲欧美另类图片| 中文字幕精品亚洲无线码一区| 麻豆成人av在线观看| 欧美性感艳星| 亚洲av中文字字幕乱码综合| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 中文亚洲av片在线观看爽| 一夜夜www| 亚洲在线观看片| 久久久久久九九精品二区国产| 亚洲色图av天堂| 国产精品久久久久久久久免 | 亚洲精品国产精品久久久不卡| 人妻久久中文字幕网| 又粗又爽又猛毛片免费看| 日韩高清综合在线| 精品福利观看| 叶爱在线成人免费视频播放| 午夜免费男女啪啪视频观看 | 在线观看舔阴道视频| 99热只有精品国产| 丰满人妻熟妇乱又伦精品不卡| av天堂在线播放| 久久久久国内视频| 怎么达到女性高潮| 欧美绝顶高潮抽搐喷水| 欧美日韩福利视频一区二区| 99久国产av精品| 婷婷亚洲欧美| 97人妻精品一区二区三区麻豆| 国产精品免费一区二区三区在线| 俺也久久电影网| 啦啦啦观看免费观看视频高清| 日韩免费av在线播放| 中文字幕熟女人妻在线| 国产日本99.免费观看| 国产成+人综合+亚洲专区| 一级a爱片免费观看的视频| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 精品国产亚洲在线| 国产精品三级大全| xxxwww97欧美| 久久精品国产亚洲av涩爱 | 午夜日韩欧美国产| 少妇的逼好多水| 中文在线观看免费www的网站| 亚洲成人中文字幕在线播放| 国产单亲对白刺激| 国产午夜精品久久久久久一区二区三区 | 日韩欧美在线二视频| 国产午夜精品久久久久久一区二区三区 | 啦啦啦韩国在线观看视频| 99久久精品热视频| 亚洲不卡免费看| 午夜福利18| 亚洲国产色片| 日韩中文字幕欧美一区二区| 国产成人影院久久av| 99久久九九国产精品国产免费| 18禁在线播放成人免费| 国产亚洲精品综合一区在线观看| 哪里可以看免费的av片| 99热精品在线国产| 亚洲精品美女久久久久99蜜臀| 国产高清videossex| eeuss影院久久| 在线看三级毛片| 国产一区二区三区视频了| 国模一区二区三区四区视频| 国产精品99久久99久久久不卡| 欧美日韩乱码在线| 国产亚洲精品一区二区www| 久久精品综合一区二区三区| 精品久久久久久久人妻蜜臀av| 在线免费观看的www视频| 日本在线视频免费播放| 久久精品夜夜夜夜夜久久蜜豆| 美女高潮的动态| 日本免费一区二区三区高清不卡| 午夜精品久久久久久毛片777| 亚洲电影在线观看av| 18美女黄网站色大片免费观看| 伊人久久精品亚洲午夜| 亚洲七黄色美女视频| netflix在线观看网站| 久久香蕉国产精品| 成人欧美大片| 欧美日韩中文字幕国产精品一区二区三区| 日本一本二区三区精品| 亚洲av美国av| 亚洲av成人精品一区久久| 日韩欧美一区二区三区在线观看| 成人18禁在线播放| 国产蜜桃级精品一区二区三区| 婷婷丁香在线五月| 97人妻精品一区二区三区麻豆| 成年免费大片在线观看| 国产真实伦视频高清在线观看 | 亚洲国产色片| 午夜福利高清视频| 国产麻豆成人av免费视频| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 人妻夜夜爽99麻豆av| 亚洲五月婷婷丁香| 色哟哟哟哟哟哟| 亚洲成av人片免费观看| 国产99白浆流出| 美女cb高潮喷水在线观看| 小说图片视频综合网站| 桃红色精品国产亚洲av| 不卡一级毛片| 别揉我奶头~嗯~啊~动态视频| 久久精品91蜜桃| 久久精品91无色码中文字幕| 麻豆成人av在线观看| 精品久久久久久久末码| 亚洲成人免费电影在线观看| 男人和女人高潮做爰伦理| 国产视频一区二区在线看| 欧美成人一区二区免费高清观看| 一级黄色大片毛片| 久久性视频一级片| 我的老师免费观看完整版| 男人和女人高潮做爰伦理| 精品国产超薄肉色丝袜足j| 美女cb高潮喷水在线观看| 日韩中文字幕欧美一区二区| 久久九九热精品免费| 成年女人看的毛片在线观看| tocl精华| 国内少妇人妻偷人精品xxx网站| 在线观看日韩欧美| 97碰自拍视频| 最新在线观看一区二区三区| 在线视频色国产色| 精品99又大又爽又粗少妇毛片 | 日韩成人在线观看一区二区三区| 欧美成人免费av一区二区三区| 午夜福利18| 欧美中文日本在线观看视频| 少妇丰满av| 熟女少妇亚洲综合色aaa.| 午夜a级毛片| 国产精品久久电影中文字幕| 黄色丝袜av网址大全| 草草在线视频免费看| 欧美另类亚洲清纯唯美| 99久久精品热视频| 亚洲人成网站高清观看| 国产一区二区在线av高清观看| 啦啦啦韩国在线观看视频| 女生性感内裤真人,穿戴方法视频| 99久久精品一区二区三区| 久久久久久久午夜电影| 亚洲中文日韩欧美视频| 日本免费a在线| 757午夜福利合集在线观看| av视频在线观看入口| 亚洲色图av天堂| 久久精品91无色码中文字幕| 日本黄大片高清| 午夜精品久久久久久毛片777| 国产亚洲精品久久久久久毛片| 制服丝袜大香蕉在线| 啦啦啦观看免费观看视频高清| 欧美极品一区二区三区四区| 色综合婷婷激情| 精华霜和精华液先用哪个| 日韩欧美精品免费久久 | 亚洲精品在线观看二区| 免费看日本二区| 午夜福利欧美成人| 中文在线观看免费www的网站| 美女高潮喷水抽搐中文字幕| 日日摸夜夜添夜夜添小说| 90打野战视频偷拍视频| 伊人久久精品亚洲午夜| 成人无遮挡网站| 日韩中文字幕欧美一区二区| 日韩有码中文字幕| 色在线成人网| svipshipincom国产片| 在线免费观看的www视频| 啦啦啦韩国在线观看视频| 淫秽高清视频在线观看| 天天一区二区日本电影三级| 日韩欧美精品免费久久 | 女生性感内裤真人,穿戴方法视频| 中文字幕高清在线视频| а√天堂www在线а√下载| 国产精品99久久99久久久不卡| 一本久久中文字幕| 亚洲 国产 在线| 欧美丝袜亚洲另类 | 久久久久久人人人人人| 日本撒尿小便嘘嘘汇集6| 亚洲成人中文字幕在线播放| 一本精品99久久精品77| 成年女人永久免费观看视频| 香蕉久久夜色| 中亚洲国语对白在线视频| 精品一区二区三区视频在线观看免费| 午夜激情欧美在线| 啦啦啦免费观看视频1| 99久久成人亚洲精品观看| 亚洲精品色激情综合| 国产伦精品一区二区三区四那| 真实男女啪啪啪动态图| 亚洲av美国av| 18+在线观看网站| 18禁美女被吸乳视频| 欧美日韩精品网址| 国产精品嫩草影院av在线观看 | 中文资源天堂在线| 母亲3免费完整高清在线观看| 成人特级黄色片久久久久久久| 精品熟女少妇八av免费久了| 亚洲国产精品999在线| 亚洲国产精品久久男人天堂| 国产国拍精品亚洲av在线观看 | 亚洲精品国产精品久久久不卡| 搡老熟女国产l中国老女人| 免费人成视频x8x8入口观看| 精品一区二区三区人妻视频| 久久久久性生活片| 国产精品av视频在线免费观看| 国产综合懂色| 一个人免费在线观看的高清视频| 757午夜福利合集在线观看| 国产成人影院久久av| 日韩免费av在线播放| 国产单亲对白刺激| 国产91精品成人一区二区三区| 久久人妻av系列| 黄片小视频在线播放| 最新在线观看一区二区三区| 级片在线观看| 村上凉子中文字幕在线| 国产成人av教育| 精品国产亚洲在线| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 老汉色av国产亚洲站长工具| 国内久久婷婷六月综合欲色啪| 午夜福利欧美成人| 亚洲久久久久久中文字幕| 精品乱码久久久久久99久播| 亚洲国产精品久久男人天堂| 亚洲专区中文字幕在线| 国产成人系列免费观看| 熟女人妻精品中文字幕| 久久香蕉精品热| 久久天躁狠狠躁夜夜2o2o| 欧美日韩黄片免| 夜夜夜夜夜久久久久| www.色视频.com| a级一级毛片免费在线观看| 男女之事视频高清在线观看| 日韩中文字幕欧美一区二区| aaaaa片日本免费| 国产伦精品一区二区三区视频9 | 国产av在哪里看| 亚洲国产中文字幕在线视频| 一二三四社区在线视频社区8| 国产精品亚洲美女久久久| 欧美xxxx黑人xx丫x性爽| 精品电影一区二区在线| 精品99又大又爽又粗少妇毛片 | 宅男免费午夜| 韩国av一区二区三区四区| 久久精品91无色码中文字幕| 欧美+亚洲+日韩+国产| 人妻久久中文字幕网| 热99在线观看视频| 久久伊人香网站| 国产精品久久久人人做人人爽| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 欧美精品啪啪一区二区三区| 国产淫片久久久久久久久 | 国产精品免费一区二区三区在线| 中文资源天堂在线| a在线观看视频网站| 国产av在哪里看| 国产在视频线在精品| 露出奶头的视频| 黄色日韩在线| 午夜激情福利司机影院| 国产精品久久久人人做人人爽| 亚洲中文日韩欧美视频| 18美女黄网站色大片免费观看| 男人舔奶头视频| 成熟少妇高潮喷水视频| 亚洲精品在线观看二区| 久久精品国产自在天天线| 最新在线观看一区二区三区| 国产男靠女视频免费网站| 别揉我奶头~嗯~啊~动态视频| 身体一侧抽搐| 精品一区二区三区视频在线 | 欧美日韩乱码在线| 亚洲最大成人中文| 99国产精品一区二区蜜桃av| 婷婷精品国产亚洲av| 中国美女看黄片| 日本黄色视频三级网站网址| 久久午夜亚洲精品久久| ponron亚洲| 可以在线观看的亚洲视频| 波多野结衣高清无吗| 97超级碰碰碰精品色视频在线观看| 色在线成人网| 欧美极品一区二区三区四区| 久久99热这里只有精品18| 日韩欧美国产在线观看| 天堂动漫精品| 国产激情偷乱视频一区二区| 天天添夜夜摸| 国产精品乱码一区二三区的特点| 免费人成视频x8x8入口观看| 久久国产精品影院| 亚洲最大成人手机在线| 亚洲专区中文字幕在线| 色视频www国产| 婷婷六月久久综合丁香| 亚洲成av人片在线播放无| a在线观看视频网站| 嫩草影院精品99| 99视频精品全部免费 在线| 日本黄色片子视频| 亚洲熟妇中文字幕五十中出| 成年版毛片免费区| 国产中年淑女户外野战色| 免费电影在线观看免费观看| 国产精品亚洲美女久久久| 国产精品嫩草影院av在线观看 | 欧美最新免费一区二区三区 | 校园春色视频在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲专区国产一区二区| 在线观看日韩欧美| 人人妻人人澡欧美一区二区| 可以在线观看的亚洲视频| 最近最新中文字幕大全电影3| 美女黄网站色视频| 国产成人a区在线观看| 久久精品综合一区二区三区| 亚洲国产日韩欧美精品在线观看 | 欧美丝袜亚洲另类 | 日本 av在线| 亚洲一区高清亚洲精品| 亚洲国产精品sss在线观看| 国产精品 欧美亚洲| 中文字幕人妻丝袜一区二区| 久久久久久久久中文| 黄片大片在线免费观看| av国产免费在线观看| 一区二区三区国产精品乱码| 日本一二三区视频观看| 男人舔女人下体高潮全视频| 亚洲国产欧美网| 久久性视频一级片| 最新中文字幕久久久久| 99久久99久久久精品蜜桃| 午夜精品在线福利| 婷婷六月久久综合丁香| 欧美乱色亚洲激情| 亚洲成人久久爱视频| 亚洲,欧美精品.| 白带黄色成豆腐渣| 国产亚洲欧美98| 有码 亚洲区| 国产亚洲精品久久久久久毛片| 精品午夜福利视频在线观看一区| 亚洲国产欧美人成| 欧美区成人在线视频| 人妻久久中文字幕网|