Function to unwrap phase having 2*pi jumps just after phase extraction from interferogram

This is the function to unwrap 2*pi jumps on extracted phase from interferogram. This function implemented TV (Total Variation) minimization algorithm (ANTONIN CHAMBOLLE, “An Algorithm for Total Variation Minimization and Applications,” J. Math. Imaging Vis., vol. 20, no. 1/2, pp. 89–97, Jan. 2004.) to minimize HN (High Noise) capable to stretch unwrapped phase obtained by using simple unwrap method. The function is based on paper, "H. Y. H. Huang, L. Tian, Z. Zhang, Y. Liu, Z. Chen, and G. Barbastathis, “Path-independent phase unwrapping using phase gradient and total-variation (TV) denoising,” Opt. Express, vol. 20, no. 13, p. 14075, Jun. 2012.".
Unlikely to other advanced phase unwrapping technique.
Use of this function is really simple as "unwrapped_phase = phase_unwrap_TV_min(wrapped_phase,residue_add_check)" and relatively fast. residue_add_check can be 'yes' or 'no'. If 'yes', HN (High Noise) originally shown in wrapped_phase is clearly distinguishable on output unwrapped phase. Sometimes, there are still 2*pi jumps on output but these jumps are easily dealt with simple unwrap method such as unwrap command in Matlab. if 'No', HN and possible 2*pi jumps becomes blurry.
I recommend to use "unwrapped_phase = phase_unwrap_TV_min(wrapped_phase,'yes')".
I appreciate authors of the paper which is really nice.