uw_grid_wrapped.m

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

function []=uw_grid_wrapped(ph_in,xy_in,pix_size,prefilt_win,goldfilt_flag,lowfilt_flag)
%UW_GRID_WRAPPED resample unwrapped phase to a grid and filter
%
%   Andy Hooper, June 2006
%

fprintf('Resampling phase to grid...\n')

if nargin<2
    error('not enough arguments')
end
if nargin<3
    pix_size=200
end
if nargin<4
    prefilt_win=32
end
if nargin<5
    goldfilt_flag='y'
end
if nargin<6
    lowfilt_flag='y'
end

[n_ps,n_ifg]=size(ph_in);
 
xy_in(:,1)=[1:n_ps]';

grid_x_min=min(xy_in(:,2));
grid_y_min=min(xy_in(:,3));

grid_ij(:,1)=ceil((xy_in(:,3)-grid_y_min+1e-3)/pix_size);
grid_ij(grid_ij(:,1)==max(grid_ij(:,1)),1)=max(grid_ij(:,1))-1;
grid_ij(:,2)=ceil((xy_in(:,2)-grid_x_min+1e-3)/pix_size);
grid_ij(grid_ij(:,2)==max(grid_ij(:,2)),2)=max(grid_ij(:,2))-1;

n_i=max(grid_ij(:,1));
n_j=max(grid_ij(:,2));

ph_grid=zeros(n_i,n_j,n_ifg,'single');

for i=1:n_ps
    ph_grid(grid_ij(i,1),grid_ij(i,2),:)=ph_grid(grid_ij(i,1),grid_ij(i,2),:)+shiftdim(ph_in(i,:),-1);
end

nzix=sum(ph_grid~=0,3)>0;
n_ps=sum(nzix(:));
[nz_i,nz_j]=find(ph_grid(:,:,1)~=0);
xy=[[1:n_ps]',(nz_j-0.5)*pix_size,(nz_i-0.5)*pix_size];

ph=zeros(n_ps,n_ifg,'single');
ph_lowpass=ph;

for i=1:n_ifg
    ph_this=ph_grid(:,:,i);
    if strcmpi(goldfilt_flag,'y') | strcmpi(lowfilt_flag,'y')
        [ph_this_gold,ph_this_low]=wrap_filt(ph_this,prefilt_win,0.8);
        ph_lowpass(:,i)=ph_this_low(nzix);
    end
    if strcmpi(goldfilt_flag,'y')
        ph(:,i)=ph_this_gold(nzix);
    else
        ph(:,i)=ph_this(nzix);
    end
end

save('uw_grid','ph','ph_lowpass','xy','nzix','grid_x_min','grid_y_min','n_i','n_j','n_ifg','n_ps','grid_ij','pix_size')