Sampling points with blue noise spectral characteristics are best suited for antialiasing. Compared to Poisson disk sampling, Lloyd's relaxation scheme yields samples that approximate these superior frequency characteristics much better. We complete the elementary algorithm by new techniques that allow to generate hierarchical point sets for adaptive sampling that can be efficiently tiled over the sampling plane in an aperiodic way while preserving the blue noise characteristics.