查找軌跡的全部駐留點。駐留點是移動對象長期停留的點。

語法

SETOF trstaypoint ST_StayPoint(trajectory traj, interval resample_int, float dr, interval dt, float thres);
TABLE(point geometry, length interval) ST_stayPointSpatialClustered(trajectory traj, interval resample_int, float dr, interval dt, float thres, float merger)

參數

參數名稱 描述
traj 軌跡對象。
resample_int 軌跡重采樣的時間間隔。重采樣可以讓軌跡的采樣點更加均勻。時間間隔越小，駐留點的計算越精確，但計算更慢。當不需要重采樣時，可以將此值設置為很大，例如'100000 day'。更多信息，請參見ST_Resample。
dr 空間容差。按照軌跡的SRID所對應的單位進行計算。例如，對SRID=0的軌跡，按照歐式距離計算，對SRID=4326的軌跡，按照米計算。
dt 時間容差。應該略長于單次駐留的時間。
thres 當一個點周邊（空間距離小于dr，時間差小于dt）的軌跡采樣點數量超過thres時，認為其是駐點。當連續多個點都被識別為駐點時，會找到周邊點數最多的點作為最終的駐點。
merger 將距離小于merger的兩個不同時間的駐點合并為一個。
返回值類型trstaypoint是包含3列的復合類型，三列分別如下：
列名 列類型 說明
point geometry 駐點的位置。
startt timestamp 開始駐留的時間（按第一個被識別為駐點的點計算）。
endt timestamp 結束駐留的時間（按最后一個被識別為駐點的點計算）。

描述

本函數識別一條軌跡的駐點。例如，如果我們認為在500米范圍內停留30分鐘以上算作駐留，則可以選取resample_int為5分鐘，dr為250米，dt為40分鐘，而thres為7（30分鐘處以5分鐘，再加1）進行識別。

對于有多次往返線路的軌跡，可以使用ST_stayPointSpatialClustered函數，將空間上相近的駐點進行聚合，并得到總的駐點和停留時間。這些駐點可能對應著移動對象經常停留的地點，如加油站、停車場等。

示例

SELECT ST_AsText((ST_StayPoint(traj, '12 hour', 3, '3 day', 3)).point),(ST_StayPoint(traj, '12 hour', 3, '3 day', 3)).*
FROM (SELECT ST_MakeTrajectory('LINESTRING(1 1, 2 2, 3 3, 10 10, 11 11, 12 12, 13 13, 3 3, 2 2, 2 2, 3 3, 8 8,10 10,9 9,8 8)') AS traj) tj;
   st_astext    |                   point                    |       startt        |        endt
----------------+--------------------------------------------+---------------------+---------------------
 POINT(1 1)     | 0101000000000000000000F03F000000000000F03F | 2000-01-01 00:00:00 | 2000-01-03 00:00:00
 POINT(11 11)   | 010100000000000000000026400000000000002640 | 2000-01-04 00:00:00 | 2000-01-07 00:00:00
 POINT(3 3)     | 010100000000000000000008400000000000000840 | 2000-01-08 00:00:00 | 2000-01-11 00:00:00
 POINT(9.5 9.5) | 010100000000000000000023400000000000002340 | 2000-01-12 00:00:00 | 2000-01-15 00:00:00
(4 rows)

SELECT ST_AsText((ST_stayPointSpatialClustered(traj, '12 hour', 3, '3 day', 3, 3)).point), (ST_stayPointSpatialClustered(traj, '12 hour', 3, '3 day', 3, 3)).*
FROM (SELECT ST_MakeTrajectory('LINESTRING(1 1, 2 2, 3 3, 10 10, 11 11, 12 12, 13 13, 3 3, 2 2, 2 2, 3 3, 8 8,10 10,9 9,8 8)') AS traj) tj;
     st_astext      |                   point                    | length
--------------------+--------------------------------------------+--------
 POINT(10.25 10.25) | 010100000000000000008024400000000000802440 | 6 days
 POINT(2 2)         | 010100000000000000000000400000000000000040 | 5 days
(2 rows)