Pore scale modeling plays a key role in fluid flow through porous media and associated macroscale constitutive relationships. The polyhedral shape and effective local pore size within granular material microstructure are computed in this study by means of the Euclidean Distance Transform (EDT), a local maxima search (non-maximum suppression), and a segmentation process. Various synthetic packed particles are simulated and employed as comparative models during the computation of pore size distribution (PSD). Reconstructed un-sheared and sheared Ottawa 20–30 sand samples are used to compute PSD for non-trivial and non-spherical models.