# Patent

## Title
Movement predicting method

## Abstract

A movement predicting method is disclosed. The method utilizes at least one phone communicating in and moving between neighboring cells of base stations to predict population movement in a prediction interval. Firstly, obtain the traffic volumes in the cells generated by the phone and the handover information generated by the phone moving between the cells in a day. Next, calculate the traffic volumes and the handover information to obtain a movement probability for the population moving between the cells and an average residence time that the population stays in the region of each cell in the different periods of the day. Finally, according to the data obtained, predict an appearance probability that the population appears in each region at the end point of the prediction interval.

## Background
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to a predicting method, particularly to a population movement predicting method.
2. Description of the Related Art
Population movement is an important social phenomenon, whereby population spreads globally, and whereby civilization propagates and exchanges, and whereby peoples assimilate and fuse. Population movement also correlates with propagation of infectious diseases. Based on the data collected in the outbreak of dengue fever in 2007-2008 in Rio de Janeiro the second largest city of Brazil, Brazilian researchers pointed out that population movement is a key factor of spreading dengue virus. The research urges a new understanding of the epidemiological dynamics of dengue fever.
A U.S. Pat. No. 5,572,221 disclosed a technology to predict the future population movements via comparing the past movement tracks of mobile terminals with the history records. However, the prior art has the following disadvantages: 1. The technology may infringe right of privacy because it needs to collect movement tracks of a great number of individuals and. 2. The technology is limited in usage because it cannot predict the movements occurring in a specified unit time, such as the morning or afternoon rush hours. 3. The technology cannot predict the position distribution after a specified time point because it cannot estimate the residence time of individuals.
Accordingly, the present invention proposes a movement predicting method to overcome the abovementioned problems.
SUMMARY OF THE INVENTION
The primary objective of the present invention is to provide a movement predicting method, which uses the traffic volumes and handover information of the cells of switching machines in a communication network to establish a population movement model for predicting propagation of infectious diseases or tendencies of traffic flows.
To achieve the abovementioned objective, the present invention proposes a movement predicting method, which predicts the movements of people over a prediction interval, using at least one phone moving between neighboring cells covered by a plurality of base stations. Firstly, obtain the traffic volumes generated by the phones in the cells and the handover information generated by the phones moving between the cells in one day. Next, calculate the traffic volumes and the handover information to obtain the movement probabilities that people move from one cell to the neighboring cells and the average residence time that people stay in one region covered by a cell, for different periods in one day. Next, predict the appearance probabilities that people appear in each region at the end point of the prediction interval according to the start point of the prediction interval, the movement probability, and the average residence time.
Below, the embodiments are described in detail in cooperation with the drawings to make easily understood the technical contents, characteristics and efficacies of the present invention.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 schematically shows a communication network according to one embodiment of the present invention;
FIG. 2 shows a flowchart of a movement predicting method according to one embodiment of the present invention;
FIG. 3 shows a flowchart of calculating traffic volumes and handover information to obtain movement probabilities according to one embodiment of the present invention;
FIG. 4 shows a flowchart of predicting appearance probabilities according to one embodiment of the present invention;
FIG. 5 schematically shows that population moves between four cells according to one embodiment of the present invention; and
FIG. 6 schematically shows freeway sections covered by several cells according to one embodiment of the present invention.



DETAILED DESCRIPTION OF THE INVENTION
Refer to FIG. 1, wherein a telecommunication network contains a plurality neighboring cells covered by a plurality of base stations 12. The phone inside a cell communicates and transceives data via the corresponding base station 12. The rear end of base stations 12 is a Mobile Switching Center (MSC) 10. The base station 12 records the traffic volume and handover information of each user into the MSC 10.
Refer to FIG. 2. The present invention predicts movement of at least one phone (or the user of the phone) communicating in moving between a plurality of neighboring cells covered by a plurality of base stations 12 over a prediction interval. Below, the present invention is exemplified with predicting the movements of phones' users. In Step S10, obtain the traffic volume of a phone in each cell during one day and the handover information of the phone moving between the cells. Next, in Step S12, calculate the traffic volume and the handover information to obtain the movement probabilities that the user moves from a cell to a neighboring cell in different communication periods during one day and the average residence time that the user stays in a cell in different communication periods during one day. Next, in Step S14, predict the appearance probability that the user appears in each region at the end point of the prediction interval according to the start point of the prediction interval, the movement probability, and the average residence time.
Refer to FIG. 3. Step S12 further comprises the following steps. In Step S16, divide the traffic volume and the handover information to obtain the traffic volume and the handover information at each unit time Δt, and perform statistics to obtain a sub-traffic volume ρi(τ) in a cell i during each communication period τ and a sub-handover information λi,j(τ) that the user moves from a cell i to a neighboring cell j during each communication period τ. Next, in Step S18, according to sub-handover information λi,j(τ), obtain a movement probability pi,j(τ) that the user moves from a cell i to a neighboring cell j during each communication period τ, and an incoming sub-handover information λi,in(τ) that the user moves into a cell i during each communication period τ. Furthermore, an outgoing sub-handover information λi,out(τ) that the user moves out of a cell i during each communication period τ is also obtained in Step S18. pi,j(τ), λi,in(τ) and λi,out(τ) are respectively expressed by Equations (1)-(3):




  
    
      
        
          
            p
            
              i
              ,
              j
            
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            
              λ
              
                i
                ,
                j
              
            
            ⁡
            
              (
              τ
              )
            
          
          
            
              ∑
              
                j
                ,
                
                  j
                  ≠
                  i
                
              
            
            ⁢
            
              
                λ
                
                  i
                  ,
                  j
                
              
              ⁡
              
                (
                τ
                )
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          
            λ
            
              i
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            ∑
            
              j
              ,
              
                j
                ≠
                i
              
            
          
          ⁢
          
            
              λ
              
                j
                ,
                i
              
            
            ⁡
            
              (
              τ
              )
            
          
        
      
    
    
      
        (
        2
        )
      
    
  
  
    
      
        
          
            λ
            
              i
              ,
              out
            
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            ∑
            
              j
              ,
              
                j
                ≠
                i
              
            
          
          ⁢
          
            
              λ
              
                i
                ,
                j
              
            
            ⁡
            
              (
              τ
              )
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




According to the Little's law, N=λR, wherein N is the average number of persons in a cell during a unit time, and wherein λ is the number of the pieces of the handover information that the users comes into a cell during a unit time, and wherein R is the average residence time that the users stay in a cell. Suppose that the average conversation time is tc and that the average interval between two conversations is ta. Thus, the conversation probability of each user is tc/ta. Therefore, Equation (4) and (5) are obtained:




  
    
      
        
          
            ρ
            i
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            
              
                N
                ⁡
                
                  (
                  
                    1
                    
                      t
                      a
                    
                  
                  )
                
              
              ⁢
              
                t
                c
              
            
            ⇒
            N
          
          =
          
            
              
                ρ
                i
              
              ⁡
              
                (
                τ
                )
              
            
            ⁢
            
              (
              
                
                  t
                  a
                
                
                  t
                  c
                
              
              )
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            λ
            
              i
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            
              
                
                  t
                  c
                
                
                  t
                  a
                
              
              ⁢
              λ
            
            ⇒
            λ
          
          =
          
            
              
                λ
                
                  i
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
              ⁡
              
                (
                τ
                )
              
            
            ⁢
            
              
                t
                a
              
              
                t
                c
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




Combine the Little's law with Equations (4) and (5) to obtain the average residence time Ri(τ) expressed by Equation (6):




  
    
      
        
          
            R
            i
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          R
          =
          
            
              
                
                  (
                  
                    
                      t
                      a
                    
                    
                      t
                      c
                    
                  
                  )
                
                ⁢
                
                  
                    ρ
                    i
                  
                  ⁡
                  
                    (
                    τ
                    )
                  
                
              
              
                
                  (
                  
                    
                      t
                      a
                    
                    
                      t
                      c
                    
                  
                  )
                
                ⁢
                
                  
                    λ
                    
                      i
                      ,
                      
                        i
                        ⁢
                        
                          
                        
                        ⁢
                        n
                      
                    
                  
                  ⁡
                  
                    (
                    τ
                    )
                  
                
              
            
            =
            
              
                
                  ρ
                  i
                
                ⁡
                
                  (
                  τ
                  )
                
              
              
                
                  λ
                  
                    i
                    ,
                    
                      i
                      ⁢
                      
                        
                      
                      ⁢
                      n
                    
                  
                
                ⁡
                
                  (
                  τ
                  )
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




Next, in Step S20, substitute the sub-traffic volume ρi(τ) and the incoming sub-handover information λi,in(τ) into Equation (6) to obtain the average residence time Ri(τ). If the unit time or the average conversation time is very short, or if a more precise average residence time Ri(τ) is desired, take the outgoing sub-handover information λi,out(τ) into consideration and express Ri(τ) with Equation (7):




  
    
      
        
          
            R
            i
          
          ⁡
          
            (
            τ
            )
          
        
        =
        
          
            2
            ⁢
            
              
            
            ⁢
            
              
                ρ
                i
              
              ⁡
              
                (
                τ
                )
              
            
          
          
            
              
                λ
                
                  i
                  ,
                  
                    i
                    ⁢
                    
                      
                    
                    ⁢
                    n
                  
                
              
              ⁡
              
                (
                τ
                )
              
            
            +
            
              
                λ
                
                  i
                  ,
                  out
                
              
              ⁡
              
                (
                τ
                )
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




Refer to FIG. 2 and FIG. 4. Step S14 further comprises the following steps. In Step S22, find out the cell where the user is at the start time point, and define the cell as the start cell. Next, in Step S24, establish a plurality of movement tracks of the user moving between the start cell and the neighboring cells over a prediction interval according to the average residence time, and multiply the movement probabilities corresponding to the movement tracks to obtain a track probability. Next, in Step S26, define the cell where the user is at the end time point as the end cell, and sum up the track probabilities of all the movement tracks corresponding to the end cell to obtain the appearance probability that the user appears in the end cell.
Refer to FIG. 5, wherein Cell 3 neighbors Cell 1 and Cell 4, and wherein Cell 2 neighbors Cell 1 and Cell 4. Suppose that a unit time has fifteen minutes. Thus, there are totally 96 pieces of unit times one day, and performing statistics can obtain 96 pieces of movement probabilities pi,j(τ) and 96 pieces of average residence time Ri(τ) for the unit time 0:00-0:15, 0:15-0:30, . . . 23:45-24:00. Suppose a user enters Cell 1 at 0:01. Cell 1 functions as the start cell herein. Suppose the residence time R1(0:01) in Cell 1 is 20 minutes in the interval 0:00-0:15. It is 0:21 20 minutes later. During the unit time 0:15-0:30, the user may go to Cell 2 or Cell 3. According to statistics, the movement probability from Cell 1 to Cell 2 p1,2(0:21) and the movement probability from Cell 1 to Cell 3 p1,3(0:21) are respectively 1⁄3 and 2⁄3 for the unit time 0:15-0:30. According to statistics, the residence time in Cell 2 R2(0:21) and the residence time in Cell 3 R3(0:21) are both 20 minutes for the unit time 0:15-0:30. It is 0:41 20 minutes later. During the unit time 0:30-0:45, the user in Cell 2 may go to Cell 1 or Cell 4. According to statistics, the movement probability from Cell 2 to Cell 1 p2,1(0:41) and the movement probability from Cell 2 to Cell 4 p2,4(0:41) are respectively 1⁄3 and 2⁄3 for the unit time 0:30-0:45. The user entering Cell 3 during the unit time 0:15-0:30 may go to Cell 1 or Cell 4. According to statistics, the movement probability from Cell 3 to Cell 1 p3,1(0:41) and the movement probability from Cell 3 to Cell 4 p3,4(0:41) are respectively 1⁄3 and 2⁄3 for the unit time 0:30-0:45. Therefore, a user entering Cell 1 at 0:01 may be at the end Cell 1 or the end Cell 4 forty minutes later. The movement tracks to the end Cell 1 include Track 1→2→1 and Track 1→3→1. The movement tracks to the end Cell 4 include Track 1→2→4 and Track 1→3→4. Thus, the appearance probability of the user in Cell 1 forty minutes later is p1,2(0:21)×p2,1(0:41)+p1,3(0:21)×p3,1(0:41)=1⁄3×1⁄3+2⁄3×1⁄3=1⁄3; the appearance probability of the user in Cell 1 forty minutes later is p1,2(0:21)×p2,4(0:41)+p1,3(0:21)×p3,4(0:41)=1⁄3×2⁄3+2⁄3×2⁄3=2⁄3.
In conclusion, the present invention can be used to establish a population movement model to predict propagation of infectious fever and traffic status, which correlate with population movement. For example, the present invention can predict the average vehicle speed on a trunk highway. Refer to FIG. 6 schematically showing several freeway sections. Base stations for the areas where a freeway passes are generally constructed along the freeway. Suppose that three sections of a freeway are respectively covered by Cell 1, Cell 2 and Cell 3. The coverage distance Di of a freeway section covered by a cell can be worked out from the map information and base station construction data or obtained via measurement along the freeway. The method of the present invention can obtain the average residence time Ri(τ) of the users staying in a cell within a unit time, whereby can be instantly worked out the moving speed Vi at which users move in a cell i during a communication period τ. The moving speed is defined as the average vehicle speed and expressed by Equation (8):




  
    
      
        
          V
          i
        
        =
        
          
            D
            i
          
          
            
              R
              i
            
            ⁡
            
              (
              τ
              )
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




In conclusion, the present invention uses the traffic volume and handover information in base stations to construct population movement model and predict the population movement-related phenomena.
The embodiments described above are only to exemplify the present invention but not to limit the scope of the present invention. Any equivalent modification or variation according to the spirit of the present invention is to be also included within the scope of the present invention.

## Claims

1. A movement predicting method, which uses at least one phone communicating in and moving between neighboring cells covered by a plurality of base stations to predict movement of users of said phone in a prediction interval, comprising steps:
obtaining traffic volume generated by said phone in each said cell and handover information generated by said phone moving between said cells during one day with said base stations and storing said traffic volume and said handover information into a Mobile Switching Center (MSC);
calculating said traffic volume and said handover information to obtain a movement probability that said user moves from one said cell to a neighboring said cell and an average residence time that said user stays in one region covered by one said cell for each communication period in one day with said MSC; and
predicting an appearance probability that said user appears in each said region at an end point of said prediction interval according to a start point of said prediction interval, said movement probability, and said average residence time with said MSC.

2. The movement predicting method according to claim 1, wherein said step of calculating said traffic volume and said handover information to obtain said movement probability and said average residence time further comprises steps performed with said MSC:
dividing said traffic volume and said handover information to obtain said traffic volume and said handover information at each unit time Δt, and performing statistics to obtain a sub-traffic volume ρi(τ) in one said cell i during each said communication period τ and a sub-handover information λi,j(τ) that said user moves from one said cell i to a neighboring said cell j during each said communication period τ;
using said sub-handover information τi,j(τ) to obtain said movement probability pi,j(τ) that said user moves from one said cell i to a neighboring said cell j during each said communication period τ and an incoming sub-handover information λi,in(τ) that said user moves into one said cell i during each said communication period τ, wherein





  
    
      
        p
        
          i
          ,
          j
        
      
      ⁡
      
        (
        τ
        )
      
    
    =
    
      
        
          
            
              λ
              
                i
                ,
                j
              
            
            ⁡
            
              (
              τ
              )
            
          
          
            
              ∑
              
                j
                ,
                
                  j
                  ≠
                  i
                
              
            
            ⁢
            
              
                λ
                
                  i
                  ,
                  j
                
              
              ⁡
              
                (
                τ
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        
          
            λ
            
              i
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
          ⁡
          
            (
            τ
            )
          
        
      
      =
      
        
          ∑
          
            j
            ,
            
              j
              ≠
              i
            
          
        
        ⁢
        
          
            λ
            
              j
              ,
              i
            
          
          ⁡
          
            (
            τ
            )
          
        
      
    
  
  ;




and
using said sub-traffic volume ρi(τ) and said incoming sub-handover information λi,in(τ) to obtain said average residence time.

3. The movement predicting method according to claim 2, wherein said average residence time is calculated with an equation:




  
    
      R
      i
    
    ⁡
    
      (
      τ
      )
    
  
  =
  
    
      
        
          ρ
          i
        
        ⁡
        
          (
          τ
          )
        
      
      
        
          λ
          
            i
            ,
            
              i
              ⁢
              
                
              
              ⁢
              n
            
          
        
        ⁡
        
          (
          τ
          )
        
      
    
    .

4. The movement predicting method according to claim 2, wherein an outgoing sub-handover information λi,out(τ) that said user moves out of one said cell i during each communication period τ is also obtained in said step of using said sub-handover information λi,j(τ) to obtain said movement probability and said incoming sub-handover information, and wherein




  
    
      λ
      
        i
        ,
        out
      
    
    ⁡
    
      (
      τ
      )
    
  
  =
  
    
      ∑
      
        j
        ,
        
          j
          ≠
          i
        
      
    
    ⁢
    
      
        
          λ
          
            i
            ,
            j
          
        
        ⁡
        
          (
          τ
          )
        
      
      .

5. The movement predicting method according to claim 4, wherein said average residence time is calculated according to said sub-traffic volume ρi(τ), said incoming sub-handover information λi,in(τ), and said outgoing sub-handover information λi,out(τ).

6. The movement predicting method according to claim 5, wherein said average residence time is calculated with an equation:




  
    
      R
      i
    
    ⁡
    
      (
      τ
      )
    
  
  =
  
    
      
        2
        ⁢
        
          
        
        ⁢
        
          
            ρ
            i
          
          ⁡
          
            (
            τ
            )
          
        
      
      
        
          
            λ
            
              i
              ,
              
                i
                ⁢
                
                  
                
                ⁢
                n
              
            
          
          ⁡
          
            (
            τ
            )
          
        
        +
        
          
            λ
            
              i
              ,
              out
            
          
          ⁡
          
            (
            τ
            )
          
        
      
    
    .

7. The movement predicting method according to claim 1, wherein said step of predicting said appearance probability at said end point according to said start point, said movement probability and said average residence time further comprises steps performed with said MSC:
finding out said cell used as a start cell where said user is at said start time point;
establishing a plurality of movement tracks of said user moving between said start cell and neighboring said cells during said prediction interval according to said average residence time, and multiplying said movement probabilities corresponding to said movement tracks to obtain track probabilities; and
defining said cell where said user is at said end time point as an end cell, and summing up said track probabilities of all said movement tracks corresponding to said end cell to obtain said appearance probability that said user appears in said end cell.

8. The movement predicting method according to claim 1, wherein a moving speed of said user in one said cell i during said communication period τ is




  
    
      D
      i
    
    
      
        R
        i
      
      ⁡
      
        (
        τ
        )
      
    
  
  ,




and wherein Di is a coverage distance covered by said cell i, and wherein Ri(τ) is said average residence time that said user stays in said cell i.