Error in processing ALOS-2 images with STAMPS

[橋本学]

I am trying to process ALOS-2/PALSAR-2 images to obtain time series of surface deformation. I used gamma_resampled_setup to convert resampled SLC's processed with Gamma and obtained interferograms. When I input these results into StaMPS (v4.01b), I encountered the error around steps 5 and 6 as attached below.

I would appreciate for any comments.

Best regards,

Manabu Hashimoto, 

DPRI, Kyoto University

 

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

STAMPS: ################ Step 5 ################

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

STAMPS: Directory is PATCH_9

PS_CORRECT_PHASE: Starting

Correcting phase for look angle error...

GETPARM: small_baseline_flag='n'

PS_CORRECT_PHASE: Finished

STAMPS: Directory is INSAR_20160418

PS_MERGE_PATCHES: Starting

Merging patches...

GETPARM: merge_resample_size='30'

GETPARM: merge_standard_dev=Inf 

   Processing PATCH_1

Operands to the || and && operators must be convertible to logical scalar values.

Error in ps_merge_patches (line 121)

    elseif (grid_size ~=0) && (ix_no_ps~=1)

Error in stamps (line 473)

            ps_merge_patches

 

>> stamps(5)

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

STAMPS: ####### StaMPS/MTI Version 4.0b6 #######

STAMPS: #######  Beta version, Jun 2018  #######

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

STAMPS: Will process current directory only

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

STAMPS: ################ Step 5 ################

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

STAMPS: Directory is PATCH_1

PS_CORRECT_PHASE: Starting

Correcting phase for look angle error...

GETPARM: small_baseline_flag='n'

PS_CORRECT_PHASE: Finished

Estimating noise standard deviation (degrees)...

  1 18-Apr-2016 53.85

  2 02-May-2016 44.56

  3 16-May-2016 30.51

  4 13-Jun-2016 37.70

  5 27-Jun-2016 31.91

  6 11-Jul-2016 35.89

  7 25-Jul-2016 31.63

  8 08-Aug-2016 30.73

  9 05-Sep-2016 29.47

 10 19-Sep-2016 35.24

 11 03-Oct-2016 25.28

 12 17-Oct-2016 24.63

 13 31-Oct-2016 23.82

 14 14-Nov-2016 29.48

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

STAMPS: ################ Step 6 ################

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

STAMPS: Directory is PATCH_1

PS_UNWRAP: Starting

Phase-unwrapping...

GETPARM: small_baseline_flag='n'

GETPARM: unwrap_patch_phase='n'

GETPARM: scla_deramp='y'

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.238404 

PS_UNWRAP: n_trial_wraps=0.146528

Resampling phase to grid...

   Number of interferograms  : 13

   Number of points per ifg  : 310404

Error using uw_grid_wrapped (line 84)

Minimum dimension of the resampled grid (18 pixels) is less than prefilter window size (32)

Error in uw_3d (line 155)

uw_grid_wrapped(ph,xy,options.grid_size,options.prefilt_win,options.goldfilt_flag,options.lowfilt_flag,options.gold_alpha,options.ph_uw_predef);

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 503)

        ps_unwrap

[David Bekaert]

Error in step 5 is a logical error,  i.e. use of & versus && in the code, you could try and change that.

I though this was already fixed on github.

Step 6 is related to unwrap grid size at 200 m leaving only 18 pixels which the grid filter is 32 pixels.

You could decrease the unwrap grid size and or the filter size.

--
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 view this discussion on the web visit https://groups.google.com/d/msgid/mainsar/418789cb-7329-4a98-bb30-ea9eb4318b10%40googlegroups.com.