A novel synthetical method for phase unwrapping is proposed in this paper. First an effective model to eliminate local branch points without influencing normal points is established. And then, by calculating the mean values of the unwrapped phase along multiple certain directions, the goal of path independence can be achieved. The feasibilities and the efficiencies of the approach are demonstrated by employing numerical simulations and experiments. It is shown that the proposed method also exhibits strong stability and high efficiency even in presence of serious noises and relatively large under sampled areas.