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

    基于特征點篩選的漁船軌跡壓縮方法及優(yōu)化

    2023-08-22 11:06:12馬蘇上王澤宇劉明劍趙文彬
    漁業(yè)現(xiàn)代化 2023年4期
    關(guān)鍵詞:壓縮率漁船軌跡

    馬蘇上,王澤宇,劉明劍,2,李 洋,趙文彬

    (1 大連海洋大學(xué)信息工程學(xué)院,遼寧 大連 116023;2 設(shè)施漁業(yè)教育部重點試驗室(大連海洋大學(xué)),遼寧 大連 116023)

    近年來漁民捕撈經(jīng)濟在海洋經(jīng)濟中日益重要[1],2021年海洋捕撈產(chǎn)值為2 303.72億元,漁船作為漁民捕撈的工具,加強對漁船的行為監(jiān)管,是保障漁業(yè)健康有序地發(fā)展的重要手段[2-3]。漁船在作業(yè)過程中呈現(xiàn)出撈區(qū)域繁多、分布廣、持續(xù)作業(yè)時間長等特點[4],這就需要在監(jiān)管過程中,對漁船不同作業(yè)特點的軌跡進行差異化分析。目前常用的船舶軌跡跟蹤方式是使用船舶自動識別系統(tǒng)[5](Automatic Identification System,AIS)對漁船行駛軌跡進行記錄與存儲,但在此過程中,AIS系統(tǒng)會產(chǎn)生大量的軌跡數(shù)據(jù)[6],不僅消耗系統(tǒng)存儲與計算資源,還存在大量的冗雜和錯誤信息[7],不利于對船舶軌跡的有效分析。如海上漁船由于捕撈作業(yè)會產(chǎn)生大量的靜態(tài)信息,此時僅需保留幾條有效的信息,就可以有效的反映漁船的作業(yè)狀態(tài),因此需要采用軌跡壓縮方法對漁船軌跡數(shù)據(jù)進行篩選,從而為漁船行為有效分析提供高精度的軌跡數(shù)據(jù)支撐。

    針對上述需求,國內(nèi)外的研究學(xué)者開展了關(guān)于漁船軌跡壓縮方面的相關(guān)研究,目前船舶軌跡壓縮方法大致分為兩類[8]:在線軌跡壓縮方法;離線軌跡壓縮方法。在線壓縮算法為局部處理算法,處理速度較快。其中,Keoghe等[9]采用滑動窗口方式實現(xiàn)了在線壓縮;文獻[10-12]通過在滑動窗口中加入閾值方式來提升壓縮效率。但以上方法沒有考慮部分特征點對整條軌跡曲線存在軌跡信息量不均衡現(xiàn)象,當(dāng)刪除一些極值點后易造成軌跡失真,不利于漁船軌跡分析;離線軌跡壓縮方法,則是對于全局?jǐn)?shù)據(jù)進行分析處理,保證了所有軌跡點都需進行判別[13],更適于對漁船軌跡進行分析。Mancini等[14]提出了均勻采樣方法,實現(xiàn)軌跡壓縮;Douglas等[15]提出的Douglas-Peucker(簡稱DP)算法,將軌跡的起始點相連接作為基線,根據(jù)其余軌跡點到直線的垂直歐式距離來進行軌跡的壓縮。在DP算法的基礎(chǔ)上,文獻[16-17]通過修改數(shù)據(jù)集、去除數(shù)據(jù)冗余信息來實現(xiàn)高效漁船軌跡壓縮;文獻[18-19]考慮了船舶軌跡的壓縮形狀,有效保留了船舶航行中作業(yè)的重要特征信息,減少了算法的時間復(fù)雜度。文獻[20-21]著重研究閾值對于軌跡壓縮的影響,根據(jù)不同情況改變閾值,實現(xiàn)更精確的壓縮。文獻[22-24]通過引入時間維度信息來對漁船航向速度等動態(tài)數(shù)據(jù)進行壓縮,能夠保留更多的重要動態(tài)行為軌跡點。

    盡管上述方法對于離線軌跡壓縮的處理趨于完善,但目前離線壓縮算法無法妥善處理和反映船舶的這些動態(tài)變化特征信息,導(dǎo)致其算法運行的時間復(fù)雜度較高,且軌跡壓縮的誤差較大,致使軌跡壓縮后存在失真的問題。針對上述問題,本研究提出了FMDP漁船軌跡壓縮方法,在保證壓縮率的前提下最大程度保留了漁船軌跡的特征信息。

    1 算法建立和實現(xiàn)

    1.1 FMDP軌跡壓縮方法

    設(shè)漁船在行駛過程中,AIS軌跡點為p(x,y),由軌跡點p所組成的軌跡點集為T。漁船在某一位置長時間停留或移動距離極短的點被稱之為停留點,停留點集為ts。

    FMDP軌跡壓縮方法:首先,將原始軌跡T進行停留點篩選出并計入停留點集ts(ts?T),同時,針對漁船行為停留時間較長,造成的差異化停留與轉(zhuǎn)向狀態(tài)特征,設(shè)計壓縮軌跡權(quán)重機制,對停留點集ts再次篩選;其次,使用加入角度閾值的DP算法進行軌跡壓縮處理完成特征點篩選,組成特征點集tr(tr?T);最后,取ts和tr的交集為整合點集tA,將tA按照時間順序的軌跡點串聯(lián)得到新的壓縮軌跡。

    1.2 留點篩選

    1.2.1 軌跡點間距離

    停留是移動對象的基本運動狀態(tài),代表著軌跡中具有重要特征的點序列[25]。當(dāng)漁船停留時可能會因為外界因素或者自身因素發(fā)生偏移現(xiàn)象,導(dǎo)致停留點并不是固定的點。為了可以更準(zhǔn)確判別軌跡點是否屬于停留點,進行如下設(shè)定:將某條漁船的原始漁船軌跡點集T依次分成n個船速計算區(qū)域,其中每個區(qū)域包含五百個AIS軌跡數(shù)據(jù)點。考慮任意相鄰的軌跡點之間距離并不屬于直線距離,因此本研究將軌跡點的經(jīng)緯度坐標(biāo)轉(zhuǎn)換為笛卡爾坐標(biāo),進行距離度量。軌跡點距離d是兩個相鄰的軌跡點的恒向線[26],其距離公式如下:

    (1)

    式中:Δx和Δy分別為兩個相鄰軌跡點的經(jīng)緯度差值[27],(°);c為結(jié)合后續(xù)計算所設(shè)計的參數(shù)系數(shù);r為地球的半徑,km;a是用來將經(jīng)度轉(zhuǎn)換的等量維度系數(shù)。

    1.2.2 停留點閾值設(shè)定

    依據(jù)不同區(qū)域行駛距離與所用的時間,求取出不同區(qū)域的平均速度作為該區(qū)域漁船行駛速度,并篩選出最慢的行駛速度vmin,并根據(jù)AIS數(shù)據(jù)發(fā)送的最短時間間隔與速度vmin,求出最短的距離l0,作為停留點閾值:

    l0=vmin×tmin

    (2)

    式中:l0為最短的距離,km;vmin為行駛速度,km/s;tmin為停留時間,s。

    1.2.3 停留點篩選

    隨后將點集T中的軌跡點pi到pi+1的距離為di,pi+1到pi+2所組成線段距離為di+1,兩線段的距離之和為di+di+1,其距離之和根據(jù)公式(1)將經(jīng)緯度轉(zhuǎn)換為笛卡爾坐標(biāo)系再進行計算。將距離之和與閾值l0進行對比,如果di小于l0則取中間的點pi+1作為停留點,之后繼續(xù)計算軌跡點pi+1到pi+2,pi+2到pi+3的距離之和與閾值l0進行判定,直至將軌跡點集T中所有符合條件的停留點全部篩選出為止,并將這些停留點組成停留點點集ts。

    1.3 特征點篩選

    1.3.1 原始DP算法的軌跡特征點篩選處理

    (3)

    式中:x與y為通過投影將經(jīng)緯度轉(zhuǎn)化到同一平面坐標(biāo)的經(jīng)度與緯度。k和b是笛卡爾坐標(biāo)系中的常數(shù)。

    1.3.2 加入了角度閾值的軌跡壓縮特征點篩選

    傳統(tǒng)的DP軌跡壓縮方法中,僅考慮了距離特性而忽略了船舶的轉(zhuǎn)向變化特性,無法保證將全部特征點進行篩選,導(dǎo)致壓縮后的軌跡不能完整反應(yīng)出原始軌跡特征。因此為保證船舶的轉(zhuǎn)向信息最大保留,在DP算法原有距離閾值的基礎(chǔ)上加入了角度閾值篩選特征點方法。

    軌跡點的轉(zhuǎn)向角度由連續(xù)的3個AIS點構(gòu)成,3個AIS點依次相連,所構(gòu)成線段的角度即為轉(zhuǎn)向角度。設(shè)定軌跡點的角度閾值,當(dāng)軌跡點滿足角度閾值條件,但未滿足原始DP算法的特征點選擇條件時,仍可作為特征點計入特征點集TR。

    角度閾值的選取與軌跡壓縮率有著密不可分的聯(lián)系。為了避免選取不合適角度閾值后導(dǎo)致算法對數(shù)據(jù)的壓縮率效率降低,將采取真實漁船軌跡數(shù)據(jù)進行測算來獲得每條船在壓縮時的最佳角度閾值,通過對比不同角度下的角度軌跡壓縮率,判斷出最佳的角度閾值,計算不同角度的范圍壓縮點的剩余點數(shù)目Va與原始軌跡點數(shù)O結(jié)合求出角度壓縮率Ca,其公式如下:

    (4)

    以2018年10月1日至7日大連世想海洋科技有限公司提供的某條漁船軌跡數(shù)據(jù)為例進行角度閾值選取。由于該船在此時間段AIS數(shù)據(jù)龐大,微小的角度變化就會造成壓縮率巨大變化,因此該算法分兩階段篩選角度閾值。第一階段篩選大致范圍,規(guī)定每10°一個區(qū)間,根據(jù)峰值挑選出合適的角度范圍,如圖1所示,當(dāng)這條漁船轉(zhuǎn)向角在40°至60°以及130°至150°時存在有最大的角度壓縮CA,因此需要對于以上兩個區(qū)間范圍繼續(xù)通過第二階段進行篩選。

    圖1 第一階段軌跡壓縮率

    第二階段將挑選出的角度范圍再次細分進行二次篩選,以確保壓縮率最大化,規(guī)定每2°一個區(qū)間,同樣根據(jù)峰值篩選出最佳角度范圍,完成最終的角度閾值確立,如圖2所示。

    注:(a)圖為40°至60°區(qū)間內(nèi)的軌跡壓縮率,(b)圖為130°至150°區(qū)間內(nèi)的軌跡壓縮率

    1.4 軌跡點集結(jié)合處理

    1.5 復(fù)雜軌跡壓縮權(quán)重機制

    考慮到漁船的作業(yè)類型不同,部分捕撈作業(yè)花費時間較長,如捕撈金槍魚和石斑魚等魚類,需要在同一地點持續(xù)2~3 h的捕撈作業(yè)。僅依靠上述軌跡壓縮方法進行壓縮時,由于大量停留點被壓縮導(dǎo)致軌跡信息嚴(yán)重缺失。針對上述情況,本研究設(shè)計了壓縮軌跡權(quán)重機制,其目的在于通過對停留點進行再次篩選,避免因為坐標(biāo)信息重復(fù)或相近而導(dǎo)致停留點被大量刪除。從而保留重要的軌跡特征信息。其具體方法如下:加入權(quán)重的軌跡壓縮,是在進行FMDP算法的停留點篩選后,將出現(xiàn)的多個停留點連續(xù)區(qū)間視為停留點區(qū)間。在每個停留點區(qū)間,設(shè)停留點數(shù)為ns,停留時間為ts,停留閾值為β=ns/ts,對β進行向下取整為δ,那么δ為該停留點區(qū)間應(yīng)保留的停留點個數(shù)。

    在每個停留點區(qū)間,由于漁船長時間的停留出現(xiàn)多個相同的坐標(biāo)點集,每個點集均是相同的坐標(biāo)點所組成。設(shè)在此區(qū)間中坐標(biāo)點集的個數(shù)為m,將坐標(biāo)點集個數(shù)m與停留點個數(shù)δ相比較:

    (1)當(dāng)坐標(biāo)點集個數(shù)m小于停留點個數(shù)δ時,只保留該區(qū)間的首尾停留點。

    (2)當(dāng)坐標(biāo)點集個數(shù)m大于等于停留點個數(shù)δ時,并將停留點數(shù)目ns與δ作比,得出比值γ,保留γ倍數(shù)序號的停留點。其公式如下所示:

    (5)

    加入權(quán)重機制的FMDP算法(FMDP algorithm with Weight mechanism added,后文稱為“FMDPW算法”)雖然提高了算法的復(fù)雜度,但處理停留時間較長的漁船作業(yè)行為更加精準(zhǔn);反之,處理停留時間較短的漁船作業(yè)行為時,未加入權(quán)重機制的FMDP算法更加快速準(zhǔn)確,提高了軌跡壓縮的效率,完成了對漁船不同作業(yè)類型軌跡的高效壓縮。其算法流程由圖3所示,圖3中i代表軌跡點從隊首到隊尾的序號;di+1代表由經(jīng)緯度坐標(biāo)系轉(zhuǎn)化為笛卡爾直角坐標(biāo)系后的第i點與第i+1點的距離。

    圖3 算法流程

    2 試驗驗證與算法對比

    2.1 參數(shù)設(shè)定

    為了對本研究提出的FMDP軌跡壓縮方法的準(zhǔn)確性和執(zhí)行的效率進行驗證。采用 Python 語言編寫,選取來自大連世想海洋科技有限公司的原始試驗數(shù)據(jù)集起止時間為2018年1月20日到2018年12月20日,在壓縮效率方面與原始DP算法性能進行對比評估。

    評估指標(biāo):以壓縮效率作為試驗結(jié)果參考,對FMDP與DP算法進行評測,評測公式如下所示:

    (6)

    式中:Cp為壓縮率;Vp為剩余軌跡點數(shù);O為原始軌跡點數(shù)。

    2.2 壓縮效率對比試驗

    對漁船作業(yè)行為進行有效監(jiān)控的前提是完成漁船軌跡數(shù)據(jù)信息的高效壓縮處理,如通過處理壓縮軌跡中轉(zhuǎn)向與停留的特征信息,可以較為準(zhǔn)確判定漁船是否有越界、騙保等違規(guī)行為;通過數(shù)據(jù)壓縮可以刪除原始軌跡中冗雜、錯誤的軌跡點。但是由于大量數(shù)據(jù)被壓縮難免出現(xiàn)重要軌跡點缺失的問題,例如:魷魚釣漁船在進行捕撈作業(yè)時長時間處于漂流狀態(tài),會產(chǎn)生的大量漂流軌跡點,之前軌跡壓縮方法會因為數(shù)據(jù)重復(fù)將大量的停留作業(yè)軌跡點刪除,導(dǎo)致漁船作業(yè)信息丟失。因此為了實現(xiàn)漁船在不同作業(yè)情況均能高效地進行軌跡壓縮,本研究依據(jù)漁船軌跡估計的復(fù)雜度在算法中加入權(quán)重機制,以此實現(xiàn)較為完整的保留漁船軌跡的特征信息。

    試驗一:FMDP算法與DP算法對比相同漁船在不同日期的軌跡壓縮情況,通過軌跡的壓縮效果,來測試FMDP算法的壓縮性能。

    (1)軌跡壓縮效果情況對比

    在數(shù)據(jù)集中選取一組較為復(fù)雜軌跡數(shù)據(jù),使用FMDP算法與DP算法進行壓縮效果對比評估。如圖4所示。通過計算FMDP算法與DP算法的軌跡壓縮率Cp均超過98%,可知DP算法與FMDP算法均能實現(xiàn)對軌跡的有效壓縮,但是通過對比壓縮情況可以得出,當(dāng)漁船作出停留轉(zhuǎn)向等行為出現(xiàn)造軌跡折返時,DP算法在進行軌跡壓縮時可能會刪除速度變化較大的動態(tài)信息點,進而導(dǎo)致軌跡的時間特征的丟失[30],而FMDP算法會保留漁船轉(zhuǎn)向、停留的重要特征信息(圖4c中畫圈處),在確保高壓縮率的基礎(chǔ)上,保留了更多的漁船軌跡的特征信息。

    圖4 12月7日軌跡壓縮對比情況

    (2)壓縮率對比試驗

    為了進一步對比FMDP算法與DP算法的壓縮效果,進行了壓縮率的對比試驗,如表1所示,DP算法的壓縮率Cp均值為98.0%,FMDP的壓縮率Cp均值為97.5%,最低為95.51%。雖然FMDP算法在壓縮率上略低于DP算法,但依據(jù)圖4得出結(jié)論,DP算法所刪除的點集,大多屬于轉(zhuǎn)向點與停留點被刪除,無法反映軌跡的時間特征。

    表1 兩種壓縮算法的試驗結(jié)果

    對于軌跡的整體規(guī)劃不理想。FMDP算法則是在確保壓縮率Cp平均達到97.5%的基礎(chǔ)上,保證了壓縮后軌跡具備更多特征信息。此外,通過對比兩種算法的軌跡壓縮消耗的時間可以看出,FMDP算法處理軌跡壓縮花費時間與DP算法所消耗的時間相近,可見算法處理數(shù)據(jù)的壓縮時間不會對試驗結(jié)果評估造成影響。

    試驗二:為判斷在FMDP算法中加入權(quán)重機制,是否在處理所有的軌跡時都具有必要性。本試驗將待壓縮軌跡分為A,B兩組:A組為復(fù)雜漁船軌跡;B組中為簡單漁船軌跡。進行壓縮效率的計算,其具體結(jié)果如表2所示。

    表2 權(quán)重對壓縮率影響的試驗結(jié)果部分樣例

    通過對比表2中的A、B兩組數(shù)據(jù)中的壓縮率變化可以得出以下結(jié)論:FMDPW算法的壓縮率雖然要低于FMDP的壓縮率,但保留了原始軌跡中更多的特征信息,且原始軌跡復(fù)雜程度與保留的特征信息數(shù)量成正比。當(dāng)軌跡復(fù)雜需要保留更多特征信息時,優(yōu)先考慮使用FMDPW算法處理原始軌跡。為進一步證明權(quán)重機制的必要性,直觀反應(yīng)軌跡壓縮后的變化情況。試驗篩選出A、B兩組中有代表性的數(shù)據(jù)進行數(shù)據(jù)可視化,其對比結(jié)果如下圖5、圖6所示。

    圖5 A組復(fù)雜軌跡壓縮率對比圖

    圖6 B組簡單軌跡壓縮率對比圖

    通過對比圖5中的數(shù)據(jù)可以看出A組數(shù)據(jù)在加入權(quán)重機制處理后,更完整地保留了原始軌跡的變化特征,可以得出結(jié)論:FMDPW算法可以有效地從船舶軌跡數(shù)據(jù)中篩選出體現(xiàn)原始軌跡特征的子集,保留原始軌跡的特征信息。雖然FMDPW算法保留了原始軌跡中更多的特征信息,但導(dǎo)致了漁船軌跡的壓縮率和壓縮效率降低。并且根據(jù)圖6可以看出,簡單的漁船軌跡在FMDPW算法壓縮后軌跡并沒有明顯變化,因此在處理簡單的漁船軌跡時,無需在算法中加入壓縮權(quán)重機制。根據(jù)原始軌跡的復(fù)雜程度選擇是否進行權(quán)重考慮,可以最大程度上保證軌跡壓縮的高效性與準(zhǔn)確性。

    通過試驗對比與壓縮率對比,可以得出FMDP算法在處理漁船軌跡時能精準(zhǔn)保留漁船的特征信息。同時FMDP算法權(quán)重機制的考慮有效避免停留點過多造成壓縮軌跡失真的問題。在保證高壓縮率的情況下保留軌跡特征信息。

    3 結(jié)論

    本研究基于原始DP算法設(shè)計了一種特征點分離的多閾值DP算法,并加入權(quán)重機制的優(yōu)化研究。

    本研究基于原始算法首次提出了將軌跡中的停留點進行分離處理的方法,并通過權(quán)重機制以及加入角度閾值改進算法,最大程度保留了漁船軌跡的特征信息,為后續(xù)計算捕撈魚量,統(tǒng)計漁場捕撈強度提供數(shù)據(jù)支持,更能有效完成漁船作業(yè)行為分析,但漁船停留點閾值上選取需要進一步驗證。后續(xù)將會研究算法的最佳距離閾值與尺度變化間的定量關(guān)系,將著重考慮閾值的影響。

    猜你喜歡
    壓縮率漁船軌跡
    漁船
    千舟競發(fā)
    廉政瞭望(2021年15期)2021-08-23 04:55:22
    軌跡
    軌跡
    國內(nèi)新型遠洋金槍魚圍網(wǎng)漁船首航
    商周刊(2018年19期)2018-12-06 09:49:41
    水密封連接器尾部接電纜的優(yōu)化設(shè)計
    纏繞墊片產(chǎn)品質(zhì)量控制研究
    軌跡
    進化的軌跡(一)——進化,無盡的適應(yīng)
    中國三峽(2017年2期)2017-06-09 08:15:29
    多載波通信系統(tǒng)中CQI無損壓縮法研究
    亚洲av美国av| 在线看三级毛片| 欧洲精品卡2卡3卡4卡5卡区| 欧美性长视频在线观看| 成人三级黄色视频| 亚洲欧美一区二区三区黑人| 性欧美人与动物交配| 波多野结衣巨乳人妻| 国产一级毛片七仙女欲春2| 久久精品国产亚洲av香蕉五月| 午夜福利免费观看在线| 桃色一区二区三区在线观看| 亚洲专区字幕在线| 国产av在哪里看| 久久久久久久午夜电影| 成人手机av| 伦理电影免费视频| 国产午夜福利久久久久久| 婷婷六月久久综合丁香| 一级片免费观看大全| 欧美一级毛片孕妇| 不卡av一区二区三区| www.自偷自拍.com| 欧美性猛交黑人性爽| 一级作爱视频免费观看| 波多野结衣巨乳人妻| 一进一出抽搐gif免费好疼| 欧美黄色片欧美黄色片| 午夜免费激情av| 丰满人妻熟妇乱又伦精品不卡| 久久久久国产精品人妻aⅴ院| 亚洲乱码一区二区免费版| 欧美日韩中文字幕国产精品一区二区三区| 99热6这里只有精品| 日日爽夜夜爽网站| 人人妻,人人澡人人爽秒播| 亚洲国产欧美人成| 国产亚洲欧美98| 色综合婷婷激情| 91字幕亚洲| 黄色视频不卡| 国产精品美女特级片免费视频播放器 | 日韩高清综合在线| 国产成人av激情在线播放| 在线观看免费视频日本深夜| 亚洲18禁久久av| 日本一二三区视频观看| 怎么达到女性高潮| 51午夜福利影视在线观看| 午夜激情福利司机影院| 无限看片的www在线观看| 亚洲精品一区av在线观看| 国产一区二区三区在线臀色熟女| 又爽又黄无遮挡网站| 好看av亚洲va欧美ⅴa在| 97碰自拍视频| 亚洲欧美一区二区三区黑人| 欧美另类亚洲清纯唯美| 久久香蕉精品热| 香蕉av资源在线| 久久精品成人免费网站| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 黄片小视频在线播放| 久久久久九九精品影院| 欧美最黄视频在线播放免费| 欧美日本亚洲视频在线播放| 亚洲黑人精品在线| 人妻丰满熟妇av一区二区三区| 老司机靠b影院| 老汉色av国产亚洲站长工具| x7x7x7水蜜桃| 中文字幕av在线有码专区| 亚洲 国产 在线| 久久久久久大精品| 人妻久久中文字幕网| 岛国在线观看网站| 免费电影在线观看免费观看| 国产精品美女特级片免费视频播放器 | 久久精品国产综合久久久| 日本五十路高清| 99国产综合亚洲精品| 亚洲成a人片在线一区二区| 国产成+人综合+亚洲专区| 亚洲国产精品sss在线观看| 国产黄色小视频在线观看| 久久久久国产一级毛片高清牌| 亚洲欧美日韩东京热| 99re在线观看精品视频| 最近最新免费中文字幕在线| 精品久久久久久久久久免费视频| 此物有八面人人有两片| 亚洲黑人精品在线| 成人18禁高潮啪啪吃奶动态图| 成人特级黄色片久久久久久久| 欧美黄色片欧美黄色片| 国产视频内射| av免费在线观看网站| 国产精品免费一区二区三区在线| 中文字幕精品亚洲无线码一区| 9191精品国产免费久久| 午夜福利18| 久久精品国产亚洲av高清一级| 特大巨黑吊av在线直播| 亚洲午夜精品一区,二区,三区| 一进一出好大好爽视频| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 天天躁狠狠躁夜夜躁狠狠躁| 丝袜人妻中文字幕| 淫妇啪啪啪对白视频| 成人av一区二区三区在线看| 天天添夜夜摸| 成人18禁在线播放| 成人国产一区最新在线观看| 国内毛片毛片毛片毛片毛片| 正在播放国产对白刺激| 一进一出抽搐gif免费好疼| 黄色视频,在线免费观看| 少妇的丰满在线观看| 亚洲色图av天堂| 日韩欧美一区二区三区在线观看| 91国产中文字幕| 可以在线观看毛片的网站| 亚洲av成人不卡在线观看播放网| 国产精品 国内视频| 观看免费一级毛片| 2021天堂中文幕一二区在线观| 51午夜福利影视在线观看| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| 亚洲精品国产精品久久久不卡| x7x7x7水蜜桃| 亚洲成人免费电影在线观看| 女生性感内裤真人,穿戴方法视频| 精品无人区乱码1区二区| 亚洲专区字幕在线| 亚洲18禁久久av| 老司机午夜十八禁免费视频| 亚洲人成网站在线播放欧美日韩| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 老汉色∧v一级毛片| 国产一区二区在线av高清观看| x7x7x7水蜜桃| 日韩中文字幕欧美一区二区| 亚洲狠狠婷婷综合久久图片| 熟女少妇亚洲综合色aaa.| 久久香蕉精品热| 国产精品日韩av在线免费观看| 国产精品亚洲av一区麻豆| 白带黄色成豆腐渣| 99精品欧美一区二区三区四区| 国产麻豆成人av免费视频| 天天一区二区日本电影三级| 欧美黑人欧美精品刺激| 狂野欧美激情性xxxx| 三级男女做爰猛烈吃奶摸视频| 亚洲熟妇中文字幕五十中出| 亚洲av电影不卡..在线观看| 国产精品综合久久久久久久免费| 99久久无色码亚洲精品果冻| 精品高清国产在线一区| 香蕉av资源在线| 一进一出抽搐动态| 精品少妇一区二区三区视频日本电影| 国产精品自产拍在线观看55亚洲| 国产精品久久电影中文字幕| 久久久久亚洲av毛片大全| 两个人视频免费观看高清| 91大片在线观看| 国内少妇人妻偷人精品xxx网站 | svipshipincom国产片| 人妻夜夜爽99麻豆av| 国产精品精品国产色婷婷| 成人18禁在线播放| 亚洲乱码一区二区免费版| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 免费av毛片视频| 久99久视频精品免费| 熟妇人妻久久中文字幕3abv| 1024香蕉在线观看| 成人国语在线视频| 亚洲成人免费电影在线观看| 最近最新中文字幕大全免费视频| 久久久久久久久免费视频了| 色在线成人网| 两性夫妻黄色片| 国产熟女xx| 免费电影在线观看免费观看| 99久久综合精品五月天人人| 色哟哟哟哟哟哟| 91av网站免费观看| 欧美一级a爱片免费观看看 | bbb黄色大片| 99久久99久久久精品蜜桃| 国内揄拍国产精品人妻在线| e午夜精品久久久久久久| 欧美日韩精品网址| 黄色视频,在线免费观看| 日韩欧美国产在线观看| 少妇人妻一区二区三区视频| 给我免费播放毛片高清在线观看| 亚洲一区中文字幕在线| 一个人免费在线观看的高清视频| 夜夜躁狠狠躁天天躁| 亚洲欧美日韩无卡精品| 亚洲人成电影免费在线| 午夜两性在线视频| 日韩av在线大香蕉| 欧美日本视频| 变态另类丝袜制服| 最新美女视频免费是黄的| 国产亚洲精品第一综合不卡| 最近在线观看免费完整版| 免费在线观看视频国产中文字幕亚洲| 国产一区二区在线av高清观看| 99热这里只有精品一区 | 成人三级黄色视频| 精品无人区乱码1区二区| 国产精品一区二区三区四区免费观看 | 国产激情偷乱视频一区二区| 黄色成人免费大全| 黑人操中国人逼视频| 久久欧美精品欧美久久欧美| 欧美一区二区国产精品久久精品 | 亚洲精品国产一区二区精华液| 手机成人av网站| 一区福利在线观看| 免费看美女性在线毛片视频| 国产av又大| 国产精品 欧美亚洲| 黄色成人免费大全| 亚洲美女视频黄频| a级毛片在线看网站| 久久精品国产99精品国产亚洲性色| 日韩欧美在线乱码| 久久久国产欧美日韩av| 国产亚洲欧美在线一区二区| 午夜激情福利司机影院| 亚洲av片天天在线观看| 国产精品影院久久| 少妇裸体淫交视频免费看高清 | 手机成人av网站| 男女做爰动态图高潮gif福利片| 一个人观看的视频www高清免费观看 | 99在线人妻在线中文字幕| 三级毛片av免费| 国产激情欧美一区二区| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 三级男女做爰猛烈吃奶摸视频| АⅤ资源中文在线天堂| 亚洲一区高清亚洲精品| 久久久久久九九精品二区国产 | 好看av亚洲va欧美ⅴa在| 草草在线视频免费看| 国内精品久久久久久久电影| 成人av一区二区三区在线看| 一本一本综合久久| 在线观看免费视频日本深夜| 欧美成人性av电影在线观看| 欧美乱色亚洲激情| or卡值多少钱| 成人欧美大片| 天堂√8在线中文| 精品久久久久久久久久久久久| 久久久国产欧美日韩av| 两人在一起打扑克的视频| 妹子高潮喷水视频| 午夜免费激情av| 亚洲午夜精品一区,二区,三区| 天天添夜夜摸| 亚洲国产欧美网| 欧美黄色片欧美黄色片| 国产又黄又爽又无遮挡在线| 亚洲美女视频黄频| www.精华液| 免费看a级黄色片| 好男人在线观看高清免费视频| 亚洲色图av天堂| 在线观看美女被高潮喷水网站 | 欧美黑人欧美精品刺激| 91大片在线观看| 国产主播在线观看一区二区| 欧美成人免费av一区二区三区| 成人国产一区最新在线观看| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 两性夫妻黄色片| 久热爱精品视频在线9| 少妇裸体淫交视频免费看高清 | 国产成年人精品一区二区| 日韩欧美在线二视频| 丁香欧美五月| 久久久国产成人精品二区| 久久久久久九九精品二区国产 | 最好的美女福利视频网| 国产日本99.免费观看| 黄色 视频免费看| 久久精品成人免费网站| 麻豆成人av在线观看| 69av精品久久久久久| 黑人操中国人逼视频| 久久香蕉国产精品| 国产1区2区3区精品| 欧美黑人欧美精品刺激| 国产成人av激情在线播放| 久久天堂一区二区三区四区| 国产人伦9x9x在线观看| 少妇被粗大的猛进出69影院| 久久伊人香网站| 国产一区在线观看成人免费| 欧美3d第一页| 村上凉子中文字幕在线| 夜夜夜夜夜久久久久| 国内少妇人妻偷人精品xxx网站 | 亚洲成av人片在线播放无| 麻豆国产97在线/欧美 | 国产91精品成人一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲成人中文字幕在线播放| 久久草成人影院| 人人妻人人澡欧美一区二区| 日本撒尿小便嘘嘘汇集6| 久久久久久人人人人人| 一级作爱视频免费观看| 老熟妇仑乱视频hdxx| 成人手机av| 中文亚洲av片在线观看爽| 久久天躁狠狠躁夜夜2o2o| 欧美不卡视频在线免费观看 | 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 精品欧美一区二区三区在线| 丁香欧美五月| 熟女少妇亚洲综合色aaa.| 啦啦啦免费观看视频1| 成人午夜高清在线视频| 亚洲午夜精品一区,二区,三区| 黄频高清免费视频| 亚洲欧美日韩东京热| 一夜夜www| 亚洲自偷自拍图片 自拍| or卡值多少钱| 国产精品亚洲美女久久久| 俄罗斯特黄特色一大片| 一本综合久久免费| 日日干狠狠操夜夜爽| 超碰成人久久| 最近最新中文字幕大全免费视频| 白带黄色成豆腐渣| 免费看十八禁软件| 在线观看www视频免费| a在线观看视频网站| 精品少妇一区二区三区视频日本电影| 亚洲熟妇熟女久久| 波多野结衣巨乳人妻| 久久精品综合一区二区三区| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲| 欧美乱色亚洲激情| 99国产综合亚洲精品| 91九色精品人成在线观看| 中文字幕熟女人妻在线| 香蕉久久夜色| 18禁美女被吸乳视频| 搞女人的毛片| 啪啪无遮挡十八禁网站| 国产精品免费视频内射| 禁无遮挡网站| 老熟妇仑乱视频hdxx| 国产成人aa在线观看| 伊人久久大香线蕉亚洲五| 午夜老司机福利片| 人人妻人人澡欧美一区二区| 婷婷丁香在线五月| 在线永久观看黄色视频| 欧美久久黑人一区二区| 黄色视频不卡| 韩国av一区二区三区四区| 黄色视频,在线免费观看| 好看av亚洲va欧美ⅴa在| 少妇裸体淫交视频免费看高清 | 亚洲精品国产精品久久久不卡| 国产av麻豆久久久久久久| 日韩三级视频一区二区三区| 国产激情偷乱视频一区二区| 久久中文字幕一级| 激情在线观看视频在线高清| 久久久久国内视频| 国产一区二区三区视频了| 制服诱惑二区| 此物有八面人人有两片| 精品久久蜜臀av无| 久久婷婷人人爽人人干人人爱| svipshipincom国产片| 香蕉av资源在线| 亚洲人成网站高清观看| 久久久久久久午夜电影| 欧美绝顶高潮抽搐喷水| 免费人成视频x8x8入口观看| 黄色成人免费大全| tocl精华| 悠悠久久av| 美女 人体艺术 gogo| 婷婷精品国产亚洲av在线| 夜夜夜夜夜久久久久| 精华霜和精华液先用哪个| 色av中文字幕| 黄片小视频在线播放| 国产探花在线观看一区二区| 九九热线精品视视频播放| 亚洲一区二区三区不卡视频| 99久久无色码亚洲精品果冻| 在线观看www视频免费| 成人av在线播放网站| 两个人的视频大全免费| x7x7x7水蜜桃| 在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 国产精品一区二区三区四区久久| 天堂av国产一区二区熟女人妻 | 国产av麻豆久久久久久久| 亚洲人成77777在线视频| 搡老熟女国产l中国老女人| 一个人免费在线观看电影 | 少妇被粗大的猛进出69影院| 50天的宝宝边吃奶边哭怎么回事| 一个人免费在线观看的高清视频| 757午夜福利合集在线观看| 久久精品91无色码中文字幕| 成人av一区二区三区在线看| 国产69精品久久久久777片 | 一本久久中文字幕| 午夜福利高清视频| 欧美3d第一页| 免费av毛片视频| 中文字幕久久专区| 香蕉丝袜av| 亚洲国产精品合色在线| 在线观看免费日韩欧美大片| 此物有八面人人有两片| 成人亚洲精品av一区二区| 午夜免费观看网址| 精品久久久久久久人妻蜜臀av| 成在线人永久免费视频| 国产黄片美女视频| 一本一本综合久久| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 女生性感内裤真人,穿戴方法视频| 男男h啪啪无遮挡| 黄频高清免费视频| 午夜两性在线视频| 美女 人体艺术 gogo| 精华霜和精华液先用哪个| 日韩欧美国产在线观看| 久久久久久久午夜电影| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 999久久久国产精品视频| 夜夜夜夜夜久久久久| 欧美黑人精品巨大| 国内精品久久久久精免费| 天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成av人片免费观看| 99国产极品粉嫩在线观看| 91麻豆av在线| 亚洲,欧美精品.| 一区二区三区国产精品乱码| 亚洲七黄色美女视频| 亚洲专区国产一区二区| 午夜精品久久久久久毛片777| 村上凉子中文字幕在线| 一本精品99久久精品77| 一个人免费在线观看的高清视频| 成年版毛片免费区| 成人av在线播放网站| 天天躁狠狠躁夜夜躁狠狠躁| 最新美女视频免费是黄的| 欧美性猛交黑人性爽| 日韩精品中文字幕看吧| 好男人在线观看高清免费视频| 亚洲精品一区av在线观看| 日本三级黄在线观看| 变态另类丝袜制服| 亚洲欧美日韩东京热| 91国产中文字幕| 亚洲真实伦在线观看| 亚洲男人的天堂狠狠| 亚洲精品色激情综合| 韩国av一区二区三区四区| 精品久久久久久久末码| 久久久久久九九精品二区国产 | 国产私拍福利视频在线观看| 精品高清国产在线一区| 亚洲国产欧美一区二区综合| 亚洲专区中文字幕在线| 美女扒开内裤让男人捅视频| 99热只有精品国产| 国产午夜精品论理片| 亚洲精品久久国产高清桃花| www.自偷自拍.com| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉国产精品| 欧美日韩精品网址| 国产麻豆成人av免费视频| 午夜a级毛片| 亚洲 欧美一区二区三区| 一夜夜www| 亚洲中文字幕一区二区三区有码在线看 | 色av中文字幕| 国产黄a三级三级三级人| 欧美绝顶高潮抽搐喷水| 两性夫妻黄色片| 久久久国产欧美日韩av| 国产在线精品亚洲第一网站| 久久中文字幕人妻熟女| 在线观看免费视频日本深夜| 免费在线观看亚洲国产| 国产蜜桃级精品一区二区三区| 精品久久久久久久人妻蜜臀av| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 在线播放国产精品三级| 亚洲性夜色夜夜综合| 黑人巨大精品欧美一区二区mp4| 日本一区二区免费在线视频| 亚洲片人在线观看| 看黄色毛片网站| 最近最新免费中文字幕在线| 嫩草影视91久久| 亚洲第一电影网av| 亚洲精品中文字幕一二三四区| 国产99白浆流出| 老司机午夜福利在线观看视频| 国产高清视频在线观看网站| 特级一级黄色大片| 好男人在线观看高清免费视频| 五月伊人婷婷丁香| 成人一区二区视频在线观看| 欧美在线一区亚洲| 嫩草影视91久久| 正在播放国产对白刺激| av在线天堂中文字幕| 欧美成人午夜精品| 国产精品亚洲av一区麻豆| 欧美黑人欧美精品刺激| 亚洲国产精品合色在线| 女警被强在线播放| 国产免费av片在线观看野外av| 亚洲午夜理论影院| 99riav亚洲国产免费| 午夜影院日韩av| xxx96com| 岛国在线免费视频观看| 大型黄色视频在线免费观看| 亚洲国产精品久久男人天堂| a级毛片a级免费在线| 全区人妻精品视频| 性欧美人与动物交配| 国产av不卡久久| 日本熟妇午夜| 99精品久久久久人妻精品| 久久这里只有精品中国| netflix在线观看网站| 久久性视频一级片| 国产又色又爽无遮挡免费看| 一级毛片精品| 亚洲国产欧美一区二区综合| 婷婷精品国产亚洲av在线| 日本一本二区三区精品| 日本a在线网址| 12—13女人毛片做爰片一| 日本在线视频免费播放| 久久中文字幕一级| 悠悠久久av| 欧美国产日韩亚洲一区| 香蕉丝袜av| 午夜视频精品福利| 国产久久久一区二区三区| 极品教师在线免费播放| 人妻夜夜爽99麻豆av| 国产久久久一区二区三区| 99国产综合亚洲精品| 我的老师免费观看完整版| 久久伊人香网站| 精品国产乱码久久久久久男人| 成人av一区二区三区在线看| 国产av一区在线观看免费| 一个人免费在线观看电影 | 淫秽高清视频在线观看| 欧美三级亚洲精品| 美女扒开内裤让男人捅视频| 国产精品自产拍在线观看55亚洲| 日本黄大片高清| 亚洲国产精品久久男人天堂| 欧美日韩亚洲综合一区二区三区_| 999久久久精品免费观看国产| 国产激情久久老熟女| 99国产综合亚洲精品| 老司机靠b影院| 在线播放国产精品三级| 18禁美女被吸乳视频| 嫩草影视91久久| 91老司机精品| a级毛片在线看网站| 日韩欧美精品v在线| 亚洲av中文字字幕乱码综合| 日韩大码丰满熟妇| 亚洲成人中文字幕在线播放| 亚洲国产欧美网|