Mining the Spatial-Temporal Characteristics of Population Flow and Material Flow Using Big Data

Population ﬂow and material ﬂow are the core elements of the “space of ﬂows.” Studying the temporal and spatial distribution characteristics of these two types of ﬂows (TToF) can provide us a greater understanding of the research area. Most of the existing works related to TToF only focus on exploring the diﬀerence characteristics of one of the members of TToF in a certain time or space scale in the research area. Diﬀerent from these related works, the spatial-temporal characteristics of the population ﬂow and material ﬂow in Taiwan Province and the spatial-temporal autocorrelation of Taiwan’s expressway network are explored by means of multimembership and layer-by-layer reﬁnement. The research work carried out in this paper includes the following: (1) studying the diﬀerentiated characteristics of the TToF in diﬀerent time units; (2) studying the spatial diﬀerences among each type of the TToF under diﬀerent scales; (3) dividing both the population ﬂow and material ﬂow into two subtypes and then analyzing the temporal variation characteristics of the four subtypes of ﬂows; and (4) studying the global and local spatial-temporal autocorrelation of Taiwan’s expressway network. The results show the following. (1) The spatial-temporal diﬀerentiation characteristics of the TToF are obvious in diﬀerent time units and on diﬀerent scales. (2) The contribution of the population ﬂow to the TToF in ﬂow quantities is far greater than that of the material ﬂow. (3) The population ﬂow and material ﬂow are dominated by the “minority population ﬂow” and “small-scale material ﬂow,” respectively. (4) Meanwhile, in Taiwan’s expressway network, there is a signiﬁcant spatial-temporal positive correlation mainly reﬂected in the spatial ﬁrst-order adjacent road sections.


Introduction
With the rapid development of information technology, the influence of a fixed geographical location on production and organization has gradually weakened. Power and function are organized in the flows. e dominant form of space is no longer the traditional space ("space of places") but rather is a new space ("space of flows"). is new space transcends the boundaries of time and space, playing a guiding role in the development of regions' spatial patterns in the future. With the help of the interaction of various "flows" in the space, it extends from its interior to many fields of human production and life. e core elements of the "space of flows" include population flows, material flows, capital flows, and information flows. Among these flows, information flows, material flows, and population flows are often mentioned and discussed. However, the data associated with the information flows are very difficult to obtain because it usually involves the private information of units or individuals. According to the literature we collected, we determine that the current studies on flows are more focused on population flows and material flows.
Population mobility is an extremely complex research topic. Many researchers can only choose some more specific subtopics for research, such as studying the impact of population mobility on regional economic development [1][2][3], exploring the relationship between population mobility and urbanization development [4][5][6], and mining the spatial pattern and formation mechanism of population mobility [7][8][9][10]. Based on the interprovincial migration data during the period 1985-2010, literature [7] used bicomponent trend mapping technique with every five years as a time interval to analyze the intensity and change pattern of in-migration and out-migration from 1985 to 2010. Literature [8] analyzed the spatial-temporal characteristics and influencing factors of population changes in Beijing-Tianjin-Hebei (BTH) region from the two aspects of population distribution patterns and population flow in long and short periods. eir work enlightens us that in order to understand a problem or phenomenon more thoroughly, we can analyze it from many different perspectives. e process of population mobility is often accompanied by the flow of information, material, capital, and technology. e mobility of a population among different regions leads to the dynamic adjustment of relevant resources, such as the labor force, talent, technology, and so on [11]. e impact of this dynamic adjustment is often complex. For example, a labor force inflow greatly supports and speeds up the infrastructure construction in the destination areas. In contrast, a labor force outflow also causes a large number of old people and children to be left behind and brings a series of social problems. As another example, the mobility of talents from backward areas to developed areas will further increase the gap in the comprehensive strength and competitiveness of two types of areas. In contrast, the gap between the two types of areas may gradually become narrow. Furthermore, long-term or regular large-scale population migration will not only affect the decision making of relevant government departments but also stimulate the construction of the infrastructure (such as housing, medical care, education, and so on) closely related to the floating population in the relevant regions [12]. e material flow is the most important flow in an economic system. It influences regional economic growth by affecting both supply and demand. It is not only conducive to the deepening of the regional division of labor and specialized production and the effective input of elements, but it can also expand output demand [13]. Currently, management scholars are committed to studying the patterns of material flows [14,15] and discussing the location planning of material flow centers [16,17]. Economists are committed to mining the relationship between the material flow and economic development [18,19]. From the geographical perspective, the related research on material flows mainly focuses on the discovery of the spatial patterns and the evolutionary development of a material flow enterprise or the regional logistics chain [20][21][22][23][24]. e method of data envelopment analysis (DEA) and exploring statistics has been used by literature [23] to analyze the whole evolution trend, structure, and spatial pattern's evolutionary process and its evolutionary mechanism of Yangtze River delta city logistics development efficiency from 2002 to 2011 based on the input-output index of city logistics. Literature [24] takes the major 22 ports along the Yangtze River as the research object and develops a comprehensive evaluation index system of port logistics development. en, it analyzes the spatiotemporal evolution characteristics and mechanism of port logistics development along the Yangtze River, by using the methods of Gini coefficient and standard deviational ellipse. eir work has extended the application of some classical correlation analysis methods.
To sum up, population mobility is a complex and important research topic. Its process is often accompanied by the flow of information, material, capital, and technology. e analysis of population mobility is of great significance to the development of economy and the formulation of policies. Most of the existing related literatures study population mobility from a certain dimension or local subproblem with specific research purposes but lack multidimensional perspective of exploration and research. Material flow is an important manifestation of economic activities. e study of material flow is not only the need of theoretical research on industrial economy but also the practical need of the development of logistics industry. At present, most of the literatures focus on the research of logistics mode and logistics technology. e research of material flow needs to be further expanded and deepened. is paper is committed to providing multidimensional analysis methods for the study of population flow and material flow. Taiwan Province and its electronic toll collection (ETC) data are taken as the research area and fundamental experiment data, respectively.
Taiwan Province, with its superior geographical position, is located at the southeast gateway of China as an important channel to the Pacific Ocean. e study of the population flow and material flow in Taiwan Province can help us to better understand the province. At present, the main purpose of relevant studies on the population flow in Taiwan Province is to determine the long-term migration trend of the population on and between the islands [25][26][27], while the research on the material flow focuses on its logistics policy and competitiveness [28][29][30].
Different from the starting point of the existing related work [25][26][27][28][29][30], this paper focuses on the research of the spatial-temporal differentiation characteristics of the population flow and material flow (referred to "two types of flows, TToF") in Taiwan Province. In addition, different from those literatures which analyze the spatial-temporal differentiation characteristics of TToF in a specific geographical area and only in a certain time scale and space scale, we adopted the multimembership and layer-by-layer refinement methods for conducting our analysis work. e research work carried out in this paper mainly includes the following: (1) studying the differentiation characteristics of the TToF in different time units (hourly, daily, and monthly), (2) studying the spatial differences of each type of the TToF under different scales (detection sites, road sections, counties, and cities), and (3) studying the global and local spatial-temporal autocorrelation of Taiwan expressway network. Our experiments and analysis are mainly based on the ETC data of Taiwan Province in 2017. Our work would help to provide a better understanding of the travel habits of 2 Mathematical Problems in Engineering the people and the mode of the material flow in Taiwan Province.

Experimental Data and Methods
In this paper, we take Taiwan Province as the research area and the ETC data of its expressway as the fundamental experimental data. en, the spatial-temporal characteristics of the population flow and material flow in Taiwan Province and the spatial-temporal autocorrelation of Taiwan's expressway network are explored by means of multimembership and layer-by-layer refinement.  [31]. Figure 1 shows the administrative map of Taiwan Province. e experimental data used in this paper are the ETC data of Taiwan expressway in 2017. e data cover 16 counties and cities in Taiwan except Chiayi City, Hualien County, and Taitung County.

Experimental Data.
e ETC data of the expressway in Taiwan Province are our fundamental experimental data. e original data are recorded by the detection sites on the expressway and can be obtained from the official website. Considering the research purpose of this paper, it is necessary to further process the original data.

Data Sources.
At present, there are 331 detection sites along Taiwan expressway. ey divide Taiwan expressway into 481 detection areas. e location of each detection site and the number of detection sites in each county and city are shown in Figure 2 and Table 1, respectively.
In this paper, we take the M05A files from the ETC data of Taiwan expressway in 2017 as the fundamental experimental data. e data in this file are recorded every 5 minutes, and thus 253,017,045 records are obtained. Each record contains six fields, such as TimeInterval (report generation time), GantryFrom (no. of upstream detection sites), GantryTo (no. of downstream detection sites), VehicleType (vehicle's type), SpaceMeanSpeed (average speed), and VehicleCount (traffic volume for each vehicle type in the current time interval). e vehicle types are as follows: 31 for minibuses, 32 for pickup trucks, 41 for buses, 42 for trucks, and 5 for coupling vehicles.

Data Preprocessing.
To find the clear spatial-temporal differences of the population flow and material flow in Taiwan Province, we first obtain the latitude and longitude information of the 331 detection sites and import the data in the relevant record tables into an SQL Server 2008 database after merging them on a monthly basis. After that, extra fields such as "hour, week, and county/city of starting detection site and ending detection site" are added to each record. e detailed steps of data preprocessing are shown in Table 2. e times generated by the records can be used to analyze the temporal differences of the traffic volume, while the spatial differences in the traffic volume can be analyzed by combining the longitude and latitude information of detecting sites with the administrative divisions of Taiwan Province.
After analyzing the spatial-temporal differences of the traffic volume, this paper further studies the spatial-temporal autocorrelation of Taiwan's expressway network. Considering that the overall distribution of the phenomenon is described by the global analysis index, it is not only difficult to detect the correlation patterns of different regions, but it may also conceal some local heterogeneity by using only one single value to reflect the overall situation. erefore, this paper analyzes both the global and the local spatial-temporal autocorrelation of Taiwan's expressway network. e multimembership and layer-by-layer refinement analysis method is adopted to conduct our research work. Our research framework is shown in Figure 3. Seen from Figure 3, the analysis idea of multimembership and layer-by-layer refinement mainly embodies the following aspects. (1) Classifying the experimental data into two parts according to the actual purpose of vehicles, one for the analysis on population flow and the other for material flow.
(2) Both of the population flow and material flow are divided into two subtypes, and then four kinds of flows are obtained.
(3) For the four kinds of flows, we analyzed their differences in the temporal characteristics in different time units (hourly, daily, and monthly). (4) For the TToF, we studied their spatial differences under different scales (detection sites, road sections, counties, and cities). (5) For the TToF, we explored not only their global spatial-temporal autocorrelation but also local spatial-temporal autocorrelation.

Spatial and Temporal Difference Analysis Method.
To make the analysis results more accurate and make the best use of the experimental data, according to the actual purposes of vehicles, the records corresponding to vehicle types 31 and 41 in the experimental data are regarded as population flow data while the records corresponding to vehicle types 5, 32, and 42 are regarded as material flow data. In addition, to achieve a better understanding of Taiwanese travel habits and Mathematical Problems in Engineering the material flow patterns of Taiwan, we further subdivide the TToF into four subtypes: (1) considering the records corresponding to vehicle type 31 as minority population flow; (2) considering the records corresponding to vehicle type 41 as mass population flow data; (3) considering the records corresponding to vehicle type 32 as small-scale material flow data; and (4) considering the records corresponding to vehicle types 5 and 42 as large-scale material flow data.
We seek to determine the spatial-temporal differentiation characteristics of the TToF in Taiwan Province. By adopting the means of multimembership and layer-by-layer refinement methods, we not only study the temporal differences of different types of flows in different units (monthly, daily, and hourly) but also study the spatial differences of different types of flows in different units (detection sites, road sections, and counties or cities).

Spatial-Temporal Autocorrelation Analysis Method.
In this paper, the ST-ACF (spatial-temporal autocorrelation function) proposed by Pfeifer et al. [32,33] and the CCF (cross correlation function) proposed by Box and Jenkins [34] were used to analyze the global spatial-temporal autocorrelation and local spatial-temporal autocorrelation of Taiwan's expressway network, respectively. e ST-ACF can be used to measure the correlation between samples with an arbitrary spatial delay and time delay for the two aspects of space and time, respectively, while the CCF can be used to measure the correlation between two regional samples over a specified time delay. Some scholars have proved that CCF can be used to measure the spatial-temporal autocorrelation among road network sections [35]. e calculation formula for ST-ACF [32] is as follows: where ρ h0 (s) is the spatial-temporal autocorrelation coefficient; h and s stand for the space delay and time delay, respectively; r h (s) stands for the sample covariance with a spatial delay of h and a time delay of s; σ h (0) represents the sample variance with a spatial delay of h and a time delay of 0; W h represents the spatial weight matrix whose spatial delay is h (when h � 0, this weight matrix is an identity matrix); and Z(t) and Z(t + s) stand for the specific property values of research unit Z at times t and (t + s), respectively. In our experiment, Z(t) stands for the traffic Mathematical Problems in Engineering volume at time t of road section Z and s � 5 since the ETC system generates a record every five minutes. Let X and Y be two time series. en, their spatialtemporal correlation coefficients under the specified time delay s and based on the CCF function [34] can be calculated as follows: where X and Y are, respectively, the average values of X and Y and σ X and σ Y are the variances of the two respective sequences.
To calculate the values of the ST-ACF and CCF, the spatial weight matrix needs to be substituted into the formulas. In this paper, the spatial first-order adjacency matrix and spatial second-order adjacency matrix of each road section are used as the weight matrix following the practice of literature [36]. e first-order and second-order adjacency relations of each section in the road network are defined as follows. According to the above definition, the first-order and second-order adjacency matrices of the road network shown in Figure 4 en, the first-order and second-order adjacency matrices of Taiwan's expressway network which consists of 331 detection sites and 481 road sections can be generated. e two matrices will be denoted as W 1 and W 2 , respectively.

Analysis and Discussion
Different kinds of flows are classified by the vehicle type in each record. By this way, two broad types (population flow and material flow) and four subtypes (minority population flow, mass population flow, small-scale material flow, and large-scale material flow) are obtained. en, both the differentiated characteristics of each kind of flow in different time units and the spatial differences among each type of the TToF under different scales are analyzed. Moreover, we also discussed the global and local spatial-temporal autocorrelation of Taiwan's expressway network.

Analysis of the Differences in the Temporal
Characteristics of the Total Circulation of the TToF. As mentioned above, the differentiated characteristics of the TToF are analyzed in different time units (monthly, daily, and hourly).

Differences in the Monthly Characteristics of the TToF.
Since each data table is built monthly during data preprocessing, the total circulation of the TToF in each month can be obtained by querying the statistics based on the vehicle types. For the different days of each month, in order to more reasonably reflect the difference in the monthly flows, we first calculate the daily average flow of each month and then use the min-max standardization method to normalize the relevant data.
e results are shown in Figure 5. Figure 5 shows the following. (1) August is the most active period of Taiwan's material flow, and its tide period appears in June. According to the data provided by the official website of http://www.tianqihoubao.com/lishi, in June 2017, there are more than 20 rainy days on average in the major counties and cities of Taiwan Province, and there are even several days of heavy rain and rainstorms. Bad weather affected the traffic volume on the expressway, which has greatly reduced the material flow. (2) In January and February, the population flow is the most active. e main reason is that there are many important festivals such as New Year's Day, the Spring Festival, Valentine's Day, the Lantern Festival, and Peace Memorial Day in January and February.
ey are bound to stimulate the growth of the population flow. Influenced by bad weather, June is also a low population flow period. (3) Regarding the overall change trends, the change trend of the monthly mean flow of material flow is much gentler than that of the population flow. (4) e contribution of the population flow to total circulation is far greater than that of the material flow, which shows that its monthly change trend is very consistent with that of the total circulation of the TToF. e correlation between the blue curve and green curve is approximately 98.74%, while the correlation between the yellow curve and green curve is only approximately 34.89%.

Differences in the Daily Characteristics of the TToF.
By querying the statistics according to the new field "week" added in the preprocessing process and the vehicle types, the Process the ETC data of the expressway in Taiwan Province (the detailed process is shown in Table 2  Step 1 Obtain the longitudes and latitudes of the 331 detection sites and the counties and cities where they are located in from Google maps.
Step 2 For the convenience of the classification and statistics, use the command "copy * .csv reslutSave.csv" to merge the record data in the M05A file by month to form twelve csv files.
Step 3 Import the twelve CSV files into an SQL server database to obtain the twelve record tables.
Step 4 Additionally, for the convenience of the classification and statistics, add fields such as [hour, week, and administrative region where the origin and destination sites are located] to each record. e values of these added fields are set according to the information from Step1.
Step  Figure 6. e following can be seen from Figure 6. (1) On weekdays, the daily material flow shows an overall upward trend, and the peak appears on Friday. (2) On weekdays, the daily population flow decreases from Monday to ursday, but increases significantly on Friday. (3) e population flow on weekends is significantly higher than that on weekdays. One reason for this phenomenon is that Taiwan's infrastructure, such as its roads and corresponding services, is perfect. In contrast, the material flow drops sharply on weekends. Selectively avoiding the peak period of public travel can save transportation costs for logistics companies. (4) From the perspective of the total TToF, the total flow on weekends is much higher than that on weekdays because the contribution of the population flow to the total flow is much greater than that of the material flow. e percentage of the population flow to the total flow on weekdays is approximately 70.10% while that on weekends is approximately 81.08%.

Differences in the Hourly Characteristics of the TToF.
By querying the statistics according to the new field "hour" added in the preprocessing process and the vehicle types, the total circulation volume of the TToF on different hours of a day can be obtained. After normalizing the total population flow, the total material flow, and the sum of the two in each hour, the hourly change trend of the TToF is obtained, as shown in Figure 7. e following can be seen from Figure 7. (1) In a day, the most active periods of the material flow are 9:00-11:00 and 14:00-15:00 with peak points at 10:00 and 15:00, respectively. (2) e most active time periods of the population flow are 8:00-11:00 and 16:00-18:00 with the peak points at  4) On the whole, Taiwan's population flow and material flow increase rapidly starting at 6: 00 and reach the first peak at 10:00. ere is a slight decline from 10:00 to 12:00, and then it rises again. e second and third peaks of the TToF are at 15:00 and 17:00, respectively. After 17:00, the both components of the TToF rapidly decrease. We can find that the hourly change trends of the TToF coincide with Taiwanese work and living habits.

Analysis of the Differences in the Temporal Characteristics of the Minority Population Flow and Mass Population Flow.
To study Taiwanese travel habits more deeply and make full use of the data in each field of the original data, we further divide the population flow into the "minority population flow" (vehicle type 31) and the "mass population flow" (vehicle type 41) according to vehicle types in the original data. en, we analyze the time characteristics of the two kinds of population flows and find the differences between them. e processing methods of the two kinds of population flow data are consistent with those of the TToF data in the first subsection of this section.

Differences in the Monthly Characteristics of the Two Kinds of Population Flows.
e monthly mean flow and the monthly change curves of the two kinds of population flows are normalized and shown in Figure 8. e following can be seen from Figure 8. (1) January to February and July to August are the peak periods of the minority population flow with the troughs in June and November. (2) In contrast to the minority population flow, January to February and July to August are the troughs of the mass population flow, March to May and October to December are the active periods, and the highest peak time appears in November. e main reason for this phenomenon is that January to February and July to August are the months with the lowest and highest average temperatures in Taiwan, respectively. During these time periods, people tend to travel by car. On the whole, the monthly mean change trends of the mass population flow is gentler than that of the minority population flow. e contribution of the minority population flow to the total population flow is far greater than that of the mass population flow.

Differences in the Daily Characteristics of the Two Kinds of Population Flows.
e daily change trends of the mean minority population flow and mass population flow in a week are shown in Figure 9.    Mathematical Problems in Engineering e following is shown in Figure 9. (1) e minority population flow shows a slight downward trend from Monday to ursday, and it gradually increases in the next three days, reaching its peak on Sunday. (2) e mass population flow shows an overall upward trend from Monday to Saturday and has peaks on Saturday and declines slightly on Sunday. Combining (1) and (2), on weekdays, Taiwanese prefer to take public transport. It also reflects Taiwan's perfect transportation infrastructure and services. In addition, on the whole, the circulation volume of these two kinds of population flows on the weekends is higher than that on weekdays.

Differences in the Hourly Characteristics of the Two
Kinds of Population Flows. Figure 10 shows the hourly change trends of the mean minority population flow and mass population flow. Figure 10 shows us the following. (1) e minority population flow and mass population flow in each hour show double peak characteristics. (2) e two peaks of the two kinds of population flows are not the same. e first peak of the mass population flow appears at 8:00 while the first peak of minority population flow appears at 10:00. 17:00 is the other peak of the two kinds of population flows.

Analysis of the Differences in the Temporal Characteristics of the Small-Scale Material Flow and Large-Scale Material
Flow. Correspondingly, in order to further understand the pattern of the material flow in Taiwan Province, we divide the material flow into the small-scale material flow (vehicle type 32) and large-scale material flow (vehicle types 5 and 42) according to the vehicle types in the original data. en, the temporal characteristics of the two kinds of material flows are analyzed to find the differences between them. e data preprocessing of these two kinds of material flows is also consistent with that of the TToF in the first subsection of this section.

Differences in the Monthly Characteristics of the Two Kinds of Material Flows.
e monthly mean flows and the monthly change curves of the two kinds of material flows are normalized and shown in Figure 11. e statistical results show the following ( Figure 11). (1) March, June, August, and November are the active periods of the large-scale material flow while its troughs are from January to February. (2) e active periods of the small-scale material flow appear in January and August, and June is its trough. We think one of the main reasons for (1) and (2) is that the discount season in Taiwan varies with different products. (a) It is a good   time to buy shoes or clothes in January, February, July, and August. e transportation of these goods is more suitable for small-scale flows. (b) March, April, and September are the offseasons for shopping and the maximum discount for goods may be 20%. October to December is the peak of the discount season, and most businesses can offer a discount like "give 100 for every 1000." During this period, people's willingness to buy all kinds of goods (especially large commodities) is significantly higher than that of other months. Buying these large commodities will naturally result in a larger large-scale material flow. In general, the monthly mean flow change trends of the large-scale material flow and small-scale material flow are different. e varying amplitude of the small-scale material flow is relatively small. e contribution of the small-scale material flow to the total material flow is greater than that of the large-scale material flow.

Differences in the Daily Characteristics of the Two Kinds of Material Flows.
e daily change trends of the mean flows of the two kinds of material flows are shown in Figure 12. Figure 12 reveals the following. (1) From Monday to Friday, the large-scale material flow shows an upward trend, and the growth rate from Monday to Tuesday is obvious.
ursday is its most active period. (2) e small-scale material flow shows a slight downward trend from Monday to Wednesday with a significant increase from ursday to Friday. (3) e peaks of the two kinds of material flows are staggered. (4) On weekends, the small-scale material flow and the large-scale material flow are significantly less than those on weekdays. Especially, the large-scale material flow is rapidly declining on weekends.

Differences in the Hourly Characteristics of the Two Kinds of Material Flows.
e hourly change trends of the two kinds of material flows are shown in Figure 13.
(1) e large-scale material flow has three characteristic peaks with peak points at 6:00, 10:00-11:00, and 14:00. (2) e small-scale material flow shows the characteristic of double peaks. e peak points appear at 9:00-10:00 and 16: 00-17:00, which are staggered with the peaks of the largescale material flow. (3) On the whole, 6:00-17:00 is the active period for the two kinds of material flows and 00:00-5:00 and 18:00-23:00 are the relatively silent periods for material flows over the course of a whole day.

Analysis of the Differences in the Spatial Characteristics.
In this section, we computed and found the top 10 road sections for each kind of flow. en, the spatial distribution characteristics of those top 10 road sections are analyzed.  2h 3h 4h 5h 6h 7h 8h 9h 10h 11h 12h 13h 14h 15h 16h 17h 18h 19h 20h 21h

Analysis of the Differences among County and City
Units. Furthermore, this paper takes the administrative region (city or county) of each detection site as a unit to explore the differences in the population flow and logistics in each county or city. First, we calculate the outflow, inflow, and total flow (sum of the inflow and outflow) of the population flow and material flow in each city or county, respectively (Tables 3 and 4). en, the natural breakpoint classification method [37] is adopted to classify the total population flow and total material flow in each city or county. e classification results are shown in Figures 15(a) and 15(b). e classification results in Figure 15 show the following.

Analysis of the Influence Factors of the Intensity of the
TToF. Here, we explore the possible factors affecting the intensity of the TToF in each county and city. Based on the availability of data, 23 indexes that may impact the intensity are preliminarily selected from "the statistical yearbook of Taiwan Province 2016" [38]. en, adopting the method of literature [39], the correlation coefficient analysis and variation coefficient analysis are carried out for each index to delete the redundant indexes. Finally, 12 indexes are obtained that are used in the impact factor analysis. Finally, the statistical product and service solutions (SPSS) software is used to analyze the correlation between the total population  flow, the total material flow, and the relevant index data of each county and city. e correlation matrices we get are shown in Tables 5 and 6. e following can be seen from the correlation matrix in Tables 5 and 6   result for the TToF. (c) e correlation between the intensity of the material flow and STC is slightly higher than that with the BTC. It confirms that the proportion of the small-scale material flow with respect to the total material flow is higher than that of the large-scale material flow. (3) From a comprehensive point of view, different industries and structures will also have great impacts on the intensity of the TToF, which is shown by the correlations between the   intensity of the TToF and the NPBPI, NPBSI and NPBTI being significantly different. e contribution of the second industry to the TToF is greater than that of the third industry, and it is further greater than that of the first industry.

Spatial and Temporal Autocorrelation Analysis.
To further explore the spatial-temporal correlation of Taiwan's expressway network, ST-ACF and CCF are adopted to analyze the global spatial-temporal autocorrelation and local spatial-temporal autocorrelation, respectively.

Global Spatial-Temporal Autocorrelation Analysis.
In this paper, the spatial first-order adjacency matrix and the second-order adjacency matrix of each road section are used as the spatial weight matrix to calculate the ST-ACF index of all road sections. Meanwhile, to better reflect the temporal and spatial differences in the spatial-temporal correlation of Taiwan's expressway network, we analyze the changes of ST-ACF index in the early peak (8:00-11:00), the middle peak (11:00-16:00), and the late peak (16:00-18:00) periods. e calculation results of the ST-ACF index are shown in Figure 16. e following can be seen from Figure 16. (1) e ST-ACF value changes periodically over the whole day ( Figure 16(d)), and the change period is consistent with the time length of the analysis (every 5 minutes, 288 intervals in the whole day). (2) e spatial-temporal correlation is negatively correlated with the spatial distance, which is reflected by the ST-ACF value in Figure 16(d) decreasing as the spatial order increasing. (a) When the spatial order is first order, the ST-ACF value is greater than 0 and the degree of confidence is greater than 95%, showing a significant positive correlation. (b) e spatial-temporal correlation is not significant when the spatial order is second order. is shows that the spatial-temporal correlation of Taiwan's expressway network is mainly reflected in the direct adjacent road sections, but not in the indirect connected road sections. (3) Compared with the early peak and the late peak, the ST-ACF value in the middle peak is significantly higher, and it changes faster. As a whole, the different traffic patterns in different periods for Taiwan's expressway network show that the spatial-temporal autocorrelation has the characteristics of time dynamics and spatial heterogeneity.

Local Spatial-Temporal Autocorrelation Analysis.
e global spatial-temporal autocorrelation analysis describes Taiwan's expressway network as a whole, but it may cover up some local heterogeneity. erefore, it is necessary to analyze the local spatial-temporal correlation. Considering the fact that the core areas of the TToF are located in the northern part of Taiwan Province and considering the ranking of the intensity of the TToF in each detection section, five consecutive and intersecting road sections (as shown in Figure 17(a)) are selected from the top 50 sections from the northern part of Taiwan Province for the local spatial-temporal autocorrelation analysis. en, using the method in literature [36], the CCF values of different periods (early peak, middle peak, and late peak) in the spatial firstorder adjacency matrix are calculated. e results are shown in Figure 17.
(1) e CCF values in different road sections show obvious heterogeneity. (a) In the middle peak, with the spatial first-order adjacency matrix and time delay s being zero, the CCF values of the five road sections vary from −0.2 to 0.35. is means that there are significant differences in the correlation between some road sections and their spatial first-order adjacent road sections. (b) As shown in Figure 17(b), road sections 1, 2, and 3 show a strong positive correlation (CCF values are greater than zero). Road section 4 has a CCF close to zero, which means that it has a weak correlation with its directly adjacent road   Table 5 represent the total flow of population, subtotal of buses (vehicles), and subtotal of buses (vehicles). TMFC, BTC, and STC in Table 6 represent the total flow of material, subtotal of trucks (vehicles), and subtotal of trucks (vehicles), respectively. * * When the confidence (double test) is 0.01, the correlation is significant. * When the confidence (double test) is 0.05, the correlation is significant. e first row of data of each index is Pearson correlation, and the second row is significance (double tail). Sample capacity n is 16. sections. Meanwhile, road section 5 shows a negative correlation. (c) e correlation of each road section varies with time. For example, in the late peak period, the correlation of road sections 3 and 5 gradually changes from negative to positive, the correlation of road sections 1 and 2 changes from positive to negative, and section 4 does not change much from beginning to end. (2) e correlation of the same road section also changes in different periods. For example, the correlation of section 3 shows a positive correlation in the early peak period and a negative correlation in the middle and late peak periods. It shows that the autocorrelation structure also changes with time.

Conclusion
e population flow and material flow are fundamental elements of space of flows. e main shortcoming of most of the existing works related to TToF is that they ignore the differences of spatiotemporal characteristics in different scales and the correlation between them. To grasp the spatial-temporal differences of TToF in specific research area deeply and accurately, more detailed work is needed. In this paper, we study the spatial-temporal characteristics of the population flow and material flow in Taiwan Province and analyze the spatial-temporal autocorrelation of Taiwan's expressway network based on the ETC data of Taiwan Province in 2017. e population flow and material flow are distinguished by vehicle types in the experimental data. To fully explore the potential value of the data, we adopt a multimembership and layer-by-layer refinement method to carry out this research to discover Taiwanese travel habits and the material flow patterns of Taiwan Province. e thoughts of multimembership and layer-by-layer refinement are mainly embodied in four aspects. (1) e population flow is divided into two subtypes named the "minority population flow" and the "mass population flow," while the material flow is divided into the "small-scale material flow" and the "large-scale material flow" correspondingly. (2) e differences among the time characteristics of these kinds of flows are studied in monthly, daily, and hourly units, respectively. (3) e spatial differences of the different types of flows are studied under different scales (road sections, cities, and counties). (4) e ST-ACF and CCF are adopted to study the global and local spatial-temporal autocorrelation of Taiwan's expressway network, respectively. e results show the following. (1) Both the different temporal characteristics and the spatial characteristics of the TToF in Taiwan Province are obvious. (2) e flow intensity of the TToF is largely influenced by some objective factors, such as the population size (especially the number of employed  0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23   population), industrial structure, and so on. (3) ere is a significant positive spatial-temporal correlation in Taiwan's expressway network, but this correlation is mainly reflected in the spatial first-order adjacent road sections. e value of our work is mainly embodied in four aspects: (1) constructing a research framework for studying the TToF by adopting the multimembership and layer-by-layer refinement analysis idea, (2) helping decision makers and some other researchers to further understand Taiwanese's travel habits and the material flow patterns of Taiwan Province, (3) providing some basis for decision makers to conduct traffic planning and adjust industrial layout, and (4) showing another novel use of the ETC data. However, there are also some limitations in this paper. (1) e annual gross domestic product (GDP) of tourism is relatively considerable in Taiwan Province, and so the minority flow generated by the tourism industry will inevitably affect the total population flow of Taiwan Province. However, limited by the availability of data, factors such as the number of tourists in each county and city and their tourism incomes are not taken into account in the process of the factor analysis of the intensity of the TToF.
(2) e current work of this paper mainly focuses on the analysis of the flow "intensity" of each type of flow and fails to properly excavate the value of the flow "direction." In the follow-up work, the starting sites and ending sites and their administrative regions could be used to construct the population flow and material flow correlated network of site-to-site and county (city)-tocounty (city) networks. en, the social network analysis method could be used to evaluate the whole network and each node in the network.

Data Availability
Our data have been uploaded to Dryad and assigned with a DOI. e data are approved by Dryad curators, and the following temporary link can be used for obtaining the data: https://datadryad.org/stash/share/OhiOnewzuI59cxMXgR WEIma4TqS_05hUulUQGRlOCkQ. e data can also be downloaded from the following website: http://tisvcloud. freeway.gov.tw/history/TDCS/M05A/.

Conflicts of Interest
e authors declare that they have no conflicts of interest.