Next Article in Journal
Fast Registration of Point Cloud Based on Custom Semantic Extraction
Previous Article in Journal
Research on Ground Object Classification Method of High Resolution Remote-Sensing Images Based on Improved DeeplabV3+
Previous Article in Special Issue
A New Wavelet-Based Privatization Mechanism for Probability Distributions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pedestrian Flow Prediction and Route Recommendation with Business Events †

1
School of Computer Science and Engineering, University of Electronic Science and Technology of China, Qingshuihe Campus, Chengdu 611731, China
2
College of Computer Science and Electronic Engineering, Hunan University, Changsha 410082, China
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in Gu, J.; Song, C.; Jiang, W.; Shi, L.; Gong, H.; Liu, M. Pedestrian Flow Prediction with Business Events. In proceedings of the 2019 IEEE International Conference on Mobile Ad-Hoc and Sensor Networks (MSN), Shenzhen, China, 11–13 December 2019; pp. 43–48. https://doi.org/10.1109/MSN48538.2019.00022.
Sensors 2022, 22(19), 7478; https://doi.org/10.3390/s22197478
Submission received: 5 September 2022 / Revised: 26 September 2022 / Accepted: 29 September 2022 / Published: 2 October 2022
(This article belongs to the Special Issue Internet of Things, Big Data and Smart Systems)

Abstract

:
Due to the potential economic benefits, pedestrian flow is considered an essential indication of public spaces. Pedestrian flow prediction is designed to assist operators in making decisions (such as shopping center owners). Operators hold certain events, such as sales promotions, to attract surrounding pedestrians; we refer to this type of event as a business event. Business events attract pedestrian flows, which means business opportunities for the merchants. Moreover, their placement will affect the distributions of the pedestrian flows. However, deciding which route is chosen for a specified event is difficult. To the best of our knowledge, we are the first to consider business events when predicting pedestrian flow. In this paper, we investigate two problems: one is pedestrian flow prediction with business events, and the other is route recommendation for business events. First, we propose an Attraction-Based Matrix Factorization model (ABMF) to efficiently predict the pedestrian flow with business events, which introduces the attraction index of different categories to pedestrians in matrix factorization. Second, we leverage the Skip-gram mode to learn the latent representations and improve the pair-wise ranking loss to a flow-aware-based method (SG-FWARP), which aims to learn events’ latent representations for route recommendation. Compared with other state-of-the-art methods, the experimental results show ABMF can predict pedestrian flow matrix with a similarity of over 0.9 compared with the ground truth, and SG-FWARP can recommend routes for business events with high accuracy.

1. Introduction

The rapid progress of urbanization brings prosperity to businesses. A large number of human beings visit the business district, which forms dynamic pedestrian flows between Points of interests (POIs). Mobile crowdsourcing data from location-based social network services (LBSNs) provide information on individuals’ preferences for routes and locations [1,2,3]. Pedestrian flow prediction provides crucial information to operators for decision-making. For instance, by utilizing crowdsourced trajectories, pedestrian flow analysis enhances the structure of technological show events [4], and pedestrian flow prediction in vast road networks is used to benefit the German outdoor advertising market [5]. The pedestrian flow provides the operators purchasing power and has the potential ability to benefit them. Additionally, the operators frequently organize events to draw the nearby pedestrian flow in order to gain greater commercial benefits. Researchers studied the prediction with different types of events in real life. For example, the work [6] predicts the crowdedness of POIs in order to enhance personalized trip recommendations by avoiding crowds. The work [7] forecasts attendance of the activity in event-based social networks (e.g., Douban).
However, events that occur along the route are rarely studied by academics, such as sales promotions, and we refer to this novel type of event as business events. Unlike previous studies on pedestrian flows, we focus on pedestrian flow prediction under business events. Business events are frequently hosted in areas with high pedestrian flows because event organizers want to draw attention to their products by gathering pedestrian flow. Figure 1 depicts the effect of a business event on pedestrian flow: there are two paths from v 1 to v 4 , both of which pass through v 2 and v 3 . Two business events are associated with the two routes v 1 v 3 and v 3 v 4 , respectively. These events attract a greater number of pedestrians. As business events are scheduled to be placed, predicting pedestrian flow with business events is a big challenge. The challenge is that different business events on different routes will result in varied distributions of the overall pedestrian flow.
In this paper, by utilizing mobile crowdsourcing data, we study how business events affect pedestrian flows. Based on historical pedestrian flows, we present an Attraction- Based Matrix Factorization model (ABMF) to predict pedestrian flows with business events accurately, which introduces the attraction index of different categories to pedestrians in matrix factorization. In some cases, the operator only knows which event he plans to hold but does not know which route fits the event. Thus, it is essential to recommend routes for a business event. So, we investigate which routes fit for holding an event, and recommend top-N routes for a specified business event. The goal of top-N route recommendation is to recommend for each event N routes that are most appropriate for it. To achieve this goal, by leveraging the Skip-gram model [8], we learn the latent representations for a route to capture the contextual information in sequences. By considering the flow factors, we improve the pair-wise ranking loss to a flow-aware-based method, which aims to learn events’ latent representations for route recommendation. The main contributions of this paper are:
  • As far as we know, we are the first to take business events into account when predicting pedestrian flow.
  • To predict pedestrian flow based on business events, we present the ABMF model, and its highlight lies in that it introduces the attraction index of different categories to pedestrians in matrix factorization.
  • To recommend routes for event placement, we learn route representations based on the Skip-gram model and consider the flow factors to recommend top-N routes for events.
  • We compare the performance of the proposed two methods with state-of-the-art solutions on a simulation dataset and real-world datasets, and experiment results reveal that our algorithms outperform the baselines.
This paper is organized as follows: Section 2 reviews related work. Section 3 shows the data model and problem definition. We detail the ABMF method and the parameter estimation in Section 4. Then, we introduce the extension of route recommendation for events in Section 5. We report the experimental results in Section 6 and conclude our paper in Section 7.

2. Related Work

In this section, we review the related work from the following aspects:
(1) Pedestrian Flow Analysis: Pedestrian flow analysis, which is used to extract information from stored trajectories via pedestrian flow modeling, simulation, and optimization, has garnered increased attention. In the literature on pedestrian modeling, a number of modeling techniques have already been proposed at microscopic and macroscopic scales, including trajectory clustering for extracting similar trajectories from a dataset [9]. An analysis of pedestrian flow in complicated scenes computes a representation of the main pedestrian flows [10]. The work [11] predicts the influence of urban environment layout on the spatial distribution of pedestrian flows. Based on the pedestrian trajectories, the work [4] shows a technique for real pedestrian track analysis during an actual exhibition. There are numerous approaches for simulating pedestrian traffic. The work [12] simulates the pedestrian flow by considering the interaction among pedestrians. Kaminka et al. [13] discuss a model of urban pedestrian flow based on agents.
(2) Pedestrian Flow Prediction: Predicting pedestrian flow is a fundamental topic in urban computing [14]. Previous research focuses on predicting people’s movements using their past location data [15,16]. In recent years, researchers have developed methods to forecast the city’s crowd flow in a variety of situations, including taxi and bike movements [17]. Zhou et al. [18] propose MOHER to predict the potential crowd flow in a certain mode, which uses the LTSM module to predict the sequential flow. To predict the crowd flow over the entire city, Zhang et al. [19] develop an end-to-end ST-ResNet structure based on unique spatio-temporal data features. Some studies concentrate on predicting pedestrian flow: Eravci et al. [20] forecast the scale of pedestrian flows to provide commercial advice; Ma et al. [21] propose an approach that automatically predicts crowd density in the short term, which presents a prediction algorithm using v-support vector regression (vSVR). Duan et al. [22] introduce a complementary attention gated network to predict pedestrian trajectory, which captures both frequent and peculiar modals in spatial and temporal patterns.
(3) Urban Computing and Business Recommendation: The rich human mobility data generated in urban spaces reflects a city’s underlying problems, helping urban planners to formulate plans effectively. Big data-driven urban computing is a popular field with many valuable applications. Li et al. [23] select locations for ambulance stations. Bao et al. [24] plan bike lanes relaying on the sharing-bikes’ trajectories. Liu et al. [25] aim to select billboard locations based on the large-scale taxi trajectories. The recommendation technique is employed in business category selection, store location selection, etc. Zhao et al. [26] utilize the data from location-based social networks to recommend new business categories in a partitioned business district, which mines the business opportunities and guides the planners to open new commercial shops in certain categories in a specific district. Some researchers employ the business data to support the business owners in LBSNs, designing the zone recommendation system [27], business prediction system [28], and retail allocation system [29]. Lu et al. [30] apply user and business properties for personalized business recommendation.
(4) Sequential Modeling and Embedding Learning: Modeling the sequential pattern is an important technique in the recommendation field. Many works apply the Markov chain in the successive check-ins to observe the sequential pattern. The studies in [31] learn the transitive pattern of categories in sequential check-ins. Zhang et al. [32] learn the transitive probability from the additive Markov chain to recommend POIs. With the success of deep learning, the neural network has been utilized to model the check-in sequences. The work [33] applies the word2vec framework to model the check-in sequences to capture the sequential contexts. The prior works [31,33] inspire us to capture the sequential pattern from user travel trajectories to improve business event recommendation. The word2vec framework [8] is a neural language model to learn latent representations of words effectively. The main idea is to observe a word’s contextual relations in sentences, which perform better than the viewpoint of word similarity and transitivity in sentences. With the success of the framework in capturing the contextual correlations of items, the word2vec framework is widely used for user modeling [34] and POI modeling [35], etc. For example, Liu et al. [33] adopt a Skip-Gram model in POI recommendation by capturing the sequential POI check-ins and performs better than the Markov chain model. However, the work [33,36] neglects the crowd flow influence, so we propose a flow-aware model to recommend routes for the business event. Compared to these studies, the main difference lies in our proposed methods naturally focusing on predicting the pedestrian flow with business events and recommending routes for business events. This prediction assists operators in making decisions on where to place business events, increasing possible economic gains.

3. Preliminaries

The motivation to predict pedestrian flow with business events is as follows: The business events that are held along the route will attract pedestrian flow, which could have advantages and is the driving force behind the prediction of pedestrian flow with business events. This section explains the data model and defines the problem of predicting pedestrian flow under business activities. The most common notations used in this paper are listed in Table 1.

3.1. Data Model

Assume there are n locations (or POIs) and k business events in a region, which is represented by a fully connected graph G = ( V , E ) , where V represents the set of locations and E represents the set of routes. We assume that, at most, one business event will occur on the route between two locations. We treat several business events that take place between two locations as a single event with a mixed preference, which covers the preference of every business event on this route. We describe a pedestrian flow matrix as follows:
Definition 1.
A pedestrian flow matrix is a n-order matrix denoted as M, with the element m i j referring to the pedestrian flow between locations v i and v j under a set of business events.
The pedestrian flow history dataset is defined as M A , M A = { ( M 1 , A 1 , R 1 ) , ( M 2 , A 2 , R 2 ) , ( M i , A i , R i ) , } , where M i , A i and R i denote the pedestrian flow matrix, the set of business events and the position matrix of each business event during the ith time window ( i [ 1 , n ] ), respectively. Here, R i and M i are of equal size. The attraction of business events in A i and their interactions influence the pedestrian flow matrix M i . It should be noted that all M A data is collected either during the week or on the weekends (holiday).

3.2. Problem Definition

Figure 2a depicts a heat map of crowds for a district. The business event, i.e., a festival booth in Figure 2b, is taken on the route (The terms “route” and “edge” are used interchangeably in this paper.), and the business event in Figure 2c is taken on a route between the POIs in the campus. The green arrows in Figure 2b,c reflect the pedestrian flow on the original route, while the blue signs represent the event on the route. The problem of pedestrian flow prediction with business events is defined as follows:
Problem 1
(Pedestrian flow prediction with business events.). Given the historical data M A with ( s 1 ) time windows, the set of business events A s at the s t h time window, and the position matrix R s of A s , we aim to predict the pedestrian flow matrix M s at the s t h time window.
Business events have an impact on the pedestrian flow matrix. Typically, business events attract pedestrians to visit. Figure 2b shows how visitors are drawn to the events. The phenomenon demonstrates that business events can attract the pedestrian flow for the routes they place. As a result, pedestrian traffic modeling will help operators make decisions about where to arrange business events. The challenge is that different business events placed on different routes result in varying distributions of the overall pedestrian flow.
Figure 2. An example of a business event: (a) Heat map of crowds; (b) Festival booth; (c) Business event in campus.
Figure 2. An example of a business event: (a) Heat map of crowds; (b) Festival booth; (c) Business event in campus.
Sensors 22 07478 g002

4. Predicting Pedestrian Flow

This section presents an Attraction-Based Matrix Factorization (ABMF) model to overcome the aforementioned challenge. Figure 3 outlines the structure of pedestrian flow prediction with business events, including the optimization and prediction procedures.

4.1. ABMF Model

Unlike the traditional matrix factorization [37], we construct a location–location matrix, where each m i j represents the pedestrian flow between location v i to location v j . Let X R z × n and Y R z × n be the latent location feature matrix and the latent factor feature matrix [38], with column vectors X i and Y i representing the z-dimensional location-specific and factor-specific feature vectors of location v i . In particular, the vector of location feature represents the location’s preference for corresponding properties, and the vector of factor feature captures the properties level of location v i , where each property level indicates the ability to attract visitors. For instance, X R 3 × 2 and Y R 3 × 2 , X 1 = ( 0.8 , 0.6 , 0.1 ) T , Y 1 = ( 200 , 60 , 30 ) T . The elements in X 1 represent the preference for three categories of location 1, which correspond to the degree of belonging to a category. The elements in Y 1 represent the ability to attract visitors for corresponding categories of location 1. Therefore, the location 1 can attract 0.8 × 200 + 0.6 × 60 + 0.1 × 30 = 199 visitors. Similarly, we also construct two matrices for business events. We denote H R z × | C | and F R z × | C | as the latent business event feature matrix and the latent factor feature matrix of business event, where the column vectors H j and F j represent the z-dimensional business event-specific and factor-specific feature vectors of business event a j . C denotes a set of categories C = { c 1 , c 2 , , c | C | } and | C | is the amount of categories of business events. The business event feature vector, in particular, reflects the preference of the business event for corresponding properties, whereas the factor feature vector captures the property level of the business event a i , which is the attractive crowd index for each property. We create an undirected graph G = ( V , E ) , where V = { v 1 , v 2 , , v | V | } signifies the set of the locations and e i j represents the edge (route) between v i and v j . Let A = { a 1 , a 2 , , a | A | } denote the set of all the business events. If the visitors are visiting at v i , and then they are attracted to visit v j , the pedestrian flow from v i to v j is formed. When a business event a k ( a k A ) takes place on e i j , the pedestrian flow on e i j includes four aspects: (1) the pedestrian flow from v j to v i ; (2) the pedestrian flow from v i to v j ; (3) the pedestrian flow that a k attracts from all the locations; (4) the pedestrian flow that a k attracts from other business events, excluding the pedestrian flow that other business events attract from a k . As a result, we integrate X , Y , H , F together to model the pedestrian flow between v i and v j .
The current set of a business event is denoted as A , where A A . The predicted pedestrian flow m ^ i j between location v i and location v j is defined as follows:
m ^ i j = X i T Y j + X j T Y i + I i j a k s V 1 | V | X s T F f ( a k ) + I i j a k ( a m A / a k 1 | A / a k | H f ( a m ) T F f ( a k ) a m A / a k 1 | A / a m | H f ( a k ) T F f ( a m ) )
where the function f ( a k ) converts a k into its category index. The indicator function, I i j a k , is equal to 1 if the a k takes place on e i j and to 0 otherwise. A / a k represents the business events set excluding a k , while A / a m represents the set of business events excluding a m . 1 V is the attraction influence probability of a k across all POIs. The attraction influence probability of a k on other business events is 1 | A / a k | , while the attraction influence probability of other business events on a k is 1 | A / a m | . Note that m i i = 0 ( i [ 1 , n ] ) . If there is no event between locations v i and location v j , m i j is represented as: m ^ i j = X i T Y j + X j T Y i , as shown in Figure 4.
Inspired by [38], the elements in the pedestrian flow matrix are supposed to be drawn from a Gaussian distribution with the mean as stated in Equation (1) and the precision as λ m . We apply zero-mean spherical Gaussian priors on location and factor feature vectors with precisions of λ x and λ y , respectively. Similarly, we apply zero-mean spherical Gaussian priors with the precision of λ h and λ f to business event and factor feature vectors. Therefore, using a Maximum-a-Posteriori (MAP) estimation, we derive the objective function for X, Y, H, and F in Equation (2), where · F signifies the Frobenius norm to avoid overfitting.
When m ^ i j is substituted for Equation (1), the ABMF model is learned by the following objective function, where m ^ i j is the estimated value as indicated in Equation (1).
L o s s = min X , Y , H , F λ m 2 i = 1 n j = 1 n ( m i j m ^ i j ) 2 + λ x 2 X F 2 + λ y 2 Y F 2 + λ h 2 H F 2 + λ f 2 F F 2

4.2. Parameter Estimation

We utilize stochastic gradient descent (SGD) to learn the objective function. For each training instance, we update the associated parameters as well as the gradient’s ascending direction to learn the model. The gradients of L o s s with regard to the latent factor features are as follows:
L o s s X i = λ m j = 1 n ( m i j m ^ i j ) ( Y j Y i ) + I i j a k F f ( a k ) | V | + λ x X i
L o s s Y j = λ m i = 1 n ( m i j m ^ i j ) ( X i X j ) + λ y Y j
L o s s H c = λ m i = 1 n j = 1 n ( m i j m ^ i j ) I i j a k F ˜ + λ h H c
L o s s F r = λ m i = 1 n j = 1 n ( m i j m ^ i j ) I i j a k H ˜ + λ f F r
We define the following variables to simplify the formula:
F ˜ = | D 1 | | A / a k | F f ( a k ) a m A / a k 1 | A / a m | F f ( a m )
H ˜ = 1 | V | s = 1 | V | X s + a m A / a k 1 | A / a k | H f ( a m ) | D 2 | | A | 1 H f ( a k )
where D 1 is the business events set where a m D 1 and satisfies a m A / a k & f ( a m ) = c ; D 2 is the business events set where a m D 2 and satisfies a m A / a k & f ( a m ) = r .
The following are the details of Algorithm 1: the inputs consist of the abstracted graph G = ( V , E ) , the historical data M A , the set of all the business events A, regularization coefficients, the dimension z, business events category number | C | and the learning rate η . In lines 1–7, for each pair of historical data in M A , it randomly initializes X , Y , H , F . Then, each column vector of X , Y , H ,   and   F is simultaneously updated as follows based on the gradients and step size η :
X i X i η L o s s X i
Y j Y j η L o s s Y j
H c H c η L o s s H c
F r F r η L o s s F r
After that, the algorithm computes the learned matrices X i , Y i , H i , F i for each pair of ( M i , A i ) . In lines 8–12, after training all the samples, for each X a b i ( i [ 1 , n ] ) , the mean value is calculated excluding the maximal X a b i and the minimum X a b i , because we want to lessen the impact of the maximum and the minimum on the mean value. Then, we obtain the final optimized X and apply the same method to calculate Y , H , F .
Algorithm 1: The ABMF Optimization
Sensors 22 07478 i001
When making predictions, the operator plans to allocate a set of business events A s in the region with the position matrix R s . Based on the learnt matrices, we use Equation (1) to predict each element of the pedestrian flow matrix.

5. Extension: Route Recommendation for Business Events

In this section, we make an extension of the pedestrian flow prediction with business events. In Problem 1, we define the problem of predicting pedestrian flow under a set of business events whose position is known in advance. In some cases, the operator only knows which event he plans to hold but doesnot know which route fits the event. Thus, it is essential to recommend routes for a business event. This is why we make the extension.
In the following subsections, we first introduce the concepts and problems. In order to recommend routes for business event and attract more pedestrian flows, we then apply a Skip-gram-based model to represent each route like representing word vector in Section 5.2. Finally, we improve the pair-wise ranking loss to a flow-aware-based method by considering the flow factors in Section 5.3, which aims to learn events’ latent representations for top-N route recommendation by minimizing the F-WARP loss.

5.1. Concepts and Problem Definition

The set of business events is denoted by A = { a 1 , a 2 , . . . } . The set of all routes in a travel graph is denoted by E = { e 1 , e 2 , . . . } . The travel trajectories of users are defined as follows:
Definition 2
(Travel trajectories). The set of travel trajectories is denoted by S = { S 1 , S 2 , . . . } , where S i refers to the ith travel trajectory and S i = { e 1 , e 2 , . . . , e i , . . . } , e i denotes a route in a travel trajectory. Note that S i is the subset of E.
Definition 3
(Target route and context route). In a travel trajectory S u , the chosen e i is the target route, and other routes in S u are context routes.
Problem 2
(Routes recommendation for an event). Given a business event a and users’ travel trajectories set E, the task is to recommend routes to hold this event, which aims to attract more people to visit this event so as to bring more benefit to the merchant.
For example, the operators plan to hold a sales promotion and will be recommended routes for this promotion. These routes not only have a high preference degree for this promotion but also bring pedestrian flow to increase potential benefits.

5.2. Learning Route Representations

The workflow to recommend routes for events is presented as flows: We first learn routes’ latent vector representations from context routes in the trajectories. Then, we capture events’ flow-aware preference for routes and learn events’ latent vector representations. At last, we recommend top-N routes for an event by the inner product of the latent factor vectors of an event and routes.
Based on the Skip-Gram model [8], we propose an embedding model to learn route representations, which captures routes’ contextual information from trajectory sequences. We learn the representations of context routes from e i + k to e i k given a target route e i . Here, k is set to control the window size of context, as shown in Figure 5. In our scenario, the collected route corpus is fed into the Skip-gram model. Compared with the word corpus, each route and each travel trajectory correspond to a word and a sentence, respectively. We aggregate all travel trajectories to construct a route corpus. We learn the route embeddings by maximizing the function below:
L m = S u S 1 | S u | e i S u k z k , z 0 log P ( e i + z | e i )
where S is a set containing all travel trajectories. L m aims to maximize the context route’s conditional occurrence likelihood for all travel trajectories.
In addition, we formulate the probability P ( e i + z | e i ) using a softmax function (The bold type in this paper refers to a vector.):
P ( e i + z | e i ) = exp ( v i + z · v i ) e x E exp ( v x · v i )
where D refers to the dimension of the latent space; v i R D × 1 and v i + z R D × 1 refer to the latent representations of the target route e i and one of its context routes e i + z , respectively. It is hard to directly optimize Equation (14) as the size of set E is extremely large. In this work, a negative sampling [8] method is utilized to promote optimization efficiency. For each route e i E , a set of K routes that not occur in e i ’s context window are sampled. We reformulate the objective function as a new form that is easier to minimize,
L m = S u S 1 | S u | e i S u k z k , c 0 ( log σ ( v i + z · v i ) + k = 1 K E k log σ ( v k · v i )
where σ ( · ) refers to the commonly used sigmoid function. The noise distribution P n ( r ) is applied to sample the K negative routes, this distribution is the unigram distribution raised to the 3/4rd power [8]. E ( · ) means calculating the expected value for all negative samples used. We utilize the backpropagation algorithm to fit the Skip-gram model.

5.3. Flow-Aware Recommendation Model

In the above subsection, routes’ latent representations are learned by utilizing the context patterns in a global view, so the various flow characteristics on the routes are not considered. In order to address this problem, we design a flow-aware preference learning model for top-N route recommendations as events’ route selection. It is seen that the events are usually taken on the routes, so the flow size on the route is crucial for events. We define f a , e to denote the pedestrian flow size when an event a is taken on the route e. Note that the larger flow size f on the route e, the more confident that event a fit for the route e. Thus, the events’ preference rankings of the routes can be obtained. For example, f a , e > f a , e indicates the route e ranked higher than the route e for a given event a. In this paper, f a , e = 0 shows a is not fit for e. It can also refer to the fact that a has not been taken on e.
In order to design a flow-aware preference learning model for top-N route recommendation, with the help of pair-wise ranking loss, we utilize the Weighted Approximately Ranked Pairwise (WARP) loss [39] to learn events’ latent representations. By leveraging the precision@N measure, the WARP loss measures the pair-wise violations relying upon routes’ positions in the ranking list. For each event a, its arranged route set and un-arranged route set are denoted as Z a + and Z a , respectively. We define the WARP loss as follows,
L w a r p = a A e Z a + L [ r a n k ( ω ^ a , e ) ]
where r a n k ( ω ^ a , e ) is the rank of an arranged route e Z a + in a’s flow-aware ranking list of routes. r a n k ( ω ^ a , e ) is estimated by e Z a I ( ω ^ a , e ω ^ a , e ) , and I ( · ) is an indicator function. To optimize the WARP loss function, we supplant the discrete indicator function by a continuous margin function [33]: m a x ( 0 , 1 ω ^ a , e + ω ^ a , e ) . L ( · ) is applied to convert a ranking order to a loss value, and L ( s ) = i = 1 s 1 i . ω ^ a , e represents an event a’s preference for a route e, and it can be predicted by our factorization model: ω ^ a , e = a v e , where a R D × 1 is the latent vector of an event a, and v e refers to e’s latent vector that obtained by the Skip-gram model in the Section 5.2. To accommodate the Skip-gram model, we set the dimension of events’ latent representations is the same as routes’ latent representations.
For most events, the number of routes un-arranged is substantially more than the arranged ones. We need to approximate the r a n k function efficiently. For each event a, we sample an un-arranged route e randomly when an arranged route e is given, until the sampled routes do not conform with the margin function. ω ^ a , e is estimated by | Z a | 1 N , where · is floor function to obtain an integer, and N is the quantity of sampling times, | · | is the set’s cardinality.
We need to catch events’ flow-aware preference for routes, so we consider the flow size on the routes. Thus, we propose the F-WARP loss to enhance the WARP loss for more reasonable recommendations. For each pair of positive and negative routes ( e , e ) , we consider adding a weight θ e , e , which is defined as: θ e , e = 1 + γ · ( f a , e f a , e ) , where γ is to control the difference of increasement. If the difference is larger, this pair of routes does not conform to the margin function, so this pair of routes contribute much to the total loss. When we take the weight of the route pair into account, we reformulate the loss function as follows:
L f w a r p = a A e Z a + L [ e Z a { Z a + \ e } m a x ( 0 , θ e , e · ( 1 ω ^ a , e + ω ^ a , e ) ) ] + λ a A a F 2
where λ manages the degree of regularization and aims to avoid over-fitting. In the negative sampling routes process. It is not essential to bind the negative routes to un-arranged routes; the arranged routes of the event with a small flow size can likewise be sampled as negative cases.
We employ stochastic gradient descent (SGD) to learn the latent factors of events. For each arranged route, we sample a negative route to update the event’s latent factors in each iteration. The gradient of L f w a r p about the kth latent factor of a is calculated as follows:
L f w a r p a k = L ( | E | 2 N ) θ e , e ( v e , k v e , k ) + 2 λ a k
The update rule for the latent factor is:
a k a k η L f w a r p a k
where η is the learning rate.
After the latent representations of events have been learned, we evaluate the suitability of an event a for a route e by computing the inner product of a and v e . We recommend the routes for an event by ordering the candidate routes in a descending sort of the predicted scores, and then we select top-ranked N routes for the recommendation. We define the flow-aware route recommendation for business events based on the Skip-gram model as SG-FWARP.

6. Experiments

This section analyzes the pedestrian flow before comparing the performance of the ABMF and SG-FWARP with existing approaches.

6.1. The Experimental Setup

(1) Datasets. In the experiment, we describe three datasets as follows:
  • The first is a surveillance video dataset of the New York Grand Central Station [40], which records people coming and passing by.
  • The second one is a simulation dataset, which simulates pedestrian flow under different business event sets by using the social force model [41]. In the simulation dataset, the number of locations ranges from 10–60, and there are 8 types of business events, such as promotions, exhibitions of cloths, games and so on, each type of business event has a preference to attract visitors.
  • The third one is Foursquare [2,42], a publicly crowdsourced large-scale check-in LBSNs dataset is free available. Foursquare of New York (NYC) has 1,385,223 check-ins, whereas Foursquare of Tokyo (TKY) has 573,703 check-ins. We obtain the check-in flows and identify the business events using the dataset’s visit counts for the routes.
(2) Experimental Settings. We choose 80 % of the datasets to train the ABMF model, and the remaining 20 % to test. We use cross-validation in training SG-FWARP. λ m is set to 0.0001 , whereas λ x , λ y , λ h and λ f are all set to 0.01 . z is assigned to different values based on the simulation data’s scale. η is set as 0.001 first, and we decrease its value when iterations are sufficient. λ is set to 0.01 . We recommend top-10 routes for a business event. We select the regularization coefficient by experimental value. The process is as follows: we first determine the learning rate; then we determine the order of magnitude of λ m (such as 0.0001 or 0.01 , etc). Finally, we fine-tune it further to obtain a proper value of a regularization coefficient.

6.2. Metrics

  • The ABMF model is evaluated according to matrix similarity and root-mean-square error (RMSE), which are defined as:
    s i m ( M , N ) = i j ( m i j M ¯ ) ( n i j N ¯ ) ( i j ( m i j M ¯ ) 2 ) ( i j ( n i j N ¯ ) 2 )
    R M S E = 1 n i = 1 n j = 1 n ( m i j n i j ) ) 2
    where s i m ( M , N ) represents the similarity degree between M and N adjusted to 0 1 . Both M and N are n-order matrices that represent the ground truth and predicted pedestrian flow matrix, respectively. The average values of matrix M and matrix N are represented by M ¯ and N ¯ , respectively.
  • To evaluate the performance of route recommendation for events, we use precision@N, which refers to the ratio of the successfully predicted routes to the top-N recommendations. We use Mean Reciprocal Rank (MRR) as another metric. This ranking metric measures the recommendation accuracy by finding out how far the first successfully predicted route is from the top of the recommendation list. MRR is defined as follows:
    M R R = 1 A i = 1 | A | 1 r i
    where | A | is the size of the event set. For the ith event, r i refers to the ranking position of the first route in the recommendation list in the ground-truth result. Note that all experiments are conducted 10 times for latent factor models, and we report the averaged results.

6.3. Baseline methods

  • The ABMF model is used to predict pedestrian flow based on a set of business events. We compare the prediction results of the ABMF model with the other three prediction methods. The baseline algorithms include vSVR [21], PMF [37,43], and MOHER [18] (1) vSVR is an application of SVM (Support Vector Machine) to regression problems, which can be applied to predict the pedestrian flow density, and then covered to the number of individuals. (2) PMF is an effective method that is frequently used as a baseline in current work [44], which factors the pedestrian flow matrix into two feature matrices. (3) MOHER is to predict the potential crowd flow in a certain mode, which uses the LTSM module to predict the sequential flow. However, the baseline algorithms cannot reflect the diverse influence of business events, because they fuse the influence of business events into the pedestrian flow.
  • We compare SG-FWARP with three top-N recommendation methods. (1) The first is WRMF [45], which is the weighted regularized matrix factorization model designed to handle implicit feedback data (i.e., an event taken on a route or not) for top-N recommendations. This method is used as a baseline in the latest work [46]. (2) The second one is WARP-MF [47]. This is a pairwise ranking method that utilizes matrix factorization to minimize the basic WARP loss. The latent factors of events and routes are learned by randomly sampling the positive and negative route pairs. (3) The third one is similarity pairwise ranking matrix factorization (SPRMF) [48]. This method uses a new penalty to eliminate the differences in the scores between popular and personalized items based on their similarity.

6.4. Analysis of Pedestrian Flow

A 5-min video collected from New York Grand Central Station is used to analyze the pedestrian flows. We use KLT keypoint tracker [49] to extract trajectories. Figure 6a shows the clustered main coarse-grained trajectories. Then, the extracted trajectories are shown in Figure 6b. A red ellipse represents the attractive region. The pedestrians congregate around the attractive region, indicating that something draws pedestrians to visit. As shown in Figure 6c, the horizontal axis represents the indexes of the clustered routes extracted from the video in Figure 6b. The vertical axis reflects the number of pedestrian flows on each route during the video. According to the statistical results, the number of pedestrian flow on route 1 is higher than that of other routes, indicating that the attractive region draws in the pedestrians for route 1. As a result, it makes sense to introduce the attractive crowd index for each property in the ABMF model.

6.5. Experimental Results

We present the experimental results of the two problems separately.

6.5.1. Experimental Results for the ABMF

Based on the simulation pedestrian dataset with business events, the pedestrian flow matrix without/with business events are mapped on the heat map. Figure 7a shows the original pedestrian flow matrix with five locations, and Figure 7b shows the pedestrian flow matrix of five locations when arranging two business events on e 14 and e 54 , respectively. It is obvious that the pedestrian flow values on route e 14 and e 54 in Figure 7b are larger than the corresponding values in Figure 7a, indicating that the scheduled business events draw additional pedestrian flow for the associated routes.
Figure 8 depicts the performance of various algorithms with various scales of simulation datasets. In Figure 8a, the number of locations corresponds to the scale of the simulation dataset, and each business event set has the same categories under all the simulation datasets. The dimension of latent features is set differently with the different scales of the datasets. The number of locations in the dataset ranges from 10 to 60, and the corresponding dimension is set from 5 to 30. Overall, the RMSE of the ABMF is lower than that of the other three methods. The smaller RMSE is, the more accurate the algorithm is. The reason is that ABMF considers the ability of categories to attract pedestrian flow and refines the source of the flow of a route. The performance of PMF is better than that of vSVR because PMF considers potential factors in matrix factorization. In Figure 8b, the similarity between the predicted matrix by the ABMF and the ground truth matrix is higher than that of the other three algorithms overall. Regarding the metric of similarity, its range is 0 to 1, and the greater the similarity, the closer the predicted result is to the ground truth. As the ABMF takes into account the entire flow sources of a route with a business event, the ABMF is utilized for pedestrian flow prediction with the diversity of business events set. The result shows that the ABMF outperforms vSVR, PMF, and MOHER, because the three methods fuse the influence of business in the model, while ABMF models the attraction of business events to the pedestrian flow in the design. MOHER performs better than vSVR and PMF because MOHER uses LSTM to mix different flow modes and then predict the pedestrian flow. The performance of MOHER is not as good as ABMF, because MOHER is more suitable for the mix of different transportation modes, and only predicting pedestrian flow can not play a full role. In Figure 9, ABMF performs well on large-scale real datasets, and its performance is superior to vSVR, PMF, and MOHER. The analysis in Figure 9 with the Foursquare dataset is consistent with that in Figure 8.
In this experiment, the simulation data is associated with 20 locations and 8 kinds of business events, and the dimension of location latent factor z is set to different values in order to evaluate the metrics of similarity and RMSE. Figure 10 displays the performance of various algorithms with varying dimensions of the latent factor. As illustrated in Figure 10b, the similarity of both four methods rises with the dimension increasing. In Figure 10a, the RMSE of the ABMF is lower than that of the other three methods with the dimension increasing, indicating that a larger dimension model is more accurate in predicting pedestrian flow, and the ABMF model performs the best of the four algorithms. When the dimension is set to 30, the similarity of the ABMF remains close to 1, indicating that the ABMF performs well in terms of accuracy and stability. The analysis in Figure 11 under the Foursquare dataset is similar to that in Figure 10. Figure 11 shows that when the dimension increases, the performance of ABMF in terms of RMSE and similarity changes little, indicating the solid performance of the ABMF model. The increase in the potential feature dimension plays a promotive role in the performance of ABMF.
Figure 12 displays the influence of business event diversity on the routes in the simulation dataset. There are four different types of business events in the business events set; each business event is associated with a type. We select three routes from the simulation dataset, and we place 0 , 1 , 2 , 3 or 4 business events on each route to obtain the changed pedestrian flow matrix, respectively. Each time, the business events are placed on the three routes simultaneously. The results show that the pedestrian flow size of each route rises as the type of business events increases.

6.5.2. Experimental Results for SG-FWARP

We first analyze how the context window size impacts the performance of SG-FWARP in Table 2. The dimensionality of latent factor vectors D is set as 300 in the experiments. The size of the context window varies from 1 to 6. From Table 2, we observe that precision and MRR increase first with the increasing context window size, then decrease after arriving at a larger value at the window size of 3 or 4. The reason is as follows: a larger window size represents the target route’s context and could be more comprehensively modeled. Compared with the large-scale text corpus where the Skip-gram is typically used, the route corpus we used is sparse, so a small context window size is proper to model the route’s context influence in our experiments.
Next, we report the impact of dimensionality on the latent vectors in Figure 13. In the SG-FWARP method, we set the dimensionality from 50 to 400 with 50 as the increment. Figure 13 shows that a higher latent vector dimensionality brings higher precision and MRR. This trend reflects that higher latent vector dimensionality helps to capture more latent factors and more accurately represent routes and events. However, we find that precision and MRR are stable after achieving a certain threshold (e.g., 300). As a consequence, an optimal latent vector dimensionality can be obtained from empirical value, so we can obtain a high recommendation quality with proper computational overheads. In the following comparison experiment, we set 300 as the default dimensionality of a latent vector.
In the comparison experiment, we compare the performance of our proposed method SG-FWARP with three baseline methods. By using cross-validation, for the methods WRMF, WARP-MF, and SPRMF, we set the dimension of latent factor vector, γ , and regularization coefficient λ as 200, 1, and 0.01 , respectively. For SG-FWARP, the latent factor vector dimensionality, γ , and regularization coefficient are set to 300, 1, and 0.01 , respectively. In addition, the optimal context window size is set for different datasets.
Table 3 shows the experimental results on two datasets when evaluating top-10 route recommendations as events’ route recommendations. We can see that the SG-FWARP model outperforms WARP-MF, WRMF, and SPRMF, because SG-FWARP captures the semantic information of sequence routes. The main difference is that semantic information (i.e., categories of routes) is considered to learn the latent preferences, thus improving the performance. The other three algorithms WARP-MF, WRMF, and SPRMF fuse the flow factor in the matrix factorization, while SG-FWARP model introduces the flow factor into the WARP loss to improve the reasonability in flow-aware route recommendation for business events.
Overall, SG-FWARP is superior to the three baseline methods by two main design mechanisms. First, we apply the Skip-gram model to learn the routes’ latent representations to capture the context (before and after the target route) of routes. Second, based on the learned route latent representations, the F-warp loss function combines the flow size of the routes with a pairwise ranking algorithm to learn users’ latent representations for flow-aware route recommendation.

6.6. Discussion

Overall, the problem of route recommendation for business events happens before the pedestrian flow prediction. Route recommendation for business events is to assist operators in making decisions when the operators have no idea where to place the business event. We explained the experimental results in the above subsections in detail. We first verify the performance of ABMF in predicting pedestrian flow, then we show the experimental results to demonstrate the feasibility of recommending routes for business events. However, suppose the operators apply SG-FWARP every time they make a decision. In this case, the problem of pedestrian flow prediction with business events and the problem of route recommendation for business events will affect each other in the long term. Below, we conclude our paper and present future work.

7. Conclusions

In this paper, we first present the ABMF model to predict the pedestrian flow with business events, which is essential because pedestrian flow can bring potential profits to operators; the highlight of the ABMF model is that it introduces the attraction index of different categories to pedestrians in matrix factorization. Then, we learn route representations based on the Skip-gram model, and we improve the basic pair-wise ranking loss to a flow-aware-based method by considering the flow factors, which aims to learn events’ latent representations for route recommendation. To the best of our knowledge, we are the first to predict pedestrian flow by considering business events. We evaluate the performance of the ABMF and SG-FWARP; both are superior to the baseline methods, which confirms the ABMF and SG-FWARP models are better and more applicable to pedestrian flow prediction with business events and route recommendation for business events. The limit of this paper is that we did not consider the interaction effect between route recommendation for the events and pedestrian flow prediction; the two problems can affect each other in the long run. In the future, we will consider the joint optimization of the two problems and find the optimal solutions under the interaction effect by utilizing Markov Process. Furthermore, we will also consider the time-aware factors to improve the ABMF and SG-FWARP models.

Author Contributions

Conceptualization, J.G. and C.S.; methodology, J.G., C.S. and W.J.; software, J.G. and Z.R.; validation, J.G. and Z.R.; formal analysis, J.G., L.L. and M.L.; investigation, J.G. and C.S.; writing—original draft preparation, J.G.; writing—review and editing, J.G., L.L., W.J., M.L. and C.S.; visualization, J.G. and Z.R.; supervision, C.S.; project administration, C.S.; funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China under Grant 2021YFB3101302 and 2021YFB3101303; the National Natural Science Foundation of China under Grant No. 62020106013 and No. 61872061; the Science and Technology Achievements Transformation Demonstration Project of Sichuan Province of China No. 2018CC0094; and the Fundamental Research Funds for the Central Universities No. ZYGX2019J075.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the first author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, Z.; Ma, H.; Guo, B.; Yang, Z. Crowdsensing 2.0. Commun. ACM 2021, 64, 76–80. [Google Scholar] [CrossRef]
  2. Gu, J.; Song, C.; Jiang, W.; Wang, X.; Liu, M. Enhancing personalized trip recommendation with attractive routes. In Proceedings of the AAAI Conference on Artificial Intelligence, New York, NY, USA, 7–12 February 2020; Volume 34, pp. 662–669. [Google Scholar]
  3. Estrada, R.; Mizouni, R.; Otrok, H.; Mourad, A. Task coalition formation for Mobile CrowdSensing based on workers’ routes preferences. Veh. Commun. 2021, 31, 100376. [Google Scholar] [CrossRef]
  4. Asahara, A.; Sato, N.; Nomiya, M. Pedestrian-Flow Analysis System for Improving Layout of Exhibitions. In Proceedings of the 14th International Symposium on Spatial and Temporal Databases, Hong Kong, China, 26–28 August 2015. [Google Scholar]
  5. May, M.; Scheider, S.; Schulz, D.; Hecker, D. Pedestrian flow prediction in extensive road networks using biased observational data. In Proceedings of the 16th ACM Sigspatial International Conference on Advances in Geographic Information Systems, Irvine, CA, USA, 5–7 November 2008. [Google Scholar]
  6. Wang, X.; Leckie, C.; Chan, J.; Lim, K.H.; Vaithianathan, T. Improving Personalized Trip Recommendation by Avoiding Crowds. In Proceedings of the 25th ACM International on Conference on Information and Knowledge Management, Indianapolis, IN, USA, 24–28 October 2016. [Google Scholar]
  7. Rong, D.; Yu, Z.; Tao, M.; Wang, Z.; Guo, B. Predicting activity attendance in event-based social networks: Content, context and social influence. In Proceedings of the 2014 ACM International Joint Conference on Pervasive and Ubiquitous Computing, Seattle, WA, USA, 13–17 September 2014. [Google Scholar]
  8. Mikolov, T.; Sutskever, I.; Chen, K.; Corrado, G.; Dean, J. Distributed representations of words and phrases and their compositionality. Adv. Neural Inf. Process. Syst. 2013, 26, 3111–3119. [Google Scholar]
  9. Sahaleh, S.; Bierlaire, M.; Farooq, B.; Danalet, A.; Hänseler, F. Scenario analysis of pedestrian flow in public spaces. In Proceedings of the STRC 12th Swiss Transport Research Conference, Ascona, Switzerland, 2–4 May 2012. [Google Scholar]
  10. Widhalm, P.; Brandle, N. Learning Major Pedestrian Flows in Crowded Scenes. In Proceedings of the 2010 20th International Conference on Pattern Recognition, Istanbul, Turkey, 23–26 August 2010. [Google Scholar]
  11. Cooper, C.H.; Harvey, I.; Orford, S.; Chiaradia, A.J. Using multiple hybrid spatial design network analysis to predict longitudinal effect of a major city centre redevelopment on pedestrian flows. Transportation 2021, 48, 643–672. [Google Scholar]
  12. Tanaka, R.; Takahashi, H. Multi-agent simulation approach of pedestrian flow with group walking models. ICIC Exp. Lett. B Appl. 2020, 11, 363–371. [Google Scholar]
  13. Kaminka, G.A.; Fridman, N. Simulating Urban Pedestrian Crowds of Different Cultures. ACM Trans. Intell. Syst. Technol. 2018, 9, 1–27. [Google Scholar] [CrossRef]
  14. Zheng, Y.; Capra, L.; Wolfson, O.; Yang, H. Urban Computing: Concepts, Methodologies, and Applications. ACM Trans. Intell. Syst. Technol. 2014, 5, 1–55. [Google Scholar] [CrossRef]
  15. Fan, Z.; Song, X.; Adachi, R.; Adachi, R. Citymomentum: An online approach for crowd behavior prediction at a citywide level. In Proceedings of the 2015 ACM International Joint Conference on Pervasive and Ubiquitous Computing, Osaka, Japan, 7–11 September 2015. [Google Scholar]
  16. Song, X.; Zhang, Q.; Sekimoto, Y.; Shibasaki, R. Prediction of human emergency behavior and their mobility following arge-scale disaster. In Proceedings of the 20th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, New York, NY, USA, 24–27 August 2014. [Google Scholar]
  17. Hoang, M.X.; Zheng, Y.; Singh, A.K. FCCF: Forecasting citywide crowd flows based on big data. In Proceedings of the 24th ACM Sigspatial International Conference on Advances in Geographic Information Systems, Burlingame, CA, USA, 31 October–3 November 2016. [Google Scholar]
  18. Zhou, Q.; Gu, J.; Lu, X.; Zhuang, F.; Zhao, Y.; Wang, Q.; Zhang, X. Modeling heterogeneous relations across multiple modes for potential crowd flow prediction. In Proceedings of the AAAI Conference on Artificial Intelligence, Virtual, 2–9 February 2021; Volume 35, pp. 4723–4731. [Google Scholar]
  19. Zhang, J.; Zheng, Y.; Qi, D. Deep Spatio-Temporal Residual Networks for Citywide Crowd Flows Prediction. In Proceedings of the 31st AAAI Conference on Artificial Intelligence, San Francisco, CA, USA, 4–9 February 2016. [Google Scholar]
  20. Eravci, B.; Bulut, N.; Etemoglu, C.; Ferhatosmanoǧlu, H. Location Recommendations for New Businesses Using Check-in Data. In Proceedings of the 2016 IEEE 16th International Conference on Data Mining Workshops (ICDMW), Barcelona, Spain, 12–15 December 2016. [Google Scholar]
  21. Ma, Y.; Bai, G. Short term prediction of crowd density using v-SVR. In Proceedings of the 2010 IEEE Youth Conference on Information, Computing and Telecommunications, Beijing, China, 28–30 November 2010. [Google Scholar]
  22. Duan, J.; Wang, L.; Long, C.; Zhou, S.; Zheng, F.; Shi, L.; Hua, G. Complementary Attention Gated Network for Pedestrian Trajectory Prediction. In Proceedings of the AAAI Conference on Artificial Intelligence, Virtual, 22 February–1 March 2022. [Google Scholar]
  23. Li, Y.; Zheng, Y.; Ji, S.; Wang, W.; Hou, U.L.; Gong, Z. Location selection for ambulance stations: A data-driven approach. In Proceedings of the 23rd SIGSPATIAL International Conference on Advances in Geographic Information Systems, Seattle, WA, USA, 3–6 November 2015. [Google Scholar]
  24. Bao, J.; He, T.; Ruan, S.; Li, Y.; Zheng, Y. Planning Bike Lanes based on Sharing-Bikes’ Trajectories. In Proceedings of the 23rd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Halifax, NS, Canada, 13–17 August 2017. [Google Scholar]
  25. Liu, D.; Weng, D.; Li, Y.; Bao, J.; Zheng, Y.; Qu, H.; Wu, Y. SmartAdP: Visual Analytics of Large-scale Taxi Trajectories for Selecting Billboard Locations. IEEE Trans. Vis. Comput. Graph. 2017, 23, 1–10. [Google Scholar] [CrossRef]
  26. Zhao, S.; King, I.; Lyu, M.R.; Zeng, J.; Yuan, M. Mining Business Opportunities from Location-based Social Networks. In Proceedings of the 40th International ACM SIGIR Conference on Research and Development in Information Retrieval, Tokyo, Japan, 7–11 August 2017. [Google Scholar]
  27. Lin, J.; Oentaryo, R.J.; Lim, E.P.; Vu, C.; Vu, A.; Kwee, A.T.; Prasetyo, P.K. A Business Zone Recommender System Based on Facebook and Urban Planning Data. In European Conference on Information Retrieval; Springer: Cham, Switzerland, 2016; pp. 641–647. [Google Scholar]
  28. Sonosy, O.A.; Rady, S.; Badr, N.L.; Hashem, M. Exploiting location based social networks in business predictions. In Proceedings of the 2015 11th International Conference on Innovations in Information Technology (IIT), Dubai, United Arab Emirates, 1–3 November 2015. [Google Scholar]
  29. Hu, B. Recommendation in Location-based Social Networks. Ph.D. Thesis, Simon Fraser University, Vancouver, BC, Canada, 2014. [Google Scholar]
  30. Lu, K.; Zhang, Y.; Zhang, L.; Wang, S. Exploiting User and Business Attributes for Personalized Business Recommendation. In Proceedings of the 38th International ACM SIGIR Conference on Research and Development in Information Retrieval, Santiago, Chile, 9–13 August 2015. [Google Scholar]
  31. Ye, J.; Zhu, Z.; Cheng, H. What’s Your Next Move: User Activity Prediction in Location-based Social Networks. In Proceedings of the 2013 SIAM International Conference on Data Mining; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2013. [Google Scholar]
  32. Zhang, J.D.; Li, Y.; Li, Y. LORE: Exploiting sequential influence for location recommendations. In Proceedings of the 22nd ACM SIGSPATIAL International Conference on Advances in Geographic Information Systems, Dallas, TX, USA, 4–7 November 2014. [Google Scholar]
  33. Liu, X.; Liu, Y.; Li, X. Exploring the context of locations for personalized location recommendations. In Proceedings of the Twenty-Fifth International Joint Conference on Artificial Intelligence, New York, NY, USA, 9–15 July 2016. [Google Scholar]
  34. Tang, D.; Qin, B.; Yang, Y.; Yang, Y. User modeling with neural network for review rating prediction. In Proceedings of the Twenty-Fourth International Joint Conference on Artificial Intelligence, Buenos Aires, Argentina, 25–31 July 2015. [Google Scholar]
  35. Shanshan, F.; Gao, C.; Bo, A.; Yeow Meng, C. POI2Vec: Geographical Latent Representation for Predicting Future Visitors. In Proceedings of the Thirty-First AAAI Conference on Artificial Intelligence, San Francisco, CA, USA, 4–9 February 2017. [Google Scholar]
  36. Zhao, S.; Zhao, T.; King, I.; Lyu, M.R. Geo-Teaser: Geo-Temporal Sequential Embedding Rank for Point-of-interest Recommendation. In Proceedings of the 26th International Conference on World Wide Web Companion, Perth, Australia, 3–7 April 2017. [Google Scholar]
  37. Salakhutdinov, R.; Mnih, A. Probabilistic Matrix Factorization. In Proceedings of the 20th International Conference on Neural Information Processing Systems, Vancouver, BC, Canada, 3–6 December 2007. [Google Scholar]
  38. Li, H.; Hong, R.; Zhu, S.; Yong, G. Point-of-Interest Recommender Systems: A Separate-Space Perspective. In Proceedings of the 2015 IEEE International Conference on Data Mining, Atlantic City, NJ, USA, 14–17 November 2015. [Google Scholar]
  39. Weston, J.; Bengio, S.; Usunier, N. Large scale image annotation: Learning to rank with joint word-image embeddings. Mach. Learn. 2010, 81, 21–35. [Google Scholar] [CrossRef]
  40. Zhou, B.; Wang, X.; Tang, X. Understanding Collective Crowd Behaviors:Learning a Mixture Model of Dynamic Pedestrian-Agents. In Proceedings of the 2012 IEEE Conference on Computer Vision and Pattern Recognition, Providence, RI, USA, 16–21 June 2012. [Google Scholar]
  41. Cheng, L.; Yarlagadda, R.; Fookes, C.; Yarlagadda, P.K. A review of pedestrian group dynamics and methodologies in modelling pedestrian group behaviours. World J. Mech. Eng. 2014, 1, 1–13. [Google Scholar]
  42. Lu, E.H.C.; Chen, C.Y.; Tseng, V.S. Personalized trip recommendation with multiple constraints by mining user check-in behaviors. In Proceedings of the 20th International Conference on Advances in Geographic Information Systems, Redondo Beach, CA, USA, 6–9 November 2012. [Google Scholar]
  43. Pujahari, A.; Sisodia, D.S. Pair-wise preference relation based probabilistic matrix factorization for collaborative filtering in recommender system. Knowl.-Based Syst. 2020, 196, 105798. [Google Scholar] [CrossRef]
  44. Zhang, K.; Wang, Z.; Liang, J.; Zhao, X. A Bayesian matrix factorization model for dynamic user embedding in recommender system. Front. Comput. Sci. 2022, 16, 165346. [Google Scholar] [CrossRef]
  45. Hu, Y.; Koren, Y.; Volinsky, C. Collaborative filtering for implicit feedback datasets. In Proceedings of the 2008 Eighth IEEE International Conference on Data Mining, Pisa, Italy, 15–19 December 2008; pp. 263–272. [Google Scholar]
  46. Fan, R.; Chen, J.; Zhang, J.; Lian, D.; Chen, E. Improving Implicit Alternating Least Squares with Ring-based Regularization. In Proceedings of the 45th International ACM SIGIR Conference on Research and Development in Information Retrieval, Madrid, Spain, 11–15 July 2022; pp. 102–111. [Google Scholar]
  47. Zhao, S.; King, I.; Lyu, M.R. Geo-pairwise ranking matrix factorization model for point-of-interest recommendation. In Proceedings of the International Conference on Neural Information Processing; Springer: Berlin/Heidelberg, Germany, 2017; pp. 368–377. [Google Scholar]
  48. Liu, J.; Yang, Z.; Li, T.; Wu, D.; Wang, R. SPR: Similarity pairwise ranking for personalized recommendation. Knowl. Based Syst. 2022, 239, 107828. [Google Scholar] [CrossRef]
  49. Tomasi, C.; Kanade, T. Detection and Tracking of Point Features. Int. J. Comput. Vis. 1991, 9, 137–154. [Google Scholar] [CrossRef]
Figure 1. Business event increases the pedestrian flow of a route.
Figure 1. Business event increases the pedestrian flow of a route.
Sensors 22 07478 g001
Figure 3. Framework of pedestrian flow prediction with business events.
Figure 3. Framework of pedestrian flow prediction with business events.
Sensors 22 07478 g003
Figure 4. Vector representations for pedestrian flow matrix M.
Figure 4. Vector representations for pedestrian flow matrix M.
Sensors 22 07478 g004
Figure 5. An instance of the context window of a trajectory. The window includes 4 routes visited before and after the ith route.
Figure 5. An instance of the context window of a trajectory. The window includes 4 routes visited before and after the ith route.
Sensors 22 07478 g005
Figure 6. The statistics in New York Grand Central Station. (a) Trajectory clustering; (b) The extracted trajectories; (c) The pedestrian flow number.
Figure 6. The statistics in New York Grand Central Station. (a) Trajectory clustering; (b) The extracted trajectories; (c) The pedestrian flow number.
Sensors 22 07478 g006
Figure 7. Heat map of pedestrian flow matrix: (a) The original matrix; (b) The changed matrix.
Figure 7. Heat map of pedestrian flow matrix: (a) The original matrix; (b) The changed matrix.
Sensors 22 07478 g007
Figure 8. Performance comparison with different scale of simulation dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Figure 8. Performance comparison with different scale of simulation dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Sensors 22 07478 g008
Figure 9. Performance comparison with different scale of Foursquare(NYC) dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Figure 9. Performance comparison with different scale of Foursquare(NYC) dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Sensors 22 07478 g009
Figure 10. Performance with different dimensions on simulation dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Figure 10. Performance with different dimensions on simulation dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Sensors 22 07478 g010
Figure 11. Performance with different dimensions on Foursquare (NYC) dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Figure 11. Performance with different dimensions on Foursquare (NYC) dataset: (a) Performance on RMSE; (b) Performance on Similarity.
Sensors 22 07478 g011
Figure 12. The impact of diversity events on the route in the simulation dataset.
Figure 12. The impact of diversity events on the route in the simulation dataset.
Sensors 22 07478 g012
Figure 13. Performance of different latent factor vector dimensionality (top-10 recommendations): (a) Precision; (b) MRR.
Figure 13. Performance of different latent factor vector dimensionality (top-10 recommendations): (a) Precision; (b) MRR.
Sensors 22 07478 g013
Table 1. Frequently-used notations.
Table 1. Frequently-used notations.
NotationInterpretation
V / E the set of locations/routes, v i V , | V | = n , e i j E
Athe business events set, a i A
A the current business events set
M A the set of pedestrian flow history data
M i the ith pedestrian flow matrix
A i the ith the business events set, A i A
R i the position matrix of A i
Cthe set of categories, c i C
Xthe latent location feature matrix
Hthe latent business event feature matrix
Ylatent factor feature matrix of location
Flatent factor feature matrix of business event
f ( a i ) the function to map a i to its category index
λ x , λ y , λ h , λ f the regularization parameters of X, Y, H and F
λ m the reciprocal of the variance
η the learning rate
sthe predicted pedestrian flow matrix index
zthe latent feature dimension
S i the ith trajectory, S i S
v / a the vector of route/event
f a , e the flow size with event a on the route e
L m loss function of learning route representations
L w a r p loss function of flow-aware WARP loss
Table 2. Performance of different context window size on SG-FWARP (top-10 recommendations).
Table 2. Performance of different context window size on SG-FWARP (top-10 recommendations).
DatasetMetricContext Window Size
123456
Foursquare (NYC)Precision@100.20240.23420.25480.27460.23830.2262
MRR@100.56560.61980.68410.70710.61080.5026
Foursquare (TKY)Precision@100.09830.14240.16630.15690.13220.1172
MRR@100.43180.46550.56450.51420.47310.4128
Table 3. Performance comparison on top-10 route recommendation.
Table 3. Performance comparison on top-10 route recommendation.
DatasetAlgorithmPrecisionMRR
Foursquare NYCWRMF0.13920.3712
WARP-MF0.12530.3517
SPRMF0.19270.5461
SG-FWARP0.24360.6814
Foursquare TKYWRMF0.09120.2310
WARP-MF0.08360.2150
SPRMF0.12610.3952
SG-FWARP0.16320.5147
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gu, J.; Song, C.; Ren, Z.; Lu, L.; Jiang, W.; Liu, M. Pedestrian Flow Prediction and Route Recommendation with Business Events. Sensors 2022, 22, 7478. https://doi.org/10.3390/s22197478

AMA Style

Gu J, Song C, Ren Z, Lu L, Jiang W, Liu M. Pedestrian Flow Prediction and Route Recommendation with Business Events. Sensors. 2022; 22(19):7478. https://doi.org/10.3390/s22197478

Chicago/Turabian Style

Gu, Jiqing, Chao Song, Zheng Ren, Li Lu, Wenjun Jiang, and Ming Liu. 2022. "Pedestrian Flow Prediction and Route Recommendation with Business Events" Sensors 22, no. 19: 7478. https://doi.org/10.3390/s22197478

APA Style

Gu, J., Song, C., Ren, Z., Lu, L., Jiang, W., & Liu, M. (2022). Pedestrian Flow Prediction and Route Recommendation with Business Events. Sensors, 22(19), 7478. https://doi.org/10.3390/s22197478

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop