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

    干雪深度反演的同極化相位差模型

    2021-08-14 01:36:34宋依娜肖鵬峰張學(xué)良
    測(cè)繪學(xué)報(bào) 2021年7期
    關(guān)鍵詞:冰粒雪深積雪

    宋依娜,肖鵬峰,2,3,張學(xué)良,卓 越,馬 威

    1.南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023;2.自然資源部國(guó)土衛(wèi)星遙感應(yīng)用重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023;3.江蘇省地理信息技術(shù)重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023

    積雪是冰凍圈和水圈重要的組成部分,在全球變化研究中具有重要地位[1]。一方面,積雪作為一種特殊的地表覆蓋物,通過(guò)反射大部分太陽(yáng)輻射以及在消融過(guò)程中吸收大量能量,對(duì)全球氣候系統(tǒng)產(chǎn)生了顯著的降溫作用[2-3]。同時(shí),積雪融水作為全球水資源的重要組成部分,為全球1/6人口提供了水源[4],并進(jìn)一步產(chǎn)生了水力發(fā)電、農(nóng)業(yè)灌溉等其他用途[5]。另一方面,積雪具有時(shí)空不穩(wěn)定性,過(guò)量或非適時(shí)的積雪也會(huì)引發(fā)自然災(zāi)害,給人類(lèi)生命和財(cái)產(chǎn)安全帶來(lái)巨大的威脅和損失[6]。而積雪深度作為積雪的重要因子,表征局部氣候環(huán)境特征與水資源條件,是水文預(yù)測(cè)、氣候模擬、雪災(zāi)監(jiān)測(cè)與評(píng)價(jià)等模型的重要輸入?yún)?shù)[7]。獲取高精度雪深空間分布信息可為水資源管理、氣候變化研究和防災(zāi)減災(zāi)工程等提供科學(xué)支撐。

    合成孔徑雷達(dá)(synthetic aperture radar,SAR)因其有較高的空間分辨率[8]且對(duì)積雪特性有著較高的敏感度等優(yōu)勢(shì)[9-10],在流域尺度的雪深反演中具有較強(qiáng)的適用性[11-12]。極化SAR通過(guò)回波獲取關(guān)于積雪的散射特征[13],被廣泛應(yīng)用于積雪相關(guān)研究,如干濕雪的分類(lèi)、雪濕度和雪密度等參數(shù)的反演[14-16],但是目前的研究對(duì)極化SAR的極化相位信息的分析較少[17]。極化SAR計(jì)算同一天線、不同極化方式的極化相干,可獲取極化相干系數(shù)和極化相位信息,包括同極化相位差(co-polarized phase difference,CPD)。積雪為各向異性介質(zhì),導(dǎo)致HH和VV極化信號(hào)在雪層中具有不同的傳播速度,因此產(chǎn)生CPD。文獻(xiàn)[18]根據(jù)此原理建立積雪CPD模型,定量分析了CPD與積雪各向異性結(jié)構(gòu)、雪深的關(guān)系。但是,根據(jù)CPD模型反演干雪深度需要冰粒的各向異性結(jié)構(gòu)(以冰粒軸間比表示)和積雪密度兩個(gè)參數(shù),其中冰粒軸間比參數(shù)比較抽象,難以通過(guò)野外觀測(cè)獲得。而且,目前國(guó)際上基于CPD模型進(jìn)行雪深反演的研究主要使用X波段,暫未見(jiàn)使用C波段極化SAR數(shù)據(jù)的研究報(bào)道。文獻(xiàn)[19]假設(shè)研究區(qū)積雪各向異性介電常數(shù)不變,使用一個(gè)實(shí)測(cè)站點(diǎn)處的積雪密度0.07 g/cm3和假定軸間比1.5作為雪深反演算法的輸入,利用全極化TerraSAR-X數(shù)據(jù)反演整個(gè)影像的雪深,并以該實(shí)測(cè)站點(diǎn)的雪深進(jìn)行驗(yàn)證,獲得較高的反演精度。但是,由于積雪的各向異性、密度空間分布不均,僅使用一個(gè)實(shí)測(cè)點(diǎn)的積雪密度和假定軸間比反演整個(gè)影像的雪深,將給反演結(jié)果帶來(lái)較大誤差。文獻(xiàn)[20]根據(jù)實(shí)測(cè)點(diǎn)處CPD的正負(fù),結(jié)合該點(diǎn)處實(shí)測(cè)雪深、雪密度等參數(shù),代入CPD正演模型,求得冰粒軸間比,再借助積雪密度反演算法[21]獲取全局雪深分布,但此反演方法較為復(fù)雜,且積雪密度反演算法可能會(huì)帶來(lái)新的誤差。

    本文基于CPD模型提出了新的雪深反演方法,將反演所需要的參數(shù)減少為遙感獲取的CPD數(shù)據(jù)以及實(shí)測(cè)雪深數(shù)據(jù),根據(jù)全部實(shí)測(cè)點(diǎn)的雪深與CPD的關(guān)系,擬合得到最適合研究區(qū)積雪特性的反演模型系數(shù),用它們反演研究區(qū)的雪深。為了驗(yàn)證反演方法在C波段、不同積雪條件下的應(yīng)用能力,使用GF-3全極化數(shù)據(jù)計(jì)算CPD,以新疆阿爾泰山南坡克蘭河上游為研究區(qū),根據(jù)實(shí)測(cè)雪情差異劃分深雪區(qū)與淺雪區(qū)進(jìn)行雪深反演,分析反演結(jié)果的差異,以期為山區(qū)雪深反演提供新思路。

    1 研究區(qū)與數(shù)據(jù)

    1.1 研究區(qū)概況

    克蘭河流域位于新疆北部,阿爾泰山南麓,在北溫帶大陸性氣候條件影響下,所形成的積雪為典型的“干寒型”積雪[22],年平均雪深在40 cm以上??颂m河屬融雪徑流補(bǔ)給為主的河流,年平均徑流量6.0×108m3[23]。積雪融水對(duì)當(dāng)?shù)剞r(nóng)業(yè)灌溉及畜牧業(yè)發(fā)展具有重要影響,研究該流域的積雪深度具有重要意義。

    研究區(qū)范圍如圖1所示,位于克蘭河流域上游東支小東溝源頭區(qū)域。研究區(qū)西南部為戈壁,地形起伏較為平緩;東北部是阿爾泰山南麓侵蝕中山山地,海拔較高,地形起伏較大,是典型的山區(qū)環(huán)境[24]。

    1.2 研究數(shù)據(jù)

    (1)SAR數(shù)據(jù)。CPD為HH和VV極化通道進(jìn)行極化相干得到復(fù)相干系數(shù)的相位角,SAR衛(wèi)星必須具有HH和VV極化通道極化相干成像的能力才能計(jì)算CPD,目前常見(jiàn)的可以計(jì)算CPD的SAR衛(wèi)星有Radarsat-2、TerraSAR-X/TanDEM-X、ALOS-PALSAR等。研究表明C波段是CPD對(duì)雪深變化仍有較高敏感性的頻率下限[25],本文研究選取我國(guó)的GF-3衛(wèi)星數(shù)據(jù),探討基于CPD模型使用C波段極化SAR數(shù)據(jù)反演雪深的能力。GF-3衛(wèi)星于2016年8月10日發(fā)射,最高分辨率為1 m,設(shè)計(jì)有12種成像模式[26]。綜合考慮GF-3衛(wèi)星數(shù)據(jù)的成像時(shí)間和圖像覆蓋區(qū)域、成像質(zhì)量、同步觀測(cè)試驗(yàn)可行性等,選取2018年1月17日獲取的阿勒泰地區(qū)GF-3衛(wèi)星全極化條帶數(shù)據(jù)作為實(shí)驗(yàn)數(shù)據(jù),影像范圍為47°49′—48°07′N(xiāo),88°03′—88°24′E。研究所用的SAR影像相關(guān)參數(shù)見(jiàn)表1。

    (2)野外實(shí)測(cè)數(shù)據(jù)。CPD反演雪深的可行性分析、模型參數(shù)的確定與反演結(jié)果的精度驗(yàn)證均需地面實(shí)測(cè)數(shù)據(jù)的支持,因此于2018年1月17日對(duì)研究區(qū)進(jìn)行星地同步積雪觀測(cè),設(shè)計(jì)了積雪樣帶觀測(cè)和積雪屬性剖面觀測(cè)相結(jié)合的觀測(cè)方案。此時(shí)為研究區(qū)的積雪穩(wěn)定期,觀測(cè)期間最高氣溫在0℃以下,未發(fā)生降水。積雪樣帶觀測(cè)共計(jì)44個(gè)樣點(diǎn),主要測(cè)量雪深、積雪密度、雪層溫度等參數(shù)。積雪剖面觀測(cè)共計(jì)8個(gè)剖面點(diǎn),主要測(cè)量雪深、雪層密度、雪層溫度、雪層介電常數(shù)、液態(tài)水含量等參數(shù)。雪深直接由鋼制直尺測(cè)量,取觀測(cè)點(diǎn)周?chē)? m的范圍內(nèi)3次測(cè)量的平均值,測(cè)量精度為0.1 cm;積雪密度由鋼制直尺、托板和自制的量雪筒測(cè)量,測(cè)量精度為1 g/cm3;自積雪表層向下每10 cm劃分雪層,測(cè)得積雪參數(shù)代表該雪層積雪參數(shù),其中雪層溫度利用經(jīng)過(guò)標(biāo)定的針式溫度計(jì)進(jìn)行測(cè)量,測(cè)量精度為0.1℃,雪層介電常數(shù)、密度和液態(tài)水含量通過(guò)雪特性分析儀SnowFork測(cè)量3次取平均值得到[27]。

    表1 所用的GF-3影像的成像參數(shù)Tab.1 Parameters of the used GF-3 image

    雪深測(cè)量結(jié)果如圖2所示,可見(jiàn)研究區(qū)雪深空間分布不均,西南部戈壁區(qū)域積雪較淺,東北部山地區(qū)域積雪較深。這是因?yàn)閬?lái)自西南方向的暖濕氣團(tuán)在東移過(guò)程與阿爾泰山成近正交角,地形對(duì)氣流的抬升作用顯著,有利于山區(qū)降水,使得山麓的降水多于河谷平原地區(qū)[28]。為了探究模型在不同積雪條件下的適用性,將研究區(qū)分為深雪區(qū)與淺雪區(qū)兩個(gè)區(qū)域,分別進(jìn)行反演與分析。淺雪區(qū)的雪深觀測(cè)點(diǎn)16個(gè),平均雪深27.5 cm;深雪區(qū)的雪深觀測(cè)點(diǎn)36個(gè),平均雪深72.6 cm。

    深雪區(qū)觀測(cè)剖面5個(gè),平均積雪密度0.23 g/cm3,平均體積含水量0.64%;淺雪區(qū)觀測(cè)剖面3個(gè),平均積雪密度0.17 g/cm3,平均體積含水量0.27%。雪層含水量都接近于0,可認(rèn)為積雪為干雪,滿足CPD模型條件。分層積雪密度與液態(tài)水含量的觀測(cè)結(jié)果如圖3所示??梢钥闯鰷\雪區(qū)密度與含水量分層垂直差異小,深雪區(qū)分層積雪密度存在明顯的垂直差異,雪層密度呈中間大、底部和頂部小的特征。深雪區(qū)底層積雪的含水量高于表層積雪,垂直分布表現(xiàn)出較大的差異。圖4為研究區(qū)雪粒圖像,可以看出深雪區(qū)較淺雪區(qū)雪粒垂直分異更明顯,積雪重結(jié)晶變質(zhì)作用更強(qiáng)。

    圖1 研究區(qū)位置Fig.1 Location of the study area

    圖2 觀測(cè)雪深點(diǎn)的空間分布Fig.2 Position of the measured snow depth points

    圖3 積雪剖面觀測(cè)結(jié)果Fig.3 Distribution of snow depth,snow density,and snow wetness of the measured snow samples

    圖4 各雪層的雪粒觀測(cè)結(jié)果Fig.4 Observed results of snow particle of each layer

    (3)其他數(shù)據(jù)。根據(jù)阿勒泰國(guó)家基準(zhǔn)氣候站提供的氣溫和雪深數(shù)據(jù),在星地同步觀測(cè)前,1月13日前后發(fā)生了一次較大的降雪,自1月13日至2月底無(wú)較大降雪天氣過(guò)程,平均雪深約為14.5 cm。此外,采用國(guó)家基礎(chǔ)地理信息中心發(fā)布的全球30 m地表覆蓋數(shù)據(jù)GlobeLand 30掩膜森林等復(fù)雜散射機(jī)理的區(qū)域,使其不參與雪深反演;采用30 m分辨率的SRTM DEM數(shù)據(jù)作為數(shù)字高程數(shù)據(jù)。

    2 本文研究方法

    2.1 積雪各向異性結(jié)構(gòu)

    積雪的物理特性與其微觀幾何結(jié)構(gòu)有關(guān)。在理想狀態(tài)下,干雪可以看作是冰粒和空氣的二相混合物[29]。為了模擬雪的各向異性結(jié)構(gòu),假定干雪冰粒為規(guī)則的橢球體,以橢球體中心為原點(diǎn),建立三維笛卡兒坐標(biāo)系(x,y,z),z軸為冰粒晶體對(duì)稱(chēng)軸,平行于重力方向,x和y平面平行于平坦地表。ax、ay和az是橢球體在x、y、z方向上的半軸長(zhǎng),假設(shè)ax=ay,冰粒軸間比r=ax/az,表示冰粒的形狀。r=1時(shí),為各向同性隨機(jī)結(jié)構(gòu)的新降雪;r>1時(shí),積雪在自身重力壓實(shí)作用下迅速密實(shí)化[30],逐漸轉(zhuǎn)變?yōu)楦飨虍愋?、圓盤(pán)狀水平排列結(jié)構(gòu)的雪粒[31],r越大,各向異性越強(qiáng);r<1時(shí),積雪受溫度梯度作用發(fā)生重結(jié)晶和變質(zhì)作用,形成各向異性、針狀垂直排列的變質(zhì)雪粒[32-33],r越趨近于0,各向異性越強(qiáng)。由新降雪到發(fā)生重結(jié)晶與變質(zhì)作用變成陳雪的過(guò)程,伴隨冰粒晶體形狀及軸間比r的演變?nèi)鐖D5所示。

    圖5 積雪冰粒演化過(guò)程Fig.5 Evolution of snow particle shape

    2.2 同極化相位差模型

    當(dāng)電磁波穿透干雪層時(shí),積雪粒子的大小遠(yuǎn)小于電磁波波長(zhǎng),干雪可認(rèn)為是均質(zhì)的、體散射可忽略不計(jì)的介質(zhì)[34-35],因此后向散射信號(hào)主要來(lái)源于下墊面,因雪層折射率而延遲。文獻(xiàn)[36]基于Maxwell-Garnett理論將冰川冰的各向異性結(jié)構(gòu)與觀測(cè)到的負(fù)相位差聯(lián)系起來(lái),建立了CPD模型。文獻(xiàn)[37]進(jìn)一步基于干雪微結(jié)構(gòu)及其微波極化特性,將CPD模型與雪的微觀結(jié)構(gòu)聯(lián)系起來(lái)。

    如圖6所示,當(dāng)電磁波以θ入射角穿透一個(gè)深度為SD的干雪層,冰粒在電磁波極化電場(chǎng)作用下極化,依據(jù)Laplace方程及其邊界條件,計(jì)算出3個(gè)方向的退極化因子Ni,i∈{x,y,z}[38]

    (1)

    Nx=Ny=0.5(1-Nz)

    (2)

    圖6 極化信號(hào)穿透雪層產(chǎn)生相位差示意圖(改編自文獻(xiàn)[37])Fig.6 Electromagnetic wave penetrates the snow layer and results in the phase difference (adapted from reference[37])

    根據(jù)Maxwell-Garnett等效介質(zhì)理論計(jì)算冰粒各向異性介電常數(shù)εeff,i,i∈{x,y,z}

    (3)

    式中,Ni為退極化因子;εice為冰的等效介電常數(shù);εair為空氣的等效介電常數(shù);fvol為積雪與冰的密度比,故各向異性介電常數(shù)可由積雪密度ρsnow、冰粒軸間比r計(jì)算得到。

    因?yàn)檎凵渎蕁與介電常數(shù)εe存在εe=n2的關(guān)系,故HH極化和VV極化的折射率nH和nV與3個(gè)軸向上的等效介電常數(shù)εeff,i、入射角θ滿足以下關(guān)系

    (4)

    (5)

    HH和VV極化信號(hào)在雪層中的折射率不同,導(dǎo)致后向散射信號(hào)所傳播的路徑和距離不同,從而產(chǎn)生了相位差異。利用波長(zhǎng)與相位的關(guān)系,以及極化信號(hào)傳播路徑和雪層的幾何關(guān)系,得到不同雪深引起的同極化相位差

    (6)

    2.3 雪深反演方法

    根據(jù)式(6)可計(jì)算厚度為SD的積雪導(dǎo)致的相位差,由此可以推導(dǎo)雪深反演模型為

    (7)

    根據(jù)nV、nH與εeff,i的關(guān)系、εeff,i與積雪密度ρsnow、Ni的關(guān)系以及Ni與冰粒軸間比r的關(guān)系可知,在確定SAR數(shù)據(jù)計(jì)算CPD之后,輸入積雪密度和積雪冰粒軸間比參數(shù)即可得到雪深。

    本文研究在文獻(xiàn)[19]的基礎(chǔ)上進(jìn)一步深入,通過(guò)模擬CPD與雪深的關(guān)系可知,在假定積雪各向異性介電常數(shù)不變的情況下,CPD與雪深呈線性關(guān)系,如圖7所示。當(dāng)研究區(qū)有多個(gè)實(shí)測(cè)點(diǎn)時(shí),可根據(jù)實(shí)測(cè)點(diǎn)雪深數(shù)據(jù)與從SAR數(shù)據(jù)計(jì)算的CPD擬合得到最適合研究區(qū)積雪特性的線性模型經(jīng)驗(yàn)系數(shù)a、b,從而得到該區(qū)域全部像元的雪深。即

    (8)

    (9)

    式中,γc為絕對(duì)相干系數(shù);SVV和SHH分別為VV和HH極化單視復(fù)數(shù)影像;〈〉為濾波器。本文借助干涉雷達(dá)數(shù)據(jù)處理的軟件Gamma計(jì)算CPD,包括數(shù)據(jù)預(yù)處理、極化相干、空間濾波、地理編碼等步驟??紤]實(shí)測(cè)樣點(diǎn)的空間分布,設(shè)置地理編碼CPD的空間分辨率為10 m。

    圖7 積雪各向異性介電常恒定下CPD與雪深的關(guān)系Fig.7 Relationship between CPD and snow depth based on constant snow anisotropic relative permittivity

    3 反演結(jié)果分析

    3.1 反演精度分析

    深雪區(qū)與淺雪區(qū)積雪屬性差異大,區(qū)域內(nèi)部積雪屬性差異小。為獲得更精確的反演結(jié)果,分區(qū)域進(jìn)行雪深反演。對(duì)于同一實(shí)測(cè)點(diǎn),實(shí)測(cè)雪深可以看作一個(gè)常數(shù)值,而不同濾波窗口大小計(jì)算得到的CPD不同,則反演得到的擬合關(guān)系、反演精度不同。本文使用留P法進(jìn)行交叉驗(yàn)證,以有效避免過(guò)擬合和欠擬合狀態(tài),評(píng)價(jià)不同濾波窗口大小下深雪區(qū)、淺雪區(qū)的反演模型優(yōu)劣。為了充分利用所有實(shí)測(cè)點(diǎn)數(shù)據(jù),P取1,即每次假定一個(gè)點(diǎn)的雪深未知,使用其余點(diǎn)處實(shí)測(cè)雪深值與CPD擬合得到模型系數(shù),反演得到該點(diǎn)雪深,然后使用相關(guān)系數(shù)R和均方根誤差RMSE來(lái)分析所有點(diǎn)上雪深實(shí)測(cè)值與反演結(jié)果的誤差,如圖8所示。

    結(jié)果表明,淺雪區(qū)的最優(yōu)濾波窗口為59×59像元,相關(guān)系數(shù)R為0.83,均方根誤差RMSE為2.72 cm;深雪區(qū)的最優(yōu)濾波窗口為33×33像元,相關(guān)系數(shù)R為0.54,均方根誤差RMSE為11.69 cm。隨著濾波窗口大小的增加,淺雪區(qū)和深雪區(qū)的R均呈現(xiàn)出先上升再下降的趨勢(shì),RMSE呈現(xiàn)先大幅上漲波動(dòng)、然后減小、最后穩(wěn)定的趨勢(shì)。在濾波窗口25×25像元以下的結(jié)果,R與RMSE出現(xiàn)較大波動(dòng),這是因?yàn)闉V波窗口較小則不能有效抑制相干斑噪聲,相位的波動(dòng)使相干性降低;當(dāng)濾波窗口過(guò)大,則會(huì)平滑CPD的變化,削弱圖像的細(xì)節(jié),也導(dǎo)致相干性降低。

    圖8 反演結(jié)果R和RMSE與濾波窗口大小的關(guān)系Fig.8 Relationship between R,RMSE and the window size of filter

    不同區(qū)域最優(yōu)濾波窗口大小下雪深的空間分布如圖9所示。圖9(a)顯示深雪區(qū)反演得到的雪深自西向東逐漸減少,南部和北部的雪深小于中部區(qū)域。圖9(c)為淺雪區(qū)雪深反演的結(jié)果,可見(jiàn)除西北部反演雪深明顯大于其他區(qū)域以外,其他區(qū)域反演雪深呈相間分布。掩膜森林、水體等具有復(fù)雜的散射機(jī)制及受人類(lèi)活動(dòng)影響嚴(yán)重的人造地表區(qū)域的雪深結(jié)果如圖9(b)、(d)所示。

    3.2 反演不確定性分析

    地形通過(guò)影響水熱條件從而影響積雪空間分布,地形起伏帶來(lái)雷達(dá)構(gòu)像幾何關(guān)系的變化從而影響雷達(dá)散射機(jī)制,對(duì)雪深的反演結(jié)果產(chǎn)生影響[39]。為分析地形對(duì)反演精度的影響,選取坡度作為變量與反演誤差進(jìn)行相關(guān)性分析。如圖10所示,雪深反演誤差與坡度顯著相關(guān),隨著坡度的增加,雪深的反演誤差呈現(xiàn)增加趨勢(shì),且淺雪區(qū)的相關(guān)性強(qiáng)于深雪區(qū)。由淺雪區(qū)實(shí)測(cè)點(diǎn)坡度分布范圍可知,淺雪區(qū)的地形起伏較深雪區(qū)更大,由于構(gòu)建CPD模型時(shí)假定冰粒的z軸平行于重力方向,復(fù)雜的地形會(huì)導(dǎo)致極化信號(hào)穿透雪層的幾何關(guān)系發(fā)生改變,影響模型的適用性,同時(shí)電磁波回波由于受到地形起伏干擾作用會(huì)引入過(guò)多的噪聲,從而增大反演誤差。

    圖9 最優(yōu)濾波窗口下的雪深反演結(jié)果Fig.9 Spatial distribution ofretrieved snow depth with optimal window size of filter

    CPD模型以雪層均質(zhì)、忽略體散射為基礎(chǔ),然而嚴(yán)格來(lái)說(shuō),自然界中的媒介質(zhì)都是不均勻的,不可避免會(huì)引入誤差。由實(shí)測(cè)結(jié)果可知,深雪區(qū)積雪密度較淺雪區(qū)大,由圖4實(shí)測(cè)雪粒圖像也可以看出深雪區(qū)經(jīng)過(guò)重力壓實(shí)與變質(zhì)作用導(dǎo)致粒子粗糙度更大,雪層均質(zhì)性更差,極化信號(hào)穿透雪層時(shí)體散射增強(qiáng),此時(shí),后向散射信號(hào)不僅僅取決于下墊面,也由雪層中不同的散射體決定,有效散射中心向雪層上表面移動(dòng),因此同極化相干性受到體去相關(guān)的影響而降低[18];根據(jù)電磁波在介質(zhì)中的穿透深度公式[29],針對(duì)本文所用的GF-3數(shù)據(jù),對(duì)不同濕度下的積雪穿透深度進(jìn)行了模擬,結(jié)果如圖11所示??梢钥闯鯟波段穿透深度隨著積雪濕度增加而顯著下降。由圖3的剖面雪深與濕度的關(guān)系可以看出,由于深雪區(qū)較強(qiáng)的土壤熱傳導(dǎo)作用導(dǎo)致積雪底層含水量較大,部分底層含水量超過(guò)3%,水的介電常數(shù)在很大程度上改變了雪的介電特性,大幅降低了極化信號(hào)對(duì)雪層的穿透能力,導(dǎo)致極化信號(hào)不能到達(dá)下墊面,因此深雪區(qū)雪深反演精度低于淺雪區(qū)。

    圖10 反演誤差與坡度的關(guān)系Fig.10 Correlation between slope and inversion errors

    此外,為了探討不同入射角觀測(cè)幾何條件下CPD模型反演能力的差異,使用入射角發(fā)生改變時(shí)單位雪深引起的CPD變化來(lái)表征,單位雪深引起的CPD變化越大則代表CPD對(duì)該入射角下的電磁波越敏感,反演能力越強(qiáng)。模擬CPD與不同入射角之間的關(guān)系如圖12所示。CPD隨SAR電磁波入射角的增加而增加,單位雪深引起的CPD變化也隨著SAR電磁波入射的增加而增加,但二者是非線性關(guān)系。當(dāng)入射角增大時(shí),HH極化與VV極化的折射率差增大,所引起的CPD變化也增大,即CPD對(duì)入射角較大的SAR電磁波更加敏感。故不同入射角觀測(cè)幾何條件下,本文方法的反演能力存在差異,它更適合具有大入射角的SAR衛(wèi)星。

    圖11 C波段極化信號(hào)積雪穿透深度模擬Fig.11 Simulation the penetration depth of C band polarimetric signal

    圖12 模擬CPD與不同入射角之間的關(guān)系Fig.12 Simulation the relationship between CPD and incidence angle

    3.3 反演方法比較

    為驗(yàn)證本文所提出的反演方法的優(yōu)勢(shì),與已有方法進(jìn)行對(duì)比。根據(jù)文獻(xiàn)[19]提出的雪深反演方法,使用研究區(qū)平均積雪密度和假定軸間比1.5代入反演模型,求得雪深;根據(jù)文獻(xiàn)[20]提出的雪深反演方法,結(jié)合不同濾波窗口下的實(shí)測(cè)點(diǎn)處CPD、實(shí)測(cè)雪深、雪密度等參數(shù)求得軸間比,再結(jié)合該點(diǎn)處的實(shí)測(cè)積雪密度,輸入模型反演該點(diǎn)處的雪深。使用留P法(P=1)進(jìn)行交叉驗(yàn)證來(lái)評(píng)價(jià)3種反演方法優(yōu)劣,比較結(jié)果見(jiàn)表2。根據(jù)文獻(xiàn)[19]提出的雪深反演方法得到掩膜后的深雪區(qū)、淺雪區(qū)的雪深分布結(jié)果分別如圖13(a)、13(c)所示,根據(jù)文獻(xiàn)[20]提出的雪深反演方法得到掩膜后的深雪區(qū)、淺雪區(qū)的雪深分布結(jié)果分別如圖13(b)、13(d)所示。

    表2 本文方法與已有方法的雪深反演精度比較Tab.2 Comparison of the proposed method with the existing methods

    圖13 已有雪深反演方法的反演結(jié)果Fig.13 Spatial distribution of retrieved snow depth based on the Patil’s and Majumdar’s methods

    可以看出,對(duì)于本文研究所用GF-3數(shù)據(jù),本文方法比另外兩種方法精度更高,且將反演所需要的參數(shù)減少為遙感獲取的CPD數(shù)據(jù)和進(jìn)行模型訓(xùn)練的實(shí)測(cè)雪深數(shù)據(jù)。比較結(jié)果證明了本文方法具有較好的反演積雪深度的能力。

    4 結(jié) 論

    本文以新疆阿勒泰克蘭河上游地區(qū)為研究區(qū),在獲取C波段全極化GF-3數(shù)據(jù)及地面同步觀測(cè)數(shù)據(jù)的基礎(chǔ)上,基于CPD模型構(gòu)建線性的半經(jīng)驗(yàn)雪深反演模型,根據(jù)實(shí)測(cè)積雪屬性差異,分深雪區(qū)、淺雪區(qū)反演雪深,探討了反演方法在不同積雪條件下的適用性,對(duì)反演不確定性進(jìn)行了分析并與已有方法進(jìn)行比較,主要結(jié)論包括:

    (1)在積雪各向異性介電常數(shù)視為恒定的理想情況下,CPD僅是雪深的函數(shù),可根據(jù)全部實(shí)測(cè)點(diǎn)的雪深數(shù)據(jù)與從SAR數(shù)據(jù)獲取的CPD數(shù)據(jù),擬合得到最適合研究區(qū)積雪特性的模型經(jīng)驗(yàn)系數(shù),反演研究區(qū)雪深。反演精度的高低與計(jì)算CPD過(guò)程中使用的濾波器的窗口大小有關(guān),淺雪區(qū)的最優(yōu)濾波窗口為59×59像元,此時(shí)R為0.83,RMSE為2.72 cm,深雪區(qū)的最優(yōu)濾波窗口為33×33像元,此時(shí)R為0.54,RMSE為11.69 cm。

    (2)雪深反演誤差與坡度顯著相關(guān),隨著坡度的增加,雪深的反演誤差呈現(xiàn)增加的趨勢(shì),雪深反演不確定性受雪層變質(zhì)程度、含水量及衛(wèi)星入射角觀測(cè)幾何條件影響,反演方法對(duì)于干燥、雪層變質(zhì)結(jié)晶程度低、均質(zhì)的積雪,以及具有大入射角的SAR衛(wèi)星有更好的適用性。

    (3)對(duì)比已有基于CPD模型的雪深反演方法,本文所提出的反演方法精度更高,并且將反演所需要的參數(shù)減少為遙感獲取的CPD數(shù)據(jù)和進(jìn)行模型訓(xùn)練的實(shí)測(cè)雪深數(shù)據(jù),能夠可靠地反演積雪深度。

    猜你喜歡
    冰粒雪深積雪
    冰雹中也有“年輪”
    我們
    一種基于頻率與相位算法的雷達(dá)雪深探測(cè)方法
    高原冬季雪深與重慶夏季降水的年際關(guān)系研究
    大糧積雪 誰(shuí)解老將廉頗心
    炎黃地理(2017年10期)2018-01-31 02:15:05
    積雪
    2000~2014年西藏高原積雪覆蓋時(shí)空變化
    我國(guó)冰粒降水天氣的觀測(cè)特征統(tǒng)計(jì)分析
    鐵路防災(zāi)雪深圖像采集的設(shè)計(jì)和實(shí)現(xiàn)
    近30年西藏雪深時(shí)空變化及其對(duì)氣候變化的響應(yīng)
    国产在线男女| 国产麻豆成人av免费视频| 久久久久性生活片| 一本精品99久久精品77| 黄色配什么色好看| 日韩欧美国产一区二区入口| 天堂√8在线中文| 最后的刺客免费高清国语| 欧美3d第一页| 午夜日韩欧美国产| 亚洲 国产 在线| 国产精品人妻久久久影院| 日本 欧美在线| 国产精品久久久久久久久免| 人妻制服诱惑在线中文字幕| 国产三级中文精品| 搡老岳熟女国产| 日本 欧美在线| 亚洲无线在线观看| а√天堂www在线а√下载| 国产欧美日韩精品亚洲av| 露出奶头的视频| 日韩欧美免费精品| 欧美黑人巨大hd| 成人av在线播放网站| 一级黄片播放器| 婷婷丁香在线五月| av黄色大香蕉| 日韩大尺度精品在线看网址| 在线国产一区二区在线| 看免费成人av毛片| 精品一区二区三区视频在线观看免费| 国产黄a三级三级三级人| 给我免费播放毛片高清在线观看| 日韩大尺度精品在线看网址| 日韩欧美国产一区二区入口| 亚洲av第一区精品v没综合| 国产视频内射| 看片在线看免费视频| 日韩欧美 国产精品| 波多野结衣高清作品| 久久午夜亚洲精品久久| 又黄又爽又刺激的免费视频.| 小说图片视频综合网站| 日日撸夜夜添| 最近最新免费中文字幕在线| 老司机福利观看| 欧美色欧美亚洲另类二区| 此物有八面人人有两片| 日韩国内少妇激情av| 中文在线观看免费www的网站| 村上凉子中文字幕在线| 久久人人精品亚洲av| 搡老岳熟女国产| 亚州av有码| 国产黄a三级三级三级人| 制服丝袜大香蕉在线| 成人无遮挡网站| 麻豆成人午夜福利视频| 欧美丝袜亚洲另类 | 免费一级毛片在线播放高清视频| 亚洲无线观看免费| 免费av不卡在线播放| 91午夜精品亚洲一区二区三区 | 欧美中文日本在线观看视频| 国产伦在线观看视频一区| 午夜影院日韩av| 天美传媒精品一区二区| 欧美另类亚洲清纯唯美| 国产高清不卡午夜福利| 在线免费观看的www视频| 国产精华一区二区三区| 窝窝影院91人妻| 午夜福利高清视频| 乱人视频在线观看| 亚洲欧美精品综合久久99| 亚洲最大成人中文| 久久99热这里只有精品18| 亚洲五月天丁香| 一本精品99久久精品77| 一进一出抽搐gif免费好疼| 亚洲内射少妇av| 三级毛片av免费| 99热精品在线国产| 国产av在哪里看| 极品教师在线免费播放| 亚洲欧美精品综合久久99| 色综合色国产| 久久精品夜夜夜夜夜久久蜜豆| 女人十人毛片免费观看3o分钟| 一级黄片播放器| 国产伦在线观看视频一区| 亚洲精品成人久久久久久| 欧美zozozo另类| 国产精品久久电影中文字幕| 国产一区二区三区视频了| 国产黄a三级三级三级人| 午夜激情福利司机影院| 天天一区二区日本电影三级| 久久精品国产清高在天天线| 久99久视频精品免费| 色综合婷婷激情| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 又爽又黄a免费视频| 国产精品久久久久久久电影| 亚洲精品粉嫩美女一区| 国产亚洲欧美98| 亚洲最大成人中文| 色吧在线观看| 99在线视频只有这里精品首页| 99久久精品一区二区三区| 日韩欧美国产一区二区入口| 嫩草影院入口| 国产在线精品亚洲第一网站| 97碰自拍视频| 国产伦人伦偷精品视频| av黄色大香蕉| 悠悠久久av| 色综合站精品国产| 精品日产1卡2卡| 三级男女做爰猛烈吃奶摸视频| 丰满人妻一区二区三区视频av| 在现免费观看毛片| 搡老熟女国产l中国老女人| eeuss影院久久| 欧美成人一区二区免费高清观看| 麻豆一二三区av精品| 久久久午夜欧美精品| 69人妻影院| 黄色丝袜av网址大全| 亚洲一区二区三区色噜噜| 成人特级av手机在线观看| 啦啦啦啦在线视频资源| 久久香蕉精品热| 色哟哟·www| 两人在一起打扑克的视频| 一个人看视频在线观看www免费| 精品国产三级普通话版| 少妇猛男粗大的猛烈进出视频 | 亚洲午夜理论影院| 九色国产91popny在线| 国产伦一二天堂av在线观看| 国产爱豆传媒在线观看| 国内揄拍国产精品人妻在线| a级毛片免费高清观看在线播放| 男人的好看免费观看在线视频| 欧美成人免费av一区二区三区| 国产在线男女| 动漫黄色视频在线观看| 草草在线视频免费看| 一个人免费在线观看电影| x7x7x7水蜜桃| 国内久久婷婷六月综合欲色啪| 国产精品98久久久久久宅男小说| 有码 亚洲区| 欧美丝袜亚洲另类 | 最近中文字幕高清免费大全6 | 精华霜和精华液先用哪个| 97热精品久久久久久| 联通29元200g的流量卡| 18禁黄网站禁片免费观看直播| 成人无遮挡网站| 日日摸夜夜添夜夜添小说| 亚洲美女黄片视频| 亚洲欧美清纯卡通| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 日韩亚洲欧美综合| 又黄又爽又免费观看的视频| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久久久久久丰满 | 亚洲欧美清纯卡通| 色5月婷婷丁香| 99久久精品热视频| 国产日本99.免费观看| 国产精品98久久久久久宅男小说| 黄色视频,在线免费观看| 午夜老司机福利剧场| 精品久久久久久成人av| 伊人久久精品亚洲午夜| 午夜福利在线观看吧| 欧美zozozo另类| 女人被狂操c到高潮| 伦精品一区二区三区| 日本免费一区二区三区高清不卡| 九色成人免费人妻av| 俺也久久电影网| 毛片女人毛片| 51国产日韩欧美| 亚洲男人的天堂狠狠| 久久久久久久久久成人| 亚洲精品一区av在线观看| 亚洲av日韩精品久久久久久密| 亚洲四区av| 国产色爽女视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 精品久久久久久久久亚洲 | 久久久久精品国产欧美久久久| 白带黄色成豆腐渣| 少妇被粗大猛烈的视频| 久久午夜亚洲精品久久| 欧美人与善性xxx| 99热网站在线观看| xxxwww97欧美| 国产老妇女一区| 我要看日韩黄色一级片| 99热这里只有是精品50| av国产免费在线观看| 国产成人影院久久av| 少妇的逼水好多| eeuss影院久久| 精品欧美国产一区二区三| 91狼人影院| 国产免费av片在线观看野外av| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 欧美性猛交黑人性爽| 日韩欧美在线二视频| 一区二区三区高清视频在线| 最后的刺客免费高清国语| 中文资源天堂在线| 亚洲av成人av| 人妻久久中文字幕网| 我的老师免费观看完整版| 成人午夜高清在线视频| 亚洲午夜理论影院| 在线免费观看的www视频| 欧美3d第一页| 两个人的视频大全免费| 在线观看午夜福利视频| 一个人看视频在线观看www免费| 久久精品影院6| 亚洲欧美精品综合久久99| 一区二区三区高清视频在线| 免费观看的影片在线观看| 国产精品亚洲美女久久久| 亚洲av美国av| 国产精品久久久久久久久免| 欧美又色又爽又黄视频| 欧美色视频一区免费| 白带黄色成豆腐渣| 一个人免费在线观看电影| 亚洲三级黄色毛片| 99久久中文字幕三级久久日本| 女人被狂操c到高潮| 1000部很黄的大片| 一本久久中文字幕| h日本视频在线播放| 听说在线观看完整版免费高清| 久久久久久久亚洲中文字幕| 婷婷六月久久综合丁香| 91精品国产九色| 日本免费一区二区三区高清不卡| 成年人黄色毛片网站| 国产精品乱码一区二三区的特点| 嫩草影院精品99| 日本一本二区三区精品| 婷婷六月久久综合丁香| 女的被弄到高潮叫床怎么办 | 国产69精品久久久久777片| 不卡一级毛片| 狠狠狠狠99中文字幕| 可以在线观看毛片的网站| 国产精品久久久久久精品电影| 又爽又黄无遮挡网站| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久末码| 又黄又爽又免费观看的视频| 欧美性感艳星| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 久久人妻av系列| 亚洲成人中文字幕在线播放| 欧美3d第一页| 成人国产一区最新在线观看| 国产人妻一区二区三区在| 麻豆成人av在线观看| 国产精品女同一区二区软件 | 一a级毛片在线观看| 亚洲av美国av| 欧美最新免费一区二区三区| 久久人人精品亚洲av| 国产高清有码在线观看视频| 久久久久免费精品人妻一区二区| 九九热线精品视视频播放| av在线老鸭窝| 久久婷婷人人爽人人干人人爱| 两人在一起打扑克的视频| 一进一出抽搐gif免费好疼| 一区二区三区免费毛片| 黄色配什么色好看| 久久久精品欧美日韩精品| 淫妇啪啪啪对白视频| 美女免费视频网站| 干丝袜人妻中文字幕| 日韩一本色道免费dvd| 欧美xxxx黑人xx丫x性爽| 伊人久久精品亚洲午夜| 搡老妇女老女人老熟妇| 久久久久久久久久成人| av在线天堂中文字幕| 国产探花在线观看一区二区| 久久精品影院6| 老师上课跳d突然被开到最大视频| 少妇被粗大猛烈的视频| 成熟少妇高潮喷水视频| 一级a爱片免费观看的视频| 久久久久性生活片| 国产男靠女视频免费网站| 精品一区二区三区视频在线| 国产免费一级a男人的天堂| 欧美xxxx黑人xx丫x性爽| 欧美日韩瑟瑟在线播放| 免费大片18禁| 国产精品不卡视频一区二区| 久久久精品欧美日韩精品| 日韩一区二区视频免费看| 日韩欧美三级三区| 亚洲在线自拍视频| 亚洲自偷自拍三级| 日韩大尺度精品在线看网址| 99在线人妻在线中文字幕| 亚洲人成网站高清观看| 国产男人的电影天堂91| 天堂网av新在线| 丰满人妻一区二区三区视频av| 日本三级黄在线观看| 亚洲美女黄片视频| 午夜激情福利司机影院| 91久久精品电影网| 桃红色精品国产亚洲av| 在线a可以看的网站| 丝袜美腿在线中文| 人妻制服诱惑在线中文字幕| а√天堂www在线а√下载| 91在线观看av| a级毛片免费高清观看在线播放| 成人三级黄色视频| 亚洲四区av| 成人高潮视频无遮挡免费网站| 亚洲av日韩精品久久久久久密| 久久人妻av系列| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲欧美98| 亚洲一区二区三区色噜噜| 中国美女看黄片| 人妻制服诱惑在线中文字幕| 黄色视频,在线免费观看| 99精品在免费线老司机午夜| 精品午夜福利在线看| 亚洲成人久久爱视频| 天堂av国产一区二区熟女人妻| 一级黄片播放器| 亚洲精华国产精华液的使用体验 | 亚洲欧美日韩卡通动漫| 免费人成在线观看视频色| 永久网站在线| 久久精品综合一区二区三区| 蜜桃亚洲精品一区二区三区| 国产精品国产高清国产av| 亚洲一级一片aⅴ在线观看| 两人在一起打扑克的视频| 精品不卡国产一区二区三区| 成人特级av手机在线观看| 日韩欧美在线二视频| 亚洲精品亚洲一区二区| 日本 欧美在线| 99久久精品一区二区三区| 久久这里只有精品中国| 亚洲欧美激情综合另类| 中文亚洲av片在线观看爽| 国产精品精品国产色婷婷| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 国产亚洲av嫩草精品影院| 天天一区二区日本电影三级| 色播亚洲综合网| 老熟妇乱子伦视频在线观看| 亚洲精品久久国产高清桃花| 动漫黄色视频在线观看| 人妻少妇偷人精品九色| 欧美又色又爽又黄视频| 国产精品国产三级国产av玫瑰| 91精品国产九色| 欧美bdsm另类| 日韩欧美在线二视频| 黄色女人牲交| 校园春色视频在线观看| 久久精品国产清高在天天线| 国产色婷婷99| 久久久久久国产a免费观看| 中文字幕av在线有码专区| 露出奶头的视频| 在线观看舔阴道视频| 制服丝袜大香蕉在线| 日本熟妇午夜| 丝袜美腿在线中文| 久久香蕉精品热| 国产亚洲av嫩草精品影院| 国产女主播在线喷水免费视频网站 | 国产私拍福利视频在线观看| 成人美女网站在线观看视频| 久久精品综合一区二区三区| 欧美+日韩+精品| 岛国在线免费视频观看| 一级黄色大片毛片| 久久久色成人| 在线观看av片永久免费下载| 日本黄色片子视频| 少妇丰满av| 99久久成人亚洲精品观看| 国产精品免费一区二区三区在线| 免费在线观看日本一区| 国产精品精品国产色婷婷| 伊人久久精品亚洲午夜| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 亚洲人与动物交配视频| 久久人人精品亚洲av| 成熟少妇高潮喷水视频| 久久中文看片网| 中文字幕久久专区| or卡值多少钱| 看免费成人av毛片| 超碰av人人做人人爽久久| 国产欧美日韩精品一区二区| 精品久久久久久,| 久久精品夜夜夜夜夜久久蜜豆| 三级毛片av免费| 成人毛片a级毛片在线播放| 一夜夜www| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 乱码一卡2卡4卡精品| 村上凉子中文字幕在线| 午夜爱爱视频在线播放| 国产毛片a区久久久久| 舔av片在线| 国产精品野战在线观看| 免费看美女性在线毛片视频| 欧美不卡视频在线免费观看| 色精品久久人妻99蜜桃| 国产一级毛片七仙女欲春2| 国产精品美女特级片免费视频播放器| 亚洲美女黄片视频| 精品久久久久久成人av| 波多野结衣高清作品| 在线国产一区二区在线| 69av精品久久久久久| 一级黄色大片毛片| 国产亚洲av嫩草精品影院| 波野结衣二区三区在线| 国产精品女同一区二区软件 | 三级毛片av免费| 亚洲色图av天堂| 日本免费a在线| 黄色配什么色好看| videossex国产| 身体一侧抽搐| 99热6这里只有精品| 五月伊人婷婷丁香| 午夜久久久久精精品| 亚洲18禁久久av| 亚洲成av人片在线播放无| 男女之事视频高清在线观看| 一级黄色大片毛片| 国产人妻一区二区三区在| 国产午夜精品久久久久久一区二区三区 | 成熟少妇高潮喷水视频| 成人国产麻豆网| 久久久久久久久久久丰满 | 成人亚洲精品av一区二区| 在线观看一区二区三区| 亚洲天堂国产精品一区在线| 九九热线精品视视频播放| 国产伦一二天堂av在线观看| 久久天躁狠狠躁夜夜2o2o| 成年女人看的毛片在线观看| 午夜爱爱视频在线播放| 日韩中字成人| 中文字幕人妻熟人妻熟丝袜美| 天堂√8在线中文| 在线国产一区二区在线| а√天堂www在线а√下载| 又紧又爽又黄一区二区| 国产欧美日韩一区二区精品| 亚洲乱码一区二区免费版| 特大巨黑吊av在线直播| 婷婷六月久久综合丁香| 99精品在免费线老司机午夜| 偷拍熟女少妇极品色| 网址你懂的国产日韩在线| 高清毛片免费观看视频网站| 我的老师免费观看完整版| 亚洲成a人片在线一区二区| 一级a爱片免费观看的视频| 久久人妻av系列| 女生性感内裤真人,穿戴方法视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久午夜福利片| 国产女主播在线喷水免费视频网站 | 欧美三级亚洲精品| 久久久久性生活片| 日韩精品有码人妻一区| 久久精品国产鲁丝片午夜精品 | 免费人成在线观看视频色| 我要搜黄色片| 网址你懂的国产日韩在线| 欧美黑人巨大hd| 亚洲成人免费电影在线观看| 一级黄片播放器| avwww免费| 国产又黄又爽又无遮挡在线| 18禁在线播放成人免费| h日本视频在线播放| 亚洲成人久久爱视频| 少妇人妻一区二区三区视频| 国产综合懂色| 舔av片在线| 成年版毛片免费区| 乱系列少妇在线播放| 亚洲欧美日韩无卡精品| 又黄又爽又刺激的免费视频.| 国产成人福利小说| 熟妇人妻久久中文字幕3abv| 午夜免费激情av| 午夜福利在线在线| 国内久久婷婷六月综合欲色啪| 长腿黑丝高跟| 少妇高潮的动态图| 久久天躁狠狠躁夜夜2o2o| 此物有八面人人有两片| 热99re8久久精品国产| 亚洲真实伦在线观看| 久久精品国产亚洲网站| 日本色播在线视频| 亚洲在线观看片| 国产欧美日韩精品亚洲av| 窝窝影院91人妻| 欧美精品国产亚洲| 成人欧美大片| 国产高清激情床上av| 哪里可以看免费的av片| 久久久久精品国产欧美久久久| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 国产真实伦视频高清在线观看 | 三级男女做爰猛烈吃奶摸视频| 午夜福利在线观看吧| 婷婷精品国产亚洲av| 两个人视频免费观看高清| 欧美绝顶高潮抽搐喷水| 极品教师在线免费播放| 亚洲在线观看片| 美女黄网站色视频| 亚洲国产精品久久男人天堂| 国产一区二区在线观看日韩| 久久久久久久久久久丰满 | 91久久精品国产一区二区成人| 成人高潮视频无遮挡免费网站| 久久久久久久精品吃奶| 国产乱人伦免费视频| 97热精品久久久久久| 国产视频一区二区在线看| 免费观看的影片在线观看| 国产一区二区三区视频了| 国内精品宾馆在线| 日本黄大片高清| 国产麻豆成人av免费视频| 波多野结衣高清无吗| 99热6这里只有精品| 婷婷色综合大香蕉| 日本三级黄在线观看| 国产男人的电影天堂91| 欧美不卡视频在线免费观看| 99精品久久久久人妻精品| 久久精品国产亚洲av香蕉五月| 国产激情偷乱视频一区二区| 欧美成人一区二区免费高清观看| 亚洲精品亚洲一区二区| 国产在视频线在精品| 国产精品精品国产色婷婷| 在线观看免费视频日本深夜| 男人狂女人下面高潮的视频| 亚洲熟妇中文字幕五十中出| 精品福利观看| 狂野欧美白嫩少妇大欣赏| 我的女老师完整版在线观看| 中文资源天堂在线| 最近在线观看免费完整版| 久久精品国产自在天天线| 麻豆一二三区av精品| 亚洲五月天丁香| 亚洲自拍偷在线| 国产探花在线观看一区二区| 精品久久久久久久末码| 免费看a级黄色片| 免费观看精品视频网站| 99久久无色码亚洲精品果冻| 亚洲精品一区av在线观看| 啦啦啦观看免费观看视频高清| 欧美激情在线99| 精品人妻1区二区| 老女人水多毛片| а√天堂www在线а√下载| 午夜福利高清视频| 国产大屁股一区二区在线视频| 美女xxoo啪啪120秒动态图| 99久久久亚洲精品蜜臀av| 精品久久久久久成人av| 欧美性感艳星| 午夜福利在线观看免费完整高清在 | 波多野结衣高清作品| 白带黄色成豆腐渣|