• <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
    亚洲精品一二三| 一级a做视频免费观看| 欧美最新免费一区二区三区| 亚洲三级黄色毛片| 中国美白少妇内射xxxbb| 91久久精品电影网| 2018国产大陆天天弄谢| 亚洲国产色片| 人妻 亚洲 视频| 永久网站在线| 一边亲一边摸免费视频| 国产成人aa在线观看| 国产深夜福利视频在线观看| 不卡视频在线观看欧美| freevideosex欧美| 最近的中文字幕免费完整| 建设人人有责人人尽责人人享有的| 青春草视频在线免费观看| 久久国产精品大桥未久av | 一区二区三区四区激情视频| 亚洲成人一二三区av| 亚洲人成网站在线观看播放| 日日啪夜夜撸| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品电影小说| 高清毛片免费看| 少妇高潮的动态图| 久久精品国产自在天天线| 偷拍熟女少妇极品色| 一级黄片播放器| 久久这里有精品视频免费| 亚洲国产精品一区二区三区在线| 少妇的逼好多水| 国内精品宾馆在线| av专区在线播放| 三级国产精品欧美在线观看| 大片电影免费在线观看免费| 欧美精品亚洲一区二区| 精品一品国产午夜福利视频| 久久热精品热| 美女xxoo啪啪120秒动态图| 韩国av在线不卡| a级一级毛片免费在线观看| 中文字幕人妻丝袜制服| 一本—道久久a久久精品蜜桃钙片| 国产精品久久久久久av不卡| 丝袜喷水一区| 交换朋友夫妻互换小说| 国产亚洲午夜精品一区二区久久| 国产精品99久久99久久久不卡 | 亚洲精品aⅴ在线观看| 一本—道久久a久久精品蜜桃钙片| 一级毛片我不卡| 久久精品国产a三级三级三级| 日韩熟女老妇一区二区性免费视频| 少妇被粗大猛烈的视频| 精品久久久精品久久久| 亚洲三级黄色毛片| 一区二区三区精品91| 久久ye,这里只有精品| 97精品久久久久久久久久精品| 久久精品熟女亚洲av麻豆精品| 在线看a的网站| 国产精品人妻久久久影院| 亚洲欧美清纯卡通| 亚洲丝袜综合中文字幕| 午夜日本视频在线| 久久99热这里只频精品6学生| 国产精品久久久久久精品古装| 麻豆成人av视频| 中文欧美无线码| 日韩成人av中文字幕在线观看| 在线观看www视频免费| 国产精品国产av在线观看| 亚洲av中文av极速乱| 丰满乱子伦码专区| 久久精品久久久久久久性| 亚洲国产精品999| 免费人妻精品一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 99热这里只有精品一区| 女人久久www免费人成看片| 日韩精品有码人妻一区| 日韩免费高清中文字幕av| 国产亚洲av片在线观看秒播厂| 亚洲av二区三区四区| 亚洲怡红院男人天堂| 亚洲国产最新在线播放| 精品一区二区三卡| 国产黄色免费在线视频| 高清不卡的av网站| 精品亚洲成国产av| 在线观看av片永久免费下载| 久久人人爽av亚洲精品天堂| 亚洲av在线观看美女高潮| 久久久久久久久久人人人人人人| 国产一区亚洲一区在线观看| 秋霞在线观看毛片| 国产成人精品福利久久| 亚州av有码| 老司机影院成人| av专区在线播放| 亚洲真实伦在线观看| 亚洲精品久久午夜乱码| 在线观看人妻少妇| 亚洲精品视频女| 亚洲av电影在线观看一区二区三区| 久久精品久久久久久久性| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 国产免费一区二区三区四区乱码| 欧美97在线视频| 男女啪啪激烈高潮av片| 一区二区三区四区激情视频| 国产免费又黄又爽又色| 只有这里有精品99| 日日摸夜夜添夜夜爱| 99久久人妻综合| av.在线天堂| 2018国产大陆天天弄谢| av女优亚洲男人天堂| 久久久精品免费免费高清| 这个男人来自地球电影免费观看 | 99热全是精品| 三级经典国产精品| 2018国产大陆天天弄谢| 亚洲国产欧美日韩在线播放 | 男人爽女人下面视频在线观看| 99久久精品一区二区三区| 免费高清在线观看视频在线观看| 热re99久久国产66热| 日本91视频免费播放| 两个人的视频大全免费| 免费观看性生交大片5| 黑人巨大精品欧美一区二区蜜桃 | 看十八女毛片水多多多| 亚洲内射少妇av| 建设人人有责人人尽责人人享有的| 免费大片黄手机在线观看| 色婷婷av一区二区三区视频| 亚洲精品第二区| 内射极品少妇av片p| 国内揄拍国产精品人妻在线| 麻豆成人午夜福利视频| 亚洲中文av在线| 十分钟在线观看高清视频www | 最近中文字幕2019免费版| 国产一区亚洲一区在线观看| 老司机影院成人| 久久久久久久亚洲中文字幕| 国产精品久久久久成人av| 男人爽女人下面视频在线观看| 日韩 亚洲 欧美在线| 免费看日本二区| 精品久久久噜噜| 欧美+日韩+精品| 一级毛片黄色毛片免费观看视频| 国产永久视频网站| 久久久久久伊人网av| 免费大片18禁| 热99国产精品久久久久久7| 男人狂女人下面高潮的视频| 美女中出高潮动态图| 色94色欧美一区二区| 亚洲真实伦在线观看| 最近最新中文字幕免费大全7| 成人影院久久| 国产精品久久久久久精品电影小说| 亚洲不卡免费看| 亚洲欧美一区二区三区黑人 | 久久这里有精品视频免费| 成人黄色视频免费在线看| 国模一区二区三区四区视频| 日韩成人伦理影院| 国产欧美另类精品又又久久亚洲欧美| 日韩一区二区三区影片| 亚洲av二区三区四区| 一边亲一边摸免费视频| 亚洲图色成人| 国产精品偷伦视频观看了| 精品少妇内射三级| 在线播放无遮挡| 亚洲欧美成人综合另类久久久| 热99国产精品久久久久久7| 色5月婷婷丁香| 国产欧美日韩综合在线一区二区 | 色94色欧美一区二区| 国产国拍精品亚洲av在线观看| 欧美丝袜亚洲另类| 纵有疾风起免费观看全集完整版| 六月丁香七月| 国产伦在线观看视频一区| 亚洲综合色惰| 美女视频免费永久观看网站| 欧美精品人与动牲交sv欧美| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 中文欧美无线码| 中文字幕精品免费在线观看视频 | 久久精品熟女亚洲av麻豆精品| 亚洲精品日韩av片在线观看| 国产成人aa在线观看| av有码第一页| 国产视频首页在线观看| 亚洲国产毛片av蜜桃av| 少妇丰满av| 久久久久网色| .国产精品久久| 亚洲精品国产色婷婷电影| 日日啪夜夜爽| 深夜a级毛片| 国内精品宾馆在线| 国产黄片美女视频| 亚洲av电影在线观看一区二区三区| 日日啪夜夜爽| 看非洲黑人一级黄片| 在现免费观看毛片| 色网站视频免费| 搡老乐熟女国产| 亚洲国产精品国产精品| 亚洲精品视频女| 99九九在线精品视频 | 国产免费视频播放在线视频| 日韩一区二区视频免费看| 3wmmmm亚洲av在线观看| 国产乱来视频区| 欧美日韩综合久久久久久| 国产日韩一区二区三区精品不卡 | 亚洲欧美成人精品一区二区| 日韩大片免费观看网站| 国产91av在线免费观看| 如日韩欧美国产精品一区二区三区 | 日韩一本色道免费dvd| 91精品伊人久久大香线蕉| 成人午夜精彩视频在线观看| 精品久久久久久电影网| 欧美成人精品欧美一级黄| 国产女主播在线喷水免费视频网站| 高清不卡的av网站| 久久久久久久久久成人| 一级毛片电影观看| 美女中出高潮动态图| 亚洲av成人精品一二三区| 秋霞在线观看毛片| 永久网站在线| 一本久久精品| 波野结衣二区三区在线| 嫩草影院新地址| 日韩精品免费视频一区二区三区 | 黑丝袜美女国产一区| 丰满迷人的少妇在线观看| 又黄又爽又刺激的免费视频.| 深夜a级毛片| 欧美激情极品国产一区二区三区 | 亚洲熟女精品中文字幕| 韩国av在线不卡| 欧美xxⅹ黑人| 午夜福利影视在线免费观看| 免费av不卡在线播放| 人人澡人人妻人| 97在线人人人人妻| 精品熟女少妇av免费看| 日本wwww免费看| 精品久久久久久久久亚洲| 成人黄色视频免费在线看| 日韩免费高清中文字幕av| 久热这里只有精品99| 午夜av观看不卡| 性色avwww在线观看| 国产成人精品一,二区| 91久久精品电影网| 最近中文字幕高清免费大全6| 亚洲第一av免费看| 亚洲av国产av综合av卡| 国产精品秋霞免费鲁丝片| 我要看黄色一级片免费的| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看| 简卡轻食公司| 国产成人一区二区在线| 草草在线视频免费看| 简卡轻食公司| 久久久久久久久久人人人人人人| 午夜av观看不卡| 欧美 亚洲 国产 日韩一| 国产精品99久久久久久久久| 97超视频在线观看视频| 男人狂女人下面高潮的视频| 永久网站在线| 久久97久久精品| tube8黄色片| 亚洲熟女精品中文字幕| 一本久久精品| 99re6热这里在线精品视频| 久久久国产一区二区| 国产精品伦人一区二区| 寂寞人妻少妇视频99o| 天堂8中文在线网| 久久精品国产鲁丝片午夜精品| 这个男人来自地球电影免费观看 | 精华霜和精华液先用哪个| 亚洲精品国产av成人精品| 日韩在线高清观看一区二区三区| 日韩不卡一区二区三区视频在线| 国产成人aa在线观看| 大陆偷拍与自拍| 国产亚洲91精品色在线| 亚洲av不卡在线观看| 日韩欧美 国产精品| 亚洲欧美日韩卡通动漫| 精品一品国产午夜福利视频| 免费播放大片免费观看视频在线观看| 亚洲国产色片| 国产在线免费精品| 一边亲一边摸免费视频| 亚洲av综合色区一区| 夜夜看夜夜爽夜夜摸| av在线播放精品| 午夜免费观看性视频| 菩萨蛮人人尽说江南好唐韦庄| 国产乱来视频区| 精品熟女少妇av免费看| 亚洲一级一片aⅴ在线观看| 免费看av在线观看网站| 伊人久久精品亚洲午夜| 一区二区av电影网| 午夜老司机福利剧场| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 99久久综合免费| 中文字幕久久专区| 最黄视频免费看| 王馨瑶露胸无遮挡在线观看| 老女人水多毛片| 日韩欧美 国产精品| 99九九在线精品视频 | 欧美国产精品一级二级三级 | 亚洲国产av新网站| 97精品久久久久久久久久精品| 黄色日韩在线| 国产精品三级大全| 亚洲精品一区蜜桃| 自拍偷自拍亚洲精品老妇| 久久午夜综合久久蜜桃| 午夜福利视频精品| 久久99热6这里只有精品| 免费少妇av软件| 插逼视频在线观看| 国产视频首页在线观看| 中文字幕制服av| 国产精品一区二区在线观看99| 国产国拍精品亚洲av在线观看| 美女视频免费永久观看网站| 熟女电影av网| av一本久久久久| 三上悠亚av全集在线观看 | 91精品伊人久久大香线蕉| 纯流量卡能插随身wifi吗| 精品国产一区二区三区久久久樱花| 国产男人的电影天堂91| 国产精品久久久久久久久免| 国产成人精品久久久久久| 成人特级av手机在线观看| 人妻系列 视频| 欧美激情国产日韩精品一区| 日韩一本色道免费dvd| √禁漫天堂资源中文www| 欧美精品亚洲一区二区| 久久久久国产网址| 精华霜和精华液先用哪个| 一本久久精品| 99久久精品一区二区三区| 日韩中字成人| 女性被躁到高潮视频| 精品人妻熟女毛片av久久网站| av免费在线看不卡| 亚洲精品第二区| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 少妇 在线观看| 亚洲在久久综合| 亚洲国产欧美在线一区| 免费看av在线观看网站| 天堂中文最新版在线下载| 99九九线精品视频在线观看视频| 91成人精品电影| 亚洲精品自拍成人| 一本一本综合久久| 日日摸夜夜添夜夜添av毛片| 亚洲内射少妇av| 国产高清三级在线| 最黄视频免费看| 少妇人妻精品综合一区二区| 国产精品一区二区在线观看99| 午夜激情久久久久久久| 欧美成人午夜免费资源| 制服丝袜香蕉在线| 亚洲欧美日韩卡通动漫| 91成人精品电影| 性高湖久久久久久久久免费观看| 最新中文字幕久久久久| 春色校园在线视频观看| 黄色欧美视频在线观看| 汤姆久久久久久久影院中文字幕| 免费观看性生交大片5| 亚洲国产精品国产精品| 国产永久视频网站| 嫩草影院新地址| 男女无遮挡免费网站观看| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说| 人妻 亚洲 视频| 菩萨蛮人人尽说江南好唐韦庄| 韩国av在线不卡| 欧美xxⅹ黑人| 日日爽夜夜爽网站| 欧美变态另类bdsm刘玥| 99久久中文字幕三级久久日本| 国产视频内射| 永久网站在线| 91精品伊人久久大香线蕉| 国产精品不卡视频一区二区| 亚州av有码| 国产69精品久久久久777片| 人人妻人人澡人人爽人人夜夜| av天堂久久9| 久久精品国产鲁丝片午夜精品| 日本av免费视频播放| 成年人午夜在线观看视频| 久久女婷五月综合色啪小说| 噜噜噜噜噜久久久久久91| 久久久久久久国产电影| 天天操日日干夜夜撸| 在线观看美女被高潮喷水网站| 国产乱人偷精品视频| 91久久精品电影网| 中文字幕av电影在线播放| 一个人免费看片子| 国产真实伦视频高清在线观看| 黄色视频在线播放观看不卡| .国产精品久久| 国产亚洲最大av| 国产日韩欧美亚洲二区| 国产精品国产三级专区第一集| 午夜福利影视在线免费观看| 在线播放无遮挡| 国产国拍精品亚洲av在线观看| 国产深夜福利视频在线观看| 国产 一区精品| 丝袜脚勾引网站| 久久国内精品自在自线图片| 大片电影免费在线观看免费| 乱码一卡2卡4卡精品| 亚洲怡红院男人天堂| 亚洲av在线观看美女高潮| 久久精品久久久久久噜噜老黄| 亚洲国产成人一精品久久久| 一级,二级,三级黄色视频| 人妻少妇偷人精品九色| 亚洲精品aⅴ在线观看| 麻豆成人午夜福利视频| 成人毛片a级毛片在线播放| 国产老妇伦熟女老妇高清| 全区人妻精品视频| 精品国产国语对白av| 99精国产麻豆久久婷婷| 亚洲欧美清纯卡通| 欧美区成人在线视频| 一本久久精品| 久久久久久久久久人人人人人人| 热re99久久国产66热| 老司机影院毛片| 美女内射精品一级片tv| 波野结衣二区三区在线| videossex国产| 国产av精品麻豆| 99久久综合免费| 成人免费观看视频高清| 99久久综合免费| 夫妻性生交免费视频一级片| 久久6这里有精品| 亚洲人成网站在线观看播放| 国产亚洲91精品色在线| 久久久久久人妻| 亚洲伊人久久精品综合| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线| 亚洲欧美日韩另类电影网站| 亚洲激情五月婷婷啪啪| 国产亚洲精品久久久com| 国产av码专区亚洲av| av国产久精品久网站免费入址| 免费av不卡在线播放| 精品卡一卡二卡四卡免费| 波野结衣二区三区在线| 亚洲欧美清纯卡通| 午夜福利在线观看免费完整高清在| 99九九在线精品视频 | 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 婷婷色综合大香蕉| 久久久久久人妻| 日韩不卡一区二区三区视频在线| 亚洲av不卡在线观看| 国产极品粉嫩免费观看在线 | 亚洲精品成人av观看孕妇| freevideosex欧美| 国产无遮挡羞羞视频在线观看| 亚洲精品日韩av片在线观看| 免费黄网站久久成人精品| 日韩伦理黄色片| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| av视频免费观看在线观看| 在线观看美女被高潮喷水网站| 大香蕉久久网| 热re99久久国产66热| 免费观看av网站的网址| 曰老女人黄片| 在线观看免费日韩欧美大片 | 亚洲国产毛片av蜜桃av| 国产成人91sexporn| 亚洲第一av免费看| 亚洲欧洲精品一区二区精品久久久 | 亚洲怡红院男人天堂| 内地一区二区视频在线| 人人妻人人澡人人看| 精品亚洲成国产av| 日韩av免费高清视频| 少妇被粗大的猛进出69影院 | 亚洲无线观看免费| 国产黄片视频在线免费观看| 国产白丝娇喘喷水9色精品| 久久午夜福利片| av福利片在线观看| 99久国产av精品国产电影| 国内少妇人妻偷人精品xxx网站| 秋霞在线观看毛片| 成人黄色视频免费在线看| 另类精品久久| 少妇人妻一区二区三区视频| 亚洲精品中文字幕在线视频 | 乱码一卡2卡4卡精品| 丰满人妻一区二区三区视频av| 国产高清三级在线| 亚洲精品一二三| 精华霜和精华液先用哪个| 特大巨黑吊av在线直播| 国产极品粉嫩免费观看在线 | 卡戴珊不雅视频在线播放| 国产精品一区二区三区四区免费观看| 免费观看无遮挡的男女| 黑丝袜美女国产一区| 一本久久精品| 日本-黄色视频高清免费观看| 在现免费观看毛片| 欧美精品一区二区大全| 国产精品久久久久久精品电影小说| 寂寞人妻少妇视频99o| 老司机影院成人| 美女大奶头黄色视频| 中文乱码字字幕精品一区二区三区| 成人亚洲欧美一区二区av| 一区二区三区乱码不卡18| 插逼视频在线观看| av不卡在线播放| 久久99一区二区三区| 色婷婷久久久亚洲欧美| 三级国产精品欧美在线观看| 国产精品久久久久久精品电影小说| 国产永久视频网站| 色视频www国产| 看非洲黑人一级黄片| 又爽又黄a免费视频| 国产成人精品婷婷| 国产精品不卡视频一区二区| 91在线精品国自产拍蜜月| 亚洲经典国产精华液单| 国产精品嫩草影院av在线观看| 男女啪啪激烈高潮av片| 热re99久久精品国产66热6| 成人国产麻豆网| 免费少妇av软件| 久久ye,这里只有精品| 国产精品99久久久久久久久| 亚洲成人手机| 人妻夜夜爽99麻豆av| 国产精品国产三级国产专区5o| av福利片在线| 丁香六月天网| 日本猛色少妇xxxxx猛交久久| 久久av网站| 久久午夜福利片| 国产成人精品无人区| 国产视频内射| 成人国产麻豆网| 一级爰片在线观看| 中国国产av一级| 少妇高潮的动态图| 一区二区三区乱码不卡18| 国产成人精品久久久久久| 熟女人妻精品中文字幕| 亚洲国产精品一区二区三区在线| 国内精品宾馆在线| 纵有疾风起免费观看全集完整版| 国产片特级美女逼逼视频| 亚洲激情五月婷婷啪啪| 久久精品久久久久久久性| 少妇人妻 视频| 新久久久久国产一级毛片| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 国产高清有码在线观看视频| 美女脱内裤让男人舔精品视频| 国产精品99久久99久久久不卡 |