matlab 缠绕相位镜像对称操作

一、算法原理

1、镜像对称

   对缠绕相位进行快速傅里叶变换，需要缠绕相位数据具有隐含的周期性。而这种隐含的周期性可以通过对缠绕相位进行镜像对
称操作实现。

网上可查到的FFT变换最小二乘相位解缠代码，写的五花八门，难理解，算法效率低不说，基本全是错的，还是各种收费。因此，本文免费给出缠绕相位镜像对称的操作。

图 1 为对缠绕相位进行镜像对称操作的示意图。

图1 缠绕相位进行镜像对称操作的示意图

   从上图可以看出，做镜像对称操作就是分别以图像的边界为对称轴，进行数据的翻转，最后再进行二次翻转，具体的镜像对称操作的公式如下：

  由上式可以得到一个图像大小几乎为原缠绕相位四倍的区域，该图像中每一个像元的相位数值都与原缠绕相位息息相关。并且，进行镜像对称操作之后，该区域已经具有隐含的周期性，具备进行 FFT 变换的条件。需要注意的是，镜像对称操作不但适用于缠绕相位，而且还要应用于真实相位。

2、参考文献

[1] 邓晓龙. InSAR相位解缠算法的研究[D]. 哈尔滨工业大学, 2013.

二、代码实现

clc;
clear;
close all;

%% --------------------------------初始相位--------------------------------
phi0 = peaks(200)*3;
figure(1),
surfl(phi0);
%% -------------------------------包裹相位---------------------------------
phi = angle(exp(1i*phi0));
figure(2);
imagesc(phi);
title('原始缠绕相位');
axis on
%% ------------------------------相位解包裹--------------------------------

[m,n] = size(phi);
% 镜像对称（直接用计算公式）
% varphi = zeros(2*m,2*n);
% varphi(1:m,1:n) = phi;                        % 左上
% varphi(m+1:end,1:n) = phi(m:-1:1,:);          % 左下
% varphi(1:m,n+1:2*(n-1)) = phi(:,n-1:-1:2);    % 右下
% varphi(m+1:end,n+1:end) = phi(m:-1:1,n:-1:1); % 右下
% 镜像对称（调用函数）
varphi = zeros(2*m,2*n);
varphi(1:m,1:n) = phi;                 % 不变
varphi(m+1:end,1:n) = flipud(phi);     % 从上向下翻转
varphi(1:m,n+1:end) = fliplr(phi);     % 从左向右翻转
varphi(m+1:end,n+1:end) = rot90(phi,2);% 旋转90°

figure(3)
imagesc(varphi)
title('镜像对称');
axis on

三、结果展示