Grain-size distribution in consolidated volcaniclastic rocks can be determined by image analisys of field images, or of polished slabs or of thin sections. In order to convert 2D data, derived by the processing of a binarized (B/W) image tough ImageJ or NHI image software, in 3D data, you have to apply a stereological correction. You can do this trough CSD Correction (Higgins) or other software. Now, assuming that all the particles have the same density, you can put volume fraction data on SFT software (Wohletz) for determining grain-size parameters. This method has been commonly used for determining crystal size distribution on various kind of rocks. You can find an example of an application on consolidated volcaniclastic rocks and some explanation in Mundula et al. (2013- Journal of Volcanology and Geothermal Research).