We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent b2de417 commit 9464532Copy full SHA for 9464532
1 file changed
oceanstream/echodata/denoise/impulse_noise.py
@@ -172,6 +172,9 @@ def impulse_noise_mask(
172
range_values = range_values.compute()
173
elif range_values.dims:
174
range_dim = range_values.dims[0]
175
+ # Materialise to avoid dask median NaN issues
176
+ if hasattr(range_values, "compute"):
177
+ range_values = range_values.compute()
178
else:
179
range_dim = range_coord
180
@@ -184,6 +187,8 @@ def impulse_noise_mask(
184
187
window_m = float(s.rstrip("mM "))
185
188
dz_arr = range_values.diff(range_dim)
186
189
dz = float(dz_arr.median()) if dz_arr.size > 0 else 0.1
190
+ if np.isnan(dz) or dz == 0:
191
+ dz = 0.1
192
bin_sz = max(1, int(round(window_m / abs(dz))))
193
194
bin_sz = max(1, int(s))
0 commit comments