This paper proposes a novel method for modeling human retinal cone distribution. It is based on Blue-noise sampling algorithms that share interesting properties with the sampling performed by the mosaic formed by cone photoreceptors in the retina. Here we present the method together with a series of examples of various real retinal patches. The same samples have also been created with alternative algorithms and compared with plots of the center of the inner segments of cone photoreceptors from imaged retinas. Results are evaluated with different distance measure used in the field, like nearest-neighbor analysis and pair correlation function. The proposed method can describe features of a human retinal cone distribution with a certain degree of similarity to the available data and can be efficiently used for modeling local patches of retina.