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

    快速LZD匹配異源DEM空洞修復(fù)應(yīng)用研究

    2023-12-11 07:11:58李小虎
    計算機工程與應(yīng)用 2023年23期
    關(guān)鍵詞:格網(wǎng)線程空洞

    李小虎,蘆 穎,田 壯

    中國電子科技集團公司第十五研究所,北京 100080

    地形數(shù)據(jù)信息處理技術(shù)是數(shù)字地形分析(digital terrain analyse,DTA)的基礎(chǔ),一直以來作為測繪學(xué)、地貌學(xué)和計算機圖形學(xué)的研究熱點,被廣泛應(yīng)用于道路規(guī)劃、精準農(nóng)業(yè)、數(shù)字戰(zhàn)場、實景三維等眾多領(lǐng)域[1-2]。隨著全球信息化、數(shù)字化的快速發(fā)展,作為一種重要的地理信息資源,數(shù)字高程模型(digital elevation model,DEM)產(chǎn)品的數(shù)據(jù)采集方式更加多元,可根據(jù)用戶的個性化需求定制不同精度、形態(tài)復(fù)雜度、參考框架和覆蓋范圍的數(shù)據(jù)產(chǎn)品[3-4]。其中,機載激光雷達(airborne light detection and ranging,ALiDAR)憑借高密度、高精度的三維數(shù)據(jù)采集優(yōu)勢已成為獲取高質(zhì)量DEM的主要數(shù)據(jù)來源[5-7]。然而,受云霧遮擋、高程地形(陡峭的斜坡和深谷)和信號干擾等因素對立體成像的影響[8],導(dǎo)致構(gòu)建的ALiDAR航測點云在某些區(qū)域散布較多的數(shù)據(jù)空洞-像素高程值為零,從而嚴重影響了地表數(shù)字化模擬產(chǎn)品的可用性[9-10]。

    如何有效填補上述DEM 數(shù)據(jù)空缺,是學(xué)者們進行技術(shù)性研究任務(wù)的關(guān)鍵難點。目前主流方法有兩種[11-13]:一是利用空缺區(qū)域鄰域信息的直接插值法;二是結(jié)合三維數(shù)據(jù)融合與空間插值的綜合法。改進的插值方法有:基于高斯權(quán)值-向量基(Gauss weight-vector base,GW-VB)的插值方法,在明確稀疏與稠密點云的劃定界限后,能夠有效提升離散點間的平滑處理效果[14];等高線間隔重差值法(contour interval redifference,CIR),該方法能夠較好地捕捉到總體的地形趨勢,但是由于缺乏詳細的描述特征,常應(yīng)用于水文區(qū)域的空缺填充[15];高精度曲面模型(high accuracy surface modelling,HASM)插值法:只需借助空缺周圍數(shù)據(jù)進行填補,能夠快速填補位于簡單地形表面的小尺度的數(shù)據(jù)空洞[16]。改進后的自適應(yīng)Kriging 插值法:能夠克服離散點分布不均勻條件下所引入的建模誤差,通過最優(yōu)無偏估計提高插值精度[17]。

    插值法在空洞范圍較小且平坦的數(shù)據(jù)修復(fù)中表現(xiàn)較為理想,對于山脊、山谷等大面積區(qū)域缺失的情況下,常采用三維數(shù)據(jù)融合與空間插值綜合法,以滿足高程產(chǎn)品的質(zhì)量要求。其中,岳國棟等[18]提出了干系數(shù)與垂直基線高度加權(quán)融合的方法,利用一組主輔影像去除側(cè)視異常區(qū)域;李世明[19]提出了多源點云漸進配準和拉普拉斯融合算法,在面向城市三維空間精細建模中完善了局部復(fù)雜的結(jié)構(gòu)表面;基于泊松融合[20]的無縫修復(fù)法,通過統(tǒng)計分析極大地增加了填充可用候選集,且具備復(fù)制紋理的效果;陳雪榮[21]探索出了優(yōu)化后的徑向基神經(jīng)網(wǎng)絡(luò)(radial basis function,RBF)方法,可有效去除大尺度噪聲,并用網(wǎng)絡(luò)預(yù)測值代替實測值有效解決非線性空洞修補;余婷婷等[22]提出的將隨機森林回歸算法(random forest regression,RFR)應(yīng)用在DEM 空洞修補精度校正上,適用于處理各種地形異質(zhì)程度的高程數(shù)據(jù),泛化能力強。上述方法雖在去除大尺度數(shù)據(jù)空洞方面做了大量工作,但當數(shù)據(jù)分辨率較高且空洞覆蓋面積較大時,海量的訓(xùn)練樣本導(dǎo)致模型不易收斂,從而影響了修復(fù)效率。

    本文提出了一種基于異源輔助數(shù)據(jù)的快速最小高差三維匹配算法,配合統(tǒng)一計算架構(gòu)(compute unified device architecture,CUDA)實現(xiàn)GPU 多核策略并行計算,提高LZD 空間匹配模型求解和“臺階”融合算法對空洞邊緣離散點平滑處理的效率和效果。修復(fù)后得到的DEM數(shù)據(jù)完整度高、過程可靠,利用該方法形成的軟件成果已在工程應(yīng)用中取得了較好的實效,具有一定的推廣價值。

    1 方法原理

    1.1 LZD表面匹配原理

    為了進行快速、高精度的空洞修復(fù),首先需要統(tǒng)一輔助數(shù)據(jù)與待修補數(shù)據(jù)間的坐標基準,將兩個結(jié)構(gòu)、位姿差異的三維空間進行表面匹配。目前廣泛使用的方式有最近點迭代算法(iterative closest point,ICP)和LZD 算法等[23]。鑒于待匹配兩數(shù)據(jù)規(guī)則格網(wǎng)的地形表述特點,選用LZD 算法進行待轉(zhuǎn)換空間映射模型的計算,如圖1所示。

    圖1 空間坐標關(guān)系示意圖Fig.1 Spatial coordinate relationship diagram

    依據(jù)兩組DEM 數(shù)據(jù)組織形式,設(shè)待匹配模型為Zc=f(Xc,Yc),其中Xc、Yc是待修復(fù)表面投影至二維平面的坐標值;同理,標準模型為Zs=f(Xs,Ys)。待修復(fù)表面上一點Pc=[Xc,Yc,Zc]T與基準數(shù)據(jù)上對應(yīng)點Ps=[Xs,Ys,Zs]T間的幾何關(guān)系可進行如下表示:

    其中,T=[tx ty tz]T為Pc相對Ps的空間三個維度平移矩陣,s為比例縮放系數(shù),R是含有參數(shù)θ、?、ψ的表示旋轉(zhuǎn)偏移的單位正交矩陣。由于DEM格網(wǎng)地理坐標作為一種相對坐標,不存在偏扭、彎曲等復(fù)雜情況,根據(jù)布爾莎[24]仿射變換模型得到其簡化表達式如下:

    式中,θ、?、ψ分別為繞O-Zs、O-Ys、O-Xs軸的滾轉(zhuǎn)角、偏航角和俯仰角。將式(2)帶入式(1)中可建立待匹配坐標系O-XcYcZc向基準坐標系O-XsYsZs進行空間轉(zhuǎn)換的矩陣模型:

    其中,樣本個數(shù)i=m×n,m、n分別為待匹配柵格的橫、縱坐標像素量。依據(jù)最小二乘(ordinary least squares,OLS)迭代估計方法,取Z方向高差平方和最小值建立目標方程。則最小高差極值的計算公式為:

    存在平面同名點時,格網(wǎng)值的權(quán)重因子wsi=1,否則wsi=0。當基準與待匹配數(shù)據(jù)分辨率不一致時,利用雙線性插值法(bilinear interpolation,BI)進行重采樣,構(gòu)造j個離散的虛擬像素pfj(xj,yj,zj),其中zj=0。

    將局部虛擬空間近似為連續(xù)二次線性曲面:

    通過4 個臨近參考點確定橫縱兩個維度的曲線參數(shù)a0、a1、b0、b1,將pfj帶入式(5)計算等距插值zj=hfj(xj,yj),將數(shù)據(jù)修改為相同尺寸。

    利用最小二乘迭代求解(4)中E接近限差時的7個約束方程參數(shù),從而構(gòu)建三維匹配的空間坐標轉(zhuǎn)換模型。

    1.2 臺階融合

    受不同時相、傳感幾何輻射差異等因素的影響,空洞銜接邊緣會產(chǎn)生不可避免的高程值突變[25],針對存在的“臺階”式拼接縫,利用圖像融合技術(shù)對其進行平滑處理。

    經(jīng)臺階融合算法處理得到的高程坐標Zfusion,是根據(jù)其自身每個格網(wǎng)處的高程值Zdem與同名點高程值Zref按如下規(guī)則加權(quán)平均得到:

    式中,pdem表示待處理DEM 數(shù)據(jù)點的權(quán)重,與其對應(yīng)的同名參考點權(quán)值之和為1。結(jié)合離散圖像相鄰區(qū)域像素值變化小于區(qū)域內(nèi)部的特點,可利用柵格距離衰減規(guī)律對每個格網(wǎng)處的pdem進行計算。

    若取窗口大小為n×n,中心點權(quán)重p(xc,yc)=0.5,其余8個方向n2-1 個點的權(quán)重和之為0.5,則待修復(fù)數(shù)據(jù)鄰域內(nèi)各點的高程梯度倒數(shù)絕對值為:

    式中,像素偏移量i、j∈[0,n-1],且i、j不同時為0,若Zdem(xc+i,yc+j)=Zdem(xc,yc),則gdem(xc,yc,i,j)=0。

    以待修復(fù)點作為掩模中心進行加權(quán)鄰域求梯度平均,可得到待處理高程權(quán)值系數(shù):

    將pdem帶入式(6),經(jīng)臺階融合輸出空洞填充邊緣的高程坐標信息,使格網(wǎng)像素能夠自然、平滑地銜接在一起。

    實現(xiàn)過程如下:

    (1)設(shè)定每個格網(wǎng)點初始化的權(quán)值p0如下:

    i、j分別為DEM數(shù)據(jù)的行號和列號,有效初值為215,為確定待修復(fù)點的高程權(quán)值分布,需要對圖像順序、逆序兩次進行計算。

    (2)定義一個偏移量offset:?。ǜ鞣较颍?quán)值重新計算的格子數(shù)默認為20。

    (3)以柵格平面的左上角點作為原點,行號、列號依次遞增,遍歷尋找空洞邊緣點,取該像元鄰域內(nèi)左、左上、上、右上方向上各offset∈[1,(n-1)/2]個點并根據(jù)梯度計算依次設(shè)置權(quán)重。

    (4)更改原點為像平面右下角,行號、列號依次遞減,在左下、下、右下、右方向上逆序重復(fù)上述操作,更新上一步重新計算的像元的權(quán)重pdem=min{p(xc,yc),pdem} 。

    (5)取滿足pdem(x,y)<1 的所有點,利用加權(quán)平均的方法重新計算高程值,以消除待修復(fù)數(shù)據(jù)存在的臺階問題。

    2 快速DEM空洞修復(fù)方案

    考慮到大尺度數(shù)據(jù)空洞修復(fù)時存在的實時性問題,首先,基于LZD 原則將表面匹配中三個維度的殘差分析簡化為取Z 方向極小值,進而利用GPU 并行處理機制,將同名點搜索和格網(wǎng)點像元值的填補過程分解至多個線程塊同步進行,以縮短處理過程的總體運算周期。

    2.1 匹配模型參數(shù)獲取

    為了得到最優(yōu)的表面匹配結(jié)果,通過最小二乘迭代逐步優(yōu)化模型數(shù)值,確定各參數(shù)近限差的增量矩陣dX。對變換陣Trans的7 參選取初始值X0=[R0,t0,s0]T,其中令Rψ0、R?0、Rθ0、tx0、ty0、tz0=0,s0=1 表示空間兩坐標系重合,無相對旋轉(zhuǎn)、平移,默認等比例縮放。

    遍歷待匹配樣本中的每個高程點Pci=[Xci,Yci,Zci]T,帶入式(3)進行Trans變換后,在標準DEM中尋找格網(wǎng)坐標滿足xci=xsi、yci=ysi的同名點Psi。

    當同名點不存在時,本算法中使用了雙線性插值構(gòu)造虛擬同名點的高程坐標,數(shù)據(jù)量小時也可考慮反距離加權(quán)平方(inverse distance weight,IDW)算法。

    給定一個足夠小的步長delta,當前點在X、Y方向的一階導(dǎo)數(shù)計算公式如下:

    式中,待匹配所有點滿足Zc=f(Xc,Yc)。對誤差函數(shù)e(x,y)進行一階泰勒展開,得到線性化誤差函數(shù):

    每次迭代各點誤差函數(shù)的過程就是求最優(yōu)增量參數(shù)矩陣dX*的過程:

    經(jīng)系數(shù)反推,可建立增量方程組B=AdX,設(shè)計觀測系數(shù)矩陣Ai=[dzxdzy-1]T,得到Bi:

    將Ai*Bi的結(jié)果存儲到A矩陣的第i行,Zci-Zsi的值存儲到B矩陣的第i行。

    遍歷結(jié)束后,使用奇異值分解(singular value decomposition,SVD)方法解超定方程以計算dX增量:

    并將dX*累加到變換矩陣中X=X+dX*,如果dX小于迭代限差Vx或達到最大迭代次數(shù)Va時,跳出迭代過程完成參數(shù)計算。

    在匹配過程中,迭代終止條件由各參數(shù)相鄰迭代值之差來判定,即七個轉(zhuǎn)換參數(shù)的增量,七個轉(zhuǎn)換參數(shù)的限差值和設(shè)定的最大迭代次數(shù)由各個觀測值的精度決定。

    匹配模型參數(shù)獲取的具體流程如圖2所示。

    圖2 算法流程圖Fig.2 Algorithm flow chart

    2.2 空洞修復(fù)任務(wù)的GPU并行方案設(shè)計

    引入最小高差算法雖然避免了三個維度最小距離計算造成的冗余問題,但在進行同名點搜索、鄰域加權(quán)時,窮舉式的距離查找過程和大規(guī)模的DEM 數(shù)據(jù)量依然影響了空間匹配的算法效率。

    因此,加入GPU并行處理技術(shù),通過PCIe總線協(xié)助CPU 進行運算控制操作,充分發(fā)揮CPU 的復(fù)雜邏輯控制能力和GPU在處理圖像高密集計算和超長流水線任務(wù)的結(jié)構(gòu)優(yōu)勢,形成的CPU-GPU異構(gòu)模型如圖3所示。

    圖3 CPU-GPU異構(gòu)模型Fig.3 CPU-GPU heterogeneous model

    NVIDA 推出的標準運算平臺CUDA,代碼層面更易編寫,在性能、成本維護方面表現(xiàn)突出。本文利用其并行計算特性,提出一種多核多線程并發(fā)的DEM 匹配修復(fù)方案。其中,并行化編程最需考慮的線程同步和數(shù)據(jù)歸并問題,利用CUDA平臺線程調(diào)度和內(nèi)存的訪問管理原則,將空洞修復(fù)計算過程以核-格網(wǎng)-線程快-線程(KERNAL-GRID-BLOCK-THREAD,KGBT)的三級組織架構(gòu)進行調(diào)度配置。

    構(gòu)造順序的處理邏輯步驟如下:

    (1)獲取所有待匹配格網(wǎng)像元同名點坐標值;

    (2)迭代計算殘差極值并求解空間映射模型的特征向量;

    (3)利用加權(quán)融合算法重新計算空洞邊緣高程坐標。

    獲取同名點坐標和加權(quán)融合過程需要進行大量重復(fù)計算,因此分別映射至knl_search、knl_restore 兩個核函數(shù)中,并采用單指令多數(shù)據(jù)流(single instruction multiple data,SIMD)以異步的形式執(zhí)行相同的算法指令:

    首先,以程序初始化和影像數(shù)據(jù)GeoTiff(tiff),Erdas Imagine Images(img),ASCII DEM(dem)等格式的解析作為起始模塊。利用GDAL(geospatial data abstraction library)庫讀取元數(shù)據(jù)和高程信息數(shù)組取至主機端內(nèi)存。

    針對基準模型大小遠大于待匹配區(qū)域的現(xiàn)象,利用OpenCv庫對存儲DEM的數(shù)據(jù)模型文件進行裁剪處理,調(diào)用RasterIO 將大圖分割為512×512 大小的塊后,將柵格數(shù)組與元數(shù)據(jù)一同拷貝至全局顯存(gloabl memory,GM),完成程序的初始化。

    設(shè)置臨近點計算的輻射距離為0.2 m,采用CUDA編程模型并調(diào)用核函數(shù)knl_search,以線程id 作為檢索項,將遍歷像點的任務(wù)映射成若干并發(fā)執(zhí)行的子任務(wù)。在某一線程網(wǎng)格(Grid)的多個線程塊(Block)中對區(qū)域分割后的子圖進行同名點搜索,利用異步流迭代方法將同名點對應(yīng)的坐標集合存儲至當前Block(32×32)的共享存儲器中,供塊內(nèi)線程共用,減少訪問延遲以提高讀寫效率。在knl_search計算模型中,對于每個線程,分別搜索當前待匹配格網(wǎng)距離相近的所有點,最近距離小于等于閾值時確定為同名點對。同名點不存在時,設(shè)置橫縱縮放比例和幾何中心利用BI 插值函數(shù)構(gòu)造虛擬點,resize調(diào)整數(shù)據(jù)至相同尺寸。待所有樣本完成查詢操作后,輸出參數(shù)控制所有點對信息進行坐標反算。

    迭代求解映射參數(shù)向量的過程,是采用奇異值分解(singular value decomposition,SVD)對稠密矩陣進行降維,其計算代價大,但由于運算先后關(guān)聯(lián)性較強,內(nèi)部邏輯耦合度高,不適合進行并行化改造,因此該矩陣運算仍在CPU串行域中進行。

    待模型數(shù)據(jù)解析完備,針對圖中多個待修復(fù)區(qū)域進行的8鄰域權(quán)值計算,并將多個空洞分配給多個線程塊block同步進行“臺階”融合。正、逆遍歷所有待融合點,計算中心點p(xc,yc,zc)與鄰域各點高程權(quán)重關(guān)系。對拼接邊緣像素進行平滑處理,以函數(shù)返回值形式合并修復(fù)后數(shù)據(jù)至CPU內(nèi)存,將結(jié)果模型寫入輸出的DEM文件中,完成空洞修復(fù)運算。

    任務(wù)基本設(shè)置:

    (1)分配2個順序進行的核函數(shù)(kernel),進行同名點插值計算和臺階像素融合。將knl_search函數(shù)設(shè)置為逐點搜索,n個點對應(yīng)n個線程。等待圖像所有點遍歷結(jié)束,通知主線程迭代求解單應(yīng)性矩陣;對于knl_restore函數(shù),若相鄰Z坐標缺失個數(shù)>1,即判定為待修復(fù)空洞,采用流水線作業(yè)實行多個空洞邊緣并發(fā)修復(fù)。

    (2)knl_search 的m個子圖分別對應(yīng)m個格網(wǎng),將格網(wǎng)劃分為多個pixel 的線程塊(block),每個線程塊分別對應(yīng)256個CUDA線程。

    (3)子圖中第k個空洞映射至knl_restore 的第k個線程塊中。

    采用程序并行化關(guān)鍵解決的幾點問題:

    (1)在大樣本點集查詢?nèi)蝿?wù)中,采用K近鄰法(Knearest neighbor,KNN)查詢單元格間的拓撲關(guān)系,通過兩步法確定法向最小距離點。KNN算法的搜索過程可有以下三種常用方式:窮舉法、kd樹(k-dimensional-tree)和球樹(ball-tree),根據(jù)預(yù)測平均時間隨樣本數(shù)量變化曲線知[26],窮舉的方式更適用于容量較大的DEM數(shù)據(jù)。

    (2)設(shè)Q為大小為u的n維參考點集,R為大小為v的查詢點集,通過線性掃描對k個最近鄰域點進行表決。表面匹配算法只關(guān)心距離最小點的坐標值,在此基礎(chǔ)上,對v個長度為k點集的分別進行插入排序。這種暴力實現(xiàn)的方式非常耗費運算資源,而結(jié)果集合的搜索與有序表的生成在匹配點之間不存在任何依賴關(guān)系,因此,可充分利用GPU 編程將遍歷格網(wǎng)數(shù)據(jù)進行同名點查尋詢的任務(wù)映射成若干并發(fā)執(zhí)行的子任務(wù),提高算法效率。

    (3)融合算法是根據(jù)所有點鄰域加權(quán)取均值,進行Z坐標平滑。這一過程在程序中,遍歷待修復(fù)DEM 所有點獲取權(quán)重是獨立且高重復(fù)性的,因此將任務(wù)分配至GPU運算單元重疊執(zhí)行,在并行層面分別對所有待處理點高程進行重算,將循環(huán)遍歷的過程原子化,顯著降低了內(nèi)存占用和運算時間。

    CPU并行架構(gòu)下,設(shè)3維參考點集Q的樣本數(shù)為待匹配柵格大小i×j,查詢點集為v,由于核心較少導(dǎo)致能夠同時開啟的線程數(shù)目有限,可得鄰域搜索的時間復(fù)雜度為O(ijuv)。結(jié)合圖形處理單元加速運行該算法,每個樣本單獨開啟一個核函數(shù)同步處理,在相同數(shù)據(jù)量的情況下,時間復(fù)雜度減少至O(uv)。然后,對于搜索到的臨近點集,進行v個長度為k的插入排序,該操作的時間復(fù)雜度為O(vk2),使用GPU 并開啟v個線程并行執(zhí)行可將時間復(fù)雜度降低為O(k2)。在邊緣平滑算法的優(yōu)化中,逐點高程值重算被簡化為同步進行的窗口(n×n)卷積,計算時間縮短樣本個數(shù)倍,時間復(fù)雜度為O(n2)。從算法理論角度分析該策略具備明顯的優(yōu)化性能,而系統(tǒng)實際運行的總時長,不僅需要對算法執(zhí)行指令的條數(shù)進行評估,還涉及到了線程工作量后的步驟復(fù)雜度。

    圖4 為GPU 多核架構(gòu)與多線程作業(yè)下的具體實現(xiàn)流程。

    圖4 GPU多核多線程空洞修復(fù)程序流程Fig.4 GPU multi-core multi-thread hole repair program

    3 實驗內(nèi)容

    3.1 實驗數(shù)據(jù)及平臺

    本文選取的基準數(shù)據(jù)來源于奮進號航天飛機的雷達地形測繪(shuttle radar topography mission,SRTM)數(shù)據(jù),全圖數(shù)據(jù)包含1 200×1 200個采樣點高程值,水平和垂直精度分別為±20 m和±16 m,格網(wǎng)間距90 m(3″),投影基準采用WGS84 坐標系統(tǒng),經(jīng)緯范圍(102.4 E~104 E,29 N~30.3 N)。待修復(fù)數(shù)據(jù)均為中國境域采集的SRTM-DEM 30數(shù)據(jù),格網(wǎng)間距30 m(1″)。

    分別對待修復(fù)數(shù)據(jù)中海拔較高的青藏高原(區(qū)域1)和地勢相對平緩的長江三角洲沖積平原(區(qū)域2)的山脊數(shù)據(jù)缺失和條狀云層遮擋進行算法驗證,其中區(qū)域1空洞占比4.5%,連通空洞個數(shù)為8;區(qū)域二占比2.2%,空洞個數(shù)為4。

    實驗環(huán)境:軟件平臺采用Visual Studio2010和C/C++語言進行軟件設(shè)計與編譯,通過觀測模型的收斂速度、修復(fù)總耗時、加速比等,在搭建的三個異構(gòu)混合平臺上依次進行對照組實驗,以驗證本文方法可用性和實現(xiàn)效率。工作站所用處理器(CPU)及顯卡(GPU)型號說明和預(yù)先設(shè)置匹配模型收斂條件分別如表1、表2所示。

    表1 CPU/GPU參數(shù)Table 1 CPU/GPU parameters

    表2 空間映射模型7參增量限差Table 2 Seven parameters for spatial mapping model

    根據(jù)兩個DEM(或DSM 等)的元數(shù)據(jù)信息,設(shè)定Geotiff格式數(shù)據(jù)點平面的坐標。以標準DEM的左上角點作為原點Os(0,0),Xs(行)方向每行間隔為經(jīng)度方向分辨率×111 000 m,Ys(列)方向間隔為緯度方向的分辨率×111 000 m。待匹配DEM 數(shù)據(jù)設(shè)定為Os(0,0)-XcYc原理同標準數(shù)據(jù)。

    在修復(fù)前進行數(shù)據(jù)準備:

    (1)將待匹配DEM數(shù)據(jù)用GDAL庫讀入內(nèi)存,包括元數(shù)據(jù)信息和高程數(shù)組;

    (2)將標準DEM中對應(yīng)待匹配DEM區(qū)域(略大)提取為算法中的標準DEM;

    (3)根據(jù)精度需求進行異源數(shù)據(jù)等間隔采樣;

    (4)切分地形區(qū)域為若干(m)子圖;

    3.2 實驗結(jié)果與分析

    LZD 算法的模型解算精度取決于各參數(shù)限差和迭代次數(shù),鑒于兩組地區(qū)地理面積相近,此處均采用最大迭代次數(shù)300,采樣大小500×500 pixel。當所有參數(shù)增量滿足迭代限差時,算法終止,并確定模型參數(shù)向量。

    分別取待匹配DEM 中青藏高原山區(qū)、長江三角洲沖積平原兩組典型的數(shù)據(jù)模型作為對照組進行格網(wǎng)數(shù)據(jù)空洞修復(fù),以驗證該方法在不考慮實現(xiàn)策略時,不同地勢條件下算法的可行性。其中,青藏高原地區(qū)地勢變化較大,受傳感范圍限制,在數(shù)據(jù)的采集和重構(gòu)過程中易丟失高海拔山脊處的三維信息;長江三角洲平原位于長江入???,該地區(qū)平均坡度較小,但由于厚云遮擋導(dǎo)致內(nèi)含兩條狹長的數(shù)據(jù)空缺。兩組待修復(fù)數(shù)據(jù)大小均為26 MB,全國基準DEM 大小為5 GB,經(jīng)模型求解和高程坐標估值計算后,得到修復(fù)后的結(jié)果如圖5所示。

    圖5 漏洞修補對比圖Fig.5 Comparative diagram of vulnerability patching

    由圖5可知,對于兩地區(qū)分布的不規(guī)則空洞的修復(fù)效果表現(xiàn)為:均能夠較好地恢復(fù)局部紋理特征、幾何性質(zhì)穩(wěn)定,在源數(shù)據(jù)的基礎(chǔ)上保持現(xiàn)有分辨率,并擴大了數(shù)據(jù)的空間覆蓋,實現(xiàn)了影像質(zhì)量的有效提升。

    在盡可能保證精度的前提下,取區(qū)域1的數(shù)據(jù)進行計算性能測試。分別采用CPU多線程和CUDA多線程調(diào)度策略,設(shè)置大小相同的數(shù)據(jù)區(qū)域(區(qū)域1)和采集密度(500×500 pixel),依次記錄模型的匹配耗時與空洞修復(fù)總耗時。配備的三組硬件實驗環(huán)境為,平臺1:CPU為Xeon CPU E3-1230V6,顯卡采用NVIDA GeFouce GTX 980;平臺2:CPU 為Xeon CPU E3-1230V6,顯卡采用NVIDA GeFouce GTX 1080 Ti;平臺3:CPU 為Intel Core i7-7700,顯卡采用NVIDA GeFouce GTX 1080 Ti。在多種主機設(shè)備組合配置下,對照并行算法的執(zhí)行效果。

    實驗結(jié)果表明,部分任務(wù)調(diào)用CPU-GPU 異構(gòu)的CUDA 編程模型進行實現(xiàn)后,修復(fù)總時長<0.86 s,具體運行效率如表3所示。

    表3 不同策略下區(qū)域1修復(fù)耗時與加速比Table 3 Time-consuming and speed-up ratio of Region 1 repair under different strategies

    采用GPU 技術(shù)協(xié)助CPU 處理重復(fù)計算,縮短了點對最優(yōu)距離計算和排序時間,然而隨像素尺寸增加數(shù)據(jù)在CPU-GPU 間拷貝的開銷增大,擴容后的多線程CPU策略(平臺2)較平臺1 總耗時明顯縮減,且匹配耗時占比由20.4%增加至26.1%,串行優(yōu)勢在充足內(nèi)存容量條件下更為顯著,同時,并行計算效率也提升了3.6倍。平臺1與平臺2采用相同CPU和相異的圖形處理器,利用并行異構(gòu)策略,使總執(zhí)行加速比分別提升至4.0 和4.5,優(yōu)化后的匹配、融合耗時較CPU 多線程最多可減少71.082%和76.006%。在平臺3上運行異構(gòu)程序,算法總體修復(fù)的執(zhí)行效率相較于單線程CPU 能夠達到9.7 倍以上的加速比,具有較好的修復(fù)性能。

    為了檢驗本文算法在不同地形場景下的匹配精度,分別將待修復(fù)數(shù)據(jù)中的區(qū)域1和區(qū)域2與全國基準數(shù)據(jù)進行LZD 三維匹配。其中,青藏山區(qū)和沖積平原數(shù)據(jù)大小為各取50 個點作為測試樣本,與基準數(shù)據(jù)對照計算Z坐標估值殘差,得到的最小參數(shù)增量、誤差矩陣和Z方向精度均方誤差(單位:m)如表4所示。

    表4 迭代最終參數(shù)增量與誤差Table 4 Iterate over final parameter increment and error

    以區(qū)域1的山脊為例,可看到模型計算的高程值在空洞邊緣存在“臺階”形斷層,為完善修復(fù)效果,對該空洞區(qū)域進行了加權(quán)融合。取空洞邊緣點Z坐標值得加權(quán)平均,進行DEM高程值重算,再次賦值后達到奇異點去除和邊緣平滑的效果,空洞平滑修補前后對比如圖6所示。

    圖6 空洞平滑修補前后對比圖Fig.6 Comparative diagram of cavity before and after smooth repair

    如圖6(a)所示,融合前的高海拔數(shù)據(jù)缺失面,邊沿分布著非連續(xù)的像素條帶,在數(shù)據(jù)準確度和視覺效果層面都存在明顯缺陷。經(jīng)加權(quán)融合糾正后的該區(qū)域見圖6(b),保留了主要地形特征,未見明顯修復(fù)痕跡,達到趨近于基準數(shù)據(jù)的完善要求。

    4 結(jié)束語

    機載雷達(airborne radar,AR)在高海拔、崎嶇地勢和極端氣候時出現(xiàn)的測繪盲區(qū),導(dǎo)致衍生出的數(shù)字高程數(shù)據(jù)(DEM)存在數(shù)據(jù)空洞等的問題,現(xiàn)有的空間插值和表面匹配算法為該空缺修復(fù)研究提供了有力支撐。為了進一步滿足準確、高效的應(yīng)用需求,本文提出了一種改進的快速LZD三維匹配算法結(jié)合“臺階”融合的修復(fù)方法。借助異源輔助數(shù)據(jù),在搭建的不同平臺下,分別對2 組差異明顯的非結(jié)構(gòu)化地形樣本進行了方法驗證。實驗結(jié)果表明本文采用的方法,修復(fù)后的圖像信息豐富、空洞邊緣平滑,誤差控制在±37.725 m,符合地形數(shù)據(jù)的處理需求。在模型精度滿足固定限差的條件下匹配效率成倍提升,較CPU 多線程架構(gòu)運行耗時可減少76%,加速比最優(yōu)可達9.7倍,在大尺度數(shù)字高程模型修復(fù)領(lǐng)域具有明顯優(yōu)勢。

    猜你喜歡
    格網(wǎng)線程空洞
    實時電離層格網(wǎng)數(shù)據(jù)精度評估
    淺談linux多線程協(xié)作
    空洞的眼神
    基于空間信息格網(wǎng)與BP神經(jīng)網(wǎng)絡(luò)的災(zāi)損快速評估系統(tǒng)
    用事實說話勝過空洞的說教——以教育類報道為例
    新聞傳播(2015年20期)2015-07-18 11:06:46
    平均Helmert空間重力異常格網(wǎng)構(gòu)制方法
    基于位置服務(wù)的地理格網(wǎng)編碼設(shè)計
    測繪通報(2013年2期)2013-12-11 07:27:50
    臭氧層空洞也是幫兇
    Linux線程實現(xiàn)技術(shù)研究
    么移動中間件線程池并發(fā)機制優(yōu)化改進
    国产精品一区www在线观看| 欧美亚洲日本最大视频资源| 亚洲图色成人| 久久毛片免费看一区二区三区| 亚洲精品美女久久av网站| 亚洲精品亚洲一区二区| 国产亚洲午夜精品一区二区久久| 美女cb高潮喷水在线观看| 美女视频免费永久观看网站| 日韩视频在线欧美| 你懂的网址亚洲精品在线观看| 国产黄频视频在线观看| 在线观看美女被高潮喷水网站| 国产精品久久久久成人av| 亚洲精品中文字幕在线视频| 久久毛片免费看一区二区三区| 两个人的视频大全免费| 夫妻午夜视频| 麻豆乱淫一区二区| 亚洲国产精品999| 天天躁夜夜躁狠狠久久av| 大码成人一级视频| 曰老女人黄片| 成人午夜精彩视频在线观看| 日本黄大片高清| 99久久中文字幕三级久久日本| 免费观看无遮挡的男女| 乱人伦中国视频| 亚洲精品美女久久av网站| 午夜免费男女啪啪视频观看| 午夜免费观看性视频| 亚洲国产精品一区二区三区在线| 日日撸夜夜添| 欧美 亚洲 国产 日韩一| 精品99又大又爽又粗少妇毛片| 亚洲久久久国产精品| 免费高清在线观看视频在线观看| 国产精品欧美亚洲77777| 国产精品熟女久久久久浪| 亚洲精品乱码久久久v下载方式| 久久久久国产网址| 亚洲国产av影院在线观看| 黄色毛片三级朝国网站| 十八禁高潮呻吟视频| 啦啦啦在线观看免费高清www| 日韩欧美精品免费久久| 又黄又爽又刺激的免费视频.| 老司机影院毛片| av免费观看日本| a级毛片在线看网站| 中文字幕精品免费在线观看视频 | 天堂中文最新版在线下载| 亚洲中文av在线| 日韩熟女老妇一区二区性免费视频| 伊人久久精品亚洲午夜| 黄色配什么色好看| 国产av码专区亚洲av| 亚洲人成77777在线视频| 老司机影院毛片| 男女啪啪激烈高潮av片| 欧美精品人与动牲交sv欧美| 亚洲av日韩在线播放| 日韩一本色道免费dvd| av福利片在线| 亚洲精品乱久久久久久| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 成人18禁高潮啪啪吃奶动态图 | 国产成人a∨麻豆精品| 午夜福利影视在线免费观看| 一本一本综合久久| 国产黄色视频一区二区在线观看| 蜜桃国产av成人99| 欧美亚洲日本最大视频资源| 国产深夜福利视频在线观看| videos熟女内射| 亚洲人成网站在线观看播放| 欧美激情极品国产一区二区三区 | 永久网站在线| 草草在线视频免费看| 啦啦啦视频在线资源免费观看| 亚洲av不卡在线观看| 九草在线视频观看| 国产精品国产三级国产av玫瑰| 婷婷色综合www| 国产一区有黄有色的免费视频| 国产精品久久久久久久电影| 国产欧美另类精品又又久久亚洲欧美| 亚洲五月色婷婷综合| 久久久欧美国产精品| 18在线观看网站| 黑人欧美特级aaaaaa片| 亚洲国产精品国产精品| 91久久精品电影网| 九色成人免费人妻av| 老司机影院成人| 久久狼人影院| 亚洲美女黄色视频免费看| 久久久久久久精品精品| 中文字幕制服av| 国产一区二区三区综合在线观看 | 国产免费视频播放在线视频| 亚洲人成网站在线播| .国产精品久久| 色吧在线观看| 老司机影院成人| av视频免费观看在线观看| 亚洲av国产av综合av卡| 免费观看无遮挡的男女| 美女中出高潮动态图| 精品人妻偷拍中文字幕| 丰满迷人的少妇在线观看| 999精品在线视频| 日本黄色日本黄色录像| 草草在线视频免费看| 美女内射精品一级片tv| 秋霞伦理黄片| 国产一区二区三区av在线| 黄色毛片三级朝国网站| 国产亚洲精品第一综合不卡 | 久久久国产一区二区| 少妇 在线观看| 亚洲,欧美,日韩| 久久精品久久久久久噜噜老黄| 成人无遮挡网站| 又黄又爽又刺激的免费视频.| 高清午夜精品一区二区三区| 久久免费观看电影| 人人妻人人爽人人添夜夜欢视频| 狂野欧美白嫩少妇大欣赏| av不卡在线播放| 51国产日韩欧美| 999精品在线视频| 成年av动漫网址| 精品亚洲成国产av| 午夜视频国产福利| 久久久国产精品麻豆| 国产成人精品无人区| 日韩中字成人| 大片电影免费在线观看免费| 国产黄片视频在线免费观看| 女性被躁到高潮视频| 自线自在国产av| 欧美激情 高清一区二区三区| 国产高清不卡午夜福利| 亚洲精品自拍成人| 毛片一级片免费看久久久久| tube8黄色片| 寂寞人妻少妇视频99o| 多毛熟女@视频| 天堂俺去俺来也www色官网| 91成人精品电影| 全区人妻精品视频| 亚洲情色 制服丝袜| 中文字幕免费在线视频6| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久| 欧美 日韩 精品 国产| 日本免费在线观看一区| 黄片播放在线免费| 一级二级三级毛片免费看| 成年av动漫网址| av有码第一页| 大香蕉久久网| 五月天丁香电影| 国产免费视频播放在线视频| 丝袜脚勾引网站| 高清在线视频一区二区三区| 欧美激情国产日韩精品一区| 亚洲一区二区三区欧美精品| 亚洲国产av新网站| 国产成人免费无遮挡视频| 国产 精品1| 中文字幕人妻熟人妻熟丝袜美| 不卡视频在线观看欧美| 91午夜精品亚洲一区二区三区| 久久久久久伊人网av| 免费观看无遮挡的男女| 久久99热6这里只有精品| 在线观看一区二区三区激情| 色94色欧美一区二区| 精品久久久久久电影网| 精品国产国语对白av| 久久久久久久久久久丰满| 看免费成人av毛片| 日韩伦理黄色片| 一级毛片 在线播放| 国产免费福利视频在线观看| 狠狠婷婷综合久久久久久88av| 久久久亚洲精品成人影院| 色哟哟·www| a 毛片基地| 日韩成人伦理影院| 制服诱惑二区| 精品久久久久久久久av| 亚洲四区av| 国产亚洲一区二区精品| 欧美变态另类bdsm刘玥| 精品国产露脸久久av麻豆| 欧美亚洲 丝袜 人妻 在线| 国产成人91sexporn| 亚洲欧美日韩卡通动漫| 国产黄片视频在线免费观看| 国产成人免费无遮挡视频| 国产精品久久久久久av不卡| 久久午夜福利片| 午夜91福利影院| 插阴视频在线观看视频| 高清不卡的av网站| 日日摸夜夜添夜夜爱| 99九九在线精品视频| 亚洲av免费高清在线观看| 亚洲精品日本国产第一区| 久久人人爽人人片av| 欧美日韩亚洲高清精品| 国产乱人偷精品视频| 一个人看视频在线观看www免费| freevideosex欧美| 日韩成人av中文字幕在线观看| 少妇 在线观看| 日韩电影二区| 这个男人来自地球电影免费观看 | av一本久久久久| 女性被躁到高潮视频| 国产亚洲精品第一综合不卡 | 纵有疾风起免费观看全集完整版| 成人国语在线视频| 一级毛片我不卡| 一区二区日韩欧美中文字幕 | 国产精品免费大片| 老女人水多毛片| 视频在线观看一区二区三区| videossex国产| 尾随美女入室| 一级a做视频免费观看| 18禁动态无遮挡网站| 欧美精品一区二区大全| 亚洲综合色网址| 久久国内精品自在自线图片| 亚洲av在线观看美女高潮| 全区人妻精品视频| 久久ye,这里只有精品| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 久久 成人 亚洲| av视频免费观看在线观看| www.色视频.com| 亚洲欧美一区二区三区国产| 夫妻性生交免费视频一级片| 久久ye,这里只有精品| 日本91视频免费播放| 免费人妻精品一区二区三区视频| 欧美 亚洲 国产 日韩一| 亚洲性久久影院| 亚洲综合色网址| 亚洲精品国产av蜜桃| 美女大奶头黄色视频| 久久午夜福利片| 色5月婷婷丁香| 五月玫瑰六月丁香| av黄色大香蕉| 少妇丰满av| 国产日韩欧美视频二区| 国产成人精品久久久久久| av国产精品久久久久影院| 十八禁高潮呻吟视频| 久久99精品国语久久久| 欧美精品亚洲一区二区| 免费观看av网站的网址| 好男人视频免费观看在线| 国产毛片在线视频| 欧美少妇被猛烈插入视频| 高清毛片免费看| 91国产中文字幕| 精品一区二区三区视频在线| 国产精品一区二区在线观看99| 美女xxoo啪啪120秒动态图| 欧美精品高潮呻吟av久久| 日韩伦理黄色片| 久久精品国产亚洲网站| 中文天堂在线官网| 能在线免费看毛片的网站| 一个人看视频在线观看www免费| 有码 亚洲区| 久久精品国产鲁丝片午夜精品| 免费人妻精品一区二区三区视频| 黄色视频在线播放观看不卡| 国产精品久久久久久久电影| 热re99久久国产66热| 日本vs欧美在线观看视频| 国产精品99久久99久久久不卡 | 久久久久精品性色| 91精品三级在线观看| 日本色播在线视频| 国产成人精品在线电影| 精品熟女少妇av免费看| 寂寞人妻少妇视频99o| 有码 亚洲区| 亚洲av免费高清在线观看| 国产亚洲av片在线观看秒播厂| 美女大奶头黄色视频| 欧美日韩视频高清一区二区三区二| 亚洲国产欧美日韩在线播放| 精品酒店卫生间| 亚州av有码| 少妇人妻精品综合一区二区| 黄色配什么色好看| 欧美精品国产亚洲| 国产一区二区三区av在线| 亚洲欧美日韩卡通动漫| 韩国高清视频一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 3wmmmm亚洲av在线观看| 成人无遮挡网站| 久久婷婷青草| 亚洲天堂av无毛| 亚洲国产精品成人久久小说| 伦精品一区二区三区| 一区二区av电影网| 国产精品久久久久成人av| 另类亚洲欧美激情| av免费观看日本| 好男人视频免费观看在线| 日韩电影二区| 久久99一区二区三区| 成人国语在线视频| 国产乱来视频区| 亚洲av二区三区四区| 精品一区在线观看国产| 国产日韩欧美亚洲二区| 国产成人精品无人区| 色视频在线一区二区三区| 亚洲精品一二三| 97在线视频观看| 黄色欧美视频在线观看| 精品久久蜜臀av无| 在线观看一区二区三区激情| 亚洲一级一片aⅴ在线观看| 亚洲精品一区蜜桃| 日本vs欧美在线观看视频| 高清毛片免费看| 精品国产乱码久久久久久小说| 亚洲精品456在线播放app| 在线观看免费日韩欧美大片 | 又大又黄又爽视频免费| 国产深夜福利视频在线观看| 草草在线视频免费看| 午夜精品国产一区二区电影| 国产熟女午夜一区二区三区 | 中国三级夫妇交换| 一级爰片在线观看| 亚洲内射少妇av| 男女无遮挡免费网站观看| 国产一区亚洲一区在线观看| 国产高清三级在线| 久久久国产一区二区| 日本欧美视频一区| 国产成人精品婷婷| 亚洲少妇的诱惑av| 狠狠婷婷综合久久久久久88av| av免费在线看不卡| 水蜜桃什么品种好| 久久久久久久久久久免费av| 男男h啪啪无遮挡| 国产欧美亚洲国产| 丝袜美足系列| 美女脱内裤让男人舔精品视频| 国产免费一区二区三区四区乱码| 色94色欧美一区二区| 国产成人a∨麻豆精品| 性高湖久久久久久久久免费观看| 日日摸夜夜添夜夜爱| 少妇高潮的动态图| 亚洲色图综合在线观看| 伦理电影大哥的女人| 男男h啪啪无遮挡| videosex国产| 最近手机中文字幕大全| 国产成人精品无人区| 日韩亚洲欧美综合| 成人毛片a级毛片在线播放| 亚洲欧美中文字幕日韩二区| 欧美精品人与动牲交sv欧美| 九色亚洲精品在线播放| 亚洲av国产av综合av卡| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 热re99久久精品国产66热6| 久久99精品国语久久久| 欧美成人精品欧美一级黄| 久久精品国产自在天天线| 国产爽快片一区二区三区| 精品国产一区二区久久| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲丝袜综合中文字幕| 中国国产av一级| 亚洲,欧美,日韩| 久久青草综合色| 精品酒店卫生间| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 国产老妇伦熟女老妇高清| 伊人亚洲综合成人网| 国产成人aa在线观看| 最新的欧美精品一区二区| 韩国高清视频一区二区三区| 日韩中字成人| videos熟女内射| 大码成人一级视频| 久久久欧美国产精品| 9色porny在线观看| 成年人免费黄色播放视频| 亚洲欧美日韩卡通动漫| 99视频精品全部免费 在线| 大香蕉97超碰在线| 亚洲精品自拍成人| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 国产色婷婷99| 国产视频首页在线观看| 久久鲁丝午夜福利片| 黄色毛片三级朝国网站| 日本免费在线观看一区| 丝袜脚勾引网站| 蜜臀久久99精品久久宅男| 999精品在线视频| 成人国产麻豆网| 丁香六月天网| 免费久久久久久久精品成人欧美视频 | 一区二区av电影网| 乱码一卡2卡4卡精品| 国产欧美另类精品又又久久亚洲欧美| av视频免费观看在线观看| 一级二级三级毛片免费看| 免费大片黄手机在线观看| 99久久综合免费| 九九在线视频观看精品| 熟女人妻精品中文字幕| 国产精品无大码| 久久久久久久国产电影| 国产成人免费观看mmmm| 一区二区日韩欧美中文字幕 | 丰满少妇做爰视频| 天天操日日干夜夜撸| 日韩一区二区三区影片| 又大又黄又爽视频免费| 精品久久久噜噜| 亚洲欧洲精品一区二区精品久久久 | 汤姆久久久久久久影院中文字幕| 女的被弄到高潮叫床怎么办| 一本—道久久a久久精品蜜桃钙片| 国产精品麻豆人妻色哟哟久久| 在线观看三级黄色| 亚洲精品一二三| 日日爽夜夜爽网站| 午夜福利在线观看免费完整高清在| a级毛片黄视频| 亚洲精品,欧美精品| 久久久久久久国产电影| 在线观看三级黄色| 午夜福利,免费看| 青青草视频在线视频观看| 亚洲精品日韩av片在线观看| 亚洲精品一区蜜桃| 交换朋友夫妻互换小说| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 亚洲精品久久午夜乱码| a级片在线免费高清观看视频| 亚洲内射少妇av| 国产一区有黄有色的免费视频| 亚洲av成人精品一二三区| 国产精品久久久久久久久免| 街头女战士在线观看网站| 欧美日韩av久久| 在线看a的网站| 亚洲av免费高清在线观看| 全区人妻精品视频| 一本大道久久a久久精品| 在线看a的网站| 欧美精品亚洲一区二区| 男的添女的下面高潮视频| 女性被躁到高潮视频| 国产免费一级a男人的天堂| 国产免费一区二区三区四区乱码| 老司机影院毛片| 人妻少妇偷人精品九色| 国产精品偷伦视频观看了| 三级国产精品片| 亚洲精品视频女| 在线亚洲精品国产二区图片欧美 | 丰满饥渴人妻一区二区三| 视频区图区小说| 黄色欧美视频在线观看| 色视频在线一区二区三区| 人体艺术视频欧美日本| 99热全是精品| 国产成人91sexporn| 两个人免费观看高清视频| 天美传媒精品一区二区| 久久人妻熟女aⅴ| 麻豆乱淫一区二区| 亚洲av成人精品一二三区| 久久狼人影院| 一区二区三区四区激情视频| 精品熟女少妇av免费看| 久久精品久久久久久久性| 九九久久精品国产亚洲av麻豆| 男女无遮挡免费网站观看| 久久久久久久精品精品| 亚洲欧美一区二区三区黑人 | 久久久久精品性色| 欧美日韩视频精品一区| 亚洲av成人精品一区久久| 欧美性感艳星| 欧美日韩一区二区视频在线观看视频在线| 麻豆精品久久久久久蜜桃| 亚洲国产精品专区欧美| 日韩一区二区视频免费看| 高清不卡的av网站| 欧美精品一区二区大全| 中文字幕免费在线视频6| 一边亲一边摸免费视频| 国产在线免费精品| 国产精品一区www在线观看| 精品国产乱码久久久久久小说| 成人综合一区亚洲| 国产av国产精品国产| 永久免费av网站大全| 女的被弄到高潮叫床怎么办| 国产国拍精品亚洲av在线观看| 日韩,欧美,国产一区二区三区| 我要看黄色一级片免费的| 欧美精品一区二区免费开放| 欧美精品高潮呻吟av久久| 国产成人av激情在线播放 | av卡一久久| 中文天堂在线官网| 内地一区二区视频在线| 亚洲精品国产色婷婷电影| 老司机影院毛片| 一边摸一边做爽爽视频免费| av免费在线看不卡| 久久久久久久久久成人| 五月天丁香电影| 国产一区二区在线观看日韩| 久久这里有精品视频免费| 啦啦啦中文免费视频观看日本| 成年美女黄网站色视频大全免费 | 女人精品久久久久毛片| 国产成人一区二区在线| av不卡在线播放| 免费观看a级毛片全部| 国产日韩欧美在线精品| 国精品久久久久久国模美| 人成视频在线观看免费观看| 自线自在国产av| 国产片内射在线| 99视频精品全部免费 在线| 久久精品久久久久久久性| 大码成人一级视频| 黄色一级大片看看| 亚洲欧美中文字幕日韩二区| 国产伦精品一区二区三区视频9| 午夜免费男女啪啪视频观看| 国产在视频线精品| 成人手机av| 精品一区二区三卡| 国产精品久久久久久av不卡| 亚洲精品日韩在线中文字幕| 在线观看人妻少妇| 免费观看性生交大片5| 男的添女的下面高潮视频| av线在线观看网站| 9色porny在线观看| 久久99一区二区三区| 蜜桃久久精品国产亚洲av| 精品久久久久久久久av| 国产高清不卡午夜福利| 妹子高潮喷水视频| 在线观看三级黄色| 久久久久久久久久久免费av| 国产精品国产三级国产专区5o| 日韩中文字幕视频在线看片| 成年人免费黄色播放视频| 狂野欧美激情性bbbbbb| 成人国产av品久久久| 国产精品女同一区二区软件| av.在线天堂| 午夜精品国产一区二区电影| 成人影院久久| 亚洲欧美一区二区三区黑人 | 国产69精品久久久久777片| 精品午夜福利在线看| 性高湖久久久久久久久免费观看| 麻豆乱淫一区二区| 在线观看三级黄色| 国产白丝娇喘喷水9色精品| 麻豆精品久久久久久蜜桃| 啦啦啦在线观看免费高清www| 欧美激情极品国产一区二区三区 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av.av天堂| 日本黄色日本黄色录像| 最近的中文字幕免费完整| 各种免费的搞黄视频| 少妇高潮的动态图| 久久99热6这里只有精品| 少妇的逼好多水| 人体艺术视频欧美日本| 久久久久久久大尺度免费视频| 国产亚洲最大av| 亚洲av日韩在线播放| 久久精品国产鲁丝片午夜精品|