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

    兩種被動微波遙感混合像元分解方法比較

    2013-09-25 14:27:54顧玲嘉任瑞治
    吉林大學學報(地球科學版) 2013年6期
    關鍵詞:亮溫陸地被動

    顧玲嘉,趙 凱,任瑞治,孫 健,

    1.吉林大學電子科學與工程學院,長春 130012

    2.中國科學院東北地理與農(nóng)業(yè)生態(tài)研究所,長春 130012

    0 引言

    被動微波遙感技術(shù)在全球遙感觀測中有著廣泛的應用,它能夠穿透云層獲得地表一定深度的信息[1]。由于星載微波輻射計無法裝載大型天線,被動微波數(shù)據(jù)的空間分辨率都較低[2]。通常使用的星載被動微波輻射計,如SSM/I、AMSR-E和FY3BMWRI的空間分辨率都是幾十千米[3-6]。對于觀測地區(qū)的數(shù)據(jù),低空間分辨率導致獲得的微波數(shù)據(jù)是其像元內(nèi)多類地物的混合結(jié)果。

    經(jīng)前期調(diào)研[7],針對星載被動微波數(shù)據(jù)空間分辨率低導致的混合像元分問題,國內(nèi)外研究者基于圖像處理技術(shù)解決的方案主要有2類:一種是圖像超分辨率增強方法,另一種是混合像元分解方法。2類方法的實現(xiàn)都可以提高被動微波地物反演參數(shù),但是處理的方法和最終獲得的數(shù)據(jù)有所不同。

    被動微波數(shù)據(jù)超分辨率增強方法可以直接提高亮溫的分辨率。其主要的理論依據(jù)是衛(wèi)星的觀測數(shù)據(jù)腳印之間有很多部分的重疊,利用這些重疊觀測可以獲得更高分辨率的亮溫數(shù)據(jù),這種方法也是目前國內(nèi)外研究者多數(shù)采用的方法。與之相比,混合像元分解方法對原有被動微波數(shù)據(jù)處理后,獲得的是原混合像元中的各組分亮溫,其組分分類根據(jù)具體的應用有所不同,組分的獲得需要借助于高于被動微波空間分辨率的地物分類數(shù)據(jù)。將組分亮溫數(shù)據(jù)代入到地物參數(shù)反演中,可以提高直接采用混合數(shù)據(jù)進行反演的結(jié)果。例如若直接使用水陸邊界處的被動微波混合數(shù)據(jù)進行海水反演,由于空間分辨率低,該混合像元數(shù)據(jù)是陸地和海洋的綜合結(jié)果,相對于海洋表面,陸地表面呈現(xiàn)高發(fā)射率,導致反演精度減低。如果通過微波混合像元分解,將該混合像元內(nèi)的海水組分亮溫提取出來并進行海水參數(shù)的反演,其反演數(shù)據(jù)結(jié)果的精度肯定是要高于前者。

    基于地物分類的被動微波混合像元分解方法取決于微波像元內(nèi)地物的分類。目前國內(nèi)外采用的被動微波混合方法主要針對沿海地區(qū)水陸分界比較明顯的區(qū)域進行分解,像元內(nèi)的地物分類主要分為海洋和陸地,代表方法是Bellerby等[8]結(jié)合海岸線數(shù)據(jù)利用相鄰像元信息的混合像元分解方法。Bennartz[9]采用周圍完全是陸地類型的亮溫數(shù)據(jù)平均值來代替陸地組分亮溫[10-13],然后代入建立的陸地水體混合像元分解模型,求得天線足印中海洋的組分亮溫。MaaB等[14]研究并提出了一種類似于Bennartz的混合像元分解方法,兩者的主要差異在于應用領域。MaaB應用該方法來提高海冰參數(shù),而Bennartz利用該方法實現(xiàn)大氣參數(shù)的反演[15]。在Bellerby提出的微波混合像元分解算法的基礎上,筆者[16]曾提出基于水陸分類信息的被動微波混合像元分解方法,并對2010年中國吉林省永吉縣洪澇災害進行了有效評估。

    上述被動微波混合像元分解方法都是利用已知地物分類的先驗知識,建立被動微波混合像元分解模型,通過求解方程組獲得被動微波混合像元中水體和陸地的組分亮溫,其有效性的驗證都是根據(jù)具體的研究地區(qū)展開,在各自研究的地區(qū)內(nèi)可以保證其組分亮溫求解的精度。其中,Bennartz提出的基于陸地組分亮溫的混合像元分解模型相對簡單,其陸地組分亮溫是由周圍完全是陸地類型的亮溫數(shù)據(jù)平均值計算得到;因此,這種模型的適用范圍較小。筆者提出基于地物分類信息的被動微波混合像元分解方法,并將其與Bellerby提出的分解方法進行對比研究??紤]到實際測試場地中較難獲得不同地物的被動微波組分亮溫數(shù)據(jù),筆者通過模擬具有不同組分亮溫分布的觀測地區(qū),設計具有不同組分亮溫分布的觀測地區(qū),對2種分解方法進行評價,進一步分析2種分解方法的適用范圍和影響其求解精度的主要因素,為建立具有更高求解精度和適用范圍的被動微波混合像元分解方法提供理論基礎。

    1 基于地物分類信息的被動微波混合像元分解方法

    Bellerby假設相鄰9個像元區(qū)域中陸地和海洋組分的亮溫值相同,這樣會使計算結(jié)果產(chǎn)生一定誤差[16];筆者在 Bellerby算法[17]的基礎上,針對其算法的局限性進行改進。已知被動微波混合像元亮溫值TB,混合像元中水體分布比例矩陣PW和陸地分布比例矩陣PL,考慮到鄰近區(qū)域內(nèi)水體被動微波組分亮溫變化不大的特征,可對m×n窗口尺寸內(nèi)的被動微波混合數(shù)據(jù)進行分解,通過構(gòu)建方程組求解水體組分亮溫TW和陸地組分亮溫TL:

    令L(i,j)為水體陸地分類結(jié)果。當點(i,j)是陸地上的點時,L(i,j)為1;當點(i,j)是水體的點時,L(i,j)為0。水體陸地分類結(jié)果來自于觀測地區(qū)土地利用分類數(shù)據(jù)。則P(i,j)為天線增益函數(shù)ˉG與水體陸地分類結(jié)果的卷積:

    此外,還需設計一個目標函數(shù)作為判斷可能解優(yōu)劣的標準。選取F作為目標函數(shù),ξ為事先給定的閾值,其取值取決于對解的精度要求。目標函數(shù)F定義為

    該混合像元分解模型中,m×n方程組中有m×n+1個未知量,即TL中的m×n個量和TW。這種方程數(shù)小于未知量數(shù)的方程組屬于欠定方程組,有無窮多組解,求得最優(yōu)解是該算法需要解決的問題。

    筆者提出在整個解空間中搜索最優(yōu)解的計算方法,對不同類型的組分亮溫值設定其對應的取值區(qū)間。具體采用k-means聚類算法對所有的被動微波亮溫數(shù)據(jù)進行統(tǒng)計分類,聚類后的中心值代表不同組分的綜合初值,即CL和CW。根據(jù)陸地被動微波混合像元的統(tǒng)計結(jié)果如均值、最大值、最小值、方差,以及混合像元中不同組分的比例,計算陸地亮溫的變化閾值RL,確定TL(i,j)初值選取范圍是[CLRL,CL+RL];同理,計算水體亮溫的變化閾值RW,確定TW初值選取范圍是[CW-RW,CW+RW]。由于水體亮溫應小于陸地亮溫,需滿足CW+RW<CL-RL,對初值選取范圍進行重新設定。

    根據(jù)目標函數(shù)調(diào)整搜索窗口尺寸m×n,采用帶約束條件的非負最小二乘法迭代運算求得最優(yōu)解。由于初值的選取范圍考慮到不同地物的亮溫變化空間,且存在多個條件限制方程組的求解范圍,因此采用筆者提出的方法可以得到較好的最優(yōu)解。

    2 仿真實驗與分析

    仿真實驗中,首先設計星載被動微波天線增益函數(shù),將其與具有不同模擬組分亮溫值分布的觀測地區(qū)進行卷積,獲得該地區(qū)的微波混合像元亮溫;然后,分別采用Bellerby和筆者提出的分解方法對混合像元進行分解;最后,將其各自分解后得到的組分亮溫與原始組分亮溫進行比較,進而驗證2種被動微波混合像元分解方法的有效性。

    2.1 微波天線增益函數(shù)

    通常,星載微波天線增益函數(shù)近似滿足二維高斯分布。實際的天線增益方向圖是一個包括旁瓣的復雜三維結(jié)構(gòu),Drusch等[18]研究表明天線函數(shù)可以近似為高斯分布函數(shù)。因此,天線增益函數(shù)的公式定義為

    其中:g0和g1是描述高斯分布的常數(shù);x和y是當前位置到測量中心點的水平和垂直距離;a和b表示-3dB橢圓的半軸,可以根據(jù)被動微波數(shù)據(jù)中不同頻率對應的足印來確定;橢圓半徑r定義為

    -3dB點定義為從中心點到天線增益降低到最大值一半的位置,Drusch等[18]將式(6)代入式(5),推導得到

    本研究使用-3dB點作為橢圓半徑,則微波天線增益函數(shù)為

    2.2 試驗樣區(qū)設計

    假設觀測區(qū)域為75km×75km的范圍,組分亮溫網(wǎng)格的空間分辨率為5km×5km,則共包含有15×15個網(wǎng)格。令被動微波數(shù)據(jù)的采樣分辨率為25km×25km,考慮到相對于其他地表類型,水體類型的被動微波遙感數(shù)據(jù)亮溫值變化的動態(tài)范圍較小,根據(jù)被動微波遙感數(shù)據(jù)AMSR-E和MWRI的主要下墊面輻射亮溫統(tǒng)計特性,分別采用如下方案進行樣地組分亮溫設計。

    方案一:假設觀測地區(qū)內(nèi)所有陸地組分亮溫相同,為260K,所有水體組分亮溫相同,為120K。水體組分亮溫在區(qū)域呈十字分布,令水體比例逐漸增加,構(gòu)成被動微波組分亮溫。具體組分亮溫分布比例如圖1a、b所示,其中藍色的圖像位置代表水體,紅色的圖像位置代表陸地,都呈十字分布。由于定義的天線增益函數(shù)中橢圓的長短半軸均為25km,所以用每5×5網(wǎng)格的組分亮溫與天線增益函數(shù)作卷積計算,計算后得到的模擬被動微波混合亮溫數(shù)據(jù)為3像元×3像元,具體值如圖1c、d所示。由圖1c、d可見,水體組分亮溫比例的增加使混合像元亮溫數(shù)據(jù)的十字中心位置更接近水體組分的亮溫值,即表現(xiàn)為藍色圖像。

    分別采用Bellerby和筆者提出的方法對微波混合像元進行分解,筆者提出的模型采用的窗口尺寸與Bellerby方法一致,為3像元×3像元,其各組分亮溫的初值選取由地物分類矩陣中各地物的比例動態(tài)估計。將2種模型分解后得到的組分亮溫與原始組分亮溫進行比較。實驗結(jié)果表明:在方案一條件下,2種方法求得的組分亮溫與原組分亮溫的誤差為0K,即滿足該條件下,2種方法都可以完全獲得混合像元的組分亮溫。

    方案二:假設觀測地區(qū)內(nèi)所有陸地組分亮溫相同,為260K;所有水體組分亮溫相同,為120K。水體組分亮溫位置在區(qū)域呈隨機分布。令水體比例逐漸增加,構(gòu)成被動微波組分亮溫。其具體分布如圖2a、b所示。計算后得到的被動微波混合亮溫數(shù)據(jù)如圖2c、d所示。從圖2c、d可見,水體組分亮溫的隨機分布使得混合像元亮溫數(shù)據(jù)的分布也較隨機,與方案一中混合像元亮溫呈十字分布有所不同。

    實驗結(jié)果表明:在方案二條件下,采用Bellerby方法和筆者提出的方法求得的組分亮溫與原組分亮溫的誤差也為0K。

    方案一和方案二的結(jié)論證明:對于地物類型僅為水體和陸地2類的觀測地區(qū),若模型求解區(qū)域?qū)挠^測地區(qū)中陸地亮溫組分相同、水體亮溫組分相同,只要已知地物分類信息,無論水體的位置和比例如何變化,采用這2種混合像元分解模型均可以從被動混合像元中準確地獲得對應的水體和陸地組分亮溫。

    方案三:假設當觀測地區(qū)沒有水體時,陸地組分亮溫每列相同,且從第1列到第15列由246K到260K遞增變化,變化精度為1K(圖3a)。當觀測地區(qū)含有水體時,根據(jù)相鄰微波亮溫數(shù)據(jù)中水體組分亮溫變化不大的特點,仍假設觀測區(qū)域內(nèi)水體組分亮溫相同,為120K,且水體組分亮溫位置與方案二中的分布位置和比例變化相同,即呈隨機分布并逐漸增加,構(gòu)成如圖3b、c所示的被動微波組分亮溫。將亮溫組分數(shù)據(jù)和天線增益方程進行卷積計算,得到其對應的被動微波混合亮溫如圖3d、e、f所示。圖3d中,由于水體組分亮溫的比例為0,因此構(gòu)成的混合像元亮溫圖仍是陸地分布。隨著水體組分亮溫比例增加,其對應的混合像元亮溫圖3e、f由黃色向藍色變化,即水體比例的增加使得構(gòu)成的混合像元亮溫數(shù)據(jù)值降低。

    觀測地區(qū)陸地和水體組分亮溫的數(shù)據(jù)統(tǒng)計結(jié)果和采用2種方法計算得到的組分亮溫與實際組分亮溫差值的統(tǒng)計結(jié)果見表1和表2。

    由表1和表2可見,當水體比例由0逐漸增加時:對于陸地組分亮溫,Bellerby方法得到的組分亮溫絕對誤差的平均值約為3.1K,而筆者提出方法獲得的組分亮溫絕對誤差的平均值約為0.1K;對于水體組分亮溫,Bellerby方法得到的組分亮溫絕對誤差的平均值約為5.9K,筆者提出方法獲得的組分亮溫據(jù)對誤差的平均值約為0.2K。當觀測地區(qū)陸地組分亮溫變化時,由于Bellerby采用的混合像元分解方法假設3×3混合像元內(nèi)所有陸地組分亮溫相同,所有水體組分亮溫也相同,因此采用最小二乘法求得的陸地組分亮溫精度與3×3混合像元區(qū)域即75km×75km范圍內(nèi)的組分亮溫變化有關。與之相比,筆者采用的混合像元分解模型假設3×3混合像元中每個混合像元內(nèi)的陸地組分亮溫相同,即每25km×25km范圍內(nèi)的陸地組分亮溫一致,同時假設所有混合像元中的水體組分亮溫相同,因此采用該方法求得的陸地組分亮溫的誤差主要來自于25km×25km范圍(即5×5個組分亮溫網(wǎng)格)內(nèi)陸地組分亮溫的變化。因此,筆者提出的方法能得到更加準確的分解結(jié)果。

    圖1 方案一組分亮溫分布和混合像元亮溫Fig.1 The distribution of component brightness temperature and mixed data in the first project

    圖2 方案二中組分亮溫分布和混合像元亮溫Fig.2 The distribution of component brightness temperature and mixed data in the second project

    圖3 方案三中組分亮溫分布和混合像元亮溫Fig.3 The distribution of component brightness temperature and mixed data in the third project

    表1 方案三陸地組分亮溫數(shù)據(jù)統(tǒng)計Table1 Statistics results of land component brightness temperature in the third project

    表2 方案三水體組分亮溫數(shù)據(jù)統(tǒng)計Table2 Statistics results of water component brightness temperature in the third project

    3 結(jié)論與建議

    1)由于陸地、水體的亮度溫度差異很大,且水體亮度溫度比較穩(wěn)定,因此對于地物類型僅為水體和陸地2類的觀測地區(qū),若模型求解區(qū)域?qū)挠^測地區(qū)中陸地、水體亮溫組分相同,只要已知地物分類信息,無論水體的位置和比例如何變化,采用Bellerby和筆者提出的混合像元分解模型均可以從被動混合像元中準確地獲得實際的水體組分和陸地組分亮溫。

    2)當觀測地區(qū)陸地組分亮溫變化時,筆者提出的方法進行混合像元分解后獲得的組分亮溫的精度要高于Bellerby提出的方法。筆者提出的微波混合像元分解方法中需要對于欠定方程組進行求解,其中各組分亮溫的初值變化范圍和條件對于求解的精度尤為重要。欠定方程組有無窮多組解,實驗結(jié)果表明,如果定義好求解條件,可以從多個解中找到最優(yōu)解。

    3)與水體組分亮溫相比,陸地上不同類型地表的亮度溫度差異較大,而且觀測區(qū)域內(nèi)陸地類型的微波輻射亮度溫度也不完全一致,這是導致目前提出的微波混合像元分解方法求解精度降低的主要原因。通過對實驗結(jié)果進行分析,如果對微波混合像元中陸地類型進行再次分類,即將陸地根據(jù)地表分類數(shù)據(jù)劃分為不同地表類型,進一步求得微波混合像元中多類地物的比例,建立基于多類地物分類的被動微波混合像元分解模型,通過合理定義每種類型地物組分亮溫范圍和設定求解約束條件,將有可能進一步提高被動微波混合像元求解的精度,這也是未來需要研究的工作。

    長春凈月潭遙感實驗站工作人員對本文仿真實驗方案設計提供了幫助和支持,謹致謝忱!

    (References):

    [1]Thomas J J.Passive Microwave Remote Sensing for Land Applications[J].Earth and Environmental Science,2008(1):9-18.

    [2]尹紅剛,張德海.使用BG反演算法細分重建微波輻射計圖像[J].遙感技術(shù)與應用,2006,21(2):120-124.Yin Honggang,Zhang Dehai.Subdivided Reconstruction of Microwave Radiometer Images with BG Inversion Method[J].Remote Sensing Technology and Application,2006,21(2):120-124.

    [3]Derksen C,Walker A,Ledrew E,et al.Combining SMMR and SSM/I Data for Time Series Analysis of Central North American Snow Water Equivalent[J].Journal of Hydrometeorology,2003,4(2):304-316.

    [4]Paloscia S,Macelloni G,Pampaloni P,et al.Global Scale Monitoring of Soil and Vegetation by Using AMSR-E Multi-Temporal Data [C ]//International Geoscience and Remote Sensing Symposium(IGARSS).Denver:IEEE,2006:1609-1612.

    [5]楊虎,施建成.FY-3微波成像儀地表參數(shù)反演研究[J].遙感技術(shù)與應用,2005,20(1):194-200.Yang Hu,Shi Jiancheng.On the Estimation of Land Surface Parameters by Using FY-3AMicrowave Radiometer Imager(MWRI)[J].Remote Sensing Technology and Application,2005,20(1):194-200.

    [6]Gu Lingjia,Zhao Kai,Zhang Shuwen,et al.Comparative Analysis of Microwave Brightness Temperature Data in Northeast China Using AMSR-E and MWRI Products[J].Chinese Geographical Science,2011,21(1):84-93.

    [7]顧玲嘉,趙凱,孫健,等.被動微波遙感數(shù)據(jù)超分辨率增強與混合像元分解研究綜述[J].遙感技術(shù)與應用,2012,27(1):25-31.Gu Lingjia,Zhao Kai,Sun Jian,et al.Review on the Super Resolution Enhancement and Unmixing Method for Passive Microwave Remote Sensing Data[J].Remote Sensing Technology and Application.2012,27(1):25-31.

    [8]Bellerby T,Taberner M,Wilmshurst A.Passive Microwave Retrieval of Ocean Surface Windspeeds for British Coastal Waters[J].IEEE Transactions on Geoscience and Remote Sensing,1999,37(5):2578-2584.

    [9]Bennartz R.On the Use of SSM/I Measurements in Coastal Regions[J].Journal of Atmospheric and Oceanic Technology,1999,16(4):417-431.

    [10]Bennartz R,Paape K,Gutman G G.Global Land Surface Characterization Using Combined Passive Microwave and Visible/Near-Infrared Satellite Data[J].Eos,Transactions,American Geophysical U-nion,1997,78(1):59-60.

    [11]Grody N C.Surface Identification Using Satellite Microwave Radiometers[J].IEEE Transactions on Geoscience and Remote Sensing,1988,26(6):850-859.

    [12]Jones A S,Vonder Haar T H.Retrieval of Microwave Surface Emittance Over Land Using Coincident Microwave and Infrared Satellite Measurements[J].Journal of Geophysical Research Atmospheres,1997,102(D12):13609-13626.

    [13]Prigent C,Rossow B W,Matthews E.Microwave Land Surface Emissivities Estimated from SSM/I Observations[J].Journal of Geophysical Research,1997,102(D18):21867-21890.

    [14]MaaB N,Kaleschke L.Improving Passive Microwave Sea Ice Concentration Algorithms for Coastal Areas:Applications to the Baltic Sea[J].Tellus Series A:Dynamic Meteorology and Oceanography,2010,62(4):393-410.

    [15]Wessel P,Smith W H F.A Global,Self-Consistent,Hierarchical,High-Resolution Shoreline Database[J].J Geophys Res:B Solid Earth,1996,101(B4):8741-8743.

    [16]Gu LingJia,Zhao Kai,Zhang Shuang,et al.An AMSR-E Data Unmixing Method for Monitoring Flood and Waterlogging Disaster[J].Chinese Geographical Science,2011,21(6):666-675.

    [17]Bellerby T,Taberner M,Wilmshurst A.Retrieval of Land and Sea Brightness Temperatures from Mixed Coastal Pixels in Passive Microwave Data[J].IEEE Transactions on Geoscience and Remote Sensing,1998,36(6):1844-1851.

    [18]Drusch M,Wood E F,Lindau R.The Impact of the SSM/I Antenna Gain Function on Land Surface Parameter Retrieval[J].Geophysical Research Letter,1999,26(23):3481-3484.

    猜你喜歡
    亮溫陸地被動
    誰在推著陸地跑
    霰譜分布特征對強對流云高頻微波亮溫影響的模擬研究
    新聞語篇中被動化的認知話語分析
    基于南太平洋的AMSR2 L1R亮溫數(shù)據(jù)質(zhì)量評估
    海洋通報(2020年2期)2020-09-04 09:22:34
    陸地開來“宙斯盾”
    主動句都能轉(zhuǎn)換成被動句嗎
    第五課 拒絕被動
    趣味(語文)(2019年5期)2019-09-02 01:52:44
    爬爬爬,以水中沖向陸地
    2014年2月12日新疆于田MS7.3地震熱紅外亮溫異常分析
    中國地震(2015年1期)2015-11-08 11:11:32
    2011年盈江5.8級地震熱輻射亮溫異常分析
    地震研究(2014年3期)2014-02-27 09:30:49
    99热网站在线观看| 汤姆久久久久久久影院中文字幕| 亚洲少妇的诱惑av| 亚洲成人一二三区av| 麻豆精品久久久久久蜜桃| 国产 精品1| 国产 精品1| 悠悠久久av| 男女无遮挡免费网站观看| videosex国产| 一区在线观看完整版| 男女边吃奶边做爰视频| 一级毛片我不卡| 亚洲专区中文字幕在线 | 亚洲国产av新网站| 黑人欧美特级aaaaaa片| 国产精品免费视频内射| 成年动漫av网址| 亚洲欧美精品自产自拍| 亚洲免费av在线视频| 国产免费视频播放在线视频| 老汉色∧v一级毛片| 成年动漫av网址| 亚洲精品美女久久av网站| 亚洲欧美精品自产自拍| 电影成人av| 精品久久久久久电影网| 亚洲国产av新网站| 一边亲一边摸免费视频| 国产高清国产精品国产三级| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 国产乱人偷精品视频| av在线app专区| 精品少妇黑人巨大在线播放| 国产成人精品在线电影| 国产女主播在线喷水免费视频网站| av福利片在线| 日本一区二区免费在线视频| 亚洲欧美成人综合另类久久久| 纯流量卡能插随身wifi吗| 一边摸一边做爽爽视频免费| 麻豆乱淫一区二区| 青春草国产在线视频| 国产精品久久久久久精品古装| 97人妻天天添夜夜摸| 免费黄色在线免费观看| 国产精品国产三级国产专区5o| 国产无遮挡羞羞视频在线观看| 国产精品熟女久久久久浪| 9热在线视频观看99| 久久精品久久久久久久性| 国产精品久久久久久精品古装| 国产成人系列免费观看| 国产毛片在线视频| 亚洲欧美一区二区三区黑人| 中文精品一卡2卡3卡4更新| 免费在线观看视频国产中文字幕亚洲 | 久久性视频一级片| 亚洲情色 制服丝袜| 国产精品一二三区在线看| 亚洲成人国产一区在线观看 | 丰满乱子伦码专区| 亚洲五月色婷婷综合| 另类亚洲欧美激情| 在线观看免费日韩欧美大片| 91老司机精品| 中文字幕制服av| 在线天堂中文资源库| 欧美日韩一级在线毛片| 日日摸夜夜添夜夜爱| 亚洲国产成人一精品久久久| 国产一区二区激情短视频 | 久久99精品国语久久久| 青青草视频在线视频观看| 热re99久久精品国产66热6| 国产成人一区二区在线| 欧美日韩视频精品一区| 国产av国产精品国产| 女的被弄到高潮叫床怎么办| 波野结衣二区三区在线| 最新在线观看一区二区三区 | 国产乱来视频区| 久久亚洲国产成人精品v| 免费高清在线观看视频在线观看| xxxhd国产人妻xxx| 18禁观看日本| 人体艺术视频欧美日本| 国产一区亚洲一区在线观看| 久久久久精品性色| av视频免费观看在线观看| 国产黄色免费在线视频| 国产高清国产精品国产三级| 丝袜在线中文字幕| 久久久久久久久久久免费av| 成人国产av品久久久| 90打野战视频偷拍视频| 国产成人免费观看mmmm| 婷婷色综合www| 一个人免费看片子| 高清黄色对白视频在线免费看| 高清av免费在线| 亚洲精品成人av观看孕妇| a级片在线免费高清观看视频| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 国产免费福利视频在线观看| 国产亚洲欧美精品永久| 日韩一区二区三区影片| 啦啦啦 在线观看视频| 少妇人妻 视频| 丰满少妇做爰视频| 一区福利在线观看| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| 亚洲少妇的诱惑av| 母亲3免费完整高清在线观看| 国产视频首页在线观看| 久久久久久免费高清国产稀缺| 精品人妻熟女毛片av久久网站| 亚洲欧洲国产日韩| 老司机深夜福利视频在线观看 | 无遮挡黄片免费观看| 制服丝袜香蕉在线| 亚洲精品国产一区二区精华液| 97精品久久久久久久久久精品| 国产不卡av网站在线观看| 中文欧美无线码| 亚洲精品av麻豆狂野| 嫩草影视91久久| 精品国产露脸久久av麻豆| 高清av免费在线| 韩国高清视频一区二区三区| 久久久久久久久久久免费av| 又粗又硬又长又爽又黄的视频| 国产精品免费视频内射| 捣出白浆h1v1| 国产97色在线日韩免费| 国产精品99久久99久久久不卡 | av国产久精品久网站免费入址| 波多野结衣av一区二区av| 精品福利永久在线观看| 久久久久久久精品精品| 91成人精品电影| 中文天堂在线官网| 久久99热这里只频精品6学生| 免费黄网站久久成人精品| 精品国产露脸久久av麻豆| 午夜久久久在线观看| 丰满饥渴人妻一区二区三| 伦理电影大哥的女人| 久久久久国产一级毛片高清牌| 亚洲欧美一区二区三区国产| 黄色一级大片看看| 色综合欧美亚洲国产小说| 建设人人有责人人尽责人人享有的| 黄色视频不卡| 成人国产av品久久久| 一级片免费观看大全| 国产精品人妻久久久影院| 嫩草影院入口| 卡戴珊不雅视频在线播放| 97精品久久久久久久久久精品| 日日啪夜夜爽| 亚洲国产精品999| 亚洲人成77777在线视频| 高清黄色对白视频在线免费看| 久久精品久久久久久噜噜老黄| 黑人猛操日本美女一级片| 9191精品国产免费久久| 飞空精品影院首页| 美女午夜性视频免费| 欧美人与善性xxx| 亚洲七黄色美女视频| 亚洲国产av新网站| 丰满迷人的少妇在线观看| 免费在线观看视频国产中文字幕亚洲 | 在线看a的网站| 最近最新中文字幕免费大全7| 亚洲欧洲精品一区二区精品久久久 | 精品国产露脸久久av麻豆| 午夜av观看不卡| 黄片播放在线免费| 亚洲专区中文字幕在线 | 日韩欧美精品免费久久| 极品少妇高潮喷水抽搐| 成人亚洲精品一区在线观看| 最黄视频免费看| 99久久综合免费| 自线自在国产av| videosex国产| 99热网站在线观看| 色婷婷久久久亚洲欧美| 午夜福利一区二区在线看| 母亲3免费完整高清在线观看| 男女国产视频网站| 韩国av在线不卡| 亚洲成av片中文字幕在线观看| 18在线观看网站| 欧美日韩视频高清一区二区三区二| 亚洲精品中文字幕在线视频| 亚洲男人天堂网一区| 国产黄色免费在线视频| 欧美另类一区| 欧美国产精品va在线观看不卡| 天天躁夜夜躁狠狠躁躁| 免费黄网站久久成人精品| 国产毛片在线视频| 麻豆av在线久日| 国产精品 欧美亚洲| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 午夜福利在线免费观看网站| 精品国产一区二区久久| 欧美日韩视频高清一区二区三区二| 另类精品久久| 9热在线视频观看99| 亚洲久久久国产精品| 国产成人91sexporn| 黄色视频在线播放观看不卡| 人人妻人人爽人人添夜夜欢视频| 最新的欧美精品一区二区| 精品国产国语对白av| 久久久久久人妻| 老司机影院毛片| 国产亚洲最大av| 日日啪夜夜爽| 欧美日韩成人在线一区二区| 色婷婷久久久亚洲欧美| 国产免费现黄频在线看| 久久久久精品国产欧美久久久 | 久久综合国产亚洲精品| 国产色婷婷99| 成人国产麻豆网| 久热爱精品视频在线9| 人人妻人人澡人人爽人人夜夜| 国产人伦9x9x在线观看| 巨乳人妻的诱惑在线观看| 成人黄色视频免费在线看| 成年人午夜在线观看视频| 国产成人精品福利久久| 成人影院久久| 免费少妇av软件| 美女主播在线视频| 九色亚洲精品在线播放| 香蕉丝袜av| 国产av国产精品国产| av视频免费观看在线观看| 性少妇av在线| 亚洲国产欧美日韩在线播放| 久久免费观看电影| 欧美人与性动交α欧美软件| 久久久久精品人妻al黑| www.熟女人妻精品国产| tube8黄色片| 亚洲天堂av无毛| 9热在线视频观看99| 欧美最新免费一区二区三区| 免费观看性生交大片5| 80岁老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆 | 亚洲成色77777| 一区福利在线观看| 波多野结衣av一区二区av| 久久99一区二区三区| 男女高潮啪啪啪动态图| 亚洲国产中文字幕在线视频| 男男h啪啪无遮挡| 男女午夜视频在线观看| 日韩av在线免费看完整版不卡| 成人国产av品久久久| 国产成人免费无遮挡视频| 在线 av 中文字幕| 丁香六月天网| 一级片免费观看大全| 中文字幕最新亚洲高清| 国产爽快片一区二区三区| 热99久久久久精品小说推荐| 99热国产这里只有精品6| av国产精品久久久久影院| 日韩制服丝袜自拍偷拍| 满18在线观看网站| 久久久精品94久久精品| 国产亚洲av高清不卡| 欧美老熟妇乱子伦牲交| 日韩大码丰满熟妇| 久久精品久久精品一区二区三区| 人人妻,人人澡人人爽秒播 | 亚洲国产精品999| 日本wwww免费看| 9热在线视频观看99| 精品人妻熟女毛片av久久网站| 免费不卡黄色视频| 精品国产超薄肉色丝袜足j| 最近中文字幕2019免费版| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡 | 国产精品女同一区二区软件| 亚洲美女搞黄在线观看| 亚洲七黄色美女视频| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网 | www.自偷自拍.com| 精品免费久久久久久久清纯 | 日韩精品有码人妻一区| 高清欧美精品videossex| 精品一品国产午夜福利视频| 免费人妻精品一区二区三区视频| 欧美 亚洲 国产 日韩一| 性高湖久久久久久久久免费观看| 天美传媒精品一区二区| a级片在线免费高清观看视频| www.精华液| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区久久| 免费黄网站久久成人精品| 成人黄色视频免费在线看| 麻豆精品久久久久久蜜桃| 亚洲国产欧美一区二区综合| 国产精品久久久久久精品电影小说| 最近中文字幕2019免费版| 在线观看免费高清a一片| 人人妻人人爽人人添夜夜欢视频| 久久ye,这里只有精品| 美女中出高潮动态图| 亚洲精品乱久久久久久| 男女午夜视频在线观看| 亚洲免费av在线视频| 亚洲精品一区蜜桃| 色婷婷av一区二区三区视频| 高清在线视频一区二区三区| 国产成人91sexporn| 久久人人97超碰香蕉20202| 日本色播在线视频| 美女视频免费永久观看网站| 国产熟女午夜一区二区三区| 国产精品久久久久久久久免| 国产熟女午夜一区二区三区| 在线亚洲精品国产二区图片欧美| 免费在线观看黄色视频的| 妹子高潮喷水视频| 免费观看a级毛片全部| 久久久久久免费高清国产稀缺| 欧美亚洲 丝袜 人妻 在线| 老司机影院成人| 精品免费久久久久久久清纯 | 精品国产一区二区三区四区第35| 亚洲综合色网址| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕色久视频| 波多野结衣av一区二区av| 日日撸夜夜添| tube8黄色片| 亚洲成国产人片在线观看| 精品亚洲乱码少妇综合久久| 一级毛片 在线播放| 韩国精品一区二区三区| 久久久精品国产亚洲av高清涩受| 国产乱人偷精品视频| 亚洲国产av新网站| 精品亚洲成国产av| 在线观看三级黄色| 亚洲国产看品久久| 精品少妇内射三级| 亚洲av电影在线进入| 丝袜美足系列| 精品国产乱码久久久久久男人| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人av在线免费| av卡一久久| 90打野战视频偷拍视频| 美国免费a级毛片| 精品国产露脸久久av麻豆| 国语对白做爰xxxⅹ性视频网站| 成人亚洲欧美一区二区av| 久久精品国产综合久久久| 亚洲av电影在线进入| netflix在线观看网站| 波野结衣二区三区在线| 中文字幕精品免费在线观看视频| 国产成人精品久久久久久| 免费日韩欧美在线观看| 国产精品久久久久久人妻精品电影 | 又黄又粗又硬又大视频| 国产xxxxx性猛交| 国产一区二区三区av在线| av电影中文网址| 精品人妻熟女毛片av久久网站| 七月丁香在线播放| 免费日韩欧美在线观看| 国产毛片在线视频| 国产精品免费大片| 午夜免费男女啪啪视频观看| 我的亚洲天堂| 90打野战视频偷拍视频| 国产精品久久久久成人av| 男女免费视频国产| 亚洲精品第二区| 成人手机av| 欧美成人精品欧美一级黄| 日本一区二区免费在线视频| 九九爱精品视频在线观看| 精品一区二区免费观看| 一级爰片在线观看| 男女高潮啪啪啪动态图| 99re6热这里在线精品视频| 日韩成人av中文字幕在线观看| 在线精品无人区一区二区三| av卡一久久| 人体艺术视频欧美日本| 美女国产高潮福利片在线看| 亚洲一区二区三区欧美精品| 亚洲精品在线美女| 少妇人妻 视频| videos熟女内射| 久久久久视频综合| av视频免费观看在线观看| 在线观看国产h片| 国产精品二区激情视频| 午夜福利视频精品| 精品免费久久久久久久清纯 | 99久久人妻综合| 中文字幕人妻熟女乱码| 黑人巨大精品欧美一区二区蜜桃| 日本一区二区免费在线视频| 日日啪夜夜爽| 自拍欧美九色日韩亚洲蝌蚪91| 日韩电影二区| 搡老乐熟女国产| 91精品三级在线观看| 日本爱情动作片www.在线观看| 欧美日本中文国产一区发布| 亚洲欧美一区二区三区久久| 国产福利在线免费观看视频| 亚洲精品成人av观看孕妇| 91国产中文字幕| 中文字幕亚洲精品专区| 菩萨蛮人人尽说江南好唐韦庄| 国产野战对白在线观看| 午夜激情久久久久久久| 午夜福利,免费看| 久久国产亚洲av麻豆专区| 欧美日韩一级在线毛片| 伊人久久国产一区二区| 免费久久久久久久精品成人欧美视频| 黄片小视频在线播放| 日日摸夜夜添夜夜爱| 丝袜美腿诱惑在线| 最近的中文字幕免费完整| 日本av免费视频播放| 自线自在国产av| 中文字幕色久视频| 亚洲精品一二三| 国产av国产精品国产| 国产av精品麻豆| 在线 av 中文字幕| 黄片无遮挡物在线观看| 99久久人妻综合| 自拍欧美九色日韩亚洲蝌蚪91| av视频免费观看在线观看| 日韩成人av中文字幕在线观看| 亚洲成色77777| 热99国产精品久久久久久7| 老鸭窝网址在线观看| 妹子高潮喷水视频| 七月丁香在线播放| 另类亚洲欧美激情| 91精品国产国语对白视频| 欧美日韩一区二区视频在线观看视频在线| 最近2019中文字幕mv第一页| 欧美日韩一级在线毛片| 免费少妇av软件| 久久影院123| 精品酒店卫生间| 精品久久久精品久久久| 1024视频免费在线观看| 国产爽快片一区二区三区| 大码成人一级视频| 一级爰片在线观看| 国产一区二区 视频在线| 十八禁网站网址无遮挡| 夜夜骑夜夜射夜夜干| 午夜日本视频在线| 午夜精品国产一区二区电影| 国产一区亚洲一区在线观看| 欧美黄色片欧美黄色片| 亚洲国产欧美网| 午夜福利网站1000一区二区三区| 激情视频va一区二区三区| 成人国产麻豆网| 97在线人人人人妻| av片东京热男人的天堂| 国产精品一二三区在线看| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 午夜激情av网站| 亚洲国产日韩一区二区| 色视频在线一区二区三区| 考比视频在线观看| 亚洲国产欧美在线一区| 熟妇人妻不卡中文字幕| 欧美日韩成人在线一区二区| 麻豆乱淫一区二区| 麻豆av在线久日| 亚洲精品第二区| 18在线观看网站| 亚洲欧美日韩另类电影网站| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 国产免费福利视频在线观看| 老司机影院成人| 欧美精品一区二区大全| 水蜜桃什么品种好| 久久天躁狠狠躁夜夜2o2o | 超碰成人久久| 51午夜福利影视在线观看| 久久精品亚洲熟妇少妇任你| 亚洲国产精品一区二区三区在线| 免费日韩欧美在线观看| 国产精品久久久久久人妻精品电影 | 狠狠精品人妻久久久久久综合| 国产男人的电影天堂91| 国产在线视频一区二区| 一边摸一边抽搐一进一出视频| 欧美人与性动交α欧美精品济南到| 天天影视国产精品| av国产精品久久久久影院| 亚洲av日韩在线播放| 婷婷色av中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 91精品三级在线观看| 日本欧美国产在线视频| 美女主播在线视频| 午夜福利乱码中文字幕| 热re99久久国产66热| 日韩av不卡免费在线播放| 久久久久久人妻| 欧美另类一区| 欧美日韩综合久久久久久| 亚洲五月色婷婷综合| 99热网站在线观看| 亚洲男人天堂网一区| 欧美人与性动交α欧美软件| 哪个播放器可以免费观看大片| 在线观看一区二区三区激情| 色视频在线一区二区三区| 久久精品国产亚洲av高清一级| 婷婷色av中文字幕| 国产亚洲欧美精品永久| 女性生殖器流出的白浆| 侵犯人妻中文字幕一二三四区| 日本vs欧美在线观看视频| 宅男免费午夜| 2021少妇久久久久久久久久久| 亚洲熟女精品中文字幕| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜爱| 国产精品 国内视频| 亚洲精品一区蜜桃| 黑人欧美特级aaaaaa片| 国产黄频视频在线观看| 高清在线视频一区二区三区| 9色porny在线观看| 国产av码专区亚洲av| 久久久久久免费高清国产稀缺| xxx大片免费视频| 国产成人精品福利久久| 侵犯人妻中文字幕一二三四区| 亚洲av男天堂| 在线精品无人区一区二区三| av女优亚洲男人天堂| 久久精品久久久久久噜噜老黄| 午夜福利一区二区在线看| 一区二区三区乱码不卡18| 香蕉国产在线看| 成年动漫av网址| 少妇猛男粗大的猛烈进出视频| 欧美精品一区二区大全| 亚洲精品一区蜜桃| 精品一区二区三区av网在线观看 | 好男人视频免费观看在线| 超色免费av| 欧美在线黄色| 啦啦啦视频在线资源免费观看| 国产精品 国内视频| 精品亚洲成a人片在线观看| 黄色怎么调成土黄色| 纵有疾风起免费观看全集完整版| 丝瓜视频免费看黄片| 欧美久久黑人一区二区| 永久免费av网站大全| 亚洲伊人色综图| 亚洲av电影在线观看一区二区三区| 男女之事视频高清在线观看 | 国产色婷婷99| 午夜日韩欧美国产| 91aial.com中文字幕在线观看| 精品视频人人做人人爽| 亚洲欧美一区二区三区久久| 男人操女人黄网站| 两个人免费观看高清视频| 日韩免费高清中文字幕av| 一区在线观看完整版| 亚洲人成网站在线观看播放| 最近手机中文字幕大全| 欧美日韩福利视频一区二区| 亚洲人成网站在线观看播放| 高清欧美精品videossex| 日韩精品有码人妻一区| 国产伦理片在线播放av一区| 亚洲国产精品成人久久小说| 99热网站在线观看| 国产欧美亚洲国产| 麻豆精品久久久久久蜜桃|