Abstract. The problem of separating touching or overlapping objects is classical in imaging. Many solutions have been proposed in 2D. While similar, the problem in 3D has differentiating features. Contrary to the 2D case, apparent overlap due to projection effects does not exist, but real or apparent interpenetration can occur only due to either physical particle fusion or partial volume effects. Often the ability to separate objects logically is sufficient, however sometimes finding the orientation of tangent separating plane is useful. In this article, we propose a method for separating 3D touching objects and estimate a precise separating plane based on power watershed. Power watershed is used in two steps, first to obtain individual object identification, and in a second step to compute a probability field allowing sub-voxel accuracy in the plane fitting procedure. We show that our procedure is much more precise than a simple segmentation approach. We illustrate this in an application involving the shearing of a sample of sand grains imaged in various configurations by micro-CT tomography. Our technique allows the measurement of the orientation of the contacts between grains, a quantity that theoretical developments in soil mechanics explicitly use, but which has up until now been difficult to measure from experiments.