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

    基于二維小波變換的亮溫降尺度方法研究

    2022-08-08 05:45:32何勤業(yè)盧涵宇袁詠儀盧天健
    關(guān)鍵詞:亮溫土壤水分波段

    何勤業(yè),盧涵宇*,袁詠儀,盧天健

    (1.貴州大學(xué) 大數(shù)據(jù)與信息工程學(xué)院, 貴州 貴陽(yáng) 550025;2.貴州六盤水三力達(dá)科技有限公司, 貴州 六盤水 532001)

    0 引言

    土壤水分是氣候系統(tǒng)中的一個(gè)關(guān)鍵狀態(tài)變量,對(duì)理解陸地-大氣能量、水分和碳循環(huán)通量具有重要意義[1]。土壤水分在洪水和干旱預(yù)報(bào)、氣候變化、天氣預(yù)報(bào)、水資源管理和農(nóng)業(yè)生產(chǎn)等方面有著廣泛的應(yīng)用,對(duì)人類社會(huì)和環(huán)境有著重要的影響[2],因此,精確度和高空間分辨率的土壤水分全球測(cè)繪對(duì)于滿足這些應(yīng)用的需求至關(guān)重要。

    目前最廣泛應(yīng)用的測(cè)量地表土壤水分的方法是被動(dòng)微波遙感,它能全天候全天時(shí)觀測(cè)數(shù)據(jù)。由于微波遙感的1.41 GHz(L波段)具有更好的穿透大氣和植被覆蓋能力及受射頻干擾(radio grequency interference,RFI)影響較小特點(diǎn),因此在大多數(shù)與土壤水分相關(guān)的研究以該頻段為主要研究對(duì)象。第一個(gè)搭載L波段微波輻射計(jì)的衛(wèi)星是在2009年11月2日由歐空局發(fā)射的土壤水分和海洋鹽度(soil moisture and ocean salinity,SMOS)衛(wèi)星。該衛(wèi)星可以提供地表的雙極化和大范圍入射角(5°~60°)的亮溫?cái)?shù)據(jù),同時(shí)這也是第一顆監(jiān)測(cè)全球土壤水分和海水鹽度的衛(wèi)星[3];但是L波段輻射計(jì)的空間分辨率只有40 km,并不足以滿足水文氣象、生態(tài)學(xué)、水資源管理等應(yīng)用超過空間分辨率為10 km的使用需求,所以美國(guó)國(guó)家航空航天局(national aeronautics and space administration,NASA)于2015年1月31日在加利福尼亞州發(fā)射土壤水分主動(dòng)-被動(dòng)(soil moisture active passive,SMAP)任務(wù)[4]。該衛(wèi)星搭載著空間分辨率為3 km的L波段雷達(dá)與空間分辨率為40 km的L波段輻射計(jì),可以通過幾種主被動(dòng)結(jié)合的降尺度方法,提供空間分辨率為9 km的全球土壤水分[4],但是,在2015年7月7日SMAP雷達(dá)故障,無法提供數(shù)據(jù)。而為了能繼續(xù)提供高分辨率數(shù)據(jù),SMAP使用了Backus-Gilbert(BG)optimal interpolation對(duì)原始的SMAP Level 1B亮溫產(chǎn)品(L1B_TB)進(jìn)行插值處理從而得到空間分辨率為9 km的亮溫?cái)?shù)據(jù)[5]。

    目前土壤水分的降尺度方法可以大致分為基于經(jīng)驗(yàn)?zāi)P偷姆椒?、基于?shù)理模型的方法和基于微波遙感數(shù)據(jù)融合的方法。其中基于經(jīng)驗(yàn)?zāi)P偷姆椒ㄒ全@取的高分辨率的地表參數(shù),再利用獲取的地表參數(shù)建立與微波土壤水分的模型,并將模型參數(shù)用于微波土壤水分的降尺度處理;但降尺度結(jié)果的準(zhǔn)確程度主要取決于此類經(jīng)驗(yàn)?zāi)P偷木_性,而且為了能更好的構(gòu)建經(jīng)驗(yàn)?zāi)P瓦€需要獲取大量的實(shí)地?cái)?shù)據(jù),因此限制了經(jīng)驗(yàn)?zāi)P驮谳^大區(qū)域內(nèi)的應(yīng)用[6]。而基于數(shù)理模型的方法主要是通過數(shù)理統(tǒng)計(jì)模型(如譜分析、多重分形和地學(xué)統(tǒng)計(jì)等空間插值方法)對(duì)微波數(shù)據(jù)進(jìn)行降尺度處理,其缺陷在于這些數(shù)理模型往往只是對(duì)單一的微波數(shù)據(jù)進(jìn)行空間上的處理,缺少額外的輔助數(shù)據(jù)作為支撐,并不能確保降尺度的精確程度[7]。本文使用的是基于微波遙感數(shù)據(jù)融合的方法,這種方法是通過將2種不同分辨率的微波數(shù)據(jù)進(jìn)行融合以達(dá)到降尺度的目的。其中文獻(xiàn)[8]利用傅里葉變換的方法將主被動(dòng)微波亮溫?cái)?shù)據(jù)進(jìn)行結(jié)合,但是SMAP的主動(dòng)微波雷達(dá)已經(jīng)損壞,而且傅里葉變換在分析微波信號(hào)這類非平穩(wěn)信號(hào)時(shí)存在局限性,會(huì)導(dǎo)致降尺度結(jié)果出現(xiàn)較大的誤差。文獻(xiàn)[9]使用時(shí)間序列觀測(cè)回歸分析降尺度法研究被動(dòng)微波中L波段與S波段亮溫的線性關(guān)系,將S波段作為輔助數(shù)據(jù)對(duì)L波段的亮溫進(jìn)行降尺度處理,但是目前還沒有搭載S波段的衛(wèi)星,而且S波段的分辨率較低遠(yuǎn)遠(yuǎn)不及10 km,無法滿足土壤水分的應(yīng)用需求。本文在上述2種方法的基礎(chǔ)上進(jìn)行了改進(jìn),首先是在選用微波遙感數(shù)據(jù)時(shí)選用了被動(dòng)微波的X波段作為L(zhǎng)波段亮溫降尺度的輔助波段,這是因?yàn)閄波段的分辨率為10 km,可以滿足土壤水分的應(yīng)用需求,而且X波段的亮溫?cái)?shù)據(jù)可以直接通過AMSR2得到;其次在對(duì)X波段與L波段進(jìn)行融合時(shí)使用了二維離散小波變換(2-D DWT)的方法,這種方法在分析非平穩(wěn)信號(hào)時(shí)沒有傅里葉分析的局限,可以更加精準(zhǔn)的對(duì)微波信號(hào)進(jìn)行分析與變換。具體步驟是通過2-D DWT分辨率為9 km的X波段亮溫與分辨率為36 km的L波段亮溫進(jìn)行分解與融合,從而得到降尺度后的L波段亮溫?cái)?shù)據(jù)。實(shí)驗(yàn)結(jié)果表明降尺度后的L波段亮溫可以滿足反演土壤水分的需求,且反演的土壤水分的精度與SMAP 9 km亮溫反演的土壤水分近似,具有極高的適用性。

    本文主要通過結(jié)合L波段與X波段亮溫的特性來實(shí)現(xiàn)亮溫降尺度的研究,并且將小波變換技術(shù)引入到了亮溫降尺度和土壤水分反演這2個(gè)領(lǐng)域,從而拓展了亮溫降尺度和土壤水分反演的研究方向。

    1 實(shí)驗(yàn)區(qū)域與數(shù)據(jù)

    本研究使用的L波段的亮溫?cái)?shù)據(jù)采用了SMAP-3級(jí)輻射計(jì)(SPL3SMP版本7)中的亮溫?cái)?shù)據(jù),SMAP數(shù)據(jù)于“美國(guó)國(guó)家冰雪數(shù)據(jù)中心”(https://nsidc.org/data)進(jìn)行下載。X波段的亮溫?cái)?shù)據(jù)則使用了AMSR2中10.65 GHz的L3亮溫?cái)?shù)據(jù),其數(shù)據(jù)來源于“日本宇宙航空研究開發(fā)機(jī)構(gòu)”(https://gportal.jaxa.jp/gpr/),并且選擇了位于帕里站點(diǎn)獲取的2016年1月1日到12月31日的土壤水分?jǐn)?shù)據(jù)作為驗(yàn)證數(shù)據(jù),該數(shù)據(jù)來源于“國(guó)家青藏高原科學(xué)數(shù)據(jù)中心”(http://data.tpdc.ac.cn)。

    1.1 實(shí)驗(yàn)區(qū)域

    本次研究區(qū)域選擇為帕里區(qū)域。該區(qū)域的土地覆蓋以草地、裸露地面為主,地形平坦,平均年降水量小于400 mm,且降雨大多集中在5—10月。在帕里區(qū)域內(nèi)一共部署了25個(gè)觀測(cè)站點(diǎn),在本研究中只用到了其中的8個(gè)觀測(cè)站點(diǎn),并分為了A、B共2個(gè)9 km網(wǎng)格,其中黑色三角點(diǎn)表示土壤水分的觀測(cè)站點(diǎn)[10]。帕里研究區(qū)土壤水分采樣點(diǎn)分布如圖1所示。

    圖1 帕里研究區(qū)土壤水分采樣點(diǎn)分布圖

    1.2 SMAP數(shù)據(jù)

    土壤水分主動(dòng)被動(dòng)(SMAP)衛(wèi)星于2015年1月31日發(fā)射,運(yùn)行周期為2~3 d,上升軌道時(shí)間為18:00,下降軌道時(shí)間為上午06:00[4]。SMAP衛(wèi)星分別搭載了L波段(1.41 GHz)輻射計(jì)和L波段(1.26 GHz)雷達(dá),入射角度為40°[4],可以提供空間分辨率為9、36 km的亮溫?cái)?shù)據(jù)。而空間分辨率為9 km的亮溫?cái)?shù)據(jù)在L波段雷達(dá)故障后就無法獲取,因此基于對(duì)BG 最佳插值的應(yīng)用,將SMAP的原始觀測(cè)數(shù)據(jù)進(jìn)行了修正從而獲得了空間分辨率為9 km的亮溫?cái)?shù)據(jù)。

    在本研究中使用了SMAP衛(wèi)星在帕里實(shí)驗(yàn)區(qū)域的分辨率36、9 km的亮溫?cái)?shù)據(jù),以及用亮溫反演土壤水分需要的輔助數(shù)據(jù),為了方便小波變換的應(yīng)用選取的范圍為16×16的9 km網(wǎng)格,具體的經(jīng)緯度為27.259 5°~28.534 2°、88.879 7°~90.373 4°。

    1.3 AMSR2數(shù)據(jù)

    GCOM-W由Japan Aerospace Exploration Agency(JAXA)在2012年5月18日發(fā)射,軌道位于地球上方約700 km處,上升軌道時(shí)間為13:30,下降軌道時(shí)間為01:30[11]。GCOM-W衛(wèi)星攜帶的無源微波輻射計(jì)Advanced Microwave Scanning Radiometer 2(AMSR2)是由AMSR-E改進(jìn)而來,AMSR2的天線旋轉(zhuǎn)周期為1.5 s[11],可以獲取1 450 km范圍內(nèi)的數(shù)據(jù),使得AMSR2能以2 d為一周期獲取全球99%以上的數(shù)據(jù)[11]。AMSR2的入射角為55°,有6個(gè)不同頻率的微波波段(6.9、10.7、18.7、23.8、36.5、89 GHz),可以提供空間分辨率為10、25 km的雙極化亮度溫度。本研究中使用了10.65 GHz(X波段)采集的空間分辨率為10 km亮度溫度作為降尺度的輔助數(shù)據(jù)。

    1.4 現(xiàn)場(chǎng)測(cè)量

    土壤水分是評(píng)價(jià)亮溫降尺度結(jié)果的一個(gè)重要標(biāo)準(zhǔn),因此本研究使用了帕里區(qū)域中按照布設(shè)的8個(gè)觀測(cè)站點(diǎn)(如圖1黑色三角點(diǎn)所示)獲取的土壤水分。同時(shí)為了方便進(jìn)行土壤水分的對(duì)比驗(yàn)證還對(duì)這8個(gè)觀測(cè)站點(diǎn)根據(jù)其所在的網(wǎng)格進(jìn)行了平均計(jì)算,從而得到這2個(gè)9 km網(wǎng)格的土壤水分。

    2 研究方法

    雖然L波段對(duì)于土壤水分含量的探測(cè)更加精準(zhǔn),但是SMAP的L波段的空間分辨率比AMSR2的X波段低,因此為了提高L波段的空間分辨率,使得L波段更加泛用,本文通過二維離散小波變換(2-D DWT)將L波段與X波段的亮溫?cái)?shù)據(jù)融合,從而得到的同時(shí)擁有L波段亮溫?cái)?shù)據(jù)特征和X波段空間分辨率的亮溫?cái)?shù)據(jù)。同時(shí)為了對(duì)這種降尺度方法的結(jié)果進(jìn)行對(duì)比還加入了SMAP 9 km亮溫?cái)?shù)據(jù)(通過Backus-Gilbert(BG)optimal interpolation獲取)。為了方便將降尺度后的亮溫?cái)?shù)據(jù)與SMAP的分辨率為9 km的亮溫?cái)?shù)據(jù)進(jìn)行對(duì)比,還需要對(duì)AMSR2的X波段亮溫?cái)?shù)據(jù)進(jìn)行插值,得到X波段分辨率為9 km的亮溫?cái)?shù)據(jù)。

    2.1 Backus-Gilbert(BG)最佳插值

    Backus-Gilbert(BG)最佳插值基于Backus-Gilbert理論,是一種對(duì)成像微波輻射計(jì)數(shù)據(jù)進(jìn)行插值的技術(shù),因?yàn)楸3峙c微波輻射計(jì)相關(guān)的天線增益函數(shù)的空間分辨率,所以其插值過程稱為最優(yōu)的。SMAP 9 km的亮溫?cái)?shù)據(jù)是通過對(duì)SMAP 1B級(jí)亮溫產(chǎn)品(L1B_TB)中的天線溫度(TA)的測(cè)量值進(jìn)行Backus-Gilbert(BG)最佳插值得到的,其插值生成的TA數(shù)據(jù)在經(jīng)過校準(zhǔn)后,就能在9 km網(wǎng)格上生成SMAP 1C級(jí)亮溫產(chǎn)品(L1C_TB_E)。

    2.2 二維離散小波變換(2-D DWT)

    小波變換作為一種分析工具,可以對(duì)不同空間分辨率下的亮溫?cái)?shù)據(jù)、土壤水分?jǐn)?shù)據(jù)等地球物理分布變量數(shù)據(jù)進(jìn)行特性變化的分析。

    連續(xù)小波變換(CWT)可描述為:假設(shè)存在一個(gè)能量有限的信號(hào)空間為L(zhǎng)2(R),且f(t)∈L2(R),可以用小波基將f(t)展開[12],得到f(t)的連續(xù)小波變換Wf(a,τ),

    (1)

    式中:a是伸縮因子;τ是平移因子;函數(shù)ψa,τ(t)為小波基函數(shù),其逆變換為

    (2)

    離散小波變換(DWT)是通過對(duì)CWT進(jìn)行二進(jìn)離散化實(shí)現(xiàn)的,可以將信號(hào)分解低頻與高頻2個(gè)部分到離散小波基的方法一般是將平移因子a和伸縮因子τ變成冪級(jí)數(shù)結(jié)構(gòu)[13],即

    (3)

    式中:a0≠1;τ0是常數(shù);j,k∈Z對(duì)應(yīng)的離散小波基為

    (4)

    最終f(t)的離散小波變換為

    (5)

    2-D DWT就是將二維數(shù)據(jù)分別變成行數(shù)據(jù)與列數(shù)據(jù)后再進(jìn)行一維離散小波變換。二維離散小波變換主要通過Mallat前向金字塔算法的擴(kuò)展實(shí)現(xiàn),二維小波變換原理如圖2所示。每進(jìn)行一次一維小波變換,原始數(shù)據(jù)就會(huì)分成高頻與低頻2個(gè)分量,因此二維數(shù)據(jù)在進(jìn)行一次二維小波變換后,原始數(shù)據(jù)將會(huì)被分成4個(gè)分量,分別是原始二維數(shù)據(jù)的低頻分量LL和高頻分量HL、LH、HH。

    圖2 二維小波變換原理圖

    其中HL包含表示原始數(shù)據(jù)在水平方向上的高頻分量,LH表示原始數(shù)據(jù)在垂直方向上的高頻分量,HH表示原始數(shù)據(jù)在對(duì)角處的高頻分量,這些高頻分量表示的是原始信號(hào)的邊緣信息,而LL表示原始數(shù)據(jù)的低頻分量,是原始數(shù)據(jù)的中心部分,代表了原始數(shù)據(jù)的主要信息。而為了能夠獲得更加精確的小波系數(shù),可以將第一次二維小波變換后獲得的近似值(低頻分量)LL作為下一次二維小波變換的輸入數(shù)據(jù),最多可進(jìn)行二維小波變換p=log2N次,N為原始數(shù)據(jù)的長(zhǎng)度。在實(shí)際應(yīng)用中,分解的級(jí)數(shù)也與小波基的選取有關(guān)。

    本研究使用了db5小波作為小波基,并進(jìn)行了2次2-D DWT,最后又使用了小波逆變換重構(gòu)出了新的亮溫?cái)?shù)據(jù)。之所以使用db5小波作為小波基是因?yàn)閐b5小波具有較好的正則性,且消失矩較大,比較光滑,具有較好的頻域局化能力,同時(shí)還不會(huì)使得時(shí)域的緊支撐性過低從而增大計(jì)算量使得實(shí)時(shí)性變差。而進(jìn)行2次二維小波變換的原因是為了在不失真的情況下獲得更準(zhǔn)確的小波系數(shù)。其中二維小波的重構(gòu)過程為:先對(duì)二維小波變換后獲得的高低頻分量的每一列進(jìn)行一維小波逆變換,再對(duì)逆變換后的結(jié)果的每一行進(jìn)行一維小波逆變換,即可得到重構(gòu)的數(shù)據(jù)。

    本研究的技術(shù)路線如圖3所示。

    圖3 技術(shù)路線

    2.3 土壤水分反演

    由于沒有實(shí)地觀測(cè)的分辨率為9 km亮溫?cái)?shù)據(jù),無法對(duì)降尺度后的L波段的亮溫?cái)?shù)據(jù)進(jìn)行直接有效的驗(yàn)證。因此用單通道檢索算法(single-channel algorithm,SCA)和地表輔助數(shù)據(jù)[14]。該算法主要是利用對(duì)土壤水分變換敏感的亮溫和輔助數(shù)據(jù)(地表溫度、地表土壤質(zhì)地、地表覆蓋情況、植被含水量)來實(shí)現(xiàn)對(duì)土壤水分的反演。

    SCA算法主要是基于零階近似模型ω-τ構(gòu)建的,其模型為

    (6)

    其實(shí)現(xiàn)過程具體分為以下5個(gè)步驟:

    (7)

    ② 通過地表覆蓋分類圖確定對(duì)應(yīng)的植被類型和參數(shù)b。

    ③ 根據(jù)地表覆蓋分類圖確定的植被類型獲取植被單次反照率ω和植被光學(xué)厚度τ,通過公式(8)、(9)計(jì)算出植被透過率γP和植被含水量VWC,用來校正植被影響

    (8)

    τ=b·VWC。

    (9)

    ④ 在將植被的影響去除后,使用Choudhury模型[15]對(duì)土壤粗糙度進(jìn)行校正,獲取光滑表面的發(fā)射率。

    ⑤ 最后用菲涅爾方程計(jì)算土壤介電常數(shù)后,使用Dobson模型[16]估算出土壤水分含量。

    該算法因?yàn)槭褂昧舜罅康妮o助數(shù)據(jù)來去除土壤溫度、大氣溫度、植被和地表粗糙度的影響,所以反演得到的土壤水分?jǐn)?shù)據(jù)精度很高,同時(shí)也是目前SMAP衛(wèi)星的官方的土壤水分反演算法。在本研究中的土壤水分反演用到的地表輔助數(shù)據(jù)均由SMAP數(shù)據(jù)集得到,并根據(jù)實(shí)驗(yàn)區(qū)域內(nèi)的植被分布選取單次散射反照率ω和植被參數(shù)b,且b=0.13,ω=0.05。在將亮溫?cái)?shù)據(jù)反演得到土壤水分?jǐn)?shù)據(jù)后,就可以與實(shí)地探測(cè)土壤水分?jǐn)?shù)據(jù)進(jìn)行對(duì)比,根據(jù)土壤水分?jǐn)?shù)據(jù)的對(duì)比的結(jié)果就可以確定降尺度后亮溫?cái)?shù)據(jù)的精度。

    3 實(shí)驗(yàn)結(jié)果與分析

    本次實(shí)驗(yàn)為了在獲取分辨率為9 km的L波段亮溫?cái)?shù)據(jù)的同時(shí)證明加入輔助數(shù)據(jù)對(duì)L波段亮溫?cái)?shù)據(jù)降尺度處理的重要性,以及驗(yàn)證使用2-D DWT對(duì)亮溫降尺度的有效性,因此計(jì)劃借助AMSR2的X波段亮溫?cái)?shù)據(jù)作為輔助數(shù)據(jù),將SMAP 36 km的亮溫?cái)?shù)據(jù)降尺度到9 km,并與SMAP 9 km的亮溫?cái)?shù)據(jù)進(jìn)行比較。由于AMSR2的X波段只有10 km的亮溫?cái)?shù)據(jù),因此需要使用線性插值的方法將其改變?yōu)? km的亮溫?cái)?shù)據(jù),以此滿足小波變換降尺度方法的要求。

    本次研究選取的時(shí)間段為2016-01-01—12-31日,共計(jì)366 d;但是,由于SMAP和AMSR2的運(yùn)行周期不一樣,因此為了保證在實(shí)驗(yàn)區(qū)域內(nèi)SMAP和AMSR2都能提供亮溫?cái)?shù)據(jù),在上午時(shí)間段選取了55 d作為研究樣本。

    3.1 亮溫降尺度結(jié)果

    本研究中使用了二維離散小波變換對(duì)L波段亮溫進(jìn)行了降尺度,其降尺度的亮溫圖像如圖4所示,顯示了在2016-06-22—09-26,期間,以32 d為間隔選取了4 d的L波段亮溫在上午時(shí)間段V極化和H極化的降尺度結(jié)果。其中第一行是SMAP 36 km原始的L波段亮溫圖像,第二行是AMSR2的9 km X波段亮溫圖像,第三行是二維離散小波變換降尺度的L波段亮溫圖像, 第四行是為了驗(yàn)證算法加入的SMAP 9 km的L波段亮溫圖像。通過比較可以發(fā)現(xiàn),第三行的降尺度圖像基本上可以呈現(xiàn)出SMAP 36 km亮溫圖像的空間分布情況與AMSR2 9 km亮溫圖像的空間結(jié)構(gòu),并與SMAP 9 km的亮溫圖像相比,二維離散小波變換的降尺度圖像的異質(zhì)性更高,與原始觀測(cè)的36 km亮溫圖像的關(guān)系更緊密,空間上的細(xì)節(jié)更多。

    (a)V極化

    其亮溫降尺度的圖像分析結(jié)果見表1,分別計(jì)算了SMAP 9 km的亮溫圖像與二維離散小波變換的降尺度亮溫圖像在不同極化下的信息熵、平均梯度、標(biāo)準(zhǔn)差、空間頻率。結(jié)果表明通過2-D DWT得來的降尺度圖像在不同極化的條件計(jì)算得到的標(biāo)準(zhǔn)差和空間頻率均高于SMAP 9 km的亮溫圖像,說明通過2-D DWT降尺度得到的亮溫圖像的有用信息更多,圖像空間的總體活躍度更高。而在信息熵和平均梯度的計(jì)算上2-D DWT圖像和SMAP 9 km圖像各有優(yōu)劣,總體持平,說明這2種亮溫圖像在平均信息量和圖像對(duì)比度上有著近似的特征。這些情況表明二維離散小波變換的亮溫圖像有著更多的細(xì)節(jié)表達(dá)能力,圖像信息更多,圖像空間的活躍程度更高。

    表1 亮溫降尺度圖像分析結(jié)果

    3.2 土壤水分降尺度結(jié)果和驗(yàn)證

    為了驗(yàn)證降尺度的結(jié)果,需要利用SCA算法與SMAP提供的輔助數(shù)據(jù),對(duì)2-D DWT降尺度后的9 km亮溫?cái)?shù)據(jù)和SMAP提供的9 km亮溫?cái)?shù)據(jù)進(jìn)行表面土壤水分反演。

    反演的土壤水分和站點(diǎn)土壤水分的長(zhǎng)時(shí)間序列圖與亮溫反演土壤水分的散點(diǎn)圖如圖5所示。其中區(qū)域A和區(qū)域B分別是選取的2個(gè)不同的觀測(cè)區(qū)域。為了確保土壤水分反演的精度,只反演了地表溫度大于0 ℃的時(shí)期。從長(zhǎng)時(shí)間序列圖中可知這2種不同的亮溫反演的土壤水分與站點(diǎn)觀測(cè)的土壤水分有很高的相似度,特別是V極化下反演的土壤水分相似度更高。從散點(diǎn)圖可知,用降尺度后的亮溫反演的土壤水分和用SMAP 9 km的亮溫反演的土壤水分極為相似,相關(guān)性最高為0.984 8,最低為0.970 3,表明兩者之間相似度極高,而且偏差值最高僅為0.005 6,表明兩者差值極低。其均方根誤差最高為0.008 3 cm3/cm3,表明兩者間的離散度非常低。以上情況說明用2-D DWT進(jìn)行降尺度的亮溫反演的土壤水分與SMAP 9 km亮溫反演的土壤水分有著極高的相似度。

    (a)V極化

    反演土壤水分與站點(diǎn)土壤水分的對(duì)比統(tǒng)計(jì)結(jié)果見表2。將2種亮溫在不同區(qū)域、不同極化的情況下反演的土壤水與相對(duì)應(yīng)的站點(diǎn)觀測(cè)的土壤水進(jìn)行了對(duì)比,計(jì)算了反演的土壤水分與站點(diǎn)土壤水分的相關(guān)系數(shù)(R)、偏差值(BIAS)、均方根誤差(RMSE)、無偏均方根誤差(ubRMSE)。比較結(jié)果表明,反演的土壤水分偏差值均在0.2以內(nèi),且相關(guān)系數(shù)均大于0.68,均方根誤差在0.2 cm3/cm3以內(nèi),無偏均方根誤差最高僅為0.021 9 cm3/cm3,表明反演的土壤水分與站點(diǎn)觀測(cè)的土壤水分的有著相似的性能,且在V極化情況下能達(dá)到最優(yōu)。

    表2 反演土壤水分與站點(diǎn)土壤水分對(duì)比統(tǒng)計(jì)結(jié)果

    4 結(jié)語(yǔ)

    本文通過二維離散小波變換借助X波段的高分辨率亮溫對(duì)L波段亮溫進(jìn)行降尺度處理,并與SMAP 9 km的亮溫?cái)?shù)據(jù)進(jìn)行了對(duì)比。其結(jié)果表明通過二維離散小波變換得來的亮溫圖像在整體上要優(yōu)于SMAP 9 km的亮溫圖像,且與SMAP 36 km的亮溫圖像的聯(lián)系更加緊密,并且在反演的土壤水分的結(jié)果對(duì)比中兩者結(jié)果十分近似,且都與站點(diǎn)觀測(cè)的土壤水分有很高的相關(guān)度,無偏均方根誤差都在0.02左右。這些數(shù)據(jù)說明在對(duì)L波段亮溫進(jìn)行降尺度時(shí),可以使用二維離散小波變換對(duì)L波段亮溫進(jìn)行降尺度處理,而且精度與SMAP 9 km的亮溫?cái)?shù)據(jù)相比略有提高。本文的研究只限于帕里區(qū)域,且在反演土壤水分時(shí)只在2個(gè)像元上進(jìn)行土壤水分反演結(jié)果驗(yàn)證,因此今后的工作需要將該方法改進(jìn),爭(zhēng)取應(yīng)用到全球區(qū)域的亮溫降尺度,為亮溫降尺度方法提供一種新的參考思路。

    猜你喜歡
    亮溫土壤水分波段
    春日暖陽(yáng)
    霰譜分布特征對(duì)強(qiáng)對(duì)流云高頻微波亮溫影響的模擬研究
    基于南太平洋的AMSR2 L1R亮溫?cái)?shù)據(jù)質(zhì)量評(píng)估
    西藏高原土壤水分遙感監(jiān)測(cè)方法研究
    M87的多波段輻射過程及其能譜擬合
    2014年2月12日新疆于田MS7.3地震熱紅外亮溫異常分析
    日常維護(hù)對(duì)L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    不同覆蓋措施對(duì)棗園土壤水分和溫度的影響
    L波段雷達(dá)磁控管的使用與維護(hù)
    河南科技(2014年18期)2014-02-27 14:14:53
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    亚洲成人国产一区在线观看| 久久久久久久久免费视频了| 国产精华一区二区三区| 深夜精品福利| 嫩草影视91久久| 91九色精品人成在线观看| 欧美在线一区亚洲| 欧美成人性av电影在线观看| 国产av不卡久久| 宅男免费午夜| 国产av一区二区精品久久| 午夜a级毛片| 精品国产一区二区三区四区第35| 国产精品久久久av美女十八| 2021天堂中文幕一二区在线观 | 午夜久久久久精精品| 亚洲天堂国产精品一区在线| 人成视频在线观看免费观看| 成人永久免费在线观看视频| 又黄又粗又硬又大视频| 久久婷婷成人综合色麻豆| 亚洲欧美精品综合一区二区三区| 色哟哟哟哟哟哟| 在线视频色国产色| 色在线成人网| 成人免费观看视频高清| 亚洲国产欧美一区二区综合| 久久香蕉激情| 国产不卡一卡二| 亚洲欧美日韩高清在线视频| 黄色毛片三级朝国网站| 欧美日本视频| 亚洲国产精品合色在线| 精品第一国产精品| 亚洲av熟女| 狠狠狠狠99中文字幕| 免费人成视频x8x8入口观看| 国产成人精品久久二区二区91| 日本a在线网址| 国产成人精品无人区| 高清毛片免费观看视频网站| 国产成人影院久久av| 高清毛片免费观看视频网站| a在线观看视频网站| 午夜免费激情av| 99riav亚洲国产免费| 一本久久中文字幕| 18禁裸乳无遮挡免费网站照片 | 大香蕉久久成人网| 高清毛片免费观看视频网站| 中文字幕最新亚洲高清| 欧美不卡视频在线免费观看 | 一进一出抽搐gif免费好疼| 午夜免费成人在线视频| 听说在线观看完整版免费高清| 国产亚洲av嫩草精品影院| 亚洲欧美精品综合久久99| 国产99久久九九免费精品| 曰老女人黄片| 1024香蕉在线观看| 亚洲精品久久国产高清桃花| 精品日产1卡2卡| 亚洲国产精品合色在线| 国产男靠女视频免费网站| 日韩欧美三级三区| 国产成人影院久久av| 婷婷丁香在线五月| 天天一区二区日本电影三级| 免费看a级黄色片| 久久久久久久久久黄片| tocl精华| 中文在线观看免费www的网站 | 国产1区2区3区精品| 变态另类成人亚洲欧美熟女| 香蕉av资源在线| 亚洲精品久久国产高清桃花| 精品一区二区三区av网在线观看| ponron亚洲| 亚洲熟妇熟女久久| 18禁黄网站禁片午夜丰满| 日韩大尺度精品在线看网址| 人人妻,人人澡人人爽秒播| 欧美性猛交╳xxx乱大交人| xxxwww97欧美| 不卡av一区二区三区| 美女大奶头视频| 国产高清有码在线观看视频 | 久久欧美精品欧美久久欧美| 99热只有精品国产| 亚洲精品粉嫩美女一区| 国产成年人精品一区二区| 亚洲精品色激情综合| 亚洲狠狠婷婷综合久久图片| 国内毛片毛片毛片毛片毛片| 亚洲av成人av| 国产三级在线视频| 可以在线观看毛片的网站| 中文字幕人妻熟女乱码| 香蕉国产在线看| 制服诱惑二区| 人人妻人人澡人人看| 久久久久国产精品人妻aⅴ院| 99热这里只有精品一区 | 久久国产亚洲av麻豆专区| 精品第一国产精品| 久久青草综合色| 两个人视频免费观看高清| 午夜福利高清视频| 久久午夜综合久久蜜桃| 十八禁人妻一区二区| 午夜免费激情av| av视频在线观看入口| videosex国产| 成人亚洲精品一区在线观看| 91成年电影在线观看| 叶爱在线成人免费视频播放| 一进一出好大好爽视频| 十八禁网站免费在线| 国产av不卡久久| 亚洲国产看品久久| 国产高清激情床上av| 一二三四社区在线视频社区8| 欧美性猛交╳xxx乱大交人| 日本 欧美在线| 欧美色视频一区免费| 美女 人体艺术 gogo| 精品第一国产精品| 黄色 视频免费看| 1024香蕉在线观看| 18禁裸乳无遮挡免费网站照片 | 一本一本综合久久| 亚洲男人的天堂狠狠| 国产伦一二天堂av在线观看| 国产又爽黄色视频| 宅男免费午夜| 欧美成人一区二区免费高清观看 | 久热这里只有精品99| 两个人视频免费观看高清| 久久草成人影院| 久久 成人 亚洲| 悠悠久久av| 亚洲成人国产一区在线观看| 国产精华一区二区三区| 啦啦啦观看免费观看视频高清| av免费在线观看网站| 国产一区在线观看成人免费| 18禁美女被吸乳视频| 99热6这里只有精品| 国产在线精品亚洲第一网站| 国产高清视频在线播放一区| 一二三四社区在线视频社区8| 亚洲九九香蕉| aaaaa片日本免费| 亚洲天堂国产精品一区在线| 麻豆久久精品国产亚洲av| 久久久久久大精品| av片东京热男人的天堂| 国产av一区在线观看免费| 久久久国产成人精品二区| 午夜激情福利司机影院| 国产成人精品久久二区二区免费| 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| 嫁个100分男人电影在线观看| xxxwww97欧美| 男女那种视频在线观看| xxx96com| 欧美黑人精品巨大| 国产一级毛片七仙女欲春2 | 亚洲一区二区三区色噜噜| 亚洲国产毛片av蜜桃av| 久久中文看片网| 俄罗斯特黄特色一大片| 成人免费观看视频高清| 亚洲av第一区精品v没综合| 日韩高清综合在线| av片东京热男人的天堂| 91大片在线观看| 真人一进一出gif抽搐免费| 国产又黄又爽又无遮挡在线| 搡老熟女国产l中国老女人| 精品一区二区三区av网在线观看| 成人永久免费在线观看视频| 日韩精品免费视频一区二区三区| 久久久久九九精品影院| 亚洲国产日韩欧美精品在线观看 | 欧美激情 高清一区二区三区| 亚洲精华国产精华精| 国产成人啪精品午夜网站| 亚洲最大成人中文| 亚洲精品美女久久久久99蜜臀| 欧美日韩福利视频一区二区| 香蕉久久夜色| 草草在线视频免费看| 久久精品成人免费网站| 午夜久久久久精精品| 久久国产精品影院| 天堂√8在线中文| 久99久视频精品免费| 欧美黄色淫秽网站| 女性被躁到高潮视频| 久久婷婷人人爽人人干人人爱| 麻豆av在线久日| 人人澡人人妻人| 黑人操中国人逼视频| 男女做爰动态图高潮gif福利片| 午夜免费激情av| 欧美大码av| 国产一区二区三区视频了| 国产黄a三级三级三级人| 国内精品久久久久久久电影| 色婷婷久久久亚洲欧美| 欧美激情极品国产一区二区三区| 桃色一区二区三区在线观看| 麻豆久久精品国产亚洲av| 色综合站精品国产| 亚洲av电影不卡..在线观看| 亚洲人成网站在线播放欧美日韩| 久久中文看片网| 欧美日韩福利视频一区二区| 日本熟妇午夜| 亚洲久久久国产精品| 99久久国产精品久久久| 欧美日本视频| www日本黄色视频网| av超薄肉色丝袜交足视频| 亚洲天堂国产精品一区在线| 久久久久亚洲av毛片大全| 久久人妻福利社区极品人妻图片| 午夜福利成人在线免费观看| 国产一级毛片七仙女欲春2 | 人妻久久中文字幕网| 国产精品98久久久久久宅男小说| 亚洲国产欧洲综合997久久, | 亚洲精品在线观看二区| 啪啪无遮挡十八禁网站| 窝窝影院91人妻| 欧美亚洲日本最大视频资源| 国产蜜桃级精品一区二区三区| 欧美日韩福利视频一区二区| 国产午夜福利久久久久久| 国产高清有码在线观看视频 | 亚洲专区字幕在线| 啦啦啦 在线观看视频| 精品久久久久久成人av| 欧美日本亚洲视频在线播放| 国产麻豆成人av免费视频| 欧美精品亚洲一区二区| 搡老妇女老女人老熟妇| 中文资源天堂在线| 天堂动漫精品| 级片在线观看| 美女午夜性视频免费| bbb黄色大片| 变态另类成人亚洲欧美熟女| cao死你这个sao货| 99国产精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 国产片内射在线| 欧美av亚洲av综合av国产av| 非洲黑人性xxxx精品又粗又长| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 在线国产一区二区在线| 女生性感内裤真人,穿戴方法视频| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 欧美激情高清一区二区三区| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 国产av又大| 天天一区二区日本电影三级| 亚洲精品粉嫩美女一区| 国产久久久一区二区三区| 亚洲精品国产区一区二| 制服丝袜大香蕉在线| 国产97色在线日韩免费| 变态另类成人亚洲欧美熟女| 久久久精品欧美日韩精品| 亚洲男人的天堂狠狠| 午夜视频精品福利| 成人亚洲精品一区在线观看| 丰满人妻熟妇乱又伦精品不卡| 12—13女人毛片做爰片一| 俄罗斯特黄特色一大片| 男女那种视频在线观看| 国产男靠女视频免费网站| 久久精品成人免费网站| 美国免费a级毛片| 欧美激情 高清一区二区三区| 在线看三级毛片| 老熟妇乱子伦视频在线观看| 1024香蕉在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲av片天天在线观看| 亚洲欧美日韩无卡精品| av视频在线观看入口| 啪啪无遮挡十八禁网站| 亚洲国产毛片av蜜桃av| 黄色视频,在线免费观看| 欧美人与性动交α欧美精品济南到| 国产精品美女特级片免费视频播放器 | 亚洲五月天丁香| 特大巨黑吊av在线直播 | 午夜福利在线在线| 国产精品香港三级国产av潘金莲| av有码第一页| а√天堂www在线а√下载| 在线看三级毛片| 久久热在线av| 手机成人av网站| 搡老岳熟女国产| 久久久国产精品麻豆| 欧美成人一区二区免费高清观看 | 国内精品久久久久精免费| 亚洲精品一区av在线观看| 欧美日韩乱码在线| 日日爽夜夜爽网站| 久热这里只有精品99| 99热6这里只有精品| 亚洲三区欧美一区| 亚洲精品中文字幕一二三四区| 看免费av毛片| 国产蜜桃级精品一区二区三区| 国内精品久久久久久久电影| 国产亚洲精品av在线| 亚洲激情在线av| 成年免费大片在线观看| 悠悠久久av| 18禁美女被吸乳视频| 高潮久久久久久久久久久不卡| 亚洲精品在线观看二区| 1024手机看黄色片| 一级黄色大片毛片| 国产一区二区三区在线臀色熟女| 国产激情偷乱视频一区二区| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 色在线成人网| 欧美成狂野欧美在线观看| 欧美乱色亚洲激情| 欧美性猛交黑人性爽| 午夜久久久久精精品| 在线视频色国产色| 久久午夜亚洲精品久久| 不卡av一区二区三区| 亚洲片人在线观看| 特大巨黑吊av在线直播 | 日本 欧美在线| 国产亚洲精品av在线| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 亚洲 欧美一区二区三区| 日日摸夜夜添夜夜添小说| 999精品在线视频| 免费电影在线观看免费观看| 真人一进一出gif抽搐免费| 国产精品98久久久久久宅男小说| 女性生殖器流出的白浆| 精品无人区乱码1区二区| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 国产欧美日韩一区二区三| 国产精品 国内视频| 黑人操中国人逼视频| 国产精品一区二区三区四区久久 | svipshipincom国产片| 麻豆久久精品国产亚洲av| av超薄肉色丝袜交足视频| av电影中文网址| 黄色视频不卡| 国内精品久久久久精免费| 搡老岳熟女国产| 很黄的视频免费| 亚洲精品国产区一区二| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频| 免费高清视频大片| 老司机福利观看| 大型黄色视频在线免费观看| 午夜福利在线观看吧| 日本撒尿小便嘘嘘汇集6| 亚洲国产高清在线一区二区三 | 午夜成年电影在线免费观看| 亚洲全国av大片| 黄片大片在线免费观看| 国产精品国产高清国产av| 国产精品久久久人人做人人爽| 免费看日本二区| 一级黄色大片毛片| 免费搜索国产男女视频| 国产亚洲精品一区二区www| 亚洲精品美女久久av网站| 午夜成年电影在线免费观看| av中文乱码字幕在线| 男人操女人黄网站| 国产亚洲精品久久久久久毛片| 欧美成人性av电影在线观看| 黑人操中国人逼视频| 久久久久久久久久黄片| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 免费观看人在逋| 免费看日本二区| 亚洲全国av大片| 国产色视频综合| 777久久人妻少妇嫩草av网站| 国产黄色小视频在线观看| 成人永久免费在线观看视频| 99久久无色码亚洲精品果冻| 在线观看免费视频日本深夜| 嫁个100分男人电影在线观看| x7x7x7水蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一码二码三码区别大吗| 亚洲欧美精品综合一区二区三区| 成人手机av| 少妇粗大呻吟视频| 成年免费大片在线观看| 日韩精品免费视频一区二区三区| xxx96com| 欧美性猛交黑人性爽| netflix在线观看网站| 99国产综合亚洲精品| 亚洲最大成人中文| 精品久久久久久久末码| 日本成人三级电影网站| 99精品久久久久人妻精品| 国产精品影院久久| 久久精品国产99精品国产亚洲性色| 国产精品久久久av美女十八| 国产高清videossex| 亚洲久久久国产精品| 国产1区2区3区精品| 久久中文字幕人妻熟女| 国产黄片美女视频| 亚洲色图 男人天堂 中文字幕| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 精品久久久久久久人妻蜜臀av| 亚洲精品久久成人aⅴ小说| 中文字幕人妻熟女乱码| 亚洲专区国产一区二区| 亚洲av电影在线进入| 久久国产精品人妻蜜桃| 在线观看一区二区三区| 成年免费大片在线观看| 欧美+亚洲+日韩+国产| 最近在线观看免费完整版| 又紧又爽又黄一区二区| 天堂√8在线中文| 日本一本二区三区精品| 国语自产精品视频在线第100页| 欧美绝顶高潮抽搐喷水| 国产欧美日韩精品亚洲av| 50天的宝宝边吃奶边哭怎么回事| 欧美黑人欧美精品刺激| 日本一本二区三区精品| 久久婷婷人人爽人人干人人爱| 中出人妻视频一区二区| 欧美日本亚洲视频在线播放| 欧美色欧美亚洲另类二区| 99国产极品粉嫩在线观看| 可以在线观看的亚洲视频| 国产乱人伦免费视频| 欧美日韩乱码在线| 中文字幕精品免费在线观看视频| 91成年电影在线观看| 成人一区二区视频在线观看| 国产精品免费一区二区三区在线| 首页视频小说图片口味搜索| 伊人久久大香线蕉亚洲五| 欧美日韩一级在线毛片| 国产亚洲av高清不卡| 日本 欧美在线| 亚洲午夜精品一区,二区,三区| 美女高潮到喷水免费观看| 精品第一国产精品| 国产精品久久久人人做人人爽| 午夜福利一区二区在线看| 欧美色欧美亚洲另类二区| 999久久久国产精品视频| 国产精品久久久久久人妻精品电影| 99久久99久久久精品蜜桃| 变态另类丝袜制服| 精品久久久久久久久久免费视频| 国产精品香港三级国产av潘金莲| 久久性视频一级片| 手机成人av网站| 满18在线观看网站| 91九色精品人成在线观看| 亚洲狠狠婷婷综合久久图片| 成人18禁在线播放| 国产精品久久久人人做人人爽| 99久久综合精品五月天人人| 变态另类成人亚洲欧美熟女| 搡老妇女老女人老熟妇| 在线十欧美十亚洲十日本专区| 成在线人永久免费视频| 男男h啪啪无遮挡| 色综合站精品国产| 国产成人av激情在线播放| 午夜a级毛片| 啦啦啦韩国在线观看视频| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 国产欧美日韩一区二区精品| 搡老岳熟女国产| 精品一区二区三区av网在线观看| 制服丝袜大香蕉在线| 亚洲人成网站在线播放欧美日韩| 麻豆成人午夜福利视频| 亚洲精品在线美女| 美女 人体艺术 gogo| av在线播放免费不卡| 免费在线观看日本一区| av视频在线观看入口| 国产成人欧美| 成熟少妇高潮喷水视频| 亚洲精品粉嫩美女一区| 老司机在亚洲福利影院| 午夜精品久久久久久毛片777| 97人妻精品一区二区三区麻豆 | 久久香蕉精品热| 97超级碰碰碰精品色视频在线观看| 他把我摸到了高潮在线观看| 亚洲国产日韩欧美精品在线观看 | 日本五十路高清| 伦理电影免费视频| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| 两人在一起打扑克的视频| 18禁国产床啪视频网站| 国产亚洲精品一区二区www| a在线观看视频网站| 男女午夜视频在线观看| 制服人妻中文乱码| 亚洲狠狠婷婷综合久久图片| 熟女电影av网| aaaaa片日本免费| 国产91精品成人一区二区三区| 欧美乱码精品一区二区三区| www.www免费av| 午夜影院日韩av| 人人澡人人妻人| 亚洲精品中文字幕一二三四区| 老司机靠b影院| 一卡2卡三卡四卡精品乱码亚洲| 亚洲久久久国产精品| 少妇熟女aⅴ在线视频| 国产亚洲av嫩草精品影院| 国产男靠女视频免费网站| xxx96com| 久久精品国产亚洲av高清一级| 成人av一区二区三区在线看| 国产成人啪精品午夜网站| 91成年电影在线观看| 久久香蕉精品热| 亚洲成国产人片在线观看| 亚洲av美国av| 亚洲国产精品合色在线| 人人妻人人澡人人看| 午夜免费观看网址| 国产成+人综合+亚洲专区| 真人一进一出gif抽搐免费| 久久狼人影院| 国产1区2区3区精品| 免费看十八禁软件| 色在线成人网| 免费在线观看日本一区| 日本 欧美在线| 人人妻人人澡欧美一区二区| 国产成人av激情在线播放| 99国产精品一区二区蜜桃av| 精品电影一区二区在线| 日本在线视频免费播放| 叶爱在线成人免费视频播放| 老汉色∧v一级毛片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲性夜色夜夜综合| 亚洲精品国产区一区二| 欧美中文综合在线视频| 免费高清在线观看日韩| 国产99白浆流出| 午夜久久久在线观看| 亚洲人成77777在线视频| 中文亚洲av片在线观看爽| www国产在线视频色| 久久香蕉国产精品| 在线观看免费视频日本深夜| 亚洲第一欧美日韩一区二区三区| 久久久久久国产a免费观看| 欧美成人午夜精品| 精品午夜福利视频在线观看一区| 国产黄a三级三级三级人| 日本成人三级电影网站| 国产v大片淫在线免费观看| 亚洲精品粉嫩美女一区| 国产精品亚洲一级av第二区| 91麻豆av在线| 精品久久久久久久毛片微露脸| 亚洲三区欧美一区| 午夜福利成人在线免费观看| 天天添夜夜摸| 悠悠久久av| 国产人伦9x9x在线观看| 日日干狠狠操夜夜爽| 国产精品久久久久久亚洲av鲁大| 久久久久久久久免费视频了| 男女视频在线观看网站免费 | 免费女性裸体啪啪无遮挡网站| 亚洲精品国产区一区二| 欧美国产日韩亚洲一区| 久久中文字幕人妻熟女|