I'm looking for a formula (isn't everyone?).

this/these formula/e evaluate/s a position on the surface of a sphere, such that...

... the entire sphere is covered.

... the density of the sphere is equal at every point.

... the input is (prefferably but not neccesarily) in the range [0,1).

I have no idea where to start, any answer or even simply a start in the right direction would be wildly apreciated.

TIA

James

What follows requires two (rather than one) uniformly distributed random variables?

Define Sphericla coordinates as follows: the sphere has radius r; the 'longitude' is denoted by q (and takes the range [0,2p]); the latitude is denoted by j and is in [-p/2,p/2].

For a given r and over all q, if j lies in the range [j,j+dj] then we are describing a 'ring' of the sphere with area

2prcos(j)rdj = 2pr²cos(j)dj

Therefore the area covered for j in [-p/2,j] is (integrating over j):

A=2pr²[sin(j)+1]

As an aside, if we set j = p/2, we get A=4pr² which is the familiar formula for the area of a sphere.

Now if we want a mapping between j and x (in [0,1]) such that the area is proportional to dx, then

4fr² x=2pr²[sin(j)+1]

j = sin^-1(2x-1)

Finally, we observe that the area of a small piece of the sphere (described by r, [q,q+dq], [j,j+dj] is proportional to dq, therefore if y is in [0,1] our second mapping is:

q = 2py

If you really must have just a single uniformly distributed variable, z, in [0,1] as an 'input' then I suppose you could construct the x, y as being the odd and even digits of z after the decimal point, but if you were ever going to write some sort of a computer program, it would make no sense.

Andre

That's really cool.

You seemed to brush over the proof that this has an uniform density. I'm sure it does and if you don't reply i'll go away beleiving it to be true, maybe even try to myself. It would just be satisfying to see a proof.

Thanks again for your time, it's much appreciated.

James

We constructed x, y such that

(i) x maps onto j and y maps onto q. There is no dependence of q on x etc.
(ii) The area bounded by x in [0,k] is proportional to k.

Alternatively, having constructed x, y we can work the other way and compute the area between (x,y) and (x+dx,y+dy). We know that the area between (r,q, j) and (r,q+dq,j+ dj) is given by:

dA=r²cos(j)djdq

and that

j = sin^-1(2x-1)

q = 2py

implies

cos(j)dj = 2dx

dq = 2pdy

substituting:

dA=r² 2dx2pdy=4pr²dxdy

The lack of this dependence of this expression on x, y, q, j is one and the same as saying the area described on the sphere is proportional to the area described in the unit square.

As an aside, you could use this mapping to map the earth on a unit square. Such an area-preserving transformation is known as (guess what?) an ëqual area projection". If you search in Google for "Wolfram equal area projection" you'll see other ways of doing it. I can't see the mapping we described above.

Andre

excellent.

thanks again.