Uniformly distributed in Poisson-Voronoi cell When I read research article: here, I found following set-up:
Figure illustration of one realization of Big node (BN) positions from the homogeneous Poisson point process (PPP) and of small nodes (SN) uniformly distributed in the corresponding Poisson-Voronoi cells.
Each SN connects to its closest BN such that the coverage area of a BN is its Poisson-Voronoi cell. The SNs are assumed to be uniformly distributed in the Poisson-Voronoi cell of their serving BN.
However, I do not really understand the meaning of ``The SNs are assumed to be uniformly distributed in the Poisson-Voronoi cell of their serving BN". 
Can someone elaborate this explanation much simpler way? 
I am not sure if someone can help me of coding this one (matlab) to get such figure. 

 A: By "uniform", they mean that the point does not tend to exist more often in some regions (of each Voronoi cell) than others.  In Cartesian coordinates for a rectangular cell, this means that the x and y coordinates are independent uniform random variables, say X and Y. 
It's relatively straight forward to place a point uniformly on a simple geometric object, such as a rectangle, disk, or triangle. Those links are for pages that detail how to simulate Poisson point processes, but to do that, each point is placed uniformly. 
But it is a bit more complicated to do it for a general Voronoi cell. Here are two possible methods.
1) The simplest and crudest method is to bound a Voronoi cell with a rectangle or disk, then place/sample the point uniformly on the rectangle or disk. If it doesn't randomly land inside the rectangle or disk, then start again. Crude, slightly inefficient, but simple. 
2) The more elegant way is to first partition (ie divide) the Voronoi cells into triangles. (For a Voronoi cell, this is easy. Each side corresponds to one side of a triangle (ie two points), where the third point is the original Poisson point corresponding to the Voronoi cell.) Given the n triangles that form the Voronoi cell, randomly choose a triangle based on the areas of the triangle (ie the probabilities are the areas renormalized by the total area of the cell). In MATLAB, you would do this:
cdfTri=cumsum(areaTri)/areaCell; %create triangle CDF 
indexTri=find(rand(1)<=cdfTri,1); %use CDF to randomly choose a triangle

Then for that single triangle, uniformly place a point on it (there are expressions for this; see below). Repeat for all (bounded) Voronoi cells. 
The topic of placing a single point uniformly on a general triangle is discussed in this StackExchange post; see that for suggestions of proof.  For the formulas, it cites the paper "Shape distributions" by Osada, Funkhouser, Chazelle and Dobkin, where no proof is given. Discussions on this problem can also be found here and here. 
For Method 2, I have uploaded some code in both MATLAB and Python (using SciPy and Numpy). The main random placement functions are located respectively in the files funVoronoiUniform.m and funVoronoiUniform.py. 
I've also included two files that empirically test those functions respectively, namely  VoronoiUniformTest.m or VoronoiUniformTest.py, by empirically estimating the centroids of the Voronoi cells. To test the functions, run VoronoiUniformTest.m or VoronoiUniformTest.py with the aforementioned files in the same directory/folder. 
https://github.com/hpaulkeeler/voronoi_uniform
