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

    遙感影像區(qū)域面積快速計算并行算法研究

    2021-07-20 04:49:23楊靜靜馬駿
    計算機時代 2021年6期
    關鍵詞:并行計算遙感影像

    楊靜靜 馬駿

    摘 ?要: 為了有效地解決遙感應用中感興趣區(qū)域內有效遙感影像面積統(tǒng)計效率的問題,提出一種基于MPI(Message Passing Interface)的并行計算環(huán)境,采用經典數(shù)學定理鞋帶公式(Shoelace Formula)及計算機圖形學中向量積法計算多邊形面積。使用三角分割法對多邊形進行切割并行,判斷多邊形之間拓撲關系,獲取交集頂點表,使用鞋帶公式并行計算交多邊形面積。以1.96407的加速比有效提高遙感影像有效面積統(tǒng)計效率,加快網頁加載速度,提高用戶體驗度。

    關鍵詞: 遙感影像; 交多邊形面積; MPI; 并行計算; 鞋帶公式; 向量積

    中圖分類號:TP312 ? ? ? ? ?文獻標識碼:A ? ? 文章編號:1006-8228(2020)06-08-05

    Abstract: In order to effectively solve the problem of statistical efficiency of effective remote sensing image area in the region of interest in remote sensing applications, this paper proposes a parallel computing environment based on MPI (Message Passing Interface), which uses the classic mathematical theorem Shoelace Formula and Cross Product Algorithm in computer graphics to calculate the polygon area. Triangular segmentation is used to cut the polygons in parallel, determine the topological relationship between the polygons, obtain the intersection vertex data, and Shoelace Formula is used to calculate the area of the intersection polygons in parallel. With an acceleration ratio of 1.96407, it effectively improves the statistical efficiency of the effective area of remote sensing images, speeds up the loading of web pages, and improves user experience.

    Key words: remote sensing image; intersection polygon area; MPI; parallel computing; Shoelace Formula; Cross Product

    0 引言

    隨著衛(wèi)星遙感技術的迅速發(fā)展,采集的衛(wèi)星遙感影像爆炸式冪指數(shù)增加。面對巨量的遙感影像,滿足大范圍檢索條件要求的數(shù)據趨于上萬條,對判斷遙感影像間拓撲關系及計算遙感影像交并面積的計算效率提出挑戰(zhàn)。

    每一個遙感影像都有專屬的坐標位置屬性信息,故將計算多個遙感影像交并面積問題轉換為計算多個簡單多邊形交并面積問題。計算多邊形交并面積,首先需要獲取多邊形交并集頂點。判斷任意多邊形間拓撲關系不僅是計算幾何與計算機圖形學中的基本問題,也是GIS疊加分析的理論基礎。因此,對這一問題的研究無論是在理論上,還是在實踐上都有重要意義[4]。國內外針對多邊形交并判斷及交并面積計算問題,分別提出各類針對不同應用環(huán)境的多個多邊形拓撲關系判斷算法。如Shamos算法[5]、O' Rourke算法[6]用于處理凸多邊形,Weiler-Atherton算法[7]、Vatti算法[8]、Greiner-Hormann算法[9]、Sutherland-Hodgman算法[10]和M.Rivero提出的算法[11]能處理任意簡單多邊形的情況。針對這些經典算法在不同生活領域的應用,分別存在不同角度及不同層面的不足。對此,許多相關領域專家分別提出的不同的改進方法。

    齊東洲等提出優(yōu)化數(shù)據結構提高算法執(zhí)行效率的計算交點、將交點插入多邊形頂點序列、遍歷三個步驟的高效多邊形布爾計算方法[12]。陳濤提出了適用于凹多邊形裁剪的Cyrus-Beck改進算法[13]。劉勇奎等[14]、侯寶明等[15]和魏許青[16]等分別在Weiler算法的基礎上改進算法;劉勇奎等提出的算法采用單鏈表數(shù)據結構,減少了計算交點的時間;魏許青利用算法中的多邊形鏈表求出交/并集頂點表,求出交/并多邊形面積;侯寶明等采用接縫技術的算法大大簡化了運算過程。王結臣等提出基于掃描線和梯形分割思想的復雜多邊形裁剪算法[17];彭杰等提出基于交點排序思想的多邊形裁剪算法[18]等。結合以上多邊形裁切算法計算多邊形交/并面積的研究已經逐漸擴展。陳占龍等實現(xiàn)多核環(huán)境下基于Hilbert曲線空間劃分的多邊形并行合并[19];范俊甫等基于OGC簡單要素模型和Vatti算法,在集群MPI環(huán)境下對圖層級多邊形疊加合并[20];趙斯思等實現(xiàn)基于GPU加速的多邊形裁剪[21];高藝等提出基于GPU的柵格化多邊形相交面積算法GPURAS,并結合蒙特卡洛方法和遮擋查詢技術算法優(yōu)勢互補,提高計算效率[22]。上述研究學者主要在改進經典算法、經典算法結合優(yōu)勢互補方向深入研究,而在基于計算幾何中數(shù)學定理的MPI并行環(huán)境下計算多邊形交/并面積的研究較少。

    本文提出了一種基于MPI并行計算環(huán)境,鞋帶公式及向量積法計算多邊形面積。用三角分割法對多邊形進行切割并行,判斷多邊形之間拓撲關系,獲取交集頂點表,使用鞋帶公式并行計算交多邊形面積的并行算法研究。本文算法在實際應用中有效提高數(shù)據執(zhí)行效率,縮短程序執(zhí)行時間,提高用戶體驗度。在遙感影像有效區(qū)域面積統(tǒng)計研究方面具有一定的參考價值。

    1 多邊形交集面積計算算法

    1.1 多邊形相交判斷及交點坐標獲取

    為了計算方便借坐標原點為輔助點,視多邊形A為實體多邊形,則多邊形B為切割多邊形(首先要判斷多邊形A、B點排列方向順時針輸入,保證多邊形面積為正)。坐標原點與多邊形任意相鄰的兩個頂點構成一個三角形(如圖1多邊形三角分割),而三角形的面積可由三個頂點構成的兩個平面向量的外積求得。多邊形面積的計算和原點的選取沒有關系,通常選取原點坐標點o(0,0)或者多邊形的第一個點,這里選取原點坐標點。以切割多邊形B任意兩點與原點o(0,0)組成的三角形ocd,每條邊為直線向量切割實體多邊形A任意兩點與原點組成的三角形oab,用向量積法判斷多邊形A、B每條邊之間拓撲關系,獲取交點坐標集。向量積求三角形面積如圖2所示。

    1.2 多邊形面積計算

    1.2.1 鞋帶公式

    鞋帶公式[1-3],又稱“鞋帶算法”(Shoelace Formula,也稱為高斯面積公式),是一種數(shù)學算法,其用于計算簡單多邊形的面積,這些多邊形的頂點由平面中的笛卡爾坐標表示。對這些坐標進行叉乘計算,找到包圍多邊形的區(qū)域,并從周圍的多邊形中找到其中的多邊形區(qū)域。根據其交叉乘的過程像系鞋帶一樣,故稱為鞋帶公式。

    1.2.2 多邊形面積計算

    本文用鞋帶公式計算多邊形面積。 依據1.1節(jié)中獲取的交點坐標Pi,循環(huán)遍歷鞋帶算法計算兩三角形有向交面積之和Sum。

    傳統(tǒng)經典算法通過循環(huán)遍歷切割多邊形,向量積法多次循環(huán)判斷切割三角形兩兩之間拓撲關系判斷及獲取多邊形交點坐標,調用鞋帶算法計算交多邊形面積,即計算三角形之間有向交面積并統(tǒng)計所有切割三角形交面積之和作為多邊形有效交面積。傳統(tǒng)多邊形交面積計算程序數(shù)據計算串行執(zhí)行耗時較長,在遇到多邊形頂多較多的情況時會降低算法的執(zhí)行效率,影響算法應用效果。

    2 基于MPI并行計算交多邊形面積

    本文針對上述算法存在的缺陷,提出基于MPI信息傳遞接口并行環(huán)境計算多邊形面積,對數(shù)據分塊進行多線程并行計算,減少數(shù)據計算執(zhí)行等待時間,提高算法的執(zhí)行效率。

    2.1 MPI消息傳遞接口

    消息傳遞接口(Message Passing Interface,簡稱MPI),是一種跨語言的通訊協(xié)議,用于編寫并行計算機,支持點對點和廣播通信[24]。目前用于編寫并行計算機較為流行。MPI是一個并行計算消息傳遞接口標準,由MPI論壇(MPI Forum)推出,制定該標準的目的是提高并行程序的可移植性和開發(fā)效率[25]。MPI現(xiàn)在已經成為產業(yè)界廣泛支持的并行計算標準[26]。MPI主要消息傳遞接口如下。

    ⑴ Mpi_Init(); 初始化MPI環(huán)境,建立多個MPI進程之間的聯(lián)系,為后續(xù)通信做準備。

    ⑵ Mpi_Finalize(); 退出MPI環(huán)境。

    ⑶ Mpi_Comm_rank(MPI_COMM_WORLD,&rank);獲取當前進程在指定通信域中的編號,返回整型的錯誤值,將自身與其他程序區(qū)分。MPI_COMM_WORLD類型的通信域,標識參與計算的MPI進程組;&rank返回調用進程中的標識號。

    ⑷ Mpi_Comm_size(MPI_COMM_WORLD,&size);獲取指定通信域的進程個數(shù),確定自身完成任務比例。

    ⑸ MPI_Bcast(&n1,1,MPI_INT,0, MPI_COMM_WORLD);廣播數(shù)據集。

    ⑹ MPI_Reduce(&myResult,&result,1,MPI_DOUBLE, MPI_SUM,0,MPI_COMM_WORLD); 數(shù)據歸約,由進程0進行歸約,把每個進程計算出來的myResult進行相加(MPI_SUM),賦給result。

    ⑺ Mpi_send(buf,counter,datatype,dest,tag,comm):buf發(fā)送緩沖區(qū)的起始地址,可以是數(shù)組或結構指針;count:非負整數(shù),發(fā)送的數(shù)據個數(shù);datatype:發(fā)送數(shù)據的數(shù)據類型;dest:整型,目的的進程號;tag:整型,消息標志;comm:MPI進程組所在的通信域。

    ⑻ Mpi_recv(buf,count,datatype,source,tag,comm,status):source:整型,接收數(shù)據的來源,即發(fā)送數(shù)據進程的進程號; status:MPI_Status結構指針,返回狀態(tài)信息。

    2.2 基于MPI并行計算多邊形面積

    據算法使用背景本文使用組通信并行編程模型,利用消息傳遞方法實現(xiàn)任務間的并行。通過指定通信因子和組來對進程進行一種邏輯上劃分,通訊因子定義了進程組內或組間通訊的上下文。Foster方法下并行化設計步驟。

    ⑴ 劃分:數(shù)據分發(fā),0號進程讀入整個坐標數(shù)據向量,將獲取的數(shù)據分為comm_size塊(comm_size個核),并將數(shù)據發(fā)送給需要分量的其他進程。不同進程對收到的數(shù)據塊進行計算順序相鄰的兩個頂點與原點組成的三角形面積。

    ⑵ 通信:各線程內坐標數(shù)據塊的臨界坐標需要相互間通信,并發(fā)地執(zhí)行對數(shù)據的相交判斷及面積計算,期間多機有相互間的網絡通信開銷。

    ⑶ 聚集:多任務運行時,0號進程收集向量的所有分量,然后由0號進程聚集降低原始任務之間本需要的相互通信開銷。

    ⑷ 映射:聚合過的任務主要通過for循環(huán)部分通過調度對任務進行并行化計算,循環(huán)中的迭代需要按照迭代次數(shù)及核數(shù)進行輪播分配。

    如圖3組通信任務分配所示。源數(shù)據指定特定通信組,通信組內由多個處理機組成,通信組內通過MPI_Bcast廣播數(shù)據組,根據組內處理機數(shù)量對源數(shù)據進行任務分塊,組內處理機之間通過網絡傳送消息實現(xiàn)相互通信,實現(xiàn)任務間并行。MPI_Reduce規(guī)約各個進程任務塊計算結果。

    本文主要用到上一節(jié)中前6個消息傳遞接口。

    3 實驗與分析

    在Windows7 64位系統(tǒng)環(huán)境,8G運行內存,同樣機型4臺,一臺作為主進程(0號進程,負責數(shù)據讀入、分發(fā)、聚合),其他3臺機作為并行節(jié)點。

    實驗主要采用梯度驗證法設置一些測試梯度,以提高測試準確性。通常采用加速比S和效率E來測量并行計算效率,如等式⑸和等式⑹。

    其中,Ts為串行算法執(zhí)行耗費時間,Tp為并行算法執(zhí)行耗費時間,P為算法并行節(jié)點數(shù)。

    為了更加準確的驗證實驗結論,本文分兩種情況進行實驗驗證,一是節(jié)點數(shù)一定改變多邊形頂點數(shù);二是多邊形頂點數(shù)一定改變節(jié)點數(shù)。首先我們驗證節(jié)點數(shù)一定為2,多邊形頂點數(shù)分別為3,10,20,40,實驗結果如表1節(jié)點數(shù)一定改變多邊形頂點個數(shù)和圖4節(jié)點數(shù)一定改變多邊形頂點個數(shù)所示。

    當并行節(jié)點數(shù)一定時,改變串行、并行頂點數(shù)量,發(fā)現(xiàn)隨著頂點數(shù)的增加計算時間增大,當數(shù)據過小時CPU消費的時間會影響并行計算效率,并行效果不明顯,在頂點數(shù)一定時串行和并行計算時加速比和效率是呈線性增加趨勢。

    再次驗證頂點數(shù)一定為20,節(jié)點數(shù)分別為1,2,3,4,實驗結果如表2多邊形頂點個數(shù)一定改變節(jié)點數(shù)和圖5多邊形頂點個數(shù)一定改變節(jié)點數(shù)所示。

    當頂點數(shù)一定時,改變并行節(jié)點數(shù)量,發(fā)現(xiàn)隨著節(jié)點數(shù)的增加,計算時間增大,當數(shù)據過小時,CPU消費的時間會影響并行計算效率。節(jié)點數(shù)繼續(xù)增加,加速比增長緩慢且有下降趨勢,因為線程啟動需要開銷一定的時間。故隨著節(jié)點數(shù)的增加,并行程序的效率越來越低。

    綜上所述,本研究中當節(jié)點數(shù)為2時,并行計算效率最高,達到相對最佳并行計算節(jié)點。

    4 結論

    本研究主要為了解決多邊形頂點足夠多時,多邊形分割及三角形向量面積計算需要通過多層循環(huán)遍歷,程序串行執(zhí)行效率低下問題。本文提出基于MPI并行計算多邊形交面積,提高多邊形面積計算效率。在多邊形頂點增多時,本文并行計算多邊形面積效率明顯比串行計算速度快,相比串行方法本文并行計算方法在提高計算效率方面具有一定的研究價值。但是本文實驗利用消息傳遞方法實現(xiàn)多處理機處理任務并行,該方法需要網絡環(huán)境才能實現(xiàn)多處理機的并行,本次實驗僅在一種網絡環(huán)境下進行實驗,接下來可以嘗試換一種網絡帶寬進行并行測試,來驗證并行實驗的有效性。

    參考文獻(References):

    [1] Dahlke, Karl. "Shoelace Formula" (http://www.mathrefer-ence.com/la-det,shoe.html). Retrieved 9 June 2008.

    [2] Shoelace Theorem (http://www.artofproblemsolving.com/wiki/index.php?title=Shoelace_Theorem), Art of Problem Solving Wiki.

    [3] Younhee Lee, ?Woong Lim. Shoelace Formula: Connect-ing the Area of a Polygon with Vector Cross Product[J]. The Mathematics Teacher,2017.110(8):631-636

    [4] 朱雅音,王化文,萬豐等. 確定兩個任意簡單多邊形交_并_差的算法[J]. 計算機研究與發(fā)展, 2003.4:69-76

    [5] Preparata. Computational Geometry: An Introduction[J].?texts & monographs in computer science,1986.47(176).

    [6] O' Rourke. Computational Geometry in C[J]. Cambridge:Cam-bridge University Press,1985.37:128-129

    [7] Kevin Weiler , Peter Atherton. Hidden surface removal?using polygon area sorting[C]. Proceedings of the 4th Annual Conference on Computer Graphic and Interactive Techniques,1977.11:214-222

    [8] Bala R. Vatti. ?A generic solution to polygon clipping[J].?Communication of the ACM,1992.35:56-63

    [9] Hoffmann CM , The problems of accuracy and robustness?in geometric computations[J]. IEEE Computer,1989.22:31-42

    [10] 蘇誠,韓俊剛.Sutherland-Hodgman裁剪算法的改進[J].西安郵電學院學報,2013.18(3):80-82

    [11] Rivero, F R Feito. Boolean operations on general planarpolygons[M].Computers& Graphics,2000.24(6):881-898

    [12] 齊東洲,吳敏.高效的多邊形布爾計算方法[J].計算機應用,2014.S2:85-89

    [13] 陳濤.適用于凹多邊形的Cyrus-Beck改進算法[J].計算機科學,2006.33(12):217-220

    [14] 劉永奎,高云,黃有群.一個有效的多邊形裁剪算法[J].軟件學報,2003.14(4):845-856

    [15] 侯寶明,劉雪娜.任意多邊形區(qū)域交的有效算法[J],計算機輔助工程,2009.18(2):73-76

    [16] 魏許青.計算多邊形交集_并集面積的算法[J].計算機工程與科學,2007.12:89-90

    [17] 王結臣,沈定濤,陳焱明,等.一種有效的復雜多邊形裁剪算法[J].武漢大學學報(信息科學版),2010.35(3):369-372

    [18] 彭杰,劉南,唐遠彬等.一種基于交點排序的高效多邊形裁剪算法[J].浙江大學學報(理學版),2012.39(1):107-111

    [19] 陳占龍,吳亮,劉煥煥.多核環(huán)境下Hilbert曲線劃分簡單要素多邊形合并算法[J].計算機應用研究,2012.29(7):2747-2750

    [20] 范俊甫,馬廷,周成虎等.基于集群MPI的圖層級多邊形并行合并算法[J].地球信息科學學報,2014.16(4): 15-21

    [21] 趙斯思,周成虎.GPU加速的多邊形疊加分析[J]. 地理科學進展,2013.32(1):114-120

    [22] 高藝,羅健欣,裘杭萍等.基于GPU的任意多邊形相交面積計算方法_高藝[J].測繪工程,2017.26(12):58-62

    [23] Meister, A. L. F., Generalia de genesi figurarum planarumt inde pendentibus earum affectionibus, https://books.google.com/books?id=fOE_AAAAcAAJ,1769,Nov.Com. G?tt. (in Latin), 1:144

    [24] https://baike.baidu.com/item/MPI/15277241?fr=aladdin

    [25] 李久楷,朱俊,寧交賢.MPI并行計算性能的研究[J].四川大學學報(自然科學版),2009.33(5):496-499

    [26] 申煥,石曉春,邱宏華.基于MPI的海量遙感影像并行處理技術探析[J].全球定位系統(tǒng),2012.37(6):73-76

    猜你喜歡
    并行計算遙感影像
    基于自適應線程束的GPU并行粒子群優(yōu)化算法
    遙感影像資料在海圖制圖中的應用研究
    航海(2016年6期)2017-01-09 11:28:24
    云計算中MapReduce分布式并行處理框架的研究與搭建
    矩陣向量相乘的并行算法分析
    大經貿(2016年9期)2016-11-16 16:25:33
    并行硬件簡介
    基于GPU的超聲場仿真成像平臺
    遙感數(shù)字圖像處理課程實驗綜述
    科技視界(2016年13期)2016-06-13 12:01:14
    基于Matlab的遙感圖像IHS小波融合算法的并行化設計
    科技視界(2016年11期)2016-05-23 08:13:35
    臥龍湖濕地遙感調查和演變分析
    高分遙感影像中道路信息提取方法綜述
    免费观看a级毛片全部| 激情视频va一区二区三区| freevideosex欧美| 狠狠婷婷综合久久久久久88av| 巨乳人妻的诱惑在线观看| 18禁观看日本| 久久人人97超碰香蕉20202| 亚洲av在线观看美女高潮| www.自偷自拍.com| 国产深夜福利视频在线观看| 女人被躁到高潮嗷嗷叫费观| 在线亚洲精品国产二区图片欧美| 天天躁夜夜躁狠狠躁躁| 飞空精品影院首页| 免费看av在线观看网站| 天天躁夜夜躁狠狠躁躁| videos熟女内射| 国产一区有黄有色的免费视频| 国产成人精品福利久久| 久久精品国产自在天天线| 美女xxoo啪啪120秒动态图| 午夜激情av网站| 午夜福利视频在线观看免费| 国产免费现黄频在线看| 国产熟女午夜一区二区三区| 久久久久久久久免费视频了| 亚洲成人一二三区av| 欧美变态另类bdsm刘玥| 电影成人av| 在线观看三级黄色| 欧美国产精品va在线观看不卡| 日日啪夜夜爽| 亚洲四区av| 伊人久久大香线蕉亚洲五| 中文字幕人妻丝袜一区二区 | 亚洲天堂av无毛| www.精华液| 亚洲一区中文字幕在线| 成人午夜精彩视频在线观看| 视频在线观看一区二区三区| 成人18禁高潮啪啪吃奶动态图| 9热在线视频观看99| 天天躁狠狠躁夜夜躁狠狠躁| 最近最新中文字幕大全免费视频 | 两性夫妻黄色片| 丰满乱子伦码专区| 水蜜桃什么品种好| 欧美变态另类bdsm刘玥| 日本av免费视频播放| 在线观看国产h片| 丝袜脚勾引网站| 日本欧美视频一区| 电影成人av| 亚洲欧洲日产国产| 亚洲三区欧美一区| 亚洲精品美女久久av网站| 波多野结衣一区麻豆| 9色porny在线观看| 久久久久久人人人人人| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 日日摸夜夜添夜夜爱| 丝袜在线中文字幕| 亚洲精品美女久久久久99蜜臀 | 久久精品国产亚洲av高清一级| 亚洲av.av天堂| 精品卡一卡二卡四卡免费| 狠狠精品人妻久久久久久综合| 亚洲精品美女久久av网站| 一区在线观看完整版| 国产综合精华液| 老司机影院成人| 在线观看免费高清a一片| 亚洲av日韩在线播放| 叶爱在线成人免费视频播放| 亚洲精品视频女| 国产极品粉嫩免费观看在线| 美女国产视频在线观看| www.熟女人妻精品国产| 高清不卡的av网站| 午夜日韩欧美国产| 黄色毛片三级朝国网站| 久久综合国产亚洲精品| 午夜日韩欧美国产| 美女福利国产在线| 亚洲国产欧美日韩在线播放| 午夜免费观看性视频| 亚洲国产毛片av蜜桃av| kizo精华| 欧美精品一区二区大全| 少妇人妻久久综合中文| 一区二区三区四区激情视频| 黄色视频在线播放观看不卡| 亚洲国产成人一精品久久久| 韩国高清视频一区二区三区| 波多野结衣av一区二区av| 人妻一区二区av| 久久国产精品大桥未久av| 日本午夜av视频| 中国国产av一级| 五月天丁香电影| 精品国产超薄肉色丝袜足j| 亚洲少妇的诱惑av| 欧美人与善性xxx| 欧美xxⅹ黑人| 这个男人来自地球电影免费观看 | 国产 一区精品| 一级毛片电影观看| 侵犯人妻中文字幕一二三四区| 国产免费又黄又爽又色| 国产av码专区亚洲av| 美女国产视频在线观看| 热re99久久国产66热| 亚洲国产最新在线播放| 夜夜骑夜夜射夜夜干| 国产欧美日韩一区二区三区在线| 赤兔流量卡办理| 久久精品国产亚洲av高清一级| 日韩精品有码人妻一区| 黑人猛操日本美女一级片| 老司机影院成人| 久久久久久久亚洲中文字幕| 国产男女超爽视频在线观看| 久久这里只有精品19| 精品第一国产精品| 高清不卡的av网站| 妹子高潮喷水视频| 久久久久久久亚洲中文字幕| 人体艺术视频欧美日本| 99香蕉大伊视频| 欧美国产精品一级二级三级| 欧美 亚洲 国产 日韩一| 国产无遮挡羞羞视频在线观看| 老汉色∧v一级毛片| 久久久欧美国产精品| 精品人妻熟女毛片av久久网站| 国产在线免费精品| 中文字幕人妻丝袜一区二区 | 久久婷婷青草| 亚洲国产日韩一区二区| 汤姆久久久久久久影院中文字幕| 色吧在线观看| 国产av国产精品国产| 秋霞在线观看毛片| 日本欧美视频一区| 国产男女超爽视频在线观看| 亚洲第一区二区三区不卡| 久久99蜜桃精品久久| 日韩中文字幕视频在线看片| 免费观看无遮挡的男女| 少妇熟女欧美另类| 国产淫语在线视频| 天天躁夜夜躁狠狠久久av| 有码 亚洲区| 中国国产av一级| 男人爽女人下面视频在线观看| 亚洲精品久久成人aⅴ小说| 日韩av在线免费看完整版不卡| 亚洲精品久久午夜乱码| 午夜影院在线不卡| 欧美日韩精品成人综合77777| 色播在线永久视频| 亚洲三级黄色毛片| 日韩伦理黄色片| 99九九在线精品视频| 在线观看一区二区三区激情| 国产精品无大码| 91精品国产国语对白视频| 久久久a久久爽久久v久久| 母亲3免费完整高清在线观看 | 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区 | 妹子高潮喷水视频| 久久99一区二区三区| 久久人人爽av亚洲精品天堂| 女人高潮潮喷娇喘18禁视频| 国产极品天堂在线| 啦啦啦视频在线资源免费观看| 国产精品久久久av美女十八| 国产 一区精品| 日韩免费高清中文字幕av| 在线观看免费视频网站a站| 久久99蜜桃精品久久| 王馨瑶露胸无遮挡在线观看| 18禁国产床啪视频网站| 最近最新中文字幕免费大全7| 国产精品99久久99久久久不卡 | 欧美人与性动交α欧美精品济南到 | 丝袜脚勾引网站| 不卡视频在线观看欧美| 男女边吃奶边做爰视频| 亚洲av在线观看美女高潮| 中文天堂在线官网| 亚洲国产最新在线播放| 国产探花极品一区二区| 日韩av免费高清视频| 另类精品久久| 国产精品 国内视频| 亚洲精品av麻豆狂野| 亚洲国产毛片av蜜桃av| 看非洲黑人一级黄片| 超碰成人久久| 九草在线视频观看| 男女高潮啪啪啪动态图| 涩涩av久久男人的天堂| 青草久久国产| av不卡在线播放| 亚洲久久久国产精品| 亚洲av.av天堂| 人妻系列 视频| 成人亚洲精品一区在线观看| 欧美另类一区| 国产成人精品一,二区| 只有这里有精品99| 老司机亚洲免费影院| 不卡视频在线观看欧美| 男女国产视频网站| 国产又爽黄色视频| 人人妻人人澡人人看| 中文乱码字字幕精品一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 777久久人妻少妇嫩草av网站| 岛国毛片在线播放| 激情五月婷婷亚洲| tube8黄色片| 91国产中文字幕| 一本久久精品| 欧美黄色片欧美黄色片| 国产精品欧美亚洲77777| 国产爽快片一区二区三区| 伊人亚洲综合成人网| 国产精品久久久久久av不卡| a级片在线免费高清观看视频| 极品人妻少妇av视频| 国产成人aa在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲日产国产| 国产精品久久久av美女十八| 欧美人与性动交α欧美精品济南到 | 天天躁夜夜躁狠狠躁躁| 成年人免费黄色播放视频| 亚洲精品自拍成人| 日本vs欧美在线观看视频| 毛片一级片免费看久久久久| 欧美成人午夜免费资源| 90打野战视频偷拍视频| 只有这里有精品99| 男女高潮啪啪啪动态图| av女优亚洲男人天堂| 在线观看一区二区三区激情| 国产熟女欧美一区二区| 国产成人免费观看mmmm| 又黄又粗又硬又大视频| 亚洲国产看品久久| 国产成人精品久久二区二区91 | 成人黄色视频免费在线看| 黄色怎么调成土黄色| 国产一级毛片在线| 国产成人精品福利久久| 一区二区av电影网| 99久久综合免费| 免费少妇av软件| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区三区| 久久人人爽人人片av| 国产av精品麻豆| 色视频在线一区二区三区| 国产精品香港三级国产av潘金莲 | 国产精品国产三级专区第一集| 天天躁日日躁夜夜躁夜夜| 亚洲第一青青草原| 国产精品欧美亚洲77777| 丁香六月天网| 国产成人91sexporn| 老汉色∧v一级毛片| 男女啪啪激烈高潮av片| www日本在线高清视频| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 欧美黄色片欧美黄色片| 久久久久人妻精品一区果冻| 巨乳人妻的诱惑在线观看| 人人妻人人澡人人看| 韩国高清视频一区二区三区| 成人毛片a级毛片在线播放| 99久久综合免费| 99国产精品免费福利视频| 国产亚洲欧美精品永久| 久久久久网色| 久久国产精品大桥未久av| 午夜免费鲁丝| 亚洲第一青青草原| 久久久久人妻精品一区果冻| 国产精品国产av在线观看| 国产精品人妻久久久影院| 久久久久久久国产电影| 中文字幕人妻熟女乱码| 妹子高潮喷水视频| 久久精品亚洲av国产电影网| 久久精品夜色国产| 精品福利永久在线观看| 欧美国产精品va在线观看不卡| 久久女婷五月综合色啪小说| 亚洲精品av麻豆狂野| 熟女av电影| 啦啦啦在线观看免费高清www| 伊人久久大香线蕉亚洲五| 日本av免费视频播放| 一级a爱视频在线免费观看| 丝袜美腿诱惑在线| 亚洲婷婷狠狠爱综合网| 亚洲中文av在线| 色视频在线一区二区三区| 国产精品国产av在线观看| 午夜激情久久久久久久| 久久久久久久亚洲中文字幕| 少妇的丰满在线观看| 九九爱精品视频在线观看| 日本色播在线视频| 亚洲av综合色区一区| 精品人妻偷拍中文字幕| 777久久人妻少妇嫩草av网站| 日本-黄色视频高清免费观看| 午夜福利视频精品| 美女高潮到喷水免费观看| 在线观看美女被高潮喷水网站| 国产精品国产av在线观看| 免费久久久久久久精品成人欧美视频| 国产又爽黄色视频| av网站在线播放免费| 国产日韩欧美亚洲二区| 色吧在线观看| tube8黄色片| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 亚洲综合色惰| 国产亚洲av片在线观看秒播厂| 欧美日韩一级在线毛片| 天天躁夜夜躁狠狠躁躁| 一区二区av电影网| 伊人久久国产一区二区| av免费在线看不卡| 国产亚洲av片在线观看秒播厂| 日韩不卡一区二区三区视频在线| 伊人久久大香线蕉亚洲五| 日韩中字成人| 纯流量卡能插随身wifi吗| 国产 一区精品| 热re99久久精品国产66热6| 麻豆av在线久日| 最近2019中文字幕mv第一页| 国产精品 国内视频| 国产精品久久久久久精品电影小说| 午夜福利在线观看免费完整高清在| 亚洲综合精品二区| 国产有黄有色有爽视频| 久久国内精品自在自线图片| 高清黄色对白视频在线免费看| 午夜91福利影院| 国产白丝娇喘喷水9色精品| 亚洲精品久久午夜乱码| 精品酒店卫生间| 亚洲精品久久午夜乱码| 国产一区二区激情短视频 | 久久精品熟女亚洲av麻豆精品| 考比视频在线观看| 久久狼人影院| 只有这里有精品99| 亚洲精品av麻豆狂野| 精品99又大又爽又粗少妇毛片| 国产白丝娇喘喷水9色精品| 日日撸夜夜添| 妹子高潮喷水视频| 色94色欧美一区二区| 亚洲av福利一区| 国产亚洲一区二区精品| 日本免费在线观看一区| 成人漫画全彩无遮挡| 涩涩av久久男人的天堂| 国产极品天堂在线| 国产在视频线精品| 尾随美女入室| 亚洲av福利一区| 18禁动态无遮挡网站| a级毛片在线看网站| 久久精品国产综合久久久| 国产成人一区二区在线| 大片电影免费在线观看免费| 亚洲国产成人一精品久久久| 母亲3免费完整高清在线观看 | 少妇人妻精品综合一区二区| 2022亚洲国产成人精品| 久久久精品区二区三区| 波多野结衣一区麻豆| 18禁动态无遮挡网站| 免费在线观看黄色视频的| 国产爽快片一区二区三区| 国产精品秋霞免费鲁丝片| 日日啪夜夜爽| 又粗又硬又长又爽又黄的视频| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲日产国产| 欧美黄色片欧美黄色片| 啦啦啦视频在线资源免费观看| 久久精品人人爽人人爽视色| 午夜91福利影院| 中文字幕人妻丝袜制服| 国产老妇伦熟女老妇高清| 母亲3免费完整高清在线观看 | 寂寞人妻少妇视频99o| 国产在视频线精品| 美女大奶头黄色视频| 97人妻天天添夜夜摸| 一级片免费观看大全| 日韩一区二区三区影片| 99香蕉大伊视频| 国产深夜福利视频在线观看| 日韩中文字幕视频在线看片| av视频免费观看在线观看| 成人18禁高潮啪啪吃奶动态图| 99re6热这里在线精品视频| av在线播放精品| 亚洲经典国产精华液单| 99久国产av精品国产电影| 男的添女的下面高潮视频| 亚洲欧美精品综合一区二区三区 | 欧美国产精品va在线观看不卡| 人成视频在线观看免费观看| 秋霞在线观看毛片| 最黄视频免费看| 伊人久久大香线蕉亚洲五| 七月丁香在线播放| 日韩一区二区视频免费看| 一级毛片电影观看| 国产在线视频一区二区| 一区二区三区乱码不卡18| 最近最新中文字幕大全免费视频 | 2021少妇久久久久久久久久久| 欧美精品一区二区大全| 在线观看www视频免费| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 国产av一区二区精品久久| 青春草亚洲视频在线观看| 高清视频免费观看一区二区| 色播在线永久视频| 国产女主播在线喷水免费视频网站| 亚洲精品成人av观看孕妇| 中文欧美无线码| 女人被躁到高潮嗷嗷叫费观| 国产精品亚洲av一区麻豆 | 最黄视频免费看| 日韩制服丝袜自拍偷拍| 一级片'在线观看视频| 午夜福利一区二区在线看| 国产综合精华液| 久久久国产精品麻豆| 欧美 亚洲 国产 日韩一| 1024视频免费在线观看| 纯流量卡能插随身wifi吗| 女人精品久久久久毛片| 成人午夜精彩视频在线观看| 久久青草综合色| 成人亚洲精品一区在线观看| 丝袜在线中文字幕| 日韩伦理黄色片| 午夜福利视频在线观看免费| 亚洲精品国产色婷婷电影| 久久久久国产网址| 国产在视频线精品| 婷婷色麻豆天堂久久| 久热久热在线精品观看| 18禁裸乳无遮挡动漫免费视频| 成年女人在线观看亚洲视频| 精品少妇久久久久久888优播| 日韩熟女老妇一区二区性免费视频| 人妻人人澡人人爽人人| 视频在线观看一区二区三区| 免费黄网站久久成人精品| 亚洲国产欧美在线一区| 国产成人一区二区在线| 天天躁日日躁夜夜躁夜夜| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 七月丁香在线播放| 18+在线观看网站| 最近中文字幕高清免费大全6| 老司机影院毛片| 水蜜桃什么品种好| 99热国产这里只有精品6| 欧美精品高潮呻吟av久久| av天堂久久9| 久久久国产精品麻豆| 亚洲国产色片| 如日韩欧美国产精品一区二区三区| 女人高潮潮喷娇喘18禁视频| 精品一区在线观看国产| 国产成人精品福利久久| 亚洲国产精品成人久久小说| 国产精品无大码| 国产精品蜜桃在线观看| av免费观看日本| 18禁动态无遮挡网站| 少妇的丰满在线观看| 国产午夜精品一二区理论片| 卡戴珊不雅视频在线播放| 国产成人精品久久二区二区91 | 春色校园在线视频观看| 日韩av在线免费看完整版不卡| 亚洲国产成人一精品久久久| 性少妇av在线| tube8黄色片| 亚洲国产欧美在线一区| 91精品三级在线观看| 婷婷色综合www| 母亲3免费完整高清在线观看 | 国产极品天堂在线| 亚洲少妇的诱惑av| 成人影院久久| 国产精品一区二区在线观看99| a 毛片基地| 69精品国产乱码久久久| 久久久久久久亚洲中文字幕| 欧美日韩视频高清一区二区三区二| 精品国产露脸久久av麻豆| 久久人人爽人人片av| 亚洲国产av影院在线观看| 国产av一区二区精品久久| 久久韩国三级中文字幕| a级毛片在线看网站| 久久久久久人妻| 色哟哟·www| 伦理电影大哥的女人| 亚洲精品日本国产第一区| 亚洲天堂av无毛| 国产又爽黄色视频| av国产精品久久久久影院| 黄频高清免费视频| 亚洲精品中文字幕在线视频| 国产精品亚洲av一区麻豆 | 如日韩欧美国产精品一区二区三区| 午夜福利视频在线观看免费| 在线观看美女被高潮喷水网站| av.在线天堂| 国产精品久久久久久精品古装| 一级毛片 在线播放| 亚洲国产精品999| 咕卡用的链子| 中文欧美无线码| 在线观看一区二区三区激情| 亚洲精品久久午夜乱码| 人成视频在线观看免费观看| 免费黄网站久久成人精品| 一级毛片电影观看| 最近中文字幕2019免费版| 日韩中文字幕欧美一区二区 | 亚洲人成77777在线视频| 午夜激情av网站| 久久久精品国产亚洲av高清涩受| 免费看不卡的av| 国产免费一区二区三区四区乱码| 亚洲人成电影观看| 国产精品麻豆人妻色哟哟久久| 日韩成人av中文字幕在线观看| 成年美女黄网站色视频大全免费| 欧美 日韩 精品 国产| 欧美精品国产亚洲| 女性生殖器流出的白浆| 色婷婷久久久亚洲欧美| 一级,二级,三级黄色视频| 激情视频va一区二区三区| 婷婷色综合大香蕉| 精品国产超薄肉色丝袜足j| 五月伊人婷婷丁香| 丝袜在线中文字幕| 最黄视频免费看| 一级爰片在线观看| 久久久精品94久久精品| 久久亚洲国产成人精品v| 日本av手机在线免费观看| 久久久国产精品麻豆| 日本欧美国产在线视频| 国产黄色免费在线视频| 999精品在线视频| 精品久久久久久电影网| 自线自在国产av| 成人国产av品久久久| www.熟女人妻精品国产| 欧美+日韩+精品| 国产精品久久久久久av不卡| 久久人妻熟女aⅴ| 最近2019中文字幕mv第一页| 香蕉国产在线看| 亚洲av在线观看美女高潮| 久久久久国产一级毛片高清牌| 考比视频在线观看| xxxhd国产人妻xxx| 人体艺术视频欧美日本| 黄色视频在线播放观看不卡| 国产一区二区三区av在线| 久久精品国产鲁丝片午夜精品| av网站在线播放免费| 久久av网站| 大片免费播放器 马上看| 男女免费视频国产| 三级国产精品片| 日韩精品有码人妻一区| 考比视频在线观看| 国产精品久久久久成人av| 五月伊人婷婷丁香| 国产不卡av网站在线观看| 在现免费观看毛片|