Samples are biased

[jeffkaufman]

In generating samples, makedots.py does:

for i in range(pop):
  while True:
    samplepoint = make_ogr_point( uniform(ll,rr), uniform(bb,tt) )
    if geom.Intersects( samplepoint ):
      break
  x, y = merc.LatLonToMeters( samplepoint.GetY(), samplepoint.GetX() )

This samples uniformly in lat-lng space, but that's not actually quite uniform over the surface of the globe since that's a sphere. I think it's too minor to be worth fixing, but it was still bugging me.

[robert-chiniquy]

Nice observation. If you can come up with a fix for that, go for it.