The Scientific World Journal

Volume 2013, Article ID 187479, 10 pages

http://dx.doi.org/10.1155/2013/187479

## Bayesian Estimation-Based Pedestrian Tracking in Microcells

^{1}The Cybermedia Center, Osaka University, Toyonaka 560-0043, Japan^{2}Graduate School of Engineering, Osaka University, Suita 565-0871, Japan^{3}Graduate School of Information Science and Technology, Osaka University, Suita 565-0871, Japan

Received 8 July 2013; Accepted 21 August 2013

Academic Editors: S. H. Rubin and S. Sun

Copyright © 2013 Yoshiaki Taniguchi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We consider a pedestrian tracking system where sensor nodes are placed only at specific points so that the monitoring region is divided into multiple smaller regions referred to as microcells. In the proposed pedestrian tracking system, sensor nodes composed of pairs of binary sensors can detect pedestrian arrival and departure events. In this paper, we focus on pedestrian tracking in microcells. First, we investigate actual pedestrian trajectories in a microcell on the basis of observations using video sequences, after which we prepare a pedestrian mobility model. Next, we propose a method for pedestrian tracking in microcells based on the developed pedestrian mobility model. In the proposed method, we extend the Bayesian estimation to account for time-series information to estimate the correspondence between pedestrian arrival and departure events. Through simulations, we show that the tracking success ratio of the proposed method is increased by 35.8% compared to a combinatorial optimization-based tracking method.

#### 1. Introduction

Wireless sensor network technology has attracted considerable attention in recent years. Pedestrian tracking is one of the most promising applications of wireless sensor networks. For example, the ability to evaluate the trajectories of pedestrians in shopping malls or event sites would allow vendors and event organizers to arrange goods or exhibitions more effectively so that they attract the attention of more people. A variety of sensors have been considered for implementation of pedestrian tracking by using wireless sensor networks, such as cameras [1, 2], laser range scanners [3], and binary sensors [4–12].

Binary sensors, such as infrared and pressure sensors, are among the most simple, inexpensive, and energy-efficient sensors. Even though they can detect only the presence or absence of pedestrians within the sensing range, they have found widespread adoption for various applications. A number of pedestrian tracking methods using a binary sensor network have been proposed [4–12]. However, the above-mentioned studies have relied on the assumption that sensor nodes are deployed uniformly so that the sensing region covers the entire monitoring region, which requires a large number of sensor nodes.

In this paper, we consider a pedestrian tracking system where sensor nodes are placed only at specific points in the monitoring region, so that the monitoring region is divided into multiple smaller regions referred to as *microcells*. Each sensor node is composed of a pair of binary sensors coupled with a wireless communication device, and its function is to detect pedestrian arrival and departure events [13–15]. We refer to the locations at which sensor nodes are placed as *gates*. Figure 1 shows a schematic representation of the proposed tracking system implemented in a university campus. In the tracking system, sensor information is collected by a tracking server through a wireless network, and the tracking server estimates pedestrian trajectories based on sensor information.

In this paper, we focus on pedestrian tracking in a single microcell. First, we investigate actual pedestrian trajectories in a microcell based on video sequences, after which we develop a pedestrian mobility model. Next, we show that pedestrians move approximately in a straight line and that pedestrian velocities in a microcell follow a normal distribution. These results suggest that we can estimate pedestrian trajectories in a microcell by observing pedestrian arrival and departure events at specific points (i.e., gates) and estimating the correspondence between these events. Next, we propose a novel pedestrian tracking method based on the developed pedestrian mobility model. In the proposed method, we extend the Bayesian estimation to account for time-series information to estimate the correspondence between pedestrian arrival and departure events. We evaluate the performance of proposed method in a variety of situations by using actual and artificial trajectories through simulation experiments.

The rest of this paper is organized as follows. In Section 2, we present the observation results of actual pedestrian trajectories and develop a pedestrian mobility model in a microcell. In Section 3, we propose a Bayesian estimation-based pedestrian tracking method. We evaluate the performance of proposed method in Section 4. Finally, we conclude this paper with possible directions of future research in Section 5.

#### 2. Modeling Pedestrian Mobility in a Microcell

In this section, we investigate actual pedestrian trajectories in a microcell by analyzing video sequences and develop a pedestrian mobility model.

##### 2.1. Observation Overview

We took a video of actual pedestrians walking on a road to investigate actual pedestrian trajectories in a microcell. The video was taken from the 7th floor of the building of the Cybermedia Center, Osaka University, Japan. Figure 2 shows a sketch and a snapshot of the observation environment, which corresponds to the gray microcell and its surrounding environment in Figure 1. In the observation environment, there are a road and two buildings on both sides: the Cybermedia Center at the near end, and the dining hall is at the opposite side. To go to the dining hall from the road, pedestrians have to move to avoid the guard rails. The duration of the recorded video was 600 s, and the sampling rate was 500 ms. In the video sequences, a 6 6 m square region was chosen as a microcell (indicated by the quadrangular region enclosed by red lines in Figure 2). For the duration of the video, 125 pedestrians passed through the microcell. Therefore, the overall arrival rate was 125/600 = 0.208 pedestrians/s. Below, we refer to this rate as *cell arrival rate* and denote it as . We manually obtained the locations of all pedestrians in all frames. Then, we obtained the trajectory of each pedestrian by interlinking the pedestrians’ locations in ascending order in the observation period. To convert pedestrian trajectories from the camera coordinate system to a world coordinate system in overhead view, we used the projective transformation method [16]. We also obtained pedestrian velocities between two successive locations for each pedestrian.

##### 2.2. Pedestrian Trajectories in a Microcell

Figure 3(a) shows the actual pedestrian trajectories, while Figure 3(b) shows linearly approximated trajectories drawn from the arrival location of each pedestrian to the corresponding departure location. Figure 4 shows the cumulative distribution function of the distance error obtained by comparing the actual locations with the corresponding approximated locations in the sampling interval. As shown in Figure 4, the distance error is within 0.30 m in 90% of the cases. Since 0.30 m is a relatively small error compared to the size of pedestrians, the linear approximation is considered sufficient for obtaining pedestrian trajectories in a microcell. These results indicate that we can estimate pedestrian trajectories by obtaining only the pedestrian arrival and departure locations.

##### 2.3. Pedestrian Velocities in a Microcell

Figure 5 shows the cumulative distribution function of the pedestrian velocities obtained in the sampling interval. From the observation results, the mean pedestrian velocity and the variance were calculated to be 1.35 m/s and m^{2}/s^{2}, respectively. In the same figure, we also show a normal distribution [17] which approximates the observation results. The experimentally observed pedestrian velocities are found to be well approximated by the normal distribution. There are also other distributions which share similar trends with the observed cumulative distribution function, such as the log-normal distribution [18], the gamma distribution [19], and the Cauchy distribution [20]. We fitted these distributions to the observed cumulative distribution function and calculated the mean square error for each case, in which the normal distribution yielded the smallest error. Therefore, we considered that pedestrian velocities in a microcell follow a normal distribution.

We use the normal distribution to model pedestrian velocities. Therefore, the probability density function of pedestrian velocities is given by where is pedestrian velocity. Therefore, the probability density function of the pedestrian transit time required for a pedestrian to cover a distance can be written as follows: where is the pedestrian transit time. The details of the derivation of (2) from (1) are presented in the appendix.

##### 2.4. Pedestrian Gate Transition and Gate Arrival Probabilities

To investigate the relationships between the pedestrian arrival and departure locations, we divide the microcell into equal-sized grids and obtain segments (gates) at the boundary. Below, the th gate is denoted as , and the set of gates in the microcell is denoted as . Figure 6 shows the relationship between microcell and gates for . If a pedestrian passes through the microcell as shown in Figure 6, the arrival and departure events are detected at gates and , respectively.

We refer to the probability that a pedestrian arrives at gate and departs from gate as *gate transition probability* and denote it by . The gate transition probability is calculated from observation as follows:
where indicates the number of pedestrians who arrive at gate and depart from gate and indicates the total number of pedestrians in the video sequences. The set of gate transition probabilities is denoted by . In addition, we refer to the probability that a pedestrian arrives at gate as *gate arrival probability* and denote it by . The gate arrival probability is calculated from observation as follows:
where indicates the number of pedestrians who arrive at gate . The set of gate arrival probabilities is denoted by .

Figures 7 and 8, respectively, show the distributions of gate transition and gate arrival probabilities when is set to 80, that is, when the width of each gate is 0.3 m. The gate transition probabilities are symmetrical with respect to the line , and most gate transition probabilities are zero. These results indicate that there are few bidirectional paths which are frequently used and that the gate transition probabilities do not follow a uniform distribution in real-world environments. From Figure 8, we can see that the distribution of gate arrival probabilities is also nonuniform. The gate arrival probability is zero at 40% of the gates and takes a high value at some adjacent gates. Table 1 summarizes the parameters and their values in the observed environment.

##### 2.5. Summary of Pedestrian Mobility in a Microcell

Through these experiments we found that (1) pedestrians move approximately along straight lines, (2) pedestrian velocities follow a normal distribution, and (3) gate transition and gate arrival probabilities do not follow a normal distribution. We conjecture that (1) stems from the fact that pedestrians aim to reach their destination along the shortest path. Furthermore, we assume that (2) is a general characteristic of pedestrian movement. Finally, (3) appears to depend on the environment, for example, the locations and configuration of facilities. Although our experiment provides only a single example, we expect that this new knowledge will be applicable in other similar environments.

#### 3. Proposed Pedestrian Tracking Method

Considering the model presented in the previous section, in this section we propose a Bayesian estimation-based pedestrian tracking method.

##### 3.1. Overview

First, we present an overview of the tracking system. In this paper, we focus on pedestrian tracking in a microcell constituting a small part of the entire monitoring region (Figure 1). We should note here that the microcell is assumed to be sufficiently small area so that pedestrian trajectories can be approximated linearly. The border of the microcell is divided into segments (gates), and a sensor node is placed at each gate. Each sensor node is composed of a pair of binary sensors with a wireless communication device, and its function is to detect pedestrian arrival and departure events. Here, note that we consider that an arrival/departure event is detected at the location of the sensor node. We denote an arrival event detected by the sensor node at gate at time by and a departure event detected by the sensor node at gate at time by . Sensor information about arrival/departure events and the timing of events is collected by a tracking server, which estimates the pedestrian trajectories based on this information.

Since each arrival and departure event is observed independently at each gate, matching arrival and departure events are required. In this study, matching is performed when the tracking server obtains information on a departure event. When such an event is detected, the tracking server must have more than one arrival event as a candidate for matching the departure event. To select the appropriate arrival event from a set of candidate arrival events, we propose a matching method based on the Bayesian estimation [21] which is extended to account for time-series information.

##### 3.2. Matching Likelihood

Before presenting the proposed tracking method, we first derive the likelihood that departure event corresponds to arrival event . We refer to this as *matching likelihood* and denote it by . We can calculate the matching likelihood by using the following theorem.

Theorem 1. *When the tracking system is in a steady state, the matching likelihood is written as follows:
**
where is the distance between gates and .*

*Proof. *Let be an arrival event detected at gate and a departure event detected at gate . According to the Bayes theorem, we can obtain the relationship between the conditional and marginal probabilities of stochastic events and as
By taking time into account, (6) can be extended as follows:
In a steady state, we assume that the occurrence rate of arrival and departure events is independent of time. As a result, (7) becomes
Let be the velocity of a pedestrian passing between two gates. Also in (8) can be denoted by the product of the probability distribution of the transit time between gates and and the probability of pedestrian transit from gate to gate in the steady state. Consequently, (8) can be rewritten as follows by using (6):
In (9), denotes the probability that a pedestrian moves from to in the steady state, that is, gate transition probability . Therefore, we obtain the matching likelihood as in (5).

##### 3.3. Bayesian Estimation-Based Tracking Method

We now propose a tracking method using the results described in the previous section. We assume that the distribution of gate transition probabilities , mean velocity , and variance of velocity is estimated a priori by pre-learning.

The tracking server maintains a set of candidate arrival events . When the tracking server obtains information about an arrival event, it adds the arrival event to the set of candidate arrival events for future matching. Conversely, when the tracking server obtains information about a departure event , it starts matching the departure event to arrival events in the set of candidate arrival events . Here, we denote the th candidate arrival event in the set by . The tracking server first calculates the matching likelihood for each candidate arrival event using (5). Next, it selects the candidate arrival event with the highest matching likelihood as the arrival event corresponding to departure event : After matching, the tracking server estimates a line from the location of gate to the location of gate as the pedestrian trajectory.

Since there is a limitation on the transit time in a microcell, on the tracking server does not need to contain all arrival events. In the proposed method, candidate arrival events which have occurred prior to a certain period of time are removed from . In addition, the candidate arrival event is removed from after matching depending on its *matching reliability *, which is defined as follows:
After matching, if the matching reliability is higher than a predefined reliability threshold, that is, , candidate arrival event is removed from .

Figure 9 shows an example of pedestrian tracking in a microcell. In this example, a pedestrian arrives at the microcell at gate at time , after which another pedestrian arrives at gate at time . Therefore, the tracking server maintains a set of two candidate arrival events: . After a certain amount of time, a pedestrian departs from the microcell at gate at time . At this stage, there are two possible pedestrian trajectories to gate : trajectory 1 from gate or trajectory 2 from gate . Among these candidate trajectories, the matching likelihood for the trajectory from gate is higher than that from gate . Therefore, trajectory 1 is selected as the tracking result. If the matching reliability is higher than the reliability threshold, candidate arrival event is removed from .

#### 4. Performance Evaluation

In this section, we evaluate the proposed method. Although our method utilizes linearly approximated pedestrian trajectories, which are slightly different from actual pedestrian trajectories, the difference is negligible, as shown in Section 2. Therefore, in this section we focus on the estimation accuracy of the correspondence between arrival and departure events to evaluate the performance of the proposed method. We define *tracking success ratio* as the ratio of the number of successful matches to the total number of matches and use this ratio as a performance evaluation index.

##### 4.1. Comparison with a Combinatorial Optimization-Based Method Using Actual Pedestrian Trajectories

First, we evaluate the proposed method using actual pedestrian trajectories obtained in Section 2. To evaluate the basic performance of the proposed method, we use the parameters and their values listed in Table 1. The results are compared with those obtained using a combinatorial optimization-based method [22] (referred to as the compared method below), which uses the mean pedestrian velocity for matching arrival and departure events. In the compared method, the objective function is written as where is a set of arrival and departure events in a specific time period. The combination of pairs of arrival and departure events for which the objective function is minimized is selected as the matching result. Although the accuracy of matching increases if the distribution of velocities is taken into account, the compared method has a larger search space and higher complexity.

We calculate the tracking success ratio for both the proposed method and the compared method using the data set of actual pedestrian trajectories. The tracking success ratio of the proposed method is 1, while that of the compared method reaches a maximum of 0.736, amounting to an improvement of 35.8% in tracking success ratio with the proposed method. It should, however, be noted that the tracking success ratio in this evaluation reaches 1 since the cell arrival rate in the actual data set is relatively low, and we assume that the parameters of the proposed method are known. In the following sections, we evaluate the proposed method in a variety of situations.

##### 4.2. Effect of Reliability Threshold

In this section, we present a simulation using artificial pedestrian trajectories to investigate the performance of the proposed method in a variety of situations. In the simulation, the size of the microcell is set to m and the number of gates is set to 80. A pedestrian arrives at gate in accordance with a Poisson process with a gate arrival rate . The gate arrival rate is the product of the cell arrival rate and gate arrival probability . The departure gate of the pedestrian is determined according to gate transition probabilities , where we assume that the pedestrian moves to the departure gate along a straight line. The pedestrian velocities follow a normal distribution. In this simulation, we use the parameter settings shown in Table 1 for the mean velocity , variance of velocity , gate transition probabilities , and gate-arrival probability . The observation duration is set to 600 s, and the cell arrival rate is set to 1, 2, or 3 pedestrians/s. The correspondence between arrival and departure events is estimated using the proposed method. The following results are averaged over 100 iterations.

Figure 10 shows a plot of the relation between tracking success ratio and the reliability threshold , where we can see that there is an optimal value for the reliability threshold which maximizes the tracking success ratio. The optimal reliability threshold was between 0.8 and 0.9 in this evaluation. The reason why there is an optimal value is as follows. When the reliability threshold is too low, candidate arrival events are easily deleted from the set of candidate arrival events even in the case of a mismatch. As a result, the number of mismatches increases and the tracking success ratio decreases. In contrast, when the reliability threshold is too high, more candidate arrival events remain in the set of candidate arrival events, increasing the possibility for a mismatch.

In addition, as shown in Figure 10, the tracking success ratio is higher for lower cell arrival rates. This is because the number of candidate arrival events for matching a departure event is higher for higher cell arrival rates, and as a result it becomes difficult to estimate the correct arrival event.

##### 4.3. Effect of Distribution of Gate Transition Probabilities

Then, we investigate how the distribution of gate transition probabilities affects the tracking success ratio. We used two types of gate transition probability distributions: a uniform distribution and an observation-based distribution (Figure 7). In this section, the reliability threshold is set to 0.9.

Figure 11 shows the relationship between the tracking success ratio, the pedestrian cell arrival rate, and the distribution of gate transition probabilities. The tracking success ratio is lower for higher cell arrival rates, and this trend is independent of the distribution of transition probabilities. However, when we used the observed distribution of gate transition probabilities, the tracking success ratio remains much higher for different cell arrival rates as compared to the results obtained with the uniform distribution of gate transition probabilities. From this result, we conclude that the proposed method is more suitable for situations where pedestrian transition is nonuniform.

##### 4.4. Effect of Learning Period

In the previous sections, we used the same values for the parameters in the mobility model in the simulation and for the parameters in the proposed method. In practical situations, the parameters of the proposed method are obtained by prelearning. In this section, we investigate the effect of the duration of the learning period on the tracking success ratio of the proposed method. In this simulation, the observation period is set to 3600 s and the cell arrival rate is set to 2 pedestrians/s. We calculate the distribution of gate transit probabilities using the beginning (i.e., the learning period) of the observation period and then evaluate the tracking success ratio of the entire observation period.

Figure 12 shows the tracking success ratio plotted against the learning period. To show that there is an upper limit to the tracking success ratio, we also show the tracking success ratio when the entire observed period of 3600 s is used as a learning period (Figure 12). The tracking success ratio clearly increases together with the learning period since with longer learning periods we can obtain more precise parameter values. When the learning period is set to 600 s, the tracking success ratio is at the upper limit. In addition, the tracking success ratio is almost at the upper limit when the learning period is set to 300 s. Therefore, the learning period can be set to a value between these values.

#### 5. Conclusions and Future Work

In this paper, first we presented a model of pedestrian mobility in a microcell on the basis of observation of actual pedestrian trajectories. We showed that pedestrians move along approximately straight lines and that pedestrian velocities follow a normal distribution. Based on these results, we proposed a novel method for pedestrian tracking in a microcell. In the proposed method, we extend Bayesian estimation to account for time-series information in order to estimate the correspondence between pedestrian arrival and departure events. Through simulations, we evaluated the basic performance of the proposed method in a variety of situations and demonstrated that the tracking success ratio of the proposed method is improved by 35.8% compared to a combinatorial optimization-based tracking method.

In this paper, we focused on pedestrian tracking in a microcell. In future work, we plan to develop an inter-microcell pedestrian tracking method by using results obtained in this paper. We also plan to evaluate the proposed method using actual pedestrian trajectories under different environments in terms of location, season, and time in order to investigate the feasibility of the method in real-world environments.

#### Appendix

Theorem A.1. *Assume that pedestrian velocity follows a normal distribution with a mean value and variance . The probability density function of the distribution of pedestrian transit time necessary for a pedestrian to cover a distance can be written as*

*Proof. *Since pedestrian velocity follows a normal distribution, the probability density function of the distribution of pedestrian velocities can be written as

Let and be the cumulative distribution functions of (A.1) and (A.2), respectively, where and are random variables. The following relationships exist among the probability density functions and the cumulative distribution functions:
Here, can be written as follows:
By differentiating (A.4) with respect to the transit time , we obtain
By substituting (A.2) into (A.5), we obtain (A.1).

#### Acknowledgment

This work was partly supported by the KAKENHI 25330103 of the JSPS in Japan.

#### References

- A. Yilmaz, O. Javed, and M. Shah, “Object tracking: a survey,”
*ACM Computing Surveys*, vol. 38, no. 4, 2006. View at Publisher · View at Google Scholar · View at Scopus - K. Heath and L. Guibas, “Multi-person tracking from sparse 3D trajectories in a camera sensor network,” in
*Proceedings of the 2nd ACM/IEEE International Conference on Distributed Smart Cameras (ICDSC '08)*, September 2008. View at Publisher · View at Google Scholar · View at Scopus - S. Gidel, P. Checchin, C. Blanc, T. Chateau, and L. Trassoudaine, “Pedestrian detection and tracking in an urban environment using a multilayer laser scanner,”
*IEEE Transactions on Intelligent Transportation Systems*, vol. 11, no. 3, pp. 579–588, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Schiff and K. Goldberg, “Automated intruder tracking using particle filtering and a network of binary motion sensors,” in
*Proceedings of the IEEE International Conference on Automation Science and Engineering (CASE '06)*, pp. 580–587, October 2006. View at Publisher · View at Google Scholar · View at Scopus - H. Morishita, R. Fukui, and T. Sato, “High resolution pressure sensor distributed floor for future human-robot symbiosis environments,” in
*Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IRDS '02)*, pp. 1246–1251, October 2002. View at Scopus - T. Murakita, T. Ikeda, and H. Ishiguro, “Human tracking using floor sensors based on the markov chain monte carlo method,” in
*Proceedings of the 17th International Conference on Pattern Recognition (ICPR '04)*, pp. 917–920, August 2004. View at Publisher · View at Google Scholar · View at Scopus - J. Aslam, Z. Butler, F. Constantin, V. Crespi, G. Cybenko, and D. Rus, “Tracking a moving object with a binary sensor network,” in
*Proceedings of the 1st International Conference on Embedded Networked Sensor Systems (SenSys '03)*, pp. 150–161, November 2003. View at Scopus - L. Lazos, R. Poovendran, and J. A. Ritcey, “Probabilistic detection of mobile targets in heterogeneous sensor networks,” in
*Proceedings of the 6th International Symposium on Information Processing in Sensor Networks (IPSN '07)*, pp. 519–528, April 2007. View at Publisher · View at Google Scholar · View at Scopus - X. Liu, G. Zhao, and X. Ma, “Target localization and tracking in noisy binary sensor networks with known spatial topology,” in
*Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP '07)*, pp. 1029–1032, April 2007. View at Publisher · View at Google Scholar · View at Scopus - Z. Wang, E. Bulut, and B. K. Szymanski, “Distributed energy-efficient target tracking with binary sensor networks,”
*ACM Transactions on Sensor Networks*, vol. 6, no. 4, pp. 1–32, 2010. View at Publisher · View at Google Scholar · View at Scopus - Y. Busnel, L. Querzoni, R. Baldoni, M. Bertier, and A. Kermarrec, “Analysis of deterministic tracking of multiple objects using a binary sensor network,”
*ACM Transactions on Sensor Networks*, vol. 8, no. 1, pp. 1–27, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Singh, R. Kumar, U. Madhow, S. Suri, and R. Cagley, “Multiple-target tracking with binary proximity sensors,”
*ACM Transactions on Sensor Networks*, vol. 8, no. 1, pp. 1–26, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Fujii, Y. Taniguchi, G. Hasegawa, and H. Nakano, “Monte Carlobased bi-directional pedestrian counting method for compound-eye sensor systems,”
*Journal of Emerging Trends in Computing and Information Sciences*, vol. 4, pp. 48–57, 2013. View at Google Scholar - “Pyroelectric sensor,” http://www.eco-compteur.com/?wpid=15036.
- “Acoustic slab sensor,” http://www.eco-compteur.com/?wpid=15035.
- G. Sparr, “Projective invariants for affine shapes of point configurations,” in
*Proceeding of the DARPA-ESPRIT Workshop on Applications of Invariants in Computer Vision*, pp. 151–170, March 1991. - J. K. Patel and C. B. Read,
*Handbook of the Normal Distribution*, Marcel Dekker, New York, NY, USA, 1996. - C. C. Heyde, “On a property of the lognormal distribution,”
*Journal of the Royal Statistical Society*, vol. 25, pp. 392–393, 1963. View at Google Scholar - H. C. S. Thom, “A note on the gamma distribution,”
*Monthly Weather Review*, vol. 86, pp. 117–122, 1958. View at Google Scholar - M. V. Menon, “A characterization of the Cauchy distribution,”
*The Annals of Mathematical Statistics*, vol. 33, pp. 1267–1271, 1962. View at Google Scholar - J. M. Bernardo and A. F. M. Smith,
*Bayesian Theory*, Wiley & Sons, New York, NY, USA, 2000. - S. Aihara, M. Sasabe, and H. Nakano, “Mobility model based on incoming and outgoing nodes to an area,” in
*Proceedings of the Mobile Mapping Technologies Symposium (MMT '07)*, May 2007.