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

    基于多源DEM數(shù)據(jù)的都江堰灌區(qū)地形因子提取及對比分析

    2024-03-01 08:30:58肖文全曹依帆馮丹禹
    四川水利 2024年1期
    關(guān)鍵詞:柵格坡度粗糙度

    肖文全,曹依帆,秦 濤,羅 尚,翟 星,馮丹禹,趙 輝

    (四川水發(fā)勘測設(shè)計(jì)研究有限公司,成都 610072)

    0 引言

    20世紀(jì)50年代,數(shù)字高程模型(DEM)的概念被提出,它被定義為一種以數(shù)字?jǐn)?shù)組形式顯示地面標(biāo)高的實(shí)體模型[1]。近年來,隨著計(jì)算機(jī)技術(shù)、地理信息技術(shù)(GIS)和遙感技術(shù)(RS)的迅速發(fā)展,DEM數(shù)據(jù)的種類逐漸增多,分辨率也逐步提高?;贒EM數(shù)據(jù)提取的地形因子可以有效地描述地貌的形態(tài)特征,從而推動(dòng)地貌形態(tài)分類的研究進(jìn)展[2]。因此,DEM數(shù)據(jù)在眾多基礎(chǔ)地理信息數(shù)據(jù)中占有極為重要的地位。針對DEM的研究也成為當(dāng)下的熱點(diǎn),高藝偉等[3]基于ASTER GDEM2數(shù)據(jù)提取陜南地區(qū)地形起伏度,并采用均值變點(diǎn)法確定最佳鄰域分析單元;韓萍等[4]基于多源DEM數(shù)據(jù),提取青島市部分地區(qū)小流域水文特征并進(jìn)行對比分析;胡最等[5]依據(jù)1∶25萬DEM數(shù)據(jù),開展湖南省地貌形態(tài)特征分類研究。雖然各專家學(xué)者已對各式DEM展開不同探討,但利用多源DEM數(shù)據(jù)進(jìn)行地形因子提取和對比分析的研究還略顯不足,且已有研究表明,不同分辨率的DEM數(shù)據(jù)所表達(dá)的地形信息的容量和精度等均存在顯著差異[6]。鑒于此,針對同一研究區(qū),利用不同分辨率DEM數(shù)據(jù)提取地形因子并進(jìn)行定量分析具有一定研究價(jià)值。本文以都江堰灌區(qū)作為研究區(qū),首先,基于三種不同分辨率的DEM數(shù)據(jù),提取高程、坡度、坡向、地形起伏度、地表粗糙度、地表切割深度等六個(gè)地形因子;接著,對提取結(jié)果進(jìn)行定性和定量分析;最后,利用皮爾遜相關(guān)性分析篩選出研究區(qū)最佳地形因子。研究結(jié)果可為都江堰灌區(qū)的地貌類型精細(xì)化劃分、災(zāi)害風(fēng)險(xiǎn)評估及土壤侵蝕等提供依據(jù)和參考。

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

    1.1 研究區(qū)概況

    都江堰灌區(qū)位于四川盆地西部,介于東經(jīng)103°29′~105°24′、北緯29°24′~31°29′之間。從地形地貌條件來看,灌區(qū)地貌類型以丘陵為主,丘陵面積占總面積的比例高達(dá)71.1%[7];海拔高度介于210~1 004 m,以龍泉山脈為分界,具有典型的西北高、東南低特征。從氣候水文方面來看,灌區(qū)屬亞熱帶濕潤氣候類型,全年氣候適宜,年均降雨量900~1 240 mm,降雨量受季節(jié)影響較大,呈現(xiàn)春旱、夏洪、秋澇、冬干的特點(diǎn)[8]。目前,都江堰灌區(qū)實(shí)際農(nóng)田灌溉面積達(dá)75.37萬hm2,涉及成都、綿陽等7市40縣(市、區(qū))。研究區(qū)概況如圖1所示。

    圖1 研究區(qū)概況

    1.2 數(shù)據(jù)來源

    本文利用的三種DEM數(shù)據(jù)分別為ASTR GDEM、SRTM DEM及ALOS DEM。其中,ALOS DEM空間分辨率為12.5 m,來源于NASA地球數(shù)據(jù)中心(https://search.asf.alaska.edu/#/);ASTR GDEM空間分辨率為30 m,SRTM DEM空間分辨率為90 m,均來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/)。基礎(chǔ)數(shù)據(jù)的準(zhǔn)確性是提取地形因子的關(guān)鍵,基于GIS軟件平臺,對三種原始DEM數(shù)據(jù)進(jìn)行鑲嵌、裁剪及投影變換等預(yù)處理,以獲取研究區(qū)同一坐標(biāo)系下準(zhǔn)確的基礎(chǔ)數(shù)據(jù)。

    2 研究方法

    2.1 歸一化

    不同的地形因子具有不同的原始值,并且對量綱和量綱單位也具有不同的要求,為了使各地形因子具有可比性,需要采取適當(dāng)?shù)氖侄螌Ω鞯匦我蜃舆M(jìn)行處理,從而將各地形因子的值限定在一定的范圍內(nèi)[9]。因此,為了便于開展地形因子間的相關(guān)性分析,本文對各地形因子的原始值進(jìn)行歸一化處理,使處理結(jié)果全部介于[0,1]。

    (1)

    式中,X為歸一化值;x為各地形因子原始值;xmin為各地形因子的最小值;xmax為各地形因子的最大值。

    2.2 皮爾遜相關(guān)性分析

    皮爾遜相關(guān)系數(shù)又叫皮爾遜積矩相關(guān)系數(shù),它是一種可以準(zhǔn)確量算兩個(gè)變量之間的相關(guān)程度的經(jīng)典統(tǒng)計(jì)學(xué)方法,被廣泛地應(yīng)用于國民消費(fèi)結(jié)構(gòu)[10]及災(zāi)害敏感性評估[11]等方面。以變量x和變量y為例,那么x與y之間的相關(guān)系數(shù)計(jì)算公式為:

    (2)

    式中,R表示變量x和變量y之間的相關(guān)系數(shù);n表示變量x和變量y的個(gè)數(shù)。相關(guān)系數(shù)與具體的相關(guān)程度如表1所示。

    表1 皮爾遜相關(guān)系數(shù)與相關(guān)程度

    3 地形因子提取與分級

    3.1 高程

    高程是地貌形態(tài)的最基本要素,DEM所提供的原始數(shù)據(jù)就是高程數(shù)據(jù)[12]。對三種高程數(shù)據(jù)按分辨率從高到低排序,獲得研究區(qū)的高程值分別介于210~1 004 m、192~1 014 m及254~1 042 m,考慮前人對高程的分級標(biāo)準(zhǔn)并結(jié)合研究區(qū)實(shí)際高程分布情況,將高程值分為<400 m、400~600 m、600~800 m、800~1 000 m及>1 000 m等五個(gè)等級,如圖2所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    3.2 坡度

    坡度定量反映了地表的陡緩程度,是常見的微觀地形因子之一[13]。利用Arc GIS 10.2軟件的表面分析工具,從三種不同DEM數(shù)據(jù)中提取坡度因子,按分辨率從高到低排序,最終提取結(jié)果為0~69.36°、0~71.28°、0~47.97°。依據(jù)曾佩楓等[1]對坡度的等級劃分,本文將坡度分為<3°、3°~5°、5°~10°、10°~15°、15°~20°、20°~25°、25°~30°及>30°等八個(gè)等級,如圖3所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    3.3 坡向

    坡向是地表任意一點(diǎn)切平面的法線矢量在水平面的投影與過該點(diǎn)正北方向的夾角,可通俗地理解為自高而低的方向[1]。與坡度的提取及排序方式相似,將研究區(qū)坡向劃分為平坦-1、北(0~22.5°及337.5°~360°)、東北(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°)等共九個(gè)類別,如圖4所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    3.4 地形起伏度

    地形起伏度(Rf)是指某一特定區(qū)域內(nèi)最高點(diǎn)高程與最低點(diǎn)高程的差值,能夠較為直觀地反映地表形態(tài)[14],是地貌類型劃分研究中最常用的因子之一,其具體計(jì)算公式為:

    Rf=Hmax-Hmin

    (3)

    最佳分析窗口的確定是影響地形起伏度準(zhǔn)確性的最關(guān)鍵因素,基于前人[2,15-16]研究成果,對于ALOSDEM、ASTR GDEM及SRTM DEM,分別采取13×13窗口,15×15窗口,17×17窗口為最佳分析窗口,其面積分別為0.026 4 km2、0.202 5 km2及2.34 km2。利用Arc GIS 10.2軟件平臺的焦點(diǎn)統(tǒng)計(jì)工具,從三種DEM數(shù)據(jù)中提取地形起伏度因子,按分辨率從大到小排序,最終結(jié)果為0~216 m、0~363 m、0~509 m。依據(jù)以往研究成果的數(shù)字格局地貌研究分類方法,將三種提取結(jié)果劃分為五類,分別為<30 m(平原)、30~50 m(臺地)、50~200 m(丘陵)、200~500 m(小起伏山地)及>500 m(中起伏度山地),如圖5所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    3.5 地表粗糙度

    地表粗糙度(R)是指某一固定區(qū)域內(nèi)地球表面積與其投影面積的比值,是反映地表形態(tài)的一個(gè)重要宏觀因子[2]。在提取的坡度數(shù)據(jù)基礎(chǔ)上,利用Arc GIS 10.2軟件的柵格計(jì)算器工具,獲取研究區(qū)地表粗糙度數(shù)據(jù)。結(jié)果分別為1~2.87、1~3.12及1~1.49。參照賈麗娜[17]的分級標(biāo)準(zhǔn),將三種分辨率的地表粗糙度數(shù)據(jù)分為1~1.03、1.03~1.08、1.08~1.14、1.14~1.21、1.21~1.29、1.29~1.39、1.39~1.57、1.57~1.93及>1.93等共九個(gè)等級,如圖6所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    3.6 地表切割深度

    地表切割深度(D)定義為某一特定區(qū)域內(nèi)平均高程與最小高程的差值[18]。在提取地表切割深度的過程中,最佳分析窗口的選擇與提取地形起伏度的窗口相同,提取結(jié)果分別為0~112.23、0~204.94及0~285.27。參照梁宏艷等[2]的劃分標(biāo)準(zhǔn),將三種提取結(jié)果分為三個(gè)等級,分別為<30 m(微切割山地)、30~100 m(淺切割山地)及100~500 m(中等切割山地),如圖7所示。

    (a)12.5 m分辨率 (b)30 m分辨率 (c)90 m分辨率

    4 結(jié)果分析

    4.1 地形因子對比分析

    為了更直觀地分析三種不同分辨率地形因子的等級劃分結(jié)果,對各地形因子不同等級的柵格數(shù)量、面積、柵格數(shù)量占比及面積占比進(jìn)行統(tǒng)計(jì)。不同分辨率下各等級高程數(shù)據(jù)的柵格數(shù)量占比如表2所示。

    表2 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(高程)

    由表2可知,對于30 m和90 m分辨率的高程因子來說,各等級所占比例接近,且400~600 m等級的占比最高,達(dá)70%以上。而對于12.5 m分辨率的高程因子,400~600 m所占比例卻不足50%,更多的區(qū)域被劃分為<400 m等級。

    不同分辨率下各等級坡度數(shù)據(jù)的柵格數(shù)量占比如表3所示。

    表3 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(坡度)

    由表3可知,三種分辨率的坡度數(shù)據(jù)在柵格數(shù)量占比上未表現(xiàn)出明顯規(guī)律。除<3°分級的差異較大外,其余各等級所占比例相差不大。

    不同分辨率下各等級坡向數(shù)據(jù)的柵格數(shù)量占比如表4所示。

    表4 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(坡向)

    由表4可知,三種分辨率的坡向數(shù)據(jù),除在平面所占比例相差較大外,其余八個(gè)分級所占比例接近,且均在10%上下波動(dòng)。

    不同分辨率下各等級地形起伏度數(shù)據(jù)的柵格數(shù)量占比如表5所示。

    表5 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(地形起伏度)

    由表5可知,30 m和90 m分辨率的地形起伏度數(shù)據(jù)在對平原和丘陵的劃分結(jié)果相似,而12.5 m和90 m分辨率的地形起伏度數(shù)據(jù)對臺地的劃分結(jié)果相似。此外,90 m分辨率的地形起伏度數(shù)據(jù)中,有極少數(shù)區(qū)域被劃分為中起伏度山地。

    不同分辨率下各等級地表粗糙度數(shù)據(jù)的柵格數(shù)量占比如表6所示。

    表6 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(地表粗糙度)

    由表6可知,三種分辨率的地表粗糙度數(shù)據(jù)中,各等級柵格數(shù)量所占比例差別不大,90%以上的柵格屬于1~1.08之間。

    不同分辨率下各等級地表切割深度數(shù)據(jù)的柵格數(shù)量占比如表7所示。

    表7 研究區(qū)不同分辨率下各等級柵格數(shù)量占比(地表切割深度)

    由表7可知,對于30 m分辨率和90 m分辨率的地表切割深度劃分結(jié)果,各級占比差異不大。其中,70%以上的區(qū)域被劃分為微切割山地,20%以上的區(qū)域被劃分為淺切割山地,僅有極少區(qū)域被劃分為中等切割山地。而對于12.5 m分辨率的地形切割深度分級結(jié)果來說,僅有微切割山地和淺切割山地兩種類型,且微切割山地的占比高達(dá)94.78%。

    4.2 最佳地形因子確定

    為了確定研究區(qū)的最佳地形因子,利用皮爾遜相關(guān)系數(shù)對各因子進(jìn)行相關(guān)性分析,具體步驟為:①利用Arc GIS 10.2軟件的隨機(jī)選點(diǎn)功能,生成2 000個(gè)獨(dú)立的隨機(jī)點(diǎn);②提取各地形因子的原始值至隨機(jī)點(diǎn);③利用公式(1),將提取的原始值進(jìn)行歸一化處理;④基于SPSS 25軟件,對歸一化過后的數(shù)據(jù)進(jìn)行皮爾遜相關(guān)分析。各因子的相關(guān)性分析結(jié)果如表8所示。

    表8 地形因子相關(guān)性分析結(jié)果

    由表8的分析結(jié)果并結(jié)合表1可知,坡向與其他地形因子間的相關(guān)系數(shù)均小于0.2,即極低相關(guān)性,表明坡向與其他地形因子相比,對研究區(qū)地形特征差異表達(dá)不明顯,應(yīng)予以剔除。除此之外,切割深度與起伏度、切割深度與坡度的相關(guān)系數(shù)均大于0.8,即極強(qiáng)相關(guān),表明三者表達(dá)的地形信息有所重復(fù)??紤]到起伏度和坡度已被大量研究確定為地形劃分的主要因子,因此,選擇保留起伏度和坡度而刪除切割深度。最終,研究區(qū)最佳地貌類型劃分的地形因子為高程、坡度、地形起伏度及地表粗糙度。

    5 結(jié)論與討論

    以都江堰灌區(qū)為研究區(qū),首先,基于三種不同分辨率的DEM數(shù)據(jù),提取高程、坡度、坡向、地形起伏度、地表粗糙度及地表切割深度等六個(gè)地形因子;接著,對提取結(jié)果進(jìn)行等級劃分和統(tǒng)計(jì);然后,利用皮爾遜相關(guān)性分析確定研究區(qū)地貌分類的最佳地形因子。研究結(jié)果如下:

    (1)由不同分辨率提取的研究區(qū)地形因子結(jié)果存在明顯差異。對高程而言,90%以上的柵格分布于<600 m等級;對坡度而言,除<3°等級相差較大外,其余等級差距較小;對坡向而言,平面等級的占比相差在10%以內(nèi),其余等級差距不大;對地形起伏度而言,12.5 m分辨率DEM的提取結(jié)果與其他兩種提取結(jié)果差異較大,所有的區(qū)域都被劃分為平原、臺地及丘陵。呈現(xiàn)出分辨率越低,提取結(jié)果種類越多的趨勢;對地表粗糙度而言,各等級柵格數(shù)量所占比例接近;對地表切割深度而言,12.5 m分辨率DEM的提取結(jié)果僅劃分出微地形切割和淺地形切割兩種,且微地形切割的比例高達(dá)95%。

    (2)由皮爾遜相關(guān)性分析結(jié)果可知,都江堰灌區(qū)的最佳地形因子為高程、坡度、地形起伏度及地表粗糙度。

    猜你喜歡
    柵格坡度粗糙度
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    基于無人機(jī)影像的巖體結(jié)構(gòu)面粗糙度獲取
    甘肅科技(2020年20期)2020-04-13 00:30:18
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    冷沖模磨削表面粗糙度的加工試驗(yàn)與應(yīng)用
    模具制造(2019年4期)2019-06-24 03:36:48
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    鋼材銹蝕率與表面三維粗糙度參數(shù)的關(guān)系
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    CT和MR對人上脛腓關(guān)節(jié)面坡度的比較研究
    不同剖面形狀的柵格壁對柵格翼氣動(dòng)特性的影響
    男人的好看免费观看在线视频| 色精品久久人妻99蜜桃| 日韩欧美免费精品| 丝袜美腿在线中文| avwww免费| 精品国产美女av久久久久小说| 亚洲熟妇中文字幕五十中出| 69人妻影院| 国产精品久久视频播放| 亚洲欧美激情综合另类| 一级黄片播放器| 两人在一起打扑克的视频| 久久久精品大字幕| h日本视频在线播放| 一进一出抽搐gif免费好疼| 国产真人三级小视频在线观看| 亚洲精品色激情综合| 亚洲av一区综合| 亚洲成av人片免费观看| 国产精品 欧美亚洲| 9191精品国产免费久久| 99久久久亚洲精品蜜臀av| 在线观看一区二区三区| 观看美女的网站| 久久午夜亚洲精品久久| 亚洲av电影不卡..在线观看| 在线播放无遮挡| h日本视频在线播放| 国内少妇人妻偷人精品xxx网站| 美女大奶头视频| 91麻豆精品激情在线观看国产| АⅤ资源中文在线天堂| 一区二区三区免费毛片| 亚洲国产中文字幕在线视频| 亚洲午夜理论影院| 97碰自拍视频| 久久精品影院6| 一个人观看的视频www高清免费观看| 别揉我奶头~嗯~啊~动态视频| 国产v大片淫在线免费观看| 啦啦啦观看免费观看视频高清| 国产精品亚洲美女久久久| 人妻夜夜爽99麻豆av| avwww免费| 亚洲性夜色夜夜综合| 国产成人av教育| 国产精品野战在线观看| 三级国产精品欧美在线观看| 免费电影在线观看免费观看| 欧美zozozo另类| 国产欧美日韩一区二区精品| 国产精品美女特级片免费视频播放器| 国产 一区 欧美 日韩| 免费av不卡在线播放| 日韩 欧美 亚洲 中文字幕| 亚洲精品国产精品久久久不卡| 黄片大片在线免费观看| 国产精品久久久久久人妻精品电影| 少妇裸体淫交视频免费看高清| 国产一区二区三区在线臀色熟女| 高清毛片免费观看视频网站| www国产在线视频色| 国产探花在线观看一区二区| 精品一区二区三区av网在线观看| 亚洲va日本ⅴa欧美va伊人久久| 一进一出抽搐gif免费好疼| 久久婷婷人人爽人人干人人爱| 成人特级av手机在线观看| 黄色视频,在线免费观看| 国产精品自产拍在线观看55亚洲| 国产高清三级在线| 国产亚洲欧美98| 人妻久久中文字幕网| 欧美大码av| 国产在线精品亚洲第一网站| 亚洲第一电影网av| 啦啦啦韩国在线观看视频| 韩国av一区二区三区四区| 国内精品久久久久精免费| netflix在线观看网站| 免费看光身美女| 波野结衣二区三区在线 | 亚洲av成人精品一区久久| 国产在视频线在精品| 在线播放国产精品三级| 日本三级黄在线观看| 免费看日本二区| 国产69精品久久久久777片| 日本在线视频免费播放| 久久精品国产亚洲av涩爱 | 国产精品98久久久久久宅男小说| 久久精品人妻少妇| 18禁国产床啪视频网站| 亚洲精品影视一区二区三区av| 97超级碰碰碰精品色视频在线观看| 国产aⅴ精品一区二区三区波| 免费无遮挡裸体视频| 无遮挡黄片免费观看| 久久久国产精品麻豆| 成年人黄色毛片网站| 久久伊人香网站| 亚洲最大成人中文| 日本五十路高清| 国产一区二区在线av高清观看| 中亚洲国语对白在线视频| 国产v大片淫在线免费观看| 18+在线观看网站| 午夜福利高清视频| 男女床上黄色一级片免费看| 日韩av在线大香蕉| 真人一进一出gif抽搐免费| 内地一区二区视频在线| 在线播放无遮挡| 亚洲 国产 在线| 波多野结衣高清作品| 亚洲精品影视一区二区三区av| 在线观看日韩欧美| 久久久精品大字幕| 中文资源天堂在线| 国产精品99久久99久久久不卡| 熟妇人妻久久中文字幕3abv| 丁香六月欧美| 国产色爽女视频免费观看| 99国产精品一区二区三区| 网址你懂的国产日韩在线| 欧美三级亚洲精品| 此物有八面人人有两片| 国产91精品成人一区二区三区| 最近视频中文字幕2019在线8| 噜噜噜噜噜久久久久久91| 免费搜索国产男女视频| 精品熟女少妇八av免费久了| 国产精品久久久久久久久免 | 五月玫瑰六月丁香| 90打野战视频偷拍视频| 欧美日韩乱码在线| 99在线视频只有这里精品首页| 久久草成人影院| 男女那种视频在线观看| 亚洲成av人片在线播放无| 小说图片视频综合网站| 午夜免费成人在线视频| 中文字幕人成人乱码亚洲影| 国产黄a三级三级三级人| 九九热线精品视视频播放| 精品久久久久久久末码| 久久久久精品国产欧美久久久| 黄色成人免费大全| 久久久成人免费电影| 久久久久久久精品吃奶| 亚洲第一欧美日韩一区二区三区| 熟妇人妻久久中文字幕3abv| 亚洲精品日韩av片在线观看 | 欧美黄色淫秽网站| 成人精品一区二区免费| 欧美三级亚洲精品| 每晚都被弄得嗷嗷叫到高潮| 国产视频一区二区在线看| 欧美在线黄色| 国产高清有码在线观看视频| 国产单亲对白刺激| 搞女人的毛片| 日韩成人在线观看一区二区三区| 美女黄网站色视频| 一边摸一边抽搐一进一小说| 国产aⅴ精品一区二区三区波| 国产视频一区二区在线看| 搡老岳熟女国产| 狠狠狠狠99中文字幕| 精品午夜福利视频在线观看一区| 国产毛片a区久久久久| 国产不卡一卡二| 18禁在线播放成人免费| 黄色片一级片一级黄色片| 一级毛片女人18水好多| 999久久久精品免费观看国产| 岛国视频午夜一区免费看| 久久九九热精品免费| 日韩欧美国产在线观看| 亚洲一区二区三区色噜噜| 一进一出好大好爽视频| 夜夜爽天天搞| 99久久成人亚洲精品观看| 国产黄片美女视频| 亚洲美女视频黄频| 国产精品久久久久久人妻精品电影| 免费在线观看亚洲国产| 91在线观看av| 久久人人精品亚洲av| 国产伦人伦偷精品视频| 国产69精品久久久久777片| 亚洲va日本ⅴa欧美va伊人久久| 啦啦啦免费观看视频1| 免费人成视频x8x8入口观看| 日韩欧美在线二视频| 欧美激情久久久久久爽电影| 久久亚洲精品不卡| 十八禁人妻一区二区| 天堂网av新在线| 精品电影一区二区在线| 国产三级黄色录像| 桃色一区二区三区在线观看| 亚洲无线观看免费| 亚洲国产精品999在线| 久久久国产精品麻豆| 欧美色欧美亚洲另类二区| 国产亚洲精品一区二区www| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| 国内久久婷婷六月综合欲色啪| 午夜福利高清视频| 嫩草影院精品99| 人妻久久中文字幕网| 两个人的视频大全免费| 久久6这里有精品| 日韩大尺度精品在线看网址| 一级a爱片免费观看的视频| 日韩亚洲欧美综合| 又黄又爽又免费观看的视频| 九色成人免费人妻av| 久久草成人影院| 亚洲成人久久性| 丝袜美腿在线中文| 成人精品一区二区免费| www.色视频.com| 非洲黑人性xxxx精品又粗又长| 精品福利观看| 色综合站精品国产| 国产精品永久免费网站| 神马国产精品三级电影在线观看| 成人无遮挡网站| 久久精品国产综合久久久| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 国产视频一区二区在线看| 制服丝袜大香蕉在线| 精品一区二区三区av网在线观看| 亚洲成人中文字幕在线播放| 男人的好看免费观看在线视频| 欧美绝顶高潮抽搐喷水| 国产成人福利小说| 1024手机看黄色片| 亚洲国产欧美网| 久久99热这里只有精品18| 亚洲精品影视一区二区三区av| 国产色婷婷99| 精品国产三级普通话版| 三级毛片av免费| 日韩中文字幕欧美一区二区| 日日干狠狠操夜夜爽| 午夜精品一区二区三区免费看| 婷婷精品国产亚洲av在线| 国产精品免费一区二区三区在线| 日本成人三级电影网站| 亚洲熟妇中文字幕五十中出| 精品无人区乱码1区二区| 欧美日韩精品网址| 宅男免费午夜| 无人区码免费观看不卡| 日韩精品中文字幕看吧| 99热只有精品国产| 熟女人妻精品中文字幕| 亚洲国产欧美网| 日本与韩国留学比较| 母亲3免费完整高清在线观看| 人人妻人人看人人澡| 久久精品国产清高在天天线| 99久国产av精品| 日韩精品青青久久久久久| 十八禁人妻一区二区| 99久久九九国产精品国产免费| 丰满乱子伦码专区| 又黄又爽又刺激的免费视频.| 免费观看a级毛片全部| 最近的中文字幕免费完整| 欧美三级亚洲精品| 男女视频在线观看网站免费| 欧美人与善性xxx| 2018国产大陆天天弄谢| 国产久久久一区二区三区| 只有这里有精品99| 在线免费十八禁| 色网站视频免费| 两个人视频免费观看高清| 国产午夜精品一二区理论片| 亚洲18禁久久av| 内地一区二区视频在线| 国产精品无大码| 亚洲成人久久爱视频| 精品熟女少妇av免费看| 日日摸夜夜添夜夜爱| 国产伦理片在线播放av一区| 国产国拍精品亚洲av在线观看| av黄色大香蕉| 婷婷色综合大香蕉| 国产亚洲一区二区精品| 久久久亚洲精品成人影院| 18禁裸乳无遮挡免费网站照片| 亚洲在久久综合| 噜噜噜噜噜久久久久久91| 免费人成在线观看视频色| 日韩大片免费观看网站| 蜜桃久久精品国产亚洲av| 国产午夜精品久久久久久一区二区三区| 午夜爱爱视频在线播放| 国产午夜福利久久久久久| 青春草亚洲视频在线观看| 狂野欧美白嫩少妇大欣赏| 一二三四中文在线观看免费高清| 最近2019中文字幕mv第一页| 午夜免费男女啪啪视频观看| 亚洲精品日韩av片在线观看| 欧美xxxx黑人xx丫x性爽| 国产在线男女| 精品酒店卫生间| av在线蜜桃| 男女国产视频网站| 亚洲高清免费不卡视频| 人妻少妇偷人精品九色| 久久久久国产网址| 在线免费观看的www视频| 美女高潮的动态| 国产精品.久久久| 亚洲自偷自拍三级| 中文字幕免费在线视频6| 久久热精品热| 性色avwww在线观看| 人妻一区二区av| 一级a做视频免费观看| 色尼玛亚洲综合影院| 日韩精品青青久久久久久| 午夜免费激情av| 成人午夜精彩视频在线观看| 亚洲av男天堂| 久热久热在线精品观看| 永久免费av网站大全| 大片免费播放器 马上看| 看黄色毛片网站| 日韩精品有码人妻一区| 毛片女人毛片| 国产探花极品一区二区| 亚洲18禁久久av| 亚洲欧美一区二区三区国产| www.av在线官网国产| 777米奇影视久久| 干丝袜人妻中文字幕| 国产av码专区亚洲av| 91在线精品国自产拍蜜月| 国产午夜精品一二区理论片| 黑人高潮一二区| 国产精品蜜桃在线观看| 亚洲精品自拍成人| 国产在视频线精品| 搞女人的毛片| 免费av毛片视频| 欧美 日韩 精品 国产| 亚洲精品视频女| 日韩欧美精品免费久久| 特级一级黄色大片| 中文精品一卡2卡3卡4更新| 欧美成人午夜免费资源| 日韩精品青青久久久久久| 国产在视频线精品| 成人高潮视频无遮挡免费网站| 日韩强制内射视频| 国产色婷婷99| 亚洲一级一片aⅴ在线观看| 欧美日韩亚洲高清精品| 男的添女的下面高潮视频| 欧美性感艳星| 尤物成人国产欧美一区二区三区| 国产精品福利在线免费观看| 2021少妇久久久久久久久久久| 久久久久久久久久成人| 91av网一区二区| 亚洲国产精品成人综合色| 亚洲av中文字字幕乱码综合| 国产黄a三级三级三级人| 亚洲精品一区蜜桃| 七月丁香在线播放| 天堂影院成人在线观看| 亚洲一级一片aⅴ在线观看| 男女那种视频在线观看| xxx大片免费视频| 青春草视频在线免费观看| 午夜激情久久久久久久| 精品久久久噜噜| 亚洲一级一片aⅴ在线观看| 亚洲天堂国产精品一区在线| 国产成人免费观看mmmm| 十八禁国产超污无遮挡网站| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品50| 国产av不卡久久| 成人欧美大片| 狂野欧美激情性xxxx在线观看| 亚洲国产日韩欧美精品在线观看| 欧美日韩视频高清一区二区三区二| 亚洲在久久综合| 欧美激情在线99| av在线亚洲专区| 搡老妇女老女人老熟妇| 免费看光身美女| 丝袜美腿在线中文| 大陆偷拍与自拍| 国产亚洲av嫩草精品影院| av在线蜜桃| 日韩 亚洲 欧美在线| 国产极品天堂在线| 免费看日本二区| 免费播放大片免费观看视频在线观看| 亚洲国产精品国产精品| 我的女老师完整版在线观看| 午夜精品一区二区三区免费看| 国产精品人妻久久久影院| 成人一区二区视频在线观看| 亚洲国产欧美人成| 国产成人午夜福利电影在线观看| 国产高清国产精品国产三级 | 麻豆av噜噜一区二区三区| 免费观看a级毛片全部| 成年av动漫网址| 亚洲综合色惰| 夜夜爽夜夜爽视频| 激情五月婷婷亚洲| 欧美人与善性xxx| 一二三四中文在线观看免费高清| 国产男女超爽视频在线观看| 亚洲va在线va天堂va国产| 国产精品一区二区性色av| 最近的中文字幕免费完整| 国产精品福利在线免费观看| 欧美极品一区二区三区四区| 一级毛片 在线播放| 大又大粗又爽又黄少妇毛片口| 麻豆国产97在线/欧美| 亚洲欧美成人精品一区二区| 国产av国产精品国产| 国产男女超爽视频在线观看| 婷婷色综合www| 我的老师免费观看完整版| 免费黄色在线免费观看| 日本色播在线视频| 国产在线一区二区三区精| 最近最新中文字幕免费大全7| 大又大粗又爽又黄少妇毛片口| 久久久精品94久久精品| 国产欧美日韩精品一区二区| 久久久精品欧美日韩精品| 欧美人与善性xxx| 天天躁日日操中文字幕| 夜夜看夜夜爽夜夜摸| 欧美97在线视频| 菩萨蛮人人尽说江南好唐韦庄| 99热这里只有是精品50| 最近视频中文字幕2019在线8| 黄色配什么色好看| 国产在视频线精品| 亚洲精品乱码久久久久久按摩| 天堂俺去俺来也www色官网 | 成人二区视频| 日日撸夜夜添| 日韩一区二区视频免费看| 亚洲av一区综合| 免费av不卡在线播放| 精品欧美国产一区二区三| 亚洲精品乱码久久久久久按摩| 国产一级毛片在线| 亚洲综合色惰| 亚洲人成网站在线观看播放| 69人妻影院| 国产精品久久久久久久电影| 精品人妻熟女av久视频| eeuss影院久久| 午夜久久久久精精品| 久久久久久国产a免费观看| 精品少妇黑人巨大在线播放| 国产在视频线精品| 国产极品天堂在线| 99久久人妻综合| 国产亚洲91精品色在线| 亚洲无线观看免费| 最后的刺客免费高清国语| 日韩人妻高清精品专区| 如何舔出高潮| 亚洲熟妇中文字幕五十中出| 特级一级黄色大片| 国产在视频线精品| 亚洲真实伦在线观看| 久久久久久久国产电影| 一级片'在线观看视频| 久久久精品94久久精品| 99久久九九国产精品国产免费| 99久国产av精品国产电影| 神马国产精品三级电影在线观看| 亚洲图色成人| 久久久精品94久久精品| h日本视频在线播放| 久久久久精品性色| 国产一区二区在线观看日韩| 久久久久久久亚洲中文字幕| 偷拍熟女少妇极品色| 精品不卡国产一区二区三区| 99热全是精品| 亚洲综合色惰| 91精品一卡2卡3卡4卡| 18禁在线无遮挡免费观看视频| ponron亚洲| 亚洲av不卡在线观看| 国产男人的电影天堂91| 十八禁网站网址无遮挡 | 91aial.com中文字幕在线观看| 极品少妇高潮喷水抽搐| 嫩草影院精品99| .国产精品久久| 国产麻豆成人av免费视频| 成人国产麻豆网| 午夜免费观看性视频| 欧美zozozo另类| 永久免费av网站大全| 极品教师在线视频| 天堂√8在线中文| 十八禁网站网址无遮挡 | 亚洲精品日韩av片在线观看| 精品久久久噜噜| 日本爱情动作片www.在线观看| 天堂中文最新版在线下载 | av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 全区人妻精品视频| 国产精品国产三级国产专区5o| 亚洲av成人精品一二三区| 日韩欧美一区视频在线观看 | 国内少妇人妻偷人精品xxx网站| 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 国产成人freesex在线| 老司机影院成人| 日本午夜av视频| 黄色欧美视频在线观看| 搡老妇女老女人老熟妇| 永久免费av网站大全| 午夜精品国产一区二区电影 | 日韩av在线大香蕉| 日韩av不卡免费在线播放| 久久精品国产亚洲av天美| 综合色av麻豆| 精品不卡国产一区二区三区| 成年免费大片在线观看| 99久久中文字幕三级久久日本| 久久久久免费精品人妻一区二区| 午夜免费激情av| 亚洲欧美日韩无卡精品| 欧美97在线视频| 亚洲aⅴ乱码一区二区在线播放| 永久免费av网站大全| 精品久久久久久成人av| 久久久久久国产a免费观看| 亚洲成色77777| 久久久久性生活片| 亚洲精品色激情综合| 不卡视频在线观看欧美| 国产人妻一区二区三区在| 亚洲美女视频黄频| 91精品一卡2卡3卡4卡| 春色校园在线视频观看| av免费观看日本| 免费观看精品视频网站| 丝袜美腿在线中文| 亚洲精品日韩在线中文字幕| 老女人水多毛片| 欧美一区二区亚洲| 我的老师免费观看完整版| 久久久久久久大尺度免费视频| 男女边吃奶边做爰视频| 精品不卡国产一区二区三区| 女人被狂操c到高潮| 高清av免费在线| 看十八女毛片水多多多| 久久人人爽人人片av| 亚洲av男天堂| 亚洲性久久影院| 国内精品宾馆在线| 少妇被粗大猛烈的视频| 午夜免费男女啪啪视频观看| 午夜福利视频1000在线观看| 91aial.com中文字幕在线观看| 午夜久久久久精精品| 18禁在线播放成人免费| 永久免费av网站大全| 别揉我奶头 嗯啊视频| 国产精品精品国产色婷婷| 亚洲av日韩在线播放| 成人亚洲精品av一区二区| 国产亚洲最大av| 夜夜看夜夜爽夜夜摸| 中文字幕亚洲精品专区| 噜噜噜噜噜久久久久久91| 欧美三级亚洲精品| 久久久久久久大尺度免费视频| 综合色av麻豆| 免费观看av网站的网址| 99久国产av精品| 国产亚洲最大av| 国产成人精品久久久久久| 国产v大片淫在线免费观看| a级毛片免费高清观看在线播放| 99视频精品全部免费 在线| 女的被弄到高潮叫床怎么办| 淫秽高清视频在线观看| 亚洲国产高清在线一区二区三| 亚洲精品国产av蜜桃| 日韩强制内射视频| 日本-黄色视频高清免费观看| 一级av片app|