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

    地形分析中柵格中間層構(gòu)建的尺度適應(yīng)性

    2017-11-29 08:22:59張立朝張合朝
    測繪通報(bào) 2017年1期
    關(guān)鍵詞:中間層柵格切線

    張立朝,趙 鵬,張合朝

    (1. 河南省測繪地理信息局, 河南 鄭州 450003; 2. 信息工程大學(xué)地理空間信息學(xué)院,河南 鄭州 450052; 3. 河南省科源測繪中心,河南 鄭州 450003)

    地形分析中柵格中間層構(gòu)建的尺度適應(yīng)性

    張立朝1,趙 鵬2,張合朝3

    (1. 河南省測繪地理信息局, 河南 鄭州 450003; 2. 信息工程大學(xué)地理空間信息學(xué)院,河南 鄭州 450052; 3. 河南省科源測繪中心,河南 鄭州 450003)

    柵格中間層是指在地形分析中由DEM、DSM、矢量數(shù)據(jù)等原始數(shù)據(jù),經(jīng)過一系列柵格運(yùn)算獲取且可以重用的基礎(chǔ)柵格層。本文系統(tǒng)研究了柵格中間層構(gòu)建的尺度適應(yīng)性,定義了范圍、比例尺、分辨率和分析4種尺度,并確定了范圍、比例尺和分辨率尺度的確定原則;特別對于分析尺度,通過地貌因子計(jì)算進(jìn)行了切線、剖面曲率柵格層分析尺度適宜性試驗(yàn),分析了曲率計(jì)算中適宜的算法和分析尺度,保證了地形基礎(chǔ)因子?xùn)鸥裰虚g層具有滿足分析要求的數(shù)值精度和匹配地形起伏的形態(tài)精度。

    柵格中間層;尺度適應(yīng)性;地形分析;數(shù)值精度;形態(tài)精度

    尺度常常被定義為在研究對象或現(xiàn)象時(shí)所采用的空間或時(shí)間單位,或現(xiàn)象(過程)在空間和時(shí)間上所涉及的范圍和發(fā)生的頻率。在現(xiàn)有的地形分析尺度適應(yīng)性研究成果中,涵蓋了DEM地形分析的各個(gè)方面:①數(shù)據(jù)源尺度,如地形圖比例尺、遙感影像比例尺等對DEM地形分析的影響;②DEM內(nèi)插方法對地形分析的影響;③地形分析方法對地形因子和地形特征提取的影響分析。

    由于柵格數(shù)據(jù)在通過離散方式表達(dá)連續(xù)變化的地形表面過程中存在著尺度依賴性,因此,在柵格中間層的構(gòu)建過程中,應(yīng)合理確定柵格中間層的覆蓋范圍、相應(yīng)比例尺、柵格單元分辨率等尺度問題,以保證同一級別的柵格中間層在進(jìn)行疊置、模糊綜合評判分析時(shí),不會因?yàn)槟撤N尺度不同造成誤差甚至導(dǎo)致錯(cuò)誤分析結(jié)果;另外在進(jìn)行鄰域分析時(shí)還會遇到分析尺度大小影響分析結(jié)果的情況,也是在柵格中間層構(gòu)建過程判斷算法適應(yīng)性的關(guān)鍵問題。

    1 柵格中間層構(gòu)建的尺度

    圖1 柵格中間層構(gòu)建尺度關(guān)系

    本文定義了4種柵格中間層的構(gòu)建尺度,分別是范圍尺度、比例尺尺度、分辨率和分析尺度。地形分析尺度是非常復(fù)雜的問題,一種統(tǒng)一的尺度不可能適合所有分析對象和分析算法,應(yīng)該在諸多尺度問題的研究中找到一個(gè)相對平衡的統(tǒng)一尺度范圍,盡可能滿足柵格中間層總體的精度、算法適應(yīng)性和對各類分析對象的適應(yīng)性。上述4種柵格中間層構(gòu)建尺度的結(jié)構(gòu)關(guān)系如圖1所示。

    (1) 范圍尺度。柵格中間層數(shù)據(jù)在構(gòu)建時(shí)以地形表面為基礎(chǔ),是對地形基礎(chǔ)因子、地形特征及地形影響因素對專題應(yīng)用影響程度的數(shù)字化表達(dá),本身就具有范圍尺度含義。柵格中間層的范圍尺度直接決定著其他尺度和柵格中間層的應(yīng)用目的。

    (2) 比例尺尺度。柵格中間層的比例尺尺度體現(xiàn)在地物數(shù)據(jù)柵格化的矢量數(shù)據(jù)比例尺,柵格中間層的分辨率尺度與比例尺尺度是相互關(guān)聯(lián)的關(guān)系。

    (3) 分辨率。在柵格中間層的應(yīng)用中,各層級中間層將會進(jìn)行大量的疊置分析、模糊操作,如果分辨率不同,會導(dǎo)致格網(wǎng)計(jì)算的不匹配,就會帶來大量的內(nèi)插計(jì)算,致使冗余和效率低下。因此,如果能根據(jù)分析級別確定滿足分析精度要求的統(tǒng)一的格網(wǎng)分辨率,利于柵格中間層的構(gòu)建和應(yīng)用。

    (4) 分析尺度。分析尺度主要是在計(jì)算地形因子和地形特征提取時(shí),分析窗口大小的尺度,分析窗口一般分為3×3、5×5、7×7等矩形分析窗口,不同大小的分析窗口往往跟描述地形表面的曲面次數(shù)有關(guān),選擇高次曲面相應(yīng)的分析窗口也應(yīng)更大,應(yīng)根據(jù)具體算法的適應(yīng)性確定分析尺度的大小。

    2 柵格中間層構(gòu)建的適宜尺度分析

    由于尺度分析的復(fù)雜性,在確定柵格中間層構(gòu)建適宜尺度時(shí),應(yīng)把握從整體著眼、方便應(yīng)用,盡量貼近現(xiàn)有標(biāo)準(zhǔn)的原則,在幾個(gè)尺度之間找到恰當(dāng)?shù)钠胶恻c(diǎn),達(dá)到在柵格中間層整體構(gòu)建中現(xiàn)有尺度的矢量數(shù)據(jù)和柵格數(shù)據(jù)能直接參與,滿足大多數(shù)分析任務(wù)精度和速度要求的目的,同時(shí)不應(yīng)過于追求細(xì)節(jié),執(zhí)著于某一特定分析對象或某個(gè)分析算法的尺度細(xì)節(jié)要求。

    范圍尺度、比例尺尺度和分辨率3種尺度之間存在相互關(guān)聯(lián)的關(guān)系,因此在尺度適應(yīng)性研究中,可根據(jù)空間分析任務(wù)級別及各尺度相互之間的關(guān)系分別確定上述3種尺度;分析尺度則應(yīng)根據(jù)具體的柵格數(shù)據(jù)精度、分析算法適應(yīng)性確定。

    3 柵格中間層的尺度適應(yīng)性試驗(yàn)

    本文分別以3×3、5×5、7×7局部擬合窗口生成切線曲率柵格圖和剖面曲率柵格圖,DEM數(shù)據(jù)采用比例尺尺度為1∶5萬時(shí)25 m×25 m的適宜分辨率,通過曲率計(jì)算值的精度、與地形表面形態(tài)匹配情況兩方面比較驗(yàn)證得出適宜的曲率柵格中間層構(gòu)建中的分析尺度。

    3.1 3×3局部擬合窗口

    3×3局部擬合窗口如圖2所示。局部曲面擬合計(jì)算采用以下公式

    (1)

    圖2 3×3局部擬合窗口

    式中

    (2)

    (3)

    (4)

    (5)

    (6)

    移動(dòng)DEM 3×3局部擬合窗口,可順序求出除邊緣一行和一列以外的所有網(wǎng)格點(diǎn)對應(yīng)的多項(xiàng)式系數(shù),代入下式

    可求出試驗(yàn)區(qū)域所有格網(wǎng)點(diǎn)切線曲率的值,將這些切線曲率的值存入新的柵格圖層,得到切線曲率柵格層。同樣將多項(xiàng)式系數(shù)代入下式

    可求出試驗(yàn)區(qū)域所有格網(wǎng)點(diǎn)剖面曲率的值,將這些剖面曲率的值存入新的柵格圖層,得到剖面曲率柵格層。將切線曲率柵格層與DEM柵格暈渲圖對比,剖面曲率柵格層與坡度變率柵格層對比,如圖3所示(切線曲率柵格層與剖面曲率柵格層采用二值化顯示)。

    圖3 3×3局部擬合窗口曲率計(jì)算結(jié)果比較驗(yàn)證

    從圖3可以看出,利用3×3局部擬合窗口生成的切線曲率柵格、剖面曲率柵格與相應(yīng)的DEM柵格暈渲、坡度變率柵格在形態(tài)上相似,但明顯的噪聲點(diǎn)很多,二值化邊界模糊。因此,如果用于地形特征提取,由于大量的噪聲點(diǎn)存在,不能很好地反映較大區(qū)域的地形形態(tài)特征,對于線性地形特征,提取的地形特征線將會界限模糊甚至根本無法體現(xiàn)線性形態(tài),無法用于更高層級的柵格中間層構(gòu)建。

    3.2 5×5局部擬合窗口

    5×5局部擬合窗口如圖4所示。如果需要不提高擬合多項(xiàng)式次數(shù)而增加擬合窗口格網(wǎng)點(diǎn)個(gè)數(shù),就要對擬合多項(xiàng)式進(jìn)行變形,本文利用泰勒級數(shù)公式求出擬合多項(xiàng)式fx±nwn=0,1,2,3,…,n的展開式。通過展開式可以得到fx的一階和二階差分公式,然后采用該點(diǎn)周圍5×5范圍網(wǎng)格點(diǎn)的值對局部二次曲面進(jìn)行最小二乘法擬合,計(jì)算曲面上任意一點(diǎn)的曲率。隨著格網(wǎng)個(gè)數(shù)的增加,展開式就會越來越復(fù)雜,不利于曲率求解,因此下面推導(dǎo)的5×5和7×7局部擬合窗口曲率計(jì)算公式在保證實(shí)際應(yīng)用效果的前提下作了一定量的化簡。

    局部曲面擬合計(jì)算同樣采用式(1),其各個(gè)參數(shù)計(jì)算如下

    r=-z1+16z2-30z3+16z4-z5-z21+16z22-30z23+16z24-z25/96w2+

    -z6+16z7-30z8+16z9-z10-z11+16z12-30z13+16z14-z15-z16+16z17-30z18+16z19-z20/72w2

    7

    t=-z1+16z6-30z11+16z16-z21-z5+16z10-30z15+16z20-z25/96w2+

    -z2+16z7-30z12+16z17-z22-z3+16z8-30z13+16z18-z23-z4+16z9-30z14+16z19-z24/72w2

    8

    s=z1-8z2+8z4-z5/144w2--z6-8z7+8z9-z10/18w2+z16-8z17+18z19-z20/18w2-

    z21-8z22+8z24-z25/144w2

    (9)

    p=z1+z6+z11+z16+z21-8z2+z7+z12+z17+z22+8z4+z9+z14+z19+z24-z5+z10+z15+z20+z25/60w

    (10)

    q=z1+z2+z3+z4+z5-8z6+z7+z8+z9+z10+8z16+z17+z18+z19+z20-z21+z22+z23+z24+z25/60w

    (11)

    圖4 5×5局部擬合窗口

    移動(dòng)DEM 5×5局部擬合窗口,可順序求出除邊緣兩行和兩列點(diǎn)以外的所有網(wǎng)格點(diǎn)對應(yīng)的多項(xiàng)式系數(shù)。與3×3局部擬合窗口相同,將切線曲率柵格層與DEM柵格暈渲對比,剖面曲率柵格層與坡度變率柵格層對比,如圖5所示。

    從圖5可以看出,利用5×5局部擬合窗口生成的切線曲率柵格、剖面曲率柵格與相應(yīng)的DEM柵格暈渲、坡度變率柵格在形態(tài)上相似程度很高,無明顯的噪聲點(diǎn),二值化邊界清晰。說明用5×5局部擬合窗口計(jì)算曲率精度提高,對噪聲的壓制能力增強(qiáng)。因此,用于地形特征提取,能很好地反映較大區(qū)域的地形形態(tài)特征,尤其對于線性地形特征,提取的地形特征線將會線形清晰,可以用于更高層級的柵格中間層構(gòu)建,利用5×5局部擬合窗口計(jì)算曲率是一種適宜的分析尺度。

    圖5 5×5局部擬合窗口曲率計(jì)算結(jié)果比較驗(yàn)證

    3.3 7×7局部擬合窗口

    同樣,移動(dòng)DEM7×7局部擬合窗口,可順序求出除邊緣3行和3列點(diǎn)以外的所有網(wǎng)格點(diǎn)對應(yīng)的多項(xiàng)式系數(shù)。與3×3、5×5局部擬合窗口相同,將切線曲率柵格層與DEM柵格暈渲對比,剖面曲率柵格層與坡度變率柵格層對比,如圖6所示。從圖6可以看出,利用7×7局部擬合窗口生成的切線曲率柵格、剖面曲率柵格噪聲點(diǎn)較5×5局部擬合窗口更少,二值化邊界更加清晰,說明用7×7局部擬合窗口計(jì)算曲率精度更高,對噪聲的壓制能力進(jìn)一步增強(qiáng)。但是與相應(yīng)的DEM柵格暈渲、坡度變率柵格在形態(tài)上相似程度下降,擬合曲面與DEM地形曲面之間匹配性減弱,不利于地形特征的提取。

    圖6 7×7局部擬合窗口曲率計(jì)算結(jié)果比較驗(yàn)證

    因此,太大的局部擬合窗口并不適用于曲率的計(jì)算和曲率柵格層的構(gòu)建,為了使曲率計(jì)算精度較高,又同時(shí)保持?jǐn)M合曲面與DEM地形曲面之間的匹配性,根據(jù)以上試驗(yàn)情況與結(jié)論,推薦使用5×5局部擬合窗口。

    4 結(jié)束語

    本文在曲率柵格中間層適宜分析尺度試驗(yàn)中,得到了與劉學(xué)軍等對于坡度分析尺度研究相似的試驗(yàn)結(jié)果,也印證了坡度分析尺度的可行性,相應(yīng)可推斷在微觀地貌因子的柵格層構(gòu)建中,精度和分析尺度的適應(yīng)性存在一定的共性。

    地物因子的獲取主要是通過矢量地物數(shù)據(jù)柵格化或直接從其他原始數(shù)據(jù)處理結(jié)果中得到,其分析尺度相對簡單。地貌因子的獲取比較復(fù)雜,尤其是微觀地貌因子計(jì)算涉及曲面微分幾何計(jì)算,需要考慮算法的計(jì)算精度、擬合曲面次數(shù)、分析尺度的適宜性等問題,要選擇合適的計(jì)算方法、擬合曲面和分析窗口大小,使計(jì)算結(jié)果達(dá)到分析任務(wù)要求的數(shù)值精度,保證在參與數(shù)值計(jì)算時(shí)得到結(jié)果的準(zhǔn)確性;同時(shí)還要使計(jì)算結(jié)果達(dá)到匹配地表起伏的形態(tài)精度,保證在參與地形特征提取時(shí)得到結(jié)果的形態(tài)相似性。對于宏觀地貌因子的分析尺度,現(xiàn)階段的研究成果中沒有太多可借鑒的結(jié)論。另外,由于地形分析任務(wù)的級別和類型不同,對分析范圍的要求差別也不一樣,不能像微觀地貌因子一樣確定統(tǒng)一的分析尺度,應(yīng)根據(jù)分析任務(wù)具體情況確定相應(yīng)的分析窗口大小。

    [1] 趙敏,林行剛,趙乃國,等. 一種DSM模型下飛行障礙物尺度的研究[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,43(9):1249-1250.

    [2] 吳勇,湯國安,楊昕.小波派生多尺度DEM的精度分析[J].測繪通報(bào),2007(4):38-45.

    [3] 李雙成,蔡云龍.地理尺度轉(zhuǎn)換若干問題的初步探討[J].地理研究,2005,25(1):11-18.

    [4] DUNGAN J I,PERRY J N,DALE M R T,et al.A Balanced View of Scale in Spatial Statistical Analysis[J].Ecography,2002(25):626-640.

    [5] 李霖,應(yīng)申.空間尺度基礎(chǔ)性問題研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2005,30(3):200-202.[6] 湯國安,劉學(xué)軍,房亮,等.DEM 及數(shù)字地形分析中尺度問題研究綜述[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版),2006,31(12):1060-1065.

    [7] 劉學(xué)軍,盧華興,仁政,等.論DEM地形分析中的尺度問題[J].地理研究,2007,26(3):433-442.

    [8] FULLER I C,LARGE A R G,CHARLTON M E,et al.Reachscale Sediment Transfers:An Evaluation of Two Morphological Budgeting Approaches[J].Earth Surface Processes and Landforms,2003,28(8):889-903.

    [9] 張景雄.空間信息的尺度、不確定性與融合[M].武漢:武漢大學(xué)出版社,2008:8-14.

    [10] 申成磊,馬勁松.柵格尺度對DEM分析中地形因子的影響[J].現(xiàn)代測繪,2007,30(6):3-6.

    [11] ANDR S A,F(xiàn)RDRIC C,YANN G,et al.Interpolation of Digital Elevation Models Using AMLE and Related Methods[J].IEEE Transactions on Geoscience and Remote Sensing,2002,40(2):314-325.

    [12] 劉學(xué)軍,龔健雅,周啟鳴,等.基于DEM坡度坡向算法精度的分析研究[J].測繪學(xué)報(bào),2004,33(3):258-263.

    [13] WARREN S D,HOHMANN M G,AUERSWALD K,et al.An Evaluation of Methods to Determine Slope Using Digital Elevation Data[J].Catena,2004(58):215-233.

    [14] 劉學(xué)軍,張平,朱瑩.DEM 坡度計(jì)算的適宜窗口分析[J].測繪學(xué)報(bào),2009,38(3):264-270.

    [15] MAMILLAN R A,SHARY P A.Landfoms and Landform Elements in Geomorphometry[J].Developments in Soil Science,2009(33):233-234.

    [16] 秦承志,呼雪梅.柵格數(shù)字地形分析中的尺度問題研究方法[J].地理研究,2014,33(2):270-283.

    ScaleAdaptabilityofRasterIntermediateLayerinTerrainAnalysis

    ZHANG Lichao1,ZHAO Peng2,ZHANG Hechao3

    (1. Information Center of Henan Administration of Surveying, Mapping and Geoinformation, Zhengzhou 450003,China; 2. Institute of Geospatial Information, Information Engineering University, Zhengzhou 450052,China; 3. Keyuan Surveying and Mapping Center of Henan Province, Zhengzhou 450003,China)

    Raster intermediate layer is fundamental and reusable raster layer for terrain analysis created by the DEM, DSM, vector data, which is obtained through a series of raster calculation. The scale adaptability of raster intermediate layer is studied. The range, scale, resolution, analysis are defined. The principle of range, scale, resolution are determined. Particularly, the paper takes the experiment of scale adaptability of the tangent, curvature raster intermediate layer aiming at the analysis scale, and puts forward the adaptive algorithm and analysis scale. The result shows that raster intermediate layer of terrain foundation factors created from this solution meets the need of numerical precision and morphology precision for terrain analysis.

    raster intermediate layer;scale adaptability;terrain analysis;numerical precision;morphology precision

    P208

    A

    0494-0911(2017)01-0130-04

    張立朝,趙鵬,張合朝.地形分析中柵格中間層構(gòu)建的尺度適應(yīng)性[J].測繪通報(bào),2017(1):130-133.

    10.13474/j.cnki.11-2246.2017.0029.

    2016-03-03

    國家自然科學(xué)基金(41201390);礦山空間信息技術(shù)國家測繪地理信息局重點(diǎn)實(shí)驗(yàn)室開放基金(KLM201411);信息工程大學(xué)地理空間信息學(xué)院自立科研課題(Y1106)

    張立朝(1981—),男,工程師,主要從事地理國情監(jiān)測及地理信息服務(wù)。E-mail:13838205603@163.com

    猜你喜歡
    中間層柵格切線
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    圓錐曲線的切線方程及其推廣的結(jié)論
    切線在手,函數(shù)無憂
    過圓錐曲線上一點(diǎn)作切線的新方法
    鎳基高溫合金TLP擴(kuò)散焊中間層材料研究進(jìn)展
    焊接(2016年8期)2016-02-27 13:05:10
    B含量對IC10合金TLP焊接用中間層材料及接頭組織的影響
    焊接(2016年6期)2016-02-27 13:04:55
    不同剖面形狀的柵格壁對柵格翼氣動(dòng)特性的影響
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計(jì)
    社會中間層建設(shè)與活動(dòng)機(jī)制網(wǎng)研究
    動(dòng)態(tài)柵格劃分的光線追蹤場景繪制
    亚洲av熟女| 国产成+人综合+亚洲专区| 我的老师免费观看完整版| 国产免费av片在线观看野外av| 天天添夜夜摸| 午夜激情福利司机影院| 色综合婷婷激情| 在线观看www视频免费| 午夜a级毛片| 久久中文字幕人妻熟女| 亚洲免费av在线视频| 天堂av国产一区二区熟女人妻 | 精华霜和精华液先用哪个| 国产午夜精品久久久久久| 男人的好看免费观看在线视频 | 波多野结衣高清作品| 国产精品久久久久久人妻精品电影| 国产69精品久久久久777片 | 亚洲成人精品中文字幕电影| 精品国产乱码久久久久久男人| 午夜激情福利司机影院| 国产午夜精品论理片| 亚洲 国产 在线| 十八禁网站免费在线| a级毛片在线看网站| 久久久久性生活片| 变态另类成人亚洲欧美熟女| 在线观看一区二区三区| 母亲3免费完整高清在线观看| xxxwww97欧美| 亚洲成av人片在线播放无| 国产午夜福利久久久久久| 国产精品av视频在线免费观看| av视频在线观看入口| 久久精品国产99精品国产亚洲性色| 好看av亚洲va欧美ⅴa在| 日韩欧美 国产精品| 亚洲精品色激情综合| 亚洲精品久久国产高清桃花| 九色国产91popny在线| 久久国产精品影院| 欧美黄色片欧美黄色片| 日韩大码丰满熟妇| 91老司机精品| 国产亚洲精品久久久久久毛片| 两个人免费观看高清视频| 免费人成视频x8x8入口观看| 国产片内射在线| 成年版毛片免费区| avwww免费| 露出奶头的视频| 欧美色视频一区免费| 久久精品夜夜夜夜夜久久蜜豆 | 女人高潮潮喷娇喘18禁视频| 中文字幕人成人乱码亚洲影| 他把我摸到了高潮在线观看| 亚洲性夜色夜夜综合| cao死你这个sao货| 亚洲国产欧美一区二区综合| 非洲黑人性xxxx精品又粗又长| 欧美中文综合在线视频| 国产又黄又爽又无遮挡在线| 国产精品国产高清国产av| 一进一出抽搐gif免费好疼| 国产精品久久久久久亚洲av鲁大| 国产三级在线视频| 国内揄拍国产精品人妻在线| 老汉色av国产亚洲站长工具| 欧美性猛交╳xxx乱大交人| 午夜免费激情av| 欧美色欧美亚洲另类二区| 亚洲国产精品久久男人天堂| 国产av麻豆久久久久久久| 国产亚洲精品综合一区在线观看 | 国产精品美女特级片免费视频播放器 | 一边摸一边做爽爽视频免费| 国产高清videossex| 一进一出抽搐动态| 国产亚洲av嫩草精品影院| 色精品久久人妻99蜜桃| 又大又爽又粗| 无人区码免费观看不卡| 久久久精品大字幕| 成人特级黄色片久久久久久久| 欧美日韩黄片免| 91大片在线观看| 老司机深夜福利视频在线观看| 亚洲成人久久性| 啦啦啦免费观看视频1| 国产午夜精品久久久久久| 亚洲精品久久国产高清桃花| 制服诱惑二区| 精品欧美一区二区三区在线| 亚洲国产日韩欧美精品在线观看 | 国产黄色小视频在线观看| 久久人妻福利社区极品人妻图片| 国产精品美女特级片免费视频播放器 | 操出白浆在线播放| 久久人妻av系列| 日韩中文字幕欧美一区二区| 久久精品aⅴ一区二区三区四区| 少妇被粗大的猛进出69影院| 长腿黑丝高跟| 无遮挡黄片免费观看| 亚洲国产精品999在线| 99精品欧美一区二区三区四区| 欧美性猛交黑人性爽| 国产精品久久视频播放| 国产成人av教育| 老司机福利观看| 在线国产一区二区在线| 国产99白浆流出| 国产私拍福利视频在线观看| 99久久无色码亚洲精品果冻| 黄色a级毛片大全视频| 老汉色∧v一级毛片| 国产黄a三级三级三级人| 午夜老司机福利片| 国产精品一区二区三区四区免费观看 | 无限看片的www在线观看| 亚洲成人中文字幕在线播放| 两性夫妻黄色片| 看片在线看免费视频| 午夜影院日韩av| 成人午夜高清在线视频| 欧美一级a爱片免费观看看 | 国产人伦9x9x在线观看| 日韩欧美国产一区二区入口| 中亚洲国语对白在线视频| 最近最新中文字幕大全电影3| 三级国产精品欧美在线观看 | 99在线视频只有这里精品首页| 精品久久久久久久久久免费视频| 国产乱人伦免费视频| 2021天堂中文幕一二区在线观| 国产精品免费一区二区三区在线| 成人午夜高清在线视频| 国产真实乱freesex| 女人被狂操c到高潮| 中文字幕最新亚洲高清| 欧美日韩国产亚洲二区| 久久婷婷成人综合色麻豆| 草草在线视频免费看| a级毛片a级免费在线| 日本一二三区视频观看| 午夜福利成人在线免费观看| 成人精品一区二区免费| 午夜激情av网站| 男人舔女人的私密视频| 91麻豆精品激情在线观看国产| 免费观看人在逋| 精品无人区乱码1区二区| 一级毛片女人18水好多| 全区人妻精品视频| 18禁黄网站禁片免费观看直播| 波多野结衣高清作品| 欧美色视频一区免费| 精品久久久久久成人av| 午夜精品一区二区三区免费看| 亚洲专区字幕在线| 人妻夜夜爽99麻豆av| 久久亚洲精品不卡| 91av网站免费观看| 国产亚洲av嫩草精品影院| 人妻久久中文字幕网| 免费观看精品视频网站| 欧美精品亚洲一区二区| 男人舔奶头视频| 亚洲激情在线av| 久久精品国产亚洲av高清一级| 欧美色视频一区免费| 成人18禁高潮啪啪吃奶动态图| 老司机在亚洲福利影院| 亚洲国产看品久久| 亚洲成人久久爱视频| 又黄又爽又免费观看的视频| 色综合站精品国产| 国产又黄又爽又无遮挡在线| 丁香欧美五月| 欧美乱妇无乱码| 好看av亚洲va欧美ⅴa在| 在线永久观看黄色视频| 国产区一区二久久| 久久久久久免费高清国产稀缺| av有码第一页| 后天国语完整版免费观看| 性欧美人与动物交配| 日韩欧美免费精品| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 国产探花在线观看一区二区| 精品一区二区三区av网在线观看| 日日爽夜夜爽网站| 欧美乱码精品一区二区三区| 小说图片视频综合网站| 欧美 亚洲 国产 日韩一| 亚洲人与动物交配视频| 久久久久国产一级毛片高清牌| 亚洲精品国产一区二区精华液| 久久久久久久久中文| 亚洲电影在线观看av| 日韩精品免费视频一区二区三区| 听说在线观看完整版免费高清| bbb黄色大片| 50天的宝宝边吃奶边哭怎么回事| 国产视频一区二区在线看| 巨乳人妻的诱惑在线观看| 久久香蕉国产精品| 亚洲精品久久国产高清桃花| 熟女少妇亚洲综合色aaa.| 欧美日韩瑟瑟在线播放| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人中文| 欧美在线黄色| 麻豆av在线久日| 91老司机精品| 91字幕亚洲| 嫩草影院精品99| 婷婷丁香在线五月| 久久精品91蜜桃| a级毛片a级免费在线| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 免费在线观看黄色视频的| 成人欧美大片| 巨乳人妻的诱惑在线观看| 国产成人影院久久av| 91老司机精品| 999久久久国产精品视频| 极品教师在线免费播放| 很黄的视频免费| 亚洲自偷自拍图片 自拍| 久久婷婷人人爽人人干人人爱| tocl精华| 亚洲av片天天在线观看| 国产一区二区在线观看日韩 | 女人被狂操c到高潮| 国产激情欧美一区二区| 老司机靠b影院| 五月玫瑰六月丁香| 丝袜人妻中文字幕| 久久久久精品国产欧美久久久| 国产黄片美女视频| 一进一出好大好爽视频| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 黄色成人免费大全| 国产高清激情床上av| 亚洲精华国产精华精| 天堂影院成人在线观看| 国产精品一区二区精品视频观看| 国产精品久久久人人做人人爽| 毛片女人毛片| 色尼玛亚洲综合影院| 在线观看www视频免费| 免费高清视频大片| 一本一本综合久久| 亚洲在线自拍视频| 在线观看66精品国产| 久久这里只有精品中国| 母亲3免费完整高清在线观看| 日本在线视频免费播放| 国产精品一区二区精品视频观看| 久久这里只有精品中国| 午夜福利在线观看吧| 蜜桃久久精品国产亚洲av| 午夜激情av网站| 日韩欧美一区二区三区在线观看| 精品久久久久久久末码| 国产成人av激情在线播放| 三级国产精品欧美在线观看 | а√天堂www在线а√下载| 国产久久久一区二区三区| 麻豆久久精品国产亚洲av| 久久久水蜜桃国产精品网| 一级片免费观看大全| 欧洲精品卡2卡3卡4卡5卡区| 一级黄色大片毛片| 十八禁网站免费在线| 免费在线观看亚洲国产| 亚洲欧美精品综合一区二区三区| 国产精品,欧美在线| 亚洲成人久久爱视频| 欧美乱色亚洲激情| 亚洲美女视频黄频| 在线观看免费午夜福利视频| 中文资源天堂在线| 免费电影在线观看免费观看| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 波多野结衣高清作品| 99热6这里只有精品| 老司机在亚洲福利影院| 亚洲国产欧洲综合997久久,| 国产成人一区二区三区免费视频网站| 91九色精品人成在线观看| 婷婷精品国产亚洲av在线| 日韩欧美免费精品| 狂野欧美激情性xxxx| 又黄又爽又免费观看的视频| 久久天堂一区二区三区四区| tocl精华| 91九色精品人成在线观看| 中文字幕人成人乱码亚洲影| 桃红色精品国产亚洲av| 九色国产91popny在线| 精品电影一区二区在线| 一夜夜www| 欧美久久黑人一区二区| 亚洲精品在线观看二区| av超薄肉色丝袜交足视频| 最近最新中文字幕大全电影3| 波多野结衣高清无吗| 男女视频在线观看网站免费 | 免费在线观看完整版高清| av片东京热男人的天堂| 91字幕亚洲| 久久久国产成人免费| 老司机福利观看| 我的老师免费观看完整版| 久久久久久久精品吃奶| 色av中文字幕| 真人一进一出gif抽搐免费| 成人永久免费在线观看视频| 久久这里只有精品中国| 热99re8久久精品国产| 成人高潮视频无遮挡免费网站| 在线视频色国产色| 亚洲国产精品成人综合色| 国产成人av激情在线播放| 国产男靠女视频免费网站| 一边摸一边抽搐一进一小说| 麻豆一二三区av精品| av欧美777| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 18禁国产床啪视频网站| 亚洲人与动物交配视频| 精品福利观看| 亚洲av成人av| www.精华液| 欧美一区二区精品小视频在线| 日本a在线网址| 国产久久久一区二区三区| 国产不卡一卡二| 亚洲欧美精品综合久久99| 亚洲美女黄片视频| 国产单亲对白刺激| 国产一区二区在线观看日韩 | 在线免费观看的www视频| 成人18禁在线播放| 午夜免费观看网址| 国产aⅴ精品一区二区三区波| 国产欧美日韩一区二区三| 法律面前人人平等表现在哪些方面| 国产又色又爽无遮挡免费看| 欧美另类亚洲清纯唯美| 亚洲av成人不卡在线观看播放网| 欧美zozozo另类| 一级毛片精品| 亚洲美女黄片视频| 99国产精品一区二区蜜桃av| 亚洲国产精品久久男人天堂| 久久香蕉精品热| 国产精品av久久久久免费| 亚洲欧美激情综合另类| av免费在线观看网站| 淫秽高清视频在线观看| svipshipincom国产片| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 欧美日本视频| 国产精品免费一区二区三区在线| 最新在线观看一区二区三区| 动漫黄色视频在线观看| 一个人免费在线观看电影 | 麻豆一二三区av精品| 少妇被粗大的猛进出69影院| 国产视频内射| 动漫黄色视频在线观看| 欧美在线一区亚洲| 欧美日本视频| 久久精品影院6| 精品久久久久久久久久久久久| 久久国产乱子伦精品免费另类| 制服诱惑二区| 欧美精品亚洲一区二区| 久久久久久久久久黄片| 90打野战视频偷拍视频| 婷婷丁香在线五月| 97人妻精品一区二区三区麻豆| 国产视频内射| 亚洲电影在线观看av| 欧美黑人欧美精品刺激| 精品久久久久久久久久久久久| 亚洲美女黄片视频| 一进一出抽搐动态| 日本成人三级电影网站| 国产高清有码在线观看视频 | 久久精品91蜜桃| 这个男人来自地球电影免费观看| 99国产精品一区二区三区| 黄色视频,在线免费观看| 国产激情偷乱视频一区二区| 亚洲,欧美精品.| 99热这里只有精品一区 | 午夜视频精品福利| a级毛片在线看网站| 好看av亚洲va欧美ⅴa在| 国产欧美日韩精品亚洲av| 久久伊人香网站| 日韩成人在线观看一区二区三区| 日本一区二区免费在线视频| 动漫黄色视频在线观看| 国产黄色小视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 91麻豆精品激情在线观看国产| 国产单亲对白刺激| 久久精品人妻少妇| 日韩高清综合在线| 无限看片的www在线观看| 亚洲中文av在线| 久久热在线av| 国产主播在线观看一区二区| 欧美黄色淫秽网站| 国产一级毛片七仙女欲春2| 精品久久久久久久末码| 日韩精品中文字幕看吧| 少妇裸体淫交视频免费看高清 | 国产午夜福利久久久久久| 日日干狠狠操夜夜爽| 91在线观看av| 亚洲国产欧美人成| 国内久久婷婷六月综合欲色啪| 婷婷精品国产亚洲av在线| 这个男人来自地球电影免费观看| 免费在线观看视频国产中文字幕亚洲| 亚洲国产欧美人成| www.熟女人妻精品国产| 婷婷精品国产亚洲av在线| 99在线视频只有这里精品首页| 黑人巨大精品欧美一区二区mp4| 欧美性猛交╳xxx乱大交人| 国产精品电影一区二区三区| 亚洲片人在线观看| 午夜福利在线观看吧| 国产午夜精品论理片| 亚洲中文av在线| 制服人妻中文乱码| 少妇人妻一区二区三区视频| 午夜激情福利司机影院| 日韩欧美 国产精品| 麻豆av在线久日| 国产免费男女视频| av有码第一页| 伦理电影免费视频| 少妇人妻一区二区三区视频| 麻豆国产av国片精品| 成年版毛片免费区| 欧美日韩乱码在线| 亚洲国产欧美网| 午夜亚洲福利在线播放| 久久午夜综合久久蜜桃| 美女黄网站色视频| 久久伊人香网站| 丁香欧美五月| 亚洲成a人片在线一区二区| 亚洲精品国产一区二区精华液| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 日本熟妇午夜| 嫩草影院精品99| 黄色毛片三级朝国网站| 日韩欧美一区二区三区在线观看| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| √禁漫天堂资源中文www| 欧美黑人精品巨大| 老司机福利观看| 成人三级黄色视频| 夜夜夜夜夜久久久久| 国产精品影院久久| 亚洲人与动物交配视频| 老汉色∧v一级毛片| 18美女黄网站色大片免费观看| 香蕉久久夜色| a级毛片a级免费在线| 亚洲欧美日韩高清在线视频| 18禁国产床啪视频网站| 成熟少妇高潮喷水视频| 精品久久久久久,| 99久久精品国产亚洲精品| 午夜福利高清视频| 最近在线观看免费完整版| 久久精品91蜜桃| 十八禁人妻一区二区| 欧美日韩精品网址| 大型av网站在线播放| 亚洲专区字幕在线| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| 听说在线观看完整版免费高清| 国语自产精品视频在线第100页| 久久九九热精品免费| av天堂在线播放| www.精华液| 国产黄a三级三级三级人| 不卡av一区二区三区| 久久性视频一级片| 精品久久久久久,| 久久热在线av| 我要搜黄色片| 国产午夜福利久久久久久| 久久久久久亚洲精品国产蜜桃av| 丰满人妻一区二区三区视频av | 免费观看人在逋| 一边摸一边做爽爽视频免费| 日本撒尿小便嘘嘘汇集6| 在线永久观看黄色视频| 国产亚洲精品综合一区在线观看 | 国产爱豆传媒在线观看 | 五月玫瑰六月丁香| 免费在线观看视频国产中文字幕亚洲| 给我免费播放毛片高清在线观看| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 在线观看美女被高潮喷水网站 | а√天堂www在线а√下载| 欧美日韩黄片免| 狂野欧美激情性xxxx| 亚洲18禁久久av| 不卡av一区二区三区| 国产av又大| 搡老岳熟女国产| 一级毛片女人18水好多| 99国产精品一区二区蜜桃av| 999久久久国产精品视频| 午夜成年电影在线免费观看| 露出奶头的视频| 亚洲专区国产一区二区| 亚洲欧美激情综合另类| 美女黄网站色视频| 欧美绝顶高潮抽搐喷水| 好男人在线观看高清免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 国内精品久久久久精免费| av视频在线观看入口| 国产爱豆传媒在线观看 | 日韩三级视频一区二区三区| 757午夜福利合集在线观看| 午夜精品久久久久久毛片777| 国产又色又爽无遮挡免费看| 看片在线看免费视频| 中亚洲国语对白在线视频| 午夜福利在线在线| 国产久久久一区二区三区| or卡值多少钱| 亚洲成人免费电影在线观看| 少妇的丰满在线观看| 长腿黑丝高跟| 国产黄片美女视频| 国产日本99.免费观看| 久久中文看片网| 欧美 亚洲 国产 日韩一| 麻豆国产av国片精品| 午夜a级毛片| 亚洲av片天天在线观看| 日本一二三区视频观看| 巨乳人妻的诱惑在线观看| 国产av麻豆久久久久久久| 午夜福利在线观看吧| 2021天堂中文幕一二区在线观| 亚洲在线自拍视频| 亚洲aⅴ乱码一区二区在线播放 | 国产成人啪精品午夜网站| 露出奶头的视频| 亚洲中文字幕日韩| 精品欧美一区二区三区在线| 久久久久久国产a免费观看| 国产精品1区2区在线观看.| e午夜精品久久久久久久| 欧美黄色片欧美黄色片| 欧美激情久久久久久爽电影| 热99re8久久精品国产| 国产精品,欧美在线| 亚洲最大成人中文| 国产成人一区二区三区免费视频网站| 亚洲成人国产一区在线观看| 天天一区二区日本电影三级| 亚洲人成电影免费在线| 午夜免费观看网址| 99国产精品一区二区三区| 可以在线观看的亚洲视频| 91麻豆精品激情在线观看国产| 老司机深夜福利视频在线观看| 国产熟女午夜一区二区三区| 午夜a级毛片| 韩国av一区二区三区四区| 欧美色视频一区免费| 99久久精品热视频| 久久久久免费精品人妻一区二区| 99在线视频只有这里精品首页| 久久精品国产亚洲av香蕉五月| 亚洲一区二区三区不卡视频| 1024手机看黄色片| 少妇被粗大的猛进出69影院| 超碰成人久久| 一进一出抽搐动态| 岛国视频午夜一区免费看| 白带黄色成豆腐渣| 国产日本99.免费观看| 国产av不卡久久|