PC-TC/formatPOIsWithStayPoints.m at master · YunlongWang/PC-TC · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
%
% Copyright (c) 2018, Yunlong Wang
% All rights reserved. Please read the "license.txt" for license terms.
%
% This code is part of PC-TC(v0.1)
%
% Contact Info: [email protected]
%

function POIs = formatPOIsWithStayPoints(data_original,stayPoints,data_with_clusters,IDX)

%% format POIs
maxNum = max(IDX);
latitudeEachCluster = zeros(maxNum,1);
longitudeEachCluster = zeros(maxNum,1);
radiusEachCluster = zeros(maxNum,1);
hitEach=zeros(maxNum,1);
durationEach=zeros(maxNum,1);
hitDays=zeros(maxNum,1);


for i=1:maxNum
   latitudeEachCluster(i) = mean(stayPoints.lat(IDX==i));
   longitudeEachCluster(i) = mean(stayPoints.lng(IDX==i));
   maxLat = max(stayPoints.lat(IDX==i));
   minLat = min(stayPoints.lat(IDX==i));
   if maxLat-latitudeEachCluster(i)>latitudeEachCluster(i)-minLat
      radiusLat = maxLat-latitudeEachCluster(i);
   else
      radiusLat = latitudeEachCluster(i)-minLat;
   end
   maxLng = max(stayPoints.lng(IDX==i));
   minLng = min(stayPoints.lng(IDX==i));
   if maxLng-longitudeEachCluster(i)>longitudeEachCluster(i)-minLng
      radiusLng = maxLng-longitudeEachCluster(i);
   else
      radiusLng = longitudeEachCluster(i)-minLng;
   end

   if radiusLat >radiusLng
       distInDegree = radiusLat;
   else
       distInDegree = radiusLng;
   end

   radiusEachCluster(i) = 40075017 * distance(latitudeEachCluster(i),longitudeEachCluster(i),latitudeEachCluster(i)+distInDegree,longitudeEachCluster(i))/360;

   tempIndex = find(data_with_clusters.clusters_final==i);
   hitEach(i)=size(tempIndex,1);
   durationEach(i)=sum(data_with_clusters.timeSpan(tempIndex));
   dayLabels = data_with_clusters.day_labels(tempIndex);
   hitDays(i) = size(unique(dayLabels),1);
end

durationPerHitDay = durationEach./double(hitDays);

POIs_center_latitude = latitudeEachCluster;
POIs_center_longitude = longitudeEachCluster;
POIs_radius = radiusEachCluster;
POIs_visit_times = hitEach;
POIs_coverage_days = ones(maxNum,1)*max(data_original.day_labels);% not D!!! because D may not contain all the day labels
POIs_average_duration = durationEach./POIs_visit_times;

popularPlacesClusterIndex = (1:maxNum)';
clusters = table(popularPlacesClusterIndex,POIs_center_latitude,POIs_center_longitude,POIs_radius,POIs_visit_times,POIs_coverage_days,POIs_average_duration);

POIs = clusters;