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

    中國西南地區(qū)地形起伏度的最佳分析尺度確定

    2018-03-29 03:07:50靜,
    水土保持通報 2018年1期
    關(guān)鍵詞:變點(diǎn)度值西南地區(qū)

    鐘 靜, 盧 濤

    (1.中國科學(xué)院 成都生物研究所, 中國科學(xué)院 山地生態(tài)恢復(fù)與生物資源利用重點(diǎn)實驗室,生態(tài)恢復(fù)與生物多樣性保育四川省重點(diǎn)實驗室, 四川 成都 610041; 2.中國科學(xué)院大學(xué), 北京 100049)

    地形起伏度(relief amplitude)指某一確定面積內(nèi)最高點(diǎn)和最低點(diǎn)的高差[1-2]。其是受地質(zhì)構(gòu)造與地表剝蝕相互作用而成[3-4],是定量描述地貌形態(tài)特征,劃分地貌單元的重要依據(jù)[5-6]。在傳統(tǒng)研究方法中,由于數(shù)據(jù)量大和計算繁瑣方面的限制,地形起伏度僅在地貌地形圖的繪制等方面發(fā)揮作用[7]。隨著區(qū)域人口、經(jīng)濟(jì)的發(fā)展,生態(tài)環(huán)境問題逐漸引起人們的關(guān)注,地形起伏度作為生態(tài)環(huán)境評價的重要指標(biāo),日益受到國內(nèi)外學(xué)者的廣泛應(yīng)用。近年來,隨著數(shù)字高程模型(digital elevation model, DEM)數(shù)據(jù)庫的建立和計算機(jī)的廣泛使用,以DEM為基本信息源對區(qū)域地形起伏度的研究逐漸興盛起來[8]。

    提取地形起伏度的關(guān)鍵在于確定最佳統(tǒng)計單元[9],即確定一定的面積,使得該面積既能反映區(qū)域地貌的完整性又能普遍代表該區(qū)域的地貌特征[1]。最佳統(tǒng)計單元的確定具有較大的不確定性和明顯的尺度效應(yīng)[10-13],隨著統(tǒng)計單元大小的改變,區(qū)域內(nèi)高差值隨之改變,使得單元內(nèi)地形起伏度值也發(fā)生改變,進(jìn)而直接影響整個研究區(qū)域地形起伏度的提取精度[14]。此外,由于不同研究區(qū)地形特征各異,各個領(lǐng)域研究目的不盡相同,對數(shù)據(jù)精度要求也不同。因此,研究所采用的基礎(chǔ)數(shù)據(jù)、提取方法、最佳統(tǒng)計單元的確定以及等級劃分都存在差異[15]。關(guān)于最佳統(tǒng)計單元的確定,主要方法有樣地判別法、最大高差法、人工目測法和均值變點(diǎn)分析法等。如涂漢明等[5]采用樣地判別法,認(rèn)為中國存在2,10,16,20,22 km2共5種不同規(guī)模的地形起伏度最佳統(tǒng)計單元。王玲等[16]采用新疆1∶25萬DEM 數(shù)據(jù),基于3×3,5×5,7×7,…,30×30的窗口序列,通過最大高差法和人工目測法確定新疆地形起伏度的最佳統(tǒng)計單元為2.56 km2。王讓虎等[17]采用均值變點(diǎn)分析法,確定中國東北地區(qū)地形起伏度提取的最佳窗口為2.62 km2。陳學(xué)兄等[18]同樣采用窗口分析法和均值變點(diǎn)分析法,確定0.899 km2為陜西省地形起伏度的最佳統(tǒng)計單元??梢?,窗口分析法和均值變點(diǎn)分析法已成為確定地形起伏度最佳統(tǒng)計單元的主流研究方法。

    目前,中國地形起伏度的研究區(qū)域多集中于北部及中東部地區(qū),而對地形起伏較大,結(jié)構(gòu)復(fù)雜的區(qū)域,研究相對較少[19]。中國西南地區(qū)作為地貌類型多樣的獨(dú)特地理單元,其復(fù)雜的地形特征極大地影響著區(qū)域的生態(tài)格局。明確西南地區(qū)地形結(jié)構(gòu)特點(diǎn)和地貌結(jié)構(gòu)劃分,對區(qū)域環(huán)境保護(hù)和發(fā)展具有重要意義。本文擬在總結(jié)概括國內(nèi)外關(guān)于地形起伏度的研究方法和技術(shù)思路的前提下,確定以中國西南地區(qū)分辨率為30 m×30 m的DEM數(shù)據(jù)為基礎(chǔ)數(shù)據(jù),通過地貌形態(tài)特征劃分地貌類型中典型區(qū)域,并利用Python模塊實現(xiàn)窗口的提取以及最佳統(tǒng)計單元的確定,從而實現(xiàn)西南地區(qū)地形起伏度的提取,其中如何兼顧各地貌類型的特點(diǎn)并確定出唯一適宜整個西南地區(qū)地形起伏度提取的最佳統(tǒng)計窗口成為本文研究重點(diǎn)。

    1 研究區(qū)概況

    本文所指的中國西南地區(qū),包括云南、四川、貴州、廣西、重慶、青海的南部地區(qū)以及西藏的中東部地區(qū),地理坐標(biāo)為東經(jīng)83.862°—112.069°,北緯21.135°—36.482°,總面積為2.55×106km2。西南地區(qū)橫跨中國地形一、二、三級階梯,地形結(jié)構(gòu)復(fù)雜,地貌類型多樣且內(nèi)部空間差異較大[19],既有平均海拔高達(dá)4 500 m的青藏高原,又有地勢平坦的成都平原,還有山高谷深、山河相間的橫斷山區(qū),以及石灰?guī)r廣布,地表崎嶇不平的喀斯特地貌區(qū)。西南地區(qū)生態(tài)系統(tǒng)復(fù)雜多樣,生物多樣性和水資源等自然資源豐富,是中國長江流域、珠江流域以及西南諸河的源頭及上游區(qū),對于維系中國華東、華南乃至東南亞的生態(tài)及社會經(jīng)濟(jì)安全起著重要的作用,是中國重要的生態(tài)屏障區(qū)。

    2 研究材料與方法

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

    文中采用的DEM數(shù)據(jù)為ASTER GDEMv2(30 m),數(shù)據(jù)來源于中國科學(xué)院計算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.Gscloud.cn)。原始數(shù)據(jù)覆蓋研究區(qū)范圍,通過ArcGIS軟件進(jìn)行拼接投影等前期處理后,采用西南地區(qū)的矢量圖形對DEM數(shù)據(jù)進(jìn)行裁剪,從而得到整個西南地區(qū)的DEM數(shù)據(jù)。

    2.2 地貌類型劃分

    西南地區(qū)地形復(fù)雜,涵蓋高原、山地、丘陵、平原等多種地貌類型。為了更為全面地掌握整個西南地區(qū)地形起伏度狀況,以及各類地貌類型間最佳統(tǒng)計面積的差異,本文分別針對各代表性地貌類型以及整個西南地區(qū),確定了各自地形起伏度所對應(yīng)的最佳統(tǒng)計單元。由于地貌發(fā)育的階段性和特殊性,自然區(qū)劃往往無法給出各種地貌劃分的明確界線。因此,本文在參考已有文獻(xiàn)和中國地形圖的基礎(chǔ)上,確定出高原、山地、丘陵以及平原作為研究區(qū)4種代表性地貌單元,以此兼顧區(qū)域地貌完整性與普適性。

    2.3 地形起伏度提取

    采用窗口分析方法提取各典型地貌區(qū)的地形起伏度值。如圖1所示,窗口分析法的基本原理是:基于DEM柵格數(shù)據(jù),確定具有固定分析半徑n×n的窗口值,通過設(shè)定起始窗口,移動步距以及終止窗口分別計算各窗口內(nèi)最大值和最小值的差值作為目標(biāo)窗口的地形起伏度值。最后,通過窗口步距進(jìn)行逐一移動,從而實現(xiàn)各個窗口下目標(biāo)柵格的地形起伏度值獲取。

    圖1 窗口分析示意圖

    本文根據(jù)區(qū)域地貌類型,確定了2×2(0.036 km2)為起始窗口,移動步距為1,終止窗口為183×183(30.14 km2)的矩形分析窗口,分別對青藏高原、橫斷山區(qū)、喀斯特地貌和成都平原4種典型地貌進(jìn)行地形起伏度提取。依托ArcGIS軟件的空間分析功能,通過Spatial Analyst工具中的柵格鄰域計算工具Neighborhood Statistics,選擇矩形分析窗口并依次計算窗口n×n(n=2,3,4,…,181,182,183)內(nèi)最大高程值和最小高程值,進(jìn)而利用Raster Calculator計算最大值層和最小值層的差值即可得到地勢起伏度層。其中,計算所得每個柵格的值表示以該窗口為目標(biāo)窗口的確定鄰域內(nèi)的地形起伏度值。最后,計算n×n窗口下柵格起伏度值的平均值,所得值作為該分析窗口下的地形起伏度值。由于本文研究區(qū)域范圍較廣,采用30 m分辨率的DEM數(shù)據(jù),并且窗口設(shè)置數(shù)量較多,窗口移動步距較短,使得計算過程繁復(fù),數(shù)據(jù)處理緩慢并產(chǎn)生大量緩存文件。因此,本文通過ArcGIS中的Python模塊,通過編譯程序,實現(xiàn)地形起伏度的自動提取以減少數(shù)據(jù)誤差,提高數(shù)據(jù)精度。

    2.4 最佳統(tǒng)計單元確定

    地形起伏度提取的關(guān)鍵在于最佳統(tǒng)計單元的確定。已有研究表明,地勢起伏度隨統(tǒng)計面積的變化曲線呈logarithmic曲線[20-21]。即隨著統(tǒng)計單元面積的擴(kuò)大,地形起伏度值隨之呈現(xiàn)出先快速增大,經(jīng)過唯一拐點(diǎn)之后增速明顯減緩并呈現(xiàn)緩慢增大的過程。該曲線上存在的唯一“拐點(diǎn)”所對應(yīng)的統(tǒng)計單元面積即為最佳統(tǒng)計單元。均值變點(diǎn)分析法對有且只有一個拐點(diǎn)的檢驗最為有效[16],該方法已逐漸發(fā)展為類似研究的最主流方法。本文以統(tǒng)計窗口面積為橫坐標(biāo),各窗口對應(yīng)的平均起伏度值為縱坐標(biāo)繪制窗口面積與地形起伏度對應(yīng)關(guān)系的擬合曲線,通過均值變點(diǎn)分析法尋找拐點(diǎn)?;居嬎氵^程為:根據(jù)研究數(shù)據(jù)設(shè)定以漸變窗口下的平均地形起伏度值作為非線性系統(tǒng)輸出數(shù)據(jù),即樣本序列{xm,m=1,2,…,N}。其中,每一個m將樣本量分為2段,分別是m1,m2,m3,…,mn-1與mn,mn+1,…,mN。該方法主要由以下3個步驟完成:

    (1) 計算各段樣本的算數(shù)平均值(Xm1)和(Xm2)以及各段的統(tǒng)計量Sm。

    (1)

    式中:Sm——2段樣本離差平方和的差值;N——總樣本數(shù)。

    (2) 計算總樣本量的離差平方和。

    (2)

    (3) 計算期望值E(S-Sm)。變點(diǎn)的存在會使總樣本的統(tǒng)計量S與樣本分段后的統(tǒng)計量Sm之間的差距增大[20],最大期望值所對應(yīng)的窗口大小即為最佳統(tǒng)計窗口。

    根據(jù)均值變點(diǎn)分析法的計算原理,同樣通過編譯Python程序運(yùn)行計算過程,從而確定各典型地貌單元和整個西南地區(qū)地形起伏度提取的最佳統(tǒng)計面積。

    2.5 地形起伏度空間自相關(guān)分析

    空間自相關(guān)是指一定范圍研究區(qū)域內(nèi),研究對象與空間位置之間存在的相關(guān)關(guān)系。該值是檢驗?zāi)骋豢臻g領(lǐng)域的屬性值與相鄰空間上的屬性值是否相關(guān)的重要指標(biāo)[22]??臻g自相關(guān)分為正、負(fù)兩類,正相關(guān)表示某單元的屬性值與其相鄰空間單元的屬性值具有相同的變化趨勢,表示空間格局存在集聚性,負(fù)相關(guān)則相反[23]。

    運(yùn)用該模型將各個坡面采樣點(diǎn)的 7Be含量轉(zhuǎn)化為表層土壤的侵蝕或沉積量可以說明坡面土壤侵蝕的空間分布特征。

    本文以全局自相關(guān)指數(shù)Moran’s I對屬性值在整個區(qū)域的空間特征進(jìn)行描述,其值介于-1~1之間[23-24]。所得值大于0,表明研究對象存在空間正相關(guān);等于零表示不存在空間相關(guān)性;而小于零則表示存在負(fù)相關(guān)關(guān)系[25-26]。計算公式如下:

    (3)

    由于全局空間自相關(guān)不能對聚集或異常發(fā)生的空間位置進(jìn)行確定。因此本研究進(jìn)一步采用空間關(guān)聯(lián)局域指標(biāo)LISA(local indicators of spatial association)指數(shù)分析局部空間自相關(guān)性。其本質(zhì)是將Moran’s指數(shù)分解到一定的區(qū)域單元,以分析要素屬性與周邊屬性值的相關(guān)關(guān)系[27]。其關(guān)聯(lián)模式可分為“高—高”關(guān)聯(lián)、“低—低”關(guān)聯(lián)2種空間正相關(guān),以及“低—高”、“高—低”關(guān)聯(lián)2種空間負(fù)相關(guān)。其中“高—高”表示某一區(qū)域的屬性值與其周圍區(qū)域的屬性值都較高;“高—低”表示某一區(qū)域的屬性值較高而周圍區(qū)域?qū)傩灾递^低;“低—低”與“低—高”則反之。將整個西南地區(qū)縣級行政單元與地形起伏度值進(jìn)行疊加,采用ArcGIS的zone statistic工具進(jìn)行各縣域地形起伏度的分區(qū)統(tǒng)計并完成基于縣域的地形起伏度賦值,進(jìn)而通過空間分析工具實現(xiàn)西南地區(qū)地形起伏度的空間自相分析。

    3 結(jié)果與分析

    3.1 地形起伏度擬合曲線

    通過Python模塊編程,實現(xiàn)高原(青藏高原區(qū))、山地(橫斷山區(qū))、丘陵(喀斯特地貌區(qū))、平原(成都平原區(qū))以及整個西南地區(qū)2×2,3×3,…,182×182,183×183窗口下的地形起伏度提取。得到窗口面積及對應(yīng)的平均地形起伏度值如表1所示。

    表1 研究區(qū)網(wǎng)格單元大小與平均地形起伏度的對應(yīng)關(guān)系

    從表1可見,4個典型地貌類型區(qū)以及整個西南地區(qū)地形起伏度值均呈現(xiàn)隨窗口對應(yīng)的柵格面積增加而增加的特征。4個地貌區(qū)在同樣的統(tǒng)計面積下計算出的地形起伏度值差異較大。如同為50×50的窗口下,高原(青藏高原區(qū))計算所得的地形起伏度值為120.83 m,山區(qū)(橫斷山區(qū))為555.02 m,丘陵(喀斯特地貌區(qū))地形起伏度值為257.11 m,平原(成都平原區(qū))為47.27 m,不同地貌類型地形起伏度值差值高達(dá)507.75 m。而整個西南地區(qū)由于面積遠(yuǎn)遠(yuǎn)大于其他區(qū)域,在同樣的50×50窗口下所得的地形起伏度值高達(dá)349.99 m。可見,由于地貌類型以及研究區(qū)面積大小的不同,最佳統(tǒng)計尺度也存在明顯差異。

    利用Origin軟件對表1中各窗口對應(yīng)的統(tǒng)計面積和平均地形起伏度進(jìn)行對數(shù)方程擬合,分別得到青藏高原區(qū)、橫斷山區(qū)、喀斯特地貌區(qū)和成都平原區(qū)4種地貌類型以及整個西南地區(qū)的地形起伏度擬合曲線(圖2)。從圖2中可以看出,所有區(qū)域的擬合曲線均擬合效果較好,且通過統(tǒng)計檢驗。

    圖2 統(tǒng)計單元與平均地形起伏度對應(yīng)關(guān)系擬合曲線

    3.2 最佳統(tǒng)計單元確定

    由圖2可見,4個地貌區(qū)的地形起伏度值在0~2 km2均呈現(xiàn)出快速增加的趨勢,曲線較陡;在2~3 km2起伏度值增速減緩,曲線隨之由陡變緩。說明在2~3 km2存在曲線由陡變緩的拐點(diǎn)?;诰底凕c(diǎn)分析法,采用Python模塊編程分別計算4種地貌類型的典型區(qū)域以及整個西南地區(qū)的最佳統(tǒng)計窗口值,得出青藏高原最佳統(tǒng)計窗口為50×50(2.25 km2);橫斷山區(qū)最佳統(tǒng)計窗口為54×54(2.62 km2);喀斯特地貌區(qū)最佳統(tǒng)計單元為51×51(2.34 km2);成都平原最佳統(tǒng)計窗口為49×49(2.16 km2);整個西南地區(qū)最佳統(tǒng)計窗口為52×52(2.43 km2)。其中橫斷山區(qū)與成都平原之間的窗口面積相差最大(0.46 km2),其他區(qū)域存在差異但相對較小。

    3.3 西南地區(qū)地形起伏度

    通過計算52×52(2.43 km2)窗口下西南地區(qū)地形起伏度值,得到整個西南地區(qū)地形起伏度的專題圖層。在國內(nèi)外地貌傳統(tǒng)分類的基礎(chǔ)上,采用地形起伏度作為第一分級指標(biāo),將中國地形起伏度以30,70,200,500,1 000,2 500 m為間隔,劃分為7個等級。本研究根據(jù)這7個等級將地形起伏度分為平坦、緩起伏、微起伏、小起伏、中起伏、大起伏和極大起伏區(qū)域,進(jìn)而繪制西南地區(qū)地形起伏度圖(圖3),并計算得到各級起伏度的面積(表2)。

    圖3 中國西南地區(qū)地形起伏度分級

    根據(jù)圖3和表2可知,整個西南地區(qū)地形起伏度主要以中、小起伏為主,低起伏次之。其中,小起伏度面積最大,占總面積38.68%;中起伏度面積占比23.58%,二者占總面積的比例之和高達(dá)62.26%;低起伏區(qū)面積比為21.15%,僅次于中起伏度地區(qū);緩起伏區(qū)和平坦區(qū)面積占總研究區(qū)面積的比例分別為8.15%和5.80%;只有2 900.85 km2的區(qū)域為極大起伏區(qū),僅占總面積的0.11%??梢?,整個西南地區(qū)地形起伏度差異較大,除了極少部分區(qū)域存在極大起伏和平坦區(qū)外,大部分地區(qū)地形起伏度呈現(xiàn)出中、小起伏的特征。

    從分布格局而言,西南地區(qū)地形起伏度的分布規(guī)律為:平坦區(qū)(0~30 m)和緩起伏區(qū)域(30~70 m)分布較為集中,主要分布在西北部區(qū)域和東北部區(qū)域;微起伏區(qū)(70~200 m)主要分布在東南部部分區(qū)域;大起伏度區(qū)(1 000~2 500 m)以及極大起伏區(qū)(2 500~6 903 m)則主要分布在研究區(qū)域的中部地區(qū);而小起伏(200~500 m)和中起伏區(qū)(500~1 000 m)主要分布在東南部區(qū)域。從整體地勢起伏情況來看,西南地區(qū)呈現(xiàn)出中部高,東南部次之,西北部和東北部較低的特征。這與西南地區(qū)地貌類型分布情況相適應(yīng):研究區(qū)西北部為青藏高原區(qū),東緣為中起伏和大起伏區(qū)域交錯分布,可見高原邊緣地形起伏較大,而高原內(nèi)部地形起伏度較小,多為30~70 m的緩起伏與高原邊緣起伏較大而高原面較為平坦的特點(diǎn)相符合;中部橫斷山區(qū)橫跨中國地形第一、二階梯,山河相間,縱列分布,山高谷深,地形起伏度大,成為整個研究區(qū)內(nèi)大、中起伏和極大起伏的主要分布區(qū);中東部地區(qū)為云貴高原,由于該區(qū)域內(nèi)部廣泛分布著崎嶇不平的喀斯特地貌,使得該區(qū)地形起伏主要以中小起伏為主;東南部涵蓋兩廣丘陵區(qū)域,而東北部則為四川盆地以及盆地內(nèi)部的成都平原,與圖示低起伏、微起伏區(qū)和平坦區(qū)分布區(qū)域相符合。

    表2 中國西南地區(qū)地形起伏度分級

    3.4 地形起伏度空間自相關(guān)分析

    全局空間相關(guān)性分析結(jié)果表明,整個西南地區(qū)地形起伏度的Moran’I值為0.603(p<0.001)??梢娢髂系貐^(qū)的地形起伏度分布在整體上存在著集聚或者分散格局。為進(jìn)一步了解西南地區(qū)地形起伏度的空間異質(zhì)性并明確集聚區(qū)域所在的空間位置,本文繪制了中國西南地區(qū)地形起伏度的LISA集聚圖(圖4)。從圖4可以看出,呈現(xiàn)空間正相關(guān)性的“高—高”區(qū)域主要集中在青藏高原南緣、橫斷山區(qū)、川西北高山高原區(qū)所在的各區(qū)縣。主要原因在于這些區(qū)域內(nèi)部的地形復(fù)雜,高差各異,地形起伏度值都較大。而“低—低”聚集區(qū)域主要集中在四川盆地和廣西山地丘陵性盆地所在的縣域單元。主要原因是區(qū)域內(nèi)部的四川盆地西部為地勢平坦的平原,中部為方山丘陵,東部為低山丘陵,整體地形起伏較?。欢鴱V西山地中丘陵平地所占地面積較多,整體地形起伏相對也不大,區(qū)域內(nèi)地形起伏度值表現(xiàn)為低值與低值的集聚。呈現(xiàn)空間負(fù)相關(guān)性的“低—高”聚集只存在于川西甘南的碌曲縣、若爾蓋縣及紅原縣的少部分地區(qū),表明該區(qū)域內(nèi)部地形起伏度之間存在較大的空間異質(zhì)性。其余區(qū)縣單元地形起伏度的相關(guān)性都不顯著??傮w而言,西南地區(qū)的中部山區(qū)以及成都平原區(qū)、廣西丘陵區(qū)各縣域單元之間存在著明顯的聚集格局。

    圖4 中國西南地區(qū)地形起伏度的LISA 集聚

    4 討論和結(jié)論

    (1) 通過Python編程實現(xiàn)地形起伏度的自動提取和最佳統(tǒng)計單元的確定,操作簡單高效且減少大量不必要的緩存文件。以往的常規(guī)算法通常產(chǎn)生大量緩存文件,嚴(yán)重限制了窗口數(shù)量的設(shè)置[28,20]。通過減少窗口數(shù)量或者增加窗口步距,雖然能夠減少計算量,但因此也造成了結(jié)果不夠精確[19,29-31]。近年來,Python模塊雖逐漸應(yīng)用于地形起伏度的自動提取,但其在最佳窗口確定上應(yīng)用非常少見[17]。本研究通過Python編程,完成地形起伏度自動提取的同時,還通過均值變點(diǎn)分析法實現(xiàn)最佳窗口的判定。該方法的采用,高效低耗且有效避免由于人為操作產(chǎn)生的誤差,從而有效保證了地形起伏度提取精確度。

    (2) 影響地形起伏度提取最佳適宜尺度的因素較多,窗口大小及數(shù)量的設(shè)置,地貌類型的差異以及研究區(qū)面積大小等是影響區(qū)域地形起伏度最佳適宜尺度的主要因素。曹偉超等[32]計算得出中國西南地區(qū)最佳統(tǒng)計單元為2.34 km2,與本文所得結(jié)果2.43 km2相近。但該文的分析方法存在2個不足,一是窗口數(shù)量僅為41×41,數(shù)據(jù)量較?。欢枪拯c(diǎn)的判定是在均值變點(diǎn)分析法的基礎(chǔ)上通過人為主觀確定。上述2個方面都將影響研究結(jié)果的精確度和客觀性。而本文采用Python編程,較好地避免了窗口數(shù)量較大帶來的計算困擾以及人為判定拐點(diǎn)所導(dǎo)致的誤差,因此能夠相對客觀地求取最佳窗口。值得特別指出的是,韓海輝等[29]認(rèn)為青藏高原的最佳統(tǒng)計面積為1.17 km2,該結(jié)果與本文所得的2.25 km2相差較大。究其原因,除了該文采用的最大窗口為32×32,與本文183×183相去甚遠(yuǎn)以外,兩篇文章的研究區(qū)范圍也存在很大差異,韓海輝以整個青藏高原為研究區(qū),而本文只關(guān)注青藏高原中東部地區(qū)(多為高原面)。

    (3) 西南地區(qū)地形起伏度以中小起伏為主,低起伏次之。其中,中部的橫斷山區(qū)地形起伏度最大,東南部的喀斯特地貌區(qū)次之,最低值為成都平原。地形起伏度的空間自相關(guān)分析結(jié)果中,“高高”與“低低”正相關(guān)的集聚地區(qū)也與實際地形起伏度分布的區(qū)域相符合。為了盡可能保證最佳分析尺度的可信度,以往的研究或通過采取不同分辨率的基礎(chǔ)數(shù)據(jù)以兼顧原始數(shù)據(jù)的尺度效應(yīng)[21,33],或設(shè)置不同窗口類型以兼顧?quán)徲蝾愋蛯ψ罴逊治龀叨鹊挠绊慬18],或采取多種樣地類型以兼顧地貌類型的多樣性[34]。但就某一個研究案例而言,很少同時兼顧多種因素,以至影響最佳分析尺度的精確度。本研究基于面積廣闊的中國西南地區(qū),采用30 m分辨率的原數(shù)據(jù)并設(shè)置步長為1,終止窗口為183×183的窗口數(shù)據(jù),有效考慮了與分析尺度相關(guān)的多種影響因素,從而提高了最佳分析尺度確定的精確度。

    [1] 劉振東,涂漢明.中國地勢起伏度統(tǒng)計單元的初步研究[J].熱帶地理,1989,9(1):31-38.

    [2] Deffontaines B, Lee J C, Angelier J, et al. New geomorphic data on the active Taiwan orogen-a multisource approach[J]. Journal of Geophysical Research-Solid Earth, 1994,99(B10):20243-20266.

    [3] Kuhni A, Pfiffner O A. The relief of the Swiss Alps and adjacent areas and its relation to lithology and structure: Topographic analysis from a 250 m DEM [J]. Geomorphology, 2001,41(4):285-307.

    [4] 張會平,楊農(nóng),張岳橋,等.岷江水系流域地貌特征及其構(gòu)造指示意義[J].第四紀(jì)研究,2006,26(1):126-135.

    [5] 涂漢明,劉振東.中國地勢起伏度研究[J].測繪學(xué)報,1991,20(4):311-319.

    [6] 哈凱,丁慶龍,門明新,等.山地丘陵區(qū)土地利用分布及其與地形因子關(guān)系:以河北省懷來縣為例[J].地理研究,2015,34(5):909-921.

    [7] 游珍,楊艷昭,姜魯光,等.基于DEM數(shù)據(jù)的瀾滄江—湄公河流域地形起伏度研究[J].云南大學(xué)學(xué)報:自然科學(xué)版,2012,34(4):393-400.

    [8] 朱紅春,陳楠,劉海英,等.自1∶10000比例尺DEM提取地形起伏度:以陜北黃土高原的試驗為例[J].測繪科學(xué),2005,30(4):86-88.

    [9] 范建容,張子瑜,李立華.四川省山地類型界定與山區(qū)類型劃分[J].地理研究,2015,34(1):65-73.

    [10] Prima O D A, Echigo A, Yokoyama R, et al. Supervised landform classification of Northeast Honshu from DEM-derived thematic maps[J]. Geomorphology, 2006,78(3/4):373-386.

    [11] Zevenbergen L W, Thorne C R. Quantitative-analysis of land surface-topography[J]. Earth Surface Processes and Landforms, 1987,12(1):47-56.

    [12] 江曉波.中國山地范圍界定的初步意見[J].山地學(xué)報,2008,26(2):129-136.

    [13] Liu Chun, Sun Weiwei, Wu Hangbin. Determination of complexity factor and its relationship with accuracy of representation for DEM terrain[J]. Geo-spatial Information Science, 2010,13(4):249-256.

    [14] 郎玲玲,程維明,朱啟疆,等.多尺度DEM提取地勢起伏度的對比分析:以福建低山丘陵區(qū)為例[J].地球信息科學(xué),2007,9(6):1-6.

    [15] Feng Zhiming, Tang Yan, Yang Yanzhao, et al. The relief degree of land surface in China and its correlation with population distribution [J]. Acta Geographica Sinica, 2007,62(2):1073-1082.

    [16] 王玲,同小娟.基于變點(diǎn)分析的地形起伏度研究[J].地理與地理信息科學(xué),2007,23(6):65-67.

    [17] 王讓虎,張樹文,蒲羅曼,等.基于ASTER GDEM和均值變點(diǎn)分析的中國東北地形起伏度研究[J].干旱區(qū)資源與環(huán)境,2016,30(6):49-54.

    [18] 陳學(xué)兄,張小軍,常慶瑞.陜西省地形起伏度最佳計算單元研究[J].水土保持通報,2016,36(3):265-270.

    [19] 陳田田,彭立,劉邵權(quán),等.基于GIS的橫斷山區(qū)地形起伏度與人口和經(jīng)濟(jì)的關(guān)系[J].中國科學(xué)院大學(xué)學(xué)報,2016,33(4):505-512.

    [20] 張軍,李曉東,陳春艷,等.新疆地勢起伏度的分析研究[J].蘭州大學(xué)學(xué)報:自然科學(xué)版,2008,44(S1):10-13.

    [21] 張偉,李愛農(nóng).基于DEM的中國地形起伏度適宜計算尺度研究[J].地理與地理信息科學(xué),2012,28(4):8-12.

    [22] Cliff A D, Ord J K. Spatial Processes: Models & Applications[M]. London: Taylor & Francis, 1981.

    [23] 孟斌,王勁峰,張文宗,等.基于空間分析方法的中國區(qū)域差異研究[J].地理科學(xué),2005,25(4):11-18.

    [24] 谷建立,張海濤,陳家贏,等.基于DEM的縣域土地利用空間自相關(guān)格局分析[J].農(nóng)業(yè)工程學(xué)報,2012,28(23):216-224.

    [25] 陳彥光.基于Moran統(tǒng)計量的空間自相關(guān)理論發(fā)展和方法改進(jìn)[J].地理研究,2009,28(6):1449-1463.

    [26] Li Hongfei, Calder C A, Cressie N. Beyond Moran's I:Testing for spatial dependence based on the spatial autoregressive model[J]. Geographical Analysis, 2007,39(4):357-375.

    [27] Anselin L. Local indicators of spatial association-LISA[J]. Geographical Analysis, 1995,27(2):93-115.

    [28] 王康,何俊仕,于德浩,等.采用ArcGIS平臺的地勢起伏度自動提取技術(shù)研究[J].沈陽理工大學(xué)學(xué)報,2013,32(2):63-67.

    [29] 韓海輝,高婷,易歡,等.基于變點(diǎn)分析法提取地勢起伏度:以青藏高原為例[J].地理科學(xué),2012,32(1):101-104.

    [30] 范鵬宇,朱校娟,郭嘯川.基于DEM的地形起伏度最佳計算尺度的研究:以仙居縣為例[J].地礦測繪,2016,32(4):36-38.

    [31] 王巖,劉少峰.基于DEM的青海貴德地區(qū)地形起伏度的研究[J].地質(zhì)通報,2008,27(4):2117-2121.

    [32] 曹偉超,陶和平,孔博,等.基于DEM數(shù)據(jù)分割的西南地區(qū)地貌形態(tài)自動識別研究[J].中國水土保持,2011(3):38-41.

    [33] 張錦明,游雄.地形起伏度最佳分析區(qū)域研究[J].測繪科學(xué)技術(shù)學(xué)報,2011,28(5):369-373.

    [34] 劉新華,楊勤科,湯國安.中國地形起伏度的提取及在水土流失定量評價中的應(yīng)用[J].水土保持通報,2001,21(1):57-59.

    猜你喜歡
    變點(diǎn)度值西南地區(qū)
    基于MCI指數(shù)的西南地區(qū)近60年夏季干旱特征
    探討公路項目路基連續(xù)壓實質(zhì)量檢測技術(shù)
    西南地區(qū)干濕演變特征及其對ENSO事件的響應(yīng)
    回歸模型參數(shù)的變點(diǎn)檢測方法研究
    正態(tài)分布序列均值變點(diǎn)檢測的貝葉斯方法
    基于二元分割的多變點(diǎn)估計
    獨(dú)立二項分布序列變點(diǎn)的識別方法
    山西隊?wèi)蛳蛭髂系貐^(qū)的傳播
    中華戲曲(2019年2期)2019-02-06 06:54:16
    無線傳輸中短碼長噴泉碼的度分布優(yōu)化算法*
    Завершено строительство крупнейшего комплексного транспортного узла на юго
    --западе Китая
    中亞信息(2016年7期)2016-10-20 01:41:30
    真实男女啪啪啪动态图| 国产探花极品一区二区| 欧美成人免费av一区二区三区| 春色校园在线视频观看| 欧美性猛交黑人性爽| 欧美日韩精品成人综合77777| 欧美区成人在线视频| 日韩在线高清观看一区二区三区| 欧美成人免费av一区二区三区| 能在线免费观看的黄片| 亚洲欧美日韩卡通动漫| 天堂√8在线中文| 亚洲欧美精品专区久久| 变态另类成人亚洲欧美熟女| 精品久久久噜噜| 国产精品一二三区在线看| 最后的刺客免费高清国语| 精品久久国产蜜桃| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久亚洲av鲁大| 欧美最新免费一区二区三区| 欧美日韩国产亚洲二区| 国产乱人偷精品视频| 干丝袜人妻中文字幕| 99久国产av精品国产电影| 超碰av人人做人人爽久久| 国产亚洲91精品色在线| 色吧在线观看| 精品免费久久久久久久清纯| 啦啦啦韩国在线观看视频| 极品教师在线视频| 在现免费观看毛片| or卡值多少钱| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品国产成人久久av| 久久精品夜夜夜夜夜久久蜜豆| 久久久午夜欧美精品| 国产午夜福利久久久久久| www.av在线官网国产| 久久亚洲国产成人精品v| 真实男女啪啪啪动态图| 啦啦啦啦在线视频资源| 亚洲欧美中文字幕日韩二区| 97热精品久久久久久| or卡值多少钱| 国产成人精品婷婷| 男女下面进入的视频免费午夜| 久久热精品热| 亚洲第一区二区三区不卡| 少妇熟女aⅴ在线视频| 亚洲,欧美,日韩| 国产午夜精品一二区理论片| 99热这里只有是精品50| 中文字幕免费在线视频6| 亚洲成a人片在线一区二区| 国产老妇伦熟女老妇高清| 国产精品久久久久久av不卡| 欧美性猛交黑人性爽| 亚洲国产欧美在线一区| 婷婷六月久久综合丁香| 久久99蜜桃精品久久| 国产精品一区二区性色av| 不卡视频在线观看欧美| kizo精华| 中文资源天堂在线| 国产成人a区在线观看| 欧美色欧美亚洲另类二区| 99riav亚洲国产免费| 偷拍熟女少妇极品色| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| av在线播放精品| 激情 狠狠 欧美| 亚洲国产欧美在线一区| 亚洲av免费在线观看| 成人国产麻豆网| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久亚洲av鲁大| 三级国产精品欧美在线观看| 一级黄色大片毛片| 亚洲成人av在线免费| 国产精品99久久久久久久久| 精品久久久久久久人妻蜜臀av| 国产在线精品亚洲第一网站| 干丝袜人妻中文字幕| 三级毛片av免费| 日本撒尿小便嘘嘘汇集6| 麻豆av噜噜一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产精品久久视频播放| 最近视频中文字幕2019在线8| 97在线视频观看| 中国国产av一级| 久久久精品94久久精品| 国产私拍福利视频在线观看| 免费无遮挡裸体视频| 欧美激情国产日韩精品一区| 国产在线精品亚洲第一网站| 极品教师在线视频| 青春草国产在线视频 | 久久精品国产99精品国产亚洲性色| 国产视频内射| 韩国av在线不卡| 国产人妻一区二区三区在| 国产真实乱freesex| 在线免费十八禁| 在线免费十八禁| 国产亚洲欧美98| 欧美成人精品欧美一级黄| 校园春色视频在线观看| 国产av一区在线观看免费| 午夜久久久久精精品| 久久人妻av系列| 亚洲三级黄色毛片| 欧美潮喷喷水| 亚洲成av人片在线播放无| 国产伦精品一区二区三区四那| 国产成人午夜福利电影在线观看| 成人午夜精彩视频在线观看| 国产欧美日韩精品一区二区| 亚洲av中文字字幕乱码综合| 国产精品野战在线观看| 国产伦精品一区二区三区四那| 日韩欧美在线乱码| 成人午夜精彩视频在线观看| 亚洲国产日韩欧美精品在线观看| 在线播放无遮挡| 日韩,欧美,国产一区二区三区 | 尾随美女入室| 成人永久免费在线观看视频| 国产精品一区二区三区四区免费观看| 国产伦在线观看视频一区| 不卡视频在线观看欧美| 少妇高潮的动态图| 少妇高潮的动态图| 夫妻性生交免费视频一级片| 久久久久久久久久久免费av| 全区人妻精品视频| av.在线天堂| 久久热精品热| 少妇人妻一区二区三区视频| 热99在线观看视频| 高清日韩中文字幕在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 黄色日韩在线| 九九在线视频观看精品| 99久久久亚洲精品蜜臀av| 亚洲国产高清在线一区二区三| 22中文网久久字幕| 欧美bdsm另类| av又黄又爽大尺度在线免费看 | 欧美一区二区精品小视频在线| 国产精品福利在线免费观看| 国产一区二区在线观看日韩| av免费在线看不卡| 国产精品三级大全| 偷拍熟女少妇极品色| 好男人视频免费观看在线| 久久久国产成人免费| 99热全是精品| 国产精品久久久久久久久免| 99久久人妻综合| 啦啦啦韩国在线观看视频| 少妇被粗大猛烈的视频| 少妇被粗大猛烈的视频| 在线免费观看的www视频| 色哟哟·www| 丝袜美腿在线中文| 亚洲国产色片| 日本黄色片子视频| 亚洲图色成人| 啦啦啦观看免费观看视频高清| 久久亚洲国产成人精品v| 三级经典国产精品| 夜夜看夜夜爽夜夜摸| 26uuu在线亚洲综合色| 天美传媒精品一区二区| 日本黄色片子视频| 久久久a久久爽久久v久久| 亚洲丝袜综合中文字幕| 国产精品乱码一区二三区的特点| 人体艺术视频欧美日本| 亚洲成人精品中文字幕电影| 99热只有精品国产| 免费看美女性在线毛片视频| 在线观看午夜福利视频| 热99re8久久精品国产| 别揉我奶头 嗯啊视频| 日韩欧美在线乱码| 成人亚洲欧美一区二区av| 国产乱人视频| 国产蜜桃级精品一区二区三区| 成年女人永久免费观看视频| 最后的刺客免费高清国语| 国产单亲对白刺激| 欧美一级a爱片免费观看看| 尾随美女入室| 日韩欧美 国产精品| 久久精品夜夜夜夜夜久久蜜豆| 成人漫画全彩无遮挡| 欧美xxxx黑人xx丫x性爽| 中文字幕精品亚洲无线码一区| 老女人水多毛片| 国产午夜精品论理片| 国产极品精品免费视频能看的| 一区福利在线观看| 国产成人福利小说| videossex国产| 亚洲欧美精品综合久久99| 天堂av国产一区二区熟女人妻| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 99热只有精品国产| 日日干狠狠操夜夜爽| 99在线视频只有这里精品首页| 九九爱精品视频在线观看| 欧美人与善性xxx| 日韩欧美 国产精品| 99热网站在线观看| 一夜夜www| 欧美xxxx性猛交bbbb| 日韩大尺度精品在线看网址| 内射极品少妇av片p| 日韩 亚洲 欧美在线| 国产精品免费一区二区三区在线| 久久精品国产亚洲av香蕉五月| 日韩三级伦理在线观看| 精品久久久久久久久av| 91久久精品国产一区二区三区| 日本av手机在线免费观看| 99久久无色码亚洲精品果冻| 最近的中文字幕免费完整| 九九爱精品视频在线观看| 欧美一级a爱片免费观看看| 熟妇人妻久久中文字幕3abv| 国产熟女欧美一区二区| 中文字幕av成人在线电影| 草草在线视频免费看| 国产av不卡久久| 91久久精品国产一区二区成人| 日韩欧美精品免费久久| 最近2019中文字幕mv第一页| 久久99热6这里只有精品| 国产亚洲av片在线观看秒播厂 | 亚洲人成网站高清观看| 国产在线男女| 三级经典国产精品| 久久欧美精品欧美久久欧美| 日韩欧美一区二区三区在线观看| 综合色丁香网| 久久九九热精品免费| 如何舔出高潮| 国产欧美日韩精品一区二区| 免费av不卡在线播放| 亚洲成人久久性| 国产高清激情床上av| 国产黄a三级三级三级人| 我要看日韩黄色一级片| 国产亚洲av嫩草精品影院| 国产亚洲欧美98| 国产午夜精品一二区理论片| 国产黄片美女视频| 只有这里有精品99| 久久中文看片网| 青青草视频在线视频观看| 亚洲最大成人手机在线| 国产中年淑女户外野战色| 十八禁国产超污无遮挡网站| 岛国毛片在线播放| 22中文网久久字幕| 深夜精品福利| 成人午夜高清在线视频| 国产精品一区www在线观看| 亚洲久久久久久中文字幕| 亚洲丝袜综合中文字幕| 久久久久久久久久久免费av| 在线免费十八禁| 99久久成人亚洲精品观看| 久久久久久久久中文| av天堂中文字幕网| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲av香蕉五月| 亚洲综合色惰| 国产精品日韩av在线免费观看| 国产精品久久视频播放| 久久久久性生活片| 国产色爽女视频免费观看| 老女人水多毛片| 国产av一区在线观看免费| 日日摸夜夜添夜夜添av毛片| 小说图片视频综合网站| av国产免费在线观看| av在线天堂中文字幕| 欧美xxxx性猛交bbbb| 啦啦啦观看免费观看视频高清| 中文亚洲av片在线观看爽| 国产精品一及| 国产成人a∨麻豆精品| 极品教师在线视频| 伊人久久精品亚洲午夜| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 51国产日韩欧美| 麻豆一二三区av精品| 成人毛片a级毛片在线播放| 变态另类丝袜制服| 日本欧美国产在线视频| 国产 一区 欧美 日韩| 日本黄色视频三级网站网址| 99热只有精品国产| 最近手机中文字幕大全| 久久精品国产鲁丝片午夜精品| 波野结衣二区三区在线| 免费电影在线观看免费观看| 好男人视频免费观看在线| 久久久久久久久久久丰满| 国语自产精品视频在线第100页| 国产成人91sexporn| 一边摸一边抽搐一进一小说| 别揉我奶头 嗯啊视频| 在线观看免费视频日本深夜| 精品一区二区免费观看| 亚洲五月天丁香| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇中文字幕五十中出| 免费大片18禁| 一本一本综合久久| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩国产亚洲二区| 成人毛片a级毛片在线播放| 菩萨蛮人人尽说江南好唐韦庄 | 成人毛片60女人毛片免费| 搡女人真爽免费视频火全软件| 卡戴珊不雅视频在线播放| 亚洲美女搞黄在线观看| 看十八女毛片水多多多| 菩萨蛮人人尽说江南好唐韦庄 | 精品一区二区三区人妻视频| 婷婷六月久久综合丁香| 校园人妻丝袜中文字幕| 亚洲人与动物交配视频| 黄色欧美视频在线观看| 人人妻人人澡人人爽人人夜夜 | 国语自产精品视频在线第100页| 国产精品福利在线免费观看| АⅤ资源中文在线天堂| 禁无遮挡网站| 亚洲不卡免费看| 91精品一卡2卡3卡4卡| 12—13女人毛片做爰片一| 久久精品久久久久久噜噜老黄 | 亚洲欧洲国产日韩| 深爱激情五月婷婷| 中文字幕精品亚洲无线码一区| 青春草视频在线免费观看| av福利片在线观看| 精华霜和精华液先用哪个| 一级二级三级毛片免费看| av在线播放精品| 热99re8久久精品国产| 级片在线观看| 你懂的网址亚洲精品在线观看 | 直男gayav资源| 国产白丝娇喘喷水9色精品| 成人三级黄色视频| 亚洲无线在线观看| 久久九九热精品免费| 午夜久久久久精精品| 免费av不卡在线播放| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 狂野欧美白嫩少妇大欣赏| 老师上课跳d突然被开到最大视频| 网址你懂的国产日韩在线| av国产免费在线观看| 69人妻影院| 在线播放国产精品三级| 国产一级毛片七仙女欲春2| 国产精品免费一区二区三区在线| 69人妻影院| 老司机福利观看| 不卡视频在线观看欧美| 久久人人精品亚洲av| 爱豆传媒免费全集在线观看| 日韩高清综合在线| 午夜精品国产一区二区电影 | 国产精品一区www在线观看| 亚洲av中文字字幕乱码综合| 亚洲欧美精品专区久久| 搡老妇女老女人老熟妇| 日本一二三区视频观看| 国产精品蜜桃在线观看 | 91精品一卡2卡3卡4卡| 日本免费一区二区三区高清不卡| 亚洲欧美日韩无卡精品| av在线天堂中文字幕| 亚洲一区高清亚洲精品| 成人永久免费在线观看视频| 欧美+亚洲+日韩+国产| 一级毛片电影观看 | 欧美最新免费一区二区三区| 在线免费观看的www视频| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 久久人人精品亚洲av| 青春草亚洲视频在线观看| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 黄片无遮挡物在线观看| 日本免费a在线| 亚洲成a人片在线一区二区| 色5月婷婷丁香| 成人三级黄色视频| av在线播放精品| 九九爱精品视频在线观看| 色综合色国产| 国产成人91sexporn| 能在线免费看毛片的网站| 久久人妻av系列| 夜夜看夜夜爽夜夜摸| 国产成年人精品一区二区| 1024手机看黄色片| 国产精品电影一区二区三区| av在线蜜桃| 一级毛片电影观看 | 黄色日韩在线| 日韩精品有码人妻一区| av在线播放精品| 成熟少妇高潮喷水视频| 黄色一级大片看看| 成人国产麻豆网| 一区二区三区免费毛片| 国产大屁股一区二区在线视频| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 国产在线精品亚洲第一网站| 三级国产精品欧美在线观看| 变态另类丝袜制服| 久99久视频精品免费| 亚洲欧美成人精品一区二区| 国产真实伦视频高清在线观看| 日本av手机在线免费观看| 日韩制服骚丝袜av| 国产成人a∨麻豆精品| 又爽又黄a免费视频| 欧美日韩综合久久久久久| 亚洲av男天堂| 在线免费十八禁| 欧美高清成人免费视频www| 中文欧美无线码| 99国产精品一区二区蜜桃av| 中文在线观看免费www的网站| 国产精品一区二区性色av| 激情 狠狠 欧美| 亚洲精品456在线播放app| 免费一级毛片在线播放高清视频| 日产精品乱码卡一卡2卡三| а√天堂www在线а√下载| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 五月伊人婷婷丁香| 99在线人妻在线中文字幕| 亚洲成a人片在线一区二区| 看免费成人av毛片| 免费av毛片视频| 国内精品宾馆在线| 日韩成人伦理影院| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 国产av麻豆久久久久久久| 国产三级在线视频| 亚洲五月天丁香| 国产久久久一区二区三区| 婷婷亚洲欧美| 国产乱人视频| 一个人免费在线观看电影| 99热这里只有是精品50| 免费看a级黄色片| 99久久无色码亚洲精品果冻| 丝袜喷水一区| 神马国产精品三级电影在线观看| 一夜夜www| 日本一本二区三区精品| 国产精品久久久久久亚洲av鲁大| 91在线精品国自产拍蜜月| av黄色大香蕉| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 欧洲精品卡2卡3卡4卡5卡区| 国产精品一区二区性色av| 久久精品综合一区二区三区| av黄色大香蕉| 亚洲国产精品合色在线| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 乱码一卡2卡4卡精品| 欧美精品一区二区大全| 一级毛片电影观看 | 最近最新中文字幕大全电影3| 又粗又爽又猛毛片免费看| 最近中文字幕高清免费大全6| 十八禁国产超污无遮挡网站| 黄片无遮挡物在线观看| 中国美女看黄片| 国产午夜精品一二区理论片| 男人舔女人下体高潮全视频| 久久久色成人| 中文资源天堂在线| 国产高清激情床上av| 成人性生交大片免费视频hd| 亚洲成人av在线免费| а√天堂www在线а√下载| 好男人在线观看高清免费视频| 能在线免费看毛片的网站| 亚洲最大成人手机在线| 亚洲人成网站在线观看播放| 国产真实伦视频高清在线观看| 午夜福利在线观看吧| 免费av观看视频| 国产国拍精品亚洲av在线观看| 亚洲av电影不卡..在线观看| 国产精品久久久久久亚洲av鲁大| 最后的刺客免费高清国语| 欧美高清成人免费视频www| 变态另类丝袜制服| 丝袜喷水一区| 99久久中文字幕三级久久日本| 午夜精品国产一区二区电影 | 蜜臀久久99精品久久宅男| 国产精品久久久久久精品电影小说 | 日本黄色片子视频| 国产精品久久久久久久电影| 日韩强制内射视频| 一区福利在线观看| 少妇人妻精品综合一区二区 | 成人毛片60女人毛片免费| 亚洲国产高清在线一区二区三| 日韩强制内射视频| 日韩欧美国产在线观看| 黄片无遮挡物在线观看| 一本久久中文字幕| 看片在线看免费视频| 99热只有精品国产| 一级毛片我不卡| 婷婷色av中文字幕| 91午夜精品亚洲一区二区三区| 国内精品宾馆在线| 人人妻人人澡人人爽人人夜夜 | 一级毛片aaaaaa免费看小| 午夜a级毛片| 99在线视频只有这里精品首页| 黄色一级大片看看| 在线观看66精品国产| a级毛片a级免费在线| 国产在视频线在精品| 深夜a级毛片| 日日干狠狠操夜夜爽| 在线国产一区二区在线| 美女国产视频在线观看| 国产69精品久久久久777片| 亚洲中文字幕一区二区三区有码在线看| 少妇裸体淫交视频免费看高清| 国产日韩欧美在线精品| 亚洲综合色惰| 久久久精品大字幕| 亚洲精品国产成人久久av| 国产日韩欧美在线精品| 女人十人毛片免费观看3o分钟| 国产熟女欧美一区二区| 久久99热6这里只有精品| 欧美bdsm另类| 乱码一卡2卡4卡精品| 黄色视频,在线免费观看| 不卡视频在线观看欧美| av在线播放精品| 一级黄色大片毛片| 亚洲国产色片| avwww免费| 真实男女啪啪啪动态图| 国产精品久久久久久精品电影| kizo精华| 69人妻影院| 人妻制服诱惑在线中文字幕| 国产成人影院久久av| 国产爱豆传媒在线观看| 99久久成人亚洲精品观看| 精品久久久久久久久av| 国产亚洲av片在线观看秒播厂 | 三级国产精品欧美在线观看| 99久久九九国产精品国产免费| 国产高清激情床上av| 一本久久精品| 国内久久婷婷六月综合欲色啪| 大型黄色视频在线免费观看| 亚洲va在线va天堂va国产| 国产精品一区二区在线观看99 | 日韩欧美一区二区三区在线观看| 蜜桃久久精品国产亚洲av| 22中文网久久字幕| 久久久色成人| 不卡视频在线观看欧美| 国产精品日韩av在线免费观看| 直男gayav资源| 日韩一本色道免费dvd| 能在线免费看毛片的网站| 天天一区二区日本电影三级| 一边亲一边摸免费视频| 久久久久性生活片| 两个人的视频大全免费| 嫩草影院入口| 国产色婷婷99| 深夜a级毛片|