Based on the modulation ordering, a novel phase unwrapping algorithm for measurement of photoelastic fringe orders in the whole field of a model is put forward. The path for phase unwrapping is determined by the intensity modulation in the pixels on the boundary of the phase unwrapped areas. The algorithm is very superior, and in the worst case the error, if any, is limited to local minimum areas.