get stars from multiple nside 4096 pixels rather than 256 (smaller area -> more computationally feasible)
randomize positions so that the stellar density is homogeneous (density varies by healpix in galactic coord but the catalog is store in healpix in icrs)