People flow estimation system, estimation parameter calculation device, people flow estimation device, estimation parameter calculation method, people flow estimation method, estimation parameter calculation program, and people flow estimation program.

Abstract

To provide a people flow estimation system capable of capturing the flow of people regardless of the presence or absence of transactions, as well as, capable of achieving high accuracy in the estimated flow of people.SOLUTION: A people flow estimation system 1 is constituted of: a stay point extraction unit 31, a statistical correction unit 34, and an estimation parameter calculation unit 12 that, using spatial base station data 11 acquired by a base station and relating to a communication terminal device located in a predetermined space and spatial trajectory spatial data 51 relating to a trajectory of the communication terminal device moving in the space, calculate an estimation parameter relating to a ratio of the number of communication terminal devices indicated by the spatial base station data to the number of communication terminal devices staying in the space indicated by the spatial trajectory data; a visited location estimation unit 32 that calculates the number of the communication terminal devices staying in a feature existing in the space, using the spatial trajectory data 51; and an estimation unit 33 that estimates the number of the communication terminal devices staying in the feature, using the calculated number of the communication terminal devices and the estimation parameter.SELECTED DRAWING: Figure 2

Description

[Technical Field]

[0001] The present invention relates to a people flow estimation system, an estimation parameter calculation device, a people flow estimation device, an estimation parameter calculation method, a people flow estimation method, an estimation parameter calculation program, and a people flow estimation program. [Background technology]

[0002] Services are being offered that use portable communication terminal devices such as smartphones to grasp the state of fluid human movement (hereinafter also referred to as "people flow"). For example, such services provide information on the state of congestion in stores and the population of urban areas by time of day, and can be used for individual behavior and urban planning. Such a service for grasping people flow is described, for example, in Patent Document 1.

[0003] Patent Document 1 describes an estimation system that acquires multiple pieces of terminal device location information and transaction information indicating the details of transactions conducted by the terminal device owners, calculates the number of terminal device owners who visited a facility where a transaction indicated in the transaction information was conducted, and generates statistical information regarding the number of people who visited a specified facility based on the number of terminal devices and the number of people who are determined to have visited the specified facility. Because the estimation system described in Patent Document 1 generates statistical information based on the number of people who conducted some kind of transaction at the facility, there is a concern that the accuracy of the estimation may be low if there is no consistent relationship between the number of people who conducted a transaction and the number of people who visited the facility but did not conduct a transaction. Furthermore, an estimation system that estimates statistical data based on transaction information can only be applied to facilities where the terminal device owners conduct transactions, limiting its scope of application. [Prior art documents] [Patent documents]

[0004] [Patent Document 1] Patent Publication No. 2021-43646 Summary of the Invention [Problem to be solved by the invention]

[0005] The present invention has been made in consideration of the above points, and aims to provide a people flow estimation system, an estimation parameter calculation device, a people flow estimation device, an estimation parameter calculation method, a people flow estimation method, an estimation parameter calculation program, and a people flow estimation program that can grasp people flow regardless of whether a transaction is taking place or not, and that estimates people flow with high accuracy. [Means for solving the problem]

[0006] In order to achieve the above object, one aspect of the present invention is a people flow estimation system comprising: an estimation parameter calculation device that uses spatial base station data acquired by a base station relating to communication terminal devices within a predetermined space and spatial trajectory data relating to trajectories of communication terminal devices moving within the space to calculate an estimation parameter relating to a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices remaining in the space indicated by the spatial trajectory data; a feature terminal number calculation unit that uses the spatial trajectory data to calculate the number of communication terminal devices remaining in a feature present in the space; and an estimation unit that uses the number of communication terminal devices calculated by the feature terminal number calculation unit and the estimation parameter to estimate the number of communication terminal devices present in the feature.

[0007] One aspect of the present invention is an estimation unit that estimates the number of communication terminal devices present in a feature using a calculated number of communication terminal devices remaining in a feature present in a space and an estimation parameter related to a ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining in the space, Using spatial base station data relating to communication terminal devices within a predetermined space, which are acquired by a base station, and spatial trajectory data relating to the trajectories of communication terminal devices moving within the space, a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices staying within the space indicated by the spatial trajectory data is used. Therefore, the above The estimation parameter calculation device includes an estimation parameter calculation unit that calculates estimation parameters.

[0008] One aspect of the present invention is a people flow estimation device that includes an estimation parameter acquisition unit that acquires estimation parameters related to the ratio between the number of communication terminal devices indicated by spatial base station data related to communication terminal devices within a specified space acquired by a base station and the number of communication terminal devices remaining in the space indicated by spatial trajectory data related to trajectories of communication terminal devices moving within the specified space; a feature terminal number calculation unit that uses the spatial trajectory data to calculate the number of communication terminal devices remaining within a feature present in the space; and an estimation unit that uses the number of communication terminal devices calculated by the feature terminal number calculation unit and the estimation parameters to estimate the number of communication terminal devices present within the feature.

[0009] One aspect of the present invention is an estimation unit that estimates the number of communication terminal devices present in a feature using a calculated number of communication terminal devices remaining in a feature present in a space and an estimation parameter related to a ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining in the space, Using spatial base station data relating to communication terminal devices within a predetermined space, which are acquired by a base station, and spatial trajectory data relating to the trajectories of communication terminal devices moving within the space, a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices staying within the space indicated by the spatial trajectory data is used. Therefore, the above The method includes calculating the estimation parameters.

[0010] One aspect of the present invention is By computer, The people flow estimation method includes: acquiring an estimation parameter related to the ratio between the number of communication terminal devices indicated by spatial base station data related to communication terminal devices within a specified space acquired by a base station and the number of communication terminal devices remaining within the specified space indicated by spatial trajectory data related to trajectories of communication terminal devices moving within the specified space; calculating the number of communication terminal devices remaining within a feature present in the space using the spatial trajectory data; and estimating the number of communication terminal devices present within the feature using the calculated number of communication terminal devices remaining within the feature and the estimation parameter.

[0011] One aspect of the present invention is an estimation unit that estimates the number of communication terminal devices present in a feature using a calculated number of communication terminal devices remaining in a feature present in a space and estimation parameters related to a ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining in the space,Using spatial base station data relating to communication terminal devices within a predetermined space, which are acquired by a base station, and spatial trajectory data relating to the trajectories of communication terminal devices moving within the space, a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices staying within the space indicated by the spatial trajectory data is used. Therefore, the above This is an estimation parameter calculation program that enables a computer to realize the function of calculating estimation parameters.

[0012] A people flow estimation program according to one embodiment of the present invention causes a computer to implement the following functions: a function to acquire estimation parameters relating to the ratio between the number of communication terminal devices indicated by spatial base station data relating to communication terminal devices within a specified space acquired by a base station and the number of communication terminal devices remaining within the specified space indicated by spatial trajectory data relating to the trajectories of communication terminal devices moving within the specified space; a function to calculate the number of communication terminal devices remaining within a feature present in the space using the spatial trajectory data; and a function to estimate the number of communication terminal devices present within the feature using the calculated number of communication terminal devices remaining within the feature and the estimation parameters. [Effects of the Invention]

[0013] According to the above aspects, it is possible to provide a people flow estimation system, an estimation parameter calculation device, a people flow estimation device, an estimation parameter calculation method, a people flow estimation method, an estimation parameter calculation program, and a people flow estimation program that can grasp people flow regardless of whether or not a transaction is taking place, and that estimates people flow with high accuracy. [Brief explanation of the drawings]

[0014] [Figure 1] 1A, 1B, 1C, and 1D are diagrams each illustrating an outline of one embodiment of the present invention. [Figure 2] FIG. 1 is a functional block diagram for explaining a people flow estimation system according to an embodiment of the present invention. [Figure 3] 3 is a diagram for explaining a process performed in a stay point extraction unit shown in FIG. 2. FIG. [Figure 4]10A, 10B, 10C, and 10D are diagrams for explaining the process of extracting a stagnation point. [Figure 5] 3A, 3B, and 3C are diagrams illustrating the processing performed in the statistical correction unit shown in FIG. 2. [Figure 6] FIG. 10 is a diagram for explaining a process for determining at which feature a user of a communication terminal device is staying. [Figure 7] 1 is a flowchart illustrating a people flow estimation method according to an embodiment of the present invention. [Figure 8] FIG. 8 is a schematic diagram for explaining each process of the flowchart shown in FIG. 7. DETAILED DESCRIPTION OF THE INVENTION

[0015] An embodiment of the present invention will be described below. Note that the people flow estimation system, estimation parameter calculation device, people flow estimation device, estimation parameter calculation method, people flow estimation method, estimation parameter calculation program, and people flow estimation program described in this embodiment are merely examples, and aspects of the present invention are not limited to the configurations, processing procedures, etc. described herein.

[0016] [overview] 1(a), 1(b), 1(c), and 1(d) are diagrams for explaining an overview of this embodiment. FIG. 1(a) shows a mesh space M, and FIG. 1(b) shows a building B within the mesh space M. FIG. 1(c) shows a movement trajectory t of a communication terminal device moving within the mesh space M. FIG. 1(d) shows a communication terminal device staying within building B. In this specification, "staying" refers to the movement trajectory of a communication terminal staying within a predetermined range (range threshold) for a predetermined time (time threshold) or more. In this embodiment, the time threshold may be set to, for example, 10 minutes or 15 minutes. Furthermore, in this embodiment, all of the movement trajectories of a communication terminal device, excluding partial trajectories determined as "staying," are considered to be trajectories on which the communication terminal device has "moved." Furthermore, all communication terminal devices present within the mesh space, including both "staying" and "moving" terminals, are considered to be "staying" communication terminal devices. The range of the mesh space M may be, for example, about 500 m×500 m, and may be a rectangular section when viewed from above, like the mesh space M shown in FIG. 1(a) and the like.

[0017] The number of communication terminal devices residing within the mesh space M shown in FIG. 1(a) can be obtained, for example, by a base station communicating with communication terminal devices within the mesh space. The base station generates communication information for each corresponding area, including communication (call and data transmission / reception) of the communication terminal devices and communication information for periodically tracking the locations of the communication terminal devices. The history of the generated information (base station communication history (CDR)) includes location information of the communication terminal devices within the area. The location information of the communication terminal devices included in the base station communication history will hereinafter be referred to as "base station data." It is believed that the number of users (people) residing in the mesh space M approximately matches the number of communication terminal devices. Therefore, the number of users of communication terminal devices residing within the mesh space M can be accurately obtained from the base station data.

[0018] Furthermore, as shown in FIG. 1(c), the movement trajectory of the communication terminal is represented by location information and time information associated with the location information. In this specification, a pair of location information and time information will be referred to as "trajectory data." The trajectory data is acquired, for example, by a Global Navigation Satellite System (GNSS). Positioning of the communication terminal device by the GNSS may be performed, for example, when an application that performs positioning of the communication terminal device is used in the communication terminal device.

[0019] Let's consider estimating the number of communication terminal owners (hereinafter simply referred to as "owners") residing within Building B in a mesh space. The location information of communication terminals obtained from a base station relates to communication terminals within a range of approximately 40 to 70 meters from the base station, making it difficult to identify the number of communication terminals residing within a specific feature. Note that a feature here refers to any object on the ground, regardless of whether it is man-made or natural. On the other hand, while GNSS positioning can determine the number of communication terminals moving within Building B and their movement trajectories, this requires a sufficient amount of positioning data for each communication terminal. Therefore, measuring communication terminals within Building B using trajectory data is difficult to apply to any communication terminal other than those that acquire positioning data relatively frequently, and it is difficult to collect trajectory data for a sufficient number of communication terminals to obtain highly accurate statistical results.

[0020] From the above, in this embodiment, base station data (spatial base station data) for the entire mesh space M, which can accurately observe communication terminal devices within the mesh space M, and trajectory data (spatial trajectory data) for the entire mesh space M, which can acquire more trajectory data than within building B, are acquired. The communication terminal devices whose number is determined by the spatial trajectory data are a portion of the communication terminal devices whose number is determined by the spatial base station data. In this embodiment, the number obtained by the spatial base station data is calculated relative to the number of staying owners obtained by the spatial trajectory data, and estimation parameters are created from this ratio. Then, assuming that the created estimation parameters are also applicable within building B, the number of people staying in building B is estimated by applying the estimation parameters to the number of staying people based on the feature trajectory data acquired within building B.

[0021] In the above process, it is possible to associate the spatial base station data with the planning information of the communication terminal device. Therefore, according to this embodiment, it is possible to identify attributes of the user, such as gender, age, occupation, and residential area, estimated based on the spatial base station data, and acquire estimation parameters (attribute estimation parameters) for each user attribute. Furthermore, according to this embodiment, it is possible to acquire map data corresponding to the mesh space M and acquire information related to the attributes of building B, such as commercial facilities, office buildings, and public facilities. Then, it is possible to correct the attribute estimation parameters according to the attributes of building B.

[0022] Furthermore, in this embodiment, in order to increase the number of feature trajectory data acquired and improve the accuracy of the estimation results, trajectory data from a communication terminal device that has acquired a large number of data is used to correct trajectory data from a communication terminal device that has acquired a small number of data. This type of processing increases the number of feature trajectory data that can be used to calculate estimation parameters, making it possible to improve the accuracy of the estimated values. This type of aspect will be described below.

[0023] [People Flow Estimation System] FIG. 2 is a functional block diagram illustrating a people flow estimation system according to this embodiment. The people flow estimation system 1 shown in FIG. 2 includes an estimation parameter calculation device and a people flow estimation device. The stagnation point extraction unit 31, the visited destination estimation unit 32, the estimation parameter calculation unit 12, and the statistical correction unit 34 of the people flow estimation system 1 constitute the estimation parameter calculation device. The estimation unit 33 corresponds to the people flow estimation device. Each function may be a program that works in conjunction with computer hardware and other software to calculate estimation parameters and / or estimate people flow based on the estimation parameters. In other words, the program that constitutes the people flow estimation system may function to both calculate people flow estimation parameters and estimate people flow.

[0024] The spatial base station data 11, spatial trajectory data 51, high-density trajectory data 52, staying point data 53, visited place estimation result data 54, in-building occupancy data 55, low-density trajectory data 56, and corrected staying point data 57 are all data acquired or generated by the estimation parameter calculation device or people flow estimation device, or a memory unit that stores data temporarily or for a long period of time. Such data or memory unit may be stored inside or externally to the people flow estimation system, or may even be located on a network.

[0025] Furthermore, the people flow estimation system 1 shown in FIG. 2 is capable of calculating estimation parameters and estimating people flow in parallel. However, this embodiment is not limited to this configuration, and estimation parameters may be calculated in advance and used to estimate people flow in other situations. In such cases, for example, by matching the conditions used when calculating the estimation parameters with the conditions used when estimating people flow, it is possible to increase the accuracy of the estimated number of users. Possible conditions include, for example, conditions related to the area or date and time corresponding to the mesh space, and conditions related to building attributes.

[0026] [Estimation parameter calculation device] As described above, the estimation parameter calculation device includes a stay point extraction unit 31, a visited place estimation unit 32, an estimation parameter calculation unit 12, and a statistical correction unit 34. As shown in FIG. 2 , the estimation parameter calculation unit 12 acquires spatial base station data of a communication terminal device located within a mesh space, which is a predetermined space, acquired by a base station, and spatial trajectory data relating to the trajectory of the communication terminal device moving within the mesh space. The spatial trajectory data is data in which positioning data obtained by GNSS is arranged in chronological order, pairing the position of the communication terminal device with the time at which the position was acquired. In this embodiment, the spatial trajectory data is created using positioning data obtained by a GPS (Global Positioning System). The spatial base station data and the spatial trajectory data are temporarily saved as spatial base station data 11 and spatial trajectory data 51, respectively. In this embodiment, the spatial trajectory data 51 includes high-density trajectory data, which has a relatively large number of positioning data points within a predetermined time range, and low-density trajectory data, which has a relatively small number of positioning data points. The spatial trajectory data is divided into high-density trajectory data 52 and low-density trajectory data 56 and stored.

[0027] High-density trajectory data 52 refers to data containing a predetermined number of positioning data points within a predetermined time range, while low-density trajectory data 56 refers to data containing fewer than a predetermined number of positioning data points within the same predetermined time range. The high-density trajectory data 52 allows for accurate estimation of the location where a communication terminal device stayed in a mesh space. On the other hand, the low-density trajectory data 56 contains a small number of positioning data points, making it difficult to estimate the location where a communication terminal device stayed in a mesh space. The people flow estimation system 1 includes a stay point extraction unit 31 that calculates the proportion of stay time of a communication terminal device in a mesh space based on the high-density trajectory data, and a statistical correction unit 34 that corrects the stay time of a communication terminal device in a mesh space that has obtained low-density trajectory data, using the proportion of stay time calculated by the stay point extraction unit 31. In other words, the stay point extraction unit 31 and the statistical correction unit 34 constitute a stay correction unit in the people flow estimation system 1.

[0028] (Retention point extraction part) Next, the stay point extraction unit will be described. FIG. 3 is a diagram illustrating the processing performed by the stay point extraction unit 31. The stay point extraction unit 31 extracts stay points based on high-density trajectory data and creates or updates spatial trajectory data each time a communication terminal device acquires positioning data. The spatial trajectory data may include, for example, positioning data acquired by a single communication terminal device in a mesh space over the course of a day. The stay point extraction unit 31 divides such spatial trajectory data into a stay point portion and a moving portion. This division is performed, for example, for a single communication terminal device, and positioning data within a predetermined window width from the core point is sequentially acquired. Moving portions t1 and t2 are formed by consecutive positioning data. Furthermore, a set of positioning data within a window set around a core point, which is a point midway between the moving portions (movement trajectories) t1 and t2, is referred to as a "neighborhood point." If the elapsed time of a set of neighboring points whose distance range is within a threshold is equal to or greater than a stop time threshold, the set of positioning data is designated as a "stop cluster."

[0029] By repeating the above process, multiple stop clusters are generated in the moving portions t1 and t2. The stop clusters connect overlapping positioning data, and if the communication terminal device stays in the connected stop cluster for a stay time threshold or longer, the range of the connected clusters is designated as a "stay point" P. In addition, positioning data that is a distance greater than a noise threshold, which is a preset distance from the moving portions t1 and t2 or the stay point P, is removed as a noise point Np. The noise threshold may be, for example, approximately 100 m to 300 m. Note that this process can be realized, for example, by a stay point extraction (D-star) technique.

[0030] 4(a), 4(b), 4(c), and 4(d) are diagrams for explaining the above processing. FIG. 4(a) is a diagram showing positioning data p including core point pc acquired first within the range of window width 41. The positioning data shown in FIG. 4(a) is acquired in order from left to right in the figure, as shown on the time axis. Note that the window width is a concept set in relation to the positioning order, and does not take time or distance into consideration.

[0031] FIG. 4(b) is a diagram illustrating the process of removing noise points Np. Note that FIG. 4(b) is a schematic diagram illustrating the process, showing the order in which noise points Np and positioning data p are acquired. The noise points Np exist within the range of the window width. As shown in FIG. 4(b), during the acquisition of consecutive positioning data p or core points pc, points whose positioning order is close to the core point pc are extracted as "neighboring points." In this case, if positioning data is acquired that is farther away from the core point pc than the noise threshold, the acquired positioning data is determined to be a noise point Np. The noise points Np are deleted when stop clusters are combined, as described below.

[0032] 4(c) is a diagram illustrating a process of creating a stop cluster by applying a time threshold and a range threshold. As shown in FIG. 4(c), when the elapsed time from the acquisition time of the first acquired positioning data pf to the acquisition time of the last acquired positioning data pl (including noise point Np) within the range threshold reaches the time threshold, the positioning data pf, pl, and the positioning data acquired between the acquisition of the positioning data pf and the acquisition of the positioning data pl are treated as a stop cluster 45 (clustered).

[0033] 4(d) is a diagram for explaining the joining of multiple stop clusters. As shown in FIG. 4(d), multiple stop clusters 45 and 46 are created by the process shown in FIG. 4(c), and if the positioning data pd included in stop cluster 45 is also included in stop cluster 46, stop clusters 45 and 46 are joined. If the time interval between the acquisition time of the first acquired positioning data and the acquisition time of the last acquired positioning data among the positioning data included in the joined stop clusters exceeds the stay threshold, the range corresponding to the stop cluster becomes stay point P. Stay point extraction unit 31 repeats the above process to extract stay points of users of communication terminal devices moving within the mesh space.

[0034] (Statistical correction section) The statistical correction unit 34 corrects the low-density trajectory data to increase the number of extracted stagnation points and perform processing to improve the accuracy of the estimated values. Figures 5(a), 5(b), and 5(c) are diagrams for explaining the processing performed by the statistical correction unit 34. Figure 5(a) is a diagram for explaining the correction processing of low-density trajectory data. Figure 5(b) is a diagram for explaining the processing of calculating the proportion of stagnation time in the observable time, which is used in the processing shown in Figure 5(a). Figure 5(c) is a diagram showing the stagnation time at each feature obtained as a result of the processing shown in Figure 5(a). As shown in Figure 5(b), the statistical correction unit 34 corrects the high-density trajectory data d of multiple communication terminal devices. gb1 From d gbn The statistical correction unit 34 then calculates the relationship between the observable time and the proportion of the staying time. Here, the observable time refers to the time interval from the time when a specific communication terminal device first acquires positioning data in the mesh space to the time when it last acquires positioning data. The proportion of the staying time refers to the ratio of the time that the communication terminal device was at the staying point to the observable time. In the example shown in FIG. 5(b), when the observable time is 1 hour, the proportion of the staying time is 0.9, and when the observable time is 5 hours, the proportion of the staying time is 0.6.

[0035] The statistical correction unit 34 calculates the density of the plurality of low-density spatial trajectory data d shown in FIG. ga1 From d gan In the example shown in FIG. 5(a), for example, if the time interval between successively acquired positioning data p1, p2, and p3 exceeds a stay threshold, the positions indicated by the positioning data p1, p2, and p3 may be regarded as stay points. Alternatively, positioning data may be interpolated between the positioning data p1, p2, and p3, and all of them may be regarded as stay points. The statistical correction unit 34 determines the proportion of the stay time according to the observation time from when the positioning data p1 to p3 is acquired. Then, the stay time is evenly allocated to the positions indicated by the positioning data p1, p2, and p3 that are regarded as stay points. The corrected positioning data d shown in FIG. 5(a) ga1 ´ to d gan5(c), the positions indicated by the positioning data p1, p2, and p3 are associated with the features. As a result, the statistical correction unit 34 can calculate the stay time at each of the features, such as the A Building, the TA Building, and the home improvement store, as shown in FIG.

[0036] In the example shown in FIG. 5(c), the observable time is 5 hours, of which 0.6 is the dwell time. Therefore, the total dwell time is 3 hours, and the dwell time is allocated equally to each feature, with 1 hour allocated to each feature. The process of associating dwell points with feature data will be described later. This process of correcting low-density trajectory data based on the proportion of dwell time in high-density trajectory data is based on the idea that low-density trajectory data is obtained by thinning out the positioning data from high-density trajectory data. The statistical correction unit 34 considers positioning data located along the movement portion indicated by the low-density trajectory data as dwell points and allocates the dwell time of the communication terminal device to each dwell point equally. The corrected dwell point data 57 is accumulated together with the dwell point data 53 and input to the visited place estimation unit 32 and the estimation parameter calculation unit 12.

[0037] (Visited Place Estimation Department) The visited place estimation unit 32 has a function of associating the stay points indicated by the stay point data 53 and the corrected stay point data 57 with features. FIG. 6 is a diagram for explaining the process of associating stay point P with features, i.e., determining at which feature a user is staying. FIG. 6 shows stay point P and buildings Ba, Bb, and Bc as features. The visited place estimation unit 32 associates stay point P with a building that satisfies one of the following three conditions: (1) a building located in an area that includes stay point P; (2) a building whose distance d from stay point P is equal to or less than a predetermined distance threshold; or (3) a building that is at least partially included in a range A of a circle with a radius R centered on stay point P. If there are features that meet multiple conditions, the association is prioritized in the order of (1), (2), and (3). The user of the communication terminal device located at stay point P is determined to be staying in the associated building. The visited place estimation unit 32 generates visited place estimation result data 54 that identifies the building in which the user of the communication terminal device stayed, using data indicating the stay point P in the mesh space and one of buildings Ba, Bb, or Bc associated with this stay point P. In addition, in this embodiment, the user staying in a building will hereinafter also be referred to as "visiting" the building.

[0038] In the above condition (1), the relationship between the position of the stay point P and the area occupied by the building may be determined from map information acquired by the visited place estimation unit 32. When the stay point P is expressed by latitude and longitude, the map information preferably includes latitude and longitude information. The distance d and radius R in conditions (2) and (3) may be, for example, the distance between the center of gravity of the stay point P and the center of gravity of the building. The visited place estimation unit 32 sends visited place estimation result data 54 to the estimation unit 33.

[0039] (Estimation parameter calculation section) The estimation parameter calculation unit 12 receives input of the spatial base station data 11 and visited place estimation data 54 (i.e., spatial trajectory data 51 associated with features). The estimation parameter calculation unit 12 then calculates the number of users present in the mesh space using the spatial base station data 11, and calculates the number of users staying in the mesh space using the visited place estimation data 54. Hereinafter, in this specification, the number of users calculated using the spatial base station data 11 will be referred to as the mobile population, and the number of users staying in the mesh space calculated using the spatial trajectory data 51 will be referred to as the GPS population. The calculations of the mobile population and the GPS population are performed within the same time range, which may be, for example, a day or every hour.

[0040] Here, we will explain an example of how to count the mobile population in a mesh space, the GPS population, and the GPS population within a feature. The mobile population in a mesh space is calculated, for example, from the base station communication history of base stations whose communication management ranges overlap with the mesh space, and corresponds to the number of communication terminal devices that communicated. The GPS population in a communication terminal mesh space can be obtained from the identification information and stay points of the communication terminal devices included in the positioning data. As mentioned above, the GPS population within a feature (building) is calculated by comparing the stay points acquired based on the positioning data with map information and associating the stay points with buildings that satisfy any of the above-mentioned (1) to (3). The number of users who stayed in a building, i.e., visited a building, within a specified time range is calculated as follows.

[0041] [Table 1] Table 1 above shows an example of the number of users staying in building Ba in Figure 6 between 9:30 and 11:45. Table 1 shows an example in which one user stayed between 9:30 and 10:00, one user stayed between 10:00 and 11:00, and one user stayed between 11:00 and 11:45. In this embodiment, users are counted every hour, and times less than one hour are expressed as a decimal and multiplied by the number of users staying.

[0042] Next, the estimation parameter calculation unit 12 calculates estimation parameters by dividing the mobile population in the mesh space by the GPS population associated with buildings, etc. The estimation parameters indicate the proportion of users whose positions have been measured by GPS, that is, users of communication terminal devices whose positioning information has been acquired multiple times within a short distance, out of all users extracted from the base station communication history. This embodiment is based on the premise that the estimation parameters generated in the mesh space can also be applied within features.

[0043] In addition, the base station communication history used to calculate the mobile population is associated with the user's contract information. Furthermore, the application requesting the acquisition of positioning data stores information including the user's registration information. In this embodiment, estimation parameters may be calculated for each user's attributes based on the user's contract information and registration information. Table 2 shows an example in which the mobile population and GPS population of a mesh space are acquired by gender and age group, and estimation parameters for users in their 20s are calculated by gender. According to Table 2, the mobile population of men in their 20s in the mesh space is 1,200, the GPS population is 100, and the mobile population of women in their 20s is 100 and the GPS population is 100. The estimation parameter for men in their 20s is 12, and the estimation parameter for women in their 20s is 6.

[0044] [Table 2]

[0045] Furthermore, the estimation unit 33 may correct the estimation parameters based on the attributes of the building, which is a feature. Tables 3 and 4 show examples of correcting the estimation parameters for a male in his twenties in Table 2 based on the use as an attribute of building Ba shown in FIG. 6. The example shown in Table 3 is an estimation parameter obtained by correcting the estimation parameter 12 shown in Table 2. In the example in Table 3, the estimation parameters are corrected based on the use of the building, whether it is a weekday or a holiday, and the time of day. For example, if a company that uses building Ba as an office uses a flextime system and many employees start work after 11:00, the estimation parameters for the morning are corrected to be smaller (from 12 to 10). Also, on holidays, fewer users visit building Ba, which is an office, than on weekdays. By taking these conditions into consideration, the estimation parameters for holidays are corrected to be smaller (from 12 to 3).

[0046] [Table 3]

[0047] [Table 4]

[0048] Table 4 shows an example of the estimated parameters for men in their 20s adjusted based on the attributes of offices, commercial facilities, and residential areas. If Building Ba is a commercial facility that houses restaurants, and the time of day is noon, it is expected that many users will visit the restaurants. For this reason, the estimated parameters are adjusted to increase the mobile population visiting Building Ba (from 12 to 20). Also, if Building Ba is a residential area such as an apartment building, and the time of day is a weekday morning, it is expected that many users living in the residential area will head to their workplaces and will not stay at the residential area. For this reason, the estimated parameters are adjusted to decrease the mobile population visiting Building Ba (from 12 to 2).

[0049] Furthermore, in this embodiment, the estimation parameters may be corrected based on the attributes of the application that requested the acquisition of positioning data. For example, if it is known that the application has a large male user base, the estimation parameters may be corrected so that the male mobile population ratio is increased. However, this embodiment is not limited to using positioning data requested by a single application, and spatial trajectory data and stay points may be extracted using all requested positioning data regardless of the application.

[0050] (Estimation Department) The estimation unit 33 acquires visited place estimation result data 54, i.e., the GPS population of people who visited the building. The estimation unit 33 also inputs estimation parameters from the estimation parameter calculation unit 12. The visited place estimation result data 54 is then multiplied by the estimation parameters to estimate the mobile population of users staying in building Ba. The estimated mobile population is accumulated as data 55 on the number of users in the building. Table 5 is a table for explaining how the GPS population in building Ba is multiplied by the estimation parameters to calculate the estimated mobile population of users visiting building Ba. According to Table 5, when the GPS population of men in their twenties who visited building Ba is 10, the mobile population of users who visited the feature is estimated to be 120.

[0051] [Table 5]

[0052] [Person flow estimation method] 7 and 8 are diagrams illustrating the people flow estimation method of this embodiment. FIG. 7 is a flowchart illustrating the people flow estimation method, and FIG. 8 is a schematic diagram illustrating each process shown in the flowchart. The people flow estimation method includes calculating estimation parameters and estimating the number of communication terminal devices present within a feature using the number of communication terminal devices within the feature and the estimation parameters. As shown in the flowchart in FIG. 7, the people flow estimation method first acquires base station data in a mesh space (step S701). In parallel with this, the people flow estimation method acquires spatial trajectory data of communication terminal devices within the same mesh space (step S702). Then, it is determined whether the spatial trajectory data is high-density trajectory data, in which the number of positioning data acquired within a predetermined time is equal to or greater than a predetermined threshold (step S703). If the positioning data is high-density trajectory data (step S703: YES), the method proceeds to calculation of estimation parameters (step S705).

[0053] On the other hand, if the positioning data is low-density trajectory data (step S703: NO), this data is corrected (step S704). The corrected low-density trajectory data, together with the high-density trajectory data, proceeds to the calculation of estimated parameters (step S705). The estimated parameters calculated using the high-density trajectory data and the corrected low-density trajectory data are used in the estimation process of the user within the feature (step S706). After the above process, it is determined whether or not the process has ended (step S707). If an instruction to end the process is received, the process ends (step S707: YES), and the acquisition of spatial base station data and spatial trajectory data continues until an instruction to end is received (step S707: NO).

[0054] The process of calculating the estimation parameters and the process of estimating within a feature will be explained using Figure 8. The process of calculating the estimation parameters is performed by the stay point extraction unit 31, statistical correction unit 57, and estimation parameter calculation unit 12 shown in Figure 2. The process of estimating users within a feature is performed by the visited place estimation unit 32 and estimation unit 33. The people flow estimation method of this embodiment calculates the number of users, i.e., the mobile population, from spatial base station data 11 within a mesh. The number of users, i.e., the GPS population, is calculated from spatial trajectory data in the mesh. The GPS population is calculated by counting the number of communication terminal devices staying within the mesh space, i.e., the number of users who visited the mesh space.

[0055] The estimation parameters are calculated by dividing the mobile population by the GPS population. The estimation parameters may also be calculated for each age group and gender of users. In the example shown in Figure 8, the estimation parameters are calculated for teenage males, teenage females, and males in their twenties. The age group and gender of users can be obtained from the contract information of the communication device and the registration information of the application that acquires the positioning data. As shown in Figure 8, some of the communication devices detected by the base station data are also detected by the positioning data, so the estimation parameters will be a value greater than or equal to 1.

[0056] In addition, the people flow estimation method uses, for example, spatial trajectory data d of multiple users. m1 ···d mn By associating each of the staying points with the features, the feature trajectory data d g1 ···d gn Generate multiple feature trajectory data d g1 ···d gn These data constitute the visited place estimation result data 54 in FIG. 2. g1 ···d gn By aggregating the data, it becomes possible to count the number of users who visited building Ba, for example. Feature trajectory data is aggregated, for example, every hour. In the example shown in FIG. 8, the number of visitors to building Ba from 2:00 PM to 10:00 PM on January 30, 2021, is aggregated every hour.

[0057] The people flow estimation method of the present embodiment can be executed by a computer, and the execution of the method by a computer is performed by a program that causes the computer to realize each of the above functions.

[0058] As described above, this embodiment estimates people flow using base station data and positioning data, making it possible to grasp people flow regardless of whether a transaction is occurring. Furthermore, because the base station data and positioning data are obtained by measurement, the accuracy of the estimated people flow can be improved. Since spatial trajectory data generated by positioning data is associated with features, it becomes possible to estimate the GPS population of users on a feature-by-feature basis. Furthermore, because the GPS population of users who visited a feature is multiplied by an estimation parameter based on the base station data and positioning data, it is possible to estimate the mobile population of users on a feature-by-feature basis and improve the accuracy of the estimate.

[0059] The information on people flow obtained by this embodiment is useful information for, for example, new store opening plans and urban transportation plans. When opening a new store, it is possible to open the new store in a location with more visitors. Furthermore, if there are many stores visited by men in their 20s but no stores visited by men in their 50s, it is possible to open a store that caters to the preferences of men in their 50s. Furthermore, if visitors gather at Station A around 11 p.m. on weekdays, the people flow data is useful for determining the number of taxis to dispatch during this time period.

[0060] It should be noted that this embodiment is not limited to the above-described embodiment. For example, in the above-described embodiment, the features are buildings, but the features may be other features such as parks. [Explanation of symbols]

[0061] 1. People flow estimation system 11 Spatial base station data 12 Estimation parameter calculation section 31 Retention point extraction part 32 Visited Place Estimation Department 33 Estimation Department 34 Statistical correction section 51 Spatial trajectory data 52 High-density GPS trajectory data 53 Retention point data 54 Visited place estimation result data 55 Building occupancy data 56 Low-density GPS track data 57 Corrected retention point data

Claims

1. an estimation parameter calculation device that uses spatial base station data relating to communication terminal devices within a predetermined space, acquired by a base station, and spatial trajectory data relating to trajectories of communication terminal devices moving within the space, to calculate an estimation parameter relating to a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices staying within the space indicated by the spatial trajectory data; a people flow estimation device including: a feature terminal number calculation unit that calculates the number of communication terminal devices staying in a feature that exists in the space using the spatial trajectory data; and an estimation unit that estimates the number of communication terminal devices that exist in the feature using the number of communication terminal devices calculated by the feature terminal number calculation unit and the estimation parameters; A people flow estimation system equipped with the following.

2. An estimation parameter calculation device used in a people flow estimation device, including an estimation unit that estimates the number of communication terminal devices present within a feature using a calculated number of communication terminal devices remaining within a feature present within a space and estimation parameters related to the ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining within the space, an estimation parameter calculation device comprising an estimation parameter calculation unit that uses spatial base station data relating to communication terminal devices within a specified space, acquired by the base station, and spatial trajectory data relating to trajectories of communication terminal devices moving within the space, to calculate the estimation parameters based on a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices remaining within the space indicated by the spatial trajectory data.

3. 3. The estimation parameter calculation device according to claim 2, further comprising a stay correction unit that calculates a proportion of a stay time of a communication terminal device in the space based on positioning data of a communication terminal device having a relatively large number of pieces of positioning data within a predetermined time range, and corrects the number of stay points, which are ranges in which a communication terminal device having a relatively small number of pieces of positioning data is stopped on the trajectory within the space, using the proportion of the stay time.

4. 4. The estimation parameter calculation device according to claim 2, wherein the estimation parameter calculation unit calculates the estimation parameters for each attribute of a subscriber of a communication terminal device indicated by the spatial base station data and a communication terminal device staying in the space.

5. an estimation parameter acquisition unit that acquires an estimation parameter related to the ratio between the number of communication terminal devices indicated by spatial base station data related to communication terminal devices within a predetermined space acquired by a base station and the number of communication terminal devices staying within the space indicated by spatial trajectory data related to trajectories of communication terminal devices moving within the predetermined space; a feature terminal number calculation unit that calculates the number of communication terminal devices staying within a feature that exists in the space using the spatial trajectory data; A people flow estimation device comprising an estimation unit that estimates the number of communication terminal devices present within the feature using the number of communication terminal devices calculated by the feature terminal number calculation unit and the estimation parameters.

6. 6. The people flow estimation device according to claim 5, wherein a time threshold, a range threshold, and a distance threshold are set in advance, and when the spatial trajectory data indicates that the communication terminal device is in a range within the range threshold for a time equal to or greater than the time threshold, the feature terminal number calculation unit creates a stagnation point for the communication terminal device within the predetermined range, and associates the stagnation point with any one of a feature in a range including the stagnation point, a feature whose distance to the stagnation point is equal to or less than the distance threshold, and a feature at least a part of which is included in a range of a circle centered on the stagnation point and whose radius is equal to or less than the diameter threshold.

7. The people flow estimation device according to claim 5 or 6, wherein the estimation unit corrects the estimation parameters based on attributes of the features.

8. An estimation parameter calculation method for use in a people flow estimation device, comprising: an estimation unit that estimates the number of communication terminal devices present within a feature using a calculated number of communication terminal devices remaining within a feature present within a space and estimation parameters related to the ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining within the space, The computer An estimation parameter calculation method, comprising: using spatial base station data relating to communication terminal devices within a specified space, acquired by the base station, and spatial trajectory data relating to trajectories of communication terminal devices moving within the space, and calculating the estimation parameter based on a ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices remaining within the space indicated by the spatial trajectory data.

9. By a computer, Acquiring an estimation parameter related to the ratio between the number of communication terminal devices indicated by spatial base station data related to communication terminal devices within a predetermined space acquired by a base station and the number of communication terminal devices staying within the space indicated by spatial trajectory data related to trajectories of communication terminal devices moving within the predetermined space; Calculating the number of communication terminal devices staying within a feature present in the space using the spatial trajectory data; and estimating the number of communication terminal devices present within the feature using the calculated number of communication terminal devices staying within the feature and the estimation parameters.

10. An estimation parameter calculation program for use in a people flow estimation device, comprising: an estimation unit that estimates the number of communication terminal devices present within a feature using a calculated number of communication terminal devices remaining within a feature present within a space and estimation parameters related to the ratio between the number of communication terminal devices acquired by a base station in the space and the number of communication terminal devices remaining within the space, An estimation parameter calculation program that causes a computer to realize a function of calculating the estimation parameter by using spatial base station data related to communication terminal devices within a specified space, acquired by the base station, and spatial trajectory data related to the trajectories of communication terminal devices moving within the space, and by using the ratio between the number of communication terminal devices indicated by the spatial base station data and the number of communication terminal devices remaining within the space indicated by the spatial trajectory data.

11. a function of acquiring an estimation parameter relating to the ratio between the number of communication terminal devices indicated by spatial base station data relating to communication terminal devices within a predetermined space acquired by a base station and the number of communication terminal devices remaining within the space indicated by spatial trajectory data relating to trajectories of communication terminal devices moving within the predetermined space; a function of calculating the number of communication terminal devices staying within a feature existing in the space using the spatial trajectory data; A people flow estimation program that causes a computer to realize a function of estimating the number of communication terminal devices present within the feature using the calculated number of communication terminal devices remaining within the feature and the estimation parameters.

Images