Fix events.tsv onset off by one in bids_preproc - #315
Open
arnodelorme wants to merge 1 commit into
Open
Conversation
|
Claude encountered an error after 0s —— View job I'll analyze this and get back to you. |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
🤖
Problem
bids_preprocrewrote the derivativeevents.tsvwithev_time = EEG['times'][e['latency']] / 1000.0. Event latencies are 1-based samples (latency 1.0 = first sample, EEGLAB convention) whileEEG['times']is a 0-based array, so:times[1]= 7.8125 ms at 128 Hz instead of 0.latency == pnts) raisedIndexErrorand was dropped.sample_data/eeglab_data.set, or anything after resampling) raisedIndexErrorand was dropped. Oneeglab_data.setthe old code kept 0 of 154 events.The dropped events also masked a second failure:
pop_loadsetyieldsduration=Nonefor empty MATLAB fields, soe.get('duration', 0.0) / sratewould have raisedTypeErroronce an event got past the onset lookup.Fix
(latency - 1) / srate, the same conversion already used byplugins/EEG_BIDS/pop_exportbids.py. Notimes[...]lookup, so noIndexErrorand no dropped events. Duration staysduration / srate(EEGLAB stores duration in samples);Noneis treated as 0._write_events_tsv(EEG, fpath)so it can be tested on real data without running the pipeline.grep -rn "times\[" src/eegprep/plugins/EEG_BIDS/shows this was the only such lookup underEEG_BIDS.Epoched data
bids_preprocwritesevents.tsvafter the optionalpop_epochstage, so epoched datasets are exported too.pop_epochproduces 1-based latencies concatenated across epochs (... + 1 + pnts * epoch_index), so(latency - 1) / srategives the onset relative to the start of the concatenated epoched array, consistent withpop_exportbids. No behaviour beyond that is introduced; the helper docstring states this.Verification
New
TestWriteEventsTsv(2 tests, realsample_data/eeglab_data.set, no torch/MATLAB): all 154 events written with onset(latency-1)/srate(fractional latencies included, duration 0, type preserved), and latency 1 / latencypntsmap to 0.0 s and(pnts-1)/srate.