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

    載體軌跡停留信息提取的核密度法及其可視化

    2016-11-07 02:37:36向隆剛邵曉天
    測繪學(xué)報 2016年9期
    關(guān)鍵詞:鄰域時空軌跡

    向隆剛,邵曉天

    1. 武漢大學(xué)測繪遙感信息工程國家重點實驗室,湖北 武漢430079; 2. 地球空間信息技術(shù)協(xié)同創(chuàng)新中心,湖北 武漢 430079

    ?

    載體軌跡停留信息提取的核密度法及其可視化

    向隆剛1,2,邵曉天1,2

    1. 武漢大學(xué)測繪遙感信息工程國家重點實驗室,湖北 武漢430079; 2. 地球空間信息技術(shù)協(xié)同創(chuàng)新中心,湖北 武漢 430079

    軌跡停留蘊含重要語義信息,其有效提取是開展軌跡Stop/Move模型分析的前提。本文首先依據(jù)核密度思想,通過累計鄰域點時空貢獻(xiàn)來定義軌跡點的停留指數(shù),在此基礎(chǔ)上設(shè)計了停留指數(shù)圖,以圖形方式直觀表達(dá)軌跡點的時空聚集程度變化。進(jìn)一步針對源于停留指數(shù)的潛在停留段,提出了一種基于潛在停留段時空臨近關(guān)系的逐級合并算法,以自動發(fā)現(xiàn)和提取停留。試驗表明,該算法兼顧停留識別的完整性和準(zhǔn)確性,可以有效識別復(fù)雜多樣的軌跡停留,即使面對噪聲嚴(yán)重的軌跡,停留提取的正確率依然較高。

    核密度;軌跡;停留指數(shù);停留;時空貢獻(xiàn)

    隨著GPS、無線通信和移動設(shè)備的結(jié)合與發(fā)展,時空軌跡數(shù)據(jù)獲取技術(shù)突飛猛進(jìn),相關(guān)行業(yè)已經(jīng)累積了海量的時空軌跡數(shù)據(jù)。如何分析挖掘這些時空軌跡數(shù)據(jù),從中提取信息與知識,已逐漸受到國內(nèi)外學(xué)者的重視,成為空間信息領(lǐng)域的一大研究熱點。軌跡Stop/Move模型是近年來提出的一種較為新穎的軌跡分析方法[1-2],該方法首先將軌跡點序列劃分為Stop/Move序列,并賦予相關(guān)語義,在此基礎(chǔ)上可以開展語義更為明確的分析挖掘研究[3-4]。Stop,即停留,指在移動對象的軌跡中靜止或在一定范圍內(nèi)緩慢移動的部分,是軌跡中具有重要語義特征的點序列。一方面,停留提取是支持面向軌跡數(shù)據(jù)的知識發(fā)現(xiàn)與行為識別的重要手段,如從候鳥軌跡中發(fā)現(xiàn)遷徙過程中的停歇地,從漁船軌跡中分析出海過程中的捕魚行為;另一方面,從軌跡停留出發(fā),研究大量軌跡在停留-移動序列上的相似性,可以進(jìn)一步發(fā)掘移動對象的時空移動模式[5-6],從而為城市交通、居民出行提供決策支持。以出租車軌跡數(shù)據(jù)為例,在停留提取的基礎(chǔ)上,分析出租車司機的午餐行為在空間和時間上的分布,可以輔助政府設(shè)計合理的出租車就餐點,從而提高出租車的服務(wù)效率,并方便居民出行;又如微軟的GeoLife項目[7],通過收集游客軌跡數(shù)據(jù),并分析挖掘其中關(guān)于停留-移動的頻繁模式,在旅游線路推薦和朋友圈發(fā)現(xiàn)等方面進(jìn)行了一些有益嘗試。為此,國內(nèi)外學(xué)者提出了多種停留提取方法,大致可以分為集成地理背景信息的方法[8-9]、基于軌跡導(dǎo)出特征的方法[10-11]、密度聚類法[12-14]以及逐級合并方法[15]。

    集成地理背景信息的方法是指在軌跡停留提取過程中,從地理背景信息出發(fā),通過檢測移動對象相對于不同地理實體(POI、ROI)的時空關(guān)系來識別軌跡停留。具體來說,這類方法通過計算軌跡在地理實體內(nèi)部或其鄰域范圍內(nèi)的累積連續(xù)停留時間,來判斷軌跡在該處是否存在停留。該類方法適用于在已有地理背景數(shù)據(jù)環(huán)境下查詢或挖掘軌跡與地理實體間聯(lián)系的應(yīng)用,但無法從地理背景信息缺失的區(qū)域提取出停留[9]。

    基于軌跡導(dǎo)出特征的方法通過在軌跡中諸如速度、加速度、點密度、轉(zhuǎn)向角、信號缺失等一個或多個特征上設(shè)置規(guī)則,已從軌跡數(shù)據(jù)中提取出停留。如文獻(xiàn)[10]通過設(shè)置累積轉(zhuǎn)向角閾值來發(fā)現(xiàn)漁船軌跡中方向頻繁變化的捕魚點;文獻(xiàn)[11]綜合利用信號在缺失時間、轉(zhuǎn)向角多上的多個規(guī)則來挖掘居民出行中的活動地點。由于軌跡停留識別的最佳特征選取往往取決于特定的數(shù)據(jù)或應(yīng)用,使得該類方法的通用性受到一定程度的制約。

    密度聚類法通過在傳統(tǒng)的空間密度聚類方法中引入時間維信息,從軌跡中提取時空兩方面同時聚集的停留。如基于DBSCAN算法的改進(jìn)算法CB-SMoT[12]、TrajDBSCAN[13],該類方法針對軌跡的時空特征,重新定義了如鄰居、核心點、密度可達(dá)、密度相連等DBSCAN算法中的基本概念;又如文獻(xiàn)[14]對OPTICS中的可達(dá)距離在軌跡數(shù)據(jù)上進(jìn)行重新定義,設(shè)計了適用于軌跡的停留提取算法T-OPTICS。該類方法面向軌跡數(shù)據(jù),對傳統(tǒng)密度聚類算法進(jìn)行了相應(yīng)擴展,但在處理帶長距離漂移的軌跡數(shù)據(jù)時,其停留識別的效果不佳,易將單個停留誤識別為多個停留[14-15]。

    為了克服漂移噪聲對停留識別的影響,文獻(xiàn)[15]針對居民出行軌跡提出了逐級合并的停留識別方法,該方法首先使用簡單的速度判別準(zhǔn)則,將軌跡初步分為停留/移動類型子軌跡,并根據(jù)有意義的停留與移動在持續(xù)時間、跨越距離上的限制,動態(tài)更新子軌跡的停留/移動標(biāo)簽,并通過逐級合并相同類型相鄰的子軌跡來優(yōu)化識別結(jié)果。該方法雖然在一定程度上克服了噪聲對停留識別的影響,但是其簡單合并策略易于導(dǎo)致不同停留的誤合并。

    通過對現(xiàn)有停留提取研究的總結(jié),筆者發(fā)現(xiàn):①現(xiàn)有研究往往側(cè)重于停留的定義或停留識別準(zhǔn)則的選取,對表征停留本身的時空聚集程度及其可視化方法卻較少研究;②現(xiàn)有研究較少考慮含大量漂移噪聲的室內(nèi)停留的優(yōu)化處理問題,對軌跡停留的復(fù)雜多樣性考慮不足。為此,本文基于核密度思想提出停留指數(shù),并據(jù)此設(shè)計停留指數(shù)圖,直觀表達(dá)軌跡點的時空聚集程度變化,在此基礎(chǔ)上,本文開發(fā)停留提取算法,并采用自采軌跡和GeoLife軌跡[16]開展試驗,結(jié)果表明該算法在停留提取方面具有較高的正確率,能夠識別復(fù)雜多樣的軌跡停留。

    1 軌跡及其停留

    定義1(軌跡):traj=(tid,),其中,pi=(xi,yi,ti),0≤i≤N,且?i

    tid是軌跡的唯一標(biāo)識,pi=(xi,yi,ti)是構(gòu)成軌跡的1個時空點,其中:xi和yi是二維歐氏空間中的坐標(biāo)值(通常為經(jīng)緯度),ti是時間全序軸上的時刻值,xi、yi和ti均來自于實數(shù)域。p0和pN是兩個特殊時空點,即始點和終點。

    軌跡在其生命周期內(nèi)并不是一直移動,可能出于某一目的在某一地點停留一段時間(如就餐和加油等)。在該段時間內(nèi),軌跡的空間位置保持不變或在較小范圍內(nèi)變化;此后,軌跡繼續(xù)在空間中移動,并到達(dá)下一地點停留,或者抵達(dá)終點。停留,即Stop,是軌跡點序列中的一個連續(xù)子序列:

    定義2(停留):stop=(tid,sid,),其中,0

    由定義2可知,停留不是由絕對速度定義而來,而是在時空兩維中,由空間范圍和時間長度兩者共同界定的概念。由于受到多種因素影響,軌跡停留呈現(xiàn)出復(fù)雜多樣的特征:從發(fā)生地點來看,可以是室內(nèi),也可以是室外;從外在表現(xiàn)來看,停留可以是絕對靜止、小范圍內(nèi)移動,或者兩者的結(jié)合;而從數(shù)據(jù)承載上來看,停留可以是單個軌跡點,多個或者大量連續(xù)的軌跡點。

    2 停留指數(shù)及其可視化

    核密度分析法[18-19]作為一種基于密度的空間點模式分析技術(shù),是研究空間點聚集性的有效方法。該方法通過對核密度的計算,得到核密度在空間上的分布,并據(jù)此提取具有重要語義的聚集模式,如交通事件、犯罪熱點等[20-21]。通常來說,移動對象接近停留點時速度逐步降低,停留期間,速度為零或者很低,而離開停留點后,速度將逐步提高,因而,從數(shù)據(jù)承載方面來看,停留表現(xiàn)為大量軌跡點在局部小范圍內(nèi)的聚集,據(jù)此,本文從核密度思想出發(fā),提出停留指數(shù),以反映軌跡在其生命周期內(nèi)時空聚集程度的變化,在此基礎(chǔ)上以可視化方式直觀揭示軌跡的“潛在”停留,并進(jìn)一步發(fā)展出軌跡停留的自動提取算法。

    2.1停留指數(shù)

    核密度法通過累計目標(biāo)點鄰域內(nèi)其他點對中心點的空間權(quán)重,估算目標(biāo)點處的密度,見式(1)

    (1)

    式中,n是落在目標(biāo)點s距離為h的鄰域內(nèi)的點數(shù);函數(shù)K是核函數(shù),其值隨著s與點Ci的距離的增大而遞減。K函數(shù)的性質(zhì)體現(xiàn)了空間影響隨距離衰減的地理學(xué)第一定律[22-23],是核密度思想的核心。

    由于記錄方式或信號質(zhì)量等原因,軌跡數(shù)據(jù)的采樣頻率并非是恒等不變的。因而,在將核密度法應(yīng)用于軌跡點停留指數(shù)的計算時,除了考慮鄰域內(nèi)的點數(shù)和空間距離之外,應(yīng)進(jìn)一步將點在鄰域內(nèi)的逗留時間納入進(jìn)來,以反映軌跡在目標(biāo)點領(lǐng)域內(nèi)(即空間維信息)的逗留時間(即時間維信息)。據(jù)此,軌跡點的停留指數(shù)定義如下

    (2)

    2.1.1逗留時間計算

    一種簡單直觀的計算ST(Pj)的方法是取Pj到Pj+1的時間采樣間隔,即ST(Pj)=interval(Pj,Pj+1)=Pj+1·tj+1-Pj·tj。由于軌跡是現(xiàn)實世界中對象的連續(xù)移動在計算機世界中的離散化表達(dá),上述方法將人為隔斷跨越鄰域的相鄰軌跡點之間的時空聯(lián)系,勢必造成邊緣點時空貢獻(xiàn)的誤計算。

    設(shè)想圖1(a)所示的情形,pj到pj+1的子軌跡跨越了N(pi,τ)邊界,而在實際情況中,pj和pj+1可能間隔很久(如經(jīng)過長距離隧道所造成的信號丟失),故以時間采樣間隔來計算鄰域點的時空貢獻(xiàn)將導(dǎo)致點pi停留指數(shù)的計算結(jié)果虛高??紤]到軌跡移動的連續(xù)性,在計算點pj的逗留時間時,僅需考慮子軌跡pj→pj+1落在N(pi,τ)內(nèi)的部分(在圖1中由實線表示)所對應(yīng)的時間,故停留時間的計算公式如下

    (3)

    式中,edin為子軌跡pj→pj+1在鄰域內(nèi)的長度。此處,筆者假設(shè)軌跡對象在兩個采樣點之間勻速運動。

    圖1 逗留時間計算與空間權(quán)重校正Fig.1 Calculation of stop time (ST) and its spatial weight correction

    2.1.2空間權(quán)重校正

    在更復(fù)雜的軌跡移動中,軌跡可能多次進(jìn)出鄰域,而每次進(jìn)出時引入的鄰域點對目標(biāo)點的時空貢獻(xiàn)顯然是不相同的。圖1(b)示意了一條軌跡連續(xù)兩次進(jìn)出N(pi,τ)的過程,按照時間順序是:enter1→leave1→enter2→leave2,形成兩段落在N(pi,τ)內(nèi)的子軌跡,即stay1和stay2。不同于stay1,stay2在鄰域內(nèi)不存在到目標(biāo)點pi的直接通路,故在計算鄰域點時空貢獻(xiàn)時,應(yīng)對stay1與stay2加以區(qū)分。不難得出,stay2上的鄰域點,設(shè)為pj,對目標(biāo)點pi的時空貢獻(xiàn)不僅取決于pj到pi的歐氏距離,還應(yīng)考慮軌跡在pi與pj之間的移動距離。據(jù)此,筆者對鄰域點pj時空貢獻(xiàn)的空間權(quán)重進(jìn)行校正,引入如下所示的校正因子

    (4)

    式中,n=min(i,j);m=max(i,j);TD(pn,pm)為軌跡從pn到pm的移動距離;tdout為其中不在鄰域內(nèi)的長度,稱作脫離距離。當(dāng)軌跡從pn移動到pm時未曾脫離鄰域,即脫離距離為0,那么SWP(Pj)=1,這意味著無需對pj時空貢獻(xiàn)進(jìn)行校正。

    校正因子依據(jù)鄰域點的脫離距離對其時空貢獻(xiàn)進(jìn)行矯正:對于噪聲導(dǎo)致的暫時脫離鄰域的情形,由于脫離距離較短,其對時空貢獻(xiàn)的影響可忽略不計;而對于離開-返回形成的脫離鄰域的情形,由于脫離距離通常較長,其時空貢獻(xiàn)將被調(diào)至較低水平。

    (5)

    需要說明的是,此處去掉常系數(shù)是為了更好地反映時空貢獻(xiàn)蘊含的語義,即保證目標(biāo)點pi本身在對于pi停留指數(shù)的時空貢獻(xiàn)中:空間權(quán)重為1,逗留時間則依賴于點pi+1的空間分布情況,具體由式(3)決定。

    由定義(2)可知,停留是移動對象在一定時長內(nèi)保持靜止或在小范圍內(nèi)移動的狀態(tài),相應(yīng)于該定義,停留指數(shù)在小范圍內(nèi)(領(lǐng)域半徑)計算移動對象累積的逗留時間。式(5)之所以采用高斯函數(shù),其原因在于為不同鄰域點的逗留時間賦予不同的空間衰減權(quán),從而反映出軌跡點的不同時空聚集程度,即鄰域點的空間分布越集中于目標(biāo)點附近,目標(biāo)點的停留指數(shù)則越高,這也符合人們對于停留的直觀認(rèn)識。不同于CB-SMoT等在擴展DBSCAN算法時先空間,后時間的做法,停留指數(shù)在累積每一個鄰域點的貢獻(xiàn)時,同時考慮其時空特征(即空間距離和逗留時間),不僅可以表征大量軌跡點聚集所形成的停留,而且能夠反映軌跡點較少(甚至僅有一個軌跡點)但持續(xù)時間較長的停留。后者一般對應(yīng)于移動對象在停留期間GPS設(shè)備降低采樣頻率、不采樣、甚至關(guān)機的情形,而CB-SMoT等方法在處理由此形成的停留時往往無能為力。

    由式(5)可知,停留指數(shù)僅由移動對象在鄰域內(nèi)的逗留時間決定,與速度并不發(fā)生直接的關(guān)系。設(shè)想如下情形:飛機在機場上空盤旋時,雖然速度很高,但始終徘徊于機場附近,故由機場這一鄰域決定的停留指數(shù)仍然較大。但在一般情況下,停留指數(shù)與速度存在一種間接的反比關(guān)系,比如移動對象通過鄰域的速度越大,所花費的時間越少,因而其停留指數(shù)越小,反之亦然。

    2.2停留指數(shù)圖

    以軌跡點為橫軸,停留指數(shù)為縱軸,即可繪出軌跡的停留指數(shù)圖。圖2是一條從常州自駕到溧陽的軌跡的停留指數(shù)圖,從上到下分別是τ=50 m和τ=1000 m的停留指數(shù)圖,以及相應(yīng)于τ=1000 m的熱力圖。為了顯示方便,圖2所示的停留指數(shù)圖在縱軸進(jìn)行了最大值截斷。借助于圖形化表達(dá)的停留指數(shù)圖,用戶不必解譯晦澀難懂的軌跡數(shù)據(jù),即可獲知軌跡在其生命周期中的移動-停留情況。以圖2為例,不難得出此次自駕之旅有3次明顯的停留行為,且停留③持續(xù)了較長時間。此外,同停留①相比,參與停留②的軌跡點雖少,但其停留指數(shù)較高。

    需要指出的是,停留是與空間尺度相關(guān)的概念,以長距離自駕為例,大尺度上可將在城鎮(zhèn)的一次落腳休整看作一個大停留,而縮放到小尺度之后,則可將這次落腳細(xì)分為多個小停留,如就餐、購物和住店等。不難得出,停留指數(shù)通過鄰域半徑τ建立起同空間尺度的關(guān)系,因而在不同大小的鄰域半徑下,可以觀察到不同尺度下的停留情況。繼續(xù)以圖2為例,當(dāng)τ增大到1000 m時,在τ=50 m時觀察到的前2小個停留被連接成1個大停留,即在大尺度下可將在常州市的2次停留看作1個停留。

    3 基于停留指數(shù)的停留提取

    從軌跡停留指數(shù)的可視化結(jié)果不難看出,實際停留對應(yīng)的子軌跡具有較高的停留指數(shù)值,這啟發(fā)筆者試著從停留指數(shù)發(fā)現(xiàn)并提取停留,為此,擬從軌跡中搜索高停留指數(shù)的軌跡點序列形成的子軌跡,即潛在停留段。

    圖2 軌跡停留指數(shù)圖樣例Fig.2 Example of a trajectory and its stop index plot

    定義3(潛在停留段):給定軌跡Traj和鄰域半徑τ,潛在停留段S=pi→pi+1→…pj是traj中的子軌跡,j≥i,且滿足:

    (1) ?k,i≥k≥j,SI(pk)≥ΔSI。

    (2) 若存在點pi-1或點pj+1,則SI(Pi-1)<ΔSI,SI(pj+1)<ΔSI。其中,ΔSI是判別潛在停留段的閾值。

    潛在停留段可以由大量高停留指數(shù)的軌跡點組成,如圖2中的③,也可以由幾個、甚至一個高停留指數(shù)的軌跡點組成,如圖2中的②。潛在停留段判別閾值ΔSI是提取潛在停留段的關(guān)鍵,若設(shè)置得太小,則會提取出較多的偽停留,而若設(shè)置得太大,則會忽略掉一些較短的停留。為此,本文給出在鄰域半徑τ下,潛在停留段判別閾值的推薦值的計算公式

    (6)

    式中,v是靜止判別閾值,為軌跡對象以速度v作勻速直線運動通過鄰域所需的最長時間。在式(6)中,靜止判別閾值v的選取需要考慮移動對象的速度,一般來說,移動對象的速度越快,v應(yīng)設(shè)置得越大,但與此同時,τ也應(yīng)相應(yīng)增大,以匹配移動對象的速度能力。考慮到人的正常步行速度在5 km/h左右[24],本文將靜止判別閾值設(shè)定為v=0.6 m/s,實際上,這也是大部分相關(guān)研究所使用的靜止判別標(biāo)準(zhǔn)[15]。

    在GPS信號良好時,一個潛在停留段即對應(yīng)一個實際停留。但是,由于室內(nèi)遮擋或高樓反射等原因,GPS的定位精度變差,產(chǎn)生漂移噪聲,一個完整的停留將對應(yīng)多個潛在停留段。如果處理不當(dāng),很可能將該完整停留誤識別成多個停留,這也是現(xiàn)有停留識別方法都難以有效解決的問題[12-14]。文獻(xiàn)[14]提出的逐級合并停留識別方法雖然在克服噪聲方面效果較好,但簡單的合并規(guī)則也導(dǎo)致了誤合并的發(fā)生[15]。圖3展示了一條有2次實際停留行為的軌跡,以及TrajDBSCAN算法和文獻(xiàn)[14]算法在各自文章的推薦參數(shù)下對停留的識別結(jié)果,其中,TrajDBSCAN法將B處的停留識別成了5個停留,而逐級合并算法則將2個停留合并成1個停留。

    圖3 不同算法的停留識別結(jié)果對比Fig.3 Comparison of stop results extracted by different algorithms

    在研究和總結(jié)現(xiàn)有方法對于停留識別的基礎(chǔ)上,本文針對源于停留指數(shù)的潛在停留段,提出了一種基于時空鄰近關(guān)系的逐級合并方法,據(jù)此發(fā)展出基于停留指數(shù)的停留識別算法。首先給出停留段概念,它或者對應(yīng)潛在停留段,或者由潛在停留段合并而來,然后定義停留段的時空相鄰關(guān)系,如下:

    定義4(時空相鄰?fù)A舳?:對于軌跡traj的兩個停留段Sa與Sb,設(shè)Sa發(fā)生在Sb之前,若Sa與Sb滿足下面兩個條件之一,那么Sa與Sb時空相鄰。

    (1) 若Sa與Sb是兩個連續(xù)的停留段,即中間不存在其他停留段,且Itv(Sa,Sb)≤MinMov。

    (2) 若Sa與Sb的凸包重疊,且Itv(Sa,Sb)≤ω·MinMov,ω>1。

    其中,Itv(Sa,Sb)表示停留段Sa與Sb之間的時間間隔;MinMov表示一個有意義的移動應(yīng)持續(xù)的最短時間。

    當(dāng)Sa與Sb時空相鄰時,可以將Sa與Sb,以及Sa與Sb之間的軌跡點,合并成一個新停留段。在定義4的2個條件中:條件1用于合并因短時脫離形成的兩個連續(xù)停留段,而條件2用于合并被噪聲數(shù)據(jù)隔開的停留段。此外,條件2的兩個輸入停留段并不要求是連續(xù)的,從而放寬了時間間隔方面的要求,本文統(tǒng)一取ω=2。

    停留的提取可以從定位潛在停留段開始,通過逐級合并時空相鄰的停留段,最終得到完整的停留。在逐級合并中,時空相鄰小停留段被合并成大停留段,大停留段則繼承小停留段的時間范圍,并將小停留段共同的凸包作為自身的凸包,以參與下一級合并。因此,停留段的逐級合并,實質(zhì)上是停留段時空范圍的擴張過程。圖4示意了4個潛在停留段的合并過程,按照時間順序分別為:S1→S2→S3→S4。

    在第1級合并中,由于S1和S3凸包重疊,使得S1、S2、S3被成功合并,其結(jié)果凸包與S4的凸包重疊,從而觸發(fā)第2級合并,最終使得4個潛在停留段被合并成1個停留段。需要指出的是,在合并停留段時,小停留段之間的軌跡點(即停留指數(shù)值較低的非停留段部分)并不參與大停留段凸包合并,這是因為非停留段部分可能包含噪聲數(shù)據(jù),若將其加入合并,勢必將導(dǎo)致大停留段的空間范圍的假性擴張,從而有可能觸發(fā)誤合并。

    綜上所述,不難得出基于停留指數(shù)的軌跡停留提取算法,其偽代碼如下:

    算法:基于停留指數(shù)的軌跡停留提取算法

    輸入:traj,τ,MinMov

    輸出:停留段集合SC

    1.for軌跡中的每一個點

    2. 根據(jù)τ計算Traj中軌跡點的停留指數(shù)

    3.endfor

    4. 據(jù)定義(3)提取所有潛在停留段,標(biāo)記為未訪問,并加入到SC

    5.while在SC中存在未被訪問的停留段

    6. 從SC中取出一個尚未被訪問的停留段S,并標(biāo)記為已訪問

    7.if在SC中存在與S時空相鄰的停留段 ∥據(jù)定義(4)

    8. 將S與其時空相鄰的所有停留段合并成大停留段S‘

    9. 將S與其時空相鄰的所有停留段的共同凸包設(shè)為S‘的凸包

    10. 將與S時空相鄰的所有停留段標(biāo)記為已訪問

    11. 將S‘加入到SC,并標(biāo)記為未訪問

    12.endif

    13.endwhile

    14. 輸出停留段集合SC

    圖4 潛在停留段合并Fig.4 Merge of potential stops

    4 軌跡停留提取試驗

    為了驗證基于停留指數(shù)的軌跡停留提取方法的有效性,本節(jié)分別選用自采和GeoLife兩組居民出行的軌跡,對日常出行中感興趣的停留開展提取試驗。試驗軌跡涵蓋了步行、公共交通、駕車等多種出行方式,其中自采軌跡為由GarmineTrex20設(shè)備采集的8條軌跡數(shù)據(jù),該設(shè)備內(nèi)置GLONASS和GPS雙星系統(tǒng),室外定位精度約為15m,軌跡點采樣間隔被固定為3s;GeoLife軌跡是微軟亞洲研究院發(fā)布的免費軌跡數(shù)據(jù)源,由于其數(shù)據(jù)采集設(shè)備規(guī)格參差不齊,筆者從中選取了15條出行信息豐富,停留數(shù)量較多、且?guī)в性肼暤能壽E,其采樣間隔為5s。

    4.1參數(shù)設(shè)置

    停留提取算法的運行共依賴于2個參數(shù),即鄰域半徑τ和移動最短持續(xù)時間MinMov,需要指出的是,停留是一個與空間尺度相關(guān)的概念(見2.2節(jié)),因而,提取算法的參數(shù)也應(yīng)與用戶所關(guān)注的空間尺度相適應(yīng)。考慮到GPS設(shè)備的定位誤差一般在15m左右,在城市范圍內(nèi)研究停留提取問題時,可依據(jù)軌跡數(shù)據(jù)的噪聲情況,將參數(shù)τ設(shè)置為30~90m之間的值,不宜過大,也不宜過小。過大時,潛在停留段的邊界將會向兩端膨脹,從而可能引發(fā)相鄰?fù)A舻腻e誤合并;過小時,潛在停留段的邊界會收縮,甚至消失,從而可能導(dǎo)致停留未被提取。參數(shù)MinMov的設(shè)置同樣與軌跡數(shù)據(jù)的噪聲情況有關(guān):當(dāng)GPS信號質(zhì)量較差時,單個停留易被分割成多潛在停留段,應(yīng)適當(dāng)增大MinMov值,以增強潛在停留段的合并能力,反之則需減小MinMov值。

    考慮到試驗數(shù)據(jù)均為居民在城市范圍內(nèi)的出行數(shù)據(jù),本文選擇τ=50m,作為停留提取的尺度參數(shù),以保證識別結(jié)果在邊界范圍上的合理性。同時,試驗數(shù)據(jù)存在較多室內(nèi)停留,致使GPS信號質(zhì)量較差,數(shù)據(jù)噪聲較大,本文采用MinMov=300s作為停留合并的閾值,以確保停留識別的完整性。此外,根據(jù)定義(2),停留應(yīng)持續(xù)一定的時間,其值設(shè)置實際上是一個應(yīng)用相關(guān)的問題。本文將其設(shè)置為180s,即表示僅對持續(xù)時間超過此值的停留感興趣,意在剔除那些因為短時減速慢行形成的,但語義蘊含較弱的停留。

    4.2結(jié)果分析

    針對自采軌跡和GeoLife軌跡的試驗結(jié)果分別見表1和表2所示,其中,真實停留數(shù)是通過人工檢驗得到的實際停留的數(shù)目,發(fā)現(xiàn)停留數(shù)是通過算法輸出的停留數(shù),正確識別停留數(shù)為輸出的停留中真實停留的數(shù)目,合并停留數(shù)為通過合并手段得到的停留數(shù),未識別停留數(shù)為真實停留中未被識別的停留數(shù)。以第13條GeoLife軌跡為例,其停留指數(shù)圖見圖5所示,相應(yīng)的停留提取結(jié)果表2第14行:人工檢出6個停留,算法輸出6個停留,其中有3個為合并而來,而在這6個輸出停留中,5個是真實停留,意味著有1個真實停留未被檢出,且多檢出1個偽停留。

    表1 自采軌跡的停留識別結(jié)果

    表2 GeoLife軌跡的停留識別結(jié)果

    圖5 GeoLife樣例軌跡的停留指數(shù)圖Fig.5 Example of GeoLife trajectory and its stop index plot

    從表1和表2可以得出,盡管試驗軌跡帶有不少含漂移噪聲的室內(nèi)停留,本文方法對于真實停留的正確檢測率仍然超過了90%。以圖5所示的第13條GeoLife軌跡為例,⑤為帶漂移噪聲的長時室內(nèi)停留,借助于時空相鄰?fù)A舳蔚闹鸺壓喜?,算法成功提取出該停留,且其完整性與實際基本相符。自采數(shù)據(jù)總共僅有6處輸出停留是通過合并得到的,原因在于自采數(shù)據(jù)的定位精度和采樣頻率較高,一些室內(nèi)停留甚至無需通過合并即可被完整的提取出來;而對于GeoLife數(shù)據(jù)來說,大量的漂移噪聲將實際停留分割成多個高停留指數(shù)的潛在停留段,有將近4成的輸出停留是通過合并得來的。

    從表1和表2也可以看出,本文方法主要存在兩方面的問題:一是誤識別,即輸出偽停留;二是未識別,即有真實停留未被檢出。從本質(zhì)上來看,這兩個問題的原因是相同的,即由于飄移噪聲的影響,屬于一個實際停留的多個停留段的時間間隔過大或者凸包不重疊,導(dǎo)致實際停留未能被正確合并出來,使得這些停留段或者被誤識別為多個小停留(即發(fā)生誤識別,如圖5中位于I區(qū)的⑥,經(jīng)過合并后仍然被分割成兩個凸包不重疊的子停留),或者因短于180 s而被去除(即導(dǎo)致未識別,如圖5中位于II區(qū)的①)。

    5 總結(jié)與展望

    本文在分析與總結(jié)現(xiàn)有的軌跡停留識別方法的基礎(chǔ)之上,提出了基于核密度思想、顧及時空鄰近性,表征軌跡點時空聚集程度的數(shù)值指標(biāo):停留指數(shù),并進(jìn)一步發(fā)展出停留指數(shù)圖。停留指數(shù)取決于軌跡在鄰域內(nèi)的逗留時間,符合人們對于停留的直觀認(rèn)識;停留指數(shù)圖通過鄰域半徑建立起同空間尺度的關(guān)系,以可視化的方式表達(dá)不同尺度下的軌跡停留情況。源于停留指數(shù),本文首先定義潛在停留段,提出了逐級合并的軌跡停留識別算法。該算法充分分析了潛在停留段之間的時空鄰近關(guān)系,能夠識別復(fù)雜多樣的軌跡停留:GPS信號良好的停留可直接由潛在停留段得來,而GPS信號微弱的停留則可通過合并多個潛在停留段而來。基于自采軌跡和GeoLife軌跡的試驗表明,本文方法可以直觀展示軌跡在其生命周期內(nèi)時空聚集程度的變化,且能更進(jìn)一步、有效提取包括單點/多點和室內(nèi)/室外在內(nèi)的多類型、含噪聲的軌跡停留。

    本文算法的不足是在處理帶有嚴(yán)重漂移噪聲的長時室內(nèi)停留時,仍有可能發(fā)生誤識別和未識別的情形,后續(xù)研究將結(jié)合地理空間上下文信息,包括路網(wǎng)匹配信息和地理逆編碼信息,以進(jìn)一步合并由噪聲隔開的停留段。

    [1]SPACCAPIETRA S, PARENT C, DAMIANI M L, et al. A Conceptual View on Trajectories[J]. Data and Knowledge Engineering, 2008, 65(1): 126-146.

    [2]YAN Zhixian, CHAKRABORTY D, PARENT C, et al. Semantic Trajectories: Mobility Data Computation and Annotation[J]. ACM Transactions on Intelligent Systems and Technology, 2013, 4(3): 49.

    [3]PARENT C, SPACCAPIETRA S, RENSO C, et al. Semantic Trajectories Modeling and Analysis[J]. ACM Computing Surveys, 2013, 45(4): 42.

    [4]向隆剛, 吳濤, 龔健雅. 面向地理空間信息的軌跡模型及時空模式查詢[J]. 測繪學(xué)報, 2014, 43(9): 982-988. DOI: 10.13485/j.cnki.11-2089.2014.0121.

    XIANG Longgang, WU Tao, GONG Jianya. A Geo-spatial Information Oriented Trajectory Model and Spatio-temporal Pattern Querying[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(9): 982-988. DOI: 10.13485/j.cnki.11-2089.2014.0121.

    [5]齊凌艷, 陳榮國, 溫馨. 基于語義軌跡停留點的位置服務(wù)匹配與應(yīng)用研究[J]. 地球信息科學(xué)學(xué)報, 2014, 16(5): 720-726.

    QI Lingyan, CHEN Rongguo, WEN Xin. Research on the LBS Matching Based on Stay Point of the Semantic Trajectory[J]. Journal of Geo-Information Science, 2014, 16(5): 720-726.

    [6]VIEIRA M, BAKALOV P, TSOTRAS V J. Querying Trajectories Using Flexible Patterns[C]∥Proceedings of the 13th International Conference on Extending Database Technology. New York: ACM, 2010: 406-417.

    [7]GIANNOTTI F, NANNI M, PEDRESCHI D, et al. Trajectory Pattern Mining[C]∥Proceedings of the 13th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. New York: ACM, 2007: 330-339.

    [8]ZHENG Yu, ZHANG Lizhu, XIE Xing, et al. Mining Interesting Locations and Travel Sequences from GPS Trajectories[C]∥Proceedings of the 18th International Conference on World Wide Web. New York: ACM, 2009: 791-800.

    [9]ALVARES L, BOGORNY V, KUIJPERS B, et al. A Model for Enriching Trajectories with Semantic Geographical Information[C]∥Proceedings of the 15th Annual ACM International Symposium on Advances in Geographic Information Systems. New York: ACM, 2007.

    [10]ROCHA J M R, TIMES V C, OLIVEIRA G, et al. DB-SMoT: A Direction-based Spatio-temporal Clustering Method[C]∥Proceedings of the 5th IEEE International Conference Intelligent Systems. London: IEEE, 2010: 114-119.

    [11]STOPHER P, JIANG Q, FITZGERALD C. Processing GPS Data from Travel Surveys[C]∥Proceedings of the 2nd International Colloquium on the Behavioral Foundations of Integrated Land-use and Transportation Models: Frameworks, Models and Applications. Toronto: [s.n.], 2005.

    [12]TIETBOHL A, BOGORNY V, KUIJPERS B, et al. A Clustering-based Approach for Discovering Interesting Places in Trajectories[C]∥Proceedings of the 2008 ACM Symposium on Applied Computing. New York: ACM, 2008: 863-868.

    [13]TRAN L H, NGUYEN Q V H, NGOC H, et al. Robust and Hierarchical Stop Discovery in Sparse and Diverse Trajectories[EB/OL]. (2012-06-15)[2015-09-24]. http:∥cn.bing.com/academic/profile?id=1894965895&v=paper_preview&mkt=zh-cn

    [14]ZIMMERMANN M, KIRSTE T, SPILIOPOULOU M. Finding Stops in Error-Prone Trajectories of Moving Objects with Time-based Clustering[M]∥TAVANGARIAN D, KIRSTE T, TIMMERMANN D, et al. Intelligent Interactive Assistance and Mobile Multimedia Computing. Berlin: Springer, 2009, 53: 275-286.

    [15]張治華. 基于GPS軌跡的出行信息提取研究[D]. 上海: 華東師范大學(xué), 2010.

    ZHANG Zhihua. Deriving Trip Information from GPS Trajectories[D]. Shanghai: East China Normal University, 2010.

    [16]Microsoft Research. GeoLife GPS Trajectories[EB/OL]. [2015-06-10]. http:∥research.microsoft.com/en-us/downloads/b16d359d-d164-469e-9fd4-daa38f2b2e13.

    [17]向隆剛, 龔健雅, 吳濤, 等. 一種面向Stop/Move抽象的軌跡時空關(guān)系[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2014, 39(8): 956-962.

    XIANG Longgang, GONG Jianya, WU Tao, et al. Spatio-temporal Trajectory Relationships Based on Stop/Move Abstraction[J]. Geomatics and Information Science of Wuhan University, 2014, 39(8): 956-962.

    [18]PARZEN E. On Estimation of a Probability Density Function and Mode[J]. The Annals of Mathematical Statistics, 1962, 33(8): 1065-1076.

    [19]王遠(yuǎn)飛, 何洪林. 空間數(shù)據(jù)分析方法[M]. 北京: 科學(xué)出版社, 2007: 57-93.

    WANG Yuanfei, HE Honglin. Spatial Data Analysis Methods[M]. Beijing: Science Press, 2007: 57-93.

    [20]禹文豪, 艾廷華. 核密度估計法支持下的網(wǎng)絡(luò)空間POI點可視化與分析[J]. 測繪學(xué)報, 2015, 44(1): 82-90. DOI: 10.11947/j.AGCS.2015.20130538.

    YU Wenhao, AI Tinghua. The Visualization and Analysis of POI Features under Network Space Supported by Kernel Density Estimation[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(1): 82-90. DOI: 10.11947/j.AGCS.2015.20130538.

    [21]SCHABENBERGER O, GOTWAY C A. Statistical Methods for Spatial Data Analysis[M]. Boca Raton: Chapman & Hall /CRC, 2005.

    [22]MILLER H J. Tobler’s First Law and Spatial Analysis[J]. Annals of the Association of American Geographers, 2004, 94(2): 284-289.

    [23]李小文, 曹春香, 常超一. 地理學(xué)第一定律與時空鄰近度的提出[J]. 自然雜志, 2007, 29(2): 69-71.

    LI Xiaowen, CAO Chunxiang, CHANG Chaoyi. The First Law of Geography and Spatial-temporal Proximity[J]. Chinese Journal of Nature, 2007, 29(2): 69-71.

    [24]CHANDRA S, BHARTI A K. Speed Distribution Curves for Pedestrians during Walking and Crossing[J]. Procedia-social and Behavioral Sciences, 2013(104): 660-667.

    (責(zé)任編輯:張艷玲)

    Visualization and Extraction of Trajectory Stops Based on Kernel-density

    XIANG Longgang1,2,SHAO Xiaotian1,2

    1. State Key Laboratory of Information Engineering in Surveying,Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China; 2. Collaborative Innovation Center of Geospatial Technology, Wuhan 430079, China

    Trajectory stops imply important semantic information, and the extraction of trajectory stops is the premise to carry out advanced Stop/Move analysis. This paper, based on the idea of kernel density, firstly introduces the concept of stop index, which is derived by cumulating spatio-temporal contribution of neighboring points, and further designs stop index graph to intuitively visualize the evolution of spatio-clustering degree during a trajectory. Importantly, stops index and its graph are related to spatial scale through neighboring radius, which then can be exploited to analyze trajectory stops under multiple scales. In addition, this paper introduces stop sequence rooted from stop index, and proposes an algorithm for the automatic extraction of trajectory stops by progressively merging stop sequences. According to the algorithm, a stop under strong GPS signal exactly corresponds to a stop sequence, while a stop under weak GPS signal could be derived by merging multiple stop sequences. Experiments based on own-acquired and GeoLife trajectories show that the algorithm has achieves the balance between the completeness and accuracy of stop extraction, and could effectively discover and extract complex and diverse trajectory stops. Even facing trajectories with serious drift noises, the algorithm still achieves a high rate of accuracy on stop extraction.

    kernel density; trajectory; stop index; stop; spatio-temporal contribution

    The National Natural Science Foundation of China (Nos.41471374;41001296)

    XIANG Longgang(1976—),male, PhD, associate professor, majors in the fields of spatial database, trajectory data analysis and virtual globe technologies.

    SHAO Xiaotian

    向隆剛,邵曉天.載體軌跡停留信息提取的核密度法及其可視化[J].測繪學(xué)報,2016,45(9):1122-1131.

    10.11947/j.AGCS.2016.20150347.

    XIANG Longgang,SHAO Xiaotian.Visualization and Extraction of Trajectory Stops Based on Kernel-density[J]. Acta Geodaetica et Cartographica Sinica,2016,45(9):1122-1131. DOI:10.11947/j.AGCS.2016.20150347.

    P208

    A

    1001-1595(2016)09-1122-10

    國家自然科學(xué)基金資助項目(41471374;41001296)

    2015-07-08

    向隆剛(1976—),男,博士,副教授,主要從事空間數(shù)據(jù)庫、軌跡數(shù)據(jù)處理分析與虛擬地球技術(shù)研究。

    E-mail: geoxlg@whu.edu.cn

    邵曉天

    E-mail: shaoxiaotian@whu.edu.cn

    修回日期: 2016-04-05

    猜你喜歡
    鄰域時空軌跡
    跨越時空的相遇
    鏡中的時空穿梭
    軌跡
    軌跡
    稀疏圖平方圖的染色數(shù)上界
    玩一次時空大“穿越”
    基于鄰域競賽的多目標(biāo)優(yōu)化算法
    軌跡
    進(jìn)化的軌跡(一)——進(jìn)化,無盡的適應(yīng)
    中國三峽(2017年2期)2017-06-09 08:15:29
    關(guān)于-型鄰域空間
    日本猛色少妇xxxxx猛交久久| 91精品伊人久久大香线蕉| 亚洲成人久久爱视频| 别揉我奶头 嗯啊视频| 免费av观看视频| 我的老师免费观看完整版| 少妇被粗大猛烈的视频| 日本爱情动作片www.在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲天堂国产精品一区在线| 国内精品一区二区在线观看| 深夜a级毛片| 有码 亚洲区| 日韩亚洲欧美综合| 国产一区有黄有色的免费视频 | 中文资源天堂在线| 国产精品福利在线免费观看| 一个人免费在线观看电影| 亚洲精品成人久久久久久| 国产成人精品婷婷| 男插女下体视频免费在线播放| 色尼玛亚洲综合影院| 久久鲁丝午夜福利片| 日本与韩国留学比较| 国产高清国产精品国产三级 | 亚洲自拍偷在线| 国产视频首页在线观看| 欧美3d第一页| 日本熟妇午夜| 国产精品av视频在线免费观看| 免费大片18禁| 中文字幕亚洲精品专区| 黄色一级大片看看| 久久精品国产自在天天线| 男女国产视频网站| 午夜激情欧美在线| 赤兔流量卡办理| 国产爱豆传媒在线观看| 乱人视频在线观看| 99热这里只有是精品在线观看| 夜夜爽夜夜爽视频| 午夜福利网站1000一区二区三区| 深夜a级毛片| 国产爱豆传媒在线观看| 亚洲精品久久久久久婷婷小说 | 天堂√8在线中文| 有码 亚洲区| 亚洲av一区综合| 18禁裸乳无遮挡免费网站照片| 亚洲国产色片| 少妇人妻精品综合一区二区| 麻豆一二三区av精品| 午夜福利高清视频| 99九九线精品视频在线观看视频| 嘟嘟电影网在线观看| 欧美激情国产日韩精品一区| 久久精品国产自在天天线| 亚洲欧美日韩卡通动漫| 老师上课跳d突然被开到最大视频| 老司机影院毛片| 国产精品av视频在线免费观看| 天天躁日日操中文字幕| av黄色大香蕉| 身体一侧抽搐| 国产精品爽爽va在线观看网站| 少妇人妻精品综合一区二区| 国产成人freesex在线| 亚洲欧美日韩无卡精品| 少妇裸体淫交视频免费看高清| 国产69精品久久久久777片| 色5月婷婷丁香| 亚洲精品,欧美精品| 久久久久久国产a免费观看| 久久国内精品自在自线图片| 啦啦啦韩国在线观看视频| 欧美一区二区精品小视频在线| 九色成人免费人妻av| 在线免费十八禁| 国产高清三级在线| 久久精品91蜜桃| 国产亚洲91精品色在线| 热99re8久久精品国产| 国产成人午夜福利电影在线观看| 亚洲欧美日韩无卡精品| av卡一久久| 色尼玛亚洲综合影院| 国产av码专区亚洲av| 久久婷婷人人爽人人干人人爱| 午夜福利在线观看吧| 日韩欧美国产在线观看| 欧美精品一区二区大全| 最近视频中文字幕2019在线8| 免费黄色在线免费观看| 成人高潮视频无遮挡免费网站| 日韩中字成人| 亚洲五月天丁香| 国产在线一区二区三区精 | 一个人免费在线观看电影| 老司机福利观看| 99久久九九国产精品国产免费| 久久久成人免费电影| 日韩欧美精品免费久久| 国产 一区 欧美 日韩| 亚洲第一区二区三区不卡| 天堂√8在线中文| 国产麻豆成人av免费视频| 秋霞伦理黄片| 亚洲精品成人久久久久久| 少妇猛男粗大的猛烈进出视频 | 亚洲av福利一区| 在线免费观看的www视频| 国产av不卡久久| 午夜亚洲福利在线播放| 国产黄片视频在线免费观看| 久久午夜福利片| 精品无人区乱码1区二区| 麻豆av噜噜一区二区三区| 国产综合懂色| www.av在线官网国产| 久久婷婷人人爽人人干人人爱| 男插女下体视频免费在线播放| 联通29元200g的流量卡| 精品人妻视频免费看| av视频在线观看入口| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看 | 欧美+日韩+精品| 日本午夜av视频| 18禁在线无遮挡免费观看视频| 观看免费一级毛片| 亚洲真实伦在线观看| 日日干狠狠操夜夜爽| av黄色大香蕉| 只有这里有精品99| 国语对白做爰xxxⅹ性视频网站| 一级爰片在线观看| 村上凉子中文字幕在线| 亚洲va在线va天堂va国产| 久久久a久久爽久久v久久| 成人毛片a级毛片在线播放| 久久久久久久午夜电影| 亚洲不卡免费看| 国产午夜福利久久久久久| 成年女人永久免费观看视频| 国产精品久久久久久久久免| 午夜精品国产一区二区电影 | 国产一级毛片在线| 黄色日韩在线| or卡值多少钱| 国产综合懂色| 亚洲精品,欧美精品| 亚洲国产精品成人综合色| 老师上课跳d突然被开到最大视频| 国产精品综合久久久久久久免费| 日日摸夜夜添夜夜爱| 久久人妻av系列| 日韩人妻高清精品专区| 亚洲五月天丁香| 青春草视频在线免费观看| 2022亚洲国产成人精品| 亚洲国产最新在线播放| 亚洲无线观看免费| 久久久久久久久久黄片| av女优亚洲男人天堂| 97超碰精品成人国产| 26uuu在线亚洲综合色| 亚洲va在线va天堂va国产| 午夜视频国产福利| 美女国产视频在线观看| 色视频www国产| av福利片在线观看| 国产精品精品国产色婷婷| 99热这里只有精品一区| 婷婷色麻豆天堂久久 | 亚洲欧美精品自产自拍| 我要搜黄色片| 国产探花在线观看一区二区| 国产 一区精品| 亚洲综合精品二区| 国产亚洲精品av在线| 久久婷婷人人爽人人干人人爱| 简卡轻食公司| 亚洲国产精品合色在线| 亚洲国产欧洲综合997久久,| 老司机福利观看| 大话2 男鬼变身卡| av在线播放精品| 成年av动漫网址| 两个人的视频大全免费| 日韩av在线大香蕉| 日日摸夜夜添夜夜添av毛片| 亚洲一级一片aⅴ在线观看| 国产中年淑女户外野战色| 黄色一级大片看看| 我要看日韩黄色一级片| 日韩欧美国产在线观看| 亚洲精品日韩av片在线观看| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放 | 少妇的逼水好多| eeuss影院久久| 亚洲国产色片| 欧美激情久久久久久爽电影| 国产精品爽爽va在线观看网站| 日本三级黄在线观看| 国产爱豆传媒在线观看| 久久久久九九精品影院| 九草在线视频观看| 国产精品精品国产色婷婷| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品亚洲一区二区| 毛片女人毛片| 女人久久www免费人成看片 | 在线a可以看的网站| a级毛片免费高清观看在线播放| 中文字幕制服av| 纵有疾风起免费观看全集完整版 | 日本一二三区视频观看| 国产毛片a区久久久久| 少妇的逼好多水| 久久久久久久国产电影| 国产三级中文精品| 男人和女人高潮做爰伦理| 国产免费男女视频| 久久久久久久午夜电影| 精品久久久噜噜| 97人妻精品一区二区三区麻豆| 国产精品熟女久久久久浪| 国产精品乱码一区二三区的特点| 中文字幕熟女人妻在线| 精品国内亚洲2022精品成人| 色吧在线观看| 精品99又大又爽又粗少妇毛片| 国内少妇人妻偷人精品xxx网站| 99久久无色码亚洲精品果冻| 亚洲在久久综合| 天堂√8在线中文| 精品久久久久久成人av| 深夜a级毛片| 久久久国产成人免费| 大香蕉久久网| 99九九线精品视频在线观看视频| 国产久久久一区二区三区| 女人久久www免费人成看片 | a级毛色黄片| 欧美激情国产日韩精品一区| 国产 一区精品| 嫩草影院精品99| 国产真实乱freesex| 亚洲图色成人| 99久久九九国产精品国产免费| 看非洲黑人一级黄片| 丝袜美腿在线中文| 亚洲图色成人| 精品国产露脸久久av麻豆 | 我的老师免费观看完整版| 国产一区二区亚洲精品在线观看| 久久久色成人| 久久精品熟女亚洲av麻豆精品 | 欧美激情久久久久久爽电影| 免费av毛片视频| 男人和女人高潮做爰伦理| 久久久久国产网址| 偷拍熟女少妇极品色| 联通29元200g的流量卡| 99热网站在线观看| 精品久久久久久成人av| 婷婷色综合大香蕉| 欧美3d第一页| 人人妻人人澡欧美一区二区| 午夜激情福利司机影院| 国内精品美女久久久久久| 欧美成人a在线观看| 国产亚洲精品av在线| 国产亚洲91精品色在线| 国产综合懂色| 免费av毛片视频| 91精品国产九色| 男的添女的下面高潮视频| 丰满少妇做爰视频| 亚洲经典国产精华液单| 一夜夜www| 看黄色毛片网站| 美女黄网站色视频| videos熟女内射| 黄片wwwwww| 日韩一区二区三区影片| 国产三级中文精品| 午夜激情欧美在线| 国产熟女欧美一区二区| 麻豆国产97在线/欧美| 亚洲欧美精品专区久久| 国产成人精品婷婷| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 五月玫瑰六月丁香| 国产av码专区亚洲av| 麻豆国产97在线/欧美| 99热这里只有是精品在线观看| 一个人看的www免费观看视频| 亚洲国产欧洲综合997久久,| 日韩一区二区三区影片| 看十八女毛片水多多多| 国产精品,欧美在线| 老司机影院成人| 精品久久久噜噜| 午夜日本视频在线| 综合色丁香网| 亚洲人成网站在线播| 99视频精品全部免费 在线| 少妇的逼好多水| 欧美丝袜亚洲另类| 色5月婷婷丁香| 日本wwww免费看| 亚洲av免费高清在线观看| 秋霞伦理黄片| www日本黄色视频网| 能在线免费看毛片的网站| 成人综合一区亚洲| 国产三级中文精品| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 少妇被粗大猛烈的视频| 中文字幕亚洲精品专区| 久久精品久久久久久噜噜老黄 | 汤姆久久久久久久影院中文字幕 | 国产一级毛片在线| 美女大奶头视频| 久久精品国产亚洲网站| 亚洲激情五月婷婷啪啪| 国产精品久久视频播放| 亚洲成av人片在线播放无| 国产成人福利小说| 观看免费一级毛片| 少妇人妻一区二区三区视频| 一二三四中文在线观看免费高清| 国产一区有黄有色的免费视频 | 久久久久久国产a免费观看| 国产成人a区在线观看| 国产私拍福利视频在线观看| 最近最新中文字幕大全电影3| 久久婷婷人人爽人人干人人爱| 欧美极品一区二区三区四区| 久久这里只有精品中国| 高清在线视频一区二区三区 | 免费看日本二区| 2021天堂中文幕一二区在线观| 少妇熟女欧美另类| 男人舔奶头视频| 欧美xxxx性猛交bbbb| 国产高潮美女av| 精品一区二区三区视频在线| 日本免费一区二区三区高清不卡| 免费黄色在线免费观看| 2021天堂中文幕一二区在线观| 99热全是精品| 亚洲欧美精品专区久久| 91精品国产九色| 久久国产乱子免费精品| 久久久久久久国产电影| 久久99热这里只有精品18| 亚洲美女搞黄在线观看| www.色视频.com| 国产av一区在线观看免费| 久久这里只有精品中国| 亚洲天堂国产精品一区在线| 国产欧美另类精品又又久久亚洲欧美| 免费观看精品视频网站| 成人特级av手机在线观看| 午夜视频国产福利| 18禁在线播放成人免费| 日韩,欧美,国产一区二区三区 | 亚洲欧美精品自产自拍| 亚洲欧美中文字幕日韩二区| 国产三级中文精品| 美女高潮的动态| 国产精品无大码| 成人二区视频| 少妇被粗大猛烈的视频| 国产av一区在线观看免费| 亚洲国产精品专区欧美| 最近最新中文字幕大全电影3| 美女黄网站色视频| 一级毛片aaaaaa免费看小| 变态另类丝袜制服| 亚洲成av人片在线播放无| 听说在线观看完整版免费高清| 久久亚洲国产成人精品v| 听说在线观看完整版免费高清| 中国美白少妇内射xxxbb| 国产国拍精品亚洲av在线观看| 亚洲av电影不卡..在线观看| 国产亚洲5aaaaa淫片| 丝袜美腿在线中文| 免费一级毛片在线播放高清视频| 久久精品夜夜夜夜夜久久蜜豆| 国产精华一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 97超视频在线观看视频| 美女高潮的动态| 高清午夜精品一区二区三区| 乱人视频在线观看| 老师上课跳d突然被开到最大视频| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久com| 欧美+日韩+精品| 亚洲欧洲日产国产| 99久久成人亚洲精品观看| 亚洲成人精品中文字幕电影| 一本一本综合久久| 国产成人freesex在线| 国产黄a三级三级三级人| 国产av一区在线观看免费| 国产一区二区三区av在线| 简卡轻食公司| 夫妻性生交免费视频一级片| 国产黄片视频在线免费观看| 国产又色又爽无遮挡免| 欧美一区二区亚洲| 亚洲色图av天堂| 丰满少妇做爰视频| 亚洲精品aⅴ在线观看| 狂野欧美白嫩少妇大欣赏| 汤姆久久久久久久影院中文字幕 | 99久久人妻综合| 99久久中文字幕三级久久日本| 成人美女网站在线观看视频| 国产亚洲最大av| 国产免费男女视频| videossex国产| 精品熟女少妇av免费看| 午夜激情福利司机影院| www日本黄色视频网| 直男gayav资源| 亚洲精品色激情综合| 免费一级毛片在线播放高清视频| 综合色丁香网| 国产高清不卡午夜福利| 国产精品久久久久久精品电影小说 | 2021天堂中文幕一二区在线观| 蜜桃亚洲精品一区二区三区| 亚洲国产精品专区欧美| 一级毛片aaaaaa免费看小| 我要看日韩黄色一级片| 色综合色国产| 亚洲av免费高清在线观看| 又粗又爽又猛毛片免费看| 精品无人区乱码1区二区| 97人妻精品一区二区三区麻豆| 色尼玛亚洲综合影院| 亚洲国产欧美在线一区| 97超碰精品成人国产| 狠狠狠狠99中文字幕| 在线观看66精品国产| 国产精品熟女久久久久浪| 亚洲精品日韩在线中文字幕| 18禁动态无遮挡网站| 精品一区二区免费观看| 精品人妻熟女av久视频| 亚洲乱码一区二区免费版| 内地一区二区视频在线| 亚洲18禁久久av| 免费av观看视频| 亚洲av熟女| 高清午夜精品一区二区三区| 少妇人妻精品综合一区二区| 少妇丰满av| 国产亚洲av片在线观看秒播厂 | 免费观看性生交大片5| 亚洲成人av在线免费| 丰满少妇做爰视频| 男的添女的下面高潮视频| 免费看av在线观看网站| 九九久久精品国产亚洲av麻豆| 国产精品1区2区在线观看.| 一个人观看的视频www高清免费观看| 熟女电影av网| 久久精品夜色国产| 久久精品国产自在天天线| videos熟女内射| 亚洲国产欧洲综合997久久,| 人人妻人人看人人澡| 日日摸夜夜添夜夜爱| 成年女人看的毛片在线观看| 国产视频内射| 久久久成人免费电影| 成人欧美大片| 久久6这里有精品| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲5aaaaa淫片| 精品99又大又爽又粗少妇毛片| 国产精品人妻久久久影院| 亚洲成av人片在线播放无| 国产成人a∨麻豆精品| 亚洲av日韩在线播放| 久久精品久久久久久久性| 国产精品久久久久久精品电影| 国产精品野战在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美精品自产自拍| 欧美丝袜亚洲另类| 色综合亚洲欧美另类图片| 国产精品野战在线观看| 日韩av在线大香蕉| 男的添女的下面高潮视频| 99久国产av精品| 亚洲精品色激情综合| 尾随美女入室| 一边摸一边抽搐一进一小说| 婷婷色麻豆天堂久久 | 麻豆国产97在线/欧美| 国产精品乱码一区二三区的特点| 亚洲av免费高清在线观看| 国产午夜精品久久久久久一区二区三区| 尤物成人国产欧美一区二区三区| a级毛片免费高清观看在线播放| 久久久久国产网址| 麻豆成人午夜福利视频| 婷婷色综合大香蕉| 大又大粗又爽又黄少妇毛片口| 天美传媒精品一区二区| 国产成人精品久久久久久| a级毛色黄片| 国产成人免费观看mmmm| 中国美白少妇内射xxxbb| 丰满少妇做爰视频| 免费看a级黄色片| 中文在线观看免费www的网站| 国产视频首页在线观看| 在线a可以看的网站| 国产中年淑女户外野战色| 婷婷六月久久综合丁香| 人人妻人人看人人澡| 国产精品久久视频播放| 日日啪夜夜撸| 尤物成人国产欧美一区二区三区| 中文字幕熟女人妻在线| 精品国产露脸久久av麻豆 | 99久久人妻综合| 欧美性猛交黑人性爽| 老司机影院毛片| 久久人人爽人人片av| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲91精品色在线| 成人高潮视频无遮挡免费网站| 99在线视频只有这里精品首页| 免费观看a级毛片全部| 欧美97在线视频| 久久6这里有精品| 男的添女的下面高潮视频| 免费黄网站久久成人精品| 久久精品综合一区二区三区| av福利片在线观看| 麻豆精品久久久久久蜜桃| 日韩大片免费观看网站 | 国产一区二区在线av高清观看| 欧美不卡视频在线免费观看| 欧美xxxx性猛交bbbb| 精品一区二区免费观看| 中文资源天堂在线| 免费观看的影片在线观看| 亚洲精品日韩av片在线观看| 99热精品在线国产| 性插视频无遮挡在线免费观看| 日日啪夜夜撸| 国产高清三级在线| 91在线精品国自产拍蜜月| 我的老师免费观看完整版| 精品久久久久久成人av| 国产一区二区亚洲精品在线观看| 麻豆成人av视频| 国产免费视频播放在线视频 | 人妻夜夜爽99麻豆av| 久久人人爽人人片av| 亚洲成人中文字幕在线播放| 免费黄网站久久成人精品| 亚洲av男天堂| 久久久久久久亚洲中文字幕| 国产成人福利小说| 欧美日本视频| av在线播放精品| 最近中文字幕2019免费版| 国产亚洲精品久久久com| 日本免费一区二区三区高清不卡| 国产亚洲一区二区精品| 亚洲成人中文字幕在线播放| 高清在线视频一区二区三区 | 午夜激情福利司机影院| 亚洲乱码一区二区免费版| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 亚洲人成网站在线观看播放| 青青草视频在线视频观看| 亚洲成色77777| 国产老妇伦熟女老妇高清| 中文欧美无线码| 内地一区二区视频在线| 亚洲自偷自拍三级| 免费av观看视频| 丝袜美腿在线中文| 国产精品久久久久久久电影| 少妇熟女欧美另类| 亚洲自偷自拍三级| 国产白丝娇喘喷水9色精品| 精品无人区乱码1区二区| 国产成年人精品一区二区| 少妇高潮的动态图| 麻豆乱淫一区二区| 亚洲熟妇中文字幕五十中出| 黄色日韩在线| 国产单亲对白刺激| 欧美一区二区国产精品久久精品| 欧美日韩一区二区视频在线观看视频在线 | 久久99热这里只有精品18| 国产精品99久久久久久久久|