Discrete cosine transform (DCT) is introduced to solve discrete Poisson equation which has Neumann boundary condition.Phase unwrapping results are obtained under least-squares method meaning.