Xjyuk

The phrase "to the numbers born"?

Is it okay to consider publishing in my first year of PhD?

Is it a good practice to use a static variable in a Test Class and use that in the actual class instead of Test.isRunningTest()?

Can we generate random numbers using irrational numbers like π and e?

How to obtain a position of last non-zero element

Why not take a picture of a closer black hole?

The difference between dialogue marks

Dropping list elements from nested list after evaluation

What to do when moving next to a bird sanctuary with a loosely-domesticated cat?

How to notate time signature switching consistently every measure

Can you cast a spell on someone in the Ethereal Plane, if you are on the Material Plane and have the True Seeing spell active?

Is Cinnamon a desktop environment or a window manager? (Or both?)

Are there any other methods to apply to solving simultaneous equations?

For what reasons would an animal species NOT cross a *horizontal* land bridge?

Keeping a retro style to sci-fi spaceships?

Unitary representations of finite groups over finite fields

Is an up-to-date browser secure on an out-of-date OS?

Match Roman Numerals

Why isn't the circumferential light around the M87 black hole's event horizon symmetric?

Can an undergraduate be advised by a professor who is very far away?

What information about me do stores get via my credit card?

Why can't devices on different VLANs, but on the same subnet, communicate?

Why doesn't UInt have a toDouble()?

Did any laptop computers have a built-in 5 1/4 inch floppy drive?

Is there a way to generate a uniformly distributed point on a sphere from a fixed amount of random real numbers?

4

1

$begingroup$

The obvious solution of Lattitude & Longitude doesn't work because it generates points more densely near the poles, and the other thing I came up with (Pick a random point in the unit cube, if it's in the sphere map it to the surface, and restart if it's outside) doesn't always find a point within a fixed number of tries.

geometry

share|cite|improve this question

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  So what you want is a uniform distribution. It would be helpful to state this explicitly.
  $endgroup$
  – robjohn♦
  3 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Distribute longitude uniformly and the sine of the latitude uniformly. Then the distribution of points on the sphere will be uniform.
  $endgroup$
  – robjohn♦
  3 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @robjohn thank you, you're right that I forgot to specify that.
  $endgroup$
  – The Zach Man
  1 hour ago
  

  
  
  
  

add a comment | 

4

1

$begingroup$

The obvious solution of Lattitude & Longitude doesn't work because it generates points more densely near the poles, and the other thing I came up with (Pick a random point in the unit cube, if it's in the sphere map it to the surface, and restart if it's outside) doesn't always find a point within a fixed number of tries.

geometry

share|cite|improve this question

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

$endgroup$

• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  So what you want is a uniform distribution. It would be helpful to state this explicitly.
  $endgroup$
  – robjohn♦
  3 hours ago
  

  
  
  
  


• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  $begingroup$
  Distribute longitude uniformly and the sine of the latitude uniformly. Then the distribution of points on the sphere will be uniform.
  $endgroup$
  – robjohn♦
  3 hours ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  $begingroup$
  @robjohn thank you, you're right that I forgot to specify that.
  $endgroup$
  – The Zach Man
  1 hour ago
  

  
  
  
  

add a comment | 

$begingroup$

The obvious solution of Lattitude & Longitude doesn't work because it generates points more densely near the poles, and the other thing I came up with (Pick a random point in the unit cube, if it's in the sphere map it to the surface, and restart if it's outside) doesn't always find a point within a fixed number of tries.

geometry

share|cite|improve this question

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

$endgroup$

share|cite|improve this question

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

share|cite|improve this question

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

edited 56 mins ago

robjohn♦

271k27313642

edited 56 mins ago

robjohn♦

271k27313642

asked 3 hours ago

The Zach ManThe Zach Man

1107

asked 3 hours ago

The Zach ManThe Zach Man

1107

2 Answers
2

active

oldest

votes

5

$begingroup$

The Lambert cylindrical equal area projection maps the sphere to a cylinder, area to equal area. It is easy to generate a uniform distribution on a cylinder. Simply map it back to the sphere.

For $(u_1,u_2)$ uniform on $[0,1]^2$, either

$mathrmlat=arcsin(2u_1-1),mathrmlon=2pi u_2$

or

$z=2u_1-1,x=sqrt1-z^2cos(2pi u_2),y=sqrt1-z^2sin(2pi u_2)$

share|cite|improve this answer

answered 1 hour ago

robjohn♦robjohn

271k27313642

$endgroup$

add a comment | 

3

$begingroup$

Your method, even though it doesn't finish in a fixed number of times, is a reasonable way to do it. Each trial succeeds with probability $fracpi6$, which is better than $frac12$: the average number of trials is less than $2$.

Another standard method is to use the normal distribution. Generate $x, y, z$ independently from a standard normal distribution, then take the point $(x,y,z)$ and divide it by $sqrtx^2+y^2+z^2$ as you did for points inside the cube. The multivariate normal distribution is rotationally symmetric, so this will get you evenly distributed points on the sphere.

(The Box–Muller transform is one way to generate normally distributed random numbers, and some versions of it do not use rejection sampling, so they can be done with a "fixed amount" of randomness.)

share|cite|improve this answer

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

$endgroup$

add a comment | 

Your Answer

StackExchange.ifUsing("editor", function ()
return StackExchange.using("mathjaxEditing", function ()
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
);
);
, "mathjax-editing");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "69"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3184449%2fis-there-a-way-to-generate-a-uniformly-distributed-point-on-a-sphere-from-a-fixe%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

5

$begingroup$

The Lambert cylindrical equal area projection maps the sphere to a cylinder, area to equal area. It is easy to generate a uniform distribution on a cylinder. Simply map it back to the sphere.

For $(u_1,u_2)$ uniform on $[0,1]^2$, either

$mathrmlat=arcsin(2u_1-1),mathrmlon=2pi u_2$

or

$z=2u_1-1,x=sqrt1-z^2cos(2pi u_2),y=sqrt1-z^2sin(2pi u_2)$

share|cite|improve this answer

answered 1 hour ago

robjohn♦robjohn

271k27313642

$endgroup$

add a comment | 

5

$begingroup$

The Lambert cylindrical equal area projection maps the sphere to a cylinder, area to equal area. It is easy to generate a uniform distribution on a cylinder. Simply map it back to the sphere.

For $(u_1,u_2)$ uniform on $[0,1]^2$, either

$mathrmlat=arcsin(2u_1-1),mathrmlon=2pi u_2$

or

$z=2u_1-1,x=sqrt1-z^2cos(2pi u_2),y=sqrt1-z^2sin(2pi u_2)$

share|cite|improve this answer

answered 1 hour ago

robjohn♦robjohn

271k27313642

$endgroup$

add a comment | 

$begingroup$

The Lambert cylindrical equal area projection maps the sphere to a cylinder, area to equal area. It is easy to generate a uniform distribution on a cylinder. Simply map it back to the sphere.

For $(u_1,u_2)$ uniform on $[0,1]^2$, either

$mathrmlat=arcsin(2u_1-1),mathrmlon=2pi u_2$

or

$z=2u_1-1,x=sqrt1-z^2cos(2pi u_2),y=sqrt1-z^2sin(2pi u_2)$

share|cite|improve this answer

answered 1 hour ago

robjohn♦robjohn

271k27313642

$endgroup$

share|cite|improve this answer

answered 1 hour ago

robjohn♦robjohn

271k27313642

answered 1 hour ago

robjohn♦robjohn

271k27313642

answered 1 hour ago

robjohn♦robjohn

271k27313642

3

$begingroup$

Your method, even though it doesn't finish in a fixed number of times, is a reasonable way to do it. Each trial succeeds with probability $fracpi6$, which is better than $frac12$: the average number of trials is less than $2$.

Another standard method is to use the normal distribution. Generate $x, y, z$ independently from a standard normal distribution, then take the point $(x,y,z)$ and divide it by $sqrtx^2+y^2+z^2$ as you did for points inside the cube. The multivariate normal distribution is rotationally symmetric, so this will get you evenly distributed points on the sphere.

(The Box–Muller transform is one way to generate normally distributed random numbers, and some versions of it do not use rejection sampling, so they can be done with a "fixed amount" of randomness.)

share|cite|improve this answer

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

$endgroup$

add a comment | 

3

$begingroup$

Your method, even though it doesn't finish in a fixed number of times, is a reasonable way to do it. Each trial succeeds with probability $fracpi6$, which is better than $frac12$: the average number of trials is less than $2$.

Another standard method is to use the normal distribution. Generate $x, y, z$ independently from a standard normal distribution, then take the point $(x,y,z)$ and divide it by $sqrtx^2+y^2+z^2$ as you did for points inside the cube. The multivariate normal distribution is rotationally symmetric, so this will get you evenly distributed points on the sphere.

(The Box–Muller transform is one way to generate normally distributed random numbers, and some versions of it do not use rejection sampling, so they can be done with a "fixed amount" of randomness.)

share|cite|improve this answer

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

$endgroup$

add a comment | 

$begingroup$

Your method, even though it doesn't finish in a fixed number of times, is a reasonable way to do it. Each trial succeeds with probability $fracpi6$, which is better than $frac12$: the average number of trials is less than $2$.

Another standard method is to use the normal distribution. Generate $x, y, z$ independently from a standard normal distribution, then take the point $(x,y,z)$ and divide it by $sqrtx^2+y^2+z^2$ as you did for points inside the cube. The multivariate normal distribution is rotationally symmetric, so this will get you evenly distributed points on the sphere.

(The Box–Muller transform is one way to generate normally distributed random numbers, and some versions of it do not use rejection sampling, so they can be done with a "fixed amount" of randomness.)

share|cite|improve this answer

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

$endgroup$

share|cite|improve this answer

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107

answered 3 hours ago

Misha LavrovMisha Lavrov

49k757107