Contact rolling resistance is the most widely used method to incorporate particle shape effects in the discrete element method (DEM). The main reason for this is that such approach allows for using spherical particles hence offering substantial computational benefits compared to non-spherical DEM models. This paper shows how rolling resistance parameters for 3D DEM models can be easily calibrated with 2D sand grain images.