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

    基于ArcGIS的SRTM缺失數(shù)據(jù)修復(fù)處理方法

    2011-04-18 08:23:30徐良彭光雄沈蔚
    城市勘測(cè) 2011年1期
    關(guān)鍵詞:空洞標(biāo)準(zhǔn)差插值

    徐良,彭光雄,沈蔚

    (1.常州市測(cè)繪院,江蘇常州 213000; 2.中南大學(xué)地學(xué)與環(huán)境工程學(xué)院,湖南長(zhǎng)沙 410083;3.大洋生物資源開發(fā)與利用上海市高校重點(diǎn)實(shí)驗(yàn)室,上海海洋大學(xué),上海 21306;4.有色金屬成礦預(yù)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,湖南長(zhǎng)沙 410083)

    企業(yè)風(fēng)采

    基于ArcGIS的SRTM缺失數(shù)據(jù)修復(fù)處理方法

    徐良1?,彭光雄2,4,沈蔚3

    (1.常州市測(cè)繪院,江蘇常州 213000; 2.中南大學(xué)地學(xué)與環(huán)境工程學(xué)院,湖南長(zhǎng)沙 410083;3.大洋生物資源開發(fā)與利用上海市高校重點(diǎn)實(shí)驗(yàn)室,上海海洋大學(xué),上海 21306;4.有色金屬成礦預(yù)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室,湖南長(zhǎng)沙 410083)

    利用柵矢轉(zhuǎn)換和掩膜技術(shù)及ArcGIS的空間插值模塊,構(gòu)建了一種SRTM缺失數(shù)據(jù)的修復(fù)處理方法。通過(guò)不同條件下的多種樣本試驗(yàn),探討了外界因素對(duì)插值精度和運(yùn)算效率的影響。結(jié)果表明樣條函數(shù)法是適宜于SRTM缺失數(shù)據(jù)修復(fù)的最佳插值方法。當(dāng)人工模擬的空洞窗口大小為4×4像元時(shí),其插值的誤差平均值和標(biāo)準(zhǔn)差分別為0.90%和0.75%。隨著地形起伏程度的增加,插值誤差的平均值和標(biāo)準(zhǔn)差也隨之加大。同時(shí)插值的精度會(huì)隨著空洞窗口的增大而降低,窗口大小從4×4像元增大至12×12像元時(shí),誤差平均值從0.92%增加至3.84%,標(biāo)準(zhǔn)差則從0.75%增加至3.26%。試驗(yàn)表明6×6像元的窗口大小是保障插值誤差在SRTM數(shù)據(jù)精度范圍內(nèi)的最大窗口。

    SRTM;數(shù)據(jù)修復(fù);空間插值;ArcGIS

    1 引 言

    美國(guó)等航天飛機(jī)雷達(dá)地形測(cè)量(Shuttle Radar Topography Mission,SRTM)的數(shù)字高程數(shù)據(jù)是地理空間建模的重要基礎(chǔ)數(shù)據(jù),在地學(xué)研究中得到廣泛的應(yīng)用[1~5]。然而由于雷達(dá)信號(hào)受到干擾、或者發(fā)生鏡面反射、雷達(dá)陰影或回波滯后等情況導(dǎo)致SRTM高程數(shù)據(jù)在水域、高山區(qū)和峽谷地區(qū)存在有小塊的數(shù)據(jù)空缺點(diǎn)、空白區(qū)。國(guó)內(nèi)外對(duì)SRTM數(shù)據(jù)空洞的填補(bǔ)方法展開了廣泛的研究[6~10]。主要的方法有以下3種類型:一是利用相應(yīng)的其他高分辨率數(shù)字高程數(shù)據(jù)直接填充;二是利用空間插值方法進(jìn)行填充[6,7];三是利用其他輔助數(shù)據(jù)源進(jìn)行數(shù)據(jù)融合后填充[8~10]。其中空間插值方法仍是當(dāng)前研究的主流,主要包括幾何法、統(tǒng)計(jì)法、空間統(tǒng)計(jì)法、函數(shù)法、隨機(jī)模擬法、物理模型模擬法和綜合法等[11]。還有的學(xué)者將分形和形態(tài)學(xué)濾波等方法應(yīng)用到地形數(shù)據(jù)內(nèi)插中來(lái)[12~15]。修復(fù)SRTM缺失數(shù)據(jù)的內(nèi)插方法多種多樣,研究者針對(duì)某一種特定的修復(fù)方法,往往只能自行編寫程序,方能實(shí)現(xiàn)SRTM缺失數(shù)據(jù)的修復(fù)。即使這種修復(fù)方法本身是優(yōu)秀的,卻需要廣大的應(yīng)用人員去編寫繁雜的軟件程序來(lái)實(shí)現(xiàn)其功能,則其實(shí)際應(yīng)用的可操作性將大大降低。因此如果能基于一種主流或通用的軟件工具,通過(guò)某種有效的流程計(jì)算即可實(shí)現(xiàn)上述目的,則對(duì)于地學(xué)研究人員更加廣泛高效地利用SRTM數(shù)據(jù)具有重要的意義。

    本文以主流的GIS平臺(tái)ArcGIS 9.3為工具軟件,通過(guò)構(gòu)建合理有效的SRTM缺失數(shù)據(jù)修復(fù)處理流程,基于其自身的插值工具,在比較多種插值方法計(jì)算性能的基礎(chǔ)上,優(yōu)選出了相對(duì)最為優(yōu)良的插值方法,為SRTM數(shù)據(jù)的廣大用戶填充修復(fù)缺失數(shù)據(jù)提供了一種方便實(shí)用而有效的解決方案。

    2 SRTM缺失數(shù)據(jù)的處理方法

    2.1 數(shù)據(jù)來(lái)源

    SRTM項(xiàng)目于2000年由美國(guó)、意大利和德國(guó)等多家國(guó)際機(jī)構(gòu)共同發(fā)起承擔(dān),其任務(wù)是用于獲取全球范圍高精度的數(shù)字地形。SRTM產(chǎn)品包括3種不同分辨率的數(shù)字高程模型(Digital Elevation Model,DEM)數(shù)據(jù):SRTM1數(shù)據(jù)覆蓋范圍僅僅包括美國(guó)大陸,其空間分辨率為1弧秒(30 m);SRTM3數(shù)據(jù)覆蓋全球,空間分辨率為3弧秒(90 m);SRTM30數(shù)據(jù)同樣覆蓋全球,空間分辨率是30弧秒(1 000 m)[16]。針對(duì)數(shù)據(jù)缺失問(wèn)題,F(xiàn)illed finished版的SRTM產(chǎn)品雖然已經(jīng)對(duì)少于16個(gè)連續(xù)數(shù)據(jù)點(diǎn)的缺失區(qū)域進(jìn)行了填補(bǔ),但是其余面積較大的數(shù)據(jù)缺失區(qū)域仍沒有進(jìn)行處理。因此SRTM的數(shù)據(jù)缺失現(xiàn)象依然存在(如圖1a),必須進(jìn)行修復(fù)處理才能進(jìn)行實(shí)際應(yīng)用。本文從GLCF網(wǎng)站下載了中國(guó)境內(nèi)的一幅SRTM3數(shù)據(jù)作為試驗(yàn)數(shù)據(jù),產(chǎn)品編號(hào)為SRTM_ffB03_p124r039[17],空間分辨率為90m,標(biāo)稱絕對(duì)高程精度±16 m,絕對(duì)平面精度±20 m[16]。

    2.2 基于ArcGIS的空間內(nèi)插算法比較

    ArcGIS是當(dāng)前主流的GIS軟件平臺(tái),在地學(xué)分析中得到廣泛的應(yīng)用。反距離加權(quán)方法(Inverse Distanceto a Power,IDW)、克里金法(Kriging)、最近距離法(Nearest Neighbor,NN)、樣條函數(shù)法(Spline)等是ArcGIS軟件中幾種較為經(jīng)典的空間插值方法。對(duì)于不同的空間插值方法而言,沒有絕對(duì)最優(yōu)的空間內(nèi)插方法,只有特定條件下的最優(yōu)方法[11]。根據(jù)ArcGIS空間插值模塊的自身特點(diǎn)和SRTM數(shù)據(jù)的內(nèi)在特征,通過(guò)比較分析和檢驗(yàn)多種不同條件下的插值試驗(yàn)結(jié)果,最后選擇出適宜于SRTM缺失數(shù)據(jù)修復(fù)的最優(yōu)空間內(nèi)插方法。

    選擇SRTM數(shù)據(jù)中地貌類型豐富且沒有數(shù)據(jù)空洞的區(qū)域切割后作為測(cè)試數(shù)據(jù)。在測(cè)試數(shù)據(jù)中選擇不同地貌類型所在位置區(qū)域設(shè)置16個(gè)4×4像元窗口大小的檢測(cè)樣本(如圖1b);在同一地貌區(qū)域設(shè)置1×1像元至12×12像元共12個(gè)不同窗口大小的檢測(cè)樣本(如圖1c)。

    圖1 SRTM數(shù)據(jù)空洞及測(cè)試樣本分布示意圖

    插值方法的性能特點(diǎn)和試驗(yàn)參數(shù) 表1

    以ArcGIS 9.3為插值工具軟件,利用IDW、Kriging、NN和Spline等4種插值方法對(duì)4×4像元窗口大小的16個(gè)樣本和不同窗口大小的12個(gè)樣本進(jìn)行插值比較分析。插值方法的性能特點(diǎn)[11]和所采用的插值試驗(yàn)參數(shù)見表1所示。通過(guò)這組插值試驗(yàn),比較分析了相同條件下不同插值方法對(duì)插值結(jié)果的影響、地形起伏大小對(duì)插值結(jié)果的影響、缺失數(shù)據(jù)的面積大小對(duì)插值結(jié)果的影響以及不同插值方法的計(jì)算效率等4種

    2.3 SRTM缺失數(shù)據(jù)修復(fù)流程與處理

    很多文獻(xiàn)中提到的基于柵格數(shù)據(jù)的空間插值一般只描述了插值方法的工作原理,對(duì)于具體的插值過(guò)程往往很少進(jìn)行說(shuō)明[10,18]。而實(shí)際情況中插值過(guò)程的不同,會(huì)導(dǎo)致插值結(jié)果的顯著差異,特別是利用現(xiàn)有的插值工具軟件進(jìn)行柵格數(shù)據(jù)的空間插值處理。以基于ArcGIS 9.3的SRTM空洞數(shù)據(jù)插值為例說(shuō)明如下。ArcGIS 9.3的Geostatistical Analyst模塊既可以對(duì)柵格數(shù)據(jù)也可以對(duì)矢量點(diǎn)數(shù)據(jù)進(jìn)行空間插值,而Interploation模塊則僅能對(duì)矢量點(diǎn)數(shù)據(jù)進(jìn)行空間插值。用Geostatistical Analyst模塊對(duì)帶有人工空洞的SRTM柵格數(shù)據(jù)直接進(jìn)行插值,所得的結(jié)果如圖2a所示,可見空洞區(qū)域并沒有得到合理的填充。由于空洞區(qū)域的像評(píng)價(jià)指標(biāo),探討了外界因素對(duì)插值的影響,從而優(yōu)選出適宜于SRTM缺失數(shù)據(jù)修復(fù)的最佳插值方法。試驗(yàn)選取完好的SRTM數(shù)據(jù)在各個(gè)樣本區(qū)域通過(guò)人工填充負(fù)值(-32 768,跟真實(shí)的缺失數(shù)據(jù)值一致)的方式模擬SRTM缺失數(shù)據(jù)的空洞,并由此獲得用于計(jì)算插值精度的真實(shí)值。插值結(jié)果的精度用誤差百分比來(lái)衡量,其計(jì)算方法如公式1所示。元本身是有值的(-32 768),只不過(guò)是缺失的無(wú)效值而已,ArcGIS的插值模塊對(duì)此無(wú)法進(jìn)行判別而仍當(dāng)作正常值進(jìn)行處理,從而導(dǎo)致空洞區(qū)域的像元仍然保留插值前的值。因此利用ArcGIS軟件進(jìn)行SRTM柵格數(shù)據(jù)直接插值來(lái)填充空洞區(qū)域的方法是行不通的。

    解決上述問(wèn)題的思路是對(duì)SRTM空洞區(qū)的缺失數(shù)據(jù)先進(jìn)行一個(gè)預(yù)處理,使得ArcGIS僅對(duì)有效數(shù)據(jù)進(jìn)行插值而對(duì)缺失數(shù)據(jù)不插值。首先,利用Conversion Tools模塊的Raster to Point工具將SRTM柵格數(shù)據(jù)轉(zhuǎn)為Shape格式的點(diǎn)矢量文件;其次,通過(guò)Shape點(diǎn)矢量文件的屬性表對(duì)所有值為-32 768(缺失值)的點(diǎn)進(jìn)行選取并刪除;最后,對(duì)不包含缺失值的Shape點(diǎn)矢量文件進(jìn)行插值。

    圖2b是沒有去除空值點(diǎn)的矢量點(diǎn)插值的結(jié)果,其實(shí)質(zhì)是跟柵格數(shù)據(jù)直接插值是相同的,從而所得到結(jié)果也是不合理的。這是由于每個(gè)點(diǎn)位上該點(diǎn)自身的權(quán)重是最大的,兩點(diǎn)之間的距離為90 m(由SRTM數(shù)據(jù)的空間分辨率決定),當(dāng)插值輸出像元大小為90 m時(shí),則插值的結(jié)果其實(shí)就是將每個(gè)點(diǎn)自身的值(-32 768)賦給該點(diǎn),從而在圖2b中形成一條值為-32 768的封閉等值線。圖2c是去除空值點(diǎn)后的矢量點(diǎn)插值的結(jié)果,由于消除了缺失值點(diǎn)的影響,則在圖2c中形成多條相互連通的等值線。

    基于ArcGIS的SRTM缺失數(shù)據(jù)修復(fù)方法的總體流程如圖3所示。利用Conversion Tools模塊的Raster to Point工具實(shí)現(xiàn)柵矢轉(zhuǎn)換;利用Attributes屬性表剔除缺失點(diǎn)數(shù)據(jù);利用Interploation模塊將矢量點(diǎn)數(shù)據(jù)插值生產(chǎn)柵格數(shù)據(jù);利用Reclass模塊的Reclassify工具以-32 768為閾值提取空洞掩膜Mask(0,1)和Mask(1,0),如圖4b所示;利用Math模塊的Times(乘)和Plus (加)運(yùn)算進(jìn)行掩膜處理,如式(2)所示;最終得到SRTM的缺失數(shù)據(jù)的修復(fù)結(jié)果,如圖4c所示。

    圖2 柵格數(shù)據(jù)和Shape點(diǎn)矢量數(shù)據(jù)插值結(jié)果比較

    圖3 基于ArcGIS的SRTM缺失數(shù)據(jù)修復(fù)流程

    圖4 SRTM缺失數(shù)據(jù)修復(fù)前后比較

    式中,SRTM是修復(fù)后的數(shù)據(jù),SRTM1是帶有空洞的數(shù)據(jù),SRTM2是插值后的數(shù)據(jù),Mask(0,1)是空洞賦值為0的掩膜,Mask(0,1)是空洞賦值為1的掩膜。

    3 修復(fù)處理結(jié)果分析與方法評(píng)價(jià)

    3.1 修復(fù)處理結(jié)果分析

    對(duì)不同地形條件下的16個(gè)4×4像元窗口大小的人工模擬空洞經(jīng)IDW、Kriging、NN和Spline等4種方法的插試驗(yàn)值進(jìn)行了統(tǒng)計(jì),結(jié)果如表2和圖5所示。圖5a顯示的是16個(gè)樣本四種不同方法插值精度的比較曲線,誤差百分比的平均值由高到低依次為IDW、 NN、Krging和Spline,分別為1.32%、1.13%、1.07%和0.90%,各自的標(biāo)準(zhǔn)差則分別為0.90%、0.81%、0.62%和0.75%。誤差的平均值反映出各種插值方法所得結(jié)果與真實(shí)值的接近程度,標(biāo)準(zhǔn)差則體現(xiàn)了各種方法的穩(wěn)定性。綜合比較分析誤差的平均值和標(biāo)準(zhǔn)差的情況,可以看出Spline是這4種插值方法中最佳的一種。圖5b顯示的是插值精度隨地形起伏變化的情況,Spline法插值的結(jié)果表明隨著地形起伏程度的增加,其誤差平均值和標(biāo)準(zhǔn)差也隨之加大。誤差平均值由0.56%增加至1.21%,標(biāo)準(zhǔn)差則由0.41%增加至1.08%,說(shuō)明地形起伏越小插值的精度越高,穩(wěn)定性也越好。

    空值區(qū)域?yàn)?×4窗口的插值統(tǒng)計(jì)結(jié)果 表2

    圖5 4×4窗口時(shí)的插值統(tǒng)計(jì)曲線

    圖6 不同方法插值運(yùn)算效率曲線

    圖7 插值精度隨窗口大小的變化曲線

    插值精度是反映插值方法性能的重要指標(biāo)之一,而另一重要指標(biāo)則是運(yùn)算效率,特別是對(duì)于大規(guī)模的插值計(jì)算尤其重要。通過(guò)截取不同像元總量的SRTM數(shù)據(jù),基于同臺(tái)電腦,利用上述4種方法進(jìn)行插值試驗(yàn),得到如圖6所示的統(tǒng)計(jì)結(jié)果??梢娤嗤闆r下插值的運(yùn)算效率由低至高依次為Krging、Spline、IDW、和NN。其中Krging的運(yùn)算量相對(duì)而言非常大、效率十分低下,而Spline和IDW相差不大。結(jié)合插值的精度和運(yùn)算效率,則可進(jìn)一步判斷出Spline是最適宜于STRM數(shù)據(jù)的插值方法。

    在相同地形條件下,人工模擬出1×1像元至12× 12像元共12個(gè)不同窗口SRTM數(shù)據(jù)空洞(如圖1c),并以Spline方法為例測(cè)試空洞面積大小對(duì)插值精度的影響,統(tǒng)計(jì)結(jié)果如表3和圖7所示??梢?,插值的誤差百分比隨著窗口的增大而逐漸增加,且以誤差百分比的最大值增加趨勢(shì)尤其顯著。誤差百分比的平均值和標(biāo)準(zhǔn)差的增加趨勢(shì)基本一致,隨著窗口增大在插值精度降低的同時(shí)穩(wěn)定性也隨之變小。例如當(dāng)窗口大小從4×4增大至12×12時(shí),誤差百分比的平均值從0.92%增加至3.84%,標(biāo)準(zhǔn)差則從0.75%增加至3.26%。測(cè)試窗口所在區(qū)域的平均海拔為1 102 m,則該地區(qū)4×4窗口插值的平均精度約為±10.2 m,12×12窗口插值的平均精度約為±42.3 m。按SRTM3數(shù)據(jù)的標(biāo)稱絕對(duì)高程精度±16 m來(lái)計(jì)算,則該地區(qū)相應(yīng)的誤差百分比約為1.6%,對(duì)應(yīng)的窗口大小在6×6左右。超出6×6窗口大小插值得到的結(jié)果將會(huì)引起較大的失真。當(dāng)然,這種用海拔和數(shù)據(jù)的標(biāo)稱絕對(duì)高程精度所衡量的適宜窗口大小,只是相對(duì)意義上的衡量方法,它會(huì)隨著海拔高度、地形起伏程度以及空洞區(qū)域的形狀和位置等因素的變化而會(huì)有所差異。

    不同窗口大小的Spline插值統(tǒng)計(jì)結(jié)果 表3

    3.2 基于ArcGIS的修復(fù)方法評(píng)價(jià)

    基于ArcGIS的修復(fù)方法,僅僅利用ArcGIS軟件本身而不需要編寫內(nèi)插程序就能對(duì)SRTM的缺失數(shù)據(jù)進(jìn)行修復(fù)處理,且操作簡(jiǎn)單快捷,修復(fù)精度高,具有良好的實(shí)用效果和廣闊的應(yīng)用前景。通過(guò)上述插值精度和運(yùn)算效率的比較,表明Spline法是基于ArcGIS修復(fù)SRTM缺失數(shù)據(jù)的最佳空間內(nèi)插算法。該方法可修復(fù)缺失數(shù)據(jù)窗口為任意大小的SRTM數(shù)據(jù),當(dāng)且僅當(dāng)SRTM的缺失數(shù)據(jù)的窗口小于或等于6×6像元時(shí),其誤差百分比可達(dá)1.4%左右,基本可達(dá)到SRTM數(shù)據(jù)±16 m的標(biāo)稱絕對(duì)高程精度。

    4 結(jié)論與討論

    利用ArcGIS軟件平臺(tái),通過(guò)構(gòu)建合理的處理流程,探討了一種SRTM缺失數(shù)據(jù)的修復(fù)方法。主要得到以下結(jié)論:

    (1)利用ArcGIS軟件進(jìn)行SRTM缺失數(shù)據(jù)的空間插值時(shí),必須先要進(jìn)行柵矢轉(zhuǎn)換和空值點(diǎn)的剔除處理,才能保證插值結(jié)果的有效性;

    (2)通過(guò)對(duì)ArcGIS軟件現(xiàn)有的幾種插值方法的測(cè)試比較,發(fā)現(xiàn)樣條函數(shù)法的插值結(jié)果最好,運(yùn)算效率也較為理想,是最適宜于SRTM缺失數(shù)據(jù)修復(fù)的插值方法;

    (3)地形起伏對(duì)插值的精度影響較大,隨著地形起伏程度的增加,插值誤差的平均值和標(biāo)準(zhǔn)差也隨之加大;

    (4)SRTM缺失數(shù)據(jù)區(qū)域的面積大小和形狀也對(duì)插值的精度也存在很大的影響,插值的精度會(huì)隨著空洞窗口的增大而降低。初步的測(cè)試結(jié)果表明空洞窗口大小在6×6像元以下時(shí)基本能保障插值填補(bǔ)數(shù)據(jù)的精度,而超過(guò)此窗口大小時(shí)將會(huì)產(chǎn)生較大的失真;

    (5)測(cè)試樣本的數(shù)量和分布區(qū)域有待更進(jìn)一步的提高和完善,6×6像元的精度保障窗口大小只是一個(gè)初步的判斷,實(shí)際結(jié)果會(huì)因?yàn)榈匦纹鸱潭纫约翱斩磪^(qū)域的形狀和位置等因素的變化而會(huì)有所差異。

    [1] 游松財(cái),孫朝陽(yáng).中國(guó)流域的SRTM30數(shù)據(jù)提取與計(jì)算[J].地球信息科學(xué)學(xué)報(bào),2009,11(02):189~194

    [2] 龍恩,程維明,周成虎等.基于Srtm-DEM與遙感的長(zhǎng)白山基本地貌類型提取方法[J].山地學(xué)報(bào),2009,25(5): 557~565

    [3] 彭光雄,宮阿都,崔偉宏等.多時(shí)相影像的典型區(qū)農(nóng)作物識(shí)別分類方法對(duì)比研究.地球信息科學(xué)學(xué)報(bào),2009,11 (02):225~230

    [4] 趙尚民,程維明,柴慧霞等.基于遙感與SRTM的青藏高原冰緣地貌信息提取方法—以1∶100萬(wàn)標(biāo)準(zhǔn)分幅拉薩幅(H46)為例[J].地理研究,2007,26(6):1175~1187

    [5] Oky D A P,Ayako E,Ryuzo Yo,et al.Supervised landform classification of North-east Honshu from DEM-derived thematic maps[J].Geomorphology,2006,78:373~386

    [6] 游松財(cái),孫朝陽(yáng).中國(guó)區(qū)域SRTM90m數(shù)字高程數(shù)據(jù)空值區(qū)域的填補(bǔ)方法比較[J].地理科學(xué)進(jìn)展,2005,24(6): 88 92

    [7] 王濤,毋河海.SRTM高程數(shù)據(jù)中空缺單元的內(nèi)插填補(bǔ)[J].測(cè)繪科學(xué),2006,31(3):76 77

    [8] Kaab A.Combination of SRTM3 and Repeat ASTER Data for Deriving Alpine Glacier Flow velocities in the Bhutan Himalaya [J].Remote Sensing of Environment,2005(94):463~671

    [9] 凌峰,王乘,張秋文.基于ASTER數(shù)據(jù)和空間誤差分析的SRTM無(wú)效區(qū)域填充[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,34(12),108~111

    [10] 闞璦珂,朱利東,張瑞軍等.基于數(shù)據(jù)融合的SRTM數(shù)據(jù)空洞填補(bǔ)方法[J].地理空間信息,2007,5(3):62~64

    [11] 李新,程國(guó)棟,盧玲.空間內(nèi)插方法比較[J].地球科學(xué)進(jìn)展,2000,15(3):260~265

    [12] 張山山.分形方法在地形數(shù)據(jù)內(nèi)插中的應(yīng)用.西南交通大學(xué)學(xué)報(bào),2000,35(2):141~14

    [13] Weidner U.,F(xiàn)orstner W.Towards Automatic Building Extraction from High Resolution Digital Elevation Models[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 1995,50,38~49

    [14] Kilian J.,Haala N.,Englich M.Capturing and evaluation of airborne laser scanner data[J].International Archives of Photogrammetry and Remote Sensing,Vienna,1996,32,(B3):383~388

    [15] Hug C.,Wehr A.Detecting and Identifying Topographic Objects in Imaging Laser Altimeter Data[J].International Archives of Photogrammetry and Remote Sensing,1997,32:19~26

    [16] USGS(2004),Shuttle Radar Topography Mission,1 Arc Second scene SRTM_u03_n008e004,Unfilled Unfinished 2. 0,Global Land Cover Facility,University of Maryland,College Park,Maryland,F(xiàn)ebruary 2000

    [17] http://glcf.umiacs.umd.edu/data/srtm/index.shtml.

    [18] 黃文捷.利用插值算法填補(bǔ)SRTM3 DEM數(shù)據(jù)空洞的比較分析.江西測(cè)繪.2007,69(3):25~28

    Method Study on SRTM Void Data Recovery Based on ArcGIS

    Xu Liang1,Peng GuangXiong2,4,Shen Wei3
    (1.School of Geoscience and Environmental Engineering,Central South University,Changsha 410083,China; 2.The Key Laboratory of Shanghai Education Commission for Oceanic Fisheries Resources Exploitation,Shanghai Ocean University,Shanghai 201306,China; 3.Key Laboratory of Metallogenic Prediction of Nonferrous Metals,Ministry of Education,Changsha 410083,China)

    The digital elevation data of Shuttle Radar Topography Mission(SRTM)is an important foundation data for geospatial modeling.However the void data of SRTM caused by the deficiencies of radar produce great inconvenience for its application.An approach and process flow of SRTM void data recovery is constructed using raster to vector conversion and mask processing technology based on spatial interpolation in ArcGIS software.To ensure the validity of interpolation results,the most important step of this process flow is deleting points of void data.The interpolation accuracy and operations efficiency of Inverse Distance to a Power(IDW),Kriging,Nearest Neighbor(NN)and Spline are compared and analyzed based on a variety of test samples under different conditions.16 test samples of 4×4 pixel window size distribute on different landforms regions,while12 test samples of different pixel window size from 1×1 to12×12 locate on the same landforms region.The hollows of SRTM void data are simulated artificially through fill some negative value to the SRTM data in the areas of test samples.With this way to obtain the true value used to estimate interpolation accuracy.The mean of interpolation error is in decreasing order is IDW,NN,Krging and Spline,whose value is 1.32%,1.13%,1.07%and 0. 90%respectively.Their standard deviation value are 0.90%、0.81%、0.62%and 75%respectively.Results show that Spline interpolation method is the optimal one applicable to recover void data of SRTM.When the artificial simulation window size of empty area is 4×4 pixels,the mean and standard deviation of interpolation error is 0.90%and 0.75%respectively.With the levels of undulating topography increase,the mean and standard deviation of interpolation error will increase accordingly.On the other hand,interpolation accuracy will descend with the accretion of window size.When the window size of empty area increase form 4×4 pixels to12×12 pixels,the mean of interpolation error increased from 0.92% to 3.84%,standard deviation increased from 0.75%to 3.26%.6×6 pixels is the largest window size of empty area which can restrict the interpolation error within the data precision of SRTM.Experimental results show that the approach proposed in this paper is a good way to recover SRTM void data only using ArcGIS software.

    SRTM;data recovery;spatial interpolation;ArcGIS

    1672-8262(2011)01-5-06

    N945.1

    A

    2010—12—21

    徐良(1977—),男,工程師,主要從事城市測(cè)繪與GIS、城市三維等方面的工作。

    中國(guó)博士后科學(xué)基金(20080430586);大洋生物資源開發(fā)與利用重點(diǎn)實(shí)驗(yàn)室開放基金(B-8208-07-0001-1)。

    猜你喜歡
    空洞標(biāo)準(zhǔn)差插值
    用Pro-Kin Line平衡反饋訓(xùn)練儀對(duì)早期帕金森病患者進(jìn)行治療對(duì)其動(dòng)態(tài)平衡功能的影響
    基于Sinc插值與相關(guān)譜的縱橫波速度比掃描方法
    空洞的眼神
    一種改進(jìn)FFT多譜線插值諧波分析方法
    基于四項(xiàng)最低旁瓣Nuttall窗的插值FFT諧波分析
    用事實(shí)說(shuō)話勝過(guò)空洞的說(shuō)教——以教育類報(bào)道為例
    新聞傳播(2015年20期)2015-07-18 11:06:46
    對(duì)于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
    Blackman-Harris窗的插值FFT諧波分析與應(yīng)用
    臭氧層空洞也是幫兇
    醫(yī)學(xué)科技論文中有效數(shù)字的確定
    久久精品国产亚洲网站| 蜜桃亚洲精品一区二区三区| 99久久精品一区二区三区| 亚洲欧美日韩东京热| 国产精品国产三级专区第一集| 日本wwww免费看| 成人特级av手机在线观看| 熟女人妻精品中文字幕| 波多野结衣巨乳人妻| 九九在线视频观看精品| 亚洲av男天堂| 少妇裸体淫交视频免费看高清| 人妻少妇偷人精品九色| 午夜视频国产福利| 国产精品国产三级专区第一集| 又粗又硬又长又爽又黄的视频| 亚洲精品影视一区二区三区av| 秋霞在线观看毛片| 中文字幕av在线有码专区| 午夜福利在线观看免费完整高清在| www日本黄色视频网| АⅤ资源中文在线天堂| 成人性生交大片免费视频hd| 99热这里只有是精品50| 免费观看在线日韩| 爱豆传媒免费全集在线观看| 久久韩国三级中文字幕| 级片在线观看| 国产成人精品久久久久久| 色综合亚洲欧美另类图片| 久久99热这里只频精品6学生 | 久久久久久久久久久丰满| 中文资源天堂在线| 又粗又爽又猛毛片免费看| 亚洲自偷自拍三级| 欧美xxxx性猛交bbbb| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站| 日韩欧美三级三区| 亚洲精品国产av成人精品| 欧美另类亚洲清纯唯美| 成人高潮视频无遮挡免费网站| 婷婷色麻豆天堂久久 | 看免费成人av毛片| 亚洲精品成人久久久久久| 一级二级三级毛片免费看| 人妻少妇偷人精品九色| 天堂网av新在线| 国产熟女欧美一区二区| 午夜a级毛片| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 联通29元200g的流量卡| 国产老妇女一区| 男人舔女人下体高潮全视频| 偷拍熟女少妇极品色| 亚洲成人久久爱视频| 亚洲在线观看片| 在线观看66精品国产| 少妇人妻一区二区三区视频| 日韩强制内射视频| 国产真实乱freesex| 国产在视频线在精品| 国产毛片a区久久久久| 亚洲欧美精品专区久久| 国产精品日韩av在线免费观看| 亚洲欧美清纯卡通| 高清在线视频一区二区三区 | 嫩草影院精品99| 最近2019中文字幕mv第一页| 精品国内亚洲2022精品成人| 1000部很黄的大片| 国产成人福利小说| 亚洲国产色片| 欧美97在线视频| 国产免费一级a男人的天堂| 欧美精品国产亚洲| 国产老妇伦熟女老妇高清| 日韩视频在线欧美| 91久久精品电影网| 床上黄色一级片| 亚洲欧美一区二区三区国产| 国产成人午夜福利电影在线观看| 久久久午夜欧美精品| 成人漫画全彩无遮挡| 免费看av在线观看网站| 内射极品少妇av片p| 婷婷色综合大香蕉| 久久6这里有精品| 精品99又大又爽又粗少妇毛片| 国语自产精品视频在线第100页| 婷婷六月久久综合丁香| 亚洲婷婷狠狠爱综合网| 精品久久久久久久久久久久久| av天堂中文字幕网| 国产私拍福利视频在线观看| 午夜福利视频1000在线观看| 日韩视频在线欧美| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| a级毛色黄片| 搞女人的毛片| 欧美日韩精品成人综合77777| 69av精品久久久久久| 高清毛片免费看| 老司机影院毛片| 国产精品国产三级国产av玫瑰| 国产片特级美女逼逼视频| a级一级毛片免费在线观看| 亚洲精品乱码久久久久久按摩| 青春草亚洲视频在线观看| 久久精品久久久久久久性| 两性午夜刺激爽爽歪歪视频在线观看| 91av网一区二区| 一区二区三区四区激情视频| 欧美另类亚洲清纯唯美| 国产av码专区亚洲av| av又黄又爽大尺度在线免费看 | 精品国产露脸久久av麻豆 | av在线天堂中文字幕| 99热全是精品| 日日摸夜夜添夜夜添av毛片| av免费观看日本| 精品久久久久久久末码| 男女那种视频在线观看| 欧美一区二区精品小视频在线| 国产亚洲最大av| 麻豆av噜噜一区二区三区| 少妇熟女aⅴ在线视频| 婷婷六月久久综合丁香| 毛片女人毛片| 看片在线看免费视频| 人人妻人人澡欧美一区二区| 亚洲av熟女| 级片在线观看| 成人毛片60女人毛片免费| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 亚洲av免费在线观看| 亚洲欧美中文字幕日韩二区| 亚洲国产日韩欧美精品在线观看| 婷婷色综合大香蕉| 午夜老司机福利剧场| 国产精品电影一区二区三区| 欧美日本亚洲视频在线播放| 一个人免费在线观看电影| 女的被弄到高潮叫床怎么办| av福利片在线观看| 国产亚洲91精品色在线| 少妇人妻精品综合一区二区| 欧美激情国产日韩精品一区| 热99在线观看视频| 日本欧美国产在线视频| 亚洲三级黄色毛片| 久久99蜜桃精品久久| 毛片一级片免费看久久久久| av在线老鸭窝| 又粗又爽又猛毛片免费看| 91在线精品国自产拍蜜月| 日韩欧美国产在线观看| 国内少妇人妻偷人精品xxx网站| 视频中文字幕在线观看| 在线播放国产精品三级| 中文天堂在线官网| 青青草视频在线视频观看| 亚洲精品乱码久久久v下载方式| 免费观看在线日韩| 亚洲最大成人av| 蜜桃久久精品国产亚洲av| 天堂中文最新版在线下载 | 蜜臀久久99精品久久宅男| 午夜亚洲福利在线播放| videos熟女内射| 又粗又爽又猛毛片免费看| 亚洲欧美日韩卡通动漫| 国产精品麻豆人妻色哟哟久久 | 99九九线精品视频在线观看视频| 蜜桃久久精品国产亚洲av| 精品一区二区三区人妻视频| 舔av片在线| 别揉我奶头 嗯啊视频| 国产女主播在线喷水免费视频网站 | 免费播放大片免费观看视频在线观看 | 直男gayav资源| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 亚洲国产精品国产精品| 日韩人妻高清精品专区| 亚洲欧美日韩卡通动漫| 啦啦啦观看免费观看视频高清| 国产 一区精品| 18禁裸乳无遮挡免费网站照片| 国产又色又爽无遮挡免| 精品久久久久久久久亚洲| 男人的好看免费观看在线视频| 91久久精品国产一区二区成人| 国产色爽女视频免费观看| 国产精品国产三级专区第一集| 性插视频无遮挡在线免费观看| 99热精品在线国产| av国产免费在线观看| 日本一二三区视频观看| 久久亚洲国产成人精品v| 亚洲伊人久久精品综合 | 欧美三级亚洲精品| 欧美一区二区亚洲| 最近最新中文字幕大全电影3| 亚洲av中文av极速乱| 久久韩国三级中文字幕| 欧美潮喷喷水| 在线观看一区二区三区| 国产片特级美女逼逼视频| 看免费成人av毛片| 九色成人免费人妻av| 中文字幕亚洲精品专区| 成人毛片60女人毛片免费| 两个人视频免费观看高清| 国产精品电影一区二区三区| 成人国产麻豆网| 国产美女午夜福利| 三级经典国产精品| 色综合站精品国产| 国产av一区在线观看免费| 国产精品国产三级专区第一集| 一级毛片我不卡| 国产精品不卡视频一区二区| 亚洲综合色惰| 桃色一区二区三区在线观看| 日本黄色视频三级网站网址| 视频中文字幕在线观看| 成人综合一区亚洲| 黄片wwwwww| 国产美女午夜福利| 国产老妇女一区| 一级毛片久久久久久久久女| 成人综合一区亚洲| 少妇丰满av| 午夜激情欧美在线| 欧美区成人在线视频| av在线天堂中文字幕| 国产成人精品一,二区| 91aial.com中文字幕在线观看| 亚洲av免费高清在线观看| 国产成人freesex在线| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 波多野结衣高清无吗| 国产在视频线精品| 欧美+日韩+精品| 亚洲国产欧洲综合997久久,| 亚洲久久久久久中文字幕| 亚洲精品自拍成人| 性色avwww在线观看| 亚洲成av人片在线播放无| 少妇人妻一区二区三区视频| 免费大片18禁| 97人妻精品一区二区三区麻豆| av在线播放精品| 久久精品久久久久久久性| 久久精品国产亚洲av涩爱| 九九久久精品国产亚洲av麻豆| 国产一区二区在线av高清观看| 一本久久精品| 级片在线观看| 大香蕉久久网| 日韩三级伦理在线观看| 水蜜桃什么品种好| 亚洲色图av天堂| 高清av免费在线| 极品教师在线视频| 夜夜看夜夜爽夜夜摸| 在线观看美女被高潮喷水网站| 天堂影院成人在线观看| 成人毛片a级毛片在线播放| 激情 狠狠 欧美| 韩国高清视频一区二区三区| 国产又色又爽无遮挡免| 亚洲av成人精品一二三区| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站 | 中文字幕久久专区| 在线观看一区二区三区| 一级黄片播放器| 身体一侧抽搐| av卡一久久| 亚洲成av人片在线播放无| 中文精品一卡2卡3卡4更新| 国产精品永久免费网站| 日日撸夜夜添| 亚洲av电影不卡..在线观看| 久久99精品国语久久久| 亚洲欧美中文字幕日韩二区| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 国产伦一二天堂av在线观看| 国国产精品蜜臀av免费| 国产大屁股一区二区在线视频| 亚洲精品日韩av片在线观看| 熟女人妻精品中文字幕| 亚洲欧美精品自产自拍| 丰满少妇做爰视频| 99久久九九国产精品国产免费| 舔av片在线| 国产精品野战在线观看| 亚洲av福利一区| 又爽又黄a免费视频| 黄片wwwwww| 麻豆av噜噜一区二区三区| 久久热精品热| 欧美色视频一区免费| 你懂的网址亚洲精品在线观看 | 永久免费av网站大全| 午夜免费激情av| 网址你懂的国产日韩在线| 综合色丁香网| 日韩成人伦理影院| 91久久精品国产一区二区三区| 国产高清三级在线| 99久久精品热视频| 精品久久久久久成人av| 婷婷色综合大香蕉| 女人被狂操c到高潮| 久久欧美精品欧美久久欧美| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 美女cb高潮喷水在线观看| www.av在线官网国产| 午夜亚洲福利在线播放| 中文资源天堂在线| 我要看日韩黄色一级片| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| 中文字幕久久专区| 亚洲国产欧美在线一区| 又黄又爽又刺激的免费视频.| 国产伦精品一区二区三区视频9| 蜜桃亚洲精品一区二区三区| 国产精品无大码| 国产精品一二三区在线看| 日韩av在线免费看完整版不卡| 欧美色视频一区免费| 国产一级毛片七仙女欲春2| 春色校园在线视频观看| 精品国产一区二区三区久久久樱花 | 精品99又大又爽又粗少妇毛片| 亚洲国产精品sss在线观看| 桃色一区二区三区在线观看| 日本五十路高清| 熟女人妻精品中文字幕| 欧美一区二区精品小视频在线| 国产老妇伦熟女老妇高清| 亚洲中文字幕一区二区三区有码在线看| 97热精品久久久久久| 网址你懂的国产日韩在线| 91久久精品国产一区二区成人| 亚洲国产精品合色在线| 最近中文字幕高清免费大全6| 国产免费男女视频| 综合色丁香网| 亚洲在久久综合| 深爱激情五月婷婷| 久久久久久久久久成人| 少妇裸体淫交视频免费看高清| 国产成人一区二区在线| 超碰av人人做人人爽久久| 91狼人影院| 我要搜黄色片| 精品少妇黑人巨大在线播放 | 国产熟女欧美一区二区| 午夜免费激情av| 日本免费在线观看一区| 亚洲欧美清纯卡通| 免费观看在线日韩| 精品欧美国产一区二区三| 麻豆久久精品国产亚洲av| 老司机影院成人| av在线蜜桃| 久久久久网色| 久久久久久九九精品二区国产| 麻豆成人午夜福利视频| 国产乱人视频| 秋霞在线观看毛片| 日本与韩国留学比较| 色视频www国产| 只有这里有精品99| 男人狂女人下面高潮的视频| 五月伊人婷婷丁香| 午夜福利成人在线免费观看| 国内揄拍国产精品人妻在线| 久久精品国产亚洲av天美| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站| 国产精品国产三级国产av玫瑰| 亚洲在线自拍视频| 国产成人福利小说| 日韩欧美在线乱码| 欧美一区二区亚洲| 看免费成人av毛片| 亚洲五月天丁香| 九九在线视频观看精品| 99久国产av精品国产电影| 1000部很黄的大片| 午夜激情福利司机影院| 91aial.com中文字幕在线观看| 极品教师在线视频| 亚洲成人中文字幕在线播放| 中文字幕精品亚洲无线码一区| 两性午夜刺激爽爽歪歪视频在线观看| 看十八女毛片水多多多| 国产成年人精品一区二区| 热99re8久久精品国产| 国内精品一区二区在线观看| 色综合站精品国产| 韩国高清视频一区二区三区| 日本黄大片高清| 狠狠狠狠99中文字幕| 观看免费一级毛片| 91av网一区二区| 乱系列少妇在线播放| 色网站视频免费| 久久精品久久精品一区二区三区| 国产激情偷乱视频一区二区| 99久久精品热视频| 国产精品伦人一区二区| 亚洲av中文av极速乱| 国内精品宾馆在线| 国产黄片美女视频| 亚洲av一区综合| 亚洲人成网站在线观看播放| 色吧在线观看| 中文字幕亚洲精品专区| 日本与韩国留学比较| 国产黄色视频一区二区在线观看 | 噜噜噜噜噜久久久久久91| 午夜精品一区二区三区免费看| 亚洲综合色惰| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 中文字幕精品亚洲无线码一区| 国产免费福利视频在线观看| 国模一区二区三区四区视频| 国产免费一级a男人的天堂| 51国产日韩欧美| 国产精品福利在线免费观看| 午夜福利成人在线免费观看| 日韩大片免费观看网站 | 国产爱豆传媒在线观看| 国产麻豆成人av免费视频| 人人妻人人澡欧美一区二区| 欧美高清性xxxxhd video| 最近最新中文字幕免费大全7| 免费人成在线观看视频色| 久久久久久大精品| 天天一区二区日本电影三级| 亚洲欧美清纯卡通| 久久欧美精品欧美久久欧美| 成人亚洲精品av一区二区| 我的女老师完整版在线观看| 成人午夜精彩视频在线观看| 日本免费a在线| 国产精品1区2区在线观看.| 在现免费观看毛片| 日产精品乱码卡一卡2卡三| 国产一级毛片在线| 亚洲一级一片aⅴ在线观看| av专区在线播放| 午夜亚洲福利在线播放| 国产精品国产高清国产av| 久久久久九九精品影院| 天天躁夜夜躁狠狠久久av| 成人二区视频| 丰满人妻一区二区三区视频av| 日韩一区二区视频免费看| 在线观看美女被高潮喷水网站| 黄色日韩在线| 水蜜桃什么品种好| 国产一区二区在线观看日韩| 亚洲在线自拍视频| 青青草视频在线视频观看| 夫妻性生交免费视频一级片| 毛片女人毛片| 亚洲欧美精品专区久久| 亚洲天堂国产精品一区在线| 九九热线精品视视频播放| 欧美激情久久久久久爽电影| 国内精品宾馆在线| 村上凉子中文字幕在线| 乱码一卡2卡4卡精品| 丰满乱子伦码专区| 国产精品国产高清国产av| 国产精品久久视频播放| 国内少妇人妻偷人精品xxx网站| 免费看日本二区| 国产精品久久久久久精品电影| av在线老鸭窝| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 午夜激情欧美在线| 亚洲怡红院男人天堂| 国产精品野战在线观看| 色播亚洲综合网| 午夜福利网站1000一区二区三区| 成人二区视频| 亚洲精品乱久久久久久| 国产成人一区二区在线| 亚洲欧美清纯卡通| 久久99蜜桃精品久久| 国产在视频线在精品| 亚洲精品国产成人久久av| 亚洲色图av天堂| 国产精品久久久久久精品电影小说 | 亚洲四区av| 国产三级在线视频| 久久久a久久爽久久v久久| 最近视频中文字幕2019在线8| 亚洲经典国产精华液单| 午夜爱爱视频在线播放| 午夜a级毛片| 成人欧美大片| 久久韩国三级中文字幕| 在线观看av片永久免费下载| 欧美成人精品欧美一级黄| 日韩三级伦理在线观看| 精品少妇黑人巨大在线播放 | 亚洲乱码一区二区免费版| 国产精品蜜桃在线观看| 色噜噜av男人的天堂激情| 免费黄网站久久成人精品| 久久久成人免费电影| 亚洲av成人av| 日韩欧美国产在线观看| 97人妻精品一区二区三区麻豆| 欧美另类亚洲清纯唯美| 精品国产三级普通话版| 亚洲不卡免费看| 日本wwww免费看| 精品99又大又爽又粗少妇毛片| 草草在线视频免费看| 老女人水多毛片| 国产午夜精品久久久久久一区二区三区| 在现免费观看毛片| 边亲边吃奶的免费视频| 国产精品熟女久久久久浪| 免费观看在线日韩| 免费观看的影片在线观看| 两个人的视频大全免费| 乱系列少妇在线播放| 国产极品天堂在线| 亚洲一区高清亚洲精品| 婷婷色综合大香蕉| 色尼玛亚洲综合影院| 在线免费观看的www视频| 亚洲三级黄色毛片| 噜噜噜噜噜久久久久久91| 国产免费男女视频| 色播亚洲综合网| 中文亚洲av片在线观看爽| av女优亚洲男人天堂| 天天躁夜夜躁狠狠久久av| 日本黄色片子视频| 少妇丰满av| 国产亚洲av片在线观看秒播厂 | 寂寞人妻少妇视频99o| 最近的中文字幕免费完整| 亚洲乱码一区二区免费版| 亚洲国产成人一精品久久久| 99久久精品热视频| 国产精品电影一区二区三区| av卡一久久| 国产亚洲5aaaaa淫片| 亚洲在线自拍视频| 日日摸夜夜添夜夜爱| 日韩av不卡免费在线播放| 国产在视频线在精品| 视频中文字幕在线观看| 纵有疾风起免费观看全集完整版 | 秋霞伦理黄片| 亚洲五月天丁香| 国产三级在线视频| 亚洲三级黄色毛片| 亚洲五月天丁香| 亚洲欧美精品专区久久| 国产探花在线观看一区二区| 六月丁香七月| 亚洲乱码一区二区免费版| 中文字幕久久专区| 亚洲婷婷狠狠爱综合网| 亚洲欧美精品专区久久| 青春草视频在线免费观看| 国产精品久久久久久av不卡| 久久99精品国语久久久| 日韩视频在线欧美| av卡一久久| 村上凉子中文字幕在线| 久久精品熟女亚洲av麻豆精品 | 午夜爱爱视频在线播放| 国产精品久久久久久精品电影| 久久久久久久久久成人| 国产人妻一区二区三区在| www.色视频.com| 欧美日韩国产亚洲二区| 97热精品久久久久久| 午夜福利高清视频| 色视频www国产| 久久久久国产网址| 免费av毛片视频| 一级爰片在线观看| 欧美一区二区国产精品久久精品| 伦精品一区二区三区| av视频在线观看入口| 欧美区成人在线视频| 久久精品夜色国产| 最新中文字幕久久久久| 国产一区有黄有色的免费视频 |