• <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麻豆狂野| 18在线观看网站| 午夜福利一区二区在线看| 两个人免费观看高清视频| 一进一出好大好爽视频| 大香蕉久久网| 视频区欧美日本亚洲| 男女之事视频高清在线观看| 亚洲av电影在线进入| 一级,二级,三级黄色视频| 男女免费视频国产| 亚洲一码二码三码区别大吗| 99国产精品99久久久久| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美在线二视频 | 久久青草综合色| 窝窝影院91人妻| 在线十欧美十亚洲十日本专区| av电影中文网址| 丝瓜视频免费看黄片| 欧美黑人欧美精品刺激| 成人黄色视频免费在线看| 夫妻午夜视频| 一本大道久久a久久精品| 女人被躁到高潮嗷嗷叫费观| 成人亚洲精品一区在线观看| 免费在线观看黄色视频的| 人妻丰满熟妇av一区二区三区 | 成在线人永久免费视频| 最新在线观看一区二区三区| 757午夜福利合集在线观看| 亚洲美女黄片视频| 国产成人精品无人区| 精品熟女少妇八av免费久了| 免费一级毛片在线播放高清视频 | 老熟妇乱子伦视频在线观看| 一区二区日韩欧美中文字幕| 午夜福利在线观看吧| 极品人妻少妇av视频| 美女 人体艺术 gogo| 少妇猛男粗大的猛烈进出视频| 久久久久久免费高清国产稀缺| 日韩精品免费视频一区二区三区| 性色av乱码一区二区三区2| 十分钟在线观看高清视频www| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 美女高潮喷水抽搐中文字幕| 欧美日韩视频精品一区| 国产aⅴ精品一区二区三区波| 老熟妇乱子伦视频在线观看| av视频免费观看在线观看| 在线看a的网站| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 嫩草影视91久久| 成人免费观看视频高清| 欧美人与性动交α欧美精品济南到| 亚洲欧洲精品一区二区精品久久久| 国产1区2区3区精品| 极品人妻少妇av视频| 亚洲中文字幕日韩| 色综合婷婷激情| 亚洲avbb在线观看| 国产高清国产精品国产三级| 成人18禁在线播放| 久久午夜综合久久蜜桃| 夫妻午夜视频| 国产精品久久电影中文字幕 | 美女 人体艺术 gogo| 大香蕉久久网| 亚洲精品一二三| 少妇被粗大的猛进出69影院| 久久久水蜜桃国产精品网| 男女免费视频国产| 午夜精品在线福利| 最新在线观看一区二区三区| 国产精品免费大片| 久久久久精品国产欧美久久久| 两个人看的免费小视频| 真人做人爱边吃奶动态| netflix在线观看网站| 亚洲自偷自拍图片 自拍| 欧美精品亚洲一区二区| 韩国av一区二区三区四区| 亚洲 欧美一区二区三区| 丁香六月欧美| 国产精品 欧美亚洲| 狂野欧美激情性xxxx| 国产成人欧美在线观看 | 亚洲av成人不卡在线观看播放网| 欧美av亚洲av综合av国产av| 亚洲精品美女久久久久99蜜臀| 久久99一区二区三区| 99热网站在线观看| 欧美久久黑人一区二区| 国产亚洲欧美在线一区二区| 日韩欧美国产一区二区入口| 国产又爽黄色视频| av视频免费观看在线观看| 亚洲第一青青草原| 黄频高清免费视频| 在线观看日韩欧美| 丰满人妻熟妇乱又伦精品不卡| 成人精品一区二区免费| 欧美日韩精品网址| 十八禁网站免费在线| 中文字幕人妻熟女乱码| 久久精品国产综合久久久| 欧美av亚洲av综合av国产av| 欧美午夜高清在线| 久久九九热精品免费| 亚洲精品久久午夜乱码| 中文字幕制服av| 成在线人永久免费视频| 免费av中文字幕在线| 老司机深夜福利视频在线观看| 国产精品久久久av美女十八| 天天操日日干夜夜撸| 校园春色视频在线观看| 婷婷成人精品国产| 老司机影院毛片| 久久精品国产99精品国产亚洲性色 | 97人妻天天添夜夜摸| 在线免费观看的www视频| 成人手机av| 一级毛片精品| 天天躁夜夜躁狠狠躁躁| 交换朋友夫妻互换小说| 一进一出抽搐动态| 手机成人av网站| a级片在线免费高清观看视频| 19禁男女啪啪无遮挡网站| av福利片在线| 亚洲精品国产色婷婷电影| 久久久久久久久免费视频了| 精品亚洲成国产av| 亚洲一区二区三区欧美精品| 久久人妻av系列| 欧美亚洲日本最大视频资源| 一夜夜www| 黑人猛操日本美女一级片| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 日本vs欧美在线观看视频| 久久久久国产精品人妻aⅴ院 | 精品国内亚洲2022精品成人 | 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| av中文乱码字幕在线| av网站免费在线观看视频| 高清在线国产一区| 亚洲成人免费av在线播放| 国产精品 国内视频| 看免费av毛片| 超碰97精品在线观看| 亚洲久久久国产精品| 国产精品免费大片| 亚洲第一欧美日韩一区二区三区| 国产精品 欧美亚洲| 脱女人内裤的视频| 超碰成人久久| 日韩制服丝袜自拍偷拍| 精品乱码久久久久久99久播| 国产伦人伦偷精品视频| 操出白浆在线播放| 亚洲欧美日韩高清在线视频| 动漫黄色视频在线观看| 欧美午夜高清在线| 欧美日韩成人在线一区二区| 老司机午夜福利在线观看视频| 久久 成人 亚洲| 国产亚洲精品第一综合不卡| 中文欧美无线码| 国产在视频线精品| 人人妻人人爽人人添夜夜欢视频| 欧美精品一区二区免费开放| 日韩欧美在线二视频 | 国产不卡一卡二| 国产一区二区三区在线臀色熟女 | 精品国产一区二区久久| 国产在线精品亚洲第一网站| 成人特级黄色片久久久久久久| 久久狼人影院| 国产av又大| 国产亚洲精品一区二区www | 精品卡一卡二卡四卡免费| 天天影视国产精品| 亚洲色图 男人天堂 中文字幕| 精品卡一卡二卡四卡免费| 99re在线观看精品视频| 精品一区二区三区四区五区乱码| 一级黄色大片毛片| 精品久久久久久久毛片微露脸| 精品久久蜜臀av无| 激情视频va一区二区三区| 国产精品自产拍在线观看55亚洲 | 亚洲熟妇中文字幕五十中出 | 久久精品国产综合久久久| 日韩欧美一区二区三区在线观看 | 国产视频一区二区在线看| 久久精品人人爽人人爽视色| 黄片播放在线免费| 久久狼人影院| 三级毛片av免费| 欧美不卡视频在线免费观看 | 国产色视频综合| 亚洲一区二区三区不卡视频| 精品人妻熟女毛片av久久网站| 99热网站在线观看| 中文字幕制服av| 夜夜爽天天搞| 国产在视频线精品| 日韩熟女老妇一区二区性免费视频| 亚洲成国产人片在线观看| 亚洲中文字幕日韩| 满18在线观看网站| 王馨瑶露胸无遮挡在线观看| 黄色女人牲交| 大码成人一级视频| 天天躁夜夜躁狠狠躁躁| 在线观看免费日韩欧美大片| 757午夜福利合集在线观看| 亚洲精品美女久久av网站| 熟女少妇亚洲综合色aaa.| 午夜精品国产一区二区电影| 亚洲五月天丁香| 真人做人爱边吃奶动态| 天堂动漫精品| 国内久久婷婷六月综合欲色啪| 少妇粗大呻吟视频| 欧美亚洲 丝袜 人妻 在线| 美女福利国产在线| 欧美老熟妇乱子伦牲交| 成年女人毛片免费观看观看9 | 在线看a的网站| 夜夜躁狠狠躁天天躁| 一a级毛片在线观看| 三级毛片av免费| 麻豆乱淫一区二区| 亚洲一码二码三码区别大吗| 国产在视频线精品| av片东京热男人的天堂| 欧美日韩黄片免| 亚洲av熟女| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 成熟少妇高潮喷水视频| 免费黄频网站在线观看国产| 免费在线观看黄色视频的| 久热爱精品视频在线9| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 无遮挡黄片免费观看| 国产精品成人在线| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯 | 性色av乱码一区二区三区2| 亚洲色图综合在线观看| avwww免费| 视频区图区小说| 一级黄色大片毛片| 老鸭窝网址在线观看| 91在线观看av| 国产免费现黄频在线看| 9色porny在线观看| 国产又色又爽无遮挡免费看| 女人被狂操c到高潮| 国产高清视频在线播放一区| 中文欧美无线码| 手机成人av网站| 黑丝袜美女国产一区| av天堂久久9| 国产日韩欧美亚洲二区| 成人黄色视频免费在线看| 欧美激情久久久久久爽电影 | 国产av一区二区精品久久| 久久这里只有精品19| 超碰97精品在线观看| 最近最新免费中文字幕在线| www日本在线高清视频| 成人免费观看视频高清| 国产片内射在线| 巨乳人妻的诱惑在线观看| 亚洲人成伊人成综合网2020| 亚洲专区中文字幕在线| 一区在线观看完整版| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 精品高清国产在线一区| 亚洲av美国av| 亚洲情色 制服丝袜| 99re在线观看精品视频| 99热网站在线观看| ponron亚洲| 中文欧美无线码| 777久久人妻少妇嫩草av网站| 精品人妻熟女毛片av久久网站| 亚洲国产中文字幕在线视频| av在线播放免费不卡| 国产97色在线日韩免费| 国产野战对白在线观看| 99精品久久久久人妻精品| 天堂俺去俺来也www色官网| 黄频高清免费视频| 超碰成人久久| 亚洲精品中文字幕一二三四区| 真人做人爱边吃奶动态| 日韩熟女老妇一区二区性免费视频| 中文字幕最新亚洲高清| 成年动漫av网址| 可以免费在线观看a视频的电影网站| 国产精品久久久久久精品古装| 最近最新中文字幕大全免费视频| 欧美乱色亚洲激情| 捣出白浆h1v1| 国精品久久久久久国模美| 欧美在线黄色| 欧美另类亚洲清纯唯美| 一a级毛片在线观看| 国产单亲对白刺激| 国产精品久久久av美女十八| 欧美一级毛片孕妇| 国产精品久久久av美女十八| 中出人妻视频一区二区| 国产在线观看jvid| 中文字幕色久视频| 国产一区在线观看成人免费| 校园春色视频在线观看| 国产xxxxx性猛交| 最新在线观看一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 99久久精品国产亚洲精品| 午夜视频精品福利| 如日韩欧美国产精品一区二区三区| 国产精品98久久久久久宅男小说| 交换朋友夫妻互换小说| 久久精品亚洲熟妇少妇任你| 亚洲七黄色美女视频| 又黄又粗又硬又大视频| 极品教师在线免费播放| 国产在视频线精品| 亚洲精品中文字幕一二三四区| 一级毛片精品| 久久久国产精品麻豆| 国产精品免费大片| 欧美日韩亚洲高清精品| 亚洲中文av在线| 女人被狂操c到高潮| 亚洲va日本ⅴa欧美va伊人久久| 国产一区有黄有色的免费视频| 国产成人精品久久二区二区91| 免费av中文字幕在线| 欧美亚洲日本最大视频资源| 国产有黄有色有爽视频| 可以免费在线观看a视频的电影网站| 亚洲五月色婷婷综合| 极品教师在线免费播放| 丰满迷人的少妇在线观看| 99久久人妻综合| 国产成人影院久久av| 精品国产国语对白av| 黄色女人牲交| 国产成人免费观看mmmm| 在线观看www视频免费| 亚洲人成77777在线视频| svipshipincom国产片| 亚洲av成人不卡在线观看播放网| av网站免费在线观看视频| 亚洲在线自拍视频| 国产成人系列免费观看| 午夜精品在线福利| 色老头精品视频在线观看| 成年人免费黄色播放视频| 午夜久久久在线观看| 久久婷婷成人综合色麻豆| 满18在线观看网站| a级毛片黄视频| 亚洲av欧美aⅴ国产| 99国产精品一区二区三区| 亚洲七黄色美女视频| 亚洲一区中文字幕在线| 国产精品 国内视频| 一级黄色大片毛片| 成年人黄色毛片网站| 黄色 视频免费看| 在线观看日韩欧美| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 乱人伦中国视频| 国产欧美日韩综合在线一区二区| 欧美 日韩 精品 国产| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品第一综合不卡| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 久久久久视频综合| 亚洲中文字幕日韩| 国产精品.久久久| 国内久久婷婷六月综合欲色啪| 国产深夜福利视频在线观看| 日本a在线网址| 热re99久久国产66热| 午夜福利乱码中文字幕| 国产麻豆69| 在线观看www视频免费| 一夜夜www| 看黄色毛片网站| 91精品三级在线观看| 国产男女内射视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 91九色精品人成在线观看| 建设人人有责人人尽责人人享有的| 亚洲精品美女久久av网站| 男女高潮啪啪啪动态图| 国产男靠女视频免费网站| 中文欧美无线码| 成年人午夜在线观看视频| 欧美黑人欧美精品刺激| 校园春色视频在线观看| 国产91精品成人一区二区三区| 欧美大码av| 美女高潮喷水抽搐中文字幕| 国产精品电影一区二区三区 | 黄色丝袜av网址大全| 一本一本久久a久久精品综合妖精| 人人澡人人妻人| 美女高潮喷水抽搐中文字幕| 一区二区三区国产精品乱码| 中文字幕另类日韩欧美亚洲嫩草| 亚洲va日本ⅴa欧美va伊人久久| 男女免费视频国产| 亚洲欧美精品综合一区二区三区| 1024视频免费在线观看| 12—13女人毛片做爰片一| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看 | 国产aⅴ精品一区二区三区波| 宅男免费午夜| 成人精品一区二区免费| 亚洲国产精品一区二区三区在线| 水蜜桃什么品种好| 麻豆国产av国片精品| 91字幕亚洲| 少妇的丰满在线观看| 两个人免费观看高清视频| 黄片播放在线免费| 久久国产精品大桥未久av| 女性被躁到高潮视频| 男女下面插进去视频免费观看| 真人做人爱边吃奶动态| 亚洲av美国av| 亚洲午夜理论影院| 99国产精品一区二区蜜桃av | 99香蕉大伊视频| 亚洲少妇的诱惑av| 欧美另类亚洲清纯唯美| 亚洲第一青青草原| 亚洲中文日韩欧美视频| 日本a在线网址| 久久香蕉精品热| 精品熟女少妇八av免费久了| 亚洲午夜精品一区,二区,三区| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 国产精品九九99| 精品久久久久久久久久免费视频 | 色94色欧美一区二区| 国产成+人综合+亚洲专区| 99国产精品一区二区蜜桃av | 视频区图区小说| 村上凉子中文字幕在线| 伦理电影免费视频| 亚洲精品美女久久av网站| 亚洲一区二区三区不卡视频| 又大又爽又粗| 亚洲专区中文字幕在线| 校园春色视频在线观看| 国产一区二区三区综合在线观看| netflix在线观看网站| 精品福利永久在线观看| 女人被躁到高潮嗷嗷叫费观| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 精品久久久精品久久久| 午夜福利在线免费观看网站| 在线国产一区二区在线| 中亚洲国语对白在线视频| 久久婷婷成人综合色麻豆| 亚洲一区高清亚洲精品| 亚洲视频免费观看视频| 亚洲av成人一区二区三| 深夜精品福利| 欧美人与性动交α欧美软件| 国产成人影院久久av| 丝瓜视频免费看黄片| 老司机午夜福利在线观看视频| 黄色女人牲交| 看黄色毛片网站| 电影成人av| 三上悠亚av全集在线观看| 婷婷丁香在线五月| 一区二区日韩欧美中文字幕| 国产99久久九九免费精品| 久久人妻福利社区极品人妻图片| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区三区综合在线观看| 午夜福利影视在线免费观看| 麻豆av在线久日| 日韩制服丝袜自拍偷拍| 色尼玛亚洲综合影院| 日韩三级视频一区二区三区| 不卡av一区二区三区| 国产精华一区二区三区| av天堂在线播放| 亚洲精品国产色婷婷电影| 亚洲精品粉嫩美女一区| 9热在线视频观看99| 亚洲人成电影免费在线| 国产不卡一卡二| 久久青草综合色| 日韩欧美免费精品| 精品一区二区三区视频在线观看免费 | 国产av精品麻豆| 97人妻天天添夜夜摸| 一进一出抽搐动态| 久久精品国产亚洲av香蕉五月 | 日本a在线网址| 不卡一级毛片| 99久久国产精品久久久| 少妇被粗大的猛进出69影院| 丝袜美足系列| 午夜福利在线观看吧| 亚洲伊人色综图| 欧美中文综合在线视频| 大香蕉久久网| 国产精品一区二区精品视频观看| a在线观看视频网站| videos熟女内射| 亚洲情色 制服丝袜| 99精国产麻豆久久婷婷| 亚洲欧美日韩另类电影网站| 99久久国产精品久久久| 精品国产美女av久久久久小说| 91精品国产国语对白视频| 亚洲av熟女| 精品电影一区二区在线| 欧洲精品卡2卡3卡4卡5卡区| 精品久久久久久久毛片微露脸| 国产精品99久久99久久久不卡| 男人操女人黄网站| 日本欧美视频一区| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 国产精品综合久久久久久久免费 | cao死你这个sao货| 伊人久久大香线蕉亚洲五| 老司机福利观看| 岛国毛片在线播放| 黄色成人免费大全| 久久久精品国产亚洲av高清涩受| 一区二区日韩欧美中文字幕| 91在线观看av| 午夜精品在线福利| avwww免费| 亚洲国产欧美网| 国产精华一区二区三区| 国产xxxxx性猛交| 看片在线看免费视频| 在线观看午夜福利视频| 涩涩av久久男人的天堂| 久久热在线av| 男女免费视频国产| 老汉色∧v一级毛片| www.精华液| 国产精品自产拍在线观看55亚洲 | x7x7x7水蜜桃| 午夜福利,免费看| 国产精品久久久久成人av| 国产精品成人在线| 欧美精品高潮呻吟av久久| 黑人猛操日本美女一级片| 国产免费现黄频在线看| 亚洲中文字幕日韩| 亚洲综合色网址| 国产免费现黄频在线看| 亚洲中文字幕日韩| 成人三级做爰电影| 午夜福利免费观看在线| 在线观看www视频免费| 老司机午夜福利在线观看视频| 亚洲精品美女久久久久99蜜臀| 欧美性长视频在线观看| 99久久国产精品久久久| 在线免费观看的www视频| 亚洲国产精品合色在线| 国产极品粉嫩免费观看在线| 日韩欧美一区二区三区在线观看 | 啦啦啦 在线观看视频| 久久国产精品影院| 十分钟在线观看高清视频www| 高清视频免费观看一区二区| 久久精品亚洲精品国产色婷小说| 黄片大片在线免费观看| 免费少妇av软件| 免费观看精品视频网站| 亚洲精华国产精华精| 国产一区有黄有色的免费视频|