Here, we present a fast algorithm for two-dimensional (2D) phase unwrapping which behaves as a recursive linear filter. This linear behavior allows us to easily find its frequency response and stability conditions. Previously, we published a robust to noise recursive 2D phase unwrapping system with smoothing capabilities. But our previous approach was rather heuristic in the sense that not general 2D theory was given. Here an improved and better understood version of our previous 2D recursive phase unwrapper is presented. In addition, a full characterization of it is shown in terms of its frequency response and stability. The objective here is to extend our previous unwrapping algorithm and give a more solid theoretical foundation to it.