Poisson point process

Poisson distribution of point counts

A Poisson point process is characterized via the Poisson distribution. The Poisson distribution is the probability distribution of a random variable ${\textstyle N}$ (called a Poisson random variable) such that the probability that ${\displaystyle \textstyle N}$ equals ${\displaystyle \textstyle n}$ is given by:

${\displaystyle \Pr\{N=n\}={\frac {\Lambda ^{n}}{n!}}e^{-\Lambda }}$

where ${\textstyle n!}$ denotes factorial and the parameter ${\textstyle \Lambda }$ determines the shape of the distribution. (In fact, ${\textstyle \Lambda }$ equals the expected value of ${\textstyle N}$.)

By definition, a Poisson point process has the property that the number of points in a bounded region of the process's underlying space is a Poisson-distributed random variable.^[36]

Complete independence

Consider a collection of disjoint and bounded subregions of the underlying space. By definition, the number of points of a Poisson point process in each bounded subregion will be completely independent of all the others.

This property is known under several names such as complete randomness, complete independence,^[38] or independent scattering^[39][40] and is common to all Poisson point processes. In other words, there is a lack of interaction between different regions and the points in general,^[41] which motivates the Poisson process being sometimes called a purely or completely random process.^[38]

Interpreted as a counting process

The homogeneous Poisson point process, when considered on the positive half-line, can be defined as a counting process, a type of stochastic process, which can be denoted as ${\textstyle \{N(t),t\geq 0\}}$.^[29][32] A counting process represents the total number of occurrences or events that have happened up to and including time ${\textstyle t}$. A counting process is a homogeneous Poisson counting process with rate ${\textstyle \lambda >0}$ if it has the following three properties:^[29][32]

• ${\textstyle N(0)=0;}$
• has independent increments; and
• the number of events (or points) in any interval of length ${\textstyle t}$ is a Poisson random variable with parameter (or mean) ${\textstyle \lambda t}$.

The last property implies:

${\displaystyle \operatorname {E} [N(t)]=\lambda t.}$

In other words, the probability of the random variable ${\textstyle N(t)}$ being equal to ${\textstyle n}$ is given by:

${\displaystyle \Pr\{N(t)=n\}={\frac {(\lambda t)^{n}}{n!}}e^{-\lambda t}.}$

The Poisson counting process can also be defined by stating that the time differences between events of the counting process are exponential variables with mean ${\textstyle 1/\lambda }$.^[45] The time differences between the events or arrivals are known as interarrival ^[46] or interoccurrence times.^[45]

Interpreted as a point process on the real line

Interpreted as a point process, a Poisson point process can be defined on the real line by considering the number of points of the process in the interval ${\textstyle (a,b]}$. For the homogeneous Poisson point process on the real line with parameter ${\textstyle \lambda >0}$, the probability of this random number of points, written here as ${\textstyle N(a,b]}$, being equal to some counting number ${\textstyle n}$ is given by:^[47]

${\displaystyle \Pr\{N(a,b]=n\}={\frac {[\lambda (b-a)]^{n}}{n!}}e^{-\lambda (b-a)},}$

For some positive integer ${\textstyle k}$, the homogeneous Poisson point process has the finite-dimensional distribution given by:^[47]

${\displaystyle \Pr\{N(a_{i},b_{i}]=n_{i},i=1,\dots ,k\}=\prod _{i=1}^{k}{\frac {[\lambda (b_{i}-a_{i})]^{n_{i}}}{n_{i}!}}e^{-\lambda (b_{i}-a_{i})},}$

where the real numbers ${\textstyle a_{i}<b_{i}\leq a_{i+1}}$.

In other words, ${\textstyle N(a,b]}$ is a Poisson random variable with mean ${\textstyle \lambda (b-a)}$, where ${\textstyle a\leq b}$. Furthermore, the number of points in any two disjoint intervals, say, ${\textstyle (a_{1},b_{1}]}$ and ${\textstyle (a_{2},b_{2}]}$ are independent of each other, and this extends to any finite number of disjoint intervals.^[47] In the queueing theory context, one can consider a point existing (in an interval) as an event, but this is different to the word event in the probability theory sense.^[b] It follows that ${\textstyle \lambda }$ is the expected number of arrivals that occur per unit of time.^[32]

Spatial Poisson point process

A spatial Poisson process is a Poisson point process defined in the plane ${\displaystyle \textstyle \mathbb {R} ^{2}}$.^[57][64] For its mathematical definition, one first considers a bounded, open or closed (or more precisely, Borel measurable) region ${\textstyle B}$ of the plane. The number of points of a point process ${\displaystyle \textstyle N}$ existing in this region ${\displaystyle \textstyle B\subset \mathbb {R} ^{2}}$ is a random variable, denoted by ${\displaystyle \textstyle N(B)}$. If the points belong to a homogeneous Poisson process with parameter ${\displaystyle \textstyle \lambda >0}$, then the probability of ${\displaystyle \textstyle n}$ points existing in ${\displaystyle \textstyle B}$ is given by:

${\displaystyle \Pr\{N(B)=n\}={\frac {(\lambda |B|)^{n}}{n!}}e^{-\lambda |B|}}$

where ${\displaystyle \textstyle |B|}$ denotes the area of ${\displaystyle \textstyle B}$.

For some finite integer ${\displaystyle \textstyle k\geq 1}$, we can give the finite-dimensional distribution of the homogeneous Poisson point process by first considering a collection of disjoint, bounded Borel (measurable) sets ${\displaystyle \textstyle B_{1},\dots ,B_{k}}$. The number of points of the point process ${\displaystyle \textstyle N}$ existing in ${\displaystyle \textstyle B_{i}}$ can be written as ${\displaystyle \textstyle N(B_{i})}$. Then the homogeneous Poisson point process with parameter ${\displaystyle \textstyle \lambda >0}$ has the finite-dimensional distribution:^[65]

${\displaystyle \Pr\{N(B_{i})=n_{i},i=1,\dots ,k\}=\prod _{i=1}^{k}{\frac {(\lambda |B_{i}|)^{n_{i}}}{n_{i}!}}e^{-\lambda |B_{i}|}.}$

Applications

The spatial Poisson point process features prominently in spatial statistics,^[21][22] stochastic geometry, and continuum percolation theory.^[23] This point process is applied in various physical sciences such as a model developed for alpha particles being detected. In recent years, it has been frequently used to model seemingly disordered spatial configurations of certain wireless communication networks.^[17][18][19] For example, models for cellular or mobile phone networks have been developed where it is assumed the phone network transmitters, known as base stations, are positioned according to a homogeneous Poisson point process.

Defined in higher dimensions

The previous homogeneous Poisson point process immediately extends to higher dimensions by replacing the notion of area with (high dimensional) volume. For some bounded region ${\displaystyle \textstyle B}$ of Euclidean space ${\displaystyle \textstyle \mathbb {R} ^{d}}$, if the points form a homogeneous Poisson process with parameter ${\displaystyle \textstyle \lambda >0}$, then the probability of ${\displaystyle \textstyle n}$ points existing in ${\displaystyle \textstyle B\subset \mathbb {R} ^{d}}$ is given by:

${\displaystyle \Pr\{N(B)=n\}={\frac {(\lambda |B|)^{n}}{n!}}e^{-\lambda |B|}}$

where ${\displaystyle \textstyle |B|}$ now denotes the ${\displaystyle \textstyle d}$-dimensional volume of ${\displaystyle \textstyle B}$. Furthermore, for a collection of disjoint, bounded Borel sets ${\displaystyle \textstyle B_{1},\dots ,B_{k}\subset \mathbb {R} ^{d}}$, let ${\displaystyle \textstyle N(B_{i})}$ denote the number of points of ${\displaystyle \textstyle N}$ existing in ${\displaystyle \textstyle B_{i}}$. Then the corresponding homogeneous Poisson point process with parameter ${\displaystyle \textstyle \lambda >0}$ has the finite-dimensional distribution:^[67]

${\displaystyle \Pr\{N(B_{i})=n_{i},i=1,\dots ,k\}=\prod _{i=1}^{k}{\frac {(\lambda |B_{i}|)^{n_{i}}}{n_{i}!}}e^{-\lambda |B_{i}|}.}$

Homogeneous Poisson point processes do not depend on the position of the underlying space through its parameter ${\displaystyle \textstyle \lambda }$, which implies it is both a stationary process (invariant to translation) and an isotropic (invariant to rotation) stochastic process.^[25] Similarly to the one-dimensional case, the homogeneous point process is restricted to some bounded subset of ${\textstyle \mathbb {R} ^{d}}$, then depending on some definitions of stationarity, the process is no longer stationary.^[25][52]

Points are uniformly distributed

If the homogeneous point process is defined on the real line as a mathematical model for occurrences of some phenomenon, then it has the characteristic that the positions of these occurrences or events on the real line (often interpreted as time) will be uniformly distributed. More specifically, if an event occurs (according to this process) in an interval ${\displaystyle \textstyle (a,b]}$ where ${\displaystyle \textstyle a\leq b}$, then its location will be a uniform random variable defined on that interval.^[65] Furthermore, the homogeneous point process is sometimes called the uniform Poisson point process (see Terminology). This uniformity property extends to higher dimensions in the Cartesian coordinate, but not in, for example, polar coordinates.^[68][69]

Defined on the real line

On the real line, the inhomogeneous or non-homogeneous Poisson point process has mean measure given by a one-dimensional integral. For two real numbers ${\displaystyle \textstyle a}$ and ${\displaystyle \textstyle b}$, where ${\displaystyle \textstyle a\leq b}$, denote by ${\displaystyle \textstyle N(a,b]}$ the number points of an inhomogeneous Poisson process with intensity function ${\displaystyle \textstyle \lambda (t)}$ occurring in the interval ${\displaystyle \textstyle (a,b]}$. The probability of ${\displaystyle \textstyle n}$ points existing in the above interval ${\displaystyle \textstyle (a,b]}$ is given by:

${\displaystyle \Pr\{N(a,b]=n\}={\frac {[\Lambda (a,b)]^{n}}{n!}}e^{-\Lambda (a,b)}.}$

where the mean or intensity measure is:

${\displaystyle \Lambda (a,b)=\int _{a}^{b}\lambda (t)\,\mathrm {d} t,}$

which means that the random variable ${\displaystyle \textstyle N(a,b]}$ is a Poisson random variable with mean ${\displaystyle \textstyle \operatorname {E} [N(a,b]]=\Lambda (a,b)}$.

A feature of the one-dimension setting, is that an inhomogeneous Poisson process can be transformed into a homogeneous by a monotone transformation or mapping, which is achieved with the inverse of ${\displaystyle \textstyle \Lambda }$.^[70][71]

Spatial Poisson process

An inhomogeneous Poisson process defined in the plane ${\displaystyle \textstyle \mathbb {R} ^{2}}$ is called a spatial Poisson process^[16] It is defined with intensity function and its intensity measure is obtained performing a surface integral of its intensity function over some region.^[20][74] For example, its intensity function (as a function of Cartesian coordinates ${\textstyle x}$ and ${\displaystyle \textstyle y}$) can be

${\displaystyle \lambda (x,y)=e^{-(x^{2}+y^{2})},}$

so the corresponding intensity measure is given by the surface integral

${\displaystyle \Lambda (B)=\int _{B}e^{-(x^{2}+y^{2})}\,\mathrm {d} x\,\mathrm {d} y,}$

where ${\textstyle B}$ is some bounded region in the plane ${\textstyle \mathbb {R} ^{2}}$.

In higher dimensions

In the plane, ${\textstyle \Lambda (B)}$ corresponds to a surface integral while in ${\textstyle \mathbb {R} ^{d}}$ the integral becomes a (${\textstyle d}$-dimensional) volume integral.

Applications

When the real line is interpreted as time, the inhomogeneous process is used in the fields of counting processes and in queueing theory.^[72][75] Examples of phenomena which have been represented by or appear as an inhomogeneous Poisson point process include:

• Goals being scored in a soccer game.^[76]
• Defects in a circuit board^[77]

In the plane, the Poisson point process is important in the related disciplines of stochastic geometry^[1][33] and spatial statistics.^[21][22] The intensity measure of this point process is dependent on the location of underlying space, which means it can be used to model phenomena with a density that varies over some region. In other words, the phenomena can be represented as points that have a location-dependent density.^[20] This processes has been used in various disciplines and uses include the study of salmon and sea lice in the oceans,^[78] forestry,^[6] and search problems.^[79]

Interpretation of the intensity function

The Poisson intensity function ${\textstyle \lambda (x)}$ has an interpretation, considered intuitive,^[20] with the volume element ${\textstyle \mathrm {d} x}$ in the infinitesimal sense: ${\textstyle \lambda (x)\,\mathrm {d} x}$ is the infinitesimal probability of a point of a Poisson point process existing in a region of space with volume ${\textstyle \mathrm {d} x}$ located at ${\textstyle x}$.^[20]

For example, given a homogeneous Poisson point process on the real line, the probability of finding a single point of the process in a small interval of width ${\textstyle \delta }$ is approximately ${\textstyle \lambda \delta }$. In fact, such intuition is how the Poisson point process is sometimes introduced and its distribution derived.^[80][41][81]

Simple point process

If a Poisson point process has an intensity measure that is a locally finite and diffuse (or non-atomic), then it is a simple point process. For a simple point process, the probability of a point existing at a single point or location in the underlying (state) space is either zero or one. This implies that, with probability one, no two (or more) points of a Poisson point process coincide in location in the underlying space.^[82][18][83]

Step 1: Number of points

The number of points ${\textstyle N}$ in the window, denoted here by ${\textstyle W}$, needs to be simulated, which is done by using a (pseudo)-random number generating function capable of simulating Poisson random variables.

Step 2: Positioning of points

The second stage requires randomly placing the ${\displaystyle \textstyle N}$ points in the window ${\displaystyle \textstyle W}$.

Poisson distribution

Despite its name, the Poisson point process was neither discovered nor studied by its namesake. It is cited as an example of Stigler's law of eponymy.^[2][3] The name arises from the process's inherent relation to the Poisson distribution, derived by Poisson as a limiting case of the binomial distribution.^[88] It describes the probability of the sum of ${\displaystyle \textstyle n}$ Bernoulli trials with probability ${\displaystyle \textstyle p}$, often likened to the number of heads (or tails) after ${\displaystyle \textstyle n}$ biased coin flips with the probability of a head (or tail) occurring being ${\displaystyle \textstyle p}$. For some positive constant ${\displaystyle \textstyle \Lambda >0}$, as ${\displaystyle \textstyle n}$ increases towards infinity and ${\displaystyle \textstyle p}$ decreases towards zero such that the product ${\displaystyle \textstyle np=\Lambda }$ is fixed, the Poisson distribution more closely approximates that of the binomial.^[89]

In 1841, Poisson derived the Poisson distribution by studying the binomial distribution in the limit as ${\displaystyle \textstyle p}$ goes to zero and ${\displaystyle \textstyle n}$ goes to infinity. The distribution appears only once in Poisson's work,^[90] and the result was not well known during his time. Over the following years, others used the distribution without citing Poisson, including Philipp Ludwig von Seidel and Ernst Abbe.^[91][2] At the end of the 19th century, Ladislaus Bortkiewicz revived interest in the distribution by citing Poisson and using real data on the number of deaths from horse kicks in the Prussian army.^[88][92]

Discovery

There are a number of claims for early uses or discoveries of the Poisson point process.^[2][3] For example, John Michell in 1767, a decade before Poisson was born, was interested in the probability a star being within a certain region of another star under the erroneous assumption that the stars were "scattered by mere chance", and studied an example consisting of the six brightest stars in the Pleiades, without deriving the Poisson distribution. This work inspired Simon Newcomb to study the problem and to calculate the Poisson distribution as an approximation for the binomial distribution in 1860.^[3]

At the beginning of the 20th century the Poisson process (in one dimension) would arise independently in different situations.^[2][3] In Sweden 1903, Filip Lundberg published a thesis containing work, now considered fundamental and pioneering, where he proposed to model insurance claims with a homogeneous Poisson process.^[93][94]

In Denmark A.K. Erlang derived the Poisson distribution in 1909 when developing a mathematical model for the number of incoming phone calls in a finite time interval. Erlang unaware of Poisson's earlier work and assumed that the number phone calls arriving in each interval of time were independent of each other. He then found the limiting case, which is effectively recasting the Poisson distribution as a limit of the binomial distribution.^[2]

In 1910 Ernest Rutherford and Hans Geiger published experimental results on counting alpha particles. Their experimental work had mathematical contributions from Harry Bateman, who derived Poisson probabilities as a solution to a family of differential equations, though the solution had been derived earlier, resulting in the independent discovery of the Poisson process.^[2] After this time, there were many studies and applications of the Poisson process, but its early history is complicated, which has been explained by the various applications of the process in numerous fields by biologists, ecologists, engineers and various physical scientists.^[2]

Early applications

The years after 1909 led to a number of studies and applications of the Poisson point process, however, its early history is complex, which has been explained by the various applications of the process in numerous fields by biologists, ecologists, engineers and others working in the physical sciences. The early results were published in different languages and in different settings, with no standard terminology and notation used.^[2] For example, in 1922 Swedish chemist and Nobel Laureate Theodor Svedberg proposed a model in which a spatial Poisson point process is the underlying process to study how plants are distributed in plant communities.^[95] A number of mathematicians started studying the process in the early 1930s, and important contributions were made by Andrey Kolmogorov, William Feller and Aleksandr Khinchin,^[2] among others.^[96] In the field of teletraffic engineering, mathematicians and statisticians studied and used Poisson and other point processes.^[97]

History of terms

The Swede Conny Palm in his 1943 dissertation studied the Poisson and other point processes in the one-dimensional setting by examining them in terms of the statistical or stochastic dependence between the points in time.^[98][97] In his work exists the first known recorded use of the term point processes as Punktprozesse in German.^[98][3]

It is believed^[2] that William Feller was the first in print to refer to it as the Poisson process in a 1940 paper. Although the Swede Ove Lundberg used the term Poisson process in his 1940 PhD dissertation,^[3] in which Feller was acknowledged as an influence,^[99] it has been claimed that Feller coined the term before 1940.^[89] It has been remarked that both Feller and Lundberg used the term as though it were well-known, implying it was already in spoken use by then.^[3] Feller worked from 1936 to 1939 alongside Harald Cramér at Stockholm University, where Lundberg was a PhD student under Cramér who did not use the term Poisson process in a book by him, finished in 1936, but did in subsequent editions, which his has led to the speculation that the term Poisson process was coined sometime between 1936 and 1939 at the Stockholm University.^[3]

Laplace functionals

For a Poisson point process ${\displaystyle \textstyle N}$ with intensity measure ${\displaystyle \textstyle \Lambda }$ on some space ${\displaystyle X}$, the Laplace functional is given by:^[18]

${\displaystyle L_{N}(f)=\mathbb {E} e^{-\int _{X}f(x)\,N(\mathrm {d} x)}=e^{-\int _{X}(1-e^{-f(x)})\Lambda (\mathrm {d} x)},}$

One version of Campbell's theorem involves the Laplace functional of the Poisson point process.

Probability generating functionals

The probability generating function of non-negative integer-valued random variable leads to the probability generating functional being defined analogously with respect to any non-negative bounded function ${\displaystyle \textstyle v}$ on ${\displaystyle \textstyle \mathbb {R} ^{d}}$ such that ${\displaystyle \textstyle 0\leq v(x)\leq 1}$. For a point process ${\displaystyle \textstyle {N}}$ the probability generating functional is defined as:^[122]

${\displaystyle G(v)=\operatorname {E} \left[\prod _{x\in N}v(x)\right]}$

where the product is performed for all the points in ${\textstyle N}$. If the intensity measure ${\displaystyle \textstyle \Lambda }$ of ${\displaystyle \textstyle {N}}$ is locally finite, then the ${\textstyle G}$ is well-defined for any measurable function ${\displaystyle \textstyle u}$ on ${\displaystyle \textstyle \mathbb {R} ^{d}}$. For a Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$ the generating functional is given by:

${\displaystyle G(v)=e^{-\int _{\mathbb {R} ^{d}}[1-v(x)]\,\Lambda (\mathrm {d} x)},}$

which in the homogeneous case reduces to

${\displaystyle G(v)=e^{-\lambda \int _{\mathbb {R} ^{d}}[1-v(x)]\,\mathrm {d} x}.}$

Moment measure

For a general Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$ the first moment measure is its intensity measure:^[18][19]

${\displaystyle M^{1}(B)=\Lambda (B),}$

which for a homogeneous Poisson point process with constant intensity ${\displaystyle \textstyle \lambda }$ means:

${\displaystyle M^{1}(B)=\lambda |B|,}$

where ${\displaystyle \textstyle |B|}$ is the length, area or volume (or more generally, the Lebesgue measure) of ${\displaystyle \textstyle B}$.

The Mecke equation

The Mecke equation characterizes the Poisson point process. Let ${\displaystyle \mathbb {N} _{\sigma }}$ be the space of all ${\displaystyle \sigma }$-finite measures on some general space ${\displaystyle {\mathcal {Q}}}$. A point process ${\displaystyle \eta }$ with intensity ${\displaystyle \lambda }$ on ${\displaystyle {\mathcal {Q}}}$ is a Poisson point process if and only if for all measurable functions ${\displaystyle f:{\mathcal {Q}}\times \mathbb {N} _{\sigma }\to \mathbb {R} _{+}}$ the following holds

${\displaystyle E\left[\int f(x,\eta )\eta (\mathrm {d} x)\right]=\int E\left[f(x,\eta +\delta _{x})\right]\lambda (\mathrm {d} x)}$

For further details see.^[123]

Factorial moment measure

For a general Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$ the ${\displaystyle \textstyle n}$-th factorial moment measure is given by the expression:^[124]

${\displaystyle M^{(n)}(B_{1}\times \cdots \times B_{n})=\prod _{i=1}^{n}[\Lambda (B_{i})],}$

where ${\displaystyle \textstyle \Lambda }$ is the intensity measure or first moment measure of ${\displaystyle \textstyle {N}}$, which for some Borel set ${\displaystyle \textstyle B}$ is given by

${\displaystyle \Lambda (B)=M^{1}(B)=\operatorname {E} [N(B)].}$

For a homogeneous Poisson point process the ${\displaystyle \textstyle n}$-th factorial moment measure is simply:^[18][19]

${\displaystyle M^{(n)}(B_{1}\times \cdots \times B_{n})=\lambda ^{n}\prod _{i=1}^{n}|B_{i}|,}$

where ${\displaystyle \textstyle |B_{i}|}$ is the length, area, or volume (or more generally, the Lebesgue measure) of ${\displaystyle \textstyle B_{i}}$. Furthermore, the ${\displaystyle \textstyle n}$-th factorial moment density is:^[124]

${\displaystyle \mu ^{(n)}(x_{1},\dots ,x_{n})=\lambda ^{n}.}$

Rényi’s theorem

Simple point processes are completely characterized by their void probabilities.^[126] In other words, complete information of a simple point process is captured entirely in its void probabilities, and two simple point processes have the same void probabilities if and if only if they are the same point processes. The case for Poisson process is sometimes known as Rényi’s theorem, which is named after Alfréd Rényi who discovered the result for the case of a homogeneous point process in one dimension.^[127]

In one form^[127] Rényi’s theorem says that, if ${\displaystyle \textstyle \Lambda }$ is a diffuse (or non-atomic) Radon measure on ${\displaystyle \textstyle \mathbb {R} ^{d}}$ and ${\displaystyle \textstyle N}$ is a locally finite simple point process on ${\displaystyle \textstyle \mathbb {R} ^{d}}$ such that for any set ${\displaystyle \textstyle A}$ being a finite union of rectangles there holds true:

${\displaystyle \Pr\{N(A)=0\}=v(A)=e^{-\Lambda (A)}}$,

then ${\displaystyle \textstyle N}$ is a Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$.

Thinning

For the Poisson process, the independent ${\displaystyle \textstyle p(x)}$-thinning operations results in another Poisson point process. More specifically, a ${\displaystyle \textstyle p(x)}$-thinning operation applied to a Poisson point process with intensity measure ${\displaystyle \textstyle \Lambda }$ gives a point process of removed points that is also Poisson point process ${\displaystyle \textstyle {N}_{p}}$ with intensity measure ${\displaystyle \textstyle \Lambda _{p}}$, which for a bounded Borel set ${\displaystyle \textstyle B}$ is given by:

${\displaystyle \Lambda _{p}(B)=\int _{B}p(x)\,\Lambda (\mathrm {d} x)}$

This thinning result of the Poisson point process is sometimes known as Prekopa's theorem.^[129] Furthermore, after randomly thinning a Poisson point process, the kept or remaining points also form a Poisson point process, which has the intensity measure

${\displaystyle \Lambda _{p}(B)=\int _{B}(1-p(x))\,\Lambda (\mathrm {d} x).}$

The two separate Poisson point processes formed respectively from the removed and kept points are stochastically independent of each other.^[128] In other words, if a region is known to contain ${\displaystyle \textstyle n}$ kept points (from the original Poisson point process), then this will have no influence on the random number of removed points in the same region. This ability to randomly create two independent Poisson point processes from one is sometimes known as splitting^[130][131] the Poisson point process.

Superposition

If there is a countable collection of point processes ${\displaystyle \textstyle N_{1},N_{2},\dots }$, then their superposition, or, in set theory language, their union, which is^[132]

${\displaystyle N=\bigcup _{i=1}^{\infty }N_{i},}$

also forms a point process. In other words, any points located in any of the point processes ${\displaystyle \textstyle N_{1},N_{2}\dots }$ will also be located in the superposition of these point processes ${\displaystyle \textstyle {N}}$.

Clustering

The operation clustering is performed when each point ${\displaystyle \textstyle x}$ of some point process ${\displaystyle \textstyle {N}}$ is replaced by another (possibly different) point process. If the original process ${\displaystyle \textstyle {N}}$ is a Poisson point process, then the resulting process ${\displaystyle \textstyle {N}_{c}}$ is called a Poisson cluster point process.

Random displacement

A mathematical model may require randomly moving points of a point process to other locations on the underlying mathematical space, which gives rise to a point process operation known as displacement ^[134] or translation.^[135] The Poisson point process has been used to model, for example, the movement of plants between generations, owing to the displacement theorem,^[134] which loosely says that the random independent displacement of points of a Poisson point process (on the same underlying space) forms another Poisson point process.

Mapping

Another property that is considered useful is the ability to map a Poisson point process from one underlying space to another space.^[136]

Clumping heuristic

One method for approximating random events or phenomena with Poisson processes is called the clumping heuristic.^[140] The general heuristic or principle involves using the Poisson point process (or Poisson distribution) to approximate events, which are considered rare or unlikely, of some stochastic process. In some cases these rare events are close to being independent, hence a Poisson point process can be used. When the events are not independent, but tend to occur in clusters or clumps, then if these clumps are suitably defined such that they are approximately independent of each other, then the number of clumps occurring will be close to a Poisson random variable ^[139] and the locations of the clumps will be close to a Poisson process.^[140]

Stein's method

Stein's method is a mathematical technique originally developed for approximating random variables such as Gaussian and Poisson variables, which has also been applied to point processes. Stein's method can be used to derive upper bounds on probability metrics, which give way to quantify how different two random mathematical objects vary stochastically.^[138][141] Upperbounds on probability metrics such as total variation and Wasserstein distance have been derived.^[138]

Researchers have applied Stein's method to Poisson point processes in a number of ways,^[138] such as using Palm calculus.^[108] Techniques based on Stein's method have been developed to factor into the upper bounds the effects of certain point process operations such as thinning and superposition.^[142][143] Stein's method has also been used to derive upper bounds on metrics of Poisson and other processes such as the Cox point process, which is a Poisson process with a random intensity measure.^[138]

Poisson-type random measures

The Poisson-type random measures (PT) are a family of three random counting measures which are closed under restriction to a subspace, i.e. closed under Point process operation#Thinning. These random measures are examples of the mixed binomial process and share the distributional self-similarity property of the Poisson random measure. They are the only members of the canonical non-negative power series family of distributions to possess this property and include the Poisson distribution, negative binomial distribution, and binomial distribution. The Poisson random measure is independent on disjoint subspaces, whereas the other PT random measures (negative binomial and binomial) have positive and negative covariances. The PT random measures are discussed^[146] and include the Poisson random measure, negative binomial random measure, and binomial random measure.

Poisson point processes on more general spaces

For mathematical models the Poisson point process is often defined in Euclidean space,^[1][36] but has been generalized to more abstract spaces and plays a fundamental role in the study of random measures,^[147][148] which requires an understanding of mathematical fields such as probability theory, measure theory and topology.^[149]

In general, the concept of distance is of practical interest for applications, while topological structure is needed for Palm distributions, meaning that point processes are usually defined on mathematical spaces with metrics.^[150] Furthermore, a realization of a point process can be considered as a counting measure, so points processes are types of random measures known as random counting measures.^[115] In this context, the Poisson and other point processes have been studied on a locally compact second countable Hausdorff space.^[151]