stamps.m

StaMps最新测试版

function stamps(start_step,end_step,patches_flag,est_gamma_parm)
%STAMPS Stanford Method for Persistent Scatterers
%   STAMPS(START_STEP,END_STEP,PATCHES_FLAG,EST_GAMMA_FLAG) Default is to run all steps.
%   A subset of steps may be selected with START_STEP and/or END_STEP
%   STEP 1 = Initial load of data
%   STEP 2 = Estimate gamma 
%   STEP 3 = Select PS pixels
%   STEP 4 = Weed out adjacent pixels
%   STEP 5 = Correct wrapped phase for spatially-uncorrelated look angle error and merge patches
%   STEP 6 = Unwrap phase
%   STEP 7 = Calculate spatially correlated look angle (DEM) error 
%   STEP 8 = Filter spatially correlated noise 
%   
%   PATCHES_FLAG Default 'y'. Set to 'n' to process all data as one patch
%
%   EST_GAMMA_PARM is an optional parameter passed to PS_EST_GAMMA_QUICK
%
%   If current directory is a single patch, stamps only operates in the
%   current directory, but if current directory contains many patches,
%   stamps operates on them all.
%
%   Andy Hooper, June 2006
%
%   =================================================================
%   07/2006 AH: END_STEP added
%   09/2006 AH: ps_load removed (obsolete)
%   09/2006 AH: small baselines added 
%   11/2006 AH: patches added
%   01/2007 AH: calculate spatially correlated look angle error added
%   =================================================================

nfill=40;
fillstr=[repmat('#',1,nfill),'\n'];
skipstr='\n';
msgstr=fillstr;

msgstr(round(nfill)/2-11:round(nfill/2)+12)='   StaMPS Version 3.0.3 ';
fprintf([skipstr,fillstr,msgstr]);
msgstr(round(nfill)/2-11:round(nfill/2)+12)=' Andy Hooper, Jul 2008  ';
fprintf([msgstr,fillstr,skipstr]);


quick_est_gamma_flag=getparm('quick_est_gamma_flag');
unwrap_method=getparm('unwrap_method');
unwrap_prefilter_flag=getparm('unwrap_prefilter_flag');
small_baseline_flag=getparm('small_baseline_flag');

if nargin<1 
    start_step=1;
end

if nargin<2
    end_step=8;
end

if nargin<3
    if start_step<6
        patches_flag='y';
    else
        patches_flag='n';
    end
end

if nargin<4
    est_gamma_parm=0;
end

if strcmpi(patches_flag,'y')
    if exist('patch.list','file')
        fid=fopen('patch.list');
        i=0;
        while 1
            nextline=fgetl(fid);
            if ischar(nextline)
                i=i+1;
                patchdir(i).name=nextline;
            else
                break
            end
        end
    else
        patchdir=dir('PATCH*');
    end
    if isempty(patchdir)
        patches_flag='n';
    else
        ps_parms_default
        patches_flag='y';
    end
end

if ~strcmpi(patches_flag,'y')
    patchdir(1).name='.';
    fprintf('Processing current directory only...\n')
else
    fprintf('Processing patch directories...\n')
end

currdir=pwd;

nfill=40;
fillstr=[repmat('#',1,nfill),'\n'];
msgstr=fillstr;

for i=1:length(patchdir)
    cd(patchdir(i).name)
    fprintf('%s\n',pwd)
    
    if start_step==1
        msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 1 ';
        fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
        if strcmpi(small_baseline_flag,'y')
            sb_load_initial
        else
            ps_load_initial;
        end
    elseif start_step <=4
        setpsver(1)
    end

    if start_step<=2 & end_step >=2 
        msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 2 ';
        fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
        if strcmpi(quick_est_gamma_flag,'y')
            ps_est_gamma_quick(est_gamma_parm);
        else
            ps_est_gamma(est_gamma_parm);
        end
    end

    if start_step<=3 & end_step >=3 
        msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 3 ';
        fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
        if strcmpi(quick_est_gamma_flag,'y')
            ps_select;
        else
            ps_select(1);
        end
    end

    if start_step<=4 & end_step >=4 
        msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 4 ';
        fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
        if strcmpi(small_baseline_flag,'y')
            ps_weed(0,1);
        else
            ps_weed;
        end
    end

    if start_step<=5 & end_step >=5 
        msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 5 ';
        fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
        ps_correct_phase;
        ps_calc_ifg_std;
    end

    
    cd(currdir)
end

if start_step<=5 & end_step >=5 
    if patches_flag=='y'
        ps_merge_patches
    end
end


if start_step<=6 & end_step >=6 
    msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 6 ';
    fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
    ps_unwrap
    if strcmpi(small_baseline_flag,'y')
        sb_invert_uw
    end
end

if start_step<=7 & end_step >=7 
    msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 7 ';
    fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
    if strcmpi(small_baseline_flag,'y')
        %ps_aps_from_topo(1,1)
        ps_calc_scla(1)
    end
    %ps_aps_from_topo(0,1)
    ps_calc_scla
    %ps_smooth_scla
end

if start_step<=8 & end_step >=8
    msgstr(round(nfill)/2-7:round(nfill/2)+7)=' StaMPS Step 8 ';
    fprintf([skipstr,fillstr,msgstr,fillstr,skipstr]);
    ps_scn_filt
end