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

    一種計(jì)算植物冠層光分布的并行仿真系統(tǒng)

    2018-03-13 22:26:49謝文昊
    現(xiàn)代電子技術(shù) 2018年6期
    關(guān)鍵詞:客戶機(jī)

    謝文昊

    摘 要: 光照對(duì)植物冠層生長(zhǎng)、繁殖的影響是巨大的,但在仿真植物冠層光照分布的過(guò)程中由于需要處理的數(shù)據(jù)量很大,如果利用單處理器計(jì)算機(jī)計(jì)算,會(huì)導(dǎo)致計(jì)算耗時(shí)過(guò)長(zhǎng),甚至無(wú)法得到高精度的仿真結(jié)果。為了提高計(jì)算速度,將用于計(jì)算植物冠層光分布的深度緩存排序算法進(jìn)行并行化處理,移植到大型多處理器計(jì)算機(jī)上進(jìn)行計(jì)算,并通過(guò)在普通計(jì)算機(jī)上開發(fā)客戶端對(duì)仿真計(jì)算進(jìn)行控制。設(shè)計(jì)實(shí)現(xiàn)一套基于客戶機(jī)?服務(wù)器模式的仿真系統(tǒng)。應(yīng)用該系統(tǒng)對(duì)不同生長(zhǎng)時(shí)期的水稻進(jìn)行光照仿真,得到了效果圖。從而得出結(jié)論該系統(tǒng)在分析植物生長(zhǎng)趨勢(shì)模擬全天冠層光照變化方面有著很大的應(yīng)用前景。

    關(guān)鍵詞: 植物冠層; 光照分布; 并行仿真; 客戶機(jī)?服務(wù)器模式; MPI并行模式; 深度緩存排序算法

    中圖分類號(hào): TN212?34; TP181 文獻(xiàn)標(biāo)識(shí)碼: A 文章編號(hào): 1004?373X(2018)06?0150?04

    Abstract: Illumination has a great influence on the growth and propagation of plant canopy. As a large amount of data needs to be processed during the simulation process of plant canopy illumination distribution, the time consumption of calculation can be very long, and the high?precision simulation results cannot even be obtained if the single processor computer is used for calculation. To improve the calculation speed, the parallel processing is performed for the Z?buffer sorting algorithm which is used for calculation of plant canopy illumination distribution and then transplanted to the large multiprocessor computer for calculation. The simulation calculation is controlled by developing the client on the ordinary computer. A simulation system based on client?server mode was designed and implemented. The system was applied to illumination simulation for paddies in different grown periods and the effect diagram was obtained. It can be concluded that the system has a great application prospect in plant growth tendency analysis and full?day canopy illumination variation simulation.

    Keywords: plant canopy; illumination distribution; parallel simulation; client?server mode; MPI parallel mode; Z?buffer sorting algorithm

    0 引 言

    植物生長(zhǎng)過(guò)程離不開光照的影響,因此計(jì)算光照在植物冠層內(nèi)的分布情況就很有必要[1]。隨著計(jì)算機(jī)圖形學(xué)的發(fā)展和應(yīng)用,目前常用的模擬方法有:輻射度方法[2]、光線跟蹤方法、深度緩存排序算法等[3]。其中深度緩存排序算法主要是先分別計(jì)算植物冠層在各方向光線直射下由于相互遮擋而產(chǎn)生的葉片對(duì)光輻射的截獲情況再進(jìn)行多方向疊加從而得到仿真結(jié)果。該方法雖然忽略了光的反射透射作用,但對(duì)于主要表現(xiàn)為光線直射部分的光合有效輻射分布的模擬有著很高的精度,因此被廣泛使用。另外,各方向上的光分布計(jì)算相互獨(dú)立,也很適合算法并行化處理。

    近年來(lái),計(jì)算機(jī)硬件技術(shù)不斷創(chuàng)新發(fā)展,高性能的多處理器計(jì)算機(jī)已經(jīng)成為科學(xué)計(jì)算的主要工具之一[4]。MPI并行模式也已作為一種普遍流行的并行編程規(guī)范被人們廣泛接受。并行計(jì)算是相對(duì)于串行計(jì)算而言的,并行計(jì)算是指在同一時(shí)間間隔內(nèi)增加操作系統(tǒng)進(jìn)程數(shù)量,利用多臺(tái)計(jì)算機(jī)共同實(shí)現(xiàn)同一個(gè)任務(wù)。MPI(Message Passing Interface)可以將一個(gè)大型任務(wù)分解成多個(gè)進(jìn)程,再把這些進(jìn)程分配給多臺(tái)計(jì)算機(jī)進(jìn)行并行處理,通過(guò)計(jì)算結(jié)果的綜合來(lái)得到最終結(jié)果[5]。本文將傳統(tǒng)串行實(shí)現(xiàn)的用于計(jì)算植物冠層光分布的深度緩存排序算法并行化處理,并移植到高性能多處理器計(jì)算機(jī)上運(yùn)行,建立服務(wù)器程序。因?yàn)槠胀ㄓ脩艉茈y直接操作高性能計(jì)算機(jī),所以相應(yīng)地設(shè)計(jì)開發(fā)了安裝在普通電腦上用于控制服務(wù)器的客戶端。從而建立了一套客戶端?服務(wù)器模式的植物冠層光分布的并行仿真計(jì)算系統(tǒng)。

    1 并行算法

    1.1 深度緩存排序算法原理

    深度緩存(Z?buffer)[6]最早由Catmull提出,是一種典型的、最簡(jiǎn)單的圖像空間面消隱算法,但所需的存儲(chǔ)容量較大。該算法的原理是,先將待處理的景物表面上的采樣點(diǎn)變換到圖像空間(屏幕坐標(biāo)系),計(jì)算其深度值,并根據(jù)采樣點(diǎn)在屏幕上的投影位置,將其深度值與已存儲(chǔ)在Z緩存器中相應(yīng)像素處的原可見點(diǎn)的深度值進(jìn)行比較。如果新的采樣點(diǎn)的深度(即Z軸坐標(biāo)值)大于原可見點(diǎn)的深度,表明新的采樣點(diǎn)遮住了原可見點(diǎn),則用該采樣點(diǎn)處的顏色值更新幀緩存器中相應(yīng)像素的顏色值,同時(shí)用其深度值更新Z緩存器中的深度值;否則,不做更改。深度緩存算法是在為物件進(jìn)行著色時(shí),執(zhí)行“隱藏面消除”工作的一項(xiàng)技術(shù),所以隱藏物件背后的部分就不會(huì)被顯示出來(lái)。 在3D環(huán)境中每個(gè)像素中會(huì)利用一組數(shù)據(jù)資料來(lái)定義像素在顯示時(shí)的縱深度。Z?buffer所用的位數(shù)越高,則代表該顯示卡所提供的物件縱深感也越精確。endprint

    當(dāng)忽略光的反射、透射作用,只考慮光的直射影響時(shí),植株結(jié)構(gòu)相互遮擋是植物冠層內(nèi)的光分布差異情況產(chǎn)生的主要原因,所以是可以應(yīng)用深度緩存算法進(jìn)行計(jì)算的。文獻(xiàn)[7]中給出了在計(jì)算植物冠層光分布時(shí)的原模型,如圖1所示。先確定一塊矩形的植物冠層區(qū)域?yàn)閷?shí)際觀測(cè)區(qū)域,假設(shè)該實(shí)測(cè)區(qū)域周圍存在著與實(shí)測(cè)區(qū)域冠層結(jié)構(gòu)完全相同的植物冠層,這樣假設(shè)是為了避免其周圍區(qū)域(稱其為復(fù)制域)的冠層結(jié)構(gòu)因?yàn)檎趽醵鴮?duì)區(qū)域內(nèi)的光分布產(chǎn)生影響,見圖1。

    此外,為了方便投影,該方法還要在三維空間坐標(biāo)系中對(duì)葉面進(jìn)行劃分,將其分割成很小的空間三角形面元,在如圖1所示的光線照射下來(lái)進(jìn)行投影[7]。

    該算法在投影過(guò)程中先在冠層正上方建立光照平面,將它作為全部入射光線的來(lái)源,并將光照平面進(jìn)行二維網(wǎng)格化[8],在光照平面二維網(wǎng)格化的基礎(chǔ)上假定每個(gè)網(wǎng)格上有一條入射光線。接下來(lái),由光線方向?qū)⑷~片面元投影到上述的光照平面上。對(duì)于類似于圖1中的A這樣的投影不在投影平面的小面元,可以用復(fù)制域中和A位置相同的面元B的投影來(lái)代替它。最后使用Z?buffer深度排序法計(jì)算投影后落在同一個(gè)小網(wǎng)格的小面元序列的深度緩存值,以便找到距離光照平面最短的小面元。這樣可以通過(guò)此面元光照的分布情況,進(jìn)而統(tǒng)計(jì)出光照在植物冠層中的三維分布[9]。

    上述投影過(guò)程中,面元的投影算法可以通過(guò)下式得到。設(shè)實(shí)測(cè)區(qū)域內(nèi)任意一點(diǎn)[(x,y,z)]投影后的坐標(biāo)就可以表示為:

    [px=x+(ytop-y)?cos(SunAzimuth)tan(SunH)pz=z+(ytop-y)?sin(SunAzimuth)tan(SunH)]

    式中:采用右手坐標(biāo)系,x軸向正南,y軸向正上,z軸向正西;[px],[pz]為[x,z]方向上投影后的坐標(biāo)值;[ytop]為實(shí)測(cè)區(qū)域高度;SunAzimuth為太陽(yáng)方位角,以正南為0°,向西為正,向東為負(fù);SunH為太陽(yáng)高度角,即太陽(yáng)與地平面的夾角。

    1.2 算法的并行化

    仿真野外條件下植物冠層光照的分布情況,就要在各方向上先分別對(duì)植物冠層進(jìn)行深度緩存排序來(lái)計(jì)算單一方向光入射時(shí)植物冠層光分布結(jié)果,再對(duì)各方向計(jì)算結(jié)果進(jìn)行累加。單處理器計(jì)算機(jī)只能依次串行的計(jì)算每個(gè)方向上的光照情況再將結(jié)果相加,工作量大、耗時(shí)長(zhǎng)。另外又受到實(shí)際內(nèi)存限制常常得不到令人滿意的高精度仿真結(jié)果。而多處理器高性能計(jì)算機(jī)則可以同時(shí)并發(fā)地執(zhí)行多項(xiàng)任務(wù),所以并行化算法將每個(gè)方向上的光照分布計(jì)算工作平均分布到高性能計(jì)算機(jī)不同的處理器中并發(fā)執(zhí)行,并將結(jié)果匯總到共用的數(shù)據(jù)容器中,下面把像這樣的負(fù)責(zé)某一方向上深度緩存排序計(jì)算的處理器稱為計(jì)算節(jié)點(diǎn)。如圖2所示,每一個(gè)計(jì)算節(jié)點(diǎn)上運(yùn)行著一個(gè)計(jì)算光分布的進(jìn)程,彼此相互獨(dú)立平行。利于讀寫鎖維護(hù)共用數(shù)據(jù)容器的正確使用。當(dāng)最后一個(gè)計(jì)算節(jié)點(diǎn)上光分布運(yùn)算進(jìn)程結(jié)束時(shí),共用數(shù)據(jù)容器中的光分布結(jié)果便是自然條件下多方向光同時(shí)入射植物冠層時(shí)產(chǎn)生的光輻射分布情況[10]。

    2 系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)

    植物冠層并行仿真系統(tǒng)采用客戶機(jī)?服務(wù)器系統(tǒng)結(jié)構(gòu)模式,并基于面向?qū)ο蟮乃枷脒M(jìn)行設(shè)計(jì),使用C++語(yǔ)言進(jìn)行實(shí)現(xiàn),并結(jié)合OpenGL做三維可視化顯示。整個(gè)系統(tǒng)采用模塊化的設(shè)計(jì)思想,將系統(tǒng)不同功能封裝在不同的邏輯模塊中。其中主要的模塊有:并行運(yùn)算模塊、服務(wù)器監(jiān)聽模塊、客戶機(jī)人機(jī)交互模塊和三維顯示模塊。

    2.1 并行計(jì)算模塊的設(shè)計(jì)和實(shí)現(xiàn)

    該模塊運(yùn)行在聯(lián)想深騰1800高性能多處理器計(jì)算機(jī)上,操作系統(tǒng)為Redhat5。算法采用MPI+C并行編程技術(shù),完成計(jì)算植物冠層光分布的深度緩存排序算法的并行化工作。其中主要類為:

    1) Light類:系統(tǒng)光源類。主要建模實(shí)現(xiàn)了單方向上的平行光線。

    2) Sky類:天空類。將天空半球的漫反射光線劃分為多方向上的平行入射光線,并為每一方向上的光線提供初始輻射強(qiáng)度。

    3) Triangle類:葉片面元類,記錄葉片面元坐標(biāo)信息,用于計(jì)算深度緩存排序。

    4) Calculatefork類:節(jié)點(diǎn)計(jì)算控制類。用于控制計(jì)算節(jié)點(diǎn)中光分布模擬的算法流程。

    5) Communicate類:節(jié)點(diǎn)通信控制類。負(fù)責(zé)不同計(jì)算節(jié)點(diǎn)之間的通信控制,管理讀寫鎖等工作。

    各類的具體關(guān)系見圖3。

    2.2 服務(wù)器監(jiān)聽模塊的設(shè)計(jì)和實(shí)現(xiàn)

    該模塊同樣運(yùn)行于高性能計(jì)算機(jī)上,主要負(fù)責(zé)監(jiān)聽客戶端發(fā)來(lái)的計(jì)算請(qǐng)求,根據(jù)請(qǐng)求建立并行運(yùn)算服務(wù),檢測(cè)用戶上傳數(shù)據(jù)是否合理,并在運(yùn)算過(guò)程中將實(shí)時(shí)信息傳給遠(yuǎn)程客戶端。模塊采用Socket套接字編程技術(shù)進(jìn)行網(wǎng)絡(luò)進(jìn)程間的通信。使用TCP/IP協(xié)議三步握手式建立連接方法,確保數(shù)據(jù)在本地服務(wù)器和遠(yuǎn)程客戶機(jī)之間的無(wú)誤傳輸。采用多進(jìn)程技術(shù)保證服務(wù)器可以同時(shí)響應(yīng)多個(gè)計(jì)算請(qǐng)求。主進(jìn)程一直等待監(jiān)聽,如有客戶端傳來(lái)的計(jì)算請(qǐng)求則開啟新的子進(jìn)程完成計(jì)算??偟膩?lái)說(shuō),該部分提供了并行計(jì)算模塊和遠(yuǎn)程客戶端之間的通信控制服務(wù)。其主要算法流圖如圖4所示。

    2.3 客戶機(jī)人機(jī)交互模塊的設(shè)計(jì)和實(shí)現(xiàn)

    該模塊運(yùn)行于普通個(gè)人電腦的Windows操作系統(tǒng)下,在VC++環(huán)境下采用MFC實(shí)現(xiàn)。它主要為用戶提供一個(gè)人機(jī)交互的GUI平臺(tái),接收用戶設(shè)計(jì)的運(yùn)算參數(shù)和數(shù)據(jù)文件并進(jìn)行上傳。向用戶顯示運(yùn)算過(guò)程信息,并在運(yùn)算結(jié)束后將服務(wù)器指定位置下的運(yùn)算結(jié)果下載到本地,供科學(xué)分析和三維顯示模塊生成仿真效果圖使用。從而達(dá)到用客戶端GUI控制服務(wù)器并行運(yùn)算的目的。

    2.4 三維顯示模塊的設(shè)計(jì)和實(shí)現(xiàn)endprint

    該模塊負(fù)責(zé)根據(jù)下載到客戶端的光分布運(yùn)算結(jié)果數(shù)據(jù)建立仿真效果圖并用OpenGL方法來(lái)實(shí)現(xiàn)三維渲染,完成系統(tǒng)所需的可視化處理和輸出。為分析人員提供直觀感性的仿真結(jié)果。植物冠層光分布的并行仿真計(jì)算系統(tǒng)有客戶端、服務(wù)器兩大部分,并由服務(wù)器監(jiān)聽模塊、人機(jī)交互模塊、并行計(jì)算模塊以及三維顯示模塊四個(gè)主要模塊組成。它們的相互關(guān)系如圖5所示。

    以上四部分相互協(xié)作完成了系統(tǒng)的主要數(shù)據(jù)輸入,并行仿真計(jì)算和結(jié)果輸出等功能。

    3 實(shí) 驗(yàn)

    植物冠層光分布的并行仿真計(jì)算系統(tǒng)服務(wù)器程序的兩部分在Linux系統(tǒng)下分別使用Inter編譯器和GUN編譯器編譯成功,客戶端程序在Windows上用Visual Studio 2005開發(fā)調(diào)試成功,圖6給出了仿真計(jì)算不同時(shí)期水稻光照分布后客戶端顯示的效果截圖。

    通過(guò)圖6可以清楚地看到水稻不同時(shí)期光照的效果圖,并且能夠了解服務(wù)器端運(yùn)行信息,使用戶感覺(jué)計(jì)算就像在本地運(yùn)行一樣。同時(shí)由于每個(gè)計(jì)算節(jié)點(diǎn)任務(wù)量相對(duì)固定,均為單一方向光分布計(jì)算。所以該系統(tǒng)仿真計(jì)算耗時(shí)不隨天空光照方向數(shù)增多而增加。服務(wù)器返回結(jié)果顯示,每次計(jì)算耗時(shí)僅為十多秒,這極大地提高了計(jì)算效率,是普通串行算法無(wú)法達(dá)到的。

    另外系統(tǒng)客戶端不但能生成冠層的三維結(jié)構(gòu),還能應(yīng)用圖像控制標(biāo)簽欄對(duì)仿真效果圖進(jìn)行縮放、3D旋轉(zhuǎn)、顏色渲染等操作,這方便了對(duì)三維冠層的觀察操作,提供了一種良好的交互方式。同時(shí)系統(tǒng)還能通過(guò)任務(wù)標(biāo)簽欄管理計(jì)算結(jié)果文件。使科研人員可以進(jìn)一步查閱數(shù)據(jù)分析仿真結(jié)果。

    4 結(jié) 論

    綜上,植物冠層光分布的并行仿真計(jì)算系統(tǒng)將傳統(tǒng)的深度緩存排序算法進(jìn)行了并行化處理,并基于面向?qū)ο蟮脑O(shè)計(jì)思想進(jìn)行實(shí)現(xiàn),系統(tǒng)主要部分實(shí)現(xiàn)模塊化管理,使它擁有良好的擴(kuò)展性和維護(hù)性。目前系統(tǒng)基本功能已經(jīng)實(shí)現(xiàn),接下來(lái)還可以在客戶端根據(jù)用戶多種分析需要進(jìn)一步添加更多的計(jì)算結(jié)果數(shù)據(jù)顯示分析方法,使其不斷完善。而在服務(wù)器端,目前系統(tǒng)實(shí)現(xiàn)了深度緩存排序算法的并行化,在其并行化的基礎(chǔ)上還可以將其他光照分布算法做同樣處理并移植到大型機(jī)去進(jìn)行仿真計(jì)算,隨著系統(tǒng)的不斷成熟和發(fā)展,系統(tǒng)客戶端可以廣泛提供給沒(méi)有條件購(gòu)買大型機(jī)的科研教學(xué)用戶使用,使更多人可以方便快速地進(jìn)行仿真計(jì)算。相信在不遠(yuǎn)的將來(lái),該系統(tǒng)一定會(huì)成為科研和教學(xué)中十分有用的冠層光輻射分布計(jì)算的仿真模擬平臺(tái)。

    參考文獻(xiàn)

    [1] VERHOEF W. Light scattering by leaf layers with application to canopy reflectance modeling: the SAIL model [J]. Remote sensing of environment, 1984, 16(2): 125?141.

    [2] GORAL C M, TORRANCE K E, GREENBERG D P. Modeling the interaction of light between diffuse surfaces [J]. ACM SIGGRAPH computer graphics, 2008, 18(3): 213?222.

    [3] CHELLE M, ANDRIEU B. The nested radiosity model for the distribution of light within plant canopies [J]. Ecological modelling, 1998, 111(1): 75?91.

    [4] 谷傳才,關(guān)楠,于金銘,等.多處理器混合關(guān)鍵性系統(tǒng)中的劃分調(diào)度策略[J].軟件學(xué)報(bào),2014,25(2):284?297.

    GU Chuancai, GUAN Nan, YU Jinming, et al. The division scheduling strategy of multiprocessor hybrid critical system [J]. Journal of software, 2014, 25(2): 284?297.

    [5] 盧云娥,黃宗宇,李超陽(yáng),等.基于微機(jī)集群系統(tǒng)的MPI并行計(jì)算[J].電子設(shè)計(jì)工程,2011,19(5):78?81.

    LU Yune, HUANG Zongyu, LI Chaoyang, et al. Parallel computation of MPI based on microcomputer cluster system [J]. Electronic design engineering, 2011, 19(5): 78?81.

    [6] 郭惠,符紅光,羅東輝.三維動(dòng)態(tài)幾何中直線消隱的實(shí)現(xiàn)[J].計(jì)算機(jī)應(yīng)用,2007,27(3):663?665.

    GUO Hui, FU Hongguang, LUO Donghui. The realization of linear cancellation in three?dimensional dynamic geometry [J]. Computer application, 2007, 27(3): 663?665.

    [7] 肖強(qiáng),勞彩蓮,王春霞.植物冠層光輻射分布的快速深度緩存算法[J].計(jì)算機(jī)工程,2010,36(15):288?290.

    XIAO Qiang, LAO Cailian, WANG Chunxia. The rapid depth cache algorithm of the photoemission distribution of plant canopy [J]. Computer engineering, 2010, 36(15): 288?290.

    [8] 李靜,王文成,吳恩華.基于空盒自適應(yīng)生成的動(dòng)態(tài)場(chǎng)景光線跟蹤計(jì)算[J].計(jì)算機(jī)學(xué)報(bào),2009,32(6):1172?1182.

    LI Jing, WANG Wencheng, WU Enhua. The dynamic scene light tracking calculation based on the adaptive generation of the empty box [J]. Chinese journal of computers, 2009, 32(6): 1172?1182.

    [9] 黃帆,勞彩蓮,肖翠霞.基于光線跟蹤的冠層光分布模型參數(shù)研究[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2013,18(6):96?101.

    HUANG Fan, LAO Cailian, XIAO Cuixia. Parameters in the model of canopy light distribution with ray tracing approach [J]. Journal of China Agricultural University, 2013, 18(6): 96?101.

    [10] 王昊鵬,趙凱,宋鳳斌,等.植物冠層光照的建模與虛擬現(xiàn)實(shí)仿真研究[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究,2010,26(1):68?74.

    WANG Haopeng, ZHAO Kai, SONG Fengbin, et al. Study on the modeling of canopy illumination and virtual reality simulation [J]. System sciences and comprehensive studies in agriculture, 2010, 26(1): 68?74.endprint

    猜你喜歡
    客戶機(jī)
    Reducing the global cancer burden with gastrointestinal screening: China’s 30 years practice
    虛擬技術(shù)在計(jì)算機(jī)實(shí)驗(yàn)教學(xué)示范中心的應(yīng)用研究
    卷宗(2018年7期)2018-06-06 05:27:16
    隔山亦能打牛,本本巧變遠(yuǎn)控利器
    基于雙TLB的二進(jìn)制翻譯訪存性能優(yōu)化
    用移動(dòng)設(shè)備完成環(huán)境監(jiān)測(cè)
    售后服務(wù)管理系統(tǒng)的設(shè)計(jì)
    瘦客戶機(jī):安全與便捷的選擇
    惠普核心代理:瘦客戶機(jī)正在從邊緣走向主流
    升騰瘦客戶機(jī)借神碼翱翔“云端”
    基于Web數(shù)據(jù)提高訪問(wèn)速度的方法
    亚洲精品美女久久av网站| 天天躁夜夜躁狠狠躁躁| 美女中出高潮动态图| 亚洲国产日韩一区二区| 日韩伦理黄色片| 久久久久久久久免费视频了| 成人二区视频| 老汉色av国产亚洲站长工具| 香蕉国产在线看| 黑人猛操日本美女一级片| 99九九在线精品视频| 精品国产超薄肉色丝袜足j| 久久毛片免费看一区二区三区| 国产免费福利视频在线观看| 一本—道久久a久久精品蜜桃钙片| 久久 成人 亚洲| 秋霞在线观看毛片| 国产毛片在线视频| 只有这里有精品99| 国产日韩一区二区三区精品不卡| 精品少妇一区二区三区视频日本电影 | 日韩一本色道免费dvd| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频精品| 亚洲国产精品国产精品| 国产有黄有色有爽视频| 女人精品久久久久毛片| 丰满少妇做爰视频| 日本欧美视频一区| 电影成人av| 亚洲av福利一区| 高清欧美精品videossex| 97在线人人人人妻| 久久精品国产亚洲av高清一级| 一级毛片 在线播放| 好男人视频免费观看在线| 久久久久久久国产电影| 久久国产精品大桥未久av| 国产国语露脸激情在线看| 久久人人爽av亚洲精品天堂| 午夜老司机福利剧场| 亚洲四区av| 在线观看www视频免费| 热re99久久精品国产66热6| 亚洲精华国产精华液的使用体验| 视频在线观看一区二区三区| 天美传媒精品一区二区| 99久久人妻综合| 成人免费观看视频高清| 亚洲成国产人片在线观看| 一区二区三区乱码不卡18| 老熟女久久久| 肉色欧美久久久久久久蜜桃| 成人国产av品久久久| 国产免费一区二区三区四区乱码| 国产野战对白在线观看| 亚洲av欧美aⅴ国产| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲综合精品二区| tube8黄色片| 国产av码专区亚洲av| 日本vs欧美在线观看视频| 亚洲婷婷狠狠爱综合网| 午夜影院在线不卡| 欧美精品一区二区大全| 丰满迷人的少妇在线观看| 青春草国产在线视频| 国产精品国产av在线观看| 久久免费观看电影| 两个人看的免费小视频| 精品99又大又爽又粗少妇毛片| 国产成人精品福利久久| 亚洲av成人精品一二三区| 91精品伊人久久大香线蕉| 久久久久久人妻| 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 久久精品夜色国产| 九草在线视频观看| 国产精品亚洲av一区麻豆 | 精品亚洲成国产av| 最近最新中文字幕大全免费视频 | 国产熟女午夜一区二区三区| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 女人久久www免费人成看片| 26uuu在线亚洲综合色| 日日爽夜夜爽网站| 亚洲精品一区蜜桃| 日本-黄色视频高清免费观看| 又粗又硬又长又爽又黄的视频| 久久久久网色| 国产成人一区二区在线| 国产一区二区在线观看av| 午夜老司机福利剧场| 国产成人精品一,二区| 亚洲一码二码三码区别大吗| 9色porny在线观看| 日韩制服丝袜自拍偷拍| 纵有疾风起免费观看全集完整版| 人人妻人人添人人爽欧美一区卜| av在线app专区| 亚洲欧美精品自产自拍| 亚洲国产精品一区三区| 纯流量卡能插随身wifi吗| 成人亚洲精品一区在线观看| 2018国产大陆天天弄谢| 在线观看www视频免费| 岛国毛片在线播放| 中文精品一卡2卡3卡4更新| 丝袜脚勾引网站| 蜜桃在线观看..| 国产精品无大码| 亚洲,一卡二卡三卡| 老熟女久久久| 亚洲av日韩在线播放| 国产精品 欧美亚洲| 国产精品一二三区在线看| 久久久久国产网址| 十分钟在线观看高清视频www| 久久久久网色| 美女主播在线视频| 我要看黄色一级片免费的| 精品国产一区二区三区四区第35| 丝袜人妻中文字幕| 黄色怎么调成土黄色| 女的被弄到高潮叫床怎么办| 久久99一区二区三区| av一本久久久久| 老汉色∧v一级毛片| 午夜福利在线免费观看网站| 免费久久久久久久精品成人欧美视频| 晚上一个人看的免费电影| 日韩熟女老妇一区二区性免费视频| 国产精品无大码| 欧美国产精品一级二级三级| 99久国产av精品国产电影| 午夜福利一区二区在线看| 大香蕉久久成人网| 99久久综合免费| 久久久国产欧美日韩av| 男的添女的下面高潮视频| 久久热在线av| 亚洲av在线观看美女高潮| 高清黄色对白视频在线免费看| 亚洲av电影在线观看一区二区三区| 人妻系列 视频| 老汉色av国产亚洲站长工具| 日韩av免费高清视频| 国产精品国产三级国产专区5o| 国产成人精品福利久久| 多毛熟女@视频| 亚洲av综合色区一区| 夫妻性生交免费视频一级片| 国产人伦9x9x在线观看 | 一本色道久久久久久精品综合| 美女国产高潮福利片在线看| www.av在线官网国产| 一级毛片 在线播放| 国产精品免费视频内射| 免费高清在线观看日韩| 欧美精品av麻豆av| 老女人水多毛片| 美女大奶头黄色视频| 亚洲国产精品一区三区| 极品少妇高潮喷水抽搐| 国产精品偷伦视频观看了| 午夜福利视频精品| 另类亚洲欧美激情| 精品一区二区三区四区五区乱码 | 中文字幕色久视频| 麻豆精品久久久久久蜜桃| 成年人免费黄色播放视频| 亚洲伊人久久精品综合| 美女国产高潮福利片在线看| 女性生殖器流出的白浆| 中文精品一卡2卡3卡4更新| 麻豆乱淫一区二区| 狠狠婷婷综合久久久久久88av| 色哟哟·www| 日韩一区二区三区影片| 黄频高清免费视频| 免费观看性生交大片5| 国产亚洲精品第一综合不卡| 91aial.com中文字幕在线观看| 精品久久久精品久久久| 国产又色又爽无遮挡免| 色吧在线观看| 国产精品蜜桃在线观看| 久久免费观看电影| 中文字幕精品免费在线观看视频| 久久av网站| 极品人妻少妇av视频| 曰老女人黄片| 夫妻性生交免费视频一级片| av天堂久久9| 日韩制服丝袜自拍偷拍| 男女边摸边吃奶| av天堂久久9| av线在线观看网站| 久久精品亚洲av国产电影网| www日本在线高清视频| 久久久欧美国产精品| 免费播放大片免费观看视频在线观看| 亚洲欧美一区二区三区国产| 99国产综合亚洲精品| 一边亲一边摸免费视频| 久久人人爽人人片av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜激情av网站| 国产一区二区三区综合在线观看| 午夜激情久久久久久久| 丝袜脚勾引网站| 亚洲精华国产精华液的使用体验| 边亲边吃奶的免费视频| 国产亚洲午夜精品一区二区久久| av国产久精品久网站免费入址| 在线观看www视频免费| 久久av网站| 国产精品一国产av| 看免费av毛片| 中文字幕制服av| av有码第一页| 日韩伦理黄色片| 毛片一级片免费看久久久久| 精品99又大又爽又粗少妇毛片| av有码第一页| 日韩一卡2卡3卡4卡2021年| 久久97久久精品| 久久精品国产a三级三级三级| 成人亚洲欧美一区二区av| 熟女少妇亚洲综合色aaa.| 久久久久久久大尺度免费视频| 最近最新中文字幕大全免费视频 | 搡老乐熟女国产| 国产野战对白在线观看| 亚洲欧美成人综合另类久久久| 国产精品免费视频内射| 满18在线观看网站| 夜夜骑夜夜射夜夜干| 亚洲在久久综合| 高清av免费在线| 亚洲精品国产av蜜桃| 久久久久久久久久久久大奶| 日韩免费高清中文字幕av| 国产成人一区二区在线| 伦理电影免费视频| 国产淫语在线视频| 免费不卡的大黄色大毛片视频在线观看| 在线观看www视频免费| 国产福利在线免费观看视频| av电影中文网址| 在线精品无人区一区二区三| 久久影院123| 天天影视国产精品| 亚洲伊人久久精品综合| 久久久久久人妻| 只有这里有精品99| 侵犯人妻中文字幕一二三四区| 国产熟女欧美一区二区| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 日韩中文字幕视频在线看片| 老熟女久久久| 日韩精品有码人妻一区| 欧美成人精品欧美一级黄| 中文天堂在线官网| 久久久久网色| 国产免费福利视频在线观看| 春色校园在线视频观看| 久久国产精品大桥未久av| 天美传媒精品一区二区| 一本久久精品| 久久国产精品大桥未久av| 人人妻人人添人人爽欧美一区卜| 尾随美女入室| 中文字幕最新亚洲高清| 考比视频在线观看| 国产精品 欧美亚洲| tube8黄色片| 搡女人真爽免费视频火全软件| 精品国产超薄肉色丝袜足j| 中文天堂在线官网| av女优亚洲男人天堂| 国产精品av久久久久免费| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 欧美精品国产亚洲| 国产日韩欧美亚洲二区| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦啦在线视频资源| 久久99精品国语久久久| 三级国产精品片| 国产福利在线免费观看视频| 色婷婷久久久亚洲欧美| 在线看a的网站| 亚洲欧美一区二区三区久久| 亚洲成av片中文字幕在线观看 | 人人妻人人添人人爽欧美一区卜| 91aial.com中文字幕在线观看| 美国免费a级毛片| 日韩欧美一区视频在线观看| 亚洲欧洲国产日韩| 亚洲av免费高清在线观看| 欧美老熟妇乱子伦牲交| 日本wwww免费看| 一级黄片播放器| 亚洲成av片中文字幕在线观看 | 亚洲,一卡二卡三卡| 十八禁高潮呻吟视频| 精品国产乱码久久久久久男人| 久久影院123| 考比视频在线观看| 国产片内射在线| 久久狼人影院| 有码 亚洲区| 波多野结衣av一区二区av| 国产极品粉嫩免费观看在线| 天天躁夜夜躁狠狠躁躁| 亚洲欧美成人综合另类久久久| 久热久热在线精品观看| 69精品国产乱码久久久| 免费看av在线观看网站| 亚洲熟女精品中文字幕| 婷婷成人精品国产| 五月伊人婷婷丁香| 中文字幕制服av| 亚洲天堂av无毛| 国产日韩欧美视频二区| 日产精品乱码卡一卡2卡三| 国产精品99久久99久久久不卡 | 有码 亚洲区| 日韩视频在线欧美| 有码 亚洲区| 欧美97在线视频| 中文字幕色久视频| tube8黄色片| 99久久中文字幕三级久久日本| 乱人伦中国视频| 人人妻人人澡人人看| 晚上一个人看的免费电影| 午夜av观看不卡| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 在现免费观看毛片| 久久久久久久久久久久大奶| 免费黄色在线免费观看| 18在线观看网站| 午夜激情av网站| 美女大奶头黄色视频| 9热在线视频观看99| 午夜福利网站1000一区二区三区| 国产精品久久久久久精品电影小说| 中文字幕人妻丝袜制服| 国产老妇伦熟女老妇高清| 久久久久久久久久人人人人人人| 国产精品 欧美亚洲| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 久久热在线av| 日韩制服丝袜自拍偷拍| 美女大奶头黄色视频| 日日爽夜夜爽网站| 亚洲精品久久成人aⅴ小说| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 国产亚洲最大av| 男人舔女人的私密视频| 母亲3免费完整高清在线观看 | 久久精品国产亚洲av天美| 国产黄色免费在线视频| 性色avwww在线观看| 国产毛片在线视频| 午夜激情久久久久久久| 99久久人妻综合| 欧美精品一区二区免费开放| 超碰成人久久| 老汉色av国产亚洲站长工具| 久久国产亚洲av麻豆专区| 美女高潮到喷水免费观看| 国产免费视频播放在线视频| 叶爱在线成人免费视频播放| 免费在线观看视频国产中文字幕亚洲 | 9色porny在线观看| 日韩,欧美,国产一区二区三区| 啦啦啦在线免费观看视频4| 岛国毛片在线播放| 午夜老司机福利剧场| 国产视频首页在线观看| 久久精品aⅴ一区二区三区四区 | 欧美97在线视频| 久久久国产一区二区| 午夜免费观看性视频| 免费av中文字幕在线| 亚洲精华国产精华液的使用体验| 精品少妇久久久久久888优播| www.av在线官网国产| 国产成人a∨麻豆精品| 国产男女超爽视频在线观看| 不卡av一区二区三区| 天天操日日干夜夜撸| 一级a爱视频在线免费观看| 国产精品国产三级专区第一集| 在线观看www视频免费| 成人黄色视频免费在线看| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 久久婷婷青草| 黄片播放在线免费| 90打野战视频偷拍视频| 春色校园在线视频观看| 精品人妻偷拍中文字幕| 一区二区三区激情视频| 欧美日韩精品成人综合77777| 亚洲图色成人| 大码成人一级视频| 国产av一区二区精品久久| 成人18禁高潮啪啪吃奶动态图| 国产成人一区二区在线| 热re99久久国产66热| 中文字幕另类日韩欧美亚洲嫩草| 777久久人妻少妇嫩草av网站| 久久免费观看电影| 亚洲精品中文字幕在线视频| 极品人妻少妇av视频| 国产精品久久久av美女十八| 精品一区在线观看国产| 国产一区二区三区综合在线观看| 亚洲国产欧美网| 久久女婷五月综合色啪小说| 欧美国产精品va在线观看不卡| 精品福利永久在线观看| 青春草国产在线视频| 少妇被粗大猛烈的视频| 新久久久久国产一级毛片| 欧美日韩一级在线毛片| 亚洲国产色片| 99热全是精品| 亚洲一码二码三码区别大吗| 久久国产亚洲av麻豆专区| 午夜免费观看性视频| 爱豆传媒免费全集在线观看| 熟妇人妻不卡中文字幕| 91国产中文字幕| 人妻少妇偷人精品九色| 亚洲av欧美aⅴ国产| 在线观看人妻少妇| 国产一区二区在线观看av| 午夜精品国产一区二区电影| 青草久久国产| 一级黄片播放器| 欧美精品人与动牲交sv欧美| 亚洲四区av| 熟女电影av网| 国产一区二区 视频在线| 日本欧美视频一区| 99热网站在线观看| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 欧美日韩国产mv在线观看视频| 蜜桃国产av成人99| 欧美97在线视频| 亚洲欧美成人综合另类久久久| 丝袜在线中文字幕| a级毛片在线看网站| 一区二区三区乱码不卡18| 欧美激情 高清一区二区三区| 女人精品久久久久毛片| 香蕉精品网在线| 91午夜精品亚洲一区二区三区| 一本大道久久a久久精品| 久久久久久久久久人人人人人人| 在线观看免费高清a一片| 国产欧美亚洲国产| 国产1区2区3区精品| 看免费成人av毛片| 日韩av免费高清视频| 九九爱精品视频在线观看| 一区在线观看完整版| 性色avwww在线观看| 人体艺术视频欧美日本| 最近手机中文字幕大全| 大陆偷拍与自拍| 嫩草影院入口| 精品国产一区二区三区久久久樱花| 在线观看免费高清a一片| 18+在线观看网站| 国产精品一区二区在线不卡| 国产成人一区二区在线| 青春草亚洲视频在线观看| 久久这里有精品视频免费| 中国三级夫妇交换| 十八禁高潮呻吟视频| 久久久亚洲精品成人影院| 在线免费观看不下载黄p国产| 两性夫妻黄色片| 一区二区三区激情视频| 国产精品成人在线| 欧美成人精品欧美一级黄| 成人国产麻豆网| 国产野战对白在线观看| 亚洲伊人久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 丝袜美足系列| 欧美激情高清一区二区三区 | 又粗又硬又长又爽又黄的视频| 成人午夜精彩视频在线观看| 亚洲精品国产色婷婷电影| 热99久久久久精品小说推荐| 亚洲一级一片aⅴ在线观看| 极品少妇高潮喷水抽搐| 色吧在线观看| 青春草视频在线免费观看| 欧美另类一区| 国产成人一区二区在线| 天堂8中文在线网| 欧美人与善性xxx| 人人妻人人添人人爽欧美一区卜| 日韩熟女老妇一区二区性免费视频| 国产在视频线精品| 中文天堂在线官网| 人妻系列 视频| 看非洲黑人一级黄片| 日本91视频免费播放| 高清av免费在线| 1024香蕉在线观看| 免费av中文字幕在线| 99re6热这里在线精品视频| 嫩草影院入口| 丝袜美腿诱惑在线| 国产精品熟女久久久久浪| 亚洲精品久久久久久婷婷小说| 欧美日韩亚洲高清精品| 免费在线观看黄色视频的| 老女人水多毛片| 777米奇影视久久| kizo精华| 啦啦啦在线免费观看视频4| 精品国产乱码久久久久久男人| a级毛片黄视频| 欧美亚洲 丝袜 人妻 在线| 天天躁夜夜躁狠狠躁躁| av免费在线看不卡| 天天影视国产精品| 亚洲欧洲国产日韩| 国产免费又黄又爽又色| 国产高清国产精品国产三级| 欧美黄色片欧美黄色片| av卡一久久| 制服丝袜香蕉在线| 青春草国产在线视频| 视频区图区小说| 欧美97在线视频| 人人澡人人妻人| 色网站视频免费| 黄色视频在线播放观看不卡| 亚洲熟女精品中文字幕| 精品亚洲成国产av| 男女高潮啪啪啪动态图| 捣出白浆h1v1| 中文字幕精品免费在线观看视频| freevideosex欧美| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 男人添女人高潮全过程视频| 人人澡人人妻人| 日韩熟女老妇一区二区性免费视频| 最近中文字幕高清免费大全6| 99国产精品免费福利视频| 中文字幕av电影在线播放| 日本免费在线观看一区| 欧美日韩国产mv在线观看视频| 欧美黄色片欧美黄色片| 青青草视频在线视频观看| 国产成人精品福利久久| 一级毛片我不卡| 免费观看在线日韩| 亚洲国产色片| 男人操女人黄网站| 国产精品麻豆人妻色哟哟久久| 91久久精品国产一区二区三区| videosex国产| 欧美日韩精品网址| 在线观看免费日韩欧美大片| 亚洲成人一二三区av| 日韩不卡一区二区三区视频在线| 老熟女久久久| 国产精品人妻久久久影院| 一区二区日韩欧美中文字幕| 久久人人爽av亚洲精品天堂| 久久久国产精品麻豆| 高清欧美精品videossex| 91午夜精品亚洲一区二区三区| 超碰97精品在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩视频在线欧美| 两个人看的免费小视频| 亚洲成色77777| a级毛片黄视频| 亚洲美女黄色视频免费看| 国产在线一区二区三区精| 久久久国产一区二区| 亚洲成av片中文字幕在线观看 | 新久久久久国产一级毛片| 国产精品三级大全| 亚洲成av片中文字幕在线观看 | 一级爰片在线观看| 男女高潮啪啪啪动态图| 成人亚洲精品一区在线观看| 精品国产超薄肉色丝袜足j| 亚洲欧洲日产国产|