Methods for distributing a random variable by spatial interpolation with statistical corrections. The method includes assigning a numerical value of the random variable at each vertex of an array of equilateral triangles formed in a planar coordinate frame and defining a plurality of test points at respective spatial locations in the planar coordinate frame that are bounded by the array of equilateral triangles. A numerical value of the random variable is distributed at each of the test points by spatial interpolation from one or more of the numerical values of the random variable assigned at each vertex of the array of equilateral triangles. The method further includes adjusting the numerical value of the random variable distributed at each of the test points with a respective correction factor.