uw_unwrap_from_grid.m

来自「StaMps最新测试版」· M 代码 · 共 26 行

function [ph_uw]=uw_unwrap_from_grid(ph,xy,pix_size)
%UW_UNWRAP_FROM_GRID unwrap PS from unwrapped gridded ifgs 
%
%   Andy Hooper, June 2007
%
fprintf('Unwrapping from grid...\n')

uw=load('uw_grid');
uu=load('uw_phaseuw');

% grid_ij=zeros(size(xy,1),2);
% grid_ij(:,1)=ceil((xy(:,3)-uw.grid_y_min+1e-3)/pix_size);
% grid_ij(grid_ij(:,1)>uw.n_i,1)=uw.n_i;
% grid_ij(:,2)=ceil((xy(:,2)-uw.grid_x_min+1e-3)/pix_size);
% grid_ij(grid_ij(:,2)>uw.n_j,2)=uw.n_j;

[n_ps,n_ifg]=size(ph);
gridix=zeros(size(uw.nzix));
gridix(uw.nzix)=[1:uw.n_ps];
    
ph_uw=zeros(n_ps,n_ifg);

for i=1:n_ps
    ph_uw_pix=uu.ph_uw(gridix(uw.grid_ij(i,1),uw.grid_ij(i,2)),:);
    ph_uw(i,:)=ph_uw_pix+angle(ph(i,:).*exp(-j*ph_uw_pix));
end