Assume you have a Poisson point process of constant intensity $\lambda$ in the Euclidean plane. From this point process we construct the Delaunay triangulation (or the Voronoi tessellation for that matter). 


It is known [Stoyan et all] that the expected degree of a typical vertex has degree 6. Moreover, there are several results for the average area of a typical triangle, edge lengths and so on. However, all the results seem to be for constant intensity. My question is: Is it known how the expected degree changes as we change the intensity?


More specifically, assume we have a rotationally invariant intensity $\rho(r)$ where $r$ is the distance to the origin and assume that $\rho(r)\to\infty$ as $r$ increases. How, does the expected degree of a node depends on its distance to the origin. The intuition is that the degree is going to increase as $r$ increases but is there any known result or reference in this direction? This is not my research area so I will appreciate any help or comment!



Expected Degree of a vertex in Delaunay Triangulations
The expected degree for the Delaunay triangulation on hyperbolic space will depend on the density of the points. With low enough density points, you should get arbitrarily high degree. You should be able to get the expected degree as high as you want for a point distribution in the plane by taking the conformal representation of the hyperbolic plane as a varying metric in a unit disc (the Poincare disc), and placing the points with density proportional to this metric. Then the density goes to infinity at the boundary of the circle.


Note that since the Poincare disc model takes circles in the hyperbolic plane to circles in the disc, the Delaunay triangulation of a set of points (which is characterized by the fact that the circumcircle of every triangle does not contain another point) is the same in the disc and the hyperbolic plane. Thus, the expected average degree of the Delaunay triangulation should be the same in both of these scenarios. 


The above construction satisfies the conditions stated in your question, but I don't know if you'll be satisfied. Did you want the points to be distributed on the entire plane? If so, I expect it's impossible because there's no way to put a conformal metric on all of $\mathbb{R}^2$ that gives the hyperbolic plane (although note that this isn't quite a proof). 


