A coupling method incorporating digital image processing and discrete element method for modeling of geomaterials