Imaging multipole gravity anomaly sources by 3D probability tomography