Skip to content

Commit

Permalink
Make frequency filtering robust to missing data
Browse files Browse the repository at this point in the history
  • Loading branch information
carlosgjs committed Jul 18, 2023
1 parent c03b95b commit 74f569d
Showing 1 changed file with 4 additions and 1 deletion.
5 changes: 4 additions & 1 deletion src/noisepy/seis/S1_fft_cc_MPI.py
Original file line number Diff line number Diff line change
Expand Up @@ -291,8 +291,11 @@ def _filter_channel_data(
tuples: List[Tuple[Channel, ChannelData]], samp_freq: int
) -> List[Tuple[Channel, ChannelData]]:
frequencies = set(t[1].sampling_rate for t in tuples)
frequencies = list(filter(lambda f: f >= samp_freq, frequencies))
if len(frequencies) == 0:
return []
closest_freq = min(
filter(lambda f: f >= samp_freq, frequencies),
frequencies,
key=lambda f: max(f - samp_freq, 0),
)
filtered_tuples = list(filter(lambda tup: tup[1].sampling_rate == closest_freq, tuples))
Expand Down

0 comments on commit 74f569d

Please sign in to comment.