Unwrapping error

[Tom]

Hello,

   I used small baseline in the new version of StaMPS and encountered an error in step 6. Does anybody have idea about this problem?

>> stamps(6,6)

STAMPS: ########################################
STAMPS: ####### StaMPS/MTI Version 4.0b6 #######
STAMPS: #######  Beta version, Jun 2018  #######
STAMPS: ########################################

STAMPS: Will process current directory only

STAMPS: ########################################
STAMPS: ################ Step 6 ################
STAMPS: ########################################
STAMPS: Directory is SMALL_BASELINES

PS_UNWRAP: Starting
Phase-unwrapping...
GETPARM: small_baseline_flag='y'
GETPARM: unwrap_patch_phase='n'
GETPARM: scla_deramp='n'
GETPARM: subtr_tropo='n'
GETPARM: tropo_method='a_l'
GETPARM: drop_ifg_index=
GETPARM: unwrap_hold_good_values='n'
PS_UNWRAP: Code to hold good values skipped
GETPARM: unwrap_time_win=730
GETPARM: unwrap_method='3D'
GETPARM: unwrap_grid_size=200
GETPARM: unwrap_gold_n_win=32
GETPARM: unwrap_prefilter_flag='y'
GETPARM: unwrap_gold_alpha=0.8
GETPARM: unwrap_la_error_flag='y'
GETPARM: unwrap_spatial_cost_func_flag='n'
GETPARM: max_topo_err=20
GETPARM: lambda=0.0554658
PS_UNWRAP: n_trial_wraps=0.186141
Resampling phase to grid...
   Number of interferograms  : 113
   Number of points per ifg  : 2477123
   Number of resampled points: 38121
Your variables are reaching 2GB limit, revert to save -v7.3
This will be slower but avoids matlab not saving the data
Interpolating grid...
   Number of unique edges in grid: 86646
Unwrapping in time-space...
   Estimating look angle error (elapsed time=170s)
   Smoothing in time (elapsed time=174s)
   Performing complex smoothing on 630 edges (elapsed time=188s)
Unwrapping in space...
   Processing IFG 1 of 113
The logical indices contain a true value outside of the array bounds.

Error in uw_stat_costs (line 167)
    ph_uw(:,i1)=ifguw(uw.nzix);

Error in uw_3d (line 163)
uw_stat_costs(options.unwrap_method,options.variance);

Error in ps_unwrap (line 235)
    [ph_uw_some,msd_some]=uw_3d(ph_w(:,unwrap_ifg_index),ps.xy,day,ifgday_ix(unwrap_ifg_index,:),ps.bperp(unwrap_ifg_index),options);
   
Error in stamps (line 504)
        ps_unwrap

[David Bekaert]

The error shows that one of the value of nzix is outside the matrix size.

What are the values of uw.nzix,  min(uw.nzix), max(uw.nzix) and size(ifguw).

Can you provide more background on how the data was prepared.

Thanks,

D.

--
You received this message because you are subscribed to the Google Groups "MAINSAR" group.
To unsubscribe from this group and stop receiving emails from it, send an email to mainsar+u...@googlegroups.com.
To post to this group, send email to mai...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Tom]

Hi David,

   Thank you. I dont know how did this problem occur. When finished running step 6, I just changed the parameter 'merge_resample_size' in step 5 and rerun steps 5 and  6 and this error occur. So I just deleted the unwrapping results phuw*.mat then rerun steps 5 and 6, and this error disappeared.

在 2018年10月29日星期一 UTC+1上午1:53:24，David Bekaert写道：