Error: Boolean index shape mismatch (163 vs 287)

[LITIDisma]

Hello, i loaded a mapped reference line into QGIS, along with 84 transects. In the configuration, set settings['max_dist_ref'] = 500. When running the script, an error appears in the reject_outliers function, specifically:

File "C:\Coastsat\CoastSat-master\Atacama.py", line 397, in
cross_distance = SDS_transects.reject_outliers(cross_distance,output,settings_outliers)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "C:\Coastsat\CoastSat-master\coastsat\SDS_transects.py", line 456, in reject_outliers
chainage[~np.array(idx_kept)] = np.nan
~~~~~~~~^^^^^^^^^^^^^^^^^^^^^
IndexError: boolean index did not match indexed array along dimension 0; dimension is 163 but corresponding boolean dimension is 287

I'm not sure what’s causing this mismatch or what to check to resolve it. Any guidance would be greatly appreciated.

Thank you very much in advance.

[kvos]

can you print a few of the variables, like chainage and idx_kept? and their length? then I can troubleshoot

[LITIDisma]

Excuse me, what should I do to print those variables?

[mwu752]

I believe the bug is caused by the .pkl data file, which still stores all 287 retrieved coastlines. After running SDS_tools.remove_duplicates(output) and SDS_tools.remove_inaccurate_georef(output, X), there are 124 (287-163) coastlines are removed, but these deletions are not actually deleted in the original .pkl file. I encountered the same issue

[kvos]

thanks for resolving this one @mwu752