Returns N (approximatively) evenly spread points over the unit disc (aka. Vogel's method)

README.md

Golden dot spiral on the disc

This code sample demonstrates how to generates a set of points that covers a disc almost evenly, for any number of points. The trick is to use a Fermat's spiral and the golden angle.

42 points	256 points	999 points

Such points sets are useful for anything that requires sampling a disc.

spiral-disc.py

import numpy

def get_spiral_disc(n):
  golden_angle = numpy.pi * (3 - numpy.sqrt(5))
  theta = golden_angle * numpy.arange(n)
  radius = numpy.sqrt(numpy.arange(n) / float(n))
  return radius[:,None] * numpy.stack([numpy.cos(theta), numpy.sin(theta)], axis = 1)