Analysis of the Cycling Flow Between Origin and Destination for Dockless Shared Bicycles Based on Singular Value Decomposition

: Recently, an increasing number of cities have deployed bicycle-sharing systems to solve the ﬁrst / last mile connection problem, generating a large quantity of data. In this paper, singular value decomposition (SVD) was used to extract the main features of the cycling ﬂow from the origin and destination (OD) data of shared bicycles in Beijing. The results show that (1) pairs of OD ﬂow clusters can be derived from the pairs of vectors after SVD, and each pair of clusters represents a small part of an area with dockless shared bicycles; (2) the spatial clusters derived from the top vectors of SVD are highly coincident with the hot spot areas in the heatmap of shared bicycles; (3) approximately 30% of the study area accounts for nearly 80% of bike riding; (4) nearly 70% of the clustered area derived from the top 1000 vectors of SVD is associated with subway stations; and (5) the types of point of interest (POI) di ﬀ er between the origin area and destination area for the clustered area of the top 1000 vectors.


Introduction
The public bike-sharing (PBS) system has received increasing attention as a potential way to improve the first and last mile connections to other modes of transit and lessen the environmental impact of transportation activities [1][2][3].There are two main types of shared bicycle systems, namely docked and dockless bike.The docked bike-sharing system only allows users to ride bikes from one dock to another and the bikes should be locked in place.In recent years, dockless shared bicycles have become increasingly popular among urban residents due to their advantages of being unlocked by smart phones and parked anywhere within the service area.Regular spatial patterns may be concealed under the seemingly disordered movement of the bikes.Additionally, commuting regularities are usually linked to the spatial distribution or the transfer routing of bikes [4].
Many methods have been proposed to identify the regular pattern of bike movements, such as the fast-greedy algorithm [5], the Poisson mixture model, the Markov chain [6], the expectation maximization algorithm [7], and geographical business intelligence [8].Usage patterns and impact factors have been compared between a city's downtown and suburban district [9].These studies show that spatial analysis of the bicycle-sharing system is conducive to research on urban planning, urban traffic management, bicycle layout optimization, and urban population migration.However, there is limited research on the cycling flow between origin and destination (OD) pairs for dockless shared bicycles.Research on cycling flow could present the spatial distribution and movements of sharing bicycles on a smaller scale [10].In this paper, the OD flow is defined as the cycling flow between a pair of OD points.The analysis of the OD flows of taxi traces [11], mobile phones [12,13], and Integrated Circuit (IC) cards [14,15] can help to explore human mobility and location characteristics and provide necessary data for traffic planning, road construction, as well as traffic control and management [16,17].An analysis of cycling flow between OD pairs is required for the spatial distribution and movements for dockless shared bicycles.
Many studies have focused on finding clusters of origin and destination points by considering the attributes and spatial distributions of the OD points [18].Among them, spatial statistics-based methods such as Moran's I [19], hierarchical clustering methods such as K-nearest-neighbors [20], and density-based clustering methods including Density-Based Spatial Clustering of Applications with Noise (DBSCAN) [21] can be used for clustering origin and destination points.Some studies have proposed many methods for extracting OD clusters using long-distance OD data, such as taxi trips.Bicycles are not suitable for long-distance commuting due to the restrictions of the bicycle itself.Some methods include, for example, the simple line clustering method [22], based on the OD line to find spatial linkage.However, this research attempts to construct an OD matrix to specify the travel demand from origin points to destination points microscopically from the perspective of short-distance commuting.
Usually, an OD matrix is a large sparse matrix [23].A dimensionality reduction method, which is called feature extraction, is necessary on such a large matrix to simplify the data without losing too much information [24].Many methods are used to mine information from the OD matrix, such as affinity propagation, which is used to cluster the dimensionality-reduced matrix [25].The tensor decomposition model is also widely used to handle taxi OD data from multiple dimensions [26].The main purpose of the above operation is to simplify the dataset and extract the main feature from a matrix.These existing methods have good performance in OD clustering and spatial distribution analysis.In order to further analyze the characteristics of shared bicycle flow microscopically, a simple and fast algorithm called singular value decomposition (SVD) is widely used in data mining to easily extract information in pairs [27,28].
SVD, a factorization of a real or complex matrix, is commonly used for dimensionality reduction [25], and has been widely used in many applications, such as recommender systems [29], feature recognition [28], and machine fault diagnosis [30].These studies show that SVD can find the main components from matrix-shaped data.In particular, SVD has been used to extract features from the OD matrix of smart card records of Shenzhen, China [31].Currently, little research exists on applying SVD directly to OD data mining and extraction.Theoretically, SVD is suitable for discovering the main clusters of cycling flow between short-distance OD pairs from the singular values.
This study aimed to analyze the characteristics of cycling flow between OD pairs for dockless shared bicycles.SVD was applied to extract the features of the paired clusters from the OD matrix.The remainder of this paper is organized into the following sections.Section 2 introduces the research area and the data and preprocessing methods, followed by the SVD methodology and the results in Section 3. Section 4 analyzes the results of SVD and discusses the spatial features of cycling flows.Finally, Section 5 draws conclusions and presents future research prospects.

Research Area and Data
As the capital of China, Beijing took the lead in deploying bicycle-sharing projects in 2016.With the joint efforts of the government and shared bicycle companies, the bicycle-sharing system of Beijing began to develop in the first quarter of 2017.A large number of smart phone users prefer to use dockless shared bicycles in Beijing, China.In this study, the research area is the area within the 6th Ring Road of Beijing (Figure 1).In all, 3,104,089 rows of bicycle-sharing data were collected for Beijing from 10 May 2017 to 24 May 2017 (with missing data on 17 May 2017).The records came from Mobike, a leading bicyclesharing company in China.Some sample data are shown in Table 1.Each record is generated from a passenger ride, which contains a unique number to identify the record (order_id), a unique user number (user_id), a unique bicycle number (bike_id), a time stamp when the user unlocked the bicycle and started cycling (start_time), and the location coordinates where the user started (geohashed_start_loc) and stopped (geohashed_end_loc) cycling.A geo-hash decoding algorithm was applied to the last two columns to decompress the latitude and longitude coordinates from the encrypted string.Excluding some POI categories that were not related to bicycle travel, such as Auto Service, we finally collected more than 900,000 POIs in 13 major categories.Table 2 shows the names of the different POI categories and basic descriptions.In all, 3,104,089 rows of bicycle-sharing data were collected for Beijing from 10 May 2017 to 24 May 2017 (with missing data on 17 May 2017).The records came from Mobike, a leading bicycle-sharing company in China.Some sample data are shown in Table 1.Each record is generated from a passenger ride, which contains a unique number to identify the record (order_id), a unique user number (user_id), a unique bicycle number (bike_id), a time stamp when the user unlocked the bicycle and started cycling (start_time), and the location coordinates where the user started (geohashed_start_loc) and stopped (geohashed_end_loc) cycling.A geo-hash decoding algorithm was applied to the last two columns to decompress the latitude and longitude coordinates from the encrypted string.In addition, we collected point of interest (POI) information within the 6th Ring Road in Beijing.Excluding some POI categories that were not related to bicycle travel, such as Auto Service, we finally collected more than 900,000 POIs in 13 major categories.Table 2 shows the names of the different POI categories and basic descriptions.

OD Matrix Construction
The research area had to be rasterized to construct the OD matrix.The research area was divided using a certain size of fishnet, and the number of OD points of shared bicycles in each grid was counted to construct the OD matrix.To define the size of the fishnet, the network distances between each pair of OD points were calculated based on the road network in Beijing.Figure 2 shows the frequency distribution histogram of network distance.The results indicate that more than 87% of bicycle trips were longer than 500 meters.Considering the spatial scale and computational complexity, this study implements a latitude and longitude network with a side length of 0.005 degrees (approximately 430 meters in Beijing) as the rasterization unit.

OD Matrix Construction
The research area had to be rasterized to construct the OD matrix.The research area was divided using a certain size of fishnet, and the number of OD points of shared bicycles in each grid was counted to construct the OD matrix.To define the size of the fishnet, the network distances between each pair of OD points were calculated based on the road network in Beijing.Figure 2 shows the frequency distribution histogram of network distance.The results indicate that more than 87% of bicycle trips were longer than 500 meters.Considering the spatial scale and computational complexity, this study implements a latitude and longitude network with a side length of 0.005 degrees (approximately 430 meters in Beijing) as the rasterization unit.A total of 12,250 grids covered the extent of the research area after rasterization.The grids were numbered from bottom to top and from left to right.The average uses of shared bicycles on weekdays and the weekend were counted every ten minutes for each grid, as shown in Figure 3.The data trend on the weekend was quite different from that on weekdays, and different colors were used to distinguish the usage of shared bicycles on the weekend and weekdays.The plot shows three peaks on weekdays during which bicycle usage rises sharply.The morning rush hour lasts from 7 a.m. to 9 a.m. on weekdays.The weekday evening rush hour starts at 5 p.m., and after reaching its maximum at 6 p.m., the line falls continuously until midnight.A small peak at approximately 12 p.m. indicates a noon rush hour on weekdays.By contrast, no obvious peaks are observed on the weekend.A total of 12,250 grids covered the extent of the research area after rasterization.The grids were numbered from bottom to top and from left to right.The average uses of shared bicycles on weekdays and the weekend were counted every ten minutes for each grid, as shown in Figure 3.The data trend on the weekend was quite different from that on weekdays, and different colors were used to distinguish the usage of shared bicycles on the weekend and weekdays.The plot shows three peaks on weekdays during which bicycle usage rises sharply.The morning rush hour lasts from 7 a.m. to 9 a.m. on weekdays.The weekday evening rush hour starts at 5 p.m., and after reaching its maximum at 6 p.m., the line falls continuously until midnight.A small peak at approximately 12 p.m. indicates a noon rush hour on weekdays.By contrast, no obvious peaks are observed on the weekend.To study the travel characteristics of shared bicycles during different time periods, two OD matrices were constructed for the morning and evening rush hours on weekdays.The shape of the OD matrix was 12,250 × 12,250, where the longitudinal dimension represents each origin grid and the transverse is the destination grid (Figure 4).The elements M(m,n) in the matrix represent the number of bicycles that migrate from the grid O(m) to the grid D(n).

Singular Value Decomposition
In linear algebra, SVD is a factorization of a matrix and has been used extensively for dimensionality reduction in pattern recognition and information retrieval applications [32].SVD is a generalization of the eigendecomposition of a positive semidefinite normal matrix to any matrix via an extension of polar decomposition [33].Figure 5 is a schematic diagram of SVD.To study the travel characteristics of shared bicycles during different time periods, two OD matrices were constructed for the morning and evening rush hours on weekdays.The shape of the OD matrix was 12,250 × 12,250, where the longitudinal dimension represents each origin grid and the transverse is the destination grid (Figure 4).The elements M(m,n) in the matrix represent the number of bicycles that migrate from the grid O(m) to the grid D(n).To study the travel characteristics of shared bicycles during different time periods, two OD matrices were constructed for the morning and evening rush hours on weekdays.The shape of the OD matrix was 12,250 × 12,250, where the longitudinal dimension represents each origin grid and the transverse is the destination grid (Figure 4).The elements M(m,n) in the matrix represent the number of bicycles that migrate from the grid O(m) to the grid D(n).

Singular Value Decomposition
In linear algebra, SVD is a factorization of a matrix and has been used extensively for dimensionality reduction in pattern recognition and information retrieval applications [32].SVD is a generalization of the eigendecomposition of a positive semidefinite normal matrix to any matrix via an extension of polar decomposition [33].Figure 5 is a schematic diagram of SVD.

Singular Value Decomposition
In linear algebra, SVD is a factorization of a matrix and has been used extensively for dimensionality reduction in pattern recognition and information retrieval applications [32].SVD is a generalization of the eigendecomposition of a positive semidefinite normal matrix to any matrix via an extension of polar decomposition [33].Figure 5 is a schematic diagram of SVD.Suppose M is an m × n matrix composed of real numbers.Then, there exists a factorization, called the SVD of M, of the following form [34]: where U is an m × m unitary matrix, Σ is a diagonal m × n matrix with non-negative real numbers called singular values, V is an n × n unitary matrix, and  is the transpose of V.Each pair of column vectors  ⃗ and  ⃗ in matrices U and V are singular vectors corresponding to the diagonal entries and singular values σi in the diagonal matrix Σ, respectively.The singular values σi of the diagonal matrix Σ are listed in descending order.The larger the singular value σi, the more information can be provided by the value in the matrix Σ and the corresponding vectors.The cumulative proportion is calculated to measure the total cumulative contribution rate of the anterior components of σi.The cumulative proportion is calculated by the following formula [35]: where CP(k) describes the total accumulative contribution rate of the anterior k components and σi is the singular value in diagonal matrix Σ.

SVD Results
After applying SVD to the OD matrices, three matrices U, Σ, and V of size 12,250 × 12,250 are generated.The eigenvalues in the diagonal matrix Σ calculated by SVD are sorted in descending order.Figure 6 shows the singular values in Σ of the OD matrices obtained from morning rush hour and evening rush hour.The results show that the singular value decreases sharply and approaches zero as the ordinal rises.The cumulative proportion of the top 2000 singular values from the two OD matrices exceeded 85% (morning rush hour, 92.78%; evening rush hour, 87.95%).Generally, these top singular values and singular vectors can accurately represent the information of the original matrix.Each singular value in the diagonal matrix Σ corresponds to a pair of singular vectors in matrices U and V. Sorting by singular values, the top-ranked vectors contain more information.We visualized Suppose M is an m × n matrix composed of real numbers.Then, there exists a factorization, called the SVD of M, of the following form [34]: where U is an m × m unitary matrix, Σ is a diagonal m × n matrix with non-negative real numbers called singular values, V is an n × n unitary matrix, and V T is the transpose of V.Each pair of column vectors → u i and → v i in matrices U and V are singular vectors corresponding to the diagonal entries and singular values σ i in the diagonal matrix Σ, respectively.The singular values σ i of the diagonal matrix Σ are listed in descending order.The larger the singular value σ i , the more information can be provided by the value in the matrix Σ and the corresponding vectors.The cumulative proportion is calculated to measure the total cumulative contribution rate of the anterior components of σ i .The cumulative proportion is calculated by the following formula [35]: where CP(k) describes the total accumulative contribution rate of the anterior k components and σ i is the singular value in diagonal matrix Σ.

SVD Results
After applying SVD to the OD matrices, three matrices U, Σ, and V of size 12,250 × 12,250 are generated.The eigenvalues in the diagonal matrix Σ calculated by SVD are sorted in descending order.Figure 6 shows the singular values in Σ of the OD matrices obtained from morning rush hour and evening rush hour.The results show that the singular value decreases sharply and approaches zero as the ordinal rises.The cumulative proportion of the top 2000 singular values from the two OD matrices exceeded 85% (morning rush hour, 92.78%; evening rush hour, 87.95%).Generally, these top singular values and singular vectors can accurately represent the information of the original matrix.Suppose M is an m × n matrix composed of real numbers.Then, there exists a factorization, called the SVD of M, of the following form [34]: where U is an m × m unitary matrix, Σ is a diagonal m × n matrix with non-negative real numbers called singular values, V is an n × n unitary matrix, and  is the transpose of V.Each pair of column vectors  ⃗ and  ⃗ in matrices U and V are singular vectors corresponding to the diagonal entries and singular values σi in the diagonal matrix Σ, respectively.The singular values σi of the diagonal matrix Σ are listed in descending order.The larger the singular value σi, the more information can be provided by the value in the matrix Σ and the corresponding vectors.The cumulative proportion is calculated to measure the total cumulative contribution rate of the anterior components of σi.The cumulative proportion is calculated by the following formula [35]: where CP(k) describes the total accumulative contribution rate of the anterior k components and σi is the singular value in diagonal matrix Σ.

SVD Results
After applying SVD to the OD matrices, three matrices U, Σ, and V of size 12,250 × 12,250 are generated.The eigenvalues in the diagonal matrix Σ calculated by SVD are sorted in descending order.Figure 6 shows the singular values in Σ of the OD matrices obtained from morning rush hour and evening rush hour.The results show that the singular value decreases sharply and approaches zero as the ordinal rises.The cumulative proportion of the top 2000 singular values from the two OD matrices exceeded 85% (morning rush hour, 92.78%; evening rush hour, 87.95%).Generally, these top singular values and singular vectors can accurately represent the information of the original matrix.Each singular value in the diagonal matrix Σ corresponds to a pair of singular vectors in matrices U and V. Sorting by singular values, the top-ranked vectors contain more information.We visualized

Discussion of OD Flows
To analyze the spatial distribution of the above six vectors, we remapped the vector values to the spatial grid according to the grid number, as is shown in Figure 7.The results show that OD flow clusters can be found after setting zero as null.Each cluster on the map is derived from a twinning vector, and some close clusters may overlap.Specifically, the twelve clusters in Figure 8 are associated with the top six pairs of vectors in Figure 7.These pairs of clusters can divide the space into small pieces, and each pair can represent an OD flow from cluster  ⃗ to cluster  ⃗.Actually, each cluster marked by a vector symbol in Figure 8 is the reflection of a vector in Figure 7.

Discussion of OD Flows
To analyze the spatial distribution of the above six vectors, we remapped the vector values to the spatial grid according to the grid number, as is shown in Figure 7.The results show that OD flow clusters can be found after setting zero as null.Each cluster on the map is derived from a twinning vector, and some close clusters may overlap.Specifically, the twelve clusters in Figure 8 are associated with the top six pairs of vectors in Figure 7.These pairs of clusters can divide the space into small pieces, and each pair can represent an OD flow from cluster → u i to cluster → v i .Actually, each cluster marked by a vector symbol in Figure 8 is the reflection of a vector in Figure 7.
Each pair of singular vectors represents a small part of the area in Beijing that accounts for a large portion of shared bicycle usage.Figure 9 shows two pairs of singular vectors during the morning rush hour.As mentioned previously, each vector represents a part of the area.In Figure 9a, the red area with red stripes is an origin area in the morning and covers a large part of a residential area.The area with blue stripes is the destination.This chart reflects that people who live near subway stations ride shared bicycles to the subway stations.In Figure 9b, the origin areas cover two subway stations, while the destination areas are composed of the central business district (CBD) of Beijing.This chart indicates that people might be cycling from subway stations to the CBD.Additional visualization results show that the migration relationship revolves around the subway station, especially for top-ranked vectors.
The weight of each vector is related to the order of the singular values, so the top vectors can extract more important information from the OD matrix.In this case, the top singular vectors reflect the hot spots of shared bicycle travel.Figure 10 shows an overlay of the spatial reflection of the top 50 singular vectors and a heatmap of shared bicycle OD points obtained using the kernel density during the morning rush hour.The results show that the origin areas identified by SVD and the heatmap overlap substantially; thus, SVD can extract the main OD hot spot areas.
The top 20, 50, 100, 200, 500, and 1000 singular vectors are overlaid from Figure 10 with the heatmap of the origin points for the morning rush hour in Figure 11.The results show that the greater the number of singular vectors added to the heatmap, the larger the areas covered.The covered areas are highly coincident with the hot spots of the heatmap.Similarly, the destinations during the morning rush hour as well as the origins and destinations during the afternoon rush hour have analogous properties.Figure 12 shows that as more singular vectors are added, the coverage rates of both the OD areas increase.This upward trend is consistent with that of the cumulative proportion, which is calculated by Equation ( 2).The top 1000 vectors are used here.As the cumulative proportion approaches 80%, the area ratio of the map covered by the OD areas is approximately 30%.The results show that nearly 80% of shared bicycle travel is distributed in approximately 30% of the area within the 6th Ring Road of Beijing.
Among the top 1000 SVD vectors, nearly 70% of the clustered areas are related to subway stations, which indicates that people often ride shared bicycles between subway stations and other places as part of a short-distance commute.This fact is in line with the original intention of using shared bicycles to solve the first/last mile transportation problem.The results verify our previous research on the effects of free-floating shared bicycles on urban public transportation [36].Each pair of singular vectors represents a small part of the area in Beijing that accounts for a large portion of shared bicycle usage.Figure 9 shows two pairs of singular vectors during the morning rush hour.As mentioned previously, each vector represents a part of the area.In Figure 9a, the red area with red stripes is an origin area in the morning and covers a large part of a residential area.The area with blue stripes is the destination.This chart reflects that people who live near subway The weight of each vector is related to the order of the singular values, so the top vectors can extract more important information from the OD matrix.In this case, the top singular vectors reflect the hot spots of shared bicycle travel.Figure 10 shows an overlay of the spatial reflection of the top 50 singular vectors and a heatmap of shared bicycle OD points obtained using the kernel density during the morning rush hour.The results show that the origin areas identified by SVD and the heatmap overlap substantially; thus, SVD can extract the main OD hot spot areas.The top 20, 50, 100, 200, 500, and 1000 singular vectors are overlaid from Figure 10 with the heatmap of the origin points for the morning rush hour in Figure 11.The results show that the greater the number of singular vectors added to the heatmap, the larger the areas covered.The covered areas are highly coincident with the hot spots of the heatmap.Similarly, the destinations during the morning rush hour as well as the origins and destinations during the afternoon rush hour have both the OD areas increase.This upward trend is consistent with that of the cumulative proportion, which is calculated by Equation ( 2).The top 1000 vectors are used here.As the cumulative proportion approaches 80%, the area ratio of the map covered by the OD areas is approximately 30%.The results show that nearly 80% of shared bicycle travel is distributed in approximately 30% of the area within the 6th Ring Road of Beijing.Among the top 1000 SVD vectors, nearly 70% of the clustered areas are related to subway stations, which indicates that people often ride shared bicycles between subway stations and other places as part of a short-distance commute.This fact is in line with the original intention of using shared bicycles to solve the first/last mile transportation problem.The results verify our previous research on the effects of free-floating shared bicycles on urban public transportation [36].
Moreover, the proportions of various POIs in the areas covered by the top 1000 vectors during the morning rush hour and evening rush hour were calculated, as shown in Figure 13 The results indicate that people who are related to those places have a contrasting migration trend during the morning and evening rush hours, which is consistent with the tidal characteristics of shared bicycles reported in another article [36].The results indicate that people who are related to those places have a contrasting migration trend during the morning and evening rush hours, which is consistent with the tidal characteristics of shared bicycles reported in another article [36].

Conclusions
In this paper, SVD was applied to extract the main components from the OD matrix of dockless shared bicycles in Beijing.SVD can be used to find clusters of OD areas for bicycle riding in pairs.After SVD, one pair of clusters can be derived from a vector on both the OD map, and each pair of vectors represents a small part of the area in Beijing that accounts for a large portion of bicycle usage.The clusters derived from the top vectors are highly coincident with the hot spot areas in the heatmap of shared bicycle OD points.More clusters derived from the top singular vectors correspond to a higher area coverage rate of shared bicycle travel.The results show that nearly 80% of shared bicycle travel is distributed in approximately 30% of the study area in Beijing.In addition, more POIs related to Food and Beverages, Shopping, Daily Life Service, and Enterprises are located in the areas covered by the top 1000 vectors during morning rush hour.The proportion of various POIs differs between the OD areas of shared bicycle travel.
In further research, the OD matrices of shared bicycle travel could be constructed with different dimensions, including time, riding distance, or duration.Other matrix decomposition methods could be applied to explore potential space-time regulations.

Conclusions
In this paper, SVD was applied to extract the main components from the OD matrix of dockless shared bicycles in Beijing.SVD can be used to find clusters of OD areas for bicycle riding in pairs.After SVD, one pair of clusters can be derived from a vector on both the OD map, and each pair of vectors represents a small part of the area in Beijing that accounts for a large portion of bicycle usage.The clusters derived from the top vectors are highly coincident with the hot spot areas in the heatmap of shared bicycle OD points.More clusters derived from the top singular vectors correspond to a higher area coverage rate of shared bicycle travel.The results show that nearly 80% of shared bicycle travel is distributed in approximately 30% of the study area in Beijing.In addition, more POIs related to Food and Beverages, Shopping, Daily Life Service, and Enterprises are located in the areas covered by the top 1000 vectors during morning rush hour.The proportion of various POIs differs between the OD areas of shared bicycle travel.
In further research, the OD matrices of shared bicycle travel could be constructed with different dimensions, including time, riding distance, or duration.Other matrix decomposition methods could be applied to explore potential space-time regulations.

Figure 3 .
Figure 3. Average uses of shared bicycles on weekdays and the weekend.

Figure 3 .
Figure 3. Average uses of shared bicycles on weekdays and the weekend.

14 Figure 3 .
Figure 3. Average uses of shared bicycles on weekdays and the weekend.

Figure 5 .
Figure 5. Visualization of the matrices in singular value decomposition (SVD).

Figure 6 .
Figure 6.Singular values of both morning and evening rush hours.

Figure 5 .
Figure 5. Visualization of the matrices in singular value decomposition (SVD).

Figure 6 .
Figure 6.Singular values of both morning and evening rush hours.

Figure 6 .
Figure 6.Singular values of both morning and evening rush hours.

Figure 7 .
Figure 7. Sample vectors from the SVD results during the morning rush hour.

Figure 7 .
Figure 7. Sample vectors from the SVD results during the morning rush hour.

14 Figure 8 .
Figure 8. Spatialization of the top six pairs of vectors.

Figure 8 . 14 Figure 9 .
Figure 8. Spatialization of the top six pairs of vectors.

Figure 10 .
Figure 10.Overlay of the top 50 singular vectors and heatmap of origin and destination (OD) points.

Figure 11 .
Figure 11.Overlay of the top 20-1000 singular vectors and heatmap of origin points.Figure 11.Overlay of the top 20-1000 singular vectors and heatmap of origin points.

Figure 11 . 14 Figure 12 .
Figure 11.Overlay of the top 20-1000 singular vectors and heatmap of origin points.Figure 11.Overlay of the top 20-1000 singular vectors and heatmap of origin points.ISPRS Int.J. Geo-Inf.2019, 8, x FOR PEER REVIEW 11 of 14 . The results show that the most common POIs in this area are Food and Beverages (No. 1), Shopping (No. 2), Daily Life Service (No. 3), and Enterprises (No. 13).In addition, the POIs of Daily Life Service (No. 3), Accommodation Service (No. 6), and Commercial House (No.8) have a higher proportion in the origin areas than in the destination areas during the morning rush hour.By contrast, the proportion of POIs related to Enterprises (No. 13) in the destination areas is larger than the value in the origin areas during the morning rush hour.That is, people tend to unlock bicycles from the areas related to Daily Life Service (No. 3), Accommodation Service (No. 6).and Commercial House (No.8) to start riding to the areas of Enterprises (No. 13) in the morning.Conversely, the trend of POI between the origin and destination areas in the evening is opposite to the trend in the morning in the areas related to Daily Life Service (No. 3), Accommodation Service (No. 6), Commercial House (No.8), and Enterprises (No. 13).

Figure 12 .
Figure 12.Spatial coverage rate and cumulative proportion for the top 1000 singular vectors.Moreover, the proportions of various POIs in the areas covered by the top 1000 vectors during the morning rush hour and evening rush hour were calculated, as shown in Figure 13.The results show that the most common POIs in this area are Food and Beverages (No. 1), Shopping (No. 2), Daily Life Service (No. 3), and Enterprises (No. 13).In addition, the POIs of Daily Life Service (No. 3), Accommodation Service (No. 6), and Commercial House (No.8) have a higher proportion in the origin areas than in the destination areas during the morning rush hour.By contrast, the proportion of POIs related to Enterprises (No. 13) in the destination areas is larger than the value in the origin areas during the morning rush hour.That is, people tend to unlock bicycles from the areas related to Daily Life Service (No. 3), Accommodation Service (No. 6).and Commercial House (No.8) to start riding to the areas of Enterprises (No. 13) in the morning.Conversely, the trend of POI between the origin and destination areas in the evening is opposite to the trend in the morning in the areas related to Daily Life Service (No. 3), Accommodation Service (No. 6), Commercial House (No.8), and Enterprises (No. 13).The results indicate that people who are related to those places have a contrasting migration trend during the morning and evening rush hours, which is consistent with the tidal characteristics of shared bicycles reported in another article[36].

Table 2 .
Point of interest (POI) types and descriptions.

Table 2 .
Point of interest (POI) types and descriptions.