An MAP algorithm with edge-preserving prior for muon tomography