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

    基于Landsat 8數(shù)據(jù)的蕪湖市地表溫度反演

    2018-05-14 08:59:49湯志戴照福
    安徽農(nóng)業(yè)科學 2018年20期
    關鍵詞:透射率水汽波段

    湯志 戴照福

    摘要 以空間分辨率100 m的Landsat 8 TIRS熱紅外傳感器第10波段數(shù)據(jù)為基礎,采用對大氣透射率估算方程修正后的單窗算法,反演了蕪湖市地表溫度,并使用MODIS地溫產(chǎn)品數(shù)據(jù)為標準對其進行精度評價。評價結果表明,修正后的算法所得到的反演結果具有較高的數(shù)據(jù)精度,能夠更加準確細致地表示出地表溫度空間分布的細節(jié)信息。

    關鍵詞 單窗算法;地表溫度;Landsat 8;蕪湖市

    中圖分類號 P237 文獻標識碼 A 文章編號 0517-6611(2018)20-0047-04

    Abstract Based on the tenth band data of the Landsat 8 TIRS thermal infrared senor with a spatial resolution of 100 m, we used the monowindow algorithm modified by the atmospheric transmittance estimation equation to retrieve the land surface temperature of Wuhu City and used the MODIS land surface temperature product data as the standard to evaluate its accuracy. The evaluation results showed as followed: the retrieved results obtained by the modified algorithm had a higher data accuracy and the detailed information of land surface temperature spatial distribution could be expressed more accurately.

    Key words Monowindow algorithm;Land surface temperature;Landsat 8;Wuhu

    地表溫度在地理學、地球物理學和地球生物化學等領域研究中是一個十分重要的參量,在地-氣相互作用如潛熱傳輸和顯熱傳輸中是個很大的影響因素,地表溫度的量化對氣溫[1]、植被[2]、土壤含水量[3]和生態(tài)系統(tǒng)[4]等研究有重要意義。

    要了解區(qū)域內(nèi)部廣泛連續(xù)的地表溫度時空分布狀況,目前較多是通過熱紅外遙感手段來獲取。國內(nèi)外學者在這方面已經(jīng)做了大量研究,Price[5]使用NOAA/AVHRR遙感數(shù)據(jù)首次將海面溫度遙感分裂窗算法應用到陸地表面溫度反演中。Jiménez-Muoz等[6]在2003年提出了一種只需要大氣透射率和地表反射率的單通道算法。覃志豪[7]針對TM6衛(wèi)星數(shù)據(jù)提出了一種只需要1個熱紅外波段即可進行地表溫度反演的算法,稱為單窗算法,該算法所需參數(shù)少,且只需用1個熱紅外波段而被人廣泛采用。

    基于MODIS數(shù)據(jù)和分裂窗算法反演地表溫度的方法已經(jīng)十分成熟,雖然其反演結果精度高,但該數(shù)據(jù)空間分辨率較低,為1 km,無法細致地描述地表溫度的空間分布,將逐漸不能滿足相關領域的研究和發(fā)展。Landsat 8 TIRS傳感器中的2個熱紅外波段(Band10和Band11)數(shù)據(jù)空間分辨率均為100 m,但由于Band11波段的參數(shù)目前還不穩(wěn)定,可能會造成較大的反演誤差,故筆者選取Band10波段,采用單窗算法進行地表溫度的反演,并對結果進行了精度驗證和分析比較,以期為高空間分辨率熱紅外波段數(shù)據(jù)的地表溫度反演提

    供參考。

    1 原理與方法

    1.1 單窗算法

    覃志豪[7]根據(jù)地表熱量輻射傳輸方程,推導出一種針對TM熱紅外波段數(shù)據(jù)的地表溫度反演算法(Monowindow Algorithm),該算法只需地表輻射率、大氣透射率和大氣平均作用溫度這3個參數(shù)即可計算地表溫度,算法不僅較容易實現(xiàn),而且保持了較高的精度,公式為:

    式中,TS為地表溫度;a、b為系數(shù);τ為大氣透射率;ε為地表輻射率;Tb為遙感器接收到的輻射亮度值所轉(zhuǎn)化而來的亮度溫度;Ta為大氣平均作用溫度。

    1.2 a和b的取值

    Rozenstein等[8]利用LOWTRAN大氣輻射傳輸軟件模擬了不同溫度下TIRS傳感器第10和第11波段的反演回歸系數(shù)a和b(表1)。根據(jù)表1,先估算出研究區(qū)氣溫,從而可查詢到該溫度在表中相應的回歸系數(shù)。溫度區(qū)間越短,相關系數(shù)平方(r2)越大,估計值標準誤差(SEE)越小,回歸系數(shù)的精度越高,故表1中不同溫度區(qū)間都含有估計出的溫度值時,應選擇較短溫度區(qū)間所對應的回歸系數(shù)。該研究選取溫度為10~40 ℃的a、b值,即a10=-62.806,b10=0.434。

    1.3 大氣透射率的計算

    因CO2和N2等氣體的影響在同一景遙感影像中可認為是不變的,所以大氣透射率值主要是由大氣水汽含量值所決定的[8]。大氣的水汽含量不是一個定值,為了獲得更加精確的大氣透射率,可以通過MODTRAN模擬出TIRS第10波段的大氣透射率和水汽含量之間的變化關系,建立擬合方程,為提高相關系數(shù)平方(R2),降低擬合結果的均方根誤差(RMSE),將大氣水汽含量分成了2個范圍分別進行擬合[9],擬合結果見表2。

    在表2的擬合方程中,γ表示大氣透射率,w表示大氣水汽含量。在沒有實測氣象數(shù)據(jù)的情況下,在該景影像的可見光波段進行FLAASH大氣校正后由模型生成的記事本文件中,可以查詢到大氣校正模型對該景影像的大氣水汽含量的估測值,而在有實測氣象數(shù)據(jù)的情況下,可首先利用相對濕度和氣溫與絕對水汽壓之間的關系式[式(4)]求出研究區(qū)的絕對水汽壓[10],再由絕對水汽壓求得該區(qū)域的大氣水汽含量[式(5)][11]。

    通過查詢國家氣象科學數(shù)據(jù)共享服務平臺的中國地面氣候資料日值數(shù)據(jù)得到研究區(qū)當日RH=0.53,T0=15 ℃,計算得出絕對水氣壓e=1.48 kPa,并按照公式(5)計算出大氣水汽含量w=3.11 g/cm2,故采用大氣水汽含量在3.0~6.0 g/cm2的擬合方程計算大氣透射率,計算得出τ=0.63。

    1.4 地表輻射率的確定

    地表輻射率主要由地表結構和遙感器波段區(qū)間所決定,在單窗算法下,由于只使用了1個熱紅外波段,所以在估計地表輻射率時只需考慮地表結構的影響。地球表面類型多樣,結構復雜,因受傳感器分辨率的限制,目前Landsat 8衛(wèi)星數(shù)據(jù)還不能將地表結構完全展現(xiàn)出來,Qin等[12]提出的混合像元法將地表劃分成水體、建筑、裸土和植被4種類型,認為對于TM數(shù)據(jù)的地表溫度遙感反演來說,上述劃分已經(jīng)可以滿足地溫反演對地表輻射率的精度要求。

    2 蕪湖市地表溫度反演實例

    2.1 研究區(qū)概況

    蕪湖市地處長江三角洲西南部,北鄰長江,南倚皖南山區(qū),區(qū)域經(jīng)緯度為117°57′40″~118°43′23″E,30°38′35″~31°28′26″N,市區(qū)總面積為6 026 km2,是國家區(qū)域中心城市、長江三角洲城市群城市以及南京都市圈城市。全市北低南高,地貌類型多樣,丘陵山地主要集中于西南部,東北部地勢平坦,河流縱橫交錯,土地覆蓋類型主要為植被、建筑等(圖1)。

    2.2 數(shù)據(jù)來源及預處理

    2.2.1 遙感數(shù)據(jù)。

    該研究選取的數(shù)據(jù)為2013年4月7日10:39過境的Landsat 8 OLI/TIRS C1 Level-1數(shù)據(jù),數(shù)據(jù)來源于USGS網(wǎng)站。其中,Landsat 8 OLI陸地成像儀包含多光譜波段,空間分辨率為30 m,用于地表發(fā)射率的計算,TIRS包含2個熱紅外波段(第10和第11波段),空間分辨率為100 m,該研究選用第10波段反演地表溫度。數(shù)據(jù)預處理包括影像拼接、裁剪、輻射定標和大氣校正等。因覆蓋研究區(qū)需要2幅影像,所以首先需要進行影像拼接,拼接時使用直方圖匹配來消除2幅圖像的色差,再裁剪出研究區(qū)的影像,對OLI數(shù)據(jù)進行輻射定標和大氣校正,對TIRS數(shù)據(jù)按照公式(12)將原始DN值轉(zhuǎn)換為熱輻射亮度值用于亮度溫度的計算。

    2.2.2 驗證數(shù)據(jù)。

    由于缺少地面實測數(shù)據(jù),該研究使用同一天MODIS衛(wèi)星MODLT1D傳感器的地表溫度產(chǎn)品數(shù)據(jù)來進行地表溫度反演結果的精度驗證,數(shù)據(jù)來源于地理空間數(shù)據(jù)云網(wǎng)站。該數(shù)據(jù)的DN值是按照16位圖像進行顯示的,范圍在0~655 35,故需要將其DN值乘以0.02轉(zhuǎn)換成通常所用的絕對溫度,數(shù)據(jù)如圖2所示。

    2.3 反演結果

    將所求得的地表輻射率、大氣透射率和大氣平均作用溫度等各參數(shù)帶入公式(1)中,反演出研究區(qū)的地表溫度(圖3)。圖3-a為反演過程中生成的研究區(qū)域亮度溫度圖像,圖3-b為研究區(qū)單窗算法反演地表溫度結果圖像。

    2.3.1 精度驗證。

    通過對研究區(qū)隨機布點,生成15個用于地表溫度反演結果精度驗證的樣本點,統(tǒng)計各個點處的反演溫度(Ts)和MODIS數(shù)據(jù)中的地溫(Tm),使用這兩者之間差值的絕對值(| Ts-Tm |)來表示Landsat 8 TIRSBand10波段地表溫度反演結果與MODIS地溫的最終誤差。

    由表3可知,研究區(qū)地表溫度反演結果與MODIS地溫平均誤差為0.81 ℃,說明地表溫度反演結果的總體精度較高,大部分樣點處的地表溫度誤差都在1 ℃以下,少數(shù)樣點地溫在1 ℃以上。經(jīng)對比發(fā)現(xiàn),地表溫度較高(25 ℃左右)的樣本點誤差相對較大(大于1 ℃);地表溫度較低(20 ℃左右)的樣本點誤差相對較小(小于1 ℃)。宋挺等[15]在其所做的Landsat 8 數(shù)據(jù)地表溫度反演算法對比研究中認為,原始單窗算法反演Landsat 8數(shù)據(jù)得出的地溫平均值與MODIS數(shù)據(jù)的地溫平均值相差了2.11 ℃。由此可見,基于Landsat 8數(shù)據(jù)使用單窗算法反演地表溫度過程中,根據(jù)研究區(qū)的大氣水汽含量值對原始算法中的大氣透射率參數(shù)進行修正,可以顯著提高反演結果的數(shù)據(jù)精度。

    2.3.2 反演結果。

    通過對數(shù)據(jù)進行統(tǒng)計并結合圖3可以得出,研究區(qū)亮度溫度最低為19.81 ℃,最高為30.47 ℃,平均為26.58 ℃,地表溫度最低為18.74 ℃,最高為31.52 ℃,平均為23.81 ℃。亮度溫度的空間分布與地表溫度的空間分布基本一致,亮度溫度值整體高于地表溫度值,這是由于高空中的遙感衛(wèi)星傳感器所接收到的輻射不僅包括來自地表的直接輻射,還包括大氣輻射的地面反射以及大氣多次散射等各種輻射,這就使得到達遙感器處的輻射值高于地表輻射值,在圖像中的表現(xiàn)即為亮度溫度高于地表溫度。在圖3-b中,地表溫度較高的區(qū)域主要分布在道路、裸地以及人口聚集的城區(qū),地表溫度較低的區(qū)域主要分布在河流、湖泊等水體以及植被茂密的山體處,反演結果的空間分布符合研究區(qū)實際情況。

    3 結論與討論

    基于Landsat 8 TIRS Band10熱紅外波段數(shù)據(jù),使用對大氣透射率修正后的單窗算法反演了蕪湖市地表溫度,并利用MODIS地溫產(chǎn)品數(shù)據(jù)進行精度評價,得出平均誤差為0.81 ℃,說明該算法反演地表溫度的結果精度較好,能夠反映研究區(qū)地表溫度分布的實際情況。單窗算法所得到的研究區(qū)地表溫度分布與MODIS地溫數(shù)據(jù)中的地表溫度分布基本保持一致,且Landsat 8因其熱紅外波段數(shù)據(jù)空間分辨率較高,故能夠比MODIS地溫數(shù)據(jù)更為精確細致地表示出地表溫度分布狀況,能夠為城市土壤環(huán)境管理、作物種植、農(nóng)業(yè)規(guī)劃、森林火災監(jiān)測等一系列與地表溫度相關的領域在現(xiàn)實中的應用提供更為精確的指導建議和輔助決策。

    該研究的不足之處在于,由于單窗算法只使用1個熱紅外波段數(shù)據(jù)進行地溫反演,傳感器所成圖像原始亮度值直接影響到反演的結果,所以算法對熱紅外傳感器的性能和大氣輻射傳輸過程的模擬精度要求很高,也就是說該算法中大氣水汽含量和大氣平均作用溫度的準確估算對保證反演結果的精度至關重要,如何更為精確地量化該算法所需要的參數(shù)以及如何利用Landsat 8 2個熱紅外波段數(shù)據(jù)做劈窗算法反演地表溫度,均是今后需要思考的方向。

    參考文獻

    [1] 姜會飛,廖樹華,葉爾克江,等.地面溫度與氣溫關系的統(tǒng)計分析[J].中國農(nóng)業(yè)氣象,2004,25(3):1-4.

    [2] 許娟,楊武年,任娟.基于Landsat數(shù)據(jù)地表特征參數(shù)與地表溫度關系[J].測繪與空間地理信息,2015(1):56-59.

    [3] 劉斐,陳軍,慕軍營,等.土壤溫度檢測及其與含水率關系研究[J].干旱地區(qū)農(nóng)業(yè)研究,2013,31(3):95-99.

    [4] 李琪,薛紅喜,王云龍,等.土壤溫度和水分對克氏針茅草原生態(tài)系統(tǒng)碳通量的影響初探[J].農(nóng)業(yè)環(huán)境科學學報,2011,30(3):605-610.

    [5] PRICE J C.Land surface temperature measurements from the split window channels of the NOAA 7 Advanced Very High Resolution Radiometer[J].Journal of geophysical research,1984,89(D5):7231-7237.

    [6] JIMNEZMUN~OZ J C,SOBRINO J A.A generalized singlechannel method for retrieving land surface temperature from remote sensing data[J].Journal of geophysical research,2003,108(D22):4688-4697.

    [7] 覃志豪.用陸地衛(wèi)星TM6數(shù)據(jù)演算地表溫度的單窗算法[J].地理學報,2001,56(4):456-466.

    [8] ROZENSTEIN O,QIN Z H,DERIMIAN Y,et al.Derivationof land surface temperature for Landsat8 TIRS using a split window algorithm[J].Sensors,2014,14(4):5768-5780.

    [9] 蔣大林,匡鴻海,曹曉峰,等.基于Landsat8的地表溫度反演算法研究:以滇池流域為例[J].遙感技術與應用,2015,30(3):448-454.

    [10] 喬建民.不同土地覆被類型對城市地表溫度的影響:以山東省龍口市為例[J].寧波農(nóng)業(yè)科技,2012(3):21-25.

    [11] 楊景梅,邱金桓.用地面濕度參量計算我國整層大氣可降水量及有效水汽含量方法的研究[J].大氣科學,2002,26(1):9-22.

    [12] 覃志豪,李文娟,徐斌,等.陸地衛(wèi)星TM6波段范圍內(nèi)地表比輻射率的估計[J].國土資源遙感,2004,16(3):28-32.

    [13] USGS.Using the USGS Landsat 8 Product [EB/OL].(2013-12-04)[2018-02-20].http://landsat.usgs.gov/Landsat8_ Using _ Product.php.

    [14] 覃志豪,LI W J,ZHANG M H,等.單窗算法的大氣參數(shù)估計方法[J].國土資源遙感,2003,15(2):37-43.

    [15] 宋挺,段崢,劉軍志,等.Landsat8數(shù)據(jù)地表溫度反演算法對比[J].遙感學報,2015,19(3):451-464.

    猜你喜歡
    透射率水汽波段
    春日暖陽
    青藏高原上空平流層水汽的時空演變特征
    純色太陽鏡鏡片耐日光輻照性能試驗研究
    1979~2011年間平流層溫度及平流層水汽的演變趨勢
    M87的多波段輻射過程及其能譜擬合
    深圳“5·11”特大暴雨過程的水汽輸送特征分析
    日常維護對L 波段雷達的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    石墨烯光柵太赫茲透射特性的研究
    玻璃長波透射率對普通雙層中空玻璃熱特性的影響
    晶格常數(shù)對左手材料的能量透射率及負折射率的影響
    計算物理(2014年1期)2014-03-11 17:00:48
    国产精品久久久人人做人人爽| 欧美性长视频在线观看| 精品熟女少妇八av免费久了| 满18在线观看网站| 极品少妇高潮喷水抽搐| 五月天丁香电影| 青青草视频在线视频观看| 国产精品亚洲av一区麻豆| 视频在线观看一区二区三区| 免费在线观看完整版高清| 精品亚洲成a人片在线观看| a级毛片黄视频| 丰满人妻熟妇乱又伦精品不卡| 免费看十八禁软件| 欧美黑人欧美精品刺激| 久久久国产欧美日韩av| 国产日韩欧美在线精品| 久久久久网色| 亚洲五月婷婷丁香| 亚洲av第一区精品v没综合| 他把我摸到了高潮在线观看 | 一二三四在线观看免费中文在| 9色porny在线观看| www.精华液| 国产亚洲午夜精品一区二区久久| 亚洲国产av新网站| 国产精品秋霞免费鲁丝片| 99香蕉大伊视频| 这个男人来自地球电影免费观看| 高清黄色对白视频在线免费看| 色94色欧美一区二区| 日本精品一区二区三区蜜桃| 精品人妻在线不人妻| 亚洲国产毛片av蜜桃av| 亚洲精品国产区一区二| 久久久久久久国产电影| 一级黄色大片毛片| 在线观看舔阴道视频| 久久 成人 亚洲| 久久热在线av| 国产成人啪精品午夜网站| 亚洲成a人片在线一区二区| 亚洲色图 男人天堂 中文字幕| 久久99一区二区三区| 一个人免费在线观看的高清视频| 日韩欧美一区视频在线观看| 不卡一级毛片| 曰老女人黄片| 亚洲美女黄片视频| 一本一本久久a久久精品综合妖精| videosex国产| av网站免费在线观看视频| 亚洲一码二码三码区别大吗| tocl精华| 一本—道久久a久久精品蜜桃钙片| av网站在线播放免费| 女警被强在线播放| 香蕉国产在线看| 91国产中文字幕| 热99re8久久精品国产| 精品一区二区三区四区五区乱码| 中文字幕色久视频| 怎么达到女性高潮| 精品国产乱子伦一区二区三区| 国产精品电影一区二区三区 | 青草久久国产| 日韩有码中文字幕| 午夜福利在线观看吧| 亚洲性夜色夜夜综合| 国产精品九九99| 人妻 亚洲 视频| 精品人妻熟女毛片av久久网站| 成年人免费黄色播放视频| 女人久久www免费人成看片| 国产精品.久久久| 亚洲精品中文字幕一二三四区 | 精品卡一卡二卡四卡免费| 欧美日韩国产mv在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲成人国产一区在线观看| 午夜福利视频在线观看免费| 18禁观看日本| 国产日韩欧美视频二区| 丰满少妇做爰视频| 青青草视频在线视频观看| 午夜福利视频精品| 久久人妻av系列| 成人国语在线视频| 视频区图区小说| 成人手机av| 国产精品 国内视频| 最近最新中文字幕大全免费视频| 亚洲性夜色夜夜综合| 亚洲精品美女久久av网站| 在线观看免费视频网站a站| 国产野战对白在线观看| 狠狠精品人妻久久久久久综合| 亚洲五月色婷婷综合| www日本在线高清视频| 国产av国产精品国产| 国产在线精品亚洲第一网站| 欧美成人免费av一区二区三区 | 亚洲视频免费观看视频| 在线观看一区二区三区激情| 嫩草影视91久久| 麻豆乱淫一区二区| 最近最新中文字幕大全电影3 | 人人妻人人添人人爽欧美一区卜| 91麻豆av在线| av片东京热男人的天堂| 国产熟女午夜一区二区三区| 高清毛片免费观看视频网站 | bbb黄色大片| 国产av又大| 最新的欧美精品一区二区| 精品福利永久在线观看| 天堂中文最新版在线下载| 人人妻人人澡人人爽人人夜夜| 亚洲av电影在线进入| 啦啦啦 在线观看视频| 中文字幕人妻丝袜一区二区| 亚洲,欧美精品.| 日本黄色日本黄色录像| 国产精品久久久av美女十八| 亚洲全国av大片| 欧美亚洲日本最大视频资源| 国产亚洲一区二区精品| 国产精品1区2区在线观看. | 人人妻人人澡人人看| bbb黄色大片| 午夜福利在线免费观看网站| 午夜福利一区二区在线看| 国产主播在线观看一区二区| 国产淫语在线视频| 国产淫语在线视频| 国产真人三级小视频在线观看| 美女午夜性视频免费| 久久精品熟女亚洲av麻豆精品| 免费黄频网站在线观看国产| 精品亚洲成国产av| 欧美av亚洲av综合av国产av| 一个人免费看片子| videos熟女内射| 国产亚洲精品第一综合不卡| 日本五十路高清| videos熟女内射| 丁香六月天网| 日韩 欧美 亚洲 中文字幕| 欧美精品高潮呻吟av久久| 少妇精品久久久久久久| 精品亚洲成国产av| 少妇精品久久久久久久| 精品视频人人做人人爽| 99re在线观看精品视频| 一个人免费看片子| 丁香六月天网| 天天躁狠狠躁夜夜躁狠狠躁| 无遮挡黄片免费观看| 国产又爽黄色视频| 久久狼人影院| 巨乳人妻的诱惑在线观看| 久久久国产欧美日韩av| 国产日韩欧美在线精品| 久久狼人影院| av视频免费观看在线观看| 精品人妻在线不人妻| 国产不卡av网站在线观看| 亚洲av美国av| 国产在线视频一区二区| 岛国毛片在线播放| 露出奶头的视频| 黄色视频,在线免费观看| 精品久久久久久久毛片微露脸| 性高湖久久久久久久久免费观看| 亚洲精品粉嫩美女一区| 日韩三级视频一区二区三区| 19禁男女啪啪无遮挡网站| 女人久久www免费人成看片| 精品国产一区二区三区久久久樱花| 欧美日韩亚洲综合一区二区三区_| 香蕉丝袜av| 国产成人精品无人区| 国产精品 国内视频| 两性夫妻黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产熟女午夜一区二区三区| 欧美日韩一级在线毛片| 三级毛片av免费| 狂野欧美激情性xxxx| 久久久久精品国产欧美久久久| 青青草视频在线视频观看| 亚洲欧美日韩高清在线视频 | 巨乳人妻的诱惑在线观看| 久久精品国产亚洲av香蕉五月 | 国产精品麻豆人妻色哟哟久久| videos熟女内射| 亚洲精品久久成人aⅴ小说| 国产不卡av网站在线观看| 欧美激情极品国产一区二区三区| 国产精品av久久久久免费| 男女免费视频国产| 久久天躁狠狠躁夜夜2o2o| 欧美亚洲 丝袜 人妻 在线| 久9热在线精品视频| 老司机午夜十八禁免费视频| 精品久久蜜臀av无| av一本久久久久| 精品国内亚洲2022精品成人 | 欧美乱妇无乱码| 久久这里只有精品19| 一级片'在线观看视频| av天堂久久9| 91精品国产国语对白视频| 欧美 亚洲 国产 日韩一| 一级黄色大片毛片| 亚洲精品在线观看二区| 久久精品国产亚洲av香蕉五月 | 久久久精品免费免费高清| 国产精品久久久人人做人人爽| 老司机影院毛片| 最近最新免费中文字幕在线| 中文欧美无线码| aaaaa片日本免费| www.精华液| 久久 成人 亚洲| 老司机在亚洲福利影院| 午夜精品国产一区二区电影| av在线播放免费不卡| 中文字幕人妻丝袜一区二区| 国产男女超爽视频在线观看| 黄频高清免费视频| 日韩欧美免费精品| 18在线观看网站| 亚洲国产av影院在线观看| 亚洲avbb在线观看| 亚洲 欧美一区二区三区| 一区二区三区国产精品乱码| 一区二区三区乱码不卡18| 一区二区三区乱码不卡18| 日韩免费高清中文字幕av| 国产精品久久久久久精品电影小说| 97人妻天天添夜夜摸| 国产视频一区二区在线看| 女同久久另类99精品国产91| 天堂俺去俺来也www色官网| 每晚都被弄得嗷嗷叫到高潮| 自线自在国产av| 欧美av亚洲av综合av国产av| 久久中文看片网| 男女之事视频高清在线观看| 最近最新免费中文字幕在线| 国产精品秋霞免费鲁丝片| 十分钟在线观看高清视频www| 国产高清国产精品国产三级| 美国免费a级毛片| 精品国产超薄肉色丝袜足j| 韩国精品一区二区三区| 国产av一区二区精品久久| 国产免费现黄频在线看| 久久精品成人免费网站| 成人手机av| 国产av又大| 一本—道久久a久久精品蜜桃钙片| 国产高清videossex| xxxhd国产人妻xxx| 日日摸夜夜添夜夜添小说| 亚洲色图av天堂| 免费在线观看影片大全网站| 国产精品麻豆人妻色哟哟久久| 捣出白浆h1v1| 免费女性裸体啪啪无遮挡网站| 久久九九热精品免费| 狂野欧美激情性xxxx| 最近最新中文字幕大全电影3 | 免费黄频网站在线观看国产| 亚洲第一青青草原| videosex国产| 久久性视频一级片| 精品一区二区三区四区五区乱码| 国产精品欧美亚洲77777| 亚洲色图综合在线观看| 国产精品电影一区二区三区 | 制服诱惑二区| 777米奇影视久久| 高潮久久久久久久久久久不卡| 在线观看舔阴道视频| 亚洲中文字幕日韩| 涩涩av久久男人的天堂| 欧美在线黄色| 欧美老熟妇乱子伦牲交| 91麻豆精品激情在线观看国产 | 999久久久国产精品视频| 极品少妇高潮喷水抽搐| 免费少妇av软件| 久久久水蜜桃国产精品网| 亚洲av日韩在线播放| 国产极品粉嫩免费观看在线| 国产片内射在线| 久久精品亚洲精品国产色婷小说| 黄色 视频免费看| 国产主播在线观看一区二区| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费午夜福利视频| 久久久久精品国产欧美久久久| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品古装| 一区二区日韩欧美中文字幕| 欧美精品一区二区大全| 成人国产av品久久久| 日本vs欧美在线观看视频| 久久久久国产一级毛片高清牌| 午夜91福利影院| 欧美亚洲日本最大视频资源| 人妻 亚洲 视频| 久久人妻熟女aⅴ| 操美女的视频在线观看| 热re99久久国产66热| 亚洲伊人色综图| 亚洲人成77777在线视频| 91麻豆精品激情在线观看国产 | 久久人妻福利社区极品人妻图片| 纵有疾风起免费观看全集完整版| 桃花免费在线播放| 91九色精品人成在线观看| 水蜜桃什么品种好| 久久精品aⅴ一区二区三区四区| 天天操日日干夜夜撸| 一级片'在线观看视频| 丰满少妇做爰视频| 亚洲全国av大片| 五月开心婷婷网| 纵有疾风起免费观看全集完整版| 老司机午夜福利在线观看视频 | 国产日韩欧美亚洲二区| 国产男女超爽视频在线观看| 国产国语露脸激情在线看| 日日爽夜夜爽网站| 日韩免费av在线播放| av有码第一页| 国产国语露脸激情在线看| 成人精品一区二区免费| 丝瓜视频免费看黄片| 啦啦啦 在线观看视频| 中文字幕色久视频| 蜜桃国产av成人99| 18禁美女被吸乳视频| 91老司机精品| 久久国产精品人妻蜜桃| 老汉色∧v一级毛片| 老鸭窝网址在线观看| 无限看片的www在线观看| 99在线人妻在线中文字幕 | 天堂中文最新版在线下载| 免费一级毛片在线播放高清视频 | av网站免费在线观看视频| 我的亚洲天堂| 日韩精品免费视频一区二区三区| 亚洲成a人片在线一区二区| 欧美 日韩 精品 国产| 亚洲精品自拍成人| 国产免费视频播放在线视频| 亚洲av美国av| 久久香蕉激情| 侵犯人妻中文字幕一二三四区| 日韩欧美免费精品| 99精品在免费线老司机午夜| 国产不卡一卡二| 性少妇av在线| 国产精品久久久久久精品古装| av超薄肉色丝袜交足视频| 一区二区三区国产精品乱码| 久久av网站| 亚洲综合色网址| 久久人人97超碰香蕉20202| 最新的欧美精品一区二区| 午夜福利影视在线免费观看| 精品少妇黑人巨大在线播放| 老司机福利观看| 丝瓜视频免费看黄片| 波多野结衣一区麻豆| 亚洲欧美日韩另类电影网站| 日本a在线网址| 丁香六月欧美| 1024香蕉在线观看| 国产成人欧美在线观看 | 日本精品一区二区三区蜜桃| 国产免费福利视频在线观看| 国产高清videossex| 999精品在线视频| 亚洲九九香蕉| 久久ye,这里只有精品| 每晚都被弄得嗷嗷叫到高潮| 日韩熟女老妇一区二区性免费视频| 久久久精品国产亚洲av高清涩受| 一本久久精品| 国产精品香港三级国产av潘金莲| 亚洲色图 男人天堂 中文字幕| 一级毛片精品| 国产成人系列免费观看| 亚洲五月婷婷丁香| 看免费av毛片| 国产亚洲精品第一综合不卡| 亚洲精品中文字幕一二三四区 | 老汉色∧v一级毛片| 久久久水蜜桃国产精品网| 欧美日韩福利视频一区二区| 日日摸夜夜添夜夜添小说| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 在线av久久热| 欧美日韩av久久| 老汉色av国产亚洲站长工具| 久久久久久亚洲精品国产蜜桃av| 成年版毛片免费区| 国产av又大| 精品少妇久久久久久888优播| 精品人妻1区二区| 亚洲国产精品一区二区三区在线| 多毛熟女@视频| 飞空精品影院首页| 欧美精品一区二区大全| 窝窝影院91人妻| 国产精品麻豆人妻色哟哟久久| 纵有疾风起免费观看全集完整版| 国产成人一区二区三区免费视频网站| 欧美人与性动交α欧美软件| 久久毛片免费看一区二区三区| 无人区码免费观看不卡 | a级片在线免费高清观看视频| 人人妻人人澡人人看| 欧美日韩av久久| av有码第一页| 亚洲综合色网址| 国产欧美日韩一区二区三区在线| 日韩欧美免费精品| 久久热在线av| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻在线不人妻| 亚洲人成电影免费在线| 天天躁日日躁夜夜躁夜夜| 欧美在线黄色| 欧美成狂野欧美在线观看| 黑人操中国人逼视频| 日韩视频一区二区在线观看| 国产极品粉嫩免费观看在线| 香蕉久久夜色| kizo精华| 一区二区av电影网| 在线观看舔阴道视频| 日韩欧美一区二区三区在线观看 | 十八禁高潮呻吟视频| 一级片免费观看大全| 国产99久久九九免费精品| 免费少妇av软件| 免费人妻精品一区二区三区视频| 女人爽到高潮嗷嗷叫在线视频| 热99久久久久精品小说推荐| 操美女的视频在线观看| 变态另类成人亚洲欧美熟女 | 黄片小视频在线播放| 最黄视频免费看| 国产老妇伦熟女老妇高清| 男女床上黄色一级片免费看| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 久久中文字幕人妻熟女| 国产精品欧美亚洲77777| 成人av一区二区三区在线看| 亚洲av第一区精品v没综合| 成人影院久久| 嫩草影视91久久| 在线看a的网站| 免费在线观看日本一区| 搡老岳熟女国产| 国产成人免费观看mmmm| 黄网站色视频无遮挡免费观看| 精品免费久久久久久久清纯 | 国产成人av激情在线播放| 搡老岳熟女国产| 最黄视频免费看| 少妇 在线观看| 免费女性裸体啪啪无遮挡网站| 看免费av毛片| 国产淫语在线视频| 国产国语露脸激情在线看| 手机成人av网站| 国产欧美亚洲国产| 国产精品 欧美亚洲| 女人高潮潮喷娇喘18禁视频| 国产黄频视频在线观看| 我要看黄色一级片免费的| 国产一区二区激情短视频| 久久影院123| 午夜福利欧美成人| 老司机午夜十八禁免费视频| 亚洲三区欧美一区| 电影成人av| 午夜福利影视在线免费观看| 18禁美女被吸乳视频| 色婷婷久久久亚洲欧美| 天天添夜夜摸| 成人精品一区二区免费| 午夜福利在线免费观看网站| 成人18禁高潮啪啪吃奶动态图| 下体分泌物呈黄色| 天堂动漫精品| 久9热在线精品视频| 少妇猛男粗大的猛烈进出视频| 国产野战对白在线观看| 精品久久久久久电影网| 悠悠久久av| 久久久国产成人免费| 成人特级黄色片久久久久久久 | 黑人巨大精品欧美一区二区mp4| 在线看a的网站| 757午夜福利合集在线观看| 人人妻人人澡人人爽人人夜夜| 男女高潮啪啪啪动态图| 丁香欧美五月| 精品午夜福利视频在线观看一区 | 黄片小视频在线播放| 午夜激情久久久久久久| 狠狠精品人妻久久久久久综合| 免费看十八禁软件| 老司机深夜福利视频在线观看| 国产三级黄色录像| 免费av中文字幕在线| 国产不卡一卡二| 美女扒开内裤让男人捅视频| 国产高清激情床上av| 欧美一级毛片孕妇| 午夜福利欧美成人| √禁漫天堂资源中文www| 午夜福利欧美成人| 大片电影免费在线观看免费| 日韩精品免费视频一区二区三区| 欧美另类亚洲清纯唯美| 99re在线观看精品视频| 亚洲精品乱久久久久久| videos熟女内射| 精品一区二区三区四区五区乱码| 欧美 日韩 精品 国产| 桃红色精品国产亚洲av| 一级a爱视频在线免费观看| 国产免费福利视频在线观看| 中文字幕高清在线视频| 十八禁网站网址无遮挡| 久久人妻av系列| 日韩欧美一区视频在线观看| 中文字幕色久视频| 免费看a级黄色片| 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区视频在线观看免费 | 香蕉久久夜色| 女人久久www免费人成看片| 50天的宝宝边吃奶边哭怎么回事| 久热这里只有精品99| 精品午夜福利视频在线观看一区 | 在线天堂中文资源库| 国产一区二区 视频在线| 国产精品国产av在线观看| 欧美另类亚洲清纯唯美| 亚洲国产毛片av蜜桃av| 99久久精品国产亚洲精品| 男人舔女人的私密视频| 欧美日韩视频精品一区| 丰满饥渴人妻一区二区三| 午夜福利在线观看吧| 麻豆乱淫一区二区| 中文欧美无线码| 婷婷成人精品国产| 精品熟女少妇八av免费久了| 久久久久久免费高清国产稀缺| 这个男人来自地球电影免费观看| 黄色视频在线播放观看不卡| 大型av网站在线播放| 国产成+人综合+亚洲专区| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 精品人妻熟女毛片av久久网站| 成人18禁高潮啪啪吃奶动态图| kizo精华| 国产精品一区二区在线不卡| 欧美午夜高清在线| 欧美大码av| 人人妻人人爽人人添夜夜欢视频| 久久午夜综合久久蜜桃| kizo精华| 国产欧美日韩精品亚洲av| av天堂久久9| 国产一卡二卡三卡精品| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 久久久久久久久免费视频了| 精品人妻熟女毛片av久久网站| 久久99热这里只频精品6学生| www.999成人在线观看| 久久久久久久久免费视频了| 亚洲第一青青草原| 亚洲欧洲精品一区二区精品久久久| h视频一区二区三区| 午夜久久久在线观看| 国产精品香港三级国产av潘金莲| 2018国产大陆天天弄谢| 亚洲国产看品久久| 欧美日韩亚洲国产一区二区在线观看 | 国产黄频视频在线观看| 欧美日韩黄片免| 国产成人一区二区三区免费视频网站| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看 | 久久久精品国产亚洲av高清涩受| 国产无遮挡羞羞视频在线观看|