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

    阜平縣地表溫度遙感反演及其空間特征研究

    2021-04-06 04:08:28李春波任磊于蒙李會(huì)愷焦尚斌
    城市地質(zhì) 2021年1期
    關(guān)鍵詞:阜平縣坡向海拔

    李春波 任磊 于蒙 李會(huì)愷 焦尚斌

    摘 要:以阜平縣為研究區(qū)域,基于Landsat 8遙感影像和ASTGTM2 DEM數(shù)據(jù),利用輻射傳輸方程法反演地表溫度,運(yùn)用統(tǒng)計(jì)學(xué)方法分析地表溫度與歸一化植被指數(shù)(NDVI)、海拔和坡向之間的關(guān)系。結(jié)果表明:阜平縣地表溫度整體呈西高東低分布,地表溫度與植被覆蓋、地形要素之間的關(guān)系密切,地表溫度高值區(qū)域主要分布在研究區(qū)西部和東北部海拔較高且NDVI值較高區(qū)域,低值區(qū)域主要分布在東南部海拔較低且NDVI值較低區(qū)域。地表溫度隨海拔和NDVI的升高逐漸減小,呈現(xiàn)明顯的負(fù)線性相關(guān);不同坡向的地表溫度分布具有一定的差異,陽(yáng)坡的地表溫度高于陰坡。

    關(guān)鍵詞:地表溫度反演;輻射傳輸方程法;影響因素;阜平縣

    Abstract: Taking Fuping County as the research area, based on Landsat 8 remote sensing images and ASTGTM2 DEM data, this paper applies radiation transfer equation method to invert the surface temperature of the study area, and statistical methods to analyze the relationship between land surface temperature and normalized vegetation index, elevation and aspect. The results show that the land surface temperature in Fuping County is high in the west and low in the east. The surface temperature and vegetation coverage and topographic elements are closely related. The high surface temperature areas are mainly distributed in the west and northeast of the study area with higher elevations and higher NDVI values. The low surface temperature areas are mainly distributed in the southeastern areas with lower elevations and lower NDVI values. The statistic results show that the surface temperature gradually decreases with elevation and NDVI, showing a significant negative correlation; the surface temperature distribution of different aspects varies, with the surface temperature of the sunny slope higher than that of the shaded slope.

    Keywords: surface temperature inversion; radiation transfer equation method; influencing factors; Fuping County

    地表溫度(Land Surface Temperature,LST)是影響自然環(huán)境形成與變化的最活躍因素之一,同時(shí)也是地表能量平衡的關(guān)鍵參數(shù)之一,研究地表溫度的分布特征與變化具有重要的意義(史新等,2018)。目前,基于遙感技術(shù)反演地表溫度已廣泛應(yīng)用于城市熱島效應(yīng)研究、地?zé)豳Y源預(yù)測(cè)、區(qū)域干旱監(jiān)測(cè)等諸多領(lǐng)域(李學(xué)敏等,2019;何雪琴等,2019;李紅軍等,2006),部分學(xué)者還進(jìn)行了地表溫度影響因素的研究。伍顯等(2018)基于Landsat 8數(shù)據(jù)和DEM數(shù)據(jù)研究了貴州思雅河流域地表溫度與土地利用、地形要素的關(guān)系,孫常峰等(2014)基于Landsat 5數(shù)據(jù)和ASTGTM2 DEM數(shù)據(jù)研究了泰山地表溫度與地形因子、歸一化水汽指數(shù)(Normalized Different Moisture Index,NDMI)、歸一化植被指數(shù)(Normalized Vegetation Index,NDVI)的關(guān)系,熊俊楠等(2019)基于Landsat 8數(shù)據(jù)和SRTM DEM數(shù)據(jù)研究了青藏高原桑珠孜區(qū)地表溫度與地形因子、NDMI、NDVI之間的關(guān)系。這些研究表明,地表溫度與地形、植被覆蓋等因素之間存在一定的關(guān)系,但不同研究給出的結(jié)論略有不同,這說(shuō)明在不同地區(qū)地表溫度的影響因素對(duì)地表溫度的影響程度存在差異。

    本文以阜平縣為研究區(qū)域,基于Landsat 8遙感影像反演研究區(qū)的地表溫度、NDVI,基于ASTGTM2 DEM數(shù)據(jù)計(jì)算坡向,運(yùn)用統(tǒng)計(jì)學(xué)方法分析地表溫度與NDVI、海拔、坡向之間的關(guān)系,為下一步開展阜平縣地?zé)豳Y源預(yù)測(cè)提供依據(jù)。

    1 研究區(qū)概況與數(shù)據(jù)

    1.1 研究區(qū)概況

    阜平縣位于保定市西部,面積約2500 km2,地理坐標(biāo)為東經(jīng) 113°46′~114°32′、北偉 38°39′~39°7′(圖1)。該縣地處太行山中北部東麓,為全山區(qū)縣,縣域內(nèi)地形復(fù)雜,山巒綿亙、溝壑縱橫,地勢(shì)由東南向西北逐漸升高,海拔高度一般為200~2200 m(圖2)。氣候?qū)倥瘻貛О霛駶?rùn)半干旱大陸性季風(fēng)氣候,四季分明,年平均氣溫為12.6 ℃,年平均降水量630 mm左右,年平均蒸發(fā)量2024 mm左右,無(wú)霜期約225 d??h域生態(tài)環(huán)境整體較好,森林覆蓋率達(dá)39.47%,植被覆蓋度達(dá)80%以上。阜平縣地?zé)豳Y源豐富,其中城南莊溫泉、溫塘寺溫泉、下堡溫泉等目前已處于開發(fā)利用階段。

    1.2 數(shù)據(jù)及預(yù)處理

    研究數(shù)據(jù)包括Landsat 8 OLI/TIRS遙感影像和ASTGTM2 DEM數(shù)據(jù),均來(lái)源于地理空間數(shù)據(jù)云平臺(tái)(http://www.gscloud.cn)。其中Landsat 8遙感影像1景,標(biāo)識(shí)號(hào)為L(zhǎng)C81240332015225LGN01,條帶號(hào)124、行編號(hào)33,平均云量為0.09%,獲取時(shí)間為2015年8月13日;ASTGTM2 DEM數(shù)據(jù)4景,條帶號(hào)113~114,行編號(hào)38~39,空間分辨率為30 m。

    Landsat 8數(shù)據(jù)預(yù)處理主要包括輻射定標(biāo)、大氣校正、影像裁剪等,ASTGTM2數(shù)據(jù)預(yù)處理主要包括數(shù)據(jù)鑲嵌、投影、裁剪等,主要用于地表溫度反演以及計(jì)算NDVI、坡向等。

    2 地表溫度反演

    基于遙感技術(shù)反演地表溫度的方法大致可分為輻射傳輸方程法、單窗算法、單通道算法、劈窗算法等。侯宇初等(2019)基于Landsat 8數(shù)據(jù)在大興安嶺地區(qū)對(duì)比分析了輻射傳輸方程法、單窗算法和劈窗算法的地表溫度反演精度,認(rèn)為輻射傳輸方程法反演的結(jié)果精度更高;史新等(2018)利用Landsat 8數(shù)據(jù)在廣東三河壩流域?qū)Ρ攘溯椛鋫鬏敺匠谭?、單窗算法和單通道算?種算法,認(rèn)為3種算法的反演結(jié)果線性擬合程度類似,空間分布一致;Yu等(2014)分析了輻射傳輸方程法、分裂窗算法和單通道算法3種算法的精度,認(rèn)為輻射傳輸方程法的精度高于另外2種算法。結(jié)合上述研究,本次研究采用輻射傳輸方程法反演阜平縣地表溫度。

    輻射傳輸方程法的基本原理是將衛(wèi)星傳感器接收到的總輻射分為3部分,分別是地表熱輻射經(jīng)大氣衰減后到達(dá)傳感器的輻射、大氣上行輻射和大氣下行輻射經(jīng)地表反射后通過(guò)大氣傳輸?shù)絺鞲衅鞯妮椛洌ㄠ嚂螅?014)。在假設(shè)地表為朗伯體的條件下,輻射傳輸方程可簡(jiǎn)化為:

    式中:Lλ為傳感器接收到的波長(zhǎng)為λ的熱紅外輻射強(qiáng)度,Bλ (Ts )為溫度Ts時(shí)黑體的的輻射強(qiáng)度,ελ是地表比輻射率,τλ為大氣透過(guò)率,、? 分別為大氣上行輻射強(qiáng)度和大氣下行輻射強(qiáng)度。

    普朗克定律闡明了物體熱輻射強(qiáng)度與波長(zhǎng)、溫度之間的關(guān)系,根據(jù)該定律,可通過(guò)黑體在熱紅外波段的輻射強(qiáng)度反演地表溫度:

    式中:K1、K2為熱紅外傳感器的定標(biāo)參數(shù),對(duì)于Landsat 8 TIRS band 10波段,K1 = 774.89 W·m-2·μm-1·sr-1,K2=1321.08 K。B(Ts )可由公式(1)推導(dǎo)出:

    由公式(2)、公式(3)可知,利用輻射傳輸方程法反演地表溫度必須的參數(shù)為地表比輻射率ε和大氣透過(guò)率τ、大氣上行輻射強(qiáng)度L↑和大氣下行輻射強(qiáng)度L↓等大氣剖面參數(shù)。對(duì)于Landsat 8數(shù)據(jù)可在NASA提供的網(wǎng)站(http://atmcorr.gsfc.nasa.gov)輸入成像時(shí)間以及中心經(jīng)緯度等信息獲取大氣剖面參數(shù)。地表比輻射率的計(jì)算可根據(jù)公式(4)計(jì)算(Sobrino et al.,2004):

    3 結(jié)果分析

    阜平縣地表溫度為20.22~55.07 ℃(圖3),低值區(qū)域(小于30 ℃)占比10.57%,高值區(qū)域(大于40 ℃)占比12.47%。地表溫度低值區(qū)域主要分布在研究區(qū)西部和東北部海拔較高、地形起伏較大且NDVI值較高的區(qū)域,另外在研究區(qū)東南部的水體區(qū)域地表溫度值也較低;地表溫度高值區(qū)域主要分布在研究區(qū)東南部海拔較低且地形相對(duì)較為平坦、城鎮(zhèn)分布較為集中、NDVI值較低區(qū)域,地表溫度與地形、植被覆蓋之間的關(guān)系較為明顯(圖4)。

    為進(jìn)一步分析地表溫度與植被、地形之間的關(guān)系,參考相關(guān)研究,選取NDVI、海拔和坡向3個(gè)地表溫度影響因素,分別運(yùn)用隨機(jī)點(diǎn)統(tǒng)計(jì)、分區(qū)統(tǒng)計(jì)等方法,分析地表溫度與NDVI、海拔和坡向之間的關(guān)系。

    3.1 隨機(jī)點(diǎn)統(tǒng)計(jì)

    在研究區(qū)范圍內(nèi)隨機(jī)選取3000個(gè)點(diǎn),在ArcGIS 10.2軟件中分別提取地表溫度(LST)、NDVI、海拔高度和坡向的數(shù)值,生成散點(diǎn)圖,并進(jìn)行相關(guān)性分析。

    相關(guān)研究表明,NDVI為負(fù)值的區(qū)域一般為水體,在分析LST與NDVI之間的關(guān)系時(shí)應(yīng)該去除,因此本次研究?jī)H分析NDVI為正值的數(shù)據(jù)。隨機(jī)點(diǎn)分析結(jié)果顯示,LST隨NDVI的增大而逐漸減小(圖5),LST與NDVI之間呈現(xiàn)明顯的負(fù)相關(guān)。LST與地形之間的關(guān)系較為復(fù)雜,LST與海拔高程呈現(xiàn)明顯的負(fù)線性相關(guān)(圖6),但在海拔300 m以下,LST隨海拔的升高呈現(xiàn)增加的趨勢(shì),這是因?yàn)檠芯繀^(qū)海拔高程最低的位置位于研究區(qū)東南部的水體區(qū)域,而水體的LST小于周邊區(qū)域且海拔高程也低于周邊造成的。LST與坡向的關(guān)系與一般認(rèn)知相同,即陽(yáng)坡的地表溫度比陰坡的地表溫度稍高(圖7)。

    3.2 分區(qū)統(tǒng)計(jì)

    為進(jìn)一步分析LST與NDVI、海拔和坡向之間的關(guān)系,運(yùn)用分區(qū)統(tǒng)計(jì)的方法,將NDVI、海拔、坡度和坡向的數(shù)據(jù)按照一定的增量劃分為若干區(qū)間,統(tǒng)計(jì)區(qū)間內(nèi)LST的平均值,并進(jìn)行相關(guān)性分析。其中,NDVI數(shù)據(jù)按照每增加0.1劃分1個(gè)區(qū)間,共劃分10個(gè)區(qū)間;海拔數(shù)據(jù)按照每增加100 m劃分1個(gè)區(qū)間,共劃分20個(gè)區(qū)間;坡向以正北方向?yàn)?,按照順時(shí)針?lè)较驅(qū)⑵孪驍?shù)據(jù)劃分為平面、北坡(337.5°~22.5°)、東北坡(22.5°~67.5°)、東坡(67.5°~112.5°)、東南坡(112.5°~157.5°)、南坡(157.5°~202.5°)、西南坡(202.5°~247.5°)、西坡(247.5°~292.5°)、西北坡(292.5°~337.5°)9大類(圖8),一般而言,朝向南面的坡稱為陽(yáng)坡,朝向背面的坡稱為陰坡。

    結(jié)果表明,分區(qū)統(tǒng)計(jì)與隨機(jī)點(diǎn)統(tǒng)計(jì)結(jié)果具有一定的相似性。以NDVI增加0.1為區(qū)間,LST平均值隨NDVI的增大逐漸減小,表明LST的平均值與NDVI呈負(fù)線性相關(guān),且相關(guān)系數(shù)R2達(dá)0.99(圖9)。以海拔高度每增加100 m為區(qū)間,LST平均值與海拔呈負(fù)線性相關(guān),相關(guān)系數(shù)R2為0.99(圖10)。不同坡向的LST平均值差異較大,陽(yáng)坡的地表溫度高于陰坡的地表溫度,最大相差約2 ℃,其他坡度的差異相對(duì)較小,表明對(duì)于山區(qū)陰陽(yáng)坡向的不同也是影響地表溫度分布的一個(gè)重要因素(圖11)。

    4 結(jié)論

    (1)阜平縣地表溫度為20.22~55.07 ℃,低值區(qū)域主要分布在研究區(qū)西部和東北部海拔較高、地形起伏較大且NDVI值較高的區(qū)域,高值區(qū)域主要分布在研究區(qū)東南部海拔較低且地形相對(duì)較為平坦、城鎮(zhèn)分布較為集中、NDVI值較低區(qū)域。

    (2)阜平縣地表溫度與植被覆蓋、地形要素之間的關(guān)系較為密切。地表溫度與NDVI呈負(fù)相關(guān),地表溫度隨NDVI的增加逐漸減小;地表溫度與海拔之間呈負(fù)相關(guān),當(dāng)海拔高度大于300 m時(shí),每升高100 m,地表溫度下降約0.84 ℃;不同坡向的地表溫度分布具有一定的差異,陽(yáng)坡的地表溫度高于陰坡。

    (3)本次研究?jī)H分析了植被覆蓋和地形與地表溫度之間的關(guān)系,實(shí)際上地表溫度是在太陽(yáng)輻射、季節(jié)變化、地表覆蓋、人類活動(dòng)等多種因素下共同作用的結(jié)果,地表溫度與其他影響因素之間的關(guān)系還有待進(jìn)一步研究。

    參考文獻(xiàn):

    鄧書斌,2014. ENVI遙感圖像處理方法[M]. 2版. 北京:高等教育出版社.

    何雪琴,鄧清海,高宗軍,等,2019. 遙感技術(shù)在干熱巖選區(qū)中的應(yīng)用:以魯東地區(qū)為例[J]. 山東國(guó)土資源,35(6):67-73.

    侯宇初,張冬有,2019. 基于Landsat8遙感影像的地表溫度反演方法對(duì)比研究[J]. 中國(guó)農(nóng)學(xué)通報(bào),35(10):142-147.

    李紅軍,鄭力,雷玉平,等,2006. 植被指數(shù)-地表溫度特征空間研究及其在旱情監(jiān)測(cè)中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào),22(11):170-174.

    李學(xué)敏,文力,王儷璇,2019. 基于Landsat遙感數(shù)據(jù)的城市熱島效應(yīng)及影響因素:以襄陽(yáng)市為例[J]. 湖北農(nóng)業(yè)科學(xué),58(16):41-48.

    孫常峰,孔繁華,尹海偉,等,2014. 山區(qū)夏季地表溫度的影響因素:以泰山為例[J].? 生態(tài)學(xué)報(bào),34(12):3394-3404.

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

    史新,周買春,劉振華,等,2018. 基于Landsat 8數(shù)據(jù)的3種地表溫度反演算法在三河壩流域的對(duì)比分析[J]. 遙感技術(shù)與應(yīng)用,33(3):465-475.

    伍顯,安裕倫,趙海兵,等,2018. 基于Landsat8的地表溫度與地形要素的關(guān)系研究[J]. 江蘇農(nóng)業(yè)科學(xué),46(18):266-270.

    熊俊楠,李偉,程維明,等,2019. 高原地區(qū)LST空間分異特征及影響因素研究:以桑珠孜區(qū)為例[J]. 國(guó)土資源遙感,31(2):164-171.

    SOBRINO J A, JIM?NEZ-MU?OZ J C, PAOLINI L, 2004. Land surface temperature retrieval from LANDSAT TM 5[J]. Remote Sensing of Environment, 90(4):434-440.

    YU Xiaolei, GUO Xulin, WU Zhaocong, 2014. Land Surface Temperature Retrieval from Landsat 8 TIRS-Comparison between Radiative Transfer Equation-based Method,Split Window Algorithm and Single Channel Method[J]. Remote sensing, 6(10):9829-9852.

    猜你喜歡
    阜平縣坡向海拔
    創(chuàng)建楓橋式人民法庭 打造多元解紛新格局——阜平縣人民法院訴源治理工作紀(jì)實(shí)
    襪子也瘋狂
    ——河北省阜平縣龍泉關(guān)學(xué)校襪子娃娃社團(tuán)活動(dòng)
    高海拔的甘南
    河北省阜平縣脫貧成效及經(jīng)驗(yàn)淺析
    DEM地表坡向變率的向量幾何計(jì)算法
    南美三大高海拔球場(chǎng)
    海拔4000米的父愛
    阜平縣栗園莊金(銀)礦化帶地質(zhì)特征及成因淺析
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    亚洲激情五月婷婷啪啪| 亚洲精品视频女| 国产精品不卡视频一区二区| 人妻少妇偷人精品九色| 亚洲av中文av极速乱| 高清在线视频一区二区三区| 女性被躁到高潮视频| 国产黄片视频在线免费观看| 久久久精品94久久精品| 中文字幕最新亚洲高清| 亚洲国产欧美在线一区| 亚洲国产精品999| 国产午夜精品一二区理论片| 久久狼人影院| av在线老鸭窝| 亚洲精品乱码久久久v下载方式| 蜜桃国产av成人99| 日韩强制内射视频| 亚洲av中文av极速乱| 国产精品一区二区三区四区免费观看| 丝袜喷水一区| 亚洲成人手机| 日本猛色少妇xxxxx猛交久久| 日日撸夜夜添| 高清午夜精品一区二区三区| 大陆偷拍与自拍| 亚洲精品一区蜜桃| 热99国产精品久久久久久7| 在线观看国产h片| 亚洲内射少妇av| 2022亚洲国产成人精品| 亚洲成色77777| 久久久久久久亚洲中文字幕| 青春草国产在线视频| 女性被躁到高潮视频| 国产av国产精品国产| 日韩制服骚丝袜av| 免费高清在线观看视频在线观看| .国产精品久久| 80岁老熟妇乱子伦牲交| 成年av动漫网址| 欧美精品亚洲一区二区| 国产成人a∨麻豆精品| 精品一区在线观看国产| 搡女人真爽免费视频火全软件| 22中文网久久字幕| 国产精品一区二区在线观看99| 考比视频在线观看| 精品亚洲乱码少妇综合久久| 国产日韩欧美视频二区| 午夜福利影视在线免费观看| 国产成人av激情在线播放 | 婷婷色麻豆天堂久久| 免费看不卡的av| 狠狠精品人妻久久久久久综合| a 毛片基地| av卡一久久| 99九九线精品视频在线观看视频| 国产熟女欧美一区二区| 午夜福利,免费看| 欧美激情 高清一区二区三区| 在线观看人妻少妇| 久久精品国产亚洲网站| 老司机影院成人| 最近中文字幕高清免费大全6| 一本色道久久久久久精品综合| 又大又黄又爽视频免费| 欧美三级亚洲精品| 在线观看免费视频网站a站| 蜜桃久久精品国产亚洲av| 午夜激情福利司机影院| 亚洲av免费高清在线观看| 观看av在线不卡| 日本色播在线视频| 国产 一区精品| 午夜福利在线观看免费完整高清在| 欧美激情国产日韩精品一区| 午夜免费男女啪啪视频观看| 亚洲人成网站在线播| 国产黄色视频一区二区在线观看| 国产精品久久久久久av不卡| 亚洲情色 制服丝袜| 多毛熟女@视频| 自拍欧美九色日韩亚洲蝌蚪91| 菩萨蛮人人尽说江南好唐韦庄| 国产精品偷伦视频观看了| 亚洲人成77777在线视频| 国产探花极品一区二区| 欧美变态另类bdsm刘玥| 亚洲美女搞黄在线观看| 老熟女久久久| 性色av一级| 黄色毛片三级朝国网站| 男女国产视频网站| 中文字幕最新亚洲高清| 欧美国产精品一级二级三级| 2018国产大陆天天弄谢| 大片电影免费在线观看免费| 九九久久精品国产亚洲av麻豆| 国产黄频视频在线观看| 晚上一个人看的免费电影| 99久久精品国产国产毛片| 观看美女的网站| 精品久久国产蜜桃| 国产精品国产三级国产av玫瑰| 精品久久久久久久久亚洲| 日本午夜av视频| 人人妻人人澡人人看| xxx大片免费视频| 亚洲成人手机| 水蜜桃什么品种好| 午夜激情福利司机影院| 女的被弄到高潮叫床怎么办| videos熟女内射| 少妇被粗大的猛进出69影院 | 少妇 在线观看| 欧美 日韩 精品 国产| av黄色大香蕉| 午夜视频国产福利| 九草在线视频观看| 高清黄色对白视频在线免费看| av在线观看视频网站免费| 草草在线视频免费看| 日本猛色少妇xxxxx猛交久久| 人人妻人人添人人爽欧美一区卜| 丰满少妇做爰视频| 少妇熟女欧美另类| 久久久久久久久久久丰满| 考比视频在线观看| 女性生殖器流出的白浆| 中文天堂在线官网| 日韩中字成人| 777米奇影视久久| 综合色丁香网| 99久久精品国产国产毛片| 久久精品久久久久久久性| 国产亚洲精品第一综合不卡 | 赤兔流量卡办理| 国产成人精品在线电影| 欧美精品人与动牲交sv欧美| 七月丁香在线播放| 一级二级三级毛片免费看| 91精品伊人久久大香线蕉| 国产片特级美女逼逼视频| 美女视频免费永久观看网站| av在线app专区| 在线播放无遮挡| 欧美日韩一区二区视频在线观看视频在线| 国产精品国产三级专区第一集| 国产一区二区在线观看av| 欧美国产精品一级二级三级| av卡一久久| 日韩一区二区三区影片| 国产亚洲av片在线观看秒播厂| 亚洲精品一二三| 夜夜看夜夜爽夜夜摸| 日本av免费视频播放| 中文精品一卡2卡3卡4更新| 精品少妇久久久久久888优播| 老司机影院成人| 在线亚洲精品国产二区图片欧美 | 最新中文字幕久久久久| 国产男人的电影天堂91| 黄色毛片三级朝国网站| 久久久久国产精品人妻一区二区| 国产日韩欧美亚洲二区| 高清视频免费观看一区二区| 在线看a的网站| 日韩精品有码人妻一区| 欧美另类一区| 午夜91福利影院| 免费看av在线观看网站| 日韩电影二区| 亚洲人成77777在线视频| 亚洲精品,欧美精品| 欧美精品国产亚洲| 高清不卡的av网站| 少妇人妻精品综合一区二区| 久久精品国产鲁丝片午夜精品| 少妇的逼水好多| 在线天堂最新版资源| 成人黄色视频免费在线看| 一级a做视频免费观看| 亚洲高清免费不卡视频| 国产亚洲一区二区精品| 日本欧美视频一区| 99国产精品免费福利视频| 制服丝袜香蕉在线| 搡女人真爽免费视频火全软件| 在线观看国产h片| 女性被躁到高潮视频| 丝瓜视频免费看黄片| 日韩成人伦理影院| 国产成人精品无人区| 国产伦精品一区二区三区视频9| 青春草亚洲视频在线观看| 中文字幕精品免费在线观看视频 | 自线自在国产av| 亚洲精品乱码久久久v下载方式| 欧美老熟妇乱子伦牲交| 男女免费视频国产| 视频区图区小说| 国产成人av激情在线播放 | 国产色爽女视频免费观看| 国产在线免费精品| 亚洲精品视频女| av天堂久久9| 韩国高清视频一区二区三区| av视频免费观看在线观看| 国产男女内射视频| 亚洲一级一片aⅴ在线观看| 国产成人免费观看mmmm| 久久久久久伊人网av| 亚洲综合色惰| 女人精品久久久久毛片| 色婷婷av一区二区三区视频| www.色视频.com| 亚洲综合色惰| 成人国产av品久久久| 亚洲国产精品专区欧美| 黑人高潮一二区| av一本久久久久| 亚洲,一卡二卡三卡| 亚洲欧美清纯卡通| 韩国av在线不卡| 亚洲图色成人| 久久精品熟女亚洲av麻豆精品| 欧美老熟妇乱子伦牲交| 国产精品一区二区三区四区免费观看| 你懂的网址亚洲精品在线观看| 91午夜精品亚洲一区二区三区| 美女主播在线视频| 亚洲精品色激情综合| 制服人妻中文乱码| 好男人视频免费观看在线| 欧美精品亚洲一区二区| 中文字幕制服av| 青春草视频在线免费观看| 国产探花极品一区二区| 亚洲第一区二区三区不卡| 免费人妻精品一区二区三区视频| 人人妻人人澡人人看| 色哟哟·www| 91精品国产国语对白视频| 久久久久精品性色| 女的被弄到高潮叫床怎么办| 黄色欧美视频在线观看| 欧美人与善性xxx| 狠狠婷婷综合久久久久久88av| 大香蕉久久成人网| 亚洲国产精品专区欧美| 日韩强制内射视频| 一级毛片我不卡| 建设人人有责人人尽责人人享有的| 97在线人人人人妻| 热99久久久久精品小说推荐| 欧美最新免费一区二区三区| 亚洲精品日韩在线中文字幕| 日本wwww免费看| 久久久久人妻精品一区果冻| 一级毛片我不卡| 热re99久久精品国产66热6| 日本91视频免费播放| 伦精品一区二区三区| 三级国产精品欧美在线观看| 国产精品一国产av| 高清av免费在线| 91精品伊人久久大香线蕉| 免费黄色在线免费观看| 精品人妻一区二区三区麻豆| 97超碰精品成人国产| av.在线天堂| 3wmmmm亚洲av在线观看| 免费久久久久久久精品成人欧美视频 | 麻豆乱淫一区二区| 欧美精品高潮呻吟av久久| 不卡视频在线观看欧美| 全区人妻精品视频| 成人手机av| 精品久久久久久久久av| 大片免费播放器 马上看| 中文字幕亚洲精品专区| av专区在线播放| 国产探花极品一区二区| a级毛片在线看网站| 成人国产av品久久久| 久久久国产精品麻豆| 最新的欧美精品一区二区| 亚洲精品色激情综合| 美女福利国产在线| 欧美3d第一页| 成人午夜精彩视频在线观看| 国产有黄有色有爽视频| 26uuu在线亚洲综合色| 伦理电影免费视频| 国产精品嫩草影院av在线观看| 亚洲美女搞黄在线观看| 欧美人与善性xxx| 日韩av在线免费看完整版不卡| 亚洲一级一片aⅴ在线观看| 国产免费视频播放在线视频| 亚洲欧美色中文字幕在线| 欧美日韩综合久久久久久| 夜夜看夜夜爽夜夜摸| 亚洲欧洲国产日韩| 校园人妻丝袜中文字幕| 一级黄片播放器| 中文字幕久久专区| 久久久久久久久久人人人人人人| 免费观看a级毛片全部| 久久久精品94久久精品| 免费看av在线观看网站| 肉色欧美久久久久久久蜜桃| 中文字幕最新亚洲高清| 国产精品久久久久久精品古装| 寂寞人妻少妇视频99o| 国产 一区精品| a级毛片免费高清观看在线播放| 插阴视频在线观看视频| 嘟嘟电影网在线观看| 另类亚洲欧美激情| 亚洲欧美一区二区三区黑人 | 妹子高潮喷水视频| 久久狼人影院| 高清午夜精品一区二区三区| 日韩中字成人| 丰满迷人的少妇在线观看| a级毛片黄视频| 精品亚洲成a人片在线观看| 免费大片黄手机在线观看| 久久久午夜欧美精品| 国产国拍精品亚洲av在线观看| 一区二区三区免费毛片| 欧美 日韩 精品 国产| 国产精品.久久久| 只有这里有精品99| 亚洲av福利一区| 日日爽夜夜爽网站| 亚洲欧美日韩另类电影网站| 免费黄频网站在线观看国产| 欧美+日韩+精品| 天天操日日干夜夜撸| 国产一区二区在线观看av| 黄片播放在线免费| 下体分泌物呈黄色| av在线app专区| 看免费成人av毛片| 亚洲综合色惰| 久久国内精品自在自线图片| 亚洲精品久久成人aⅴ小说 | 午夜久久久在线观看| 少妇人妻精品综合一区二区| 亚洲,欧美,日韩| 男女啪啪激烈高潮av片| 国产一区亚洲一区在线观看| 美女内射精品一级片tv| 国产精品久久久久久久久免| 少妇被粗大的猛进出69影院 | 如何舔出高潮| 国产成人午夜福利电影在线观看| 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡动漫免费视频| 99九九在线精品视频| 晚上一个人看的免费电影| kizo精华| 男的添女的下面高潮视频| 成人国语在线视频| 99久久综合免费| 久久久久久久大尺度免费视频| 欧美老熟妇乱子伦牲交| 男的添女的下面高潮视频| 亚洲一级一片aⅴ在线观看| 毛片一级片免费看久久久久| 热re99久久精品国产66热6| 欧美变态另类bdsm刘玥| 视频在线观看一区二区三区| 国产精品欧美亚洲77777| 在线免费观看不下载黄p国产| 亚洲精品一二三| 制服诱惑二区| 啦啦啦视频在线资源免费观看| 国产欧美日韩一区二区三区在线 | 色视频在线一区二区三区| 香蕉精品网在线| www.av在线官网国产| 少妇人妻久久综合中文| 美女cb高潮喷水在线观看| 精品一区二区三卡| 成人毛片a级毛片在线播放| 一级毛片电影观看| 性色av一级| 伦理电影大哥的女人| 99热网站在线观看| 国产精品熟女久久久久浪| 美女中出高潮动态图| 亚洲精品日韩在线中文字幕| 国产综合精华液| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 中文精品一卡2卡3卡4更新| 九色亚洲精品在线播放| 亚洲国产精品999| 日韩伦理黄色片| 国产在视频线精品| 夜夜爽夜夜爽视频| 国产黄色免费在线视频| 日韩成人av中文字幕在线观看| 久久久久久久大尺度免费视频| 黑丝袜美女国产一区| 久久人人爽人人爽人人片va| 一级黄片播放器| 久久久久国产网址| 七月丁香在线播放| 人妻少妇偷人精品九色| 国产精品不卡视频一区二区| 少妇被粗大猛烈的视频| 亚洲欧美清纯卡通| 少妇的逼水好多| 国产熟女午夜一区二区三区 | 国产精品.久久久| 黑人高潮一二区| 韩国av在线不卡| 日韩熟女老妇一区二区性免费视频| 国产精品国产三级国产av玫瑰| 黄色视频在线播放观看不卡| 久久久久久久久久成人| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 精品酒店卫生间| 国产日韩一区二区三区精品不卡 | 国产高清三级在线| 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 涩涩av久久男人的天堂| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 人妻制服诱惑在线中文字幕| 中国国产av一级| 考比视频在线观看| 国产欧美另类精品又又久久亚洲欧美| 免费黄网站久久成人精品| 欧美精品一区二区大全| 亚洲精品av麻豆狂野| 亚洲三级黄色毛片| av免费观看日本| 高清视频免费观看一区二区| 久久精品久久久久久噜噜老黄| 午夜福利影视在线免费观看| 另类亚洲欧美激情| 精品久久久久久电影网| 九色亚洲精品在线播放| av有码第一页| 免费久久久久久久精品成人欧美视频 | 人体艺术视频欧美日本| 91久久精品国产一区二区三区| 成年av动漫网址| a级毛片在线看网站| 中文字幕av电影在线播放| 久久狼人影院| 少妇丰满av| 熟女人妻精品中文字幕| 日韩熟女老妇一区二区性免费视频| 大片电影免费在线观看免费| 如何舔出高潮| 国产免费一级a男人的天堂| 精品少妇久久久久久888优播| 欧美精品高潮呻吟av久久| 18禁在线播放成人免费| 美女视频免费永久观看网站| 国产成人精品久久久久久| 色视频在线一区二区三区| 午夜影院在线不卡| 一区在线观看完整版| 男人爽女人下面视频在线观看| 欧美日本中文国产一区发布| 大片电影免费在线观看免费| 嘟嘟电影网在线观看| 久久久久久久久久久久大奶| 久久韩国三级中文字幕| 成人午夜精彩视频在线观看| 国产欧美日韩综合在线一区二区| 男人爽女人下面视频在线观看| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 乱人伦中国视频| 最新的欧美精品一区二区| 春色校园在线视频观看| 日本黄色日本黄色录像| 男女国产视频网站| 在线观看人妻少妇| 日本色播在线视频| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 精品久久久久久久久av| 国产在线一区二区三区精| 免费人妻精品一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 边亲边吃奶的免费视频| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 我的女老师完整版在线观看| 美女视频免费永久观看网站| 欧美精品人与动牲交sv欧美| 亚洲婷婷狠狠爱综合网| 国产毛片在线视频| 啦啦啦在线观看免费高清www| 狠狠精品人妻久久久久久综合| 男女边摸边吃奶| 大香蕉97超碰在线| 久久精品人人爽人人爽视色| 国产亚洲午夜精品一区二区久久| 久久精品国产鲁丝片午夜精品| 一区二区三区精品91| 日韩亚洲欧美综合| 波野结衣二区三区在线| 午夜免费鲁丝| 少妇的逼水好多| 在线观看三级黄色| 久久久久久久久久久丰满| 欧美日韩视频精品一区| 亚洲五月色婷婷综合| 两个人的视频大全免费| 欧美日韩国产mv在线观看视频| 国产在线视频一区二区| 在线观看国产h片| 久久精品熟女亚洲av麻豆精品| 桃花免费在线播放| 亚洲av不卡在线观看| 国产黄频视频在线观看| 黄色视频在线播放观看不卡| 亚洲第一区二区三区不卡| videossex国产| 国产精品一区二区三区四区免费观看| 国产精品一区二区在线不卡| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级国产av玫瑰| 免费黄网站久久成人精品| 一区在线观看完整版| 久久人人爽人人片av| 纵有疾风起免费观看全集完整版| 国语对白做爰xxxⅹ性视频网站| av有码第一页| 黄色怎么调成土黄色| 日日摸夜夜添夜夜添av毛片| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 最近最新中文字幕免费大全7| 熟妇人妻不卡中文字幕| av福利片在线| 999精品在线视频| 街头女战士在线观看网站| 免费黄色在线免费观看| 边亲边吃奶的免费视频| 2021少妇久久久久久久久久久| 成人亚洲欧美一区二区av| 啦啦啦视频在线资源免费观看| 欧美激情国产日韩精品一区| 亚洲精品国产色婷婷电影| 秋霞伦理黄片| 高清av免费在线| 一级爰片在线观看| 成人国语在线视频| 一个人免费看片子| 桃花免费在线播放| 十八禁高潮呻吟视频| 亚洲怡红院男人天堂| 亚洲人成网站在线播| 3wmmmm亚洲av在线观看| 久热久热在线精品观看| 免费观看a级毛片全部| 另类精品久久| 黄片无遮挡物在线观看| 国产69精品久久久久777片| 亚洲精品,欧美精品| 三级国产精品片| 18在线观看网站| 少妇熟女欧美另类| 18在线观看网站| 国产乱来视频区| 只有这里有精品99| 三级国产精品片| av在线观看视频网站免费| 国产精品国产三级国产专区5o| 久久久午夜欧美精品| 免费黄色在线免费观看| 日韩av在线免费看完整版不卡| 韩国高清视频一区二区三区| av免费观看日本| 亚洲av在线观看美女高潮| 中国三级夫妇交换| 日本91视频免费播放| 国产精品国产三级国产av玫瑰| 老女人水多毛片| 日韩中文字幕视频在线看片| 亚洲第一av免费看| 久久人人爽人人爽人人片va| 日韩一区二区视频免费看| 午夜久久久在线观看| 久久久精品免费免费高清| 伦理电影免费视频| 亚洲精华国产精华液的使用体验| 久久午夜综合久久蜜桃| 伦精品一区二区三区| 免费观看性生交大片5| 精品酒店卫生间| 久久午夜福利片| av.在线天堂| 亚洲国产精品成人久久小说| 日韩av不卡免费在线播放| www.av在线官网国产| 91在线精品国自产拍蜜月| 亚洲中文av在线| videossex国产| 肉色欧美久久久久久久蜜桃| 国产老妇伦熟女老妇高清|