STAMPS step 2 stops at ps_est_gamma_quick.m line 225. Input data from ISCE, Sentinel-1.

[alexis.ap...@gmail.com]

Hello,

Could you help with the following problem?

Stamps step 2 stops mysteriously at line 225 ps_est_gamma_quick.m like it needs keyboard input. Output is copied below.

However it can't continue from there no way. 

Sentinel-1 stack is prepared with ISCE (topsStack and prepStackToStaMPS from contrib)

I didn't notice any error in the previous steps of the process (ISCE topsStack runs, make_single_master_stack_isce, mt_prep_isce,). 

When I created with "mt_prep_isce" 35 patches (7 range x 5 azimuth) instead of 6 (3 range x 2 azimuth), for some Paches step 2 was completed ok.

 

step 2 output from matlab:

STAMPS: ########################################

STAMPS: ################ Step 2 ################

STAMPS: ########################################

STAMPS: Directory is PATCH_2

PS_EST_GAMMA_QUICK: Starting

PS_EST_GAMMA_QUICK: Estimating gamma for candidate pixels

GETPARM: filter_grid_size=50 

GETPARM: filter_weighting='P-square'

GETPARM: clap_win=32 

GETPARM: clap_low_pass_wavelength=800 

GETPARM: clap_alpha=1 

GETPARM: clap_beta=0.3 

GETPARM: max_topo_err=20 

GETPARM: lambda=0.056 

GETPARM: gamma_change_convergence=0.005 

GETPARM: gamma_max_iterations=3 

GETPARM: small_baseline_flag='n'

Found inc angle file 

PS_EST_GAMMA_QUICK: n_trial_wraps=0.003267

PS_EST_GAMMA_QUICK: Initialising random distribution...

PS_EST_GAMMA_QUICK: 12499048 PS candidates to process

PS_EST_GAMMA_QUICK: iteration #1

PS_EST_GAMMA_QUICK: Calculating patch phases...

K>> 

Thanks!

Alex

[ma1...@bristol.ac.uk]

Hi

Im having this same problem. Did you solved it?

[Alex]

Hello,

I haven't found the cause of this error and I don't know under which circumstances it is happening.

I did some changes to the input though that probably affected the processing avoiding the error after that. 

I cropped the input stack to a smaller area using the very useful David Bekaert's script "run_SLCcropStack.csh" from ISCE contrib folder and I also tried to make the patches more squared (than rectangle) changing the range and azimuth analogy in mt_prep_isce params.

Alex

[ma1...@bristol.ac.uk]

Thanks for replying

The problem seems to be with the clap_filt but I don't know where to look at and check if it is my data.

STAMPS: ########################################

STAMPS: ################ Step 2 ################

STAMPS: ########################################

STAMPS: Directory is PATCH_1

PS_EST_GAMMA_QUICK: Starting

PS_EST_GAMMA_QUICK: Estimating gamma for candidate pixels

GETPARM: filter_grid_size=50 

GETPARM: filter_weighting='P-square'

GETPARM: clap_win=32 

GETPARM: clap_low_pass_wavelength=800 

GETPARM: clap_alpha=1 

GETPARM: clap_beta=0.3 

GETPARM: max_topo_err=20 

GETPARM: lambda=0.0554658 

GETPARM: gamma_change_convergence=0.005 

GETPARM: gamma_max_iterations=3 

GETPARM: small_baseline_flag='n'

Found look angle file 

PS_EST_GAMMA_QUICK: n_trial_wraps=0.295776

PS_EST_GAMMA_QUICK: Initialising random distribution...

Warning: P-file /opt/matlab-R2018a/toolbox/matlab/datafun/hist.p is older than

MATLAB file /opt/matlab-R2018a/toolbox/matlab/datafun/hist.m.

/opt/matlab-R2018a/toolbox/matlab/datafun/hist.p may be obsolete and may need

to be regenerated.

Type "help pcode" for information about generating P-files. 

> In ps_est_gamma_quick (line 174)

  In stamps (line 326) 

PS_EST_GAMMA_QUICK: 151055 PS candidates to process

PS_EST_GAMMA_QUICK: iteration #1

PS_EST_GAMMA_QUICK: Calculating patch phases...

At this point it hangs and suddenly I am able to use the command line and this message appears.

Operation terminated by user during clap_filt (line 88)

In ps_est_gamma_quick (line 225)
        ph_filt(:,:,i)=clap_filt(ph_grid(:,:,i),clap_alpha,clap_beta,n_win*0.75,n_win*0.25,low_pass);
        
In stamps (line 326)
                    ps_est_gamma_quick(est_gamma_parm);

On Monday, 25 March 2019 18:22:51 UTC, Alex wrote:

[weilin tang]

Hello, have you resolved your problem? I also encountered the same problem as you