最小二乘法实现位相解包裹

采用最小二乘法，通过MATLAB编程语言实现位相解包裹。

% *******************************************
clear all
close all
clc
%% *************初始相位**************
N = 512;
G = 2;
phi0 = peaks（N）*G; %模拟初始相位
figure（1）
surf（phi0‘FaceColor‘‘interp‘ ‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight left axis tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘Initial Phase‘‘FontSize‘14）
set（figure（1）‘name‘‘Initial Phase 3D‘‘Numbertitle‘‘off‘）;
phi = angle（exp（j*phi0））;         %包裹相位
figure（2）;
imshow（phi[]）;
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;%title（‘Wrapped Phase‘‘FontSize‘14）
set（figure（2）‘name‘‘Wrapped Phase‘‘Numbertitle‘‘off‘）;
axis on
%% *************相位解包裹**************
[mn] = size（phi）;
phidx=zeros（mn）;
phidy=zeros（mn）;
phidx（1:m-1:）= angle（exp（j*（phi（2:m:）-phi（1:m-1:））））;
phidy（:1:n-1）= angle（exp（j*（phi（:2:n）-phi（:1:n-1））））;
%********************对包裹相位求二阶偏微分**************
Rou3 = zeros（mn）;
Rou3dx = zeros（mn）;
Rou3dy = zeros（mn）;
Rou3dx（1:m-1:） = phidx（2:m:）-phidx（1:m-1:）;
Rou3dy（:1:n-1） = phidy（:2:n）-phidy（:1:n-1）;
Rou3 = Rou3dx + Rou3dy;
figure（3）;
surf（Rou3‘FaceColor‘‘interp‘ ‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight left axis tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘lou3‘‘FontSize‘14）
set（figure（3）‘name‘‘R（xy） 3D‘‘Numbertitle‘‘off‘）;
figure（4）;
imshow（Rou3[]）;
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;
set（figure（4）‘name‘‘R（xy） 2D‘‘Numbertitle‘‘off‘）;
%% ***********************DCT求解泊松方程********************
tic
PP3 = dct2（Rou3）;
for ii=1:m
    for jj=1:n
        k1=2*cos（（ii-1）*pi/（m））;
        k2=2*cos（（jj-1）*pi/（n））;
        KK = k1+k2-4;
        PH3（iijj） = PP3（iijj）/KK;
    end
end
PH3（11） = -（PH3（12） + PH3（21） - PP3（11））/2;
phi3 = idct2（PH3）;
toc
phi3 = phi3（1:m1:n）; %解包裹出的相位


figure（5）;
surf（phi3‘FaceColor‘‘interp‘ ‘EdgeColor‘‘none‘‘FaceLighting‘‘phong‘）;
camlight left axis tight
xlabel（‘X/Pixels‘‘FontSize‘14）;ylabel（‘Y/Pixels‘‘FontSize‘14）;zlabel（‘Phase/Radians‘‘FontSize‘14）;%title（‘BLS Phase Unwrapping‘‘FontSize‘14）
set（figure（5）‘name‘‘BLS Phase Unwrapping‘‘Numbertitle‘‘off‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2360  2018-04-16 15:47  zuixiaoerchengfa.m