Hi,
My problem is that I want to use SSP for blink detection, but no epochs can be found. This is due to one large amplitude artefact in my dataset, which increases the blink detection threshold, as it is defined as (max(eog) - min(eog)) / 4 by default.
The artefact is annotated as bad, and when I use the function find_eog_events directly with reject_by_annotation=True, it works perfectly.
When using compute_proj_eog, however, I cannot specify reject_by_annotation. The function find_eog_events is called by compute_proj_eog with reject_by_annotation set to the default, False.
I think this could be easily fixed by offering reject_by_annotation as an argument for the function compute_proj_eog.
I am wondering, is there a reason this argument does not exist for the SSP functions? Is there something I am generally doing wrong?
If not, could this be useful for others and shall I open a pull request with the changes?