• <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)
    久久人妻熟女aⅴ| 黄色毛片三级朝国网站 | 麻豆精品久久久久久蜜桃| 日日啪夜夜爽| 精品少妇黑人巨大在线播放| 日韩人妻高清精品专区| 国产亚洲最大av| 最后的刺客免费高清国语| 亚洲国产欧美日韩在线播放 | 啦啦啦视频在线资源免费观看| 三级国产精品片| 国产伦精品一区二区三区四那| 精品熟女少妇av免费看| 大码成人一级视频| 热re99久久国产66热| 亚洲国产色片| 国产精品久久久久久精品古装| 亚洲三级黄色毛片| 插逼视频在线观看| 国产熟女午夜一区二区三区 | 丰满乱子伦码专区| 国产免费视频播放在线视频| 黑人巨大精品欧美一区二区蜜桃 | 久久婷婷青草| 久久久亚洲精品成人影院| 777米奇影视久久| 一本久久精品| 国产日韩欧美在线精品| 男的添女的下面高潮视频| 国产一区二区在线观看av| 久久 成人 亚洲| 亚洲av综合色区一区| 在线播放无遮挡| 亚洲精品乱码久久久久久按摩| 日韩av免费高清视频| 国内精品宾馆在线| 免费黄频网站在线观看国产| h日本视频在线播放| 看非洲黑人一级黄片| 黄色怎么调成土黄色| 国产精品国产三级国产专区5o| 人妻 亚洲 视频| 三级国产精品欧美在线观看| 一本—道久久a久久精品蜜桃钙片| 99九九在线精品视频 | 高清av免费在线| 精品久久久精品久久久| 国产色爽女视频免费观看| 美女脱内裤让男人舔精品视频| 少妇的逼好多水| 校园人妻丝袜中文字幕| 国语对白做爰xxxⅹ性视频网站| 精品少妇内射三级| 欧美精品一区二区大全| 男女边摸边吃奶| 亚洲精品乱码久久久久久按摩| 精品国产一区二区久久| av福利片在线| 边亲边吃奶的免费视频| 伊人久久国产一区二区| 国内精品宾馆在线| 97精品久久久久久久久久精品| 特大巨黑吊av在线直播| 80岁老熟妇乱子伦牲交| 波野结衣二区三区在线| 久久99一区二区三区| 在线观看av片永久免费下载| 在线观看免费日韩欧美大片 | 成人毛片60女人毛片免费| 日韩中字成人| 亚洲内射少妇av| 一二三四中文在线观看免费高清| 国产片特级美女逼逼视频| 国产乱人偷精品视频| 国产在线一区二区三区精| 波野结衣二区三区在线| 又粗又硬又长又爽又黄的视频| 中文乱码字字幕精品一区二区三区| 夜夜爽夜夜爽视频| 日日啪夜夜撸| 美女xxoo啪啪120秒动态图| 欧美丝袜亚洲另类| 免费av不卡在线播放| 精品一区二区免费观看| 国产精品久久久久久精品电影小说| 亚洲三级黄色毛片| 街头女战士在线观看网站| 天堂8中文在线网| 国产精品伦人一区二区| 老司机影院毛片| 亚洲国产精品专区欧美| 国产精品久久久久久精品电影小说| 人妻 亚洲 视频| 汤姆久久久久久久影院中文字幕| 卡戴珊不雅视频在线播放| 丝瓜视频免费看黄片| 国产av码专区亚洲av| 黑人高潮一二区| 在现免费观看毛片| 美女中出高潮动态图| 天天躁夜夜躁狠狠久久av| 91久久精品国产一区二区三区| 少妇熟女欧美另类| 国产精品.久久久| 欧美bdsm另类| 午夜免费鲁丝| 久久6这里有精品| 亚洲av国产av综合av卡| 深夜a级毛片| 哪个播放器可以免费观看大片| 看十八女毛片水多多多| 成人二区视频| 免费少妇av软件| 成人综合一区亚洲| 欧美高清成人免费视频www| 亚洲欧美一区二区三区黑人 | av天堂久久9| 王馨瑶露胸无遮挡在线观看| 日韩免费高清中文字幕av| 日韩 亚洲 欧美在线| 新久久久久国产一级毛片| 18+在线观看网站| 久久这里有精品视频免费| 另类亚洲欧美激情| 国产精品一区二区性色av| 哪个播放器可以免费观看大片| 一级a做视频免费观看| 六月丁香七月| 日韩中文字幕视频在线看片| 国产伦精品一区二区三区视频9| 亚洲精品aⅴ在线观看| 汤姆久久久久久久影院中文字幕| 国产亚洲91精品色在线| 男男h啪啪无遮挡| 亚洲国产成人一精品久久久| 99re6热这里在线精品视频| 天天躁夜夜躁狠狠久久av| 丰满人妻一区二区三区视频av| 国产精品久久久久久av不卡| 五月天丁香电影| 精品一区二区三区视频在线| 丝袜喷水一区| 99久久人妻综合| 80岁老熟妇乱子伦牲交| 天天躁夜夜躁狠狠久久av| 亚洲欧洲日产国产| 亚洲精品国产av成人精品| 又爽又黄a免费视频| 亚洲成色77777| 欧美日韩在线观看h| 国产欧美另类精品又又久久亚洲欧美| 在线观看国产h片| av不卡在线播放| 精品少妇内射三级| 国产一级毛片在线| 国产伦在线观看视频一区| 熟女电影av网| 美女国产视频在线观看| 亚洲国产精品一区二区三区在线| 日产精品乱码卡一卡2卡三| 七月丁香在线播放| av天堂久久9| 啦啦啦视频在线资源免费观看| 欧美日韩av久久| 乱系列少妇在线播放| 国产黄色视频一区二区在线观看| 亚洲欧美日韩另类电影网站| 亚洲国产精品国产精品| 国产成人一区二区在线| 欧美老熟妇乱子伦牲交| 久久久久久久久久久免费av| 国产 一区精品| 色婷婷av一区二区三区视频| 久久青草综合色| 丝袜在线中文字幕| 亚洲欧美一区二区三区国产| 国产成人精品婷婷| 成年人午夜在线观看视频| 日韩一区二区视频免费看| 五月开心婷婷网| 日韩伦理黄色片| av国产精品久久久久影院| 国产在线一区二区三区精| 色网站视频免费| 亚洲经典国产精华液单| 看免费成人av毛片| 婷婷色综合大香蕉| 狂野欧美激情性xxxx在线观看| 国产乱来视频区| 亚洲久久久国产精品| 你懂的网址亚洲精品在线观看| 免费观看在线日韩| 亚洲欧洲国产日韩| 国产精品女同一区二区软件| 色视频在线一区二区三区| 少妇 在线观看| 国产精品人妻久久久久久| 美女大奶头黄色视频| 久久久精品94久久精品| 建设人人有责人人尽责人人享有的| 久久久欧美国产精品| 亚洲精品乱码久久久久久按摩| 国产精品一区二区在线不卡| 久久久久国产精品人妻一区二区| 久热这里只有精品99| 国产永久视频网站| 亚洲美女视频黄频| 欧美日韩视频高清一区二区三区二| 人妻人人澡人人爽人人| 中文字幕免费在线视频6| 波野结衣二区三区在线| 国内少妇人妻偷人精品xxx网站| 亚洲内射少妇av| 热99国产精品久久久久久7| 新久久久久国产一级毛片| 人妻制服诱惑在线中文字幕| 亚洲精品日韩在线中文字幕| 精品一区二区免费观看| 精品久久久久久久久亚洲| 多毛熟女@视频| 啦啦啦中文免费视频观看日本| 国产日韩欧美视频二区| 日韩电影二区| 日韩电影二区| 高清黄色对白视频在线免费看 | 国产av国产精品国产| 欧美丝袜亚洲另类| 日韩熟女老妇一区二区性免费视频| 亚洲怡红院男人天堂| 夜夜骑夜夜射夜夜干| 国产成人aa在线观看| 国产精品免费大片| 国产中年淑女户外野战色| 国产69精品久久久久777片| 欧美变态另类bdsm刘玥| 国产黄色免费在线视频| 青青草视频在线视频观看| 狠狠精品人妻久久久久久综合| 高清午夜精品一区二区三区| 国产亚洲欧美精品永久| 丝袜脚勾引网站| 久热这里只有精品99| 国产精品成人在线| 亚洲精品国产成人久久av| 国产精品99久久99久久久不卡 | 春色校园在线视频观看| 中文字幕人妻熟人妻熟丝袜美| h日本视频在线播放| 日韩精品有码人妻一区| 嫩草影院入口| 精品人妻一区二区三区麻豆| 色网站视频免费| 亚洲综合色惰| 中文字幕精品免费在线观看视频 | 99热全是精品| 成人综合一区亚洲| 成人亚洲欧美一区二区av| 国产精品成人在线| 国精品久久久久久国模美| 国产乱人偷精品视频| 精品国产一区二区三区久久久樱花| 男女边摸边吃奶| 国产av码专区亚洲av| 国产极品粉嫩免费观看在线 | 国产精品人妻久久久久久| 久久精品久久久久久久性| 国产毛片在线视频| 日本vs欧美在线观看视频 | 久久久久久久久久久丰满| 久热久热在线精品观看| 日本黄色日本黄色录像| 草草在线视频免费看| 一区二区三区乱码不卡18| 建设人人有责人人尽责人人享有的| 精品视频人人做人人爽| 69精品国产乱码久久久| 免费观看的影片在线观看| 成年女人在线观看亚洲视频| 亚洲精品亚洲一区二区| 久久99热这里只频精品6学生| 国产深夜福利视频在线观看| 欧美三级亚洲精品| 高清av免费在线| 精品人妻熟女av久视频| 午夜视频国产福利| 日韩中字成人| av卡一久久| 日本91视频免费播放| 亚洲精品自拍成人| 久久ye,这里只有精品| 一二三四中文在线观看免费高清| 国产片特级美女逼逼视频| 国产精品一二三区在线看| 精品视频人人做人人爽| 大话2 男鬼变身卡| 97超视频在线观看视频| 下体分泌物呈黄色| 秋霞伦理黄片| 国产美女午夜福利| 久久毛片免费看一区二区三区| 国产一区二区在线观看日韩| 国模一区二区三区四区视频| 午夜视频国产福利| 五月伊人婷婷丁香| 日本欧美国产在线视频| 美女中出高潮动态图| 汤姆久久久久久久影院中文字幕| 久久青草综合色| 一级毛片aaaaaa免费看小| 欧美bdsm另类| 亚洲国产精品成人久久小说| 亚洲精品中文字幕在线视频 | 久久av网站| av网站免费在线观看视频| 观看美女的网站| 秋霞在线观看毛片| 51国产日韩欧美| 欧美老熟妇乱子伦牲交| 亚洲欧洲国产日韩| 亚洲欧美清纯卡通| 青春草视频在线免费观看| 久久热精品热| 9色porny在线观看| 亚洲精品国产av蜜桃| 国产av精品麻豆| 国产欧美另类精品又又久久亚洲欧美| 热99国产精品久久久久久7| 男女无遮挡免费网站观看| 久久久久久久久久久久大奶| 激情五月婷婷亚洲| 亚洲精品中文字幕在线视频 | 久久99一区二区三区| 天堂8中文在线网| 午夜精品国产一区二区电影| 天天躁夜夜躁狠狠久久av| 亚洲成人手机| 中文资源天堂在线| av免费在线看不卡| 免费在线观看成人毛片| av.在线天堂| 色哟哟·www| 一区二区三区乱码不卡18| 日韩在线高清观看一区二区三区| 日本vs欧美在线观看视频 | 免费观看a级毛片全部| 一区二区三区免费毛片| 97超碰精品成人国产| 一本色道久久久久久精品综合| 国产一级毛片在线| 亚洲av国产av综合av卡| 老司机亚洲免费影院| 日韩成人伦理影院| 国产 一区精品| 中文乱码字字幕精品一区二区三区| 高清视频免费观看一区二区| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 亚洲国产精品成人久久小说| 久久精品国产亚洲av涩爱| 十分钟在线观看高清视频www | 亚洲人成网站在线观看播放| 亚洲伊人久久精品综合| 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 热re99久久国产66热| 另类亚洲欧美激情| 亚洲精品色激情综合| 亚洲精品中文字幕在线视频 | 9色porny在线观看| av线在线观看网站| 一区二区三区精品91| 国产成人a∨麻豆精品| 欧美+日韩+精品| 日本猛色少妇xxxxx猛交久久| 中文在线观看免费www的网站| 国产 一区精品| 国产美女午夜福利| 欧美区成人在线视频| 嫩草影院新地址| 亚洲第一av免费看| 精品国产一区二区久久| 少妇 在线观看| 波野结衣二区三区在线| 午夜免费鲁丝| av专区在线播放| 最近最新中文字幕免费大全7| 下体分泌物呈黄色| 成人亚洲欧美一区二区av| 天美传媒精品一区二区| 在线观看人妻少妇| 哪个播放器可以免费观看大片| 日本欧美国产在线视频| 国产成人freesex在线| 亚洲自偷自拍三级| 丰满饥渴人妻一区二区三| 夫妻性生交免费视频一级片| 国产国拍精品亚洲av在线观看| 晚上一个人看的免费电影| 一级毛片黄色毛片免费观看视频| 视频中文字幕在线观看| 乱人伦中国视频| 久久久亚洲精品成人影院| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线 | 边亲边吃奶的免费视频| 视频中文字幕在线观看| 美女中出高潮动态图| 亚洲国产精品专区欧美| 国产精品伦人一区二区| 国产精品.久久久| 欧美激情极品国产一区二区三区 | 日日摸夜夜添夜夜爱| 日本欧美国产在线视频| 一本久久精品| 内射极品少妇av片p| 搡老乐熟女国产| 日韩视频在线欧美| 免费大片黄手机在线观看| 秋霞伦理黄片| 男女边摸边吃奶| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 国产亚洲午夜精品一区二区久久| 国产深夜福利视频在线观看| av网站免费在线观看视频| 美女xxoo啪啪120秒动态图| 久久久a久久爽久久v久久| 自拍偷自拍亚洲精品老妇| 色婷婷久久久亚洲欧美| 十八禁网站网址无遮挡 | 简卡轻食公司| 最近中文字幕2019免费版| 一边亲一边摸免费视频| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 亚洲av男天堂| 国产精品一二三区在线看| 男女免费视频国产| 99久久中文字幕三级久久日本| av在线观看视频网站免费| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 十分钟在线观看高清视频www | 国产伦精品一区二区三区四那| 国产精品人妻久久久影院| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 欧美精品亚洲一区二区| 最新的欧美精品一区二区| 国产中年淑女户外野战色| 亚洲第一区二区三区不卡| 久久人妻熟女aⅴ| 国产av国产精品国产| 精品久久久噜噜| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| av播播在线观看一区| 国产美女午夜福利| 永久网站在线| 99九九在线精品视频 | 9色porny在线观看| kizo精华| 国产日韩一区二区三区精品不卡 | 国产精品伦人一区二区| 国产精品国产三级专区第一集| 乱系列少妇在线播放| 国产女主播在线喷水免费视频网站| 日本爱情动作片www.在线观看| 久久久国产精品麻豆| av卡一久久| 高清黄色对白视频在线免费看 | 搡女人真爽免费视频火全软件| 久久久久久久亚洲中文字幕| 中国国产av一级| 三级经典国产精品| 精品一区二区三区视频在线| 91aial.com中文字幕在线观看| 亚洲精品乱码久久久v下载方式| xxx大片免费视频| 国产精品人妻久久久久久| 最近中文字幕高清免费大全6| 国产av码专区亚洲av| a 毛片基地| 亚洲美女黄色视频免费看| 亚洲av在线观看美女高潮| 国产成人精品婷婷| 久久精品国产亚洲av涩爱| 一区二区三区乱码不卡18| 91午夜精品亚洲一区二区三区| 中文字幕亚洲精品专区| 国产成人精品婷婷| 伊人久久精品亚洲午夜| 久久6这里有精品| 女人久久www免费人成看片| av国产精品久久久久影院| 亚洲精品456在线播放app| 精品久久久噜噜| 亚洲内射少妇av| 免费人成在线观看视频色| 99精国产麻豆久久婷婷| 我的老师免费观看完整版| 中文字幕精品免费在线观看视频 | 插阴视频在线观看视频| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 伊人久久精品亚洲午夜| 在线观看免费高清a一片| av天堂中文字幕网| 99久国产av精品国产电影| 成人毛片60女人毛片免费| 国产精品一区www在线观看| 亚洲精品自拍成人| 日本黄色片子视频| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| 黑人高潮一二区| 日韩电影二区| 久久午夜综合久久蜜桃| 欧美+日韩+精品| 少妇猛男粗大的猛烈进出视频| 人体艺术视频欧美日本| 亚洲av免费高清在线观看| 寂寞人妻少妇视频99o| 亚洲欧美中文字幕日韩二区| 男女边摸边吃奶| 国产美女午夜福利| 最新的欧美精品一区二区| 国产精品一二三区在线看| 在线观看免费日韩欧美大片 | 国产精品久久久久久久久免| 日韩av在线免费看完整版不卡| 久久人人爽人人片av| 亚洲熟女精品中文字幕| 自线自在国产av| 国产精品人妻久久久影院| 午夜视频国产福利| 久久青草综合色| 九九在线视频观看精品| 丝袜脚勾引网站| 免费高清在线观看视频在线观看| 91精品国产九色| 七月丁香在线播放| 夜夜爽夜夜爽视频| 精品久久国产蜜桃| 伦精品一区二区三区| 久久6这里有精品| 精品国产露脸久久av麻豆| 下体分泌物呈黄色| 亚洲真实伦在线观看| 日本黄色片子视频| 日韩伦理黄色片| 国产日韩欧美亚洲二区| 少妇的逼好多水| 欧美+日韩+精品| 精品国产国语对白av| 两个人的视频大全免费| 久久精品久久久久久久性| 国产精品秋霞免费鲁丝片| 国产日韩欧美在线精品| 国产视频首页在线观看| 国产伦在线观看视频一区| 毛片一级片免费看久久久久| 国产免费福利视频在线观看| 亚洲精品第二区| 99热国产这里只有精品6| 国产色爽女视频免费观看| 色哟哟·www| 国产精品99久久久久久久久| 午夜福利视频精品| 亚洲av日韩在线播放| 国产日韩欧美视频二区| 日韩中字成人| 在线观看人妻少妇| √禁漫天堂资源中文www| 人妻一区二区av| 少妇裸体淫交视频免费看高清| 欧美日韩综合久久久久久| 成年人午夜在线观看视频| 2022亚洲国产成人精品| 18禁动态无遮挡网站| 免费不卡的大黄色大毛片视频在线观看| 久久久精品免费免费高清| 校园人妻丝袜中文字幕| 人体艺术视频欧美日本| 亚洲欧美中文字幕日韩二区| 婷婷色麻豆天堂久久| 婷婷色av中文字幕| 国产成人91sexporn| 黄色欧美视频在线观看| 久久午夜综合久久蜜桃| 大香蕉久久网| 插阴视频在线观看视频| 一级,二级,三级黄色视频| 国产精品蜜桃在线观看| 美女国产视频在线观看| 久久久久久久亚洲中文字幕| 国产亚洲最大av| 久久精品国产亚洲av涩爱| 免费观看无遮挡的男女| 熟女电影av网| 少妇人妻精品综合一区二区| 国产真实伦视频高清在线观看| 丝瓜视频免费看黄片| 黑人高潮一二区| 在线播放无遮挡| 国产免费视频播放在线视频| av.在线天堂| 日韩一区二区三区影片| 欧美日韩精品成人综合77777| av黄色大香蕉| 多毛熟女@视频| 久久鲁丝午夜福利片| 午夜福利影视在线免费观看| 国产永久视频网站| 亚洲精品中文字幕在线视频 |