“…In the first approach, three interferograms with the same size of the original one, I(x,y,0), I(x,y,π/2), and I(x,y,π) are generated by extracting the values of the original interferogram periodically at locations with the same phase shift, and using linear interpolation to allocate the empty pixels. Interpolation is performed by averaging the values of adjacent pixels in a similar way as is done in Refs [24][25][26]. Note that the interpolation is more accurate for the interferogram with a π/2 shift, as we have information in two pixels of the unit cell.…”