Constrained inverse volume rendering for planetary nebulae