SARvey - ValueError: attempt to get argmax of an empty sequence after all PSCs are removed in Step 1

[epehlivanli]

2025-06-10 18:01:54,535 - INFO - SARvey version: 1.2.0 - Strawberry Pie, 2025-02-19_01, Run: MTInSAR
2025-06-10 18:01:54,536 - INFO -                               Parameter           value         default
2025-06-10 18:01:54,537 - INFO -                               _________           _____         _______
2025-06-10 18:01:54,537 - INFO -                              input_path         inputs/         inputs/
2025-06-10 18:01:54,537 - INFO -                             output_path /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs <---   outputs/
2025-06-10 18:01:54,537 - INFO -                               num_cores              10 <---         50
2025-06-10 18:01:54,537 - INFO -                             num_patches               1               1
2025-06-10 18:01:54,537 - INFO -               apply_temporal_unwrapping            True            True
2025-06-10 18:01:54,537 - INFO -               spatial_unwrapping_method            puma            puma
2025-06-10 18:01:54,537 - INFO -                           logging_level            INFO            INFO
2025-06-10 18:01:54,537 - INFO -                            logfile_path       logfiles/       logfiles/
2025-06-10 18:01:54,537 - INFO - 
2025-06-10 18:01:54,541 - INFO -     ---------------------------------------------------------------------------------
2025-06-10 18:01:54,541 - INFO -                          STEP 1:     CONSISTENCY CHECK
2025-06-10 18:01:54,541 - INFO -     ---------------------------------------------------------------------------------
2025-06-10 18:01:54,541 - INFO -                               Parameter           value         default
2025-06-10 18:01:54,541 - INFO -                               _________           _____         _______
2025-06-10 18:01:54,541 - INFO -                            coherence_p1             0.7 <---        0.9
2025-06-10 18:01:54,541 - INFO -                               grid_size            None <---        200
2025-06-10 18:01:54,541 - INFO -                            mask_p1_file            None <---           
2025-06-10 18:01:54,541 - INFO -                  num_nearest_neighbours              50 <---         30
2025-06-10 18:01:54,541 - INFO -                          max_arc_length          999999 <---       None
2025-06-10 18:01:54,541 - INFO -                          velocity_bound             0.1             0.1
2025-06-10 18:01:54,541 - INFO -                         dem_error_bound           250.0 <---      100.0
2025-06-10 18:01:54,541 - INFO -                num_optimization_samples             100             100
2025-06-10 18:01:54,541 - INFO -                arc_unwrapping_coherence             0.4 <---        0.6
2025-06-10 18:01:54,541 - INFO -                             min_num_arc               3               3
2025-06-10 18:01:54,541 - INFO - 
2025-06-10 18:01:54,623 - INFO - reading box None from file: /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs/temporal_coherence.h5 ...
2025-06-10 18:01:56,133 - INFO - No mask for area of interest given.
2025-06-10 18:01:58,057 - INFO - reading box None from file: /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs/ifg_stack.h5 ...
2025-06-10 18:02:01,960 - INFO - write data to /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs/p1_ifg_wr.h5...
2025-06-10 18:02:02,055 - INFO - create distance matrix between all points...
2025-06-10 18:02:19,185 - INFO - Triangulate points with 50-nearest neighbours.
2025-06-10 18:02:21,712 - DEBUG - time used: 00 mins 2.5 secs.
2025-06-10 18:02:21,713 - INFO - Triangulate points with global delaunay.
2025-06-10 18:02:22,306 - INFO - remove arcs with length > 999999.
2025-06-10 18:02:23,712 - INFO - retrieve arcs from adjacency matrix.
2025-06-10 18:02:24,687 - INFO - no. arcs:	614730
2025-06-10 18:02:54,263 - INFO - ifg arc observations created.
2025-06-10 18:02:54,264 - INFO - write data to /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs/point_network.h5...
2025-06-10 18:03:00,230 - INFO - ########## TEMPORAL UNWRAPPING: AMBIGUITY FUNCTION ##########
2025-06-10 18:03:00,230 - INFO - start parallel processing with 10 cores.
2025-06-10 18:03:00,295 - INFO - Finished temporal unwrapping.
2025-06-10 18:03:00,295 - DEBUG - time used: 00 mins 0.1 secs.
2025-06-10 18:03:00,296 - INFO - write data to /scratch/08590/emrhnp/MiamiTsxSMDT36/Grove_subset/Grove_newcon/outputs/point_network_parameter.h5...
2025-06-10 18:09:06,873 - INFO - Detect points with low quality arcs (mean): < 0.4
2025-06-10 18:09:06,874 - INFO - Removal of points whose arcs are incoherent in average.
2025-06-10 18:13:06,807 - INFO - Detected 21095 point(s) with mean coherence of all connected arcs < 0.4 
2025-06-10 18:13:06,807 - INFO - Removal of low quality arcs: < 0.4
2025-06-10 18:13:06,809 - INFO - Removed 614730 arc(s)
2025-06-10 18:13:06,813 - INFO - Removal of arcs and PSC that cannot be tested.
2025-06-10 18:13:06,863 - INFO - Detected 21095 point(s) with less than 3 arcs
2025-06-10 18:13:06,863 - INFO - Remove 21095 point(s)
2025-06-10 18:13:07,213 - INFO - Removed 0 arc(s) connected to the removed point(s)
2025-06-10 18:13:08,191 - INFO - ########## NOISY POINT REMOVAL BASED ON ARC PARAMETERS ##########
2025-06-10 18:13:08,191 - INFO - Selection of the reference PSC

While running SARvey with standard config.json parameters on a small area subset (new construction area Grove), I get a crash in Step 1 (Consistency Check) when all PSCs were removed, leading to an unhandled ValueError in the subsequent "noisy point removal" step.

From the log:

Detected 21095 point(s) with mean coherence of all connected arcs < 0.4 
Removed 21095 point(s)
Detected 21095 point(s) with less than 3 arcs
Remove 21095 point(s)

Then:

ValueError: attempt to get argmax of an empty sequence

This occurs at:

spatial_ref_idx = np.where(design_mat[np.argmax(net_par_obj.gamma), :] != 0)[0][0]

Plotted p1_ts.h5:

Plotted temporal_coherence.h5:

To run:

wget http://149.165.154.65:/data/HDF5EOS/Emirhan/Grove_subset/Grove_newcon.tar.gz
tar -xzf Grove_newcon.tar.gz
cd Grove_newcon
sarvey -f config.json 0 2